Reconciling the many faces of lateral gene transfer

Molecular and evolutionary basis of O-antigenic polysaccharide-driven phage sensitivity in environmental pseudomonads

IMPORTANCE

The application of plant-beneficial microorganisms to protect crop plants is a promising alternative to the usage of chemicals. However, biocontrol research often faces difficulties in implementing this approach due to the inconsistency of the bacterial inoculant to establish itself within the root microbiome. Beneficial bacterial inoculants can be decimated by the presence of their natural predators, notably bacteriophages (also called phages). Thus, it is important to gain knowledge regarding the mechanisms behind phage-bacteria interactions to overcome this challenge. Here, we evidence that the major long O-antigenic polysaccharide (O-PS, O-antigen) of the widely used model plant-beneficial bacterium Pseudomonas protegens CHA0 is the receptor of its natural predator, the phage ΦGP100. We examined the distribution of the gene cluster directing the synthesis of this O-PS and identified signatures of horizontal gene acquisitions. Altogether, our study highlights the importance of bacterial cell surface structure variation in the complex interplay between phages and their Pseudomonas hosts.

Dissecting the HGT network of carbon metabolic genes in soil-borne microbiota

The microbiota inhabiting soil plays a significant role in essential life-supporting element cycles. Here, we investigated the occurrence of horizontal gene transfer (HGT) and established the HGT network of carbon metabolic genes in 764 soil-borne microbiota genomes. Our study sheds light on the crucial role of HGT components in microbiological diversification that could have far-reaching implications in understanding how these microbial communities adapt to changing environments, ultimately impacting agricultural practices. In the overall HGT network of carbon metabolic genes in soil-borne microbiota, a total of 6,770 nodes and 3,812 edges are present. Among these nodes, phyla Proteobacteria, Actinobacteriota, Bacteroidota, and Firmicutes are predominant. Regarding specific classes, Actinobacteria, Gammaproteobacteria, Alphaproteobacteria, Bacteroidia, Actinomycetia, Betaproteobacteria, and Clostridia are dominant. The Kyoto Encyclopedia of Genes and Genomes (KEGG) functional assignments of glycosyltransferase (18.5%), glycolysis/gluconeogenesis (8.8%), carbohydrate-related transporter (7.9%), fatty acid biosynthesis (6.5%), benzoate degradation (3.1%) and butanoate metabolism (3.0%) are primarily identified. Glycosyltransferase involved in cell wall biosynthesis, glycosylation, and primary/secondary metabolism (with 363 HGT entries), ranks first overwhelmingly in the list of most frequently identified carbon metabolic HGT enzymes, followed by pimeloyl-ACP methyl ester carboxylesterase, alcohol dehydrogenase, and 3-oxoacyl-ACP reductase. Such HGT events mainly occur in the peripheral functions of the carbon metabolic pathway instead of the core section. The inter-microbe HGT genetic traits in soil-borne microbiota genetic sequences that we recognized, as well as their involvement in the metabolism and regulation processes of carbon organic, suggest a pervasive and substantial effect of HGT on the evolution of microbes.

Reconstructing horizontal gene flow network to understand prokaryotic evolution

Horizontal gene transfer (HGT) is a major source of phenotypic innovation and a mechanism of niche adaptation in prokaryotes. Quantification of HGT is critical to decipher its myriad roles in microbial evolution and adaptation. Advances in genome sequencing and bioinformatics have augmented our ability to understand the microbial world, particularly the direct or indirect influence of HGT on diverse life forms. Methods for detecting HGT can be classified into phylogenetic-based and parametric or composition-based approaches. Here, we exploited the complementary strengths of both the approaches to construct a high confidence horizontal gene flow network. Our network is unique in its ability to detect the transfer of native genes of a genome to genomes from other taxa, thus establishing donor and recipient organisms (taxa), rather than through a post hoc analysis as is the practice with several other approaches. The scale-free horizontal gene flow network presented here provides new insights into modes of transfer for the exchange of genetic information and also illuminates differential gene flow across phyla.

Reconstructing horizontal gene flow network to understand prokaryotic evolution

Horizontal gene transfer (HGT) is a major source of phenotypic innovation and a mechanism of niche adaptation in prokaryotes. Quantification of HGT is critical to decipher its myriad roles in microbial evolution and adaptation. Advances in genome sequencing and bioinformatics have augmented our ability to understand the microbial world, particularly the direct or indirect influence of HGT on diverse life forms. Methods for detecting HGT can be classified into phylogenetic-based and parametric or composition-based approaches. Here, we exploited the complementary strengths of both the approaches to construct a high confidence horizontal gene flow network. Our network is unique in its ability to detect the transfer of native genes of a genome to genomes from other taxa, thus establishing donor and recipient organisms (taxa), rather than through a post hoc analysis as is the practice with several other approaches. The scale-free horizontal gene flow network presented here provides new insights into modes of transfer for the exchange of genetic information and also illuminates differential gene flow across phyla.

The microbial ecology of Escherichia coli in the vertebrate gut

Escherichia coli has a rich history as biology's 'rock star', driving advances across many fields. In the wild, E. coli resides innocuously in the gut of humans and animals but is also a versatile pathogen commonly associated with intestinal and extraintestinal infections and antimicrobial resistance-including large foodborne outbreaks such as the one that swept across Europe in 2011, killing 54 individuals and causing approximately 4000 infections and 900 cases of haemolytic uraemic syndrome. Given that most E. coli are harmless gut colonizers, an important ecological question plaguing microbiologists is what makes E. coli an occasionally devastating pathogen? To address this question requires an enhanced understanding of the ecology of the organism as a commensal. Here, we review how our knowledge of the ecology and within-host diversity of this organism in the vertebrate gut has progressed in the 137 years since E. coli was first described. We also review current approaches to the study of within-host bacterial diversity. In closing, we discuss some of the outstanding questions yet to be addressed and prospects for future research.

Bioinformatic and Mechanistic Analysis of the Palmerolide PKS-NRPS Biosynthetic Pathway From the Microbiome of an Antarctic Ascidian

Complex interactions exist between microbiomes and their hosts. Increasingly, defensive metabolites that have been attributed to host biosynthetic capability are now being recognized as products of host-associated microbes. These unique metabolites often have bioactivity targets in human disease and can be purposed as pharmaceuticals. Polyketides are a complex family of natural products that often serve as defensive metabolites for competitive or pro-survival purposes for the producing organism, while demonstrating bioactivity in human diseases as cholesterol lowering agents, anti-infectives, and anti-tumor agents. Marine invertebrates and microbes are a rich source of polyketides. Palmerolide A, a polyketide isolated from the Antarctic ascidian Synoicum adareanum, is a vacuolar-ATPase inhibitor with potent bioactivity against melanoma cell lines. The biosynthetic gene clusters (BGCs) responsible for production of secondary metabolites are encoded in the genomes of the producers as discrete genomic elements. A candidate palmerolide BGC was identified from a S. adareanum microbiome-metagenome based on a high degree of congruence with a chemical structure-based retrobiosynthetic prediction. Protein family homology analysis, conserved domain searches, active site and motif identification were used to identify and propose the function of the ∼75 kbp trans-acyltransferase (AT) polyketide synthase-non-ribosomal synthase (PKS-NRPS) domains responsible for the stepwise synthesis of palmerolide A. Though PKS systems often act in a predictable co-linear sequence, this BGC includes multiple trans-acting enzymatic domains, a non-canonical condensation termination domain, a bacterial luciferase-like monooxygenase (LLM), and is found in multiple copies within the metagenome-assembled genome (MAG). Detailed inspection of the five highly similar pal BGC copies suggests the potential for biosynthesis of other members of the palmerolide chemical family. This is the first delineation of a biosynthetic gene cluster from an Antarctic microbial species, recently proposed as Candidatus Synoicihabitans palmerolidicus. These findings have relevance for fundamental knowledge of PKS combinatorial biosynthesis and could enhance drug development efforts of palmerolide A through heterologous gene expression.

A review of the taxonomy, genetics, and biology of the genus Escherichia and the type species Escherichia coli

Historically, bacteriologists have relied heavily on biochemical and structural phenotypes for bacterial taxonomic classification. However, advances in comparative genomics have led to greater insights into the remarkable genetic diversity within the microbial world, and even within well-accepted species such as Escherichia coli. The extraordinary genetic diversity in E. coli recapitulates the evolutionary radiation of this species in exploiting a wide range of niches (i.e., ecotypes), including the gastrointestinal system of diverse vertebrate hosts as well as non-host natural environments (soil, natural waters, wastewater), which drives the adaptation, natural selection, and evolution of intragenotypic conspecific specialism as a strategy for survival. Over the last few years, there has been increasing evidence that many E. coli strains are very host (or niche)-specific. While biochemical and phylogenetic evidence support the classification of E. coli as a distinct species, the vast genomic (diverse pan-genome and intragenotypic variability), phenotypic (e.g., metabolic pathways), and ecotypic (host-/niche-specificity) diversity, comparable to the diversity observed in known species complexes, suggest that E. coli is better represented as a complex. Herein we review the taxonomic classification of the genus Escherichia and discuss how phenotype, genotype, and ecotype recapitulate our understanding of the biology of this remarkable bacterium.

A Model to Investigate the Impact of Farm Practice on Antimicrobial Resistance in UK Dairy Farms

The ecological and human health impact of antibiotic use and the related antimicrobial resistance (AMR) in animal husbandry is poorly understood. In many countries, there has been considerable pressure to reduce overall antibiotic use in agriculture or to cease or minimise use of human critical antibiotics. However, a more nuanced approach would consider the differential impact of use of different antibiotic classes; for example, it is not known whether reduced use of bacteriostatic or bacteriolytic classes of antibiotics would be of greater value. We have developed an ordinary differential equation model to investigate the effects of farm practice on the spread and persistence of AMR in the dairy slurry tank environment. We model the chemical fate of bacteriolytic and bacteriostatic antibiotics within the slurry and their effect on a population of bacteria, which are capable of resistance to both types of antibiotic. Through our analysis, we find that changing the rate at which a slurry tank is emptied may delay the proliferation of multidrug-resistant bacteria by up to five years depending on conditions. This finding has implications for farming practice and the policies that influence waste management practices. We also find that, within our model, the development of multidrug resistance is particularly sensitive to the use of bacteriolytic antibiotics, rather than bacteriostatic antibiotics, and this may be cause for controlling the usage of bacteriolytic antibiotics in agriculture.

Young JP, Andersen SU. Symbiosis genes show a unique pattern of introgression and selection within a Rhizobium leguminosarum species complex

Rhizobia supply legumes with fixed nitrogen using a set of symbiosis genes. These can cross rhizobium species boundaries, but it is unclear how many other genes show similar mobility. Here, we investigate inter-species introgression using de novo assembly of 196 Rhizobium leguminosarum sv. trifolii genomes. The 196 strains constituted a five-species complex, and we calculated introgression scores based on gene-tree traversal to identify 171 genes that frequently cross species boundaries. Rather than relying on the gene order of a single reference strain, we clustered the introgressing genes into four blocks based on population structure-corrected linkage disequilibrium patterns. The two largest blocks comprised 125 genes and included the symbiosis genes, a smaller block contained 43 mainly chromosomal genes, and the last block consisted of three genes with variable genomic location. All introgression events were likely mediated by conjugation, but only the genes in the symbiosis linkage blocks displayed overrepresentation of distinct, high-frequency haplotypes. The three genes in the last block were core genes essential for symbiosis that had, in some cases, been mobilized on symbiosis plasmids. Inter-species introgression is thus not limited to symbiosis genes and plasmids, but other cases are infrequent and show distinct selection signatures.

Detecting horizontal gene transfer: a probabilistic approach

BACKGROUND

Horizontal gene transfer (HGT) is the event of a DNA sequence being transferred between species not by inheritance. HGT is a crucial factor in prokaryotic evolution and is a significant source for genomic novelty resulting in antibiotic resistance or the outbreak of virulent strains. Detection of HGT and the mechanisms responsible and enabling it, is hence of prime importance.Existing algorithms rely on a strong phylogenetic signal distinguishing the transferred sequence from its recipient genome. Closely related species pose an even greater challenge as most genes are very similar and therefore, the phylogenetic signal is weak anyhow. Notwithstanding, the importance of detecting HGT between such organisms is extremely high for the role of HGT in the emergence of new highly virulent strains.

RESULTS

In a recent work we devised a novel technique that relies on loss of synteny around a gene as a witness for HGT. We used a novel heuristic for synteny measurement, SI (Syntent Index), and the technique was tested on both simulated and real data and was found to provide a greater sensitivity than other HGT techniques. This synteny-based approach suffers low specificity, in particular more closely related species. Here we devise an adaptive approach to cope with this by varying the criteria according to species distance. The new approach is doubly adaptive as it also considers the lengths of the genes being transferred. In particular, we use Chernoff bound to decree HGT both in simulations and real bacterial genomes taken from EggNog database.

CONCLUSIONS

Here we show empirically that this approach is more conservative than the previous χ² based approach and provides a lower false positive rate, especially for closely related species and under wide range of genome parameters.

Comparative in silico analysis of ftsZ gene from different bacteria reveals the preference for core set of codons in coding sequence structuring and secondary structural elements determination

The deluge of sequence information in the recent times provide us with an excellent opportunity to compare organisms on a large genomic scale. In this study we have tried to decipher the variation in the gene organization and structuring of a vital bacterial gene called ftsZ which codes for an integral component of the bacterial cell division, the FtsZ protein. FtsZ is homologous to tubulin protein and has been found to be ubiquitous in eubacteria. FtsZ is showing increasing promise as a target for antibacterial drug discovery. Our study of ftsZ protein from 143 different bacterial species spanning a wider range of morphological and physiological type demonstrates that the ftsZ gene of about ninety three percent of the organisms show relatively biased codon usage profile and significant GC deviation from their genomic GC content. Comparative codon usage analysis of ftsZ and a core housekeeping gene rpoB demonstrated that codon usage pattern of ftsZ CDS is shaped by natural selection to a large extent and mimics that of a housekeeping gene. We have also detected a tendency among the different organisms to utilize a core set of codons in structuring the ftsZ coding sequence. We observed that the compositional frequency of the amino acid serine in the FtsZ protein appears to be a indicator of the bacterial lifestyle. Our meticulous analysis of the ftsZ gene linked with the corresponding FtsZ protein show that there is a bias towards the use of specific synonymous codons particularly in the helix and strand regions of the multi-domain FtsZ protein. Overall our findings suggest that in an indispensable and vital protein such as FtsZ, there is an inherent tendency to maintain form for optimized performance in spite of the extrinsic variability in coding features.

Alternative characterizations of Fitch's xenology relation

Horizontal gene transfer (HGT) is an important factor for the evolution of prokaryotes as well as eukaryotes. According to Walter M. Fitch, two genes are xenologs if they are separated by at least one HGT. This concept is formalized through Fitch relations, which are defined as binary relations that comprise all pairs (x, y) of genes x and y for which y has been horizontally transferred at least once since it diverged from the last common ancestor of x and y. This definition, in particular, preserves the directional character of the transfer. Fitch relations are characterized by a small set of forbidden induced subgraphs on three vertices and can be recognized in linear time. The mathematical characterization of Fitch relations is crucial to understand whether putative xenology relations are at least to some extent "biologically feasible". In this contribution, we provide two novel characterizations of Fitch relations. In particular, these results allow us directly to reconstruct gene trees (together with the location of the horizontal transfer events) that explain the underlying Fitch relation. As a biological side result, we can conclude that the phylogenetic signal to infer these gene trees is entirely contained in those pairs of genes x and y for which no directional transfer has been taken place in the common history of y and the last common ancestor of x and y. In other words, non-HGT events provide the essential information about the gene trees. In addition, we utilize the new characterizations to present an alternative, short and elegant proof of the characterization theorem established by Geiß et al. (J Math Bio 77(5), 2018).

Horizontal gene transfer in human-associated microorganisms inferred by phylogenetic reconstruction and reconciliation

Horizontal gene transfer (HGT) is widespread in the evolution of prokaryotes, especially those associated with the human body. Here, we implemented large-scale gene-species phylogenetic tree reconstructions and reconciliations to identify putative HGT-derived genes in the reference genomes of microbiota isolated from six major human body sites by the NIH Human Microbiome Project. Comparisons with a control group representing microbial genomes from diverse natural environments indicated that HGT activity increased significantly in the genomes of human microbiota, which is confirmatory of previous findings. Roughly, more than half of total genes in the genomes of human-associated microbiota were transferred (donated or received) by HGT. Up to 60% of the detected HGTs occurred either prior to the colonization of the human body or involved bacteria residing in different body sites. The latter could suggest 'genetic crosstalk' and movement of bacterial genes within the human body via hitherto poorly understood mechanisms. We also observed that HGT activity increased significantly among closely-related microorganisms and especially when they were united by physical proximity, suggesting that the 'phylogenetic effect' can significantly boost HGT activity. Finally, we identified several core and widespread genes least influenced by HGT that could become useful markers for building robust 'trees of life' and address several outstanding technical challenges to improve the phylogeny-based genome-wide HGT detection method for future applications.

Bipartite Network Analysis of Gene Sharings in the Microbial World

Extensive microbial gene flows affect how we understand virology, microbiology, medical sciences, genetic modification, and evolutionary biology. Phylogenies only provide a narrow view of these gene flows: plasmids and viruses, lacking core genes, cannot be attached to cellular life on phylogenetic trees. Yet viruses and plasmids have a major impact on cellular evolution, affecting both the gene content and the dynamics of microbial communities. Using bipartite graphs that connect up to 149,000 clusters of homologous genes with 8,217 related and unrelated genomes, we can in particular show patterns of gene sharing that do not map neatly with the organismal phylogeny. Homologous genes are recycled by lateral gene transfer, and multiple copies of homologous genes are carried by otherwise completely unrelated (and possibly nested) genomes, that is, viruses, plasmids and prokaryotes. When a homologous gene is present on at least one plasmid or virus and at least one chromosome, a process of "gene externalization," affected by a postprocessed selected functional bias, takes place, especially in Bacteria. Bipartite graphs give us a view of vertical and horizontal gene flow beyond classic taxonomy on a single very large, analytically tractable, graph that goes beyond the cellular Web of Life.

Potential of the Burkholderia cepacia Complex to Produce 4-Hydroxy-3-Methyl-2-Alkyquinolines

A few Burkholderia species, especially Burkholderia pseudomallei, Burkholderia thailandensis, Burkholderia ambifaria, and Burkholderia cepacia, are known to produce and release various 4-hydroxy-3-methyl-2-alkylquinolines (HMAQs), a family of molecules analogous to the 4-hydroxy-2-alkylquinolines [aka 2-n-alkyl-4(1H)-quinolones] of Pseudomonas aeruginosa, which include the Pseudomonas quinolone signal (PQS). However, while these exoproducts play several roles in P. aeruginosa virulence and survival, the available literature is very limited on their distribution and function in Burkholderia. In this perspective article, we studied the distribution of the hmqABCDEFG operon, which encodes the enzymes involved in the biosynthesis of HMAQs, in the Burkholderia cepacia complex (Bcc) group. Based on the available sequence data, about one third of Bcc species carry a homolog of the hmqABCDEFG, and not all sequenced strains in a given species possess this operon. Looking at the synteny of genes surrounding the hmqABCDEFG operon, we found that for some species, the operon seems to have been deleted or replaced by other genes. Finally, we review the literature on the possible function of HMAQs. Understanding the Hmq system may provide clues concerning their functions in Bcc.

Genomics insights into ecotype formation of ammonia-oxidizing archaea in the deep ocean

Various lineages of ammonia-oxidizing archaea (AOA) are present in deep waters, but the mechanisms that determine ecotype formation are obscure. We studied 18 high-quality genomes of the marine group I AOA lineages (alpha, gamma and delta) from the Mariana and Ogasawara trenches. The genomes of alpha AOA resembled each other, while those of gamma and delta lineages were more divergent and had even undergone insertion of some phage genes. The instability of the gamma and delta AOA genomes could be partially due to the loss of DNA polymerase B (polB) and methyladenine DNA glycosylase (tag) genes responsible for the repair of point mutations. The alpha AOA genomes harbour genes encoding a thrombospondin-like outer membrane structure that probably serves as a barrier to gene flow. Moreover, the gamma and alpha AOA lineages rely on vitamin B₁₂ -independent MetE and B₁₂ -dependent MetH, respectively, for methionine synthesis. The delta AOA genome contains genes involved in uptake of sugar and peptide perhaps for heterotrophic lifestyle. Our study provides insights into co-occurrence of cladogenesis and anagenesis in the formation of AOA ecotypes that perform differently in nitrogen and carbon cycling in dark oceans.

Genome-wide analysis of horizontally acquired genes in the genus Mycobacterium

Horizontal gene transfer (HGT) was attributed as a major driving force for the innovation and evolution of prokaryotic genomes. Previously, multiple research endeavors were undertaken to decipher HGT in different bacterial lineages. The genus Mycobacterium houses some of the most deadly human pathogens; however, the impact of HGT in Mycobacterium has never been addressed in a systematic way. Previous initiatives to explore the genomic imprints of HGTs in Mycobacterium were focused on few selected species, specifically among the members of Mycobacterium tuberculosis complex. Considering the recent availability of a large number of genomes, the current study was initiated to decipher the probable events of HGTs among 109 completely sequenced Mycobacterium species. Our comprehensive phylogenetic analysis with more than 9,000 families of Mycobacterium proteins allowed us to list several instances of gene transfers spread across the Mycobacterium phylogeny. Moreover, by examining the topology of gene phylogenies here, we identified the species most likely to donate and receive these genes and provided a detailed overview of the putative functions these genes may be involved in. Our study suggested that horizontally acquired foreign genes had played an enduring role in the evolution of Mycobacterium genomes and have contributed to their metabolic versatility and pathogenicity.

"Candidatus Nitrosotenuis aquarius," an Ammonia-Oxidizing Archaeon from a Freshwater Aquarium Biofilter

Ammonia is a metabolic waste product excreted by aquatic organisms that causes toxicity when it accumulates. Aquaria and aquaculture systems therefore use biological filters that promote the growth of nitrifiers to convert ammonia to nitrate. Ammonia-oxidizing bacteria (AOB) have been isolated from aquarium biofilters and are available as commercial supplements, but recent evidence suggests that ammonia-oxidizing archaea (AOA) are abundant in aquarium biofilters. In this study, we report the cultivation and closed genome sequence of the novel AOA representative "Candidatus Nitrosotenuis aquarius," which was enriched from a freshwater aquarium biofilter. "Ca Nitrosotenuis aquarius" oxidizes ammonia stoichiometrically to nitrite with a concomitant increase in thaumarchaeotal cells and a generation time of 34.9 h. "Ca Nitrosotenuis aquarius" has an optimal growth temperature of 33°C, tolerates up to 3 mM NH₄Cl, and grows optimally at 0.05% salinity. Transmission electron microscopy revealed that "Ca Nitrosotenuis aquarius" cells are rod shaped, with a diameter of ∼0.4 μm and length ranging from 0.6 to 3.6 μm. In addition, these cells possess surface layers (S-layers) and multiple proteinaceous appendages. Phylogenetically, "Ca Nitrosotenuis aquarius" belongs to the group I.1a Thaumarchaeota, clustering with environmental sequences from freshwater aquarium biofilters, aquaculture systems, and wastewater treatment plants. The complete 1.70-Mbp genome contains genes involved in ammonia oxidation, bicarbonate assimilation, flagellum synthesis, chemotaxis, S-layer production, defense, and protein glycosylation. Incubations with differential inhibitors indicate that "Ca Nitrosotenuis aquarius"-like AOA contribute to ammonia oxidation within the aquarium biofilter from which it originated.IMPORTANCE Nitrification is a critical process for preventing ammonia toxicity in engineered biofilter environments. This work describes the cultivation and complete genome sequence of a novel AOA representative enriched from a freshwater aquarium biofilter. In addition, despite the common belief in the aquarium industry that AOB mediate ammonia oxidation, the present study suggests an in situ role for "Ca Nitrosotenuis aquarius"-like AOA in freshwater aquarium biofilters.

A Novel Strategy for Detecting Recent Horizontal Gene Transfer and Its Application to Rhizobium Strains

Recent horizontal gene transfer (HGT) is crucial for enabling microbes to rapidly adapt to their novel environments without relying upon rare beneficial mutations that arise spontaneously. For several years now, computational approaches have been developed to detect HGT, but they typically lack the sensitivity and ability to detect recent HGT events. Here we introduce a novel strategy, named RecentHGT. The number of genes undergoing recent HGT between two bacterial genomes was estimated by a new algorithm derived from the expectation-maximization algorithm and is based on the theoretical sequence-similarity distribution of orthologous genes. We tested the proposed strategy by applying it to a set of 10 Rhizobium genomes, and detected several large-scale recent HGT events. We also found that our strategy was more sensitive than other available HGT detection methods. These HGT events were mainly mediated by symbiotic plasmids. Our new strategy can provide clear evidence of recent HGT events and thus it brings us closer to the goal of detecting these potentially adaptive evolution processes in rhizobia as well as pathogens.

Pan-Genomic Approaches in Lactobacillus reuteri as a Porcine Probiotic: Investigation of Host Adaptation and Antipathogenic Activity

After the introduction of a ban on the use of antibiotic growth promoters (AGPs) for livestock, reuterin-producing Lactobacillus reuteri is getting attention as an alternative to AGPs. In this study, we investigated genetic features of L. reuteri associated with host specificity and antipathogenic effect. We isolated 104 L. reuteri strains from porcine feces, and 16 strains, composed of eight strains exhibiting the higher antipathogenic effect (group HS) and eight strains exhibiting the lower effect (group LS), were selected for genomic comparison. We generated draft genomes of the 16 isolates and investigated their pan-genome together with the 26 National Center for Biotechnology Information-registered genomes. L. reuteri genomes organized six clades with multi-locus sequence analysis, and the clade IV includes the 16 isolates. First, we identified six L. reuteri clade IV-specific genes including three hypothetical protein-coding genes. The three annotated genes encode transposases and cell surface proteins, indicating that these genes are the result of adaptation to the host gastrointestinal epithelia and that these host-specific traits were acquired by horizontal gene transfer. We also identified differences between groups HS and LS in the pdu-cbi-cob-hem gene cluster, which is essential for reuterin and cobalamin synthesis, and six genes specific to group HS are revealed. While the strains of group HS possessed all genes of this cluster, LS strains have lost many genes of the cluster. This study provides a deeper understanding of the relationship between probiotic properties and genomic features of L. reuteri.

A Preliminary List of Horizontally Transferred Genes in Prokaryotes Determined by Tree Reconstruction and Reconciliation

Genome-wide global detection of genes involved in horizontal gene transfer (HGT) remains an active area of research in medical microbiology and evolutionary genomics. Utilizing the explicit evolutionary method of comparing topologies of a total of 154,805 orthologous gene trees against corresponding 16S rRNA “reference” trees, we previously detected a total of 660,894 candidate HGT events in 2,472 completely-sequenced prokaryotic genomes. Here, we report an HGT-index for each individual gene-reference tree pair reconciliation, representing the total number of detected HGT events on the gene tree divided by the total number of genomes (taxa) member of that tree. HGT-index is thus a simple measure indicating the sensitivity of prokaryotic genes to participate (or not participate) in HGT. Our preliminary list provides HGT-indices for a total of 69,365 genes (detected in >10 and <50% available prokaryotic genomes) that are involved in a wide range of biological processes such as metabolism, information, and bacterial response to environment. Identification of horizontally-derived genes is important to combat antibiotic resistance and is a step forward toward reconstructions of improved phylogenies describing the history of life. Our effort is thus expected to benefit ongoing research in the fields of clinical microbiology and evolutionary biology.

A novel alignment-free method for detection of lateral genetic transfer based on TF-IDF

Lateral genetic transfer (LGT) plays an important role in the evolution of microbes. Existing computational methods for detecting genomic regions of putative lateral origin scale poorly to large data. Here, we propose a novel method based on TF-IDF (Term Frequency-Inverse Document Frequency) statistics to detect not only regions of lateral origin, but also their origin and direction of transfer, in sets of hierarchically structured nucleotide or protein sequences. This approach is based on the frequency distributions of k-mers in the sequences. If a set of contiguous k-mers appears sufficiently more frequently in another phyletic group than in its own, we infer that they have been transferred from the first group to the second. We performed rigorous tests of TF-IDF using simulated and empirical datasets. With the simulated data, we tested our method under different parameter settings for sequence length, substitution rate between and within groups and post-LGT, deletion rate, length of transferred region and k size, and found that we can detect LGT events with high precision and recall. Our method performs better than an established method, ALFY, which has high recall but low precision. Our method is efficient, with runtime increasing approximately linearly with sequence length.

Horizontal Gene Transfer Contributes to the Evolution of Arthropod Herbivory

Within animals, evolutionary transition toward herbivory is severely limited by the hostile characteristics of plants. Arthropods have nonetheless counteracted many nutritional and defensive barriers imposed by plants and are currently considered as the most successful animal herbivores in terrestrial ecosystems. We gather a body of evidence showing that genomes of various plant feeding insects and mites possess genes whose presence can only be explained by horizontal gene transfer (HGT). HGT is the asexual transmission of genetic information between reproductively isolated species. Although HGT is known to have great adaptive significance in prokaryotes, its impact on eukaryotic evolution remains obscure. Here, we show that laterally transferred genes into arthropods underpin many adaptations to phytophagy, including efficient assimilation and detoxification of plant produced metabolites. Horizontally acquired genes and the traits they encode often functionally diversify within arthropod recipients, enabling the colonization of more host plant species and organs. We demonstrate that HGT can drive metazoan evolution by uncovering its prominent role in the adaptations of arthropods to exploit plants.

Bioprospecting of Plant Growth Promoting Bacilli and Related Genera Prevalent in Soils of Pristine Sacred Groves: Biochemical and Molecular Approach

Bacillus spp. and related genera native to soils of the pristine sacred groves from Meghalaya, India were characterized using biochemical and 16S rRNA gene analysis which revealed dominance of Bacillus, Paenibacillus, Lysinibacillus and Viridibacillus in the groves. Biochemical estimation was carried out for in vitro testing of plant growth promoting traits present in these isolates. PCR screening were performed for plant growth-promoting related genes involved in the biosynthesis of acid phosphatase (AcPho), indolepyruvate decarboxylase (ipdC), 1-aminocyclopropane-1-carboxylate deaminase (accd) and siderophore biosynthesis protein (asbA). 76% of the sacred grove isolates gave an amplified fragment for AcPho. Three of the isolates gave an amplified fragment for IpdC gene. Apart from 2 isolates, all the other isolates including the reference strains were positive for the amplification of the accd gene indicating their potential to produce ACC deaminase enzyme. 42% of the isolates gave an amplified fragment for asbA gene indicating the potential ability of these isolates to produce the catechol type siderophore, petrobactin. Overall findings indicated multiple PGP genetic traits present in these isolates which suggested that these isolates are capable of expressing multiple PGP traits. Phylogenetic and sequence analysis of accd and asbA genes from the isolates revealed that asbA genes from Paenibacillus taichungiensis SG3 and Paenibacillus tylopili SG24 indicated the occurrence of intergeneric horizontal transfer between Paenibacillus and Bacillus.

Horizontal Gene Transfer and the History of Life

Microbes acquire DNA from a variety of sources. The last decades, which have seen the development of genome sequencing, have revealed that horizontal gene transfer has been a major evolutionary force that has constantly reshaped genomes throughout evolution. However, because the history of life must ultimately be deduced from gene phylogenies, the lack of methods to account for horizontal gene transfer has thrown into confusion the very concept of the tree of life. As a result, many questions remain open, but emerging methodological developments promise to use information conveyed by horizontal gene transfer that remains unexploited today.

HGTree: database of horizontally transferred genes determined by tree reconciliation

The HGTree database provides putative genome-wide horizontal gene transfer (HGT) information for 2472 completely sequenced prokaryotic genomes. This task is accomplished by reconstructing approximate maximum likelihood phylogenetic trees for each orthologous gene and corresponding 16S rRNA reference species sets and then reconciling the two trees under parsimony framework. The tree reconciliation method is generally considered to be a reliable way to detect HGT events but its practical use has remained limited because the method is computationally intensive and conceptually challenging. In this regard, HGTree (http://hgtree.snu.ac.kr) represents a useful addition to the biological community and enables quick and easy retrieval of information for HGT-acquired genes to better understand microbial taxonomy and evolution. The database is freely available and can be easily scaled and updated to keep pace with the rapid rise in genomic information.

Massive horizontal gene transfer, strictly vertical inheritance and ancient duplications differentially shape the evolution of Bacillus cereus enterotoxin operons hbl, cytK and nhe

Background

Bacillus cereus sensu lato comprises eight closely related species including the human pathogens Bacillus anthracis and Bacillus cereus. Within B. cereus sensu lato, chromosomally and plasmid-encoded toxins exist. While plasmid-mediated horizontal gene transfer of the emetic toxin, anthrax and insecticidal toxins is known, evolution of enterotoxin genes within the group has not been studied.

Results

We report draft genome assemblies of 25 strains, a phylogenetic network of 142 strains based on ANI derived from genome sequences and a phylogeny based on whole-genome SNP analysis. The data clearly support subdivision of B. cereus sensu lato into seven phylogenetic groups. While group I, V and VII represent B. pseudomycoides, B. toyonensis and B. cytotoxicus, which are distinguishable at species level (ANI border ≥ 96 %), strains ascribed to the other five species do not match phylogenic groups. The chromosomal enterotoxin operons nheABC and hblCDAB are abundant within B. cereus both isolated from infections and from the environment. While the duplicated hbl variant hbl_a is present in 22 % of all strains investigated, duplication of nheABC is extremely rare (0.02 %) and appears to be phylogenetically unstable. Distribution of toxin genes was matched to a master tree based on seven concatenated housekeeping genes, which depicts species relationships in B. cereus sensu lato as accurately as whole-genome comparisons. Comparison to the phylogeny of enterotoxin genes uncovered ample evidence for horizontal transfer of hbl, cytK and plcR, as well as frequent deletion of both toxins and duplication of hbl. No evidence for nhe deletion was found and stable horizontal transfer of nhe is rare. Therefore, evolution of B. cereus enterotoxin operons is shaped unexpectedly different for yet unknown reasons.

Conclusions

Frequent exchange of the pathogenicity factors hbl, cytK and plcR in B. cereus sensu lato appears to be an important mechanism of B. cereus virulence evolution, including so-called probiotic or non-pathogenic species, which might have consequences for risk assessment procedures. In contrast, exclusively vertical inheritance of nhe was observed, and since nhe-negative strains appear to be extremely rare, we suggest that fitness loss may be associated with deletion or horizontal transfer of the nhe operon.

Electronic supplementary material

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

Analysis of three genomes within the thermophilic bacterial species Caldanaerobacter subterraneus with a focus on carbon monoxide dehydrogenase evolution and hydrolase diversity

Background

The Caldanaerobacter subterraneus species includes thermophilic fermentative bacteria able to grow on carbohydrates substrates with acetate and L-alanine as the main products. In this study, comprehensive analysis of three genomes of C. subterraneus subspecies was carried in order to identify genes encoding key metabolic enzymes and to document the genomic basis for the evolution of these organisms.

Methods

Average nucleotide identity and in silico DNA relatedness were estimated for the studied C. subterraneus genomes. Genome synteny was evaluated using R2CAT software. Protein conservation was analyzed using mGenome Subtractor. Horizontal gene transfer was predicted through the GOHTAM pipeline (using tetranucleotide composition) and phylogenetic analyses (by maximum likelihood). Hydrolases were identified through the MEROPS and CAZy platforms.

Results

The three genomes of C. subterraneus showed high similarity, although there are substantial differences in their gene composition and organization. Each subspecies possesses a gene cluster encoding a carbon monoxide dehydrogenase (CODH) and an energy converting hydrogenase (ECH). The CODH gene is associated with an operon that resembles the Escherichia coli hydrogenase hyc/hyf operons, a novel genetic context distinct from that found in archetypical hydrogenogenic carboxydotrophs. Apart from the CODH-associated hydrogenase, these bacteria also contain other hydrogenases, encoded by ech and hyd genes. An Mbx ferredoxin:NADP oxidoreductase homolog similar to that originally described in the archaeon Pyrococcus furiosus was uniquely encoded in the C. subterraneus subsp. yonseiensis genome. Compositional analysis demonstrated that some genes of the CODH-ECH and mbx operons present distinct sequence patterns in relation to the majority of the other genes of each genome. Phylogenetic reconstructions of the genes from these operons and those from the ech operon are incongruent to the species tree. Notably, the cooS gene of C. subterraneus subsp. pacificus and its homologs in C. subterraneus subsp. tengcongensis and C. subterraneus subsp. yonseiensis form distinct clades. The strains have diverse hydrolytic enzymes and they appear to be proteolytic and glycolytic. Divergent glycosidases from 14 families, among them amylases, chitinases, alpha-glucosidases, beta-glucosidases, and cellulases, were identified. Each of the three genomes also contains around 100 proteases from 50 subfamilies, as well about ten different esterases.

Conclusions

Genomic information suggests that multiple horizontal gene transfers conferred the adaptation of C. subterraneus subspecies to extreme niches throughout the carbon monoxide utilization and hydrogen production. The variety of hydrolases found in their genomes indicate the versatility of the species in obtaining energy and carbon from diverse substrates, therefore these organisms constitute a remarkable resource of enzymes with biotechnological potential.

Electronic supplementary material

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

Detecting Horizontal Gene Transfer between Closely Related Taxa

Horizontal gene transfer (HGT), the transfer of genetic material between organisms, is crucial for genetic innovation and the evolution of genome architecture. Existing HGT detection algorithms rely on a strong phylogenetic signal distinguishing the transferred sequence from ancestral (vertically derived) genes in its recipient genome. Detecting HGT between closely related species or strains is challenging, as the phylogenetic signal is usually weak and the nucleotide composition is normally nearly identical. Nevertheless, there is a great importance in detecting HGT between congeneric species or strains, especially in clinical microbiology, where understanding the emergence of new virulent and drug-resistant strains is crucial, and often time-sensitive.

We developed a novel, self-contained technique named Near HGT, based on the synteny index, to measure the divergence of a gene from its native genomic environment and used it to identify candidate HGT events between closely related strains. The method confirms candidate transferred genes based on the constant relative mutability (CRM). Using CRM, the algorithm assigns a confidence score based on “unusual” sequence divergence. A gene exhibiting exceptional deviations according to both synteny and mutability criteria, is considered a validated HGT product. We first employed the technique to a set of three E. coli strains and detected several highly probable horizontally acquired genes. We then compared the method to existing HGT detection tools using a larger strain data set.

When combined with additional approaches our new algorithm provides richer picture and brings us closer to the goal of detecting all newly acquired genes in a particular strain.

The transfer of genetic material between organisms, usually denoted as horizontal (or lateral) gene transfer (HGT or LGT), is a prime mechanism in microbial evolution and responsible for genetic innovation and the evolution of genome architecture. Detecting HGT between closely related species or strains is imperative as drug-resistant pathogenic strains most often acquire their virulence from closely related bacteria. The proposed method combines two evolutionary signals that were not employed in the past for this task. One is the synteny index (SI), measuring the loss of synteny in an organism, and the other is a novel concept—constant relative mutability (CRM), maintaining that genes preserve their relative evolution rate along linages (although the latter ones may each change).

We show both in simulation and real biological data that the method is sound and, in the cases examined, provides stronger sensitivity than existing methods. We therefore believe this novel approach represents a significant advance, for the first time enabling the detection of previously ignored HGT events that will bring us closer to the goal of detecting all newly acquired genes in a particular strain.

Availability: The method is publicly available at http://research.haifa.ac.il/~ssagi/software/nearHGT.zip

Transfer of energy pathway genes in microbial enhanced biological phosphorus removal communities

BACKGROUND

Lateral gene transfer (LGT) is an important evolutionary process in microbial evolution. In sewage treatment plants, LGT of antibiotic resistance and xenobiotic degradation-related proteins has been suggested, but the role of LGT outside these processes is unknown. Microbial communities involved in Enhanced Biological Phosphorus Removal (EBPR) have been used to treat wastewater in the last 50 years and may provide insights into adaptation to an engineered environment. We introduce two different types of analysis to identify LGT in EBPR sewage communities, based on identifying assembled sequences with more than one strong taxonomic match, and on unusual phylogenetic patterns. We applied these methods to investigate the role of LGT in six energy-related metabolic pathways.

RESULTS

The analyses identified overlapping but non-identical sets of transferred enzymes. All of these were homologous with sequences from known mobile genetic elements, and many were also in close proximity to transposases and integrases in the EBPR data set. The taxonomic method had higher sensitivity than the phylogenetic method, identifying more potential LGTs. Both analyses identified the putative transfer of five enzymes within an Australian community, two in a Danish community, and none in a US-derived culture.

CONCLUSIONS

Our methods were able to identify sequences with unusual phylogenetic or compositional properties as candidate LGT events. The association of these candidates with known mobile elements supports the hypothesis of transfer. The results of our analysis strongly suggest that LGT has influenced the development of functionally important energy-related pathways in EBPR systems, but transfers may be unique to each community due to different operating conditions or taxonomic composition.

Inferring horizontal gene transfer

Horizontal or Lateral Gene Transfer (HGT or LGT) is the transmission of portions of genomic DNA between organisms through a process decoupled from vertical inheritance. In the presence of HGT events, different fragments of the genome are the result of different evolutionary histories. This can therefore complicate the investigations of evolutionary relatedness of lineages and species. Also, as HGT can bring into genomes radically different genotypes from distant lineages, or even new genes bearing new functions, it is a major source of phenotypic innovation and a mechanism of niche adaptation. For example, of particular relevance to human health is the lateral transfer of antibiotic resistance and pathogenicity determinants, leading to the emergence of pathogenic lineages. Computational identification of HGT events relies upon the investigation of sequence composition or evolutionary history of genes. Sequence composition-based ("parametric") methods search for deviations from the genomic average, whereas evolutionary history-based ("phylogenetic") approaches identify genes whose evolutionary history significantly differs from that of the host species. The evaluation and benchmarking of HGT inference methods typically rely upon simulated genomes, for which the true history is known. On real data, different methods tend to infer different HGT events, and as a result it can be difficult to ascertain all but simple and clear-cut HGT events.

Comparative genomics of closely related Salmonella enterica serovar Typhi strains reveals genome dynamics and the acquisition of novel pathogenic elements

Background

Typhoid fever is an infectious disease of global importance that is caused by Salmonella enterica subsp. enterica serovar Typhi (S. Typhi). This disease causes an estimated 200,000 deaths per year and remains a serious global health threat. S. Typhi is strictly a human pathogen, and some recovered individuals become long-term carriers who continue to shed the bacteria in their faeces, thus becoming main reservoirs of infection.

Results

A comparative genomics analysis combined with a phylogenomic analysis revealed that the strains from the outbreak and carrier were closely related with microvariations and possibly derived from a common ancestor. Additionally, the comparative genomics analysis with all of the other completely sequenced S. Typhi genomes revealed that strains BL196 and CR0044 exhibit unusual genomic variations despite S. Typhi being generally regarded as highly clonal. The two genomes shared distinct chromosomal architectures and uncommon genome features; notably, the presence of a ~10 kb novel genomic island containing uncharacterised virulence-related genes, and zot in particular. Variations were also detected in the T6SS system and genes that were related to SPI-10, insertion sequences, CRISPRs and nsSNPs among the studied genomes. Interestingly, the carrier strain CR0044 harboured far more genetic polymorphisms (83% mutant nsSNPs) compared with the closely related BL196 outbreak strain. Notably, the two highly related virulence-determinant genes, rpoS and tviE, were mutated in strains BL196 and CR0044, respectively, which revealed that the mutation in rpoS is stabilising, while that in tviE is destabilising. These microvariations provide novel insight into the optimisation of genes by the pathogens. However, the sporadic strain was found to be far more conserved compared with the others.

Conclusions

The uncommon genomic variations in the two closely related BL196 and CR0044 strains suggests that S. Typhi is more diverse than previously thought. Our study has demonstrated that the pathogen is continually acquiring new genes through horizontal gene transfer in the process of host adaptation, providing novel insight into its unusual genomic dynamics. The understanding of these strains and virulence factors, and particularly the strain that is associated with the large outbreak and the less studied asymptomatic Typhi carrier in the population, will have important impact on disease control.

Electronic supplementary material

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

HGTector: an automated method facilitating genome-wide discovery of putative horizontal gene transfers

Background

First pass methods based on BLAST match are commonly used as an initial step to separate the different phylogenetic histories of genes in microbial genomes, and target putative horizontal gene transfer (HGT) events. This will continue to be necessary given the rapid growth of genomic data and the technical difficulties in conducting large-scale explicit phylogenetic analyses. However, these methods often produce misleading results due to their inability to resolve indirect phylogenetic links and their vulnerability to stochastic events.

Results

A new computational method of rapid, exhaustive and genome-wide detection of HGT was developed, featuring the systematic analysis of BLAST hit distribution patterns in the context of a priori defined hierarchical evolutionary categories. Genes that fall beyond a series of statistically determined thresholds are identified as not adhering to the typical vertical history of the organisms in question, but instead having a putative horizontal origin. Tests on simulated genomic data suggest that this approach effectively targets atypically distributed genes that are highly likely to be HGT-derived, and exhibits robust performance compared to conventional BLAST-based approaches. This method was further tested on real genomic datasets, including Rickettsia genomes, and was compared to previous studies. Results show consistency with currently employed categories of HGT prediction methods. In-depth analysis of both simulated and real genomic data suggests that the method is notably insensitive to stochastic events such as gene loss, rate variation and database error, which are common challenges to the current methodology. An automated pipeline was created to implement this approach and was made publicly available at: https://github.com/DittmarLab/HGTector. The program is versatile, easily deployed, has a low requirement for computational resources.

Conclusions

HGTector is an effective tool for initial or standalone large-scale discovery of candidate HGT-derived genes.

Electronic supplementary material

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

Recombination blurs phylogenetic groups routine assignment in Escherichia coli: setting the record straight

The characterization of population structures plays a main role for understanding outbreaks and the dynamics of bacterial spreading. In Escherichia coli, the widely used combination of multiplex-PCR scheme together with goeBURST has some limitations. The purpose of this study is to show that the combination of different phylogenetic approaches based on concatenated sequences of MLST genes results in a more precise assignment of E. coli phylogenetic groups, complete understanding of population structure and reconstruction of ancestral clones. A collection of 80 Escherichia coli strains of different origins was analyzed following the Clermont and Doumith's multiplex-PCR schemes. Doumith's multiplex-PCR showed only 1.7% of misassignment, whereas Clermont's-2000 protocol reached 14.0%, although the discrepancies reached 30% and 38.7% respectively when recombinant C, F and E phylogroups were considered. Therefore, correct phylogroup attribution is highly variable and depends on the clonal composition of the sample. As far as population structure of these E. coli strains, including 48 E. coli genomes from GenBank, goeBURST provides a quite dispersed population structure; whereas NeighborNet approach reveals a complex population structure. MLST-based eBURST can infer different founder genotypes, for instance ST23/ST88 could be detected as the founder genotypes for STC23; however, phylogenetic reconstructions might suggest ST410 as the ancestor clone and several evolutionary trajectories with different founders. To improve our routine understanding of E. coli molecular epidemiology, we propose a strategy based on three successive steps; first, to discriminate three main groups A/B1/C, D/F/E and B2 following Doumith's protocol; second, visualization of population structure based on MLST genes according to goeBURST, using NeighborNet to establish more complex relationships among STs; and third, to perform, a cost-free characterization of evolutionary trajectories in variants emerging along the clonal expansion using parsimony methods of phylogenetic analysis.

Horizontal gene transfer in eukaryotic plant pathogens

Gene transfer has been identified as a prevalent and pervasive phenomenon and an important source of genomic innovation in bacteria. The role of gene transfer in microbial eukaryotes seems to be of a reduced magnitude but in some cases can drive important evolutionary innovations, such as new functions that underpin the colonization of different niches. The aim of this review is to summarize published cases that support the hypothesis that horizontal gene transfer (HGT) has played a role in the evolution of phytopathogenic traits in fungi and oomycetes. Our survey of the literature identifies 46 proposed cases of transfer of genes that have a putative or experimentally demonstrable phytopathogenic function. When considering the life-cycle steps through which a pathogen must progress, the majority of the HGTs identified are associated with invading, degrading, and manipulating the host. Taken together, these data suggest HGT has played a role in shaping how fungi and oomycetes colonize plant hosts.

Genomics in Coccidioides: insights into evolution, ecology, and pathogenesis

Coccidioides immitis and C. posadasii, the causative agents of the mammalian disease coccidioidomycosis, are dimorphic fungal pathogens distributed throughout desert-like environments in North and South America. Coccidioides spp. are members of the Onygenales, a diverse group of pathogenic and nonpathogenic fungi. Recently, full genomes have been published for Coccidioides and a number of other Onygenales species. Phylogenomic comparisons and additional studies in Coccidioides population genomics and gene expression have shed light on the ecology and pathogenesis of Coccidioides and the other medically important species in this clade. Observed patterns of gene family expansion/contraction and evidence of gene flow have provided insight to the evolution of Coccidioides and greatly broadened our understanding of the diversity and sources of genetic variation found in fungi. In the future, expansion of the number of sequenced isolates from all populations will allow deeper insight into the evolutionary processes that have shaped this unique human pathogen. In addition, deep sequencing of isolates from a single Coccidioides population and pairing of those data with phenotype information on growth and pathogenicity for genome-wide association analysis will allow researchers to find genes responsible for any phenotype, virulence included, that shows variation in the population.

Genome-scale comparative analysis of gene fusions, gene fissions, and the fungal tree of life

During the course of evolution genes undergo both fusion and fission by which ORFs are joined or separated. These processes can amend gene function and represent an important factor in the evolution of protein interaction networks. Gene fusions have been suggested to be useful characters for identifying evolutionary relationships because they constitute synapomorphies or cladistic characters. To investigate the fidelity of gene-fusion characters, we developed an approach for identifying differentially distributed gene fusions among whole-genome datasets: fdfBLAST. Applying this tool to the Fungi, we identified 63 gene fusions present in two or more genomes. Using a combination of phylogenetic and comparative genomic analyses, we then investigated the evolution of these genes across 115 fungal genomes, testing each gene fusion for evidence of homoplasy, including gene fission, convergence, and horizontal gene transfer. These analyses demonstrated 110 gene-fission events. We then identified a minimum of three mechanisms that drive gene fission: separation, degeneration, and duplication. These data suggest that gene fission plays an important and hitherto underestimated role in gene evolution. Gene fusions therefore are highly labile characters, and their use for polarizing evolutionary relationships, without reference to gene and species phylogenies, is limited. Accounting for these considerable sources of homoplasy, we identified fusion characters that provide support for multiple nodes in the phylogeny of the Fungi, including relationships within the deeply derived flagellum-forming fungi (i.e., the chytrids).

Man and his spaceships: Vehicles for extraterrestrial colonization?

The resiliency and adaptive ability of microbial life in real time on Earth relies heavily upon horizontal gene transfer. Based on that knowledge, how likely is earth based microbial life to colonize extraterrestrial targets such as Mars? To address this question, we consider manned and unmanned space exploration, the resident microbiota that is likely to inhabit those vehicles, the adaptive potential of that microbiota in an extraterrestrial setting especially with regards to mobile genetic elements, and the likelihood that Mars like environments could initiate and sustain colonization.

Occurrence of horizontal gene transfer of P(IB)-type ATPase genes among bacteria isolated from the uranium rich deposit of Domiasiat in North East India

Uranium (U) tolerant aerobic heterotrophs were isolated from the subsurface soils of one of the pre-mined U-rich deposits at Domiasiat located in the north-eastern part of India. On screening of genomic DNA from 62 isolates exhibiting superior U and heavy metal tolerance, 32 isolates were found to be positive for P(IB)-type ATPase genes. Phylogenetic incongruence and anomalous DNA base compositions revealed the acquisition of P(IB)-type ATPase genes by six isolates through horizontal gene transfer (HGT). Three of these instances of HGT appeared to have occurred at inter-phylum level and the other three instances indicated to have taken place at intra-phylum level. This study provides an insight into one of the possible survival strategies that bacteria might employ to adapt to environments rich in uranium and heavy metals.

Microbial stowaways: inimitable survivors or hopeless pioneers?

The resiliency of prokaryotic life has provided colonization across the globe and in the recesses of Earth's most extreme environments. Horizontal gene transfer provides access to a global bank of genetic resources that creates diversity and allows real-time adaptive potential to the clonal prokaryotic world. We assess the likelihood that this Earth-based strategy could provide survival and adaptive potential, in the case of microbial stowaways off Earth.

Detection of horizontal transfer of individual genes by anomalous oligomer frequencies

Background

Understanding the history of life requires that we understand the transfer of genetic material across phylogenetic boundaries. Detecting genes that were acquired by means other than vertical descent is a basic step in that process. Detection by discordant phylogenies is computationally expensive and not always definitive. Many have used easily computed compositional features as an alternative procedure. However, different compositional methods produce different predictions, and the effectiveness of any method is not well established.

Results

The ability of octamer frequency comparisons to detect genes artificially seeded in cyanobacterial genomes was markedly increased by using as a training set those genes that are highly conserved over all bacteria. Using a subset of octamer frequencies in such tests also increased effectiveness, but this depended on the specific target genome and the source of the contaminating genes. The presence of high frequency octamers and the GC content of the contaminating genes were important considerations. A method comprising best practices from these tests was devised, the Core Gene Similarity (CGS) method, and it performed better than simple octamer frequency analysis, codon bias, or GC contrasts in detecting seeded genes or naturally occurring transposons. From a comparison of predictions with phylogenetic trees, it appears that the effectiveness of the method is confined to horizontal transfer events that have occurred recently in evolutionary time.

Conclusions

The CGS method may be an improvement over existing surrogate methods to detect genes of foreign origin.