MyMedR

Displaying all 9 publications

Abstract:

Sort:

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

Harrisson K, Pavlova A, Gan HM, Lee YP, Austin CM, Sunnucks P

Heredity (Edinb), 2016 Jun;116(6):506-15.
PMID: 26883183 DOI: 10.1038/hdy.2016.8

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.
Fulltext Purifying selection and genetic drift shaped Pleistocene evolution of the mitochondrial genome in an endangered Australian freshwater fish

Pavlova A, Gan HM, Lee YP, Austin CM, Gilligan DM, Lintermans M, et al.

Heredity (Edinb), 2017 05;118(5):466-476.
PMID: 28051058 DOI: 10.1038/hdy.2016.120

Genetic variation in mitochondrial genes could underlie metabolic adaptations because mitochondrially encoded proteins are directly involved in a pathway supplying energy to metabolism. Macquarie perch from river basins exposed to different climates differ in size and growth rate, suggesting potential presence of adaptive metabolic differences. We used complete mitochondrial genome sequences to build a phylogeny, estimate lineage divergence times and identify signatures of purifying and positive selection acting on mitochondrial genes for 25 Macquarie perch from three basins: Murray-Darling Basin (MDB), Hawkesbury-Nepean Basin (HNB) and Shoalhaven Basin (SB). Phylogenetic analysis resolved basin-level clades, supporting incipient speciation previously inferred from differentiation in allozymes, microsatellites and mitochondrial control region. The estimated time of lineage divergence suggested an early- to mid-Pleistocene split between SB and the common ancestor of HNB+MDB, followed by mid-to-late Pleistocene splitting between HNB and MDB. These divergence estimates are more recent than previous ones. Our analyses suggested that evolutionary drivers differed between inland MDB and coastal HNB. In the cooler and more climatically variable MDB, mitogenomes evolved under strong purifying selection, whereas in the warmer and more climatically stable HNB, purifying selection was relaxed. Evidence for relaxed selection in the HNB includes elevated transfer RNA and 16S ribosomal RNA polymorphism, presence of potentially mildly deleterious mutations and a codon (ATP6113) displaying signatures of positive selection (ratio of nonsynonymous to synonymous substitution rates (dN/dS) >1, radical change of an amino-acid property and phylogenetic conservation across the Percichthyidae). In addition, the difference could be because of stronger genetic drift in the smaller and historically more subdivided HNB with low per-population effective population sizes.
Labile sex chromosomes in the Australian freshwater fish family Percichthyidae

Pavlova A, Harrisson KA, Turakulov R, Lee YP, Ingram BA, Gilligan D, et al.

Mol Ecol Resour, 2021 Dec 04.
PMID: 34863023 DOI: 10.1111/1755-0998.13569

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.
Climate-driven mitochondrial selection: A test in Australian songbirds

Lamb AM, Gan HM, Greening C, Joseph L, Lee YP, Morán-Ordóñez A, et al.

Mol Ecol, 2018 02;27(4):898-918.
PMID: 29334409 DOI: 10.1111/mec.14488

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.
Fulltext De novo genome assembly and annotation of Australia's largest freshwater fish, the Murray cod (Maccullochella peelii), from Illumina and Nanopore sequencing read

Austin CM, Tan MH, Harrisson KA, Lee YP, Croft LJ, Sunnucks P, et al.

Gigascience, 2017 08 01;6(8):1-6.
PMID: 28873963 DOI: 10.1093/gigascience/gix063

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.
Fulltext Population mitogenomics provides insights into evolutionary history, source of invasions and diversifying selection in the House Crow (Corvus splendens)

Krzemińska U, Morales HE, Greening C, Nyári ÁS, Wilson R, Song BK, et al.

Heredity (Edinb), 2018 04;120(4):296-309.
PMID: 29180719 DOI: 10.1038/s41437-017-0020-7

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.
Accelerated differentiation of neo-W nuclear-encoded mitochondrial genes between two climate-associated bird lineages signals potential co-evolution with mitogenomes

Low GW, Pavlova A, Gan HM, Ko MC, Sadanandan KR, Lee YP, et al.

Heredity (Edinb), 2024 Aug 22.
PMID: 39174672 DOI: 10.1038/s41437-024-00718-w

There is considerable evidence for mitochondrial-nuclear co-adaptation as a key evolutionary driver. Hypotheses regarding the roles of sex-linkage have emphasized Z-linked nuclear genes with mitochondrial function (N-mt genes), whereas it remains contentious whether the perfect co-inheritance of W genes with mitogenomes could hinder or facilitate co-adaptation. Young (neo-) sex chromosomes that possess relatively many N-mt genes compared to older chromosomes provide unprecedented hypothesis-testing opportunities. Eastern Yellow Robin (EYR) lineages in coastal and inland habitats with different climates are diverged in mitogenomes, and in a ~ 15.4 Mb nuclear region enriched with N-mt genes, in contrast with otherwise-similar nuclear genomes. This nuclear region maps to passerine chromosome 1A, previously found to be neo-sex in the inland EYR genome. To compare sex-linked Chr1A-derived genes between lineages, we assembled and annotated the coastal EYR genome. We found that: (i) the coastal lineage shares a similar neo-sex system with the inland lineage, (ii) neo-W and neo-Z N-mt genes are not more diverged between lineages than are comparable non-N-mt genes, and showed little evidence for broad positive selection, (iii) however, W-linked N-mt genes are more diverged between lineages than are their Z-linked gametologs. The latter effect was ~7 times stronger for N-mt than non-N-mt genes, suggesting that W-linked N-mt genes might have diverged between lineages under environmental selection through co-evolution with mitogenomes. Finally, we identify a candidate gene driver for divergent selection, NDUFA12. Our data represent a rare example suggesting a possible role for W-associated mitochondrial-nuclear interactions in climate-associated adaptation and lineage differentiation.
Signatures of polygenic adaptation associated with climate across the range of a threatened fish species with high genetic connectivity

Harrisson KA, Amish SJ, Pavlova A, Narum SR, Telonis-Scott M, Rourke ML, et al.

Mol Ecol, 2017 Nov;26(22):6253-6269.
PMID: 28977721 DOI: 10.1111/mec.14368

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.
Fulltext Chromosome-length genome assembly and linkage map of a critically endangered Australian bird: the helmeted honeyeater

Robledo-Ruiz DA, Gan HM, Kaur P, Dudchenko O, Weisz D, Khan R, et al.

Gigascience, 2022 Mar 29;11.
PMID: 35348671 DOI: 10.1093/gigascience/giac025

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.