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.

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.

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.

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.

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.].

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Variability of rRNA Operon Copy Number and Growth Rate Dynamics of Bacillus Isolated from an Extremely Oligotrophic Aquatic Ecosystem.

The ribosomal RNA (rrn) operon is a key suite of genes related to the production of protein synthesis machinery and thus to bacterial growth physiology. Experimental evidence has suggested an intrinsic relationship between the number of copies of this operon and environmental resource availability, especially the availability of phosphorus (P), because bacteria that live in oligotrophic ecosystems usually have few rrn operons and a slow growth rate. The Cuatro Ciénegas Basin (CCB) is a complex aquatic ecosystem that contains an unusually high microbial diversity that is able to persist under highly oligotrophic conditions. These environmental conditions impose a variety of strong selective pressures that shape the genome dynamics of their inhabitants. The genus Bacillus is one of the most abundant cultivable bacterial groups in the CCB and usually possesses a relatively large number of rrn operon copies (6-15 copies). The main goal of this study was to analyze the variation in the number of rrn operon copies of Bacillus in the CCB and to assess their growth-related properties as well as their stoichiometric balance (N and P content). We defined 18 phylogenetic groups within the Bacilli clade and documented a range of from six to 14 copies of the rrn operon. The growth dynamic of these Bacilli was heterogeneous and did not show a direct relation to the number of operon copies. Physiologically, our results were not consistent with the Growth Rate Hypothesis, since the copies of the rrn operon were decoupled from growth rate. However, we speculate that the diversity of the growth properties of these Bacilli as well as the low P content of their cells in an ample range of rrn copy number is an adaptive response to oligotrophy of the CCB and could represent an ecological mechanism that allows these taxa to coexist. These findings increase the knowledge of the variability in the number of copies of the rrn operon in the genus Bacillus and give insights about the physiology of this bacterial group under extreme oligotrophic conditions.

Antagonism influences assembly of a Bacillus guild in a local community and is depicted as a food-chain network.

Understanding the principles that govern community assemblages is a central goal of ecology. There is limited experimental evidence in natural settings showing that microbial assembly in communities are influenced by antagonistic interactions. We, therefore, analyzed antagonism among bacterial isolates from a taxonomically related bacterial guild obtained from five sites in sediments from a fresh water system. We hypothesized that if antagonistic interactions acted as a shaping force of the community assembly, then the frequency of resistance to antagonism among bacterial isolates originating from a given site would be higher than the resistance to conspecifics originating from a different assemblage. Antagonism assays were conducted between 78 thermoresistant isolates, of which 72 were Bacillus spp. Sensitive, resistant and antagonistic isolates co-occurred at each site, but the within-site frequency of resistance observed was higher than that observed when assessed across-sites. We found that antagonism results from bacteriocin-like substances aimed at the exclusion of conspecifics. More than 6000 interactions were scored and described by a directed network with hierarchical structure that exhibited properties that resembled a food chain, where the different Bacillus taxonomic groups occupied specific positions. For some tested interacting pairs, the unidirectional interaction could be explained by competition that inhibited growth or completely excluded one of the pair members. This is the first report on the prevalence and specificity of Bacillus interactions in a natural setting and provides evidence for the influence of bacterial antagonist interactions in the assemblage of a taxonomically related guild in local communities.