Virulence evolution of the human pathogen Neisseria meningitidis by recombination in the core and accessory genome

Background

Neisseria meningitidis is a naturally transformable, facultative pathogen colonizing the human nasopharynx. Here, we analyze on a genome-wide level the impact of recombination on gene-complement diversity and virulence evolution in N. meningitidis. We combined comparative genome hybridization using microarrays (mCGH) and multilocus sequence typing (MLST) of 29 meningococcal isolates with computational comparison of a subset of seven meningococcal genome sequences.

Principal Findings

We found that lateral gene transfer of minimal mobile elements as well as prophages are major forces shaping meningococcal population structure. Extensive gene content comparison revealed novel associations of virulence with genetic elements besides the recently discovered meningococcal disease associated (MDA) island. In particular, we identified an association of virulence with a recently described canonical genomic island termed IHT-E and a differential distribution of genes encoding RTX toxin- and two-partner secretion systems among hyperinvasive and non-hyperinvasive lineages. By computationally screening also the core genome for signs of recombination, we provided evidence that about 40% of the meningococcal core genes are affected by recombination primarily within metabolic genes as well as genes involved in DNA replication and repair. By comparison with the results of previous mCGH studies, our data indicated that genetic structuring as revealed by mCGH is stable over time and highly similar for isolates from different geographic origins.

Conclusions

Recombination comprising lateral transfer of entire genes as well as homologous intragenic recombination has a profound impact on meningococcal population structure and genome composition. Our data support the hypothesis that meningococcal virulence is polygenic in nature and that differences in metabolism might contribute to virulence.     

Inter-genomic displacement via lateral gene transfer of bacterial <Emphasis Type="Italic">trp</Emphasis> operons in an overall context of vertical genealogy

Background  

The growing conviction that lateral gene transfer plays a significant role in prokaryote genealogy opens up a need for comprehensive evaluations of gene-enzyme systems on a case-by-case basis. Genes of tryptophan biosynthesis are frequently organized as whole-pathway operons, an attribute that is expected to facilitate multi-gene transfer in a single step. We have asked whether events of lateral gene transfer are sufficient to have obscured our ability to track the vertical genealogy that underpins tryptophan biosynthesis.     

Screening synteny blocks in pairwise genome comparisons through integer programming

Background  

It is difficult to accurately interpret chromosomal correspondences such as true orthology and paralogy due to significant divergence of genomes from a common ancestor. Analyses are particularly problematic among lineages that have repeatedly experienced whole genome duplication (WGD) events. To compare multiple "subgenomes" derived from genome duplications, we need to relax the traditional requirements of "one-to-one" syntenic matchings of genomic regions in order to reflect "one-to-many" or more generally "many-to-many" matchings. However this relaxation may result in the identification of synteny blocks that are derived from ancient shared WGDs that are not of interest. For many downstream analyses, we need to eliminate weak, low scoring alignments from pairwise genome comparisons. Our goal is to objectively select subset of synteny blocks whose total scores are maximized while respecting the duplication history of the genomes in comparison. We call this "quota-based" screening of synteny blocks in order to appropriately fill a quota of syntenic relationships within one genome or between two genomes having WGD events.     

Score-based prediction of genomic islands in prokaryotic genomes using hidden Markov models

Background  

Horizontal gene transfer (HGT) is considered a strong evolutionary force shaping the content of microbial genomes in a substantial manner. It is the difference in speed enabling the rapid adaptation to changing environmental demands that distinguishes HGT from gene genesis, duplications or mutations. For a precise characterization, algorithms are needed that identify transfer events with high reliability. Frequently, the transferred pieces of DNA have a considerable length, comprise several genes and are called genomic islands (GIs) or more specifically pathogenicity or symbiotic islands.     

Global gene expression analysis reveals reduced abundance of putative microRNA targets in human prostate tumours

Background

The availability of hundreds of bacterial genomes allowed a comparative genomic study of the Type VI Secretion System (T6SS), recently discovered as being involved in pathogenesis. By combining comparative and phylogenetic approaches using more than 500 prokaryotic genomes, we characterized the global T6SS genetic structure in terms of conservation, evolution and genomic organization.

Results

This genome wide analysis allowed the identification of a set of 13 proteins constituting the T6SS protein core and a set of conserved accessory proteins. 176 T6SS loci (encompassing 92 different bacteria) were identified and their comparison revealed that T6SS-encoded genes have a specific conserved genetic organization. Phylogenetic reconstruction based on the core genes showed that lateral transfer of the T6SS is probably its major way of dissemination among pathogenic and non-pathogenic bacteria. Furthermore, the sequence analysis of the VgrG proteins, proposed to be exported in a T6SS-dependent way, confirmed that some C-terminal regions possess domains showing similarities with adhesins or proteins with enzymatic functions.

Conclusion

The core of T6SS is composed of 13 proteins, conserved in both pathogenic and non-pathogenic bacteria. Subclasses of T6SS differ in regulatory and accessory protein content suggesting that T6SS has evolved to adapt to various microenvironments and specialized functions. Based on these results, new functional hypotheses concerning the assembly and function of T6SS proteins are proposed.     

A genomic analysis of the archaeal system Ignicoccus hospitalis-Nanoarchaeum equitans

Background

The relationship between the hyperthermophiles Ignicoccus hospitalis and Nanoarchaeum equitans is the only known example of a specific association between two species of Archaea. Little is known about the mechanisms that enable this relationship.

Results

We sequenced the complete genome of I. hospitalis and found it to be the smallest among independent, free-living organisms. A comparative genomic reconstruction suggests that the I. hospitalis lineage has lost most of the genes associated with a heterotrophic metabolism that is characteristic of most of the Crenarchaeota. A streamlined genome is also suggested by a low frequency of paralogs and fragmentation of many operons. However, this process appears to be partially balanced by lateral gene transfer from archaeal and bacterial sources.

Conclusions

A combination of genomic and cellular features suggests highly efficient adaptation to the low energy yield of sulfur-hydrogen respiration and efficient inorganic carbon and nitrogen assimilation. Evidence of lateral gene exchange between N. equitans and I. hospitalis indicates that the relationship has impacted both genomes. This association is the simplest symbiotic system known to date and a unique model for studying mechanisms of interspecific relationships at the genomic and metabolic levels.     

Horizontal Gene Transfers in prokaryotes show differential preferences for metabolic and translational genes

Background  

Horizontal gene transfer (HGT) is an important process, which contributes in bacterial pathogenesis and drug resistance. A number of methods have been proposed for detection of horizontal gene transfer. One successful approach to the detection of HGT events is due to Novichkov et al. (J. Bacteriology 186, 6575–85), who rely on comparing phylogenetic distances within a gene family with genomic distances of the source organisms. Building on their approach, we introduce outlier detection in the correlation between those two sets of distances. This approach is designed to detect horizontal transfers of core set of genes present in many bacteria. The principle behind method allows detection of xenologous gene displacements as well as acquisition of novel genes.     

JContextExplorer: a tree-based approach to facilitate cross-species genomic context comparison

Background

Cross-species comparisons of gene neighborhoods (also called genomic contexts) in microbes may provide insight into determining functionally related or co-regulated sets of genes, suggest annotations of previously un-annotated genes, and help to identify horizontal gene transfer events across microbial species. Existing tools to investigate genomic contexts, however, lack features for dynamically comparing and exploring genomic regions from multiple species. As DNA sequencing technologies improve and the number of whole sequenced microbial genomes increases, a user-friendly genome context comparison platform designed for use by a broad range of users promises to satisfy a growing need in the biological community.

Results

Here we present JContextExplorer: a tool that organizes genomic contexts into branching diagrams. We implement several alternative context-comparison and tree rendering algorithms, and allow for easy transitioning between different clustering algorithms. To facilitate genomic context analysis, our tool implements GUI features, such as text search filtering, point-and-click interrogation of individual contexts, and genomic visualization via a multi-genome browser. We demonstrate a use case of our tool by attempting to resolve annotation ambiguities between two highly homologous yet functionally distinct genes in a set of 22 alpha and gamma proteobacteria.

Conclusions

JContextExplorer should enable a broad range of users to analyze and explore genomic contexts. The program has been tested on Windows, Mac, and Linux operating systems, and is implemented both as an executable JAR file and java WebStart. Program executables, source code, and documentation is available at http://www.bme.ucdavis.edu/facciotti/resources_data/software/.     

Variation suggestive of horizontal gene transfer at a lipopolysaccharide (<Emphasis Type="Italic">lps</Emphasis>) biosynthetic locus in <Emphasis Type="Italic">Xanthomonas oryzae</Emphasis> pv. <Emphasis Type="Italic">oryzae</Emphasis>, the bacterial leaf blight pathogen of rice

Background  

In animal pathogenic bacteria, horizontal gene transfer events (HGT) have been frequently observed in genomic regions that encode functions involved in biosynthesis of the outer membrane located lipopolysaccharide (LPS). As a result, different strains of the same pathogen can have substantially different lps biosynthetic gene clusters. Since LPS is highly antigenic, the variation at lps loci is attributed to be of advantage in evading the host immune system. Although LPS has been suggested as a potentiator of plant defense responses, interstrain variation at lps biosynthetic gene clusters has not been reported for any plant pathogenic bacterium.     

Comparison of eukaryotic phylogenetic profiling approaches using species tree aware methods

Background  

Phylogenetic profiling encompasses an important set of methodologies for in silico high throughput inference of functional relationships between genes. The simplest profiles represent the distribution of gene presence-absence in a set of species as a sequence of 0's and 1's, and it is assumed that functionally related genes will have more similar profiles. The methodology has been successfully used in numerous studies of prokaryotic genomes, although its application in eukaryotes appears problematic, with reported low accuracy due to the complex genomic organization within this domain of life. Recently some groups have proposed an alternative approach based on the correlation of homologous gene group sizes, taking into account all potentially informative genetic events leading to a change in group size, regardless of whether they result in a de novo group gain or total gene group loss.     

Genomic and evolutionary comparisons of diazotrophic and pathogenic bacteria of the order Rhizobiales

Background

Species belonging to the Rhizobiales are intriguing and extensively researched for including both bacteria with the ability to fix nitrogen when in symbiosis with leguminous plants and pathogenic bacteria to animals and plants. Similarities between the strategies adopted by pathogenic and symbiotic Rhizobiales have been described, as well as high variability related to events of horizontal gene transfer. Although it is well known that chromosomal rearrangements, mutations and horizontal gene transfer influence the dynamics of bacterial genomes, in Rhizobiales, the scenario that determine pathogenic or symbiotic lifestyle are not clear and there are very few studies of comparative genomic between these classes of prokaryotic microorganisms trying to delineate the evolutionary characterization of symbiosis and pathogenesis.

Results

Non-symbiotic nitrogen-fixing bacteria and bacteria involved in bioremediation closer to symbionts and pathogens in study may assist in the origin and ancestry genes and the gene flow occurring in Rhizobiales. The genomic comparisons of 19 species of Rhizobiales, including nitrogen-fixing, bioremediators and pathogens resulted in 33 common clusters to biological nitrogen fixation and pathogenesis, 15 clusters exclusive to all nitrogen-fixing bacteria and bacteria involved in bioremediation, 13 clusters found in only some nitrogen-fixing and bioremediation bacteria, 01 cluster exclusive to some symbionts, and 01 cluster found only in some pathogens analyzed. In BBH performed to all strains studied, 77 common genes were obtained, 17 of which were related to biological nitrogen fixation and pathogenesis. Phylogenetic reconstructions for Fix, Nif, Nod, Vir, and Trb showed possible horizontal gene transfer events, grouping species of different phenotypes.

Conclusions

The presence of symbiotic and virulence genes in both pathogens and symbionts does not seem to be the only determinant factor for lifestyle evolution in these microorganisms, although they may act in common stages of host infection. The phylogenetic analysis for many distinct operons involved in these processes emphasizes the relevance of horizontal gene transfer events in the symbiotic and pathogenic similarity.     

The obesity and inflammatory marker haptoglobin attracts monocytes via interaction with chemokine (C-C motif) receptor 2 (CCR2)

Background

Phylogeographic reconstruction of some bacterial populations is hindered by low diversity coupled with high levels of lateral gene transfer. A comparison of recombination levels and diversity at seven housekeeping genes for eleven bacterial species, most of which are commonly cited as having high levels of lateral gene transfer shows that the relative contributions of homologous recombination versus mutation for Burkholderia pseudomallei is over two times higher than for Streptococcus pneumoniae and is thus the highest value yet reported in bacteria. Despite the potential for homologous recombination to increase diversity, B. pseudomallei exhibits a relative lack of diversity at these loci. In these situations, whole genome genotyping of orthologous shared single nucleotide polymorphism loci, discovered using next generation sequencing technologies, can provide very large data sets capable of estimating core phylogenetic relationships. We compared and searched 43 whole genome sequences of B. pseudomallei and its closest relatives for single nucleotide polymorphisms in orthologous shared regions to use in phylogenetic reconstruction.

Results

Bayesian phylogenetic analyses of >14,000 single nucleotide polymorphisms yielded completely resolved trees for these 43 strains with high levels of statistical support. These results enable a better understanding of a separate analysis of population differentiation among >1,700 B. pseudomallei isolates as defined by sequence data from seven housekeeping genes. We analyzed this larger data set for population structure and allele sharing that can be attributed to lateral gene transfer. Our results suggest that despite an almost panmictic population, we can detect two distinct populations of B. pseudomallei that conform to biogeographic patterns found in many plant and animal species. That is, separation along Wallace's Line, a biogeographic boundary between Southeast Asia and Australia.

Conclusion

We describe an Australian origin for B. pseudomallei, characterized by a single introduction event into Southeast Asia during a recent glacial period, and variable levels of lateral gene transfer within populations. These patterns provide insights into mechanisms of genetic diversification in B. pseudomallei and its closest relatives, and provide a framework for integrating the traditionally separate fields of population genetics and phylogenetics for other bacterial species with high levels of lateral gene transfer.     

ACID: annotation of cassette and integron data

Background  

Although integrons and their associated gene cassettes are present in ~10% of bacteria and can represent up to 3% of the genome in which they are found, very few have been properly identified and annotated in public databases. These genetic elements have been overlooked in comparison to other vectors that facilitate lateral gene transfer between microorganisms.     

Codivergence of mycoviruses with their hosts

Background

The associations between pathogens and their hosts are complex and can result from any combination of evolutionary events such as codivergence, switching, and duplication of the pathogen. Mycoviruses are RNA viruses which infect fungi and for which natural vectors are so far unknown. Thus, lateral transfer might be improbable and codivergence their dominant mode of evolution. Accordingly, mycoviruses are a suitable target for statistical tests of virus-host codivergence, but inference of mycovirus phylogenies might be difficult because of low sequence similarity even within families.

Methodology

We analyzed here the evolutionary dynamics of all mycovirus families by comparing virus and host phylogenies. Additionally, we assessed the sensitivity of the co-phylogenetic tests to the settings for inferring virus trees from their genome sequences and approximate, taxonomy-based host trees.

Conclusions

While sequence alignment filtering modes affected branch support, the overall results of the co-phylogenetic tests were significantly influenced only by the number of viruses sampled per family. The trees of the two largest families, Partitiviridae and Totiviridae, were significantly more similar to those of their hosts than expected by chance, and most individual host-virus links had a significant positive impact on the global fit, indicating that codivergence is the dominant mode of virus diversification. However, in this regard mycoviruses did not differ from closely related viruses sampled from non-fungus hosts. The remaining virus families were either dominated by other evolutionary modes or lacked an apparent overall pattern. As this negative result might be caused by insufficient taxon sampling, the most parsimonious hypothesis still is that host-parasite evolution is basically the same in all mycovirus families. This is the first study of mycovirus-host codivergence, and the results shed light not only on how mycovirus biology affects their co-phylogenetic relationships, but also on their presumable host range itself.     

Evolution of glutamate dehydrogenase genes: evidence for lateral gene transfer within and between prokaryotes and eukaryotes

Background

Lateral gene transfer can introduce genes with novel functions into genomes or replace genes with functionally similar orthologs or paralogs. Here we present a study of the occurrence of the latter gene replacement phenomenon in the four gene families encoding different classes of glutamate dehydrogenase (GDH), to evaluate and compare the patterns and rates of lateral gene transfer (LGT) in prokaryotes and eukaryotes.

Results

We extend the taxon sampling of gdh genes with nine new eukaryotic sequences and examine the phylogenetic distribution pattern of the various GDH classes in combination with maximum likelihood phylogenetic analyses. The distribution pattern analyses indicate that LGT has played a significant role in the evolution of the four gdh gene families. Indeed, a number of gene transfer events are identified by phylogenetic analyses, including numerous prokaryotic intra-domain transfers, some prokaryotic inter-domain transfers and several inter-domain transfers between prokaryotes and microbial eukaryotes (protists).

Conclusion

LGT has apparently affected eukaryotes and prokaryotes to a similar extent within the gdh gene families. In the absence of indications that the evolution of the gdh gene families is radically different from other families, these results suggest that gene transfer might be an important evolutionary mechanism in microbial eukaryote genome evolution.

     

Different patterns of evolution for duplicated DNA repair genes in bacteria of the Xanthomonadales group

Background

DNA repair genes encode proteins that protect organisms against genetic damage generated by environmental agents and by-products of cell metabolism. The importance of these genes in life maintenance is supported by their high conservation, and the presence of duplications of such genes may be easily traced, especially in prokaryotic genomes.

Results

The genome sequences of two Xanthomonas species were used as the basis for phylogenetic analyses of genes related to DNA repair that were found duplicated. Although 16S rRNA phylogenetic analyses confirm their classification at the basis of the gamma proteobacteria subdivision, differences were found in the origin of the various genes investigated. Except for lexA, detected as a recent duplication, most of the genes in more than one copy are represented by two highly divergent orthologs. Basically, one of such duplications is frequently positioned close to other gamma proteobacteria, but the second is often positioned close to unrelated bacteria. These orthologs may have occurred from old duplication events, followed by extensive gene loss, or were originated from lateral gene transfer (LGT), as is the case of the uvrD homolog.

Conclusions

Duplications of DNA repair related genes may result in redundancy and also improve the organisms' responses to environmental challenges. Most of such duplications, in Xanthomonas, seem to have arisen from old events and possibly enlarge both functional and evolutionary genome potentiality.

     

Chromosome-breakage genomic instability and chromothripsis in breast cancer

Background

Chromosomal breakage followed by faulty DNA repair leads to gene amplifications and deletions in cancers. However, the mere assessment of the extent of genomic changes, amplifications and deletions may reduce the complexity of genomic data observed by array comparative genomic hybridization (array CGH). We present here a novel approach to array CGH data analysis, which focuses on putative breakpoints responsible for rearrangements within the genome.

Results

We performed array comparative genomic hybridization in 29 primary tumors from high risk patients with breast cancer. The specimens were flow sorted according to ploidy to increase tumor cell purity prior to array CGH. We describe the number of chromosomal breaks as well as the patterns of breaks on individual chromosomes in each tumor. There were differences in chromosomal breakage patterns between the 3 clinical subtypes of breast cancers, although the highest density of breaks occurred at chromosome 17 in all subtypes, suggesting a particular proclivity of this chromosome for breaks. We also observed chromothripsis affecting various chromosomes in 41% of high risk breast cancers.

Conclusions

Our results provide a new insight into the genomic complexity of breast cancer. Genomic instability dependent on chromosomal breakage events is not stochastic, targeting some chromosomes clearly more than others. We report a much higher percentage of chromothripsis than described previously in other cancers and this suggests that massive genomic rearrangements occurring in a single catastrophic event may shape many breast cancer genomes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-579) contains supplementary material, which is available to authorized users.     

Chlamydiae has contributed at least 55 genes to Plantae with predominantly plastid functions

Background

The photosynthetic organelle (plastid) originated via primary endosymbiosis in which a phagotrophic protist captured and harnessed a cyanobacterium. The plastid was inherited by the common ancestor of the red, green (including land plants), and glaucophyte algae (together, the Plantae). Despite the critical importance of primary plastid endosymbiosis, its ancient derivation has left behind very few “footprints” of early key events in organelle genesis.

Methodology/Principal Findings

To gain insights into this process, we conducted an in-depth phylogenomic analysis of genomic data (nuclear proteins) from 17 Plantae species to identify genes of a surprising provenance in these taxa, Chlamydiae bacteria. Previous studies show that Chlamydiae contributed many genes (at least 21 in one study) to Plantae that primarily have plastid functions and were postulated to have played a fundamental role in organelle evolution. Using our comprehensive approach, we identify at least 55 Chlamydiae-derived genes in algae and plants, of which 67% (37/55) are putatively plastid targeted and at least 3 have mitochondrial functions. The remainder of the proteins does not contain a bioinformatically predicted organelle import signal although one has an N-terminal extension in comparison to the Chlamydiae homolog. Our data suggest that environmental Chlamydiae were significant contributors to early Plantae genomes that extend beyond plastid metabolism. The chlamydial gene distribution and protein tree topologies provide evidence for both endosymbiotic gene transfer and a horizontal gene transfer ratchet driven by recurrent endoparasitism as explanations for gene origin.

Conclusions/Significance

Our findings paint a more complex picture of gene origin than can easily be explained by endosymbiotic gene transfer from an organelle-like point source. These data significantly extend the genomic impact of Chlamydiae on Plantae and show that about one-half (30/55) of the transferred genes are most closely related to sequences emanating from the genome of the only environmental isolate that is currently available. This strain, Candidatus Protochlamydia amoebophila UWE25 is an endosymbiont of Acanthamoeba and likely represents the type of endoparasite that contributed the genes to Plantae.     

Ankle ligament reconstruction after wide resection of the osteosarcoma of the distal fibula: a case report

Background

Restoration of the lateral ankle after distal fibulectomy is a difficult reconstructive procedure. Many surgical techniques have been proposed. This report shows another fibular reconstructive option with promising outcome.

Case presentation

We report the case of a 30-year-old woman who presented with a solitary mass located in the lateral aspect of the ankle. The mass had grown rapidly for 2 months and caused increasing pain. Physical examination showed a 3.0 cm diameter tender, nonmobile hard mass in the lateral malleolus. Radiographs showed an osteolytic lesion involving the lateral cortex at the distal fibula. After incisional biopsy, pathologic examination found a well-differentiated intramedullary osteosarcoma. Neoadjuvant chemotherapy with doxorubicin was provided for 3 months prior to definitive surgical treatment. Magnetic resonance imaging showed persistent tumor in the biopsy site. After distal fibulectomy and wide resection, split tibialis posterior tendon transfer to the remaining peroneus brevis restored the stability of the ankle. The pain resolved within 3 months. The ankle was stable and no recurrence of the cancer was found at a 7 year follow-up.

Conclusion

Reconstruction following distal fibulectomy and surrounding soft tissue resection responds favorably to split tibialis posterior transfer to the remaining peroneus brevis suggesting that this technique can provide a good and functional outcome.