Inferring the expression variability of human transposable element-derived exons by linear model analysis of deep RNA sequencing data

Regulation of human interferon signaling by transposon exonization

Innate immune signaling is essential for clearing pathogens and damaged cells, and must be tightly regulated to avoid excessive inflammation or autoimmunity. Here, we found that the alternative splicing of exons derived from transposable elements is a key mechanism controlling immune signaling in human cells. By analyzing long-read transcriptome datasets, we identified numerous transposon exonization events predicted to generate functional protein variants of immune genes, including the type I interferon receptor IFNAR2. We demonstrated that the transposon-derived isoform of IFNAR2 is more highly expressed than the canonical isoform in almost all tissues, and functions as a decoy receptor that potently inhibits interferon signaling including in cells infected with SARS-CoV-2. Our findings uncover a primate-specific axis controlling interferon signaling and show how a transposon exonization event can be co-opted for immune regulation.

Widespread splicing of repetitive element loci into coding regions of gene transcripts

We performed a thorough characterization of expressed repetitive element loci (RE) in the human orbitofrontal cortex (OFC) using directional RNA sequencing data. Considering only sequencing reads that map uniquely onto the human genome, we discovered that the overwhelming majority of intronic and exonic RE are expressed in the same orientation as the gene in which they reside. Our mapping approach enabled the identification of novel differentially expressed RE transcripts between the OFC and peripheral blood lymphocytes. Further analysis revealed that RE are extensively spliced into coding regions of gene transcripts yielding thousands of novel mRNA variants with altered coding potential. Lower frequency splicing of RE into untranslated regions of gene transcripts was also observed. The same pattern of RE splicing in the brain was also detected for Drosophila, zebrafish, mouse, rat, dog and rabbit. RE splicing occurs largely at canonical GT-AG splice junctions with LINE and SINE elements forming the most RE splice junctions in the human OFC. This type of splicing usually gives rise to a minor splice variant of the endogenous gene and in silico analysis suggests that RE splicing has the potential to introduce novel open reading frames. Reanalysis of previously published sequencing data performed in the mouse cerebellum revealed that thousands of RE splice variants are associated with translating ribosomes. Our results demonstrate that RE expression is more complex than previously envisioned and raise the possibility that RE splicing might generate functional protein isoforms.

De novo sequencing, assembly and analysis of eight different transcriptomes from the Malayan pangolin

Pangolins are scale-covered mammals, containing eight endangered species. Maintaining pangolins in captivity is a significant challenge, in part because little is known about their genetics. Here we provide the first large-scale sequencing of the critically endangered Manis javanica transcriptomes from eight different organs using Illumina HiSeq technology, yielding ~75 Giga bases and 89,754 unigenes. We found some unigenes involved in the insect hormone biosynthesis pathway and also 747 lipids metabolism-related unigenes that may be insightful to understand the lipid metabolism system in pangolins. Comparative analysis between M. javanica and other mammals revealed many pangolin-specific genes significantly over-represented in stress-related processes, cell proliferation and external stimulus, probably reflecting the traits and adaptations of the analyzed pregnant female M. javanica. Our study provides an invaluable resource for future functional works that may be highly relevant for the conservation of pangolins.

Evolutionary impact of transposable elements on genomic diversity and lineage-specific innovation in vertebrates

Since their discovery, a growing body of evidence has emerged demonstrating that transposable elements are important drivers of species diversity. These mobile elements exhibit a great variety in structure, size and mechanisms of transposition, making them important putative actors in organism evolution. The vertebrates represent a highly diverse and successful lineage that has adapted to a wide range of different environments. These animals also possess a rich repertoire of transposable elements, with highly diverse content between lineages and even between species. Here, we review how transposable elements are driving genomic diversity and lineage-specific innovation within vertebrates. We discuss the large differences in TE content between different vertebrate groups and then go on to look at how they affect organisms at a variety of levels: from the structure of chromosomes to their involvement in the regulation of gene expression, as well as in the formation and evolution of non-coding RNAs and protein-coding genes. In the process of doing this, we highlight how transposable elements have been involved in the evolution of some of the key innovations observed within the vertebrate lineage, driving the group's diversity and success.

Transcriptome profiling of the cancer and adjacent nontumor tissues from cervical squamous cell carcinoma patients by RNA sequencing

Cervical cancer is the third most common cancer and the fourth leading cause of cancer deaths among women in the world. The discovery of vital diagnostic and therapeutic markers against cervical squamous cell carcinoma (CSCC) would broaden our understanding on the molecular basis of CSCC. In this study, we thoroughly analyzed the transcriptome of CSCC and matched adjacent nontumor (ATN) tissue. RNA sequencing was performed to screen the differentially expressed genes (DEGs) of three pairs of CSCC and ATN tissues. Functional enrichment analysis was used to uncover the biological functions of DEGs. Protein interaction network was carried out to reveal interaction of DEGs. Quantitative real-time PCR was conducted to validate the expression of DEGs. Immunohistochemistry was used to detect the relationship between clinicopathological parameters of CSCC and DEGs. There were a total of 347 significantly common DEGs in the three paired examples, including 104 consistent upregulated and 148 consistent downregulated DEGs. The 347 DEGs were categorized into 73 functional categories by Gene Ontology (GO) analysis. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis suggested six significantly signal pathways. The protein interaction network uncovered three important DEGs, including retinol dehydrogenase 12 (RDH12), ubiquitin D (UBD), and serum amyloid A1 (SAA1). We found that RDH12 expression was decreased in 74.5 % of CSCC tissues. RDH12 expression was negatively associated with tumor size and depth of cervical invasion. The UBD was overexpressed in 61.7 % of CSCC tissues and was positively related with tumor size and lymphatic metastasis. The SAA1 protein was overexpressed in 57.4 % of CSCC tissues and was positively related with clinicopathological parameters of tumor size, lymphatic metastasis, and depth of cervical invasion. The RDH12, UBD, and SAA1 genes might participate in the progression of CSCC.