Kousuke Tanimoto

Massive transcriptional start site analysis of human genes in hypoxia cells

Combining our full-length cDNA method and the massively parallel sequencing technology, we developed a simple method to collect precise positional information of transcriptional start sites (TSSs) together with digital information of the gene-expression levels in a high throughput manner. We applied this method to observe gene-expression changes in a colon cancer cell line cultured in normoxic and hypoxic conditions. We generated more than 100 million 36-base TSS-tag sequences and revealed comprehensive features of hypoxia responsive alterations in the transcriptional landscape of the human genome. The features include presence of inducible ‘hot regions’ in 54 genomic regions, 220 novel hypoxia inducible promoters that may drive non-protein-coding transcripts, 191 hypoxia responsive alternative promoters and detailed views of 120 novel as well as known hypoxia responsive genes. We further analyzed hypoxic response of different cells using additional 60 million TSS-tags and found that the degree of the gene-expression changes were different among cell lines, possibly reflecting cellular robustness against hypoxia. The novel dynamic figure of the human gene transcriptome will deepen our understanding of the transcriptional program of the human genome as well as bringing new insights into the biology of cancer cells in hypoxia.

Genome-wide characterization of transcriptional start sites in humans by integrative transcriptome analysis

We performed a genome-wide analysis of transcriptional start sites (TSSs) in human genes by multifaceted use of a massively parallel sequencer. By analyzing 800 million sequences that were obtained from various types of transcriptome analyses, we characterized 140 million TSS tags in 12 human cell types. Despite the large number of TSS clusters (TSCs), the number of TSCs was observed to decrease sharply with increasing expression levels. Highly expressed TSCs exhibited several characteristic features: Nucleosome-seq analysis revealed highly ordered nucleosome structures, ChIP-seq analysis detected clear RNA polymerase II binding signals in their surrounding regions, evaluations of previously sequenced and newly shotgun-sequenced complete cDNA sequences showed that they encode preferable transcripts for protein translation, and RNA-seq analysis of polysome-incorporated RNAs yielded direct evidence that those transcripts are actually translated into proteins. We also demonstrate that integrative interpretation of transcriptome data is essential for the selection of putative alternative promoter TSCs, two of which also have protein consequences. Furthermore, discriminative chromatin features that separate TSCs at different expression levels were found for both genic TSCs and intergenic TSCs. The collected integrative information should provide a useful basis for future biological characterization of TSCs.

The Tumor-Suppressive miR-497-195 Cluster Targets Multiple Cell-Cycle Regulators in Hepatocellular Carcinoma

MicroRNAs (miRNAs) are key post-transcriptional regulators of gene expression and commonly deregulated in carcinogenesis. To explore functionally crucial tumor-suppressive (TS)-miRNAs in hepatocellular carcinoma (HCC), we performed integrative function- and expression-based screenings of TS-miRNAs in six HCC cell lines. The screenings identified seven miRNAs, which showed growth-suppressive activities through the overexpression of each miRNA and were endogenously downregulated in HCC cell lines. Further expression analyses using a large panel of HCC cell lines and primary tumors demonstrated four miRNAs, miR-101, -195, -378 and -497, as candidate TS-miRNAs frequently silenced in HCCs. Among them, two clustered miRNAs miR-195 and miR-497 showed significant growth-suppressive activity with induction of G1 arrest. Comprehensive exploration of their targets using Argonute2-immunoprecipitation-deep-sequencing (Ago2-IP-seq) and genome-wide expression profiling after their overexpression followed by pathway analysis, revealed a significant enrichment of cell cycle regulators. Among the candidates, we successfully identified CCNE1, CDC25A, CCND3, CDK4, and BTRC as direct targets for miR-497 and miR-195. Moreover, target genes frequently upregulated in HCC in a tumor-specific manner, such as CDK6, CCNE1, CDC25A and CDK4, showed an inverse correlation in the expression of miR-195 and miR-497, and their targets. These results suggest the molecular pathway regulating cell cycle progression to be integrally altered by downregulation of miR-195 and miR-497 expression, leading to the aberrant cell proliferation in hepatocarcinogenesis.

Genome-wide identification and annotation of HIF-1α binding sites in two cell lines using massively parallel sequencing

We identified 531 and 616 putative HIF-1α target sites by ChIP-Seq in the cancerous cell line DLD-1 and the non-cancerous cell line TIG-3, respectively. We also examined the positions and expression levels of transcriptional start sites (TSSs) in these cell lines using our TSS-Seq method. We observed that 121 and 48 genes in DLD-1 and TIG-3 cells, respectively, had HIF-1α binding sites in proximal regions of the previously reported TSSs that were up-regulated at the transcriptional level. In addition, 193 and 123 of the HIF-1α target sites, respectively, were located in proximal regions of previously uncharacterized TSSs, namely, TSSs of putative alternative promoters of protein-coding genes or promoters of putative non-protein-coding transcripts. The hypoxic response of DLD-1 cells was more significant than that of TIG-3 cells with respect to both the number of target sites and the degree of induced changes in transcript expression. The Nucleosome-Seq and ChIP-Seq analyses of histone modifications revealed that the chromatin formed an open structure in regions surrounding the HIF-1α binding sites, but this event occurred prior to the actual binding of HIF-1α. Different cellular histories may be encoded by chromatin structures and determine the activation of specific genes in response to hypoxic shock.

Screening for possible miRNA–mRNA associations in a colon cancer cell line

MicroRNAs (miRNAs) are small non-coding RNAs mediating the regulation of gene expression in various biological contexts, including carcinogenesis. Here, we screened putative associations between 34, 45, and 103 miRNAs and 164, 391, and 81 mRNAs via Argonaute1 (Ago1) or Ago2 immunoprecipitation (IP) experiments in a colon cancer cell line. We used a combination of RIP Seq analysis. RNAs that were co-immunoprecipitated with Ago1 or Ago2 were used for massively parallel small RNA and mRNA sequencing. The detected miRNAs and mRNAs were further associated with one another based on in silico target predictions. Analysis of the putative associations indicated that, although Ago1 and Ago2 shared a similar repertory of miRNAs, the mRNAs possibly regulated by those miRNAs seemed different. The mRNAs detected with Ago1 IP were indicated to be frequently associated with genes having constitutive cellular functions, regulated by a smaller number of miRNAs, and appeared to receive more stringent translational regulation. In contrast, putative miRNA-mRNA associations detected with Ago2 IP appeared to be related to signal transduction genes, which had a larger number of possible miRNA binding sites. We then conducted a similar analysis using the colon cancer cells cultured under hypoxia and identified potential hypoxia-induced miRNA-mRNA associations, which included several well-characterized cancer-related genes as novel putative miRNA targets.

MiR-634 activates the mitochondrial apoptosis pathway and enhances chemotherapy-induced cytotoxicity

Some tumor-suppressing miRNAs target multiple oncogenes concurrently and therefore may be useful as cancer therapeutic agents. Further, such miRNAs may be useful to address chemotherapeutic resistance in cancer, which remains a primary clinical challenge in need of solutions. Thus, cytoprotective processes upregulated in cancer cells that are resistant to chemotherapy are a logical target for investigation. Here, we report that overexpression of miR-634 activates the mitochondrial apoptotic pathway by direct concurrent targeting of genes associated with mitochondrial homeostasis, antiapoptosis, antioxidant ability, and autophagy. In particular, we show how enforced expression of miR-634 enhanced chemotherapy-induced cytotoxicity in a model of esophageal squamous cell carcinoma, where resistance to chemotherapy remains clinically problematic. Our findings illustrate how reversing miR-634-mediated cytoprotective processes may offer a broadly useful approach to improving cancer therapy.

Characterization of STAT6 Target Genes in Human B Cells and Lung Epithelial Cells

Using ChIP Seq, we identified 556 and 467 putative STAT6 target sites in the Burkitts lymphoma cell line Ramos and in the normal lung epithelial cell line BEAS2B, respectively. We also examined the positions and expression of transcriptional start sites (TSSs) in these cells using our TSS Seq method. We observed that 44 and 132 genes in Ramos and BEAS2B, respectively, had STAT6 binding sites in proximal regions of their previously reported TSSs that were up-regulated at the transcriptional level. In addition, 406 and 109 of the STAT6 target sites in Ramos and BEAS2B, respectively, were located in proximal regions of previously uncharacterized TSSs. The target genes identified in Ramos and BEAS2B cells in this study and in Th2 cells in previous studies rarely overlapped and differed in their identity. Interestingly, ChIP Seq analyses of histone modifications and RNA polymerase II revealed that chromatin formed an active structure in regions surrounding the STAT6 binding sites; this event also frequently occurred in different cell types, although neither STAT6 binding nor TSS induction was observed. The rough landscape of STAT6-responsive sites was found to be shaped by chromatin structure, but distinct cellular responses were mainly mediated by distinct sets of transcription factors.

CREB is activated by ER stress and modulates the unfolded protein response by regulating the expression of IRE1α and PERK.

Living cells are frequently exposed to various stresses. Hypoxic conditions induce endoplasmic reticulum (ER) stress, and activate the unfolded protein response (UPR) to maintain homeostasis. We previously reported that CREB has an important role in the proper response to prolonged hypoxia. To further understand the role of CREB in the hypoxic response, CREB stable knock-down (CREB-KD) cells were established from breast cancer MDA-MB231 cells and analyzed. CREB was activated by ER stress, and activation of CREB and the UPR pathway occurred in a coordinated manner in response to different stimuli, including ER stress-inducing chemicals, prolonged hypoxia, and oxygen-glucose deprivation (OGD). Depletion of CREB decreased the expression of IRE1α and PERK, two critical UPR signaling molecules. Promoter analysis and a chromatin immunoprecipitation assay indicated that CREB binds to the promoter region of these genes and regulates their expression. ER stress induced by hypoxia was reduced in CREB-KD cells, leading to reduced tumor metastasis to the lung. Finally, OGD strongly activated the UPR and induced cell death in control cells, whereas the UPR was moderately activated in CREB-KD cells, which were more resistant to cell death. This study demonstrates a new role for CREB as a regulator of ER stress, which is required to properly respond to stressful conditions, such as hypoxia.

Simultaneous Detection of Both Single Nucleotide Variations and Copy Number Alterations by Next-Generation Sequencing in Gorlin Syndrome

Gorlin syndrome (GS) is an autosomal dominant disorder that predisposes affected individuals to developmental defects and tumorigenesis, and caused mainly by heterozygous germline PTCH1 mutations. Despite exhaustive analysis, PTCH1 mutations are often unidentifiable in some patients; the failure to detect mutations is presumably because of mutations occurred in other causative genes or outside of analyzed regions of PTCH1, or copy number alterations (CNAs). In this study, we subjected a cohort of GS-affected individuals from six unrelated families to next-generation sequencing (NGS) analysis for the combined screening of causative alterations in Hedgehog signaling pathway-related genes. Specific single nucleotide variations (SNVs) of PTCH1 causing inferred amino acid changes were identified in four families (seven affected individuals), whereas CNAs within or around PTCH1 were found in two families in whom possible causative SNVs were not detected. Through a targeted resequencing of all coding exons, as well as simultaneous evaluation of copy number status using the alignment map files obtained via NGS, we found that GS phenotypes could be explained by PTCH1 mutations or deletions in all affected patients. Because it is advisable to evaluate CNAs of candidate causative genes in point mutation-negative cases, NGS methodology appears to be useful for improving molecular diagnosis through the simultaneous detection of both SNVs and CNAs in the targeted genes/regions.

Comprehensive investigation of CASK mutations and other genetic etiologies in 41 patients with intellectual disability and microcephaly with pontine and cerebellar hypoplasia (MICPCH)

The CASK gene (Xp11.4) is highly expressed in the mammalian nervous system and plays several roles in neural development and synaptic function. Loss-of-function mutations of CASK are associated with intellectual disability and microcephaly with pontine and cerebellar hypoplasia (MICPCH), especially in females. Here, we present a comprehensive investigation of 41 MICPCH patients, analyzed by mutational search of CASK and screening of candidate genes using an SNP array, targeted resequencing and whole-exome sequencing (WES). In total, we identified causative or candidate genomic aberrations in 37 of the 41 cases (90.2%). CASK aberrations including a rare mosaic mutation in a male patient, were found in 32 cases, and a mutation in ITPR1, another known gene in which mutations are causative for MICPCH, was found in one case. We also found aberrations involving genes other than CASK, such as HDAC2, MARCKS, and possibly HS3ST5, which may be associated with MICPCH. Moreover, the targeted resequencing screening detected heterozygous variants in RELN in two cases, of uncertain pathogenicity, and WES analysis suggested that concurrent mutations of both DYNC1H1 and DCTN1 in one case could lead to MICPCH. Our results not only identified the etiology of MICPCH in nearly all the investigated patients but also suggest that MICPCH is a genetically heterogeneous condition, in which CASK inactivating mutations appear to account for the majority of cases.