Peregrine: A rapid and unbiased method to produce strand-specific RNA-Seq libraries from small quantities of starting material

Divergent Heat Stress Responses in Bactrocera tryoni and Ceratitis capitata

Invasive Tephritid fruit flies rank among the most destructive agricultural and horticultural pests worldwide. Heat treatment is commonly employed as a post-harvest method to exterminate fruit flies in fruits or vegetables. These pest species exhibit distinct tolerance to heat treatments, suggesting that the molecular pathways affected by heat may differ among species. In this study, the Queensland fruit fly (Qfly), Bactrocera tryoni, was utilised as a model investigate its molecular response to heat stress through heat bioassays. RNA samples from flies before and after heat treatment were extracted and sequenced to identify genes with significant changes in expression. These findings were compared to another serious Tephritid fruit fly species, the Mediterranean fruit fly (Medfly), Ceratitis capitata, under similar heat treatment conditions. The analysis reveals only three common genes: heat shock protein 70 (HSP70), HSP68, and 14-3-3 zeta protein. However, despite these shared genes, their expression patterns differ between Qfly and Medfly. This suggests that these genes might play different roles in the heat responses of each species and could be regulated differently. This study presents the first evidence of differing molecular responses to heat between Qfly and Medfly, potentially linked to their varied origins, habitats, and genetic backgrounds. These findings offer new insights into Tephritid fruit fly responses to heat at the molecular level, which may help refine post-harvest strategies to control these pests in the future.

Cost-Effective Whole Transcriptome Sequencing Landscape and Diagnostic Potential Biomarkers in Active Tuberculosis

Tuberculosis (TB) is a prevalent and severe infectious disease that poses a significant threat to human health. However, it is frequently disregarded as there are not enough quick and accurate ways to diagnose tuberculosis. Here, we develop a strategy for tuberculosis detection to address the challenges, including an experimental strategy, namely, Double Adapter Directional Capture sequencing (DADCSeq), an easily operated and low-cost whole transcriptome sequencing method, and a computational method to identify hub differentially expressed genes as well as the diagnosis of TB based on whole transcriptome data using DADCSeq on peripheral blood mononuclear cells (PBMCs) from active TB and latent TB or healthy control. Applying our approach to create a robust and stable TB multi-mRNA risk probability model (TBMMRP) that can accurately distinguish active and latent TB patients, including active TB and healthy controls in clinical cohorts, this diagnostic biomarker was successfully validated by several independent cross-platform cohorts with favorable performance in differentiating active TB from latent TB or active TB from healthy controls and further demonstrated superior or similar diagnostic accuracy compared to previous diagnostic markers. Overall, we develop a low-cost and effective strategy for tuberculosis diagnosis; as the clinical cohort increases, we can expand to different disease kinds and learn new features through our disease diagnosis strategy.

Untacking small RNA profiling and RNA fragment footprinting: Approaches and challenges in library construction

Small RNAs (sRNAs) with sizes ranging from 15 to 50 nucleotides (nt) are critical regulators of gene expression control. Prior studies have shown that sRNAs are involved in a broad range of biological processes, such as organ development, tumorigenesis, and epigenomic regulation; however, emerging evidence unveils a hidden layer of diversity and complexity of endogenously encoded sRNAs profile in eukaryotic organisms, including novel types of sRNAs and the previously unknown post-transcriptional RNA modifications. This underscores the importance for accurate, unbiased detection of sRNAs in various cellular contexts. A multitude of high-throughput methods based on next-generation sequencing (NGS) are developed to decipher the sRNA expression and their modifications. Nonetheless, distinct from mRNA sequencing, the data from sRNA sequencing suffer frequent inconsistencies and high variations emanating from the adapter contaminations and RNA modifications, which overall skew the sRNA libraries. Here, we summarize the sRNA-sequencing approaches, and discuss the considerations and challenges for the strategies and methods of sRNA library construction. The pros and cons of sRNA sequencing have significant implications for implementing RNA fragment footprinting approaches, including CLIP-seq and Ribo-seq. We envision that this review can inspire novel improvements in small RNA sequencing and RNA fragment footprinting in future. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA Processing > Processing of Small RNAs Regulatory RNAs/RNAi/Riboswitches > Biogenesis of Effector Small RNAs.

Non-additive gene interactions underpin molecular and phenotypic responses in honey bee larvae exposed to imidacloprid and thymol

Understanding the cumulative risk of chemical mixtures at environmentally realistic concentrations is a key challenge in honey bee ecotoxicology. Ecotoxicogenomics, including transcriptomics, measures responses in individual organisms at the molecular level which can provide insights into the mechanisms underlying phenotypic responses induced by one or more stressors and link impacts on individuals to populations. Here, fifth instar honey bee larvae were sampled from a previously reported field experiment exploring the phenotypic impacts of environmentally realistic chronic exposures of the pesticide imidacloprid (5 μg.kg^-1 for six weeks) and the acaricide thymol (250 g.kg^-1 applied via Apiguard gel in-hive for four weeks), both separately and in combination. RNA-seq was used to discover individual and interactive chemical effects on larval gene expression and to uncover molecular mechanisms linked to reported adult and colony phenotypes. The separate and combined treatments had distinct gene expression profiles which represented differentially affected signaling and metabolic pathways. The molecular signature of the mixture was characterised by additive interactions in canonical stress responses associated with oxidative stress and detoxification, and non-additive interactions in secondary responses including developmental, neurological, and immune pathways. Novel emergent impacts on eye development genes correlated with long-term defects in visual learning performance as adults. This is consistent with these chemicals working through independent modes of action that combine to impact common downstream pathways, and highlights the importance of establishing mechanistic links between molecular and phenotypic responses when predicting effects of chemical mixtures on ecologically relevant population outcomes.

Molecular Response of the Mediterranean Fruit Fly (Diptera: Tephritidae) to Heat

Tephritid fruit flies are highly successful invaders and some-such as the Mediterranean fruit fly, Ceratitis capitata (Wiedemann)-are able to adapt to a large range of crops. Biosecurity controls require that shipments of produce are ensured to be pest-free, which is increasingly difficult due to the ban of key pesticides. Instead, stress-based strategies including controlled atmosphere, temperature, and irradiation can be used to eradicate flies inside products. However, unlike pesticide science, we do not yet have a robust scientific approach to measure cost-effectively whether a sufficiently lethal stress has been delivered and understand what this stress does to the biology of the pest. The latter is crucial as it would enable a combination of stresses targeting multiple molecular pathways and thus allow for lower doses of each to achieve higher lethality and reduce the development of resistance. Using heat as an example, this is the first study investigating the molecular stress response to heat in Tephritidae. Using a novel setup delivering measured doses of heat on C. capitata larvae and a high-density 11 timepoint gene expression experiment, we identified key components of lethal heat-stress response. While unraveling the complete molecular mechanism of fruit fly response to lethal stress would be a long-term project, this work curates and develops 31 potential biomarkers to assess whether sufficient lethal stress has been delivered. Further, as these protocols are straightforward and less expensive than other-omic approaches, our studies and approach will assist other researchers working on stress response.

Oxidative stress delays development and alters gene expression in the agricultural pest moth, Helicoverpa armigera

Stress is a widespread phenomenon that all organisms must endure. Common in nature is oxidative stress, which can interrupt cell homeostasis to cause cell damage and may be derived from respiration or from environmental exposure through diet. As a result of the routine exposure from respiration, many organisms can mitigate the effects of oxidative stress, but less is known about responses to oxidative stress from other sources. Helicoverpa armigera is a major agricultural pest moth that causes significant damage to crops worldwide. Here, we examined the effects of oxidative stress on H. armigera by chronically exposing individuals to paraquat-a free radical producer-and measuring changes in development (weight, developmental rate, lifespan), and gene expression. We found that oxidative stress strongly affected development in H. armigera, with stressed samples spending more time as caterpillars than control samples (>24 vs. ~15 days, respectively) and therefore living longer overall. We found 1,618 up- and 761 down-regulated genes, respectively, in stressed versus control samples. In the up-regulated gene set, was an over-representation of biological processes related to cuticle and chitin development, glycine metabolism, and oxidation-reduction. Oxidative stress clearly impacts physiology and biochemistry in H. armigera and the interesting finding of an extended lifespan in stressed individuals could demonstrate hormesis, the phenomenon whereby toxic compounds can actually be beneficial at low doses. Collectively, our findings provide new insights into physiological and gene expression responses to oxidative stress in invertebrates.

Oxidative stress damages rRNA inside the ribosome and differentially affects the catalytic center

Intracellular levels of reactive oxygen species (ROS) increase as a consequence of oxidative stress and represent a major source of damage to biomolecules. Due to its high cellular abundance RNA is more frequently the target for oxidative damage than DNA. Nevertheless the functional consequences of damage on stable RNA are poorly understood. Using a genome-wide approach, based on 8-oxo-guanosine immunoprecipitation, we present evidence that the most abundant non-coding RNA in a cell, the ribosomal RNA (rRNA), is target for oxidative nucleobase damage by ROS. Subjecting ribosomes to oxidative stress, we demonstrate that oxidized 23S rRNA inhibits the ribosome during protein biosynthesis. Placing single oxidized nucleobases at specific position within the ribosome's catalytic center by atomic mutagenesis resulted in markedly different functional outcomes. While some active site nucleobases tolerated oxidative damage well, oxidation at others had detrimental effects on protein synthesis by inhibiting different sub-steps of the ribosomal elongation cycle. Our data provide molecular insight into the biological consequences of RNA oxidation in one of the most central cellular enzymes and reveal mechanistic insight on the role of individual active site nucleobases during translation.

Stenotrophomonas maltophilia Differential Gene Expression in Synthetic Cystic Fibrosis Sputum Reveals Shared and Cystic Fibrosis Strain-Specific Responses to the Sputum Environment

Stenotrophomonas maltophilia is a Gram-negative opportunistic pathogen that can infect the lungs of people with cystic fibrosis (CF). The highly viscous mucus in the CF lung, expectorated as sputum, serves as the primary nutrient source for microbes colonizing this site and induces virulence-associated phenotypes and gene expression in several CF pathogens. Here, we characterized the transcriptional responses of three S. maltophilia strains during exposure to synthetic CF sputum medium (SCFM2) to gain insight into how this organism interacts with the host in the CF lung. These efforts led to the identification of 881 transcripts differentially expressed by all three strains, many of which reflect the metabolic pathways used by S. maltophilia in sputum, as well as altered stress responses. The latter correlated with increased resistance to peroxide exposure after pregrowth in SCFM2 for two of the strains. We also compared the SCFM2 transcriptomes of two S. maltophilia CF isolates to that of the acute infection strain, S. maltophilia K279a, allowing us to identify CF isolate-specific signatures in differential gene expression. The expression of genes from the accessory genomes was also differentially altered in response to SCFM2. Finally, a number of biofilm-associated genes were differentially induced in SCFM2, particularly in K279a, which corresponded to increased aggregation and biofilm formation in this strain relative to both CF strains. Collectively, this work details the response of S. maltophilia to an environment that mimics important aspects of the CF lung, identifying potential survival strategies and metabolic pathways used by S. maltophilia during infections.IMPORTANCE Stenotrophomonas maltophilia is an important infecting bacterium in the airways of people with cystic fibrosis (CF). However, compared to the other CF pathogens, S. maltophilia has been relatively understudied. The significance of our research is to provide insight into the global transcriptomic changes of S. maltophilia in response to a medium that was designed to mimic important aspects of the CF lung. This study elucidates the overall metabolic changes that occur when S. maltophilia encounters the CF lung and generates a road map of candidate genes to test using in vitro and in vivo models of CF.

Dual Gene Expression Analysis Identifies Factors Associated with Staphylococcus aureus Virulence in Diabetic Mice

Staphylococcus aureus is a major human pathogen of the skin. The global burden of diabetes is high, with S. aureus being a major complication of diabetic wound infections. We investigated how the diabetic environment influences S. aureus skin infection and observed an increased susceptibility to infection in mouse models of both type I and type II diabetes. A dual gene expression approach was taken to investigate transcriptional alterations in both the host and bacterium after infection. While analysis of the host response revealed only minor changes between infected control and diabetic mice, we observed that S. aureus isolated from diabetic mice had significant increases in the levels of genes associated with translation and posttranslational modification and chaperones and reductions in the levels of genes associated with amino acid transport and metabolism. One family of genes upregulated in S. aureus isolated from diabetic lesions encoded the Clp proteases, associated with the misfolded protein response. The Clp proteases were found to be partially glucose regulated as well as influencing the hemolytic activity of S. aureus Strains lacking the Clp proteases ClpX, ClpC, and ClpP were significantly attenuated in our animal model of skin infection, with significant reductions observed in dermonecrosis and bacterial burden. In particular, mutations in clpP and clpX were significantly attenuated and remained attenuated in both normal and diabetic mice. Our data suggest that the diabetic environment also causes changes to occur in invading pathogens, and one of these virulence determinants is the Clp protease system.

Flanders hapavirus in western North America

Flanders virus (FLAV; family Rhabdoviridae) is a mosquito-borne hapavirus with no known pathology that is frequently isolated during arbovirus surveillance programs. Here, we document the presence of FLAV in Culex tarsalis mosquitoes and a Canada goose (Branta canadensis) collected in western North America, outside of the currently recognized range of FLAV. Until now, FLAV-like viruses detected in the western United States were assumed to be Hart Park virus (HPV, family Rhabdoviridae), a closely related congener. A re-examination of archived viral isolates revealed that FLAV was circulating in California as early as 1963. FLAV also was isolated in Nebraska, Colorado, South Dakota, North Dakota, and Saskatchewan, Canada. Phylogenetic analysis of the U1 pseudogene for 117 taxa and eight nuclear genes for 15 taxa demonstrated no distinct clustering between western FLAV isolates. Assuming the range of FLAV has been expanding west, these results indicate that FLAV likely spread west following multiple invasion events. However, it remains to be determined if the detection of FLAV in western North America is due to expansion or is a result of enhanced arbovirus surveillance or diagnostic techniques. Currently, the impact of FLAV infection remains unknown.

Extraction and biomolecular analysis of dermal interstitial fluid collected with hollow microneedles

Dermal interstitial fluid (ISF) is an underutilized information-rich biofluid potentially useful in health status monitoring applications whose contents remain challenging to characterize. Here, we present a facile microneedle approach for dermal ISF extraction with minimal pain and no blistering for human subjects and rats. Extracted ISF volumes were sufficient for determining transcriptome, and proteome signatures. We noted similar profiles in ISF, serum, and plasma samples, suggesting that ISF can be a proxy for direct blood sampling. Dynamic changes in RNA-seq were recorded in ISF from induced hypoxia conditions. Finally, we report the first isolation and characterization, to our knowledge, of exosomes from dermal ISF. The ISF exosome concentration is 12-13 times more enriched when compared to plasma and serum and represents a previously unexplored biofluid for exosome isolation. This minimally invasive extraction approach can enable mechanistic studies of ISF and demonstrates the potential of ISF for real-time health monitoring applications.

Template-primer binding affinity and RNase H cleavage specificity contribute to the strand transfer efficiency of HIV-1 reverse transcriptase

During reverse transcription of the HIV-1 genome, two strand-transfer events occur. Both events rely on the RNase H cleavage activity of reverse transcriptases (RTs) and template homology. Using a panel of mutants of HIV-1BH10 (group M/subtype B) and HIV-1ESP49 (group O) RTs and in vitro assays, we demonstrate that there is a strong correlation between RT minus-strand transfer efficiency and template-primer binding affinity. The highest strand transfer efficiencies were obtained with HIV-1ESP49 RT mutants containing the substitutions K358R/A359G/S360A, alone or in combination with V148I or T355A/Q357M. These HIV-1ESP49 RT mutants had been previously engineered to increase their DNA polymerase activity at high temperatures. Now, we found that RTs containing RNase H-inactivating mutations (D443N or E478Q) were devoid of strand transfer activity, whereas enzymes containing F61A or L92P had very low strand transfer activity. The strand transfer defect produced by L92P was attributed to a loss of template-primer binding affinity and, more specifically, to the higher dissociation rate constants (koff) shown by RTs bearing this substitution. Although L92P also deleteriously affected the RT's nontemplated nucleotide addition activity, neither nontemplated nucleotide addition activity nor the RT's clamp activities contributed to increased template switching when all tested mutant and WT RTs were considered. Interestingly, our results also revealed an association between efficient strand transfer and the generation of secondary cleavages in the donor RNA, consistent with the creation of invasion sites. Exposure of the elongated DNA at these sites facilitate acceptor (RNA or DNA) binding and promote template switching.

RNA Sequencing and Analysis in Microorganisms for Metabolic Network Reconstruction

There is a strict interplay between metabolic networks and transcriptional regulation in bacteria; indeed, the transcriptome regulation, affecting the expression of large gene sets, can be used to predict the likely "on" or "off" state of metabolic genes as a function of environmental factors. Up to date, many bacterial transcriptomes have been studied by RNAseq, hundreds of experiments have been performed, and Giga bases of sequences have been produced. All this transcriptional information could potentially be integrated into metabolic networks in order to obtain a more comprehensive view of their regulation and to increase their prediction power.To get high-quality transcriptomic data, to be integrated into metabolic networks, it is paramount to clearly know how to produce highly informative RNA sequencing libraries and how to manage RNA sequencing data.In this chapter, we will get across the main steps of an RNAseq experiment: from removal of ribosomal RNAs, to strand-specific library preparation, till data analysis and integration. We will try to share our experience and know-how, to give you a precise protocol to follow, and some useful recommendations or tips and tricks to adopt in order to go straightforward toward a successful RNAseq experiment.

Cancer transcriptome profiling at the juncture of clinical translation

Methodological breakthroughs over the past four decades have repeatedly revolutionized transcriptome profiling. Using RNA sequencing (RNA-seq), it has now become possible to sequence and quantify the transcriptional outputs of individual cells or thousands of samples. These transcriptomes provide a link between cellular phenotypes and their molecular underpinnings, such as mutations. In the context of cancer, this link represents an opportunity to dissect the complexity and heterogeneity of tumours and to discover new biomarkers or therapeutic strategies. Here, we review the rationale, methodology and translational impact of transcriptome profiling in cancer.

De Novo RNA Sequencing and Transcriptome Analysis of Monascus purpureus and Analysis of Key Genes Involved in Monacolin K Biosynthesis

Monascus purpureus is an important medicinal and edible microbial resource. To facilitate biological, biochemical, and molecular research on medicinal components of M. purpureus, we investigated the M. purpureus transcriptome by RNA sequencing (RNA-seq). An RNA-seq library was created using RNA extracted from a mixed sample of M. purpureus expressing different levels of monacolin K output. In total 29,713 unigenes were assembled from more than 60 million high-quality short reads. A BLAST search revealed hits for 21,331 unigenes in at least one of the protein or nucleotide databases used in this study. The 22,365 unigenes were categorized into 48 functional groups based on Gene Ontology classification. Owing to the economic and medicinal importance of M. purpureus, most studies on this organism have focused on the pharmacological activity of chemical components and the molecular function of genes involved in their biogenesis. In this study, we performed quantitative real-time PCR to detect the expression of genes related to monacolin K (mokA-mokI) at different phases (2, 5, 8, and 12 days) of M. purpureus M1 and M1-36. Our study found that mokF modulates monacolin K biogenesis in M. purpureus. Nine genes were suggested to be associated with the monacolin K biosynthesis. Studies on these genes could provide useful information on secondary metabolic processes in M. purpureus. These results indicate a detailed resource through genetic engineering of monacolin K biosynthesis in M. purpureus and related species.

A Rapid Spin Column-Based Method to Enrich Pathogen Transcripts from Eukaryotic Host Cells Prior to Sequencing

When analyzing pathogen transcriptomes during the infection of host cells, the signal-to-background (pathogen-to-host) ratio of nucleic acids (NA) in infected samples is very small. Despite the advancements in next-generation sequencing, the minute amount of pathogen NA makes standard RNA-seq library preps inadequate for effective gene-level analysis of the pathogen in cases with low bacterial loads. In order to provide a more complete picture of the pathogen transcriptome during an infection, we developed a novel pathogen enrichment technique, which can enrich for transcripts from any cultivable bacteria or virus, using common, readily available laboratory equipment and reagents. To evenly enrich for pathogen transcripts, we generate biotinylated pathogen-targeted capture probes in an enzymatic process using the entire genome of the pathogen as a template. The capture probes are hybridized to a strand-specific cDNA library generated from an RNA sample. The biotinylated probes are captured on a monomeric avidin resin in a miniature spin column, and enriched pathogen-specific cDNA is eluted following a series of washes. To test this method, we performed an in vitro time-course infection using Klebsiella pneumoniae to infect murine macrophage cells. K. pneumoniae transcript enrichment efficiency was evaluated using RNA-seq. Bacterial transcripts were enriched up to ~400-fold, and allowed the recovery of transcripts from ~2000–3600 genes not observed in untreated control samples. These additional transcripts revealed interesting aspects of K. pneumoniae biology including the expression of putative virulence factors and the expression of several genes responsible for antibiotic resistance even in the absence of drugs.

Novel Viruses Isolated from Mosquitoes in Pantanal, Brazil

Genomic sequences are described from five novel viruses and divergent strains of Brejeira and Guaico Culex viruses from mosquitoes collected in Pantanal, Brazil, in 2010.

RNA-Seq methods for transcriptome analysis

Deep sequencing has been revolutionizing biology and medicine in recent years, providing single base-level precision for our understanding of nucleic acid sequences in high throughput fashion. Sequencing of RNA, or RNA-Seq, is now a common method to analyze gene expression and to uncover novel RNA species. Aspects of RNA biogenesis and metabolism can be interrogated with specialized methods for cDNA library preparation. In this study, we review current RNA-Seq methods for general analysis of gene expression and several specific applications, including isoform and gene fusion detection, digital gene expression profiling, targeted sequencing and single-cell analysis. In addition, we discuss approaches to examine aspects of RNA in the cell, technical challenges of existing RNA-Seq methods, and future directions. WIREs RNA 2017, 8:e1364. doi: 10.1002/wrna.1364 For further resources related to this article, please visit the WIREs website.

Meeting report: discussions and preliminary findings on extracellular RNA measurement methods from laboratories in the NIH Extracellular RNA Communication Consortium

Extracellular RNAs (exRNAs) have been identified in all tested biofluids and have been associated with a variety of extracellular vesicles, ribonucleoprotein complexes and lipoprotein complexes. Much of the interest in exRNAs lies in the fact that they may serve as signalling molecules between cells, their potential to serve as biomarkers for prediction and diagnosis of disease and the possibility that exRNAs or the extracellular particles that carry them might be used for therapeutic purposes. Among the most significant bottlenecks to progress in this field is the lack of robust and standardized methods for collection and processing of biofluids, separation of different types of exRNA-containing particles and isolation and analysis of exRNAs. The Sample and Assay Standards Working Group of the Extracellular RNA Communication Consortium is a group of laboratories funded by the U.S. National Institutes of Health to develop such methods. In our first joint endeavour, we held a series of conference calls and in-person meetings to survey the methods used among our members, placed them in the context of the current literature and used our findings to identify areas in which the identification of robust methodologies would promote rapid advancements in the exRNA field.

The use of exome capture RNA-seq for highly degraded RNA with application to clinical cancer sequencing

RNA-seq by poly(A) selection is currently the most common protocol for whole transcriptome sequencing as it provides a broad, detailed, and accurate view of the RNA landscape. Unfortunately, the utility of poly(A) libraries is greatly limited when the input RNA is degraded, which is the norm for research tissues and clinical samples, especially when specimens are formalin-fixed. To facilitate the use of RNA sequencing beyond cell lines and in the clinical setting, we developed an exome-capture transcriptome protocol with greatly improved performance on degraded RNA. Capture transcriptome libraries enable measuring absolute and differential gene expression, calling genetic variants, and detecting gene fusions. Through validation against gold-standard poly(A) and Ribo-Zero libraries from intact RNA, we show that capture RNA-seq provides accurate and unbiased estimates of RNA abundance, uniform transcript coverage, and broad dynamic range. Unlike poly(A) selection and Ribo-Zero depletion, capture libraries retain these qualities regardless of RNA quality and provide excellent data from clinical specimens including formalin-fixed paraffin-embedded (FFPE) blocks. Systematic improvements across key applications of RNA-seq are shown on a cohort of prostate cancer patients and a set of clinical FFPE samples. Further, we demonstrate the utility of capture RNA-seq libraries in a patient with a highly malignant solitary fibrous tumor (SFT) enrolled in our clinical sequencing program called MI-ONCOSEQ. Capture transcriptome profiling from FFPE revealed two oncogenic fusions: the pathognomonic NAB2-STAT6 inversion and a therapeutically actionable BRAF fusion, which may drive this specific cancer's aggressive phenotype.

Hypothesis: Artifacts, Including Spurious Chimeric RNAs with a Short Homologous Sequence, Caused by Consecutive Reverse Transcriptions and Endogenous Random Primers

Recent RNA-sequencing technology and associated bioinformatics have led to identification of tens of thousands of putative human chimeric RNAs, i.e. RNAs containing sequences from two different genes, most of which are derived from neighboring genes on the same chromosome. In this essay, we redefine "two neighboring genes" as those producing individual transcripts, and point out two known mechanisms for chimeric RNA formation, i.e. transcription from a fusion gene or trans-splicing of two RNAs. By our definition, most putative RNA chimeras derived from canonically-defined neighboring genes may either be technical artifacts or be cis-splicing products of 5'- or 3'-extended RNA of either partner that is redefined herein as an unannotated gene, whereas trans-splicing events are rare in human cells. Therefore, most authentic chimeric RNAs result from fusion genes, about 1,000 of which have been identified hitherto. We propose a hypothesis of "consecutive reverse transcriptions (RTs)", i.e. another RT reaction following the previous one, for how most spurious chimeric RNAs, especially those containing a short homologous sequence, may be generated during RT, especially in RNA-sequencing wherein RNAs are fragmented. We also point out that RNA samples contain numerous RNA and DNA shreds that can serve as endogenous random primers for RT and ensuing polymerase chain reactions (PCR), creating artifacts in RT-PCR.

Transcriptomic Analysis of Yersinia enterocolitica Biovar 1B Infecting Murine Macrophages Reveals New Mechanisms of Extracellular and Intracellular Survival

Yersinia enterocolitica is typically considered an extracellular pathogen; however, during the course of an infection, a significant number of bacteria are stably maintained within host cell vacuoles. Little is known about this population and the role it plays during an infection. To address this question and to elucidate the spatially and temporally dynamic gene expression patterns of Y. enterocolitica biovar 1B through the course of an in vitro infection, transcriptome sequencing and differential gene expression analysis of bacteria infecting murine macrophage cells were performed under four distinct conditions. Bacteria were first grown in a nutrient-rich medium at 26 °C to establish a baseline of gene expression that is unrelated to infection. The transcriptomes of these bacteria were then compared to bacteria grown in a conditioned cell culture medium at 37 °C to identify genes that were differentially expressed in response to the increased temperature and medium but not in response to host cells. Infections were then performed, and the transcriptomes of bacteria found on the extracellular surface and intracellular compartments were analyzed individually. The upregulated genes revealed potential roles for a variety of systems in promoting intracellular virulence, including the Ysa type III secretion system, the Yts2 type II secretion system, and the Tad pilus. It was further determined that mutants of each of these systems had decreased virulence while infecting macrophages. Overall, these results reveal the complete set of genes expressed by Y. enterocolitica in response to infection and provide the groundwork for future virulence studies.

The rotary zone thermal cycler: a low-power system enabling automated rapid PCR

Advances in molecular biology, microfluidics, and laboratory automation continue to expand the accessibility and applicability of these methods beyond the confines of conventional, centralized laboratory facilities and into point of use roles in clinical, military, forensic, and field-deployed applications. As a result, there is a growing need to adapt the unit operations of molecular biology (e.g., aliquoting, centrifuging, mixing, and thermal cycling) to compact, portable, low-power, and automation-ready formats. Here we present one such adaptation, the rotary zone thermal cycler (RZTC), a novel wheel-based device capable of cycling up to four different fixed-temperature blocks into contact with a stationary 4-microliter capillary-bound sample to realize 1-3 second transitions with steady state heater power of less than 10 W. We demonstrate the utility of the RZTC for DNA amplification as part of a highly integrated rotary zone PCR (rzPCR) system that uses low-volume valves and syringe-based fluid handling to automate sample loading and unloading, thermal cycling, and between-run cleaning functionalities in a compact, modular form factor. In addition to characterizing the performance of the RZTC and the efficacy of different online cleaning protocols, we present preliminary results for rapid single-plex PCR, multiplex short tandem repeat (STR) amplification, and second strand cDNA synthesis.

A solvent replenishment solution for managing evaporation of biochemical reactions in air-matrix digital microfluidics devices

Digital microfluidics (DMF) is a powerful technique for sample preparation and analysis for a broad range of biological and chemical applications. In many cases, it is desirable to carry out DMF on an open surface, such that the matrix surrounding the droplets is ambient air. However, the utility of the air-matrix DMF format has been severely limited by problems with droplet evaporation, especially when the droplet-based biochemical reactions require high temperatures for long periods of time. We present a simple solution for managing evaporation in air-matrix DMF: just-in-time replenishment of the reaction volume using droplets of solvent. We demonstrate that this solution enables DMF-mediated execution of several different biochemical reactions (RNA fragmentation, first-strand cDNA synthesis, and PCR) over a range of temperatures (4-95 °C) and incubation times (up to 1 h or more) without use of oil, humidifying chambers, or off-chip heating modules. Reaction volumes and temperatures were maintained roughly constant over the course of each experiment, such that the reaction kinetics and products generated by the air-matrix DMF device were comparable to those of conventional benchscale reactions. This simple yet effective solution for evaporation management is an important advance in developing air-matrix DMF for a wide variety of new, high-impact applications, particularly in the biomedical sciences.

Nhumirim virus, a novel flavivirus isolated from mosquitoes from the Pantanal, Brazil

We describe the isolation of a novel flavivirus, isolated from a pool of mosquitoes identified as Culex (Culex) chidesteri collected in 2010 in the Pantanal region of west-central Brazil. The virus is herein designated Nhumirim virus (NHUV) after the name of the ranch from which the mosquito pool was collected. Flavivirus RNA was detected by real-time RT-PCR of homogenized mosquitoes and from the corresponding C6/36 culture supernatant. Based on full-genome sequencing, the virus isolate was genetically distinct from but most closely related to Barkedji virus (BJV), a newly described flavivirus from Senegal. Phylogenetic analysis demonstrated that NHUV grouped with mosquito-borne flaviviruses forming a clade with BJV. This clade may be genetically intermediate between the Culex-borne flaviviruses amplified by birds and the insect-only flaviviruses.

An optimized kit-free method for making strand-specific deep sequencing libraries from RNA fragments

Deep sequencing of strand-specific cDNA libraries is now a ubiquitous tool for identifying and quantifying RNAs in diverse sample types. The accuracy of conclusions drawn from these analyses depends on precise and quantitative conversion of the RNA sample into a DNA library suitable for sequencing. Here, we describe an optimized method of preparing strand-specific RNA deep sequencing libraries from small RNAs and variably sized RNA fragments obtained from ribonucleoprotein particle footprinting experiments or fragmentation of long RNAs. Our approach works across a wide range of input amounts (400 pg to 200 ng), is easy to follow and produces a library in 2–3 days at relatively low reagent cost, all while giving the user complete control over every step. Because all enzymatic reactions were optimized and driven to apparent completion, sequence diversity and species abundance in the input sample are well preserved.

Characterization of a novel insect-specific flavivirus from Brazil: potential for inhibition of infection of arthropod cells with medically important flaviviruses

In the past decade, there has been an upsurge in the number of newly described insect-specific flaviviruses isolated pan-globally. We recently described the isolation of a novel flavivirus (tentatively designated 'Nhumirim virus'; NHUV) that represents an example of a unique subset of apparently insect-specific viruses that phylogenetically affiliate with dual-host mosquito-borne flaviviruses despite appearing to be limited to replication in mosquito cells. We characterized the in vitro growth potential and 3' untranslated region (UTR) sequence homology with alternative flaviviruses, and evaluated the virus's capacity to suppress replication of representative Culex spp.-vectored pathogenic flaviviruses in mosquito cells. Only mosquito cell lines were found to support NHUV replication, further reinforcing the insect-specific phenotype of this virus. Analysis of the sequence and predicted RNA secondary structures of the 3' UTR indicated NHUV to be most similar to viruses within the yellow fever serogroup and Japanese encephalitis serogroup, and viruses in the tick-borne flavivirus clade. NHUV was found to share the fewest conserved sequence elements when compared with traditional insect-specific flaviviruses. This suggests that, despite apparently being insect specific, this virus probably diverged from an ancestral mosquito-borne flavivirus. Co-infection experiments indicated that prior or concurrent infection of mosquito cells with NHUV resulted in a significant reduction in virus production of West Nile virus (WNV), St Louis encephalitis virus (SLEV) and Japanese encephalitis virus. The inhibitory effect was most effective against WNV and SLEV with over a 10(6)-fold and 10(4)-fold reduction in peak titres, respectively.

Capture and Amplification by Tailing and Switching (CATS). An ultrasensitive ligation-independent method for generation of DNA libraries for deep sequencing from picogram amounts of DNA and RNA

Massive parallel sequencing (MPS) technologies have paved the way into new areas of research including individualized medicine. However, sequencing of trace amounts of DNA or RNA still remains a major challenge, especially for degraded nucleic acids like circulating DNA. This together with high cost and time requirements impedes many important applications of MPS in medicine and fundamental science. We have established a fast, cheap and highly efficient protocol called ‘Capture and Amplification by Tailing and Switching’ (CATS) to directly generate ready-to-sequence libraries for MPS from nanogram and picogram quantities of both DNA and RNA. Furthermore, those DNA libraries are strand-specific, can be prepared within 2–3 h and do not require preliminary sample amplification steps. To exemplify the capacity of the technique, we have generated and sequenced DNA libraries from hundred-picogram amounts of circulating nucleic acids isolated from human blood plasma, one nanogram of mRNA-enriched total RNA from cultured cells and few nanograms of bisulfite-converted DNA. The approach for DNA library preparation from minimal and fragmented input described here will find broad application in diverse research areas such as translational medicine including therapy monitoring, prediction, prognosis and early detection of various human disorders and will permit high-throughput DNA sequencing from previously inaccessible material such as minute forensic and archeological samples.

Analysis of plant microbe interactions in the era of next generation sequencing technologies

Next generation sequencing (NGS) technologies have impressively accelerated research in biological science during the last years by enabling the production of large volumes of sequence data to a drastically lower price per base, compared to traditional sequencing methods. The recent and ongoing developments in the field allow addressing research questions in plant-microbe biology that were not conceivable just a few years ago. The present review provides an overview of NGS technologies and their usefulness for the analysis of microorganisms that live in association with plants. Possible limitations of the different sequencing systems, in particular sources of errors and bias, are critically discussed and methods are disclosed that help to overcome these shortcomings. A focus will be on the application of NGS methods in metagenomic studies, including the analysis of microbial communities by amplicon sequencing, which can be considered as a targeted metagenomic approach. Different applications of NGS technologies are exemplified by selected research articles that address the biology of the plant associated microbiota to demonstrate the worth of the new methods.

RIPiT-Seq: a high-throughput approach for footprinting RNA:protein complexes

Development of high-throughput approaches to map the RNA interaction sites of individual RNA binding proteins (RBPs) transcriptome-wide is rapidly transforming our understanding of post-transcriptional gene regulatory mechanisms. Here we describe a ribonucleoprotein (RNP) footprinting approach we recently developed for identifying occupancy sites of both individual RBPs and multi-subunit RNP complexes. RNA:protein immunoprecipitation in tandem (RIPiT) yields highly specific RNA footprints of cellular RNPs isolated via two sequential purifications; the resulting RNA footprints can then be identified by high-throughput sequencing (Seq). RIPiT-Seq is broadly applicable to all RBPs regardless of their RNA binding mode and thus provides a means to map the RNA binding sites of RBPs with poor inherent ultraviolet (UV) crosslinkability. Further, among current high-throughput approaches, RIPiT has the unique capacity to differentiate binding sites of RNPs with overlapping protein composition. It is therefore particularly suited for studying dynamic RNP assemblages whose composition evolves as gene expression proceeds.

Use of a capture-based pathogen transcript enrichment strategy for RNA-Seq analysis of the Francisella tularensis LVS transcriptome during infection of murine macrophages

Francisella tularensis is a zoonotic intracellular pathogen that is capable of causing potentially fatal human infections. Like all successful bacterial pathogens, F. tularensis rapidly responds to changes in its environment during infection of host cells, and upon encountering different microenvironments within those cells. This ability to appropriately respond to the challenges of infection requires rapid and global shifts in gene expression patterns. In this study, we use a novel pathogen transcript enrichment strategy and whole transcriptome sequencing (RNA-Seq) to perform a detailed characterization of the rapid and global shifts in F. tularensis LVS gene expression during infection of murine macrophages. We performed differential gene expression analysis on all bacterial genes at two key stages of infection: phagosomal escape, and cytosolic replication. By comparing the F. tularensis transcriptome at these two stages of infection to that of the bacteria grown in culture, we were able to identify sets of genes that are differentially expressed over the course of infection. This analysis revealed the temporally dynamic expression of a number of known and putative transcriptional regulators and virulence factors, providing insight into their role during infection. In addition, we identified several F. tularensis genes that are significantly up-regulated during infection but had not been previously identified as virulence factors. These unknown genes may make attractive therapeutic or vaccine targets.