Lateral gene transfer and the origins of prokaryotic groups

The Transmission and Antibiotic Resistance Variation in a Multiple Drug Resistance Clade of Vibrio cholerae Circulating in Multiple Countries in Asia

Vibrio cholerae has caused massive outbreaks and even trans-continental epidemics. In 2008 and 2010, at least 3 remarkable cholera outbreaks occurred in Hainan, Anhui and Jiangsu provinces of China. To address the possible transmissions and the relationships to the 7^th pandemic strains of those 3 outbreaks, we sequenced the whole genomes of the outbreak isolates and compared with the global isolates from the 7^th pandemic. The three outbreaks in this study were caused by a cluster of V. cholerae in clade 3.B which is parallel to the clade 3.C that was transmitted from Nepal to Haiti and caused an outbreak in 2010. Pan-genome analysis provided additional evolution information on the mobile element and acquired multiple antibiotic resistance genes. We suggested that clade 3.B should be monitored because the multiple antibiotic resistant characteristics of this clade and the ‘amplifier’ function of China in the global transmission of current Cholera pandemic. We also show that dedicated whole genome sequencing analysis provided more information than the previous techniques and should be applied in the disease surveillance networks.

Comparative Genomics of Core Metabolism Genes of Cellulolytic and Non-cellulolytic Clostridium Species

Microbial production of fuels such as ethanol, butanol, hydrogen (H₂), and methane (CH₄) from waste biomass has the potential to provide sustainable energy systems that can displace fossil fuel consumption. Screening for microbial diversity and genome sequencing of a wide-range of microorganisms can identify organisms with natural abilities to synthesize these alternative fuels and/or other biotechnological applications. Clostridium species are the most widely studied strict anaerobes capable of fermentative synthesis of ethanol, butanol, or hydrogen directly from waste biomass. Clostridium termitidis CT1112 is a mesophilic, cellulolytic species capable of direct cellulose fermentation to ethanol and organic acids, with concomitant synthesis of H₂ and CO₂. On the basis of 16S ribosomal RNA (rRNA) and chaperonin 60 (cpn60) gene sequence data, phylogenetic analyses revealed a close relationship between C. termitidis and C. cellobioparum. Comparative bioinformatic analyses of the C. termitidis genome with 18 cellulolytic and 10 non-cellulolytic Clostridium species confirmed this relationship, and further revealed that the majority of core metabolic pathway genes in C. termitidis and C. cellobioparum share more than 90% amino acid sequence identity. The gene loci and corresponding amino acid sequences of the encoded enzymes for each pathway were correlated by percentage identity, higher score (better alignment), and lowest e-value (most significant "hit"). In addition, the function of each enzyme was proposed by conserved domain analysis. In this chapter we discuss the comparative analysis of metabolic pathways involved in synthesis of various useful products by cellulolytic and non-cellulolytic biofuel and solvent producing Clostridium species. This study has generated valuable information concerning the core metabolism genes and pathways of C. termitidis CT1112, which is helpful in developing metabolic engineering strategies to enhance its natural capacity for better industrial applications.

A Phylogenomic Analysis of the Bacterial Phylum Fibrobacteres

The Fibrobacteres has been recognized as a bacterial phylum for over a decade, but little is known about the group beyond its environmental distribution, and characterization of its sole cultured representative genus, Fibrobacter, after which the phylum was named. Based on these incomplete data, it is thought that cellulose hydrolysis, anaerobic metabolism, and lack of motility are unifying features of the phylum. There are also contradicting views as to whether an uncultured sister lineage, candidate phylum TG3, should be included in the Fibrobacteres. Recently, chitin-degrading cultured representatives of TG3 were isolated from a hypersaline soda lake, and the genome of one species, Chitinivibrio alkaliphilus, sequenced and described in detail. Here, we performed a comparative analysis of Fibrobacter succinogenes, C. alkaliphilus and eight near or substantially complete Fibrobacteres/TG3 genomes of environmental populations recovered from termite gut, anaerobic digester, and sheep rumen metagenomes. We propose that TG3 should be amalgamated with the Fibrobacteres phylum based on robust monophyly of the two lineages and shared character traits. Polymer hydrolysis, using a distinctive set of glycoside hydrolases and binding domains, appears to be a prominent feature of members of the Fibrobacteres. Not all members of this phylum are strictly anaerobic as some termite gut Fibrobacteres have respiratory chains adapted to the microaerophilic conditions found in this habitat. Contrary to expectations, flagella-based motility is predicted to be an ancestral and common trait in this phylum and has only recently been lost in F. succinogenes and its relatives based on phylogenetic distribution of flagellar genes. Our findings extend current understanding of the Fibrobacteres and provide an improved basis for further investigation of this phylum.

Respiratory quinones in Archaea: phylogenetic distribution and application as biomarkers in the marine environment

The distribution of respiratory quinone electron carriers among cultivated organisms provides clues on both the taxonomy of their producers and the redox processes these are mediating. Our study of the quinone inventories of 25 archaeal species belonging to the phyla Eury-, Cren- and Thaumarchaeota facilitates their use as chemotaxonomic markers for ecologically important archaeal clades. Saturated and monounsaturated menaquinones with six isoprenoid units forming the alkyl chain may serve as chemotaxonomic markers for Thaumarchaeota. Other diagnostic biomarkers are thiophene-bearing quinones for Sulfolobales and methanophenazines as functional quinone analogues of the Methanosarcinales. The ubiquity of saturated menaquinones in the Archaea in comparison to Bacteria suggests that these compounds may represent an ancestral and diagnostic feature of the Archaea. Overlap between quinone compositions of distinct thermophilic and halophilic archaea and bacteria may indicate lateral gene transfer. The biomarker potential of thaumarchaeal quinones was exemplarily demonstrated on a water column profile of the Black Sea. Both, thaumarchaeal quinones and membrane lipids showed similar distributions with maxima at the chemocline. Quinone distributions indicate that Thaumarchaeota dominate respiratory activity at a narrow interval in the chemocline, while they contribute only 9% to the microbial biomass at this depth, as determined by membrane lipid analysis.

Horizontally transferred genes cluster spatially and metabolically

Background

Genomic uptake of DNA by prokaryotes often encompasses more than a single gene. In many cases, several horizontally transferred genes may be acquired together. Accordingly, we expect that horizontally transferred genes cluster spatially in the genome more often than expected if transfers were independent. Further, genes that depend on each other functionally may be unlikely to have beneficial fitness effects when taken up individually by a foreign genome. Hence, we also expect the co-acquisition of functionally related genes, resulting in the clustering of horizontally transferred genes in functional networks.

Results

Analysing spatial and metabolic clustering of recent horizontal (or lateral) gene transfers among 21 γ-proteobacteria, we confirm both predictions. When comparing two datasets of predicted transfers that differ in their expected false-positive rate, we find that the more stringent dataset shows a stronger enrichment of clustered pairs.

Conclusions

The enrichment of interdependent metabolic genes among predicted transfers supports a biologically significant role of horizontally transferred genes in metabolic adaptation. Our results further suggest that spatial and metabolic clustering may be used as a benchmark for methods that predict recent horizontal gene transfers.

Reviewers

This article was reviewed by Peter Gogarten in collaboration with Luiz Thiberio Rangel, and by Yuri Wolf.

Plant growth-promoting bacteria as inoculants in agricultural soils

Plant-microbe interactions in the rhizosphere are the determinants of plant health, productivity and soil fertility. Plant growth-promoting bacteria (PGPB) are bacteria that can enhance plant growth and protect plants from disease and abiotic stresses through a wide variety of mechanisms; those that establish close associations with plants, such as the endophytes, could be more successful in plant growth promotion. Several important bacterial characteristics, such as biological nitrogen fixation, phosphate solubilization, ACC deaminase activity, and production of siderophores and phytohormones, can be assessed as plant growth promotion (PGP) traits. Bacterial inoculants can contribute to increase agronomic efficiency by reducing production costs and environmental pollution, once the use of chemical fertilizers can be reduced or eliminated if the inoculants are efficient. For bacterial inoculants to obtain success in improving plant growth and productivity, several processes involved can influence the efficiency of inoculation, as for example the exudation by plant roots, the bacterial colonization in the roots, and soil health. This review presents an overview of the importance of soil-plant-microbe interactions to the development of efficient inoculants, once PGPB are extensively studied microorganisms, representing a very diverse group of easily accessible beneficial bacteria.

Evolutionary history influences the salinity preference of bacterial taxa in wetland soils

Salinity is a major driver of bacterial community composition across the globe. Despite growing recognition that different bacterial species are present or active at different salinities, the mechanisms by which salinity structures community composition remain unclear. We tested the hypothesis that these patterns reflect ecological coherence in the salinity preferences of phylogenetic groups using a reciprocal transplant experiment of fresh- and saltwater wetland soils. The salinity of both the origin and host environments affected community composition (16S rRNA gene sequences) and activity (CO2 and CH4 production, and extracellular enzyme activity). These changes in community composition and activity rates were strongly correlated, which suggests the effect of environment on function could be mediated, at least in part, by microbial community composition. Based on their distribution across treatments, each phylotype was categorized as having a salinity preference (freshwater, saltwater, or none) and phylogenetic analyses revealed a significant influence of evolutionary history on these groupings. This finding was corroborated by examining the salinity preferences of high-level taxonomic groups. For instance, we found that the majority of α- and γ-proteobacteria in these wetland soils preferred saltwater, while many β-proteobacteria prefer freshwater. Overall, our results indicate the effect of salinity on bacterial community composition results from phylogenetically-clustered salinity preferences.

Dynamic evolution of Geranium mitochondrial genomes through multiple horizontal and intracellular gene transfers

The exchange of genetic material between cellular organelles through intracellular gene transfer (IGT) or between species by horizontal gene transfer (HGT) has played an important role in plant mitochondrial genome evolution. The mitochondrial genomes of Geraniaceae display a number of unusual phenomena including highly accelerated rates of synonymous substitutions, extensive gene loss and reduction in RNA editing. Mitochondrial DNA sequences assembled for 17 species of Geranium revealed substantial reduction in gene and intron content relative to the ancestor of the Geranium lineage. Comparative analyses of nuclear transcriptome data suggest that a number of these sequences have been functionally relocated to the nucleus via IGT. Evidence for rampant HGT was detected in several Geranium species containing foreign organellar DNA from diverse eudicots, including many transfers from parasitic plants. One lineage has experienced multiple, independent HGT episodes, many of which occurred within the past 5.5 Myr. Both duplicative and recapture HGT were documented in Geranium lineages. The mitochondrial genome of Geranium brycei contains at least four independent HGT tracts that are absent in its nearest relative. Furthermore, G. brycei mitochondria carry two copies of the cox1 gene that differ in intron content, providing insight into contrasting hypotheses on cox1 intron evolution.

A fresh look at the evolution and diversification of photochemical reaction centers

In this review, I reexamine the origin and diversification of photochemical reaction centers based on the known phylogenetic relations of the core subunits, and with the aid of sequence and structural alignments. I show, for example, that the protein folds at the C-terminus of the D1 and D2 subunits of Photosystem II, which are essential for the coordination of the water-oxidizing complex, were already in place in the most ancestral Type II reaction center subunit. I then evaluate the evolution of reaction centers in the context of the rise and expansion of the different groups of bacteria based on recent large-scale phylogenetic analyses. I find that the Heliobacteriaceae family of Firmicutes appears to be the earliest branching of the known groups of phototrophic bacteria; however, the origin of photochemical reaction centers and chlorophyll synthesis cannot be placed in this group. Moreover, it becomes evident that the Acidobacteria and the Proteobacteria shared a more recent common phototrophic ancestor, and this is also likely for the Chloroflexi and the Cyanobacteria. Finally, I argue that the discrepancies among the phylogenies of the reaction center proteins, chlorophyll synthesis enzymes, and the species tree of bacteria are best explained if both types of photochemical reaction centers evolved before the diversification of the known phyla of phototrophic bacteria. The primordial phototrophic ancestor must have had both Type I and Type II reaction centers.

Insights into thermoadaptation and the evolution of mesophily from the bacterial phylum Thermotogae

Thermophiles are extremophiles that grow optimally at temperatures >45 °C. To survive and maintain function of their biological molecules, they have a suite of characteristics not found in organisms that grow at moderate temperature (mesophiles). At the cellular level, thermophiles have mechanisms for maintaining their membranes, nucleic acids, and other cellular structures. At the protein level, each of their proteins remains stable and retains activity at temperatures that would denature their mesophilic homologs. Conversely, cellular structures and proteins from thermophiles may not function optimally at moderate temperatures. These differences between thermophiles and mesophiles presumably present a barrier for evolutionary transitioning between the 2 lifestyles. Therefore, studying closely related thermophiles and mesophiles can help us determine how such lifestyle transitions may happen. The bacterial phylum Thermotogae contains hyperthermophiles, thermophiles, mesophiles, and organisms with temperature ranges wide enough to span both thermophilic and mesophilic temperatures. Genomic, proteomic, and physiological differences noted between other bacterial thermophiles and mesophiles are evident within the Thermotogae. We argue that the Thermotogae is an ideal group of organisms for understanding of the response to fluctuating temperature and of long-term evolutionary adaptation to a different growth temperature range.

Knowledge-based discovery for designing CRISPR-CAS systems against invading mobilomes in thermophiles

Clustered regularly interspaced short palindromic repeats (CRISPRs) are direct features of the prokaryotic genomes involved in resistance to their bacterial viruses and phages. Herein, we have identified CRISPR loci together with CRISPR-associated sequences (CAS) genes to reveal their immunity against genome invaders in the thermophilic archaea and bacteria. Genomic survey of this study implied that genomic distribution of CRISPR-CAS systems was varied from strain to strain, which was determined by the degree of invading mobiloms. Direct repeats found to be equal in some extent in many thermopiles, but their spacers were differed in each strain. Phylogenetic analyses of CAS superfamily revealed that genes cmr, csh, csx11, HD domain, devR were belonged to the subtypes of cas gene family. The members in cas gene family of thermophiles were functionally diverged within closely related genomes and may contribute to develop several defense strategies. Nevertheless, genome dynamics, geological variation and host defense mechanism were contributed to share their molecular functions across the thermophiles. A thermophilic archaean, Thermococcus gammotolerans and thermophilic bacteria, Petrotoga mobilis and Thermotoga lettingae have shown superoperons-like appearance to cluster cas genes, which were typically evolved for their defense pathways. A cmr operon was identified with a specific promoter in a thermophilic archaean, Caldivirga maquilingensis. Overall, we concluded that knowledge-based genomic survey and phylogeny-based functional assignment have suggested for designing a reliable genetic regulatory circuit naturally from CRISPR-CAS systems, acquired defense pathways, to thermophiles in future synthetic biology.

Structure and Evolution of the Archaeal Lipid Synthesis Enzyme sn-Glycerol-1-phosphate Dehydrogenase

One of the most critical events in the origins of cellular life was the development of lipid membranes. Archaea use isoprenoid chains linked via ether bonds to sn-glycerol 1-phosphate (G1P), whereas bacteria and eukaryotes use fatty acids attached via ester bonds to enantiomeric sn-glycerol 3-phosphate. NAD(P)H-dependent G1P dehydrogenase (G1PDH) forms G1P and has been proposed to have played a crucial role in the speciation of the Archaea. We present here, to our knowledge, the first structures of archaeal G1PDH from the hyperthermophilic methanogen Methanocaldococcus jannaschii with bound substrate dihydroxyacetone phosphate, product G1P, NADPH, and Zn(2+) cofactor. We also biochemically characterized the enzyme with respect to pH optimum, cation specificity, and kinetic parameters for dihydroxyacetone phosphate and NAD(P)H. The structures provide key evidence for the reaction mechanism in the stereospecific addition for the NAD(P)H-based pro-R hydrogen transfer and the coordination of the Zn(2+) cofactor during catalysis. Structure-based phylogenetic analyses also provide insight into the origins of G1PDH.

A Post-Genomic View of the Ecophysiology, Catabolism and Biotechnological Relevance of Sulphate-Reducing Prokaryotes

Dissimilatory sulphate reduction is the unifying and defining trait of sulphate-reducing prokaryotes (SRP). In their predominant habitats, sulphate-rich marine sediments, SRP have long been recognized to be major players in the carbon and sulphur cycles. Other, more recently appreciated, ecophysiological roles include activity in the deep biosphere, symbiotic relations, syntrophic associations, human microbiome/health and long-distance electron transfer. SRP include a high diversity of organisms, with large nutritional versatility and broad metabolic capacities, including anaerobic degradation of aromatic compounds and hydrocarbons. Elucidation of novel catabolic capacities as well as progress in the understanding of metabolic and regulatory networks, energy metabolism, evolutionary processes and adaptation to changing environmental conditions has greatly benefited from genomics, functional OMICS approaches and advances in genetic accessibility and biochemical studies. Important biotechnological roles of SRP range from (i) wastewater and off gas treatment, (ii) bioremediation of metals and hydrocarbons and (iii) bioelectrochemistry, to undesired impacts such as (iv) souring in oil reservoirs and other environments, and (v) corrosion of iron and concrete. Here we review recent advances in our understanding of SRPs focusing mainly on works published after 2000. The wealth of publications in this period, covering many diverse areas, is a testimony to the large environmental, biogeochemical and technological relevance of these organisms and how much the field has progressed in these years, although many important questions and applications remain to be explored.

CRISPR Content Correlates with the Pathogenic Potential of Escherichia coli

Guide RNA molecules (crRNA) produced from clustered regularly interspaced short palindromic repeat (CRISPR) arrays, altogether with effector proteins (Cas) encoded by cognate cas (CRISPR associated) genes, mount an interference mechanism (CRISPR-Cas) that limits acquisition of foreign DNA in Bacteria and Archaea. The specificity of this action is provided by the repeat intervening spacer carried in the crRNA, which upon hybridization with complementary sequences enables their degradation by a Cas endonuclease. Moreover, CRISPR arrays are dynamic landscapes that may gain new spacers from infecting elements or lose them for example during genome replication. Thus, the spacer content of a strain determines the diversity of sequences that can be targeted by the corresponding CRISPR-Cas system reflecting its functionality. Most Escherichia coli strains possess either type I-E or I-F CRISPR-Cas systems. To evaluate their impact on the pathogenicity of the species, we inferred the pathotype and pathogenic potential of 126 strains of this and other closely related species and analyzed their repeat content. Our results revealed a negative correlation between the number of I-E CRISPR units in this system and the presence of pathogenicity traits: the median number of repeats was 2.5-fold higher for commensal isolates (with 29.5 units, range 0-53) than for pathogenic ones (12.0, range 0-42). Moreover, the higher the number of virulence factors within a strain, the lower the repeat content. Additionally, pathogenic strains of distinct ecological niches (i.e., intestinal or extraintestinal) differ in repeat counts. Altogether, these findings support an evolutionary connection between CRISPR and pathogenicity in E. coli.

A single exposure of sediment sulphate-reducing bacteria to oxytetracycline concentrations relevant to aquaculture enduringly disturbed their activity, abundance and community structure

AIM

Although feed medicated with antibiotics is widely used in animal production to prevent and treat bacterial infections, the effect of these drugs on nontarget anaerobic bacteria is unknown. We aimed to clarify whether a single exposure of sulphate-reducing bacteria (SRB) from a tilapia pond to oxytetracycline (OTC) concentrations relevant to aquaculture impacts their function, abundance and community structure.

METHODS AND RESULTS

To demonstrate changes in SO4(2-) content, SRB abundance, dsrB copy number and SRB diversity, sediment mesocosms were spiked with 5, 25, 50 and 75 mg OTC kg(-1) and examined for 30 days by means of ion chromatography, qPCR, cultivation and fluorescent in situ hybridization (FISH). On day 3, we measured higher SO4(2-) concentrations (ca. two-fold) and a reduction in dsrB copy numbers of approximately 50% in the treatments compared to the controls. After 30 days, a subtle yet measurable enrichment of bacteria from the order Desulfovibrionales occurred in mesocosms receiving ≥ 50 mg OTC kg(-1), notwithstanding that SRB counts decreased two orders of magnitude. OTC was dynamically and reversibly converted into 4-epioxytetracycline and other related compounds in a dose-dependent manner during the experiment.

CONCLUSIONS

A single exposure to rather high OTC concentrations triggered functional and structural changes in a SRB community that manifested quickly and persisted for a month.

SIGNIFICANCE AND IMPACT OF THE STUDY

This study improves our limited knowledge on the ecotoxicology of antibiotics in anaerobic environments.

Bacterial gene import and mesophilic adaptation in archaea

It is widely believed that the archaeal ancestor was hyperthermophilic, but during archaeal evolution, several lineages - including haloarchaea and their sister methanogens, the Thaumarchaeota, and the uncultured Marine Group II and Marine Group III Euryarchaeota (MGII/III) - independently adapted to lower temperatures. Recent phylogenomic studies suggest that the ancestors of these lineages were recipients of massive horizontal gene transfer from bacteria. Many of the acquired genes, which are often involved in metabolism and cell envelope biogenesis, were convergently acquired by distant mesophilic archaea. In this Opinion article, we explore the intriguing hypothesis that the import of these bacterial genes was crucial for the adaptation of archaea to mesophilic lifestyles.

Fungal metabolic gene clusters-caravans traveling across genomes and environments

Metabolic gene clusters (MGCs), physically co-localized genes participating in the same metabolic pathway, are signature features of fungal genomes. MGCs are most often observed in specialized metabolism, having evolved in individual fungal lineages in response to specific ecological needs, such as the utilization of uncommon nutrients (e.g., galactose and allantoin) or the production of secondary metabolic antimicrobial compounds and virulence factors (e.g., aflatoxin and melanin). A flurry of recent studies has shown that several MGCs, whose functions are often associated with fungal virulence as well as with the evolutionary arms race between fungi and their competitors, have experienced horizontal gene transfer (HGT). In this review, after briefly introducing HGT as a source of gene innovation, we examine the evidence for HGT's involvement on the evolution of MGCs and, more generally of fungal metabolism, enumerate the molecular mechanisms that mediate such transfers and the ecological circumstances that favor them, as well as discuss the types of evidence required for inferring the presence of HGT in MGCs. The currently available examples indicate that transfers of entire MGCs have taken place between closely related fungal species as well as distant ones and that they sometimes involve large chromosomal segments. These results suggest that the HGT-mediated acquisition of novel metabolism is an ongoing and successful ecological strategy for many fungal species.

Revisiting bacterial phylogeny: Natural and experimental evidence for horizontal gene transfer of 16S rRNA

Current methods used for phylogenetic classification of prokaryotes largely rely on the sequences of 16S rRNA genes that are ubiquitously present in the cell. Theoretical basis of this methodology is based on the assumption that 16S rRNA genes are only vertically inherited and are thus indigenous to each species. However, microbial genomic analysis has revealed the existence of prokaryotic species containing two types of rRNA (rrn) operons of seemingly different origins. It has also been reported that some bacteria contain 16S rRNA that are mosaics of sequences from multiple species. This suggests that horizontal gene transfer (HGT) occurred for 16S rRNA genes. In addition, a recent HGT experiment mimicking the natural HGT process has shown that a wide range of foreign 16S rRNA genes can be transferred into Escherichia coli, including those from different phylogenetic classes (with a minimum sequence identity of 80.9% to the Escherichia coli 16S rRNA gene). Thus, in contrast to the complexity hypothesis that states informational genes are rarely horizontally transferred between species, 16S rRNA is occasionally amenable to HGT. Results of the current method for rapid identification and classification of prokaryotes based on the 16S rRNA gene should thus be carefully analyzed and interpreted.

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.

Use of genome sequencing to assess nucleotide structure variation of Staphylococcus aureus strains cultured in spaceflight on Shenzhou-X, under simulated microgravity and on the ground

The extreme environment of space could affect microbial behavior and may increase the risk of infectious disease during spaceflight. However, the molecular genetic changes of methicillin-resistant Staphylococcus aureus (MRSA) in response to the spaceflight environment have not been fully clarified. In the present study, we determined the draft genome sequences for an ancestral S. aureus strain (LCT-SAO) isolated from a clinical sample and three derivative strains, LCT-SAS, LCT-SAM and LCT-SAG, cultured in parallel during the spaceflight Shenzhou-X, under simulated microgravity and on the ground, respectively. To evaluate the impact of short-term spaceflight on the MRSA strains, comparative genomic analysis was implemented. Genome-based mapping of toxin genes and antibiotic resistance genes confirmed that these strains have the conventional pathogenicity and resistance to drugs, as none of the strains showed significant changes in these regions after culturing in the three different environments; this result suggests that spaceflight may not change bacterial virulence or drug resistance. Thirty-nine strain-specific sequence variants (SVs) were identified throughout the genomes, and the three derivatives exhibited almost the same mutation rates. Fifty-nine percent of SVs were located in the intergenic regions of the genomes, indicating that S. aureus may have an extremely robust repair mechanism responsible for recognizing and repairing DNA replication mismatches. It is noteworthy that strain LCT-SAS, cultured in space, presented the most unique SVs (n=9) and shared the fewest SVs with LCT-SAM (n=5) and LCT-SAG (n=4). Furthermore, we identified 10 potential deletion regions and 2 potential insertion regions, with LCT-SAS appearing more fragile than other strains by this measure. These results suggest that the environment of space is inherently complicated, with multiple variables, and cannot be simulated in a simple manner. Our results represent the first analysis of nucleotide structure variation of S. aureus strains in a spaceflight environment and also provide a valuable insight for understanding the mutation strategies of MRSA on earth.

Evolution of the F0F1 ATP synthase complex in light of the patchy distribution of different bioenergetic pathways across prokaryotes

Bacteria and archaea are characterized by an amazing metabolic diversity, which allows them to persist in diverse and often extreme habitats. Apart from oxygenic photosynthesis and oxidative phosphorylation, well-studied processes from chloroplasts and mitochondria of plants and animals, prokaryotes utilize various chemo- or lithotrophic modes, such as anoxygenic photosynthesis, iron oxidation and reduction, sulfate reduction, and methanogenesis. Most bioenergetic pathways have a similar general structure, with an electron transport chain composed of protein complexes acting as electron donors and acceptors, as well as a central cytochrome complex, mobile electron carriers, and an ATP synthase. While each pathway has been studied in considerable detail in isolation, not much is known about their relative evolutionary relationships. Wanting to address how this metabolic diversity evolved, we mapped the distribution of nine bioenergetic modes on a phylogenetic tree based on 16S rRNA sequences from 272 species representing the full diversity of prokaryotic lineages. This highlights the patchy distribution of many pathways across different lineages, and suggests either up to 26 independent origins or 17 horizontal gene transfer events. Next, we used comparative genomics and phylogenetic analysis of all subunits of the F0F1 ATP synthase, common to most bacterial lineages regardless of their bioenergetic mode. Our results indicate an ancient origin of this protein complex, and no clustering based on bioenergetic mode, which suggests that no special modifications are needed for the ATP synthase to work with different electron transport chains. Moreover, examination of the ATP synthase genetic locus indicates various gene rearrangements in the different bacterial lineages, ancient duplications of atpI and of the beta subunit of the F0 subcomplex, as well as more recent stochastic lineage-specific and species-specific duplications of all subunits. We discuss the implications of the overall pattern of conservation and flexibility of the F0F1 ATP synthase genetic locus.

The futalosine pathway played an important role in menaquinone biosynthesis during early prokaryote evolution

Menaquinone (MK) is an important component of the electron-transfer system in prokaryotes. One of its precursors, 1,4-dihydroxy-2-naphthoate, can be synthesized from chorismate by the classical MK pathway. Interestingly, in some bacteria, chorismate can also be converted to 1,4-dihydroxy-6-naphthoate by four enzymes encoded by mqnABCD in an alternative futalosine pathway. In this study, six crucial enzymes belonging to these two independent nonhomologous pathways were identified in the predicted proteomes of prokaryotes representing a broad phylogenetic distribution. Although the classical MK pathway was found in 32.1% of the proteomes, more than twice the proportion containing the futalosine pathway, the latter was found in a broader taxonomic range of organisms (18 of 31 phyla). The prokaryotes equipped with the classical MK pathway were almost all aerobic or facultatively anaerobic, but those with the futalosine pathway were not only aerobic or facultatively anaerobic but also anaerobic. Phylogenies of enzymes of the classical MK pathway indicated that its genes in archaea were probably acquired by an ancient horizontal gene transfer from bacterial donors. Therefore, the organization of the futalosine pathway likely predated that of the classical MK pathway in the evolutionary history of prokaryotes.

Genomic heterogeneity and ecological speciation within one subspecies of Bacillus subtilis

Closely related bacterial genomes usually differ in gene content, suggesting that nearly every strain in nature may be ecologically unique. We have tested this hypothesis by sequencing the genomes of extremely close relatives within a recognized taxon and analyzing the genomes for evidence of ecological distinctness. We compared the genomes of four Death Valley isolates plus the laboratory strain W23, all previously classified as Bacillus subtilis subsp. spizizenii and hypothesized through multilocus analysis to be members of the same ecotype (an ecologically homogeneous population), named putative ecotype 15 (PE15). These strains showed a history of positive selection on amino acid sequences in 38 genes. Each of the strains was under a different regimen of positive selection, suggesting that each strain is ecologically unique and represents a distinct ecological speciation event. The rate of speciation appears to be much faster than can be resolved with multilocus sequencing. Each PE15 strain contained unique genes known to confer a function for bacteria. Remarkably, no unique gene conferred a metabolic system or subsystem function that was not already present in all the PE15 strains sampled. Thus, the origin of ecotypes within this clade shows no evidence of qualitative divergence in the set of resources utilized. Ecotype formation within this clade is consistent with the nanoniche model of bacterial speciation, in which ecotypes use the same set of resources but in different proportions, and genetic cohesion extends beyond a single ecotype to the set of ecotypes utilizing the same resources.

Phylogenetic study of Geitlerinema and Microcystis (Cyanobacteria) using PC-IGS and 16S-23S ITS as markers: investigation of horizontal gene transfer

Selection of genes that have not been horizontally transferred for prokaryote phylogenetic inferences is regarded as a challenging task. The markers internal transcribed spacer of ribosomal genes (16S-23S ITS) and phycocyanin intergenic spacer (PC-IGS), based on the operons of ribosomal and phycocyanin genes respectively, are among the most used markers in cyanobacteria. The region of the ribosomal genes has been considered stable, whereas the phycocyanin operon may have undergone horizontal transfer. To investigate the occurrence of horizontal transfer of PC-IGS, phylogenetic trees of Geitlerinema and Microcystis strains were generated using PC-IGS and 16S-23S ITS and compared. Phylogenetic trees based on the two markers were mostly congruent for Geitlerinema and Microcystis, indicating a common evolutionary history among ribosomal and phycocyanin genes with no evidence for horizontal transfer of PC-IGS. Thus, PC-IGS is a suitable marker, along with 16S-23S ITS for phylogenetic studies of cyanobacteria.

Bacterial origin of a diverse family of UDP-glycosyltransferase genes in the Tetranychus urticae genome

UDP-glycosyltransferases (UGTs) catalyze the conjugation of a variety of small lipophilic molecules with uridine diphosphate (UDP) sugars, altering them into more water-soluble metabolites. Thereby, UGTs play an important role in the detoxification of xenobiotics and in the regulation of endobiotics. Recently, the genome sequence was reported for the two-spotted spider mite, Tetranychus urticae, a polyphagous herbivore damaging a number of agricultural crops. Although various gene families implicated in xenobiotic metabolism have been documented in T. urticae, UGTs so far have not. We identified 80 UGT genes in the T. urticae genome, the largest number of UGT genes in a metazoan species reported so far. Phylogenetic analysis revealed that lineage-specific gene expansions increased the diversity of the T. urticae UGT repertoire. Genomic distribution, intron-exon structure and structural motifs in the T. urticae UGTs were also described. In addition, expression profiling after host-plant shifts and in acaricide resistant lines supported an important role for UGT genes in xenobiotic metabolism. Expanded searches of UGTs in other arachnid species (Subphylum Chelicerata), including a spider, a scorpion, two ticks and two predatory mites, unexpectedly revealed the complete absence of UGT genes. However, a centipede (Subphylum Myriapoda) and a water flea and a crayfish (Subphylum Crustacea) contain UGT genes in their genomes similar to insect UGTs, suggesting that the UGT gene family might have been lost early in the Chelicerata lineage and subsequently re-gained in the tetranychid mites. Sequence similarity of T. urticae UGTs and bacterial UGTs and their phylogenetic reconstruction suggest that spider mites acquired UGT genes from bacteria by horizontal gene transfer. Our findings show a unique evolutionary history of the T. urticae UGT gene family among other arthropods and provide important clues to its functions in relation to detoxification and thereby host adaptation.

The extent of genome flux and its role in the differentiation of bacterial lineages

Horizontal gene transfer (HGT) and gene loss are key processes in bacterial evolution. However, the role of gene gain and loss in the emergence and maintenance of ecologically differentiated bacterial populations remains an open question. Here, we use whole-genome sequence data to quantify gene gain and loss for 27 lineages of the plant-associated bacterium Pseudomonas syringae. We apply an extensive error-control procedure that accounts for errors in draft genome data and greatly improves the accuracy of patterns of gene occurrence among these genomes. We demonstrate a history of extensive genome fluctuation for this species and show that individual lineages could have acquired thousands of genes in the same period in which a 1% amino acid divergence accrues in the core genome. Elucidating the dynamics of genome fluctuation reveals the rapid turnover of gained genes, such that the majority of recently gained genes are quickly lost. Despite high observed rates of fluctuation, a phylogeny inferred from patterns of gene occurrence is similar to a phylogeny based on amino acid replacements within the core genome. Furthermore, the core genome phylogeny suggests that P. syringae should be considered a number of distinct species, with levels of divergence at least equivalent to those between recognized bacterial species. Gained genes are transferred from a variety of sources, reflecting the depth and diversity of the potential gene pool available via HGT. Overall, our results provide further insights into the evolutionary dynamics of genome fluctuation and implicate HGT as a major factor contributing to the diversification of P. syringae lineages.

Horizontal gene transfer and the rock record: comparative genomics of phylogenetically distant bacteria that induce wrinkle structure formation in modern sediments

Wrinkle structures are sedimentary features that are produced primarily through the trapping and binding of siliciclastic sediments by mat-forming micro-organisms. Wrinkle structures and related sedimentary structures in the rock record are commonly interpreted to represent the stabilizing influence of cyanobacteria on sediments because cyanobacteria are known to produce similar textures and structures in modern tidal flat settings. However, other extant bacteria such as filamentous representatives of the family Beggiatoaceae can also interact with sediments to produce sedimentary features that morphologically resemble many of those associated with cyanobacteria-dominated mats. While Beggiatoa spp. and cyanobacteria are metabolically and phylogenetically distant, genomic analyses show that the two groups share hundreds of homologous genes, likely as the result of horizontal gene transfer. The comparative genomics results described here suggest that some horizontally transferred genes may code for phenotypic traits such as filament formation, chemotaxis, and the production of extracellular polymeric substances that potentially underlie the similar biostabilizing influences of these organisms on sediments. We suggest that the ecological utility of certain basic life modes such as the construction of mats and biofilms, coupled with the lateral mobility of genes in the microbial world, introduces an element of uncertainty into the inference of specific phylogenetic origins from gross morphological features preserved in the ancient rock record.

Bacterial genome instability

Bacterial genomes are remarkably stable from one generation to the next but are plastic on an evolutionary time scale, substantially shaped by horizontal gene transfer, genome rearrangement, and the activities of mobile DNA elements. This implies the existence of a delicate balance between the maintenance of genome stability and the tolerance of genome instability. In this review, we describe the specialized genetic elements and the endogenous processes that contribute to genome instability. We then discuss the consequences of genome instability at the physiological level, where cells have harnessed instability to mediate phase and antigenic variation, and at the evolutionary level, where horizontal gene transfer has played an important role. Indeed, this ability to share DNA sequences has played a major part in the evolution of life on Earth. The evolutionary plasticity of bacterial genomes, coupled with the vast numbers of bacteria on the planet, substantially limits our ability to control disease.

Constraint and opportunity in genome innovation

The development of rigorous molecular taxonomy pioneered by Carl Woese has freed evolution science to explore numerous cellular activities that lead to genome change in evolution. These activities include symbiogenesis, inter- and intracellular horizontal DNA transfer, incorporation of DNA from infectious agents, and natural genetic engineering, especially the activity of mobile elements. This article reviews documented examples of all these processes and proposes experiments to extend our understanding of cell-mediated genome change.

Composition of the DNA-uptake complex of Vibrio cholerae

Natural competence for transformation is a developmental program that allows certain bacteria to take up free extracellular DNA from the environment and integrate this DNA into their genome. Thereby, natural transformation acts as mode of horizontal gene transfer and impacts bacterial evolution. The number of genes induced upon competence induction varies significantly between organisms. However, all of the naturally competent bacteria possess competence genes that encode so-called DNA-uptake machineries. Some components of these multi-protein complexes resemble subunits of type IV pili and type II secretion systems. However, knowledge on the mechanistic aspects of such DNA-uptake complexes is still very limited. Here, we discuss some new findings regarding the DNA-uptake machinery of the naturally transformable human pathogen Vibrio cholerae. The potential of this organism to initiate the competence program was discovered less than a decade ago. However, recent studies have provided new insight into both the regulatory pathways of competence induction and into the DNA uptake dynamics.

Common Virulence Factors and Tissue Targets of Entomopathogenic Bacteria for Biological Control of Lepidopteran Pests

This review focuses on common insecticidal virulence factors from entomopathogenic bacteria with special emphasis on two insect pathogenic bacteria Photorhabdus (Proteobacteria: Enterobacteriaceae) and Bacillus (Firmicutes: Bacillaceae). Insect pathogenic bacteria of diverse taxonomic groups and phylogenetic origin have been shown to have striking similarities in the virulence factors they produce. It has been suggested that the detection of phage elements surrounding toxin genes, horizontal and lateral gene transfer events, and plasmid shuffling occurrences may be some of the reasons that virulence factor genes have so many analogs throughout the bacterial kingdom. Comparison of virulence factors of Photorhabdus, and Bacillus, two bacteria with dissimilar life styles opens the possibility of re-examining newly discovered toxins for novel tissue targets. For example, nematodes residing in the hemolymph may release bacteria with virulence factors targeting neurons or neuromuscular junctions. The first section of this review focuses on toxins and their context in agriculture. The second describes the mode of action of toxins from common entomopathogens and the third draws comparisons between Gram positive and Gram negative bacteria. The fourth section reviews the implications of the nervous system in biocontrol.

Strain-level genomic variation in natural populations of Lebetimonas from an erupting deep-sea volcano

Chemolithoautotrophic Epsilonproteobacteria are ubiquitous in sulfidic, oxygen-poor habitats, including hydrothermal vents, marine oxygen minimum zones, marine sediments and sulfidic caves and have a significant role in cycling carbon, hydrogen, nitrogen and sulfur in these environments. The isolation of diverse strains of Epsilonproteobacteria and the sequencing of their genomes have revealed that this group has the metabolic potential to occupy a wide range of niches, particularly at dynamic deep-sea hydrothermal vents. We expand on this body of work by examining the population genomics of six strains of Lebetimonas, a vent-endemic, thermophilic, hydrogen-oxidizing Epsilonproteobacterium, from a single seamount in the Mariana Arc. Using Lebetimonas as a model for anaerobic, moderately thermophilic organisms in the warm, anoxic subseafloor environment, we show that genomic content is highly conserved and that recombination is limited between closely related strains. The Lebetimonas genomes are shaped by mobile genetic elements and gene loss as well as the acquisition of novel functional genes by horizontal gene transfer, which provide the potential for adaptation and microbial speciation in the deep sea. In addition, these Lebetimonas genomes contain two operons of nitrogenase genes with different evolutionary origins. Lebetimonas expressed nifH during growth with nitrogen gas as the sole nitrogen source, thus providing the first evidence of nitrogen fixation in any Epsilonproteobacteria from deep-sea hydrothermal vents. In this study, we provide a comparative overview of the genomic potential within the Nautiliaceae as well as among more distantly related hydrothermal vent Epsilonproteobacteria to broaden our understanding of microbial adaptation and diversity in the deep sea.

Algal endosymbionts as vectors of horizontal gene transfer in photosynthetic eukaryotes

Photosynthesis in eukaryotes occurs in the plastid, an organelle that is derived from a single cyanobacterial primary endosymbiosis in the common ancestor of the supergroup Plantae (or Archaeplastida) that includes green, red, and glaucophyte algae and plants. However a variety of other phytoplankton such as the chlorophyll c-containing diatoms, dinoflagellates, and haptophytes contain a red alga-derived plastid that traces its origin to secondary or tertiary (eukaryote engulfs eukaryote) endosymbiosis. The hypothesis of Plantae monophyly has only recently been substantiated, however the extent and role of endosymbiotic and horizontal gene transfer (EGT and HGT) in algal genome evolution still remain to be fully understood. What is becoming clear from analysis of complete genome data is that algal gene complements can no longer be considered essentially eukaryotic in provenance; i.e., with the expected addition of several hundred cyanobacterial genes derived from EGT and a similar number derived from the mitochondrial ancestor. For example, we now know that foreign cells such as Chlamydiae and other prokaryotes have made significant contributions to plastid functions in Plantae. Perhaps more surprising is the recent finding of extensive bacterium-derived HGT in the nuclear genome of the unicellular red alga Porphyridium purpureum that does not relate to plastid functions. These non-endosymbiont gene transfers not only shaped the evolutionary history of Plantae but also were propagated via secondary endosymbiosis to a multitude of other phytoplankton. Here we discuss the idea that Plantae (in particular red algae) are one of the major players in eukaryote genome evolution by virtue of their ability to act as "sinks" and "sources" of foreign genes through HGT and endosymbiosis, respectively. This hypothesis recognizes the often under-appreciated Rhodophyta as major sources of genetic novelty among photosynthetic eukaryotes.

Functional genes to assess nitrogen cycling and aromatic hydrocarbon degradation: primers and processing matter

Targeting sequencing to genes involved in key environmental processes, i.e., ecofunctional genes, provides an opportunity to sample nature's gene guilds to greater depth and help link community structure to process-level outcomes. Vastly different approaches have been implemented for sequence processing and, ultimately, for taxonomic placement of these gene reads. The overall quality of next generation sequence analysis of functional genes is dependent on multiple steps and assumptions of unknown diversity. To illustrate current issues surrounding amplicon read processing we provide examples for three ecofunctional gene groups. A combination of in silico, environmental and cultured strain sequences was used to test new primers targeting the dioxin and dibenzofuran degrading genes dxnA1, dbfA1, and carAa. The majority of obtained environmental sequences were classified into novel sequence clusters, illustrating the discovery value of the approach. For the nitrite reductase step in denitrification, the well-known nirK primers exhibited deficiencies in reference database coverage, illustrating the need to refine primer-binding sites and/or to design multiple primers, while nirS primers exhibited bias against five phyla. Amino acid-based OTU clustering of these two N-cycle genes from soil samples yielded only 114 unique nirK and 45 unique nirS genus-level groupings, likely a reflection of constricted primer coverage. Finally, supervised and non-supervised OTU analysis methods were compared using the nifH gene of nitrogen fixation, with generally similar outcomes, but the clustering (non-supervised) method yielded higher diversity estimates and stronger site-based differences. High throughput amplicon sequencing can provide inexpensive and rapid access to nature's related sequences by circumventing the culturing barrier, but each unique gene requires individual considerations in terms of primer design and sequence processing and classification.

Evolutionary diversification and characterization of the eubacterial gene family encoding DXR type II, an alternative isoprenoid biosynthetic enzyme

BACKGROUND

Isoprenoids constitute a vast family of natural compounds performing diverse and essential functions in all domains of life. In most eubacteria, isoprenoids are synthesized through the methylerythritol 4-phosphate (MEP) pathway. The production of MEP is usually catalyzed by deoxyxylulose 5-phosphate reductoisomerase (DXR-I) but a few organisms use an alternative DXR-like enzyme (DXR-II).

RESULTS

Searches through 1498 bacterial complete proteomes detected 130 sequences with similarity to DXR-II. Phylogenetic analysis identified three well-resolved clades: the DXR-II family (clustering 53 sequences including eleven experimentally verified as functional enzymes able to produce MEP), and two previously uncharacterized NAD(P)-dependent oxidoreductase families (designated DLO1 and DLO2 for DXR-II-like oxidoreductases 1 and 2). Our analyses identified amino acid changes critical for the acquisition of DXR-II biochemical function through type-I functional divergence, two of them mapping onto key residues for DXR-II activity. DXR-II showed a markedly discontinuous distribution, which was verified at several levels: taxonomic (being predominantly found in Alphaproteobacteria and Firmicutes), metabolic (being mostly found in bacteria with complete functional MEP pathways with or without DXR-I), and phenotypic (as no biological/phenotypic property was found to be preferentially distributed among DXR-II-containing strains, apart from pathogenicity in animals). By performing a thorough comparative sequence analysis of GC content, 3:1 dinucleotide frequencies, codon usage and codon adaptation indexes (CAI) between DXR-II sequences and their corresponding genomes, we examined the role of horizontal gene transfer (HGT), as opposed to an scenario of massive gene loss, in the evolutionary origin and diversification of the DXR-II subfamily in bacteria.

CONCLUSIONS

Our analyses support a single origin of the DXR-II family through functional divergence, in which constitutes an exceptional model of acquisition and maintenance of redundant gene functions between non-homologous genes as a result of convergent evolution. Subsequently, although old episodic events of HGT could not be excluded, the results supported a prevalent role of gene loss in explaining the distribution of DXR-II in specific pathogenic eubacteria. Our results highlight the importance of the functional characterization of evolutionary shortcuts in isoprenoid biosynthesis for screening specific antibacterial drugs and for regulating the production of isoprenoids of human interest.

Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences

Profiling phylogenetic marker genes, such as the 16S rRNA gene, is a key tool for studies of microbial communities but does not provide direct evidence of a community’s functional capabilities. Here we describe PICRUSt (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States), a computational approach to predict the functional composition of a metagenome using marker gene data and a database of reference genomes. PICRUSt uses an extended ancestral-state reconstruction algorithm to predict which gene families are present and then combines gene families to estimate the composite metagenome. Using 16S information, PICRUSt recaptures key findings from the Human Microbiome Project and accurately predicts the abundance of gene families in host-associated and environmental communities, with quantifiable uncertainty. Our results demonstrate that phylogeny and function are sufficiently linked that this ‘predictive metagenomic’ approach should provide useful insights into the thousands of uncultivated microbial communities for which only marker gene surveys are currently available.

Integrating metagenomic and amplicon databases to resolve the phylogenetic and ecological diversity of the Chlamydiae

In the era of metagenomics and amplicon sequencing, comprehensive analyses of available sequence data remain a challenge. Here we describe an approach exploiting metagenomic and amplicon data sets from public databases to elucidate phylogenetic diversity of defined microbial taxa. We investigated the phylum Chlamydiae whose known members are obligate intracellular bacteria that represent important pathogens of humans and animals, as well as symbionts of protists. Despite their medical relevance, our knowledge about chlamydial diversity is still scarce. Most of the nine known families are represented by only a few isolates, while previous clone library-based surveys suggested the existence of yet uncharacterized members of this phylum. Here we identified more than 22,000 high quality, non-redundant chlamydial 16S rRNA gene sequences in diverse databases, as well as 1900 putative chlamydial protein-encoding genes. Even when applying the most conservative approach, clustering of chlamydial 16S rRNA gene sequences into operational taxonomic units revealed an unexpectedly high species, genus and family-level diversity within the Chlamydiae, including 181 putative families. These in silico findings were verified experimentally in one Antarctic sample, which contained a high diversity of novel Chlamydiae. In our analysis, the Rhabdochlamydiaceae, whose known members infect arthropods, represents the most diverse and species-rich chlamydial family, followed by the protist-associated Parachlamydiaceae, and a putative new family (PCF8) with unknown host specificity. Available information on the origin of metagenomic samples indicated that marine environments contain the majority of the newly discovered chlamydial lineages, highlighting this environment as an important chlamydial reservoir.

Bacteria-human somatic cell lateral gene transfer is enriched in cancer samples

There are 10× more bacterial cells in our bodies from the microbiome than human cells. Viral DNA is known to integrate in the human genome, but the integration of bacterial DNA has not been described. Using publicly available sequence data from the human genome project, the 1000 Genomes Project, and The Cancer Genome Atlas (TCGA), we examined bacterial DNA integration into the human somatic genome. Here we present evidence that bacterial DNA integrates into the human somatic genome through an RNA intermediate, and that such integrations are detected more frequently in (a) tumors than normal samples, (b) RNA than DNA samples, and (c) the mitochondrial genome than the nuclear genome. Hundreds of thousands of paired reads support random integration of Acinetobacter-like DNA in the human mitochondrial genome in acute myeloid leukemia samples. Numerous read pairs across multiple stomach adenocarcinoma samples support specific integration of Pseudomonas-like DNA in the 5′-UTR and 3′-UTR of four proto-oncogenes that are up-regulated in their transcription, consistent with conversion to an oncogene. These data support our hypothesis that bacterial integrations occur in the human somatic genome and may play a role in carcinogenesis. We anticipate that the application of our approach to additional cancer genome projects will lead to the more frequent detection of bacterial DNA integrations in tumors that are in close proximity to the human microbiome.

There are 10× more bacterial cells in the human body than there are human cells that are part of the human microbiome. Many of those bacteria are in constant, intimate contact with human cells. We sought to establish if bacterial cells insert their own DNA into the human genome. Such random mutations could cause disease in the same manner that mutagens like UV rays from the sun or chemicals in cigarettes induce mutations. We detected the integration of bacterial DNA in the human genome more readily in tumors than normal samples. In particular, extensive amounts of DNA with similarity to Acinetobacter DNA were fused to human mitochondrial DNA in acute myeloid leukemia samples. We also identified specific integrations of DNA with similarity to Pseudomonas DNA near the untranslated regulatory regions of four proto-oncogenes. This supports our hypothesis that bacterial integrations occur in the human somatic genome that may potentially play a role in carcinogenesis. Further study in this area may provide new avenues for cancer prevention.

Microdiversity of extracellular enzyme genes among sequenced prokaryotic genomes

Understanding the relationship between prokaryotic traits and phylogeny is important for predicting and modeling ecological processes. Microbial extracellular enzymes have a pivotal role in nutrient cycling and the decomposition of organic matter, yet little is known about the phylogenetic distribution of genes encoding these enzymes. In this study, we analyzed 3058 annotated prokaryotic genomes to determine which taxa have the genetic potential to produce alkaline phosphatase, chitinase and β-N-acetyl-glucosaminidase enzymes. We then evaluated the relationship between the genetic potential for enzyme production and 16S rRNA phylogeny using the consenTRAIT algorithm, which calculated the phylogenetic depth and corresponding 16S rRNA sequence identity of clades of potential enzyme producers. Nearly half (49.2%) of the genomes analyzed were found to be capable of extracellular enzyme production, and these were non-randomly distributed across most prokaryotic phyla. On average, clades of potential enzyme-producing organisms had a maximum phylogenetic depth of 0.008004-0.009780, though individual clades varied broadly in both size and depth. These values correspond to a minimum 16S rRNA sequence identity of 98.04-98.40%. The distribution pattern we found is an indication of microdiversity, the occurrence of ecologically or physiologically distinct populations within phylogenetically related groups. Additionally, we found positive correlations among the genes encoding different extracellular enzymes. Our results suggest that the capacity to produce extracellular enzymes varies at relatively fine-scale phylogenetic resolution. This variation is consistent with other traits that require a small number of genes and provides insight into the relationship between taxonomy and traits that may be useful for predicting ecological function.

Quantifying homologous replacement of loci between haloarchaeal species

In vitro studies of the haloarchaeal genus Haloferax have demonstrated their ability to frequently exchange DNA between species, whereas rates of homologous recombination estimated from natural populations in the genus Halorubrum are high enough to maintain random association of alleles between five loci. To quantify the effects of gene transfer and recombination of commonly held (relaxed core) genes during the evolution of the class Halobacteria (haloarchaea), we reconstructed the history of 21 genomes representing all major groups. Using a novel algorithm and a concatenated ribosomal protein phylogeny as a reference, we created a directed horizontal genetic transfer (HGT) network of contemporary and ancestral genomes. Gene order analysis revealed that 90% of testable HGTs were by direct homologous replacement, rather than nonhomologous integration followed by a loss. Network analysis revealed an inverse log-linear relationship between HGT frequency and ribosomal protein evolutionary distance that is maintained across the deepest divergences in Halobacteria. We use this mathematical relationship to estimate the total transfers and amino acid substitutions delivered by HGTs in each genome, providing a measure of chimerism. For the relaxed core genes of each genome, we conservatively estimate that 11–20% of their evolution occurred in other haloarchaea. Our findings are unexpected, because the transfer and homologous recombination of relaxed core genes between members of the class Halobacteria disrupts the coevolution of genes; however, the generation of new combinations of divergent but functionally related genes may lead to adaptive phenotypes not available through cumulative mutations and recombination within a single population.

Phylogenetic conservatism of functional traits in microorganisms

A central question in biology is how biodiversity influences ecosystem functioning. Underlying this is the relationship between organismal phylogeny and the presence of specific functional traits. The relationship is complicated by gene loss and convergent evolution, resulting in the polyphyletic distribution of many traits. In microorganisms, lateral gene transfer can further distort the linkage between phylogeny and the presence of specific functional traits. To identify the phylogenetic conservation of specific traits in microorganisms, we developed a new phylogenetic metric-consenTRAIT-to estimate the clade depth where organisms share a trait. We then analyzed the distribution of 89 functional traits across a broad range of Bacteria and Archaea using genotypic and phenotypic data. A total of 93% of the traits were significantly non-randomly distributed, which suggested that vertical inheritance was generally important for the phylogenetic dispersion of functional traits in microorganisms. Further, traits in microbes were associated with a continuum of trait depths (τD), ranging from a few deep to many shallow clades (average τD: 0.101-0.0011 rRNA sequence dissimilarity). Next, we demonstrated that the dispersion and the depth of clades that contain a trait is correlated with the trait's complexity. Specifically, complex traits encoded by many genes like photosynthesis and methanogenesis were found in a few deep clusters, whereas the ability to use simple carbon substrates was highly phylogenetically dispersed. On the basis of these results, we propose a framework for predicting the phylogenetic conservatism of functional traits depending on the complexity of the trait. This framework enables predicting how variation in microbial composition may affect microbially-mediated ecosystem processes as well as linking phylogenetic and trait-based patterns of biogeography.

Multiple Rieske/cytb complexes in a single organism

Most organisms contain a single Rieske/cytb complex. This enzyme can be integrated in any respiratory or photosynthetic electron transfer chain that is quinone-based and sufficiently energy rich to allow for the turnover of three enzymes - a quinol reductase, a Rieske/cytb complex and a terminal oxidase. Despite this universal usability of the enzyme a variety of phylogenetically distant organisms have multiple copies thereof and no reason for this redundancy is obvious. In this review we present an overview of the distribution of multiple copies among species and describe their properties from the scarce experimental results, analysis of their amino acid sequences and genomic context. We discuss the predicted redox properties of the Rieske cluster in relation to the nature of the pool quinone. It appears that acidophilic iron-oxidizing bacteria specialized one of their two copies for reverse electron transfer, archaeal Thermoprotei adapted their three copies to the interaction with different oxidases and several, phylogenetically unrelated species imported a second complex with a putative heme ci that may confer some yet to be determined properties to the complex. These hypothesis and all the more the so far completely unexplained cases call for further studies and we put forward a number of suggestions for future research that we hope to be stimulating for the field. This article is part of a Special Issue entitled: Respiratory complex III and related bc complexes.

Being Aquifex aeolicus: Untangling a hyperthermophile's checkered past

Lateral gene transfer (LGT) is an important factor contributing to the evolution of prokaryotic genomes. The Aquificae are a hyperthermophilic bacterial group whose genes show affiliations to many other lineages, including the hyperthermophilic Thermotogae, the Proteobacteria, and the Archaea. Previous phylogenomic analyses focused on Aquifex aeolicus identified Thermotogae and Aquificae either as successive early branches or sisters in a rooted bacterial phylogeny, but many phylogenies and cellular traits have suggested a stronger affiliation with the Epsilonproteobacteria. Different scenarios for the evolution of the Aquificae yield different phylogenetic predictions. Here, we outline these scenarios and consider the fit of the available data, including three sequenced Aquificae genomes, to different sets of predictions. Evidence from phylogenetic profiles and trees suggests that the Epsilonproteobacteria have the strongest affinities with the three Aquificae analyzed. However, this pattern is shown by only a minority of encoded proteins, and the Archaea, many lineages of thermophilic bacteria, and members of genus Clostridium and class Deltaproteobacteria also show strong connections to the Aquificae. The phylogenetic affiliations of different functional subsystems showed strong biases: Most but not all genes implicated in the core translational apparatus tended to group Aquificae with Thermotogae, whereas a wide range of metabolic and cellular processes strongly supported the link between Aquificae and Epsilonproteobacteria. Depending on which sets of genes are privileged, either Thermotogae or Epsilonproteobacteria is the most plausible adjacent lineage to the Aquificae. Both scenarios require massive sharing of genes to explain the history of this enigmatic group, whose history is further complicated by specific affinities of different members of Aquificae to different partner lineages.