Transcriptome profiling analysis of underlying regulation of growing follicle development in the chicken

Characterisation of runs of homozygosity and inbreeding coefficients in the red-brown Korean native chickens

OBJECTIVE

The analysis of runs of homozygosity (ROH) has been applied to assess the level of inbreeding and identify selection signatures in various livestock species. The objectives of this study were to characterize the ROH pattern, estimate the rate of inbreeding, and identify signatures of selection in the red-brown Korean native chickens.

METHODS

The Illumina 60K single nucleotide polymorphism chip data of 651 chickens was used in the analysis. Runs of homozygosity were analysed using the PLINK v1.9 software. Inbreeding coefficients were estimated using the GCTA software and their correlations were examined. Genomic regions with high levels of ROH were explored to identify selection signatures.

RESULTS

A total of 32,176 ROH segments were detected in this study. The majority of the ROH segments were shorter than 4 Mb. The average ROH inbreeding coefficients (FROH) varied with the length of ROH segments. The means of inbreeding coefficients calculated from different methods were also variable. The correlations between different inbreeding coefficients were positive and highly variable (r = 0.18-1). Five ROH islands harbouring important quantitative trait loci were identified.

CONCLUSION

This study assessed the level of inbreeding and patterns of homozygosity in Red-brown native Korean chickens. The results of this study suggest that the level of recent inbreeding is low which indicates substantial progress in the conservation of red-brown Korean native chickens. Additionally, Candidate genomic regions associated with important production traits were detected in homozygous regions.

Transcriptome profiling reveals SLC5A5 regulates chicken ovarian follicle granulosa cell proliferation, apoptosis, and steroid hormone synthesis

The egg-laying performance of hens holds significant economic importance within the poultry industry. Broody inheritance of the parent stock of chickens can result in poor options for the improvement of egg production, and is a phenomenon influenced by multiple genetic factors. However, few studies have been conducted to delineate the molecular mechanism of ovarian regression in brooding chickens. Here, we explored the pivotal genes responsible for the regulation of ovarian follicles in laying hens, using RNA-sequencing analysis on the small ovarian follicles from broody and laying chickens. Sequencing data analysis revealed the differential expression of 200 genes, with a predominant enrichment in biological processes related to cell activation and metabolism. Among these genes, we focused on solute carrier family 5 member 5 (SLC5A5), which exhibited markedly higher RNA expression levels in follicles from laying compared with broody chickens. Subsequent cellular function studies with knockdown of SLC5A5 in chicken ovarian follicle granulosa cells (GCs) led to the down-regulation of genes associated with cell proliferation and steroid hormone synthesis, and concurrent promotion of gene expression linked to apoptosis. These findings indicated that SLC5A5 deficiency led to the inhibition of proliferation, steroid hormone synthesis and secretion, and promotion of apoptosis in chicken GCs. Our study demonstrated a pivotal role for SLC5A5 in the development and function of chicken GCs, shedding light on its potential significance in the broader context of chicken ovarian follicle development, and providing a prospective target to improve the egg-laying performance of chickens via molecular marker-assisted breeding technology.

Identification of genes related to growth and amino acid metabolism from the transcriptome profile of the liver of growing laying hens

The growing period is a critical period for the growth and development of hens and affects their production performance during the laying period. During the early stage of growing, bone and muscle growth is rapid, making it necessary to provide sufficient amino acids (AA) to support the growth and development of laying hens. In this experiment, RNA-Sequencing (RNA-Seq) was applied to compare the liver tissues from 6- to 12-wk-old growing laying hens to identify candidate genes related to growth and AA transport and metabolism. In the liver tissues, 596 differentially expressed genes (DEG) were identified, of which 424 genes were up-regulated and 172 were down-regulated. Through enrichment analysis and DEGs analysis, some DEGs and pathways related to AA transport and metabolism were identified. Additionally, there were significantly increased activities in the liver of glutamate dehydrogenase (GDH), glutamic oxaloacetic transaminase (GOT), and glutamate pyruvate transaminase (GPT). Meanwhile, the level of serum insulin-like growth factor binding protein-5 (IGFBP-5) significantly elevated, and insulin-like growth factor-1 (IGF-1) levels significantly reduced at 12 wk compared to 6 wk. The AA contents in the breast muscle were not significantly altered, while the levels of the free AA in the serum underwent significant changes. This study discovered that the transport and metabolism of AAs in growing laying hens at different ages changed, which influenced the growth and development of growing laying hens. This contributes to future research on the mechanisms of growth and AA metabolism during the growing period of laying hens.

Transcriptome analysis of ovarian tissues highlights genes controlling energy homeostasis and oxidative stress as potential drivers of heterosis for egg number and clutch size in crossbred laying hens

Heterosis is the major benefit of crossbreeding and has been exploited in laying hens breeding for a long time. This genetic phenomenon has been linked to various modes of nonadditive gene action. However, the molecular mechanism of heterosis for egg production in laying hens has not been fully elucidated. To fill this research gap, we sequenced mRNAs and lncRNAs of the ovary stroma containing prehierarchical follicles in White Leghorn, Rhode Island Red chickens as well as their reciprocal crossbreds that demonstrated heterosis for egg number and clutch size. We further delineated the modes of mRNAs and lncRNAs expression to identify their potential functions in the observed heterosis. Results showed that dominance was the principal mode of nonadditive expression exhibited by mRNAs and lncRNAs in the prehierarchical follicles of crossbred hens. Specifically, low-parent dominance was the main mode of mRNA expression, while high-parent dominance was the predominant mode of lncRNA expression. Important pathways enriched by genes that showed higher expression in crossbreds compared to either one or both parental lines were cell adhesion molecules, tyrosine and purine metabolism. In contrast, ECM-receptor interaction, focal adhesion, PPAR signaling, and ferroptosis were enriched in genes with lower expression in the crossbred. Protein network interaction identified nonadditively expressed genes including apolipoprotein B (APOB), transferrin, acyl-CoA synthetase medium-chain family member (APOBEC) 3, APOBEC1 complementation factor, and cathepsin S as hub genes. Among these potential hub genes, APOB was the only gene with underdominance expression common to the 2 reciprocal crossbred lines, and has been linked to oxidative stress. LncRNAs with nonadditive expression in the crossbred hens targeted natriuretic peptide receptor 1, epidermal differentiation protein beta, spermatogenesis-associated gene 22, sperm-associated antigen 16, melanocortin 2 receptor, dolichol kinase, glycine amiinotransferase, and prolactin releasing hormone receptor. In conclusion, genes with nonadditive expression in the crossbred may play crucial roles in follicle growth and atresia by improving follicle competence and increasing oxidative stress, respectively. These 2 phenomena could underpin heterosis for egg production in crossbred laying hens.

Dynamic mRNA expression during chicken ovarian follicle development

Ovarian follicle development is a complex and well-orchestrated biological process of great economic significance for poultry production. Specifically, understanding the molecular mechanisms underlying follicular development is essential for high-efficiency follicular development can benefit the entire industry. In addition, domestic egg-laying hens often spontaneously develop ovarian cancer, providing an opportunity to study the genetic, biochemical, and environmental risk factors associated with the development of this cancer. Here, we provide high-quality RNA sequencing data for chicken follicular granulosa cells across 10 developmental stages, which resulted in a total of 204.57 Gb of clean sequencing data (6.82 Gb on average per sample). We also performed gene expression, time-series, and functional enrichment analyses across the 10 developmental stages. Our study revealed that SWF (small while follicle), F1 (F1 hierarchical follicles), and POFs (postovulatory follicles) best represent the transcriptional changes associated with the prehierarchical, preovulatory, and postovulatory stages, respectively. We found that the preovulatory stage F1 showed the greatest divergence in gene expression from the POF stage. Our research lays a foundation for further elucidation of egg-laying performance of chicken and human ovarian disease.

Transcriptome profiling analysis of uterus during chicken laying periods

The avian eggshell is formed in the uterus. Changes in uterine function may have a significant effect on eggshell quality. To identify the vital genes impacting uterine functional maintenance in the chicken, uteri in three different periods (22W, 31W, 51W) were selected for RNA sequencing and bioinformatics analysis. In our study, 520, 706 and 736 differentially expressed genes (DEGs) were respectively detected in the W31 vs W22 group, W51 vs W31 group and W51 vs W22 group. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated DEGs were enriched in the extracellular matrix, extracellular region part, extracellular region, extracellular matrix structural constituent, ECM receptor interaction, collagen-containing extracellular matrix and collagen trimer in the uterus (P < 0.05). Protein-protein interaction analysis revealed that FN1, LOX, THBS2, COL1A1, COL1A2, COL5A1, COL5A2, POSTN, MMP13, VANGL2, RAD54B, SPP1, SDC1, BTC, ANGPTL3 might be key candidate genes for uterine functional maintenance in chicken. This study discovered dominant genes and pathways which enhanced our knowledge of chicken uterine functional maintenance.

Influence of microbiota inoculum as a substitute for antibiotic growth promoter during the initial laying phase on productivity performance, egg quality, and the morphology of reproductive organs in laying hens

Background and Aim

Antibiotics that increase growth have long been employed as a component of chicken growth. Long-term, unchecked usage may lead to microbial imbalance, resistance, and immune system suppression. Probiotics are a suitable and secure feed additive that may be provided as a solution. The objective of this research was to ascertain the effects of dietary multistrain probiotics (Lactobacillus acidophilus, Bifidobacterium spp., and Lactobacillus plantarum) on the morphology (length and weight) of reproductive organs and productivity performance of laying hens during the early stage of laying.

Materials and Methods

One hundred ISA Brown commercial layer chicks of the same body weight (BW) that were 5 days old were divided into five treatments, each with four replicates and four chicks in each duplicate. There were five different dietary interventions: (T1) 100% base feed; (T2) base feed with 2.5 g of antibiotic growth promoter/kg feed; (T3) base feed plus probiotics; (T4) base feed at 1 mL/kg with probiotics; and (T5) base feed with probiotics, 3 mL/kg feed, 5 mL/kg of feed. The parameters observed were performance, internal and exterior egg quality, and the morphology (length and weight) of laying hens' reproductive organs.

Results

Probiotic supplementation (L. acidophilus, Bifidobacterium, and L. plantarum) significantly affected the BW, feed intake, egg weight, yolk index, albumin index, Haugh unit, egg height, egg width, and morphology (length and weight) of laying hens' reproductive organs compared to the control group (basic feed). In addition, there was no discernible difference between treatment groups in theeggshell weight and thickness variables across all treatment groups.

Conclusion

When laying hens were between 17 and 21 weeks old, during the early laying period, microbiota inoculum supplements (L. acidophilus, Bifidobacterium, and L. plantarum) increased growth, the quality of the internal and external layers' eggs, and the morphology of the laying hens' reproductive organs.

1α,25-Dihydroxyvitamin D3 Improves Follicular Development and Steroid Hormone Biosynthesis by Regulating Vitamin D Receptor in the Layers Model

1α,25-Dihydroxyvitamin D₃ (VitD₃) is the active form of vitamin D, and it regulates gene expression and protein synthesis in mammalian follicle development. However, the function of VitD₃ in the follicular development of layers remains unclear. This study investigated, through in vivo and in vitro experiments, the effects of VitD₃ on follicle development and steroid hormone biosynthesis in young layers. In vivo, ninety 18-week-old Hy-Line Brown laying hens were randomly divided into three groups for different treatments of VitD₃ (0, 10, and 100 μg/kg). VitD₃ supplementation promoted follicle development, increasing the number of small yellow follicles (SYFs) and large yellow follicles (LYFs) and the thickness of the granulosa layer (GL) of SYFs. Transcriptome analysis revealed that VitD₃ supplementation altered gene expression in the ovarian steroidogenesis, cholesterol metabolism, and glycerolipid metabolism signaling pathways. Steroid hormone-targeted metabolomics profiling identified 20 steroid hormones altered by VitD₃ treatment, with 5 being significantly different among the groups. In vitro, it was found that VitD₃ increased cell proliferation, promoted cell-cycle progression, regulated the expression of cell-cycle-related genes, and inhibited the apoptosis of granulosa cells from pre-hierarchical follicles (phGCs) and theca cells from prehierarchical follicles (phTCs). In addition, the steroid hormone biosynthesis-related genes, estradiol (E2) and progesterone (P4) concentrations, and vitamin D receptor (VDR) expression level was significantly altered by VitD₃. Our findings identified that VitD₃ altered the gene expression related to steroid metabolism and the production of testosterone, estradiol, and progesterone in the pre-hierarchical follicles (PHFs), resulting in positive effects on poultry follicular development.

miR-140-3p promotes follicle granulosa cell proliferation and steroid hormone synthesis via targeting AMH in chickens

Granulosa cells (GCs) are the ovary's most critical cells since they undergo cell differentiation and hormone synthesis changes closely associated with follicle development. While micro RNA 140-3p (miRNA-140-3p) has an apparent cell signaling role, particularly in cell proliferation, its biological role in chicken ovarian follicle growth and development remains elusive. This study explored miR-140-3p's effects on chicken GC proliferation and steroid hormone synthesis. MiR-140-3p dramatically increased GC proliferation, prevented apoptosis, increased progesterone synthesis, and enhanced gene expression related to steroid hormone synthesis. In addition, the anti-Müllerian hormone (AMH) gene was identified as a direct miR-140-3p target. MiR-140-3p abundance correlated negatively with AMH mRNA and protein levels in GCs. Our findings show that miR-140-3p influences chicken GC proliferation and steroid hormone synthesis by suppressing AMH expression.

Untargeted Metabolomics Revealed Potential Biomarkers of Small Yellow Follicles of Chickens during Sexual Maturation

Sexual maturation provides economically important traits in poultry production. Research on the initiation mechanism of sexual maturity is of great significance for breeding high-yield laying hens. However, the underlying mechanisms are not fully clear. Here, one hundred and fifty Chahua No. 2 laying hens (the CH2 group, which has precocious puberty) and one hundred and fifty Wu Liang Shan black-bone laying hens (the WLS group, a late-maturing chicken breed) with similar weights and ages were randomly selected. ELISA was used to determine the secretion levels of luteinizing hormone (LH), estradiol (E2), and progesterone (P4) in 150-day-old serum and small yellow follicle (SYF) tissues. A histology examination, immunohistochemistry, and quantitative real-time PCR (qPCR) were used to explore the molecular mechanism of how some genes related to oxidative stress affect sexual maturation. The results showed that the secretion levels of LH, E2, and P4 in the CH2 group serum and SYF were higher than those in the WLS group. The results of the real-time PCR of all genes showed that the expression levels of cytochrome P450 family 11 subfamily A member 1, steroidogenic acute regulatory protein, follicle-stimulating hormone receptor, and cytochrome P450 family 19 subfamily A member 1 in the CH2 group were significantly higher than those in the WLS groups (p < 0.001). Untargeted metabolomics combined with multivariate statistical analysis was used to identify biomarkers of SYF tissues in the CH2 and WLS groups. A trajectory analysis of the principal component analysis (PCA) results showed that the samples within the group were clustered and that the samples were dispersed between the CH2 and the WLS groups, indicating that the results of the measured data were reliable and could be used for further research. Further analysis showed that a total of 319 metabolites in small yellow follicles of the CH2 and WLS groups were identified, among which 54 downregulated differential metabolites were identified. These 54 metabolites were found as potential CH2 biomarkers compared with WLS at 150 days, and the different expressions of L-arginine, L-prolinamide, (R)-4-hydroxymandelate, glutathione, and homovanillic acid were more significant. Twenty metabolic pathways were found when significantly differential metabolites were queried in the KEGG database. According to the impact values of the metabolic pathways, eighteen differential metabolites belonged to the mTOR signaling pathway, glutathione metabolism, ABC transporters, the cell ferroptosis pathway, and D-arginine and D-ornithine metabolism. Interestingly, we identified that the cell ferroptosis pathway played an important role in chicken follicle selection for the first time. The histology and immunohistochemistry of SYF showed that the number of granulosa cells increased in the CH2 groups and the expression levels of glutathione peroxidase 4, tumor protein p53, ribosomal protein S6 kinase, and sterol regulatory element binding protein 1 in the granulosa cell layer were upregulated in the CH2 group at the time of sexual maturation. Furthermore, we also speculated that the antioxidant system may play an indispensable role in regulating sexual maturity in chickens. Overall, our findings suggest differentially expressed metabolites and metabolic pathways between CH2 and WLS chickens, providing new insights into the initiation mechanism of sexual maturation.

Integration of transcriptome sequencing and whole genome resequencing reveal candidate genes in egg production of upright and pendulous-comb chickens

Egg production performance plays an important role in the poultry industry across the world. Previous studies have shown a great difference in egg production performance between pendulous-comb (PC) and upright-comb (UC) chickens. However, there are no reports to identify potential candidate genes for egg production in PC and UC chickens. In the present study, 1,606 laying chickens were raised, and the egg laid by individual chicken was collected for 100 d. Moreover, the expression level of estrogen and progesterone hormones was measured at the start-laying and peak-laying periods of hens. Besides, 4 PC and 4 UC chickens were selected at 217 d of age to perform transcriptome sequencing (RNA-seq) and whole genome resequencing (WGS) to screen the potential candidate genes of egg production. The results showed that PC chicken demonstrated better egg production performance (P < 0.05) and higher estrogen and progesterone hormone expression levels than UC chicken (P < 0.05). RNA-seq analysis showed that 341 upregulated and 1,036 downregulated differentially expressed genes (DEGs) were identified in the ovary tissues of PC and UC chickens. These DEGs were mainly enriched in protein-related, lipid-related, and nucleic acids-related biological processes including ribosome, peptide biosynthetic process, lipid transport terms, and catalytic activity acting on RNA which can significantly affect egg production in chicken. The enrichment results of WGS analysis were consistent with RNA-seq. Further, joint analysis of WGS and RNA-seq data was utilized to screen 30 genes and CAMK1D, CLSTN2, MAST2, PIK3C2G, TBC1D1, STK3, ADGRB3, and PPARGC1A were identified as potential candidate genes for egg production in PC and UC chickens. In summary, our study provides a wealth of information for a better understanding of the genetic and molecular mechanism for the future breeding of PC and UC chickens for egg production.

Interacting Networks of the Hypothalamic-Pituitary-Ovarian Axis Regulate Layer Hens Performance

Egg production is a vital biological and economic trait for poultry breeding. The 'hypothalamic-pituitary-ovarian (HPO) axis' determines the egg production, which affects the layer hens industry income. At the organism level, the HPO axis is influenced by the factors related to metabolic and nutritional status, environment, and genetics, whereas at the cellular and molecular levels, the HPO axis is influenced by the factors related to endocrine and metabolic regulation, cytokines, key genes, signaling pathways, post-transcriptional processing, and epigenetic modifications. MiRNAs and lncRNAs play a critical role in follicle selection and development, atresia, and ovulation in layer hens; in particular, miRNA is known to affect the development and atresia of follicles by regulating apoptosis and autophagy of granulosa cells. The current review elaborates on the regulation of the HPO axis and its role in the laying performance of hens at the organism, cellular, and molecular levels. In addition, this review provides an overview of the interactive network regulation mechanism of the HPO axis in layer hens, as well as comprehensive knowledge for successfully utilizing their genetic resources.

MicroRNAs with non-additive expression in the ovary of hybrid hens target genes enriched in key reproductive pathways that may influence heterosis for egg laying traits

Heterosis has been extensively exploited in chicken breeding to improve laying traits in commercial hybrid stock. However, the molecular mechanisms underlying it remains elusive. This study characterizes the miRNAome in the pre-hierarchical follicles of purebred and hybrid laying hens, and investigate the functions of miRNAs with non-additive expression in the pre-hierarchical follicles as they modulate heterosis for egg number and clutch size. To achieve that aim, White Leghorn and Rhode Island Red chicken lines were reciprocally crossed to generate hybrids. The crossbreds demonstrated heterosis for egg number and clutch size, and pre-hierarchical follicles from 4 birds of each genotype were collected at 53 weeks of age. Mode of miRNA expression was characterized after miRNA sequencing. A total of 50 miRNAs including 30 novel ones, were found to exhibit non-additive expression. Dominance was the predominant mode of expression exhibited by majority of the miRNAs. Functional analysis of target genes of the known miRNAs with non-additive expression revealed Gene Ontology terms related to regulation of transcription, metabolic processes and gene expression. KEGG and REACTOME pathways including hedgehog, cellular senescence, wnt, TGF-β, progesterone-mediated oocyte maturation, oocyte meiosis, GnRH signaling, signal transduction and generic transcription, which can be linked to primordial follicle activation, growth and ovulation, were significantly enriched by target genes of miRNAs with non-additive expression. Majority of the genes enriched in these biological pathways were targeted by gga-miR-19a, gga-miR-19b, gga-miR-375, gga-miR-135a, and gga-miR-7 and 7b, thus, revealing their synergistic roles in enhancing processes that could influence heterosis for egg number and clutch size in hybrid hens.

Grade follicles transcriptional profiling analysis in different laying stages in chicken

During follicular development, a series of key events such as follicular recruitment and selection are crucially governed by strict complex regulation. However, its molecular mechanisms remain obscure. To identify the dominant genes controlling chicken follicular development, the small white follicle (SWF), the small yellow follicle (SYF), and the large yellow follicle (LYF) in different laying stages (W22, W31, W51) were collected for RNA sequencing and bioinformatics analysis. There were 1866, 1211, and 1515 differentially expressed genes (DEGs) between SWF and SYF in W22, W31, and W51, respectively. 4021, 2295, and 2902 DEGs were respectively identified between SYF and LYF in W22, W31, and W51. 5618, 4016, and 4809 DEGs were respectively identified between SWF and LYF in W22, W31, and W51. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that extracellular matrix, extracellular region, extracellular region part, ECM-receptor interaction, collagen extracellular matrix, and collagen trimer were significantly enriched (P < 0.05). Protein–protein interaction analysis revealed that COL4A2, COL1A2, COL4A1, COL5A2, COL12A1, ELN, ALB, and MMP10 might be key candidate genes for follicular development in chicken. The current study identified dominant genes and pathways contributing to our understanding of chicken follicular development. 

Transcriptome analyses of potential regulators of pre- and post-ovulatory follicles in the pigeon (Columba livia)

To identify the dominant genes controlling follicular maturation, ovulation and regression for pigeon, we used RNA-seq to explore the gene expression profiles of pre- and post-ovulatory follicles of pigeon. We obtained total of 4.73million (96% of the raw data) high-quality clean reads, which could be aligned with 20282 genes. Gene expression profile analysis identified 1461 differentially expressed genes (DEGs) between the pre- (P4) and post-ovulatory follicles (P5). Of these, 843 genes were upregulated, and 618 genes were down-regulated. Furthermore, many DEGs were significantly enriched in some pathways closely related to follicle maturation, ovulation and regression, such as ECM-receptor interaction, vascular smooth muscle contraction, progesterone-mediated oocyte maturation, phagosome. Importantly, the DGEs in ECM-receptor interaction pathway included COL1A1 , COL1A2 , COL4A1 , COL4A2 , ITGA11 , ITGB3 and SDC3 , in the progesterone-mediated oocyte maturation pathway involved CDK1 , CDC25A , CCNB3 , CDC20 and Plk1 , and in the vascular smooth muscle contraction covered CALD1 , KCNMA1 , KCNMB1 , CACNA1 , ACTA2 , MYH10 , MYL3 , MYL6 , MYL9 , closely related to promoting follicular maturation and ovulation in pre-ovulatory follicles. Moreover, it seems that the lysosomal cathepsin family has a decisive role in the regression of early stage of post-ovulatory follicle. Taken together, these data enrich the research of molecular mechanisms of pigeon follicular activities at the transcriptional level and provide novel insight of breeding-related physiology for birds.

Differential proteomic analysis demonstrates follicle fluid participate immune reaction and protein translation in yak

Background

Ovarian follicle fluid (FF) as a microenvironment surrounding oocyte plays critical roles in physio-biochemical processes of follicle development and oocyte maturation. It is hypothesized that proteins in yak FF participate in the physio-biochemical pathways. The primary aims of this study were to find differentially expressed proteins (DEPs) between mature and immature FF, and to elucidating functions of the mature and immature FF in yak.

Results

The mature and immature FF samples were obtained from three healthy yaks that were nonpregnant, aged from four to five years, and free from any anatomical reproductive disorders. The FF samples were subjected to mass spectrometry with the isobaric tags for relative and absolute quantification (iTRAQ). The FF samples went through correlation analysis, principle component analysis, and expression pattern analysis based on quantification of the identified proteins. Four hundred sixty-three DEPs between mature and immature FF were identified. The DEPs between the mature and immature FF samples underwent gene ontology (GO), Kyoto encyclopedia of genes and genomes (KEGG), and protein-protein interaction (PPI) analysis. The DEPs highly expressed in the mature FF mainly took parts in the complement and coagulation cascades, defense response, acute-phase response, response to other organism pathways to avoid invasion of exogenous microorganisms. The complement activation pathway contains eight DEPs, namely C2, C5, C6, C7, C9, C4BPA, CFH, and MBL2. The three DEPs, CATHL4, CHGA, and PGLYRP1, take parts in defense response pathway to prevent invasion of exogenetic microorganism. The coagulation cascades pathway involves many coagulation factors, such as F7, F13A1, FGA, FGB, FGG, KLKB1, KNG1, MASP1, SERPINA1, and SERPIND1. While the DEPs highly expressed in the immature FF participated in protein translation, peptide biosynthetic process, DNA conformation change, and DNA geometric change pathways to facilitate follicle development. The translation pathway contains many ribosomal proteins, such as RPL3, RPL5, RPS3, RPS6, and other translation factors, such as EIF3J, EIF4G2, ETF1, MOV10, and NARS. The DNA conformation change and DNA geometric change involve nine DEPs, DDX1, G3BP1, HMGB1, HMGB2, HMGB3, MCM3, MCM5, MCM6, and RUVBL2. Furthermore, the expressed levels of the main DEPs, C2 and SERPIND1, were confirmed by western blot.

Conclusions

The differential proteomics revealed the up-regulated DEPs in mature FF take parts in immunoreaction to prevent invasion of microorganisms and the up-regulated DEPs in immature FF participate in protein synthesis, which may improve our knowledge of the follicular microenvironment and its biological roles for reproductive processes in yak. The DEPs, C2 and SERPIND1, can be considered as protein markers for mature yak follicle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-021-03097-0.

Dynamic transcriptome and chromatin architecture in granulosa cells during chicken folliculogenesis

Folliculogenesis is a complex biological process involving a central oocyte and its surrounding somatic cells. Three-dimensional chromatin architecture is an important transcription regulator; however, little is known about its dynamics and role in transcriptional regulation of granulosa cells during chicken folliculogenesis. We investigate the transcriptomic dynamics of chicken granulosa cells over ten follicular stages and assess the chromatin architecture dynamics and how it influences gene expression in granulosa cells at three key stages: the prehierarchical small white follicles, the first largest preovulatory follicles, and the postovulatory follicles. Our results demonstrate the consistency between the global reprogramming of chromatin architecture and the transcriptomic divergence during folliculogenesis, providing ample evidence for compartmentalization rearrangement, variable organization of topologically associating domains, and rewiring of the long-range interaction between promoter and enhancers. These results provide key insights into avian reproductive biology and provide a foundational dataset for the future in-depth functional characterization of granulosa cells.

The domestic chicken Gallus gallus domesticus is a classic model for the study of folliculogenesis. Here the authors integrate multi-omics analyses characterizing the dynamic transcriptome and chromatin architecture in granulosa cells during chicken folliculogenesis.

Systematic Analysis of Long Noncoding RNA and mRNA in Granulosa Cells during the Hen Ovulatory Cycle

Simple Summary

Chicken is one of the most economically important farm poultry, and providing many food products, such as meat and eggs for human consumption. However, follicle transcriptome studies in chickens with timepoints relating to changes in luteinizing hormone level remain unknown. In this study, the largest preovulatory follicle of chicken underwent the early, middle, and terminal stages of ovulatory cycle. Our work provides a comprehensive analysis of lncRNAs and mRNAs in chicken granulosa cells during the ovulatory cycle. A total of 12,479 mRNAs and 7528 lncRNAs were identified among the three stages. Thousands of lncRNAs were annotated, and the most differentially abundant genes were detected in the luteinizing hormone surge stage. Functional features of the lncRNAs and mRNAs at each stage were revealed, which was also associated with the changes in serum luteinizing hormone level. Especially, genes related to oxidative stress, steroids regulation, and inflammatory process were enriched in the luteinizing hormone surge stage, The comprehensive data generated in this study provides the foundation for future investigations to improve the reproductive performance of chickens and explore the mechanisms responsible for female ovarian diseases.

Abstract

Long non-coding RNAs (lncRNAs) and mRNAs are temporally expressed during chicken follicle development. However, follicle transcriptome studies in chickens with timepoints relating to changes in luteinizing hormone (LH) levels are rare. In this study, gene expression in Rohman layers was investigated at three distinct stages of the ovulatory cycle: zeitgeber time 0 (ZT0, 9:00 a.m.), zeitgeber time 12 (ZT12, 9:00 p.m.), and zeitgeber time 20 (ZT20, 5:00 a.m.) representing the early, middle, and LH surge stages, respectively, of the ovulatory cycle. Gene expression profiles were explored during follicle development at ZT0, ZT12, and ZT20 using Ribo-Zero RNA sequencing. The three stages were separated into two major stages, including the pre-LH surge and the LH surge stages. A total of 12,479 mRNAs and 7528 lncRNAs were identified among the three stages, and 4531, 523 differentially expressed genes (DEGs) and 2367, 211 differentially expressed lncRNAs (DELs) were identified in the ZT20 vs. ZT12, and ZT12 vs. ZT0, comparisons. Functional enrichment analysis revealed that genes involved in cell proliferation and metabolism processes (lipid-related) were mainly enriched in the ZT0 and ZT12 stages, respectively, and genes related to oxidative stress, steroids regulation, and inflammatory process were enriched in the ZT20 stage. These findings provide the basis for further investigation of the specific genetic and molecular functions of follicle development in chickens.

Identification of reliable reference genes for expression studies in the magnum of laying hens housed in cage and cage-free systems

Stress factors during poultry production can evoke changes in gene transcription and protein synthesis in the hen oviduct and could affect the internal and external egg quality. Studies of relative gene expression require the identification of the most stable reference genes for the quantitative polymerase chain reaction (qPCR) to investigate the reproductive tissues' response in laying hens kept in different production systems. The objective of this study was to determine the most stable reference genes of the magnum tissues of laying hens housed in two different production systems. Hy-Line Brown hens were reared under the same sanitary conditions until 15 weeks of age. Later on, they were transferred into two different production systems, conventional cage (CC) and cage free (CF), up to 82 weeks of age. At 50 and 60 weeks, a total of six hens from each production system were euthanized, and magnum samples were collected. The qPCR was used to determine the RNA transcription level of five reference genes, ACTB, 18S, GAPDH, MSX2 and HMBS. These genes were evaluated for transcript expression in magnum tissues by NormFinder, BestKeeper, geNorm and RefFinder software. The results indicated that the most stable gene in the CF housing system was HMBS in three of the algorithms and in the CC housing system was the 18S, and the best combination of reference genes was HMBS/GAPDH in CF and 18S/HMBS in CC. In conclusion, HMBS, 18S and GAPDH genes could be used together as reference genes for the normalization of the magnum tissues transcript expression of hens in CF and CC housing systems.

TGF-β1-induced collagen promotes chicken ovarian follicle development via an intercellular cooperative pattern

Follicle development is a complex process under strict regulation of diverse hormones and cytokines including transforming growth factor β (TGF-β) superfamily members. TGF-β is pivotal for the regulation of ovarian functions under physiological and pathological conditions. In this study, effect of TGF-β1 on chicken follicle development was examined through investigating the accumulation and action of collagen, an indispensable member of the extracellular matrix (ECM) involved in this process. The granulosa cells (GCs) and theca cells (TCs) were separated from growing follicles of the laying chicken for treatment of TGF-β1 and analysis of expression of ECM components and key proteins in intracellular signaling pathways. Results showed that collagen was mainly distributed in the follicular theca layer and was produced with the formation of the granulosa layer during ovarian development. Collagen accumulation increased with follicle growth and treatment of GCs with TGF-β1 elicited an increased expression of collagen. After production from GCs, collagen was transferred to the neighboring TCs to promote cell proliferation and inhibit apoptosis. Treatment of collagen remarkably increased expression of p-ERK, mitogen-activated protein kinase (MAPK), and p-MAPK, but treatment with hydroxylase inhibitor (to break collagen structure) reversed these alterations. In conclusion, during follicle growth collagen was secreted by GCs under TGF-β1 stimulation and was subsequently collaboratively transferred to neighboring TCs to increase cell proliferation and thus to promote follicle development via an intercellular cooperative pattern during development of chicken growing follicles.

Transcriptome Analysis of Long Noncoding RNAs and mRNAs in Granulosa Cells of Jinghai Yellow Chickens Illuminated With Red Light

Jinghai Yellow chickens are a new indigenous breed with a dual purpose in China, but their egg laying performance is limited. Compared with white light (WL), exposure to red light (RL) can improve the egg laying performance of hens. Herein, to elucidate the molecular mechanism by which RL affects the egg laying performance, RNA sequencing was used to analyze long noncoding RNAs (lncRNAs) and mRNAs from granulosa cells of small yellow follicles from Jinghai Yellow chickens in RL and WL groups. A total of 12,466 lncRNAs were identified among the assembled transcripts, of which 168 lncRNAs were significantly different between the RL and WL groups (101 downregulated and 67 upregulated). Additionally, 1182 differentially expressed mRNAs were identified (958 downregulated and 224 upregulated). Integrated network analysis demonstrated that numerous differential mRNAs were involved in follicular development through steroid hormone synthesis, oocyte meiosis, and the PI3K-Akt signaling pathway. The impact of lncRNAs on cis and trans target mRNAs indicates that some lncRNAs play important roles in follicular development of small yellow follicles. The results provide a starting point for studies aimed at understanding the molecular mechanisms by which monochromatic light affects follicular development and egg production in hens.