Using Pan RNA-Seq Analysis to Reveal the Ubiquitous Existence of 5' and 3' End Small RNAs

Expression of Mitochondrial Long Non-Coding RNAs, MDL1 and MDL1AS, Are Good Prognostic and/or Diagnostic Biomarkers for Several Cancers, Including Colorectal Cancer

Non-coding RNAs provide new opportunities to identify biomarkers that properly classify cancer patients. Here, we study the biomarker status of the mitochondrial long non-coding RNAs, MDL1 and MDL1AS. Expression of these genes was studied in public transcriptomic databases. In addition, a cohort of 69 locally advanced rectal cancer (LARC) patients with a follow-up of more than 5 years was used to determine the prognostic value of these markers. Furthermore, cell lines of colorectal (HCT116) and breast (MDA-MB-231) carcinoma were employed to study the effects of downregulating MDL1AS in vitro. Expression of MDL1AS (but not MDL1) was significantly different in tumor cells than in the surrounding tissue in a tumor-type-specific context. Both MDL1 and MDL1AS were accurate biomarkers for the 5-year survival of LARC patients (p = 0.040 and p = 0.007, respectively) with promising areas under the curve in the ROC analyses (0.820 and 0.930, respectively). MDL1AS downregulation reduced mitochondrial respiration in both cell lines. Furthermore, this downregulation produced a decrease in growth and migration on colorectal cells, but the reverse effects on breast cancer cells. In summary, MDL1 and MDL1AS can be used as reliable prognostic biomarkers of LARC, and MDL1AS expression provides relevant information on the diagnosis of different cancers.

Emerging functions of mitochondria-encoded noncoding RNAs

Mitochondria, organelles that harbor their own circular genomes, are critical for energy production and homeostasis maintenance in eukaryotic cells. Recent studies discovered hundreds of mitochondria-encoded noncoding RNAs (mt-ncRNAs), including novel subtypes of mitochondria-encoded circular RNAs (mecciRNAs) and mitochondria-encoded double-stranded RNAs (mt-dsRNAs). Here, we discuss the emerging field of mt-ncRNAs by reviewing their expression patterns, biogenesis, metabolism, regulatory roles, and functional mechanisms. Many mt-ncRNAs have regulatory roles in cellular physiology, and some are associated with, or even act as, causal factors in human diseases. We also highlight developments in technologies and methodologies and further insights into future perspectives and challenges in studying these noncoding RNAs, as well as their potential biomedical applications.

The first discovery of Tc1 transposons in yeast

Background

Identification of transposons without close homologs is still a difficult task. IS630/Tc1/mariner transposons, classified into a superfamily, are probably the most widespread DNA transposons in nature. Tc1/mariner transposons have been discovered in animals, plants, and filamentous fungi, however, not in yeast.

Results

In the present study, we report the discovery of two intact Tc1 transposons in yeast and filamentous fungi, respectively. The first one, named Tc1-OP1 (DD40E), represents Tc1 transposons in Ogataea parapolymorpha. The second one, named Tc1-MP1 (DD34E), represents Tc1 transposons in the Rhizopodaceae and Mucoraceae families. As a homolog of Tc1-OP1 and Tc1-MP1, IS630-AB1 (DD34E) was discovered as an IS630 transposon in Acinetobacter spp.

Conclusion

Tc1-OP1 is not only the first reported Tc1 transposon in yeast, but also the first reported nonclassical Tc1 transposon. Tc1-OP1 is the largest of IS630/Tc1/mariner transposons reported to date and significantly different from others. Notably, Tc1-OP1 encodes a serine-rich domain and a transposase, extending the current knowledge of Tc1 transposons. The phylogenetic relationships of Tc1-OP1, Tc1-MP1 and IS630-AB1 indicated that these transposons had evolved from a common ancestor. Tc1-OP1, Tc1-MP1 and IS630-AB1 can be used as reference sequences to facilitate the identification of IS630/Tc1/mariner transposons. More Tc1/mariner transposons will be identified in yeast, following our discovery.

Complete Mitogenome Analysis of Five Leafhopper Species of Idiocerini (Hemiptera: Cicadellidae)

Insect mitochondrial genomes (mitogenomes) are of great interest in exploring molecular evolution, phylogenetics, and biogeography. So far, only 12 mitogenomes of the leafhopper tribe Idiocerini have been released in GenBank, although the tribe comprises 488 known species including some agricultural, forestry, and horticultural pests. In order to compare and analyze the mitochondrial genome structure of Idiocerini and even the selective pressure of 13 protein-coding genes (PCGs) of the family Cicadellidae, the complete mitogenomes of five species including Nabicerus dentimus, Sahlbergotettix salicicola, Podulmorinus opacus, Podulmorinus consimilis, and a new species of a new genus were determined by next-generation sequencing. The size of the newly determined mitogenomes ranged from 14,733 bp to 15,044 bp, comprising the standard set of 13 PCGs, 22 transfer RNA genes, two ribosomal RNA genes, and a long non-coding control region (CR). The extent of purifying selection presented different pictures in the tribe and the family. The less pronounced genes (0.5 < dN/dS < 1) were nad5 and nad4l in Idiocerin, whereas in the family Cicadellidae including the sequences of Idiocerin, nad1-nad6 and cox1 genes were less pronounced. The codon encoding leucine was the most common in all species, and the codon encoding serine 1 was the most common in all species except for P. opacus. Interestingly, in P. opacus, another of the most common codons is that encoding serine 2. Among the 17 examined species of the Idiocerini, 14 species contained the tandem repeats, and 11 species of them contained the motif "TTATA". These findings will promote research on the structure and evolution of the mitochondrial genome and highlight the need for more mitogenomes in Cicadellidae.

Full-length chloroplast genome of Dongxiang wild rice reveals small single-copy region switching

BACKGROUND

Plant chloroplast DNA (cpDNA) typically has a circular structure, including a large single-copy region (LSC), a small single-copy region (SSC) and two inverted repeats (IR1 and IR2). The organization of these four elementary regions LSC-IR1-SSC-IR2 is highly conserved across all plant cpDNAs. Very few structural variations (SVs) occurring at the elementary-region level have been reported.

RESULTS

In the present study, we assembled the full-length cpDNA of Dongxiang wild rice line 159 (DXWR159). Using the long PacBio subreads, we discovered a large inversion of SSC and a large duplication of IR in DXWR159 cpDNAs. Significantly, we reported for the first time forward and reverse SSCs of cpDNAs in similar proportions and named the frequent inversion of a whole SSC as SSC switching.

CONCLUSIONS

Our study helps researchers to correctly assemble the chloroplast genomes. Our recombination model explained the formation of large SVs in cpDNAs and provided insights into a novel scientific question that if there are common mechanisms in the formation or translocation of all kinds of transposon-like elements (TLEs). We propose that: (1) large inversion is the most accepted mutation type of SVs in cpDNAs; (2) SSC switching ubiquitous occurs in plant cpDNAs; and (3) further investigation of molecular mechanism underlying SSC switching may reveal new driving forces for large SVs.

Precise annotation of Drosophila mitochondrial genomes leads to insights into AT-rich regions

In the present study, we performed precise annotation of Drosophila melanogaster, D. simulans, D. grimshawi, Bactrocera oleae mitochondrial (mt) genomes using pan RNA-seq analysis. Using PacBio cDNA-seq data from D. simulans, we precisely annotated the Transcription Initiation Sites (TISs) of the mt Heavy and Light strands in Drosophila mt genomes and reported that the polyA(+) and polyA(-) motifs in the CRs are associated with TISs. The discovery of the conserved polyA(+) and polyA(-) motifs provides insights into many polyA and polyT sequences in CRs of insect mt genomes, leading to reveal the mt transcription and its regulation in invertebrates. Notably, we propose that: (1) polyA/polyT motifs in CRs function as signals to initiate mtDNA transcription; (2) the duplication, recombination or mutation of these polyA/polyT sequences formed the AT-rich regions during evolution; and (3) since CRs of many invertebrate species still contain many polyA/polyT sequences, there is a high probability that several TISs and TTSs exist in invertebrate mt genomes.

How the Replication and Transcription Complex Functions in Jumping Transcription of SARS-CoV-2

Background: Coronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Although unprecedented efforts are underway to develop therapeutic strategies against this disease, scientists have acquired only a little knowledge regarding the structures and functions of the CoV replication and transcription complex (RTC). Ascertaining all the RTC components and the arrangement of them is an indispensably step for the eventual determination of its global structure, leading to completely understanding all of its functions at the molecular level.

Results: The main results include: 1) hairpins containing the canonical and non-canonical NSP15 cleavage motifs are canonical and non-canonical transcription regulatory sequence (TRS) hairpins; 2) TRS hairpins can be used to identify recombination regions in CoV genomes; 3) RNA methylation participates in the determination of the local RNA structures in CoVs by affecting the formation of base pairing; and 4) The eventual determination of the CoV RTC global structure needs to consider METTL3 in the experimental design.

Conclusions: In the present study, we proposed the theoretical arrangement of NSP12-15 and METTL3 in the global RTC structure and constructed a model to answer how the RTC functions in the jumping transcription of CoVs. As the most important finding, TRS hairpins were reported for the first time to interpret NSP15 cleavage, RNA methylation of CoVs and their association at the molecular level. Our findings enrich fundamental knowledge in the field of gene expression and its regulation, providing a crucial basis for future studies.

Complete Mitochondrial Genomes of Five Racerunners (Lacertidae: Eremias) and Comparison with Other Lacertids: Insights into the Structure and Evolution of the Control Region

Comparative studies on mitochondrial genomes (mitogenomes) as well as the structure and evolution of the mitochondrial control region are few in the Lacertidae family. Here, the complete mitogenomes of five individuals of Eremias scripta (2 individuals), Eremias nikolskii, Eremias szczerbaki, and Eremias yarkandensis were determined using next-generation sequencing and were compared with other lacertids available in GenBank. The circular mitogenomes comprised the standard set of 13 protein-coding genes (PCGs), 22 transfer RNA genes, 2 ribosomal RNA genes and a long non-coding control region (CR). The extent of purifying selection was less pronounced for the COIII and ND2 genes in comparison with the rest of the PCGs. The codons encoding Leucine (CUN), Threonine, and Isolecucine were the three most frequently present. The secondary structure of rRNA of Lacertidae (herein, E. scripta KZL15 as an example) comprised four domains and 28 helices for 12S rRNA, with six domains and 50 helices for 16S rRNA. Five types and twenty-one subtypes of CR in Lacertidae were described by following the criteria of the presence and position of tandem repeats (TR), termination-associated sequence 1 (TAS1), termination-associated sequence 2 (TAS2), conserved sequence block 1 (CBS1), conserved sequence block 2 (CSB2), and conserved sequence block 3 (CSB3). The compositions of conserved structural elements in four genera, Acanthodactylus, Darevskia, Eremias, and Takydromus, were further explored in detail. The base composition of TAS2 – TATACATTAT in Lacertidae was updated. In addition, the motif “TAGCGGCTTTTTTG” of tandem repeats in Eremias and the motif ”GCGGCTT” in Takydromus were presented. Nucleotide lengths between CSB2 and CSB3 remained 35 bp in Eremias and Darevskia. The phylogenetic analyses of Lacertidae recovered the higher-level relationships among the three subfamilies and corroborated a hard polytomy in the Lacertinae phylogeny. The phylogenetic position of E. nikolskii challenged the monophyly of the subgenus Pareremias within Eremias. Some mismatches between the types of CR and their phylogeny demonstrated the complicated evolutionary signals of CR such as convergent evolution. These findings will promote research on the structure and evolution of the CR and highlight the need for more mitogenomes in Lacertidae.

Mitochondria Encoded Non-coding RNAs in Cell Physiology

Mitochondria are the powerhouses of mammalian cells, which participate in series of metabolic processes and cellular events. Mitochondria have their own genomes, and it is generally acknowledged that human mitochondrial genome encodes 13 proteins, 2 rRNAs and 22 tRNAs. However, the complexity of mitochondria derived transcripts is just starting to be envisaged. Currently, there are at least 8 lncRNAs, some dsRNAs, various small RNAs, and hundreds of circRNAs known to be generated from mitochondrial genome. These non-coding RNAs either translocate into cytosol/nucleus or reside in mitochondria to play various biological functions. Here we present an overview of regulatory non-coding RNAs encoded by the mammalian mitochondria genome. For overall understandings of non-coding RNAs in mitochondrial function, a brief summarization of nuclear-encoded non-coding RNAs in mitochondria is also included. We discuss about roles of these non-coding RNAs in cellular physiology and the communication between mitochondria and the nucleus.

A Negative Feedback Model to Explain Regulation of SARS-CoV-2 Replication and Transcription

BACKGROUND

Coronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Although a preliminary understanding of the replication and transcription of SARS-CoV-2 has recently emerged, their regulation remains unknown.

RESULTS

By comprehensive analysis of genome sequence and protein structure data, we propose a negative feedback model to explain the regulation of CoV replication and transcription, providing a molecular basis of the "leader-to-body fusion" model. The key step leading to the proposal of our model was that the transcription regulatory sequence (TRS) motifs were identified as the cleavage sites of nsp15, a nidoviral RNA uridylate-specific endoribonuclease (NendoU). According to this model, nsp15 regulates the synthesis of subgenomic RNAs (sgRNAs), and genomic RNAs (gRNAs) by cleaving TRSs. The expression level of nsp15 controls the relative proportions of sgRNAs and gRNAs, which in turn change the expression level of nsp15 to reach equilibrium between the CoV replication and transcription.

CONCLUSION

The replication and transcription of CoVs are regulated by a negative feedback mechanism that influences the persistence of CoVs in hosts. Our findings enrich fundamental knowledge in the field of gene expression and its regulation, and provide new clues for future studies. One important clue is that nsp15 may be an important and ideal target for the development of drugs (e.g., uridine derivatives) against CoVs.

Complete mitochondrial genome of a blue-tailed skink Plestiodon capito (Reptilia, Squamata, Scincidae) and comparison with other Scincidae lizards

Vertebrate mitochondrial genomes (mitogenomes) are valuable for studying phylogeny, evolutionary genetics and genomics. To date, however, compared to other vertebrate groups, our knowledge about the mitogenomes of skinks (the family Scincidae), even of reptile, has been relatively limited. In the present study, we determined the complete mitogenome of a blue-tailed skink Plestiodon capito for the first time, and compared it with other skinks available in GenBank. The circular genome is 17,344 bp long, showing a typical vertebrate pattern with 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and one control region (CR). The gene organization, nucleotide composition, and codon usage are similar to those from skinks previously published. Twelve out of 13 PCGs initiates with canonical start codon (ATG), while COX1 starts with GTG. The codon usage analysis revealed a preferential use of the Leu^CUN (Leu1), Pro, and Thr codons with the A/U ending. All tRNAs in P. capito were predicted to fold into typical clover-leaf secondary structure, except tRNA-Ser ^AGY. The secondary structures of 12S rRNA and 16S rRNA comprises 34 helices and 56 helices, respectively. The alignment of the Plesitodon species CRs exhibited high genetic variability and rich A + T content. Besides, variable types and numbers of tandem repeat units were also identified in the CR of Plestiodon. Phylogenetic analyses recovered P. capito as the sister species to P. tunganus; monophyly of the Scincidae is well supported. Our results will help to better understand structure and evolution of the mitochondrial DNA control region in reptiles as well as the evolutionary status of P. capito, and to lay foundation for further phylogenetic study of skinks in a mitogenomic framework.

Full-length genome sequence of segmented RNA virus from ticks was obtained using small RNA sequencing data

BACKGROUND

In 2014, a novel tick-borne virus of the Flaviviridae family was first reported in the Mogiana region of Brazil and named the Mogiana tick virus (MGTV). Thereafter, the Jingmen tick virus (JMTV), Kindia tick virus (KITV), and Guangxi tick virus (GXTV)-evolutionarily related to MGTV-were reported.

RESULTS

In the present study, we used small RNA sequencing (sRNA-seq) to detect viruses in ticks and discovered a new MGTV strain in Amblyomma testudinarium ticks collected in China's Yunnan Province in 2016. We obtained the full-length genome sequence of this MGTV strain Yunnan2016 (GenBank: MT080097, MT080098, MT080099 and MT080100) and recommended it for its inclusion in the NCBI RefSeq database for future studies on MGTV, JMTV, KITV and GXTV. Phylogenetic analysis showed that MGTV, JMTV, KITV and GXTV are monophyletic and belong to a MGTV group. Furthermore, this MGTV group of viruses may be phylogenetically related to geographical regions that were formerly part of the supercontinents Gondwana and Laurasia.

CONCLUSIONS

To the best of our knowledge, this is the first study in which 5' and 3' sRNAs were used to generate full-length genome sequences of, but not limited to, RNA viruses. We also demonstrated the feasibility of using the sRNA-seq based method for the detection of viruses in pooled two and even possible one small ticks. MGTV may preserve the characteristic of ancient RNA viruses, which can be used to study the origin and evolution of RNA viruses. In addition, MGTV can be used as novel species for studies in phylogeography.

Precise annotation of tick mitochondrial genomes reveals multiple copy number variation of short tandem repeats and one transposon-like element

Background

In the present study, we used long-PCR amplification coupled with Next-Generation Sequencing (NGS) to obtain complete mitochondrial (mt) genomes of individual ticks and unprecedently performed precise annotation of these mt genomes. We aimed to: (1) develop a simple, cost-effective and accurate method for the study of extremely high AT-content mt genomes within an individual animal (e.g. Dermacentor silvarum) containing miniscule DNA; (2) provide a high-quality reference genome for D. silvarum with precise annotation and also for future studies of other tick mt genomes; and (3) detect and analyze mt DNA variation within an individual tick.

Results

These annotations were confirmed by the PacBio full-length transcriptome data to cover both entire strands of the mitochondrial genomes without any gaps or overlaps. Moreover, two new and important findings were reported for the first time, contributing fundamental knowledge to mt biology. The first was the discovery of a transposon-like element that may eventually reveal much about mechanisms of gene rearrangements in mt genomes. Another finding was that Copy Number Variation (CNV) of Short Tandem Repeats (STRs) account for mitochondrial sequence diversity (heterogeneity) within an individual tick, insect, mouse or human, whereas SNPs were not detected. The CNV of STRs in the protein-coding genes resulted in frameshift mutations in the proteins, which can cause deleterious effects. Mitochondria containing these deleterious STR mutations accumulate in cells and can produce deleterious proteins.

Conclusions

We proposed that the accumulation of CNV of STRs in mitochondria may cause aging or diseases. Future tests of the CNV of STRs hypothesis help to ultimately reveal the genetic basis of mitochondrial DNA variation and its consequences (e.g., aging and diseases) in animals. Our study will lead to the reconsideration of the importance of STRs and a unified study of CNV of STRs with longer and shorter repeat units (particularly polynucleotides) in both nuclear and mt genomes.

Precise annotation of human, chimpanzee, rhesus macaque and mouse mitochondrial genomes leads to insight into mitochondrial transcription in mammals

In the present study, we applied our 'precise annotation' to the mitochondrial (mt) genomes of human, chimpanzee, rhesus macaque and mouse using 5' and 3' end small RNAs. Our new annotations updated previous annotations. In particular, our new annotations led to two important novel findings: (1) the identification of five Conserved Sequence Blocks (CSB1, CSB2, CSB3, LSP and HSP) in the control regions; and (2) the annotation of Transcription Initiation and novel Transcription Termination Sites. Based on these annotations, we proposed a novel model of mt transcription which can account for the mt transcription and its regulation in mammals. According to our model, Transcription Termination Sites function as switches to regulate the production of short, long primary transcripts and uninterrupted transcription, rather than simply terminate the mt transcription. Moreover, the expression levels of mitochondrial transcription termination factors control the proportions of rRNAs, mRNAs and lncRNAs in total mt RNA. Our findings point to the existence of many other, as yet unidentified, Transcription Termination Sites in mammals.

Using high-resolution annotation of insect mitochondrial DNA to decipher tandem repeats in the control region

In this study, we used a small RNA sequencing (sRNA-seq) based method to annotate the mitochondrial genome of the insect Erthesina fullo Thunberg at 1 bp resolution. The high-resolution annotations cover both entire strands of the mitochondrial genome without any gaps or overlaps. Most of the new annotations were consistent with the previous annotations which had been obtained using PacBio full-length transcripts. Two important findings were that animals transcribe both entire strands of mitochondrial genomes and the tandem repeats in the control region of the E. fullo mitochondrial genome contains the repeated Transcription Initiation Sites (TISs) of the heavy strand. In addition, we found that the copy numbers of tandem repeats showed a great diversity within an individual, suggesting that mitochondrial DNA recombination occurs in an individual. In conclusion, the sRNA-seq based method uses 5' and 3' end small RNAs to annotate nuclear non-coding and mitochondrial genes at 1 bp resolution, and can be used to identify new steady RNAs, particularly long non-coding RNAs (lncRNAs). The high-resolution annotations of mitochondrial genomes can also be used to study the molecular phylogenetics and evolution of animals or to investigate mitochondrial gene transcription, RNA processing, RNA maturation and several other related topics. The complete mitochondrial genome sequence of E. fullo with the new annotations using the sRNA-seq based method is available at the NCBI GenBank database under the accession number MK374364. We publish our theories, methods, the high quality sRNA-seq and RNA-seq data (SRA: SRP174926) for extensive use.