TherMos: Estimating protein-DNA binding energies from in vivo binding profiles

Accurately characterizing transcription factor (TF)-DNA affinity is a central goal of regulatory genomics. Although thermodynamics provides the most natural language for describing the continuous range of TF-DNA affinity, traditional motif discovery algorithms focus instead on classification paradigms that aim to discriminate 'bound' and 'unbound' sequences. Moreover, these algorithms do not directly model the distribution of tags in ChIP-seq data. Here, we present a new algorithm named Thermodynamic Modeling of ChIP-seq (TherMos), which directly estimates a positionspecific binding energy matrix (PSEM) from ChIPseq/exo tag profiles. In cross-validation tests on seven genome-wide TF-DNA binding profiles, one of which we generated via ChIP-seq on a complex developing tissue, TherMos predicted quantitative TF-DNA binding with greater accuracy than five well-known algorithms. We experimentally validated TherMos binding energy models for Klf4 and Esrrb, using a novel protocol to measure PSEMs in vitro. Strikingly, our measurements revealed strong nonadditivity at multiple positions within the two PSEMs. Among the algorithms tested, only TherMos was able to model the entire binding energy landscape of Klf4 and Esrrb. Our study reveals new insights into the energetics of TF-DNA binding in vivo and provides an accurate first-principles approach to binding energy inference from ChIP-seq and ChIP-exo data. © 2013 The Author(s).Link_to_subscribed_fulltex

Deciphering the Sox-Oct partner code by quantitative cooperativity measurements

Several Sox-Oct transcription factor (TF) combinations have been shown to cooperate on diverse enhancers to determine cell fates. Here, we developed a method to quantify biochemically the Sox-Oct cooperation and assessed the pairing of the high-mobility group (HMG) domains of 11 Sox TFs with Oct4 on a series of composite DNA elements. This way, we clustered Sox proteins according to their dimerization preferences illustrating that Sox HMG domains evolved different propensities to cooperate with Oct4. Sox2, Sox14, Sox21 and Sox15 strongly cooperate on the canonical element but compete with Oct4 on a recently discovered compressed element. Sry also cooperates on the canonical element but binds additively to the compressed element. In contrast, Sox17 and Sox4 cooperate more strongly on the compressed than on the canonical element. Sox5 and Sox18 show some cooperation on both elements, whereas Sox8 and Sox9 compete on both elements. Testing rationally mutated Sox proteins combined with structural modeling highlights critical amino acids for differential Sox-Oct4 partnerships and demonstrates that the cooperativity correlates with the efficiency in producing induced pluripotent stem cells. Our results suggest selective Sox-Oct partnerships in genome regulation and provide a toolset to study protein cooperation on DNA

Association of MC1R Variants and host phenotypes with melanoma risk in CDKN2A mutation carriers: a GenoMEL study

<p><b>Background</b> Carrying the cyclin-dependent kinase inhibitor 2A (CDKN2A) germline mutations is associated with a high risk for melanoma. Penetrance of CDKN2A mutations is modified by pigmentation characteristics, nevus phenotypes, and some variants of the melanocortin-1 receptor gene (MC1R), which is known to have a role in the pigmentation process. However, investigation of the associations of both MC1R variants and host phenotypes with melanoma risk has been limited.</p> <p><b>Methods</b> We included 815 CDKN2A mutation carriers (473 affected, and 342 unaffected, with melanoma) from 186 families from 15 centers in Europe, North America, and Australia who participated in the Melanoma Genetics Consortium. In this family-based study, we assessed the associations of the four most frequent MC1R variants (V60L, V92M, R151C, and R160W) and the number of variants (1, ≥2 variants), alone or jointly with the host phenotypes (hair color, propensity to sunburn, and number of nevi), with melanoma risk in CDKN2A mutation carriers. These associations were estimated and tested using generalized estimating equations. All statistical tests were two-sided.</p> <p><b>Results</b> Carrying any one of the four most frequent MC1R variants (V60L, V92M, R151C, R160W) in CDKN2A mutation carriers was associated with a statistically significantly increased risk for melanoma across all continents (1.24 × 10−6 ≤ P ≤ .0007). A consistent pattern of increase in melanoma risk was also associated with increase in number of MC1R variants. The risk of melanoma associated with at least two MC1R variants was 2.6-fold higher than the risk associated with only one variant (odds ratio = 5.83 [95% confidence interval = 3.60 to 9.46] vs 2.25 [95% confidence interval = 1.44 to 3.52]; Ptrend = 1.86 × 10−8). The joint analysis of MC1R variants and host phenotypes showed statistically significant associations of melanoma risk, together with MC1R variants (.0001 ≤ P ≤ .04), hair color (.006 ≤ P ≤ .06), and number of nevi (6.9 × 10−6 ≤ P ≤ .02).</p> <p><b>Conclusion</b> Results show that MC1R variants, hair color, and number of nevi were jointly associated with melanoma risk in CDKN2A mutation carriers. This joint association may have important consequences for risk assessments in familial settings.</p&gt

Pathway-Based Analysis of a Melanoma Genome-Wide Association Study: Analysis of Genes Related to Tumour-Immunosuppression

Systemic immunosuppression is a risk factor for melanoma, and sunburn-induced immunosuppression is thought to be causal. Genes in immunosuppression pathways are therefore candidate melanoma-susceptibility genes. If variants within these genes individually have a small effect on disease risk, the association may be undetected in genome-wide association (GWA) studies due to low power to reach a high significance level. Pathway-based approaches have been suggested as a method of incorporating a priori knowledge into the analysis of GWA studies. In this study, the association of 1113 single nucleotide polymorphisms (SNPs) in 43 genes (39 genomic regions) related to immunosuppression have been analysed using a gene-set approach in 1539 melanoma cases and 3917 controls from the GenoMEL consortium GWA study. The association between melanoma susceptibility and the whole set of tumour-immunosuppression genes, and also predefined functional subgroups of genes, was considered. The analysis was based on a measure formed by summing the evidence from the most significant SNP in each gene, and significance was evaluated empirically by case-control label permutation. An association was found between melanoma and the complete set of genes (pemp = 0.002), as well as the subgroups related to the generation of tolerogenic dendritic cells (pemp = 0.006) and secretion of suppressive factors (pemp = 0.0004), thus providing preliminary evidence of involvement of tumour-immunosuppression gene polymorphisms in melanoma susceptibility. The analysis was repeated on a second phase of the GenoMEL study, which showed no evidence of an association. As one of the first attempts to replicate a pathway-level association, our results suggest that low power and heterogeneity may present challenges

Loss of the scavenger mRNA decapping enzyme DCPS causes syndromic intellectual disability with neuromuscular defects

mRNA decay is an essential and active process that allows cells to continuously adapt gene expression to internal and environmental cues. There are two mRNA degradation pathways: 3' to 5' and 5' to 3'. The DCPS protein is the scavenger mRNA decapping enzyme which functions in the last step of the 3' end mRNA decay pathway. We have identified a DCPS pathogenic mutation in a large family with three affected individuals presenting with a novel recessive syndrome consisting of craniofacial anomalies, intellectual disability and neuromuscular defects. Using patient's primary cells, we show that this homozygous splice mutation results in a DCPS loss-of-function allele. Diagnostic biochemical analyses using various m(7)G cap derivatives as substrates reveal no DCPS enzymatic activity in patient's cells. Our results implicate DCPS and more generally RNA catabolism, as a critical cellular process for neurological development, normal cognition and organismal homeostasis in human

The crystal structure of the Sox4 HMG domain-DNA complex suggests a mechanism for positional interdependence in DNA recognition

10.1042/BJ20111768Biochemical Journal443139-47BIJO

A unique Oct4 interface is crucial for reprogramming to pluripotency

Terminally differentiated cells can be reprogrammed to pluripotency by the forced expression of Oct4, Sox2, Klf4 and c-Myc. However, it remains unknown how this leads to the multitude of epigenetic changes observed during the reprogramming process. Interestingly, Oct4 is the only factor that cannot be replaced by other members of the same family to induce pluripotency. To understand the unique role of Oct4 in reprogramming, we determined the structure of its POU domain bound to DNA. We show that the linker between the two DNA-binding domains is structured as an α-helix and exposed to the protein's surface, in contrast to the unstructured linker of Oct1. Point mutations in this α-helix alter or abolish the reprogramming activity of Oct4, but do not affect its other fundamental properties. On the basis of mass spectrometry studies of the interactome of wild-type and mutant Oct4, we propose that the linker functions as a protein-protein interaction interface and plays a crucial role during reprogramming by recruiting key epigenetic players to Oct4 target genes. Thus, we provide molecular insights to explain how Oct4 contributes to the reprogramming process

A variant in FTO shows association with melanoma risk not due to BMI

We report the results of an association study of melanoma that is based on the genome-wide imputation of the genotypes of 1,353 cases and 3,566 controls of European origin conducted by the GenoMEL consortium. This revealed an association between several SNPs in intron 8 of the FTO gene, including rs16953002, which replicated using 12,313 cases and 55,667 controls of European ancestry from Europe, the USA and Australia (combined P = 3.6 7 10(-12), per-allele odds ratio for allele A = 1.16). In addition to identifying a new melanoma-susceptibility locus, this is to our knowledge the first study to identify and replicate an association with SNPs in FTO not related to body mass index (BMI). These SNPs are not in intron 1 (the BMI-related region) and exhibit no association with BMI. This suggests FTO's function may be broader than the existing paradigm that FTO variants influence multiple traits only through their associations with BMI and obesity

Genome-wide association study identifies novel loci predisposing to cutaneous melanoma

We performed a multistage genome-wide association study of melanoma. In a discovery cohort of 1804 melanoma cases and 1026 controls, we identified loci at chromosomes 15q13.1 (HERC2/OCA2 region) and 16q24.3 (MC1R) regions that reached genome-wide significance within this study and also found strong evidence for genetic effects on susceptibility to melanoma from markers on chromosome 9p21.3 in the p16/ARF region and on chromosome 1q21.3 (ARNT/LASS2/ANXA9 region). The most significant single-nucleotide polymorphisms (SNPs) in the 15q13.1 locus (rs1129038 and rs12913832) lie within a genomic region that has profound effects on eye and skin color; notably, 50% of variability in eye color is associated with variation in the SNP rs12913832. Because eye and skin colors vary across European populations, we further evaluated the associations of the significant SNPs after carefully adjusting for European substructure. We also evaluated the top 10 most significant SNPs by using data from three other genome-wide scans. Additional in silico data provided replication of the findings from the most significant region on chromosome 1q21.3 rs7412746 (P = 6 x 10(-10)). Together, these data identified several candidate genes for additional studies to identify causal variants predisposing to increased risk for developing melanoma