Identification of rare microbial colonizers of plastic materials incubated in a coral reef environment.

Plastic waste accumulation in marine environments has complex, unintended impacts on ecology that cross levels of community organization. To measure succession in polyolefin-colonizing marine bacterial communities, an in situ time-series experiment was conducted in the oligotrophic coastal waters of the Bermuda Platform. Our goals were to identify polyolefin colonizing taxa and isolate bacterial cultures for future studies of the biochemistry of microbe-plastic interactions. HDPE, LDPE, PP, and glass coupons were incubated in surface seawater for 11 weeks and sampled at two-week intervals. 16S rDNA sequencing and ATR-FTIR/HIM were used to assess biofilm community structure and chemical changes in polymer surfaces. The dominant colonizing taxa were previously reported cosmopolitan colonizers of surfaces in marine environments, which were highly similar among the different plastic types. However, significant differences in rare community composition were observed between plastic types, potentially indicating specific interactions based on surface chemistry. Unexpectedly, a major transition in community composition occurred in all material treatments between days 42 and 56 (p < 0.01). Before the transition, Alteromonadaceae, Marinomonadaceae, Saccharospirillaceae, Vibrionaceae, Thalassospiraceae, and Flavobacteriaceae were the dominant colonizers. Following the transition, the relative abundance of these taxa declined, while Hyphomonadaceae, Rhodobacteraceae and Saprospiraceae increased. Over the course of the incubation, 8,641 colonizing taxa were observed, of which 25 were significantly enriched on specific polyolefins. Seven enriched taxa from families known to include hydrocarbon degraders (Hyphomonadaceae, Parvularculaceae and Rhodobacteraceae) and one n-alkane degrader (Ketobacter sp.). The ASVs that exhibited associations with specific polyolefins are targets of ongoing investigations aimed at retrieving plastic-degrading microbes in culture.

Sulforaphane Bioavailability in Healthy Subjects Fed a Single Serving of Fresh Broccoli Microgreens.

Cruciferous vegetable consumption is associated with numerous health benefits attributed to the phytochemical sulforaphane (SFN) that exerts antioxidant and chemopreventive properties, among other bioactive compounds. Broccoli sprouts, rich in SFN precursor glucoraphanin (GRN), have been investigated in numerous clinical trials. Broccoli microgreens are similarly rich in GRN but have remained largely unexplored. The goal of this study was to examine SFN bioavailability and the microbiome profile in subjects fed a single serving of fresh broccoli microgreens. Eleven subjects participated in a broccoli microgreens feeding study. Broccoli microgreens GRN and SFN contents and stability were measured. Urine and stool SFN metabolite profiles and microbiome composition were examined. Broccoli microgreens had similar GRN content to values previously reported for broccoli sprouts, which was stable over time. Urine SFN metabolite profiles in broccoli microgreens-fed subjects were similar to those reported previously in broccoli sprouts-fed subjects, including the detection of SFN-nitriles. We also reported the detection of SFN metabolites in stool samples for the first time. A single serving of broccoli microgreens did not significantly alter microbiome composition. We showed in this study that broccoli microgreens are a significant source of SFN. Our work provides the foundation for future studies to establish the health benefits of broccoli microgreens consumption.

Woronichinia naegeliana: A common nontoxigenic component of temperate freshwater cyanobacterial blooms with 30% of its genome in transposons.

Monitoring in the U.S. state of Washington across the period 2007-2019 showed that Woronichinia has been present in many lakes state-wide. This cyanobacterium was commonly dominant or sub-dominant in cyanobacterial blooms in the wet temperate region west of the Cascade Mountains. In these lakes, Woronichinia often co-existed with Microcystis, Dolichospermum and Aphanizomenon flos-aquae and the cyanotoxin microcystin has often been present in those blooms, although it has not been known whether Woronichinia is a toxin producer. We report the first complete genome of Woronichinia naegeliana WA131, assembled from the metagenome of a sample collected from Wiser Lake, Washington, in 2018. The genome contains no genes for cyanotoxin biosynthesis or taste-and-odor compounds, but there are biosynthetic gene clusters for other bioactive peptides, including anabaenopeptins, cyanopeptolins, microginins and ribosomally produced, post-translationally modified peptides. Genes for photosynthesis, nutrient acquisition, vitamin synthesis and buoyancy that are typical of bloom-forming cyanobacteria are present, although nitrate and nitrite reductase genes are conspicuously absent. However, the 7.9 Mbp genome is 3-4 Mbp larger than those of the above-mentioned frequently co-existing cyanobacteria. The increased genome size is largely due to an extraordinary number of insertion sequence elements (transposons), which account for 30.3% of the genome and many of which are present in multiple copies. The genome contains a relatively large number of pseudogenes, 97% of which are transposase genes. W. naegeliana WA131 thus seems to be able to limit the potentially deleterious effects of high rates of recombination and transposition to the mobilome fraction of its genome.

Differences in the regulatory strategies of marine oligotrophs and copiotrophs reflect differences in motility.

Aquatic bacteria frequently are divided into lifestyle categories oligotroph or copiotroph. Oligotrophs have proportionately fewer transcriptional regulatory genes than copiotrophs and are generally non-motile/chemotactic. We hypothesized that the absence of chemotaxis/motility in oligotrophs prevents them from occupying nutrient patches long enough to benefit from transcriptional regulation. We first confirmed that marine oligotrophs are generally reduced in genes for transcriptional regulation and motility/chemotaxis. Next, using a non-motile oligotroph (Ca. Pelagibacter st. HTCC7211), a motile copiotroph (Alteromonas macleodii st. HOT1A3), and [14 C]l-alanine, we confirmed that l-alanine catabolism is not transcriptionally regulated in HTCC7211 but is in HOT1A3. We then found that HOT1A3 took 2.5-4 min to initiate l-alanine oxidation at patch l-alanine concentrations, compared to <30 s for HTCC7211. By modelling cell trajectories, we predicted that, in most scenarios, non-motile cells spend <2 min in patches, compared to >4 min for chemotactic/motile cells. Thus, the time necessary for transcriptional regulation to initiate prevents transcriptional regulation from being beneficial for non-motile oligotrophs. This is supported by a mechanistic model we developed, which predicted that HTCC7211 cells with transcriptional regulation of l-alanine metabolism would produce 12% of their standing ATP stock upon encountering an l-alanine patch, compared to 880% in HTCC7211 cells without transcriptional regulation.

Age and micronutrient effects on the microbiome in a mouse model of zinc depletion and supplementation.

Older adult populations are at risk for zinc deficiency, which may predispose them to immune dysfunction and age-related chronic inflammation that drives myriad diseases and disorders. Recent work also implicates the gut microbiome in the onset and severity of age-related inflammation, indicating that dietary zinc status and the gut microbiome may interact to impact age-related host immunity. We hypothesize that age-related alterations in the gut microbiome contribute to the demonstrated zinc deficits in host zinc levels and increased inflammation. We tested this hypothesis with a multifactor two-part study design in a C57BL/6 mouse model. The two studies included young (2 month old) and aged (24 month old) mice fed either (1) a zinc adequate or zinc supplemented diet, or (2) a zinc adequate or marginal zinc deficient diet, respectively. Overall microbiome composition did not significantly change with zinc status; beta diversity was driven almost exclusively by age effects. Microbiome differences due to age are evident at all taxonomic levels, with more than half of all taxonomic units significantly different. Furthermore, we found 150 out of 186 genera were significantly different between the two age groups, with Bacteriodes and Parabacteroides being the primary taxa of young and old mice, respectively. These data suggest that modulating individual micronutrient concentrations does not lead to comprehensive microbiome shifts, but rather affects specific components of the gut microbiome. However, a phylogenetic agglomeration technique (ClaaTU) revealed phylogenetic clades that respond to modulation of dietary zinc status and inflammation state in an age-dependent manner. Collectively, these results suggest that a complex interplay exists between host age, gut microbiome composition, and dietary zinc status.

Genome sequence of freshwater nontoxigenic Limnoraphis associated with microcystin-producing blooms.

A sample from a 2019 cyanobacterial bloom in a freshwater reservoir in eastern Oregon, USA, was used to produce a metagenome from which the complete, circular 7.3 Mbp genome of Limnoraphis sp. WC205 was assembled. The Limnoraphis sp. WC205 genome contains gas vesicle genes, genes for N2-fixation and genes for both phycocyanin- and phycoerythrin-containing phycobilisomes. Limnoraphis was present in Willow Creek Reservoir throughout the summer and fall, coexisting with various other cyanobacteria in blooms that were associated with microcystin. The absence of cyanotoxin genes from the Limnoraphis sp. WC205 genome showed this cyanobacterium to be non-toxigenic, although it is predicted to produce cyanobactins closely related to Microcystis aeruginosa microcyclamides. DNA sequence corresponding to the Microcystis mcyG gene identified Microcystis as the microcystin producer in this lake.

Multi-Stimuli-Responsive Janus Hollow Polydopamine Nanotubes.

A tubular-shaped Janus nanoparticle based on polydopamine that responds to near-infrared, magnetic, and pH stimuli is reported. The robust tubular polydopamine structure was obtained by optimizing the halloysite template-to-dopamine ratio during synthesis. The inner and outer surfaces of the tube were exposed at different steps of the template-sonication--etching process, enabling the differential surface modification of these surfaces. Poly(ethylene glycol) (PEG) and poly(N-isopropylacrylamide) (PNIPAM) were grafted to the outer and inner surface of the nanotube, respectively. The PEG-coated surface limited aggregation of the nanoparticles at elevated temperatures. The PNIPAM-coated interior enhanced doxorubicin loading and endowed the nanoparticle with temperature-responsive behavior. The deposition of precipitated Fe3O4 nanoparticles further modified the nanoparticles. The resulting magnetic Janus nanoparticles responded to pH, temperature, and magnetic fields. Temperature changes could be induced by near-infrared laser, and all three stimuli were found to influence release rates of adsorbed doxorubicin from the nanoparticles. The interaction of the stimuli on release kinetics was elucidated using a linear mixed model; reduced pH and NIR irradiation enhanced release while applying a static magnetic field retarded release. Furthermore, the mechanism was shifted toward Fickian behavior by applying a static magnetic field and low pH conditions. However, NIR irradiation only shifted the behavior toward Fickian behavior at low pH.

7-epi-cylindrospermopsin and microcystin producers among diverse Anabaena/Dolichospermum/Aphanizomenon CyanoHABs in Oregon, USA.

Several genomes of Nostocales ADA clade members from the US Pacific Northwest were recently sequenced. Biosynthetic genes for microcystin, cylindrospermopsin or anatoxin-a were present in 7 of the 15 Dolichospermum/Anabaena strains and none of the 5 Aphanizomenon flos-aquae (AFA) strains. Toxin analyses (ELISA and LC-MS/MS) were conducted to quantitate and identify microcystin (MC) and cylindrospermopsin (CYN) congeners/analogs in samples dominated by Dolichospermum spp. of known genome sequence. MC-LR was the main congener produced by Dolichospermum spp. from Junipers Reservoir, Lake Billy Chinook and Odell Lake, while a congener provisionally identified as [Dha7]MC-HtyR was produced by a Dolichospermum sp. in Detroit Reservoir. A second Dolichospermum sp. from Detroit Reservoir was found to produce 7-epi-CYN, with 7-deoxy-CYN also present, but no CYN. The monitoring history of each of these lakes indicates the capacity for high levels of cyanotoxins during periods when Dolichospermum spp. are the dominant cyanobacteria. The diversity of ADA strains found in the US Pacific NW emphasizes the importance of these cyanobacteria as potentially toxic HAB formers in this temperate climatic region. Our results linking congener and genetic identity add data points that will help guide development of improved tools for predicting congener specificity from cyanotoxin gene sequences.

Metabolism of key atmospheric volatile organic compounds by the marine heterotrophic bacterium Pelagibacter HTCC1062 (SAR11).

Plants and phytoplankton are natural sources of the volatile organic compounds (VOCs) acetone and isoprene, which are reactive and can alter atmospheric chemistry. In earlier research we reported that, when co-cultured with a diatom, the marine bacterium Pelagibacter (strain HTCC1062; 'SAR11 clade') reduced the concentration of compounds tentatively identified as acetone and isoprene. In this study, experiments with Pelagibacter monocultures confirmed that these cells are capable of metabolizing acetone and isoprene at rates similar to bacterial communities in seawater and high enough to consume substantial fractions of the total marine acetone and isoprene budgets if extrapolated to global SAR11 populations. Homologues of an acetone/cyclohexanone monooxygenase were identified in the HTCC1062 genome and in the genomes of a wide variety of other abundant marine taxa, and were expressed at substantial levels (c. 10-4 of transcripts) across TARA oceans metatranscriptomes from ocean surface samples. The HTCC1062 genome lacks the canonical isoprene degradation pathway, suggesting an unknown alternative biochemical pathway is used by these cells for isoprene uptake. Fosmidomycin, an inhibitor of bacterial isoprenoid biosynthesis, blocked HTCC1062 growth, but the cells were rescued when isoprene was added to the culture, indicating SAR11 cells may be capable of synthesizing isoprenoid compounds from exogenous isoprene.

Interplay between Cruciferous Vegetables and the Gut Microbiome: A Multi-Omic Approach.

Brassica vegetables contain a multitude of bioactive compounds that prevent and suppress cancer and promote health. Evidence suggests that the gut microbiome may be essential in the production of these compounds; however, the relationship between specific microbes and the abundance of metabolites produced during cruciferous vegetable digestion are still unclear. We utilized an ex vivo human fecal incubation model with in vitro digested broccoli sprouts (Broc), Brussels sprouts (Brus), a combination of the two vegetables (Combo), or a negative control (NC) to investigate microbial metabolites of cruciferous vegetables. We conducted untargeted metabolomics on the fecal cultures by LC-MS/MS and completed 16S rRNA gene sequencing. We identified 72 microbial genera in our samples, 29 of which were significantly differentially abundant between treatment groups. A total of 4499 metabolomic features were found to be significantly different between treatment groups (q ≤ 0.05, fold change > 2). Chemical enrichment analysis revealed 45 classes of compounds to be significantly enriched by brassicas, including long-chain fatty acids, coumaric acids, and peptides. Multi-block PLS-DA and a filtering method were used to identify microbe−metabolite interactions. We identified 373 metabolites from brassica, which had strong relationships with microbes, such as members of the family Clostridiaceae and genus Intestinibacter, that may be microbially derived.

Bacillus velezensis AP183 Inhibits Staphylococcus aureus Biofilm Formation and Proliferation in Murine and Bovine Disease Models.

The increasing frequency of S. aureus antimicrobial resistance has spurred interest in identifying alternative therapeutants. We investigated the S. aureus-inhibitory capacity of B. velezensis strains in mouse and bovine models. Among multiple B. velezensis strains that inhibited S. aureus growth in vitro, B. velezensis AP183 provided the most potent inhibition of S. aureus proliferation and bioluminescence in a mouse cutaneous wound (P = 0.02). Histology revealed abundant Gram-positive cocci in control wounds that were reduced in B. velezensis AP183-treated tissues. Experiments were then conducted to evaluate the ability of B. velezensis AP183 to prevent S. aureus biofilm formation on a tracheostomy tube substrate. B. velezensis AP183 could form a biofilm on a tracheostomy tube inner cannula substrate, and that this biofilm was antagonistic to S. aureus colonization. B. velezensis AP183 was also observed to inhibit the growth of S. aureus isolates originated from bovine mastitis cases. To evaluate the inflammatory response of mammary tissue to intramammary inoculation with B. velezensis AP183, we used high dose and low dose inocula in dairy cows. At the high dose, a significant increase in somatic cell count (SCC) and clinical mastitis was observed at all post-inoculation time points (P < 0.01), which resolved quickly compared to S. aureus-induced mastitis; in contrast, the lower dose of B. velezensis AP183 resulted in a slight increase of SCC and no clinical mastitis. In a subsequent experiment, all mammary quarters in four cows were induced to have grade 1 clinical mastitis by intramammary inoculation of a S. aureus mastitis isolate; following mastitis induction, eight quarters were treated with B. velezensis AP183 and milk samples were collected from pretreatment and post-treatment samples for 9 days. In groups treated with B. velezensis AP183, SCC and abundance of S. aureus decreased with significant reductions in S. aureus after 3 days post-inoculation with AP183 (P = 0.04). A milk microbiome analysis revealed significant reductions in S. aureus relative abundance in the AP183-treated group by 8 days post-inoculation (P = 0.02). These data indicate that B. velezensis AP183 can inhibit S. aureus biofilm formation and its proliferation in murine and bovine disease models.

Diversification of the Type VI Secretion System in Agrobacteria.

The type VI secretion system (T6SS) is used by many Gram-negative bacteria to deploy toxic effectors for interbacterial competition. This system provides a competitive advantage in planta to agrobacteria, a diverse group with phytopathogenic members capable of genetically transforming plants. To inform on the ecology and evolution of agrobacteria, we revealed processes that diversify their effector gene collections. From genome sequences of diverse strains, we identified T6SS loci, functionally validated associated effector genes for toxicity, and predicted genes homologous to those that encode proteins known to interact with effectors. The gene loci were analyzed in a phylogenetic framework, and results show that strains of some species-level groups have different patterns of T6SS expression and are enriched in specific sets of T6SS loci. Findings also demonstrate that the modularity of T6SS loci and their associated genes engenders dynamicity, promoting reshuffling of entire loci, fragments therein, and domains to swap toxic effector genes across species. However, diversification is constrained by the need to maintain specific combinations of gene subtypes, congruent with observations that certain genes function together to regulate T6SS loading and activation. Data are consistent with a scenario where species can acquire unique T6SS loci that are then reshuffled across the genus in a restricted manner to generate new combinations of effector genes. IMPORTANCE The T6SS is used by several taxa of Gram-negative bacteria to secrete toxic effector proteins to attack others. Diversification of effector collections shapes bacterial interactions and impacts the health of hosts and ecosystems in which bacteria reside. We uncovered the diversity of T6SS loci across a genus of plant-associated bacteria and show that diversification is driven by the acquisition of new loci and reshuffling among species. However, linkages between specific subtypes of genes need to be maintained to ensure that proteins whose interactions are necessary to activate the T6SS remain together. Results reveal how organization of gene loci and domain structure of genes provides flexibility to diversify under the constraints imposed by the system. Findings inform on the evolution of a mechanism that influences bacterial communities.

Complete genomes derived by directly sequencing freshwater bloom populations emphasize the significance of the genus level ADA clade within the Nostocales.

The genome sequences of 16 Nostocales cyanobacteria have been determined. Most of them are complete or near-complete genome sequences derived by long-read metagenome sequencing of recent harmful algal blooms (HABs) in freshwater lakes without the potential bias of culture isolation. The genomes are all members of the recently recognized ADA clade (Driscoll et al., Harmful Algae, 77:93, 2018), which we argue represents a genus. We identify 10 putative species-level branches within the clade, on the basis of 91-gene phylogenomic and average nucleotide identity analyses. The assembled genomes each correspond to a single morphotype in the original sample, but distinct genomes from different HABs in some cases correspond to similar morphotypes. We present data indicating that the ADA clade is a highly significant component of current cyanobacterial HABs, including members assigned to the prevalent Dolichospermum and Aphanizomenon genera, as well as Cuspidothrix and Anabaena. In general, currently used genus and species names within the ADA clade are not monophyletic. We infer that the morphological characters routinely used in taxonomic assignments are not reliable for discriminating species within the ADA clade. Taxonomic revisions will be needed to create a genus with a single name (we recommend Anabaena) and to adopt species names that do not depend on morphological traits that lack sufficient discrimination and specificity, while recognizing the utility of some easily observable and distinct morphologies.

Comparative genomics of the ADA clade within the Nostocales.

The ADA clade of Nostocales cyanobacteria, a group that is prominent in current harmful algal bloom events, now includes over 40 genome sequences with the recent addition of sixteen novel sequenced genomes (Dreher et al., Harmful Algae, 2021). Fourteen genomes are complete (closed), enabling highly detailed assessments of gene content and genome architecture. ADA genomes contain 5 rRNA operons, genes expected to support a photoautotrophic and diazotrophic lifestyle, and a varied array of genes for the synthesis of bioactive secondary metabolites. Genes for the production of the taste-and-odor compound geosmin and the four major classes of cyanotoxins - anatoxin-a, cylindrospermopsin, microcystin and saxitoxin - are represented in members of the ADA clade. Notably, the gene array for the synthesis of cylindrospermopsin by Dolichospermum sp. DET69 was located on a plasmid, raising the possibility of facile horizontal transmission. However, genes supporting independent conjugative transfer of this plasmid are lacking. Further, analysis of genomic loci containing this and other cyanotoxin gene arrays shows evidence that these arrays have long-term stability and do not appear to be genomic islands easily capable of horizontal transmission to other cells. There is considerable diversity in the gene complements of individual ADA genomes, including the variable presence of physiologically important genes: genomes in three species-level subclades lack the gas vesicle genes that facilitate a planktonic lifestyle, and, surprisingly, the genome of Cuspidothrix issatschenkoi CHARLIE-1, a reported diazotroph, lacks the genes for nitrogen fixation. Notably, phylogenetically related genomes possess limited synteny, indicating a prominent role for chromosome rearrangements during ADA strain evolution. The genomes contain abundant insertion sequences and repetitive transposase genes, which could be the main drivers of genome rearrangement through active transposition and homologous recombination. No prophages were found, and no evidence of viral infection was observed in the bloom population samples from which the genomes discussed here were derived. Phages thus seem to have a limited influence on ADA evolution.

Unexpected distribution of the 4-formylaminooxyvinylglycine (FVG) biosynthetic pathway in Pseudomonas and beyond.

The biological herbicide and antibiotic 4-formylaminooxyvinylglycine (FVG) was originally isolated from several rhizosphere-associated strains of Pseudomonas fluorescens. Biosynthesis of FVG is dependent on the gvg biosynthetic gene cluster in P. fluorescens. In this investigation, we used comparative genomics to identify strains with the genetic potential to produce FVG due to presence of a gvg gene cluster. These strains primarily belong to two groups of Pseudomonas, P. fluorescens and P. syringae, however, a few strains with the gvg cluster were found outside of Pseudomonas. Mass spectrometry confirmed that all tested strains of the P. fluorescens species group produced FVG. However, P. syringae strains did not produce FVG under standard conditions. Several lines of evidence regarding the transmission of the gvg cluster including a robust phylogenetic analysis suggest that it was introduced multiple times through horizontal gene transfer within the Pseudomonas lineage as well as in select lineages of Thiomonas, Burkholderia and Pantoea. Together, these data broaden our understanding of the evolution and diversity of FVG biosynthesis. In the course of this investigation, additional gene clusters containing only a subset of the genes required to produce FVG were identified in a broad range of bacteria, including many non-pseudomonads.

Using a lake sediment record to infer the long-term history of cyanobacteria and the recent rise of an anatoxin producing Dolichospermum sp.

Lakes that experience recurrent toxic cyanobacterial harmful algae blooms (cyanoHABS) are often subject to cultural eutrophication, where landscape development and upland activities increase the nutrient inputs to the water column and fuel cyanoHABS. Few studies have focused on the response of a lake to nutrient inputs for which the natural geomorphic setting predisposes a nutrient-rich water column to already support abundant cyanobacteria. Here, we present a sediment core record from a lake surrounded by parkland that experiences recurrent cyanoHABs which produce dangerous levels of the neurotoxin, anatoxin-a, impacting the recreational use of the lake and park. Using photoautotrophic pigments in the sediment record, we establish cyanobacteria have long been part of the diverse and abundant phytoplankton community within the lake. Despite this long record, shotgun metagenome and other DNA analyses of the sediment record suggest that the current anatoxin-a producer Dolichospermum sp. WA102 only emerged to dominate the cyanobacterial community in the mid-1990s. A period of lakeshore farming that finished in the 1950s-1960s and possibly the stocking of rainbow trout fry (1970-2016) coincide with a progressive shift in primary production, together with a change in bacterial communities. Based on the history of the lake and contemporary ecology of Dolichospermum, we propose that the legacy of nutrient inputs and changes in nutrient cycling within the lake has encouraged the development of an ecosystem where the toxin producing Dolichospermum sp WA102 is highly competitive. Understanding the historical presence of cyanobacteria in the lake provides a context for current-day management strategies of cyanoHABs.

Phylogenetic Integration Reveals the Zebrafish Core Microbiome and Its Sensitivity to Environmental Exposures.

Zebrafish are increasingly used to study how environmental exposures impact vertebrate gut microbes. However, we understand little about which microbial taxa are common to the zebrafish gut across studies and facilities. Here, we define the zebrafish core gut microbiome to resolve microbiota that are both relatively robust to study or facility effects and likely to drive proper microbiome assembly and functioning due to their conservation. To do so, we integrated publicly available gut microbiome 16S gene sequence data from eight studies into a phylogeny and identified monophyletic clades of gut bacteria that are unexpectedly prevalent across individuals. Doing so revealed 585 core clades of bacteria in the zebrafish gut, including clades within Aeromonas, Pseudomonas, Cetobacterium, Shewanella, Chitinibacter, Fluviicola, Flectobacillus, and Paucibacter. We then applied linear regression to discern which of these core clades are sensitive to an array of different environmental exposures. We found that 200 core clades were insensitive to any exposure we assessed, while 134 core clades were sensitive to more than two exposures. Overall, our analysis defines the zebrafish core gut microbiome and its sensitivity to exposure, which helps future studies to assess the robustness of their results and prioritize taxa for empirical assessments of how gut microbiota mediate the effects of exposure on the zebrafish host.

Ancient co-option of an amino acid ABC transporter locus in Pseudomonas syringae for host signal-dependent virulence gene regulation.

Pathogenic bacteria frequently acquire virulence traits via horizontal gene transfer, yet additional evolutionary innovations may be necessary to integrate newly acquired genes into existing regulatory pathways. The plant bacterial pathogen Pseudomonas syringae relies on a horizontally acquired type III secretion system (T3SS) to cause disease. T3SS-encoding genes are induced by plant-derived metabolites, yet how this regulation occurs, and how it evolved, is poorly understood. Here we report that the two-component system AauS-AauR and substrate-binding protein AatJ, proteins encoded by an acidic amino acid-transport (aat) and -utilization (aau) locus in P. syringae, directly regulate T3SS-encoding genes in response to host aspartate and glutamate signals. Mutants of P. syringae strain DC3000 lacking aauS, aauR or aatJ expressed lower levels of T3SS genes in response to aspartate and glutamate, and had decreased T3SS deployment and virulence during infection of Arabidopsis. We identified an AauR-binding motif (Rbm) upstream of genes encoding T3SS regulators HrpR and HrpS, and demonstrated that this Rbm is required for maximal T3SS deployment and virulence of DC3000. The Rbm upstream of hrpRS is conserved in all P. syringae strains with a canonical T3SS, suggesting AauR regulation of hrpRS is ancient. Consistent with a model of conserved function, an aauR deletion mutant of P. syringae strain B728a, a bean pathogen, had decreased T3SS expression and growth in host plants. Together, our data suggest that, upon acquisition of T3SS-encoding genes, a strain ancestral to P. syringae co-opted an existing AatJ-AauS-AauR pathway to regulate T3SS deployment in response to specific host metabolite signals.

Unexpected conservation and global transmission of agrobacterial virulence plasmids.

The accelerated evolution and spread of pathogens are threats to host species. Agrobacteria require an oncogenic Ti or Ri plasmid to transfer genes into plants and cause disease. We developed a strategy to characterize virulence plasmids and applied it to analyze hundreds of strains collected between 1927 and 2017, on six continents and from more than 50 host species. In consideration of prior evidence for prolific recombination, it was surprising that oncogenic plasmids are descended from a few conserved lineages. Characterization of a hierarchy of features that promote or constrain plasticity allowed inference of the evolutionary history across the plasmid lineages. We uncovered epidemiological patterns that highlight the importance of plasmid transmission in pathogen diversification as well as in long-term persistence and the global spread of disease.

Pan-tissue transcriptome analysis of long noncoding RNAs in the American beaver Castor canadensis.

BACKGROUND: Long noncoding RNAs (lncRNAs) have roles in gene regulation, epigenetics, and molecular scaffolding and it is hypothesized that they underlie some mammalian evolutionary adaptations. However, for many mammalian species, the absence of a genome assembly precludes the comprehensive identification of lncRNAs. The genome of the American beaver (Castor canadensis) has recently been sequenced, setting the stage for the systematic identification of beaver lncRNAs and the characterization of their expression in various tissues. The objective of this study was to discover and profile polyadenylated lncRNAs in the beaver using high-throughput short-read sequencing of RNA from sixteen beaver tissues and to annotate the resulting lncRNAs based on their potential for orthology with known lncRNAs in other species. RESULTS: Using de novo transcriptome assembly, we found 9528 potential lncRNA contigs and 187 high-confidence lncRNA contigs. Of the high-confidence lncRNA contigs, 147 have no known orthologs (and thus are putative novel lncRNAs) and 40 have mammalian orthologs. The novel lncRNAs mapped to the Oregon State University (OSU) reference beaver genome with greater than 90% sequence identity. While the novel lncRNAs were on average shorter than their annotated counterparts, they were similar to the annotated lncRNAs in terms of the relationships between contig length and minimum free energy (MFE) and between coverage and contig length. We identified beaver orthologs of known lncRNAs such as XIST, MEG3, TINCR, and NIPBL-DT. We profiled the expression of the 187 high-confidence lncRNAs across 16 beaver tissues (whole blood, brain, lung, liver, heart, stomach, intestine, skeletal muscle, kidney, spleen, ovary, placenta, castor gland, tail, toe-webbing, and tongue) and identified both tissue-specific and ubiquitous lncRNAs. CONCLUSIONS: To our knowledge this is the first report of systematic identification of lncRNAs and their expression atlas in beaver. LncRNAs-both novel and those with known orthologs-are expressed in each of the beaver tissues that we analyzed. For some beaver lncRNAs with known orthologs, the tissue-specific expression patterns were phylogenetically conserved. The lncRNA sequence data files and raw sequence files are available via the web supplement and the NCBI Sequence Read Archive, respectively.