Genome-wide transcriptional profiling and functional analysis of long noncoding RNAs and mRNAs in chicken macrophages associated with the infection of avian pathogenic E. coli.

BACKGROUND: Avian pathogenic E. coli (APEC) can cause localized or systemic infections, collectively known as avian colibacillosis, resulting in huge economic losses to poultry industry globally per year. In addition, increasing evidence indicates that long non-coding RNAs (lncRNAs) play a critical role in regulating host inflammation in response to bacterial infection. However, the role of lncRNAs in the host response to APEC infection remains unclear. RESULTS: Here, we found 816 differentially expressed (DE) lncRNAs and 1,798 DE mRNAs in APEC infected chicken macrophages by RNAseq. The identified DE lncRNA-mRNAs were involved in Toll like receptor signaling pathway, VEGF signaling pathway, fatty acid metabolism, phosphatidylinositol signaling system, and other types of O-glycan biosynthesis. Furthermore, we found the novel lncRNA TCONS_00007391 as an important immune regulator in APEC infection was able to regulate the inflammatory response by directly targeting CD86. CONCLUSION: These findings provided a better understanding of host response to APEC infection and also offered the potential drug targets for therapy development against APEC infection.

DNA methylome and transcriptome identified Key genes and pathways involved in Speckled Eggshell formation in aged laying hens.

BACKGROUND: The quality of poultry eggshells is closely related to the profitability of egg production. Eggshell speckles reflect an important quality trait that influences egg appearance and customer preference. However, the mechanism of speckle formation remains poorly understood. In this study, we systematically compared serum immune and antioxidant indices of hens laying speckled and normal eggs. Transcriptome and methylome analyses were used to elucidate the mechanism of eggshell speckle formation. RESULTS: The results showed that seven differentially expressed genes (DEGs) were identified between the normal and speckle groups. Gene set enrichment analysis (GSEA) revealed that the expressed genes were mainly enriched in the calcium signaling pathway, focal adhesion, and MAPK signaling pathway. Additionally, 282 differentially methylated genes (DMGs) were detected, of which 15 genes were associated with aging, including ARNTL, CAV1, and GCLC. Pathway analysis showed that the DMGs were associated with T cell-mediated immunity, response to oxidative stress, and cellular response to DNA damage stimulus. Integrative analysis of transcriptome and DNA methylation data identified BFSP2 as the only overlapping gene, which was expressed at low levels and hypomethylated in the speckle group. CONCLUSIONS: Overall, these results indicate that aging- and immune-related genes and pathways play a crucial role in the formation of speckled eggshells, providing useful information for improving eggshell quality.

Genome-wide association study exploring the genetic architecture of eggshell speckles in laying hens.

BACKGROUND: Eggshell speckle phenotype is an important trait in poultry production because they affect eggshell quality. However, the genetic architecture of speckled eggshells remains unclear. In this study, we determined the heritability of eggshell speckles and conducted a genome-wide association study (GWAS) on purebred Rhode Island Red (RIR) hens at 28 weeks to detect potential genomic loci and candidate genes associated with eggshell speckles. RESULTS: The heritability of eggshell speckles was 0.35 at 28 weeks, and the speckle level is not related to other eggshell quality traits in terms of phenotypic correlation. We detected 311 SNPs (6 significantly, and 305 suggestively associated) and 39 candidate genes associated with eggshell speckles. Based on the pathway analysis, the 39 candidate genes were mainly involved in alpha-linolenic acid metabolism, linoleic acid metabolism, ether lipid metabolism, GnRH signaling pathway, vascular smooth muscle contraction, and MAPK signaling pathway. Ultimately, ten genes, LOC423226, SPTBN5, EHD4, LOC77155, TYRO3, ITPKA, DLL4, PLA2G4B, PLA2G4EL5, and PLA2G4EL6 were considered the most promising genes associated with eggshell speckles that were implicated in immunoregulation, calcium transport, and phospholipid metabolism, while its function in laying hens requires further studies. CONCLUSIONS: This study provides new insights into understanding the genetic basis of eggshell speckles and has practical application value for the genetic improvement of eggshell quality.

A ∼ 4.1 kb deletion in IRX1 gene upstream is completely associated with rumplessness in Piao chicken.

Piao chicken, a Chinese indigenous rumpless chicken breed, lacks pygostyle, caudal vertebra, uropygial gland and tail feathers. The rumplessness in Piao chicken presents an autosomal dominant inheritance pattern. However, the molecular genetic mechanisms underlying the rumplessness in Piao chicken remains unclear. In this study, whole-genome resequencing was performed for 146 individuals from 10 chicken breeds, including 9 tailed chicken breeds and Piao rumpless breed. Tailbone CT scan for Piao chickens and WL chickens, revealed that some Piao chicken tails were normal in number, and for a few Piao chickens tail length and tail bone numbers were between the rumpless and the normal tailed chickens. The results showed that the rumpless phenotype has not been completely fixed in Piao chicken breed. Using selection signature analysis and structural variation detection, we found a 4174 bp deletion located in the upstream region of IRX1 gene on chromosome 2 related to rumpless phenotype. Structural variation genotyping showed that the deletion was present in all 32 rumpless Piao chickens (del/del, wild/del) and absent from all 112 tailed chickens included in the dataset for the other 9 breeds and 2 tailed Piao chickens (wild/wild). In summary, all rumpless Piao chickens tested here carry this deletion mutation, to show a complete linkage association with rumplessness trait. We suggested that the 4174 bp deletion could be causative for rumpless phenotype in Piao chicken since this is the only mutation to show the complete linkage disequilibrium with rumplessness on whole genome level across all of 146 chickens from the 10 breeds. This study could facilitate a better understanding of the genetic characteristics of Piao chicken.

Genomic and transcriptomic analyses reveal genetic adaptation to cold conditions in the chickens.

Under the pressure of natural and artificial selection, domestic animals, including chickens, have evolved unique mechanisms of genetic adaptations such as high-altitude adaptation, hot and arid climate adaptation, and desert adaptation. Here, we investigated the genetic basis of cold tolerance in chicken by integrating whole-genome and transcriptome sequencing technologies. Genome-wide comparative analyses of 118 chickens living in different latitudes showed 46 genes and several pathways that may be involved in cold adaptation. The results of the functional enrichment analysis of differentially expressed genes proved the important role of metabolic pathways and immune-related pathways in cold tolerance in chickens. The subsequent integration of whole genome and transcriptome sequencing technology further identified six genes - dnah5 (dynein axonemal heavy chain 5), ptgs2 (prostaglandin-endoperoxide synthase 2), inhba (inhibin beta A subunit), irx2 (iroquois homeobox 2), ensgalg00000054917, and ensgalg00000046652 - requiring more detailed studies. In addition, we also discovered different allele frequency distributions of five SNPs (single nucleotide polymorphisms) within ptgs2 and nine SNPs within dnah5 in chickens in different latitudes, suggesting strong selective pressure of these two genes in chickens. We provide a novel insight into the genetic adaptation in chickens to cold environments, and provide a reference for evaluating and developing adaptive chicken breeds in cold environments.

Systematic Selection Signature Analysis of Chinese Gamecocks Based on Genomic and Transcriptomic Data.

Selection pressures driven by natural causes or human interference are key factors causing genome variants and signatures of selection in specific regions of the genome. Gamecocks were bred for cockfighting, presenting pea-combs, larger body sizes, stronger limbs, and higher levels of aggression than other chickens. In this study, we aimed to explore the genomic differences between Chinese gamecocks and commercial, indigenous, foreign, and cultivated breeds by detecting the regions or sites under natural or artificial selection using genome-wide association studies (GWAS), genome-wide selective sweeps based on the genetic differentiation index (FST), and transcriptome analyses. Ten genes were identified using GWAS and FST: gga-mir-6608-1, SOX5, DGKB, ISPD, IGF2BP1, AGMO, MEOX2, GIP, DLG5, and KCNMA1. The ten candidate genes were mainly associated with muscle and skeletal development, glucose metabolism, and the pea-comb phenotype. Enrichment analysis results showed that the differentially expressed genes between the Luxi (LX) gamecock and Rhode Island Red (RIR) chicken were mainly related to muscle development and neuroactive-related pathways. This study will help to understand the genetic basis and evolution of Chinese gamecocks and support the further use of gamecocks as an excellent breeding material from a genetic perspective.

Gene expression profiling of RIP2-knockdown in HD11 macrophages - elucidation of potential pathways (gene network) when challenged with avian pathogenic E.coli (APEC).

BACKGROUND: Receptor interacting serine/threonine kinase 2 (RIP2), ubiquitous in many tissue/cell types, is the key regulator of immune and inflammatory responses for many diseases, including avian pathogenic E. coli (APEC), which causes a wide variety of localized or systemic infections. However, the molecular mechanisms by which RIP2 drives its transcriptional program to affect immune and inflammatory response upon APEC infection remains poorly understood. RESULTS: In this study, RNA-seq and bioinformatics analyses were used to detect gene expression and new direct/indirect RIP2 targets in the treatments of wild type HD11 cells (WT), RIP2 knockdown cells (shRIP2), APEC stimulation cells (APEC), and RIP2 knockdown cells combined with APEC infection (shRIP2 + APEC). The results revealed that a total of 4691 and 2605 differentially expressed genes (DEGs) were screened in shRIP2 + APEC vs. APEC and shRIP2 vs. WT, respectively. Functional annotation analysis showed that apoptosis, MAPK, p53, Toll-like receptor, and Nod-like receptor signaling pathways were involved in APEC-induced RIP2 knockdown HD11 cells. By analyzing the enriched pathway and gene networks, we identified that several DEGs, including HSP90AB1, BID, and CASP9 were targeted by RIP2 upon APEC infection. CONCLUSION: As a whole, this study can not only provide data support for constructing gene networks of RIP2 knockdown with APEC challenge but also provide new ideas for improving the immune and inflammatory response.

Analysis of miRNA Expression Profiling of RIP2 Knockdown in Chicken HD11 Cells When Infected with Avian Pathogenic E. coli (APEC).

Colibacillosis is an acute and chronic avian disease caused by avian pathogenic E. coli (APEC). Previous studies have demonstrated that RIP2 plays a significant role in APEC infection. Moreover, increasing evidence indicates that microRNAs (miRNAs) are involved in host-pathogen interactions and the immune response. However, the role of miRNAs in the host against APEC infection remains unclear. Herein, we attempted to reveal new miRNAs potentially involved in the regulation of the immune and inflammatory response against APEC infection, with a particular focus on those possibly correlated with RIP2 expression, via miRNA-seq, RT-qPCR, Western blotting, dual-luciferase reporter assay, and CCK-8. The results showed that a total of 93 and 148 differentially expressed (DE) miRNAs were identified in the knockdown of RIP2 cells following APEC infection (shRIP2+APEC) vs. knockdown of RIP2 cells (shRIP2) and shRIP2 vs. wild-type cells (WT), respectively. Among those identified DE miRNAs, the biological function of gga-miR-455-5p was investigated. It was found that gga-miR-455-5p regulated by RIP2 was involved in the immune and inflammatory response against APEC infection via targeting of IRF2 to modulate the expression of type I interferons. Additionally, RIP2 could directly regulate the production of the type I interferons. Altogether, these findings highlighted the crucial role of miRNAs, especially gga-miR-455-5p, in host defense against APEC infection.

Transcriptional Regulation of RIP2 Gene by NFIB Is Associated with Cellular Immune and Inflammatory Response to APEC Infection.

Avian pathogenic E. coli (APEC) can cause localized or systemic infection, resulting in large economic losses per year, and impact health of humans. Previous studies showed that RIP2 (receptor interacting serine/threonine kinase 2) and its signaling pathway played an important role in immune response against APEC infection. In this study, chicken HD11 cells were used as an in vitro model to investigate the function of chicken RIP2 and the transcription factor binding to the RIP2 core promoter region via gene overexpression, RNA interference, RT-qPCR, Western blotting, dual luciferase reporter assay, CHIP-PCR, CCK-8, and flow cytometry assay following APEC stimulation. Results showed that APEC stimulation promoted RIP2 expression and cells apoptosis, and inhibited cells viability. Knockdown of RIP2 significantly improved cell viability and suppressed the apoptosis of APEC-stimulated cells. Transcription factor NFIB (Nuclear factor I B) and GATA1 (globin transcription factor 1) binding site was identified in the core promoter region of RIP2 from -2300 bp to -1839 bp. However, only NFIB was confirmed to be bound to the core promoter of RIP2. Overexpression of NFIB exacerbated cell injuries with significant reduction in cell viability and increased cell apoptosis and inflammatory cytokines levels, whereas opposite results were observed in NFIB inhibition treatment group. Moreover, RIP2 was up-regulated by NFIB overexpression, and RIP2 silence mitigated the effect of NFIB overexpression in cell apoptosis, inflammation, and activation of NFκB signaling pathways. This study demonstrated that NFIB overexpression accelerated APEC-induced apoptosis and inflammation via up-regulation of RIP2 mediated downstream pathways in chicken HD11 cells.

Identification of candidate genomic regions for chicken egg number traits based on genome-wide association study.

BACKGROUND: Since the domestication of chicken, various breeds have been developed for food production, entertainment, and so on. Compared to indigenous chicken breeds which generally do not show elite production performance, commercial breeds or lines are selected intensely for meat or egg production. In the present study, in order to understand the molecular mechanisms underlying the dramatic differences of egg number between commercial egg-type chickens and indigenous chickens, we performed a genome-wide association study (GWAS) in a mixed linear model. RESULTS: We obtained 148 single nucleotide polymorphisms (SNPs) associated with egg number traits (57 significantly, 91 suggestively). Among them, 4 SNPs overlapped with previously reported quantitative trait loci (QTL), including 2 for egg production and 2 for reproductive traits. Furthermore, we identified 32 candidate genes based on the function of the screened genes. These genes were found to be mainly involved in regulating hormones, playing a role in the formation, growth, and development of follicles, and in the development of the reproductive system. Some genes such as NELL2 (neural EGFL like 2), KITLG (KIT ligand), GHRHR (Growth hormone releasing hormone receptor), NCOA1 (Nuclear receptor coactivator 1), ITPR1 (inositol 1, 4, 5-trisphosphate receptor type 1), GAMT (guanidinoacetate N-methyltransferase), and CAMK4 (calcium/calmodulin-dependent protein kinase IV) deserve our attention and further study since they have been reported to be closely related to egg production, egg number and reproductive traits. In addition, the most significant genomic region obtained in this study was located at 48.61-48.84 Mb on GGA5. In this region, we have repeatedly identified four genes, in which YY1 (YY1 transcription factor) and WDR25 (WD repeat domain 25) have been shown to be related to oocytes and reproductive tissues, respectively, which implies that this region may be a candidate region underlying egg number traits. CONCLUSION: Our study utilized the genomic information from various chicken breeds or populations differed in the average annual egg number to understand the molecular genetic mechanisms involved in egg number traits. We identified a series of SNPs, candidate genes, or genomic regions that associated with egg number, which could help us in developing the egg production trait in chickens.

Alternative polyadenylation coordinates embryonic development, sexual dimorphism and longitudinal growth in Xenopus tropicalis.

RNA alternative polyadenylation contributes to the complexity of information transfer from genome to phenome, thus amplifying gene function. Here, we report the first X. tropicalis resource with 127,914 alternative polyadenylation (APA) sites derived from embryos and adults. Overall, APA networks play central roles in coordinating the maternal-zygotic transition (MZT) in embryos, sexual dimorphism in adults and longitudinal growth from embryos to adults. APA sites coordinate reprogramming in embryos before the MZT, but developmental events after the MZT due to zygotic genome activation. The APA transcriptomes of young adults are more variable than growing adults and male frog APA transcriptomes are more divergent than females. The APA profiles of young females were similar to embryos before the MZT. Enriched pathways in developing embryos were distinct across the MZT and noticeably segregated from adults. Briefly, our results suggest that the minimal functional units in genomes are alternative transcripts as opposed to genes.

Evolution of cis- and trans-regulatory divergence in the chicken genome between two contrasting breeds analyzed using three tissue types at one-day-old.

BACKGROUND: Gene expression variation is a key underlying factor influencing phenotypic variation, and can occur via cis- or trans-regulation. To understand the role of cis- and trans-regulatory variation on population divergence in chicken, we developed reciprocal crosses of two chicken breeds, White Leghorn and Cornish Game, which exhibit major differences in body size and reproductive traits, and used them to determine the degree of cis versus trans variation in the brain, liver, and muscle tissue of male and female 1-day-old specimens. RESULTS: We provided an overview of how transcriptomes are regulated in hybrid progenies of two contrasting breeds based on allele specific expression analysis. Compared with cis-regulatory divergence, trans-acting genes were more extensive in the chicken genome. In addition, considerable compensatory cis- and trans-regulatory changes exist in the chicken genome. Most importantly, stronger purifying selection was observed on genes regulated by trans-variations than in genes regulated by the cis elements. CONCLUSIONS: We present a pipeline to explore allele-specific expression in hybrid progenies of inbred lines without a specific reference genome. Our research is the first study to describe the regulatory divergence between two contrasting breeds. The results suggest that artificial selection associated with domestication in chicken could have acted more on trans-regulatory divergence than on cis-regulatory divergence.

A genome-wide association study explores the genetic determinism of host resistance to Salmonella pullorum infection in chickens.

BACKGROUND: Salmonella infection is a serious concern in poultry farming because of its impact on both economic loss and human health. Chicks aged 20 days or less are extremely vulnerable to Salmonella pullorum (SP), which causes high mortality. Furthermore, an outbreak of SP infection can result in a considerable number of carriers that become potential transmitters, thus, threatening fellow chickens and offspring. In this study, we conducted a genome-wide association study (GWAS) to detect potential genomic loci and candidate genes associated with two disease-related traits: death and carrier state. METHODS: In total, 818 birds were phenotyped for death and carrier state traits through a SP challenge experiment, and genotyped by using a 600 K high-density single nucleotide polymorphism (SNP) array. A GWAS using a single-marker linear mixed model was performed with the GEMMA software. RNA-sequencing on spleen samples was carried out for further identification of candidate genes. RESULTS: We detected a region that was located between 33.48 and 34.03 Mb on chicken chromosome 4 and was significantly associated with death, with the most significant SNP (rs314483802) accounting for 11.73% of the phenotypic variation. Two candidate genes, FBXW7 and LRBA, were identified as the most promising genes involved in resistance to SP. The expression levels of FBXW7 and LRBA were significantly downregulated after SP infection, which suggests that they may have a role in controlling SP infections. Two other significant loci and related genes (TRAF3 and gga-mir-489) were associated with carrier state, which indicates a different polygenic determinism compared with that of death. In addition, genomic inbreeding coefficients showed no correlation with resistance to SP within each breed in our study. CONCLUSIONS: The results of this GWAS with a carefully organized Salmonella challenge experiment represent an important milestone in understanding the genetics of infectious disease resistance, offer a theoretical basis for breeding SP-resistant chicken lines using marker-assisted selection, and provide new information for salmonellosis research in humans and other animals.

Genome-wide association study revealed genomic regions related to white/red earlobe color trait in the Rhode Island Red chickens.

BACKGROUND: Earlobe color is a naturally and artificially selected trait in chicken. As a head furnishing trait, it has been selected as a breed characteristic. Research has demonstrated that white/red earlobe color was related to at least three loci and sex-linked. However, there has been little work to date to identify the specific genomic regions and genes response to earlobe color in Rhode Island Red chickens. Currently, it is possible to identify the genomic regions responsible for white/red earlobe in Rhode Island Red chicken to eliminate this gap in knowledge by using genome-wide association (GWA) analysis. RESULTS: In the present study, genome-wide association (GWA) analysis was conducted to explore the candidate genomic regions response to chicken earlobe color phenotype. Hens with red dominant and white dominant earlobe was used for case-control analysis by Illumina 600 K SNP arrays. The GWA results showed that 2.38 Mb genomic region (50.13 to 52.51 Mb) with 282 SNPs on chromosome Z were significantly correlated to earlobe color, including sixteen known genes and seven anonymous genes. The sixteen genes were PAM, SLCO4C1, ST8SIA4, FAM174A, CHD1, RGMB, RIOK2, LIX1, LNPEP, SHB, RNF38, TRIM14, NANS, CLTA, GNE, and CPLX1. CONCLUSIONS: The study has revealed the white/red earlobe trait is polygenic and sex-linked in Rhode Island Red chickens. In the genome significant ~2.38 Mb region, twenty-three genes were found and some of them could play critical roles in the formation of white/red earlobe color, especially gene SLCO4C1. Taken together, the candidate genes findings herein can help elucidate the genomic architecture of response to white/red earlobe and provide a new insight on mechanisms underlying earlobe color in Rhode Island Red chickens and other breeds.

Chicken gga-miR-130a targets HOXA3 and MDFIC and inhibits Marek's disease lymphoma cell proliferation and migration.

Marek's disease (MD) is an infectious disease of chickens caused by MD virus (MDV), which is a herpesvirus that initiates tumor formation. Studies have indicated that microRNAs (miRNAs) are linked with the development of cancers or tumors. Previously, gga-miR-130a was discovered downregulated in MDV-infected tissues. Here, we aimed to explore the further function of gga-miR-130a in MD. The expression of gga-miR-130a in MDV-infected and uninfected spleens was detected by quantitative real-time PCR (qRT-PCR). Subsequently, proliferation and migration assays of MDV-transformed lymphoid cells (MSB1) were carried out by transfecting gga-miR-130a. The target genes of gga-miR-130a were predicted using TargetScan and miRDB and clustered through Gene Ontology analysis. The target genes were validated by western blot, qRT-PCR, and a dual luciferase reporter assay. Our results show that the expression of gga-miR-130a was reduced in MDV-infected spleens. Gga-miR-130a showed an inhibitory effect on MSB1 cell proliferation and migration. Two target genes, homeobox A3 (HOXA3) and MyoD family inhibitor domain containing (MDFIC), were predicted and clustered to cell proliferation. Results indicate that gga-miR-130a regulates HOXA3 and MDFIC at the protein level but not at the mRNA level. Moreover, the gga-miR-130a binding sites of two target genes have been confirmed. We conclude that gga-miR-130a can arrest MSB1 cell proliferation and migration, and target HOXA3 and MDFIC, which are both involved in the regulation of cell proliferation. Collectively, gga-miR-130a plays a critical role in the tumorigenesis associated with chicken MD.

An EAV-HP insertion in 5' Flanking region of SLCO1B3 causes blue eggshell in the chicken.

The genetic determination of eggshell coloration has not been determined in birds. Here we report that the blue eggshell is caused by an EAV-HP insertion that promotes the expression of SLCO1B3 gene in the uterus (shell gland) of the oviduct in chicken. In this study, the genetic map location of the blue eggshell gene was refined by linkage analysis in an F(2) chicken population, and four candidate genes within the refined interval were subsequently tested for their expression levels in the shell gland of the uterus from blue-shelled and non-blue-shelled hens. SLCO1B3 gene was found to be the only one expressed in the uterus of blue-shelled hens but not in that of non-blue-shelled hens. Results from a pyrosequencing analysis showed that only the allele of SLCO1B3 from blue-shelled chickens was expressed in the uterus of heterozygous hens (O*LC/O*N). SLCO1B3 gene belongs to the organic anion transporting polypeptide (OATP) family; and the OATPs, functioning as membrane transporters, have been reported for the transportation of amphipathic organic compounds, including bile salt in mammals. We subsequently resequenced the whole genomic region of SLCO1B3 and discovered an EAV-HP insertion in the 5' flanking region of SLCO1B3. The EAV-HP insertion was found closely associated with blue eggshell phenotype following complete Mendelian segregation. In situ hybridization also demonstrated that the blue eggshell is associated with ectopic expression of SLCO1B3 in shell glands of uterus. Our finding strongly suggests that the EAV-HP insertion is the causative mutation for the blue eggshell phenotype. The insertion was also found in another Chinese blue-shelled breed and an American blue-shelled breed. In addition, we found that the insertion site in the blue-shelled chickens from Araucana is different from that in Chinese breeds, which implied independent integration events in the blue-shelled chickens from the two continents, providing a parallel evolutionary example at the molecular level.

Genome-wide identification of copy number variations between two chicken lines that differ in genetic resistance to Marek's disease.

BACKGROUND: Copy number variation (CNV) is a major source of genome polymorphism that directly contributes to phenotypic variation such as resistance to infectious diseases. Lines 63 and 72 are two highly inbred experimental chicken lines that differ greatly in susceptibility to Marek's disease (MD), and have been used extensively in efforts to identify the genetic and molecular basis for genetic resistance to MD. Using next generation sequencing, we present a genome-wide assessment of CNVs that are potentially associated with genetic resistance to MD. METHODS: Three chickens randomly selected from each line were sequenced to an average depth of 20×. Two popular software, CNVnator and Pindel, were used to call genomic CNVs separately. The results were combined to obtain a union set of genomic CNVs in the two chicken lines. RESULTS: A total of 5,680 CNV regions (CNVRs) were identified after merging the two datasets, of which 1,546 and 1,866 were specific to the MD resistant or susceptible line, respectively. Over half of the line-specific CNVRs were shared by 2 or more chickens, reflecting the reduced diversity in both inbred lines. The CNVRs fixed in the susceptible lines were significantly enriched in genes involved in MAPK signaling pathway. We also found 67 CNVRs overlapping with 62 genes previously shown to be strong candidates of the underlying genes responsible for the susceptibility to MD. CONCLUSIONS: Our findings provide new insights into the genetic architecture of the two chicken lines and additional evidence that MAPK signaling pathway may play an important role in host response to MD virus infection. The rich source of line-specific CNVs is valuable for future disease-related association studies in the two chicken lines.

Genome-wide association study revealed a promising region and candidate genes for eggshell quality in an F2 resource population.

BACKGROUND: Eggshell is subject to quality loss with aging process of laying hens, and damaged eggshells result in economic losses of eggs. However, the genetic architecture underlying the dynamic eggshell quality remains elusive. Here, we measured eggshell quality traits, including eggshell weight (ESW), eggshell thickness (EST) and eggshell strength (ESS) at 11 time points from onset of laying to 72 weeks of age and conducted comprehensive genome-wide association studies (GWAS) in 1534 F2 hens derived from reciprocal crosses between White Leghorn (WL) and Dongxiang chickens (DX). RESULTS: ESWs at all ages exhibited moderate SNP-based heritability estimates (0.30 ~ 0.46), while the estimates for EST (0.21 ~ 0.31) and ESS (0.20 ~ 0.27) were relatively low. Eleven independent univariate genome-wide screens for each trait totally identified 1059, 1026 and 1356 significant associations with ESW, EST and ESS, respectively. Most significant loci were in a region spanning from 57.3 to 71.4 Mb of chromosome 1 (GGA1), which together account for 8.4 ~ 16.5% of the phenotypic variance for ESW from 32 to 72 weeks of age, 4.1 ~ 6.9% and 2.95 ~ 16.1% for EST and ESS from 40 to 72 weeks of age. According to linkage disequilibrium (LD) and conditional analysis, the significant SNPs in this region were in extremely strong linkage disequilibrium status. Ultimately, two missense SNPs in GGA1 and one in GGA4 were considered as promising loci on three independent genes including ITPR2, PIK3C2G, and NCAPG. The homozygotes of advantageously effective alleles on PIK3C2G and ITPR2 possessed the best eggshell quality and could partly counteract the negative effect of aging process. NCAPG had certain effect on eggshell quality for young hens. CONCLUSIONS: Identification of the promising region as well as potential candidate genes will greatly advance our understanding of the genetic basis underlying dynamic eggshell quality and has the practical significance in breeding program for the improvement of eggshell quality, especially at the later part of laying cycle.

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

BACKGROUND: As a major economic trait in chickens, egg weight (EW) receives widespread interests in breeding, production and consumption. However, limited information is available for underlying genetic architecture of longitudinal trend in EW. Herein, we measured EWs at nine time points from onset of laying to 60 week of age, and conducted comprehensive genome-wide association studies (GWAS) in 1,534 F2 hens derived from reciprocal crosses between White Leghorn and Dongxiang chickens. RESULTS: Egg weights at all ages except the first egg weight (FEW) exhibited high SNP-based heritability estimates (0.47~0.60). Strong pair-wise genetic correlations (0.77~1.00) were found among all EWs. Nine separate univariate genome-wide screens suggested 73 signals showing significant associations with longitudinal EWs. After multivariate and conditional analyses, four variants on three chromosomes remained independent contributions. The minor alleles at two loci exerted consistent and positive substitution effects on EWs, and other two were negative. The four loci together accounted for 3.84 % of the phenotypic variance for FEW and 7.29~11.06 % for EWs from 32 to 60 week of age. We obtained five candidate genes, of which NCAPG harbors a non-synonymous SNP (rs14491030) causing a valine-to-alanine amino-acid substitution. Genome partitioning analysis indicated a strong linear correlation between the variance explained by each chromosome and its length, which provided evidence that EW follows a highly polygenic nature of inheritance. CONCLUSIONS: Identification of significant genetic causes that together implicate EWs at different ages will greatly advance our understanding of the genetic basis behind longitudinal EWs, and would be helpful to illuminate the future breeding direction on how to select desired egg size.

Genome-wide association studies for feed intake and efficiency in two laying periods of chickens.

BACKGROUND: Feed contributes to over 60 % of the total production costs in the poultry industry. Increasing feed costs prompt geneticists to include feed intake and efficiency as selection goals in breeding programs. In the present study, we used an F2 chicken population in a genome-wide association study (GWAS) to detect potential genetic variants and candidate genes associated with daily feed intake (FI) and feed efficiency, including residual feed intake (RFI) and feed conversion ratio (FCR). METHODS: A total of 1534 F2 hens from a White Leghorn and Dongxiang reciprocal cross were phenotyped for feed intake and efficiency between 37 and 40 weeks (FI1, RFI1, and FCR1) and between 57 and 60 weeks (FI2, RFI2, and FCR2), and genotyped using the chicken 600 K single nucleotide polymorphism (SNP) genotyping array. Univariate, bivariate, and conditional genome-wide association studies (GWAS) were performed with GEMMA, a genome-wide efficient mixed model association algorithm. The statistical significance threshold for association was inferred by the simpleM method. RESULTS: We identified eight genomic regions that each contained at least one genetic variant that showed a significant association with FI. Genomic regions on Gallus gallus (GGA) chromosome 4 coincided with known quantitative trait loci (QTL) that affect feed intake of layers. Of particular interest, eight SNPs on GGA1 in the region between 169.23 and 171.55 Mb were consistently associated with FI in both univariate and bivariate GWAS, which explained 3.72 and 2.57 % of the phenotypic variance of FI1 and FI2, respectively. The CAB39L gene can be considered as a promising candidate for FI1. For RFI, a haplotype block on GGA27 harbored a significant SNP associated with RFI2. The major allele of rs315135692 was favorable for a lower RFI, with a phenotypic difference of 3.35 g/day between opposite homozygous genotypes. Strong signals on GGA1 were detected in the bivariate GWAS for FCR. CONCLUSIONS: The results demonstrated the polygenic nature of feed intake. GWAS identified novel variants and confirmed a QTL that was previously reported for feed intake in chickens. Genetic variants associated with feed efficiency may be used in genomic breeding programs to select more efficient layers.