Matthew E. Johnson

Genetic Susceptibility to Type 2 Diabetes and Obesity: Follow-Up of Findings from Genome-Wide Association Studies

Elucidating the underlying genetic variations influencing various complex diseases is one of the major challenges currently facing clinical genetic research. Although these variations are often difficult to uncover, approaches such as genome-wide association studies (GWASs) have been successful at finding statistically significant associations between specific genomic loci and disease susceptibility. GWAS has been especially successful in elucidating genetic variants that influence type 2 diabetes (T2D) and obesity/body mass index (BMI). Specifically, several GWASs have confirmed that a variant in transcription factor 7-like 2 (TCF7L2) confers risk for T2D, while a variant in fat mass and obesity-associated protein (FTO) confers risk for obesity/BMI; indeed both of these signals are considered the most statistically associated loci discovered for these respective traits to date. The discovery of these two key loci in this context has been invaluable for providing novel insight into mechanisms of heritability and disease pathogenesis. As follow-up studies of TCF7L2 and FTO have typically lead the way in how to follow up a GWAS discovery, we outline what has been learned from such investigations and how they have implications for the myriad of other loci that have been subsequently reported in this disease context.

DNA binding by FOXP3 domain-swapped dimer suggests mechanisms of long-range chromosomal interactions

FOXP3 is a lineage-specific transcription factor that is required for regulatory T cell development and function. In this study, we determined the crystal structure of the FOXP3 forkhead domain bound to DNA. The structure reveals that FOXP3 can form a stable domain-swapped dimer to bridge DNA in the absence of cofactors, suggesting that FOXP3 may play a role in long-range gene interactions. To test this hypothesis, we used circular chromosome conformation capture coupled with high throughput sequencing (4C-seq) to analyze FOXP3-dependent genomic contacts around a known FOXP3-bound locus, Ptpn22. Our studies reveal that FOXP3 induces significant changes in the chromatin contacts between the Ptpn22 locus and other Foxp3-regulated genes, reflecting a mechanism by which FOXP3 reorganizes the genome architecture to coordinate the expression of its target genes. Our results suggest that FOXP3 mediates long-range chromatin interactions as part of its mechanisms to regulate specific gene expression in regulatory T cells.

Two novel type 2 diabetes loci revealed through integration of TCF7L2 DNA occupancy and SNP association data

Background The transcription factor 7-like 2 (TCF7L2) locus is strongly implicated in the pathogenesis of type 2 diabetes (T2D). We previously mapped the genomic regions bound by TCF7L2 using ChIP (chromatin immunoprecipitation)-seq in the colorectal carcinoma cell line, HCT116, revealing an unexpected highly significant over-representation of genome-wide association studies (GWAS) loci associated primarily with endocrine (in particular T2D) and cardiovascular traits. Methods In order to further explore if this observed phenomenon occurs in other cell lines, we carried out ChIP-seq in HepG2 cells and leveraged ENCODE data for five additional cell lines. Given that only a minority of the predicted genetic component to most complex traits has been identified to date, plus our GWAS-related observations with respect to TCF7L2 occupancy, we investigated if restricting association analyses to the genes yielded from this approach, in order to reduce the constraints of multiple testing, could reveal novel T2D loci. Results We found strong evidence for the continued enrichment of endocrine and cardiovascular GWAS categories, with additional support for cancer. When investigating all the known GWAS loci bound by TCF7L2 in the shortest gene list, derived from HCT116, the coronary artery disease-associated variant, rs46522 at the UBE2Z-GIP-ATP5G1-SNF8 locus, yielded significant association with T2D within DIAGRAM. Furthermore, when we analyzed tag-SNPs (single nucleotide polymorphisms) in genes not previously implicated by GWAS but bound by TCF7L2 within 5 kb, we observed a significant association of rs4780476 within CPPED1 in DIAGRAM. Conclusions ChIP-seq data generated with this GWAS-implicated transcription factor provided a biologically plausible method to limit multiple testing in the assessment of genome-wide genotyping data to uncover two novel T2D-associated loci.

Characterization of the transcriptional machinery bound across the widely presumed type 2 diabetes causal variant, rs7903146, within TCF7L2.

Resolving the underlying functional mechanism to a given genetic association has proven extremely challenging. However, the strongest associated type 2 diabetes (T2D) locus reported to date, TCF7L2, presents an opportunity for translational analyses, as many studies in multiple ethnicities strongly point to SNP rs7903146 in intron 3 as being the causal variant within this gene. We carried out oligo pull-down combined with mass spectrophotometry (MS) to elucidate the specific transcriptional machinery across this SNP using protein extracts from HCT116 cells. We observed that poly (ADP-ribose) polymerase 1 (PARP-1) is by far the most abundant binding factor. Pursuing the possibility of a feedback mechanism, we observed that PARP-1, along with the next most abundant binding proteins, DNA topoisomerase I and ATP-dependent RNA helicase A, dimerize with the TCF7L2 protein and with each other. We uncovered further evidence of a feedback mechanism using a luciferase reporter approach, including observing expression differences between alleles for rs7903146. We also found that there was an allelic difference in the MS results for proteins with less abundant binding, namely X-ray repair cross-complementing 5 and RPA/p70. Our results point to a protein complex binding across rs7903146 within TCF7L2 and suggests a possible mechanism by which this locus confers its T2D risk.

Genome-wide analyses of ChIP-Seq derived FOXA2 DNA occupancy in liver points to genetic networks underpinning multiple complex traits.

BACKGROUND Forkhead Box A2 (FOXA2) exerts an influence on glucose homeostasis via activity in the liver. In addition, a key genome-wide association study (GWAS) recently demonstrated that genetic variation, namely rs6048205, at the FOXA2 locus is robustly associated with fasting glucose levels. Our hypothesis was that this DNA-binding protein regulates the expression of a set of molecular pathways critical to endocrine traits. METHODS Drawing on our laboratory and bioinformatic experience with chromatin immunoprecipitation followed by massively parallel sequencing, we analyzed our existing FOXA2 chromatin immunoprecipitation followed by massively parallel sequencing data generated in human liver, using the algorithm hypergeometric optimization of motif enrichment, to gain insight into its global genomic binding pattern from a disease perspective. RESULTS We performed a pathway analysis of the gene list using the gene set enrichment analysis algorithm, which yielded a number of significant annotations. Motivated by the fact that the FOXA2 locus has been implicated by GWAS, we cross-referenced the occupancy sites with the National Institutes of Health GWAS catalog and found strong evidence for the enrichment of loci implicated in endocrine, neuropsychiatric, cardiovascular, and cancer trait categories, but interestingly there was no evidence for enrichment for inflammation related traits. Intriguingly, a FOXA2 occupancy site coincided with rs6048205, suggesting that this variant confers its effect, at least partially, via a perturbation of a FOXA2 feedback mechanism. CONCLUSION Our data strongly suggest that FOXA2 is acting as a master regulator of key pathways that are enriched for loci implicated by GWAS for most trait categories, with the clear exception of inflammation, suggesting that this factor exerts its effect in this context via noninflammatory processes.

The Loss of TET2 Promotes CD8+ T Cell Memory Differentiation

T cell differentiation requires appropriate regulation of DNA methylation. In this article, we demonstrate that the methylcytosine dioxygenase ten-eleven translocation (TET)2 regulates CD8+ T cell differentiation. In a murine model of acute viral infection, TET2 loss promotes early acquisition of a memory CD8+ T cell fate in a cell-intrinsic manner without disrupting Ag-driven cell expansion or effector function. Upon secondary recall, TET2-deficient memory CD8+ T cells demonstrate superior pathogen control. Genome-wide methylation analysis identified a number of differentially methylated regions in TET2-deficient versus wild-type CD8+ T cells. These differentially methylated regions did not occur at the loci of differentially expressed memory markers; rather, several hypermethylated regions were identified in known transcriptional regulators of CD8+ T cell memory fate. Together, these data demonstrate that TET2 is an important regulator of CD8+ T cell fate decisions.

Genome-scale Capture C promoter interaction analysis implicates novel effector genes at GWAS loci for bone mineral density

Osteoporosis is a devastating disease with an essential genetic component. Genome wide association studies (GWAS) have discovered genetic variants robustly associated with bone mineral density (BMD), however they only report genomic signals and not necessarily the precise localization of culprit effector genes. Therefore, we sought to carry out physical and direct ‘variant to gene mapping’ in a relevant primary human cell type. We developed ‘SPATIaL-seq’ (genome-Scale, Promoter-focused Analysis of chromaTIn Looping), a massively parallel, high resolution Capture-C based method to simultaneously characterize the genome-wide interactions of all human promoters. By intersecting our SPATIaL-seq and ATAC-seq data from human mesenchymal progenitor cell -derived osteoblasts, we observed consistent contacts between candidate causal variants and putative target gene promoters in open chromatin for ~30% of the 110 BMD loci investigated. Knockdown of two novel implicated genes, ING3 at ‘CPED1-WNT16’ and EPDR1 at ‘STARD3NL’, had pronounced inhibitory effects on osteoblastogenesis. Our approach therefore aids target discovery in osteoporosis and can be applied to other common genetic diseases.

Genetic And Epigenetic Fine Mapping Of Complex Trait Associated Loci In The Human Liver

Deciphering the impact of genetic variation on gene regulation is fundamental to understanding common, complex human diseases. Although histone modifications are important markers of gene regulatory regions of the genome, any specific histone modification has not been assayed in more than a few individuals in the human liver. As a result, the impacts of genetic variation that direct histone modification states in the liver are poorly understood. Here, we generate the most comprehensive genome-wide dataset of two epigenetic marks, H3K4me3 and H3K27ac, and annotate thousands of putative regulatory elements in the human liver. We integrate these findings with genome-wide gene expression data collected from the same human liver tissues and high-resolution promoter-focused chromatin interaction maps collected from human liver-derived HepG2 cells. We demonstrate widespread functional consequences of natural genetic variation on putative regulatory element activity and gene expression levels. Leveraging these extensive datasets, we fine-map a total of 77 GWAS loci that have been associated with at least one complex phenotype. Our results contribute to the repertoire of genes and regulatory mechanisms governing complex disease development and further the basic understanding of genetic and epigenetic regulation of gene expression in the human liver tissue.

Leveraging epigenomics and contactomics data to investigate SNP pairs in GWAS

Although Genome Wide Association Studies (GWAS) have led to many valuable insights into the genetic bases of common diseases over the past decade, the issue of missing heritability has surfaced, as the discovered main effect genetic variants found to date do not account for much of a trait’s predicted genetic component. We present a workflow, integrating epigenomics and topologically associating domain data, aimed at discovering trait-associated SNP pairs from GWAS where neither SNP achieved independent genome-wide significance. Each analyzed SNP pair consists of one SNP in a putative active enhancer and another SNP in a putative physically interacting gene promoter in a trait-relevant tissue. As a proof-of-principle case study, we used this approach to identify focused collections of SNP pairs that we analyzed in three independent Type 2 diabetes (T2D) GWAS. This approach led us to discover 35 significant SNP pairs, encompassing both novel signals and signals for which we have found orthogonal support from other sources. Nine of these pairs are consistent with eQTL results, two are consistent with our own capture C experiments, and seven involve signals supported by recent T2D literature.

Functional Characterization of the Morpheus Gene Family

DATA ACCESS The cDNA sequences reported in this paper have been deposited in the GenBank database (accession numbers): KF175165-KF175225 and BACs accession numbers that are used in this study: AC148621, AC190226, AC097327, AC097332, AC190226, AC097333, AC145401, AC187943, AC166855, AC166597, AC167295, AC235773, AC202644, and AC234805. ABSTRACT The burst of segmental duplications during human and great ape evolution focuses on a set of “core” duplicons encoding great-ape-specific gene families. Characterization of these gene families is complicated by their high copy number, incomplete sequence, and polymorphic nature. We investigate the structure, transcriptional diversity, and protein localization of the nuclear pore complex-interacting protein (NPIP) or Morpheus gene family. The corresponding core, LCRA, encodes one of the most rapidly evolving genes in the human genome; LCRA has expanded to ~20 copies from a single ancestral locus in Old World monkey and is associated with most of the recurrent chromosome 16 microdeletions implicated in autism and mental retardation. Phylogenetic analysis and cDNA sequencing suggest two distinct subfamilies or subtypes, NPIPA and NPIPB. The latter expanded recently within the great apes due to a series of structural changes within the canonical gene structure. Among Old World monkey, we observe a testis-specific pattern of expression that contrasts with the ubiquitous pattern observed among human tissues. This change in the expression profile coincides with the structural events that reshaped the structure and organization of the gene family. Most of the expressed human copies are capable of producing an open reading frame. Immunofluorescence analyses of the morpheus genes showed a primary localization to both the nucleus and its periphery. We show that morpheus genes may be upregulated upon pI:C treatment and find evidence of human autoantibodies produced against the NPIPB protein, raising the possibility that morpheus genes may be related to immune- or autoimmune-related function.