CEQer: a graphical tool for copy number and allelic imbalance detection from whole-exome sequencing data

Recurrent somatic mutations of FAT family cadherins induce an aggressive phenotype and poor prognosis in anaplastic large cell lymphoma

BACKGROUND

Anaplastic Large Cell Lymphoma (ALCL) is a rare and aggressive T-cell lymphoma, classified into ALK-positive and ALK-negative subtypes, based on the presence of chromosomal translocations involving the ALK gene. The current standard of treatment for ALCL is polychemotherapy, with a high overall survival rate. However, a subset of patients does not respond to or develops resistance to these therapies, posing a serious challenge for clinicians. Recent targeted treatments such as ALK kinase inhibitors and anti-CD30 antibody-drug conjugates have shown promise but, for a fraction of patients, the prognosis is still unsatisfactory.

METHODS

We investigated the genetic landscape of ALK + ALCL by whole-exome sequencing; recurring mutations were characterized in vitro and in vivo using transduced ALCL cellular models.

RESULTS

Recurrent mutations in FAT family genes and the transcription factor RUNX1T1 were found. These mutations induced changes in ALCL cells morphology, growth, and migration, shedding light on potential factors contributing to treatment resistance. In particular, FAT4 silencing in ALCL cells activated the β-catenin and YAP1 pathways, which play crucial roles in tumor growth, and conferred resistance to chemotherapy. Furthermore, STAT1 and STAT3 were hyper-activated in these cells. Gene expression profiling showed global changes in pathways related to cell adhesion, cytoskeletal organization, and oncogenic signaling. Notably, FAT mutations associated with poor outcome in patients.

CONCLUSIONS

These findings provide novel insights into the molecular portrait of ALCL, that could help improve treatment strategies and the prognosis for ALCL patients.

Late relapse of chronic myeloid leukemia after allogeneic bone marrow transplantation points to KANSARL (KANSL1::ARL17A) alteration: a case report with insights on the molecular landscape

Chronic myeloid leukemia is a myeloproliferative neoplasm characterized by the presence of the Philadelphia chromosome and the consequent BCR::ABL1 oncoprotein. In the era before the introduction of tyrosine kinase inhibitors (TKIs), the only potentially curative treatment was allogeneic hematopoietic stem cell transplantation (HSCT). Here, we present the case of a patient affected by CML who experienced a relapse 20 years after allogeneic HSCT. Following relapse, the patient was treated with imatinib and bosutinib, resulting in a deep molecular response and successfully discontinued treatment. Additional analysis including whole-exome sequencing and RNA sequencing provided some insights on the molecular mechanisms of the relapse: the identification of the fusion transcript KANSL1::ARL17A (KANSARL), a cancer predisposition fusion gene, could justify a condition of genomic instability which may be associated with the onset and/or probably the late relapse of his CML.

Idiopathic erythrocytosis: a germline disease?

Polycythemia Vera (PV) is typically caused by V617F or exon 12 JAK2 mutations. Little is known about Polycythemia cases where no JAK2 variants can be detected, and no other causes identified. This condition is defined as idiopathic erythrocytosis (IE). We evaluated clinical-laboratory parameters of a cohort of 56 IE patients and we determined their molecular profile at diagnosis with paired blood/buccal-DNA exome-sequencing coupled with a high-depth targeted OncoPanel to identify a possible underling germline or somatic cause. We demonstrated that most of our cohort (40/56: 71.4%) showed no evidence of clonal hematopoiesis, suggesting that IE is, in large part, a germline disorder. We identified 20 low mutation burden somatic variants (Variant allelic fraction, VAF, < 10%) in only 14 (25%) patients, principally involving DNMT3A and TET2. Only 2 patients presented high mutation burden somatic variants, involving DNMT3A, TET2, ASXL1 and WT1. We identified recurrent germline variants in 42 (75%) patients occurring mainly in JAK/STAT, Hypoxia and Iron metabolism pathways, among them: JAK3-V722I and HIF1A-P582S; a high fraction of patients (48.2%) resulted also mutated in homeostatic iron regulatory gene HFE-H63D or C282Y. By generating cellular models, we showed that JAK3-V722I causes activation of the JAK-STAT5 axis and upregulation of EPAS1/HIF2A, while HIF1A-P582S causes suppression of hepcidin mRNA synthesis, suggesting a major role for these variants in the onset of IE.

Detection of DNA copy number alterations by matrix-assisted laser desorption/ionization time-of-flight mass spectrometric analysis of single nucleotide polymorphisms

OBJECTIVES

Copy number alterations (CNAs) are frequently found in malignant tissues. Different approaches have been used for CNA detection. However, it is not easy to detect a large panel of CNA targets in heterogenous tumors.

METHODS

We have developed a CNAs detection approach through quantitatively analyzed allelic imbalance by allelotyping single nucleotide polymorphisms (SNPs) by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS). Furthermore, the copy number changes were quantified by real-competitive PCR (rcPCR) to distinguish loss of heterozygosity (LOH) and genomic amplification. The approach was used to validate the CNA regions detected by next generation sequencing (NGS) in early-stage lung carcinoma.

RESULTS

CNAs were detected in heterogeneous DNA samples where tumor DNA is present at only 10% through the SNP based allelotyping. In addition, two different types of CNAs (loss of heterozygosity and chromosome amplification) were able to be distinguished quantitatively by rcPCR. Validation on a total of 41 SNPs from the selected CNA regions showed that copy number changes did occur, and the tissues from early-stage lung carcinoma were distinguished from normal.

CONCLUSIONS

CNA detection by MALDI-TOF MS can be used for validating potentially interesting genomic regions identified from next generation sequencing, and for detecting CNAs in tumor tissues consisting of a mixture of neoplastic and normal cells.

ETNK1 mutations induce a mutator phenotype that can be reverted with phosphoethanolamine

Recurrent somatic mutations in ETNK1 (Ethanolamine-Kinase-1) were identified in several myeloid malignancies and are responsible for a reduced enzymatic activity. Here, we demonstrate in primary leukemic cells and in cell lines that mutated ETNK1 causes a significant increase in mitochondrial activity, ROS production, and Histone H2AX phosphorylation, ultimately driving the increased accumulation of new mutations. We also show that phosphoethanolamine, the metabolic product of ETNK1, negatively controls mitochondrial activity through a direct competition with succinate at mitochondrial complex II. Hence, reduced intracellular phosphoethanolamine causes mitochondria hyperactivation, ROS production, and DNA damage. Treatment with phosphoethanolamine is able to counteract complex II hyperactivation and to restore a normal phenotype.

Integrated Genomic, Functional, and Prognostic Characterization of Atypical Chronic Myeloid Leukemia

Supplemental Digital Content is available for this article.

Atypical chronic myeloid leukemia (aCML) is a BCR-ABL1-negative clonal disorder, which belongs to the myelodysplastic/myeloproliferative group. This disease is characterized by recurrent somatic mutations in SETBP1, ASXL1 and ETNK1 genes, as well as high genetic heterogeneity, thus posing a great therapeutic challenge. To provide a comprehensive genomic characterization of aCML we applied a high-throughput sequencing strategy to 43 aCML samples, including both whole-exome and RNA-sequencing data. Our dataset identifies ASXL1, SETBP1, and ETNK1 as the most frequently mutated genes with a total of 43.2%, 29.7 and 16.2%, respectively. We characterized the clonal architecture of 7 aCML patients by means of colony assays and targeted resequencing. The results indicate that ETNK1 variants occur early in the clonal evolution history of aCML, while SETBP1 mutations often represent a late event. The presence of actionable mutations conferred both ex vivo and in vivo sensitivity to specific inhibitors with evidence of strong in vitro synergism in case of multiple targeting. In one patient, a clinical response was obtained. Stratification based on RNA-sequencing identified two different populations in terms of overall survival, and differential gene expression analysis identified 38 significantly overexpressed genes in the worse outcome group. Three genes correctly classified patients for overall survival.

Identification of genetic polymorphisms modulating nausea and vomiting in two series of opioid-treated cancer patients

Nausea and vomiting are often associated with opioid analgesia in cancer patients; however, only a subset of patients develop such side effects. Here, we tested the hypothesis that the occurrence of nausea and vomiting is modulated by the genetic background of the patients. Whole exome sequencing of DNA pools from patients with either low (n = 937) or high (n = 557) nausea and vomiting intensity, recruited in the European Pharmacogenetic Opioid Study, revealed a preliminary association of 53 polymorphisms. PCR-based genotyping of 45 of these polymorphisms in the individual patients of the same series confirmed the association for six SNPs in AIM1L, CLCC1, MUC16, PDE3A, POM121L2, and ZNF165 genes. Genotyping of the same 45 polymorphisms in 264 patients of the Italian CERP study, also treated with opioids for cancer pain, instead confirmed the association for two SNPs in ZNF568 and PDE3A genes. Only one SNP, rs12305038 in PDE3A, was confirmed in both series, although with opposite effects of the minor allele on the investigated phenotype. Overall, our findings suggest that genetic factors are indeed associated with nausea and vomiting in opioid-treated cancer patients, but the role of individual polymorphisms may be weak.

Deep sequencing of a recurrent oligodendroglioma and the derived xenografts reveals new insights into the evolution of human oligodendroglioma and candidate driver genes

We previously reported the establishment of a rare xenograft derived from a recurrent oligodendroglioma with 1p/19q codeletion. Here, we analyzed in detail the exome sequencing datasets from the recurrent oligodendroglioma (WHO grade III, recurrent O²⁰¹⁰) and the first-generation xenograft (xenograft¹). Somatic SNVs and small InDels (n = 80) with potential effects at the protein level in recurrent O²⁰¹⁰ included variants in IDH1 (NM_005896:c.395G>A; p. Arg132His), FUBP1 (NM_003902:c.1307_1310delTAGA; p.Ile436fs), and CIC (NM_015125:c.4421T>G; p.Val1474Gly). All but 2 of these 80 variants were also present in xenograft¹, along with 7 new variants. Deep sequencing of the 87 SNVs and InDels in the original tumor (WHO grade III, primary O²⁰⁰⁵) and in a second-generation xenograft (xenograft²) revealed that only 11 variants, including IDH1 (NM_005896:c.395G>A; p. Arg132His), PSKH1 (NM_006742.2:c.650G>A; p.Arg217Gln), and SNX12 (NM_001256188:c.470G>A; p.Arg157His), along with a variant in the TERT promoter (C250T, NM_198253.2: c.-146G>A), were already present in primary O²⁰⁰⁵. Allele frequencies of the 11 variants were calculated to assess their potential as putative driver genes. A missense change in NDST4 (NM_022569:c.2392C>G; p.Leu798Val) on 4q exhibited an increasing allele frequency (~ 20%, primary O²⁰⁰⁵, 80%, recurrent O²⁰¹⁰ and 100%, xenograft¹), consistent with a selection event. Sequencing of NDST4 in a cohort of 15 oligodendrogliomas, however, revealed no additional cases with potential protein disrupting variants. Our analysis illuminated a tumor evolutionary series of events, which included 1p/19q codeletion, IDH1 R132H, and TERT C250T as early events, followed by loss of function of NDST4 and mutations in FUBP1 and CIC as late events.

Lorlatinib Treatment Elicits Multiple On- and Off-Target Mechanisms of Resistance in ALK-Driven Cancer

: Targeted therapy changed the standard of care in ALK-dependent tumors. However, resistance remains a major challenge. Lorlatinib is a third-generation ALK inhibitor that inhibits most ALK mutants resistant to current ALK inhibitors. In this study, we utilize lorlatinib-resistant anaplastic large cell lymphoma (ALCL), non-small cell lung cancer (NSCLC), and neuroblastoma cell lines in vitro and in vivo to investigate the acquisition of resistance and its underlying mechanisms. ALCL cells acquired compound ALK mutations G1202R/G1269A and C1156F/L1198F in vitro at high drug concentrations. ALCL xenografts selected in vivo showed recurrent N1178H (5/10 mice) and G1269A (4/10 mice) mutations. Interestingly, intracellular localization of NPM/ALKN^1178H skewed toward the cytoplasm in human cells, possibly mimicking overexpression. RNA sequencing of resistant cells showed significant alteration of PI3K/AKT and RAS/MAPK pathways. Functional validation by small-molecule inhibitors confirmed the involvement of these pathways in resistance to lorlatinib. NSCLC cells exposed in vitro to lorlatinib acquired hyperactivation of EGFR, which was blocked by erlotinib to restore sensitivity to lorlatinib. In neuroblastoma, whole-exome sequencing and proteomic profiling of lorlatinib-resistant cells revealed a truncating NF1 mutation and hyperactivation of EGFR and ErbB4. These data provide an extensive characterization of resistance mechanisms that may arise in different ALK-positive cancers following lorlatinib treatment. SIGNIFICANCE: High-throughput genomic, transcriptomic, and proteomic profiling reveals various mechanisms by which multiple tumor types acquire resistance to the third-generation ALK inhibitor lorlatinib.

OncoScore: a novel, Internet-based tool to assess the oncogenic potential of genes

The complicated, evolving landscape of cancer mutations poses a formidable challenge to identify cancer genes among the large lists of mutations typically generated in NGS experiments. The ability to prioritize these variants is therefore of paramount importance. To address this issue we developed OncoScore, a text-mining tool that ranks genes according to their association with cancer, based on available biomedical literature. Receiver operating characteristic curve and the area under the curve (AUC) metrics on manually curated datasets confirmed the excellent discriminating capability of OncoScore (OncoScore cut-off threshold = 21.09; AUC = 90.3%, 95% CI: 88.1-92.5%), indicating that OncoScore provides useful results in cases where an efficient prioritization of cancer-associated genes is needed.

Identification of genomic deletions causing inherited retinal degenerations by coverage analysis of whole exome sequencing data

BACKGROUND

Inherited retinal degenerations (IRDs) are a common cause of visual disturbance with a high clinical and genetic heterogeneity. Recent sequencing techniques such as whole exome sequencing (WES) contribute to the discovery of novel genes. The aim of the current study was to use WES data to identify large deletions that include at least one exon in known IRD genes.

METHODS

Patients diagnosed with IRDs underwent a comprehensive ophthalmic evaluation. WES was performed using the NimbleGen V2 paired-end kit and HiSeq 2000. An analysis of exon coverage data was performed on 60 WES samples. Exonic deletions were verified by 'PCR walking' analysis.

RESULTS

We analysed data obtained from 60 WES samples of index patients with IRDs. By calculating the average coverage for all exons in the human genome, we were able to identify homozygous and hemizygous deletions of at least one exon in six families (10%), including a single-exon deletion in EYS, deletions of three consecutive exons in MYO7A and NPHP4, deletions of four and eight consecutive exons in RPGR and a multigene deletion on the X-chromosome, including CHM. By using PCR-walking analysis, we were able to identify the borders of five of the deletions and to screen our set of patients for these deletions.

CONCLUSIONS

We performed here a comprehensive analysis of WES data as a tool for identifying large genomic deletions in patients with IRDs. Our analysis indicates that large deletions are relatively frequent (about 10% of our WES cohort) and should be screened when analysing WES data.

Intragenic duplication in the PKHD1 gene in autosomal recessive polycystic kidney disease

Background

In the present study, we report on a couple who underwent prenatal genetic diagnosis for autosomal recessive polycystic kidney disease (ARPKD).

Case presentation

This healthy couple had previously had a healthy boy but had experienced two consecutive neonatal deaths due to respiratory distress resulting from pulmonary hypoplasia caused by oligohydramnios. The woman consulted our facility after she realized she was pregnant again. We promptly performed a carrier test for the PKHD1 gene by target exome sequencing of samples from the couple. A pathogenic mutation was identified only in the paternal allele (c.9008C>T, p.S3003F). The mutation was confirmed by Sanger sequencing of the DNA from formalin-fixed, paraffin-embedded, kidney tissue of the second neonate patient and was not found in the healthy sibling. We then performed haplotype analyses using microsatellite markers scattered throughout the PKHD1 gene. DNA from the amniocentesis was determined to belong to a carrier, and the couple decided to continue with the pregnancy, obtaining a healthy newborn. Subsequent detailed examination of the exome data suggested higher read depth at exons 45 and 46. Multiplex ligation-dependent probe amplification allowed identification of duplication of these two exons. This case suggests the potential usefulness of target exome sequencing in the prenatal diagnosis of the PKHD1 gene in ARPKD.

Conclusions

This is the first report of intragenic duplication in the PKHD1 gene in ARPKD.

Electronic supplementary material

The online version of this article (doi:10.1186/s12881-015-0245-3) contains supplementary material, which is available to authorized users.

Preimplantation genetic diagnosis/screening by comprehensive molecular testing

Although embryo screening by preimplantation genetic diagnosis (PGD) has become the standard technique for the treatment of recurrent pregnancy loss in couples with a balanced gross chromosomal rearrangement, the implantation and pregnancy rates of PGD using conventional fluorescence in situ hybridization (FISH) remain suboptimal. Comprehensive molecular testing, such as array comparative genomic hybridization and next-generation sequencing, can improve these rates, but amplification bias in the whole genome amplification method remains an obstacle to accurate diagnosis. Recent advances in amplification procedures combined with improvements in the microarray platform and analytical method have overcome the amplification bias, and the data accuracy of the comprehensive PGD method has reached the level of clinical laboratory testing. Currently, comprehensive PGD is also applied to recurrent pregnancy loss due to recurrent fetal aneuploidy or infertility with recurrent implantation failure, known as preimplantation genetic screening. However, there are still numerous problems to be solved, including misdiagnosis due to somatic mosaicism, cell cycle-related background noise, and difficulty in diagnosis of polyploidy. The technology for comprehensive PGD also requires further improvement.

Treatment Efficacy and Resistance Mechanisms Using the Second-Generation ALK Inhibitor AP26113 in Human NPM-ALK-Positive Anaplastic Large Cell Lymphoma

ALK is a tyrosine kinase receptor involved in a broad range of solid and hematologic tumors. Among 70% to 80% of ALK(+) anaplastic large cell lymphomas (ALCL) are caused by the aberrant oncogenic fusion protein NPM-ALK. Crizotinib was the first clinically relevant ALK inhibitor, now approved for the treatment of late-stage and metastatic cases of lung cancer. However, patients frequently develop drug resistance to Crizotinib, mainly due to the appearance of point mutations located in the ALK kinase domain. Fortunately, other inhibitors are available and in clinical trial, suggesting the potential for second-line therapies to overcome Crizotinib resistance. This study focuses on the ongoing phase I/II trial small-molecule tyrosine kinase inhibitor (TKI) AP26113, by Ariad Pharmaceuticals, which targets both ALK and EGFR. Two NPM-ALK(+) human cell lines, KARPAS-299 and SUP-M2, were grown in the presence of increasing concentrations of AP26113, and eight lines were selected that demonstrated resistance. All lines show IC50 values higher (130 to 1,000-fold) than the parental line. Mechanistically, KARPAS-299 populations resistant to AP26113 show NPM-ALK overexpression, whereas SUP-M2-resistant cells harbor several point mutations spanning the entire ALK kinase domain. In particular, amino acid substitutions: L1196M, S1206C, the double F1174V+L1198F and L1122V+L1196M mutations were identified. The knowledge of the possible appearance of new clinically relevant mechanisms of drug resistance is a useful tool for the management of new TKI-resistant cases.

IMPLICATIONS

This work defines reliable ALCL model systems of AP26113 resistance and provides a valuable tool in the management of all cases of relapse upon NPM-ALK-targeted therapy.

YMAP: a pipeline for visualization of copy number variation and loss of heterozygosity in eukaryotic pathogens

The design of effective antimicrobial therapies for serious eukaryotic pathogens requires a clear understanding of their highly variable genomes. To facilitate analysis of copy number variations, single nucleotide polymorphisms and loss of heterozygosity events in these pathogens, we developed a pipeline for analyzing diverse genome-scale datasets from microarray, deep sequencing, and restriction site associated DNA sequence experiments for clinical and laboratory strains of Candida albicans, the most prevalent human fungal pathogen. The Y_MAP pipeline (http://lovelace.cs.umn.edu/Ymap/) automatically illustrates genome-wide information in a single intuitive figure and is readily modified for the analysis of other pathogens with small genomes.

Recurrent ETNK1 mutations in atypical chronic myeloid leukemia

Despite the recent identification of recurrent SETBP1 mutations in atypical chronic myeloid leukemia (aCML), a complete description of the somatic lesions responsible for the onset of this disorder is still lacking. To find additional somatic abnormalities in aCML, we performed whole-exome sequencing on 15 aCML cases. In 2 cases (13.3%), we identified somatic missense mutations in the ETNK1 gene. Targeted resequencing on 515 hematological clonal disorders revealed the presence of ETNK1 variants in 6 (8.8%) of 68 aCML and 2 (2.6%) of 77 chronic myelomonocytic leukemia samples. These mutations clustered in a small region of the kinase domain, encoding for H243Y and N244S (1/8 H243Y; 7/8 N244S). They were all heterozygous and present in the dominant clone. The intracellular phosphoethanolamine/phosphocholine ratio was, on average, 5.2-fold lower in ETNK1-mutated samples (P < .05). Similar results were obtained using myeloid TF1 cells transduced with ETNK1 wild type, ETNK1-N244S, and ETNK1-H243Y, where the intracellular phosphoethanolamine/phosphocholine ratio was significantly lower in ETNK1-N244S (0.76 ± 0.07) and ETNK1-H243Y (0.37 ± 0.02) than in ETNK1-WT (1.37 ± 0.32; P = .01 and P = .0008, respectively), suggesting that ETNK1 mutations may inhibit the catalytic activity of the enzyme. In summary, our study shows for the first time the evidence of recurrent somatic ETNK1 mutations in the context of myeloproliferative/myelodysplastic disorders.

Exome sequencing identifies frequent mutation of MLL2 in non-small cell lung carcinoma from Chinese patients

Lung cancer is the most common cause of cancer mortality worldwide, with an estimated 1.4 million deaths each year. Here we report whole-exome sequencing of nine tumor/normal tissue pairs from Chinese patients with non-small cell lung carcinoma (NSCLC). This allows us to identify a number of significantly mutated genes in NSCLC, which were highly enriched in DNA damage repair, NF-κB pathway, JAK/STAT signaling and chromatin modification. Notably, we identify a histone-lysine methyltransferase gene, namely, MLL2, as one of the most significantly mutated genes in our screen. In a following validation study, we identify deleterious mutations of MLL2 in 12 out of 105 (11.4%) NSCLC patients. Additionally, reduced or lost expression of MLL2 was commonly observed in tumor tissues as compared with paired adjacent non-tumor tissues regardless of mutation status. Together, our study defines the landscape of somatic mutations in Chinese NSCLC and supports the role of MLL2 mutation in the pathogenesis of the disease.