Chris J. Brauer

Riverscape genomics of a threatened fish across a hydroclimatically heterogeneous river basin.

Understanding how natural selection generates and maintains adaptive genetic diversity in heterogeneous environments is key to predicting the evolutionary response of populations to rapid environmental change. Detecting selection in complex spatial environments remains challenging, especially for threatened species where the effects of strong genetic drift may overwhelm signatures of selection. We carried out a basinwide riverscape genomic analysis in the threatened southern pygmy perch (Nannoperca australis), an ecological specialist with low dispersal potential. High‐resolution environmental data and 5162 high‐quality filtered SNPs were used to clarify spatial population structure and to assess footprints of selection associated with a steep hydroclimatic gradient and with human disturbance across the naturally and anthropogenically fragmented Murray–Darling Basin (Australia). Our approach included FST outlier tests to define neutral loci, and a combination of spatially explicit genotype–environment association analyses to identify candidate adaptive loci while controlling for the effects of landscape structure and shared population history. We found low levels of genetic diversity and strong neutral population structure consistent with expectations based on spatial stream hierarchy and life history. In contrast, variables related to precipitation and temperature appeared as the most important environmental surrogates for putatively adaptive genetic variation at both regional and local scales. Human disturbance also influenced the variation in candidate loci for adaptation, but only at a local scale. Our study contributes to understanding of adaptive evolution along naturally and anthropogenically fragmented ecosystems. It also offers a tangible example of the potential contributions of landscape genomics for informing in situ and ex situ conservation management of biodiversity.

Catchment-Scale Conservation Units Identified for the Threatened Yarra Pygmy Perch (Nannoperca obscura) in Highly Modified River Systems

Habitat fragmentation caused by human activities alters metapopulation dynamics and decreases biological connectivity through reduced migration and gene flow, leading to lowered levels of population genetic diversity and to local extinctions. The threatened Yarra pygmy perch, Nannoperca obscura, is a poor disperser found in small, isolated populations in wetlands and streams of southeastern Australia. Modifications to natural flow regimes in anthropogenically-impacted river systems have recently reduced the amount of habitat for this species and likely further limited its opportunity to disperse. We employed highly resolving microsatellite DNA markers to assess genetic variation, population structure and the spatial scale that dispersal takes place across the distribution of this freshwater fish and used this information to identify conservation units for management. The levels of genetic variation found for N. obscura are amongst the lowest reported for a fish species (mean heterozygosity of 0.318 and mean allelic richness of 1.92). We identified very strong population genetic structure, nil to little evidence of recent migration among demes and a minimum of 11 units for conservation management, hierarchically nested within four major genetic lineages. A combination of spatial analytical methods revealed hierarchical genetic structure corresponding with catchment boundaries and also demonstrated significant isolation by riverine distance. Our findings have implications for the national recovery plan of this species by demonstrating that N. obscura populations should be managed at a catchment level and highlighting the need to restore habitat and avoid further alteration of the natural hydrology.

A novel holistic framework for genetic-based captive-breeding and reintroduction programs

Research in reintroduction biology has provided a greater understanding of the often limited success of species reintroductions and highlighted the need for scientifically rigorous approaches in reintroduction programs. We examined the recent genetic-based captive-breeding and reintroduction literature to showcase the underuse of the genetic data gathered. We devised a framework that takes full advantage of the genetic data through assessment of the genetic makeup of populations before (past component of the framework), during (present component), and after (future component) captive-breeding and reintroduction events to understand their conservation potential and maximize their success. We empirically applied our framework to two small fishes: Yarra pygmy perch (Nannoperca obscura) and southern pygmy perch (Nannoperca australis). Each of these species has a locally adapted and geographically isolated lineage that is endemic to the highly threatened lower Murray-Darling Basin in Australia. These two populations were rescued during Australias recent decade-long Millennium Drought, when their persistence became entirely dependent on captive-breeding and subsequent reintroduction efforts. Using historical demographic analyses, we found differences and similarities between the species in the genetic impacts of past natural and anthropogenic events that occurred in situ, such as European settlement (past component). Subsequently, successful maintenance of genetic diversity in captivity-despite skewed brooder contribution to offspring-was achieved through carefully managed genetic-based breeding (present component). Finally, genetic monitoring revealed the survival and recruitment of released captive-bred offspring in the wild (future component). Our holistic framework often requires no additional data collection to that typically gathered in genetic-based breeding programs, is applicable to a wide range of species, advances the genetic considerations of reintroduction programs, and is expected to improve with the use of next-generation sequencing technology.

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

Distinguishing the relative influence of historic (i.e. natural) versus anthropogenic factors in metapopulation structure is an important but often overlooked step in management programs of threatened species. Biotas in freshwater wetlands and floodplains, such as those in the Murray–Darling Basin (MDB)—one of Australia’s most impacted ecosystems, are particularly susceptible to anthropogenic fragmentation. Here we present a comprehensive multilocus assessment of genetic variation in the threatened southern pygmy perch Nannoperca australis (578 individuals; 45 localities; microsatellite, allozyme and mitochondrial DNA datasets), an ecological specialist with low dispersal potential. We assess patterns of spatial structure and genetic diversity in populations spanning the highly fragmented MDB and test whether recent anthropogenic modification has disrupted range-wide connectivity. We detected strong and hierarchical population structure, very low genetic diversity and lack of contemporary gene flow across the MDB. In contrast, the apparent absence of pronounced or long-term phylogeographic structure suggests that observed population divergences generally do not reflect deeply historic natural fragmentation. Coalescent-based analyses supported this inference, revealing that divergence times between populations from the upper and lower MDB fall into the period of European settlement. It appears that the observed contemporary isolation of populations is partly explained by the severe modification of the MDB post-dating the onset of European settlement. Our integrated approach substantially improves the interpretation of how fragmentation impacts present-day biodiversity. It also provides novel contributions for risk-assessing management actions in the context of captive breeding and translocations of small freshwater fishes, a group of increasing global conservation concern.

Multi-generational evaluation of genetic diversity and parentage in captive southern pygmy perch ( Nannoperca australis )

Maintaining genetic diversity within captive breeding populations is a key challenge for conservation managers. We applied a multi-generational genetic approach to the captive breeding program of an endangered Australian freshwater fish, the southern pygmy perch (Nannoperca australis). During previous work, fish from the lower Murray-Darling Basin were rescued before drought exacerbated by irrigation resulted in local extinction. This endemic lineage of the species was captive-bred in genetically designed groups, and equal numbers of F1 individuals were reintroduced to the wild with the return of favourable habitat. Here, we implemented a contingency plan by continuing the genetic-based captive breeding in the event that a self-sustaining wild population was not established. F1 individuals were available as putative breeders from the subset of groups that produced an excess of fish in the original restoration program. We used microsatellite-based parentage analyses of these F1 fish to form breeding groups that minimized inbreeding. We assessed their subsequent parental contribution to F2 individuals and the maintenance of genetic diversity. We found skewed parental contribution to F2 individuals, yet minimal loss of genetic diversity from their parents. However, the diversity was substantially less than that of the original rescued population. We attribute this to the unavoidable use of F1 individuals from a limited number of the original breeding groups. Alternative genetic sources for supplementation or reintroduction should be assessed to determine their suitability. The genetic fate of the captive-bred population highlights the strong need to integrate DNA-based tools for monitoring and adaptive management of captive breeding programs.

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.

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.

Comparative ecological transcriptomics and the contribution of gene expression to the evolutionary potential of a threatened fish

Understanding whether small populations with low genetic diversity can respond to rapid environmental change via phenotypic plasticity is an outstanding research question in biology. RNA sequencing (RNA‐seq) has recently provided the opportunity to examine variation in gene expression, a surrogate for phenotypic variation, in nonmodel species. We used a comparative RNA‐seq approach to assess expression variation within and among adaptively divergent populations of a threatened freshwater fish, Nannoperca australis, found across a steep hydroclimatic gradient in the Murray–Darling Basin, Australia. These populations evolved under contrasting selective environments (e.g., dry/hot lowland; wet/cold upland) and represent opposite ends of the species’ spectrum of genetic diversity and population size. We tested the hypothesis that environmental variation among isolated populations has driven the evolution of divergent expression at ecologically important genes using differential expression (DE) analysis and an anova‐based comparative phylogenetic expression variance and evolution model framework based on 27,425 de novo assembled transcripts. Additionally, we tested whether gene expression variance within populations was correlated with levels of standing genetic diversity. We identified 290 DE candidate transcripts, 33 transcripts with evidence for high expression plasticity, and 50 candidates for divergent selection on gene expression after accounting for phylogenetic structure. Variance in gene expression appeared unrelated to levels of genetic diversity. Functional annotation of the candidate transcripts revealed that variation in water quality is an important factor influencing expression variation for N. australis. Our findings suggest that gene expression variation can contribute to the evolutionary potential of small populations.

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.

On the roles of landscape heterogeneity and environmental variation in determining population genomic structure in a dendritic system

Dispersal and natural selection are key evolutionary processes shaping the distribution of phenotypic and genetic diversity. For species inhabiting complex spatial environments however, it is unclear how the balance between gene flow and selection may be influenced by landscape heterogeneity and environmental variation. Here, we evaluated the effects of dendritic landscape structure and the selective forces of hydroclimatic variation on population genomic parameters for the Murray River rainbowfish, Melanotaenia fluviatilis across the Murray–Darling Basin, Australia. We genotyped 249 rainbowfish at 17,503 high‐quality SNP loci and integrated these with models of network connectivity and high‐resolution environmental data within a riverscape genomics framework. We tested competing models of gene flow before using multivariate genotype–environment association (GEA) analysis to test for signals of adaptive divergence associated with hydroclimatic variation. Patterns of neutral genetic variation were consistent with expectations based on the stream hierarchy model and M. fluviatilis’ moderate dispersal ability. Models incorporating dendritic network structure suggested that landscape heterogeneity is a more important factor determining connectivity and gene flow than waterway distance. Extending these results, we also introduce a novel approach to controlling for the unique effects of dendritic network structure in GEA analyses of populations of aquatic species. We identified 146 candidate loci potentially underlying a polygenic adaptive response to seasonal fluctuations in stream flow and variation in the relative timing of temperature and precipitation extremes. Our findings underscore an emerging predominant role for seasonal variation in hydroclimatic conditions driving local adaptation and are relevant for informing proactive conservation management.