Genome-wide association study of seedling nitrogen-use efficiency-associated traits in common wheat(Triticum aestivum L.)

Nitrogen(N)fertilizer application is essential for crop-plant growth and development.Identifying genetic loci associated with N-use efficiency(NUE)could increase wheat yields and reduce environmental pollution caused by overfertilization.We subjected a panel of 389 wheat accessions to N and chlorate(a nitrate analog)treatments to identify quantitative trait loci(QTL)controlling NUE-associated traits at the wheat seedling stage.Genotyping the panel with a 660K single-nucleotide polymorphism(SNP)array,we identified 397 SNPs associated with N-sensitivity index and chlorate inhibition rate.These SNPs were merged into 49 QTL,of which eight were multi-environment stable QTL and 27 were located near previously reported QTL.A set of 135 candidate genes near the 49 QTL included TaBOX(F-box family protein)and TaERF(ethylene-responsive transcription factor).A Tabox mutant was more sensitive to low-N stress than the wild-type plant.We developed two functional markers for Hap 1,the favorable allele of TaBOX.

Genetic Diversity,Population Structure,and Genome-Wide Association Study of Seven Agronomic Traits in 273 Diverse Upload Cotton Accessions

Upland cotton(Gossypium hirsutum)is the most important plant producing natural fibers for the textile industry.In this study,we first investigated the phenotypic variation of seven agronomic traits of 273 diverse cotton accessions in the years 2017 and 2018,which were from 18 geographical regions.We found large variations among the traits in different geographical regions and only half of the traits in either years 2017 or 2018 followed a normal distribution.We then genotyped the collection with 81,612 high quality SNPs.Phylogenetic tree and population structure revealed a diverse genetic structure of the core collection,and geographical diversification was an important factor,but account for part of the variances of genetic diversification.We then performed genome-wide association study for the seven traits in the years 2017 and 2018,and the average values of each trait in the two years,respectively.We identified a total of 19 significant marker-trait associations and found that Pollen Ole e 1 allergen/extension could be the candidate gene associated with the fall-off cotton bolls from the last three branches.In addition,large variations were observed for the heritability of traits in the years 2017 and 2018.These results provide new potential candidate genes for further functional validation,which could be useful for genetic improvement and breeding of new cotton cultivars with better agronomic performances.

Genetic dissection of rice appearance quality and cooked rice elongation by genome-wide association study

Appearance and cooked rice elongation are key quality traits of rice. Although some QTL for these traits have been identified, understanding of the genetic relationship between them remains limited. In the present study, large phenotypic variation was observed in 760 accessions from the 3 K Rice Genomes Project for both appearance quality and cooked rice elongation. Most component traits of appearance quality and cooked rice elongation showed significant pairwise correlations, but a low correlation was found between appearance quality and cooked rice elongation. A genome-wide association study identified 74 QTL distributed on all 12 chromosomes for grain length, grain width, length to width ratio, degree of endosperm with chalkiness, rice elongation difference, and elongation index. Thirteen regions containing QTL stably expressed in multiple environments and/or exerting pleiotropic effects on multiple traits were detected. By gene-based association analysis and haplotype analysis, 46 candidate genes, including five cloned genes, and 49 favorable alleles were identified for these 13 QTL. The effect of the candidate gene Wx on rice elongation difference was validated by a transgenic strategy. These results shed light on the genetic bases of appearance quality and cooked rice elongation and provide gene resources for improving rice quality by molecular breeding.

Genome-wide association study and genomic prediction of Fusarium ear rot resistance in tropical maize germplasm

Fusarium ear rot(FER)is a destructive maize fungal disease worldwide.In this study,three tropical maize populations consisting of 874 inbred lines were used to perform genomewide association study(GWAS)and genomic prediction(GP)analyses of FER resistance.Broad phenotypic variation and high heritability for FER were observed,although it was highly influenced by large genotype-by-environment interactions.In the 874 inbred lines,GWAS with general linear model(GLM)identified 3034 single-nucleotide polymorphisms(SNPs)significantly associated with FER resistance at the P-value threshold of 1×10^(-5),the average phenotypic variation explained(PVE)by these associations was 3%with a range from 2.33%to 6.92%,and 49 of these associations had PVE values greater than 5%.The GWAS analysis with mixed linear model(MLM)identified 19 significantly associated SNPs at the P-value threshold of 1×10^(-4),the average PVE of these associations was 1.60%with a range from 1.39%to 2.04%.Within each of the three populations,the number of significantly associated SNPs identified by GLM and MLM ranged from 25 to 41,and from 5 to 22,respectively.Overlapping SNP associations across populations were rare.A few stable genomic regions conferring FER resistance were identified,which located in bins 3.04/05,7.02/04,9.00/01,9.04,9.06/07,and 10.03/04.The genomic regions in bins 9.00/01 and 9.04 are new.GP produced moderate accuracies with genome-wide markers,and relatively high accuracies with SNP associations detected from GWAS.Moderate prediction accuracies were observed when the training and validation sets were closely related.These results implied that FER resistance in maize is controlled by minor QTL with small effects,and highly influenced by the genetic background of the populations studied.Genomic selection(GS)by incorporating SNP associations detected from GWAS is a promising tool for improving FER resistance in maize.

Genetic dissection of the grain-filling rate and related traits through linkage analysis and genome-wide association study in bread wheat

Wheat grain yield is generally sink-limited during grain filling.The grain-filling rate(GFR)plays a vital role but is poorly studied due to the difficulty of phenotype surveys.This study explored the grain-filling traits in a recombinant inbred population and wheat collection using two highly saturated genetic maps for linkage analysis and genome-wide association study(GWAS).Seventeen stable additive quantitative trait loci(QTLs)were identified on chromosomes 1B,4B,and 5A.The linkage interval between IWB19555 and IWB56078 showed pleiotropic effects on GFR_(1),GFR_(max),kernel length(KL),kernel width(KW),kernel thickness(KT),and thousand kernel weight(TKW),with the phenotypic variation explained(PVE)ranging from 13.38%(KW)to 33.69%(TKW).198 significant marker-trait associations(MTAs)were distributed across most chromosomes except for 3D and 4D.The major associated sites for GFR included IWB44469(11.27%),IWB8156(12.56%)and IWB24812(14.46%).Linkage analysis suggested that IWB35850,identified through GWAS,was located in approximately the same region as QGFR_(max)2B.3-11,where two high-confidence candidate genes were present.Two important grain weight(GW)-related QTLs colocalized with grain-filling QTLs.The findings contribute to understanding the genetic architecture of the GFR and provide a basic approach to predict candidate genes for grain yield trait QTLs.

Genome-wide association study and metabolic pathway prediction of barrenness in maize as a response to high planting density

Increasing the planting density is one way to enhance grain production in maize.However,high planting density brings about growth and developmental defects such as barrenness,which is the major factor limiting grain yield.In this study,the barrenness was characterized in an association panel comprising 280 inbred lines under normal(67500 plants ha–1,ND)and high(120000 plants ha–1,HD)planting densities in 2017 and 2018.The population was genotyped using 776254 single nucleotide polymorphism(SNP)markers with criteria of minor allele frequency>5%and<20%missing data.A genome-wide association study(GWAS)was conducted for barrenness under ND and HD,as well as the barrenness ratio(HD/ND),by applying a Mixed Linear Model that controls both population structure and relative kinship(Q+K).In total,20 SNPs located in nine genes were significantly(P<6.44×10–8)associated with barrenness under the different planting densities.Among them,seven SNPs for barrenness at ND and HD were located in two genes,four of which were common under both ND and HD.In addition,13 SNPs for the barrenness ratio were located in seven genes.A complementary pathway analysis indicated that the metabolic pathways of amino acids,such as glutamate and arginine,and the mitogen-activated protein kinase(MAPK)signaling pathway might play important roles in tolerance to high planting density.These results provide insights into the genetic basis of high planting density tolerance and will facilitate high yield maize breeding.

A genome-wide association study identifies novel genetic loci that modify pharmacokinetic-pharmacodynamic responses to clopidogrel

OBJECTIVE Genetic variants in the pharmacokinetic(PK)mechanism are the main underlying factors that modify the antiplatelet efficacy of clopidogrel.Hence,joint analysis of genetic variants that modify pharmacodynamic(PD)and PK responses to clopidogrel should be effective for identifying the genetic variants affecting the antiplatelet response to the drug.METHODS A genome-wide association study was conducted to identify new genetic loci that modify PD responses to clopidogrel and its active metabolite H4 in 115 Chinese patients with coronary heart disease(CHD).RESULTS We identified novel variants in two transporter genes(rs12456693 in SLC14A2 and rs2487032 in ABCA1)and in N6AMT1(rs2254638)associated with clopidogrel-treated P2Y12reaction unit(PRU)and plasma H4 concentration.The associations between these single nucleotide polymorphisms(SNPs)and PK parameters of clopidogrel and H4 were observed in 31 additional CHD patients(P<0.05).The new variants,together with CYP2C19*2 and clinical factors,dramatically improved the predictability of PRU variability to 37.7%compared with the published value of approximately 20%.The function of these SNPs on the activation of clopidogrel was validated in 32 liver S9 fractions,and the N6AMT1 rs2254638 T variant was found to be associated with decreased formation of H4(P=0.0386).Meanwhile,N6AMT1 rs2254638 was further identified to exert a marginal risk effect for MACE in an independent CHD patient cohort(OR:1.428,95%CI:0.978-2.086,P=0.0653,FDR=0.4726).In conclusion,we systematically identified new genetic variants as risk factors for the reduced efficacy of clopidogrel.CONCLUSION Our study findings enhanced the understanding of the absorption and metabolic mechanisms that influence PD responses to clopidogrel treatment.

Genetic Dissection of Grain Size Traits Through Genome-Wide Association Study Based on Genic Markers in Rice

Grain size plays a significant role in rice,starting from affecting yield to consumer preference,which is the driving force for deep investigation and improvement of grain size characters.Quantitative inheritance makes these traits complex to breed on account of several alleles contributing to the complete trait expression.We employed genome-wide association study in an association panel of 88 rice genotypes using 142 new candidate gene based SSR(cgSSR)markers,derived from yield-related candidate genes,with the efficient mixed-model association coupled mixed linear model for dissecting complete genetic control of grain size traits.A total of 10 significant associations were identified for four grain size-related characters(grain weight,grain length,grain width,and length-width ratio).Among the identified associations,seven marker trait associations explain more than 10%of the phenotypic variation,indicating major putative QTLs for respective traits.The allelic variations at genes OsBC1L4,SHO1 and OsD2 showed association between 1000-grain weight and grain width,1000-grain weight and grain length,and grain width and length-width ratio,respectively.The cgSSR markers,associated with corresponding traits,can be utilized for direct allelic selection,while other significantly associated cgSSRs may be utilized for allelic accumulation in the breeding programs or grain size improvement.The new cgSSR markers associated with grain size related characters have a significant impact on practical plant breeding to increase the number of causative alleles for these traits through marker aided rice breeding programs.

Genome-Wide Association Study of Nitrogen Use Efficiency and Agronomic Traits in Upland Rice

Genome-wide association study(GWAS)was performed for 16 agronomic traits including nitrogen use efficiency(NUE)and yield-related components using a panel of 190 mainly japonica rice varieties and a set of 38390 single nucleotide polymorphism(SNP)markers.This panel was evaluated under rainfed upland conditions in Madagascar in two consecutive cropping seasons with two contrasted nitrogen input levels.Using another set of five grain traits,we identified previously known genes(GW5,GS3,Awn1 and Glabrous1),thus validating the pertinence and accuracy of our datasets for GWAS.A total of 369 significant associations were detected between SNPs and agronomic traits,gathered into 46 distinct haplotype groups and 28 isolated markers.Few association signals were identified for the complex quantitative trait NUE,however,larger number of quantitative trait loci(QTLs)were detected for its component traits,with 10 and 2 association signals for nitrogen utilization efficiency and nitrogen uptake efficiency,respectively.Several detected association signals co-localized with genes involved in nitrogen transport or nitrogen remobilization within 100 kb.The present study thus confirmed the potential of GWAS to identify candidate genes and new loci associated with agronomic traits.However,because of the quantitative and complex nature of NUE-related traits,GWAS might have not captured a large number of QTLs with limited effects.

Integrated linkage mapping and genome-wide association study to dissect the genetic basis of zinc deficiency tolerance in maize at seedling stage

Zinc(Zn)deficiency is the most widespread micronutrient deficiency,affecting yield and quality of crops worldwide.Identifying genes associated with Zn-deficiency tolerance in maize is a basis for elucidating its genetic mechanism.A K22×CI7 recombinant inbred population consisting of 210 lines and an association panel of 508 lines were used to identify genetic loci influencing Zn-deficiency tolerance.Under-Zn and-Zn/CK conditions,15 quantitative trait loci(QTL)were detected,each explaining 5.7%-12.6%of phenotypic variation.Sixty-one significant single-nucleotide polymorphisms(SNPs)were identified at P<10^(-5)by genome-wide association study(GWAS),accounting for 5%-14%of phenotypic variation.Among respectively 198 and 183 candidate genes identified within the QTL regions and the 100-kb regions flanking these significant SNPs,12 were associated with Zn-deficiency tolerance.Among these candidate genes,four genes associated with hormone signaling in response to Zn-deficiency stress were co-localized with QTL or SNPs,including the genes involved in the auxin(ZmARF7),and ethylene(ZmETR5,ZmESR14,and ZmEIN2)signaling pathways.Three candidate genes were identified as being responsible for Zn transport,including ZmNAS3 detected by GWAS,ZmVIT and ZmYSL11 detected by QTL mapping.Expression of ZmYSL11 was up-regulated in Zn-deficient shoots.Four candidate genes that displayed different expression patterns in response to Zn deficiency were detected in the regions overlapping peak GWAS signals,and the haplotypes for each candidate gene were further analyzed.

Genome-wide association study reveals novel loci for adult type 1 diabetes in a 5-year nested case-control study

BACKGROUND Type 1 diabetes(T1D)is a severe and prevalent metabolic disease.Due to its high heredity,an increasing number of genome-wide association studies have been performed,most of which were from hospital-based case-control studies with a relatively small sample size.The association of single nucleotide polymorphisms(SNPs)and T1D has been less studied and is less understood in natural cohorts.AIM To investigate the significant variants of T1D,which could be potential biomarkers for T1D prediction or even therapy.METHODS A genome-wide association study(GWAS)of adult T1D was performed in a nested case-control study(785 cases vs 804 controls)from a larger 5-year cohort study in Suzhou,China.Potential harmful or protective SNPs were evaluated for T1D.Subsequent expression and splicing quantitative trait loci(eQTL and sQTL)analyses were carried out to identify target genes modulated by these SNPs.RESULTS A harmful SNP for T1D,rs3117017[odds ratio(OR)=3.202,95%confidence interval(CI):2.296-4.466,P=9.33×10-4]and three protective SNPs rs55846421(0.113,0.081-0.156,1.76×10-9),rs75836320(0.283,0.205-0.392,1.07×10-4),rs362071(0.568,0.495-0.651,1.66×10-4)were identified.Twenty-two genes were further identified as potential candidates for T1D onset.CONCLUSION We identified a potential genetic basis of T1D,both protective and harmful,using a GWAS in a larger nested case-control study of a Chinese population.

A Hybrid Algorithm for Disease Association Study

Single nucleotide polymorphisms (SNPs), which are the most common form of DNA variations, have great potential as a medical diagnostic tool. However, compared to the number of SNPs involved, the available training data sets generally have a fairly small sample size, which is a main challenge to traditional data analysis methods. This paper proposed an improved univariate marginal distribution algorithm (UMDA) named multi-population UMDA (MPUMDA) for disease association study. In order to illustrate the effectiveness of our algorithm, we compared it with some current known methods, and the results showed that our method is potentially interesting as an alter-native tool in disease association study.

Pathway-based analysis of genome-wide association study of circadian phenotypes

Sleepiness affects normal social life, which attracts more and more attention. Circadian phenotypes contribute to obvious individual differences in susceptibility to sleepiness. We aimed to identify candidate single nucleotide polymorphisms(SNPs) which may cause circadian phenotypes, elucidate the potential mechanisms, and generate corresponding SNP-gene-pathways. A genome-wide association studies(GWAS) dataset of circadian phenotypes was utilized in the study. Then, the Identify Candidate Causal SNPs and Pathways analysis was employed to the GWAS dataset after quality control filters. Furthermore, genotype-phenotype association analysis was performed with HapMap database. Four SNPs in three different genes were determined to correlate with usual weekday bedtime,totally providing seven hypothetical mechanisms. Eleven SNPs in six genes were identified to correlate with usual weekday sleep duration, which provided six hypothetical pathways. Our results demonstrated that fifteen candidate SNPs in eight genes played vital roles in six hypothetical pathways implicated in usual weekday bedtime and six potential pathways involved in usual weekday sleep duration.

Genome-Wide Association Study of Cooked Rice Textural Attributes and Starch Physicochemical Properties in indica Rice

Rice cooking and eating qualities(CEQ)are mainly determined by cooked rice textural parameters and starch physicochemical properties.However,the genetic bases of grain texture and starch properties in rice have not been fully understood.We conducted a genome-wide association study for apparent amylose content(AAC),starch pasting viscosities,and cooked rice textural parameters using 279 indica rice accessions from the 3000 Rice Genome Project.We identified 26 QTLs in the whole population and detected single nucleotide polymorphisms(SNPs)with the lowest P-value at the Waxy(Wx)locus for all traits except pasting temperature and cohesiveness.Additionally,we detected significant SNPs at the SUBSTANDARD STARCH GRAIN6(SSG6)locus for AAC,setback(SB),hardness,adhesiveness,chewiness(CHEW),gumminess(GUM),and resilience.We subsequently divided the population using a SNP adjacent to the Waxy locus,and identified 23 QTLs and 12 QTLs in two sub-panels,WxT and WxA,respectively.In these sub-panels,SSG6 was also identified to be associated with pasting parameters,including peak viscosity,hot paste viscosity,cold paste viscosity,and consistency viscosity.Furthermore,a candidate gene encoding monosaccharide transporter 5(OsMST5)was identified to be associated with AAC,breakdown,SB,CHEW,and GUM.In total,39 QTLs were co-localized with known genes or previously reported QTLs.These identified genes and QTLs provide valuable information for genetic manipulation to improve rice CEQ.

Genome-wide association with transcriptomics reveals a shade-tolerance gene network in soybean

Shade tolerance is essential for soybeans in inter/relay cropping systems.A genome-wide association study(GWAS)integrated with transcriptome sequencing was performed to identify genes and construct a genetic network governing the trait in a set of recombinant inbred lines derived from two soybean parents with contrasting shade tolerance.An improved GWAS procedure,restricted two-stage multi-locus genome-wide association study based on gene/allele sequence markers(GASM-RTM-GWAS),identified 140 genes and their alleles associated with shade-tolerance index(STI),146 with relative pith cell length(RCL),and nine with both.Annotation of these genes by biological categories allowed the construction of a protein–protein interaction network by 187 genes,of which half were differentially expressed under shading and non-shading conditions as well as at different growth stages.From the identified genes,three ones jointly identified for both traits by both GWAS and transcriptome and two genes with maximum links were chosen as beginners for entrance into the network.Altogether,both STI and RCL gene systems worked for shade-tolerance with genes interacted each other,this confirmed that shadetolerance is regulated by more than single group of interacted genes,involving multiple biological functions as a gene network.

Causal associations between inflammatory bowel disease and anxiety:A bidirectional Mendelian randomization study

BACKGROUND Anxiety is common in patients with inflammatory bowel disease(IBD),including those with ulcerative colitis(UC)and Crohn’s disease(CD);however,the causal relationship between IBD and anxiety remains unknown.AIM To investigate the causal relationship between IBD and anxiety by using bidirectional Mendelian randomization analysis.METHODS Single nucleotide polymorphisms retrieved from genome-wide association studies(GWAS)of the European population were identified as genetic instrument variants.GWAS statistics for individuals with UC(6968 patients and 20464 controls;adults)and CD(5956 patients and 14927 controls;adults)were obtained from the International IBD Genetics Consortium.GWAS statistics for individuals with anxiety were obtained from the Psychiatric Genomics Consortium(2565 patients and 14745 controls;adults)and FinnGen project(20992 patients and 197800 controls;adults),respectively.Inverse-variance weighted was applied to assess the causal relationship,and the results were strengthened by heterogeneity,pleiotropy and leave-one-out analyses.RESULTS Genetic susceptibility to UC was associated with an increased risk of anxiety[odds ratio:1.071(95%confidence interval:1.009-1.135),P=0.023],while genetic susceptibility to CD was not associated with anxiety.Genetic susceptibility to anxiety was not associated with UC or CD.No heterogeneity or pleiotropy was observed,and the leave-one-out analysis excluded the potential influence of a particular variant.CONCLUSION This study revealed that genetic susceptibility to UC was significantly associated with anxiety and highlighted the importance of early screening for anxiety in patients with UC.

Genome-wide and candidate gene association studies identify BnPAP17 as conferring the utilization of organic phosphorus in oilseed rape

Phosphorus(P)is essential for living plants,and P deficiency is one of the key factors limiting the yield in rapeseed production worldwide.As the most important organ for plants,root morphology traits(RMTs)play a key role in P absorption.To investigate the genetic variability of RMT under low P availability,we dissected the genetic structure of RMTs by genome-wide association studies(GWAS),linkage mapping and candidate gene association studies(CGAS).A total of 52 suggestive loci were associated with RMTs under P stress conditions in 405 oilseed rape accessions.The purple acid phosphatase gene BnPAP17 was found to control the lateral root number(LRN)and root dry weight(RDW)under low P stress.The expression of BnPAP17 was increased in shoot tissue in P-efficient cultivars compared to root tissue and P-inefficient cultivars in response to low P stress.Moreover,the haplotype of BnPAP17^(Hap3)was detected for the selective breeding of P efficiency in oilseed rape.Over-expression of the BnPAP17^(Hap3)could promote the shoot and root growth with enhanced tolerance to low P stress and organic phosphorus(Po)utilization in oilseed rape.Collectively,these findings increase our understanding of the mechanisms underlying BnPAP17-mediated low P stress tolerance in oilseed rape.

Multi-locus association study of schizophrenia susceptibility genes with a posterior probability method

Schizophrenia is a serious neuropsychiatric illness affecting about 1% of the world’s population. It is considered a complex inheritance disorder. A number of genes are involved in combination in the etiology of the disorder. Evidence implicates the altered dopaminergic trans- mission in schizophrenia. In the present study, in order to identify susceptibility genes for schizophrenia in dopaminergic metabolism, we analyzed 59 single nucleotide polymorphisms (SNPs) in 24 genes of the dopaminergic pathway among 82 unrelated patients with schizophre- nia and 108 matched normal controls. Considering that traditional single-locus association stud- ies ignore the multigenic nature of complex diseases and do not take into account possible in- teractions between susceptibility genes, we proposed a multi-locus analysis method, using the posterior probability of morbidity as a measure of absolute disease risk for a multi-locus genotype combination, and developed an algorithm based on perturbation and average to detect the sus- ceptibility multi-locus genotype combinations, as well as to repress noise and avoid false positive results at our best. A three-locus SNP genotype combination involved in the interactions of COMT and ALDH3B1 genes was detected to be significantly susceptible to schizophrenia.

A genome-wide association study of facial morphology identifies novel genetic loci in Han Chinese

The human face is a heritable surface with many complex sensory organs.In recent years,many genetic loci associated with facial features have been reported in different populations,yet there is a lack of studies on the Han Chinese population.Here,we report a genome-wide association study of 3D normal human faces of 2,659 Han Chinese with autosegment phenotypes of facial morphology.We identify singlenucleotide polymorphisms(SNPs)encompassing four genomic regions showing significant associations with different facial regions,including SNPs in DENND1 B associated with the chin,SNPs among PISRT1 associated with eyes,SNPs between DCHS2 and SFRP2 associated with the nose,and SNPs in VPS13 B associated with the nose.We replicate 24 SNPs from previously reported genetic loci in different populations,whose candidate genes are DCHS2,SUPT3 H,HOXD1,SOX9,PAX3,and EDAR.These results provide a more comprehensive understanding of the genetic basis of variation in human facial morphology.

Genome-wide association study of the backfat thickness trait in two pig populations

Backfat thickness is a good predictor of carcass lean content,an economically important trait,and a main breeding target in pig improvement.In this study,the candidate genes and genomic regions associated with the tenth rib backfat thickness trait were identified in two independent pig populations,using a genome-wide association study of porcine 60K SNP genotype data applying the compressed mixed linear model(CMLM)statistical method.For each population,30 most significant single-nucleotide polymorphisms(SNPs)were selected and SNP annotation implemented using Sus scrofa Build 10.2.In the first population,25 significant SNPs were distributed on seven chromosomes,and SNPs on SSC1 and SSC7 showed great significance for fat deposition.The most significant SNP(ALGA0006623)was located on SSC1,upstream of the MC4R gene.In the second population,27 significant SNPs were recognized by annotation,and 12 SNPs on SSC12 were related to fat deposition.Two haplotype blocks,M1GA0016251-MARC0075799 and ALGA0065251-MARC0014203-M1GA0016298-ALGA0065308,were detected in significant regions where the PIPNC1 and GH1 genes were identified as contributing to fat metabolism.The results indicated that genetic mechanism regulating backfat thickness is complex,and that genome-wide associations can be affected by populations with different genetic backgrounds.