Horizontal gene transfer in evolution: facts and challenges

Evidence Supporting the Uptake and Genomic Incorporation of Environmental DNA in the "Ancient Asexual" Bdelloid Rotifer Philodina roseola

Increasing evidence suggests that bdelloid rotifers regularly undergo horizontal gene transfer, apparently as a surrogate mechanism of genetic exchange in the absence of true sexual reproduction, in part because of their ability to withstand desiccation. We provide empirical support for this latter hypothesis using the bdelloid Philodina roseola, which we demonstrate to readily internalize environmental DNA in contrast to a representative monogonont rotifer (Brachionus rubens), which, like other monogononts, is facultative sexual and cannot withstand desiccation. In addition, environmental DNA that was more similar to the host DNA was retained more often and for a longer period of time. Indirect evidence (increased variance in the reproductive output of the untreated F1 generation) suggests that environmental DNA can be incorporated into the genome during desiccation and is thus heritable. Our observed fitness effects agree with sexual theory and also occurred when the animals were desiccated in groups (thereby acting as DNA donors), but not individually, indicating the mechanism could occur in nature. Thus, although DNA uptake and its genomic incorporation appears proximally related to anhydrobiosis in bdelloids, it might also facilitate accidental genetic exchange with closely related taxa, thereby maintaining higher levels of genetic diversity than is otherwise expected for this group of "ancient asexuals".

Escherichia coli out in the cold: Dissemination of human-derived bacteria into the Antarctic microbiome

Discharge of untreated sewage into Antarctic environments presents a risk of introducing non-native microorganisms, but until now, adverse consequences have not been conclusively identified. Here we show that sewage disposal introduces human derived Escherichia coli carrying mobile genetic elements and virulence traits with the potential to affect the diversity and evolution of native Antarctic microbial communities. We compared E. coli recovered from environmental and animal sources in Antarctica to a reference collection of E. coli from humans and non-Antarctic animals. The distribution of phylogenetic groups and frequency of 11 virulence factors amongst the Antarctic isolates were characteristic of E. coli strains more commonly associated with humans. The rapidly emerging E. coli ST131 and ST95 clones were found amongst the Antarctic isolates, and ST95 was the predominant E. coli recovered from Weddell seals. Class 1 integrons were found in 15% of the Antarctic E. coli with 4 of 5 identified gene cassette arrays containing antibiotic resistance genes matching those common in clinical contexts. Disposing untreated sewage into the Antarctic environment does disseminate non-native microorganisms, but the extent of this impact and implications for Antarctic ecosystem health are, as yet, poorly understood.

Exploring lateral genetic transfer among microbial genomes using TF-IDF

Many microbes can acquire genetic material from their environment and incorporate it into their genome, a process known as lateral genetic transfer (LGT). Computational approaches have been developed to detect genomic regions of lateral origin, but typically lack sensitivity, ability to distinguish donor from recipient, and scalability to very large datasets. To address these issues we have introduced an alignment-free method based on ideas from document analysis, term frequency-inverse document frequency (TF-IDF). Here we examine the performance of TF-IDF on three empirical datasets: 27 genomes of Escherichia coli and Shigella, 110 genomes of enteric bacteria, and 143 genomes across 12 bacterial and three archaeal phyla. We investigate the effect of k-mer size, gap size and delineation of groups on the inference of genomic regions of lateral origin, finding an interplay among these parameters and sequence divergence. Because TF-IDF identifies donor groups and delineates regions of lateral origin within recipient genomes, aggregating these regions by gene enables us to explore, for the first time, the mosaic nature of lateral genes including the multiplicity of biological sources, ancestry of transfer and over-writing by subsequent transfers. We carry out Gene Ontology enrichment tests to investigate which biological processes are potentially affected by LGT.

Xenogeneic Silencing and Its Impact on Bacterial Genomes

The H-NS (heat-stable nucleoid structuring) protein affects both nucleoid compaction and global gene regulation. H-NS appears to act primarily as a silencer of AT-rich genetic material acquired by horizontal gene transfer. As such, it is key in the regulation of most genes involved in virulence and in adaptation to new environmental niches. Here we review recent progress in understanding the biochemistry of H-NS and how xenogeneic silencing affects bacterial evolution. We highlight the strengths and weaknesses of some of the models proposed in H-NS-mediated nucleoprotein complex formation. Based on recent single-molecule studies, we also propose a novel mode of DNA compaction by H-NS termed intrabridging to explain over two decades of observations of the H-NS molecule.

A functional difference between native and horizontally acquired genes in bdelloid rotifers

The form of RNA processing known as SL trans-splicing involves the transfer of a short conserved sequence, the spliced leader (SL), from a noncoding SL RNA to the 5' ends of mRNA molecules. SL trans-splicing occurs in several animal taxa, including bdelloid rotifers (Rotifera, Bdelloidea). One striking feature of these aquatic microinvertebrates is the large proportion of foreign genes, i.e. those acquired by horizontal gene transfer from other organisms, in their genomes. However, whether such foreign genes behave similarly to native genes has not been tested in bdelloids or any other animal. We therefore used a combination of experimental and computational methods to examine whether transcripts of foreign genes in bdelloids were SL trans-spliced, like their native counterparts. We found that many foreign transcripts contain SLs, use similar splice acceptor sequences to native genes, and are able to undergo alternative trans-splicing. However, a significantly lower proportion of foreign mRNAs contains SL sequences than native transcripts. This demonstrates a novel functional difference between foreign and native genes in bdelloids and suggests that SL trans-splicing is not essential for the expression of foreign genes, but is acquired during their domestication.

Impact of genomics on the understanding of microbial evolution and classification: the importance of Darwin's views on classification

Analyses of genome sequences, by some approaches, suggest that the widespread occurrence of horizontal gene transfers (HGTs) in prokaryotes disguises their evolutionary relationships and have led to questioning of the Darwinian model of evolution for prokaryotes. These inferences are critically examined in the light of comparative genome analysis, characteristic synapomorphies, phylogenetic trees and Darwin's views on examining evolutionary relationships. Genome sequences are enabling discovery of numerous molecular markers (synapomorphies) such as conserved signature indels (CSIs) and conserved signature proteins (CSPs), which are distinctive characteristics of different prokaryotic taxa. Based on these molecular markers, exhibiting high degree of specificity and predictive ability, numerous prokaryotic taxa of different ranks, currently identified based on the 16S rRNA gene trees, can now be reliably demarcated in molecular terms. Within all studied groups, multiple CSIs and CSPs have been identified for successive nested clades providing reliable information regarding their hierarchical relationships and these inferences are not affected by HGTs. These results strongly support Darwin's views on evolution and classification and supplement the current phylogenetic framework based on 16S rRNA in important respects. The identified molecular markers provide important means for developing novel diagnostics, therapeutics and for functional studies providing important insights regarding prokaryotic taxa.

Horizontal gene transfer events reshape the global landscape of arm race between viruses and homo sapiens

Horizontal gene transfer (HGT) drives the evolution of recipient organism particularly if it provides a novel function which enhances the fitness or its adaption to the environment. Virus-host co-evolution is attractive for studying co-evolutionary processes, since viruses strictly replicate inside of the host cells and thus their evolution is inexorably tangled with host biology. HGT, as a mechanism of co-evolution between human and viruses, has been widely documented, however, the roles HGT play during the interaction between human and viruses are still in their infancy. In this study, we performed a comprehensive analysis on the genes horizontally transferred between viruses and their corresponding human hosts. Our study suggests that the HGT genes in human are predominantly enriched in immune related GO terms while viral HGT genes are tend to be encoded by viruses which promote the invasion of immune system of hosts. Based on our results, it gives us a hint about the evolution trajectory of HGT events. Overall, our study suggests that the HGT between human and viruses are highly relevant to immune interaction and probably reshaped the arm race between hosts and viruses.

Recent events dominate interdomain lateral gene transfers between prokaryotes and eukaryotes and, with the exception of endosymbiotic gene transfers, few ancient transfer events persist

While there is compelling evidence for the impact of endosymbiotic gene transfer (EGT; transfer from either mitochondrion or chloroplast to the nucleus) on genome evolution in eukaryotes, the role of interdomain transfer from bacteria and/or archaea (i.e. prokaryotes) is less clear. Lateral gene transfers (LGTs) have been argued to be potential sources of phylogenetic information, particularly for reconstructing deep nodes that are difficult to recover with traditional phylogenetic methods. We sought to identify interdomain LGTs by using a phylogenomic pipeline that generated 13 465 single gene trees and included up to 487 eukaryotes, 303 bacteria and 118 archaea. Our goals include searching for LGTs that unite major eukaryotic clades, and describing the relative contributions of LGT and EGT across the eukaryotic tree of life. Given the difficulties in interpreting single gene trees that aim to capture the approximately 1.8 billion years of eukaryotic evolution, we focus on presence-absence data to identify interdomain transfer events. Specifically, we identify 1138 genes found only in prokaryotes and representatives of three or fewer major clades of eukaryotes (e.g. Amoebozoa, Archaeplastida, Excavata, Opisthokonta, SAR and orphan lineages). The majority of these genes have phylogenetic patterns that are consistent with recent interdomain LGTs and, with the notable exception of EGTs involving photosynthetic eukaryotes, we detect few ancient interdomain LGTs. These analyses suggest that LGTs have probably occurred throughout the history of eukaryotes, but that ancient events are not maintained unless they are associated with endosymbiotic gene transfer among photosynthetic lineages.

The Geobacillus Pan-Genome: Implications for the Evolution of the Genus

The genus Geobacillus is comprised of a diverse group of spore-forming Gram-positive thermophilic bacterial species and is well known for both its ecological diversity and as a source of novel thermostable enzymes. Although the mechanisms underlying the thermophilicity of the organism and the thermostability of its macromolecules are reasonably well understood, relatively little is known of the evolutionary mechanisms, which underlie the structural and functional properties of members of this genus. In this study, we have compared 29 Geobacillus genomes, with a specific focus on the elements, which comprise the conserved core and flexible genomes. Based on comparisons of conserved core and flexible genomes, we present evidence of habitat delineation with specific Geobacillus genomes linked to specific niches. Our analysis revealed that Geobacillus and Anoxybacillus share a high proportion of genes. Moreover, the results strongly suggest that horizontal gene transfer is a major factor deriving the evolution of Geobacillus from Bacillus, with genetic contributions from other phylogenetically distant taxa.

Evolutionary change and phylogenetic relationships in light of horizontal gene transfer

Horizontal gene transfer has, over the past 25 years, become a part of evolutionary thinking. In the present paper I discuss horizontal gene transfer (HGT) in relation to contingency, natural selection, evolutionary change speed and the Tree-of-Life endeavour, with the aim of contributing to the understanding of the role of HGT in evolutionary processes. In addition, the challenges that HGT imposes on the current view of evolution are emphasized.

Comparisons of the transferability of plasmids pCAR1, pB10, R388, and NAH7 among Pseudomonas putida at different cell densities

The transferability of plasmids pCAR1, pB10, R388, and NAH7 was compared using the same donor-recipient system at different cell density combinations in liquid or on a solid surface. pCAR1 was efficiently transferred in liquid, whereas the other plasmids were preferentially transferred on a solid surface. Difference of liquid or solid affected the transfer frequency especially at lower cell densities.

Horizontal transfer of DNA methylation patterns into bacterial chromosomes

Horizontal gene transfer (HGT) is the non-inherited acquisition of novel DNA sequences. HGT is common and important in bacteria because it enables the rapid generation of new phenotypes such as antibiotic resistance. Here we show that in vivo and in vitro DNA methylation patterns can be horizontally transferred into bacterial chromosomes to program cell phenotypes. The experiments were performed using a synthetic system in Escherichia coli where different DNA methylation patterns within the cis-regulatory sequence of the agn43 gene turn on or off a fluorescent reporter (CFP). With this system we demonstrated that DNA methylation patterns not only accompany the horizontal transfer of genes into the bacterial cytoplasm but can be transferred into chromosomes by: (i) bacteriophage P1 transduction; and (ii) transformation of extracellular synthetic DNA. We also modified the experimental system by replacing CFP with the SgrS small RNA, which regulates glucose and methyl α-D-glucoside uptake, and showed that horizontally acquired DNA methylation patterns can increase or decrease cell fitness. That is, horizontally acquired DNA methylation patterns can result in the selection for and against cells that have HGT. Findings from these proof-of-concept experiments have applications in synthetic biology and potentially broad implications for bacterial adaptation and evolution.

Abiotic Gene Transfer: Rare or Rampant?

Phylogenetic studies reveal that horizontal gene transfer (HGT) plays a prominent role in evolution and genetic variability of life. Five biotic mechanisms of HGT among prokaryotic organisms have been extensively characterized: conjugation, competence, transduction, gene transfer agent particles, and transitory fusion with recombination, but it is not known whether they can account for all natural HGT. It is even less clear how HGT could have occurred before any of these mechanisms had developed. Here, we consider contemporary conditions and experiments on microorganisms to estimate possible roles of abiotic HGT-currently and throughout evolution. Candidate mechanisms include freeze-and-thaw, microbeads-agitation, and electroporation-based transformation, and we posit that these laboratory techniques have analogues in nature acting as mechanisms of abiotic HGT: freeze-and-thaw cycles in polar waters, agitation by sand at foreshores and riverbeds, and lightning-triggered electroporation in near-surface aqueous habitats. We derive conservative order-of-magnitude estimates for rates of microorganisms subjected to freeze-and-thaw cycles, sand agitation, and lightning-triggered electroporation, at 10²⁴, 10¹⁹, and 10¹⁷ per year, respectively. Considering the yield of viable transformants, which is by far the highest in electroporation, we argue this may still favor lightning-triggered transformation over the other two mechanisms. Electroporation-based gene transfer also appears to be the most general of these abiotic candidates, and perhaps even of all known HGT mechanisms. Future studies should provide improved estimates of gene transfer rates and cell viability, currently and in the past, but to assess the importance of abiotic HGT in nature will likely require substantial progress-also in knowledge of biotic HGT.

Biochemistry and physiology within the framework of the extended synthesis of evolutionary biology

Functional biologists, like Claude Bernard, ask "How?", meaning that they investigate the mechanisms underlying the emergence of biological functions (proximal causes), while evolutionary biologists, like Charles Darwin, asks "Why?", meaning that they search the causes of adaptation, survival and evolution (remote causes). Are these divergent views on what is life? The epistemological role of functional biology (molecular biology, but also biochemistry, physiology, cell biology and so forth) appears essential, for its capacity to identify several mechanisms of natural selection of new characters, individuals and populations. Nevertheless, several issues remain unsolved, such as orphan metabolic activities, i.e., adaptive functions still missing the identification of the underlying genes and proteins, and orphan genes, i.e., genes that bear no signature of evolutionary history, yet provide an organism with improved adaptation to environmental changes. In the framework of the Extended Synthesis, we suggest that the adaptive roles of any known function/structure are reappraised in terms of their capacity to warrant constancy of the internal environment (homeostasis), a concept that encompasses both proximal and remote causes.

Deep sequencing revealed molecular signature of horizontal gene transfer of plant like transcripts in the mosquito Anopheles culicifacies: an evolutionary puzzle

In prokaryotes, horizontal gene transfer (HGT) has been regarded as an important evolutionary drive to acquire and retain beneficial genes for their survival in diverse ecologies. However, in eukaryotes, the functional role of HGTs remains questionable, although current genomic tools are providing increased evidence of acquisition of novel traits within non-mating metazoan species. Here, we provide another transcriptomic evidence for the acquisition of massive plant genes in the mosquito, Anopheles culicifacies. Our multiple experimental validations including genomic PCR, RT-PCR, real-time PCR, immuno-blotting and immuno-florescence microscopy, confirmed that plant like transcripts (PLTs) are of mosquito origin and may encode functional proteins. A comprehensive molecular analysis of the PLTs and ongoing metagenomic analysis of salivary microbiome provide initial clues that mosquitoes may have survival benefits through the acquisition of nuclear as well as chloroplast encoded plant genes. Our findings of PLTs further support the similar questionable observation of HGTs in other higher organisms, which is still a controversial and debatable issue in the community of evolutionists. We believe future understanding of the underlying mechanism of the feeding associated molecular responses may shed new insights in the functional role of PLTs in the mosquito.

Inferring orthologous gene regulatory networks using interspecies data fusion

Motivation: The ability to jointly learn gene regulatory networks (GRNs) in, or leverage GRNs between related species would allow the vast amount of legacy data obtained in model organisms to inform the GRNs of more complex, or economically or medically relevant counterparts. Examples include transferring information from Arabidopsis thaliana into related crop species for food security purposes, or from mice into humans for medical applications. Here we develop two related Bayesian approaches to network inference that allow GRNs to be jointly inferred in, or leveraged between, several related species: in one framework, network information is directly propagated between species; in the second hierarchical approach, network information is propagated via an unobserved ‘hypernetwork’. In both frameworks, information about network similarity is captured via graph kernels, with the networks additionally informed by species-specific time series gene expression data, when available, using Gaussian processes to model the dynamics of gene expression.

Results: Results on in silico benchmarks demonstrate that joint inference, and leveraging of known networks between species, offers better accuracy than standalone inference. The direct propagation of network information via the non-hierarchical framework is more appropriate when there are relatively few species, while the hierarchical approach is better suited when there are many species. Both methods are robust to small amounts of mislabelling of orthologues. Finally, the use of Saccharomyces cerevisiae data and networks to inform inference of networks in the budding yeast Schizosaccharomyces pombe predicts a novel role in cell cycle regulation for Gas1 (SPAC19B12.02c), a 1,3-beta-glucanosyltransferase.

Availability and implementation: MATLAB code is available from http://go.warwick.ac.uk/systemsbiology/software/.

Contact:d.l.wild@warwick.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.

Horizontal functional gene transfer from bacteria to fishes

Invertebrates can acquire functional genes via horizontal gene transfer (HGT) from bacteria but fishes are not known to do so. We provide the first reliable evidence of one HGT event from marine bacteria to fishes. The HGT appears to have occurred after emergence of the teleosts. The transferred gene is expressed and regulated developmentally. Its successful integration and expression may change the genetic and metabolic repertoire of fishes. In addition, this gene contains conserved domains and similar tertiary structures in fishes and their putative donor bacteria. Thus, it may function similarly in both groups. Evolutionary analyses indicate that it evolved under purifying selection, further indicating its conserved function. We document the first likely case of HGT of functional gene from prokaryote to fishes. This discovery certifies that HGT can influence vertebrate evolution.

A reconstruction problem for a class of phylogenetic networks with lateral gene transfers

Background

Lateral, or Horizontal, Gene Transfers are a type of asymmetric evolutionary events where genetic material is transferred from one species to another. In this paper we consider LGT networks, a general model of phylogenetic networks with lateral gene transfers which consist, roughly, of a principal rooted tree with its leaves labelled on a set of taxa, and a set of extra secondary arcs between nodes in this tree representing lateral gene transfers. An LGT network gives rise in a natural way to a principal phylogenetic subtree and a set of secondary phylogenetic subtrees, which, roughly, represent, respectively, the main line of evolution of most genes and the secondary lines of evolution through lateral gene transfers.

Results

We introduce a set of simple conditions on an LGT network that guarantee that its principal and secondary phylogenetic subtrees are pairwise different and that these subtrees determine, up to isomorphism, the LGT network. We then give an algorithm that, given a set of pairwise different phylogenetic trees \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_0,T_1,\ldots ,T_k$$\end{document}T0,T1,…,Tk on the same set of taxa, outputs, when it exists, the LGT network that satisfies these conditions and such that its principal phylogenetic tree is \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_0$$\end{document}T0 and its secondary phylogenetic trees are \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_1,\ldots ,T_k$$\end{document}T1,…,Tk.

Electronic supplementary material

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

Recombination and horizontal transfer of nodulation and ACC deaminase (acdS) genes within Alpha- and Betaproteobacteria nodulating legumes of the Cape Fynbos biome

The goal of this work is to study the evolution and the degree of horizontal gene transfer (HGT) within rhizobial genera of both Alphaproteobacteria (Mesorhizobium, Rhizobium) and Betaproteobacteria (Burkholderia), originating from South African Fynbos legumes. By using a phylogenetic approach and comparing multiple chromosomal and symbiosis genes, we revealed conclusive evidence of high degrees of horizontal transfer of nodulation genes among closely related species of both groups of rhizobia, but also among species with distant genetic backgrounds (Rhizobium and Mesorhizobium), underscoring the importance of lateral transfer of symbiosis traits as an important evolutionary force among rhizobia of the Cape Fynbos biome. The extensive exchange of symbiosis genes in the Fynbos is in contrast with a lack of significant events of HGT among Burkholderia symbionts from the South American Cerrado and Caatinga biome. Furthermore, homologous recombination among selected housekeeping genes had a substantial impact on sequence evolution within Burkholderia and Mesorhizobium. Finally, phylogenetic analyses of the non-symbiosis acdS gene in Mesorhizobium, a gene often located on symbiosis islands, revealed distinct relationships compared to the chromosomal and symbiosis genes, suggesting a different evolutionary history and independent events of gene transfer. The observed events of HGT and incongruence between different genes necessitate caution in interpreting topologies from individual data types.

Correlation between the green-island phenotype and Wolbachia infections during the evolutionary diversification of Gracillariidae leaf-mining moths

Internally feeding herbivorous insects such as leaf miners have developed the ability to manipulate the physiology of their host plants in a way to best meet their metabolic needs and compensate for variation in food nutritional composition. For instance, some leaf miners can induce green-islands on yellow leaves in autumn, which are characterized by photosynthetically active green patches in otherwise senescing leaves. It has been shown that endosymbionts, and most likely bacteria of the genus Wolbachia, play an important role in green-island induction in the apple leaf-mining moth Phyllonorycter blancardella. However, it is currently not known how widespread is this moth-Wolbachia-plant interaction. Here, we studied the co-occurrence between Wolbachia and the green-island phenotype in 133 moth specimens belonging to 74 species of Lepidoptera including 60 Gracillariidae leaf miners. Using a combination of molecular phylogenies and ecological data (occurrence of green-islands), we show that the acquisitions of the green-island phenotype and Wolbachia infections have been associated through the evolutionary diversification of Gracillariidae. We also found intraspecific variability in both green-island formation and Wolbachia infection, with some species being able to form green-islands without being infected by Wolbachia. In addition, Wolbachia variants belonging to both A and B supergroups were found to be associated with green-island phenotype suggesting several independent origins of green-island induction. This study opens new prospects and raises new questions about the ecology and evolution of the tripartite association between Wolbachia, leaf miners, and their host plants.

The evolutionary dynamics of tRNA-gene copy number and codon-use in E. coli

Background

The introduction of foreign DNA by Lateral Gene Transfer (LGT) can quickly and drastically alter genome composition. Problems can arise if the genes introduced by LGT use codons that are not suited to the host’s translational machinery. Here we investigate compensatory adaptation of E. coli in response to the introduction of large volumes of codons that are rarely used by the host genome.

Results

We analyze genome sequences from the E. coli/Shigella complex, and find that certain tRNA genes are present in multiple copies in two pathogenic Shigella and O157:H7 subgroups of E. coli. Furthermore, we show that the codons that correspond to these multi-copy number tRNA genes are enriched in the high copy number Selfish Genetic Elements (SGE’s) in Shigella and laterally introduced genes in O157:H7. We analyze the duplicate copies and find evidence for the selective retention of tRNA genes introduced by LGT in response to the changed codon content of the genome.

Conclusion

These data support a model where the relatively rapid influx of LGT genes and SGE’s introduces a large number of genes maladapted to the host’s translational machinery. Under these conditions, it becomes advantageous for the host to retain tRNA genes that are required for the incorporation of amino acids at these codons. Subsequently, the increased number of copies of these specific tRNA genes adjusts the cellular tRNA pool to the demands set by global shifts in codon usage.

Electronic supplementary material

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

Antibiotics, Antibiotic Resistance Genes, and Bacterial Community Composition in Fresh Water Aquaculture Environment in China

Environmental antibiotic resistance has drawn increasing attention due to its great threat to human health. In this study, we investigated concentrations of antibiotics (tetracyclines, sulfonamides and (fluoro)quinolones) and abundances of antibiotic resistance genes (ARGs), including tetracycline resistance genes, sulfonamide resistance genes, and plasmid-mediated quinolone resistance genes, and analyzed bacterial community composition in aquaculture environment in Guangdong, China. The concentrations of sulfametoxydiazine, sulfamethazine, sulfamethoxazole, oxytetracycline, chlorotetracycline, doxycycline, ciprofloxacin, norfloxacin, and enrofloxacin were as high as 446 μg kg(-1) and 98.6 ng L(-1) in sediment and water samples, respectively. The relative abundances (ARG copies/16S ribosomal RNA (rRNA) gene copies) of ARGs (sul1, sul2, sul3, tetM, tetO, tetW, tetS, tetQ, tetX, tetB/P, qepA, oqxA, oqxB, aac(6')-Ib, and qnrS) were as high as 2.8 × 10(-2). The dominant phyla were Proteobacteria, Bacteroidetes, and Firmicutes in sediment samples and Proteobacteria, Actinobacteria and Bacteroidetes in water samples. The genera associated with pathogens were also observed, such as Acinetobacter, Arcobacter, and Clostridium. This study comprehensively investigated antibiotics, ARGs, and bacterial community composition in aquaculture environment in China. The results indicated that fish ponds are reservoirs of ARGs and the presence of potential resistant and pathogen-associated taxonomic groups in fish ponds might imply the potential risk to human health.

Ecological Specialization of Two Photobiont-Specific Maritime Cyanolichen Species of the Genus Lichina

All fungi in the class Lichinomycetes are lichen-forming and exclusively associate with cyanobacteria. Two closely related maritime species of the genus Lichina (L. confinis and L. pygmaea) show similar distribution ranges in the Northeast Atlantic, commonly co-occurring at the same rocky shores but occupying different littoral zones. By means of 16S rRNA and phycocyanin operon markers we studied a) the phylogenetic relationships of cyanobionts associated with these species, b) the match of divergence times between both symbionts, and c) whether Lichina species differ in photobiont association and in how geography and ecology affect selectivity. The cyanobionts studied are closely related to both marine and freshwater strains of the genus Rivularia. We found evidence of a high specificity to particular cyanobiont lineages in both species: Lichina pygmaea and L. confinis incorporate specific lineages of Rivularia that do not overlap at the haplotype nor the OTU levels. Dating divergences of the fungal and cyanobacterial partners revealed an asynchronous origin of both lineages. Within each fungal species, selectivity varied across the studied area, influenced by environmental conditions (both atmospheric and marine), although patterns were highly correlated between both lichen taxa. Ecological speciation due to the differential association of photobionts to each littoral zone is suspected to have occurred in marine Lichina.

Pre_GI: a global map of ontological links between horizontally transferred genomic islands in bacterial and archaeal genomes

The Predicted Genomic Islands database (Pre_GI) is a comprehensive repository of prokaryotic genomic islands (islands, GIs) freely accessible at http://pregi.bi.up.ac.za/index.php. Pre_GI, Version 2015, catalogues 26 744 islands identified in 2407 bacterial/archaeal chromosomes and plasmids. It provides an easy-to-use interface which allows users the ability to query against the database with a variety of fields, parameters and associations. Pre_GI is constructed to be a web-resource for the analysis of ontological roads between islands and cartographic analysis of the global fluxes of mobile genetic elements through bacterial and archaeal taxonomic borders. Comparison of newly identified islands against Pre_GI presents an alternative avenue to identify their ontology, origin and relative time of acquisition. Pre_GI aims to aid research on horizontal transfer events and materials through providing data and tools for holistic investigation of migration of genes through ecological niches and taxonomic boundaries.

Reconstructible phylogenetic networks: do not distinguish the indistinguishable

Phylogenetic networks represent the evolution of organisms that have undergone reticulate events, such as recombination, hybrid speciation or lateral gene transfer. An important way to interpret a phylogenetic network is in terms of the trees it displays, which represent all the possible histories of the characters carried by the organisms in the network. Interestingly, however, different networks may display exactly the same set of trees, an observation that poses a problem for network reconstruction: from the perspective of many inference methods such networks are "indistinguishable". This is true for all methods that evaluate a phylogenetic network solely on the basis of how well the displayed trees fit the available data, including all methods based on input data consisting of clades, triples, quartets, or trees with any number of taxa, and also sequence-based approaches such as popular formalisations of maximum parsimony and maximum likelihood for networks. This identifiability problem is partially solved by accounting for branch lengths, although this merely reduces the frequency of the problem. Here we propose that network inference methods should only attempt to reconstruct what they can uniquely identify. To this end, we introduce a novel definition of what constitutes a uniquely reconstructible network. For any given set of indistinguishable networks, we define a canonical network that, under mild assumptions, is unique and thus representative of the entire set. Given data that underwent reticulate evolution, only the canonical form of the underlying phylogenetic network can be uniquely reconstructed. While on the methodological side this will imply a drastic reduction of the solution space in network inference, for the study of reticulate evolution this is a fundamental limitation that will require an important change of perspective when interpreting phylogenetic networks.

Inter- and intra-domain horizontal gene transfer, gain-loss asymmetry and positive selection mark the evolutionary history of the CBM14 family

Protein-carbohydrate interactions are ubiquitous in nature and at the core of many physiological processes of profound importance to health and disease. Specificity in protein-carbohydrate interactions is conferred by carbohydrate-binding modules (CBMs) that can accurately discriminate among the multitude of saccharides found in nature, thus targeting proteins to their particular substrates. Family 14 carbohydrate-binding modules (CBM14s), more specifically, are short modules that bind explicitly to chitin, the second most abundant carbohydrate in nature. Although considerable effort has been placed in elucidating protein-carbohydrate interactions at the molecular level for biological and biotechnological applications, in contrast the evolutionary relationships among these modules are minimally understood. Using the CBM14 family as an example, here we describe one of the first global molecular evolutionary analyses of a CBM family across all domains of life, with an emphasis on its origin, taxonomic distribution and pattern of diversification as a result of gene and module duplication, and positive selection. Our genome-wide searches recovered an impressive number of CBM14s from diverse lineages across nearly all domains of life. However, their highly disseminated distribution in taxa outside the Opisthokonta group strongly suggests a later evolutionary origin and elevated rates of inter- and intra-domain horizontal gene transfer. Moreover, accelerated rates of asymmetric gains and losses reveal a dynamic mode of birth-and-death evolution, whereas positive selection acting on paralogous CBM14-containing proteins suggest changes in substrate specificity and an increase in the functional promiscuity of this ancient CBM family. The importance of these results is discussed.

The evolutionary fate of the horizontally transferred agrobacterial mikimopine synthase gene in the genera Nicotiana and Linaria

Few cases of spontaneously horizontally transferred bacterial genes into plant genomes have been described to date. The occurrence of horizontally transferred genes from the T-DNA of Agrobacterium rhizogenes into the plant genome has been reported in the genus Nicotiana and in the species Linaria vulgaris. Here we compare patterns of evolution in one of these genes (a gene encoding mikimopine synthase, mis) following three different events of horizontal gene transfer (HGT). As this gene plays an important role in Agrobacterium, and there are known cases showing that genes from pathogens can acquire plant protection function, we hypothesised that in at least some of the studied species we will find signs of selective pressures influencing mis sequence. The mikimopine synthase (mis) gene evolved in a different manner in the branch leading to Nicotiana tabacum and N. tomentosiformis, in the branch leading to N. glauca and in the genus Linaria. Our analyses of the genus Linaria suggest that the mis gene began to degenerate soon after the HGT. In contrast, in the case of N. glauca, the mis gene evolved under significant selective pressures. This suggests a possible role of mikimopine synthase in current N. glauca and its ancestor(s). In N. tabacum and N. tomentosiformis, the mis gene has a common frameshift mutation that disrupted its open reading frame. Interestingly, our results suggest that in spite of the frameshift, the mis gene could evolve under selective pressures. This sequence may still have some regulatory role at the RNA level as suggested by coverage of this sequence by small RNAs in N. tabacum.

Natural transformation and genome evolution in Streptococcus pneumoniae

Streptococcus pneumoniae is a frequent colonizer of the human nasopharynx that has the potential to cause severe infections such as pneumonia, bacteremia and meningitis. Despite considerable efforts to reduce the burden of pneumococcal disease, it continues to be a major public health problem. After the Second World War, antimicrobial therapy was introduced to fight pneumococcal infections, followed by the first effective vaccines more than half a century later. These clinical interventions generated a selection pressure that drove the evolution of vaccine-escape mutants and strains that were highly resistant against antibiotics. The remarkable ability of S. pneumoniae to acquire drug resistance and evade vaccine pressure is due to its recombination-mediated genetic plasticity. S. pneumoniae is competent for natural genetic transformation, a property that enables the pneumococcus to acquire new traits by taking up naked DNA from the environment and incorporating it into its genome through homologous recombination. In the present paper, we review current knowledge on pneumococcal transformation, and discuss how the pneumococcus uses this mechanism to adapt and survive under adverse and fluctuating conditions.

Integron associated mobile genes: Just a collection of plug in apps or essential components of cell network hardware?

Lateral gene transfer (LGT) impacts on the evolution of prokaryotes in both the short and long-term. The short-term impacts of mobilized genes are a concern to humans since LGT explains the global rise of multi drug resistant pathogens seen in the past 70 years. However, LGT has been a feature of prokaryotes from the earliest days of their existence and the concept of a bifurcating tree of life is not entirely applicable to prokaryotes since most genes in extant prokaryotic genomes have probably been acquired from other lineages. Successful transfer and maintenance of a gene in a new host is understandable if it acts independently of cell networks and confers an advantage. Antibiotic resistance provides an example of this whereby a gene can be advantageous in virtually any cell across broad species backgrounds. In a longer evolutionary context however laterally transferred genes can be assimilated into even essential cell networks. How this happens is not well understood and we discuss recent work that identifies a mobile gene, unique to a cell lineage, which is detrimental to the cell when lost. We also present some additional data and believe our emerging model will be helpful in understanding how mobile genes integrate into cell networks.

Functional phylogenomics analysis of bacteria and archaea using consistent genome annotation with UniFam

BACKGROUND

Phylogenetic studies have provided detailed knowledge on the evolutionary mechanisms of genes and species in Bacteria and Archaea. However, the evolution of cellular functions, represented by metabolic pathways and biological processes, has not been systematically characterized. Many clades in the prokaryotic tree of life have now been covered by sequenced genomes in GenBank. This enables a large-scale functional phylogenomics study of many computationally inferred cellular functions across all sequenced prokaryotes.

RESULTS

A total of 14,727 GenBank prokaryotic genomes were re-annotated using a new protein family database, UniFam, to obtain consistent functional annotations for accurate comparison. The functional profile of a genome was represented by the biological process Gene Ontology (GO) terms in its annotation. The GO term enrichment analysis differentiated the functional profiles between selected archaeal taxa. 706 prokaryotic metabolic pathways were inferred from these genomes using Pathway Tools and MetaCyc. The consistency between the distribution of metabolic pathways in the genomes and the phylogenetic tree of the genomes was measured using parsimony scores and retention indices. The ancestral functional profiles at the internal nodes of the phylogenetic tree were reconstructed to track the gains and losses of metabolic pathways in evolutionary history.

CONCLUSIONS

Our functional phylogenomics analysis shows divergent functional profiles of taxa and clades. Such function-phylogeny correlation stems from a set of clade-specific cellular functions with low parsimony scores. On the other hand, many cellular functions are sparsely dispersed across many clades with high parsimony scores. These different types of cellular functions have distinct evolutionary patterns reconstructed from the prokaryotic tree.

Examination of exhaustive cloning attempts reveals that C. elegans piRNAs, transposons, and repeat sequences are efficiently cloned in yeast, but not in bacteria

Genome sequencing requires insertion of random fragments of the sequenced organism’s DNA into a unicellular host, most often Escherichia coli bacteria. This manipulation was found in the past to be analogous to naturally occurring horizontal gene transfer, and moreover has proved valuable to understanding toxicity of foreign genetic elements to E. coli. Sequencing of the Caenorhabditis elegans genome was similarly achieved via DNA transformation into E. coli. However, numerous attempts have proven a significant percentage of the genome unclonable using bacteria, although clonable via yeast. We examined the genomic segments that were not clonable in bacteria but were clonable in yeast, and observed that, in line with previous hypotheses, such sequences are more repetitive on average compared with the entire C. elegans genome. In addition, we found that these gap-sequences encode significantly more for DNA transposons. Surprisingly, we discovered that although the vast majority of the C. elegans genome is clonable in bacteria (77.5%), almost all the thousands of sequences that encode for PIWI-interacting small RNAs, or 21U-RNAs (91.6%) were only clonable in yeast. These results might help understanding why most piRNAs in C. elegans are physically clustered on particular loci on chromosome IV. In worms and in a large number of other organisms, piRNAs serve to distinguish “Self” from “Non-Self” sequences, and thus to protect the integrity of the genome against foreign genetic elements, such as transposons. We discuss the possible implications of these discoveries.

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.

The utility of metabolomics in natural product and biomarker characterization

BACKGROUND

Metabolomics is a well-established rapidly developing research field involving quantitative and qualitative metabolite assessment within biological systems. Recent improvements in metabolomics technologies reveal the unequivocal value of metabolomics tools in natural products discovery, gene-function analysis, systems biology and diagnostic platforms.

SCOPE OF REVIEW

We review here some of the prominent metabolomics methodologies employed in data acquisition and analysis of natural products and disease-related biomarkers.

MAJOR CONCLUSIONS

This review demonstrates that metabolomics represents a highly adaptable technology with diverse applications ranging from environmental toxicology to disease diagnosis. Metabolomic analysis is shown to provide a unique snapshot of the functional genetic status of an organism by examining its biochemical profile, with relevance toward resolving phylogenetic associations involving horizontal gene transfer and distinguishing subgroups of genera possessing high genetic homology, as well as an increasing role in both elucidating biosynthetic transformations of natural products and detecting preclinical biomarkers of numerous disease states.

GENERAL SIGNIFICANCE

This review expands the interest in multiplatform combinatorial metabolomic analysis. The applications reviewed range from phylogenetic assignment, biosynthetic transformations of natural products, and the detection of preclinical biomarkers.

Potentially novel copper resistance genes in copper-enriched activated sludge revealed by metagenomic analysis

In this study, we utilized the Illumina high-throughput metagenomic approach to investigate diversity and abundance of both microbial community and copper resistance genes (CuRGs) in activated sludge (AS) which was enriched under copper selective stress up to 800 mg/L. The raw datasets (~3.5 Gb for each sample, i.e., the copper-enriched AS and the control AS) were merged and normalized for the BLAST analyses against the SILVA SSU rRNA gene database and self-constructed copper resistance protein database (CuRD). Also, the raw metagenomic sequences were assembled into contigs and analyzed based on Open Reading Frames (ORFs) to identify potentially novel copper resistance genes. Among the different resistance systems for copper detoxification under the high copper stress condition, the Cus system was the most enriched system. The results also indicated that genes encoding multi-copper oxidase played a more important role than those encoding efflux proteins. More significantly, several potentially novel copper resistance ORFs were identified by Pfam search and phylogenic analysis. This study demonstrated a new understanding of microbial-mediated copper resistance under high copper stress using high-throughput shotgun sequencing technique.

Causes, consequences and solutions of phylogenetic incongruence

Phylogenetic analysis is used to recover the evolutionary history of species, genes or proteins. Understanding phylogenetic relationships between organisms is a prerequisite of almost any evolutionary study, as contemporary species all share a common history through their ancestry. Moreover, it is important because of its wide applications that include understanding genome organization, epidemiological investigations, predicting protein functions, and deciding the genes to be analyzed in comparative studies. Despite immense progress in recent years, phylogenetic reconstruction involves many challenges that create uncertainty with respect to the true evolutionary relationships of the species or genes analyzed. One of the most notable difficulties is the widespread occurrence of incongruence among methods and also among individual genes or different genomic regions. Presence of widespread incongruence inhibits successful revealing of evolutionary relationships and applications of phylogenetic analysis. In this article, I concisely review the effect of various factors that cause incongruence in molecular phylogenies, the advances in the field that resolved some factors, and explore unresolved factors that cause incongruence along with possible ways for tackling them.

Dynamics of soil bacterial communities in response to repeated application of manure containing sulfadiazine

Large amounts of manure have been applied to arable soils as fertilizer worldwide. Manure is often contaminated with veterinary antibiotics which enter the soil together with antibiotic resistant bacteria. However, little information is available regarding the main responders of bacterial communities in soil affected by repeated inputs of antibiotics via manure. In this study, a microcosm experiment was performed with two concentrations of the antibiotic sulfadiazine (SDZ) which were applied together with manure at three different time points over a period of 133 days. Samples were taken 3 and 60 days after each manure application. The effects of SDZ on soil bacterial communities were explored by barcoded pyrosequencing of 16S rRNA gene fragments amplified from total community DNA. Samples with high concentration of SDZ were analyzed on day 193 only. Repeated inputs of SDZ, especially at a high concentration, caused pronounced changes in bacterial community compositions. By comparison with the initial soil, we could observe an increase of the disturbance and a decrease of the stability of soil bacterial communities as a result of SDZ manure application compared to the manure treatment without SDZ. The number of taxa significantly affected by the presence of SDZ increased with the times of manure application and was highest during the treatment with high SDZ-concentration. Numerous taxa, known to harbor also human pathogens, such as Devosia, Shinella, Stenotrophomonas, Clostridium, Peptostreptococcus, Leifsonia, Gemmatimonas, were enriched in the soil when SDZ was present while the abundance of bacteria which typically contribute to high soil quality belonging to the genera Pseudomonas and Lysobacter, Hydrogenophaga, and Adhaeribacter decreased in response to the repeated application of manure and SDZ.

Phylogeny and evolution of RNA structure

Darwin's conviction that all living beings on Earth are related and the graph of relatedness is tree-shaped has been essentially confirmed by phylogenetic reconstruction first from morphology and later from data obtained by molecular sequencing. Limitations of the phylogenetic tree concept were recognized as more and more sequence information became available. The other path-breaking idea of Darwin, natural selection of fitter variants in populations, is cast into simple mathematical form and extended to mutation-selection dynamics. In this form the theory is directly applicable to RNA evolution in vitro and to virus evolution. Phylogeny and population dynamics of RNA provide complementary insights into evolution and the interplay between the two concepts will be pursued throughout this chapter. The two strategies for understanding evolution are ultimately related through the central paradigm of structural biology: sequence ⇒ structure ⇒ function. We elaborate on the state of the art in modeling both phylogeny and evolution of RNA driven by reproduction and mutation. Thereby the focus will be laid on models for phylogenetic sequence evolution as well as evolution and design of RNA structures with selected examples and notes on simulation methods. In the perspectives an attempt is made to combine molecular structure, population dynamics, and phylogeny in modeling evolution.

Horizontal gene transfer in the acquisition of novel traits by metazoans

Horizontal gene transfer is accepted as an important evolutionary force modulating the evolution of prokaryote genomes. However, it is thought that horizontal gene transfer plays only a minor role in metazoan evolution. In this paper, I critically review the rising evidence on horizontally transferred genes and on the acquisition of novel traits in metazoans. In particular, I discuss suspected examples in sponges, cnidarians, rotifers, nematodes, molluscs and arthropods which suggest that horizontal gene transfer in metazoans is not simply a curiosity. In addition, I stress the scarcity of studies in vertebrates and other animal groups and the importance of forthcoming studies to understand the importance and extent of horizontal gene transfer in animals.

Intestinal microbiome of poultry and its interaction with host and diet

The gastrointestinal (GI) tract of poultry is densely populated with microorganisms which closely and intensively interact with the host and ingested feed. The gut microbiome benefits the host by providing nutrients from otherwise poorly utilized dietary substrates and modulating the development and function of the digestive and immune system. In return, the host provides a permissive habitat and nutrients for bacterial colonization and growth. Gut microbiome can be affected by diet, and different dietary interventions are used by poultry producers to enhance bird growth and reduce risk of enteric infection by pathogens. There also exist extensive interactions among members of the gut microbiome. A comprehensive understanding of these interactions will help develop new dietary or managerial interventions that can enhance bird growth, maximize host feed utilization, and protect birds from enteric diseases caused by pathogenic bacteria.

Evolution and dynamics of megaplasmids with genome sizes larger than 100 kb in the Bacillus cereus group

BACKGROUND

Plasmids play a crucial role in the evolution of bacterial genomes by mediating horizontal gene transfer. However, the origin and evolution of most plasmids remains unclear, especially for megaplasmids. Strains of the Bacillus cereus group contain up to 13 plasmids with genome sizes ranging from 2 kb to 600 kb, and thus can be used to study plasmid dynamics and evolution.

RESULTS

This work studied the origin and evolution of 31 B. cereus group megaplasmids (>100 kb) focusing on the most conserved regions on plasmids, minireplicons. Sixty-five putative minireplicons were identified and classified to six types on the basis of proteins that are essential for replication. Twenty-nine of the 31 megaplasmids contained two or more minireplicons. Phylogenetic analysis of the protein sequences showed that different minireplicons on the same megaplasmid have different evolutionary histories. Therefore, we speculated that these megaplasmids are the results of fusion of smaller plasmids. All plasmids of a bacterial strain must be compatible. In megaplasmids of the B. cereus group, individual minireplicons of different megaplasmids in the same strain belong to different types or subtypes. Thus, the subtypes of each minireplicon they contain may determine the incompatibilities of megaplasmids. A broader analysis of all 1285 bacterial plasmids with putative known minireplicons whose complete genome sequences were available from GenBank revealed that 34% (443 plasmids) of the plasmids have two or more minireplicons. This indicates that plasmid fusion events are general among bacterial plasmids.

CONCLUSIONS

Megaplasmids of B. cereus group are fusion of smaller plasmids, and the fusion of plasmids likely occurs frequently in the B. cereus group and in other bacterial taxa. Plasmid fusion may be one of the major mechanisms for formation of novel megaplasmids in the evolution of bacteria.

Evolution of bacterial-like phosphoprotein phosphatases in photosynthetic eukaryotes features ancestral mitochondrial or archaeal origin and possible lateral gene transfer

Protein phosphorylation is a reversible regulatory process catalyzed by the opposing reactions of protein kinases and phosphatases, which are central to the proper functioning of the cell. Dysfunction of members in either the protein kinase or phosphatase family can have wide-ranging deleterious effects in both metazoans and plants alike. Previously, three bacterial-like phosphoprotein phosphatase classes were uncovered in eukaryotes and named according to the bacterial sequences with which they have the greatest similarity: Shewanella-like (SLP), Rhizobiales-like (RLPH), and ApaH-like (ALPH) phosphatases. Utilizing the wealth of data resulting from recently sequenced complete eukaryotic genomes, we conducted database searching by hidden Markov models, multiple sequence alignment, and phylogenetic tree inference with Bayesian and maximum likelihood methods to elucidate the pattern of evolution of eukaryotic bacterial-like phosphoprotein phosphatase sequences, which are predominantly distributed in photosynthetic eukaryotes. We uncovered a pattern of ancestral mitochondrial (SLP and RLPH) or archaeal (ALPH) gene entry into eukaryotes, supplemented by possible instances of lateral gene transfer between bacteria and eukaryotes. In addition to the previously known green algal and plant SLP1 and SLP2 protein forms, a more ancestral third form (SLP3) was found in green algae. Data from in silico subcellular localization predictions revealed class-specific differences in plants likely to result in distinct functions, and for SLP sequences, distinctive and possibly functionally significant differences between plants and nonphotosynthetic eukaryotes. Conserved carboxyl-terminal sequence motifs with class-specific patterns of residue substitutions, most prominent in photosynthetic organisms, raise the possibility of complex interactions with regulatory proteins.

Computational approaches to species phylogeny inference and gene tree reconciliation

An intricate relation exists between gene trees and species phylogenies, due to evolutionary processes that act on the genes within and across the branches of the species phylogeny. From an analytical perspective, gene trees serve as character states for inferring accurate species phylogenies, and species phylogenies serve as a backdrop against which gene trees are contrasted for elucidating evolutionary processes and parameters. In a 1997 paper, Maddison discussed this relation, reviewed the signatures left by three major evolutionary processes on the gene trees, and surveyed parsimony and likelihood criteria for utilizing these signatures to elucidate computationally this relation. Here, I review progress that has been made in developing computational methods for analyses under these two criteria, and survey remaining challenges.

Parallel evolution of genome structure and transcriptional landscape in the Epsilonproteobacteria

Background

Gene reshuffling, point mutations and horizontal gene transfer contribute to bacterial genome variation, but require the genome to rewire its transcriptional circuitry to ensure that inserted, mutated or reshuffled genes are transcribed at appropriate levels. The genomes of Epsilonproteobacteria display very low synteny, due to high levels of reshuffling and reorganisation of gene order, but still share a significant number of gene orthologs allowing comparison. Here we present the primary transcriptome of the pathogenic Epsilonproteobacterium Campylobacter jejuni, and have used this for comparative and predictive transcriptomics in the Epsilonproteobacteria.

Results

Differential RNA-sequencing using 454 sequencing technology was used to determine the primary transcriptome of C. jejuni NCTC 11168, which consists of 992 transcription start sites (TSS), which included 29 putative non-coding and stable RNAs, 266 intragenic (internal) TSS, and 206 antisense TSS. Several previously unknown features were identified in the C. jejuni transcriptional landscape, like leaderless mRNAs and potential leader peptides upstream of amino acid biosynthesis genes. A cross-species comparison of the primary transcriptomes of C. jejuni and the related Epsilonproteobacterium Helicobacter pylori highlighted a lack of conservation of operon organisation, position of intragenic and antisense promoters or leaderless mRNAs. Predictive comparisons using 40 other Epsilonproteobacterial genomes suggests that this lack of conservation of transcriptional features is common to all Epsilonproteobacterial genomes, and is associated with the absence of genome synteny in this subdivision of the Proteobacteria.

Conclusions

Both the genomes and transcriptomes of Epsilonproteobacteria are highly variable, both at the genome level by combining and division of multicistronic operons, but also on the gene level by generation or deletion of promoter sequences and 5′ untranslated regions. Regulatory features may have evolved after these species split from a common ancestor, with transcriptome rewiring compensating for changes introduced by genomic reshuffling and horizontal gene transfer.