ABCG2 regulatory single-nucleotide polymorphisms alter in vivo enhancer activity and expression.

OBJECTIVES: The expression and activity of the breast cancer resistance protein (ABCG2) contributes toward the pharmacokinetics of endogenous and xenobiotic substrates. The effect of genetic variation on the activity of cis-regulatory elements and nuclear response elements in the ABCG2 locus and their contribution toward ABCG2 expression have not been investigated systematically. In this study, the effect of genetic variation on the in vitro and in vivo enhancer activity of six previously identified liver enhancers in the ABCG2 locus was examined. METHODS: Reference and variant liver enhancers were tested for their ability to alter luciferase activity in vitro in HepG2 and HEK293T cell lines and in vivo using a hydrodynamic tail vein assay. Positive in vivo single-nucleotide polymorphisms (SNPs) were tested for association with gene expression and for altered protein binding in electrophoretic mobility shift assays. RESULTS: Multiple SNPs were found to alter enhancer activity in vitro. Four of these variants (rs9999111, rs12508471, ABCG2RE1*2, and rs149713212) decreased and one (rs2725263) increased enhancer activity in vivo. In addition, rs9999111 and rs12508471 were associated with ABCG2 expression in lymphoblastoid cell lines, lymphocytes, and T cells, and showed increased HepG2 nuclear protein binding. CONCLUSION: This study identifies SNPs within regulatory regions of the ABCG2 locus that alter enhancer activity in vitro and in vivo. Several of these SNPs correlate with tissue-specific ABCG2 expression and alter DNA/protein binding. These SNPs could contribute toward reported tissue-specific variability in ABCG2 expression and may influence the correlation between ABCG2 expression and disease risk or the pharmacokinetics and pharmacodynamics of breast cancer resistance protein substrates.

In Vivo Hepatic Enhancer Elements in the Human ABCG2 Locus.

ABCG2 encodes the mitoxantrone resistance protein (MXR; breast cancer resistance protein), an ATP-binding cassette (ABC) efflux membrane transporter. Computational analysis of the ∼300 kb region of DNA surrounding ABCG2 (chr4:88911376-89220011, hg19) identified 30 regions with potential cis-regulatory capabilities. These putative regulatory regions were tested for their enhancer and suppressor activity in a human liver cell line using luciferase reporter assays. The in vitro enhancer and suppressor assays identified four regions that decreased gene expression and five regions that increased expression >1.6-fold. Four of five human hepatic in vitro enhancers were confirmed as in vivo liver enhancers using the mouse hydrodynamic tail vein injection assay. Two of the in vivo liver enhancers (ABCG2RE1 and ABCG2RE9) responded to 17ß-estradiol or rifampin in human cell lines, and ABCG2RE9 had ChIP-seq evidence to support the binding of several transcription factors and the transcriptional coactivator p300 in human hepatocytes. This study identified genomic regions surrounding human ABCG2 that can function as regulatory elements, some with the capacity to alter gene expression upon environmental stimulus. The results from this research will drive future investigations of interindividual variation in ABCG2 expression and function that contribute to differences in drug response.

Genomic Characterization of Metformin Hepatic Response.

Metformin is used as a first-line therapy for type 2 diabetes (T2D) and prescribed for numerous other diseases. However, its mechanism of action in the liver has yet to be characterized in a systematic manner. To comprehensively identify genes and regulatory elements associated with metformin treatment, we carried out RNA-seq and ChIP-seq (H3K27ac, H3K27me3) on primary human hepatocytes from the same donor treated with vehicle control, metformin or metformin and compound C, an AMP-activated protein kinase (AMPK) inhibitor (allowing to identify AMPK-independent pathways). We identified thousands of metformin responsive AMPK-dependent and AMPK-independent differentially expressed genes and regulatory elements. We functionally validated several elements for metformin-induced promoter and enhancer activity. These include an enhancer in an ataxia telangiectasia mutated (ATM) intron that has SNPs in linkage disequilibrium with a metformin treatment response GWAS lead SNP (rs11212617) that showed increased enhancer activity for the associated haplotype. Expression quantitative trait locus (eQTL) liver analysis and CRISPR activation suggest that this enhancer could be regulating ATM, which has a known role in AMPK activation, and potentially also EXPH5 and DDX10, its neighboring genes. Using ChIP-seq and siRNA knockdown, we further show that activating transcription factor 3 (ATF3), our top metformin upregulated AMPK-dependent gene, could have an important role in gluconeogenesis repression. Our findings provide a genome-wide representation of metformin hepatic response, highlight important sequences that could be associated with interindividual variability in glycemic response to metformin and identify novel T2D treatment candidates.

Systematic dissection of coding exons at single nucleotide resolution supports an additional role in cell-specific transcriptional regulation.

In addition to their protein coding function, exons can also serve as transcriptional enhancers. Mutations in these exonic-enhancers (eExons) could alter both protein function and transcription. However, the functional consequence of eExon mutations is not well known. Here, using massively parallel reporter assays, we dissect the enhancer activity of three liver eExons (SORL1 exon 17, TRAF3IP2 exon 2, PPARG exon 6) at single nucleotide resolution in the mouse liver. We find that both synonymous and non-synonymous mutations have similar effects on enhancer activity and many of the deleterious mutation clusters overlap known liver-associated transcription factor binding sites. Carrying a similar massively parallel reporter assay in HeLa cells with these three eExons found differences in their mutation profiles compared to the liver, suggesting that enhancers could have distinct operating profiles in different tissues. Our results demonstrate that eExon mutations could lead to multiple phenotypes by disrupting both the protein sequence and enhancer activity and that enhancers can have distinct mutation profiles in different cell types.

Genome-wide discovery of drug-dependent human liver regulatory elements.

Inter-individual variation in gene regulatory elements is hypothesized to play a causative role in adverse drug reactions and reduced drug activity. However, relatively little is known about the location and function of drug-dependent elements. To uncover drug-associated elements in a genome-wide manner, we performed RNA-seq and ChIP-seq using antibodies against the pregnane X receptor (PXR) and three active regulatory marks (p300, H3K4me1, H3K27ac) on primary human hepatocytes treated with rifampin or vehicle control. Rifampin and PXR were chosen since they are part of the CYP3A4 pathway, which is known to account for the metabolism of more than 50% of all prescribed drugs. We selected 227 proximal promoters for genes with rifampin-dependent expression or nearby PXR/p300 occupancy sites and assayed their ability to induce luciferase in rifampin-treated HepG2 cells, finding only 10 (4.4%) that exhibited drug-dependent activity. As this result suggested a role for distal enhancer modules, we searched more broadly to identify 1,297 genomic regions bearing a conditional PXR occupancy as well as all three active regulatory marks. These regions are enriched near genes that function in the metabolism of xenobiotics, specifically members of the cytochrome P450 family. We performed enhancer assays in rifampin-treated HepG2 cells for 42 of these sequences as well as 7 sequences that overlap linkage-disequilibrium blocks defined by lead SNPs from pharmacogenomic GWAS studies, revealing 15/42 and 4/7 to be functional enhancers, respectively. A common African haplotype in one of these enhancers in the GSTA locus was found to exhibit potential rifampin hypersensitivity. Combined, our results further suggest that enhancers are the predominant targets of rifampin-induced PXR activation, provide a genome-wide catalog of PXR targets and serve as a model for the identification of drug-responsive regulatory elements.

Genome-wide identification of signaling center enhancers in the developing limb.

The limb is widely used as a model developmental system and changes to gene expression patterns in its signaling centers, notably the zone of polarizing activity (ZPA) and the apical ectodermal ridge (AER), are known to cause limb malformations and evolutionary differences in limb morphology. Although several genes that define these limb signaling centers have been described, the identification of regulatory elements that are active within these centers has been limited. By dissecting mouse E11.5 limbs that fluorescently mark the ZPA or AER, followed by fluorescence-activated cell sorting and low-cell H3K27ac ChIP-seq, we identified thousands of specific signaling-center enhancers. Our ChIP-seq datasets show strong correlation with ZPA- and AER-expressed genes, previously characterized functional ZPA and AER enhancers and enrichment for relevant biological terms related to limb development and malformation for the neighboring genes. Using transgenic assays, we show that several of these sequences function as ZPA and AER enhancers. Our results identify novel ZPA and AER enhancers that could be important regulators of genes involved in the establishment of these specialized regions and the patterning of tetrapod limbs.

A multidisciplinary intervention to reduce antibiotic duration in lower respiratory tract infections.

OBJECTIVES: Prolonged antibiotic courses are common in patients with lower respiratory tract infections (LRTIs) and contribute to antibiotic resistance and side effects. This study describes a multidisciplinary intervention to reduce antibiotic duration in LRTI patients. METHODS: This was a prospective before-and-after intervention study conducted from November 2011 to December 2012. Antibiotic duration was recorded for 6 months for all LRTI admissions (pneumonia, exacerbation of chronic obstructive pulmonary disease, exacerbation of asthma, and other LRTIs), followed by the introduction of an intervention intended to reduce the duration of antibiotic treatment. The intervention incorporated an antibiotic duration based on the CURB65 score, automatic stop dates and pharmacist feedback to prescribers. RESULTS: Two hundred and eighty-one patients were included in the pre-intervention group and 221 in the post-intervention group. The intervention resulted in a reduction in the duration of antibiotic treatment from 8.3 to 6.8 days (P < 0.001, 18.1% relative reduction). The rate of antibiotic-related adverse effects reduced from 31% to 19% (P = 0.03, 39.3% relative reduction). There was no increase in mortality or length of stay CONCLUSIONS: A simple intervention can significantly reduce antibiotic duration and antibiotic-related side effects.

Sequence signatures extracted from proximal promoters can be used to predict distal enhancers.

BACKGROUND: Gene expression is controlled by proximal promoters and distal regulatory elements such as enhancers. While the activity of some promoters can be invariant across tissues, enhancers tend to be highly tissue-specific. RESULTS: We compiled sets of tissue-specific promoters based on gene expression profiles of 79 human tissues and cell types. Putative transcription factor binding sites within each set of sequences were used to train a support vector machine classifier capable of distinguishing tissue-specific promoters from control sequences. We obtained reliable classifiers for 92% of the tissues, with an area under the receiver operating characteristic curve between 60% (for subthalamic nucleus promoters) and 98% (for heart promoters). We next used these classifiers to identify tissue-specific enhancers, scanning distal non-coding sequences in the loci of the 200 most highly and lowly expressed genes. Thirty percent of reliable classifiers produced consistent enhancer predictions, with significantly higher densities in the loci of the most highly expressed compared to lowly expressed genes. Liver enhancer predictions were assessed in vivo using the hydrodynamic tail vein injection assay. Fifty-eight percent of the predictions yielded significant enhancer activity in the mouse liver, whereas a control set of five sequences was completely negative. CONCLUSIONS: We conclude that promoters of tissue-specific genes often contain unambiguous tissue-specific signatures that can be learned and used for the de novo prediction of enhancers.

Massively parallel decoding of mammalian regulatory sequences supports a flexible organizational model.

Despite continual progress in the cataloging of vertebrate regulatory elements, little is known about their organization and regulatory architecture. Here we describe a massively parallel experiment to systematically test the impact of copy number, spacing, combination and order of transcription factor binding sites on gene expression. A complex library of ∼5,000 synthetic regulatory elements containing patterns from 12 liver-specific transcription factor binding sites was assayed in mice and in HepG2 cells. We find that certain transcription factors act as direct drivers of gene expression in homotypic clusters of binding sites, independent of spacing between sites, whereas others function only synergistically. Heterotypic enhancers are stronger than their homotypic analogs and favor specific transcription factor binding site combinations, mimicking putative native enhancers. Exhaustive testing of binding site permutations suggests that there is flexibility in binding site order. Our findings provide quantitative support for a flexible model of regulatory element activity and suggest a framework for the design of synthetic tissue-specific enhancers.

Depression in patients with idiopathic pulmonary fibrosis.

Depression carries enormous global morbidity and is 1.5-7 times likelier to occur in individuals with chronic illness than in the general population. Idiopathic pulmonary fibrosis (IPF) has a rising incidence with a severe impact on quality of life. An indication of the prevalence of depression in this group is therefore of paramount interest. A prospective study was performed. A total of 118 participants with IPF who attended the interstitial lung disease clinic in Ninewells Hospitals, Dundee, Scotland, from May 2010 to September 2011 were recruited. Informed consent was obtained. The male to female ratio was 60:58. The Wakefield Self-assessment of Depression Inventory was used (scores ≥15 denote a depressed state). Pulmonary function tests were measured to correlate disease severity with depression scores. Of them, 58 patients had significant depressive symptoms scoring ≥15; only nine were taking antidepressant medication. The mean depression score of female participants was 15.0 ± 0.77 (SD 5.9), compared with a mean male score of 13.1 ± 0.99 (SD 7.5). Disease severity, age, duration since diagnosis and number of co-morbidities were not significantly correlated with depression. The study population had a high prevalence of depressive symptoms. Medical therapy for pulmonary fibrosis is limited and therefore palliation of symptoms and pulmonary rehabilitation form the main strategy for management. Depression should be actively screened in patients with IPF.

A compact, in vivo screen of all 6-mers reveals drivers of tissue-specific expression and guides synthetic regulatory element design.

BACKGROUND: Large-scale annotation efforts have improved our ability to coarsely predict regulatory elements throughout vertebrate genomes. However, it is unclear how complex spatiotemporal patterns of gene expression driven by these elements emerge from the activity of short, transcription factor binding sequences. RESULTS: We describe a comprehensive promoter extension assay in which the regulatory potential of all 6 base-pair (bp) sequences was tested in the context of a minimal promoter. To enable this large-scale screen, we developed algorithms that use a reverse-complement aware decomposition of the de Bruijn graph to design a library of DNA oligomers incorporating every 6-bp sequence exactly once. Our library multiplexes all 4,096 unique 6-mers into 184 double-stranded 15-bp oligomers, which is sufficiently compact for in vivo testing. We injected each multiplexed construct into zebrafish embryos and scored GFP expression in 15 tissues at two developmental time points. Twenty-seven constructs produced consistent expression patterns, with the majority doing so in only one tissue. Functional sequences are enriched near biologically relevant genes, match motifs for developmental transcription factors, and are required for enhancer activity. By concatenating tissue-specific functional sequences, we generated completely synthetic enhancers for the notochord, epidermis, spinal cord, forebrain and otic lateral line, and show that short regulatory sequences do not always function modularly. CONCLUSIONS: This work introduces a unique in vivo catalog of short, functional regulatory sequences and demonstrates several important principles of regulatory element organization. Furthermore, we provide resources for designing compact, reverse-complement aware k-mer libraries.

Peripheral nervous system genes expressed in central neurons induce growth on inhibitory substrates.

Trauma to the spinal cord and brain can result in irreparable loss of function. This failure of recovery is in part due to inhibition of axon regeneration by myelin and chondroitin sulfate proteoglycans (CSPGs). Peripheral nervous system (PNS) neurons exhibit increased regenerative ability compared to central nervous system neurons, even in the presence of inhibitory environments. Previously, we identified over a thousand genes differentially expressed in PNS neurons relative to CNS neurons. These genes represent intrinsic differences that may account for the PNS's enhanced regenerative ability. Cerebellar neurons were transfected with cDNAs for each of these PNS genes to assess their ability to enhance neurite growth on inhibitory (CSPG) or permissive (laminin) substrates. Using high content analysis, we evaluated the phenotypic profile of each neuron to extract meaningful data for over 1100 genes. Several known growth associated proteins potentiated neurite growth on laminin. Most interestingly, novel genes were identified that promoted neurite growth on CSPGs (GPX3, EIF2B5, RBMX). Bioinformatic approaches also uncovered a number of novel gene families that altered neurite growth of CNS neurons.

Pharmacogene regulatory elements: from discovery to applications.

Regulatory elements play an important role in the variability of individual responses to drug treatment. This has been established through studies on three classes of elements that regulate RNA and protein abundance: promoters, enhancers and microRNAs. Each of these elements, and genetic variants within them, are being characterized at an exponential pace by next-generation sequencing (NGS) technologies. In this review, we outline examples of how each class of element affects drug response via regulation of drug targets, transporters and enzymes. We also discuss the impact of NGS technologies such as chromatin immunoprecipitation sequencing (ChIP-Seq) and RNA sequencing (RNA-Seq), and the ramifications of new techniques such as high-throughput chromosome capture (Hi-C), chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) and massively parallel reporter assays (MPRA). NGS approaches are generating data faster than they can be analyzed, and new methods will be required to prioritize laboratory results before they are ready for the clinic. However, there is no doubt that these approaches will bring about a systems-level understanding of the interplay between genetic variants and drug response. An understanding of the importance of regulatory variants in pharmacogenomics will facilitate the identification of responders versus non-responders, the prevention of adverse effects and the optimization of therapies for individual patients.

A chemical genetic approach identifies piperazine antipsychotics as promoters of CNS neurite growth on inhibitory substrates.

Injury to the central nervous system (CNS) can result in lifelong loss of function due in part to the regenerative failure of CNS neurons. Inhibitory proteins derived from myelin and the astroglial scar are major barriers for the successful regeneration of injured CNS neurons. Previously, we described the identification of a novel compound, F05, which promotes neurite growth from neurons challenged with inhibitory substrates in vitro, and promotes axonal regeneration in vivo (Usher et al., 2010). To identify additional regeneration-promoting compounds, we used F05-induced gene expression profiles to query the Broad Institute Connectivity Map, a gene expression database of cells treated with >1300 compounds. Despite no shared chemical similarity, F05-induced changes in gene expression were remarkably similar to those seen with a group of piperazine phenothiazine antipsychotics (PhAPs). In contrast to antipsychotics of other structural classes, PhAPs promoted neurite growth of CNS neurons challenged with two different glial derived inhibitory substrates. Our pharmacological studies suggest a mechanism whereby PhAPs promote growth through antagonism of calmodulin signaling, independent of dopamine receptor antagonism. These findings shed light on mechanisms underlying neurite-inhibitory signaling, and suggest that clinically approved antipsychotic compounds may be repurposed for use in CNS injured patients.

Massively parallel functional dissection of mammalian enhancers in vivo.

The functional consequences of genetic variation in mammalian regulatory elements are poorly understood. We report the in vivo dissection of three mammalian enhancers at single-nucleotide resolution through a massively parallel reporter assay. For each enhancer, we synthesized a library of >100,000 mutant haplotypes with 2-3% divergence from the wild-type sequence. Each haplotype was linked to a unique sequence tag embedded within a transcriptional cassette. We introduced each enhancer library into mouse liver and measured the relative activities of individual haplotypes en masse by sequencing the transcribed tags. Linear regression analysis yielded highly reproducible estimates of the effect of every possible single-nucleotide change on enhancer activity. The functional consequence of most mutations was modest, with ∼22% affecting activity by >1.2-fold and ∼3% by >2-fold. Several, but not all, positions with higher effects showed evidence for purifying selection, or co-localized with known liver-associated transcription factor binding sites, demonstrating the value of empirical high-resolution functional analysis.

BioAssay Ontology (BAO): a semantic description of bioassays and high-throughput screening results.

BACKGROUND: High-throughput screening (HTS) is one of the main strategies to identify novel entry points for the development of small molecule chemical probes and drugs and is now commonly accessible to public sector research. Large amounts of data generated in HTS campaigns are submitted to public repositories such as PubChem, which is growing at an exponential rate. The diversity and quantity of available HTS assays and screening results pose enormous challenges to organizing, standardizing, integrating, and analyzing the datasets and thus to maximize the scientific and ultimately the public health impact of the huge investments made to implement public sector HTS capabilities. Novel approaches to organize, standardize and access HTS data are required to address these challenges. RESULTS: We developed the first ontology to describe HTS experiments and screening results using expressive description logic. The BioAssay Ontology (BAO) serves as a foundation for the standardization of HTS assays and data and as a semantic knowledge model. In this paper we show important examples of formalizing HTS domain knowledge and we point out the advantages of this approach. The ontology is available online at the NCBO bioportal http://bioportal.bioontology.org/ontologies/44531. CONCLUSIONS: After a large manual curation effort, we loaded BAO-mapped data triples into a RDF database store and used a reasoner in several case studies to demonstrate the benefits of formalized domain knowledge representation in BAO. The examples illustrate semantic querying capabilities where BAO enables the retrieval of inferred search results that are relevant to a given query, but are not explicitly defined. BAO thus opens new functionality for annotating, querying, and analyzing HTS datasets and the potential for discovering new knowledge by means of inference.

Transcriptional profiling of intrinsic PNS factors in the postnatal mouse.

Neurons in the peripheral nervous system (PNS) display a higher capacity to regenerate after injury than those in the central nervous system, suggesting cell specific transcriptional modules underlying axon growth and inhibition. We report a systems biology based search for PNS specific transcription factors (TFs). Messenger RNAs enriched in dorsal root ganglion (DRG) neurons compared to cerebellar granule neurons (CGNs) were identified using subtractive hybridization and DNA microarray approaches. Network and transcription factor binding site enrichment analyses were used to further identify TFs that may be differentially active. Combining these techniques, we identified 32 TFs likely to be enriched and/or active in the PNS. Twenty-five of these TFs were then tested for an ability to promote CNS neurite outgrowth in an overexpression screen. Real-time PCR and immunohistochemical studies confirmed that one representative TF, STAT3, is intrinsic to PNS neurons, and that constitutively active STAT3 is sufficient to promote CGN neurite outgrowth.

High content screening of cortical neurons identifies novel regulators of axon growth.

Neurons in the central nervous system lose their intrinsic capacity for axon regeneration as they mature, and it is widely hypothesized that changes in gene expression are responsible. Testing this hypothesis and identifying the relevant genes has been challenging because hundreds to thousands of genes are developmentally regulated in CNS neurons, but only a small subset are likely relevant to axon growth. Here we used automated high content analysis (HCA) methods to functionally test 743 plasmids encoding developmentally regulated genes in neurite outgrowth assays using postnatal cortical neurons. We identified both growth inhibitors (Ephexin, Aldolase A, Solute Carrier 2A3, and Chimerin), and growth enhancers (Doublecortin, Doublecortin-like, Kruppel-like Factor 6, and CaM-Kinase II gamma), some of which regulate established growth mechanisms like microtubule dynamics and small GTPase signaling. Interestingly, with only one exception the growth-suppressing genes were developmentally upregulated, and the growth-enhancing genes downregulated. These data provide important support for the hypothesis that developmental changes in gene expression control neurite outgrowth, and identify potential new gene targets to promote neurite outgrowth.

A modifier locus on chromosome 5 contributes to L1 cell adhesion molecule X-linked hydrocephalus in mice.

Humans with L1 cell adhesion molecule (L1CAM) mutations exhibit X-linked hydrocephalus, as well as other severe neurological disorders. L1-6D mutant mice, which are homozygous for a deletion that removes the sixth immunoglobulin-like domain of L1cam, seldom display hydrocephalus on the 129/Sv background. However, the same L1-6D mutation produces severe hydrocephalus on the C57BL/6J background. To begin to understand how L1cam deficiencies result in hydrocephalus and to identify modifier loci that contribute to X-linked hydrocephalus by genetically interacting with L1cam, we conducted a genome-wide scan on F2 L1-6D mice, bred from L1-6D 129S2/SvPasCrlf and C57BL/6J mice. Linkage studies, utilizing chi-square tests and quantitative trait loci mapping techniques, were performed. Candidate modifier loci were further investigated in an extension study. Linkage was confirmed for a locus on chromosome 5, which we named L1cam hydrocephalus modifier 1 (L1hydro1), p = 4.04 X 10(-11).

EST Express: PHP/MySQL based automated annotation of ESTs from expression libraries.

BACKGROUND: Several biological techniques result in the acquisition of functional sets of cDNAs that must be sequenced and analyzed. The emergence of redundant databases such as UniGene and centralized annotation engines such as Entrez Gene has allowed the development of software that can analyze a great number of sequences in a matter of seconds. RESULTS: We have developed "EST Express", a suite of analytical tools that identify and annotate ESTs originating from specific mRNA populations. The software consists of a user-friendly GUI powered by PHP and MySQL that allows for online collaboration between researchers and continuity with UniGene, Entrez Gene and RefSeq. Two key features of the software include a novel, simplified Entrez Gene parser and tools to manage cDNA library sequencing projects. We have tested the software on a large data set (2,016 samples) produced by subtractive hybridization. CONCLUSION: EST Express is an open-source, cross-platform web server application that imports sequences from cDNA libraries, such as those generated through subtractive hybridization or yeast two-hybrid screens. It then provides several layers of annotation based on Entrez Gene and RefSeq to allow the user to highlight useful genes and manage cDNA library projects.