Molecular indexing enables quantitative targeted RNA sequencing and reveals poor efficiencies in standard library preparations.

We present a simple molecular indexing method for quantitative targeted RNA sequencing, in which mRNAs of interest are selectively captured from complex cDNA libraries and sequenced to determine their absolute concentrations. cDNA fragments are individually labeled so that each molecule can be tracked from the original sample through the library preparation and sequencing process. Multiple copies of cDNA fragments of identical sequence become distinct through labeling, and replicate clones created during PCR amplification steps can be identified and assigned to their distinct parent molecules. Selective capture enables efficient use of sequencing for deep sampling and for the absolute quantitation of rare or transient transcripts that would otherwise escape detection by standard sequencing methods. We have also constructed a set of synthetic barcoded RNA molecules, which can be introduced as controls into the sample preparation mix and used to monitor the efficiency of library construction. The quantitative targeted sequencing revealed extremely low efficiency in standard library preparations, which were further confirmed by using synthetic barcoded RNA molecules. This finding shows that standard library preparation methods result in the loss of rare transcripts and highlights the need for monitoring library efficiency and for developing more efficient sample preparation methods.

Digital encoding of cellular mRNAs enabling precise and absolute gene expression measurement by single-molecule counting.

We present a new approach for the sensitive detection and accurate quantitation of messenger ribonucleic acid (mRNA) gene transcripts in single cells. First, the entire population of mRNAs is encoded with molecular barcodes during reverse transcription. After amplification of the gene targets of interest, molecular barcodes are counted by sequencing or scored on a simple hybridization detector to reveal the number of molecules in the starting sample. Since absolute quantities are measured, calibration to standards is unnecessary, and many of the relative quantitation challenges such as polymerase chain reaction (PCR) bias are avoided. We apply the method to gene expression analysis of minute sample quantities and demonstrate precise measurements with sensitivity down to sub single-cell levels. The method is an easy, single-tube, end point assay utilizing standard thermal cyclers and PCR reagents. Accurate and precise measurements are obtained without any need for cycle-to-cycle intensity-based real-time monitoring or physical partitioning into multiple reactions (e.g., digital PCR). Further, since all mRNA molecules are encoded with molecular barcodes, amplification can be used to generate more material for multiple measurements and technical replicates can be carried out on limited samples. The method is particularly useful for small sample quantities, such as single-cell experiments. Digital encoding of cellular content preserves true abundance levels and overcomes distortions introduced by amplification.

Counting individual DNA molecules by the stochastic attachment of diverse labels.

We implement a unique strategy for single molecule counting termed stochastic labeling, where random attachment of a diverse set of labels converts a population of identical DNA molecules into a population of distinct DNA molecules suitable for threshold detection. The conceptual framework for stochastic labeling is developed and experimentally demonstrated by determining the absolute and relative number of selected genes after stochastically labeling approximately 360,000 different fragments of the human genome. The approach does not require the physical separation of molecules and takes advantage of highly parallel methods such as microarray and sequencing technologies to simultaneously count absolute numbers of multiple targets. Stochastic labeling should be particularly useful for determining the absolute numbers of RNA or DNA molecules in single cells.

Expression profiling. Combinatorial labeling of single cells for gene expression cytometry.

We present a technically simple approach for gene expression cytometry combining next-generation sequencing with stochastic barcoding of single cells. A combinatorial library of beads bearing cell- and molecular-barcoding capture probes is used to uniquely label transcripts and reconstruct the digital gene expression profile of thousands of individual cells in a single experiment without the need for robotics or automation. We applied the technology to dissect the human hematopoietic system and to characterize heterogeneous response to in vitro stimulation. High sensitivity is demonstrated by detection of low-abundance transcripts and rare cells. Under current implementation, the technique can analyze a few thousand cells simultaneously and can readily scale to 10,000s or 100,000s of cells.

Construction of cDNA libraries for highly efficient DNA sequencing from the 3' end of expressed genes.

The majority of expressed sequence tag (EST) sequences available today have been derived from the 5' ends of cDNA clones. Obtaining high-quality DNA sequences from the 3' ends of oligo(dT)-primed cDNA on a large scale has been difficult because of slippage of the DNA polymerase enzyme used in direct PCR and cycle sequencing. With the completion of whole genome sequencing for more and more organisms, mRNA 3'-UTR sequences can be particularly useful for clustering large numbers of ESTs for the effective discrimination of individual genes and gene families. We have identified a flaw in the widely used oligo(dT) primers for cDNA synthesis, and here we describe an improved priming approach to effectively synthesize cDNA devoid of homopolymeric nucleotide stretches from mRNA poly(A) tails to enable highly efficient and reliable DNA sequence determination from 3' mRNA ends. Using this method, we produced a rat lung cDNA library and successfully sequenced the 3' ends of 98% of all attempted clones.

Full-length cDNA synthesis for long-distance RT-PCR of large mRNA transcripts.

Reverse transcription of mRNA often leads to the synthesis of partial, non-full-length cDNAs. Methods to facilitate reverse transcription across RNA regions of secondary structure, as well as enzyme modifications to eliminate RNase H activities inherent to reverse transcriptase enzymes, have been previously reported. However, because all reverse transcriptases have high error rates of polymerization, the misincorporation of nucleotides can also cause the reverse transcriptase to stumble. Hence, even in the absence of RNA secondary structure and RNase H activity, the synthesis of full-length cDNA from long mRNA transcripts still remains a challenge. We describe here the coupling of a 3'-->5' exonuclease function during reverse transcription. The incorporation of a proofreading activity, when used in conjunction with denaturant buffers and RNase H-deficient reverse transcriptases, can successfully generate full-length cDNAs of up to 15 kb.

High resolution discovery and confirmation of copy number variants in 90 Yoruba Nigerians.

BACKGROUND: Copy number variants (CNVs) account for a large proportion of genetic variation in the genome. The initial discoveries of long (> 100 kb) CNVs in normal healthy individuals were made on BAC arrays and low resolution oligonucleotide arrays. Subsequent studies that used higher resolution microarrays and SNP genotyping arrays detected the presence of large numbers of CNVs that are < 100 kb, with median lengths of approximately 10 kb. More recently, whole genome sequencing of individuals has revealed an abundance of shorter CNVs with lengths < 1 kb. RESULTS: We used custom high density oligonucleotide arrays in whole-genome scans at approximately 200-bp resolution, and followed up with a localized CNV typing array at resolutions as close as 10 bp, to confirm regions from the initial genome scans, and to detect the occurrence of sample-level events at shorter CNV regions identified in recent whole-genome sequencing studies. We surveyed 90 Yoruba Nigerians from the HapMap Project, and uncovered approximately 2,700 potentially novel CNVs not previously reported in the literature having a median length of approximately 3 kb. We generated sample-level event calls in the 90 Yoruba at nearly 9,000 regions, including approximately 2,500 regions having a median length of just approximately 200 bp that represent the union of CNVs independently discovered through whole-genome sequencing of two individuals of Western European descent. Event frequencies were noticeably higher at shorter regions < 1 kb compared to longer CNVs (> 1 kb). CONCLUSIONS: As new shorter CNVs are discovered through whole-genome sequencing, high resolution microarrays offer a cost-effective means to detect the occurrence of events at these regions in large numbers of individuals in order to gain biological insights beyond the initial discovery.

Serine 171, a conserved residue in the gamma-aminobutyric acid type A (GABAA) receptor gamma2 subunit, mediates subunit interaction and cell surface localization.

Serine 171 in the GABA(A) receptor gamma2 subunit is highly conserved in the ligand-gated ion channel superfamily. In this paper, we report that mutating serine 171 within gamma2 to glycine or cysteine prevents the interaction of gamma2 with alpha2 and beta1 when these subunits are co-expressed in human embryo kidney 293 cells, resulting in intracellular retention of gamma2. Structure analysis based on a three-dimensional homology model of gamma2 (Ernst, M., Brauchart, D., Boresch, S., and Sieghart, W. (2003) Neuroscience 119, 933-943) reveals that serine 171 may play a critical role in the formation and stabilization of an exposed turn structure that is part of the subunit interaction site. Mutation of serine 171 in the gamma2 subunit could therefore result in alteration of the structure of the subunit interaction site, preventing correct subunit assembly.

Circular rapid amplification of cDNA ends for high-throughput extension cloning of partial genes.

The rapid amplification of cDNA ends (RACE) procedure is a widely used PCR-based method to clone the cDNA ends of mRNA transcripts. Current RACE methods often produce a high background of nonspecific PCR products, which can exclude the identification of the target cDNA of interest. We describe here an improved RACE procedure using circular cDNA templates and demonstrate the successful extension cloning of 4406 cDNAs.

Cloning and characterization of a GABAA receptor gamma2 subunit variant.

We have cloned a novel gamma-aminobutyric acid type A (GABAA) receptor gamma2 subunit variant named gamma2XL. gamma2XL contains an alternatively spliced exon, resulting in the addition of 40 amino acids to the N-terminal extracellular domain between Ser171 and Tyr172. We show that gamma2XL failed to localize to the cell surface when it was coexpressed with the alpha2 and beta1 subunits in human embryonic kidney 293 cells. Expression of gamma2XL in 293 cells suppressed GABAA receptor binding in a dose-dependent manner by preventing GABAA receptor cell-surface localization. We also generated a gamma2 mutant with Ser171 and Tyr172 converted to glycine and threonine, respectively. We demonstrate that this mutant has a significantly lower affinity for the alpha2 and beta1 subunits and failed to reach the cell surface when coexpressed with these subunits. Together, our results indicate that Ser171 and Tyr172 in the gamma2 subunit constitute a critical motif. When this motif is disrupted by insertion of the alternative exon, access of the gamma2 subunit to the cell surface is prevented.

PCR isolation and cloning of novel splice variant mRNAs from known drug target genes.

Alternative splicing of pre-mRNAs is an important mechanism for the generation of vertebrate protein diversity. Unfortunately, the contribution of alternative splicing to protein diversity is currently not well understood because many full-length mRNA sequences have yet to be identified. In this report, we describe the use of RT-PCR to identify and clone 279 novel alternatively spliced mRNAs from 114 well-known drug target genes. Our findings demonstrate the existence of many novel alternatively spliced mRNA transcripts and suggest that many more genes undergo functionally significant alternative splicing than previously thought.

E. coli selection of human genes encoding secreted and membrane proteins based on cDNA fusions to a leaderless beta-lactamase reporter.

Although several signal peptide-trapping methods have been devised and used to detect signal sequences, none have relied on using E.coli to identify eukaryotic proteins with signal peptides. Here, we describe a system for selecting human secreted and membrane proteins in E. coli followed by the direct validation of secretion in human cells. The method is based on cDNA fusions to a leaderless beta-lactamase reporter gene to isolate clones encoding signal peptides of human genes. We found that beta-lactamase fusion proteins carrying a eukaryotic signal peptide at its N-terminus were able to direct their export into the periplasm in E. coli to confer survival upon challenge with carbenicillin. When libraries constructed from 5' end-enriched cDNAs fused to beta-lactamase were screened in E.coli, approximately 0.5%-1% of the cDNAs are selected, and over half of the surviving clones were found to encode for secreted fusion proteins when tested in human cells. These clones were sequenced and shown to represent human genes encoding signal peptides of secreted and membrane proteins. We conclude that this is an efficient and effective strategy to easily enrich cDNA libraries for the identification of novel genes likely to encode secreted enzymes, growth factors, and receptors.