ABSTRACT
Large panels of comprehensively characterized human cancer models, including the Cancer Cell Line Encyclopedia (CCLE), have provided a rigorous framework with which to study genetic variants, candidate targets, and small-molecule and biological therapeutics and to identify new marker-driven cancer dependencies. To improve our understanding of the molecular features that contribute to cancer phenotypes, including drug responses, here we have expanded the characterizations of cancer cell lines to include genetic, RNA splicing, DNA methylation, histone H3 modification, microRNA expression and reverse-phase protein array data for 1,072 cell lines from individuals of various lineages and ethnicities. Integration of these data with functional characterizations such as drug-sensitivity, short hairpin RNA knockdown and CRISPR-Cas9 knockout data reveals potential targets for cancer drugs and associated biomarkers. Together, this dataset and an accompanying public data portal provide a resource for the acceleration of cancer research using model cancer cell lines.
Subject(s)
Cell Line, Tumor , Neoplasms/genetics , Neoplasms/pathology , Antineoplastic Agents/pharmacology , Biomarkers, Tumor , DNA Methylation , Drug Resistance, Neoplasm , Ethnicity/genetics , Gene Editing , Histones/metabolism , Humans , MicroRNAs/genetics , Molecular Targeted Therapy , Neoplasms/metabolism , Protein Array Analysis , RNA SplicingABSTRACT
BACKGROUND: Somatic alterations in the cancer genome, some of which are associated with changes in gene expression, have been characterized in multiple studies across diverse cancer types. However, less is known about germline variants that influence tumor biology by shaping the cancer transcriptome. METHODS: We performed expression quantitative trait loci (eQTL) analyses using multi-dimensional data from The Cancer Genome Atlas to explore the role of germline variation in mediating the cancer transcriptome. After accounting for associations between somatic alterations and gene expression, we determined the contribution of inherited variants to the cancer transcriptome relative to that of somatic variants. Finally, we performed an interaction analysis using estimates of tumor cellularity to identify cell type-restricted eQTLs. RESULTS: The proportion of genes with at least one eQTL varied between cancer types, ranging between 0.8% in melanoma to 28.5% in thyroid cancer and was correlated more strongly with intratumor heterogeneity than with somatic alteration rates. Although contributions to variance in gene expression was low for most genes, some eQTLs accounted for more than 30% of expression of proximal genes. We identified cell type-restricted eQTLs in genes known to be cancer drivers including LPP and EZH2 that were associated with disease-specific mortality in TCGA but not associated with disease risk in published GWAS. Together, our results highlight the need to consider germline variation in interpreting cancer biology beyond risk prediction.
Subject(s)
Genome-Wide Association Study , Melanoma , Genetic Predisposition to Disease , Genome-Wide Association Study/methods , Humans , Polymorphism, Single Nucleotide , TranscriptomeABSTRACT
Cell autonomous cancer dependencies are now routinely identified using CRISPR loss-of-function viability screens. However, a bias exists that makes it difficult to assess the true essentiality of genes located in amplicons, since the entire amplified region can exhibit lethal scores. These false-positive hits can either be discarded from further analysis, which in cancer models can represent a significant number of hits, or methods can be developed to rescue the true-positives within amplified regions. We propose two methods to rescue true positive hits in amplified regions by correcting for this copy number artefact. The Local Drop Out (LDO) method uses the relative lethality scores within genomic regions to assess true essentiality and does not require additional orthogonal data (e.g. copy number value). LDO is meant to be used in screens covering a dense region of the genome (e.g. a whole chromosome or the whole genome). The General Additive Model (GAM) method models the screening data as a function of the known copy number values and removes the systematic effect from the measured lethality. GAM does not require the same density as LDO, but does require prior knowledge of the copy number values. Both methods have been developed with single sample experiments in mind so that the correction can be applied even in smaller screens. Here we demonstrate the efficacy of both methods at removing the copy number effect and rescuing hits from some of the amplified regions. We estimate a 70-80% decrease of false positive hits with either method in regions of high copy number compared to no correction.
Subject(s)
Clustered Regularly Interspaced Short Palindromic Repeats , DNA Copy Number Variations/genetics , Neoplasms/genetics , Artifacts , Astrocytoma/genetics , Astrocytoma/pathology , Brain Neoplasms/genetics , Brain Neoplasms/pathology , Cell Line, Tumor , Cell Proliferation , Datasets as Topic , False Positive Reactions , Genomics , Humans , Models, Theoretical , Neoplasms/pathology , Stomach Neoplasms/genetics , Stomach Neoplasms/pathologyABSTRACT
Genomic structural variants (SVs) are abundant in humans, differing from other forms of variation in extent, origin and functional impact. Despite progress in SV characterization, the nucleotide resolution architecture of most SVs remains unknown. We constructed a map of unbalanced SVs (that is, copy number variants) based on whole genome DNA sequencing data from 185 human genomes, integrating evidence from complementary SV discovery approaches with extensive experimental validations. Our map encompassed 22,025 deletions and 6,000 additional SVs, including insertions and tandem duplications. Most SVs (53%) were mapped to nucleotide resolution, which facilitated analysing their origin and functional impact. We examined numerous whole and partial gene deletions with a genotyping approach and observed a depletion of gene disruptions amongst high frequency deletions. Furthermore, we observed differences in the size spectra of SVs originating from distinct formation mechanisms, and constructed a map of SV hotspots formed by common mechanisms. Our analytical framework and SV map serves as a resource for sequencing-based association studies.
Subject(s)
DNA Copy Number Variations/genetics , Genetics, Population , Genome, Human/genetics , Genomics , Gene Duplication/genetics , Genetic Predisposition to Disease/genetics , Genotype , Humans , Mutagenesis, Insertional/genetics , Reproducibility of Results , Sequence Analysis, DNA , Sequence Deletion/geneticsABSTRACT
Despite great progress in identifying genetic variants that influence human disease, most inherited risk remains unexplained. A more complete understanding requires genome-wide studies that fully examine less common alleles in populations with a wide range of ancestry. To inform the design and interpretation of such studies, we genotyped 1.6 million common single nucleotide polymorphisms (SNPs) in 1,184 reference individuals from 11 global populations, and sequenced ten 100-kilobase regions in 692 of these individuals. This integrated data set of common and rare alleles, called 'HapMap 3', includes both SNPs and copy number polymorphisms (CNPs). We characterized population-specific differences among low-frequency variants, measured the improvement in imputation accuracy afforded by the larger reference panel, especially in imputing SNPs with a minor allele frequency of Subject(s)
DNA Copy Number Variations
, Genome, Human
, Polymorphism, Single Nucleotide
, Population Groups/genetics
, Human Genome Project
, Humans
ABSTRACT
Genetic variation among individual humans occurs on many different scales, ranging from gross alterations in the human karyotype to single nucleotide changes. Here we explore variation on an intermediate scale--particularly insertions, deletions and inversions affecting from a few thousand to a few million base pairs. We employed a clone-based method to interrogate this intermediate structural variation in eight individuals of diverse geographic ancestry. Our analysis provides a comprehensive overview of the normal pattern of structural variation present in these genomes, refining the location of 1,695 structural variants. We find that 50% were seen in more than one individual and that nearly half lay outside regions of the genome previously described as structurally variant. We discover 525 new insertion sequences that are not present in the human reference genome and show that many of these are variable in copy number between individuals. Complete sequencing of 261 structural variants reveals considerable locus complexity and provides insights into the different mutational processes that have shaped the human genome. These data provide the first high-resolution sequence map of human structural variation--a standard for genotyping platforms and a prelude to future individual genome sequencing projects.
Subject(s)
Genetic Variation/genetics , Genome, Human/genetics , Physical Chromosome Mapping , Sequence Analysis, DNA , Chromosome Inversion/genetics , Euchromatin/genetics , Gene Deletion , Geography , Haplotypes , Humans , Mutagenesis, Insertional/genetics , Polymorphism, Single Nucleotide/genetics , Racial Groups/genetics , Reproducibility of ResultsABSTRACT
Mutations in the genes encoding isocitrate dehydrogenase 1 and 2 (IDH1/2) occur in a variety of tumor types, resulting in production of the proposed oncometabolite, 2-hydroxyglutarate (2-HG). How mutant IDH and 2-HG alter signaling pathways to promote cancer, however, remains unclear. Additionally, there exist relatively few cell lines with IDH mutations. To examine the effect of endogenous IDH mutations and 2-HG, we created a panel of isogenic epithelial cell lines with either wild-type IDH1/2 or clinically relevant IDH1/2 mutations. Differences were noted in the ability of IDH mutations to cause robust 2-HG accumulation. IDH1/2 mutants that produce high levels of 2-HG cause an epithelial-mesenchymal transition (EMT)-like phenotype, characterized by changes in EMT-related gene expression and cellular morphology. 2-HG is sufficient to recapitulate aspects of this phenotype in the absence of an IDH mutation. In the cells types examined, mutant IDH-induced EMT is dependent on up-regulation of the transcription factor ZEB1 and down-regulation of the miR-200 family of microRNAs. Furthermore, sustained knockdown of IDH1 in IDH1 R132H mutant cells is sufficient to reverse many characteristics of EMT, demonstrating that continued expression of mutant IDH is required to maintain this phenotype. These results suggest mutant IDH proteins can reversibly deregulate discrete signaling pathways that contribute to tumorigenesis.
Subject(s)
Epithelial-Mesenchymal Transition , Gene Expression Regulation, Enzymologic , Gene Expression Regulation, Neoplastic , Homeodomain Proteins/metabolism , Isocitrate Dehydrogenase/biosynthesis , MicroRNAs/biosynthesis , Mutation, Missense , Neoplasm Proteins/metabolism , Neoplasms/metabolism , RNA, Neoplasm/biosynthesis , Transcription Factors/metabolism , Amino Acid Substitution , Cell Line, Tumor , Glutarates/metabolism , Homeodomain Proteins/genetics , Humans , Isocitrate Dehydrogenase/genetics , MicroRNAs/genetics , Neoplasm Proteins/genetics , Neoplasms/genetics , Neoplasms/pathology , RNA, Neoplasm/genetics , Transcription Factors/genetics , Up-Regulation/genetics , Zinc Finger E-box-Binding Homeobox 1ABSTRACT
High coverage whole genome sequencing provides near complete information about genetic variation. However, other technologies can be more efficient in some settings by (a) reducing redundant coverage within samples and (b) exploiting patterns of genetic variation across samples. To characterize as many samples as possible, many genetic studies therefore employ lower coverage sequencing or SNP array genotyping coupled to statistical imputation. To compare these approaches individually and in conjunction, we developed a statistical framework to estimate genotypes jointly from sequence reads, array intensities, and imputation. In European samples, we find similar sensitivity (89%) and specificity (99.6%) from imputation with either 1× sequencing or 1 M SNP arrays. Sensitivity is increased, particularly for low-frequency polymorphisms (MAF < 5%), when low coverage sequence reads are added to dense genome-wide SNP arrays--the converse, however, is not true. At sites where sequence reads and array intensities produce different sample genotypes, joint analysis reduces genotype errors and identifies novel error modes. Our joint framework informs the use of next-generation sequencing in genome wide association studies and supports development of improved methods for genotype calling.
Subject(s)
Genomics/methods , Oligonucleotide Array Sequence Analysis/methods , Polymorphism, Single Nucleotide , Sequence Analysis, DNA/methods , Algorithms , Cluster Analysis , Databases, Genetic , Genome-Wide Association Study , Genotype , Humans , Sensitivity and Specificity , White PeopleABSTRACT
Investigators have linked rare copy number variation (CNVs) to neuropsychiatric diseases, such as schizophrenia. One hypothesis is that CNV events cause disease by affecting genes with specific brain functions. Under these circumstances, we expect that CNV events in cases should impact brain-function genes more frequently than those events in controls. Previous publications have applied "pathway" analyses to genes within neuropsychiatric case CNVs to show enrichment for brain-functions. While such analyses have been suggestive, they often have not rigorously compared the rates of CNVs impacting genes with brain function in cases to controls, and therefore do not address important confounders such as the large size of brain genes and overall differences in rates and sizes of CNVs. To demonstrate the potential impact of confounders, we genotyped rare CNV events in 2,415 unaffected controls with Affymetrix 6.0; we then applied standard pathway analyses using four sets of brain-function genes and observed an apparently highly significant enrichment for each set. The enrichment is simply driven by the large size of brain-function genes. Instead, we propose a case-control statistical test, cnv-enrichment-test, to compare the rate of CNVs impacting specific gene sets in cases versus controls. With simulations, we demonstrate that cnv-enrichment-test is robust to case-control differences in CNV size, CNV rate, and systematic differences in gene size. Finally, we apply cnv-enrichment-test to rare CNV events published by the International Schizophrenia Consortium (ISC). This approach reveals nominal evidence of case-association in neuronal-activity and the learning gene sets, but not the other two examined gene sets. The neuronal-activity genes have been associated in a separate set of schizophrenia cases and controls; however, testing in independent samples is necessary to definitively confirm this association. Our method is implemented in the PLINK software package.
Subject(s)
Brain/physiopathology , DNA Copy Number Variations/genetics , Genes/genetics , Genetic Predisposition to Disease , Schizophrenia/genetics , Schizophrenia/physiopathology , Case-Control Studies , Computer Simulation , Databases, Genetic , Gene Deletion , Genome, Human/genetics , Humans , Models, Genetic , Risk AssessmentABSTRACT
As we move forward from the current generation of genome-wide association (GWA) studies, additional cohorts of different ancestries will be studied to increase power, fine map association signals, and generalize association results to additional populations. Knowledge of genetic ancestry as well as population substructure will become increasingly important for GWA studies in populations of unknown ancestry. Here we propose genotyping pooled DNA samples using genome-wide SNP arrays as a viable option to efficiently and inexpensively estimate admixture proportion and identify ancestry informative markers (AIMs) in populations of unknown origin. We constructed DNA pools from African American, Native Hawaiian, Latina, and Jamaican samples and genotyped them using the Affymetrix 6.0 array. Aided by individual genotype data from the African American cohort, we established quality control filters to remove poorly performing SNPs and estimated allele frequencies for the remaining SNPs in each panel. We then applied a regression-based method to estimate the proportion of admixture in each cohort using the allele frequencies estimated from pooling and populations from the International HapMap Consortium as reference panels, and identified AIMs unique to each population. In this study, we demonstrated that genotyping pooled DNA samples yields estimates of admixture proportion that are both consistent with our knowledge of population history and similar to those obtained by genotyping known AIMs. Furthermore, through validation by individual genotyping, we demonstrated that pooling is quite effective for identifying SNPs with large allele frequency differences (i.e., AIMs) and that these AIMs are able to differentiate two closely related populations (HapMap JPT and CHB).
Subject(s)
Gene Pool , Genetics, Population/methods , Genome, Human/genetics , Phylogeny , Asian People/genetics , Gene Frequency/genetics , Genetic Markers , Genotype , Humans , Principal Component Analysis , Quality Control , Reproducibility of ResultsABSTRACT
The considerable uncertainty regarding cancer risks associated with inherited mutations of BRCA2 is due to unknown factors. To investigate whether common genetic variants modify penetrance for BRCA2 mutation carriers, we undertook a two-staged genome-wide association study in BRCA2 mutation carriers. In stage 1 using the Affymetrix 6.0 platform, 592,163 filtered SNPs genotyped were available on 899 young (<40 years) affected and 804 unaffected carriers of European ancestry. Associations were evaluated using a survival-based score test adjusted for familial correlations and stratified by country of the study and BRCA2*6174delT mutation status. The genomic inflation factor (λ) was 1.011. The stage 1 association analysis revealed multiple variants associated with breast cancer risk: 3 SNPs had p-values<10(-5) and 39 SNPs had p-values<10(-4). These variants included several previously associated with sporadic breast cancer risk and two novel loci on chromosome 20 (rs311499) and chromosome 10 (rs16917302). The chromosome 10 locus was in ZNF365, which contains another variant that has recently been associated with breast cancer in an independent study of unselected cases. In stage 2, the top 85 loci from stage 1 were genotyped in 1,264 cases and 1,222 controls. Hazard ratios (HR) and 95% confidence intervals (CI) for stage 1 and 2 were combined and estimated using a retrospective likelihood approach, stratified by country of residence and the most common mutation, BRCA2*6174delT. The combined per allele HR of the minor allele for the novel loci rs16917302 was 0.75 (95% CI 0.66-0.86, ) and for rs311499 was 0.72 (95% CI 0.61-0.85, ). FGFR2 rs2981575 had the strongest association with breast cancer risk (per allele HR = 1.28, 95% CI 1.18-1.39, ). These results indicate that SNPs that modify BRCA2 penetrance identified by an agnostic approach thus far are limited to variants that also modify risk of sporadic BRCA2 wild-type breast cancer.
Subject(s)
BRCA2 Protein/genetics , Breast Neoplasms/genetics , Genome-Wide Association Study/methods , Polymorphism, Single Nucleotide , Adult , Chromosomes, Human, Pair 10 , Chromosomes, Human, Pair 20 , DNA-Binding Proteins/genetics , Female , Gene Frequency , Genetic Predisposition to Disease/genetics , Haplotypes , Heterozygote , Humans , Linkage Disequilibrium , Middle Aged , Mutation , Penetrance , Receptor, Fibroblast Growth Factor, Type 2/genetics , Risk Factors , Transcription Factors/genetics , White People/geneticsABSTRACT
Although androgen deprivation treatment often effectively decreases prostate cancer, incurable metastatic castration-resistant prostate cancer (CRPC) eventually occurs. It is important to understand how CRPC metastasis progresses, which is not clearly defined. The loss of PTEN, a phosphatase to dephosphorylate phosphatidylinositol 3,4,5-trisphosphate in the PI3K pathway, occurs in up to 70% to 80% of CRPC. We generated a mouse androgen-independent prostate cancer cell line (PKO) from PTEN null and Hi-Myc transgenic mice in C57BL/6 background. We confirmed that this PKO cell line has an activated PI3K pathway and can metastasize into the femur and tibia of immunodeficient nude and immunocompetent C57BL/6 mice. In vitro, we found that androgen deprivation significantly enhanced PKO cell migration/invasion via the p110ß isoform-depended PAK1-MAPK activation. Inhibition of the p110ß-PAK1 axis significantly decreased prostate cancer cell migration/invasion. Of note, our analysis using clinical samples showed that PAK1 is more activated in CRPC than in advanced prostate cancer; high PAK1/phosphorylated-PAK1 levels are associated with decreased survival rates in patients with CRPC. All the information suggests that this cell line reflects the characteristics of CRPC cells and can be applied to dissect the mechanism of CRPC initiation and progression. This study also shows that PAK1 is a potential target for CRPC treatment. IMPLICATIONS: This study uses a newly generated PTEN null prostate cancer cell line to define a critical functional role of p110ß-PAK1 in CRPC migration/invasion. This study also shows that the p110ß-PAK1 axis can potentially be a therapeutic target in CRPC metastasis.
Subject(s)
Prostatic Neoplasms, Castration-Resistant , Animals , Humans , Male , Mice , Androgen Antagonists , Androgens/therapeutic use , Cell Line, Tumor , Mice, Inbred C57BL , Mice, Transgenic , p21-Activated Kinases/genetics , Phosphatidylinositol 3-Kinases/metabolism , Prostatic Neoplasms, Castration-Resistant/metabolism , PTEN Phosphohydrolase/genetics , PTEN Phosphohydrolase/metabolism , Receptors, Androgen/metabolismABSTRACT
For a majority of patients with non-small cell lung cancer with EGFR mutations, treatment with EGFR inhibitors (EGFRi) induces a clinical response. Despite this initial reduction in tumor size, residual disease persists that leads to disease relapse. Elucidating the preexisting biological differences between sensitive cells and surviving drug-tolerant persister cells and deciphering how drug-tolerant cells evolve in response to treatment could help identify strategies to improve the efficacy of EGFRi. In this study, we tracked the origins and clonal evolution of drug-tolerant cells at a high resolution by using an expressed barcoding system coupled with single-cell RNA sequencing. This platform enabled longitudinal profiling of gene expression and drug sensitivity in response to EGFRi across a large number of clones. Drug-tolerant cells had higher expression of key survival pathways such as YAP and EMT at baseline and could also differentially adapt their gene expression following EGFRi treatment compared with sensitive cells. In addition, drug combinations targeting common downstream components (MAPK) or orthogonal factors (chemotherapy) showed greater efficacy than EGFRi alone, which is attributable to broader targeting of the heterogeneous EGFRi-tolerance mechanisms present in tumors. Overall, this approach facilitates thorough examination of clonal evolution in response to therapy that could inform the development of improved diagnostic approaches and treatment strategies for targeting drug-tolerant cells. SIGNIFICANCE: The evolution and heterogeneity of EGFR inhibitor tolerance are identified in a large number of clones at enhanced cellular and temporal resolution using an expressed barcode technology coupled with single-cell RNA sequencing.
Subject(s)
Carcinoma, Non-Small-Cell Lung , Lung Neoplasms , Humans , Carcinoma, Non-Small-Cell Lung/drug therapy , Carcinoma, Non-Small-Cell Lung/genetics , Carcinoma, Non-Small-Cell Lung/pathology , Lung Neoplasms/drug therapy , Lung Neoplasms/genetics , Lung Neoplasms/pathology , ErbB Receptors/genetics , ErbB Receptors/metabolism , Neoplasm Recurrence, Local , Drug ToleranceABSTRACT
A common outcome of androgen deprivation in prostate cancer therapy is disease relapse and progression to castration-resistant prostate cancer (CRPC) via multiple mechanisms. To gain insight into the recent clinical findings that highlighted genomic alterations leading to hyperactivation of PI3K, we examined the roles of the commonly expressed p110 catalytic isoforms of PI3K in a murine model of Pten-null invasive CRPC. While blocking p110α had negligible effects in the development of Pten-null invasive CRPC, either genetic or pharmacologic perturbation of p110ß dramatically slowed CRPC initiation and progression. Once fully established, CRPC tumors became partially resistant to p110ß inhibition, indicating the acquisition of new dependencies. Driven by our genomic analyses highlighting potential roles for the p110ß/RAC/PAK1 and ß-catenin pathways in CRPC, we found that combining p110ß with RAC/PAK1 or tankyrase inhibitors significantly reduced the growth of murine and human CRPC organoids in vitro and in vivo. Because p110ß activity is dispensable for most physiologic processes, our studies support novel therapeutic strategies both for preventing disease progression into CRPC and for treating CRPC. IMPLICATIONS: This work establishes p110ß as a promising target for preventing the progression of primary PTEN-deficient prostate tumors to CRPC, and for treating established CRPC in combination with RAC/PAK1 or tankyrase inhibitors.
Subject(s)
Prostatic Neoplasms, Castration-Resistant , Tankyrases , Androgen Antagonists , Animals , Humans , Male , Mice , PTEN Phosphohydrolase/genetics , Phosphatidylinositol 3-Kinases , Prostate , Prostatic Neoplasms, Castration-Resistant/drug therapy , Prostatic Neoplasms, Castration-Resistant/geneticsABSTRACT
Genome-wide genotyping of a cohort using pools rather than individual samples has long been proposed as a cost-saving alternative for performing genome-wide association (GWA) studies. However, successful disease gene mapping using pooled genotyping has thus far been limited to detecting common variants with large effect sizes, which tend not to exist for many complex common diseases or traits. Therefore, for DNA pooling to be a viable strategy for conducting GWA studies, it is important to determine whether commonly used genome-wide SNP array platforms such as the Affymetrix 6.0 array can reliably detect common variants of small effect sizes using pooled DNA. Taking obesity and age at menarche as examples of human complex traits, we assessed the feasibility of genome-wide genotyping of pooled DNA as a single-stage design for phenotype association. By individually genotyping the top associations identified by pooling, we obtained a 14- to 16-fold enrichment of SNPs nominally associated with the phenotype, but we likely missed the top true associations. In addition, we assessed whether genotyping pooled DNA can serve as an inexpensive screen as the second stage of a multi-stage design with a large number of samples by comparing the most cost-effective 3-stage designs with 80% power to detect common variants with genotypic relative risk of 1.1, with and without pooling. Given the current state of the specific technology we employed and the associated genotyping costs, we showed through simulation that a design involving pooling would be 1.07 times more expensive than a design without pooling. Thus, while a significant amount of information exists within the data from pooled DNA, our analysis does not support genotyping pooled DNA as a means to efficiently identify common variants contributing small effects to phenotypes of interest. While our conclusions were based on the specific technology and study design we employed, the approach presented here will be useful for evaluating the utility of other or future genome-wide genotyping platforms in pooled DNA studies.
Subject(s)
Genome-Wide Association Study/methods , Oligonucleotide Array Sequence Analysis/methods , Sequence Analysis, DNA/methods , Adolescent , Child , Cohort Studies , Computer Simulation , Female , Genetic Variation , Genome-Wide Association Study/economics , Humans , Male , Menarche/genetics , Obesity/genetics , Oligonucleotide Array Sequence Analysis/economics , Polymorphism, Single NucleotideABSTRACT
Three founder mutations in BRCA1 and BRCA2 contribute to the risk of hereditary breast and ovarian cancer in Ashkenazi Jews (AJ). They are observed at increased frequency in the AJ compared to other BRCA mutations in Caucasian non-Jews (CNJ). Several authors have proposed that elevated allele frequencies in the surrounding genomic regions reflect adaptive or balancing selection. Such proposals predict long-range linkage disequilibrium (LD) resulting from a selective sweep, although genetic drift in a founder population may also act to create long-distance LD. To date, few studies have used the tools of statistical genomics to examine the likelihood of long-range LD at a deleterious locus in a population that faced a genetic bottleneck. We studied the genotypes of hundreds of women from a large international consortium of BRCA1 and BRCA2 mutation carriers and found that AJ women exhibited long-range haplotypes compared to CNJ women. More than 50% of the AJ chromosomes with the BRCA1 185delAG mutation share an identical 2.1 Mb haplotype and nearly 16% of AJ chromosomes carrying the BRCA2 6174delT mutation share a 1.4 Mb haplotype. Simulations based on the best inference of Ashkenazi population demography indicate that long-range haplotypes are expected in the context of a genome-wide survey. Our results are consistent with the hypothesis that a local bottleneck effect from population size constriction events could by chance have resulted in the large haplotype blocks observed at high frequency in the BRCA1 and BRCA2 regions of Ashkenazi Jews.
Subject(s)
Arthritis/genetics , BRCA1 Protein/genetics , BRCA2 Protein/genetics , Deafness/genetics , Haplotypes/genetics , Polychondritis, Relapsing/genetics , Base Sequence , Computer Simulation , Female , Founder Effect , Genotype , Heterozygote , Humans , Jews/genetics , Sequence DeletionABSTRACT
BACKGROUND: Autism spectrum disorder is a heritable developmental disorder in which chromosomal abnormalities are thought to play a role. METHODS: As a first component of a genomewide association study of families from the Autism Genetic Resource Exchange (AGRE), we used two novel algorithms to search for recurrent copy-number variations in genotype data from 751 multiplex families with autism. Specific recurrent de novo events were further evaluated in clinical-testing data from Children's Hospital Boston and in a large population study in Iceland. RESULTS: Among the AGRE families, we observed five instances of a de novo deletion of 593 kb on chromosome 16p11.2. Using comparative genomic hybridization, we observed the identical deletion in 5 of 512 children referred to Children's Hospital Boston for developmental delay, mental retardation, or suspected autism spectrum disorder, as well as in 3 of 299 persons with autism in an Icelandic population; the deletion was also carried by 2 of 18,834 unscreened Icelandic control subjects. The reciprocal duplication of this region occurred in 7 affected persons in AGRE families and 4 of the 512 children from Children's Hospital Boston. The duplication also appeared to be a high-penetrance risk factor. CONCLUSIONS: We have identified a novel, recurrent microdeletion and a reciprocal microduplication that carry substantial susceptibility to autism and appear to account for approximately 1% of cases. We did not identify other regions with similar aggregations of large de novo mutations.
Subject(s)
Autistic Disorder/genetics , Chromosome Aberrations , Chromosome Deletion , Chromosomes, Human, Pair 16/genetics , Genetic Predisposition to Disease , Child , Chromosomes, Human, Pair 15/genetics , DNA Mutational Analysis , Developmental Disabilities/genetics , Female , Genotype , Humans , Intellectual Disability/genetics , Male , Phenotype , Sequence Analysis, DNA/methodsABSTRACT
There is a compelling need for new therapeutic strategies for glioblastoma multiforme (GBM). Preclinical target and therapeutic discovery for GBMs is primarily conducted using cell lines grown in serum-containing media, such as U-87 MG, which do not reflect the gene expression profiles of tumors found in GBM patients. To address this lack of representative models, we sought to develop a panel of patient-derived GBM models and characterize their genomic features, using RNA sequencing (RNA-seq) and growth characteristics, both when grown as neurospheres in culture, and grown orthotopically as xenografts in mice. When we compared these with commonly used GBM cell lines in the Cancer Cell Line Encyclopedia (CCLE), we found these patient-derived models to have greater diversity in gene expression and to better correspond to GBMs directly sequenced from patient tumor samples. We also evaluated the potential of these models for targeted therapy, by using the genomic characterization to identify small molecules that inhibit the growth of distinct subsets of GBMs, paving the way for precision medicines for GBM.
ABSTRACT
Advanced ovarian cancers are a leading cause of cancer-related death in women and are currently treated with surgery and chemotherapy. This standard of care is often temporarily successful but exhibits a high rate of relapse, after which, treatment options are few. Here we investigate whether biomarker-guided use of multiple targeted therapies, including small molecules and antibody-drug conjugates, is a viable alternative. A panel of patient-derived ovarian cancer xenografts (PDX), similar in genetics and chemotherapy responsiveness to human tumors, was exposed to 21 monotherapies and combination therapies. Three monotherapies and one combination were found to be active in different subsets of PDX. Analysis of gene expression data identified biomarkers associated with responsiveness to each of the three targeted therapies, none of which directly inhibits an oncogenic driver. While no single treatment had as high a response rate as chemotherapy, nearly 90% of PDXs were eligible for and responded to at least one biomarker-guided treatment, including tumors resistant to standard chemotherapy. The distribution of biomarker positivity in The Cancer Genome Atlas data suggests the potential for a similar precision approach in human patients. SIGNIFICANCE: This study exploits a panel of patient-derived xenografts to demonstrate that most ovarian tumors can be matched to effective biomarker-guided treatments.
Subject(s)
Antineoplastic Combined Chemotherapy Protocols/pharmacology , Biomarkers, Tumor/genetics , Ovarian Neoplasms/drug therapy , Xenograft Model Antitumor Assays/methods , Antineoplastic Agents/pharmacology , Carcinoma, Ovarian Epithelial/drug therapy , Carcinoma, Ovarian Epithelial/genetics , Carcinoma, Ovarian Epithelial/mortality , Carcinoma, Ovarian Epithelial/pathology , Drug Resistance, Neoplasm/drug effects , Drug Resistance, Neoplasm/genetics , Female , Gene Expression Regulation, Neoplastic , Humans , Kaplan-Meier Estimate , Molecular Targeted Therapy/methods , Ovarian Neoplasms/genetics , Ovarian Neoplasms/mortality , Ovarian Neoplasms/pathology , Precision Medicine , Proof of Concept StudyABSTRACT
Interferons (IFNs) are cytokines that play a critical role in limiting infectious and malignant diseases 1-4 . Emerging data suggest that the strength and duration of IFN signaling can differentially impact cancer therapies, including immune checkpoint blockade 5-7 . Here, we characterize the output of IFN signaling, specifically IFN-stimulated gene (ISG) signatures, in primary tumors from The Cancer Genome Atlas. While immune infiltration correlates with the ISG signature in some primary tumors, the existence of ISG signature-positive tumors without evident infiltration of IFN-producing immune cells suggests that cancer cells per se can be a source of IFN production. Consistent with this hypothesis, analysis of patient-derived tumor xenografts propagated in immune-deficient mice shows evidence of ISG-positive tumors that correlates with expression of human type I and III IFNs derived from the cancer cells. Mechanistic studies using cell line models from the Cancer Cell Line Encyclopedia that harbor ISG signatures demonstrate that this is a by-product of a STING-dependent pathway resulting in chronic tumor-derived IFN production. This imposes a transcriptional state on the tumor, poising it to respond to the aberrant accumulation of double-stranded RNA (dsRNA) due to increased sensor levels (MDA5, RIG-I and PKR). By interrogating our functional short-hairpin RNA screen dataset across 398 cancer cell lines, we show that this ISG transcriptional state creates a novel genetic vulnerability. ISG signature-positive cancer cells are sensitive to the loss of ADAR, a dsRNA-editing enzyme that is also an ISG. A genome-wide CRISPR genetic suppressor screen reveals that the entire type I IFN pathway and the dsRNA-activated kinase, PKR, are required for the lethality induced by ADAR depletion. Therefore, tumor-derived IFN resulting in chronic signaling creates a cellular state primed to respond to dsRNA accumulation, rendering ISG-positive tumors susceptible to ADAR loss.