Adapting to environmental changes using specialized paralogs

Unraveling radiation resistance strategies in two bacterial strains from the high background radiation area of Chavara-Neendakara: A comprehensive whole genome analysis

This paper reports the results of gamma irradiation experiments and whole genome sequencing (WGS) performed on vegetative cells of two radiation resistant bacterial strains, Metabacillus halosaccharovorans (VITHBRA001) and Bacillus paralicheniformis (VITHBRA024) (D10 values 2.32 kGy and 1.42 kGy, respectively), inhabiting the top-ranking high background radiation area (HBRA) of Chavara-Neendakara placer deposit (Kerala, India). The present investigation has been carried out in the context that information on strategies of bacteria having mid-range resistance for gamma radiation is inadequate. WGS, annotation, COG and KEGG analyses and manual curation of genes helped us address the possible pathways involved in the major domains of radiation resistance, involving recombination repair, base excision repair, nucleotide excision repair and mismatch repair, and the antioxidant genes, which the candidate could activate to survive under ionizing radiation. Additionally, with the help of these data, we could compare the candidate strains with that of the extremely radiation resistant model bacterium Deinococccus radiodurans, so as to find the commonalities existing in their strategies of resistance on the one hand, and also the rationale behind the difference in D10, on the other. Genomic analysis of VITHBRA001 and VITHBRA024 has further helped us ascertain the difference in capability of radiation resistance between the two strains. Significantly, the genes such as uvsE (NER), frnE (protein protection), ppk1 and ppx (non-enzymatic metabolite production) and those for carotenoid biosynthesis, are endogenous to VITHBRA001, but absent in VITHBRA024, which could explain the former's better radiation resistance. Further, this is the first-time study performed on any bacterial population inhabiting an HBRA. This study also brings forward the two species whose radiation resistance has not been reported thus far, and add to the knowledge on radiation resistant capabilities of the phylum Firmicutes which are abundantly observed in extreme environment.

Substrate specificity and ecological significance of PstS homologs in phosphorus uptake in marine Synechococcus sp. WH8102

Phosphorus, a vital macronutrient, often limits primary productivity in marine environments. Marine Synechococcus strains, including WH8102, rely on high-affinity phosphate-binding proteins (PstS) to scavenge inorganic phosphate in oligotrophic oceans. However, WH8102 possesses three distinct PstS homologs whose substrate specificity and ecological roles are unclear. The three PstS homologs were heterologously expressed and purified to investigate their substrate specificity and binding kinetics. Our study revealed that all three PstS homologs exhibited a high degree of specificity for phosphate but differed in phosphate binding affinities. Notably, PstS1b displayed nearly 10-fold higher binding affinity (KD = 0.44 µM) compared to PstS1a (KD = 3.3 μM) and PstS2 (KD = 4.3 μM). Structural modeling suggested a single amino acid variation in the binding pocket of PstS1b (threonine instead of serine in PstS1a and PstS2) likely contributed to its higher Pi affinity. Genome context data, together with the protein biophysical data, suggest distinct ecological roles for the three PstS homologs. We propose that PstS1b may be involved in scavenging inorganic phosphorus in oligotrophic conditions and that PstS1a may be involved in transporting recycled phosphate derived from organic phosphate cleavage. The role of PstS2 is less clear, but it may be involved in phosphate uptake when environmental phosphate concentrations are transiently higher. The conservation of three distinct PstS homologs in Synechococcus clade III strains likely reflects distinct adaptations for P acquisition under varying oligotrophic conditions.IMPORTANCEPhosphorus is an essential macronutrient that plays a key role in marine primary productivity and biogeochemistry. However, intense competition for bioavailable phosphorus in the marine environment limits growth and productivity of ecologically important cyanobacteria. In oligotrophic oceans, marine Synechococcus strains, like WH8102, utilize high-affinity phosphate-binding proteins (PstS) to scavenge inorganic phosphate. However, WH8102 possesses three distinct PstS homologs, with unclear substrate specificity and ecological roles, creating a knowledge gap in understanding phosphorus acquisition mechanisms in picocyanobacteria. Through genomic, functional, biophysical, and structural analysis, our study unravels the ecological functions of these homologs. Our findings enhance our understanding of cyanobacterial nutritional uptake strategies and shed light on the crucial role of these conserved nutrient uptake systems in adaptation to specific niches, which ultimately underpins the success of marine Synechococcus across a diverse array of marine ecosystems.

XRE transcription factors conserved in Caulobacter and φCbK modulate adhesin development and phage production

The xenobiotic response element (XRE) family of transcription factors (TFs), which are commonly encoded by bacteria and bacteriophage, regulate diverse features of bacterial cell physiology and impact phage infection dynamics. Through a pangenome analysis of Caulobacter species isolated from soil and aquatic ecosystems, we uncovered an apparent radiation of a paralogous XRE TF gene cluster, several of which have established functions in the regulation of holdfast adhesin development and biofilm formation in C. crescentus. We further discovered related XRE TFs throughout the class Alphaproteobacteria and its phages, including the φCbK Caulophage, suggesting that members of this cluster impact host-phage interactions. Here we show that a closely related group of XRE transcription factors encoded by both C. crescentus and φCbK can physically interact and function to control the transcription of a common gene set, influencing processes including holdfast development and the production of φCbK virions. The φCbK-encoded XRE paralog, tgrL, is highly expressed at the earliest stages of infection and can directly inhibit transcription of host genes including hfiA, a potent holdfast inhibitor, and gafYZ, an activator of prophage-like gene transfer agents (GTAs). XRE proteins encoded from the C. crescentus chromosome also directly repress gafYZ transcription, revealing a functionally redundant set of host regulators that may protect against spurious production of GTA particles and inadvertent cell lysis. Deleting the C. crescentus XRE transcription factors reduced φCbK burst size, while overexpressing these host genes or φCbK tgrL rescued this burst defect. We conclude that this XRE TF gene cluster, shared by C. crescentus and φCbK, plays an important role in adhesion regulation under phage-free conditions, and influences host-phage dynamics during infection.

Radiation impacts gene redundancy and biofilm regulation of cryoconite microbiomes in Northern Hemisphere glaciers

BACKGROUND

Glaciers harbor diverse microorganisms adapted to extreme conditions with high radiation, fluctuating temperature, and low nutrient availability. In glacial ecosystems, cryoconite granules are hotspots of microbial metabolic activity and could influences the biogeochemical cycle on glacier surface. Climate change could influence glacier dynamics by changing regional meteorological factors (e.g., radiation, precipitation, temperature, wind, and evaporation). Moreover, meteorological factors not only influence glacier dynamics but also directly or indirectly influence cryoconite microbiomes. However, the relationship of the meteorological factors and cryoconite microbiome are poorly understood.

RESULTS

Here, we collected 88 metagenomes from 26 glaciers distributed in the Northern Hemisphere with corresponding public meteorological data to reveal the relationship between meteorological factors and variation of cryoconite microbiome. Our results showed significant differences in taxonomic and genomic characteristics between cryoconite generalists and specialists. Additionally, we found that the biogeography of both generalists and specialists was influenced by solar radiation. Specialists with smaller genome size and lower gene redundancy were more abundant under high radiation stress, implying that streamlined genomes are more adapted to high radiation conditions. Network analysis revealed that biofilm regulation is a ubiquitous function in response to radiation stress, and hub genes were associated with the formation and dispersion of biofilms.

CONCLUSION

These findings enhance our understanding of glacier cryoconite microbiome variation on a hemispheric scale and indicate the response mechanisms to radiation stress, which will support forecasts of the ecological consequences of future climate change. Video Abstract.

XRE Transcription Factors Conserved in Caulobacter and φCbK Modulate Adhesin Development and Phage Production

Upon infection, transcriptional shifts in both a host bacterium and its invading phage determine host and viral fitness. The xenobiotic response element (XRE) family of transcription factors (TFs), which are commonly encoded by bacteria and phages, regulate diverse features of bacterial cell physiology and impact phage infection dynamics. Through a pangenome analysis of Caulobacter species isolated from soil and aquatic ecosystems, we uncovered an apparent radiation of a paralogous XRE TF gene cluster, several of which have established functions in the regulation of holdfast adhesin development and biofilm formation in C. crescentus. We further discovered related XRE TFs across the class Alphaproteobacteria and its phages, including the φCbK Caulophage, suggesting that members of this gene cluster impact host-phage interactions. Here we show that that a closely related group of XRE proteins, encoded by both C. crescentus and φCbK, can form heteromeric associations and control the transcription of a common gene set, influencing processes including holdfast development and the production of φCbK virions. The φCbK XRE paralog, tgrL, is highly expressed at the earliest stages of infection and can directly repress transcription of hfiA, a potent holdfast inhibitor, and gafYZ, a transcriptional activator of prophage-like gene transfer agents (GTAs) encoded on the C. crescentus chromosome. XRE proteins encoded from the C. crescentus chromosome also directly repress gafYZ transcription, revealing a functionally redundant set of host regulators that may protect against spurious production of GTA particles and inadvertent cell lysis. Deleting host XRE transcription factors reduced φCbK burst size, while overexpressing these genes or φCbK tgrL rescued this burst defect. We conclude that an XRE TF gene cluster, shared by C. crescentus and φCbK, plays an important role in adhesion regulation under phage-free conditions, and influences host-phage dynamics during infection.

Differences in gene expression patterns between cultured and natural Haloquadratum walsbyi ecotypes

Solar crystallizer ponds are characterized by high population density with a relatively simple community structure in terms of species composition. The microbial community in the solar saltern of Santa Pola (Alicante, Spain), is largely dominated by the hyperhalophilic square archaeon Haloquadratum walsbyi. Here we studied metatranscriptomes retrieved from a crystallizer pond during the winter of 2012 and summer of 2014 and compared Hqr. walsbyi’s transcription patterns with that of the cultured strain Hqr. walsbyi HBSQ001. Significant differences were found between natural and the cultured grown strain in the distribution of transcript levels per gene. This likely reflects the adaptation of the cultured strain to the relative homogeneous growth conditions while the natural species, which is represented by multiple ecotypes, is adapted to heterogeneous environmental conditions and challenges of nutrient competition, viral attack, and other stressors. An important consequence of this study is that expression patterns obtained under artificial cultivation conditions cannot be directly extrapolated to gene expression under natural conditions. Moreover, we found 195 significantly differential expressed genes between the seasons, with 140 genes being higher expressed in winter and mainly encode proteins involved in energy and carbon source acquiring processes, and in stress responses.

Biochemical and genetic characterization comparison of four extradiol dioxygenases in Rhizorhabdus wittichii RW1

Rhizorhabdus (previously Sphingomonas) wittichii RW1 uses a diverse array of aromatic organic compounds as energy and carbon sources, including some extremely recalcitrant compounds such as dibenzo-p-dioxin and dibenzofuran. Extradiol dioxygenases play a key role in the metabolism of dibenzofuran (DBF), dibenzo-p-dioxin (DBD), PCBs, and various other aromatic compounds. In this study, a detailed kinetic analysis of four extradiol dioxygenases identified in R. wittichii RW1 (DbfB, Edo2, Edo3, and Edo4) showed all of them to be typical 2,3dihydroxybiphenyl (DHB) dioxygenases with DHB as preferred substrate (k_cat/K_m values of 0.13-188 (µM ^-1 s^-1)) and only slightly lower activity against trihydroxybiphenyl (THB) whereas monocyclic substrates were, to different extents, poor substrates due to high k_m values. All extradiol dioxygenases analyzed were subject to mechanism-based inactivation by 2,2`,3-trihydroxybiphenylether (THBE) the intermediate of DBD degradation. However, Edo4 was superior as reflected by the relatively high partition ratio and the comparably low efficiency of inactivation. Significant differences were observed with respect to their inactivation by 3-chlorocatechol. The absence of any significant mechanism-based inactivation makes Edo3 a perfect candidate for being recruited for chlorobiphenyl degradation where inactivation of extradiol dioxygenases by this intermediate creates significant metabolic problems. KEY POINTS: • Characterization of additional extradiol dioxygenases encoded by RW1 • Identification of differences in 2,2`,3-trihydroxybiphenylether transformation • Identification of differences in inhibition by 3-chlorocatechol.

The Pseudomonas aeruginosa DksA1 protein is involved in H2O2 tolerance and within-macrophages survival and can be replaced by DksA2

In Gram-negative pathogens, the stringent response regulator DksA controls the expression of hundreds of genes, including virulence-related genes. Interestingly, Pseudomonas aeruginosa has two functional DksA paralogs: DksA1 is constitutively expressed and has a zinc-finger motif, while DksA2 is expressed only under zinc starvation conditions and does not contain zinc. DksA1 stimulates the production of virulence factors in vitro and is required for full pathogenicity in vivo. DksA2 can replace these DksA1 functions. Here, the role of dksA paralogs in P. aeruginosa tolerance to H₂O₂-induced oxidative stress has been investigated. The P. aeruginosa dksA1 dksA2 mutant showed impaired H₂O₂ tolerance in planktonic and biofilm-growing cultures and increased susceptibility to macrophages-mediated killing compared to the wild type. Complementation with either dksA1 or dksA2 genes restored the wild type phenotypes. The DksA-dependent tolerance to oxidative stress involves, at least in part, the positive transcriptional control of both katA and katE catalase-encoding genes. These data support the hypothesis that DksA1 and DksA2 are eco-paralogs with indistinguishable function but optimal activity under different environmental conditions, and highlight their mutual contribution to P. aeruginosa virulence.

Linkages of Various Calcium Sources on Immune Performance, Diarrhea Rate, Intestinal Barrier, and Post-gut Microbial Structure and Function in Piglets

The purpose of this experiment was to investigate the effects of different sources of calcium on immune performance, diarrhea rate, intestinal barrier, and post-intestinal flora structure and function in weaned piglets. A total of 1,000 weaned piglets were randomly assigned to five groups 10 replicate pens per treatment, 20 piglets per pen and fed calcium carbonate, calcium citrate, multiple calcium, and organic trace minerals of different concentrations of acidifier diets. The results of the study showed that the replacement of calcium carbonate with calcium citrate and multiple calcium had almost no significant effect on immune indexes (IL-1β, IL-6, IL-10, TNF-α) of piglets compared with the control group (p > 0.05). The five groups did not show a change in the diarrhea rate and diarrhea index (p > 0.05). The diet containing multiple calcium dramatically decreased the TP compared to the C and L diet (p < 0.05). No significant difference in HDL was noted in the five groups (p > 0.05). However, the concentration of LDL in blood in the multiple calcium group was significantly higher than that in groups L and D (p < 0.05). Moreover, the concentration of Glu in blood in the multiple calcium group was significantly higher than that in group C (p < 0.05). Compared with the control group, calcium citrate plus organic trace minerals diet markedly increased UCG-005 abundance in the colon (p < 0.05). In addition, the relative abundance of Prevotellaceae_NK3B31_group had an upward trend in the colon of the M group compared to the D group (p = 0.070). Meanwhile, calcium citrate plus organic trace minerals diet markedly increased Clostridium_sensu_stricto_1 abundance in the colon (p < 0.05). Metagenomic predictions by PICRUSt suggested that the colonic and fecal microbiota was mainly involved in carbohydrate metabolism, amino acid metabolism, energy metabolism, and metabolism of cofactors and vitamins.

Different trajectories of polyploidization shape the genomic landscape of the Brettanomyces bruxellensis yeast species

Polyploidization events are observed across the tree of life and occur in many fungi, plant, and animal species. During evolution, polyploidy is thought to be an important source of speciation and tumorigenesis. However, the origin of polyploid populations is not always clear, and little is known about the precise nature and structure of their complex genome. Using a long-read sequencing strategy, we sequenced 71 strains from the Brettanomyces bruxellensis yeast species, which is found in anthropized environments (e.g., beer, contaminant of wine, kombucha, and ethanol production) and characterized by several polyploid subpopulations. To reconstruct the polyploid genomes, we phased them by using different strategies and found that each subpopulation had a unique polyploidization history with distinct trajectories. The polyploid genomes contain either genetically closely related (with a genetic divergence <1%) or diverged copies (>3%), indicating auto- as well as allopolyploidization events. These latest events have occurred independently with a specific and unique donor in each of the polyploid subpopulations and exclude the known Brettanomyces sister species as possible donors. Finally, loss of heterozygosity events has shaped the structure of these polyploid genomes and underline their dynamics. Overall, our study highlights the multiplicity of the trajectories leading to polyploid genomes within the same species.

Multiple Copies of flhDC in Paraburkholderia unamae Regulate Flagellar Gene Expression, Motility, and Biofilm Formation

FlhDC is a heterohexameric complex that acts as a master regulator of flagellar biosynthesis genes in numerous bacteria. Previous studies have identified a single flhDC operon encoding this complex. However, we found that two flhDC loci are present throughout Paraburkholderia, and two additional flhC copies are also present in Paraburkholderia unamae. Systematic deletion analysis in P. unamae of the different flhDC copies showed that one of the operons, flhDC1, plays the predominant role, with deletion of its genes resulting in a severe inhibition of motility and biofilm formation. Expression analysis using promoter-lacZ fusions and real-time quantitative PCR support the primary role of flhDC1 in flagellar gene regulation, with flhDC2 a secondary contributor. Phylogenetic analysis shows the presence of the flhDC1 and flhDC2 operons throughout Paraburkholderia. In contrast, Burkholderia and other bacteria only carry the copy syntenous with flhDC2. The variations in impact each copy of flhDC has on downstream processes indicate that regulation of FlhDC in P. unamae, and likely other Paraburkholderia species, is regulated at least in part by the presence of multiple copies of these genes. IMPORTANCE Motility is important in the colonization of plant roots by beneficial and pathogenic bacteria, with flagella playing essential roles in host cell adhesion, entrance, and biofilm formation. Flagellar biosynthesis is energetically expensive. Its complex regulation by the FlhDC master regulator is well studied in peritrichous flagella expressing enterics. We report the unique presence throughout Paraburkholderia of multiple copies of flhDC. In P. unamae, the flhDC1 copy showed higher expression and a greater effect on swim motility, flagellar development, and regulation of downstream genes, than the flhDC2 copy that is syntenous to flhDC in Escherichia coli and pathogenic Burkholderia spp. The flhDC genes have evolved differently in these plant-growth-promoting bacteria, giving an additional layer of complexity in gene regulation by FlhDC.

The two Pseudomonas aeruginosa DksA stringent response proteins are largely interchangeable at the whole transcriptome level and in the control of virulence-related traits

The stringent response regulator DksA plays a key role in Gram negative bacteria adaptation to challenging environments. Intriguingly, the plant and human pathogen Pseudomonas aeruginosa is unique as it expresses two functional DksA paralogs: DksA1 and DksA2. However, the role of DksA2 in P. aeruginosa adaptive strategies has been poorly investigated so far. Here, RNA-Seq analysis and phenotypic assays showed that P. aeruginosa DksA1 and DksA2 proteins are largely interchangeable. Relative to wild type P. aeruginosa, transcription of 1779 genes was altered in a dksA1 dksA2 double mutant, and the wild type expression level of ≥90% of these genes was restored by in trans complementation with either dksA1 or dksA2. Interestingly, the expression of a small sub-set of genes seems to be preferentially or exclusively complemented by either dksA1 or dksA2. In addition, evidence has been provided that the DksA-dependent regulation of virulence genes expression is independent and hierarchically dominant over two major P. aeruginosa regulatory circuits, i.e., quorum sensing and cyclic-di-GMP signalling systems. Our findings support the prominent role of both DksA paralogs in P. aeruginosa environmental adaptation.

New insights into hydroxyectoine synthesis and its transcriptional regulation in the broad-salt growing halophilic bacterium Chromohalobacter salexigens

Elucidating the mechanisms controlling the synthesis of hydroxyectoine is important to design novel genetic engineering strategies for optimizing the production of this biotechnologically relevant compatible solute. The genome of the halophilic bacterium Chromohalobacter salexigens carries two ectoine hydroxylase genes, namely ectD and ectE, whose encoded proteins share the characteristic consensus motif of ectoine hydroxylases but showed only a 51.9% identity between them. In this work, we have shown that ectE encodes a secondary functional ectoine hydroxylase and that the hydroxyectoine synthesis mediated by this enzyme contributes to C.␣salexigens thermoprotection. The evolutionary pattern of EctD and EctE and related proteins suggests that they may have arisen from duplication of an ancestral gene preceding the directional divergence that gave origin to the orders Oceanospirillales and Alteromonadales. Osmoregulated expression of ectD at exponential phase, as well as the thermoregulated expression of ectD at the stationary phase, seemed to be dependent on the general stress factor RpoS. In contrast, expression of ectE was always RpoS-dependent regardless of the growth phase and osmotic or heat stress conditions tested. The data presented here suggest that the AraC-GlxA-like EctZ transcriptional regulator, whose encoding gene lies upstream of ectD, plays a dual function under exponential growth as both a transcriptional activator of osmoregulated ectD expression and a repressor of ectE transcription, privileging the synthesis of the main ectoine hydroxylase EctD. Inactivation of ectZ resulted in a higher amount of the total ectoines pool at the expenses of a higher accumulation of ectoine, with maintenance of the hydroxyectoine levels. In addition to the transcriptional control, our results suggest a strong post-transcriptional regulation of hydroxyectoine synthesis. Data on the accumulation of ectoine and hydroxyectoine in rpoS and ectZ strains pave the way for using these genetic backgrounds for metabolic engineering for hydroxyectoine production.

A Bifan Motif Shaped by ArsR1, ArsR2, and Their Cognate Promoters Frames Arsenic Tolerance of Pseudomonas putida

Prokaryotic tolerance to inorganic arsenic is a widespread trait habitually determined by operons encoding an As (III)-responsive repressor (ArsR), an As (V)-reductase (ArsC), and an As (III)-export pump (ArsB), often accompanied by other complementary genes. Enigmatically, the genomes of many environmental bacteria typically contain two or more copies of this basic genetic device arsRBC. To shed some light on the logic of such apparently unnecessary duplication(s) we have inspected the regulation—together and by separate—of the two ars clusters borne by the soil bacterium Pseudomonas putida strain KT2440, in particular the cross talk between the two repressors ArsR1/ArsR2 and the respective promoters. DNase I footprinting and gel retardation analyses of Pars1 and Pars2 with their matching regulators revealed non-identical binding sequences and interaction patterns for each of the systems. However, in vitro transcription experiments exposed that the repressors could downregulate each other’s promoters, albeit within a different set of parameters. The regulatory frame that emerges from these data corresponds to a particular type of bifan motif where all key interactions have a negative sign. The distinct regulatory architecture that stems from coexistence of various ArsR variants in the same cells could enter an adaptive advantage that favors the maintenance of the two proteins as separate repressors.

Comparative genomics and molecular adaptational analysis of Arthrobacter from Sikkim Himalaya provided insights into its survivability under multiple high-altitude stress

Arthrobacter is a dominant aerobic bacterium under the class Actinobacteria, known for its nutritionally versatile nature and wide prevalence in stressful environments. In the current study representative two strains of Arthrobacter, ERGS1:01 and ERGS4:06, with efficient survivability under high altitude stress conditions were selected for comparative genomic studies with their mesophilic counterparts. Physiological analysis and genome insights supported the survival of these strains under multiple high-altitude stress conditions. Molecular cold-adaptation and substitution analysis of the studied strains supported the incidence of more cold-adapted proteins for functionality at low temperatures. Studied strains preferred amino acids like serine, asparagine, lysine, tryptophan for favoring increased flexibility supporting their broad temperature survivability. To the best of our knowledge, this is the first molecular cold adaptation analysis performed for the genus Arthrobacter and has revealed that 'aromaticity', one of the cold-adaptor indicators, should be carefully considered while evaluating cold adaptation strategies in psychrotrophic/psychrophilic bacteria.

Evolutionary history of dimethylsulfoniopropionate (DMSP) demethylation enzyme DmdA in marine bacteria

Dimethylsulfoniopropionate (DMSP), an osmolyte produced by oceanic phytoplankton and bacteria, is primarily degraded by bacteria belonging to the Roseobacter lineage and other marine Alphaproteobacteria via DMSP-dependent demethylase A protein (DmdA). To date, the evolutionary history of DmdA gene family is unclear. Some studies indicate a common ancestry between DmdA and GcvT gene families and a co-evolution between Roseobacter and the DMSP-producing-phytoplankton around 250 million years ago (Mya). In this work, we analyzed the evolution of DmdA under three possible evolutionary scenarios: (1) a recent common ancestor of DmdA and GcvT, (2) a coevolution between Roseobacter and the DMSP-producing-phytoplankton, and (3) an enzymatic adaptation for utilizing DMSP in marine bacteria prior to Roseobacter origin. Our analyses indicate that DmdA is a new gene family originated from GcvT genes by duplication and functional divergence driven by positive selection before a coevolution between Roseobacter and phytoplankton. Our data suggest that Roseobacter acquired dmdA by horizontal gene transfer prior to an environment with higher DMSP. Here, we propose that the ancestor that carried the DMSP demethylation pathway genes evolved in the Archean, and was exposed to a higher concentration of DMSP in a sulfur-rich atmosphere and anoxic ocean, compared to recent Roseobacter eco-orthologs (orthologs performing the same function under different conditions), which should be adapted to lower concentrations of DMSP.

Hydrogen and Carbon Monoxide-Utilizing Kyrpidia spormannii Species From Pantelleria Island, Italy

Volcanic and geothermal areas are hot and often acidic environments that emit geothermal gasses, including H₂, CO and CO₂. Geothermal gasses mix with air, creating conditions where thermoacidophilic aerobic H₂- and CO-oxidizing microorganisms could thrive. Here, we describe the isolation of two Kyrpidia spormannii strains, which can grow autotrophically by oxidizing H₂ and CO with oxygen. These strains, FAVT5 and COOX1, were isolated from the geothermal soils of the Favara Grande on Pantelleria Island, Italy. Extended physiology studies were performed with K. spormannii FAVT5, and showed that this strain grows optimally at 55°C and pH 5.0. The highest growth rate is obtained using H₂ as energy source (μ_max 0.19 ± 0.02 h^–1, doubling time 3.6 h). K. spormannii FAVT5 can additionally grow on a variety of organic substrates, including some alcohols, volatile fatty acids and amino acids. The genome of each strain encodes for two O₂-tolerant hydrogenases belonging to [NiFe] group 2a hydrogenases and transcriptome studies using K. spormannii FAVT5 showed that both hydrogenases are expressed under H₂ limiting conditions. So far no Firmicutes except K. spormannii FAVT5 have been reported to exhibit a high affinity for H₂, with a K_s of 327 ± 24 nM. The genomes of each strain encode for one putative CO dehydrogenase, belonging to Form II aerobic CO dehydrogenases. The genomic potential and physiological properties of these Kyrpidia strains seem to be quite well adapted to thrive in the harsh environmental volcanic conditions.

ArsH protects Pseudomonas putida from oxidative damage caused by exposure to arsenic

The two As resistance arsRBC operons of Pseudomonas putida KT2440 are followed by a downstream gene called arsH that encodes an NADPH-dependent flavin mononucleotide reductase. In this work, we show that the arsH1 and (to a lesser extent) arsH2 genes of P. putida KT2440 strengthened its tolerance to both inorganic As(V) and As(III) and relieved the oxidative stress undergone by cells exposed to either oxyanion. Furthermore, overexpression of arsH1 and arsH2 endowed P. putida with a high tolerance to the oxidative stress caused by diamide (a drainer of metabolic NADPH) in the absence of any arsenic. To examine whether the activity of ArsH was linked to a direct action on the arsenic compounds tested, arsH1 and arsH2 genes were expressed in Escherichia coli, which has an endogenous arsRBC operon but lacks an arsH ortholog. The resulting clones both deployed a lower production of reactive oxygen species (ROS) when exposed to As salts and had a superior endurance to physiological redox insults. These results suggest that besides the claimed direct action on organoarsenicals, ArsH contributes to relieve toxicity of As species by mediating reduction of ROS produced in vivo upon exposure to the oxyanion, e.g. by generating FMNH2 to fuel ROS-quenching activities.

Genomic insights into the Acidobacteria reveal strategies for their success in terrestrial environments

Members of the phylum Acidobacteria are abundant and ubiquitous across soils. We performed a large-scale comparative genome analysis spanning subdivisions 1, 3, 4, 6, 8 and 23 (n = 24) with the goal to identify features to help explain their prevalence in soils and understand their ecophysiology. Our analysis revealed that bacteriophage integration events along with transposable and mobile elements influenced the structure and plasticity of these genomes. Low- and high-affinity respiratory oxygen reductases were detected in multiple genomes, suggesting the capacity for growing across different oxygen gradients. Among many genomes, the capacity to use a diverse collection of carbohydrates, as well as inorganic and organic nitrogen sources (such as via extracellular peptidases), was detected - both advantageous traits in environments with fluctuating nutrient environments. We also identified multiple soil acidobacteria with the potential to scavenge atmospheric concentrations of H2 , now encompassing mesophilic soil strains within the subdivision 1 and 3, in addition to a previously identified thermophilic strain in subdivision 4. This large-scale acidobacteria genome analysis reveal traits that provide genomic, physiological and metabolic versatility, presumably allowing flexibility and versatility in the challenging and fluctuating soil environment.

The Rapid Evolution of an Ohnolog Contributes to the Ecological Specialization of Incipient Yeast Species

Identifying the molecular changes that lead to ecological specialization during speciation is one of the major goals of molecular evolution. One question that remains to be thoroughly investigated is whether ecological specialization derives strictly from adaptive changes and their associated trade-offs, or from conditionally neutral mutations that accumulate under relaxed selection. We used whole-genome sequencing, genome annotation and computational analyses to identify genes that have rapidly diverged between two incipient species of Saccharomyces paradoxus that occupy different climatic regions along a south-west to north-east gradient. As candidate loci for ecological specialization, we identified genes that show signatures of adaptation and accelerated rates of amino acid substitutions, causing asymmetric evolution between lineages. This set of genes includes a glycyl-tRNA-synthetase, GRS2, which is known to be transcriptionally induced under heat stress in the model and sister species S. cerevisiae. Molecular modelling, expression analysis and fitness assays suggest that the accelerated evolution of this gene in the Northern lineage may be caused by relaxed selection. GRS2 arose during the whole-genome duplication (WGD) that occurred 100 million years ago in the yeast lineage. While its ohnolog GRS1 has been preserved in all post-WGD species, GRS2 has frequently been lost and is evolving rapidly, suggesting that the fate of this ohnolog is still to be resolved. Our results suggest that the asymmetric evolution of GRS2 between the two incipient S. paradoxus species contributes to their restricted climatic distributions and thus that ecological specialization derives at least partly from relaxed selection rather than a molecular trade-off resulting from adaptive evolution.

Comprehensive Molecular Characterization of Bacterial Communities in Feces of Pet Birds Using 16S Marker Sequencing

Birds and other animals live and evolve in close contact with millions of microorganisms (microbiota). While the avian microbiota has been well characterized in domestic poultry, the microbiota of other bird species has been less investigated. The aim of this study was to describe the fecal bacterial communities of pet birds. Pooled fecal samples from 22 flocks representing over 150 individual birds of three different species (Melopsittacus undulatus or budgerigars, Nymphicus hollandicus or cockatiels, and Serinus canaria or domestic canaries) were used for analysis using the 16S rRNA gene sequencing in the MiSeq platform (Illumina). Firmicutes was the most abundant phylum (median 88.4 %; range 12.9-98.4 %) followed by other low-abundant phyla such as Proteobacteria (median 2.3 %; 0.1-85.3 %) and Actinobacteria (median 1.7 %; 0-18.3 %). Lactobacillaceae (mostly Lactobacillus spp.) was the most abundant family (median 78.1 %; 1.4-97.5 %), especially in budgerigars and canaries, and it deserves attention because of the ascribed beneficial properties of lactic acid bacteria. Importantly, feces from birds contain intestinal, urinary, and reproductive-associated microbiota thus posing a serious problem to study one anatomical region at a time. Other groups of interest include the family Clostridiaceae that showed very low abundance (overall median <0.1 %) with the exception of two samples from cockatiels (14 and 45.9 %) and one sample from budgerigars (19.9 %). Analysis of UniFrac metrics showed that overall, the microbial communities from the 22 flocks tended to cluster together for each bird species, meaning each species shed distinctive bacterial communities in feces. This descriptive analysis provides insight into the fecal microbiota of pet birds.

Arsenic resistance strategy in Pantoea sp. IMH: Organization, function and evolution of ars genes

Pantoea sp. IMH is the only bacterium found in genus Pantoea with a high As resistance capacity, but its molecular mechanism is unknown. Herein, the organization, function, and evolution of ars genes in IMH are studied starting with analysis of the whole genome. Two ars systems - ars1 (arsR1B1C1H1) and ars2 (arsR2B2C2H2) - with low sequence homology and two arsC-like genes, were found in the IMH genome. Both ars1 and ars2 are involved in the As resistance, where ars1 is the major contributor at 15 °C and ars2 at 30 °C. The difference in the behavior of these two ars systems is attributed to the disparate activities of their arsR promoters at different temperatures. Sequence analysis based on concatenated ArsRBC indicates that ars1 and ars2 clusters may be acquired from Franconibacter helveticus LMG23732 and Serratia marcescens (plasmid R478), respectively, by horizontal gene transfer (HGT). Nevertheless, two arsC-like genes, probably arising from the duplication of arsC, do not contribute to the As resistance. Our results indicate that Pantoea sp. IMH acquired two different As resistance genetic systems by HGT, allowing the colonization of changing ecosystems, and highlighting the flexible adaptation of microorganisms to resist As.

Gene-Family Extension Measures and Correlations

The existence of multiple copies of genes is a well-known phenomenon. A gene family is a set of sufficiently similar genes, formed by gene duplication. In earlier works conducted on a limited number of completely sequenced and annotated genomes it was found that size of gene family and size of genome are positively correlated. Additionally, it was found that several atypical microbes deviated from the observed general trend. In this study, we reexamined these associations on a larger dataset consisting of 1484 prokaryotic genomes and using several ranking approaches. We applied ranking methods in such a way that genomes with lower numbers of gene copies would have lower rank. Until now only simple ranking methods were used; we applied the Kemeny optimal aggregation approach as well. Regression and correlation analysis were utilized in order to accurately quantify and characterize the relationships between measures of paralog indices and genome size. In addition, boxplot analysis was employed as a method for outlier detection. We found that, in general, all paralog indexes positively correlate with an increase of genome size. As expected, different groups of atypical prokaryotic genomes were found for different types of paralog quantities. Mycoplasmataceae and Halobacteria appeared to be among the most interesting candidates for further research of evolution through gene duplication.

Comparative analysis of genome-wide Mlo gene family in Cajanus cajan and Phaseolus vulgaris

The Mlo gene was discovered in barley because the mutant 'mlo' allele conferred broad-spectrum, non-race-specific resistance to powdery mildew caused by Blumeria graminis f. sp. hordei. The Mlo genes also play important roles in growth and development of plants, and in responses to biotic and abiotic stresses. The Mlo gene family has been characterized in several crop species, but only a single legume species, soybean (Glycine max L.), has been investigated so far. The present report describes in silico identification of 18 CcMlo and 20 PvMlo genes in the important legume crops Cajanus cajan (L.) Millsp. and Phaseolus vulgaris L., respectively. In silico analysis of gene organization, protein properties and conserved domains revealed that the C. cajan and P. vulgaris Mlo gene paralogs are more divergent from each other than from their orthologous pairs. The comparative phylogenetic analysis classified CcMlo and PvMlo genes into three major clades. A comparative analysis of CcMlo and PvMlo proteins with the G. max Mlo proteins indicated close association of one CcMlo, one PvMlo with two GmMlo genes, indicating that there was no further expansion of the Mlo gene family after the separation of these species. Thus, most of the diploid species of eudicots might be expected to contain 15-20 Mlo genes. The genes CcMlo12 and 14, and PvMlo11 and 12 are predicted to participate in powdery mildew resistance. If this prediction were verified, these genes could be targeted by TILLING or CRISPR to isolate powdery mildew resistant mutants.

Comparative Genomics between Two Xenorhabdus bovienii Strains Highlights Differential Evolutionary Scenarios within an Entomopathogenic Bacterial Species

Bacteria of the genus Xenorhabdus are symbionts of soil entomopathogenic nematodes of the genus Steinernema. This symbiotic association constitutes an insecticidal complex active against a wide range of insect pests. Within Xenorhabdus bovienii species, the X. bovienii CS03 strain (Xb CS03) is nonvirulent when directly injected into lepidopteran insects, and displays a low virulence when associated with its Steinernema symbiont. The genome of Xb CS03 was sequenced and compared with the genome of a virulent strain, X. bovienii SS-2004 (Xb SS-2004). The genome size and content widely differed between the two strains. Indeed, Xb CS03 had a large genome containing several specific loci involved in the inhibition of competitors, including a few NRPS-PKS loci (nonribosomal peptide synthetases and polyketide synthases) producing antimicrobial molecules. Consistently, Xb CS03 had a greater antimicrobial activity than Xb SS-2004. The Xb CS03 strain contained more pseudogenes than Xb SS-2004. Decay of genes involved in the host invasion and exploitation (toxins, invasins, or extracellular enzymes) was particularly important in Xb CS03. This may provide an explanation for the nonvirulence of the strain when injected into an insect host. We suggest that Xb CS03 and Xb SS-2004 followed divergent evolutionary scenarios to cope with their peculiar life cycle. The fitness strategy of Xb CS03 would involve competitor inhibition, whereas Xb SS-2004 would quickly and efficiently kill the insect host. Hence, Xenorhabdus strains would have widely divergent host exploitation strategies, which impact their genome structure.

Flagella of halophilic archaea: differences in supramolecular organization

Archaeal flagella are similar functionally to bacterial flagella, but structurally they are completely different. Helical archaeal flagellar filaments are formed of protein subunits called flagellins (archaellins). Notwithstanding progress in studies of archaeal flagella achieved in recent years, many problems in this area are still unsolved. In this review, we analyze the formation of these supramolecular structures by the example of flagellar filaments of halophilic archaea. Recent data on the structure of the flagellar filaments demonstrate that their supramolecular organization differs considerably in different haloarchaeal species.

The two paralogue phoN (phosphinothricin acetyl transferase) genes of Pseudomonas putida encode functionally different proteins

Phosphinothricin (PPT) is a non-specific inhibitor of glutamine synthetase that has been employed as herbicide for selection of transgenic plants expressing cognate resistance genes. While the soil bacterium Pseudomonas putida KT2440 has been generally considered PPT-sensitive, inspection of its genome sequence reveals the presence of two highly similar open reading frames (PP_1924 and PP_4846) encoding acetylases with a potential to cause tolerance to the herbicide. To explore this possibility, each of these genes (named phoN1 and phoN2) was separately cloned and their activities examined in vivo and in vitro. Genetic and biochemical evidence indicated that phoN1 encodes a bona fide PPT-acetyl transferase, the expression of which suffices to make P. putida tolerant to high concentrations of the herbicide. In contrast, PhoN2 does not act on PPT but displays instead activity against methionine sulfoximine (MetSox), another glutamine synthetase inhibitor. When the geometry of the substrate-binding site of PhoN1 was grafted with the equivalent residues of the predicted PhoN2 structure, the resulting protein increased significantly MetSox resistance of the expression host concomitantly with the loss of activity on PPT. These observations uncover intricate biochemical and genetic interactions among soil microorganisms and how they can be perturbed by exposure to generic herbicides in soil.

Ail and PagC-related proteins in the entomopathogenic bacteria of Photorhabdus genus

Among pathogenic Enterobacteriaceae, the proteins of the Ail/OmpX/PagC family form a steadily growing family of outer membrane proteins with diverse biological properties, potentially involved in virulence such as human serum resistance, adhesion and entry into eukaryotic culture cells. We studied the proteins Ail/OmpX/PagC in the bacterial Photorhabdus genus. The Photorhabdus bacteria form symbiotic complexes with nematodes of Heterorhabditis species, associations which are pathogenic to insect larvae. Our phylogenetic analysis indicated that in Photorhabdus asymbiotica and Photorhabdus luminescens only Ail and PagC proteins are encoded. The genomic analysis revealed that the Photorhabdus ail and pagC genes were present in a unique copy, except two ail paralogs from P. luminescens. These genes, referred to as ail1_Pl and ail2_Pl, probably resulted from a recent tandem duplication. Surprisingly, only ail1_Pl expression was directly controlled by PhoPQ and low external Mg²⁺ conditions. In P. luminescens, the magnesium-sensing two-component regulatory system PhoPQ regulates the outer membrane barrier and is required for pathogenicity against insects. In order to characterize Ail functions in Photorhabdus, we showed that only ail2_Pl and pagC_Pl had the ability, when expressed into Escherichia coli, to confer resistance to complement in human serum. However no effect in resistance to antimicrobial peptides was found. Thus, the role of Ail and PagC proteins in Photorhabdus life cycle is discussed.

Unusual genome complexity in Lactobacillus salivarius JCM1046

Background

Lactobacillus salivarius strains are increasingly being exploited for their probiotic properties in humans and animals. Dissemination of antibiotic resistance genes among species with food or probiotic-association is undesirable and is often mediated by plasmids or integrative and conjugative elements. L. salivarius strains typically have multireplicon genomes including circular megaplasmids that encode strain-specific traits for intestinal survival and probiotic activity. Linear plasmids are less common in lactobacilli and show a very limited distribution in L. salivarius. Here we present experimental evidence that supports an unusually complex multireplicon genome structure in the porcine isolate L. salivarius JCM1046.

Results

JCM1046 harbours a 1.83 Mb chromosome, and four plasmids which constitute 20% of the genome. In addition to the known 219 kb repA-type megaplasmid pMP1046A, we identified and experimentally validated the topology of three additional replicons, the circular pMP1046B (129 kb), a linear plasmid pLMP1046 (101 kb) and pCTN1046 (33 kb) harbouring a conjugative transposon. pMP1046B harbours both plasmid-associated replication genes and paralogues of chromosomally encoded housekeeping and information-processing related genes, thus qualifying it as a putative chromid. pLMP1046 shares limited sequence homology or gene synteny with other L. salivarius plasmids, and its putative replication-associated protein is homologous to the RepA/E proteins found in the large circular megaplasmids of L. salivarius. Plasmid pCTN1046 harbours a single copy of an integrated conjugative transposon (Tn6224) which appears to be functionally intact and includes the tetracycline resistance gene tetM.

Conclusion

Experimental validation of sequence assemblies and plasmid topology resolved the complex genome architecture of L. salivarius JCM1046. A high-coverage draft genome sequence would not have elucidated the genome complexity in this strain. Given the expanding use of L. salivarius as a probiotic, it is important to determine the genotypic and phenotypic organization of L. salivarius strains. The identification of Tn6224-like elements in this species has implications for strain selection for probiotic applications.

Electronic supplementary material

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

Functional coexistence of twin arsenic resistance systems in Pseudomonas putida KT2440

The genome of the soil bacterium Pseudomonas putida KT2440 bears two virtually identical arsRBCH operons putatively encoding resistance to inorganic arsenic species. Single and double chromosomal deletions in each of these ars clusters of this bacterium were tested for arsenic sensitivity and found that the contribution of each operon to the resistance to the metalloid was not additive, as either cluster sufficed to endow cells with high-level resistance. However, otherwise identical traits linked to each of the ars sites diverged when temperature was decreased. Growth of the various mutants at 15°C (instead of the standard 30°C for P. putida) uncovered that ars2 affords a much higher resistance to As (III) than the ars1 counterpart. Reverse transcription polymerase chain reaction of arsB1 and arsB2 genes as well as lacZ fusions to the Pars1 and Pars2 promoters traced the difference to variations in transcription of the corresponding gene sets at each temperature. Functional redundancy may thus be selected as a stable condition - rather than just as transient state - if it affords one key activity to be expressed under a wider range of physicochemical settings. This seems to provide a straightforward solution to regulatory problems in environmental bacteria that thrive under changing scenarios.

Safety in numbers: multiple occurrences of highly similar homologs among Azotobacter vinelandii carbohydrate metabolism proteins probably confer adaptive benefits

BACKGROUND

Gene duplication and horizontal gene transfer are common processes in bacterial and archaeal genomes, and are generally assumed to result in either diversification or loss of the redundant gene copies. However, a recent analysis of the genome of the soil bacterium Azotobacter vinelandii DJ revealed an abundance of highly similar homologs among carbohydrate metabolism genes. In many cases these multiple genes did not appear to be the result of recent duplications, or to function only as a means of stimulating expression by increasing gene dosage, as the homologs were located in varying functional genetic contexts. Based on these initial findings we here report in-depth bioinformatic analyses focusing specifically on highly similar intra-genome homologs, or synologs, among carbohydrate metabolism genes, as well as an analysis of the general occurrence of very similar synologs in prokaryotes.

RESULTS

Approximately 900 bacterial and archaeal genomes were analysed for the occurrence of synologs, both in general and among carbohydrate metabolism genes specifically. This showed that large numbers of highly similar synologs among carbohydrate metabolism genes are very rare in bacterial and archaeal genomes, and that the A. vinelandii DJ genome contains an unusually large amount of such synologs. The majority of these synologs were found to be non-tandemly organized and localized in varying but metabolically relevant genomic contexts. The same observation was made for other genomes harbouring high levels of such synologs. It was also shown that highly similar synologs generally constitute a very small fraction of the protein-coding genes in prokaryotic genomes. The overall synolog fraction of the A. vinelandii DJ genome was well above the data set average, but not nearly as remarkable as the levels observed when only carbohydrate metabolism synologs were considered.

CONCLUSIONS

Large numbers of highly similar synologs are rare in bacterial and archaeal genomes, both in general and among carbohydrate metabolism genes. However, A. vinelandii and several other soil bacteria harbour large numbers of highly similar carbohydrate metabolism synologs which seem not to result from recent duplication or transfer events. These genes may confer adaptive benefits with respect to certain lifestyles and environmental factors, most likely due to increased regulatory flexibility and/or increased gene dosage.

Expression of the six chromate ion transporter homologues of Burkholderia xenovorans LB400

The chromate ion transporter (CHR) superfamily comprises transporters that confer chromate resistance by extruding toxic chromate ions from cytoplasm. Burkholderia xenovorans strain LB400 has been reported to encode six CHR homologues in its multireplicon genome. We found that strain LB400 displays chromate-inducible resistance to chromate. Susceptibility tests of Escherichia coli strains transformed with cloned B. xenovorans chr genes indicated that the six genes confer chromate resistance, although under different growth conditions, and suggested that expression of chr genes is regulated by sulfate. Expression of chr genes was measured by quantitative reverse transcription-PCR (RT-qPCR) from total RNA of B. xenovorans LB400 grown under different concentrations of sulfate and exposed or not to chromate. The chr homologues displayed distinct expression levels, but showed no significant differences in transcription under the various sulfate concentrations tested, indicating that sulfate does not regulate chr gene expression in B. xenovorans. The chrA2 gene, encoded in the megaplasmid, was the only chr gene whose expression was induced by chromate and it was shown to constitute the chromate-responsive chrBACF operon. These data suggest that this determinant is mainly responsible for the B. xenovorans LB400 chromate resistance phenotype.

Haloarcula marismortui archaellin genes as ecoparalogs

The genome of haloarchaeon Haloarcula marismortui contains two archaellin genes-flaA2 and flaB. Earlier we isolated and characterized two H. marismortui strains in that archaella consisting of FlaA2 archaellin (with a minor FlaB fraction) or of FlaB only. Both the FlaA2 and FlaB strains were motile and produced functional helical archaella. Thus, it may seem that the FlaA2 archaellin is redundant. In this study we investigated the biological roles of archaellin redundancy and demonstrated that FlaA2 archaellin is better adapted to more severe conditions of high temperature/low salinity, while FlaB has an advantage with increasing salinity. We used the thermodynamic data and bioinformatics sequence analysis to demonstrate that archaella formed by FlaA2 are more stable than those formed by FlaB. Our combined data indicate that the monomer FlaA2 archaellin is more flexible and leads to more compact and stable formation of filamentous structures. The difference in response to environmental stress indicates that FlaA2 and FlaB replace each other under different environmental conditions and can be considered as ecoparalogs.

The opportunistic pathogen Propionibacterium acnes: insights into typing, human disease, clonal diversification and CAMP factor evolution

We previously described a Multilocus Sequence Typing (MLST) scheme based on eight genes that facilitates population genetic and evolutionary analysis of P. acnes. While MLST is a portable method for unambiguous typing of bacteria, it is expensive and labour intensive. Against this background, we now describe a refined version of this scheme based on two housekeeping (aroE; guaA) and two putative virulence (tly; camp2) genes (MLST₄) that correctly predicted the phylogroup (IA₁, IA₂, IB, IC, II, III), clonal complex (CC) and sequence type (ST) (novel or described) status for 91% isolates (n = 372) via cross-referencing of the four gene allelic profiles to the full eight gene versions available in the MLST database (http://pubmlst.org/pacnes/). Even in the small number of cases where specific STs were not completely resolved, the MLST₄ method still correctly determined phylogroup and CC membership. Examination of nucleotide changes within all the MLST loci provides evidence that point mutations generate new alleles approximately 1.5 times as frequently as recombination; although the latter still plays an important role in the bacterium's evolution. The secreted/cell-associated ‘virulence’ factors tly and camp2 show no clear evidence of episodic or pervasive positive selection and have diversified at a rate similar to housekeeping loci. The co-evolution of these genes with the core genome might also indicate a role in commensal/normal existence constraining their diversity and preventing their loss from the P. acnes population. The possibility that members of the expanded CAMP factor protein family, including camp2, may have been lost from other propionibacteria, but not P. acnes, would further argue for a possible role in niche/host adaption leading to their retention within the genome. These evolutionary insights may prove important for discussions surrounding camp2 as an immunotherapy target for acne, and the effect such treatments may have on commensal lineages.

Microbial and viral metagenomes of a subtropical freshwater reservoir subject to climatic disturbances

Extreme climatic activities, such as typhoons, are widely known to disrupt our natural environment. In particular, studies have revealed that typhoon-induced perturbations can result in several long-term effects on various ecosystems. In this study, we have conducted a 2-year metagenomic survey to investigate the microbial and viral community dynamics associated with environmental changes and seasonal variations in an enclosed freshwater reservoir subject to episodic typhoons. We found that the microbial community structure and the associated metagenomes continuously changed, where microbial richness increased after typhoon events and decreased during winter. Among the environmental factors that influenced changes in the microbial community, precipitation was considered to be the most significant. Similarly, the viral community regularly showed higher relative abundances and diversity during summer in comparison to winter, with major variations happening in several viral families including Siphoviridae, Myoviridae, Podoviridae and Microviridae. Interestingly, we also found that the precipitation level was associated with the terrestrial viral abundance in the reservoir. In contrast to the dynamic microbial community (L-divergence 0.73 ± 0.25), we found that microbial metabolic profiles were relatively less divergent (L-divergence 0.24 ± 0.04) at the finest metabolic resolution. This study provides for the first time a glimpse at the microbial and viral community dynamics of a subtropical freshwater ecosystem, adding a comprehensive set of new knowledge to aquatic environments.

Comparative genomics of Blattabacterium cuenoti: the frozen legacy of an ancient endosymbiont genome

Many insect species have established long-term symbiotic relationships with intracellular bacteria. Symbiosis with bacteria has provided insects with novel ecological capabilities, which have allowed them colonize previously unexplored niches. Despite its importance to the understanding of the emergence of biological complexity, the evolution of symbiotic relationships remains hitherto a mystery in evolutionary biology. In this study, we contribute to the investigation of the evolutionary leaps enabled by mutualistic symbioses by sequencing the genome of Blattabacterium cuenoti, primary endosymbiont of the omnivorous cockroach Blatta orientalis, and one of the most ancient symbiotic associations. We perform comparative analyses between the Blattabacterium cuenoti genome and that of previously sequenced endosymbionts, namely those from the omnivorous hosts the Blattella germanica (Blattelidae) and Periplaneta americana (Blattidae), and the endosymbionts harbored by two wood-feeding hosts, the subsocial cockroach Cryptocercus punctulatus (Cryptocercidae) and the termite Mastotermes darwiniensis (Termitidae). Our study shows a remarkable evolutionary stasis of this symbiotic system throughout the evolutionary history of cockroaches and the deepest branching termite M. darwiniensis, in terms of not only chromosome architecture but also gene content, as revealed by the striking conservation of the Blattabacterium core genome. Importantly, the architecture of central metabolic network inferred from the endosymbiont genomes was established very early in Blattabacterium evolutionary history and could be an outcome of the essential role played by this endosymbiont in the host's nitrogen economy.

The ecology of bacterial genes and the survival of the new

Much of the observed variation among closely related bacterial genomes is attributable to gains and losses of genes that are acquired horizontally as well as to gene duplications and larger amplifications. The genomic flexibility that results from these mechanisms certainly contributes to the ability of bacteria to survive and adapt in varying environmental challenges. However, the duplicability and transferability of individual genes imply that natural selection should operate, not only at the organismal level, but also at the level of the gene. Genes can be considered semiautonomous entities that possess specific functional niches and evolutionary dynamics. The evolution of bacterial genes should respond both to selective pressures that favor competition, mostly among orthologs or paralogs that may occupy the same functional niches, and cooperation, with the majority of other genes coexisting in a given genome. The relative importance of either type of selection is likely to vary among different types of genes, based on the functional niches they cover and on the tightness of their association with specific organismal lineages. The frequent availability of new functional niches caused by environmental changes and biotic evolution should enable the constant diversification of gene families and the survival of new lineages of genes.

Evolutionary constraints and expression analysis of gene duplications in Rhodobacter sphaeroides 2.4.1

Background

Gene duplication is a major force that contributes to the evolution of new metabolic functions in all organisms. Rhodobacter sphaeroides 2.4.1 is a bacterium that displays a wide degree of metabolic versatility and genome complexity and therefore is a fitting model for the study of gene duplications in bacteria. A comprehensive analysis of 234 duplicate gene-pairs in R. sphaeroides was performed using structural constraint and expression analysis.

Results

The results revealed that most gene-pairs in in-paralogs are maintained under negative selection (ω ≤ 0.3), but the strength of selection differed among in-paralog gene-pairs. Although in-paralogs located on different replicons are maintained under purifying selection, the duplicated genes distributed between the primary chromosome (CI) and the second chromosome (CII) are relatively less selectively constrained than the gene-pairs located within each chromosome. The mRNA expression patterns of duplicate gene-pairs were examined through microarray analysis of this organism grown under seven different growth conditions. Results revealed that ~62% of paralogs have similar expression patterns (cosine ≥ 0.90) over all of these growth conditions, while only ~7% of paralogs are very different in their expression patterns (cosine < 0.50).

Conclusions

The overall findings of the study suggest that only a small proportion of paralogs contribute to the metabolic diversity and the evolution of novel metabolic functions in R. sphaeroides. In addition, the lack of relationships between structural constraints and gene-pair expression suggests that patterns of gene-pair expression are likely associated with conservation or divergence of gene-pair promoter regions and other coregulation mechanisms.

Polyclonality of concurrent natural populations of Alteromonas macleodii

We have analyzed a natural population of the marine bacterium, Alteromonas macleodii, from a single sample of seawater to evaluate the genomic diversity present. We performed full genome sequencing of four isolates and 161 metagenomic fosmid clones, all of which were assigned to A. macleodii by sequence similarity. Out of the four strain genomes, A. macleodii deep ecotype (AltDE1) represented a different genome, whereas AltDE2 and AltDE3 were identical to the previously described AltDE. Although the core genome (~80%) had an average nucleotide identity of 98.51%, both AltDE and AltDE1 contained flexible genomic islands (fGIs), that is, genomic islands present in both genomes in the same genomic context but having different gene content. Some of the fGIs encode cell surface receptors known to be phage recognition targets, such as the O-chain of the lipopolysaccharide, whereas others have genes involved in physiological traits (e.g., nutrient transport, degradation, and metal resistance) denoting microniche specialization. The presence in metagenomic fosmids of genomic fragments differing from the sequenced strain genomes, together with the presence of new fGIs, indicates that there are at least two more A. macleodii clones present. The availability of three or more sequences overlapping the same genomic region also allowed us to estimate the frequency and distribution of recombination events among these different clones, indicating that these clustered near the genomic islands. The results indicate that this natural A. macleodii population has multiple clones with a potential for different phage susceptibility and exploitation of resources, within a seemingly unstructured habitat.

Biological consequences of ancient gene acquisition and duplication in the large genome of Candidatus Solibacter usitatus Ellin6076

Members of the bacterial phylum Acidobacteria are widespread in soils and sediments worldwide, and are abundant in many soils. Acidobacteria are challenging to culture in vitro, and many basic features of their biology and functional roles in the soil have not been determined. Candidatus Solibacter usitatus strain Ellin6076 has a 9.9 Mb genome that is approximately 2–5 times as large as the other sequenced Acidobacteria genomes. Bacterial genome sizes typically range from 0.5 to 10 Mb and are influenced by gene duplication, horizontal gene transfer, gene loss and other evolutionary processes. Our comparative genome analyses indicate that the Ellin6076 large genome has arisen by horizontal gene transfer via ancient bacteriophage and/or plasmid-mediated transduction, and widespread small-scale gene duplications, resulting in an increased number of paralogs. Low amino acid sequence identities among functional group members, and lack of conserved gene order and orientation in regions containing similar groups of paralogs, suggest that most of the paralogs are not the result of recent duplication events. The genome sizes of additional cultured Acidobacteria strains were estimated using pulsed-field gel electrophoresis to determine the prevalence of the large genome trait within the phylum. Members of subdivision 3 had larger genomes than those of subdivision 1, but none were as large as the Ellin6076 genome. The large genome of Ellin6076 may not be typical of the phylum, and encodes traits that could provide a selective metabolic, defensive and regulatory advantage in the soil environment.

Shift in ribonucleotide reductase gene expression in Pseudomonas aeruginosa during infection

The roles of different ribonucleotide reductases (RNRs) in bacterial pathogenesis have not been studied systematically. In this work we analyzed the importance of the different Pseudomonas aeruginosa RNRs in pathogenesis using the Drosophila melanogaster host-pathogen interaction model. P. aeruginosa codes for three different RNRs with different environmental requirements. Class II and III RNR chromosomal mutants exhibited reduced virulence in this model. Translational reporter fusions of RNR gene nrdA, nrdJ, or nrdD to the green fluorescent protein were constructed to measure the expression of each class during the infection process. Analysis of the P. aeruginosa infection by flow cytometry revealed increased expression of nrdJ and nrdD and decreased nrdA expression during the infection process. Expression of each RNR class fits with the pathogenicities of the chromosomal deletion mutants. An extended understanding of the pathogenicity and physiology of P. aeruginosa will be important for the development of novel drugs against infections in cystic fibrosis patients.

Loss of genetic redundancy in reductive genome evolution

Biological systems evolved to be functionally robust in uncertain environments, but also highly adaptable. Such robustness is partly achieved by genetic redundancy, where the failure of a specific component through mutation or environmental challenge can be compensated by duplicate components capable of performing, to a limited extent, the same function. Highly variable environments require very robust systems. Conversely, predictable environments should not place a high selective value on robustness. Here we test this hypothesis by investigating the evolutionary dynamics of genetic redundancy in extremely reduced genomes, found mostly in intracellular parasites and endosymbionts. By combining data analysis with simulations of genome evolution we show that in the extensive gene loss suffered by reduced genomes there is a selective drive to keep the diversity of protein families while sacrificing paralogy. We show that this is not a by-product of the known drivers of genome reduction and that there is very limited convergence to a common core of families, indicating that the repertoire of protein families in reduced genomes is the result of historical contingency and niche-specific adaptations. We propose that our observations reflect a loss of genetic redundancy due to a decreased selection for robustness in a predictable environment.

Adaptation of cells to new environments

The evolutionary success of an organism is a testament to its inherent capacity to keep pace with environmental conditions that change over short and long periods. Mechanisms underlying adaptive processes are being investigated with renewed interest and excitement. This revival is partly fueled by powerful technologies that can probe molecular phenomena at a systems scale. Such studies provide spectacular insight into the mechanisms of adaptation, including rewiring of regulatory networks via natural selection of horizontal gene transfers, gene duplication, deletion, readjustment of kinetic parameters, and myriad other genetic reorganizational events. Here, we discuss advances in prokaryotic systems biology from the perspective of evolutionary principles that have shaped regulatory networks for dynamic adaptation to environmental change.

Ribonucleotide reduction - horizontal transfer of a required function spans all three domains

Background

Ribonucleotide reduction is the only de novo pathway for synthesis of deoxyribonucleotides, the building blocks of DNA. The reaction is catalysed by ribonucleotide reductases (RNRs), an ancient enzyme family comprised of three classes. Each class has distinct operational constraints, and are broadly distributed across organisms from all three domains, though few class I RNRs have been identified in archaeal genomes, and classes II and III likewise appear rare across eukaryotes. In this study, we examine whether this distribution is best explained by presence of all three classes in the Last Universal Common Ancestor (LUCA), or by horizontal gene transfer (HGT) of RNR genes. We also examine to what extent environmental factors may have impacted the distribution of RNR classes.

Results

Our phylogenies show that the Last Eukaryotic Common Ancestor (LECA) possessed a class I RNR, but that the eukaryotic class I enzymes are not directly descended from class I RNRs in Archaea. Instead, our results indicate that archaeal class I RNR genes have been independently transferred from bacteria on two occasions. While LECA possessed a class I RNR, our trees indicate that this is ultimately bacterial in origin. We also find convincing evidence that eukaryotic class I RNR has been transferred to the Bacteroidetes, providing a stunning example of HGT from eukaryotes back to Bacteria. Based on our phylogenies and available genetic and genomic evidence, class II and III RNRs in eukaryotes also appear to have been transferred from Bacteria, with subsequent within-domain transfer between distantly-related eukaryotes. Under the three-domains hypothesis the RNR present in the last common ancestor of Archaea and eukaryotes appears, through a process of elimination, to have been a dimeric class II RNR, though limited sampling of eukaryotes precludes a firm conclusion as the data may be equally well accounted for by HGT.

Conclusions

Horizontal gene transfer has clearly played an important role in the evolution of the RNR repertoire of organisms from all three domains of life. Our results clearly show that class I RNRs have spread to Archaea and eukaryotes via transfers from the bacterial domain, indicating that class I likely evolved in the Bacteria. However, against the backdrop of ongoing transfers, it is harder to establish whether class II or III RNRs were present in the LUCA, despite the fact that ribonucleotide reduction is an essential cellular reaction and was pivotal to the transition from RNA to DNA genomes. Instead, a general pattern of ongoing horizontal transmission emerges wherein environmental and enzyme operational constraints, especially the presence or absence of oxygen, are likely to be major determinants of the RNR repertoire of genomes.

An empirical strategy for characterizing bacterial proteomes across species in the absence of genomic sequences

Global protein identification through current proteomics methods typically depends on the availability of sequenced genomes. In spite of increasingly high throughput sequencing technologies, this information is not available for every microorganism and rarely available for entire microbial communities. Nevertheless, the protein-level homology that exists between related bacteria makes it possible to extract biological information from the proteome of an organism or microbial community by using the genomic sequences of a near neighbor organism. Here, we demonstrate a trans-organism search strategy for determining the extent to which near-neighbor genome sequences can be applied to identify proteins in unsequenced environmental isolates. In proof of concept testing, we found that within a CLUSTAL W distance of 0.089, near-neighbor genomes successfully identified a high percentage of proteins within an organism. Application of this strategy to characterize environmental bacterial isolates lacking sequenced genomes, but having 16S rDNA sequence similarity to Shewanella resulted in the identification of 300-500 proteins in each strain. The majority of identified pathways mapped to core processes, as well as to processes unique to the Shewanellae, in particular to the presence of c-type cytochromes. Examples of core functional categories include energy metabolism, protein and nucleotide synthesis and cofactor biosynthesis, allowing classification of bacteria by observation of conserved processes. Additionally, within these core functionalities, we observed proteins involved in the alternative lactate utilization pathway, recently described in Shewanella.