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.

Publisher Correction: The tuatara genome reveals ancient features of amniote evolution.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Multiple hybridization events and repeated evolution of homoeologue expression bias in parthenogenetic, polyploid New Zealand stick insects.

During hybrid speciation, homoeologues combine in a single genome. Homoeologue expression bias (HEB) occurs when one homoeologue has higher gene expression than another. HEB has been well characterized in plants but rarely investigated in animals, especially invertebrates. Consequently, we have little idea as to the role that HEB plays in allopolyploid invertebrate genomes. If HEB is constrained by features of the parental genomes, then we predict repeated evolution of similar HEB patterns among hybrid genomes formed from the same parental lineages. To address this, we reconstructed the history of hybridization between the New Zealand stick insect genera Acanthoxyla and Clitarchus using a high-quality genome assembly from Clitarchus hookeri to call variants and phase alleles. These analyses revealed the formation of three independent diploid and triploid hybrid lineages between these genera. RNA sequencing revealed a similar magnitude and direction of HEB among these hybrid lineages, and we observed that many enriched functions and pathways were also shared among lineages, consistent with repeated evolution due to parental genome constraints. In most hybrid lineages, a slight majority of the genes involved in mitochondrial function showed HEB towards the maternal homoeologues, consistent with only weak effects of mitonuclear incompatibility. We also observed a proteasome functional enrichment in most lineages and hypothesize this may result from the need to maintain proteostasis in hybrid genomes. Reference bias was a pervasive problem, and we caution against relying on HEB estimates from a single parental reference genome.

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.

Weak premating isolation between Clitarchus stick insect species despite divergent male and female genital morphology.

Documenting natural hybrid systems builds our understanding of mate choice, reproductive isolation and speciation. The stick insect species Clitarchus hookeri and C. tepaki differ in their genital morphology and hybridize along a narrow peninsula in northern New Zealand. We utilize three lines of evidence to understand the role of premating isolation and species boundaries: (a) genetic differentiation using microsatellites and mitochondrial DNA; (b) variation in 3D surface topology of male claspers and 2D morphometrics of female opercular organs; and (c) behavioural reproductive isolation among parental and hybrid populations through mating crosses. The genetic data show introgression between the parental species and formation of a genetically variable hybrid swarm. Similarly, the male and female morphometric data show genital divergence between the parental species as well as increased variation within the hybrid populations. This genital divergence has not resulted in reproductive isolation between species, instead weak perimating isolation has enabled the formation of a hybrid swarm. Behavioural analysis demonstrates that the entire mating process influences the degree of reproductive isolation between species undergoing secondary contact. Mechanical isolation may appear strong, whereas perimating isolation is weak.

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.

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.

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.

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.

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.

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.

Divergent transcriptional responses to low temperature among populations of alpine and lowland species of New Zealand stick insects (Micrarchus).

In widespread and genetically structured populations, temperature variation may lead to among-population differentiation of thermal biology. The New Zealand stick insect genus Micrarchus contains four species that inhabit different thermal environments, two of which are geographically widespread. RNA-Seq and quantitative PCR were used to investigate the transcriptional responses to cold shock among lowland and alpine species to identify cold-responsive transcripts that differ between the species and to determine whether there is intraspecific geographical variation in gene expression. We also used mitochondrial DNA, nuclear 28S ribosomal DNA and transcriptome-wide SNPs to determine phylogeographic structure and the potential for differences in genetic backgrounds to contribute to variation in gene expression. RNA-Seq identified 2160 unigenes that were differentially expressed as a result of low-temperature exposure across three populations from two species (M. hystriculeus and M. nov. sp. 2), with a majority (68% ± 20%) being population specific. This extensive geographical variation is consistent across years and is likely a result of background genetic differences among populations caused by genetic drift and possibly local adaptation. Responses to cold shock shared among alpine M. nov. sp. 2 populations included the enrichment of cuticular structure-associated transcripts, suggesting that cuticle modification may have accompanied colonization of low-temperature alpine environments and the development of a more cold-hardy phenotype.

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.

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.

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.

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.

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.

Reconciling phylogeography and ecological niche models for New Zealand beetles: Looking beyond glacial refugia.

Mitochondrial DNA (cox1) sequence data and recently developed coalescent phylogeography models were used to construct geo-spatial histories for the New Zealand fungus beetles Epistranus lawsoni and Pristoderus bakewelli (Zopheridae). These methods utilize continuous-time Markov chains and Bayesian stochastic search variable selection incorporated in BEAST to identify historical dispersal patterns via ancestral state reconstruction. Ecological niche models (ENMs) were incorporated to reconstruct the potential geographic distribution of each species during the Last Glacial Maximum (LGM). Coalescent analyses suggest a North Island origin for E. lawsoni, with gene flow predominately north-south between adjacent regions. ENMs for E. lawsoni indicated glacial refugia in coastal regions of both main islands, consistent with phylogenetic patterns but at odds with the coalescent dates, which implicate much older topographic events. Dispersal matrices revealed patterns of gene flow consistent with projected refugia, suggesting long-term South Island survival with population vicariance around the Southern Alps. Phylogeographic relationships are more ambiguous for P. bakewelli, although long-term survival on both main islands is evident. Divergence dates for both species are consistent with the topographic evolution of New Zealand over the last 10Ma, whereas the signature of the LGM is less apparent in the time-scaled phylogeny.

Phylogenetic analysis of New Zealand earthworms (Oligochaeta: Megascolecidae) reveals ancient clades and cryptic taxonomic diversity.

We have constructed the first ever phylogeny for the New Zealand earthworm fauna (Megascolecinae and Acanthodrilinae) including representatives from other major continental regions. Bayesian and maximum likelihood phylogenetic trees were constructed from 427 base pairs from the mitochondrial large subunit (16S) rRNA gene and 661 base pairs from the nuclear large subunit (28S) rRNA gene. Within the Acanthodrilinae we were able to identify a number of well-supported clades that were restricted to continental landmasses. Estimates of nodal support for these major clades were generally high, but relationships among clades were poorly resolved. The phylogenetic analyses revealed several independent lineages in New Zealand, some of which had a comparable phylogenetic depth to monophyletic groups sampled from Madagascar, Africa, North America and Australia. These results are consistent with at least some of these clades having inhabited New Zealand since rifting from Gondwana in the Late Cretaceous. Within the New Zealand Acanthodrilinae, major clades tended to be restricted to specific regions of New Zealand, with the central North Island and Cook Strait representing major biogeographic boundaries. Our field surveys of New Zealand and subsequent identification has also revealed extensive cryptic taxonomic diversity with approximately 48 new species sampled in addition to the 199 species recognized by previous authors. Our results indicate that further survey and taxonomic work is required to establish a foundation for future biogeographic and ecological research on this vitally important component of the New Zealand biota.

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.