Amy Webb

Profiling Solute Carrier Transporters in the Human Blood–Brain Barrier

The neuroprotective function of the blood–brain barrier (BBB) presents a major challenge for drug delivery to the central nervous system (CNS). Critical to this function, BBB membrane transporters include the ATP‐binding cassette (ABC) transporters, which limit drug penetration across the BBB, and the less‐well‐studied solute carrier (SLC) transporters. In this work, expression profiling of 359 SLC transporters, comparative expression analysis with kidney and liver, and immunoassays in brain microvessels (BMVs) identified previously unknown transporters at the human BBB.

Whole transcriptome RNA-Seq allelic expression in human brain

BackgroundMeasuring allelic RNA expression ratios is a powerful approach for detecting cis-acting regulatory variants, RNA editing, loss of heterozygosity in cancer, copy number variation, and allele-specific epigenetic gene silencing. Whole transcriptome RNA sequencing (RNA-Seq) has emerged as a genome-wide tool for identifying allelic expression imbalance (AEI), but numerous factors bias allelic RNA ratio measurements. Here, we compare RNA-Seq allelic ratios measured in nine different human brain regions with a highly sensitive and accurate SNaPshot measure of allelic RNA ratios, identifying factors affecting reliable allelic ratio measurement. Accounting for these factors, we subsequently surveyed the variability of RNA editing across brain regions and across individuals.ResultsWe find that RNA-Seq allelic ratios from standard alignment methods correlate poorly with SNaPshot, but applying alternative alignment strategies and correcting for observed biases significantly improves correlations. Deploying these methods on a transcriptome-wide basis in nine brain regions from a single individual, we identified genes with AEI across all regions (SLC1A3, NHP2L1) and many others with region-specific AEI. In dorsolateral prefrontal cortex (DLPFC) tissues from 14 individuals, we found evidence for frequent regulatory variants affecting RNA expression in tens to hundreds of genes, depending on stringency for assigning AEI. Further, we find that the extent and variability of RNA editing is similar across brain regions and across individuals.ConclusionsThese results identify critical factors affecting allelic ratios measured by RNA-Seq and provide a foundation for using this technology to screen allelic RNA expression on a transcriptome-wide basis. Using this technology as a screening tool reveals tens to hundreds of genes harboring frequent functional variants affecting RNA expression in the human brain. With respect to RNA editing, the similarities within and between individuals leads us to conclude that this post-transcriptional process is under heavy regulatory influence to maintain an optimal degree of editing for normal biological function.

Transcriptomic variation of pharmacogenes in multiple human tissues and lymphoblastoid cell lines

Variation in the expression level and activity of genes involved in drug disposition and action (‘pharmacogenes’) can affect drug response and toxicity, especially when in tissues of pharmacological importance. Previous studies have relied primarily on microarrays to understand gene expression differences, or have focused on a single tissue or small number of samples. The goal of this study was to use RNA-sequencing (RNA-seq) to determine the expression levels and alternative splicing of 389 Pharmacogenomics Research Network pharmacogenes across four tissues (liver, kidney, heart and adipose) and lymphoblastoid cell lines, which are used widely in pharmacogenomics studies. Analysis of RNA-seq data from 139 different individuals across the 5 tissues (20–45 individuals per tissue type) revealed substantial variation in both expression levels and splicing across samples and tissue types. Comparison with GTEx data yielded a consistent picture. This in-depth exploration also revealed 183 splicing events in pharmacogenes that were previously not annotated. Overall, this study serves as a rich resource for the research community to inform biomarker and drug discovery and use.

Allele-Selective Transcriptome Recruitment to Polysomes Primed for Translation: Protein- Coding and Noncoding RNAs, and RNA Isoforms

mRNA translation into proteins is highly regulated, but the role of mRNA isoforms, noncoding RNAs (ncRNAs), and genetic variants remains poorly understood. mRNA levels on polysomes have been shown to correlate well with expressed protein levels, pointing to polysomal loading as a critical factor. To study regulation and genetic factors of protein translation we measured levels and allelic ratios of mRNAs and ncRNAs (including microRNAs) in lymphoblast cell lines (LCL) and in polysomal fractions. We first used targeted assays to measure polysomal loading of mRNA alleles, confirming reported genetic effects on translation of OPRM1 and NAT1, and detecting no effect of rs1045642 (3435C>T) in ABCB1 (MDR1) on polysomal loading while supporting previous results showing increased mRNA turnover of the 3435T allele. Use of high-throughput sequencing of complete transcript profiles (RNA-Seq) in three LCLs revealed significant differences in polysomal loading of individual RNA classes and isoforms. Correlated polysomal distribution between protein-coding and non-coding RNAs suggests interactions between them. Allele-selective polysome recruitment revealed strong genetic influence for multiple RNAs, attributable either to differential expression of RNA isoforms or to differential loading onto polysomes, the latter defining a direct genetic effect on translation. Genes identified by different allelic RNA ratios between cytosol and polysomes were enriched with published expression quantitative trait loci (eQTLs) affecting RNA functions, and associations with clinical phenotypes. Polysomal RNA-Seq combined with allelic ratio analysis provides a powerful approach to study polysomal RNA recruitment and regulatory variants affecting protein translation.

Expression of mRNA transcripts encoding membrane transporters detected with whole transcriptome sequencing of human brain and liver.

Background Membrane transporters control the influx and efflux of endogenous and xenobiotic substrates, including nutrients and drugs, across cellular membranes. Objective Whole transcriptome sequencing enables simultaneous analysis of overall and allele-specific mRNA expression, and the detection of multiple RNA isoforms. Methods Here we characterize variation in RNA transcripts emanating from gene loci encoding transporters based on RNAseq data from 10 human brains (including cocaine overdose and normal brain tissues) and 12 normal livers. Results mRNA expression was detected in 65% of transporter genes in either tissue, with many genes generating multiple mRNA transcripts. Single-nucleotide polymorphisms within transporters with previous evidence for pharmacogenomics impact were detected. We also identified noncoding RNAs in the vicinity of transporter genes with potential regulatory functions. Conclusion The results obtained with RNAseq provide detailed information on transporter mRNA expression at the molecular level, affording new avenues for the study of membrane transport, with relevance to drug efficacy and toxicity.

Replication and validation of genetic polymorphisms associated with survival after allogeneic blood or marrow transplant

Multiple candidate gene-association studies of non-HLA single-nucleotide polymorphisms (SNPs) and outcomes after blood or marrow transplant (BMT) have been conducted. We identified 70 publications reporting 45 SNPs in 36 genes significantly associated with disease-related mortality, progression-free survival, transplant-related mortality, and/or overall survival after BMT. Replication and validation of these SNP associations were performed using DISCOVeRY-BMT (Determining the Influence of Susceptibility COnveying Variants Related to one-Year mortality after BMT), a well-powered genome-wide association study consisting of 2 cohorts, totaling 2888 BMT recipients with acute myeloid leukemia, acute lymphoblastic leukemia, or myelodysplastic syndrome, and their HLA-matched unrelated donors, reported to the Center for International Blood and Marrow Transplant Research. Gene-based tests were used to assess the aggregate effect of SNPs on outcome. None of the previously reported significant SNPs replicated at P < .05 in DISCOVeRY-BMT. Validation analyses showed association with one previously reported donor SNP at P < .05 and survival; more associations would be anticipated by chance alone. No gene-based tests were significant at P < .05. Functional annotation with publicly available data shows these candidate SNPs most likely do not have biochemical function; only 13% of candidate SNPs correlate with gene expression or are predicted to impact transcription factor binding. Of these, half do not impact the candidate gene of interest; the other half correlate with expression of multiple genes. These findings emphasize the peril of pursing candidate approaches and the importance of adequately powered tests of unbiased genome-wide associations with BMT clinical outcomes given the ultimate goal of improving patient outcomes.

Allelic mRNA expression imbalance in C-type lectins reveals a frequent regulatory SNP in the human surfactant protein A (SP-A) gene

Genetic variation in C-type lectins influences infectious disease susceptibility but remains poorly understood. We used allelic mRNA expression imbalance (AEI) technology for surfactant protein (SP)-A1, SP-A2, SP-D, dendritic cell-specific ICAM-3-grabbing non-integrin (DC-SIGN), macrophage mannose receptor (MRC1) and Dectin-1, expressed in human macrophages and/or lung tissues. Frequent AEI, an indicator of regulatory polymorphisms, was observed in SP-A2, SP-D and DC-SIGN. AEI was measured for SP-A2 in 38 lung tissues using four marker single-nucleotide polymorphisms (SNPs) and was confirmed by next-generation sequencing of one lung RNA sample. Genomic DNA at the SP-A2 DNA locus was sequenced by Ion Torrent technology in 16 samples. Correlation analysis of genotypes with AEI identified a haplotype block, and, specifically, the intronic SNP rs1650232 (30% minor allele frequency); the only variant consistently associated with an approximately twofold change in mRNA allelic expression. Previously shown to alter a NAGNAG splice acceptor site with likely effects on SP-A2 expression, rs1650232 generates an alternative splice variant with three additional bases at the start of exon 3. Validated as a regulatory variant, rs1650232 is in partial linkage disequilibrium with known SP-A2 marker SNPs previously associated with risk for respiratory diseases including tuberculosis. Applying functional DNA variants in clinical association studies, rather than marker SNPs, will advance our understanding of genetic susceptibility to infectious diseases.

Expression and splicing of ABC and SLC transporters in the human blood-brain barrier measured with RNAseq

ABSTRACT The blood‐brain barrier (BBB) expresses numerous membrane transporters that supply needed nutrients to the central nervous system (CNS), consisting mostly of solute carriers (SLC transporters), or remove unwanted substrates via extrusion pumps through the action of ATP binding cassette (ABC) transporters. Previous work has identified many BBB transporters using hybridization arrays or qRT‐PCR, using targeted probes. Here we have performed next‐generation sequencing of the transcriptome (RNAseq) extracted from cerebral cortex tissues and brain microvessel endothelial cells (BMEC) obtained from two donors. The same RNA samples had previously been measured for transporter expression using qRT‐PCR (Geier et al., 2013), yielding similar expression levels for overlapping mRNAs (R = 0.66, p < 0.001). RNAseq confirms a number of transporters highly enriched in BMECs (e.g., ABCB1, ABCG2, SLCO2B1, and SLC47A1), but also detects novel BMEC transporters. Multiple splice isoforms detected by RNAseq are either robustly enriched or depleted in BMECs, indicating differential RNA processing in the BBB. The Complete RNAseq data are publically available (GSE94064).

Molecular profiling of locally-advanced rectal adenocarcinoma using microRNA expression (Review)

Treatment for locally-advanced rectal cancer (LARC) typically consists of neoadjuvant chemoradiation followed by total mesorectal excision. Recently, there has been growing interest in non-operative management for patients who are medically-inoperable or wish to avoid surgical morbidity and permanent colostomy. Approximately 50% of patients who receive pre-operative neoadjuvant chemoradiation develop some degree of pathologic response. Approximately 10-20% of patients are found to have a complete pathologic response, a finding which has frequently been shown to predict better clinical outcomes, including local-regional control, distant metastasis and survival. Many recent studies have evaluated the role of molecular biomarkers in predicting response to neoadjuvant therapy. MicroRNAs (miRNAs) are an emerging class of biomarkers that have the potential to predict which patients are most likely to benefit from pre-operative therapy and from a selective surgical approach. Here, we review the published literature on microRNAs as prognostic and predictive biomarkers in rectal cancer after pre-operative therapy. In the future, the development of prospectively validated miRNA signatures will allow clinical implementation of miRNAs as prognostic and predictive signatures in LARC.

Genome-Wide Prioritization and Transcriptomics Reveal Novel Signatures Associated with Thiazide Diuretics Blood Pressure Response

Background— Thiazide diuretics are among the most commonly prescribed antihypertensives. However, <50% of thiazide-treated patients achieve blood pressure (BP) control. Herein, we used different omics (genomics and transcriptomics) to identify novel biomarkers of thiazide diuretics BP response. Methods and Results— Genome-wide analysis included 228 white hypertensives with BP determined at baseline and after 9 weeks of hydrochlorothiazide. Single-nucleotide polymorphisms with P <5×10−5 were prioritized according to their biological function, using RegulomeDB, haploreg, and Genome-Wide Annotation of Variants. The results from the prioritization approach revealed rs10995 as the most likely functional single-nucleotide polymorphism, among single-nucleotide polymorphisms tested, that has been associated with hydrochlorothiazide BP response. The rs10995 G-allele was associated with better BP response to hydrochlorothiazide versus noncarriers (&Dgr; systolic BP/&Dgr; diastolic BP: −12.3/−8.2 versus −6.8/−3.5 mm Hg, respectively, &Dgr; systolic BP P=3×10−4, &Dgr; diastolic BP P=5×10−5). This association was replicated in independent participants treated with chlorthalidone. In addition, rs10995 G-allele was associated with increased mRNA expression of VASP (vasodilator-stimulated phosphoprotein). Moreover, baseline expression of the VASP mRNA was significantly higher in 25 good responders to hydrochlorothiazide compared with 25 poor responders (P=0.01). This finding was replicated in independent participants treated with chlorthalidone (P=0.04). Last, allelic-specific expression analysis revealed a significant but modest imbalance with rs10995 and rs10156, a single-nucleotide polymorphism in high linkage disequilibrium (r2=0.7) with rs10995, which both could contribute to the observed genetic effects by affecting VASP mRNA expression. Conclusions— This study highlights the strength of using different omics to identify novel biomarkers of drug response and suggests VASP as a potential determinant of thiazide diuretics BP response. Clinical Trial Registration— URL: http://www.clinicaltrials.gov. Unique identifiers: NCT00246519 and NCT01203852.