BS Seeker: precise mapping for bisulfite sequencing

<p>Abstract</p> <p>Background</p> <p>Bisulfite 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.</p> <p>Results</p> <p>We tested our aligner on synthetic data, a bisulfite-converted <it>Arabidopsis </it>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.</p> <p>Conclusions</p> <p>BS 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 <url>http://pellegrini.mcdb.ucla.edu/BS_Seeker/BS_Seeker.html</url>.</p

Tissue-Specific Transcriptomes Reveal Gene Expression Trajectories in Two Maturing Skin Epithelial Layers in Zebrafish Embryos.

Epithelial cells are the building blocks of many organs, including skin. The vertebrate skin initially consists of two epithelial layers, the outer periderm and inner basal cell layers, which have distinct properties, functions, and fates. The embryonic periderm ultimately disappears during development, whereas basal cells proliferate to form the mature, stratified epidermis. Although much is known about mechanisms of homeostasis in mature skin, relatively little is known about the two cell types in pre-stratification skin. To define the similarities and distinctions between periderm and basal skin epithelial cells, we purified them from zebrafish at early development stages and deeply profiled their gene expression. These analyses identified groups of genes whose tissue enrichment changed at each stage, defining gene flow dynamics of maturing vertebrate epithelia. At each of 52 and 72 hr post-fertilization (hpf), more than 60% of genes enriched in skin cells were similarly expressed in both layers, indicating that they were common epithelial genes, but many others were enriched in one layer or the other. Both expected and novel genes were enriched in periderm and basal cell layers. Genes encoding extracellular matrix, junctional, cytoskeletal, and signaling proteins were prominent among those distinguishing the two epithelial cell types. In situ hybridization and BAC transgenes confirmed our expression data and provided new tools to study zebrafish skin. Collectively, these data provide a resource for studying common and distinguishing features of maturing epithelia

ATXR5 and ATXR6 are H3K27 monomethyltransferases required for chromatin structure and gene silencing.

Constitutive heterochromatin in Arabidopsis thaliana is marked by repressive chromatin modifications, including DNA methylation, histone H3 dimethylation at Lys9 (H3K9me2) and monomethylation at Lys27 (H3K27me1). The enzymes catalyzing DNA methylation and H3K9me2 have been identified; alterations in these proteins lead to reactivation of silenced heterochromatic elements. The enzymes responsible for heterochromatic H3K27me1, in contrast, remain unknown. Here we show that the divergent SET-domain proteins ARABIDOPSIS TRITHORAX-RELATED PROTEIN 5 (ATXR5) and ATXR6 have H3K27 monomethyltransferase activity, and atxr5 atxr6 double mutants have reduced H3K27me1 in vivo and show partial heterochromatin decondensation. Mutations in atxr5 and atxr6 also lead to transcriptional activation of repressed heterochromatic elements. Notably, H3K9me2 and DNA methylation are unaffected in double mutants. These results indicate that ATXR5 and ATXR6 form a new class of H3K27 methyltransferases and that H3K27me1 represents a previously uncharacterized pathway required for transcriptional repression in Arabidopsis

Detecting coordinated regulation of multi-protein complexes using logic analysis of gene expression

<p>Abstract</p> <p>Background</p> <p>Many of the functional units in cells are multi-protein complexes such as RNA polymerase, the ribosome, and the proteasome. For such units to work together, one might expect a high level of regulation to enable co-appearance or repression of sets of complexes at the required time. However, this type of coordinated regulation between whole complexes is difficult to detect by existing methods for analyzing mRNA co-expression. We propose a new methodology that is able to detect such higher order relationships.</p> <p>Results</p> <p>We detect coordinated regulation of multiple protein complexes using <it>logic analysis </it>of gene expression data. Specifically, we identify gene triplets composed of genes whose expression profiles are found to be related by various types of logic functions. In order to focus on complexes, we associate the members of a gene triplet with the distinct protein complexes to which they belong. In this way, we identify complexes related by specific kinds of regulatory relationships. For example, we may find that the transcription of complex C is increased only if the transcription of both complex A AND complex B is repressed. We identify hundreds of examples of coordinated regulation among complexes under various stress conditions. Many of these examples involve the ribosome. Some of our examples have been previously identified in the literature, while others are novel. One notable example is the relationship between the transcription of the ribosome, RNA polymerase and mannosyltransferase II, which is involved in N-linked glycan processing in the Golgi.</p> <p>Conclusions</p> <p>The analysis proposed here focuses on relationships among triplets of genes that are not evident when genes are examined in a pairwise fashion as in typical clustering methods. By grouping gene triplets, we are able to decipher coordinated regulation among sets of three complexes. Moreover, using all triplets that involve coordinated regulation with the ribosome, we derive a large network involving this essential cellular complex. In this network we find that all multi-protein complexes that belong to the same functional class are regulated in the same direction as a group (either induced or repressed).</p

Genome-wide analysis of mono-, di- and trimethylation of histone H3 lysine 4 in Arabidopsis thaliana

Analysis of the genome-wide distribution patterns of histone H3 lysine4 methylation in Arabidopsis thaliana seedlings shows that it has widespread roles in regulating gene expression

Modeling the regulatory network of histone acetylation in Saccharomyces cerevisiae

Acetylation of histones plays an important role in regulating transcription. Histone acetylation is mediated partly by the recruitment of specific histone acetyltransferases (HATs) and deacetylases (HDACs) to genomic loci by transcription factors, resulting in modulation of gene expression. Although several specific interactions between transcription factors and HATs and HDACs have been elaborated in Saccharomyces cerevisiae, the full regulatory network remains uncharacterized. We have utilized a linear regression of optimized sigmoidal functions to correlate transcription factor binding patterns to the acetylation profiles of 11 lysines in the four core histones measured at all S. cerevisiae promoters. The resulting associations are combined with large-scale protein–protein interaction data sets to generate a comprehensive model that relates recruitment of specific HDACs and HATs to transcription factors and their target genes and the resulting effects on individual lysines. This model provides a broad and detailed view of the regulatory network, describing which transcription factors are most significant in regulating acetylation of specific lysines at defined promoters. We validate the model, both computationally and experimentally, to demonstrate that it yields accurate predictions of these regulatory mechanisms

Propionibacterium acnes bacteriophages display limited genetic diversity and broad killing activity against bacterial skin isolates.

UnlabelledInvestigation of the human microbiome has revealed diverse and complex microbial communities at distinct anatomic sites. The microbiome of the human sebaceous follicle provides a tractable model in which to study its dominant bacterial inhabitant, Propionibacterium acnes, which is thought to contribute to the pathogenesis of the human disease acne. To explore the diversity of the bacteriophages that infect P. acnes, 11 P. acnes phages were isolated from the sebaceous follicles of donors with healthy skin or acne and their genomes were sequenced. Comparative genomic analysis of the P. acnes phage population, which spans a 30-year temporal period and a broad geographic range, reveals striking similarity in terms of genome length, percent GC content, nucleotide identity (&gt;85%), and gene content. This was unexpected, given the far-ranging diversity observed in virtually all other phage populations. Although the P. acnes phages display a broad host range against clinical isolates of P. acnes, two bacterial isolates were resistant to many of these phages. Moreover, the patterns of phage resistance correlate closely with the presence of clustered regularly interspaced short palindromic repeat elements in the bacteria that target a specific subset of phages, conferring a system of prokaryotic innate immunity. The limited diversity of the P. acnes bacteriophages, which may relate to the unique evolutionary constraints imposed by the lipid-rich anaerobic environment in which their bacterial hosts reside, points to the potential utility of phage-based antimicrobial therapy for acne.ImportancePropionibacterium acnes is a dominant member of the skin microflora and has also been implicated in the pathogenesis of acne; however, little is known about the bacteriophages that coexist with and infect this bacterium. Here we present the novel genome sequences of 11 P. acnes phages, thereby substantially increasing the amount of available genomic information about this phage population. Surprisingly, we find that, unlike other well-studied bacteriophages, P. acnes phages are highly homogeneous and show a striking lack of genetic diversity, which is perhaps related to their unique and restricted habitat. They also share a broad ability to kill clinical isolates of P. acnes; phage resistance is not prevalent, but when detected, it appears to be conferred by chromosomally encoded immunity elements within the host genome. We believe that these phages display numerous features that would make them ideal candidates for the development of a phage-based therapy for acne

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

MTHFD1 controls DNA methylation in Arabidopsis.

DNA methylation is an epigenetic mechanism that has important functions in transcriptional silencing and is associated with repressive histone methylation (H3K9me). To further investigate silencing mechanisms, we screened a mutagenized Arabidopsis thaliana population for expression of SDCpro-GFP, redundantly controlled by DNA methyltransferases DRM2 and CMT3. Here, we identify the hypomorphic mutant mthfd1-1, carrying a mutation (R175Q) in the cytoplasmic bifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase (MTHFD1). Decreased levels of oxidized tetrahydrofolates in mthfd1-1 and lethality of loss-of-function demonstrate the essential enzymatic role of MTHFD1 in Arabidopsis. Accumulation of homocysteine and S-adenosylhomocysteine, genome-wide DNA hypomethylation, loss of H3K9me and transposon derepression indicate that S-adenosylmethionine-dependent transmethylation is inhibited in mthfd1-1. Comparative analysis of DNA methylation revealed that the CMT3 and CMT2 pathways involving positive feedback with H3K9me are mostly affected. Our work highlights the sensitivity of epigenetic networks to one-carbon metabolism due to their common S-adenosylmethionine-dependent transmethylation and has implications for human MTHFD1-associated diseases