Identification of trunk mutations in gastric carcinoma: a case study

Background

Intratumor heterogeneity (ITH) poses an urgent challenge for cancer precision medicine because it can cause drug resistance against cancer target therapy and immunotherapy. The search for trunk mutations that are present in all cancer cells is therefore critical for each patient.

Case presentation

In this study, we aimed to evaluate the efficiency of multiregional sequencing for the identification of trunk mutations present in all regions of a tumor as a case study. We applied multiregional whole-exome sequencing (WES) to investigate the genetic heterogeneity and homogeneity of a case of gastric carcinoma. Approximately 83% of common missense mutations present in two samples and approximately 89% of common missense mutations present in three samples were trunk mutations. Notably, trunk mutations appeared to have higher variant allele frequencies (VAFs) than non-trunk mutations.

Conclusions

Our results indicate that small-scale multiregional sampling and subsequent screening of low VAF somatic mutations might be a cost-effective strategy for identifying the majority of trunk mutations in gastric carcinoma.

     

Genetic evolution of nevus of Ota reveals clonal heterogeneity acquiring BAP1 and TP53 mutations

Melanoma presents molecular alterations based on its anatomical location and exposure to environmental factors. Due to its intrinsic genetic heterogeneity, a simple snapshot of a tumor's genetic alterations does not reflect the tumor clonal complexity or specific gene–gene cooperation. Here, we studied the genetic alterations and clonal evolution of a unique patient with a Nevus of Ota that developed into a recurring uveal‐like dermal melanoma. The Nevus of Ota and ulterior lesions contained GNAQ mutations were c‐KIT positive, and tumors showed an increased RAS pathway activity during progression. Whole‐exome sequencing of these lesions revealed the acquisition of BAP1 and TP53 mutations during tumor evolution, thereby unmasking clonal heterogeneity and allowing the identification of cooperating genes within the same tumor. Our results highlight the importance of studying tumor genetic evolution to identify cooperating mechanisms and delineate effective therapies.     

X-Linked chronic granulomatous disease: mutations in the CYBB gene encoding the gp91-phox component of respiratory-burst oxidase.

Chronic granulomatous disease (CGD) is a hereditary disorder of host defense due to absent or decreased activity of phagocyte NADPH oxidase. The X-linked form of the disease derives from defects in the CYBB gene, which encodes the 91-kD glycoprotein component (termed "gp91-phox") of the oxidase. We have identified the mutations in the CYBB gene responsible for X-linked CGD in 131 consecutive independent kindreds. Screening by SSCP analysis identified mutations in 124 of the kindreds, and sequencing of all exons and intron boundary regions revealed the other seven mutations. We detected 103 different specific mutations; no single mutation appeared in more than seven independent kindreds. The types of mutations included large and small deletions (11%), frameshifts (24%), nonsense mutations (23%), missense mutations (23%), splice-region mutations (17%), and regulatory-region mutations (2%). The distribution of mutations within the CYBB gene exhibited great heterogeneity, with no apparent mutational hot spots. Evaluation of 87 available mothers revealed X-linked carrier status in all but 10. The heterogeneity of mutations and the lack of any predominant genotype indicate that the disease represents many different mutational events, without a founder effect, as is expected for a disorder with a previously lethal phenotype.     

Evolutionary dynamics of intratumor heterogeneity

Intraneoplastic diversity in human tumors is a widespread phenomenon of critical importance for tumor progression and the response to therapeutic intervention. Insights into the evolutionary events that control tumor heterogeneity would be a major breakthrough in our comprehension of cancer development and could lead to more effective prevention methods and therapies. In this paper, we design an evolutionary mathematical framework to study the dynamics of heterogeneity over time. We consider specific situations arising during tumorigenesis, such as the emergence of positively selected mutations ("drivers") and the accumulation of neutral variation ("passengers"). We perform exact computer simulations of the emergence of diverse tumor cell clones over time, and derive analytical estimates for the extent of heterogeneity within a population of cancer cells. Our methods contribute to a quantitative understanding of tumor heterogeneity and the impact of heritable alterations on this tumor trait.     

SciClone: Inferring Clonal Architecture and Tracking the Spatial and Temporal Patterns of Tumor Evolution

The sensitivity of massively-parallel sequencing has confirmed that most cancers are oligoclonal, with subpopulations of neoplastic cells harboring distinct mutations. A fine resolution view of this clonal architecture provides insight into tumor heterogeneity, evolution, and treatment response, all of which may have clinical implications. Single tumor analysis already contributes to understanding these phenomena. However, cryptic subclones are frequently revealed by additional patient samples (e.g., collected at relapse or following treatment), indicating that accurately characterizing a tumor requires analyzing multiple samples from the same patient. To address this need, we present SciClone, a computational method that identifies the number and genetic composition of subclones by analyzing the variant allele frequencies of somatic mutations. We use it to detect subclones in acute myeloid leukemia and breast cancer samples that, though present at disease onset, are not evident from a single primary tumor sample. By doing so, we can track tumor evolution and identify the spatial origins of cells resisting therapy.     

Detection of mosaic pattern of mitochondrial DNA alterations in different populations of cells from the same endometrial tumor

Somatic mitochondrial DNA (mtDNA) alterations including point mutations and microsatellite instability (MSI) have been frequently detected in human cancers. To further explore the extensiveness of mtDNA alterations, we have analyzed the occurrence of somatic mtDNA mutations in different populations of endometrial cancer cells from the same tumor tissues as compared with adjacent non-tumor cells. Laser-captured micro-dissection was used to harvest endometrial cancer cells from separated areas of the same tumor and adjacent normal cells. Total DNA isolated from micro-dissected cells was PCR amplified and analyzed for mtDNA alterations by polyacrylamide gel electrophoresis and DNA sequencing. Multiple mtDNA alterations were detected in different portions of the same tumor. Different populations of endometrial cancer cells carried different patterns of mtDNA mutations. Interestingly, unlike previous reports, most mutations were found to be heteroplasmic. We have demonstrated the occurrence of hyper-variability of mtDNA alterations in a single piece of tumor tissue. Our observations support the hypothesis that the accumulation of mtDNA alterations is random and expands independently. The data presented here showed the heterogeneity of cancer cells in terms of mtDNA alterations in endometrial cancer.     

eTumorMetastasis: A Network-based Algorithm Predicts Clinical Outcomes Using Whole-exome Sequencing Data of Cancer Patients

Continual reduction in sequencing cost is expanding the accessibility of genome sequencing data for routine clinical applications. However, the lack of methods to construct machine learning-based predictive models using these datasets has become a crucial bottleneck for the application of sequencing technology in clinics. Here, we develop a new algorithm, eTumorMetastasis, which transforms tumor functional mutations into network-based profiles and identifies network operational gene (NOG) signatures. NOG signatures model the tipping point at which a tumor cell shifts from a state that doesn’t favor recurrence to one that does. We show that NOG signatures derived from genomic mutations of tumor founding clones (i.e., the ‘most recent common ancestor’ of the cells within a tumor) significantly distinguish the recurred and non-recurred breast tumors as well as outperform the most popular genomic test (i.e., Oncotype DX). These results imply that mutations of the tumor founding clones are associated with tumor recurrence and can be used to predict clinical outcomes. As such, predictive tools could be used in clinics to guide treatment routes. Finally, the concepts underlying the eTumorMetastasis pave the way for the application of genome sequencing in predictions for other complex genetic diseases. eTumorMetastasis pseudocode and related data used in this study are available at https://github.com/WangEdwinLab/eTumorMetastasis.     

Structure-Based Analysis Reveals Cancer Missense Mutations Target Protein Interaction Interfaces

Recently it has been shown that cancer mutations selectively target protein-protein interactions. We hypothesized that mutations affecting distinct protein interactions involving established cancer genes could contribute to tumor heterogeneity, and that novel mechanistic insights might be gained into tumorigenesis by investigating protein interactions under positive selection in cancer. To identify protein interactions under positive selection in cancer, we mapped over 1.2 million nonsynonymous somatic cancer mutations onto 4,896 experimentally determined protein structures and analyzed their spatial distribution. In total, 20% of mutations on the surface of known cancer genes perturbed protein-protein interactions (PPIs), and this enrichment for PPI interfaces was observed for both tumor suppressors (Odds Ratio 1.28, P-value < 10⁻⁴) and oncogenes (Odds Ratio 1.17, P-value < 10⁻³). To study this further, we constructed a bipartite network representing structurally resolved PPIs from all available human complexes in the Protein Data Bank (2,864 proteins, 3,072 PPIs). Analysis of frequently mutated cancer genes within this network revealed that tumor-suppressors, but not oncogenes, are significantly enriched with functional mutations in homo-oligomerization regions (Odds Ratio 3.68, P-Value < 10⁻⁸). We present two important examples, TP53 and beta-2-microglobulin, for which the patterns of somatic mutations at interfaces provide insights into specifically perturbed biological circuits. In patients with TP53 mutations, patient survival correlated with the specific interactions that were perturbed. Moreover, we investigated mutations at the interface of protein-nucleotide interactions and observed an unexpected number of missense mutations but not silent mutations occurring within DNA and RNA binding sites. Finally, we provide a resource of 3,072 PPI interfaces ranked according to their mutation rates. Analysis of this list highlights 282 novel candidate cancer genes that encode proteins participating in interactions that are perturbed recurrently across tumors. In summary, mutation of specific protein interactions is an important contributor to tumor heterogeneity and may have important implications for clinical outcomes.     

Evolution under Spatially Heterogeneous Selection in Solid Tumors

Spatial genetic and phenotypic diversity within solid tumors has been well documented. Nevertheless, how this heterogeneity affects temporal dynamics of tumorigenesis has not been rigorously examined because solid tumors do not evolve as the standard population genetic model due to the spatial constraint. We therefore, propose a neutral spatial (NS) model whereby the mutation accumulation increases toward the periphery; the genealogical relationship is spatially determined and the selection efficacy is blunted (due to kin competition). In this model, neutral mutations are accrued and spatially distributed in manners different from those of advantageous mutations. Importantly, the distinctions could be blurred in the conventional model. To test the NS model, we performed a three-dimensional multiple microsampling of two hepatocellular carcinomas. Whole-genome sequencing (WGS) revealed a 2-fold increase in mutations going from the center to the periphery. The operation of natural selection can then be tested by examining the spatially determined clonal relationships and the clonal sizes. Due to limited migration, only the expansion of highly advantageous clones can sweep through a large part of the tumor to reveal the selective advantages. Hence, even multiregional sampling can only reveal a fraction of fitness differences in solid tumors. Our results suggest that the NS patterns are crucial for testing the influence of natural selection during tumorigenesis, especially for small solid tumors.     

Intra-tumor heterogeneity of MLH1 promoter methylation revealed by deep single molecule bisulfite sequencing

A single tumor may contain cells with different somatic mutations. By characterizing this genetic heterogeneity within tumors, advances have been made in the prognosis, treatment and understanding of tumorigenesis. In contrast, the extent of epigenetic intra-tumor heterogeneity and how it influences tumor biology is under-explored. We have characterized epigenetic heterogeneity within individual tumors using next-generation sequencing. We used deep single molecule bisulfite sequencing and sample-specific DNA barcodes to determine the spectrum of MLH1 promoter methylation across an average of 1000 molecules in each of 33 individual samples in parallel, including endometrial cancer, matched blood and normal endometrium. This first glimpse, deep into each tumor, revealed unexpectedly heterogeneous patterns of methylation at the MLH1 promoter within a subset of endometrial tumors. This high-resolution analysis allowed us to measure the clonality of methylation in individual tumors and gain insight into the accumulation of aberrant promoter methylation on both alleles during tumorigenesis.     

Sequence Variation in the Gene Encoding the Nonstructural 3 Protein of Hepatitis C Virus: Evidence for Immune Selection

To determine whether the persistent nature of hepatitis C infection is related to the emergence of antigenic variants driven by immune selection, we examined the sequence heterogeneity in a portion of the hepatitis C virus (HCV) nonstructural 3 (NS3) gene of a patient infected over the course of more than 2 years. By PCR amplification, cloning, and sequencing, we observed several variable and conserved regions in the NS3 segment of the HCV genome. All variable regions had higher ratios of nonsynonymous/synonymous mutations and encompassed immunodominant epitopes, and their locations were not essential to maintain the known function of HCV RNA helicase. In contrast, the regions that are critical for HCV RNA helicase activity were found to be conserved with lower heterogeneity or lower ratios of nonsynonymous/synonymous mutations, and none except one of these regions was encoded within immunodominant epitopes. Our results are consistent with immune selection of viral variants at the epitope and molecular levels that may enable HCV to evade host defenses over time. Plotting the relatedness of sequence variants revealed a star topology suggesting that a wild-type HCV sequence is maintained, unlike HIV. Received: 2 November 2000 / Accepted: 1 October 2001     

Whole mitochondrial DNA variations in hippocampal surgical specimens and blood samples with high-throughput sequencing: A case of mesial temporal lobe epilepsy with hippocampal sclerosis

Introduction

Hippocampal sclerosis is the most common lesion in patients with mesial temporal lobe epilepsy. Recently, there has been growing evidence on the involvement of mitochondria also in sporadic forms of epilepsy. In addition, it has been increasingly argued that mitochondrial dysfunction has an important role in epileptogenesis and seizure generation in temporal lobe epilepsy. Although mtDNA polymorphisms have been identified as potential risk factors for neurological diseases, the link between homoplasmy and heteroplasmy within tissues is not clear. We investigated whether mitochondrial DNA (mtDNA) polymorphisms are involved in a case report of a patient with mesial temporal lobe epilepsy-hippocampal sclerosis (MTLE-HS).

Design

We report the whole genome mtDNA deep sequencing results and clinical features of a 36-year-old woman with MTLE-HS. We used pyrosequencing technology to sequence a whole mitochondrial genome isolated from six different regions of her brain and blood. To assess the possible role of mitochondrial DNA variations in affected tissues, we compared all specimens from different regions of the hippocampus and blood.

Results

In total, 35 homoplasmic and 18 heteroplasmic variations have been detected in 6 different regions of the hippocampus and in blood samples. While the samples did not display any difference in homoplasmic variations, it has been shown that hippocampus regions contain more heteroplasmic variations than blood. The number of heteroplasmic variations was highest in the CA2 region of the brain and accumulated in ND2, ND4 and ND5 genes. Also, dentate and subiculum regions of the hippocampus had similar heteroplasmic variation profiles.

Discussion

We present a new rare example of parallel mutation at 16223 position. Our case suggests that defects in mitochondrial function might be underlying the pathogenesis of seizures in temporal lobe epilepsy.     

Identity by descent mapping of founder mutations in cancer using high-resolution tumor SNP data

Dense genotype data can be used to detect chromosome fragments inherited from a common ancestor in apparently unrelated individuals. A disease-causing mutation inherited from a common founder may thus be detected by searching for a common haplotype signature in a sample population of patients. We present here FounderTracker, a computational method for the genome-wide detection of founder mutations in cancer using dense tumor SNP profiles. Our method is based on two assumptions. First, the wild-type allele frequently undergoes loss of heterozygosity (LOH) in the tumors of germline mutation carriers. Second, the overlap between the ancestral chromosome fragments inherited from a common founder will define a minimal haplotype conserved in each patient carrying the founder mutation. Our approach thus relies on the detection of haplotypes with significant identity by descent (IBD) sharing within recurrent regions of LOH to highlight genomic loci likely to harbor a founder mutation. We validated this approach by analyzing two real cancer data sets in which we successfully identified founder mutations of well-characterized tumor suppressor genes. We then used simulated data to evaluate the ability of our method to detect IBD tracts as a function of their size and frequency. We show that FounderTracker can detect haplotypes of low prevalence with high power and specificity, significantly outperforming existing methods. FounderTracker is thus a powerful tool for discovering unknown founder mutations that may explain part of the "missing" heritability in cancer. This method is freely available and can be used online at the FounderTracker website.     

Aphidicolin resistance in herpes simplex virus type I reveals features of the DNA polymerase dNTP binding site

We describe the mapping and sequencing of mutations within the DNA polymerase gene of herpes simplex virus type 1 which confer resistance to aphidicolin, a DNA polymerase inhibitor. The mutations occur near two regions which are highly conserved among DNA polymerases related to the herpes simplex enzyme. They also occur near other herpes simplex mutations which affect the interactions between the polymerase and deoxyribonucleoside triphosphate substrates. Consequently, we argue in favor of the idea that the aphidicolin binding site overlaps the substrate binding site and that the near-by conserved regions are functionally required for substrate binding. Our mutants also exhibit abnormal sensitivity to another DNA polymerase inhibitor, phosphonoacetic acid. This drug is thought to bind as an analogue of pyrophosphate. A second-site mutation which suppresses the hypersensitivity of one mutant to phosphonoacetic acid (but not its aphidicolin resistance) is described. This second mutation may represent a new class of mutations, which specifically affects pyrophosphate, but not substrate, binding.     

Exome Sequencing of Cell-Free DNA from Metastatic Cancer Patients Identifies Clinically Actionable Mutations Distinct from Primary Disease

The identification of the molecular drivers of cancer by sequencing is the backbone of precision medicine and the basis of personalized therapy; however, biopsies of primary tumors provide only a snapshot of the evolution of the disease and may miss potential therapeutic targets, especially in the metastatic setting. A liquid biopsy, in the form of cell-free DNA (cfDNA) sequencing, has the potential to capture the inter- and intra-tumoral heterogeneity present in metastatic disease, and, through serial blood draws, track the evolution of the tumor genome.In order to determine the clinical utility of cfDNA sequencing we performed whole-exome sequencing on cfDNA and tumor DNA from two patients with metastatic disease; only minor modifications to our sequencing and analysis pipelines were required for sequencing and mutation calling of cfDNA. The first patient had metastatic sarcoma and 47 of 48 mutations present in the primary tumor were also found in the cell-free DNA. The second patient had metastatic breast cancer and sequencing identified an ESR1 mutation in the cfDNA and metastatic site, but not in the primary tumor. This likely explains tumor progression on Anastrozole. Significant heterogeneity between the primary and metastatic tumors, with cfDNA reflecting the metastases, suggested separation from the primary lesion early in tumor evolution. This is best illustrated by an activating PIK3CA mutation (H1047R) which was clonal in the primary tumor, but completely absent from either the metastasis or cfDNA. Here we show that cfDNA sequencing supplies clinically actionable information with minimal risks compared to metastatic biopsies. This study demonstrates the utility of whole-exome sequencing of cell-free DNA from patients with metastatic disease. cfDNA sequencing identified an ESR1 mutation, potentially explaining a patient’s resistance to aromatase inhibition, and gave insight into how metastatic lesions differ from the primary tumor.     

A case of interdigitating dendritic cell sarcoma studied by whole-exome sequencing

Interdigitating dendritic cell sarcoma (IDCS) is an aggressive neoplasm and is an extremely rare disease, with a challenging diagnosis. Etiology of IDCS is also unknown and most studies with only case reports. In our case, immunohistochemistry showed that the tumor cells were positive for S100, CD45, and CD68, but negative for CD1a and CD21. This study aimed to investigate the causative factors of IDCS by sequencing the protein-coding regions of IDCS. We performed whole-exome sequencing with genomic DNA from blood and sarcoma tissue of the IDCS patient using the Illumina Hiseq 2500 platform. After that, we conducted Sanger sequencing for validation of sarcoma-specific variants and gene ontology analysis using DAVID bioinformatics resources. Through comparing sequencing data of sarcoma with normal blood, we obtained 15 nonsynonymous single nucleotide polymorphisms (SNPs) as sarcoma-specific variants. Although the 15 SNPs were not validated by Sanger sequencing due to tumor heterogeneity and low sensitivity of Sanger sequencing, we examined the function of the genes in which each SNP is located. Based on previous studies and gene ontology database, we found that POLQ encoding DNA polymerase theta enzyme and FNIP1 encoding tumor suppressor folliculin-interacting protein might have contributed to the IDCS. Our study provides potential causative genetic factors of IDCS and plays a role in advancing the understanding of IDCS pathogenesis.     

Loss-of-function mutations in PTPN11 cause metachondromatosis, but not Ollier disease or Maffucci syndrome

Metachondromatosis (MC) is a rare, autosomal dominant, incompletely penetrant combined exostosis and enchondromatosis tumor syndrome. MC is clinically distinct from other multiple exostosis or multiple enchondromatosis syndromes and is unlinked to EXT1 and EXT2, the genes responsible for autosomal dominant multiple osteochondromas (MO). To identify a gene for MC, we performed linkage analysis with high-density SNP arrays in a single family, used a targeted array to capture exons and promoter sequences from the linked interval in 16 participants from 11 MC families, and sequenced the captured DNA using high-throughput parallel sequencing technologies. DNA capture and parallel sequencing identified heterozygous putative loss-of-function mutations in PTPN11 in 4 of the 11 families. Sanger sequence analysis of PTPN11 coding regions in a total of 17 MC families identified mutations in 10 of them (5 frameshift, 2 nonsense, and 3 splice-site mutations). Copy number analysis of sequencing reads from a second targeted capture that included the entire PTPN11 gene identified an additional family with a 15 kb deletion spanning exon 7 of PTPN11. Microdissected MC lesions from two patients with PTPN11 mutations demonstrated loss-of-heterozygosity for the wild-type allele. We next sequenced PTPN11 in DNA samples from 54 patients with the multiple enchondromatosis disorders Ollier disease or Maffucci syndrome, but found no coding sequence PTPN11 mutations. We conclude that heterozygous loss-of-function mutations in PTPN11 are a frequent cause of MC, that lesions in patients with MC appear to arise following a "second hit," that MC may be locus heterogeneous since 1 familial and 5 sporadically occurring cases lacked obvious disease-causing PTPN11 mutations, and that PTPN11 mutations are not a common cause of Ollier disease or Maffucci syndrome.