Global Composition of the Bacteriophage Community in Honey Bees

The microbial communities in animal digestive systems are critical for host development and health. They stimulate the immune system during development, synthesize important chemical compounds like hormones, aid in digestion, competitively exclude pathogens, etc. Compared to the bacterial and fungal components of the microbiome, we know little about the temporal and spatial dynamics of bacteriophage communities in animal digestive systems. Recently, the bacteriophages of the honey bee gut were characterized in two European bee populations. Most of the bacteriophages described in these two reports were novel, harbored many metabolic genes in their genomes, and had a community structure that suggests coevolution with their bacterial hosts. To describe the conservation of bacteriophages in bees and begin to understand their role in the bee microbiome, we sequenced the virome of Apis mellifera from Austin, TX, and compared bacteriophage compositions among three locations around the world. We found that most bacteriophages from Austin are novel, sharing no sequence similarity with anything in public repositories. However, many bacteriophages are shared among the three bee viromes, indicating specialization of bacteriophages in the bee gut. Our study, along with the two previous bee virome studies, shows that the bee gut bacteriophage community is simple compared to that of many animals, consisting of several hundred types of bacteriophages that primarily infect four of the dominant bacterial phylotypes in the bee gut. IMPORTANCE Viruses that infect bacteria (bacteriophages) are abundant in the microbial communities that live on and in plants and animals. However, our knowledge of the structure, dynamics, and function of these viral communities lags far behind our knowledge of their bacterial hosts. We sequenced the first bacteriophage community of honey bees from the United States and compared the U.S. honey bee bacteriophage community to those of samples from Europe. Our work is an important characterization of an economically critical insect species and shows how bacteriophage communities can contain highly conserved individuals and be highly variable in composition across a wide geographic range.

Priming With Rhinovirus Protects Mice Against a Lethal Pulmonary Coronavirus Infection

Rhinoviruses (RV) have been shown to inhibit subsequent infection by heterologous respiratory viruses, including influenza viruses and severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2). To better understand the mechanisms whereby RV protects against pulmonary coronavirus infection, we used a native murine virus, mouse hepatitis virus strain 1 (MHV-1), that causes severe disease in the lungs of infected mice. We found that priming of the respiratory tract with RV completely prevented mortality and reduced morbidity of a lethal MHV-1 infection. Replication of MHV-1 was reduced in RV-primed mouse lungs although expression of antiviral type I interferon, IFN-β, was more robust in mice infected with MHV-1 alone. We further showed that signaling through the type I interferon receptor was required for survival of mice given a non-lethal dose of MHV-1. RV-primed mice had reduced pulmonary inflammation and hemorrhage and influx of leukocytes, especially neutrophils, in the airways upon MHV-1 infection. Although MHV-1 replication was reduced in RV-primed mice, RV did not inhibit MHV-1 replication in coinfected lung epithelial cells in vitro. In summary, RV-mediated priming in the respiratory tract reduces viral replication, inflammation, and tissue damage, and prevents mortality of a pulmonary coronavirus infection in mice. These results contribute to our understanding of how distinct respiratory viruses interact with the host to affect disease pathogenesis, which is a critical step in understanding how respiratory viral coinfections impact human health.

Nuku, a family of primate retrocopies derived from KU70

The gene encoding the ubiquitous DNA repair protein, Ku70p, has undergone extensive copy number expansion during primate evolution. Gene duplications of KU70 have the hallmark of long interspersed element-1 mediated retrotransposition with evidence of target-site duplications, the poly-A tails, and the absence of introns. Evolutionary analysis of this expanded family of KU70-derived “NUKU” retrocopies reveals that these genes are both ancient and also actively being created in extant primate species. NUKU retrocopies show evidence of functional divergence away from KU70, as evinced by their altered pattern of tissue expression and possible tissue-specific translation. Molecular modeling predicted that amino acid changes in Nuku2p at the interaction interface with Ku80p would prevent the assembly of the Ku heterodimer. The lack of Nuku2p-Ku80p interaction was confirmed by yeast two-hybrid assay, which contrasts the robust interaction of Ku70p-Ku80p. While several NUKU retrocopies appear to have been degraded by mutation, NUKU2 shows evidence of positive natural selection, suggesting that this retrocopy is undergoing neofunctionalization. Although Nuku proteins do not appear to antagonize retrovirus transduction in cell culture, the observed expansion and rapid evolution of NUKUs could be being driven by alternative selective pressures related to infectious disease or an undefined role in primate physiology.

Plasmodium falciparum DHFR and DHPS Mutations Are Associated With HIV-1 Co-Infection and a Novel DHPS Mutation I504T Is Identified in Western Kenya

Antifolate resistance is significant in Kenya and presumed to result from extensive use and cross-resistance between antifolate antimalarials and antibiotics, including cotrimoxazole/Bactrim used for HIV-1 chemotherapy. However, little is known about antifolate-resistant malaria in the context of newly diagnosed HIV-1 co-infection prior to administration of HIV-1 chemotherapy. Blood samples from a cross-sectional study of asymptomatic adult Kenyans enrolled during voluntary HIV testing were analyzed by PCR for Plasmodium spp. More than 95% of volunteers with identifiable parasite species (132 HIV-1 co-infected) were infected with Plasmodium falciparum alone or P. falciparum with Plasmodium ovale and/or Plasmodium malariae. Deep sequencing was used to screen for mutations in P. falciparum dihydrofolate reductase (dhfr) (N51I, C59R, S108N, I164L) and dihydropteroate synthase (dhps) (S436H, A437G, K540E, A581G) from 1133 volunteers. Individual mutations in DHPS but not DHFR correlated with HIV-1 status. DHFR haplotype diversity was significantly different among volunteers by gender and HIV-1 status. DHPS haplotype diversity by HIV-1 status was significantly different between volunteers paired by age and gender, indicating that patterns of resistance were independent of these variables. Molecular simulations for a novel DHPS mutation (I504T) suggested that the mutated protein has increased affinity for the endogenous ligand DHPPP and decreased affinity for drug binding. A sub-group of monoclonal infections revealed that age and parasitemia were not correlated and enabled identification of a rare septuple-mutant haplotype (IRNL-HGEA). In our study, adult Kenyans newly diagnosed with HIV-1 infection were predominantly infected with moderately resistant P. falciparum, with patterns of infecting parasite genotypes significantly associated with HIV-1 status. Together with the discovery of DHPS I504T, these data indicate that antifolate resistance continues to evolve in Kenya. Further, they highlight the need to understand the effects of associated mutations on both fitness and resistance of P. falciparum in the context of HIV-1 co-infection to better inform treatment for asymptomatic malaria.

Organic Farming Sharpens Plant Defenses in the Field

Plants deploy a variety of chemical and physical defenses to protect themselves against herbivores and pathogens. Organic farming seeks to enhance these responses by improving soil quality, ultimately altering bottom up regulation of plant defenses. While laboratory studies suggest this approach is effective, it remains unclear whether organic agriculture encourages more-active plant defenses under real-world conditions. Working on the farms of cooperating growers, we examined gene expression in the leaves of two potato (Solanum tuberosum) varieties, grown on organic vs. conventional farms. For one variety, Norkotah, we found significantly heightened initiation of genes associated with plant-defense pathways in plants grown in organic vs. conventional fields. Organic Norkotah fields exhibited lower levels of nitrate in soil and of nitrogen in plant foliage, alongside differences in communities of soil bacteria, suggesting possible links between soil management and observed differences in plant defenses. Additionally, numbers of predatory and phloem-feeding insects were higher in organic than conventional fields. A second potato variety, Alturas, which is generally grown using fewer inputs and in poorer-quality soils, exhibited lower overall herbivore and predator numbers, few differences in soil ecology, and no differences in gene-activity in organic and conventional farming systems. Altogether, our results suggest that organic farming has the potential to increase plants' resistance to herbivores, possibly facilitating reduced need for insecticide applications. These benefits appear to be mediated by plant variety and/or farming context.

ΦX174 Attenuation by Whole-Genome Codon Deoptimization

Natural selection acting on synonymous mutations in protein-coding genes influences genome composition and evolution. In viruses, introducing synonymous mutations in genes encoding structural proteins can drastically reduce viral growth, providing a means to generate potent, live-attenuated vaccine candidates. However, an improved understanding of what compositional features are under selection and how combinations of synonymous mutations affect viral growth is needed to predictably attenuate viruses and make them resistant to reversion. We systematically recoded all nonoverlapping genes of the bacteriophage ΦX174 with codons rarely used in its Escherichia coli host. The fitness of recombinant viruses decreases as additional deoptimizing mutations are made to the genome, although not always linearly, and not consistently across genes. Combining deoptimizing mutations may reduce viral fitness more or less than expected from the effect size of the constituent mutations and we point out difficulties in untangling correlated compositional features. We test our model by optimizing the same genes and find that the relationship between codon usage and fitness does not hold for optimization, suggesting that wild-type ΦX174 is at a fitness optimum. This work highlights the need to better understand how selection acts on patterns of synonymous codon usage across the genome and provides a convenient system to investigate the genetic determinants of virulence.

One Hundred Mitochondrial Genomes of Cicadas

Mitochondrial genomes can provide valuable information on the biology and evolutionary histories of their host organisms. Here, we present and characterize the complete coding regions of 107 mitochondrial genomes (mitogenomes) of cicadas (Insecta: Hemiptera: Auchenorrhyncha: Cicadoidea), representing 31 genera, 61 species, and 83 populations. We show that all cicada mitogenomes retain the organization and gene contents thought to be ancestral in insects, with some variability among cicada clades in the length of a region between the genes nad2 and cox1, which encodes 3 tRNAs. Phylogenetic analyses using these mitogenomes recapitulate a recent 5-gene classification of cicadas into families and subfamilies, but also identify a species that falls outside of the established taxonomic framework. While protein-coding genes are under strong purifying selection, tests of relative evolutionary rates reveal significant variation in evolutionary rates across taxa, highlighting the dynamic nature of mitochondrial genome evolution in cicadas. These data will serve as a useful reference for future research into the systematics, ecology, and evolution of the superfamily Cicadoidea.

Saturation Mutagenesis Genome Engineering of Infective ΦX174 Bacteriophage via Unamplified Oligo Pools and Golden Gate Assembly

Here we present a novel protocol for the construction of saturation single-site-and massive multisite-mutant libraries of a bacteriophage. We segmented the ΦX174 genome into 14 nontoxic and nonreplicative fragments compatible with Golden Gate assembly. We next used nicking mutagenesis with oligonucleotides prepared from unamplified oligo pools with individual segments as templates to prepare near-comprehensive single-site mutagenesis libraries of genes encoding the F capsid protein (421 amino acids scanned) and G spike protein (172 amino acids scanned). Libraries possessed greater than 99% of all 11 860 programmed mutations. Golden Gate cloning was then used to assemble the complete ΦX174 mutant genome and generate libraries of infective viruses. This protocol will enable reverse genetics experiments for studying viral evolution and, with some modifications, can be applied for engineering therapeutically relevant bacteriophages with larger genomes.

Organic Soils Control Beetle Survival While Competitors Limit Aphid Population Growth

Soil chemistry and microbial diversity can impact the vigor and nutritive qualities of plants, as well as plants' ability to deploy anti-herbivore defenses. Soil qualities often vary dramatically on organic versus conventional farms, reflecting the many differences in soil management practices between these farming systems. We examined soil-mediated effects on herbivore performance by growing potato plants (Solanum tuberosum L.) in soils collected from organic or conventional commercial farm fields, and then exposing these plants to herbivory by green peach aphids (Myzus persicae Sulzer, Hemiptera: Aphididae) and/or Colorado potato beetles (Leptinotarsa decemlineata Say, Coleoptera: Chrysomelidae). Responses of the two potato pests varied dramatically. Survivorship of Colorado potato beetles was almost 3× higher on plants grown in organic than in conventional soils, but was unaffected by the presence of aphids. In contrast, aphid colony growth was twice as rapid when aphids were reared alone rather than with Colorado potato beetles, but was unaffected by soil type. We saw no obvious differences in soil nutrients when comparing organic and conventional soils. However, we saw a higher diversity of bacteria in organic soils, and potato plants grown in this soil had a lower carbon concentration in foliar tissue. In summary, the herbivore species differed in their susceptibility to soil- versus competitor-mediated effects, and these differences may be driven by microbe-mediated changes in host plant quality. Our results suggest that soil-mediated effects on pest growth can depend on herbivore species and community composition, and that soil management strategies that promote plant health may also increase host quality for pests.

Comparative genomics of Bifidobacterium species isolated from marmosets and humans

The genus Bifidobacterium is purported to have beneficial consequences for human health and is a major component of many gastrointestinal probiotics. Although species of Bifidobacterium are generally at low relative frequency in the adult human gastrointestinal tract, they can constitute high proportions of the gastrointestinal communities of adult marmosets. To identify genes that might be important for the maintenance of Bifidobacterium in adult marmosets, ten strains of Bifidobacterium were isolated from the feces of seven adult marmosets, and their genomes were sequenced. There were six B. reuteri strains, two B. callitrichos strains, one B. myosotis sp. nov. and one B. tissieri sp. nov. among our isolates. Phylogenetic analysis showed that three of the four species we isolated were most closely related to B. bifidum, B. breve and B. longum, which are species found in high abundance in human infants. There were 1357 genes that were shared by at least one strain of B. reuteri, B. callitrichos, B. breve, and B. longum, and 987 genes that were found in all strains of the four species. There were 106 genes found in B. reuteri and B. callitrichos but not in human bifidobacteria, and several of these genes were involved in nutrient uptake. These pathways for nutrient uptake appeared to be specific to Bifidobacterium from New World monkeys. Additionally, the distribution of Bifidobacterium in fecal samples from captive adult marmosets constituted as much as 80% of the gut microbiome, although this was variable between individuals and colonies. We suggest that nutrient transporters may be important for the maintenance of Bifidobacterium during adulthood in marmosets.

A Rapid Method for Sequencing Double-Stranded RNAs Purified from Yeasts and the Identification of a Potent K1 Killer Toxin Isolated from Saccharomyces cerevisiae

Mycoviruses infect a large number of diverse fungal species, but considering their prevalence, relatively few high-quality genome sequences have been determined. Many mycoviruses have linear double-stranded RNA genomes, which makes it technically challenging to ascertain their nucleotide sequence using conventional sequencing methods. Different specialist methodologies have been developed for the extraction of double-stranded RNAs from fungi and the subsequent synthesis of cDNAs for cloning and sequencing. However, these methods are often labor-intensive, time-consuming, and can require several days to produce cDNAs from double-stranded RNAs. Here, we describe a comprehensive method for the rapid extraction and sequencing of dsRNAs derived from yeasts, using short-read next generation sequencing. This method optimizes the extraction of high-quality double-stranded RNAs from yeasts and 3′ polyadenylation for the initiation of cDNA synthesis for next-generation sequencing. We have used this method to determine the sequence of two mycoviruses and a double-stranded RNA satellite present within a single strain of the model yeast Saccharomyces cerevisiae. The quality and depth of coverage was sufficient to detect fixed and polymorphic mutations within viral populations extracted from a clonal yeast population. This method was also able to identify two fixed mutations within the alpha-domain of a variant K1 killer toxin encoded on a satellite double-stranded RNA. Relative to the canonical K1 toxin, these newly reported mutations increased the cytotoxicity of the K1 toxin against a specific species of yeast.

Multiple origins of interdependent endosymbiotic complexes in a genus of cicadas

Bacterial endosymbionts that provide nutrients to hosts often have genomes that are extremely stable in structure and gene content. In contrast, the genome of the endosymbiont Hodgkinia cicadicola has fractured into multiple distinct lineages in some species of the cicada genus Tettigades To better understand the frequency, timing, and outcomes of Hodgkinia lineage splitting throughout this cicada genus, we sampled cicadas over three field seasons in Chile and performed genomics and microscopy on representative samples. We found that a single ancestral Hodgkinia lineage has split at least six independent times in Tettigades over the last 4 million years, resulting in complexes of between two and six distinct Hodgkinia lineages per host. Individual genomes in these symbiotic complexes differ dramatically in relative abundance, genome size, organization, and gene content. Each Hodgkinia lineage retains a small set of core genes involved in genetic information processing, but the high level of gene loss experienced by all genomes suggests that extensive sharing of gene products among symbiont cells must occur. In total, Hodgkinia complexes that consist of multiple lineages encode nearly complete sets of genes present on the ancestral single lineage and presumably perform the same functions as symbionts that have not undergone splitting. However, differences in the timing of the splits, along with dissimilar gene loss patterns on the resulting genomes, have led to very different outcomes of lineage splitting in extant cicadas.

Selecting among three basic fitness landscape models: Additive, multiplicative and stickbreaking

Fitness landscapes map genotypes to organismal fitness. Their topographies depend on how mutational effects interact - epistasis - andare important for understanding evolutionary processes such as speciation, the rate of adaptation, the advantage of recombination, and the predictability versus stochasticity of evolution. The growing amount of data has made it possible to better test landscape models empirically. We argue that this endeavor will benefit from the development and use of meaningful basic models against which to compare more complex models. Here we develop statistical and computational methods for fitting fitness data from mutation combinatorial networks to three simple models: additive, multiplicative and stickbreaking. We employ a Bayesian framework for doing model selection. Using simulations, we demonstrate that our methods work and we explore their statistical performance: bias, error, and the power to discriminate among models. We then illustrate our approach and its flexibility by analyzing several previously published datasets. An R-package that implements our methods is available in the CRAN repository under the name Stickbreaker.

Biochemical association of metabolic profile and microbiome in chronic pressure ulcer wounds

Chronic, non-healing wounds contribute significantly to the suffering of patients with co-morbidities in the clinical population with mild to severely compromised immune systems. Normal wound healing proceeds through a well-described process. However, in chronic wounds this process seems to become dysregulated at the transition between resolution of inflammation and re-epithelialization. Bioburden in the form of colonizing bacteria is a major contributor to the delayed headlining in chronic wounds such as pressure ulcers. However how the microbiome influences the wound metabolic landscape is unknown. Here, we have used a Systems Biology approach to determine the biochemical associations between the taxonomic and metabolomic profiles of wounds colonized by bacteria. Pressure ulcer biopsies were harvested from primary chronic wounds and bisected into top and bottom sections prior to analysis of microbiome by pyrosequencing and analysis of metabolome using 1H nuclear magnetic resonance (NMR) spectroscopy. Bacterial taxonomy revealed that wounds were colonized predominantly by three main phyla, but differed significantly at the genus level. While taxonomic profiles demonstrated significant variability between wounds, metabolic profiles shared significant similarity based on the depth of the wound biopsy. Biochemical association between taxonomy and metabolic landscape indicated significant wound-to-wound similarity in metabolite enrichment sets and metabolic pathway impacts, especially with regard to amino acid metabolism. To our knowledge, this is the first demonstration of a statistically robust correlation between bacterial colonization and metabolic landscape within the chronic wound environment.

Sympatric speciation in a bacterial endosymbiont results in two genomes with the functionality of one

Mutualisms that become evolutionarily stable give rise to organismal interdependencies. Some insects have developed intracellular associations with communities of bacteria, where the interdependencies are manifest in patterns of complementary gene loss and retention among members of the symbiosis. Here, using comparative genomics and microscopy, we show that a three-member symbiotic community has become a four-way assemblage through a novel bacterial lineage-splitting event. In some but not all cicada species of the genus Tettigades, the endosymbiont Candidatus Hodgkinia cicadicola has split into two new cytologically distinct but metabolically interdependent species. Although these new bacterial genomes are partitioned into discrete cell types, the intergenome patterns of gene loss and retention are almost perfectly complementary. These results defy easy classification: they show genomic patterns consistent with those observed after both speciation and whole-genome duplication. We suggest that our results highlight the potential power of nonadaptive forces in shaping organismal complexity.

Dynamic recruitment of amino acid transporters to the insect/symbiont interface

Symbiosis is well known to influence bacterial symbiont genome evolution and has recently been shown to shape eukaryotic host genomes. Intriguing patterns of host genome evolution, including remarkable numbers of gene duplications, have been observed in the pea aphid, a sap-feeding insect that relies on a bacterial endosymbiont for amino acid provisioning. Previously, we proposed that gene duplication has been important for the evolution of symbiosis based on aphid-specific gene duplication in amino acid transporters (AATs), with some paralogs highly expressed in the cells housing symbionts (bacteriocytes). Here, we use a comparative approach to test the role of gene duplication in enabling recruitment of AATs to bacteriocytes. Using genomic and transcriptomic data, we annotate AATs from sap-feeding and non sap-feeding insects and find that, like aphids, AAT gene families have undergone independent large-scale gene duplications in three of four additional sap-feeding insects. RNA-seq differential expression data indicate that, like aphids, the sap-feeding citrus mealybug possesses several lineage-specific bacteriocyte-enriched paralogs. Further, differential expression data combined with quantitative PCR support independent evolution of bacteriocyte enrichment in sap-feeding insect AATs. Although these data indicate that gene duplication is not necessary to initiate host/symbiont amino acid exchange, they support a role for gene duplication in enabling AATs to mediate novel host/symbiont interactions broadly in the sap-feeding suborder Sternorrhyncha. In combination with recent studies on other symbiotic systems, gene duplication is emerging as a general pattern in host genome evolution.

An AT mutational bias in the tiny GC-rich endosymbiont genome of Hodgkinia

The fractional guanine + cytosine (GC) contents of sequenced bacterial genomes range from 13% to 75%. Despite several decades of research aimed at understanding this wide variation, the forces controlling GC content are not well understood. Recent work has suggested that a universal adenine + thymine (AT) mutational bias exists in all bacteria and that the elevated GC contents found in some bacterial genomes is due to genome-wide selection for increased GC content. These results are generally consistent with the low GC contents observed in most strict endosymbiotic bacterial genomes, where the loss of DNA repair mechanisms combined with the population genetic effects of small effective population sizes and decreased recombination should lower the efficacy of selection and shift the equilibrium GC content in the mutationally favored AT direction. Surprisingly, the two smallest bacterial genomes, Candidatus Hodgkinia cicadicola (144 kb) and Candidatus Tremblaya princeps (139 kb), have the unusual combination of highly reduced genomes and elevated GC contents, raising the possibility that these bacteria may be exceptions to the otherwise apparent universal bacterial AT mutational bias. Here, using population genomic data generated from the Hodgkinia genome project, we show that Hodgkinia has a clear AT mutational bias. These results provide further evidence that an AT mutational bias is universal in bacteria, even in strict endosymbionts with elevated genomic GC contents.