Quantitative whole-genome analysis of DNA-protein interactions by in vivo methylase protection in E. coli

Bacmethy: A novel and convenient tool for investigating bacterial DNA methylation pattern and their transcriptional regulation effects

DNA methylation serves as the primary mode of epigenetic regulation in prokaryotes, particularly through transcriptional regulation. With the rapid implementation of third-generation sequencing technology, we are currently experiencing a golden age of bacterial epigenomics. However, there has been a lack of comprehensive research exploring the versatility and consequential impact of bacterial DNA methylome on cellular and physiological functions. There is a critical need for a user-friendly bioinformatics tool that can effectively characterize DNA methylation modification features and predict the regulation patterns. To address this gap, the current study introduces Bacmethy, an innovative tool that utilizes SMRT-seq data and offers a range of analytical modules. First, the tool classifies methylation sites in the genome, highlighting the distinct regulations present under varying modification fractions and location enrichment. Furthermore, this tool enables us to identify regulatory region methylation and potential cis and trans interactions between methylation sites and regulatory effectors. Using benchmark data sets and our data, we show that our tool facilitates the understanding of the distinctive traits of DNA methylation modifications and predicts transcriptional regulation effects on important physiological and pathological functions. Bacmethy code is freely available, and the Docker image is downloadable. Bacmethy has been made available as a user-friendly web server interface at https://bacmethy.med.sustech.edu.cn.

Differential adenine methylation analysis reveals increased variability in 6mA in the absence of methyl-directed mismatch repair

IMPORTANCE

Methylation greatly influences the bacterial genome by guiding DNA repair and regulating pathogenic and stress-response phenotypes. But, the rate of epigenetic changes and their consequences on molecular phenotypes are underexplored. Through a detailed characterization of genome-wide adenine methylation in a commonly used laboratory strain of Escherichia coli, we reveal that mismatch repair deficient populations experience an increase in epimutations resulting in a genome-wide reduction of 6mA methylation in a manner consistent with genetic drift. Our findings highlight how methylation patterns evolve and the constraints on epigenetic evolution due to post-replicative DNA repair, contributing to a deeper understanding of bacterial genome evolution and how epimutations may introduce semi-permanent variation that can influence adaptation.

Marine DNA methylation patterns are associated with microbial community composition and inform virus-host dynamics

BACKGROUND

DNA methylation in prokaryotes is involved in many different cellular processes including cell cycle regulation and defense against viruses. To date, most prokaryotic methylation systems have been studied in culturable microorganisms, resulting in a limited understanding of DNA methylation from a microbial ecology perspective. Here, we analyze the distribution patterns of several microbial epigenetics marks in the ocean microbiome through genome-centric metagenomics across all domains of life.

RESULTS

We reconstructed 15,056 viral, 252 prokaryotic, 56 giant viral, and 6 eukaryotic metagenome-assembled genomes from northwest Pacific Ocean seawater samples using short- and long-read sequencing approaches. These metagenome-derived genomes mostly represented novel taxa, and recruited a majority of reads. Thanks to single-molecule real-time (SMRT) sequencing technology, base modification could also be detected for these genomes. This showed that DNA methylation can readily be detected across dominant oceanic bacterial, archaeal, and viral populations, and microbial epigenetic changes correlate with population differentiation. Furthermore, our genome-wide epigenetic analysis of Pelagibacter suggests that GANTC, a DNA methyltransferase target motif, is related to the cell cycle and is affected by environmental conditions. Yet, the presence of this motif also partitions the phylogeny of the Pelagibacter phages, possibly hinting at a competitive co-evolutionary history and multiple effects of a single methylation mark.

CONCLUSIONS

Overall, this study elucidates that DNA methylation patterns are associated with ecological changes and virus-host dynamics in the ocean microbiome. Video Abstract.

A GCDGC-specific DNA (cytosine-5) methyltransferase that methylates the GCWGC sequence on both strands and the GCSGC sequence on one strand

5-Methylcytosine is one of the major epigenetic marks of DNA in living organisms. Some bacterial species possess DNA methyltransferases that modify cytosines on both strands to produce fully-methylated sites or on either strand to produce hemi-methylated sites. In this study, we characterized a DNA methyltransferase that produces two sequences with different methylation patterns: one methylated on both strands and another on one strand. M.BatI is the orphan DNA methyltransferase of Bacillus anthracis coded in one of the prophages on the chromosome. Analysis of M.BatI modified DNA by bisulfite sequencing revealed that the enzyme methylates the first cytosine in sequences of 5ʹ-GCAGC-3ʹ, 5ʹ-GCTGC-3ʹ, and 5ʹ-GCGGC-3ʹ, but not of 5ʹ-GCCGC-3ʹ. This resulted in the production of fully-methylated 5ʹ-GCWGC-3ʹ and hemi-methylated 5ʹ-GCSGC-3ʹ. M.BatI also showed toxicity when expressed in E. coli, which was caused by a mechanism other than DNA modification activity. Homologs of M.BatI were found in other Bacillus species on different prophage like regions, suggesting the spread of the gene by several different phages. The discovery of the DNA methyltransferase with unique modification target specificity suggested unrevealed diversity of target sequences of bacterial cytosine DNA methyltransferase.

Roles of adenine methylation and genetic mutations in adaptation to different temperatures in Serratia marcescens

Epigenetic modifications can contribute to adaptation, but the relative contributions of genetic and epigenetic variation are unknown. Previous studies on the role of epigenetic changes in adaptation in eukaryotes have nearly exclusively focused on cytosine methylation (m5C), while prokaryotes exhibit a richer system of methyltransferases targetting adenines (m6A) or cytosines (m4C, m5C). DNA methylation in prokaryotes has many roles, but its potential role in adaptation still needs further investigation. We collected phenotypic, genetic, and epigenetic data using single molecule real-time sequencing of clones of the bacterium Serratia marcescens that had undergone experimental evolution in contrasting temperatures to investigate the relationship between environment and genetic, epigenetic, and phenotypic changes. The genomic distribution of GATC motifs, which were the main target for m6A methylation, and of variable m6A epiloci pointed to a potential link between m6A methylation and regulation of gene expression in S. marcescens. Evolved strains, while genetically homogeneous, exhibited many polymorphic m6A epiloci. There was no strong support for a genetic control of methylation changes in our experiment, and no clear evidence of parallel environmentally induced or environmentally selected methylation changes at specific epiloci was found. Both genetic and epigenetic variants were associated with some phenotypic traits. Overall, our results suggest that both genetic and adenine methylation changes have the potential to contribute to phenotypic adaptation in S. marcescens, but that any environmentally induced epigenetic change occurring in our experiment would probably have been quite labile.

Comparison of long-read sequencing technologies in interrogating bacteria and fly genomes

The newest generation of DNA sequencing technology is highlighted by the ability to generate sequence reads hundreds of kilobases in length. Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) have pioneered competitive long read platforms, with more recent work focused on improving sequencing throughput and per-base accuracy. We used whole-genome sequencing data produced by three PacBio protocols (Sequel II CLR, Sequel II HiFi, RS II) and two ONT protocols (Rapid Sequencing and Ligation Sequencing) to compare assemblies of the bacteria Escherichia coli and the fruit fly Drosophila ananassae. In both organisms tested, Sequel II assemblies had the highest consensus accuracy, even after accounting for differences in sequencing throughput. ONT and PacBio CLR had the longest reads sequenced compared to PacBio RS II and HiFi, and genome contiguity was highest when assembling these datasets. ONT Rapid Sequencing libraries had the fewest chimeric reads in addition to superior quantification of E. coli plasmids versus ligation-based libraries. The quality of assemblies can be enhanced by adopting hybrid approaches using Illumina libraries for bacterial genome assembly or polishing eukaryotic genome assemblies, and an ONT-Illumina hybrid approach would be more cost-effective for many users. Genome-wide DNA methylation could be detected using both technologies, however ONT libraries enabled the identification of a broader range of known E. coli methyltransferase recognition motifs in addition to undocumented D. ananassae motifs. The ideal choice of long read technology may depend on several factors including the question or hypothesis under examination. No single technology outperformed others in all metrics examined.

Dynamic landscape of protein occupancy across the Escherichia coli chromosome

Free-living bacteria adapt to environmental change by reprogramming gene expression through precise interactions of hundreds of DNA-binding proteins. A predictive understanding of bacterial physiology requires us to globally monitor all such protein-DNA interactions across a range of environmental and genetic perturbations. Here, we show that such global observations are possible using an optimized version of in vivo protein occupancy display technology (in vivo protein occupancy display-high resolution, IPOD-HR) and present a pilot application to Escherichia coli. We observe that the E. coli protein-DNA interactome organizes into 2 distinct prototypic features: (1) highly dynamic condition-dependent transcription factor (TF) occupancy; and (2) robust kilobase scale occupancy by nucleoid factors, forming silencing domains analogous to eukaryotic heterochromatin. We show that occupancy dynamics across a range of conditions can rapidly reveal the global transcriptional regulatory organization of a bacterium. Beyond discovery of previously hidden regulatory logic, we show that these observations can be utilized to computationally determine sequence specificity models for the majority of active TFs. Our study demonstrates that global observations of protein occupancy combined with statistical inference can rapidly and systematically reveal the transcriptional regulatory and structural features of a bacterial genome. This capacity is particularly crucial for non-model bacteria that are not amenable to routine genetic manipulation.

A Thermosensitive, Phase-Variable Epigenetic Switch: pap Revisited

It has been more than a decade since the last comprehensive review of the phase-variable uropathogen-associated pyelonephritis-associated pilus (pap) genetic switch. Since then, important data have come to light, including additional factors that regulate pap expression, better characterization of H-NS regulation, the structure of the Lrp octamer in complex with pap regulatory DNA, the temperature-insensitive phenotype of a mutant lacking the acetyltransferase RimJ, evidence that key components of the regulatory machinery are acetylated, and new insights into the role of DNA binding by key regulators in shaping both the physical structure and regulatory state of the papI and papBA promoters. This review revisits pap, integrating these newer observations with older ones to produce a new model for the concerted behavior of this virulence-regulatory region.

DNA Adenine Methyltransferase (Dam) Overexpression Impairs Photorhabdus luminescens Motility and Virulence

Dam, the most described bacterial DNA-methyltransferase, is widespread in gamma-proteobacteria. Dam DNA methylation can play a role in various genes expression and is involved in pathogenicity of several bacterial species. The purpose of this study was to determine the role played by the dam ortholog identified in the entomopathogenic bacterium Photorhabdus luminescens. Complementation assays of an Escherichia coli dam mutant showed the restoration of the DNA methylation state of the parental strain. Overexpression of dam in P. luminescens did not impair growth ability in vitro. In contrast, compared to a control strain harboring an empty plasmid, a significant decrease in motility was observed in the dam-overexpressing strain. A transcriptome analysis revealed the differential expression of 208 genes between the two strains. In particular, the downregulation of flagellar genes was observed in the dam-overexpressing strain. In the closely related bacterium Xenorhabdus nematophila, dam overexpression also impaired motility. In addition, the dam-overexpressing P. luminescens strain showed a delayed virulence compared to that of the control strain after injection in larvae of the lepidopteran Spodoptera littoralis. These results reveal that Dam plays a major role during P. luminescens insect infection.

Independent Evolution of Winner Traits without Whole Genome Duplication in Dekkera Yeasts

Dekkera yeasts have often been considered as alternative sources of ethanol production that could compete with S. cerevisiae. The two lineages of yeasts independently evolved traits that include high glucose and ethanol tolerance, aerobic fermentation, and a rapid ethanol fermentation rate. The Saccharomyces yeasts attained these traits mainly through whole genome duplication approximately 100 million years ago (Mya). However, the Dekkera yeasts, which were separated from S. cerevisiae approximately 200 Mya, did not undergo whole genome duplication (WGD) but still occupy a niche similar to S. cerevisiae. Upon analysis of two Dekkera yeasts and five closely related non-WGD yeasts, we found that a massive loss of cis-regulatory elements occurred in an ancestor of the Dekkera yeasts, which led to improved mitochondrial functions similar to the S. cerevisiae yeasts. The evolutionary analysis indicated that genes involved in the transcription and translation process exhibited faster evolution in the Dekkera yeasts. We detected 90 positively selected genes, suggesting that the Dekkera yeasts evolved an efficient translation system to facilitate adaptive evolution. Moreover, we identified that 12 vacuolar H+-ATPase (V-ATPase) function genes that were under positive selection, which assists in developing tolerance to high alcohol and high sugar stress. We also revealed that the enzyme PGK1 is responsible for the increased rate of glycolysis in the Dekkera yeasts. These results provide important insights to understand the independent adaptive evolution of the Dekkera yeasts and provide tools for genetic modification promoting industrial usage.

A role for the bacterial GATC methylome in antibiotic stress survival

Antibiotic resistance is an increasingly serious public health threat¹. Understanding pathways allowing bacteria to survive antibiotic stress may unveil new therapeutic targets^2–8. We explore the role of the bacterial epigenome in antibiotic stress survival using classical genetic tools and single-molecule real-time sequencing to characterize genomic methylation kinetics. We find that Escherichia coli survival under antibiotic pressure is severely compromised without adenine methylation at GATC sites. While the adenine methylome remains stable during drug stress, without GATC methylation, methyl-dependent mismatch repair (MMR) is deleterious, and fueled by the drug-induced error-prone polymerase PolIV, overwhelms cells with toxic DNA breaks. In multiple E. coli strains, including pathogenic and drug-resistant clinical isolates, DNA adenine methyltransferase deficiency potentiates antibiotics from the β-lactam and quinolone classes. This work indicates that the GATC methylome provides structural support for bacterial survival during antibiotics stress and suggests targeting bacterial DNA methylation as a viable approach to enhancing antibiotic activity.

The Epigenomic Landscape of Prokaryotes

DNA methylation acts in concert with restriction enzymes to protect the integrity of prokaryotic genomes. Studies in a limited number of organisms suggest that methylation also contributes to prokaryotic genome regulation, but the prevalence and properties of such non-restriction-associated methylation systems remain poorly understood. Here, we used single molecule, real-time sequencing to map DNA modifications including m6A, m4C, and m5C across the genomes of 230 diverse bacterial and archaeal species. We observed DNA methylation in nearly all (93%) organisms examined, and identified a total of 834 distinct reproducibly methylated motifs. This data enabled annotation of the DNA binding specificities of 620 DNA Methyltransferases (MTases), doubling known specificities for previously hard to study Type I, IIG and III MTases, and revealing their extraordinary diversity. Strikingly, 48% of organisms harbor active Type II MTases with no apparent cognate restriction enzyme. These active ‘orphan’ MTases are present in diverse bacterial and archaeal phyla and show motif specificities and methylation patterns consistent with functions in gene regulation and DNA replication. Our results reveal the pervasive presence of DNA methylation throughout the prokaryotic kingdoms, as well as the diversity of sequence specificities and potential functions of DNA methylation systems.

DNA methylation is a chemical modification of DNA present in many prokaryotic genomes. The best-known role of DNA methylation is as a component of restriction-modification systems. In these systems, restriction enzymes target foreign DNA for cleavage, while DNA methylation protects the host genome from destruction. Studies in a handful of organisms show that DNA methylation may also act independently of restriction systems and function in genome regulation. However, a lack of technologies has limited the study of DNA methylation to a small number of organisms, and the broader patterns and functions of DNA methylation remain unknown. Here we use SMRT-sequencing to determine the genome wide DNA methylation patterns of more than 200 diverse bacteria and archaea. We show that DNA methylation is pervasive and present in more than 90% of studied organisms. Analysis of this data enabled annotation of the specific DNA binding sites of more than 600 restriction systems, revealing their extraordinary diversity. Strikingly, we observed widespread DNA methylation in the absence of restriction systems. Analyses of these patterns reveal that they are conserved through evolution, and likely function in genome regulation. Thus DNA methylation may play a far wider function in prokaryotic genome biology than was previously supposed.

Bacterial DNA Methylation and Methylomes

Formation of C5-methylcytosine, N4-methylcytosine, and N6-methyladenine in bacterial genomes is postreplicative and involves transfer of a methyl group from S-adenosyl-methionine to a base embedded in a specific DNA sequence context. Most bacterial DNA methyltransferases belong to restriction-modification systems; in addition, "solitary" or "orphan" DNA methyltransferases are frequently found in the genomes of bacteria and phage. Base methylation can affect the interaction of DNA-binding proteins with their cognate sites, either by a direct effect (e.g., steric hindrance) or by changes in DNA topology. In both Alphaproteobacteria and Gammaproteobacteria, the roles of DNA base methylation are especially well known for N6-methyladenine, including control of chromosome replication, nucleoid segregation, postreplicative correction of DNA mismatches, cell cycle-coupled transcription, formation of bacterial cell lineages, and regulation of bacterial virulence. Technical procedures that permit genome-wide analysis of DNA methylation are nowadays expanding our knowledge of the extent, evolution, and physiological significance of bacterial DNA methylation.

DNA Methylation

The DNA of Escherichia coli contains 19,120 6-methyladenines and 12,045 5-methylcytosines in addition to the four regular bases, and these are formed by the postreplicative action of three DNA methyltransferases. The majority of the methylated bases are formed by the Dam and Dcm methyltransferases encoded by the dam (DNA adenine methyltransferase) and dcm (DNA cytosine methyltransferase) genes. Although not essential, Dam methylation is important for strand discrimination during the repair of replication errors, controlling the frequency of initiation of chromosome replication at oriC, and the regulation of transcription initiation at promoters containing GATC sequences. In contrast, there is no known function for Dcm methylation, although Dcm recognition sites constitute sequence motifs for Very Short Patch repair of T/G base mismatches. In certain bacteria (e.g., Vibrio cholerae, Caulobacter crescentus) adenine methylation is essential, and, in C. crescentus, it is important for temporal gene expression, which, in turn, is required for coordinating chromosome initiation, replication, and division. In practical terms, Dam and Dcm methylation can inhibit restriction enzyme cleavage, decrease transformation frequency in certain bacteria, and decrease the stability of short direct repeats and are necessary for site-directed mutagenesis and to probe eukaryotic structure and function.

Single molecule-level detection and long read-based phasing of epigenetic variations in bacterial methylomes

Beyond its role in host defense, bacterial DNA methylation also plays important roles in the regulation of gene expression, virulence and antibiotic resistance. Bacterial cells in a clonal population can generate epigenetic heterogeneity to increase population-level phenotypic plasticity. Single molecule, real-time (SMRT) sequencing enables the detection of N6-methyladenine and N4-methylcytosine, two major types of DNA modifications comprising the bacterial methylome. However, existing SMRT sequencing-based methods for studying bacterial methylomes rely on a population-level consensus that lacks the single-cell resolution required to observe epigenetic heterogeneity. Here, we present SMALR (single-molecule modification analysis of long reads), a novel framework for single molecule-level detection and phasing of DNA methylation. Using seven bacterial strains, we show that SMALR yields significantly improved resolution and reveals distinct types of epigenetic heterogeneity. SMALR is a powerful new tool that enables de novo detection of epigenetic heterogeneity and empowers investigation of its functions in bacterial populations.

Bacterial DNA methylation is involved in many processes, from host defense to antibiotic resistance, however current methods for examining methylated genomes lack single-cell resolution. Here Beaulaurier et al. present Single Molecule Modification Analysis of Long Reads, a new tool for de novo detection of epigenetic heterogeneity.

Genome-wide methylation patterns in Salmonella enterica Subsp. enterica Serovars

The methylation of DNA bases plays an important role in numerous biological processes including development, gene expression, and DNA replication. Salmonella is an important foodborne pathogen, and methylation in Salmonella is implicated in virulence. Using single molecule real-time (SMRT) DNA-sequencing, we sequenced and assembled the complete genomes of eleven Salmonella enterica isolates from nine different serovars, and analysed the whole-genome methylation patterns of each genome. We describe 16 distinct N⁶-methyladenine (m6A) methylated motifs, one N⁴-methylcytosine (m4C) motif, and one combined m6A-m4C motif. Eight of these motifs are novel, i.e., they have not been previously described. We also identified the methyltransferases (MTases) associated with 13 of the motifs. Some motifs are conserved across all Salmonella serovars tested, while others were found only in a subset of serovars. Eight of the nine serovars contained a unique methylated motif that was not found in any other serovar (most of these motifs were part of Type I restriction modification systems), indicating the high diversity of methylation patterns present in Salmonella.

Microarray-based detection of Salmonella enterica serovar Enteritidis genes involved in chicken reproductive tract colonization

Salmonella enterica serovar Enteritidis has developed the potential to contaminate table eggs internally, by colonization of the chicken reproductive tract and internalization in the forming egg. The serotype Enteritidis has developed mechanisms to colonize the chicken oviduct more successfully than other serotypes. Until now, the strategies exploited by Salmonella Enteritidis to do so have remained largely unknown. For that reason, a microarray-based transposon library screen was used to identify genes that are essential for the persistence of Salmonella Enteritidis inside primary chicken oviduct gland cells in vitro and inside the reproductive tract in vivo. A total of 81 genes with a potential role in persistence in both the oviduct cells and the oviduct tissue were identified. Major groups of importance include the Salmonella pathogenicity islands 1 and 2, genes involved in stress responses, cell wall, and lipopolysaccharide structure, and the region-of-difference genomic islands 9, 21, and 40.

A genome-wide screen identifies Salmonella Enteritidis lipopolysaccharide biosynthesis and the HtrA heat shock protein as crucial factors involved in egg white persistence at chicken body temperature

Eggs contaminated with Salmonella Enteritidis are an important source of human foodborne Salmonella infections. Salmonella Enteritidis is able to contaminate egg white during formation of the egg within the chicken oviduct, and it has developed strategies to withstand the antimicrobial properties of egg white to survive in this hostile environment. The mechanisms involved in the persistence of Salmonella Enteritidis in egg white are likely to be complex. To address this issue, a microarray-based transposon library screen was performed to identify genes necessary for survival of Salmonella Enteritidis in egg white at chicken body temperature. The majority of identified genes belonged to the lipopolysaccharide biosynthesis pathway. Additionally, we provide evidence that the serine protease/heat shock protein (HtrA) appears essential for the survival of Salmonella Enteritidis in egg white at chicken body temperature.

Genetic evidence for the involvement of the S-layer protein gene sap and the sporulation genes spo0A, spo0B, and spo0F in Phage AP50c infection of Bacillus anthracis

In order to better characterize the Bacillus anthracis typing phage AP50c, we designed a genetic screen to identify its bacterial receptor. Insertions of the transposon mariner or targeted deletions of the structural gene for the S-layer protein Sap and the sporulation genes spo0A, spo0B, and spo0F in B. anthracis Sterne resulted in phage resistance with concomitant defects in phage adsorption and infectivity. Electron microscopy of bacteria incubated with AP50c revealed phage particles associated with the surface of bacilli of the Sterne strain but not with the surfaces of Δsap, Δspo0A, Δspo0B, or Δspo0F mutants. The amount of Sap in the S layer of each of the spo0 mutant strains was substantially reduced compared to that of the parent strain, and incubation of AP50c with purified recombinant Sap led to a substantial reduction in phage activity. Phylogenetic analysis based on whole-genome sequences of B. cereus sensu lato strains revealed several closely related B. cereus and B. thuringiensis strains that carry sap genes with very high similarities to the sap gene of B. anthracis. Complementation of the Δsap mutant in trans with the wild-type B. anthracis sap or the sap gene from either of two different B. cereus strains that are sensitive to AP50c infection restored phage sensitivity, and electron microscopy confirmed attachment of phage particles to the surface of each of the complemented strains. Based on these data, we postulate that Sap is involved in AP50c infectivity, most likely acting as the phage receptor, and that the spo0 genes may regulate synthesis of Sap and/or formation of the S layer.

Exploring the roles of DNA methylation in the metal-reducing bacterium Shewanella oneidensis MR-1

We performed whole-genome analyses of DNA methylation in Shewanella oneidensis MR-1 to examine its possible role in regulating gene expression and other cellular processes. Single-molecule real-time (SMRT) sequencing revealed extensive methylation of adenine (N6mA) throughout the genome. These methylated bases were located in five sequence motifs, including three novel targets for type I restriction/modification enzymes. The sequence motifs targeted by putative methyltranferases were determined via SMRT sequencing of gene knockout mutants. In addition, we found that S. oneidensis MR-1 cultures grown under various culture conditions displayed different DNA methylation patterns. However, the small number of differentially methylated sites could not be directly linked to the much larger number of differentially expressed genes under these conditions, suggesting that DNA methylation is not a major regulator of gene expression in S. oneidensis MR-1. The enrichment of methylated GATC motifs in the origin of replication indicates that DNA methylation may regulate genome replication in a manner similar to that seen in Escherichia coli. Furthermore, comparative analyses suggest that many Gammaproteobacteria, including all members of the Shewanellaceae family, may also utilize DNA methylation to regulate genome replication.

Impaired cardiovascular function caused by different stressors elicits a common pathological and transcriptional response in zebrafish embryos

Zebrafish embryos have been widely used to study the genes and processes needed for normal vertebrate heart development. We recently observed that exposure to 2,3,7,8-tetra-chlorodibenzo-p-dioxin (TCDD) or retinoic acid (RA) produces very similar signs of heart failure in developing zebrafish via divergent molecular pathways. The fact that diverse stressors and mutations cause severe pericardial edema and circulatory collapse in developing zebrafish has been largely unexplored. We hypothesized that unrelated chemicals can trigger a common pathological response leading to the same end-stage heart failure. To test this hypothesis, we compared the effects of TCDD, RA, carbaryl, valproic acid, and morpholino oligonucleotide (MO) knockdown of TBX5 on the developing heart in zebrafish embryos. These model stressors have all been previously reported to affect zebrafish heart development, and elicited very similar signs of embryonic heart failure. Microarray analysis showed that one cluster of 92 transcripts affected by these different treatments was significantly downregulated by all treatments. This gene cluster is composed of transcripts required for chromosome assembly, DNA replication, and cell cycle progression. We refer to this cluster as the cell cycle gene cluster (CCGC). Immunohistochemistry revealed that downregulation of the CCGC precedes a halt in cardiomyocyte proliferation in the hearts of zebrafish exposed to any of the treatments. Previous work has shown that the initial response to TCDD is a decrease in cardiac output. Since this precedes the signs of edema, heart failure, and fall in CCGC expression, we postulated that any factor that decreases cardiac output will produce the same syndrome of heart failure responses. To test this, we used MO knockdown of cardiac troponin T2 (TNNT2) to specifically block contractility. The TNNT2-MO produced exactly the same signs of cardiotoxicity as the other treatments, including downregulation of the signature CCGC. Our results indicate that agents altering cardiac output can have amplified consequences during specific periods in development.

Characterization of undermethylated sites in Vibrio cholerae

The activities of DNA methyltransferases are important for a variety of cellular functions in bacteria. In this study, we developed a modified high-throughput technique called methyl homopolymer tail mediated sequencing (methyl HTM-seq) to identify the undermethylated sites in the Vibrio cholerae genome for the two DNA methyltransferases, Dam, an adenine methyltransferase, and VchM, a cytosine methyltransferase, during growth in rich medium in vitro. Many of the undermethylated sites occurred in intergenic regions, and for most of these sites, we identified the transcription factors responsible for undermethylation. This confirmed the presence of previously hypothesized DNA-protein interactions for these transcription factors and provided insight into the biological state of these cells during growth in vitro. DNA adenine methylation has previously been shown to mediate heritable epigenetic switches in gene regulation. However, none of the undermethylated Dam sites tested showed evidence of regulation by this mechanism. This study is the first to identify undermethylated adenines and cytosines genomewide in a bacterium using second-generation sequencing technology.

Proximal recognition sites facilitate intrasite hopping by DNA adenine methyltransferase: mechanistic exploration of epigenetic gene regulation

The methylation of adenine in palindromic 5'-GATC-3' sites by Escherichia coli Dam supports diverse roles, including the essential regulation of virulence genes in several human pathogens. As a result of a unique hopping mechanism, Dam methylates both strands of the same site prior to fully dissociating from the DNA, a process referred to as intrasite processivity. The application of a DpnI restriction endonuclease-based assay allowed the direct interrogation of this mechanism with a variety of DNA substrates. Intrasite processivity is disrupted when the DNA flanking a single GATC site is longer than 400 bp on either side. Interestingly, the introduction of a second GATC site within this flanking DNA reinstates intrasite methylation of both sites. Our results show that intrasite methylation occurs only when GATC sites are clustered, as is found in gene segments both known and postulated to undergo in vivo epigenetic regulation by Dam methylation. We propose a model for intrasite methylation in which Dam bound to flanking DNA is an obligate intermediate. Our results provide insights into how intrasite processivity, which appears to be context-dependent, may contribute to the diverse biological roles that are carried out by Dam.

Accurate prediction of gene expression by integration of DNA sequence statistics with detailed modeling of transcription regulation

Gene regulation involves a hierarchy of events that extend from specific protein-DNA interactions to the combinatorial assembly of nucleoprotein complexes. The effects of DNA sequence on these processes have typically been studied based either on its quantitative connection with single-domain binding free energies or on empirical rules that combine different DNA motifs to predict gene expression trends on a genomic scale. The middle-point approach that quantitatively bridges these two extremes, however, remains largely unexplored. Here, we provide an integrated approach to accurately predict gene expression from statistical sequence information in combination with detailed biophysical modeling of transcription regulation by multidomain binding on multiple DNA sites. For the regulation of the prototypical lac operon, this approach predicts within 0.3-fold accuracy transcriptional activity over a 10,000-fold range from DNA sequence statistics for different intracellular conditions.

Epigenetic regulation of the bacterial cell cycle

N(6)-methyl-adenines can serve as epigenetic signals for interactions between regulatory DNA sequences and regulatory proteins that control cellular functions, such as the initiation of chromosome replication or the expression of specific genes. Several of these genes encode master regulators of the bacterial cell cycle. DNA adenine methylation is mediated by Dam in gamma-proteobacteria and by CcrM in alpha-proteobacteria. A major difference between them is that CcrM is cell cycle regulated, while Dam is active throughout the cell cycle. In alpha-proteobacteria, GANTC sites can remain hemi-methylated for a significant period of the cell cycle, depending on their location on the chromosome. In gamma-proteobacteria, most GATC sites are only transiently hemi-methylated, except regulatory GATC sites that are protected from Dam methylation by specific DNA-binding proteins.

Escherichia coli DNA adenine methyltransferase: intrasite processivity and substrate-induced dimerization and activation

Methylation of GATC sites in Escherichia coli by DNA adenine methyltransferase (EcoDam) is essential for proper DNA replication timing, gene regulation, and mismatch repair. The low cellular concentration of EcoDam and the high number of GATC sites in the genome (approximately 20000) support the reliance on methylation efficiency-enhancing strategies such as extensive intersite processivity. Here, we present evidence that EcoDam has evolved other unique mechanisms of activation not commonly observed with restriction-modification methyltransferases. EcoDam dimerizes on short, synthetic DNA, resulting in enhanced catalysis; however, dimerization is not observed on large genomic DNA where the potential for intersite processive methylation precludes any dimerization-dependent activation. An activated form of the enzyme is apparent on large genomic DNA and can also be achieved with high concentrations of short, synthetic substrates. We suggest that this activation is inherent on polymeric DNA where either multiple GATC sites are available for methylation or the partitioning of the enzyme onto nonspecific DNA is favored. Unlike other restriction-modification methyltransferases, EcoDam carries out intrasite processive catalysis whereby the enzyme-DNA complex methylates both strands of an unmethylated GATC site prior to dissociation from the DNA. This occurs with short 21 bp oligonucleotides and is highly dependent upon salt concentrations. Kinetic modeling which invokes enzyme activation by both dimerization and excess substrate provides mechanistic insights into key regulatory checkpoints for an enzyme involved in multiple, diverse biological pathways.

Protein occupancy landscape of a bacterial genome

Protein-DNA interactions are fundamental to core biological processes, including transcription, DNA replication, and chromosomal organization. We have developed in vivo protein occupancy display (IPOD), a technology that reveals protein occupancy across an entire bacterial chromosome at the resolution of individual binding sites. Application to Escherichia coli reveals thousands of protein occupancy peaks, highly enriched within and in close proximity to noncoding regulatory regions. In addition, we discovered extensive (>1 kilobase) protein occupancy domains (EPODs), some of which are localized to highly expressed genes, enriched in RNA-polymerase occupancy. However, the majority are localized to transcriptionally silent loci dominated by conserved hypothetical ORFs. These regions are highly enriched in both predicted and experimentally determined binding sites of nucleoid proteins and exhibit extreme biophysical characteristics such as high intrinsic curvature. Our observations implicate these transcriptionally silent EPODs as the elusive organizing centers, long proposed to topologically isolate chromosomal domains.

DNA Methylation

The DNA of Escherichia coli contains 19,120 6-methyladenines and 12,045 5-methylcytosines in addition to the four regular bases, and these are formed by the postreplicative action of three DNA methyltransferases. The majority of the methylated bases are formed by the Dam and Dcmmethyltransferases encoded by the dam (DNA adenine methyltransferase) and dcm (DNA cytosine methyltransferase) genes. Although not essential, Dam methylation is important for strand discrimination during repair of replication errors, controlling the frequency of initiation of chromosome replication at oriC, and regulation of transcription initiation at promoters containing GATC sequences. In contrast, there is no known function for Dcm methylation, although Dcm recognition sites constitute sequence motifs for Very Short Patch repair of T/G base mismatches. In certain bacteria (e.g., Vibrio cholera and Caulobactercrescentus) adenine methylation is essential, and in C.crescentus it is important for temporal gene expression which, in turn, is required for coordination of chromosome initiation, replication, and division. In practical terms, Dam and Dcm methylation can inhibit restriction enzyme cleavage,decrease transformation frequency in certain bacteria,and decrease the stability of short direct repeats andare necessary for site-directed mutagenesis and to probe eukaryotic structure and function.

Escherichia coli DNA adenine methyltransferase: the structural basis of processive catalysis and indirect read-out

We have investigated the structural basis of processive GATC methylation by the Escherichia coli DNA adenine methyltransferase, which is critical in chromosome replication and mismatch repair. We determined the contribution of the orthologically conserved phosphate interactions involving residues Arg(95), Asn(126), Asn(132), Arg(116), and Lys(139), which directly contact the DNA outside the cognate recognition site (GATC) to processive catalysis, and that of residue Arg(137), which is not conserved and contacts the DNA backbone within the GATC sequence. Alanine substitutions at the conserved positions have large impacts on processivity yet do not impact k(cat)/K(m)(DNA) or DNA affinity (K(D)(DNA)). However, these mutants cause large preferences for GATC sites varying in flanking sequences when considering the pre-steady state efficiency constant k(chem)/K(D)(DNA). These changes occur mainly at the level of the methylation rate constant, which results in the observed decreases in processive catalysis. Thus, processivity and catalytic efficiency (k(cat)/K(m)(DNA)) are uncoupled in these mutants. These results reveal that the binding energy involved in DNA recognition contributes to the assembly of the active site rather than tight binding. Furthermore, the conserved residues (Arg(95), Asn(126), Asn(132), and Arg(116)) repress the modulation of the response of the enzyme to flanking sequence effects. Processivity impacted mutants do not show substrate-induced dimerization as is observed for the wild type enzyme. This study describes the structural means by which an enzyme that does not completely enclose its substrate has evolved to achieve processive catalysis, and how interactions with DNA flanking the recognition site alter this processivity.

Development of a mariner-based transposon and identification of Listeria monocytogenes determinants, including the peptidyl-prolyl isomerase PrsA2, that contribute to its hemolytic phenotype

Listeriolysin O (LLO) is a pore-forming toxin that mediates phagosomal escape and cell-to-cell spread of the intracellular pathogen Listeria monocytogenes. In order to identify factors that control the production, activity, or secretion of this essential virulence factor, we constructed a Himar1 mariner transposon delivery system and screened 50,000 mutants for a hypohemolytic phenotype on blood agar plates. Approximately 200 hypohemolytic mutants were identified, and the 51 most prominent mutants were screened ex vivo for intracellular growth defects. Eight mutants with a phenotype were identified, and they contained insertions in the following genes: lmo0964 (similar to yjbH), lmo1268 (clpX), lmo1401 (similar to ymdB), lmo1575 (similar to ytqI), lmo1695 (mprF), lmo1821 (similar to prpC), lmo2219 (prsA2), and lmo2460 (similar to cggR). Some of these genes are involved in previously unexplored areas of research with L. monocytogenes: the genes yjbH and clpX regulate the disulfide stress response in Bacillus subtilis, and the prpC phosphatase has been implicated in virulence in other gram-positive pathogens. Here we demonstrate that prsA2, an extracytoplasmic peptidyl-prolyl cis/trans isomerase, is critical for virulence and contributes to the folding of LLO and to the activity of another virulence factor, the broad-range phospholipase C (PC-PLC). Furthermore, although it has been shown that prsA2 expression is linked to PrfA, the master virulence transcription factor in L. monocytogenes pathogenesis, we demonstrate that prsA2 is not directly controlled by PrfA. Finally, we show that PrsA2 is involved in flagellum-based motility, indicating that this factor likely serves a broad physiological role.

Roles of DNA adenine methylation in host-pathogen interactions: mismatch repair, transcriptional regulation, and more

The DNA adenine methyltransferase (Dam methylase) of Gammaproteobacteria and the cell cycle-regulated methyltransferase (CcrM) methylase of Alphaproteobacteria catalyze an identical reaction (methylation of adenosine moieties using S-adenosyl-methionine as a methyl donor) at similar DNA targets (GATC and GANTC, respectively). Dam and CcrM are of independent evolutionary origin. Each may have evolved from an ancestral restriction-modification system that lost its restriction component, leaving an 'orphan' methylase devoted solely to epigenetic genome modification. The formation of 6-methyladenine reduces the thermodynamic stability of DNA and changes DNA curvature. As a consequence, the methylation state of specific adenosine moieties can affect DNA-protein interactions. Well-known examples include binding of the replication initiation complex to the methylated oriC, recognition of hemimethylated GATCs in newly replicated DNA by the MutHLS mismatch repair complex, and discrimination of methylation states in promoters and regulatory DNA motifs by RNA polymerase and transcription factors. In recent years, Dam and CcrM have been shown to play roles in host-pathogen interactions. These roles are diverse and have only partially been understood. Especially intriguing is the evidence that Dam methylation regulates virulence genes in Escherichia coli, Salmonella, and Yersinia at the posttranscriptional level.

Competitive Lrp and Dam assembly at the pap regulatory region: implications for mechanisms of epigenetic regulation

Escherichia coli DNA adenine methyltransferase (Dam) and Leucine-responsive regulatory protein (Lrp) are key regulators of the pap operon, which codes for the pilus proteins necessary for uropathogenic E. coli cellular adhesion. The pap operon is regulated by a phase variation mechanism in which the methylation states of two GATC sites in the pap regulatory region and the binding position of Lrp determine whether the pilus genes are expressed. The post-replicative reassembly of Dam, Lrp, and the local regulator PapI onto a hemimethylated pap intermediate is a critical step of the phase variation switching mechanism and is not well understood. We show that Lrp, in the presence and in the absence of PapI and nonspecific DNA, specifically protects pap regulatory GATC sites from Dam methylation when allowed to compete with Dam for assembly on unmethylated and hemimethylated pap DNA. The methylation protection is dependent upon the concentration of Lrp and does not occur with non-regulatory GATC sites. Our data suggest that only at low Lrp concentrations will Dam compete effectively for binding and methylation of the proximal GATC site, leading to a phase switch resulting in the expression of pili.

Modulation of Escherichia coli DNA methyltransferase activity by biologically derived GATC-flanking sequences

Escherichia coli DNA adenine methyltransferase (EcoDam) methylates the N-6 position of the adenine in the sequence 5'-GATC-3' and plays vital roles in gene regulation, mismatch repair, and DNA replication. It remains unclear how the small number of critical GATC sites involved in the regulation of replication and gene expression are differentially methylated, whereas the approximately 20,000 GATCs important for mismatch repair and dispersed throughout the genome are extensively methylated. Our prior work, limited to the pap regulon, showed that methylation efficiency is controlled by sequences immediately flanking the GATC sites. We extend these studies to include GATC sites involved in diverse gene regulatory and DNA replication pathways as well as sites previously shown to undergo differential in vivo methylation but whose function remains to be assigned. EcoDam shows no change in affinity with variations in flanking sequences derived from these sources, but methylation kinetics varied 12-fold. A-tracts immediately adjacent to the GATC site contribute significantly to these differences in methylation kinetics. Interestingly, only when the poly(A) is located 5' of the GATC are the changes in methylation kinetics revealed. Preferential methylation is obscured when two GATC sites are positioned on the same DNA molecule, unless both sites are surrounded by large amounts of nonspecific DNA. Thus, facilitated diffusion and sequences immediately flanking target sites contribute to higher order specificity for EcoDam; we suggest that the diverse biological roles of the enzyme are in part regulated by these two factors, which may be important for other enzymes that sequence-specifically modify DNA.

A graph-based approach to systematically reconstruct human transcriptional regulatory modules

MOTIVATION

A major challenge in studying gene regulation is to systematically reconstruct transcription regulatory modules, which are defined as sets of genes that are regulated by a common set of transcription factors. A commonly used approach for transcription module reconstruction is to derive coexpression clusters from a microarray dataset. However, such results often contain false positives because genes from many transcription modules may be simultaneously perturbed upon a given type of conditions. In this study, we propose and validate that genes, which form a coexpression cluster in multiple microarray datasets across diverse conditions, are more likely to form a transcription module. However, identifying genes coexpressed in a subset of many microarray datasets is not a trivial computational problem.

RESULTS

We propose a graph-based data-mining approach to efficiently and systematically identify frequent coexpression clusters. Given m microarray datasets, we model each microarray dataset as a coexpression graph, and search for vertex sets which are frequently densely connected across [theta m] datasets (0 < or = theta < or = 1). For this novel graph-mining problem, we designed two techniques to narrow down the search space: (1) partition the input graphs into (overlapping) groups sharing common properties; (2) summarize the vertex neighbor information from the partitioned datasets onto the 'Neighbor Association Summary Graph's for effective mining. We applied our method to 105 human microarray datasets, and identified a large number of potential transcription modules, activated under different subsets of conditions. Validation by ChIP-chip data demonstrated that the likelihood of a coexpression cluster being a transcription module increases significantly with its recurrence. Our method opens a new way to exploit the vast amount of existing microarray data accumulation for gene regulation study. Furthermore, the algorithm is applicable to other biological networks for approximate network module mining.

AVAILABILITY

http://zhoulab.usc.edu/NeMo/.

A comprehensive genetic characterization of bacterial motility

We have developed a powerful experimental framework that combines competitive selection and microarray-based genetic footprinting to comprehensively reveal the genetic basis of bacterial behaviors. Application of this method to Escherichia coli motility identifies 95% of the known flagellar and chemotaxis genes, and reveals three dozen novel loci that, to varying degrees and through diverse mechanisms, affect motility. To probe the network context in which these genes function, we developed a method that uncovers genome-wide epistatic interactions through comprehensive analyses of double-mutant phenotypes. This allows us to place the novel genes within the context of signaling and regulatory networks, including the Rcs phosphorelay pathway and the cyclic di-GMP second-messenger system. This unifying framework enables sensitive and comprehensive genetic characterization of complex behaviors across the microbial biosphere.

Bacteria thrive in a limitless range of extreme environments, accompanied by exotic metabolisms and sophisticated behaviors. However, our modern molecular understanding of bacteria comes from studies of a limited range of phenotypes in a handful of model organisms such as E. coli and Bacillus subtilis. With the availability of thousands of sequenced bacterial genomes, there is now an urgent need for methods that rapidly and comprehensively reveal the genetic basis of phenotypes across the microbial biosphere. To this end, we have developed a genome-wide experimental framework that quantifies the degree to which every gene in the genome contributes to a phenotype of interest, and reveals the organization of genes within regulatory networks and signaling pathways. We show here that the application of this methodology to E. coli swimming and surface motility reveals essentially all the previously known components of flagellar-mediated chemotaxis on the time scale of weeks. Remarkably, we also identify three dozen additional novel loci that operate through diverse mechanisms to affect a behavior that was assumed to be completely characterized. The speed, ease, and broad applicability of this framework should greatly accelerate the global analysis of a wide range of uncharacterized bacterial behaviors.

DNA adenine methylation changes dramatically during establishment of symbiosis

The DNA adenine methylation status on specific 5'-GANTC-3' sites and its change during the establishment of plant-microbe interactions was demonstrated in several species of alpha-proteobacteria. Restriction landmark genome scanning (RLGS), which is a high-resolution two dimensional DNA electrophoresis method, was used to monitor the genomewide change in methylation. In the case of Mesorhizobium loti MAFF303099, real RLGS images obtained with the restriction enzyme MboI, which digests at GATC sites, almost perfectly matched the virtual RLGS images generated based on genome sequences. However, only a few spots were observed when the restriction enzyme HinfI was used, suggesting that most GANTC (HinfI) sites were tightly methylated and specific sites were unmethylated. DNA gel blot analysis with the cloned specifically unmethylated regions (SUMs) showed that some SUMs were methylated differentially in bacteroids compared to free-living bacteria. SUMs have also been identified in other symbiotic and parasitic bacteria. These results suggest that DNA adenine methylation may contribute to the establishment and/or maintenance of symbiotic and parasitic relationships.

Altered levels of Salmonella DNA adenine methylase are associated with defects in gene expression, motility, flagellar synthesis, and bile resistance in the pathogenic strain 14028 but not in the laboratory strain LT2

Comparative genomic analysis has revealed limited strain diversity between Salmonella pathogenic and nonpathogenic isolates. Thus, some of the relative virulence and host-immune response disparities may be credited to differential gene regulation rather than gross differences in genomic content. Here we show that altered levels of Salmonella DNA adenine methylase (Dam) resulted in acute defects in virulence-associated gene expression, motility, flagellin synthesis, and bile resistance in the Salmonella pathogenic strain 14028 but not in avirulent laboratory strain LT2. The defects in motility exhibited by 14028 in response to altered Dam levels was not dependent on the presence of the regulatory protein, RpoS. The transitioning between flagellar types (phase variation) was also differentially regulated in 14028 versus LT2 in response to dam levels, resulting in distinct differences in flagellin expression states. These data suggest that differential gene regulation may contribute to the relative virulence disparities observed between Salmonella serovars that are closely related at the DNA level.

Epigenetic gene regulation in the bacterial world

Like many eukaryotes, bacteria make widespread use of postreplicative DNA methylation for the epigenetic control of DNA-protein interactions. Unlike eukaryotes, however, bacteria use DNA adenine methylation (rather than DNA cytosine methylation) as an epigenetic signal. DNA adenine methylation plays roles in the virulence of diverse pathogens of humans and livestock animals, including pathogenic Escherichia coli, Salmonella, Vibrio, Yersinia, Haemophilus, and Brucella. In Alphaproteobacteria, methylation of adenine at GANTC sites by the CcrM methylase regulates the cell cycle and couples gene transcription to DNA replication. In Gammaproteobacteria, adenine methylation at GATC sites by the Dam methylase provides signals for DNA replication, chromosome segregation, mismatch repair, packaging of bacteriophage genomes, transposase activity, and regulation of gene expression. Transcriptional repression by Dam methylation appears to be more common than transcriptional activation. Certain promoters are active only during the hemimethylation interval that follows DNA replication; repression is restored when the newly synthesized DNA strand is methylated. In the E. coli genome, however, methylation of specific GATC sites can be blocked by cognate DNA binding proteins. Blockage of GATC methylation beyond cell division permits transmission of DNA methylation patterns to daughter cells and can give rise to distinct epigenetic states, each propagated by a positive feedback loop. Switching between alternative DNA methylation patterns can split clonal bacterial populations into epigenetic lineages in a manner reminiscent of eukaryotic cell differentiation. Inheritance of self-propagating DNA methylation patterns governs phase variation in the E. coli pap operon, the agn43 gene, and other loci encoding virulence-related cell surface functions.

GATC flanking sequences regulate Dam activity: evidence for how Dam specificity may influence pap expression

Escherichia coli DNA adenine methyltransferase (Dam) plays essential roles in DNA replication, mismatch repair and gene regulation. The differential methylation by Dam of the two GATC sequences in the pap promoter regulates the expression of pili genes necessary for uropathogenic E.coli cellular adhesion. Dam processively methylates GATC sites in various DNA substrates, yet the two pap GATC sites are not processively methylated. We previously proposed that the flanking sequences surrounding the two pap GATC sites contribute to the enzyme's distributive methylation. We show here that replacement of the poorly methylated pap GATC sites with sites predicted to be processively methylated indeed results in an increase in Dam processivity. The increased processivity is due to a change in the methyltransfer kinetics and not the binding efficiency of Dam. A competition experiment in which the flanking sequences of only one pap GATC site were altered demonstrates that the GATC flanking sequences directly regulate the enzyme's catalytic efficiency. The GATC flanking sequences in Dam-regulated promoters in E.coli and other bacteria are similar to those in the pap promoter. Gene regulation from some of these promoters involves mechanisms and proteins that are quite different from those in the pap operon. Further, GATC sequences previously identified to remain unmethylated within the E.coli genome, but whose function remains largely unassigned, are flanked by sequences predicted to be poorly methylated. We conclude that the GATC flanking sequences may be critical for expression of pap and other Dam-regulated genes by affecting the activity of Dam at such sites and, thus, its processivity. A model is proposed, illustrating how the sequences flanking the GATC sites in Dam-regulated promoters may contribute to this epigenetic mechanism of gene expression, and how flanking sequences contribute to the diverse biological roles of Dam.

Rewiring of the Yeast Transcriptional Network Through the Evolution of Motif Usage

Recent experiments revealed large-scale differences in the transcription programs of related species, yet little is known about the genetic basis underlying the evolution of gene expression and its contribution to phenotypic diversity. Here we describe a large-scale modulation of the yeast transcription program that is connected to the emergence of the capacity for rapid anaerobic growth. Genes coding for mitochondrial and cytoplasmic ribosomal proteins display a strongly correlated expression pattern in Candida albicans, but this correlation is lost in the fermentative yeast Saccharomyces cerevisiae. We provide evidence that this change in gene expression is connected to the loss of a specific cis-regulatory element from dozens of genes following the apparent whole-genome duplication event. Our results shed new light on the genetic mechanisms underlying the large-scale evolution of transcriptional networks.

LcrV synthesis is altered by DNA adenine methylase overproduction in Yersinia pseudotuberculosis and is required to confer immunity in vaccinated hosts

Yersinia pseudotuberculosis mutants that overproduce the DNA adenine methylase (DamOP Yersinia) are attenuated, confer robust protective immune responses, and synthesize or secrete several Yersinia outer proteins (Yops) under conditions that are nonpermissive for synthesis and secretion in wild-type strains. To understand the molecular basis of immunity elicited by DamOP Yersinia, we investigated the effects of Dam overproduction on the synthesis and localization of a principal Yersinia immunogen, LcrV, a low-calcium-responsive virulence factor involved in Yop synthesis, localization, and suppression of host inflammatory activities. Dam overproduction relaxed the stringent temperature and calcium regulation of LcrV synthesis. Moreover, the LcrV-dependent synthesis and localization of the actin cytotoxin, YopE, were shown to be relaxed in DamOP cells, suggesting that the synthesis and localization of Yops can occur via both LcrV-dependent and -independent mechanisms. Last, the immunity conferred by DamOP Yersinia was strictly dependent on the presence of LcrV, which may result from its role (i) as an immunogen, (ii) as an immunomodulator of host anti-inflammatory activities, or (iii) in the altered synthesis and localization of Yops that could contribute to immunogen repertoire expansion.

Phylogenetic analysis of 5'-noncoding regions from the ABA-responsive rab16/17 gene family of sorghum, maize and rice provides insight into the composition, organization and function of cis-regulatory modules

Phylogenetic analysis of sequences from gene families and homologous genes from species of varying divergence can be used to identify conserved noncoding regulatory elements. In this study, phylogenetic analysis of 5'-noncoding sequences was optimized using rab17, a well-characterized ABA-responsive gene from maize, and five additional rab16/17 homologs from sorghum and rice. Conserved 5'-noncoding sequences among the maize, sorghum, and rice rab16/17 homologs were identified with the aid of the software program FootPrinter and by screening for known transcription-factor-binding sites. Searches for 7 of 8 (7/8)bp sequence matches within aligned 5'-noncoding segments of the rab genes identified many of the cis-elements previously characterized by biochemical analysis in maize rab17 plus several additional putative regulatory elements. Differences in the composition of conserved noncoding sequences among rab16/17 genes were related to variation in rab gene mRNA levels in different tissues and to response to ABA treatment using qRT-PCR. Absence of a GRA-like element in the promoter of sorghum dhn2 relative to maize rab17 was correlated with an approximately 85-fold reduction of dhn2 RNA in sorghum shoots. Overall, we conclude that phylogenetic analysis of gene families among rice, sorghum, and maize will help identify regulatory sequences in the noncoding regions of genes and contribute to our understanding of grass gene regulatory networks.

PEI - a potent, but not harmless, mucosal immuno-stimulator of mixed T-helper cell response and FasL-mediated cell death in mice

Polyethyleneimine (PEI) is one of the most effective gene delivery systems available today. However, very little is known about its ability to stimulate a systemic immune response and the molecular mechanisms thereof. However, this information is vital for the future development of new gene delivery systems. Here we address this issue by studying gene expression profiles from spleen lymphocytes after in vivo immunization of mice with PEI formulated with a reporter plasmid (PEI+) or the formulation alone (PEI-). PEI- was found to provoke the activation of genes with important immunostimulatory functions, but without the necessary costimulatory signals. PEI+ resulted in: a mixed Th1/Th2 response; activation of both CD8(+) and CD4(+) T cells, with a larger effect on CD4(+); and FasL-mediated antigen-induced cell death. A comparison of the immune responses of PEI+ with that of the clinically used tetanus toxoid-aluminum phosphate vaccine showed that the DNA vaccine provoked a stronger immune response as compared to the protein vaccine. However, many genes involved in other cellular responses such as apoptosis, stress responses and oncogenesis were activated in PEI+, supporting the theory of immunostimulation by danger genes, but also pointing toward possible adverse reactions such as Alzheimer's disease.

Role of SeqA and Dam in Escherichia coli gene expression: a global/microarray analysis

High-density oligonucleotide arrays were used to monitor global transcription patterns in Escherichia coli with various levels of Dam and SeqA proteins. Cells lacking Dam methyltransferase showed a modest increase in transcription of the genes belonging to the SOS regulon. Bacteria devoid of the SeqA protein, which preferentially binds hemimethylated DNA, were found to have a transcriptional profile almost identical to WT bacteria overexpressing Dam methyltransferase. The latter two strains differed from WT in two ways. First, the origin proximal genes were transcribed with increased frequency due to increased gene dosage. Second, chromosomal domains of high transcriptional activity alternate with regions of low activity, and our results indicate that the activity in each domain is modulated in the same way by SeqA deficiency or Dam overproduction. We suggest that the methylation status of the cell is an important factor in forming and/or maintaining chromosome structure.

Genome-wide analysis of deoxyadenosine methyltransferase-mediated control of gene expression in Escherichia coli

Deoxyadenosine methyltransferase (Dam) methylates the deoxyadenine residues in 5'-GATC-3' sequences and is important in many cellular processes in Escherichia coli. We performed a computational analysis of the entire E. coli genome and confirmed that GATC sequences are distributed unevenly in regulatory regions, which suggests that Dam might regulate gene transcription. To test this, a high-density DNA microarray of 4097 E. coli genes was constructed and used to assess the gene expression profiles of the wild type and the dam-16::kam mutant strain grown under four different conditions. We also used two-dimensional electrophoretic analysis of the proteome to assess the protein profiles. The expression of a large number of genes was affected by the dam deficiency. Genes involved in aerobic respiration, stress and SOS responses, amino acid metabolism and nucleotide metabolism were expressed at higher levels in the mutant cells, especially in aerobic conditions. In contrast, transcription of genes participating in anaerobic respiration, flagella biosynthesis, chemotaxis and motility was decreased in the dam mutant strain under both aerobic and low aerobic conditions. Thus, Dam-controlled genes are involved in adjusting the metabolic and respiratory pathways and bacterial motility to suit particular environmental conditions. The promoters of most of these Dam-controlled genes were also found to contain GATC sequences that overlap with recognition sites for two global regulators, fumarate nitrate reduction (Fnr) and catabolite activator protein (CRP). We propose that Dam-mediated methylation plays an important role in the global regulation of genes, particularly those with Fnr and CRP binding sites.

DNA adenine methylase is essential for viability and plays a role in the pathogenesis of Yersinia pseudotuberculosis and Vibrio cholerae

Salmonella strains that lack or overproduce DNA adenine methylase (Dam) elicit a protective immune response to different Salmonella species. To generate vaccines against other bacterial pathogens, the dam genes of Yersinia pseudotuberculosis and Vibrio cholerae were disrupted but found to be essential for viability. Overproduction of Dam significantly attenuated the virulence of these two pathogens, leading to, in Yersinia, the ectopic secretion of virulence proteins (Yersinia outer proteins) and a fully protective immune response in vaccinated hosts. Dysregulation of Dam activity may provide a means for the development of vaccines against varied bacterial pathogens.

From topology to dynamics in biochemical networks

Abstract formulations of the regulation of gene expression as random Boolean switching networks have been studied extensively over the past three decades. These models have been developed to make statistical predictions of the types of dynamics observed in biological networks based on network topology and interaction bias, p. For values of mean connectivity chosen to correspond to real biological networks, these models predict disordered dynamics. However, chaotic dynamics seems to be absent from the functioning of a normal cell. While these models use a fixed number of inputs for each element in the network, recent experimental evidence suggests that several biological networks have distributions in connectivity. We therefore study randomly constructed Boolean networks with distributions in the number of inputs, K, to each element. We study three distributions: delta function, Poisson, and power law (scale free). We analytically show that the critical value of the interaction bias parameter, p, above which steady state behavior is observed, is independent of the distribution in the limit of the number of elements N--> infinity. We also study these networks numerically. Using three different measures (types of attractors, fraction of elements that are active, and length of period), we show that finite, scale-free networks are more ordered than either the Poisson or delta function networks below the critical point. Thus the topology of scale-free biochemical networks, characterized by a wide distribution in the number of inputs per element, may provide a source of order in living cells. (c) 2001 American Institute of Physics.