Uncovering novel bacterial and archaeal diversity: genomic insights from metagenome-assembled genomes in Cuatro Cienegas, Coahuila.

A comprehensive study was conducted in the Cuatro Ciénegas Basin (CCB) in Coahuila, Mexico, which is known for its diversity of microorganisms and unique physicochemical properties. The study focused on the "Archaean Domes" (AD) site in the CCB, which is characterized by an abundance of hypersaline, non-lithifying microbial mats. In AD, we analyzed the small domes and circular structures using metagenome assembly genomes (MAGs) with the aim of expanding our understanding of the prokaryotic tree of life by uncovering previously unreported lineages, as well as analyzing the diversity of bacteria and archaea in the CCB. A total of 325 MAGs were identified, including 48 Archaea and 277 Bacteria. Remarkably, 22 archaea and 104 bacteria could not be classified even at the genus level, highlighting the remarkable novel diversity of the CCB. Besides, AD site exhibited significant diversity at the phylum level, with Proteobacteria being the most abundant, followed by Desulfobacteria, Spirochaetes, Bacteroidetes, Nanoarchaeota, Halobacteriota, Cyanobacteria, Planctomycetota, Verrucomicrobiota, Actinomycetes and Chloroflexi. In Archaea, the monophyletic groups of MAGs belonged to the Archaeoglobi, Aenigmarchaeota, Candidate Nanoarchaeota, and Halobacteriota. Among Bacteria, monophyletic groups were also identified, including Spirochaetes, Proteobacteria, Planctomycetes, Actinobacteria, Verrucomicrobia, Bacteroidetes, Candidate Bipolaricaulota, Desulfobacteria, and Cyanobacteria. These monophyletic groups were possibly influenced by geographic isolation, as well as the extreme and fluctuating environmental conditions in the pond AD, such as stoichiometric imbalance of C:N:P of 122:42:1, fluctuating pH (5-9.8) and high salinity (5.28% to saturation).

Comprehensive comparative analysis of the periodontal pathogen Porphyromonas gingivalis: exploring the pan-genome, the reconstruction of the gene regulatory network and genome-scale metabolic network.

Porphyromonas gingivalis is a nonmotile, obligate anaerobic, Gram-negative bacterium known for its association with periodontal disease and its involvement in systemic diseases such as atherosclerosis, cardiovascular disease, colon cancer, and Alzheimer's disease. This bacterium produces several virulence factors, including capsules, fimbriae, lipopolysaccharides, proteolytic enzymes, and hemagglutinins. A comparative genomic analysis revealed the open pangenome of P. gingivalis and identified complete type IV secretion systems in strain KCOM2805 and almost complete type VI secretion systems in strains KCOM2798 and ATCC49417, which is a new discovery as previous studies did not find the proteins involved in secretion systems IV and VI. Conservation of some virulence factors between different strains was observed, regardless of their genetic diversity and origin. In addition, we performed for the first time a reconstruction analysis of the gene regulatory network, identifying transcription factors and proteins involved in the regulatory mechanisms of bacterial pathogenesis. In particular, QseB regulates the expression of hemagglutinin and arginine deaminase, while Rex may suppress the release of gingipain through interactions with PorV and the formatum/nitrate transporter. Our study highlights the central role of conserved virulence factors and regulatory pathways, particularly QseB and Rex, in P. gingivalis and provides insights into potential therapeutic targets.

Campylobacter prophage diversity reveals pervasive recombination between prophages from different Campylobacter species.

IMPORTANCE: Prophages play an important role in shaping the genetic diversity and evolution of their hosts. Acquisition or loss of prophages can lead to genomic variations, including changes in the bacterial phenotype promoted by recombination events, genetic repertoire exchanges and dissemination of virulence factors, and antibiotic resistance. By studying prophages in Campylobacter species, scientists can gain insights into the evolutionary patterns, pathogenicity mechanisms, epidemiology, and population dynamics of these species. This has implications for public health, antibiotic resistance surveillance, and the development of targeted therapeutic approaches.

Genetic insights into the microevolutionary dynamics and early introductions of human monkeypox virus in Mexico.

The recent global outbreak of mpox, caused by monkeypox virus (MPV) emerged in Europe in 2022 and rapidly spread to over 40 countries. The Americas are currently facing the highest impact, reporting over 50,000 cases by early 2023. In this study, we analyzed 880 MPV isolates worldwide to gain insights into the evolutionary patterns and initial introduction events of the virus in Mexico. We found that MPV entered Mexico on multiple occasions, from the United Kingdom, Portugal, and Canada, and subsequently spread locally in different regions of Mexico. Additionally, we show that MPV has an open pangenome, highlighting the role of gene turnover in shaping its genomic diversity, rather than single-nucleotide polymorphisms (SNPs), which do not contribute significantly to genome diversity. Although the genome contains multiple SNPs in coding regions, these remain under purifying selection, suggesting their evolutionary conservation. One notable exception is amino acid position 63 of the protein encoded by the Cop-A4L gene, which is intricately related to viral maturity, which we found to be under strong positive selection. Ancestral state reconstruction indicated that the ancestral state at position 63 corresponds to the amino acid valine, which is present only in isolates of clade I. However, the isolates from the current outbreak contained threonine at position 63. Our findings contribute new information about the evolution of monkeypox virus.

HADEG: A curated hydrocarbon aerobic degradation enzymes and genes database.

Databases of genes and enzymes involved in hydrocarbon degradation have been previously reported. However, these databases specialize on only a specific group of hydrocarbons and/or are constructed partly based on enzyme sequences with putative functions indicated by in silico research, with no experimental evidence. Here, we present a curated database of Hydrocarbon Aerobic Degradation Enzymes and Genes (HADEG) containing proteins and genes involved in alkane, alkene, aromatic, and plastic aerobic degradation and biosurfactant production based solely on experimental evidence, which are present in bacteria, and fungi. HADEG includes 259 proteins for petroleum hydrocarbon degradation, 160 for plastic degradation, and 32 for biosurfactant production. This database will help identify and predict hydrocarbon degradation genes/pathways and biosurfactant production in genomes.

Emergence of Raoultella ornithinolytica in human infections from different hospitals in Ecuador with OXA-48-producing resistance.

Purpose: The purpose of this study was to highlight the clinical and molecular features of 13 Raoultella ornithinolytica strains isolated from clinical environments in Ecuador, and to perform comparative genomics with previously published genomes of Raoultella spp. As Raoultella is primarily found in environmental, clinical settings, we focused our work on identifying mechanisms of resistance that can provide this bacterium an advantage to establish and persist in hospital environments. Methods: We analyzed 13 strains of Raoultella ornithinolytica isolated from patients with healthcare associated infections (HAI) in three hospitals in Quito and one in Santo Domingo de Los Tsáchilas, Ecuador, between November 2017 and April 2018. These isolates were subjected to phenotypic antimicrobial susceptibility testing, end-point polymerase chain reaction (PCR) to detect the presence of carbapenemases and whole-genome sequencing. Results: Polymerase chain reaction revealed that seven isolates were positive isolates for blaOXA-48 and one for blaKPC-2 gene. Of the seven strains that presented the blaOXA-48 gene, six harbored it on an IncFII plasmid, one was inserted into the bacterial chromosome. The blaKPC gene was detected in an IncM2/IncR plasmid. From the bioinformatics analysis, nine genomes had the gene blaOXA-48, originating from Ecuador. Moreover, all R. ornithinolytica strains contained the ORN-1 gene, which confers resistance for ß-lactams, such as penicillins and cephalosporins. Comparative genome analysis of the strains showed that the pangenome of R. ornithinolytica is considered an open pangenome, with 27.77% of core genes, which could be explained by the fact that the antibiotic resistance genes in the ancestral reconstruction are relatively new, suggesting that this genome is constantly incorporating new genes. Conclusion: These results reveal the genome plasticity of R. ornithinolytica, particularly in acquiring antibiotic-resistance genes. The genomic surveillance and infectious control of these uncommon species are important since they may contribute to the burden of antimicrobial resistance and human health.

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.

Identification of Evolutionary Trajectories Shared across Human Betacoronaviruses.

Comparing the evolution of distantly related viruses can provide insights into common adaptive processes related to shared ecological niches. Phylogenetic approaches, coupled with other molecular evolution tools, can help identify mutations informative on adaptation, although the structural contextualization of these to functional sites of proteins may help gain insight into their biological properties. Two zoonotic betacoronaviruses capable of sustained human-to-human transmission have caused pandemics in recent times (SARS-CoV-1 and SARS-CoV-2), although a third virus (MERS-CoV) is responsible for sporadic outbreaks linked to animal infections. Moreover, two other betacoronaviruses have circulated endemically in humans for decades (HKU1 and OC43). To search for evidence of adaptive convergence between established and emerging betacoronaviruses capable of sustained human-to-human transmission (HKU1, OC43, SARS-CoV-1, and SARS-CoV-2), we developed a methodological pipeline to classify shared nonsynonymous mutations as putatively denoting homoplasy (repeated mutations that do not share direct common ancestry) or stepwise evolution (sequential mutations leading towards a novel genotype). In parallel, we look for evidence of positive selection and draw upon protein structure data to identify potential biological implications. We find 30 candidate mutations, from which 4 (codon sites 18121 [nsp14/residue 28], 21623 [spike/21], 21635 [spike/25], and 23948 [spike/796]; SARS-CoV-2 genome numbering) further display evolution under positive selection and proximity to functional protein regions. Our findings shed light on potential mechanisms underlying betacoronavirus adaptation to the human host and pinpoint common mutational pathways that may occur during establishment of human endemicity.

Evidence of natural selection and dominance of SARS-CoV-2 variant Lambda (C.37) over variants of concern in Cusco, Peru.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineage C.37 (Lambda) has spread rapidly in Peru and other Latin American countries. However, most studies in Peru have focused on Lima, the capital city, without knowing the dynamics of the spread of the variant in other departments. Cusco, Peru, is one of the most popular departments in the country for tourists, so the introduction of new variants of SARS-CoV-2 might occur despite closure of the borders. Therefore, in this work, we analyzed the variants circulating in Cusco. The aim of this work was to better understand the distribution of SARS-CoV-2 lineages circulating in Cusco and to characterize the genomes of these strains. To this end, 46 SARS-CoV-2 genomes from vaccinated and unvaccinated patients were sequenced in the first half of 2021. The genomes were analyzed using phylogenetic and natural selection methods. Phylogenetic trees from Cusco showed dominance of the Lambda lineage over the variants of concern (VOCs), and there was no clustering of variants by district. Natural selection analysis revealed mutations, mainly in the spike protein, at positions 75, 246, 247, 707, 769, and 1020. In addition, we found that unvaccinated patients accumulated more new mutations than did vaccinated patients, and these included the F101Y mutation in ORF7a, E419A in NSP3, a deletion in S (21,618-22,501), and a deletion in ORF3a (25,437-26,122).

Identification of evolutionary trajectories shared across human betacoronaviruses.

Comparing the evolution of distantly related viruses can provide insights into common adaptive processes related to shared ecological niches. Phylogenetic approaches, coupled with other molecular evolution tools, can help identify mutations informative on adaptation, whilst the structural contextualization of these to functional sites of proteins may help gain insight into their biological properties. Two zoonotic betacoronaviruses capable of sustained human-to-human transmission have caused pandemics in recent times (SARS-CoV-1 and SARS-CoV-2), whilst a third virus (MERS-CoV) is responsible for sporadic outbreaks linked to animal infections. Moreover, two other betacoronaviruses have circulated endemically in humans for decades (HKU1 and OC43). To search for evidence of adaptive convergence between established and emerging betacoronaviruses capable of sustained human-to-human transmission (HKU1, OC43, SARS-CoV-1 and SARS-CoV-2), we developed a methodological pipeline to classify shared non-synonymous mutations as putatively denoting homoplasy (repeated mutations that do not share direct common ancestry) or stepwise evolution (sequential mutations leading towards a novel genotype). In parallel, we look for evidence of positive selection, and draw upon protein structure data to identify potential biological implications. We find 30 mutations, with four of these [codon sites 18121 (nsp14/residue 28), 21623 (spike/21), 21635 (spike/25) and 23948 (spike/796); SARS-CoV-2 genome numbering] displaying evolution under positive selection and proximity to functional protein regions. Our findings shed light on potential mechanisms underlying betacoronavirus adaptation to the human host and pinpoint common mutational pathways that may occur during establishment of human endemicity.

Genome Evolution and Early Introductions of the SARS-CoV-2 Omicron Variant in Mexico.

A new variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), named Omicron (Pango lineage designation B.1.1.529), was first reported to the World Health Organization by South African health authorities on 24 November 2021. The Omicron variant possesses numerous mutations associated with increased transmissibility and immune escape properties. In November 2021, Mexican authorities reported Omicron's presence in the country. In this study, we infer the first introductory events of Omicron and the impact that human mobility has had on the spread of the virus. We also evaluated the adaptive evolutionary processes in Mexican SARS-CoV-2 genomes during the first month of the circulation of Omicron. We inferred 160 introduction events of Omicron in Mexico since its first detection in South Africa; subsequently, after the first introductions there was an evident increase in the prevalence of SARS-CoV-2 during January. This higher prevalence of the novel variant resulted in a peak of reported cases; on average 6 weeks after, a higher mobility trend was reported. During the peak of cases in the country from January to February 2022, the Omicron BA.1.1 sub-lineage dominated, followed by the BA.1 and BA.15 sub-lineages. Additionally, we identified the presence of diversifying natural selection in the genomes of Omicron and found six non-synonymous mutations in the receptor binding domain of the spike protein, all of them related to evasion of the immune response. In contrast, the other proteins in the genome are highly conserved; however, we identified homoplasic mutations in non-structural proteins, indicating a parallel evolution.

Importance of Microbiome of Fecal Samples Obtained from Adolescents with Different Weight Conditions on Resistance Gene Transfer.

Antimicrobial resistance (AMR) is a relevant public health problem worldwide, and microbiome bacteria may contribute to the horizontal gene transfer associated with antimicrobial resistance. The microbiome of fecal samples from Mexican adolescents were analyzed and correlated with eating habits, and the presence of AMR genes on bacteria in the microbiome was evaluated. Fecal samples from adolescents were collected and processed to extract genomic DNA. An Illumina HiSeq 1500 system was used to determine resistance genes and the microbiome of adolescents through the amplification of gene resistance and the V3-V4 regions of RNA, respectively. Analysis of the microbiome from fecal samples taken from 18 obese, overweight, and normal-weight adolescents revealed that the Firmicutes was the most frequent phylum, followed by Bacteroidetes, Actinobacteria, Proteobacteria and Verrucomicrobia. The following species were detected as the most frequent in the samples: F. prausnitzii, P. cori, B. adolescentis, E. coli and A. muciniphila. The presence of Bacteroides, Prevotella and Ruminococcus was used to establish the enterotype; enterotype 1 was more common in women and enterotype 2 was more common in men. Twenty-nine AMR genes were found for ß-lactamases, fluoroquinolones, aminoglycosides, macrolide, lincosamides, streptogramin (MLS), tetracyclines and sulfonamides. The presence of microorganisms in fecal samples that harbor AMR genes that work against antimicrobials frequently used for the treatment of microbial infections such as b-lactams, macrolides, aminoglycosides, MLS, and tetracyclines is of great concern, as these organisms may be an important reservoir for horizontal AMR gene transfer.

Emergence and widespread circulation of a recombinant SARS-CoV-2 lineage in North America.

Although recombination is a feature of coronavirus evolution, previously detected recombinant lineages of SARS-CoV-2 have shown limited circulation thus far. Here, we present a detailed phylogenetic analysis of four SARS-CoV-2 lineages to investigate the possibility of virus recombination among them. Our analyses reveal well-supported phylogenetic differences between the Orf1ab region encoding viral non-structural proteins and the rest of the genome, including Spike (S) protein and remaining reading frames. By accounting for several deletions in NSP6, Orf3a, and S, we conclude that the B.1.628 major cluster, now designated as lineage XB, originated from a recombination event between viruses of B.1.631 and B.1.634 lineages. This scenario is supported by the spatiotemporal distribution of these lineages across the USA and Mexico during 2021, suggesting that the recombination event originated in this geographical region. This event raises important questions regarding the role and potential effects of recombination on SARS-CoV-2 evolution.

The Alpha Variant (B.1.1.7) of SARS-CoV-2 Failed to Become Dominant in Mexico.

During the coronavirus disease 2019 (COVID-19) pandemic, the emergence and rapid increase of the B.1.1.7 (Alpha) lineage of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), first identified in the United Kingdom in September 2020, was well documented in different areas of the world and became a global public health concern because of its increased transmissibility. The B.1.1.7 lineage was first detected in Mexico during December 2020, showing a slow progressive increase in its circulation frequency, which reached its maximum in May 2021 but never became predominant. In this work, we analyzed the patterns of diversity and distribution of this lineage in Mexico using phylogenetic and haplotype network analyses. Despite the reported increase in transmissibility of the B.1.1.7 lineage, in most Mexican states, it did not displace cocirculating lineages, such as B.1.1.519, which dominated the country from February to May 2021. Our results show that the states with the highest prevalence of B.1.1.7 were those at the Mexico-U.S. border. An apparent pattern of dispersion of this lineage from the northern states of Mexico toward the center or the southeast was observed in the largest transmission chains, indicating possible independent introduction events from the United States. However, other entry points cannot be excluded, as shown by multiple introduction events. Local transmission led to a few successful haplotypes with a localized distribution and specific mutations indicating sustained community transmission. IMPORTANCE The emergence and rapid increase of the B.1.1.7 (Alpha) lineage of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) throughout the world were due to its increased transmissibility. However, it did not displace cocirculating lineages in most of Mexico, particularly B.1.1.519, which dominated the country from February to May 2021. In this work, we analyzed the distribution of B.1.1.7 in Mexico using phylogenetic and haplotype network analyses. Our results show that the states with the highest prevalence of B.1.1.7 (around 30%) were those at the Mexico-U.S. border, which also exhibited the highest lineage diversity, indicating possible introduction events from the United States. Also, several haplotypes were identified with a localized distribution and specific mutations, indicating that sustained community transmission occurred in the country.

Genomic characterization of two bacteriophages (vB_EcoS-phiEc3 and vB_EcoS-phiEc4) belonging to the genus Kagunavirus with lytic activity against uropathogenic Escherichia coli.

In this study, the genomes of two lytic bacteriophages, vB_EcoS-phiEc3 and vB_EcoS-phiEc4, were sequenced and characterized using bioinformatics approaches. Whole-genome analysis showed that both phages belonged to the Kagunavirus genus, Guernseyvirinae subfamily and Siphoviridae family. Moreover, their genomes had 45, 288 bp and 44,540 bp, and G + C content of 48.42% and 50.04%, respectively. The genome of vB_EcoS-phiEc3 harbored 80 protein coding sequences (CDSs), whereas vB_EcoS-phiEc4 harbored 75 CDSs. Among them, 50 CDSs in vB_EcoS-phiEc3 and 44 CDSs in vB_EcoS-phiEc4 were considered as functional genes. Their lytic activity against multidrug-resistant uropathogenic Escherichia coli (UPEC) strains, as well as the absence of antibiotic resistance genes, lysogenic and virulence genes, enable vB_EcoS-phiEc3 and vB_EcoS-phiEc4 as a safe therapy option against UPEC infections.

Metagenomic analysis and antimicrobial activity of two fermented milk kefir samples.

In recent years, the fermented milk product kefir has been intensively studied because of its health benefits. Here, we evaluated the microbial consortia of two kefir samples, from Escarcega, Campeche, and Campeche (México). We considered a functional comparison between both samples, including fungal and bacterial inhibition; second, we applied shotgun metagenomics to assess the structure and functional diversity of the communities of microorganisms. These two samples exhibited antagonisms against bacterial and fungal pathogens. Bioactive polyketides and nonribosomal peptides were identified by LC-HRMS analysis. We also observed a high bacterial diversity and an abundance of Actinobacteria in both kefir samples, and a greater abundance of Saccharomyces species in kefir of Escarcega than in the Campeche kefir. When the prophage compositions were evaluated, the Campeche sample showed a higher diversity of prophage sequences. In Escarcega, we observed a prevalence of prophage families that infect Enterobacteria and Lactobacillus. The sequences associated with secondary metabolites, such as plipastatin, fengycin, and bacillaene, and also bacteriocins like helveticin and zoocin, were also found in different proportions, with greater diversity in the Escarcega sample. The analyses described in this work open the opportunity to understand the microbial diversity in kefir samples from two distant localities.

OXA-48 Carbapenemase in Klebsiella pneumoniae Sequence Type 307 in Ecuador.

Antibiotic resistance is on the rise, leading to an increase in morbidity and mortality due to infectious diseases. Klebsiella pneumoniae is a Gram-negative bacterium that causes bronchopneumonia, abscesses, urinary tract infection, osteomyelitis, and a wide variety of infections. The ubiquity of this microorganism confounds with the great increase in antibiotic resistance and have bred great concern worldwide. K. pneumoniae sequence type (ST) 307 is a widespread emerging clone associated with hospital-acquired infections, although sporadic community infections have also been reported. The aim of our study is to describe the first case of Klebsiella pneumoniae (ST) 307 harboring the blaOXA-48-like gene in Ecuador. We characterized a new plasmid that carry OXA-48 and could be the source of future outbreaks. The strain was recovered from a patient with cancer previously admitted in a Ukrainian hospital, suggesting that this mechanism of resistance could be imported. These findings highlight the importance of programs based on active molecular surveillance for the intercontinental spread of multidrug-resistant microorganisms with emergent carbapenemases.

Metagenomics of Atacama Lithobiontic Extremophile Life Unveils Highlights on Fungal Communities, Biogeochemical Cycles and Carbohydrate-Active Enzymes.

Halites, which are typically found in various Atacama locations, are evaporitic rocks that are considered as micro-scaled salterns. Both structural and functional metagenomic analyses of halite nodules were performed. Structural analyses indicated that the halite microbiota is mainly composed of NaCl-adapted microorganisms. In addition, halites appear to harbor a limited diversity of fungal families together with a biodiverse collection of protozoa. Functional analysis indicated that the halite microbiome possesses the capacity to make an extensive contribution to carbon, nitrogen, and sulfur cycles, but possess a limited capacity to fix nitrogen. The halite metagenome also contains a vast repertory of carbohydrate active enzymes (CAZY) with glycosyl transferases being the most abundant class present, followed by glycosyl hydrolases (GH). Amylases were also present in high abundance, with GH also being identified. Thus, the halite microbiota is a potential useful source of novel enzymes that could have biotechnological applicability. This is the first metagenomic report of fungi and protozoa as endolithobionts of halite nodules, as well as the first attempt to describe the repertoire of CAZY in this community. In addition, we present a comprehensive functional metagenomic analysis of the metabolic capacities of the halite microbiota, providing evidence for the first time on the sulfur cycle in Atacama halites.

A Review on Viral Metagenomics in Extreme Environments.

Viruses are the most abundant biological entities in the biosphere, and have the ability to infect Bacteria, Archaea, and Eukaryotes. The virome is estimated to be at least ten times more abundant than the microbiome with 107 viruses per milliliter and 109 viral particles per gram in marine waters and sediments or soils, respectively. Viruses represent a largely unexplored genetic diversity, having an important role in the genomic plasticity of their hosts. Moreover, they also play a significant role in the dynamics of microbial populations. In recent years, metagenomic approaches have gained increasing popularity in the study of environmental viromes, offering the possibility of extending our knowledge related to both virus diversity and their functional characterization. Extreme environments represent an interesting source of both microbiota and their virome due to their particular physicochemical conditions, such as very high or very low temperatures and >1 atm hydrostatic pressures, among others. Despite the fact that some progress has been made in our understanding of the ecology of the microbiota in these habitats, few metagenomic studies have described the viromes present in extreme ecosystems. Thus, limited advances have been made in our understanding of the virus community structure in extremophilic ecosystems, as well as in their biotechnological potential. In this review, we critically analyze recent progress in metagenomic based approaches to explore the viromes in extreme environments and we discuss the potential for new discoveries, as well as methodological challenges and perspectives.

Extremophile deep-sea viral communities from hydrothermal vents: Structural and functional analysis.

Ten publicly available metagenomic data sets from hydrothermal vents were analyzed to determine the taxonomic structure of the viral communities present, as well as their potential metabolic functions. The type of natural selection on two auxiliary metabolic genes was also analyzed. The structure of the virome in the hydrothermal vents was quite different in comparison with the viruses present in sediments, with specific populations being present in greater abundance in the plume samples when compared with the sediment samples. ssDNA genomes such as Circoviridae and Microviridae were predominantly present in the sediment samples, with Caudovirales which are dsDNA being present in the vent samples. Genes potentially encoding enzymes that participate in carbon, nitrogen and sulfur metabolic pathways were found in greater abundance, than those involved in the oxygen cycle, in the hydrothermal vents. Functional profiling of the viromes, resulted in the discovery of genes encoding proteins involved in bacteriophage capsids, DNA synthesis, nucleotide synthesis, DNA repair, as well as viral auxiliary metabolic genes such as cytitidyltransferase and ribonucleotide reductase. These auxiliary metabolic genes participate in the synthesis of phospholipids and nucleotides respectively and are likely to contribute to enhancing the fitness of their bacterial hosts within the hydrothermal vent communities. Finally, evolutionary analysis suggested that these auxiliary metabolic genes are highly conserved and evolve under purifying selection, and are thus maintained in their genome.