RNA-DNA sequence differences in Saccharomyces cerevisiae

Metabolic stress in space: ROS-induced mutations in mice hint at a new path to cancer

Long-duration spaceflight beyond Earth's magnetosphere poses serious health risks, including muscle atrophy, bone loss, liver and kidney damage, and the Spaceflight-Associated Neuro-ocular Syndrome (SANS). RNA-seq of mice aboard the International Space Station (ISS) for 37 days revealed extraordinary hypermutation in tissue-specific genes, with guanine base conversion predominating, potentially contributing to spaceflight-associated health risks. Our results suggest that the genome-wide accelerated mutation that we measured, seemingly independent of radiation dose, was induced by oxidative damage from higher atmospheric carbon dioxide (CO2) levels and increased reactive oxygen species (ROS) on the ISS. This accelerated mutation, faster via RNA transcription than replication and more numerous than by radiation alone, unveils novel hotspots in the mammalian proteome. Notably, these hotspots correlate with commonly mutated genes across various human cancers, highlighting the ISS as a crucial platform for studying accelerated mutation, genome instability, and the induction of disease-causing mutations in model organisms. Our results suggest that metabolic processes can contribute to somatic mutation, and thus may play a role in the development of cancer. A metabolic link to genetic instability potentially has far-reaching implications for various diseases, with implications for human health on Earth and in space.

RNA-DNA differences in variant calls from cattle tissues result in erroneous eQTLs

BACKGROUND

Association testing between molecular phenotypes and genomic variants can help to understand how genotype affects phenotype. RNA sequencing provides access to molecular phenotypes such as gene expression and alternative splicing while DNA sequencing or microarray genotyping are the prevailing options to obtain genomic variants.

RESULTS

We genotype variants for 74 male Braunvieh cattle from both DNA (~ 13-fold coverage) and deep total RNA sequencing from testis, vas deferens, and epididymis tissue (~ 250 million reads per tissue). We show that RNA sequencing can be used to identify approximately 40% of variants (7-10 million) called from DNA sequencing, with over 80% precision. Within highly expressed coding regions, over 92% of expected variants were called with nearly 98% precision. Allele-specific expression and putative post-transcriptional modifications negatively impact variant genotyping accuracy from RNA sequencing and contribute to RNA-DNA differences. Variants called from RNA sequencing detect roughly 75% of eGenes identified using variants called from DNA sequencing, demonstrating a nearly 2-fold enrichment of eQTL variants. We observe a moderate-to-strong correlation in nominal association p-values (Spearman ρ2 ~ 0.6), although only 9% of eGenes have the same top associated variant.

CONCLUSIONS

We find hundreds of thousands of RNA-DNA differences in variants called from RNA and DNA sequencing on the same individuals. We identify several highly significant eQTL when using RNA sequencing variant genotypes which are not found with DNA sequencing variant genotypes, suggesting that using RNA sequencing variant genotypes for association testing results in an increased number of false positives. Our findings demonstrate that caution must be exercised beyond filtering for variant quality or imputation accuracy when analysing or imputing variants called from RNA sequencing.

Repurposing RNA sequencing for discovery of RNA modifications in clinical cohorts

The study of RNA modifications in large clinical cohorts can reveal relationships between the epitranscriptome and human diseases, although this is especially challenging. We developed ModTect (https://github.com/ktan8/ModTect), a statistical framework to identify RNA modifications de novo by standard RNA-sequencing with deletion and mis-incorporation signals. We show that ModTect can identify both known (N 1-methyladenosine) and previously unknown types of mRNA modifications (N 2,N 2-dimethylguanosine) at nucleotide-resolution. Applying ModTect to 11,371 patient samples and 934 cell lines across 33 cancer types, we show that the epitranscriptome was dysregulated in patients across multiple cancer types and was additionally associated with cancer progression and survival outcomes. Some types of RNA modification were also more disrupted than others in patients with cancer. Moreover, RNA modifications contribute to multiple types of RNA-DNA sequence differences, which unexpectedly escape detection by Sanger sequencing. ModTect can thus be used to discover associations between RNA modifications and clinical outcomes in patient cohorts.

RNA abasic sites in yeast and human cells

RNA abasic sites and the mechanisms involved in their regulation are mostly unknown; in contrast, DNA abasic sites are well-studied. We found surprisingly that, in yeast and human cells, RNA abasic sites are prevalent. When a base is lost from RNA, the remaining ribose is found as a closed-ring or an open-ring sugar with a reactive C1' aldehyde group. Using primary amine-based reagents that react with the aldehyde group, we uncovered evidence for abasic sites in nascent RNA, messenger RNA, and ribosomal RNA from yeast and human cells. Mass spectroscopic analysis confirmed the presence of RNA abasic sites. The RNA abasic sites were found to be coupled to R-loops. We show that human methylpurine DNA glycosylase cleaves N-glycosidic bonds on RNA and that human apurinic/apyrimidinic endonuclease 1 incises RNA abasic sites in RNA-DNA hybrids. Our results reveal that, in yeast and human cells, there are RNA abasic sites, and we identify a glycosylase that generates these sites and an AP endonuclease that processes them.

Gene fragmentation and RNA editing without borders: eccentric mitochondrial genomes of diplonemids

Diplonemids are highly abundant heterotrophic marine protists. Previous studies showed that their strikingly bloated mitochondrial genome is unique because of systematic gene fragmentation and manifold RNA editing. Here we report a comparative study of mitochondrial genome architecture, gene structure and RNA editing of six recently isolated, phylogenetically diverse diplonemid species. Mitochondrial gene fragmentation and modes of RNA editing, which include cytidine-to-uridine (C-to-U) and adenosine-to-inosine (A-to-I) substitutions and 3' uridine additions (U-appendage), are conserved across diplonemids. Yet as we show here, all these features have been pushed to their extremes in the Hemistasiidae lineage. For example, Namystynia karyoxenos has its genes fragmented into more than twice as many modules than other diplonemids, with modules as short as four nucleotides. Furthermore, we detected in this group multiple A-appendage and guanosine-to-adenosine (G-to-A) substitution editing events not observed before in diplonemids and found very rarely elsewhere. With >1,000 sites, C-to-U and A-to-I editing in Namystynia is nearly 10 times more frequent than in other diplonemids. The editing density of 12% in coding regions makes Namystynia's the most extensively edited transcriptome described so far. Diplonemid mitochondrial genome architecture, gene structure and post-transcriptional processes display such high complexity that they challenge all other currently known systems.

Systematic Nucleotide Exchange Analysis of ESTs From the Human Cancer Genome Project Report: Origins of 347 Unknown ESTs Indicate Putative Transcription of Non-Coding Genomic Regions

Expressed sequence tags (ESTs) provide an imprint of cellular RNA diversity irrespectively of sequence homology with template genomes. NCBI databases include many unknown RNAs from various normal and cancer cells. These are usually ignored assuming sequencing artefacts or contamination due to their lack of sequence homology with template DNA. Here, we report genomic origins of 347 ESTs previously assumed artefacts/unknown, from the FAPESP/LICR Human Cancer Genome Project. EST template detection uses systematic nucleotide exchange analyses called swinger transformations. Systematic nucleotide exchanges replace systematically particular nucleotides with different nucleotides. Among 347 unknown ESTs, 51 ESTs match mitogenome transcription, 17 and 2 ESTs are from nuclear chromosome non-coding regions, and uncharacterized nuclear genes. Identified ESTs mapped on 205 protein-coding genes, 10 genes had swinger RNAs in several biosamples. Whole cell transcriptome searches for 17 ESTs mapping on non-coding regions confirmed their transcription. The 10 swinger-transcribed genes identified more than once associate with cancer induction and progression, suggesting swinger transformation occurs mainly in highly transcribed genes. Swinger transformation is a unique method to identify noncanonical RNAs obtained from NGS, which identifies putative ncRNA transcribed regions. Results suggest that swinger transcription occurs in highly active genes in normal and genetically unstable cancer cells.

Transcripts with systematic nucleotide deletion of 1-12 nucleotide in human mitochondrion suggest potential non-canonical transcription

Raw transcriptomic data contain numerous RNA reads whose homology with template DNA doesn't match canonical transcription. Transcriptome analyses usually ignore such noncanonical RNA reads. Here, analyses search for noncanonical mitochondrial RNAs systematically deleting 1 to 12 nucleotides after each transcribed nucleotide triplet, producing deletion-RNAs (delRNAs). We detected delRNAs in the human whole cell and purified mitochondrial transcriptomes, and in Genbank's human EST database corresponding to systematic deletions of 1 to 12 nucleotides after each transcribed trinucleotide. DelRNAs detected in both transcriptomes mapped along with 55.63% of the EST delRNAs. A bias exists for delRNAs covering identical mitogenomic regions in both transcriptomic and EST datasets. Among 227 delRNAs detected in these 3 datasets, 81.1% and 8.4% of delRNAs were mapped on mitochondrial coding and hypervariable region 2 of dloop. Del-transcription analyses of GenBank's EST database confirm observations from whole cell and purified mitochondrial transcriptomes, eliminating the possibility that detected delRNAs are false positives matches, cytosolic DNA/RNA nuclear contamination or sequencing artefacts. These detected delRNAs are enriched in frameshift-inducing homopolymers and are poor in frameshift-preventing circular code codons (a set of 20 codons which regulate reading frame detection, over- and underrepresented in coding and other frames of genes, respectively) suggesting a motif-based regulation of non-canonical transcription. These findings show that rare non-canonical transcripts exist. Such non canonical del-transcription does increases mitochondrial coding potential and non-coding regulation of intracellular mechanisms, and could explain the dark DNA conundrum.

A-to-I mRNA editing in fungi: occurrence, function, and evolution

A-to-I RNA editing is an important post-transcriptional modification that converts adenosine (A) to inosine (I) in RNA molecules via hydrolytic deamination. Although editing of mRNAs catalyzed by adenosine deaminases acting on RNA (ADARs) is an evolutionarily conserved mechanism in metazoans, organisms outside the animal kingdom lacking ADAR orthologs were thought to lack A-to-I mRNA editing. However, recent discoveries of genome-wide A-to-I mRNA editing during the sexual stage of the wheat scab fungus Fusarium graminearum, model filamentous fungus Neurospora crassa, Sordaria macrospora, and an early diverging filamentous ascomycete Pyronema confluens indicated that A-to-I mRNA editing is likely an evolutionarily conserved feature in filamentous ascomycetes. More importantly, A-to-I mRNA editing has been demonstrated to play crucial roles in different sexual developmental processes and display distinct tissue- or development-specific regulation. Contrary to that in animals, the majority of fungal RNA editing events are non-synonymous editing, which were shown to be generally advantageous and favored by positive selection. Many non-synonymous editing sites are conserved among different fungi and have potential functional and evolutionary importance. Here, we review the recent findings about the occurrence, regulation, function, and evolution of A-to-I mRNA editing in fungi.

SMARTcleaner: identify and clean off-target signals in SMART ChIP-seq analysis

Background

Noises and artifacts may arise in several steps of the next-generation sequencing (NGS) process. Recently, an NGS library preparation method called SMART, or Switching Mechanism At the 5′ end of the RNA Transcript, is introduced to prepare ChIP-seq (chromatin immunoprecipitation and deep sequencing) libraries from small amount of DNA material, using the DNA SMART ChIP-seq Kit. The protocol adds Ts to the 3′ end of DNA templates, which is subsequently recognized and used by SMART poly(dA) primers for reverse transcription and then addition of PCR primers and sequencing adapters. The poly(dA) primers, however, can anneal to poly(T) sequences in a genome and amplify DNA fragments that are not enriched in the immunoprecipitated DNA templates. This off-target amplification results in false signals in the ChIP-seq data.

Results

Here, we show that the off-target ChIP-seq reads derived from false amplification of poly(T/A) genomic sequences have unique and strand-specific features. Accordingly, we develop a tool (called “SMARTcleaner”) that can exploit these features to remove SMART ChIP-seq artifacts. Application of SMARTcleaner to several SMART ChIP-seq datasets demonstrates that it can remove reads from off-target amplification effectively, leading to significantly improved ChIP-seq peaks and results.

Conclusions

SMARTcleaner could identify and clean the false signals in SMART-based ChIP-seq libraries, leading to improvement in peak calling, and downstream data analysis and interpretation.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2577-4) contains supplementary material, which is available to authorized users.

RNA editing in bacteria: occurrence, regulation and significance

DNA harbors the blueprint for life. However, the instructions stored in the DNA could be altered at the RNA level before they are executed. One of these processes is RNA editing, which was shown to modify RNA sequences in many organisms. The most abundant modification is the deamination of adenosine (A) into inosine (I). In turn, inosine can be identified as a guanosine (G) by the ribosome and other cellular machineries such as reverse transcriptase. In multicellular organisms, enzymes from the ADAR (adenosine deaminase acting on RNA) family mediate RNA editing in mRNA, whereas enzymes from the ADAT family mediate A-to-I editing on tRNAs. In bacteria however, until recently, only one editing site was described, in tRNA^Arg, but never in mRNA. The tRNA site was shown to be modified by tadA (tRNA specific adenosine deaminase) which is believed to be the ancestral enzyme for the RNA editing family of enzymes. In our recent work, we have shown for the first time, editing on multiple sites in bacterial mRNAs and identified tadA as the enzyme responsible for this editing activity. Focusing on one of the identified targets - the self-killing toxin hokB, we found that editing is physiologically regulated and that it increases protein activity. Here we discuss possible modes of regulation on hokB editing, potential roles of RNA editing in bacteria, possible implications, and future research directions.

Immunoprecipitation of RNA:DNA Hybrids from Budding Yeast

During transcription, the nascent transcript behind an elongating RNA polymerase (RNAP) can invade the DNA duplex and hybridize with the complementary DNA template strand, generating a three-stranded "R-loop" structure, composed of an RNA:DNA duplex and an unpaired non-template DNA strand. R-loops can be strongly associated with actively transcribed loci by all RNAPs including the mitochondrial RNA polymerase (mtRNAP). In this chapter, we describe two protocols for the detection of RNA:DNA hybrids in living budding yeast cells, one that uses conventional chromatin immunoprecipitation (ChIP-qPCR) and one that uses DNA:RNA immunoprecipitation (DRIP-qPCR). Both protocols make use of the S9.6 antibody, which is believed to recognize the intermediate A/B helical RNA:DNA duplex conformation, with no sequence specificity.

RNA editing in bacteria recodes multiple proteins and regulates an evolutionarily conserved toxin-antitoxin system

Adenosine (A) to inosine (I) RNA editing is widespread in eukaryotes. In prokaryotes, however, A-to-I RNA editing was only reported to occur in tRNAs but not in protein-coding genes. By comparing DNA and RNA sequences of Escherichia coli, we show for the first time that A-to-I editing occurs also in prokaryotic mRNAs and has the potential to affect the translated proteins and cell physiology. We found 15 novel A-to-I editing events, of which 12 occurred within known protein-coding genes where they always recode a tyrosine (TAC) into a cysteine (TGC) codon. Furthermore, we identified the tRNA-specific adenosine deaminase (tadA) as the editing enzyme of all these editing sites, thus making it the first identified RNA editing enzyme that modifies both tRNAs and mRNAs. Interestingly, several of the editing targets are self-killing toxins that belong to evolutionarily conserved toxin-antitoxin pairs. We focused on hokB, a toxin that confers antibiotic tolerance by growth inhibition, as it demonstrated the highest level of such mRNA editing. We identified a correlated mutation pattern between the edited and a DNA hard-coded Cys residue positions in the toxin and demonstrated that RNA editing occurs in hokB in two additional bacterial species. Thus, not only the toxin is evolutionarily conserved but also the editing itself within the toxin is. Finally, we found that RNA editing in hokB increases as a function of cell density and enhances its toxicity. Our work thus demonstrates the occurrence, regulation, and functional consequences of RNA editing in bacteria.

Reviewing evidence for systematic transcriptional deletions, nucleotide exchanges, and expanded codons, and peptide clusters in human mitochondria

Polymerization sometimes transforms sequences by (a) systematic deletions of mono-, dinucleotides after trinucleotides, or (b) 23 systematic nucleotide exchanges (9 symmetric, X<>Y, e.g. G<>T, 14 asymmetric, X > Y > Z > X, e.g. A > G > T > A), producing del- and swinger RNAs. Some peptides correspond to del- and swinger RNA translations, also according to tetracodons, codons expanded by a silent nucleotide. Here new analyzes assume different proteolytic patterns, partially alleviating false negative peptide detection biases, expanding noncanonical mitoproteome profiles. Mito-genomic, -transcriptomic and -proteomic evidence for noncanonical transcriptions and translations are reviewed and clusters of del- and swinger peptides (also along tetracodons) are described. Noncanonical peptide clusters indicate regulated expression of cryptically encoded mitochondrial protein coding genes. These candidate noncanonical proteins don't resemble known proteins.