Fumihito Miura

A large-scale full-length cDNA analysis to explore the budding yeast transcriptome

We performed a large-scale cDNA analysis to explore the transcriptome of the budding yeast Saccharomyces cerevisiae. We sequenced two cDNA libraries, one from the cells exponentially growing in a minimal medium and the other from meiotic cells. Both libraries were generated by using a vector-capping method that allows the accurate mapping of transcription start sites (TSSs). Consequently, we identified 11,575 TSSs associated with 3,638 annotated genomic features, including 3,599 ORFs, to suggest that most yeast genes have two or more TSSs. In addition, we identified 45 previously undescribed introns, including those affecting current ORF annotations and those spliced alternatively. Furthermore, the analysis revealed 667 transcription units in the intergenic regions and transcripts derived from antisense strands of 367 known features. We also found that 348 ORFs carry TSSs in their 3′-halves to generate sense transcripts starting from inside the ORFs. These results indicate that the budding yeast transcriptome is considerably more complex than previously thought, and it shares many recently revealed characteristics with the transcriptomes of mammals and other higher eukaryotes. Thus, the genome-wide active transcription that generates novel classes of transcripts appears to be an intrinsic feature of the eukaryotic cells. The budding yeast will serve as a versatile model for the studies on these aspects of transcriptome, and the full-length cDNA clones can function as an invaluable resource in such studies.

Amplification-free whole-genome bisulfite sequencing by post-bisulfite adaptor tagging

DNA methylation plays a key role in epigenetic regulation of eukaryotic genomes. Hence the genome-wide distribution of 5-methylcytosine, or the methylome, has been attracting intense attention. In recent years, whole-genome bisulfite sequencing (WGBS) has enabled methylome analysis at single-base resolution. However, WGBS typically requires microgram quantities of DNA as well as global PCR amplification, thereby precluding its application to samples of limited amounts. This is presumably because bisulfite treatment of adaptor-tagged templates, which is inherent to current WGBS methods, leads to substantial DNA fragmentation. To circumvent the bisulfite-induced loss of intact sequencing templates, we conceived an alternative method termed Post-Bisulfite Adaptor Tagging (PBAT) wherein bisulfite treatment precedes adaptor tagging by two rounds of random primer extension. The PBAT method can generate a substantial number of unamplified reads from as little as subnanogram quantities of DNA. It requires only 100 ng of DNA for amplification-free WGBS of mammalian genomes. Thus, the PBAT method will enable various novel applications that would not otherwise be possible, thereby contributing to the rapidly growing field of epigenomics.

High-resolution DNA methylome analysis of primordial germ cells identifies gender-specific reprogramming in mice

Dynamic epigenetic reprogramming occurs during mammalian germ cell development, although the targets of this process, including DNA demethylation and de novo methylation, remain poorly understood. We performed genome-wide DNA methylation analysis in male and female mouse primordial germ cells at embryonic days 10.5, 13.5, and 16.5 by whole-genome shotgun bisulfite sequencing. Our high-resolution DNA methylome maps demonstrated gender-specific differences in CpG methylation at genome-wide and gene-specific levels during fetal germline progression. There was extensive intra- and intergenic hypomethylation with erasure of methylation marks at imprinted, X-linked, or germline-specific genes during gonadal sex determination and partial methylation at particular retrotransposons. Following global demethylation and sex determination, CpG sites switched to de novo methylation in males, but the X-linked genes appeared resistant to the wave of de novo methylation. Significant differential methylation at a subset of imprinted loci was identified in both genders, and non-CpG methylation occurred only in male gonocytes. Our data establish the basis for future studies on the role of epigenetic modifications in germline development and other biological processes.

Mouse Oocyte Methylomes at Base Resolution Reveal Genome-Wide Accumulation of Non-CpG Methylation and Role of DNA Methyltransferases

DNA methylation is an epigenetic modification that plays a crucial role in normal mammalian development, retrotransposon silencing, and cellular reprogramming. Although methylation mainly occurs on the cytosine in a CG site, non-CG methylation is prevalent in pluripotent stem cells, brain, and oocytes. We previously identified non-CG methylation in several CG-rich regions in mouse germinal vesicle oocytes (GVOs), but the overall distribution of non-CG methylation and the enzymes responsible for this modification are unknown. Using amplification-free whole-genome bisulfite sequencing, which can be used with minute amounts of DNA, we constructed the base-resolution methylome maps of GVOs, non-growing oocytes (NGOs), and mutant GVOs lacking the DNA methyltransferase Dnmt1, Dnmt3a, Dnmt3b, or Dnmt3L. We found that nearly two-thirds of all methylcytosines occur in a non-CG context in GVOs. The distribution of non-CG methylation closely resembled that of CG methylation throughout the genome and showed clear enrichment in gene bodies. Compared to NGOs, GVOs were over four times more methylated at non-CG sites, indicating that non-CG methylation accumulates during oocyte growth. Lack of Dnmt3a or Dnmt3L resulted in a global reduction in both CG and non-CG methylation, showing that non-CG methylation depends on the Dnmt3a-Dnmt3L complex. Dnmt3b was dispensable. Of note, lack of Dnmt1 resulted in a slight decrease in CG methylation, suggesting that this maintenance enzyme plays a role in non-dividing oocytes. Dnmt1 may act on CG sites that remain hemimethylated in the de novo methylation process. Our results provide a basis for understanding the mechanisms and significance of non-CG methylation in mammalian oocytes.

Translation Initiation: A Regulatory Role for poly(A) Tracts in Front of the AUG Codon in Saccharomyces cerevisiae

The 5′-UTR serves as the loading dock for ribosomes during translation initiation and is the key site for translation regulation. Many genes in the yeast Saccharomyces cerevisiae contain poly(A) tracts in their 5′-UTRs. We studied these pre-AUG poly(A) tracts in a set of 3274 recently identified 5′-UTRs in the yeast to characterize their effect on in vivo protein abundance, ribosomal density, and protein synthesis rate in the yeast. The protein abundance and the protein synthesis rate increase with the length of the poly(A), but exhibit a dramatic decrease when the poly(A) length is ≥12. The ribosomal density also reaches the lowest level when the poly(A) length is ≥12. This supports the hypothesis that a pre-AUG poly(A) tract can bind to translation initiation factors to enhance translation initiation, but a long (≥12) pre-AUG poly(A) tract will bind to Pab1p, whose binding size is 12 consecutive A residues in yeast, resulting in repression of translation. The hypothesis explains why a long pre-AUG poly(A) leads to more efficient translation initiation than a short one when PABP is absent, and why pre-AUG poly(A) is short in the early genes but long in the late genes of vaccinia virus.

Highly sensitive targeted methylome sequencing by post-bisulfite adaptor tagging

The current gold standard method for methylome analysis is whole-genome bisulfite sequencing (WGBS), but its cost is substantial, especially for the purpose of multi-sample comparison of large methylomes. Shotgun bisulfite sequencing of target-enriched DNA, or targeted methylome sequencing (TMS), can be a flexible, cost-effective alternative to WGBS. However, the current TMS protocol requires a considerable amount of input DNA and hence is hardly applicable to samples of limited quantity. Here we report a method to overcome this limitation by using post-bisulfite adaptor tagging (PBAT), in which adaptor tagging is conducted after bisulfite treatment to circumvent bisulfite-induced loss of intact sequencing templates, thereby enabling TMS of a 100-fold smaller amount of input DNA with far fewer cycles of polymerase chain reaction than in the current protocol. We thus expect that the PBAT-mediated TMS will serve as an invaluable method in epigenomics.

DNA methylation and gene expression in Mimulus guttatus.

BackgroundThe presence of methyl groups on cytosine nucleotides across an organism’s genome (methylation) is a major regulator of genome stability, crossing over, and gene regulation. The capacity for DNA methylation to be altered by environmental conditions, and potentially passed between generations, makes it a prime candidate for transgenerational epigenetic inheritance. Here we conduct the first analysis of the Mimulus guttatus methylome, with a focus on the relationship between DNA methylation and gene expression.ResultsWe present a whole genome methylome for the inbred line Iron Mountain 62 (IM62). DNA methylation varies across chromosomes, genomic regions, and genes. We develop a model that predicts gene expression based on DNA methylation (R2 = 0.2). Post hoc analysis of this model confirms prior relationships, and identifies novel relationships between methylation and gene expression. Additionally, we find that DNA methylation is significantly depleted near gene transcriptional start sites, which may explain the recently discovered elevated rate of recombination in these same regions.ConclusionsThe establishment here of a reference methylome will be a useful resource for the continued advancement of M. guttatus as a model system. Using a model-based approach, we demonstrate that methylation patterns are an important predictor of variation in gene expression. This model provides a novel approach for differential methylation analysis that generates distinct and testable hypotheses regarding gene expression.

Differential display analysis of mutants for the transcription factor Pdr1p regulating multidrug resistance in the budding yeast

The transcription factor Pdr1p recognizes Pdr1p/Pdr3p‐response element (PDRE) to activate genes involved in multidrug resistance of the budding yeast. To identify novel targets of Pdr1p, we compared transcriptomes among the yeast cells bearing wild, disrupted and gain‐of‐function alleles of PDR1 using a high‐throughput fluorescent differential display PCR. Consequently, we identified 20 transcripts apparently regulated by Pdr1p, which are derived from well‐known target genes as well as those that have never been described in the context of drug resistance. Intriguingly, among the latter, a previously unrecognized gene bearing a small putative open reading frame preceded by a functional PDRE was found.

Unexpected Complexity of the Budding Yeast Transcriptome

The genome of the budding yeast Saccharomyces cerevisiae was sequenced over a decade ago and has been annotated to encode ∼6,000 genes. However, recent high throughput studies using tiling array hybridization and cDNA sequencing have revealed an unexpectedly large number of previously undescribed transcripts. They largely lack protein‐coding capacity and are transcribed from both strands of intragenic and intergenic regions in the genome. Accordingly, pervasive transcription leading to a plethora of noncoding RNAs, which was first revealed for mammalian genomes to attract intense attentions, is likely an intrinsic feature of eukaryotic genomes. Although it is not clear what fraction of these transcription events are functional, some were shown to induce transcriptional interference or histone modifications to regulate gene expression. The budding yeast may serve as an excellent model to study pervasive transcription and noncoding RNAs.

Software updates in the Illumina HiSeq platform affect whole-genome bisulfite sequencing.

BackgroundMethylation of cytosine in genomic DNA is a well-characterized epigenetic modification involved in many cellular processes and diseases. Whole-genome bisulfite sequencing (WGBS), such as MethylC-seq and post-bisulfite adaptor tagging sequencing (PBAT-seq), uses the power of high-throughput DNA sequencers and provides genome-wide DNA methylation profiles at single-base resolution. However, the accuracy and consistency of WGBS outputs in relation to the operating conditions of high-throughput sequencers have not been explored.ResultsWe have used the Illumina HiSeq platform for our PBAT-based WGBS, and found that different versions of HiSeq Control Software (HCS) and Real-Time Analysis (RTA) installed on the system provided different global CpG methylation levels (approximately 5% overall difference) for the same libraries. This problem was reproduced multiple times with different WGBS libraries and likely to be associated with the low sequence diversity of bisulfite-converted DNA. We found that HCS was the major determinant in the observed differences. To determine which version of HCS is most suitable for WGBS, we used substrates with predetermined CpG methylation levels, and found that HCS v2.0.5 is the best among the examined versions. HCS v2.0.12 showed the poorest performance and provided artificially lower CpG methylation levels when 5-methylcytosine is read as guanine (first read of PBAT-seq and second read of MethylC-seq). In addition, paired-end sequencing of low diversity libraries using HCS v2.2.38 or the latest HCS v2.2.58 was greatly affected by cluster densities.ConclusionsSoftware updates in the Illumina HiSeq platform can affect the outputs from low-diversity sequencing libraries such as WGBS libraries. More recent versions are not necessarily the better, and HCS v2.0.5 is currently the best for WGBS among the examined HCS versions. Thus, together with other experimental conditions, special care has to be taken on this point when CpG methylation levels are to be compared between different samples by WGBS.