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

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

Genetic diversity of bitter taste receptor gene family in Sichuan domestic and Tibetan chicken populations

The sense of bitter taste plays a critical role in animals as it can help them to avoid intake of toxic and harmful substances. Previous research had revealed that chicken has only three bitter taste receptor genes (Tas2r1, Tas2r2 and Tas2r7). To better understand the genetic polymorphisms and importance of bitter taste receptor genes (Tas2rs) in chicken, here, we sequenced Tas2rs of 30 Sichuan domestic chickens and 30 Tibetan chickens. Thirteen single-nucleotide polymorphisms (SNPs) including three nonsynonymous mutations (m.359G >C, m.503C >A and m.583A >G) were detected in Tas2r1 (m. is the abbreviation for mutation); three SNPs were detected in Tas2r2, but none of them were missense mutation; eight SNPs were detected in Tas2r7 including six nonsynonymous substitutions (m.178G >A, m.421A >C, m.787C >T, m.832G >T, m.907A >T and m.943G >A). Tajima’s D neutral test indicates that there is no population expansion in both populations, and the size of the population is relatively stable. All the three networks indicate that red jungle fowls share haplotypes with domestic chickens. In addition, we found that haplotypes H1 and HE1 were positively associated with high-altitude adaptation, whereas haplotypes H4 and HE4 showed a negative correlation with high-altitude adaptation in Tas2rs. Although, chicken has only three Tas2rs, our results showed that both Sichuan domestic chickens and Tibetan chickens have abundant haplotypes in Tas2rs, especially in Tas2r7, which might help chickens to recognize a wide variety of bitter-tasting compounds.     

Cold adaptation in the Asian tiger mosquito's native range precedes its invasion success in temperate regions

Adaptation to environmental conditions within the native range of exotic species can condition the invasion success of these species outside their range. The striking success of the Asian tiger mosquito, Aedes albopictus, to invade temperate regions has been attributed to the winter survival of diapause eggs in cold environments. In this study, we evaluate genetic polymorphisms (SNPs) and wing morphometric variation among three biogeographical regions of the native range of A. albopictus. Reconstructed demographic histories of populations show an initial expansion in Southeast Asia and suggest that marine regression during late Pleistocene and climate warming after the last glacial period favored expansion of populations in southern and northern regions, respectively. Searching for genomic signatures of selection, we identified significantly differentiated SNPs among which several are located in or within 20 kb distance from candidate genes for cold adaptation. These genes involve cellular and metabolic processes and several of them have been shown to be differentially expressed under diapausing conditions. The three biogeographical regions also differ for wing size and shape, and wing size increases with latitude supporting Bergmann's rule. Adaptive genetic and morphometric variation observed along the climatic gradient of A. albopictus native range suggests that colonization of northern latitudes promoted adaptation to cold environments prior to its worldwide invasion.     

Relationship of runs of homozygosity with adaptive and production traits in a paternal broiler line

Genomic regions under high selective pressure present specific runs of homozygosity (ROH), which provide valuable information on the genetic mechanisms underlying the adaptation to environment imposed challenges. In broiler chickens, the adaptation to conventional production systems in tropical environments lead the animals with favorable genotypes to be naturally selected, increasing the frequency of these alleles in the next generations. In this study, ~1400 chickens from a paternal broiler line were genotyped with the 600 K Affymetrix® Axiom® high-density (HD) genotyping array for estimation of linkage disequilibrium (LD), effective population size (N_e), inbreeding and ROH. The average LD between adjacent single nucleotide polymorphisms (SNPs) in all autosomes was 0.37, and the LD decay was higher in microchromosomes followed by intermediate and macrochromosomes. The N_e of the ancestral population was high and declined over time maintaining a sufficient number of animals to keep the inbreeding coefficient of this population at low levels. The ROH analysis revealed genomic regions that harbor genes associated with homeostasis maintenance and immune system mechanisms, which may have been selected in response to heat stress. Our results give a comprehensive insight into the relationship between shared ROH regions and putative regions related to survival and production traits in a paternal broiler line selected for over 20 years. These findings contribute to the understanding of the effects of environmental and artificial selection in shaping the distribution of functional variants in the chicken genome.     

Genome resequencing and bioinformatic analysis of SNP containing candidate genes in the autoimmune vitiligo Smyth line chicken model

Background

The Smyth line (SL) chicken is the only animal model for autoimmune vitiligo that spontaneously displays all clinical and biological manifestations of the human disorder. To understand the genetic components underlying the susceptibility to develop SL vitiligo (SLV), whole genome resequencing analysis was performed in SLV chickens compared with non-vitiliginous parental Brown line (BL) chickens, which maintain a very low incidence rate of vitiligo.

Results

Illumina sequencing technology and reference based assembly on Red Jungle Fowl genome sequences were used. Results of genome resequencing of pooled DNA of each 10 BL and SL chickens reached 5.1x and 7.0x coverage, respectively. The total number of SNPs was 4.8 and 5.5 million in BL and SL genome, respectively. Through a series of filtering processes, a total of ~1 million unique SNPs were found in the SL alone. Eventually of the 156 reliable marker SNPs, which can induce non-synonymous-, frameshift-, nonsense-, and no-start mutations in amino acid sequences in proteins, 139 genes were chosen for further analysis. Of these, 14 randomly chosen SNPs were examined for SNP verification by PCR and Sanger sequencing to detect SNP positions in 20 BL and 70 SL chickens. The results of the analysis of the 14 SNPs clearly showed differential frequencies of nucleotide bases in the SNP positions between BL and SL chickens. Bioinformatic analysis showed that the 156 most reliable marker SNPs included genes involved in dermatological diseases/conditions such as ADAMTS13, ASPM, ATP6V0A2, BRCA2, COL12A1, GRM5, LRP2, OBSCN, PLAU, RNF168, STAB2, and XIRP1. Intermolecular gene network analysis revealed that candidate genes identified in SLV play a role in networks centered on protein kinases (MAPK, ERK1/2, PKC, PRKDC), phosphatase (PPP1CA), ubiquitinylation (UBC) and amyloid production (APP).

Conclusions

Various potential genetic markers showing amino acid changes and potential roles in vitiligo development were identified in the SLV chicken through genome resequencing. The genetic markers and bioinformatic interpretations of amino acid mutations found in SLV chickens may provide insight into the genetic component responsible for the onset and the progression of autoimmune vitiligo and serve as valuable markers to develop diagnostic tools to detect vitiligo susceptibility.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-707) contains supplementary material, which is available to authorized users.     

Whole genome sequencing of Ethiopian highlanders reveals conserved hypoxia tolerance genes

Background

Although it has long been proposed that genetic factors contribute to adaptation to high altitude, such factors remain largely unverified. Recent advances in high-throughput sequencing have made it feasible to analyze genome-wide patterns of genetic variation in human populations. Since traditionally such studies surveyed only a small fraction of the genome, interpretation of the results was limited.

Results

We report here the results of the first whole genome resequencing-based analysis identifying genes that likely modulate high altitude adaptation in native Ethiopians residing at 3,500 m above sea level on Bale Plateau or Chennek field in Ethiopia. Using cross-population tests of selection, we identify regions with a significant loss of diversity, indicative of a selective sweep. We focus on a 208 kbp gene-rich region on chromosome 19, which is significant in both of the Ethiopian subpopulations sampled. This region contains eight protein-coding genes and spans 135 SNPs. To elucidate its potential role in hypoxia tolerance, we experimentally tested whether individual genes from the region affect hypoxia tolerance in Drosophila. Three genes significantly impact survival rates in low oxygen: cic, an ortholog of human CIC, Hsl, an ortholog of human LIPE, and Paf-AHα, an ortholog of human PAFAH1B3.

Conclusions

Our study reveals evolutionarily conserved genes that modulate hypoxia tolerance. In addition, we show that many of our results would likely be unattainable using data from exome sequencing or microarray studies. This highlights the importance of whole genome sequencing for investigating adaptation by natural selection.     

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.     

Identification of candidate genes related to highland adaptation from multiple Chinese local chicken breeds by whole genome sequencing analysis

Understanding the genetic mechanism of highland adaptation is of great importance for breeding improvement of Tibetan chickens (TBC). Some studies of TBC have identified some candidate genes and pathways from multiple subgroups, but the related genetic mechanisms remain largely unknown. Different genetic backgrounds and the independent genetic basis of highland adaptation make it difficult to identity the selective region of highland adaptation with all TBC samples. In this study, we conducted pre-analysis in a large-scale population to select a TBC subgroup with the purest and highest level of highland-specific lineage for the further analysis. Finally, the 37 samples from a TBC subgroup and 19 Lahsa White chickens were used to represent the highland group for further analysis with 80 samples from five Chinese local lowland breeds as controls. Population structure analysis revealed that highland adaptation significantly affected population stratification in Chinese local chicken breeds. Genome-wide selection signal analysis identified 201 candidate genes associated with highland adaptation of TBC, and these genes were significantly enriched in calcium signaling, vascular smooth muscle contraction and the cellular response to oxidative stress pathways. Additionally, we identified a narrow 1.76 kb region containing an overlapping region between HBZ and an active enhancer, and our identified region showed a highly significant signal. The highland group selected the haplotype with high activity to improve the oxygen-carrying capacity, thus being adapted to a hypoxic environment. We also found that STX2 was significantly selected in the highland group, thus potentially reducing the oxidative stress caused by hypoxia, and that STX2 exhibited the opposite effects on highland adaptation and reproductive traits. Our findings advance our understanding of extreme environment adaptation of highland chickens, and provide some variants and genes beneficial to TBC genetic breeding improvement.     

Polymorphisms of the interleukin-15 gene and their associations with fatness and muscle fiber traits in chickens

Interleukin-15 (IL-15) is a cytokine that has been proposed to modulate skeletal muscle and adipose tissue mass. In the present study, an F₂ resource population of Gushi chickens crossed with Anka broilers was used to investigate the genetic effects of the chicken IL-15 gene. Two single nucleotide polymorphisms (SNPs) (g.31224G>A and g.31266T>G) were identified in exon 5 of the IL-15 gene by means of polymerase chain reaction?Crestriction fragment length polymorphism (PCR-RFLP) and DNA sequencing. Associations between the two SNPs and chicken fatness and muscle fiber traits were determined using linkage disequilibrium, haplotype construction, and association analysis. Both of the SNPs were associated with abdominal fat weight, leg muscle fiber diameter, and leg muscle fiber density (p?<?0.05). Haplotypes of the two linked SNPs were associated with abdominal fat weight, fat thickness under the skin, and leg muscle fiber diameter (p?<?0.05). The results suggested that the IL-15 gene might be associated with the causative mutation or the quantitative trait locus (QTL) controlling the fatness traits and muscle fiber traits in chickens.     

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

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