The persistence and stabilization of auxiliary genes in the human skin virome.

BACKGROUND: The human skin contains a diverse microbiome that provides protective functions against environmental pathogens. Studies have demonstrated that bacteriophages modulate bacterial community composition and facilitate the transfer of host-specific genes, potentially influencing host cellular functions. However, little is known about the human skin virome and its role in human health. Especially, how viral-host relationships influence skin microbiome structure and function is poorly understood. RESULTS: Population dynamics and genetic diversity of bacteriophage communities in viral metagenomic data collected from three anatomical skin locations from 60 subjects at five different time points revealed that cutaneous bacteriophage populations are mainly composed of tailed Caudovirales phages that carry auxiliary genes to help improve metabolic remodeling to increase bacterial host fitness through antimicrobial resistance. Sequence variation in the MRSA associated antimicrobial resistance gene, erm(C) was evaluated using targeted sequencing to further confirm the presence of antimicrobial resistance genes in the human virome and to demonstrate how functionality of such genes may influence persistence and in turn stabilization of bacterial host and their functions. CONCLUSIONS: This large temporal study of human skin associated viruses indicates that the human skin virome is associated with auxiliary metabolic genes and antimicrobial resistance genes to help increase bacterial host fitness.

Genome Sequence of Feline Papillomavirus Strain P20 Assembled from Metagenomic Data from the Skin of a House Cat Owner.

A feline papillomavirus genome was assembled from metagenomic sequencing data collected from the skin of a house cat owner. The circular genome of strain P20 is 8,069 bp in length, has a GC content of 54.38%, and displays genome organization typical of feline papillomaviruses. The genome exhibits approximately 75% sequence similarity to other feline papillomavirus genomes.

Association analyses of host genetics, root-colonizing microbes, and plant phenotypes under different nitrogen conditions in maize.

The root-associated microbiome (rhizobiome) affects plant health, stress tolerance, and nutrient use efficiency. However, it remains unclear to what extent the composition of the rhizobiome is governed by intraspecific variation in host plant genetics in the field and the degree to which host plant selection can reshape the composition of the rhizobiome. Here, we quantify the rhizosphere microbial communities associated with a replicated diversity panel of 230 maize (Zea mays L.) genotypes grown in agronomically relevant conditions under high N (+N) and low N (-N) treatments. We analyze the maize rhizobiome in terms of 150 abundant and consistently reproducible microbial groups and we show that the abundance of many root-associated microbes is explainable by natural genetic variation in the host plant, with a greater proportion of microbial variance attributable to plant genetic variation in -N conditions. Population genetic approaches identify signatures of purifying selection in the maize genome associated with the abundance of several groups of microbes in the maize rhizobiome. Genome-wide association study was conducted using the abundance of microbial groups as rhizobiome traits, and n=622 plant loci were identified that are linked to the abundance of n=104 microbial groups in the maize rhizosphere. In 62/104 cases, which is more than expected by chance, the abundance of these same microbial groups was correlated with variation in plant vigor indicators derived from high throughput phenotyping of the same field experiment. We provide comprehensive datasets about the three-way interaction of host genetics, microbe abundance, and plant performance under two N treatments to facilitate targeted experiments toward harnessing the full potential of root-associated microbial symbionts in maize production.

Improving Taxonomic Delimitation of Fungal Species in the Age of Genomics and Phenomics.

Species concepts have long provided a source of debate among biologists. These lively debates have been important for reaching consensus on how to communicate across scientific disciplines and for advancing innovative strategies to study evolution, population biology, ecology, natural history, and disease epidemiology. Species concepts are also important for evaluating variability and diversity among communities, understanding biogeographical distributions, and identifying causal agents of disease across animal and plant hosts. While there have been many attempts to address the concept of species in the fungi, there are several concepts that have made taxonomic delimitation especially challenging. In this review we discuss these major challenges and describe methodological approaches that show promise for resolving ambiguity in fungal taxonomy by improving discrimination of genetic and functional traits. We highlight the relevance of eco-evolutionary theory used in conjunction with integrative taxonomy approaches to improve the understanding of interactions between environment, ecology, and evolution that give rise to distinct species boundaries. Beyond recent advances in genomic and phenomic methods, bioinformatics tools and modeling approaches enable researchers to test hypothesis and expand our knowledge of fungal biodiversity. Looking to the future, the pairing of integrative taxonomy approaches with multi-locus genomic sequencing and phenomic techniques, such as transcriptomics and proteomics, holds great potential to resolve many unknowns in fungal taxonomic classification.

The application of the skin virome for human identification.

The use of skin virome offers a unique approach for human identification purposes in instances where a viable and statistically relevant human DNA profile is unavailable. The skin virome may act as an alternative DNA profile and/or an additional form of probative genetic material. To date, no study has attempted to investigate the human virome over a time series across various physical locations of the body to identify its diagnostic potential as a tool for human identification. For this study, we set out to evaluate the stability, diversity, and individualization of the human skin virome. An additional goal was to identify putative viral signatures that can be used in conjunction with traditional forensic STR loci. In order to accomplish this, human viral metagenomes were collected and sequenced from 42 individuals at three anatomical locations (left hand, right hand, and scalp) across multiple collection periods over a 6-month window of time. Assembly dependent and independent bioinformatic approaches, along with a database centered assessment of viral identification, resulted in three sets of stable putative viral markers. In total, with the three sets combined, we identified 59 viral biomarker regions, consisting of viral species and uncharacterized viral genome assemblies, that were stable over the sampling period. Additionally, we found the abundance profiles of these 59 viral biomarkers, based on presence or absence, to be significantly different across subjects (P < 0.001). Here we demonstrate that not only is the human virome applicable to be used for human identification, but we have identified many viral signatures that can putatively be used for forensic applications, thus providing a foundation to the novel field of forensic virology.

Global Distribution and Richness of Armillaria and Related Species Inferred From Public Databases and Amplicon Sequencing Datasets.

Armillaria is a globally distributed fungal genus most notably composed of economically important plant pathogens that are found predominantly in forest and agronomic systems. The genus sensu lato has more recently received attention for its role in woody plant decomposition and in mycorrhizal symbiosis with specific plants. Previous phylogenetic analyses suggest that around 50 species are recognized globally. Despite this previous work, no studies have analyzed the global species richness and distribution of the genus using data derived from fungal community sequencing datasets or barcoding initiatives. To assess the global diversity and species richness of Armillaria, we mined publicly available sequencing datasets derived from numerous primer regions for the ribosomal operon, as well as ITS sequences deposited on Genbank, and clustered them akin to metabarcoding studies. Our estimates reveal that species richness ranges from 50 to 60 species, depending on whether the ITS1 or ITS2 marker is used. Eastern Asia represents the biogeographic region with the highest species richness. We also assess the overlap of species across geographic regions and propose some hypotheses regarding the drivers of variability in species diversity and richness between different biogeographic regions.

Transcriptomics Reveals the Putative Mycoparasitic Strategy of the Mushroom Entoloma abortivum on Species of the Mushroom Genus Armillaria.

During mycoparasitism, a fungus-the host-is parasitized by another fungus-the mycoparasite. The genetic underpinnings of these relationships have been best characterized in ascomycete fungi. However, within basidiomycete fungi, there are rare instances of mushroom-forming species parasitizing the reproductive structures, or sporocarps, of other mushroom-forming species, which have been rarely investigated on a genetic level. One of the most enigmatic of these occurs between Entoloma abortivum and species of Armillaria, where hyphae of E. abortivum are hypothesized to disrupt the development of Armillaria sporocarps, resulting in the formation of carpophoroids. However, it remains unknown whether carpophoroids are the direct result of a mycoparasitic relationship. To address the nature of this unique interaction, we analyzed gene expression of field-collected Armillaria and E. abortivum sporocarps and carpophoroids. Transcripts in the carpophoroids are primarily from E. abortivum, supporting the hypothesis that this species is parasitizing Armillaria. Most notably, we identified differentially upregulated E. abortivum ß-trefoil-type lectins in the carpophoroid, which we hypothesize bind to Armillaria cell wall galactomannoproteins, thereby mediating recognition between the mycoparasite and the host. The most differentially upregulated E. abortivum transcripts in the carpophoroid code for oxalate decarboxylases-enzymes that degrade oxalic acid. Oxalic acid is a virulence factor in many plant pathogens, including Armillaria species; however, E. abortivum has evolved a sophisticated strategy to overcome this defense mechanism. The number of gene models and genes that code for carbohydrate-active enzymes in the E. abortivum transcriptome was reduced compared to other closely related species, perhaps as a result of the specialized nature of this interaction. IMPORTANCE By studying fungi that parasitize other fungi, we can understand the basic biology of these unique interactions. Studies focused on the genetic mechanisms regulating mycoparasitism between host and parasite have thus far concentrated on a single fungal lineage within the Ascomycota. The work presented here expands our understanding of mycoparasitic relationships to the Basidiomycota and represents the first transcriptomic study to our knowledge that examines fungal-fungal relationships in their natural setting. The results presented here suggest that even distantly related mycoparasites utilize similar mechanisms to parasitize their host. Given that species of the mushroom-forming pathogen Armillaria cause plant root-rot diseases in many agroecosystems, an enhanced understanding of this interaction may contribute to better control of these diseases through biocontrol applications.

Symbiotic nitrogen fixation in the reproductive structures of a basidiomycete fungus.

Nitrogen (N) fixation is a driving force for the formation of symbiotic associations between N2-fixing bacteria and eukaryotes.1 Limited examples of these associations are known in fungi, and none with sexual structures of non-lichenized species.2-6 The basidiomycete Guyanagaster necrorhizus is a sequestrate fungus endemic to the Guiana Shield.7 Like the root rot-causing species in its sister genera Armillaria and Desarmillaria, G. necrorhizus sporocarps fruit from roots of decaying trees (Figures 1A-1C),8 and genome sequencing is consistent with observations that G. necrorhizus is a white-rotting decomposer. This species also represents the first documentation of an arthropod-dispersed sequestrate fungus. Numerous species of distantly related wood-feeding termites, which scavenge for N-rich food, feed on the mature spore-bearing tissue, or gleba, of G. necrorhizus. During feeding, mature spores adhere to termites for subsequent dispersal.9 Using chemical assays, isotope analysis, and high-throughput sequencing, we show that the sporocarps harbor actively N2-fixing Enterobacteriaceae species and that the N content within fungal tissue increases with maturation. Untargeted proteomic profiling suggests that ATP generation in the gleba is accomplished via fermentation. The use of fermentation-an anaerobic process-indicates that the sporocarp environment is anoxic, likely an adaptation to protect the oxygen-sensitive nitrogenase enzyme. Sporocarps also have a thick outer covering, possibly to limit oxygen diffusion. The enriched N content within mature sporocarps may offer a dietary inducement for termites in exchange for spore dispersal. These results show that the flexible metabolic capacity of fungi may facilitate N2-fixing associations, as well as higher-level organismal associations.

Rhizosphere Microbiomes in a Historical Maize-Soybean Rotation System Respond to Host Species and Nitrogen Fertilization at the Genus and Subgenus Levels.

Root-associated microbes are key players in plant health, disease resistance, and nitrogen (N) use efficiency. It remains largely unclear how the interplay of biological and environmental factors affects rhizobiome dynamics in agricultural systems. In this study, we quantified the composition of rhizosphere and bulk soil microbial communities associated with maize (Zea mays L.) and soybean (Glycine max L.) in a long-term crop rotation study under conventional fertilization and low-N regimes. Over two growing seasons, we evaluated the effects of environmental conditions and several treatment factors on the abundance of rhizosphere- and soil-colonizing microbial taxa. Time of sampling, host plant species, and N fertilization had major effects on microbiomes, while no effect of crop rotation was observed. Using variance partitioning as well as 16S sequence information, we further defined a set of 82 microbial genera and functional taxonomic groups at the subgenus level that show distinct responses to treatment factors. We identified taxa that are highly specific to either maize or soybean rhizospheres, as well as taxa that are sensitive to N fertilization in plant rhizospheres and bulk soil. This study provides insights to harness the full potential of soil microbes in maize and soybean agricultural systems through plant breeding and field management. IMPORTANCE Plant roots are colonized by large numbers of microbes, some of which may help the plant acquire nutrients and fight diseases. Our study contributes to a better understanding of root-colonizing microbes in the widespread and economically important maize-soybean crop rotation system. The long-term goal of this research is to optimize crop plant varieties and field management to create the best possible conditions for beneficial plant-microbe interactions to occur. These beneficial microbes may be harnessed to sustainably reduce dependency on pesticides and industrial fertilizer. We identify groups of microbes specific to the maize or to the soybean host and microbes that are sensitive to nitrogen fertilization. These microbes represent candidates that may be influenced through plant breeding or field management, and future research will be directed toward elucidating their roles in plant health and nitrogen usage.

Complete Genome Sequence of Geobacter sp. Strain FeAm09, a Moderately Acidophilic Soil Bacterium.

A moderately acidophilic Geobacter sp. strain, FeAm09, was isolated from forest soil. The complete genome sequence is 4,099,068 bp with an average GC content of 61.1%. No plasmids were detected. The genome contains a total of 3,843 genes and 3,608 protein-coding genes, including genes supporting iron and nitrogen biogeochemical cycling.

Sequence-based classification and identification of Fungi.

Fungal taxonomy and ecology have been revolutionized by the application of molecular methods and both have increasing connections to genomics and functional biology. However, data streams from traditional specimen- and culture-based systematics are not yet fully integrated with those from metagenomic and metatranscriptomic studies, which limits understanding of the taxonomic diversity and metabolic properties of fungal communities. This article reviews current resources, needs, and opportunities for sequence-based classification and identification (SBCI) in fungi as well as related efforts in prokaryotes. To realize the full potential of fungal SBCI it will be necessary to make advances in multiple areas. Improvements in sequencing methods, including long-read and single-cell technologies, will empower fungal molecular ecologists to look beyond ITS and current shotgun metagenomics approaches. Data quality and accessibility will be enhanced by attention to data and metadata standards and rigorous enforcement of policies for deposition of data and workflows. Taxonomic communities will need to develop best practices for molecular characterization in their focal clades, while also contributing to globally useful datasets including ITS. Changes to nomenclatural rules are needed to enable validPUBLICation of sequence-based taxon descriptions. Finally, cultural shifts are necessary to promote adoption of SBCI and to accord professional credit to individuals who contribute to community resources.

Genome Sequence of Fusarium Isolate MYA-4552 from the Midgut of Anoplophora glabripennis, an Invasive, Wood-Boring Beetle.

The Fusarium solani species complex (FSSC) is a clade of environmentally ubiquitous fungi that includes plant, animal, and insect associates. Here, we report the draft genome sequence of the undescribed species FSSC 6 (isolate MYA-4552), housed in the gut of the wood-boring cerambycid beetle Anoplophora glabripennis.

The khmer software package: enabling efficient nucleotide sequence analysis.

The khmer package is a freely available software library for working efficiently with fixed length DNA words, or k-mers. khmer provides implementations of a probabilistic k-mer counting data structure, a compressible De Bruijn graph representation, De Bruijn graph partitioning, and digital normalization. khmer is implemented in C++ and Python, and is freely available under the BSD license at  https://github.com/dib-lab/khmer/.

Genome of an arbuscular mycorrhizal fungus provides insight into the oldest plant symbiosis.

The mutualistic symbiosis involving Glomeromycota, a distinctive phylum of early diverging Fungi, is widely hypothesized to have promoted the evolution of land plants during the middle Paleozoic. These arbuscular mycorrhizal fungi (AMF) perform vital functions in the phosphorus cycle that are fundamental to sustainable crop plant productivity. The unusual biological features of AMF have long fascinated evolutionary biologists. The coenocytic hyphae host a community of hundreds of nuclei and reproduce clonally through large multinucleated spores. It has been suggested that the AMF maintain a stable assemblage of several different genomes during the life cycle, but this genomic organization has been questioned. Here we introduce the 153-Mb haploid genome of Rhizophagus irregularis and its repertoire of 28,232 genes. The observed low level of genome polymorphism (0.43 SNP per kb) is not consistent with the occurrence of multiple, highly diverged genomes. The expansion of mating-related genes suggests the existence of cryptic sex-related processes. A comparison of gene categories confirms that R. irregularis is close to the Mucoromycotina. The AMF obligate biotrophy is not explained by genome erosion or any related loss of metabolic complexity in central metabolism, but is marked by a lack of genes encoding plant cell wall-degrading enzymes and of genes involved in toxin and thiamine synthesis. A battery of mycorrhiza-induced secreted proteins is expressed in symbiotic tissues. The present comprehensive repertoire of R. irregularis genes provides a basis for future research on symbiosis-related mechanisms in Glomeromycota.

Metagenomic profiling reveals lignocellulose degrading system in a microbial community associated with a wood-feeding beetle.

The Asian longhorned beetle (Anoplophoraglabripennis) is an invasive, wood-boring pest that thrives in the heartwood of deciduous tree species. A large impediment faced by A. glabripennis as it feeds on woody tissue is lignin, a highly recalcitrant biopolymer that reduces access to sugars and other nutrients locked in cellulose and hemicellulose. We previously demonstrated that lignin, cellulose, and hemicellulose are actively deconstructed in the beetle gut and that the gut harbors an assemblage of microbes hypothesized to make significant contributions to these processes. While lignin degrading mechanisms have been well characterized in pure cultures of white rot basidiomycetes, little is known about such processes in microbial communities associated with wood-feeding insects. The goals of this study were to develop a taxonomic and functional profile of a gut community derived from an invasive population of larval A. glabripennis collected from infested host trees and to identify genes that could be relevant for the digestion of woody tissue and nutrient acquisition. To accomplish this goal, we taxonomically and functionally characterized the A. glabripennis midgut microbiota through amplicon and shotgun metagenome sequencing and conducted a large-scale comparison with the metagenomes from a variety of other herbivore-associated communities. This analysis distinguished the A. glabripennis larval gut metagenome from the gut communities of other herbivores, including previously sequenced termite hindgut metagenomes. Genes encoding enzymes were identified in the A. glabripennis gut metagenome that could have key roles in woody tissue digestion including candidate lignin degrading genes (laccases, dye-decolorizing peroxidases, novel peroxidases and ß-etherases), 36 families of glycoside hydrolases (such as cellulases and xylanases), and genes that could facilitate nutrient recovery, essential nutrient synthesis, and detoxification. This community could serve as a reservoir of novel enzymes to enhance industrial cellulosic biofuels production or targets for novel control methods for this invasive and highly destructive insect.