Search

The extraordinary diversity of viruses infecting bacteria and archaea is now primarily studied through metagenomics. While metagenomes enable high-throughput exploration of the viral sequence space, metagenome-derived sequences lack key information compared to isolated viruses, in particular host association. Different computational approaches are available to predict the host(s) of uncultivated viruses based on their genome sequences, but thus far individual approaches are limited either in precision or in recall, i.e., for a number of viruses they yield erroneous predictions or no prediction at all. Here, we describe iPHoP, a two-step framework that integrates multiple methods to reliably predict host taxonomy at the genus rank for a broad range of viruses infecting bacteria and archaea, while retaining a low false discovery rate. Based on a large dataset of metagenome-derived virus genomes from the IMG/VR database, we illustrate how iPHoP can provide extensive host prediction and guide further characterization of uncultivated viruses. The extraordinary diversity of viruses infecting bacteria and archaea is now primarily studied through metagenomics, but information regarding viral hosts is lost. This study presents iPHoP (integrated Phage-Host Prediction), a new tool that leverages multiple complementary approaches to maximize host prediction for metagenome-assembled virus genomes. [ABSTRACT FROM AUTHOR]

A universal taxonomy of viruses is essential for a comprehensive view of the virus world and for communicating the complicated evolutionary relationships among viruses. However, there are major differences in the conceptualisation and approaches to virus classification and nomenclature among virologists, clinicians, agronomists, and other interested parties. Here, we provide recommendations to guide the construction of a coherent and comprehensive virus taxonomy, based on expert scientific consensus. Firstly, assignments of viruses should be congruent with the best attainable reconstruction of their evolutionary histories, i.e., taxa should be monophyletic. This fundamental principle for classification of viruses is currently included in the International Committee on Taxonomy of Viruses (ICTV) code only for the rank of species. Secondly, phenotypic and ecological properties of viruses may inform, but not override, evolutionary relatedness in the placement of ranks. Thirdly, alternative classifications that consider phenotypic attributes, such as being vector-borne (e.g., "arboviruses"), infecting a certain type of host (e.g., "mycoviruses," "bacteriophages") or displaying specific pathogenicity (e.g., "human immunodeficiency viruses"), may serve important clinical and regulatory purposes but often create polyphyletic categories that do not reflect evolutionary relationships. Nevertheless, such classifications ought to be maintained if they serve the needs of specific communities or play a practical clinical or regulatory role. However, they should not be considered or called taxonomies. Finally, while an evolution-based framework enables viruses discovered by metagenomics to be incorporated into the ICTV taxonomy, there are essential requirements for quality control of the sequence data used for these assignments. Combined, these four principles will enable future development and expansion of virus taxonomy as the true evolutionary diversity of viruses becomes apparent. Transforming an existing phenotypic classification of viruses into one based on evolutionary relationships that can accommodate the vast number of viruses characterized in metagenomics and environmental studies is an ongoing challenge. This Consensus View explains how such a taxonomy can be expanded to encapsulate viral diversity and to recognize independent biological origins of different virus groups. [ABSTRACT FROM AUTHOR]

Background The HIV-1 subtype B epidemic amongst men who have sex with men (MSM) is resurgent in many countries despite the widespread use of effective combination antiretroviral therapy (cART). In this combined mathematical and phylogenetic study of observational data, we aimed to find out the extent to which the resurgent epidemic is the result of newly introduced strains or of growth of already circulating strains. Methods and Findings As of November 2011, the ATHENA observational HIV cohort of all patients in care in the Netherlands since 1996 included HIV-1 subtype Bpolymerase sequences from 5,852 patients. Patients who were diagnosed between 1981 and 1995 were included in the cohort if they were still alive in 1996. The ten most similar sequences to each ATHENA sequence were selected from the Los Alamos HIV Sequence Database, and a phylogenetic tree was created of a total of 8,320 sequences. Large transmission clusters that included ≥ 10 ATHENA sequences were selected, with a local support value ≥ 0.9 and median pairwise patristic distance below the fifth percentile of distances in the whole tree. Time-varying reproduction numbers of the large MSM-majority clusters were estimated through mathematical modeling. We identified 106 large transmission clusters, including 3,061 (52%) ATHENA and 652 Los Alamos sequences. Half of the HIV sequences from MSM registered in the cohort in the Netherlands (2,128 of 4,288) were included in 91 large MSMmajority clusters. Strikingly, at least 54 (59%) of these 91 MSM-majority clusters were already circulating before 1996, when cART was introduced, and have persisted to the present. Overall, 1,226 (35%) of the 3,460 diagnoses among MSM since 1996 were found in these 54 long-standing clusters. The reproduction numbers of all large MSM-majority clusters were around the epidemic threshold value of one over the whole study period. A tendency towards higher numbers was visible in recent years, especially in the more recently introduced clusters. The mean age of MSM at diagnosis increased by 0.45 years/ year within clusters, but new clusters appeared with lower mean age. Major strengths of this study are the high proportion of HIV-positive MSM with a sequence in this study and the combined application of phylogenetic and modeling approaches. Main limitations are the assumption that the sampled population is representative of the overall HIV-positive population and the assumption that the diagnosis interval distribution is similar between clusters. Conclusions The resurgent HIV epidemic amongst MSM in the Netherlands is driven by several large, persistent, self-sustaining, and, in many cases, growing sub-epidemics shifting towards new generations of MSM. Many of the sub-epidemics have been present since the early epidemic, to which new sub-epidemics are being added., Introduction The HIV-1 epidemic amongst men who have sex with men (MSM) is resurgent in many Western countries [1-5]. This may seem paradoxical, as HIV-positive men are being diagnosed at [...]

Peronospora effusa (previously known as P. farinosa f. sp. spinaciae, and here referred to as Pfs) is an obligate biotrophic oomycete that causes downy mildew on spinach (Spinacia oleracea). To combat this destructive many disease resistant cultivars have been bred and used. However, new Pfs races rapidly break the employed resistance genes. To get insight into the gene repertoire of Pfs and identify infection-related genes, the genome of the first reference race, Pfs1, was sequenced, assembled, and annotated. Due to the obligate biotrophic nature of this pathogen, material for DNA isolation can only be collected from infected spinach leaves that, however, also contain many other microorganisms. The obtained sequences can, therefore, be considered a metagenome. To filter and obtain Pfs sequences we utilized the CAT tool to taxonomically annotate ORFs residing on long sequences of a genome pre-assembly. This study is the first to show that CAT filtering performs well on eukaryotic contigs. Based on the taxonomy, determined on multiple ORFs, contaminating long sequences and corresponding reads were removed from the metagenome. Filtered reads were re-assembled to provide a clean and improved Pfs genome sequence of 32.4 Mbp consisting of 8,635 scaffolds. Transcript sequencing of a range of infection time points aided the prediction of a total of 13,277 gene models, including 99 RxLR(-like) effector, and 14 putative Crinkler genes. Comparative analysis identified common features in the predicted secretomes of different obligate biotrophic oomycetes, regardless of their phylogenetic distance. Their secretomes are generally smaller, compared to hemi-biotrophic and necrotrophic oomycete species. We observe a reduction in proteins involved in cell wall degradation, in Nep1-like proteins (NLPs), proteins with PAN/apple domains, and host translocated effectors. The genome of Pfs1 will be instrumental in studying downy mildew virulence and for understanding the molecular adaptations by which new isolates break spinach resistance. [ABSTRACT FROM AUTHOR]

Teredinidae are a family of highly adapted wood-feeding and wood-boring bivalves, commonly known as shipworms, whose evolution is linked to the acquisition of cellulolytic gammaproteobacterial symbionts harbored in bacteriocytes within the gills. In the present work we applied metagenomics to characterize microbiomes of the gills and digestive tract of Neoteredo reynei, a mangrove-adapted shipworm species found over a large range of the Brazilian coast. Comparative metagenomics grouped the gill symbiont community of different N. reynei specimens, indicating closely related bacterial types are shared. Similarly, the intestine and digestive gland communities were related, yet were more diverse than and showed no overlap with the gill community. Annotation of assembled metagenomic contigs revealed that the gill symbiotic community of N. reynei encodes a plethora of plant cell wall polysaccharides degrading glycoside hydrolase encoding genes, and Biosynthetic Gene Clusters (BGCs). In contrast, the digestive tract microbiomes seem to play little role in wood digestion and secondary metabolites biosynthesis. Metagenome binning recovered the nearly complete genome sequences of two symbiotic Teredinibacter strains from the gills, a representative of Teredinibacter turnerae “clade I” strain, and a yet to be cultivated Teredinibacter sp. type. These Teredinibacter genomes, as well as un-binned gill-derived gammaproteobacteria contigs, also include an endo-β-1,4-xylanase/acetylxylan esterase multi-catalytic carbohydrate-active enzyme, and a trans-acyltransferase polyketide synthase (trans-AT PKS) gene cluster with the gene cassette for generating β-branching on complex polyketides. Finally, we use multivariate analyses to show that the secondary metabolome from the genomes of Teredinibacter representatives, including genomes binned from N. reynei gills’ metagenomes presented herein, stands out within the Cellvibrionaceae family by size, and enrichments for polyketide, nonribosomal peptide and hybrid BGCs. Results presented here add to the growing characterization of shipworm symbiotic microbiomes and indicate that the N. reynei gill gammaproteobacterial community is a prolific source of biotechnologically relevant enzymes for wood-digestion and bioactive compounds production. [ABSTRACT FROM AUTHOR]

The Salmonella enterica serovars Enteritidis, Dublin, and Gallinarum are closely related but differ in virulence and host range. To identify the genetic elements responsible for these differences and to better understand how these serovars are evolving, we sequenced the genomes of Enteritidis strain LK5 and Dublin strain SARB12 and compared these genomes to the publicly available Enteritidis P125109, Dublin CT 02021853 and Dublin SD3246 genome sequences. We also compared the publicly available Gallinarum genome sequences from biotype Gallinarum 287/91 and Pullorum RKS5078. Using bioinformatic approaches, we identified single nucleotide polymorphisms, insertions, deletions, and differences in prophage and pseudogene content between strains belonging to the same serovar. Through our analysis we also identified several prophage cargo genes and pseudogenes that affect virulence and may contribute to a host-specific, systemic lifestyle. These results strongly argue that the Enteritidis, Dublin and Gallinarum serovars of Salmonella enterica evolve by acquiring new genes through horizontal gene transfer, followed by the formation of pseudogenes. The loss of genes necessary for a gastrointestinal lifestyle ultimately leads to a systemic lifestyle and niche exclusion in the host-specific serovars. [ABSTRACT FROM AUTHOR]

Progress in extending the life expectancy of cystic fibrosis (CF) patients remains jeopardized by the increasing incidence of fungal respiratory infections. Pseudallescheria boydii (P. boydii), an emerging pathogen of humans, is a filamentous fungus frequently isolated from the respiratory secretions of CF patients. It is commonly believed that infection by this fungus occurs through inhalation of airborne conidia, but the mechanisms allowing the adherence of Pseudallescheria to the host epithelial cells and its escape from the host immune defenses remain largely unknown. Given that the cell wall orchestrates all these processes, we were interested in studying its dynamic changes in conidia as function of the age of cultures. We found that the surface hydrophobicity and electronegative charge of conidia increased with the age of culture. Melanin that can influence the cell surface properties, was extracted from conidia and estimated using UV-visible spectrophotometry. Cells were also directly examined and compared using electron paramagnetic resonance (EPR) that determines the production of free radicals. Consistent with the increased amount of melanin, the EPR signal intensity decreased suggesting polymerization of melanin. These results were confirmed by flow cytometry after studying the effect of melanin polymerization on the surface accessibility of mannose-containing glycoconjugates to fluorescent concanavalin A. In the absence of melanin, conidia showed a marked increase in fluorescence intensity as the age of culture increased. Using atomic force microscopy, we were unable to find rodlet-forming hydrophobins, molecules that can also affect conidial surface properties. In conclusion, the changes in surface properties and biochemical composition of the conidial wall with the age of culture highlight the process of conidial maturation. Mannose-containing glycoconjugates that are involved in immune recognition, are progressively masked by polymerization of melanin, an antioxidant that is commonly thought to allow fungal escape from the host immune defenses. [ABSTRACT FROM AUTHOR]

As part of a virus discovery investigation using a metagenomic approach, a highly divergent novel Human papillomavirus type was identified in pooled convenience nasal/oropharyngeal swab samples collected from patients with febrile respiratory illness. Phylogenetic analysis of the whole genome and the L1 gene reveals that the new HPV identified in this study clusters with previously described gamma papillomaviruses, sharing only 61.1% (whole genome) and 63.1% (L1) sequence identity with its closest relative in the Papillomavirus episteme (PAVE) database. This new virus was named HPV_SD2 pending official classification. The complete genome of HPV-SD2 is 7,299 bp long (36.3% G/C) and contains 7 open reading frames (L2, L1, E6, E7, E1, E2 and E4) and a non-coding long control region (LCR) between L1 and E6. The metagenomic procedures, coupled with the bioinformatic methods described herein are well suited to detect small circular genomes such as those of human papillomaviruses. [ABSTRACT FROM AUTHOR]

Azole compounds are the primary therapy for patients with diseases caused by Aspergillus fumigatus. However, prolonged treatment may cause resistance to develop, which is associated with treatment failure. The azole target cyp51A is a hotspot for mutations that confer phenotypic resistance, but in an increasing number of resistant isolates the underlying mechanism remains unknown. Here, we report the discovery of a novel resistance mechanism, caused by a mutation in the CCAAT-binding transcription factor complex subunit HapE. From one patient, fourA. fumigatus isolates were serially collected. The last two isolates developed an azole resistant phenotype during prolonged azole therapy. Because the resistant isolates contained a wild type cyp51A gene and the isolates were isogenic, the complete genomes of the last susceptible isolate and the first resistant isolate (taken 17 weeks apart) were sequenced using Illumina technology to identify the resistance conferring mutation. By comparing the genome sequences to each other as well as to two A. fumigatus reference genomes, several potential non-synonymous mutations in protein-coding regions were identified, six of which could be confirmed by PCR and Sanger sequencing. Subsequent sexual crossing experiments showed that resistant progeny always contained a P88L substitution in HapE, while the presence of the other five mutations did not correlate with resistance in the progeny. Cloning the mutated hapE gene into the azole susceptible akuBKU80 strain showed that the HapE P88L mutation by itself could confer the resistant phenotype. This is the first time that whole genome sequencing and sexual crossing strategies have been used to find the genetic basis of a trait of interest in A. fumigatus. The discovery may help understand alternate pathways for azole resistance in A. fumigatus with implications for the molecular diagnosis of resistance and drug discovery. [ABSTRACT FROM AUTHOR]

The endemic marine sponge Arenosclera brasiliensis (Porifera, Demospongiae, Haplosclerida) is a known source of secondary metabolites such as arenosclerins A-C. In the present study, we established the composition of the A. brasiliensis microbiome and the metabolic pathways associated with this community. We used 454 shotgun pyrosequencing to generate approximately 640,000 high-quality sponge-derived sequences (,150 Mb). Clustering analysis including sponge, seawater and twenty-three other metagenomes derived from marine animal microbiomes shows that A. brasiliensis contains a specific microbiome. Fourteen bacterial phyla (including Proteobacteria, Cyanobacteria, Actinobacteria, Bacteroidetes, Firmicutes and Cloroflexi) were consistently found in the A. brasiliensis metagenomes. The A. brasiliensis microbiome is enriched for Betaproteobacteria (e.g., Burkholderia) and Gammaproteobacteria (e.g., Pseudomonas and Alteromonas) compared with the surrounding planktonic microbial communities. Functional analysis based on Rapid Annotation using Subsystem Technology (RAST) indicated that the A. brasiliensis microbiome is enriched for sequences associated with membrane transport and one-carbon metabolism. In addition, there was an overrepresentation of sequences associated with aerobic and anaerobic metabolism as well as the synthesis and degradation of secondary metabolites. This study represents the first analysis of sponge-associated microbial communities via shotgun pyrosequencing, a strategy commonly applied in similar analyses in other marine invertebrate hosts, such as corals and algae. We demonstrate that A. brasiliensis has a unique microbiome that is distinct from that of the surrounding planktonic microbes and from other marine organisms, indicating a species-specific microbiome. [ABSTRACT FROM AUTHOR]

The 7th cholera pandemic reached Latin America in 1991, spreading from Peru to virtually all Latin American countries. During the late epidemic period, a strain that failed to ferment sucrose dominated cholera outbreaks in the Northern Brazilian Amazon region. In order to understand the genomic characteristics and the determinants of this altered sucrose fermenting phenotype, the genome of the strain IEC224 was sequenced. This paper reports a broad genomic study of this strain, showing its correlation with the major epidemic lineage. The potentially mobile genomic regions are shown to possess GC content deviation, and harbor the main V. cholera virulence genes. A novel bioinformatic approach was applied in order to identify the putative functions of hypothetical proteins, and was compared with the automatic annotation by RAST. The genome of a large bacteriophage was found to be integrated to the IEC224's alanine aminopeptidase gene. The presence of this phage is shown to be a common characteristic of the El Tor strains from the Latin American epidemic, as well as its putative ancestor from Angola. The defective sucrose fermenting phenotype is shown to be due to a single nucleotide insertion in the V. cholerae sucrose-specific transportation gene. This frame-shift mutation truncated a membrane protein, altering its structural pore-like conformation. Further, the identification of a common bacteriophage reinforces both the monophyletic and African-Origin hypotheses for the main causative agent of the 1991 Latin America cholera epidemics. [ABSTRACT FROM AUTHOR]

Multiple factors drive the progression from healthy mucosa towards sporadic colorectal carcinomas and accumulating evidence associates intestinal bacteria with disease initiation and progression. Therefore, the aim of this study was to provide a first high-resolution map of colonic dysbiosis that is associated with human colorectal cancer (CRC). To this purpose, the microbiomes colonizing colon tumor tissue and adjacent non-malignant mucosa were compared by deep rRNA sequencing. The results revealed striking differences in microbial colonization patterns between these two sites. Although inter-individual colonization in CRC patients was variable, tumors consistently formed a niche for Coriobacteria and other proposed probiotic bacterial species, while potentially pathogenic Enterobacteria were underrepresented in tumor tissue. As the intestinal microbiota is generally stable during adult life, these findings suggest that CRC-associated physiological and metabolic changes recruit tumor-foraging commensal-like bacteria. These microbes thus have an apparent competitive advantage in the tumor microenvironment and thereby seem to replace pathogenic bacteria that may be implicated in CRC etiology. This first glimpse of the CRC microbiome provides an important step towards full understanding of the dynamic interplay between intestinal microbial ecology and sporadic CRC, which may provide important leads towards novel microbiome-related diagnostic tools and therapeutic interventions. [ABSTRACT FROM AUTHOR]

The mouse is a valuable model for unravelling the role of hepcidin in iron homeostasis, however, such studies still report hepcidin mRNA levels as a surrogate marker for bioactive hepcidin in its pivotal function to block ferroportin-mediated iron transport. Here, we aimed to assess bioactive mouse Hepcidin-1 (Hep-1) and its paralogue Hepcidin-2 (Hep-2) at the peptide level. To this purpose, fourier transform ion cyclotron resonance (FTICR) and tandem-MS was used for hepcidin identification, after which a time-of-flight (TOF) MS-based methodology was exploited to routinely determine Hep-1 and -2 levels in mouse serum and urine. This method was biologically validated by hepcidin assessment in: i) 3 mouse strains (C57Bl/6; DBA/2 and BABL/c) upon stimulation with intravenous iron and LPS, ii) homozygous Hfe knock out, homozygous transferrin receptor 2 (Y245X) mutated mice and double affected mice, and iii) mice treated with a sublethal hepatotoxic dose of paracetamol. The results showed that detection of Hep-1 was restricted to serum, whereas Hep-2 and its presumed isoforms were predominantly present in urine. Elevations in serum Hep-1 and urine Hep-2 upon intravenous iron or LPS were only moderate and varied considerably between mouse strains. Serum Hep-1 was decreased in all three hemochromatosis models, being lowest in the double affected mice. Serum Hep-1 levels correlated with liver hepcidin-1 gene expression, while acute liver damage by paracetamol depleted Hep-1 from serum. Furthermore, serum Hep-1 appeared to be an excellent indicator of splenic iron accumulation. In conclusion, Hep-1 and Hep-2 peptide responses in experimental mouse agree with the known biology of hepcidin mRNA regulators, and their measurement can now be implemented in experimental mouse models to provide novel insights in post-transcriptional regulation, hepcidin function, and kinetics. [ABSTRACT FROM AUTHOR]

Progress in extending the life expectancy of cystic fibrosis (CF) patients remains jeopardized by the increasing incidence of fungal respiratory infections. Pseudallescheria boydii (P. boydii), an emerging pathogen of humans, is a filamentous fungus frequently isolated from the respiratory secretions of CF patients. It is commonly believed that infection by this fungus occurs through inhalation of airborne conidia, but the mechanisms allowing the adherence of Pseudallescheria to the host epithelial cells and its escape from the host immune defenses remain largely unknown. Given that the cell wall orchestrates all these processes, we were interested in studying its dynamic changes in conidia as function of the age of cultures. We found that the surface hydrophobicity and electronegative charge of conidia increased with the age of culture. Melanin that can influence the cell surface properties, was extracted from conidia and estimated using UV-visible spectrophotometry. Cells were also directly examined and compared using electron paramagnetic resonance (EPR) that determines the production of free radicals. Consistent with the increased amount of melanin, the EPR signal intensity decreased suggesting polymerization of melanin. These results were confirmed by flow cytometry after studying the effect of melanin polymerization on the surface accessibility of mannose-containing glycoconjugates to fluorescent concanavalin A. In the absence of melanin, conidia showed a marked increase in fluorescence intensity as the age of culture increased. Using atomic force microscopy, we were unable to find rodlet-forming hydrophobins, molecules that can also affect conidial surface properties. In conclusion, the changes in surface properties and biochemical composition of the conidial wall with the age of culture highlight the process of conidial maturation. Mannose-containing glycoconjugates that are involved in immune recognition, are progressively masked by polymerization of melanin, an antioxidant that is commonly thought to allow fungal escape from the host immune defenses. [ABSTRACT FROM AUTHOR]

Heterozygous mutations in p63 are associated with split hand/foot malformations (SHFM), orofacial clefting, and ectodermal abnormalities. Elucidation of the p63 gene network that includes target genes and regulatory elements may reveal new genes for other malformation disorders. We performed genome-wide DNA-binding profiling by chromatin immunoprecipitation (ChIP), followed by deep sequencing (ChIP-seq) in primary human keratinocytes, and identified potential target genes and regulatory elements controlled by p63. We show that p63 binds to an enhancer element in the SHFM1 locus on chromosome 7q and that this element controls expression of DLX6 and possibly DLX5, both of which are important for limb development. A unique micro-deletion including this enhancer element, but not the DLX5/DLX6 genes, was identified in a patient with SHFM. Our study strongly indicates disruption of a non-coding cis-regulatory element located more than 250 kb from the DLX5/DLX6 genes as a novel disease mechanism in SHFM1. These data provide a proof-of-concept that the catalogue of p63 binding sites identified in this study may be of relevance to the studies of SHFM and other congenital malformations that resemble the p63-associated phenotypes., Competing Interests: The authors have declared that no competing interests exist.