Fine-scale biogeographical and temporal diversification processes of peacock swallowtails (PapiliosubgenusAchillides) in the Indo-Australian Archipelago

A comprehensive phylogeny and revised taxonomy illuminate the origin and diversification of the global radiation of Papilio (Lepidoptera: Papilionidae)

The swallowtail genus Papilio (Lepidoptera: Papilionidae) is species rich, distributed worldwide, and has broad morphological habits and ecological niches. Because of its elevated species richness, it has been historically difficult to reconstruct a densely sampled phylogeny for this clade. Here we provide a taxonomic working list for the genus, resulting in 235 Papilio species, and assemble a molecular dataset of seven gene fragments representing ca. 80% of the currently described diversity. Phylogenetic analyses reconstructed a robust tree with highly supported relationships within subgenera, although a few nodes in the early history of the Old World Papilio remain unresolved. Contrasting with previous results, we found that Papilio alexanor is sister to all Old World Papilio and that the subgenus Eleppone is no longer monotypic. The latter includes the recently described Fijian Papilio natewa with the Australian Papilio anactus and is sister to subgenus Araminta (formerly included in subgenus Menelaides) occurring in Southeast Asia. Our phylogeny also includes rarely studied (P. antimachus, P. benguetana) or endangered species (P. buddha, P. chikae). Taxonomic changes resulting from this study are elucidated. Molecular dating and biogeographic analyses indicate that Papilio originated ca. 30 million years ago (Oligocene), in a northern region centered on Beringia. A rapid early Miocene radiation in the Paleotropics is revealed within Old World Papilio, potentially explaining their low early branch support. Most subgenera originated in the early to middle Miocene followed by synchronous southward biogeographic dispersals and repeated local extirpations in northern latitudes. This study provides a comprehensive phylogenetic framework for Papilio with clarification of subgeneric systematics and species taxonomic changes enumerated, which will facilitate further studies to address questions on their ecology and evolutionary biology using this model clade.

Molecular phylogeny and biogeography of the genus Symbrenthia (Lepidoptera, Nymphalidae) correlates with the past geography of the Oriental region

Jesters, butterflies in the genus Symbrenthia Hübner, 1819, comprise 14 species mainly distributed in the Oriental region. Although this genus has attracted the attention of many researchers in the past, its taxonomy and biogeographic history remain unclear. In this study, we investigate phylogenetic and biogeographic relationships inferred from one mitochondrial (COI) and two nuclear genes (ArgKin, wingless), using both likelihood and Bayesian approaches. With the exception of S. hippalus, which we find to be either sister to Mynes Boisduval, 1832 or sister to Symbrenthia + Mynes + Araschnia, all species of Symbrenthia form a single monophyletic group. We describe a new genus Mynbrenthia Fric & Rindos gen. nov. to accommodate the taxon hippalus. The genus Symbrenthia splits into four sub-groups, "Brensymthia" (with S. niphanda and S. sinoides), "hypselis" (with S. hypselis, S. brabira, S. leoparda and S. doni), "intricata" (with S. intricata and S. hypatia) and "hippoclus" group (including S. platena and a complex of S. hippoclus and S. lilaea). The genus probably originated in Sundaland or continental Asia during the Eocene. The history of the genus Symbrenthia was more influenced by dispersal events and then by subsequent vicariances. Whereas the "hypselis" group colonised the Indo-Australian Archipelago from the Asian continent, the "hippoclus" group dispersed to continental Asia from the Indo-Australian Archipelago.

Admixture of evolutionary rates across a butterfly hybrid zone

Hybridization is a major evolutionary force that can erode genetic differentiation between species, whereas reproductive isolation maintains such differentiation. In studying a hybrid zone between the swallowtail butterflies Papilio syfanius and Papilio maackii (Lepidoptera: Papilionidae), we made the unexpected discovery that genomic substitution rates are unequal between the parental species. This phenomenon creates a novel process in hybridization, where genomic regions most affected by gene flow evolve at similar rates between species, while genomic regions with strong reproductive isolation evolve at species-specific rates. Thus, hybridization mixes evolutionary rates in a way similar to its effect on genetic ancestry. Using coalescent theory, we show that the rate-mixing process provides distinct information about levels of gene flow across different parts of genomes, and the degree of rate-mixing can be predicted quantitatively from relative sequence divergence ([Formula: see text]) between the hybridizing species at equilibrium. Overall, we demonstrate that reproductive isolation maintains not only genomic differentiation, but also the rate at which differentiation accumulates. Thus, asymmetric rates of evolution provide an additional signature of loci involved in reproductive isolation.

Are the Yellow and Red Marked Club-Tail Losaria coon the Same Species?

Losaria coon (Fabricius, 1793) is currently comprised of ten subspecies, which were originally described under two names, Papilio coon and P. doubledayi before 1909, when they were combined as one species. The main difference between them is the colour of abdomen and hindwing subterminal spots-yellow in coon and red in doubledayi. Wing morphology, male and female genitalia, and molecular evidence (DNA barcodes) were analysed for multiple subspecies of L. coon and three other Losaria species-rhodifer, neptunus, and palu. Our molecular data support the separation of L. coon and L. doubledayi stat. rev. as two distinct species, with L. rhodifer positioned between them in phylogenetic analyses. Wing morphology and genitalic structures also confirm the molecular conclusions. Our findings divide L. coon into two species occupying different geographic ranges: with L. coon restricted to southern Sumatra, Java, and Bawean Island, while L. doubledayi occurs widely in regions from North India to northern Sumatra, including Hainan and Nicobar Islands. Hence, future conservation efforts must reassess the status and threat factors of the two species to form updated strategies.

A glide over the Indo-Australian geological maze: repeated transgressions of Lydekker’s and Wallace’s Lines in archdukes, barons and dukes (Nymphalidae: Limenitidinae: Adoliadini)

Archdukes, barons, counts, dukes and marquises are forest-dwelling butterflies found in mainland Asia and most islands of the Indo-Australian archipelago west of Wallace’s Line, with only a few species occurring as far east as the Bismarck Archipelago. This pattern is unusual among butterfly groups of the region, which often present more widespread geographical ranges bearing little signature of Lydekker’s and Wallace’s Lines. Using a molecular multimarker matrix, we infer the first dated phylogeny for this clade and estimate its biogeographical history. We recover the Oriental genus Euthalia as polyphyletic, although other genera are monophyletic. The clade originated in continental Indomalaya in the late Oligocene ~24 Mya, when the Indo-Australian archipelago was at one of the most dynamic stages of its orogeny. Multiple independent colonization events towards the Lesser Sunda Islands, Moluccas, Australia and New Guinea suggest the relative permeability of Lydekker’s and Wallace’s Lines to these butterflies. Colonization of Melanesia took place twice, probably before the recent formation of Sulawesi. The study of Indo-Australian Adoliadini provides additional evidence that biogeographical barriers long thought to prevent exchange between the Asian and Australian biotas are, in fact, permeable especially to vagile insect lineages in the region.

Molecular phylogeny and historical biogeography of Parnara butterflies (Lepidoptera: Hesperiidae)

The butterfly genus Parnara (Hesperiinae: Baorini), of which some are major pests of economic crops (e.g., rice, wild rice stems and sugarcane), currently consists of 10 species and several subspecies and has a highly disjunct distribution in Australia, Africa, and Asia. We determined the systematic relationships and biogeographical history of the genus by reconstructing the phylogeny based on eight genes and 101 specimens representing all 10 recognized species. Four species delimitation methods (ABGD, bPTP, GMYC and BPP) were also employed to assess the taxonomic status of each species. Based on these results and analyses, we recognize 11 extant species in the genus. The status of the taxon P. naso poutieri (Boisduval, 1833) from Madagascar is revised as a distinct species, Parnara poutieri (Boisduval, 1833) stat. rev. The subspecies P. guttata mangala (Moore, 1866) syn. nov. is synonymized with P. guttata guttata (Bremer & Grey, 1853), while P. bada (Moore, 1878) is provisionally treated as a complex of two species, namely P. bada and P. apostata (Snellen, 1886). The monophyly of Parnara is strongly supported, with the following relationships: P. amalia + ((P. monasi + (P. poutieri + P. naso)) + ((P. kawazoei + P. bada complex) + (P. ganga + (P. ogasawarensis + (P. guttata + P. batta))))). Divergence time and ancestral range estimates indicate that the common ancestor of Parnara originated in an implausible area of Australia, Africa, and Oriental region in the mid-Oligocene and then differentiated in the late Miocene-late Pliocene. Dispersal and range expansion have played an important role in diversification of the genus in Asia and Afica. Relatively stable geotectonic plates at the time when most extant lineages appeared during the late Miocene-early Pliocene might have been the factor responsible for the relatively constant low dynamic rate of diversification within the group.

Taxon cycle predictions supported by model-based inference in Indo-Pacific trap-jaw ants (Hymenoptera: Formicidae: Odontomachus)

Nonequilibrium dynamics and non-neutral processes, such as trait-dependent dispersal, are often missing from quantitative island biogeography models despite their potential explanatory value. One of the most influential nonequilibrium models is the taxon cycle, but it has been difficult to test its validity as a general biogeographical framework. Here, we test predictions of the taxon cycle model using six expected phylogenetic patterns and a time-calibrated phylogeny of Indo-Pacific Odontomachus (Hymenoptera: Formicidae: Ponerinae), one of the ant genera that E.O. Wilson used when first proposing the hypothesis. We used model-based inference and a newly developed trait-dependent dispersal model to jointly estimate ancestral biogeography, ecology (habitat preferences for forest interiors, vs. "marginal" habitats, such as savannahs, shorelines, disturbed areas) and the linkage between ecology and dispersal rates. We found strong evidence that habitat shifts from forest interior to open and disturbed habitats increased macroevolutionary dispersal rate. In addition, lineages occupying open and disturbed habitats can give rise to both island endemics re-occupying only forest interiors and taxa that re-expand geographical ranges. The phylogenetic predictions outlined in this study can be used in future work to evaluate the relative weights of neutral (e.g., geographical distance and area) and non-neutral (e.g., trait-dependent dispersal) processes in historical biogeography and community ecology.

The analysis of genetic variation in the mitochondrial genome and its application for the identification of Papilio species

Mitochondrial DNA (mtDNA) markers are ideal for evolutionary studies, including phylogeography, population genetics, phylogeny, etc. However, different mitochondrial genes always own different evolutionary rate. In this study, we analysed the genetic variation across the 16 complete mtDNA from 13 species in the genus Papilio and recognized the best DNA barcoding for Papilio species. The mitochondrial gene arrangement for each species shares a similar order, similar to the typical Papilionidae species, which indicated the relatively conservative state of gene arrangement in Papilio. The sliding window of genetic diversity showed that there was a significant difference in the genetic diversity of each gene in the mitochondrial genome of Papilio. The relatively mean clock rate of the ND1 was broadly lower than the other genes in mitochondrial genome of Papilio; while the ATP8 owns the largest values of mean clock rate. Those results suggested that the rate of evolution of each gene is not balanced and all mitochondrial genes except ND1 and ATP8 could act as barcoding for the identification of Papilio species. The phylogenetic analyses of complete mtDNA data for 13 Papilio species divided them into five major branches, which keep the same topological structure with previous studies.

An ant genus-group (Prenolepis) illuminates the biogeography and drivers of insect diversification in the Indo-Pacific

The Malay Archipelago and the tropical South Pacific (hereafter the Indo-Pacific region) are considered biodiversity hotspots, yet a general understanding of the origins and diversification of species-rich groups in the region remains elusive. We aimed to test hypotheses for the evolutionary processes driving insect species diversity in the Indo-Pacific using a higher-level and comprehensive phylogenetic hypothesis for an ant clade consisting of seven genera. We estimated divergence times and reconstructed the biogeographical history of ant species in the Prenolepis genus-group (Formicidae: Formicinae: Lasiini). We used a fossil-calibrated phylogeny to infer ancestral geographical ranges utilizing a biogeographic model that includes founder-event speciation. Ancestral state reconstructions of the ants' ecological preferences, and diversification rates were estimated for selected Indo-Pacific clades. Overall, we report that faunal interchange between Asia and Australia has occurred since at least 20-25 Ma, and early dispersal to the Fijian Basin happened during the early and mid-Miocene (ca. 10-20 Ma). Differences in diversification rates across Indo-Pacific clades may be related to ecological preference breadth, which in turn may have facilitated geographical range expansions. Ancient dispersal routes suggested by our results agree with the palaeogeography of the region. For this particular group of ants, the rapid orogenesis in New Guinea and possibly subsequent ecological shifts may have promoted their rapid diversification and widespread distribution across the Indo-Pacific.

Islands conserve high species richness and areas of endemism of Hormaphidinae aphids

Patterns of biodiversity and endemism are important and form the foundation for biogeography and conservation studies. Hormaphidinae is an aphid group mainly distributed in the Sino-Japanese, Oriental, and Nearctic zoogeographic realms. To infer the areas of endemism of Hormaphidinae aphids in the Sino-Japanese and Oriental realms, we employed a geographical distribution dataset covering all 225 species in subfamily Hormaphidinae. In total, 1,245 distribution occurrence records for all species were analyzed in addition to the number of species in a certain grid cell to calculate species richness. Two methods (parsimony analysis of endemicity (PAE) and the use of an optimization criterion-NDM/VNDM) using a total of 6 grid sizes were applied to detect possible areas of endemism and to assess the probable effects of the 2 different methods and 6 grid sizes on the results of the patterns of the areas of endemism. The results revealed that islands presented most of the areas with high species richness and endemic species, particularly Japan, Taiwan Island, Java Island, the Malaysian Peninsula, southeast Himalaya, and the Hengduan Mountains. Most of these areas of endemism were located on islands or a peninsula. Islands were therefore shown to play an essential role in facilitating the formation of high species richness and endemism. Different grid sizes directly determined the regions of the areas of endemism, with small grid sizes tending to detect small and discrete areas of endemism and large grid sizes identifying continuous areas of endemism. In terms of the two methods, NDM/VNDM was found to identify more areas of endemism than PAE.

EF-1α DNA Sequences Indicate Multiple Origins of Introduced Populations of Essigella californica (Hemiptera: Aphididae)

Aphids in the pine-feeding Nearctic genus Essigella (Sternorrhyncha, Aphididae, Lachninae) have been introduced in Europe, North Africa, Oceania, and South America. Mitochondrial, nuclear, and endosymbiont DNA sequences of 12 introduced populations from three continents confirm they all belong to Essigella californica (Essig, 1909). Intron sequence variation of the nuclear gene EF-1α has revealed the existence of four distinct groups. Group I gathers one population from China, where the species is newly reported, and several from Europe (France and Italy); Group II is represented by one population from Argentina; Group III includes two populations from Southern Australia with one from New Zealand; and Group IV corresponds to five populations from Eastern and South-Eastern Australia. These results indicate that introduced populations of E. californica have at least four source populations. They also show that intron variation of EF-1α can be a method to discriminate populations of asexually reproducing aphids.

Advances using molecular data in insect systematics

The size of molecular datasets has been growing exponentially since the mid 1980s, and new technologies have now dramatically increased the slope of this increase. New datasets include genomes, transcriptomes, and hybrid capture data, producing hundreds or thousands of loci. With these datasets, we are approaching a consensus on the higher level insect phylogeny. Huge datasets can produce new challenges in interpreting branch support, and new opportunities in developing better models and more sophisticated partitioning schemes. Dating analyses are improving as we recognize the importance of careful fossil calibration selection. With thousands of genes now available, coalescent methods have come of age. Barcode libraries continue to expand, and new methods are being developed for incorporating them into phylogenies with tens of thousands of individuals.

Macroevolution of hyperdiverse flightless beetles reflects the complex geological history of the Sunda Arc

The Sunda Arc forms an almost continuous chain of islands and thus a potential dispersal corridor between mainland Southeast Asia and Melanesia. However, the Sunda Islands have rather different geological histories, which might have had an important impact on actual dispersal routes and community assembly. Here, we reveal the biogeographical history of hyperdiverse and flightless Trigonopterus weevils. Different approaches to ancestral area reconstruction suggest a complex east to west range expansion. Out of New Guinea, Trigonopterus repeatedly reached the Moluccas and Sulawesi transgressing Lydekker's Line. Sulawesi repeatedly acted as colonization hub for different segments of the Sunda Arc. West Java, East Java and Bali are recognized as distinct biogeographic areas. The timing and diversification of species largely coincides with the geological chronology of island emergence. Colonization was not inhibited by traditional biogeographical boundaries such as Wallace's Line. Rather, colonization patterns support distance dependent dispersal and island age limiting dispersal.

Diversification of clearwing butterflies with the rise of the Andes

AIM

Despite the greatest butterfly diversity on Earth occurring in the Neotropical Andes and Amazonia, there is still keen debate about the origins of this exceptional biota. A densely sampled calibrated phylogeny for a widespread butterfly subtribe, Oleriina (Nymphalidae: Ithomiini) was used to estimate the origin, colonization history and diversification of this species-rich group.

LOCATION

Neotropics.

METHODS

Ancestral elevation and biogeographical ranges were reconstructed using data generated from detailed range maps and applying the dispersal-extinction-cladogenesis model using stratified palaeogeographical time slice matrices. The pattern of diversification through time was examined by comparing constant and variable rate models. We also tested the hypothesis that a change in elevation is associated with speciation.

RESULTS

The Oleriina likely originated in the Andes in the Early to Middle Miocene and rapidly diversified to include four genera all of which also originated in the Andes. These clades, together with four species groups, experienced varying spatial and temporal patterns of diversification. An overall early burst and decreasing diversification rate is identified, and this pattern is reflected for most subclades.

MAIN CONCLUSIONS

Changes in the palaeogeological landscape, particularly the prolonged uplift of the Andes, had a profound impact on the diversification of the subtribe. The Oleriina mostly remained within the Andes and vicariant speciation resulted in some instances. Dynamic dispersal occurred with the disappearance of geological barriers such as the Acre System and the subtribe exploited newly available habitats. Our results confirm the role of the Andean uplift in the evolution of Neotropical biodiversity.

Phylogeny and Historical Biogeography of Asian Pterourus Butterflies (Lepidoptera: Papilionidae): A Case of Intercontinental Dispersal from North America to East Asia

The phylogenetic status of the well-known Asian butterflies often known as Agehana (a species group, often treated as a genus or a subgenus, within Papilio sensu lato) has long remained unresolved. Only two species are included, and one of them especially, Papilio maraho, is not only rare but near-threatened, being monophagous on its vulnerable hostplant, Sassafras randaiense (Lauraceae). Although the natural history and population conservation of "Agehana" has received much attention, the biogeographic origin of this group still remains enigmatic. To clarify these two questions, a total of 86 species representatives within Papilionidae were sampled, and four genes (concatenated length 3842 bp) were used to reconstruct their phylogenetic relationships and historical scenarios. Surprisingly, "Agehana" fell within the American Papilio subgenus Pterourus and not as previously suggested, phylogenetically close to the Asian Papilio subgenus Chilasa. We therefore formally synonymize Agehana with Pterourus. Dating and biogeographic analysis allow us to infer an intercontinental dispersal of an American ancestor of Asian Pterourus in the early Miocene, which was coincident with historical paleo-land bridge connections, resulting in the present "East Asia-America" disjunction distribution. We emphasize that species exchange between East Asia and America seems to be a quite frequent occurrence in butterflies during the Oligocene to Miocene climatic optima.

Deciphering the evolution of birdwing butterflies 150 years after Alfred Russel Wallace

One hundred and fifty years after Alfred Wallace studied the geographical variation and species diversity of butterflies in the Indomalayan-Australasian Archipelago, the processes responsible for their biogeographical pattern remain equivocal. We analysed the macroevolutionary mechanisms accounting for the temporal and geographical diversification of the charismatic birdwing butterflies (Papilionidae), a major focus of Wallace's pioneering work. Bayesian phylogenetics and dating analyses of the birdwings were conducted using mitochondrial and nuclear genes. The combination of maximum likelihood analyses to estimate biogeographical history and diversification rates reveals that diversity-dependence processes drove the radiation of birdwings, and that speciation was often associated with founder-events colonizing new islands, especially in Wallacea. Palaeo-environment diversification models also suggest that high extinction rates occurred during periods of elevated sea level and global warming. We demonstrated a pattern of spatio-temporal habitat dynamics that continuously created or erased habitats suitable for birdwing biodiversity. Since birdwings were extinction-prone during the Miocene (warmer temperatures and elevated sea levels), the cooling period after the mid-Miocene climatic optimum fostered birdwing diversification due to the release of extinction. This also suggests that current global changes may represent a serious conservation threat to this flagship group.

Biogeography of Anurans from the Poorly Known and Threatened Coastal Sandplains of Eastern Brazil

The east coast of Brazil comprises an extensive area inserted in the Tropical Atlantic Domain and is represented by sandy plains of beach ridges commonly known as Restingas. The coastal environments are unique and house a rich amphibian fauna, the geographical distribution patterns of which are incipient. Biogeographical studies can explain the current distributional patterns and provide the identification of natural biogeographical units. These areas are important in elucidating the evolutionary history of the taxa and the areas where they occur. The aim of this study was to seek natural biogeographical units in the Brazilian sandy plains of beach ridges by means of distribution data of amphibians and to test the main predictions of the vicariance model to explain the patterns found. We revised and georeferenced data on the geographical distribution of 63 anuran species. We performed a search for latitudinal distribution patterns along the sandy coastal plains of Brazil using the non-metric multidimensional scaling method (NMDS) and the biotic element analysis to identify natural biogeographical units. The results showed a monotonic variation in anuran species composition along the latitudinal gradient with a break in the clinal pattern from 23°S to 25°S latitude (states of Rio de Janeiro to São Paulo). The major predictions of the vicariance model were corroborated by the detection of four biotic elements with significantly clustered distribution and by the presence of congeneric species distributed in distinct biotic elements. The results support the hypothesis that vicariance could be one of the factors responsible for the distribution patterns of the anuran communities along the sandy coastal plains of eastern Brazil. The results of the clusters are also congruent with the predictions of paleoclimatic models made for the Last Glacial Maximum of the Pleistocene, such as the presence of historical forest refugia and biogeographical patterns already detected for amphibians in the Atlantic Rainforest.

The Scirtothrips dorsalis Species Complex: Endemism and Invasion in a Global Pest

Invasive arthropods pose unique management challenges in various environments, the first of which is correct identification. This apparently mundane task is particularly difficult if multiple species are morphologically indistinguishable but accurate identification can be determined with DNA barcoding provided an adequate reference set is available. Scirtothrips dorsalis is a highly polyphagous plant pest with a rapidly expanding global distribution and this species, as currently recognized, may be comprised of cryptic species. Here we report the development of a comprehensive DNA barcode library for S. dorsalis and seven nuclear markers via next-generation sequencing for identification use within the complex. We also report the delimitation of nine cryptic species and two morphologically distinguishable species comprising the S. dorsalis species complex using histogram analysis of DNA barcodes, Bayesian phylogenetics, and the multi-species coalescent. One member of the complex, here designated the South Asia 1 cryptic species, is highly invasive, polyphagous, and likely the species implicated in tospovirus transmission. Two other species, South Asia 2, and East Asia 1 are also highly polyphagous and appear to be at an earlier stage of global invasion. The remaining members of the complex are regionally endemic, varying in their pest status and degree of polyphagy. In addition to patterns of invasion and endemism, our results provide a framework both for identifying members of the complex based on their DNA barcode, and for future species delimiting efforts.

Critiquing blind dating: the dangers of over-confident date estimates in comparative genomics

Phylogenomic advances provide more rigorous estimates for the timing of evolutionary divergences than previously available (e.g., Bayesian relaxed-clock estimates with soft fossil constraints). However, because many family-level clades and higher, as well as model species within those clades, have not been included in phylogenomic studies, the literature presents temporal estimates likely harboring substantial errors. Blindly using such dates can substantially retard scientific advancement. We suggest a way forward by conducting analyses that minimize prior assumptions and use large datasets, and demonstrate how using such a phylogenomic approach can lead to significantly more parsimonious conclusions without a good fossil record. We suggest that such an approach calls for research into the biological causes of conflict between molecular and fossil signatures.

Australasian sky islands act as a diversity pump facilitating peripheral speciation and complex reversal from narrow endemic to widespread ecological supertramp

The Australasian archipelago is biologically extremely diverse as a result of a highly puzzling geological and biological evolution. Unveiling the underlying mechanisms has never been more attainable as molecular phylogenetic and geological methods improve, and has become a research priority considering increasing human-mediated loss of biodiversity. However, studies of finer scaled evolutionary patterns remain rare particularly for megadiverse Melanesian biota. While oceanic islands have received some attention in the region, likewise insular mountain blocks that serve as species pumps remain understudied, even though Australasia, for example, features some of the most spectacular tropical alpine habitats in the World. Here, we sequenced almost 2 kb of mitochondrial DNA from the widespread diving beetle Rhantus suturalis from across Australasia and the Indomalayan Archipelago, including remote New Guinean highlands. Based on expert taxonomy with a multigene phylogenetic backbone study, and combining molecular phylogenetics, phylogeography, divergence time estimation, and historical demography, we recover comparably low geographic signal, but complex phylogenetic relationships and population structure within R. suturalis. Four narrowly endemic New Guinea highland species are subordinated and two populations (New Guinea, New Zealand) seem to constitute cases of ongoing speciation. We reveal repeated colonization of remote mountain chains where haplotypes out of a core clade of very widespread haplotypes syntopically might occur with well-isolated ones. These results are corroborated by a Pleistocene origin approximately 2.4 Ma ago, followed by a sudden demographic expansion 600,000 years ago that may have been initiated through climatic adaptations. This study is a snapshot of the early stages of lineage diversification by peripatric speciation in Australasia, and supports New Guinea sky islands as cradles of evolution, in line with geological evidence suggesting very recent origin of high altitudes in the region.