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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.