Community assembly of the human piercing microbiome.

Predicting how biological communities respond to disturbance requires understanding the forces that govern their assembly. We propose using human skin piercings as a model system for studying community assembly after rapid environmental change. Local skin sterilization provides a 'clean slate' within the novel ecological niche created by the piercing. Stochastic assembly processes can dominate skin microbiomes due to the influence of environmental exposure on local dispersal, but deterministic processes might play a greater role within occluded skin piercings if piercing habitats impose strong selection pressures on colonizing species. Here we explore the human ear-piercing microbiome and demonstrate that community assembly is predominantly stochastic but becomes significantly more deterministic with time, producing increasingly diverse and ecologically complex communities. We also observed changes in two dominant and medically relevant antagonists (Cutibacterium acnes and Staphylococcus epidermidis), consistent with competitive exclusion induced by a transition from sebaceous to moist environments. By exploiting this common yet uniquely human practice, we show that skin piercings are not just culturally significant but also represent ecosystem engineering on the human body. The novel habitats and communities that skin piercings produce may provide general insights into biological responses to environmental disturbances with implications for both ecosystem and human health.

A new multiplex qPCR assay to detect and differentiate big cat species in the illegal wildlife trade.

All species of big cats, including tigers, cheetahs, leopards, lions, snow leopards, and jaguars, are protected under the Convention on the International Trade in Endangered Species (CITES). This is due in large part to population declines resulting from anthropogenic factors, especially poaching and the unregulated and illegal trade in pelts, bones, teeth and other products that are derived from these iconic species. To enhance and scale up monitoring for big cat products in this trade, we created a rapid multiplex qPCR test that can identify and differentiate DNA from tiger (Panthera tigris), cheetah (Acinonyx jubatus), leopard (Panthera pardus), lion (Panthera leo), snow leopard (Panthera uncia), and jaguar (Panthera onca) in wildlife products using melt curve analysis to identify each species by its unique melt peak temperature. Our results showed high PCR efficiency (> 90%), sensitivity (detection limit of 5 copies of DNA per PCR reaction) and specificity (no cross amplification between each of the 6 big cat species). When paired with a rapid (< 1 h) DNA extraction protocol that amplifies DNA from bone, teeth, and preserved skin, total test time is less than three hours. This test can be used as a screening method to improve our understanding of the scale and scope of the illegal trade in big cats and aid in the enforcement of international regulations that govern the trade in wildlife and wildlife products, both ultimately benefiting the conservation of these species worldwide.

Estimating biodiversity across the tree of life on Mount Everest's southern flank with environmental DNA.

Species composition in high-alpine ecosystems is a useful indicator for monitoring climatic and environmental changes at the upper limits of habitable environments. We used environmental DNA (eDNA) analysis to document the breadth of high-alpine biodiversity present on Earth's highest mountain, Mt. Everest (8,849 m a.s.l.) in Nepal's Khumbu region. In April-May 2019, we collected eDNA from ten ponds and streams between 4,500 m and 5,500 m. Using multiple sequencing and bioinformatic approaches, we identified taxa from 36 phyla and 187 potential orders across the Tree of Life in Mt. Everest's high-alpine and aeolian ecosystem. These organisms, all recorded above 4,500 m-an elevational belt comprising <3% of Earth's land surface-represents ∼16% of global taxonomic order estimates. Our eDNA inventory will aid future high-Himalayan biomonitoring and retrospective molecular studies to assess changes over time as climate-driven warming, glacial melt, and anthropogenic influences reshape this rapidly transforming world-renowned ecosystem.

Measuring protected-area effectiveness using vertebrate distributions from leech iDNA.

Protected areas are key to meeting biodiversity conservation goals, but direct measures of effectiveness have proven difficult to obtain. We address this challenge by using environmental DNA from leech-ingested bloodmeals to estimate spatially-resolved vertebrate occupancies across the 677 km2 Ailaoshan reserve in Yunnan, China. From 30,468 leeches collected by 163 park rangers across 172 patrol areas, we identify 86 vertebrate species, including amphibians, mammals, birds and squamates. Multi-species occupancy modelling shows that species richness increases with elevation and distance to reserve edge. Most large mammals (e.g. sambar, black bear, serow, tufted deer) follow this pattern; the exceptions are the three domestic mammal species (cows, sheep, goats) and muntjak deer, which are more common at lower elevations. Vertebrate occupancies are a direct measure of conservation outcomes that can help guide protected-area management and improve the contributions that protected areas make towards global biodiversity goals. Here, we show the feasibility of using invertebrate-derived DNA to estimate spatially-resolved vertebrate occupancies across entire protected areas.

Promoting equity and inclusion with student-driven initiatives.

Representation is crucial in building more inclusive communities in science, technology, engineering, mathematics, and medicine (STEMM) fields. STEMM Diversity is a student-driven initiative that was founded to promote equity, diversity, and inclusion (EDI) at McGill University. Here, we discuss the lessons learned while developing STEMM Diversity that can help guide others to develop their own student-driven initiatives.

Transgenes of genetically modified animals detected non-invasively via environmental DNA.

We demonstrate that simple, non-invasive environmental DNA (eDNA) methods can detect transgenes of genetically modified (GM) animals from terrestrial and aquatic sources in invertebrate and vertebrate systems. We detected transgenic fragments between 82-234 bp through targeted PCR amplification of environmental DNA extracted from food media of GM fruit flies (Drosophila melanogaster), feces, urine, and saliva of GM laboratory mice (Mus musculus), and aquarium water of GM tetra fish (Gymnocorymbus ternetzi). With rapidly growing accessibility of genome-editing technologies such as CRISPR, the prevalence and diversity of GM animals will increase dramatically. GM animals have already been released into the wild with more releases planned in the future. eDNA methods have the potential to address the critical need for sensitive, accurate, and cost-effective detection and monitoring of GM animals and their transgenes in nature.

Resistance, resilience, and functional redundancy of freshwater bacterioplankton communities facing a gradient of agricultural stressors in a mesocosm experiment.

Agricultural pollution with fertilizers and pesticides is a common disturbance to freshwater biodiversity. Bacterioplankton communities are at the base of aquatic food webs, but their responses to these potentially interacting stressors are rarely explored. To test the extent of resistance and resilience in bacterioplankton communities faced with agricultural stressors, we exposed freshwater mesocosms to single and combined gradients of two commonly used pesticides: the herbicide glyphosate (0-15 mg/L) and the neonicotinoid insecticide imidacloprid (0-60 µg/L), in high or low nutrient backgrounds. Over the 43-day experiment, we tracked variation in bacterial density with flow cytometry, carbon substrate use with Biolog EcoPlates, and taxonomic diversity and composition with environmental 16S rRNA gene amplicon sequencing. We show that only glyphosate (at the highest dose, 15 mg/L), but not imidacloprid, nutrients, or their interactions measurably changed community structure, favouring members of the Proteobacteria including the genus Agrobacterium. However, no change in carbon substrate use was detected throughout, suggesting functional redundancy despite taxonomic changes. We further show that communities are resilient at broad, but not fine taxonomic levels: 24 days after glyphosate application the precise amplicon sequence variants do not return, and tend to be replaced by phylogenetically close taxa. We conclude that high doses of glyphosate - but still within commonly acceptable regulatory guidelines - alter freshwater bacterioplankton by favouring a subset of higher taxonomic units (i.e., genus to phylum) that transiently thrive in the presence of glyphosate. Longer-term impacts of glyphosate at finer taxonomic resolution merit further investigation.

Community rescue in experimental phytoplankton communities facing severe herbicide pollution.

Community rescue occurs when ecological or evolutionary processes restore positive growth in a highly stressful environment that was lethal to the community in its ancestral form, thus averting biomass collapse in a deteriorating environment. Laboratory evidence suggests that community rescue is most likely in high-biomass communities that have previously experienced moderate doses of sublethal stress. We assessed this result under more natural conditions, in a mesocosm experiment with phytoplankton communities exposed to the ubiquitous herbicide glyphosate. We tested whether community biomass and prior herbicide exposure would facilitate community rescue after severe contamination. We found that prior exposure to glyphosate was a very strong predictor of the rescue outcome, while high community biomass was not. Furthermore, although glyphosate had negative effects on diversity, it did not influence community composition significantly, suggesting a modest role for genus sorting in this rescue process. Our results expand the scope of community rescue theory to complex ecosystems and confirm that prior stress exposure is a key predictor of rescue.

Predator-induced collapse of niche structure and species coexistence.

Biological invasions are both a pressing environmental challenge and an opportunity to investigate fundamental ecological processes, such as the role of top predators in regulating biodiversity and food-web structure. In whole-ecosystem manipulations of small Caribbean islands on which brown anole lizards (Anolis sagrei) were the native top predator, we experimentally staged invasions by competitors (green anoles, Anolis smaragdinus) and/or new top predators (curly-tailed lizards, Leiocephalus carinatus). We show that curly-tailed lizards destabilized the coexistence of competing prey species, contrary to the classic idea of keystone predation. Fear-driven avoidance of predators collapsed the spatial and dietary niche structure that otherwise stabilized coexistence, which intensified interspecific competition within predator-free refuges and contributed to the extinction of green-anole populations on two islands. Moreover, whereas adding either green anoles or curly-tailed lizards lengthened food chains on the islands, adding both species reversed this effect-in part because the apex predators were trophic omnivores. Our results underscore the importance of top-down control in ecological communities, but show that its outcomes depend on prey behaviour, spatial structure, and omnivory. Diversity-enhancing effects of top predators cannot be assumed, and non-consumptive effects of predation risk may be a widespread constraint on species coexistence.

An efficient and robust laboratory workflow and tetrapod database for larger scale environmental DNA studies.

BACKGROUND: The use of environmental DNA for species detection via metabarcoding is growing rapidly. We present a co-designed lab workflow and bioinformatic pipeline to mitigate the 2 most important risks of environmental DNA use: sample contamination and taxonomic misassignment. These risks arise from the need for polymerase chain reaction (PCR) amplification to detect the trace amounts of DNA combined with the necessity of using short target regions due to DNA degradation. FINDINGS: Our high-throughput workflow minimizes these risks via a 4-step strategy: (i) technical replication with 2 PCR replicates and 2 extraction replicates; (ii) using multi-markers (12S,16S,CytB); (iii) a "twin-tagging," 2-step PCR protocol; and (iv) use of the probabilistic taxonomic assignment method PROTAX, which can account for incomplete reference databases. Because annotation errors in the reference sequences can result in taxonomic misassignment, we supply a protocol for curating sequence datasets. For some taxonomic groups and some markers, curation resulted in >50% of sequences being deleted from public reference databases, owing to (i) limited overlap between our target amplicon and reference sequences, (ii) mislabelling of reference sequences, and (iii) redundancy. Finally, we provide a bioinformatic pipeline to process amplicons and conduct PROTAX assignment and tested it on an invertebrate-derived DNA dataset from 1,532 leeches from Sabah, Malaysia. Twin-tagging allowed us to detect and exclude sequences with non-matching tags. The smallest DNA fragment (16S) amplified most frequently for all samples but was less powerful for discriminating at species rank. Using a stringent and lax acceptance criterion we found 162 (stringent) and 190 (lax) vertebrate detections of 95 (stringent) and 109 (lax) leech samples. CONCLUSIONS: Our metabarcoding workflow should help research groups increase the robustness of their results and therefore facilitate wider use of environmental and invertebrate-derived DNA, which is turning into a valuable source of ecological and conservation information on tetrapods.

One-locus-several-primers: A strategy to improve the taxonomic and haplotypic coverage in diet metabarcoding studies.

In diet metabarcoding analyses, insufficient taxonomic coverage of PCR primer sets generates false negatives that may dramatically distort biodiversity estimates. In this paper, we investigated the taxonomic coverage and complementarity of three cytochrome c oxidase subunit I gene (COI) primer sets based on in silico analyses and we conducted an in vivo evaluation using fecal and spider web samples from different invertivores, environments, and geographic locations. Our results underline the lack of predictability of both the coverage and complementarity of individual primer sets: (a) sharp discrepancies exist observed between in silico and in vivo analyses (to the detriment of in silico analyses); (b) both coverage and complementarity depend greatly on the predator and on the taxonomic level at which preys are considered; (c) primer sets' complementarity is the greatest at fine taxonomic levels (molecular operational taxonomic units [MOTUs] and variants). We then formalized the "one-locus-several-primer-sets" (OLSP) strategy, that is, the use of several primer sets that target the same locus (here the first part of the COI gene) and the same group of taxa (here invertebrates). The proximal aim of the OLSP strategy is to minimize false negatives by increasing total coverage through multiple primer sets. We illustrate that the OLSP strategy is especially relevant from this perspective since distinct variants within the same MOTUs were not equally detected across all primer sets. Furthermore, the OLSP strategy produces largely overlapping and comparable sequences, which cannot be achieved when targeting different loci. This facilitates the use of haplotypic diversity information contained within metabarcoding datasets, for example, for phylogeography and finer analyses of prey-predator interactions.

Linking a mutation to survival in wild mice.

Adaptive evolution in new or changing environments can be difficult to predict because the functional connections between genotype, phenotype, and fitness are complex. Here, we make these explicit connections by combining field and laboratory experiments in wild mice. We first directly estimate natural selection on pigmentation traits and an underlying pigment locus, Agouti, by using experimental enclosures of mice on different soil colors. Next, we show how a mutation in Agouti associated with survival causes lighter coat color through changes in its protein binding properties. Together, our findings demonstrate how a sequence variant alters phenotype and then reveal the ensuing ecological consequences that drive changes in population allele frequency, thereby illuminating the process of evolution by natural selection.

Keystone Genes.

The keystone species concept is used in ecology to describe individual species with disproportionately large effects on their communities. We extend this idea to the level of genes with disproportionately large effects on ecological processes. Such 'keystone genes' (KGs) would underlie traits involved in species interactions or causing critical biotic and/or abiotic changes that influence emergent community and ecosystem properties. We propose a general framework for how KGs could be identified, while keeping KGs under the umbrella of 'ecologically important genes' (EIGs) that also include categories such as 'foundation genes', 'ecosystem engineering genes', and more. Although likely rare, KGs and other EIGs could dominate certain ecological processes; thus, their discovery and study are relevant for understanding eco-evolutionary dynamics.

Carrion fly-derived DNA metabarcoding is an effective tool for mammal surveys: Evidence from a known tropical mammal community.

Metabarcoding of vertebrate DNA derived from carrion flies has been proposed as a promising tool for biodiversity monitoring. To evaluate its efficacy, we conducted metabarcoding surveys of carrion flies on Barro Colorado Island (BCI), Panama, which has a well-known mammal community, and compared our results against diurnal transect counts and camera trapping. We collected 1,084 flies in 29 sampling days, conducted metabarcoding with mammal-specific (16S) and vertebrate-specific (12S) primers, and sequenced amplicons on Illumina MiSeq. For taxonomic assignment, we compared blast with the new program protax, and we found that protax improved species identifications. We detected 20 mammal, four bird, and one lizard species from carrion fly metabarcoding, all but one of which are known from BCI. Fly metabarcoding detected more mammal species than concurrent transect counts (29 sampling days, 13 species) and concurrent camera trapping (84 sampling days, 17 species), and detected 67% of the number of mammal species documented by 8 years of transect counts and camera trapping combined, although fly metabarcoding missed several abundant species. This study demonstrates that carrion fly metabarcoding is a powerful tool for mammal biodiversity surveys and has the potential to detect a broader range of species than more commonly used methods.

Spider Web DNA: A New Spin on Noninvasive Genetics of Predator and Prey.

Noninvasive genetic sampling enables biomonitoring without the need to directly observe or disturb target organisms. This paper describes a novel and promising source of noninvasive spider and insect DNA from spider webs. Using black widow spiders (Latrodectus spp.) fed with house crickets (Acheta domesticus), we successfully extracted, amplified, and sequenced mitochondrial DNA from spider web samples that identified both spider and prey to species. Detectability of spider DNA did not differ between assays with amplicon sizes from 135 to 497 base pairs. Spider and prey DNA remained detectable at least 88 days after living organisms were no longer present on the web. Spider web DNA as a proof-of-concept may open doors to other practical applications in conservation research, pest management, biogeography studies, and biodiversity assessments.