Juergen Brosius

Biases in small RNA deep sequencing data

High-throughput RNA sequencing (RNA-seq) is considered a powerful tool for novel gene discovery and fine-tuned transcriptional profiling. The digital nature of RNA-seq is also believed to simplify meta-analysis and to reduce background noise associated with hybridization-based approaches. The development of multiplex sequencing enables efficient and economic parallel analysis of gene expression. In addition, RNA-seq is of particular value when low RNA expression or modest changes between samples are monitored. However, recent data uncovered severe bias in the sequencing of small non-protein coding RNA (small RNA-seq or sRNA-seq), such that the expression levels of some RNAs appeared to be artificially enhanced and others diminished or even undetectable. The use of different adapters and barcodes during ligation as well as complex RNA structures and modifications drastically influence cDNA synthesis efficacies and exemplify sources of bias in deep sequencing. In addition, variable specific RNA G/C-content is associated with unequal polymerase chain reaction amplification efficiencies. Given the central importance of RNA-seq to molecular biology and personalized medicine, we review recent findings that challenge small non-protein coding RNA-seq data and suggest approaches and precautions to overcome or minimize bias.

A global view of the nonprotein-coding transcriptome in Plasmodium falciparum

Nonprotein-coding RNAs (npcRNAs) represent an important class of regulatory molecules that act in many cellular pathways. Here, we describe the experimental identification and validation of the small npcRNA transcriptome of the human malaria parasite Plasmodium falciparum. We identified 630 novel npcRNA candidates. Based on sequence and structural motifs, 43 of them belong to the C/D and H/ACA-box subclasses of small nucleolar RNAs (snoRNAs) and small Cajal body-specific RNAs (scaRNAs). We further observed the exonization of a functional H/ACA snoRNA gene, which might contribute to the regulation of ribosomal protein L7a gene expression. Some of the small npcRNA candidates are from telomeric and subtelomeric repetitive regions, suggesting their potential involvement in maintaining telomeric integrity and subtelomeric gene silencing. We also detected 328 cis-encoded antisense npcRNAs (asRNAs) complementary to P. falciparum protein-coding genes of a wide range of biochemical pathways, including determinants of virulence and pathology. All cis-encoded asRNA genes tested exhibit lifecycle-specific expression profiles. For all but one of the respective sense–antisense pairs, we deduced concordant patterns of expression. Our findings have important implications for a better understanding of gene regulatory mechanisms in P. falciparum, revealing an extended and sophisticated npcRNA network that may control the expression of housekeeping genes and virulence factors.

Experimental identification and characterization of 97 novel npcRNA candidates in Salmonella enterica serovar Typhi

We experimentally identified and characterized 97 novel, non-protein-coding RNA candidates (npcRNAs) from the human pathogen Salmonella enterica serovar Typhi (hereafter referred to as S. typhi). Three were specific to S. typhi, 22 were restricted to Salmonella species and 33 were differentially expressed during S. typhi growth. We also identified Salmonella Pathogenicity Island-derived npcRNAs that might be involved in regulatory mechanisms of virulence, antibiotic resistance and pathogenic specificity of S. typhi. An in-depth characterization of S. typhi StyR-3 npcRNA showed that it specifically interacts with RamR, the transcriptional repressor of the ramA gene, which is involved in the multidrug resistance (MDR) of Salmonella. StyR-3 interfered with RamR–DNA binding activity and thus potentially plays a role in regulating ramA gene expression, resulting in the MDR phenotype. Our study also revealed a large number of cis-encoded antisense npcRNA candidates, supporting previous observations of global sense–antisense regulatory networks in bacteria. Finally, at least six of the npcRNA candidates interacted with the S. typhi Hfq protein, supporting an important role of Hfq in npcRNA networks. This study points to novel functional npcRNA candidates potentially involved in various regulatory roles including the pathogenicity of S. typhi.

The genome of a Mesozoic paleovirus reveals the evolution of hepatitis B viruses

Paleovirology involves the identification of ancient endogenous viral elements within eukaryotic genomes. The evolutionary origins of the reverse-transcribing hepatitis B viruses, however, remain elusive, due to the small number of endogenized sequences present in host genomes. Here we report a comprehensively dated genomic record of hepatitis B virus endogenizations that spans bird evolution from >82 to <12.1 million years ago. The oldest virus relic extends over a 99% complete hepatitis B virus genome sequence and constitutes the first discovery of a Mesozoic paleovirus genome. We show that Hepadnaviridae are >63 million years older than previously known and provide direct evidence for coexistence of hepatitis B viruses and birds during the Mesozoic and Cenozoic Eras. Finally, phylogenetic analyses and distribution of hepatitis B virus relics suggest that birds potentially are the ancestral hosts of Hepadnaviridae and mammalian hepatitis B viruses probably emerged after a bird-mammal host switch. Our study reveals previously undiscovered and multi-faceted insights into prehistoric hepatitis B virus evolution and provides valuable resources for future studies, such as in-vitro resurrection of Mesozoic hepadnaviruses.

Multiple Lineages of Ancient CR1 Retroposons Shaped the Early Genome Evolution of Amniotes

Chicken repeat 1 (CR1) retroposons are long interspersed elements (LINEs) that are ubiquitous within amniote genomes and constitute the most abundant family of transposed elements in birds, crocodilians, turtles, and snakes. They are also present in mammalian genomes, where they reside as numerous relics of ancient retroposition events. Yet, despite their relevance for understanding amniote genome evolution, the diversity and evolution of CR1 elements has never been studied on an amniote-wide level. We reconstruct the temporal and quantitative activity of CR1 subfamilies via presence/absence analyses across crocodilian phylogeny and comparative analyses of 12 crocodilian genomes, revealing relative genomic stasis of retroposition during genome evolution of extant Crocodylia. Our large-scale phylogenetic analysis of amniote CR1 subfamilies suggests the presence of at least seven ancient CR1 lineages in the amniote ancestor; and amniote-wide analyses of CR1 successions and quantities reveal differential retention (presence of ancient relics or recent activity) of these CR1 lineages across amniote genome evolution. Interestingly, birds and lepidosaurs retained the fewest ancient CR1 lineages among amniotes and also exhibit smaller genome sizes. Our study is the first to analyze CR1 evolution in a genome-wide and amniote-wide context and the data strongly suggest that the ancestral amniote genome contained myriad CR1 elements from multiple ancient lineages, and remnants of these are still detectable in the relatively stable genomes of crocodilians and turtles. Early mammalian genome evolution was thus characterized by a drastic shift from CR1 prevalence to dominance and hyperactivity of L2 LINEs in monotremes and L1 LINEs in therians.

Genomic impact of eukaryotic transposable elements

The third international conference on the genomic impact of eukaryotic transposable elements (TEs) was held 24 to 28 February 2012 at the Asilomar Conference Center, Pacific Grove, CA, USA. Sponsored in part by the National Institutes of Health grant 5 P41 LM006252, the goal of the conference was to bring together researchers from around the world who study the impact and mechanisms of TEs using multiple computational and experimental approaches. The meeting drew close to 170 attendees and included invited floor presentations on the biology of TEs and their genomic impact, as well as numerous talks contributed by young scientists. The workshop talks were devoted to computational analysis of TEs with additional time for discussion of unresolved issues. Also, there was ample opportunity for poster presentations and informal evening discussions. The success of the meeting reflects the important role of Repbase in comparative genomic studies, and emphasizes the need for close interactions between experimental and computational biologists in the years to come.

Isolation and Posttranscriptional Modification Analysis of Native BC1 RNA from Mouse Brain

We present a simple and general affinity method, based on size fractionation and nucleic acid complementarity, to isolate sufficient amounts of native RNA molecules for further physicochemical studies, such as modification state of nucleotides. In the case presented here, we purified four micrograms of dendritic neuronal BC1 RNA from 130 grams of mouse brain (initially yielding a total of 200 mg RNA). Directly combined liquid chromatography-electrospray ionization mass spectrometry (LC/MS) analysis revealed no base or sugar backbone modifications in native BC1 RNA, despite earlier indications that C-54 could be methylated in vitro (Cm5, position 54).

Alternative processing as evolutionary mechanism for the origin of novel nonprotein coding RNAs.

The evolution of new genes can ensue through either gene duplication and the neofunctionalization of one of the copies or the formation of a de novo gene from hitherto nonfunctional, neutrally evolving intergenic or intronic genomic sequences. Only very rarely are entire genes created de novo. Mostly, nonfunctional sequences are coopted as novel parts of existing genes, such as in the process of exonization whereby introns become exons through changes in splicing. Here, we report a case in which a novel nonprotein coding RNA evolved by intron-sequence recruitment into its structure. cDNAs derived from rat brain small RNAs, revealed a novel small nucleolar RNA (snoRNA) originating from one of the Snord115 copies in the rat Prader–Willi syndrome locus. We suggest that a single-point substitution in the Snord115 region led to the expression of a longer snoRNA variant, designated as L-Snord115. Cell culture and footprinting experiments confirmed that a single nucleotide substitution at Snord115 position 67 destabilized the kink-turn motif within the canonical snoRNA, while distal intronic sequences provided an alternate D-box region. The exapted sequence displays putative base pairing to 28S rRNA and mRNA targets.

De-novo emergence and template switching of SINE retroposons during the early evolution of passerine birds

Passeriformes (“perching birds” or passerines) make up more than half of all extant bird species. Here, we resolve their deep phylogenetic relationships using presence/absence patterns of short interspersed elements (SINEs), a group of retroposons which is abundant in mammalian genomes but considered largely inactive in avian genomes. The resultant retroposon-based phylogeny provides a powerful and independent corroboration of previous indications derived from sequence-based analyses. Notably, SINE activity began in the common ancestor of Eupasseres (passerines excl. the New Zealand wrens Acanthisittidae) and ceased before the rapid diversification of oscine passerines (songbirds). Furthermore, we find evidence for very recent SINE activity within suboscine passerines, following the emergence of a SINE via acquisition of a different tRNA head as we suggest through template switching. We propose that the early evolution of passerines was unusual among birds in that it was accompanied by activity of SINEs. Their genomic and transcriptomic impact warrants further study in the light of the massive diversification of passerines.

Characterization of novel non-protein-coding RNA from V. cholerae O1 El Tor

Abstract Introduction A plethora of bacterial small RNA (npcRNA) has progressively gained momentum in revising our apprehension on their engagement as key player in modulating general stress responses and bacterial virulence. Vibrio cholerae serogroup O1, biotype El Tor causes cholera, which is an acute dehydrating, watery diarrhoeal disease, with epidemic and pandemic potential. We have identified a total of 224 novel npcRNA candidates from V. cholerae El Tor by experimental RNomics strategy. Objective: To characterize the expression profile and functionality of identified novel npcRNA from V. cholerae O1 El Tor. Methods Differential expression of selected npcRNA candidates was studied using Northern blot analysis under various stress conditions in the background wild type strain as well as Hfq protein knockout strain. Interaction of the selected npcRNA candidates with Hfq protein was also evaluated. Results & Discussion Two intergenic npcRNA genes were incorporated as targets in a multiplex PCR that could be a crucial tool in molecular epidemiological study of V. cholerae . Differential expression of 9 selected npcRNA candidates under a series of stress conditions was demonstrated. Interestingly, one of the candidate, Vc_npcR_3853 was demonstrated to have the potential in negatively regulating the expression of VC0092 and VC0304, which encodeLexA repressor and guanosine-5′-triphosphate, 3′-diphosphate pyrophosphatase, respectively. Conclusion Collectively, differential expression of those npcRNA candidates under wild type and Hfq knockout background as well as their interaction with Hfq protein set the stage towards discerning the functional role of these novel npcRNAs in orchestrating pathophysiology of V. cholerae El Tor.