A comprehensive characterization of regulatory elements in the chicken genome across tissues will have substantial impacts on both fundamental and applied research. Here, we systematically identified and characterized regulatory elements in the chicken genome by integrating 377 genome-wide sequencing datasets from 23 adult tissues. In total, we annotated 1.57 million regulatory elements, representing 15 distinct chromatin states, and predicted about 1.2 million enhancer-gene pairs and 7662 super-enhancers. This functional annotation of the chicken genome should have wide utility on identifying regulatory elements accounting for gene regulation underlying domestication, selection, and complex trait regulation, which we explored. In short, this comprehensive atlas of regulatory elements provides the scientific community with a valuable resource for chicken genetics and genomics.

Objective: Follicle selection is an important process in chicken egg laying. Among several small yellow (SY) follicles, the one exhibiting the highest expression of follicle stimulation hormone receptor (FSHR) will be selected to become a hierarchal follicle. The role of lncRNA, miRNA and other non-coding RNA in chicken follicle selection is unclear. Methods: In this study, the whole transcriptome sequencing of SY follicles with different expression levels of FSHR in Jining Bairi hens was performed, and the expression of 30 randomly selected mRNAs, lncRNAs and miRNAs was validated by quantitative real-time polymerase chain reaction. Preliminary studies and bioinformatics analysis were performed on the selected mRNA, lncRNA, miRNA and their target genes. The effect of identified gene was examined in the granulosa cells of chicken follicles. Results: Integrated transcriptomic analysis on chicken SY follicles differing in FSHR expression revealed 467 differentially expressed mRNA genes, 134 differentially expressed lncRNA genes and 34 differentially expressed miRNA genes, and sosondowah ankyrin repeat domain family member A (SOWAHA) was the common target gene of three miRNAs and one lncRNA. SOWAHA was mainly expressed in small white (SW) and SY follicles and was affected by follicle stimulation hormone (FSH) treatment in the granulosa cells. Knockdown of SOWAHA inhibited the expression of Wnt family member 4 (Wnt4) and steroidogenic acute regulatory protein (StAR) in the granulosa cells of prehierarchal follicles, while stimulated Wnt4 in hierarchal follicles. Overexpression of SOWAHA increased the expression of Wnt4 in the granulosa cells of prehierarchal follicles, decreased that of StAR and cytochrome P450 family 11 subfamily A member 1 in the granulosa cells of hierarchal follicles and inhibited the proliferation of granulosa cells. Conclusion: Integrated analysis of chicken SY follicle transcriptomes identified SOWAHA as a network gene that is affected by FSH in granulosa cells of ovarian follicles. SOWAHA affected the expression of genes involved in chicken follicle selection and inhibited the proliferation of granulosa cells, suggesting an inhibitory role in chicken follicle selection.

Chicken is a valuable model for understanding fundamental biology, vertebrate evolution and diseases, as well as a major source of nutrient-dense and lean-protein-enriched food globally. Although it is the first non-mammalian amniote genome to be sequenced, the chicken genome still lacks a systematic characterization of functional impacts of genetic variants. Here, through integrating 7,015 RNA-Seq and 2,869 whole-genome sequence data, the Chicken Genotype- Tissue Expression (ChickenGTEx) project presents the pilot reference of regulatory variants in 28 chicken tissue transcriptomes, including millions of regulatory effects on primary expression (including protein-coding genes, lncRNA and exon) and post-transcriptional modifications (alternative splicing and 3’ untranslated region alternative polyadenylation). We explored the tissue-sharing and context-specificity of these regulatory variants, their underlying molecular mechanisms of action, and their utility in interpreting adaptation and genome-wide associations of 108 chicken complex traits. Finally, we illustrated shared and lineage-specific features of gene regulation between chickens and mammals, and demonstrated how the ChickenGTEx resource can further assist with translating genetic findings across species. One-Sentence Summary: The ChickenGTEx provides a multi-tissue reference of regulatory variants for chicken genetics and genomics, functional genomics, precision breeding, veterinary medicine, vertebrate evolution and even human biomedicine.

The signaling pathway of the wingless-type mouse mammary tumor virus integration site (Wnt) plays an important role in ovarian and follicular development. In our previous study, WNT4 was shown to be involved in the selection and development of chicken follicles by upregulating the expression of follicle-stimulating hormone receptors (FSHR), stimulating the proliferation of follicular granulosa cells, and increasing the secretion of steroidal hormones. FSH also stimulates the expression of WNT4. To further explore the molecular mechanism by which FSH upregulates WNT4 and characterize the cis-elements regulating WNT4 transcription, in this study, we determined the critical regulatory regions affecting chicken WNT4 transcription. We then identified a single-nucleotide polymorphism (SNP) in this region, and finally analyzed the associations of the SNP with chicken production traits. The results showed that the 5′ regulatory region from −3354 to −2689 of WNT4 had the strongest activity and greatest response to FSH stimulation, and we identified one SNP site in this segment, −3015 (G > C), as affecting the binding of NFAT5 (nuclear factor of activated T cells 5) and respones to FSH stimulation. When G was replaced with C at this site, it eliminated the NFAT5 binding. The mRNA level of WNT4 in small yellow follicles of chickens with genotype GG was significantly higher than that of the other two genotypes. Moreover, this locus was found to be significantly associated with comb length in hens. Individuals with the genotype CC had longer combs. Collectively, these data suggested that SNP−3015 (G > C) is involved in the regulation of WNT4 gene expression by responding FSH and affecting the binding of NFAT5 and that it is associated with chicken comb length. The current results provide a reference for further revealing the response mechanism between WNT and FSH.

Follicle selection is an important step in the laying process of chicken, which is closely related to the laying performance and fecundity of hens. Follicle selection mainly depends on the regulation of follicle-stimulating hormone (FSH) secreted by pituitary gland and the expression of follicle stimulation hormone receptor. To uncover the role of FSH in chicken follicle selection, in this study, we analyzed the changes in the mRNA transcriptome profiles of FSH-treated chicken granulosa cells from prehierarchical follicles by long-read sequencing Oxford Nanopore Technologies (ONT) approach. Among the 10,764 genes detected, 31 differentially expressed (DE) transcripts of 28 DE genes were significantly upregulated by FSH treatment. These DE transcripts (DETs) were mainly related to the steroid biosynthetic process by GO analysis and enriched in pathways of ovarian steroidogenesis and aldosterone synthesis and secretion by KEGG analysis. Among these genes, the mRNA and protein expression of TNF receptor associated factor 7 (TRAF7) was upregulated after FSH treatment. Further study revealed that TRAF7 stimulated the mRNA expression of steroidogenic enzymes steroidogenic acute regulatory protein (StAR) and cytochrome P450 family 11 subfamily A member 1 (CYP11A1) genes and the proliferation of granulosa cells. This is the first study to investigate differences in chicken prehierarchical follicular granulosa cells before and after FSH treatment by using ONT transcriptome sequencing, which provides a reference for a more comprehensive understanding of the molecular mechanism of follicle selection in chicken.