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.

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.

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.

Deep segregation in the open ocean: Macaronesia as an evolutionary hotspot for low dispersal marine invertebrates.

Diversification and speciation of terrestrial organisms are anticipated in oceanic islands such as Macaronesia, a group of Atlantic islands that have remained unconnected to continental landmasses. Hitherto, the diversification of marine organisms in oceanic islands, especially those with low vagility, has received little direct empirical analysis using molecular markers. Here, we focus on such a case study, through applying a multilocus molecular approach to investigate the diversity and evolution of a group that lacks a planktonic larval stage, the isopod genus Dynamene, in Macaronesia and Northeast Atlantic. Sequences of two mitochondrial (cytochrome c oxidase subunit I and 16S rRNA) and two nuclear (18S rRNA and 28S rRNA) loci were obtained from specimens of Dynamene edwardsi (Lucas, 1849), Dynamene magnitorata Holdich, 1968 and Dynamene bidentata (Adams, 1800) collected along the Northeast Atlantic and Macaronesia. Although no major phylogeographic structure was detected in D. bidentata and D. magnitorata, from five to nine deeply divergent lineages were evident within D. edwardsi. The divergent lineages displayed genetic distances comparable to those found among established species of peracarids. D. edwardsi exhibits a long, rich and complex phylogeographic history in Macaronesia, where the geodynamics of the islands possibly associated with founder effects and subsequent lack of gene flow among populations confounds patterns based on geographic proximity of targeted populations. Our findings collectively suggest a much larger role of oceanic islands in the diversification of marine invertebrates than previously anticipated. The work provides insights into the origins and dynamics of ongoing geographic segregation and associated deep divergence among sister evolutionary lineages in Macaronesia.

Acidity promotes degradation of multi-species environmental DNA in lotic mesocosms.

Accurate quantification of biodiversity is fundamental to understanding ecosystem function and for environmental assessment. Molecular methods using environmental DNA (eDNA) offer a non-invasive, rapid, and cost-effective alternative to traditional biodiversity assessments, which require high levels of expertise. While eDNA analyses are increasingly being utilized, there remains considerable uncertainty regarding the dynamics of multispecies eDNA, especially in variable systems such as rivers. Here, we utilize four sets of upland stream mesocosms, across an acid-base gradient, to assess the temporal and environmental degradation of multispecies eDNA. Sampling included water column and biofilm sampling over time with eDNA quantified using qPCR. Our findings show that the persistence of lotic multispecies eDNA, sampled from water and biofilm, decays to non-detectable levels within 2 days and that acidic environments accelerate the degradation process. Collectively, the results provide the basis for a predictive framework for the relationship between lotic eDNA degradation dynamics in spatio-temporally dynamic river ecosystems.

Stepping stones to isolation: Impacts of a changing climate on the connectivity of fragmented fish populations.

In the marine environment, understanding the biophysical mechanisms that drive variability in larval dispersal and population connectivity is essential for estimating the potential impacts of climate change on the resilience and genetic structure of populations. Species whose populations are small, isolated and discontinuous in distribution will differ fundamentally in their response and resilience to environmental stress, compared with species that are broadly distributed, abundant and frequently exchange conspecifics. Here, we use an individual-based modelling approach, combined with a population genetics projection model, to consider the impacts of a warming climate on the population connectivity of two contrasting Antarctic fish species, Notothenia rossii and Champsocephalus gunnari. Focussing on the Scotia Sea region, sea surface temperatures are predicted to increase significantly by the end of the 21st century, resulting in reduced planktonic duration and increased egg and larval mortality. With shorter planktonic durations, the results of our study predict reduced dispersal of both species across the Scotia Sea, from Antarctic Peninsula sites to islands in the north and east, and increased dispersal among neighbouring sites, such as around the Antarctic Peninsula. Increased mortality modified the magnitude of population connectivity but had little effect on the overall patterns. Whilst the predicted changes in connectivity had little impact on the projected regional population genetic structure of N. rossii, which remained broadly genetically homogeneous within distances of ~1,500 km, the genetic isolation of C. gunnari populations in the northern Scotia Sea was predicted to increase with rising sea temperatures. Our study highlights the potential for increased isolation of island populations in a warming world, with implications for the resilience of populations and their ability to adapt to ongoing environmental change, a matter of high relevance to fisheries and ecosystem-level management.

Performance of amplicon and shotgun sequencing for accurate biomass estimation in invertebrate community samples.

New applications of DNA and RNA sequencing are expanding the field of biodiversity discovery and ecological monitoring, yet questions remain regarding precision and efficiency. Due to primer bias, the ability of metabarcoding to accurately depict biomass of different taxa from bulk communities remains unclear, while PCR-free whole mitochondrial genome (mitogenome) sequencing may provide a more reliable alternative. Here, we used a set of documented mock communities comprising 13 species of freshwater macroinvertebrates of estimated individual biomass, to compare the detection efficiency of COI metabarcoding (three different amplicons) and shotgun mitogenome sequencing. Additionally, we used individual COI barcoding and de novo mitochondrial genome sequencing, to provide reference sequences for OTU assignment and metagenome mapping (mitogenome skimming), respectively. We found that, even though both methods occasionally failed to recover very low abundance species, metabarcoding was less consistent, by failing to recover some species with higher abundances, probably due to primer bias. Shotgun sequencing results provided highly significant correlations between read number and biomass in all but one species. Conversely, the read-biomass relationships obtained from metabarcoding varied across amplicons. Specifically, we found significant relationships for eight of 13 (amplicons B1FR-450 bp, FF130R-130 bp) or four of 13 (amplicon FFFR, 658 bp) species. Combining the results of all three COI amplicons (multiamplicon approach) improved the read-biomass correlations for some of the species. Overall, mitogenomic sequencing yielded more informative predictions of biomass content from bulk macroinvertebrate communities than metabarcoding. However, for large-scale ecological studies, metabarcoding currently remains the most commonly used approach for diversity assessment.

Range-wide genomic data synthesis reveals transatlantic vicariance and secondary contact in Atlantic cod.

Recent advances in genetic and genomic analysis have greatly improved our understanding of spatial population structure in marine species. However, studies addressing phylogeographic patterns at oceanic spatial scales remain rare. In Atlantic cod (Gadus morhua), existing range-wide examinations suggest significant transatlantic divergence, although the fine-scale contemporary distribution of populations and potential for secondary contact are largely unresolved. Here, we explore transatlantic phylogeography in Atlantic cod using a data-synthesis approach, integrating multiple genome-wide single-nucleotide polymorphism (SNP) datasets representative of different regions to create a single range-wide dataset containing 1,494 individuals from 54 locations and genotyped at 796 common loci. Our analysis highlights significant transatlantic divergence and supports the hypothesis of westward post-glacial colonization of Greenland from the East Atlantic. Accordingly, our analysis suggests the presence of transatlantic secondary contact off eastern North America and supports existing perspectives on the phylogeographic history of Atlantic cod with an unprecedented combination of genetic and geographic resolution. Moreover, we demonstrate the utility of integrating distinct SNP databases of high comparability.

Plio-Pleistocene phylogeography of the Southeast Asian Blue Panchax killifish, Aplocheilus panchax.

The complex climatic and geological history of Southeast Asia has shaped this region's high biodiversity. In particular, sea level fluctuations associated with repeated glacial cycles during the Pleistocene both facilitated, and limited, connectivity between populations. In this study, we used data from two mitochondrial and three anonymous nuclear markers to determine whether a fresh/brackish water killifish, Aplocheilus panchax, Hamilton, 1822, could be used to further understand how climatic oscillations and associated sea level fluctuations have shaped the distribution of biota within this region, and whether such patterns show evidence of isolation within palaeodrainage basins. Our analyses revealed three major mitochondrial clades within A. panchax. The basal divergence of A. panchax mitochondrial lineages was approximately 3.5 Ma, whilst the subsequent divergence timings of these clades occurred early Pleistocene (~2.6 Ma), proceeding through the Pleistocene. Continuous phylogeographic analysis showed a clear west-east dispersal followed by rapid radiation across Southeast Asia. Individuals from Krabi, just north of the Isthmus of Kra, were more closely related to the Indian lineages, providing further evidence for a freshwater faunal disjunction at the Isthmus of Kra biogeographic barrier. Our results suggest that Sulawesi, across the Wallace Line, was colonised relatively recently (~30 ka). Nuclear DNA is less geographically structured, although Mantel tests indicated that nuclear genetic distances were correlated with geographic proximity. Overall, these results imply that recent gene flow, as opposed to historical isolation, has been the key factor determining patterns of nuclear genetic variation in A. panchax, however, some evidence of historical isolation is retained within the mitochondrial genome. Our study further validates the existence of a major biogeographic boundary at the Kra Isthmus, and also demonstrates the use of widely distributed fresh/brackishwater species in phylogeographic studies, and their ability to disperse across major marine barriers in relatively recent time periods.

An ?Aukward' Tale: A Genetic Approach to Discover the Whereabouts of the Last Great Auks.

One hundred and seventy-three years ago, the last two Great Auks, Pinguinusimpennis, ever reliably seen were killed. Their internal organs can be found in the collections of the Natural History Museum of Denmark, but the location of their skins has remained a mystery. In 1999, Great Auk expert Errol Fuller proposed a list of five potential candidate skins in museums around the world. Here we take a palaeogenomic approach to test which-if any-of Fuller's candidate skins likely belong to either of the two birds. Using mitochondrial genomes from the five candidate birds (housed in museums in Bremen, Brussels, Kiel, Los Angeles, and Oldenburg) and the organs of the last two known individuals, we partially solve the mystery that has been on Great Auk scholars' minds for generations and make new suggestions as to the whereabouts of the still-missing skin from these two birds.

Annual time-series analysis of aqueous eDNA reveals ecologically relevant dynamics of lake ecosystem biodiversity.

The use of environmental DNA (eDNA) in biodiversity assessments offers a step-change in sensitivity, throughput and simultaneous measures of ecosystem diversity and function. There remains, however, a need to examine eDNA persistence in the wild through simultaneous temporal measures of eDNA and biota. Here, we use metabarcoding of two markers of different lengths, derived from an annual time series of aqueous lake eDNA to examine temporal shifts in ecosystem biodiversity and in an ecologically important group of macroinvertebrates (Diptera: Chironomidae). The analyses allow different levels of detection and validation of taxon richness and community composition (ß-diversity) through time, with shorter eDNA fragments dominating the eDNA community. Comparisons between eDNA, community DNA, taxonomy and UK species abundance data further show significant relationships between diversity estimates derived across the disparate methodologies. Our results reveal the temporal dynamics of eDNA and validate the utility of eDNA metabarcoding for tracking seasonal diversity at the ecosystem scale.

Recommendations for developing and applying genetic tools to assess and manage biological invasions in marine ecosystems.

The European Union's Marine Strategy Framework Directive (MSFD) aims to adopt integrated ecosystem management approaches to achieve or maintain "Good Environmental Status" for marine waters, habitats and resources, including mitigation of the negative effects of non-indigenous species (NIS). The Directive further seeks to promote broadly standardized monitoring efforts and assessment of temporal trends in marine ecosystem condition, incorporating metrics describing the distribution and impacts of NIS. Accomplishing these goals will require application of advanced tools for NIS surveillance and risk assessment, particularly given known challenges associated with surveying and monitoring with traditional methods. In the past decade, a host of methods based on nucleic acids (DNA and RNA) analysis have been developed or advanced that promise to dramatically enhance capacity in assessing and managing NIS. However, ensuring that these rapidly evolving approaches remain accessible and responsive to the needs of resource managers remains a challenge. This paper provides recommendations for future development of these genetic tools for assessment and management of NIS in marine systems, within the context of the explicit requirements of the MSFD. Issues considered include technological innovation, methodological standardization, data sharing and collaboration, and the critical importance of shared foundational resources, particularly integrated taxonomic expertise. Though the recommendations offered here are not exhaustive, they provide a basis for future intentional (and international) collaborative development of a genetic toolkit for NIS research, capable of fulfilling the immediate and long term goals of marine ecosystem and resource conservation.

Does density influence relative growth performance of farm, wild and F1 hybrid Atlantic salmon in semi-natural and hatchery common garden conditions?

The conditions encountered by Atlantic salmon, Salmo salar L., in aquaculture are markedly different from the natural environment. Typically, farmed salmon experience much higher densities than wild individuals, and may therefore have adapted to living in high densities. Previous studies have demonstrated that farmed salmon typically outgrow wild salmon by large ratios in the hatchery, but these differences are much less pronounced in the wild. Such divergence in growth may be explained partly by the offspring of wild salmon experiencing higher stress and thus lower growth when compared under high-density farming conditions. Here, growth of farmed, wild and F1 hybrid salmon was studied at contrasting densities within a hatchery and semi-natural environment. Farmed salmon significantly outgrew hybrid and wild salmon in all treatments. Importantly, however, the reaction norms were similar across treatments for all groups. Thus, this study was unable to find evidence that the offspring of farmed salmon have adapted more readily to higher fish densities than wild salmon as a result of domestication. It is suggested that the substantially higher growth rate of farmed salmon observed in the hatchery compared with wild individuals may not solely be caused by differences in their ability to grow in high-density hatchery scenarios.

A common garden design reveals population-specific variability in potential impacts of hybridization between populations of farmed and wild Atlantic salmon, Salmo salar L.

Released individuals can have negative impacts on native populations through various mechanisms, including competition, disease transfer and introduction of maladapted gene complexes. Previous studies indicate that the level of farmed Atlantic salmon introgression in native populations is population specific. However, few studies have explored the potential role of population diversity or river characteristics, such as temperature, on the consequences of hybridization. We compared freshwater growth of multiple families derived from two farmed, five wild and two F1 hybrid salmon populations at three contrasting temperatures (7°C, 12°C and 16°C) in a common garden experiment. As expected, farmed salmon outgrew wild salmon at all temperatures, with hybrids displaying intermediate growth. However, differences in growth were population specific and some wild populations performed better than others relative to the hybrid and farmed populations at certain temperatures. Therefore, the competitive balance between farmed and wild salmon may depend both on the thermal profile of the river and on the genetic characteristics of the respective farmed and wild strains. While limited to F1 hybridization, this study shows the merits in adopting a more complex spatially resolved approach to risk management of local populations.

The biogeography of the atlantic salmon (Salmo salar) gut microbiome.

Although understood in many vertebrate systems, the natural diversity of host-associated microbiota has been little studied in teleosts. For migratory fishes, successful exploitation of multiple habitats may affect and be affected by the composition of the intestinal microbiome. We collected 96 Salmo salar from across the Atlantic encompassing both freshwater and marine phases. Dramatic differences between environmental and gut bacterial communities were observed. Furthermore, community composition was not significantly impacted by geography. Instead life-cycle stage strongly defined both the diversity and identity of microbial assemblages in the gut, with evidence for community destabilisation in migratory phases. Mycoplasmataceae phylotypes were abundantly recovered in all life-cycle stages. Patterns of Mycoplasmataceae phylotype recruitment to the intestinal microbial community among sites and life-cycle stages support a dual role for deterministic and stochastic processes in defining the composition of the S. salar gut microbiome.

Oceanography and life history predict contrasting genetic population structure in two Antarctic fish species.

Understanding the key drivers of population connectivity in the marine environment is essential for the effective management of natural resources. Although several different approaches to evaluating connectivity have been used, they are rarely integrated quantitatively. Here, we use a 'seascape genetics' approach, by combining oceanographic modelling and microsatellite analyses, to understand the dominant influences on the population genetic structure of two Antarctic fishes with contrasting life histories, Champsocephalus gunnari and Notothenia rossii. The close accord between the model projections and empirical genetic structure demonstrated that passive dispersal during the planktonic early life stages is the dominant influence on patterns and extent of genetic structuring in both species. The shorter planktonic phase of C. gunnari restricts direct transport of larvae between distant populations, leading to stronger regional differentiation. By contrast, geographic distance did not affect differentiation in N. rossii, whose longer larval period promotes long-distance dispersal. Interannual variability in oceanographic flows strongly influenced the projected genetic structure, suggesting that shifts in circulation patterns due to climate change are likely to impact future genetic connectivity and opportunities for local adaptation, resilience and recovery from perturbations. Further development of realistic climate models is required to fully assess such potential impacts.

Borneo and Indochina are major evolutionary hotspots for Southeast Asian biodiversity.

Tropical Southeast (SE) Asia harbors extraordinary species richness and in its entirety comprises four of the Earth's 34 biodiversity hotspots. Here, we examine the assembly of the SE Asian biota through time and space. We conduct meta-analyses of geological, climatic, and biological (including 61 phylogenetic) data sets to test which areas have been the sources of long-term biological diversity in SE Asia, particularly in the pre-Miocene, Miocene, and Plio-Pleistocene, and whether the respective biota have been dominated by in situ diversification, immigration and/or emigration, or equilibrium dynamics. We identify Borneo and Indochina, in particular, as major "evolutionary hotspots" for a diverse range of fauna and flora. Although most of the region's biodiversity is a result of both the accumulation of immigrants and in situ diversification, within-area diversification and subsequent emigration have been the predominant signals characterizing Indochina and Borneo's biota since at least the early Miocene. In contrast, colonization events are comparatively rare from younger volcanically active emergent islands such as Java, which show increased levels of immigration events. Few dispersal events were observed across the major biogeographic barrier of Wallace's Line. Accelerated efforts to conserve Borneo's flora and fauna in particular, currently housing the highest levels of SE Asian plant and mammal species richness, are critically required.