Clonal evolution in tyrosine kinase inhibitor-resistance: lessons from in vitro-models.

Introduction: Resistance in anti-cancer treatment is a result of clonal evolution and clonal selection. In chronic myeloid leukemia (CML), the hematopoietic neoplasm is predominantly caused by the formation of the BCR::ABL1 kinase. Evidently, treatment with tyrosine kinase inhibitors (TKIs) is tremendously successful. It has become the role model of targeted therapy. However, therapy resistance to TKIs leads to loss of molecular remission in about 25% of CML patients being partially due to BCR::ABL1 kinase mutations, while for the remaining cases, various other mechanisms are discussed. Methods: Here, we established an in vitro-TKI resistance model against the TKIs imatinib and nilotinib and performed exome sequencing. Results: In this model, acquired sequence variants in NRAS, KRAS, PTPN11, and PDGFRB were identified in TKI resistance. The well-known pathogenic NRAS p.(Gln61Lys) variant provided a strong benefit for CML cells under TKI exposure visible by increased cell number (6.2-fold, p < 0.001) and decreased apoptosis (-25%, p < 0.001), proving the functionality of our approach. The transfection of PTPN11 p.(Tyr279Cys) led to increased cell number (1.7-fold, p = 0.03) and proliferation (2.0-fold, p < 0.001) under imatinib treatment. Discussion: Our data demonstrate that our in vitro-model can be used to study the effect of specific variants on TKI resistance and to identify new driver mutations and genes playing a role in TKI resistance. The established pipeline can be used to study candidates acquired in TKI-resistant patients, thereby providing new options for the development of new therapy strategies to overcome resistance.

CRISPR/Cas9-induced knockout reveals the role of ABCB1 in the response to temozolomide, carmustine and lomustine in glioblastoma multiforme.

Glioblastoma multiforme (GBM) is the most common malignant brain tumor with limited therapeutic options. Besides surgery, chemotherapy using temozolomide, carmustine or lomustine is the main pillar of therapy. However, therapy success is limited and prognosis still is very poor. One restraining factor is drug resistance caused by drug transporters of the ATP-binding cassette family, e.g. ABCB1 and ABCG2, located at the blood-brain barrier and on tumor cells. The active efflux of xenobiotics including drugs, e.g. temozolomide, leads to low intracellular drug concentrations and subsequently insufficient anti-tumor effects. Nevertheless, the role of efflux transporters in GBM is controversially discussed. In the present study, we analyzed the role of ABCB1 and ABCG2 in GBM cells showing that ABCB1, but marginally ABCG2, is relevant. Applying a CRISPR/Cas9-derived ABCB1 knockout, the response to temozolomide was significantly augmented demonstrated by decreased cell number (p < 0.001) and proliferation rate (p = 0.04), while apoptosis was increased (p = 0.04). For carmustine, a decrease of cells in G1-phase was detected pointing to cell cycle arrest in the ABCB1 knockout (p = 0.006). For lomustine, however, loss of ABCB1 did not alter the response to the treatment. Overall, this study shows that ABCB1 is involved in the active transport of temozolomide out of the tumor cells diminishing the response to temozolomide. Interestingly, loss of ABCB1 also affected the response to the lipophilic drug carmustine. These findings show that ABCB1 is not only relevant at the blood-brain barrier, but also in the tumor cells diminishing success of chemotherapy.

Correction: Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease.

[This corrects the article DOI: 10.1371/journal.pgen.1009679.].

Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease.

Numerous genetic studies have established a role for rare genomic variants in Congenital Heart Disease (CHD) at the copy number variation (CNV) and de novo variant (DNV) level. To identify novel haploinsufficient CHD disease genes, we performed an integrative analysis of CNVs and DNVs identified in probands with CHD including cases with sporadic thoracic aortic aneurysm. We assembled CNV data from 7,958 cases and 14,082 controls and performed a gene-wise analysis of the burden of rare genomic deletions in cases versus controls. In addition, we performed variation rate testing for DNVs identified in 2,489 parent-offspring trios. Our analysis revealed 21 genes which were significantly affected by rare CNVs and/or DNVs in probands. Fourteen of these genes have previously been associated with CHD while the remaining genes (FEZ1, MYO16, ARID1B, NALCN, WAC, KDM5B and WHSC1) have only been associated in small cases series or show new associations with CHD. In addition, a systems level analysis revealed affected protein-protein interaction networks involved in Notch signaling pathway, heart morphogenesis, DNA repair and cilia/centrosome function. Taken together, this approach highlights the importance of re-analyzing existing datasets to strengthen disease association and identify novel disease genes and pathways.

Genome sequencing in families with congenital limb malformations.

The extensive clinical and genetic heterogeneity of congenital limb malformation calls for comprehensive genome-wide analysis of genetic variation. Genome sequencing (GS) has the potential to identify all genetic variants. Here we aim to determine the diagnostic potential of GS as a comprehensive one-test-for-all strategy in a cohort of undiagnosed patients with congenital limb malformations. We collected 69 cases (64 trios, 1 duo, 5 singletons) with congenital limb malformations with no molecular diagnosis after standard clinical genetic testing and performed genome sequencing. We also developed a framework to identify potential noncoding pathogenic variants. We identified likely pathogenic/disease-associated variants in 12 cases (17.4%) including four in known disease genes, and one repeat expansion in HOXD13. In three unrelated cases with ectrodactyly, we identified likely pathogenic variants in UBA2, establishing it as a novel disease gene. In addition, we found two complex structural variants (3%). We also identified likely causative variants in three novel high confidence candidate genes. We were not able to identify any noncoding variants. GS is a powerful strategy to identify all types of genomic variants associated with congenital limb malformation, including repeat expansions and complex structural variants missed by standard diagnostic approaches. In this cohort, no causative noncoding SNVs could be identified.

Molecular characterization of Burkitt lymphoma in the breast or ovary.

Breast and ovary have been described as rare but typical sites of presentation of Burkitt lymphoma (BL) in females, particularly after puberty. We revised a historic series of 44 lymphomas of the breast or the ovary in women diagnosed between 1973 and 2014 as BL. Fluorescence in situ hybridization (FISH) was applied to all, and array-based copy number analysis as well as expression profiling to a subset of those cases. Of the 42 cases evaluable for FISH, 19 cases showed an IG-MYC translocation but only 9 of those fulfilled the criteria of the current WHO classification for the diagnosis of BL. Those nine cases resembled BL of other sites with regard to molecular features. Our findings along with literature data suggest that breast and ovarian BL (1) seem to be rarer than hitherto assumed, (2) share typical molecular features with other BL, and (3) predominantly affect women in the fertile age.

ZFP36L1 plays an ambiguous role in the regulation of cell expansion and negatively regulates CDKN1A in chronic myeloid leukemia cells.

The mRNA-destabilizing proteins ZFP36L1 and ZFP36L2 are described as mediators of quiescence and play a pivotal role in hematopoietic malignancies. Both genes are mainly classified as tumor suppressor genes as they posttranscriptionally downregulate the expression of oncogenes and contribute to cellular quiescence. Here, we analyzed the role of ZFP36L1 and ZFP36L2 in chronic myeloid leukemia (CML). We found ZFP36L1 and ZFP36L2 expression to be deregulated in patients with CML. By use of in vitro models of tyrosine kinase inhibitor resistance, an increase in ZFP36L1 and ZFP36L2 expression was detected during the development of imatinib resistance. CRISPR/Cas9-derived knockout of ZFP36L1, but not of ZFP36L2, in imatinib-sensitive cells led to decreased proliferation rates in response to tyrosine kinase inhibitor treatment. This effect was also observed in untreated ZFP36L1 knockout cells, albeit to a lower extent. Genomewide gene expression analyses of ZFP36L1 knockout cells revealed differential expression of cell cycle regulators, in particular upregulation of the cell cycle inhibitor CDKN1A. In addition, the 3' untranslated region of CDKN1A was proven to be a direct target of ZFP36L1. This indicates that tumor suppressor genes can also be targeted by ZFP36L1. Hence, ZFP36L1 cannot unambiguously be regarded as a tumor suppressor gene.

Cancer incidence and spectrum among children with genetically confirmed Beckwith-Wiedemann spectrum in Germany: a retrospective cohort study.

BACKGROUND: Beckwith-Wiedemann syndrome (BWS) is a cancer predisposition syndrome caused by defects on chromosome 11p15.5. The quantitative cancer risks in BWS patients depend on the underlying (epi)genotype but have not yet been assessed in a population-based manner. METHODS: We identified a group of 321 individuals with a molecularly confirmed diagnosis of BWS and analysed the cancer incidence up to age 15 years and cancer spectrum by matching their data with the German Childhood Cancer Registry. RESULTS: We observed 13 cases of cancer in the entire BWS cohort vs 0.4 expected. This corresponds to a 33-fold increased risk (standardised incidence ratio (SIR) = 32.6; 95% confidence interval = 17.3-55.7). The specific cancers included hepatoblastoma (n = 6); nephroblastoma (n = 4); astrocytoma (n = 1); neuroblastoma (n = 1) and adrenocortical carcinoma (n = 1). The cancer SIR was highest in patients with a paternal uniparental disomy of 11p15.5 (UPDpat). A high cancer risk remained when cases of cancer diagnosed prior to the BWS diagnosis were excluded. CONCLUSIONS: This study confirms an increased cancer risk in children with BWS. Our findings suggest that the highest cancer risk is associated with UPDpat. We were unable to confirm an excessive cancer risk in patients with IC1 gain of methylation (IC1-GOM) and this finding requires further investigation.

Decoding the DNA Methylome of Mantle Cell Lymphoma in the Light of the Entire B Cell Lineage.

We analyzed the in silico purified DNA methylation signatures of 82 mantle cell lymphomas (MCL) in comparison with cell subpopulations spanning the entire B cell lineage. We identified two MCL subgroups, respectively carrying epigenetic imprints of germinal-center-inexperienced and germinal-center-experienced B cells, and we found that DNA methylation profiles during lymphomagenesis are largely influenced by the methylation dynamics in normal B cells. An integrative epigenomic approach revealed 10,504 differentially methylated regions in regulatory elements marked by H3K27ac in MCL primary cases, including a distant enhancer showing de novo looping to the MCL oncogene SOX11. Finally, we observed that the magnitude of DNA methylation changes per case is highly variable and serves as an independent prognostic factor for MCL outcome.

Phenotypic spectrum and extent of DNA methylation defects associated with multilocus imprinting disturbances.

AIM: To characterize the genotypic and phenotypic extent of multilocus imprinting disturbances (MLID). MATERIALS & METHODS: We analyzed 37 patients with imprinting disorders (explorative cohort) for DNA methylation changes using the Infinium HumanMethylation450 BeadChip. For validation, three independent cohorts with imprinting disorders or cardinal features thereof were analyzed (84 patients with imprinting disorders, 52 with growth disorder, 81 with developmental delay). RESULTS: In the explorative cohort 21 individuals showed array-based MLID with each one displaying an Angelman or Temple syndrome phenotype, respectively. Epimutations in ZDBF2 and FAM50B were associated with severe MLID regarding number of affected regions. By targeted analysis we identified methylation changes of ZDBF2 and FAM50B also in the three validation cohorts. CONCLUSION: We corroborate epimutations in ZDBF2 and FAM50B as frequent changes in MLID whereas these rarely occur in other patients with cardinal features of imprinting disorders. Moreover, we show cell lineage specific differences in the genomic extent of FAM50B epimutation.

Identification of ZCCHC8 as fusion partner of ROS1 in a case of congenital glioblastoma multiforme with a t(6;12)(q21;q24.3).

Congenital gliobastoma multiforme (GBM) is rare and little is known about the molecular defects underlying the initiation and progression of this tumor type. We present a case of congenital GBM analyzed by conventional cytogenetics, fluorescence in situ hybridization, array comparative genomic hybridization and next generation sequencing. On cytogenetic analysis we detected a reciprocal translocation t(6;12)(q21;q24.3). By sequencing, the translocation was shown to form a fusion between the 5' region of ZCCHC8 and the 3' region of ROS1. RT-PCR analyses confirmed the existence of an in-frame fusion transcript with ZCCHC8 exons 1-3 joined to ROS1 exons 36-43. In addition to the ZCCHC8-ROS1 fusion, we detected a deletion in the short arm of chromosome 9, including homozygous loss of the CDKN2A/2B locus in 9p21.3 and heterozygous deletion of the HAUS6 gene in 9p22.1. The latter encodes a protein involved in faithful chromosome segregation by regulating microtubule nucleation and its deletion might be associated with the marked subclonal changes of ploidy observed in the tumor. This report adds the ZCCHC8-ROS1 fusion as oncogenic driver in GBM and supports the role of ROS1 activation in the pathogenesis of a subset of GBM. © 2016 Wiley Periodicals, Inc.

A summary of molecular genetic findings in fructose-1,6-bisphosphatase deficiency with a focus on a common long-range deletion and the role of MLPA analysis.

BACKGROUND: Fructose-1,6-bisphosphatase deficiency is a rare inborn error of metabolism affecting gluconeogenesis with only sporadic reports on its molecular genetic basis. RESULTS: We report our experience with mutation analysis in 14 patients (13 families) with fructose-1,6-bisphosphatase deficiency using conventional Sanger sequencing and multiplex ligation-dependent probe amplification analysis, and we provide a mutation update for the fructose bisphosphatase-1 gene (FBP1). Mutations were found on both chromosomes in all of our 14 patients including 5 novel mutations. Among the novel mutations is a 5412-bp deletion (c.-24-26_170 + 5192del) including the entire coding sequence of exon 2 of FBP1 that was repeatedly found in patients from Turkey and Armenia which may explain earlier poorly defined findings in patients from this area. This deletion can be detected with specific primers by generation of a junction fragment and by MLPA and SNP array assays. MLPA analysis was able to detect copy number variations in two further patients, one heterozygous for a deletion within exon 8, another heterozygous for a novel deletion of the entire FBP1 gene. CONCLUSIONS: Based on our update for the FBP1 gene, currently listing 35 mutations worldwide, and knowledge of PCR conditions that allow simple detection of a common FBP1 deletion in the Armenian and Turkish population, molecular genetic diagnosis has become easier in FBP1 deficiency. Furthermore, MLPA analysis may plays a useful role in patients with this disorder.

Genotype/phenotype correlation in a female patient with 21q22.3 and 12p13.33 duplications.

Many chromosomal rearrangements that lead to copy-number gains or losses have been shown to cause distinctive and recognizable clinical phenotypes. Conventional cytogenetic analysis can detect many, but not all, rearrangements depending on its power of resolution. The wide use of whole-genome array-based comparative genomic hybridization (array-CGH) techniques has allowed the detection of novel syndromes and to establish genotype-phenotype correlations by delineating at high resolution the regions involved in specific chromosomal aberrations. We report on a two and half-year-old female patient with intellectual disability and distinctive phenotypic features resulting from a de novo duplication of about 0.3 Mb in 21q22.3 associated with duplication of about 0.3 Mb in 12p13.33. The patient's chromosomal abnormalities were identified at the cytogenetic molecular level, using SNP array analysis, while GTG banding technique revealed a normal karyotype. Clinical findings of the patient were compared with Down syndrome and 12p duplication syndrome. This study suggests that an area of contiguous genes on the distal part of chromosome 21 (21q22.3) contribute to the Down syndrome phenotype and indicates that genes in the distal region of 12p (12p13.33) account for many facial characteristics and hypotonia of trisomy 12p syndrome.

DNA methylome analysis in Burkitt and follicular lymphomas identifies differentially methylated regions linked to somatic mutation and transcriptional control.

Although Burkitt lymphomas and follicular lymphomas both have features of germinal center B cells, they are biologically and clinically quite distinct. Here we performed whole-genome bisulfite, genome and transcriptome sequencing in 13 IG-MYC translocation-positive Burkitt lymphoma, nine BCL2 translocation-positive follicular lymphoma and four normal germinal center B cell samples. Comparison of Burkitt and follicular lymphoma samples showed differential methylation of intragenic regions that strongly correlated with expression of associated genes, for example, genes active in germinal center dark-zone and light-zone B cells. Integrative pathway analyses of regions differentially methylated in Burkitt and follicular lymphomas implicated DNA methylation as cooperating with somatic mutation of sphingosine phosphate signaling, as well as the TCF3-ID3 and SWI/SNF complexes, in a large fraction of Burkitt lymphomas. Taken together, our results demonstrate a tight connection between somatic mutation, DNA methylation and transcriptional control in key B cell pathways deregulated differentially in Burkitt lymphoma and other germinal center B cell lymphomas.

New insights into the genetics of X-linked dystonia-parkinsonism (XDP, DYT3).

X-linked recessive dystonia-parkinsonism is a rare movement disorder that is highly prevalent in Panay Island in the Philippines. Earlier studies identified seven different genetic alterations within a 427-kb disease locus on the X chromosome; however, the exact disease-causing variant among these is still not unequivocally determined. To further investigate the genetic cause of this disease, we sequenced all previously reported genetic alterations in 166 patients and 473 Filipino controls. Singly occurring variants in our ethnically matched controls would have allowed us to define these as polymorphisms, but none were found. Instead, we identified five patients carrying none of the disease-associated variants, and one male control carrying all of them. In parallel, we searched for novel single-nucleotide variants using next-generation sequencing. We did not identify any shared variants in coding regions of the X chromosome. However, by validating intergenic variants discovered via genome sequencing, we were able to define the boundaries of the disease-specific haplotype and narrow the disease locus to a 294-kb region that includes four known genes. Using microarray-based analyses, we ruled out the presence of disease-linked copy number variants within the implicated region. Finally, we utilized in silico analysis and detected no strong evidence of regulatory regions surrounding the disease-associated variants. In conclusion, our finding of disease-specific variants occurring in complete linkage disequilibrium raises new insights and intriguing questions about the origin of the disease haplotype, the existence of phenocopies and of reduced penetrance, and the causative genetic alteration in XDP.

A familial disorder of altered DNA-methylation.

BACKGROUND: In a subset of imprinting disorders caused by epimutations, multiple imprinted loci are affected. Familial occurrence of multilocus imprinting disorders is rare. PURPOSE/OBJECTIVE: We have investigated the clinical and molecular features of a familial DNA-methylation disorder. METHODS: Tissues of affected individuals and blood samples of family members were investigated by conventional and molecular karyotyping. Sanger sequencing and RT-PCR of imprinting-associated genes (NLRP2, NLRP7, ZFP57, KHDC3L, DNMT1o), exome sequencing and locus-specific, array-based and genome-wide technologies to determine DNA-methylation were performed. RESULTS: In three offspring of a healthy couple, we observed prenatal onset of severe growth retardation and dysmorphism associated with altered DNA-methylation at paternally and maternally imprinted loci. Array-based analyses in various tissues of the offspring identified the DNA-methylation of 2.1% of the genes in the genome to be recurrently altered. Despite significant enrichment of imprinted genes (OR 9.49), altered DNA-methylation predominately (90.2%) affected genes not known to be imprinted. Sequencing of genes known to cause comparable conditions and exome sequencing in affected individuals and their ancestors did not unambiguously point to a causative gene. CONCLUSIONS: The family presented herein suggests the existence of a familial disorder of DNA-methylation affecting imprinted but also not imprinted gene loci potentially caused by a maternal effect mutation in a hitherto not identified gene.

A recurrent 11q aberration pattern characterizes a subset of MYC-negative high-grade B-cell lymphomas resembling Burkitt lymphoma.

The genetic hallmark of Burkitt lymphoma (BL) is the t(8;14)(q24;q32) and its variants leading to activation of the MYC oncogene. It is a matter of debate whether true BL without MYC translocation exists. Here, we identified 59 lymphomas concordantly called BL by 2 gene expression classifiers among 753 B-cell lymphomas. Only 2 (3%) of these 59 molecular BL lacked a MYC translocation, which both shared a peculiar pattern of chromosome 11q aberration characterized by interstitial gains including 11q23.2-q23.3 and telomeric losses of 11q24.1-qter. We extended our analysis to 17 MYC-negative high-grade B-cell lymphomas with a similar 11q aberration and showed this aberration to be recurrently associated with morphologic and clinical features of BL. The minimal region of gain was defined by high-level amplifications in 11q23.3 and associated with overexpression of genes including PAFAH1B2 on a transcriptional and protein level. The recurrent region of loss contained a focal homozygous deletion in 11q24.2-q24.3 including the ETS1 gene, which was shown to be mutated in 4 of 16 investigated cases. These findings indicate the existence of a molecularly distinct subset of B-cell lymphomas reminiscent of BL, which is characterized by deregulation of genes in 11q.

Recurrent mutation of JAK3 in T-cell prolymphocytic leukemia.

T-cell prolymphocytic leukemia (T-PLL) is an aggressive post-thymic T-cell malignancy characterized by the recurrent inv(14)(q11q32)/t(14;14)(q11;q32) or t(X;14)(q28;q11) leading to activation of either the TCL1 or MTCP1 gene, respectively. However, these primary genetic events are insufficient to drive leukemogenesis. Recently, activating mutations in JAK3 have been identified in other T-cell malignancies. Since JAK3 is essential for T-cell maturation, we analyzed a cohort of 32 T-PLL patients for mutational hot spots in the JAK3 gene using a step-wise screening approach. We identified 14 mutations in 11 of 32 patients (34%). The most frequently detected mutation in our cohort was M511I (seen in 57% of cases) previously described as an activating change in other T-cell malignancies. Three patients carried two mutations in JAK3. In two patients M511I and R657Q were simultaneously detected and in another patient V674F and V678L. In the latter case we could demonstrate that the mutations were on the same allele in cis. Protein modeling and homology analyses of mutations present in other members of the JAK family suggested that these mutations likely activate JAK3, possibly by disrupting the activation loop and the interface between N and C lobes, increasing the accessibility of the catalytic loop. In addition, four of the 21 patients lacking a JAK3 point mutation presented an aberrant karyotype involving the chromosomal band 19p13 harboring the JAK3 locus. The finding of recurrent activating JAK3 mutations in patients with T-PLL could enable the use of JAK3 inhibitors to treat patients with this unfavorable malignancy who otherwise have a very poor prognosis.

Genetic characteristics of the human hepatic stellate cell line LX-2.

The human hepatic cell line LX-2 has been described as tool to study mechanisms of hepatic fibrogenesis and the testing of antifibrotic compounds. It was originally generated by immortalisation with the Simian Vacuolating Virus 40 (SV40) transforming (T) antigen and subsequent propagation in low serum conditions. Although this immortalized line is used in an increasing number of studies, detailed genetic characterisation has been lacking. We here have performed genetic characterisation of the LX-2 cell line and established a single-locus short tandem repeat (STR) profile for the cell line and characterized the LX-2 karyotype by several cytogenetic and molecular cytogenetic techniques. Spectral karyotyping (SKY) revealed a complex karyotype with a set of aberrations consistently present in the metaphases analyses which might serve as cytogenetic markers. In addition, various subclonal and single cell aberrations were detected. Our study provides criteria for genetic authentication of LX-2 and offers insights into the genotype changes which might underlie part of its phenotypic features.