Junko Mizushima-Sugano

Diverse transcriptional initiation revealed by fine, large‐scale mapping of mRNA start sites

Determination of the mRNA start site is the first step in identifying the promoter region, which is of key importance for transcriptional regulation of gene expression. The ‘oligo‐capping’ method enabled us to introduce a sequence tag to the first base of an mRNA by replacing the cap structure of the mRNA. Using cDNA libraries made from oligo‐capped mRNAs, we could identify the transcriptional start site of an individual mRNA just by sequencing the 5′‐end of the cDNA. The fine mapping of transcriptional start sites was performed for 5880 mRNAs in 276 human genes. Contrary to our expectations, the majority of the genes showed a diverse distribution of transcriptional start sites. They were distributed over 61.7 bp with a standard deviation of 19.5. Our finding may reflect the dynamic nature of transcriptional initiation events of human genes in vivo.

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.

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.

Selective inhibition of transcription of the E. coli tufB operon by guanosine-5′-diphosphate-3′-diphosphate

SummaryWe have studied the effect of guanosine-5′-diphosphate-3′-diphosphate (ppGpp) on the transcription of the E. coli tufB and recA operons in a cell-free system containing of purified RNA polymerase holoenzyme. The transcription of the tufB operon which is under stringent control, was markedly inhibited by 0.5 mM ppGpp, and the extent of this inhibition was found to be greatly influenced by the Mg2+ and K+ concentrations in the reaction mixture. Maximal inhibition was obtained in the presence of 2 mM Mg2+ and 80–120 mM K+, whereas at higher concentrations of Mg2+ or lower concentrations of K+, practically no inhibition was observed. In contrast, transcription of the recA operon which is not subject to stringent control, was little affected by ppGpp at any of Mg2+ and K+ concentrations tested. The nucleotide inhibited initiation of transcription of tufB, while the rate of RNA chain elongation was not greatly inhibited in the presence of ppGpp.

Sequencing and phasing cancer mutations in lung cancers using a long-read portable sequencer

Abstract Here, we employed cDNA amplicon sequencing using a long-read portable sequencer, MinION, to characterize various types of mutations in cancer-related genes, namely, EGFR, KRAS, NRAS and NF1. For homozygous SNVs, the precision and recall rates were 87.5% and 91.3%, respectively. For previously reported hotspot mutations, the precision and recall rates reached 100%. The precise junctions of EML4-ALK, CCDC6-RET and five other gene fusions were also detected. Taking advantages of long-read sequencing, we conducted phasing of EGFR mutations and elucidated the mutational allelic backgrounds of anti-tumor drug-sensitive and resistant mutations, which could provide useful information for selecting therapeutic approaches. In the H1975 cells, 72% of the reads harbored both L858R and T790M mutations, and 22% of the reads harbored neither mutation. To ensure that the clinical requirements can be met in potentially low cancer cell populations, we further conducted a serial dilution analysis of the template for EGFR mutations. Several percent of the mutant alleles could be detected depending on the yield and quality of the sequencing data. Finally, we characterized the mutation genotypes in eight clinical samples. This method could be a convenient long-read sequencing-based analytical approach and thus may change the current approaches used for cancer genome sequencing.

Intrinsic Promoter Activities of Primary DNA Sequences in the Human Genome

Abstract In order to understand an overview of promoter activities intrinsic to primary DNA sequences in the human genome within a particular cell type, we carried out systematic quantitative luciferase assays of DNA fragments corresponding to putative promoters for 472 human genes which are expressed in HEK (human embryonic kidney epithelial) 293 cells. We observed the promoter activities of them were distributed in a bimodal manner; putative promoters belonging to the first group (with strong promoter activities) were designated as P1 and the latter (with weak promoter activities) as P2. The frequencies of the TATA-boxes, the CpG islands, and the overall G + C-contents were significantly different between these two populations, indicating there are two separate groups of promoters. Interestingly, similar analysis using 251 randomly isolated genomic DNA fragments showed that P2-type promoter occasionally occurs within the human genome. Furthermore, 35 DNA fragments corresponding to putative promoters of non-protein-coding transcripts (ncRNAs) shared similar features with the P2 in both promoter activities and sequence compositions. At least, a part of ncRNAs, which have been massively identified by full-length cDNA projects with no functional relevance inferred, may have originated from those sporadic promoter activities of primary DNA sequences inherent to the human genome.