Mapping main, epistatic and sex-specific QTL for body composition in a chicken population divergently selected for low or high growth rate

New insights into genetics underlying of plumage color

Plumage color can be considered as a social signal in chickens and a breeding identification tool among breeders. The relationship between plumage color and trait groups of immunity, growth and fertility is still a controversial issue. This research aimed to determine the genome-wide additive and epistatic variants affecting plumage color variation in chickens using the chicken Illumina 60k high-density SNP array. Two scenarios of genome-wide additive association studies using all SNPs and independent SNPs were carried out. To perform epistatic association analysis, the LD pruning approach was used to reduce the complexity of the analysis. We detected seven novel significant loci using all of the SNPs in the model and 14 SNPs using the LD pruning approach associated with plumage color. Moreover, 89 significantly associated SNP-SNP interactions (P-value <10^-6 ) distributed in 25 chromosomes were identified, indicating that all of the signals together putatively influence the quantitative variation of plumage color. By annotating genes relevant to top SNPs, we have distinguished 18 potential candidate genes comprising HNF4beta, CKMT1B, TBC1D22A, RPL8, CACNA2D1, FZD4, SGMS1, IRF8, OPTN, LOC420362, TRABD, OvoDA1, DAD1, USP6, RBM12B, MIR1772, MIR1709 and MIR6696 and also 89 putative gene-gene combinations responsible for plumage color variation in chickens. Furthermore, several KEGG pathways including metabolic pathway, cytokine-cytokine receptor interaction, focal adhesion, melanogenesis, glycosaminoglycan biosynthesis-keratan sulfate and sphingolipid metabolism were enriched in the gene-set analysis. The results indicated that plumage color is a highly polygenic trait which, in turn, can be affected by multiple coding genes, regulatory genes and gene-gene epistasis interactions. In addition to genes with additive effects, epistatic genes with tiny individual effect sizes but significant effects in a pair have the potential to control plumage coloration in chickens.

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.

Genome-wide association study reveals novel loci associated with fat-deposition and meat-quality traits in Pekin ducks

Meat-quality traits play an essential role in meat poultry production. To determine the genetic mechanisms of meat quality in Pekin ducks, we performed a large-scale GWAS to identify quantitative trait loci affecting meat quality in Pekin ducks. We measured 10 traits in 542 Pekin ducks and genotyped each duck using genotyping-by-sequencing. The genetic parameters (genomic heritability, genetic correlation) for 10 meat-quality related traits were evaluated. Based on the large genotype-phenotype dataset, we performed GWASs for all of these traits. A total of 33 significant QTL (P < 3.03 × 10^-5 ) across 13 chromosomes were identified by loci-based analysis. Some newly identified candidate genes were discovered for fat-deposition and meat-quality traits, including PAG1 for body weight and eviscerated weight, INTU and NUP35 for abdominal fat weight and ratio, NUP3 and ARHGDIB for skin fat weight and ratio, GOLGA5 for breast muscle toughness and breast tenderness, and CTDSPL and PKP1 for breast muscle thickness. The current study is the first systematic report regarding duck meat quality.

Genetic Architecture of Innate Fear Behavior in Chickens

The genetic architecture of innate fear behavior in chickens is poorly understood. Here, we performed quantitative trait loci (QTL) analysis of innate responses to tonic immobility (TI) and open field (OF) fears in 242 newly hatched chicks of an F₂ population between the native Japanese Nagoya breed and the White Leghorn breed using 881 single nucleotide polymorphism markers obtained by restriction site-associated DNA sequencing. At genome-wide 5% significance levels, four QTL for TI traits were revealed on chromosomes 1-3 and 24. Two of these loci had sex-specific effects on the traits. For OF traits, three QTL were revealed on chromosomes 2, 4 and 7. The TI and OF QTL identified showed no overlaps in genomic regions and different modes of inheritance. The three TI QTL and one OF QTL exerted antagonistic effects on the traits. The results demonstrated that context-dependent QTL underlie the variations in innate TI and OF behaviors.

Association of the thyroid hormone responsive spot 14 alpha gene with growth-related traits in Korean native chicken

Objective

Thyroid hormone responsive spot 14 alpha (THRSP) has been used to investigate the regulation of de novo lipogenesis because the variation of THRSP mRNA content in the tissue affects directly the ability of that tissue to synthetize lipids. Also, this gene responds to thyroid hormone stimulation and high level of carbohydrate feeding or insulin-injection. This study was carried out to investigate variations within THRSP and their effects on body and carcass weights in Korean native chicken (KNC).

Methods

A total of 585 chickens which represent the five lines of KNC (Black, Gray-Brown, Red-Brown, White, and Yellow-Brown) were reared and body weight data were recorded every two weeks from hatch until 20 weeks of age. Polymerase chain reaction-restriction fragment length polymorphism, DNA chips for Agilent 2100 Bioanalyzer, and Fluidigm Genotyping Technology, were applied to genotype selected markers. A linear mixed-effect model was used to access association between these single nucleotide polymorphism (SNP) markers and growth-related traits.

Results

A total of 30 polymorphisms were investigated in THRSP. Of these, nine SNPs for loci were selected to perform association analyses. Significant associations were detected between g.-49G>T SNP with body weight at 20 weeks of age (BW20), g.451T>C SNP with growth at 10 to 12 weeks of age (GR10-12), and g.1432A>C SNP with growth at 14 to 16 weeks trait (GR14-16) and body weight at 18 weeks of age (BW18). Moreover, diplotype of the THRSP gene significantly affected body weight at 12 weeks of age (BW12) and GR10-12 traits. Diplotype of ht1/ht2 was favorable for BW12 and GR10-12 traits.

Conclusion

These results suggest that THRSP can be regarded as a candidate gene for growth traits in KNC.

QTL Fine Mapping for Sex Determination Region in Bighead Carp (Hypophthalmichthys nobilis) and Comparison with Silver Carp (Hypophthalmichthys molitrix)

Bighead carp (Hypophthalmichthys nobilis) and silver carp (Hypophthalmichthys molitrix) are genetically close aquaculture fish in the Cyprinidae, which have been confirmed to hold XX/XY sex determination. However, genomic locations of potential sex-related loci in these two fishes are still unknown. In this study, a high-resolution genetic linkage map was constructed by using 2976 SNP and 924 microsatellite markers in a F1 full-sib family of bighead carp, the length of which spanned 2022.34 cM with an average inter-marker distance of 0.52 cM. Comparative genomics revealed a high level of genomic synteny between bighead carp and zebrafish as well as grass carp. QTL fine mapping for sex trait was performed based on this linkage map of bighead carp and an unpublished linkage map of silver carp. A map distance of 3.863 cM (69.787-73.650 cM) on LG19 of bighead carp and 4.705 cM (79.096-83.801 cM) on LG21 of silver carp was significantly associated with sex phenotypes, and these two LGs are homologous between two fish species. Fourteen markers harboring in these regions were in strong linkage disequilibrium with the sex phenotype variance explained (PVE) varying from 89 to 100%. Two common markers were mapped on the QTL regions of bighead carp and silver carp, suggesting that these two carp species may have similar genetic bases for sex determination. Eleven potentially sex-related genes were identified within or near the sex QTL markers in two species. This study provided insights into elucidating mechanisms and evolution of sex determination in cyprinid fishes.

Quantitative trait loci mapping for the shear force value in breast muscle of F2 chickens

The shear force value is one of the major traits that determine meat quality. In the present study, we performed QTL analysis for chicken breast muscle shear force value at 7 wk of age using 545 single nucleotide polymorphism (SNP) markers developed via restriction-site associated DNA sequencing (RAD-seq). An F2 resource family was generated by mating Oh-Shamo, a native Japanese chicken breed, and the White Plymouth Rock chicken breed. A total of 215 F2 birds were produced. Simple interval mapping revealed one significant main-effect QTL between 6.28 and 8.10 Mb SNPs on the chromosome Z with a logarithm of odds score of 5.53 at the genome-wide 5% level. At this QTL, the confidence interval, phenotypic variance explained, and additive effect were 26 cM, 12.24%, and -0.31 in males and -0.34 in females, respectively. No QTL with epistatic interaction effects were detected. To our knowledge, this is the first report on a QTL affecting the shear force value in the chicken breast muscle, using SNP markers derived from RAD-seq.

Quantitative trait loci for morphometric and mineral composition traits of the tibia bone in a broiler × layer cross

Many economic losses occur in the poultry industry due to leg fragility. Knowing the genomic regions that influence traits associated with the growth and composition of the leg's bone can help to improve the selection process leading to increased leg resistance to fracture. The present study aimed to map quantitative trait loci (QTL) for mineral composition and morphometric traits of the tibia in 478 animals from an F2 broiler × layer cross. The measurement of weight, length and width of Tibia was carried out at 42 days of age. Ash, dry matter, levels of calcium (Ca), phosphorus (P), magnesium (Mg), Zinc (Zn) and Calcium:Phosphorus (Ca:P) ratio were also recorded. The population was genotyped for 128 microsatellite markers and one single nucleotide polymorphism, covering 2630 cM of the chicken genome. A likelihood ratio test was performed to find QTLs. Additive and dominance effects of the QTLs were included in the model. In the chromosomes 2 (GGA2), 6 (GGA6), 8 (GGA8), 24 (GGA24) and 26 (GGA26) some suggestive QTLs (P&lt;0.00276) were mapped for tibia weight (GGA2 and GGA26), ash percentage (GGA2 and GGA6), dry matter percentage (GGA2), Ca (GGA8 and GGA24) and Ca:P ratio (GGA8), many of which are close to genes already identified as good candidates for those traits. The suggestive QTL on GGA2 has a pleiotropic effect on ash percentage, dry matter and bone weight, whereas in the GGA8 there seems to be two QTLs, one for Ca and another for Ca:P ratio. Thus, this study identified at least five genomic regions, in different chromosomes, that can be targeted for further research to identify potential mutations influencing the development and composition of leg bones in Gallus gallus.

Mapping of Quantitative Trait Loci for Growth and Carcass-Related Traits in Chickens Using a Restriction-Site Associated DNA Sequencing Method

In the present study, quantitative trait loci (QTLs) analysis was performed to identify the chromosomal positions of growth and carcass-related trait QTLs using 319 F₂ chickens obtained from intercrosses of an Oh-Shamo male and four White Plymouth Rock females. Body weight was measured weekly until the birds were 7 weeks old. Carcass-related traits were also measured at this timepoint. A genetic linkage map was constructed using 545 single nucleotide polymorphism (SNP) markers that were developed using a restriction-site associated DNA sequencing method. The linkage map included the 23 autosomes and the Z chromosome. Using simple interval QTL mapping, we were able to identify 10 significant and suggestive main-effect QTLs for growth and carcass-related traits present on chromosomes 1, 2, 3, 5, 8, 19, 24, and Z. These loci explained 5.60–16.52% of the phenotypic variances. The chromosomal positions of the 10 QTLs overlapped with those of previously reported QTLs, whereas the targeted traits varied. Our QTLs will aid future breeding programs in improving growth and meat yield of chickens (e.g., via marker-assisted selection), particularly in the Japanese brand chicken industry.

A genome-wide association study reveals novel genomic regions and positional candidate genes for fat deposition in broiler chickens

Background

Excess fat content in chickens has a negative impact on poultry production. The discovery of QTL associated with fat deposition in the carcass allows the identification of positional candidate genes (PCGs) that might regulate fat deposition and be useful for selection against excess fat content in chicken’s carcass. This study aimed to estimate genomic heritability coefficients and to identify QTLs and PCGs for abdominal fat (ABF) and skin (SKIN) traits in a broiler chicken population, originated from the White Plymouth Rock and White Cornish breeds.

Results

ABF and SKIN are moderately heritable traits in our broiler population with estimates ranging from 0.23 to 0.33. Using a high density SNP panel (355,027 informative SNPs), we detected nine unique QTLs that were associated with these fat traits. Among these, four QTL were novel, while five have been previously reported in the literature. Thirteen PCGs were identified that might regulate fat deposition in these QTL regions: JDP2, PLCG1, HNF4A, FITM2, ADIPOR1, PTPN11, MVK, APOA1, APOA4, APOA5, ENSGALG00000000477, ENSGALG00000000483, and ENSGALG00000005043. We used sequence information from founder animals to detect 4843 SNPs in the 13 PCGs. Among those, two were classified as potentially deleterious and two as high impact SNPs.

Conclusions

This study generated novel results that can contribute to a better understanding of fat deposition in chickens. The use of high density array of SNPs increases genome coverage and improves QTL resolution than would have been achieved with low density. The identified PCGs were involved in many biological processes that regulate lipid storage. The SNPs identified in the PCGs, especially those predicted as potentially deleterious and high impact, may affect fat deposition. Validation should be undertaken before using these SNPs for selection against carcass fat accumulation and to improve feed efficiency in broiler chicken production.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4779-6) contains supplementary material, which is available to authorized users.

Epistatic effects between pairs of the growth hormone secretagogue receptor 1a, growth hormone, growth hormone receptor, non-SMC condensin I complex, subunit G and stearoyl-CoA desaturase genes on carcass, price-related and fatty acid composition traits in Japanese Black cattle

Growth hormone secretagogue receptor 1a (GHSR1a), growth hormone (GH), growth hormone receptor (GHR), non-SMC condensin I complex, subunit G (NCAPG) and stearoyl-CoA desaturase (SCD), are known to play important roles in growth and lipid metabolisms. Single and epistatic effects of the five genes on carcass, price-related and fatty acid (FA) composition traits were analyzed in a commercial Japanese Black cattle population of Ibaraki Prefecture. A total of 650 steers and 116 heifers for carcass and price-related traits, and 158 steers for FA composition traits were used in this study. Epistatic effects between pairs of the five genes were found in several traits. Alleles showing strain-specific differences in the five genes had significant single and epistatic effects in some traits. The data suggest that a TG-repeat polymorphism of the GHSR1a.5'UTR-(TG)_n locus plays a central role in gene-gene epistatic interaction of FA composition traits in the adipose tissue of Japanese Black cattle.

Transcriptional analysis of abdominal fat in chickens divergently selected on bodyweight at two ages reveals novel mechanisms controlling adiposity: validating visceral adipose tissue as a dynamic endocrine and metabolic organ

Background

Decades of intensive genetic selection in the domestic chicken (Gallus gallus domesticus) have enabled the remarkable rapid growth of today’s broiler (meat-type) chickens. However, this enhanced growth rate was accompanied by several unfavorable traits (i.e., increased visceral fatness, leg weakness, and disorders of metabolism and reproduction). The present descriptive analysis of the abdominal fat transcriptome aimed to identify functional genes and biological pathways that likely contribute to an extreme difference in visceral fatness of divergently selected broiler chickens.

Methods

We used the Del-Mar 14 K Chicken Integrated Systems microarray to take time-course snapshots of global gene transcription in abdominal fat of juvenile [1-11 weeks of age (wk)] chickens divergently selected on bodyweight at two ages (8 and 36 wk). Further, a RNA sequencing analysis was completed on the same abdominal fat samples taken from high-growth (HG) and low-growth (LG) cockerels at 7 wk, the age with the greatest divergence in body weight (3.2-fold) and visceral fatness (19.6-fold).

Results

Time-course microarray analysis revealed 312 differentially expressed genes (FDR ≤ 0.05) as the main effect of genotype (HG versus LG), 718 genes in the interaction of age and genotype, and 2918 genes as the main effect of age. The RNA sequencing analysis identified 2410 differentially expressed genes in abdominal fat of HG versus LG chickens at 7 wk. The HG chickens are fatter and over-express numerous genes that support higher rates of visceral adipogenesis and lipogenesis. In abdominal fat of LG chickens, we found higher expression of many genes involved in hemostasis, energy catabolism and endocrine signaling, which likely contribute to their leaner phenotype and slower growth. Many transcription factors and their direct target genes identified in HG and LG chickens could be involved in their divergence in adiposity and growth rate.

Conclusions

The present analyses of the visceral fat transcriptome in chickens divergently selected for a large difference in growth rate and abdominal fatness clearly demonstrate that abdominal fat is a very dynamic metabolic and endocrine organ in the chicken. The HG chickens overexpress many transcription factors and their direct target genes, which should enhance in situ lipogenesis and ultimately adiposity. Our observation of enhanced expression of hemostasis and endocrine-signaling genes in diminished abdominal fat of LG cockerels provides insight into genetic mechanisms involved in divergence of abdominal fatness and somatic growth in avian and perhaps mammalian species, including humans.

Electronic supplementary material

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

Genome-Wide Association Study for Muscle Fat Content and Abdominal Fat Traits in Common Carp (Cyprinus carpio)

Muscle fat content is an important phenotypic trait in fish, as it affects the nutritional, technical and sensory qualities of flesh. To identify loci and candidate genes associated with muscle fat content and abdominal fat traits, we performed a genome-wide association study (GWAS) using the common carp 250 K SNP assay in a common carp F₂ resource population. A total of 18 loci surpassing the genome-wide suggestive significance level were detected for 4 traits: fat content in dorsal muscle (MFdo), fat content in abdominal muscle (MFab), abdominal fat weight (AbFW), and AbFW as a percentage of eviscerated weight (AbFP). Among them, one SNP (carp089419) affecting both AbFW and AbFP reached the genome-wide significance level. Ten of those loci were harbored in or near known genes. Furthermore, relative expressions of 5 genes related to MFdo were compared using dorsal muscle samples with high and low phenotypic values. The results showed that 4 genes were differentially expressed between the high and low phenotypic groups. These genes are, therefore, prospective candidate genes for muscle fat content: ankyrin repeat domain 10a (ankrd10a), tetratricopeptide repeat, ankyrin repeat and coiled-coil containing 2 (tanc2), and four jointed box 1 (fjx1) and choline kinase alpha (chka). These results offer valuable insights into the complex genetic basis of fat metabolism and deposition.

Genome-wide association study of 8 carcass traits in Jinghai Yellow chickens using specific-locus amplified fragment sequencing technology

Carcass traits are important to the commercial chicken industry, and understanding the genetics of these traits will be useful in the development of commercially viable varieties of chickens. We conducted a genome-wide association study based on 8 carcass trait phenotypes in a population of 400 43-week-old Jinghai Yellow chickens. Specific-locus amplified fragment sequencing technology was used to identify 90,961 single nucleotide polymorphisms (SNP) distributed among 29 chromosomes and the mitochondrial genome. SNP that were significantly associated with phenotypic traits were identified by a simple general linear model. Fifteen SNP attained genome-wide significance (P < 1.87E−6) and were associated with 5 of the 8 carcass traits; only one SNP was significantly associated with 2 traits (foot weight and wing weight). Twelve genes were associated with these 15 SNP. A region of chromosome 4 between 75.5 and 76.1 Mb was associated with carcass weight, foot weight, and wing weight. An 84-kb region on chromosome 3 (51.2 Mb) was associated with eviscerated weight and semi-eviscerated weight.

Genome-Wide Linkage Analysis Identifies Loci for Physical Appearance Traits in Chickens

Physical appearance traits, such as feather-crested head, comb size and type, beard, wattles size, and feathered feet, are used to distinguish between breeds of chicken and also may be associated with economic traits. In this study, a genome-wide linkage analysis was used to identify candidate regions and genes for physical appearance traits and to potentially provide further knowledge of the molecular mechanisms that underlie these traits. The linkage analysis was conducted with an F2 population derived from Beijing-You chickens and a commercial broiler line. Single-nucleotide polymorphisms were analyzed using the Illumina 60K Chicken SNP Beadchip. The data were used to map quantitative trait loci and genes for six physical appearance traits. A 10-cM/0.51-Mb region (0.0−10.0 cM/0.00−0.51 Mb) with 1% genome-wide significant level on LGE22C19W28_E50C23 linkage group (LGE22) for crest trait was identified, which is likely very closely linked to the HOXC8. A QTL with 5% chromosome-wide significant level for comb weight, which partly overlaps with a region identified in a previous study, was identified at 74 cM/25.55 Mb on chicken (Gallus gallus; GG) chromosome 3 (i.e., GGA3). For beard and wattles traits, an identical region 11 cM/2.23 Mb (0.0−11.0 cM/0.00−2.23 Mb) including WNT3 and GH genes on GGA27 was identified. Two QTL with 1% genome-wide significant level for feathered feet trait, one 9-cM/2.80-Mb (48.0-57.0/13.40-16.20 Mb) region on GGA13, and another 12-cM/1.45-Mb (41.0−53.0 cM/11.37−12.82 Mb) region on GGA15 were identified. These candidate regions and genes provide important genetic information for the physical appearance traits in chicken.

Detecting signatures of selection in nine distinct lines of broiler chickens

Modern commercial chickens have been bred for one of two specific purposes: meat production (broilers) or egg production (layers). This has led to large phenotypic changes, so that the genomic signatures of selection may be detectable using statistical techniques. Genetic differentiation between nine distinct broiler lines was calculated using Weir and Cockerham's pairwise FST estimator for 11 003 genome-wide markers to identify regions showing evidence of differential selection across lines. Differentiation measures were averaged into overlapping sliding windows for each line, and a permutation approach was used to determine the significance of each window. A total of 51 regions were found to show significant differentiation between the lines. Several lines were consistently found to share significant regions, suggesting that the pattern of line divergence is related to selection for broiler traits. The majority of the 51 regions contain QTL relating to broiler traits, but only five of them were found to be significantly enriched for broiler QTL, including a region on chromosome 27 containing 39 broiler QTL and 114 genes. Additionally, a number of these regions have been identified by other selection mapping studies. This study has identified a large number of potential selection signatures, and further tests with higher-density marker data may narrow these regions down to individual genes.

A genome scan revealed significant associations of growth traits with a major QTL and GHR2 in tilapia

Growth is an important trait in animal breeding. However, the genetic effects underpinning fish growth variability are still poorly understood. QTL mapping and analysis of candidate genes are effective methods to address this issue. We conducted a genome-wide QTL analysis for growth in tilapia. A total of 10, 7 and 8 significant QTLs were identified for body weight, total length and standard length at 140 dph, respectively. The majority of these QTLs were sex-specific. One major QTL for growth traits was identified in the sex-determining locus in LG1, explaining 71.7%, 67.2% and 64.9% of the phenotypic variation (PV) of body weight, total length and standard length, respectively. In addition, a candidate gene GHR2 in a QTL was significantly associated with body weight, explaining 13.1% of PV. Real-time qPCR revealed that different genotypes at the GHR2 locus influenced the IGF-1 expression level. The markers located in the major QTL for growth traits could be used in marker-assisted selection of tilapia. The associations between GHR2 variants and growth traits suggest that the GHR2 gene should be an important gene that explains the difference in growth among tilapia species.

Re-sequencing data for refining candidate genes and polymorphisms in QTL regions affecting adiposity in chicken

In this study, we propose an approach aiming at fine-mapping adiposity QTL in chicken, integrating whole genome re-sequencing data. First, two QTL regions for adiposity were identified by performing a classical linkage analysis on 1362 offspring in 11 sire families obtained by crossing two meat-type chicken lines divergently selected for abdominal fat weight. Those regions, located on chromosome 7 and 19, contained a total of 77 and 84 genes, respectively. Then, SNPs and indels in these regions were identified by re-sequencing sires. Considering issues related to polymorphism annotations for regulatory regions, we focused on the 120 and 104 polymorphisms having an impact on protein sequence, and located in coding regions of 35 and 42 genes situated in the two QTL regions. Subsequently, a filter was applied on SNPs considering their potential impact on the protein function based on conservation criteria. For the two regions, we identified 42 and 34 functional polymorphisms carried by 18 and 24 genes, and likely to deeply impact protein, including 3 coding indels and 4 nonsense SNPs. Finally, using gene functional annotation, a short list of 17 and 4 polymorphisms in 6 and 4 functional genes has been defined. Even if we cannot exclude that the causal polymorphisms may be located in regulatory regions, this strategy gives a complete overview of the candidate polymorphisms in coding regions and prioritize them on conservation- and functional-based arguments.

Changes in variance explained by top SNP windows over generations for three traits in broiler chicken

The purpose of this study was to determine if the set of genomic regions inferred as accounting for the majority of genetic variation in quantitative traits remain stable over multiple generations of selection. The data set contained phenotypes for five generations of broiler chicken for body weight, breast meat, and leg score. The population consisted of 294,632 animals over five generations and also included genotypes of 41,036 single nucleotide polymorphism (SNP) for 4,866 animals, after quality control. The SNP effects were calculated by a GWAS type analysis using single step genomic BLUP approach for generations 1-3, 2-4, 3-5, and 1-5. Variances were calculated for windows of 20 SNP. The top ten windows for each trait that explained the largest fraction of the genetic variance across generations were examined. Across generations, the top 10 windows explained more than 0.5% but less than 1% of the total variance. Also, the pattern of the windows was not consistent across generations. The windows that explained the greatest variance changed greatly among the combinations of generations, with a few exceptions. In many cases, a window identified as top for one combination, explained less than 0.1% for the other combinations. We conclude that identification of top SNP windows for a population may have little predictive power for genetic selection in the following generations for the traits here evaluated.

Overview of Genomic Insights into Chicken Growth Traits Based on Genome-Wide Association Study and microRNA Regulation

Over the two past decades, a significant number of studies have observed animal growth traits to examine animal genetic mechanisms due to their ease of measurement and high heritability. Chicken which has a significant impact on fundamental biology is a major source of protein worldwide, making it an ideal model for examining animal growth trait development. The genetic mechanisms of chicken growth traits have been studied using quantitative trait loci mapping through genome-scan and candidate gene approaches, genome-wide association studies (GWAS), comparative genomic strategies, microRNA (miRNA) regulation of growth development analysis, and epigenomic analysis. This review focuses on chicken GWAS and miRNA regulation of growth traits. Several recently published GWAS reports showed that most genome-wide significant single nucleotide polymorphisms are located on chromosomes 1 and 4 in chickens. Chicken growth, particularly skeletal muscle growth and development, is greatly regulated by miRNA. Using dwarf and normal chickens, let-7b was found to be involved in determining chicken dwarf phenotypes by regulating growth hormone receptor gene expression.

Genome-wide interval mapping using SNPs identifies new QTL for growth, body composition and several physiological variables in an F2 intercross between fat and lean chicken lines

BACKGROUND

For decades, genetic improvement based on measuring growth and body composition traits has been successfully applied in the production of meat-type chickens. However, this conventional approach is hindered by antagonistic genetic correlations between some traits and the high cost of measuring body composition traits. Marker-assisted selection should overcome these problems by selecting loci that have effects on either one trait only or on more than one trait but with a favorable genetic correlation. In the present study, identification of such loci was done by genotyping an F2 intercross between fat and lean lines divergently selected for abdominal fatness genotyped with a medium-density genetic map (120 microsatellites and 1302 single nucleotide polymorphisms). Genome scan linkage analyses were performed for growth (body weight at 1, 3, 5, and 7 weeks, and shank length and diameter at 9 weeks), body composition at 9 weeks (abdominal fat weight and percentage, breast muscle weight and percentage, and thigh weight and percentage), and for several physiological measurements at 7 weeks in the fasting state, i.e. body temperature and plasma levels of IGF-I, NEFA and glucose. Interval mapping analyses were performed with the QTLMap software, including single-trait analyses with single and multiple QTL on the same chromosome.

RESULTS

Sixty-seven QTL were detected, most of which had never been described before. Of these 67 QTL, 47 were detected by single-QTL analyses and 20 by multiple-QTL analyses, which underlines the importance of using different statistical models. Close analysis of the genes located in the defined intervals identified several relevant functional candidates, such as ACACA for abdominal fatness, GHSR and GAS1 for breast muscle weight, DCRX and ASPSCR1 for plasma glucose content, and ChEBP for shank diameter.

CONCLUSIONS

The medium-density genetic map enabled us to genotype new regions of the chicken genome (including micro-chromosomes) that influenced the traits investigated. With this marker density, confidence intervals were sufficiently small (14 cM on average) to search for candidate genes. Altogether, this new information provides a valuable starting point for the identification of causative genes responsible for important QTL controlling growth, body composition and metabolic traits in the broiler chicken.

Comparative analysis of quantitative trait loci for body weight, growth rate and growth curve parameters from 3 to 72 weeks of age in female chickens of a broiler-layer cross

Background

Comparisons of quantitative trait loci (QTL) for growth and parameters of growth curves assist in understanding the genetics and ultimately the physiology of growth. Records of body weight at 3, 6, 12, 24, 48 and 72 weeks of age and growth rate between successive age intervals of about 500 F₂ female chickens of the Roslin broiler-layer cross were available for analysis. These data were analysed to detect and compare QTL for body weight, growth rate and parameters of the Gompertz growth function.

Results

Over 50 QTL were identified for body weight at specific ages and most were also detected in the nearest preceding and/or subsequent growth stage. The sum of the significant and suggestive additive effects for bodyweight at specific ages accounted for 23-43% of the phenotypic variation. A single QTL for body weight on chromosome 4 at 48 weeks of age had the largest additive effect (550.4 ± 68.0 g, 11.5% of the phenotypic variation) and a QTL at a similar position accounted 14.5% of the phenotypic variation at 12 weeks of age. Age specific QTL for growth rate were detected suggesting that there are specific genes that affect developmental processes during the different stages of growth. Relatively few QTL influencing Gompertz growth curve parameters were detected and overlapped with loci affecting growth rate. Dominance effects were generally not significant but from 12 weeks of age they exceeded the additive effect in a few cases. No evidence for epistatic QTL pairs was found.

Conclusions

The results confirm the location for body weight and body weight gain during growth that were identified in previous studies and were consistent with QTL for the parameters of the Gompertz growth function. Chromosome 4 explained a relatively large proportion of the observed growth variation across the different ages, and also harboured most of the detected QTL for Gompertz parameters, confirming its importance in controlling growth. Very few QTL were detected for body weight or gain at 48 and 72 weeks of age, probably reflecting the effect of differences in reproduction and random environmental effects.

Genetic dissection of growth traits in a Chinese indigenous × commercial broiler chicken cross

Background

In China, consumers often prefer indigenous broiler chickens over commercial breeds, as they have characteristic meat qualities requested within traditional culinary customs. However, the growth-rate of these indigenous breeds is slower than that of the commercial broilers, which means they have not yet reached their full economic value. Therefore, combining the valuable meat quality of the native chickens with the efficiency of the commercial broilers is of interest. In this study, we generated an F₂ intercross between the slow growing native broiler breed, Huiyang Beard chicken, and the fast growing commercial broiler breed, High Quality chicken Line A, and used it to map loci explaining the difference in growth rate between these breeds.

Results

A genome scan to identify main-effect loci affecting 24 growth-related traits revealed nine distinct QTL on six chromosomes. Many QTL were pleiotropic and conformed to the correlation patterns observed between phenotypes. Most of the mapped QTL were found in locations where growth QTL have been reported in other populations, although the effects were greater in this population. A genome scan for pairs of interacting loci identified a number of additional QTL in 10 other genomic regions. The epistatic pairs explained 6–8% of the residual phenotypic variance. Seven of the 10 epistatic QTL mapped in regions containing candidate genes in the ubiquitin mediated proteolysis pathway, suggesting the importance of this pathway in the regulation of growth in this chicken population.

Conclusions

The main-effect QTL detected using a standard one-dimensional genome scan accounted for a significant fraction of the observed phenotypic variance in this population. Furthermore, genes in known pathways present interesting candidates for further exploration. This study has thus located several QTL regions as promising candidates for further study, which will increase our understanding of the genetic mechanisms underlying growth-related traits in chickens.

Polymorphisms in FGFBP1 and FGFBP2 genes associated with carcass and meat quality traits in chickens

In the past, the focus of broiler breeding programs on yield and carcass traits improvement led to problems related to meat quality. Awareness of public concern for quality resulted in inclusion of meat quality traits in the evaluation process. Nevertheless, few genes associated with meat quality attributes are known. Previous studies mapped quantitative trait loci for weight at 35 and 42 days in a region of GGA4 flanked by the microsatellite markers, MCW0240 and LEI0063. In this region, there are 2 fibroblast growth factor binding protein (FGFBP) genes that play an important role in embryogenesis, cellular differentiation, and proliferation in chickens. The objective of this study was to identify and associate single nucleotide polymorphisms (SNPs) in FGFBP1 and FGFBP2 with performance, carcass, and meat quality in experimental and commercial chicken populations. In the commercial population, SNP g.2014G>A in FGFBP1 was associated with decreased carcass weight (P < 0.05), and SNP g.651G>A in FGFBP2 was associated with thawing loss and meat redness content (P < 0.05). Four haplotypes were constructed based on 2 SNPs and were associated with breast weight, thawing loss, and meat redness content. The diplotypes were associated with thawing loss, lightness, and redness content. The SNPs evaluated in the present study may be used as markers in poultry breeding programs to aid in improving growth and meat quality traits.

A genome-wide scan of selective sweeps in two broiler chicken lines divergently selected for abdominal fat content

Background

Genomic regions controlling abdominal fatness (AF) were studied in the Northeast Agricultural University broiler line divergently selected for AF. In this study, the chicken 60KSNP chip and extended haplotype homozygosity (EHH) test were used to detect genome-wide signatures of AF.

Results

A total of 5357 and 5593 core regions were detected in the lean and fat lines, and 51 and 57 reached a significant level (P<0.01), respectively. A number of genes in the significant core regions, including RB1, BBS7, MAOA, MAOB, EHBP1, LRP2BP, LRP1B, MYO7A, MYO9A and PRPSAP1, were detected. These genes may be important for AF deposition in chickens.

Conclusions

We provide a genome-wide map of selection signatures in the chicken genome, and make a contribution to the better understanding the mechanisms of selection for AF content in chickens. The selection for low AF in commercial breeding using this information will accelerate the breeding progress.

Quantitative trait loci segregating in crosses between New Hampshire and White Leghorn chicken lines: III. Fat deposition and intramuscular fat content

In this study, a genome scan was performed to detect genomic loci that affect fat deposition in white adipose tissues and muscles in 278 F (2) males of reciprocal crosses between the genetically and phenotypically extreme inbred chicken lines New Hampshire (NHI) and White Leghorn (WL77). Genome-wide highly significant quantitative trait loci (QTL) influencing fat deposition in white adipose tissues were found on GGA2 and 4. The peak QTL positions for different visceral and subcutaneous white adipose tissues were located between 41.4 and 112.4 Mb on GGA2 and between 76.2 and 78.7 Mb on GGA4, which explained 4.2-10.4% and 4.3-11.6% respectively of the phenotypic F (2) variances. Contrary to our expectations, the QTL allele descending from the lean line WL77 on GGA4 led to increased fat deposition. We suggest a transgressive action of the obesity allele only if it is not in the genetic background of the line WL77. Additional highly significant loci for subcutaneous adipose tissue mass were identified on GGA12 and 15. For intramuscular fat content, a suggestive QTL was located on GGA14. The analysed crosses provide a valuable resource for further fine mapping of fatness genes and subsequent gene discovery.

Whole genome QTL mapping for growth, meat quality and breast meat yield traits in turkey

Background

The turkey (Meleagris gallopavo) is an important agricultural species and is the second largest contributor to the world's poultry meat production. Demand of turkey meat is increasing very rapidly. Genetic markers linked to genes affecting quantitative traits can increase the selection response of animal breeding programs. The use of these molecular markers for the identification of quantitative trait loci, and subsequently fine-mapping of quantitative trait loci regions, allows for pinpointing of genes that underlie such economically important traits.

Results

The quantitative trait loci analyses of the growth curve, body weight, breast yield and the meat quality traits showed putative quantitative trait loci on 21 of the 27 turkey chromosomes covered by the linkage map. Forty-five quantitative trait loci were detected across all traits and these were found in 29 different regions on 21 chromosomes. Out of the 45 quantitative trait loci, twelve showed significant (p < 0.01) evidence of linkage while the remaining 33 showed suggestive evidence (p < 0.05) of linkage with different growth, growth curve, meat quality and breast yield traits.

Conclusion

A large number of quantitative trait loci were detected across the turkey genome, which affected growth, breast yield and meat quality traits. Pleiotropic effects or close linkages between quantitative trait loci were suggested for several of the chromosomal regions. The comparative analysis regarding the location of quantitative trait loci on different turkey, and on the syntenic chicken chromosomes, along with their phenotypic associations, revealed signs of functional conservation between these species.

Mapping quantitative trait loci for egg production traits in an F2 intercross of Oh-Shamo and White Leghorn chickens

We performed quantitative trait locus (QTL) analyses for egg production traits, including age at first egg (AFE) and egg production rates (EPR) measured every 4 weeks from 22 to 62 weeks of hen age, in a population of 421 F(2) hens derived from an intercross between the Oh-Shamo (Japanese Large Game) and White Leghorn breeds of chickens. Simple interval mapping revealed a main-effect QTL for AFE on chromosome 1 and four main-effect QTL for EPR on chromosomes 1 and 11 (three on chromosome 1 and one on chromosome 11) at the genome-wide 5% levels. Among the three EPR QTL on chromosome 1, two were identified at the early stage of egg laying (26-34 weeks of hen age) and the remaining one was discovered at the late stage (54-58 weeks). The alleles at the two EPR QTL derived from the Oh-Shamo breed unexpectedly increased the trait values, irrespective of the Oh-Shamo being inferior to the White Leghorn in the trait. This suggests that the Oh-Shamo, one of the indigenous Japanese breeds, is an untapped resource that is important for further improvement of current elite commercial laying chickens. In addition, six epistatic QTL were identified on chromosomes 2, 4, 7, 8, 17 and 19, where none of the above main-effect QTL were located. This is the first example of detection of epistatic QTL affecting egg production traits. The main and epistatic QTL identified accounted for 4-8% of the phenotypic variance. The total contribution of all QTL detected for each trait to the phenotypic and genetic variances ranged from 4.1% to 16.9% and from 11.5% to 58.5%, respectively.