Planning and implementing genetic rescue of an endangered freshwater fish population in a regulated river, where low flow reduces breeding opportunities and may trigger inbreeding depression.

Augmenting depleted genetic diversity can improve the fitness and evolutionary potential of wildlife populations, but developing effective management approaches requires genetically monitored test cases. One such case is the small, isolated and inbred Cotter River population of an endangered Australian freshwater fish, the Macquarie perch Macquaria australasica, which over 3 years (2017-2019) received 71 translocated migrants from a closely related, genetically more diverse population. We used genetic monitoring to test whether immigrants bred, interbred with local fish and augmented population genetic diversity. We also investigated whether levels of river flow affected recruitment, inbreeding depression and juvenile dispersal. Fish length was used to estimate the age, birth year cohort and growth of 524 individuals born between 2016 and 2020 under variable flow conditions. DArT genome-wide genotypes were used to assess individual ancestry, heterozygosity, short-term effective population size and identify parent-offspring and full-sibling families. Of 442 individuals born after translocations commenced, only two (0.45%) were of mixed ancestry; these were half-sibs with one translocated parent in common. Numbers of breeders and genetic diversity for five birth year cohorts of the Cotter River fish were low, especially in low-flow years. Additionally, individuals born in the year of lowest flow evidently suffered from inbreeding depression for juvenile growth. The year of highest flow was associated with the largest number of breeders, lowest inbreeding in the offspring and greatest juvenile dispersal distances. Genetic diversity decreased in the upstream direction, flagging restricted access of breeders to the most upstream breeding sites, exacerbated by low river flow. Our results suggest that the effectiveness of translocations could be increased by focussing on upstream sites and moving more individuals per year; using riverine sources should be considered. Our results indicate that river flow sufficient to facilitate fish movement through the system would increase the number of breeders, promote individuals' growth, reduce inbreeding depression and promote genetic rescue.

Easy-to-use R functions to separate reduced-representation genomic datasets into sex-linked and autosomal loci, and conduct sex assignment.

Identifying sex-linked markers in genomic datasets is important because their presence in supposedly neutral autosomal datasets can result in incorrect estimates of genetic diversity, population structure and parentage. However, detecting sex-linked loci can be challenging, and available scripts neglect some categories of sex-linked variation. Here, we present new R functions to (1) identify and separate sex-linked loci in ZW and XY sex determination systems and (2) infer the genetic sex of individuals based on these loci. We tested these functions on genomic data for two bird and one mammal species and compared the biological inferences made before and after removing sex-linked loci using our function. We found that our function identified autosomal loci with ≥98.8% accuracy and sex-linked loci with an average accuracy of 87.8%. We showed that standard filters, such as low read depth and call rate, failed to remove up to 54.7% of sex-linked loci. This led to (i) overestimation of population FIS by up to 24%, and the number of private alleles by up to 8%; (ii) wrongly inferring significant sex differences in heterozygosity; (iii) obscuring genetic population structure and (iv) inferring ~11% fewer correct parentages. We discuss how failure to remove sex-linked markers can lead to incorrect biological inferences (e.g. sex-biased dispersal and cryptic population structure) and misleading management recommendations. For reduced-representation datasets with at least 15 known-sex individuals of each sex, our functions offer convenient resources to remove sex-linked loci and to sex the remaining individuals (freely available at https://github.com/drobledoruiz/conservation_genomics).

Chromosome-length genome assembly and linkage map of a critically endangered Australian bird: the helmeted honeyeater.

BACKGROUND: The helmeted honeyeater (Lichenostomus melanops cassidix) is a Critically Endangered bird endemic to Victoria, Australia. To aid its conservation, the population is the subject of genetic rescue. To understand, monitor, and modulate the effects of genetic rescue on the helmeted honeyeater genome, a chromosome-length genome and a high-density linkage map are required. RESULTS: We used a combination of Illumina, Oxford Nanopore, and Hi-C sequencing technologies to assemble a chromosome-length genome of the helmeted honeyeater, comprising 906 scaffolds, with length of 1.1 Gb and scaffold N50 of 63.8 Mb. Annotation comprised 57,181 gene models. Using a pedigree of 257 birds and 53,111 single-nucleotide polymorphisms, we obtained high-density linkage and recombination maps for 25 autosomes and Z chromosome. The total sex-averaged linkage map was 1,347 cM long, with the male map being 6.7% longer than the female map. Recombination maps revealed sexually dimorphic recombination rates (overall higher in males), with average recombination rate of 1.8 cM/Mb. Comparative analyses revealed high synteny of the helmeted honeyeater genome with that of 3 passerine species (e.g., 32 Hi-C scaffolds mapped to 30 zebra finch autosomes and Z chromosome). The genome assembly and linkage map suggest that the helmeted honeyeater exhibits a fission of chromosome 1A into 2 chromosomes relative to zebra finch. PSMC analysis showed a ∼15-fold decline in effective population size to ∼60,000 from mid- to late Pleistocene. CONCLUSIONS: The annotated chromosome-length genome and high-density linkage map provide rich resources for evolutionary studies and will be fundamental in guiding conservation efforts for the helmeted honeyeater.

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.

Using multiple sources during reintroduction of a locally extinct population benefits survival and reproduction of an endangered freshwater fish.

Through using different sources, population reintroductions can create genetically diverse populations at low risk of harmful inbreeding and well equipped for adaptation to future environments. Genetic variation from one source can mask locally nonoptimal alleles from another, thereby enhancing adaptive potential and population persistence. We assessed the outcomes in survival, growth and reproduction of using two differentiated sources (genetically diverse Yarra and moderately diverse Dartmouth) for translocations and stocking to reintroduce the endangered Australian freshwater Macquarie perch Macquaria australasica into the Ovens River. For stocking, same- and different-population parents ("cross-types") were used during hatchery production. Genetic samples and data on individual fish were collected over three years of monitoring the Ovens. We genetically assigned Ovens fish to their broodstock parents and tested whether cross-type and genetic dissimilarity between parents are associated with offspring survival, and whether cross-type and parental dissimilarity or individual genetic diversity are associated with somatic growth rates of stocked fish. We genetically identified translocated fish and assessed local recruit ancestry. Of 296 Ovens fish, 31.1% were inferred to be stocked, 1.3% translocated and 67.6% locally born. Cross-type strongly predicted survival of stocked offspring: those with two Yarra parents had the highest survival, followed by offspring with two-population, then Dartmouth, ancestry. Of the Ovens recruits, 59.5% had Yarra, 33.5% two-population and 7.0% Dartmouth ancestry, despite 67% of stocked and 98% of translocated fish originating from Dartmouth. Offspring with two Yarra parents grew faster than offspring of Dartmouth or two-population ancestry. Although Dartmouth fish appear to be less fit in the Ovens compared to Yarra fish, possibly due to deleterious variation or genetic or plastic maladaptation, they contribute to the reintroduced population through local interbreeding with Yarra fish and relatively high survival of stocked offspring of two-population ancestry. Thus, combining compatible stocks is likely to benefit restoration of other wildlife populations.

Evolutionary history and genetic connectivity across highly fragmented populations of an endangered daisy.

Conservation management can be aided by knowledge of genetic diversity and evolutionary history, so that ecological and evolutionary processes can be preserved. The Button Wrinklewort daisy (Rutidosis leptorrhynchoides) was a common component of grassy ecosystems in south-eastern Australia. It is now endangered due to extensive habitat loss and the impacts of livestock grazing, and is currently restricted to a few small populations in two regions >500 km apart, one in Victoria, the other in the Australian Capital Territory and nearby New South Wales (ACT/NSW). Using a genome-wide SNP dataset, we assessed patterns of genetic structure and genetic differentiation of 12 natural diploid populations. We estimated intrapopulation genetic diversity to scope sources for genetic management. Bayesian clustering and principal coordinate analyses showed strong population genetic differentiation between the two regions, and substantial substructure within ACT/NSW. A coalescent tree-building approach implemented in SNAPP indicated evolutionary divergence between the two distant regions. Among the populations screened, the last two known remaining Victorian populations had the highest genetic diversity, despite having among the lowest recent census sizes. A maximum likelihood population tree method implemented in TreeMix suggested little or no recent gene flow except potentially between very close neighbours. Populations that were more genetically distinctive had lower genetic diversity, suggesting that drift in isolation is likely driving population differentiation though loss of diversity, hence re-establishing gene flow among them is desirable. These results provide background knowledge for evidence-based conservation and support genetic rescue within and between regions to elevate genetic diversity and alleviate inbreeding.

Integrating Mitochondrial Aerobic Metabolism into Ecology and Evolution.

Biologists have long appreciated the critical role that energy turnover plays in understanding variation in performance and fitness among individuals. Whole-organism metabolic studies have provided key insights into fundamental ecological and evolutionary processes. However, constraints operating at subcellular levels, such as those operating within the mitochondria, can also play important roles in optimizing metabolism over different energetic demands and time scales. Herein, we explore how mitochondrial aerobic metabolism influences different aspects of organismal performance, such as through changing adenosine triphosphate (ATP) and reactive oxygen species (ROS) production. We consider how such insights have advanced our understanding of the mechanisms underpinning key ecological and evolutionary processes, from variation in life-history traits to adaptation to changing thermal conditions, and we highlight key areas for future research.

Accounting for cryptic population substructure enhances detection of inbreeding depression with genomic inbreeding coefficients: an example from a critically endangered marsupial.

Characterizing inbreeding depression in wildlife populations can be critical to their conservation. Coefficients of individual inbreeding can be estimated from genome-wide marker data. The degree to which sensitivity of inbreeding coefficients to population genetic substructure alters estimates of inbreeding depression in wild populations is not well understood. Using generalized linear models, we tested the power of two frequently used inbreeding coefficients that are calculated from genome-wide SNP markers, FH and F^III , to predict four fitness traits estimated over two decades in an isolated population of the critically endangered Leadbeater's possum. FH estimates inbreeding as excess observed homozygotes relative to equilibrium expectations, whereas F^III quantifies allelic similarity between the gametes that formed an individual, and upweights rare homozygotes. We estimated FH and F^III from 1,575 genome-wide SNP loci in individuals with fitness trait data (N = 179-237 per trait), and computed revised coefficients, FHby group and F^IIIby group , adjusted for population genetic substructure by calculating them separately within two different genetic groups of individuals identified in the population. Using FH or F^III in the models, inbreeding depression was detected for survival to sexual maturity, longevity and whether individuals bred during their lifetime. F^IIIby group (but not FHby group ) additionally revealed significant inbreeding depression for lifetime reproductive output (total offspring assigned to each individual). Estimates of numbers of lethal equivalents indicated substantial inbreeding load, but differing between inbreeding estimators. Inbreeding depression, declining population size, and low and declining genetic diversity suggest that genetic rescue may assist in preventing extinction of this unique Leadbeater's possum population.

Genomic evidence of neo-sex chromosomes in the eastern yellow robin.

BACKGROUND: Understanding sex-biased natural selection can be enhanced by access to well-annotated chromosomes including ones inherited in sex-specific fashion. The eastern yellow robin (EYR) is an endemic Australian songbird inferred to have experienced climate-driven sex-biased selection and is a prominent model for studying mitochondrial-nuclear interactions in the wild. However, the lack of an EYR reference genome containing both sex chromosomes (in birds, a female bearing Z and W chromosomes) limits efforts to understand the mechanisms of these processes. Here, we assemble the genome for a female EYR and use low-depth (10×) genome resequencing data from 19 individuals of known sex to identify chromosome fragments with sex-specific inheritance. FINDINGS: MaSuRCA hybrid assembly using Nanopore and Illumina reads generated a 1.22-Gb EYR genome in 20,702 scaffolds (94.2% BUSCO completeness). Scaffolds were tested for W-linked (female-only) inheritance using a k-mer approach, and for Z-linked inheritance using median read-depth test in male and female reads (read-depths must indicate haploid female and diploid male representation). This resulted in 2,372 W-linked scaffolds (total length: 97,872,282 bp, N50: 81,931 bp) and 586 Z-linked scaffolds (total length: 121,817,358 bp, N50: 551,641 bp). Anchoring of the sex-linked EYR scaffolds to the reference genome of a female zebra finch revealed 2 categories of sex-linked genomic regions. First, 653 W-linked scaffolds (25.7 Mb) were anchored to the W sex chromosome and 215 Z-linked scaffolds (74.4 Mb) to the Z. Second, 1,138 W-linked scaffolds (70.9 Mb) and 179 Z-linked scaffolds (51.0 Mb) were anchored to a large section (coordinates ∼5 to ∼60 Mb) of zebra finch chromosome 1A. The first ∼5 Mb and last ∼14 Mb of the reference chromosome 1A had only autosomally behaving EYR scaffolds mapping to them. CONCLUSIONS: We report a female (W chromosome-containing) EYR genome and provide genomic evidence for a neo-sex (neo-W and neo-Z) chromosome system in the EYR, involving most of a large chromosome (1A) previously only reported to be autosomal in passerines.

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.

Pleistocene-dated biogeographic barriers drove divergence within the Australo-Papuan region in a sex-specific manner: an example in a widespread Australian songbird.

Understanding how environmental change has shaped species evolution can inform predictions of how future climate change might continue to do so. Research of widespread biological systems spanning multiple climates that have been subject to environmental change can yield generalizable inferences about the neutral and adaptive processes driving lineage divergence during periods of environmental change. We contribute to the growing body of multi-locus phylogeographic studies investigating the effect of Pleistocene climate change on species evolution by focusing on a widespread Australo-Papuan songbird with several mitochondrial lineages that diverged during the Pleistocene, the grey shrike-thrush (Colluricincla harmonica). We employed multi-locus phylogenetic, population genetic and coalescent analyses to (1) assess whether nuclear genetic diversity suggests a history congruent with that based on phenotypically defined subspecies ranges, mitochondrial clade boundaries and putative biogeographical barriers, (2) estimate genetic diversity within and genetic differentiation and gene flow among regional populations and (3) estimate population divergence times. The five currently recognized subspecies of grey shrike-thrush are genetically differentiated in nuclear and mitochondrial genomes, but connected by low levels of gene flow. Divergences among these populations are concordant with recognized historical biogeographical barriers and date to the Pleistocene. Discordance in the order of population divergence events based on mitochondrial and nuclear genomes suggests a history of sex-biased gene flow and/or mitochondrial introgression at secondary contacts. This study demonstrates that climate change can impact sexes with different dispersal biology in different ways. Incongruence between population and mitochondrial trees calls for a genome-wide investigation into dispersal, mitochondrial introgression and mitonuclear evolution.

Concordant divergence of mitogenomes and a mitonuclear gene cluster in bird lineages inhabiting different climates.

Metabolic processes in eukaryotic cells depend on interactions between mitochondrial and nuclear gene products (mitonuclear interactions). These interactions could have a direct role in population divergence. Here, we study mitonuclear co-evolution in a widespread bird that experienced population divergence followed by bidirectional mitochondrial introgression into different nuclear backgrounds. Using >60,000 single nucleotide polymorphisms, we quantify patterns of nuclear genetic differentiation between populations that occupy areas with different climates and harbour deeply divergent mitochondrial lineages despite ongoing nuclear gene flow. We find that strong genetic differentiation and sequence divergence in a region of ~15.4 megabases on chromosome 1A mirror the geographic pattern of mitochondrial DNA divergence. This result is seen in two different transects representing populations with different nuclear backgrounds. The chromosome 1A region is enriched for genes performing mitochondrial functions (N-mt genes). Molecular signatures of selective sweeps in this region alongside those in the mitochondrial genome suggest a history of adaptive mitonuclear co-introgression. Moreover, evidence for large linkage disequilibrium blocks in this genomic region suggests that low recombination could facilitate functional interactions between co-evolved nuclear alleles. Our results are consistent with mitonuclear co-evolution as an important mechanism for population divergence and local adaptation.

Climate-driven mitochondrial selection: A test in Australian songbirds.

Diversifying selection between populations that inhabit different environments can promote lineage divergence within species and ultimately drive speciation. The mitochondrial genome (mitogenome) encodes essential proteins of the oxidative phosphorylation (OXPHOS) system and can be a strong target for climate-driven selection (i.e., associated with inhabiting different climates). We investigated whether Pleistocene climate changes drove mitochondrial selection and evolution within Australian birds. First, using phylogeographic analyses of the mitochondrial ND2 gene for 17 songbird species, we identified mitochondrial clades (mitolineages). Second, using distance-based redundancy analyses, we tested whether climate predicts variation in intraspecific genetic divergence beyond that explained by geographic distances and geographic position. Third, we analysed 41 complete mitogenome sequences representing each mitolineage of 17 species using codon models in a phylogenetic framework and a biochemical approach to identify signals of selection on OXPHOS protein-coding genes and test for parallel selection in mitolineages of different species existing in similar climates. Of 17 species examined, 13 had multiple mitolineages (range: 2-6). Climate was a significant predictor of mitochondrial variation in eight species. At least two amino acid replacements in OXPHOS complex I could have evolved under positive selection in specific mitolineages of two species. Protein homology modelling showed one of these to be in the loop region of the ND6 protein channel and the other in the functionally critical helix HL region of ND5. These findings call for direct tests of the functional and evolutionary significance of mitochondrial protein candidates for climate-associated selection.

Population mitogenomics provides insights into evolutionary history, source of invasions and diversifying selection in the House Crow (Corvus splendens).

The House Crow (Corvus splendens) is a useful study system for investigating the genetic basis of adaptations underpinning successful range expansion. The species originates from the Indian subcontinent, but has successfully spread through a variety of thermal environments across Asia, Africa and Europe. Here, population mitogenomics was used to investigate the colonisation history and to test for signals of molecular selection on the mitochondrial genome. We sequenced the mitogenomes of 89 House Crows spanning four native and five invasive populations. A Bayesian dated phylogeny, based on the 13 mitochondrial protein-coding genes, supports a mid-Pleistocene (~630,000 years ago) divergence between the most distant genetic lineages. Phylogeographic patterns suggest that northern South Asia is the likely centre of origin for the species. Codon-based analyses of selection and assessments of changes in amino acid properties provide evidence of positive selection on the ND2 and ND5 genes against a background of purifying selection across the mitogenome. Protein homology modelling suggests that four amino acid substitutions inferred to be under positive selection may modulate coupling efficiency and proton translocation mediated by OXPHOS complex I. The identified substitutions are found within native House Crow lineages and ecological niche modelling predicts suitable climatic areas for the establishment of crow populations within the invasive range. Mitogenomic patterns in the invasive range of the species are more strongly associated with introduction history than climate. We speculate that invasions of the House Crow have been facilitated by standing genetic variation that accumulated due to diversifying selection within the native range.

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.

Perpendicular axes of differentiation generated by mitochondrial introgression.

Differential introgression of mitochondrial vs. nuclear DNA generates discordant patterns of geographic variation and can promote population divergence and speciation. We examined a potential case of mitochondrial introgression leading to two perpendicular axes of differentiation. The Eastern Yellow Robin Eopsaltria australis, a widespread Australian bird, shows a deep mitochondrial split that is perpendicular to north-south nuclear DNA and plumage colour differentiation. We propose a scenario to explain this pattern: (i) first, both nuclear and mitochondrial genomes differentiated in concert during north-south population divergence; (ii) later, their histories disconnected after two mitochondrial introgression events resulting in a deep mitochondrial split perpendicular to the nuclear DNA structure. We explored this scenario by coalescent modelling of ten mitochondrial genes and 400 nuclear DNA loci. Initial mitochondrial and nuclear genome divergences were estimated to have occurred in the early Pleistocene, consistent with the proposed scenario. Subsequent climatic transitions may have driven later mitochondrial introgression. We consider neutral introgression unlikely and instead propose that the evidence is more consistent with adaptive mitochondrial introgression and selection against incompatible mitochondrial-nuclear combinations. This likely generated an axis of coastal-inland mitochondrial differentiation in the face of nuclear gene flow, perpendicular to the initial north-south axis of differentiation (reflected in genomewide nuclear DNA and colour variation).