Congjiao Sun

Genome-Wide Characterization of Insertion and Deletion Variation in Chicken Using Next Generation Sequencing

Insertion and deletion (INDEL) is one of the main events contributing to genetic and phenotypic diversity, which receives less attention than SNP and large structural variation. To gain a better knowledge of INDEL variation in chicken genome, we applied next generation sequencing on 12 diverse chicken breeds at an average effective depth of 8.6. Over 1.3 million non-redundant short INDELs (1–49 bp) were obtained, the vast majority (92.48%) of which were novel. Follow-up validation assays confirmed that most (88.00%) of the randomly selected INDELs represent true variations. The majority (95.76%) of INDELs were less than 10 bp. Both the detected number and affected bases were larger for deletions than insertions. In total, INDELs covered 3.8 Mbp, corresponding to 0.36% of the chicken genome. The average genomic INDEL density was estimated as 0.49 per kb. INDELs were ubiquitous and distributed in a non-uniform fashion across chromosomes, with lower INDEL density in micro-chromosomes than in others, and some functional regions like exons and UTRs were prone to less INDELs than introns and intergenic regions. Nearly 620,253 INDELs fell in genic regions, 1,765 (0.28%) of which located in exons, spanning 1,358 (7.56%) unique Ensembl genes. Many of them are associated with economically important traits and some are the homologues of human disease-related genes. We demonstrate that sequencing multiple individuals at a medium depth offers a promising way for reliable identification of INDELs. The coding INDELs are valuable candidates for further elucidation of the association between genotypes and phenotypes. The chicken INDELs revealed by our study can be useful for future studies, including development of INDEL markers, construction of high density linkage map, INDEL arrays design, and hopefully, molecular breeding programs in chicken.

Gut metagenomic analysis reveals prominent roles of Lactobacillus and cecal microbiota in chicken feed efficiency

Interactions between the host and gut microbiota can affect gut metabolism. In this study, the individual performances of 252 hens were recorded to evaluate feed efficiency. Hens with contrasting feed efficiencies (14 birds per group) were selected to investigate their duodenal, cecal and fecal microbial composition by sequencing the 16S rRNA gene V4 region. The results showed that the microbial community in the cecum was quite different from those in the duodenum and feces. The highest biodiversity and all differentially abundant taxa between the different efficiency groups were observed in the cecal microbial community with false discovery rate (FDR) <0.05. Of these differentially abundant cecal microbes, Lactobacillus accounted for a greater proportion than the others. The abundances of Lactobacillus and Akkermansia were significantly higher while that of Faecalibacterium was lower (FDR < 0.05) in the better feed efficiency (BFE) group. Phylogenetic investigation of communities by reconstruction of unobserved states (PICRUSt) analysis revealed that the functions relating to glycometabolism and amino acid metabolism were enriched in the cecal microbiota of the BFE group. These results indicated the prominent role of cecal microbiota in the feed efficiency of chickens and suggested plausible uses of Lactobacillus to improve the feed efficiency of host. Supplementary information The online version of this article (doi:10.1038/srep45308) contains supplementary material, which is available to authorized users.

Identification of Promising Mutants Associated with Egg Production Traits Revealed by Genome-Wide Association Study

Egg number (EN), egg laying rate (LR) and age at first egg (AFE) are important production traits related to egg production in poultry industry. To better understand the knowledge of genetic architecture of dynamic EN during the whole laying cycle and provide the precise positions of associated variants for EN, LR and AFE, laying records from 21 to 72 weeks of age were collected individually for 1,534 F2 hens produced by reciprocal crosses between White Leghorn and Dongxiang Blue-shelled chicken, and their genotypes were assayed by chicken 600 K Affymetrix high density genotyping arrays. Subsequently, pedigree and SNP-based genetic parameters were estimated and a genome-wide association study (GWAS) was conducted on EN, LR and AFE. The heritability estimates were similar between pedigree and SNP-based estimates varying from 0.17 to 0.36. In the GWA analysis, we identified nine genome-wide significant loci associated with EN of the laying periods from 21 to 26 weeks, 27 to 36 weeks and 37 to 72 weeks. Analysis of GTF2A1 and CLSPN suggested that they influenced the function of ovary and uterus, and may be considered as relevant candidates. The identified SNP rs314448799 for accumulative EN from 21 to 40 weeks on chromosome 5 created phenotypic differences of 6.86 eggs between two homozygous genotypes, which could be potentially applied to the molecular breeding for EN selection. Moreover, our finding showed that LR was a moderate polygenic trait. The suggestive significant region on chromosome 16 for AFE suggested the relationship between sex maturity and immune in the current population. The present study comprehensively evaluates the role of genetic variants in the development of egg laying. The findings will be helpful to investigation of causative genes function and future marker-assisted selection and genomic selection in chickens.

Promising Loci and Genes for Yolk and Ovary Weight in Chickens Revealed by a Genome-Wide Association Study.

Because it serves as the cytoplasm of the oocyte and provides a large amount of reserves, the egg yolk has biological significance for developing embryos. The ovary and its hierarchy of follicles are the main reproductive organs responsible for yolk deposition in chickens. However, the genetic architecture underlying the yolk and ovarian follicle weights remains elusive. Here, we measured the yolk weight (YW) at 11 age points from onset of egg laying to 72 weeks of age and measured the follicle weight (FW) and ovary weight (OW) at 73 weeks as part of a comprehensive genome-wide association study (GWAS) in 1,534 F2 hens derived from reciprocal crosses between White Leghorn (WL) and Dongxiang chickens (DX). For all ages, YWs exhibited moderate single nucleotide polymorphism (SNP)-based heritability estimates (0.25–0.38), while the estimates for FW (0.16) and OW (0.20) were relatively low. Independent univariate genome-wide screens for each trait identified 12, 3, and 31 novel significant associations with YW, FW, and OW, respectively. A list of candidate genes such as ZAR1, STARD13, ACER1b, ACSBG2, and DHRS12 were identified for having a plausible function in yolk and follicle development. These genes are important to the initiation of embryogenesis, lipid transport, lipoprotein synthesis, lipid droplet promotion, and steroid hormone metabolism, respectively. Our study provides for the first time a genome-wide association (GWA) analysis for follicle and ovary weight. Identification of the promising loci as well as potential candidate genes will greatly advance our understanding of the genetic basis underlying dynamic yolk weight and ovarian follicle development and has practical significance in breeding programs for the alteration of yolk weight at different age points.

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.

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.

Polymorphisms in Ion Transport Genes Are Associated with Eggshell Mechanical Property

Eggshell mechanical property traits such as eggshell breaking strength (ESS), eggshell thickness (EST) and eggshell weight (ESW) are most common and important indexes to evaluate eggshell quality in poultry industry. Uterine ion transporters involve in eggshell formation and might be associated with eggshell mechanical property traits. In this study, 99 SNPs in 15 ion transport genes were selected to genotype 976 pedigreed hens of Rhode Island Red. ESS, EST and ESW were measured for each bird at 55 weeks of age. The association study showed that 14 SNPs in 8 genes were significantly related (p < 0.05) with at least one trait, and their contributions to phenotypic variance ranged from 0.23% to 4.14%. Both ATP2A3 and SLC4A5 had a significant effect on all the three traits. Strong linkage disequilibrium (LD) was detected among SNPs in four genes: ATP2A3, ITPR1, SLC8A3, SCNN1a. The significant effects of those diplotypes on eggshell mechanical property traits were found, and their contributions to phenotypic variance ranged from 0.50% to 0.70%. It was concluded that the identified SNPs and diplotypes in this study were potential markers influencing the eggshell mechanical properties, which could contribute to the genetic improvement of eggshell quality.

Genetic architecture dissection by genome-wide association analysis reveals avian eggshell ultrastructure traits

The ultrastructure of an eggshell is considered the major determinant of eggshell quality, which has biological and economic significance for the avian and poultry industries. However, the interrelationships and genome-wide architecture of eggshell ultrastructure remain to be elucidated. Herein, we measured eggshell thickness (EST), effective layer thickness (ET), mammillary layer thickness (MT), and mammillary density (MD) and conducted genome-wide association studies in 927 F2 hens. The SNP-based heritabilities of eggshell ultrastructure traits were estimated to be 0.39, 0.36, 0.17 and 0.19 for EST, ET, MT and MD, respectively, and a total of 719, 784, 1 and 10 genome-wide significant SNPs were associated with EST, ET, MT and MD, respectively. ABCC9, ITPR2, KCNJ8 and WNK1, which are involved in ion transport, were suggested to be the key genes regulating EST and ET. ITM2C and KNDC1 likely affect MT and MD, respectively. Additionally, there were linear relationships between the chromosome lengths and the variance explained per chromosome for EST (R2 = 0.57) and ET (R2 = 0.67). In conclusion, the interrelationships and genetic architecture of eggshell ultrastructure traits revealed in this study are valuable for our understanding of the avian eggshell and contribute to research on a variety of other calcified shells.

Back Cover: Divergent Proteome Patterns of Egg Albumen from Domestic Chicken, Duck, Goose, Turkey, Quail and Pigeon

To uncover a diversity of genetic and biological unknowns, a comprehensive and comparative proteomic analysis is performed on egg albumen of domestic chicken, duck, goose, turkey, quail, and pigeon with tandem mass tags quantification technology. In this study, a total of 148, 138, 150, 162, 183, and 179 proteins are identified in egg albumen of the above six species, respectively. Venn plots, PCA, and cluster analysis all reveal the highest similarity of protein composition between duck and goose (≈75%). Additionally, the six species have 52 proteins detected in common in the egg albumen. As revealed by GO and pathway analyses, the plausible functions of these highly conserved proteins are to provide a secure environment and prevent the early death of embryonic cells. Species‐specific proteins such as haptoglobin in pigeon, serpin‐like protein HMSD in duck, and ovodefensin in chicken are also screened and are likely associated with species‐dependent biological processes. Furthermore, Enzyme Code analysis indicated egg albumen have abundant enzyme activity, with hydrolases accounting for more than half of the total enzymes. This study is the first to provide the proteome profiles of egg albumen for the major poultry species, which will be instructive for the understanding of species‐specific biological problems with egg albumen.

iTRAQ-Based Quantitative Proteomics Identifies Potential Regulatory Proteins Involved in Chicken Eggshell Brownness.

Brown eggs are popular in many countries and consumers regard eggshell brownness as an important indicator of egg quality. However, the potential regulatory proteins and detailed molecular mechanisms regulating eggshell brownness have yet to be clearly defined. In the present study, we performed quantitative proteomics analysis with iTRAQ technology in the shell gland epithelium of hens laying dark and light brown eggs to investigate the candidate proteins and molecular mechanisms underlying variation in chicken eggshell brownness. The results indicated 147 differentially expressed proteins between these two groups, among which 65 and 82 proteins were significantly up-regulated in the light and dark groups, respectively. Functional analysis indicated that in the light group, the down-regulated iron-sulfur cluster assembly protein (Iba57) would decrease the synthesis of protoporphyrin IX; furthermore, the up-regulated protein solute carrier family 25 (mitochondrial carrier; adenine nucleotide translocator), member 5 (SLC25A5) and down-regulated translocator protein (TSPO) would lead to increased amounts of protoporphyrin IX transported into the mitochondria matrix to form heme with iron, which is supplied by ovotransferrin protein (TF). In other words, chickens from the light group produce less protoporphyrin IX, which is mainly used for heme synthesis. Therefore, the exported protoporphyrin IX available for eggshell deposition and brownness is reduced in the light group. The current study provides valuable information to elucidate variation of chicken eggshell brownness, and demonstrates the feasibility and sensitivity of iTRAQ-based quantitative proteomics analysis in providing useful insights into the molecular mechanisms underlying brown eggshell pigmentation.