An integrated spatio-temporal view of riverine biodiversity using environmental DNA metabarcoding.

Anthropogenically forced changes in global freshwater biodiversity demand more efficient monitoring approaches. Consequently, environmental DNA (eDNA) analysis is enabling ecosystem-scale biodiversity assessment, yet the appropriate spatio-temporal resolution of robust biodiversity assessment remains ambiguous. Here, using intensive, spatio-temporal eDNA sampling across space (five rivers in Europe and North America, with an upper range of 20-35 km between samples), time (19 timepoints between 2017 and 2018) and environmental conditions (river flow, pH, conductivity, temperature and rainfall), we characterise the resolution at which information on diversity across the animal kingdom can be gathered from rivers using eDNA. In space, beta diversity was mainly dictated by turnover, on a scale of tens of kilometres, highlighting that diversity measures are not confounded by eDNA from upstream. Fish communities showed nested assemblages along some rivers, coinciding with habitat use. Across time, seasonal life history events, including salmon and eel migration, were detected. Finally, effects of environmental conditions were taxon-specific, reflecting habitat filtering of communities rather than effects on DNA molecules. We conclude that riverine eDNA metabarcoding can measure biodiversity at spatio-temporal scales relevant to species and community ecology, demonstrating its utility in delivering insights into river community ecology during a time of environmental change.

Comparative population genomics of manta rays has global implications for management.

Understanding population connectivity and genetic diversity is of fundamental importance to conservation. However, in globally threatened marine megafauna, challenges remain due to their elusive nature and wide-ranging distributions. As overexploitation continues to threaten biodiversity across the globe, such knowledge gaps compromise both the suitability and effectiveness of management actions. Here, we use a comparative framework to investigate genetic differentiation and diversity of manta rays, one of the most iconic yet vulnerable groups of elasmobranchs on the planet. Despite their recent divergence, we show how oceanic manta rays (Mobula birostris) display significantly higher heterozygosity than reef manta rays (Mobula alfredi) and that M. birostris populations display higher connectivity worldwide. Through inferring modes of colonization, we reveal how both contemporary and historical forces have likely influenced these patterns, with important implications for population management. Our findings highlight the potential for fisheries to disrupt population dynamics at both local and global scales and therefore have direct relevance for international conservation of marine species.

Managing human-mediated range shifts: understanding spatial, temporal and genetic variation in marine non-native species.

The use of molecular tools to manage natural resources is increasingly common. However, DNA-based methods are seldom used to understand the spatial and temporal dynamics of species' range shifts. This is important when managing range shifting species such as non-native species (NNS), which can have negative impacts on biotic communities. Here, we investigated the ascidian NNS Ciona robusta, Clavelina lepadiformis, Microcosmus squamiger and Styela plicata using a combined methodological approach. We first conducted non-molecular biodiversity surveys for these NNS along the South African coastline, and compared the results with historical surveys. We detected no consistent change in range size across species, with some displaying range stability and others showing range shifts. We then sequenced a section of cytochrome c oxidase subunit I (COI) from tissue samples and found genetic differences along the coastline but no change over recent times. Finally, we found that environmental DNA metabarcoding data showed broad congruence with both the biodiversity survey and the COI datasets, but failed to capture the complete incidence of all NNS. Overall, we demonstrated how a combined methodological approach can effectively detect spatial and temporal variation in genetic composition and range size, which is key for managing both thriving NNS and threatened species. This article is part of the theme issue 'Species' ranges in the face of changing environments (part I)'.

Domestication-induced reduction in eye size revealed in multiple common garden experiments: The case of Atlantic salmon (Salmo salar L.).

Domestication leads to changes in traits that are under directional selection in breeding programmes, though unintentional changes in nonproduction traits can also arise. In offspring of escaping fish and any hybrid progeny, such unintentionally altered traits may reduce fitness in the wild. Atlantic salmon breeding programmes were established in the early 1970s, resulting in genetic changes in multiple traits. However, the impact of domestication on eye size has not been studied. We measured body size corrected eye size in 4000 salmon from six common garden experiments conducted under artificial and natural conditions, in freshwater and saltwater environments, in two countries. Within these common gardens, offspring of domesticated and wild parents were crossed to produce 11 strains, with varying genetic backgrounds (wild, domesticated, F1 hybrids, F2 hybrids and backcrosses). Size-adjusted eye size was influenced by both genetic and environmental factors. Domesticated fish reared under artificial conditions had smaller adjusted eye size when compared to wild fish reared under identical conditions, in both the freshwater and marine environments, and in both Irish and Norwegian experiments. However, in parr that had been introduced into a river environment shortly after hatching and sampled at the end of their first summer, differences in adjusted eye size observed among genetic groups were of a reduced magnitude and were nonsignificant in 2-year-old sea migrating smolts sampled in the river immediately prior to sea entry. Collectively, our findings could suggest that where natural selection is present, individuals with reduced eye size are maladapted and consequently have reduced fitness, building on our understanding of the mechanisms that underlie a well-documented reduction in the fitness of the progeny of domesticated salmon, including hybrid progeny, in the wild.

Environmental DNA provides higher resolution assessment of riverine biodiversity and ecosystem function via spatio-temporal nestedness and turnover partitioning.

Rapidly assessing biodiversity is essential for environmental monitoring; however, traditional approaches are limited in the scope needed for most ecological systems. Environmental DNA (eDNA) based assessment offers enhanced scope for assessing biodiversity, while also increasing sampling efficiency and reducing processing time, compared to traditional methods. Here we investigated the effects of landuse and seasonality on headwater community richness and functional diversity, via spatio-temporal dynamics, using both eDNA and traditional sampling. We found that eDNA provided greater resolution in assessing biodiversity dynamics in time and space, compared to traditional sampling. Community richness was seasonally linked, peaking in spring and summer, with temporal turnover having a greater effect on community composition compared to localized nestedness. Overall, our assessment of ecosystem function shows that community formation is driven by regional resource availability, implying regional management requirements should be considered. Our findings show that eDNA based ecological assessment is a powerful, rapid and effective assessment strategy that enables complex spatio-temporal studies of community diversity and ecosystem function, previously infeasible using traditional methods.

Animals, protists and bacteria share marine biogeographic patterns.

Over millennia, ecological and evolutionary mechanisms have shaped macroecological patterns across the tree of life. Research describing these patterns at both regional and global scales has traditionally focused on the study of metazoan species. Consequently, there is a limited understanding of cross-phylum biogeographic structuring and an escalating need to understand the macroecology of both microscopic and macroscopic organisms. Here we used environmental DNA (eDNA) metabarcoding to explore the biodiversity of marine metazoans, protists and bacteria along an extensive and highly heterogeneous coastline. Our results showed remarkably consistent biogeographic structure across the kingdoms of life despite billions of years of evolution. Analyses investigating the drivers of these patterns for each taxonomic kingdom found that environmental conditions (such as temperature) and, to a lesser extent, anthropogenic stressors (such as fishing pressure and pollution) explained some of the observed variation. Additionally, metazoans displayed biogeographic patterns that suggested regional biotic homogenization. Against the backdrop of global pervasive anthropogenic environmental change, our work highlights the importance of considering multiple domains of life to understand the maintenance and drivers of biodiversity patterns across broad taxonomic, ecological and geographical scales.

Molecular characterization of a marine turtle tumor epizootic, profiling external, internal and postsurgical regrowth tumors.

Sea turtle populations are under threat from an epizootic tumor disease (animal epidemic) known as fibropapillomatosis. Fibropapillomatosis continues to spread geographically, with prevalence of the disease also growing at many longer-affected sites globally. However, we do not yet understand the precise environmental, mutational and viral events driving fibropapillomatosis tumor formation and progression.Here we perform transcriptomic and immunohistochemical profiling of five fibropapillomatosis tumor types: external new, established and postsurgical regrowth tumors, and internal lung and kidney tumors. We reveal that internal tumors are molecularly distinct from the more common external tumors. However, they have a small number of conserved potentially therapeutically targetable molecular vulnerabilities in common, such as the MAPK, Wnt, TGFß and TNF oncogenic signaling pathways. These conserved oncogenic drivers recapitulate remarkably well the core pan-cancer drivers responsible for human cancers. Fibropapillomatosis has been considered benign, but metastatic-related transcriptional signatures are strongly activated in kidney and established external tumors. Tumors in turtles with poor outcomes (died/euthanized) have genes associated with apoptosis and immune function suppressed, with these genes providing putative predictive biomarkers.Together, these results offer an improved understanding of fibropapillomatosis tumorigenesis and provide insights into the origins, inter-tumor relationships, and therapeutic treatment for this wildlife epizootic.

Disentangling the effects of sex, life history and genetic background in Atlantic salmon: growth, heart and liver under common garden conditions.

Livestock domestication has long been a part of agriculture, estimated to have first occurred approximately 10 000 years ago. Despite the plethora of traits studied, there is little understanding of the possible impacts domestication has had on internal organs, which are key determinants of survival. Moreover, the genetic basis of observed associated changes in artificial environments is still puzzling. Here we examine impacts of captivity on two organs in Atlantic salmon (Salar salar) that have been domesticated for approximately 50 years: heart and liver, in addition to growth. We studied multiple families of wild, domesticated, F1 and F2 hybrid, and backcrossed strains of S. salar in replicated common garden tanks during the freshwater and marine stages of development. Heart and liver weight were investigated, along with heart morphology metrics examined in just the wild, domesticated and F1 hybrid strains (heart height and width). Growth was positively linked with the proportion of the domesticated strain, and recombination in F2 hybrids (and the potential disruption of co-adapted gene complexes) did not influence growth. Despite the influence of domestication on growth, we found no evidence for domestication-driven divergence in heart or liver morphology. However, sexual dimorphism was detected in heart morphology, and after controlling for body size, females exhibited significantly larger heart weight and heart width when compared with males. Wild females also had an increased heart height when compared with wild males, and this was not observed in any other strain. Females sampled in saltwater showed significantly larger heart height with rounder hearts, than saltwater males. Collectively, these results demonstrate an additive basis of growth and, despite a strong influence of domestication on growth, no clear evidence of changes in heart or liver morphology associated with domestication was identified.

Phylogenomics and species delimitation for effective conservation of manta and devil rays.

Practical biodiversity conservation relies on delineation of biologically meaningful units. Manta and devil rays (Mobulidae) are threatened worldwide, yet morphological similarities and a succession of recent taxonomic changes impede the development of an effective conservation strategy. Here, we generate genome-wide single nucleotide polymorphism (SNP) data from a geographically and taxonomically representative set of manta and devil ray samples to reconstruct phylogenetic relationships and evaluate species boundaries under the general lineage concept. We show that nominal species units supported by alternative data sources constitute independently evolving lineages, and find robust evidence for a putative new species of manta ray in the Gulf of Mexico. Additionally, we uncover substantial incomplete lineage sorting indicating that rapid speciation together with standing variation in ancestral populations has driven phylogenetic uncertainty within Mobulidae. Finally, we detect cryptic diversity in geographically distinct populations, demonstrating that management below the species level may be warranted in certain species. Overall, our study provides a framework for molecular genetic species delimitation that is relevant to wide-ranging taxa of conservation concern, and highlights the potential for genomic data to support effective management, conservation and law enforcement strategies.

Executing multi-taxa eDNA ecological assessment via traditional metrics and interactive networks.

Current approaches to ecological assessment are limited by the traditional morpho-taxonomic methods presently employed and the inability to meet increasing demands for rapid assessments. Advancements in high throughput sequencing now enable rapid high-resolution ecological assessment using environmental DNA (eDNA). Here we test the ability of using eDNA-based ecological assessment methods against traditional assessment of two key indicator groups (diatoms and macroinvertebrates) and show how eDNA across multiple gene regions (COI, rbcL, 12S and 18S) can be used to infer interactive networks that link to ecological assessment criteria. We compared results between taxonomic and eDNA based assessments and found significant positive associations between macroinvertebrate (p < 0.001 R2 = 0.645) and diatom (p = 0.015, R2 = 0.222) assessment metrics. We further assessed the ability of eDNA based assessment to identify environmentally sensitive genera and found an order of magnitude greater potential for 18S, versus COI or rbcL, to determine environmental filtering of ecologically assessed communities. Lastly, we compared the ability of traditional metrics against co-occurrence network properties of our combined 18S, COI and rbcL indicator genera to infer habitat quality measures currently used by managers. We found that transitivity (network connectivity), linkage density and cohesion were significantly associated with habitat modification scores (HMS), whereas network properties were inconsistent with linking to the habitat quality score (HQS) metric. The incorporation of multi-marker eDNA network assessment opens up a means for finer scale ecological assessment, currently limited using traditional methods. While utilization of eDNA-based assessment is recommended, direct comparisons with traditional approaches are difficult as the methods are intrinsically different and should be treated as such with regards to future research. Overall, our findings show that eDNA can be used for effective ecological assessment while offering a wider range of scope and application compared to traditional assessment methods.

Epistatic regulation of growth in Atlantic salmon revealed: a QTL study performed on the domesticated-wild interface.

BACKGROUND: Quantitative traits are typically considered to be under additive genetic control. Although there are indications that non-additive factors have the potential to contribute to trait variation, experimental demonstration remains scarce. Here, we investigated the genetic basis of growth in Atlantic salmon by exploiting the high level of genetic diversity and trait expression among domesticated, hybrid and wild populations. RESULTS: After rearing fish in common-garden experiments under aquaculture conditions, we performed a variance component analysis in four mapping populations totaling ~ 7000 individuals from six wild, two domesticated and three F1 wild/domesticated hybrid strains. Across the four independent datasets, genome-wide significant quantitative trait loci (QTLs) associated with weight and length were detected on a total of 18 chromosomes, reflecting the polygenic nature of growth. Significant QTLs correlated with both length and weight were detected on chromosomes 2, 6 and 9 in multiple datasets. Significantly, epistatic QTLs were detected in all datasets. DISCUSSION: The observed interactions demonstrated that the phenotypic effect of inheriting an allele deviated between half-sib families. Gene-by-gene interactions were also suggested, where the combined effect of two loci resulted in a genetic effect upon phenotypic variance, while no genetic effect was detected when the two loci were considered separately. To our knowledge, this is the first documentation of epistasis in a quantitative trait in Atlantic salmon. These novel results are of relevance for breeding programs, and for predicting the evolutionary consequences of domestication-introgression in wild populations.

Publisher Correction: Detection of introduced and resident marine species using environmental DNA metabarcoding of sediment and water.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Comparative genetic stock structure in three species of commercially exploited Indo-Malay Carangidae (Teleosteii, Perciformes).

We examine genetic structuring in three commercially important species of the teleost family Carangidae from Malaysian waters: yellowtail scad Atule mate, bigeye scad Selar crumenophthalmus and yellowstripe scad Selaroides leptolepis, from the Indo-Malay Archipelago. In view of their distribution across contrasting habitats, we tested the hypothesis that pelagic species display less genetic divergence compared with demersal species, due to their potential to undertake long-distance migrations in oceanic waters. To evaluate population genetic structure, we sequenced two mitochondrial (mt)DNA [650 bp of cytochrome oxidase I (coI), 450 bp of control region (CR)] and one nuclear gene (910 bp of rag1) in each species. One hundred and eighty samples from four geographical regions within the Indo-Malay Archipelago including a population of yellowtail from Kuwait were examined. Findings revealed that the extent of genetic structuring among populations in the semi-pelagic and pelagic, yellowtail and bigeye were lower than demersal yellowstripe, consistent with the hypothesis that pelagic species display less genetic divergence compared with demersal species. The yellowtail phylogeny identified three distinct clades with bootstrap values of 86%-99% in mtDNA and 63%-67% in rag1. However, in bigeye, three clades were also observed from mtDNA data while only one clade was identified in rag1 dataset. In yellowstripe, the mtDNA tree was split into three closely related clades and two clades in rag1 tree with bootstraps value of 73%-99% and 56% respectively. However, no geographic structure appears in both mtDNA and rag1 datasets. Hierarchical molecular variance analysis (AMOVA), pair wise FST comparisons and the nearest-neighbour statistic (Snn ) showed significant genetic differences among Kuwait and Indo-Malay yellowtail. Within the Indo-Malay Archipelago itself, two distinct mitochondrial lineages were detected in yellowtail suggesting potential cryptic species. Findings suggests varying degrees of genetic structuring, key information relevant to management of exploited stocks, though more rapidly evolving genetic markers should be used in future to better delimit the nature and dynamics of putative stock boundaries.

Demographic reconstruction from ancient DNA supports rapid extinction of the great auk.

The great auk was once abundant and distributed across the North Atlantic. It is now extinct, having been heavily exploited for its eggs, meat, and feathers. We investigated the impact of human hunting on its demise by integrating genetic data, GPS-based ocean current data, and analyses of population viability. We sequenced complete mitochondrial genomes of 41 individuals from across the species' geographic range and reconstructed population structure and population dynamics throughout the Holocene. Taken together, our data do not provide any evidence that great auks were at risk of extinction prior to the onset of intensive human hunting in the early 16th century. In addition, our population viability analyses reveal that even if the great auk had not been under threat by environmental change, human hunting alone could have been sufficient to cause its extinction. Our results emphasise the vulnerability of even abundant and widespread species to intense and localised exploitation.

Author Correction: Gene-associated markers provide tools for tackling illegal fishing and false eco-certification.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Environmental DNA size sorting and degradation experiment indicates the state of Daphnia magna mitochondrial and nuclear eDNA is subcellular.

Environmental DNA analysis has emerged as a key component of biodiversity and environmental monitoring. However, the state and fate of eDNA in natural environments is still poorly understood for many ecological systems. Here we assess the state and fate of eDNA derived from the water flea, Daphnia magna, using a full factorial mesocosm experiment. We measured the quantity and degradation of eDNA over a two month period across a range of filters differing in pore size (0, 0.2, 1 and 10 µm), which spans the range of eDNA source material including subcellular, cellular and tissue. We also used two primer sets targeting mitochondrial (COI) and nuclear (18S) genomic regions. Our findings demonstrated that eDNA was most prevalent in the effluent water, but also reliably detected on the 0.2 µm filter, suggesting subcellular material is the predominate state of eDNA. Temporal eDNA quantity dynamics followed an exponential decay function over the course of 6-17 days, demonstrating a predictable decline in eDNA concentration. Nuclear eDNA was more abundant than mitochondrial eDNA, which may be a result of greater primer affinity, or indicate greater availability of nuclear eDNA gene targets in the environment. In contrast to two previous size-sorting experiments, which utilizing fish eDNA, our findings suggest that the state of invertebrate eDNA is much smaller than previously suspected. Overall, our data suggest that the detection of eDNA greatly depends on our knowledge of the state and fate of eDNA, which differ among species, and likely across environmental conditions.

Detection of introduced and resident marine species using environmental DNA metabarcoding of sediment and water.

Environmental DNA (eDNA) surveys are increasingly being used for biodiversity monitoring, principally because they are sensitive and can provide high resolution community composition data. Despite considerable progress in recent years, eDNA studies examining how different environmental sample types can affect species detectability remain rare. Comparisons of environmental samples are especially important for providing best practice guidance on early detection and subsequent mitigation of non-indigenous species. Here we used eDNA metabarcoding of COI (cytochrome c oxidase subunit I) and 18S (nuclear small subunit ribosomal DNA) genes to compare community composition between sediment and water samples in artificial coastal sites across the United Kingdom. We first detected markedly different communities and a consistently greater number of distinct operational taxonomic units in sediment compared to water. We then compared our eDNA datasets with previously published rapid assessment biodiversity surveys and found excellent concordance among the different survey techniques. Finally, our eDNA surveys detected many non-indigenous species, including several newly introduced species, highlighting the utility of eDNA metabarcoding for both early detection and temporal / spatial monitoring of non-indigenous species. We conclude that careful consideration on environmental sample type is needed when conducting eDNA surveys, especially for studies assessing community change.

Evolutionary drivers of kype size in Atlantic salmon (Salmo salar): domestication, age and genetics.

The diversity of reproduction and associated mating patterns in Atlantic salmon (Salmo salar) has long captivated evolutionary biologists. Salmo salar exhibit strategies involving migration, bold mating behaviours and radical morphological and physiological change. One such radical change is the elongation and curvature of the lower jaw in sexually mature males into a hook-like appendage called the kype. The kype is a secondary sexual characteristic used in mating hierarchies and a prime candidate for sexual selection. As one of the core global aquaculture fish species, however, mate choice, and thus sexual selection, has been replaced by industrial artificial fertilization seeking to develop more commercially viable strains. Removal of mate choice provides a unique opportunity to examine the kype over successive generations in the absence of sexual selection. Here we use a large-scale common-garden experiment, incorporating six experimental strains (wild, farmed and wild × farmed hybrids), experiencing one to three sea winters, to assess the impact of age and genetic background. After controlling for allometry, fork length-adjusted kype height (AKH) was significantly reduced in the domesticated strain in comparison to two wild strains. Furthermore, genetic variation at a locus on linkage group SSA1 was associated with kype height, and a locus on linkage group SSA23 was associated with fork length-adjusted kype length (AKL). The reduction in fork length-AKH in domesticated salmon suggests that the kype is of importance in mate choice and that it has decreased due to relaxation of sexual selection. Fork length-AKL showed an increase in domesticated individuals, highlighting that it may not be an important cue in mate choice. These results give us insight into the evolutionary significance of the kype, as well as implications of genetic induced phenotypic change caused by domesticated individuals escaping into the natural environment.