Genome-wide association study dissects genetic architecture underlying longitudinal egg weights in chickens

Genome-wide association study identified candidate genes for egg production traits in the Longyan Shan-ma duck

Egg production is an important economic trait in layer ducks and understanding the genetics basis is important for their breeding. In this study, a genome-wide association study (GWAS) for egg production traits in 303 female Longyan Shan-ma ducks was performed based on a genotyping-by-sequencing strategy. Sixty-two single nucleotide polymorphisms (SNPs) associated with egg weight traits were identified (P < 9.48 × 10-5), including 8 SNPs at 5% linkage disequilibrium (LD)-based Bonferroni-corrected genome-wide significance level (P < 4.74 × 10-6). One hundred and nineteen SNPs were associated with egg number traits (P < 9.48 × 10-5), including 13 SNPs with 5% LD-based Bonferroni-corrected genome-wide significance (P < 4.74 × 10-6). These SNPs annotated 146 target genes which contained known candidate genes for egg production traits, such as prolactin and prolactin releasing hormone receptor. This study identified that these associated genes were significantly enriched in egg production-related pathways (P < 0.05), such as the oxytocin signaling, MAPK signaling, and calcium signaling pathways. It was notable that 18 genes were differentially expressed in ovarian tissues between higher and lower egg production in Shan-ma ducks. The identified potential candidate genes and pathways provide insight into the genetic basis underlying the egg production trait of layer ducks.

Genome-wide association analysis of eggshell color of an F2 generation population reveals candidate genes in chickens

Eggshell color is an important visual characteristic that affects consumer preferences for eggs. Eggshell color, which has moderate to high heritability, can be effectively enhanced through molecular marker selection. Various studies have been conducted on eggshell color at specific time points. However, few longitudinal data are available on eggshell color. Therefore, the objective of this study was to investigate eggshell color using the Commission International de L'Eclairage L*a*b* system with multiple measurements at different ages (age at the first egg and at 32, 36, 40, 44, 48, 52, 56, 60, 66, and 72 weeks) within the same individuals from an F2 resource population produced by crossing White Leghorn and Dongxiang Blue chicken. Using an Affymetrix 600 single nucleotide polymorphism (SNP) array, we estimated the genetic parameters of the eggshell color trait, performed genome-wide association studies (GWASs), and screened for the potential candidate genes. The results showed that pink-shelled eggs displayed a significant negative correlation between L* values and both a* and b* values. Genetic heritability based on SNPs showed that the heritability of L*, a*, and b* values ranged from 0.32 to 0.82 for pink-shelled eggs, indicating a moderate to high level of genetic control. The genetic correlations at each time point were mostly above 0.5. The major-effect regions affecting the pink eggshell color were identified in the 10.3-13.0 Mb interval on Gallus gallus chromosome 20, and candidate genes were selected, including SLC35C2, PCIF1, and SLC12A5. Minor effect polygenic regions were identified on chromosomes 1, 6, 9, 12, and 15, revealing 11 candidate genes, including MTMR3 and SLC35E4. Members of the solute carrier family play an important role in influencing eggshell color. Overall, our findings provide valuable insights into the phenotypic and genetic aspects underlying the variation in eggshell color. Using GWAS analysis, we identified multiple quantitative trait loci (QTLs) for pink eggshell color, including a major QTL on chromosome 20. Genetic variants associated with eggshell color may be used in genomic breeding programs.

New insights into the genetic loci related to egg weight and age at first egg traits in broiler breeder

Egg weight (EW) and age at first egg (AFE) are economically important traits in breeder chicken production. The genetic basis of these traits, however, is far from understood, especially for broiler breeders. In this study, genetic parameter estimation, genome-wide association analysis, meta-analysis, and selective sweep analysis were carried out to identify genetic loci associated with EW and AFE in 6,842 broiler breeders. The study found that the heritability of EW ranged from 0.42 to 0.44, while the heritability of AFE was estimated at 0.33 in the maternal line. Meta-analysis and selective sweep analysis identified two colocalized regions on GGA4 that significantly influenced EW at 32 wk (EW32W) and at 43 wk (EW43W) with both paternal and maternal lines. The genes AR, YIPF6, and STARD8 were located within the significant region (GGA4: 366.86-575.50 kb), potentially affecting EW through the regulation of follicle development, cell proliferation, and lipid transfer etc. The promising genes LCORL and NCAPG were positioned within the significant region (GGA4:75.35-75.42 Mb), potentially influencing EW through pleiotropic effects on growth and development. Additionally, 3 significant regions were associated with AFE on chromosomes GGA7, GGA19, and GGA27. All of these factors affected the AFE by influencing ovarian development. In our study, the genomic information from both paternal and maternal lines was used to identify genetic regions associated with EW and AFE. Two genomic regions and eight genes were identified as the most likely candidates affecting EW and AFE. These findings contribute to a better understanding of the genetic basis of egg production traits in broiler breeders and provide new insights into future technology development.

Utilizing Two Populations Derived from Tropical Maize for Genome-Wide Association Analysis of Banded Leaf and Sheath Blight Resistance

Banded leaf and sheath blight (BLSB) in maize is a soil-borne fungal disease caused by Rhizoctonia solani Kühn, resulting in significant yield losses. Investigating the genes responsible for regulating resistance to BLSB is crucial for yield enhancement. In this study, a multiparent maize population was developed, comprising two recombinant inbred line (RIL) populations totaling 442 F8RILs. The populations were generated by crossing two tropical inbred lines, CML444 and NK40-1, known for their BLSB resistance, as female parents, with the high-yielding but BLSB-susceptible inbred line Ye107 serving as the common male parent. Subsequently, we utilized 562,212 high-quality single nucleotide polymorphisms (SNPs) generated through genotyping-by-sequencing (GBS) for a comprehensive genome-wide association study (GWAS) aimed at identifying genes responsible for BLSB resistance. The objectives of this study were to (1) identify SNPs associated with BLSB resistance through genome-wide association analyses, (2) explore candidate genes regulating BLSB resistance in maize, and (3) investigate pathways involved in BLSB resistance and discover key candidate genes through Gene Ontology (GO) analysis. The GWAS analysis revealed nineteen SNPs significantly associated with BLSB that were consistently identified across four environments in the GWAS, with phenotypic variation explained (PVE) ranging from 2.48% to 11.71%. Screening a 40 kb region upstream and downstream of the significant SNPs revealed several potential candidate genes. By integrating information from maize GDB and the NCBI, we identified five novel candidate genes, namely, Zm00001d009723, Zm00001d009975, Zm00001d009566, Zm00001d009567, located on chromosome 8, and Zm00001d026376, on chromosome 10, related to BLSB resistance. These candidate genes exhibit association with various aspects, including maize cell membrane proteins and cell immune proteins, as well as connections to cell metabolism, transport, transcriptional regulation, and structural proteins. These proteins and biochemical processes play crucial roles in maize defense against BLSB. When Rhizoctonia solani invades maize plants, it induces the expression of genes encoding specific proteins and regulates corresponding metabolic pathways to thwart the invasion of this fungus. The present study significantly contributes to our understanding of the genetic basis of BLSB resistance in maize, offering valuable insights into novel candidate genes that could be instrumental in future breeding efforts to develop maize varieties with enhanced BLSB resistance.

Genetic Architecture of Abdominal Fat Deposition Revealed by a Genome-Wide Association Study in the Laying Chicken

Fat has a high energy density, and excessive fatness has been recognized as a problem for egg production and the welfare of chickens. The identification of a genetic polymorphism controlling fat deposition would be helpful to select against excessive fatness in the laying hen. This study aimed to estimate genomic heritability and identify the genetic architecture of abdominal fat deposition in a population of chickens from a Dongxiang blue-shelled local breed crossbred with the White Leghorn. A genome-wide association study was conducted on abdominal fat percentage, egg production and body weights using a sample of 1534 hens genotyped with a 600 K Chicken Genotyping Array. The analysis yielded a heritability estimate of 0.19 ± 0.04 for abdominal fat percentage; 0.56 ± 0.04 for body weight at 72 weeks; 0.11 ± 0.03 for egg production; and 0.24 ± 0.04 for body weight gain. The genetic correlation of abdominal fat percentage with egg production between 60 and 72 weeks of age was -0.35 ± 0.18. This implies a potential trade-off between these two traits related to the allocation of resources. Strong positive genetic correlations were found between fat deposition and weight traits. A promising locus close to COL12A1 on chromosome 3, associated with abdominal fat percent, was found in the present study. Another region located around HTR2A on chromosome 1, where allele substitution was predicted to be associated with body weight gain, accounted for 2.9% of phenotypic variance. Another region located on chromosome 1, but close to SOX5, was associated with egg production. These results may be used to influence the balanced genetic selection for laying hens.

Disentangling clustering configuration intricacies for divergently selected chicken breeds

Divergently selected chicken breeds are of great interest not only from an economic point of view, but also in terms of sustaining diversity of the global poultry gene pool. In this regard, it is essential to evaluate the classification (clustering) of varied chicken breeds using methods and models based on phenotypic and genotypic breed differences. It is also important to implement new mathematical indicators and approaches. Accordingly, we set the objectives to test and improve clustering algorithms and models to discriminate between various chicken breeds. A representative portion of the global chicken gene pool including 39 different breeds was examined in terms of an integral performance index, i.e., specific egg mass yield relative to body weight of females. The generated dataset was evaluated within the traditional, phenotypic and genotypic classification/clustering models using the k-means method, inflection points clustering, and admixture analysis. The latter embraced SNP genotype datasets including a specific one focused on the performance-associated NCAPG-LCORL locus. The k-means and inflection points analyses showed certain discrepancies between the tested models/submodels and flaws in the produced cluster configurations. On the other hand, 11 core breeds were identified that were shared between the examined models and demonstrated more adequate clustering and admixture patterns. These findings will lay the foundation for future research to improve methods for clustering as well as genome- and phenome-wide association/mediation analyses.

Genome-wide association and selective sweep analyses reveal genetic loci for FCR of egg production traits in ducks

Background

As a major economic trait in poultry, egg production efficiency attracts widespread interest in breeding and production. However, limited information is available about the underlying genetic architecture of egg production traits in ducks. In this paper, we analyzed six egg production-related traits in 352 F₂ ducks derived from reciprocal crosses between mallard and Pekin ducks.

Results

Feed conversation ratio (FCR) was positively correlated with feed intake but negatively correlated with egg-related traits, including egg weight and egg production, both phenotypically and genetically. Estimates of pedigree-based heritability were higher than 0.2 for all traits investigated, except hip-width. Based on whole-genome sequencing data, we conducted genome-wide association studies to identify genomic regions associated with these traits. In total, 11 genomic regions were associated with FCR. No genomic regions were identified as significantly associated with hip-width, total feed intake, average daily feed intake, and total egg production. Analysis of selective sweeps between mallard and Pekin ducks confirmed three of these genomic regions on chromosomes 13, 3 and 6. Within these three regions, variants in candidate genes that were in linkage disequilibrium with the GWAS leader single nucleotide polymorphisms (SNPs) (Chr13:2,196,728, P = 7.05 × 10^–14; Chr3:76,991,524, P = 1.06 × 10^–12; Chr6:20,356,803, P = 1.14 × 10^–10) were detected. Thus, we identified 31 potential candidate genes associated with FCR, among which the strongest candidates are those that are highly expressed in tissues involved in reproduction and nervous system functions of ducks: CNTN4, CRBR, GPR63, KLHL32, FHL5, TRNT1, MANEA, NDUFAF4, and SCD.

Conclusions

For the first time, we report the identification of genomic regions that are associated with FCR in ducks and our results illustrate the genomic changes that occurred during their domestication and are involved in egg production efficiency.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00684-5.

Identification and Functional Annotation of Genes Related to Bone Stability in Laying Hens Using Random Forests

Skeletal disorders, including fractures and osteoporosis, in laying hens cause major welfare and economic problems. Although genetics have been shown to play a key role in bone integrity, little is yet known about the underlying genetic architecture of the traits. This study aimed to identify genes associated with bone breaking strength and bone mineral density of the tibiotarsus and the humerus in laying hens. Potentially informative single nucleotide polymorphisms (SNP) were identified using Random Forests classification. We then searched for genes known to be related to bone stability in close proximity to the SNPs and identified 16 potential candidates. Some of them had human orthologues. Based on our findings, we can support the assumption that multiple genes determine bone strength, with each of them having a rather small effect, as illustrated by our SNP effect estimates. Furthermore, the enrichment analysis showed that some of these candidates are involved in metabolic pathways critical for bone integrity. In conclusion, the identified candidates represent genes that may play a role in the bone integrity of chickens. Although further studies are needed to determine causality, the genes reported here are promising in terms of alleviating bone disorders in laying hens.

Hierarchical mixed-model expedites genome-wide longitudinal association analysis

A hierarchical random regression model (Hi-RRM) was extended into a genome-wide association analysis for longitudinal data, which significantly reduced the dimensionality of repeated measurements. The Hi-RRM first modeled the phenotypic trajectory of each individual using a RRM and then associated phenotypic regressions with genetic markers using a multivariate mixed model (mvLMM). By spectral decomposition of genomic relationship and regression covariance matrices, the mvLMM was transformed into a multiple linear regression, which improved computing efficiency while implementing mvLMM associations in efficient mixed-model association expedited (EMMAX). Compared with the existing RRM-based association analyses, the statistical utility of Hi-RRM was demonstrated by simulation experiments. The method proposed here was also applied to find the quantitative trait nucleotides controlling the growth pattern of egg weights in poultry data.

Detection of loci exhibiting pleiotropic effects on body weight and egg number in female broilers

The objective of the present study was to discover the genetic variants, functional candidate genes, biological processes and molecular functions underlying the negative genetic correlation observed between body weight (BW) and egg number (EN) traits in female broilers. To this end, first a bivariate genome-wide association and second stepwise conditional-joint analyses were performed using 2586 female broilers and 240 k autosomal SNPs. The aforementioned analyses resulted in a total number of 49 independent cross-phenotype (CP) significant SNPs with 35 independent markers showing antagonistic action i.e., positive effects on one trait and negative effects on the other trait. A number of 33 independent CP SNPs were located within 26 and 14 protein coding and long non-coding RNA genes, respectively. Furthermore, 26 independent markers were situated within 44 reported QTLs, most of them related to growth traits. Investigation of the functional role of protein coding genes via pathway and gene ontology analyses highlighted four candidates (CPEB3, ACVR1, MAST2 and CACNA1H) as most plausible pleiotropic genes for the traits under study. Three candidates (CPEB3, MAST2 and CACNA1H) were associated with antagonistic pleiotropy, while ACVR1 with synergistic pleiotropic action. Current results provide a novel insight into the biological mechanism of the genetic trade-off between growth and reproduction, in broilers.

Design and evaluation of a custom 50K Infinium SNP array for egg-type chickens

With the development of molecular genetics and high-throughput sequencing technology, genotyping arrays consisting of large numbers of SNP have raised great interest in animal and plant research. However, the application of commercial chicken 600K SNP arrays has varied in different populations of egg-type chickens. Moreover, their genotyping cost is too high for large-scale population applications. Herein, we independently developed a custom Illumina 50K BeadChip, named PhenoixChip-I, for egg-type chickens based on SNP from 479 sequenced individuals in 7 lines. We filtered and selected SNP with stringent criteria, such as high polymorphism, genome coverage, design score, and priorities. Finally, a total of 43,681 effective SNP successfully genotyped were included on our custom array. Approximately 14K SNP were previously reported to be associated with important economic traits in egg-type chickens. Subsequently, we verified the applicability and efficiency of the PhenoixChip-I SNP array from many aspects, including evaluating its use scientific research (population structure analysis and genome-wide association study) and the poultry breeding industry (genomic selection). The findings in our study will play a crucial role in accelerating the genetic improvement of egg-type chickens.

Deciphering the mode of action and position of genetic variants impacting on egg number in broiler breeders

Background

Aim of the present study was first to identify genetic variants associated with egg number (EN) in female broilers, second to describe the mode of their gene action (additive and/or dominant) and third to provide a list with implicated candidate genes for the trait. A number of 2586 female broilers genotyped with the high density (~ 600 k) SNP array and with records on EN (mean = 132.4 eggs, SD = 29.8 eggs) were used. Data were analyzed with application of additive and dominant multi-locus mixed models.

Results

A number of 7 additive, 4 dominant and 6 additive plus dominant marker-trait significant associations were detected. A total number of 57 positional candidate genes were detected within 50 kb downstream and upstream flanking regions of the 17 significant markers. Functional enrichment analysis pinpointed two genes (BHLHE40 and CRTC1) to be involved in the ‘entrainment of circadian clock by photoperiod’ biological process. Gene prioritization analysis of the positional candidate genes identified 10 top ranked genes (GDF15, BHLHE40, JUND, GDF3, COMP, ITPR1, ELF3, ELL, CRLF1 and IFI30). Seven prioritized genes (GDF15, BHLHE40, JUND, GDF3, COMP, ELF3, CRTC1) have documented functional relevance to reproduction, while two more prioritized genes (ITPR1 and ELL) are reported to be related to egg quality in chickens.

Conclusions

Present results have shown that detailed exploration of phenotype-marker associations can disclose the mode of action of genetic variants and help in identifying causative genes associated with reproductive traits in the species.

Endocrine and genetic factors affecting egg laying performance in chickens: a review

1. Egg-laying performance reflects the overall reproductive performance of breeding hens. The genetic traits for egg-laying performance have low or medium heritability, and, depending on the period involved, usually ranges from 0.16 to 0.64. Egg-laying in chickens is regulated by a combination of environmental, endocrine and genetic factors. 2. The main endocrine factors that regulate egg-laying are gonadotropin-releasing hormone (GnRH), prolactin (PRL), follicle-stimulating hormone (FSH) and luteinising hormone (LH). 3. In the last three decades, many studies have explored this aspect at a molecular genetic level. Recent studies identified 31 reproductive hormone-based candidate genes that were significantly associated with egg-laying performance. With the development of genome-sequencing technology, 64 new candidate genes and 108 single nucleotide polymorphisms (SNPs) related to egg-laying performance have been found using genome-wide association studies (GWAS), providing novel insights into the molecular genetic mechanisms governing egg production. At the same time, microRNAs that regulate genes responsible for egg-laying in chickens were reviewed. 4. Research on endocrinological and genetic factors affecting egg-laying performance will greatly improve the reproductive performance of chickens and promote the protection, development, and utilisation of poultry. This review summarises studies on the endocrine and genetic factors of egg-laying performance in chickens from 1972 to 2019.

Multiple association analysis of loci and candidate genes that regulate body size at three growth stages in Simmental beef cattle

Background

Body size traits as one of the main breeding selection criteria was widely used to monitor cattle growth and to evaluate the selection response. In this study, body size was defined as body height (BH), body length (BL), hip height (HH), heart size (HS), abdominal size (AS), and cannon bone size (CS). We performed genome-wide association studies (GWAS) of these traits over the course of three growth stages (6, 12 and 18 months after birth) using three statistical models, single-trait GWAS, multi-trait GWAS and LONG-GWAS. The Illumina Bovine HD 770 K BeadChip was used to identify genomic single nucleotide polymorphisms (SNPs) in 1217 individuals.

Results

In total, 19, 29, and 10 significant SNPs were identified by the three models, respectively. Among these, 21 genes were promising candidate genes, including SOX2, SNRPD1, RASGEF1B, EFNA5, PTBP1, SNX9, SV2C, PKDCC, SYNDIG1, AKR1E2, and PRIM2 identified by single-trait analysis; SLC37A1, LAP3, PCDH7, MANEA, and LHCGR identified by multi-trait analysis; and P2RY1, MPZL1, LINGO2, CMIP, and WSCD1 identified by LONG-GWAS.

Conclusions

Multiple association analysis was performed for six growth traits at each growth stage. These findings offer valuable insights for the further investigation of potential genetic mechanism of growth traits in Simmental beef cattle.

Genome-wide association study provides insights into the genetic architecture of bone size and mass in chickens

Bone size is an important trait for chickens because of its association with osteoporosis in layers and meat production in broilers. Here, we employed high density genotyping platforms to detect candidate genes for bone traits. Estimates of the narrow heritabilities ranged from 0.37 ± 0.04 for shank length to 0.59 ± 0.04 for tibia length. The dominance heritability was 0.12 ± 0.04 for shank length. Using a linear mixed model approach, we identified a promising locus within NCAPG on chromosome 4, which was associated with tibia length and mass, femur length and area, and shank length. In addition, three other loci were associated with bone size or mass at a Bonferroni-corrected genome-wide significance threshold of 1%. One region on chicken chromosome 1 between 168.38 and 171.82 Mb harbored HTR2A, LPAR6, CAB39L, and TRPC4. A second region that accounted for 2.2% of the phenotypic variance was located around WNT9A on chromosome 2, where allele substitution was predicted to be associated with tibia length. Four candidate genes identified on chromosome 27 comprising SPOP, NGFR, GIP, and HOXB3 were associated with tibia length and mass, femur length and area, and shank length. Genome partitioning analysis indicated that the variance explained by each chromosome was proportional to its length.

Genetic architecture related to contour feathers density in an F2 resource population via a genome-wide association study

The density of contour feathers is an important trait as it is closely related to heat dissipation in birds. Thus, identification of the major genes that control this trait will be useful to improve heat tolerance in chicken. So far, no GWAS study for the density of contour feathers in birds has been previously published; therefore, this study was aimed to identify genomic regions controlling the density of contour feathers. A total of 1252 hens were genotyped, using the 600 K Affymetrix Axiom Chicken Genotyping Array. The association analyses were performed using the GenABEL package in the R program. In brief, 146 significant SNP markers were mainly located on chromosome 1 and were identified to associate with the density of contour feathers in the current GWAS analysis. Moreover, we identified several within/nearby candidate genes (SUCLA2, DNAJC15, DHRS12, MLNR, and RB1) that are either directly or indirectly involved in the genetic control of the density of contour feathers in chicken. This study laid the foundation for studying the mechanism that underlies the density of chicken feathers. Furthermore, it is feasible to shear the back feathers of live chickens and measure the density of the feathers to improve heat tolerance in breeding practice.

Genome-wide association analysis of egg production performance in chickens across the whole laying period

BACKGROUND

Egg production is the most economically-important trait in layers as it directly influences benefits of the poultry industry. To better understand the genetic architecture of egg production, we measured traits including age at first egg (AFE), weekly egg number (EN) from onset of laying eggs to 80 weeks which was divided into five stage (EN1: from onset of laying eggs to 23 weeks, EN2: from 23 to 37 weeks, EN3: from 37 to 50 weeks, EN4: from 50 to 61 weeks, EN5: from 61 to 80 weeks) based on egg production curve and total egg number across the whole laying period (Total-EN). Then we performed genome-wide association studies (GWAS) in 1078 Rhode Island Red hens using a linear mixed model.

RESULTS

Estimates of pedigree and SNP-based genetic parameter showed that AFE and EN1 exhibited high heritability (0.51 ± 0.09, 0.53 ± 0.08), while the h2 for EN in other stages varied from low (0.07 ± 0.04) to moderate (0.24 ± 0.07) magnitude. Subsequently, seven univariate GWAS for AFE and ENs were carried out independently, from which a total of 161 candidate SNPs located on GGA1, GGA2, GGA5, GGA6, GGA9 and GGA24 were identified. Thirteen SNP located on GGA6 were associated with AFE and an interesting gene PRLHR that may affect AFE through regulating oxytocin secretion in chickens. Sixteen genome-wide significant SNPs associated with EN3 were in a strong linkage disequilibrium (LD) region spanning from 117.87 Mb to 118.36 Mb on GGA1 and the most significant SNP (rs315777735) accounted for 3.57% of phenotypic variance. Genes POLA1, PDK3, PRDX4 and APOO identified by annotating sixteen genome-wide significant SNPs can be considered as candidates associated with EN3. Unfortunately, our study did not find any candidate gene for the total egg number.

CONCLUSIONS

Findings in our study could provide promising genes and SNP markers to improve egg production performance based on marker-assisted breeding selection, while further functional validation is still needed in other populations.

Discovery and characterization of functional modules associated with body weight in broilers

Aim of the present study was to investigate whether body weight (BW) in broilers is associated with functional modular genes. To this end, first a GWAS for BW was conducted using 6,598 broilers and the high density SNP array. The next step was to search for positional candidate genes and QTLs within strong LD genomic regions around the significant SNPs. Using all positional candidate genes, a network was then constructed and community structure analysis was performed. Finally, functional enrichment analysis was applied to infer the functional relevance of modular genes. A total number of 645 positional candidate genes were identified in strong LD genomic regions around 11 genome-wide significant markers. 428 of the positional candidate genes were located within growth related QTLs. Community structure analysis detected 5 modules while functional enrichment analysis showed that 52 modular genes participated in developmental processes such as skeletal system development. An additional number of 14 modular genes (GABRG1, NGF, APOBEC2, STAT5B, STAT3, SMAD4, MED1, CACNB1, SLAIN2, LEMD2, ZC3H18, TMEM132D, FRYL and SGCB) were also identified as related to body weight. Taken together, current results suggested a total number of 66 genes as most plausible functional candidates for the trait examined.

Bidirectional Selection for Body Weight on Standing Genetic Variation in a Chicken Model

Experimental populations of model organisms provide valuable opportunities to unravel the genomic impact of selection in a controlled system. The Virginia body weight chicken lines represent a unique resource to investigate signatures of selection in a system where long-term, single-trait, bidirectional selection has been carried out for more than 60 generations. At 55 generations of divergent selection, earlier analyses of pooled genome resequencing data from these lines revealed that 14.2% of the genome showed extreme differentiation between the selected lines, contained within 395 genomic regions. Here, we report more detailed analyses of these data exploring the regions displaying within- and between-line genomic signatures of the bidirectional selection applied in these lines. Despite the strict selection regime for opposite extremes in body weight, this did not result in opposite genomic signatures between the lines. The lines often displayed a duality of the sweep signatures, where an extended region of homozygosity in one line, in contrast to mosaic pattern of heterozygosity in the other line. These haplotype mosaics consisted of short, distinct haploblocks of variable between-line divergence, likely the results of a complex demographic history involving bottlenecks, introgressions and moderate inbreeding. We demonstrate this using the example of complex haplotype mosaicism in the growth1 QTL. These mosaics represent the standing genetic variation available at the onset of selection in the founder population. Selection on standing genetic variation can thus result in different signatures depending on the intensity and direction of selection.

Identification of potential genomic regions and candidate genes for egg albumen quality by a genome-wide association study

Albumen quality is a leading economic trait in the chicken industry. Major studies have paid attention to genetic architecture underlying albumen quality. However, the putative quantitative trait locus (QTL) for this trait is still unclear. In this genome-wide association study, we used an F2 resource population to study longitudinal albumen quality. Seven single-nucleotide polymorphism (SNP) loci were found to be significantly (p<8.43×10-7) related to albumen quality by univariate analysis, while 11 SNPs were significantly (p<8.43×10-7) associated with albumen quality by multivariate analysis. A QTL on GGA4 had a pervasive function on albumen quality, including a SNP at the missense of NCAPG, and a SNP at the intergenic region of FGFPB1. It was further found that the putative QTLs at GGA1, GGA2, and GGA7 had the strongest effects on albumen height (AH) at 32 weeks, Haugh units (HU) at 44 weeks, and AH at 55 weeks. Moreover, novel SNPs on GGA5 and GGA3 were associated with AH and HU at 32, 44, and 48 weeks of age. These results confirmed the regions for egg weight that were detected in a previous study and were similar with QTL for albumen quality. These results showed that GGA4 had the strongest effect on albumen quality. Only a few significant loci were detected for most characteristics probably reflecting the attributes of a pleiotropic gene and a minor-polygene in quantitative traits.

Genetic architecture and candidate genes detected for chicken internal organ weight with a 600 K single nucleotide polymorphism array

Objective

Internal organs indirectly affect economic performance and well-being of animals. Study of internal organs during later layer period will allow full utilization of layer hens. Hence, we conducted a genome-wide association study (GWAS) to identify potential quantitative trait loci or genes that potentially contribute to internal organ weight.

Methods

A total of 1,512 chickens originating from White Leghorn and Dongxiang Blue-Shelled chickens were genotyped using high-density Affymetrix 600 K single nucleotide polymorphism (SNP) array. We conducted a GWAS, linkage disequilibrium analysis, and heritability estimated based on SNP information by using GEMMA, Haploview and GCTA software.

Results

Our results displayed that internal organ weights show moderate to high (0.283 to 0.640) heritability. Variance partitioned across chromosomes and chromosome lengths had a linear relationship for liver weight and gizzard weight (R² = 0.493, 0.753). A total of 23 highly significant SNPs that associated with all internal organ weights were mainly located on Gallus gallus autosome (GGA) 1 and GGA4. Six SNPs on GGA2 affected heart weight. After the final analysis, five top SNPs were in or near genes 5-Hydroxytryptamine receptor 2A, general transcription factor IIF polypeptide 2, WD repeat and FYVE domain containing 2, non-SMC condensin I complex subunit G, and sonic hedgehog, which were considered as candidate genes having a pervasive role in internal organ weights.

Conclusion

Our findings provide an understanding of the underlying genetic architecture of internal organs and are beneficial in the selection of chickens.

Genome-wide association study on chicken carcass traits using sequence data imputed from SNP array

Chicken carcass traits are economically important for the chicken industry. Detecting which genes affect chicken carcass traits is of great benefit to the genetic improvement of this important agricultural species. To investigate the genetic mechanism of carcass traits in chickens, we carried out a genome-wide association study (GWAS). A total of 435 Chinese indigenous chickens were phenotyped for carcass weight (CW), eviscerated weight with giblets (EWG), and eviscerated weight (EW) after slaughter at 91 days and were genotyped using a 600-K single nucleotide polymorphism (SNP) genotyping array. Twenty-four birds were selected for sequencing, and the 600 K SNP panel data were imputed to sequence data with the 24 birds as the reference. Univariate GWASs were performed with GEMMA software using the whole genome sequence data imputed from SNP chip data. Finally, 3, 25, and 63 suggestively significant SNPs were identified to be associated with carcass weight (CW), eviscerated weight with giblets (EWG), and eviscerated weight (EW), respectively. Six candidate genes, RNF219, SCEL, MYCBP2, ETS1, APLP2, and PRDM10 were detected. SCEL and MYCBP2 were potentially associated with these three traits, RNF219 and APLP2 were potentially associated with EWG and EW, and ETS1 and PRDM10 were only potentially associated with EWG and EW, respectively. Compared with forefathers' research, 10 reported QTLs associated with CW were located within a 5-Mb distance near the SNPs with P value lower than 1×10^-5. This study enriched the knowledge of the genetic mechanisms of chicken carcass traits.

Genome-Wide Association Analysis of Age-Dependent Egg Weights in Chickens

Egg weight (EW) is an economically-important trait and displays a consecutive increase with the hen's age. Because extremely large eggs cause a range of problems in the poultry industry, we performed a genome-wide association study (GWAS) in order to identify genomic variations that are associated with EW. We utilized the Affymetrix 600 K high density SNP array in a population of 1,078 hens at seven time points from day at first egg to 80 weeks age (EW80). Results reveal that a 90 Kb genomic region (169.42 Mb ~ 169.51 Mb) in GGA1 is significantly associated with EW36 and is also potentially associated with egg weight at 28, 56, and 66 week of age. The leading SNP could account for 3.66% of the phenotypic variation, while two promising genes (DLEU7 and MIR15A) can be mapped to this narrow significant region and may affect EW in a pleiotropic manner. In addition, one gene (CECR2 on GGA1) and two genes (MEIS1 and SPRED2 on GGA3), which involved in the processes of embryogenesis and organogenesis, were also considered to be candidates related to first egg weight (FEW) and EW56, respectively. Findings in our study could provide worthy theoretical basis to generate eggs of ideal size based on marker assisted breeding selection.

A genome-wide study to identify genes responsible for oviduct development in chickens

Molecular genetic tools provide a method for improving the breeding selection of chickens (Gallus gallus). Although some studies have identified genes affecting egg quality, little is known about the genes responsible for oviduct development. To address this issue, here we used a genome-wide association (GWA) study to detect genes or genomic regions that are related to oviduct development in a chicken F₂ resource population by employing high-density 600 K single-nucleotide polymorphism (SNP) arrays. For oviduct length and weight, which exhibited moderate heritability estimates of 0.35 and 0.39, respectively, chromosome 1 (GGA1) explained 9.45% of the genetic variance, while GGA4 to GGA8 and GGA11 explained over 1% of the variance. Independent univariate genome-wide screens for oviduct length and weight detected 69 significant SNPs on GGA1 and 49 suggestive SNPs on GGA1, GGA4, and GGA8. One hundred and fourteen suggestive SNPs were associated with oviduct length, while 73 SNPs were associated with oviduct weight. The significant genomic regions affecting oviduct weight ranged from 167.79–174.29 Mb on GGA1, 73.16–75.70 Mb on GGA4, and 4.88–4.92 Mb on GGA8. The genes CKAP2, CCKAR, NCAPG, IGFBP3, and GORAB were shown to have potential roles in oviduct development. These genes are involved in cell survival, appetite, and growth control. Our results represent the first GWA analysis of genes controlling oviduct weight and length. The identification of genomic loci and potential candidate genes affecting oviduct development greatly increase our understanding of the genetic basis underlying oviduct development, which could have an impact on the selection of egg quality.

Genetic Architecture and Candidate Genes Identified for Follicle Number in Chicken

Follicular development has a major impact on reproductive performance. Most previous researchers focused on molecular mechanisms of follicular development. The genetic architecture underlying the number of follicle, however, has yet not to be thoroughly defined in chicken. Here we report a genome-wide association study for the genetic architecture determining the numbers of follicles in a large F₂ resource population. The results showed heritability were low to moderate (0.05-0.28) for number of pre-ovulatory follicles (POF), small yellow follicles (SYF) and atresia follicles (AF). The highly significant SNPs associated with SYF were mainly located on GGA17 and GGA28. Only four significant SNPs were identified for POF on GGA1. The variance partitioned across chromosomes and chromosome lengths had a linear relationship for SYF (R² = 0.58). The enriched genes created by the closest correspondent significant SNPs were found to be involved in biological pathways related to cell proliferation, cell cycle and cell survival. Two promising candidate genes, AMH and RGS3, were suggested to be prognostic biomarkers for SYF. In conclusion, this study offers the first evidence of genetic variance and positional candidate genes which influence the number of SYF in chicken. These identified informative SNPs may facilitate selection for an improved reproductive performance of laying hens.

Annotating long intergenic non-coding RNAs under artificial selection during chicken domestication

Background

Numerous biological functions of long intergenic non-coding RNAs (lincRNAs) have been identified. However, the contribution of lincRNAs to the domestication process has remained elusive. Following domestication from their wild ancestors, animals display substantial changes in many phenotypic traits. Therefore, it is possible that diverse molecular drivers play important roles in this process.

Results

We analyzed 821 transcriptomes in this study and annotated 4754 lincRNA genes in the chicken genome. Our population genomic analysis indicates that 419 lincRNAs potentially evolved during artificial selection related to the domestication of chicken, while a comparative transcriptomic analysis identified 68 lincRNAs that were differentially expressed under different conditions. We also found 47 lincRNAs linked to special phenotypes.

Conclusions

Our study provides a comprehensive view of the genome-wide landscape of lincRNAs in chicken. This will promote a better understanding of the roles of lincRNAs in domestication, and the genetic mechanisms associated with the artificial selection of domestic animals.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1036-6) contains supplementary material, which is available to authorized users.

Genetic architecture of bone quality variation in layer chickens revealed by a genome-wide association study

Skeletal problems in layer chickens are gaining attention due to animal welfare and economic losses in the egg industry. The genetic improvement of bone traits has been proposed as a potential solution to these issues; however, genetic architecture is not well understood. We conducted a genome-wide association study (GWAS) on bone quality using a sample of 1534 hens genotyped with a 600 K Chicken Genotyping Array. Using a linear mixed model approach, a novel locus close to GSG1L, associated with femur bone mineral density (BMD), was uncovered in this study. In addition, nine SNPs in genes were associated with bone quality. Three of these genes, RANKL, ADAMTS and SOST, were known to be associated with osteoporosis in humans, which makes them good candidate genes for osteoporosis in chickens. Genomic partitioning analysis supports the fact that common variants contribute to the variations of bone quality. We have identified several strong candidate genes and genomic regions associated with bone traits measured in end-of-lay cage layers, which accounted for 1.3–7.7% of the phenotypic variance. These SNPs could provide the relevant information to help elucidate which genes affect bone quality in chicken.

Genetic Mapping of Quantitative Trait Loci for Egg Production and Egg Quality Traits in Chickens: a Review

Chickens display a wide spectrum of phenotypic variations in quantitative traits such as egg-related traits. Quantitative trait locus (QTL) analysis is a statistical method used to understand the relationship between phenotypic (trait measurements) and genotypic data (molecular markers). We have performed QTL analyses for egg-related traits using an original resource population based on the Japanese Large Game (Oh-Shamo) and the White Leghorn breeds of chickens. In this article, we summarize the results of our extensive QTL analyses for 11 and 66 traits for egg production and egg quality, respectively. We reveal that at least 30 QTL regions on 17 different chromosomes affect phenotypic variation in egg-related traits. Each locus had an age-specific effect on traits, and a variety in effects was also apparent, such as additive, dominance, and epistatic-interaction effects. Although genome-wide association study (GWAS) is suitable for gene-level resolution mapping of GWAS loci with additive effects, QTL mapping studies enable us to comprehensively understand genetic control, such as chromosomal regions, genetic contribution to phenotypic variance, mode of inheritance, and age-specificity of both common and rare alleles. QTL analyses also describe the relationship between genotypes and phenotypes in experimental populations. Accumulation of QTL information, including GWAS loci, is also useful for studies of population genomics approached without phenotypic data in order to validate the identified genomic signatures of positive selection. The combination of QTL studies and next-generation sequencing techniques with uncharacterized genetic resources will enhance current understanding of the relationship between genotypes and phenotypes in livestock animals.

Genome-Wide Detection of Selective Signatures in Chicken through High Density SNPs

Chicken is recognized as an excellent model for studies of genetic mechanism of phenotypic and genomic evolution, with large effective population size and strong human-driven selection. In the present study, we performed Extended Haplotype Homozygosity (EHH) tests to identify significant core regions employing 600K SNP Chicken chip in an F2 population of 1,534 hens, which was derived from reciprocal crosses between White Leghorn and Dongxiang chicken. Results indicated that a total of 49,151 core regions with an average length of 9.79 Kb were identified, which occupied approximately 52.15% of genome across all autosomes, and 806 significant core regions attracted us mostly. Genes in candidate regions may experience positive selection and were considered to have possible influence on beneficial economic traits. A panel of genes including AASDHPPT, GDPD5, PAR3, SOX6, GPC1 and a signal pathway of AKT1 were detected with the most extreme P-values. Further enrichment analyses indicated that these genes were associated with immune function, sensory organ development and neurogenesis, and may have experienced positive selection in chicken. Moreover, some of core regions exactly overlapped with genes excavated in our previous GWAS, suggesting that these genes have undergone positive selection may affect egg production. Findings in our study could draw a comparatively integrate genome-wide map of selection signature in the chicken genome, and would be worthy for explicating the genetic mechanisms of phenotypic diversity in poultry breeding.

Genome-Wide Association Studies for Comb Traits in Chickens

The comb, as a secondary sexual character, is an important trait in chicken. Indicators of comb length (CL), comb height (CH), and comb weight (CW) are often selected in production. DNA-based marker-assisted selection could help chicken breeders to accelerate genetic improvement for comb or related economic characters by early selection. Although a number of quantitative trait loci (QTL) and candidate genes have been identified with advances in molecular genetics, candidate genes underlying comb traits are limited. The aim of the study was to use genome-wide association (GWA) studies by 600 K Affymetrix chicken SNP arrays to detect genes that are related to comb, using an F2 resource population. For all comb characters, comb exhibited high SNP-based heritability estimates (0.61-0.69). Chromosome 1 explained 20.80% genetic variance, while chromosome 4 explained 6.89%. Independent univariate genome-wide screens for each character identified 127, 197, and 268 novel significant SNPs with CL, CH, and CW, respectively. Three candidate genes, VPS36, AR, and WNT11B, were determined to have a plausible function in all comb characters. These genes are important to the initiation of follicle development, gonadal growth, and dermal development, respectively. The current study provides the first GWA analysis for comb traits. Identification of the genetic basis as well as promising candidate genes will help us understand the underlying genetic architecture of comb development and has practical significance in breeding programs for the selection of comb as an index for sexual maturity or reproduction.

The association of very low-density lipoprotein receptor (VLDLR) haplotypes with egg production indicates VLDLR is a candidate gene for modulating egg production

The very low-density lipoprotein receptor (VLDLR) transports egg yolk precursors into oocytes. However, our knowledge of the distribution patterns of VLDLR variants among breeds and their relationship to egg production is still incomplete. In this study, eight single nucleotide polymorphisms (SNPs) that account for 87% of all VLDLR variants were genotyped in Nick Chick (NC, n=91), Lohmann Brown (LohB, n=50) and Lueyang (LY, n=381) chickens, the latter being an Chinese indigenous breed. Egg production by NC and LY chickens was recorded from 17 to 50 weeks. Only four similar haplotypes were found in NC and LohB, of which two accounted for 100% of all NC haplotypes and 92.5% of LohB haplotypes. In contrast, there was considerable haplotypic diversity in LY. Comparison of egg production in LY showed that hens with NC-like haplotypes had a significantly higher production (p < 0.05) than those without the haplotypes. However, VLDLR expression was not significantly different between the haplotypes. These findings indicate a divergence in the distribution of VLDLR haplotypes between selected and non-selected breeds and suggest that the near fixation of VLDLR variants in NC and LohB is compatible with signature of selection. These data also support VLDLR as a candidate gene for modulating egg production.