Swinging back from the brink? Polygamous mating strategies revealed for an iconic threatened freshwater fish.

Catastrophic fish death events are increasing in frequency and severity globally. A series of major recent fish deaths in the semi-arid lower Darling-Baaka river system (LDBR) of Australia are emblematic of these issues with tens of millions of native fish perishing. In 2018-2019 there was a major death event for Australia's largest freshwater fish, Murray cod (Maccullochella peelii). To aid the recovery and guide restoration activities of local Murray cod populations, it is essential to gather information on the mating strategies and effective population size following the fish death event. After the fish deaths, we collected larvae during the 2020 and 2021 breeding seasons and used single nucleotide polymorphisms (SNPs) to provide insight mating strategies and to estimate effective population size. Larvae were detected in both years along the entire length of the LDBR. Sixteen percent of the inferred breeding individuals were found to contribute to multiple pairings, confirming a complex and polygamous mating system. A high frequency of polygamy was evident both within and between years with 100 % polygamy identified among parents that produced offspring in both 2020 and 2021 and 95 % polygamy identified among parents involved in multiple spawning events within years. Post-larval Murray cod samples collected between 2016 and 2021 were co-analysed to further inform kinship patterns. Again, monogamy was rare with no confirmed cases of the same male-female pair contributing to multiple breeding events within or between seasons. Effective population size based on Murray cod collected after the fish death event was estimated at 721.6 (CI 471-1486), though this has likely declined following a subsequent catastrophic fish death event in the LDBR in March 2023. Our data provide insight into the variability of Murray cod mating strategies, and we anticipate that this knowledge will assist in planning conservation actions to ultimately help recover a species in crisis.

Labile sex chromosomes in the Australian freshwater fish family Percichthyidae.

Sex-specific ecology has management implications, but rapid sex-chromosome turnover in fishes hinders sex-marker development for monomorphic species. We used annotated genomes and reduced-representation sequencing data for two Australian percichthyids, Macquarie perch Macquaria australasica and golden perch M. ambigua, and whole genome resequencing for 50 Macquarie perch of each sex, to identify sex-linked loci and develop an affordable sexing assay. In silico pool-seq tests of 1,492,004 Macquarie perch SNPs revealed that a 275-kb scaffold was enriched for gametologous loci. Within this scaffold, 22 loci were sex-linked in a predominantly XY system, with females being homozygous for the X-linked allele at all 22, and males having the Y-linked allele at >7. Seven XY-gametologous loci (all males, but no females, are heterozygous or homozygous for the male-specific allele) were within a 146-bp region. A PCR-RFLP sexing assay targeting one Y-linked SNP, tested in 66 known-sex Macquarie perch and two of each sex of three confamilial species, plus amplicon sequencing of 400 bp encompassing the 146-bp region, revealed that the few sex-linked positions differ between species and between Macquarie perch populations. This indicates sex-chromosome lability in Percichthyidae, supported by nonhomologous scaffolds containing sex-linked loci for Macquarie- and golden perches. The present resources facilitate genomic research in Percichthyidae, including formulation of hypotheses about candidate genes of interest such as transcription factor SOX1b that occurs in the 275-kb scaffold ~38 kb downstream of the 146-bp region containing seven XY-gametologous loci. Sex-linked markers will be useful for determining genetic sex in some populations and studying sex chromosome turnover.

A novel framework for evaluating in situ breeding management strategies in endangered populations.

Conservation breeding management aims to reduce inbreeding and maximize the retention of genetic diversity in endangered populations. However, breeding management of wild populations is still rare, and there is a need for approaches that provide data-driven evidence of the likelihood of success of alternative in situ strategies. Here, we provide an analytical framework that uses in silico simulations to evaluate, for real wild populations, (i) the degree of population-level inbreeding avoidance, (ii) the genetic quality of mating pairs, and (iii) the potential genetic benefits of implementing two breeding management strategies. The proposed strategies aim to improve the genetic quality of breeding pairs by splitting detrimental pairs and allowing the members to re-pair in different ways. We apply the framework to the wild population of the Critically Endangered helmeted honeyeater by combining genomic data and field observations to estimate the inbreeding (i.e., pair-kinship) and genetic quality (i.e., Mate Suitability Index) of all mating pairs for seven consecutive breeding seasons. We found no evidence of population-level inbreeding avoidance and that ~91.6% of breeding pairs were detrimental to the genetic health of the population. Furthermore, the framework revealed that neither proposed management strategy would significantly improve the genetic quality or reduce inbreeding of the mating pairs in this population. Our results demonstrate the usefulness of our analytical framework for testing the efficacy of different in situ breeding management strategies and for making evidence-based management decisions.

Lifetime Fitness Costs of Inbreeding and Being Inbred in a Critically Endangered Bird.

Reduced fitness as a result of inbreeding is a major threat facing many species of conservation concern [1-4]. However, few case studies for assessing the magnitude of inbreeding depression in the wild means that its relative importance as a risk factor for population persistence remains under-appreciated [5]. The increasing availability and affordability of genomic technologies provide new opportunities to address knowledge gaps around the magnitude and manifestation of inbreeding depression in wild populations [6-12]. Here, we combine over three decades of individual lifetime reproductive data and genomic data to estimate the relative lifetime and short-term fitness costs of both being inbred and engaging in inbreeding in the last wild population (<250 individuals remaining) of an iconic and critically endangered bird: the helmeted honeyeater Lichenostomus melanops cassidix. The magnitude of inbreeding depression was substantial: the mean predicted lifetime reproductive success of the most inbred (homozygosity = 0.82) individuals was on average 87%-90% lower than that of the least inbred (homozygosity = 0.75). For individual reproductive events and lifetime measures, we provide rare empirical evidence that pairing with a genetically dissimilar individual can reduce fitness costs associated with being an inbred individual. By comparing lifetime and short-term fitness measures, we demonstrate how short-term measures of reproductive success that are associated with only weak signatures of inbreeding depression can still underlie stronger lifetime effects. Our study represents a valuable case study, highlighting the critical importance of inbreeding depression as a factor influencing the immediate viability of populations in threatened species management.

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.

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.

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.

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 (N(e) < 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.

Species- and sex-specific connectivity effects of habitat fragmentation in a suite of woodland birds.

Loss of functional connectivity following habitat loss and fragmentation could drive species declines. A comprehensive understanding of fragmentation effects on functional connectivity of an ecological assemblage requires investigation of multiple species with different mobilities, at different spatial scales, for each sex, and in different landscapes. Based on published data on mobility and ecological responses to fragmentation of 10 woodland-dependent birds, and using simulation studies, we predicted that (1) fragmentation would impede dispersal and gene flow of eight "decliners" (species that disappear from suitable patches when landscape-level tree cover falls below species-specific thresholds), but not of two "tolerant" species (whose occurrence in suitable habitat patches is independent of landscape tree cover); and that fragmentation effects would be stronger (2) in the least mobile species, (3) in the more philopatric sex, and (4) in the more fragmented region. We tested these predictions by evaluating spatially explicit isolation-by-landscape-resistance models of gene flow in fragmented landscapes across a 50 x 170 km study area in central Victoria, Australia, using individual and population genetic distances. To account for sex-biased dispersal and potential scale- and configuration-specific effects, we fitted models specific to sex and geographic zones. As predicted, four of the least mobile decliners showed evidence of reduced genetic connectivity. The responses were strongly sex specific, but in opposite directions in the two most sedentary species. Both tolerant species and (unexpectedly) four of the more mobile decliners showed no reduction in gene flow. This is unlikely to be due to time lags because more mobile species develop genetic signatures of fragmentation faster than do less mobile ones. Weaker genetic effects were observed in the geographic zone with more aggregated vegetation, consistent with gene flow being unimpeded by landscape structure. Our results indicate that for all but the most sedentary species in our system, the movement of the more dispersive sex (females in most cases) maintains overall genetic connectivity across fragmented landscapes in the study area, despite some small-scale effects on the more philopatric sex for some species. Nevertheless, to improve population viability for the less mobile bird species, structural landscape connectivity must be increased.

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.

Does reduced mobility through fragmented landscapes explain patch extinction patterns for three honeyeaters?

Habitat loss and associated fragmentation are major drivers of biodiversity decline, and understanding how they affect population processes (e.g. dispersal) is an important conservation goal. In a large-scale test employing 10 × 10 km units of replication, three species of Australian birds, the fuscous honeyeater, yellow-tufted honeyeater and white-plumed honeyeater, responded differently to fragmentation. The fuscous and yellow-tufted honeyeaters are 'decliners' that disappeared from suitable habitat in landscapes where levels of tree-cover fell below critical thresholds of 17 and 8%, respectively. The white-plumed honeyeater is a 'tolerant' species whose likelihood of occurrence in suitable habitat was independent of landscape-level tree-cover. To determine whether the absence of the two decliner species in low tree-cover landscapes can be explained by reduced genetic connectivity, we looked for signatures of reduced mobility and gene flow in response to fragmentation across agricultural landscapes in the Box-Ironbark region of north-central Victoria, Australia. We compared patterns of genetic diversity and population structure at the regional scale and across twelve 100 km(2) landscapes with different tree-cover extents. We used genetic data to test landscape models predicting reduced dispersal through the agricultural matrix. We tested for evidence of sex-biased dispersal and sex-specific responses to fragmentation. Reduced connectivity may have contributed to the disappearance of the yellow-tufted honey-eater from low tree-cover landscapes, as evidenced by male bias and increased relatedness among males in low tree-cover landscapes and signals of reduced gene flow and mobility through the agricultural matrix. We found no evidence for negative effects of fragmentation on gene flow in the other decliner, the fuscous honeyeater, suggesting that undetected pressures act on this species. As expected, there was no evidence for decreased movement through fragmented landscapes for the tolerant white-plumed honeyeater. We demonstrated effects of habitat loss and fragmentation (stronger patterns of genetic differentiation, increased relatedness among males) on the yellow-tufted honeyeater above the threshold at which probability of occurrence dropped. Increasing extent and structural connectivity of habitat should be an appropriate management action for this species and other relatively sedentary woodland specialist species for which it can be taken as representative.

Limited population structure, genetic drift and bottlenecks characterise an endangered bird species in a dynamic, fire-prone ecosystem.

Fire is a major disturbance process in many ecosystems world-wide, resulting in spatially and temporally dynamic landscapes. For populations occupying such environments, fire-induced landscape change is likely to influence population processes, and genetic patterns and structure among populations. The Mallee Emu-wren Stipiturus mallee is an endangered passerine whose global distribution is confined to fire-prone, semi-arid mallee shrublands in south-eastern Australia. This species, with poor capacity for dispersal, has undergone a precipitous reduction in distribution and numbers in recent decades. We used genetic analyses of 11 length-variable, nuclear loci to examine population structure and processes within this species, across its global range. Populations of the Mallee Emu-wren exhibited a low to moderate level of genetic diversity, and evidence of bottlenecks and genetic drift. Bayesian clustering methods revealed weak genetic population structure across the species' range. The direct effects of large fires, together with associated changes in the spatial and temporal patterns of suitable habitat, have the potential to cause population bottlenecks, serial local extinctions and subsequent recolonisation, all of which may interact to erode and homogenise genetic diversity in this species. Movement among temporally and spatially shifting habitat, appears to maintain long-term genetic connectivity. A plausible explanation for the observed genetic patterns is that, following extensive fires, recolonisation exceeds in-situ survival as the primary driver of population recovery in this species. These findings suggest that dynamic, fire-dominated landscapes can drive genetic homogenisation of populations of species with low-mobility and specialised habitat that otherwise would be expected to show strongly structured populations. Such effects must be considered when formulating management actions to conserve species in fire-prone systems.

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.