Profound Tissue Specificity in Proliferation Control Underlies Cancer Drivers and Aneuploidy Patterns.

Genomics has provided a detailed structural description of the cancer genome. Identifying oncogenic drivers that work primarily through dosage changes is a current challenge. Unrestrained proliferation is a critical hallmark of cancer. We constructed modular, barcoded libraries of human open reading frames (ORFs) and performed screens for proliferation regulators in multiple cell types. Approximately 10% of genes regulate proliferation, with most performing in an unexpectedly highly tissue-specific manner. Proliferation drivers in a given cell type showed specific enrichment in somatic copy number changes (SCNAs) from cognate tumors and helped predict aneuploidy patterns in those tumors, implying that tissue-type-specific genetic network architectures underlie SCNA and driver selection in different cancers. In vivo screening confirmed these results. We report a substantial contribution to the catalog of SCNA-associated cancer drivers, identifying 147 amplified and 107 deleted genes as potential drivers, and derive insights about the genetic network architecture of aneuploidy in tumors.

Integrated loss- and gain-of-function screens define a core network governing human embryonic stem cell behavior.

Understanding the genetic control of human embryonic stem cell function is foundational for developmental biology and regenerative medicine. Here we describe an integrated genome-scale loss- and gain-of-function screening approach to identify genetic networks governing embryonic stem cell proliferation and differentiation into the three germ layers. We identified a deep link between pluripotency maintenance and survival by showing that genetic alterations that cause pluripotency dissolution simultaneously increase apoptosis resistance. We discovered that the chromatin-modifying complex SAGA and in particular its subunit TADA2B are central regulators of pluripotency, survival, growth, and lineage specification. Joint analysis of all screens revealed that genetic alterations that broadly inhibit differentiation across multiple germ layers drive proliferation and survival under pluripotency-maintaining conditions and coincide with known cancer drivers. Our results show the power of integrated multilayer genetic screening for the robust mapping of complex genetic networks.

Functional genomics reveals that tumors with activating phosphoinositide 3-kinase mutations are dependent on accelerated protein turnover.

Activating mutations in the phosphoinositide 3-kinase (PI3K) signaling pathway are frequently identified in cancer. To identify pathways that support PI3K oncogenesis, we performed a genome-wide RNAi screen in isogenic cell lines harboring wild-type or mutant PIK3CA to search for PI3K synthetic-lethal (SL) genes. A combined analysis of these results with a meta-analysis of two other large-scale RNAi screening data sets in PI3K mutant cancer cell lines converged on ribosomal protein translation and proteasomal protein degradation as critical nononcogene dependencies for PI3K-driven tumors. Genetic or pharmacologic inhibition of either pathway alone, but not together, selectively killed PI3K mutant tumor cells in an mTOR-dependent manner. The expression of ribosomal and proteasomal components was significantly up-regulated in primary human colorectal tumors harboring PI3K pathway activation. Importantly, a PI3K SL gene signature containing the top hits of the SL genes identified in our meta-analysis robustly predicted overall patient survival in colorectal cancer, especially among patients with tumors with an activated PI3K pathway. These results suggest that disruption of protein turnover homeostasis via ribosome or proteasome inhibition may be a novel treatment strategy for PI3K mutant human tumors.

Chromosome evolution screens recapitulate tissue-specific tumor aneuploidy patterns.

Whole chromosome and arm-level copy number alterations occur at high frequencies in tumors, but their selective advantages, if any, are poorly understood. Here, utilizing unbiased whole chromosome genetic screens combined with in vitro evolution to generate arm- and subarm-level events, we iteratively selected the fittest karyotypes from aneuploidized human renal and mammary epithelial cells. Proliferation-based karyotype selection in these epithelial lines modeled tissue-specific tumor aneuploidy patterns in patient cohorts in the absence of driver mutations. Hi-C-based translocation mapping revealed that arm-level events usually emerged in multiples of two via centromeric translocations and occurred more frequently in tetraploids than diploids, contributing to the increased diversity in evolving tetraploid populations. Isogenic clonal lineages enabled elucidation of pro-tumorigenic mechanisms associated with common copy number alterations, revealing Notch signaling potentiation as a driver of 1q gain in breast cancer. We propose that intrinsic, tissue-specific proliferative effects underlie tumor copy number patterns in cancer.

A GATA4-regulated secretory program suppresses tumors through recruitment of cytotoxic CD8 T cells.

The GATA4 transcription factor acts as a master regulator of development of multiple tissues. GATA4 also acts in a distinct capacity to control a stress-inducible pro-inflammatory secretory program that is associated with senescence, a potent tumor suppression mechanism, but also operates in non-senescent contexts such as tumorigenesis. This secretory pathway is composed of chemokines, cytokines, growth factors, and proteases. Since GATA4 is deleted or epigenetically silenced in cancer, here we examine the role of GATA4 in tumorigenesis in mouse models through both loss-of-function and overexpression experiments. We find that GATA4 promotes non-cell autonomous tumor suppression in multiple model systems. Mechanistically, we show that Gata4-dependent tumor suppression requires cytotoxic CD8 T cells and partially requires the secreted chemokine CCL2. Analysis of transcriptome data in human tumors reveals reduced lymphocyte infiltration in GATA4-deficient tumors, consistent with our murine data. Notably, activation of the GATA4-dependent secretory program combined with an anti-PD-1 antibody robustly abrogates tumor growth in vivo.

Lack of association of Klotho gene variants with valvular and vascular calcification in Caucasians: a candidate gene study of the Framingham Offspring Cohort.

BACKGROUND: Valvular and vascular calcification are important early aging phenotypes and represent risk factors for cardiovascular morbidity and mortality. Klotho is a gene primarily expressed in the kidney that has an important role in calcium-phosphate homeostasis. The functional KL-VS variant of Klotho has been associated with aging and cardiovascular disease in human studies, but its role in valvular and vascular calcification remains unknown. We performed a candidate gene study in the Framingham Offspring Cohort to evaluate the effect of KL-VS variant of the Klotho gene on valvular calcification. METHODS: We analyzed the Klotho KL-VS genotype (rs9536314) from the Affymetrix 550K genome-wide dataset, distributed by dbGAP, on 1389 cases and 2139 controls from the Framingham Heart Study Offspring Cohort. Allele and genotype frequencies were compared between cases and controls. Valvular calcification was defined as presence of calcification on the mitral annulus or the aortic valve as determined by echocardiography. A sensitivity analysis of coronary artery calcification by electron beam computed tomography was performed on 1363 patients. RESULTS: The frequency of the TT versus the TG allele was not different between the cases and the controls (39 versus 41%). The KL-VS variant of Klotho was not associated with valvular or vascular calcification, despite adequate power to detect association (86% for odds ratios ≥1.2). In sensitivity analyses, no association (P > 0.001) between other common variants of Klotho, ß-Klotho or fibroblast growth factor-23 and the end points of valvular or vascular calcification was observed. CONCLUSIONS: In our adequately powered candidate gene study, we did not observe an association with the functional KL-VS variant of Klotho and presence of valvular or vascular calcification. Future studies aimed at combining cohorts with echocardiographic phenotypes need to be conducted to identify genetic variants associated with valvular calcification.

Increased FOG-2 in failing myocardium disrupts thyroid hormone-dependent SERCA2 gene transcription.

Reduced expression of sarcoplasmic reticulum calcium ATPase (SERCA)2 and other genes in the adult cardiac gene program has raised consideration of an impaired responsiveness to thyroid hormone (T3) that develops in the advanced failing heart. Here, we show that human and murine cardiomyopathy hearts have increased expression of friend of GATA (FOG)-2, a cardiac nuclear hormone receptor corepressor protein. Cardiac-specific overexpression of FOG-2 in transgenic mice led to depressed cardiac function, activation of the fetal gene program, congestive heart failure, and early death. SERCA2 transcript and protein levels were reduced in FOG-2 transgenic hearts, and FOG-2 overexpression impaired T3-mediated SERCA2 expression in cultured cardiomyocytes. FOG-2 physically interacts with thyroid hormone receptor-alpha1 and abrogated even high levels of T3-mediated SERCA2 promoter activity. These results demonstrate that SERCA2 is an important target of FOG-2 and that increased FOG-2 expression may contribute to a decline in cardiac function in end-stage heart failure by impaired T3 signaling.

Tumor aneuploidy correlates with markers of immune evasion and with reduced response to immunotherapy.

Immunotherapies based on immune checkpoint blockade are highly effective in a subset of patients. An ongoing challenge is the identification of biomarkers that predict which patients will benefit from these therapies. Aneuploidy, also known as somatic copy number alterations (SCNAs), is widespread in cancer and is posited to drive tumorigenesis. Analyzing 12 human cancer types, we find that, for most, highly aneuploid tumors show reduced expression of markers of cytotoxic infiltrating immune cells, especially CD8+ T cells, and increased expression of cell proliferation markers. Different types of SCNAs predict the proliferation and immune signatures, implying distinct underlying mechanisms. Using published data from two clinical trials of immune checkpoint blockade therapy for metastatic melanoma, we found that tumor aneuploidy inversely correlates with patient survival. Together with other tumor characteristics such as tumor mutational load, aneuploidy may thus help identify patients most likely to respond to immunotherapy.

Patterns of resistance and incomplete response to docetaxel by gene expression profiling in breast cancer patients.

PURPOSE: Chemotherapy for operable breast cancer decreases the risk of death. Docetaxel is one of the most active agents in breast cancer, but resistance or incomplete response is frequent. PATIENTS AND METHODS: Core biopsies from 24 patients were obtained before treatment with neoadjuvant docetaxel (four cycles, 100 mg/m(2) every 3 weeks), and response was assessed after chemotherapy. After 3 months of neoadjuvant chemotherapy, surgical specimens (n = 13) were obtained, and laser capture microdissection (LCM; n = 8) was performed to enrich for tumor cells. From each core, surgical, and LCM specimen, sufficient total RNA (3 to 6 microg) was extracted for cDNA array analysis using the Affymetrix HgU95-Av2 GeneChip (Affymetrix, Santa Clara, CA). RESULTS: From the initial core biopsies, differential patterns of expression of 92 genes correlated with docetaxel response (P = .001). However, the molecular patterns of the residual cancers after 3 months of docetaxel treatment were strikingly similar, independent of initial sensitivity or resistance. This relative genetic homogeneity after treatment was observed in both LCM and non-LCM surgical specimens. The residual tumor after treatment in tumors that were initially sensitive indicates selection of a residual and resistant subpopulation of cells. The gene expression pattern was populated by genes involved in cell cycle arrest at G(2)M (eg, mitotic cyclins and cdc2) and survival pathways involving the mammalian target of rapamycin. CONCLUSION: A specific and consistent gene expression pattern was found in residual tumors after docetaxel treatment. These profiles provide therapeutic targets that could lead to improved treatment.

Coordinated upregulation of oxidative pathways and downregulation of lipid biosynthesis underlie obesity resistance in perilipin knockout mice: a microarray gene expression profile.

Obesity is a major risk factor for diabetes and heart disease. We previously reported that the inactivation of the gene for perilipin (plin), an adipocyte lipid droplet surface protein, produced lean and obesity-resistant mice. To dissect the underlying mechanisms involved, we used oligonucleotide microarrays to analyze the gene-expression profile of white adipose tissue (WAT), liver, heart, skeletal muscle, and kidney of plin(-/-) and plin(+/+) mice. As compared with wild-type littermates, the WAT of plin(-/-) mice had 270 and 543 transcripts that were significantly up- or downregulated. There was a coordinated upregulation of genes involved in beta-oxidation, the Krebs cycle, and the electron transport chain concomitant with a downregulation of genes involved in lipid biosynthesis. There was also a significant downregulation of the stearoyl CoA desaturase-1 gene, which has been associated with obesity resistance. Thus, in response to the constitutive activation of lipolysis associated with absence of perilipin, WAT activated pathways to rid itself of the products of lipolysis and activated pathways of energy expenditure that contribute to the observed obesity resistance. The biochemical pathways involved in obesity resistance in plin(-/-) mice identified in this study may represent potential targets for the treatment of obesity.

Gene expression profiling for the prediction of therapeutic response to docetaxel in patients with breast cancer.

BACKGROUND: Systemic chemotherapy for operable breast cancer substantially decreases the risk of death. Patients often have de novo resistance or incomplete response to docetaxel, one of the most active agents in this disease. We postulated that gene expression profiles of the primary breast cancer can predict the response to docetaxel. METHODS: We took core biopsy samples from primary breast tumours in 24 patients before treatment and then assessed tumour response to neoadjuvant docetaxel (four cycles, 100 mg/m2 daily for 3 weeks) by cDNA analysis of RNA extracted from biopsy samples using HgU95-Av2 GeneChip. FINDINGS: From the core biopsy samples, we extracted sufficient total RNA (3-6 microg) for cDNA array analysis using HgU95-Av2 GeneChip. Differential patterns of expression of 92 genes correlated with docetaxel response (p=0.001). Sensitive tumours had higher expression of genes involved in cell cycle, cytoskeleton, adhesion, protein transport, protein modification, transcription, and stress or apoptosis; whereas resistant tumours showed increased expression of some transcriptional and signal transduction genes. In leave-one-out cross-validation analysis, ten of 11 sensitive tumours (90% specificity) and 11 of 13 resistant tumours (85% sensitivity) were correctly classified, with an accuracy of 88%. This 92-gene predictor had positive and negative predictive values of 92% and 83%, respectively. Correlation between RNA expression measured by the arrays and semiquantitative RT-PCR was also ascertained, and our results were validated in an independent set of six patients. INTERPRETATION: If validated, these molecular profiles could allow development of a clinical test for docetaxel sensitivity, thus reducing unnecessary treatment for women with breast cancer.

Notch pathway activation contributes to inhibition of C2C12 myoblast differentiation by ethanol.

The loss of muscle mass in alcoholic myopathy may reflect alcohol inhibition of myogenic cell differentiation into myotubes. Here, using a high content imaging system we show that ethanol inhibits C2C12 myoblast differentiation by reducing myogenic fusion, creating smaller and less complex myotubes compared with controls. Ethanol administration during C2C12 differentiation reduced MyoD and myogenin expression, and microarray analysis identified ethanol activation of the Notch signaling pathway target genes Hes1 and Hey1. A reporter plasmid regulated by the Hes1 proximal promoter was activated by alcohol treatment in C2C12 cells. Treatment of differentiating C2C12 cells with a gamma secretase inhibitor (GSI) abrogated induction of Hes1. On a morphological level GSI treatment completely rescued myogenic fusion defects and partially restored other myotube parameters in response to alcohol. We conclude that alcohol inhibits C2C12 myoblast differentiation and the inhibition of myogenic fusion is mediated by Notch pathway activation.

Formin homology 2 domain containing 3 variants associated with hypertrophic cardiomyopathy.

BACKGROUND: Incomplete penetrance and variable expression of hypertrophic cardiomyopathy (HCM) is well appreciated. Common genetic polymorphisms variants that may affect HCM penetrance and expression have been predicted but are not well established. METHODS AND RESULTS: We performed a case-control genomewide association study to identify common HCM-associated genetic polymorphisms and then asked whether such common variants were more represented in HCM or could explain the heterogeneity of HCM phenotypes. We identified an intronic FHOD3 variant (rs516514) associated with HCM (odds ratio, 2.45; 95% confidence interval, 1.76-3.41; P=1.25×10(-7)) and validated this finding in an independent cohort. Next, we tested FHOD3-V1151I (rs2303510), a nonsynonymous variant in partial linkage disequilibrium with rs516514, and we detected an even stronger association with HCM (P=1.76×10(-9)). Although HCM patients were more likely to carry these, FHOD3 allele subjects homozygous for FHOD3-1151I had similar HCM phenotypes as carriers of the V1151 allele. FHOD3 expression is increased in the setting of HCM, and both alleles of FHOD3-V1151I were detected in HCM myectomy tissue. Previously, FHOD3 was found to be required for formation of the sarcomere, and here we demonstrate that its fly homolog fhos is required for normal adult heart systolic contraction. CONCLUSIONS: Here we demonstrate the association of a common nonsynonymous FHOD3 genetic variant with HCM. This discovery further strengthens the potential role of gene mutations and polymorphisms that alter the amino acid sequence of sarcomere proteins and HCM.

Application of gene network analysis techniques identifies AXIN1/PDIA2 and endoglin haplotypes associated with bicuspid aortic valve.

Bicuspid Aortic Valve (BAV) is a highly heritable congenital heart defect. The low frequency of BAV (1% of general population) limits our ability to perform genome-wide association studies. We present the application of four a priori SNP selection techniques, reducing the multiple-testing penalty by restricting analysis to SNPs relevant to BAV in a genome-wide SNP dataset from a cohort of 68 BAV probands and 830 control subjects. Two knowledge-based approaches, CANDID and STRING, were used to systematically identify BAV genes, and their SNPs, from the published literature, microarray expression studies and a genome scan. We additionally tested Functionally Interpolating SNPs (fitSNPs) present on the array; the fourth consisted of SNPs selected by Random Forests, a machine learning approach. These approaches reduced the multiple testing penalty by lowering the fraction of the genome probed to 0.19% of the total, while increasing the likelihood of studying SNPs within relevant BAV genes and pathways. Three loci were identified by CANDID, STRING, and fitSNPS. A haplotype within the AXIN1-PDIA2 locus (p-value of 2.926x10(-06)) and a haplotype within the Endoglin gene (p-value of 5.881x10(-04)) were found to be strongly associated with BAV. The Random Forests approach identified a SNP on chromosome 3 in association with BAV (p-value 5.061x10(-06)). The results presented here support an important role for genetic variants in BAV and provide support for additional studies in well-powered cohorts. Further, these studies demonstrate that leveraging existing expression and genomic data in the context of GWAS studies can identify biologically relevant genes and pathways associated with a congenital heart defect.