Genetic mapping of the sex-linked barring gene in the chicken

Genomic insights into shank and eggshell color in Italian local chickens

Eggshell and shank color in poultry is an intriguing topic of research due to the roles in selection, breed recognition, and environmental adaptation. This study delves into the genomics foundations of shank and eggshell pigmentation in Italian local chickens through genome-wide association studies analysis to uncover the mechanisms governing these phenotypes. To this purpose, 483 animals from 20 local breeds (n = 466) and 2 commercial lines (n = 17) were considered and evaluated for shank and eggshell color. All animals were genotyped using the Affymetrix Axiom 600 K Chicken Genotyping Array. As regards shank color, the most interesting locus was detected on chromosome Z, close to the TYRP1 gene, known to play a key role in avian pigmentation. Additionally, several novel loci and genes associated with shank pigmentation, skin pigmentation, UV protection, and melanocyte regulation were identified (e.g., MTAP, CDKN2A, CDKN2B). In eggshell, fewer significant loci were identified, including SLC7A11 and MITF on chromosomes 4 and 12, respectively, associated with melanocyte processes and pigment synthesis. This comprehensive study shed light on the genetic architecture underlying shank and eggshell color in Italian native chicken breeds, contributing to a better understanding of this phenomenon which plays a role in breed identification and conservation, and has ecological and economic implications.

Genome-wide association study revealed the genomic regions associated with skin pigmentation in an Ogye x White Leghorn F2 chicken population

Skin color in chickens is an economically important trait that determines the first impression of a consumer toward a broiler and can ultimately affect consumer choice in the market. Therefore, identification of genomic regions associated with skin color is crucial for increasing the sales value of chickens. Although previous studies have attempted to reveal the genetic markers associated with the skin coloration in chickens, most were limited to investigations of candidate genes, such as melanin-related genes, and focused on case/control studies based on a single or small population. In this study, we performed a genome-wide association study (GWAS) on 770 F2 intercrosses produced by an experimental population of 2 chicken breeds, namely Ogye and White Leghorns, with different skin colors. The GWAS demonstrated that the L* value among the 3 skin color traits is highly heritable, and the genomic regions located on 2 chromosomes (20 and Z) were detected to harbor SNPs significantly associated with the skin color trait, accounting for most of the total genetic variance. Particular genomic regions spanning a ∼2.94 Mb region on GGA Z and a ∼3.58 Mb region on GGA 20 were significantly associated with skin color traits, and in these regions, certain candidate genes, including MTAP, FEM1C, GNAS, and EDN3, were found. Our findings could help elucidate the genetic mechanisms underlying chicken skin pigmentation. Furthermore, the candidate genes can be used to provide a valuable breeding strategy for the selection of specific chicken breeds with ideal skin coloration.

A genome-wide association study for eumelanin pigmentation in chicken plumage using a computer vision approach

Chicken plumage colouration is an important trait related to productivity in poultry industry. Therefore, the genetic basis for pigmentation in chicken plumage is an area of great interest. However, the colour trait is generally regarded as a qualitative trait and representing colour variations is difficult. In this study, we developed a method to quantify and classify colour using an F₂ population crossed from two pure lines: White Leghorn and the Korean indigenous breed Yeonsan Ogye. Using red, green, and blue values in the cropped body region, we identified significant genomic regions on chromosomes 33:3 160 480-7 447 197 and Z:78 748 287-79 173 793. Furthermore, we identified two potential candidate genes (PMEL and MTAP) that might have significant effects on melanin-based plumage pigmentation. Our study presents a new phenotyping method using a computer vision approach and provides new insights into the genetic basis of melanin-based feather colouration in chickens.

Embryonic transcriptome unravels mechanisms and pathways underlying embryonic development with respect to muscle growth, egg production, and plumage formation in native and broiler chickens

Background: Muscle development, egg production, and plumage colors are different between native and broiler chickens. The study was designed to investigate why improved Aseel (PD4) is colorful, stronger, and grew slowly compared with the control broiler (CB).

Methods: A microarray was conducted using the 7th-day embryo (7EB) and 18th-day thigh muscle (18TM) of improved Aseel and broiler, respectively. Also, we have selected 24 Gallus gallus candidate reference genes from NCBI, and total RNA was isolated from the broiler, improved Aseel embryo tissues, and their expression profiles were studied by real-time quantitative PCR (qPCR). Furthermore, microarray data were validated with qPCR using improved Aseel and broiler embryo tissues.

Results: In the differential transcripts screening, all the transcripts obtained by microarray of slow and fast growth groups were screened by fold change ≥ 1 and false discovery rate (FDR) ≤ 0.05. In total, 8,069 transcripts were differentially expressed between the 7EB and 18TM of PD4 compared to the CB. A further analysis showed that a high number of transcripts are differentially regulated in the 7EB of PD4 (6,896) and fewer transcripts are differentially regulated (1,173) in the 18TM of PD4 compared to the CB. On the 7th^- and 18th-day PD4 embryos, 3,890, 3,006, 745, and 428 transcripts were up- and downregulated, respectively. The commonly up- and downregulated transcripts are 91 and 44 between the 7th- and 18th-day of embryos. In addition, the best housekeeping gene was identified. Furthermore, we validated the differentially expressed genes (DEGs) related to muscle growth, myostatin signaling and development, and fatty acid metabolism genes in PD4 and CB embryo tissues by qPCR, and the results correlated with microarray expression data.

Conclusion: Our study identified DEGs that regulate the myostatin signaling and differentiation pathway; glycolysis and gluconeogenesis; fatty acid metabolism; Jak-STAT, mTOR, and TGF-β signaling pathways; tryptophan metabolism; and PI3K-Akt signaling pathways in PD4. The results revealed that the gene expression architecture is present in the improved Aseel exhibiting embryo growth that will help improve muscle development, differentiation, egg production, protein synthesis, and plumage formation in PD4 native chickens. Our findings may be used as a model for improving the growth in Aseel as well as optimizing the growth in the broiler.

Population Genomics of Variegated Toad-Headed Lizard Phrynocephalus versicolor and Its Adaptation to the Colorful Sand of the Gobi Desert

The variegated toad-headed agama, Phrynocephalus versicolor, lives in the arid landscape of the Chinese Gobi Desert. We analyzed populations from three different locations which vary in substrate color and altitude: Heishankou (HSK), Guazhou County (GZ), and Ejin Banner (EJN). The substrate color is either light-yellow (GZ-y), yellow (EJN-y), or black (HSK-b); the corresponding lizard population colors largely match their substrate in the degree of melanism. We assembled the P. versicolor genome and sequenced over 90 individuals from the three different populations. Genetic divergence between populations corresponds to their geographic distribution. We inferred the genetic relationships among these populations and used selection scans and differential expression to identify genes that show signatures of selection. Slc2a11 and akap12, among other genes, are highly differentiated and may be responsible for pigment adaptation to substrate color in P. versicolor.

Transcript expression profiling of fibromelanosis-related genes in black-bone chickens

1. The aim of the present study was to identify differentially expressed genes (DEGs) and metabolic pathways involved in this phenotype. Fibromelanosis is the most striking feature of black-bone chickens, such as the Silkie and Dongxiang indigenous breeds. Due to the accumulation of eumelanin in connective tissues, fibromelanosis manifests as black colouration of the skin, muscles, gut, and periosteum. Studies on fibromelanosis can provide useful information pertaining to human diseases and offer commercial value to the poultry industry. However, the genetic basis of fibromelanosis remains unclear.2. Digital gene expression analysis was performed on black and white skin samples collected from the HW1 black-bone chicken line to detect differences in genome-wide expression patterns. A total of >30 billion bp were sequenced, and 2,707,926,466 bp and 2,948,782,964 bp of clean data obtained for creation of libraries for black and white skin, respectively. In total, 252 DEGs from 15,508 mapped genes were identified with 83 up-regulated in white skin and 169 up-regulated in black skin.3. Gene ontology analysis highlighted that genes from the extracellular region and associated components were abundant among the DEGs. Pathway analysis revealed that many DEGs were linked to amino acid metabolism and the immune system. qRT-PCR validation using 14 genes showed good conformity with the sequence analysis of fibromelanosis-related genes.4. The results showed that L-dopachrometautomerase precursor (DCT), tyrosine aminotransferase (TAT), 4-hydroxyphenylpyruvate dioxygenase (HPD) from the tyrosine metabolism pathway, coagulation factor II (F2), fibrinogen beta chain (FGB), plasminogen (PLG) and complement component 7 (C7) from the complement and coagulation cascades were important genes in the fibromelanosis process in black-bone chickens. These candidate genes require further correlation analysis and functional verification.

Genetic basis of chicken plumage color in artificial population of complex epistasis

Chicken plumage color, the genetic basis of which is often affected by epistasis, has long interested scientists. In the current study, a population of complex epistasis was constructed by crossing dominant White Leghorn chickens with recessive white feather chickens. Through a genome-wide association study, we identified single nucleotide polymorphisms and genes significantly associated with white and colored plumage in hens at different developmental stages. Interestingly, white plumage in adulthood was associated with the recessive white feather gene (TYR), whereas white feathers at birth stage were associated with the dominant white feather gene (PMEL), indicating age-related roles for these genes. TYR was shown to exert an epistatic effect on PMEL in adult hens. Additionally, TYR had an epistatic effect on barred plumage, while barred plumage had an epistatic effect on black plumage. TYR had no epistatic effect on the yellow plumage. We confirmed that the barred plumage gene is CDKN2A, as reported in previous studies. Golgb1 and REEP3, which play important roles in the Golgi network and affect the formation of feather pigments, are important candidate genes for yellow plumage. The candidate genes for black plumage are CAMKK1 and IFT22. Further research is warranted to elucidate the molecular mechanisms underlying these traits.

Consistency of dark skeletal muscles in Thai native black-bone chickens (Gallus gallus domesticus)

Black-bone chickens (Gallus gallus domesticus) have become economically valuable, particularly in Southeast Asia as a consequence of popular traditional Chinese medical practices. Chickens with whole body organ darkness are considered to have higher value and are, therefore, more often requested. This research study aimed to investigate the darkness in 34 skeletal muscles of 10 Thai black-bone chickens (five males and five females). The evaluation of muscle darkness was done on two levels: (i) a color chart was employed at the macroanatomical level and (ii) by using melanin pigment to evaluate the structure at the microanatomy level. The results revealed that the accumulation of melanin pigment in the muscle tissue was observed in the endomysium, perimysium and epimysium. With respect to the results of the color chart test, iliotibialis lateralis pars preacetabularis, gastrocnemius, fibularis longus and puboischiofemoralis pars medialis showed the highest degree of darkness, while serratus profundus, pectoralis, iliotibialis cranialis, flexor cruris lateralis, and flexor cruris medialis appeared to be the least dark. In addition, we found that the highest and lowest amounts of melanin pigment was noted in the flexor carpi ulnaris and pectoralis (p < 0.05), respectively; however, there was no significant difference (p > 0.05) observed between the sexes. These results reveal that the 34 specified muscles of black-bone chickens showed uneven distribution of darkness due to the differing accumulations of melanin pigments of each muscle.This information may provide background knowledge for a better understanding of melanin accumulation and lead to breeding improvements in Thai black-bone chickens.

Identification and Differential Abundance of Mitochondrial Genome Encoding Small RNAs (mitosRNA) in Breast Muscles of Modern Broilers and Unselected Chicken Breed

Background: Although small non-coding RNAs are mostly encoded by the nuclear genome, thousands of small non-coding RNAs encoded by the mitochondrial genome, termed as mitosRNAs were recently reported in human, mouse and trout. In this study, we first identified chicken mitosRNAs in breast muscle using small RNA sequencing method and the differential abundance was analyzed between modern pedigree male (PeM) broilers (characterized by rapid growth and large muscle mass) and the foundational Barred Plymouth Rock (BPR) chickens (characterized by slow growth and small muscle mass). Methods: Small RNA sequencing was performed with total RNAs extracted from breast muscles of PeM and BPR (n = 6 per group) using the 1 × 50 bp single end read method of Illumina sequencing. Raw reads were processed by quality assessment, adapter trimming, and alignment to the chicken mitochondrial genome (GenBank Accession: X52392.1) using the NGen program. Further statistical analyses were performed using the JMP Genomics 8. Differentially expressed (DE) mitosRNAs between PeM and BPR were confirmed by quantitative PCR. Results: Totals of 183,416 unique small RNA sequences were identified as potential chicken mitosRNAs. After stringent filtering processes, 117 mitosRNAs showing >100 raw read counts were abundantly produced from all 37 mitochondrial genes (except D-loop region) and the length of mitosRNAs ranged from 22 to 46 nucleotides. Of those, abundance of 44 mitosRNAs were significantly altered in breast muscles of PeM compared to those of BPR: all mitosRNAs were higher in PeM breast except those produced from 16S-rRNA gene. Possibly, the higher mitosRNAs abundance in PeM breast may be due to a higher mitochondrial content compared to BPR. Our data demonstrate that in addition to 37 known mitochondrial genes, the mitochondrial genome also encodes abundant mitosRNAs, that may play an important regulatory role in muscle growth via mitochondrial gene expression control.

RNA sequencing for global gene expression associated with muscle growth in a single male modern broiler line compared to a foundational Barred Plymouth Rock chicken line

Background

Modern broiler chickens exhibit very rapid growth and high feed efficiency compared to unselected chicken breeds. The improved production efficiency in modern broiler chickens was achieved by the intensive genetic selection for meat production. This study was designed to investigate the genetic alterations accumulated in modern broiler breeder lines during selective breeding conducted over several decades.

Methods

To identify genes important in determining muscle growth and feed efficiency in broilers, RNA sequencing (RNAseq) was conducted with breast muscle in modern pedigree male (PeM) broilers (n = 6 per group), and with an unselected foundation broiler line (Barred Plymouth Rock; BPR). The RNAseq analysis was carried out using Ilumina Hiseq (2 x 100 bp paired end read) and raw reads were assembled with the galgal4 reference chicken genome. With normalized RPM values, genes showing >10 average read counts were chosen and genes showing <0.05 p-value and >1.3 fold change were considered as differentially expressed (DE) between PeM and BPR. DE genes were subjected to Ingenuity Pathway Analysis (IPA) for bioinformatic functional interpretation.

Results

The results indicate that 2,464 DE genes were identified in the comparison between PeM and BPR. Interestingly, the expression of genes encoding mitochondrial proteins in chicken are significantly biased towards the BPR group, suggesting a lowered mitochondrial content in PeM chicken muscles compared to BPR chicken. This result is inconsistent with more slow muscle fibers bearing a lower mitochondrial content in the PeM. The molecular, cellular and physiological functions of DE genes in the comparison between PeM and BPR include organismal injury, carbohydrate metabolism, cell growth/proliferation, and skeletal muscle system development, indicating that cellular mechanisms in modern broiler lines are tightly associated with rapid growth and differential muscle fiber contents compared to the unselected BPR line. Particularly, PDGF (platelet derived growth factor) signaling and NFE2L2 (nuclear factor, erythroid 2-like 2; also known as NRF2) mediated oxidative stress response pathways appear to be activated in modern broiler compared to the foundational BPR line. Upstream and network analyses revealed that the MSTN (myostatin) –FST (follistatin) interactions and inhibition of AR (androgen receptor) were predicted to be effective regulatory factors for DE genes in modern broiler line. PRKAG3 (protein kinase, AMP-activated, gamma 3 non-catalytic subunit) and LIPE (lipase E) are predicted as core regulatory factors for myogenic development, nutrient and lipid metabolism.

Conclusion

The highly upregulated genes in PeM may represent phenotypes of subclinical myopathy commonly observed in the commercial broiler breast tissue, that can lead to muscle hardening, named as woody breast. By investigating global gene expression in a highly selected pedigree broiler line and a foundational breed (Barred Plymouth Rock), the results provide insight into cellular mechanisms that regulate muscle growth, fiber composition and feed efficiency.

Electronic supplementary material

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

A genome-wide association study identifies novel single nucleotide polymorphisms associated with dermal shank pigmentation in chickens

Shank color of domestic chickens varies from black to blue, green, yellow, or white, which is controlled by the combination of melanin and xanthophylls in dermis and epidermis. Dermal shank pigmentation of chickens is determined by sex-linked inhibitor of dermal melanin (Id), which is located on the distal end of the long arm of Z chromosome, through controlling dermal melanin pigmentation. Although previous studies have focused on the identification of Id and the linear relationship with barring and recessive white skin, no causal mutations have yet been identified in relation to the mutant dermal pigment inhibiting allele at the Id locus. In this study, we first used the 600K Affymetrix Axiom HD genotyping array, which includes ~580,961 SNP of which 26,642 SNP were on the Z chromosome to perform a genome-wide association study on pure lines of 19 Tibetan hens with dermal pigmentation shank and 21 Tibetan hens with yellow shank to refine the Id location. Association analysis was conducted by the PLINK software using the standard chi-squared test, and then Bonferroni correction was used to adjust multiple testing. The genome-wide study revealed that 3 SNP located at 78.5 to 79.2 Mb on the Z chromosome in the current assembly of chicken genome (galGal4) were significantly associated with dermal shank pigmentation of chickens, but none of them were located in known genes. The interval we refined was partly converged with previous results, suggesting that the Id gene is in or near our refined genome region. However, the genomic context of this region was complex. There were only 15 SNP markers developed by the genotyping array within the interval region, in which only 1 SNP marker passed quality control. Additionally, there were about 5.8-Mb gaps on both sides of the refined interval. The follow-up replication studies may be needed to further confirm the functional significance for these newly identified SNP.

Gene duplication of endothelin 3 is closely correlated with the hyperpigmentation of the internal organs (Fibromelanosis) in silky chickens

During early development in vertebrates, pluripotent cells are generated from the neural crest and migrate according to their presumptive fate. In birds and mammals, one of the progeny cells, melanoblasts, generally migrate through a dorsolateral route of the trunk region and differentiate to melanocytes. However, Silky is an exceptional chicken in which numerous melanoblasts travel via a ventral pathway and disperse into internal organs. Finally, these ectopic melanocytes induce heavy dermal and visceral melanization known as Fibromelanosis (Fm). To identify the genetic basis of this phenotype, we confirmed the mode of inheritance of Fm as autosomal dominant and then performed linkage analysis with microsatellite markers and sequence-tagged site markers. Using 85 backcross progeny from crossing Black Minorca chickens (BM-C) with F(1) individuals between White Silky (WS) and BM-C Fm was located on 10.2-11.7 Mb of chicken chromosome 20. In addition, we noticed a DNA marker that all Silky chickens and the F(1) individuals showed heterozygous genotyping patterns, suggesting gene duplication in the Fm region. By quantitative real-time PCR assay, Silky line-specific gene duplication was detected as an ~130-kb interval. It contained five genes including endothelin 3 (EDN3), which encoded a potent mitogen for melanoblasts/melanocytes. EDN3 with another three of these duplicated genes in Silky chickens expressed almost twofold of those in BM-C. Present results strongly suggest that the increase of the expression levels resulting from the gene duplication in the Fm region is the trigger of hypermelanization in internal organs of Silky chickens.

Poultry genome sequences: progress and outstanding challenges

The first build of the chicken genome sequence appeared in March, 2004 - the first genome sequence of any animal agriculture species. That sequence was done primarily by whole genome shotgun Sanger sequencing, along with the use of an extensive BAC contig-based physical map to assemble the sequence contigs and scaffolds and align them to the known chicken chromosomes and linkage groups. Subsequent sequencing and mapping efforts have improved upon that first build, and efforts continue in search of missing and/or unassembled sequence, primarily on the smaller microchromosomes and the sex chromosomes. In the past year, a draft turkey genome sequence of similar quality has been obtained at a much lower cost primarily due to the development of 'next-generation' sequencing techniques. However, assembly and alignment of the sequence contigs and scaffolds still depended on a detailed BAC contig map of the turkey genome that also utilized comparison to the existing chicken sequence. These 2 land fowl (Galliformes) genomes show a remarkable level of similarity, despite an estimated 30-40 million years of separate evolution since their last common ancestor. Among the advantages offered by these sequences are routine re-sequencing of commercial and research lines to identify the genetic correlates of phenotypic change (for example, selective sweeps), a much improved understanding of poultry diversity and linkage disequilibrium, and access to high-density SNP typing and association analysis, detailed transcriptomic and proteomic studies, and the use of genome-wide marker- assisted selection to enhance genetic gain in commercial stocks.