Evolution and phylogeny of the New Zealand cicada genus Kikihia Dugdale (Homoptera: Auchenorrhyncha: Cicadidae) with special reference to the origin of the Kermadec and Norfolk Islands' species

Aim Determine the phylogeny and dispersal patterns of the cicada genus Kikihia in New Zealand and the origin of the Norfolk, Kermadec, and Chatham Island cicadas. Location New Zealand, Norfolk Island, Kermadec Islands and Chatham Island. Methods DNA sequences from 16 species and four soon to be described species of cicadas from New Zealand and Norfolk Island (Australia) were examined. A total of 1401 base pairs were analysed from whole genome extraction of three mitochondrial genes (cytochrome oxidase subunit II, ATPase6 and ATPase8). These DNA sequences were aligned and analysed using standard likelihood approaches to phylogenetic analysis. Dates of divergences between clades were determined using a molecular clock based on Bayesian statistics. Results Most species in the genus Kikihia diverged between 3 and 5 million years ago (Ma) coincident with a period of rapid mountain building in New Zealand. Cicada species on the Kermadec and Norfolk Islands invaded recently from New Zealand and are closely related to the New Zealand North Island species Kikihia cutora. Main conclusions Speciation in the genus Kikihia was likely due in large part to the appearance of new habitats associated with the rise of the Southern Alps, starting c. 5 Ma. Dispersal of Kikihia species within mainland New Zealand probably occurred gradually rather than through long‐distance jumps. However, invasion of Norfolk, the Kermadecs and Chatham Islands had to have occurred through long‐distance dispersal.     

Biogeography and phylogeny of the New Zealand cicada genera (Hemiptera: Cicadidae) based on nuclear and mitochondrial DNA data

Aim Determine the geographical and temporal origins of New Zealand cicadas. Location New Zealand, eastern Australia and New Caledonia. Methods DNA sequences from 14 species of cicadas from New Zealand, Australia, and New Caledonia were examined. A total of 4628 bp were analysed from whole genome extraction of four mitochondrial genes (cytochrome oxidase subunits I and II, and ribosomal 12S and 16S subunits) and one nuclear gene (elongation factor‐1 alpha). These DNA sequences were aligned and analysed using standard phylogenetic methods based primarily on the maximum likelihood optimality criterion. Dates of divergences between clades were determined using several molecular clock methods. Results New Zealand cicadas form two well‐defined clades. One clade groups with Australian taxa, the other with New Caledonian taxa. The molecular clock analyses indicate that New Zealand genera diverged from the Australian and New Caledonian genera within the last 11.6 Myr. Main conclusions New Zealand was likely colonized by two or more invasions. One NZ lineage has its closest relatives in Australia and the other in New Caledonia. These invasions occurred well after New Zealand became isolated from other land masses, therefore cicadas must have crossed large bodies of water to reach New Zealand.     

Differentiating between hypotheses of lineage sorting and introgression in New Zealand alpine cicadas (Maoricicada Dugdale)

Lineage sorting and introgression can lead to incongruence among gene phylogenies, complicating the inference of species trees for large groups of taxa that have recently and rapidly radiated. In addition, it can be difficult to determine which of these processes is responsible for this incongruence. We explore these issues with the radiation of New Zealand alpine cicadas of the genus Maoricicada Dugdale. Gene trees were estimated from four putative independent loci: mitochondrial DNA (2274 nucleotides), elongation factor 1-alpha (1275 nucleotides), period (1709 nucleotides), and calmodulin (678 nucleotides). We reconstructed phylogenies using maximum likelihood and Bayesian methods from 44 individuals representing the 19 species and subspecies of Maoricicada and two outgroups. Species-level relationships were reconstructed using a novel extension of gene tree parsimony, whereby gene trees were weighted by their Bayesian posterior probabilities. The inferred gene trees show marked incongruence in the placement of some taxa, especially the enigmatic forest and scrub dwelling species, M. iolanthe. Using the species tree estimated by gene tree parsimony, we simulated coalescent gene trees in order to test the null hypothesis that the nonrandom placement of M. iolanthe among gene trees has arisen by chance. Under the assumptions of constant population size, known generation time, and panmixia, we were able to reject this null hypothesis. Furthermore, because the two alternative placements of M. iolanthe are in each case with species that share a similar song structure, we conclude that it is more likely that an ancient introgression event rather than lineage sorting has caused this incongruence.     

Surviving glacial ages within the Biotic Gap: phylogeography of the New Zealand cicada Maoricicada campbelli

Aim New Zealand is an ideal location in which to investigate the roles of landscape and climate change on speciation and biogeography. An earlier study of the widespread endemic cicada Maoricicada campbelli (Myers) found two phylogeographically distinguishable major clades – northern South Island plus North Island (northern‐SI + NI) and Otago. These two clades appeared to have diverged on either side of an area of the South Island known as the Biotic Gap. We sampled more intensively to test competing theories for this divergence. We aimed to discover if M. campbelli had survived within the Biotic Gap during recent glacial maxima, and if predicted areas of secondary contact between the two major clades existed. Location New Zealand. Methods We analysed mitochondrial DNA sequences (1520 bp; 212 individuals; 91 populations) using phylogenetic (maximum likelihood, Bayesian), population genetic (analysis of molecular variance) and molecular dating methods (Bayesian relaxed clock with improved priors). Results We found strong geographical structuring of genetic variation. Our dating analyses suggest that M. campbelli originated 1.83–2.58 Ma, and split into the two major clades 1.45–2.09 Ma. The main subclades in the northern‐SI + NI clade arose almost simultaneously at 0.69–1.03 Ma. Most subclades are supported by long internal branches and began to diversify 0.40–0.78 Ma. We found four narrow areas of secondary contact between the two major clades. We also found a difference between calling songs of the Otago vs. northern‐SI + NI clades. Main conclusions Phylogeographical patterns within M. campbelli indicate an early Pleistocene split into two major clades, followed by late Pleistocene range expansion and in situ population differentiation of subclades. The northern‐SI + NI clade diversified so rapidly that the main subclade relationships cannot be resolved, and we now have little evidence for a disjunction across the Biotic Gap. Structure within the main subclades indicates rapid divergence after a common bottlenecking event, perhaps attributable to an extremely cold glacial maximum at c. 0.43 Ma. Clade structure and dating analyses indicate that M. campbelli survived in many refugia during recent glacial maxima, including within the Biotic Gap. The narrow overlap between the two major clades is attributed to recent contact during the current interglacial and slow gene diffusion. The two major clades appear to be in the early stages of speciation based on genetic and behavioural differences.     

Phylogeography of the New Zealand cicada Maoricicada campbelli based on mitochondrial DNA sequences: ancient clades associated with cenozoic environmental change

New Zealand's isolation, its well-studied rapidly changing landscape, and its many examples of rampant speciation make it an excellent location for studying the process of genetic differentiation. Using 1520 base pairs of mitochondrial DNA from the cytochrome oxidase subunit I, ATPase subunits 6 and 8 and tRNA(Asp) genes, we detected two well-differentiated, parapatrically distributed clades within the widespread New Zealand cicada species Maoricicada campbelli that may prove to represent two species. The situation that we uncovered is unusual in that an ancient lineage with low genetic diversity is surrounded on three sides by two recently diverged lineages. Using a relaxed molecular clock model coupled with Bayesian statistics, we dated the earliest divergence within M. campbelli at 2.3 +/- 0.55 million years. Our data suggest that geological and climatological events of the late Pliocene divided a once-widespread species into northern and southern components and that near the middle of the Pleistocene the northern lineage began moving south eventually reaching the southern clade. The southern clade seems to have moved northward to only a limited extent. We discovered five potential zones of secondary contact through mountain passes that will be examined in future work. We predict that, as in North American periodical cicadas, contact between these highly differentiated lineages will exist but will not involve gene flow.     

New Zealand cicadas of the genus Maoricicada (Homoptera: Tibicinidae)

Abstract

The genus Maoricicada Dugdale, 1972 includes 19 taxa, of which 15 are described here. Ten from the South Island mountains are described as new: M. alticola, M. clamitans, M. mangu gourlayi, M. m. multicostata, M. m. celer, M. nigra frigida, M. otagoensis otagoensis, M. o. maceweni, M. phaeoptera, and M. tenuis. The distributions of most taxa are mapped, their songs are represented by sonograms, and their known periods of emergence are recorded. Evidence of hybridisation is unknown. Most species are characteristic of open, unforested sites; M. nigra nigra and M. oromelaena occur close to the summer snowline, a zone not known to be inhabited by cicadas elsewhere in the world.     

Diversification of Chionochloa (Poaceae) and biogeography of the New Zealand Southern Alps

Aim We test hypotheses regarding the origin of diversity and patterns of species richness in and around the New Zealand Southern Alps with 25 species of Chionochloa (Poaceae, Danthonioideae). Location New Zealand. Methods We inferred a well‐resolved and mostly robustly supported chloroplast phylogeny based on multiple DNA sequence markers (trnT–L–F, rpl16, trnD–psbM, atpB–rbcL, matK and ndhF), sampling 92% of the recognized species and 82% of the subspecific taxa. Nuclear ribosomal internal transcribed spacer sequences were also sampled, but proved uninformative. Biogeographic reconstruction and character optimization were done using both parsimony and likelihood approaches, and molecular dating used relaxed clock approaches. Results Most of the species diversity in Chionochloa stemmed from a common ancestor in the southern South Island with subsequent dispersal between areas. One clade of apparently cryptic taxa diversified within the central South Island ‘endemism gap’, persisting there throughout at least the latter half of the Pleistocene. Exclusively alpine and other habitat specialist species originated independently, the former relatively recently (between 7.6 Ma and the present). Main conclusions The phylogeny of Chionochloa and other published phylogenies of New Zealand plant groups demonstrate that the higher degree of endemism in the north and south of the New Zealand South Island relative to a central endemism gap cannot be explained by Alpine Fault displacement. Furthermore, our results suggest that if extinctions resulting from glaciations played a role in the origin of the central endemism gap, their impact was less than might be presumed on the basis of the distribution of taxa as they are currently defined. The diversification of Chionochloa and a number of New Zealand plant groups, such as Ranunculus, was contemporaneous with the initiation of the uplift of the Southern Alps. In contrast to patterns of diversifications within the alpine regions typical of the hyperdiverse Andes, exclusively alpine species in New Zealand arose independently from ancestors distributed in more lowland areas. Similarly, habitat specialists in Chionochloa arose independently from more generalist ancestors. Thus, although diversification in these groups may have been stimulated by mountain building and Pleistocene climatic oscillations, cladogenesis did not occur within the high alpine habitat itself.     

Biogeography and evolution of Dolichandra (Bignonieae,Bignoniaceae)

Little information on evolutionary relationships of Neotropical organisms or on the factors that have shaped the diversity currently encountered in this region is available. However, it is clear that biotic interactions and abiotic aspects have played important roles for species diversification in the region. This study focuses on Dolichandra (Bignonieae, Bignoniaceae), a clade of Neotropical lianas that is distributed broadly across different habitats and with diverse pollination and dispersal systems. We used sequences from two plastid DNA markers (ndhF and rpl32‐trnL) and one nuclear gene (PepC) to infer phylogenetic relationships in Dolichandra using parsimony and Bayesian approaches. We then used this phylogenetic framework as basis to study the biogeographic history, reconstruct the evolution of morphological characters and test the impact of morphology and environment on the diversification of the genus. More specifically, we: (1) time‐calibrate the phylogenetic tree of Dolichandra; (2) estimate the ancestral areas of the various lineages; (3) estimate the ancestral states of discrete and continuous morphological traits; (4) test for phylogenetic signal in environmental and phenotypic data; and (5) test whether morphological characters and/or niche evolution are correlated with cladogenesis. All Dolichandra spp. are monophyletic in the combined molecular phylogeny; relationships among species are generally well resolved, although poorly supported in some instances. The genus is inferred to have originated 36.43–26.23 Mya, possibly in eastern South America. Ancestral state reconstructions of continuous and discrete floral characters inferred a mixed morphology as the ancestral condition for the group. Phylogenetic signal differed between perianth and sexual whorls and gradual evolution was recovered for all traits except style length and anther length. Environmental variables showed no phylogenetic signal and a pattern of variation that was not correlated with branch length, suggesting that environmental transitions were concomitant with speciation. Dispersal is inferred to be the main driver of the differential distribution observed among species. In addition, climatic preferences and floral characters seem to have been important reproductive barriers in Dolichandra. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 179 , 403–420.     

Conflicting patterns of DNA barcoding and taxonomy in the cicada genus Tettigettalna from southern Europe (Hemiptera: Cicadidae)

DNA barcodes have great potential to assist in species identification, especially when high taxonomical expertise is required. We investigated the utility of the 5′ mitochondrial cytochrome c oxidase I (COI) region to discriminate between 13 European cicada species. These included all nine species currently recognized under the genus Tettigettalna, from which seven are endemic to the southern Iberian Peninsula. These cicadas have species‐specific male calling songs but are morphologically very similar. Mean COI divergence between congeners ranged from 0.4% to 10.6%, but this gene was proven insufficient to determine species limits within genus Tettigettalna because a barcoding gap was absent for several of its species, that is, the highest intraspecific distance exceeded the lowest interspecific distance. The genetic data conflicted with current taxonomic classification for T. argentata and T. mariae. Neighbour‐joining and Bayesian analyses revealed that T. argentata is geographically structured (clades North and South) and might constitute a species complex together with T. aneabi and T. mariae. The latter diverges very little from the southern clade of T. argentata and shares with it its most common haplotype. T. mariae is often in sympatry with T. argentata but it remains unclear whether introgression or incomplete lineage sorting may be responsible for the sharing of haplotypes. T. helianthemi and T. defauti also show high intraspecific variation that might signal hidden cryptic diversity. These taxonomic conflicts must be re‐evaluated with further studies using additional genes and extensive morphological and acoustic analyses.     

Giants and dwarfs: molecular phylogenies reveal multiple origins of annual spurges within Euphorbia subg. Esula

Euphorbia (Euphorbiaceae) comprises over 2150 species and is thus the second-largest genus of flowering plants. In Europe, it is represented by more than 100 species with highest diversity in the Mediterranean area; the majority of taxa belong to subgenus Esula Pers., including about 500 taxa. The few available phylogenetic studies yielded contrasting results regarding the monophyly of subg. Esula, and the phylogenetic relationships among its constituents remain poorly understood. We have sampled DNA sequences from the nuclear ribosomal internal transcribed spacer (ITS) and the plastid trnT-trnF region from about 100, predominantly European taxa of subg. Esula in order to infer its phylogenetic history. The plastid data support monophyly of subg. Esula whereas the ITS phylogeny, which is generally less resolved, is indecisive in this respect. Although some major clades have partly incongruent positions in the ITS and plastid phylogenies, the taxonomic content of the major terminal clades is congruent in both trees. As traditional sectional delimitations are largely not corroborated, an improved classification is proposed. Character state reconstruction illustrates that the annual life form developed independently several times in different clades of subgenus Esula from perennial ancestors, and that several morphological traits used in previous classifications of Euphorbia developed in parallel in different lineages.     

Specificity of cicada calling songs in the genus Tibicina (Hemiptera: Cicadidae)

Abstract.  Males of Tibicina cicada species produce a sustained and monotonous calling song by tymbal activity. This acoustic signal constitutes the first step in pair formation, attracting females at long range, and is involved in male–male interactions. The specificity of this signal was investigated for the first time for seven species and one subspecies of Tibicina occurring in France. This analysis was achieved by describing tymbal anatomy, tymbal mechanism and calling song structure. Male calling songs are emitted following the same general scheme: tymbals are activated alternately and the successive buckling of the sclerotized ribs that they bear produces a regular succession of groups of pulses. The structural and mechanical properties shared by Tibicina species and subspecies lead to a considerable uniformity of the signal shape. Nevertheless, a principal component analysis applied to eight temporal and three frequency parameters revealed differences between the signals of the species studied. In particular, calling songs differed in groups of pulse rate and/or in peak of the second frequency band (carrier frequency). These acoustic differences are probably linked to differences in the numbers of tymbal ribs and body size. Groups of pulse rate and/or peak of the second frequency band could encode specific information. However, Tibicina calling songs may not act as distinct specific-mate recognition systems and may not play a leading role in the mating isolation process; rather, they might merely belong to a complex set of specific spatial, ecological, ethological and morphological characters that ensure syngamy.     

Untangling the assassin's web: Phylogeny and classification of the spider-associated Emesine complex (Hemiptera: Reduviidae)

Web-building spiders are formidable predators, yet assassin bugs in the Emesine Complex (Hemiptera: Reduviidae: Emesinae, Saicinae, and Visayanocorinae) prey on spiders. The Emesine Complex comprises >1000 species and these web-associated predatory strategies may have driven their diversification. However, lack of natural history data and a robust phylogenetic framework currently preclude tests of this hypothesis. We combine Sanger (207 taxa, 3865 bp) and high-throughput sequencing data (15 taxa, 381 loci) to generate the first taxon- and data-rich phylogeny for this group. We discover rampant paraphyly among subfamilies and tribes, necessitating revisions to the classification. We use ancestral character state reconstructions for 40 morphological characters to identify diagnostic features for a revised classification. Our new classification treats Saicinae Stål and Visayanocorinae Miller as junior synonyms of Emesinae Amyot and Serville, synonymizes the emesine tribes Ploiariolini Van Duzee and Metapterini Stål with Emesini Amyot and Serville, and recognises six tribes within Emesinae (Collartidini Wygodzinsky, Emesini, Leistarchini Stål, Oncerotrachelini trib.n. , Saicini Stål stat.n. , and Visayanocorini Miller stat.n. ). We show that a pretarsal structure putatively involved in web-associated behaviours evolved in the last common ancestor of Emesini, the most species-rich clade within Emesinae, suggesting that web-associations could be widespread in Emesinae.     

The diversification of the genus Monodelphis and the chronology of Didelphidae (Didelphimorphia)

To evaluate the monophyletic status of the genus Monodelphis, and its species complexes, we used a 9.3‐kb multimarker alignment to build a phylogenetic tree based on the largest taxon sampling for this didelphid genus to date. Furthermore, using a Bayesian framework and six calibration points, we inferred the divergence times for the major Monodelphis lineages and their current geographical distribution to perform an ancestral state reconstruction for geographical areas. Our results indicate the monophyletic nature of Monodelphis and suggest ‘kunsi’ as a new species complex that includes Monodelphis kunsi and an undescribed species. Monodelphis is further subdivided into three clades showing a common vicariance pattern, with each major clade consisting of a northern South American lineage joined with an Atlantic Forest lineage. This geographic consistency suggests a vicariant event that might have been related to a warm period at the Oligocene/Miocene border, according to our time results. © 2015 The Linnean Society of London     

Phylogenetic relationships within the New World subfamily Larreoideae (Zygophyllaceae) confirm polyphyly of the disjunct genus Bulnesia

The Larreoideae subfamily is the major representative of the family Zygophyllaceae in South America, where several of its members are common to dominant in arid regions of the Southern Cone. However, there are currently no phylogenetic analyses of the subfamily that may help to understand its origin and diversification. Additionally, there are taxonomic discrepancies around Bulnesia Gay (1845), one of its more important genera. Accordingly, we performed a phylogenetic analysis combining chloroplast (rbcL and trnL-F) and nuclear (ITS) DNA sequences. Bayesian and Parsimony analyses were performed to highlight the intergeneric relationships within Larreoideae. All genera with the exception of Bulnesia are monophyletic and we propose to redefine Bulnesia, dividing it in two genera. Furthermore, other taxonomic issues of the remaining genera are solved. This study represents the first approximation to clarify the phylogenetic relationships amongst all Larreoideae genera, producing a phylogenetic framework that can be used in future macro-ecological studies.     

Support for vicariant origins of the New Zealand Onychophora

Aim The distribution of Onychophora across the southern continents has long been considered the result of vicariance events. However, it has recently been hypothesized that New Zealand was completely inundated during the late Oligocene (25–22 Ma) and therefore that the entire biota is the result of long-distance dispersal. We tested this assumption using phylogenetic and molecular dating of DNA sequence data from Onychophora. Location New Zealand, Australia, South Africa, Chile (South America). Methods We obtained DNA sequence data from the nuclear genes 28S and 18S rRNA to reconstruct relationships among species of Peripatopsidae (Onychophora). We performed molecular dating under a Bayesian relaxed clock model with a range of prior distributions using the rifting of South America and South Africa as a calibration. Results Our phylogenetic trees revealed that the New Zealand genera Ooperipatellus and Peripatoides, together with selected Australian genera (Euperipatoides, Phallocephale and an undescribed genus from Tasmania), form a monophyletic group that is the sister group to genera from Chile (Metaperipatus) and South Africa (Peripatopsis and Opisthopatus). The relaxed clock dating analyses yielded mean divergence times from 71.3 to 78.9 Ma for the split of the New Zealand Peripatoides from their Australian sister taxa. The 0.95 Bayesian posterior intervals were very broad and ranged from 24.5 to 137.6 Ma depending on the prior assumptions. The mean divergence of the New Zealand species of Ooperipatellus from the Australian species Ooperipatellus insignis was estimated at between 39.9 and 46.2 Ma, with posterior intervals ranging from 9.5 to 91.6 Ma. Main conclusions The age of Peripatoides is consistent with long-term survival in New Zealand and implies that New Zealand was not completely submerged during the Oligocene. Ooperipatellus is less informative on the question of continuous land in the New Zealand region because we cannot exclude a post-Oligocene divergence. The great age of Peripatoides is consistent with a vicariant origin of this genus resulting from the rifting of New Zealand from the eastern margin of Gondwana and supports the assumptions of previous authors who considered the Onychophora to be a relict component of the New Zealand biota.     

Historical biogeography of the livebearing fish genus Poeciliopsis (Poeciliidae: Cyprinodontiformes)

To assess the historical biogeography of freshwater topminnows in the genus Poeciliopsis, we examined sequence variation in two mitochondrial genes, cytochrome b (1140 bp) and NADH subunit 2 (1047 bp). This widespread fish genus is distributed from Arizona to western Colombia, and nearly half of its 21 named species have distributions that border on the geologically active Trans-Mexican Volcanic Belt (TMVB), a region that defines the uplifted plateau (Mesa Central) of Mexico. We used the parametric bootstrap method to test the hypothesis that a single vicariant event associated with the TMVB was responsible for divergence of taxa found to the north and south of this boundary. Because the single-event hypothesis was rejected as highly unlikely, we hypothesize that at least two geological events were responsible for divergence of these species. The first (8-16 million years ago) separated ancestral populations that were distributed across the present TMVB region. A second event (2.8-6.4 million years ago) was associated with northward dispersal and subsequent vicariance of two independent southern lineages across the TMVB. The geological history of this tectonically and volcanically active region is discussed and systematic implications for the genus are outlined.     

Phylogeny,biogeography, and morphological ancestral character reconstruction in the Mediterranean genus Fumana (Cistaceae)

Fumana is a diverse genus of the Cistaceae family, consisting of 21 currently accepted species. In this study, nuclear (ITS) and plastid (matK, trnT‐L) molecular markers were used to reconstruct the phylogeny and to estimate divergence times, including 19 species of Fumana. Phylogenetic analyses (Bayesian Inference, Maximum Parsimony and Maximum Likelihood) confirmed the monophyly of Fumana and did not support the infrageneric divisions previously established. The results support four main clades that group species that differ in vegetative and reproductive characters. Given the impossibility to define morphological characters common to all species within the clades, our proposal is to reject infrageneric divisions. Molecular dating and ancestral area analyses provide evidence for a Miocene diversification of the genus in the north‐western Mediterranean. Ancestral state reconstructions revealed ancestral character states for some traits related to xeric and arid habitats, suggesting a preadaptation to the Mediterranean climate.     

A new genus and species of tettigarctid cicada from the early Miocene of New Zealand: Paratettigarcta zealandica (Hemiptera,Auchenorrhyncha, Tettigarctidae)

A new genus and species of primitive cicada (Hemiptera: Tettigarctidae) is described from the early Miocene of southern New Zealand. Paratettigarcta zealandica gen. et sp. n. is the first cicada (Cicadoidea) fossil from New Zealand and exhibits wing venation patterns typical for the subfamily Tettigarctinae. It differs from other fossil taxa and the extant genus Tettigarcta in the early divergence of CuA2 from the nodal line in the forewing, its parallel-sided subcostal cell, the early bifurcation of vein M and long apical cells of the hindwing, and in wing pigmentation patterns.