Wing pattern variation and DNA barcodes defy taxonomic splitting in the New Zealand Pimelea Looper Notoreas perornata (Walker) (Lepidoptera: Geometridae: Larentiinae): the importance of populations as conservation units

The endemic Notoreas perornata (Walker, 1863) complex (Lepidoptera: Geometridae: Larentiinae) from the North Island and northern South Island of New Zealand is reviewed. Larvae feed on Pimelea spp. (Thymelaeaceae), frequently in highly fragmented and threatened shrubland habitats. Allopatric populations tend to differ in size and wing pattern characteristics, but not in genitalia; moreover extensive variation renders recognition of subspecies / allopatric species based on any species concept problematic. A mitochondrial DNA gene tree is not congruent with morphology and indicates rapid recent divergence that has not settled into diagnosable lineages. Based on our results, we synonymise Notoreas simplex Hudson, 1898 with N. perornata (Walker, 1863), and retain N. perornata as a single, highly diverse but monotypic species. All known populations are illustrated to display variation. For conservation purposes, we recommend the continued recognition within the species of 10 populations or groups of populations that appear to be on the way to diverging at subspecific level based on morphological and/or DNA data. The conservation status of all these populations is reviewed. One conservation unit, comprising the populations from Westland, has not been seen since 1998 and is feared possibly extinct.

Fast-tracking bespoke DNA reference database generation from museum collections for biomonitoring and conservation

Despite recent advances in high-throughput DNA sequencing technologies, a lack of locally relevant DNA reference databases limits the potential for DNA-based monitoring of biodiversity for conservation and biosecurity applications. Museums and national collections represent a compelling source of authoritatively identified genetic material for DNA database development, yet obtaining DNA barcodes from long-stored specimens may be difficult due to sample degradation. Here we demonstrate a sensitive and efficient laboratory and bioinformatic process for generating DNA barcodes from hundreds of invertebrate specimens simultaneously via the Illumina MiSeq system. Using this process, we recovered full-length (334) or partial (105) COI barcodes from 439 of 450 (98%) national collection-held invertebrate specimens. This included full-length barcodes from 146 specimens which produced low-yield DNA and no visible PCR bands, and which produced as little as a single sequence per specimen, demonstrating high sensitivity of the process. In many cases, the identity of the most abundant sequences per specimen were not the correct barcodes, necessitating the development of a taxonomy-informed process for identifying correct sequences among the sequencing output. The recovery of only partial barcodes for some taxa indicates a need to refine certain PCR primers. Nonetheless, our approach represents a highly sensitive, accurate and efficient method for targeted reference database generation, providing a foundation for DNA-based assessments and monitoring of biodiversity.

An exploration of assembly strategies and quality metrics on the accuracy of the rewarewa (Knightia excelsa) genome

We used long read sequencing data generated from Knightia excelsa, a nectar-producing Proteaceae tree endemic to Aotearoa (New Zealand), to explore how sequencing data type, volume and workflows can impact final assembly accuracy and chromosome reconstruction. Establishing a high-quality genome for this species has specific cultural importance to Māori and commercial importance to honey producers in Aotearoa. Assemblies were produced by five long read assemblers using data subsampled based on read lengths, two polishing strategies and two Hi-C mapping methods. Our results from subsampling the data by read length showed that each assembler tested performed differently depending on the coverage and the read length of the data. Subsampling highlighted that input data with longer read lengths but perhaps lower coverage constructed more contiguous, kmers and gene-complete assemblies than short read length input data with higher coverage. The final genome assembly was constructed into 14 pseudochromosomes using an initial flye long read assembly, a racon/medaka/pilon combined polishing strategy, salsa2 and allhic scaffolding, juicebox curation, and Macadamia linkage map validation. We highlighted the importance of developing assembly workflows based on the volume and read length of sequencing data and established a robust set of quality metrics for generating high-quality assemblies. Scaffolding analyses highlighted that problems found in the initial assemblies could not be resolved accurately by Hi-C data and that assembly scaffolding was more successful when the underlying contig assembly was of higher accuracy. These findings provide insight into how quality assessment tools can be implemented throughout genome assembly pipelines to inform the de novo reconstruction of a high-quality genome assembly for nonmodel organisms.

Divergent Gene Expression Following Duplication of Meiotic Genes in the Stick Insect Clitarchus hookeri

Some animal groups, such as stick insects (Phasmatodea), have repeatedly evolved alternative reproductive strategies, including parthenogenesis. Genomic studies have found modification of the genes underlying meiosis exists in some of these animals. Here we examine the evolution of copy number, evolutionary rate, and gene expression in candidate meiotic genes of the New Zealand geographic parthenogenetic stick insect Clitarchus hookeri. We characterized 101 genes from a de novo transcriptome assembly from female and male gonads that have homology with meiotic genes from other arthropods. For each gene we determined copy number, the pattern of gene duplication relative to other arthropod orthologs, and the potential for meiosis-specific expression. There are five genes duplicated in C. hookeri, including one also duplicated in the stick insect Timema cristinae, that are not or are uncommonly duplicated in other arthropods. These included two sister chromatid cohesion associated genes (SA2 and SCC2), a recombination gene (HOP1), an RNA-silencing gene (AGO2) and a cell-cycle regulation gene (WEE1). Interestingly, WEE1 and SA2 are also duplicated in the cyclical parthenogenetic aphid Acyrthosiphon pisum and Daphnia duplex, respectively, indicating possible roles in the evolution of reproductive mode. Three of these genes (SA2, SCC2, and WEE1) have one copy displaying gonad-specific expression. All genes, with the exception of WEE1, have significantly different nonsynonymous/synonymous ratios between the gene duplicates, indicative of a shift in evolutionary constraints following duplication. These results suggest that stick insects may have evolved genes with novel functions in gamete production by gene duplication.

The tuatara genome reveals ancient features of amniote evolution

The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2-is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which-at approximately 5 Gb-is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing.

A new leaf-mining moth from New Zealand, Sabulopteryxbotanica sp. nov. (Lepidoptera, Gracillariidae, Gracillariinae), feeding on the rare endemic shrub Teucriumparvifolium (Lamiaceae), with a revised checklist of New Zealand Gracillariidae

Sabulopteryxbotanica Hoare & Patrick, sp. nov. (Lepidoptera, Gracillariidae, Gracillariinae) is described as a new species from New Zealand. It is regarded as endemic, and represents the first record of its genus from the southern hemisphere. Though diverging in some morphological features from previously described species, it is placed in genus Sabulopteryx Triberti, based on wing venation, abdominal characters, male and female genitalia and hostplant choice; this placement is supported by phylogenetic analysis based on the COI mitochondrial gene. The life history is described: the larva is an underside leaf-miner on the endemic divaricating shrub Teucriumparvifolium (Lamiaceae), and exits the mine to pupate in a cocoon in a folded leaf of the host plant. The remarkable history of the discovery and rediscovery of this moth is discussed: for many years it was only known from a single sap-feeding larva found in a leaf-mine in a pressed herbarium specimen of the host. The adult was discovered by BHP in Christchurch Botanic Gardens in 2013. Most distribution records of the moth come from a recent search for mines and cocoons on herbarium specimens of T.parvifolium. Sabulopteryxbotanica has high conservation status, and is regarded as 'Nationally Vulnerable' according to the New Zealand Department of Conservation threat classification system, based on the rarity and declining status of its host plant. However, the presence of apparently thriving populations of S.botanica on cultivated plants of T.parvifolium, especially at the type locality, Christchurch Botanic Gardens, suggests that encouraging cultivation of the plant could greatly improve the conservation status of the moth. A revised checklist of New Zealand Gracillariidae is presented, assigning all species to the currently recognised subfamilies. The Australian Macarostolaida (Meyrick, 1880) is newly recorded from New Zealand (Auckland), where it is established on Eucalyptus.

Estimating the biodiversity of terrestrial invertebrates on a forested island using DNA barcodes and metabarcoding data

Invertebrates are a major component of terrestrial ecosystems, however, estimating their biodiversity is challenging. We compiled an inventory of invertebrate biodiversity along an elevation gradient on the temperate forested island of Hauturu, New Zealand, by DNA barcoding of specimens obtained from leaf litter samples and pitfall traps. We compared the barcodes and biodiversity estimates from this data set with those from a parallel DNA metabarcoding analysis of soil from the same locations, and with pre-existing sequences in reference databases, before exploring the use of combined data sets as a basis for estimating total invertebrate biodiversity. We obtained 1,282 28S and 1,610 COI barcodes from a total of 1,947 invertebrate specimens, which were clustered into 247 (28S) and 366 (COI) OTUs, of which ≤ 10% were represented in GenBank. Coleoptera were most abundant (730 sequenced specimens), followed by Hymenoptera, Diptera, Lepidoptera, and Amphipoda. The most abundant OTU from both the 28S (153 sequences) and COI (140 sequences) data sets was an undescribed beetle from the family Salpingidae. Based on the occurrences of COI OTUs along the elevation gradient, we estimated there are ~1,000 arthropod species (excluding mites) on Hauturu, including 770 insects, of which 344 are beetles. A DNA metabarcoding analysis of soil DNA from the same sites resulted in the identification of similar numbers of OTUs in most invertebrate groups compared with the DNA barcoding, but less than 10% of the DNA barcoding COI OTUs were also detected by the metabarcoding analysis of soil DNA. A mark-recapture analysis based on the overlap between these data sets estimated the presence of approximately 6,800 arthropod species (excluding mites) on the island, including ~3,900 insects. Estimates of New Zealand-wide biodiversity for selected arthropod groups based on matching of the COI DNA barcodes with pre-existing reference sequences suggested over 13,200 insect species are present, including 4,000 Coleoptera, 2,200 Diptera, and 2,700 Hymenoptera species, and 1,000 arachnid species (excluding mites). These results confirm that metabarcoding analyses of soil DNA tends to recover different components of terrestrial invertebrate biodiversity compared to traditional invertebrate sampling, but the combined methods provide a novel basis for estimating invertebrate biodiversity.

The genus Mecodema Blanchard 1853 (Coleoptera: Carabidae: Broscini) from the North Island, New Zealand

Mecodema (Coleoptera: Carabidae: Broscini) is a hyperdiverse endemic New Zealand genus of ground beetles with only a few geographically widespread species found throughout the two main islands, as well as many offshore islands. Using specimens from a number of private and institutional collections, in addition new specimens were acquired by extensive pitfall trapping, we describe or redescribe all of the known North Island Mecodema species. Additionally, we redescribe three South Island species from the former genus Metaglymma, as morphological evidence shows that these species are nested within Mecodema. Species descriptions are formed by using 128 morphological characters, which include external characters, as well as both male and female internal structures. There are four new combinations: Mecodema antarctica comb. n., M. aberrans comb. n., M. moniliferum comb. n. and M. tibiale comb. n. We synonymise M. occiputale under Mecodema curvidens, and M. sulcatum under Mecodema oblongum, and reinstate M. scitulum Broun (northwest Hunua Range, Auckland). Twenty four new species are described: Mecodema argentum sp. n., M. atuanui sp. n., M. dunnorum sp. n., M. genesispotini sp. n., M. godzilla sp. n., M. jacinda sp. n., M. kipjac sp. n., M. kokoroiho sp. n., M. mohi sp. n., M. ngaiatonga sp. n., M. ngaitahuhu sp. n., M. papake sp. n., M. perexiguus sp. n., M. rusticulus sp. n., M. temata sp. n., M. teparawhau sp. n., M. teroroa sp. n., M. tewhara sp. n., M. tuhoe sp. n., M. undecimus sp. n., M. wharekahika sp. n., M. xylanthrax sp. n., M. yconomus sp. n., M. zonula sp. n. North Island regional species endemism is very high in Northland (15/16 endemic species), with species becoming more widespread in the southern regions, e.g. Wellington only has two endemic species from a total of eight species. This research increases the total number of described Mecodema species to 102, and will allow a modern taxonomic framework for completion of the revision of the South Island species.

New Zealand Tree and Giant Wētā (Orthoptera) Transcriptomics Reveal Divergent Selection Patterns in Metabolic Loci

Exposure to low temperatures requires an organism to overcome physiological challenges. New Zealand wētā belonging to the genera Hemideina and Deinacrida are found across a wide range of thermal environments and therefore subject to varying selective pressures. Here we assess the selection pressures across the wētā phylogeny, with a particular emphasis on identifying genes under positive or diversifying selection. We used RNA-seq to generate transcriptomes for all 18 Deinacrida and Hemideina species. A total of 755 orthologous genes were identified using a bidirectional best-hit approach, with the resulting gene set encompassing a diverse range of functional classes. Analysis of ortholog ratios of synonymous to nonsynonymous amino acid changes found 83 genes that are under positive selection for at least one codon. A wide variety of Gene Ontology terms, enzymes, and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways are represented among these genes. In particular, enzymes involved in oxidative phosphorylation, melanin synthesis, and free-radical scavenging are represented, consistent with physiological and metabolic changes that are associated with adaptation to alpine environments. Structural alignment of the transcripts with the most codons under positive selection revealed that the majority of sites are surface residues, and therefore have the potential to influence the thermostability of the enzyme, with the exception of prophenoloxidase where two residues near the active site are under selection. These proteins provide interesting candidates for further analysis of protein evolution.

Evolutionary history of Polyneoptera and its implications for our understanding of early winged insects

Polyneoptera represents one of the major lineages of winged insects, comprising around 40,000 extant species in 10 traditional orders, including grasshoppers, roaches, and stoneflies. Many important aspects of polyneopteran evolution, such as their phylogenetic relationships, changes in their external appearance, their habitat preferences, and social behavior, are unresolved and are a major enigma in entomology. These ambiguities also have direct consequences for our understanding of the evolution of winged insects in general; for example, with respect to the ancestral habitats of adults and juveniles. We addressed these issues with a large-scale phylogenomic analysis and used the reconstructed phylogenetic relationships to trace the evolution of 112 characters associated with the external appearance and the lifestyle of winged insects. Our inferences suggest that the last common ancestors of Polyneoptera and of the winged insects were terrestrial throughout their lives, implying that wings did not evolve in an aquatic environment. The appearance of the first polyneopteran insect was mainly characterized by ancestral traits such as long segmented abdominal appendages and biting mouthparts held below the head capsule. This ancestor lived in association with the ground, which led to various specializations including hardened forewings and unique tarsal attachment structures. However, within Polyneoptera, several groups switched separately to a life on plants. In contrast to a previous hypothesis, we found that social behavior was not part of the polyneopteran ground plan. In other traits, such as the biting mouthparts, Polyneoptera shows a high degree of evolutionary conservatism unique among the major lineages of winged insects.

Systematics of the New Zealand Weevil Etheophanus Broun (Curculionidae: Molytinae)

Etheophanus Broun is considered a molytine based on the form of the pharyngeal plate, presence of a small spiculum relictum in the male, and presence of a pair of small internal apodemes on the antero-lateral corners of the 5th abdominal ventrite of the female. Examination of primary type specimens and newer material confirm one new species Etheophanus kuscheli sp. n. and two synonomies (Etheophanus nitidellus Broun, 1923 [= Etheophanus obscurus Broun, 1923] and Etheophanus striatus Broun, 1910 [=Etheophanus punctiventris Broun, 1914]). Generic and species diagnoses, a key to the species, and lectotype designations for three species are included. Phylogenetic reconstructions based on a combined analysis of the nuclear 28S rRNA and mitochondrial cytochrome c oxidase subunit I genes confirmed the status of E. kuscheli and a species complex, the E. nitidellus/E. optandus clade distributed in the southern portion of the South Island. The relationship E. pinguis [northern North Island] (E. striatus [southern North Island, northern South Island] (E. kuscheli [northwestern South Island] (E. nitidellus, E. optandus [southwestern North Island]) corresponds to geographic patterns found in other beetle lineages. Etheophanus striatus is composed of three lineages, one widespread in the north and south islands and two allopatric populations in the northwest South Island. The E. nitidellus/E. optandus complex includes four distinct lineages, one restricted to Fiordland, the other three sympatric in the region affected by the Haast Corridor.

Assembling large genomes: analysis of the stick insect (Clitarchus hookeri) genome reveals a high repeat content and sex-biased genes associated with reproduction

BACKGROUND

Stick insects (Phasmatodea) have a high incidence of parthenogenesis and other alternative reproductive strategies, yet the genetic basis of reproduction is poorly understood. Phasmatodea includes nearly 3000 species, yet only the genome of Timema cristinae has been published to date. Clitarchus hookeri is a geographical parthenogenetic stick insect distributed across New Zealand. Sexual reproduction dominates in northern habitats but is replaced by parthenogenesis in the south. Here, we present a de novo genome assembly of a female C. hookeri and use it to detect candidate genes associated with gamete production and development in females and males. We also explore the factors underlying large genome size in stick insects.

RESULTS

The C. hookeri genome assembly was 4.2 Gb, similar to the flow cytometry estimate, making it the second largest insect genome sequenced and assembled to date. Like the large genome of Locusta migratoria, the genome of C. hookeri is also highly repetitive and the predicted gene models are much longer than those from most other sequenced insect genomes, largely due to longer introns. Miniature inverted repeat transposable elements (MITEs), absent in the much smaller T. cristinae genome, is the most abundant repeat type in the C. hookeri genome assembly. Mapping RNA-Seq reads from female and male gonadal transcriptomes onto the genome assembly resulted in the identification of 39,940 gene loci, 15.8% and 37.6% of which showed female-biased and male-biased expression, respectively. The genes that were over-expressed in females were mostly associated with molecular transportation, developmental process, oocyte growth and reproductive process; whereas, the male-biased genes were enriched in rhythmic process, molecular transducer activity and synapse. Several genes involved in the juvenile hormone synthesis pathway were also identified.

CONCLUSIONS

The evolution of large insect genomes such as L. migratoria and C. hookeri genomes is most likely due to the accumulation of repetitive regions and intron elongation. MITEs contributed significantly to the growth of C. hookeri genome size yet are surprisingly absent from the T. cristinae genome. Sex-biased genes identified from gonadal tissues, including genes involved in juvenile hormone synthesis, provide interesting candidates for the further study of flexible reproduction in stick insects.

Positive selection and comparative molecular evolution of reproductive proteins from New Zealand tree weta (Orthoptera, Hemideina)

Animal reproductive proteins, especially those in the seminal fluid, have been shown to have higher levels of divergence than non-reproductive proteins and are often evolving adaptively. Seminal fluid proteins have been implicated in the formation of reproductive barriers between diverging lineages, and hence represent interesting candidates underlying speciation. RNA-seq was used to generate the first male reproductive transcriptome for the New Zealand tree weta species Hemideina thoracica and H. crassidens. We identified 865 putative reproductive associated proteins across both species, encompassing a diverse range of functional classes. Candidate gene sequencing of nine genes across three Hemideina, and two Deinacrida species suggests that H. thoracica has the highest levels of intraspecific genetic diversity. Non-monophyly was observed in the majority of sequenced genes indicating that either gene flow may be occurring between the species, or that reciprocal monophyly at these loci has yet to be attained. Evidence for positive selection was found for one lectin-related reproductive protein, with an overall omega of 7.65 and one site in particular being under strong positive selection. This candidate gene represents the first step in the identification of proteins underlying the evolutionary basis of weta reproduction and speciation.

Analysis of the genome of the New Zealand giant collembolan (Holacanthella duospinosa) sheds light on hexapod evolution

BACKGROUND

The New Zealand collembolan genus Holacanthella contains the largest species of springtails (Collembola) in the world. Using Illumina technology we have sequenced and assembled a draft genome and transcriptome from Holacanthella duospinosa (Salmon). We have used this annotated assembly to investigate the genetic basis of a range of traits critical to the evolution of the Hexapoda, the phylogenetic position of H. duospinosa and potential horizontal gene transfer events.

RESULTS

Our genome assembly was ~375 Mbp in size with a scaffold N50 of ~230 Kbp and sequencing coverage of ~180×. DNA elements, LTRs and simple repeats and LINEs formed the largest components and SINEs were very rare. Phylogenomics (370,877 amino acids) placed H. duospinosa within the Neanuridae. We recovered orthologs of the conserved sex determination genes thought to play a role in sex determination. Analysis of CpG content suggested the absence of DNA methylation, and consistent with this we were unable to detect orthologs of the DNA methyltransferase enzymes. The small subunit rRNA gene contained a possible retrotransposon. The Hox gene complex was broken over two scaffolds. For chemosensory ability, at least 15 and 18 ionotropic glutamate and gustatory receptors were identified, respectively. However, we were unable to identify any odorant receptors or their obligate co-receptor Orco. Twenty-three chitinase-like genes were identified from the assembly. Members of this multigene family may play roles in the digestion of fungal cell walls, a common food source for these saproxylic organisms. We also detected 59 and 96 genes that blasted to bacteria and fungi, respectively, but were located on scaffolds that otherwise contained arthropod genes.

CONCLUSIONS

The genome of H. duospinosa contains some unusual features including a Hox complex broken over two scaffolds, in a different manner to other arthropod species, a lack of odorant receptor genes and an apparent lack of environmentally responsive DNA methylation, unlike many other arthropods. Our detection of candidate horizontal gene transfer candidates confirms that this phenomenon is occurring across Collembola. These findings allow us to narrow down the regions of the arthropod phylogeny where key innovations have occurred that have facilitated the evolutionary success of Hexapoda.

Applications of phylogenetics to solve practical problems in insect conservation

Phylogenetic approaches have much promise for the setting of conservation priorities and resource allocation. There has been significant development of analytical methods for the measurement of phylogenetic diversity within and among ecological communities as a way of setting conservation priorities. Application of these tools to insects has been low as has been the uptake by conservation managers. A critical reason for the lack of uptake includes the scarcity of detailed phylogenetic and species distribution data from much of insect diversity. Environmental DNA technologies offer a means for the high throughout collection of phylogenetic data across landscapes for conservation planning.

Building strong relationships between conservation genetics and primary industry leads to mutually beneficial genomic advances

Several reviews in the past decade have heralded the benefits of embracing high-throughput sequencing technologies to inform conservation policy and the management of threatened species, but few have offered practical advice on how to expedite the transition from conservation genetics to conservation genomics. Here, we argue that an effective and efficient way to navigate this transition is to capitalize on emerging synergies between conservation genetics and primary industry (e.g., agriculture, fisheries, forestry and horticulture). Here, we demonstrate how building strong relationships between conservation geneticists and primary industry scientists is leading to mutually-beneficial outcomes for both disciplines. Based on our collective experience as collaborative New Zealand-based scientists, we also provide insight for forging these cross-sector relationships.

De Novo Transcriptome Analysis of the Common New Zealand Stick Insect Clitarchus hookeri (Phasmatodea) Reveals Genes Involved in Olfaction, Digestion and Sexual Reproduction

Phasmatodea, more commonly known as stick insects, have been poorly studied at the molecular level for several key traits, such as components of the sensory system and regulators of reproduction and development, impeding a deeper understanding of their functional biology. Here, we employ de novo transcriptome analysis to identify genes with primary functions related to female odour reception, digestion, and male sexual traits in the New Zealand common stick insect Clitarchus hookeri (White). The female olfactory gene repertoire revealed ten odorant binding proteins with three recently duplicated, 12 chemosensory proteins, 16 odorant receptors, and 17 ionotropic receptors. The majority of these olfactory genes were over-expressed in female antennae and have the inferred function of odorant reception. Others that were predominantly expressed in male terminalia (n = 3) and female midgut (n = 1) suggest they have a role in sexual reproduction and digestion, respectively. Over-represented transcripts in the midgut were enriched with digestive enzyme gene families. Clitarchus hookeri is likely to harbour nine members of an endogenous cellulase family (glycoside hydrolase family 9), two of which appear to be specific to the C. hookeri lineage. All of these cellulase sequences fall into four main phasmid clades and show gene duplication events occurred early in the diversification of Phasmatodea. In addition, C. hookeri genome is likely to express γ-proteobacteria pectinase transcripts that have recently been shown to be the result of horizontal transfer. We also predicted 711 male terminalia-enriched transcripts that are candidate accessory gland proteins, 28 of which were annotated to have molecular functions of peptidase activity and peptidase inhibitor activity, two groups being widely reported to regulate female reproduction through proteolytic cascades. Our study has yielded new insights into the genetic basis of odour detection, nutrient digestion, and male sexual traits in stick insects. The C. hookeri reference transcriptome, together with identified gene families, provides a comprehensive resource for studying the evolution of sensory perception, digestive systems, and reproductive success in phasmids.

Evaluating a multigene environmental DNA approach for biodiversity assessment

Background

There is an increasing demand for rapid biodiversity assessment tools that have a broad taxonomic coverage. Here we evaluate a suite of environmental DNA (eDNA) markers coupled with next generation sequencing (NGS) that span the tree of life, comparing them with traditional biodiversity monitoring tools within ten 20×20 meter plots along a 700 meter elevational gradient.

Results

From six eDNA datasets (one from each of 16S, 18S, ITS, trnL and two from COI) we identified sequences from 109 NCBI taxonomy-defined phyla or equivalent, ranging from 31 to 60 for a given eDNA marker. Estimates of alpha and gamma diversity were sensitive to the number of sequence reads, whereas beta diversity estimates were less sensitive. The average within-plot beta diversity was lower than between plots for all markers. The soil beta diversity of COI and 18S markers showed the strongest response to the elevational variation of the eDNA markers (COI: r=0.49, p<0.001; 18S: r=0.48, p<0.001). Furthermore pairwise beta diversities for these two markers were strongly correlated with those calculated from traditional vegetation and invertebrate biodiversity measures.

Conclusions

Using a soil-based eDNA approach, we demonstrate that standard phylogenetic markers are capable of recovering sequences from a broad diversity of eukaryotes, in addition to prokaryotes by 16S. The COI and 18S eDNA markers are the best proxies for aboveground biodiversity based on the high correlation between the pairwise beta diversities of these markers and those obtained using traditional methods.

Electronic supplementary material

The online version of this article (doi:10.1186/s13742-015-0086-1) contains supplementary material, which is available to authorized users.

Single origin of the Mascarene stick insects: ancient radiation on sunken islands?

BACKGROUND

The study of islands as model systems plays a key role in understanding many evolutionary processes. Knowledge of the historical events leading to present-day island communities is pivotal for exploring fundamental mechanisms of speciation and adaptation. The remote Mascarene archipelago (Mauritius, Réunion, Rodrigues), considered to be the product of an age-progressive trend of north-to-south volcanic activity in the Indian Ocean, hosts a remarkably diverse, endemic and threatened concentration of flora and fauna that has traditionally been considered to be biogeographically related to Madagascar and Africa. To explore the evolutionary diversity of the Mascarene stick insects (Phasmatodea), we constructed a global phylogeny from approximately 2.4 kb of mitochondrial and nuclear sequence data of more than 120 species representing all major phasmatodean lineages.

RESULTS

Based on the obtained time-calibrated molecular tree we demonstrate that the current phasmid community of the Mascarene archipelago, which consists of members of four presumably unrelated traditional subfamilies, is the result of a single ancient dispersal event from Australasia and started radiating between 16-29 million years ago, significantly predating the age of Mauritius (8-10 million years).

CONCLUSIONS

We propose that the Mascarene stick insects diversified on landmasses now eroded away, presumably to the north of Mauritius. In consequence, ancient islands have probably persisted in the Indian Ocean until the emergence of Mauritius and not only served as stepping stones for colonisation events during sea-level lowstands, but as long-lasting cradles of evolution. These ancient landmasses most likely allowed for adaptive speciation and served as significant sources of diversity that contributed to the biomes of the Mascarene archipelago and the megadiverse Madagascar.

Mate detection and seasonal variation in stick insect mating behaviour (Phamatodea: Clitarchus hookeri)

For animals that exhibit a scramble competition mating system, sexual selection pressures on mate searching ability are expected to be strong. Scramble competition mating systems evolve when populations provide females with equal accessibility to all male competitors, yet sex ratio and population density influences mating systems and varies seasonally. The stick insect species,Clitarchus hookeri, is frequently found in copula, yet very little is known about it’s mating behaviour. We preformed behavioural tests and assayed antennal sensory morphology to determine whether males used chemosensory cues to detect females. Through natural field observations we found populations to be significantly male-biased earlier in the season, while later, populations began to display equal sex ratios. With increasing female availability mating pair proportions steadily increased, while copulation duration declined. These results supportC. hookerias a scramble competitor, and demonstrate males may alter their behaviour in response to the seasonal variation in female density.

Revision of the stick insect genus Clitarchus Stål (Phasmatodea: Phasmatidae): new synonymies and two new species from northern New Zealand

We describe two new species of Clitarchus Stål from Northland, New Zealand. Clitarchus rakauwhakanekeneke sp. nov. is described from the Poor Knights Islands and Clitarchus tepaki sp. nov. is described from the Te Paki / North Cape area and the Karikari Peninsula at the northernmost tip of New Zealand. Two new synonymies are proposed including Clitarchus multidentatus Brunner (syn. nov.) and Clitarchus tuberculatus Salmon (syn. nov.) as synonyms of Clitarchus hookeri (White). Clitarchus magnus Brunner, recorded from Thailand, is transferred to Ramulus Saussure and given the replacement name Ramulus changmaiense nom. nov. The holotype of C. multidentatus was recorded as being collected from New Caledonia; however we believe this is a labelling error and the specimen was collected from New Zealand. These taxonomic changes render Clitarchus endemic to New Zealand and consisting of three species; C. hookeri, C. rakauwhakanekeneke and C. tepaki. Keys to the adult males and females of Clitarchus species are given in addition to notes on host plants, ecology and geographic distributions.

A living fossil tale of Pangaean biogeography

The current distributions of widespread groups of terrestrial animals and plants are supposedly the result of a mixture of either vicariance owing to continental split or more recent trans-oceanic dispersal. For organisms exhibiting a vicariant biogeographic pattern-achieving their current distribution by riding on the plates of former supercontinents-this view is largely inspired by the belief that Pangaea lacked geographical or ecological barriers, or that extinctions and dispersal would have erased any biogeographic signal since the early Mesozoic. We here present a time-calibrated molecular phylogeny of Onychophora (velvet worms), an ancient and exclusively terrestrial panarthropod group distributed throughout former Pangaean landmasses. Our data not only demonstrate that trans-oceanic dispersal does not need be invoked to explain contemporary distributions, but also reveal that the early diversification of the group pre-dates the break-up of Pangaea, maintaining regionalization even in landmasses that have remained contiguous throughout the history of the group. These results corroborate a growing body of evidence from palaeontology, palaeogeography and palaeoclimatic modelling depicting ancient biogeographic regionalization over the continuous landmass of Pangaea.

Inclusion of chloroplast genes that have undergone expansion misleads phylogenetic reconstruction in the Chlorophyta

PREMISE OF THE STUDY

Chlorophytes comprise a substantial proportion of green plant diversity. However, sister-group relationships and circumscription of the classes Chlorophyceae, Trebouxiophyceae, and Ulvophyceae have been problematic to resolve. Some analyses support a sister relationship between the trebouxiophycean Leptosira and chlorophyceans, potentially altering the circumscription of two classes, also supported by a shared fragmentation in the chloroplast gene rpoB. We sought to determine whether the latter is a synapomorphy or whether the supporting analyses are vulnerable to systematic bias.

METHODS

We sequenced a portion of rpoB spanning the fragmented region in strains for which it had not previously been sampled: four Chlorophyceae, six counterclockwise (CCW) group (ulvophyceans and trebouxiophyceans) and one streptophyte. We then explored the effect of subsampling proteins and taxa on phylogenetic reconstruction from a data set of 41 chloroplast proteins.

KEY RESULTS

None of the CCW or streptophyte strains possessed the split in rpoB, including inferred near relatives of Leptosira, but it was found in all chlorophycean strains. We reconstructed alternative phylogenies (Leptosira + Chlorophyceae and Leptosira + Chlorellales) using two different protein groups (Rpo and Rps), both subject to coding-region expansion. A conserved region of RpoB remained suitable for analysis of more recent divergences.

CONCLUSIONS

The Rps sequences can explain earlier findings linking Leptosira with the Chlorophyceae and should be excluded from phylogenetic analyses attempting to resolve deep nodes because their expansion violates the assumptions of substitution models. We reaffirm that Leptosira is a trebouxiophycean and that fragmentation of rpoB has occurred at least twice in chlorophyte evolution.

Positive selection in glycolysis among Australasian stick insects

BACKGROUND

The glycolytic pathway is central to cellular energy production. Selection on individual enzymes within glycolysis, particularly phosphoglucose isomerase (Pgi), has been associated with metabolic performance in numerous organisms. Nonetheless, how whole energy-producing pathways evolve to allow organisms to thrive in different environments and adopt new lifestyles remains little explored. The Lanceocercata radiation of Australasian stick insects includes transitions from tropical to temperate climates, lowland to alpine habitats, and winged to wingless forms. This permits a broad investigation to determine which steps within glycolysis and what sites within enzymes are the targets of positive selection. To address these questions we obtained transcript sequences from seven core glycolysis enzymes, including two Pgi paralogues, from 29 Lanceocercata species.

RESULTS

Using maximum likelihood methods a signature of positive selection was inferred in two core glycolysis enzymes. Pgi and Glyceraldehyde 3-phosphate dehydrogenase (Gaphd) genes both encode enzymes linking glycolysis to the pentose phosphate pathway. Positive selection among Pgi paralogues and orthologues predominately targets amino acids with residues exposed to the protein's surface, where changes in physical properties may alter enzyme performance.

CONCLUSION

Our results suggest that, for Lancerocercata stick insects, adaptation to new stressful lifestyles requires a balance between maintaining cellular energy production, efficiently exploiting different energy storage pools and compensating for stress-induced oxidative damage.

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.

Phylogeographic history of the New Zealand stick insect Niveaphasma annulata (Phasmatodea) estimated from mitochondrial and nuclear loci

We have assessed the utility of a single-copy nuclear locus and mitochondrial DNA (mtDNA) in a phylogeographic study of the New Zealand stick insect Niveaphasma annulata (Hutton). We amplified sequences from the mitochondrial cytochrome oxidase subunit I (COI) gene and the single-copy nuclear gene elongation factor-1alpha (EF1alpha) from 97 individuals. Allelic phase at the EF1alpha locus was determined using Denaturing Gradient Gel Electrophoresis. Phylogenetic analyses showed broad congruence between the geographic distribution of three major COI clades and EF1alpha alleles, which suggested that the phylogenetic patterns reflect population history rather than lineage sorting. However, the geographic boundaries of these clades were not always in exact agreement between the two loci. Our data indicate that Niveaphasma annulata was most likely separated into a number of refugia during Pleistocene glacial advances. Subsequent to glacial retreat these refugial populations have expanded and now form a number of zones of secondary contact. We contrast these patterns with those observed from other New Zealand taxa. Our study offers compelling evidence for the use of nuclear genes alongside mtDNA for future phylogeographic studies.

Extreme convergence in stick insect evolution: phylogenetic placement of the Lord Howe Island tree lobster

The 'tree lobsters' are an enigmatic group of robust, ground-dwelling stick insects (order Phasmatodea) from the subfamily Eurycanthinae, distributed in New Guinea, New Caledonia and associated islands. Its most famous member is the Lord Howe Island stick insect Dryococelus australis (Montrouzier), which was believed to have become extinct but was rediscovered in 2001 and is considered to be one of the rarest insects in the world. To resolve the evolutionary position of Dryococelus, we constructed a phylogeny from approximately 2.4 kb of mitochondrial and nuclear sequence data from representatives of all major phasmatodean lineages. Our data placed Dryococelus and the New Caledonian tree lobsters outside the New Guinean Eurycanthinae as members of an unrelated Australasian stick insect clade, the Lanceocercata. These results suggest a convergent origin of the 'tree lobster' body form. Our reanalysis of tree lobster characters provides additional support for our hypothesis of convergent evolution. We conclude that the phenotypic traits leading to the traditional classification are convergent adaptations to ground-living behaviour. Our molecular dating analyses indicate an ancient divergence (more than 22 Myr ago) between Dryococelus and its Australian relatives. Hence, Dryococelus represents a long-standing separate evolutionary lineage within the stick insects and must be regarded as a key taxon to protect with respect to phasmatodean diversity.

DNA barcoding of the endemic New Zealand leafroller moth genera, Ctenopseustis and Planotortrix

Molecular techniques such as DNA barcoding have become popular in assisting species identification especially for cryptic species complexes. We have analysed data from a 468-bp region of the mitochondrial cytochrome oxidase subunit I (COI) gene from 200 specimens of 12 species of endemic New Zealand leafroller moths (Tortricidae) from the genera Planotortrix and Ctenopseustis to assess whether the DNA barcoding region can distinguish these species. Among the 200 sequences analysed, 72 haplotypes were recovered, with each genus forming a separate major clade. Maximum likelihood phylogenetic methods were used to test whether species fell into reciprocally monophyletic clades. The optimal phylogeny showed that four species within the genus Ctenopseustis (C. obliquana, C. herana, C. filicis and C. fraterna) and three within Planotortrix (P. octo, P. excessana and P. avicenniae) are polyphyletic. Shimodaira-Hasegawa tests rejected a null hypothesis of monophyly for the species C. obliquana, C. herana, P. octo and P. excessana. Comparisons of within and between species levels of sequence divergence for the same set of seven species showed cases where maximum levels of within-species divergence were greater than some levels of between-species divergence. DNA barcoding using this region of the COI gene is able to distinguish the two genera and some species within each genus; however, many species cannot be identified using this method. Finally, we discuss the possible reasons for this polyphyly, including incomplete lineage sorting, introgression, horizontal gene transfer and incorrect taxonomy.

Multistate Characters and Diet Shifts: Evolution of Erotylidae (Coleoptera)

The dominance of angiosperms has played a direct role in the diversification of insects, especially Coleoptera. The shift to angiosperm feeding from other diets is likely to have increased the rate of speciation in Phytophaga. However, Phytophaga is only one of many hyperdiverse lineages of beetles and studies of host-shift proliferation have been somewhat limited to groups that primitively feed on plants. We have studied the diet-diverse beetle family Erotylidae (Cucujoidea) to determine if diet is correlated with high diversification rates and morphological evolution by first reconstructing ancestral diets and then testing for associations between diet and species number and diet and ovipositor type. A Bayesian phylogenetic analysis of morphological data that was previously published in Leschen (2003, Pages 1-108 in Fauna of New Zealand, 47; 53 terminal taxa and 1 outgroup, 120 adult characters and 1 diet character) yielded results that are similar to the parsimony analyses of Leschen (2003). Ancestral state reconstructions based on Bayesian and parsimony inference were largely congruent and both reconstructed microfungal feeding (the diet of the outgroup Biphyllidae) at the root of the Erotylidae tree. Shifts among microfungal, saprophagous, and phytophagous diets were most frequent. The largest numbers of species are contained in lineages that are macrofungal feeders (subfamily Erotylinae) and phytophagous (derived Languriinae), although the Bayesian posterior predictive tests of character state correlation were unable to detect any significant associations. Ovipositor morphology correlated with diet (i.e., acute forms were associated with phytophagy and unspecialized forms were associated with a mixture of diets). Although there is a general trend to increased species number associated with the shift from microfungal feeding to phytophagy (based on character mapping and mainly restricted to shifts in Languriinae), there is a large radiation of taxa feeding on macrofungi. Cycad feeding is scattered in more deeply diverged taxa and may have preceded the evolution of angiosperm feeding in some groups. Preliminary analysis of diet mapped onto higher beetle phylogenies suggests that about half of the major Coleoptera lineages may have had fungus-feeding ancestors. We discuss the roles of stochastic models and prior distributions of the reconstruction of ancestral character states in the context of the current data.

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.