Mee J. Kim

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.

Coding exons function as tissue-specific enhancers of nearby genes

Enhancers are essential gene regulatory elements whose alteration can lead to morphological differences between species, developmental abnormalities, and human disease. Current strategies to identify enhancers focus primarily on noncoding sequences and tend to exclude protein coding sequences. Here, we analyzed 25 available ChIP-seq data sets that identify enhancers in an unbiased manner (H3K4me1, H3K27ac, and EP300) for peaks that overlap exons. We find that, on average, 7% of all ChIP-seq peaks overlap coding exons (after excluding for peaks that overlap with first exons). By using mouse and zebrafish enhancer assays, we demonstrate that several of these exonic enhancer (eExons) candidates can function as enhancers of their neighboring genes and that the exonic sequence is necessary for enhancer activity. Using ChIP, 3C, and DNA FISH, we further show that one of these exonic limb enhancers, Dync1i1 exon 15, has active enhancer marks and physically interacts with Dlx5/6 promoter regions 900 kb away. In addition, its removal by chromosomal abnormalities in humans could cause split hand and foot malformation 1 (SHFM1), a disorder associated with DLX5/6. These results demonstrate that DNA sequences can have a dual function, operating as coding exons in one tissue and enhancers of nearby gene(s) in another tissue, suggesting that phenotypes resulting from coding mutations could be caused not only by protein alteration but also by disrupting the regulation of another gene.

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.

Identification and characterization of novel polymorphisms in the basal promoter of the human transporter, MATE1.

Objectives Human multidrug and toxin extrusion member 1, MATE1 (SLC47A1), plays an important role in the renal and biliary excretion of endogenous and exogenous organic cations including many therapeutic drugs. In this study, we characterized the transcriptional effects of five polymorphic variants and six common haplotypes in the basal promoter region of MATE1 that were identified in 272 DNA samples from ethnically diverse US populations. Methods We measured luciferase activities of the six common promoter haplotypes of MATE1 using in-vitro and in-vivo reporter assays. Results Haplotypes that contain the most common variant (mean allele frequency in four ethnic groups: 0.322), g.–66T>C, showed a significant decrease in reporter activities compared to the reference. Two transcription factors, activating protein-1 (AP-1) and activating protein-2 repressor (AP-2rep), were predicted to bind to the promoter in the region of g.–66T>C. Results from electrophoretic mobility shift assays showed that the g.–66T allele, exhibited greater binding to AP-1 than the g.–66C allele. AP-2rep inhibited the binding of AP-1 to the MATE1 basal promoter region, and the effect was considerably greater for the g.–66T>C. These data suggest that the reduced transcriptional activity of g.–66T>C results from a reduction in the binding potency of the transcriptional activator, AP-1, and an enhanced binding potency of the repressor, AP-2rep to the MATE1 basal promoter region. Consistent with the reporter assays, MATE1 mRNA expression levels were significantly lower in kidney samples from individuals who were homozygous or heterozygous for g.–66T>C in comparison with samples from individuals who were homozygous for the g.–66T allele. Conclusion Our study suggests that the rate of transcription of MATE1 is regulated by AP-1 and AP-2rep and that a common promoter variant, g.–66T>C may affect the expression level of MATE1 in human kidney, and ultimately result in variation in drug disposition and response.

Functional Characterization of Liver Enhancers That Regulate Drug-Associated Transporters

Little is known about how genetic variations in enhancers influence drug response. In this study, we investigated whether nucleotide variations in enhancers that regulate drug transporters can alter their expression levels. Using comparative genomics and liver‐specific transcription factor binding site (TFBS) analyses, we identified evolutionary conserved regions (ECRs) surrounding nine liver membrane transporters that interact with commonly used pharmaceuticals. The top 50 ECRs were screened for enhancer activity in vivo, of which five—located around ABCB11, SLC10A1, SLCO1B1, SLCO1A2, and SLC47A1—exhibited significant enhancer activity. Common variants identified in a large ethnically diverse cohort (n = 272) were assayed for differential enhancer activity, and three variants were found to have significant effects on reporter activity as compared with the reference allele. In addition, one variant was associated with reduced SLCO1A2 mRNA expression levels in human liver tissues, and another was associated with increased methotrexate (MTX) clearance in patients. This work provides a general model for the rapid characterization of liver enhancers and identifies associations between enhancer variants and drug response.

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.

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

BackgroundGene 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.ResultsWe 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.ConclusionsWe 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.

The Hydrodynamic Tail Vein Assay as a Tool for the Study of Liver Promoters and Enhancers

The hydrodynamic tail vein injection is a technique that is used to deliver nucleic acids into live mice. Delivery through this method results in the in vivo transfection of foreign DNA primarily in the liver. Here, we describe the use of this technique to test for regulatory activity of liver promoters and enhancers, using a dual luciferase reporter system as the readable/measureable output and how this application can be used for pharmacogenomic studies.

Functional characterization of SIM1-associated enhancers

Haploinsufficiency of the single-minded homology 1 (SIM1) gene in humans and mice leads to severe obesity, suggesting that altered expression of SIM1, by way of regulatory elements such as enhancers, could predispose individuals to obesity. Here, we identified transcriptional enhancers that could regulate SIM1, using comparative genomics coupled with zebrafish and mouse transgenic enhancer assays. Owing to the dual role of Sim1 in hypothalamic development and in adult energy homeostasis, the enhancer activity of these sequences was annotated from embryonic to adult age. Of the seventeen tested sequences, two SIM1 candidate enhancers (SCE2 and SCE8) were found to have brain-enhancer activity in zebrafish. Both SCE2 and SCE8 also exhibited embryonic brain-enhancer expression in mice, and time course analysis of SCE2 activity showed overlapping expression with Sim1 from embryonic to adult age, notably in the hypothalamus in adult mice. Using a deletion series, we identified the critical region in SCE2 that is needed for enhancer activity in the developing brain. Sequencing this region in obese and lean cohorts revealed a higher prevalence of single nucleotide polymorphisms (SNPs) that were unique to obese individuals, with one variant reducing developmental-enhancer activity in zebrafish. In summary, we have characterized two brain enhancers in the SIM1 locus and identified a set of obesity-specific SNPs within one of them, which may predispose individuals to obesity.

Genetic Association Study of Adiposity and Melanocortin-4 Receptor (MC4R) Common Variants: Replication and Functional Characterization of Non-Coding Regions

Common genetic variants 3′ of MC4R within two large linkage disequilibrium (LD) blocks spanning 288 kb have been associated with common and rare forms of obesity. This large association region has not been refined and the relevant DNA segments within the association region have not been identified. In this study, we investigated whether common variants in the MC4R gene region were associated with adiposity-related traits in a biracial population-based study. Single nucleotide polymorphisms (SNPs) in the MC4R region were genotyped with a custom array and a genome-wide array and associations between SNPs and five adiposity-related traits were determined using race-stratified linear regression. Previously reported associations between lower BMI and the minor alleles of rs2229616/Val103Ile and rs52820871/Ile251Leu were replicated in white female participants. Among white participants, rs11152221 in a proximal 3′ LD block (closer to MC4R) was significantly associated with multiple adiposity traits, but SNPs in a distal 3′ LD block (farther from MC4R) were not. In a case-control study of severe obesity, rs11152221 was significantly associated. The association results directed our follow-up studies to the proximal LD block downstream of MC4R. By considering nucleotide conservation, the significance of association, and proximity to the MC4R gene, we identified a candidate MC4R regulatory region. This candidate region was sequenced in 20 individuals from a study of severe obesity in an attempt to identify additional variants, and the candidate region was tested for enhancer activity using in vivo enhancer assays in zebrafish and mice. Novel variants were not identified by sequencing and the candidate region did not drive reporter gene expression in zebrafish or mice. The identification of a putative insulator in this region could help to explain the challenges faced in this study and others to link SNPs associated with adiposity to altered MC4R expression.