Plio-Pleistocene diversification and biogeographic barriers in southern Australia reflected in the phylogeography of a widespread and common lizard species.

Palaeoclimatic events and biogeographical processes since the mid-Tertiary have played an important role in shaping the evolution and distribution of Australian fauna. However, their impacts on fauna in southern and arid zone regions of Australia are not well understood. Here we investigate the phylogeography of an Australian scincid lizard, Tiliqua rugosa, across southern Australia using mitochondrial DNA (mtDNA) and 11 nuclear DNA markers (nuDNA), including nine anonymous nuclear loci. Phylogenetic analyses revealed three major mtDNA lineages within T. rugosa, geographically localised north and south of the Murray River in southern Australia, and west of the Nullarbor Plain. Molecular variance and population analyses of both mtDNA and nuDNA haplotypes revealed significant variation among the three populations, although potential introgression of nuDNA markers was also detected for the Northern and Southern population. Coalescent times for major mtDNA lineages coincide with an aridification phase, which commenced after the early Pliocene and increased in intensity during the Late Pliocene-Pleistocene. Species distribution modelling and a phylogeographic diffusion model suggest that the range of T. rugosa may have contracted during the Last Glacial Maximum and the locations of optimal habitat appear to coincide with the geographic origin of several distinct mtDNA lineages. Overall, our analyses suggest that Plio-Pleistocene climatic changes and biogeographic barriers associated with the Nullarbor Plain and Murray River have played a key role in shaping the present-day distribution of genetic diversity in T. rugosa and many additional ground-dwelling animals distributed across southern Australia.

Multi-locus sequence data illuminate demographic drivers of Pleistocene speciation in semi-arid southern Australian birds (Cinclosoma spp.).

BACKGROUND: During the Pleistocene, shifts of species distributions and their isolation in disjunct refugia led to varied outcomes in how taxa diversified. Some species diverged, others did not. Here, we begin to address another facet of the role of the Pleistocene in generating today's diversity. We ask which processes contributed to divergence in semi-arid southern Australian birds. We isolated 11 autosomal nuclear loci and one mitochondrial locus from a total of 29 specimens of the sister species pair, Chestnut Quail-thrush Cinclosoma castanotum and Copperback Quail-thrush C. clarum. RESULTS: A population clustering analysis confirmed the location of the current species boundary as a well-known biogeographical barrier in southern Australia, the Eyrean Barrier. Coalescent-based analyses placed the time of species divergence to the Middle Pleistocene. Gene flow between the species since divergence has been low. The analyses suggest the effective population size of the ancestor was 54 to 178 times smaller than populations since divergence. This contrasts with recent multi-locus studies in some other Australian birds (butcherbirds, ducks) where a lack of phenotypic divergence was accompanied by larger historical population sizes. Post-divergence population size histories of C. clarum and C. castanotum were inferred using the extended Bayesian skyline model. The population size of C. clarum increased substantially during the late Pleistocene and continued to increase through the Last Glacial Maximum and Holocene. The timing of this expansion across its vast range is broadly concordant with that documented in several other Australian birds. In contrast, effective population size of C. castanotum was much more constrained and may reflect its smaller range and more restricted habitat east of the Eyrean Barrier compared with that available to C. clarum to the west. CONCLUSIONS: Our results contribute to awareness of increased population sizes, following significant contractions, as having been important in shaping diversity in Australian arid and semi-arid zones. Further, we improve knowledge of the role of Pleistocene climatic shifts in areas of the planet that were not glaciated at that time but which still experienced that period's cyclical climatic fluctuations.

Evolution of arid zone birds in Australia: leapfrog distribution patterns and mesic-arid connections in quail-thrush (Cinclosoma, Cinclosomatidae).

The quail-thrush, Cinclosoma, include between five and seven species distributed broadly across arid and semi-arid inland Australia, mesic forests of south-eastern Australia and New Guinea. It has been suggested that the arid zone species of quail-thrush arose from forest ancestors as Australia changed from a warm wet climate to a cooler drier climate since the late-Miocene. We generated multilocus (mitochondrial ND2 and eight nuclear loci) gene and species trees with complete taxon sampling of Cinclosoma to investigate evolutionary relationships and species status of some taxa. Topologies reconstructed in congruent, highly-resolved gene trees and species trees that supported the recognition of seven species. Ancestral state reconstruction and divergence time estimates suggest that arid-adapted taxa radiated in parallel with a drying climate and changing habitat. A 'leapfrog' distribution in phenotypes of arid zone taxa was likely a result of ancestral retention of inconspicuous (or camouflaged) plumage patterns. A specimen-based report from 1968 of hybridization between non-sister taxa Cinclosoma castanotum and Cinclosoma marginatum was verified using molecular analysis on specimens collected at the same locality 40 years later. We discuss the implications of hybridization to the evolution of this species group.

Robbery in progress: Historical museum collections bring to light a mitochondrial capture within a bird species widespread across southern Australia, the Copperback Quail-thrush Cinclosoma clarum.

We surveyed mitochondrial, autosomal, and Z chromosome diversity within and between the Copperback Quail-thrush Cinclosoma clarum and Chestnut Quail-thrush C. castanotum, which together span the arid and semi-arid zones of southern Australia, and primarily from specimens held in museum collections. We affirm the recent taxonomic separation of the two species and then focus on diversity within the more widespread of the two species, C. clarum. To guide further study of the system and what it offers to understanding the genomics of the differentiation and speciation processes, we develop and present a hypothesis to explain mitonuclear discordance that emerged in ourdata. Following a period of historical allopatry, secondary contact has resulted in an eastern mitochondrial genome replacing the western mitochondrial genome in western populations. This is predicted under a population-level invasion in the opposite direction, that of the western population invading the range of the eastern one. Mitochondrial captures can be driven by neutral, demographic processes, or adaptive mechanisms, and we favor the hypothesized capture being driven by neutral means. We cannot fully reject the adaptive process but suggest how these alternatives may be further tested. We acknowledge an alternative hypothesis, which finds some support in phenotypic data published elsewhere, namely that outcomes of secondary contact have been more complex than our current genomic data suggest. Discriminating and reconciling these two alternative hypotheses, which may not be mutually exclusive, could be tested with closer sampling at levels of population, individual, and nucleotide than has so far been possible. This would be further aided by knowledge of the genetic basis to phenotypic variation described elsewhere.

Where and when does a ring start and end? Testing the ring-species hypothesis in a species complex of Australian parrots.

Speciation, despite ongoing gene flow can be studied directly in nature in ring species that comprise two reproductively isolated populations connected by a chain or ring of intergrading populations. We applied three tiers of spatio-temporal analysis (phylogeny/historical biogeography, phylogeography and landscape/population genetics) to the data from mitochondrial and nuclear genomes of eastern Australian parrots of the Crimson Rosella Platycercus elegans complex to understand the history and present genetic structure of the ring they have long been considered to form. A ring speciation hypothesis does not explain the patterns we have observed in our data (e.g. multiple genetic discontinuities, discordance in genotypic and phenotypic assignments where terminal differentiates meet). However, we cannot reject that a continuous circular distribution has been involved in the group's history or indeed that one was formed through secondary contact at the 'ring's' east and west; however, we reject a simple ring-species hypothesis as traditionally applied, with secondary contact only at its east. We discuss alternative models involving historical allopatry of populations. We suggest that population expansion shown by population genetics parameters in one of these isolates was accompanied by geographical range expansion, secondary contact and hybridization on the eastern and western sides of the ring. Pleistocene landscape and sea-level and habitat changes then established the birds' current distributions and range disjunctions. Populations now show idiosyncratic patterns of selection and drift. We suggest that selection and drift now drive evolution in different populations within what has been considered the ring.

Combined mitochondrial and nuclear data enhance resolution of a rapid radiation of Australian rainbow skinks (Scincidae: Carlia).

In contrast to low diversity seen in many Australian rainforest squamate genera, sclerophyll adapted groups--such as Carlia--show signs of faster diversification. Here we expand upon a previous single-locus mitochondrial DNA phylogenetic study of Carlia which described a major polytomy at an intermediate level of divergence. With additional mtDNA data, two nuclear intron loci and comprehensive taxonomic coverage, we provide support, congruent across loci, for the existence of three major clades. In doing so we recognise three genera for the 'Carlia group of skinks': clade 1, Carlia Gray, 1845; 2, Lygisaurus De Vis, 1884 (includes all species formerly known as Lygisaurus, and also includes C. parrhasius); 3, Liburnascincus Wells and Wellington, 1984 (consisting of three boulder-dwelling species). Remaining regions of low bootstrap and posterior probability support are associated with short internodes and apparent conflict among loci, as inferred by Partition Branch Support. Likelihood-based diversification-rate analysis rejects constant rate models, and indicates that Carlia underwent a period of relatively rapid diversification early in the evolution of the group, a rate 3-4 times faster than subsequent rates, and faster than comparable wet forest skinks.

A multilocus perspective on refugial isolation and divergence in rainforest skinks (Carlia).

To explore the evolutionary consequences of climate-induced fluctuations in distribution of rainforest habitat we contrasted demographic histories of divergence among three lineages of Australian rainforest endemic skinks. The red-throated rainbow skink, Carlia rubrigularis, consists of morphologically indistinguishable northern and southern mitochondrial DNA (mtDNA) lineages that are partially reproductively isolated at their parapatric boundary. The third lineage (C. rhomboidalis) inhabits rainforests just to the south of C. rubrigularis, has blue, rather than red-throated males, and for mtDNA is more closely related to southern C. rubrigularis than is northern C. rubrigularis. Multigene coalescent analyses supported more recent divergence between morphologically distinct lineages than between morphologically conservative lineages. There was effectively no migration and therefore stronger isolation between southern C. rubrigularis and C. rhomboidalis, and low unidirectional migration between morphologically conservative lineages of C. rubrigularis. We found little or no evidence for strong differences in effective population size, and hence different contributions of genetic drift in the demographic history of the three lineages. Overall the results suggest contrasting responses to long-term fluctuations in rainforest habitats, leading to varying opportunities for speciation.

A species assemblage approach to comparative phylogeography of birds in southern Australia.

We present a novel approach to investigating the divergence history of biomes and their component species using single-locus data prior to investing in multilocus data. We use coalescent-based hierarchical approximate Bayesian computation (HABC) methods (MsBayes) to estimate the number and timing of discrete divergences across a putative barrier and to assign species to their appropriate period of co-divergence. We then apply a coalescent-based full Bayesian model of divergence (IMa) to suites of species shown to have simultaneously diverged. The full Bayesian model results in reduced credibility intervals around divergence times and allows other parameters associated with divergence to be summarized across species assemblages. We apply this approach to 10 bird species that are wholly or patchily discontinuous in semi-arid habitats between Australia's southwest (SW) and southeast (SE) mesic zones. There was substantial support for up to three discrete periods of divergence. HABC indicates that two species wholly restricted to more mesic habitats diverged earliest, between 594,382 and 3,417,699 years ago, three species from semi-arid habitats diverged between 0 and 1,508,049 years ago, and four diverged more recently, between 0 and 396,843 years ago. Eight species were assigned to three periods of co-divergence with confidence. For full Bayesian analyses, we accounted for uncertainty in the two remaining species by analyzing all possible suites of species. Estimates of divergence times from full Bayesian divergence models ranged between 429,105 and 2,006,355; 67,172 and 663,837; and 24,607 and 171,085 for the earliest, middle, and most recent periods of co-divergence, respectively. This single-locus approach uses the power of multitaxa coalescent analyses as an efficient means of generating a foundation for further, targeted research using multilocus and genomic tools applied to an understudied biome.

Evidence for differential assortative female preference in association with refugial isolation of rainbow skinks in Australia's tropical rainforests.

BACKGROUND: Divergence driven by female preference can give rise to pre-mating isolation more rapidly than post-mating isolation can evolve through the accumulation of allelic incompatibilities. Moreover pre-mating isolation may be more effective at maintaining morphological differentiation between divergent populations. In the context of Australian rainforest endemic skinks that were historically subjected to refugial isolation, this study examined the following predictions: 1) that assortative female preference is associated with more recent divergence of southern C. rubrigularis (S-RED) and C. rhomboidalis (BLUE), but not with deeply divergent S-RED and northern C. rubrigularis (N-RED); and 2) that upon secondary contact, morphological differentiation is maintained between S-RED and BLUE, whereas N-RED and S-RED remain morphogically indistinguishable. PRINCIPAL FINDINGS: Female preference trials found no evidence for assortative female preference between N-RED and S-RED, supporting a previous genetic hybrid zone study which inferred post-mating but no pre-mating isolation. In contrast there is evidence for assortative female preference between S-RED and BLUE, with BLUE females preferring to associate with BLUE males, but S-RED females showing no preference. Multi-locus coalescent analyses, used to estimate post-divergence gene-flow between proximally located S-RED and BLUE populations, rejected zero gene-flow from BLUE to S-RED and thus RED and BLUE have maintained morphological differentiation despite secondary contact. Morphometric analyses confirmed a lack of morphological divergence between N-RED and S-RED and established that BLUE is morphologically divergent from RED in traits other than throat colour. CONCLUSIONS/SIGNIFICANCE: Long-term isolation has been sufficient to generate post-mating isolation but no morphological divergence between N-RED and S-RED. In contrast, greater morphological differentiation is associated with evidence for assortative female preference between more recently diverged S-RED and BLUE. Combined with previous estimates of lineage-wide gene flow, these results are consistent with the suggestion that assortative female preference is more effective than post-mating isolation in maintaining morphological differentiation between divergent populations.

Ancient DNA, pig domestication, and the spread of the Neolithic into Europe.

The Neolithic Revolution began 11,000 years ago in the Near East and preceded a westward migration into Europe of distinctive cultural groups and their agricultural economies, including domesticated animals and plants. Despite decades of research, no consensus has emerged about the extent of admixture between the indigenous and exotic populations or the degree to which the appearance of specific components of the "Neolithic cultural package" in Europe reflects truly independent development. Here, through the use of mitochondrial DNA from 323 modern and 221 ancient pig specimens sampled across western Eurasia, we demonstrate that domestic pigs of Near Eastern ancestry were definitely introduced into Europe during the Neolithic (potentially along two separate routes), reaching the Paris Basin by at least the early 4th millennium B.C. Local European wild boar were also domesticated by this time, possibly as a direct consequence of the introduction of Near Eastern domestic pigs. Once domesticated, European pigs rapidly replaced the introduced domestic pigs of Near Eastern origin throughout Europe. Domestic pigs formed a key component of the Neolithic Revolution, and this detailed genetic record of their origins reveals a complex set of interactions and processes during the spread of early farmers into Europe.

Phylogeny and ancient DNA of Sus provides insights into neolithic expansion in Island Southeast Asia and Oceania.

Human settlement of Oceania marked the culmination of a global colonization process that began when humans first left Africa at least 90,000 years ago. The precise origins and dispersal routes of the Austronesian peoples and the associated Lapita culture remain contentious, and numerous disparate models of dispersal (based primarily on linguistic, genetic, and archeological data) have been proposed. Here, through the use of mtDNA from 781 modern and ancient Sus specimens, we provide evidence for an early human-mediated translocation of the Sulawesi warty pig (Sus celebensis) to Flores and Timor and two later separate human-mediated dispersals of domestic pig (Sus scrofa) through Island Southeast Asia into Oceania. Of the later dispersal routes, one is unequivocally associated with the Neolithic (Lapita) and later Polynesian migrations and links modern and archeological Javan, Sumatran, Wallacean, and Oceanic pigs with mainland Southeast Asian S. scrofa. Archeological and genetic evidence shows these pigs were certainly introduced to islands east of the Wallace Line, including New Guinea, and that so-called "wild" pigs within this region are most likely feral descendants of domestic pigs introduced by early agriculturalists. The other later pig dispersal links mainland East Asian pigs to western Micronesia, Taiwan, and the Philippines. These results provide important data with which to test current models for human dispersal in the region.

The alpha-globin gene family of an Australian marsupial, Macropus eugenii: the long evolutionary history of the theta-globin gene and its functional status in mammals.

Comparative evolutionary analyses of gene families among divergent lineages can provide information on the order and timing of major gene duplication events and evolution of gene function. Here we investigate the evolutionary history of the alpha-globin gene family in mammals by isolating and characterizing alpha-like globin genes from an Australian marsupial, the tammar wallaby, Macropus eugenii. Sequence and phylogenetic analyses indicate that the tammar alpha-globin family consists of at least four genes including a single adult-expressed gene (alpha), two embryonic/neonatally expressed genes (zeta and zeta'), and theta-globin, each orthologous to the respective alpha-, zeta-, and theta-globin genes of eutherian mammals. The results suggest that the theta-globin lineage arose by duplication of an ancestral adult alpha-globin gene and had already evolved an unusual promoter region, atypical of all known alpha-globin gene promoters, prior to the divergence of the marsupial and eutherian lineages. Evolutionary analyses, using a maximum likelihood approach, indicate that theta-globin, has evolved under strong selective constraints in both marsupials and the lineage leading to human theta-globin, suggesting a long-term functional status. Overall, our results indicate that at least a four-gene cluster consisting of three alpha-like and one beta-like globin genes linked in the order 5'-zeta-alpha-theta-omega-3' existed in the common ancestor of marsupials and eutherians. However, results are inconclusive as to whether the two tammar zeta-globin genes arose by duplication prior to the radiation of the marsupial and eutherian lineages, with maintenance of exon sequences by gene conversion, or more recently within marsupials.