Holly M. Bik

PhyloSift: phylogenetic analysis of genomes and metagenomes

Like all organisms on the planet, environmental microbes are subject to the forces of molecular evolution. Metagenomic sequencing provides a means to access the DNA sequence of uncultured microbes. By combining DNA sequencing of microbial communities with evolutionary modeling and phylogenetic analysis we might obtain new insights into microbiology and also provide a basis for practical tools such as forensic pathogen detection. In this work we present an approach to leverage phylogenetic analysis of metagenomic sequence data to conduct several types of analysis. First, we present a method to conduct phylogeny-driven Bayesian hypothesis tests for the presence of an organism in a sample. Second, we present a means to compare community structure across a collection of many samples and develop direct associations between the abundance of certain organisms and sample metadata. Third, we apply new tools to analyze the phylogenetic diversity of microbial communities and again demonstrate how this can be associated to sample metadata. These analyses are implemented in an open source software pipeline called PhyloSift. As a pipeline, PhyloSift incorporates several other programs including LAST, HMMER, and pplacer to automate phylogenetic analysis of protein coding and RNA sequences in metagenomic datasets generated by modern sequencing platforms (e.g., Illumina, 454).

Sequencing our way towards understanding global eukaryotic biodiversity

Microscopic eukaryotes are abundant, diverse and fill critical ecological roles across every ecosystem on Earth, yet there is a well-recognized gap in understanding of their global biodiversity. Fundamental advances in DNA sequencing and bioinformatics now allow accurate en masse biodiversity assessments of microscopic eukaryotes from environmental samples. Despite a promising outlook, the field of eukaryotic marker gene surveys faces significant challenges: how to generate data that are most useful to the community, especially in the face of evolving sequencing technologies and bioinformatics pipelines, and how to incorporate an expanding number of target genes.

An Introduction to Social Media for Scientists

Online social media tools can be some of the most rewarding and informative resources for scientists—IF you know how to use them.

Metagenetic community analysis of microbial eukaryotes illuminates biogeographic patterns in deep-sea and shallow water sediments.

Microbial eukaryotes (nematodes, protists, fungi, etc., loosely referred to as meiofauna) are ubiquitous in marine sediments and probably play pivotal roles in maintaining ecosystem function. Although the deep-sea benthos represents one of the worlds largest habitats, we lack a firm understanding of the biodiversity and community interactions amongst meiobenthic organisms in this ecosystem. Within this vast environment, key questions concerning the historical genetic structure of species remain a mystery, yet have profound implications for our understanding of global biodiversity and how we perceive and mitigate the impact of environmental change and anthropogenic disturbance. Using a metagenetic approach, we present an assessment of microbial eukaryote communities across depth (shallow water to abyssal) and ocean basins (deep-sea Pacific and Atlantic). Within the 12 sites examined, our results suggest that some taxa can maintain eurybathic ranges and cosmopolitan deep-sea distributions, but the majority of species appear to be regionally restricted. For Operationally Clustered Taxonomic Units (OCTUs) reporting wide distributions, there appears to be a taxonomic bias towards a small subset of taxa in most phyla; such bias may be driven by specific life history traits amongst these organisms. In addition, low genetic divergence between geographically disparate deep-sea sites suggests either a shorter coalescence time between deep-sea regions or slower rates of evolution across this vast oceanic ecosystem. While high-throughput studies allow for broad assessment of genetic patterns across microbial eukaryote communities, intragenomic variation in rRNA gene copies and the patchy coverage of reference databases currently present substantial challenges for robust taxonomic interpretations of eukaryotic data sets.

Dramatic Shifts in Benthic Microbial Eukaryote Communities following the Deepwater Horizon Oil Spill

Benthic habitats harbour a significant (yet unexplored) diversity of microscopic eukaryote taxa, including metazoan phyla, protists, algae and fungi. These groups are thought to underpin ecosystem functioning across diverse marine environments. Coastal marine habitats in the Gulf of Mexico experienced visible, heavy impacts following the Deepwater Horizon oil spill in 2010, yet our scant knowledge of prior eukaryotic biodiversity has precluded a thorough assessment of this disturbance. Using a marker gene and morphological approach, we present an intensive evaluation of microbial eukaryote communities prior to and following oiling around heavily impacted shorelines. Our results show significant changes in community structure, with pre-spill assemblages of diverse Metazoa giving way to dominant fungal communities in post-spill sediments. Post-spill fungal taxa exhibit low richness and are characterized by an abundance of known hydrocarbon-degrading genera, compared to prior communities that contained smaller and more diverse fungal assemblages. Comparative taxonomic data from nematodes further suggests drastic impacts; while pre-spill samples exhibit high richness and evenness of genera, post-spill communities contain mainly predatory and scavenger taxa alongside an abundance of juveniles. Based on this community analysis, our data suggest considerable (hidden) initial impacts across Gulf beaches may be ongoing, despite the disappearance of visible surface oil in the region.

Microbiota of the indoor environment: a meta-analysis

BackgroundAs modern humans, we spend the majority of our time in indoor environments. Consequently, environmental exposure to microorganisms has important implications for human health, and a better understanding of the ecological drivers and processes that impact indoor microbial assemblages will be key for expanding our knowledge of the built environment. In the present investigation, we combined recent studies examining the microbiota of the built environment in order to identify unifying community patterns and the relative importance of indoor environmental factors. Ultimately, the present meta-analysis focused on studies of bacteria and archaea due to the limited number of high-throughput fungal studies from the indoor environment. We combined 16S ribosomal RNA (rRNA) gene datasets from 16 surveys of indoor environments conducted worldwide, additionally including 7 other studies representing putative environmental sources of microbial taxa (outdoor air, soil, and the human body).ResultsCombined analysis of subsets of studies that shared specific experimental protocols or indoor habitats revealed community patterns indicative of consistent source environments and environmental filtering. Additionally, we were able to identify several consistent sources for indoor microorganisms, particularly outdoor air and skin, mirroring what has been shown in individual studies. Technical variation across studies had a strong effect on comparisons of microbial community assemblages, with differences in experimental protocols limiting our ability to extensively explore the importance of, for example, sampling locality, building function and use, or environmental substrate in structuring indoor microbial communities.ConclusionsWe present a snapshot of an important scientific field in its early stages, where studies have tended to focus on heavy sampling in a few geographic areas. From the practical perspective, this endeavor reinforces the importance of negative “kit” controls in microbiome studies. From the perspective of understanding mechanistic processes in the built environment, this meta-analysis confirms that broad factors, such as geography and building type, structure indoor microbes. However, this exercise suggests that individual studies with common sampling techniques may be more appropriate to explore the relative importance of subtle indoor environmental factors on the indoor microbiome.

Low endemism, continued deep-shallow interchanges, and evidence for cosmopolitan distributions in free-living marine nematodes (order Enoplida).

BackgroundNematodes represent the most abundant benthic metazoa in one of the largest habitats on earth, the deep sea. Characterizing major patterns of biodiversity within this dominant group is a critical step towards understanding evolutionary patterns across this vast ecosystem. The present study has aimed to place deep-sea nematode species into a phylogenetic framework, investigate relationships between shallow water and deep-sea taxa, and elucidate phylogeographic patterns amongst the deep-sea fauna.ResultsMolecular data (18 S and 28 S rRNA) confirms a high diversity amongst deep-sea Enoplids. There is no evidence for endemic deep-sea lineages in Maximum Likelihood or Bayesian phylogenies, and Enoplids do not cluster according to depth or geographic location. Tree topologies suggest frequent interchanges between deep-sea and shallow water habitats, as well as a mixture of early radiations and more recently derived lineages amongst deep-sea taxa. This study also provides convincing evidence of cosmopolitan marine species, recovering a subset of Oncholaimid nematodes with identical gene sequences (18 S, 28 S and cox1) at trans-Atlantic sample sites.ConclusionsThe complex clade structures recovered within the Enoplida support a high global species richness for marine nematodes, with phylogeographic patterns suggesting the existence of closely related, globally distributed species complexes in the deep sea. True cosmopolitan species may additionally exist within this group, potentially driven by specific life history traits of Enoplids. Although this investigation aimed to intensively sample nematodes from the order Enoplida, specimens were only identified down to genus (at best) and our sampling regime focused on an infinitesimal small fraction of the deep-sea floor. Future nematode studies should incorporate an extended sample set covering a wide depth range (shelf, bathyal, and abyssal sites), utilize additional genetic loci (e.g. mtDNA) that are informative at the species level, and apply high-throughput sequencing methods to fully assay community diversity. Finally, further molecular studies are needed to determine whether phylogeographic patterns observed in Enoplids are common across other ubiquitous marine groups (e.g. Chromadorida, Monhysterida).

Environmental DNA metabarcoding: transforming how we survey animal and plant communities

The genomic revolution has fundamentally changed how we survey biodiversity on earth. High‐throughput sequencing (“HTS”) platforms now enable the rapid sequencing of DNA from diverse kinds of environmental samples (termed “environmental DNA” or “eDNA”). Coupling HTS with our ability to associate sequences from eDNA with a taxonomic name is called “eDNA metabarcoding” and offers a powerful molecular tool capable of noninvasively surveying species richness from many ecosystems. Here, we review the use of eDNA metabarcoding for surveying animal and plant richness, and the challenges in using eDNA approaches to estimate relative abundance. We highlight eDNA applications in freshwater, marine and terrestrial environments, and in this broad context, we distill what is known about the ability of different eDNA sample types to approximate richness in space and across time. We provide guiding questions for study design and discuss the eDNA metabarcoding workflow with a focus on primers and library preparation methods. We additionally discuss important criteria for consideration of bioinformatic filtering of data sets, with recommendations for increasing transparency. Finally, looking to the future, we discuss emerging applications of eDNA metabarcoding in ecology, conservation, invasion biology, biomonitoring, and how eDNA metabarcoding can empower citizen science and biodiversity education.

Intra-genomic variation in the ribosomal repeats of nematodes.

Ribosomal loci represent a major tool for investigating environmental diversity and community structure via high-throughput marker gene studies of eukaryotes (e.g. 18S rRNA). Since the estimation of species’ abundance is a major goal of environmental studies (by counting numbers of sequences), understanding the patterns of rRNA copy number across species will be critical for informing such high-throughput approaches. Such knowledge is critical, given that ribosomal RNA genes exist within multi-copy repeated arrays in a genome. Here we measured the repeat copy number for six nematode species by mapping the sequences from whole genome shotgun libraries against reference sequences for their rRNA repeat. This revealed a 6-fold variation in repeat copy number amongst taxa investigated, with levels of intragenomic variation ranging from 56 to 323 copies of the rRNA array. By applying the same approach to four C. elegans mutation accumulation lines propagated by repeated bottlenecking for an average of ~400 generations, we find on average a 2-fold increase in repeat copy number (rate of increase in rRNA estimated at 0.0285-0.3414 copies per generation), suggesting that rRNA repeat copy number is subject to selection. Within each Caenorhabditis species, the majority of intragenomic variation found across the rRNA repeat was observed within gene regions (18S, 28S, 5.8S), suggesting that such intragenomic variation is not a product of selection for rRNA coding function. We find that the dramatic variation in repeat copy number among these six nematode genomes would limit the use of rRNA in estimates of organismal abundance. In addition, the unique pattern of variation within a single genome was uncorrelated with patterns of divergence between species, reflecting a strong signature of natural selection for rRNA function. A better understanding of the factors that control or affect copy number in these arrays, as well as their rates and patterns of evolution, will be critical for informing estimates of global biodiversity.

Next-generation phenomics for the Tree of Life

The phenotype represents a critical interface between the genome and the environment in which organisms live and evolve. Phenotypic characters also are a rich source of biodiversity data for tree building, and they enable scientists to reconstruct the evolutionary history of organisms, including most fossil taxa, for which genetic data are unavailable. Therefore, phenotypic data are necessary for building a comprehensive Tree of Life. In contrast to recent advances in molecular sequencing, which has become faster and cheaper through recent technological advances, phenotypic data collection remains often prohibitively slow and expensive. The next-generation phenomics project is a collaborative, multidisciplinary effort to leverage advances in image analysis, crowdsourcing, and natural language processing to develop and implement novel approaches for discovering and scoring the phenome, the collection of phentotypic characters for a species. This research represents a new approach to data collection that has the potential to transform phylogenetics research and to enable rapid advances in constructing the Tree of Life. Our goal is to assemble large phenomic datasets built using new methods and to provide the public and scientific community with tools for phenomic data assembly that will enable rapid and automated study of phenotypes across the Tree of Life.