Search

The human gut microbiome is promising therapeutic target, but development of interventions is hampered by limited understanding of the microbial ecosystem. Therefore, recent years have seen a surge in the engineering of inference algorithms seeking to unravel rules of ecological interactions from metagenomic data. Research groups score algorithmic performance in a variety of different ways, however, there exists no unified framework to score and rank each inference approach. The machine learning field presents a useful solution to this issue: a unified set of validation data and accompanying scoring metric. Here, we present MTIST: a platform for benchmarking microbial ecosystem inference tools. We use a generalized Lotka-Volterra framework to simulate microbial abundances over time, akin to what would be obtained by quantitative metagenomic sequencing studies or lab experiments, to generate a massivein silicotraining dataset (MTIST) for algorithmic validation, as well as an “ecological sign” score (ES score) to rate them. MTIST comprises 24,570 time series of microbial abundance data packaged into 648 datasets. Together, the MTIST dataset and the ES score serve as a platform to develop and compare microbiome ecosystem inference approaches.

The microbial populations in the gut microbiome have recently been associated with COVID-19 disease severity. However, a causal impact of the gut microbiome on COVID-19 patient health has not been established. Here we provide evidence that gut microbiome dysbiosis is associated with translocation of bacteria into the blood during COVID-19, causing life-threatening secondary infections. Antibiotics and other treatments during COVID-19 can potentially confound microbiome associations. We therefore first demonstrate that the gut microbiome is directly affected by SARS-CoV-2 infection in a dose-dependent manner in a mouse model, causally linking viral infection and gut microbiome dysbiosis. Comparison with stool samples collected from 97 COVID-19 patients at two different clinical sites also revealed substantial gut microbiome dysbiosis, paralleling our observations in the animal model. Specifically, we observed blooms of opportunistic pathogenic bacterial genera known to include antimicrobial-resistant species in hospitalized COVID-19 patients. Analysis of blood culture results testing for secondary microbial bloodstream infections with paired microbiome data obtained from these patients suggest that bacteria translocate from the gut into the systemic circulation of COVID-19 patients. These results are consistent with a direct role for gut microbiome dysbiosis in enabling dangerous secondary infections during COVID 19.

The microbial populations in the gut microbiome have recently been associated with COVID-19 disease severity. However, a causal impact of the gut microbiome on COVID-19 patient health has not been established. Here we provide evidence that gut microbiome dysbiosis is associated with translocation of bacteria into the blood during COVID-19, causing life-threatening secondary infections. Antibiotics and other treatments during COVID-19 can potentially confound microbiome associations. We therefore first demonstrate in a mouse model that SARS-CoV-2 infection can induce gut microbiome dysbiosis, which correlated with alterations to Paneth cells and goblet cells, and markers of barrier permeability. Comparison with stool samples collected from 96 COVID-19 patients at two different clinical sites also revealed substantial gut microbiome dysbiosis, paralleling our observations in the animal model. Specifically, we observed blooms of opportunistic pathogenic bacterial genera known to include antimicrobial-resistant species in hospitalized COVID-19 patients. Analysis of blood culture results testing for secondary microbial bloodstream infections with paired microbiome data obtained from these patients indicates that bacteria may translocate from the gut into the systemic circulation of COVID-19 patients. These results are consistent with a direct role for gut microbiome dysbiosis in enabling dangerous secondary infections during COVID-19.

Although protocols exist for the recovery of ancient DNA from land snail and marine bivalve shells, marine conch shells have yet to be studied from a paleogenomic perspective. We first present reference assemblies for both a 623.7 Mbp nuclear genome and a 15.4 kbp mitochondrial genome for Strombus pugilis, the West Indian fighting conch. We next detail a method to extract and sequence DNA from conch shells and apply it to conch from Bocas del Toro, Panama across three time periods: recently-eaten and discarded (n = 3), Late Holocene (984-1258 before present [BP]) archaeological midden (n = 5), and mid-Holocene (5711-7187 BP) paleontological fossil coral reef (n = 5). These results are compared to control DNA extracted from live-caught tissue and fresh shells (n = 5). Using high-throughput sequencing, we were able to obtain S. pugilis nuclear sequence reads from shells across all age periods: up to 92.5 thousand filtered reads per sample in live-caught shell material, 4.57 thousand for modern discarded shells, 12.1 thousand reads for archaeological shells, and 114 reads in paleontological shells. We confirmed authenticity of the ancient DNA recovered from the archaeological and paleontological shells based on 5.7× higher average frequency of deamination-driven misincorporations and 15% shorter average read lengths compared to the modern shells. Reads also mapped to the S. pugilis mitochondrial genome for all but the paleontological shells, with consistent ratios of mitochondrial to nuclear mapped reads across sample types. Our methods can be applied to diverse archaeological sites to facilitate reconstructions of the long-term impacts of human behaviour on mollusc evolutionary biology.

Although protocols exist for the recovery of ancient DNA from land snail and marine bivalve shells, marine conch shells have yet to be studied from a paleogenomic perspective. We first present reference assemblies for both a 623.7 Mbp nuclear genome and a 15.4 kbp mitochondrial genome for Strombus pugilis, the West Indian fighting conch. We next detail a method to extract and sequence DNA from conch shells and apply it to conch from Bocas del Toro, Panama across three time periods: recently-eaten and discarded (n=3), Late Holocene (984-1258 BP) archaeological midden (n=5), and a mid-Holocene (5711-7187 BP) paleontological fossil coral reef (n=5). These results are compared to control DNA extracted from live-caught tissue and fresh shells (n=5). Using high-throughput sequencing, we were able to obtain S. pugilis nuclear sequence reads from shells across all age periods: up to 92.5 thousand filtered reads per sample in live-caught shell material, 4.57 thousand for modern discarded shells, 12.1 thousand reads for archaeological shells, and 114 reads in paleontological shells. We confirmed authenticity of the ancient DNA recovered from the archaeological and paleontological shells based on 5.7x higher average frequency of deamination-driven misincorporations and 15% shorter average read lengths compared to the modern shells. Reads also mapped to the S. pugilis mitochondrial genome for all but the paleontological shells, with consistent ratios of mitochondrial to nuclear mapped reads across sample types. Our methods can be applied to diverse archaeological sites to facilitate reconstructions of the long-term impacts of human behavior on mollusc evolutionary biology.

High-quality genomic resources facilitate population-level and species-level comparisons to answer questions about behavioral ecology, morphological and physiological adaptations, as well as the evolution of genomic architecture. Squamate reptiles (lizards and snakes) are particularly diverse in characteristics that have intrigued evolutionary biologists, but high-quality genomic resources for squamates are relatively sparse. Lizards in the genus Sceloporus have a long history as important ecological, evolutionary, and physiological models, making them a valuable target for the development of genomic resources. We present a high-quality chromosome-level reference genome assembly, SceUnd1.0, (utilizing 10X Genomics Chromium, HiC, and PacBio data) and tissue/developmental stage transcriptomes for the Eastern Fence Lizard, Sceloporus undulatus. We performed synteny analysis with other available squamate chromosome-level assemblies to identify broad patterns of chromosome evolution including the fusion of micro- and macrochromosomes in S. undulatus. Using this new S. undulatus genome assembly we conducted reference-based assemblies for 34 other Sceloporus species to improve draft nuclear genomes assemblies from 1% coverage to 43% coverage on average. Across these species, typically >90% of reads mapped for species within 20 million years divergence from S. undulatus, this dropped to 75% reads mapped for species at 35 million years divergence. Finally we use RNAseq and whole genome resequencing data to compare the three assemblies as references, each representing an increased level of sequencing, cost and assembly efforts: Supernova Assembly with data from10X Genomics Chromium library; HiRise Assembly that added data from HiC library; and PBJelly Assembly that added data from PacBio sequencing. We found that the Supernova Assembly contained the full genome and was a suitable reference for RNAseq, but the chromosome-level scaffolds provided by the addition of the HiC data allowed the reference to be used for other whole genome analysis, including synteny and whole genome association mapping analyses. The addition of PacBio data provided negligible gains. Overall, these new genomic resources provide valuable tools for advanced molecular analysis of an organism that has become a model in physiology and evolutionary ecology.

The Holocene arrival of humans on Madagascar precipitated major changes to the island’s biodiversity. The now-extinct, endemic “subfossil” megafauna of Madagascar were likely hunted by the island’s early human inhabitants. Perhaps in part due to preferential hunting of larger prey, no surviving species on Madagascar is larger than 10 kg. Outside of Madagascar, size-selective hunting pressure has resulted in the phyletic dwarfism of many still-living species across a diversity of phyla. On Madagascar, some subfossil bones of extant lemurs are considerably larger than those of the modern members of their species, but relatively large distances between the subfossil localities and modern samples that have been compared to date makes it impossible to reject the possibility that these size differences more simply reflect pre-existing ecogeographic variation. Here, we used high-resolution 3D scan data to conduct comparative morphological analyses of subfossil and modern skeletal remains of one of the larger extant lemurs, Verreaux’s sifakas (Propithecus verreauxi) from subfossil and modern sites ∼10 km adjacent: Taolambiby (bones dated to 725-560 – 1075-955 cal. years before present) and Beza Mahafaly Special Reserve, respectively. We found that the average subfossil sifaka bone (n=12) is 9% and significantly larger than that of modern sifakas (n=31 individuals; permutation test; p=0.037). When restricting the analysis to the single element and side with the largest representation in the subfossil sample (n=4 right distal femora), the average subfossil bone is 10% larger (p=0.046). While we cannot yet conclude whether this size difference reflects evolutionary change or an archaeological aggregation/taphonomic process, if this is a case of phyletic dwarfism in response to human size-selective harvesting pressures then the estimated rate of evolutionary change is slightly higher than that previously calculated for other archaeological cases of this phenomenon.

We have identified a fixed nonsynonymous sequence difference between humans (Val381; derived variant) and Neandertals (Ala381; ancestral variant) in the ligand-binding domain of the aryl hydrocarbon receptor (AHR) gene. In an exome sequence analysis of four Neandertal and Denisovan individuals compared with nine modern humans, there are only 90 total nucleotide sites genome-wide for which archaic hominins are fixed for the ancestral nonsynonymous variant and the modern humans are fixed for the derived variant. Of those sites, only 27, including Val381 in the AHR, also have no reported variability in the human dbSNP database, further suggesting that this highly conserved functional variant is a rare event. Functional analysis of the amino acid variant Ala381 within the AHR carried by Neandertals and nonhuman primates indicate enhanced polycyclic aromatic hydrocarbon (PAH) binding, DNA binding capacity, and AHR mediated transcriptional activity compared with the human AHR. Also relative to human AHR, the Neandertal AHR exhibited 150-1000 times greater sensitivity to induction of Cyp1a1 and Cyp1b1 expression by PAHs (e.g., benzo(a)pyrene). The resulting CYP1A1/CYP1B1 enzymes are responsible for PAH first pass metabolism, which can result in the generation of toxic intermediates and perhaps AHR-associated toxicities. In contrast, the human AHR retains the ancestral sensitivity observed in primates to nontoxic endogenous AHR ligands (e.g., indole, indoxyl sulfate). Our findings reveal that a functionally significant change in the AHR occurred uniquely in humans, relative to other primates, that would attenuate the response to many environmental pollutants, including chemicals present in smoke from fire use during cooking.

Noninvasive sampling is an important development in population genetic monitoring of wild animals. Particularly, the collection of environmental DNA (eDNA) which can be collected without needing to encounter the target animal facilitates the genetic analysis of endangered species. One method that has been applied to these sample types is target capture and enrichment which overcomes the issue of high proportions of exogenous (nonhost) DNA from these lower quality samples. We tested whether target capture of mitochondrial DNA from sampled feeding traces of the aye‐aye, an endangered lemur species would yield mitochondrial DNA sequences for population genetic monitoring. We sampled gnawed wood where aye‐ayes excavate wood‐boring insect larvae from trees. We designed RNA probes complementary to the aye‐aye's mitochondrial genome and used these to isolate aye‐aye DNA from other nontarget DNA in these samples. We successfully retrieved six near‐complete mitochondrial genomes from two sites within the aye‐aye's geographic range that had not been sampled previously. Our method demonstrates the application of next‐generation molecular techniques to species of conservation concern. This method can likely be applied to alternative foraged remains to sample endangered species other than aye‐ayes.

Non-invasive sampling is an important development in population genetic monitoring of wild animals. Particularly, the collection of environmental DNA (eDNA) which can be collected without needing to encounter the target animal, facilitates the genetic analysis of cryptic and threatened species. One method that has been applied to these types of sample is target capture and enrichment which overcomes the issue of high proportions of exogenous (non-host) DNA from these lower quality samples. We tested whether target capture of mitochondrial DNA from sampled feeding traces of wild aye-ayes would yield mitochondrial DNA sequences for population genetic monitoring. We sampled gnawed wood from feeding traces where aye-ayes excavate wood-boring insect larvae from trees. We designed RNA probes complementary to the aye-aye’s mitochondrial genome and used these to isolate aye-aye DNA from other non-target DNA in these samples. We successfully retrieved six near-complete mitochondrial genomes from two sites within the aye-aye’s geographic range that had not been sampled previously. This method can likely be applied to alternative foraged remains to sample species other than aye-ayes. Our method demonstrates the application to next-generation molecular techniques to species of conservation concern.

A genome-wide association study (GWAS) identifies regions of the genome that likely affect the variable state of a phenotype of interest. These regions can then be studied with population genetic methods to make inferences about the evolutionary history of the trait. There are increasing opportunities to use GWAS results-even from clinically motivated studies-for tests of classic anthropological hypotheses. One such example, presented here as a case study for this approach, involves tooth development variation related to dental crowding. Specifically, more than 10% of humans fail to develop one or more permanent third molars (M3 agenesis). M3 presence/absence variation within human populations has a significant genetic component (heritability estimate h 2 = 0.47). The evolutionary significance of M3 agenesis has a long history of anthropological speculation. First, the modern frequency of M3 agenesis could reflect a relaxation of selection pressure to retain larger and more teeth following the origins of cooking and other food-softening behaviors (i.e., the genetic drift hypothesis or, classically, the "probable mutation effect"). Alternatively, commensurate with increasing hominin brain size and facial shortening, M3 agenesis may have conferred an adaptive fitness advantage if it reduced the risk of M3 impaction and potential health complications (i.e., the positive selection hypothesis). A recent GWAS identified 70 genetic loci that may play a role in human M3 presence/absence variation. To begin evaluating the contrasting evolutionary scenarios for M3 agenesis, we used the integrated haplotype score (iHS) statistic to test whether those 70 genetic regions are enriched for genomic signatures of recent positive selection. None of our findings are inconsistent with the null hypothesis of genetic drift to explain the high prevalence of human M3 agenesis. This result might suggest that M3 impaction rates for modern humans do not accurately retrodict those of the preagricultural past. Alternatively, the absence of support for the positive selection hypothesis could reflect a lack of power; this analysis should be repeated following the completion of more comprehensive GWAS analyses for human M3 agenesis.

Due to our intensive subsistence and habitat-modification strategies-including broad-spectrum harvesting and predation, widespread landscape burning, settlement construction, and translocation of other species-humans have major roles as ecological actors who influence fundamental trophic interactions. Here we review how the long-term history of human-environment interaction has shaped the evolutionary biology of diverse non-human, non-domesticated species. Clear examples of anthropogenic effects on non-human morphological evolution have been documented in modern studies of substantial changes to body size or other major traits in terrestrial and aquatic vertebrates, invertebrates, and plants in response to selective human harvesting, urbanized habitats, and human-mediated translocation. Meanwhile, archaeological records of harvested marine invertebrates and terrestrial vertebrates suggest that similar processes extend considerably into prehistory, perhaps to 50,000 yr BP or earlier. These results are consistent with palaeoenvironmental and other records that demonstrate long-term human habitat modification and intensive harvesting practices. Thus, while considerable attention has been focused on recent human impacts on 'natural' habitats, integrated evidence from modern biology and archaeology suggests a deep history of human entanglement with our ecosystems including substantial effects on the evolutionary biology of non-human taxa. The number and magnitude of such effects will probably increase given the continued intensification of anthropogenic activities and ecosystem impacts, including climate change and direct genetic modification.

The Eurasian sympatry of Neandertals and anatomically modern humans – beginning at least 45,000 years ago and lasting for more than 5,000 years – has long sparked anthropological interest into the factors that potentially contributed to Neandertal extinction. Among many different hypotheses, the “differential pathogen resistance” extinction model posits that Neandertals were disproportionately affected by exposure to novel infectious diseases that were transmitted during the period of spatiotemporal sympatry with modern humans. Comparisons of new archaic hominin paleogenome sequences with modern human genomes have confirmed a history of genetic admixture – and thus direct contact – between humans and Neandertals. Analyses of these data have also shown that Neandertal nuclear genome genetic diversity was likely considerably lower than that of the Eurasian anatomically modern humans with whom they came into contact, perhaps leaving Neandertal innate immune systems relatively more susceptible to novel pathogens. In this study, we compared levels of genetic diversity in genes for which genetic variation is hypothesized to benefit pathogen defense among Neandertals and African, European, and Asian modern humans, using available exome sequencing data (six chromosomes per population). We observed that Neandertals had only 31-39% as many nonsynonymous (amino acid changing) polymorphisms across 73 innate immune system genes compared to modern human populations. We also found that Neandertal genetic diversity was relatively low in an unbiased set of balancing selection candidate genes for primates – genes with the highest 1% genetic diversity genome-wide in non-human apes. In contrast, Neandertals had similar to higher levels of genetic diversity than humans in 13 major histocompatibility complex (MHC) genes. Thus, while Neandertals may have been relatively more susceptible to some novel pathogens and differential pathogen resistance could be considered as one potential contributing factor in their extinction, this model does have limitations.