Identification and functional analysis of long non-coding RNAs in mouse cleavage stage embryonic development based on single cell transcriptome data

Identification and functional prediction of lncRNAs associated with intramuscular lipid deposition in Guangling donkeys

Many studies have shown that long non-coding RNAs (lncRNAs) play key regulatory roles in various biological processes. However, the importance and molecular regulatory mechanisms of lncRNAs in donkey intramuscular fat deposition remain to be further investigated. In this study, we used published transcriptomic data from the longissimus dorsi muscle of Guangling donkeys to identify lncRNAs and obtained 196 novel lncRNAs. Compared with the coding genes, the novel lncRNAs and the known lncRNAs exhibited some typical features, such as shorter transcript length and smaller exons. A total of 272 coding genes and 52 lncRNAs were differentially expressed between the longissimus dorsi muscles of the low-fat and high-fat groups. The differentially expressed genes were found to be involved in various biological processes related to lipid metabolism. The potential target genes of differentially expressed lncRNAs were predicted by cis and trans. Functional analysis of lncRNA targets showed that some lncRNAs may act on potential target genes involved in lipid metabolism processes and regulate lipid deposition in the longissimus dorsi muscle. This study provides valuable information for further investigation of the molecular mechanisms of lipid deposition traits in donkeys, which may improve meat traits and facilitate the selection process of donkeys in future breeding.

Moringa oleifera leaf polysaccharide regulates fecal microbiota and colonic transcriptome in calves

This study investigated the effects of Moringa oleifera polysaccharide on growth performance indicators, serum biochemical indicators, immune organ indicators, colonic morphology, colonic microbiomics and colonic transcriptomics in newborn calves. 21 newborn calves were randomly divided into three groups of 7 calves per treatment group: control group (no Moringa oleifera polysaccharide addition); low-dose group (Moringa oleifera polysaccharide 0.5 g/kg); and high-dose group (Moringa oleifera polysaccharide 1 g/kg). This trial used gavage to feed MOP to calves. The test lasted 8 weeks. Calves were humanely electroshocked on the last day of the trial and slaughtered afterwards. Thymus, spleen, blood and colonic contents were collected for further testing. The results of this trial showed that MOP significantly increased the body weight of newborn calves and reduced the rate of calf diarrhea, thus promoting calf growth. Fecal scores showed a linear decrease with the addition of MOP. In terms of serum biochemistry, feeding MOP significantly increased serum ALB levels in a linear fashion. In terms of serum antioxidants, feeding MOP linearly increased CAT and T-AOC levels and decreased MDA concentrations, and in terms of serum immunity, feeding MOP linearly increased IgA, IgG, and IgM levels. At the same time, MOP regulated the abundance of Firmicutes and Bacteroidetes in the intestinal tract of calves, which reduced the occurrence of diarrhea. In addition, moringa polysaccharide could regulate genes related to inflammatory signaling pathways such as MAPK signaling pathway, TGF-beta signaling pathway, PI3K-Akt signaling pathway and TNF signaling pathway in calves' intestine to reduce the occurrence of intestinal inflammation. In conclusion, MOP can be used as a novel ruminant additive for the prevention of enteritis in calves.

Genome-Wide Identification and Involvement in Response to Biotic and Abiotic Stresses of lncRNAs in Turbot (Scophthalmus maximus)

Long non-coding RNAs (lncRNAs) play crucial roles in a variety of biological processes, including stress response. However, the number, characteristics and stress-related expression of lncRNAs in turbot are still largely unknown. In this study, a total of 12,999 lncRNAs were identified at the genome-wide level of turbot for the first time using 24 RNA-seq datasets. Sequence characteristic analyses of transcripts showed that lncRNA transcripts were shorter in average length, lower in average GC content and in average expression level as compared to the coding genes. Expression pattern analyses of lncRNAs in 12 distinct tissues showed that lncRNAs, especially lincRNA, exhibited stronger tissue-specific expression than coding genes. Moreover, 612, 1351, 1060, 875, 420 and 1689 differentially expressed (DE) lncRNAs under Vibrio anguillarum, Enteromyxum scophthalmi, and Megalocytivirus infection and heat, oxygen, and salinity stress conditions were identified, respectively. Among them, 151 and 62 lncRNAs showed differential expression under various abiotic and biotic stresses, respectively, and 11 lncRNAs differentially expressed under both abiotic and biotic stresses were selected as comprehensive stress-responsive lncRNA candidates. Furthermore, expression pattern analysis and qPCR validation both verified the comprehensive stress-responsive functions of these 11 lncRNAs. In addition, 497 significantly co-expressed target genes (correlation coefficient (R) > 0.7 and q-value < 0.05) for these 11 comprehensive stress-responsive lncRNA candidates were identified. Finally, GO and KEGG enrichment analyses indicated that these target genes were enriched mainly in molecular function, such as cytokine activity and active transmembrane transporter activity, in biological processes, such as response to stimulus and immune response, and in pathways, such as protein families: signaling and cellular processes, transporters and metabolism. These findings not only provide valuable reference resources for further research on the molecular basis and function of lncRNAs in turbot but also help to accelerate the progress of molecularly selective breeding of stress-resistant turbot strains or varieties.

Identification and Functional Analysis of Transcriptome Profiles, Long Non-Coding RNAs, Single-Nucleotide Polymorphisms, and Alternative Splicing from the Oocyte to the Preimplantation Stage of Sheep by Single-Cell RNA Sequencing

Numerous dynamic and complicated processes characterize development from the oocyte to the embryo. However, given the importance of functional transcriptome profiles, long non-coding RNAs, single-nucleotide polymorphisms, and alternative splicing during embryonic development, the effect that these features have on the blastomeres of 2-, 4-, 8-, 16-cell, and morula stages of development has not been studied. Here, we carried out experiments to identify and functionally analyze the transcriptome profiles, long non-coding RNAs, single-nucleotide polymorphisms (SNPs), and alternative splicing (AS) of cells from sheep from the oocyte to the blastocyst developmental stages. We found between the oocyte and zygote groups significantly down-regulated genes and the second-largest change in gene expression occurred between the 8- and 16-cell stages. We used various methods to construct a profile to characterize cellular and molecular features and systematically analyze the related GO and KEGG profile of cells of all stages from the oocyte to the blastocyst. This large-scale, single-cell atlas provides key cellular information and will likely assist clinical studies in improving preimplantation genetic diagnosis.

Genome-wide characterization of lncRNAs and mRNAs in muscles with differential intramuscular fat contents

Meat quality is one of the most important traits in pig production. Long non-coding RNAs (lncRNAs) have been involved in diverse biological processes such as muscle development through regulating gene expression. However, studies on lncRNAs lag behind and a comparatively small number of lncRNAs have been identified in pigs. Also, the effects of lncRNAs on meat quality remain to be characterized. Here, we analyzed lncRNAs in longissimus thoracis (LT) and semitendinosus (ST) muscles, being different in meat quality, with RNA-sequencing technology. A total of 500 differentially expressed lncRNAs (DELs) and 2,094 protein-coding genes (DEGs) were identified. Through KEGG analysis on DELs, we first made clear that fat deposition might be the main reason resulting in the differential phenotype of LT and ST, for which cGMP–PKG and VEGF signaling pathways were the most important ones. In total, forty-one key DELs and 50 DEGs involved in the differential fat deposition were then characterized. One of the key genes, cAMP-response element binding protein 1, was selected to confirm its role in porcine adipogenesis with molecular biology methods and found that it promotes the differentiation of porcine preadipocytes, consistent with its higher expression level and intramuscular fat contents in LT than that in ST muscle. Furthermore, through integrated analysis of DELs and DEGs, transcription factors important for differential fat deposition were characterized among which BCL6 has the most target DEGs while MEF2A was targeted by the most DELs. The results provide candidate genes crucial for meat quality, which will contribute to improving meat quality with molecular-breeding strategies.

Research Progress of Totipotent Stem Cells

Totipotent stem cells (TSCs), can develop into complete organisms, used in biological fields such as regenerative medicine, mammalian breeding, and conservation. However, cells from early-stage embryos cultured are hard to self-renew and maintain developmental totipotency, which becomes a key factor limiting the research of TSCs. Fortunately, a break-through in the study of induced pluripotent stem cells returning to their totipotent state has been made, resulting in the establishment of multiple TSCs and igniting a new wave of stem cell research. Furthermore, the blastocyst-like structures can be generated by the established TSCs, which lays a foundation for synthetic embryos in vitro. In this review, we summarize the totipotent stage of the early embryos, the establishment and cultivation of TSCs, and the developmental ability exploration of TSCs to promote further research of TSCs.

Integrative Analysis of lncRNA-mRNA Co-expression Provides Novel Insights Into the Regulation of Developmental Transitions in Female Varroa destructor

Varroa destructor is a major pathogenic driver of the Western honeybee colony losses globally. Understanding the developmental regulation of V. destructor is critical to develop effective control measures. Development is a complex biological process regulated by numerous genes and long non-coding RNAs (lncRNAs); however, the underlying regulation of lncRNAs in the development of V. destructor remains unknown. In this study, we analyzed the RNA sequencing (RNA-Seq) data derived from the four stages of female V. destructor in the reproductive phase (i.e., egg, protonymph, deutonymph, and adult). The identified differentially expressed mRNAs and lncRNAs exhibited a stage-specific pattern during developmental transitions. Further functional enrichment established that fat digestion and absorption, ATP-binding cassette (ABC) transporters, mitogen-activated protein kinase (MAPK) signaling pathway, and ubiquitin-proteasome pathway play key roles in the maturation of female V. destructor. Moreover, the lncRNAs and mRNAs of some pivotal genes were significantly upregulated at the deutonymph stage, such as cuticle protein 65/6.4/63/38 and mucin 5AC, suggesting that deutonymph is the key stage of metamorphosis development and pathogen resistance acquisition for female V. destructor. Our study provides novel insights into a foundational understanding of V. destructor biology.

Long Noncoding RNAs: An Overview

Recently an explosion in the discovery of long noncoding RNAs (lncRNAs) was obtained by high-throughput sequencing. Genome-wide transcriptome analyses, in conjugation with research for epigenetic modifications of chromatins, identified a novel type of non-protein coding transcripts longer than 200 nucleotides named lncRNAs . They are gradually emerging as functional and critical participants in many physiological processes. Here we gave an overview of the characteristics, biological functions, and working mechanism for this new class of noncoding RNA molecules.

Dynamics of Known Long Non-Coding RNAs during the Maternal-to-Zygotic Transition in Rabbit

The control of pre-implantation development in mammals undergoes a maternal-to-zygotic transition (MZT) after fertilization. The transition involves maternal clearance and zygotic genome activation remodeling the terminal differentiated gamete to confer totipotency. In the study, we first determined the profile of long non-coding RNAs (lncRNAs) of mature rabbit oocyte, 2-cell, 4-cell, 8-cell, and morula embryos using RNA-seq. A total of 2673 known rabbit lncRNAs were identified. The lncRNAs exhibited dynamic expression patterns during pre-implantation development. Moreover, 107 differentially expressed lncRNAs (DE lncRNAs) were detected between mature oocyte and 2-cell embryo, while 419 DE lncRNAs were detected between 8-cell embryo and morula, consistent with the occurrence of minor and major zygotic genome activation (ZGA) wave of rabbit pre-implanted embryo. This study then predicted the potential target genes of DE lncRNAs based on the trans-regulation mechanism of lncRNAs. The GO and KEGG analyses showed that lncRNAs with stage-specific expression patterns promoted embryo cleavage and synchronic development by regulating gene transcription and translation, intracellular metabolism and organelle organization, and intercellular signaling transduction. The correlation analysis between mRNAs and lncRNAs identified that lncRNAs ENSOCUG00000034943 and ENSOCUG00000036338 may play a vital role in the late-period pre-implantation development by regulating ILF2 gene. This study also found that the sequential degradation of maternal lncRNAs occurred through maternal and zygotic pathways. Furthermore, the function analysis of the late-degraded lncRNAs suggested that these lncRNAs may play a role in the mRNA degradation in embryos via mRNA surveillance pathway. Therefore, this work provides a global view of known lncRNAs in rabbit pre-implantation development and highlights the role of lncRNAs in embryogenesis regulation.

Long noncoding RNA RP11-909N17.2 presages a poor prognosis of non-small cell lung cancer

BACKGROUND

Long non-coding RNAs (lncRNAs) were detected extraordinarily expressed in various tumors and could combine with microRNAs (miRNAs) to play important role in tumor cells. This study is to explore the role of lncRNA RP11-909N17.2 in NSCLC and discuss in what way it functions in NSCLC.

METHODS

120 NSCLC patients were enlisted in this study. Expression levels of lncRNA RP11-909N17.2 and miR-767-3p were detected and the correlation between lncRNA RP11-909N17.2 expression and the clinical data characteristics was analyzed. Prognosis potential of lncRNA RP11-909N17.2 was inferred with Kaplan-Meier and multivariate Cox regression assays. Biological functions of NSCLC cells were accessed by cell counting Kit-8, transwell migration and invasion assay. Mechanism of RP11-909N17.2 action on NSCLC cells was investigated by luciferase activity assay with wide-type or mutation.

RESULTS

LncRNA RP11-909N17.2 has an ascendant expression while miR-767-3p has descended one in NSCLC tissue specimens and cells. Over-expression of lncRNA RP11-909N17.2 can shorten the overall survival period of NSCLC patients when compared with low expression. Knockdown of lncRNA RP11-909N17.2 suppressed biology function of NSCLC cell including proliferation, migration, and invasion.

CONCLUSION

LncRNA RP11-909N17.2 can be developed into a prognostic index for NSCLC. LncRNA RP11-909N17.2 plays a promoting role in NSCLC cells possibly by binding miR-767-3p as a sponge.

Building Pluripotency Identity in the Early Embryo and Derived Stem Cells

The fusion of two highly differentiated cells, an oocyte with a spermatozoon, gives rise to the zygote, a single totipotent cell, which has the capability to develop into a complete, fully functional organism. Then, as development proceeds, a series of programmed cell divisions occur whereby the arising cells progressively acquire their own cellular and molecular identity, and totipotency narrows until when pluripotency is achieved. The path towards pluripotency involves transcriptome modulation, remodeling of the chromatin epigenetic landscape to which external modulators contribute. Both human and mouse embryos are a source of different types of pluripotent stem cells whose characteristics can be captured and maintained in vitro. The main aim of this review is to address the cellular properties and the molecular signature of the emerging cells during mouse and human early development, highlighting similarities and differences between the two species and between the embryos and their cognate stem cells.

Detection of genomic regions associated malformations in newborn piglets: a machine-learning approach

Background

A significant proportion of perinatal losses in pigs occurs due to congenital malformations. The purpose of this study is the identification of genomic loci associated with fetal malformations in piglets.

Methods

The malformations were divided into two groups: associated with limb defects (piglet splay leg) and associated with other congenital anomalies found in newborn piglets. 148 Landrace and 170 Large White piglets were selected for the study. A genome-wide association study based on the gradient boosting machine algorithm was performed to identify markers associated with congenital anomalies and piglet splay leg.

Results

Forty-nine SNPs (23 SNPs in Landrace pigs and 26 SNPs in Large White) were associated with congenital anomalies, 22 of which were localized in genes. A total of 156 SNPs (28 SNPs in Landrace; 128 in Large White) were identified for piglet splay leg, of which 79 SNPs were localized in genes. We have demonstrated that the gradient boosting machine algorithm can identify SNPs and their combinations associated with significant selection indicators of studied malformations and productive characteristics.

Data availability

Genotyping and phenotyping data are available at http://www.compubioverne.group/data-and-software/.

A systematic evaluation of bioinformatics tools for identification of long noncoding RNAs

High-throughput RNA sequencing unveiled the complexity of transcriptome and significantly increased the records of long noncoding RNAs (lncRNAs), which were reported to participate in a variety of biological processes. Identification of lncRNAs is a key step in lncRNA analysis, and a bunch of bioinformatics tools have been developed for this purpose in recent years. While these tools allow us to identify lncRNA more efficiently and accurately, they may produce inconsistent results, making selection a confusing issue. We compared the performance of 41 analysis models based on 14 software packages and different data sets, including high-quality data and low-quality data from 33 species. In addition, computational efficiency, robustness, and joint prediction of the models were explored. As a practical guidance, key points for lncRNA identification under different situations were summarized. In this investigation, no one of these models could be superior to others under all test conditions. The performance of a model relied to a great extent on the source of transcripts and the quality of assemblies. As general references, FEELnc_all_cl, CPC, and CPAT_mouse work well in most species while COME, CNCI, and lncScore are good choices for model organisms. Since these tools are sensitive to different factors such as the species involved and the quality of assembly, researchers must carefully select the appropriate tool based on the actual data. Alternatively, our test suggests that joint prediction could behave better than any single model if proper models were chosen. All scripts/data used in this research can be accessed at http://bioinfo.ihb.ac.cn/elit.

Sperm Originated Chromatin Imprints and LincRNAs in Organismal Development and Cancer

Importance of sperm-derived transcripts and chromatin imprints in organismal development is poorly investigated. Here using an integrative approach, we show that human sperm transcripts are equally important as oocyte. Sperm-specific and sperm-oocyte common transcripts carry distinct chromatin structures at their promoters correlating with corresponding transcript levels in sperm. Interestingly, sperm-specific H3K4me3 patterns at the lincRNA promoters are not maintained in the germ layers and somatic tissues. However, bivalent chromatin at the sperm-specific protein-coding gene promoters is maintained throughout the development. Sperm-specific transcripts reach their peak expression during zygotic genome activation, whereas sperm-oocyte common transcripts are present during early preimplantation development but decline at the onset of zygotic genome activation. Additionally, there is an inverse correlation between sperm-specific and sperm-oocyte lincRNAs throughout the development. Sperm-lincRNAs also show aberrant activation in tumors. Overall, our observations indicate that sperm transcripts carrying chromatin imprints may play an important role in human development and cancer.

•

Sp-lincRNAs carry distinct chromatin structures correlating with transcript levels

•

Sp-lincRNAs are active during ZGA in preimplantation developmental stages

•

SpOc-lincRNAs are active in pre-ZGA and decline at the onset of ZGA

•

Sp-lincRNAs are silent post implantation but show aberrant cancer-specific activation

Long Noncoding RNA HOTTIP Serves as an Independent Predictive Biomarker for the Prognosis of Patients with Clear Cell Renal Cell Carcinoma

Several studies have indicated that HOXA transcript at the distal tip (HOTTIP) play important roles in the tumorigenesis and development of various cancers. We aim to investigate the expression and prognostic value of HOTTIP in clear cell renal cell carcinoma (ccRCC). A systematic review of PubMed, Embase, Medline, and Web of Science databases was performed to select eligible literatures relevant to the correlation between HOTTIP expression and clinical outcome of different cancers. The association between the HOTTIP level and overall survival (OS), lymph node metastasis (LNM), or clinical stage was subsequently analyzed. Survival analyses were performed in a large cohort of more than 500 patients with ccRCC from The Cancer Genome Atlas (TCGA) using bioinformatic methods. Seventeen studies with a total of 1594 patients with thirteen kinds of carcinomas were included in this analysis. The result showed that high HOTTIP expression could predict worse outcome in cancer patients, with the pooled hazard ratio (HR) of 2.34 (95% confidence interval (CI) 1.96–2.79, p < 0.0001). The result also showed that elevated HOTTIP expression was correlated with more LNM (OR = 2.61, 95% CI 1.91-3.58, p < 0.0001) and advanced clinical stage (OR = 3.57, 95% CI 2.58-4.93, p < 0.0001). We further validated that ccRCC patients with higher HOTTIP expression tend to have unsatisfactory outcomes both in the entire TCGA dataset and different clinical stratums, like age, grade, and stage. The tumor of those patients was associated with a larger size, easier to metastasis, advanced clinical stage, and a higher pathological grade. These findings suggested that increased HOTTIP expression might act as a novel prognostic marker for ccRCC patients.

Identification and Comparative Analysis of Long Non-Coding RNA in the Skeletal Muscle of Two Dezhou Donkey Strains

Long non-coding RNA (lncRNA) has been extensively studied in many livestock. However, compared with other livestock breeds, there is less research regarding donkey lncRNA function. It has been reported that lncRNA plays an important role in the timing control of development, aging, and death of livestock. Therefore, the study of donkey skeletal muscle lncRNA is of great significance for exploring donkey meat production performance. In this study, RNA-Seq was used to perform high-throughput sequencing of skeletal muscle (longissimus dorsi and gluteus) of two Dezhou donkey strains (SanFen and WuTou). The differentially expressed lncRNAs were screened between different strains and tissues. Then candidate genes for conjoint analysis were screened based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Finally, the accuracy of the sequencing data was verified by real-time quantitative polymerase chain reaction (RT-qPCR). Herein, we identified 3869 novel lncRNAs and 73 differentially expressed lncRNAs. Through the comparison between groups, the specific expression of lncRNAs were found in different strains and muscle tissues. Importantly, we constructed the lncRNA-miRNA-mRNA interaction network and found three important candidate lncRNAs (MSTRG.9787.1, MSTRG.3144.1, and MSTRG.9886.1) and four candidate genes (ACTN1, CDON, FMOD, and BMPR1B). Analysis of the KEGG pathway indicates that the TGF-β signaling pathway plays a pivotal role in the growth and development of skeletal muscle in Dezhou donkey strains.

Long noncoding RNA 2193 regulates meiosis through global epigenetic modification and cytoskeleton organization in pig oocytes

Long noncoding RNAs (lncRNAs) regulate a variety of physiological and pathological processes. However, the biological function of lncRNAs in mammalian germ cells remains largely unexplored. Here we identified one novel lncRNA (lncRNA2193) from single-cell RNA sequencing performed on porcine oocytes and investigated its function in oocyte meiosis. During in vitro maturation (IVM), from germinal vesicle (GV, 0 hr), GV breakdown (GVBD, 24 hr), to metaphase II stage (MII, 44 hr), the transcriptional abundance of lncRNA2193 remained stable and high. LncRNA2193 interference by small interfering RNA microinjection into porcine GV oocytes could significantly inhibit rates of GVBD and the first polar body extrusion, but enhance the rates of oocytes with a nuclear abnormality. Moreover, lncRNA2193 knockdown disturbed cytoskeletal organization (F-actin and spindle), and decreased DNA 5-methylcytosine (5mC) and histone trimethylation (H3K4me3, H3K9me3, H3K27me3, and H3K36me3) levels. The lncRNA2193 downregulation induced a decrease of 5mC level could be partially due to the reduction of DNA methyltransferase 3A and 3B, and the elevation of 5mC-hydroxylase ten-11 translocation 2 (TET2). After parthenogenetic activation of MII oocytes, parthenotes exhibited higher fragmentation but lower cleavage rates in the lncRNA2193 downregulated group. However, lncRNA2193 interference performed on mature MII oocytes and parthenotes at 1-cell stage did not affect the cleavage and blasctocyst rates of pathenotes. Taken together, lncRNA2193 plays an important role in porcine oocyte maturation, providing more insights for relevant investigations on mammalian germ cells.

Cell lineage-specific transcriptome analysis for interpreting cell fate specification of proembryos

In Arabidopsis, a zygote undergoes asymmetrical cell division that establishes the first two distinct cell types of early proembryos, apical and basal cells. However, the genome-wide transcriptional activities that guide divergence of apical and basal cell development remain unknown. Here, we present a comprehensive transcriptome analysis of apical and basal cell lineages, uncovering distinct molecular pathways during cell lineage specification. Selective deletion of inherited transcripts and specific de novo transcription contribute to the establishment of cell lineage-specific pathways for cell fate specification. Embryo-related pathways have been specifically activated in apical cell lineage since 1-cell embryo stage, but quick transcriptome remodeling toward suspensor-specific pathways are found in basal cell lineage. Furthermore, long noncoding RNAs and alternative splicing isoforms may be involved in cell lineage specification. This work also provides a valuable lineage-specific transcriptome resource to elucidate the molecular pathways for divergence of apical and basal cell lineages at genome-wide scale.

Asymmetric division of the Arabidopsis zygote produces apical and basal cells that mainly develop into embryo and suspensor, respectively. Here, Zhou et al. show that de novo transcription and selective RNA turnover establish distinct apical and basal transcriptomes as early as the 1-cell stage of embryo development.

Long noncoding RNA TUG1 regulates prostate cancer cell proliferation, invasion and migration via the Nrf2 signaling axis

BACKGROUND

Long noncoding RNAs (lncRNAs) have been identified to modulate the development and progression of prostate cancer (PCa) via the regulation of their target genes. However, the biological function underlying the effect of lncRNA TUG1 in PCa remains unclear.

METHODS

Reverse transcription-quantitative polymerase chain reaction (qRT-PCR) and Western blotting analysis were used to assess the mRNA expression of TUG1 and protein expression levels of Nrf2 pathway members, respectively. The migration, invasion, and proliferation abilities of cells were assessed by the wound-healing, Transwell migration/invasion, and CCK8 assays, respectively.

RESULTS

TUG1 was strikingly upregulated in PCa cells compared with non-tumorigenic human prostate epithelial cells. The LncTar Web Server, which is a bioinformatics tool, was used to predict the target association between TUG1 and Nrf2. Moreover, the expression of TUG1 showed a strikingly positive correlation with that of Nrf2 in TCGA PCa RNA-Seq data (r = 0.26,P = 4.63E-09). Subsequently, inhibition of TUG1 using siRNA resulted in deceased proliferation, migration, and invasion of PCa cells; however, these effects were reversed by treatment with oltipraz (an activator of Nrf2). Finally, we evaluated the Nrf2 pathway to reveal the underlying mechanism of TUG1 in PCa cells, and found that TUG1 knockdown decreased the protein expression of Nrf2 downstream members (e.g., HO-1, FTH1, and NQO1).

CONCLUSIONS

LncRNA TUG1 plays an oncogenic role in human PCa cells by promoting the cell proliferation and invasion in PCa cell lines, at least partly via the Nrf2 signaling pathway.

A Single Cell but Many Different Transcripts: A Journey into the World of Long Non-Coding RNAs

In late 2012 it was evidenced that most of the human genome is transcribed but only a small percentage of the transcripts are translated. This observation supported the importance of non-coding RNAs and it was confirmed in several organisms. The most abundant non-translated transcripts are long non-coding RNAs (lncRNAs). In contrast to protein-coding RNAs, they show a more cell-specific expression. To understand the function of lncRNAs, it is fundamental to investigate in which cells they are preferentially expressed and to detect their subcellular localization. Recent improvements of techniques that localize single RNA molecules in tissues like single-cell RNA sequencing and fluorescence amplification methods have given a considerable boost in the knowledge of the lncRNA functions. In recent years, single-cell transcription variability was associated with non-coding RNA expression, revealing this class of RNAs as important transcripts in the cell lineage specification. The purpose of this review is to collect updated information about lncRNA classification and new findings on their function derived from single-cell analysis. We also retained useful for all researchers to describe the methods available for single-cell analysis and the databases collecting single-cell and lncRNA data. Tables are included to schematize, describe, and compare exposed concepts.

Long noncoding RNA CCDC144NL-AS1 knockdown induces naïve-like state conversion of human pluripotent stem cells

BACKGROUND

Human naïve pluripotency state cells can be derived from direct isolation of inner cell mass or primed-to-naïve resetting of human embryonic stem cells (hESCs) through different combinations of transcription factors, small molecular inhibitors, and growth factors. Long noncoding RNAs (lncRNAs) have been identified to be crucial in diverse biological processes, including pluripotency regulatory circuit of mouse pluripotent stem cells (PSCs), but few are involved in human PSCs' regulation of pluripotency and naïve pluripotency derivation. This study initially planned to discover more lncRNAs possibly playing significant roles in the regulation of human PSCs' pluripotency, but accidently identified a lncRNA whose knockdown in human PSCs induced naïve-like pluripotency conversion.

METHODS

Candidate lncRNAs tightly correlated with human pluripotency were screened from 55 RNA-seq data containing human ESC, human induced pluripotent stem cell (iPSC), and somatic tissue samples. Then loss-of-function experiments in human PSCs were performed to investigate the function of these candidate lncRNAs. The naïve-like pluripotency conversion caused by CCDC144NL-AS1 knockdown (KD) was characterized by quantitative real-time PCR, immunofluorescence staining, western blotting, differentiation of hESCs in vitro and in vivo, RNA-seq, and chromatin immunoprecipitation. Finally, the signaling pathways in CCDC144NL-AS1-KD human PSCs were examined through western blotting and analysis of RNA-seq data.

RESULTS

The results indicated that knockdown of CCDC144NL-AS1 induces naïve-like state conversion of human PSCs in the absence of additional transcription factors or small molecular inhibitors. CCDC144NL-AS1-KD human PSCs reveal naïve-like pluripotency features, such as elevated expression of naïve pluripotency-associated genes, increased developmental capacity, analogous transcriptional profiles to human naïve PSCs, and global reduction of repressive chromatin modification marks. Furthermore, CCDC144NL-AS1-KD human PSCs display inhibition of MAPK (ERK), accumulation of active β-catenin, and upregulation of some LIF/STAT3 target genes, and all of these are concordant with previously reported traits of human naïve PSCs.

CONCLUSIONS

Our study unveils an unexpected role of a lncRNA, CCDC144NL-AS1, in the naïve-like state conversion of human PSCs, providing a new perspective to further understand the regulation process of human early pluripotency states conversion. It is suggested that CCDC144NL-AS1 can be potentially valuable for future research on deriving higher quality naïve state human PSCs and promoting their therapeutic applications.

The oocyte-to-embryo transition in mouse: past, present, and future

The oocyte-to-embryo transition (OET) arguably initiates with formation of a primordial follicle and culminates with reprogramming of gene expression during the course of zygotic genome activation. This transition results in converting a highly differentiated cell, i.e. oocyte, to undifferentiated cells, i.e. initial blastomeres of a preimplantation embryo. A plethora of changes occur during the OET and include, but are not limited to, changes in transcription, chromatin structure, and protein synthesis; accumulation of macromolecules and organelles that will comprise the oocyte's maternal contribution to the early embryo; sequential acquisition of meiotic and developmental competence to name but a few. This review will focus on transcriptional and post-transcriptional changes that occur during OET in mouse because such changes are likely the major driving force for OET. We often take a historical and personal perspective, and highlight how advances in experimental methods often catalyzed conceptual advances in understanding the molecular bases for OET. We also point out questions that remain open and therefore represent topics of interest for future investigation.

Non-coding RNA MFI2-AS1 promotes colorectal cancer cell proliferation, migration and invasion through miR-574-5p/MYCBP axis

Objective

Long non‐coding RNAs (lncRNAs) and microRNAs (miRNAs) play essential roles in the tumour progression. LncRNAs mostly act as competing endogenous RNAs (ceRNAs) by sponging miRNAs. This study aimed to study the association of a novel lncRNA MFI2‐AS1 with miR‐574‐5p/MYCBP axis in the development of colorectal cancer (CRC).

Methods

Ninety‐four CRC tissues and paired adjacent non‐tumour tissues were included in our study. The relative expression level of MFI2‐AS1 was detected, and its relationship with clinico‐pathological factors was analysed. Then, the CRC cells lines (LoVo and RKO) were transfected with MFI2‐AS1 siRNA, miR‐574‐5p mimics and inhibitors. Cell proliferation, migration, invasion, cell cycle distribution and DNA damage in response to different transfection conditions were examined. Dual‐luciferase reporter assay was performed to identify the target interactions between MFI2‐AS1 and miR‐574‐5p, miR‐574‐5p and MYCBP.

Results

LncRNA MFI2‐AS1 and MYCBP were up‐regulated in CRC tissues when compared with adjacent non‐tumour tissues. The expression levels of MFI2‐AS1 were significantly associated with tumour histological grade, lymph and distant metastasis, TNM stage and vascular invasion. Both MFI2‐AS1 siRNA and miR‐574‐5p mimics inhibited proliferation, migration and invasion in LoVo and RKO cells. The transfection of miR‐574‐5p inhibitor showed MFI2‐AS1 siRNA‐induced changes in CRC cells. Dual‐luciferase reporter assay revealed target interactions between MFI2‐AS1 and miR‐574‐5p, miR‐574‐5p and MYCBP.

Conclusions

These findings suggested that lncRNA MFI2‐AS1 and MYCBP have promoting effects in CRC tissues. LncRNA MFI2‐AS1 promoted CRC cell proliferation, migration and invasion through activating MYCBP and by sponging miR‐574‐5p.

Identification and Functional Prediction of Long Intergenic Non-coding RNAs Related to Subcutaneous Adipose Development in Pigs

An increasing number of studies have shown that long intergenic non-coding RNAs (lincRNAs) are a very important class of non-coding RNAs that plays a vital role in many biological processes. Adipose tissue is an important place for storing energy, but few studies on lincRNAs were related to pig subcutaneous fat development. Here, we used published RNA-seq data from subcutaneous adipose tissue of Italian Large White pigs and identified 252 putative lincRNAs, wherein 34 were unannotated. These lincRNAs had relatively shorter length, lower number of exons, and lower expression level compared with protein-coding transcripts. Gene ontology and pathway analysis indicated that the adjacent genes of lincRNAs were involved in lipid metabolism. In addition, differentially expressed lincRNAs (DELs) between low and high backfat thickness pigs were identified. Through the detection of quantitative trait locus (QTL), DELs were mainly located in QTLs related to adipose development. Based on the expression correlation of DEL genes and their differentially expressed potential target genes, we constructed a co-expression network and a potential pathway of DEL's effect on lipid metabolism. Our study identified and analyzed lincRNAs in subcutaneous adipose tissue, and results suggested that lincRNAs may be involved in the regulation of subcutaneous fat development. Our findings provided new insights into the biological function of porcine lincRNAs.

Long non-coding RNA repertoire and open chromatin regions constitute midbrain dopaminergic neuron - specific molecular signatures

Midbrain dopaminergic (DA) neurons are involved in diverse neurological functions, including control of movements, emotions or reward. In turn, their dysfunctions cause severe clinical manifestations in humans, such as the appearance of motor and cognitive symptoms in Parkinson’s Disease. The physiology and pathophysiology of these neurons are widely studied, mostly with respect to molecular mechanisms implicating protein-coding genes. In contrast, the contribution of non-coding elements of the genome to DA neuron function is poorly investigated. In this study, we isolated DA neurons from E14.5 ventral mesencephalons in mice, and used RNA-seq and ATAC-seq to establish and describe repertoires of long non-coding RNAs (lncRNAs) and putative DNA regulatory regions specific to this neuronal population. We identified 1,294 lncRNAs constituting the repertoire of DA neurons, among which 939 were novel. Most of them were not found in hindbrain serotonergic (5-HT) neurons, indicating a high degree of cell-specificity. This feature was also observed regarding open chromatin regions, as 39% of the ATAC-seq peaks from the DA repertoire were not detected in the 5-HT neurons. Our work provides for the first time DA-specific catalogues of non-coding elements of the genome that will undoubtedly participate in deepening our knowledge regarding DA neuronal development and dysfunctions.

Long noncoding RNAs of single hematopoietic stem and progenitor cells in healthy and dysplastic human bone marrow

Long noncoding RNAs (lncRNAs) are regulators of cell differentiation and development. The lncRNA transcriptome in human hematopoietic stem and progenitor cells is not comprehensively defined. We investigated lncRNAs in 979 human bone marrow-derived CD34⁺ cells by single cell RNA sequencing followed by de novo transcriptome reconstruction. We identified 3,173 lncRNAs in total, among which 2,365 were previously unknown, and we characterized lncRNA stem, differentiation, and maturation signatures. lncRNA expression exhibited high cell-to-cell variation, which was only apparent in single cell analysis. lncRNA expression followed a lineage-specific and highly dynamic pattern during early hematopoiesis. lncRNAs in hematopoietic cells closely correlated with protein-coding genes of known functions in the regulation of hematopoiesis and cell fate decisions, and the potential regulatory roles of lncRNAs in hematopoiesis were imputed by projection from protein-coding genes with a “guilt-by-association” approach. We characterized lncRNAs preferentially expressed in hematopoietic stem cells and in various downstream differentiated lineage progenitors. We also profiled lncRNA expression in single cells from patients with myelodysplastic syndromes and in aneuploid cells in particular. Our study provides a global view of lncRNAs in human hematopoietic stem and progenitor cells. We observed a highly ordered pattern of lncRNA expression and participation in regulation of early hematopoiesis, and coordinate aberrant messenger RNA and lncRNA transcriptomes in dysplastic hematopoiesis. (Registered at clinicaltrials.gov with identifiers: 00001620, 00001397)

Differential expression profiles of long non‑coding RNAs during the mouse pronuclear stage under normal gravity and simulated microgravity

Pronuclear migration, which is the initial stage of embryonic development and the marker of zygote formation, is a crucial process during mammalian preimplantation embryonic development. Recent studies have revealed that long non-coding RNAs (lncRNAs) serve an important role in early embryonic development. However, the functional regulation of lncRNAs in this process has yet to be elucidated, largely due to the difficulty of assessing gene expression alterations during the very short time in which pronuclear migration occurs. It has previously been reported that migration of the pronucleus of a zygote can be obstructed by simulated microgravity. To investigate pronuclear migration in mice, a rotary cell culture system was employed, which generates simulated microgravity, in order to interfere with murine pronuclear migration. Subsequently, lncRNA sequencing was performed to investigate the mechanism underlying this process. In the present study, a comprehensive analysis of lncRNA profile during the mouse pronuclear stage was conducted, in which 3,307 lncRNAs were identified based on single-cell RNA sequencing data. Furthermore, 52 lncRNAs were identified that were significantly differentially expressed. Subsequently, 10 lncRNAs were selected for validation by reverse transcription-quantitative polymerase chain reaction, in which the same relative expression pattern was observed. The results revealed that 12 lncRNAs (lnc006745, lnc007956, lnc013100, lnc013782, lnc017097, lnc019869, lnc025838, lnc027046, lnc005454, lnc007956, lnc019410 and lnc019607), with tubulin β 4B class IVb or actinin α 4 as target genes, may be associated with the expression of microtubule and microfilament proteins. Binding association was confirmed using a dual-luciferase reporter assay. Finally, Gene Ontology analysis revealed that the target genes of the differentially expressed lncRNAs participated in cellular processes associated with protein transport, binding, catalytic activity, membrane-bounded organelle, protein complex and the cortical cytoskeleton. These findings suggested that these lncRNAs may be associated with migration of the mouse pronucleus.

Upregulation of long non-coding RNA RAB1A-2 induces FGF1 expression worsening lung cancer prognosis

The chromosomal locations of lncRNAs (long non-coding RNAs, lncRNAs) infer their biological functions in cancer. Lnc-RAB1A-2, a Ras-related protein Rab-1A (RAB1A) upstream lncRNA, was chosen for assessment of its impact on lung cancer prognosis in a case-based analysis and investigation of its biological function though a series of functional assays. Lnc-RAB1A-2 was significantly upregulated in 276 lung cancer tissues compared with corresponding non-tumor tissues, and its expression level was significantly correlated with clinical stage and metastasis status in lung cancer patients. Patients with high expression levels of this lncRNA had a shorter median survival time (16.0 months vs. 23.0 months, P = 0.011 in southern samples; 8.0 months vs. 19.0 months, P = 0.020 in eastern samples; 13.0 months vs. 19.0 months, P = 0.002 in merged samples) and a higher risk of death than those with lower levels (HR = 1.52; 95% CI = 1.01-2.26, in merged samples). Additionally, overexpression of lnc-RAB1A-2 significantly promoted lung cancer cell proliferation in vitro and in vivo. Further analyses using digital gene expression tag profiling revealed that lnc-RAB1A-2 could affect the expression of fibroblast growth factor 1 (FGF1), a gene involved in the PI3K/AKT/mTOR pathway that is largely activated by RAB1A. FGF1 was confirmed to be a down-stream gene of lnc-RAB1A-2. Collectively, our study demonstrated that lnc-RAB1A-2 is associated with poor lung cancer prognosis by promoting lung cancer development.

Long non-coding RNAs potentially function synergistically in the cellular reprogramming of SCNT embryos

BACKGROUND

Long non-coding RNAs (lncRNAs), a type of epigenetic regulator, are thought to play important roles in embryonic development in mice, and several developmental defects are associated with epigenetic modification disorders. The most dramatic epigenetic reprogramming event occurs during somatic cell nuclear transfer (SCNT) when the expression profile of a differentiated cell is abolished, and a newly embryo-specific expression profile is established. However, the molecular mechanism underlying somatic reprogramming remains unclear, and the dynamics and functions of lncRNAs in this process have not yet been illustrated, resulting in inefficient reprogramming.

RESULTS

In this study, 63 single-cell RNA-seq libraries were first generated and sequenced. A total of 7009 mouse polyadenylation lncRNAs (including 5204 novel lncRNAs) were obtained, and a comprehensive analysis of in vivo and SCNT mouse pre-implantation embryo lncRNAs was further performed based on our single-cell RNA sequencing data. Expression profile analysis revealed that lncRNAs were expressed in a developmental stage-specific manner during mouse early-stage embryonic development, whereas a more temporal and spatially specific expression pattern was identified in mouse SCNT embryos with changes in the state of chromatin during somatic cell reprogramming, leading to incomplete zygotic genome activation, oocyte to embryo transition and 2-cell to 4-cell transition. No obvious differences between other stages and mouse NTC or NTM embryos at the same stage were observed. Gene oncology (GO) enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and weighted gene co-expression network analysis (WGCNA) of lncRNAs and their association with known protein-coding genes suggested that several lncRNAs and their associated with known protein-coding genes might be involved in mouse embryonic development and cell reprogramming.

CONCLUSIONS

This is a novel report on the expression landscapes of lncRNAs of mouse NT embryos by scRNA-seq analysis. This study will provide insight into the molecular mechanism underlying the involvement of lncRNAs in mouse pre-implantation embryonic development and epigenetic reprogramming in mammalian species after SCNT-based cloning.

Preliminary RNA-Seq Analysis of Long Non-Coding RNAs Expressed in Human Term Placenta

Development of particular structures and proper functioning of the placenta are under the influence of sophisticated pathways, controlled by the expression of substantial genes that are additionally regulated by long non-coding RNAs (lncRNAs). To date, the expression profile of lncRNA in human term placenta has not been fully established. This study was conducted to characterize the lncRNA expression profile in human term placenta and to verify whether there are differences in the transcriptomic profile between the sex of the fetus and pregnancy multiplicity. RNA-Seq data were used to profile, quantify, and classify lncRNAs in human term placenta. The applied methodology enabled detection of the expression of 4463 isoforms from 2899 annotated lncRNA loci, plus 990 putative lncRNA transcripts from 607 intergenic regions. Those placentally expressed lncRNAs displayed features such as shorter transcript length, longer exon length, fewer exons, and lower expression levels compared to messenger RNAs (mRNAs). Among all placental transcripts, 175,268 were classified as mRNAs and 15,819 as lncRNAs, and 56,727 variants were discovered within unannotated regions. Five differentially expressed lncRNAs (HAND2-AS1, XIST, RP1-97J1.2, AC010084.1, TTTY15) were identified by a sex-bias comparison. Splicing events were detected within 37 genes and 4 lncRNA loci. Functional analysis of cis-related potential targets for lncRNAs identified 2021 enriched genes. It is presumed that the obtained data will expand the current knowledge of lncRNAs in placenta and human non-coding catalogs, making them more contemporary and specific.

LncRNA TUG1 promoted viability and associated with gemcitabine resistant in pancreatic ductal adenocarcinoma

OBJECTIVE

To investigate the underlying mechanism of lncRNA TUG1 in pancreatic ductal adenocarcinoma (PDAC).

METHODS

The expression of TUG1 was defined by qRT-PCR. The apoptotic cells were detected by flow cytometry assay. The cell migration and invasion were measured by scratch assay and Transwell assay. The level of ERK pathway was detected using Western blot.

RESULTS

Compared with normal tissues and cells, the expression of TUG1 was up-regulated in pancreatic cancer tissue and cells. Meanwhile, knockdown of TUG1 could promote PDAC cells apoptosis and inhibit PDAC cells viability, migration and invasion. In addition, overexpression of TUG1 enhanced the gemcitabine chemoresistance of PDAC cells. Surprisingly, gemcitabine combined with SCH772984 (a suppressor of ERK pathway) could reverse the drug resistance resulted from overexpression of TUG1.

CONCLUSION

TUG1 promoted the viability of PDAC cells and enhanced its resistance of gemcitabine.

Long non-coding RNA exchange during the oocyte-to-embryo transition in mice

The oocyte-to-embryo transition (OET) transforms a differentiated gamete into pluripotent blastomeres. The accompanying maternal-zygotic RNA exchange involves remodeling of the long non-coding RNA (lncRNA) pool. Here, we used next generation sequencing and de novo transcript assembly to define the core population of 1,600 lncRNAs expressed during the OET (lncRNAs). Relative to mRNAs, OET lncRNAs were less expressed and had shorter transcripts, mainly due to fewer exons and shorter 5′ terminal exons. Approximately half of OET lncRNA promoters originated in retrotransposons suggesting their recent emergence. Except for a small group of ubiquitous lncRNAs, maternal and zygotic lncRNAs formed two distinct populations. The bulk of maternal lncRNAs was degraded before the zygotic genome activation. Interestingly, maternal lncRNAs seemed to undergo cytoplasmic polyadenylation observed for dormant mRNAs. We also identified lncRNAs giving rise to trans-acting short interfering RNAs, which represent a novel lncRNA category. Altogether, we defined the core OET lncRNA transcriptome and characterized its remodeling during early development. Our results are consistent with the notion that rapidly evolving lncRNAs constitute signatures of cells-of-origin while a minority plays an active role in control of gene expression across OET. Our data presented here provide an excellent source for further OET lncRNA studies.

Long non-coding RNAs involved in the regulatory network during porcine pre-implantation embryonic development and iPSC induction

Long non-coding RNAs (lncRNA) play a key role in the orchestration of transcriptional regulation during development and many other cellular processes. The importance of the regulatory co-expression network was highlighted in the identification of the mechanism of these processes in humans and mice. However, elucidation of the properties of porcine lncRNAs involved in the regulatory network during pre-implantation embryonic development and fibroblast reprogramming to induced pluripotent stem cell (iPSC) has been limited to date. Using a weighted gene co-expression network analysis, we constructed the regulatory network and determined that the novel lncRNAs were functionally involved in key events of embryonic development during the pre-implantation period; moreover, reprogramming could be delineated by a small number of potentially functional modules of co-expressed genes. These findings indicate that lncRNAs may be involved in the transcriptional regulation of zygotic genome activation, first lineage segregation and somatic reprogramming to pluripotency. Furthermore, we performed a conservation and synteny analysis with the significant lncRNAs involved in these vital events and validated the results via experimental assays. In summary, the current findings provide a valuable resource to dissect the protein coding gene and lncRNA regulatory networks that underlie the progressive development of embryos and somatic reprogramming.

Cross-species inference of long non-coding RNAs greatly expands the ruminant transcriptome

Background

mRNA-like long non-coding RNAs (lncRNAs) are a significant component of mammalian transcriptomes, although most are expressed only at low levels, with high tissue-specificity and/or at specific developmental stages. Thus, in many cases lncRNA detection by RNA-sequencing (RNA-seq) is compromised by stochastic sampling. To account for this and create a catalogue of ruminant lncRNAs, we compared de novo assembled lncRNAs derived from large RNA-seq datasets in transcriptional atlas projects for sheep and goats with previous lncRNAs assembled in cattle and human. We then combined the novel lncRNAs with the sheep transcriptional atlas to identify co-regulated sets of protein-coding and non-coding loci.

Results

Few lncRNAs could be reproducibly assembled from a single dataset, even with deep sequencing of the same tissues from multiple animals. Furthermore, there was little sequence overlap between lncRNAs that were assembled from pooled RNA-seq data. We combined positional conservation (synteny) with cross-species mapping of candidate lncRNAs to identify a consensus set of ruminant lncRNAs and then used the RNA-seq data to demonstrate detectable and reproducible expression in each species. In sheep, 20 to 30% of lncRNAs were located close to protein-coding genes with which they are strongly co-expressed, which is consistent with the evolutionary origin of some ncRNAs in enhancer sequences. Nevertheless, most of the lncRNAs are not co-expressed with neighbouring protein-coding genes.

Conclusions

Alongside substantially expanding the ruminant lncRNA repertoire, the outcomes of our analysis demonstrate that stochastic sampling can be partly overcome by combining RNA-seq datasets from related species. This has practical implications for the future discovery of lncRNAs in other species.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0391-0) contains supplementary material, which is available to authorized users.

Identification and Functional Analysis of Long Intergenic Non-coding RNAs Underlying Intramuscular Fat Content in Pigs

Intramuscular fat (IMF) content is an important trait that can affect pork quality. Previous studies have identified many genes that can regulate IMF. Long intergenic non-coding RNAs (lincRNAs) are emerging as key regulators in various biological processes. However, lincRNAs related to IMF in pig are largely unknown, and the mechanisms by which they regulate IMF are yet to be elucidated. Here we reconstructed 105,687 transcripts and identified 1,032 lincRNAs in pig longissimus dorsi muscle (LDM) of four stages with different IMF contents based on published RNA-seq. These lincRNAs show typical characteristics such as shorter length and lower expression compared with protein-coding genes. Combined with methylation data, we found that both the promoter and genebody methylation of lincRNAs can negatively regulate lincRNA expression. We found that lincRNAs exhibit high correlation with their protein-coding neighbors in expression. Co-expression network analysis resulted in eight stage-specific modules, gene ontology and pathway analysis of them suggested that some lincRNAs were involved in IMF-related processes, such as fatty acid metabolism and peroxisome proliferator-activated receptor signaling pathway. Furthermore, we identified hub lincRNAs and found six of them may play important roles in IMF development. This work detailed some lincRNAs which may affect of IMF development in pig, and facilitated future research on these lincRNAs and molecular assisted breeding for pig.

Identification and functional analysis of long non-coding RNAs in human and mouse early embryos based on single-cell transcriptome data

Epigenetics regulations have an important role in fertilization and proper embryonic development, and several human diseases are associated with epigenetic modification disorders, such as Rett syndrome, Beckwith-Wiedemann syndrome and Angelman syndrome. However, the dynamics and functions of long non-coding RNAs (lncRNAs), one type of epigenetic regulators, in human pre-implantation development have not yet been demonstrated. In this study, a comprehensive analysis of human and mouse early-stage embryonic lncRNAs was performed based on public single-cell RNA sequencing data. Expression profile analysis revealed that lncRNAs are expressed in a developmental stage-specific manner during human early-stage embryonic development, whereas a more temporal-specific expression pattern was identified in mouse embryos. Weighted gene co-expression network analysis suggested that lncRNAs involved in human early-stage embryonic development are associated with several important functions and processes, such as oocyte maturation, zygotic genome activation and mitochondrial functions. We also found that the network of lncRNAs involved in zygotic genome activation was highly preservative between human and mouse embryos, whereas in other stages no strong correlation between human and mouse embryo was observed. This study provides insight into the molecular mechanism underlying lncRNA involvement in human pre-implantation embryonic development.

Long non-coding RNA TCONS_00041960 enhances osteogenesis and inhibits adipogenesis of rat bone marrow mesenchymal stem cell by targeting miR-204-5p and miR-125a-3p

A growing number of long non‐coding RNAs (lncRNAs) have been found to be involved in diverse biological processes such as cell cycle regulation, embryonic development, and cell differentiation. However, limited knowledge is available concerning the underlying mechanisms of lncRNA functions. In this study, we found down‐regulation of TCONS_00041960 during adipogenic and osteogenic differentiation of glucocorticoid‐treated bone marrow mesenchymal stem cells (BMSCs). Furthermore, up‐regulation of TCONS_00041960 promoted expression of osteogenic genes Runx2, osterix, and osteocalcin, and anti‐adipogenic gene glucocorticoid‐induced leucine zipper (GILZ). Conversely, expression of adipocyte‐specific markers was decreased in the presence of over‐expressed TCONS_00041960. Mechanistically, we determined that TCONS_00041960 as a competing endogenous RNA interacted with miR‐204‐5p and miR‐125a‐3p to regulate Runx2 and GILZ, respectively. Overall, we identified a new TCONS_00041960‐miR‐204‐5p/miR‐125a‐3p‐Runx2/GILZ axis involved in regulation of adipogenic and osteogenic differentiation of glucocorticoid‐treated BMSCs.

Understanding the Role of lncRNAs in Nervous System Development

The diversity of lncRNAs has expanded within mammals in tandem with the evolution of increased brain complexity, suggesting that lncRNAs play an integral role in this process. In this chapter, we will highlight the identification and characterization of lncRNAs in nervous system development. We discuss the potential role of lncRNAs in nervous system and brain evolution, along with efforts to create comprehensive catalogues that analyze spatial and temporal changes in lncRNA expression during nervous system development. Additionally, we focus on recent endeavors that attempt to assign function to lncRNAs during nervous system development. We highlight discrepancies that have been observed between in vitro and in vivo studies of lncRNA function and the challenges facing researchers in conducting mechanistic analyses of lncRNAs in the developing nervous system. Altogether, this chapter highlights the emerging role of lncRNAs in the developing brain and sheds light on novel, RNA-mediated mechanisms by which nervous system development is controlled.

Noncoding RNAs: Master Regulators of Inflammatory Signaling

Inflammatory signaling underlies many diseases, from arthritis to cancer. Our understanding of inflammation has thus far been limited to the world of proteins, because we are only just beginning to understand the role that noncoding RNAs (ncRNA) might play. It is now clear that ncRNA do not constitute transcriptional 'noise' but instead harbor physiological functions in controlling signaling pathways. In this review, we cover the newly discovered mechanisms and functions of ncRNAs in the regulation of inflammatory signaling. We also describe advances in experimental techniques allowing this field of research to take root. These findings have opened new avenues for putative therapeutic intervention in inflammatory diseases, which may be seen translated into clinical outcomes in the future.

Long and small noncoding RNAs during oocyte-to-embryo transition in mammals

Oocyte-to-embryo transition is a process during which an oocyte ovulates, is fertilized, and becomes a developing embryo. It involves the first major genome reprogramming event in life of an organism where gene expression, which gave rise to a differentiated oocyte, is remodeled in order to establish totipotency in blastomeres of an early embryo. This remodeling involves replacement of maternal RNAs with zygotic RNAs through maternal RNA degradation and zygotic genome activation. This review is focused on expression and function of long noncoding RNAs (lncRNAs) and small RNAs during oocyte-to-embryo transition in mammals. LncRNAs are an assorted rapidly evolving collection of RNAs, which have no apparent protein-coding capacity. Their biogenesis is similar to mRNAs including transcriptional control and post-transcriptional processing. Diverse molecular and biological roles were assigned to lncRNAs although most of them probably did not acquire a detectable biological role. Since some lncRNAs serve as precursors for small noncoding regulatory RNAs in RNA silencing pathways, both types of noncoding RNA are reviewed together.

Genome-wide analysis and functional prediction of long non-coding RNAs in mouse uterus during the implantation window

Establishment of the receptive uterus is a crucial step for embryo implantation. In this study, the expression profiles and characterization of long non-coding RNAs (lncRNAs) in pregnant mouse uteri on day 4, day 5 at implantation sites and inter-implantation sites were conducted using RNA-seq. A total of 7,764 putative lncRNA transcripts were identified, including 6,179 known lncRNA transcripts and 1,585 novel lncRNA transcripts. Bioinformatics analysis of the cis and trans lncRNA targets showed that the differentially expressed lncRNAs were mainly involved in tissue remodelling, immune response and metabolism-related processes, indicating that lncRNAs could be involved in the regulation of embryo implantation. We also discovered that differentially expressed lncRNAs might regulate multiple signalling pathways that play an important role in the regulation of embryo implantation. In addition, nine known lncRNAs and four novel lncRNAs were randomly selected and validated by qRT-PCR. The expression of Tug1, Neat1, Gas5, Malat1, H19 and Rmst were significantly regulated in the mouse uterus during the implantation window. Our results are the first to systematically identify lncRNAs in the mouse uterus and provide a catalogue of lncRNAs for further understanding their functions in pregnant mouse uteri during the implantation window.

Single cell transcriptomics of noncoding RNAs and their cell-specificity

Recent developments of single cell transcriptome profiling methods have led to the realization that many seemingly homogeneous cells have surprising levels of expression variability. The biological implications of the high degree of variability is unclear but one possibility is that many genes are restricted in expression to small lineages of cells, suggesting the existence of many more cell types than previously estimated. Noncoding RNA (ncRNA) are thought to be key parts of gene regulatory processes and their single cell expression patterns may help to dissect the biological function of single cell variability. Technology for measuring ncRNA in single cell is still in development and most of the current single cell datasets have reliable measurements for only long noncoding RNA (lncRNA). Most works report that lncRNAs show lineage-specific restricted expression patterns, which suggest that they might determine, at least in part, lineage fates and cell subtypes. However, evidence is still inconclusive as to whether lncRNAs and other ncRNAs are more lineage-specific than protein-coding genes. Nevertheless, measurement of ncRNAs in single cells will be important for studies of cell types and single cell function. WIREs RNA 2017, 8:e1433. doi: 10.1002/wrna.1433 For further resources related to this article, please visit the WIREs website.

Single-cell transcriptomics reveals gene signatures and alterations associated with aging in distinct neural stem/progenitor cell subpopulations

Aging associated cognitive decline has been linked to dampened neural stem/progenitor cells (NSC/NPCs) activities manifested by decreased proliferation, reduced propensity to produce neurons, and increased differentiation into astrocytes. While gene transcription changes objectively reveal molecular alterations of cells undergoing various biological processes, the search for molecular mechanisms underlying aging of NSC/NPCs has been confronted by the enormous heterogeneity in cellular compositions of the brain and the complex cellular microenvironment where NSC/NPCs reside. Moreover, brain NSC/NPCs themselves are not a homogenous population, making it even more difficult to uncover NSC/NPC sub-type specific aging mechanisms. Here, using both population-based and single cell transcriptome analyses of young and aged mouse forebrain ependymal and subependymal regions and comprehensive “big-data” processing, we report that NSC/NPCs reside in a rather inflammatory environment in aged brain, which likely contributes to the differentiation bias towards astrocytes versus neurons. Moreover, single cell transcriptome analyses revealed that different aged NSC/NPC subpopulations, while all have reduced cell proliferation, use different gene transcription programs to regulate age-dependent decline in cell cycle. Interestingly, changes in cell proliferation capacity are not influenced by inflammatory cytokines, but likely result from cell intrinsic mechanisms. The Erk/Mapk pathway appears to be critically involved in regulating age-dependent changes in the capacity for NSC/NPCs to undergo clonal expansion. Together this study is the first example of using population and single cell based transcriptome analyses to unveil the molecular interplay between different NSC/NPCs and their microenvironment in the context of the aging brain.

Identification of long non-coding RNA in the horse transcriptome

Background

Efforts to resolve the transcribed sequences in the equine genome have focused on protein-coding RNA. The transcription of the intergenic regions, although detected via total RNA sequencing (RNA-seq), has yet to be characterized in the horse. The most recent equine transcriptome based on RNA-seq from several tissues was a prime opportunity to obtain a concurrent long non-coding RNA (lncRNA) database.

Results

This lncRNA database has a breadth of eight tissues and a depth of over 20 million reads for select tissues, providing the deepest and most expansive equine lncRNA database. Utilizing the intergenic reads and three categories of novel genes from a previously published equine transcriptome pipeline, we better describe these groups by annotating the lncRNA candidates. These lncRNA candidates were filtered using an approach adapted from human lncRNA annotation, which removes transcripts based on size, expression, protein-coding capability and distance to the start or stop of annotated protein-coding transcripts.

Conclusion

Our equine lncRNA database has 20,800 transcripts that demonstrate characteristics unique to lncRNA including low expression, low exon diversity and low levels of sequence conservation. These candidate lncRNA will serve as a baseline lncRNA annotation and begin to describe the RNA-seq reads assigned to the intergenic space in the horse.

Electronic supplementary material

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

Identification of long non-coding transcripts with feature selection: a comparative study

Background

The unveiling of long non-coding RNAs as important gene regulators in many biological contexts has increased the demand for efficient and robust computational methods to identify novel long non-coding RNAs from transcripts assembled with high throughput RNA-seq data. Several classes of sequence-based features have been proposed to distinguish between coding and non-coding transcripts. Among them, open reading frame, conservation scores, nucleotide arrangements, and RNA secondary structure have been used with success in literature to recognize intergenic long non-coding RNAs, a particular subclass of non-coding RNAs.

Results

In this paper we perform a systematic assessment of a wide collection of features extracted from sequence data. We use most of the features proposed in the literature, and we include, as a novel set of features, the occurrence of repeats contained in transposable elements. The aim is to detect signatures (groups of features) able to distinguish long non-coding transcripts from other classes, both protein-coding and non-coding. We evaluate different feature selection algorithms, test for signature stability, and evaluate the prediction ability of a signature with a machine learning algorithm. The study reveals different signatures in human, mouse, and zebrafish, highlighting that some features are shared among species, while others tend to be species-specific. Compared to coding potential tools and similar supervised approaches, including novel signatures, such as those identified here, in a machine learning algorithm improves the prediction performance, in terms of area under precision and recall curve, by 1 to 24%, depending on the species and on the signature.

Conclusions

Understanding which features are best suited for the prediction of long non-coding RNAs allows for the development of more effective automatic annotation pipelines especially relevant for poorly annotated genomes, such as zebrafish. We provide a web tool that recognizes novel long non-coding RNAs with the obtained signatures from fasta and gtf formats. The tool is available at the following url: http://www.bioinformatics-sannio.org/software/.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1594-z) contains supplementary material, which is available to authorized users.

Localisation of RNAs and proteins in nucleolar precursor bodies of early mouse embryos

Early embryos of all mammalian species contain morphologically distinct but transcriptionally silent nucleoli called the nucleolar precursor bodies (NPBs), which, unlike normal nucleoli, have been poorly studied at the biochemical level. To bridge this gap, here we examined the occurrence of RNA and proteins in early mouse embryos with two fluorochromes - an RNA-binding dye pyronin Y (PY) and the protein-binding dye fluorescein-5'-isothiocyanate (FITC). The staining patterns of zygotic NPBs were then compared with those of nucleolus-like bodies (NLBs) in fully grown surrounded nucleolus (SN)-type oocytes, which are morphologically similar to NPBs. We show that both entities contain proteins, but unlike NLBs, NPBs are significantly impoverished for RNA. Detectable amounts of RNA appear on the NPB surface only after resumption of rDNA transcription and includes pre-rRNAs and 28S rRNA as evidenced by fluorescence in situ hybridisation with specific oligonucleotide probes. Immunocytochemical assays demonstrate that zygotic NPBs contain rRNA processing factors fibrillarin, nucleophosmin and nucleolin, while UBF (the RNA polymerase I transcription factor) and ribosomal proteins RPL26 and RPS10 are not detectable. Based on the results obtained and data in the contemporary literature, we suggest a scheme of NPB assembly and maturation to normal nucleoli that assumes utilisation of maternally derived nucleolar proteins but of nascent rRNAs.

Comprehensive analysis of long non-coding RNAs highlights their spatio-temporal expression patterns and evolutional conservation in Sus scrofa

Despite modest sequence conservation and rapid evolution, long non-coding RNAs (lncRNAs) appear to be conserved in expression pattern and function. However, analysis of lncRNAs across tissues and developmental stages remains largely uncharacterized in mammals. Here, we systematically investigated the lncRNAs of the Guizhou miniature pig (Sus scrofa), which was widely used as biomedical model. We performed RNA sequencing across 9 organs and 3 developmental skeletal muscle, and developed a filtering pipeline to identify 10,813 lncRNAs (9,075 novel). Conservation patterns analysis revealed that 57% of pig lncRNAs showed homology to humans and mice based on genome alignment. 5,455 lncRNAs exhibited typical hallmarks of regulatory molecules, such as high spatio-temporal specificity. Notably, conserved lncRNAs exhibited higher tissue specificity than pig-specific lncRNAs and were significantly enriched in testis and ovary. Weighted co-expression network analysis revealed a set of conserved lncRNAs that are likely involved in postnatal muscle development. Based on the high degree of similarity in the structure, organization, and dynamic expression of pig lncRNAs compared with human and mouse lncRNAs, we propose that these lncRNAs play an important role in organ physiology and development in mammals. Our results provide a resource for studying animal evolution, morphological complexity, breeding, and biomedical research.