Transcriptomic Analysis of circRNAs and mRNAs Reveals a Complex Regulatory Network That Participate in Follicular Development in Chickens

Multi-tissue transcriptome profiling linked the association between tissue-specific circRNAs and the heterosis for feed intake and efficiency in chicken

Heterosis has been widely utilized in chickens. The nonadditive inheritance of genes contributes to this biological phenomenon. However, the role of circRNAs played in the heterosis is poorly determined. In this study, we observed divergent heterosis for residual feed intake (RFI) between 2 crossbreds derived from a reciprocal cross between White Leghorns and Beijing You chickens. Then, circRNA landscape for 120 samples covering the hypothalamus, liver, duodenum mucosa and ovary were profiled to elucidate the regulatory mechanisms of heterosis. We detected that a small proportion of circRNAs (7.83-20.35%) were additively and non-additively expressed, in which non-additivity was a major inheritance of circRNAs in the crossbreds. Tissue-specific expression of circRNAs was prevalent across 4 tissues. Weighted gene co-expression network analysis revealed circRNA-mRNA co-expression modules associated with feed intake and RFI in the hypothalamus and liver, and the co-expressed genes were enriched in oxidative phosphorylation pathway. We further identified 8 nonadditive circRNAs highly correlated with 16 nonadditive genes regulating negative heterosis for RFI in the 2 tissues. Circ-ITSN2 was validated in the liver tissue for its significantly positive correlation with PGPEP1L. Moreover, the bioinformatic analysis indicated that candidate circRNAs might be functioned by binding the microRNAs and interacting with the RNA binding proteins. The integration of multi-tissue transcriptome firstly linked the association between tissue-specific circRNAs and the heterosis for feed intake and efficiency in chicken, which provide novel insights into the molecular mechanism underlying heterosis for feed efficiency. The validated circRNAs can act as potential biomarkers for predicting RFI and its heterosis.

A Network of Circular RNA and MicroRNA Sequencing Provides Insights into Pigment Deposition of Changshun Blue Eggshell Chickens

Eggshell color plays important biological roles and attracts the attention of both egg retailers and researchers. However, whether non-coding RNAs are involved in pigment deposition among different eggshell colors remains unknown. In this study, RNA sequencing was used to analyse the uterine gland transcriptome (CircRNA and miRNA) of Changshun chicken blue-shell hens producing four different eggshell color eggs including dark blue PK(DB) and light blue (LB), dark brown and greenish (between blue and pink, DP) and pink (p). We found that miR-192-x, targeting SLC16a7, was expressed in DB, DP, and LB groups compared with the PK group, which indicates that miR-192-x may play a role in the blue eggshell color. KEGG and GO analyses showed that the "metabolic pathways" with targeted genes such BLVRA and HMOX1 were detected in dark and light blue color eggshell chickens, which confirms the different ratios of biliverdin and HO-1 involved in the deposition of blue color. As annotated by connectivity analysis, RASGRF1 and RASGRF2, belonging to the RASGRF family, are involved in the Ras signaling pathway, which plays an important role in cell growth, differentiation, metastasis and apoptosis. Our findings enrich the database of circRNA, miRNAs and genes for chicken uterine tissue, which will be useful in accelerating molecular selection for blue eggshell color layers.

Mining key circRNA-associated-ceRNA networks for milk fat metabolism in cows with varying milk fat percentages

BACKGROUND

Cow milk fat is an essential indicator for evaluating and measuring milk quality and cow performance. Growing research has identified the molecular functions of circular RNAs (circRNAs) necessary for mammary gland development and lactation in mammals.

METHOD

The present study analyzed circRNA expression profiling data in mammary epithelial cells (MECs) from cows with highly variable milk fat percentage (MFP) using differential expression analysis and weighted gene co-expression network analysis (WGCNA).

RESULTS

A total of 309 differentially expressed circRNAs (DE-circRNAs) were identified in the high and low MFP groups. WGCNA analysis revealed that the pink module was significantly associated with MFP (r = - 0.85, P = 0.007). Parental genes of circRNAs in this module were enriched mainly in lipid metabolism-related signaling pathways, such as focal adhesion, ECM-receptor interaction, adherens junction and AMPK. Finally, six DE-circRNAs were screened from the pink module: circ_0010571, circ_0007797, circ_0002746, circ_0003052, circ_0004319, and circ_0012840. Among them, circ_0002746, circ_0003052, circ_0004319, and circ_0012840 had circular structures and were highly expressed in mammary tissues. Subcellular localization revealed that these four DE-circRNAs may play a regulatory role in the mammary glands of dairy cows, mainly as competitive endogenous RNAs (ceRNAs). Seven hub target genes (GNB1, GNG2, PLCB1, PLCG1, ATP6V0C, NDUFS4, and PIGH) were obtained by constructing the regulatory network of their ceRNAs and then analyzed by CytoHubba and MCODE plugins in Cytoscape. Functional enrichment analysis revealed that these genes are crucial and most probable ceRNA regulators in milk fat metabolism.

CONCLUSIONS

Our study identified several vital circRNAs and ceRNAs affecting milk fat synthesis, providing new research ideas and a theoretical basis for cow lactation, milk quality, and breed improvement.

Analysis of circRNA-miRNA-mRNA regulatory network of embryonic gonadal development in Mulard duck

The aim of the study was to explore the regulatory mechanism of differences in embryonic gonadal development between intergeneric distance hybrid offspring Mulard ducks and parent ducks. The morphological differences gonadal tissues of Muscovy ducks, Pekin ducks and Mulard ducks at 12.5-day embryonic age were observed by sectioning and hematoxylin-eosin (HE) staining. Then followed by transcriptome sequencing to screen for gonadal development-related differentially expressed circRNAs and mRNAs to construct a competitive endogenous RNA (ceRNA) regulatory network. Finally, qRT-PCR and luciferase reporter system were used to verify the sequencing data and targeting relationship of ceRNA pairs. The results showed that the seminiferous tubule lumen of Mulard ducks was not obvious, while there were obvious seminiferous tubules and tubular structures in testis of Pekin ducks and Muscovy ducks, with number and shape indicating maturity. There were 18 upregulated circRNAs and 16 downregulated circRNAs in Mulard ducks and Pekin ducks, respectively, and 39 upregulated circRNAs and 1 downregulated circRNA in Mulard ducks and Muscovy ducks, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis found that genes involves in dorso-ventral axis formation, for example, neurogenic locus notch homolog protein 1 (NOTCH1), were significantly enriched (P < 0.05). The novel_circ_0002265-gga-miR-122-5p-PAFAH1B2 regulatory network was constructed. The qRT-PCR results showed that the sequencing results were reliable. The dual-luciferase reporter assay showed that gga-miR-122-5p exists binding site of circ_0002265 and PAFAH1B2, indicating circ_0002265-gga-miR-122-5p-PAFAH1B2 targeting relationship. In summary, the embryonic gonadal development of intergeneric hybrid Mulard ducks may be regulated by differentially expressed circRNAs and genes, such as novel_circ_0000519, novel_circ_0003537, NOTCH1, FGFR2, PAFAH1B1, and PAFAH1B2, among which circ_0002265-gga-miR-122-5p-PAFAH1B2 may participate in the targeted regulation of gonadal development in Mulard ducks. The findings of this study are helpful for analyzing the mechanism of embryonic gonadal development differences in avians.

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.

Estrogen Receptor Gene 1 (ESR1) Mediates Lipid Metabolism in Goose Hierarchical Granulosa Cells Rather than in Pre-Hierarchical Granulosa Cells

(1) Background: The role of estrogen receptor gene 1 (ESR1) in female reproduction and lipid metabolism has been extensively investigated. However, its contribution to lipid metabolism during the development of poultry follicles remains unclear. (2) Methods: This study aimed to explore the function of ESR1 via overexpressing (ESR1^ov) and interfering (ESR1^si) with its expression in pre-hierarchical granulosa cells (phGCs) and hierarchical granulosa cells (poGCs). (3) Results: We successfully cloned and obtained an 1866 bp segment of the full-length CDS region of the Sichuan white goose ESR1 gene. In phGCs of the ESR1^ov and ESR1^si groups, there were no significant changes compared to the control group. However, in poGCs, the ESR1^ov group exhibited decreased lipid deposition, triglycerides, and cholesterol compared to the control group, while the ESR1^si group showed increased lipid deposition, triglycerides, and cholesterol. The expression of APOB and PPARα was significantly reduced in the ESR1^ov group compared to the ESR1^ov-NC group. Moreover, significant changes in the expression of ACCα, DGAT1, SCD, CPT1, and ATGL were observed between the ESR1^si and ESR1^si-NC group. (4) Conclusions: These findings shed light on the function and molecular mechanism of ESR1 in lipid metabolism in goose poGCs, providing a better understanding of the physiological process of goose follicular development.

Effects of quercetin on granulosa cells from prehierarchical follicles by modulating MAPK signaling pathway in chicken

Quercetin (Que), widely found in a huge variety of plants, plays important roles in ovarian function. However, to data, there have been no reports about Que regulating granulosa cells (GCs) in prehierarchical follicles in chicken. Herein, GCs from follicles diameter from 4 to 8 mm in chicken were treated by Que in vitro culture to investigate how Que exerts its effect on follicular development. GCs treated by Que in concentrations of 10, 100, and 1,000 ng/mL were tested for cell proliferation and progesterone secretion. Eight cDNA libraries were constructed from GCs (4 samples per group) to explore transcriptome expression changes. The role of the MAPK/ERK signaling pathway was validated in this process. Treatment with 100 and 1,000 ng/mL levels of Que significantly promoted cell proliferation and progesterone secretion (P < 0.05). RNA-seq analysis data showed that 402 and 263 differentially expressed genes (DEGs) were up- and down-regulated, respectively. Functional enrichment analysis that the pathways related to follicular development included biosynthesis of amino acids, MAPK signaling pathway, and calcium signaling pathway. Notably, the function exerted in GCs of the different levels of Que was associated with the suppression of the MAPK pathway. In conclusion, our results proved that low levels of Que could promote MAPK signaling pathway, but high levels of Que inhibit MAPK signaling pathway in GCs from the prehierarchical follicles, promote cell proliferation and progesterone secretion, and benefit follicle selection.

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.

Transcriptomic analysis of ovarian follicles uncovers the crucial genes relevant to follicle selection and preovulatory hierarchy in hens

Follicle selection and preovulatory hierarchy of hen ovaries were important stages of follicle development and crucially determining egg-laying performance. The selected follicles with a higher expression level of follicle-stimulating hormone receptor (FSHR) mRNA that facilitates response to FSH, and rapidly develops into preovulatory follicles with distinctive characteristics of granulosa cells (GCs) proliferation and differentiation. Identification of the key genes involved in these developmental events is helpful for elucidation of the molecular mechanism underlying egg-laying traits in chicken and other domestic fowl. Herein, the comparative transcriptomic analysis of ovarian prehierarchical follicles before selection (BSF), follicles at selection stage (ASF), and hierarchical follicles (HF) were implemented in the Jilin Black chicken (JB) and Lohmann Brown layer (LB) with the divergences in their egg-laying performance by RNA-sequencing. The results showed that nine deferentially expressed genes (DEGs), including STMN4, FABP3, ROBO2, RSPO4, and DMRT1 were revealed between follicles BSF and ASF; and seventeen DEGs, such as SLC6A15, SLITRK3, PRKG2 and TMC3 were mined between ASF and HF. These two group DEGs being co-expressed between BSF and ASF, and between ASF and HF were compared and substantiated in the JB and LB layers, respectively. Furthermore, 10 signaling pathways, such as cAMP signaling, PPAR signaling pathway, AMPK(Adenosine 5'-monophosphate (AMP)-activated protein kinase) pathway, and estrogen signaling pathway were also identified. Moreover, the roles of two representative candidates ROBO2 and PRKG2 genes presented as downregulated mRNA expression pattern in the transcriptomic profiles were further verified in vitro. The results demonstrated that downregulation of ROBO2 or PRKG2 significantly increased the expression levels of FSHR mRNA and protein with the boosted expression of CCND1, STAR, and BCL-2, whereas remarkably inhibited the expression of Caspase-3, consequently, brought about the decrease of GC apoptosis in the ovarian follicles, but increase of GC proliferation and differentiation serving as the hallmarks for follicle selection. It indicated that ROBO2 and PRKG2 may play indispensable roles in follicle selection and preovulatory hierarchy of hen ovaries separately. Our findings provided a comparative transcriptomic evidence for clarifying the molecular mechanism of the follicle development underlying egg-laying traits in chicken.

Expression profile of long non-coding RNA in inner Mongolian cashmere goat with putative roles in hair follicles development

The hair follicle is a complex skin accessory organ, which determines hair growth. Long non-coding RNAs (lncRNAs) have been proven to play an important role in hair follicle development, but their specific mechanism is still unclear. In this study, high-throughput sequencing was used to obtain the expression profiles of lncRNA in the hair follicles of Inner Mongolian cashmere goats at different embryonic stages (45, 55, 65, and 75 days), and a total of 6,630 lncRNA were identified. According to the rules of hair follicle development, we combined miRNA and mRNA databases (published) and predicted lncRNA-miRNA, miRNA-mRNA, and lncRNA-mRNA interaction pairs in the 45 vs. 75 comparison group. We obtained 516 lncRNA-mRNA, 1,011 lncRNA-miRNA, and 7,411 miRNA-mRNA relationship pairs. Finally, target genes were analyzed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG), and it was found that they were mainly enriched in the Wnt signaling pathway and PI3K-Akt signaling pathway related to hair follicle development, indicating that lncRNA may interact with miRNA/mRNA to directly or indirectly regulate the expression of genes related to hair follicle development. Dual-luciferase reporter gene analysis showed that lncRNA MSTRG.1705.1 could bind to Chi-miR-1, while lncRNA MSTRG.11809.1 had no binding site for Chi-miR-433. In conclusion, this study aims to further analyze the molecular regulation mechanism of hair follicle development and to lay a theoretical foundation for revealing the regulation mechanism of cashmere hair follicle growth.

Construction and functional analysis of ceRNA regulatory network related to the development of secondary hair follicles in Inner Mongolia cashmere goats

Cashmere goat hair follicles are divided into primary hair follicles and secondary hair follicles. The primary hair follicles produce coarse hair, and the secondary hair follicles produce cashmere. The development of hair follicles is affected by a variety of signaling molecules and pathways. Studies have shown that non-coding RNAs are widely involved in the development of hair follicles of the goat, including small RNAs (miRNAs), long non-coding RNAs (lncRNA), and circular RNAs (circRNAs). In recent years, circRNAs, as a new type of circular closed non-coding RNAs, have attracted great attention due to their high stability. However, its regulatory effect on cashmere goat hair follicles mainly focuses on the periodic regulation of secondary hair follicles, and there is no report on the development of cashmere goat hair follicles during the fetal period. Therefore, this study was based on the circRNA, miRNA, and mRNA expression profiles obtained by whole-transcriptional sequencing of the skin tissue of the Inner Mongolia cashmere goats in the fetal period (days 45, 55, 65, and 75) and screening out the morphological changes of hair follicles at different periods. A total of 113 circRNAs related to the development of secondary hair follicles were present. According to the principle of the ceRNA regulatory network, a ceRNA regulatory network composed of 13 circRNAs, 21 miRNAs, and 110 mRNAs related to the development of secondary hair follicles was constructed. Then, qRT-PCR and Sanger sequencing identified circRNA2034, circRNA5712, circRNA888, and circRNA9127 were circRNAs. Next, the dual-luciferase reporter gene verified the targeting relationship of circRNA5712-miR-27b-3p-Dll4. In conclusion, this study constructed a ceRNA regulatory network for the development of cashmere goat secondary hair follicles, laying a foundation for the analysis of circRNAs regulating the morphogenesis and development of cashmere goat secondary hair follicles through the ceRNA mechanism.

Comprehensive analysis of differentially expressed profiles of mRNA, lncRNA, and miRNA of Yili geese ovary at different egg-laying stages

Background

The development of the ovaries is an important factor that affects egg production performance in geese. Ovarian development is regulated by genes that are expressed dynamically and stage-specifically. The transcriptome profile analysis on ovarian tissues of goose at different egg laying stages could provide an important basis for screening and identifying key genes regulating ovarian development.

Results

In this study, 4 ovary tissues at each breeding period of pre-laying (PP), laying (LP), and ceased-laying period (CP), respectively, with significant morphology difference, were used for RNA extraction and mRNAs, lncRNAs, and miRNAs comparison in Yili geese. CeRNA regulatory network was constructed for key genes screening. A total of 337, 1136, and 525 differentially expressed DE mRNAs, 466, 925, and 742 DE lncRNAs and 258, 1131 and 909 DE miRNAs were identified between PP and LP, between CP and LP, and between CP and PP groups, respectively. Functional enrichment analysis showed that the differentially expressed mRNAs and non-coding RNA target genes were mainly involved in the cell process, cytokine-cytokine receptor interaction, phagosome, calcium signaling pathway, steroid biosynthesis and ECM-receptor interaction. Differential genes and non-coding RNAs, PDGFRB, ERBB4, LHCGR, MSTRG.129094.34, MSTRG.3524.1 and gga-miR-145–5p, related to reproduction and ovarian development were highly enriched. Furthermore, lncRNA-miRNA-mRNA regulatory networks related to ovary development were constructed.

Conclusions

Our study found dramatic transcriptomic differences in ovaries of Yili geese at different egg-laying stages, and a differential lncRNA-miRNA-mRNA regulatory network related to cell proliferation, differentiation and apoptosis and involved in stromal follicle development were established and preliminarily validated, which could be regarded as a key regulatory pathway of ovarian development in Yili geese.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08774-4.

The opposite effects of VGLL1 and VGLL4 genes on granulosa cell proliferation and apoptosis of hen ovarian prehierarchical follicles

Transcription cofactors Vestigial like family (VGLL) members consisting of four homologs (VGLL1-4) are associated with cell growth and metastasis in mammals, among which VGLL1 gene has been documented to possess tumorigenic functions in various types of tumor, and VGLL4 acts as a new tumor suppressor; likewise several studies indicated that they potentially play a role in the regulation of ovary growth and function. However, the biological effects of chicken VGLL1 and VGLL4 on the proliferation, apoptosis, and steroidogenesis of the granulosa cells (GCs) during ovarian follicle development remain unknown now. This study found that VGLL1 and VGLL4 genes present divergent expression patterns of the transcripts in the GCs of various sized prehierarchical follicles (PFs) before follicle selection. Specific small interfering RNA (siRNA) was employed to elucidate the exact roles of VGLL1 and VGLL4 in regulating the PF development of the hen ovary. The results demonstrated that the mRNA expression levels of the steroidogenic-related enzyme steroidogenic acute regulatory protein (STAR) gene and the cell proliferation-related factors B-cell lymphoma-2 (BCL2), and cyclin D1 (CCND1) genes were significantly down-regulated in the cells with VGLL1 silence but remarkably up-regulated in the cells lacking VGLL4. Whereas the expression level of the cell apoptosis biomarker caspase-3 (CASP3) transcript was noticeably enhanced in the GCs without VGLL1 but significantly decreased in the GCs deprived of VGLL4. Further results showed that the siRNA-mediated silence of VGLL1 caused a significant increase in apoptosis with a reduction in the proliferation of GCs. Nevertheless, knockdown of VGLL4 resulted in a remarkable decrement in apoptosis but a memorable augment in proliferation of the GCs. Taken together, this study proved that VGLL1 promotes cell proliferation and steroidogenesis but inhibits apoptosis. In contrast, VGLL4 stimulates GC apoptosis while suppressing the GC proliferation and steroidogenesis in the hen ovarian follicles. We conluded that VGLL1 and VGLL4 affect oppositely the ovarian prehierarchical follicle development by the different regulatory manner in the GC proliferation and apoptosis of chicken ovary.

Transcriptome comparative analysis of ovarian follicles reveals the key genes and signaling pathways implicated in hen egg production

Background

Ovarian follicle development plays an important role in determination of poultry egg production. The follicles at the various developmental stages possess their own distinct molecular genetic characteristics and have different biological roles in chicken ovary development and function. In the each stage, several genes of follicle-specific expression and biological pathways are involved in the vary-sized follicular development and physiological events. Identification of the pivotal genes and signaling pathways that control the follicular development is helpful for understanding their exact regulatory functions and molecular mechanisms underlying egg-laying traits of laying hens.

Results

The comparative mRNA transcriptomic analysis of ovarian follicles at three key developmental stages including slow growing white follicles (GWF), small yellow follicles (SYF) of recruitment into the hierarchy, and differentiated large yellow follicles (LYF), was accomplished in the layers with lower and higher egg production. Totally, 137, 447, and 229 of up-regulated differentially expressed genes (DEGs), and 99, 97, and 157 of down-regulated DEGs in the GWF, SYF and LYF follicles, including VIPR1, VIPR2, ADRB2, and HSD17B1 were identified, respectively. Moreover, NDUFAB1 and GABRA1 genes, two most promising candidates potentially associated with egg-laying performance were screened out from the 13 co-expressed DEGs in the GWF, SYF and LYF samples. We further investigated the biological effects of NDUFAB1 and GABRA1 on ovarian follicular development and found that NDUFAB1 promotes follicle development by stimulating granulosa cell (GC) proliferation and decreasing cell apoptosis, increases the expression of CCND1 and BCL-2 but attenuates the expression of caspase-3, and facilitates steroidogenesis by enhancing the expression of STAR and CYP11A1. In contrast, GABRA1 inhibits GC proliferation and stimulates cell apoptosis, decreases the expression of CCND1, BCL-2, STAR, and CYP11A1 but elevates the expression of caspase-3. Furthermore, the three crucial signaling pathways such as PPAR signaling pathway, cAMP signaling pathway and neuroactive ligand-receptor interaction were significantly enriched, which may play essential roles in ovarian follicle growth, differentiation, follicle selection, and maturation.

Conclusions

The current study provided new molecular data for insight into the regulatory mechanism underlying ovarian follicle development associated with egg production in chicken.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08213-w.

A genome-wide scan to identify signatures of selection in two Iranian indigenous chicken ecotypes

Background

Various regions of the chicken genome have been under natural and artificial selection for thousands of years. The substantial diversity that exits among chickens from different geographic regions provides an excellent opportunity to investigate the genomic regions under selection which, in turn, will increase our knowledge about the mechanisms that underlie chicken diversity and adaptation. Several statistics have been developed to detect genomic regions that are under selection. In this study, we applied approaches based on differences in allele or haplotype frequencies (F_ST and hapFLK, respectively) between populations, differences in long stretches of consecutive homozygous sequences (ROH), and differences in allele frequencies within populations (composite likelihood ratio (CLR)) to identify inter- and intra-populations traces of selection in two Iranian indigenous chicken ecotypes, the Lari fighting chicken and the Khazak or creeper (short-leg) chicken.

Results

Using whole-genome resequencing data of 32 individuals from the two chicken ecotypes, approximately 11.9 million single nucleotide polymorphisms (SNPs) were detected and used in genomic analyses after quality processing. Examination of the distribution of ROH in the two populations indicated short to long ROH, ranging from 0.3 to 5.4 Mb. We found 90 genes that were detected by at least two of the four applied methods. Gene annotation of the detected putative regions under selection revealed candidate genes associated with growth (DCN, MEOX2 and CACNB1), reproduction (ESR1 and CALCR), disease resistance (S1PR1, ALPK1 and MHC-B), behavior pattern (AGMO, GNAO1 and PSEN1), and morphological traits (IHH and NHEJ1).

Conclusions

Our findings show that these two phenotypically different indigenous chicken populations have been under selection for reproduction, immune, behavioral, and morphology traits. The results illustrate that selection can play an important role in shaping signatures of differentiation across the genomic landscape of two chicken populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00664-9.

Expression Profiling and Functional Analysis of Circular RNAs in Inner Mongolian Cashmere Goat Hair Follicles

Background

Inner Mongolian cashmere goats have hair of excellent quality and high economic value, and the skin hair follicle traits of cashmere goats have a direct and important effect on cashmere yield and quality. Circular RNA has been studied in a variety of tissues and cells.

Result

In this study, high-throughput sequencing was used to obtain the expression profiles of circular RNA (circRNA) in the hair follicles of Inner Mongolian cashmere goats at different embryonic stages (45, 55, 65, and 75 days). A total of 21,784 circRNAs were identified. At the same time, the differentially expressed circRNA in the six comparison groups formed in the four stages were: d75vsd45, 59 upregulated and 33 downregulated DE circRNAs; d75vsd55, 61 upregulated and 102 downregulated DE circRNAs; d75vsd65, 32 upregulated and 33 downregulated DE circRNAs; d65vsd55, 67 upregulated and 169 downregulated DE circRNAs; d65vsd45, 96 upregulated and 63 downregulated DE circRNAs; and d55vsd45, 76 upregulated and 42 downregulated DE circRNAs. Six DE circRNA were randomly selected to verify the reliability of the sequencing results by quantitative RT-PCR. Subsequently, the circRNA corresponding host genes were analyzed by the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway. The results showed that the biological processes related to hair follicle growth and development enriched by GO mainly included hair follicle morphogenesis and cell development, and the signaling pathways related to hair follicle development included the Notch signaling pathway and NF-κB signaling pathway. We combined the DE circRNA of d75vsd45 with miRNA and mRNA databases (unpublished) to construct the regulatory network of circRNA–miRNA–mRNA, and formed a total of 102 pairs of circRNA–miRNA and 126 pairs of miRNA–mRNA interactions. The binding relationship of circRNA3236–chi-miR-27b-3p and circRNA3236–chi-miR-16b-3p was further verified by dual-luciferase reporter assays, and the results showed that circRNA3236 and chi-miR-27b-3p, and circRNA3236 and chi-miR-16b-3p have a targeted binding relationship.

Conclusion

To summarize, we established the expression profiling of circRNA in the fetal skin hair follicles of cashmere goats, and found that the host gene of circRNA may be involved in the development of hair follicles of cashmere goats. The regulatory network of circRNA–miRNA–mRNA was constructed and preliminarily verified using DE circRNAs.