Linear amplification for deep sequencing

Linear amplification for deep sequencing (LADS) is an amplification method that produces representative libraries for Illumina next-generation sequencing within 2 d. The method relies on attaching two different sequencing adapters to blunt-end repaired and A-tailed DNA fragments, wherein one of the adapters is extended with the sequence for the T7 RNA polymerase promoter. Ligated and size-selected DNA fragments are transcribed in vitro with high RNA yields. Subsequent cDNA synthesis is initiated from a primer complementary to the first adapter, ensuring that the library will only contain full-length fragments with two distinct adapters. Contrary to the severely biased representation of AT- or GC-rich fragments in standard PCR-amplified libraries, the sequence coverage in T7-amplified libraries is indistinguishable from that of nonamplified libraries. Moreover, in contrast to amplification-free methods, LADS can generate sequencing libraries from a few nanograms of DNA, which is essential for all applications in which the starting material is limited.     

Transcriptomic Analysis of the Mouse Mammary Gland Reveals New Insights for the Role of Serotonin in Lactation

Serotonin regulates numerous processes in the mammary gland. Our objective was to discover novel genes, pathways and functions which serotonin modulates during lactation. The rate limiting enzyme in the synthesis of non-neuronal serotonin is tryptophan-hydroxylase (TPH1). Therefore, we used TPH1 deficient dams (KO; serotonin deficient, n = 4) and compared them to wild-type (WT; n = 4) and rescue (RC; KO + 100 mg/kg 5-hydroxytryptophan injected daily, n = 4) dams. Mammary tissues were collected on day 10 of lactation. Total RNA extraction, amplification, library preparation and sequencing were performed following the Illumina mRNA-Seq. Overall, 97 and 204 genes (false discovery rate, FDR ≤ 0.01) exhibited a minimum of a 2-fold expression difference between WT vs. KO and WT vs. RC dams, respectively. Most differentially expressed genes were related to calcium homeostasis, apoptosis regulation, cell cycle, cell differentiation and proliferation, and the immune response. Additionally, gene set enrichment analysis using Gene Ontology and Medical Subject Headings databases revealed the alteration of several biological processes (FDR ≤ 0.01) including fat cell differentiation and lipid metabolism, regulation of extracellular signal-related kinase and mitogen-activated kinase cascades, insulin resistance, nuclear transport, membrane potential regulation, and calcium release from the endoplasmic reticulum into the cytosol. The majority of the biological processes and pathways altered in the KO dams are central for mammary gland homeostasis. Increasing peripheral serotonin in the RC dams affects specific pathways that favor lactation. Our data confirms the importance of serotonin during lactation in the mammary gland.     

Influence of RNA-Seq library construction,sampling methods,and tissue harvesting time on gene expression estimation

RNA sequencing (RNA-Seq) is popular for measuring gene expression in non-model organisms, including wild populations. While RNA-Seq can detect gene expression variation among wild-caught individuals and yield important insights into biological function, sampling methods can also affect gene expression estimates. We examined the influence of multiple technical variables on estimated gene expression in a non-model fish, the westslope cutthroat trout (Oncorhynchus clarkii lewisi), using two RNA-Seq library types: 3′ RNA-Seq (QuantSeq) and whole mRNA-Seq (NEB). We evaluated effects of dip netting versus electrofishing, and of harvesting tissue immediately versus 5 min after euthanasia on estimated gene expression in blood, gill, and muscle. We found no significant differences in gene expression between sampling methods or tissue collection times with either library type. When library types were compared using the same blood samples, 58% of genes detected by both NEB and QuantSeq showed significantly different expression between library types, and NEB detected 31% more genes than QuantSeq. Although the two library types recovered different numbers of genes and expression levels, results with NEB and QuantSeq were consistent in that neither library type showed differences in gene expression between sampling methods and tissue harvesting times. Our study suggests that researchers can safely rely on different fish sampling strategies in the field. In addition, while QuantSeq is more cost effective, NEB detects more expressed genes. Therefore, when it is crucial to detect as many genes as possible (especially low expressed genes), when alternative splicing is of interest, or when working with an organism lacking good genomic resources, whole mRNA-Seq is more powerful.     

CDNA library construction from a small amount of RNA: adaptor-ligation approach for two-round cRNA amplification using T7 and SP6 RNA polymerases

In this study, we developed a method that allows cDNA library construction from a small amount of RNA without causing serious size bias in the resulting cDNA population. For this purpose, we adopted two-round cRNA amplification by T7 and SP6 RNA polymerases. The first-round cDNAs, flanked by the promoter sequences of T7 and SP6 RNA polymerases, were synthesized from 1 microg total RNA and then subjected to two rounds of cRNA amplification. Comparison of the sizes of the first-round and the second-round cRNAs indicated that the size-bias effect of the second-round cRNA synthesis was not serious. The resultant double-stranded cDNAs were cloned into a plasmid by in vitro lambda phage recombination with an efficiency of 1.2 x 10(11) colony-forming unit/microgram of starting total RNA. Characterization of the resultant cDNA library in terms of the insert size, clone redundancy, and integrity of 3' ends of cDNAs indicated that the amplified library was comparable to a library constructed by a conventional method, although large cDNAs tend to be slightly truncated in the amplified library. This method enables the construction of a library from a small amount of RNA, and calculations suggest that the strategy would be efficient enough to use even a single cell as starting material.     

A Modified RNA-Seq Approach for Whole Genome Sequencing of RNA Viruses from Faecal and Blood Samples

To date, very large scale sequencing of many clinically important RNA viruses has been complicated by their high population molecular variation, which creates challenges for polymerase chain reaction and sequencing primer design. Many RNA viruses are also difficult or currently not possible to culture, severely limiting the amount and purity of available starting material. Here, we describe a simple, novel, high-throughput approach to Norovirus and Hepatitis C virus whole genome sequence determination based on RNA shotgun sequencing (also known as RNA-Seq). We demonstrate the effectiveness of this method by sequencing three Norovirus samples from faeces and two Hepatitis C virus samples from blood, on an Illumina MiSeq benchtop sequencer. More than 97% of reference genomes were recovered. Compared with Sanger sequencing, our method had no nucleotide differences in 14,019 nucleotides (nt) for Noroviruses (from a total of 2 Norovirus genomes obtained with Sanger sequencing), and 8 variants in 9,542 nt for Hepatitis C virus (1 variant per 1,193 nt). The three Norovirus samples had 2, 3, and 2 distinct positions called as heterozygous, while the two Hepatitis C virus samples had 117 and 131 positions called as heterozygous. To confirm that our sample and library preparation could be scaled to true high-throughput, we prepared and sequenced an additional 77 Norovirus samples in a single batch on an Illumina HiSeq 2000 sequencer, recovering >90% of the reference genome in all but one sample. No discrepancies were observed across 118,757 nt compared between Sanger and our custom RNA-Seq method in 16 samples. By generating viral genomic sequences that are not biased by primer-specific amplification or enrichment, this method offers the prospect of large-scale, affordable studies of RNA viruses which could be adapted to routine diagnostic laboratory workflows in the near future, with the potential to directly characterize within-host viral diversity.     

Chum-RNA allows preparation of a high-quality cDNA library from a single-cell quantity of mRNA without PCR amplification

Linear RNA amplification using T7 RNA polymerase is useful in genome-wide analysis of gene expression using DNA microarrays, but exponential amplification using polymerase chain reaction (PCR) is still required for cDNA library preparation from single-cell quantities of RNA. We have designed a small RNA molecule called chum-RNA that has enabled us to prepare a single-cell cDNA library after four rounds of T7-based linear amplification, without using PCR amplification. Chum-RNA drove cDNA synthesis from only 0.49 femtograms of mRNA (730 mRNA molecules) as a substrate, a quantity that corresponds to a minor population of mRNA molecules in a single mammalian cell. Analysis of the independent cDNA clone of this library (6.6 × 10⁵ cfu) suggests that 30-fold RNA amplification occurred in each round of the amplification process. The size distribution and representation of mRNAs in the resulting one-cell cDNA library retained its similarity to that of the million-cell cDNA library. The use of chum-RNA might also facilitate reactions involving other DNA/RNA modifying enzymes whose Michaelis constant (K_m) values are around 1 mM, allowing them to be activated in the presence of only small quantities of substrate.     

Method for improved Illumina sequencing library preparation using NuGEN Ovation RNA-Seq System

In this study, we tested the NuGEN Ovation RNA-Seq System for library preparation followed by next-generation sequencing on an Illumina GAIIx. The cDNA product of the NuGEN kit may have significant amounts of ssDNA with hairpin structures that are generated during the amplification process. These structures interfere with efficient downstream end repair, A-tailing, and adapter ligation, all standard steps in post-amplification sequencing library construction. We were able to increase the efficiency of sequencing library yields 4- to 6-fold or greater by treatment of NuGEN-amplified cDNA product with the single-strand endonuclease S1. These results suggest that this treatment effectively cleaves hairpin structures generated during amplification that are resistant to the standard enzyme cocktails used for the end-repair step.     

Identification of extracellular miRNA in archived serum samples by next-generation sequencing from RNA extracted using multiple methods

miRNAs act as important regulators of gene expression by promoting mRNA degradation or by attenuating protein translation. Since miRNAs are stably expressed in bodily fluids, there is growing interest in profiling these miRNAs, as it is minimally invasive and cost-effective as a diagnostic matrix. A technical hurdle in studying miRNA dynamics is the ability to reliably extract miRNA as small sample volumes and low RNA abundance create challenges for extraction and downstream applications. The purpose of this study was to develop a pipeline for the recovery of miRNA using small volumes of archived serum samples. The RNA was extracted employing several widely utilized RNA isolation kits/methods with and without addition of a carrier. The small RNA library preparation was carried out using Illumina TruSeq small RNA kit and sequencing was carried out using Illumina platform. A fraction of five microliters of total RNA was used for library preparation as quantification is below the detection limit. We were able to profile miRNA levels in serum from all the methods tested. We found out that addition of nucleic acid based carrier molecules had higher numbers of processed reads but it did not enhance the mapping of any miRBase annotated sequences. However, some of the extraction procedures offer certain advantages: RNA extracted by TRIzol seemed to align to the miRBase best; extractions using TRIzol with carrier yielded higher miRNA-to-small RNA ratios. Nuclease free glycogen can be carrier of choice for miRNA sequencing. Our findings illustrate that miRNA extraction and quantification is influenced by the choice of methodologies. Addition of nucleic acid- based carrier molecules during extraction procedure is not a good choice when assaying miRNA using sequencing. The careful selection of an extraction method permits the archived serum samples to become valuable resources for high-throughput applications.     

Scalable transcriptome preparation for massive parallel sequencing

Background

The tremendous output of massive parallel sequencing technologies requires automated robust and scalable sample preparation methods to fully exploit the new sequence capacity.

Methodology

In this study, a method for automated library preparation of RNA prior to massively parallel sequencing is presented. The automated protocol uses precipitation onto carboxylic acid paramagnetic beads for purification and size selection of both RNA and DNA. The automated sample preparation was compared to the standard manual sample preparation.

Conclusion/Significance

The automated procedure was used to generate libraries for gene expression profiling on the Illumina HiSeq 2000 platform with the capacity of 12 samples per preparation with a significantly improved throughput compared to the standard manual preparation. The data analysis shows consistent gene expression profiles in terms of sensitivity and quantification of gene expression between the two library preparation methods.     

mRNA-Seq reveals complex patterns of gene regulation and expression in the mouse skeletal muscle transcriptome associated with calorie restriction

Sarcopenia is an age-associated loss of skeletal muscle mass and strength that increases the risk of disability. Calorie restriction (CR), the consumption of fewer calories while maintaining adequate nutrition, mitigates sarcopenia and many other age-related diseases. To identify potential mechanisms by which CR preserves skeletal muscle integrity during aging, we used mRNA-Seq for deep characterization of gene regulation and mRNA abundance in skeletal muscle of old mice compared with old mice subjected to CR. mRNA-Seq revealed complex CR-associated changes in expression of mRNA isoforms, many of which occur without a change in total message abundance and thus would not be detected by methods other than mRNA-Seq. Functional annotation of differentially expressed genes reveals CR-associated upregulation of pathways involved in energy metabolism and lipid biosynthesis, and downregulation of pathways mediating protein breakdown and oxidative stress, consistent with earlier microarray-based studies. CR-associated changes not noted in previous studies involved downregulation of genes controlling actin cytoskeletal structures and muscle development. These CR-associated changes reflect generally healthier muscle, consistent with CR's mitigation of sarcopenia. mRNA-Seq generates a rich picture of the changes in gene expression associated with CR, and may facilitate identification of genes that are primary mediators of CR's effects.