Steven G. Rozen

Intrinsic Subtypes of Gastric Cancer, Based on Gene Expression Pattern, Predict Survival and Respond Differently to Chemotherapy

BACKGROUND & AIMS Gastric cancer (GC) is a heterogeneous disease comprising multiple subtypes that have distinct biological properties and effects in patients. We sought to identify new, intrinsic subtypes of GC by gene expression analysis of a large panel of GC cell lines. We tested if these subtypes might be associated with differences in patient survival times and responses to various standard-of-care cytotoxic drugs. METHODS We analyzed gene expression profiles for 37 GC cell lines to identify intrinsic GC subtypes. These subtypes were validated in primary tumors from 521 patients in 4 independent cohorts, where the subtypes were determined by either expression profiling or subtype-specific immunohistochemical markers (LGALS4, CDH17). In vitro sensitivity to 3 chemotherapy drugs (5-fluorouracil, cisplatin, oxaliplatin) was also assessed. RESULTS Unsupervised cell line analysis identified 2 major intrinsic genomic subtypes (G-INT and G-DIF) that had distinct patterns of gene expression. The intrinsic subtypes, but not subtypes based on Laurens histopathologic classification, were prognostic of survival, based on univariate and multivariate analysis in multiple patient cohorts. The G-INT cell lines were significantly more sensitive to 5-fluorouracil and oxaliplatin, but more resistant to cisplatin, than the G-DIF cell lines. In patients, intrinsic subtypes were associated with survival time following adjuvant, 5-fluorouracil-based therapy. CONCLUSIONS Intrinsic subtypes of GC, based on distinct patterns of expression, are associated with patient survival and response to chemotherapy. Classification of GC based on intrinsic subtypes might be used to determine prognosis and customize therapy.

Exome sequencing identifies distinct mutational patterns in liver fluke–related and non-infection-related bile duct cancers

The impact of different carcinogenic exposures on the specific patterns of somatic mutation in human tumors remains unclear. To address this issue, we profiled 209 cholangiocarcinomas (CCAs) from Asia and Europe, including 108 cases caused by infection with the liver fluke Opisthorchis viverrini and 101 cases caused by non–O. viverrini–related etiologies. Whole-exome sequencing (n = 15) and prevalence screening (n = 194) identified recurrent somatic mutations in BAP1 and ARID1A, neither of which, to our knowledge, has previously been reported to be mutated in CCA. Comparisons between intrahepatic O. viverrini–related and non–O. viverrini–related CCAs demonstrated statistically significant differences in mutation patterns: BAP1, IDH1 and IDH2 were more frequently mutated in non–O. viverrini CCAs, whereas TP53 mutations showed the reciprocal pattern. Functional studies demonstrated tumor suppressive functions for BAP1 and ARID1A, establishing the role of chromatin modulators in CCA pathogenesis. These findings indicate that different causative etiologies may induce distinct somatic alterations, even within the same tumor type.

Genome-Wide Mutational Signatures of Aristolochic Acid and Its Application as a Screening Tool

Genome-wide mutational signatures of the group 1 carcinogen aristolochic acid are observed in urothelial cancers and liver cancers from Asia. Carcinogen AAlert Aristolochic acid (AA) is a natural compound derived from plants in the Aristolochia genus. For centuries, Aristolochia has been used throughout Asia to treat a variety of ailments as a component of traditional Chinese medicine. In recent years, however, a more sinister side of this herb has come to light when it was linked to kidney damage and cancers of the urinary tract. Now, two studies by Poon et al. and Hoang et al. present a “molecular signature” of AA-induced DNA damage, which helps to explain the mutagenic effects of AA and may also be useful as a way to detect unsuspected AA exposure as a cause of cancer. The molecular signature seen in AA-associated tumors is characterized by a predominance of A:T-to-T:A transversions, a relatively unusual type of mutation that is infrequently seen in other types of cancer, including those caused by other carcinogens. These mutations concentrate at splice sites, causing the inappropriate inclusion or exclusion of entire exons in the resulting mRNA. The overall mutation rate is another notable feature of AA-associated cancers because it is several times higher than the rate of mutations caused by other carcinogens such as tobacco and ultraviolet light. In both studies, the authors also used the molecular signature to discover that AA was a likely cause of tumors previously attributed to other carcinogens. In one case, a urinary tract cancer that had been attributed to smoking and, in the other case, a liver cancer previously attributed to a chronic hepatitis infection were both identified as having the telltale signature of AA mutagenesis. The identification of a specific molecular signature for AA has both clinical and public health implications. For individual patients, the molecular signature could help physicians identify which tumors were caused by AA. Although this information cannot yet be used to optimize the treatment of individual patients, those who are diagnosed with AA-associated cancers could be monitored more closely for the appearance of additional tumors. Meanwhile, a better understanding of the mutagenic effects of AA should also help to strengthen public health efforts to decrease exposure to this carcinogenic herb. Aristolochic acid (AA), a natural product of Aristolochia plants found in herbal remedies and health supplements, is a group 1 carcinogen that can cause nephrotoxicity and upper urinary tract urothelial cell carcinoma (UTUC). Whole-genome and exome analysis of nine AA-associated UTUCs revealed a strikingly high somatic mutation rate (150 mutations/Mb), exceeding smoking-associated lung cancer (8 mutations/Mb) and ultraviolet radiation–associated melanoma (111 mutations/Mb). The AA-UTUC mutational signature was characterized by A:T to T:A transversions at the sequence motif A[C|T]AGG, located primarily on nontranscribed strands. AA-induced mutations were also significantly enriched at splice sites, suggesting a role for splice-site mutations in UTUC pathogenesis. RNA sequencing of AA-UTUC confirmed a general up-regulation of nonsense-mediated decay machinery components and aberrant splicing events associated with splice-site mutations. We observed a high frequency of somatic mutations in chromatin modifiers, particularly KDM6A, in AA-UTUC, demonstrated the sufficiency of AA to induce renal dysplasia in mice, and reproduced the AA mutational signature in experimentally treated human renal tubular cells. Finally, exploring other malignancies that were not known to be associated with AA, we screened 93 hepatocellular carcinoma genomes/exomes and identified AA-like mutational signatures in 11. Our study highlights an unusual genome-wide AA mutational signature and the potential use of mutation signatures as “molecular fingerprints” for interrogating high-throughput cancer genome data to infer previous carcinogen exposures.

Exome sequencing identifies highly recurrent MED12 somatic mutations in breast fibroadenoma

Fibroadenomas are the most common breast tumors in women under 30 (refs. 1,2). Exome sequencing of eight fibroadenomas with matching whole-blood samples revealed recurrent somatic mutations solely in MED12, which encodes a Mediator complex subunit. Targeted sequencing of an additional 90 fibroadenomas confirmed highly frequent MED12 exon 2 mutations (58/98, 59%) that are probably somatic, with 71% of mutations occurring in codon 44. Using laser capture microdissection, we show that MED12 fibroadenoma mutations are present in stromal but not epithelial mammary cells. Expression profiling of MED12-mutated and wild-type fibroadenomas revealed that MED12 mutations are associated with dysregulated estrogen signaling and extracellular matrix organization. The fibroadenoma MED12 mutation spectrum is nearly identical to that of previously reported MED12 lesions in uterine leiomyoma but not those of other tumors. Benign tumors of the breast and uterus, both of which are key target tissues of estrogen, may thus share a common genetic basis underpinned by highly frequent and specific MED12 mutations.

Whole-genome reconstruction and mutational signatures in gastric cancer

BackgroundGastric cancer is the second highest cause of global cancer mortality. To explore the complete repertoire of somatic alterations in gastric cancer, we combined massively parallel short read and DNA paired-end tag sequencing to present the first whole-genome analysis of two gastric adenocarcinomas, one with chromosomal instability and the other with microsatellite instability.ResultsIntegrative analysis and de novo assemblies revealed the architecture of a wild-type KRAS amplification, a common driver event in gastric cancer. We discovered three distinct mutational signatures in gastric cancer - against a genome-wide backdrop of oxidative and microsatellite instability-related mutational signatures, we identified the first exome-specific mutational signature. Further characterization of the impact of these signatures by combining sequencing data from 40 complete gastric cancer exomes and targeted screening of an additional 94 independent gastric tumors uncovered ACVR2A, RPL22 and LMAN1 as recurrently mutated genes in microsatellite instability-positive gastric cancer and PAPPA as a recurrently mutated gene in TP53 wild-type gastric cancer.ConclusionsThese results highlight how whole-genome cancer sequencing can uncover information relevant to tissue-specific carcinogenesis that would otherwise be missed from exome-sequencing data.

Extensive diversity in circadian regulation of plasma lipids and evidence for different circadian metabolic phenotypes in humans

The circadian system regulates daily rhythms in lipid metabolism and adipose tissue function. Although disruption of circadian clock function is associated with negative cardiometabolic end points, very little is known about interindividual variation in circadian-regulated metabolic pathways. Here, we used targeted lipidomics-based approaches to profile the time course of 263 lipids in blood plasma in 20 healthy individuals. Over a span of 28 h, blood was collected every 4 h and plasma lipids were analyzed by HPLC/MS. Across subjects, about 13% of lipid metabolites showed circadian variation. Rhythmicity spanned all metabolite classes examined, suggesting widespread circadian control of lipid-mediated energy storage, transport, and signaling. Intersubject agreement for lipids identified as rhythmic was only about 20%, however, and the timing of lipid rhythms ranged up to 12 h apart between individuals. Healthy subjects therefore showed substantial variation in the timing and strength of rhythms across different lipid species. Strong interindividual differences were also observed for rhythms of blood glucose and insulin, but not cortisol. Using consensus clustering with iterative feature selection, subjects clustered into different groups based on strength of rhythmicity for a subset of triglycerides and phosphatidylcholines, suggesting that there are different circadian metabolic phenotypes in the general population. These results have potential implications for lipid metabolism disorders linked to circadian clock disruption.

Genomic landscapes of breast fibroepithelial tumors

Breast fibroepithelial tumors comprise a heterogeneous spectrum of pathological entities, from benign fibroadenomas to malignant phyllodes tumors. Although MED12 mutations have been frequently found in fibroadenomas and phyllodes tumors, the landscapes of genetic alterations across the fibroepithelial tumor spectrum remain unclear. Here, by performing exome sequencing of 22 phyllodes tumors followed by targeted sequencing of 100 breast fibroepithelial tumors, we observed three distinct somatic mutation patterns. First, we frequently observed MED12 and RARA mutations in both fibroadenomas and phyllodes tumors, emphasizing the importance of these mutations in fibroepithelial tumorigenesis. Second, phyllodes tumors exhibited mutations in FLNA, SETD2 and KMT2D, suggesting a role in driving phyllodes tumor development. Third, borderline and malignant phyllodes tumors harbored additional mutations in cancer-associated genes. RARA mutations exhibited clustering in the portion of the gene encoding the ligand-binding domain, functionally suppressed RARA-mediated transcriptional activation and enhanced RARA interactions with transcriptional co-repressors. This study provides insights into the molecular pathogenesis of breast fibroepithelial tumors, with potential clinical implications.

Mutation signatures of carcinogen exposure: genome-wide detection and new opportunities for cancer prevention.

Exposure to environmental mutagens is an important cause of human cancer, and measures to reduce mutagenic and carcinogenic exposures have been highly successful at controlling cancer. Until recently, it has been possible to connect the chemical characteristics of mutagens to actual mutations observed in human tumors only indirectly. Now, next-generation sequencing technology enables us to observe in detail the DNA-sequence-level effects of well-known mutagens, such as ultraviolet radiation and tobacco smoke, as well as endogenous mutagenic processes, such as those involving activated DNA cytidine deaminases (APOBECs). We can also observe the effects of less well-known but potent mutagens, including those recently found to be present in some herbal remedies. Crucially, we can now tease apart the superimposed effects of several mutational exposures and processes and determine which ones occurred during the development of individual tumors. Here, we review advances in detecting these mutation signatures and discuss the implications for surveillance and prevention of cancer. The number of sequenced tumors from diverse cancer types and multiple geographic regions is growing explosively, and the genomes of these tumors will bear the signatures of even more diverse mutagenic exposures. Thus, we envision development of wide-ranging compendia of mutation signatures from tumors and a concerted effort to experimentally elucidate the signatures of a large number of mutagens. This information will be used to link signatures observed in tumors to the exposures responsible for them, which will offer unprecedented opportunities for prevention.

Evaluation and optimisation of indel detection workflows for ion torrent sequencing of the BRCA1 and BRCA2 genes

BackgroundThe Ion Torrent PGM is a popular benchtop sequencer that shows promise in replacing conventional Sanger sequencing as the gold standard for mutation detection. Despite the PGM’s reported high accuracy in calling single nucleotide variations, it tends to generate many false positive calls in detecting insertions and deletions (indels), which may hinder its utility for clinical genetic testing.ResultsRecently, the proprietary analytical workflow for the Ion Torrent sequencer, Torrent Suite (TS), underwent a series of upgrades. We evaluated three major upgrades of TS by calling indels in the BRCA1 and BRCA2 genes. Our analysis revealed that false negative indels could be generated by TS under both default calling parameters and parameters adjusted for maximum sensitivity. However, indel calling with the same data using the open source variant callers, GATK and SAMtools showed that false negatives could be minimised with the use of appropriate bioinformatics analysis. Furthermore, we identified two variant calling measures, Quality-by-Depth (QD) and VARiation of the Width of gaps and inserts (VARW), which substantially reduced false positive indels, including non-homopolymer associated errors without compromising sensitivity. In our best case scenario that involved the TMAP aligner and SAMtools, we achieved 100% sensitivity, 99.99% specificity and 29% False Discovery Rate (FDR) in indel calling from all 23 samples, which is a good performance for mutation screening using PGM.ConclusionsNew versions of TS, BWA and GATK have shown improvements in indel calling sensitivity and specificity over their older counterpart. However, the variant caller of TS exhibits a lower sensitivity than GATK and SAMtools. Our findings demonstrate that although indel calling from PGM sequences may appear to be noisy at first glance, proper computational indel calling analysis is able to maximize both the sensitivity and specificity at the single base level, paving the way for the usage of this technology for future clinical genetic testing.

Nucleus-Translocated ACSS2 Promotes Gene Transcription for Lysosomal Biogenesis and Autophagy

Overcoming metabolic stress is a critical step in tumor growth. Acetyl coenzyme A (acetyl-CoA) generated from glucose and acetate uptake is important for histone acetylation and gene expression. However, how acetyl-CoA is produced under nutritional stress is unclear. We demonstrate here that glucose deprivation results in AMP-activated protein kinase (AMPK)-mediated acetyl-CoA synthetase 2 (ACSS2) phosphorylation at S659, which exposed the nuclear localization signal of ACSS2 for importin α5 binding and nuclear translocation. In the nucleus, ACSS2 binds to transcription factor EB and translocates to lysosomal and autophagy gene promoter regions, where ACSS2 incorporates acetate generated from histone acetylation turnover to locally produce acetyl-CoA for histone H3 acetylation in these regions and promote lysosomal biogenesis, autophagy, cell survival, and brain tumorigenesis. In addition, ACSS2 S659 phosphorylation positively correlates with AMPK activity in glioma specimens and grades of glioma malignancy. These results underscore the significance of nuclear ACSS2-mediated histone acetylation in maintaining cell homeostasis and tumor development.