Molecular profiling of actionable mutations in refractory cancer patients has the potential to enable "precision medicine," wherein individualized therapies are guided based on genomic profiling. The molecular-screening program was intended to route participants to different candidate drugs in trials based on clinical-sequencing reports. In this screening program, we used a custom target-enrichment panel consisting of cancer-related genes to interrogate single-nucleotide variants, insertions and deletions, copy number variants, and a subset of gene fusions. From August 2014 through April 2015, 654 patients consented to participate in the program at Samsung Medical Center. Of these patients, 588 passed the quality control process for the 381-gene cancer-panel test, and 418 patients were included in the final analysis as being eligible for any anticancer treatment (127 gastric cancer, 122 colorectal cancer, 62 pancreatic/biliary tract cancer, 67 sarcoma/other cancer, and 40 genitourinary cancer patients). Of the 418 patients, 55 (12%) harbored a biomarker that guided them to a biomarker-selected clinical trial, and 184 (44%) patients harbored at least one genomic alteration that was potentially targetable. This study demonstrated that the panel-based sequencing program resulted in an increased rate of trial enrollment of metastatic cancer patients into biomarker-selected clinical trials. Given the expanding list of biomarker-selected trials, the guidance percentage to matched trials is anticipated to increase. IMPLICATIONS FOR PRACTICE: This study demonstrated that the panel-based sequencing program resulted in an increased rate of trial enrollment of metastatic cancer patients into biomarker-selected clinical trials. Given the expanding list of biomarker-selected trials, the guidance percentage to matched trials is anticipated to increase.
Purpose Histologic transformation of EGFR mutant lung adenocarcinoma (LADC) into small-cell lung cancer (SCLC) has been described as one of the major resistant mechanisms for epidermal growth factor receptor (EGFR) tyrosine kinase inhibitors (TKIs). However, the molecular pathogenesis is still unclear. Methods We investigated 21 patients with advanced EGFR-mutant LADCs that were transformed into EGFR TKI-resistant SCLCs. Among them, whole genome sequencing was applied for nine tumors acquired at various time points from four patients to reconstruct their clonal evolutionary history and to detect genetic predictors for small-cell transformation. The findings were validated by immunohistochemistry in 210 lung cancer tissues. Results We identified that EGFR TKI-resistant LADCs and SCLCs share a common clonal origin and undergo branched evolutionary trajectories. The clonal divergence of SCLC ancestors from the LADC cells occurred before the first EGFR TKI treatments, and the complete inactivation of both RB1 and TP53 were observed from the early LADC stages in sequenced tumors. We extended the findings by immunohistochemistry in the early-stage LADC tissues of 75 patients treated with EGFR TKIs; inactivation of both Rb and p53 was strikingly more frequent in the small-cell-transformed group than in the nontransformed group (82% v 3%; odds ratio, 131; 95% CI, 19.9 to 859). Among patients registered in a predefined cohort (n = 65), an EGFR mutant LADC that harbored completely inactivated Rb and p53 had a 43× greater risk of small-cell transformation (relative risk, 42.8; 95% CI, 5.88 to 311). Branch-specific mutational signature analysis revealed that apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like (APOBEC)-induced hypermutation was frequent in the branches toward small-cell transformation. Conclusion EGFR TKI-resistant SCLCs are branched out early from the LADC clones that harbor completely inactivated RB1 and TP53. The evaluation of RB1 and TP53 status in EGFR TKI-treated LADCs is informative in predicting small-cell transformation.
The Polycomb group (PcG) proteins are key conserved regulators of development, initially discovered in Drosophila and now strongly implicated in human disease. Nevertheless, differing silencing properties between the Drosophila and mammalian PcG systems have been observed. While specific DNA targeting sites for PcG proteins called Polycomb response elements (PREs) have been identified only in Drosophila, involvement of non-coding RNAs for PcG targeting has been favored in mammals. Another difference lies in the distribution patterns of PcG proteins. In mouse and human cells, PcG proteins show broad distributions, significantly overlapping with H3K27me3 domains. In contrast, only sharp peaks on PRE regions are observed for most PcG proteins in Drosophila, raising the question of how large domains of H3K27me3, up to many tens of kilobases, are formed and maintained in Drosophila. In this Extra View, we provide evidence that PcG distributions on silent chromatin in Drosophila are considerably broader than previously detected. Using BioTAP-XL, a chromatin crosslinking and tandem affinity purification approach, we find a broad, rather than PRE-limited overlap of PcG proteins with H3K27me3, suggesting a conserved spreading mechanism for PcG in flies and mammals.
Intravenous leiomyomatosis is an unusual smooth muscle neoplasm with quasi-malignant intravascular growth but a histologically banal appearance. Herein, we report expression and molecular cytogenetic analyses of a series of 12 intravenous leiomyomatosis cases to better understand the pathogenesis of intravenous leiomyomatosis. All cases were analyzed for the expression of HMGA2, MDM2, and CDK4 proteins by immunohistochemistry based on our previous finding of der(14)t(12;14)(q14.3;q24) in intravenous leiomyomatosis. Seven of 12 (58%) intravenous leiomyomatosis cases expressed HMGA2, and none expressed MDM2 or CDK4. Colocalization of hybridization signals for probes from the HMGA2 locus (12q14.3) and from 14q24 by interphase fluorescence in situ hybridization (FISH) was detected in a mean of 89.2% of nuclei in HMGA2-positive cases by immunohistochemistry, but in only 12.4% of nuclei in negative cases, indicating an association of HMGA2 expression and this chromosomal rearrangement (P=8.24 × 10(-10)). Four HMGA2-positive cases had greater than two HMGA2 hybridization signals per cell. No cases showed loss of a hybridization signal by interphase FISH for the frequently deleted region of 7q22 in uterine leiomyomata. One intravenous leiomyomatosis case analyzed by array comparative genomic hybridization revealed complex copy number variations. Finally, expression profiling was performed on three intravenous leiomyomatosis cases. Interestingly, hierarchical cluster analysis of the expression profiles revealed segregation of the intravenous leiomyomatosis cases with leiomyosarcoma rather than with myometrium, uterine leiomyoma of the usual histological type, or plexiform leiomyoma. These findings suggest that intravenous leiomyomatosis cases share some molecular cytogenetic characteristics with uterine leiomyoma, and expression profiles similar to that of leiomyosarcoma cases, further supporting their intermediate, quasi-malignant behavior.Modern Pathology advance online publication, 19 February 2016; doi:10.1038/modpathol.2016.36.
The transcription factor Wilms' tumor suppressor 1 (WT1) is key to podocyte development and viability; however, WT1 transcriptional networks in podocytes remain elusive. We provide a comprehensive analysis of the genome-wide WT1 transcriptional network in podocytes in vivo using chromatin immunoprecipitation followed by sequencing (ChIPseq) and RNA sequencing techniques. Our data show a specific role for WT1 in regulating the podocyte-specific transcriptome through binding to both promoters and enhancers of target genes. Furthermore, we inferred a podocyte transcription factor network consisting of WT1, LMX1B, TCF21, Fox-class and TEAD family transcription factors, and MAFB that uses tissue-specific enhancers to control podocyte gene expression. In addition to previously described WT1-dependent target genes, ChIPseq identified novel WT1-dependent signaling systems. These targets included components of the Hippo signaling system, underscoring the power of genome-wide transcriptional-network analyses. Together, our data elucidate a comprehensive gene regulatory network in podocytes suggesting that WT1 gene regulatory function and podocyte cell-type specification can best be understood in the context of transcription factor-regulatory element network interplay.
Dysregulation of ribosome biogenesis or translation can promote cancer, but the underlying mechanisms remain unclear. UTP18 is a component of the small subunit processome, a nucleolar multi-protein complex whose only known function is to cleave pre-ribosomal RNA to yield the 18S ribosomal RNA component of 40S ribosomal subunits. Here, we show that UTP18 also alters translation to promote stress resistance and growth, and that UTP18 is frequently gained and overexpressed in cancer. We observed that UTP18 localizes to the cytoplasm in a subset of cells, and that serum withdrawal increases cytoplasmic UTP18 localization. Cytoplasmic UTP18 associates with the translation complex and Hsp90 to upregulate the translation of IRES-containing transcripts such as HIF1a, Myc and VEGF, thereby inducing stress resistance. Hsp90 inhibition decreases cytoplasmic UTP18 and UTP18-induced increases in translation. Importantly, elevated UTP18 expression correlates with increased aggressiveness and decreased survival in numerous cancers. Enforced UTP18 overexpression promotes transformation and tumorigenesis, whereas UTP18 knockdown inhibits these processes. This stress adaptation mechanism is thus co-opted for growth by cancers, and its inhibition may represent a promising new therapeutic target.
Kidney injury molecule-1 (KIM-1)/T cell Ig and mucin domain-containing protein-1 (TIM-1) is upregulated more than other proteins after AKI, and it is highly expressed in renal damage of various etiologies. In this capacity, KIM-1/TIM-1 acts as a phosphatidylserine receptor on the surface of injured proximal tubular epithelial cells, mediating phagocytosis of apoptotic cells, and it may also act as a costimulatory molecule for immune cells. Despite recognition of KIM-1 as an important therapeutic target for kidney disease, the regulators of KIM-1 transcription in the kidney remain unknown. Using a bioinformatics approach, we identified upstream regulators of KIM-1 after AKI. In response to tubular injury in rat and human kidneys or oxidant stress in human proximal tubular epithelial cells (HPTECs), KIM-1 expression increased significantly in a manner that corresponded temporally and regionally with increased phosphorylation of checkpoint kinase 1 (Chk1) and STAT3. Both ischemic and oxidant stress resulted in a dramatic increase in reactive oxygen species that phosphorylated and activated Chk1, which subsequently bound to STAT3, phosphorylating it at S727. Furthermore, STAT3 bound to the KIM-1 promoter after ischemic and oxidant stress, and pharmacological or genetic induction of STAT3 in HPTECs increased KIM-1 mRNA and protein levels. Conversely, inhibition of STAT3 using siRNAs or dominant negative mutants reduced KIM-1 expression in a kidney cancer cell line (769-P) that expresses high basal levels of KIM-1. These observations highlight Chk1 and STAT3 as critical upstream regulators of KIM-1 expression after AKI and may suggest novel approaches for therapeutic intervention.
ChIP-seq has become the primary method for identifying in vivo protein-DNA interactions on a genome-wide scale, with nearly 800 publications involving the technique appearing in PubMed as of December 2012. Individually and in aggregate, these data are an important and information-rich resource. However, uncertainties about data quality confound their use by the wider research community. Recently, the Encyclopedia of DNA Elements (ENCODE) project developed and applied metrics to objectively measure ChIP-seq data quality. The ENCODE quality analysis was useful for flagging datasets for closer inspection, eliminating or replacing poor data, and for driving changes in experimental pipelines. There had been no similarly systematic quality analysis of the large and disparate body of published ChIP-seq profiles. Here, we report a uniform analysis of vertebrate transcription factor ChIP-seq datasets in the Gene Expression Omnibus (GEO) repository as of April 1, 2012. The majority (55%) of datasets scored as being highly successful, but a substantial minority (20%) were of apparently poor quality, and another ∼25% were of intermediate quality. We discuss how different uses of ChIP-seq data are affected by specific aspects of data quality, and we highlight exceptional instances for which the metric values should not be taken at face value. Unexpectedly, we discovered that a significant subset of control datasets (i.e., no immunoprecipitation and mock immunoprecipitation samples) display an enrichment structure similar to successful ChIP-seq data. This can, in turn, affect peak calling and data interpretation. Published datasets identified here as high-quality comprise a large group that users can draw on for large-scale integrated analysis. In the future, ChIP-seq quality assessment similar to that used here could guide experimentalists at early stages in a study, provide useful input in the publication process, and be used to stratify ChIP-seq data for different community-wide uses.
Increased levels of hypoxia and hypoxia-inducible factor 1α (HIF-1α) in human sarcomas correlate with tumor progression and radiation resistance. Prolonged antiangiogenic therapy of tumors not only delays tumor growth but may also increase hypoxia and HIF-1α activity. In our recent clinical trial, treatment with the vascular endothelial growth factor A (VEGF-A) antibody, bevacizumab, followed by a combination of bevacizumab and radiation led to near complete necrosis in nearly half of sarcomas. Gene Set Enrichment Analysis of microarrays from pretreatment biopsies found that the Gene Ontology category "Response to hypoxia" was upregulated in poor responders and that the hierarchical clustering based on 140 hypoxia-responsive genes reliably separated poor responders from good responders. The most commonly used chemotherapeutic drug for sarcomas, doxorubicin (Dox), was recently found to block HIF-1α binding to DNA at low metronomic doses. In four sarcoma cell lines, HIF-1α shRNA or Dox at low concentrations blocked HIF-1α induction of VEGF-A by 84-97% and carbonic anhydrase 9 by 83-93%. HT1080 sarcoma xenografts had increased hypoxia and/or HIF-1α activity with increasing tumor size and with anti-VEGF receptor antibody (DC101) treatment. Combining DC101 with HIF-1α shRNA or metronomic Dox had a synergistic effect in suppressing growth of HT1080 xenografts, at least in part via induction of tumor endothelial cell apoptosis. In conclusion, sarcomas respond to increased hypoxia by expressing HIF-1α target genes that may promote resistance to antiangiogenic and other therapies. HIF-1α inhibition blocks this evasive resistance and augments destruction of the tumor vasculature.
OBJECTIVE: Short bowel syndrome remains a condition of high morbidity and mortality, and current therapeutic options carry significant side effects. To identify new treatments we focused on postresection changes in microRNAs--short noncoding RNAs, which suppress target genes--and suggest a previously undiscovered role for microRNA-125a (mir-125a) in intestinal adaptation. METHODS: Rats underwent either 80% massive small bowel resection or transection and were harvested after 48 hours. Jejunum was harvested for microRNA microarrays, laser capture microdissection, and RNA and protein analysis. Mir-125a was overexpressed in intestinal epithelium-6 (crypt-derived) cells (IEC-6) and effects on proliferation and apoptosis determined using MTS and flow cytometry. Expression of potential targets of mir-125a in rat jejunum and IEC-6 cells was determined using quantitative real-time polymerase chain reaction (RNA) and Western blotting (protein). RESULTS: Resection upregulated mir-125a and mir-214 by 2.4-folds and 3.2-folds, respectively. Highest levels of expression were noted in the crypt fraction. Mir-125a overexpression induced apoptosis and resultant growth arrest in IEC-6 cells. The expression of the prosurvival Bcl-2 family member Mcl-1 was downregulated in both mir-125a-overexpressing IEC-6 cells and in jejunum of resected rats, confirming Mcl-1 as a previously undiscovered target of mir-125a. CONCLUSIONS: Upregulation of mir-125a suppresses the prosurvival protein Mcl1, producing the increase in apoptosis known to accompany the proliferative changes characteristic of intestinal adaptation. Our data highlight a potential role for microRNAs as mediators of the adaptive process and may facilitate the development of new therapeutic options for short bowel syndrome.
MicroRNAs (miRNAs) are endogenous noncoding RNA molecules that are involved in post-transcriptional gene silencing. Using global miRNA expression profiling, we found miR-21, -155, and 18a to be highly upregulated in rat kidneys following tubular injury induced by ischemia/reperfusion (I/R) or gentamicin administration. Mir-21 and -155 also showed decreased expression patterns in blood and urinary supernatants in both models of kidney injury. Furthermore, urinary levels of miR-21 increased 1.2-fold in patients with clinical diagnosis of acute kidney injury (AKI) (n = 22) as compared with healthy volunteers (n = 25) (p < 0.05), and miR-155 decreased 1.5-fold in patients with AKI (p < 0.01). We identified 29 messenger RNA core targets of these 3 miRNAs using the context likelihood of relatedness algorithm and found these predicted gene targets to be highly enriched for genes associated with apoptosis or cell proliferation. Taken together, these results suggest that miRNA-21 and -155 could potentially serve as translational biomarkers for detection of AKI and may play a critical role in the pathogenesis of kidney injury and tissue repair process.
PURPOSE: To facilitate the identification of genes associated with cataract and other ocular defects, the authors developed and validated a computational tool termed iSyTE (integrated Systems Tool for Eye gene discovery; http://bioinformatics.udel.edu/Research/iSyTE). iSyTE uses a mouse embryonic lens gene expression data set as a bioinformatics filter to select candidate genes from human or mouse genomic regions implicated in disease and to prioritize them for further mutational and functional analyses. METHODS: Microarray gene expression profiles were obtained for microdissected embryonic mouse lens at three key developmental time points in the transition from the embryonic day (E)10.5 stage of lens placode invagination to E12.5 lens primary fiber cell differentiation. Differentially regulated genes were identified by in silico comparison of lens gene expression profiles with those of whole embryo body (WB) lacking ocular tissue. RESULTS: Gene set analysis demonstrated that this strategy effectively removes highly expressed but nonspecific housekeeping genes from lens tissue expression profiles, allowing identification of less highly expressed lens disease-associated genes. Among 24 previously mapped human genomic intervals containing genes associated with isolated congenital cataract, the mutant gene is ranked within the top two iSyTE-selected candidates in approximately 88% of cases. Finally, in situ hybridization confirmed lens expression of several novel iSyTE-identified genes. CONCLUSIONS: iSyTE is a publicly available Web resource that can be used to prioritize candidate genes within mapped genomic intervals associated with congenital cataract for further investigation. Extension of this approach to other ocular tissue components will facilitate eye disease gene discovery.
Glioblastoma contains a hierarchy of stem-like cancer cells, but how this hierarchy is established is unclear. Here, we show that asymmetric Numb localization specifies glioblastoma stem-like cell (GSC) fate in a manner that does not require Notch inhibition. Numb is asymmetrically localized to CD133-hi GSCs. The predominant Numb isoform, Numb4, decreases Notch and promotes a CD133-hi, radial glial-like phenotype. However, upregulation of a novel Numb isoform, Numb4 delta 7 (Numb4d7), increases Notch and AKT activation while nevertheless maintaining CD133-hi fate specification. Numb knockdown increases Notch and promotes growth while favoring a CD133-lo, glial progenitor-like phenotype. We report the novel finding that Numb4 (but not Numb4d7) promotes SCF(Fbw7) ubiquitin ligase assembly and activation to increase Notch degradation. However, both Numb isoforms decrease epidermal growth factor receptor (EGFR) expression, thereby regulating GSC fate. Small molecule inhibition of EGFR activity phenocopies the effect of Numb on CD133 and Pax6. Clinically, homozygous NUMB deletions and low Numb mRNA expression occur primarily in a subgroup of proneural glioblastomas. Higher Numb expression is found in classical and mesenchymal glioblastomas and correlates with decreased survival. Thus, decreased Numb promotes glioblastoma growth, but the remaining Numb establishes a phenotypically diverse stem-like cell hierarchy that increases tumor aggressiveness and therapeutic resistance.
BACKGROUND: Recent withdrawals of major drugs have highlighted the critical importance of drug safety surveillance in the postmarketing phase. Limitations of spontaneous report data have led drug safety professionals to pursue alternative postmarketing surveillance approaches based on healthcare administrative claims data. These data are typically analysed by comparing the adverse event rates associated with a drug of interest to those of a single comparable reference drug. OBJECTIVE: The aim of this study was to determine whether adverse event detection can be improved by incorporating information from multiple reference drugs. We developed a pharmacological network model that implemented this approach and evaluated its performance. METHODS: We studied whether adverse event detection can be improved by incorporating information from multiple reference drugs, and describe two approaches for doing so. The first, reported previously, combines a set of related drugs into a single reference cohort. The second is a novel pharmacoepidemiological network model, which integrates multiple pair-wise comparisons across an entire set of related drugs into a unified consensus safety score for each drug. We also implemented a single reference drug approach for comparison with both multi-drug approaches. All approaches were applied within a sequential analysis framework, incorporating new information as it became available and addressing the issue of multiple testing over time. We evaluated all these approaches using statin (HMG-CoA reductase inhibitors) safety data from a large healthcare insurer in the US covering April 2000 through March 2005. RESULTS: We found that both multiple reference drug approaches offer earlier detection (6-13 months) than the single reference drug approach, without triggering additional false positives. CONCLUSIONS: Such combined approaches have the potential to be used with existing healthcare databases to improve the surveillance of therapeutics in the postmarketing phase over single-comparator methods. The proposed network approach also provides an integrated visualization framework enabling decision makers to understand the key high-level safety relationships amongst a group of related drugs.
Copy-number variation (CNV) is a major class of genomic variation with potentially important functional consequences in both normal and diseased populations. Remarkable advances in development of next-generation sequencing (NGS) platforms provide an unprecedented opportunity for accurate, high-resolution characterization of CNVs. In this unit, we give an overview of available computational tools for detection of CNVs and discuss comparative advantages and disadvantages of different approaches.
A transcriptional regulatory network represents a molecular framework in which developmental or environmental cues are transformed into differential expression of genes. Transcriptional regulation is mediated by the combinatorial interplay between cis-regulatory DNA elements and trans-acting transcription factors, and is perhaps the most important mechanism for controlling gene expression. Recent innovations, most notably the method for detecting protein-DNA interactions genome-wide, can help provide a comprehensive catalog of cis-regulatory elements and their interaction with given trans-acting factors in a given condition. A transcriptional regulatory network that integrates such information can lead to a systems-level understanding of regulatory mechanisms. In this review, we will highlight the key aspects of current knowledge on eukaryotic transcriptional regulation, especially on known transcription factors and their interacting regulatory elements. Then we will review some recent technical advances for genome-wide mapping of DNA-protein interactions based on high-throughput sequencing. Finally, we will discuss the types of biological insights that can be obtained from a network-level understanding of transcription regulation as well as future challenges in the field.
BACKGROUND: The molecular mechanisms that mediate the ischemia-reperfusion (I/R) injury in kidney are not completely understood. It is also largely unknown whether such mechanisms overlap with those governing the early development of kidney. METHODOLOGY/PRINCIPAL FINDINGS: We performed gene expression analysis to investigate the transcriptome changes during regeneration after I/R injury in the rat (0 hr, 6 hr, 24 hr, and 120 hr after reperfusion) and early development of mouse kidney (embryonic day 16 p.c. and postnatal 1 and 7 day). Pathway analysis revealed a wide spectrum of molecular functions that may participate in the regeneration and developmental processes of kidney as well as the functional association between them. While the genes associated with cell cycle, immunity, inflammation, and apoptosis were globally activated during the regeneration after I/R injury, the genes encoding various transporters and metabolic enzymes were down-regulated. We also observed that these injury-associated molecular functions largely overlap with those of early kidney development. In particular, the up-regulation of kinases and kinesins with roles in cell division was common during regeneration and early developmental kidney as validated by real-time PCR and immunohistochemistry. CONCLUSIONS: In addition to the candidate genes whose up-regulation constitutes an overlapping expression signature between kidney regeneration and development, this study lays a foundation for studying the functional relationship between two biological processes.