Phylogeography of two New Zealand lizards: McCann’s skink (Oligosoma maccanni) and the brown skink (O. zelandicum)

Oligosoma aureocola sp. nov. (Reptilia: Scincidae) from the northern Southland high country of Aotearoa/New Zealand

A species of diurnal skink from the Mataura Range and Mid Dome in central northern Southland, Aotearoa/New Zealand is described as Oligosoma aureocola sp. nov. It is a small species, coloured mid- to dark brown with smooth longitudinal stripes, and lives along rocky alpine ridges, low-stature shrublands, and tussock grasslands. This skink is a conspicuous species, easily sighted basking and foraging in talus or alpine plants such as golden spear grass (Aciphylla aurea).

Mammal-related Cryptosporidium infections in endemic reptiles of New Zealand

New Zealand's endemic reptile fauna is highly threatened and pathogens causing infectious diseases may be a significant risk to already endangered species. Here, we investigate Cryptosporidium infection in captive endemic New Zealand reptiles. We found two mammal-related Cryptosporidium species (C. hominis and C. parvum) and six subtypes from three gp60 families (Ib, Ig and IIa) in 12 individuals of captive endemic Tuatara, Otago and Grand skinks, and Jewelled and Rough geckos. Cryptosporidium serpentis was identified in two Jewelled geckos using 18S. In New Zealand, C. hominis and C. parvum are associated with infections in humans and introduced domestic animals but have also been recently found in wildlife. Our finding of Cryptosporidium infection in endemic reptiles can help inform strategies to monitor the conservation of species and manage potential introductions of pathogens to in-situ and ex-situ populations.

Mitogenomes resolve the phylogeography and divergence times within the endemic New Zealand Callaeidae (Aves: Passerida)

The biogeographical origins of the endemic birds of New Zealand (Aotearoa) are of great interest, particularly Palaeogene lineages such as Callaeidae, a passerine family characterized by brightly coloured wattles behind the beak and, in some cases, extreme sexual dimorphism in bill size and shape. Ancestral representatives of Callaeidae are thought to have split from their closest relatives outside New Zealand in the Oligocene, but little is known about the timing of divergences within the family. We present a fully dated molecular phylogeny of Callaeidae mitogenomes and discuss the biogeographical implications. Our results suggest that formation of Pliocene marine seaways, such as the Manawatu Strait, are likely to have played a significant role in the differentiation of North Island and South Island kōkako (Callaeas spp.) and saddlebacks/tīeke (Philesturnus spp.).

Spatial variation in the climatic predictors of species compositional turnover and endemism

Previous research focusing on broad-scale or geographically invariant species-environment dependencies suggest that temperature-related variables explain more of the variation in reptile distributions than precipitation. However, species-environment relationships may exhibit considerable spatial variation contingent upon the geographic nuances that vary between locations. Broad-scale, geographically invariant analyses may mask this local variation and their findings may not generalize to different locations at local scales. We assess how reptile-climatic relationships change with varying spatial scale, location, and direction. Since the spatial distributions of diversity and endemism hotspots differ for other species groups, we also assess whether reptile species turnover and endemism hotspots are influenced differently by climatic predictors. Using New Zealand reptiles as an example, the variation in species turnover, endemism and turnover in climatic variables was measured using directional moving window analyses, rotated through 360°. Correlations between the species turnover, endemism and climatic turnover results generated by each rotation of the moving window were analysed using multivariate generalized linear models applied at national, regional, and local scales. At national-scale, temperature turnover consistently exhibited the greatest influence on species turnover and endemism, but model predictive capacity was low (typically r (2) = 0.05, P < 0.001). At regional scales the relative influence of temperature and precipitation turnover varied between regions, although model predictive capacity was also generally low. Climatic turnover was considerably more predictive of species turnover and endemism at local scales (e.g., r (2) = 0.65, P < 0.001). While temperature turnover had the greatest effect in one locale (the northern North Island), there was substantial variation in the relative influence of temperature and precipitation predictors in the remaining four locales. Species turnover and endemism hotspots often occurred in different locations. Climatic predictors had a smaller influence on endemism. Our results caution against assuming that variability in temperature will always be most predictive of reptile biodiversity across different spatial scales, locations and directions. The influence of climatic turnover on the species turnover and endemism of other taxa may exhibit similar patterns of spatial variation. Such intricate variation might be discerned more readily if studies at broad scales are complemented by geographically variant, local-scale analyses.

A retrospective approach to testing the DNA barcoding method

A decade ago, DNA barcoding was proposed as a standardised method for identifying existing species and speeding the discovery of new species. Yet, despite its numerous successes across a range of taxa, its frequent failures have brought into question its accuracy as a short-cut taxonomic method. We use a retrospective approach, applying the method to the classification of New Zealand skinks as it stood in 1977 (primarily based upon morphological characters), and compare it to the current taxonomy reached using both morphological and molecular approaches. For the 1977 dataset, DNA barcoding had moderate-high success in identifying specimens (78-98%), and correctly flagging specimens that have since been confirmed as distinct taxa (77-100%). But most matching methods failed to detect the species complexes that were present in 1977. For the current dataset, there was moderate-high success in identifying specimens (53-99%). For both datasets, the capacity to discover new species was dependent on the methodological approach used. Species delimitation in New Zealand skinks was hindered by the absence of either a local or global barcoding gap, a result of recent speciation events and hybridisation. Whilst DNA barcoding is potentially useful for specimen identification and species discovery in New Zealand skinks, its error rate could hinder the progress of documenting biodiversity in this group. We suggest that integrated taxonomic approaches are more effective at discovering and describing biodiversity.

Limited, episodic diversification and contrasting phylogeography in a New Zealand cicada radiation

BACKGROUND

The New Zealand (NZ) cicada fauna contains two co-distributed lineages that independently colonized the isolated continental fragment in the Miocene. One extensively studied lineage includes 90% of the extant species (Kikihia + Maoricicada + Rhodopsalta; ca 51 spp.), while the other contains just four extant species (Amphipsalta - 3 spp. + Notopsalta - 1 sp.) and has been little studied. We examined mitochondrial and nuclear-gene phylogenies and phylogeography, Bayesian relaxed-clock divergence timing (incorporating literature-based uncertainty of molecular clock estimates) and ecological niche models of the species from the smaller radiation.

RESULTS

Mitochondrial and nuclear-gene trees supported the monophyly of Amphipsalta. Most interspecific diversification within Amphipsalta-Notopsalta occurred from the mid-Miocene to the Pliocene. However, interspecific divergence time estimates had large confidence intervals and were highly dependent on the assumed tree prior, and comparisons of uncorrected and patristic distances suggested difficulty in estimation of branch lengths. In contrast, intraspecific divergence times varied little across analyses, and all appear to have occurred during the Pleistocene. Two large-bodied forest taxa (A. cingulata, A. zelandica) showed minimal phylogeographic structure, with intraspecific diversification dating to ca. 0.16 and 0.37 Ma, respectively. Mid-Pleistocene-age phylogeographic structure was found within two smaller-bodied species (A. strepitans - 1.16 Ma, N. sericea - 1.36 Ma] inhabiting dry open habitats. Branches separating independently evolving species were long compared to intraspecific branches. Ecological niche models hindcast to the Last Glacial Maximum (LGM) matched expectations from the genetic datasets for A. zelandica and A. strepitans, suggesting that the range of A. zelandica was greatly reduced while A. strepitans refugia were more extensive. However, no LGM habitat could be reconstructed for A. cingulata and N. sericea, suggesting survival in microhabitats not detectable with our downscaled climate data.

CONCLUSIONS

Unlike the large and continuous diversification exhibited by the Kikihia-Maoricicada-Rhodopsalta clade, the contemporaneous Amphipsalta-Notopsalta lineage contains four comparatively old (early branching) species that show only recent diversification. This indicates either a long period of stasis with no speciation, or one or more bouts of extinction that have pruned the radiation. Within Amphipsalta-Notopsalta, greater population structure is found in dry-open-habitat species versus forest specialists. We attribute this difference to the fact that NZ lowland forests were repeatedly reduced in extent during glacial periods, while steep, open habitats likely became more available during late Pleistocene uplift.

Phylogeography of the endangered Otago skink, Oligosoma otagense: population structure, hybridisation and genetic diversity in captive populations

Climatic cooling and substantial tectonic activity since the late Miocene have had a pronounced influence on the evolutionary history of the fauna of New Zealand's South Island. However, many species have recently experienced dramatic range reductions due to habitat fragmentation and the introduction of mammalian predators and competitors. These anthropogenic impacts have been particularly severe in the tussock grasslands of the Otago region. The Otago skink (Oligosoma otagense), endemic to the region, is one of the most critically endangered vertebrates in New Zealand. We use mitochondrial DNA sequence data to investigate the evolutionary history of the Otago skink, examine its population genetic structure, and assess the level of genetic diversity in the individuals in the captive breeding program. Our data indicate that the Otago skink diverged from its closest relatives in the Miocene, consistent with the commencement of tectonic uplift of the Southern Alps. However, there is evidence for past introgression with the scree skink (O. waimatense) in the northern Otago-southern Canterbury region. The remnant populations in eastern Otago and western Otago are estimated to have diverged in the mid-Pliocene, with no haplotypes shared between these two regions. This divergence accounts for 95% of the genetic diversity in the species. Within both regions there is strong genetic structure among populations, although shared haplotypes are generally evident between adjacent localities. Although substantial genetic diversity is present in the captive population, all individuals originate from the eastern region and the majority had haplotypes that were not evident in the intensively managed populations at Macraes Flat. Our data indicate that eastern and western populations should continue to be regarded as separate management units. Knowledge of the genetic diversity of the breeding stock will act to inform the captive management of the Otago skink and contribute to a key recovery action for the species.

Onset of glaciation drove simultaneous vicariant isolation of Alpine insects in New Zealand

The origin of the New Zealand "beech gap," a low-diversity zone in the central South Island corresponding with a disjunction in the distribution of many taxa, has been the focus of biogeographical debate for many decades. Here, we use comparative phylogeographic analysis (COI; H3) of six alpine stonefly genera (116 individuals, 102 localities) to test a vicariant evolutionary hypothesis for the origin of this "biotic gap." We find strikingly similar phylogeographic patterns in all six genera, with the deepest genetic divergences always found between samples north and south of the beech gap. The magnitude of north-south genetic differentiation for COI is similar across all six genera (ranging from 0.074 to 0.091), with a test for simultaneous vicariance confirming that divergence is consistent with a single evolutionary event. The concordant cladogenesis detected across multiple taxa is consistent with vicariant isolation caused by the onset of glaciation in the late Pliocene. This study thus indicates an important cladogenetic role for glaciation, an abiotic evolutionary process that is more typically associated with loss of biodiversity.

A SINE of restricted gene flow across the Alpine Fault: phylogeography of the New Zealand common skink (Oligosoma nigriplantare polychroma)

New Zealand has experienced a complex climatic and geological history since the Pliocene. Thus, identifying the processes most important in having driven the evolution of New Zealand's biota has proven difficult. Here we examine the phylogeography of the New Zealand common skink (Oligosoma nigriplantare polychroma) which is distributed throughout much of New Zealand and crosses many putative biogeographical boundaries. Using mitochondrial DNA sequence data, we revealed five geographically distinct lineages that are highly differentiated (pairwise Phi(ST) 0.54-0.80). The phylogeographical pattern and inferred age of the lineages suggests Pliocene mountain building along active fault lines promoted their divergence 3.98-5.45 million years ago. A short interspersed nuclear element (SINE) polymorphism in the myosin gene intron (MYH-2) confirmed a pattern of restricted gene flow between lineages on either side of the mountain ranges associated with the Alpine Fault that runs southwest to northeast across the South Island of New Zealand. An analysis of molecular variance confirmed that approximately 40% of the genetic differentiation in O. n. polychroma is distributed across this major fault line. The straits between the main islands of New Zealand accounted for much less of the variation found within O. n. polychroma, most likely due to the repeated existence of landbridges between islands during periods of the Pleistocene that allowed migration. Overall, our findings reveal the relative roles of different climatic and geological processes, and in particular, demonstrate the importance of the Alpine Fault in the evolution of New Zealand's biota.