Environmentally-driven gene content convergence and the Bacillus phylogeny.

BACKGROUND: Members of the Bacillus genus have been isolated from a variety of environments. However, the relationship between potential metabolism and the niche from which bacteria of this genus have been isolated has not been extensively studied. The existence of a monophyletic aquatic Bacillus group, composed of members isolated from both marine and fresh water has been proposed. Here, we present a phylogenetic/phylogenomic analysis to investigate the potential relationship between the environment from which group members have been isolated and their evolutionary origin. We also carried out hierarchical clustering based on functional content to test for potential environmental effects on the genetic content of these bacteria. RESULTS: The phylogenetic reconstruction showed that Bacillus strains classified as aquatic have evolutionary origins in different lineages. Although we observed the presence of a clade consisting exclusively of aquatic Bacillus, it is not comprised of the same strains previously reported. In contrast to phylogeny, clustering based on the functional categories of the encoded proteomes resulted in groups more compatible with the environments from which the organisms were isolated. This evidence suggests a detectable environmental influence on bacterial genetic content, despite their different evolutionary origins. CONCLUSION: Our results suggest that aquatic Bacillus species have polyphyletic origins, but exhibit convergence at the gene content level.

How To Live with Phosphorus Scarcity in Soil and Sediment: Lessons from Bacteria.

UNLABELLED: Phosphorus (P) plays a fundamental role in the physiology and biochemistry of all living things. Recent evidence indicates that organisms in the oceans can break down and use P forms in different oxidation states (e.g., +5, +3, +1, and -3); however, information is lacking for organisms from soil and sediment. The Cuatro Ciénegas Basin (CCB), Mexico, is an oligotrophic ecosystem with acute P limitation, providing a great opportunity to assess the various strategies that bacteria from soil and sediment use to obtain P. We measured the activities in sediment and soil of different exoenzymes involved in P recycling and evaluated 1,163 bacterial isolates (mainly Bacillus spp.) for their ability to use six different P substrates. DNA turned out to be a preferred substrate, comparable to a more bioavailable P source, potassium phosphate. Phosphodiesterase activity, required for DNA degradation, was observed consistently in the sampled-soil and sediment communities. A capability to use phosphite (PO3 (3-)) and calcium phosphate was observed mainly in sediment isolates. Phosphonates were used at a lower frequency by both soil and sediment isolates, and phosphonatase activity was detected only in soil communities. Our results revealed that soil and sediment bacteria are able to break down and use P forms in different oxidation states and contribute to ecosystem P cycling. Different strategies for P utilization were distributed between and within the different taxonomic lineages analyzed, suggesting a dynamic movement of P utilization traits among bacteria in microbial communities. IMPORTANCE: Phosphorus (P) is an essential element for life found in molecules, such as DNA, cell walls, and in molecules for energy transfer, such as ATP. The Valley of Cuatro Ciénegas, Coahuila (Mexico), is a unique desert characterized by an extreme limitation of P and a great diversity of microbial life. How do bacteria in this valley manage to obtain P? We measured the availability of P and the enzymatic activity associated with P release in soil and sediment. Our results revealed that soil and sediment bacteria can break down and use P forms in different oxidation states and contribute to ecosystem P cycling. Even genetically related bacterial isolates exhibited different preferences for molecules, such as DNA, calcium phosphate, phosphite, and phosphonates, as substrates to obtain P, evidencing a distribution of roles for P utilization and suggesting a dynamic movement of P utilization traits among bacteria in microbial communities.

A three-species synthetic community model whose rapid response to antagonism allows the study of higher-order dynamics and emergent properties in minutes.

Microbial communities can be considered complex adaptive systems. Understanding how these systems arise from different components and how the dynamics of microbial interactions allow for species coexistence are fundamental questions in ecology. To address these questions, we built a three-species synthetic community, called BARS (Bacillota A + S + R). Each species in this community exhibits one of three ecological roles: Antagonistic, Sensitive, or Resistant, assigned in the context of a sediment community. We show that the BARS community reproduces features of complex communities and exhibits higher-order interaction (HOI) dynamics. In paired interactions, the majority of the S species (Sutcliffiella horikoshii 20a) population dies within 5 min when paired with the A species (Bacillus pumilus 145). However, an emergent property appears upon adding the third interactor, as antagonism of species A over S is not observed in the presence of the R species (Bacillus cereus 111). For the paired interaction, within the first 5 min, the surviving population of the S species acquires tolerance to species A, and species A ceases antagonism. This qualitative change reflects endogenous dynamics leading to the expression for tolerance to an antagonistic substance. The stability reached in the triple interaction exhibits a nonlinear response, highly sensitive to the density of the R species. In summary, our HOI model allows the study of the assembly dynamics of a three-species community and evaluating the immediate outcome within a 30 min frame. The BARS has features of a complex system where the paired interactions do not predict the community dynamics. The model is amenable to mechanistic dissection and to modeling how the parts integrate to achieve collective properties.

Intermunicipal travel networks of Mexico during the COVID-19 pandemic.

Human mobility networks are widely used for diverse studies in geography, sociology, and economics. In these networks, nodes usually represent places or regions and links refer to movement between them. They become essential when studying the spread of a virus, the planning of transit, or society's local and global structures. Therefore, the construction and analysis of human mobility networks are crucial for a vast number of real-life applications. This work presents a collection of networks that describe the human travel patterns between municipalities in Mexico in the 2020-2021 period. Using anonymized mobile location data, we constructed directed, weighted networks representing the volume of travels between municipalities. We analysed changes in global, local, and mesoscale network features. We observe that changes in these features are associated with factors such as COVID-19 restrictions and population size. In general, the implementation of restrictions at the start of the COVID-19 pandemic in early 2020, induced more intense changes in network features than later events, which had a less notable impact in network features. These networks will result very useful for researchers and decision-makers in the areas of transportation, infrastructure planning, epidemic control and network science at large.

A Transiently Hypersaline Microbial Mat Harbors a Diverse and Stable Archaeal Community in the Cuatro Cienegas Basin, Mexico.

Microbial mats are biologically diverse communities that are analogs to some of the earliest ecosystems on Earth. In this study, we describe a unique transiently hypersaline microbial mat uncovered in a shallow pond within the Cuatro Cienegas Basin (CCB) in northern México. The CCB is an endemism-rich site that harbors living stromatolites that have been studied to understand the conditions of the Precambrian Earth. These microbial mats form elastic domes filled with biogenic gas, and the mats have a relatively large and stable subpopulation of archaea. For this reason, this site has been termed archaean domes (AD). The AD microbial community was analyzed by metagenomics over three seasons. The mat exhibited a highly diverse prokaryotic community dominated by bacteria. Bacterial sequences are represented in 37 phyla, mainly Proteobacteria, Firmicutes, and Actinobacteria, that together comprised >50% of the sequences from the mat. Archaea represented up to 5% of the retrieved sequences, with up to 230 different archaeal species that belong to 5 phyla (Euryarchaeota, Crenarchaeota, Thaumarchaeota, Korarchaeota, and Nanoarchaeota). The archaeal taxa showed low variation despite fluctuations in water and nutrient availability. In addition, predicted functions highlight stress responses to extreme conditions present in the AD, including salinity, pH, and water/drought fluctuation. The observed complexity of the AD mat thriving in high pH and fluctuating water and salt conditions within the CCB provides an extant model of great value for evolutionary studies, as well as a suitable analog to the early Earth and Mars.

Comparing the evolutionary dynamics of predominant SARS-CoV-2 virus lineages co-circulating in Mexico.

Over 200 different SARS-CoV-2 lineages have been observed in Mexico by November 2021. To investigate lineage replacement dynamics, we applied a phylodynamic approach and explored the evolutionary trajectories of five dominant lineages that circulated during the first year of local transmission. For most lineages, peaks in sampling frequencies coincided with different epidemiological waves of infection in Mexico. Lineages B.1.1.222 and B.1.1.519 exhibited similar dynamics, constituting clades that likely originated in Mexico and persisted for >12 months. Lineages B.1.1.7, P.1 and B.1.617.2 also displayed similar dynamics, characterized by multiple introduction events leading to a few successful extended local transmission chains that persisted for several months. For the largest B.1.617.2 clades, we further explored viral lineage movements across Mexico. Many clades were located within the south region of the country, suggesting that this area played a key role in the spread of SARS-CoV-2 in Mexico.

Phylogenetic distribution and evolutionary history of bacterial DEAD-Box proteins.

DEAD-box proteins are found in all domains of life and participate in almost all cellular processes that involve RNA. The presence of DEAD and Helicase_C conserved domains distinguish these proteins. DEAD-box proteins exhibit RNA-dependent ATPase activity in vitro, and several also show RNA helicase activity. In this study, we analyzed the distribution and architecture of DEAD-box proteins among bacterial genomes to gain insight into the evolutionary pathways that have shaped their history. We identified 1,848 unique DEAD-box proteins from 563 bacterial genomes. Bacterial genomes can possess a single copy DEAD-box gene, or up to 12 copies of the gene, such as in Shewanella. The alignment of 1,208 sequences allowed us to perform a robust analysis of the hallmark motifs of DEAD-box proteins and determine the residues that occur at high frequency, some of which were previously overlooked. Bacterial DEAD-box proteins do not generally contain a conserved C-terminal domain, with the exception of some members that possess a DbpA RNA-binding domain (RBD). Phylogenetic analysis showed a separation of DbpA-RBD-containing and DbpA-RBD-lacking sequences and revealed a group of DEAD-box protein genes that expanded mainly in the Proteobacteria. Analysis of DEAD-box proteins from Firmicutes and γ-Proteobacteria, was used to deduce orthologous relationships of the well-studied DEAD-box proteins from Escherichia coli and Bacillus subtilis. These analyses suggest that DbpA-RBD is an ancestral domain that most likely emerged as a specialized domain of the RNA-dependent ATPases. Moreover, these data revealed numerous events of gene family expansion and reduction following speciation.

Phenotypic plasticity and evolution of thermal tolerance in bacteria from temperate and hot spring environments.

Phenotypic plasticity allows individuals to respond to the selective forces of a new environment, followed by adaptive evolution. We do not know to what extent phenotypic plasticity allows thermal tolerance evolution in bacteria at the border of their physiological limits. We analyzed growth and reaction norms to temperature of strains of two bacterial lineages, Bacillus cereus sensu lato and Bacillus subtilis sensu lato, that evolved in two contrasting environments, a temperate lagoon (T) and a hot spring (H). Our results showed that despite the co-occurrence of members of both lineages in the two contrasting environments, norms of reactions to temperature exhibited a similar pattern only in strains within the lineages, suggesting fixed phenotypic plasticity. Additionally, strains from the H environment showed only two to three degrees centigrade more heat tolerance than strains from the T environment. Their viability decreased at temperatures above their optimal for growth, particularly for the B. cereus lineage. However, sporulation occurred at all temperatures, consistent with the known cell population heterogeneity that allows the Bacillus to anticipate adversity. We suggest that these mesophilic strains survive in the hot-spring as spores and complete their life cycle of germination and growth during intermittent opportunities of moderate temperatures. The limited evolutionary changes towards an increase in heat tolerance in bacteria should alert us of the negative impact of climate change on all biological cycles in the planet, which at its most basic level depends on microorganisms.

Characterization of tRNA(Cys) processing in a conditional Bacillus subtilis CCase mutant reveals the participation of RNase R in its quality control.

We generated a conditional CCase mutant of Bacillus subtilis to explore the participation in vivo of the tRNA nucleotidyltransferase (CCA transferase or CCase) in the maturation of the single-copy tRNA(Cys), which lacks an encoded CCA 3' end. We observed that shorter tRNA(Cys) species, presumably lacking CCA, only accumulated when the inducible Pspac : cca was introduced into an rnr mutant strain, but not in combination with pnp. We sequenced the tRNA 3' ends produced in the various mutant tRNA(Cys) species to detect maturation and decay intermediates and observed that decay of the tRNA(Cys) occurs through the addition of poly(A) or heteropolymeric tails. A few clones corresponding to full-size tRNAs contained either CCA or other C and/or A sequences, suggesting that these are substrates for repair and/or decay. We also observed editing of tRNA(Cys) at position 21, which seems to occur preferentially in mature tRNAs. Altogether, our results provide in vivo evidence for the participation of the B. subtilis cca gene product in the maturation of tRNAs lacking CCA. We also suggest that RNase R exoRNase in B. subtilis participates in the quality control of tRNA.

Life on Phosphite: A Metagenomics Tale.

Phosphite, a species of phosphorus in a P3+ oxidation state, is believed to have played an important role in the primordial Earth. Figueroa et al. used metagenomics to uncover anaerobic bacterial communities from waste water waste sludge that sustain life from energy provided by phosphite.

Correction for Zarza et al., "Complete Genome Sequences of Two Bacillus pumilus Strains from Cuatrociénegas, Coahuila, Mexico".

[This corrects the article DOI: 10.1128/genomeA.00364-18.].

Complete Genome Sequences of Two Bacillus pumilus Strains from Cuatrociénegas, Coahuila, Mexico.

We assembled the complete genome sequences of Bacillus pumilus strains 145 and 150a from Cuatrociénegas, Mexico. We detected genes codifying for proteins potentially involved in antagonism (bacteriocins) and defense mechanisms (abortive infection bacteriophage proteins and 4-azaleucine resistance). Both strains harbored prophage sequences. Our results provide insights into understanding the establishment of microbial interactions.

The DEAD-Box RNA Helicases of Bacillus subtilis as a Model to Evaluate Genetic Compensation Among Duplicate Genes.

The presence of duplicated genes in organisms is well documented. There is increasing interest in understanding how these genes subfunctionalize and whether functional overlap can explain the fact that some of these genes are dispensable. Bacillus subtilis possesses four DEAD-box RNA helicases (DBRH) genes, cshA, cshB, deaD/yxiN, and yfmL that make a good case to study to what extent they can complement each other despite their subfunctionalization. They possess the highly conserved N-terminal catalytic domain core common to RNA helicases, but different carboxy-terminal ends. All four genes have been shown to have independent functions although all participate in rRNA assembly. None of the B. subtilis DBRH is essential for growth at 37°C, and all single deletion mutants exhibit defective growth at 18°C except for ΔdeaD/yxiN. Evaluation of double mutants did not reveal negative epistasis, suggesting that they do not have overlapping functions. The absence of any one gene distorts the expression pattern of the others, but not in a specific pattern suggestive of compensation. Overexpression of these paralogous genes in the different mutant backgrounds did not result in cross-complementation, further confirming their lack of buffering capability. Since no complementation could be observed among full sized proteins, we evaluated to what extent the superfamily 2 (SF2) helicase core of the smallest DBRH, YfmL, could be functional when hooked to each of the C-terminal end of CshA, CshB, and DeaD/YxiN. None of the different chimeras complemented the different mutants, and instead, all chimeras inhibited the growth of the ΔyfmL mutant, and other combinations were also deleterious. Our findings suggest that the long time divergence between DEAD-box RNA helicase genes has resulted in specialized activities in RNA metabolism and shows that these duplicated genes cannot buffer one another.

Complete Genome Sequence of Bacillus horikoshii Strain 20a from Cuatro Cienegas, Coahuila, Mexico.

We sequenced the Bacillus horikoshii 20a genome, isolated from sediment collected in Cuatro Cienegas, Mexico. We identified genes involved in establishing antagonistic interactions in microbial communities (antibiotic resistance and bacteriocins) and genes related to the metabolism of cyanophycin, a reserve compound and spore matrix material potentially relevant for survival in an oligotrophic environment.

Phenotypic Microdiversity and Phylogenetic Signal Analysis of Traits Related to Social Interaction in Bacillus spp. from Sediment Communities.

Understanding the relationship between phylogeny and predicted traits is important to uncover the dimension of the predictive power of a microbial composition approach. Numerous works have addressed the taxonomic composition of bacteria in communities, but little is known about trait heterogeneity in closely related bacteria that co-occur in communities. We evaluated a sample of 467 isolates from the Churince water system of the Cuatro Cienegas Basin (CCB), enriched for Bacillus spp. The 16S rRNA gene revealed a random distribution of taxonomic groups within this genus among 11 sampling sites. A subsample of 141 Bacillus spp. isolates from sediment, with seven well-represented species was chosen to evaluate the heterogeneity and the phylogenetic signal of phenotypic traits that are known to diverge within small clades, such as substrate utilization, and traits that are conserved deep in the lineage, such as prototrophy, swarming and biofilm formation. We were especially interested in evaluating social traits, such as swarming and biofilm formation, for which cooperation is needed to accomplish a multicellular behavior and for which there is little information from natural communities. The phylogenetic distribution of traits, evaluated by the Purvis and Fritz's D statistics approached a Brownian model of evolution. Analysis of the phylogenetic relatedness of the clusters of members sharing the trait using consenTRAIT algorithm, revealed more clustering and deeper phylogenetic signal for prototrophy, biofilm and swimming compared to the data obtained for substrate utilization. The explanation to the observed Brownian evolution of social traits could be either loss due to complete dispensability or to compensated trait loss due to the availability of public goods. Since many of the evaluated traits can be considered to be collective action traits, such as swarming, motility and biofilm formation, the observed microdiversity within taxonomic groups might be explained by distributed functions in structured communities.

The Genome Sequence of Avibacterium paragallinarum Strain CL Has a Large Repertoire of Insertion Sequence Elements.

The draft genome sequence of Avibacterium paragallinarum strain CL serovar C is reported here. The genome comprises 154 contigs corresponding to 2.4 Mb with 41% G+C content and many insertion sequence (IS) elements, a characteristic not previously reported in A. paragallinarum.

Draft Whole-Genome Sequence of the Type Strain Bacillus aquimaris TF12T.

Bacillus aquimaris TF12 is a Gram-positive bacteria isolated from a tidal flat of the Yellow Sea in South Korea. We report the draft whole-genome sequence of Bacillus aquimaris TF12, the type strain of a set of bacteria typically associated with marine habitats and with a potentially high biotechnology value.

Draft Whole-Genome Sequence of the Type Strain Bacillus horikoshii DSM 8719.

Members of the Bacillus genus have been extensively studied because of their ability to produce enzymes with high biotechnological value. Here, we report the draft of the whole-genome sequence of the type strain Bacillus horikoshii DSM 8719, an alkali-tolerant strain.

Microevolution Analysis of Bacillus coahuilensis Unveils Differences in Phosphorus Acquisition Strategies and Their Regulation.

Bacterial genomes undergo numerous events of gene losses and gains that generate genome variability among strains of the same species (microevolution). Our aim was to compare the genomes and relevant phenotypes of three Bacillus coahuilensis strains from two oligotrophic hydrological systems in the Cuatro Ciénegas Basin (México), to unveil the environmental challenges that this species cope with, and the microevolutionary differences in these genotypes. Since the strains were isolated from a low P environment, we placed emphasis on the search of different phosphorus acquisition strategies. The three B. coahuilensis strains exhibited similar numbers of coding DNA sequences, of which 82% (2,893) constituted the core genome, and 18% corresponded to accessory genes. Most of the genes in this last group were associated with mobile genetic elements (MGEs) or were annotated as hypothetical proteins. Ten percent of the pangenome consisted of strain-specific genes. Alignment of the three B. coahuilensis genomes indicated a high level of synteny and revealed the presence of several genomic islands. Unexpectedly, one of these islands contained genes that encode the 2-keto-3-deoxymannooctulosonic acid (Kdo) biosynthesis enzymes, a feature associated to cell walls of Gram-negative bacteria. Some microevolutionary changes were clearly associated with MGEs. Our analysis revealed inconsistencies between phenotype and genotype, which we suggest result from the impossibility to map regulatory features to genome analysis. Experimental results revealed variability in the types and numbers of auxotrophies between the strains that could not consistently be explained by in silico metabolic models. Several intraspecific differences in preferences for carbohydrate and phosphorus utilization were observed. Regarding phosphorus recycling, scavenging, and storage, variations were found between the three genomes. The three strains exhibited differences regarding alkaline phosphatase that revealed that in addition to gene gain and loss, regulation adjustment of gene expression also has contributed to the intraspecific diversity of B. coahuilensis.

Antagonistic interactions are sufficient to explain self-assemblage of bacterial communities in a homogeneous environment: a computational modeling approach.

Most of the studies in Ecology have been devoted to analyzing the effects the environment has on individuals, populations, and communities, thus neglecting the effects of biotic interactions on the system dynamics. In the present work we study the structure of bacterial communities in the oligotrophic shallow water system of Churince, Cuatro Cienegas, Mexico. Since the physicochemical conditions of this water system are homogeneous and quite stable in time, it is an excellent candidate to study how biotic factors influence the structure of bacterial communities. In a previous study, the binary antagonistic interactions of 78 bacterial strains, isolated from Churince, were experimentally determined. We employ these data to develop a computer algorithm to simulate growth experiments in a cellular grid representing the pond. Remarkably, in our model, the dynamics of all the simulated bacterial populations is determined solely by antagonistic interactions. Our results indicate that all bacterial strains (even those that are antagonized by many other bacteria) survive in the long term, and that the underlying mechanism is the formation of bacterial community patches. Patches corresponding to less antagonistic and highly susceptible strains are consistently isolated from the highly-antagonistic bacterial colonies by patches of neutral strains. These results concur with the observed features of the bacterial community structure previously reported. Finally, we study how our findings depend on factors like initial population size, differential population growth rates, homogeneous population death rates, and enhanced bacterial diffusion.