Skip to main content
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • Special Issues
    • COVID-19 Special Collection
    • Editor's Picks
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mSystems
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • Special Issues
    • COVID-19 Special Collection
    • Editor's Picks
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Ecological and Evolutionary Science

Evolutionary Timeline and Genomic Plasticity Underlying the Lifestyle Diversity in Rhizobiales

Sishuo Wang, Andrew Meade, Hon-Ming Lam, Haiwei Luo
Sarah M. Allard, Editor
Sishuo Wang
aSchool of Life Sciences and State Key Laboratory of Agrobiotechnology, The Chinese University of Hong Kong, Shatin, Hong Kong SAR
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Sishuo Wang
Andrew Meade
bSchool of Biological Sciences, University of Reading, Whiteknights, Reading, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Hon-Ming Lam
aSchool of Life Sciences and State Key Laboratory of Agrobiotechnology, The Chinese University of Hong Kong, Shatin, Hong Kong SAR
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Haiwei Luo
aSchool of Life Sciences and State Key Laboratory of Agrobiotechnology, The Chinese University of Hong Kong, Shatin, Hong Kong SAR
cShenzhen Research Institute, The Chinese University of Hong Kong, Shenzhen, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Haiwei Luo
Sarah M. Allard
University of California, San Diego
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mSystems.00438-20
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Members of the order Rhizobiales include those capable of nitrogen fixation in nodules as well as pathogens of animals and plants. This lifestyle diversity has important implications for agricultural and medical research. Leveraging large-scale genomic data, we infer that Rhizobiales originated as a free-living ancestor ∼1,500 million years ago (Mya) and that the later emergence of host-associated lifestyles broadly coincided with the rise of their eukaryotic hosts. In particular, the first nodulating lineage arose from either Azorhizobium or Bradyrhizobium 150 to 80 Mya, a time range in general concurrent with the emergence of legumes. The rates of lifestyle transitions are highly variable; nodule association is more likely to be lost than gained, whereas animal association likely represents an evolutionary dead end. We searched for statistical correlations between gene presence and lifestyle and identified genes likely contributing to the transition and adaptation to the same lifestyle in divergent lineages. Among the genes potentially promoting successful transitions to major nodulation lineages, the nod and nif clusters for nodulation and nitrogen fixation, respectively, were repeatedly acquired during each transition; the fix, dct, and phb clusters involved in energy conservation under micro-oxic conditions were present in the nonnodulating ancestors; and the secretion systems were acquired in lineage-specific patterns. Our study data suggest that increased eukaryote diversity drives lifestyle diversification of bacteria and highlight both acquired and preexisting traits facilitating the origin of host association.

IMPORTANCE Bacteria form diverse interactions with eukaryotic hosts. This is well represented by the Rhizobiales, a clade of Alphaproteobacteria strategically important for their large diversity of lifestyles with implications for agricultural and medical research. To investigate their lifestyle evolution, we compiled a comprehensive data set of genomes and lifestyle information for over 1,000 Rhizobiales genomes. We show that the origins of major host-associated lineages in Rhizobiales broadly coincided with the emergences of their host plants/animals, suggesting bacterium-host interactions as a driving force in the evolution of Rhizobiales. We further found that, in addition to gene gains, preexisting traits and recurrent losses of specific genomic traits may have played underrecognized roles in the origin of host-associated lineages, providing clues to genetic engineering of microbial agricultural inoculants and prevention of the emergence of potential plant/animal pathogens.

INTRODUCTION

Bacteria form diverse interactions with their hosts, which, from commensalism to mutualism and from parasitism to pathogenesis, have played an important role in the emergence of complex life forms and many other evolutionary and ecological processes (1–4). The diversity of lifestyles in bacteria is well epitomized by the alphaproteobacterial order Rhizobiales (5–8). A large proportion of Rhizobiales species adapt to host-associated lifestyles, many of which have important agricultural, ecological, and medical implications, making Rhizobiales an ideal system to investigate the evolution of bacterial lifestyles (5, 6). This lifestyle diversity includes rhizobia, which form nodules and establish a symbiotic relationship with legumes; plant pathogens such as the causative agents of root tumor Agrobacterium; and Bartonella and Brucella, two important alphaproteobacterial animal pathogens colonizing the erythrocytes of mammals.

In ancient times, bacteria primarily adapted to a free-living lifestyle occupying aquatic and terrestrial habitats (9). Alongside the geosphere-biosphere interactions that have occurred since the Cambrian explosion, the past 500 million years (My) have witnessed a massive diversification of animals and plants (10). This might create many opportunities for bacteria to interact with diverse eukaryotes and might facilitate the radiation of bacteria (5, 11). However, to our knowledge, the time frames of the coevolution between Rhizobiales and their hosts have never been determined. This leaves open two important issues: (i) establishing from which lineages the first nodulating organisms and animal pathogens originated and (ii) determining how the establishment of successful bacterial symbiosis is driven by their coevolution with eukaryotes.

Adaptations to diverse lifestyles suggest high genomic plasticity of Rhizobiales (12, 13). Intriguingly, species of the same lifestyle may not form monophyletic groups in Rhizobiales (6, 14). Recurrent gains and losses of genes important to specific lifestyles may make independent lifestyle transitions occur (6, 15, 16). Despite extensive genetic and physiological studies, our knowledge of these genes, however, is still limited. A genome-wide computational identification of the genes associated with each lifestyle which circumvents the low-throughput limitations of wet-lab experiments is urgently needed. By tracing the evolution of lifestyle-correlated genes during lifestyle transitions, we can better understand how rhizobia and other host-associated bacteria evolve step by step at the genomic level.

Despite the many important issues raised regarding lifestyle diversity in Rhizobiales, genomics studies have mainly focused on organisms with a given lifestyle (13, 17, 18), with the transitions between lifestyles awaiting further exploration. In the present study, we compiled comprehensive Rhizobiales genomics data sets, reconstructed a timeline of the evolutionary origins of distinct host-associated lifestyles, linked them to the evolutionary time of their hosts, and explored the functional changes underlying the transitions between lifestyles. We highlighted the roles of increased diversity of eukaryotes and genomic plasticity in the evolution of host-associated bacteria in Rhizobiales.

RESULTS AND DISCUSSION

Lifestyle diversification in Rhizobiales.The Rhizobiales taxa fall into four lifestyles, the nodule-associated, plant-associated, animal-associated, and free-living lifestyles. We inferred ancestral lifestyles of Rhizobiales using the maximum parsimony method implemented in Mesquite based upon the lifestyles of extant taxa (see Data Set S1 in the supplemental material for the complete list) and on their concatenated ribosomal protein phylogenies (see Materials and Methods). We did not use the maximum likelihood method in the main analysis because the short branches across the phylogeny can lead to overestimations of transition rates and thus to inaccurate ancestral lifestyle inferences (19), although we included a maximum likelihood analysis to assess the consistency of the results (see below). According to this procedure, the last common ancestor (LCA) of Rhizobiales was likely a free-living bacterium (Fig. 1), consistent with the fact that a vast majority of the branches of Rhizobiales that split early are represented by free-living members (Fig. 1). The nodule-associated lifestyle evolved multiple times, during which four major origins led to the formation of four well-known rhizobia genera, namely, Bradyrhizobium, Mesorhizobium, Rhizobium, and Sinorhizobium/Ensifer (nodes 2, 6, 5, and 4, respectively, in Fig. 1). Interestingly, none of the basal groups of these four rhizobia genera include nodule-associated bacteria: the basal members of Bradyrhizobium are free-living bacteria, those of Rhizobium are plant associates, and those of Sinorhizobium and Mesorhizobium take plant-associated or free-living lifestyles (Fig. 1). This suggests that the nodule-associated lifestyle in all of the four genera evolved relatively recently within each genus.

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Ancestral lifestyle reconstruction of the Rhizobiales. Ancestral lifestyles were inferred using the parsimony method in Mesquite. Branches in red, green, orange, and blue indicate nodule-associated, plant-associated, animal-associated, and free-living lifestyles, respectively. Numbered nodes represent the origins of the host-associated lifestyles (the nodule-associated, plant-associated, and animal-associated lifestyles). Zoomed-in lineages represent species that split early within the Bradyrhizobium, Mesorhizobium, Rhizobium, and Sinorhizobium. The three layers are each indicated with a distinct color surrounding the phylogeny to denote the presence/absence of key nodule-related genes (the nod, nif, and fix genes). The genes selected to represent each pathway/complex are shown in Text S1. Purple circles on the phylogeny represent nodes supported by IQ-Tree ultrafast bootstrap values of ≥95%. The outgroups are not shown.

DATA SET S1

(a) Lifestyle and genome information for all 1,264 Rhizobiales genomes downloaded from the NCBI RefSeq database (July 2017). Among them, a total of 655 genomes were analyzed in the study (column F), referred to as “the full set of genomes.” The genomes in this set were used in phylogenomic tree construction (Fig. 1), ancestral lifestyle inference (Fig. 1), lifestyle transition rate analysis (Fig. 2), and identification of lifestyle-associated genes. A set of 176 representative genomes of their OTUs (column G) were used in molecular dating (Fig. 3) and gene gain/loss analysis (Fig. 4). The reference genomes are those commonly used in academics and industry (column H) and were preferentially selected as the OTU representatives (see Materials and Methods). The source(s) of the isolation site of each genome, if available, is derived from the NCBI BioSample database, other public databases, or the literature (references are shown at the end). (b) The full list of lifestyle-correlated gene families that belong to the functional categories/pathways/complexes described in the main text. (c to f) Gene families showing significant correlation with the corresponding lifestyle with the full set of 655 genomes using BayesTraits. The details of the methods used for identification of the lifestyle-correlated gene families are provided in Materials and Methods. Only those with a q value (adjusted P value, Bonferroni correction) below 0.05 and rate differences above zero are displayed. (c) Nodule-associated (NA) lifestyle. (d) Animal-associated (AA) lifestyle. (e) Plant-associated (PA) lifestyle. (f) Free-living (FL) lifestyle. Download Data Set S1, XLSX file, 0.7 MB.
Copyright © 2020 Wang et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TEXT S1

(a) Supplemental discussion. (b) Selection of genes included in Fig. 1 and 4. (c) Calibration information. Download Text S1, DOCX file, 0.1 MB.
Copyright © 2020 Wang et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Animal associates showed two major independent origins leading to Bartonella (node 7 in Fig. 1) and Brucella (node 8 in Fig. 1). Since some of their close relatives, Ochrobactrum, for example, are increasingly recognized as opportunistic pathogens, there is a chance that their common ancestor had already adapted to an animal-associated lifestyle. We also identified two major origins of plant-associated members: one within the Methylobacterium (node 1 in Fig. 1) and the other corresponding to the LCA of the Rhizobiaceae (node 3 in Fig. 1). Species from the Rhizobiaceae further diversified into nodulating members of Rhizobium and Sinorhizobium, while others, in particular, those from Agrobacterium, plausibly retained the ancestral plant-associated lifestyle (and evolved plant pathogens later [17]) (Fig. 1). It is possible that a few host-associated lineages have a deeper origin, including Rhizobiaceae, where nodulating members exhibited mosaic distributions in phylogeny (Fig. 1), implying a more complex evolutionary history. We therefore updated the analysis by inclusion of the recently available genomes of a large number (>300) of Rhizobiaceae isolates from the roots of several nonlegume plants (8). Our updated phylogenomic tree showed that the phylogenetic positions of the major group of the nodulating Rhizobium (see Fig. S1a in the supplemental material) and of the nodulating Sinorhizobium (Fig. S1b) remained congruent with those shown in Fig. 1, suggesting the robustness of ancestral lifestyle reconstruction in these lineages.

FIG S1

The phylogeny of Rhizobium (a) and Sinorhizobium (b), including isolates sequenced from a recent study (8). The phylogenetic tree was constructed based on the concatenated alignment of ribosomal proteins. Purple circles on the phylogeny represent the nodes that are supported by IQ-Tree ultrafast bootstrap values of at least 95%. The layer adjacent to the taxon name denotes different lifestyles, with red, green, orange, and blue representing nodule-associated, plant-associated, animal-associated, and free-living lifestyle, respectively (as shown in Fig. 1). The black circles in the next layer denote those sequenced in a more recent study (8). The red star indicates the origin of nodule association inferred as described for Fig. 1. Download FIG S1, PDF file, 0.1 MB.
Copyright © 2020 Wang et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

To validate the origins of the host-associated lifestyles with the parsimony-based approach described above, we applied the maximum likelihood method implemented in BayesTraits (see Materials and Methods). Some of the host associations were estimated to have a deeper origin (nodes 1, 3, and 4 in Fig. 1), but the general pattern did not change (Fig. S2a). We further assessed the impacts of taxon sampling (Fig. S2b) and the set of genes used to build phylogeny (Fig. S2c) on lifestyle reconstruction. In both cases, the reconstructed lifestyles are in good agreement with the data in Fig. 1.

FIG S2

Ancestral lifestyle reconstruction using alternative strategies. (a) Lifestyles predicted using the maximum likelihood method with BayesTraits. (b) Lifestyles predicted using the maximum parsimony method with a broader taxon sampling, i.e., without removing genomes with completeness below 99%. (c) Lifestyles predicted using the maximum parsimony method based on the phylogeny constructed with 101 single-copy genes present in ≥95% of analyzed genomes (see Materials and Methods). Download FIG S2, PDF file, 0.4 MB.
Copyright © 2020 Wang et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

To further quantify the tendency of lifestyle transitions, we calculated transition rates between lifestyles using BayesTraits (see Materials and Methods). The transition rate from nodule association to nonnodule association was about eight times the rate calculated for the opposite direction (27.62 versus 3.30; log Bayes factor [logBF] = 72.04) (Fig. 2b), indicating that the tendency for rhizobia to lose their nodulating ability is much stronger than the gain of nodulating ability for nonrhizobia. The rate of transition to nodule associations was the highest for the plant-associated lifestyle (Fig. 2a) and was significantly higher than that determined for the animal-associated lifestyle (logBF = 3.28) but not that determined for the free-living lifestyle (logBF = 0.90). The transition from the plant-associated lifestyle to the free-living lifestyle is more likely to occur than the reverse transition, as reflected by the significantly higher transition rate (logBF = 3.48) (Fig. 2a). The rate of transition from an animal-associated lifestyle to any of the other three lifestyles was significantly lower than the rates of transition in the opposite direction and was not significantly different from zero (Fig. 2a). The pattern held when we combined all of the non-animal-associated lifestyles (logBF < 2; Fig. 2b). Nodule- and animal-associated bacteria exhibited markedly distinct patterns in lifestyle transition. Animal pathogens are often subjected to strong bottlenecks, which may lead to genetic drift and to massive losses of genes that cannot confer large advantages (2), and as a result, this could make it difficult for them to evolve into other lifestyles or even to reach an evolutionary end. On the other side, rhizobia need to dwell various habitats, including nodules, bulk soils, and rhizospheres; the free-living stage in soils and rhizosphere may provide rhizobia large population sizes, acting against the tendency of the decrease of selection effectiveness observed in animal pathogens (2). Rhizobia are also equipped with genes allowing them to thrive in all these habitats. Thus, they might easily shift back to nonrhizobia by loss of symbiosis genes, in particular, when their hosts grow in nitrogen-rich soils, as the benefit of carrying these symbiosis genes becomes vanishingly low (20).

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Rates of transitions between lifestyles. The width of each arrow is proportional to the transition rate (log transformed). Gray dashed lines denote transition rates that are not significantly different from zero (logBF < 2). The posterior distributions of transition rates derived by the reversible-jump MCMC method are shown as histograms. The means of transition rates and Z-scores are displayed adjacent to the histogram. Expressed as percentages, Z-values represent the proportions of sampled runs in the Markov chain where the transition rate is assigned to be zero. Z-values allow examination of the contingency of evolutionary transitions, with higher values signifying lifestyle transitions that are less likely to occur in evolution. Nodule-associated, plant-associated, animal-associated, and free-living lifestyles are abbreviated as NA, PA, AA, and FL, respectively. The color that represents each lifestyle is the same as that defined for Fig. 1. (a) Pairwise comparisons of transition rates between lifestyles. (b) Comparisons of transition rates between nodule-associated lifestyles (NA) and non-nodule-associated lifestyles (AA, PA, and FL) and between animal-associated lifestyles (AA) and non-animal-associated lifestyles (NA, PA, and FL).

Coevolutionary history of the associations between Rhizobiales species and their hosts.A large proportion (76%) of Rhizobiales members sampled in the present study were isolated from a host-associated environment, making it possible to explore the coevolution between Rhizobiales and their hosts. Our strategy started by estimating the time of origin of each lifestyle by the use of careful molecular clock analyses and comparing each estimated time of origin with that of their hosts which were recorded in fossils. For computational efficiency, we selected 176 representative genomes based on the operational taxonomic unit (OTU) at the 16S rRNA gene sequence identity level of 98.7% (21) (see Materials and Methods). Molecular clock analyses showed that the LCA of the Rhizobiales occurred 1,569 million years ago (Mya) (95% highest posterior density [HPD] interval, 1,667 to 1,447 Mya), which greatly predated the origin of their hosts (Fig. 3) (22–24; also see below). This lends strong support for the idea of a free-living LCA of Rhizobiales illustrated by the ancestral lifestyle reconstruction (Fig. 1; see also Text S1 in the supplemental material).

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Time tree of the Rhizobiales with selected taxa from Proteobacteria as the outgroup. Divergence time was estimated using MCMCTree on the species phylogeny shown in Fig. 1 (performed with only the 176 representative species). Nodes marked with a gray circle represent the calibration points in Proteobacteria and cyanobacteria (see Text S1 for details). Branch colors are based on the ancestral lifestyle reconstruction shown in Fig. 1. Vertical bars in green, orange, and red indicate the origin times of the primary hosts of extant Rhizobiales, namely, land plants, mammals, and legumes, respectively. The panel in the lower left corner of the figure shows a zoomed-in view of the evolutionary timeline of Bartonella, where the hosts are also shown next to the tips of the phylogeny. Node bars denote the 95% HPD interval of posterior dates (for key nodes related to the origins of host association only).

With an origin at 116 Mya (95% HPD interval, 146 to 88 Mya), nodulating Bradyrhizobium showed the earliest origin among the aforementioned four major rhizobia lineages (Fig. 3), very close to the origin time of another nodulating clade, the Azorhizobium clade, at 134 Mya (95% HPD interval, 180 to 91 Mya). Hence, the first alphaproteobacterial rhizobium was likely associated with Bradyrhizobium or Azorhizobium. The origins of nodulating lineages of Mesorhizobium, Sinorhizobium, and Rhizobium generally postdated those of Bradyrhizobium and Azorhizobium (Fig. 1), implying that the nodulating ability of the former lineages was acquired from the latter lineages or their relatives. Notably, the divergence time of the Rhizobiales lineages should be understood as a span of the posterior age estimate (indicated by the 95% HPD interval), rather than as a time point. Recent molecular dating studies suggest that legumes (the host of alphaproteobacterial rhizobia) originated 110 to 65 Mya (25) (see also references 16, 26, and 27 and Text S1). Hence, there were considerable overlaps in the origin times of legumes and alphaproteobacterial rhizobia (e.g., Azorhizobium, Bradyrhizobium, and Mesorhizobium). The origin of leguminous nodules, assuming its first appearance in the LCA of legumes (28–30), was thus roughly contemporaneous with that of nodulating alphaproteobacterial rhizobia (Fig. 3). Alternatively, certain rhizobial lineages (Azorhizobium or Bradyrhizobium) might have originated a bit earlier than legumes and evolved the nodulating ability after the emergence of legumes in independent lineages. Additionally, the uncertainties might result from methodological limits of molecular dating in organisms with ancient origins (31), such as the Rhizobiales. Nonetheless, our result remained compatible with the general hypothesis of the coevolution of alphaproteobacterial rhizobia and their hosts (Text S1).

The plant-associated members in Rhizobiaceae and Methylobacterium originated 280 Mya (95% HPD interval, 314 to 242 Mya) and 207 Mya (95% HPD interval, 241 to 171 Mya), respectively (Fig. 3), >100 million years (My) after the terrestrialization of plants (24). In support of this idea, these lineages are found in the microbial communities at the phyllosphere or endosphere of various land plants (8, 32). Note that, due to the lack of species that split early, the LCA shared by Rhizobiaceae and its sister group Martelella was ambiguous in terms of its lifestyle, and there is a chance that it had already adapted to a plant-associated lifestyle (Fig. 3). As this LCA occurred at 417 Mya (95% HPD interval, 479 to 355 Mya) and coincided with the emergence of land plants at ∼470 Mya, when scarce nutrients and water severely limited the development of early land plants (3), the primitive interaction between these Rhizobiales lineages and the earliest land plants might have contributed to the successful terrestrialization of early plants.

Animal-associated Brucella and Bartonella displayed divergent patterns in terms of the origin time (Fig. 3). Brucella originated no earlier than 20 Mya, whereas Bartonella likely originated at 343 Mya (95% HPD interval, 393 to 293 Mya), much earlier than the emergence of mammals (22), suggestive of host shifts during Bartonella evolution. Recent efforts have revealed that basal lineages Bartonella apis and Bartonella tamiae are nutritional commensals inhabiting the guts of insects (33) whose origin (23) generally coincided with the emergence of Bartonella (Fig. 3). Further, unlike Brucella, where transmission is directly mediated by contact between mammals, Bartonella is transmitted across mammals via blood-sucking arthropod vectors. The evidence provided above suggests that ancestral lineages of Bartonella might have already lived closely with arthropods and that they became arthropod-transmitted mammalian pathogens ∼150 My later, coinciding with the emergence of mammals (22) (Fig. 3).

To accommodate potential biases in the time constraints used here, we performed additional analyses with various combinations of time constraints derived from fossil records and estimates from previous studies (Fig. S3; see Text S1 for details). In general, the patterns obtained from these new analyses remained unchanged (Fig. S3). For example, the 95% HPD interval of the estimates for the origins of both nodulating Bradyrhizobium and Azorhizobium generally overlapped that of nodulating plants across all combinations of calibrations (Fig. S3c), strengthening the idea that the first rhizobial lineages in Alphaproteobacteria lay in Bradyrhizobium or Azorhizobium. In addition, note the uncertainties inherent in time estimates, which can arise from the lack of fossils from close relatives, data partitioning, and lack of gene sets (Fig. S3). In general, removal of secondary calibrations within the Proteobacteria (sets 5 to 8), decreasing the number of partitions (set 11), or using the single-copy genes identified by OrthoFinder (set 12) led to more-ancient posterior times being estimated (Fig. S3). The largest posterior age value was observed in Set 11, where the sequence data were partitioned according to the scheme recommended by ModelFinder instead of being fully partitioned as performed for the other 11 calibration sets (Fig. S3c). This is consistent with previous findings revealing that estimated ages increase as more partitions come to be used (34). Different lineages also showed different extents of variation in the posterior age across calibration sets (Fig. S3c). Considering these uncertainties, one should be cautious in offering any conclusive arguments based on the time estimates. We think that the lack of fossils in most major bacterial lineages, the ambiguity in the data used to estimate the age of cyanobacteria fossils, and the large phylogenetic distance between cyanobacteria and other bacteria are the biggest challenges in bacterial divergence time estimation. Future studies may consider using the recently developed strategy based on horizontal gene transfer (HGT) (35, 36) to better resolve the evolutionary timeline for Rhizobiales and other bacteria.

FIG S3

Different calibrations and estimated divergence times for host associations. (a) Illustration of calibration points. (b) Time constraints used in the divergence time estimate. Set 1 was used to determine the data reported in the main text and Fig. 3. (c) HPD credible intervals (95%; vertical bar) of the time estimates of the nine origins of host-associated lifestyles for 12 different sets of time constraints (see Text S1 for details). The horizontal bar in each panel denotes the range of the times of origin of the hosts. The cartoons of legumes, mosses, mice, and bees represent legumes, land plants, mammals, and insects, respectively. Download FIG S3, PDF file, 0.3 MB.
Copyright © 2020 Wang et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Genome-wide identification of genes associated with nodulating members.Gene families conserved across diverse bacterial lineages of the same lifestyle show strong signals for convergent evolution and are therefore potentially important to bacterial adaptation in their common habitats (15, 16, 37). We integrated different protein family annotations, performed BayesTraits-based analysis to search for genes significantly associated with their lifestyles using the full set of the Rhizobiales genomes (see Materials and Methods), and elaborated on their putative roles in lifestyle adaptation of a few examples of particular interest, starting from those correlated with nodule-associated Rhizobiales (Data Set S1b; see also Data Set S1c for the full list).

Rhizobia-legume symbiosis is initiated by nodule formation and invasion (38, 39). As expected, genes participating in the biosynthesis of Nod factors (nod genes), which play key roles in inducing the host plant to form infection threads (38), were among the top-ranking significantly correlated genes (Fig. 1; see also Data Set S1). Genes encoding a type III secretion system (T3SS), T4SS, and T6SS, as well as bacterial effector proteins exported through them, were also detected. Distinct from other secretion systems which simply transport proteins/compounds out of cells, T3SS, T4SS, and T6SS can transport bacterial effectors to enable them to have direct communication with the eukaryotic cytosol (40). The roles of the T4SS and T6SS in symbiosis are less well studied than those of the T3SS. Our results suggest that those secretion systems, though not universally present at the strain level, may contribute to Rhizobia-legume symbiosis potentially by affecting host specificity and/or nodule growth (40). Genes involved in hopanoid lipid synthesis and modification were also identified, supporting the view of hopanoid as a molecule promoting bacterial mobility and attachment to root surfaces (41). Moreover, we detected several Rhizobia-associated genes that may help manipulate the level of plant hormones by rhizobia, such as the 1-aminocyclopropane-1-carboxylate deaminase (acds) which facilitates nodulation by degrading the precursor of the nodule inhibitor ethylene (42), and halimadienyl-diphosphate synthase which presumably participates in biosynthesis of gibberellin (42).

To supply legumes with N in nodular environments, rhizobia are equipped with the nitrogenase-encoding nif genes, which, as expected, were identified as genes enriched in nodulating bacteria (Fig. 1; Data Set S1b). A high rate of O2 respiration is needed to fulfill the energy requirement during N2 fixation (43). However, O2 can irreversibly inactivate nitrogenase. The microaerobic environment of nodules sets up a conflict of interest between N2 fixation and other biological processes in terms of O2 concentration (39). To cope with this issue, some rhizobia recycle the H2 released from N2 fixation via hydrogenase to reduce energy loss by regenerating chemical energy and removing H2 and O2, the reversible and irreversible inhibitors of nitrogen fixation, respectively (44). Accordingly, the hyd and hup genes encoding this H2 recycling system were significantly associated with nodule-associated taxa. In addition, many O2-high-affinity cytochromes and genes involved in their synthesis or regulation were significantly correlated. The best example is the Fix system (Fig. 1), including fixABCX participating in electron transfer to nitrogenase, fixJLKT responsible for nif regulation, and fixNOQP encoding high-affinity cytochrome oxidases, which collectively ensures high efficiency of aerobic biological processes under micro-oxic environment (39). Also significantly correlated were denitrifying enzymes, which are hypothesized to participate in detoxification of the signaling module nitric oxide in Bradyrhizobium and Sinorhizobium (45). In bacteroids, N2 fixation is driven by oxidation of host-supplied C4-dicarboxylic acids (43). Accordingly, the dct genes encoding the C4-dicarboxylate transporter were enriched in rhizobia. Note, however, that the low O2 level in bacteroids may inhibit the TCA cycle due to imbalanced redox state of NADH/NAD+, thereby disrupting N2 fixation (43). Rhizobia use the storage of polyhydroxybutyrate (PHB), which could account for half of the dry weight of rhizobial cells, as a way to stabilize cellular redox conditions, deposit redundant energy, and release the inhibition of the TCA cycle at low O2 conditions (43). Accordingly, several genes (e.g., phaZ, phbB, and phbC) responsible for the synthesis and degradation of PHB were found to be significantly associated.

Rhizobia need to not only persist in nodules but also thrive in soils where they compete with other bacteria for limited resources such as iron (Fe) and phosphorus (P). To overcome the low concentrations of soluble Fe available in soils, rhizobia use siderophores to import Fe by binding it in tight complexes (46). Indeed, several genes involved in siderophore biosynthesis and transport were significantly correlated with nodule-associated bacteria. Likewise, genes comprising the phn operon encoding a C-P lyase for utilization of phosphonate, which is commonly found in soils (47), were also among Rhizobia-associated genes. In addition, several genes involved in the biosynthesis of ectoine, a protectant against osmotic stress described in many soil bacteria (48), were significantly correlated with rhizobia, implying their role in adaptation to distinct osmotic pressures across various habitats.

Genes associated with animal-associated, plant-associated, and free-living lifestyles.Using the same procedure, we identified genes significantly correlated with each of the other three lifestyles. Fewer genes were significantly correlated with animal-associated than with nodulating bacteria (Data Set S1d), likely because of the limited taxonomic distribution of the animal-associated lifestyle within the Rhizobiales. Of most interest were the VirB genes encoding the type IV secretion system (Fig. 1), known as the essential virulence factors contributing to the success of infection by Bartonella and Brucella (49). Also significantly correlated with the animal-associated bacteria were the genes encoding adhesins, which are necessary for attaching Bartonella to host cells (49). Likewise, many genes that encode transporters of amino acids and metals (e.g., Fe3+ and Mg2+) showed significant correlation, consistent with their demonstrated essentiality in the survival of Bartonella and Brucella in mammalian bloodstream (50, 51) (Data Set S1b). Genes coding for multidrug efflux systems, which may help pathogens develop antimicrobial resistance during infection (52), were also enriched in animal-associated bacteria (Data Set S1b).

In addition, genes encoding the class Ib ribonucleotide reductase (RNR), a key enzyme in synthesis of the precursors of DNA in its replication (53), were among significantly associated genes in animal pathogens. Compared with other types of RNRs, class Ib RNR is special because it uses manganese instead of iron as the cofactor and is found only in bacteria and bacteriophages (53). The use of class Ib RNR in DNA replication thus likely helps Bartonella and Brucella survive in mammalian cells where iron is rare and facilitates escape from host defense via iron limitation, making the class Ib RNR a promising target for antibacterial and antiviral drug design, particularly given its absence in eukaryotes (thereby exerting less influence on host cells) (54).

A considerable number of genes significantly correlated with plant-associated strains encoded transporters covering diverse substrates, including sugars, amino acids, and oligopeptides (Data Set S1e). Many of them participate in the transport of rhamnose, a key component of lipopolysaccharide important to bacterial attachment to plants (55). Also included were those involved in cellulose synthesis and its activation, consistent with their roles in mediating plant-bacterium interaction (56). Free-living Rhizobiales species were derived from diverse isolation sites (Data Set S1). Hence, there were few “characteristic” genes shared by free-living Rhizobiales (Data Set S1f). Of most interest was the arc3 gene encoding a protein pumping arsenite from the cell, which showed the lowest P value among all genes correlated with this lifestyle, suggesting its role in adaptation to arsenite compounds, the most prevalent environmental toxic substances (57).

Genome expansion in the origin of nodulating lineages.Lifestyle-associated genes could be acquired before (i.e., preexisting traits), during, or after a lifestyle transition. These scenarios were differentiated using ancestral genome reconstruction of 176 representative Rhizobiales genomes with BadiRate (see Materials and Methods). Overall, this analysis predicted that Rhizobiales started with ∼2,100 genes at its LCA and gradually expanded to the current genome sizes (Fig. 4a). It also predicted repeated genome expansion based on the origin of the nodule-associated lifestyle and repeated genome reduction based on the origin of the animal-associated lifestyle (Fig. 4a), consistent with an early study that analyzed only nine genomes (12).

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Ancestral genome reconstruction and gene gains/losses during the evolution of Rhizobiales. (a) Genome size evolution in the Rhizobiales. The colors of the branches in the phylogeny represent the changes of gene numbers in the genomes. Ancestral genomes were inferred with BadiRate based on gene family clustering using OrthoFinder. (b and c) The phyletic pattern and inferred gain/loss of lifestyle-correlated genes in the evolution of nodule-associated lineages (b) and animal-associated lineages (c). Branches in red and orange indicate nodulating and animal-associated lineages, respectively. Selected pathways/complexes important to the host-associated lifestyle are shown next to the trees. Genes were selected based on the criteria detailed in Text S1. Solid and open circles denote the presence and absence of the corresponding genes, respectively. Presence/absence at ancestral nodes inferred by BadiRate is indicated adjacent to the node. Pathways or complexes with ≥50% or <50% of genes present are indicated by a solid or open square adjacent to the ancestral nodes, respectively, based on the reconstruction of ancestral genome content.

The nod and nif genes appeared to be the only ones that were repeatedly gained during the lifestyle transitions to the four major nodulating lineages, possibly via HGT from rhizobial lineages that had evolved early (14), consistent with their determining role in nodulation and nitrogen fixation by rhizobia. In contrast, genes constituting the Fix system, Dct transporter system, and phn operon and those involved in PHB metabolism were already present in their nonnodulating ancestors (Fig. 4b).

The evolutionary histories of most Rhizobia-associated genes were actually a mix of these scenarios where they preexisted in some rhizobia but were not recruited until lifestyle transition in others. Examples are T3SS, T4SS, T6SS, and denitrification genes (Fig. 4b). Such distinct patterns across lineages were closely related to the ancestral lifestyle from which the corresponding rhizobial lineage evolved, implying that these genes contributed to lifestyle adaptation in a lineage-specific manner (13). For instance, the nodulating Bradyrhizobium species likely arose from a free-living ancestor (Fig. 1). Accordingly, most genes involved in bacterium-host interactions, such as T3SS, T4SS, and T6SS genes, were predicted to be absent in the nonnodulating ancestor of Bradyrhizobium (Fig. 4b). T3SS was inferred to be acquired by the ancestor of nodulating Bradyrhizobium. This is interesting because photosynthetic Bradyrhizobium lineages that split early can use T3SS to nodulate legumes in the absence of nod (58). Considering its reported roles in symbiosis and conservation in Bradyrhizobium (40), our results imply that a gain of T3SS might have contributed to the ancestral transition of Bradyrhizobium from the free-living to the nodulating lifestyle (Fig. 4b). Distinct from Bradyrhizobium, nodulating Rhizobium and Sinorhizobium likely evolved from a plant-associated ancestor (8) (Fig. 1). Plant-associated bacteria produce various cell surface polysaccharides, among which EPS (exopolysaccharide) attracts the most attention because of its high abundance in the surrounding environment of bacteria (59). EPS has been shown to contribute to both symbiosis for Rhizobium and Sinorhizobium and plant pathogenesis for Agrobacterium (60). Our results suggest the presence of EPS I (succinoglycan) synthesis genes in the LCA of Rhizobiaceae and in the shared ancestors of Rhizobium, Sinorhizobium, and Agrobacterium (Fig. 4b). We hypothesize that the role of EPS in early Rhizobiaceae lineages may be simply the attachment to the plant surface and a protective barrier (as in the case of Agrobacterium) (60) and that it was coopted later in evolution by some rhizobia as a signaling molecule during invasion and infection thread formation. Note that, in addition to EPS I, rhizobia secrete other types of EPSs. However, given their low abundance and that their synthesis genes are less well characterized in rhizobia (59), those genes were not analyzed here. The case for Mesorhizobium was similar to those for Rhizobium and Sinorhizobium, where the genes participating in EPS synthesis and T6SS were inherited from their nonnodulating ancestor (Fig. 4b). This implies a plant-associated ancestor of Mesorhizobium, although the ancestral lifestyle of Mesorhizobium could not be fully resolved (Fig. 1). The conservation of T6SS in nodulating Mesorhizobium hints at its role in symbiosis, despite there being little knowledge about it (40). Further, we employed AnGST, a method fundamentally different from BadiRate (see Materials and Methods), to repeat the ancestral reconstruction analysis and obtained a similar pattern (Fig. S4).

FIG S4

The presence/absence of lifestyle-associated genes in extant ancestral nodes (solid/white circles) and key ancestral nodes (solid/white squares; inferred by AnGST) of major nodule-associated lineages (a) and animal-associated lineages (b). Download FIG S4, PDF file, 0.5 MB.
Copyright © 2020 Wang et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Despite significant genome expansion in the origins of nodulating lineages (Fig. 4a), genes encoding the gene transfer agent (GTA), a small tailed phage-like element that consists of ∼15 genes and mediates HGT of short genomic DNA segments (61), were possibly independently lost in all of the four major nodulating lineages (Fig. S5). A previous study showed that, while GTA was retained in Azorhizobium, its expression is suppressed when members of this rhizobial lineage form bacteroids (62). Also, despite the importance of HGT via plasmids or mobile islands to the Rhizobia-legume symbiosis (12), GTAs may not play a role in this process, because they preferentially mediate transfer of chromosomal DNA (63, 64) and can pack DNA segments of ∼4 kb (61), which are much smaller than those needed for rhizobial symbiosis. It is likely that GTAs became dispensable and thus were lost during the evolution toward the rhizobial lineages. This hypothesis can be tested by introducing the full set of GTA genes into rhizobia and assessing the effects on the growth of the bacteria.

FIG S5

The phyletic pattern of genes encoding the gene transfer agent (GTA) in the Rhizobiales. The first layer surrounding the phylogeny denotes the lifestyle, in which red, green, orange, and blue represent nodule-associated, plant-associated, animal-associated, and free-living lifestyles, respectively (as described for Fig. 1). The next two layers represent the presence/absence of genes encoding RcGTA (green), a widely distributed GTA first discovered in Rhodobacter capsulatus (90, 91), and BaGTA (red), which is specific to Bartonella (92), respectively. Homologs of RcGTA and BaGTA were identified by BLASTP using reference proteins as queries with the cutoff E value 1E−3 (90). Reference sequences of RcGTA and BaGTA were obtained from references 90 and 93, respectively. Download FIG S5, PDF file, 0.4 MB.
Copyright © 2020 Wang et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Exchange of genetic materials by HGT between bacteria inhabiting the same niches is very frequent (65). Given the coexistence of nonnodulating bacteria with rhizobia in both nodular and soil environments, it is likely that key genes involved in rhizobial symbiosis (e.g., the nod, nif, and fix genes) have frequently been transferred across bacteria over time (14). Paradoxically, until now, the hundreds of known rhizobial species were discovered in some 12 genera within only the Rhizobiales and Burkholderiales (14, 66). This suggests that HGT of these symbiosis genes is insufficient for the conversion from nonrhizobia to rhizobia. A recent study revealed that genome modifications of the recipient lineages are required for newly acquired symbiosis traits to function (67). In the present study, we found that preexisting traits and gene losses might have also contributed to the origin of rhizobia. Possibly, it is a combination of different mechanisms that make diverse lineages within the Rhizobiales independently evolve to be successful symbionts of legumes, which, however, needs to be further explored.

Genome reduction and gene gain/losses in animal pathogens.The origins of animal pathogens were characterized by significant genome reductions (Fig. 4a). In particular, many genes involved in amino acid metabolism were lost during the evolution of Bartonella, leading to the incompleteness of most amino acid biosynthesis pathways in this lineage (33). To compensate for these losses, Bartonella recruited several amino acid uptake genes during and after the lifestyle transition to make use of the available resources in the host (Fig. 4c; see also Fig. S4b). In contrast, most genes gained by Brucella are involved in the uptake of metals, including nickel and cobalt (Fig. 4c; see also Fig. S4b). These genes have been proven critical to the virulence of other Gram-negative bacteria (reviewed in reference 68) and might be similarly important to Brucella (18). Other genes, e.g., those encoding T4SS, class Ib RNR (ribonucleotide reductase), and multidrug exporters, were likely present at the LCA shared by Bartonella and Brucella (Fig. 4c; see also Fig. S4b). These genes might make related lineages “preadapted” to animal association and pathogenesis. Of note, some strains belonging to Ochrobactrum, representing the closest relatives of Brucella, have evolved as opportunistic pathogens. People with an indwelling medical device are most susceptible to the bloodstream infection by Ochrobactrum (69). This could result from its ability to resist antibiotics and adhere to synthetic materials (70), consistent with the idea that the LCA of Ochrobactrum and Brucella might have already harbored relevant genes (Fig. 4c; see also Fig. S4b). We speculate that these preexisting genes might have been involved in different but related pathways in the free-living ancestors and might evolved new functions in pathogenesis during evolutionary transition to the host-associated lifestyle. For example, the T4SS in Bartonella, which exclusively functions to modulate the bacterium-eukaryotic host cell interaction, likely evolved from an ancestral conjugation system involved in the formation of stable mating junctions (71). The findings described above also suggest that to limit the risk of emergence of pathogenic bacteria, it is important to watch for (i) those with genes that potentially facilitate adaptation or conversion to a pathogenic lifestyle in genomes and (ii) those that infect related hosts or that have already been adapted to the host (72). Such examples include the aforementioned Ochrobactrum. In addition, several members of the Afipia are associated with the blood of patients (Data Set S1), highlighting its potential as a reservoir for emerging diseases (73).

Caveats and concluding remarks.Like a recent study (37), we classified host-associated Rhizobiales based on their isolation sites. However, different lifestyles could represent dynamic alternatives available to the organism at different time points (7). In our data sets, three non-nodule-associated strains (Mesorhizobium sp. strain LCM 4577, Rhizobium phaseoli Ch24-10, and Rhizobium sp. strain NXC14) possess a complete set of nodABC genes, thus likely representing strains that are capable of nodulating legumes but that happened to be isolated elsewhere. Thirteen of the 244 nodule-associated strains do not carry any nodABC genes. While some of them belong to the photosynthetic Bradyrhizobium lineage which utilizes a Nod-independent strategy for nodulation (58), others may represent nonnodulating bacteria living in nodules. Moreover, although our computational identification of lifestyle-correlated genes provides useful insights into the development of microbe-based strategies for sustainable agriculture and for prevention of plant/animal diseases, it is based purely on bioinformatics analyses. The detailed functions of many of identified genes still await exploration. The overrepresentation of agriculturally and medically important strains in sequenced genomes (8), although accounted for by BayesTraits, should also be noted.

Given the prevalence of plant-associated members, a recent study hypothesized that the LCA of Rhizobiales was already capable of colonizing the root of plants (8). Here, we revisited this idea by showing that the highly diversified lifestyles in Rhizobiales likely originated from a free-living ancestor. We also estimated the emergence of Rhizobiales some 1,500 Mya, far predating the origins of their hosts, including land plants (∼470 Mya) (24), insects (∼480 Mya), (23) and mammals (∼180 Mya) (22). Such an ancient origin of Rhizobiales indicates that the free-living lifestyle was adopted by Rhizobiales over the first half of their evolutionary trajectory. Rhizobiales evolved to live in association with early terrestrial animals and plants that emerged during the Cambrian explosion in their middle age, and at a later stage, when nodulating plants emerged, they established the legume-rhizobium symbiosis which is among the best-known symbiosis relationships.

Our results also indicate that genomic plasticity is an important feature driving the evolution of diverse lifestyles in Rhizobiales. Both the recurrent gains of genes that determine lifestyle transition and the losses of genes that do not have big advantages underlie the evolutionary adaptations of Rhizobiales to different lifestyles. Some of these lifestyle-important genes, if acquired by ancestral lineages, might make their descendants preadapted to the development of certain lifestyles. In addition, HGTs of genes important for and common to Rhizobiales of different lifestyles (e.g., T4SS genes) across rhizobia, plant associates, and animal pathogens might facilitate lifestyle diversification across the entire order (6, 13).

MATERIALS AND METHODS

Collection of Rhizobiales genome sequences and lifestyles.We retrieved genome sequences and annotations of 1,264 Rhizobiales isolates (i.e., no metagenome-assembled genome was included) from NCBI RefSeq (last accessed July 2017), estimated the completeness of each genome using CheckM v1.0.7 with the default parameters (74), and collected their lifestyle information from the BioSample and BioProject database at NCBI and from literature (see Data Set S1 in the supplemental material for the complete information). Those with a completeness value lower than 99% or without a clearly documented lifestyle information were removed. We selected these stringent criteria because genomes with low completeness may confound gene gain/loss analysis. Genomes without an identified 16S rRNA gene were also excluded. The remaining 891 strains were classified into four lifestyles based on their isolation site (37). Strains isolated from nodules are referred to as nodule associated. Those isolated from other parts of plants, the rhizosphere included, were considered plant-associated strains (37). Those isolated from animals were defined as animal-associated strains. Bacteria that were not associated with any host were classified as free-living strains; among those strains, 31%, 43%, and 12% were isolated from bulk soil, freshwater, and the ocean, respectively.

Phylogenomic tree construction.Amino acid sequences of 53 conserved ribosomal protein genes, which are less affected by HGT and thereby generally reflect well the evolutionary relationships of the genomes carrying them (75), were retrieved based on the results of reversed PSI-BLAST (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi). Sequences were aligned using MAFFT v7.222 with default parameters (76). Alignments were trimmed using trimAL v1.4.1 (77) with the parameter “-st 0.001” and were concatenated for phylogenetic analysis. For genomes with identical ribosomal protein sequences, one genome was randomly chosen as the representative, which led to retrieval of a total of 655 genomes (Data Set S1) for subsequent analyses.

The phylogenomic tree of Rhizobiales was built using IQ-Tree v1.5.3 (78) with the parameter “-s alignment -spp partition -m MFP+MERGE -rcluster 10 -mset LG,WAG,JTT -mrate R -wbtl -bb 1000” and with six genomes from other Alphaproteobacteria lineages (Caulobacter crescentus CB15, Rhodobacter sphaeroides 2.4.1, Leisingera daeponensis DSM 23529, Sphingomonas melonis FR1, Porphyrobacter cryptus DSM 12079, and Rhodospirillum rubrum F11) used as the outgroup. This procedure automatically implements a data partition model for the concatenated alignment and determines the best-fit substitution model for each partition. A total of 1,000 ultrafast bootstrap approximations (79) were performed to get the support value for each node. Those with a support value of ≥95% were considered to be well-supported nodes (79). Phylogenies were visualized with iTOL v4 (https://itol.embl.de/).

Calculation of rates of transition between lifestyles and reconstruction of ancestral lifestyles.Rates of transition between lifestyles were calculated using the reversible-jump Markov chain Monte Carlo (MCMC) method implemented in BayesTraits v3.0.1 (19), which takes into account both the topology and the branch length of the species phylogeny and is widely used in comparative phylogenetic studies. The analysis was performed across the 1,000 ultrafast bootstrap trees generated by IQ-Tree and was run for 10,000,000 iterations (1,000 stones, each sampled for 10,000 iterations) after discarding the first 1,000,000 runs as representing the burn-in. All priors were set to an exponential with a mean of 10 (80). Convergence was checked using Tracer v1.6 (http://tree.bio.ed.ac.uk/software/tracer/). Alternative hypotheses were tested using log Bayes factor (logBF) calculated as 2 × (log marginal likelihood [model I] – log marginal likelihood [model II]) (19). Values of logBF between 2.0 and 5.0, between 5.0 and 10.0, and above 10.0 were considered positive evidence, strong evidence, and very strong evidence for support, respectively (19).

Ancestral lifestyles were inferred using the parsimony model in Mesquite v2.7.5 (https://www.mesquiteproject.org/), which finds the ancestral lifestyles that minimize the number of steps of lifestyle changes based on the species phylogeny and observed lifestyle distribution, as used in many studies (5, 15). We further validated major origins of host-associated lifestyles (i.e., the nodule-, plant-, and animal-associated lifestyles) with the maximum likelihood method using the MultiState module implemented in BayesTraits. To accommodate the fast transitions across short branches, which could lead to overestimations of transition rates and thus to inaccurate ancestral state reconstruction, we estimated the kappa parameter to be 0.076 in BayesTraits, which maximizes the likelihood of the inference (see Fig. S2a in the supplemental material). Ancestral lifestyles were also inferred based on the genome tree built with a complete taxon sampling (i.e., without removing any genomes with low completeness) (Fig. S2b) and on phylogenies constructed by the use of 101 single-copy genes identified by OrthoFinder v2 (81), allowing a gene family member to be absent in at most 5% genomes (Fig. S2c).

Estimating the evolutionary timeline of the Rhizobiales.For computational efficiency, we classified all Rhizobiales genomes into 171 OTUs at the 16S rRNA gene identity level of 98.7%. For host-associated lineages of Brucella, Rhizobium, and Sinorhizobium, we applied a cutoff of 16S rRNA identity at 99.5% (82) to increase the number of the OTUs of these lineages, as using a cutoff at 98.7% (21) can capture only two OTUs of each of those lineages and may lead to inaccurate time estimations. For each OTU, a reference strain’s genome was selected if available. Otherwise, one genome was randomly chosen (see Data Set S1 for the list). Further, to better constrain the origin times of nodulating lineages in Bradyrhizobium, we included Bradyrhizobium sp. strain URHD0069 and Bradyrhizobium sp. strain URHA0002, the two basal lineages of Bradyrhizobium adapting to a free-living lifestyle (Fig. 1). Selected genomes from Alpha-, Beta-, and Gammaproteobacteria were also included as additional outgroups (11).

The amino acid sequences of 25 universally conserved single-copy genes were chosen as characters for molecular dating, as used in a previous study (9). Their amino acid sequences were aligned using MAFFT v7.222 (76) and trimmed using TrimAl v1.4.1 (77). All alignments were manually checked to avoid poorly aligned sequences. The relaxed molecular clock analysis was performed using MCMCTree from PAML v4.9i (83), a widely used MCMC-based tool for molecular dating (31). Since the use of more partitions can improve the precision of divergence time estimates (34), each gene was allowed to have its own partition, resulting in a total of 25 partitions in MCMCTree analysis. We used the approximate likelihood calculation and independent rate model implemented in MCMCTree. The first 10,000 iterations of each MCMC chain were discarded as burn-in, and the chain was run for 50,000 iterations, sampling every two iterations. Convergence was checked by repeating the analyses described above (see “Data availability” below).

MCMCTree employs soft fossil constraints for each calibration point, allowing 2.5% of the posterior probability distribution to exceed the minimum or maximum ages specified by the user (83). For the majority of the calibration points, the priors of the minimum and maximum ages were based on fossil and geologic evidence from cyanobacteria (Fig. S3; see also Text S1 in the supplemental material), as such evidence is most abundant in cyanobacteria among bacteria (84). In addition, given the potential for inaccuracy due to the use of very distantly related lineages as calibration points (31), we included two calibration points within the Proteobacteria based on other literature (see Fig. S3a for the full list of calibration points and Text S1 for the justification for each calibration point). To check whether the results were biased by the selected time constraints, we performed the same analysis with other different combinations of time constraints based upon other studies (Fig. S3).

Functional annotation of gene families.Clustering of genes and protein domains basically followed a recent study (37). For each genome, we clustered proteins into families based upon the annotation resources of COG, TIGRFAM, and KEGG Orthology. E values derived from reverse PSI-BLAST that were equal to or lower than 1e−5 were considered to represent homologs. HMMER v3.1 (85) was used to perform domain annotation for each gene based on Pfam release 30.0 with the gathering threshold and an E value cutoff of 1e−5. Furthermore, to capture the remaining genes that might not have functional annotations in the databases mentioned above, OrthoFinder v2 (81) was used with a cutoff of 1e−10 to cluster protein sequences from all analyzed genomes into homology-based orthogroups, which were used to infer the ancestral genome content (Fig. 4a).

Identification of lifestyle-associated genes.Genes significantly associated with lifestyles were identified using BayesTraits v3.0.1 (19). BayesTraits takes a binary matrix composed of two traits (e.g., nodule-associated lifestyle versus non-nodule-associated lifestyle) and determines whether the presence/absence of a gene family and the change in the lifestyle occurred dependently. It allows the identification of genomic adaptations that are associated with parallel switches in lifestyles. We applied the likelihood ratio test (LRT) to test the evolutionary association between the gene family and a lifestyle. In brief, we calculated the likelihood ratio associated with each gene family as representing twice the difference between the log likelihoods of the two models (2ΔlnL, where ΔlnL = lnLindependent – lnLdependent) against a χ2 distribution with four degrees of freedom (15, 16, 86). To determine the direction of association (i.e., whether the lifestyle is correlated with the presence or absence of the gene), we calculated the transition rate difference as ΔQ = q21 + q31 + q34 + q24 − (q12 + q13 + q43 + q42) (86). As shown in Fig. S6, we would expect a positive correlation between the lifestyle and the presence of the gene (i.e., we would expect the shaded regions to dominate) for ΔQ values of >0 (86). We identified genes as positively associated with a lifestyle if their false-discovery-rate (FDR)-adjusted (Bonferroni method) P values derived from the LRT were lower than 0.05 and ΔQ values were above zero. It is worth noting that the same gene could be phylogenetically correlated with bacteria of two or more lifestyles if it meets the criteria described above.

FIG S6

Illustration of possible combinations of two binary characters (gene and lifestyle) and of the transition rates (qXX) between these combinations. The shaded combinations dominate when the gene is positively correlated with the lifestyle (in this example, nodule association) calculated by BayesTraits. For example, when the q24 value is greater than the q42 value, that means that in the presence of the gene, the conversion from a non-nodule-associated lifestyle to a nodule-associated lifestyle occurs more easily than the opposite process. When the q31value is greater than the q13 value, that means that in the absence of the gene, it is easier for nodule-associated bacteria to lose a nodule association than for non-nodule-associated bacteria to become nodule associated. The rate difference is calculated as follows: ΔQ = q21 + q31 + q34 + q24 − (q12 + q13 + q43 + q42). Significant positive correlation between the gene and the trait should meet the following two criteria: (i) FDR-corrected P values (Bonferroni correction) should be <0.05, and (ii) ΔQ values should be >0. Download FIG S6, PDF file, 0.01 MB.
Copyright © 2020 Wang et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Ancestral genome reconstruction and inference of gene gains and losses.Evolution of genome content via gene gains and losses was reconstructed using BadiRate v1.35 (87) with the parameter “–ep CSP –anc -rmodel BDI -bmodel FR.” By employing the Sankoff parsimony method, this parameter estimates lineage-specific rates of gains and losses for each gene family while considering both the topology and branch lengths of the phylogeny (87).

In addition, we applied AnGST (88) to examine whether the general pattern shown by the BadiRate analysis was robust. AnGST differs from BadiRate in that it infers gene gains and losses based on reconciliation of topological incongruences between gene trees and species trees. Additionally, AnGST accounts for uncertainties in gene trees by incorporating the first 100 bootstrap trees generated by the IQ-Tree ultrafast bootstrap approximation. It also allows users to specify the penalty scores of HGT, gene duplication, and gene loss, which were determined to be 3, 2, and 1, respectively, by minimizing genome size flux as suggested by a previous study (88). Gene families with fewer than four members were excluded from AnGST analysis, as IQ-Tree cannot build trees for them.

Data availability.The original sequences, phylogenetic trees, molecular dating analyses, and gene gain/loss results, together with the codes generating them (89), are available at https://figshare.com/articles/Evolutionary_timeline_and_genomic_plasticity_underlying_the_lifestyle_diversity_in_Rhizobiales/9849539.

FIG S7

Gene phylogenies of symbiosis genes built with concatenated alignment of nodABCIJ using IQ-Tree. Representative genomes from each Rhizobiales OTU used in molecular dating analysis together with rhizobia from Betaproteobacteria included in a recent study (14) and four Actinomycetes genomes carrying the nod genes (94) were used in gene tree construction. The trees are rooted based on midpoint rooting (where the root is put at the midpoint of the longest path between any two tips). The support values obtained from the ultrafast bootstrapping implemented in IQ-Tree are shown next to the corresponding nodes. Branches with support values of less than 80% are condensed, resulting in polytomy. Download FIG S7, PDF file, 0.05 MB.
Copyright © 2020 Wang et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S8

Venn graph illustrating the overlap of identified genes significantly correlated with nodule associations between the representative genome and the full set of genomes. The following two methods were applied to adjust the P values to correct the inflation of type I errors in multiple comparisons for the representative genome set: the Bonferroni method (95), which is very conservative, and the Benjamini-Hochberg (BH) method (96), which is more relaxed. The results are compared with those identified using the full genome set (Bonferroni correction). Those with a rate difference above zero and an adjusted P value lower than 0.05 are displayed (see also Materials and Methods). The correlated genes for the full and representative genomes are available in Data Set S1c and at https://figshare.com/articles/Evolutionary_timeline_and_genomic_plasticity_underlying_the_lifestyle_diversity_in_Rhizobiales/9849539, respectively. Download FIG S8, PDF file, 0.02 MB.
Copyright © 2020 Wang et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

ACKNOWLEDGMENTS

We are particularly grateful to Zhichao Zhou for his valuable contributions in data preparation and collection in the initiation of the project. We thank Ying Chang for discussion on land plant origin, Hao Ding for comments on GTA, Jianhua Huang for insights into plant-bacterium interactions, Asaf Levy for the identification of lifestyle-associated genes, and Anna Skorupska for information regarding bacterial EPS. We appreciate the help of Stilianos Louca, Julie Marin, Patricia Sánchez-Baracaldo, Ying Sun, and Hao Zhang for molecular dating analysis. We are grateful to Davids Emms, Lawrence David, and Mario Dos Reis Barros for their guidance in the use of Orthofinder, AnGST, and MCMCTree, respectively. In addition, we thank Xiaoyuan Feng, Tianhua Liao, Jie Liu, Minglei Ren, Jinjin Tao, and Haowei Zhang for their helpful discussion.

This work was supported by the Shenzhen Science and Technology Committee (JCYJ20180508161811899), the Hong Kong Research Grants Council Area of Excellence Scheme (AoE/M-403/16), the impact postdoctoral fellowship scheme (IPDFS) of CUHK, and the Direct Grant of CUHK (4053257 and 3132809).

We declare that we have no competing interests.

FOOTNOTES

    • Received May 16, 2020.
    • Accepted June 19, 2020.
  • Copyright © 2020 Wang et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Toft C,
    2. Andersson S
    . 2010. Evolutionary microbial genomics: insights into bacterial host adaptation. Nat Rev Genet 11:465–475. doi:10.1038/nrg2798.
    OpenUrlCrossRefPubMedWeb of Science
  2. 2.↵
    1. Moran NA,
    2. Wernegreen JJ
    . 2000. Lifestyle evolution in symbiotic bacteria: insights from genomics. Trends Ecol Evol 15:321–328. doi:10.1016/s0169-5347(00)01902-9.
    OpenUrlCrossRefPubMedWeb of Science
  3. 3.↵
    1. Martin FM,
    2. Uroz S,
    3. Barker DG
    . 2017. Ancestral alliances: plant mutualistic symbioses with fungi and bacteria. Science 356:eaad450. doi:10.1126/science.aad4501.
    OpenUrlCrossRef
  4. 4.↵
    1. Beinart RA
    . 2019. The significance of microbial symbionts in ecosystem processes. mSystems 4:e00127-19. doi:10.1128/mSystems.00127-19.
    OpenUrlAbstract/FREE Full Text
  5. 5.↵
    1. Sachs JL,
    2. Skophammer RG,
    3. Regus JU
    . 2011. Evolutionary transitions in bacterial symbiosis. Proc Natl Acad Sci U S A 108:10800–10807. doi:10.1073/pnas.1100304108.
    OpenUrlAbstract/FREE Full Text
  6. 6.↵
    1. Carvalho FM,
    2. Souza RC,
    3. Barcellos FG,
    4. Hungria M,
    5. Vasconcelos A
    . 2010. Genomic and evolutionary comparisons of diazotrophic and pathogenic bacteria of the order Rhizobiales. BMC Microbiol 10:37. doi:10.1186/1471-2180-10-37.
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. Denison RF,
    2. Kiers ET
    . 2004. Lifestyle alternatives for rhizobia: mutualism, parasitism, and forgoing symbiosis. FEMS Microbiol Lett 237:187–193. doi:10.1016/j.femsle.2004.07.013.
    OpenUrlCrossRefPubMedWeb of Science
  8. 8.↵
    1. Garrido-Oter R,
    2. Nakano RT,
    3. Dombrowski N,
    4. Ma KW,
    5. McHardy AC,
    6. Schulze-Lefert P
    , AgBiome Team. 2018. Modular traits of the Rhizobiales root microbiota and their evolutionary relationship with symbiotic rhizobia. Cell Host Microbe 24:155–167.e5. doi:10.1016/j.chom.2018.06.006.
    OpenUrlCrossRef
  9. 9.↵
    1. Battistuzzi FU,
    2. Hedges SB
    . 2009. A major clade of prokaryotes with ancient adaptations to life on land. Mol Biol Evol 26:335–343. doi:10.1093/molbev/msn247.
    OpenUrlCrossRefPubMedWeb of Science
  10. 10.↵
    1. Sipos G,
    2. Prasanna AN,
    3. Walter MC,
    4. O'Connor E,
    5. Bálint B,
    6. Krizsán K,
    7. Kiss B,
    8. Hess J,
    9. Varga T,
    10. Slot J,
    11. Riley R,
    12. Bóka B,
    13. Rigling D,
    14. Barry K,
    15. Lee J,
    16. Mihaltcheva S,
    17. LaButti K,
    18. Lipzen A,
    19. Waldron R,
    20. Moloney NM,
    21. Sperisen C,
    22. Kredics L,
    23. Vágvölgyi C,
    24. Patrignani A,
    25. Fitzpatrick D,
    26. Nagy I,
    27. Doyle S,
    28. Anderson JB,
    29. Grigoriev IV,
    30. Güldener U,
    31. Münsterkötter M,
    32. Nagy LG
    . 2017. Genome expansion and lineage-specific genetic innovations in the forest pathogenic fungi Armillaria. Nat Ecol Evol 1:1931–1941. doi:10.1038/s41559-017-0347-8.
    OpenUrlCrossRef
  11. 11.↵
    1. Luo H,
    2. Csűros M,
    3. Hughes AL,
    4. Moran MA
    . 2013. Evolution of divergent life history strategies in marine alphaproteobacteria. mBio 4:e00373-13. doi:10.1128/mBio.00373-13.
    OpenUrlAbstract/FREE Full Text
  12. 12.↵
    1. Boussau B,
    2. Karlberg EO,
    3. Frank AC,
    4. Legault BA,
    5. Andersson S
    . 2004. Computational inference of scenarios for alpha-proteobacterial genome evolution. Proc Natl Acad Sci U S A 101:9722–9727. doi:10.1073/pnas.0400975101.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    1. Tian CF,
    2. Zhou YJ,
    3. Zhang YM,
    4. Li QQ,
    5. Zhang YZ,
    6. Li DF,
    7. Wang S,
    8. Wang J,
    9. Gilbert LB,
    10. Li YR,
    11. Chen WX
    . 2012. Comparative genomics of rhizobia nodulating soybean suggests extensive recruitment of lineage-specific genes in adaptations. Proc Natl Acad Sci U S A 109:8629–8634. doi:10.1073/pnas.1120436109.
    OpenUrlAbstract/FREE Full Text
  14. 14.↵
    1. Remigi P,
    2. Zhu J,
    3. Young JPW,
    4. Masson-Boivin C
    . 2016. Symbiosis within symbiosis: evolving nitrogen-fixing legume symbionts. Trends Microbiol 24:63–75. doi:10.1016/j.tim.2015.10.007.
    OpenUrlCrossRef
  15. 15.↵
    1. Simon M,
    2. Scheuner C,
    3. Meier-Kolthoff JP,
    4. Brinkhoff T,
    5. Wagner-Döbler I,
    6. Ulbrich M,
    7. Klenk HP,
    8. Schomburg D,
    9. Petersen J,
    10. Göker M
    . 2017. Phylogenomics of Rhodobacteraceae reveals evolutionary adaptation to marine and non-marine habitats. ISME J 11:1483–1499. doi:10.1038/ismej.2016.198.
    OpenUrlCrossRef
  16. 16.↵
    1. Zhang H,
    2. Yoshizawa S,
    3. Sun Y,
    4. Huang Y,
    5. Chu X,
    6. González JM,
    7. Pinhassi J,
    8. Luo H
    . 2019. Repeated evolutionary transitions of flavobacteria from marine to non-marine habitats. Environ Microbiol 21:648–666. doi:10.1111/1462-2920.14509.
    OpenUrlCrossRef
  17. 17.↵
    1. Lassalle F,
    2. Planel R,
    3. Penel S,
    4. Chapulliot D,
    5. Barbe V,
    6. Dubost A,
    7. Calteau A,
    8. Vallenet D,
    9. Mornico D,
    10. Bigot T,
    11. Guéguen L,
    12. Vial L,
    13. Muller D,
    14. Daubin V,
    15. Nesme X
    . 2017. Ancestral genome estimation reveals the history of ecological diversification in Agrobacterium. Genome Biol Evol 9:3413–3431. doi:10.1093/gbe/evx255.
    OpenUrlCrossRef
  18. 18.↵
    1. Wattam AR,
    2. Foster JT,
    3. Mane SP,
    4. Beckstrom-Sternberg SM,
    5. Beckstrom-Sternberg JM,
    6. Dickerman AW,
    7. Keim P,
    8. Pearson T,
    9. Shukla M,
    10. Ward DV,
    11. Williams KP,
    12. Sobral BW,
    13. Tsolis RM,
    14. Whatmore AM,
    15. O'Callaghan D
    . 2014. Comparative phylogenomics and evolution of the brucellae reveal a path to virulence. J Bacteriol 196:920–930. doi:10.1128/JB.01091-13.
    OpenUrlAbstract/FREE Full Text
  19. 19.↵
    1. Meade A,
    2. Pagel M
    . 2016. BayesTraits V3 manual. http://www.evolution.rdg.ac.uk/BayesTraitsV3.0.1/Files/BayesTraitsV3.Manual.pdf.
  20. 20.↵
    1. Weese DJ,
    2. Heath KD,
    3. Dentinger BTM,
    4. Lau JA
    . 2015. Long-term nitrogen addition causes the evolution of less-cooperative mutualists. Evolution 69:631–642. doi:10.1111/evo.12594.
    OpenUrlCrossRefPubMed
  21. 21.↵
    1. Yarza P,
    2. Yilmaz P,
    3. Pruesse E,
    4. Glöckner FO,
    5. Ludwig W,
    6. Schleifer KH,
    7. Whitman WB,
    8. Euzéby J,
    9. Amann R,
    10. Rosselló-Móra R
    . 2014. Uniting the classification of cultured and uncultured bacteria and archaea using 16S rRNA gene sequences. Nat Rev Microbiol 12:635–645. doi:10.1038/nrmicro3330.
    OpenUrlCrossRefPubMed
  22. 22.↵
    1. dos Reis M,
    2. Inoue J,
    3. Hasegawa M,
    4. Asher RJ,
    5. Donoghue PCJ,
    6. Yang Z
    . 2012. Phylogenomic datasets provide both precision and accuracy in estimating the timescale of placental mammal phylogeny. Proc Biol Sci 279:3491–3500. doi:10.1098/rspb.2012.0683.
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. Misof B,
    2. Liu S,
    3. Meusemann K,
    4. Peters RS,
    5. Donath A,
    6. Mayer C,
    7. Frandsen PB,
    8. Ware J,
    9. Flouri T,
    10. Beutel RG,
    11. Niehuis O,
    12. Petersen M,
    13. Izquierdo-Carrasco F,
    14. Wappler T,
    15. Rust J,
    16. Aberer AJ,
    17. Aspöck U,
    18. Aspöck H,
    19. Bartel D,
    20. Blanke A,
    21. Berger S,
    22. Böhm A,
    23. Buckley TR,
    24. Calcott B,
    25. Chen J,
    26. Friedrich F,
    27. Fukui M,
    28. Fujita M,
    29. Greve C,
    30. Grobe P,
    31. Gu S,
    32. Huang Y,
    33. Jermiin LS,
    34. Kawahara AY,
    35. Krogmann L,
    36. Kubiak M,
    37. Lanfear R,
    38. Letsch H,
    39. Li Y,
    40. Li Z,
    41. Li J,
    42. Lu H,
    43. Machida R,
    44. Mashimo Y,
    45. Kapli P,
    46. McKenna DD,
    47. Meng G,
    48. Nakagaki Y,
    49. Navarrete-Heredia JL,
    50. Ott M, et al
    . 2014. Phylogenomics resolves the timing and pattern of insect evolution. Science 346:763–767. doi:10.1126/science.1257570.
    OpenUrlAbstract/FREE Full Text
  24. 24.↵
    1. Morris JL,
    2. Puttick MN,
    3. Clark JW,
    4. Edwards D,
    5. Kenrick P,
    6. Pressel S,
    7. Wellman CH,
    8. Yang Z,
    9. Schneider H,
    10. Donoghue P
    . 2018. The timescale of early land plant evolution. Proc Natl Acad Sci U S A 115:E2274–E2283. doi:10.1073/pnas.1719588115.
    OpenUrlAbstract/FREE Full Text
  25. 25.↵
    1. Barba-Montoya J,
    2. dos Reis M,
    3. Schneider H,
    4. Donoghue PCJ,
    5. Yang Z
    . 2018. Constraining uncertainty in the timescale of angiosperm evolution and the veracity of a cretaceous terrestrial revolution. New Phytol 218:819–834. doi:10.1111/nph.15011.
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. Sauquet H,
    2. Von Balthazar M,
    3. Magallón S,
    4. Doyle JA,
    5. Endress PK,
    6. Bailes EJ,
    7. Barroso De Morais E,
    8. Bull-Hereñu K,
    9. Carrive L,
    10. Chartier M,
    11. Chomicki G,
    12. Coiro M,
    13. Cornette R,
    14. El Ottra JHL,
    15. Epicoco C,
    16. Foster CSP,
    17. Jabbour F,
    18. Haevermans A,
    19. Haevermans T,
    20. Hernández R,
    21. Little SA,
    22. Löfstrand S,
    23. Luna JA,
    24. Massoni J,
    25. Nadot S,
    26. Pamperl S,
    27. Prieu C,
    28. Reyes E,
    29. Dos Santos P,
    30. Schoonderwoerd KM,
    31. Sontag S,
    32. Soulebeau A,
    33. Staedler Y,
    34. Tschan GF,
    35. Wing-Sze Leung A,
    36. Schönenberger J
    . 2017. The ancestral flower of angiosperms and its early diversification. Nat Commun 8:16047. doi:10.1038/ncomms16047.
    OpenUrlCrossRef
  27. 27.↵
    1. Magallón S,
    2. Gómez-Acevedo S,
    3. Sánchez-Reyes LL,
    4. Hernández-Hernández T
    . 2015. A metacalibrated time-tree documents the early rise of flowering plant phylogenetic diversity. New Phytol 207:437–453. doi:10.1111/nph.13264.
    OpenUrlCrossRefPubMed
  28. 28.↵
    1. van Velzen R,
    2. Doyle JJ,
    3. Geurts R
    . 2019. A resurrected scenario: single gain and massive loss of nitrogen-fixing nodulation. Trends Plant Sci 24:49–57. doi:10.1016/j.tplants.2018.10.005.
    OpenUrlCrossRef
  29. 29.↵
    1. Doyle JJ
    . 2011. Phylogenetic perspectives on the origins of nodulation. Mol Plant-Microbe Interact 24:1289–1295. doi:10.1094/MPMI-05-11-0114.
    OpenUrlCrossRefPubMedWeb of Science
  30. 30.↵
    1. Griesmann M,
    2. Chang Y,
    3. Liu X,
    4. Song Y,
    5. Haberer G,
    6. Crook MB,
    7. Billault-Penneteau B,
    8. Lauressergues D,
    9. Keller J,
    10. Imanishi L,
    11. Roswanjaya YP,
    12. Kohlen W,
    13. Pujic P,
    14. Battenberg K,
    15. Alloisio N,
    16. Liang Y,
    17. Hilhorst H,
    18. Salgado MG,
    19. Hocher V,
    20. Gherbi H,
    21. Svistoonoff S,
    22. Doyle JJ,
    23. He S,
    24. Xu Y,
    25. Xu S,
    26. Qu J,
    27. Gao Q,
    28. Fang X,
    29. Fu Y,
    30. Normand P,
    31. Berry AM,
    32. Wall LG,
    33. Ané JM,
    34. Pawlowski K,
    35. Xu X,
    36. Yang H,
    37. Spannagl M,
    38. Mayer KFX,
    39. Wong GKS,
    40. Parniske M,
    41. Delaux PM,
    42. Cheng S
    . 2018. Phylogenomics reveals multiple losses of nitrogen-fixing root nodule symbiosis. Science 361:eaat1743. doi:10.1126/science.aat1743.
    OpenUrlAbstract/FREE Full Text
  31. 31.↵
    1. Dos Reis M,
    2. Donoghue PCJ,
    3. Yang Z
    . 2016. Bayesian molecular clock dating of species divergences in the genomics era. Nat Rev Genet 17:71–80. doi:10.1038/nrg.2015.8.
    OpenUrlCrossRefPubMed
  32. 32.↵
    1. Kutschera U
    . 2007. Plant-associated methylobacteria as co-evolved phytosymbionts: a hypothesis. Plant Signal Behav 2:74–78. doi:10.4161/psb.2.2.4073.
    OpenUrlCrossRefPubMed
  33. 33.↵
    1. Segers FH,
    2. Kešnerová L,
    3. Kosoy M,
    4. Engel P
    . 2017. Genomic changes associated with the evolutionary transition of an insect gut symbiont into a blood-borne pathogen. ISME J 11:1232–1244. doi:10.1038/ismej.2016.201.
    OpenUrlCrossRef
  34. 34.↵
    1. Rannala B,
    2. Yang Z
    . 2007. Inferring speciation tunes under an episodic molecular clock. Syst Biol 56:453–466. doi:10.1080/10635150701420643.
    OpenUrlCrossRefPubMedWeb of Science
  35. 35.↵
    1. Davín AA,
    2. Tannier E,
    3. Williams TA,
    4. Boussau B,
    5. Daubin V,
    6. Szöllősi GJ
    . 2018. Gene transfers can date the tree of life. Nat Ecol Evol 2:904–909. doi:10.1038/s41559-018-0525-3.
    OpenUrlCrossRef
  36. 36.↵
    1. Wolfe JM,
    2. Fournier GP
    . 2018. Horizontal gene transfer constrains the timing of methanogen evolution. Nat Ecol Evol 2:897–903. doi:10.1038/s41559-018-0513-7.
    OpenUrlCrossRef
  37. 37.↵
    1. Levy A,
    2. Salas Gonzalez I,
    3. Mittelviefhaus M,
    4. Clingenpeel S,
    5. Herrera Paredes S,
    6. Miao J,
    7. Wang K,
    8. Devescovi G,
    9. Stillman K,
    10. Monteiro F,
    11. Rangel Alvarez B,
    12. Lundberg DS,
    13. Lu TY,
    14. Lebeis S,
    15. Jin Z,
    16. McDonald M,
    17. Klein AP,
    18. Feltcher ME,
    19. Rio TG,
    20. Grant SR,
    21. Doty SL,
    22. Ley RE,
    23. Zhao B,
    24. Venturi V,
    25. Pelletier DA,
    26. Vorholt JA,
    27. Tringe SG,
    28. Woyke T,
    29. Dangl JL
    . 2017. Genomic features of bacterial adaptation to plants. Nat Genet 50:138–150. doi:10.1038/s41588-017-0012-9.
    OpenUrlCrossRef
  38. 38.↵
    1. Jones KM,
    2. Kobayashi H,
    3. Davies BW,
    4. Taga ME,
    5. Walker GC
    . 2007. How rhizobial symbionts invade plants: the Sinorhizobium - Medicago model. Nat Rev Microbiol 5:619–633. doi:10.1038/nrmicro1705.
    OpenUrlCrossRefPubMedWeb of Science
  39. 39.↵
    1. Prell J,
    2. Poole P
    . 2006. Metabolic changes of rhizobia in legume nodules. Trends Microbiol 14:161–168. doi:10.1016/j.tim.2006.02.005.
    OpenUrlCrossRefPubMedWeb of Science
  40. 40.↵
    1. Nelson MS,
    2. Sadowsky MJ
    . 2015. Secretion systems and signal exchange between nitrogen-fixing rhizobia and legumes. Front Plant Sci 6:491. doi:10.3389/fpls.2015.00491.
    OpenUrlCrossRefPubMed
  41. 41.↵
    1. Belin BJ,
    2. Busset N,
    3. Giraud E,
    4. Molinaro A,
    5. Silipo A,
    6. Newman DK
    . 2018. Hopanoid lipids: from membranes to plant-bacteria interactions. Nat Rev Microbiol 16:304–315. doi:10.1038/nrmicro.2017.173.
    OpenUrlCrossRef
  42. 42.↵
    1. Ferguson BJ,
    2. Mathesius U
    . 2014. Phytohormone regulation of legume-rhizobia interactions. J Chem Ecol 40:770–790. doi:10.1007/s10886-014-0472-7.
    OpenUrlCrossRefPubMed
  43. 43.↵
    1. Poole P,
    2. Allaway D
    . 2000. Carbon and nitrogen metabolism in rhizobium. Adv Microb Physiol 43:117–163. doi:10.1016/s0065-2911(00)43004-3.
    OpenUrlCrossRefPubMedWeb of Science
  44. 44.↵
    1. Evans HJ,
    2. Harker AR,
    3. Papen H,
    4. Russell SA,
    5. Hanus FJ,
    6. Zuber M
    . 1987. Physiology, biochemistry, and genetics of the uptake hydrogenase in Rhizobia. Annu Rev Microbiol 41:335–361. doi:10.1146/annurev.mi.41.100187.002003.
    OpenUrlCrossRefPubMedWeb of Science
  45. 45.↵
    1. Torres MJ,
    2. Rubia MI,
    3. Bedmar EJ,
    4. Delgado MJ
    . 2011. Denitrification in Sinorhizobium meliloti. Biochem Soc Trans 39:1886–1889. doi:10.1042/BST20110733.
    OpenUrlAbstract/FREE Full Text
  46. 46.↵
    1. Kraepiel AML,
    2. Bellenger JP,
    3. Wichard T,
    4. Morel F
    . 2009. Multiple roles of siderophores in free-living nitrogen-fixing bacteria. BioMetals 22:573–581. doi:10.1007/s10534-009-9222-7.
    OpenUrlCrossRefPubMed
  47. 47.↵
    1. Yu X,
    2. Doroghazi JR,
    3. Janga SC,
    4. Zhang JK,
    5. Circello B,
    6. Griffin BM,
    7. Labeda DP,
    8. Metcalf WW
    . 2013. Diversity and abundance of phosphonate biosynthetic genes in nature. Proc Natl Acad Sci U S A 110:20759–20764. doi:10.1073/pnas.1315107110.
    OpenUrlAbstract/FREE Full Text
  48. 48.↵
    1. Talibart R,
    2. Jebbar M,
    3. Gouesbet G,
    4. Himdi-Kabbab S,
    5. Wroblewski H,
    6. Blanco C,
    7. Bernard T
    . 1994. Osmoadaptation in rhizobia: ectoine-induced salt tolerance. J Bacteriol 176:5210–5217. doi:10.1128/jb.176.17.5210-5217.1994.
    OpenUrlAbstract/FREE Full Text
  49. 49.↵
    1. Ben-Tekaya H,
    2. Gorvel JP,
    3. Dehio C
    . 2013. Bartonella and Brucella-weapons and strategies for stealth attack. Cold Spring Harb Perspect Med 3:a010231. doi:10.1101/cshperspect.a010231.
    OpenUrlAbstract/FREE Full Text
  50. 50.↵
    1. Mavris M,
    2. Saenz H,
    3. Monteil M,
    4. Boulouis HJ,
    5. Dehio C,
    6. Vayssier-Taussat M
    . 2005. Characterization of genes involved in long-term bacteremia in mice by Bartonella birtlesii. Ann N Y Acad Sci 1063:312–314. doi:10.1196/annals.1355.050.
    OpenUrlCrossRefPubMed
  51. 51.↵
    1. Keriel A,
    2. Botella E,
    3. Estrach S,
    4. Bragagnolo G,
    5. Vergunst AC,
    6. Feral CC,
    7. O'Callaghan D
    . 2015. Brucella intracellular life relies on the transmembrane protein CD98 heavy chain. J Infect Dis 211:1769–1778. doi:10.1093/infdis/jiu673.
    OpenUrlCrossRefPubMed
  52. 52.↵
    1. Zwama M,
    2. Yamaguchi A,
    3. Nishino K
    . 2019. Phylogenetic and functional characterisation of the Haemophilus influenzae multidrug efflux pump AcrB. Commun Biol 2:340. doi:10.1038/s42003-019-0564-6.
    OpenUrlCrossRef
  53. 53.↵
    1. Martin JE,
    2. Imlay JA
    . 2011. The alternative aerobic ribonucleotide reductase of Escherichia coli, NrdEF, is a manganese-dependent enzyme that enables cell replication during periods of iron starvation. Mol Microbiol 80:319–334. doi:10.1111/j.1365-2958.2011.07593.x.
    OpenUrlCrossRefPubMed
  54. 54.↵
    1. Tholander F,
    2. Sjöberg BM
    . 2012. Discovery of antimicrobial ribonucleotide reductase inhibitors by screening in microwell format. Proc Natl Acad Sci U S A 109:9798–9803. doi:10.1073/pnas.1113051109.
    OpenUrlAbstract/FREE Full Text
  55. 55.↵
    1. Jofré E,
    2. Lagares A,
    3. Mori G
    . 2004. Disruption of dTDP-rhamnose biosynthesis modifies lipopolysaccharide core, exopolysaccharide production, and root colonization in Azospirillum brasilense. FEMS Microbiol Lett 231:267–275. doi:10.1016/S0378-1097(04)00003-5.
    OpenUrlCrossRefPubMedWeb of Science
  56. 56.↵
    1. Römling U,
    2. Galperin MY
    . 2015. Bacterial cellulose biosynthesis: diversity of operons, subunits, products, and functions. Trends Microbiol 23:545–557. doi:10.1016/j.tim.2015.05.005.
    OpenUrlCrossRefPubMed
  57. 57.↵
    1. Yang HC,
    2. Fu HL,
    3. Lin YF,
    4. Rosen BP
    . 2012. Pathways of arsenic uptake and efflux. Curr Top Membr 69:325–358. doi:10.1016/B978-0-12-394390-3.00012-4].
    OpenUrlCrossRefPubMed
  58. 58.↵
    1. Okazaki S,
    2. Tittabutr P,
    3. Teulet A,
    4. Thouin J,
    5. Fardoux J,
    6. Chaintreuil C,
    7. Gully D,
    8. Arrighi JF,
    9. Furuta N,
    10. Miwa H,
    11. Yasuda M,
    12. Nouwen N,
    13. Teaumroong N,
    14. Giraud E
    . 2016. Rhizobium-legume symbiosis in the absence of Nod factors: two possible scenarios with or without the T3SS. ISME J 10:64–74. doi:10.1038/ismej.2015.103.
    OpenUrlCrossRef
  59. 59.↵
    1. Skorupska A,
    2. Janczarek M,
    3. Marczak M,
    4. Mazur A,
    5. Król J
    . 2006. Rhizobial exopolysaccharides: genetic control and symbiotic functions. Microb Cell Fact 5:7. doi:10.1186/1475-2859-5-7.
    OpenUrlCrossRefPubMed
  60. 60.↵
    1. Leigh JA,
    2. Coplin DL
    . 1992. Exopolysaccharides in plant-bacterial interactions. Annu Rev Microbiol 46:307–346. doi:10.1146/annurev.mi.46.100192.001515.
    OpenUrlCrossRefPubMedWeb of Science
  61. 61.↵
    1. Lang AS,
    2. Zhaxybayeva O,
    3. Beatty JT
    . 2012. Gene transfer agents: phage-like elements of genetic exchange. Nat Rev Microbiol 10:472–482. doi:10.1038/nrmicro2802.
    OpenUrlCrossRefPubMed
  62. 62.↵
    1. Tsukada S,
    2. Aono T,
    3. Akiba N,
    4. Lee KB,
    5. Te Liu C,
    6. Toyazaki H,
    7. Oyaizu H
    . 2009. Comparative genome-wide transcriptional profiling of Azorhizobium caulinodans ORS571 grown under free-living and symbiotic conditions. Appl Environ Microbiol 75:5037–5046. doi:10.1128/AEM.00398-09.
    OpenUrlAbstract/FREE Full Text
  63. 63.↵
    1. Québatte M,
    2. Christen M,
    3. Harms A,
    4. Körner J,
    5. Christen B,
    6. Dehio C
    . 2017. Gene transfer agent promotes evolvability within the fittest subpopulation of a bacterial pathogen. Cell Syst 4:611–621.e6. doi:10.1016/j.cels.2017.05.011.
    OpenUrlCrossRef
  64. 64.↵
    1. Tomasch J,
    2. Wang H,
    3. Hall ATK,
    4. Patzelt D,
    5. Preusse M,
    6. Petersen J,
    7. Brinkmann H,
    8. Bunk B,
    9. Bhuju S,
    10. Jarek M,
    11. Geffers R,
    12. Lang AS,
    13. Wagner-Döbler I
    . 2018. Packaging of Dinoroseobacter shibae DNA into gene transfer agent particles is not random. Genome Biol Evol 10:356–369. doi:10.1093/gbe/evy005.
    OpenUrlCrossRef
  65. 65.↵
    1. Ochman H,
    2. Lawrence JG,
    3. Groisman EA
    . 2000. Lateral gene transfer and the nature of bacterial innovation. Nature 405:299–304. doi:10.1038/35012500.
    OpenUrlCrossRefPubMedWeb of Science
  66. 66.↵
    1. Masson-Boivin C,
    2. Giraud E,
    3. Perret X,
    4. Batut J
    . 2009. Establishing nitrogen-fixing symbiosis with legumes: how many rhizobium recipes? Trends Microbiol 17:458–466. doi:10.1016/j.tim.2009.07.004.
    OpenUrlCrossRefPubMedWeb of Science
  67. 67.↵
    1. Clerissi C,
    2. Touchon M,
    3. Capela D,
    4. Tang M,
    5. Cruveiller S,
    6. Genthon C,
    7. Lopez-Roques C,
    8. Parker MA,
    9. Moulin L,
    10. Masson-Boivin C,
    11. Rocha EPC
    . 2018. Parallels between experimental and natural evolution of legume symbionts. Nat Commun 9:2264. doi:10.1038/s41467-018-04778-5.
    OpenUrlCrossRef
  68. 68.↵
    1. Palmer LD,
    2. Skaar EP
    . 2016. Transition metals and virulence in bacteria. Annu Rev Genet 50:67–91. doi:10.1146/annurev-genet-120215-035146.
    OpenUrlCrossRefPubMed
  69. 69.↵
    1. Mudshingkar S,
    2. Choure A,
    3. Palewar M,
    4. Dohe V,
    5. Kagal A
    . 2013. Ochrobactrum anthropi: an unusual pathogen: are we missing them? Indian J Med Microbiol 31:306–308. doi:10.4103/0255-0857.115664.
    OpenUrlCrossRefPubMed
  70. 70.↵
    1. Alnor D,
    2. Frimodt-Møller N,
    3. Espersen F,
    4. Frederiksen W
    . 1994. Infections with the unusual human pathogens Agrobacterium species and Ochrobactrum anthropi. Clin Infect Dis 18:914–920. doi:10.1093/clinids/18.6.914.
    OpenUrlCrossRefPubMedWeb of Science
  71. 71.↵
    1. Alvarez-Martinez CE,
    2. Christie PJ
    . 2009. Biological diversity of prokaryotic type IV secretion systems. Microbiol Mol Biol Rev 73:775–808. doi:10.1128/MMBR.00023-09.
    OpenUrlAbstract/FREE Full Text
  72. 72.↵
    1. Gandon S,
    2. Hochberg ME,
    3. Holt RD,
    4. Day T
    . 2013. What limits the evolutionary emergence of pathogens? Philos Trans R Soc Lond B Biol Sci 368:20120086. doi:10.1098/rstb.2012.0086.
    OpenUrlCrossRefPubMed
  73. 73.↵
    1. Lo SC,
    2. Hung GC,
    3. Li B,
    4. Lei H,
    5. Li T,
    6. Nagamine K,
    7. Zhang J,
    8. Tsai S,
    9. Bryant R
    . 2013. Isolation of Novel Afipia septicemium and identification of previously unknown bacteria Bradyrhizobium sp. OHSU_III from blood of patients with poorly defined illnesses. PLoS One 8:e0121274. doi:10.1371/journal.pone.0076142.
    OpenUrlCrossRef
  74. 74.↵
    1. Parks DH,
    2. Imelfort M,
    3. Skennerton CT,
    4. Hugenholtz P,
    5. Tyson GW
    . 2015. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res 25:1043–1055. doi:10.1101/gr.186072.114.
    OpenUrlAbstract/FREE Full Text
  75. 75.↵
    1. Hug LA,
    2. Baker BJ,
    3. Anantharaman K,
    4. Brown CT,
    5. Probst AJ,
    6. Castelle CJ,
    7. Butterfield CN,
    8. Hernsdorf AW,
    9. Amano Y,
    10. Ise K,
    11. Suzuki Y,
    12. Dudek N,
    13. Relman DA,
    14. Finstad KM,
    15. Amundson R,
    16. Thomas BC,
    17. Banfield JF
    . 2016. A new view of the tree of life. Nat Microbiol 1:16048. doi:10.1038/nmicrobiol.2016.48.
    OpenUrlCrossRef
  76. 76.↵
    1. Katoh K,
    2. Standley DM
    . 2013. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol 30:772–780. doi:10.1093/molbev/mst010.
    OpenUrlCrossRefPubMedWeb of Science
  77. 77.↵
    1. Capella-Gutiérrez S,
    2. Silla-Martínez JM,
    3. Gabaldón T
    . 2009. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25:1972–1973. doi:10.1093/bioinformatics/btp348.
    OpenUrlCrossRefPubMedWeb of Science
  78. 78.↵
    1. Nguyen LT,
    2. Schmidt HA,
    3. Von Haeseler A,
    4. Minh BQ
    . 2015. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol 32:268–274. doi:10.1093/molbev/msu300.
    OpenUrlCrossRefPubMed
  79. 79.↵
    1. Hoang DT,
    2. Chernomor O,
    3. Von Haeseler A,
    4. Minh BQ,
    5. Vinh LS
    . 2018. UFBoot2: improving the ultrafast bootstrap approximation. Mol Biol Evol 35:518–522. doi:10.1093/molbev/msx281.
    OpenUrlCrossRefPubMed
  80. 80.↵
    1. El Baidouri F,
    2. Venditti C,
    3. Humphries S
    . 2016. Independent evolution of shape and motility allows evolutionary flexibility in Firmicutes bacteria. Nat Ecol Evol 1:9. doi:10.1038/s41559-016-0009.
    OpenUrlCrossRef
  81. 81.↵
    1. Emms DM,
    2. Kelly S
    . 2019. OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol 20:238. doi:10.1186/s13059-019-1832-y.
    OpenUrlCrossRef
  82. 82.↵
    1. Erko S,
    2. Ebers J
    . 2006. Taxonomic parameters revisited: tarnished gold standards. Microbiol Today 33:152–155.
    OpenUrl
  83. 83.↵
    1. Yang Z
    . 2007. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol 24:1586–1591. doi:10.1093/molbev/msm088.
    OpenUrlCrossRefPubMedWeb of Science
  84. 84.↵
    1. Battistuzzi FU,
    2. Feijao A,
    3. Hedges SB
    . 2004. A genomic timescale of prokaryote evolution: insights into the origin of methanogenesis, phototrophy, and the colonization of land. BMC Evol Biol 4:44. doi:10.1186/1471-2148-4-44.
    OpenUrlCrossRefPubMed
  85. 85.↵
    1. Eddy SR
    . 1998. Profile hidden Markov models. Bioinformatics 14:755–763. doi:10.1093/bioinformatics/14.9.755.
    OpenUrlCrossRefPubMedWeb of Science
  86. 86.↵
    1. Li J,
    2. Dohna HZ,
    3. Cardona CJ,
    4. Miller J,
    5. Carpenter TE
    . 2011. Emergence and genetic variation of neuraminidase stalk deletions in avian influenza viruses. PLoS One 6:e14722. doi:10.1371/journal.pone.0014722.
    OpenUrlCrossRefPubMed
  87. 87.↵
    1. Librado P,
    2. Vieira FG,
    3. Rozas J
    . 2012. BadiRate: estimating family turnover rates by likelihood-based methods. Bioinformatics 28:279–281. doi:10.1093/bioinformatics/btr623.
    OpenUrlCrossRefPubMedWeb of Science
  88. 88.↵
    1. David LA,
    2. Alm EJ
    . 2011. Rapid evolutionary innovation during an Archaean genetic expansion. Nature 469:93–96. doi:10.1038/nature09649.
    OpenUrlCrossRefPubMedWeb of Science
  89. 89.↵
    1. Goto N,
    2. Prins P,
    3. Nakao M,
    4. Bonnal R,
    5. Aerts J,
    6. Katayama T
    . 2010. BioRuby: bioinformatics software for the Ruby programming language. Bioinformatics 26:2617–2619. doi:10.1093/bioinformatics/btq475.
    OpenUrlCrossRefPubMedWeb of Science
  90. 90.↵
    1. Shakya M,
    2. Soucy SM,
    3. Zhaxybayeva O
    . 2017. Insights into origin and evolution of α-proteobacterial gene transfer agents. Virus Evol 3:vex036. doi:10.1093/ve/vex036.
    OpenUrlCrossRef
  91. 91.↵
    1. Marrs B
    . 1974. Genetic recombination in Rhodopseudomonas capsulata. Proc Natl Acad Sci U S A 71:971–973. doi:10.1073/pnas.71.3.971.
    OpenUrlAbstract/FREE Full Text
  92. 92.↵
    1. Tamarit D,
    2. Neuvonen MM,
    3. Engel P,
    4. Guy L,
    5. Andersson S
    . 2018. Origin and evolution of the Bartonella gene transfer agent. Mol Biol Evol 35:451–464. doi:10.1093/molbev/msx299.
    OpenUrlCrossRef
  93. 93.↵
    1. Guy L,
    2. Nystedt B,
    3. Toft C,
    4. Zaremba-Niedzwiedzka K,
    5. Berglund EC,
    6. Granberg F,
    7. Näslund K,
    8. Eriksson AS,
    9. Andersson S
    . 2013. A gene transfer agent and a dynamic repertoire of secretion systems hold the keys to the explosive radiation of the emerging pathogen Bartonella. PLoS Genet 9:e1003393. doi:10.1371/journal.pgen.1003393.
    OpenUrlCrossRefPubMed
  94. 94.↵
    1. Persson T,
    2. Battenberg K,
    3. Demina IV,
    4. Vigil-Stenman T,
    5. Vanden Heuvel B,
    6. Pujic P,
    7. Facciotti MT,
    8. Wilbanks EG,
    9. O'Brien A,
    10. Fournier P,
    11. Cruz Hernandez MA,
    12. Mendoza Herrera A,
    13. Médigue C,
    14. Normand P,
    15. Pawlowski K,
    16. Berry AM
    . 2015. Candidatus Frankia datiscae Dg1, the actinobacterial microsymbiont of Datisca glomerata, expresses the canonical nod genes NodABC in symbiosis with its host plant. PLoS One 10:e0127630. doi:10.1371/journal.pone.0127630.
    OpenUrlCrossRef
  95. 95.↵
    1. Dunn OJ
    . 1961. Multiple comparisons among means. J Am Stat Assoc 56:52–64. doi:10.1080/01621459.1961.10482090.
    OpenUrlCrossRefWeb of Science
  96. 96.↵
    1. Benjamini Y,
    2. Hochberg Y
    . 1995. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B 57:289–300. doi:10.1111/j.2517-6161.1995.tb02031.x.
    OpenUrlCrossRef
PreviousNext
Back to top
Download PDF
Citation Tools
Evolutionary Timeline and Genomic Plasticity Underlying the Lifestyle Diversity in Rhizobiales
Sishuo Wang, Andrew Meade, Hon-Ming Lam, Haiwei Luo
mSystems Jul 2020, 5 (4) e00438-20; DOI: 10.1128/mSystems.00438-20

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print
Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mSystems article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Evolutionary Timeline and Genomic Plasticity Underlying the Lifestyle Diversity in Rhizobiales
(Your Name) has forwarded a page to you from mSystems
(Your Name) thought you would be interested in this article in mSystems.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Evolutionary Timeline and Genomic Plasticity Underlying the Lifestyle Diversity in Rhizobiales
Sishuo Wang, Andrew Meade, Hon-Ming Lam, Haiwei Luo
mSystems Jul 2020, 5 (4) e00438-20; DOI: 10.1128/mSystems.00438-20
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS AND DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

Rhizobiales
rhizobia
lifestyle evolution
molecular clock
bacterial evolution
microbial evolution
molecular dating
symbiosis

Related Articles

Cited By...

About

  • About mSystems
  • Author Videos
  • Board of Editors
  • Policies
  • Overleaf Pilot
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Types of Articles
  • Getting Started
  • Ethics
  • Contact Us

Follow #mSystemsJ

@ASMicrobiology

       

 

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2379-5077