Discovering transcriptional modules by Bayesian data integration.

MOTIVATION: We present a method for directly inferring transcriptional modules (TMs) by integrating gene expression and transcription factor binding (ChIP-chip) data. Our model extends a hierarchical Dirichlet process mixture model to allow data fusion on a gene-by-gene basis. This encodes the intuition that co-expression and co-regulation are not necessarily equivalent and hence we do not expect all genes to group similarly in both datasets. In particular, it allows us to identify the subset of genes that share the same structure of transcriptional modules in both datasets. RESULTS: We find that by working on a gene-by-gene basis, our model is able to extract clusters with greater functional coherence than existing methods. By combining gene expression and transcription factor binding (ChIP-chip) data in this way, we are better able to determine the groups of genes that are most likely to represent underlying TMs. AVAILABILITY: If interested in the code for the work presented in this article, please contact the authors. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Accurate sodium bisulfite sequencing in plants.

DNA cytosine methylation is a conserved epigenetic modification frequently correlating with transcriptional silencing in a wide variety of eukaryotic organisms. Sodium bisulfite treatment of DNA converts unmethylated cytosine to uracil, while 5-methylated cytosine is protected. We describe techniques that ensure reliable sequencing data following sodium bisulfite conversion and to avoid common pitfalls such as amplification of unconverted DNA and inclusion of sibling clones.

Copper homeostasis in Salmonella is atypical and copper-CueP is a major periplasmic metal complex.

Salmonella enterica sv. typhimurium (S. enterica sv. Typhimurium) has two metal-transporting P(1)-type ATPases whose actions largely overlap with respect to growth in elevated copper. Mutants lacking both ATPases over-accumulate copper relative to wild-type or either single mutant. Such duplication of ATPases is unusual in bacterial copper tolerance. Both ATPases are under the control of MerR family metal-responsive transcriptional activators. Analyses of periplasmic copper complexes identified copper-CueP as one of the predominant metal pools. Expression of cueP was recently shown to be controlled by the same metal-responsive activator as one of the P(1)-type ATPase genes (copA), and copper-CueP is a further atypical feature of copper homeostasis in S. enterica sv. Typhimurium. Elevated copper is detected by a reporter construct driven by the promoter of copA in wild-type S. enterica sv. Typhimurium during infection of macrophages. Double mutants missing both ATPases also show reduced survival inside cultured macrophages. It is hypothesized that elevated copper within macrophages may have selected for specialized copper-resistance systems in pathogenic microorganism such as S. enterica sv. Typhimurium.

Co-ordination of pathogenicity island expression by the BipA GTPase in enteropathogenic Escherichia coli (EPEC).

BipA is a novel member of the ribosome binding GTPase superfamily and is widely distributed in bacteria and plants. We report here that it regulates -multiple cell surface- and virulence-associated -components in the enteropathogenic Escherichia coli (EPEC) strain E2348/69. The regulated components include bacterial flagella, the espC pathogenicity island and a type III secretion system specified by the locus of enterocyte effacement (LEE). BipA positively regulated the espC and LEE gene clusters through transcriptional control of the LEE-encoded regulator, Ler. Additionally, it affected the pattern of proteolysis of intimin, a key LEE-encoded adhesin specified by the LEE. BipA control of the LEE operated independently of the previously characterized regulators Per, integration host factor and H-NS. In contrast, it negatively regulated the flagella-mediated motility of EPEC and in a Ler-independent manner. Our results indicate that the BipA GTPase functions high up in diverse regulatory cascades to co-ordinate the expression of key pathogenicity islands and other virulence-associated factors in E. coli.

Quantitative RNA-seq analysis of the campylobacter jejuni transcriptome

Campylobacter jejuni is the most common bacterial cause of foodborne disease in the developed world. Its general physiology and biochemistry, as well as the mechanisms enabling it to colonize and cause disease in various hosts, are not well understood, and new approaches are required to understand its basic biology. High-throughput sequencing technologies provide unprecedented opportunities for functional genomic research. Recent studies have shown that direct Illumina sequencing of cDNA (RNA-seq) is a useful technique for the quantitative and qualitative examination of transcriptomes. In this study we report RNA-seq analyses of the transcriptomes of C. jejuni (NCTC11168) and its rpoN mutant. This has allowed the identification of hitherto unknown transcriptional units, and further defines the regulon that is dependent on rpoN for expression. The analysis of the NCTC11168 transcriptome was supplemented by additional proteomic analysis using liquid chromatography-MS. The transcriptomic and proteomic datasets represent an important resource for the Campylobacter research community. © 2011 SGM.

RNA processing and Arabidopsis flowering time control.

Plants control their flowering time in order to ensure that they reproduce under favourable conditions. The components involved in this complex process have been identified using a molecular genetic approach in Arabidopsis and classified into genetically separable pathways. The autonomous pathway controls the level of mRNA encoding a floral repressor, FLC, and comprises three RNA-binding proteins, FCA, FPA and FLK. FCA interacts with the 3'-end RNA-processing factor FY to autoregulate its own expression post-transcriptionally and to control FLC. Other components of the autonomous pathway, FVE and FLD, regulate FLC epigenetically. This combination of epigenetic and post-transcriptional control gives precision to the control of FLC expression and flowering time.

Dissecting Arabidopsis thaliana DICER function in small RNA processing, gene silencing and DNA methylation patterning.

Small RNAs have several important biological functions. MicroRNAs (miRNAs) and trans-acting small interfering RNAs (tasiRNAs) regulate mRNA stability and translation, and siRNAs cause post-transcriptional gene silencing of transposons, viruses and transgenes and are important in both the establishment and maintenance of cytosine DNA methylation. Here, we study the role of the four Arabidopsis thaliana DICER-LIKE genes (DCL1-DCL4) in these processes. Sequencing of small RNAs from a dcl2 dcl3 dcl4 triple mutant showed markedly reduced tasiRNA and siRNA production and indicated that DCL1, in addition to its role as the major enzyme for processing miRNAs, has a previously unknown role in the production of small RNAs from endogenous inverted repeats. DCL2, DCL3 and DCL4 showed functional redundancy in siRNA and tasiRNA production and in the establishment and maintenance of DNA methylation. Our studies also suggest that asymmetric DNA methylation can be maintained by pathways that do not require siRNAs.

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.

Role of RNA polymerase IV in plant small RNA metabolism.

In addition to the three RNA polymerases (RNAP I-III) shared by all eukaryotic organisms, plant genomes encode a fourth RNAP (RNAP IV) that appears to be specialized in the production of siRNAs. Available data support a model in which dsRNAs are generated by RNAP IV and RNA-dependent RNAP 2 (RDR2) and processed by DICER (DCL) enzymes into 21- to 24-nt siRNAs, which are associated with different ARGONAUTE (AGO) proteins for transcriptional or posttranscriptional gene silencing. However, it is not yet clear what fraction of genomic siRNA production is RNAP IV-dependent, and to what extent these siRNAs are preferentially processed by certain DCL(s) or associated with specific AGOs for distinct downstream functions. To address these questions on a genome-wide scale, we sequenced approximately 335,000 siRNAs from wild-type and RNAP IV mutant Arabidopsis plants by using 454 technology. The results show that RNAP IV is required for the production of >90% of all siRNAs, which are faithfully produced from a discrete set of genomic loci. Comparisons of these siRNAs with those accumulated in rdr2 and dcl2 dcl3 dcl4 and those associated with AGO1 and AGO4 provide important information regarding the processing, channeling, and functions of plant siRNAs. We also describe a class of RNAP IV-independent siRNAs produced from endogenous single-stranded hairpin RNA precursors.

Epigenetic remodeling of meiotic crossover frequency in Arabidopsis thaliana DNA methyltransferase mutants.

Meiosis is a specialized eukaryotic cell division that generates haploid gametes required for sexual reproduction. During meiosis, homologous chromosomes pair and undergo reciprocal genetic exchange, termed crossover (CO). Meiotic CO frequency varies along the physical length of chromosomes and is determined by hierarchical mechanisms, including epigenetic organization, for example methylation of the DNA and histones. Here we investigate the role of DNA methylation in determining patterns of CO frequency along Arabidopsis thaliana chromosomes. In A. thaliana the pericentromeric regions are repetitive, densely DNA methylated, and suppressed for both RNA polymerase-II transcription and CO frequency. DNA hypomethylated methyltransferase1 (met1) mutants show transcriptional reactivation of repetitive sequences in the pericentromeres, which we demonstrate is coupled to extensive remodeling of CO frequency. We observe elevated centromere-proximal COs in met1, coincident with pericentromeric decreases and distal increases. Importantly, total numbers of CO events are similar between wild type and met1, suggesting a role for interference and homeostasis in CO remodeling. To understand recombination distributions at a finer scale we generated CO frequency maps close to the telomere of chromosome 3 in wild type and demonstrate an elevated recombination topology in met1. Using a pollen-typing strategy we have identified an intergenic nucleosome-free CO hotspot 3a, and we demonstrate that it undergoes increased recombination activity in met1. We hypothesize that modulation of 3a activity is caused by CO remodeling driven by elevated centromeric COs. These data demonstrate how regional epigenetic organization can pattern recombination frequency along eukaryotic chromosomes.