Differentiation and phylogenetic relationships in Mycobacterium spp with special reference to the RNase P RNA gene rnpB.

The rnpB gene encodes for the RNA subunit of the catalytic ribonuclease RNase P and is present in all bacteria and has both conserved and highly variable sequence regions. Determination of rnpB in 35 Mycobacterium spp. showed species specific sequences for all species except the Mycobacterium tuberculosis complex (four species). High sequence variation was seen in the P3, P15 and P19 regions of suggested secondary structures of the corresponding RNase P RNA molecules. Phylogenetic analysis showed that rnpB gave similar tree topologies as 16S rRNA and hsp65 genes. A combined analysis of the three genes increased the number of nodes with significant support from 10 to 19. The results indicate that rnpB is useful for phylogenetic studies and is a possible target for identification and detection of Mycobacterium spp.

SPRIT: Identifying horizontal gene transfer in rooted phylogenetic trees.

BACKGROUND: Phylogenetic trees based on sequences from a set of taxa can be incongruent due to horizontal gene transfer (HGT). By identifying the HGT events, we can reconcile the gene trees and derive a taxon tree that adequately represents the species' evolutionary history. One HGT can be represented by a rooted Subtree Prune and Regraft (RSPR) operation and the number of RSPRs separating two trees corresponds to the minimum number of HGT events. Identifying the minimum number of RSPRs separating two trees is NP-hard, but the problem can be reduced to fixed parameter tractable. A number of heuristic and two exact approaches to identifying the minimum number of RSPRs have been proposed. This is the first implementation delivering an exact solution as well as the intermediate trees connecting the input trees. RESULTS: We present the SPR Identification Tool (SPRIT), a novel algorithm that solves the fixed parameter tractable minimum RSPR problem and its GPL licensed Java implementation. The algorithm can be used in two ways, exhaustive search that guarantees the minimum RSPR distance and a heuristic approach that guarantees finding a solution, but not necessarily the minimum one. We benchmarked SPRIT against other software in two different settings, small to medium sized trees i.e. five to one hundred taxa and large trees i.e. thousands of taxa. In the small to medium tree size setting with random artificial incongruence, SPRIT's heuristic mode outperforms the other software by always delivering a solution with a low overestimation of the RSPR distance. In the large tree setting SPRIT compares well to the alternatives when benchmarked on finding a minimum solution within a reasonable time. SPRIT presents both the minimum RSPR distance and the intermediate trees. CONCLUSIONS: When used in exhaustive search mode, SPRIT identifies the minimum number of RSPRs needed to reconcile two incongruent rooted trees. SPRIT also performs quick approximations of the minimum RSPR distance, which are comparable to, and often better than, purely heuristic solutions. Put together, SPRIT is an excellent tool for identification of HGT events and pinpointing which taxa have been involved in HGT.

Mosaic structure of intragenic repetitive elements in histone H1-like protein Hc2 varies within serovars of Chlamydia trachomatis.

BACKGROUND: The histone-like protein Hc2 binds DNA in Chlamydia trachomatis and is known to vary in size between 165 and 237 amino acids, which is caused by different numbers of lysine-rich pentamers. A more complex structure was seen in this study when sequences from 378 specimens covering the hctB gene, which encodes Hc2, were compared. RESULTS: This study shows that the size variation is due to different numbers of 36-amino acid long repetitive elements built up of five pentamers and one hexamer. Deletions and amino acid substitutions result in 14 variants of repetitive elements and these elements are combined into 22 configurations. A protein with similar structure has been described in Bordetella but was now also found in other genera, including Burkholderia, Herminiimonas, Minibacterium and Ralstonia.Sequence determination resulted in 41 hctB variants that formed four clades in phylogenetic analysis. Strains causing the eye disease trachoma and strains causing invasive lymphogranuloma venereum infections formed separate clades, while strains from urogenital infections were more heterogeneous. Three cases of recombination were identified. The size variation of Hc2 has previously been attributed to deletions of pentamers but we show that the structure is more complex with both duplication and deletions of 36-amino acid long elements. CONCLUSIONS: The polymorphisms in Hc2 need to be further investigated in experimental studies since DNA binding is essential for the unique biphasic life cycle of the Chlamydiacae. The high sequence variation in the corresponding hctB gene enables phylogenetic analysis and provides a suitable target for the genotyping of C. trachomatis.

Phylogeny and species delimitation in Silene sect. Arenosae (Caryophyllaceae): a new section.

A putatively monophyletic group of annual Silene species is revised taxonomically and described as the new section S. sect. Arenosae. The species of this section were previously treated as a part of a widely circumscribed and polyphyletic S. sect. Rigidulae. Silene sect. Arenosae as circumscribed here consists of nine species. Members of the section show a predominantly E Mediterranean to SW Asian distribution pattern from Turkey southward to Egypt and eastward to Iran and Pakistan, although most of the species have a limited distribution range. The species of S. sect. Arenosae are characterized by narrowly lanceolate calyx teeth, which are often highly polymorphic, and lanceolate to oblanceolate (non-spathulate) basal leaves. The provided taxonomic revision is based on morphological characters and supported by phylogenetic analyses of two nuclear loci (nrITS and an intron of the RPB2 gene) and one chloroplast locus (the intron of the rps16 gene). The species descriptions are formalized using a novel implementation of the Prometheus Description Model.

Diversifying selection and concerted evolution of a type IV secretion system in Bartonella.

We have studied the evolution of a type IV secretion system (T4SS), in Bartonella, which is thought to have changed function from conjugation to erythrocyte adherence following a recent horizontal gene transfer event. The system, called Trw, is unique among T4SSs in that genes encoding both exo- and intracellular components are located within the same duplicated fragment. This provides an opportunity to study the influence of selection on proteins involved in host-pathogen interactions. We sequenced the trw locus from several strains of Bartonella henselae and investigated its evolutionary history by comparisons to other Bartonella species. Several instances of recombination and gene conversion events where detected in the 2- to 5-fold duplicated gene fragments encompassing trwJIH, explaining the homogenization of the anchoring protein TrwI and the divergence of the minor pilus protein TrwJ. A phylogenetic analysis of the 7- to 8-fold duplicated gene coding for the major pilus protein TrwL displayed 2 distinct clades, likely representing a subfunctionalization event. The analyses of the B. henselae strains also identified a recent horizontal transfer event of almost the complete trwL region. We suggest that the switch in function of the T4SS was mediated by the duplication of the genes encoding pilus components and their diversification by combinatorial sequence shuffling within and among genomes. We suggest that the pilus proteins have evolved by diversifying selection to match a divergent set of erythrocyte surface structures, consistent with the trench warfare coevolutionary model.

Ribbon worm relationships: a phylogeny of the phylum Nemertea.

We present the most extensive phylogenetic analysis to date, to our knowledge, of higher-level nemertean relationships, based on sequence data from four different genes (the nuclear genes for nuclear large subunit rRNA (28S rRNA) and histone H3 (H3), and the mitochondrial genes for mitochondrial large subunit rRNA (16S rRNA) and cytochrome c oxidase subunit I (COI)). Well-supported clades are, in general, compatible with earlier, more limited, analyses, and current classification is largely in agreement with our results, although there are some notable exceptions. Bdellonemertea (represented by Malacobdella) is found to be a part of Monostilifera, and Polystilifera is the monophyletic sister group to Monostilifera. Cratenemertidae is the sister group to the remaining monostiliferans (including Malacobdella), a group to which we apply the new name Distromatonemertea. Heteronemertea is monophyletic and forms a clade with Hubrechtella; for this clade we introduce the name Pilidiophora. Finally, Pilidiophora and Hoplonemertea (with Malacobdella) form a monophyletic group, and we introduce the name Neonemertea to refer to this group. Palaeonemertea is found to be non-monophyletic and basal among nemerteans.

The phylogenetic position of the comb jellies (Ctenophora) and the importance of taxonomic sampling.

The transition to a vermiform body shape is one of the most important events in animal evolution, having led to the impressive radiation of Bilateria. However, the sister group of Bilateria has remained obscure. Cladistic analyses of morphology indicate that Ctenophora is the sister group of Bilateria. Previous analyses of SSU rRNA sequences have yielded conflicting results; in many studies Ctenophora forms the sister group of Cnidaria + Bilateria, but in others the ctenophores group with poriferans. Here we re-examine the SSU sequence by analyzing a dataset with 528 metazoan + outgroup sequences, including almost 120 poriferan and diploblast sequences. We use parsimony ratchet and jackknife methods, as well as Bayesian methods, to analyze the data. The results indicate strong phylogenetic signals for a cnidarian + bilaterian group and for the comb jellies to have branched off early within a group uniting all epithelial animals [(Ct,(Cn,Bi))]. We demonstrate the importance of inclusive taxonomic coverage of ribosomal sequences for resolving this problematic part of the metazoan tree: topological stability increases dramatically with the addition of taxa, and the jackknife frequencies of the internal nodes uniting the lineages [(Cn,Bi) and ((Ct,(Cn,Bi))] also increase. We consider the reconstructed topology to represent the current best hypothesis of the interrelationships of these old lineages. Some morphological features supporting alternative hypotheses are discussed in the light of this result.

A phylometagenomic exploration of oceanic alphaproteobacteria reveals mitochondrial relatives unrelated to the SAR11 clade.

BACKGROUND: According to the endosymbiont hypothesis, the mitochondrial system for aerobic respiration was derived from an ancestral Alphaproteobacterium. Phylogenetic studies indicate that the mitochondrial ancestor is most closely related to the Rickettsiales. Recently, it was suggested that Candidatus Pelagibacter ubique, a member of the SAR11 clade that is highly abundant in the oceans, is a sister taxon to the mitochondrial-Rickettsiales clade. The availability of ocean metagenome data substantially increases the sampling of Alphaproteobacteria inhabiting the oxygen-containing waters of the oceans that likely resemble the originating environment of mitochondria. METHODOLOGY/PRINCIPAL FINDINGS: We present a phylogenetic study of the origin of mitochondria that incorporates metagenome data from the Global Ocean Sampling (GOS) expedition. We identify mitochondrially related sequences in the GOS dataset that represent a rare group of Alphaproteobacteria, designated OMAC (Oceanic Mitochondria Affiliated Clade) as the closest free-living relatives to mitochondria in the oceans. In addition, our analyses reject the hypothesis that the mitochondrial system for aerobic respiration is affiliated with that of the SAR11 clade. CONCLUSIONS/SIGNIFICANCE: Our results allude to the existence of an alphaproteobacterial clade in the oxygen-rich surface waters of the oceans that represents the closest free-living relative to mitochondria identified thus far. In addition, our findings underscore the importance of expanding the taxonomic diversity in phylogenetic analyses beyond that represented by cultivated bacteria to study the origin of mitochondria.

A novel Chlamydiaceae-like bacterium found in faecal specimens from sea birds from the Bering Sea.

The family Chlamydiaceae contains several bacterial pathogens of important human and veterinary medical concern, such as Chlamydia trachomatis and Chlamydophila psittaci. Within the order Chlamydiales there are also an increasing number of chlamydia-like bacteria whose biodiversity, host range and environmental spread seem to have been largely underestimated, and which are currently being investigated for their potential medical relevance. In this study we present 16S rRNA, rnpB and ompA gene sequence data congruently indicating a novel chlamydia-like bacterium found in faecal specimens from opportunistic fish-eating sea birds, belonging to the Laridae and Alcidae families, from the Bering Sea. This novel bacterium appears to be closer to the Chlamydiaceae than other chlamydia-like bacteria and is most likely a novel genus within the Chlamydiaceae family.

Origin and evolution of the mitochondrial aminoacyl-tRNA synthetases.

Many theories favor a fusion of 2 prokaryotic genomes for the origin of the Eukaryotes, but there are disagreements on the origin, timing, and cellular structures of the cells involved. Equally controversial is the source of the nuclear genes for mitochondrial proteins, although the alpha-proteobacterial contribution to the mitochondrial genome is well established. Phylogenetic inferences show that the nuclearly encoded mitochondrial aminoacyl-tRNA synthetases (aaRSs) occupy a position in the tree that is not close to any of the currently sequenced alpha-proteobacterial genomes, despite cohesive and remarkably well-resolved alpha-proteobacterial clades in 12 of the 20 trees. Two or more alpha-proteobacterial clusters were observed in 8 cases, indicative of differential loss of paralogous genes or horizontal gene transfer. Replacement and retargeting events within the nuclear genomes of the Eukaryotes was indicated in 10 trees, 4 of which also show split alpha-proteobacterial groups. A majority of the mitochondrial aaRSs originate from within the bacterial domain, but none specifically from the alpha-Proteobacteria. For some aaRS, the endosymbiotic origin may have been erased by ongoing gene replacements on the bacterial as well as the eukaryotic side. For others that accurately resolve the alpha-proteobacterial divergence patterns, the lack of affiliation with mitochondria is more surprising. We hypothesize that the ancestral eukaryotic gene pool hosted primordial "bacterial-like" genes, to which a limited set of alpha-proteobacterial genes, mostly coding for components of the respiratory chain complexes, were added and selectively maintained.

Functional divergence and horizontal transfer of type IV secretion systems.

The type IV secretion system (TFSSs) is a multifunctional family of translocation pathways that mediate the transfer of DNA among bacteria and deliver DNA and proteins to eukaryotic cells during bacterial infections. Horizontal transmission has dominated the evolution of the TFSS, as demonstrated here by a lack of congruence between the tree topology inferred from components of the TFSS and the presumed bacterial species divergence pattern. A parsimony analysis suggests that conjugation represents the ancestral state and that the divergence from conjugation to secretion of effector molecules has occurred independently at multiple sites in the tree. The result shows that the nodes at which functional shifts have occurred coincide with those of horizontal gene transfers among distantly related bacteria. We suggest that it is the transfer between species that paved the way for the divergence of the TFSSs and discuss the general role of horizontal gene transfers for the evolution of novel gene functions.

Phylogenetic relationships and species differentiation of 39 Legionella species by sequence determination of the RNase P RNA gene rnpB.

The rnpB gene is ubiquitous in Bacteria, Archaea and Eucarya and encodes the RNA component of RNase P, an endoribonuclease P that consists of one RNA and one protein subunit (C5). In this study, partial rnpB genes were sequenced from 39 type strains and 16 additional strains of the genus Legionella. Models of the putative secondary structures of the RNase P RNA in the genus Legionella are proposed and possible interactions between RNase P RNA and C5 are discussed. The phylogenetic relationships within the genus Legionella were examined and rnpB sequences indicated six main clades that together comprised 27 of the 39 species examined. The phylogenetic relationships were further inferred by analysing combined datasets of sequences from the rnpB, mip, 16S rRNA and rpoB genes. It is concluded that rnpB is suitable for use in phylogenetic studies of closely related species and that it exhibits the potential to discriminate between Legionella species.

LDDist: a Perl module for calculating LogDet pair-wise distances for protein and nucleotide sequences.

LDDist is a Perl module implemented in C++ that allows the user to calculate LogDet pair-wise genetic distances for amino acid as well as nucleotide sequence data. It can handle site-to-site rate variation by treating a proportion of the sites as invariant and/or by assigning sites to different, presumably homogenous, rate categories. The rate-class assignments and invariant proportion can be set explicitly, or estimated by the program; the latter using either of two different capture-recapture methods. The assignment to rate categories in lieu of a phylogeny can be done using Shannon-Wiener index as a crude token for relative rate.

Phylogenetic relationships and genotyping of the genus Streptococcus by sequence determination of the RNase P RNA gene, rnpB.

The rnpB gene is universally present in bacterial species and encodes the RNA subunit of endoribonuclease P. In this study, rnpB was sequenced in 50 type strains and 29 additional strains of the genus Streptococcus. Putative secondary-structure models and possible interactions in RNase P RNA molecules are discussed. Phylogenetic relationships were studied and Bayesian, maximum-parsimony and minimum-evolution analyses supported six main clades that comprised 22 of the 50 species analysed. Phylogenetic inference was also studied for the 16S rRNA gene; it indicated a similar tree topology, but with weaker support values than for rnpB. Combined analysis of rnpB and 16S resulted in a phylogeny with significantly better support. Variability in the rnpB and 16S genes among all type strains, calculated as Shannon-Wiener information index values, was 0.45 for rnpB and 0.15 for 16S. Intraspecies proximity was assessed by principal coordinate analysis of rnpB for 32 strains of six closely related species (two clades) and showed species-specific clusters, but heterogeneity occurred in two species. It can be concluded that the rnpB gene is suitable for phylogenetic analysis of closely related taxa and has potential as a tool for species discrimination.

Bilaterian phylogeny based on analyses of a region of the sodium-potassium ATPase beta-subunit gene.

Molecular investigations of deep-level relationships within and among the animal phyla have been hampered by a lack of slowly evolving genes that are amenable to study by molecular systematists. To provide new data for use in deep-level metazoan phylogenetic studies, primers were developed to amplify a 1.3-kb region of the alpha subunit of the nuclear-encoded sodium-potassium ATPase gene from 31 bilaterians representing several phyla. Maximum parsimony, maximum likelihood, and Bayesian analyses of these sequences (combined with ATPase sequences for 23 taxa downloaded from GenBank) yield congruent trees that corroborate recent findings based on analyses of other data sets (e.g., the 18S ribosomal RNA gene). The ATPase-based trees support monophyly for several clades (including Lophotrochozoa, a form of Ecdysozoa, Vertebrata, Mollusca, Bivalvia, Gastropoda, Arachnida, Hexapoda, Coleoptera, and Diptera) but do not support monophyly for Deuterostomia, Arthropoda, or Nemertea. Parametric bootstrapping tests reject monophyly for Arthropoda and Nemertea but are unable to reject deuterostome monophyly. Overall, the sodium-potassium ATPase alpha-subunit gene appears to be useful for deep-level studies of metazoan phylogeny.