3'Pool-seq: an optimized cost-efficient and scalable method of whole-transcriptome gene expression profiling.

BACKGROUND: The advent of Next Generation Sequencing has allowed transcriptomes to be profiled with unprecedented accuracy, but the high costs of full-length mRNA sequencing have posed a limit on the accessibility and scalability of the technology. To address this, we developed 3'Pool-seq: a simple, cost-effective, and scalable RNA-seq method that focuses sequencing to the 3'-end of mRNA. We drew from aspects of SMART-seq, Drop-seq, and TruSeq to implement an easy workflow, and optimized parameters such as input RNA concentrations, tagmentation conditions, and read depth specifically for bulk-RNA. RESULTS: Thorough optimization resulted in a protocol that takes less than 12 h to perform, does not require custom sequencing primers or instrumentation, and cuts over 90% of the costs associated with TruSeq, while still achieving accurate gene expression quantification (Pearson's correlation coefficient with ERCC theoretical concentration r = 0.96) and differential gene detection (ROC analysis of 3'Pool-seq compared to TruSeq AUC = 0.921). The 3'Pool-seq dual indexing scheme was further adapted for a 96-well plate format, and ERCC spike-ins were used to correct for potential row or column pooling effects. Transcriptional profiling of troglitazone and pioglitazone treatments at multiple doses and time points in HepG2 cells was then used to show how 3'Pool-seq could distinguish the two molecules based on their molecular signatures. CONCLUSIONS: 3'Pool-seq can accurately detect gene expression at a level that is on par with TruSeq, at one tenth of the total cost. Furthermore, its unprecedented TruSeq/Nextera hybrid indexing scheme and streamlined workflow can be applied in several different formats, including 96-well plates, which allows users to thoroughly evaluate biological systems under several conditions and timepoints. Care must be taken regarding experimental design and plate layout such that potential pooling effects can be accounted for and corrected. Lastly, further studies using multiple sets of ERCC spike-ins may be used to simulate differential gene expression in a system with known ground-state values.

How Y-Family DNA polymerase IV is more accurate than Dpo4 at dCTP insertion opposite an N2-dG adduct of benzo[a]pyrene.

To bypass DNA damage, cells have Y-Family DNA polymerases (DNAPs). One Y-Family-class includes DNAP κ and DNAP IV, which accurately insert dCTP opposite N(2)-dG adducts, including from the carcinogen benzo[a]pyrene (BP). Another class includes DNAP η and DNAP V, which insert accurately opposite UV-damage, but inaccurately opposite BP-N(2)-dG. To investigate structural differences between Y-Family-classes, regions are swapped between DNAP IV (a κ/IV-class-member) and Dpo4 (a η/V-class-member); the kinetic consequences are evaluated via primer-extension studies with a BP-N(2)-dG-containing template. Four key structural elements are revealed. (1) Y-Family DNAPs have discreet non-covalent contacts between their little finger-domain (LF-Domain) and their catalytic core-domain (CC-Domain), which we call "non-covalent bridges" (NCBs). Arg37 and Arg38 in DNAP IV's CC-Domain near the active site form a non-covalent bridge (AS-NCB) by interacting with Glu251 and Asp252, respectively, in DNAP IV's LF-Domain. Without these interactions dATP/dGTP/dTTP misinsertions increase. DNAP IV's AS-NCB suppresses misinsertions better than Dpo4's equivalent AS-NCB. (2) DNAP IV also suppresses dATP/dGTP/dTTP misinsertions via a second non-covalent bridge, which is ∼8Å from the active site (Distal-NCB). Dpo4 has no Distal-NCB, rendering it inferior at dATP/dGTP/dTTP suppression. (3) dCTP insertion is facilitated by the larger minor groove opening near the active site in DNAP IV versus Dpo4, which is sensible given that Watson/Crick-like [dCTP:BP-N(2)-dG] pairing requires the BP-moiety to be in the minor groove. (4) Compared to Dpo4, DNAP IV has a smaller major groove opening, which suppresses dGTP misinsertion, implying BP-N(2)-dG bulk in the major groove during Hoogsteen syn-adduct-dG:dGTP pairing. In summary, DNAP IV has a large minor groove opening to enhance dCTP insertion, a plugged major groove opening to suppress dGTP misinsertion, and two non-covalent bridges (near and distal to the active site) to suppress dATP/dGTP/dTTP misinsertions; collectively these four structural features enhance DNAP IV's dNTP insertion fidelity opposite a BP-N(2)-dG adduct compared to Dpo4.

A method to accurately quantitate intensities of (32)P-DNA bands when multiple bands appear in a single lane of a gel is used to study dNTP insertion opposite a benzo[a]pyrene-dG adduct by Sulfolobus DNA polymerases Dpo4 and Dbh.

Quantitating relative (32)P-band intensity in gels is desired, e.g., to study primer-extension kinetics of DNA polymerases (DNAPs). Following imaging, multiple (32)P-bands are often present in lanes. Though individual bands appear by eye to be simple and well-resolved, scanning reveals they are actually skewed-Gaussian in shape and neighboring bands are overlapping, which complicates quantitation, because slower migrating bands often have considerable contributions from the trailing edges of faster migrating bands. A method is described to accurately quantitate adjacent (32)P-bands, which relies on having a standard: a simple skewed-Gaussian curve from an analogous pure, single-component band (e.g., primer alone). This single-component scan/curve is superimposed on its corresponding band in an experimentally determined scan/curve containing multiple bands (e.g., generated in a primer-extension reaction); intensity exceeding the single-component scan/curve is attributed to other components (e.g., insertion products). Relative areas/intensities are determined via pixel analysis, from which relative molarity of components is computed. Common software is used. Commonly used alternative methods (e.g., drawing boxes around bands) are shown to be less accurate. Our method was used to study kinetics of dNTP primer-extension opposite a benzo[a]pyrene-N(2)-dG-adduct with four DNAPs, including Sulfolobus solfataricus Dpo4 and Sulfolobus acidocaldarius Dbh. Vmax/Km is similar for correct dCTP insertion with Dpo4 and Dbh. Compared to Dpo4, Dbh misinsertion is slower for dATP (∼20-fold), dGTP (∼110-fold) and dTTP (∼6-fold), due to decreases in Vmax. These findings provide support that Dbh is in the same Y-Family DNAP class as eukaryotic DNAP κ and bacterial DNAP IV, which accurately bypass N(2)-dG adducts, as well as establish the scan-method described herein as an accurate method to quantitate relative intensity of overlapping bands in a single lane, whether generated from (32)P-signals or by other means (e.g., staining).

Immediate-early and delayed primary response genes are distinct in function and genomic architecture.

The transcriptional program induced by growth factor stimulation is classically described in two stages as follows: the rapid protein synthesis-independent induction of immediate-early genes, followed by the subsequent protein synthesis-dependent induction of secondary response genes. In this study, we obtained a comprehensive view of this transcriptional program. As expected, we identified both rapid and delayed gene inductions. Surprisingly, however, a large fraction of genes induced with delayed kinetics did not require protein synthesis and therefore represented delayed primary rather than secondary response genes. Of 133 genes induced within 4 h of growth factor stimulation, 49 (37%) were immediate-early genes, 58 (44%) were delayed primary response genes, and 26 (19%) were secondary response genes. Comparison of immediate-early and delayed primary response genes revealed functional and regulatory differences. Whereas many immediate-early genes encoded transcription factors, transcriptional regulators were not prevalent among the delayed primary response genes. The lag in induction of delayed primary response compared with immediate-early mRNAs was because of delays in both transcription initiation and subsequent stages of elongation and processing. Consistent with increased abundance of RNA polymerase II at their promoters, immediate-early genes were characterized by over-representation of transcription factor binding sites and high affinity TATA boxes. Immediate-early genes also had short primary transcripts with few exons, whereas delayed primary response genes more closely resembled other genes in the genome. These findings suggest that genomic features of immediate-early genes, in contrast to the delayed primary response genes, are selected for rapid induction, consistent with their regulatory functions.