Search

What a strain is and how many strains make up a natural bacterial population remain elusive concepts despite their apparent importance for assessing the role of intra-population diversity in disease emergence or response to environmental perturbations. To advance these concepts, we sequenced 138 randomly selected Salinibacter ruber isolates from two solar salterns and assessed these genomes against companion short-read metagenomes from the same samples. The distribution of genome-aggregate average nucleotide identity (ANI) values among these isolates revealed a bimodal distribution, with four-fold lower occurrence of values between 99.2% and 99.8% relative to ANI >99.8% or <99.2%, revealing a natural “gap” in the sequence space within species. Accordingly, we used this ANI gap to define genomovars and a higher ANI value of >99.99% and shared gene-content >99.0% to define strains. Using these thresholds and extrapolating from how many metagenomic reads each genomovar uniquely recruited, we estimated that –although our 138 isolates represented about 80% of the Sal. ruber population– the total population in one saltern pond is composed of 5,500 to 11,000 genomovars, the great majority of which appear to be rare in-situ. These data also revealed that the most frequently recovered isolate in lab media was often not the most abundant genomovar in-situ, suggesting that cultivation biases are significant, even in cases that cultivation procedures are thought to be robust. The methodology and ANI thresholds outlined here should represent a useful guide for future microdiversity surveys of additional microbial species.

Large-scale surveys of prokaryotic communities (metagenomes), as well as isolate genomes, have revealed that their diversity is predominantly organized in sequence-discrete units that may be equated to species. Specifically, genomes of the same species commonly show genome-aggregate average nucleotide identity (ANI) >95% among themselves and ANI <90% to members of other species, while genomes showing ANI 90%–95% are comparatively rare. However, it remains unclear if such “discontinuities” or gaps in ANI values can be observed within species and thus used to advance and standardize intra-species units. By analyzing 18,123 complete isolate genomes from 330 bacterial species with at least 10 genome representatives each and available long-read metagenomes, we show that another discontinuity exists between 99.2% and 99.8% (midpoint 99.5%) ANI in most of these species. The 99.5% ANI threshold is largely consistent with how sequence types have been defined in previous epidemiological studies but provides clusters with ~20% higher accuracy in terms of evolutionary and gene-content relatedness of the grouped genomes, while strains should be consequently defined at higher ANI values (>99.99% proposed). Collectively, our results should facilitate future micro-diversity studies across clinical or environmental settings because they provide a more natural definition of intra-species units of diversity., With funding from the Spanish government through the ‘María de Maeztu Unit of Excelence’ accreditation (CEX2021-001198-M).

Despite the importance of intra-species variants of viruses for causing disease and/or disrupting ecosystem functioning, there is no universally applicable standard to define these. A 95% whole-genome average nucleotide identity (ANI) gap is commonly used to define species, especially for bacteriophages, but whether a similar gap exists within species that can be used to define intra-species units has not been evaluated yet. Whole-genome comparisons among members of 1,016 bacteriophage species revealed a region of low frequency of pairs around 99.2-99.8% ANI, showing 3-fold or fewer pairs than expected for an even or normal distribution. This second gap is prevalent in viruses infecting various cultured or uncultured hosts, and from a variety of environments, although a few exceptions to this pattern were also observed (∼3.7% of the total species evaluated) and are likely attributed to cultivation biases. Similar results were observed for a limited set of eukaryotic viruses that are adequately sampled including SARS-CoV-2, whose ANI-based clusters matched well the WHO-defined Variants of Concern, indicating that they represent functionally and/or ecologically distinct units. The existence of sequence-discrete units appears to be predominantly driven by (high) ecological cohesiveness coupled to either recombination frequency for bacteriophages or selection and clonal evolution for other viruses such as SARS-CoV-2. These results indicate that fundamentally different underlying mechanisms could lead to similar diversity patterns. Based on these results, we propose the 99.5% ANI as a practical, standardized, and data-supported threshold for defining viral intra-species units of bacteriophages, for which we propose the term genomovars.

Current -omics methods allow the collection of a large amount of information that helps in describing the microbial diversity in nature. Here, and as a result of a culturomic approach that rendered the collection of thousands of isolates from 5 different hypersaline sites (in Spain, USA and New Zealand), we obtained 21 strains that represent two new Salinibacter species. For these species we propose the names Salinibacter pepae sp. nov. and Salinibacter grassmerensis sp. nov. (showing average nucleotide identity (ANI) values < 95.09% and 87.08% with Sal. ruber M31T, respectively). Metabolomics revealed species-specific discriminative profiles. Sal. ruber strains were distinguished by a higher percentage of polyunsaturated fatty acids and specific N-functionalized fatty acids; and Sal. altiplanensis was distinguished by an increased number of glycosylated molecules. Based on sequence characteristics and inferred phenotype of metagenome-assembled genomes (MAGs), we describe two new members of the genus Salinibacter. These species dominated in different sites and always coexisted with Sal. ruber and Sal. pepae. Based on the MAGs from three Argentinian lakes in the Pampa region of Argentina and the MAG of the Romanian lake Fără Fund, we describe the species Salinibacter pampae sp. nov. and Salinibacter abyssi sp. nov. respectively (showing ANI values 90.94% and 91.48% with Sal. ruber M31T, respectively). Sal. grassmerensis sp. nov. name was formed according to the rules of the International Code for Nomenclature of Prokaryotes (ICNP), and Sal. pepae, Sal. pampae sp. nov. and Sal. abyssi sp. nov. are proposed following the rules of the newly published Code of Nomenclature of Prokaryotes Described from Sequence Data (SeqCode). This work constitutes an example on how classification under ICNP and SeqCode can coexist, and how the official naming a cultivated organism for which the deposit in public repositories is difficult finds an intermediate solution.

Current -omics methods allow the collection of a large amount of information that helps in describing the microbial diversity in nature. Here, and as a result of a culturomic approach that rendered the collection of thousands of isolates from 5 different hypersaline sites (in Spain, USA and New Zealand), we obtained 21 strains that represent two new Salinibacter species. For these species we propose the names Salinibacter pepae sp. nov. and Salinibacter grassmerensis sp. nov. (showing average nucleotide identity (ANI) values < 95.09% and 87.08% with Sal. ruber M31T, respectively). Metabolomics revealed species-specific discriminative profiles. Sal. ruber strains were distinguished by a higher percentage of polyunsaturated fatty acids and specific N-functionalized fatty acids; and Sal. altiplanensis was distinguished by an increased number of glycosylated molecules. Based on sequence characteristics and inferred phenotype of metagenome-assembled genomes (MAGs), we describe two new members of the genus Salinibacter. These species dominated in different sites and always coexisted with Sal. ruber and Sal. pepae. Based on the MAGs from three Argentinian lakes in the Pampa region of Argentina and the MAG of the Romanian lake Fără Fund, we describe the species Salinibacter pampae sp. nov. and Salinibacter abyssi sp. nov. respectively (showing ANI values 90.94% and 91.48% with Sal. ruber M31T, respectively). Sal. grassmerensis sp. nov. name was formed according to the rules of the International Code for Nomenclature of Prokaryotes (ICNP), and Sal. pepae, Sal. pampae sp. nov. and Sal. abyssi sp. nov. are proposed following the rules of the newly published Code of Nomenclature of Prokaryotes Described from Sequence Data (SeqCode). This work constitutes an example on how classification under ICNP and SeqCode can coexist, and how the official naming a cultivated organism for which the deposit in public repositories is difficult finds an intermediate solution.

Figure S1. Geographical distribution of the six locations studied here and their inter-location distances (given in the table). The distances given in each location refer to Mallorca. The hypersaline here are: Mallorca (Es Trenc and S’Avall and separated by 3.5 km) and Santa Pola in Spain, Fără Fund in Romania, Great Salt Lake (USA), Pampa (Laguna Colorada Chica and Laguna Colorada Grande, separated by 23 Km), and Lake Grassmere (New Zealand). Table S1. Statistics of Salinibacter MAGs recovered using different bioinformatic tools. MAGs selected for further analysis and taxonomic classification are marked in red. Text S1: In addition to the taxonomic study, here we assessed the different tools SPADES and MegaHIT assemblers, and MetaBAT binning tools selecting contigs with length > 2,000 pbs or > 5,000 pbs to retrieve the best quality MAG (Sup. Table S1). We observed that different tools and combinations rendered different qualities of the MAGs. The selection of contigs with length > 5,000 pbs binning rendered MAGs with a lower genome size, lower completeness and similar contamination. From the FF metagenome (Romania), the highest quality MAG was retrieved using the assembler SPADES and MetaBAT binning tool using a minimum contig length of 2,000 pbs. From CCH and CG metagenomes (Argentina), the highest quality MAGs were retrieved using the assembler MegaHIT and MetaBAT binning tool using a minimum contig length of 2,000 pbs (Sup. Table S1). Figure S2. Genomic clustering based on Average Amino-acid Identity (AAI) of genomes included in the study. Figure S3. Pangenome hierarchical clustering based on the presence (grey) or absence (black) of orthologous genes. Table S2. Pangenome statistics between genomes belonging to the same species. Table S3: CRISPR-Cas systems found in the Salinibacter pepae and Salinibacter grassmerensis genomes. Table S4. Cell morphologies of the new isolates. All cells have been cultivated onto MA agar. for 10 days at 12ºC. The morphology was ob

Nonrecombining sex chromosomes, like the mammalian Y, often lose genes and accumulate transposable ele ments, a process termed degeneration. The correlation between suppressed recombination and degeneration is clear in animal XY systems, but the absence of recombination is confounded with other asymmetries between the X and Y. In contrast, UV sex chromosomes, like those found in bryophytes, experience symmetrical population genetic conditions. Here, we generate nearly gapless female and male chromosome-scale reference genomes of the moss Ceratodon purpureus to test for degeneration in the bryophyte UV sex chromosomes. We show that the moss sex chromosomes evolved over 300 million years ago and expanded via two chromosomal fusions. Although the sex chromosomes exhibit weaker purifying selection than autosomes, we find that suppressed recombination alone is insufficient to drive degeneration. Instead, the U and V sex chromosomes harbor thousands of broadly expressed genes, including numerous key regulators of sexual development across land plants.

Microbial communities often harbor overwhelming species and gene diversity, making it challenging to determine the important units to study this diversity. We argue that the reduced, and thus tractable, microbial diversity of manmade salterns provides an ideal system to advance this cornerstone issue. We review recent time-series genomic and metagenomic studies of the saltern-dominating bacterial and archaeal taxa to show that these taxa form persistent, sequence-discrete, species-like populations. While these populations harbor extensive intra-population gene diversity, even within a single saltern site, only a small minority of these genes appear to be functionally important during environmental perturbations. We outline an approach to detect and track such populations and their ecologically important genes that should be broadly applicable.

Metagenomic surveys have revealed that natural microbial communities are predominantly composed of sequence-discrete, species-like populations but the genetic and/or ecological processes that maintain such populations remain speculative, limiting our understanding of population speciation and adaptation to perturbations. To address this knowledge gap, we sequenced 112 Salinibacter ruber isolates and 12 companion metagenomes from four adjacent saltern ponds in Mallorca, Spain that were experimentally manipulated to dramatically alter salinity and light intensity, the two major drivers of this ecosystem. Our analyses showed that the pangenome of the local Sal. ruber population is open and similar in size (~15,000 genes) to that of randomly sampled Escherichia coli genomes. While most of the accessory (noncore) genes were isolate-specific and showed low in situ abundances based on the metagenomes compared to the core genes, indicating that they were functionally unimportant and/or transient, 3.5% of them became abundant when salinity (but not light) conditions changed and encoded for functions related to osmoregulation. Nonetheless, the ecological advantage of these genes, while significant, was apparently not strong enough to purge diversity within the population. Collectively, our results provide an explanation for how this immense intrapopulation gene diversity is maintained, which has implications for the prokaryotic species concept.

Microbial communities thriving in hypersaline brines of solar salterns are highly resistant and resilient to environmental changes, and salinity is a major factor that deterministically influences community structure. Here, we demonstrate that this resilience occurs even after rapid osmotic shocks caused by a threefold change in salinity (a reduction from 34 to 12% salts) leading to massive amounts of archaeal cell lysis. Specifically, our temporal metagenomic datasets identified two co-occurring ecotypes within the most dominant archaeal population of the brines Haloquadratum walsbyi that exhibited different salt concentration preferences. The dominant ecotype was generally more abundant and occurred in high-salt conditions (34%); the low abundance ecotype always co-occurred but was enriched at salinities around 20% or lower and carried unique gene content related to solute transport and gene regulation. Despite their apparent distinct ecological preferences, the ecotypes did not outcompete each other presumably due to weak functional differentiation between them. Further, the osmotic shock selected for a temporal increase in taxonomic and functional diversity at both the Hqr. walsbyi population and whole-community levels supporting the specialization-disturbance hypothesis, that is, the expectation that disturbance favors generalists. Altogether, our results provide new insights into how intraspecies diversity is maintained in light of substantial gene-content differences and major environmental perturbations.

Microbial communities thriving in hypersaline brines of solar salterns are highly resistant and resilient to environmental changes, and salinity is a major factor that deterministically influences community structure. Here, we demonstrate that this resilience occurs even after rapid osmotic shocks caused by a threefold change in salinity (a reduction from 34 to 12% salts) leading to massive amounts of archaeal cell lysis. Specifically, our temporal metagenomic datasets identified two co-occurring ecotypes within the most dominant archaeal population of the brines Haloquadratum walsbyi that exhibited different salt concentration preferences. The dominant ecotype was generally more abundant and occurred in high-salt conditions (34%); the low abundance ecotype always co-occurred but was enriched at salinities around 20% or lower and carried unique gene content related to solute transport and gene regulation. Despite their apparent distinct ecological preferences, the ecotypes did not outcompete each other presumably due to weak functional differentiation between them. Further, the osmotic shock selected for a temporal increase in taxonomic and functional diversity at both the Hqr. walsbyi population and whole-community levels supporting the specialization-disturbance hypothesis, that is, the expectation that disturbance favors generalists. Altogether, our results provide new insights into how intraspecies diversity is maintained in light of substantial gene-content differences and major environmental perturbations.