Swabs to genomes: a comprehensive workflow

Screening Salamanders for Symbionts

Microbial symbionts are broadly categorized by their impacts on host fitness: commensals, pathogens, and mutualists. However, recent investigations into the physiological basis of these impacts have revealed nuanced microbial influences on a wide range of host developmental, immunological, and physiological processes, including regeneration. Exploring these impacts begins with knowing which microbes are present. This methodological pipeline contains both targeted assays using PCR and culturing, as well as culture-independent approaches, to survey host salamander tissues for common and unknown microbial symbionts.

Deconstructing and Reconstructing Cheese Rind Microbiomes for Experiments in Microbial Ecology and Evolution

Cheese rind microbiomes are useful model systems for identifying the mechanisms that control microbiome diversity. Here, we describe the methods we have optimized to first deconstruct in situ cheese rind microbiome diversity and then reconstruct that diversity in laboratory environments to conduct controlled microbiome manipulations. Most cheese rind microbial species, including bacteria, yeasts, and filamentous fungi, can be easily cultured using standard lab media. Colony morphologies of taxa are diverse and can often be used to distinguish taxa at the phylum and sometimes even genus level. Through the use of cheese curd agar medium, thousands of unique community combinations or microbial interactions can be assessed. Transcriptomic experiments and transposon mutagenesis screens can pinpoint mechanisms of interactions between microbial species. Our general approach of creating a tractable synthetic microbial community from cheese can be easily applied to other fermented foods to develop other model microbiomes. © 2019 by John Wiley & Sons, Inc. Basic Protocol 1: Isolation of cheese rind microbial communities Support Protocol 1: Preparation of plate count agar with milk and salt Basic Protocol 2: Identification of cheese rind bacterial and fungal isolates using 16S and ITS sequences Basic Protocol 3: Preparation of experimental glycerol stocks of yeasts and bacteria Basic Protocol 4: Preparation of experimental glycerol stocks of filamentous fungi Basic Protocol 5: Reconstruction of cheese rind microbial communities in vitro Support Protocol 2: Preparation of lyophilized and powdered cheese curd Support Protocol 3: Preparation of 10% cheese curd agar plates and tubes Basic Protocol 6: Interaction screens using responding lawns Support Protocol 4: Preparation of liquid 2% cheese curd Basic Protocol 7: Experimental evolution Basic Protocol 8: Measuring community function: pH/acidification Basic Protocol 9: Measuring community function: Pigment production Basic Protocol 10: RNA sequencing of cheese rind biofilms.

Combining Microbial Culturing With Mathematical Modeling in an Introductory Course-Based Undergraduate Research Experience

Quantitative techniques are a critical part of contemporary biology research, but students interested in biology enter college with widely varying quantitative skills and attitudes toward mathematics. Course-based undergraduate research experiences (CUREs) may be an early way to build student competency and positive attitudes. Here we describe the design, implementation, and assessment of an introductory quantitative CURE focused on halophilic microbes. In this CURE, students culture and isolate halophilic microbes from environmental and food samples, perform growth assays, then use mathematical modeling to quantify the growth rate of strains in different salinities. To assess how the course may impact students' future academic plans and attitudes toward the use of math in biology, we used pre- and post-quarter surveys. Students who completed the course showed more positive attitudes toward science learning and an increased interest in pursuing additional quantitative biology experiences. We argue that the classroom application of microbiology methods, combined with mathematical modeling using student-generated data, provides a degree of student ownership, collaboration, iteration, and discovery that makes quantitative learning both relevant and exciting to students.

Isolation and sequence-based characterization of a koala symbiont: Lonepinella koalarum

Koalas (Phascolarctos cinereus) are highly specialized herbivorous marsupials that feed almost exclusively on Eucalyptus leaves, which are known to contain varying concentrations of many different toxic chemical compounds. The literature suggests that Lonepinella koalarum, a bacterium in the Pasteurellaceae family, can break down some of these toxic chemical compounds. Furthermore, in a previous study, we identified L. koalarum as the most predictive taxon of koala survival during antibiotic treatment. Therefore, we believe that this bacterium may be important for koala health. Here, we isolated a strain of L. koalarum from a healthy koala female and sequenced its genome using a combination of short-read and long-read sequencing. We placed the genome assembly into a phylogenetic tree based on 120 genome markers using the Genome Taxonomy Database (GTDB), which currently does not include any L. koalarum assemblies. Our genome assembly fell in the middle of a group of Haemophilus, Pasteurella and Basfia species. According to average nucleotide identity and a 16S rRNA gene tree, the closest relative of our isolate is L. koalarum strain Y17189. Then, we annotated the gene sequences and compared them to 55 closely related, publicly available genomes. Several genes that are known to be involved in carbohydrate metabolism could exclusively be found in L. koalarum relative to the other taxa in the pangenome, including glycoside hydrolase families GH2, GH31, GH32, GH43 and GH77. Among the predicted genes of L. koalarum were 79 candidates putatively involved in the degradation of plant secondary metabolites. Additionally, several genes coding for amino acid variants were found that had been shown to confer antibiotic resistance in other bacterial species against pulvomycin, beta-lactam antibiotics and the antibiotic efflux pump KpnH. In summary, this genetic characterization allows us to build hypotheses to explore the potentially beneficial role that L. koalarum might play in the koala intestinal microbiome. Characterizing and understanding beneficial symbionts at the whole genome level is important for the development of anti- and probiotic treatments for koalas, a highly threatened species due to habitat loss, wildfires, and high prevalence of Chlamydia infections.

Fungi, bacteria and oomycota opportunistically isolated from the seagrass, Zostera marina

Fungi in the marine environment are often neglected as a research topic, despite that fungi having critical roles on land as decomposers, pathogens or endophytes. Here we used culture-dependent methods to survey the fungi associated with the seagrass, Zostera marina, also obtaining bacteria and oomycete isolates in the process. A total of 108 fungi, 40 bacteria and 2 oomycetes were isolated. These isolates were then taxonomically identified using a combination of molecular and phylogenetic methods. The majority of the fungal isolates were classified as belonging to the classes Eurotiomycetes, Dothideomycetes, and Sordariomycetes. Most fungal isolates were habitat generalists like Penicillium sp. and Cladosporium sp., but we also cultured a diverse set of rare taxa including possible habitat specialists like Colletotrichum sp. which may preferentially associate with Z. marina leaf tissue. Although the bulk of bacterial isolates were identified as being from known ubiquitous marine lineages, we also obtained several Actinomycetes isolates and a Phyllobacterium sp. We identified two oomycetes, another understudied group of marine microbial eukaryotes, as Halophytophthora sp. which may be opportunistic pathogens or saprophytes of Z. marina. Overall, this study generates a culture collection of fungi which adds to knowledge of Z. marina associated fungi and highlights a need for more investigation into the functional and evolutionary roles of microbial eukaryotes associated with seagrasses.

Characterization of the Mycobiome of the Seagrass, Zostera marina, Reveals Putative Associations With Marine Chytrids

Seagrasses are globally distributed marine flowering plants that are foundation species in coastal ecosystems. Seagrass beds play essential roles as habitats and hatcheries, in nutrient cycling, and in protecting the coastline from erosion. Although many studies have focused on seagrass ecology, only a limited number have investigated their associated fungi. In terrestrial systems, fungi can have beneficial and detrimental effects on plant fitness. However, not much is known about marine fungi and even less is known about seagrass associated fungi. Here we used culture-independent sequencing of the ribosomal internal transcribed spacer (ITS) region to characterize the taxonomic diversity of fungi associated with the seagrass, Zostera marina. We sampled from two Z. marina beds in Bodega Bay over three time points to investigate fungal diversity within and between plants. Our results indicate that there are many fungal taxa for which a taxonomic assignment cannot be made living on and inside Z. marina leaves, roots and rhizomes and that these plant tissues harbor distinct fungal communities. We also identified differences in the abundances of the orders, Glomerellales, Agaricales and Malasseziales, between seagrass tissues. The most prevalent ITS amplicon sequence variants (ASVs) associated with Z. marina tissues could not initially be confidently assigned to a fungal phylum, but shared significant sequence similarity with Chytridiomycota and Aphelidomycota. To obtain a more definitive taxonomic classification of the most abundant ASV associated with Z. marina leaves, we used PCR with one primer targeting a unique region of this ASV's ITS2 and a second primer targeting fungal 28S rRNA genes to amplify part of the 28S rRNA gene region corresponding to this ASV. Sequencing and phylogenetic analysis of the resulting partial 28S rRNA gene revealed that the organism that this ASV comes from is a member of Novel Clade SW-I in the order Lobulomycetales in the phylum Chytridiomycota. This clade includes known parasites of freshwater diatoms and algae and it is possible this chytrid is directly infecting Z. marina leaf tissues. This work highlights a need for further studies focusing on marine fungi and the potential importance of these understudied communities to the larger seagrass ecosystem.

Genomes from bacteria associated with the canine oral cavity: A test case for automated genome-based taxonomic assignment

Taxonomy for bacterial isolates is commonly assigned via sequence analysis. However, the most common sequence-based approaches (e.g. 16S rRNA gene-based phylogeny or whole genome comparisons) are still labor intensive and subjective to varying degrees. Here we present a set of 33 bacterial genomes, isolated from the canine oral cavity. Taxonomy of these isolates was first assigned by PCR amplification of the 16S rRNA gene, Sanger sequencing, and taxonomy assignment using BLAST. After genome sequencing, taxonomy was revisited through a manual process using a combination of average nucleotide identity (ANI), concatenated marker gene phylogenies, and 16S rRNA gene phylogenies. This taxonomy was then compared to the automated taxonomic assignment given by the recently proposed Genome Taxonomy Database (GTDB). We found the results of all three methods to be similar (25 out of the 33 had matching genera), but the GTDB approach required fewer subjective decisions, and required far less labor. The primary differences in the non-identical taxonomic assignments involved cases where GTDB has proposed taxonomic revisions.

Automated detection of records in biological sequence databases that are inconsistent with the literature

We investigate and analyse the data quality of nucleotide sequence databases with the objective of automatic detection of data anomalies and suspicious records. Specifically, we demonstrate that the published literature associated with each data record can be used to automatically evaluate its quality, by cross-checking the consistency of the key content of the database record with the referenced publications. Focusing on GenBank, we describe a set of quality indicators based on the relevance paradigm of information retrieval (IR). Then, we use these quality indicators to train an anomaly detection algorithm to classify records as "confident" or "suspicious". Our experiments on the PubMed Central collection show assessing the coherence between the literature and database records, through our algorithms, is an effective mechanism for assisting curators to perform data cleansing. Although fewer than 0.25% of the records in our data set are known to be faulty, we would expect that there are many more in GenBank that have not yet been identified. By automated comparison with literature they can be identified with a precision of up to 10% and a recall of up to 30%, while strongly outperforming several baselines. While these results leave substantial room for improvement, they reflect both the very imbalanced nature of the data, and the limited explicitly labelled data that is available. Overall, the obtained results show promise for the development of a new kind of approach to detecting low-quality and suspicious sequence records based on literature analysis and consistency. From a practical point of view, this will greatly help curators in identifying inconsistent records in large-scale sequence databases by highlighting records that are likely to be inconsistent with the literature.

Antifungal Bacteria on Woodland Salamander Skin Exhibit High Taxonomic Diversity and Geographic Variability

Diverse bacteria inhabit amphibian skin; some of those bacteria inhibit growth of the fungal pathogen Batrachochytrium dendrobatidis Yet there has been no systematic survey of anti-B. dendrobatidis bacteria across localities, species, and elevations. This is important given geographic and taxonomic variations in amphibian susceptibility to B. dendrobatidis Our collection sites were at locations within the Appalachian Mountains where previous sampling had indicated low B. dendrobatidis prevalence. We determined the numbers and identities of anti-B. dendrobatidis bacteria on 61 Plethodon salamanders (37 P. cinereus, 15 P. glutinosus, 9 P. cylindraceus) via culturing methods and 16S rRNA gene sequencing. We sampled co-occurring species at three localities and sampled P. cinereus along an elevational gradient (700 to 1,000 meters above sea level [masl]) at one locality. We identified 50 anti-B. dendrobatidis bacterial operational taxonomic units (OTUs) and found that the degree of B. dendrobatidis inhibition was not correlated with relatedness. Five anti-B. dendrobatidis bacterial strains occurred on multiple amphibian species at multiple localities, but none were shared among all species and localities. The prevalence of anti-B. dendrobatidis bacteria was higher at Shenandoah National Park (NP), VA, with 96% (25/26) of salamanders hosting at least one anti-B. dendrobatidis bacterial species compared to 50% (7/14) at Catoctin Mountain Park (MP), MD, and 38% (8/21) at Mt. Rogers National Recreation Area (NRA), VA. At the individual level, salamanders at Shenandoah NP had more anti-B. dendrobatidis bacteria per individual (μ = 3.3) than those at Catoctin MP (μ = 0.8) and at Mt. Rogers NRA (μ = 0.4). All salamanders tested negative for B. dendrobatidis Anti-B. dendrobatidis bacterial species are diverse in central Appalachian Plethodon salamanders, and their distribution varied geographically. The antifungal bacterial species that we identified may play a protective role for these salamanders.IMPORTANCE Amphibians harbor skin bacteria that can kill an amphibian fungal pathogen, Batrachochytrium dendrobatidis Some amphibians die from B. dendrobatidis infection, whereas others do not. The bacteria that can kill B. dendrobatidis, called anti-B. dendrobatidis bacteria, are thought to influence the B. dendrobatidis infection outcome for the amphibian. Yet how anti-B. dendrobatidis bacterial species vary among amphibian species and populations is unknown. We determined the distribution of anti-B. dendrobatidis bacterial species among three salamander species (n = 61) sampled at three localities. We identified 50 unique anti-B. dendrobatidis bacterial species and found that all of the tested salamanders were negative for B. dendrobatidis Five anti-B. dendrobatidis bacterial species were commonly detected, suggesting a stable, functional association with these salamanders. The number of anti-B. dendrobatidis bacteria per individual varied among localities but not among co-occurring salamander species, demonstrating that environment is more influential than host factors in structuring the anti-B. dendrobatidis bacterial community. These anti-B. dendrobatidis bacteria may serve a protective function for their salamander hosts.

Draft Genome Sequences of Pseudomonas moraviensis UCD-KL30, Vibrio ostreicida UCD-KL16, Colwellia sp. Strain UCD-KL20, Shewanella sp. Strain UCD-KL12, and Shewanella sp. Strain UCD-KL21, Isolated from Seagrass

Here, we present the draft genome sequences for five bacterial strains. These strains were all isolated from seagrass (Zostera marina) collected from Bodega Bay, CA, as a part of an undergraduate research project focused on seagrass-associated microbes.

Draft Genome Sequence of a Multidrug-Resistant Strain of Enterococcus faecalis, PM01, Isolated from the Nest of an American Bushtit, Psaltriparius minimus

Pathogenic microorganisms associated with avian nests may detrimentally impact parental health and nest success for the nest primary users, potentially neighboring avian or terrestrial species, including humans. Here, we report the genome sequence of Enterococcus faecalis strain PM01, isolated from a failed nest of American bushtits, Psaltriparius minimus.

Draft genome sequences of eight bacteria isolated from the indoor environment: Staphylococcus capitis strain H36, S. capitis strain H65, S. cohnii strain H62, S. hominis strain H69, Microbacterium sp. strain H83, Mycobacterium iranicum strain H39, Plantibacter sp. strain H53, and Pseudomonas oryzihabitans strain H72

We report here the draft genome sequences of eight bacterial strains of the genera Staphylococcus, Microbacterium, Mycobacterium, Plantibacter, and Pseudomonas. These isolates were obtained from aerosol sampling of bathrooms of five residences in the San Francisco Bay area. Taxonomic classifications as well as the genome sequence and gene annotation of the isolates are described. As part of the "Built Environment Reference Genome" project, these isolates and associated genome data provide valuable resources for studying the microbiology of the built environment.

Reconstructing the Ancestral Relationships Between Bacterial Pathogen Genomes

Following recent developments in DNA sequencing technology, it is now possible to sequence hundreds of whole genomes from bacterial isolates at relatively low cost. Analyzing this growing wealth of genomic data in terms of ancestral relationships can reveal many interesting aspects of the evolution, ecology, and epidemiology of bacterial pathogens. However, reconstructing the ancestry of a sample of bacteria remains challenging, especially for the majority of species where recombination is frequent. Here, we review and describe the computational techniques currently available to infer ancestral relationships, including phylogenetic methods that either ignore or account for the effect of recombination, as well as model-based and model-free phylogeny-independent approaches.

Functional Annotations of Paralogs: A Blessing and a Curse

Gene duplication followed by mutation is a classic mechanism of neofunctionalization, producing gene families with functional diversity. In some cases, a single point mutation is sufficient to change the substrate specificity and/or the chemistry performed by an enzyme, making it difficult to accurately separate enzymes with identical functions from homologs with different functions. Because sequence similarity is often used as a basis for assigning functional annotations to genes, non-isofunctional gene families pose a great challenge for genome annotation pipelines. Here we describe how integrating evolutionary and functional information such as genome context, phylogeny, metabolic reconstruction and signature motifs may be required to correctly annotate multifunctional families. These integrative analyses can also lead to the discovery of novel gene functions, as hints from specific subgroups can guide the functional characterization of other members of the family. We demonstrate how careful manual curation processes using comparative genomics can disambiguate subgroups within large multifunctional families and discover their functions. We present the COG0720 protein family as a case study. We also discuss strategies to automate this process to improve the accuracy of genome functional annotation pipelines.

Genome Sequence of a Multidrug-Resistant Strain of Bacillus pumilus, CB01, Isolated from the Feces of an American Crow, Corvus brachyrhynchos

Avian species have the potential to serve as important reservoirs for the spread of pathogenic microorganisms. Here, we report the genome sequence of a drug-resistant strain of Bacillus pumilus, CB01, isolated from the feces of an American crow, Corvus brachyrhynchos.

Growth of 48 built environment bacterial isolates on board the International Space Station (ISS)

Background. While significant attention has been paid to the potential risk of pathogenic microbes aboard crewed spacecraft, the non-pathogenic microbes in these habitats have received less consideration. Preliminary work has demonstrated that the interior of the International Space Station (ISS) has a microbial community resembling those of built environments on Earth. Here we report the results of sending 48 bacterial strains, collected from built environments on Earth, for a growth experiment on the ISS. This project was a component of Project MERCCURI (Microbial Ecology Research Combining Citizen and University Researchers on ISS). Results. Of the 48 strains sent to the ISS, 45 of them showed similar growth in space and on Earth using a relative growth measurement adapted for microgravity. The vast majority of species tested in this experiment have also been found in culture-independent surveys of the ISS. Only one bacterial strain showed significantly different growth in space. Bacillus safensis JPL-MERTA-8-2 grew 60% better in space than on Earth. Conclusions. The majority of bacteria tested were not affected by conditions aboard the ISS in this experiment (e.g., microgravity, cosmic radiation). Further work on Bacillus safensis could lead to interesting insights on why this strain grew so much better in space.

Draft Genome Sequences of Two Pseudoalteromonas porphyrae Strains Isolated from Seagrass Sediment

Here, we present the draft genome sequences of Pseudoalteromonas porphyrae UCD-SED9 and UCD-SED14 (phylum Proteobacteria). These strains were isolated from sediment surrounding the roots of the seagrass, Zostera marina, collected near the UC, Davis Bodega Marine Laboratory (Bodega Bay, California). The assemblies contain 4,847,456 bp and 4,817,752 bp, respectively.

Draft Genome Sequences of Two Vibrio splendidus Strains, Isolated from Seagrass Sediment

Here, we present the draft genome sequences of Vibrio splendidus UCD-SED7 and UCD-SED10 (phylum Proteobacteria ). These strains were isolated from sediment surrounding Zostera marina roots near the UC Davis Bodega Marine Laboratory (Bodega, Bay, California). These assemblies contain 5,334,236 bp and 5,904,824 bp, respectively.

Draft Genome Sequence of Bacillus vietnamensis Strain UCD-SED5 (Phylum Firmicutes)

Here, we present the draft genome sequence of Bacillus vietnamensis UCD-SED5 (phylum Firmicutes). This strain was isolated from sediment surrounding Zostera marina roots near the UC Davis Bodega Marine Laboratory (Bodega, Bay, California) and represents the second genome of this species. The assembly consists of 4,325,707 bp, in 108 contigs.

Porphyrobacter mercurialis sp. nov., isolated from a stadium seat and emended description of the genus Porphyrobacter

A novel, Gram-negative, non-spore-forming, pleomorphic yellow-orange bacterial strain was isolated from a stadium seat. Strain Coronado^T falls within the Erythrobacteraceae family and the genus Porphyrobacter based on 16S rRNA phylogenetic analysis. This strain has Q-10 as the predominant respiratory lipoquinone, as do other members of the family. The fatty acid profile of this strain is similar to other Porphyrobacter, however Coronado^T contains predominately C18:1ω7cis and C16:0, a high percentage of the latter not being observed in any other Erythrobacteraceae. This strain is catalase-positive and oxidase-negative, can grow from 4 to 28 °C, at NaCl concentrations 0.1–1.5%, and at pH 6.0–8.0. On the basis of phenotypic and phylogenetic data presented in this study, strain Coronado^T represents a novel species in the Porphyrobacter genus for which the name Porphyrobacter mercurialis sp. nov. is proposed; the type strain is Coronado^T (=DSMZ 29971, =LMG 28700).

Draft Genome Sequence of Pseudoalteromonas tetraodonis Strain UCD-SED8 (Phylum Gammaproteobacteria)

Here, we present the draft genome sequence of Pseudoalteromonas tetraodonis UCD-SED8, a marine bacterium normally associated with the production of tetrodotoxin in pufferfish. This strain was isolated from sediment samples surrounding Zostera marina roots collected from Bodega Marine, California. The assembly consists of 4,017,727 bp contained in 35 contigs.

Phage on the stage

The resurgence of interest in bacteriophages for use in combating antibiotic resistant bacteria is coincident with an urgent call for more effective science education practices, including hands-on learning opportunities. To address this issue, a number of solutions have been proposed, including a large educational experiment, begun in 2007 by the Howard Hughes Medical Institute and currently involving over 85 colleges and universities, which has students discovering unique phages, obtaining images, and purifying phage DNA. A subset of these phage genomes is sequenced and analyzed using bioinformatics tools. Papers describing individual phage discoveries and comparative genomic studies are being published regularly. The vast majority of students in the program are in their first year of college, a critical time in capturing their interest and retaining them as science majors. This viral discovery model is being adopted and modified by a wide variety of educational institutions using a number of different bacterial hosts. In the opinion of the authors, this program and others like it represent a model accessible to virtually any undergraduate setting. And because of these programs, bacteriophage enthusiasts (academics, health professionals, biotechnology companies) can look forward to more well prepared students entering their ranks and should anticipate many more potentially useful phages discovered and characterized.

Draft Genome Sequence of Burkholderia gladioli Strain UCD-UG_CHAPALOTE (Phylum Proteobacteria)

Here, we present the draft genome of Burkholderia gladioli strain UCD-UG_CHAPALOTE. This strain is an endophyte isolated from surface sterilized seeds of an ancient Mexican landrace of corn, Chapalote. The genome contains 8,527,129 bp in 109 scaffolds.

Draft Genome Sequence of Enterobacter sp. Strain UCD-UG_FMILLET (Phylum Proteobacteria)

Here, we present the draft genome of Enterobacter sp. strain UCD-UG_FMILLET. This strain is an endophyte isolated from the roots of finger millet, an Afro-Indian cereal crop. The genome contains 4,801,411 bp in 53 scaffolds.