Adaptation by phase variation in pathogenic bacteria

Transcriptional and immunological analysis of the putative outer membrane protein and vaccine candidate TprL of Treponema pallidum

Background

An effective syphilis vaccine should elicit antibodies to Treponema pallidum subsp. pallidum (T. p. pallidum) surface antigens to induce pathogen clearance through opsonophagocytosis. Although the combination of bioinformatics, structural, and functional analyses of T. p. pallidum genes to identify putative outer membrane proteins (OMPs) resulted in a list of potential vaccine candidates, still very little is known about whether and how transcription of these genes is regulated during infection. This knowledge gap is a limitation to vaccine design, as immunity generated to an antigen that can be down-regulated or even silenced at the transcriptional level without affecting virulence would not induce clearance of the pathogen, hence allowing disease progression.

Principal findings

We report here that tp1031, the T. p. pallidum gene encoding the putative OMP and vaccine candidate TprL is differentially expressed in several T. p. pallidum strains, suggesting transcriptional regulation. Experimental identification of the tprL transcriptional start site revealed that a homopolymeric G sequence of varying length resides within the tprL promoter and that its length affects promoter activity compatible with phase variation. Conversely, in the closely related pathogen T. p. subsp. pertenue, the agent of yaws, where a naturally-occurring deletion has eliminated the tprL promoter region, elements necessary for protein synthesis, and part of the gene ORF, tprL transcription level are negligible compared to T. p. pallidum strains. Accordingly, the humoral response to TprL is absent in yaws-infected laboratory animals and patients compared to syphilis-infected subjects.

Conclusion

The ability of T. p. pallidum to stochastically vary tprL expression should be considered in any vaccine development effort that includes this antigen. The role of phase variation in contributing to T. p. pallidum antigenic diversity should be further studied.

Syphilis is still an endemic disease in many low- and middle-income countries and has been resurgent in high-income nations for almost two decades now. In endemic areas, syphilis still causes significant morbidity and mortality in patients, particularly when its causative agent, the bacterium Treponema pallidum subsp. pallidum is transmitted to the fetus during pregnancy. Although there are significant ongoing efforts to identify an effective syphilis vaccine to bring into clinical trials within the decade in the U.S., such efforts are partially hindered by the lack of knowledge on transcriptional regulation of many genes encoding vaccine candidates. Here, we start addressing this knowledge gap for the putative outer membrane protein (OMP) and vaccine candidates TprL, encoded by the tp1031 gene. As we previously reported for other putative OMP-encoding genes of the syphilis agent, tprL transcription level appears to be affected by the length of a homopolymeric sequence of guanosines (Gs) located within the gene promoter. This is a mechanism known as phase variation and often involved in altering the surface antigenic profile of a bacterial pathogen to facilitate immune evasion and/or adaptation to the host milieu.

Phase-Variable Glycosylation in Nontypeable Haemophilus influenzae

Nontypeable Haemophilus influenzae (NTHi) is a leading cause of respiratory tract infections worldwide and continues to be a global health burden. Adhesion and colonization of host cells are crucial steps in bacterial pathogenesis, and in many strains of NTHi, the interaction with the host is mediated by the high molecular weight adhesins HMW1A and HMW2A. These adhesins are N-glycoproteins that are modified by cytoplasmic glycosyltransferases HMW1C and HMW2C. Phase variation in the number of short sequence repeats in the promoters of hmw1A and hmw2A directly affects their expression. Here, we report the presence of similar variable repeat elements in the promoters of hmw1C and hmw2C in diverse NTHi isolates. In an ex vivo assay, we systematically altered the substrate and glycosyltransferase expression and showed that both of these factors quantitatively affected the site-specific efficiency of glycosylation on HMW-A. This represents a novel mechanism through which phase variation can generate diversity in the quantitative extent of site-specific post-translational modifications on antigenic surface proteins. Glycosylation occupancy was incomplete at many sites, variable between sites, and generally lower close to the C-terminus of HMW-A. We investigated the causes of this variability. As HMW-C glycosylates HMW-A in the cytoplasm, we tested how secretion affected glycosylation on HMW-A and showed that retaining HMW-A in the cytoplasm indeed increased glycosylation occupancy across the full length of the protein. Site-directed mutagenesis showed that HMW-C had no inherent preference for glycosylating asparagines in NxS or NxT sequons. This work provides key insights into factors contributing to the heterogenous modifications of NTHi HMW-A adhesins, expands knowledge of NTHi population diversity and pathogenic capability, and is relevant to vaccine design for NTHi and related pathogens.

The evolutionary advantage of heritable phenotypic heterogeneity

Phenotypic plasticity is an evolutionary driving force in diverse biological processes, including the adaptive immune system, the development of neoplasms, and the persistence of pathogens despite drug pressure. It is essential, therefore, to understand the evolutionary advantage of an allele that confers on cells the ability to express a range of phenotypes. Here, we study the fate of a new mutation that allows the expression of multiple phenotypic states, introduced into a finite population of individuals that can express only a single phenotype. We show that the advantage of such a mutation depends on the degree of phenotypic heritability between generations, called phenotypic memory. We analyze the fixation probability of the phenotypically plastic allele as a function of phenotypic memory, the variance of expressible phenotypes, the rate of environmental changes, and the population size. We find that the fate of a phenotypically plastic allele depends fundamentally on the environmental regime. In constant environments, plastic alleles are advantageous and their fixation probability increases with the degree of phenotypic memory. In periodically fluctuating environments, by contrast, there is an optimum phenotypic memory that maximizes the probability of the plastic allele's fixation. This same optimum memory also maximizes geometric mean fitness, in steady state. We interpret these results in the context of previous studies in an infinite-population framework. We also discuss the implications of our results for the design of therapies that can overcome persistence and, indirectly, drug resistance.

Transcription of TP0126, Treponema pallidum putative OmpW homolog, is regulated by the length of a homopolymeric guanosine repeat

An effective mechanism for introduction of phenotypic diversity within a bacterial population exploits changes in the length of repetitive DNA elements located within gene promoters. This phenomenon, known as phase variation, causes rapid activation or silencing of gene expression and fosters bacterial adaptation to new or changing environments. Phase variation often occurs in surface-exposed proteins, and in Treponema pallidum subsp. pallidum, the syphilis agent, it was reported to affect transcription of three putative outer membrane protein (OMP)-encoding genes. When the T. pallidum subsp. pallidum Nichols strain genome was initially annotated, the TP0126 open reading frame was predicted to include a poly(G) tract and did not appear to have a predicted signal sequence that might suggest the possibility of its being an OMP. Here we show that the initial annotation was incorrect, that this poly(G) is instead located within the TP0126 promoter, and that it varies in length in vivo during experimental syphilis. Additionally, we show that TP0126 transcription is affected by changes in the poly(G) length consistent with regulation by phase variation. In silico analysis of the TP0126 open reading frame based on the experimentally identified transcriptional start site shortens this hypothetical protein by 69 amino acids, reveals a predicted cleavable signal peptide, and suggests structural homology with the OmpW family of porins. Circular dichroism of recombinant TP0126 supports structural homology to OmpW. Together with the evidence that TP0126 is fully conserved among T. pallidum subspecies and strains, these data suggest an important role for TP0126 in T. pallidum biology and syphilis pathogenesis.

Bistability in a metabolic network underpins the de novo evolution of colony switching in Pseudomonas fluorescens

Phenotype switching is commonly observed in nature. This prevalence has allowed the elucidation of a number of underlying molecular mechanisms. However, little is known about how phenotypic switches arise and function in their early evolutionary stages. The first opportunity to provide empirical insight was delivered by an experiment in which populations of the bacterium Pseudomonas fluorescens SBW25 evolved, de novo, the ability to switch between two colony phenotypes. Here we unravel the molecular mechanism behind colony switching, revealing how a single nucleotide change in a gene enmeshed in central metabolism (carB) generates such a striking phenotype. We show that colony switching is underpinned by ON/OFF expression of capsules consisting of a colanic acid-like polymer. We use molecular genetics, biochemical analyses, and experimental evolution to establish that capsule switching results from perturbation of the pyrimidine biosynthetic pathway. Of central importance is a bifurcation point at which uracil triphosphate is partitioned towards either nucleotide metabolism or polymer production. This bifurcation marks a cell-fate decision point whereby cells with relatively high pyrimidine levels favour nucleotide metabolism (capsule OFF), while cells with lower pyrimidine levels divert resources towards polymer biosynthesis (capsule ON). This decision point is present and functional in the wild-type strain. Finally, we present a simple mathematical model demonstrating that the molecular components of the decision point are capable of producing switching. Despite its simple mutational cause, the connection between genotype and phenotype is complex and multidimensional, offering a rare glimpse of how noise in regulatory networks can provide opportunity for evolution.

The molecular basis of an experimentally evolved colony-level phenotype switch is revealed to affect a metabolic bifurcation point where nucleotides can be used for either cell division or capsule polymer synthesis.

Phenotype switching—the ability to switch rapidly between phenotypic states—is an evolutionary survival strategy commonly used by organisms in the face of unpredictable environmental conditions. However, little is known about how phenotype switches emerge and function in their early evolutionary stages. A previous study observed the evolutionary emergence of colony morphology switching in Pseudomonas fluorescens populations in response to fluctuating selection. Here we describe the underlying molecular basis of this colony switching, providing the first account of the mechanism behind a real-time evolved phenotype switch. We show that colony switching in this instance is underpinned at the cellular level by high frequency ON/OFF expression of colanic acid-like capsules in response to varying levels of a metabolite. Biochemical assays revealed that capsule switching results from mutations that reduce concentrations of intermediates in a central metabolic pathway—the pyrimidine biosynthetic pathway. Of key importance is the partitioning of these metabolic resources between polymer production (leading to capsulation) and cell division (leading to noncapsulation); this bifurcation marks a decision point whereby cells with low metabolite levels divert resources towards polymer production, increasing the likelihood of switching to the capsulated state. As a greater proportion of cells become capsulated, colony switching emerges. These findings show that, while colony switching evolved with relative ease, the underlying molecular mechanism is surprisingly complex.

Deep comparative genomics among Chlamydia trachomatis lymphogranuloma venereum isolates highlights genes potentially involved in pathoadaptation

Lymphogranuloma venereum (LGV) is a human sexually transmitted disease caused by the obligate intracellular bacterium Chlamydia trachomatis (serovars L1-L3). LGV clinical manifestations range from severe ulcerative proctitis (anorectal syndrome), primarily caused by the epidemic L2b strains, to painful inguinal lymphadenopathy (the typical LGV bubonic form). Besides potential host-related factors, the differential disease severity and tissue tropism among LGV strains is likely a function of the genetic backbone of the strains. We aimed to characterize the genetic variability among LGV strains as strain- or serovar-specific mutations may underlie phenotypic signatures, and to investigate the mutational events that occurred throughout the pathoadaptation of the epidemic L2b lineage. By analyzing 20 previously published genomes from L1, L2, L2b and L3 strains and two new genomes from L2b strains, we detected 1497 variant sites and about 100 indels, affecting 453 genes and 144 intergenic regions, with 34 genes displaying a clear overrepresentation of nonsynonymous mutations. Effectors and/or type III secretion substrates (almost all of those described in the literature) and inclusion membrane proteins showed amino acid changes that were about fivefold more frequent than silent changes. More than 120 variant sites occurred in plasmid-regulated virulence genes, and 66% yielded amino acid changes. The identified serovar-specific variant sites revealed that the L2b-specific mutations are likely associated with higher fitness and pointed out potential targets for future highly discriminatory diagnostic/typing tests. By evaluating the evolutionary pathway beyond the L2b clonal radiation, we observed that 90.2% of the intra-L2b variant sites occurring in coding regions involve nonsynonymous mutations, where CT456/tarp has been the main target. Considering the progress on C. trachomatis genetic manipulation, this study may constitute an important contribution for prioritizing study targets for functional genomics aiming to dissect the impact of the identified intra-LGV polymorphisms on virulence or tropism dissimilarities among LGV strains.

Rapid increase in pertactin-deficient Bordetella pertussis isolates, Australia

Acellular vaccines against Bordetella pertussis were introduced in Australia in 1997. By 2000, these vaccines had replaced whole-cell vaccines. During 2008–2012, a large outbreak of pertussis occurred. During this period, 30% (96/320) of B. pertussis isolates did not express the vaccine antigen pertactin (prn). Multiple mechanisms of prn inactivation were documented, including IS481 and IS1002 disruptions, a variation within a homopolymeric tract, and deletion of the prn gene. The mechanism of lack of expression of prn in 16 (17%) isolates could not be determined at the sequence level. These findings suggest that B. pertussis not expressing prn arose independently multiple times since 2008, rather than by expansion of a single prn-negative clone. All but 1 isolate had ptxA1, prn2, and ptxP3, the alleles representative of currently circulating strains in Australia. This pattern is consistent with continuing evolution of B. pertussis in response to vaccine selection pressure.

Automatic validation of computational models using pseudo-3D spatio-temporal model checking

BACKGROUND

Computational models play an increasingly important role in systems biology for generating predictions and in synthetic biology as executable prototypes/designs. For real life (clinical) applications there is a need to scale up and build more complex spatio-temporal multiscale models; these could enable investigating how changes at small scales reflect at large scales and viceversa. Results generated by computational models can be applied to real life applications only if the models have been validated first. Traditional in silico model checking techniques only capture how non-dimensional properties (e.g. concentrations) evolve over time and are suitable for small scale systems (e.g. metabolic pathways). The validation of larger scale systems (e.g. multicellular populations) additionally requires capturing how spatial patterns and their properties change over time, which are not considered by traditional non-spatial approaches.

RESULTS

We developed and implemented a methodology for the automatic validation of computational models with respect to both their spatial and temporal properties. Stochastic biological systems are represented by abstract models which assume a linear structure of time and a pseudo-3D representation of space (2D space plus a density measure). Time series data generated by such models is provided as input to parameterised image processing modules which automatically detect and analyse spatial patterns (e.g. cell) and clusters of such patterns (e.g. cellular population). For capturing how spatial and numeric properties change over time the Probabilistic Bounded Linear Spatial Temporal Logic is introduced. Given a collection of time series data and a formal spatio-temporal specification the model checker Mudi ( http://mudi.modelchecking.org ) determines probabilistically if the formal specification holds for the computational model or not. Mudi is an approximate probabilistic model checking platform which enables users to choose between frequentist and Bayesian, estimate and statistical hypothesis testing based validation approaches. We illustrate the expressivity and efficiency of our approach based on two biological case studies namely phase variation patterning in bacterial colony growth and the chemotactic aggregation of cells.

CONCLUSIONS

The formal methodology implemented in Mudi enables the validation of computational models against spatio-temporal logic properties and is a precursor to the development and validation of more complex multidimensional and multiscale models.

Phase-variable restriction/modification systems are required for Helicobacter pylori colonization

Background

One mechanism utilized by bacterial pathogens for host adaptation and immune evasion is the generation of phenotypic diversity by the phasevarion that results from the differential expression of a suite of genes regulated by the activity of a phase-variable methyltransferase within a restriction modification (RM) system. Phasevarions are active in Helicobacter pylori, however there have been no studies investigating the significance of phase-variable RM systems on host colonization.

Methods

Two mutant types incapable of phase variation were constructed; a clean deletion mutant (‘DEL’) and a mutant (‘ON’) where the homopolymeric repeat was replaced with a non-repeat synonymous sequence, resulting in expression of the full-length protein. The resulting mutants were assessed for their colonisation ability in the mouse model.

Results

Five phase-variable genes encoding either methyltransferases or members of RM systems were found in H. pylori OND79. Our mutants fell into three categories; 1, those with little effect on colonization, 2, those where expression of the full-length protein was detrimental, 3, those where both mutations were detrimental.

Conclusions

Our results demonstrated that phase-variable methyltransferases are critical to H. pylori colonization, suggesting that genome methylation and generation of epigenetic diversity is important for colonization and pathogenesis. The third category of mutants suggests that differential genome methylation status of H. pylori cell populations, achieved by the phasevarion, is essential for host adaptation. Studies of phase-variable RM mutants falling in the two other categories, not strictly required for colonization, represent a future perspective to investigate the role of phasevarion in persistence of H. pylori.

Phenotypic variation in the plant pathogenic bacterium Acidovorax citrulli

Acidovorax citrulli causes bacterial fruit blotch (BFB) of cucurbits, a disease that threatens the cucurbit industry worldwide. Despite the economic importance of BFB, little is known about pathogenicity and fitness strategies of the bacterium. We have observed the phenomenon of phenotypic variation in A. citrulli. Here we report the characterization of phenotypic variants (PVs) of two strains, M6 and 7a1, isolated from melon and watermelon, respectively. Phenotypic variation was observed following growth in rich medium, as well as upon isolation of bacteria from inoculated plants or exposure to several stresses, including heat, salt and acidic conditions. When grown on nutrient agar, all PV colonies possessed a translucent appearance, in contrast to parental strain colonies that were opaque. After 72 h, PV colonies were bigger than parental colonies, and had a fuzzy appearance relative to parental strain colonies that are relatively smooth. A. citrulli colonies are generally surrounded by haloes detectable by the naked eye. These haloes are formed by type IV pilus (T4P)-mediated twitching motility that occurs at the edge of the colony. No twitching haloes could be detected around colonies of both M6 and 7a1 PVs, and microscopy observations confirmed that indeed the PVs did not perform twitching motility. In agreement with these results, transmission electron microscopy revealed that M6 and 7a1 PVs do not produce T4P under tested conditions. PVs also differed from their parental strain in swimming motility and biofilm formation, and interestingly, all assessed variants were less virulent than their corresponding parental strains in seed transmission assays. Slight alterations could be detected in some DNA fingerprinting profiles of 7a1 variants relative to the parental strain, while no differences at all could be seen among M6 variants and parental strain, suggesting that, at least in the latter, phenotypic variation is mediated by slight genetic and/or epigenetic alterations.

Culture history and population heterogeneity as determinants of bacterial adaptation: the adaptomics of a single environmental transition

Diversity in adaptive responses is common within species and populations, especially when the heterogeneity of the frequently large populations found in environments is considered. By focusing on events in a single clonal population undergoing a single transition, we discuss how environmental cues and changes in growth rate initiate a multiplicity of adaptive pathways. Adaptation is a comprehensive process, and stochastic, regulatory, epigenetic, and mutational changes can contribute to fitness and overlap in timing and frequency. We identify culture history as a major determinant of both regulatory adaptations and microevolutionary change. Population history before a transition determines heterogeneities due to errors in translation, stochastic differences in regulation, the presence of aged, damaged, cheating, or dormant cells, and variations in intracellular metabolite or regulator concentrations. It matters whether bacteria come from dense, slow-growing, stressed, or structured states. Genotypic adaptations are history dependent due to variations in mutation supply, contingency gene changes, phase variation, lateral gene transfer, and genome amplifications. Phenotypic adaptations underpin genotypic changes in situations such as stress-induced mutagenesis or prophage induction or in biofilms to give a continuum of adaptive possibilities. Evolutionary selection additionally provides diverse adaptive outcomes in a single transition and generally does not result in single fitter types. The totality of heterogeneities in an adapting population increases the chance that at least some individuals meet immediate or future challenges. However, heterogeneity complicates the adaptomics of single transitions, and we propose that subpopulations will need to be integrated into future population biology and systems biology predictions of bacterial behavior.

Mass Production of the Beneficial Nematode Heterorhabditis bacteriophora and Its Bacterial Symbiont Photorhabdus luminescens

Entomoparasitic nematodes (EPNs) are being commercialized as a biocontrol measure for crop insect pests, as they provide advantages over common chemical insecticides. Mass production of these nematodes in liquid media has become a major challenge for commercialization. Producers are not willing to share the trade secrets of mass production and by doing so, have made culturing EPNs extremely difficult to advance existing technologies. Theoretically, mass production in liquid media is an ideal culturing method as it increases cost efficiency and nematode quantity. This paper will review current culturing methodologies and suggest basic culturing parameters for mass production. This review is focused on Heterorhabditis bacteriophora; however, this information can be useful for other nematode species.

Allelic diversity and phylogeny of homB, a novel co-virulence marker of Helicobacter pylori

Background

The homB gene is a Helicobacter pylori disease-marker candidate, strongly associated with peptic ulcer disease, while homA, its paralogue gene with 90% sequence identity, is correlated with non-ulcer dyspepsia. The HomB encoded outer membrane protein was shown to contribute to the proinflammatory properties of H. pylori and also to be involved in bacterial adherence.

This study investigated the distribution of homB and homA genes in 455 H. pylori strains from East Asian and Western countries, and carried out sequence comparison and phylogenetic analyses.

Results

Both homB and homA genes were heterogeneously distributed worldwide, with a marked difference between East Asian and Western strains.

Analysis of homB and homA sequences revealed diversity regarding the number of copies and their genomic localization, with East Asian and Western strains presenting different genotypes. Moreover, homB and homA sequence analysis suggests regulation by phase variation. It also indicates possible recombination events, leading to gene duplication or homB/homA conversion which may as well be implicated in the regulation of these genes. Phylogenetic reconstruction of homB and homA revealed clustering according to the geographic origin of strains. Allelic diversity in the middle region of the genes was observed for both homB and homA, although there was no correlation between any allele and disease. For each gene, a dominant worldwide allele was detected, suggesting that homB/homA allelic variants were independent of the geographical origin of the strain. Moreover, all alleles were demonstrated to be expressed in vivo.

Conclusion

Overall, these results suggest that homB and homA genes are good candidates to be part of the pool of H. pylori OMPs implicated in host-bacteria interface and also contributing to the generation of antigenic variability, and thus involved in H. pylori persistence.

Simple sequence repeat (SSR)-based gene diversity in Burkholderia pseudomallei and Burkholderia mallei

Pathogens Burkholderia pseudomallei (Bp) and Burkholderia mallei (Bm) contain a large number (> 12,000) of Simple Sequence Repeats (SSRs). To study the extent to which these features have contributed to the diversification of genes, we have conducted comparative studies with nineteen genomes of these bacteria. We found 210 genes with characteristic types of SSR variations. SSRs with nonamer repeat units were the most abundant, followed by hexamers and trimers. Amino acids with smaller and nonpolar R-groups are preferred to be encoded by the variant SSRs, perhaps due to their minimal impacts to protein functionality. A majority of these genes appears to code for surface or secreted proteins that may directly interact with the host factors during pathogenesis or other environmental factors. There also are others that encode diverse functions in the cytoplasm, and this protein variability may reflect an extensive involvement of phase variation in survival and adaptation of these pathogens.

Heterogeneity in tissue culture infection models: a source of novel host-pathogen interactions?

Tissue cultures have been successfully exploited to dissect cellular and molecular mechanisms of microbial infections. Most of the methods used in this model conclude with data describing host and pathogen 'average' responses. Microscopy, however, reveals that such interplay is very diverse and that both partners are composed of phenotypically heterogeneous populations. Thus, upon co-incubation in the plate assay, neither all cultured host cells are infected nor all pathogen cells inflict alterations in host physiology. Despite its obvious impact in data interpretation, the basis of this heterogeneity remains in most cases unknown. Addressing this issue is encouraging since may contribute to uncover novel interactions in the host-pathogen scenario.

Further characterization of Vibrio vulnificus rugose variants and identification of a capsular and rugose exopolysaccharide gene cluster

Capsular polysaccharide (CPS) is a major virulence factor in Vibrio vulnificus, and encapsulated strains have an opaque, smooth (OpS) colony morphology, while nonencapsulated strains have a translucent, smooth (TrS) colony morphology. Previously, we showed that OpS and TrS parental strains can yield a third colony type, rugose (R), and that the resulting strains, with the OpR and TrR phenotypes, respectively, form copious biofilms. Here we show that while OpR and TrR strains both produce three-dimensional biofilm structures that are indicative of rugose extracellular polysaccharide (rEPS) production, OpR strains also retain expression of CPS and are virulent in an iron-supplemented mouse model, while TrR strains lack CPS and are avirulent. Chlorine resistance assays further distinguished OpR and TrR isolates as exposure to 3 microg/ml NaOCl eradicated both OpS and OpR strains, while both TrS and TrR strains survived, but at rates which were significantly different from one another. Taken together, these results further emphasize the importance of CPS for virulence of V. vulnificus and establish a correlation between CPS expression and chlorine sensitivity in this organism. Using reverse transcriptase PCR, we also identified a nine-gene cluster associated with both CPS and rEPS expression in V. vulnificus, designated the wcr (capsular and rugose polysaccharide) locus, with expression occurring primarily in R variants. The latter results set the stage for characterization of functional determinants which individually or collectively contribute to expression of multiple EPS forms in this pathogen.

Functional analysis of the M.HpyAIV DNA methyltransferase of Helicobacter pylori

A large number of genes encoding restriction-modification (R-M) systems are found in the genome of the human pathogen Helicobacter pylori. R-M genes comprise approximately 10% of the strain-specific genes, but the relevance of having such an abundance of these genes is not clear. The type II methyltransferase (MTase) M.HpyAIV, which recognizes GANTC sites, was present in 60% of the H. pylori strains analyzed, whereof 69% were resistant to restriction enzyme digestion, which indicated the presence of an active MTase. H. pylori strains with an inactive M.HpyAIV phenotype contained deletions in regions of homopolymers within the gene, which resulted in premature translational stops, suggesting that M.HpyAIV may be subjected to phase variation by a slipped-strand mechanism. An M.HpyAIV gene mutant was constructed by insertional mutagenesis, and this mutant showed the same viability and ability to induce interleukin-8 in epithelial cells as the wild type in vitro but had, as expected, lost the ability to protect its self-DNA from digestion by a cognate restriction enzyme. The M.HpyAIV from H. pylori strain 26695 was overexpressed in Escherichia coli, and the protein was purified and was able to bind to DNA and protect GANTC sites from digestion in vitro. A bioinformatic analysis of the number of GANTC sites located in predicted regulatory regions of H. pylori strains 26695 and J99 resulted in a number of candidate genes. katA, a selected candidate gene, was further analyzed by quantitative real-time reverse transcription-PCR and shown to be significantly down-regulated in the M.HpyAIV gene mutant compared to the wild-type strain. This demonstrates the influence of M.HpyAIV methylation in gene expression.

Phase variation and microevolution at homopolymeric tracts in Bordetella pertussis

BACKGROUND

Bordetella pertussis, the causative agent of whooping cough, is a highly clonal pathogen of the respiratory tract. Its lack of genetic diversity, relative to many bacterial pathogens, could limit its ability to adapt to a hostile and changing host environment. This limitation might be overcome by phase variation, as observed for other mucosal pathogens. One of the most common mechanisms of phase variation is reversible expansion or contraction of homopolymeric tracts (HPTs).

RESULTS

The genomes of B. pertussis and the two closely related species, B. bronchiseptica and B. parapertussis, were screened for homopolymeric tracts longer than expected on the basis of chance, given their nucleotide compositions. Sixty-nine such HPTs were found in total among the three genomes, 74% of which were polymorphic among the three species. Nine HPTs were genotyped in a collection of 90 geographically and temporally diverse B. pertussis strains using the polymerase chain reaction/ligase detection reaction (PCR/LDR) assay. Six HPTs were polymorphic in this collection of B. pertussis strains. Of note, one of these polymorphic HPTs was found in the fimX promoter, where a single base insertion variant was present in seven strains, all of which were isolated prior to introduction of the pertussis vaccine. Transcript abundance of fimX was found to be 3.8-fold lower in strains carrying the longer allele. HPTs in three other genes, tcfA, bapC, and BP3651, varied widely in composition across the strain collection and displayed allelic polymorphism within single cultures.

CONCLUSION

Allelic polymorphism at homopolymeric tracts is common within the B. pertussis genome. Phase variability may be an important mechanism in B. pertussis for evasion of the immune system and adaptation to different niches in the human host. High sensitivity and specificity make the PCR/LDR assay a powerful tool for investigating allelic variation at HPTs. Using this method, allelic diversity and phase variation were demonstrated at several B. pertussis loci.

Population-associated differences between the phase variable LPS biosynthetic genes of Helicobacter pylori

Background

Population structures are normally determined using genes under minimal functional selection. In this study we have assessed genes that are not always essential, show differences in alleles between strains, and are involved in the directly host-selectable phenotype of LPS biosynthesis.

Results

Eight complete LPS biosynthesis genes, seven of which are associated with phase variation in some or all strains of Helicobacter pylori, have been sequenced and their divergence analyzed. The differences observed indicate that recombination within these genes largely reflects exchange between strains within the population lineages previously determined on the basis of MLST using housekeeping genes. This indicates that the differences that are used for MLST are likely to broadly associate with genes under functional selection, and differences in strain behaviour. However, instances of exchange between the subpopulations were identified, including the hpAfrica2 subpopulation. Further, there were other differences in gene complements and the chromosomal location of genes indicative of greater diversity within the population than is revealed by the available genome sequences and comparative genome hybridization studies.

Conclusion

These results indicate that the described population structure based upon MLST is broadly a good basis for studying the biology of H. pylori, but that individual alleles may not follow these associations. As a consequence, when working in unsequenced strains, it is necessary to carefully check the presence, sequence, and distribution of any individual gene of interest.

Phase variation and genomic architecture changes in Azospirillum

The plant growth-promoting rhizobacterium Azospirillum lipoferum 4B generates in vitro at high frequency a stable nonswimming phase variant designated 4V(I), which is distinguishable from the wild type by the differential absorption of dyes. The frequency of variants generated by a recA mutant of A. lipoferum 4B was increased up to 10-fold. The pleiotropic modifications characteristic of the phase variant are well documented, but the molecular processes involved are unknown. Here, the objective was to assess whether genomic rearrangements take place during phase variation of strain 4B. The random amplified polymorphic DNA (RAPD) profiles of strains 4B and 4V(I) differed. RAPD fragments observed only with the wild type were cloned, and three cosmids carrying the corresponding fragments were isolated. The three cosmids hybridized with a 750-kb plasmid and pulse-field gel electrophoresis analysis revealed that this replicon was missing in the 4V(I) genome. The same rearrangements took place during phase variation of 4BrecA. Large-scale genomic rearrangements during phase variation were demonstrated for two additional strains. In Azospirillum brasilense WN1, generation of stable variants was correlated with the disappearance of a replicon of 260 kb. For Azospirillum irakense KBC1, the variant was not stable and coincided with the formation of a new replicon, whereas the revertant recovered the parental genomic architecture. This study shows large-scale genomic rearrangements in Azospirillum strains and correlates them with phase variation.

Comparative overview of the genomic and genetic differences between the pathogenic Neisseria strains and species

The availability of complete genome sequences from multiple pathogenic Neisseria strains and species has enabled a comprehensive survey of the genomic and genetic differences occurring within these species. In this review, we describe the chromosomal rearrangements that have occurred, and the genomic islands and prophages that have been identified in the various genomes. We also describe instances where specific genes are present or absent, other instances where specific genes have been inactivated, and situations where there is variation in the version of a gene that is present. We also provide an overview of mosaic genes present in these genomes, and describe the variation systems that allow the expression of particular genes to be switched ON or OFF. We have also described the presence and location of mobile non-coding elements in the various genomes. Finally, we have reviewed the incidence and properties of various extra-chromosomal elements found within these species. The overall impression is one of genomic variability and instability, resulting in increased functional flexibility within these species.

Insights in the Pathogenic Potential of Propionibacterium acnes From Its Complete Genome

One of the long-lasting open questions in understanding acne is the role of the skin bacterium Propionibacterium acnes. Opinions within the scientific community diverge regarding the importance of this Gram-positive bacterium not only in acne but also in other P. acnes-associated diseases. The recently decoded genome gives us some clues with respect to its pathogenic potential and its strategies to survive in the harsh environment of human skin. This review focuses on the major findings of the genome sequence. Special attention has been drawn to the cell surface composition of P. acnes, its antigenic potential, and to enzymes that are potentially involved in degrading skin-derived substances.

Strain-specific differences in Neisseria gonorrhoeae associated with the phase variable gene repertoire

Background

There are several differences associated with the behaviour of the four main experimental Neisseria gonorrhoeae strains, FA1090, FA19, MS11, and F62. Although there is data concerning the gene complements of these strains, the reasons for the behavioural differences are currently unknown. Phase variation is a mechanism that occurs commonly within the Neisseria spp. and leads to switching of genes ON and OFF. This mechanism may provide a means for strains to express different combinations of genes, and differences in the strain-specific repertoire of phase variable genes may underlie the strain differences.

Results

By genome comparison of the four publicly available neisserial genomes a revised list of 64 genes was created that have the potential to be phase variable in N. gonorrhoeae, excluding the opa and pilC genes. Amplification and sequencing of the repeat-containing regions of these genes allowed determination of the presence of the potentially unstable repeats and the ON/OFF expression state of these genes. 35 of the 64 genes show differences in the composition or length of the repeats, of which 28 are likely to be associated with phase variation. Two genes were expressed differentially between strains causing disseminated infection and uncomplicated gonorrhoea. Further study of one of these in a range of clinical isolates showed this association to be due to sample size and is not maintained in a larger sample.

Conclusion

The results provide us with more evidence as to which genes identified through comparative genomics are indeed phase variable. The study indicates that there are large differences between these four N. gonorrhoeae strains in terms of gene expression during in vitro growth. It does not, however, identify any clear patterns by which previously reported behavioural differences can be correlated with the phase variable gene repertoire.

Phase variation mediated niche adaptation during prolonged experimental murine infection with Helicobacter pylori

Changes in the repeats associated with the recently redefined repertoire of 31 phase-variable genes in Helicobacter pylori were investigated following murine gastric colonization for up to one year in three unrelated H. pylori strains. Between the beginning and end of the experimental period, changes were seen in ten genes (32 %), which would alter gene expression in one or more of the three strains studied. For those genes that showed repeat length changes at the longest time points, intermediate time points showed differences between the rates of change for different functional groups of genes. Genes most likely to be associated with immediate niche fitting changed most rapidly, including phospholipase A (pldA) and LPS biosynthetic genes. Other surface proteins, which may be under adaptive immune selection, changed more slowly. Restriction-modification genes showed no particular temporal pattern. The number of genes that phase varied during adaptation to the murine gastric environment correlated inversely with their relative fitness as previously determined in this murine model of colonization. This suggests a role for these genes in determining initial fitness for colonization as well as in subsequent niche adaptation. In addition, a coding tandem repeat within a phase-variable gene which does not control actual gene expression was also investigated. This repeat was found to vary in copy number during colonization. This suggests that changes in the structures encoded by tandem repeats may also play a role in altered protein functions and/or immune evasion during H. pylori colonization.

The diversity within an expanded and redefined repertoire of phase-variable genes in Helicobacter pylori

Phase variation is a common mechanism used by pathogenic bacteria to generate intra-strain diversity that is important in niche adaptation and is strongly associated with virulence determinants. Previous analyses of the complete sequences of the Helicobacter pylori strains 26695 and J99 have identified 36 putative phase-variable genes among the two genomes through their association with homopolymeric tracts and dinucleotide repeats. Here a comparative analysis of the two genomes is reported and an updated and expanded list of 46 candidate phase-variable genes in H. pylori is described. These have been systematically investigated by PCR and sequencing for the presence of the genes, and the presence and variability in length of the repeats in strains 26695 and J99 and in a collection of unrelated H. pylori strains representative of the main global subdivisions recently suggested. This provides supportive evidence for the phase variability of 30 of the 46 candidates. Other differences in this subset of genes were observed (i) in the repeats, which can be present or absent among the strains, or stabilized in different strains and (ii) in the gene-complements of the strains. Differences between genes were not consistently correlated with the geographic population distribution of the strains. This study extends and provides new evidence for variation of this type in H. pylori, and of the high degree of diversity of the repertoire of genes which display phase-variable switching within individual strains.

Mutation rates: estimating phase variation rates when fitness differences are present and their impact on population structure

Phase variation is a mechanism of ON-OFF switching that is widely utilized by bacterial pathogens. There is currently no standardization to how the rate of phase variation is determined experimentally, and traditional methods of mutation rate estimation may not be appropriate to this process. Here, the history of mutation rate estimation is reviewed, describing the existing methods available. A new mathematical model that can be applied to this problem is also presented. This model specifically includes the confounding factors of back-mutation and the influence of fitness differences between the alternate phenotypes. These are central features of phase variation but are rarely addressed, with the result that some previously estimated phase variation rates may have been significantly overestimated. It is shown that, conversely, the model can also be used to investigate fitness differences if mutation rates are approximately known. In addition, stochastic simulations of the model are used to explore the impact of 'jackpot cultures' on the mutation rate estimation. Using the model, the impact of realistic rates and selection on population structure is investigated. In the absence of fitness differences it is predicted that there will be phenotypic stability over many generations. The rate of phenotypic change within a population is likely, therefore, to be principally determined by selection. A greater insight into the population dynamics of mutation rate processes can be gained if populations are monitored over successive time points.