Transposable elements reveal a stem cell-specific class of long noncoding RNAs

Human lncRNAs harbor conserved modules embedded in different sequence contexts

We analyzed the structure of human long non-coding RNA (lncRNAs) genes to investigate whether the non-coding transcriptome is organized in modular domains, as is the case for protein-coding genes. To this aim, we compared all known human lncRNA exons and identified 340 pairs of exons with high sequence and/or secondary structure similarity but embedded in a dissimilar sequence context. We grouped these pairs in 106 clusters based on their reciprocal similarities. These shared modules are highly conserved between humans and the four great ape species, display evidence of purifying selection and likely arose as a result of recent segmental duplications. Our analysis contributes to the understanding of the mechanisms driving the evolution of the non-coding genome and suggests additional strategies towards deciphering the functional complexity of this class of molecules.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

piR-823 tale as emerging cancer-hallmark molecular marker in different cancer types: a step-toward ncRNA-precision

PIWI-interacting RNAs (piRNAs) have received a lot of attention for their functions in cancer research. This class of short non-coding RNAs (ncRNA) has roles in genomic stability, chromatin remodeling, messenger RNA (mRNA) integrity, and genome structure. We summarized the mechanisms underlying the biogenesis and regulatory molecular functions of piRNAs. Among all piRNAs studied in cancer, this review offers a comprehensive analysis of the emerging roles of piR-823 in various types of cancer, including colorectal, gastric, liver, breast, and renal cancers, as well as multiple myeloma. piR-823 has emerged as a crucial modulator of various cancer hallmarks through regulating multiple pathways. In the current review, we analyzed several databases and conducted an extensive literature search to explore the influence of piR-823 in carcinogenesis in addition to describing the potential application of piR-823 as prognostic and diagnostic markers as well as the therapeutic potential toward ncRNA precision.

Abnormal activation of genomic LINE1 elements caused by DNA demethylation contributes to lncRNA CASC9 overexpression in esophageal squamous cell carcinoma

Long noncoding RNA (lncRNA) cancer susceptibility 9 (CASC9) has been found to be overexpressed and functions as an oncogene in many cancer types. We investigated the molecular mechanism underlying CASC9 overexpression in esophageal squamous cell carcinoma (ESCC). Transcripts containing exons 2 and 6 and exons 4 and 6 showed the highest CASC9 expression levels in ESCC, no transcripts were detected in the normal esophageal epithelial Het1A cell line. The Long Interspersed Nuclear Element-1 (LINE1 or L1) element in the genome was found to participate in the evolution of lncRNA CASC9, the antisense promoter (ASP) of L1 provides the cis-regulatory elements necessary for CASC9 activation, and the antisense chain of L1 participates in the formation of exons of CASC9. The activation of the antisense promoter was due to the aberrant hypomethylation of L1 elements. An active enhancer element was identified in the downstream region of CASC9 gene by ChIP-seq and ChIP-qPCR. The interaction between ASP and the enhancer elements was confirmed by chromosome conformation capture (3C). Thus, our results suggest that the L1 ASP activation due to aberrant hypomethylation and downstream enhancer interaction plays a key role in the overexpression of lncRNA CASC9 in ESCC.

Regulatory roles of long non-coding RNAs in short-term heat stress in adult worker bees

Long non-coding RNAs (lncRNAs) are crucial modulators of post-transcriptional gene expression regulation, cell fate determination, and disease development. However, lncRNA functions during short-term heat stress in adult worker bees are poorly understood. Here, we performed deep sequencing and bioinformatic analyses of honeybee lncRNAs. RNA interference was performed by using siRNA targeting the most highly expressed lncRNA. The silencing effect on lncRNA and the relative expression levels of seven heat shock protein (HSP) genes, were subsequently examined. Overall, 7,842 lncRNAs and 115 differentially expressed lncRNAs (DELs) were identified in adult worker bees following heat stress exposure. Structural analysis revealed that the overall expression abundance, length of transcripts, exon number, and open reading frames of lncRNAs were lower than those of mRNAs. GO analysis revealed that the target genes were mainly involved in "metabolism," "protein folding," "response to stress," and "signal transduction" pathways. KEGG analysis indicated that the "protein processing in endoplasmic reticulum" and "longevity regulating pathway-multiple species" pathways were most enriched. Quantitative real-time polymerase chain reaction (qRT-PCR) detection of the selected DELs confirmed the reliability of the sequencing data. Moreover, the siRNA experiment indicated that feeding siRNA yielded a silencing efficiency of 77.51% for lncRNA MSTRG.9645.5. Upon silencing this lncRNA, the expression levels of three HSP genes were significantly downregulated (p < 0.05), whereas those of three other HSP genes were significantly upregulated (p < 0.05). Our results provide a new perspective for understanding the regulatory mechanisms of lncRNAs in adult worker bees under short-term heat stress.

Whole transcriptome mapping reveals the lncRNA regulatory network of TFP5 treatment in diabetic nephropathy

BACKGROUND

TFP5 is a Cdk5 inhibitor peptide, which could restore insulin production. However, the role of TFP5 in diabetic nephropathy (DN) is still unclear.

OBJECTIVE

This study aims to characterize the transcriptome profiles of mRNA and lncRNA in TFP5-treated DN mice to mine key lncRNAs associated with TFP5 efficacy.

METHODS

We evaluated the role of TFP5 in DN pathology and performed RNA sequencing in C57BL/6J control mice, C57BL/6J db/db model mice, and TFP5 treatment C57BL/6J db/db model mice. The differentially expressed lncRNAs (DElncRNAs) and mRNAs (DEmRNAs) were analyzed. WGCNA was used to screen hub-gene of TFP5 in treatment of DN.

RESULTS

Our results showed that TFP5 therapy ameliorated renal tubular injury in DN mice. In addition, compared with the control group, the expression profile of lncRNAs in the model group was significantly disordered, while TFP5 alleviated the abnormal expression of lncRNAs. A total of 67 DElncRNAs shared among the three groups, 39 DElncRNAs showed a trend of increasing in the DN group and decreasing after TFP treatment, while the remaining 28 showed the opposite trend. DElncRNAs were enriched in glycosphingolipid biosynthesis signaling pathways, NF-κB signaling pathways, and complement activation signaling pathways. There were 1028 up-regulated and 1117 down-regulated DEmRNAs in the model group compared to control group, and 123 up-regulated and 153 down-regulated DEmRNAs in the TFP5 group compared to the model group. The DEmRNAs were involved in PPAR and MAPK signaling pathway. We confirmed that MSTRG.28304.1 is a key DElncRNA for TFP5 treatment of DN. TFP5 ameliorated DN maybe by inhibiting MSTRG.28304.1 through regulating the insulin resistance and PPAR signaling pathway. The qRT-PCR results confirmed the reliability of the sequencing data through verifying the expression of ENSMUST00000211209, MSTRG.31814.5, MSTRG.28304.1, and MSTRG.45642.14.

CONCLUSION

Overall, the present study provides novel insights into molecular mechanisms of TFP5 treatment in DN.

Cell-specific housekeeping role of lncRNAs in COVID-19-infected and recovered patients

A plethora of studies have demonstrated the roles of lncRNAs in modulating disease severity and outcomes during infection. However, the spatio-temporal expression of these lncRNAs is poorly understood. In this study, we used single-cell RNA-seq to understand the spatio-temporal expression dynamics of lncRNAs across healthy, SARS-CoV-2-infected, and recovered individuals and their functional role in modulating the disease and recovery. We identified 203 differentially expressed lncRNAs, including cell type-specific ones like MALAT1, NEAT1, ZFAS1, SNHG7, SNHG8, and SNHG25 modulating immune function in classical monocyte, NK T, proliferating NK, plasmablast, naive, and activated B/T cells. Interestingly, we found invariant lncRNAs (no significant change in expression across conditions) regulating essential housekeeping functions (for example, HOTAIR, NRAV, SNHG27, SNHG28, and UCA1) in infected and recovered individuals. Despite similar repeat element abundance, variant lncRNAs displayed higher Alu content, suggesting increased interactions with proximal and distal genes, crucial for immune response modulation. The comparable repeat abundance but distinct expression levels of variant and invariant lncRNAs highlight the significance of investigating the regulatory mechanisms of invariant lncRNAs. Overall, this study offers new insights into the spatio-temporal expression patterns and functional roles of lncRNAs in SARS-CoV-2-infected and recovered individuals while highlighting the importance of invariant lncRNAs in the disease context.

The Killer's Web: Interconnection between Inflammation, Epigenetics and Nutrition in Cancer

Inflammation is a key contributor to both the initiation and progression of tumors, and it can be triggered by genetic instability within tumors, as well as by lifestyle and dietary factors. The inflammatory response plays a critical role in the genetic and epigenetic reprogramming of tumor cells, as well as in the cells that comprise the tumor microenvironment. Cells in the microenvironment acquire a phenotype that promotes immune evasion, progression, and metastasis. We will review the mechanisms and pathways involved in the interaction between tumors, inflammation, and nutrition, the limitations of current therapies, and discuss potential future therapeutic approaches.

Navigating the brain and aging: exploring the impact of transposable elements from health to disease

Transposable elements (TEs) are mobile genetic elements that constitute on average 45% of mammalian genomes. Their presence and activity in genomes represent a major source of genetic variability. While this is an important driver of genome evolution, TEs can also have deleterious effects on their hosts. A growing number of studies have focused on the role of TEs in the brain, both in physiological and pathological contexts. In the brain, their activity is believed to be important for neuronal plasticity. In neurological and age-related disorders, aberrant activity of TEs may contribute to disease etiology, although this remains unclear. After providing a comprehensive overview of transposable elements and their interactions with the host, this review summarizes the current understanding of TE activity within the brain, during the aging process, and in the context of neurological and age-related conditions.

Retrotransposon-derived transcripts and their functions in immunity and disease

Retrotransposons, which account for approximately 42% of the human genome, have been increasingly recognized as "non-self" pathogen-associated molecular patterns (PAMPs) due to their virus-like sequences. In abnormal conditions such as cancer and viral infections, retrotransposons that are aberrantly expressed due to impaired epigenetic suppression display PAMPs, leading to their recognition by pattern recognition receptors (PRRs) of the innate immune system and triggering inflammation. This viral mimicry mechanism has been observed in various human diseases, including aging and autoimmune disorders. However, recent evidence suggests that retrotransposons possess highly regulated immune reactivity and play important roles in the development and function of the immune system. In this review, I discuss a wide range of retrotransposon-derived transcripts, their role as targets in immune recognition, and the diseases associated with retrotransposon activity. Furthermore, I explore the implications of chimeric transcripts formed between retrotransposons and known gene mRNAs, which have been previously underestimated, for the increase of immune-related gene isoforms and their influence on immune function. Retrotransposon-derived transcripts have profound and multifaceted effects on immune system function. The aim of this comprehensive review is to provide a better understanding of the complex relationship between retrotransposon transcripts and immune defense.

Transposable elements acquire time- and sex-specific transcriptional and epigenetic signatures along mouse fetal gonad development

Gonadal sex determination in mice is a complex and dynamic process, which is crucial for the development of functional reproductive organs. The expression of genes involved in this process is regulated by a variety of genetic and epigenetic mechanisms. Recently, there has been increasing evidence that transposable elements (TEs), which are a class of mobile genetic elements, play a significant role in regulating gene expression during embryogenesis and organ development. In this study, we aimed to investigate the involvement of TEs in the regulation of gene expression during mouse embryonic gonadal development. Through bioinformatics analysis, we aimed to identify and characterize specific TEs that operate as regulatory elements for sex-specific genes, as well as their potential mechanisms of regulation. We identified TE loci expressed in a time- and sex-specific manner along fetal gonad development that correlate positively and negatively with nearby gene expression, suggesting that their expression is integrated to the gonadal regulatory network. Moreover, chromatin accessibility and histone post-transcriptional modification analyses in differentiating supporting cells revealed that TEs are acquiring a sex-specific signature for promoter-, enhancer-, and silencer-like elements, with some of them being proximal to critical sex-determining genes. Altogether, our study introduces TEs as the new potential players in the gene regulatory network that controls gonadal development in mammals.

SINE-Associated LncRNA SAWPA Regulates Porcine Zygotic Genome Activation

In mice, retrotransposon-associated long noncoding RNAs (lncRNA) play important regulatory roles in pre-implantation development; however, it is largely unknown whether they function in the pre-implantation development in pigs. The current study aims to screen for retrotransposon-associated lncRNA in porcine early embryos and identifies a porcine 8-cell embryo-specific SINE-associated nuclear long noncoding RNA named SAWPA. SAWPA is essential for porcine embryonic development as depletion of SAWPA results in a developmental arrest at the 8-cell stage, accompanied by the inhibition of the JNK-MAPK signaling pathway. Mechanistically, SAWPA works in trans as a transcription factor for JNK through the formation of an RNA-protein complex with HNRNPA1 and MED8 binding the SINE elements upstream of JNK. Therefore, as the first functional SINE-associated long noncoding RNAs in pigs, SAWPA provides novel insights for the mechanism research on retrotransposons in mammalian pre-implantation development.

Restrictor synergizes with Symplekin and PNUTS to terminate extragenic transcription

Transcription termination pathways mitigate the detrimental consequences of unscheduled promiscuous initiation occurring at hundreds of thousands of genomic cis-regulatory elements. The Restrictor complex, composed of the Pol II-interacting protein WDR82 and the RNA-binding protein ZC3H4, suppresses processive transcription at thousands of extragenic sites in mammalian genomes. Restrictor-driven termination does not involve nascent RNA cleavage, and its interplay with other termination machineries is unclear. Here we show that efficient termination at Restrictor-controlled extragenic transcription units involves the recruitment of the protein phosphatase 1 (PP1) regulatory subunit PNUTS, a negative regulator of the SPT5 elongation factor, and Symplekin, a protein associated with RNA cleavage complexes but also involved in cleavage-independent and phosphatase-dependent termination of noncoding RNAs in yeast. PNUTS and Symplekin act synergistically with, but independently from, Restrictor to dampen processive extragenic transcription. Moreover, the presence of limiting nuclear levels of Symplekin imposes a competition for its recruitment among multiple transcription termination machineries, resulting in mutual regulatory interactions. Hence, by synergizing with Restrictor, Symplekin and PNUTS enable efficient termination of processive, long-range extragenic transcription.

Conserved structured domains in plant non-coding RNA enod40, their evolution and recruitment of sequences from transposable elements

Plant long noncoding RNA enod40 is involved in the regulation of symbiotic associations with bacteria, in particular, in nitrogen-fixing root nodules of legumes, and with fungi in phosphate-acquiring arbuscular mycorrhizae formed by various plants. The presence of enod40 genes in plants that do not form such symbioses indicates its other roles in cell physiology. The molecular mechanisms of enod40 RNA function are poorly understood. Enod40 RNAs form several structured domains, conserved to different extents. Due to relatively low sequence similarity, identification of enod40 sequences in plant genomes is not straightforward, and many enod40 genes remain unannotated even in complete genomes. Here, we used comparative structure analysis and sequence similarity searches in order to locate enod40 genes and determine enod40 RNA structures in nitrogen-fixing clade plants and in grasses. The structures combine conserved features with considerable diversity of structural elements, including insertions of structured domain modules originating from transposable elements. Remarkably, these insertions contain sequences similar to tandem repeats and several stem-loops are homologous to microRNA precursors.

The human genome contains over a million autonomous exons

Mammalian mRNA and lncRNA exons are often small compared to introns. The exon definition model predicts that exons splice autonomously, dependent on proximal exon sequence features, explaining their delineation within large introns. This model has not been examined on a genome-wide scale, however, leaving open the question of how often mRNA and lncRNA exons are autonomous. It is also unknown how frequently such exons can arise by chance. Here, we directly assayed large fragments (500-1000 bp) of the human genome by exon trapping, which detects exons spliced into a heterologous transgene, here designed with a large intron context. We define the trapped exons as "autonomous." We obtained ∼1.25 million trapped exons, including most known mRNA and well-annotated lncRNA internal exons, demonstrating that human exons are predominantly autonomous. mRNA exons are trapped with the highest efficiency. Nearly a million of the trapped exons are unannotated, most located in intergenic regions and antisense to mRNA, with depletion from the forward strand of introns. These exons are not conserved, suggesting they are nonfunctional and arose from random mutations. They are nonetheless highly enriched with known splicing promoting sequence features that delineate known exons. Novel autonomous exons are more numerous than annotated lncRNA exons, and computational models also indicate they will occur with similar frequency in any randomly generated sequence. These results show that most human coding exons splice autonomously, and provide an explanation for the existence of many unconserved lncRNAs, as well as a new annotation and inclusion levels of spliceable loci in the human genome.

Deep annotation of long noncoding RNAs by assembling RNA-seq and small RNA-seq data

Long noncoding RNAs (lncRNAs) are increasingly being recognized as modulators in various biological processes. However, due to their low expression, their systematic characterization is difficult to determine. Here, we performed transcript annotation by a newly developed computational pipeline, termed RNA-seq and small RNA-seq combined strategy (RSCS), in a wide variety of cellular contexts. Thousands of high-confidence potential novel transcripts were identified by the RSCS, and the reliability of the transcriptome was verified by analysis of transcript structure, base composition, and sequence complexity. Evidenced by the length comparison, the frequency of the core promoter and the polyadenylation signal motifs, and the locations of transcription start and end sites, the transcripts appear to be full length. Furthermore, taking advantage of our strategy, we identified a large number of endogenous retrovirus-associated lncRNAs, and a novel endogenous retrovirus-lncRNA that was functionally involved in control of Yap1 expression and essential for early embryogenesis was identified. In summary, the RSCS can generate a more complete and precise transcriptome, and our findings greatly expanded the transcriptome annotation for the mammalian community.

Transposable elements as essential elements in the control of gene expression

Interspersed repetitions called transposable elements (TEs), commonly referred to as mobile elements, make up a significant portion of the genomes of higher animals. TEs contribute in controlling the expression of genes locally and even far away at the transcriptional and post-transcriptional levels, which is one of their significant functional effects on gene function and genome evolution. There are different mechanisms through which TEs control the expression of genes. First, TEs offer cis-regulatory regions in the genome with their inherent regulatory features for their own expression, making them potential factors for controlling the expression of the host genes. Promoter and enhancer elements contain cis-regulatory sites generated from TE, which function as binding sites for a variety of trans-acting factors. Second, a significant portion of miRNAs and long non-coding RNAs (lncRNAs) have been shown to have TEs that encode for regulatory RNAs, revealing the TE origin of these RNAs. Furthermore, it was shown that TE sequences are essential for these RNAs' regulatory actions, which include binding to the target mRNA. By being a member of cis-regulatory and regulatory RNA sequences, TEs therefore play essential regulatory roles. Additionally, it has been suggested that TE-derived regulatory RNAs and cis-regulatory regions both contribute to the evolutionary novelty of gene regulation. Additionally, these regulatory systems arising from TE frequently have tissue-specific functions. The objective of this review is to discuss TE-mediated gene regulation, with a particular emphasis on the processes, contributions of various TE types, differential roles of various tissue types, based mostly on recent studies on humans.

Context-aware transcript quantification from long-read RNA-seq data with Bambu

Most approaches to transcript quantification rely on fixed reference annotations; however, the transcriptome is dynamic and depending on the context, such static annotations contain inactive isoforms for some genes, whereas they are incomplete for others. Here we present Bambu, a method that performs machine-learning-based transcript discovery to enable quantification specific to the context of interest using long-read RNA-sequencing. To identify novel transcripts, Bambu estimates the novel discovery rate, which replaces arbitrary per-sample thresholds with a single, interpretable, precision-calibrated parameter. Bambu retains the full-length and unique read counts, enabling accurate quantification in presence of inactive isoforms. Compared to existing methods for transcript discovery, Bambu achieves greater precision without sacrificing sensitivity. We show that context-aware annotations improve quantification for both novel and known transcripts. We apply Bambu to quantify isoforms from repetitive HERVH-LTR7 retrotransposons in human embryonic stem cells, demonstrating the ability for context-specific transcript expression analysis.

Comprehensive transcriptomic analysis of three varieties with different brown planthopper-resistance identifies leaf sheath lncRNAs in rice

BACKGROUND

Long non-coding RNAs (lncRNAs) have been brought great attention for their crucial roles in diverse biological processes. However, systematic identification of lncRNAs associated with specialized rice pest, brown planthopper (BPH), defense in rice remains unexplored.

RESULTS

In this study, a genome-wide high throughput sequencing analysis was performed using leaf sheaths of susceptible rice Taichung Native 1 (TN1) and resistant rice IR36 and R476 with and without BPH feeding. A total of 2283 lncRNAs were identified, of which 649 lncRNAs were differentially expressed. During BPH infestation, 84 (120 in total), 52 (70 in total) and 63 (94 in total) of differentially expressed lncRNAs were found only in TN1, IR36 and R476, respectively. Through analyzing their cis-, trans-, and target mimic-activities, not only the lncRNAs targeting resistance genes (NBS-LRR and RLKs) and transcription factors, but also the lncRNAs acting as the targets of the well-studied stress-related miRNAs (miR2118, miR528, and miR1320) in each variety were identified. Before the BPH feeding, 238 and 312 lncRNAs were found to be differentially expressed in TN1 vs. IR36 and TN1 vs. R476, respectively. Among their putative targets, the plant-pathogen interaction pathway was significantly enriched. It is speculated that the resistant rice was in a priming state by the regulation of lncRNAs. Furthermore, the lncRNAs extensively involved in response to BPH feeding were identified by Weighted Gene Co-expression Network Analysis (WGCNA), and the possible regulation networks of the key lncRNAs were constructed. These lncRNAs regulate different pathways that contribute to the basal defense and specific resistance of rice to the BPH.

CONCLUSION

In summary, we identified the specific lncRNAs targeting the well-studied stress-related miRNAs, resistance genes, and transcription factors in each variety during BPH infestation. Additionally, the possible regulating network of the lncRNAs extensively responding to BPH feeding revealed by WGCNA were constructed. These findings will provide further understanding of the regulatory roles of lncRNAs in BPH defense, and lay a foundation for functional research on the candidate lncRNAs.

The Talented LncRNAs: Meshing into Transcriptional Regulatory Networks in Cancer

As a group of diseases characterized by uncontrollable cell growth, cancer is highly multifaceted in how it overrides checkpoints controlling proliferation. Amongst the regulators of these checkpoints, long non-coding RNAs (lncRNAs) can have key roles in why natural biological processes go haywire. LncRNAs represent a large class of regulatory transcripts that can localize anywhere in cells. They were found to affect gene expression on many levels from transcription to mRNA translation and even protein stability. LncRNA participation in such control mechanisms can depend on cell context, with given transcripts sometimes acting as oncogenes or tumor suppressors. Importantly, the tissue-specificity and low expression levels of lncRNAs make them attractive therapeutic targets or biomarkers. Here, we review the various cellular processes affected by lncRNAs and outline molecular strategies they use to control gene expression, particularly in cancer and in relation to transcription factors.

The transposable element-derived transcript of LIN28B has a placental origin and is not specific to tumours

Transposable elements (TEs) are genetic elements that have evolved as crucial regulators of human development and cancer, functioning as both genes and regulatory elements. When TEs become dysregulated in cancer cells, they can serve as alternate promoters to activate oncogenes, a process known as onco-exaptation. This study aimed to explore the expression and epigenetic regulation of onco-exaptation events in early human developmental tissues. We discovered co-expression of some TEs and oncogenes in human embryonic stem cells and first trimester and term placental tissues. Previous studies identified onco-exaptation events in various cancer types, including an AluJb SINE element-LIN28B interaction in lung cancer cells, and showed that the TE-derived LIN28B transcript is associated with poor patient prognosis in hepatocellular carcinoma. This study further characterized the AluJb-LIN28B transcript and confirmed that its expression is restricted to the placenta. Targeted DNA methylation analysis revealed differential methylation of the two LIN28B promoters between placenta and healthy somatic tissues, indicating that some TE-oncogene interactions are not cancer-specific but arise from the epigenetic reactivation of developmental TE-derived regulatory events. In conclusion, our findings provide evidence that some TE-oncogene interactions are not limited to cancer and may originate from the epigenetic reactivation of TE-derived regulatory events that are involved in early development. These insights broaden our understanding of the role of TEs in gene regulation and suggest the potential importance of targeting TEs in cancer therapy beyond their conventional use as cancer-specific markers.

Long non-coding RNAs: definitions, functions, challenges and recommendations

Genes specifying long non-coding RNAs (lncRNAs) occupy a large fraction of the genomes of complex organisms. The term 'lncRNAs' encompasses RNA polymerase I (Pol I), Pol II and Pol III transcribed RNAs, and RNAs from processed introns. The various functions of lncRNAs and their many isoforms and interleaved relationships with other genes make lncRNA classification and annotation difficult. Most lncRNAs evolve more rapidly than protein-coding sequences, are cell type specific and regulate many aspects of cell differentiation and development and other physiological processes. Many lncRNAs associate with chromatin-modifying complexes, are transcribed from enhancers and nucleate phase separation of nuclear condensates and domains, indicating an intimate link between lncRNA expression and the spatial control of gene expression during development. lncRNAs also have important roles in the cytoplasm and beyond, including in the regulation of translation, metabolism and signalling. lncRNAs often have a modular structure and are rich in repeats, which are increasingly being shown to be relevant to their function. In this Consensus Statement, we address the definition and nomenclature of lncRNAs and their conservation, expression, phenotypic visibility, structure and functions. We also discuss research challenges and provide recommendations to advance the understanding of the roles of lncRNAs in development, cell biology and disease.

Genome-wide identification of endogenous retrovirus elements and their active transcription in mink genome

Mammalian endogenous retroviruses (ERVs) are ancient retroviruses that have been integrated into genomes. ERVs were believed to be inactive until the discovery of ERV transcription in the mouse genome. However, the transcription level and function of ERV elements in mammalian genomes are not well understood. In this study, we performed the first genome-wide scanning of ERV loci in the American mink (Neogale vison) genome (NeoERV) followed by transcriptomic analysis to detect actively transcribed NeoERV elements. A total of 365,791 NeoERV loci were identified, and161,205 (44%) of these loci were found to be actively transcribed based on transcriptomic data from three types of tissues (amygdala, trachea and lung). More than one third of the actively transcribed NeoERV loci were tissue-specific. Furthermore, some of the active loci were associated with host gene transcription, and the level of NeoERV transcription was positively correlated with that of host genes, specifically when active loci were located in overlapped gene regions. An in-depth analysis of the envelope protein coding env gene showed that, in general, its transcription level was higher than that of NeoERVs, which is believed to be associated with host immunity.

The porcine piRNA transcriptome response to Senecavirus a infection

Introduction

Senecavirus A (SVA) belongs to the genus Senecavirus in the family Picornaviridae. PIWI-interacting RNAs (piRNAs) are a class of small Ribonucleic Acids (RNAs) that have been found in mammalian cells in recent years. However, the expression profile of piRNAs in the host during SVA infection and their roles are poorly understood.

Methods

Here, we found the significant differential expression of 173 piRNAs in SVA-infected porcine kidney (PK-15) cells using RNA-seq and 10 significant differentially expressed (DE) piRNAs were further verified by qRT-PCR.

Results

GO annotation analysis showed that metabolism, proliferation, and differentiation were significantly activated after SVA infection. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that significant DE piRNAs were mainly enriched in AMPK pathway, Rap1 pathway, circadian rhythm and VEGF pathway. It was suggested that piRNAs may regulated antiviral immunity, intracellular homeostasis, and tumor activities during SVA infection. In addition, we found that the expression levels of the major piRNA-generating genes BMAL1 and CRY1 were significantly downregulated after SVA infection.

Discussion

This suggests that SVA may affect circadian rhythm and promote apoptosis by inhibiting the major piRNA-generating genes BMAL1 and CRY1. The piRNA transcriptome in PK-15 cells has never been reported before, and this study will further the understanding of the piRNA regulatory mechanisms underlying SVA infections.

Repetitive elements in aging and neurodegeneration

Repetitive elements (REs), such as transposable elements (TEs) and satellites, comprise much of the genome. Here, we review how TEs and (peri)centromeric satellite DNA may contribute to aging and neurodegenerative disorders, including amyotrophic lateral sclerosis (ALS). Alterations in RE expression, retrotransposition, and chromatin microenvironment may shorten lifespan, elicit neurodegeneration, and impair memory and movement. REs may cause these phenotypes via DNA damage, protein sequestration, insertional mutagenesis, and inflammation. We discuss several TE families, including gypsy, HERV-K, and HERV-W, and how TEs interact with various factors, including transactive response (TAR) DNA-binding protein 43 kDa (TDP-43) and the siRNA and piwi-interacting (pi)RNA systems. Studies of TEs in neurodegeneration have focused on Drosophila and, thus, further examination in mammals is needed. We suggest that therapeutic silencing of REs could help mitigate neurodegenerative disorders.

Transposable elements and their role in aging

Transposable elements (TEs) are an important part of eukaryotic genomes. The role of somatic transposition in aging, carcinogenesis, and other age-related diseases has been determined. This review discusses the fundamental properties of TEs and their complex interactions with cellular processes, which are crucial for understanding the diverse effects of their activity on the genetics and epigenetics of the organism. The interactions of TEs with recombination, replication, repair, and chromosomal regulation; the ability of TEs to maintain a balance between their own activity and repression, the involvement of TEs in the creation of new or alternative genes, the expression of coding/non-coding RNA, and the role in DNA damage and modification of regulatory networks are reviewed. The contribution of the derepressed TEs to age-dependent effects in individual cells/tissues in different organisms was assessed. Conflicting information about TE activity under stress as well as theories of aging mechanisms related to TEs is discussed. On the one hand, transposition activity in response to stressors can lead to organisms acquiring adaptive innovations of great importance for evolution at the population level. On the other hand, the TE expression can cause decreased longevity and stress tolerance at the individual level. The specific features of TE effects on aging processes in germline and soma and the ways of their regulation in cells are highlighted. Recent results considering somatic mutations in normal human and animal tissues are indicated, with the emphasis on their possible functional consequences. In the context of aging, the correlation between somatic TE activation and age-related changes in the number of proteins required for heterochromatin maintenance and longevity regulation was analyzed. One of the original features of this review is a discussion of not only effects based on the TEs insertions and the associated consequences for the germline cell dynamics and somatic genome, but also the differences between transposon- and retrotransposon-mediated structural genome changes and possible phenotypic characteristics associated with aging and various age-related pathologies. Based on the analysis of published data, a hypothesis about the influence of the species-specific features of number, composition, and distribution of TEs on aging dynamics of different animal genomes was formulated.

Adversary of DNA integrity: A long non-coding RNA stimulates driver oncogenic chromosomal rearrangement in human thyroid cells

The flurry of publications devoted to the functions of long non-coding RNAs (lncRNAs) published in the last decade leaves no doubt about the exceptional importance of lncRNAs in various areas including tumor biology. However, contribution of lncRNAs to the early stages of oncogenesis remains poorly understood. In this study we explored a new role for lncRNAs: stimulation of specific chromosomal rearrangements upon DNA damage. We demonstrated that lncRNA CASTL1 (ENSG00000269945) stimulates the formation of the CCDC6-RET inversion (RET/PTC1) in human thyroid cells subjected to radiation or chemical DNA damage. Facilitation of chromosomal rearrangement requires lncRNA to contain regions complementary to the introns of both CCDC6 and RET genes as deletion of these regions deprives CASTL1 of the ability to stimulate the gene fusion. We found that CASTL1 expression is elevated in tumors with CCDC6-RET fusion which is the most frequent rearrangement in papillary thyroid carcinoma. Our results open a new venue for the studies of early oncogenesis in various tumor types, especially those associated with physical or chemical DNA damage.

Full-Length Transcriptomes and Sex-Based Differentially Expressed Genes in the Brain and Ganglia of Giant River Prawn Macrobrachium rosenbergii

Macrobrachium rosenbergii is an important aquaculture prawn that exhibits sexual dimorphism in growth, with males growing much faster than females. However, the mechanisms controlling these complex traits are not well understood. The nervous system plays an important role in regulating life functions. In the present work, we applied PacBio RNA-seq to obtain and characterize the full-length transcriptomes of the brains and thoracic ganglia of female and male prawns, and we performed comparative transcriptome analysis of female and male prawns. A total of 159.1-Gb of subreads were obtained with an average length of 2175 bp and 93.2% completeness. A total of 84,627 high-quality unigenes were generated and annotated with functional databases. A total of 6367 transcript factors and 6287 LncRNAs were predicted. In total, 5287 and 6211 significantly differentially expressed genes (DEGs) were found in the brain and thoracic ganglion, respectively, and confirmed by qRT-PCR. Of the 435 genes associated with protein processing pathways in the endoplasmic reticula, 42 DEGs were detected, and 21/26 DEGs with upregulated expression in the male brain/thoracic ganglion. The DEGs in this pathway are regulated by multiple LncRNAs in polypeptide folding and misfolded protein degradation in the different organs and sexes of the prawn. Our results provide novel theories and insights for studying the nervous system, sexual control, and growth dimorphism.

RNA out of the mist

RNA has long been regarded primarily as the intermediate between genes and proteins. It was a surprise then to discover that eukaryotic genes are mosaics of mRNA sequences interrupted by large tracts of transcribed but untranslated sequences, and that multicellular organisms also express many long 'intergenic' and antisense noncoding RNAs (lncRNAs). The identification of small RNAs that regulate mRNA translation and half-life did not disturb the prevailing view that animals and plant genomes are full of evolutionary debris and that their development is mainly supervised by transcription factors. Gathering evidence to the contrary involved addressing the low conservation, expression, and genetic visibility of lncRNAs, demonstrating their cell-specific roles in cell and developmental biology, and their association with chromatin-modifying complexes and phase-separated domains. The emerging picture is that most lncRNAs are the products of genetic loci termed 'enhancers', which marshal generic effector proteins to their sites of action to control cell fate decisions during development.

Nuclear and cytoplasmic lncRNAs in root tips of the model legume Medicago truncatula under control and submergence

In this study, we aimed to identify long noncoding RNAs (lncRNAs) in root tips of the model legume Medicago truncatula using previously generated nuclear, total polyA, ribosome-associated polyA, and Riboseq RNA datasets, which might shed light on their localization and potential regulatory roles. RNA-seq data were mapped to the version 5 of the M. truncatula A17 genome and analyzed to identify genome annotated lncRNAs and putative new root tip (NRT) lncRNAs. lncRNAs were classified according to their genomic location relative to chromatin accessible regions, protein-coding genes and transposable elements (TE), finding differences between annotated lncRNAs and NRT lncRNAs, both in their genomic position as well as in the type of TEs in their vicinity. We investigated their response to submergence and found a set of regulated lncRNAs that were preferentially upregulated in the nucleus, some of which were located nearby genes of the conserved submergence upregulated gene families, and chromatin accessible regions suggesting a potential regulatory role. Finally, the accumulation of lncRNAs under submergence was validated by reverse transcription quantitative polymerase chain reaction on nuclear RNA, providing additional evidence of their localization, which could ultimately be required for their function.

Genome-wide analysis of long non-coding RNAs (lncRNAs) in tea plants (Camellia sinensis) lateral roots in response to nitrogen application

Tea (Camellia sinensis) is one of the significant cash crops in China. As a leaf crop, nitrogen supply can not only increase the number of new shoots and leaves but also improve the tenderness of the former. However, a conundrum remains in science, which is the molecular mechanism of nitrogen use efficiency, especially long non-coding RNA (lncRNA). In this study, a total of 16,452 lncRNAs were identified through high-throughput sequencing analysis of lateral roots under nitrogen stress and control conditions, of which 9,451 were differentially expressed lncRNAs (DE-lncRNAs). To figure out the potential function of nitrogen-responsive lncRNAs, co-expression clustering was employed between lncRNAs and coding genes. KEGG enrichment analysis revealed nitrogen-responsive lncRNAs may involve in many biological processes such as plant hormone signal transduction, nitrogen metabolism and protein processing in endoplasmic reticulum. The expression abundance of 12 DE-lncRNAs were further verified by RT-PCR, and their expression trends were consistent with the results of RNA-seq. This study expands the research on lncRNAs in tea plants, provides a novel perspective for the potential regulation of lncRNAs on nitrogen stress, and valuable resources for further improving the nitrogen use efficiency of tea plants.

Identification of birch lncRNAs and mRNAs responding to salt stress and characterization of functions of lncRNA

Long noncoding RNAs (lncRNAs) are important in abiotic stress tolerance. Here, we identified salt-responsive genes and lncRNAs in the roots and leaves of Betula platyphylla Suk. (birch), and characterized their lncRNAs functions. In total, 2660 mRNAs and 539 lncRNAs responding to salt treatment were identified using RNA-seq. The salt-responsive genes were substantially enriched in 'cell wall biogenesis' and 'wood development' in the roots and were enriched in 'photosynthesis' and 'response to stimulus' in the leaves. Meanwhile, the potential target genes of the salt-responsive lncRNAs in roots and leaves were both enriched in 'nitrogen compound metabolic process' and 'response to stimulus'. We further built a method for quickly identifying abiotic stress tolerance of lncRNAs, which employed transient transformation for overexpression and knock-down of the lncRNA, enabling gain- and loss-of-function analysis. Using this method, 11 randomly selected salt-responsive lncRNAs were characterized. Among them, six lncRNAs confer salt tolerance, two lncRNAs confer salt sensitivity, and the other three lncRNAs are not involved in salt tolerance. In addition, a lncRNA, LncY1, was further characterized, which improves salt tolerance by regulating two transcription factors, BpMYB96 and BpCDF3. Taken together, our results suggested that lncRNAs play important roles in the salt response of birch plants.

Characteristics of piRNAs and their comparative profiling in testes of sheep with different fertility

As a novel class of small RNAs, piRNAs are highly expressed in the animal gonads and their main known role is to inhibit transposon activity for ensuring the correctness and integrity of genome. In order to explore the characteristics of piRNAs in sheep testis and their possible regulatory roles on male reproduction, deep sequencing technology was used to sequence small RNAs and identify piRNAs in testes of sheep. The length of piRNAs in sheep testes showed a unimodal distribution between 26 and 31 nt, with a peak at 29 nt. These piRNAs exhibited obvious ping-pong signature and strand specificity. In the genome, they were mainly aligned to CDS, intron, repetitive sequence regions and unannotated regions. Furthermore, in transposon analysis, piRNAs were aligned predominantly to LINE, SINE, and LTR types of retrotransposon in sheep testes, and the piRNAs derived from each type showed obvious ping-pong signature. The piRNA clusters identified in sheep testes were mainly distributed on chromosomes 3, 7, 15, 17, 18 and 20. The results combining semen determination with pathway enrichment analysis implied that differentially expressed piRNAs between the testes of rams with different fertility might participate in spermatogenesis by regulating multiple pathways closely related to stabilization of blood-testis barrier and renewal and differentiation of spermatogonial stem cell. Taken together, the study provided new insights into the characteristics, origin and expression patterns of piRNAs in sheep testes tissue, which would help us better understand the role of piRNAs in sheep reproduction.

The dynamics of lncRNAs transcription in interspecific F1 allotriploid hybrids between Brassica species

Interspecific hybridization is the intrinsic forces behind genome evolution. Long non-coding RNAs (lncRNAs) are important for plant biological processes regulation. However, it is unclear that these non-coding fractions are impacted by interspecific hybridization. Here we examined the profiles of lncRNAs by comparing them with coding genes in Brassica napus, three accessions of Brassica rapa, and their F₁ hybrids. 6206 high-confidential lncRNAs were identified in F ₁ hybrids and their parentals, and the lncRNAs transcriptome in the F₁ hybrids was reprogrammed by the genome shock. Notably, genome-wide unbalanced of lncRNAs were observed between A_n and A_r subgenomes, ELD (Expression Level Dominance) was biased toward the A_n -genome in F₁ hybrids, and ELD of non-conserved lncRNAs was more than conserved lncRNAs. Our findings demonstrate that the reprogramed lncRNAs acts as important role in enhancing plant plasticity, leading to the acquisition of desirable traits in polyploid Brassica species.

lncRNAs in fertility: redefining the gene expression paradigm?

Comparative transcriptome approaches assume that highly or dynamically expressed genes are important. This has led to the identification of many genes critical for cellular activity and organism development. However, while testes express the highest levels of long noncoding RNAs (lncRNAs), there is scarcely any evidence for lncRNAs with significant roles in fertility. This was explained by changes in chromatin structure during spermatogenesis that lead to 'promiscuous transcription' with no functional roles for the transcripts. Recent discoveries offer novel and surprising alternatives. Here, I review the current knowledge regarding the involvement of lncRNAs in fertility, why I find gametogenesis different from other developmental processes, offer models to explain why the experimental evidence did not meet theoretical predictions, and suggest possible approaches to test the models.

Unravelling the impact of aging on the human endothelial lncRNA transcriptome

The incidence and prevalence of cardiovascular disease is highest among the elderly. There is a need to further understand the mechanisms behind endothelial cell aging in order to achieve vascular rejuvenation and minimize the onset of age-related vascular diseases. Long non-coding RNAs (lncRNAs) have been proposed to regulate numerous processes in the human genome, yet their function in vascular aging and their therapeutic potential remain largely unknown. This is primarily because the majority of studies investigating the impact of aging on lncRNA expression heavily rely on in vitro studies based on replicative senescence. Here, using a unique collection of young and aged endothelial cells isolated from native human arteries, we sought to characterize the age-related alterations in lncRNA expression profiles. We were able to detect a total of 4463 lncRNAs expressed in the human endothelium from which ∼17% (798) were altered in advanced age. One of the most affected lncRNAs in aging was the primate-specific, Prostate Cancer Associated Transcript (PCAT) 14. In our follow up analysis, using single molecule RNA FISH, we showed that PCAT14 is relatively abundant, localized almost exclusively in the nucleus of young endothelial cells, and silenced in the aged endothelium. Functionally, our studies proposed that downregulation of PCAT14 alters endothelial cell transcription profile and cell functions including endothelial cell migration, sprouting and inflammatory responses in vitro. Taken together, our data highlight that endothelial cell aging correlates with altered expression of lncRNAs, which could impair the endothelial regenerative capacity and enhance inflammatory phenotypes.

Identification of long non-coding RNA-microRNA-mRNA regulatory modules and their potential roles in drought stress response in wheat (Triticum aestivum L.)

Drought is one of the most severe abiotic stresses that influence wheat production across the globe. Understanding the molecular regulatory network of wheat in response to drought is of great importance in molecular breeding. Noncoding RNAs influence plant development and resistance to abiotic stresses by regulating gene expression. In this study, whole-transcriptome sequencing was performed on the seedlings of two wheat varieties with contrasting levels of drought tolerance under drought and control conditions to identify long noncoding RNAs (lncRNAs), micro RNAs (miRNAs), and mRNAs related to drought stress and explore the potential lncRNA-miRNA-mRNA regulatory modules in controlling wheat drought stress response. A total of 1515 differentially expressed lncRNAs (DELs), 209 differentially expressed miRNAs (DEMs), and 20462 differentially expressed genes (DEGs) were identified. Of the 20462 DEGs, 1025 were identified as potential wheat drought resistance-related DEGs. Based on the regulatory relationship and expression patterns of DELs, DEMs, and DEGs, 10 DEL-DEM-DEG regulatory modules related to wheat drought stress response were screened, and preliminary expression verification of two important candidate modules was performed. Our results revealed the possible roles of lncRNA-miRNA-mRNA modules in regulatory networks related to drought tolerance and provided useful information as valuable genomic resources in molecular breeding of wheat.

Differential Expression Pattern of linc-ROR Spliced Variants in Pluripotent and Non-Pluripotent Cell Lines

OBJECTIVE

The human large intergenic non-coding RNA-regulator of reprogramming program (linc-ROR) is known as a stem cell specific linc-RNA. linc-ROR counteracts differentiation via sequestering microRNA-145 (miR-145) that targets OCT4 transcript. Despite the research on the expression and function, the exact structure of linc-ROR transcripts is not clear. Considering the contribution of alternative splicing in transcripts structures and function, identifying different spliced variants of linc-ROR is necessary for further functional analyses. We aimed to find the alternatively spliced transcripts of linc-ROR and investigate their expression pattern in stem and cancer cell lines and during neural differentiation of NT2 cells as a model for understanding linc-ROR role in stem cell and differentiation.

MATERIALS AND METHODS

In this experimental study, linc-ROR locus was scanned for identifying novel exons. Different primer sets were used to detect new spliced variants by reverse transcription polymerase chain reaction (RT-PCR) and direct sequencing. Quantitative PCR (qPCR) and RT-PCR were employed to profile expression of linc-ROR transcripts in different cell lines and during neural differentiation of stem cells.

RESULTS

We could discover 13 novel spliced variants of linc-ROR harboring unique array of exons. Our work uncovered six novel exons, some of which were the product of exonized transposable elements. Monitoring expression profile of the linc-ROR spliced variants in a panel of pluripotent and non-pluripotent cells exhibited that all transcripts were primarily expressed in pluripotent cells. Moreover, the examined linc-ROR spliced variants showed a similar downregulation during neural differentiation of NT2 cells.

CONCLUSION

Altogether, our data showed despite the difference in the structure and composition of exons, various spliced variants of linc-ROR showed similar expression pattern in stem cells and through differentiation.

Molecular diversity and phenotypic pleiotropy of ancient genomic regulatory loci derived from human endogenous retrovirus type H (HERVH) promoter LTR7 and HERVK promoter LTR5_Hs and their contemporary impacts on pathophysiology of Modern Humans

Timelines of population-level effects of viruses on humans varied from the evolutionary scale of million years to contemporary spread of viral infections. Correspondingly, these events are exemplified by: (i) emergence of human endogenous retroviruses (HERVs) from ancient germline infections leading to stable integration of viral genomes into human chromosomes; and (ii) wide-spread viral infections reaching a global pandemic state such as the COVID-19 pandemic. Despite significant efforts, understanding of HERV’s roles in governance of genomic regulatory networks, their impacts on primate evolution and development of human-specific physiological and pathological phenotypic traits remains limited. Remarkably, present analyses revealed that expression of a dominant majority of genes (1696 of 1944 genes; 87%) constituting high-confidence down-steam regulatory targets of defined HERV loci was significantly altered in cells infected with the SARS-CoV-2 coronavirus, a pathogen causing the global COVID-19 pandemic. This study focused on defined sub-sets of DNA sequences derived from HERVs that are expressed at specific stages of human preimplantation embryogenesis and exert regulatory actions essential for self-renewal and pluripotency. Evolutionary histories of LTR7/HERVH and LTR5_Hs/HERVK were charted based on evidence of the earliest presence and expansion of highly conserved (HC) LTR sequences. Sequence conservation analyses of most recent releases 17 primate species’ genomes revealed that LTR7/HERVH have entered germlines of primates in Africa after the separation of the New World Monkey lineage, while LTR5_Hs/HERVK successfully colonized primates’ germlines after the segregation of Gibbons’ species. Subsequently, both LTR7 and LTR5_Hs undergo a marked ~ fourfold–fivefold expansion in genomes of Great Apes. Timelines of quantitative expansion of both LTR7 and LTR5_Hs loci during evolution of Great Apes appear to replicate the consensus evolutionary sequence of increasing cognitive and behavioral complexities of non-human primates, which seems particularly striking for LTR7 loci and 11 distinct LTR7 subfamilies. Consistent with previous reports, identified in this study, 351 human-specific (HS) insertions of LTR7 (175 loci) and LTR5_Hs (176 loci) regulatory sequences have been linked to genes implicated in establishment and maintenance of naïve and primed pluripotent states and preimplantation embryogenesis phenotypes. Unexpectedly, HS-LTRs manifest regulatory connectivity to genes encoding markers of 12 distinct cells’ populations of fetal gonads, as well as genes implicated in physiology and pathology of human spermatogenesis, including Y-linked spermatogenic failure, oligo- and azoospermia. Granular interrogations of genes linked with 11 distinct LTR7 subfamilies revealed that mammalian offspring survival (MOS) genes seem to remain one of consistent regulatory targets throughout ~ 30 MYA of the divergent evolution of LTR7 loci. Differential GSEA of MOS versus non-MOS genes identified clearly discernable dominant enrichment patterns of phenotypic traits affected by MOS genes linked with LTR7 (562 MOS genes) and LTR5_Hs (126 MOS genes) regulatory loci across the large panel of genomics and proteomics databases reflecting a broad spectrum of human physiological and pathological traits. GSEA of LTR7-linked MOS genes identified more than 2200 significantly enriched records of human common and rare diseases and gene signatures of 466 significantly enriched records of Human Phenotype Ontology traits, including Autosomal Dominant (92 genes) and Autosomal Recessive (93 genes) Inheritance. LTR7 regulatory elements appear linked with genes implicated in functional and morphological features of central nervous system, including synaptic transmission and protein–protein interactions at synapses, as well as gene signatures differentially regulated in cells of distinct neurodevelopmental stages and morphologically diverse cell types residing and functioning in human brain. These include Neural Stem/Precursor cells, Radial Glia cells, Bergman Glia cells, Pyramidal cells, Tanycytes, Immature neurons, Interneurons, Trigeminal neurons, GABAergic neurons, and Glutamatergic neurons. GSEA of LTR7-linked genes identified significantly enriched gene sets encoding markers of more than 80 specialized types of neurons and markers of 521 human brain regions, most prominently, subiculum and dentate gyrus. Identification and characterization of 1944 genes comprising high-confidence down-steam regulatory targets of LTR7 and/or LTR5_Hs loci validated and extended these observations by documenting marked enrichments for genes implicated in neoplasm metastasis, intellectual disability, autism, multiple cancer types, Alzheimer’s, schizophrenia, and other brain disorders. Overall, genes representing down-stream regulatory targets of ancient retroviral LTRs exert the apparently cooperative and exceedingly broad phenotypic impacts on human physiology and pathology. This is exemplified by altered expression of 93% high-confidence LTR targets in cells infected by contemporary viruses, revealing a convergence of virus-inflicted aberrations on genomic regulatory circuitry governed by ancient retroviral LTR elements and interference with human cells’ differentiation programs.

Mammalian genome innovation through transposon domestication

Since the discovery of transposons, their sheer abundance in host genomes has puzzled many. While historically viewed as largely harmless 'parasitic' DNAs during evolution, transposons are not a mere record of ancient genome invasion. Instead, nearly every element of transposon biology has been integrated into host biology. Here we review how host genome sequences introduced by transposon activities provide raw material for genome innovation and document the distinct evolutionary path of each species.

Induced pluripotent stem cells display a distinct set of MHC I-associated peptides shared by human cancers

Previous reports showed that mouse vaccination with pluripotent stem cells (PSCs) induces durable anti-tumor immune responses via T cell recognition of some elusive oncofetal epitopes. We characterize the MHC I-associated peptide (MAP) repertoire of human induced PSCs (iPSCs) using proteogenomics. Our analyses reveal a set of 46 pluripotency-associated MAPs (paMAPs) absent from the transcriptome of normal tissues and adult stem cells but expressed in PSCs and multiple adult cancers. These paMAPs derive from coding and allegedly non-coding (48%) transcripts involved in pluripotency maintenance, and their expression in The Cancer Genome Atlas samples correlates with source gene hypomethylation and genomic aberrations common across cancer types. We find that several of these paMAPs were immunogenic. However, paMAP expression in tumors coincides with activation of pathways instrumental in immune evasion (WNT, TGF-β, and CDK4/6). We propose that currently available inhibitors of these pathways could synergize with immune targeting of paMAPs for the treatment of poorly differentiated cancers.

Mutant KRAS regulates transposable element RNA and innate immunity via KRAB zinc-finger genes

RAS genes are the most frequently mutated oncogenes in cancer, yet the effects of oncogenic RAS signaling on the noncoding transcriptome remain unclear. We analyzed the transcriptomes of human airway and bronchial epithelial cells transformed with mutant KRAS to define the landscape of KRAS-regulated noncoding RNAs. We find that oncogenic KRAS signaling upregulates noncoding transcripts throughout the genome, many of which arise from transposable elements (TEs). These TE RNAs exhibit differential expression, are preferentially released in extracellular vesicles, and are regulated by KRAB zinc-finger (KZNF) genes, which are broadly downregulated in mutant KRAS cells and lung adenocarcinomas in vivo. Moreover, mutant KRAS induces an intrinsic IFN-stimulated gene (ISG) signature that is often seen across many different cancers. Our results indicate that mutant KRAS remodels the repetitive noncoding transcriptome, demonstrating the broad scope of intracellular and extracellular RNAs regulated by this oncogenic signaling pathway.

Many human cancers are driven by mutant KRAS, but its effects on noncoding RNA are unclear. Reggiardo et al. show that mutant KRAS regulates this RNA landscape by silencing KRAB zinc-finger genes that normally repress transposable element noncoding RNAs, which are preferentially released from mutant KRAS cells in extracellular vesicles.

Roles of transposable elements in the regulation of mammalian transcription

Transposable elements (TEs) comprise about half of the mammalian genome. TEs often contain sequences capable of recruiting the host transcription machinery, which they use to express their own products and promote transposition. However, the regulatory sequences carried by TEs may affect host transcription long after the TEs have lost the ability to transpose. Recent advances in genome analysis and engineering have facilitated systematic interrogation of the regulatory activities of TEs. In this Review, we discuss diverse mechanisms by which TEs contribute to transcription regulation. Notably, TEs can donate enhancer and promoter sequences that influence the expression of host genes, modify 3D chromatin architecture and give rise to novel regulatory genes, including non-coding RNAs and transcription factors. We discuss how TEs spur regulatory evolution and facilitate the emergence of genetic novelties in mammalian physiology and development. By virtue of their repetitive and interspersed nature, TEs offer unique opportunities to dissect the effects of mutation and genomic context on the function and evolution of cis-regulatory elements. We argue that TE-centric studies hold the key to unlocking general principles of transcription regulation and evolution.

Zebrafish transposable elements show extensive diversification in age, genomic distribution, and developmental expression

There is considerable interest in understanding the effect of transposable elements (TEs) on embryonic development. Studies in humans and mice are limited by the difficulty of working with mammalian embryos and by the relative scarcity of active TEs in these organisms. The zebrafish is an outstanding model for the study of vertebrate development, and over half of its genome consists of diverse TEs. However, zebrafish TEs remain poorly characterized. Here we describe the demography and genomic distribution of zebrafish TEs and their expression throughout embryogenesis using bulk and single-cell RNA sequencing data. These results reveal a highly dynamic genomic ecosystem comprising nearly 2000 distinct TE families, which vary in copy number by four orders of magnitude and span a wide range of ages. Longer retroelements tend to be retained in intergenic regions, whereas short interspersed nuclear elements (SINEs) and DNA transposons are more frequently found nearby or within genes. Locus-specific mapping of TE expression reveals extensive TE transcription during development. Although two-thirds of TE transcripts are likely driven by nearby gene promoters, we still observe stage- and tissue-specific expression patterns in self-regulated TEs. Long terminal repeat (LTR) retroelements are most transcriptionally active immediately following zygotic genome activation, whereas DNA transposons are enriched among transcripts expressed in later stages of development. Single-cell analysis reveals several endogenous retroviruses expressed in specific somatic cell lineages. Overall, our study provides a valuable resource for using zebrafish as a model to study the impact of TEs on vertebrate development.

ARID1A loss derepresses a group of human endogenous retrovirus-H loci to modulate BRD4-dependent transcription

Transposable elements (TEs) through evolutionary exaptation have become an integral part of the human genome, offering ample regulatory sequences and shaping chromatin 3D architecture. While the functional impacts of TE-derived sequences on early embryogenesis have been recognized, their roles in malignancy are only starting to emerge. Here we show that many TEs, especially the pluripotency-related human endogenous retrovirus H (HERVH), are abnormally activated in colorectal cancer (CRC) samples. Transcriptional upregulation of HERVH is associated with mutations of several tumor suppressors, particularly ARID1A. Knockout of ARID1A in CRC cells leads to increased transcription at several HERVH loci, which involves compensatory contribution by ARID1B. Suppression of HERVH in CRC cells and patient-derived organoids impairs tumor growth. Mechanistically, HERVH transcripts colocalize with nuclear BRD4 foci, modulating their dynamics and co-regulating many target genes. Altogether, we uncover a critical role for ARID1A in restraining HERVH, whose abnormal activation can promote tumorigenesis by stimulating BRD4-dependent transcription.

Here the authors show mutation of the BAF chromatin remodeler subunit ARID1A results in an ARID1B-dependent upregulation of HERVH, an ERV required for the pluripotency regulatory network. These HERVH RNAs can partition into BRD4 foci, affecting BRD4-dependent transcription. Suppression of HERVH in colorectal cancer cells and patient-derived organoids impairs tumor growth.

Expression of Extracellular Vesicle PIWI-Interacting RNAs Throughout hiPSC-Cardiomyocyte Differentiation

Extracellular Vesicles (EV) play a critical role in the regulation of regenerative processes in wounded tissues by mediating cell-to-cell communication. Multiple RNA species have been identified in EV, although their function still lacks understanding. We previously characterized the miRNA content of EV secreted over hiPSC-cardiomyocyte differentiation and found a distinct miRNA expression in hiPSC-EV driving its in vitro bioactivity. In this work, we investigated the piRNA profiles of EV derived from key stages of the hiPSC-CM differentiation and maturation, i.e., from hiPSC (hiPSC-EV), cardiac progenitors (CPC-EV), immature (CMi-EV), and mature (CMm-EV) cardiomyocytes, demonstrating that EV-piRNA expression differs greatly from the miRNA profiles we previously identified. Only four piRNA were significantly deregulated in EV, one in hiPSC-EV, and three in CPC-EV, as determined by differential expression analysis on small RNA-seq data. Our results provide a valuable source of information for further studies aiming at defining the role of piRNA in the bioactivity and therapeutic potential of EV.

Molecular genetics of idiopathic pulmonary fibrosis

Idiopathic pulmonary fibrosis (IPF) is a severe progressive interstitial lung disease with a prevalence of 2 to 29 per 100,000 of the world’s population. Aging is a significant risk factor for IPF, and the mechanisms of aging (telomere depletion, genomic instability, mitochondrial dysfunction, loss of proteostasis) are involved in the pathogenesis of IPF. The pathogenesis of IPF consists of TGF-β activation, epithelial-mesenchymal transition, and SIRT7 expression decrease. Genetic studies have shown a role of mutations and polymorphisms in mucin genes (MUC5B), in the genes responsible for the integrity of telomeres (TERC, TERC, TINF2, DKC1, RTEL1, PARN), in surfactant-related genes (SFTPC, SFTPCA, SFTPA2, ABCA3, SP-A2), immune system genes (IL1RN, TOLLIP), and haplotypes of HLA genes (DRB1*15:01, DQB1*06:02) in IPF pathogenesis. The investigation of the influence of reversible epigenetic factors on the development of the disease, which can be corrected by targeted therapy, shows promise. Among them, an association of a number of specific microRNAs and long noncoding RNAs was revealed with IPF. Therefore, dysregulation of transposons, which serve as key sources of noncoding RNA and affect mechanisms of aging, may serve as a driver for IPF development. This is due to the fact that pathological activation of transposons leads to violation of the regulation of genes, in the epigenetic control of which microRNA originating from these transposons are involved (due to the complementarity of nucleotide sequences). Analysis of the MDTE database (miRNAs derived from Transposable Elements) allowed the detection of 12 different miRNAs derived in evolution from transposons and associated with IPF (miR-31, miR-302, miR-326, miR-335, miR-340, miR-374, miR-487, miR-493, miR-495, miR-630, miR-708, miR-1343). We described the relationship of transposons with TGF-β, sirtuins and telomeres, dysfunction of which is involved in the pathogenesis of IPF. New data on IPF epigenetic mechanisms can become the basis for improving results of targeted therapy of the disease using noncoding RNAs.

Transposable Elements in Pluripotent Stem Cells and Human Disease

Transposable elements (TEs) are mobile genetic elements that can randomly integrate into other genomic sites. They have successfully replicated and now occupy around 40% of the total DNA sequence in humans. TEs in the genome have a complex relationship with the host cell, being both potentially deleterious and advantageous at the same time. Only a tiny minority of TEs are still capable of transposition, yet their fossilized sequence fragments are thought to be involved in various molecular processes, such as gene transcriptional activity, RNA stability and subcellular localization, and chromosomal architecture. TEs have also been implicated in biological processes, although it is often hard to reveal cause from correlation due to formidable technical issues in analyzing TEs. In this review, we compare and contrast two views of TE activity: one in the pluripotent state, where TEs are broadly beneficial, or at least mechanistically useful, and a second state in human disease, where TEs are uniformly considered harmful.

Meta-Analysis Suggests That Intron Retention Can Affect Quantification of Transposable Elements from RNA-Seq Data

Simple Summary

Transposable elements (TEs) are repetitive sequences comprising more than one third of the human genome with the original ability to change their location within the genome. Owing to their repetitive nature, the quantification of TEs results often challenging. RNA-seq is a useful tool for genome-wide TEs quantification, nevertheless it also presents technical issues, including low reads mappability and erroneous quantification derived from the transcription of TEs fragments embedded in canonical transcripts. Fragments derived from TEs are found within the introns of most genes, which led to the hypothesis that intron retention (IR) can affect the unbiased quantification of TEs expression. Performing meta-analysis of public RNA-seq datasets, here we observe that IR can indeed impact the quantification of TEs by increasing the number of reads mapped on intronic TE copies. Our work highlights a correlation between IR and TEs expression measurement by RNA-seq that should be taken into account to achieve reliable TEs quantification, especially in samples characterized by extensive IR, because differential IR might be confused with differential TEs expression.

Abstract

Transposable elements (TEs), also known as “jumping genes”, are repetitive sequences with the capability of changing their location within the genome. They are key players in many different biological processes in health and disease. Therefore, a reliable quantification of their expression as transcriptional units is crucial to distinguish between their independent expression and the transcription of their sequences as part of canonical transcripts. TEs quantification faces difficulties of different types, the most important one being low reads mappability due to their repetitive nature preventing an unambiguous mapping of reads originating from their sequences. A large fraction of TEs fragments localizes within introns, which led to the hypothesis that intron retention (IR) can be an additional source of bias, potentially affecting accurate TEs quantification. IR occurs when introns, normally removed from the mature transcript by the splicing machinery, are maintained in mature transcripts. IR is a widespread mechanism affecting many different genes with cell type-specific patterns. We hypothesized that, in an RNA-seq experiment, reads derived from retained introns can introduce a bias in the detection of overlapping, independent TEs RNA expression. In this study we performed meta-analysis using public RNA-seq data from lymphoblastoid cell lines and show that IR can impact TEs quantification using established tools with default parameters. Reads mapped on intronic TEs were indeed associated to the expression of TEs and influence their correct quantification as independent transcriptional units. We confirmed these results using additional independent datasets, demonstrating that this bias does not appear in samples where IR is not present and that differential TEs expression does not impact on IR quantification. We concluded that IR causes the over-quantification of intronic TEs and differential IR might be confused with differential TEs expression. Our results should be taken into account for a correct quantification of TEs expression from RNA-seq data, especially in samples in which IR is abundant.

The Role of Transposable Elements of the Human Genome in Neuronal Function and Pathology

Transposable elements (TEs) have been extensively studied for decades. In recent years, the introduction of whole-genome and whole-transcriptome approaches, as well as single-cell resolution techniques, provided a breakthrough that uncovered TE involvement in host gene expression regulation underlying multiple normal and pathological processes. Of particular interest is increased TE activity in neuronal tissue, and specifically in the hippocampus, that was repeatedly demonstrated in multiple experiments. On the other hand, numerous neuropathologies are associated with TE dysregulation. Here, we provide a comprehensive review of literature about the role of TEs in neurons published over the last three decades. The first chapter of the present review describes known mechanisms of TE interaction with host genomes in general, with the focus on mammalian and human TEs; the second chapter provides examples of TE exaptation in normal neuronal tissue, including TE involvement in neuronal differentiation and plasticity; and the last chapter lists TE-related neuropathologies. We sought to provide specific molecular mechanisms of TE involvement in neuron-specific processes whenever possible; however, in many cases, only phenomenological reports were available. This underscores the importance of further studies in this area.