Novel Chlamydiales strains isolated from a water treatment plant.

Chlamydiae are obligate intracellular bacteria infecting free-living amoebae, vertebrates and some invertebrates. Novel members are regularly discovered, and there is accumulating evidence supporting a very important diversity of chlamydiae in the environment. In this study, we investigated the presence of chlamydiae in a drinking water treatment plant. Samples were used to inoculate Acanthamoeba monolayers (Acanthamoeba co-culture), and to recover autochthonous amoebae onto non-nutritive agar. Chlamydiae were searched for by a pan-chlamydia 16S rRNA gene PCR from both Acanthamoeba co-cultures and autochthonous amoebae, and phylotypes determined by 16S rRNA gene sequencing. Autochthonous amoebae also were identified by 18S rRNA gene amplification and sequencing. From a total of 79 samples, we recovered eight chlamydial strains by Acanthamoeba co-culture, but only one of 28 amoebae harboured a chlamydia. Sequencing results and phylogenetic analysis showed our strains belonging to four distinct chlamydial lineages. Four strains, including the strain recovered within its natural host, belonged to the Parachlamydiaceae; two closely related strains belonged to the Criblamydiaceae; two distinct strains clustered with Rhabdochlamydia spp.; one strain clustered only with uncultured environmental clones. Our results confirmed the usefulness of amoeba co-culture to recover novel chlamydial strains from complex samples and demonstrated the huge diversity of chlamydiae in the environment, by identifying several new species including one representing the first strain of a new family.

Evolutionary simulations to detect functional lineage-specific genes.

MOTIVATION: Supporting the functionality of recent duplicate gene copies is usually difficult, owing to high sequence similarity between duplicate counterparts and shallow phylogenies, which hamper both the statistical and experimental inference. RESULTS: We developed an integrated evolutionary approach to identify functional duplicate gene copies and other lineage-specific genes. By repeatedly simulating neutral evolution, our method estimates the probability that an ORF was selectively conserved and is therefore likely to represent a bona fide coding region. In parallel, our method tests whether the accumulation of non-synonymous substitutions reveals signatures of selective constraint. We show that our approach has high power to identify functional lineage-specific genes using simulated and real data. For example, a coding region of average length (approximately 1400 bp), restricted to hominoids, can be predicted to be functional in approximately 94-100% of cases. Notably, the method may support functionality for instances where classical selection tests based on the ratio of non-synonymous to synonymous substitutions fail to reveal signatures of selection. Our method is available as an automated tool, ReEVOLVER, which will also be useful to systematically detect functional lineage-specific genes of closely related species on a large scale. AVAILABILITY: ReEVOLVER is available at http://www.unil.ch/cig/page7858.html.

False phylogenies on wood mice due to cryptic cytochrome-b pseudogene.

The phylogeny and phylogeography of the Old World wood mice (subgenus Sylvaemus, genus Apodemus, Muridae) are well-documented. Nevertheless, the distributions of species, such as A. fulvipectus and A. ponticus remain dubious, as well as their phylogenetic relationships with A. sylvaticus. We analysed samples of Apodemus spp. across Europe using the mitochondrial cytochrome-b gene (cyt-b) and compared the DNA and amino-acid compositions of previously published sequences. The main result stemming from this study is the presence of a well-differentiated lineage of Sylvaemus including samples of various species (A. sylvaticus, A. fulvipectus, A. ponticus) from distant locations, which were revealed to be nuclear copies of the mitochondrial cyt-b. The presence of this cryptic pseudogene in published sequences is supported by different pathways. This has led to important errors in previous molecular trees and hence to partial misinterpretations in the phylogeny of Apodemus.

Discovery of catalases in members of the Chlamydiales order.

Catalase is an important virulence factor for survival in macrophages and other phagocytic cells. In Chlamydiaceae, no catalase had been described so far. With the sequencing and annotation of the full genomes of Chlamydia-related bacteria, the presence of different catalase-encoding genes has been documented. However, their distribution in the Chlamydiales order and the functionality of these catalases remain unknown. Phylogeny of chlamydial catalases was inferred using MrBayes, maximum likelihood, and maximum parsimony algorithms, allowing the description of three clade 3 and two clade 2 catalases. Only monofunctional catalases were found (no catalase-peroxidase or Mn-catalase). All presented a conserved catalytic domain and tertiary structure. Enzymatic activity of cloned chlamydial catalases was assessed by measuring hydrogen peroxide degradation. The catalases are enzymatically active with different efficiencies. The catalase of Parachlamydia acanthamoebae is the least efficient of all (its catalytic activity was 2 logs lower than that of Pseudomonas aeruginosa). Based on the phylogenetic analysis, we hypothesize that an ancestral class 2 catalase probably was present in the common ancestor of all current Chlamydiales but was retained only in Criblamydia sequanensis and Neochlamydia hartmannellae. The catalases of class 3, present in Estrella lausannensis and Parachlamydia acanthamoebae, probably were acquired by lateral gene transfer from Rhizobiales, whereas for Waddlia chondrophila they likely originated from Legionellales or Actinomycetales. The acquisition of catalases on several occasions in the Chlamydiales suggests the importance of this enzyme for the bacteria in their host environment.

Arbuscular mycorrhizal fungi (Glomeromycota) harbour ancient fungal tubulin genes that resemble those of the chytrids (Chytridiomycota).

The genes encoding alpha- and beta-tubulins have been widely sampled in most major fungal phyla and they are useful tools for fungal phylogeny. Here, we report the first isolation of alpha-tubulin sequences from arbuscular mycorrhizal fungi (AMF). In parallel, AMF beta-tubulins were sampled and analysed to identify the presence of paralogs of this gene. The AMF alpha-tubulin amino acid phylogeny was congruent with the results previously reported for AMF beta-tubulins and showed that AMF tubulins group together at a basal position in the fungal clade and showed high sequence similarities with members of the Chytridiomycota. This is in contrast with phylogenies for other regions of the AMF genome. The amount and nature of substitutions are consistent with an ancient divergence of both orthologs and paralogs of AMF tubulins. At the amino acid level, however, AMF tubulins have hardly evolved from those of the chytrids. This is remarkable given that these two groups are ancient and the monophyletic Glomeromycota probably diverged from basal fungal ancestors at least 500 million years ago. The specific primers we designed for the AMF tubulins, together with the high molecular variation we found among the AMF species we analysed, make AMF tubulin sequences potentially useful for AMF identification purposes.

Genetic consequences of the ice ages on nurseries of the bat Myotis myotis: a mitochondrial and nuclear survey.

Analyses of mitochondrial DNA (mtDNA) control region polymorphism and of variation at 10 nuclear microsatellite loci were used to investigate the mechanisms and genetic consequences of postglacial expansion of Myotis myotis in Europe. Initial sampling consisted of 480 bats genotyped in 24 nursery colonies arranged along a transect of approximately 3000 km. The phylogeographical survey based on mtDNA sequences revealed the existence of major genetic subdivisions across this area, with several suture zones between haplogroups. Such zones of secondary contact were found in the Alps and Rhodopes, whereas other potential barriers to gene flow, like the Pyrenees, did not coincide with genetic discontinuities. Areas of population admixture increased locally the genetic diversity of colonies, which confounded the northward decrease in nucleotide diversity predicted using classical models of postglacial range expansion. However, when analyses were restricted to a subset of 15 nurseries originating from a single presumed glacial refugium, mtDNA polymorphism did indeed support a northwards decrease in diversity. Populations were also highly structured (PhiST = 0.384). Conversely, the same subset of colonies showed no significant latitudinal decrease in microsatellite diversity and much less population structure (FST = 0.010), but pairwise genetic differentiation at these nuclear markers was strongly correlated with increasing geographical distance. Together, this evidence suggests that alleles carried via male bats have maintained enough nuclear gene flow to counteract the effects of recurrent bottlenecks generally associated with recolonization processes. As females are highly philopatric, we argue that the maternally transmitted mtDNA marker better reflects the situation of past, historical gene flow, whereas current levels of gene flow are better reflected by microsatellite markers.

Low rates of X-Y recombination, not turnovers, account for homomorphic sex chromosomes in several diploid species of Palearctic green toads (Bufo viridis subgroup).

Contrasting with birds and mammals, most ectothermic vertebrates present homomorphic sex chromosomes, which might be due either to a high turnover rate or to occasional X-Y recombination. We tested these two hypotheses in a group of Palearctic green toads that diverged some 3.3 million years ago. Using sibship analyses of sex-linked markers, we show that all four species investigated share the same pair of sex chromosomes and a pattern of male heterogamety with drastically reduced X-Y recombination in males. Phylogenetic analyses of sex-linked sequences show that X and Y alleles cluster by species, not by gametolog. We conclude that X-Y homomorphy and fine-scale sequence similarity in these species do not stem from recent sex-chromosome turnovers, but from occasional X-Y recombination.

Evolution of Trypanosoma cruzi: clarifying hybridisations, mitochondrial introgressions and phylogenetic relationships between major lineages

Several different models of Trypanosoma cruzi evolution have been proposed. These models suggest that scarce events of genetic exchange occurred during the evolutionary history of this parasite. In addition, the debate has focused on the existence of one or two hybridisation events during the evolution of T. cruzi lineages. Here, we reviewed the literature and analysed available sequence data to clarify the phylogenetic relationships among these different lineages. We observed that TcI, TcIII and TcIV form a monophyletic group and that TcIII and TcIV are not, as previously suggested, TcI-TcII hybrids. Particularly, TcI and TcIII are sister groups that diverged around the same time that a widely distributed TcIV split into two clades (TcIVS and TcIVN). In addition, we collected evidence that TcIII received TcIVSkDNA by introgression on several occasions. Different demographic hypotheses (surfing and asymmetrical introgression) may explain the origin and expansion of the TcIII group. Considering these hypotheses, genetic exchange should have been relatively frequent between TcIII and TcIVS in the geographic area in which their distributions overlapped. In addition, our results support the hypothesis that two independent hybridisation events gave rise to TcV and TcVI. Consequently, TcIVS kDNA was first transferred to TcIII and later to TcV and TcVI in TcII/TcIII hybridisation events.

Genetic diversity and host plant preferences revealed by simple sequence repeat and mitochondrial markers in a population of the arbuscular mycorrhizal fungus Glomus intraradices.

Arbuscular mycorrhizal fungi (AMF) are important symbionts of plants that improve plant nutrient acquisition and promote plant diversity. Although within-species genetic differences among AMF have been shown to differentially affect plant growth, very little is actually known about the degree of genetic diversity in AMF populations. This is largely because of difficulties in isolation and cultivation of the fungi in a clean system allowing reliable genotyping to be performed. A population of the arbuscular mycorrhizal fungus Glomus intraradices growing in an in vitro cultivation system was studied using newly developed simple sequence repeat (SSR), nuclear gene intron and mitochondrial ribosomal gene intron markers. The markers revealed a strong differentiation at the nuclear and mitochondrial level among isolates. Genotypes were nonrandomly distributed among four plots showing genetic subdivisions in the field. Meanwhile, identical genotypes were found in geographically distant locations. AMF genotypes showed significant preferences to different host plant species (Glycine max, Helianthus annuus and Allium porrum) used before the fungal in vitro culture establishment. Host plants in a field could provide a heterogeneous environment favouring certain genotypes. Such preferences may partly explain within-population patterns of genetic diversity.

Contribution of efflux pumps, porins, and β-lactamases to multidrug resistance in clinical isolates of Acinetobacter baumannii.

We investigated the mechanisms of resistance to carbapenems, aminoglycosides, glycylcyclines, tetracyclines, and quinolones in 90 multiresistant clinical strains of Acinetobacter baumannii isolated from two genetically unrelated A. baumannii clones: clone PFGE-ROC-1 (53 strains producing the OXA-58 β-lactamase enzyme and 18 strains with the OXA-24 β-lactamase) and clone PFGE-HUI-1 (19 strains susceptible to carbapenems). We used real-time reverse transcriptase PCR to correlate antimicrobial resistance (MICs) with expression of genes encoding chromosomal β-lactamases (AmpC and OXA-51), porins (OmpA, CarO, Omp33, Dcap-like, OprB, Omp25, OprC, OprD, and OmpW), and proteins integral to six efflux systems (AdeABC, AdeIJK, AdeFGH, CraA, AbeM, and AmvA). Overexpression of the AdeABC system (level of expression relative to that by A. baumannii ATCC 17978, 30- to 45-fold) was significantly associated with resistance to tigecycline, minocycline, and gentamicin and other biological functions. However, hyperexpression of the AdeIJK efflux pump (level of expression relative to that by A. baumannii ATCC 17978, 8- to 10-fold) was significantly associated only with resistance to tigecycline and minocycline (to which the TetB efflux system also contributed). TetB and TetA(39) efflux pumps were detected in clinical strains and were associated with resistance to tetracyclines and doxycycline. The absence of the AdeABC system and the lack of expression of other mechanisms suggest that tigecycline-resistant strains of the PFGE-HUI-1 clone may be associated with a novel resistance-nodulation-cell efflux pump (decreased MICs in the presence of the inhibitor Phe-Arg β-naphthylamide dihydrochloride) and the TetA(39) system.

Describing ancient horizontal gene transfers at the nucleotide and gene levels by comparative pathogenicity island genometrics.

MOTIVATION: Lateral gene transfer is a major mechanism contributing to bacterial genome dynamics and pathovar emergence via pathogenicity island (PAI) spreading. However, since few of these genomic exchanges are experimentally reproducible, it is difficult to establish evolutionary scenarios for the successive PAI transmissions between bacterial genera. Methods initially developed at the gene and/or nucleotide level for genomics, i.e. comparisons of concatenated sequences, ortholog frequency, gene order or dinucleotide usage, were combined and applied here to homologous PAIs: we call this approach comparative PAI genometrics. RESULTS: YAPI, a Yersinia PAI, and related islands were compared with measure evolutionary relationships between related modules. Through use of our genometric approach designed for tracking codon usage adaptation and gene phylogeny, an ancient inter-genus PAI transfer was oriented for the first time by characterizing the genomic environment in which the ancestral island emerged and its subsequent transfers to other bacterial genera.

Ever-young sex chromosomes in European tree frogs.

Non-recombining sex chromosomes are expected to undergo evolutionary decay, ending up genetically degenerated, as has happened in birds and mammals. Why are then sex chromosomes so often homomorphic in cold-blooded vertebrates? One possible explanation is a high rate of turnover events, replacing master sex-determining genes by new ones on other chromosomes. An alternative is that X-Y similarity is maintained by occasional recombination events, occurring in sex-reversed XY females. Based on mitochondrial and nuclear gene sequences, we estimated the divergence times between European tree frogs (Hyla arborea, H. intermedia, and H. molleri) to the upper Miocene, about 5.4-7.1 million years ago. Sibship analyses of microsatellite polymorphisms revealed that all three species have the same pair of sex chromosomes, with complete absence of X-Y recombination in males. Despite this, sequences of sex-linked loci show no divergence between the X and Y chromosomes. In the phylogeny, the X and Y alleles cluster according to species, not in groups of gametologs. We conclude that sex-chromosome homomorphy in these tree frogs does not result from a recent turnover but is maintained over evolutionary timescales by occasional X-Y recombination. Seemingly young sex chromosomes may thus carry old-established sex-determining genes, a result at odds with the view that sex chromosomes necessarily decay until they are replaced. This raises intriguing perspectives regarding the evolutionary dynamics of sexually antagonistic genes and the mechanisms that control X-Y recombination.

Phylogenetic approach reveals that virus genotype largely determines HIV set-point viral load.

HIV virulence, i.e. the time of progression to AIDS, varies greatly among patients. As for other rapidly evolving pathogens of humans, it is difficult to know if this variance is controlled by the genotype of the host or that of the virus because the transmission chain is usually unknown. We apply the phylogenetic comparative approach (PCA) to estimate the heritability of a trait from one infection to the next, which indicates the control of the virus genotype over this trait. The idea is to use viral RNA sequences obtained from patients infected by HIV-1 subtype B to build a phylogeny, which approximately reflects the transmission chain. Heritability is measured statistically as the propensity for patients close in the phylogeny to exhibit similar infection trait values. The approach reveals that up to half of the variance in set-point viral load, a trait associated with virulence, can be heritable. Our estimate is significant and robust to noise in the phylogeny. We also check for the consistency of our approach by showing that a trait related to drug resistance is almost entirely heritable. Finally, we show the importance of taking into account the transmission chain when estimating correlations between infection traits. The fact that HIV virulence is, at least partially, heritable from one infection to the next has clinical and epidemiological implications. The difference between earlier studies and ours comes from the quality of our dataset and from the power of the PCA, which can be applied to large datasets and accounts for within-host evolution. The PCA opens new perspectives for approaches linking clinical data and evolutionary biology because it can be extended to study other traits or other infectious diseases.

Disease in the dark: molecular characterization of Polychromophilus murinus in temperate zone bats revealed a worldwide distribution of this malaria-like disease.

For a better understanding of the complex coevolutionary processes between hosts and parasites, accurate identification of the actors involved in the interaction is of fundamental importance. Blood parasites of the Order Haemosporidia, responsible for malaria, have become the focus of a broad range of studies in evolutionary biology. Interestingly, molecular-based studies on avian malaria have revealed much higher species diversity than previously inferred with morphology. Meanwhile, studies on bat haemosporidian have been largely neglected. In Europe, only one genus (Polychromophilus) and two species have been morphologically described. To evaluate the presence of potential cryptic species and parasite prevalence, we undertook a molecular characterization of Polychromophilus in temperate zone bats. We used a nested-PCR approach on the cytochrome b mitochondrial gene to detect the presence of parasites in 237 bats belonging to four different species and in the dipteran bat fly Nycteribia kolenatii, previously described as being the vector of Polychromophilus. Polychromophilus murinus was found in the four bat species and in the insect vector with prevalence ranging from 4% for Myotis myotis to 51% for M. daubentoni. By sequencing 682 bp, we then investigated the phylogenetic relationships of Polychromophilus to other published malarial lineages. Seven haplotypes were found, all very closely related, suggesting the presence of a single species in our samples. These haplotypes formed a well-defined clade together with Haemosporidia of tropical bats, revealing a worldwide distribution of this parasite mostly neglected by malarial studies since the 1980s.

The role of the vascular endothelium in arenavirus haemorrhagic fevers.

Viral haemorrhagic fevers (VHF) caused by arenaviruses are among the most devastating emerging human diseases. The most important pathogen among the arenaviruses is Lassa virus (LASV), the causative agent of Lassa fever that is endemic to West Africa. On the South American continent, the New World arenavirus Junin virus (JUNV), Machupo (MACV), Guanarito (GTOV), and Sabia virus (SABV) have emerged as causative agents of severe VHFs. Clinical and experimental studies on arenavirus VHF have revealed a crucial role of the endothelium in their pathogenesis. However, in contrast to other VHFs, haemorrhages are not a salient feature of Lassa fever and fatal cases do not show overt destruction of vascular tissue. The functional alteration of the vascular endothelium that precede shock and death in fatal Lassa fever may be due to more subtle direct or indirect effects of the virus on endothelial cells. Haemorrhagic disease manifestations and vascular involvement are more pronounced in the VHF caused by the South American haemorrhagic fever viruses. Recent studies on JUNV revealed perturbation of specific endothelial cell function, including expression of cell adhesion molecules, coagulation factors, and vasoactive mediators as a consequence of productive viral infection. These studies provided first possible links to some of the vascular abnormalities observed in patients, however, their relevance in vivo remains to be investigated.