Detecting and Removing Sample Contamination in Phylogenomic Data: An Example and its Implications for Cicadidae Phylogeny (Insecta: Hemiptera).

Contamination of a genetic sample with DNA from one or more nontarget species is a continuing concern of molecular phylogenetic studies, both Sanger sequencing studies and next-generation sequencing studies. We developed an automated pipeline for identifying and excluding likely cross-contaminated loci based on the detection of bimodal distributions of patristic distances across gene trees. When contamination occurs between samples within a data set, a comparison between a contaminated sample and its contaminant taxon will yield bimodal distributions with one peak close to zero patristic distance. This new method does not rely on a priori knowledge of taxon relatedness nor does it determine the causes(s) of the contamination. Exclusion of putatively contaminated loci from a data set generated for the insect family Cicadidae showed that these sequences were affecting some topological patterns and branch supports, although the effects were sometimes subtle, with some contamination-influenced relationships exhibiting strong bootstrap support. Long tip branches and outlier values for one anchored phylogenomic pipeline statistic (AvgNHomologs) were correlated with the presence of contamination. While the anchored hybrid enrichment markers used here, which target hemipteroid taxa, proved effective in resolving deep and shallow level Cicadidae relationships in aggregate, individual markers contained inadequate phylogenetic signal, in part probably due to short length. The cleaned data set, consisting of 429 loci, from 90 genera representing 44 of 56 current Cicadidae tribes, supported three of the four sampled Cicadidae subfamilies in concatenated-matrix maximum likelihood (ML) and multispecies coalescent-based species tree analyses, with the fourth subfamily weakly supported in the ML trees. No well-supported patterns from previous family-level Sanger sequencing studies of Cicadidae phylogeny were contradicted. One taxon (Aragualna plenalinea) did not fall with its current subfamily in the genetic tree, and this genus and its tribe Aragualnini is reclassified to Tibicininae following morphological re-examination. Only subtle differences were observed in trees after the removal of loci for which divergent base frequencies were detected. Greater success may be achieved by increased taxon sampling and developing a probe set targeting a more recent common ancestor and longer loci. Searches for contamination are an essential step in phylogenomic analyses of all kinds and our pipeline is an effective solution. [Auchenorrhyncha; base-composition bias; Cicadidae; Cicadoidea; Hemiptera; phylogenetic conflict.].

How the Aridification of Australia Structured the Biogeography and Influenced the Diversification of a Large Lineage of Australian Cicadas.

Over the last 30 million years, Australia's landscape has undergone dramatic cooling and drying due to the establishment of the Antarctic Circumpolar Current and change in global CO$_{2}$ levels. Studies have shown that many Australian organisms went extinct during these major cooling events, while others experienced adaptive radiations and increases in diversification rates as a result of exploiting new niches in the arid zone. Despite the many studies on diversification and biogeography in Australia, few have been continent-wide and none have focused on a group of organisms adapted to feeding on plants. We studied 162 species of cicadas in the Australian Pauropsalta complex, a large generic lineage within the tribe Cicadettini. We asked whether there were changes in the diversification rate of Pauropsalta over time and if so: 1) which clades were associated with the rate change? 2) did timing of rate shifts correspond to known periods of dramatic historical climate change, 3) did increases in diversification rate along select lineages correspond to adaptive radiations with movement into the arid zone? To address these questions, we estimated a molecular phylogeny of the Pauropsalta complex using ${\sim}$5300 bp of nucleotide sequence data distributed among five loci (one mtDNA locus and four nDNA loci). We found that this large group of cicadas did not diversify at a constant rate as they spread through Australia; instead the signature of decreasing diversification rate changed roughly around the time of the expansion of the east Antarctic ice sheets ${\sim}$16 Ma and the glaciation of the northern hemisphere ${\sim}$3 Ma. Unlike other Australian taxa, the Pauropsalta complex did not explosively radiate in response to an early invasion of the arid zone. Instead multiple groups invaded the arid zone and experienced rates of diversification similar to mesic-distributed taxa. We found evidence for relictual groups, located in pre-Mesozoic habitat, that have not diversified and continue to reside on mesic hosts in isolated "habitat islands". Future work should focus on groups of similar ages with similar distribution patterns to determine whether this tempo and pattern of diversification and biogeography is consistent with evidence from other phytophagous insects.

Inflation of Molecular Clock Rates and Dates: Molecular Phylogenetics, Biogeography, and Diversification of a Global Cicada Radiation from Australasia (Hemiptera: Cicadidae: Cicadettini).

Dated phylogenetic trees are important for studying mechanisms of diversification, and molecular clocks are important tools for studies of organisms lacking good fossil records. However, studies have begun to identify problems in molecular clock dates caused by uncertainty of the modeled molecular substitution process. Here we explore Bayesian relaxed-clock molecular dating while studying the biogeography of ca. 200 species from the global cicada tribe Cicadettini. Because the available fossils are few and uninformative, we calibrate our trees in part with a cytochrome oxidase I (COI) clock prior encompassing a range of literature estimates for arthropods. We show that tribe-level analyses calibrated solely with the COI clock recover extremely old dates that conflict with published estimates for two well-studied New Zealand subclades within Cicadettini. Additional subclade analyses suggest that COI relaxed-clock rates and maximum-likelihood branch lengths become inflated relative to EF-1[Formula: see text] intron and exon rates and branch lengths as clade age increases. We present corrected estimates derived from: (i) an extrapolated EF-1[Formula: see text] exon clock derived from COI-calibrated analysis within the largest New Zealand subclade; (ii) post hoc scaling of the tribe-level chronogram using results from subclade analyses; and (iii) exploitation of a geological calibration point associated with New Caledonia. We caution that considerable uncertainty is generated due to dependence of substitution estimates on both the taxon sample and the choice of model, including gamma category number and the choice of empirical versus estimated base frequencies. Our results suggest that diversification of the tribe Cicadettini commenced in the early- to mid-Cenozoic and continued with the development of open, arid habitats in Australia and worldwide. We find that Cicadettini is a rare example of a global terrestrial animal group with an Australasian origin, with all non-Australasian genera belonging to two distal clades. Within Australia, we show that Cicadettini is more widely distributed than any other cicada tribe, diverse in temperate, arid and monsoonal habitats, and nearly absent from rainforests. We comment on the taxonomic implications of our findings for thirteen cicada genera.

The phylogenetic utility of acetyltransferase (ARD1) and glutaminyl tRNA synthetase (QtRNA) for reconstructing Cenozoic relationships as exemplified by the large Australian cicada Pauropsalta generic complex.

The Pauropsalta generic complex is a large group of cicadas (72 described spp.; >82 undescribed spp.) endemic to Australia. No previous molecular work on deep level relationships within this complex has been conducted, but a recent morphological revision and phylogenetic analysis proposed relationships among the 11 genera. We present here the first comprehensive molecular phylogeny of the complex using five loci (1 mtDNA, 4 nDNA), two of which are from nuclear genes new to cicada systematics. We compare the molecular phylogeny to the morphological phylogeny. We evaluate the phylogenetic informativeness of the new loci to traditional cicada systematics loci to generate a baseline of performance and behavior to aid in gene choice decisions in future systematic and phylogenomic studies. Our maximum likelihood and Bayesian inference phylogenies strongly support the monophyly of most of the newly described genera; however, relationships among genera differ from the morphological phylogeny. A comparison of phylogenetic informativeness among all loci revealed that COI 3rd positions dominate the informativeness profiles relative to all other loci but exhibit some among taxon nucleotide bias. After removing COI 3rd positions, COI 1st positions dominate near the terminals, while the period intron has the most phylogenetic informativeness near the root. Among the nuclear loci, ARD1 and QtRNA have lower phylogenetic informativeness than period intron and elongation factor 1 alpha intron, but the informativeness increases at you move from the tips to the root. The increase in phylogenetic informativeness deeper in the tree suggests these loci may be useful for resolving older relationships.

Molecular phylogenetics, diversification, and systematics of Tibicen Latreille 1825 and allied cicadas of the tribe Cryptotympanini, with three new genera and emphasis on species from the USA and Canada(Hemiptera: Auchenorrhyncha: Cicadidae).

North America has a diverse cicada fauna with multiple genera from all three Cicadidae subfamilies, yet molecular phylogenetic analyses have been completed only for the well-studied periodical cicadas (Magicicada Davis). The genus Tibicen Latreille, a large group of charismatic species, is in need of such work because morphological patterns suggest multiple groups with complicated relationships to other genera in the tribe Cryptotympanini. In this paper we present a molecular phylogenetic analysis, based on mitochondrial and nuclear DNA, of 35 of the 38 extant USA species and subspecies of the genus Tibicen together with their North American tribal allies (Cornuplura Davis, Cacama Davis), selected Tibicen species from Eurasia, and representatives of other Eurasian and Pacific cryptotympanine genera. This tree shows that Tibicen contains several well-supported clades, one predominating in eastern and central North America and related to Cryptotympana Stål and Raiateana Boulard, another in western North America related to Cacama and Cornuplura, and at least two clades in Eurasia. We also present a morphological cladistic analysis of Tibicen and its close allies based on 27 characters. Character states identified in the cladistic analysis define three new genera, two for North American taxa (Hadoa gen. n. and Neotibicen gen. n.) including several Mexican species, and one for Asian species (Subsolanus gen. n.). Using relaxed molecular clocks and literature-derived mtDNA rate estimates, we estimate the timeframe of diversification of Tibicen clades and find that intergeneric divergence has occurred since the late Eocene, with most extant species within the former Tibicen originating after the mid-Miocene. We review patterns of ecology, behavior, and geography among Tibicen clades in light of the phylogenetic results and note that the study of these insects is still in its early stages. Some Mexican species formerly placed in Tibicen are here transferred to Diceroprocta, following refinement of the definition of that genus.

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.

Hybridization, mitochondrial DNA phylogeography, and prediction of the early stages of reproductive isolation: lessons from New Zealand cicadas (genus Kikihia).

One of the major tenets of the modern synthesis is that genetic differentiation among subpopulations is translated over time into genetic differentiation among species. Phylogeographic exploration is therefore essential to the study of speciation because it can reveal the presence of subpopulations that may go on to become species or that may already represent cryptic species. Acoustic species-specific mating signals provide a significant advantage for the recognition of cryptic or incipient species. Because the majority of species do not have such easily recognized premating signals, data from acoustically signaling species can serve as a valuable heuristic tool. Acoustic signals are also convenient tools for recognizing hybridization events. Here, we demonstrate that evidence of hybridization in the form of intermediate song phenotypes is present in many contact zones between species of the New Zealand grass cicadas of the Kikihia muta species complex and that recurring mitochondrial DNA (mtDNA) introgression has created misleading patterns that make it difficult to identify certain taxa using song or mtDNA alone. In one case, introgression appears to have occurred between allopatric taxa by dispersal of introgressed populations of an intermediary species ("hybridization by proxy"). We also present a comparison of mtDNA-tree- and song-based taxonomies obtained for the K. muta complex. We find that 12 mtDNA candidate species are identified using shifts in phylogenetic branching rate found by a single-threshold mixed Yule-coalescent lineage model, while only 7 candidate species are identified using songs. Results from the Yule-coalescent model are dependent on factors such as the number of modeled thresholds and the inclusion of duplicate haplotypes. Genetic distances within song species reach a maximum at about 0.028 substitutions/site when likely cases of hybridization and introgression are excluded. Large genetic breaks or "gaps" are not observed between some northern (warmer climate) song clades, possibly because climate-induced bottlenecks have been less severe. These results support ongoing calls for multimarker genetic studies as well as "integrative taxonomy" that combines information from multiple character sources, including behavior, ecology, geography, and morphology.

New genera and new species of Western Australian cicadas (Hemiptera: Cicadidae).

The cicada fauna of Western Australia is briefly reviewed. Six genera and 14 species are recorded from the State for the first time bringing the total of known species and subspecies to 105 and a list of all 105 is provided. Among the taxa here recorded are five new genera and 13 new species belonging to the tribes Macrotristriini (Illyria viridis sp. n.), Pictilini (Chrysocicada trophis sp. n.), and Cicadettini (Calipsalta gen. n., Calipsalta brunnea sp. n., C. fumosa sp. n., C. viridans sp. n., Kalarko gen. n., Kalarko ferruginosus sp. n., Ewartia adusta sp. n., Parvopsalta gen. n., Parvopsalta victoriae sp. n., Pedana gen. n., Pedana hesperia sp. n., Pegapsaltria gen. n., Pegapsaltria lutea sp. n., Pyropsalta amnica sp. n., Py. patula sp. n., and Py. rhythmica sp. n). In addition, Erempsalta hermannsburgensis (Distant, 1907) is redescribed and its presence in Western Australia (and four other States) documented for the first time. Songs are analysed for all species except two species of Pyropsalta where recordings were unavailable.

Cryptic failure of partitioned Bayesian phylogenetic analyses: lost in the land of long trees.

Partitioned Bayesian phylogenetic analyses of routine genetic data sets, constructed using MrBayes (Ronquist and Huelsenbeck 2003), can become trapped in regions of parameter space characterized by unrealistically long trees and distorted partition rate multipliers. Such analyses commonly fail to reach stationarity during hundreds of millions of generations of sampling-many times longer than most published analyses. Some data sets are so prone to this problem that paired MrBayes runs begun from different starting trees repeatedly find the same incorrect long-tree solutions and consequently pass the most commonly employed tests of stationarity, including the average standard deviation of split frequencies (ASDSF) and the potential scale reduction factor (PSRF) statistics offered by MrBayes (Gelman and Rubin 1992). In these situations, failure to reach stationarity is recognizable only in light of prior knowledge of model parameters, such as the expectation that third-codon-position sites usually evolve fastest in protein-coding genes. The conditions that lead to the long-tree problem are frequently encountered in phylogenetic studies today, and I present 6 demonstration examples from the literature. Although the effects on tree length (TL) are often dramatic, effects on topology appear to be subtle. Susceptibility to the problem is sometimes predicted by the difference between the true TL and the starting TL. In some cases, the problems described here can be avoided or reduced by manipulation of the starting TL and/or by adjustments to the prior on branch lengths. In more difficult situations, accurate branch length estimation may not be possible with Bayesian methods because of dependence of the solution on the branch length prior.

Redefinition of the cicada tribe Hemidictyini Distant, 1905, status of the tribe Iruanini Boulard, 1993 rev. stat., and the establishment of Hovanini n. tribe and Sapantangini n. tribe (Hemiptera: Cicadidae).

A recent paper by Ruschel Campos (2019) on "leaf-winged" cicadas proposed a significant reorganization of the cicada tribe Hemidictyini Distant, 1905g, including synonymization of the monogeneric tribe Lacetasini Moulds Marshall, 2018 following the results of a cladistic parsimony analysis of morphological characters. In this study, we reconsider and revise the morphological analysis of Ruschel Campos and obtain new genetic data for Hemidictya. We find that their study suffers from a limited taxon sample, inappropriate outgroup selection, and misinterpretation of genitalic characters (uncus vs. claspers). We show that Hemidictyini sensu Ruschel Campos includes members of multiple tribes and subfamilies, and we conclude that some of the taxonomic transfers by Ruschel Campos are not supported. The two most similar and leaf-like cicadas, Hemidictya Burmeister, 1835 (South America) and Hovana Distant, 1905g (Madagascar), are probably not closely related but rather an excellent example of convergent evolution. Lacetasini is not a junior synonym of the Hemidictyini but a distinct part of the Tettigomyiinae Distant, 1905g as originally classified. We return or transfer the genera Lacetas Karsch, 1890, Iruana Distant, 1905g, Bafutalna Boulard, 1993, and Murphyalna Boulard, 2012 to the Lacetasini. With the transfer of all genera of Iruanina Boulard, 1993 and Bafutalnina Boulard, 1993 to Lacetasini and with Lacetas transferred to the Iruanina, Lacetasini n. syn. becomes a subjective junior synonym of Iruanini rev. stat. in the Tettigomyiinae. We assign Hovana to Hovanini n. tribe in the Tettigomyiinae and Sapantanga Distant, 1905g to Sapantangini n. tribe in the Tibicininae Distant, 1905b. We propose that Hemidictyini sensu novo contains only the genus Hemidictya and we assign the tribe to Tibicininae with a revised diagnosis.

Glacial refugia in a maritime temperate climate: cicada (Kikihia subalpina) mtDNA phylogeography in New Zealand.

Understanding the biological significance of Pleistocene glaciations requires knowledge of the nature and extent of habitat refugia during glacial maxima. An opportunity to examine evidence of glacial forest refugia in a maritime, Southern Hemisphere setting is found in New Zealand, where the extent of Pleistocene forests remains controversial. We used the mitochondrial phylogeography of a forest-edge cicada (Kikihia subalpina) to test the hypothesis that populations of this species survived throughout South Island during the Last Glacial Maximum. We also compared mitochondrial DNA phylogeographic patterns with male song patterns that suggest allopatric divergence across Cook Strait. Cytochrome oxidase I and II sequences were analyzed using network analysis, maximum-likelihood phylogenetic estimation, Bayesian dating and Bayesian skyline plots. K. subalpina haplotypes from North Island and South Island form monophyletic clades that are concordant with song patterns. Song divergence corresponds to approximately 2% genetic divergence, and Bayesian dating suggests that the North Island and South Island population-lineages became isolated around 761,000 years BP. Almost all South Island genetic variation is found in the north of the island, consistent with refugia in Marlborough Sounds, central Nelson and northwest Nelson. All central and southern South Island and Stewart Island haplotypes are extremely similar to northern South Island haplotypes, a 'northern richness/southern purity' pattern that mirrors genetic patterns observed in many Northern Hemisphere taxa. Proposed southern South Island forest habitat fragments may have been too small to sustain populations of K. subalpina, and/or they may have harboured ecological communities with no modern-day analogues.

Steady Plio-Pleistocene diversification and a 2-million-year sympatry threshold in a New Zealand cicada radiation.

Estimation of diversification rates in evolutionary radiations requires a complete accounting of cryptic species diversity. The rapidly evolving songs of acoustically signaling insects make them good model organisms for such studies. This paper examines the timing of diversification of a large (30 taxon) group of New Zealand cicadas (genus Kikihia Dugdale). We use Bayesian relaxed-clock methods and phylogenetic trees based on nuclear and mitochondrial DNA data, and we apply alternative combinations of evolutionary rate priors and geological calibrations. The extant Kikihia taxa began to diversify near the Miocene/Pliocene boundary around the time of increased mountain-building, and both the mitochondrial and nuclear-gene trees confirm early splits of lineages currently represented by lowland forest-dwelling taxa. Most lineages originated in the Pleistocene, and sustained diversification occurred rapidly at over 0.5 lineages/my, a rate comparable to that of the Hawaiian silverswords. Diversification rate tests suggest an increase in the early to mid-Pliocene, followed by constant diversification from the Late Pliocene onward. No descendants of the many Pleistocene-age splits have evolved the ability to coexist in sympatry, and, where they do come into contact, hybrid zones have been documented based on acoustic and DNA evidence. In contrast, lineages separated in time by approximately 2Myr often overlap in distribution with no evidence of hybridization. This suggests that at least 2Myr has been required to achieve the level of divergence required for reproductive isolation.

A specialized fungal parasite (Massospora cicadina) hijacks the sexual signals of periodical cicadas (Hemiptera: Cicadidae: Magicicada).

Male periodical cicadas (Magicicada spp.) infected with conidiospore-producing ("Stage I") infections of the entomopathogenic fungus Massospora cicadina exhibit precisely timed wing-flick signaling behavior normally seen only in sexually receptive female cicadas. Male wing-flicks attract copulation attempts from conspecific males in the chorus; close contact apparently spreads the infective conidiospores. In contrast, males with "Stage II" infections that produce resting spores that wait for the next cicada generation do not produce female-specific signals. We propose that these complex fungus-induced behavioral changes, which resemble apparently independently derived changes in other cicada-Massospora systems, represent a fungus "extended phenotype" that hijacks cicadas, turning them into vehicles for fungus transmission at the expense of the cicadas' own interests.

A molecular phylogeny of the cicadas (Hemiptera: Cicadidae) with a review of tribe and subfamily classification.

A molecular phylogeny and a review of family-group classification are presented for 137 species (ca. 125 genera) of the insect family Cicadidae, the true cicadas, plus two species of hairy cicadas (Tettigarctidae) and two outgroup species from Cercopidae. Five genes, two of them mitochondrial, comprise the 4992 base-pair molecular dataset. Maximum-likelihood and Bayesian phylogenetic results are shown, including analyses to address potential base composition bias. Tettigarcta is confirmed as the sister-clade of the Cicadidae and support is found for three subfamilies identified in an earlier morphological cladistic analysis. A set of paraphyletic deep-level clades formed by African genera are together named as Tettigomyiinae n. stat. Taxonomic reassignments of genera and tribes are made where morphological examination confirms incorrect placements suggested by the molecular tree, and 11 new tribes are defined (Arenopsaltriini n. tribe, Durangonini n. tribe, Katoini n. tribe, Lacetasini n. tribe, Macrotristriini n. tribe, Malagasiini n. tribe, Nelcyndanini n. tribe, Pagiphorini n. tribe, Pictilini n. tribe, Psaltodini n. tribe, and Selymbriini n. tribe). Tribe Tacuini n. syn. is synonymized with Cryptotympanini, and Tryellina n. syn. is synonymized with an expanded Tribe Lamotialnini. Tribe Hyantiini n. syn. is synonymized with Fidicinini. Tribe Sinosenini is transferred to Cicadinae from Cicadettinae, Cicadatrini is moved to Cicadettinae from Cicadinae, and Ydiellini and Tettigomyiini are transferred to Tettigomyiinae n. stat from Cicadettinae. While the subfamily Cicadinae, historically defined by the presence of timbal covers, is weakly supported in the molecular tree, high taxonomic rank is not supported for several earlier clades based on unique morphology associated with sound production.

A new Neotibicen cicada subspecies (Hemiptera: Cicadidae) from the southeastern USA forms hybrid zones with a widespread relative despite a divergent male calling song.

A morphologically cryptic subspecies of Neotibicen similaris (Smith and Grossbeck) is described from forests of the Apalachicola region of the southeastern United States. Although the new form exhibits a highly distinctive male calling song, it hybridizes extensively where it meets populations of the nominate subspecies in parapatry, by which it is nearly surrounded. This is the first reported example of hybridization between North American nonperiodical cicadas. Acoustic and morphological characters are added to the original description of the nominate subspecies, and illustrations of complex hybrid song phenotypes are presented. The biogeography of N. similaris is discussed in light of historical changes in forest composition on the southeastern Coastal Plain.

Cicadas (Hemiptera: Cicadidae) from Camiguin of Mindanao Province and Dinagat Island in the Philippines, with a new genus and three new species.

This paper provides the first faunal checklist for the family Cicadidae (Hemiptera) from Camiguin of Mindanao Province and Dinagat Island in the Philippines, comprising ten species belonging to nine genera. Cryptotympana shillana Lee & Mohagan sp. nov., Orientopsaltria inermis (Stål, 1870), Purana crassinotata Lee, 2015, and Huechys parvula Haupt, 1924 are recorded for the first time from Camiguin. Platypleura dinagatensis Lee sp. nov., Chremistica kyoungheeae Lee, 2010, Dundubia vaginata (Fabricius, 1787), Oncotympana pallidiventris Stål, 1870, and Philipsalta nigrina Lee, Marshall & Hill sp. nov. are newly recorded from Dinagat Island. A new genus Philipsalta Lee, Marshall & Hill gen. nov. is erected. Huechysini Distant, 1905 syn. nov. is synonymized with Cicadettini Buckton, 1889. Information on geographic distributions of the Camiguin and Dinagat species is also provided.

Describing drivers of and barriers to persistence with tiotropium in patients with chronic obstructive pulmonary disease: a mixed-methods approach.

OBJECTIVE: To understand the factors influencing persistence with tiotropium in patients with chronic obstructive pulmonary disease (COPD). METHODS: Patients classified as 'persistent' or 'non-persistent' with tiotropium were identified from pharmacy dispensing records. Patients were compared for health status, beliefs and behaviours using data from questionnaires and interviews. KEY FINDINGS: Perceptions of the risks and benefits of medication, fear of worsening illness, and the GP's emphasis on the importance of the medication were key determinants of tiotropium persistence. CONCLUSIONS: Perceptions, attitudes and beliefs of patients and doctors influence persistence with tiotropium. These complex interactions need to be targeted to improve persistence with medicines in COPD.

Versatile aggressive mimicry of cicadas by an Australian predatory katydid.

BACKGROUND: In aggressive mimicry, a predator or parasite imitates a signal of another species in order to exploit the recipient of the signal. Some of the most remarkable examples of aggressive mimicry involve exploitation of a complex signal-response system by an unrelated predator species. METHODOLOGY/PRINCIPAL FINDINGS: We have found that predatory Chlorobalius leucoviridis katydids (Orthoptera: Tettigoniidae) can attract male cicadas (Hemiptera: Cicadidae) by imitating the species-specific wing-flick replies of sexually receptive female cicadas. This aggressive mimicry is accomplished both acoustically, with tegminal clicks, and visually, with synchronized body jerks. Remarkably, the katydids respond effectively to a variety of complex, species-specific Cicadettini songs, including songs of many cicada species that the predator has never encountered. CONCLUSIONS/SIGNIFICANCE: We propose that the versatility of aggressive mimicry in C. leucoviridis is accomplished by exploiting general design elements common to the songs of many acoustically signaling insects that use duets in pair-formation. Consideration of the mechanism of versatile mimicry in C. leucoviridis may illuminate processes driving the evolution of insect acoustic signals, which play a central role in reproductive isolation of populations and the formation of species.

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.