cDNA normalization by hydroxyapatite chromatography to enrich transcriptome diversity in RNA-seq applications

Construction of an experimental study and addition of adapter sequences using HiDi DNA polymerase for improving DNA normalization methods relevant to novel gene discovery

Microorganisms in the environment can be distinguished into dominant and rare microbial species based on their genes. It is difficult to obtain genetic information derived from rare microbial species (rare genes) because of the differences in relative abundance. DNA normalization is an approach that is used to obtain genetic information derived from rare microbial species from an environmental sample. This method involves the addition of adapter sequences for the amplification, denaturation, and reassociation of the DNA fragments and single-stranded DNA (ssDNA)/double-stranded DNA (dsDNA) separation. In this method, the amount of a high-copy-number of DNA fragments and a low-copy-number of DNA fragments can be equalized. Improvements in this technique are expected to provide novel genetic information or genes in rare microbial species. However, few model experimental systems have been reported to validate the DNA normalization techniques. This study is aimed to improve the DNA normalization technique used to obtain genetic information of rare genes from rare microbial species. An experimental study was constructed with two antibiotic resistance genes, whose copy numbers differed up to a million-fold. Both genes were mixed and the mixture of DNA fragments, of high- and low-copy-number, containing these genes was normalized by separating ssDNA/dsDNA fragments using hydroxyapatite. Normalized DNA fragments were introduced into Escherichia coli and DNA normalization was evaluated by counting colonies. Moreover, we improved the method to amplify a low-copy-number of DNA fragments by the addition of adapter sequences to DNA fragments using HiDi DNA polymerase to increase the efficiency of DNA normalization. This normalization method was achieved with a 100,000-fold difference. These methods allowed for quantitative evaluation of the DNA normalization efficiency. The experimental data and methods obtained in this study are expected to improve the DNA normalization efficiency to obtain novel genetic information or genes.

Comparison of normalization methods for Hi-C data

Hi-C has been predominately used to study the genome-wide interactions of genomes. In Hi-C experiments, it is believed that biases originating from different systematic deviations lead to extraneous variability among raw samples, and affect the reliability of downstream interpretations. As an important pipeline in Hi-C analysis, normalization seeks to remove the unwanted systematic biases; thus, a comparison between Hi-C normalization methods benefits their choice and the downstream analysis. In this article, a comprehensive comparison is proposed to investigate six Hi-C normalization methods in terms of multiple considerations. In light of comparison results, it has been shown that a cross-sample approach significantly outperforms individual sample methods in most considerations. The differences between these methods are analyzed, some practical recommendations are given, and the results are summarized in a table to facilitate the choice of the six normalization methods. The source code for the implementation of these methods is available at https://github.com/lhqxinghun/bioinformatics/tree/master/Hi-C/NormCompare.

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.

Single-Cell RNA-Seq by Multiple Annealing and Tailing-Based Quantitative Single-Cell RNA-Seq (MATQ-Seq)

Single-cell technologies have emerged as advanced tools to study various biological processes that demand the single cell resolution. To detect subtle heterogeneity in the transcriptome, high accuracy and sensitivity are still desired for single-cell RNA-seq. We describe here multiple annealing and dC-tailing-based quantitative single-cell RNA-seq (MATQ-seq) with ~90% capture efficiency. In addition, MATQ-seq is a total RNA assay allowing for detection of nonpolyadenylated transcripts.

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.

Capturing the 'ome': the expanding molecular toolbox for RNA and DNA library construction

All sequencing experiments and most functional genomics screens rely on the generation of libraries to comprehensively capture pools of targeted sequences. In the past decade especially, driven by the progress in the field of massively parallel sequencing, numerous studies have comprehensively assessed the impact of particular manipulations on library complexity and quality, and characterized the activities and specificities of several key enzymes used in library construction. Fortunately, careful protocol design and reagent choice can substantially mitigate many of these biases, and enable reliable representation of sequences in libraries. This review aims to guide the reader through the vast expanse of literature on the subject to promote informed library generation, independent of the application.

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.

Precise Identification of Genome-Wide Transcription Start Sites in Bacteria by 5'-Rapid Amplification of cDNA Ends (5'-RACE)

Transcription start sites are commonly used to locate promoter elements in bacterial genomes. TSS were previously studied one gene at a time, often through 5'-rapid amplification of cDNA ends (5'-RACE). This technique has now been adapted for high-throughput sequencing and can be used to precisely identify TSS in a genome-wide fashion for practically any bacterium, which greatly contributes to our understanding of gene regulatory networks in microorganisms.

Genetic Regulation of Yersinia pestis

Y. pestis exhibits dramatically different traits of pathogenicity and transmission, albeit their close genetic relationship with its ancestor-Y. pseudotuberculosis, a self-limiting gastroenteric pathogen. Y. pestis is evolved into a deadly pathogen and transmitted to mammals and/or human beings by infected flea biting or directly contacting with the infected animals. Various kinds of environmental changes are implicated into its complex life cycle and pathogenesis. Dynamic regulation of gene expression is critical for environmental adaptation or survival, primarily reflected by genetic regulation mediated by transcriptional factors and small regulatory RNAs at the transcriptional and posttranscriptional level, respectively. The effects of genetic regulation have been shown to profoundly influence Y. pestis physiology and pathogenesis such as stress resistance, biofilm formation, intracellular survival, and replication. In this chapter, we mainly summarize the progresses on popular methods of genetic regulation and on regulatory patterns and consequences of many key transcriptional and posttranscriptional regulators, with a particular emphasis on how genetic regulation influences the biofilm and virulence of Y. pestis.

Informatics for RNA Sequencing: A Web Resource for Analysis on the Cloud

Massively parallel RNA sequencing (RNA-seq) has rapidly become the assay of choice for interrogating RNA transcript abundance and diversity. This article provides a detailed introduction to fundamental RNA-seq molecular biology and informatics concepts. We make available open-access RNA-seq tutorials that cover cloud computing, tool installation, relevant file formats, reference genomes, transcriptome annotations, quality-control strategies, expression, differential expression, and alternative splicing analysis methods. These tutorials and additional training resources are accompanied by complete analysis pipelines and test datasets made available without encumbrance at www.rnaseq.wiki.

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.

Probe-Directed Degradation (PDD) for Flexible Removal of Unwanted cDNA Sequences from RNA-Seq Libraries

Most applications for RNA-seq require the depletion of abundant transcripts to gain greater coverage of the underlying transcriptome. The sequences to be targeted for depletion depend on application and species and in many cases may not be supported by commercial depletion kits. This unit describes a method for generating RNA-seq libraries that incorporates probe-directed degradation (PDD), which can deplete any unwanted sequence set, with the low-bias split-adapter method of library generation (although many other library generation methods are in principle compatible). The overall strategy is suitable for applications requiring customized sequence depletion or where faithful representation of fragment ends and lack of sequence bias is paramount. We provide guidelines to rapidly design specific probes against the target sequence, and a detailed protocol for library generation using the split-adapter method including several strategies for streamlining the technique and reducing adapter dimer content.

Selective and flexible depletion of problematic sequences from RNA-seq libraries at the cDNA stage

BACKGROUND

A major hurdle to transcriptome profiling by deep-sequencing technologies is that abundant transcripts, such as rRNAs, can overwhelm the libraries, severely reducing transcriptome-wide coverage. Methods for depletion of such unwanted sequences typically require treatment of RNA samples prior to library preparation, are costly and not suited to unusual species and applications. Here we describe Probe-Directed Degradation (PDD), an approach that employs hybridisation to DNA oligonucleotides at the single-stranded cDNA library stage and digestion with Duplex-Specific Nuclease (DSN).

RESULTS

Targeting Saccharomyces cerevisiae rRNA sequences in Illumina HiSeq libraries generated by the split adapter method we show that PDD results in efficient removal of rRNA. The probes generate extended zones of depletion as a function of library insert size and the requirements for DSN cleavage. Using intact total RNA as starting material, probes can be spaced at the minimum anticipated library size minus 20 nucleotides to achieve continuous depletion. No off-target bias is detectable when comparing PDD-treated with untreated libraries. We further provide a bioinformatics tool to design suitable PDD probe sets.

CONCLUSION

We find that PDD is a rapid procedure that results in effective and specific depletion of unwanted sequences from deep-sequencing libraries. Because PDD acts at the cDNA stage, handling of fragile RNA samples can be minimised and it should further be feasible to remediate existing libraries. Importantly, PDD preserves the original RNA fragment boundaries as is required for nucleotide-resolution footprinting or base-cleavage studies. Finally, as PDD utilises unmodified DNA oligonucleotides it can provide a low-cost option for large-scale projects, or be flexibly customised to suit different depletion targets, sample types and organisms.

Next-generation technologies and data analytical approaches for epigenomics

Epigenetics refers to the collection of heritable features that modulate the genome-environment interaction without being encoded in the actual DNA sequence. While being mitotically and sometimes even meiotically transmitted, epigenetic traits often demonstrate extensive flexibility. This allows cells to acquire diverse gene expression patterns during differentiation, but also to adapt to a changing environment. However, epigenetic alterations are not always beneficial to the organism, as they are, for example, frequently identified in human diseases such as cancer. Accurate and cost-efficient genome-scale profiling of epigenetic features is thus of major importance to pinpoint these "epimutations," for example, to monitor the epigenetic impact of environmental exposure. Over the last decade, the field of epigenetics has been revolutionized by several innovative "epigenomics" technologies exactly addressing this need. In this review, we discuss and compare widely used next-generation methods to assess DNA methylation and hydroxymethylation, noncoding RNA expression, histone modifications, and nucleosome positioning. Although recent methods are typically based on "second-generation" sequencing, we also pay attention to still commonly used array- and PCR-based methods, and look forward to the additional advantages of single-molecule sequencing. As the current bottleneck in epigenomics research is the analysis rather than generation of data, the basic difficulties and problem-solving strategies regarding data preprocessing and statistical analysis are introduced for the different technologies. Finally, we also consider the complications associated with epigenomic studies of species with yet unsequenced genomes and possible solutions.

Transcriptional responses to temperature and low oxygen stress in Atlantic salmon studied with next-generation sequencing technology

Background

Warmer seawater as a result of climate change may impose environmental challenges for Atlantic salmon aquaculture in its southernmost geographic range. Seawater temperatures above optimal level for growth may be reached in the warmest summer weeks. Caged fish can experience temperature and low oxygen saturation stress during such episodes, raising fish welfare and productivity concerns. In this work we compare the transcriptional responses in Atlantic salmon exposed to chronic high temperature (19°C) and low oxygen saturation (4-5 mg/L) stress.

Results

We used next-generation sequencing and RT-qPCR to screen for effects, and focused on growth regulation and oxidative stress in fish exposed to sub-optimal conditions. Both prolonged temperature (45 days) and low oxygen (120 days) stress had a significant negative effect on growth. The main effect of heat stress appears to be a general reduced transcriptional rate in salmon liver, while mechanisms typically associated with responses induced by chemical drugs were stimulated. Heat stress significantly down-regulated several transcripts encoding proteins involved in the protection against oxidative stress, including CuZn SOD, Mn SOD, GPx1 and GR, as well as additional stress markers HIF1A, CYP1A, MTOR and PSMC2 (RT-qPCR data). In salmon held at low oxygen concentration for four months protein ubiquitination (protein catabolism) was the most strongly affected pathway. According to the RT-qPCR data, low oxygen stress significantly up-regulated the transcriptional levels of IGFBP1B and down-regulated the levels of GR. Pathway analysis suggests that high temperature and low oxygen saturation stress affects many similar mechanisms in Atlantic salmon. Based on the gene lists, six out of the top ten predicted upstream transcriptional regulators, 1,2-dithiol-3-thione sirolimus, CD437, 5-fluorouracil, HNF4A and NFE2L2, were similar between the two treatments.

Conclusions

In conclusion, temperature and low oxygen saturation stress affect many identical mechanisms in liver cells resulting in a metabolic depression, but these effects are not necessarily mediated through altered transcription of the same genes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-817) contains supplementary material, which is available to authorized users.

Gene expression, single nucleotide variant and fusion transcript discovery in archival material from breast tumors

Advantages of RNA-Seq over array based platforms are quantitative gene expression and discovery of expressed single nucleotide variants (eSNVs) and fusion transcripts from a single platform, but the sensitivity for each of these characteristics is unknown. We measured gene expression in a set of manually degraded RNAs, nine pairs of matched fresh-frozen, and FFPE RNA isolated from breast tumor with the hybridization based, NanoString nCounter (226 gene panel) and with whole transcriptome RNA-Seq using RiboZeroGold ScriptSeq V2 library preparation kits. We performed correlation analyses of gene expression between samples and across platforms. We then specifically assessed whole transcriptome expression of lincRNA and discovery of eSNVs and fusion transcripts in the FFPE RNA-Seq data. For gene expression in the manually degraded samples, we observed Pearson correlations of >0.94 and >0.80 with NanoString and ScriptSeq protocols, respectively. Gene expression data for matched fresh-frozen and FFPE samples yielded mean Pearson correlations of 0.874 and 0.783 for NanoString (226 genes) and ScriptSeq whole transcriptome protocols respectively, p<2x10(-16). Specifically for lincRNAs, we observed superb Pearson correlation (0.988) between matched fresh-frozen and FFPE pairs. FFPE samples across NanoString and RNA-Seq platforms gave a mean Pearson correlation of 0.838. In FFPE libraries, we detected 53.4% of high confidence SNVs and 24% of high confidence fusion transcripts. Sensitivity of fusion transcript detection was not overcome by an increase in depth of sequencing up to 3-fold (increase from ~56 to ~159 million reads). Both NanoString and ScriptSeq RNA-Seq technologies yield reliable gene expression data for degraded and FFPE material. The high degree of correlation between NanoString and RNA-Seq platforms suggests discovery based whole transcriptome studies from FFPE material will produce reliable expression data. The RiboZeroGold ScriptSeq protocol performed particularly well for lincRNA expression from FFPE libraries, but detection of eSNV and fusion transcripts was less sensitive.

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.

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

Use of second generation sequencing (SGS) technologies for transcriptional profiling (RNA-Seq) has revolutionized transcriptomics, enabling measurement of RNA abundances with unprecedented specificity and sensitivity and the discovery of novel RNA species. Preparation of RNA-Seq libraries requires conversion of the RNA starting material into cDNA flanked by platform-specific adaptor sequences. Each of the published methods and commercial kits currently available for RNA-Seq library preparation suffers from at least one major drawback, including long processing times, large starting material requirements, uneven coverage, loss of strand information and high cost. We report the development of a new RNA-Seq library preparation technique that produces representative, strand-specific RNA-Seq libraries from small amounts of starting material in a fast, simple and cost-effective manner. Additionally, we have developed a new quantitative PCR-based assay for precisely determining the number of PCR cycles to perform for optimal enrichment of the final library, a key step in all SGS library preparation workflows.

Enriching pathogen transcripts from infected samples: a capture-based approach to enhanced host-pathogen RNA sequencing

To fully understand the interactions of a pathogen with its host, it is necessary to analyze the RNA transcripts of both the host and pathogen throughout the course of an infection. Although this can be accomplished relatively easily on the host side, the analysis of pathogen transcripts is complicated by the overwhelming amount of host RNA isolated from an infected sample. Even with the read depth provided by second-generation sequencing, it is extremely difficult to get enough pathogen reads for an effective gene-level analysis. In this study, we describe a novel capture-based technique and device that considerably enriches for pathogen transcripts from infected samples. This versatile method can, in principle, enrich for any pathogen in any infected sample. To test the technique's efficacy, we performed time course tissue culture infections using Rift Valley fever virus and Francisella tularensis. At each time point, RNA sequencing (RNA-Seq) was performed and the results of the treated samples were compared with untreated controls. The capture of pathogen transcripts, in all cases, led to more than an order of magnitude enrichment of pathogen reads, greatly increasing the number of genes hit, the coverage of those genes, and the depth at which each transcript was sequenced.