Systematic identification of long intergenic non-coding RNAs expressed in bovine oocytes

Integrative analysis of RNA-seq and Ribo-seq reveals that lncRNA-GRN regulates chicken follicular atresia through miR-103-3p/FBXW7 axis and encoding peptide

Follicular atresia in chickens seriously reduced the egg production and economic benefits of chickens. LncRNA plays a key role in the process of follicular atresia. In this study, RNA-seq and Ribo-seq were performed on normal and atretic follicles of Dahen broilers to screen out lncRNAs that may regulate follicle atresia, and to study the molecular mechanisms of their regulation. GRN granulin precursor (lncGRN, ID: 101748909) was highly expressed in atretic follicles with translational ability. A molecular regulatory network of lncGRN/miR-103-3p/FBXW7 was constructed through bioinformatics analysis and dual luciferase reporting. LncGRN promoted the expression of FBXW7 by adsorption of miR-103-3p, thereby inhibiting the proliferation of chicken granulosa cells (GCs), promoting apoptosis of chicken GCs and inhibiting steroid hormone synthesis thus induced follicular atresia. Meanwhile, we also found a micropeptide named GRN-122aa derived by lncGRN which can promote follicular atresia. In conclusion, our study found that lncGRN promoted follicular atresia through the lncGRN/miR-103-3p/FBXW7 axis and the translation micropeptide GRN-122aa. This study provided new insight into the post-transcriptional regulation mechanism of lncGRN suggesting that lncGRN may act as a potential to regulate chicken follicle development, and provided a theoretical argument for further improving the egg production of chickens through molecular breeding.

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.

Characterization of agouti-signaling protein (ASIP) in the bovine ovary and throughout early embryogenesis

The oocyte expresses certain genes during folliculogenesis to regulate the acquisition of oocyte competence. Oocyte competence, or oocyte quality, is directly related to the ability of the oocyte to result in a successful pregnancy following fertilization. Presently, approximately 40 % of bovine embryos will develop to the blastocyst stage in vitro. Characterization of factors regulating these processes is crucial to improve the efficiency of bovine in vitro embryo production. We demonstrated that the secreted protein, agouti-signaling protein (ASIP) is highly abundant in the bovine oocyte and aimed to characterize its spatiotemporal expression profile in the ovary and throughout early embryonic development. In addition to oocyte expression, ASIP was detected in granulosa, cumulus, and theca cells isolated from antral follicles. Both gene expression data and immunofluorescent staining indicated ASIP declines with oocyte maturation which may indicate a potential role for ASIP in the attainment of oocyte competence. Microinjection of zygotes using small interfering RNA targeting ASIP led to a 16 % reduction in the rate of development to the blastocyst stage. Additionally, we examined potential ASIP signaling mechanisms through which ASIP may function to establish oocyte developmental competence. The expression of melanocortin receptor 3 and 4 and the coreceptor attractin was detected in the oocyte and follicular cells. The addition of cortisol during in vitro maturation was found to increase significantly oocyte ASIP levels. In conclusion, these results suggest a functional role for ASIP in promoting oocyte maturation and subsequent embryonic development, potentially through signaling mechanisms involving cortisol.

Identification of novel cattle (Bos taurus) genes and biological insights of their function in pre-implantation embryo development

BACKGROUND

Appropriate regulation of genes expressed in oocytes and embryos is essential for acquisition of developmental competence in mammals. Here, we hypothesized that several genes expressed in oocytes and pre-implantation embryos remain unknown. Our goal was to reconstruct the transcriptome of oocytes (germinal vesicle and metaphase II) and pre-implantation cattle embryos (blastocysts) using short-read and long-read sequences to identify putative new genes.

RESULTS

We identified 274,342 transcript sequences and 3,033 of those loci do not match a gene present in official annotations and thus are potential new genes. Notably, 63.67% (1,931/3,033) of potential novel genes exhibited coding potential. Also noteworthy, 97.92% of the putative novel genes overlapped annotation with transposable elements. Comparative analysis of transcript abundance identified that 1,840 novel genes (recently added to the annotation) or potential new genes were differentially expressed between developmental stages (FDR < 0.01). We also determined that 522 novel or potential new genes (448 and 34, respectively) were upregulated at eight-cell embryos compared to oocytes (FDR < 0.01). In eight-cell embryos, 102 novel or putative new genes were co-expressed (|r|> 0.85, P < 1 × 10-8) with several genes annotated with gene ontology biological processes related to pluripotency maintenance and embryo development. CRISPR-Cas9 genome editing confirmed that the disruption of one of the novel genes highly expressed in eight-cell embryos reduced blastocyst development (ENSBTAG00000068261, P = 1.55 × 10-7).

CONCLUSIONS

Our results revealed several putative new genes that need careful annotation. Many of the putative new genes have dynamic regulation during pre-implantation development and are important components of gene regulatory networks involved in pluripotency and blastocyst formation.

LncRNA S100PBP promotes proliferation and steroid hormone synthesis of granulosa cells by sponging MiR-2285bc-BMPR2 in bovine†

In bovine follicular development, the proliferation of bovine granulosa cells affects follicular selection, atresia, and cystic follicle formation. When cystic follicles appear on the ovaries, granulosa cells stop proliferating, resulting in the reduction of granulosa cells layer. In our previous study, the whole transcriptome sequencing revealed that Bone morphogenetic protein receptor 2 (BMPR2) was differentially expressed between cystic and normal follicular granulosa cells. We speculated that long noncoding RNA may act as competing endogenous RNA targeting microRNAs and then regulating the expression of BMPR2 and the function of granulosa cells, thereby affecting follicular development and cyst formation. In this study, the results elucidated that long noncoding RNA S100PBP (NONBTAT011846.2) directly bound miR-2285bc, which targeted in the BMPR2 3'-UTR. miR-2285bc suppresses granulosa cells proliferation by downregulating BMPR2 expression. Furthermore, long noncoding RNA S100PBP was silenced by small interfering RNA, and long noncoding RNA S100PBP regulated BMPR2 expression by sponging miR-2285bc investigated through cross-verification. When small interfering RNA of long noncoding RNA S100PBP was transfected into granulosa cells, the results revealed similar molecular changes as those transfected with miR-2285bc mimics. Silencing long noncoding RNA S100PBP or overexpressing miR-2285bc altered the expressions of some follicular development-related genes, which could be related to follicular cyst occurrence. In conclusion, our findings support that long noncoding RNA S100PBP regulates the expression of BMPR2 through sponge miR-2285bc, promotes the proliferation of granulosa cells, inhibits their apoptosis, and increases the synthesis and secretion of follicular steroid hormones, thus promoting the development of bovine follicles.

Comprehensive transcriptomic analysis of long non-coding RNAs in bovine ovarian follicles and early embryos

Long non-coding RNAs (lncRNAs) have been the subject of numerous studies over the past decade. First thought to come from aberrant transcriptional events, lncRNAs are now considered a crucial component of the genome with roles in multiple cellular functions. However, the functional annotation and characterization of bovine lncRNAs during early development remain limited. In this comprehensive analysis, we review lncRNAs expression in bovine ovarian follicles and early embryos, based on a unique database comprising 468 microarray hybridizations from a single platform designed to target 7,724 lncRNA transcripts, of which 5,272 are intergenic (lincRNA), 958 are intronic, and 1,524 are antisense (lncNAT). Compared to translated mRNA, lncRNAs have been shown to be more tissue-specific and expressed in low copy numbers. This analysis revealed that protein-coding genes and lncRNAs are both expressed more in oocytes. Differences between the oocyte and the 2-cell embryo are also more apparent in terms of lncRNAs than mRNAs. Co-expression network analysis using WGCNA generated 25 modules with differing proportions of lncRNAs. The modules exhibiting a higher proportion of lncRNAs were found to be associated with fewer annotated mRNAs and housekeeping functions. Functional annotation of co-expressed mRNAs allowed attribution of lncRNAs to a wide array of key cellular events such as meiosis, translation initiation, immune response, and mitochondrial related functions. We thus provide evidence that lncRNAs play diverse physiological roles that are tissue-specific and associated with key cellular functions alongside mRNAs in bovine ovarian follicles and early embryos. This contributes to add lncRNAs as active molecules in the complex regulatory networks driving folliculogenesis, oogenesis and early embryogenesis all of which are necessary for reproductive success.

Transcriptome Profiling of the Liver in Nellore Cattle Phenotypically Divergent for RFI in Two Genetic Groups

The identification and selection of genetically superior animals for residual feed intake (RFI) could enhance productivity and minimize environmental impacts. The aim of this study was to use RNA-seq data to identify the differentially expressed genes (DEGs), known non-coding RNAs (ncRNAs), specific biomarkers and enriched biological processes associated with RFI of the liver in Nellore cattle in two genetic groups. In genetic group 1 (G1), 24 extreme RFI animals (12 low RFI (LRFI) versus 12 high RFI (HRFI)) were selected from a population of 60 Nellore bulls. The RNA-seq of the samples from their liver tissues was performed using an Illumina HiSeq 2000. In genetic group 2 (G2), 20 samples of liver tissue of Nellore bulls divergent for RFI (LRFI, n = 10 versus HRFI, n = 10) were selected from 83 animals. The raw data of the G2 were chosen from the ENA repository. A total of 1811 DEGs were found for the G1 and 2054 for the G2 (p-value ≤ 0.05). We detected 88 common genes in both genetic groups, of which 33 were involved in the immune response and in blocking oxidative stress. In addition, seven (B2M, ADSS, SNX2, TUBA4A, ARHGAP18, MECR, and ABCF3) possible gene biomarkers were identified through a receiver operating characteristic analysis (ROC) considering an AUC > 0.70. The B2M gene was overexpressed in the LRFI group. This gene regulates the lipid metabolism protein turnover and inhibits cell death. We also found non-coding RNAs in both groups. MIR25 was up-regulated and SNORD16 was down-regulated in the LRFI for G1. For G2, up-regulated RNase_MRP and SCARNA10 were found. We highlight MIR25 as being able to act by blocking cytotoxicity and oxidative stress and RMRP as a blocker of mitochondrial damage. The biological pathways associated with RFI of the liver in Nellore cattle in the two genetic groups were for energy metabolism, protein turnover, redox homeostasis and the immune response. The common transcripts, biomarkers and metabolic pathways found in the two genetic groups make this unprecedented work even more relevant, since the results are valid for different herds raised in different ways. The results reinforce the biological importance of these known processes but also reveal new insights into the complexity of the liver tissue transcriptome of Nellore cattle.

Identification and validation of ferroptosis-related lncRNA signatures as a novel prognostic model for glioma

Background: Ferroptosis is a newly discovered form of regulated cell death with distinct properties and recognizing functions involved in physical conditions or various diseases, including cancers. However, the relationship between gliomas and ferroptosis-related lncRNAs (FRLs) remains unclear.

Methods: We collected a total of 1850 samples from The Cancer Genome Atlas (TCGA) and Genotype Tissue Expression (GTEX) databases, including 698 tumor and 1,152 normal samples. A list of ferroptosis-related genes was downloaded from the Ferrdb website. Differentially expressed FRLs (DEFRLS) were analyzed using the “limma” package in R software. Subsequently, prognosis-related FRLs were obtained by univariate Cox analysis. Finally, a prognostic model based on the 3 FRLs was constructed using Cox regression analysis with the least absolute shrinkage and selection operator (LASSO) algorithm. The prognostic power of the model was assessed using receiver operating characteristic (ROC) curve analysis and Kaplan-Meier (K-M) survival curve analysis. In addition, we further explored the relationship of the immune landscape and somatic mutations to prognostic model characteristics. Finally, we validated the function of LINC01426 in vitro.

Results: We successfully constructed a 3-FRLs signature and classified glioma patients into high-risk and low-risk groups based on the risk score calculated from this signature. Compared with traditional clinicopathological features [age, sex, grade, isocitrate dehydrogenase (IDH) status], the prognostic accuracy of this model is more stable and stronger. Additionally, the model had stable predictive power for overall survival over a 5-year period. In addition, we found significant differences between the two groups in cellular immunity, the numbers of many immune cells, including NK cells, CD4+, CD8+ T-cells, and macrophages, and the expression of many immune-related genes. Finally, the two groups were also significantly different at the level of somatic mutations, especially in glioma prognosis-related genes such as IDH1 and ATRX, with lower mutation rates in the high-risk group leading to poorer prognosis. Finally, we found that the ferroptosis process of glioma cells was inhibited after knocking down the expression of LINC01426.

Conclusion: The proposed 3-FRL signature is a promising biomarker for predicting prognostic features in glioma patients.

Identification of Long Noncoding RNAs Involved in Eyelid Pigmentation of Hereford Cattle

Several ocular pathologies in cattle, such as ocular squamous cell carcinoma and infectious keratoconjunctivitis, have been associated with low pigmentation of the eyelids. The main objective of this study was to analyze the transcriptome of eyelid skin in Hereford cattle using strand-specific RNA sequencing technology to characterize and identify long noncoding RNAs (lncRNAs). We compared the expression of lncRNAs between pigmented and unpigmented eyelids and analyzed the interaction of lncRNAs and putative target genes to reveal the genetic basis underlying eyelid pigmentation in cattle. We predicted 4,937 putative lncRNAs mapped to the bovine reference genome, enriching the catalog of lncRNAs in Bos taurus. We found 27 differentially expressed lncRNAs between pigmented and unpigmented eyelids, suggesting their involvement in eyelid pigmentation. In addition, we revealed potential links between some significant differentially expressed lncRNAs and target mRNAs involved in the immune response and pigmentation. Overall, this study expands the catalog of lncRNAs in cattle and contributes to a better understanding of the biology of eyelid pigmentation.

Identification and Validation of Ferroptosis-Related LncRNA Signatures as a Novel Prognostic Model for Colon Cancer

Background

Ferroptosis is a newly defined form of programmed cell death that plays an important role in many cancers. However, ferroptosis-related lncRNAs (FRLs) involved in the regulation of colon cancer are not thoroughly understood. This study aimed to identify a prognostic FRL signature in colon cancer and explore its potential molecular function.

Methods

RNA-seq data and relevant clinical information were obtained from The Cancer Genome Atlas (TCGA) database, and a list of ferroptosis-related genes was extracted from the FerrDb website. Analysis of differentially expressed FRLs was performed using the 'limma' package in R software. By implementing coexpression analysis and univariate Cox analysis, we then identified prognostic FRLs. Using Cox regression analysis with the least absolute shrinkage and selection operator (LASSO) algorithm, we constructed a prognostic model based on 4 FRLs. We evaluated the prognostic power of this model using Kaplan-Meier (K-M) survival curve analysis and receiver operating characteristic (ROC) curve analysis. Moreover, the relationships between the signature and immune landscape, somatic mutation and drug sensitivity were explored. Finally, in vitro experiments were conducted to validate the functions of AP003555.1 and AC000584.1.

Results

A 4-FRL signature was constructed. Two risk groups were classified based on the risk score calculated by this signature. The signature-based risk score exhibited a more powerful capacity for survival prediction than traditional clinicopathological features in colon patients. Additionally, we observed a significant difference in immune cells, such as CD4+ and CD8+ T cells and macrophages, between the two groups. Moreover, the high-risk group exhibited lower IC50 values for certain chemotherapy drugs, such as cisplatin, docetaxel, bleomycin or axitinib. Finally, the in vitro experiments showed that ferroptosis processes were suppressed after AP003555.1 and AC000584.1 knockdown.

Conclusion

The proposed 4-FRL signature is a promising biomarker to predict clinical outcomes and therapeutic responses in colon cancer patients.

Transcriptome Analysis of mRNAs and Long Non-Coding RNAs During Subsequent Embryo Development of Porcine Cloned Zygotes After Vitrification

Cryopreservation of porcine cloned zygotes has important implications for biotechnology and biomedicine research; however, lower embryo developmental potential remains an urgent problem to be resolved. For exploring the sublethal cryodamages during embryo development, this study was designed to acquire the mRNA and long non-coding RNA (lncRNA) profiles of 2-cells, 4-cells and blastocysts derived from vitrified porcine cloned zygotes using transcriptome sequencing. We identified 167 differentially expressed (DE) mRNAs and 516 DE lncRNAs in 2-cell stage, 469 DE mRNAs and 565 lncRNAs in 4-cell stage, and 389 DE mRNAs and 816 DE lncRNAs in blastocyst stage. Functional enrichment analysis revealed that the DE mRNAs during embryo development were involved in many regulatory mechanisms related to cell cycle, cell proliferation, apoptosis, metabolism and others. Moreover, the target genes of DE lncRNAs in the three embryonic stages were also enriched in many key GO terms or pathways such as “defense response”, “linoleic acid metabolic process”, “embryonic axis specification”, “negative regulation of protein neddylation”, etc., In conclusion, the present study provided comprehensive transcriptomic data about mRNAs and lncRNAs for the vitrified porcine cloned zygotes during different developmental stages, which contributed to further understand the potential cryodamage mechanisms responsible for impaired embryo development.

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.

Integrated lncRNA and mRNA Transcriptome Analyses in the Ovary of Cynoglossus semilaevis Reveal Genes and Pathways Potentially Involved in Reproduction

Long non-coding RNAs (lncRNAs) have been reported to be involved in multiple biological processes. However, the roles of lncRNAs in the reproduction of half-smooth tongue sole (Cynoglossus semilaevis) are unclear, especially in the molecular regulatory mechanism driving ovarian development and ovulation. Thus, to explore the mRNA and lncRNA mechanisms regulating reproduction, we collected tongue sole ovaries in three stages for RNA sequencing. In stage IV vs. V, we identified 312 differentially expressed (DE) mRNAs and 58 DE lncRNAs. In stage V vs. VI, we identified 1,059 DE mRNAs and 187 DE lncRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses showed that DE mRNAs were enriched in ECM-receptor interaction, oocyte meiosis and steroid hormone biosynthesis pathways. Furthermore, we carried out gene set enrichment analysis (GSEA) to identify potential reproduction related-pathways additionally, such as fatty metabolism and retinol metabolism. Based on enrichment analysis, DE mRNAs with a potential role in reproduction were selected and classified into six categories, including signal transduction, cell growth and death, immune response, metabolism, transport and catabolism, and cell junction. The interactions of DE lncRNAs and mRNAs were predicted according to antisense, cis-, and trans-regulatory mechanisms. We constructed a competing endogenous RNA (ceRNA) network. Several lncRNAs were predicted to regulate genes related to reproduction including cyp17a1, cyp19a1, mmp14, pgr, and hsd17b1. The functional enrichment analysis of these target genes of lncRNAs revealed that they were involved in several signaling pathways, such as the TGF-beta, Wnt signaling, and MAPK signaling pathways and reproduction related-pathways such as the progesterone-mediated oocyte maturation, oocyte meiosis, and GnRH signaling pathway. RT-qPCR analysis showed that two lncRNAs (XR_522278.2 and XR_522171.2) were mainly expressed in the ovary. Dual-fluorescence in situ hybridization experiments showed that both XR_522278.2 and XR_522171.2 colocalized with their target genes cyp17a1 and cyp19a1, respectively, in the follicular cell layer. The results further demonstrated that lncRNAs might be involved in the biological processes by modulating gene expression. Taken together, this study provides lncRNA profiles in the ovary of tongue sole and further insight into the role of lncRNA involvement in regulating reproduction in tongue sole.