Targeted alveolar regeneration with Frizzled-specific agonists.

Wnt ligands oligomerize Frizzled (Fzd) and Lrp5/6 receptors to control the specification and activity of stem cells in many species. How Wnt signaling is selectively activated in different stem cell populations, often within the same organ, is not understood. In lung alveoli, we show that distinct Wnt receptors are expressed by epithelial (Fzd5/6), endothelial (Fzd4), and stromal (Fzd1) cells. Fzd5 is uniquely required for alveolar epithelial stem cell activity, whereas fibroblasts utilize distinct Fzd receptors. Using an expanded repertoire of Fzd-Lrp agonists, we could activate canonical Wnt signaling in alveolar epithelial stem cells via either Fzd5 or, unexpectedly, non-canonical Fzd6. A Fzd5 agonist (Fzd5ag) or Fzd6ag stimulated alveolar epithelial stem cell activity and promoted survival in mice after lung injury, but only Fzd6ag promoted an alveolar fate in airway-derived progenitors. Therefore, we identify a potential strategy for promoting regeneration without exacerbating fibrosis during lung injury.

A renal cell carcinoma tumorgraft platform to advance precision medicine.

Renal cell carcinoma (RCC) encompasses a heterogenous group of tumors, but representative preclinical models are lacking. We previously showed that patient-derived tumorgraft (TG) models recapitulate the biology and treatment responsiveness. Through systematic orthotopic implantation of tumor samples from 926 ethnically diverse individuals into non-obese diabetic (NOD)/severe combined immunodeficiency (SCID) mice, we generate a resource comprising 172 independently derived, stably engrafted TG lines from 148 individuals. TG lines are characterized histologically and genomically (whole-exome [n = 97] and RNA [n = 102] sequencing). The platform features a variety of histological and oncogenotypes, including TCGA clades further corroborated through orthogonal metabolomic analyses. We illustrate how it enables a deeper understanding of RCC biology; enables the development of tissue- and imaging-based molecular probes; and supports advances in drug development.

Determinants of renal cell carcinoma invasion and metastatic competence.

Metastasis is the principal cause of cancer related deaths. Tumor invasion is essential for metastatic spread. However, determinants of invasion are poorly understood. We addressed this knowledge gap by leveraging a unique attribute of kidney cancer. Renal tumors invade into large vessels forming tumor thrombi (TT) that migrate extending sometimes into the heart. Over a decade, we prospectively enrolled 83 ethnically-diverse patients undergoing surgical resection for grossly invasive tumors at UT Southwestern Kidney Cancer Program. In this study, we perform comprehensive histological analyses, integrate multi-region genomic studies, generate in vivo models, and execute functional studies to define tumor invasion and metastatic competence. We find that invasion is not always associated with the most aggressive clone. Driven by immediate early genes, invasion appears to be an opportunistic trait attained by subclones with diverse oncogenomic status in geospatial proximity to vasculature. We show that not all invasive tumors metastasize and identify determinants of metastatic competency. TT associated with metastases are characterized by higher grade, mTOR activation and a particular immune contexture. Moreover, TT grade is a better predictor of metastasis than overall tumor grade, which may have implications for clinical practice.

Retinoblastoma genetics screening and clinical management.

BACKGROUND: India accounts for 20% of the global retinoblastoma (RB) burden. However, the existing data on RB1 gene germline mutations and its influence on clinical decisions is minimally explored. METHODS: Fifty children with RB underwent complete clinical examination and appropriate multidisciplinary management. Screening of germline RB1 gene mutations was performed through next-generation sequencing and Multiplex Ligation-dependent Probe Amplification (MLPA) analysis. The mutation and non-mutation groups were compared for clinical parameters especially severity, progression and recurrence. RESULTS: Twenty-nine patients had bilateral RB (BLRB) and 21 had unilateral RB (ULRB). The genetic analysis revealed 20 RB1 variations in 29 probands, inclusive of 3 novel mutations, known 16 mutations and heterozygous whole gene deletions. The mutation detection rate (MDR) was 86.2% in BLRB and 19% in ULRB. Associations of disease recurrence (p = 0.021), progression (p = 0.000) and higher percentage of optic nerve invasion, subretinal seeds and high-risk pathological factors were observed in the mutation group. Clinical management was influenced by the presence of germline mutations, particularly while deciding on enucleation, frequency of periodic follow up and radiotherapy. CONCLUSIONS: We identified novel RB1 mutations, and our mutation detection rate was on par with the previous global studies. In our study, genetic results influenced clinical management and we suggest that it should be an essential and integral component of RB-care in India and elsewhere.

Anti-ferroptotic mechanism of IL4i1-mediated amino acid metabolism.

Interleukin-4-induced-1 (IL4i1) is an amino acid oxidase secreted from immune cells. Recent observations have suggested that IL4i1 is pro-tumorigenic via unknown mechanisms. As IL4i1 has homologs in snake venoms (L-amino acid oxidases [LAAO]), we used comparative approaches to gain insight into the mechanistic basis of how conserved amino acid oxidases regulate cell fate and function. Using mammalian expressed recombinant proteins, we found that venom LAAO kills cells via hydrogen peroxide generation. By contrast, mammalian IL4i1 is non-cytotoxic and instead elicits a cell protective gene expression program inhibiting ferroptotic redox death by generating indole-3-pyruvate (I3P) from tryptophan. I3P suppresses ferroptosis by direct free radical scavenging and through the activation of an anti-oxidative gene expression program. Thus, the pro-tumor effects of IL4i1 are likely mediated by local anti-ferroptotic pathways via aromatic amino acid metabolism, arguing that an IL4i1 inhibitor may modulate tumor cell death pathways.

Integrated genomic analysis reveals mutated ELF3 as a potential gallbladder cancer vaccine candidate.

Gallbladder cancer (GBC) is an aggressive gastrointestinal malignancy with no approved targeted therapy. Here, we analyze exomes (n = 160), transcriptomes (n = 115), and low pass whole genomes (n = 146) from 167 gallbladder cancers (GBCs) from patients in Korea, India and Chile. In addition, we also sequence samples from 39 GBC high-risk patients and detect evidence of early cancer-related genomic lesions. Among the several significantly mutated genes not previously linked to GBC are ETS domain genes ELF3 and EHF, CTNNB1, APC, NSD1, KAT8, STK11 and NFE2L2. A majority of ELF3 alterations are frame-shift mutations that result in several cancer-specific neoantigens that activate T-cells indicating that they are cancer vaccine candidates. In addition, we identify recurrent alterations in KEAP1/NFE2L2 and WNT pathway in GBC. Taken together, these define multiple targetable therapeutic interventions opportunities for GBC treatment and management.

Proteogenomics analysis unveils a TFG-RET gene fusion and druggable targets in papillary thyroid carcinomas.

Papillary thyroid cancer (PTC) is the most common type of endocrine malignancy. By RNA-seq analysis, we identify a RET rearrangement in the tumour material of a patient who does not harbour any known RAS or BRAF mutations. This new gene fusion involves exons 1-4 from the 5' end of the Trk fused Gene (TFG) fused to the 3' end of RET tyrosine kinase leading to a TFG-RET fusion which transforms immortalized human thyroid cells in a kinase-dependent manner. TFG-RET oligomerises in a PB1 domain-dependent manner and oligomerisation of TFG-RET is required for oncogenic transformation. Quantitative proteomic analysis reveals the upregulation of E3 Ubiquitin ligase HUWE1 and DUBs like USP9X and UBP7 in both tumor and metastatic lesions, which is further confirmed in additional patients. Expression of TFG-RET leads to the upregulation of HUWE1 and inhibition of HUWE1 significantly reduces RET-mediated oncogenesis.

Exome sequencing reveals a novel splice site variant in HUWE1 gene in patients with suspected Say-Meyer syndrome.

Say-Meyer syndrome is a rare and clinically heterogeneous syndrome characterized by trigonocephaly, short stature, developmental delay and hypotelorism. Nine patients with this syndrome have been reported thus far although no causative gene has yet been identified. Here, we report two siblings with clinical phenotypes of Say-Meyer syndrome with moderate to severe intellectual disability and autism spectrum disorder. Cytogenetics and array-based comparative genomic hybridization did not reveal any chromosome abnormalities or copy number alterations. Exome sequencing of the patients revealed a novel X-linked recessive splice acceptor site variant c.145-2A > G in intron 5 of HUWE1 gene in both affected siblings. RT-PCR and sequencing revealed the use of an alternate cryptic splice acceptor site downstream, which led to deletion of six nucleotides resulting loss of two amino acids p.(Cys49-Glu50del) in HUWE1 protein. Deletion of these two amino acids, which are located in a highly conserved region, is predicted to be deleterious and quite likely to affect the function of HUWE1 protein. This is the first report of a potential candidate gene mutation for Say-Meyer syndrome, which was initially described four decades ago.

S100a4 upregulation in Pik3caH1047R;Trp53R270H;MMTV-Cre-driven mammary tumors promotes metastasis.

BACKGROUND: PIK3CA mutations are frequent in human breast cancer. Pik3caH1047R mutant expression in mouse mammary gland promotes tumorigenesis. TP53 mutations co-occur with PIK3CA mutations in human breast cancers. We previously generated a conditionally activatable Pik3caH1047R;MMTV-Cre mouse model and found a few malignant sarcomatoid (spindle cell) carcinomas that had acquired spontaneous dominant-negative Trp53 mutations. METHODS: A Pik3caH1047R;Trp53R270H;MMTV-Cre double mutant mouse breast cancer model was generated. Tumors were characterized by histology, marker analysis, transcriptional profiling, single-cell RNA-seq, and bioinformatics. Cell lines were developed from mutant tumors and used to identify and confirm genes involved in metastasis. RESULTS: We found Pik3caH1047R and Trp53R270H cooperate in driving oncogenesis in mammary glands leading to a shorter latency than either alone. Double mutant mice develop multiple histologically distinct mammary tumors, including adenocarcinoma and sarcomatoid (spindle cell) carcinoma. We found some tumors to be invasive and a few metastasized to the lung and/or the lymph node. Single-cell RNA-seq analysis of the tumors identified epithelial, stromal, myeloid, and T cell groups. Expression analysis of the metastatic tumors identified S100a4 as a top candidate gene associated with metastasis. Metastatic tumors contained a much higher percentage of epithelial-mesenchymal transition (EMT)-signature positive and S100a4-expressing cells. CRISPR/CAS9-mediated knockout of S100a4 in a metastatic tumor-derived cell line disrupted its metastatic potential indicating a role for S100a4 in metastasis. CONCLUSIONS: Pik3caH1047R;Trp53R270H;MMTV-Cre mouse provides a preclinical model to mimic a subtype of human breast cancers that carry both PIK3CA and TP53 mutations. It also allows for understanding the cooperation between the two mutant genes in tumorigenesis. Our model also provides a system to study metastasis and develop therapeutic strategies for PIK3CA/TP53 double-positive cancers. S100a4 found involved in metastasis in this model can be a potential diagnostic and therapeutic target.

SCINA: A Semi-Supervised Subtyping Algorithm of Single Cells and Bulk Samples.

Advances in single-cell RNA sequencing (scRNA-Seq) have allowed for comprehensive analyses of single cell data. However, current analyses of scRNA-Seq data usually start from unsupervised clustering or visualization. These methods ignore prior knowledge of transcriptomes and the probable structures of the data. Moreover, cell identification heavily relies on subjective and possibly inaccurate human inspection afterwards. To address these analytical challenges, we developed SCINA (Semi-supervised Category Identification and Assignment), a semi-supervised model that exploits previously established gene signatures using an expectation-maximization (EM) algorithm. SCINA is applicable to scRNA-Seq and flow cytometry/CyTOF data, as well as other data of similar format. We applied SCINA to a wide range of datasets, and showed its accuracy, stability and efficiency, which exceeded most popular unsupervised approaches. SCINA discovered an intermediate stage of oligodendrocytes from mouse brain scRNA-Seq data. SCINA also detected immune cell population changes in cytometry data in a genetically-engineered mouse model. Furthermore, SCINA performed well with bulk gene expression data. Specifically, we identified a new kidney tumor clade with similarity to FH-deficient tumors (FHD), which we refer to as FHD-like tumors (FHDL). Overall, SCINA provides both methodological advances and biological insights from perspectives different from traditional analytical methods.

Actionable Activating Oncogenic ERBB2/HER2 Transmembrane and Juxtamembrane Domain Mutations.

Deregulated HER2 is a target of many approved cancer drugs. We analyzed 111,176 patient tumors and identified recurrent mutations in HER2 transmembrane domain (TMD) and juxtamembrane domain (JMD) that include G660D, R678Q, E693K, and Q709L. Using a saturation mutagenesis screen and testing of patient-derived mutations we found several activating TMD and JMD mutations. Structural modeling and analysis showed that the TMD/JMD mutations function by improving the active dimer interface or stabilizing an activating conformation. Further, we found that HER2 G660D employed asymmetric kinase dimerization for activation and signaling. Importantly, anti-HER2 antibodies and small-molecule kinase inhibitors blocked the activity of TMD/JMD mutants. Consistent with this, a G660D germline mutant lung cancer patient showed remarkable clinical response to HER2 blockade.

An Empirical Approach Leveraging Tumorgrafts to Dissect the Tumor Microenvironment in Renal Cell Carcinoma Identifies Missing Link to Prognostic Inflammatory Factors.

By leveraging tumorgraft (patient-derived xenograft) RNA-sequencing data, we developed an empirical approach, DisHet, to dissect the tumor microenvironment (eTME). We found that 65% of previously defined immune signature genes are not abundantly expressed in renal cell carcinoma (RCC) and identified 610 novel immune/stromal transcripts. Using eTME, genomics, pathology, and medical record data involving >1,000 patients, we established an inflamed pan-RCC subtype (IS) enriched for regulatory T cells, natural killer cells, TH1 cells, neutrophils, macrophages, B cells, and CD8+ T cells. IS is enriched for aggressive RCCs, including BAP1-deficient clear-cell and type 2 papillary tumors. The IS subtype correlated with systemic manifestations of inflammation such as thrombocytosis and anemia, which are enigmatic predictors of poor prognosis. Furthermore, IS was a strong predictor of poor survival. Our analyses suggest that tumor cells drive the stromal immune response. These data provide a missing link between tumor cells, the TME, and systemic factors.Significance: We undertook a novel empirical approach to dissect the renal cell carcinoma TME by leveraging tumorgrafts. The dissection and downstream analyses uncovered missing links between tumor cells, the TME, systemic manifestations of inflammation, and poor prognosis. Cancer Discov; 8(9); 1142-55. ©2018 AACR.This article is highlighted in the In This Issue feature, p. 1047.

ERK Mutations and Amplification Confer Resistance to ERK-Inhibitor Therapy.

Purpose: MAPK pathway inhibitors targeting BRAF and MEK have shown clinical efficacy in patients with RAF- and/or RAS-mutated tumors. However, acquired resistance to these agents has been an impediment to improved long-term survival in the clinic. In such cases, targeting ERK downstream of BRAF/MEK has been proposed as a potential strategy for overcoming acquired resistance. Preclinical studies suggest that ERK inhibitors are effective at inhibiting BRAF/RAS-mutated tumor growth and overcome BRAF or/and MEK inhibitor resistance. However, as observed with other MAPK pathway inhibitors, treatment with ERK inhibitors is likely to cause resistance in the clinic. Here, we aimed to model the mechanism of resistance to ERK inhibitors.Experimental Design: We tested five structurally different ATP-competitive ERK inhibitors representing three different scaffolds on BRAF/RAS-mutant cancer cell lines of different tissue types to generate resistant lines. We have used in vitro modeling, structural biology, and genomic analysis to understand the development of resistance to ERK inhibitors and the mechanisms leading to it.Results: We have identified mutations in ERK1/2, amplification and overexpression of ERK2, and overexpression of EGFR/ERBB2 as mechanisms of acquired resistance. Structural analysis of ERK showed that specific compounds that induced on-target ERK mutations were impaired in their ability to bind mutant ERK. We show that in addition to MEK inhibitors, ERBB receptor and PI3K/mTOR pathway inhibitors are effective in overcoming ERK-inhibitor resistance.Conclusions: These findings suggest that combination therapy with MEK or ERBB receptor or PI3K/mTOR and ERK inhibitors may be an effective strategy for managing the emergence of resistance in the clinic. Clin Cancer Res; 24(16); 4044-55. ©2018 AACR.

Comprehensive genomic analysis identifies pathogenic variants in maturity-onset diabetes of the young (MODY) patients in South India.

BACKGROUND: Maturity-onset diabetes of the young (MODY) is an early-onset, autosomal dominant form of non-insulin dependent diabetes. Genetic diagnosis of MODY can transform patient management. Earlier data on the genetic predisposition to MODY have come primarily from familial studies in populations of European origin. METHODS: In this study, we carried out a comprehensive genomic analysis of 289 individuals from India that included 152 clinically diagnosed MODY cases to identify variants in known MODY genes. Further, we have analyzed exome data to identify putative MODY relevant variants in genes previously not implicated in MODY. Functional validation of MODY relevant variants was also performed. RESULTS: We found MODY 3 (HNF1A; 7.2%) to be most frequently mutated followed by MODY 12 (ABCC8; 3.3%). They together account for ~ 11% of the cases. In addition to known MODY genes, we report the identification of variants in RFX6, WFS1, AKT2, NKX6-1 that may contribute to development of MODY. Functional assessment of the NKX6-1 variants showed that they are functionally impaired. CONCLUSIONS: Our findings showed HNF1A and ABCC8 to be the most frequently mutated MODY genes in south India. Further we provide evidence for additional MODY relevant genes, such as NKX6-1, and these require further validation.

Recurrent Loss of NFE2L2 Exon 2 Is a Mechanism for Nrf2 Pathway Activation in Human Cancers.

The Nrf2 pathway is frequently activated in human cancers through mutations in Nrf2 or its negative regulator KEAP1. Using a cell-line-derived gene signature for Nrf2 pathway activation, we found that some tumors show high Nrf2 activity in the absence of known mutations in the pathway. An analysis of splice variants in oncogenes revealed that such tumors express abnormal transcript variants from the NFE2L2 gene (encoding Nrf2) that lack exon 2, or exons 2 and 3, and encode Nrf2 protein isoforms missing the KEAP1 interaction domain. The Nrf2 alterations result in the loss of interaction with KEAP1, Nrf2 stabilization, induction of a Nrf2 transcriptional response, and Nrf2 pathway dependence. In all analyzed cases, transcript variants were the result of heterozygous genomic microdeletions. Thus, we identify an alternative mechanism for Nrf2 pathway activation in human tumors and elucidate its functional consequences.

Comprehensive genomic analysis of malignant pleural mesothelioma identifies recurrent mutations, gene fusions and splicing alterations.

We analyzed transcriptomes (n = 211), whole exomes (n = 99) and targeted exomes (n = 103) from 216 malignant pleural mesothelioma (MPM) tumors. Using RNA-seq data, we identified four distinct molecular subtypes: sarcomatoid, epithelioid, biphasic-epithelioid (biphasic-E) and biphasic-sarcomatoid (biphasic-S). Through exome analysis, we found BAP1, NF2, TP53, SETD2, DDX3X, ULK2, RYR2, CFAP45, SETDB1 and DDX51 to be significantly mutated (q-score ≥ 0.8) in MPMs. We identified recurrent mutations in several genes, including SF3B1 (∼2%; 4/216) and TRAF7 (∼2%; 5/216). SF3B1-mutant samples showed a splicing profile distinct from that of wild-type tumors. TRAF7 alterations occurred primarily in the WD40 domain and were, except in one case, mutually exclusive with NF2 alterations. We found recurrent gene fusions and splice alterations to be frequent mechanisms for inactivation of NF2, BAP1 and SETD2. Through integrated analyses, we identified alterations in Hippo, mTOR, histone methylation, RNA helicase and p53 signaling pathways in MPMs.

A comprehensive transcriptional portrait of human cancer cell lines.

Tumor-derived cell lines have served as vital models to advance our understanding of oncogene function and therapeutic responses. Although substantial effort has been made to define the genomic constitution of cancer cell line panels, the transcriptome remains understudied. Here we describe RNA sequencing and single-nucleotide polymorphism (SNP) array analysis of 675 human cancer cell lines. We report comprehensive analyses of transcriptome features including gene expression, mutations, gene fusions and expression of non-human sequences. Of the 2,200 gene fusions catalogued, 1,435 consist of genes not previously found in fusions, providing many leads for further investigation. We combine multiple genome and transcriptome features in a pathway-based approach to enhance prediction of response to targeted therapeutics. Our results provide a valuable resource for studies that use cancer cell lines.

Spectrum of diverse genomic alterations define non-clear cell renal carcinoma subtypes.

To further understand the molecular distinctions between kidney cancer subtypes, we analyzed exome, transcriptome and copy number alteration data from 167 primary human tumors that included renal oncocytomas and non-clear cell renal cell carcinomas (nccRCCs), consisting of papillary (pRCC), chromophobe (chRCC) and translocation (tRCC) subtypes. We identified ten significantly mutated genes in pRCC, including MET, NF2, SLC5A3, PNKD and CPQ. MET mutations occurred in 15% (10/65) of pRCC samples and included previously unreported recurrent activating mutations. In chRCC, we found TP53, PTEN, FAAH2, PDHB, PDXDC1 and ZNF765 to be significantly mutated. Gene expression analysis identified a five-gene set that enabled the molecular classification of chRCC, renal oncocytoma and pRCC. Using RNA sequencing, we identified previously unreported gene fusions, including ACTG1-MITF fusion. Ectopic expression of the ACTG1-MITF fusion led to cellular transformation and induced the expression of downstream target genes. Finally, we observed upregulation of the anti-apoptotic factor BIRC7 in MiTF-high RCC tumors, suggesting a potential therapeutic role for BIRC7 inhibitors.

Diverse modes of genomic alteration in hepatocellular carcinoma.

BACKGROUND: Hepatocellular carcinoma (HCC) is a heterogeneous disease with high mortality rate. Recent genomic studies have identified TP53, AXIN1, and CTNNB1 as the most frequently mutated genes. Lower frequency mutations have been reported in ARID1A, ARID2 and JAK1. In addition, hepatitis B virus (HBV) integrations into the human genome have been associated with HCC. RESULTS: Here, we deep-sequence 42 HCC patients with a combination of whole genome, exome and transcriptome sequencing to identify the mutational landscape of HCC using a reasonably large discovery cohort. We find frequent mutations in TP53, CTNNB1 and AXIN1, and rare but likely functional mutations in BAP1 and IDH1. Besides frequent hepatitis B virus integrations at TERT, we identify translocations at the boundaries of TERT. A novel deletion is identified in CTNNB1 in a region that is heavily mutated in multiple cancers. We also find multiple high-allelic frequency mutations in the extracellular matrix protein LAMA2. Lower expression levels of LAMA2 correlate with a proliferative signature, and predict poor survival and higher chance of cancer recurrence in HCC patients, suggesting an important role of the extracellular matrix and cell adhesion in tumor progression of a subgroup of HCC patients. CONCLUSIONS: The heterogeneous disease of HCC features diverse modes of genomic alteration. In addition to common point mutations, structural variations and methylation changes, there are several virus-associated changes, including gene disruption or activation, formation of chimeric viral-human transcripts, and DNA copy number changes. Such a multitude of genomic events likely contributes to the heterogeneous nature of HCC.

Dimerization of the kinase ARAF promotes MAPK pathway activation and cell migration.

The RAF family of kinases mediates RAS signaling, and RAF inhibitors can be effective for treating tumors with BRAF(V600E) mutant protein. However, RAF inhibitors paradoxically accelerate metastasis in RAS-mutant tumors and become ineffective in BRAF(V600E) tumors because of reactivation of downstream mitogen-activated protein kinase (MAPK) signaling. We found that the RAF isoform ARAF has an obligatory role in promoting MAPK activity and cell migration in a cell type-dependent manner. Knocking down ARAF prevented the activation of MAPK kinase 1 (MEK1) and extracellular signal-regulated kinase 1 and 2 (ERK1/2) and decreased the number of protrusions from tumor cell spheroids in three-dimensional culture that were induced by BRAF(V600E)-specific or BRAF/CRAF inhibitors (GDC-0879 and sorafenib, respectively). RAF inhibitors induced the homodimerization of ARAF and the heterodimerization of BRAF with CRAF and the scaffolding protein KSR1. In a purified protein solution, recombinant proteins of the three RAF isoforms competed for binding to MEK1. In cells in culture, overexpressing mutants of ARAF that could not homodimerize impaired the interaction between ARAF and endogenous MEK1 and thus prevented the subsequent activation of MEK1 and ERK1/2. Our findings reveal a new role for ARAF in directly activating the MAPK cascade and promoting tumor cell invasion and suggest a new therapeutic target for RAS- and RAF-mediated cancers.