Range-wide fragmentation in a threatened fish associated with post-European settlement modification in the Murray–Darling Basin, Australia

Longitudinal monitoring of neutral and adaptive genomic diversity in a reintroduction

Restoration programs in the form of ex-situ breeding combined with reintroductions are becoming critical to counteract demographic declines and species losses. Such programs are increasingly using genetic management to improve conservation outcomes. However, the lack of long-term monitoring of genetic indicators following reintroduction prevents assessments of the trajectory and persistence of reintroduced populations. We carried out an extensive monitoring program in the wild for a threatened small-bodied fish (southern pygmy perch, Nannoperca australis) to assess the long-term genomic effects of its captive breeding and reintroduction. The species was rescued prior to its extirpation from the terminal lakes of Australia's Murray-Darling Basin, and then used for genetically informed captive breeding and reintroductions. Subsequent annual or biannual monitoring of abundance, fitness, and occupancy over a period of 11 years, combined with postreintroduction genetic sampling, revealed survival and recruitment of reintroduced fish. Genomic analyses based on data from the original wild rescued, captive born, and reintroduced cohorts revealed low inbreeding and strong maintenance of neutral and candidate adaptive genomic diversity across multiple generations. An increasing trend in the effective population size of the reintroduced population was consistent with field monitoring data in demonstrating successful re-establishment of the species. This provides a rare empirical example that the adaptive potential of a locally extinct population can be maintained during genetically informed ex-situ conservation breeding and reintroduction into the wild. Strategies to improve biodiversity restoration via ex-situ conservation should include genetic-based captive breeding and longitudinal monitoring of standing genomic variation in reintroduced populations.

Variation in intraspecific demography drives localised concordance but species-wide discordance in response to past climatic change

BACKGROUND

Understanding how species biology may facilitate resilience to climate change remains a critical factor in detecting and protecting species at risk of extinction. Many studies have focused on the role of particular ecological traits in driving species responses, but less so on demographic history and levels of standing genetic variation. Additionally, spatial variation in the interaction of demographic and adaptive factors may further complicate prediction of species responses to environmental change. We used environmental and genomic datasets to reconstruct the phylogeographic histories of two ecologically similar and largely co-distributed freshwater fishes, the southern (Nannoperca australis) and Yarra (N. obscura) pygmy perches, to assess the degree of concordance in their responses to Plio-Pleistocene climatic changes. We described contemporary genetic diversity, phylogenetic histories, demographic histories, and historical species distributions across both species, and statistically evaluated the degree of concordance in co-occurring populations.

RESULTS

Marked differences in contemporary genetic diversity, historical distribution changes and historical migration were observed across the species, with a distinct lack of genetic diversity and historical range expansion suggested for N. obscura. Although several co-occurring populations within a shared climatic refugium demonstrated concordant demographic histories, idiosyncratic population size changes were found at the range edges of the more spatially restricted species. Discordant responses between species were associated with low standing genetic variation in peripheral populations. This might have hindered adaptive potential, as documented in recent demographic declines and population extinctions for the two species.

CONCLUSION

Our results highlight both the role of spatial scale in the degree of concordance in species responses to climate change, and the importance of standing genetic variation in facilitating range shifts. Even when ecological traits are similar between species, long-term genetic diversity and historical population demography may lead to discordant responses to ongoing and future climate change.

The roles of aridification and sea level changes in the diversification and persistence of freshwater fish lineages

While the influence of Pleistocene climatic changes on divergence and speciation has been well-documented across the globe, complex spatial interactions between hydrology and eustatics over longer timeframes may also determine species evolutionary trajectories. Within the Australian continent, glacial cycles were not associated with changes in ice cover and instead largely resulted in fluctuations from moist to arid conditions across the landscape. We investigated the role of hydrological and coastal topographic changes brought about by Plio-Pleistocene climatic changes on the biogeographic history of a small Australian freshwater fish, the southern pygmy perch Nannoperca australis. Using 7958 ddRAD-seq (double digest restriction-site associated DNA) loci and 45,104 filtered SNPs, we combined phylogenetic, coalescent and species distribution analyses to assess the various roles of aridification, sea level and tectonics and associated biogeographic changes across southeast Australia. Sea-level changes since the Pliocene and reduction or disappearance of large waterbodies throughout the Pleistocene were determining factors in strong divergence across the clade, including the initial formation and maintenance of a cryptic species, N. 'flindersi'. Isolated climatic refugia and fragmentation due to lack of connected waterways maintained the identity and divergence of inter- and intraspecific lineages. Our historical findings suggest that predicted increases in aridification and sea level due to anthropogenic climate change might result in markedly different demographic impacts, both spatially and across different landscape types.

Recent and rapid anthropogenic habitat fragmentation increases extinction risk for freshwater biodiversity

Anthropogenic habitat fragmentation is often implicated as driving the current global extinction crisis, particularly in freshwater ecosystems. The genetic signal of recent population isolation can be confounded by the complex spatial arrangement of dendritic river systems. Consequently, many populations may presently be managed separately based on an incorrect assumption that they have evolved in isolation. Integrating landscape genomics data with models of connectivity that account for landscape structure, we show that the cumulative effects of multiple in-stream barriers have contributed to the recent decline of a freshwater fish from the Murray-Darling Basin, Australia. In addition, individual-based eco-evolutionary simulations further demonstrate that contemporary inferences about population isolation are consistent with the 160-year time frame since construction of in-stream barriers began in the region. Our findings suggest that the impact of very recent fragmentation may be often underestimated for freshwater biodiversity. We argue that proactive conservation measures to reconnect many riverine populations are urgently needed.

Different landscape effects on the genetic structure of two broadly distributed woody legumes, Acacia salicina and A. stenophylla (Fabaceae)

Restoring degraded landscapes has primarily focused on re-establishing native plant communities. However, little is known with respect to the diversity and distribution of most key revegetation species or the environmental and anthropogenic factors that may affect their demography and genetic structure. In this study, we investigated the genetic structure of two widespread Australian legume species (Acacia salicina and Acacia stenophylla) in the Murray-Darling Basin (MDB), a large agriculturally utilized region in Australia, and assessed the impact of landscape structure on genetic differentiation. We used AFLP genetic data and sampled a total of 28 A. salicina and 30 A. stenophylla sampling locations across southeastern Australia. We specifically evaluated the importance of four landscape features: forest cover, land cover, water stream cover, and elevation. We found that both species had high genetic diversity (mean percentage of polymorphic loci, 55.1% for A. salicina versus. 64.3% for A. stenophylla) and differentiation among local sampling locations (A. salicina: ΦPT = 0.301, 30%; A. stenophylla: ΦPT = 0.235, 23%). Population structure analysis showed that both species had high levels of structure (6 clusters each) and admixture in some sampling locations, particularly A. stenophylla. Although both species have a similar geographic range, the drivers of genetic connectivity for each species were very different. Genetic variation in A. salicina seems to be mainly driven by geographic distance, while for A. stenophylla, land cover appears to be the most important factor. This suggests that for the latter species, gene flow among populations is affected by habitat fragmentation. We conclude that these largely co-occurring species require different management actions to maintain population connectivity. We recommend active management of A. stenophylla in the MDB to improve gene flow in the adversity of increasing disturbances (e.g., droughts) driven by climate change and anthropogenic factors.

Phylogenomic history of enigmatic pygmy perches: implications for biogeography, taxonomy and conservation

Pygmy perches (Percichthyidae) are a group of poorly dispersing freshwater fishes that have a puzzling biogeographic disjunction across southern Australia. Current understanding of pygmy perch phylogenetic relationships suggests past east-west migrations across a vast expanse of now arid habitat in central southern Australia, a region lacking contemporary rivers. Pygmy perches also represent a threatened group with confusing taxonomy and potentially cryptic species diversity. Here, we present the first study of the evolutionary history of pygmy perches based on genome-wide information. Data from 13 991 ddRAD loci and a concatenated sequence of 1 075 734 bp were generated for all currently described and potentially cryptic species. Phylogenetic relationships, biogeographic history and cryptic diversification were inferred using a framework that combines phylogenomics, species delimitation and estimation of divergence times. The genome-wide phylogeny clarified the biogeographic history of pygmy perches, demonstrating multiple east-west events of divergence within the group across the Australian continent. These results also resolved discordance between nuclear and mitochondrial data from a previous study. In addition, we propose three cryptic species within a southwestern species complex. The finding of potentially new species demonstrates that pygmy perches may be even more susceptible to ecological and demographic threats than previously thought. Our results have substantial implications for improving conservation legislation of pygmy perch lineages, especially in southwestern Western Australia.

De novo genome assembly and annotation of Australia's largest freshwater fish, the Murray cod (Maccullochella peelii), from Illumina and Nanopore sequencing read

One of the most iconic Australian fish is the Murray cod, Maccullochella peelii (Mitchell 1838), a freshwater species that can grow to ∼1.8 metres in length and live to age ≥48 years. The Murray cod is of a conservation concern as a result of strong population contractions, but it is also popular for recreational fishing and is of growing aquaculture interest. In this study, we report the whole genome sequence of the Murray cod to support ongoing population genetics, conservation, and management research, as well as to better understand the evolutionary ecology and history of the species. A draft Murray cod genome of 633 Mbp (N₅₀ = 109 974bp; BUSCO and CEGMA completeness of 94.2% and 91.9%, respectively) with an estimated 148 Mbp of putative repetitive sequences was assembled from the combined sequencing data of 2 fish individuals with an identical maternal lineage; 47.2 Gb of Illumina HiSeq data and 804 Mb of Nanopore data were generated from the first individual while 23.2 Gb of Illumina MiSeq data were generated from the second individual. The inclusion of Nanopore reads for scaffolding followed by subsequent gap-closing using Illumina data led to a 29% reduction in the number of scaffolds and a 55% and 54% increase in the scaffold and contig N₅₀, respectively. We also report the first transcriptome of Murray cod that was subsequently used to annotate the Murray cod genome, leading to the identification of 26 539 protein-coding genes. We present the whole genome of the Murray cod and anticipate this will be a catalyst for a range of genetic, genomic, and phylogenetic studies of the Murray cod and more generally other fish species of the Percichthydae family.

Ecological disturbance influences adaptive divergence despite high gene flow in golden perch (Macquaria ambigua): Implications for management and resilience to climate change

Populations that are adaptively divergent but maintain high gene flow may have greater resilience to environmental change as gene flow allows the spread of alleles that have already been tested elsewhere. In addition, populations naturally subjected to ecological disturbance may already hold resilience to future environmental change. Confirming this necessitates ecological genomic studies of high dispersal, generalist species. Here we perform one such study on golden perch (Macquaria ambigua) in the Murray-Darling Basin (MDB), Australia, using a genome-wide SNP data set. The MDB spans across arid to wet and temperate to subtropical environments, with low to high ecological disturbance in the form of low to high hydrological variability. We found high gene flow across the basin and three populations with low neutral differentiation. Genotype-environment association analyses detected adaptive divergence predominantly linked to an arid region with highly variable riverine flow, and candidate loci included functions related to fat storage, stress and molecular or tissue repair. The high connectivity of golden perch in the MDB will likely allow locally adaptive traits in its most arid and hydrologically variable environment to spread and be selected in localities that are predicted to become arid and hydrologically variable in future climates. High connectivity in golden perch is likely due to their generalist life history and efforts of fisheries management. Our study adds to growing evidence of adaptation in the face of gene flow and highlights the importance of considering ecological disturbance and adaptive divergence in biodiversity management.

Severe consequences of habitat fragmentation on genetic diversity of an endangered Australian freshwater fish: A call for assisted gene flow

Genetic diversity underpins the ability of populations to persist and adapt to environmental changes. Substantial empirical data show that genetic diversity rapidly deteriorates in small and isolated populations due to genetic drift, leading to reduction in adaptive potential and fitness and increase in inbreeding. Assisted gene flow (e.g. via translocations) can reverse these trends, but lack of data on fitness loss and fear of impairing population "uniqueness" often prevents managers from acting. Here, we use population genetic and riverscape genetic analyses and simulations to explore the consequences of extensive habitat loss and fragmentation on population genetic diversity and future population trajectories of an endangered Australian freshwater fish, Macquarie perch Macquaria australasica. Using guidelines to assess the risk of outbreeding depression under admixture, we develop recommendations for population management, identify populations requiring genetic rescue and/or genetic restoration and potential donor sources. We found that most remaining populations of Macquarie perch have low genetic diversity, and effective population sizes below the threshold required to retain adaptive potential. Our simulations showed that under management inaction, smaller populations of Macquarie perch will face inbreeding depression within a few decades, but regular small-scale translocations will rapidly rescue populations from inbreeding depression and increase adaptive potential through genetic restoration. Despite the lack of data on fitness loss, based on our genetic data for Macquarie perch populations, simulations and empirical results from other systems, we recommend regular and frequent translocations among remnant populations within catchments. These translocations will emulate the effect of historical gene flow and improve population persistence through decrease in demographic and genetic stochasticity. Increasing population genetic connectivity within each catchment will help to maintain large effective population sizes and maximize species adaptive potential. The approach proposed here could be readily applicable to genetic management of other threatened species to improve their adaptive potential.

Landscape genetics informs mesohabitat preference and conservation priorities for a surrogate indicator species in a highly fragmented river system

Poor dispersal species represent conservative benchmarks for biodiversity management because they provide insights into ecological processes influenced by habitat fragmentation that are less evident in more dispersive organisms. Here we used the poorly dispersive and threatened river blackfish (Gadopsis marmoratus) as a surrogate indicator system for assessing the effects of fragmentation in highly modified river basins and for prioritizing basin-wide management strategies. We combined individual, population and landscape-based approaches to analyze genetic variation in samples spanning the distribution of the species in Australia's Murray-Darling Basin, one of the world's most degraded freshwater systems. Our results indicate that G. marmoratus displays the hallmark of severe habitat fragmentation with notably scattered, small and demographically isolated populations with very low genetic diversity-a pattern found not only between regions and catchments but also between streams within catchments. By using hierarchically nested population sampling and assessing relationships between genetic uniqueness and genetic diversity across populations, we developed a spatial management framework that includes the selection of populations in need of genetic rescue. Landscape genetics provided an environmental criterion to identify associations between landscape features and ecological processes. Our results further our understanding of the impact that habitat quality and quantity has on habitat specialists with similarly low dispersal. They should also have practical applications for prioritizing both large- and small-scale conservation management actions for organisms inhabiting highly fragmented ecosystems.

swinger: a user-friendly computer program to establish captive breeding groups that minimize relatedness without pedigree information

Captive breeding programmes are often a necessity for the continued persistence of a population or species. They typically have the goal of maintaining genetic diversity and minimizing inbreeding. However, most captive breeding programmes have been based on the assumption that the founding breeders are unrelated and outbred, even though in situ anthropogenic impacts often mean these founders may have high relatedness and substantial inbreeding. In addition, polygamous group-breeding species in captivity often have uncertain pedigrees, making it difficult to select the group composition for subsequent breeding. Molecular-based estimates of relatedness and inbreeding may instead be used to select breeding groups (≥two individuals) that minimize relatedness and filter out inbred individuals. swinger constructs breeding groups based on molecular estimates of relatedness and inbreeding. The number of possible combinations of breeding groups quickly becomes intractable by hand. swinger was designed to overcome this major issue in ex situ conservation biology. The user can specify parameters within swinger to reach breeding solutions that suit the mating system of the target species and available resources. We provide evidence of the efficiency of the software with an empirical example and using simulations. The only data required are a typical molecular marker data set, such as a microsatellite or SNP data set, from which estimates of inbreeding and pairwise relatedness may be obtained. Such molecular data sets are becoming easier to gather from non-model organisms with next-generation sequencing technology. swinger is an open-source software with a user-friendly interface and is available at http://www.molecularecology.flinders.edu.au/molecular-ecology-lab/software/swinger/swinger/ and https://github.com/Yuma248/Swinger.