Individual-specific levels of CTG•CAG somatic instability are shared across multiple tissues in myotonic dystrophy type 1.

Myotonic dystrophy type 1 is a complex disease caused by a genetically unstable CTG repeat expansion in the 3'-untranslated region of the DMPK gene. Age-dependent, tissue-specific somatic instability has confounded genotype-phenotype associations, but growing evidence suggests that it also contributes directly toward disease progression. Using a well-characterized clinical cohort of DM1 patients from Costa Rica, we quantified somatic instability in blood, buccal cells, skin and skeletal muscle. Whilst skeletal muscle showed the largest expansions, modal allele lengths in skin were also very large and frequently exceeded 2000 CTG repeats. Similarly, the degree of somatic expansion in blood, muscle and skin were associated with each other. Notably, we found that the degree of somatic expansion in skin was highly predictive of that in skeletal muscle. More importantly, we established that individuals whose repeat expanded more rapidly than expected in one tissue (after correction for progenitor allele length and age) also expanded more rapidly than expected in other tissues. We also provide evidence suggesting that individuals in whom the repeat expanded more rapidly than expected in skeletal muscle have an earlier age at onset than expected (after correction for the progenitor allele length). Pyrosequencing analyses of the genomic DNA flanking the CTG repeat revealed that the degree of methylation in muscle was well predicted by the muscle modal allele length and age, but that neither methylation of the flanking DNA nor levels of DMPK sense and anti-sense transcripts could obviously explain individual- or tissue-specific patterns of somatic instability.

PI3K/AKT signaling allows for MAPK/ERK pathway independency mediating dedifferentiation-driven treatment resistance in melanoma.

BACKGROUND: Current therapeutic management of advanced melanoma patients largely depends on their BRAF mutation status. However, the vast heterogeneity of the tumors hampers the success of therapies targeting the MAPK/ERK pathway alone. Dissecting this heterogeneity will contribute to identifying key players in the oncogenic progression to tailor more effective therapies. METHODS: We performed a comprehensive molecular and phenotypic characterization of a panel of patient-derived BRAFV600E-positive melanoma cell lines. Transcriptional profiling was used to identify groups of coregulated genes whose expression relates to an increased migratory potential and a higher resistance. RESULTS: A decrease in sensitivity to MAPK/ERK pathway inhibition with vemurafenib or trametinib corresponded with an increasing quiescence and migratory properties of the cells. This was accompanied by the loss of transcriptional signatures of melanocytic differentiation, and the gain of stem cell features that conferred highly-resistant/mesenchymal-like cells with increased xenobiotic efflux capacity. Nevertheless, targeting of the implicated ABC transporters did not improve the response to vemurafenib, indicating that incomplete BRAF inhibition due to reduced drug uptake is not a main driver of resistance. Rather, indifference to MAPK/ERK pathway inhibition arose from the activation of compensatory signaling cascades. The PI3K/AKT pathway in particular showed a higher activity in mesenchymal-like cells, conferring a lower dependency on MAPK/ERK signaling and supporting stem-like properties that could be reverted by dual PI3K/mTOR inhibition with dactolisib. CONCLUSIONS: In case of MAPK/ERK independency, therapeutic focus may be shifted to the PI3K/AKT pathway to overcome late-stage resistance in melanoma tumors that have acquired a mesenchymal phenotype. Video Abstract.

Dynamic transcriptome analysis reveals signatures of paradoxical effect of vemurafenib on human dermal fibroblasts.

BACKGROUND: Vemurafenib (PLX4032) is one of the most frequently used treatments for late-stage melanoma patients with the BRAFV600E mutation; however, acquired resistance to the drug poses as a major challenge. It remains to be determined whether off-target effects of vemurafenib on normal stroma components could reshape the tumor microenvironment in a way that contributes to cancer progression and drug resistance. METHODS: By using temporally-resolved RNA- and ATAC-seq, we studied the early molecular changes induced by vemurafenib in human dermal fibroblast (HDF), a main stromal component in melanoma and other tumors with high prevalence of BRAFV600 mutations. RESULTS: Transcriptomics analyses revealed a stepwise up-regulation of proliferation signatures, together with a down-regulation of autophagy and proteolytic processes. The gene expression changes in HDF strongly correlated in an inverse way with those in BRAFV600E mutant malignant melanoma (MaMel) cell lines, consistent with the observation of a paradoxical effect of vemurafenib, leading to hyperphosphorylation of MEK1/2 and ERK1/2. The transcriptional changes in HDF were not strongly determined by alterations in chromatin accessibility; rather, an already permissive chromatin landscape seemed to facilitate the early accessibility to MAPK/ERK-regulated transcription factor binding sites. Combinatorial treatment with the MEK inhibitor trametinib did not preclude the paradoxical activation of MAPK/ERK signaling in HDF. When administered together, vemurafenib partially compensated for the reduction of cell viability and proliferation induced by trametinib. These paradoxical changes were restrained by using the third generation BRAF inhibitor PLX8394, a so-called paradox breaker compound. However, the advantageous effects on HDF during combination therapies were also lost. CONCLUSIONS: Vemurafenib induces paradoxical changes in HDF, enabled by a permissive chromatin landscape. These changes might provide an advantage during combination therapies, by compensating for the toxicity induced in stromal cells by less specific MAPK/ERK inhibitors. Our results highlight the relevance of evaluating the effects of the drugs on non-transformed stromal components, carefully considering the implications of their administration either as mono- or combination therapies. Video Abstract.

Myotonic dystrophy type 1 (DM1) clinical subtypes and CTCF site methylation status flanking the CTG expansion are mutant allele length-dependent.

Myotonic dystrophy type 1 (DM1) is a complex disease with a wide spectrum of symptoms. The exact relationship between mutant CTG repeat expansion size and clinical outcome remains unclear. DM1 congenital patients (CDM) inherit the largest expanded alleles, which are associated with abnormal and increased DNA methylation flanking the CTG repeat. However, DNA methylation at the DMPK locus remains understudied. Its relationship to DM1 clinical subtypes, expansion size and age-at-onset is not yet completely understood. Using pyrosequencing-based methylation analysis on 225 blood DNA samples from Costa Rican DM1 patients, we determined that the size of the estimated progenitor allele length (ePAL) is not only a good discriminator between CDM and non-CDM cases (with an estimated threshold at 653 CTG repeats), but also for all DM1 clinical subtypes. Secondly, increased methylation at both CTCF sites upstream and downstream of the expansion was almost exclusively present in CDM cases. Thirdly, levels of abnormal methylation were associated with clinical subtype, age and ePAL, with strong correlations between these variables. Fourthly, both ePAL and the intergenerational expansion size were significantly associated with methylation status. Finally, methylation status was associated with ePAL and maternal inheritance, with almost exclusively maternal transmission of CDM. In conclusion, increased DNA methylation at the CTCF sites flanking the DM1 expansion could be linked to ePAL, and both increased methylation and the ePAL could be considered biomarkers for the CDM phenotype.

Longitudinal increases in somatic mosaicism of the expanded CTG repeat in myotonic dystrophy type 1 are associated with variation in age-at-onset.

In myotonic dystrophy type 1 (DM1), somatic mosaicism of the (CTG)n repeat expansion is age-dependent, tissue-specific and expansion-biased. These features contribute toward variation in disease severity and confound genotype-to-phenotype analyses. To investigate how the (CTG)n repeat expansion changes over time, we collected three longitudinal blood DNA samples separated by 8-15 years and used small pool and single-molecule PCR in 43 DM1 patients. We used the lower boundary of the allele length distribution as the best estimate for the inherited progenitor allele length (ePAL), which is itself the best predictor of disease severity. Although in most patients the lower boundary of the allele length distribution was conserved over time, in many this estimate also increased with age, suggesting samples for research studies and clinical trials should be obtained as early as possible. As expected, the modal allele length increased over time, driven primarily by ePAL, age-at-sampling and the time interval. As expected, small expansions <100 repeats did not expand as rapidly as larger alleles. However, the rate of expansion of very large alleles was not obviously proportionally higher. This may, at least in part, be a result of the allele length-dependent increase in large contractions that we also observed. We also determined that individual-specific variation in the increase of modal allele length over time not accounted for by ePAL, age-at-sampling and time was inversely associated with individual-specific variation in age-at-onset not accounted for by ePAL, further highlighting somatic expansion as a therapeutic target in DM1.

Analysis of mutational dynamics at the DMPK (CTG)n locus identifies saliva as a suitable DNA sample source for genetic analysis in myotonic dystrophy type 1.

Genotype-to-phenotype correlation studies in myotonic dystrophy type 1 (DM1) have been confounded by the age-dependent, tissue-specific and expansion-biased features of somatic mosaicism of the expanded CTG repeat. Previously, we showed that by controlling for the confounding effects of somatic instability to estimate the progenitor allele CTG length in blood DNA, age at onset correlations could be significantly improved. To determine the suitability of saliva DNA as a source for genotyping, we used small pool-PCR to perform a detailed quantitative study of the somatic mutational dynamics of the CTG repeat in saliva and blood DNA from 40 DM1 patients. Notably, the modal allele length in saliva was only moderately higher in saliva and not as large as previously observed in most other tissues. The lower boundary of the allele distribution was also slightly higher in saliva than it was in blood DNA. However, the progenitor allele length estimated in blood explained more of the variation in age at onset than that estimated from saliva. Interestingly, although the modal allele length was slightly higher in saliva, the overall degree of somatic variation was typically lower than in blood DNA, revealing new insights into the tissue-specific dynamics of somatic mosaicism. These data indicate that saliva constitutes an accessible, non-invasive and suitable DNA sample source for performing genetic studies in DM1.

Candidate gene study reveals DRD1 and DRD2 as putative interacting risk factors for youth depression.

Alterations in the monoaminergic neurotransmission systems are suspected to be involved in the etiology of neuropsychiatric disorders, including depression. The role of these pathways in the risk of developing depressive symptoms during childhood or adolescence is still not completely clear. This study sought to identify putative genetic factors in genes of serotonergic and dopaminergic systems modulating the level of manifestation of depressive symptoms in children and adolescents. We analyzed 170 single nucleotide polymorphisms (SNPs) in 21 candidate dopaminergic and serotonergic genes in a non-clinical sample of 410 Costa Rican participants of ages between 7 and 18 years, assessing the severity of depressive symptoms through the Child Depression Inventory (CDI). Genotypic and haplotypic associations, as well as epistatic effects, were examined. A significant interaction effect was detected between rs1039089 in conjunction with rs877138 located upstream of the dopamine D1 receptor (DRD1) and the dopamine D2 receptor (DRD2) genes respectively, although no evidence was found for any single variant or haplotype related to a differential liability. This newly described genetic interaction among putative regulatory regions of dopamine receptors could affect the level of manifestation of depressive symptoms through an imbalance of D1-D2 heteromers and modulation of cognitive processes.

A polymorphism in the MSH3 mismatch repair gene is associated with the levels of somatic instability of the expanded CTG repeat in the blood DNA of myotonic dystrophy type 1 patients.

Somatic mosaicism of the expanded CTG repeat in myotonic dystrophy type 1 is age-dependent, tissue-specific and expansion-biased, contributing toward the tissue-specificity and progressive nature of the symptoms. Previously, using regression modelling of repeat instability we showed that variation in the rate of somatic expansion in blood DNA contributes toward variation in age of onset, directly implicating somatic expansion in the disease pathway. Here, we confirm these results using a larger more genetically homogenous Costa Rican DM1 cohort (p<0.001). Interestingly, we also provide evidence that supports subtle sex-dependent differences in repeat length-dependent age at onset and somatic mutational dynamics. Previously, we demonstrated that variation in the rate of somatic expansion was a heritable quantitative trait. Given the important role that DNA mismatch repair genes play in mediating expansions in mouse models, we tested for modifier gene effects with 13 DNA mismatch gene polymorphisms (one each in MSH2, PMS2, MSH6 and MLH1; and nine in MSH3). After correcting for allele length and age effects, we identified three polymorphisms in MSH3 that were associated with variation in somatic instability: Rs26279 (p=0.003); Rs1677658 (p=0.009); and Rs10168 (p=0.031). However, only the association with Rs26279 remained significant after multiple testing correction. Although we revealed a statistically significant association between Rs26279 and somatic instability, we did not detect an association with the age at onset. Individuals with the A/A genotype for Rs26279 tended to show a greater propensity to expand the CTG repeat than other genotypes. Interestingly, this SNP results in an amino acid change in the critical ATPase domain of MSH3 and is potentially functionally dimorphic. These data suggest that MSH3 is a key player in generating somatic variation in DM1 patients and further highlight MSH3 as a potential therapeutic target.