ABSTRACT
BACKGROUND & AIMS: Primary Biliary Cholangitis (PBC) is a chronic autoimmune liver disease characterized by destruction of intrahepatic bile ducts, portal inflammation and cirrhosis. Although rare in most populations, it is prevalent and often familial in British Columbia First Nations. We hypothesized that major genetic factors increased the risk in First Nations. METHODS: In all, 44 individuals with Primary Biliary Cholangitis and 61 unaffected relatives from 32 First Nations families participated. Family history and co-morbidities were documented. Medical records were reviewed and available biopsies were re-reviewed by our team pathologist. Genotyping was performed on DNA from 36 affected persons and 27 unaffected relatives using the Affymetrix Human Mapping 500K Array Set. MERLIN software was used to carry out multipoint parametric and nonparametric linkage analysis. Candidate genes were identified and entered into InnateDB and KEGG software to identify potential pathways affecting pathogenesis. RESULTS: In all, 34% of families were multiplex. Fifty per cent of cases and 33% of unaffected relatives reported other autoimmune disease. Three genomic regions (9q21, 17p13 and 19p13) produced LOD scores of 2.3 or greater suggestive of linkage, but no single linkage peak reached statistical significance. Candidate genes identified in the three regions suggested involvement of IL17, NFκB, IL6, JAK-STAT, IFNγ and TGFß immune signalling pathways. Specifically, four genes-ACT1, PIN1, DNMT1 and NTN1-emerged as having roles in these pathways that may influence Primary Biliary Cholangitis pathogenesis. CONCLUSIONS: Our whole genome linkage study results reflect the multifactorial nature of Primary Biliary Cholangitis, support previous studies suggesting signalling pathway involvement and identify new candidate genes for consideration.
Subject(s)
Genetic Predisposition to Disease , Indians, North American/genetics , Liver Cirrhosis, Biliary/ethnology , Liver Cirrhosis, Biliary/genetics , Adaptor Proteins, Signal Transducing , Adult , Autoimmune Diseases/ethnology , Autoimmune Diseases/genetics , British Columbia/ethnology , Case-Control Studies , DNA (Cytosine-5-)-Methyltransferase 1/genetics , Female , Humans , Lod Score , Male , Middle Aged , NIMA-Interacting Peptidylprolyl Isomerase/genetics , Netrin-1/genetics , Tumor Necrosis Factor Receptor-Associated Peptides and Proteins/geneticsABSTRACT
Characterization and quantification of tumour clonal populations over time via longitudinal sampling are essential components in understanding and predicting the response to therapeutic interventions. Computational methods for inferring tumour clonal composition from deep-targeted sequencing data are ubiquitous, however due to the lack of a ground truth biological data, evaluating their performance is difficult. In this work, we generate a benchmark data set that simulates tumour longitudinal growth and heterogeneity by in vitro mixing of cancer cell lines with known proportions. We apply four different algorithms to our ground truth data set and assess their performance in inferring clonal composition using different metrics. We also analyse the performance of these algorithms on breast tumour xenograft samples. We conclude that methods that can simultaneously analyse multiple samples while accounting for copy number alterations as a factor in allelic measurements exhibit the most accurate predictions. These results will inform future functional genomics oriented studies of model systems where time series measurements in the context of therapeutic interventions are becoming increasingly common. These studies will need computational models which accurately reflect the multi-factorial nature of allele measurement in cancer including, as we show here, segmental aneuploidies.
Subject(s)
Computer Simulation , Models, Biological , Neoplasms/etiology , Neoplasms/pathology , Algorithms , Animals , Breast Neoplasms/etiology , Breast Neoplasms/pathology , Cell Line, Tumor , Computational Biology/methods , DNA Copy Number Variations , Disease Models, Animal , Female , Heterografts , Humans , Mice , Polymorphism, Single Nucleotide , Reproducibility of Results , Exome SequencingABSTRACT
We studied the whole-genome point mutation and structural variation patterns of 133 tumors (59 high-grade serous (HGSC), 35 clear cell (CCOC), 29 endometrioid (ENOC), and 10 adult granulosa cell (GCT)) as a substrate for class discovery in ovarian cancer. Ab initio clustering of integrated point mutation and structural variation signatures identified seven subgroups both between and within histotypes. Prevalence of foldback inversions identified a prognostically significant HGSC group associated with inferior survival. This finding was recapitulated in two independent cohorts (n = 576 cases), transcending BRCA1 and BRCA2 mutation and gene expression features of HGSC. CCOC cancers grouped according to APOBEC deamination (26%) and age-related mutational signatures (40%). ENOCs were divided by cases with microsatellite instability (28%), with a distinct mismatch-repair mutation signature. Taken together, our work establishes the potency of the somatic genome, reflective of diverse DNA repair deficiencies, to stratify ovarian cancers into distinct biological strata within the major histotypes.
Subject(s)
DNA Repair/genetics , Ovarian Neoplasms/genetics , Ovarian Neoplasms/pathology , BRCA1 Protein/genetics , BRCA2 Protein/genetics , Endometriosis/complications , Endometriosis/genetics , Female , Gene Expression Regulation, Neoplastic , Genome, Human , Humans , Mutation , Ovarian Neoplasms/drug therapy , Ovarian Neoplasms/mortality , PrognosisABSTRACT
CDC-like kinase phosphorylation of serine/arginine-rich proteins is central to RNA splicing reactions. Yet, the genomic network of CDC-like kinase-dependent RNA processing events remains poorly defined. Here, we explore the connectivity of genomic CDC-like kinase splicing functions by applying graduated, short-exposure, pharmacological CDC-like kinase inhibition using a novel small molecule (T3) with very high potency, selectivity, and cell-based stability. Using RNA-Seq, we define CDC-like kinase-responsive alternative splicing events, the large majority of which monotonically increase or decrease with increasing CDC-like kinase inhibition. We show that distinct RNA-binding motifs are associated with T3 response in skipped exons. Unexpectedly, we observe dose-dependent conjoined gene transcription, which is associated with motif enrichment in the last and second exons of upstream and downstream partners, respectively. siRNA knockdown of CLK2-associated genes significantly increases conjoined gene formation. Collectively, our results reveal an unexpected role for CDC-like kinase in conjoined gene formation, via regulation of 3'-end processing and associated splicing factors.The phosphorylation of serine/arginine-rich proteins by CDC-like kinase is a central regulatory mechanism for RNA splicing reactions. Here, the authors synthesize a novel small molecule CLK inhibitor and map CLK-responsive alternative splicing events and discover an effect on conjoined gene transcription.
Subject(s)
Alternative Splicing/drug effects , Imidazoles/pharmacology , Protein Kinase Inhibitors/pharmacology , Protein Serine-Threonine Kinases/genetics , Protein-Tyrosine Kinases/genetics , Pyrimidines/pharmacology , RNA, Messenger/genetics , RNA-Binding Proteins/genetics , Exons , Gene Expression Profiling , Genome, Human , HCT116 Cells , Humans , Imidazoles/chemical synthesis , Phosphorylation/drug effects , Protein Kinase Inhibitors/chemical synthesis , Protein Serine-Threonine Kinases/antagonists & inhibitors , Protein Serine-Threonine Kinases/metabolism , Protein-Tyrosine Kinases/antagonists & inhibitors , Protein-Tyrosine Kinases/metabolism , Pyrimidines/chemical synthesis , RNA, Messenger/antagonists & inhibitors , RNA, Messenger/metabolism , RNA, Small Interfering/genetics , RNA, Small Interfering/metabolism , RNA-Binding Proteins/antagonists & inhibitors , RNA-Binding Proteins/metabolism , Structure-Activity Relationship , Transcription, GeneticABSTRACT
BACKGROUND: Follicular lymphoma (FL) is an indolent, yet incurable B cell malignancy. A subset of patients experience an increased mortality rate driven by two distinct clinical end points: histological transformation and early progression after immunochemotherapy. The nature of tumor clonal dynamics leading to these clinical end points is poorly understood, and previously determined genetic alterations do not explain the majority of transformed cases or accurately predict early progressive disease. We contend that detailed knowledge of the expansion patterns of specific cell populations plus their associated mutations would provide insight into therapeutic strategies and disease biology over the time course of FL clinical histories. METHODS AND FINDINGS: Using a combination of whole genome sequencing, targeted deep sequencing, and digital droplet PCR on matched diagnostic and relapse specimens, we deciphered the constituent clonal populations in 15 transformation cases and 6 progression cases, and measured the change in clonal population abundance over time. We observed widely divergent patterns of clonal dynamics in transformed cases relative to progressed cases. Transformation specimens were generally composed of clones that were rare or absent in diagnostic specimens, consistent with dramatic clonal expansions that came to dominate the transformation specimens. This pattern was independent of time to transformation and treatment modality. By contrast, early progression specimens were composed of clones that were already present in the diagnostic specimens and exhibited only moderate clonal dynamics, even in the presence of immunochemotherapy. Analysis of somatic mutations impacting 94 genes was undertaken in an extension cohort consisting of 395 samples from 277 patients in order to decipher disrupted biology in the two clinical end points. We found 12 genes that were more commonly mutated in transformed samples than in the preceding FL tumors, including TP53, B2M, CCND3, GNA13, S1PR2, and P2RY8. Moreover, ten genes were more commonly mutated in diagnostic specimens of patients with early progression, including TP53, BTG1, MKI67, and XBP1. CONCLUSIONS: Our results illuminate contrasting modes of evolution shaping the clinical histories of transformation and progression. They have implications for interpretation of evolutionary dynamics in the context of treatment-induced selective pressures, and indicate that transformation and progression will require different clinical management strategies.
Subject(s)
Clonal Evolution , Disease Progression , Lymphoma, Follicular/physiopathology , Clone Cells , Humans , Lymphoma, Follicular/genetics , MutationABSTRACT
We performed phylogenetic analysis of high-grade serous ovarian cancers (68 samples from seven patients), identifying constituent clones and quantifying their relative abundances at multiple intraperitoneal sites. Through whole-genome and single-nucleus sequencing, we identified evolutionary features including mutation loss, convergence of the structural genome and temporal activation of mutational processes that patterned clonal progression. We then determined the precise clonal mixtures comprising each tumor sample. The majority of sites were clonally pure or composed of clones from a single phylogenetic clade. However, each patient contained at least one site composed of polyphyletic clones. Five patients exhibited monoclonal and unidirectional seeding from the ovary to intraperitoneal sites, and two patients demonstrated polyclonal spread and reseeding. Our findings indicate that at least two distinct modes of intraperitoneal spread operate in clonal dissemination and highlight the distribution of migratory potential over clonal populations comprising high-grade serous ovarian cancers.
Subject(s)
Biomarkers, Tumor/genetics , Clone Cells/pathology , Cystadenocarcinoma, Serous/pathology , Genetic Variation/genetics , Ovarian Neoplasms/pathology , Peritoneal Neoplasms/pathology , Tumor Microenvironment/genetics , Aged , Clone Cells/metabolism , Cystadenocarcinoma, Serous/genetics , Disease Progression , Fallopian Tube Neoplasms/genetics , Fallopian Tube Neoplasms/pathology , Female , Gene Expression Profiling , Gene Expression Regulation, Neoplastic , Genome, Human , High-Throughput Nucleotide Sequencing/methods , Humans , Middle Aged , Mutation/genetics , Neoplasm Grading , Neoplasm Recurrence, Local/genetics , Neoplasm Recurrence, Local/pathology , Ovarian Neoplasms/genetics , Peritoneal Neoplasms/genetics , Phylogeny , Single-Cell Analysis/methods , Survival RateABSTRACT
The development of targeted anti-cancer therapies through the study of cancer genomes is intended to increase survival rates and decrease treatment-related toxicity. We treated a transposon-driven, functional genomic mouse model of medulloblastoma with 'humanized' in vivo therapy (microneurosurgical tumour resection followed by multi-fractionated, image-guided radiotherapy). Genetic events in recurrent murine medulloblastoma exhibit a very poor overlap with those in matched murine diagnostic samples (<5%). Whole-genome sequencing of 33 pairs of human diagnostic and post-therapy medulloblastomas demonstrated substantial genetic divergence of the dominant clone after therapy (<12% diagnostic events were retained at recurrence). In both mice and humans, the dominant clone at recurrence arose through clonal selection of a pre-existing minor clone present at diagnosis. Targeted therapy is unlikely to be effective in the absence of the target, therefore our results offer a simple, proximal, and remediable explanation for the failure of prior clinical trials of targeted therapy.
Subject(s)
Cerebellar Neoplasms/therapy , Clone Cells/drug effects , Clone Cells/metabolism , Medulloblastoma/therapy , Neoplasm Recurrence, Local/genetics , Neoplasm Recurrence, Local/pathology , Selection, Genetic/drug effects , Animals , Cerebellar Neoplasms/genetics , Cerebellar Neoplasms/pathology , Cerebellar Neoplasms/radiotherapy , Cerebellar Neoplasms/surgery , Clone Cells/pathology , Craniospinal Irradiation , DNA Mutational Analysis , Disease Models, Animal , Drosophila melanogaster/cytology , Drosophila melanogaster/genetics , Female , Genome, Human/genetics , Humans , Male , Medulloblastoma/genetics , Medulloblastoma/pathology , Medulloblastoma/radiotherapy , Medulloblastoma/surgery , Mice , Molecular Targeted Therapy/methods , Neoplasm Recurrence, Local/therapy , Radiotherapy, Image-Guided , Signal Transduction , Xenograft Model Antitumor AssaysABSTRACT
We carried out whole genome and transcriptome sequencing on four tumour/normal pairs of epithelioid sarcoma. These index cases were supplemented with whole transcriptome sequencing of three additional tumours and three cell lines. Unlike rhabdoid tumour (the other major group of SMARCB1-negative cancers), epithelioid sarcoma shows a complex genome with a higher mutational rate, comparable to that of ovarian carcinoma. Despite this mutational burden, SMARCB1 mutations remain the most frequently recurring event and are probably critical drivers of tumour formation. Several cases show heterozygous SMARCB1 mutations without inactivation of the second allele, and we explore this further in vitro. Finding CDKN2A deletions in our discovery cohort, we evaluated CDKN2A protein expression in a tissue microarray. Six out of 16 cases had lost CDKN2A in greater than or equal to 90% of cells, while the remaining cases had retained the protein. Expression analysis of epithelioid sarcoma cell lines by transcriptome sequencing shows a unique profile that does not cluster with any particular tissue type or with other SWI/SNF-aberrant lines. Evaluation of the levels of members of the SWI/SNF complex other than SMARCB1 revealed that these proteins are expressed as part of a residual complex, similarly to previously studied rhabdoid tumour lines. This residual SWI/SNF is susceptible to synthetic lethality and may therefore indicate a therapeutic opportunity.
Subject(s)
Sarcoma/genetics , Transcriptome , Blotting, Western , Cell Line, Tumor , Chromosomal Proteins, Non-Histone/genetics , DNA-Binding Proteins/genetics , Gene Knockdown Techniques , High-Throughput Nucleotide Sequencing , Humans , Immunohistochemistry , Immunoprecipitation , In Situ Hybridization, Fluorescence , Multiplex Polymerase Chain Reaction , Polymerase Chain Reaction , RNA, Small Interfering , SMARCB1 Protein , Tissue Array Analysis , Transcription Factors/genetics , TransfectionABSTRACT
Circulating tumour DNA analysis can be used to track tumour burden and analyse cancer genomes non-invasively but the extent to which it represents metastatic heterogeneity is unknown. Here we follow a patient with metastatic ER-positive and HER2-positive breast cancer receiving two lines of targeted therapy over 3 years. We characterize genomic architecture and infer clonal evolution in eight tumour biopsies and nine plasma samples collected over 1,193 days of clinical follow-up using exome and targeted amplicon sequencing. Mutation levels in the plasma samples reflect the clonal hierarchy inferred from sequencing of tumour biopsies. Serial changes in circulating levels of sub-clonal private mutations correlate with different treatment responses between metastatic sites. This comparison of biopsy and plasma samples in a single patient with metastatic breast cancer shows that circulating tumour DNA can allow real-time sampling of multifocal clonal evolution.
Subject(s)
Brain Neoplasms/genetics , Breast Neoplasms/genetics , Carcinoma, Ductal, Breast/genetics , Clonal Evolution/genetics , DNA, Neoplasm/genetics , Liver Neoplasms/genetics , Lung Neoplasms/genetics , Spinal Neoplasms/genetics , Adult , Antineoplastic Combined Chemotherapy Protocols/therapeutic use , Bayes Theorem , Brain Neoplasms/secondary , Breast Neoplasms/drug therapy , Breast Neoplasms/metabolism , Breast Neoplasms/pathology , Carcinoma, Ductal, Breast/drug therapy , Carcinoma, Ductal, Breast/metabolism , Carcinoma, Ductal, Breast/pathology , Case-Control Studies , Deoxycytidine/administration & dosage , Deoxycytidine/analogs & derivatives , Female , Humans , Lapatinib , Liver Neoplasms/secondary , Lung Neoplasms/secondary , Mutation , Neoplasm Metastasis , Quinazolines/administration & dosage , Receptor, ErbB-2/metabolism , Receptors, Estrogen/metabolism , Sequence Analysis, DNA , Spinal Neoplasms/secondary , Tamoxifen/administration & dosage , Trastuzumab/administration & dosage , GemcitabineABSTRACT
IMPORTANCE: E-cadherin (CDH1) is a cancer predisposition gene mutated in families meeting clinically defined hereditary diffuse gastric cancer (HDGC). Reliable estimates of cancer risk and spectrum in germline mutation carriers are essential for management. For families without CDH1 mutations, genetic-based risk stratification has not been possible, resulting in limited clinical options. OBJECTIVES: To derive accurate estimates of gastric and breast cancer risks in CDH1 mutation carriers and determine if germline mutations in other genes are associated with HDGC. DESIGN, SETTING, AND PARTICIPANTS: Testing for CDH1 germline mutations was performed on 183 index cases meeting clinical criteria for HDGC. Penetrance was derived from 75 mutation-positive families from within this and other cohorts, comprising 3858 probands (353 with gastric cancer and 89 with breast cancer). Germline DNA from 144 HDGC probands lacking CDH1 mutations was screened using multiplexed targeted sequencing for 55 cancer-associated genes. MAIN OUTCOMES AND MEASURES: Accurate estimates of gastric and breast cancer risks in CDH1 mutation carriers and the relative contribution of other cancer predisposition genes in familial gastric cancers. RESULTS: Thirty-one distinct pathogenic CDH1 mutations (14 novel) were identified in 34 of 183 index cases (19%). By the age of 80 years, the cumulative incidence of gastric cancer was 70% (95% CI, 59%-80%) for males and 56% (95% CI, 44%-69%) for females, and the risk of breast cancer for females was 42% (95% CI, 23%-68%). In CDH1 mutation-negative index cases, candidate mutations were identified in 16 of 144 probands (11%), including mutations within genes of high and moderate penetrance: CTNNA1, BRCA2, STK11, SDHB, PRSS1, ATM, MSR1, and PALB2. CONCLUSIONS AND RELEVANCE: This is the largest reported series of CDH1 mutation carriers, providing more precise estimates of age-associated risks of gastric and breast cancer that will improve counseling of unaffected carriers. In HDGC families lacking CDH1 mutations, testing of CTNNA1 and other tumor suppressor genes should be considered. Clinically defined HDGC families can harbor mutations in genes (ie, BRCA2) with different clinical ramifications from CDH1. Therefore, we propose that HDGC syndrome may be best defined by mutations in CDH1 and closely related genes, rather than through clinical criteria that capture families with heterogeneous susceptibility profiles.
Subject(s)
Biomarkers, Tumor/genetics , Breast Neoplasms/genetics , Cadherins/genetics , Germ-Line Mutation , Stomach Neoplasms/genetics , Adult , Age Distribution , Age Factors , Aged , Aged, 80 and over , Antigens, CD , Breast Neoplasms/epidemiology , Breast Neoplasms/pathology , Canada/epidemiology , DNA Mutational Analysis , Europe/epidemiology , Female , Genetic Predisposition to Disease , Heredity , Humans , Incidence , Male , Middle Aged , Pedigree , Penetrance , Phenotype , Predictive Value of Tests , Risk Assessment , Risk Factors , Sex Distribution , Sex Factors , Stomach Neoplasms/epidemiology , Stomach Neoplasms/pathology , Young AdultABSTRACT
OBJECTIVE: We reported association of haplotypes formed by IL-1b (IL1B)-511C/T (rs16944) and a variable number of tandem repeats (rs2234663) in intron 3 of IL-1 receptor antagonist (IL1RN) with rate of lung function decline in smoking-induced COPD. The aim of current study was to further investigate this association. METHODS: We genotyped an additional 19 polymorphisms in IL1 cluster (including IL1A, IL1B and IL1RN) in non-Hispanic whites who had the fastest (n = 268) and the slowest (n = 292) decline of FEV1% predicted in the same study. We also analyzed the association of all 21 polymorphisms with serum CRP levels. RESULTS: None of 21 polymorphisms showed significant association with rate of decline of lung function or CRP levels after adjusting for multiple comparisons. Before adjusting for multiple comparisons, only IL1RN_19327 (rs315949) showed significant association with lung function decline (P = 0.03, additive model). The frequencies of genotypes containing the IL1RN_19327A allele were 71.9% and 62.2%, respectively in the fast and slow decline groups (P = 0.02, odds ratio = 1.6, 95% confidence interval = 1.1-2.3); the IL1B_5200 (rs1143633) and rs2234663 in IL1RN were associated with serum CRP levels (P=0.04 and 0.03, respectively). CONCLUSIONS: No single marker was significantly associated with either rate of lung function decline or serum CRP levels.
Subject(s)
C-Reactive Protein/analysis , Interleukin-1/genetics , Pulmonary Disease, Chronic Obstructive/genetics , Smoking/adverse effects , Adult , Aged , Enzyme-Linked Immunosorbent Assay , Female , Genetic Predisposition to Disease/genetics , Genotype , Humans , Male , Middle Aged , Multigene Family , Polymorphism, Single Nucleotide , Pulmonary Disease, Chronic Obstructive/etiology , Respiratory Function TestsABSTRACT
Human cancers, including breast cancers, comprise clones differing in mutation content. Clones evolve dynamically in space and time following principles of Darwinian evolution, underpinning important emergent features such as drug resistance and metastasis. Human breast cancer xenoengraftment is used as a means of capturing and studying tumour biology, and breast tumour xenografts are generally assumed to be reasonable models of the originating tumours. However, the consequences and reproducibility of engraftment and propagation on the genomic clonal architecture of tumours have not been systematically examined at single-cell resolution. Here we show, using deep-genome and single-cell sequencing methods, the clonal dynamics of initial engraftment and subsequent serial propagation of primary and metastatic human breast cancers in immunodeficient mice. In all 15 cases examined, clonal selection on engraftment was observed in both primary and metastatic breast tumours, varying in degree from extreme selective engraftment of minor (<5% of starting population) clones to moderate, polyclonal engraftment. Furthermore, ongoing clonal dynamics during serial passaging is a feature of tumours experiencing modest initial selection. Through single-cell sequencing, we show that major mutation clusters estimated from tumour population sequencing relate predictably to the most abundant clonal genotypes, even in clonally complex and rapidly evolving cases. Finally, we show that similar clonal expansion patterns can emerge in independent grafts of the same starting tumour population, indicating that genomic aberrations can be reproducible determinants of evolutionary trajectories. Our results show that measurement of genomically defined clonal population dynamics will be highly informative for functional studies using patient-derived breast cancer xenoengraftment.
Subject(s)
Breast Neoplasms/genetics , Breast Neoplasms/pathology , Clone Cells/metabolism , Clone Cells/pathology , Genome, Human/genetics , Single-Cell Analysis , Xenograft Model Antitumor Assays , Animals , Breast Neoplasms/secondary , DNA Mutational Analysis , Genomics , Genotype , High-Throughput Nucleotide Sequencing , Humans , Mice , Neoplasm Transplantation , Time Factors , Transplantation, Heterologous , Xenograft Model Antitumor Assays/methodsABSTRACT
The evolution of cancer genomes within a single tumor creates mixed cell populations with divergent somatic mutational landscapes. Inference of tumor subpopulations has been disproportionately focused on the assessment of somatic point mutations, whereas computational methods targeting evolutionary dynamics of copy number alterations (CNA) and loss of heterozygosity (LOH) in whole-genome sequencing data remain underdeveloped. We present a novel probabilistic model, TITAN, to infer CNA and LOH events while accounting for mixtures of cell populations, thereby estimating the proportion of cells harboring each event. We evaluate TITAN on idealized mixtures, simulating clonal populations from whole-genome sequences taken from genomically heterogeneous ovarian tumor sites collected from the same patient. In addition, we show in 23 whole genomes of breast tumors that the inference of CNA and LOH using TITAN critically informs population structure and the nature of the evolving cancer genome. Finally, we experimentally validated subclonal predictions using fluorescence in situ hybridization (FISH) and single-cell sequencing from an ovarian cancer patient sample, thereby recapitulating the key modeling assumptions of TITAN.
Subject(s)
Algorithms , Computational Biology/methods , DNA Copy Number Variations , Models, Genetic , Neoplasms/genetics , Clone Cells/metabolism , Clone Cells/pathology , Female , Genomics/methods , Genotype , Humans , In Situ Hybridization, Fluorescence/methods , Loss of Heterozygosity , Ovarian Neoplasms/genetics , Polymorphism, Single Nucleotide , Reproducibility of Results , Sequence Analysis, DNA/methods , Triple Negative Breast Neoplasms/geneticsABSTRACT
Bladder cancer (or urothelial cell carcinoma [UCC]) is characterized by field disease (malignant alterations in surrounding mucosa) and frequent recurrences. Whole-genome, exome, and transcriptome sequencing of 38 tumors, including four metachronous tumor pairs and 20 superficial tumors, identified an APOBEC mutational signature in one-third. This was biased toward the sense strand, correlated with mean expression level, and clustered near breakpoints. A>G mutations were up to eight times more frequent on the sense strand (p<0.002) in [ACG]AT contexts. The patient-specific APOBEC signature was negatively correlated to repair-gene expression and was not related to clinicopathological parameters. Mutations in gene families and single genes were related to tumor stage, and expression of chromatin modifiers correlated with survival. Evolutionary and subclonal analyses of early/late tumor pairs showed a unitary origin, and discrete tumor clones contained mutated cancer genes. The ancestral clones contained Pik3ca/Kdm6a mutations and may reflect the field-disease mutations shared among later tumors.
Subject(s)
Carcinoma/genetics , Clonal Evolution , Point Mutation , Urinary Bladder Neoplasms/genetics , APOBEC-1 Deaminase , Carcinoma/pathology , Class I Phosphatidylinositol 3-Kinases , Cytidine Deaminase/genetics , Cytidine Deaminase/metabolism , DNA Repair , DNA, Antisense/genetics , Humans , Phosphatidylinositol 3-Kinases/genetics , Transcriptome , Urinary Bladder Neoplasms/pathologyABSTRACT
Small cell carcinoma of the ovary of hypercalcemic type (SCCOHT) is an extremely rare, aggressive cancer affecting children and young women. We identified germline and somatic inactivating mutations in the SWI/SNF chromatin-remodeling gene SMARCA4 in 75% (9/12) of SCCOHT cases in addition to SMARCA4 protein loss in 82% (14/17) of SCCOHT tumors but in only 0.4% (2/485) of other primary ovarian tumors. These data implicate SMARCA4 in SCCOHT oncogenesis.
Subject(s)
Carcinoma, Small Cell/genetics , DNA Helicases/genetics , Mutation/genetics , Nuclear Proteins/genetics , Ovarian Neoplasms/genetics , Transcription Factors/genetics , Base Sequence , Chromatin Assembly and Disassembly/genetics , Chromosome Mapping , Computational Biology , DNA, Complementary/genetics , Electrophoresis, Polyacrylamide Gel , Exome/genetics , Female , Gene Library , Humans , Immunohistochemistry , Molecular Sequence Data , Sequence Analysis, DNAABSTRACT
Alveolar macrophages play an important role in chronic obstructive pulmonary disease via production of matrix metalloproteinases (MMPs) and cathepsins as well as their inhibitors, tissue inhibitors of metalloproteinases and cystatin C. We hypothesised that expression levels of these molecules by alveolar macrophages at baseline and after stimulation would be influenced by genotype and associated with chronic obstructive pulmonary disease phenotypes. Quantitative PCR and ELISAs/gelatine zymography were used to investigate expression levels of mRNA and protein, respectively. The relationships of expression with genotype, pulmonary function and emphysema were analysed. The results showed that basal expression level of MMP12 mRNA was inversely related to the diffusing capacity of the lung for carbon monoxide/alveolar volume and to forced expiratory volume in 1 s/forced vital capacity after correction for multiple comparisons. The expression level of MMP12 protein stimulated with lipopolysaccharide was also inversely related to the diffusing capacity of the lung for carbon monoxide/alveolar volume and was positively related to the extent of emphysema. The basal expression of MMP1 mRNA was positively correlated with the extent of emphysema. Cathepsin L protein level was positively associated with forced expiratory volume in 1 s % predicted. We conclude that increased MMP12 and MMP1 expression may play a role in the pathogenesis of emphysema. Cathepsin L and MMP9 may be involved in the development of airflow limitation.
Subject(s)
Cathepsin L/genetics , Lung/metabolism , Macrophages, Alveolar/metabolism , Matrix Metalloproteinase 12/genetics , Matrix Metalloproteinase 1/genetics , Matrix Metalloproteinase 9/genetics , Pulmonary Emphysema/genetics , RNA, Messenger/analysis , Aged , Cathepsin L/metabolism , Enzyme-Linked Immunosorbent Assay , Female , Forced Expiratory Volume , Gene Expression Profiling , Humans , Lung/enzymology , Macrophages, Alveolar/enzymology , Male , Matrix Metalloproteinase 1/metabolism , Matrix Metalloproteinase 12/metabolism , Matrix Metalloproteinase 9/metabolism , Middle Aged , Pulmonary Diffusing Capacity , Pulmonary Emphysema/enzymology , Pulmonary Emphysema/physiopathology , Real-Time Polymerase Chain Reaction , Severity of Illness Index , Vital CapacityABSTRACT
[This corrects the article on p. e72162 in vol. 8.].
ABSTRACT
BACKGROUND: OVARIAN CARCINOMAS CONSIST OF AT LEAST FIVE DISTINCT DISEASES: high-grade serous, low-grade serous, clear cell, endometrioid, and mucinous. Biomarker and molecular characterization may represent a more biologically relevant basis for grouping and treating this family of tumors, rather than site of origin. Molecular characteristics have become the new standard for clinical pathology, however development of tailored type-specific therapies is hampered by a failure of basic research to recognize that model systems used to study these diseases must also be stratified. Unrelated model systems do offer value for study of biochemical processes but specific cellular context needs to be applied to assess relevant therapeutic strategies. METHODS: We have focused on the identification of clear cell carcinoma cell line models. A panel of 32 "ovarian cancer" cell lines has been classified into histotypes using a combination of mutation profiles, IHC mutation-surrogates, and a validated immunohistochemical model. All cell lines were identity verified using STR analysis. RESULTS: Many described ovarian clear cell lines have characteristic mutations (including ARID1A and PIK3CA) and an overall molecular/immuno-profile typical of primary tumors. Mutations in TP53 were present in the majority of high-grade serous cell lines. Advanced genomic analysis of bona-fide clear cell carcinoma cell lines also support copy number changes in typical biomarkers such at MET and HNF1B and a lack of any recurrent expressed re-arrangements. CONCLUSIONS: As with primary ovarian tumors, mutation status of cancer genes like ARID1A and TP53 and a general immuno-profile serve well for establishing histotype of ovarian cancer cell We describe specific biomarkers and molecular features to re-classify generic "ovarian carcinoma" cell lines into type specific categories. Our data supports the use of prototype clear cell lines, such as TOV21G and JHOC-5, and questions the use of SKOV3 and A2780 as models of high-grade serous carcinoma.
Subject(s)
Ovarian Neoplasms/genetics , Carcinoma, Ovarian Epithelial , Cell Line, Tumor , Cystadenocarcinoma, Serous/genetics , Cystadenocarcinoma, Serous/pathology , Female , Humans , Mutation , Neoplasms, Glandular and Epithelial/genetics , Neoplasms, Glandular and Epithelial/pathology , Ovarian Neoplasms/pathology , Tumor Suppressor Protein p53/geneticsABSTRACT
High-grade serous ovarian cancer (HGSC) is characterized by poor outcome, often attributed to the emergence of treatment-resistant subclones. We sought to measure the degree of genomic diversity within primary, untreated HGSCs to examine the natural state of tumour evolution prior to therapy. We performed exome sequencing, copy number analysis, targeted amplicon deep sequencing and gene expression profiling on 31 spatially and temporally separated HGSC tumour specimens (six patients), including ovarian masses, distant metastases and fallopian tube lesions. We found widespread intratumoural variation in mutation, copy number and gene expression profiles, with key driver alterations in genes present in only a subset of samples (eg PIK3CA, CTNNB1, NF1). On average, only 51.5% of mutations were present in every sample of a given case (range 10.2-91.4%), with TP53 as the only somatic mutation consistently present in all samples. Complex segmental aneuploidies, such as whole-genome doubling, were present in a subset of samples from the same individual, with divergent copy number changes segregating independently of point mutation acquisition. Reconstruction of evolutionary histories showed one patient with mixed HGSC and endometrioid histology, with common aetiologic origin in the fallopian tube and subsequent selection of different driver mutations in the histologically distinct samples. In this patient, we observed mixed cell populations in the early fallopian tube lesion, indicating that diversity arises at early stages of tumourigenesis. Our results revealed that HGSCs exhibit highly individual evolutionary trajectories and diverse genomic tapestries prior to therapy, exposing an essential biological characteristic to inform future design of personalized therapeutic solutions and investigation of drug-resistance mechanisms.
Subject(s)
Cystadenocarcinoma, Serous/genetics , DNA Mutational Analysis/methods , Gene Expression Regulation, Neoplastic , Genetic Variation/genetics , Ovarian Neoplasms/genetics , Aged , Clone Cells , Cystadenocarcinoma, Serous/secondary , Disease Progression , Drug Resistance , Female , Gene Dosage , Gene Expression Profiling , Humans , Middle Aged , Neoplasm Staging , Ovarian Neoplasms/pathology , Real-Time Polymerase Chain ReactionABSTRACT
Loss of heterozygosity (LOH) and copy number alteration (CNA) feature prominently in the somatic genomic landscape of tumors. As such, karyotypic aberrations in cancer genomes have been studied extensively to discover novel oncogenes and tumor-suppressor genes. Advances in sequencing technology have enabled the cost-effective detection of tumor genome and transcriptome mutation events at single-base-pair resolution; however, computational methods for predicting segmental regions of LOH in this context are not yet fully explored. Consequently, whole transcriptome, nucleotide-level resolution analysis of monoallelic expression patterns associated with LOH has not yet been undertaken in cancer. We developed a novel approach for inference of LOH from paired tumor/normal sequence data and applied it to a cohort of 23 triple-negative breast cancer (TNBC) genomes. Following extensive benchmarking experiments, we describe the nucleotide-resolution landscape of LOH in TNBC and assess the consequent effect of LOH on the transcriptomes of these tumors using RNA-seq-derived measurements of allele-specific expression. We show that the majority of monoallelic expression in the transcriptomes of triple-negative breast cancer can be explained by genomic regions of LOH and establish an upper bound for monoallelic expression that may be explained by other tumor-specific modifications such as epigenetics or mutations. Monoallelically expressed genes associated with LOH reveal that cell cycle, homologous recombination and actin-cytoskeletal functions are putatively disrupted by LOH in TNBC. Finally, we show how inference of LOH can be used to interpret allele frequencies of somatic mutations and postulate on temporal ordering of mutations in the evolutionary history of these tumors.