Genome-wide differential expression of long noncoding RNAs and mRNAs in ovarian follicles of two different chicken breeds

The first chicken oocyte nucleus whole transcriptomic profile defines the spectrum of maternal mRNA and non-coding RNA genes transcribed by the lampbrush chromosomes

Lampbrush chromosomes, with their unusually high rate of nascent RNA synthesis, provide a valuable model for studying mechanisms of global transcriptome up-regulation. Here, we obtained a whole-genomic profile of transcription along the entire length of all lampbrush chromosomes in the chicken karyotype. With nuclear RNA-seq, we obtained information about a wider set of transcripts, including long non-coding RNAs retained in the nucleus and stable intronic sequence RNAs. For a number of protein-coding genes, we visualized their nascent transcripts on the lateral loops of lampbrush chromosomes by RNA-FISH. The set of genes transcribed on the lampbrush chromosomes is required for basic cellular processes and is characterized by a broad expression pattern. We also present the first high-throughput transcriptome characterization of miRNAs and piRNAs in chicken oocytes at the lampbrush chromosome stage. Major targets of predicted piRNAs include CR1 and long terminal repeat (LTR) containing retrotransposable elements. Transcription of tandem repeat arrays was demonstrated by alignment against the whole telomere-to-telomere chromosome assemblies. We show that transcription of telomere-derived RNAs is initiated at adjacent LTR elements. We conclude that hypertranscription on the lateral loops of giant lampbrush chromosomes is required for synthesizing large amounts of transferred to the embryo maternal RNA for thousands of genes.

Characterization of bovine long non-coding RNAs, OOSNCR1, OOSNCR2 and OOSNCR3, and their roles in oocyte maturation and early embryonic development

In mammals, early embryogenesis relies heavily on the regulation of maternal transcripts including protein-coding and non-coding RNAs stored in oocytes. In this study, the expression of three bovine oocyte expressed long non-coding RNAs (lncRNAs), OOSNCR1, OOSNCR2, and OOSNCR3, was characterized in somatic tissues, the ovarian follicle, and throughout early embryonic development. Moreover, the functional requirement of each transcript during oocyte maturation and early embryonic development was investigated using a siRNA-mediated knockdown approach. Tissue distribution analysis revealed that OOSNCR1, OOSNCR2 and OOSNCR3 are predominantly expressed in fetal ovaries. Follicular cell expression analysis revealed that these lncRNAs are highly expressed in the oocytes, with minor expression detected in the cumulus cells (CCs) and mural granulosa cells (mGCs). The expression for all three genes was highest during oocyte maturation, decreased at fertilization, and ceased altogether by the 16-cell stage. Knockdown of OOSNCR1, OOSNCR2 and OOSNCR3 in immature oocytes was achieved by microinjection of the cumulus-enclosed germinal vesicle (GV) oocytes with siRNAs targeting these lncRNAs. Knockdown of OOSNCR1, OOSNCR2 and OOSNCR3 did not affect cumulus expansion, but oocyte survival at 12 h post-insemination was significantly reduced. In addition, knockdown of OOSNCR1, OOSNCR2 and OOSNCR3 in immature oocytes resulted in a decreased rate of blastocyst development, and reduced expression of genes associated with oocyte competency such as nucleoplasmin 2 (NPM2), growth differentiation factor 9 (GDF9), bone morphogenetic protein 15 (BMP15), and JY-1 in MII oocytes. The data herein suggest a functional requirement of OOSNCR1, OOSNCR2, and OOSNCR3 during bovine oocyte maturation and early embryogenesis.

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.

Transcriptome identification and characterization of long non-coding RNAs in the ovary of hens at four stages

Long non-coding RNAs (lncRNAs) play important roles in transcriptional and post-transcriptional regulation. LncRNAs, which are defined as non-coding RNAs more than 200 bp in length, are involved in key biological processes, such as cell proliferation and differentiation, epigenetic regulation, and gene transcriptional translation. Recent studies have shown that lncRNAs also play major regulatory roles in the reproduction of mammals. However, knowledge of the roles of lncRNAs in the chicken ovary lacking. In this study, we performed RNA-seq analyses of ovarian tissue from Hy-Line brown laying hens at four physiological stages [15, 20, 30, and 68 weeks of age (W)]. We identified 657 lncRNA transcripts that were differentially expressed during ovarian development, the number of down-regulated lncRNAs was higher than the number of up-regulated lncRNAs during development. We predicted the cis and trans target genes of the DE lncRNAs and constructed a lncRNA-mRNA interaction network, which indicated that the DE genes (DEGs) and the target genes of the DE lncRNAs are mainly involved in signaling pathways associated with ovarian development, including oocyte meiosis, calcium signaling pathways, ECM-receptor interactions, and ribosome and focal adhesion. Overall, we found that twelve lncRNAs were strongly involved in ovarian development: LNC_013443, LNC_001029, LNC_005713, LNC_016762, ENSGALT00000101857, LNC_003913, LNC_013692, LNC_012219, LNC_004140, ENSGALT00000096941, LNC_009356, and ENSGALT00000098716. In summary, our study utilized RNA-seq analysis of hen ovaries to explore key lncRNAs involved in ovarian development and function. Furthermore, the comprehensive analysis identified the target genes of these lncRNAs providing a better understanding of the mechanisms underlying ovarian development in hens and a theoretical basis for further research.

Long non-coding RNA MSTRG.5970.28 regulates proliferation and apoptosis of goose follicle granulosa cells via the miR-133a-3p/ANOS1 pathway

The development of follicles in the ovaries is a critical determinant of poultry egg production. There are existing studies on the follicular development patterns in poultry, but the specific regulatory mechanisms still need further study. In a previous study, we identified long non-coding RNA (lncRNA) MSTRG.5970.28, anosmin 1 (ANOS1), and its predicted target miR-133a-3p that may be associated with goose ovary development. However, the function of MSTRG.5970.28 in goose granulosa cells and its regulatory mechanisms affecting granulosa cell proliferation and apoptosis have not been reported. In the present study, MSTRG.5970.28 and miR-133a-3p overexpression and interference vectors were constructed. Combined with reverse-transcription real-time quantitative PCR (RT-qPCR), a dual luciferase activity assay, Cell Counting Kit-8 (CCK-8), and flow cytometric analysis, we investigated the role of the MSTRG.5970.28-miR-133a-3p-ANOS1 axis in goose follicular granulosa cells and the associated regulatory mechanisms. MSTRG.5970.28 was found to be localized in the cytoplasm and its expression was influenced by reproductive hormones. The targeting relationship among MSTRG.5970.28, ANOS1, and miR-133a-3p were verified by a dual luciferase activity assay. CCK-8 and apoptosis assays showed that MSTRG.5970.28 inhibited the proliferation and promoted apoptosis of goose granulosa cells. The regulatory role of miR-133a-3p on granulosa cell proliferation and apoptosis was opposite to MSTRG.5970.28. We found that the proliferative and apoptotic effects of granulosa cells caused by MSTRG.5970.28 overexpression were attenuated by miR-133a-3p. MSTRG.5970.28 functions as a competitive endogenous RNA that regulates ANOS1 expression by sponging miR-133a-3p and thus exerts regulatory functions in granulosa cells. In sum, the present study identified lncRNA MSTRG.5970.28 as associated with goose ovary development, which affects the expression of ANOS1 by targeting miR-133a-3p, thereby influencing the proliferation and apoptosis of goose granulosa cells.

Long Non-Coding RNA BNIP3 Inhibited the Proliferation of Bovine Intramuscular Preadipocytes via Cell Cycle

The intramuscular fat (or marbling fat) content is an essential economic trait of beef cattle and improves the flavor and palatability of meat. Several studies have highlighted the correlation between long non-coding RNAs (lncRNAs) and intramuscular fat development; however, the precise molecular mechanism remains unknown. Previously, through a high-throughput sequencing analysis, we found a lncRNA and named it a long non-coding RNA BNIP3 (lncBNIP3). The 5' RACE and 3' RACE explored 1945 bp total length of lncBNIP3, including 1621 bp of 5'RACE, and 464 bp of 3'RACE. The nucleoplasmic separation and FISH results explored the nuclear localization of lncBNIP3. Moreover, the tissue expression of lncBNIP3 was higher in the longissimus dorsi muscle, followed by intramuscular fat. Furthermore, down-regulation of lncBNIP3 increased the 5-Ethynyl-2'- deoxyuridine (EdU)-EdU-positive cells. The flow cytometry results showed that the number of cells in the S phase was significantly higher in preadipocytes transfected with si-lncBNIP3 than in the control group (si-NC). Similarly, CCK8 results showed that the number of cells after transfection of si-lncBNIP3 was significantly higher than in the control group. In addition, the mRNA expressions of proliferative marker genes CyclinB1 (CCNB1) and Proliferating Cell Nuclear Antigen (PCNA) in the si-lncBNIP3 group were significantly higher than in the control group. The Western Blot (WB) results also showed that the protein expression level of PCNA transfection of si-lncBNIP3 was significantly higher than in the control group. Similarly, the enrichment of lncBNIP3 significantly decreased the EdU-positive cells in the bovine preadipocytes. The results of flow cytometry and CCK8 assay also showed that overexpression of lncBNIP3 inhibited the proliferation of bovine preadipocytes. In addition, the overexpression of lncBNIP3 significantly inhibited the mRNA expressions of CCNB1 and PCNA. The WB results showed that the overexpression of lncBNIP3 significantly inhibited the expression of the CCNB1 protein level. To further explore the mechanism of lncBNIP3 on the proliferation of intramuscular preadipocytes, RNA-seq was performed after interference with si-lncBNIP3, and 660 differentially expressed genes (DEGs) were found, including 417 up-regulated DEGs and 243 down-regulated DEGs. The KEGG pathway analysis showed that the cell cycle was the most significant pathway for the functional enrichment of DEGs, followed by the DNA replication pathway. The RT-qPCR quantified the expression of twenty DEGs in the cell cycle. Therefore, we speculated that lncBNIP3 regulated intramuscular preadipocyte proliferation through the cell cycle and DNA replication pathways. To further confirm this hypothesis, the cell cycle inhibitor Ara-C was used to inhibit DNA replication of the S phase in intramuscular preadipocytes. Herein, Ara-C and si-lncBNIP3 were simultaneously added to the preadipocytes, and the CCK8, flow cytometry, and EdU assays were performed. The results showed that the si-lncBNIP3 could rescue the inhibitory effect of Ara-C in the bovine preadipocyte proliferation. In addition, lncBNIP3 could bind to the promoter of cell division control protein 6 (CDC6), and down-regulation of lncBNIP3 promoted the transcription activity and the expression of CDC6. Therefore, the inhibitory effect of lncBNIP3 on cell proliferation might be understood through the cell cycle pathway and CDC6 expression. This study provided a valuable lncRNA with functional roles in intramuscular fat accumulation and revealed new strategies for improving beef quality.

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.

Profiling and Functional Analysis of long non-coding RNAs in yak healthy and atretic follicles

Yak is the livestock on which people live in plateau areas, but its fecundity is low. Follicular development plays a decisive role in yak reproductive performance. As an important regulatory factor, the expression of long non-coding RNA (lncRNAs) in yak follicular development and its regulatory mechanism remains unclear. To explore the differentially expressed lncRNAs between healthy and atretic follicular in yaks. We used RNA-seq to construct lncRNA, miRNA, and mRNA expression profiles in yak atretic and healthy follicles, and the RNA sequence results were identified by qPCR. In addition, the correlation of lncRNA and targeted mRNA was also analyzed by Starbase software. Moreover, lncRNA/miRNA/mRNA networks were constructed by Cytoscape software, and the network was verified by dual-luciferase analysis. A total of 682 novel lncRNAs, 259 bta-miRNAs, and 1704 mRNAs were identified as differentially expressed between healthy and atretic follicles. Among them, 135 mRNAs were positively correlated with lncRNA expression and 97 were negatively correlated, which may be involved in the yak follicular development. In addition, pathway enrichment analysis of differentially expressed lncRNA host genes by Kyoto Genome Encyclopedia (KEGG) showed that host genes were mainly involved in hormone secretion, granulosa cell apoptosis, and follicular development. In conclusion, we identified a series of novel lncRNAs, constructed the lncRNA ceRNA regulatory network, and provided comprehensive resources for exploring the role of lncRNAs in yak ovarian follicular development.

gga-miR-449b-5p Regulates Steroid Hormone Synthesis in Laying Hen Ovarian Granulosa Cells by Targeting the IGF2BP3 Gene

MiRNAs have been found to be involved in the regulation of ovarian function as important post-transcriptional regulators, including regulators of follicular development, steroidogenesis, cell atresia, and even the development of ovarian cancer. In this study, we evaluated the regulatory role of gga-miR-449b-5p in follicular growth and steroid synthesis in ovarian granulosa cells (GCs) of laying hens through qRT-PCR, ELISAs, western blotting and dual-luciferase reporter assays, which have been described in our previous study. We demonstrated that gga-miR-449b-5p was widely expressed in granulosa and theca layers of the different-sized follicles, especially in the granulosa layer. The gga-miR-449b-5p had no significant effect on the proliferation of GCs, but could significantly regulate the expression of key steroidogenesis-related genes (StAR and CYP19A1) (p < 0.01) and the secretion of P4 and E2 (p < 0.01 and p < 0.05). Further research showed that gga-miR-449b-5p could target IGF2BP3 and downregulate the mRNA and protein expression of IGF2BP3 (p < 0.05). Therefore, this study suggests that gga-miR-449b-5p is a potent regulator of the synthesis of steroid hormones in GCs by targeting the expression of IGF2BP3 and may contribute to a better understanding of the role of functional miRNAs in laying hen ovarian development.

Female Germ Cell Development in Chickens and Humans: The Chicken Oocyte Enriched Genes Convergent and Divergent with the Human Oocyte

The development of germ cells and other physiological events in the differentiated ovary of humans are highly conserved with several mammalian species, except for the differences in timing. However, comparative knowledge on this topic is very scarce with respect to humans and lower vertebrates, such as chickens. In chickens, female germ cells enter into meiosis around embryonic day (E) 15.5 and are arrested in meiotic prophase I as primary oocytes. The oocytes arrested in meiosis I are accumulated in germ-cell cysts; shortly after hatching, they are enclosed by flattened granulosa cells in order to form primordial follicles. In humans, the process of meiotic recombination in female germ cells begins in the 10–11th week of gestation, and primordial follicles are formed at around week 20. In this review, we comprehensively elucidate both the conservation and the species-specific differences between chickens and humans with respect to germ cell, oocyte, and follicle development. Importantly, we provide functional insights into a set of chicken oocyte enriched genes (from E16 to 1 week post-hatch) that show convergent and divergent expression patterns with respect to the human oocyte (from week 11 to 26).

The Novel Competing Endogenous Long Noncoding RNA SM2 Regulates Gonadotropin Secretion in the Hu Sheep Anterior Pituitary by Targeting the Oar-miR-16b/TGF-β/SMAD2 Signaling Pathway

Pituitary gonadotropins play a pivotal role in reproduction. Long noncoding RNAs (lncRNAs) have been identified as important regulators in the hypothalamic−pituitary−ovarian (HPO) axis associated with reproduction. However, the contributions of lncRNAs to pituitary gonadotropin secretion remain largely unknown. Therefore, this work was performed to uncover the functional mechanisms of the novel lncRNA TCONS_00083279 (lncRNA SM2) and its potential targeting pathway oar-miR-16b/TGF-beta/SMAD2, which is associated with gonadotropin secretion in sheep pituitary cells. In the present study, the lncRNA SM2 showed high expression levels in the sheep pituitary gland, and it was located in both the nucleus and the cytoplasm of pituitary cells. lncRNA SM2 knockdown inhibited pituitary cell proliferation and FSH and LH secretion. The function of the lncRNA SM2 was sponged by oar-miR-16b, and this regulated the growth and gonadotropin secretion of pituitary cells by modulating SMAD2, as shown by the dual-luciferase reporter assay. FSH and LH levels were both upregulated by SMAD2 overexpression. Moreover, the levels of the lncRNA SM2, SMAD2 and TGFR1, as well as FSH and LH, in sheep pituitary cells increased significantly under gonadotropin-releasing hormone (GnRH) stimulation (p < 0.05). This work illustrates that the lncRNA SM2 regulates gonadotropin secretion in the Hu sheep anterior pituitary by targeting the oar-miR-16b/TGF-β/SMAD2 signaling pathway, providing a valuable resource for understanding the molecular mechanisms underlying sheep reproduction.

Transcription Analysis for Core Networks of lncRNAs–mRNAs: Implication for Potential Role in Sterility of Crassostrea gigas

Simple Summary

This study reveals the expression profiles of lncRNA in the gonads of the Pacific oyster Crassostrea gigas. The potential function of lncRNAs was predicted in the case of antisense and cis-regulatory mechanisms based on their physical positions and their coexpression relationships in the case of trans regulation. Sterility-related DEGs and DELs were chosen for subsequent analysis, demonstrating that trans-regulatory lncRNAs might play a vital role in the gametogenesis of C. gigas. We constructed core networks of lncRNAs–mRNAs for triploid sterile females and hermaphrodites based on pathway results, in which 28 lncRNAs and their 54 trans-regulatory genes were detected. Among 28 sterility-specific lncRNAs, MSTRG.79882.3 and MSTRG.79882.4 for triploid sterile females and MSTRG.33704.1, MSTRG.63844.1, and MSTRG.5675.1 for hermaphrodites play the most significant role.

Abstract

Long noncoding RNA (lncRNA), a type of non-protein-coding transcript, is emerging as a crucial regulator of gene expression. However, few roles of lncRNA in the reproductive process of the Pacific oyster (Crassostrea gigas) have been defined, especially in the regulatory mechanism of sterile triploids gametogenesis. To uncover the potential role of lncRNA, the gonads of diploids, sterile triploids, and partially sterile triploids underwent RNA sequencing. A total of 9618 reliable lncRNAs were identified. The target relationship between lncRNA and mRNA was predicted based on cis, trans, and antisense regulation with bioinformatic software. We chose differentially expressed lncRNAs and mRNAs when sterile triploids were compared to partially sterile triploids and diploids for subsequent functional enrichment analysis. Findings revealed that trans-regulatory lncRNAs might play a significant role in the gametogenesis of C. gigas. Combining pathway results, we constructed core networks of lncRNAs–mRNAs for triploid sterile females and hermaphrodites. Fifty-four genes related to cell division, germline-cell maintenance, and glycogen metabolism were found to be associated with sterility. A total of 28 candidate lncRNAs were predicted to trans-regulate these genes. We speculated that MSTRG.79882.3 and MSTRG.79882.4 for triploid sterile females and MSTRG.33704.1, MSTRG.63844.1, and MSTRG.5675.1 for hermaphrodites were highly important as they were predicted to regulate more sterility-specific genes than others. Our work collectively identified sterility-related lncRNAs and implicated the potential mechanism of lncRNA-mediated regulation in the gametogenesis of sterile triploid oysters.

LncRNAs in domesticated animals: from dog to livestock species

Animal genomes are pervasively transcribed into multiple RNA molecules, of which many will not be translated into proteins. One major component of this transcribed non-coding genome is the long non-coding RNAs (lncRNAs), which are defined as transcripts longer than 200 nucleotides with low coding-potential capabilities. Domestic animals constitute a unique resource for studying the genetic and epigenetic basis of phenotypic variations involving protein-coding and non-coding RNAs, such as lncRNAs. This review presents the current knowledge regarding transcriptome-based catalogues of lncRNAs in major domesticated animals (pets and livestock species), covering a broad phylogenetic scale (from dogs to chicken), and in comparison with human and mouse lncRNA catalogues. Furthermore, we describe different methods to extract known or discover novel lncRNAs and explore comparative genomics approaches to strengthen the annotation of lncRNAs. We then detail different strategies contributing to a better understanding of lncRNA functions, from genetic studies such as GWAS to molecular biology experiments and give some case examples in domestic animals. Finally, we discuss the limitations of current lncRNA annotations and suggest research directions to improve them and their functional characterisation.

LncRNA lnc_13814 promotes the cells apoptosis in granulosa cells of duck by acting as apla-miR-145-4 sponge

Follicle development is a vital factor which determines the reproductive performance of poultry. Long noncoding RNAs (lncRNAs) have been reported to maintain animal reproductive function and play key roles in ovarian development and hormone secretion. But the regulatory mechanism of lncRNAs in duck follicle development has seldom been reported. In this study, to better explore the molecular mechanism of follicle development in ducks, the follicular lncRNA was sequenced and analyzed. A total of 9,551 lncRNAs were predicted in the duck follicles. Four hundred and forty-five lncRNAs were differentially expressed between the white follicles and yellow follicles. The results of our studies showed that lnc_13814 promoted cell apoptosis in duck GCs. Furthermore, the bioinformatics analysis results demonstrated that lnc_13814 was involved in a lncRNA-miRNA-mRNA coexpression network and it was observed to sponge two follicle-related miRNAs by a luciferase activity assay. Moreover, we found that overexpression of lnc_13814 significantly increased DNA damage inducible transcript 3 (DDIT3) expression and downregulated GCs apoptosis. Finally, we found that lnc_13814 directly binds to and inhibits apla-mir-145-4; then, lnc_13814 increases the expression of DDIT3 and up-regulates GCs apoptosis. Taken together, our findings demonstrate that lncRNAs have potential effects on duck ovarian follicles and lncRNAs may represent a new approach to understand follicular development.

Comparative analysis of hypothalamus transcriptome between laying hens with different egg-laying rates

Egg-laying performance is one of the most important economic traits in the poultry industry. Commercial layers can lay one egg almost every day during their peak-laying period. However, many Chinese indigenous chicken breeds show a relatively low egg-laying rate, even during their peak-laying period. To understand what makes the difference in egg production, we compared the hypothalamus transcriptome profiles of Lushi blue-shelled-egg chickens (LBS), a Chinese indigenous breed with low egg-laying rate and Rhode Island Red chickens (RIR), a commercial layer with relatively high egg-laying rate using RNA-seq. A total of 753 differentially expressed genes (DEGs) were obtained. Of these DEGs, 38 genes were enriched in 2 Gene Ontology (GO) terms, namely reproduction term and the reproductive process term, and 6 KEGG pathways, namely Wnt signaling pathway, Oocyte meiosis, GnRH signaling pathway, Thyroid hormone signaling pathway, Thyroid hormone synthesis and MAPK signaling pathway, which have been long known to be involved in egg production regulation. To further determine the core genes from the 38 DEGs, protein-protein interaction (PPI) network, co-expression network and transcriptional regulatory network analyses were carried out. After integrated analysis and experimental validation, 4 core genes including RAC1, MRE11A, MAP7 and SOX5 were identified as the potential core genes that are responsible for the laying-rate difference between the 2 breeds. These findings paved the way for future investigating the mechanism of egg-laying regulation and enriched the chicken reproductive regulation theory.

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.

Differential actions of diacylglycerol acyltransferase (DGAT) 1 and 2 in regulating lipid metabolism and progesterone secretion of goose granulosa cells

Accumulating evidence shows that granulosa cells within both mammalian and avian ovaries have the ability to synthesize fatty acids through de novo lipogenesis and to accumulate triglycerides essential for oocyte and ovarian development. However, very little is known about the exact roles of key genes involved in the lipid metabolic pathway in granulosa cells. The goal of this study was to investigate the differential actions of diacylglycerol acyltransferase (DGAT) 1 and 2, which are recognized as the rate-limiting enzymes catalyzing the last step of triglyceride biosynthesis, in regulating lipid metabolism and steroidogenesis in granulosa cells of goose follicles at different developmental stages. It was observed that the mRNAs encoding DGAT1 and DGAT2 were ubiquitous in all examined granulosa cell layers but exhibited distinct expression profiles during follicle development. Notably, the mRNA levels of DGAT1, DGAT2, FSHR, LHR, STAR, CYP11A1, and 3βHSD remained almost constant in all except for 1-2 follicles within the 8-10 mm cohort, followed by an acute increase/decrease in the F5 follicles. At the cellular level, siRNA-mediated downregulation of DGAT1 or DGAT2 did not change the amount of lipids accumulated in both undifferentiated- and differentiated granulosa cells, while overexpression of DGAT2 promoted lipid accumulation and expression of lipogenic-related genes in these cells. Meanwhile, we found that interfering DGAT2 had no effect but interfering DGAT1 or overexpressing DGAT2 stimulated progesterone secretion in undifferentiated granulosa cells; in contrast, interference or overexpression of DGAT1/2 failed to change progesterone levels in differentiated granulosa cells but differently modulated expression of steroidogenic-related genes. Therefore, it could be concluded that DGAT1 is less efficient than DGAT2 in promoting lipid accumulation in both undifferentiated- and differentiated granulosa cells and that DGAT1 negatively while DGAT2 positively regulates progesterone production in undifferentiated granulosa cells.

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

Follicular development plays a key role in poultry reproduction, affecting clutch traits and thus egg production. Follicular growth is determined by granulosa cells (GCs), theca cells (TCs), and oocyte at the transcription, translation, and secretion levels. With the development of bioinformatic and experimental techniques, non-coding RNAs have been shown to participate in many life events. In this study, we investigated the transcriptomes of GCs and TCs in three different physiological stages: small yellow follicle (SYF), smallest hierarchical follicle (F6), and largest hierarchical follicle (F1) stages. A differential expression (DE) analysis, weighted gene co-expression network analysis (WGCNA), and bioinformatic analyses were performed. A total of 18,016 novel circular RNAs (circRNAs) were detected in GCs and TCs, 8127 of which were abundantly expressed in both cell types. and more circRNAs were differentially expressed between GCs and TCs than mRNAs. Enrichment analysis showed that the DE transcripts were mainly involved in cell growth, proliferation, differentiation, and apoptosis. In the WGCNA analysis, we identified six specific modules that were related to the different cell types in different stages of development. A series of central hub genes, including MAPK1, CITED4, SOD2, STC1, MOS, GDF9, MDH1, CAPN2, and novel_circ0004730, were incorporated into a Cytoscape network. Notably, using both DE analysis and WGCNA, ESR1 was identified as a key gene during follicular development. Our results provide valuable information on the circRNAs involved in follicle development and identify potential genes for further research to determine their roles in the regulation of different biological processes during follicle growth.

Comparative Transcriptome Analysis Suggests Key Roles for 5-Hydroxytryptamlne Receptors in Control of Goose Egg Production

To date, research on poultry egg production performance has only been conducted within inter or intra-breed groups, while those combining both inter- and intra-breed groups are lacking. Egg production performance is known to differ markedly between Sichuan white goose (Anser cygnoides) and Landes goose (Anser anser). In order to understand the mechanism of egg production performance in geese, we undertook this study. Here, 18 ovarian stromal samples from both Sichuan white goose and Landes goose at the age of 145 days (3 individuals before egg production initiation for each breed) and 730 days (3 high- and low egg production individuals during non-laying periods for each breed) were collected to reveal the genome-wide expression profiles of ovarian mRNAs and lncRNAs between these two geese breeds at different physiological stages. Briefly, 58, 347, 797, 777, and 881 differentially expressed genes (DEGs) and 56, 24, 154, 105, and 224 differentially expressed long non-coding RNAs (DElncRNAs) were found in LLD vs. HLD (low egg production Landes goose vs. high egg production Landes goose), LSC vs. HSC (low egg production Sichuan White goose vs. high egg production Sichuan white goose), YLD vs. YSC (young Landes goose vs. young Sichuan white goose), HLD vs. HSC (high egg production Landes goose vs. high egg production Sichuan white goose), and LLD vs. LSC (low egg production Landes goose vs. low egg production Sichuan white goose) groups, respectively. Functional enrichment analysis of these DEGs and DElncRNAs suggest that the "neuroactive ligand-receptor interaction pathway" is crucial for egg production, and particularly, members of the 5-hydroxytryptamine receptor (HTR) family affect egg production by regulating ovarian metabolic function. Furthermore, the big differences in the secondary structures among HTR1F and HTR1B, HTR2B, and HTR7 may lead to their different expression patterns in goose ovaries of both inter- and intra-breed groups. These results provide novel insights into the mechanisms regulating poultry egg production performance.

Genome-wide analysis of differentially expressed profiles of mRNAs, lncRNAs and circRNAs in chickens during Eimeria necatrix infection

BACKGROUND

Eimeria necatrix, the most highly pathogenic coccidian in chicken small intestines, can cause high morbidity and mortality in susceptible birds and devastating economic losses in poultry production, but the underlying molecular mechanisms in interaction between chicken and E. necatrix are not entirely revealed. Accumulating evidence shows that the long-non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) are key regulators in various infectious diseases. However, the expression profiles and roles of these two non-coding RNAs (ncRNAs) during E. necatrix infection are still unclear.

METHODS

The expression profiles of mRNAs, lncRNAs and circRNAs in mid-segments of chicken small intestines at 108 h post-infection (pi) with E. necatrix were analyzed by using the RNA-seq technique.

RESULTS

After strict filtering of raw data, we putatively identified 49,183 mRNAs, 818 lncRNAs and 4153 circRNAs. The obtained lncRNAs were classified into four types, including 228 (27.87%) intergenic, 67 (8.19%) intronic, 166 (20.29%) anti-sense and 357 (43.64%) sense-overlapping lncRNAs; of these, 571 were found to be novel. Five types were also predicted for putative circRNAs, including 180 exonic, 54 intronic, 113 antisense, 109 intergenic and 3697 sense-overlapping circRNAs. Eimeria necatrix infection significantly altered the expression of 1543 mRNAs (707 upregulated and 836 downregulated), 95 lncRNAs (49 upregulated and 46 downregulated) and 13 circRNAs (9 upregulated and 4 downregulated). Target predictions revealed that 38 aberrantly expressed lncRNAs would cis-regulate 73 mRNAs, and 1453 mRNAs could be trans-regulated by 87 differentially regulated lncRNAs. Additionally, 109 potential sponging miRNAs were also identified for 9 circRNAs. GO and KEGG enrichment analysis of target mRNAs for lncRNAs, and sponging miRNA targets and source genes for circRNAs identified associations of both lncRNAs and circRNAs with host immune defense and pathogenesis during E. necatrix infection.

CONCLUSIONS

To the best of our knowledge, the present study provides the first genome-wide analysis of mRNAs, lncRNAs and circRNAs in chicken small intestines infected with E. necatrix. The obtained data will offer novel clues for exploring the interaction mechanisms between chickens and Eimeria spp.

Dynamic expression and functional analysis of circRNA in granulosa cells during follicular development in chicken

Background

Circular RNA (circRNA) is a type of noncoding RNA involved in a variety of biological processes, especially in post-transcriptional regulation. The granulosa cells of follicles play a determining role in ovarian development. However, the function of circRNA in chicken follicles is unclear. To better understand the molecular mechanism underlying follicular development and granulosa cell function, we performed a strategy of second-generation sequencing and linear RNA depletion for granulosa cells from small yellow follicles (SYF, 5–8 mm), the smallest hierarchal follicles (F6, 9–12 mm), and the largest hierarchal follicles (F1, ~ 40 mm).

Results

We predicted a total of 11,642 circRNAs that distributed on almost all chromosomes. The majority of the splice lengths of circRNAs were 200–500 nt and mainly produced from intron and CDS regions. During follicle growth, differentially expressed (DE) circRNAs showed dynamic changes which were tissue- and stage-specific. The host genes of DE circRNAs were functionally enriched in GTPase activity and several pathways involved in reproduction. Moreover, bioinformatic prediction analysis for circRalGPS2 demonstrated that circRNAs from the same genes may share common miRNA to act as a sponge. The predicted target genes were enriched in various biological processes including cognition, cell communication, and regulation of signaling, and several pathways related to reproduction such as tight junction, oocyte meiosis, progesterone-mediated oocyte maturation, and GnRH signaling.

Conclusions

This study provides a starting point for further experimental investigations into chicken circRNAs and casts a light on the understanding of follicle development.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5462-2) contains supplementary material, which is available to authorized users.