Genome-wide evidence for positive selection and recombination in Actinobacillus pleuropneumoniae

Metagenomic surveillance and comparative genomic analysis of Chlamydia psittaci in patients with pneumonia

Chlamydia psittaci, a strictly intracellular bacterium, is an underestimated etiologic agent leading to infections in a broad range of animals and mild illness or pneumonia in humans. In this study, the metagenomes of bronchoalveolar lavage fluids from the patients with pneumonia were sequenced and highly abundant C. psittaci was found. The target-enriched metagenomic reads were recruited to reconstruct draft genomes with more than 99% completeness. Two C. psittaci strains from novel sequence types were detected and these were closely related to the animal-borne isolates derived from the lineages of ST43 and ST28, indicating the zoonotic transmissions of C. psittaci would benefit its prevalence worldwide. Comparative genomic analysis combined with public isolate genomes revealed that the pan-genome of C. psittaci possessed a more stable gene repertoire than those of other extracellular bacteria, with ~90% of the genes per genome being conserved core genes. Furthermore, the evidence for significantly positive selection was identified in 20 virulence-associated gene products, particularly bacterial membrane-embedded proteins and type three secretion machines, which may play important roles in the pathogen-host interactions. This survey uncovered novel strains of C. psittaci causing pneumonia and the evolutionary analysis characterized prominent gene candidates involved in bacterial adaptation to immune pressures. The metagenomic approach is of significance to the surveillance of difficult-to-culture intracellular pathogens and the research into molecular epidemiology and evolutionary biology of C. psittaci.

Wolbachia in Black Spiny Whiteflies and Their New Parasitoid Wasp in Japan: Evidence of the Distinct Infection Status on Aleurocanthus camelliae Cryptic Species Complex

Wolbachia, an alphaproteobacterial reproductive parasite, can cause profound mitochondrial divergence in insects, which might eventually be a part of cryptic speciation. Aleurocanthus camelliae is a cryptic species complex consisting of several morphospecies and/or haplotypes that are genetically different but morphologically indistinctive. However, little is known about the Wolbachia infection status in these tea and Citrus pests. Thus, this study aimed to profile the diversity and phenotypic characteristics of Wolbachia natural infections in the A. camelliae cryptic species complex. A monophyletic strain of Wolbachia that infected the A. camelliae cryptic species complex (wAlec) with different patterns was discovered. Whiteflies that are morphologically identical to Aleurocanthus spiniferus (Aleurocanthus cf. A. spiniferus in Eurya japonica and A. spiniferus in Citrus) were grouped into uninfected populations, whereas the fixed infection was detected in A. camelliae B1 from Theaceae. The rapid evolution of wAlec was also found to occur through a high recombination event, which produced subgroups A and B in wAlec. It may also be associated with the non-cytoplasmic incompatibility (CI) phenotype of wAlec due to undetectable CI-related genes from phage WO (WOAlec). The current discovery of a novel cryptic species of A. camelliae led to a discussion about the oscillation hypothesis, which may provide insights on cryptic speciation, particularly on how specialization and host expansion have been recorded among these species. This study also identified a parasitoid wasp belonging to the genus Eretmocerus in A. camelliae, for the first time in Japan.

Characterization of Burkholderia cepacia Complex Core Genome and the Underlying Recombination and Positive Selection

Recombination and positive selection are two key factors that play a vital role in pathogenic microorganisms’ population adaptation and diversification. The Burkholderia cepacia complex (Bcc) represents bacterial species with high similarity, which can cause severe infections among cases suffering from the chronic granulomatous disorder and cystic fibrosis (CF). At present, no genome-wide study has been carried out focusing on investigating the core genome of Bcc associated with the two evolutionary forces. The general characteristics of the core genome of Bcc species remain scarce as well. In this study, we explored the core orthologous genes of 116 Bcc strains using comparative genomic analysis and studied the two adaptive evolutionary forces: recombination and positive selection. We estimated 1005 orthogroups consisting entirely of single copy genes. These single copy orthologous genes in some Cluster of Orthologous Groups (COG) categories showed significant differences in the comparison of several evolutionary properties, and the encoding proteins were relatively simple and compact. Our findings showed that 5.8% of the core orthologous genes strongly supported recombination; in the meantime, 1.1% supported positive selection. We found that genes involved in protein synthesis as well as material transport and metabolism are favored by selection pressure. More importantly, homologous recombination contributed more genetic variation to a large number of genes and largely maintained the genetic cohesion in Bcc. This high level of recombination between Bcc species blurs their taxonomic boundaries, which leads Bcc species to be difficult or impossible to distinguish phenotypically and genotypically.

Comparative Genomics of Actinobacillus pleuropneumoniae Serotype 8 Reveals the Importance of Prophages in the Genetic Variability of the Species

Actinobacillus pleuropneumoniae is the etiologic agent of porcine pleuropneumonia. Currently, there are 18 different serotypes; the serotype 8 is the most widely distributed in the United States, Canada, United Kingdom, and southeastern Brazil. In this study, genomes of seven A. pleuropneumoniae serotype 8 clinical isolates were compared to the other genomes of twelve serotypes. The analyses of serotype 8 genomes resulted in a set of 2352 protein-coding sequences. Of these sequences, 76.6% are present in all serotypes, 18.5% are shared with some serotypes, and 4.9% were differential. This differential portion was characterized as a series of hypothetical and regulatory protein sequences: mobile element sequence. Synteny analysis demonstrated possible events of gene recombination and acquisition by horizontal gene transfer (HGT) in this species. A total of 30 sequences related to prophages were identified in the genomes. These sequences represented 0.3 to 3.5% of the genome of the strains analyzed, and 16 of them contained complete prophages. Similarity analysis between complete prophage sequences evidenced a possible HGT with species belonging to the family Pasteurellaceae. Thus, mobile genetic elements, such as prophages, are important components of the differential portion of the A. pleuropneumoniae genome and demonstrate a central role in the evolution of the species. This study represents the first study done to understand the genome of A. pleuropneumoniae serotype 8.

Identification of positively selected genes in human pathogenic treponemes: Syphilis-, yaws-, and bejel-causing strains differ in sets of genes showing adaptive evolution

BACKGROUND

Pathogenic treponemes related to Treponema pallidum are both human (causing syphilis, yaws, bejel) and animal pathogens (infections of primates, venereal spirochetosis in rabbits). A set of 11 treponemal genome sequences including those of five Treponema pallidum ssp. pallidum (TPA) strains (Nichols, DAL-1, Mexico A, SS14, Chicago), four T. p. ssp. pertenue (TPE) strains (CDC-2, Gauthier, Samoa D, Fribourg-Blanc), one T. p. ssp. endemicum (TEN) strain (Bosnia A) and one strain (Cuniculi A) of Treponema paraluisleporidarum ecovar Cuniculus (TPeC) were tested for the presence of positively selected genes.

METHODOLOGY/PRINCIPAL FINDINGS

A total of 1068 orthologous genes annotated in all 11 genomes were tested for the presence of positively selected genes using both site and branch-site models with CODEML (PAML package). Subsequent analyses with sequences obtained from 62 treponemal draft genomes were used for the identification of positively selected amino acid positions. Synthetic biotinylated peptides were designed to cover positively selected protein regions and these peptides were tested for reactivity with the patient's syphilis sera. Altogether, 22 positively selected genes were identified in the TP genomes and TPA sets of positively selected genes differed from TPE genes. While genetic variability among TPA strains was predominantly present in a number of genetic loci, genetic variability within TPE and TEN strains was distributed more equally along the chromosome. Several syphilitic sera were shown to react with some peptides derived from the protein sequences evolving under positive selection.

CONCLUSIONS/SIGNIFICANCE

The syphilis-, yaws-, and bejel-causing strains differed relative to sets of positively selected genes. Most of the positively selected chromosomal loci were identified among the TPA treponemes. The local accumulation of genetic variability suggests that the diversification of TPA strains took place predominantly in a limited number of genomic regions compared to the more dispersed genetic diversity differentiating TPE and TEN strains. The identification of positively selected sites in tpr genes and genes encoding outer membrane proteins suggests their role during infection of human and animal hosts. The driving force for adaptive evolution at these loci thus appears to be the host immune response as supported by observed reactivity of syphilitic sera with some peptides derived from protein sequences showing adaptive evolution.

The TAM: A Translocation and Assembly Module of the β-Barrel Assembly Machinery in Bacterial Outer Membranes

Assembly of proteins into the outer membrane is an essential process in the cell biology of bacteria. The integration of β-barrel proteins into the outer membrane is mediated by a system referred to as the β-barrel assembly machinery (BAM) that includes two related proteins: BamA in the BAM complex and TamA in the TAM (translocation and assembly module). Here we review what is known about the TAM in terms of its function and the structural architecture of its two subunits, TamA and TamB. By linking the energy transduction possibilities in the inner membrane to TamA in the outer membrane, the TAM provides additional capability to the β-barrel assembly machinery. Conservation of the TAM across evolutionary boundaries, and the presence of hybrid BAM/TAM complexes in some bacterial lineages, adds insight to our growing understanding of how bacterial outer membranes are built.

Genome rearrangements and selection in multi-chromosome bacteria Burkholderia spp

BACKGROUND

The genus Burkholderia consists of species that occupy remarkably diverse ecological niches. Its best known members are important pathogens, B. mallei and B. pseudomallei, which cause glanders and melioidosis, respectively. Burkholderia genomes are unusual due to their multichromosomal organization, generally comprised of 2-3 chromosomes.

RESULTS

We performed integrated genomic analysis of 127 Burkholderia strains. The pan-genome is open with the saturation to be reached between 86,000 and 88,000 genes. The reconstructed rearrangements indicate a strong avoidance of intra-replichore inversions that is likely caused by selection against the transfer of large groups of genes between the leading and the lagging strands. Translocated genes also tend to retain their position in the leading or the lagging strand, and this selection is stronger for large syntenies. Integrated reconstruction of chromosome rearrangements in the context of strains phylogeny reveals parallel rearrangements that may indicate inversion-based phase variation and integration of new genomic islands. In particular, we detected parallel inversions in the second chromosomes of B. pseudomallei with breakpoints formed by genes encoding membrane components of multidrug resistance complex, that may be linked to a phase variation mechanism. Two genomic islands, spreading horizontally between chromosomes, were detected in the B. cepacia group.

CONCLUSIONS

This study demonstrates the power of integrated analysis of pan-genomes, chromosome rearrangements, and selection regimes. Non-random inversion patterns indicate selective pressure, inversions are particularly frequent in a recent pathogen B. mallei, and, together with periods of positive selection at other branches, may indicate adaptation to new niches. One such adaptation could be a possible phase variation mechanism in B. pseudomallei.

Polymorphism analysis of the apxIA gene of Actinobacillus pleuropneumoniae serovar 5 isolated in swine herds from Brazil

The bacterium Actinobacillus pleuropneumoniae is the etiological agent of Contagious Porcine Pleuropneumonia, a disease responsible for economic losses in the swine industry worldwide. A. pleuropneumoniae is capable of producing proteinaceous exotoxins responsible for inducing hemorrhagic lesions, one of which is ApxI. Few studies have conducted an in-depth evaluation of polymorphisms of the nucleotides that make up the ApxI toxin gene. Here we analyze the polymorphisms of the apxIA gene region of A. pleuropneumoniae serovar 5 isolated from swine in different regions in Brazil and report the results of molecular sequencing and phylogenetic analysis. Analysis of the apxIA gene in 60 isolates revealed the presence of genetic diversity and variability. The polymorphisms in the nucleotide sequences determined the grouping of the Brazilian sequences and five more sequences from the GenBank database into 14 different haplotypes, which formed three main groups and revealed the presence of mutations in the nucleotide sequences. The estimation of selection pressures suggests the occurrence of genetic variations by positive selective pressure on A. pleuropneumoniae in large groups of animals in relatively small spaces. These conditions presumably favor the horizontal dissemination of apxIA gene mutations within bacterial populations with host reservoirs. As a result, the same serovar can demonstrate different antigenic capacities due to mutations in the apxIA gene. These alterations in sequences of the apxIA gene could occur in other areas of countries with intense swine production, which could lead to differences in the pathogenicity and immunogenicity of each serovar and have implications for the clinical status or diagnosis of A. pleuropneumoniae.

Genes under positive selection in the core genome of pathogenic Bacillus cereus group members

In this comparative genomics study our aim was to unravel genes under positive selection in the core genome of the Bacillus cereus group. Indeed, the members of this group share close genetic relationships but display a rather large phenotypic and ecological diversity, providing a unique opportunity for studying how genomic changes reflect ecological adaptation during the divergence of a bacterial group. For this purpose, we screened ten completely sequenced genomes of four pathogenic Bacillus species, finding that 254 out of 3093 genes have codon sites with d_N/d_S (ω) values above one. These results remained unchanged after having disentangled the confounding effects of recombination and selection signature in a Bayesian framework. The presumably adaptive nucleotide polymorphisms are distributed over a wide range of biological functions, such as antibiotic resistance, DNA repair, nutrient uptake, metabolism, cell wall assembly and spore structure. Our results indicate that adaptation to animal hosts, whether as pathogens, saprophytes or symbionts, is the major driving force in the evolution of the Bacillus cereus group. Future work should seek to understand the evolutionary dynamics of both core and accessory genes in an integrative framework to ultimately unravel the key networks involved in host adaptation.

Genome-Wide Analyses Reveal Genes Subject to Positive Selection in Pasteurella multocida

Pasteurella multocida, a Gram-negative opportunistic pathogen, has led to a broad range of diseases in mammals and birds, including fowl cholera in poultry, pneumonia and atrophic rhinitis in swine and rabbit, hemorrhagic septicemia in cattle, and bite infections in humans. In order to better interpret the genetic diversity and adaptation evolution of this pathogen, seven genomes of P. multocida strains isolated from fowls, rabbit and pigs were determined by using high-throughput sequencing approach. Together with publicly available P. multocida genomes, evolutionary features were systematically analyzed in this study. Clustering of 70,565 protein-coding genes showed that the pangenome of 33 P. multocida strains was composed of 1,602 core genes, 1,364 dispensable genes, and 1,070 strain-specific genes. Of these, we identified a full spectrum of genes related to virulence factors and revealed genetic diversity of these potential virulence markers across P. multocida strains, e.g., bcbAB, fcbC, lipA, bexDCA, ctrCD, lgtA, lgtC, lic2A involved in biogenesis of surface polysaccharides, hsf encoding autotransporter adhesin, and fhaB encoding filamentous haemagglutinin. Furthermore, based on genome-wide positive selection scanning, a total of 35 genes were subject to strong selection pressure. Extensive analyses of protein subcellular location indicated that membrane-associated genes were highly abundant among all positively selected genes. The detected amino acid sites undergoing adaptive selection were preferably located in extracellular space, perhaps associated with bacterial evasion of host immune responses. Our findings shed more light on conservation and distribution of virulence-associated genes across P. multocida strains. Meanwhile, this study provides a genetic context for future researches on the mechanism of adaptive evolution in P. multocida.

Genome-Wide Survey of Genes Under Positive Selection in Avian Pathogenic Escherichia coli Strains

The ability to obtain bacterial genomes from the same host has allowed for comparative studies that help in the understanding of the molecular evolution of specific pathotypes. Avian pathogenic Escherichia coli (APEC) is a group of extraintestinal strains responsible for causing colibacillosis in birds. APEC is also suggested to possess a role as a zoonotic agent. Despite its importance, APEC pathogenesis still has several cryptic pathogenic processes that need to be better understood. In this work, a genome-wide survey of eight APEC strains for genes with evidence of recombination revealed that ∼14% of the homologous groups evaluated present signs of recombination. Enrichment analyses revealed that nine Gene Ontology (GO) terms were significantly more represented in recombinant genes. Among these GO terms, several were noted to be ATP-related categories. The search for positive selection in these APEC genomes revealed 32 groups of homologous genes with evidence of positive selection. Among these groups, we found several related to cell metabolism, as well as several uncharacterized genes, beyond the well-known virulence factors ompC, lamB, waaW, waaL, and fliC. A GO term enrichment test showed a prevalence of terms related to bacterial cell contact with the external environment (e.g., viral entry into host cell, detection of virus, pore complex, bacterial-type flagellum filament C, and porin activity). Finally, the genes with evidence of positive selection were retrieved from genomes of non-APEC strains and tested as were done for APEC strains. The result revealed that none of the groups of genes presented evidence of positive selection, confirming that the analysis was effective in inferring positive selection for APEC and not for E. coli in general, which means that the study of the genes with evidence of positive selection identified in this study can contribute for the better understanding of APEC pathogenesis processes.

Genomic composition and dynamics among Methanomicrobiales predict adaptation to contrasting environments

Members of the order Methanomicrobiales are abundant, and sometimes dominant, hydrogenotrophic (H₂-CO₂ utilizing) methanoarchaea in a broad range of anoxic habitats. Despite their key roles in greenhouse gas emissions and waste conversion to methane, little is known about the physiological and genomic bases for their widespread distribution and abundance. In this study, we compared the genomes of nine diverse Methanomicrobiales strains, examined their pangenomes, reconstructed gene flow and identified genes putatively mediating their success across different habitats. Most strains slowly increased gene content whereas one, Methanocorpusculum labreanum, evidenced genome downsizing. Peat-dwelling Methanomicrobiales showed adaptations centered on improved transport of scarce inorganic nutrients and likely use H⁺ rather than Na⁺ transmembrane chemiosmotic gradients during energy conservation. In contrast, other Methanomicrobiales show the potential to concurrently use Na⁺ and H⁺ chemiosmotic gradients. Analyses also revealed that the Methanomicrobiales lack a canonical electron bifurcation system (MvhABGD) known to produce low potential electrons in other orders of hydrogenotrophic methanogens. Additional putative differences in anabolic metabolism suggest that the dynamics of interspecies electron transfer from Methanomicrobiales syntrophic partners can also differ considerably. Altogether, these findings suggest profound differences in electron trafficking in the Methanomicrobiales compared with other hydrogenotrophs, and warrant further functional evaluations.

Gene flow, recombination, and positive selection in Stenotrophomonas maltophilia: mechanisms underlying the diversity of the widespread opportunistic pathogen

Stenotrophomonas maltophilia is a global multidrug-resistant human opportunistic pathogen in clinical environments. Stenotrophomonas maltophilia is also ubiquitous in aqueous environments, soil, and plants. Various molecular typing methods have revealed that S. maltophilia exhibits high levels of phenotypic and genotypic diversity. However, information regarding the genomic diversity within S. maltophilia and the corresponding genetic mechanisms resulting in said diversity remain scarce. The genome sequences of 17 S. maltophilia strains were selected to investigate the mechanisms contributing to genetic diversity at the genome level. The core and large pan-genomes of the species were first estimated, resulting in a large, open pan-genome. A species phylogeny was also reconstructed based on 344 orthologous genes with one copy per genome, and the contribution of four evolutionary mechanisms to the species genome diversity was quantified: 15%-35% of the genes showed evidence for recombination, 0%-25% of the genes in one genome were likely gained, 0%-44% of the genes in some genomes were likely lost, and less than 0.3% of the genes in a genome were under positive selection pressures. We observed that, among the four main mechanisms, homologous recombination plays a key role in maintaining diversity in S. maltophilia. In this study, we provide an overview of evolution in S. maltophilia to provide a better understanding of its evolutionary dynamics and its relationship with genome diversity.

POTION: an end-to-end pipeline for positive Darwinian selection detection in genome-scale data through phylogenetic comparison of protein-coding genes

Background

Detection of genes evolving under positive Darwinian evolution in genome-scale data is nowadays a prevailing strategy in comparative genomics studies to identify genes potentially involved in adaptation processes. Despite the large number of studies aiming to detect and contextualize such gene sets, there is virtually no software available to perform this task in a general, automatic, large-scale and reliable manner. This certainly occurs due to the computational challenges involved in this task, such as the appropriate modeling of data under analysis, the computation time to perform several of the required steps when dealing with genome-scale data and the highly error-prone nature of the sequence and alignment data structures needed for genome-wide positive selection detection.

Results

We present POTION, an open source, modular and end-to-end software for genome-scale detection of positive Darwinian selection in groups of homologous coding sequences. Our software represents a key step towards genome-scale, automated detection of positive selection, from predicted coding sequences and their homology relationships to high-quality groups of positively selected genes. POTION reduces false positives through several sophisticated sequence and group filters based on numeric, phylogenetic, quality and conservation criteria to remove spurious data and through multiple hypothesis corrections, and considerably reduces computation time thanks to a parallelized design. Our software achieved a high classification performance when used to evaluate a curated dataset of Trypanosoma brucei paralogs previously surveyed for positive selection. When used to analyze predicted groups of homologous genes of 19 strains of Mycobacterium tuberculosis as a case study we demonstrated the filters implemented in POTION to remove sources of errors that commonly inflate errors in positive selection detection. A thorough literature review found no other software similar to POTION in terms of customization, scale and automation.

Conclusion

To the best of our knowledge, POTION is the first tool to allow users to construct and check hypotheses regarding the occurrence of site-based evidence of positive selection in non-curated, genome-scale data within a feasible time frame and with no human intervention after initial configuration. POTION is available at http://www.lmb.cnptia.embrapa.br/share/POTION/.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1765-0) contains supplementary material, which is available to authorized users.

Genome-wide detection of selection and other evolutionary forces

As is well known, pathogenic microbes evolve rapidly to escape from the host immune system and antibiotics. Genetic variations among microbial populations occur frequently during the long-term pathogen-host evolutionary arms race, and individual mutation beneficial for the fitness can be fixed preferentially. Many recent comparative genomics studies have pointed out the importance of selective forces in the molecular evolution of bacterial pathogens. The public availability of large-scale next-generation sequencing data and many state-of-the-art statistical methods of molecular evolution enable us to scan genome-wide alignments for evidence of positive Darwinian selection, recombination, and other evolutionary forces operating on the coding regions. In this chapter, we describe an integrative analysis pipeline and its application to tracking featured evolutionary trajectories on the genome of an animal pathogen. The evolutionary analysis of the protein-coding part of the genomes will provide a wide spectrum of genetic variations that play potential roles in adaptive evolution of bacteria.

A genome-wide identification of genes undergoing recombination and positive selection in Neisseria

Currently, there is particular interest in the molecular mechanisms of adaptive evolution in bacteria. Neisseria is a genus of gram negative bacteria, and there has recently been considerable focus on its two human pathogenic species N. meningitidis and N. gonorrhoeae. Until now, no genome-wide studies have attempted to scan for the genes related to adaptive evolution. For this reason, we selected 18 Neisseria genomes (14 N. meningitidis, 3 N. gonorrhoeae and 1 commensal N. lactamics) to conduct a comparative genome analysis to obtain a comprehensive understanding of the roles of natural selection and homologous recombination throughout the history of adaptive evolution. Among the 1012 core orthologous genes, we identified 635 genes with recombination signals and 10 genes that showed significant evidence of positive selection. Further functional analyses revealed that no functional bias was found in the recombined genes. Positively selected genes are prone to DNA processing and iron uptake, which are essential for the fundamental life cycle. Overall, the results indicate that both recombination and positive selection play crucial roles in the adaptive evolution of Neisseria genomes. The positively selected genes and the corresponding amino acid sites provide us with valuable targets for further research into the detailed mechanisms of adaptive evolution in Neisseria.

High levels of multiple infections, recombination and horizontal transmission of Wolbachia in the Andricus mukaigawae (Hymenoptera; Cynipidae) communities

Wolbachia are maternally inherited endosymbiotic bacteria of arthropods and nematodes. In arthropods, they manipulate the reproduction of their hosts to facilitate their own spread in host populations, causing cytoplasmic incompatibility, parthenogenesis induction, feminization of genetic males and male-killing. In this study, we investigated Wolbachia infection and studied wsp (Wolbachia surface protein) sequences in three wasp species associated with the unisexual galls of A. mukaigawae with the aim of determining the transmission mode and the reason for multiple infections of Wolbachia. Frequency of Wolbachia infected populations for A. mukaigawae, Synergus japonicus (inquiline), and Torymus sp. (parasitoid) was 75%, 100%, and 100%, respectively. Multiple Wolbachia infections were detected in A. mukaigawae and S. japonicus, with 5 and 8 Wolbachia strains, respectively. The two host species shared 5 Wolbachia strains and were infected by identical strains in several locations, indicating horizontal transmission of Wolbachia. The transmission potentially takes place through gall tissues, which the larvae of both wasps feed on. Furthermore, three recombination events of Wolbachia were observed: the strains W₈, W₂ and W₆ apparently have derived from W₃ and W_5a, W₆ and W₇, W₄ and W₉, respectively. W₈ and W₂ and their respective parental strains were detected in S. japonicus. W₆ was detected with only one parent (W₄) in S. japonicus; W₉ was detected in Torymus sp., suggesting horizontal transmission between hosts and parasitoids. In conclusion, our research supports earlier studies that horizontal transmission of Wolbachia, a symbiont of the Rickettsiales order, may be plant-mediated or take place between hosts and parasitoids. Our research provides novel molecular evidence for multiple recombination events of Wolbachia in gall wasp communities. We suggest that genomic recombination and potential plant-mediated horizontal transmission may be attributable to the high levels of multiple Wolbachia infections observed in A. mukaigawae and S. japonicus.

Genome evolution and phylogenomic analysis of Candidatus Kinetoplastibacterium, the betaproteobacterial endosymbionts of Strigomonas and Angomonas

It has been long known that insect-infecting trypanosomatid flagellates from the genera Angomonas and Strigomonas harbor bacterial endosymbionts (Candidatus Kinetoplastibacterium or TPE [trypanosomatid proteobacterial endosymbiont]) that supplement the host metabolism. Based on previous analyses of other bacterial endosymbiont genomes from other lineages, a stereotypical path of genome evolution in such bacteria over the duration of their association with the eukaryotic host has been characterized. In this work, we sequence and analyze the genomes of five TPEs, perform their metabolic reconstruction, do an extensive phylogenomic analyses with all available Betaproteobacteria, and compare the TPEs with their nearest betaproteobacterial relatives. We also identify a number of housekeeping and central metabolism genes that seem to have undergone positive selection. Our genome structure analyses show total synteny among the five TPEs despite millions of years of divergence, and that this lineage follows the common path of genome evolution observed in other endosymbionts of diverse ancestries. As previously suggested by cell biology and biochemistry experiments, Ca. Kinetoplastibacterium spp. preferentially maintain those genes necessary for the biosynthesis of compounds needed by their hosts. We have also shown that metabolic and informational genes related to the cooperation with the host are overrepresented amongst genes shown to be under positive selection. Finally, our phylogenomic analysis shows that, while being in the Alcaligenaceae family of Betaproteobacteria, the closest relatives of these endosymbionts are not in the genus Bordetella as previously reported, but more likely in the Taylorella genus.

ODoSE: a webserver for genome-wide calculation of adaptive divergence in prokaryotes

Quantifying patterns of adaptive divergence between taxa is a major goal in the comparative and evolutionary study of prokaryote genomes. When applied appropriately, the McDonald-Kreitman (MK) test is a powerful test of selection based on the relative frequency of non-synonymous and synonymous substitutions between species compared to non-synonymous and synonymous polymorphisms within species. The webserver ODoSE (Ortholog Direction of Selection Engine) allows the calculation of a novel extension of the MK test, the Direction of Selection (DoS) statistic, as well as the calculation of a weighted-average Neutrality Index (NI) statistic for the entire core genome, allowing for systematic analysis of the evolutionary forces shaping core genome divergence in prokaryotes. ODoSE is hosted in a Galaxy environment, which makes it easy to use and amenable to customization and is freely available at www.odose.nl.

Insight into the evolution of the histidine triad protein (HTP) family in Streptococcus

The Histidine Triad Proteins (HTPs), also known as Pht proteins in Streptococcus pneumoniae, constitute a family of surface-exposed proteins that exist in many pathogenic streptococcal species. Although many studies have revealed the importance of HTPs in streptococcal physiology and pathogenicity, little is known about their origin and evolution. In this study, after identifying all htp homologs from 105 streptococcal genomes representing 38 different species/subspecies, we analyzed their domain structures, positions in genome, and most importantly, their evolutionary histories. By further projecting this information onto the streptococcal phylogeny, we made several major findings. First, htp genes originated earlier than the Streptococcus genus and gene-loss events have occurred among three streptococcal groups, resulting in the absence of the htp gene in the Bovis, Mutans and Salivarius groups. Second, the copy number of htp genes in other groups of Streptococcus is variable, ranging from one to four functional copies. Third, both phylogenetic evidence and domain structure analyses support the division of two htp subfamilies, designated as htp I and htp II. Although present mainly in the pyogenic group and in Streptococcus suis, htp II members are distinct from htp I due to the presence of an additional leucine-rich-repeat domain at the C-terminus. Finally, htp genes exhibit a faster nucleotide substitution rate than do housekeeping genes. Specifically, the regions outside the HTP domains are under strong positive selection. This distinct evolutionary pattern likely helped Streptococcus to easily escape from recognition by host immunity.

Protein domain repetition is enriched in Streptococcal cell-surface proteins

Tandem repetition of domain in protein sequence occurs in all three domains of life. It creates protein diversity and adds functional complexity in organisms. In this work, we analyzed 52 streptococcal genomes and found 3748 proteins contained domain repeats. Proteins not harboring domain repeats are significantly enriched in cytoplasm, whereas proteins with domain repeats are significantly enriched in cytoplasmic membrane, cell wall and extracellular locations. Domain repetition occurs most frequently in S. pneumoniae and least in S. thermophilus and S. pyogenes. DUF1542 is the highest repeated domain in a single protein, followed by Rib, CW_binding_1, G5 and HemolysinCabind. 3D structures of 24 repeat-containing proteins were predicted to investigate the structural and functional effect of domain repetition. Several repeat-containing streptococcal cell surface proteins are known to be virulence-associated. Surface-associated tandem domain-containing proteins without experimental functional characterization may be potentially involved in the pathogenesis of streptococci and deserve further investigation.

Genome-wide survey of mutual homologous recombination in a highly sexual bacterial species

The nature of a species remains a fundamental and controversial question. The era of genome/metagenome sequencing has intensified the debate in prokaryotes because of extensive horizontal gene transfer. In this study, we conducted a genome-wide survey of outcrossing homologous recombination in the highly sexual bacterial species Helicobacter pylori. We conducted multiple genome alignment and analyzed the entire data set of one-to-one orthologous genes for its global strains. We detected mosaic structures due to repeated recombination events and discordant phylogenies throughout the genomes of this species. Most of these genes including the "core" set of genes and horizontally transferred genes showed at least one recombination event. Taking into account the relationship between the nucleotide diversity and the minimum number of recombination events per nucleotide, we evaluated the recombination rate in every gene. The rate appears constant across the genome, but genes with a particularly high or low recombination rate were detected. Interestingly, genes with high recombination included those for DNA transformation and for basic cellular functions, such as biosynthesis and metabolism. Several highly divergent genes with a high recombination rate included those for host interaction, such as outer membrane proteins and lipopolysaccharide synthesis. These results provide a global picture of genome-wide distribution of outcrossing homologous recombination in a bacterial species for the first time, to our knowledge, and illustrate how a species can be shaped by mutual homologous recombination.