Proteome-wide analysis of disease-associated SNPs that show allele-specific transcription factor binding

Transcription Factor-Wide Association Studies to Identify Functional SNPs in Alzheimer's Disease

Alzheimer's disease (AD) is a progressive neurodegenerative disorder with profound global impact. While genome-wide association studies (GWAS) have revealed genomic variants linked to AD, their translational impact has been limited due to challenges in interpreting the identified genetic associations. To address this challenge, we have devised a novel approach termed transcription factor-wide association studies (TF-WAS). By integrating the GWAS, expression quantitative trait loci, and transcriptome analyses, we selected 30 AD single nucleotide polymorphisms (SNPs) in noncoding regions that are likely to be functional. Using human transcription factor (TF) microarrays, we have identified 90 allele-specific TF interactions with 53 unique TFs. We then focused on several interactions involving SMAD4 and further validated them using electrophoretic mobility shift assay, luciferase, and chromatin immunoprecipitation on engineered genetic backgrounds (female cells). This approach holds promise for unraveling the intricacies of not just AD, but any complex disease with available GWAS data, providing insight into underlying molecular mechanisms and clues toward potential therapeutic targets.

High-throughput identification of functional regulatory SNPs in systemic lupus erythematosus

Genome-wide association studies implicate multiple loci in risk for systemic lupus erythematosus (SLE), but few contain exonic variants, rendering systematic identification of non-coding variants essential to decoding SLE genetics. We utilized SNP-seq and bioinformatic enrichment to interrogate 2180 single-nucleotide polymorphisms (SNPs) from 87 SLE risk loci for potential binding of transcription factors and related proteins from B cells. 52 SNPs that passed initial screening were tested by electrophoretic mobility shift and luciferase reporter assays. To validate the approach, we studied rs2297550 in detail, finding that the risk allele enhanced binding to the transcription factor Ikaros (encoded by IKZF1), thereby modulating expression of IKBKE. Correspondingly, primary cells from genotyped healthy donors bearing the risk allele expressed higher levels of the interferon / NF-κB regulator IKKε. Together, these findings define a set of likely functional non-coding lupus risk variants and identify a regulatory pathway involving rs2297550, Ikaros, and IKKε implicated by human genetics in risk for SLE.

Methods for Functional Characterization of Genetic Polymorphisms of Non-Coding Regulatory Regions of the Human Genome

Currently, numerous associations between genetic polymorphisms and various diseases have been characterized through the Genome-Wide Association Studies. Majority of the clinically significant polymorphisms are localized in non-coding regions of the genome. While modern bioinformatic resources make it possible to predict molecular mechanisms that explain influence of the non-coding polymorphisms on gene expression, such hypotheses require experimental verification. This review discusses the methods for elucidating molecular mechanisms underlying dependence of the disease pathogenesis on specific genetic variants within the non-coding sequences. A particular focus is on the methods for identification of transcription factors with binding efficiency dependent on polymorphic variations. Despite remarkable progress in bioinformatic resources enabling prediction of the impact of polymorphisms on the disease pathogenesis, there is still the need for experimental approaches to investigate this issue.

High-throughput identification of functional regulatory SNPs in systemic lupus erythematosus

Genome-wide association studies implicate multiple loci in risk for systemic lupus erythematosus (SLE), but few contain exonic variants, rendering systematic identification of non-coding variants essential to decoding SLE genetics. We utilized SNP-seq and bioinformatic enrichment to interrogate 2180 single-nucleotide polymorphisms (SNPs) from 87 SLE risk loci for potential binding of transcription factors and related proteins from B cells. 52 SNPs that passed initial screening were tested by electrophoretic mobility shift and luciferase reporter assays. To validate the approach, we studied rs2297550 in detail, finding that the risk allele enhanced binding to the transcription factor Ikaros (IKZF1), thereby modulating expression of IKBKE. Correspondingly, primary cells from genotyped healthy donors bearing the risk allele expressed higher levels of the interferon / NF-κB regulator IKKϵ. Together, these findings define a set of likely functional non-coding lupus risk variants and identify a new regulatory pathway involving rs2297550, Ikaros, and IKKϵ implicated by human genetics in risk for SLE.

Combined SNPs sequencing and allele specific proteomics capture reveal functional causality underpinning the 2p25 prostate cancer susceptibility locus

Genome wide association studies (GWASs) have identified numerous risk loci associated with prostate cancer, yet unraveling their functional significance remains elusive. Leveraging our high-throughput SNPs-seq method, we pinpointed rs4519489 within the multi-ancestry GWAS-discovered 2p25 locus as a potential functional SNP due to its significant allelic differences in protein binding. Here, we conduct a comprehensive analysis of rs4519489 and its associated gene, NOL10, employing diverse cohort data and experimental models. Clinical findings reveal a synergistic effect between rs4519489 genotype and NOL10 expression on prostate cancer prognosis and severity. Through unbiased proteomics screening, we reveal that the risk allele A of rs4519489 exhibits enhanced binding to USF1, a novel oncogenic transcription factor (TF) implicated in prostate cancer progression and prognosis, resulting in elevated NOL10 expression. Furthermore, we elucidate that NOL10 regulates cell cycle pathways, fostering prostate cancer progression. The concurrent expression of NOL10 and USF1 correlates with aggressive prostate cancer characteristics and poorer prognosis. Collectively, our study offers a robust strategy for functional SNP screening and TF identification through high-throughput SNPs-seq and unbiased proteomics, highlighting the rs4519489-USF1-NOL10 regulatory axis as a promising biomarker or therapeutic target for clinical diagnosis and treatment of prostate cancer.

E4F1 and ZNF148 are transcriptional activators of the -57A > C and wild-type TERT promoter

Point mutations within the TERT promoter are the most recurrent somatic noncoding mutations identified across different cancer types, including glioblastoma, melanoma, hepatocellular carcinoma, and bladder cancer. They are most abundant at -146C > T and -124C > T, and rarer at -57A > C, with the latter originally described as a familial case, but subsequently shown also to occur somatically. All three mutations create de novo E26-specific (ETS) binding sites and result in activation of the TERT gene, allowing cancer cells to achieve replicative immortality. Here, we used a systematic proteomics screen to identify transcription factors preferentially binding to the -146C > T, -124C > T, and -57A > C mutations. Although we confirmed binding of multiple ETS factors to the mutant -146C > T and -124C > T sequences, we identified E4F1 as a -57A > C-specific binder and ZNF148 as a TERT wild-type (WT) promoter binder that showed reduced interaction with the -124C > T allele. Both proteins are activating transcription factors that bind specifically to the -57A > C and WT (at position 124) TERT promoter sequence in corresponding cell lines, and up-regulate TERT transcription and telomerase activity. Our work describes new regulators of TERT gene expression with possible roles in cancer.

Combined CRISPRi and proteomics screening reveal a cohesin-CTCF-bound allele contributing to increased expression of RUVBL1 and prostate cancer progression

Genome-wide association studies along with expression quantitative trait locus (eQTL) mapping have identified hundreds of single-nucleotide polymorphisms (SNPs) and their target genes in prostate cancer (PCa), yet functional characterization of these risk loci remains challenging. To screen for potential regulatory SNPs, we designed a CRISPRi library containing 9,133 guide RNAs (gRNAs) to cover 2,166 candidate SNP loci implicated in PCa and identified 117 SNPs that could regulate 90 genes for PCa cell growth advantage. Among these, rs60464856 was covered by multiple gRNAs significantly depleted in screening (FDR < 0.05). Pooled SNP association analysis in the PRACTICAL and FinnGen cohorts showed significantly higher PCa risk for the rs60464856 G allele (p value = 1.2 × 10-16 and 3.2 × 10-7, respectively). Subsequent eQTL analysis revealed that the G allele is associated with increased RUVBL1 expression in multiple datasets. Further CRISPRi and xCas9 base editing confirmed that the rs60464856 G allele leads to elevated RUVBL1 expression. Furthermore, SILAC-based proteomic analysis demonstrated allelic binding of cohesin subunits at the rs60464856 region, where the HiC dataset showed consistent chromatin interactions in prostate cell lines. RUVBL1 depletion inhibited PCa cell proliferation and tumor growth in a xenograft mouse model. Gene-set enrichment analysis suggested an association of RUVBL1 expression with cell-cycle-related pathways. Increased expression of RUVBL1 and activation of cell-cycle pathways were correlated with poor PCa survival in TCGA datasets. Our CRISPRi screening prioritized about one hundred regulatory SNPs essential for prostate cell proliferation. In combination with proteomics and functional studies, we characterized the mechanistic role of rs60464856 and RUVBL1 in PCa progression.

Role of microRNAs in regulation of doxorubicin and paclitaxel responses in lung tumor cells

Lung cancer as the leading cause of cancer related mortality is always one of the main global health challenges. Despite the recent progresses in therapeutic methods, the mortality rate is still significantly high among lung cancer patients. A wide range of therapeutic methods including chemotherapy, radiotherapy, and surgery are used to treat lung cancer. Doxorubicin (DOX) and Paclitaxel (TXL) are widely used as the first-line chemotherapeutic drugs in lung cancer. However, there is a significant high percentage of DOX/TXL resistance in lung cancer patients, which leads to tumor recurrence and metastasis. Considering, the side effects of these drugs in normal tissues, it is required to clarify the molecular mechanisms of DOX/TXL resistance to introduce the efficient prognostic and therapeutic markers in lung cancer. MicroRNAs (miRNAs) have key roles in regulation of different pathophysiological processes including cell division, apoptosis, migration, and drug resistance. MiRNA deregulations are widely associated with chemo resistance in various cancers. Therefore, considering the importance of miRNAs in chemotherapy response, in the present review, we discussed the role of miRNAs in regulation of DOX/TXL response in lung cancer patients. It has been reported that miRNAs mainly induced DOX/TXL sensitivity in lung tumor cells by the regulation of signaling pathways, autophagy, transcription factors, and apoptosis. This review can be an effective step in introducing miRNAs as the non-invasive prognostic markers to predict DOX/TXL response in lung cancer patients.

Spatiotemporal and global profiling of DNA-protein interactions enables discovery of low-affinity transcription factors

Precise dissection of DNA-protein interactions is essential for elucidating the recognition basis, dynamics and gene regulation mechanism. However, global profiling of weak and dynamic DNA-protein interactions remains a long-standing challenge. Here, we establish the light-induced lysine (K) enabled crosslinking (LIKE-XL) strategy for spatiotemporal and global profiling of DNA-protein interactions. Harnessing unique abilities to capture weak and transient DNA-protein interactions, we demonstrate that LIKE-XL enables the discovery of low-affinity transcription-factor/DNA interactions via sequence-specific DNA baits, determining the binding sites for transcription factors that have been previously unknown. More importantly, we successfully decipher the dynamics of the transcription factor subproteome in response to drug treatment in a time-resolved manner, and find downstream target transcription factors from drug perturbations, providing insight into their dynamic transcriptional networks. The LIKE-XL strategy offers a complementary method to expand the DNA-protein profiling toolbox and map accurate DNA-protein interactomes that were previously inaccessible via non-covalent strategies, for better understanding of protein function in health and disease.

Combined CRISPRi and proteomics screening reveal a cohesin-CTCF-bound allele contributing to increased expression of RUVBL1 and prostate cancer progression

Genome-wide association studies along with expression quantitative trait loci (eQTL) mapping have identified hundreds of single nucleotide polymorphisms (SNPs) and their target genes in prostate cancer (PCa), yet functional characterization of these risk loci remains challenging. To screen for potential regulatory SNPs, we designed a CRISPRi library containing 9133 guide RNAs (gRNAs) to target 2,166 candidate SNP sites implicated in PCa and identified 117 SNPs that could regulate 90 genes for PCa cell growth advantage. Among these, rs60464856 was covered by multiple gRNAs significantly depleted in the screening (FDR<0.05). Pooled SNP association analysis in the PRACTICAL and FinnGen cohorts showed significantly higher PCa risk for the rs60464856 G allele (pvalue=1.2E-16 and 3.2E-7). Subsequent eQTL analysis revealed that the G allele is associated with increased RUVBL1 expression in multiple datasets. Further CRISPRi and xCas9 base editing proved the rs60464856 G allele leading to an elevated RUVBL1 expression. Furthermore, SILAC-based proteomic analysis demonstrated allelic binding of cohesin subunits at the rs60464856 region, where HiC dataset showed consistent chromatin interactions in prostate cell lines. RUVBL1 depletion inhibited PCa cell proliferation and tumor growth in xenograft mouse model. Gene set enrichment analysis suggested an association of RUVBL1 expression with cell-cycle-related pathways. An increased expression of RUVBL1 and activations of cell-cycle pathways were correlated with poor PCa survival in TCGA datasets. Together, our CRISPRi screening prioritized about one hundred regulatory SNPs essential for prostate cell proliferation. In combination with proteomics and functional studies, we characterized the mechanistic role of rs60464856 and RUVBL1 in PCa progression.

A leukemia-protective germline variant mediates chromatin module formation via transcription factor nucleation

Non-coding variants coordinate transcription factor (TF) binding and chromatin mark enrichment changes over regions spanning >100 kb. These molecularly coordinated regions are named "variable chromatin modules" (VCMs), providing a conceptual framework of how regulatory variation might shape complex traits. To better understand the molecular mechanisms underlying VCM formation, here, we mechanistically dissect a VCM-modulating noncoding variant that is associated with reduced chronic lymphocytic leukemia (CLL) predisposition and disease progression. This common, germline variant constitutes a 5-bp indel that controls the activity of an AXIN2 gene-linked VCM by creating a MEF2 binding site, which, upon binding, activates a super-enhancer-like regulatory element. This triggers a large change in TF binding activity and chromatin state at an enhancer cluster spanning >150 kb, coinciding with subtle, long-range chromatin compaction and robust AXIN2 up-regulation. Our results support a model in which the indel acts as an AXIN2 VCM-activating TF nucleation event, which modulates CLL pathology.

Functional annotation and investigation of the 10q24.33 melanoma risk locus identifies a common variant that influences transcriptional regulation of OBFC1

The 10q24.33 locus is known to be associated with susceptibility to cutaneous malignant melanoma (CMM), but the mechanisms underlying this association have been not extensively investigated. We carried out an integrative genomic analysis of 10q24.33 using epigenomic annotations and in vitro reporter gene assays to identify regulatory variants. We found two putative functional single nucleotide polymorphisms (SNPs) in an enhancer and in the promoter of OBFC1, respectively, in neural crest and CMM cells, one, rs2995264, altering enhancer activity. The minor allele G of rs2995264 correlated with lower OBFC1 expression in 470 CMM tumors and was confirmed to increase the CMM risk in a cohort of 484 CMM cases and 1801 controls of Italian origin. Hi-C and chromosome conformation capture (3C) experiments showed the interaction between the enhancer-SNP region and the promoter of OBFC1 and an isogenic model characterized by CRISPR-Cas9 deletion of the enhancer-SNP region confirmed the potential regulatory effect of rs2995264 on OBFC1 transcription. Moreover, the presence of G-rs2995264 risk allele reduced the binding affinity of the transcription factor MEOX2. Biologic investigations showed significant cell viability upon depletion of OBFC1, specifically in CMM cells that were homozygous for the protective allele. Clinically, high levels of OBFC1 expression associated with histologically favorable CMM tumors. Finally, preliminary results suggested the potential effect of decreased OBFC1 expression on telomerase activity in tumorigenic conditions. Our results support the hypothesis that reduced expression of OBFC1 gene through functional heritable DNA variation can contribute to malignant transformation of normal melanocytes.

CASCADE: high-throughput characterization of regulatory complex binding altered by non-coding variants

Non-coding DNA variants (NCVs) impact gene expression by altering binding sites for regulatory complexes. New high-throughput methods are needed to characterize the impact of NCVs on regulatory complexes. We developed CASCADE (Customizable Approach to Survey Complex Assembly at DNA Elements), an array-based high-throughput method to profile cofactor (COF) recruitment. CASCADE identifies DNA-bound transcription factor-cofactor (TF-COF) complexes in nuclear extracts and quantifies the impact of NCVs on their binding. We demonstrate CASCADE sensitivity in characterizing condition-specific recruitment of COFs p300 and RBBP5 (MLL subunit) to the CXCL10 promoter in lipopolysaccharide (LPS)-stimulated human macrophages and quantify the impact of all possible NCVs. To demonstrate applicability to NCV screens, we profile TF-COF binding to ~1,700 single-nucleotide polymorphism quantitative trait loci (SNP-QTLs) in human macrophages and identify perturbed ETS domain-containing complexes. CASCADE will facilitate high-throughput testing of molecular mechanisms of NCVs for diverse biological applications.

Bray et al. develop CASCADE, a method to profile transcription factor (TF)-cofactor (COF) complexes binding to DNA. They demonstrate the approach by profiling complex binding across the CXCL10 cytokine promoter and to ~1,700 single-nucleotide polymorphisms (SNPs). They anticipate that CASCADE can be applied to diverse biological systems to examine regulatory complex binding to DNA.

Identification of the transcription factor MAZ as a regulator of erythropoiesis

Erythropoiesis requires a combination of ubiquitous and tissue-specific transcription factors (TFs). Here, through DNA affinity purification followed by mass spectrometry, we have identified the widely expressed protein MAZ (Myc-associated zinc finger) as a TF that binds to the promoter of the erythroid-specific human α-globin gene. Genome-wide mapping in primary human erythroid cells revealed that MAZ also occupies active promoters as well as GATA1-bound enhancer elements of key erythroid genes. Consistent with an important role during erythropoiesis, knockdown of MAZ reduces α-globin expression in K562 cells and impairs differentiation in primary human erythroid cells. Genetic variants in the MAZ locus are associated with changes in clinically important human erythroid traits. Taken together, these findings reveal the zinc-finger TF MAZ to be a previously unrecognized regulator of the erythroid differentiation program.

Regulatory SNPs: Altered Transcription Factor Binding Sites Implicated in Complex Traits and Diseases

The vast majority of the genetic variants (mainly SNPs) associated with various human traits and diseases map to a noncoding part of the genome and are enriched in its regulatory compartment, suggesting that many causal variants may affect gene expression. The leading mechanism of action of these SNPs consists in the alterations in the transcription factor binding via creation or disruption of transcription factor binding sites (TFBSs) or some change in the affinity of these regulatory proteins to their cognate sites. In this review, we first focus on the history of the discovery of regulatory SNPs (rSNPs) and systematized description of the existing methodical approaches to their study. Then, we brief the recent comprehensive examples of rSNPs studied from the discovery of the changes in the TFBS sequence as a result of a nucleotide substitution to identification of its effect on the target gene expression and, eventually, to phenotype. We also describe state-of-the-art genome-wide approaches to identification of regulatory variants, including both making molecular sense of genome-wide association studies (GWAS) and the alternative approaches the primary goal of which is to determine the functionality of genetic variants. Among these approaches, special attention is paid to expression quantitative trait loci (eQTLs) analysis and the search for allele-specific events in RNA-seq (ASE events) as well as in ChIP-seq, DNase-seq, and ATAC-seq (ASB events) data.

A catalog of GWAS fine-mapping efforts in autoimmune disease

Genome-wide association studies (GWASs) have enabled unbiased identification of genetic loci contributing to common complex diseases. Because GWAS loci often harbor many variants and genes, it remains a major challenge to move from GWASs' statistical associations to the identification of causal variants and genes that underlie these association signals. Researchers have applied many statistical and functional fine-mapping strategies to prioritize genetic variants and genes as potential candidates. There is no gold standard in fine-mapping approaches, but consistent results across different approaches can improve confidence in the fine-mapping findings. Here, we combined text mining with a systematic review and formed a catalog of 85 studies with evidence of fine mapping for at least one autoimmune GWAS locus. Across all fine-mapping studies, we compiled 230 GWAS loci with allelic heterogeneity estimates and predictions of causal variants and trait-relevant genes. These 230 loci included 455 combinations of locus-by-disease association signals with 15 autoimmune diseases. Using these estimates, we assessed the probability of mediating disease risk associations across genes in GWAS loci and identified robust signals of causal disease biology. We predict that this comprehensive catalog of GWAS fine-mapping efforts in autoimmune disease will greatly help distill the plethora of information in the field and inform therapeutic strategies.

Transcription Factor AP4 Mediates Cell Fate Decisions: To Divide, Age, or Die

Simple Summary

Here, we review the literature on Activating Enhancer-Binding Protein 4 (AP4)/transcription factor AP4 (TFAP4) function and regulation and its role in cancer. Elevated expression of AP4 was detected in tumors of various organs and is associated with poor patient survival. AP4 is encoded by a Myc target gene and mediates cell fate decisions by regulating multiple processes, such as cell proliferation, epithelial-mesenchymal transition, stemness, apoptosis, and cellular senescence. Thereby, AP4 may be critical for tumor initiation and progression. In this review article, we summarize published evidence showing how AP4 functions as a transcriptional activator and repressor of a plethora of direct target genes in various physiological and pathological conditions. We also highlight the complex interactions of AP4 with c-Myc, N-Myc, p53, lncRNAs, and miRNAs in feed-back loops, which control AP4 levels and mediate AP4 functions. In the future, a better understanding of AP4 may contribute to improved prognosis and therapy of cancer.

Abstract

Activating Enhancer-Binding Protein 4 (AP4)/transcription factor AP4 (TFAP4) is a basic-helix-loop-helix-leucine-zipper transcription factor that was first identified as a protein bound to SV40 promoters more than 30 years ago. Almost 15 years later, AP4 was characterized as a target of the c-Myc transcription factor, which is the product of a prototypic oncogene that is activated in the majority of tumors. Interestingly, AP4 seems to represent a central hub downstream of c-Myc and N-Myc that mediates some of their functions, such as proliferation and epithelial-mesenchymal transition (EMT). Elevated AP4 expression is associated with progression of cancer and poor patient prognosis in multiple tumor types. Deletion of AP4 in mice points to roles of AP4 in the control of stemness, tumor initiation and adaptive immunity. Interestingly, ex vivo AP4 inactivation results in increased DNA damage, senescence, and apoptosis, which may be caused by defective cell cycle progression. Here, we will summarize the roles of AP4 as a transcriptional repressor and activator of target genes and the contribution of protein and non-coding RNAs encoded by these genes, in regulating the above mentioned processes. In addition, proteins interacting with or regulating AP4 and the cellular signaling pathways altered after AP4 dysregulation in tumor cells will be discussed.

Epigenetic fine-mapping: identification of causal mechanisms for autoimmunity

Genome-wide association studies (GWAS) have identified genetic susceptibility loci for a variety of autoimmune and inflammatory diseases. These studies confirm the fundamental genetic basis of individual autoimmune diseases, and also point to shared etiological mechanisms across the spectrum of autoimmunity. While hundreds of genetic loci have been implicated in autoimmune diseases, the translation of individual susceptibility loci into specific molecular mechanisms for individual diseases remains difficult. This review highlights recent advances in the genetics of autoimmune disease, and the emerging use of epigenetic techniques to identify pathogenic cell types and causal molecular mechanisms of autoimmunity.

SNPs in Sites for DNA Methylation, Transcription Factor Binding, and miRNA Targets Leading to Allele-Specific Gene Expression and Contributing to Complex Disease Risk: A Systematic Review

INTRODUCTION

The complex genetic diversity among human populations results from an assortment of factors acting at various sequential levels, including mutations, population migrations, genetic drift, and selection. Although there are a plethora of DNA sequence variations identified through genome-wide association studies (GWAS), the challenge remains to explain the mechanisms underlying interindividual phenotypic disparity accounting for disease susceptibility. Single nucleotide polymorphisms (SNPs) present in the sites for DNA methylation, transcription factor (TF) binding, or miRNA targets can alter the gene expression. The systematic review aimed to evaluate the complex crosstalk among SNPs, miRNAs, DNA methylation, and TFs for complex multifactorial disease risk.

METHODS

PubMed and Scopus databases were used from inception until May 15, 2019. Initially, screening of articles involved studies assessing the interaction of SNPs with TFs, DNA methylation, or miRNAs resulting in allele-specific gene expression in complex multifactorial diseases. We also included the studies which provided experimental validation of the interaction of SNPs with each of these factors. The results from various studies on multifactorial diseases were assessed.

RESULTS

A total of 11 articles for SNPs interacting with DNA methylation, 30 articles for SNPs interacting with TFs, and 11 articles for SNPs in miRNA binding sites were selected. The interactions of SNPs with epigenetic factors were found to be implicated in different types of cancers, autoimmune diseases, cardiovascular diseases, diabetes, and asthma.

CONCLUSION

The systematic review provides evidence for the interplay between genetic and epigenetic risk factors through allele-specific gene expression in various complex multifactorial diseases.

Mass spectrometry-based proteomic capture of proteins bound to the MACC1 promoter in colon cancer

MACC1 (metastasis associated in colon cancer 1) is a key driver that induces metastasis in colon cancer. However, the mechanisms by which MACC1 expression is transcriptionally regulated and the factors enriched at the MACC1 promoter remain largely unknown. The binding of proteins to specific DNA sites in the genome is a major determinant of genomic maintenance and the regulation of specific genes. The study herein utilized two methods to study the binding proteins of the MACC1 promoter region in colon cancer. Specifically, we adopted CRISPR-based chromatin affinity purification with mass spectrometry (CRISPR-ChAP-MS) and a biotin-streptavidin pulldown assay coupled with MS to identify the specific proteome bound to the MACC1 promoter in two colon cell lines with different metastatic potential. A total of 24 proteins were identified by CRISPR-ChAP-MS as binding to the MACC1 promoter, among which c-JUN was validated by ChIP-PCR. A total of 739 binding protein candidates were identified by biotin-streptavidin pulldown assays coupled with MS, of which HNF4G and PAX6 were validated and compared for their binding to the same promoter sites in the two cell lines. Our studies suggest distinctive proteomic factors associated with the MACC1 promoter in colon cells with different metastatic potential. The dynamic regulatory factors accumulated at the promoter of MACC1 may provide novel insights into the regulatory mechanisms of MACC1 transcription.

Atomic-resolution mapping of transcription factor-DNA interactions by femtosecond laser crosslinking and mass spectrometry

Transcription factors (TFs) regulate target genes by specific interactions with DNA sequences. Detecting and understanding these interactions at the molecular level is of fundamental importance in biological and clinical contexts. Crosslinking mass spectrometry is a powerful tool to assist the structure prediction of protein complexes but has been limited to the study of protein-protein and protein-RNA interactions. Here, we present a femtosecond laser-induced crosslinking mass spectrometry (fliX-MS) workflow, which allows the mapping of protein-DNA contacts at single nucleotide and up to single amino acid resolution. Applied to recombinant histone octamers, NF1, and TBP in complex with DNA, our method is highly specific for the mapping of DNA binding domains. Identified crosslinks are in close agreement with previous biochemical data on DNA binding and mostly fit known complex structures. Applying fliX-MS to cells identifies several bona fide crosslinks on DNA binding domains, paving the way for future large scale ex vivo experiments.

The Human Immunopeptidome Project: A Roadmap to Predict and Treat Immune Diseases

The science that investigates the ensembles of all peptides associated to human leukocyte antigen (HLA) molecules is termed "immunopeptidomics" and is typically driven by mass spectrometry (MS) technologies. Recent advances in MS technologies, neoantigen discovery and cancer immunotherapy have catalyzed the launch of the Human Immunopeptidome Project (HIPP) with the goal of providing a complete map of the human immunopeptidome and making the technology so robust that it will be available in every clinic. Here, we provide a long-term perspective of the field and we use this framework to explore how we think the completion of the HIPP will truly impact the society in the future. In this context, we introduce the concept of immunopeptidome-wide association studies (IWAS). We highlight the importance of large cohort studies for the future and how applying quantitative immunopeptidomics at population scale may provide a new look at individual predisposition to common immune diseases as well as responsiveness to vaccines and immunotherapies. Through this vision, we aim to provide a fresh view of the field to stimulate new discussions within the community, and present what we see as the key challenges for the future for unlocking the full potential of immunopeptidomics in this era of precision medicine.

Transcriptional Regulation Analysis of Alzheimer's Disease Based on FastNCA Algorithm

Background:

Understanding the relationship between genetic variation and gene expression is a central issue in genetics. Although many studies have identified genetic variations associated with gene expression, it is unclear how they perturb the underlying regulatory network of gene expression.

Objective:

To explore how genetic variations perturb potential transcriptional regulation networks of Alzheimer’s disease (AD) to paint a more complete picture of the complex landscape of transcription regulation.

Methods:

Fast network component analysis (FastNCA), which can capture the genetic variations in the form of single nucleotide polymorphisms (SNPs), is applied to analyse the expression activities of TFs and their regulatory strengths on TGs using microarray and RNA-seq data of AD. Then, multi-data fusion analysis was used to analyze the different TGs regulated by the same TFs in the different data by constructing the transcriptional regulatory networks of differentially expressed genes.

Results:

the common TF regulating TGs are not necessarily identical in different data, they may be involved in the same pathways that are closely related to the pathogenesis of AD, such as immune response, signal transduction and cytokine-cytokine receptor interaction pathways. Even if they are involved in different pathways, these pathways are also confirmed to have a potential link with AD.

Conclusion:

The study shows that the pathways of different TGs regulated by the same TFs in different data are all closely related to AD. Multi-data fusion analysis can form a certain complement to some extent and get more comprehensive results in the process of exploring the pathogenesis of AD.

The polymorphic variant rs1800734 influences methylation acquisition and allele-specific TFAP4 binding in the MLH1 promoter leading to differential mRNA expression

Expression of the mismatch repair gene MutL homolog 1 (MLH1) is silenced in a clinically important subgroup of sporadic colorectal cancers. These cancers exhibit hypermutability with microsatellite instability (MSI) and differ from microsatellite-stable (MSS) colorectal cancers in both prognosis and response to therapies. Loss of MLH1 is usually due to epigenetic silencing with associated promoter methylation; coding somatic mutations rarely occur. Here we use the presence of a colorectal cancer (CRC) risk variant (rs1800734) within the MLH1 promoter to investigate the poorly understood mechanisms of MLH1 promoter methylation and loss of expression. We confirm the association of rs1800734 with MSI+ but not MSS cancer risk in our own data and by meta-analysis. Using sensitive allele-specific detection methods, we demonstrate that MLH1 is the target gene for rs1800734 mediated cancer risk. In normal colon tissue, small allele-specific differences exist only in MLH1 promoter methylation, but not gene expression. In contrast, allele-specific differences in both MLH1 methylation and expression are present in MSI+ cancers. We show that MLH1 transcriptional repression is dependent on DNA methylation and can be reversed by a methylation inhibitor. The rs1800734 allele influences the rate of methylation loss and amount of re-expression. The transcription factor TFAP4 binds to the rs1800734 region but with much weaker binding to the risk than the protective allele. TFAP4 binding is absent on both alleles when promoter methylation is present. Thus we propose that TFAP4 binding shields the protective rs1800734 allele of the MLH1 promoter from BRAF induced DNA methylation more effectively than the risk allele.

The Atherosclerosis Risk Variant rs2107595 Mediates Allele-Specific Transcriptional Regulation of HDAC9 via E2F3 and Rb1

Background and Purpose- Genome-wide association studies have identified the HDAC9 (histone deacetylase 9) gene region as a major risk locus for atherosclerotic stroke and coronary artery disease in humans. Previous results suggest a role of altered HDAC9 expression levels as the underlying disease mechanism. rs2107595, the lead single nucleotide polymorphism for stroke and coronary artery disease resides in noncoding DNA and colocalizes with histone modification marks suggestive of enhancer elements. Methods- To determine the mechanisms by which genetic variation at rs2107595 regulates HDAC9 expression and thus vascular risk we employed targeted resequencing, proteome-wide search for allele-specific nuclear binding partners, chromatin immunoprecipitation, genome-editing, reporter assays, circularized chromosome conformation capture, and gain- and loss-of-function experiments in cultured human cell lines and primary immune cells. Results- Targeted resequencing of the HDAC9 locus in patients with atherosclerotic stroke and controls supported candidacy of rs2107595 as the causative single nucleotide polymorphism. A proteomic search for nuclear binding partners revealed preferential binding of the E2F3/TFDP1/Rb1 complex (E2F transcription factor 3/transcription factor Dp-1/Retinoblastoma 1) to the rs2107595 common allele, consistent with the disruption of an E2F3 consensus site by the risk allele. Gain- and loss-of-function studies showed a regulatory effect of E2F/Rb proteins on HDAC9 expression. Compared with the common allele, the rs2107595 risk allele exhibited higher transcriptional capacity in luciferase assays and was associated with higher HDAC9 mRNA levels in primary macrophages and genome-edited Jurkat cells. Circularized chromosome conformation capture revealed a genomic interaction of the rs2107595 region with the HDAC9 promoter, which was stronger for the common allele as was the in vivo interaction with E2F3 and Rb1 determined by chromatin immunoprecipitation. Gain-of-function experiments in isogenic Jurkat cells demonstrated a key role of E2F3 in mediating rs2107595-dependent transcriptional regulation of HDAC9. Conclusions- Collectively, our findings imply allele-specific transcriptional regulation of HDAC9 via E2F3 and Rb1 as a major mechanism mediating vascular risk at rs2107595.

MicroRNA-608 Promotes Apoptosis in Non-Small Cell Lung Cancer Cells Treated With Doxorubicin Through the Inhibition of TFAP4

Lung cancer is the most commonly diagnosed type of cancer and the leading cause of cancer-associated death worldwide. MicroRNAs (miRNAs) are non-coding single-stranded RNA molecules of ∼20–25 nucleotides in length. Single nucleotide polymorphisms are a class of genetic variation in the human genome, which when present in miRNA genes are associated with the risk of developing cancer. This study aimed to identify whether the miRNA (miR)-608 polymorphism rs4919510 influenced the incidence of lung cancer, and to explore the underlying mechanisms of miR-608 in the pathogenesis of the disease. A total of 37 patients with non-small cell lung cancer (NSCLC) were selected to determine the expression levels of miR-608; 96 NSCLC patients and 136 cancer-free healthy controls were recruited to determine the incidence of miR-608 rs4919510 in lung cancer patients. Additionally, the impact of miR-608 on the expression of predicted target genes, cell migration, viability, proliferation, and apoptosis was also assessed. We found that the presence of miR-608 rs4919510 did not affect the susceptibility of patients to NSCLC or the maturation of miR-608. miR-608 expression levels were found to be downregulated in NSCLC tissues. Furthermore, overexpression of miR-608 promoted doxorubicin-induced apoptosis in NSCLC cell lines A549 and HCC4006 by inhibiting the expression of transcription factor activating enhancer-binding protein 4 (TFAP4), and high expression levels of TFAP4 were observed in NSCLC tissues. Therefore, our results may provide valuable insights for the chemotherapeutical treatment of NSCLC.

Brain Banks Spur New Frontiers in Neuropsychiatric Research and Strategies for Analysis and Validation

Neuropsychiatric disorders affect hundreds of millions of patients and families worldwide. To decode the molecular framework of these diseases, many studies use human postmortem brain samples. These studies reveal brain-specific genetic and epigenetic patterns via high-throughput sequencing technologies. Identifying best practices for the collection of postmortem brain samples, analyzing such large amounts of sequencing data, and interpreting these results are critical to advance neuropsychiatry. We provide an overview of human brain banks worldwide, including progress in China, highlighting some well-known projects using human postmortem brain samples to understand molecular regulation in both normal brains and those with neuropsychiatric disorders. Finally, we discuss future research strategies, as well as state-of-the-art statistical and experimental methods that are drawn upon brain bank resources to improve our understanding of the agents of neuropsychiatric disorders.

Challenges and opportunities in stroke genetics

Stroke, ischaemic stroke and subtypes of ischaemic stroke display substantial heritability. When compared with related vascular conditions, the number of established risk loci reaching genome-wide significance for association with stroke is still in the lower range, particularly for aetiological stroke subtypes such as large artery atherosclerotic stroke or small vessel stroke. Nevertheless, for individual loci substantial progress has been made in determining the specific mechanisms mediating stroke risk. In this review, we present a roadmap for functional follow-up of common risk variants associated with stroke. First, we discuss in silico strategies for characterizing signals in non-coding regions and highlight databases providing information on quantitative trait loci for mRNA and protein expression, as well as methylation, focussing on those with presumed relevance for stroke. Next, we discuss experimental strategies for following up on non-coding risk variants and regions such as massively parallel reporter assays, proteome-wide association studies, and chromatin conformation capture (3C) assays. These and other approaches are relevant for gaining insight into the specific variants and mechanisms mediating genetic stroke risk. Finally, we discuss how genetic findings could influence clinical practice by adding to diagnostic algorithms and eventually improve treatment options for stroke.

Relationship between Multiple Sclerosis-Associated IL2RA Risk Allele Variants and Circulating T Cell Phenotypes in Healthy Genotype-Selected Controls

Single nucleotide polymorphisms (SNPs) in or near the IL2RA gene, that encodes the interleukin-2 (IL-2) receptor α (CD25), are associated with increased risk of immune-mediated diseases including multiple sclerosis (MS). We investigated how the MS-associated IL2RA SNPs rs2104286 and rs11256593 are associated with CD25 expression on T cells ex vivo by multiparameter flow cytometry in paired genotype-selected healthy controls. We observed that MS-associated IL2RA SNPs rs2104286 and rs11256593 are associated with expression of CD25 in CD4⁺ but not CD8⁺ T cells. In CD4⁺ T cells, carriers of the risk genotype had a reduced frequency of CD25⁺ T_FH1 cells (p = 0.001) and an increased frequency of CD25⁺ recent thymic emigrant cells (p = 0.006). Furthermore, carriers of the risk genotype had a reduced surface expression of CD25 in post-thymic expanded CD4⁺ T cells (CD31⁻CD45RA⁺), CD39⁺ T_Reg cells and in several non-follicular memory subsets. Our study found novel associations of MS-associated IL2RA SNPs on expression of CD25 in CD4⁺ T cell subsets. Insight into the associations of MS-associated IL2RA SNPs, as these new findings provide, offers a better understanding of CD25 variation in the immune system and can lead to new insights into how MS-associated SNPs contribute to development of MS.

High-throughput identification of noncoding functional SNPs via type IIS enzyme restriction

Genome-wide association studies (GWAS) have identified many disease-associated noncoding variants, but cannot distinguish functional single-nucleotide polymorphisms (fSNPs) from others that reside incidentally within risk loci. To address this challenge, we developed an unbiased high-throughput screen that employs type IIS enzymatic restriction to identify fSNPs that allelically modulate the binding of regulatory proteins. We coupled this approach, termed SNP-seq, with flanking restriction enhanced pulldown (FREP) to identify regulation of CD40 by three disease-associated fSNPs via four regulatory proteins, RBPJ, RSRC2 and FUBP-1/TRAP150. Applying this approach across 27 loci associated with juvenile idiopathic arthritis, we identified 148 candidate fSNPs, including two that regulate STAT4 via the regulatory proteins SATB2 and H1.2. Together, these findings establish the utility of tandem SNP-seq/FREP to bridge the gap between GWAS and disease mechanism.

A common intronic variant of PARP1 confers melanoma risk and mediates melanocyte growth via regulation of MITF

Previous genome-wide association studies have identified a melanoma-associated locus at 1q42.1 that encompasses a ∼100-kb region spanning the PARP1 gene. Expression quantitative trait locus (eQTL) analysis in multiple cell types of the melanocytic lineage consistently demonstrated that the 1q42.1 melanoma risk allele (rs3219090[G]) is correlated with higher PARP1 levels. In silico fine-mapping and functional validation identified a common intronic indel, rs144361550 (-/GGGCCC; r2 = 0.947 with rs3219090), as displaying allele-specific transcriptional activity. A proteomic screen identified RECQL as binding to rs144361550 in an allele-preferential manner. In human primary melanocytes, PARP1 promoted cell proliferation and rescued BRAFV600E-induced senescence phenotypes in a PARylation-independent manner. PARP1 also transformed TERT-immortalized melanocytes expressing BRAFV600E. PARP1-mediated senescence rescue was accompanied by transcriptional activation of the melanocyte-lineage survival oncogene MITF, highlighting a new role for PARP1 in melanomagenesis.

Allele-specific quantitative proteomics unravels molecular mechanisms modulated by cis-regulatory PPARG locus variation

Genome-wide association studies identified numerous disease risk loci. Delineating molecular mechanisms influenced by cis-regulatory variants is essential to understand gene regulation and ultimately disease pathophysiology. Combining bioinformatics and public domain chromatin information with quantitative proteomics supports prediction of cis-regulatory variants and enabled identification of allele-dependent binding of both, transcription factors and coregulators at the type 2 diabetes associated PPARG locus. We found rs7647481A nonrisk allele binding of Yin Yang 1 (YY1), confirmed by allele-specific chromatin immunoprecipitation in primary adipocytes. Quantitative proteomics also found the coregulator RING1 and YY1 binding protein (RYBP) whose mRNA levels correlate with improved insulin sensitivity in primary adipose cells carrying the rs7647481A nonrisk allele. Our findings support a concept with diverse cis-regulatory variants contributing to disease pathophysiology at one locus. Proteome-wide identification of both, transcription factors and coregulators, can profoundly improve understanding of mechanisms underlying genetic associations.

Noncoding Variants Functional Prioritization Methods Based on Predicted Regulatory Factor Binding Sites

BACKGROUNDS

With the advent of the post genomic era, the research for the genetic mechanism of the diseases has found to be increasingly depended on the studies of the genes, the gene-networks and gene-protein interaction networks. To explore gene expression and regulation, the researchers have carried out many studies on transcription factors and their binding sites (TFBSs). Based on the large amount of transcription factor binding sites predicting values in the deep learning models, further computation and analysis have been done to reveal the relationship between the gene mutation and the occurrence of the disease. It has been demonstrated that based on the deep learning methods, the performances of the prediction for the functions of the noncoding variants are outperforming than those of the conventional methods. The research on the prediction for functions of Single Nucleotide Polymorphisms (SNPs) is expected to uncover the mechanism of the gene mutation affection on traits and diseases of human beings.

RESULTS

We reviewed the conventional TFBSs identification methods from different perspectives. As for the deep learning methods to predict the TFBSs, we discussed the related problems, such as the raw data preprocessing, the structure design of the deep convolution neural network (CNN) and the model performance measure et al. And then we summarized the techniques that usually used in finding out the functional noncoding variants from de novo sequence.

CONCLUSION

Along with the rapid development of the high-throughout assays, more and more sample data and chromatin features would be conducive to improve the prediction accuracy of the deep convolution neural network for TFBSs identification. Meanwhile, getting more insights into the deep CNN framework itself has been proved useful for both the promotion on model performance and the development for more suitable design to sample data. Based on the feature values predicted by the deep CNN model, the prioritization model for functional noncoding variants would contribute to reveal the affection of gene mutation on the diseases.

PeakXus: comprehensive transcription factor binding site discovery from ChIP-Nexus and ChIP-Exo experiments

MOTIVATION

Transcription factor (TF) binding can be studied accurately in vivo with ChIP-exo and ChIP-Nexus experiments. Only fraction of TF binding mechanisms are yet fully understood and accurate knowledge of binding locations and patterns of TFs is key to understanding binding that is not explained by simple positional weight matrix models. ChIP-exo/Nexus experiments can also offer insight on the effect of single nucleotide polymorphism (SNP) at TF binding sites on expression of the target genes. This is an important mechanism of action for disease-causing SNPs at non-coding genomic regions.

RESULTS

We describe a peak caller PeakXus that is specifically designed to leverage the increased resolution of ChIP-exo/Nexus and developed with the aim of making as few assumptions of the data as possible to allow discoveries of novel binding patterns. We apply PeakXus to ChIP-Nexus and ChIP-exo experiments performed both in Homo sapiens and in Drosophila melanogaster cell lines. We show that PeakXus consistently finds more peaks overlapping with a TF-specific recognition sequence than published methods. As an application example we demonstrate how PeakXus can be coupled with unique molecular identifiers (UMIs) to measure the effect of a SNP overlapping with a TF binding site on the in vivo binding of the TF.

AVAILABILITY AND IMPLEMENTATION

Source code of PeakXus is available at https://github.com/hartonen/PeakXus

CONTACT

tuomo.hartonen@helsinki.fi or jussi.taipale@ki.se.

Genetic-epigenetic interactions in cis: a major focus in the post-GWAS era

Studies on genetic-epigenetic interactions, including the mapping of methylation quantitative trait loci (mQTLs) and haplotype-dependent allele-specific DNA methylation (hap-ASM), have become a major focus in the post-genome-wide-association-study (GWAS) era. Such maps can nominate regulatory sequence variants that underlie GWAS signals for common diseases, ranging from neuropsychiatric disorders to cancers. Conversely, mQTLs need to be filtered out when searching for non-genetic effects in epigenome-wide association studies (EWAS). Sequence variants in CCCTC-binding factor (CTCF) and transcription factor binding sites have been mechanistically linked to mQTLs and hap-ASM. Identifying these sites can point to disease-associated transcriptional pathways, with implications for targeted treatment and prevention.

Two separate effects contribute to regulatory T cell defect in systemic lupus erythematosus patients and their unaffected relatives

Forkhead box P3 (FoxP3)⁺ regulatory T cells (T_regs) are functionally deficient in systemic lupus erythematosus (SLE), characterized by reduced surface CD25 [the interleukin (IL)‐2 receptor alpha chain]. Low‐dose IL‐2 therapy is a promising current approach to correct this defect. To elucidate the origins of the SLE T_reg phenotype, we studied its role through developmentally defined regulatory T cell (T_reg) subsets in 45 SLE patients, 103 SLE‐unaffected first‐degree relatives and 61 unrelated healthy control subjects, and genetic association with the CD25‐encoding IL2RA locus. We identified two separate, uncorrelated effects contributing to T_reg CD25. (1) SLE patients and unaffected relatives remarkably shared CD25 reduction versus controls, particularly in the developmentally earliest CD4⁺FoxP3⁺CD45RO^–CD31⁺ recent thymic emigrant T_regs. This first component effect influenced the proportions of circulating CD4⁺FoxP3^highCD45RO⁺ activated T_regs. (2) In contrast, patients and unaffected relatives differed sharply in their activated T_reg CD25 state: while relatives as control subjects up‐regulated CD25 strongly in these cells during differentiation from naive T_regs, SLE patients specifically failed to do so. This CD25 up‐regulation depended upon IL2RA genetic variation and was related functionally to the proliferation of activated T_regs, but not to their circulating numbers. Both effects were found related to T cell IL‐2 production. Our results point to (1) a heritable, intrathymic mechanism responsible for reduced CD25 on early T_regs and decreased activation capacity in an extended risk population, which can be compensated by (2) functionally independent CD25 up‐regulation upon peripheral T_reg activation that is selectively deficient in patients. We expect that T_reg‐directed therapies can be monitored more effectively when taking this distinction into account.

Functional characterization of a multi-cancer risk locus on chr5p15.33 reveals regulation of TERT by ZNF148

Genome wide association studies (GWAS) have mapped multiple independent cancer susceptibility loci to chr5p15.33. Here, we show that fine-mapping of pancreatic and testicular cancer GWAS within one of these loci (Region 2 in CLPTM1L) focuses the signal to nine highly correlated SNPs. Of these, rs36115365-C associated with increased pancreatic and testicular but decreased lung cancer and melanoma risk, and exhibited preferred protein-binding and enhanced regulatory activity. Transcriptional gene silencing of this regulatory element repressed TERT expression in an allele-specific manner. Proteomic analysis identifies allele-preferred binding of Zinc finger protein 148 (ZNF148) to rs36115365-C, further supported by binding of purified recombinant ZNF148. Knockdown of ZNF148 results in reduced TERT expression, telomerase activity and telomere length. Our results indicate that the association with chr5p15.33-Region 2 may be explained by rs36115365, a variant influencing TERT expression via ZNF148 in a manner consistent with elevated TERT in carriers of the C allele.

Molecular mechanisms underlying noncoding risk variations in psychiatric genetic studies

Recent large-scale genetic approaches such as genome-wide association studies have allowed the identification of common genetic variations that contribute to risk architectures of psychiatric disorders. However, most of these susceptibility variants are located in noncoding genomic regions that usually span multiple genes. As a result, pinpointing the precise variant(s) and biological mechanisms accounting for the risk remains challenging. By reviewing recent progresses in genetics, functional genomics and neurobiology of psychiatric disorders, as well as gene expression analyses of brain tissues, here we propose a roadmap to characterize the roles of noncoding risk loci in the pathogenesis of psychiatric illnesses (that is, identifying the underlying molecular mechanisms explaining the genetic risk conferred by those genomic loci, and recognizing putative functional causative variants). This roadmap involves integration of transcriptomic data, epidemiological and bioinformatic methods, as well as in vitro and in vivo experimental approaches. These tools will promote the translation of genetic discoveries to physiological mechanisms, and ultimately guide the development of preventive, therapeutic and prognostic measures for psychiatric disorders.

Modulation of transcription factor binding and epigenetic regulation of the MLH1 CpG island and shore by polymorphism rs1800734 in colorectal cancer

The MLH1 promoter polymorphism rs1800734 is associated with MLH1 CpG island hypermethylation and expression loss in colorectal cancer (CRC). Conversely, variant rs1800734 is associated with MLH1 shore, but not island, hypomethylation in peripheral blood mononuclear cell DNA. To explore these distinct patterns, MLH1 CpG island and shore methylation was assessed in CRC cell lines stratified by rs1800734 genotype. Cell lines containing the variant A allele demonstrated MLH1 shore hypomethylation compared to wild type (GG). There was significant enrichment of transcription factor AP4 at the MLH1 promoter in GG and GA cell lines, but not the AA cell line, by chromatin immunoprecipitation studies. Preferential binding to the G allele was confirmed by sequencing in the GA cell line. The enhancer-associated histone modification H3K4me1 was enriched at the MLH1 shore; however, H3K27ac was not, indicating the shore is an inactive enhancer. These results demonstrate the role of variant rs1800734 in altering transcription factor binding as well as epigenetics at regions beyond the MLH1 CpG island in which it is located.

The dynamic DNA methylation landscape of the mutL homolog 1 shore is altered by MLH1-93G>A polymorphism in normal tissues and colorectal cancer

BACKGROUND

Colorectal cancers (CRCs) undergo distinct genetic and epigenetic alterations. Expression of mutL homolog 1 (MLH1), a mismatch repair gene that corrects DNA replication errors, is lost in up to 15% of sporadic tumours due to mutation or, more commonly, due to DNA methylation of its promoter CpG island. A single nucleotide polymorphism (SNP) in the CpG island of MLH1 (MLH1-93G>A or rs1800734) is associated with CpG island hypermethylation and decreased MLH1 expression in CRC tumours. Further, in peripheral blood mononuclear cell (PBMC) DNA of both CRC cases and non-cancer controls, the variant allele of rs1800734 is associated with hypomethylation at the MLH1 shore, a region upstream of its CpG island that is less dense in CpG sites.

RESULTS

To determine whether this genotype-epigenotype association is present in other tissue types, including colorectal tumours, we assessed DNA methylation in matched normal colorectal tissue, tumour, and PBMC DNA from 349 population-based CRC cases recruited from the Ontario Familial Colorectal Cancer Registry. Using the semi-quantitative real-time PCR-based MethyLight assay, MLH1 shore methylation was significantly higher in tumour tissue than normal colon or PBMCs (P < 0.01). When shore methylation levels were stratified by SNP genotype, normal colorectal DNA and PBMC DNA were significantly hypomethylated in association with variant SNP genotype (P < 0.05). However, this association was lost in tumour DNA. Among distinct stages of CRC, metastatic stage IV CRC tumours incurred significant hypomethylation compared to stage I-III cases, irrespective of genotype status. Shore methylation of MLH1 was not associated with MSI status or promoter CpG island hypermethylation, regardless of genotype. To confirm these results, bisulfite sequencing was performed in matched tumour and normal colorectal specimens from six CRC cases, including two cases per genotype (wildtype, heterozygous, and homozygous variant). Bisulfite sequencing results corroborated the methylation patterns found by MethyLight, with significant hypomethylation in normal colorectal tissue of variant SNP allele carriers.

CONCLUSIONS

These results indicate that the normal tissue types tested (colorectum and PBMC) experience dynamic genotype-associated epigenetic alterations at the MLH1 shore, whereas tumour DNA incurs aberrant hypermethylation compared to normal DNA.

The non-coding variant rs1800734 enhances DCLK3 expression through long-range interaction and promotes colorectal cancer progression

Genome-wide association studies have identified a great number of non-coding risk variants for colorectal cancer (CRC). To date, the majority of these variants have not been functionally studied. Identification of allele-specific transcription factor (TF) binding is of great importance to understand regulatory consequences of such variants. A recently developed proteome-wide analysis of disease-associated SNPs (PWAS) enables identification of TF-DNA interactions in an unbiased manner. Here we perform a large-scale PWAS study to comprehensively characterize TF-binding landscape that is associated with CRC, which identifies 731 allele-specific TF binding at 116 CRC risk loci. This screen identifies the A-allele of rs1800734 within the promoter region of MLH1 as perturbing the binding of TFAP4 and consequently increasing DCLK3 expression through a long-range interaction, which promotes cancer malignancy through enhancing expression of the genes related to epithelial-to-mesenchymal transition.

SDHD Promoter Mutations Ablate GABP Transcription Factor Binding in Melanoma

SDHD encodes subunit D of the succinate dehydrogenase complex, an integral membrane protein. Across cancer types, recurrent SDHD promoter mutations were reported to occur exclusively in melanomas, at a frequency of 4% to 5%. These mutations are predicted to disrupt consensus ETS transcription factor-binding sites and are correlated with both reduced SDHD gene expression and poor prognosis. However, the consequence of these mutations on SDHD expression in melanoma is still unclear. Here, we found that expression of SDHD in melanoma correlated with the expression of multiple ETS transcription factors, particularly in SDHD promoter wild-type samples. Consistent with the predicted loss of ETS transcription factor binding, we observed that recurrent hotspot mutations resulted in decreased luciferase activity in reporter assays. Furthermore, we demonstrated specific GABPA and GABPB1 binding to probes containing the wild-type promoter sequences, with binding disrupted by the SDHD hotspot promoter mutations in both quantitative mass spectrometry and band-shift experiments. Finally, using siRNA-mediated knockdown across multiple melanoma cell lines, we determined that loss of GABPA resulted in reduced SDHD expression at both RNA and protein levels. These data are consistent with a key role for GABPA/B1 as the critical ETS transcription factors deregulating SDHD expression in the context of highly recurrent promoter mutations in melanoma and warrant a detailed search for other recurrent promoter mutations that create or disrupt GABPA consensus sequences. Cancer Res; 77(7); 1649-61. ©2017 AACR.

Interleukin-2 receptor-α proximal promoter hypomethylation is associated with multiple sclerosis

Genetic studies have demonstrated association between single-nucleotide polymorphisms within the IL2RA (interleukin-2 receptor α-subunit) gene and risk of developing multiple sclerosis (MS); however, these variants do not have obvious functional consequences. DNA methylation is a source of genetic variation that could impact on autoimmune disease risk. We investigated DNA methylation of the IL2RA promoter in genomic DNA obtained from peripheral blood mononuclear cells and neural tissue using matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF) mass spectrometry. A differential methylation profile of IL2RA was identified, suggesting that IL2RA expression was regulated by DNA methylation. We extended our analysis of DNA methylation to peripheral blood mononuclear cell (PBMC) of MS cases and controls using MALDI-TOF and Illumina HumanMethylation450 arrays. Analyses of CpG sites within the proximal promoter of IL2RA in PBMC showed no differences between MS cases and controls despite an increase in IL2RA expression. In contrast, we inferred significant DNA methylation differences specific to particular leukocyte subtypes in MS cases compared with controls by deconvolution of the array data. The decrease in methylation in patients correlated with an increase in IL2RA expression in T cells from MS cases in comparison with controls. Our data suggest that differential methylation of the IL2RA promoter in T cells could be an important pathogenic mechanism in MS.

Translational bioinformatics in the era of real-time biomedical, health care and wellness data streams

Monitoring and modeling biomedical, health care and wellness data from individuals and converging data on a population scale have tremendous potential to improve understanding of the transition to the healthy state of human physiology to disease setting. Wellness monitoring devices and companion software applications capable of generating alerts and sharing data with health care providers or social networks are now available. The accessibility and clinical utility of such data for disease or wellness research are currently limited. Designing methods for streaming data capture, real-time data aggregation, machine learning, predictive analytics and visualization solutions to integrate wellness or health monitoring data elements with the electronic medical records (EMRs) maintained by health care providers permits better utilization. Integration of population-scale biomedical, health care and wellness data would help to stratify patients for active health management and to understand clinically asymptomatic patients and underlying illness trajectories. In this article, we discuss various health-monitoring devices, their ability to capture the unique state of health represented in a patient and their application in individualized diagnostics, prognosis, clinical or wellness intervention. We also discuss examples of translational bioinformatics approaches to integrating patient-generated data with existing EMRs, personal health records, patient portals and clinical data repositories. Briefly, translational bioinformatics methods, tools and resources are at the center of these advances in implementing real-time biomedical and health care analytics in the clinical setting. Furthermore, these advances are poised to play a significant role in clinical decision-making and implementation of data-driven medicine and wellness care.

A quantitative proteomics approach identifies ETV6 and IKZF1 as new regulators of an ERG-driven transcriptional network

Aberrant stem cell-like gene regulatory networks are a feature of leukaemogenesis. The ETS-related gene (ERG), an important regulator of normal haematopoiesis, is also highly expressed in T-ALL and acute myeloid leukaemia (AML). However, the transcriptional regulation of ERG in leukaemic cells remains poorly understood. In order to discover transcriptional regulators of ERG, we employed a quantitative mass spectrometry-based method to identify factors binding the 321 bp ERG +85 stem cell enhancer region in MOLT-4 T-ALL and KG-1 AML cells. Using this approach, we identified a number of known binders of the +85 enhancer in leukaemic cells along with previously unknown binders, including ETV6 and IKZF1. We confirmed that ETV6 and IKZF1 were also bound at the +85 enhancer in both leukaemic cells and in healthy human CD34+ haematopoietic stem and progenitor cells. Knockdown experiments confirmed that ETV6 and IKZF1 are transcriptional regulators not just of ERG, but also of a number of genes regulated by a densely interconnected network of seven transcription factors. At last, we show that ETV6 and IKZF1 expression levels are positively correlated with expression of a number of heptad genes in AML and high expression of all nine genes confers poorer overall prognosis.

Mining the Unknown: Assigning Function to Noncoding Single Nucleotide Polymorphisms

One of the formative goals of genetics research is to understand how genetic variation leads to phenotypic differences and human disease. Genome-wide association studies (GWASs) bring us closer to this goal by linking variation with disease faster than ever before. Despite this, GWASs alone are unable to pinpoint disease-causing single nucleotide polymorphisms (SNPs). Noncoding SNPs, which represent the majority of GWAS SNPs, present a particular challenge. To address this challenge, an array of computational tools designed to prioritize and predict the function of noncoding GWAS SNPs have been developed. However, fewer than 40% of GWAS publications from 2015 utilized these tools. We discuss several leading methods for annotating noncoding variants and how they can be integrated into research pipelines in hopes that they will be broadly applied in future GWAS analyses.

The mechanism of transactivation regulation due to polymorphic short tandem repeats (STRs) using IGF1 promoter as a model

Functional short tandem repeats (STR) are polymorphic in the population, and the number of repeats regulates the expression of nearby genes (known as expression STR, eSTR). STR in IGF1 promoter has been extensively studied for its association with IGF1 concentration in blood and various clinical traits and represents an important eSTR. We previously used an in-vitro luciferase reporter model to examine the interaction between STRs and SNPs in IGF1 promoter. Here, we further explored the mechanism how the number of repeats of the STR regulates gene transcription. An inverse correlation between the number of repeats and the extent of transactivation was found in a haplotype consisting of three promoter SNPs (C-STR-T-T). We showed that these adjacent SNPs located outside the STR were required for the STR to function as eSTR. The C allele of rs35767 provides a binding site for CCAAT/enhancer-binding-protein δ (C/EBPD), which is essential for the gradational transactivation property of eSTR and FOXA3 may also be involved. Therefore, we propose a mechanism in which the gradational transactivation by the eSTR is caused by the interaction of one or more transcriptional complexes located outside the STR, rather than by direct binding to a repeat motif of the STR.

Multiple single nucleotide polymorphisms in the first intron of the IL2RA gene affect transcription factor binding and enhancer activity

IL2RA gene encodes the alpha subunit of a high-affinity receptor for interleukin-2 which is expressed by several distinct populations of lymphocytes involved in autoimmune processes. A large number of polymorphic alleles of the IL2RA locus are associated with the development of various autoimmune diseases. With bioinformatics analysis we the dissected the first intron of the IL2RA gene and selected several single nucleotide polymorphisms (SNPs) that may influence the regulation of the IL2RA gene in cell types relevant to autoimmune pathology. We described five enhancers containing the selected SNPs that stimulated activity of the IL2RA promoter in a cell-type specific manner, and tested the effect of specific SNP alleles on activity of the respective enhancers (E1 to E5, labeled according to the distance to the promoter). The E4 enhancer with minor T variant of rs61839660 SNP demonstrated reduced activity due to disrupted binding of MEF2A/C transcription factors (TFs). Neither rs706778 nor rs706779 SNPs, both associated with a number of autoimmune diseases, had any effect on the activity of the enhancer E2. However, rare variants of several SNPs (rs139767239, rs115133228, rs12722502, rs12722635) genetically linked to either rs706778 and/or rs706779 significantly influenced the activity of E1, E3 and E5 enhancers, presumably by disrupting EBF1, GABPA and ELF1 binding sites.

RETRACTION: Illuminating the dark road from schizophrenia genetic associations to disease mechanisms

Recent large-scale genome-wide association studies (GWAS) have enabled the discovery of common genetic variations contributing to risk architectures of schizophrenia in human populations; however, the majority of GWAS-identified variants are located in large genomic regions spanning multiple genes, and recognizing the precise targets and mechanisms of these clinical associations is now the major challenge. Here, we review recent progress in schizophrenia genetics, functional genomics and related neuroscience research, and propose a functional pipeline to translate schizophrenia GWAS risk loci into disease biology and information for drug discovery. The pipeline includes identification of underlying molecular mechanisms using transcriptomic data in human brain, prioritization of putative functional causative variants by the integration of genetic epidemiological and bioinformatics methods as well as molecular approaches, and in vitro and in vivo experimental characterizations of the identified targeted species and causative variants to dissect the relevant disease biology. These approaches will accelerate progress from schizophrenia genetic studies to biological mechanisms and ultimately guide the development of prognostic, preventive and therapeutic measures.