Single-amplified genomes reveal most streamlined free-living marine bacteria.

Evolutionary adaptations of prokaryotes to the environment sometimes result in genome reduction. Our knowledge of this phenomenon among free-living bacteria remains scarce. We address the dynamics and limits of genome reduction by examining one of the most abundant bacteria in the ocean, the SAR86 clade. Despite its abundance, comparative genomics has been limited by the absence of pure cultures and the poor representation in metagenome-assembled genomes. We co-assembled multiple previously available single-amplified genomes to obtain the first complete genomes from members of the four families. All families showed a convergent evolutionary trajectory with characteristic features of streamlined genomes, most pronounced in the TMED112 family. This family has a genome size of ca. 1 Mb and only 1 bp as median intergenic distance, exceeding values found in other abundant microbes such as SAR11, OM43 and Prochlorococcus. This genomic simplification led to a reduction in the biosynthesis of essential molecules, DNA repair-related genes, and the ability to sense and respond to environmental factors, which could suggest an evolutionary dependence on other co-occurring microbes for survival (Black Queen hypothesis). Therefore, these reconstructed genomes within the SAR86 clade provide new insights into the limits of genome reduction in free-living marine bacteria.

Ecogenomics of the SAR11 clade.

Members of the SAR11 clade, despite their high abundance, are often poorly represented by metagenome-assembled genomes. This fact has hampered our knowledge about their ecology and genetic diversity. Here we examined 175 SAR11 genomes, including 47 new single-amplified genomes. The presence of the first genomes associated with subclade IV suggests that, in the same way as subclade V, they might be outside the proposed Pelagibacterales order. An expanded phylogenomic classification together with patterns of metagenomic recruitment at a global scale have allowed us to define new ecogenomic units of classification (genomospecies), appearing at different, and sometimes restricted, metagenomic data sets. We detected greater microdiversity across the water column at a single location than in samples collected from similar depth across the global ocean, suggesting little influence of biogeography. In addition, pangenome analysis revealed that the flexible genome was essential to shape genomospecies distribution. In one genomospecies preferentially found within the Mediterranean, a set of genes involved in phosphonate utilization was detected. While another, with a more cosmopolitan distribution, was unique in having an aerobic purine degradation pathway. Together, these results provide a glimpse of the enormous genomic diversity within this clade at a finer resolution than the currently defined clades.

Genome diversity of marine phages recovered from Mediterranean metagenomes: Size matters.

Marine viruses play a critical role not only in the global geochemical cycles but also in the biology and evolution of their hosts. Despite their importance, viral diversity remains underexplored mostly due to sampling and cultivation challenges. Direct sequencing approaches such as viromics has provided new insights into the marine viral world. As a complementary approach, we analysed 24 microbial metagenomes (>0.2 µm size range) obtained from six sites in the Mediterranean Sea that vary by depth, season and filter used to retrieve the fraction. Filter-size comparison showed a significant number of viral sequences that were retained on the larger-pore filters and were different from those found in the viral fraction from the same sample, indicating that some important viral information is missing using only assembly from viromes. Besides, we were able to describe 1,323 viral genomic fragments that were more than 10Kb in length, of which 36 represented complete viral genomes including some of them retrieved from a cross-assembly from different metagenomes. Host prediction based on sequence methods revealed new phage groups belonging to marine prokaryotes like SAR11, Cyanobacteria or SAR116. We also identified the first complete virophage from deep seawater and a new endemic clade of the recently discovered Marine group II Euryarchaeota virus. Furthermore, analysis of viral distribution using metagenomes and viromes indicated that most of the new phages were found exclusively in the Mediterranean Sea and some of them, mostly the ones recovered from deep metagenomes, do not recruit in any database probably indicating higher variability and endemicity in Mediterranean bathypelagic waters. Together these data provide the first detailed picture of genomic diversity, spatial and depth variations of viral communities within the Mediterranean Sea using metagenome assembly.

Novel Caudovirales associated with Marine Group I Thaumarchaeota assembled from metagenomes.

Marine Group I (MGI) Thaumarchaeota are some of the most abundant microorganisms in the deep ocean and responsible for much of the ammonia oxidation occurring in this environment. In this work, we present 35 sequences assembled from metagenomic samples of the first uncultivated Caudovirales viruses associated with Thaumarchaeota, which we designated marthavirus. Most of the sequences were obtained from cellular metagenomes confirming that they represent an important tool to study environmental viral communities due to cells retrieved while undergoing viral lysis. Metagenomic recruitment showed that this viral population is formed by very divergent entities with high intrapopulation homogeneity. However, metatranscriptomic analyses revealed the same differential expression profile with the capsid as major transcript, indicative of viruses during the lytic cycle. The cobalamine biosynthesis gene cobS, an auxiliary metabolic gene, was also highly expressed during the infection. These analyses expand our understanding of the global diversity of archaeal viruses.

Flexible genomic island conservation across freshwater and marine Methylophilaceae.

The evolutionary trajectory of Methylophilaceae includes habitat transitions from freshwater sediments to freshwater and marine pelagial that resulted in genome reduction (genome-streamlining) of the pelagic taxa. However, the extent of genetic similarities in the genomic structure and microdiversity of the two genome-streamlined pelagic lineages (freshwater "Ca. Methylopumilus" and the marine OM43 lineage) has so far never been compared. Here, we analyzed complete genomes of 91 "Ca. Methylopumilus" strains isolated from 14 lakes in Central Europe and 12 coastal marine OM43 strains. The two lineages showed a remarkable niche differentiation with clear species-specific differences in habitat preference and seasonal distribution. On the other hand, we observed a synteny preservation in their genomes by having similar locations and types of flexible genomic islands (fGIs). Three main fGIs were identified: a replacement fGI acting as phage defense, an additive fGI harboring metabolic and resistance-related functions, and a tycheposon containing nitrogen-, thiamine-, and heme-related functions. The fGIs differed in relative abundances in metagenomic datasets suggesting different levels of variability ranging from strain-specific to population-level adaptations. Moreover, variations in one gene seemed to be responsible for different growth at low substrate concentrations and a potential biogeographic separation within one species. Our study provides a first insight into genomic microdiversity of closely related taxa within the family Methylophilaceae and revealed remarkably similar dynamics involving mobile genetic elements and recombination between freshwater and marine family members.

Evolutionary pathways for deep-sea adaptation in marine planktonic Actinobacteriota.

The deep ocean, one of the largest ecosystems on earth, is dominated by microorganisms that are keystones in the regulation of biogeochemical cycles. However, the evolutionary pathways underlying the specific adaptations required (e.g., high pressure and low temperature) by this unique niche remain understudied. Here, we analyzed the first representatives belonging to the order Acidimicrobiales, a group of marine planktonic Actinobacteriota, that specifically inhabits the aphotic zone of the oceanic water column (>200 m). Compared with their epipelagic counterparts, deep-sea representatives showed the same evolution in genome architecture with higher GC content, longer intergenic spaces as well as higher nitrogen (N-ARSC) and lower carbon (C-ARSC) content in encoded amino acid residue side chains consistent with the higher nitrogen concentration and lower carbon concentration in deep waters compared to the photic zone. Metagenomic recruitment showed distribution patterns that allowed the description of different ecogenomic units within the three deep water-associated genera defined by our phylogenomic analyses (UBA3125, S20-B6 and UBA9410). The entire genus UBA3125 was found exclusively associated with oxygen minimum zones linked to the acquisition of genes involved in denitrification. Genomospecies of genus S20-B6 recruited in samples from both mesopelagic (200-1,000 m) and bathypelagic (1000-4,000 m) zones, including polar regions. Diversity in the genus UBA9410 was higher, with genomospecies widely distributed in temperate zones, others in polar regions, and the only genomospecies associated with abyssal zones (>4,000 m). At the functional level, groups beyond the epipelagic zone have a more complex transcriptional regulation including in their genomes a unique WhiB paralog. In addition, they showed higher metabolic potential for organic carbon and carbohydrate degradation as well as the ability to accumulate glycogen as a source of carbon and energy. This could compensate for energy metabolism in the absence of rhodopsins, which is only present in genomes associated with the photic zone. The abundance in deep samples of cytochrome P450 monooxygenases associated with the genomes of this order suggests an important role in remineralization of recalcitrant compounds throughout the water column.

Phosphate-related genomic islands as drivers of environmental adaptation in the streamlined marine alphaproteobacterial HIMB59.

IMPORTANCE: These results shed light on the evolutionary strategies of microbes with streamlined genomes to adapt and survive in the oligotrophic conditions that dominate the surface waters of the global ocean. At the individual level, these microbes have been subjected to evolutionary constraints that have led to a more efficient use of nutrients, removing non-essential genes named as "streamlining theory." However, at the population level, they conserve a highly diverse gene pool in flexible genomic islands resulting in polyclonal populations on the same genomic background as an evolutionary response to environmental pressures. Localization of these islands at equivalent positions in the genome facilitates horizontal transfer between clonal lineages. This high level of environmental genomic heterogeneity could explain their cosmopolitan distribution. In the case of the order HIMB59 within the class Alphaproteobacteria, two factors exert evolutionary pressure and determine this intraspecific diversity: phages and the concentration of P in the environment.

A novel and diverse group of Candidatus Patescibacteria from bathypelagic Lake Baikal revealed through long-read metagenomics.

BACKGROUND: Lake Baikal, the world's deepest freshwater lake, contains important numbers of Candidatus Patescibacteria (formerly CPR) in its deepest reaches. However, previously obtained CPR metagenome-assembled genomes recruited very poorly indicating the potential of other groups being present. Here, we have applied for the first time a long-read (PacBio CCS) metagenomic approach to analyze in depth the Ca. Patescibacteria living in the bathypelagic water column of Lake Baikal at 1600 m. RESULTS: The retrieval of nearly complete 16S rRNA genes before assembly has allowed us to detect the presence of a novel and a likely endemic group of Ca. Patescibacteria inhabiting bathypelagic Lake Baikal. This novel group seems to possess extremely high intra-clade diversity, precluding complete genomes' assembly. However, read binning and scaffolding indicate that these microbes are similar to other Ca. Patescibacteria (i.e. parasites or symbionts), although they seem to carry more anabolic pathways, likely reflecting the extremely oligotrophic habitat they inhabit. The novel bins have not been found anywhere, but one of the groups appears in small amounts in an oligotrophic and deep alpine Lake Thun. We propose this novel group be named Baikalibacteria. CONCLUSION: The recovery of 16S rRNA genes via long-read metagenomics plus the use of long-read binning to uncover highly diverse "hidden" groups of prokaryotes are key strategies to move forward in ecogenomic microbiology. The novel group possesses enormous intraclade diversity akin to what happens with Ca. Patescibacteria at the interclade level, which is remarkable in an environment that has changed little in the last 25 million years.

Flotillin-associated rhodopsin (FArhodopsin), a widespread paralog of proteorhodopsin in aquatic bacteria with streamlined genomes.

Microbial rhodopsins are found more than once in a single genome (paralogs) often have different functions. We screened a large dataset of open ocean single-amplified genomes (SAGs) for co-occurrences of multiple rhodopsin genes. Many such cases were found among Pelagibacterales (SAR11), HIMB59, and the Gammaproteobacteria Pseudothioglobus SAGs. These genomes always had a bona fide proteorhodopsin and a separate cluster of genes containing a second rhodopsin associated with a predicted flotillin coding gene and have thus been named flotillin-associated rhodopsins (FArhodopsins). Although they are members of the proteorhodopsin protein family, they form a separate clade within that family and are quite divergent from known proton-pumping proteorhodopsins. They contain either DTT, DTL, or DNI motifs in their key functional amino acids. FArhodopsins are mainly associated with the lower layers of the epipelagic zone. All marine FArhodopsins had the retinal binding lysine, but we found relatives in freshwater metagenomes lacking this key amino acid. AlphaFold predictions of marine FArhodopsins indicate that their retinal pocket might be very reduced or absent, hinting that they are retinal-less. Freshwater FArhodopsins were more diverse than marine ones, but we could not determine if there were other rhodopsins in the genome due to the lack of SAGs or isolates. Although the function of FArhodopsins could not be established, their conserved genomic context indicated involvement in the formation of membrane microdomains. The conservation of FArhodopsins in diverse and globally abundant microorganisms suggests that they may be important in the adaptation to the twilight zone of aquatic environments. IMPORTANCE Rhodopsins have been shown to play a key role in the ecology of aquatic microbes. Here, we describe a group of widespread rhodopsins in aquatic microbes associated with dim light conditions. Their characteristic genomic context found in both marine and freshwater environments indicates a novel potential involvement in membrane microstructure that could be important for the function of the coexisting proteorhodopsin proton pumps. The absence or reduction of the retinal binding pocket points to a drastically different physiological role.

Role of Relebactam in the Antibiotic Resistance Acquisition in Pseudomonas aeruginosa: In Vitro Study.

BACKGROUND: Pseudomonas aeruginosa shows resistance to several antibiotics and often develops such resistance during patient treatment. OBJECTIVE: Develop an in vitro model, using clinical isolates of P. aeruginosa, to compare the ability of the imipenem and imipenem/relebactam to generate resistant mutants to imipenem and to other antibiotics. Perform a genotypic analysis to detect how the selective pressure changes their genomes. METHODS: The antibiotics resistance was studied by microdilution assays and e-test, and the genotypic study was performed by NGS. RESULTS: The isolates acquired resistance to imipenem in an average of 6 days, and to imipenem/relebactam in 12 days (p value = 0.004). After 30 days of exposure, 75% of the isolates reached a MIC > 64 mg/L for imipenem and 37.5% for imipenem/relebactam (p value = 0.077). The 37.5% and the 12.5% imipenem/relebactam mutants developed resistance to piperacillin/tazobactam and ceftazidime, respectively, while the 87.5% and 37.5% of the imipenem mutants showed resistance to these drugs (p value = 0.003, p value = 0.015). The main biological processes altered by the SNPs were the glycosylation pathway, transcriptional regulation, histidine kinase response, porins, and efflux pumps. DISCUSSION: The addition of relebactam delays the generation of resistance to imipenem and limits the cross-resistance to other beta-lactams. The clinical relevance of this phenomenon, which has the limitation that it has been performed in vitro, should be evaluated by stewardship programs in clinical practice, as it could be useful in controlling multi-drug resistance in P. aeruginosa.

A moderately thermophilic origin of a novel family of marine group II euryarchaeota from deep ocean.

Marine group II (MGII) is the most abundant planktonic heterotrophic archaea in the ocean. The evolutionary history of MGII archaea is elusive. In this study, 13 new MGII metagenome-assembled genomes were recovered from surface to the hadal zone in Challenger Deep of the Mariana Trench; four of them from the deep ocean represent a novel group. The optimal growth temperature (OGT) of the common ancestor of MGII has been estimated to be at about 60°C and OGTs of MGIIc, MGIIb, and MGIIa at 47°C-50ºC, 37°C-44ºC, and 30°C-37ºC, respectively, suggesting the adaptation of these species to different temperatures during evolution. The estimated OGT range of MGIIc was supported by experimental measurements of cloned ß-galactosidase that showed optimal enzyme activity around 50°C. These results indicate that MGIIc may have originated from a common ancestor that lived in warm or even hot marine environment, such as hydrothermal vents.

Long-Read Metagenomics Improves the Recovery of Viral Diversity from Complex Natural Marine Samples.

The recovery of DNA from viromes is a major obstacle in the use of long-read sequencing to study their genomes. For this reason, the use of cellular metagenomes (>0.2-µm size range) emerges as an interesting complementary tool, since they contain large amounts of naturally amplified viral genomes from prelytic replication. We have applied second-generation (Illumina NextSeq; short reads) and third-generation (PacBio Sequel II; long reads) sequencing to compare the diversity and features of the viral community in a marine sample obtained from offshore waters of the western Mediterranean. We found that a major wedge of the expected marine viral diversity was directly recovered by the raw PacBio circular consensus sequencing (CCS) reads. More than 30,000 sequences were detected only in this data set, with no homologues in the long- and short-read assembly, and ca. 26,000 had no homologues in the large data set of the Global Ocean Virome 2 (GOV2), highlighting the information gap created by the assembly bias. At the level of complete viral genomes, the performance was similar in both approaches. However, the hybrid long- and short-read assembly provided the longest average length of the sequences and improved the host assignment. Although no novel major clades of viruses were found, there was an increase in the intraclade genomic diversity recovered by long reads that produced an enriched assessment of the real diversity and allowed the discovery of novel genes with biotechnological potential (e.g., endolysin genes). IMPORTANCE We explored the vast genetic diversity of environmental viruses by using a combination of cellular metagenome (as opposed to virome) sequencing using high-fidelity long-read sequences (in this case, PacBio CCS). This approach resulted in the recovery of a representative sample of the viral population, and it performed better (more phage contigs, larger average contig size) than Illumina sequencing applied to the same sample. By this approach, the many biases of assembly are avoided, as the CCS reads recovers (typically around 5 kb) complete genes and even operons, resulting in a better discovery of the viral gene diversity based on viral marker proteins. Thus, biotechnologically promising genes, such as endolysin genes, can be very efficiently searched with this approach. In addition, hybrid assembly produces more complete and longer contigs, which is particularly important for studying little-known viral groups such as the nucleocytoplasmic large DNA viruses (NCLDV).

Enhanced Recovery of Microbial Genes and Genomes From a Marine Water Column Using Long-Read Metagenomics.

Third-generation sequencing has penetrated little in metagenomics due to the high error rate and dependence for assembly on short-read designed bioinformatics. However, second-generation sequencing metagenomics (mostly Illumina) suffers from limitations, particularly in the assembly of microbes with high microdiversity and retrieval of the flexible (adaptive) fraction of prokaryotic genomes. Here, we have used a third-generation technique to study the metagenome of a well-known marine sample from the mixed epipelagic water column of the winter Mediterranean. We have compared PacBio Sequel II with the classical approach using Illumina Nextseq short reads followed by assembly to study the metagenome. Long reads allow for efficient direct retrieval of complete genes avoiding the bias of the assembly step. Besides, the application of long reads on metagenomic assembly allows for the reconstruction of much more complete metagenome-assembled genomes (MAGs), particularly from microbes with high microdiversity such as Pelagibacterales. The flexible genome of reconstructed MAGs was much more complete containing many adaptive genes (some with biotechnological potential). PacBio Sequel II CCS appears particularly suitable for cellular metagenomics due to its low error rate. For most applications of metagenomics, from community structure analysis to ecosystem functioning, long reads should be applied whenever possible. Specifically, for in silico screening of biotechnologically useful genes, or population genomics, long-read metagenomics appears presently as a very fruitful approach and can be analyzed from raw reads before a computationally demanding (and potentially artifactual) assembly step.

Phylogenomics of SAR116 Clade Reveals Two Subclades with Different Evolutionary Trajectories and an Important Role in the Ocean Sulfur Cycle.

The SAR116 clade within the class Alphaproteobacteria represents one of the most abundant groups of heterotrophic bacteria inhabiting the surface of the ocean. The small number of cultured representatives of SAR116 (only two to date) is a major bottleneck that has prevented an in-depth study at the genomic level to understand the relationship between genome diversity and its role in the marine environment. In this study, we use all publicly available genomes to provide a genomic overview of the phylogeny, metabolism, and biogeography within the SAR116 clade. This increased genomic diversity has led to the discovery of two subclades that, despite coexisting in the same environment, display different properties in their genomic makeup. One represents a novel subclade for which no pure cultures have been isolated and is composed mainly of single-amplified genomes (SAGs). Genomes within this subclade showed convergent evolutionary trajectories with more streamlined features, such as low GC content (ca. 30%), short intergenic spacers (<22 bp), and strong purifying selection (low ratio of nonsynonymous to synonymous polymorphisms [dN/dS]). Besides, they were more abundant in metagenomic databases recruiting at the deep chlorophyll maximum. Less abundant and restricted to the upper photic layers of the global ocean, the other subclade of SAR116, enriched in metagenome-assembled genomes (MAGs), included the only two pure cultures. Genomic analysis suggested that both clades have a significant role in the sulfur cycle with differences in the way both clades can metabolize dimethylsulfoniopropionate (DMSP). IMPORTANCE The SAR116 clade of Alphaproteobacteria is a ubiquitous group of heterotrophic bacteria inhabiting the surface of the ocean, but the information about their ecology and population genomic diversity is scarce due to the difficulty of getting pure culture isolates. The combination of single-cell genomics and metagenomics has become an alternative approach to study these kinds of microbes. Our results expand the understanding of the genomic diversity, distribution, and lifestyles within this clade and provide evidence of different evolutionary trajectories in the genomic makeup of the two subclades that could serve to illustrate how evolutionary pressure can drive different adaptations to the same environment. Therefore, the SAR116 clade represents an ideal model organism for the study of the evolutionary streamlining of genomes in microbes that have relatively close relatedness to each other.

Genomic Characterization of Imipenem- and Imipenem-Relebactam-Resistant Clinical Isolates of Pseudomonas aeruginosa.

Pseudomonas aeruginosa is an opportunistic human pathogen and a major cause of nosocomial infections. The global spread of carbapenem-resistant strains is growing rapidly and has become a major public health challenge. Imipenem-relebactam (I/R) is a novel carbapenem-beta-lactamase inhibitor combination that can overcome carbapenem resistance. In this study, we aimed to understand the mechanism underlying resistance to imipenem and imipenem-relebactam. For this purpose, we performed a genomic comparison of 40 new clinical P. aeruginosa strains with different antibiotic sensitivity patterns as well as the presence/absence of carbapenemases. Results indicated the presence of a reduced flexible genome (15% total) mostly represented by phages and defense mechanisms against them, showing an important role in evolution and pathogenicity. We found a high diversity of antibiotic resistance genes grouped in small clusters mobilized via integrative and conjugative elements and facilitated by the high homologous recombination detected. Ortholog genes were found in several pathogenic strains from distantly related taxa in different mobile elements with a global distribution. The microdiversity found in those strains without carbapenemases did not reveal a clear pattern that could be associated with carbapenem resistance, suggesting multiple mechanisms of resistance in the core genome. Our results provide new insight into the dynamics and high genomic plasticity by which clinical strains of P. aeruginosa acquire resistance. This knowledge can be applied to other multidrug-resistant microbes to create predictive frameworks for assessing common molecular mechanisms of antibiotic resistance and integrated into new strategies for their prevention. IMPORTANCE The growing emergence and spread of carbapenem-resistant pathogens worldwide exacerbate the clinical challenge of treating these infections. Given the importance of carbapenems for the treatment of infections caused by Pseudomonas aeruginosa, this study aimed to investigate the underlying genomic properties of the clinical isolates that exhibited resistance to imipenem and imipenem-relebactam. This information will enhance our ability to forecast traits of resistant strains and design reliable treatments against this important threat. Our results provide new insight into the dynamics and high genomic plasticity by which clinical strains of P. aeruginosa acquire resistance as well as offers a methodology that can be applied to many other opportunistic pathogens with broad antibiotic resistance.

The microbiome of the Black Sea water column analyzed by shotgun and genome centric metagenomics.

BACKGROUND: The Black Sea is the largest brackish water body in the world, although it is connected to the Mediterranean Sea and presents an upper water layer similar to some regions of the former, albeit with lower salinity and temperature. Despite its well-known hydrology and physicochemical features, this enormous water mass remains poorly studied at the microbial genomics level. RESULTS: We have sampled its different water masses and analyzed the microbiome by shotgun and genome-resolved metagenomics, generating a large number of metagenome-assembled genomes (MAGs) from them. We found various similarities with previously described Black Sea metagenomic datasets, that show remarkable stability in its microbiome. Our datasets are also comparable to other marine anoxic water columns like the Cariaco Basin. The oxic zone resembles to standard marine (e.g. Mediterranean) photic zones, with Cyanobacteria (Synechococcus but a conspicuously absent Prochlorococcus), and photoheterotrophs domination (largely again with marine relatives). The chemocline presents very different characteristics from the oxic surface with many examples of chemolithotrophic metabolism (Thioglobus) and facultatively anaerobic microbes. The euxinic anaerobic zone presents, as expected, features in common with the bottom of meromictic lakes with a massive dominance of sulfate reduction as energy-generating metabolism, a few (but detectable) methanogenesis marker genes, and a large number of "dark matter" streamlined genomes of largely unpredictable ecology. CONCLUSIONS: The Black Sea oxic zone presents many similarities to the global ocean while the redoxcline and euxinic water masses have similarities to other similar aquatic environments of marine (Cariaco Basin or other Black Sea regions) or freshwater (meromictic monimolimnion strata) origin. The MAG collection represents very well the different types of metabolisms expected in this kind of environment. We are adding critical information about this unique and important ecosystem and its microbiome.

Nasopharyngeal Microbial Communities of Patients Infected With SARS-CoV-2 That Developed COVID-19.

BACKGROUND: SARS-CoV-2 is an RNA virus causing COVID-19. The clinical characteristics and epidemiology of COVID-19 have been extensively investigated, however, only one study so far focused on the patient's nasopharynx microbiota. In this study we investigated the nasopharynx microbial community of patients that developed different severity levels of COVID-19. We performed 16S ribosomal DNA sequencing from nasopharyngeal swab samples obtained from SARS-CoV-2 positive (56) and negative (18) patients in the province of Alicante (Spain) in their first visit to the hospital. Positive SARS-CoV-2 patients were observed and later categorized in mild (symptomatic without hospitalization), moderate (hospitalization), and severe (admission to ICU). We compared the microbiota diversity and OTU composition among severity groups and built bacterial co-abundance networks for each group. RESULTS: Statistical analysis indicated differences in the nasopharyngeal microbiome of COVID19 patients. 62 OTUs were found exclusively in SARS-CoV-2 positive patients, mostly classified as members of the phylum Bacteroidota (18) and Firmicutes (25). OTUs classified as Prevotella were found to be significantly more abundant in patients that developed more severe COVID-19. Furthermore, co-abundance analysis indicated a loss of network complexity among samples from patients that later developed more severe symptoms. CONCLUSION: Our study shows that the nasopharyngeal microbiome of COVID-19 patients showed differences in the composition of specific OTUs and complexity of co-abundance networks. Taxa with differential abundances among groups could serve as biomarkers for COVID-19 severity. Nevertheless, further studies with larger sample sizes should be conducted to validate these results.

Genomes of the "Candidatus Actinomarinales" Order: Highly Streamlined Marine Epipelagic Actinobacteria.

"Candidatus Actinomarinales" was defined as a subclass of exclusively marine Actinobacteria with small cells and genomes. We have collected all the available genomes in databases to assess the diversity included in this group and analyzed it by comparative genomics. We have found the equivalent of five genera and 18 genomospecies. They have genome reduction parameters equal to those of freshwater actinobacterial "Candidatus Nanopelagicales" or marine alphaproteobacterial Pelagibacterales Genome recruitment shows that they are found only in the photic zone and mainly in surface waters, with only one genus that is found preferentially at or below the deep chlorophyll maximum. "Ca Actinomarinales" show a highly conserved core genome (80% of the gene families conserved for the whole order) with a saturation of genomic diversity of the flexible genome at the genomospecies level. We found only a flexible genomic island preserved throughout the order; it is related to the sugar decoration of the envelope and uses several tRNAs as hot spots to increase its genomic diversity. Populations had a discrete level of sequence diversity similar to other marine microbes but drastically different from the much higher levels found for Pelagibacterales Genomic analysis suggests that they are all aerobic photoheterotrophs with one type 1 rhodopsin and a heliorhodopsin. Like other actinobacteria, they possess the F420 coenzyme biosynthesis pathway, and its lower reduction potential could provide access to an increased range of redox chemical transformations. Last, sequence analysis revealed the first "Ca Actinomarinales" phages, including a prophage, with metaviromic islands related to sialic acid cleavage.IMPORTANCE Microbiology is in a new age in which sequence databases are primary sources of information about many microbes. However, in-depth analysis of environmental genomes thus retrieved is essential to substantiate the new knowledge. Here, we study 182 genomes belonging to the only known exclusively marine pelagic group of the phylum Actinobacteria The aquatic branch of this phylum is largely known from environmental sequencing studies (single-amplified genomes [SAGs] and metagenome-assembled genomes [MAGs]), and we have collected and analyzed the available information present in databases about the "Ca. Actinomarinales." They are among the most streamlined microbes to live in the epipelagic zone of the ocean, and their study is critical to obtain a proper view of the diversity of Actinobacteria and their role in aquatic ecosystems.

The Evolutionary Success of the Marine Bacterium SAR11 Analyzed through a Metagenomic Perspective.

The SAR11 clade of Alphaproteobacteria is the most abundant group of planktonic cells in the near-surface epipelagic waters of the ocean, but the mechanisms underlying its exceptional success have not been fully elucidated. Here, we applied a metagenomic approach to explore microdiversity patterns by measuring the accumulation of synonymous and nonsynonymous mutations as well as homologous recombination in populations of SAR11 from different aquatic habitats (marine epipelagic, bathypelagic, and surface freshwater). The patterns of mutation accumulation and recombination were compared to those of other groups of representative marine microbes with multiple ecological strategies that share the same marine habitat, namely, Cyanobacteria (Prochlorococcus and Synechococcus), Archaea ("Candidatus Nitrosopelagicus" and Marine Group II Thalassoarchaea), and some heterotrophic marine bacteria (Alteromonas and Erythrobacter). SAR11 populations showed widespread recombination among distantly related members, preventing divergence leading to a genetically stable population. Moreover, their high intrapopulation sequence diversity with an enrichment in synonymous replacements supports the idea of a very ancient divergence and the coexistence of multiple different clones. However, other microbes analyzed seem to follow different evolutionary dynamics where processes of diversification driven by geographic and ecological instability produce a higher number of nonsynonymous replacements and lower intrapopulation sequence diversity. Together, these data shed light on some of the evolutionary and ecological processes that lead to the large genomic diversity in SAR11. Furthermore, this approach can be applied to other similar microbes that are difficult to culture in the laboratory, but abundant in nature, to investigate the underlying dynamics of their genomic evolution.IMPORTANCE As the most abundant bacteria in oceans, the Pelagibacterales order (here SAR11) plays an important role in the global carbon cycle, but the study of the evolutionary forces driving its evolution has lagged considerably due to the inherent difficulty of obtaining pure cultures. Multiple evolutionary models have been proposed to explain the diversification of distinct lineages within a population; however, the identification of many of these patterns in natural populations remains mostly enigmatic. We have used a metagenomic approach to explore microdiversity patterns in their natural habitats. Comparison with a collection of bacterial and archaeal groups from the same environments shows that SAR11 populations have a different evolutionary regime, where multiple genotypes coexist within the same population and remain stable over time. Widespread homologous recombination could be one of the main driving factors of this homogenization.

Dysbiosis in marine aquaculture revealed through microbiome analysis: reverse ecology for environmental sustainability.

The increasing demand for products for human consumption is leading to the fast-growing expansion of numerous food sectors such as marine aquaculture (mariculture). However, excessive input of nutrients and pollutants modifies marine ecosystems. Here, we applied a metagenomic approach to investigate these perturbations in samples from marine farms of gilthead seabream cultures. Results revealed dysbiosis and functional imbalance within the net cage with a unique structure, with little interference with samples from the fish microbiota or those collected far away from the coast. Remarkably, below the cage the prokaryotic community was highly similar to the marine microbiome of photic offshore samples. We recovered 48 novel metagenome-assembled genomes. Metagenomic recruitment revealed a significant change in the microbial community which was dominated by several Proteobacteria orders (Sphingomonadales, Pseudomonadales, Caudobacterales and Rhizobiales). Genomic potential for bioremediation processes, including nitrate removal through aerobic denitrification, and degradation of aromatic compounds and other toxic products were enriched in these microbes. The detrimental side effects were the increased number of antimicrobial resistance genes and the presence of potentially emergent pathogens. Knowledge of this metabolic diversity and the microbes involved in ecological balance recovery can be used to reduce the environmental impact of these practices.