Katherine A. Harrisson

Using genomics to characterize evolutionary potential for conservation of wild populations

Genomics promises exciting advances towards the important conservation goal of maximizing evolutionary potential, notwithstanding associated challenges. Here, we explore some of the complexity of adaptation genetics and discuss the strengths and limitations of genomics as a tool for characterizing evolutionary potential in the context of conservation management. Many traits are polygenic and can be strongly influenced by minor differences in regulatory networks and by epigenetic variation not visible in DNA sequence. Much of this critical complexity is difficult to detect using methods commonly used to identify adaptive variation, and this needs appropriate consideration when planning genomic screens, and when basing management decisions on genomic data. When the genomic basis of adaptation and future threats are well understood, it may be appropriate to focus management on particular adaptive traits. For more typical conservations scenarios, we argue that screening genome‐wide variation should be a sensible approach that may provide a generalized measure of evolutionary potential that accounts for the contributions of small‐effect loci and cryptic variation and is robust to uncertainty about future change and required adaptive response(s). The best conservation outcomes should be achieved when genomic estimates of evolutionary potential are used within an adaptive management framework.

Disrupted fine‐scale population processes in fragmented landscapes despite large‐scale genetic connectivity for a widespread and common cooperative breeder: the superb fairy‐wren (Malurus cyaneus)

Understanding how habitat fragmentation affects population processes (e.g. dispersal) at different spatial scales is of critical importance to conservation. We assessed the effects of habitat fragmentation on dispersal and regional and fine-scale population structure in a currently widespread and common cooperatively breeding bird species found across south-eastern Australia, the superb fairy-wren Malurus cyaneus. Despite its relative abundance and classification as an urban tolerant species, the superb fairy-wren has declined disproportionately from low tree-cover agricultural landscapes across the Box-Ironbark region of north-central Victoria, Australia. Loss of genetic connectivity and disruption to its complex social system may be associated with the decline of this species from apparently suitable habitat in landscapes with low levels of tree cover. To assess whether reduced structural connectivity has had negative consequences for genetic connectivity in the superb fairy-wren, we used a landscape-scale approach to compare patterns of genetic diversity and gene flow at large (landscape/regional) and fine (site-level) spatial scales. In addition, using genetic distances, for each sex, we tested landscape models of decreased dispersal through treeless areas (isolation-by-resistance) while controlling for the effect of isolation-by-distance. Landscape models indicated that larger-scale gene flow across the Box-Ironbark region was constrained by distance rather than by lack of structural connectivity. Nonetheless, a pattern of isolation-by-resistance for males (the less-dispersive sex) and lower genetic diversity and higher genetic similarity within sites in low-cover fragmented landscapes indicated disruption to fine-scale gene flow mechanisms and/or mating systems. Although loss of structural connectivity did not appear to impede gene flow at larger spatial scales, fragmentation appeared to affect fine-scale population processes (e.g. local gene flow mechanisms and/or mating systems) adversely and may contribute to the decline of superb fairy-wrens in fragmented landscapes in the Box-Ironbark region.

Fine-Scale Genetic Response to Landscape Change in a Gliding Mammal

Understanding how populations respond to habitat loss is central to conserving biodiversity. Population genetic approaches enable the identification of the symptoms of population disruption in advance of population collapse. However, the spatio-temporal scales at which population disruption occurs are still too poorly known to effectively conserve biodiversity in the face of human-induced landscape change. We employed microsatellite analysis to examine genetic structure and diversity over small spatial (mostly 1-50 km) and temporal scales (20-50 years) in the squirrel glider (Petaurus norfolcensis), a gliding mammal that is commonly subjected to a loss of habitat connectivity. We identified genetically differentiated local populations over distances as little as 3 km and within 30 years of landscape change. Genetically isolated local populations experienced the loss of genetic diversity, and significantly increased mean relatedness, which suggests increased inbreeding. Where tree cover remained, genetic differentiation was less evident. This pattern was repeated in two landscapes located 750 km apart. These results lend support to other recent studies that suggest the loss of habitat connectivity can produce fine-scale population genetic change in a range of taxa. This gives rise to the prediction that many other vertebrates will experience similar genetic changes. Our results suggest the future collapse of local populations of this gliding mammal is likely unless habitat connectivity is maintained or restored. Landscape management must occur on a fine-scale to avert the erosion of biodiversity.

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.

Pleistocene divergence across a mountain range and the influence of selection on mitogenome evolution in threatened Australian freshwater cod species

Climatic differences across a taxon’s range may be associated with specific bioenergetic demands and may result in genetics-based metabolic adaptation, particularly in aquatic ectothermic organisms that rely on heat exchange with the environment to regulate key physiological processes. Extending down the east coast of Australia, the Great Dividing Range (GDR) has a strong influence on climate and the evolutionary history of freshwater fish species. Despite the GDR acting as a strong contemporary barrier to fish movement, many species, and species with shared ancestries, are found on both sides of the GDR, indicative of historical dispersal events. We sequenced complete mitogenomes from the four extant species of the freshwater cod genus Maccullochella, two of which occur on the semi-arid, inland side of the GDR, and two on the mesic coastal side. We constructed a dated phylogeny and explored the relative influences of purifying and positive selection in the evolution of mitogenome divergence among species. Results supported mid- to late-Pleistocene divergence of Maccullochella across the GDR (220–710 thousand years ago), bringing forward previously reported dates. Against a background of pervasive purifying selection, we detected potentially functionally relevant fixed amino acid differences across the GDR. Although many amino acid differences between inland and coastal species may have become fixed under relaxed purifying selection in coastal environments rather than positive selection, there was evidence of episodic positive selection acting on specific codons in the Mary River coastal lineage, which has consistently experienced the warmest and least extreme climate in the genus.

Scope for genetic rescue of an endangered subspecies though re-establishing natural gene flow with another subspecies

Genetic diversity is positively linked to the viability and evolutionary potential of species but is often compromised in threatened taxa. Genetic rescue by gene flow from a more diverse or differentiated source population of the same species can be an effective strategy for alleviating inbreeding depression and boosting evolutionary potential. The helmeted honeyeater Lichenostomus melanops cassidix is a critically endangered subspecies of the common yellow‐tufted honeyeater. Cassidix has declined to a single wild population of ~130 birds, despite being subject to intensive population management over recent decades. We assessed changes in microsatellite diversity in cassidix over the last four decades and used population viability analysis to explore whether genetic rescue through hybridization with the neighbouring Lichenostomus melanops gippslandicus subspecies constitutes a viable conservation strategy. The contemporary cassidix population is characterized by low genetic diversity and effective population size (Ne < 50), suggesting it is vulnerable to inbreeding depression and will have limited capacity to evolve to changing environments. We find that gene flow from gippslandicus to cassidix has declined substantially relative to pre‐1990 levels and argue that natural levels of gene flow between the two subspecies should be restored. Allowing gene flow (~4 migrants per generation) from gippslandicus into cassidix (i.e. genetic rescue), in combination with continued annual release of captive‐bred cassidix (i.e. demographic rescue), should lead to positive demographic and genetic outcomes. Although we consider the risk of outbreeding depression to be low, we recommend that genetic rescue be managed within the context of the captive breeding programme, with monitoring of outcomes.

Signatures of polygenic adaptation associated with climate across the range of a threatened fish species with high genetic connectivity

Adaptive differences across species’ ranges can have important implications for population persistence and conservation management decisions. Despite advances in genomic technologies, detecting adaptive variation in natural populations remains challenging. Key challenges in gene–environment association studies involve distinguishing the effects of drift from those of selection and identifying subtle signatures of polygenic adaptation. We used paired‐end restriction site‐associated DNA sequencing data (6,605 biallelic single nucleotide polymorphisms; SNPs) to examine population structure and test for signatures of adaptation across the geographic range of an iconic Australian endemic freshwater fish species, the Murray cod Maccullochella peelii. Two univariate gene–association methods identified 61 genomic regions associated with climate variation. We also tested for subtle signatures of polygenic adaptation using a multivariate method (redundancy analysis; RDA). The RDA analysis suggested that climate (temperature‐ and precipitation‐related variables) and geography had similar magnitudes of effect in shaping the distribution of SNP genotypes across the sampled range of Murray cod. Although there was poor agreement among the candidate SNPs identified by the univariate methods, the top 5% of SNPs contributing to significant RDA axes included 67% of the SNPs identified by univariate methods. We discuss the potential implications of our findings for the management of Murray cod and other species generally, particularly in relation to informing conservation actions such as translocations to improve evolutionary resilience of natural populations. Our results highlight the value of using a combination of different approaches, including polygenic methods, when testing for signatures of adaptation in landscape genomic studies.

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

Abstract 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 (N50 = 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 N50, 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.

Identifying environmental correlates of intraspecific genetic variation.

Genetic variation is critical to the persistence of populations and their capacity to adapt to environmental change. The distribution of genetic variation across a species’ range can reveal critical information that is not necessarily represented in species occurrence or abundance patterns. We identified environmental factors associated with the amount of intraspecific, individual-based genetic variation across the range of a widespread freshwater fish species, the Murray cod Maccullochella peelii. We used two different approaches to statistically quantify the relative importance of predictor variables, allowing for nonlinear relationships: a random forest model and a Bayesian approach. The latter also accounted for population history. Both approaches identified associations between homozygosity by locus and both disturbance to the natural flow regime and mean annual flow. Homozygosity by locus was negatively associated with disturbance to the natural flow regime, suggesting that river reaches with more disturbed flow regimes may support larger, more genetically diverse populations. Our findings are consistent with the hypothesis that artificially induced perennial flows in regulated channels may provide greater and more consistent habitat and reduce the frequency of population bottlenecks that can occur frequently under the highly variable and unpredictable natural flow regime of the system. Although extensive river regulation across eastern Australia has not had an overall positive effect on Murray cod numbers over the past century, regulation may not represent the primary threat to Murray cod survival. Instead, pressures other than flow regulation may be more critical to the persistence of Murray cod (for example, reduced frequency of large floods, overfishing and chemical pollution).

Genetic structure and sex-biased dispersal of a declining cooperative-breeder, the Grey-crowned Babbler, Pomatostomus temporalis, at the southern edge of its range

Abstract Loss and fragmentation of habitat can disrupt genetic exchange between populations, which is reflected in changes to the genetic structure of populations. The Grey-crowned Babbler (Pomatostomus temporalis) is a cooperatively breeding woodland bird, once common and widespread in south-eastern Australia. The species has suffered population declines of >90% across its southern distribution as a result of loss and fragmentation of habitat. We investigated patterns of genetic diversity and population structure of Grey-crowned Babblers in fragmented habitats at the southernmost extent of its range. We sampled blood from 135 individual Babblers from 39 groups stratified into six subpopulations in three regions. Genotypic data were used to estimate genetic diversity, population substructure, local relatedness and dispersal patterns. Individuals showed high heterozygosity within regions, and varying numbers of private alleles among regions suggested differences in levels of connectivity between regions. Four genetic clusters revealed population substructure consistent with treeless landscapes acting as strong barriers to gene flow. In contrast to previous studies, we identified a male-biased dispersal pattern and significant isolation-by-distance patterns for females at fine spatial scales. We recommend that conservation plans for this species incorporate opportunities to increase and enhance corridor areas to facilitate genetic exchange among subpopulations.