Phasevarions mediate random switching of gene expression in pathogenic Neisseria

Specificity of the ModA11, ModA12 and ModD1 epigenetic regulator N(6)-adenine DNA methyltransferases of Neisseria meningitidis

Phase variation (random ON/OFF switching) of gene expression is a common feature of host-adapted pathogenic bacteria. Phase variably expressed N(6)-adenine DNA methyltransferases (Mod) alter global methylation patterns resulting in changes in gene expression. These systems constitute phase variable regulons called phasevarions. Neisseria meningitidis phasevarions regulate genes including virulence factors and vaccine candidates, and alter phenotypes including antibiotic resistance. The target site recognized by these Type III N(6)-adenine DNA methyltransferases is not known. Single molecule, real-time (SMRT) methylome analysis was used to identify the recognition site for three key N. meningitidis methyltransferases: ModA11 (exemplified by M.NmeMC58I) (5'-CGY M6A: G-3'), ModA12 (exemplified by M.Nme77I, M.Nme18I and M.Nme579II) (5'-AC M6A: CC-3') and ModD1 (exemplified by M.Nme579I) (5'-CC M6A: GC-3'). Restriction inhibition assays and mutagenesis confirmed the SMRT methylome analysis. The ModA11 site is complex and atypical and is dependent on the type of pyrimidine at the central position, in combination with the bases flanking the core recognition sequence 5'-CGY M6A: G-3'. The observed efficiency of methylation in the modA11 strain (MC58) genome ranged from 4.6% at 5'-GCGC M6A: GG-3' sites, to 100% at 5'-ACGT M6A: GG-3' sites. Analysis of the distribution of modified sites in the respective genomes shows many cases of association with intergenic regions of genes with altered expression due to phasevarion switching.

Selection and Counterselection of Hia Expression Reveals a Key Role for Phase-Variable Expression of Hia in Infection Caused by Nontypeable Haemophilus influenzae

Hia is a major adhesin of nontypeable Haemophilus influenzae (NTHi) and has long been investigated as a vaccine candidate. Here we show that Hia phase variation is controlled by changes in the length of a polythymidine tract located in the hia promoter. Studies of an invasive clinical isolate (strain R2866) show that strains expressing high Hia levels are more efficiently killed by opsonophagocytosis. An opsonophagocytic assay was used to select for a subpopulation of variants that expressed a low level of Hia, which facilitated their escape from killing by anti-Hia antisera. Conversely, a subpopulation of variants expressing a high level of Hia was selected for during passaging through Chang cells. In both cases, phase variation of Hia expression corresponded directly with discrete modal changes in polythymidine tract length. In the chinchilla model of NTHi infection, we observed consistent selection for high Hia expression upon nasopharyngeal colonization, confirming the key role of phase-variable expression of Hia within a specific niche in vivo.

A comparative analysis of methylome profiles of Campylobacter jejuni sheep abortion isolate and gastroenteric strains using PacBio data

Campylobacter jejuni is a leading cause of human gastrointestinal disease and small ruminant abortions in the United States. The recent emergence of a highly virulent, tetracycline-resistant C. jejuni subsp. jejuni sheep abortion clone (clone SA) in the United States, and that strain's association with human disease, has resulted in a heightened awareness of the zoonotic potential of this organism. Pacific Biosciences' Single Molecule, Real-Time sequencing technology was used to explore the variation in the genome-wide methylation patterns of the abortifacient clone SA (IA3902) and phenotypically distinct gastrointestinal-specific C. jejuni strains (NCTC 11168 and 81-176). Several notable differences were discovered that distinguished the methylome of IA3902 from that of 11168 and 81-176: identification of motifs novel to IA3902, genome-specific hypo- and hypermethylated regions, strain level variability in genes methylated, and differences in the types of methylation motifs present in each strain. These observations suggest a possible role of methylation in the contrasting disease presentations of these three C. jejuni strains. In addition, the methylation profiles between IA3902 and a luxS mutant were explored to determine if variations in methylation patterns could be identified that might explain the role of LuxS-dependent methyl recycling in IA3902 abortifacient potential.

qDNAmod: a statistical model-based tool to reveal intercellular heterogeneity of DNA modification from SMRT sequencing data

In an isogenic cell population, phenotypic heterogeneity among individual cells is common and critical for survival of the population under different environment conditions. DNA modification is an important epigenetic factor that can regulate phenotypic heterogeneity. The single molecule real-time (SMRT) sequencing technology provides a unique platform for detecting a wide range of DNA modifications, including N6-methyladenine (6-mA), N4-methylcytosine (4-mC) and 5-methylcytosine (5-mC). Here we present qDNAmod, a novel bioinformatic tool for genome-wide quantitative profiling of intercellular heterogeneity of DNA modification from SMRT sequencing data. It is capable of estimating proportion of isogenic haploid cells, in which the same loci of the genome are differentially modified. We tested the reliability of qDNAmod with the SMRT sequencing data of Streptococcus pneumoniae strain ST556. qDNAmod detected extensive intercellular heterogeneity of DNA methylation (6-mA) in a clonal population of ST556. Subsequent biochemical analyses revealed that the recognition sequences of two type I restriction–modification (R-M) systems are responsible for the intercellular heterogeneity of DNA methylation initially identified by qDNAmod. qDNAmod thus represents a valuable tool for studying intercellular phenotypic heterogeneity from genome-wide DNA modification.

The genetics of Neisseria species

Neisseria gonorrhoeae and Neisseria meningitidis are closely related organisms that cause the sexually transmitted infection gonorrhea and serious bacterial meningitis and septicemia, respectively. Both species possess multiple mechanisms to alter the expression of surface-exposed proteins through the processes of phase and antigenic variation. This potential for wide variability in surface-exposed structures allows the organisms to always have subpopulations of divergent antigenic types to avoid immune surveillance and to contribute to functional variation. Additionally, the Neisseria are naturally competent for DNA transformation, which is their main means of genetic exchange. Although bacteriophages and plasmids are present in this genus, they are not as effective as DNA transformation for horizontal genetic exchange. There are barriers to genetic transfer, such as restriction-modification systems and CRISPR loci, that limit particular types of exchange. These host-restricted pathogens illustrate the rich complexity of genetics that can help define the similarities and differences of closely related organisms.

ModM DNA methyltransferase methylome analysis reveals a potential role for Moraxella catarrhalis phasevarions in otitis media

Moraxella catarrhalis is a significant cause of otitis media and exacerbations of chronic obstructive pulmonary disease. Here, we characterize a phase-variable DNA methyltransferase (ModM), which contains 5'-CAAC-3' repeats in its open reading frame that mediate high-frequency mutation resulting in reversible on/off switching of ModM expression. Three modM alleles have been identified (modM1-3), with modM2 being the most commonly found allele. Using single-molecule, real-time (SMRT) genome sequencing and methylome analysis, we have determined that the ModM2 methylation target is 5'-GAR(m6)AC-3', and 100% of these sites are methylated in the genome of the M. catarrhalis 25239 ModM2 on strain. Proteomic analysis of ModM2 on and off variants revealed that ModM2 regulates expression of multiple genes that have potential roles in colonization, infection, and protection against host defenses. Investigation of the distribution of modM alleles in a panel of M. catarrhalis strains, isolated from the nasopharynx of healthy children or middle ear effusions from patients with otitis media, revealed a statistically significant association of modM3 with otitis media isolates. The modulation of gene expression via the ModM phase-variable regulon (phasevarion), and the significant association of the modM3 allele with otitis media, suggests a key role for ModM phasevarions in the pathogenesis of this organism.

A high-throughput method to examine protein-nucleotide interactions identifies targets of the bacterial transcriptional regulatory protein fur

The Ferric uptake regulatory protein (Fur) is a transcriptional regulatory protein that functions to control gene transcription in response to iron in a number of pathogenic bacteria. In this study, we applied a label-free, quantitative and high-throughput analysis method, Interferometric Reflectance Imaging Sensor (IRIS), to rapidly characterize Fur-DNA interactions in vitro with predicted Fur binding sequences in the genome of Neisseria gonorrhoeae, the causative agent of the sexually transmitted disease gonorrhea. IRIS can easily be applied to examine multiple protein-protein, protein-nucleotide and nucleotide-nucleotide complexes simultaneously and demonstrated here that seventy percent of the predicted Fur boxes in promoter regions of iron-induced genes bound to Fur in vitro with a range of affinities as observed using this microarray screening technology. Combining binding data with mRNA expression levels in a gonococcal fur mutant strain allowed us to identify five new gonococcal genes under Fur-mediated direct regulation.

Phasevarions mediate epigenetic regulation of antimicrobial susceptibility in Neisseria meningitidis

Phase variation is a common feature of host-adapted bacterial pathogens such as Neisseria meningitidis. Recently, we reported that this rapid on/off switching of gene expression occurs in DNA methyltransferases, altering expression in multiple genes via changes in global methylation. In the current study, we compared MIC values of strains with ModA11, ModA12, and ModD1 phasevarions, revealing MIC differences due to ModA11 and ModA12 switching, with a ModA11_OFF strain showing 4-fold reduced susceptibilities to ceftazidime and ciprofloxacin.

Type III restriction-modification enzymes: a historical perspective

Restriction endonucleases interact with DNA at specific sites leading to cleavage of DNA. Bacterial DNA is protected from restriction endonuclease cleavage by modifying the DNA using a DNA methyltransferase. Based on their molecular structure, sequence recognition, cleavage position and cofactor requirements, restriction-modification (R-M) systems are classified into four groups. Type III R-M enzymes need to interact with two separate unmethylated DNA sequences in inversely repeated head-to-head orientations for efficient cleavage to occur at a defined location (25-27 bp downstream of one of the recognition sites). Like the Type I R-M enzymes, Type III R-M enzymes possess a sequence-specific ATPase activity for DNA cleavage. ATP hydrolysis is required for the long-distance communication between the sites before cleavage. Different models, based on 1D diffusion and/or 3D-DNA looping, exist to explain how the long-distance interaction between the two recognition sites takes place. Type III R-M systems are found in most sequenced bacteria. Genome sequencing of many pathogenic bacteria also shows the presence of a number of phase-variable Type III R-M systems, which play a role in virulence. A growing number of these enzymes are being subjected to biochemical and genetic studies, which, when combined with ongoing structural analyses, promise to provide details for mechanisms of DNA recognition and catalysis.

The complex methylome of the human gastric pathogen Helicobacter pylori

The genome of Helicobacter pylori is remarkable for its large number of restriction-modification (R-M) systems, and strain-specific diversity in R-M systems has been suggested to limit natural transformation, the major driving force of genetic diversification in H. pylori. We have determined the comprehensive methylomes of two H. pylori strains at single base resolution, using Single Molecule Real-Time (SMRT®) sequencing. For strains 26695 and J99-R3, 17 and 22 methylated sequence motifs were identified, respectively. For most motifs, almost all sites occurring in the genome were detected as methylated. Twelve novel methylation patterns corresponding to nine recognition sequences were detected (26695, 3; J99-R3, 6). Functional inactivation, correction of frameshifts as well as cloning and expression of candidate methyltransferases (MTases) permitted not only the functional characterization of multiple, yet undescribed, MTases, but also revealed novel features of both Type I and Type II R-M systems, including frameshift-mediated changes of sequence specificity and the interaction of one MTase with two alternative specificity subunits resulting in different methylation patterns. The methylomes of these well-characterized H. pylori strains will provide a valuable resource for future studies investigating the role of H. pylori R-M systems in limiting transformation as well as in gene regulation and host interaction.

The role of variable DNA tandem repeats in bacterial adaptation

DNA tandem repeats (TRs), also designated as satellite DNA, are inter- or intragenic nucleotide sequences that are repeated two or more times in a head-to-tail manner. Because TR tracts are prone to strand-slippage replication and recombination events that cause the TR copy number to increase or decrease, loci containing TRs are hypermutable. An increasing number of examples illustrate that bacteria can exploit this instability of TRs to reversibly shut down or modulate the function of specific genes, allowing them to adapt to changing environments on short evolutionary time scales without an increased overall mutation rate. In this review, we discuss the prevalence and distribution of inter- and intragenic TRs in bacteria and the mechanisms of their instability. In addition, we review evidence demonstrating a role of TR variations in bacterial adaptation strategies, ranging from immune evasion and tissue tropism to the modulation of environmental stress tolerance. Nevertheless, while bioinformatic analysis reveals that most bacterial genomes contain a few up to several dozens of intra- and intergenic TRs, only a small fraction of these have been functionally studied to date.

Phase variation leads to the misidentification of a Neisseria gonorrhoeae virulence gene

Neisseria gonorrhoeae is the causative agent of gonorrhea and an obligate pathogen of humans. The Opa proteins of these bacteria are known to mediate attachment and internalization by host cells, including neutrophils. The Opa protein repertoire of a typical N. gonorrhoeae isolate is encoded on ~11 genes distributed throughout the chromosome and is subject to stochastic changes in expression through phase variation. Together, these characteristics make Opa proteins a critical yet unpredictable aspect of any experimental investigation into the interaction of N. gonorrhoeae with host cells. The goal of this study was to identify novel virulence factors of N. gonorrhoeae by assessing the contribution of a set of uncharacterized hydrogen peroxide-induced genes to bacterial survival against neutrophil-mediated killing. To this end, a strain harboring an engineered mutation in the NGO0322 gene was identified that exhibited increased sensitivity to neutrophil-mediated killing, enhanced internalization by neutrophils, and the ability to induce high levels of neutrophil-generated reactive oxygen species. Each of these phenotypes reverted to near wild-type levels following genetic complementation of the NGO0322 mutation. However, after immunoblot analysis of Opa proteins expressed by the isogenic parent, mutant, and genetically complemented strains, it was determined that phase variation had resulted in a disparity between the Opa profiles of these strains. To determine whether Opa phase variation, rather than NGO0322 mutation, was the cause of the observed neutrophil-related phenotypes, NGO0322 function was investigated in N. gonorrhoeae strains lacking all Opa proteins or constitutively expressing the OpaD variant. In both cases, mutation of NGO0322 did not alter survival of gonococci in the presence of neutrophils. These results demonstrate the importance of controlling for the frequent and random variation in Opa protein production by N. gonorrhoeae when investigating host cell interactions.

Dual pili post-translational modifications synergize to mediate meningococcal adherence to platelet activating factor receptor on human airway cells

Pili of pathogenic Neisseria are major virulence factors associated with adhesion, twitching motility, auto-aggregation, and DNA transformation. Pili of N. meningitidis are subject to several different post-translational modifications. Among these pilin modifications, the presence of phosphorylcholine (ChoP) and a glycan on the pilin protein are phase-variable (subject to high frequency, reversible on/off switching of expression). In this study we report the location of two ChoP modifications on the C-terminus of N. meningitidis pilin. We show that the surface accessibility of ChoP on pili is affected by phase variable changes to the structure of the pilin-linked glycan. We identify for the first time that the platelet activating factor receptor (PAFr) is a key, early event receptor for meningococcal adherence to human bronchial epithelial cells and tissue, and that synergy between the pilin-linked glycan and ChoP post-translational modifications is required for pili to optimally engage PAFr to mediate adherence to human airway cells.

Neisseria meningitidis is an important human pathogen that can cause rapidly progressing, life threatening meningitis and sepsis in humans. There is no fully protective vaccine against this pathogen in current use and the key processes that dictate the transition from harmless carriage of the bacterium in the airway (the case for the vast majority of colonised hosts) to invasive disease are largely undefined. A key missing link in this organism's interaction with the human host is the identity of the receptor that is the first point of contact for the organism within the airway. In this study, we report that the receptor for this important human pathogen on airway epithelial cells is the platelet activating factor receptor (PAFr), an immunomodulatory molecule shown by others to play a role in promoting bacterial sepsis. We also show that two post-translational modifications, glycosylation and phosphorylcholine, are subject to phase-variation (high frequency, reversible switching of gene expression). They are closely associated on adjacent pilin subunits, and synergy between both are required for the efficient engagement with the PAFr. These data define a new role for these post-translational modifications in meningococcal adherence and also provide an insight into the selective pressures that underlie their phase variable expression.

Programmed heterogeneity: epigenetic mechanisms in bacteria

Contrary to the traditional view that bacterial populations are clonal, single-cell analysis reveals that phenotypic heterogeneity is common in bacteria. Formation of distinct bacterial lineages appears to be frequent during adaptation to harsh environments, including the colonization of animals by bacterial pathogens. Formation of bacterial subpopulations is often controlled by epigenetic mechanisms that generate inheritable phenotypic diversity without altering the DNA sequence. Such mechanisms are diverse, ranging from relatively simple feedback loops to complex self-perpetuating DNA methylation patterns.

Role of DNA methyltransferases in epigenetic regulation in bacteria

In prokaryotes, alteration in gene expression was observed with the modification of DNA, especially DNA methylation. Such changes are inherited from generation to generation with no alterations in the DNA sequence and represent the epigenetic signal in prokaryotes. DNA methyltransferases are enzymes involved in DNA modification and thus in epigenetic regulation of gene expression. DNA methylation not only affects the thermodynamic stability of DNA, but also changes its curvature. Methylation of specific residues on DNA can affect the protein-DNA interactions. DNA methylation in prokaryotes regulates a number of physiological processes in the bacterial cell including transcription, DNA mismatch repair and replication initiation. Significantly, many reports have suggested a role of DNA methylation in regulating the expression of a number of genes in virulence and pathogenesis thus, making DNA methlytransferases novel targets for the designing of therapeutics. Here, we summarize the current knowledge about the influence of DNA methylation on gene regulation in different bacteria, and on bacterial virulence.

Homologous recombination drives both sequence diversity and gene content variation in Neisseria meningitidis

The study of genetic and phenotypic variation is fundamental for understanding the dynamics of bacterial genome evolution and untangling the evolution and epidemiology of bacterial pathogens. Neisseria meningitidis (Nm) is among the most intriguing bacterial pathogens in genomic studies due to its dynamic population structure and complex forms of pathogenicity. Extensive genomic variation within identical clonal complexes (CCs) in Nm has been recently reported and suggested to be the result of homologous recombination, but the extent to which recombination contributes to genomic variation within identical CCs has remained unclear. In this study, we sequenced two Nm strains of identical serogroup (C) and multi-locus sequence type (ST60), and conducted a systematic analysis with an additional 34 Nm genomes. Our results revealed that all gene content variation between the two ST60 genomes was introduced by homologous recombination at the conserved flanking genes, and 94.25% or more of sequence divergence was caused by homologous recombination. Recombination was found in genes associated with virulence factors, antigenic outer membrane proteins, and vaccine targets, suggesting an important role of homologous recombination in rapidly altering the pathogenicity and antigenicity of Nm. Recombination was also evident in genes of the restriction and modification systems, which may undermine barriers to DNA exchange. In conclusion, homologous recombination can drive both gene content variation and sequence divergence in Nm. These findings shed new light on the understanding of the rapid pathoadaptive evolution of Nm and other recombinogenic bacterial pathogens.

LuxS influences Escherichia coli biofilm formation through autoinducer-2-dependent and autoinducer-2-independent modalities

Escherichia coli produces biofilms in response to the small molecule autoinducer-2 (AI-2), a product of the LuxS enzyme. LuxS is part of the activated methyl cycle and could also affect biofilm development by AI-2-independent effects on metabolism. A luxS deletion mutant of E. coli W3110 and an inducible plasmid-luxS-complemented strain were used to identify AI-2-independent phenotypes. Differential interference contrast microscopy revealed distinct surface colonization patterns. Confocal microscopy followed by quantitative image analysis determined differences in biofilm topography correlating with luxS expression; deletion mutant biofilms had a 'spreading' phenotype, whereas the complement had a 'climbing' phenotype. Addition of exogenous 4,5-dihydroxy-2,3-pentanedione (DPD), an AI-2 precursor, to the deletion mutant increased biofilm height and biomass, whereas addition of the methyl donor S-adenosyl methionine or aspartate prevented the luxS-complemented strain from producing a thick biofilm. The luxS-complemented strain autoaggregated, indicating that fimbriae production was inhibited, which was confirmed by transmission electron microscopy. DPD could not induce autoaggregation in the deletion mutant, demonstrating that fimbriation was an AI-2-independent phenotype. Carbon utilization was affected by LuxS, potentially contributing to the observed phenotypic differences. Overall, the work demonstrated that LuxS affected E. coli biofilm formation independently of AI-2 and could assist in adapting to diverse conditions.

Genome-wide mapping of methylated adenine residues in pathogenic Escherichia coli using single-molecule real-time sequencing

Single-molecule real-time (SMRT) DNA sequencing allows the systematic detection of chemical modifications such as methylation but has not previously been applied on a genome-wide scale. We used this approach to detect 49,311 putative 6-methyladenine (m6A) residues and 1,407 putative 5-methylcytosine (m5C) residues in the genome of a pathogenic Escherichia coli strain. We obtained strand-specific information for methylation sites and a quantitative assessment of the frequency of methylation at each modified position. We deduced the sequence motifs recognized by the methyltransferase enzymes present in this strain without prior knowledge of their specificity. Furthermore, we found that deletion of a phage-encoded methyltransferase-endonuclease (restriction-modification; RM) system induced global transcriptional changes and led to gene amplification, suggesting that the role of RM systems extends beyond protecting host genomes from foreign DNA.

Beyond junk-variable tandem repeats as facilitators of rapid evolution of regulatory and coding sequences

Copy Number Variations (CNVs) and Single Nucleotide Polymorphisms (SNPs) have been the major focus of most large-scale comparative genomics studies to date. Here, we discuss a third, largely ignored, type of genetic variation, namely changes in tandem repeat number. Historically, tandem repeats have been designated as non functional “junk” DNA, mostly as a result of their highly unstable nature. With the exception of tandem repeats involved in human neurodegenerative diseases, repeat variation was often believed to be neutral with no phenotypic consequences. Recent studies, however, have shown that as many as 10% to 20% of coding and regulatory sequences in eukaryotes contain an unstable repeat tract. Contrary to initial suggestions, tandem repeat variation can have useful phenotypic consequences. Examples include rapid variation in microbial cell surface, tuning of internal molecular clocks in flies and the dynamic morphological plasticity in mammals. As such, tandem repeats can be useful functional elements that facilitate evolvability and rapid adaptation.

Movement of DNA sequence recognition domains between non-orthologous proteins

Comparisons of proteins show that they evolve through the movement of domains. However, in many cases, the underlying mechanisms remain unclear. Here, we observed the movements of DNA recognition domains between non-orthologous proteins within a prokaryote genome. Restriction-modification (RM) systems, consisting of a sequence-specific DNA methyltransferase and a restriction enzyme, contribute to maintenance/evolution of genomes/epigenomes. RM systems limit horizontal gene transfer but are themselves mobile. We compared Type III RM systems in Helicobacter pylori genomes and found that target recognition domain (TRD) sequences are mobile, moving between different orthologous groups that occupy unique chromosomal locations. Sequence comparisons suggested that a likely underlying mechanism is movement through homologous recombination of similar DNA sequences that encode amino acid sequence motifs that are conserved among Type III DNA methyltransferases. Consistent with this movement, incongruence was observed between the phylogenetic trees of TRD regions and other regions in proteins. Horizontal acquisition of diverse TRD sequences was suggested by detection of homologs in other Helicobacter species and distantly related bacterial species. One of these RM systems in H. pylori was inactivated by insertion of another RM system that likely transferred from an oral bacterium. TRD movement represents a novel route for diversification of DNA-interacting proteins.

Restriction and sequence alterations affect DNA uptake sequence-dependent transformation in Neisseria meningitidis

Transformation is a complex process that involves several interactions from the binding and uptake of naked DNA to homologous recombination. Some actions affect transformation favourably whereas others act to limit it. Here, meticulous manipulation of a single type of transforming DNA allowed for quantifying the impact of three different mediators of meningococcal transformation: NlaIV restriction, homologous recombination and the DNA Uptake Sequence (DUS). In the wildtype, an inverse relationship between the transformation frequency and the number of NlaIV restriction sites in DNA was observed when the transforming DNA harboured a heterologous region for selection (ermC) but not when the transforming DNA was homologous with only a single nucleotide heterology. The influence of homologous sequence in transforming DNA was further studied using plasmids with a small interruption or larger deletions in the recombinogenic region and these alterations were found to impair transformation frequency. In contrast, a particularly potent positive driver of DNA uptake in Neisseria sp. are short DUS in the transforming DNA. However, the molecular mechanism(s) responsible for DUS specificity remains unknown. Increasing the number of DUS in the transforming DNA was here shown to exert a positive effect on transformation. Furthermore, an influence of variable placement of DUS relative to the homologous region in the donor DNA was documented for the first time. No effect of altering the orientation of DUS was observed. These observations suggest that DUS is important at an early stage in the recognition of DNA, but does not exclude the existence of more than one level of DUS specificity in the sequence of events that constitute transformation. New knowledge on the positive and negative drivers of transformation may in a larger perspective illuminate both the mechanisms and the evolutionary role(s) of one of the most conserved mechanisms in nature: homologous recombination.

Origin of the diversity in DNA recognition domains in phasevarion associated modA genes of pathogenic Neisseria and Haemophilus influenzae

Phase variable restriction-modification (R-M) systems have been identified in a range of pathogenic bacteria. In some it has been demonstrated that the random switching of the mod (DNA methyltransferase) gene mediates the coordinated expression of multiple genes and constitutes a phasevarion (phase variable regulon). ModA of Neisseria and Haemophilus influenzae contain a highly variable, DNA recognition domain (DRD) that defines the target sequence that is modified by methylation and is used to define modA alleles. 18 distinct modA alleles have been identified in H. influenzae and the pathogenic Neisseria. To determine the origin of DRD variability, the 18 modA DRDs were used to search the available databases for similar sequences. Significant matches were identified between several modA alleles and mod gene from distinct bacterial species, indicating one source of the DRD variability was via horizontal gene transfer. Comparison of DRD sequences revealed significant mosaicism, indicating exchange between the Neisseria and H. influenzae modA alleles. Regions of high inter- and intra-allele similarity indicate that some modA alleles had undergone recombination more frequently than others, generating further diversity. Furthermore, the DRD from some modA alleles, such as modA12, have been transferred en bloc to replace the DRD from different modA alleles.

Protocol for gene expression profiling using DNA microarrays in Neisseria gonorrhoeae

Gene expression profiling using DNA microarrays has become commonplace in current molecular biology practices, and has dramatically enhanced our understanding of the biology of Neisseria spp., and the interaction of these organisms with the host. With the choice of microarray platforms offered for gene expression profiling and commercially available arrays, investigators must ask several central questions to make decisions based on their research focus. Are arrays on hand for their organism and if not then would it be cost-effective to design custom arrays. Other important considerations; what types of specialized equipment for array hybridization and signal detection are required and is the specificity and sensitivity of the array adequate for your application. Here, we describe the use of a custom 12K CombiMatrix ElectraSense™ oligonucleotide microarray format for assessing global gene expression profiles in Neisseria spp.

Phasevarion mediated epigenetic gene regulation in Helicobacter pylori

Many host-adapted bacterial pathogens contain DNA methyltransferases (mod genes) that are subject to phase-variable expression (high-frequency reversible ON/OFF switching of gene expression). In Haemophilus influenzae and pathogenic Neisseria, the random switching of the modA gene, associated with a phase-variable type III restriction modification (R-M) system, controls expression of a phase-variable regulon of genes (a “phasevarion”), via differential methylation of the genome in the modA ON and OFF states. Phase-variable type III R-M systems are also found in Helicobacter pylori, suggesting that phasevarions may also exist in this key human pathogen. Phylogenetic studies on the phase-variable type III modH gene revealed that there are 17 distinct alleles in H. pylori, which differ only in their DNA recognition domain. One of the most commonly found alleles was modH5 (16% of isolates). Microarray analysis comparing the wild-type P12modH5 ON strain to a P12ΔmodH5 mutant revealed that six genes were either up- or down-regulated, and some were virulence-associated. These included flaA, which encodes a flagella protein important in motility and hopG, an outer membrane protein essential for colonization and associated with gastric cancer. This study provides the first evidence of this epigenetic mechanism of gene expression in H. pylori. Characterisation of H. pylori modH phasevarions to define stable immunological targets will be essential for vaccine development and may also contribute to understanding H. pylori pathogenesis.

Deletion of one nucleotide within the homonucleotide tract present in the hsdS gene alters the DNA sequence specificity of type I restriction-modification system NgoAV

As a result of a frameshift mutation, the hsdS locus of the NgoAV type IC restriction and modification (RM) system comprises two genes, hsdS(NgoAV1) and hsdS(NgoAV2). The specificity subunit, HsdS(NgoAV), the product of the hsdS(NgoAV1) gene, is a naturally truncated form of an archetypal specificity subunit (208 N-terminal amino acids instead of 410). The presence of a homonucleotide tract of seven guanines (poly[G]) at the 3' end of the hsdS(NgoAV1) gene makes the NgoAV system a strong candidate for phase variation, i.e., stochastic addition or reduction in the guanine number. We have constructed mutants with 6 guanines instead of 7 and demonstrated that the deletion of a single nucleotide within the 3' end of the hsdS(NgoAV1) gene restored the fusion between the hsdS(NgoAV1) and hsdS(NgoAV2) genes. We have demonstrated that such a contraction of the homonucleotide tract may occur in vivo: in a Neisseria gonorrhoeae population, a minor subpopulation of cells appeared to have only 6 guanines at the 3' end of the hsdS(NgoAV1) gene. Escherichia coli cells carrying the fused gene and expressing the NgoAVΔ RM system were able to restrict λ phage at a level comparable to that for the wild-type NgoAV system. NgoAV recognizes the quasipalindromic interrupted sequence 5'-GCA(N(8))TGC-3' and methylates both strands. NgoAVΔ recognizes DNA sequences 5'-GCA(N(7))GTCA-3' and 5'-GCA(N(7))CTCA-3', although the latter sequence is methylated only on the complementary strand within the 5'-CTCA-3' region of the second recognition target sequence.

Detection of genetically modified microorganisms in soil using the most-probable-number method with multiplex PCR and DNA dot blot

The principal objective of this study was to detect genetically modified microorganisms (GMMs) that might be accidentally released into the environment from laboratories. Two methods [plate counting and most-probable-number (MPN)] coupled with either multiplex PCR or DNA dot blots were compared using genetically modified Escherichia coli, Pseudomonas putida, and Acinetobacter oleivorans harboring an antibiotic-resistance gene with additional gfp and lacZ genes as markers. Alignments of sequences collected from databases using the Perl scripting language (Perl API) and from denaturing gradient gel electrophoresis analysis revealed that the gfp, lacZ and antibiotic-resistance genes (kanamycin, tetracycline, and ampicillin) in GMMs differed from the counterpart genes in many sequenced genomes and in soil DNA. Thus, specific multiplex PCR primer sets for detection of plasmid-based gfp and lacZ antibiotic-resistance genes could be generated. In the plate counting method, many antibiotic-resistant bacteria from a soil microcosm grew as colonies on antibiotic-containing agar plates. The multiplex PCR verification of randomly selected antibiotic-resistant colonies with specific primers proved ineffective. The MPN-multiplex PCR method and antibiotic-resistant phenotype could be successfully used to detect GMMs, although this method is quite laborious. The MPN-DNA dot blot method screened more cells at a time in a microtiter plate containing the corresponding antibiotics, and was shown to be a more efficient method for the detection of GMMs in soil using specific probes in terms of labor and accuracy.

The Composition and Metabolic Phenotype of Neisseria gonorrhoeae Biofilms

Neisseria gonorrhoeae has been shown to form biofilms during cervical infection. Thus, biofilm formation may play an important role in the infection of women. The ability of N. gonorrhoeae to form membrane blebs is crucial to biofilm formation. Blebs contain DNA and outer membrane structures, which have been shown to be major constituents of the biofilm matrix. The organism expresses a DNA thermonuclease that is involved in remodeling of the biofilm matrix. Comparison of the transcriptional profiles of gonococcal biofilms and planktonic runoff indicate that genes involved in anaerobic metabolism and oxidative stress tolerance are more highly expressed in biofilm. The expression of aniA, ccp, and norB, which encode nitrite reductase, cytochrome c peroxidase, and nitric oxide reductase respectively, is required for mature biofilm formation over glass and human cervical cells. In addition, anaerobic respiration occurs in the substratum of gonococcal biofilms and disruption of the norB gene required for anaerobic respiration, results in a severe biofilm attenuation phenotype. It has been demonstrated that accumulation of nitric oxide (NO) contributes to the phenotype of a norB mutant and can retard biofilm formation. However, NO can also enhance biofilm formation, and this is largely dependent on the concentration and donation rate or steady-state kinetics of NO. The majority of the genes involved in gonococcal oxidative stress tolerance are also required for normal biofilm formation, as mutations in the following genes result in attenuated biofilm formation over cervical cells and/or glass: oxyR, gor, prx, mntABC, trxB, and estD. Overall, biofilm formation appears to be an adaptation for coping with the environmental stresses present in the female genitourinary tract. Therefore, this review will discuss the studies, which describe the composition and metabolic phenotype of gonococcal biofilms.

Neisseria meningitidis is structured in clades associated with restriction modification systems that modulate homologous recombination

Molecular data on a limited number of chromosomal loci have shown that the population of Neisseria meningitidis (Nm), a deadly human pathogen, is structured in distinct lineages. Given that the Nm population undergoes substantial recombination, the mechanisms resulting in the evolution of these lineages, their persistence in time, and the implications for the pathogenicity of the bacterium are not yet completely understood. Based on whole-genome sequencing, we show that Nm is structured in phylogenetic clades. Through acquisition of specific genes and through insertions and rearrangements, each clade has acquired and remodeled specific genomic tracts, with the potential to impact on the commensal and virulence behavior of Nm. Despite this clear evidence of a structured population, we confirm high rates of detectable recombination throughout the whole Nm chromosome. However, gene conversion events were found to be longer within clades than between clades, suggesting a DNA cleavage mechanism associated with the phylogeny of the species. We identify 22 restriction modification systems, probably acquired by horizontal gene transfer from outside of the species/genus, whose distribution in the different strains coincides with the phylogenetic clade structure. We provide evidence that these clade-associated restriction modification systems generate a differential barrier to DNA exchange consistent with the observed population structure. These findings have general implications for the emergence of lineage structure and virulence in recombining bacterial populations, and they could provide an evolutionary framework for the population biology of a number of other bacterial species that show contradictory population structure and dynamics.

Standing genetic variation in contingency loci drives the rapid adaptation of Campylobacter jejuni to a novel host

The genome of the food-borne pathogen Campylobacter jejuni contains multiple highly mutable sites, or contingency loci. It has been suggested that standing variation at these loci is a mechanism for rapid adaptation to a novel environment, but this phenomenon has not been shown experimentally. In previous work we showed that the virulence of C. jejuni NCTC11168 increased after serial passage through a C57BL/6 IL-10^-/- mouse model of campylobacteriosis. Here we sought to determine the genetic basis of this adaptation during passage. Re-sequencing of the 1.64Mb genome to 200-500X coverage allowed us to define variation in 23 contingency loci to an unprecedented depth both before and after in vivo adaptation. Mutations in the mouse-adapted C. jejuni were largely restricted to the homopolymeric tracts of thirteen contingency loci. These changes cause significant alterations in open reading frames of genes in surface structure biosynthesis loci and in genes with only putative functions. Several loci with open reading frame changes also had altered transcript abundance. The increase in specific phases of contingency loci during in vivo passage of C. jejuni, coupled with the observed virulence increase and the lack of other types of genetic changes, is the first experimental evidence that these variable regions play a significant role in C. jejuni adaptation and virulence in a novel host.

Variable tandem repeats accelerate evolution of coding and regulatory sequences

Genotype-to-phenotype mapping commonly focuses on two major classes of mutations: single nucleotide polymorphisms (SNPs) and copy number variation (CNV). Here, we discuss an underestimated third class of genotypic variation: changes in microsatellite and minisatellite repeats. Such tandem repeats (TRs) are ubiquitous, unstable genomic elements that have historically been designated as nonfunctional "junk DNA" and are therefore mostly ignored in comparative genomics. However, as many as 10% to 20% of eukaryotic genes and promoters contain an unstable repeat tract. Mutations in these repeats often have fascinating phenotypic consequences. For example, changes in unstable repeats located in or near human genes can lead to neurodegenerative diseases such as Huntington disease. Apart from their role in disease, variable repeats also confer useful phenotypic variability, including cell surface variability, plasticity in skeletal morphology, and tuning of the circadian rhythm. As such, TRs combine characteristics of genetic and epigenetic changes that may facilitate organismal evolvability.

Low levels of antibacterial drug resistance expressed by Gram-negative bacteria isolated from poultry carcasses in New Zealand

AIM

To provide baseline data on the levels and patterns of antibacterial drug resistance expressed by Gram-negative bacteria isolated from poultry carcasses in New Zealand.

METHODS

Between July and December 2006, isolates of Escherichia coli (n=407) and Salmonella spp. (n=3) originating from carcass-rinse samples were submitted by testing laboratories affiliated to five major poultry processing plants. Isolates of Campylobacter jejuni (n=193) originating from retail poultry carcasses in 2005-2006 were retrieved from the Massey University archives. All isolates underwent disc diffusion susceptibility testing against panels of 12 (Enterobacteriaceae) and six (Campylobacter spp.) antibacterial drugs. Cephalothin-resistance in isolates of E. coli was confirmed using ETest strips, and confirmation of the resistance phenotypes for a subset of C. jejuni isolates used microbroth dilution assays. Patterns within the resistance phenotypes of the isolates were investigated using hierarchical clustering, and logistic regression modelling.

RESULTS

The majority of isolates (71.5% E. coli, 99% C. jejuni, and all three Salmonella spp. isolates) were fully susceptible to the drugs that were tested. Four (1%) E. coli isolates showed resistance to three or more drugs. The proportions of susceptible E. coli differed between the five processing plants. Resistances were detected in E. coli isolates, using disc diffusion to cephalothin (18.2%), ampicillin (4.4%), tetracycline (4.4%) and gentamicin (1.5%). There was an association between cephalothin-resistant isolates of E. coli and decreased susceptibility to gentamicin. Using ETests to ascertain the minimum inhibitory concentrations (MIC) of E. coli for cephalothin gave inconsistent results. One of 193 C. jejuni isolates was resistant to erythromycin, and microbroth dilution assays confirmed that this panel of C. jejuni was generally susceptible to antibacterial drugs.

CONCLUSIONS

The levels of resistance shown by Gram-negative bacteria isolated from chicken carcasses in New Zealand are among the lowest reported around the world. No resistance to extended-spectrum cephalosporin drugs was detected in E. coli, suggesting that CTX-M and AmpC beta-lactamases are rare or absent. Salmonella spp. are rarely isolated from poultry carcasses during routine testing in New Zealand, and the isolates identified during this study were fully susceptible to the drugs tested. A panel of C. jejuni isolates originating from retail poultry carcasses were susceptible to first-line and second-line antibacterial drugs. The use of cephalothin as a marker of resistance to first-generation cephalosporins may not be appropriate for non-type-specific E. coli of animal origin.

Diversity of DNA methyltransferases that recognize asymmetric target sequences

DNA methyltransferases (MTases) are a group of enzymes that catalyze the methyl group transfer from S-adenosyl-L-methionine in a sequence-specific manner. Orthodox Type II DNA MTases usually recognize palindromic DNA sequences and add a methyl group to the target base (either adenine or cytosine) on both strands. However, there are a number of MTases that recognize asymmetric target sequences and differ in their subunit organization. In a bacterial cell, after each round of replication, the substrate for any MTase is hemimethylated DNA, and it therefore needs only a single methylation event to restore the fully methylated state. This is in consistent with the fact that most of the DNA MTases studied exist as monomers in solution. Multiple lines of evidence suggest that some DNA MTases function as dimers. Further, functional analysis of many restriction-modification systems showed the presence of more than one or fused MTase genes. It was proposed that presence of two MTases responsible for the recognition and methylation of asymmetric sequences would protect the nascent strands generated during DNA replication from cognate restriction endonuclease. In this review, MTases recognizing asymmetric sequences have been grouped into different subgroups based on their unique properties. Detailed characterization of these unusual MTases would help in better understanding of their specific biological roles and mechanisms of action. The rapid progress made by the genome sequencing of bacteria and archaea may accelerate the identification and study of species- and strain-specific MTases of host-adapted bacteria and their roles in pathogenic mechanisms.

The phasevarion: phase variation of type III DNA methyltransferases controls coordinated switching in multiple genes

In several host-adapted pathogens, phase variation has been found to occur in genes that encode methyltransferases associated with type III restriction-modification systems. It was recently shown that in the human pathogens Haemophilus influenzae, Neisseria gonorrhoeae and Neisseria meningitidis phase variation of a type III DNA methyltransferase, encoded by members of the mod gene family, regulates the expression of multiple genes. This novel genetic system has been termed the 'phasevarion' (phase-variable regulon). The wide distribution of phase-variable mod family genes indicates that this may be a common strategy used by host-adapted bacterial pathogens to randomly switch between distinct cell types.

A novel phase variation mechanism in the meningococcus driven by a ligand-responsive repressor and differential spacing of distal promoter elements

Phase variable expression, mediated by high frequency reversible changes in the length of simple sequence repeats, facilitates adaptation of bacterial populations to changing environments and is frequently important in bacterial virulence. Here we elucidate a novel phase variable mechanism for NadA, an adhesin and invasin of Neisseria meningitidis. The NadR repressor protein binds to operators flanking the phase variable tract and contributes to the differential expression levels of phase variant promoters with different numbers of repeats likely due to different spacing between operators. We show that IHF binds between these operators, and may permit looping of the promoter, allowing interaction of NadR at operators located distally or overlapping the promoter. The 4-hydroxyphenylacetic acid, a metabolite of aromatic amino acid catabolism that is secreted in saliva, induces NadA expression by inhibiting the DNA binding activity of the repressor. When induced, only minor differences are evident between NadR-independent transcription levels of promoter phase variants and are likely due to differential RNA polymerase contacts leading to altered promoter activity. Our results suggest that NadA expression is under both stochastic and tight environmental-sensing regulatory control, both mediated by the NadR repressor, and may be induced during colonization of the oropharynx where it plays a major role in the successful adhesion and invasion of the mucosa. Hence, simple sequence repeats in promoter regions may be a strategy used by host-adapted bacterial pathogens to randomly switch between expression states that may nonetheless still be induced by appropriate niche-specific signals.

Diversification strategies, through genetic switches that randomly turn genes on and off, occur in many pathogenic bacterial populations and confer adaptive advantages to new environments and evasion of host immune responses. This is often mediated by spontaneous changes in the length of short DNA sequence repeats located in protein-coding regions or upstream regulatory regions, leading to deactivation or alteration of the associated genes. In this study we describe how a repeat sequence, distally upstream of the promoter region, alters the expression of an important adhesin of N. meningitidis. We identify the major mediator of this control, a negative regulator NadR, which binds to sequences flanking the variable repeat. Changes in the spacing between these sequences affect the ability of NadR to shut down expression from the promoter. We also identify a relevant metabolite that can block NadR activity and therefore act as a signal to induce adhesin expression. This finding sheds new light on the role of DNA-repeats identified in intergenic regions for which no role could be hypothesised, and may be a model mechanism used by bacterial pathogens for fine-tuning diversity within the host. Elucidating these mechanisms can aid in our understanding and prevention of disease.

Characterization of the NgoAXP: phase-variable type III restriction-modification system in Neisseria gonorrhoeae

Methyltransferases associated with type III restriction-modification (RM) systems are phase-variably expressed in a variety of pathogenic bacteria. NgoAXP, the type III RM system encoded by Neisseria gonorrhoeae, was characterized in this study. The cloned resngoAXP and ngoAXPmod genes were expressed in Escherichia coli strains. The restriction and modification activities of NgoAXP were confirmed in vivo by the lambda phage restriction and modification test and in vitro by the methylation of DNA substrates in the presence of [methyl-(3)H]AdoMet. As in all known type III systems, the restriction activity needed the presence of both genes, while the presence of the ngoAXPmod gene was sufficient for DNA methylation. Following its overexpression, the DNA methyltransferase M.NgoAXP was purified to apparent homogeneity using metal affinity chromatography. The specific sequence recognized by this enzyme was determined as a nonpalindromic sequence: 5'-CCACC-3', in which the adenine residue is methylated. We observed that in E. coli cells, the expression of the restriction phenotype associated with NgoAXP switched randomly. This phase variation was associated with the change in the number of pentanucleotide repeats (5'-CCAAC/G-3') present at the 5'-end of the coding region of the ngoAXPmod gene.