Microarray-based genomic surveying of gene polymorphisms in Chlamydia trachomatis

Correspondence of aCGH and long-read genome assembly for detection of copy number differences: A proof-of-concept with cichlid genomes

Copy number variation is an important source of genetic variation, yet data are often lacking due to technical limitations for detection given the current genome assemblies. Our goal is to demonstrate the extent to which an array-based platform (aCGH) can identify genomic loci that are collapsed in genome assemblies that were built with short-read technology. Taking advantage of two cichlid species for which genome assemblies based on Illumina and PacBio are available, we show that inter-species aCGH log2 hybridization ratios correlate more strongly with inferred copy number differences based on PacBio-built genome assemblies than based on Illumina-built genome assemblies. With regard to inter-species copy number differences of specific genes identified by each platform, the set identified by aCGH intersects to a greater extent with the set identified by PacBio than with the set identified by Illumina. Gene function, according to Gene Ontology analysis, did not substantially differ among platforms, and platforms converged on functions associated with adaptive phenotypes. The results of the current study further demonstrate that aCGH is an effective platform for identifying copy number variable sequences, particularly those collapsed in short read genome assemblies.

The Genomic Substrate for Adaptive Radiation: Copy Number Variation across 12 Tribes of African Cichlid Species

The initial sequencing of five cichlid genomes revealed an accumulation of genetic variation, including extensive copy number variation in cichlid lineages particularly those that have undergone dramatic evolutionary radiation. Gene duplication has the potential to generate substantial molecular substrate for the origin of evolutionary novelty. We use array-based comparative heterologous genomic hybridization to identify copy number variation events (CNVEs) for 168 samples representing 53 cichlid species including the 5 species for which full genome sequence is available. We identify an average of 50-100 CNVEs per individual. For those species represented by multiple samples, we identify 150-200 total CNVEs suggesting a substantial amount of intraspecific variation. For these species, only ∼10% of the detected CNVEs are fixed. Hierarchical clustering of species according to CNVE data recapitulates phylogenetic relationships fairly well at both the tribe and radiation level. Although CNVEs are detected on all linkage groups, they tend to cluster in "hotspots" and are likely to contain and be flanked by transposable elements. Furthermore, we show that CNVEs impact functional categories of genes with potential roles in adaptive phenotypes that could reasonably promote divergence and speciation in the cichlid clade. These data contribute to a more complete understanding of the molecular basis for adaptive natural selection, speciation, and evolutionary radiation.

Genetic Transformation of a C. trachomatis Ocular Isolate With the Functional Tryptophan Synthase Operon Confers an Indole-Rescuable Phenotype

Chlamydia trachomatis is the leading cause of preventable blindness and the most common bacterial sexually transmitted infection. Different strains are associated with ocular or urogenital infections, and a proposed mechanism that may explain this tissue tropism is the active tryptophan biosynthesis pathway encoded by the genomic trpRBA operon in urogenital strains. Here we describe genetic complementation studies that are essential to confirm the role of tryptophan synthase in the context of an ocular C. trachomatis genomic background. Ocular strain A2497 was transformed with the (urogenital) pSW2::GFP shuttle vector showing that there is no strain tropism barrier to this plasmid vector; moreover, transformation had no detrimental effect on the growth kinetics of A2497, which is important given the low transformation efficiency of C. trachomatis. A derivative of the pSW2::GFP vector was used to deliver the active tryptophan biosynthesis genes from a urogenital strain of C. trachomatis (Soton D1) to A2497 with the aim of complementing the truncated trpA gene common to most ocular strains. After confirmation of intact TrpA protein expression in the transformed A2497, the resulting transformants were cultivated in tryptophan-depleted medium with and without indole or tryptophan, showing that complementation of the truncated trpA gene by the intact and functional urogenital trpRBA operon was sufficient to bestow an indole rescuable phenotype upon A2497. This study proves that pSW2::GFP derived vectors do not conform to the cross-strain transformation barrier reported for other chlamydia shuttle vectors, suggesting these as a universal vector for transformation of all C. trachomatis strains. This vector promiscuity enabled us to test the indole rescue hypothesis by transforming ocular strain A2497 with the functional urogenital trpRBA operon, which complemented the non-functional tryptophan synthase. These data confirm that the trpRBA operon is necessary and sufficient for chlamydia to survive in tryptophan-limited environments such as the female urogenital tract.

Population-based analysis of ocular Chlamydia trachomatis in trachoma-endemic West African communities identifies genomic markers of disease severity

BACKGROUND

Chlamydia trachomatis (Ct) is the most common infectious cause of blindness and bacterial sexually transmitted infection worldwide. Ct strain-specific differences in clinical trachoma suggest that genetic polymorphisms in Ct may contribute to the observed variability in severity of clinical disease.

METHODS

Using Ct whole genome sequences obtained directly from conjunctival swabs, we studied Ct genomic diversity and associations between Ct genetic polymorphisms with ocular localization and disease severity in a treatment-naïve trachoma-endemic population in Guinea-Bissau, West Africa.

RESULTS

All Ct sequences fall within the T2 ocular clade phylogenetically. This is consistent with the presence of the characteristic deletion in trpA resulting in a truncated non-functional protein and the ocular tyrosine repeat regions present in tarP associated with ocular tissue localization. We have identified 21 Ct non-synonymous single nucleotide polymorphisms (SNPs) associated with ocular localization, including SNPs within pmpD (odds ratio, OR = 4.07, p* = 0.001) and tarP (OR = 0.34, p* = 0.009). Eight synonymous SNPs associated with disease severity were found in yjfH (rlmB) (OR = 0.13, p* = 0.037), CTA0273 (OR = 0.12, p* = 0.027), trmD (OR = 0.12, p* = 0.032), CTA0744 (OR = 0.12, p* = 0.041), glgA (OR = 0.10, p* = 0.026), alaS (OR = 0.10, p* = 0.032), pmpE (OR = 0.08, p* = 0.001) and the intergenic region CTA0744-CTA0745 (OR = 0.13, p* = 0.043).

CONCLUSIONS

This study demonstrates the extent of genomic diversity within a naturally circulating population of ocular Ct and is the first to describe novel genomic associations with disease severity. These findings direct investigation of host-pathogen interactions that may be important in ocular Ct pathogenesis and disease transmission.

Differential infection outcome of Chlamydia trachomatis in human blood monocytes and monocyte-derived dendritic cells

Background

Chlamydia trachomatis is an intracellular bacteria which consist of three biovariants; trachoma (serovars A-C), urogenital (serovars D-K) and lymphogranuloma venereum (L1-L3), causing a wide spectrum of disease in humans. Monocytes are considered to disseminate this pathogen throughout the body while dendritic cells (DCs) play an important role in mediating immune response against bacterial infection. To determine the fate of C. trachomatis within human peripheral blood monocytes and monocyte-derived DCs, these two sets of immune cells were infected with serovars Ba, D and L2, representative of the three biovariants of C. trachomatis.

Results

Our study revealed that the different serovars primarily infect monocytes and DCs in a comparable fashion, however undergo differential infection outcome, serovar L2 being the only candidate to inflict active infection. Moreover, the C. trachomatis serovars Ba and D become persistent in monocytes while the serovars predominantly suffer degradation within DCs. Effects of persistence gene Indoleamine 2, 3-dioxygenase (IDO) was not clearly evident in the differential infection outcome. The heightened levels of inflammatory cytokines secreted by the chlamydial infection in DCs compared to monocytes seemed to be instrumental for this consequence. The immune genes induced in monocytes and DCs against chlamydial infection involves a different set of Toll-like receptors, indicating that distinct intracellular signalling pathways are adopted for immune response.

Conclusion

Our results demonstrate that the host pathogen interaction in chlamydia infection is not only serovar specific but manifests cell specific features, inducing separate immune response cascade in monocytes and DCs.

Role for chlamydial inclusion membrane proteins in inclusion membrane structure and biogenesis

The chlamydial inclusion membrane is extensively modified by the insertion of type III secreted effector proteins. These inclusion membrane proteins (Incs) are exposed to the cytosol and share a common structural feature of a long, bi-lobed hydrophobic domain but little or no primary amino acid sequence similarity. Based upon secondary structural predictions, over 50 putative inclusion membrane proteins have been identified in Chlamydia trachomatis. Only a limited number of biological functions have been defined and these are not shared between chlamydial species. Here we have ectopically expressed several C. trachomatis Incs in HeLa cells and find that they induce the formation of morphologically distinct membranous vesicular compartments. Formation of these vesicles requires the bi-lobed hydrophobic domain as a minimum. No markers for various cellular organelles were observed in association with these vesicles. Lipid probes were incorporated by the Inc-induced vesicles although the lipids incorporated were dependent upon the specific Inc expressed. Co-expression of Inc pairs indicated that some colocalized in the same vesicle, others partially overlapped, and others did not associate at all. Overall, it appears that Incs may have an intrinsic ability to induce membrane formation and that individual Incs can induce membranous structures with unique properties.

Polymorphisms in inc proteins and differential expression of inc genes among Chlamydia trachomatis strains correlate with invasiveness and tropism of lymphogranuloma venereum isolates

Chlamydia trachomatis is a human bacterial pathogen that multiplies only within an intracellular membrane-bound vacuole, the inclusion. C. trachomatis includes ocular and urogenital strains, usually causing infections restricted to epithelial cells of the conjunctiva and genital mucosa, respectively, and lymphogranuloma venereum (LGV) strains, which can infect macrophages and spread into lymph nodes. However, C. trachomatis genomes display >98% identity at the DNA level. In this work, we studied whether C. trachomatis Inc proteins, which have a bilobed hydrophobic domain that may mediate their insertion in the inclusion membrane, could be a factor determining these different types of infection and tropisms. Analyses of polymorphisms and phylogeny of 48 Inc proteins from 51 strains encompassing the three disease groups showed significant amino acid differences that were mainly due to variations between Inc proteins from LGV and ocular or urogenital isolates. Studies of the evolutionary dynamics of inc genes suggested that 10 of them are likely under positive selection and indicated that most nonsilent mutations are LGV specific. Additionally, real-time quantitative PCR analyses in prototype and clinical strains covering the three disease groups identified three inc genes with LGV-specific expression. We determined the transcriptional start sites of these genes and found LGV-specific nucleotides within their promoters. Thus, subtle variations in the amino acids of a subset of Inc proteins and in the expression of inc genes may contribute to the unique tropism and invasiveness of C. trachomatis LGV strains.

Evolutionary insights into scleractinian corals using comparative genomic hybridizations

Background

Coral reefs belong to the most ecologically and economically important ecosystems on our planet. Yet, they are under steady decline worldwide due to rising sea surface temperatures, disease, and pollution. Understanding the molecular impact of these stressors on different coral species is imperative in order to predict how coral populations will respond to this continued disturbance. The use of molecular tools such as microarrays has provided deep insight into the molecular stress response of corals. Here, we have performed comparative genomic hybridizations (CGH) with different coral species to an Acropora palmata microarray platform containing 13,546 cDNA clones in order to identify potentially rapidly evolving genes and to determine the suitability of existing microarray platforms for use in gene expression studies (via heterologous hybridization).

Results

Our results showed that the current microarray platform for A. palmata is able to provide biological relevant information for a wide variety of coral species covering both the complex clade as well the robust clade. Analysis of the fraction of highly diverged genes showed a significantly higher amount of genes without annotation corroborating previous findings that point towards a higher rate of divergence for taxonomically restricted genes. Among the genes with annotation, we found many mitochondrial genes to be highly diverged in M. faveolata when compared to A. palmata, while the majority of nuclear encoded genes maintained an average divergence rate.

Conclusions

The use of present microarray platforms for transcriptional analyses in different coral species will greatly enhance the understanding of the molecular basis of stress and health and highlight evolutionary differences between scleractinian coral species. On a genomic basis, we show that cDNA arrays can be used to identify patterns of divergence. Mitochondrion-encoded genes seem to have diverged faster than nuclear encoded genes in robust corals. Accordingly, this needs to be taken into account when using mitochondrial markers for scleractinian phylogenies.

Evolution and conservation of predicted inclusion membrane proteins in chlamydiae

Chlamydia spp. are obligate intracellular pathogens that replicate within a vacuole termed the inclusion. Chlamydiae extensively modify the inclusion membrane via the insertion of chlamydial inclusion membrane proteins (Incs) which decorate the cytosolic face of the inclusion. We have assessed the overall relatedness and phylogeny of Incs in order to identify potential evolutionary trends. Despite a high degree of conservation among Incs within C. trachomatis serovars, phylogenetic analysis showed that some Incs cluster according to clinical groupings suggesting that certain Incs may contribute to tissue tropism. Bioinformatic predictions identified Incs in five chlamydial species: 55 in C. trachomatis, 68 in C. felis, 92 in C. pneumoniae, 79 in C. caviae, and 54 in C. muridarum. Inc homologues were compared between chlamydial species and 23 core Incs were identified as shared among all species. Genomic expansion of Incs was identified in C. pneumoniae, C. caviae, and C. felis but not C. trachomatis or C. muridarum.

Nucleotide and phylogenetic analyses of the Chlamydia trachomatis ompA gene indicates it is a hotspot for mutation

Background

Serovars of the human pathogen Chlamydia trachomatis occupy one of three specific tissue niches. Genomic analyses indicate that the serovars have a phylogeny congruent with their pathobiology and have an average substitution rate of less than one nucleotide per kilobase. In contrast, the gene that determines serovar specificity, ompA, has a phylogenetic association that is not congruent with tissue tropism and has a degree of nucleotide variability much higher than other genomic loci. The ompA gene encodes the major surface-exposed antigenic determinant, and the observed nucleotide diversity at the ompA locus is thought to be due to recombination and host immune selection pressure. The possible contribution of a localized increase in mutation rate, however, has not been investigated.

Results

Nucleotide diversity and phylogenetic relationships of the five constant and four variable domains of the ompA gene, as well as several loci surrounding ompA, were examined for each serovar. The loci flanking the ompA gene demonstrated that nucleotide diversity increased monotonically as ompA is approached and that their gene trees are not congruent with either ompA or tissue tropism. The variable domains of the ompA gene had a very high level of non-synonymous change, which is expected as these regions encode the surface-exposed epitopes and are under positive selection. However, the synonymous changes are clustered in the variable regions compared to the constant domains; if hitchhiking were to account for the increase in synonymous changes, these substitutions should be more evenly distributed across the gene. Recombination also cannot entirely account for this increase as the phylogenetic relationships of the constant and variable domains are congruent with each other.

Conclusions

The high number of synonymous substitutions observed within the variable domains of ompA appears to be due to an increased mutation rate within this region of the genome, whereas the increase in nucleotide substitution rate and the lack of phylogenetic congruence in the regions flanking ompA are characteristic motifs of gene conversion. Together, the increased mutation rate in the ompA gene, in conjunction with gene conversion and positive selection, results in a high degree of variability that promotes host immune evasion.

Genomic and antigenic characterization of monomeric autotransporters of Haemophilus parasuis: an ongoing process of reductive evolution

The genome of the highly pathogenic Haemophilus parasuis Nagasaki strain (serovar 5) was sequenced to 99 % completion. A genomic comparison with two other pathogenic serovar 5 H. parasuis strains identified six genes per genome (bmaA1-bmaA6) encoding β-barrel monomeric autotransporters, bmaA2 and bmaA3 being pseudogenes in at least one strain. The remaining encoded proteins were predicted to belong to the subtilisin (BmaA1 and BmaA4) and cysteine (BmaA5 and BmaA6) protease families. Allelic polymorphism was detected in other H. parasuis strains by comparative genomic hybridization using microarrays. Recombination events were observed, some of them leading to gene disruption in one of the three strains, although synteny around bmaA genes was conserved. These results suggest that bmaA genes are undergoing a process of reductive evolution. To evaluate their use as potential vaccine antigens, the products of the passenger domains of bmaA1, bmaA4, bmaA5 and bmaA6 were produced in Escherichia coli as recombinant proteins. They were detected by immunoblotting using sera of colostrum-deprived piglets recovering from a sublethal infection with H. parasuis (Nagasaki). The existence of specific antibodies after infection with H. parasuis also demonstrated in vivo expression. Using proteomics, only BmaA6 was detected in the in vitro-grown Nagasaki strain. Interestingly, the translocator domain was found in the outer membrane, while the passenger domain was located in supernatants. These results indicate that BmaA proteins could be considered as immunogen candidates to improve H. parasuis vaccines. However, their capacity to confer protective immunity needs to be studied further.

Array CGH phylogeny: how accurate are comparative genomic hybridization-based trees?

Background

Array-based Comparative Genomic Hybridization (CGH) data have been used to infer phylogenetic relationships. However, the reliability of array CGH analysis to determine evolutionary relationships has not been well established. In most CGH work, all species and strains are compared to a single reference species, whose genome was used to design the array. In the accompanying work, we critically evaluated CGH-based phylogeny using simulated competitive hybridization data. This work showed that a limited number of conditions, principally the tree topology and placement of the reference taxon in the tree, had a strong effect on the ability to recover the correct tree topology. Here, we add to our simulation study by testing the use of CGH as a phylogenetic tool with experimental CGH data from competitive hybridizations between N. crassa and other Neurospora species. In the discussion, we add to our empirical study of Neurospora by reanalyzing of data from a previous CGH phylogenetic analysis of the yeast sensu stricto complex.

Results

Array ratio data for Neurospora and related species were normalized with loess, robust spline, and linear ratio based methods, and then used to construct Neighbor-Joining and parsimony trees. These trees were compared to published phylogenetic analyses for Neurospora based on multilocus sequence analysis (MLSA). For the Neurospora dataset, the best combination of methods resulted in recovery of the MLSA tree topology less than half the time. Our reanalysis of a yeast dataset found that trees identical to established phylogeny were recovered only by pruning taxa - including the reference taxon - from the analysis.

Conclusion

Our results indicate that CGH data can be problematic for phylogenetic analysis. Success fluctuates based on the methods utilized to construct the tree and the taxa included. Selective pruning of the taxa improves the results - an impractical approach for normal phylogenetic analysis. From the more successful methods we make suggestions on the normalization and post-normalization methods that work best in estimating genetic distance between taxa.

Diverse requirements for SRC-family tyrosine kinases distinguish chlamydial species

Chlamydiae are well known for their species specificity and tissue tropism, and yet the individual species and strains show remarkable genomic synteny and share an intracellular developmental cycle unique in the microbial world. Only a relatively few chlamydial genes have been linked to specific disease or tissue tropism. Here we show that chlamydial species associated with human infections, Chlamydia trachomatis and C. pneumoniae, exhibit unique requirements for Src-family kinases throughout their developmental cycle. Utilization of Src-family kinases by C. trachomatis includes tyrosine phosphorylation of the secreted effector Tarp during the entry process, a functional role in microtubule-dependent trafficking to the microtubule organizing center, and a requirement for Src-family kinases for successful initiation of development. Nonhuman chlamydial species C. caviae and C. muridarum show none of these requirements and, instead, appear to be growth restricted by the activities of Src-family kinases. Depletion of Src-family kinases triggers a more rapid development of C. caviae with up to an 800% increase in infectious progeny production. Collectively, the results suggest that human chlamydial species have evolved requirements for tyrosine phosphorylation by Src-family kinases that are not seen in other chlamydial species. The requirement for Src-family kinases thus represents a fundamental distinction between chlamydial species that would not be readily apparent in genomic comparisons and may provide insights into chlamydial disease association and species specificity.

IMPORTANCE

Chlamydiae are well known for their species specificity and tissue tropism as well as their association with unique diseases. A paradox in the field relates to the remarkable genomic synteny shown among chlamydiae and the very few chlamydial genes linked to specific diseases. We have found that different chlamydial species exhibit unique requirements for Src-family kinases. These differing requirements for Src-family kinases would not be apparent in genomic comparisons and appear to be a previously unrecognized distinction that may provide insights to guide research in chlamydial pathogenesis.

Phylogenetic analysis of Chlamydia trachomatis Tarp and correlation with clinical phenotype

Chlamydia trachomatis is the leading cause of infectious blindness worldwide and is the most commonly reported pathogen causing sexually transmitted infections. Tarp (translocated actin recruiting phosphoprotein), a type III secreted effector that mediates actin nucleation, is central to C. trachomatis infection. The phylogenetic analysis of tarP from reference strains as well as ocular, genital, and lymphogranuloma venereum (LGV) clinical isolates demonstrated an evolutionary relationship with disease phenotype, with LGV and ocular isolates branched into clades that were separate from the urogenital isolates. The sequence analysis of Tarp indicated a high degree of variability and identified trends within clinical groupings. Tarps from LGV strains contained the highest number of tyrosine-rich repeat regions (up to nine) and the fewest (two) predicted actin binding domains. The converse was noted for Tarp proteins from ocular isolates that contained up to four actin binding domains and as few as one tyrosine-rich repeat region. The results suggest that Tarp is among the few known genes to play a role in C. trachomatis adaptations to specific niches within the host.

A critical assessment of cross-species detection of gene duplicates using comparative genomic hybridization

Background

Comparison of genomic DNA among closely related strains or species is a powerful approach for identifying variation in evolutionary processes. One potent source of genomic variation is gene duplication, which is prevalent among individuals and species. Array comparative genomic hybridization (aCGH) has been successfully utilized to detect this variation among lineages. Here, beyond the demonstration that gene duplicates among species can be quantified with aCGH, we consider the effect of sequence divergence on the ability to detect gene duplicates.

Results

Using the X chromosome genomic content difference between male D. melanogaster and female D. yakuba and D. simulans, we describe a decrease in the ability to accurately measure genomic content (copy number) for orthologs that are only 90% identical. We demonstrate that genome characteristics (e.g. chromatin environment and non-orthologous sequence similarity) can also affect the ability to accurately measure genomic content. We describe a normalization strategy and statistical criteria to be used for the identification of gene duplicates among any species group for which an array platform is available from a closely related species.

Conclusions

Array CGH can be used to effectively identify gene duplication and genome content; however, certain biases are present due to sequence divergence and other genome characteristics resulting from the divergence between lineages. Highly conserved gene duplicates will be more readily recovered by aCGH. Duplicates that have been retained for a selective advantage due to directional selection acting on many loci in one or both gene copies are likely to be under-represented. The results of this study should inform the interpretation of both previously published and future work that employs this powerful technique.

Using comparative genomic hybridization to survey genomic sequence divergence across species: a proof-of-concept from Drosophila

Background

Genome-wide analysis of sequence divergence among species offers profound insights into the evolutionary processes that shape lineages. When full-genome sequencing is not feasible for a broad comparative study, we propose the use of array-based comparative genomic hybridization (aCGH) in order to identify orthologous genes with high sequence divergence. Here we discuss experimental design, statistical power, success rate, sources of variation and potential confounding factors. We used a spotted PCR product microarray platform from Drosophila melanogaster to assess sequence divergence on a gene-by-gene basis in three fully sequenced heterologous species (D. sechellia, D. simulans, and D. yakuba). Because complete genome assemblies are available for these species this study presents a powerful test for the use of aCGH as a tool to measure sequence divergence.

Results

We found a consistent and linear relationship between hybridization ratio and sequence divergence of the sample to the platform species. At higher levels of sequence divergence (< 92% sequence identity to D. melanogaster) ~84% of features had significantly less hybridization to the array in the heterologous species than the platform species, and thus could be identified as "diverged". At lower levels of divergence (≥ 97% identity), only 13% of genes were identified as diverged. While ~40% of the variation in hybridization ratio can be accounted for by variation in sequence identity of the heterologous sample relative to D. melanogaster, other individual characteristics of the DNA sequences, such as GC content, also contribute to variation in hybridization ratio, as does technical variation.

Conclusions

Here we demonstrate that aCGH can accurately be used as a proxy to estimate genome-wide divergence, thus providing an efficient way to evaluate how evolutionary processes and genomic architecture can shape species diversity in non-model systems. Given the increased number of species for which microarray platforms are available, comparative studies can be conducted for many interesting lineages in order to identify highly diverged genes that may be the target of natural selection.

Review of microarray studies for host-intracellular pathogen interactions

Obligate intracellular bacteria are privileged soldiers on the battlefield that represent host-pathogen interactions. Microarrays are a powerful technology that can increase our knowledge about how bacteria respond to and interact with their hosts. This review summarizes the limitations inherent to host-pathogen interaction studies and essential strategies to improve microarray investigations of intracellular bacteria. We have compiled the comparative genomic and gene expression analyses of obligate intracellular bacteria currently available from microarrays. In this review we explore ways in which microarrays can be used to identify polymorphisms in different obligate intracellular bacteria such as Coxiella burnetii, Chlamydia trachomatis, Ehrlichia chaffeensis, Rickettsia prowazekii and Tropheryma whipplei. These microarray studies reveal that, while genomic content is highly conserved in obligate intracellular bacteria, genetic polymorphisms can potentially occur to increase bacterial pathogenesis. Additionally, changes in the gene expression of C. trachomatis throughout its life cycle, as well as changes in the gene expression profile of the pathogens R. prowazekii, Rickettsia rickettsii, Rickettsia typhi, T. whipplei and C. trachomatis following environmental changes, are discussed. Finally, an in vivo model of Rickettsia conorii within the skin is discussed. The gene expression analyses highlight the capacity of obligate intracellular bacteria to adapt to environmental changes and potentially to thwart the host response.

An iterative strategy combining biophysical criteria and duration hidden Markov models for structural predictions of Chlamydia trachomatis sigma66 promoters

BACKGROUND

Promoter identification is a first step in the quest to explain gene regulation in bacteria. It has been demonstrated that the initiation of bacterial transcription depends upon the stability and topology of DNA in the promoter region as well as the binding affinity between the RNA polymerase sigma-factor and promoter. However, promoter prediction algorithms to date have not explicitly used an ensemble of these factors as predictors. In addition, most promoter models have been trained on data from Escherichia coli. Although it has been shown that transcriptional mechanisms are similar among various bacteria, it is quite possible that the differences between Escherichia coli and Chlamydia trachomatis are large enough to recommend an organism-specific modeling effort.

RESULTS

Here we present an iterative stochastic model building procedure that combines such biophysical metrics as DNA stability, curvature, twist and stress-induced DNA duplex destabilization along with duration hidden Markov model parameters to model Chlamydia trachomatis sigma66 promoters from 29 experimentally verified sequences. Initially, iterative duration hidden Markov modeling of the training set sequences provides a scoring algorithm for Chlamydia trachomatis RNA polymerase sigma66/DNA binding. Subsequently, an iterative application of Stepwise Binary Logistic Regression selects multiple promoter predictors and deletes/replaces training set sequences to determine an optimal training set. The resulting model predicts the final training set with a high degree of accuracy and provides insights into the structure of the promoter region. Model based genome-wide predictions are provided so that optimal promoter candidates can be experimentally evaluated, and refined models developed. Co-predictions with three other algorithms are also supplied to enhance reliability.

CONCLUSION

This strategy and resulting model support the conjecture that DNA biophysical properties, along with RNA polymerase sigma-factor/DNA binding collaboratively, contribute to a sequence's ability to promote transcription. This work provides a baseline model that can evolve as new Chlamydia trachomatis sigma66 promoters are identified with assistance from the provided genome-wide predictions. The proposed methodology is ideal for organisms with few identified promoters and relatively small genomes.

Co-evolution of genomes and plasmids within Chlamydia trachomatis and the emergence in Sweden of a new variant strain

Background

Chlamydia trachomatis is the most common cause of sexually transmitted infections globally and the leading cause of preventable blindness in the developing world. There are two biovariants of C. trachomatis: 'trachoma', causing ocular and genital tract infections, and the invasive 'lymphogranuloma venereum' strains. Recently, a new variant of the genital tract C. trachomatis emerged in Sweden. This variant escaped routine diagnostic tests because it carries a plasmid with a deletion. Failure to detect this strain has meant it has spread rapidly across the country provoking a worldwide alert. In addition to being a key diagnostic target, the plasmid has been linked to chlamydial virulence. Analysis of chlamydial plasmids and their cognate chromosomes was undertaken to provide insights into the evolutionary relationship between chromosome and plasmid. This is essential knowledge if the plasmid is to be continued to be relied on as a key diagnostic marker, and for an understanding of the evolution of Chlamydia trachomatis.

Results

The genomes of two new C. trachomatis strains were sequenced, together with plasmids from six C. trachomatis isolates, including the new variant strain from Sweden. The plasmid from the new Swedish variant has a 377 bp deletion in the first predicted coding sequence, abolishing the site used for PCR detection, resulting in negative diagnosis. In addition, the variant plasmid has a 44 bp duplication downstream of the deletion. The region containing the second predicted coding sequence is the most highly conserved region of the plasmids investigated. Phylogenetic analysis of the plasmids and chromosomes are fully congruent. Moreover this analysis also shows that ocular and genital strains diverged from a common C. trachomatis progenitor.

Conclusion

The evolutionary pathways of the chlamydial genome and plasmid imply that inheritance of the plasmid is tightly linked with its cognate chromosome. These data suggest that the plasmid is not a highly mobile genetic element and does not transfer readily between isolates. Comparative analysis of the plasmid sequences has revealed the most conserved regions that should be used to design future plasmid based nucleic acid amplification tests, to avoid diagnostic failures.

Evidence that CT694 is a novel Chlamydia trachomatis T3S substrate capable of functioning during invasion or early cycle development

Chlamydia trachomatis is an obligate intracellular parasite, occupies a membrane-bound vacuole throughout development and is capable of manipulating the eukaryotic host by translocating effector molecules via a type III secretion system (T3SS). The infectious chlamydial elementary body (EB) is metabolically inactive yet possesses a functional T3S apparatus capable of translocating effector proteins into the host cell to facilitate invasion and other early cycle events. We present evidence here that the C. trachomatis protein CT694 represents an early cycle-associated effector protein. CT694 is secreted by the Yersinia T3SS and immunodetection studies of infected HeLa cultures indicate that CT694-specific signal accumulates directly adjacent to, but not completely overlapping with EBs during invasion. Yeast two-hybrid analyses revealed an interaction of CT694 with the repeat region and C-terminus of human AHNAK. Immunolocalization studies of CT694 ectopically expressed in HeLa cells were consistent with an interaction with endogenous AHNAK. Additionally, expression of CT694 in HeLa cells resulted in alterations in the detection of stress fibres that correlated with the ability of CT694 to interact with AHNAK. These data indicate that CT694 is a novel T3S-dependent substrate unique to C. trachomatis, and that its interaction with host proteins such as AHNAK may be important for aspects of invasion or development particular to this species.

Chlamydia trachomatis ompA variants in trachoma: what do they tell us?

Background

Trachoma, caused by Chlamydia trachomatis (Ct), is the leading infectious cause of blindness. Sequence-based analysis of the multiple strains typically present in endemic communities may be informative for epidemiology, transmission, response to treatment, and understanding the host response.

Methods

Conjunctival and nasal samples from a Gambian community were evaluated before and 2 months after mass azithromycin treatment. Samples were tested for Ct by Amplicor, with infection load determined by quantitative PCR (qPCR). ompA sequences were determined and their diversity analysed using frequency-based tests of neutrality.

Results

Ninety-five of 1,319 (7.2%) individuals from 14 villages were infected with Ct at baseline. Two genovars (A and B) and 10 distinct ompA genotypes were detected. Two genovar A variants (A1 and A2) accounted for most infections. There was an excess of rare ompA mutations, not sustained in the population. Post-treatment, 76 (5.7%) individuals had Ct infection with only three ompA genotypes present. In 12 of 14 villages, infection had cleared, while in two it increased, probably due to mass migration. Infection qPCR loads associated with infection were significantly greater for A1 than for A2. Seven individuals had concurrent ocular and nasal infection, with divergent genotypes in five.

Conclusions

The number of strains was substantially reduced after mass treatment. One common strain was associated with higher infection loads. Discordant genotypes in concurrent infection may indicate distinct infections at ocular and nasal sites. Population genetic analysis suggests the fleeting appearance of rare multiple ompA variants represents purifying selection rather than escape variants from immune pressure. Genotyping systems accessing extra-ompA variation may be more informative.

Trachoma is an important cause of blindness resulting from transmission of the bacterium Chlamydia trachomatis. One way to understand better how this infection is transmitted and how the human immune system controls it is to study the strains of bacteria associated with infection. Comparing strains before and after treatment might help us learn if someone has a new infection or the same one as before. Identifying differences between disease-causing strains should help us understand how infection leads to disease and how the human host defences work. We chose to study variation in the chlamydial gene ompA because it determines the protein MOMP, one of the leading candidates for inclusion in a vaccine to prevent trachoma. If immunity to MOMP is important in natural trachoma infections, we would expect to find evidence of this in the way the strains varied. We did not find this, but instead found that two common strains seemed to cause different types of disease. Although their MOMPs were very slightly different, this did not really explain the differences. We conclude that methods of typing strains going beyond the ompA gene will be needed to help us understand the interaction between Chlamydia and its human host.

Comparative hybridization reveals extensive genome variation in the AIDS-associated pathogen Cryptococcus neoformans

Extensive genome variation in the AIDS-associated pathogen Cryptococcus neoformans is revealed through comparative genome hybridization between strains of different mating type, molecular subtype and ploidy.

Background

Genome variability can have a profound influence on the virulence of pathogenic microbes. The availability of genome sequences for two strains of the AIDS-associated fungal pathogen Cryptococcus neoformans presented an opportunity to use comparative genome hybridization (CGH) to examine genome variability between strains of different mating type, molecular subtype, and ploidy.

Results

Initially, CGH was used to compare the approximately 100 kilobase MATa and MATα mating-type regions in serotype A and D strains to establish the relationship between the Log2 ratios of hybridization signals and sequence identity. Subsequently, we compared the genomes of the environmental isolate NIH433 (MATa) and the clinical isolate NIH12 (MATα) with a tiling array of the genome of the laboratory strain JEC21 derived from these strains. In this case, CGH identified putative recombination sites and the origins of specific segments of the JEC21 genome. Similarly, CGH analysis revealed marked variability in the genomes of strains representing the VNI, VNII, and VNB molecular subtypes of the A serotype, including disomy for chromosome 13 in two strains. Additionally, CGH identified differences in chromosome content between three strains with the hybrid AD serotype and revealed that chromosome 1 from the serotype A genome is preferentially retained in all three strains.

Conclusion

The genomes of serotypes A, D, and AD strains exhibit extensive variation that spans the range from small differences (such as regions of divergence, deletion, or amplification) to the unexpected disomy for chromosome 13 in haploid strains and preferential retention of specific chromosomes in naturally occurring diploids.

Microarray-based approach for genome-wide survey of nucleotide polymorphisms

DNA microarrays can be used to detect polymorphic loci in addition to identifying genes or regions that are absent within a genome. A survey such as this offers greater insight into the level of diversification within a species or population, which is useful in organisms that have near-identical genomic content but differ in phenotype. The identification of such variable loci can then lead to the characterization of genes linked to unique biological attributes. Here, we describe a competitive hybridization assay using DNA microarrays as a comparative genomics tool to identify nucleotide polymorphisms among closely related strains of Chlamydia trachomatis.

Comparative genomic analysis of Tropheryma whipplei strains reveals that diversity among clinical isolates is mainly related to the WiSP proteins

Background

The aim of this study was to analyze the genomic diversity of several Tropheryma whipplei strains by microarray-based comparative genomic hybridization. Fifteen clinical isolates originating from biopsy samples recovered from different countries were compared with the T. whipplei Twist strain. For each isolate, the genes were defined as either present or absent/divergent using the GACK analysis software. Genomic changes were then further characterized by PCR and sequencing.

Results

The results revealed a limited genetic variation among the T. whipplei isolates, with at most 2.24% of the probes exhibiting differential hybridization against the Twist strain. The main variation was found in genes encoding the WiSP membrane protein family. This work also demonstrated a 19.2 kb-pair deletion within the T. whipplei DIG15 strain. This deletion occurs in the same region as the previously described large genomic rearrangement between Twist and TW08/27. Thus, this can be considered as a major hot-spot for intra-specific T. whipplei differentiation. Analysis of this deleted region confirmed the role of WND domains in generating T. whipplei diversity.

Conclusion

This work provides the first comprehensive genomic comparison of several T. whipplei isolates. It reveals that clinical isolates originating from various geographic and biological sources exhibit a high conservation rate, indicating that T. whipplei rarely interacts with exogenous DNA. Remarkably, frequent inter-strain variations were dicovered that affected members of the WiSP family.

Studying bacterial genome dynamics using microarray-based comparative genomic hybridization

Genome sequencing has revealed the remarkable amount of genetic diversity that can be encountered in bacterial genomes. In particular, the comparison of genome sequences from closely related strains has uncovered significant differences in gene content, hinting at the dynamic nature of bacterial genomes. The study of these genome dynamics is crucial to leveraging genomic information because the genome sequence of a single bacterial strain may not accurately represent the genome of the species. The dynamic nature of bacterial genome content has required us to apply the concepts of comparative genomics (CG) at the species level. Although direct genome sequence comparisons are an ideal method of performing CG, one current constraint is the limited availability of multiple genome sequences from a given bacterial species. DNA microarray-based comparative genomic hybridization (MCGH), which can be used to determine the presence or absence of thousands of genes in a single hybridization experiment, provides a powerful alternative for determining genome content and has been successfully used to investigate the genome dynamics of a wide number of bacterial species. Although MCGH-based studies have already provided a new vista on bacterial genome diversity, original methods for MCGH have been limited by the absence of novel gene sequences included in the microarray. New applications of the MCGH platform not only promise to accelerate the pace of novel gene discovery but will also help provide an integrated microarray-based approach to the study of bacterial CG.

Urogenitale Chlamydieninfektionen beim Mann

Chlamydia species are responsible for a variety of infections in humans. Chlamydia trachomatis represents the most frequently found species among Chlamydia in urogenital infections of females and males. The clinically most relevant infection appears to be bacterial adnexitis of the female that may cause sterility. However, Chlamydia also indicates a remarkably high prevalence in the urogenital tract of asymptomatic males. This review summarizes the most current developments in the diagnostics and therapeutical approaches in Chlamydia infections of the male accessory sex glands.

The Rab6 effector Bicaudal D1 associates with Chlamydia trachomatis inclusions in a biovar-specific manner

Chlamydia species are obligate intracellular bacteria that replicate within a membrane-bound vacuole, the inclusion, which is trafficked to the peri-Golgi region by processes that are dependent on early chlamydial gene expression. Although neither the host nor the chlamydial proteins that regulate the intracellular trafficking have been clearly defined, several enhanced green fluorescent protein (EGFP)-tagged Rab GTPases, including Rab6, are recruited to Chlamydia trachomatis inclusions. To further characterize the association of Rab6 with C. trachomatis inclusions, we examined the intracellular localization of guanine nucleotide-binding mutants of Rab6 and demonstrated that only active GTP-bound and not inactive GDP-bound EGFP-Rab6 mutants were recruited to the inclusion, suggesting that EGFP-Rab6 interacts with the inclusion via a host Rab6 effector or a chlamydial protein that mimics a Rab6 effector. Using EGFP-tagged fusion proteins, we also demonstrated that the Rab6 effector Bicaudal D1 (BICD1) localized to C. trachomatis inclusions in a biovar-specific manner. In addition, we demonstrated that EGFP-Rab6 and its effector EGFP-BICD1 are recruited to the inclusion in a microtubule- and Golgi apparatus-independent but chlamydial gene expression-dependent mechanism. Finally, in contrast to the Rab6-dependent Golgi apparatus localization of endogenous BICD1, EGFP-BICD1 was recruited to the inclusion by a Rab6-independent mechanism. Collectively, these data demonstrate that neither Rab6 nor BICD1 is trafficked to the inclusion via a Golgi apparatus-localized intermediate, suggesting that each protein is trafficked to the C. trachomatis serovar L2 inclusion by a unique, but as-yet-undefined, mechanism.

Evolution of Chlamydia trachomatis diversity occurs by widespread interstrain recombination involving hotspots

Chlamydia trachomatis is an obligate intracellular bacterium of major public health significance, infecting over one-tenth of the world's population and causing blindness and infertility in millions. Mounting evidence supports recombination as a key source of genetic diversity among free-living bacteria. Previous research shows that intracellular bacteria such as Chlamydiaceae may also undergo recombination but whether this plays a significant evolutionary role has not been determined. Here, we examine multiple loci dispersed throughout the chromosome to determine the extent and significance of recombination among 19 laboratory reference strains and 10 present-day ocular and urogenital clinical isolates using phylogenetic reconstructions, compatibility matrices, and statistically based recombination programs. Recombination is widespread; all clinical isolates are recombinant at multiple loci with no two belonging to the same clonal lineage. Several reference strains show nonconcordant phylogenies across loci; one strain is unambiguously identified as recombinantly derived from other reference strain lineages. Frequent recombination contrasts with a low level of point substitution; novel substitutions relative to reference strains occur less than one per kilobase. Hotspots for recombination are identified downstream from ompA, which encodes the major outer membrane protein. This widespread recombination, unexpected for an intracellular bacterium, explains why strain-typing using one or two genes, such as ompA, does not correlate with clinical phenotypes. Our results do not point to specific events that are responsible for different pathogenicities but, instead, suggest a new approach to dissect the genetic basis for clinical strain pathology with implications for evolution, host cell adaptation, and emergence of new chlamydial diseases.

Ecological and evolutionary genomics of Saccharomyces cerevisiae

Saccharomyces cerevisiae, the budding yeast, is the most thoroughly studied eukaryote at the cellular, molecular, and genetic levels. Yet, until recently, we knew very little about its ecology or population and evolutionary genetics. In recent years, it has been recognized that S. cerevisiae occupies numerous habitats and that populations harbour important genetic variation. There is therefore an increasing interest in understanding the evolutionary forces acting on the yeast genome. Several researchers have used the tools of functional genomics to study natural isolates of this unicellular fungus. Here, we review some of these studies, and show not only that budding yeast is a prime model system to address fundamental molecular and cellular biology questions, but also that it is becoming a powerful model species for ecological and evolutionary genomics studies as well.

Chlamydial infections in urology

Chlamydia trachomatis is the most frequent cause for sexually transmitted diseases in European countries. The organism has an intracellular habitat with a very specific life cycle. A variety of diagnostic tests have been developed with different sensitivity and specificity. Interpretation of these tests can sometimes be difficult. Diseases caused by C. trachomatis in men comprise urethritis, prostatitis, epididymitis, infertility and reactive arthritis. Especially in prostatitis, the exact role of C. trachomatis is still under debate for the technical difficulties localizing the pathogen to the prostate. For treatment, only some antibiotics are effective because of the intracellular habitat of the pathogen. Prevention of infection comprises treatment and screening efforts.

Comparative genomic analysis of Chlamydia trachomatis oculotropic and genitotropic strains

Chlamydia trachomatis infection is an important cause of preventable blindness and sexually transmitted disease (STD) in humans. C. trachomatis exists as multiple serovariants that exhibit distinct organotropism for the eye or urogenital tract. We previously reported tissue-tropic correlations with the presence or absence of a functional tryptophan synthase and a putative GTPase-inactivating domain of the chlamydial toxin gene. This suggested that these genes may be the primary factors responsible for chlamydial disease organotropism. To test this hypothesis, the genome of an oculotropic trachoma isolate (A/HAR-13) was sequenced and compared to the genome of a genitotropic (D/UW-3) isolate. Remarkably, the genomes share 99.6% identity, supporting the conclusion that a functional tryptophan synthase enzyme and toxin might be the principal virulence factors underlying disease organotropism. Tarp (translocated actin-recruiting phosphoprotein) was identified to have variable numbers of repeat units within the N and C portions of the protein. A correlation exists between lymphogranuloma venereum serovars and the number of N-terminal repeats. Single-nucleotide polymorphism (SNP) analysis between the two genomes highlighted the minimal genetic variation. A disproportionate number of SNPs were observed within some members of the polymorphic membrane protein (pmp) autotransporter gene family that corresponded to predicted T-cell epitopes that bind HLA class I and II alleles. These results implicate Pmps as novel immune targets, which could advance future chlamydial vaccine strategies. Lastly, a novel target for PCR diagnostics was discovered that can discriminate between ocular and genital strains. This discovery will enhance epidemiological investigations in nations where both trachoma and chlamydial STD are endemic.

Large-scale genetic variation of the symbiosis-required megaplasmid pSymA revealed by comparative genomic analysis of Sinorhizobium meliloti natural strains

Background

Sinorhizobium meliloti is a soil bacterium that forms nitrogen-fixing nodules on the roots of leguminous plants such as alfalfa (Medicago sativa). This species occupies different ecological niches, being present as a free-living soil bacterium and as a symbiont of plant root nodules. The genome of the type strain Rm 1021 contains one chromosome and two megaplasmids for a total genome size of 6 Mb. We applied comparative genomic hybridisation (CGH) on an oligonucleotide microarrays to estimate genetic variation at the genomic level in four natural strains, two isolated from Italian agricultural soil and two from desert soil in the Aral Sea region.

Results

From 4.6 to 5.7 percent of the genes showed a pattern of hybridisation concordant with deletion, nucleotide divergence or ORF duplication when compared to the type strain Rm 1021. A large number of these polymorphisms were confirmed by sequencing and Southern blot. A statistically significant fraction of these variable genes was found on the pSymA megaplasmid and grouped in clusters. These variable genes were found to be mainly transposases or genes with unknown function.

Conclusion

The obtained results allow to conclude that the symbiosis-required megaplasmid pSymA can be considered the major hot-spot for intra-specific differentiation in S. meliloti.

Characterization of the genome composition of Bartonella koehlerae by microarray comparative genomic hybridization profiling

Bartonella henselae is present in a wide range of wild and domestic feline hosts and causes cat-scratch disease and bacillary angiomatosis in humans. We have estimated here the gene content of Bartonella koehlerae, a novel species isolated from cats that was recently identified as an agent of human endocarditis. The investigation was accomplished by comparative genomic hybridization (CGH) to a microarray constructed from the sequenced 1.93-Mb genome of B. henselae. Control hybridizations of labeled DNA from the human pathogen Bartonella quintana with a reduced genome of 1.58 Mb were performed to evaluate the accuracy of the array for genes with known levels of sequence divergence. Genome size estimates of B. koehlerae by pulsed-field gel electrophoresis matched that calculated by the CGH, indicating a genome of 1.7 to 1.8 Mb with few unique genes. As in B. quintana, sequences in the prophage and the genomic islands were reported absent in B. koehlerae. In addition, sequence variability was recorded in the chromosome II-like region, where B. koehlerae showed an intermediate retention pattern of both coding and noncoding sequences. Although most of the genes missing in B. koehlerae are also absent from B. quintana, its phylogenetic placement near B. henselae suggests independent deletion events, indicating that host specificity is not solely attributed to genes in the genomic islands. Rather, the results underscore the instability of the genomic islands even within bacterial populations adapted to the same host-vector system, as in the case of B. henselae and B. koehlerae.

Polymorphisms in the Chlamydia trachomatis cytotoxin locus associated with ocular and genital isolates

Chlamydia trachomatis is a strict human pathogen producing infections that cause medically important chronic inflammatory diseases, such as blinding trachoma and tubal factor infertility. Isolates exist as serotypes that fall into distinct biologic and pathological groups corresponding to differences in infection tissue tropism and invasion properties. Paradoxically, genome sequencing of several diverse strains has revealed a remarkable level of genomic synteny, suggesting that minor genetic differences determine the pathogen host- and tissue-specific infection characteristics. To better understand the genetic basis of chlamydial pathobiologic diversity, we performed comparative DNA-DNA microarray genomic hybridizations with all 15 C. trachomatis serovariants. We found there are few major genetic differences among the 15 serovars. An exception was the cytotoxin locus located in the plasticity zone, a region that exhibited significant polymorphisms among serovars. We therefore sequenced this region from all 15 serovars. The cytotoxin gene was interrupted by extensive mutations and deletions among the different serovars; however, three basic open reading frame motifs were discovered that correlated with noninvasive oculotropic, urogenitotropic, and invasive serovars. Of interest, only noninvasive genitotropic serovars possessed an intact N-terminal portion of the putative toxin gene. This region contains the UDP-glucose binding domain and the glycosyltransferase domain required for enzymatic activity of the clostridial toxin homologs, suggesting a role in urogenital infection or pathogenesis.