Lines in the land: a review of evidence for eastern Australia's major biogeographical barriers to closed forest taxa

Speciation across biomes: Rapid diversification with reproductive isolation in the Australian delicate mice

Phylogeographic studies of continental clades, especially when combined with palaeoclimate modelling, provide powerful insight into how environment drives speciation across climatic contexts. Australia, a continent characterized by disparate modern biomes and dynamic climate change, provides diverse opportunity to reconstruct the impact of past and present environments on diversification. Here, we use genomic-scale data (1310 exons and whole mitogenomes from 111 samples) to investigate Pleistocene diversification, cryptic diversity, and secondary contact in the Australian delicate mice (Hydromyini: Pseudomys), a recent radiation spanning almost all Australian environments. Across northern Australia, we find no evidence for introgression between cryptic lineages within Pseudomys delicatulus sensu lato, with palaeoclimate models supporting contraction and expansion of suitable habitat since the last glacial maximum. Despite multiple contact zones, we also find little evidence of introgression at a continental scale, with the exception of a potential hybrid zone in the mesic biome. In the arid zone, combined insights from genetic data and palaeomodels support a recent expansion in the arid specialist P. hermannsburgensis and contraction in the semi-arid P. bolami. In the face of repeated secondary contact, differences in sperm morphology and chromosomal rearrangements are potential mechanisms that maintain species boundaries in these recently diverged species. Additionally, we describe the western delicate mouse as a new species and recommend taxonomic reinstatement of the eastern delicate mouse. Overall, we show that speciation in an evolutionarily young and widespread clade has been driven by environmental change, and potentially maintained by divergence in reproductive morphology and chromosome rearrangements.

Phylogeography, hybrid zones and contemporary species boundaries in the south-eastern Australian smooth frogs (Anura: Myobatrachidae: Geocrinia)

Paleo-climatic fluctuations have driven episodic changes in species distributions, providing opportunities for populations to diverge in isolation and hybridise following secondary contact. Studies of phylogeographic diversity and patterns of gene flow across hybrid zones can provide insight into contemporary species boundaries and help to inform taxonomic and conservation inferences. Here we explore geographic diversity within the acoustically divergent yet morphologically conserved south-eastern Australian smooth frog complex and assess gene flow across a narrow hybrid zone using mitochondrial nucleotide sequences and nuclear genome-wide single nucleotide polymorphisms. Our analyses reveal the presence of an evolutionarily distinct taxon restricted to the Otway Plains and Ranges, Victoria, which forms a narrow (9-30 km wide), spatiotemporally stable (>50 years) hybrid zone with Geocrinia laevis, which we describe herein as a new species.

Where to draw the boundaries? Using landscape genomics to disentangle the scribbly gum species complex

PREMISE

Species delimitation is an integral part of evolution and ecology and is vital in conservation science. However, in some groups, species delimitation is difficult, especially where ancestral relationships inferred from morphological or genetic characters are discordant, possibly due to a complicated demographic history (e.g., recent divergences between lineages). Modern genetic techniques can take into account complex histories to distinguish species at a reasonable cost and are increasingly used in numerous applications. We focus on the scribbly gums, a group of up to five closely related and morphologically similar "species" within the eucalypts.

METHODS

Multiple populations of each recognized scribbly gum species were sampled over a wide region across climates, and genomewide scans were used to resolve species boundaries.

RESULTS

None of the taxa were completely divergent, and there were two genetically distinct entities: the inland distributed Eucalyptus rossii and a coastal conglomerate consisting of four species forming three discernible, but highly admixed groups. Divergence among taxa was likely driven by temporal vicariant processes resulting in partial separation across biogeographic barriers. High interspecific gene flow indicated separated taxa reconnected at different points in time, blurring species boundaries.

CONCLUSIONS

Our results highlight the need for genetic screening when dealing with closely related taxonomic entities, particularly those with modest morphological differences. We show that high-throughput sequencing can be effective at identifying species groupings and processes driving divergence, even in the most taxonomically complex groups, and be used as a standard practice for disentangling species complexes.

Revision of Litoria rothii (Anura: Pelodryadidae) from northern Australia

Litoria rothii is a widespread pelodryadid frog with a charismatic laughing advertisement call, distributed across the Australian Monsoon Tropics and southern New Guinea. Given its large distribution spanning well-known biogeographic barriers, variation in male advertisement calls and the prevalence of unresolved species complexes in the Australian frog fauna, we examine the genetic, morphological and acoustic diversity in the species from across its range. Our analyses reveal the presence of a previously unrecognised species in western parts of the range of L. rothii sensu lato, which we describe herein as a new species. Litoria ridibunda sp. nov. is distinguished from L. rothii on the basis of paraphyly of nuclear gene trees with L. everetti from Indonesia, colour patterns on the posterior thigh and male advertisement calls. Compared to L. rothii, the new species has a less contrasting pattern on the posterior thigh and a male advertisement call with a greater number of notes per call and a greater call duration. In particular, the magnitude of call differences between the species is highest where the ranges of the two species are in proximity in north-western Queensland. Our study further emphasises the undiagnosed diversity that remains in Australian frogs, even in relatively large, charismatic, frequently encountered species that often share human dwellings.

Divergent east-west lineages in an Australian fruit fly, (Bactrocera jarvisi), associated with the Carpentaria Basin divide

Bactrocera jarvisi is an endemic Australian fruit fly species (Diptera: Tephritidae). It occurs commonly across tropical and subtropical coastal Australia, from far-northern Western Australia, across the 'Top End' of the Northern Territory, and then down the Queensland east coast. Across this range, its distribution crosses several well documented biogeographic barriers. In order to better understand factors leading to the divergence of Australian fruit fly lineages, we carried out a population genetic study of B. jarvisi from across its range using genome-wide SNP analysis, utilising adult specimens gained from trapping and fruit rearing. Populations from the Northern Territory (NT) and Western Australia were genetically similar to each other, but divergent from the genetically uniform east-coast (= Queensland, QLD) population. Phylogenetic analysis demonstrated that the NT population derived from the QLD population. We infer a role for the Carpentaria Basin as a biogeographic barrier restricting east-west gene flow. The QLD populations were largely panmictic and recognised east-coast biogeographic barriers play no part in north-south population structuring. While the NT and QLD populations were genetically distinct, there was evidence for the historically recent translocation of flies from each region to the other. Flies reared from different host fruits collected in the same location showed no genetic divergence. While a role for the Carpentaria Basin as a barrier to gene flow for Australian fruit flies agrees with existing work on the related B. tryoni, the reason(s) for population panmixia for B. jarvisi (and B. tryoni) over the entire Queensland east coast, a linear north-south distance of >2000km, remains unknown.

Whole-genome phylogeography of the blue-faced honeyeater (Entomyzon cyanotis) and discovery and characterization of a neo-Z chromosome

Whole-genome surveys of genetic diversity and geographic variation often yield unexpected discoveries of novel structural variation, which long-read DNA sequencing can help clarify. Here, we report on whole-genome phylogeography of a bird exhibiting classic vicariant geographies across Australia and New Guinea, the blue-faced honeyeater (Entomyzon cyanotis), and the discovery and characterization of a novel neo-Z chromosome by long-read sequencing. Using short-read genome-wide SNPs, we inferred population divergence events within E. cyanotis across the Carpentarian and other biogeographic barriers during the Pleistocene (~0.3-1.7 Ma). Evidence for introgression between nonsister populations supports a hypothesis of reticulate evolution around a triad of dynamic barriers around Pleistocene Lake Carpentaria between Australia and New Guinea. During this phylogeographic survey, we discovered a large (134 Mbp) neo-Z chromosome and we explored its diversity, divergence and introgression landscape. We show that, as in some sylvioid passerine birds, a fusion occurred between chromosome 5 and the Z chromosome to form a neo-Z chromosome; and in E. cyanotis, the ancestral pseudoautosomal region (PAR) appears nonrecombinant between Z and W, along with most of the fused chromosome 5. The added recombination-suppressed portion of the neo-Z (~37.2 Mbp) displays reduced diversity and faster population genetic differentiation compared with the ancestral-Z. Yet, the new PAR (~17.4 Mbp) shows elevated diversity and reduced differentiation compared to autosomes, potentially resulting from introgression. In our case, long-read sequencing helped clarify the genomic landscape of population divergence on autosomes and sex chromosomes in a species where prior knowledge of genome structure was still incomplete.

Multilocus phylogeography, population genetics and niche evolution of Australian brown and black-tailed treecreepers (Aves: Climacteris)

The Carpentarian barrier across north-eastern Australia is a major biogeographic barrier and a generator of biodiversity within the Australian Monsoonal Tropics. Here we present a continent-wide analysis of mitochondrial (control region) and autosomal (14 anonymous loci) sequence and indel variation and niche modelling of brown and black-tailed treecreepers (Climacteris picumnus and Climacteris melanurus), a clade with a classic distribution on either side of the Carpentarian barrier. mtDNA control region sequences exhibited reciprocal monophyly and strong differentiation (Fst = 0.91), and revealed a signature of a recent selective sweep in C. picumnus. A variety of tests support an isolation-with-migration model of divergence, albeit with low levels of gene flow across the Carpentarian barrier and a divergence time between species of ~1.7–2.8 Mya. Palaeoecological niche models show that both range size as measured by available habitat and estimated historical population sizes of both species declined in the past ~600 kyr and that the area of interspecific range overlap was never historically large, perhaps decreasing opportunities for extensive gene flow. The relatively long divergence time and low opportunity for gene flow may have facilitated speciation more so than in other co-distributed bird taxa across the Australian Monsoonal Tropics.

A box on the river: The phylogenetics and phylogeography of Eucalyptus baueriana (Eucalyptus sect. Adnataria ser. Heterophloiae)

We present a phylogeographic study of the tree species Eucalyptus baueriana Schauer, which occurs in disjunct areas on the near coastal plains and ranges of the south-east Australian mainland. DArTseq data are used to build a phylogeny including E. baueriana and closely related taxa to test its monophyly, test the genetic distinctness of the three subspecies of E. baueriana, and investigate relationships between its disjunct populations. Additionally, we use population structure analysis to investigate the genetic distinctness of populations, and MaxEnt to investigate the environmental factors potentially influencing the species' distribution. We show E. baueriana is monophyletic and most closely related to three other Blue Box eucalypt species: E. conica H.Deane & Maiden, E. dalveenica T.L.Collins, R.L.Andrew & J.J.Bruhl and E. magnificata L.A.S.Johnson & K.D.Hill, with some evidence for genetic introgression between these taxa. Within E. baueriana, the deepest genetic breaks do not correspond with the subspecies classification as the two geographically restricted subspecies, together with samples of the more widespread E. baueriana subsp. baueriana from west of the Gippsland lowlands, form a south-western clade with that is sister to other populations of subsp. baueriana. The oldest genetic break in the species occurs in far eastern Gippsland (Victoria), corresponding to one of the shortest geographic disjunctions in the species' distribution. Genetic breaks in other species have been observed in this region which is broadly referred to as the southern transition zone. Both total annual rainfall and the seasonality of this rainfall are hypothesised to affect the species' distribution; gaps in its distribution are in areas of higher rainfall that support closed forest and in regions with more winter dominated rainfall.

A turn in species conservation for hairpin banksias: demonstration of oversplitting leads to better management of diversity

PREMISE

Understanding evolutionary history and classifying discrete units of organisms remain overwhelming tasks, and lags in this workload concomitantly impede an accurate documentation of biodiversity and conservation management. Rapid advances and improved accessibility of sensitive high-throughput sequencing tools are fortunately quickening the resolution of morphological complexes and   thereby improving the estimation of species diversity. The recently described and critically endangered Banksia vincentia is morphologically similar to the hairpin banksia complex (B. spinulosa s.l.), a group of eastern Australian flowering shrubs whose continuum of morphological diversity has been responsible for taxonomic controversy and possibly questionable conservation initiatives.

METHODS

To assist conservation while testing the current taxonomy of this group, we used high-throughput sequencing to infer a population-scale evolutionary scenario for a sample set that is comprehensive in its representation of morphological diversity and a 2500-km distribution.

RESULTS

Banksia spinulosa s.l. represents two clades, each with an internal genetic structure shaped through historical separation by biogeographic barriers. This structure conflicts with the existing taxonomy for the group. Corroboration between phylogeny and population statistics aligns with the hypothesis that B. collina, B. neoanglica, and B. vincentia should not be classified as species.

CONCLUSIONS

The pattern here supports how morphological diversity can be indicative of a locally expressed suite of traits rather than relationship. Oversplitting in the hairpin banksias is atypical since genomic analyses often reveal that species diversity is underestimated. However, we show that erring on overestimation can yield negative consequences, such as the disproportionate prioritization of a geographically anomalous population.

Taxonomic Revision of the Genus Miltotranes Zimmerman, 1994 (Coleoptera: Curculionidae: Molytinae), the Bowenia-Pollinating Cycad Weevils in Australia, with Description of a New Species and Implications for the Systematics of Bowenia

The Australian endemic weevils of the genus Miltotranes Zimmerman, 1994 (Curculionidae: Molytinae: Tranes group), comprising two species, M. prosternalis (Lea, 1929) and M. subopacus (Lea, 1929), are highly host-specific and the only known pollinators of Bowenia cycads, which comprise two CITES-protected species restricted to Tropical Queensland in Australia. In the present study, the taxonomy of Miltotranes is reviewed, a lectotype for the name Tranes prosternalis Lea, 1929 is designated and a new species associated with the Bowenia population in the McIlwraith Range is described as M. wilsoni sp. n. The descriptions and diagnoses of all species are supplemented with illustrations of their habitus and salient structures, and an identification key to all species and a distribution map are provided. Potential implications of the new species and of the taxonomy and biogeography of Miltotranes overall on the systematics and conservation of Bowenia are discussed.

Future-proofing the koala: synergising genomic and environmental data for effective species management

Climatic and evolutionary processes are inextricably linked to conservation. Avoiding extinction in rapidly changing environments often depends upon a species' capacity to adapt in the face of extreme selective pressures. Here, we employed exon capture and high-throughput next-generation sequencing to investigate the mechanisms underlying population structure and adaptive genetic variation in the koala (Phascolarctos cinereus), an iconic Australian marsupial that represents a unique conservation challenge because it is not uniformly threatened across its range. An examination of 250 specimens representing 91 wild source locations revealed that five major genetic clusters currently exist on a continental scale. The initial divergence of these clusters appears to have been concordant with the Mid-Brunhes Transition (∼ 430-300 kya), a major climatic reorganization that increased the amplitude of Pleistocene glacial-interglacial cycles. While signatures of polygenic selection and environmental adaptation were detected, strong evidence for repeated, climate-associated range contractions and demographic bottleneck events suggests that geographically isolated refugia may have played a more significant role in the survival of the koala through the Pleistocene glaciation than in situ adaptation. Consequently, the conservation of genome-wide genetic variation must be aligned with the protection of core koala habitat to increase the resilience of threatened populations to accelerating anthropogenic threats. Finally, we propose that the five major genetic clusters identified in this study should be accounted for in future koala conservation efforts (e.g. guiding translocations), as existing management divisions in the states of Queensland and New South Wales do not reflect historic or contemporary population structure.

Major biogeographic barriers in eastern Australia have shaped the population structure of widely distributed Eucalyptus moluccana and its putative subspecies

We have investigated the impact of recognized biogeographic barriers on genetic differentiation of grey box (Eucalyptus moluccana), a common and widespread tree species of the family Myrtaceae in eastern Australian woodlands, and its previously proposed four subspecies moluccana, pedicellata, queenslandica, and crassifolia. A range of phylogeographic analyses were conducted to examine the population genetic differentiation and subspecies genetic structure in E. moluccana in relation to biogeographic barriers. Slow evolving markers uncovering long term processes (chloroplast DNA) were used to generate a haplotype network and infer phylogeographic barriers. Additionally, fast evolving, hypervariable markers (microsatellites) were used to estimate demographic processes and genetic structure among five geographic regions (29 populations) across the entire distribution of E. moluccana. Morphological features of seedlings, such as leaf and stem traits, were assessed to evaluate population clusters and test differentiation of the putative subspecies. Haplotype network analysis revealed twenty chloroplast haplotypes with a main haplotype in a central position shared by individuals belonging to the regions containing the four putative subspecies. Microsatellite analysis detected the genetic structure between Queensland (QLD) and New South Wales (NSW) populations, consistent with the McPherson Range barrier, an east-west spur of the Great Dividing Range. The substructure was detected within QLD and NSW in line with other barriers in eastern Australia. The morphological analyses supported differentiation between QLD and NSW populations, with no difference within QLD, yet some differentiation within NSW populations. Our molecular and morphological analyses provide evidence that several geographic barriers in eastern Australia, including the Burdekin Gap and the McPherson Range have contributed to the genetic structure of E. moluccana. Genetic differentiation among E. moluccana populations supports the recognition of some but not all the four previously proposed subspecies, with crassifolia being the most differentiated.

Using 3D geometric morphometrics to aid taxonomic and ecological understanding of a recent speciation event within a small Australian marsupial (Antechinus: Dasyuridae)

Taxonomic distinction of species forms the foundation of biodiversity assessments and conservation priorities. However, traditional morphological and/or genetics-based taxonomic assessments frequently miss the opportunity of elaborating on the ecological and functional context of species diversification. Here, we used 3D geometric morphometrics of the cranium to improve taxonomic differentiation and add ecomorphological characterization of a young cryptic divergence within the carnivorous marsupial genus Antechinus. Specifically, we used 168 museum specimens to characterize the recently proposed clades A. stuartii ‘south’, A. stuartii ‘north’ and A. subtropicus. Beyond slight differences attributable to overall size (and, therefore, not necessarily diagnostic), we also found clear allometry-independent shape variation. This allowed us to define new, easily measured diagnostic traits in the palate, which differentiate the three clades. Contrary to previous suggestions, we found no support for a latitudinal gradient as causing the differentiation between the clades. However, skull shape co-varied with temperature and precipitation seasonality, suggesting that the clades may be adapted to environmental variables that are likely to be impacted by climate change. Our study demonstrates the use of 3D geometric morphometrics to improve taxonomic diagnosis of cryptic mammalian species, while providing perspectives on the adaptive origins and potential future threats of mammalian diversity.

Disparate origins for endemic bird taxa from the ‘Gondwana Rainforests’ of Central Eastern Australia

Subtropical and temperate rainforests of Central Eastern Australia are some of the largest remaining fragments of their kind globally. The biota of these rainforests appears to comprise two broad biogeographical elements: a more ancient (Miocene or older) and typically upland temperate (‘Gondwanan’) element and a younger (Plio-Pleistocene) lowland tropical element. We present the first phylogenetic synthesis of the spatiotemporal origins for the eight bird taxa endemic to Central Eastern Australian Rainforests. At least five of these eight focal taxa show Plio-Pleistocene divergences from their respective northern sister taxa, consistent with origins driven by recent expansion and contraction of lowland rainforest. In contrast, two more strictly upland species, the rufous scrub-bird (Atrichornis rufescens) and the logrunner (Orthonyx temminckii), diverged from their nearest living relatives during the Miocene, suggesting potentially longer histories of persistence and more temperate origins. Finally, we did not recover reciprocal monophyly in mitogenomes from the two extant lyrebirds, Albert’s lyrebird (Menura alberti) and the superb lyrebird (Menura novaehollandiae). The disparate divergence ages recovered among all eight taxa are consistent with the biota of the Central Eastern Australian Rainforests comprising isolates either of younger age and tropical lowland origins or of older age and temperate upland origins.

Benchmarking Taxonomic and Genetic Diversity After the Fact: Lessons Learned From the Catastrophic 2019–2020 Australian Bushfires

Environmental catastrophes are increasing in frequency and severity under climate change, and they substantially impact biodiversity. Recovery actions after catastrophes depend on prior benchmarking of biodiversity and that in turn minimally requires critical assessment of taxonomy and species-level diversity. Long-term recovery of species also requires an understanding of within-species diversity. Australia’s 2019–2020 bushfires were unprecedented in their extent and severity and impacted large portions of habitats that are not adapted to fire. Assessments of the fires’ impacts on vertebrates identified 114 species that were a high priority for management. In response, we compiled explicit information on taxonomic diversity and genetic diversity within fire-impacted vertebrates to provide to government agencies undertaking rapid conservation assessments. Here we discuss what we learned from our effort to benchmark pre-fire taxonomic and genetic diversity after the event. We identified a significant number of candidate species (genetic units that may be undescribed species), particularly in frogs and mammals. Reptiles and mammals also had high levels of intraspecific genetic structure relevant to conservation management. The first challenge was making published genetic data fit for purpose because original publications often focussed on a different question and did not provide raw sequence read data. Gaining access to analytical files and compiling appropriate individual metadata was also time-consuming. For many species, significant unpublished data was held by researchers. Identifying which data existed was challenging. For both published and unpublished data, substantial sampling gaps prevented areas of a species’ distribution being assigned to a conservation unit. Summarising sampling gaps across species revealed that many areas were poorly sampled across taxonomic groups. To resolve these issues and prepare responses to future catastrophes, we recommend that researchers embrace open data principles including providing detailed metadata. Governments need to invest in a skilled taxonomic workforce to document and describe biodiversity before an event and to assess its impacts afterward. Natural history collections should also target increasing their DNA collections based on sampling gaps and revise their collection strategies to increasingly take population-scale DNA samples in order to document within-species genetic diversity.

Koala Retrovirus in Northern Australia Shows a Mixture of Stable Endogenization and Exogenous Lineage Diversification within Fragmented Koala Populations

The koala population in northern Australia has become increasingly fragmented due to natural and man-made barriers and interventions. This situation has created a unique opportunity to study both endogenous and exogenous koala retrovirus (KoRV). To determine the impact that population isolation has had on KoRV diversity in Queensland, 272 koalas from six fragmented koala populations were profiled for their KoRV provirus across two natural biogeographical barriers (the St Lawrence Gap and the Brisbane Valley Barrier), one man-made geographical barrier (the city of Brisbane) and two translocation events (the single movement of koalas to an island and the repeated movement of koalas into a koala sanctuary). Analysis revealed that all koalas tested were KoRV-A positive, with 90 - 96% of the detected KoRV provirus from each koala representing a single, likely endogenous, KoRV-A strain. The next most abundant proviral sequence was a defective variant of the dominant KoRV-A strain, accounting for 3 - 10% of detected provirus. The remaining KoRV provirus represented expected exogenous strains of KoRV and included geographically localized patterns of KoRV-B, -C, -D, -F, -G, and -I. These results indicate that lineage diversification of exogenous KoRV is actively ongoing. In addition, comparison of KoRV provirus within known dam-sire-joey family groups from the koala sanctuary revealed that joeys consistently had KoRV proviral patterns more similar to their dams than their sires in KoRV-B, -C and -D provirus composition. Collectively, this study highlights both the consistency of endogenous KoRV and the diversity of exogenous KoRV across the fragmented koala populations in northern Australia.IMPORTANCE KoRV infection has become a permanent part of koalas in northern Australia. With KoRV presence and abundance linked to more severe chlamydial disease and neoplasia in these koalas, understanding how KoRV exists throughout an increasingly fragmented koala population is a key first step in designing conservation and management strategies. This survey of KoRV provirus in Queensland koalas indicates that endogenous KoRV provirus is ubiquitous and consistent throughout the state while exogenous KoRV provirus is diverse and distinct in fragmented koala populations. Understanding the prevalence and impact of both endogenous and exogenous KoRV will be needed to ensure a future for all koala populations.

Divergent lineages in a semi-arid mallee species, Eucalyptus behriana, correspond to a major geographic break in southeastern Australia

AIM

To infer relationships between populations of the semi-arid, mallee eucalypt, Eucalyptus behriana, to build hypotheses regarding evolution of major disjunctions in the species' distribution and to expand understanding of the biogeographical history of southeastern Australia.

LOCATION

Southeastern Australia.

TAXON

Eucalyptus behriana (Myrtaceae, Angiospermae).

METHODS

We developed a large dataset of anonymous genomic loci for 97 samples from 11 populations of E. behriana using double digest restriction site-associated DNA sequencing (ddRAD-seq), to determine genetic relationships between the populations. These relationships, along with species distribution models, were used to construct hypotheses regarding environmental processes that have driven fragmentation of the species' distribution.

RESULTS

Greatest genetic divergence was between populations on either side of the Lower Murray Basin. Populations west of the Basin showed greater genetic divergence between one another than the eastern populations. The most genetically distinct population in the east (Long Forest) was separated from others by the Great Dividing Range. A close relationship was found between the outlying northernmost population (near West Wyalong) and those in the Victorian Goldfields despite a large disjunction between them.

CONCLUSIONS

Patterns of genetic variation are consistent with a history of vicariant differentiation of disjunct populations. We infer that an early disjunction to develop in the species distribution was that across the Lower Murray Basin, an important biogeographical barrier separating many dry sclerophyll plant taxa in southeastern Australia. Additionally, our results suggest that the western populations fragmented earlier than the eastern ones. Fragmentation, both west and east of the Murray Basin, is likely tied to climatic changes associated with glacial-interglacial cycles although it remains possible that major geological events including uplift of the Mount Lofty Ranges and basalt flows in the Newer Volcanics Province also played a role.

Looks can be deceiving: speciation dynamics of co-distributed Angophora (Myrtaceae) species in a varying landscape

Understanding the mechanisms underlying species divergence remains a central goal in evolutionary biology. Landscape genetics can be a powerful tool for examining evolutionary processes. We used genome-wide scans to genotype samples from populations of eight Angophora species. Angophora is a small genus within the eucalypts comprising common and rare species in a heterogeneous landscape, making it an appropriate group to study speciation. We found A. hispida was highly differentiated from the other species. Two subspecies of A. costata (subsp. costata and subsp. euryphylla) formed a group, while the third (subsp. leiocarpa, which is only distinguished by its smooth fruits and provenance) was supported as a distinct pseudocryptic species. Other species that are morphologically distinct could not be genetically differentiated (e.g., A. floribunda and A. subvelutina). Distribution and genetic differentiation within Angophora were strongly influenced by temperature and humidity, as well as biogeographic barriers, particularly rivers and higher elevation regions. While extensive introgression was found between many populations of some species (e.g., A. bakeri and A. floribunda), others only hybridized at certain locations. Overall, our findings suggest multiple mechanisms drove evolutionary diversification in Angophora and highlight how genome-wide analyses of related species in a diverse landscape can provide insights into speciation.

Koala immunogenetics and chlamydial strain type are more directly involved in chlamydial disease progression in koalas from two south east Queensland koala populations than koala retrovirus subtypes

Chlamydial disease control is increasingly utilised as a management tool to stabilise declining koala populations, and yet we have a limited understanding of the factors that contribute to disease progression. To examine the impact of host and pathogen genetics, we selected two geographically separated south east Queensland koala populations, differentially affected by chlamydial disease, and analysed koala major histocompatibility complex (MHC) genes, circulating strains of Chlamydia pecorum and koala retrovirus (KoRV) subtypes in longitudinally sampled, well-defined clinical groups. We found that koala immunogenetics and chlamydial genotypes differed between the populations. Disease progression was associated with specific MHC alleles, and we identified two putative susceptibility (DCb 03, DBb 04) and protective (DAb 10, UC 01:01) variants. Chlamydial genotypes belonging to both Multi-Locus Sequence Typing sequence type (ST) 69 and ompA genotype F were associated with disease progression, whereas ST 281 was associated with the absence of disease. We also detected different ompA genotypes, but not different STs, when long-term infections were monitored over time. By comparison, KoRV profiles were not significantly associated with disease progression. These findings suggest that chlamydial genotypes vary in pathogenicity and that koala immunogenetics and chlamydial strains are more directly involved in disease progression than KoRV subtypes.

Integrative taxonomic investigation of Petaurus breviceps (Marsupialia: Petauridae) reveals three distinct species

The Australian sugar glider, Petaurus breviceps s.l., is widely distributed across eastern and northern Australia. Examination of historical and contemporary collections of Petaurus specimens and phylogenetic analyses have revealed considerable taxonomic diversity within the genus. We aimed to utilize an integrative taxonomic approach, combining genetic and morphological evidence, to resolve the taxonomy of Australian gliders currently recognized as Petaurus breviceps. Herein, we confirm the existence of three distinct species: P. breviceps, P. notatus comb. nov. and P. ariel comb. nov.. Petaurus breviceps and P. notatus are each represented by major mtDNA lineages in P. breviceps, while P. ariel forms a sister-lineage to P. norfolcensis and P. gracilis. Subtle morphological differences distinguish P. breviceps from the closely related P. notatus, while the morphological distinctions between P. ariel and its genetically similar sister-taxa, P. norfolcensis and P. gracilis, are more obvious. Given the purported broad geographic distribution of the taxon, P. breviceps s.l. was not listed as threatened, but dividing this taxon into three species has important conservation implications for all taxa in the group, particularly given the lamentable record for mammal extinctions in Australia. Concerted and targeted conservation efforts are necessary to preserve these distinct, newly described species.

Low-coverage genomic data resolve the population divergence and gene flow history of an Australian rain forest fig wasp

Population divergence and gene flow are key processes in evolution and ecology. Model-based analysis of genome-wide data sets allows discrimination between alternative scenarios for these processes even in nonmodel taxa. We used two complementary approaches (one based on the blockwise site frequency spectrum [bSFS], the second on the pairwise sequentially Markovian coalescent [PSMC]) to infer the divergence history of a fig wasp, Pleistodontes nigriventris. Pleistodontes nigriventris and its fig tree mutualist Ficus watkinsiana are restricted to rain forest patches along the eastern coast of Australia and are separated into The Northern population is to the north of the Southern populations by two dry forest corridors (the Burdekin and St. Lawrence Gaps). We generated whole genome sequence data for two haploid males per population and used the bSFS approach to infer the timing of divergence between northern and southern populations of P. nigriventris, and to discriminate between alternative isolation with migration (IM) and instantaneous admixture (ADM) models of postdivergence gene flow. Pleistodontes nigriventris has low genetic diversity (π = 0.0008), to our knowledge one of the lowest estimates reported for a sexually reproducing arthropod. We find strongest support for an ADM model in which the two populations diverged ca. 196 kya in the late Pleistocene, with almost 25% of northern lineages introduced from the south during an admixture event ca. 57 kya. This divergence history is highly concordant with individual population demographies inferred from each pair of haploid males using PSMC. Our analysis illustrates the inferences possible with genome-level data for small population samples of tiny, nonmodel organisms and adds to a growing body of knowledge on the population structure of Australian rain forest taxa.

Understanding Historical Demographic Processes to Inform Contemporary Conservation of an Arid zone Specialist: The Yellow-Footed Rock-Wallaby

Little genetic research has been undertaken on mammals across the vast expanse of the arid biome in Australia, despite continuing species decline and need for conservation management. Here, we evaluate the contemporary and historical genetic connectivity of the yellow-footed rock-wallaby, Petrogalexanthopusxanthopus, a threatened macropodid which inhabits rocky outcrops across the disconnected mountain range systems of the southern arid biome. We use 17 microsatellite loci together with mitochondrial control region data to determine the genetic diversity of populations and the evolutionary processes shaping contemporary population dynamics on which to base conservation recommendations. Our results indicate the highly fragmented populations have reduced diversity and limited contemporary gene flow, with most populations having been through population bottlenecks. Despite limited contemporary gene flow, the phylogeographic relationships of the mitochondrial control region indicate a lack of structure and suggests greater historical connectivity. This is an emerging outcome for mammals across this arid region. On the basis of our results, we recommend augmentation of populations of P. x.xanthopus, mixing populations from disjunct mountain range systems to reduce the chance of continued diversity loss and inbreeding depression, and therefore maximize the potential for populations to adapt and survive into the future.

How ‘Gondwanan’ is Riethia? Molecular phylogenetics elucidates the mode and tempo of diversification in Austro-Pacific Chironominae (Diptera)

Riethia Kieffer, a genus of the non-biting midge subfamily Chironominae (Diptera: Chironomidae) is distributed in Australia, New Zealand, New Caledonia and South America. This austral distribution could be due to earth history (vicariance) or from Southern Hemisphere dispersal(s). We obtained samples from each area, most intensively from throughout Australia. We included putative sister genus Pseudochironomus Malloch, many genera from tribe Tanytarsini, enigmatic taxa in Chironomini and conventional outgroups from other subfamilies. We assembled a multilocus molecular dataset for four genetic regions from 107 individuals to reconstruct the first dated molecular phylogeny for the group. Four terminal clusters corresponded to unreared (thus unassociated) larvae. Monophyly was supported for ‘core’ Riethia, Pseudochironomus, putative tribe Pseudochironomini, tribe Tanytarsini (including enigmatic Nandeva Wiedenbrug, Reiss & Fittkau) and subfamily Chironominae. All species are monophyletic except for R. cinctipes Freeman, which includes R. neocaledonica Cranston. Riethia zeylandica Freeman, previously thought to be widespread in eastern Australia, now is a New Zealand endemic with Australian specimens allocated now to several regionally restricted species. The origin of Riethia was at 60.6 Ma (‘core’) or 52.1 Ma depending on the relationship of two South American species. Both dates are before the break-up of South America and Australia. Diversification within crown group Riethia started before the Cretaceous–Paleogene boundary, with subsequent separation at 52 Ma of an Australian ‘clade I’ from its sister ‘clade II’, which comprises Australian, New Zealand and New Caledonian species. Inferred dates for species origins of New Caledonia and New Zealand taxa imply transoceanic dispersals from eastern Australia. Western Australian species diverged during the mid to late Miocene from their eastern Australian sister taxa. This correlates with the onset of drying of Australia and the separation of mesic east from west by the formation of an arid proto-Nullarbor. Taken together, the inferred tempo of diversification in the group included both older ages reflecting earth history, yet with suggested recent intra-Pacific separations due to transoceanic dispersals.

First phylogenetic assessment and taxonomic synopsis of the open-holed trapdoor spider genus Namea (Mygalomorphae: Anamidae): a highly diverse mygalomorph lineage from Australia’s tropical eastern rainforests

The tropical and subtropical rainforests of Australia’s eastern mesic zone have given rise to a complex and highly diverse biota. Numerous old endemic, niche-conserved groups persist in the montane rainforests south of Cooktown, where concepts of serial allopatric speciation resulting from the formation of xeric interzones have largely driven our biogeographic understanding of the region. Among invertebrate taxa, studies on less vagile arachnid lineages now complement extensive research on vertebrate taxa, and phylogenetic studies on mygalomorph spiders in particular are revealing significant insights about the biogeographic history of the Australian continent since the Eocene. One mygalomorph lineage entirely endemic to Australia’s tropical and subtropical eastern rainforests is the open-holed trapdoor spider genus Namea Raven, 1984 (family Anamidae). We explore, for the first time, the phylogenetic diversity and systematics of this group of spiders, with the aims of understanding patterns of rainforest diversity in Namea, of exploring the relative roles of lineage overlap versus in situ speciation in driving predicted high levels of congeneric sympatry, and of broadly reconciling morphology with evolutionary history. Original and legacy sequences were obtained for three mtDNA and four nuDNA markers from 151 specimens, including 82 specimens of Namea. We recovered a monophyletic genus Namea sister to the genus Teyl Main, 1975, and monophyletic species clades corresponding to 30 morphospecies OTUs, including 22 OTUs nested within three main species-complex lineages. Remarkable levels of sympatry for a single genus of mygalomorph spiders were revealed in rainforest habitats, with upland subtropical rainforests in south-eastern Queensland often home to multiple (up to six) congeners of usually disparate phylogenetic affinity living in direct sympatry or close parapatry, likely the result of simultaneous allopatric speciation in already co-occurring lineages, and more recent dispersal in a minority of taxa. In situ speciation, in contrast, appears to have played a relatively minor role in generating sympatric diversity within rainforest ‘islands’. At the population level, changes in the shape and spination of the male first leg relative to evolutionary history reveal subtle but consistent interspecific morphological shifts in the context of otherwise intraspecific variation, and understanding this morphological variance provides a useful framework for future taxonomic monography. Based on the phylogenetic results, we further provide a detailed taxonomic synopsis of the genus Namea, formally diagnosing three main species-complexes (the brisbanensis-complex, the dahmsi-complex and the jimna-complex), re-illustrating males of all 15 described species, and providing images of live spiders and burrows where available. In doing so, we reveal a huge undescribed diversity of Namea species from tropical and subtropical rainforest habitats, and an old endemic fauna that is beginning to shed light on more complex patterns of rainforest biogeography.

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.

Wild Origins of Macadamia Domestication Identified Through Intraspecific Chloroplast Genome Sequencing

Identifying the geographic origins of crops is important for the conservation and utilization of novel genetic variation. Even so, the origins of many food crops remain elusive. The tree nut crop macadamia has a remarkable domestication history, from subtropical rain forests in Australia through Hawaii to global cultivation all within the last century. The industry is based primarily on Macadamia integrifolia and M. integrifolia-M. tetraphylla hybrid cultivars with Hawaiian cultivars the main contributors to world production. Sequence data from the chloroplast genome assembled using a genome skimming strategy was used to determine population structure among remnant populations of the main progenitor species, M. integrifolia. Phylogenetic analysis of a 506 bp chloroplast SNP alignment from 64 wild and cultivated accessions identified phylogeographic structure and deep divergences between clades providing evidence for historical barriers to seed dispersal. High levels of variation were detected among wild accessions. Most Hawaiian cultivars, however, shared a single chlorotype that was also present at two wild sites at Mooloo and Mt Bauple from the northernmost distribution of the species in south-east Queensland. Our results provide evidence for a maternal genetic bottleneck during early macadamia domestication, and pinpoint the likely source of seed used to develop the Hawaiian cultivars. The extensive variability and structuring of M. integrifolia chloroplast genomic variation detected in this study suggests much unexploited genetic diversity is available for improvement of this recently domesticated crop.

Miocene biome turnover drove conservative body size evolution across Australian vertebrates

On deep time scales, changing climatic trends can have a predictable influence on macroevolution. From evidence of mass extinctions, we know that rapid climatic oscillations can indirectly open niche space and precipitate adaptive radiation, changing the course of ecological diversification. These dramatic shifts in the global climate, however, are rare events relative to extended periods of protracted climate change and biome turnover. It remains unclear whether during gradually changing periods, shifting habitats may instead promote non-adaptive speciation by facilitating allopatry and phenotypic conservatism. Using fossil-calibrated, species-level phylogenies for five Australian radiations comprising more than 800 species, we investigated temporal trends in biogeography and body size evolution. Here, we demonstrate that gradual Miocene cooling and aridification correlates with the restricted phenotypic diversification of multiple ecologically diverse vertebrate groups. This probably occurred as species ranges became fractured and isolated during continental biome restructuring, encouraging a shift towards conservatism in body size evolution. Our results provide further evidence that abiotic changes, not only biotic interactions, may act as selective forces influencing phenotypic macroevolution.

Dissection by genomic and plumage variation of a geographically complex hybrid zone between two Australian non-sister parrot species, Platycercus adscitus and Platycercus eximius

The study of hybrid zones advances understanding of the speciation process, and approaches incorporating genomic data are increasingly used to draw significant conclusions about the impact of hybridisation. Despite the progress made, the complex interplay of factors that can lead to substantially variable hybridisation outcomes are still not well understood, and many systems and/or groups remain comparatively poorly studied. Our study aims to broaden the literature on avian hybrid zones, investigating a potentially geographically and temporally complex putative hybrid zone between two native Australian non-sister parrot species, the pale-headed and eastern rosellas (Platycercus adscitus and Platycercus eximius, respectively). We analysed six plumage traits and >1400 RADseq loci and detected hybrid individuals and an unexpectedly complex geographic structure. The hybrid zone is larger than previously described due to either observer bias or its movement over recent decades. It comprises different subregions where genetic and plumage signals of admixture vary markedly in their concordance. Evidence of contemporary hybridisation (later generation and backcrossed individuals) both within and beyond the previously defined zone, when coupled with a lack of F1 hybrids and differential patterns of introgression among potentially diagnostic loci, indicates a lack of post-zygotic barriers to gene flow between species. Despite ongoing gene flow, species boundaries are likely maintained largely by strong pre-mating barriers. These findings are discussed in detail and future avenues for research into this system are proposed, which would be of benefit to the speciation and hybrid zone literature.

Insights into Platypus Population Structure and History from Whole-Genome Sequencing

The platypus is an egg-laying mammal which, alongside the echidna, occupies a unique place in the mammalian phylogenetic tree. Despite widespread interest in its unusual biology, little is known about its population structure or recent evolutionary history. To provide new insights into the dispersal and demographic history of this iconic species, we sequenced the genomes of 57 platypuses from across the whole species range in eastern mainland Australia and Tasmania. Using a highly improved reference genome, we called over 6.7 M SNPs, providing an informative genetic data set for population analyses. Our results show very strong population structure in the platypus, with our sampling locations corresponding to discrete groupings between which there is no evidence for recent gene flow. Genome-wide data allowed us to establish that 28 of the 57 sampled individuals had at least a third-degree relative among other samples from the same river, often taken at different times. Taking advantage of a sampled family quartet, we estimated the de novo mutation rate in the platypus at 7.0 × 10-9/bp/generation (95% CI 4.1 × 10-9-1.2 × 10-8/bp/generation). We estimated effective population sizes of ancestral populations and haplotype sharing between current groupings, and found evidence for bottlenecks and long-term population decline in multiple regions, and early divergence between populations in different regions. This study demonstrates the power of whole-genome sequencing for studying natural populations of an evolutionarily important species.

Mass turnover and recovery dynamics of a diverse Australian continental radiation

Trends in global and local climate history have been linked to observed macroevolutionary patterns across a variety of organisms. These climatic pressures may unilaterally or asymmetrically influence the evolutionary trajectory of clades. To test and compare signatures of changing global (Eocene-Oligocene boundary cooling) and continental (Miocene aridification) environments on a continental fauna, we investigated the macroevolutionary dynamics of one of Australia's most diverse endemic radiations, pygopodoid geckos. We generated a time-calibrated phylogeny (>90% taxon coverage) to test whether (i) asymmetrical pygopodoid tree shape may be the result of mass turnover deep in the group's history, and (ii) how Miocene aridification shaped trends in biome assemblages. We find evidence of mass turnover in pygopodoids following the isolation of the Australian continental plate ∼30 million years ago, and in contrast, gradual aridification is linked to elevated speciation rates in the young arid zone. Surprisingly, our results suggest that invasion of arid habitats was not an evolutionary end point. Instead, arid Australia has acted as a source for diversity, with repeated outward dispersals having facilitated diversification of this group. This pattern contrasts trends in richness and distribution of other Australian vertebrates, illustrating the profound effects historical biome changes have on macroevolutionary patterns.