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.

Transcriptomic RNAseq drug screen in cerebrocortical cultures: toward novel neurogenetic disease therapies.

Rare monogenic diseases affect millions worldwide; although over 4500 rare disease genotypes are known, disease-modifying drugs are available for only 5% of them. The sheer number of these conditions combined with their rarity precludes traditional costly drug discovery programs. An economically viable alternative is to repurpose established drugs for rare diseases. Many genetic diseases result from increased or decreased protein activity and identification of clinically approved drugs which moderate this pathogenic dosage holds therapeutic potential. To identify such agents for neurogenetic diseases, we have generated genome-wide transcriptome profiles of mouse primary cerebrocortical cultures grown in the presence of 218 blood-brain barrier (BBB) penetrant clinic-tested drugs. RNAseq and differential expression analyses were used to generate transcriptomic profiles; therapeutically relevant drug-gene interactions related to rare neurogenetic diseases identified in this fashion were further analyzed by quantitative reverse transcriptase-polymerase chain reaction, western blot and immunofluorescence. We have created a transcriptome-wide searchable database for easy access to the gene expression data resulting from the cerebrocortical drug screen (Neuron Screen) and have mined this data to identify a novel link between thyroid hormone and expression of the peripheral neuropathy associated gene Pmp22. Our results demonstrate the utility of cerebrocortical cultures for transcriptomic drug screening, and the database we have created will foster further discovery of novel links between over 200 clinic-tested BBB penetrant drugs and genes related to diverse neurologic conditions.

Genetics of CD33 in Alzheimer's disease and acute myeloid leukemia.

The CD33 single-nucleotide polymorphism (SNP) rs3865444 has been associated with the risk of Alzheimer's disease (AD). Rs3865444 is in linkage disequilibrium with rs12459419 which has been associated with efficacy of an acute myeloid leukemia (AML) chemotherapeutic agent based on a CD33 antibody. We seek to evaluate the extent to which CD33 genetics in AD and AML can inform one another and advance human disease therapy. We have previously shown that these SNPs are associated with skipping of CD33 exon 2 in brain mRNA. Here, we report that these CD33 SNPs are associated with exon 2 skipping in leukocytes from AML patients and with a novel CD33 splice variant that retains CD33 intron 1. Each copy of the minor rs12459419T allele decreases prototypic full-length CD33 expression by ∼ 25% and decreases the AD odds ratio by ∼ 0.10. These results suggest that CD33 antagonists may be useful in reducing AD risk. CD33 inhibitors may include humanized CD33 antibodies such as lintuzumab which was safe but ineffective in AML clinical trials. Here, we report that lintuzumab downregulates cell-surface CD33 by 80% in phorbol-ester differentiated U937 cells, at concentrations as low as 10 ng/ml. Overall, we propose a model wherein a modest effect on RNA splicing is sufficient to mediate the CD33 association with AD risk and suggest the potential for an anti-CD33 antibody as an AD-relevant pharmacologic agent.

HIF1α is required for survival maintenance of chronic myeloid leukemia stem cells.

Hypoxia-inducible factor-1α (HIF1α), a master transcriptional regulator of the cellular and systemic hypoxia response, is essential for the maintenance of self-renewal capacity of normal HSCs. It is still unknown whether HIF1α has a role in survival regulation of leukemia stem cells (LSCs) in chronic myeloid leukemia (CML). Using a mouse model of CML, here we report that HIF1α plays a crucial role in survival maintenance of LSCs. Deletion of HIF1α impairs the propagation of CML through impairing cell-cycle progression and inducing apoptosis of LSCs. Deletion of HIF1α results in elevated expression of p16(Ink4a) and p19(Arf) in LSCs, and knockdown of p16(Ink4a) and p19(Arf) rescues the defective colony-forming ability of HIF1α(-/-) LSCs. Compared with normal HSCs, LSCs appear to be more dependent on the HIF1α pathway. Together, these results demonstrate that HIF1α represents a critical pathway in LSCs and inhibition of the HIF1α pathway provides a therapeutic strategy for eradicating LSCs in CML.

Roadmap for the use of base editors to decipher drug mechanism of action.

CRISPR base editors are powerful tools for large-scale mutagenesis studies. This kind of approach can elucidate the mechanism of action of compounds, a key process in drug discovery. Here, we explore the utility of base editors in an early drug discovery context focusing on G-protein coupled receptors. A pooled mutagenesis screening framework was set up based on a modified version of the CRISPR-X base editor system. We determine optimized experimental conditions for mutagenesis where sgRNAs are delivered by cell transfection or viral infection over extended time periods (>14 days), resulting in high mutagenesis produced in a short region located at -4/+8 nucleotides with respect to the sgRNA match. The ß2 Adrenergic Receptor (B2AR) was targeted in this way employing a 6xCRE-mCherry reporter system to monitor its response to isoproterenol. The results of our screening indicate that residue 184 of B2AR is crucial for its activation. Based on our experience, we outline the crucial points to consider when designing and performing CRISPR-based pooled mutagenesis screening, including the typical technical hurdles encountered when studying compound pharmacology.

Activated microglia cause metabolic disruptions in developmental cortical interneurons that persist in interneurons from individuals with schizophrenia.

The mechanisms by which prenatal immune activation increase the risk for neuropsychiatric disorders are unclear. Here, we generated developmental cortical interneurons (cINs)-which are known to be affected in schizophrenia (SCZ) when matured-from induced pluripotent stem cells (iPSCs) derived from healthy controls (HCs) and individuals with SCZ and co-cultured them with or without activated microglia. Co-culture with activated microglia disturbed metabolic pathways, as indicated by unbiased transcriptome analyses, and impaired mitochondrial function, arborization, synapse formation and synaptic GABA release. Deficits in mitochondrial function and arborization were reversed by alpha lipoic acid and acetyl-L-carnitine treatments, which boost mitochondrial function. Notably, activated-microglia-conditioned medium altered metabolism in cINs and iPSCs from HCs but not in iPSCs from individuals with SCZ or in glutamatergic neurons. After removal of activated-microglia-conditioned medium, SCZ cINs but not HC cINs showed prolonged metabolic deficits, which suggests that there is an interaction between SCZ genetic backgrounds and environmental risk factors.

Cell type-specific genetic regulation of gene expression across human tissues.

The Genotype-Tissue Expression (GTEx) project has identified expression and splicing quantitative trait loci in cis (QTLs) for the majority of genes across a wide range of human tissues. However, the functional characterization of these QTLs has been limited by the heterogeneous cellular composition of GTEx tissue samples. We mapped interactions between computational estimates of cell type abundance and genotype to identify cell type-interaction QTLs for seven cell types and show that cell type-interaction expression QTLs (eQTLs) provide finer resolution to tissue specificity than bulk tissue cis-eQTLs. Analyses of genetic associations with 87 complex traits show a contribution from cell type-interaction QTLs and enables the discovery of hundreds of previously unidentified colocalized loci that are masked in bulk tissue.

Using an atlas of gene regulation across 44 human tissues to inform complex disease- and trait-associated variation.

We apply integrative approaches to expression quantitative loci (eQTLs) from 44 tissues from the Genotype-Tissue Expression project and genome-wide association study data. About 60% of known trait-associated loci are in linkage disequilibrium with a cis-eQTL, over half of which were not found in previous large-scale whole blood studies. Applying polygenic analyses to metabolic, cardiovascular, anthropometric, autoimmune, and neurodegenerative traits, we find that eQTLs are significantly enriched for trait associations in relevant pathogenic tissues and explain a substantial proportion of the heritability (40-80%). For most traits, tissue-shared eQTLs underlie a greater proportion of trait associations, although tissue-specific eQTLs have a greater contribution to some traits, such as blood pressure. By integrating information from biological pathways with eQTL target genes and applying a gene-based approach, we validate previously implicated causal genes and pathways, and propose new variant and gene associations for several complex traits, which we replicate in the UK BioBank and BioVU.

Developmental and genetic regulation of the human cortex transcriptome illuminate schizophrenia pathogenesis.

Genome-wide association studies have identified 108 schizophrenia risk loci, but biological mechanisms for individual loci are largely unknown. Using developmental, genetic and illness-based RNA sequencing expression analysis in human brain, we characterized the human brain transcriptome around these loci and found enrichment for developmentally regulated genes with novel examples of shifting isoform usage across pre- and postnatal life. We found widespread expression quantitative trait loci (eQTLs), including many with transcript specificity and previously unannotated sequence that were independently replicated. We leveraged this general eQTL database to show that 48.1% of risk variants for schizophrenia associate with nearby expression. We lastly found 237 genes significantly differentially expressed between patients and controls, which replicated in an independent dataset, implicated synaptic processes, and were strongly regulated in early development. These findings together offer genetics- and diagnosis-related targets for better modeling of schizophrenia risk. This resource is publicly available at http://eqtl.brainseq.org/phase1 .

Translating human genetics into novel treatment targets for schizophrenia.

Human genetics is a rational starting point for target identification in drug discovery, yet this approach has found little application in neuroscience. Recent large-scale analyses have begun to identify robust genetic loci for schizophrenia, providing an opportunity to derive novel drug targets. Here, we summarize a strategy for applying human genetics to neuroscience drug discovery.