Shawn J. Cokus

Shotgun bisulphite sequencing of the Arabidopsis genome reveals DNA methylation patterning

Cytosine DNA methylation is important in regulating gene expression and in silencing transposons and other repetitive sequences. Recent genomic studies in Arabidopsis thaliana have revealed that many endogenous genes are methylated either within their promoters or within their transcribed regions, and that gene methylation is highly correlated with transcription levels. However, plants have different types of methylation controlled by different genetic pathways, and detailed information on the methylation status of each cytosine in any given genome is lacking. To this end, we generated a map at single-base-pair resolution of methylated cytosines for Arabidopsis, by combining bisulphite treatment of genomic DNA with ultra-high-throughput sequencing using the Illumina 1G Genome Analyser and Solexa sequencing technology. This approach, termed BS-Seq, unlike previous microarray-based methods, allows one to sensitively measure cytosine methylation on a genome-wide scale within specific sequence contexts. Here we describe methylation on previously inaccessible components of the genome and analyse the DNA methylation sequence composition and distribution. We also describe the effect of various DNA methylation mutants on genome-wide methylation patterns, and demonstrate that our newly developed library construction and computational methods can be applied to large genomes such as that of mouse.

Genome-wide High-Resolution Mapping and Functional Analysis of DNA Methylation in Arabidopsis

Cytosine methylation is important for transposon silencing and epigenetic regulation of endogenous genes, although the extent to which this DNA modification functions to regulate the genome is still unknown. Here we report the first comprehensive DNA methylation map of an entire genome, at 35 base pair resolution, using the flowering plant Arabidopsis thaliana as a model. We find that pericentromeric heterochromatin, repetitive sequences, and regions producing small interfering RNAs are heavily methylated. Unexpectedly, over one-third of expressed genes contain methylation within transcribed regions, whereas only approximately 5% of genes show methylation within promoter regions. Interestingly, genes methylated in transcribed regions are highly expressed and constitutively active, whereas promoter-methylated genes show a greater degree of tissue-specific expression. Whole-genome tiling-array transcriptional profiling of DNA methyltransferase null mutants identified hundreds of genes and intergenic noncoding RNAs with altered expression levels, many of which may be epigenetically controlled by DNA methylation.

Genome-wide erasure of DNA methylation in mouse primordial germ cells is affected by Aid deficiency

Epigenetic reprogramming including demethylation of DNA occurs in mammalian primordial germ cells (PGCs) and in early embryos, and is important for the erasure of imprints and epimutations, and the return to pluripotency. The extent of this reprogramming and its molecular mechanisms are poorly understood. We previously showed that the cytidine deaminases AID and APOBEC1 can deaminate 5-methylcytosine in vitro and in Escherichia coli, and in the mouse are expressed in tissues in which demethylation occurs. Here we profiled DNA methylation throughout the genome by unbiased bisulphite next generation sequencing in wild-type and AID-deficient mouse PGCs at embryonic day (E)13.5. Wild-type PGCs revealed marked genome-wide erasure of methylation to a level below that of methylation deficient (Np95-/-, also called Uhrf1-/-) embryonic stem cells, with female PGCs being less methylated than male ones. By contrast, AID-deficient PGCs were up to three times more methylated than wild-type ones; this substantial difference occurred throughout the genome, with introns, intergenic regions and transposons being relatively more methylated than exons. Relative hypermethylation in AID-deficient PGCs was confirmed by analysis of individual loci in the genome. Our results reveal that erasure of DNA methylation in the germ line is a global process, hence limiting the potential for transgenerational epigenetic inheritance. AID deficiency interferes with genome-wide erasure of DNA methylation patterns, indicating that AID has a critical function in epigenetic reprogramming and potentially in restricting the inheritance of epimutations in mammals.

Conservation and divergence of methylation patterning in plants and animals

Cytosine DNA methylation is a heritable epigenetic mark present in many eukaryotic organisms. Although DNA methylation likely has a conserved role in gene silencing, the levels and patterns of DNA methylation appear to vary drastically among different organisms. Here we used shotgun genomic bisulfite sequencing (BS-Seq) to compare DNA methylation in eight diverse plant and animal genomes. We found that patterns of methylation are very similar in flowering plants with methylated cytosines detected in all sequence contexts, whereas CG methylation predominates in animals. Vertebrates have methylation throughout the genome except for CpG islands. Gene body methylation is conserved with clear preference for exons in most organisms. Furthermore, genes appear to be the major target of methylation in Ciona and honey bee. Among the eight organisms, the green alga Chlamydomonas has the most unusual pattern of methylation, having non-CG methylation enriched in exons of genes rather than in repeats and transposons. In addition, the Dnmt1 cofactor Uhrf1 has a conserved function in maintaining CG methylation in both transposons and gene bodies in the mouse, Arabidopsis, and zebrafish genomes.

Relationship between nucleosome positioning and DNA methylation

Nucleosomes compact and regulate access to DNA in the nucleus, and are composed of approximately 147 bases of DNA wrapped around a histone octamer. Here we report a genome-wide nucleosome positioning analysis of Arabidopsis thaliana using massively parallel sequencing of mononucleosomes. By combining this data with profiles of DNA methylation at single base resolution, we identified 10-base periodicities in the DNA methylation status of nucleosome-bound DNA and found that nucleosomal DNA was more highly methylated than flanking DNA. These results indicate that nucleosome positioning influences DNA methylation patterning throughout the genome and that DNA methyltransferases preferentially target nucleosome-bound DNA. We also observed similar trends in human nucleosomal DNA, indicating that the relationships between nucleosomes and DNA methyltransferases are conserved. Finally, as has been observed in animals, nucleosomes were highly enriched on exons, and preferentially positioned at intron–exon and exon–intron boundaries. RNA polymerase II (Pol II) was also enriched on exons relative to introns, consistent with the hypothesis that nucleosome positioning regulates Pol II processivity. DNA methylation is also enriched on exons, consistent with the targeting of DNA methylation to nucleosomes, and suggesting a role for DNA methylation in exon definition.

Three Acyltransferases and Nitrogen-responsive Regulator Are Implicated in Nitrogen Starvation-induced Triacylglycerol Accumulation in Chlamydomonas

Background: Nitrogen-starvation and other stresses induce triacylglycerol (TAG) accumulation in algae, but the relevant enzymes and corresponding signal transduction pathways are unknown. Results: RNA-Seq and genetic analysis revealed three acyltransferases that contribute to TAG accumulation. Conclusion: TAG synthesis results from recycling of membrane lipids and also by acylation of DAG. Significance: The genes are potential targets for manipulating TAG hyperaccumulation. Algae have recently gained attention as a potential source for biodiesel; however, much is still unknown about the biological triggers that cause the production of triacylglycerols. We used RNA-Seq as a tool for discovering genes responsible for triacylglycerol (TAG) production in Chlamydomonas and for the regulatory components that activate the pathway. Three genes encoding acyltransferases, DGAT1, DGTT1, and PDAT1, are induced by nitrogen starvation and are likely to have a role in TAG accumulation based on their patterns of expression. DGAT1 and DGTT1 also show increased mRNA abundance in other TAG-accumulating conditions (minus sulfur, minus phosphorus, minus zinc, and minus iron). Insertional mutants, pdat1-1 and pdat1-2, accumulate 25% less TAG compared with the parent strain, CC-4425, which demonstrates the relevance of the trans-acylation pathway in Chlamydomonas. The biochemical functions of DGTT1 and PDAT1 were validated by rescue of oleic acid sensitivity and restoration of TAG accumulation in a yeast strain lacking all acyltransferase activity. Time course analyses suggest than a SQUAMOSA promoter-binding protein domain transcription factor, whose mRNA increases precede that of lipid biosynthesis genes like DGAT1, is a candidate regulator of the nitrogen deficiency responses. An insertional mutant, nrr1-1, accumulates only 50% of the TAG compared with the parental strain in nitrogen-starvation conditions and is unaffected by other nutrient stresses, suggesting the specificity of this regulator for nitrogen-deprivation conditions.

Genome-Wide Association of Histone H3 Lysine Nine Methylation with CHG DNA Methylation in Arabidopsis thaliana

Methylation of histone H3 lysine 9 (H3K9) is a hallmark of transcriptional silencing in many organisms. In Arabidopsis thaliana, dimethylation of H3K9 (H3K9m2) is important in the silencing of transposons and in the control of DNA methylation. We constructed a high-resolution genome-wide map of H3K9m2 methylation by using chromatin immunoprecipitation coupled with whole genome Roche Nimblegen microarrays (ChIP-chip). We observed a very high coincidence between H3K9m2 and CHG methylation (where H is either A,T or C) throughout the genome. The coding regions of genes that are associated exclusively with methylation in a CG context did not contain H3K9m2. In addition, we observed two distinct patterns of H3K9m2. Transposons and other repeat elements present in the euchromatic arms contained small islands of H3K9m2 present at relatively low levels. In contrast, pericentromeric/centromeric regions of Arabidopsis chromosomes contained long, rarely interrupted blocks of H3K9m2 present at much higher average levels than seen in the chromosome arms. These results suggest a complex interplay between H3K9m2 and different types of DNA methylation and suggest that distinct mechanisms control H3K9m2 in different compartments of the genome.

RNA-Seq Analysis of Sulfur-Deprived Chlamydomonas Cells Reveals Aspects of Acclimation Critical for Cell Survival

Sulfur deprivation of the unicellular alga Chlamydomonas reinhardtii triggers massive changes in the levels of transcripts associated with sulfate assimilation, the synthesis and turnover of sulfur-containing metabolites, and the remodeling of the photosynthetic apparatus and cell wall. These responses are critical for survival of the organism under sulfur deprivation conditions. The Chlamydomonas reinhardtii transcriptome was characterized from nutrient-replete and sulfur-depleted wild-type and snrk2.1 mutant cells. This mutant is null for the regulatory Ser-Thr kinase SNRK2.1, which is required for acclimation of the alga to sulfur deprivation. The transcriptome analyses used microarray hybridization and RNA-seq technology. Quantitative RT-PCR evaluation of the results obtained by these techniques showed that RNA-seq reports a larger dynamic range of expression levels than do microarray hybridizations. Transcripts responsive to sulfur deprivation included those encoding proteins involved in sulfur acquisition and assimilation, synthesis of sulfur-containing metabolites, Cys degradation, and sulfur recycling. Furthermore, we noted potential modifications of cellular structures during sulfur deprivation, including the cell wall and complexes associated with the photosynthetic apparatus. Moreover, the data suggest that sulfur-deprived cells accumulate proteins with fewer sulfur-containing amino acids. Most of the sulfur deprivation responses are controlled by the SNRK2.1 protein kinase. The snrk2.1 mutant exhibits a set of unique responses during both sulfur-replete and sulfur-depleted conditions that are not observed in wild-type cells; the inability of this mutant to acclimate to S deprivation probably leads to elevated levels of singlet oxygen and severe oxidative stress, which ultimately causes cell death. The transcriptome results for wild-type and mutant cells strongly suggest the occurrence of massive changes in cellular physiology and metabolism as cells become depleted for sulfur and reveal aspects of acclimation that are likely critical for cell survival.

BS Seeker: precise mapping for bisulfite sequencing

BackgroundBisulfite sequencing using next generation sequencers yields genome-wide measurements of DNA methylation at single nucleotide resolution. Traditional aligners are not designed for mapping bisulfite-treated reads, where the unmethylated Cs are converted to Ts. We have developed BS Seeker, an approach that converts the genome to a three-letter alphabet and uses Bowtie to align bisulfite-treated reads to a reference genome. It uses sequence tags to reduce mapping ambiguity. Post-processing of the alignments removes non-unique and low-quality mappings.ResultsWe tested our aligner on synthetic data, a bisulfite-converted Arabidopsis library, and human libraries generated from two different experimental protocols. We evaluated the performance of our approach and compared it to other bisulfite aligners. The results demonstrate that among the aligners tested, BS Seeker is more versatile and faster. When mapping to the human genome, BS Seeker generates alignments significantly faster than RMAP and BSMAP. Furthermore, BS Seeker is the only alignment tool that can explicitly account for tags which are generated by certain library construction protocols.ConclusionsBS Seeker provides fast and accurate mapping of bisulfite-converted reads. It can work with BS reads generated from the two different experimental protocols, and is able to efficiently map reads to large mammalian genomes. The Python program is freely available at http://pellegrini.mcdb.ucla.edu/BS_Seeker/BS_Seeker.html.

MORC Family ATPases Required for Heterochromatin Condensation and Gene Silencing

To Silence or Not to Silence Repressed genes commonly have methylated DNA, and/or covalent histone modifications associated with silent chromatin, and/or associated small interfering (si)RNAs. All three features are components of gene-silencing systems (see the Perspective by Jacob and Martienssen). In a screen for components of DNA methylation gene-silencing systems in the flowering plant, Moissiard et al. (p. 1448, published online 3 May) identified the genes AtMoRC1 and AtMORC6, which are homologs of the mouse Microrchidia1 gene. AtMORC1 and AtMORC6 are involved in silencing transposable elements and genes corresponding to DNA-methylated loci, and yet neither gene is required for maintenance of DNA methylation. Instead, AtMoRC1 and AtMORC6 are related to proteins that remodel chromatin superstructure, and they seem to control gene-silencing through the higher-order compaction of methylated and silent chromatin. Qian et al. (p. 1445) identified an Arabidopsis gene, IDM1 (increased DNA methylation 1), that is involved in regulating DNA methylation at loci enriched for repeats and multigene families containing highly homologous genes. IDM1 protects target genes from DNA silencing and recognizes both histone H3 and methylated DNA at target loci and is able to acetylate histone H3. A conserved family of adenosine triphosphatases predicted to catalyze alterations in chromosome superstructure is required for gene silencing. Transposable elements (TEs) and DNA repeats are commonly targeted by DNA and histone methylation to achieve epigenetic gene silencing. We isolated mutations in two Arabidopsis genes, AtMORC1 and AtMORC6, which cause derepression of DNA-methylated genes and TEs but no losses of DNA or histone methylation. AtMORC1 and AtMORC6 are members of the conserved Microrchidia (MORC) adenosine triphosphatase (ATPase) family, which are predicted to catalyze alterations in chromosome superstructure. The atmorc1 and atmorc6 mutants show decondensation of pericentromeric heterochromatin, increased interaction of pericentromeric regions with the rest of the genome, and transcriptional defects that are largely restricted to loci residing in pericentromeric regions. Knockdown of the single MORC homolog in Caenorhabditis elegans also impairs transgene silencing. We propose that the MORC ATPases are conserved regulators of gene silencing in eukaryotes.