Whole transcriptome sequencing reveals that non-coding RNAs are related to embryo morphogenesis and development in rabbits

Comprehensive investigation of differentially expressed ncRNAs, mRNAs, and their ceRNA networks in the regulation of shell color formation in clam, Cyclina sinensis

Noncoding RNAs (ncRNAs) have gained significant attention in recent years due to their crucial roles in various biological processes. However, our understanding of the expression and functions of ncRNAs in Cyclina sinensis, an economically important marine bivalve, remains limited. This study aimed to address this knowledge gap by systematically identifying ncRNAs in the mantles of C. sinensis with purple and white shells. Through our analysis, we identified a differential expression of 1244 mRNAs, 196 lncRNAs, 49 circRNAs, and 23 miRNAs between purple- and white-shell clams. Functional enrichment analysis revealed the involvement of these differentially expressed ncRNAs in biomineralization and pigmentation processes. To gain further insights into the regulatory mechanisms underlying shell color formation, we established competitive endogenous RNA (ceRNA) networks. These networks allowed us to identify targeted differentially expressed miRNAs (DEMis) and genes associated with shell color formation. Based on the ceRNA networks, we obtained an up-down-up lncRNA-miRNA-mRNA network consisting of 13 upregulated lncRNAs and a circRNA-miRNA-mRNA network comprising three upregulated circRNAs (novel_circ_0004787, novel_circ_0001165, novel_circ_0000245). Through these networks, we identified and selected an upregulated novel gene (evm.TU.Hic_asm_7.988) and a downregulated sponge miRNA (hru-miR-1985) as potential contributors to shell color regulation. In summary, the present investigation offers a comprehensive analysis of ncRNA catalogs expressed in the clam mantle of C. sinensis. The findings enhance our comprehension of the molecular mechanisms governing shell coloration and offer new perspectives for selective breeding of C. sinensis in the future.

Comprehensive transcriptional analysis of early dorsal skin development in pigs

Porcine skin is similar to human skin in physiology, anatomy and histology and is often used as a model animal for human skin research. There are few studies on the transcriptome aspects of pig skin during the embryonic period. In this study, RNA sequencing was performed on the dorsal skin of Chenghua sows at embryonic day 56 (E56), embryonic day 76 (E76), embryonic day 105 (E105), and 3 days after birth (D3) to explore RNA changes in pig dorsal skin at four ages. A number of skin-related differential genes were identified by intercomparison between RNAs at four time points, and KEGG functional analysis showed that these differential genes were mainly enriched in metabolic and developmental, immune, and disease pathways, and the pathways enriched in GO analysis were highly overlapping. Collagen is an important part of the skin, with type I collagen making up the largest portion. In this study, collagen type I alpha 1 (COL1A1) and collagen type I alpha 2 (COL1A2) were significantly upregulated at four time points. In addition, lncRNA-miRNA-mRNA and miRNA-circRNA coexpression networks were constructed. The data obtained may help to explain age-related changes in transcriptional patterns during skin development and provide further references for understanding human skin development at the molecular level.

Whole-Transcriptome Sequencing of Ovary Reveals the ceRNA Regulation Network in Egg Production of Gaoyou Duck

To investigate the regulatory mechanism of the competing endogenous RNAs (ceRNAs) on the egg performance of Gaoyou ducks, full transcriptome sequencing was performed to analyze the ovarian tissues in Gaoyou ducks. The ducks were categorized into high- and low-yield groups based on the individual in-cage egg production records and the hematoxylin-eosin (HE) staining results. The differentially expressed genes (DEGs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs) were further processed by GO (gene ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) analyses. In total, 72 DEmRNAs; 23 DElncRNAs; 4 DEcircRNAs; and 5 signaling pathways, including the ovarian steroidogenesis, PI3K-Akt, hedgehog, tryptophan metabolism, and oocyte meiosis signaling pathways, were significantly enriched. These results suggest that they could be associated with the Gaoyou duck's ovarian function and affect the total egg production or double-yolked egg production. Furthermore, a coregulation network based on the related candidate ceRNAs across the high- and low-yield egg production groups was constructed. Our findings provide new insights into the mechanisms underlying the molecular regulation of related circRNA/lncRNA-miRNA-mRNA in the egg production and double-yolked egg traits of Gaoyou ducks.

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.

Circular RNA as a Novel Regulator and Promising Biomarker in Polycystic Ovary Syndrome

Polycystic ovary syndrome (PCOS) is a prevalent metabolic and reproductive disorder that causes low fertility in females. Despite its detrimental effects on women's health, care for PCOS has been impeded by its undefined pathogenesis. Thus, there is an urgent need to explore novel biomarkers and therapeutic targets for the diagnosis and treatment of PCOS. Circular RNAs (circRNAs) are a class of noncoding RNAs with covalently closed cyclic structures, present in high abundance, and show development-stage specific expression patterns. Recent studies have demonstrated that circRNAs participate in PCOS progression by modulating various biological functions, including cell proliferation, apoptosis, and steroidogenesis. In addition, circRNAs are widely present in the follicular fluid of women with PCOS, indicating their potential as diagnostic biomarkers and therapeutic targets for PCOS. This review provides the current knowledge of circRNAs in PCOS, including their regulatory functions and molecular mechanisms, and explores their potential as diagnostic biomarkers and therapeutic targets.

Non-coding RNAs from seminal plasma extracellular vesicles and success of live birth among couples undergoing fertility treatment

Background: Infertility remains a global health problem with male-factor infertility accounting for around 50% of cases. Understanding the molecular markers for the male contribution of live birth success has been limited. Here, we evaluated the expression levels of seminal plasma extracellular vesicle (spEV) non-coding RNAs (ncRNAs) in men of couples in relation with those with and without a successful live birth after infertility treatment. Method: Sperm-free spEV small RNA profiles were generated from 91 semen samples collected from male participants of couples undergoing assisted reproductive technology (ART) treatment. Couples were classified into two groups based on successful live birth (yes, n = 28) and (no, n = 63). Mapping of reads to human transcriptomes followed the order: miRNA > tRNA > piRNA > rRNA> "other" RNA > circRNA > lncRNA. Differential expression analysis of biotype-specific normalized read counts between groups were assessed using EdgeR (FDR<0.05). Result: We found a total of 12 differentially expressed spEV ncRNAs which included 10 circRNAs and two piRNAs between the live birth groups. Most (n = 8) of the identified circRNAs were downregulated in the no live birth group and targeted genes related to ontology terms such as negative reproductive system and head development, tissue morphogenesis, embryo development ending in birth or egg hatching, and vesicle-mediated transport. The differentially upregulated piRNAs overlapped with genomic regions including coding PID1 genes previously known to play a role in mitochondrion morphogenesis, signal transduction and cellular proliferation. Conclusion: This study identified novel ncRNAs profiles of spEVs differentiating men of couples with and without live birth and emphasizes the role of the male partner for ART success.

Integrated analysis of microRNAs, circular RNAs, long non-coding RNAs, and mRNAs revealed competing endogenous RNA networks involved in brown adipose tissue whitening in rabbits

BACKGROUND

The brown adipose tissue (BAT) is a target for treating obesity. BAT losses thermogenic capacity and gains a "white adipose tissue-like" phenotype ("BAT whitening") under thermoneutral environments, which is a potential factor causing a low curative effect in BAT-related obesity treatments. Circular RNAs (circRNAs) and long non-coding RNAs (lncRNAs) can act as competing endogenous RNAs (ceRNA) to mRNAs and function in various processes by sponging shared microRNAs (miRNAs). However, the roles of circRNA- and lncRNA-related ceRNA networks in regulating BAT whitening remain litter known.

RESULTS

In this study, BATs were collected from rabbits at day0 (D0), D15, D85, and 2 years (Y2). MiRNA-seq was performed to investigate miRNA changes during BAT whitening. Then, a combined analysis of circRNA-seq and whole-transcriptome sequencing was used for circRNA assembly and quantification during BAT whitening. Our data showed that 1187 miRNAs and 6204 circRNAs were expressed in the samples, and many of which were identified as significantly changed during BAT whitening. Target prediction showed that D0-selective miRNAs were significantly enriched in the Ras, MAPK, and PI3K-Akt signaling pathways, and Y2-selective miRNAs were predicted to be involved in cell proliferation. The cyclization of several circRNAs could form novel response elements of key thermogenesis miRNAs at the back-splicing junction (BSJ) sites, and in combination with a dual-luciferase reporter assay confirmed the binding between the BSJ site of novel_circ_0013792 and ocu-miR-378-5p. CircRNAs and lncRNAs have high cooperativity in sponging miRNAs during BAT whitening. Both circRNA-miRNA-mRNA and lncRNA-miRNA-mRNA triple networks were significantly involved in immune response-associated biological processes. The D15-selective circRNA might promote BAT whitening by increasing the expression of IDH2. The Y2-selective circRNA-related ceRNA network and lncRNA-related ceRNA network might regulate the formation of the WAT-like phenotype of BAT via MAPK and Ras signaling pathways, respectively.

CONCLUSIONS

Our work systematically revealed ceRNA networks during BAT whitening in rabbits and might provide new insight into BAT-based obesity treatments.

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.

Identification and characterization of circRNAs in the liver of blunt snout bream (Megalobrama amblycephala) infected with Aeromonas hydrophila

Circular RNAs (circRNAs), a class of non-coding RNAs, play an important role in regulating various biological processes. In the present study, circRNAs from the Megalobrama amblycephala liver were identified at five different time points post Aeromonas hydrophila using RNA-seq technology. A total of 250 circRNAs were identified, of which 106 were differentially expressed (DE) in ten pairwise comparisons. GO and KEGG analyses showed that the parental genes of DE circRNAs were enriched in phagocytosis, complement and coagulation cascades, and Fc gamma R-mediated phagocytosis pathways. According to ceRNA hypothesis, the interaction network of circRNAs, miRNAs and mRNAs was constructed. Moreover, WGCNA was conducted, and five specific modules significantly related to bacterial infection were identified. All the above results reveal the important role of circRNAs in immune response, which enriches the information of circRNAs in teleost, and helps to understand the immune response mechanism of M. amblycephala to A. hydrophila.

Transcriptome Analysis of In Vitro Fertilization and Parthenogenesis Activation during Early Embryonic Development in Pigs

Parthenogenesis activation (PA), as an important artificial breeding method, can stably preserve the dominant genotype of a species. However, the delayed development of PA embryos is still overly severe and largely leads to pre-implantation failure in pigs. The mechanisms underlying the deficiencies of PA embryos have not been completely understood. For further understanding of the molecular mechanism behind PA embryo failure, we performed transcriptome analysis among pig oocytes (meiosis II, MII) and early embryos at three developmental stages (zygote, morula, and blastocyst) in vitro fertilization (IVF) and PA group. Totally, 11,110 differentially expressed genes (DEGs), 4694 differentially expressed lincRNAs (DELs) were identified, and most DEGs enriched the regulation of apoptotic processes. Through cis- and trans-manner functional prediction, we found that hub lincRNAs were mostly involved in abnormal parthenogenesis embryonic development. In addition, twenty DE imprinted genes showed that some paternally imprinted genes in IVF displayed higher expression than that in PA. Notably, we identified that three DELs of imprinted genes (MEST, PLAGL1, and DIRAS3) were up regulated in IVF, and there was no significant change in PA group. Disordered expression of key genes for embryonic development might play key roles in abnormal parthenogenesis embryonic development. Our study indicates that embryos derived from different production techniques have varied in vitro development to the blastocyst stage, and they also affect the transcription level of corresponding genes, such as imprinted genes. This work will help future research on these genes and molecular-assisted breeding for pig parthenotes.

Identification and Analysis of Long Non-coding RNAs in Leuciscus waleckii Adapted to Highly Alkaline Conditions

Leuciscus waleckii is a freshwater fish that is known to inhabit the Dali Nor Lake, Inner Mongolia, China. The water in this lake has an HCO3 -/CO3 2- concentration of 54 mM (pH 9.6) and a salinity of 0.6‰. The physiological mechanisms that allow this fish to tolerate these saline/alkaline conditions have yet to be elucidated. Transcriptional component analysis has shown that the expression levels of a large number of genes involved in the pathways responsible for osmo-ionoregulation and arachidonic acid metabolism pathway expression change significantly (p < 0.05) during the regulation of acid-base balance under high alkaline stress. In this study, we investigated the role of long non-coding RNAs (lncRNAs) during adaptation to high alkaline conditions. Fish were challenged to an NaHCO3-adjusted alkalinity of 0 mM, 30 mM (pH 9.44 ± 0.08), and 50 mM (pH 9.55 ± 0.06) for 20 days in the laboratory. Gill and kidney tissues were then collected for high-throughput sequencing assays. A total of 159 million clean reads were obtained by high-throughput sequencing, and 41,248 lncRNA transcripts were identified. Of these, the mean number of exons and the mean length of the lncRNA transcripts were 4.8 and 2,079 bp, respectively. Based on the analysis of differential lncRNA transcript expression, a total of 5,244 and 6,571 lncRNA transcripts were found to be differentially expressed in the gills and kidneys, respectively. Results derived from Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the coding genes were correlated with the lncRNA expression profiles. GO analysis showed that many lncRNAs were enriched in the following processes: "transporter activity," "response to stimulus," and "binding." KEGG analysis further revealed that metabolic pathways were significantly enriched. A random selection of 16 lncRNA transcripts was tested by RT-qPCR; these results were consistent with our sequencing results. We found that a large number of genes, with the same expression profiles as those with differentially expressed lncRNAs, were associated with the regulation of acid-base balance, ion transport, and the excretion of ammonia and nitrogen. Collectively, our data indicate that lncRNA-regulated gene expression plays an important role in the process of adaptation to high alkaline conditions in L. waleckii.

Circular RNAs: Novel potential regulators in embryogenesis, female infertility, and pregnancy-related diseases

Circular RNAs (circRNAs) are endogenous noncoding RNAs with unique cyclic structures. Although they were previously considered as nonfunctional transcription byproducts, numerous studies have demonstrated that circRNAs regulate gene transcription and expression via different mechanisms. Reproductive health influences the quality of life and affects offspring propagation in women. CircRNAs have been found to modify pregnancy-related diseases, gynecologic cancers, polycystic ovary syndrome, aging, gamete, and embryo development. It's promising for circRNAs to be the novel diagnostic and therapeutic targets for multiple reproductive diseases. With the widespread application of assisted reproduction technology (ART), it has been revealed that circRNA identification contributes to estimating the quality of gametes and embryos, reflecting the success rate of ART. CRISPR-Cas9 gene editing technology has enabled the discovery of new roles of circRNAs. So far, the roles of circRNAs in the reproductive system remain poorly defined. In this review, we describe the classification and functions of circRNAs in embryogenesis and the female reproductive system diseases, revealing potential roles of circRNAs physiologically and pathologically. In so-doing, we provide ideas for developing circRNA-based therapeutic treatment and clinical application of various female reproductive system diseases.

Differentially expressed circular RNAs and the competing endogenous RNA network associated with preeclampsia

INTRODUCTION

Circular RNAs (circRNAs) are non-coding RNAs that are implicated in preeclampsia (PE) pathogenesis; however, their expression and functions in PE remain unclear. In this study, we aimed to investigate the expression of circRNAs in PE and construct a competing endogenous RNA (ceRNA) network, and analyze the associated pathways in PE pathogenesis.

METHODS

We performed circRNA sequencing to identify the differential expression profile of circRNAs in PE as compared to normal pregnancy. The circRNA candidates were validated using quantitative reverse transcription polymerase chain reaction (qRT-PCR). Subsequently, we used datasets from the GEO database to generate the interaction network between circRNAs, microRNAs (miRNAs), and mRNAs. GO and KEGG enrichment analyses were performed to understand the functional significance of the differentially expressed circRNAs in PE.

RESULTS

We identified 361 differentially expressed circRNAs (252 upregulated and 109 downregulated) in preeclamptic placentas. Within the selected 31 circRNAs, 6 of them were verified by qRT-PCR. GO and KEGG analyses revealed the potential pathways affected by these circRNAs, e.g., T cell receptor signaling and MAP kinase pathways. A total of 134 miRNAs and 199 mRNAs were revealed to be differentially expressed in PE by analyzing datasets from the GEO database. The circRNA-miRNA-mRNA network comprised 206 circRNAs, 50 miRNAs, and 38 mRNAs. KEGG analysis of the 38 mRNAs included pathways involved in AMPK and PI3K-Akt signaling.

DISCUSSION

Our results reported the differential expression profile of circRNAs and the circRNA-miRNA-mRNA network in PE, which provides potential therapeutic targets for this disease.