Genome-wide association study reveals the genetic determinism of growth traits in a Gushi-Anka F2 chicken population

Chicken growth traits are economically important, but the relevant genetic mechanisms have not yet been elucidated. Herein, we performed a genome-wide association study to identify the variants associated with growth traits. In total, 860 chickens from a Gushi-Anka F₂ resource population were phenotyped for 68 growth and carcass traits, and 768 samples were genotyped based on the genotyping-by-sequencing (GBS) method. Finally, 734 chickens and 321,314 SNPs remained after quality control and removal of the sex chromosomes, and these data were used to carry out a GWAS analysis. A total of 470 significant single-nucleotide polymorphisms (SNPs) for 43 of the 68 traits were detected and mapped on chromosomes (Chr) 1–6, -9, -10, -16, -18, -23, and -27. Of these, the significant SNPs in Chr1, -4, and -27 were found to be associated with more than 10 traits. Multiple traits shared significant SNPs, indicating that the same mutation in the region might have a large effect on multiple growth or carcass traits. Haplotype analysis revealed that SNPs within the candidate region of Chr1 presented a mosaic pattern. The significant SNPs and pathway enrichment analysis revealed that the MLNR, MED4, CAB39L, LDB2, and IGF2BP1 genes could be putative candidate genes for growth and carcass traits. The findings of this study improve our understanding of the genetic mechanisms regulating chicken growth and carcass traits and provide a theoretical basis for chicken breeding programs.Subject terms: Genome-wide association studies, Genetic linkage study, Development, DNA sequencing, Animal breeding     

Genome‐wide association study reveals novel variants for growth and egg traits in Dongxiang blue‐shelled and White Leghorn chickens

This study was designed to investigate the genetic basis of growth and egg traits in Dongxiang blue‐shelled chickens and White Leghorn chickens. In this study, we employed a reduced representation sequencing approach called genotyping by genome reducing and sequencing to detect genome‐wide SNPs in 252 Dongxiang blue‐shelled chickens and 252 White Leghorn chickens. The Dongxiang blue‐shelled chicken breed has many specific traits and is characterized by blue‐shelled eggs, black plumage, black skin, black bone and black organs. The White Leghorn chicken is an egg‐type breed with high productivity. As multibreed genome‐wide association studies (GWASs) can improve precision due to less linkage disequilibrium across breeds, a multibreed GWAS was performed with 156 575 SNPs to identify the associated variants underlying growth and egg traits within the two chicken breeds. The analysis revealed 32 SNPs exhibiting a significant genome‐wide association with growth and egg traits. Some of the significant SNPs are located in genes that are known to impact growth and egg traits, but nearly half of the significant SNPs are located in genes with unclear functions in chickens. To our knowledge, this is the first multibreed genome‐wide report for the genetics of growth and egg traits in the Dongxiang blue‐shelled and White Leghorn chickens.     

京海黄鸡体重性状全基因组关联分析

体重性状是肉鸡重要的经济性状。为了寻找可用于京海黄鸡体重性状遗传改良的分子标记及候选基因,本文以400只京海黄鸡核心群母鸡为基础,测定了0~14周龄体重,利用简化基因组测序技术(Specific-locus amplified fragment sequencing, SLAF-seq)对京海黄鸡体重性状进行全基因组关联研究(Genome-wide association stndy, GWAS),筛选与京海黄鸡体重性状相关的SNPs位点。结果共检测到100个与京海黄鸡体重相关的SNPs位点,其中15个位点效应达到全基因组显著水平(P<1.87E-06),85个位点效应达到全基因组潜在显著水平(P<3.73E-05)。通过筛选每个显著SNP周围1 Mb区域内的基因,共找到9个可能的候选基因,其中FAM124A(Family with sequence similarity 124A)、QDPR(Quinoid dihydropteridine reductase)、WDR1(WD repeat domain 1)和SLC2A9(Solute carrier family 2 (facilitated glucose transporter), member 9) 4个基因可能是影响体重性状的重要候选基因。同时还发现,4号染色体75.6~80.7 Mb区域集中了大部分与京海黄鸡中后期体重性状显著相关的SNPs位点,该区域可能是影响京海黄鸡中后期生长体重的重要候选区域。     

Genome-wide association study of shank length and diameter at different developmental stages in chicken F2 resource population

In order to find SNPs and genes affecting shank traits, we performed a GWAS in a chicken F2 population of eight half-sib families from five hatches derived from reciprocal crosses between an Arian fast-growing line and an Urmia indigenous slow-growing chicken. A total of 308 birds were genotyped using a 60K chicken SNP chip. Shank traits including shank length and diameter were measured weekly from birth to 12 weeks of age. A generalized linear model and a compressed mixed linear model (CMLM) were applied to achieve the significant regions. The value of the average genomic inflation factor (λ statistic) of the CMLM model (0.99) indicated that the CMLM was more effective than the generalized linear model in controlling the population structure. The genes surrounding significant SNPs and their biological functions were identified from NCBI, Ensembl and UniProt databases. The results indicated that 12 SNPs at 12 different ages passed the LD-adjusted 5% Bonferroni significant threshold. Two SNPs were significant for shank length and nine SNPs were significant for shank diameter. The significant SNPs were located near to or inside 11 candidate genes. The results showed that a number of significant SNPs in the middle ages were higher than the rest. The MXRA8 gene was related to the significant SNP at week 1 that promotes proliferation of growth plate chondrocytes. A unique SNP of Gga_rs16689511 located on chicken Z chromosome within the LOC101747628 gene was related to shank length at three different ages of birds (weeks 8, 9 and 11). The significant SNPs for shank diameter were found at weeks 4 and 7 (four and five SNPs respectively). The identifications of SNPs and genes here could contribute to a better understanding of the genetic control of shank traits in chicken.     

A Genome-Wide Association Study of the Maize Hypersensitive Defense Response Identifies Genes That Cluster in Related Pathways

Much remains unknown of molecular events controlling the plant hypersensitive defense response (HR), a rapid localized cell death that limits pathogen spread and is mediated by resistance (R-) genes. Genetic control of the HR is hard to quantify due to its microscopic and rapid nature. Natural modifiers of the ectopic HR phenotype induced by an aberrant auto-active R-gene (Rp1-D21), were mapped in a population of 3,381 recombinant inbred lines from the maize nested association mapping population. Joint linkage analysis was conducted to identify 32 additive but no epistatic quantitative trait loci (QTL) using a linkage map based on more than 7000 single nucleotide polymorphisms (SNPs). Genome-wide association (GWA) analysis of 26.5 million SNPs was conducted after adjusting for background QTL. GWA identified associated SNPs that colocalized with 44 candidate genes. Thirty-six of these genes colocalized within 23 of the 32 QTL identified by joint linkage analysis. The candidate genes included genes predicted to be in involved programmed cell death, defense response, ubiquitination, redox homeostasis, autophagy, calcium signalling, lignin biosynthesis and cell wall modification. Twelve of the candidate genes showed significant differential expression between isogenic lines differing for the presence of Rp1-D21. Low but significant correlations between HR-related traits and several previously-measured disease resistance traits suggested that the genetic control of these traits was substantially, though not entirely, independent. This study provides the first system-wide analysis of natural variation that modulates the HR response in plants.     

Genome-Wide Association Study Identifies Loci and Candidate Genes for Body Composition and Meat Quality Traits in Beijing-You Chickens

Body composition and meat quality traits are important economic traits of chickens. The development of high-throughput genotyping platforms and relevant statistical methods have enabled genome-wide association studies in chickens. In order to identify molecular markers and candidate genes associated with body composition and meat quality traits, genome-wide association studies were conducted using the Illumina 60 K SNP Beadchip to genotype 724 Beijing-You chickens. For each bird, a total of 16 traits were measured, including carcass weight (CW), eviscerated weight (EW), dressing percentage, breast muscle weight (BrW) and percentage (BrP), thigh muscle weight and percentage, abdominal fat weight and percentage, dry matter and intramuscular fat contents of breast and thigh muscle, ultimate pH, and shear force of the pectoralis major muscle at 100 d of age. The SNPs that were significantly associated with the phenotypic traits were identified using both simple (GLM) and compressed mixed linear (MLM) models. For nine of ten body composition traits studied, SNPs showing genome wide significance (P<2.59E−6) have been identified. A consistent region on chicken (Gallus gallus) chromosome 4 (GGA4), including seven significant SNPs and four candidate genes (LCORL, LAP3, LDB2, TAPT1), were found to be associated with CW and EW. Another 0.65 Mb region on GGA3 for BrW and BrP was identified. After measuring the mRNA content in beast muscle for five genes located in this region, the changes in GJA1 expression were found to be consistent with that of breast muscle weight across development. It is highly possible that GJA1 is a functional gene for breast muscle development in chickens. For meat quality traits, several SNPs reaching suggestive association were identified and possible candidate genes with their functions were discussed.     

Major histocompatibility complex (MHC) related cDNA probes associated with antibody response in meat-type chickens

The major histocompatibility complex (MHC) region was examined as a set of candidate genes for association between DNA markers and antibody response. Intercross F₂ families of chickens were generated from a cross between high (HC) and low (LC) Escherichia coli _i antibody lines. Restriction fragment length polymorphism (RFLP) analysis was conducted by using three MHC-related cDNA probes: chicken MHC class IV ( B-G ), chicken MHC class I ( B-F ), and human MHC-linked Tap2 . Association between RFLP bands and three antibody response traits ( E. coli , sheep red blood cells and Newcastle disease virus) were determined by two methods: by statistically analyzing each band separately and also by analyzing all bands obtained from the three probes by using multiple regression analysis to account for the multiple comparisons. The MHC class IV probe was the highest in polymorphisms but had the lowest number of bands associated with antibody response. The MHC class I probe yielded 15 polymorphic bands of which four exhibited association with antibody response traits. The Tap2 probe yielded 20 different RFLP bands of which five were associated with antibody production. Some Tap2 bands were associated with multiple antibody response traits. The multiband analysis of the three probes' bands revealed more significant effects than the analysis of each band separately. This study illustrates the efficacy of using multiple MHC region probes as candidate markers for quantitative trait loci (QTLs) controlling antibody response in chickens.     

High diversity of the chicken growth hormone gene and effects on growth and carcass traits

The chicken growth hormone (cGH) gene plays a crucial role in controlling growth and metabolism, leading to potential correlations between cGH polymorphisms and economic traits. In this study, DNA from four divergent chicken breeds were screened for single nucleotide polymorphisms (SNPs) in the cGH gene using denaturing high-performance liquid chromatography and sequencing. A total of 46 SNPs were identified, of which 4 were in the 5' untranslated region, 1 in the 3' untranslated region, 5 in exons (two of which are nonsynonymous), with the remaining 36 in introns. The nucleotide diversity in the cGH gene ( theta = 2.7 x 10(-3)) was higher than that reported for other chicken genes, even within the same breeds. The associations of five of these SNPs and their haplotypes with chicken growth and carcass traits were determined using polymerase chain reaction-restriction fragment length polymorphism analysis in a F2 resource population cross of two of the four chicken breeds (White Recessive Rock and Xinghua). This analysis shows that, among other correlations, G+1705A was significantly associated with body weight at all ages measured, shank length at three of four ages measured, and average daily gain within weeks 0 to 4. Thus, this cGH polymorphism, or another polymorphism that is in linkage disequilibrium with G+1705A, appears to correspond to a significant growth-related quantitative trait locus difference between the two breeds used to construct the resource population.     

Genetic variants in telomere‐maintenance genes are associated with ovarian cancer risk and outcome

Most ovarian cancer patients present at an advanced stage with poor prognosis. Telomeres play a critical role in protecting chromosomes stability. The associations of genetic variants in telomere maintenance genes and ovarian cancer risk and outcome are unclear. We genotyped 137 single nucleotide polymorphisms (SNPs) in telomere‐maintenance genes in 417 ovarian cancer cases and 417 matched healthy controls to evaluate their associations with cancer risk, survival and therapeutic response. False discovery rate Q‐value was calculated to account for multiple testing. Eleven SNPs from two genes showed nominally significant associations with the risks of ovarian cancer. The most significant SNP was TEP1: rs2228026 with participants carrying at least one variant allele exhibiting a 3.28‐fold (95% CI: 1.72‐6.29; P < 0.001, Q = 0.028) increased ovarian cancer risk, which remained significant after multiple testing adjusting. There was also suggested evidence for the associations of SNPs with outcome, although none of the associations had a Q < 0.05. Seven SNPs from two genes showed associations with ovarian cancer survival (P < 0.05). The strongest association was found in TNKS gene (rs10093972, hazard ratio = 1.88; 95% CI: 1.20‐2.92; P = 0.006, Q = 0.076). Five SNPs from four genes showed suggestive associations with therapeutic response (P < 0.05). In a survival tree analysis, TEP1:rs10143407 was the primary factor contributing to overall survival. Unfavourable genotype analysis showed a cumulative effect of significant SNPs on ovarian cancer risk, survival and therapeutic response. Genetic variations in telomere‐maintenance genes may be associated with ovarian cancer risk and outcome.     

丙酸通路基因多态性与猪肉质及胴体性状的关联分析

为了挖掘新的猪肉品质及胴体性状的候选基因,揭示猪肉质及胴体性状的遗传机制,文章将丙酸代谢通路作为候选通路,将通路内基因与猪肉质及胴体性状进行关联分析。实验采用37头三元杂交商品猪作为研究对象,首次针对丙酸通路中7个基因的36个SNP位点利用SNaPshot方法进行基因分型,分别用最小二乘模型及MB-MDR模型与肉质及胴体性状进行关联分析。结果发现,基因PCCB、MUT、MCEE及ACSS2上的4个SNP位点分别与肌内脂肪含量、背膘厚等性状显著相关(P<0.05),ACSS2与猪脂肪含量显著相关;MCEE及MUT与猪的背膘厚显著相关;PCCB基因与脂重显著相关。通过MB-MDR方法检测到多个SNP位点具有互作效应,并与背膘厚、水分含量、脂肪含量显著相关(P<0.05)。另外,丙酸代谢通路中的基因间的互作效应对猪肉品质有显著影响。     

Single step genome-wide association studies based on genotyping by sequence data reveals novel loci for the litter traits of domestic pigs

In this study, data genotyping by sequence (GBS) was used to perform single step GWAS (ssGWAS) to identify SNPs associated with the litter traits in domestic pigs and search for candidate genes in the region of significant SNPs. After quality control, 167,355 high-quality SNPs from 532 pigs were obtained. Phenotypic traits on 2112 gilt litters from 532 pigs were recorded including total number born (TNB), number born alive (NBA), and litter weight born alive (LWB). A single-step genomic BLUP approach (ssGBLUP) was used to implement the genome-wide association analysis at a 5% genome-wide significance level. A total of 8, 23 and 20 significant SNPs were associated with TNB, NBA, and LWB, respectively, and these significant SNPs accounted for 62.78%, 79.75%, and 58.79% of genetic variance. Furthermore, 1 (SSC14: 16314857), 4 (SSC1: 81986236, SSC1: 66599775, SSC1: 161999013, and SSC1: 267883107), and 5 (SSC9: 29030061, SSC2: 32368561, SSC5: 110375350, SSC13: 45619882 and SSC13: 45647829) significant SNPs for TNB, NBA, and LWB were inferred to be novel loci. At SSC1, the AIM1 and FOXO3 genes were found to be associated with NBA; these genes increase ovarian reproductive capacity and follicle number and decrease gonadotropin levels. The genes SLC36A4 and INTU are involved in cell growth, cytogenesis and development were found to be associated with LWB. These significant SNPs can be used as an indication for regions in the Sus scrofa genome for variability in litter traits, but further studies are expected to confirm causative mutations.     

ApoB与UCP基因间上位效应对鸡腹脂性状影响的遗传学分析

已有研究表明, 上位效应在畜禽重要复杂经济性状的表型形成过程中发挥重要作用。文章选取对肉鸡7 周龄腹脂率有显著影响的载脂蛋白B(Apoliprotein B, ApoB)基因T123G位点与解耦连蛋白(Uncoupling protein, UCP)基因C1197A位点, 在东北农业大学肉鸡高、低腹脂双向选择品系8、9、10世代内检测这两个SNPs的多态性并利用Natural and Orthogonal InterActions (NOIA)模型分析二者之间的上位效应对7周龄腹脂率的影响。结果表明, 在高脂系内这两个位点之间存在着对7周龄腹脂率有显著影响的上位效应组分(P<0.05), 并且在连续多世代选育过程中仍能持续稳定存在; 同时, 在低脂系中二者之间各上位效应组分对腹脂率均无影响(P>0.05)。此结果意味着, 至少在高脂系内, 腹脂率的遗传受到这两个基因间上位效应的影响; 同时提示两系间脂肪性状QTL或重要候选基因功能位点之间不同的遗传互作模式可能是引起这两个品系腹脂性状巨大表型差异的重要影响因素之一。     

Genome-wide association study for longevity in the Holstein cattle population

Longevity is one of the most important traits determining dairy cow profitability. In the last decades dairy cows suffered a lowering in the age at culling. With the aim to identify the genes involved in longevity, dates of birth, yields, dates of calving during lifespan and culling dates were collected for 946 culled cows which had been genotyped with the Bovine High Density panel. Using the GenABEL package in R, genome-wide association analysis was performed on three potential traits of longevity: (1) ‘days in production,’ (2) ‘days in herd,’ (3) number of calvings over lifespan.’ Five genome-wide significant single nucleotide polymorphisms (SNPs) associated with all three longevity traits were detected. Several consecutive SNPs identified on chromosomes 16 and 30 indicated the presence of two suggestive quantitative trait loci (QTL). The genes comprised in the QTL regions had biological functions related to fertility, reproductive disorders, heat stress and welfare of cows. These findings might contribute to improving breeding strategies to improve longevity.     

Identification of candidate genes and mutations in QTL regions for immune responses in chicken

There are two categories of immune responses – innate and adaptive immunity – both having polygenic backgrounds and a significant environmental component. In our study, adaptive immunity was represented by the specific antibody response toward keyhole limpet hemocyanin (KLH); innate immunity was represented by natural antibodies toward lipopolysaccharide (LPS) and lipoteichoic acid (LTA). Defining genetic bases of immune responses leads from defining quantitative trait loci (QTL) toward a single mutation responsible for variation in the phenotypic trait. The goal of the reported study was to define candidate genes and mutations for the immune traits of interest in chicken by performing an association study of SNPs located in candidate genes defined in QTL regions. Candidate genes and SNPs in QTL regions were selected in silico. SNP association was based on a custom SNP panel, GoldenGate genotyping assay (Illumina) and two statistical models: random mixed model and CAR score. The most significant SNP for immune response toward KLH was located in the JMJD6 gene located on GGA18. Four SNPs in candidate genes FOXJ1 (GGA18), EPHB1 (GGA9), PTGER4 (GGAZ) and PRKCB (GGA14) showed association with natural antibodies for LPS. A single SNP in ITGB4 (GGA18) was associated with natural antibodies for LTA. All associated SNPs mentioned above showed additive effects.     

Association of polymorphisms in oxidative stress genes with clinical outcomes for bladder cancer treated with Bacillus Calmette-Guérin

Genetic polymorphisms in oxidative stress pathway genes may contribute to carcinogenesis, disease recurrence, treatment response, and clinical outcomes. We applied a pathway-based approach to determine the effects of multiple single nucleotide polymorphisms (SNPs) within this pathway on clinical outcomes in non-muscle-invasive bladder cancer (NMIBC) patients treated with Bacillus Calmette-Guérin (BCG). We genotyped 276 SNPs in 38 genes and evaluated their associations with clinical outcomes in 421 NMIBC patients. Twenty-eight SNPs were associated with recurrence in the BCG-treated group (P<0.05). Six SNPs, including five in NEIL2 gene from the overall and BCG group remained significantly associated with recurrence after multiple comparison adjustments (q<0.1). Cumulative unfavorable genotype analysis showed that the risk of recurrence increased with increasing number of unfavorable genotypes. In the analysis of risk factors associated with progression to disease, rs3890995 in UNG, remained significant after adjustment for multiple comparison (q<0.1). These results support the hypothesis that genetic variations in host oxidative stress genes in NMIBC patients may affect response to therapy with BCG.     

Understanding genetic control of root system architecture in soybean: Insights into the genetic basis of lateral root number

Developing crops with better root systems is a promising strategy to ensure productivity in both optimum and stress environments. Root system architectural traits in 397 soybean accessions were characterized and a high‐density single nucleotide polymorphisms (SNPs)‐based genome‐wide association study was performed to identify the underlying genes associated with root structure. SNPs associated with root architectural traits specific to landraces and elite germplasm pools were detected. Four loci were detected in landraces for lateral root number (LRN) and distribution of root thickness in diameter Class I with a major locus on chromosome 16. This major loci was detected in the coding region of unknown protein, and subsequent analyses demonstrated that root traits are affected with mutated haplotypes of the gene. In elite germplasm pool, 3 significant SNPs in alanine‐glyoxalate aminotransferase, Leucine‐Rich Repeat receptor/No apical meristem, and unknown functional genes were found to govern multiple traits including root surface area and volume. However, no major loci were detected for LRN in elite germplasm. Nucleotide diversity analysis found evidence of selective sweeps around the landraces LRN gene. Soybean accessions with minor and mutated allelic variants of LRN gene were found to perform better in both water‐limited and optimal field conditions.     

A Multiple-SNP Approach for Genome-Wide Association Study of Milk Production Traits in Chinese Holstein Cattle

The multiple-SNP analysis has been studied by many researchers, in which the effects of multiple SNPs are simultaneously estimated and tested in a multiple linear regression. The multiple-SNP association analysis usually has higher power and lower false-positive rate for detecting causative SNP(s) than single marker analysis (SMA). Several methods have been proposed to simultaneously estimate and test multiple SNP effects. In this research, a fast method called MEML (Mixed model based Expectation-Maximization Lasso algorithm) was developed for simultaneously estimate of multiple SNP effects. An improved Lasso prior was assigned to SNP effects which were estimated by searching the maximum joint posterior mode. The residual polygenic effect was included in the model to absorb many tiny SNP effects, which is treated as missing data in our EM algorithm. A series of simulation experiments were conducted to validate the proposed method, and the results showed that compared with SMMA, the new method can dramatically decrease the false-positive rate. The new method was also applied to the 50k SNP-panel dataset for genome-wide association study of milk production traits in Chinese Holstein cattle. Totally, 39 significant SNPs and their nearby 25 genes were found. The number of significant SNPs is remarkably fewer than that by SMMA which found 105 significant SNPs. Among 39 significant SNPs, 8 were also found by SMMA and several well-known QTLs or genes were confirmed again; furthermore, we also got some positional candidate gene with potential function of effecting milk production traits. These novel findings in our research should be valuable for further investigation.     

Genome Wide Association Study Identifies 20 Novel Promising Genes Associated with Milk Fatty Acid Traits in Chinese Holstein

Detecting genes associated with milk fat composition could provide valuable insights into the complex genetic networks of genes underling variation in fatty acids synthesis and point towards opportunities for changing milk fat composition via selective breeding. In this study, we conducted a genome-wide association study (GWAS) for 22 milk fatty acids in 784 Chinese Holstein cows with the PLINK software. Genotypes were obtained with the Illumina BovineSNP50 Bead chip and a total of 40,604 informative, high-quality single nucleotide polymorphisms (SNPs) were used. Totally, 83 genome-wide significant SNPs and 314 suggestive significant SNPs associated with 18 milk fatty acid traits were detected. Chromosome regions that affect milk fatty acid traits were mainly observed on BTA1, 2, 5, 6, 7, 9, 13, 14, 18, 19, 20, 21, 23, 26 and 27. Of these, 146 SNPs were associated with more than one milk fatty acid trait; most of studied fatty acid traits were significant associated with multiple SNPs, especially C18:0 (105 SNPs), C18 index (93 SNPs), and C14 index (84 SNPs); Several SNPs are close to or within the DGAT1, SCD1 and FASN genes which are well-known to affect milk composition traits of dairy cattle. Combined with the previously reported QTL regions and the biological functions of the genes, 20 novel promising candidates for C10:0, C12:0, C14:0, C14:1, C14 index, C18:0, C18:1n9c, C18 index, SFA, UFA and SFA/UFA were found, which composed of HTR1B, CPM, PRKG1, MINPP1, LIPJ, LIPK, EHHADH, MOGAT1, ECHS1, STAT1, SORBS1, NFKB2, AGPAT3, CHUK, OSBPL8, PRLR, IGF1R, ACSL3, GHR and OXCT1. Our findings provide a groundwork for unraveling the key genes and causal mutations affecting milk fatty acid traits in dairy cattle.