Deciphering the Potential Coding of Human Cytomegalovirus: New Predicted Transmembrane Proteome.

CMV is a major cause of morbidity and mortality in immunocompromised individuals that will benefit from the availability of a vaccine. Despite the efforts made during the last decade, no CMV vaccine is available. An ideal CMV vaccine should elicit a broad immune response against multiple viral antigens including proteins involved in virus-cell interaction and entry. However, the therapeutic use of neutralizing antibodies targeting glycoproteins involved in viral entry achieved only partial protection against infection. In this scenario, a better understanding of the CMV proteome potentially involved in viral entry may provide novel candidates to include in new potential vaccine design. In this study, we aimed to explore the CMV genome to identify proteins with putative transmembrane domains to identify new potential viral envelope proteins. We have performed in silico analysis using the genome sequences of nine different CMV strains to predict the transmembrane domains of the encoded proteins. We have identified 77 proteins with transmembrane domains, 39 of which were present in all the strains and were highly conserved. Among the core proteins, 17 of them such as UL10, UL139 or US33A have no ascribed function and may be good candidates for further mechanistic studies.

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.

Recombination events and variability among full-length genomes of co-circulating molluscum contagiosum virus subtypes 1 and 2.

Molluscum contagiosum virus (MCV) is the sole member of the Molluscipoxvirus genus and causes a highly prevalent human disease of the skin characterized by the formation of a variable number of lesions that can persist for prolonged periods of time. Two major genotypes, subtype 1 and subtype 2, are recognized, although currently only a single complete genomic sequence corresponding to MCV subtype 1 is available. Using next-generation sequencing techniques, we report the complete genomic sequence of four new MCV isolates, including the first one derived from a subtype 2. Comparisons suggest a relatively distant evolutionary split between both MCV subtypes. Further, our data illustrate concurrent circulation of distinct viruses within a population and reveal the existence of recombination events among them. These results help identify a set of MCV genes with potentially relevant roles in molluscum contagiosum epidemiology and pathogenesis.

Classification of protein motifs based on subcellular localization uncovers evolutionary relationships at both sequence and functional levels.

BACKGROUND: Most proteins have evolved in specific cellular compartments that limit their functions and potential interactions. On the other hand, motifs define amino acid arrangements conserved between protein family members and represent powerful tools for assigning function to protein sequences. The ideal motif would identify all members of a protein family but in practice many motifs identify both family members and unrelated proteins, referred to as True Positive (TP) and False Positive (FP) sequences, respectively. RESULTS: To address the relationship between protein motifs, protein function and cellular localization, we systematically assigned subcellular localization data to motif sequences from the comprehensive PROSITE sequence motif database. Using this data we analyzed relationships between localization and function. We find that TPs and FPs have a strong tendency to localize in different compartments. When multiple localizations are considered, TPs are usually distributed between related cellular compartments. We also identified cases where FPs are concentrated in particular subcellular regions, indicating possible functional or evolutionary relationships with TP sequences of the same motif. CONCLUSIONS: Our findings suggest that the systematic examination of subcellular localization has the potential to uncover evolutionary and functional relationships between motif-containing sequences. We believe that this type of analysis complements existing motif annotations and could aid in their interpretation. Our results shed light on the evolution of cellular organelles and potentially establish the basis for new subcellular localization and function prediction algorithms.

Leveraging phylogenetic signal to unravel microbiome function and assembly rules.

Clarifying the general rules behind microbial community assembly will foster the development of microbiome-based technological solutions. Here, we study microbial community assembly through a computational analysis of phylogenetic core groups (PCGs): discrete portions of the bacterial phylogeny with high prevalence in the ecosystem under study. We first show that the existence of PCGs was a predominant feature of the varied set of microbial ecosystems studied. Then, we re-analyzed an in vitro experimental dataset using a PCG-based approach, drawing only from its community composition data and from publicly available genomic databases. Using mainly genome scale metabolic models and population dynamics modeling, we obtained ecological insights on metabolic niche structure and population dynamics comparable to those gained after canonical experimentation. Thus, leveraging phylogenetic signal to help unravel microbiome function and assembly rules offers a potential avenue to gain further insight on Earth's microbial ecosystems.

Latent antibiotic resistance genes are abundant, diverse, and mobile in human, animal, and environmental microbiomes.

BACKGROUND: Bacterial communities in humans, animals, and the external environment maintain a large collection of antibiotic resistance genes (ARGs). However, few of these ARGs are well-characterized and thus established in existing resistance gene databases. In contrast, the remaining latent ARGs are typically unknown and overlooked in most sequencing-based studies. Our view of the resistome and its diversity is therefore incomplete, which hampers our ability to assess risk for promotion and spread of yet undiscovered resistance determinants. RESULTS: A reference database consisting of both established and latent ARGs (ARGs not present in current resistance gene repositories) was created. By analyzing more than 10,000 metagenomic samples, we showed that latent ARGs were more abundant and diverse than established ARGs in all studied environments, including the human- and animal-associated microbiomes. The pan-resistomes, i.e., all ARGs present in an environment, were heavily dominated by latent ARGs. In comparison, the core-resistome, i.e., ARGs that were commonly encountered, comprised both latent and established ARGs. We identified several latent ARGs shared between environments and/or present in human pathogens. Context analysis of these genes showed that they were located on mobile genetic elements, including conjugative elements. We, furthermore, identified that wastewater microbiomes had a surprisingly large pan- and core-resistome, which makes it a potentially high-risk environment for the mobilization and promotion of latent ARGs. CONCLUSIONS: Our results show that latent ARGs are ubiquitously present in all environments and constitute a diverse reservoir from which new resistance determinants can be recruited to pathogens. Several latent ARGs already had high mobile potential and were present in human pathogens, suggesting that they may constitute emerging threats to human health. We conclude that the full resistome-including both latent and established ARGs-needs to be considered to properly assess the risks associated with antibiotic selection pressures. Video Abstract.

Extensive screening reveals previously undiscovered aminoglycoside resistance genes in human pathogens.

Antibiotic resistance is a growing threat to human health, caused in part by pathogens accumulating antibiotic resistance genes (ARGs) through horizontal gene transfer. New ARGs are typically not recognized until they have become widely disseminated, which limits our ability to reduce their spread. In this study, we use large-scale computational screening of bacterial genomes to identify previously undiscovered mobile ARGs in pathogens. From ~1 million genomes, we predict 1,071,815 genes encoding 34,053 unique aminoglycoside-modifying enzymes (AMEs). These cluster into 7,612 families (<70% amino acid identity) of which 88 are previously described. Fifty new AME families are associated with mobile genetic elements and pathogenic hosts. From these, 24 of 28 experimentally tested AMEs confer resistance to aminoglycoside(s) in Escherichia coli, with 17 providing resistance above clinical breakpoints. This study greatly expands the range of clinically relevant aminoglycoside resistance determinants and demonstrates that computational methods enable early discovery of potentially emerging ARGs.

Lactoferrin modulates gut microbiota and Toll-like receptors (TLRs) in mice with dysbiosis induced by antibiotics.

Background: Antibiotic administration can result in gut microbiota and immune system alterations that impact health. Bovine lactoferrin is a milk protein with anticancer, anti-inflammatory, antimicrobial and immune modulatory activities. The aim was to study the ability of native and iron-saturated lactoferrin to reverse the effects of clindamycin on gut microbiota and intestinal Toll-like receptor (TLR) expression in a murine model. Methods: Male C57BL/6 mice were treated with vehicle, clindamycin (Clin), native bovine lactoferrin (nLf), nLf + clindamycin (nLf_Clin), iron-saturated bovine lactoferrin (sLf) and sLf + clindamycin (sLf_Clin). Fecal samples of each group were collected, and bacterial DNA was extracted. Sequencing of 16s rRNA V4 hypervariable gene regions was conducted to assess the microbial composition. mRNA expression levels of TLRs (1-9) were determined in mouse colon by qPCR. Pearson's correlation test was carried out between bacteria showing differences in abundance among samples and TLR2, TLR8 and TLR9. Results: Beta-diversity analysis showed that the microbial community of the vehicle was different from the communities of Clin, nLf_Clin and sLf_Clin. At the family level, Bacteroidaceae, Prevotellaceae and Rikenellaceae decreased in the Clin group, and treatment with nLf or sLf reverted these effects. Clin reduced the expression of TLR2, TLR8 and TLR9 and sLf reverted the decrease in the expression of these receptors. Finally, TLR8 was positively correlated with Rikenellaceae abundance. Conclusion: In a situation of intestinal dysbiosis induced by clindamycin, lactoferrin restores the normal levels of some anti-inflammatory bacteria and TLRs and, therefore, could be a good ingredient to be added to functional foods.

Large-scale characterization of the macrolide resistome reveals high diversity and several new pathogen-associated genes.

Macrolides are broad-spectrum antibiotics used to treat a range of infections. Resistance to macrolides is often conferred by mobile resistance genes encoding Erm methyltransferases or Mph phosphotransferases. New erm and mph genes keep being discovered in clinical settings but their origins remain unknown, as is the type of macrolide resistance genes that will appear in the future. In this study, we used optimized hidden Markov models to characterize the macrolide resistome. Over 16 terabases of genomic and metagenomic data, representing a large taxonomic diversity (11 030 species) and diverse environments (1944 metagenomic samples), were searched for the presence of erm and mph genes. From this data, we predicted 28 340 macrolide resistance genes encoding 2892 unique protein sequences, which were clustered into 663 gene families (<70 % amino acid identity), of which 619 (94 %) were previously uncharacterized. This included six new resistance gene families, which were located on mobile genetic elements in pathogens. The function of ten predicted new resistance genes were experimentally validated in Escherichia coli using a growth assay. Among the ten tested genes, seven conferred increased resistance to erythromycin, with five genes additionally conferring increased resistance to azithromycin, showing that our models can be used to predict new functional resistance genes. Our analysis also showed that macrolide resistance genes have diverse origins and have transferred horizontally over large phylogenetic distances into human pathogens. This study expands the known macrolide resistome more than ten-fold, provides insights into its evolution, and demonstrates how computational screening can identify new resistance genes before they become a significant clinical problem.

Assessment of phylo-functional coherence along the bacterial phylogeny and taxonomy.

In this report we use available curated phylogenies, taxonomy, and genome annotations to assess the phylogenetic and gene content similarity associated with each different taxon and taxonomic rank. Subsequently, we employ the same data to assess the frontiers of functional coherence along the bacterial phylogeny. Our results show that within-group phylogenetic and gene content similarity of taxa in the same rank are not homogenous, and that these values show extensive overlap between ranks. Functional coherence along the 16S rRNA gene-based phylogeny was limited to 44 particular nodes presenting large variations in phylogenetic depth. For instance, the deep subtree affiliated to class Actinobacteria presented functional coherence, while the shallower family Enterobacteriaceae-affiliated subtree did not. On the other hand, functional coherence along the genome-based phylogeny delimited deep subtrees affiliated to phyla Actinobacteriota, Deinococcota, Chloroflexota, Firmicutes, and a subtree containing the rest of the bacterial phyla. The results presented here can be used to guide the exploration of results in many microbial ecology and evolution research scenarios. Moreover, we provide dedicated scripts and files that can be used to continue the exploration of functional coherence along the bacterial phylogeny employing different parameters or input data ( https://git.io/Jec5U ).

A comprehensive human minimal gut metagenome extends the host's metabolic potential.

Accumulating evidence suggests that humans could be considered as holobionts in which the gut microbiota play essential functions. Initial metagenomic studies reported a pattern of shared genes in the gut microbiome of different individuals, leading to the definition of the minimal gut metagenome as the set of microbial genes necessary for homeostasis and present in all healthy individuals. This study analyses the minimal gut metagenome of the most comprehensive dataset available, including individuals from agriculturalist and industrialist societies, also embodying highly diverse ethnic and geographical backgrounds. The outcome, based on metagenomic predictions for community composition data, resulted in a minimal metagenome comprising 3412 genes, mapping to 1856 reactions and 128 metabolic pathways predicted to occur across all individuals. These results were substantiated by the analysis of two additional datasets describing the microbial community compositions of larger Western cohorts, as well as a substantial shotgun metagenomics dataset. Subsequent analyses showed the plausible metabolic complementarity provided by the minimal gut metagenome to the human genome.

Methods for Enrichment and Sequencing of Oral Viral Assemblages: Saliva, Oral Mucosa, and Dental Plaque Viromes.

The oral cavity is a major portal of entry for human pathogens including viruses. However, metagenomics has revealed that highly personalized and time-persistent bacteriophage assemblages dominate this habitat. Most oral bacteriophages follow lysogenic life cycles, deploying complex strategies to manage bacterial homeostasis. Although bacterial dysbiosis underlies common oral pathologies such as caries and periodontitis, the cause of these bacteria replacements remains obscure, and it is theorized that bacteriophages play an important role. The enormous sensitivity of metagenomics coupled with next-generation sequencing has made technically feasible to address the putative role of bacteriophages in oral dysbiosis and represents a valuable tool to discover new human viruses.This chapter proposes a workflow that consists of a simple viral enrichment protocol, two alternative random amplification methods, and next-generation sequencing to access virome composition in three oral environments: supragingival plaque, saliva, and mucosa. These protocols circumvent some well-known sources of bias, providing genomic information about DNA and RNA viral communities with minimal contamination from human and bacterial sources.

Evaluation of bias induced by viral enrichment and random amplification protocols in metagenomic surveys of saliva DNA viruses.

BACKGROUND: Viruses are key players regulating microbial ecosystems. Exploration of viral assemblages is now possible thanks to the development of metagenomics, the most powerful tool available for studying viral ecology and discovering new viruses. Unfortunately, several sources of bias lead to the misrepresentation of certain viruses within metagenomics workflows, hindering the shift from merely descriptive studies towards quantitative comparisons of communities. Therefore, benchmark studies on virus enrichment and random amplification protocols are required to better understand the sources of bias. RESULTS: We assessed the bias introduced by viral enrichment on mock assemblages composed of seven DNA viruses, and the bias from random amplification methods on human saliva DNA viromes, using qPCR and deep sequencing, respectively. While iodixanol cushions and 0.45 µm filtration preserved the original composition of nuclease-protected viral genomes, low-force centrifugation and 0.22 µm filtration removed large viruses. Comparison of unamplified and randomly amplified saliva viromes revealed that multiple displacement amplification (MDA) induced stochastic bias from picograms of DNA template. However, the type of bias shifted to systematic using 1 ng, with only a marginal influence by amplification time. Systematic bias consisted of over-amplification of small circular genomes, and under-amplification of those with extreme GC content, a negative bias that was shared with the PCR-based sequence-independent, single-primer amplification (SISPA) method. MDA based on random priming provided by a DNA primase activity slightly outperformed those based on random hexamers and SISPA, which may reflect differences in ability to handle sequences with extreme GC content. SISPA viromes showed uneven coverage profiles, with high coverage peaks in regions with low linguistic sequence complexity. Despite misrepresentation of certain viruses after random amplification, ordination plots based on dissimilarities among contig profiles showed perfect overlapping of related amplified and unamplified saliva viromes and strong separation from unrelated saliva viromes. This result suggests that random amplification bias has a minor impact on beta diversity studies. CONCLUSIONS: Benchmark analyses of mock and natural communities of viruses improve understanding and mitigate bias in metagenomics surveys. Bias induced by random amplification methods has only a minor impact on beta diversity studies of human saliva viromes.

Genome sequence of two novel species of torque teno minivirus from the human oral cavity.

Anelloviridae is a family of circular, single-stranded DNA viruses highly prevalent among humans. We report the genome sequence of two torque teno miniviruses found in human oral mucosa samples. Genome organization, phylogenetic analysis, and pairwise comparisons reveal that they belong to novel species within the Betatorquevirus genus.