Papers by Francisco Rodriguez-Valera
MSystems, Dec 22, 2020
Candidatus Actinomarinales" was defined as a subclass of exclusively marine Actinobacteria with s... more 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.
Proceedings of the National Academy of Sciences of the United States of America, Jul 16, 2019
Data deposition: DNA sequencing data have been deposited in the GenBank database with the identif... more Data deposition: DNA sequencing data have been deposited in the GenBank database with the identifiers MK570053 to MK570059. Accession number of the Nitrosopumilus strain SW genome is CP035425.
Proceedings of the National Academy of Sciences of the United States of America, Feb 7, 2020
bioRxiv (Cold Spring Harbor Laboratory), Aug 15, 2019
Bacteriophages are the most numerous biological entities on earth and play a crucial role in shap... more Bacteriophages are the most numerous biological entities on earth and play a crucial role in shaping microbial communities. Investigating the bacteriophage community from soil samples will shed light not only on the yet largely unknown phage diversity, but also may result in novel insights into phage biology and functioning. Unfortunately, the study of soil viromes lags far behind any other ecological model system, due to the heterogeneous soil matrix that rises major technical difficulties in the extraction process. Resolving these technical challenges and establishing a standardized extraction protocol is therefore a fundamental prerequisite for replicable results and comparative virome studies. We here report the optimization of protocols for extraction of bacteriophage DNA from soil preceding metagenomic analysis such that the protocol can equally be harnessed for phage isolation. As an optimization strategy, soil samples were spiked with a viral community consisting of phages from different families (10 6 PFU/g soil): Listeria phage ΦA511 (Myovirus), Staphylococcus phage Φ2638AΔLCR (Siphovirus), and Escherichia phage ΦT7 (Podovirus). The efficacy of bacteriophage (i) elution, (ii) filtration, (iii) concentration, and (iv) DNA extraction methods was tested. Successful extraction routes were selected based on spiked phage recovery and low bacterial 16S rRNA gene contaminants. Natural agricultural soil viromes were then extracted with the optimized methods and shotgun sequenced. Our approach yielded sufficient amounts of inhibitor-free viral DNA for non-amplification dependent sequencing and low 16S rRNA gene contamination levels (≤ 0.2 ‰). Compared to previously published protocols, the number of bacterial read contamination was decreased by 65 %. In addition, 468 novel circularized soil phage genomes in size up to 235 kb were obtained from over 29,000 manually identified viral contigs, promising the discovery of a large, previously inaccessible viral diversity. .
Searching Metagenomes for New Rhodopsins
Springer eBooks, 2022
bioRxiv (Cold Spring Harbor Laboratory), Sep 12, 2019
Rhodopsins are the most abundant light-harvesting proteins. A new family of rhodopsins, heliorhod... more Rhodopsins are the most abundant light-harvesting proteins. A new family of rhodopsins, heliorhodopsins (HeRs), was recently discovered. In opposite to the known rhodopsins their N-termini face the cytoplasm. HeRs structure and function remain unknown. We present structures of two HeR-48C12 states at 1.5 Å showing its remarkable difference from all known rhodopsins. Its internal extracellular part is completely hydrophobic, while the cytoplasmic part comprises a cavity ('active site'), surrounded by charged amino acids and containing a cluster of water molecules, presumably being a primary proton acceptor from the Schiff base. At acidic pH a planar triangle molecule (acetate) is present in the 'active site' which demonstrated its ability to maintain such anions as carbonate or nitrate. Structure-based bioinformatic analysis identified 10 subfamilies of HeRs suggesting their diverse biological functions. The structures and available data suggest an enzymatic activity of HeR-48C12 subfamily and their possible involvement into fundamental redox biological processes. .
Research Square (Research Square), Nov 5, 2020
Background: The Black Sea is the largest brackish water body in the world, although it is connect... more 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 (mostly) temperature. In spite of its well-known hydrology and physico chemistry, 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 classic and genome-resolved metagenomics generating a large number of metagenome-assembled genomes (MAGs) from them. The oxic zone presents many similarities to the global ocean while the euxinic water mass has similarities to other similar aquatic environments of marine or freshwater (meromictic monimolimnion strata) origin. The MAG collection represents very well the different types of metabolisms expected in this kind of environments and includes Cyanobacteria (Synechococcus), photoheterotrophs (largely with marine relatives), facultative/microaerophilic microbes again largely marine, chemolithotrophs (N and S oxidizers) and a large number of anaerobes, mostly sulfate reducers but also a few methanogens and a large number of "dark matter" streamlined genomes of largely unpredictable ecology. Conclusions: The Black Sea presents a mixture of similarities to other water bodies. The photic zone has many microbes in common with that of the Mediterranean with the relevant exception of the absence of Prochlorococcus. The chemocline already presents very different characteristics with many examples of chemolithotrophic metabolism (Thioglobus) and facultatively anaerobic microbes. Finally 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 and a small but detectable methanogenesis.We are adding critical information about this unique and important ecosystem and its microbiome.
Microbial Rhodopsins
Springer eBooks, 2022
Long-Read Metagenomics Improves the Recovery of Viral Diversity from Complex Natural Marine Samples
MSystems, Jun 28, 2022
We explored the vast genetic diversity of environmental viruses by using a combination of cellula... more 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.
Systematic and Applied Microbiology, Jul 1, 1986
After an extensive taxonomic study for extremely halophilic rods , several groups appeared which ... more After an extensive taxonomic study for extremely halophilic rods , several groups appeared which did not fit with any previously described species of Halobacteria. In the present work two strains, Y27 (= ATCC 33960) and MA2.38 (= ATCC 33959), representatives of the rwo most distinct new groups, were studied in detail, and two new species are proposed.
Microbiome, Feb 20, 2019
Background: Polynyas in the Southern Ocean are regions of intense primary production, mainly by P... more Background: Polynyas in the Southern Ocean are regions of intense primary production, mainly by Phaeocystis antarctica. Carbon fixed by phytoplankton in the water column is transferred to higher trophic levels, and finally, to the deep ocean. However, in the Amundsen Sea, most of this organic carbon does not reach the sediment but is degraded in the water column due to high bacterial heterotrophic activity. Results: We reconstructed 12 key bacterial genomes from different phases of bloom and analyzed the expression of genes involved in organic carbon remineralization. A high correlation of gene expression between the peak and decline phases was observed in an individual genome bin-based pairwise comparison of gene expression. Polaribacter belonging to Bacteroidetes was found to be dominant in the peak phase, and its transcriptional activity was high (48.9% of the total mRNA reads). Two dominant Polaribacter bins had the potential to utilize major polymers in P. antarctica, chrysolaminarin and xylan, with a distinct set of glycosyl hydrolases. In the decline phase, Gammaproteobacteria (Ant4D3, SUP05, and SAR92), with the potential to utilize low molecular weight-dissolved organic matter (LMW-DOM) including compatible solutes, was increased. The versatility of Gammaproteobacteria may contribute to their abundance in organic carbon-rich polynya waters, while the SAR11 clade was found to be predominant in the sea ice-covered oligotrophic ocean. SAR92 clade showed transcriptional activity for utilization of both polysaccharides and LMW-DOM; this may account for their abundance both in the peak and decline phases. Ant4D3 clade was dominant in all phases of the polynya bloom, implicating the crucial roles of this clade in LMW-DOM remineralization in the Antarctic polynyas. Conclusions: Genomic reconstruction and in situ gene expression analyses revealed the unique metabolic potential of dominant bacteria of the Antarctic polynya at a finer taxonomic level. The information can be used to predict temporal community succession linked to the availability of substrates derived from the P. antarctica bloom. Global warming has resulted in compositional changes in phytoplankton from P. antarctica to diatoms, and thus, repeated parallel studies in various polynyas are required to predict global warming-related changes in carbon remineralization.
MSystems, Feb 11, 2020
The SAR11 clade is one of the most abundant bacterioplankton groups in surface waters of most of ... more The SAR11 clade is one of the most abundant bacterioplankton groups in surface waters of most of the oceans and lakes. However, only 15 SAR11 phages have been isolated thus far, and only one of them belongs to the Myoviridae family (pelagimyophages). Here, we have analyzed 26 sequences of myophages that putatively infect the SAR11 clade. They have been retrieved by mining ca. 45 Gbp aquatic assembled cellular metagenomes and viromes. Most of the myophages were obtained from the cellular fraction (0.2 m), indicating a bias against this type of virus in viromes. We have found the first myophages that putatively infect Candidatus Fonsibacter (freshwater SAR11) and another group putatively infecting bathypelagic SAR11 phylogroup Ic. The genomes have similar sizes and maintain overall synteny in spite of low average nucleotide identity values, revealing high similarity to marine cyanomyophages. Pelagimyophages recruited metagenomic reads widely from several locations but always much more from cellular metagenomes than from viromes, opposite to what happens with pelagipodophages. Comparing the genomes resulted in the identification of a hypervariable island that is related to host recognition. Interestingly, some genes in these islands could be related to host cell wall synthesis and coinfection avoidance. A cluster of curli-related proteins was widespread among the genomes, although its function is unclear. IMPORTANCE SAR11 clade members are among the most abundant bacteria on Earth. Their study is complicated by their great diversity and difficulties in being grown and manipulated in the laboratory. On the other hand, and due to their extraordinary abundance, metagenomic data sets provide enormous richness of information about these microbes. Given the major role played by phages in the lifestyle and evolution of prokaryotic cells, the contribution of several new bacteriophage genomes preying on this clade opens windows into the infection strategies and life cycle of its viruses. Such strategies could provide models of attack of large-genome phages preying on streamlined aquatic microbes.
Nature Communications, Jan 17, 2019
Smacoviridae is a family of small (~2.5 Kb) CRESS-DNA (Circular Rep Encoding Single-Stranded (ss)... more Smacoviridae is a family of small (~2.5 Kb) CRESS-DNA (Circular Rep Encoding Single-Stranded (ss) DNA) viruses. These viruses have been found in faeces, were thought to infect eukaryotes and are suspected to cause gastrointestinal disease in humans. CRISPR-Cas systems are adaptive immune systems in prokaryotes, wherein snippets of genomes from invaders are stored as spacers that are interspersed between a repeated CRISPR sequence. Here we report several spacer sequences in the faecal archaeon Candidatus Methanomassiliicoccus intestinalis matching smacoviruses, implicating the archaeon as a firm candidate for a host. This finding may be relevant to understanding the potential origin of smacovirus-associated human diseases. Our results support that CRESS-DNA viruses can infect non-eukaryotes, which would mean that smacoviruses are the viruses with the smallest genomes to infect prokaryotes known to date. A probable target strand bias suggests that, in addition to double-stranded DNA, the CRISPR-Cas system can target ssDNA.
Most marine metagenomic studies of the marine photic zone analyze only samples taken at one or tw... more Most marine metagenomic studies of the marine photic zone analyze only samples taken at one or two depths. However, when the water column is stratified, physicochemical parameters change dramatically over relatively short depth intervals. We sampled the photic water column every 15m depth at a single point of an off-shore Mediterranean site during a period of strong stratification (early autumn) to evaluate the effects of small depth increases on the microbiome. Using genomic assembly and metagenomic read recruitment, we found major shifts in the community structure over small variations of depth, with most microbes showing a distribution limited to layers approximately 30 meters thick (stenobathic). Only some representatives of the SAR11 clade and the Sphingomonadaceae appeared to be eurybathic, spanning a greater range of depths. These results were confirmed by studying a single gene (rhodopsin) for which we also found narrow depth distributions. Our results highlight the importance of considering vertical distribution as a major element when analyzing the presence of marine clades and species or comparing the microbiome present at different locations.
BMC Genomics, Jan 5, 2017
Background: To develop evolutionary models for the free living bacterium Alteromonas the genome s... more Background: To develop evolutionary models for the free living bacterium Alteromonas the genome sequences of isolates of the genus have been extensively analyzed. However, the main genetic exchange drivers in these microbes, conjugative elements (CEs), have not been considered in detail thus far. In this work, CEs have been searched in several complete Alteromonas genomes and their sequence studied to understand their role in the evolution of this genus. Six genomes are reported here for the first time. Results: We have found nine different plasmids of sizes ranging from 85 to 600 Kb, most of them were found in a single strain. Networks of gene similarity could be established among six of the plasmids that were also connected with another cluster of plasmids found in Shewanella strains. The cargo genes found in these plasmids included cassettes found before in chromosome flexible genomic islands of Alteromonas strains. We describe also the plasmids pAMCP48-600 and pAMCP49-600, the largest found in Alteromonas thus far (ca. 600 Kb) and containing all the hallmarks to be classified as chromids. We found in them some housekeeping genes and a cluster that code for an exocellular polysaccharide. They could represent the transport vectors for the previously described replacement flexible genomic islands. Integrative and conjugative elements (ICEs) were more common than plasmids and showed similar patterns of variation with cargo genes coding for components of additive flexible genomic islands. A nearly identical ICE was found in A. mediterranea MED64 and Vibrio cholera AHV1003 isolated from a human pathogen, indicating the potential exchange of these genes across phylogenetic distances exceeding the family threshold. We have seen evidence of how CEs can be vectors to transfer gene cassettes acquired in the chromosomal flexible genomic islands, both of the additive and replacement kind. These CEs showed evidence of how genetic material is exchanged among members of the same species but also (albeit less frequently) across genus and family barriers. These gradients of exchange frequency are probably one of the main drivers of species origin and maintenance in prokaryotes and also provide these taxa with large genetic diversity.
MSystems, Oct 27, 2020
The SAR11 clade of Alphaproteobacteria is the most abundant group of planktonic cells in the near... more 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.
Frontiers in Microbiology, Dec 19, 2019
We performed an incubation experiment of seawater confined in plastic bottles with samples collec... more We performed an incubation experiment of seawater confined in plastic bottles with samples collected at three depths (15, 60, and 90 m) after retrieval from a single offshore location in the Mediterranean Sea, from a late summer stratified water column. Two samples representative of each depth were collected and stored in opaque bottles after two periods of 7 h. We took advantage of the "bottle effect" to investigate changes in the natural microbial communities (abundant and rare). We recovered 94 metagenome-assembled genomes (MAGs) and 1089 metagenomic viral contigs and examined their abundance using metagenomic recruitment. We detected a significant fast growth of copiotrophic bacteria such as Alteromonas or Erythrobacter throughout the entire water column with different dynamics that we assign to "clonal," "polyclonal," or "multispecies" depending on the recruitment pattern. Results also showed a marked ecotype succession in the phototropic picocyanobacteria that were able to grow at all the depths in the absence of light, highlighting the importance of their mixotrophic potential. In addition, "wall-chain-reaction" hypothesis based on the study of phagehost dynamics showed the higher impact of viral predation on archaea in deeper waters, evidencing their prominent role during incubations. Our results provide a step forward in understanding the mechanisms underlying dynamic patterns and ecology of the marine microbiome and the importance of processing the samples immediately after collection to avoid changes in the community structure.
The ISME Journal, Jan 15, 2018
Interactions between co-existing microorganisms deeply affect the physiology of the involved orga... more Interactions between co-existing microorganisms deeply affect the physiology of the involved organisms and, ultimately, the function of the ecosystem as a whole. Copiotrophic Alteromonas are marine gammaproteobacteria that thrive during the late stages of phytoplankton blooms in the marine environment and in laboratory co-cultures with cyanobacteria such as Trichodesmium. The response of this heterotroph to the sometimes rapid and transient changes in nutrient supply when the phototroph crashes is not well understood. Here, we isolated and sequenced the strain Alteromonas macleodii str. Te101 from a laboratory culture of Trichodesmium erythraeum IMS101, yielding a chromosome of 4.63 Mb and a single plasmid of 237 kb. Increasing salinities to ≥43 ppt inhibited the growth of Trichodesmium but stimulated growth of the associated Alteromonas. We characterized the transcriptomic responses of both microorganisms and identified the complement of active transcriptional start sites in Alteromonas at single-nucleotide resolution. In replicate cultures, a similar set of genes became activated in Alteromonas when growth rates of Trichodesmium declined and mortality was high. The parallel activation of fliA, rpoS and of flagellar assembly and growth-related genes indicated that Alteromonas might have increased cell motility, growth, and multiple biosynthetic activities. Genes with the highest expression in the data set were three small RNAs (Aln1a-c) that were identified as analogs of the small RNAs CsrB-C in E. coli or RsmX-Z in pathogenic bacteria. Together with the carbon storage protein A (CsrA) homolog Te101_05290, these RNAs likely control the expression of numerous genes in responding to changes in the environment. Shengwei Hou and Mario López-Pérez are co-sharing first authors.
Applied and Environmental Microbiology, 1986
Some species of extremely halophilic archaebacteria, Halobacteriaceae, have been shown to accumul... more Some species of extremely halophilic archaebacteria, Halobacteriaceae, have been shown to accumulate large amounts of poly(P-hydroxybutyrate) under conditions of nitrogen limitation and abundant carbon source. The production of poly(P3-hydroxybutyrate), at least in large quantities, was restricted to two carbohydrate-utilizing species, Halobacterium mediterranei and H. volcanii. In addition to the nutrients in the media, the salt concentration also influenced poly(P-hydroxybutyrate) accumulation, which was greater at lower salt concentrations. The possible application of these microorganisms for the production of biodegradable plastics is discussed.
Microbiome, Feb 11, 2020
Background: Bacteriophages (phages) are the most numerous biological entities on Earth and play a... more Background: Bacteriophages (phages) are the most numerous biological entities on Earth and play a crucial role in shaping microbial communities. Investigating the bacteriophage community from soil will shed light not only on the yet largely unknown phage diversity, but may also result in novel insights towards their functioning in the global biogeochemical nutrient cycle and their significance in earthbound ecosystems. Unfortunately, information about soil viromes is rather scarce compared to aquatic environments, due to the heterogeneous soil matrix, which rises major technical difficulties in the extraction process. Resolving these technical challenges and establishing a standardized extraction protocol is, therefore, a fundamental prerequisite for replicable results and comparative virome studies. Results: We here report the optimization of protocols for the extraction of phage DNA from agricultural soil preceding metagenomic analysis such that the protocol can equally be harnessed for phage isolation. As an optimization strategy, soil samples were spiked with Listeria phage A511 (Myovirus), Staphylococcus phage 2638AΔLCR (Siphovirus) and Escherichia phage T7 (Podovirus) (each 10 6 PFU/g soil). The efficacy of phage (i) elution, (ii) filtration, (iii) concentration and (iv) DNA extraction methods was tested. Successful extraction routes were selected based on spiked phage recovery and low bacterial 16S rRNA gene contaminants. Natural agricultural soil viromes were then extracted with the optimized methods and shotgun sequenced. Our approach yielded sufficient amounts of inhibitor-free viral DNA for shotgun sequencing devoid of amplification prior library preparation, and low 16S rRNA gene contamination levels (≤ 0.2‰). Compared to previously published protocols, the number of bacterial read contamination was decreased by 65%. In addition, 379 novel putative complete soil phage genomes (≤ 235 kb) were obtained from over 13,000 manually identified viral contigs, promising the discovery of a large, previously inaccessible viral diversity. We have shown a considerably enhanced extraction of the soil phage community by protocol optimization that has proven robust in both culture-dependent as well as through viromic analyses. Our huge data set of manually curated soil viral contigs substantially increases the amount of currently available soil virome data, and provides insights into the yet largely undescribed soil viral sequence space.
Uploads
Papers by Francisco Rodriguez-Valera