Analysis of protein-altering variants in telomerase genes and their association with MUC5B common variant status in patients with idiopathic pulmonary fibrosis: a candidate gene sequencing study.

BACKGROUND: Idiopathic pulmonary fibrosis (IPF) risk has a strong genetic component. Studies have implicated variations at several loci, including TERT, surfactant genes, and a single nucleotide polymorphism at chr11p15 (rs35705950) in the intergenic region between TOLLIP and MUC5B. Patients with IPF who have risk alleles at rs35705950 have longer survival from the time of IPF diagnosis than do patients homozygous for the non-risk allele, whereas patients with shorter telomeres have shorter survival times. We aimed to assess whether rare protein-altering variants in genes regulating telomere length are enriched in patients with IPF homozygous for the non-risk alleles at rs35705950. METHODS: Between Nov 1, 2014, and Nov 1, 2016, we assessed blood samples from patients aged 40 years or older and of European ancestry with sporadic IPF from three international phase 3 clinical trials (INSPIRE, CAPACITY, ASCEND), one phase 2 study (RIFF), and US-based observational studies (Vanderbilt Clinical Interstitial Lung Disease Registry and the UCSF Interstitial Lung Disease Clinic registry cohorts) at the Broad Institute (Cambridge, MA, USA) and Human Longevity (San Diego, CA, USA). We also assessed blood samples from non-IPF controls in several clinical trials. We did whole-genome sequencing to assess telomere length and identify rare protein-altering variants, stratified by rs35705950 genotype. We also assessed rare functional variation in TERT exons and compared telomere length and disease progression across genotypes. FINDINGS: We assessed samples from 1510 patients with IPF and 1874 non-IPF controls. 30 (3%) of 1046 patients with an rs35705950 risk allele had a rare protein-altering variant in TERT compared with 34 (7%) of 464 non-risk allele carriers (odds ratio 0·40 [95% CI 0·24-0·66], p=0·00039). Subsequent analyses identified enrichment of rare protein-altering variants in PARN and RTEL1, and rare variation in TERC in patients with IPF compared with controls. We expanded our study population to provide a more accurate estimation of rare variant frequency in these four loci, and to calculate telomere length. The proportion of patients with at least one rare variant in TERT, PARN, TERC, or RTEL1 was higher in patients with IPF than in controls (149 [9%] of 1739 patients vs 205 [2%] of 8645 controls, p=2·44 × 10-8). Patients with IPF who had a variant in any of the four identified telomerase component genes had telomeres that were 3·69-16·10% shorter than patients without a variant in any of the four genes and had an earlier mean age of disease onset than patients without one or more variants (65·1 years [SD 7·8] vs 67·1 years [7·9], p=0·004). In the placebo arms of clinical trials, shorter telomeres were significantly associated with faster disease progression (1·7% predicted forced vital capacity per kb per year, p=0·002). Pirfenidone had treatment benefit regardless of telomere length (p=4·24 × 10-8 for telomere length lower than the median, p=0·0044 for telomere length greater than the median). INTERPRETATION: Rare protein-altering variants in TERT, PARN, TERC, and RTEL1 are enriched in patients with IPF compared with controls, and, in the case of TERT, particularly in individuals without a risk allele at the rs35705950 locus. This suggests that multiple genetic factors contribute to sporadic IPF, which might implicate distinct mechanisms of pathogenesis and disease progression. FUNDING: Genentech, National Institutes of Health, Francis Family Foundation, Pulmonary Fibrosis Foundation, Nina Ireland Program for Lung Health, US Department of Veterans Affairs.

Constitutive NOTCH3 Signaling Promotes the Growth of Basal Breast Cancers.

Notch ligands signal through one of four receptors on neighboring cells to mediate cell-cell communication and control cell fate, proliferation, and survival. Although aberrant Notch activation has been implicated in numerous malignancies, including breast cancer, the importance of individual receptors in distinct breast cancer subtypes and the mechanisms of receptor activation remain unclear. Using a novel antibody to detect active NOTCH3, we report here that NOTCH3 signals constitutively in a panel of basal breast cancer cell lines and in more than one third of basal tumors. Selective inhibition of individual ligands revealed that this signal does not require canonical ligand induction. A NOTCH3 antagonist antibody inhibited growth of basal lines, whereas a NOTCH3 agonist antibody enhanced the transformed phenotype in vitro and in tumor xenografts. Transcriptomic analyses generated a Notch gene signature that included Notch pathway components, the oncogene c-Myc, and the mammary stem cell regulator Id4 This signature drove clustering of breast cancer cell lines and tumors into the common subtypes and correlated with the basal classification. Our results highlight an unexpected ligand-independent induction mechanism and suggest that constitutive NOTCH3 signaling can drive an oncogenic program in a subset of basal breast cancers. Cancer Res; 77(6); 1439-52. ©2017 AACR.

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.

Identification and analysis of in vivo VEGF downstream markers link VEGF pathway activity with efficacy of anti-VEGF therapies.

PURPOSE: The aim of this study was to identify conserved pharmacodynamic and potential predictive biomarkers of response to anti-VEGF therapy using gene expression profiling in preclinical tumor models and in patients. EXPERIMENTAL DESIGN: Surrogate markers of VEGF inhibition [VEGF-dependent genes or VEGF-dependent vasculature (VDV)] were identified by profiling gene expression changes induced in response to VEGF blockade in preclinical tumor models and in human biopsies from patients treated with anti-VEGF monoclonal antibodies. The potential value of VDV genes as candidate predictive biomarkers was tested by correlating high or low VDV gene expression levels in pretreatment clinical samples with the subsequent clinical efficacy of bevacizumab (anti-VEGF)-containing therapy. RESULTS: We show that VDV genes, including direct and more distal VEGF downstream endothelial targets, enable detection of VEGF signaling inhibition in mouse tumor models and human tumor biopsies. Retrospective analyses of clinical trial data indicate that patients with higher VDV expression in pretreatment tumor samples exhibited improved clinical outcome when treated with bevacizumab-containing therapies. CONCLUSIONS: In this work, we identified surrogate markers (VDV genes) for in vivo VEGF signaling in tumors and showed clinical data supporting a correlation between pretreatment VEGF bioactivity and the subsequent efficacy of anti-VEGF therapy. We propose that VDV genes are candidate biomarkers with the potential to aid the selection of novel indications as well as patients likely to respond to anti-VEGF therapy. The data presented here define a diagnostic biomarker hypothesis based on translational research that warrants further evaluation in additional retrospective and prospective trials.

Development of a robust flow cytometry-based pharmacodynamic assay to detect phospho-protein signals for phosphatidylinositol 3-kinase inhibitors in multiple myeloma.

BACKGROUND: The phosphatidylinositol 3-kinase (PI3K) pathway plays an important role in multiple myeloma (MM), a blood cancer associated with uncontrolled proliferation of bone marrow plasma cells. This study aimed to develop a robust clinical pharmacodynamic (PD) assay to measure the on-target PD effects of the selective PI3K inhibitor GDC-0941 in MM patients. METHODS: We conducted an in vitro drug wash-out study to evaluate the feasibility of biochemical approaches in measuring the phosphorylation of S6 ribosomal protein (S6), one of the commonly used PD markers for PI3K pathway inhibition. We then developed a 7-color phospho-specific flow cytometry assay, or phospho flow assay, to measure the phosphorylation state of intracellular S6 in bone marrow aspirate (BMA) and peripheral blood (PB). Integrated mean fluorescence intensity (iMFI) was used to calculate fold changes of phosphorylation. Assay sensitivity was evaluated by comparing phospho flow with Meso Scale Discovery (MSD) and immunohistochemistry (IHC) assays. Finally, a sample handling method was developed to maintain the integrity of phospho signal during sample shipping and storage to ensure clinical application. RESULTS: The phospho flow assay provided single-cell PD monitoring of S6 phosphorylation in tumor and surrogate cells using fixed BMA and PB, assessing pathway modulation in response to GDC-0941 with sensitivity similar to that of MSD assay. The one-shot sample fixation and handling protocol herein demonstrated exceptional preservation of protein phosphorylation. In contrast, the IHC assay was less sensitive in terms of signal quantification while the biochemical approach (MSD) was less suitable to assess PD activities due to the undesirable impact associated with cell isolation on the protein phosphorylation in tumor cells. CONCLUSIONS: We developed a robust PD biomarker assay for the clinical evaluation of PI3K inhibitors in MM, allowing one to decipher the PD response in a relevant cell population. To our knowledge, this is the first report of an easily implemented clinical PD assay that incorporates an unbiased one-shot sample handling protocol, all (staining)-in-one (tube) phospho flow staining protocol, and an integrated modified data analysis for PD monitoring of kinase inhibitors in relevant cell populations in BMA and PB. The methods described here ensure a real-time, reliable and reproducible PD readout, which can provide information for dose selection as well as help to identify optimal combinations of targeted agents in early clinical trials.

A high-content cellular senescence screen identifies candidate tumor suppressors, including EPHA3.

Activation of a cellular senescence program is a common response to prolonged oncogene activation or tumor suppressor loss, providing a physiological mechanism for tumor suppression in premalignant cells. The link between senescence and tumor suppression supports the hypothesis that a loss-of-function screen measuring bona fide senescence marker activation should identify candidate tumor suppressors. Using a high-content siRNA screening assay for cell morphology and proliferation measures, we identify 12 senescence-regulating kinases and determine their senescence marker signatures, including elevation of senescence-associated ß-galactosidase, DNA damage and p53 or p16 (INK4a) expression. Consistent with our hypothesis, SNP array CGH data supports loss of gene copy number of five senescence-suppressing genes across multiple tumor samples. One such candidate is the EPHA3 receptor tyrosine kinase, a gene commonly mutated in human cancer. We demonstrate that selected intracellular EPHA3 tumor-associated point mutations decrease receptor expression level and/or receptor tyrosine kinase (RTK) activity. Our study therefore describes a new strategy to mine for novel candidate tumor suppressors and provides compelling evidence that EPHA3 mutations may promote tumorigenesis only when key senescence-inducing pathways have been inactivated.

Identification of mutant K-Ras-dependent phenotypes using a panel of isogenic cell lines.

To assess the consequences of endogenous mutant K-Ras, we analyzed the signaling and biological properties of a small panel of isogenic cell lines. These include the cancer cell lines DLD1, HCT116, and Hec1A, in which either the WT or mutant K-ras allele has been disrupted, and SW48 colorectal cancer cells and human mammary epithelial cells in which a single copy of mutant K-ras was introduced at its endogenous genomic locus. We find that single copy mutant K-Ras causes surprisingly modest activation of downstream signaling to ERK and Akt. In contrast, a negative feedback signaling loop to EGFR and N-Ras occurs in some, but not all, of these cell lines. Mutant K-Ras also had relatively minor effects on cell proliferation and cell migration but more dramatic effects on cell transformation as assessed by growth in soft agar. Surprisingly, knock-out of the wild type K-ras allele consistently increased growth in soft agar, suggesting tumor-suppressive properties of this gene under these conditions. Finally, we examined the effects of single copy mutant K-Ras on global gene expression. Although transcriptional programs triggered by mutant K-Ras were generally quite distinct in the different cell lines, there was a small number of genes that were consistently overexpressed, and these could be used to monitor K-Ras inhibition in a panel of human tumor cell lines. We conclude that there are conserved components of mutant K-Ras signaling and phenotypes but that many depend on cell context and environmental cues.

The STARD9/Kif16a kinesin associates with mitotic microtubules and regulates spindle pole assembly.

During cell division, cells form the microtubule-based mitotic spindle, a highly specialized and dynamic structure that mediates proper chromosome transmission to daughter cells. Cancer cells can show perturbed mitotic spindles and an approach in cancer treatment has been to trigger cell killing by targeting microtubule dynamics or spindle assembly. To identify and characterize proteins necessary for spindle assembly, and potential antimitotic targets, we performed a proteomic and genetic analysis of 592 mitotic microtubule copurifying proteins (MMCPs). Screening for regulators that affect both mitosis and apoptosis, we report the identification and characterization of STARD9, a kinesin-3 family member, which localizes to centrosomes and stabilizes the pericentriolar material (PCM). STARD9-depleted cells have fragmented PCM, form multipolar spindles, activate the spindle assembly checkpoint (SAC), arrest in mitosis, and undergo apoptosis. Interestingly, STARD9-depletion synergizes with the chemotherapeutic agent taxol to increase mitotic death, demonstrating that STARD9 is a mitotic kinesin and a potential antimitotic target.

Mapping the NPHP-JBTS-MKS protein network reveals ciliopathy disease genes and pathways.

Nephronophthisis (NPHP), Joubert (JBTS), and Meckel-Gruber (MKS) syndromes are autosomal-recessive ciliopathies presenting with cystic kidneys, retinal degeneration, and cerebellar/neural tube malformation. Whether defects in kidney, retinal, or neural disease primarily involve ciliary, Hedgehog, or cell polarity pathways remains unclear. Using high-confidence proteomics, we identified 850 interactors copurifying with nine NPHP/JBTS/MKS proteins and discovered three connected modules: "NPHP1-4-8" functioning at the apical surface, "NPHP5-6" at centrosomes, and "MKS" linked to Hedgehog signaling. Assays for ciliogenesis and epithelial morphogenesis in 3D renal cultures link renal cystic disease to apical organization defects, whereas ciliary and Hedgehog pathway defects lead to retinal or neural deficits. Using 38 interactors as candidates, linkage and sequencing analysis of 250 patients identified ATXN10 and TCTN2 as new NPHP-JBTS genes, and our Tctn2 mouse knockout shows neural tube and Hedgehog signaling defects. Our study further illustrates the power of linking proteomic networks and human genetics to uncover critical disease pathways.

Coordination of growth rate, cell cycle, stress response, and metabolic activity in yeast.

We studied the relationship between growth rate and genome-wide gene expression, cell cycle progression, and glucose metabolism in 36 steady-state continuous cultures limited by one of six different nutrients (glucose, ammonium, sulfate, phosphate, uracil, or leucine). The expression of more than one quarter of all yeast genes is linearly correlated with growth rate, independent of the limiting nutrient. The subset of negatively growth-correlated genes is most enriched for peroxisomal functions, whereas positively correlated genes mainly encode ribosomal functions. Many (not all) genes associated with stress response are strongly correlated with growth rate, as are genes that are periodically expressed under conditions of metabolic cycling. We confirmed a linear relationship between growth rate and the fraction of the cell population in the G0/G1 cell cycle phase, independent of limiting nutrient. Cultures limited by auxotrophic requirements wasted excess glucose, whereas those limited on phosphate, sulfate, or ammonia did not; this phenomenon (reminiscent of the "Warburg effect" in cancer cells) was confirmed in batch cultures. Using an aggregate of gene expression values, we predict (in both continuous and batch cultures) an "instantaneous growth rate." This concept is useful in interpreting the system-level connections among growth rate, metabolism, stress, and the cell cycle.