Gorges partition diversity within New Zealand flathead Galaxias populations

Understanding the landscape factors governing population connectivity in riverine ecosystems represents an ongoing challenge for freshwater biologists. We used DNA sequence analysis to test the hypothesis that major geomorphological features underpin freshwater-limited fish diversity in a tectonically dynamic region of New Zealand. Phylogeographic analysis of 101 Galaxias depressiceps cytochrome b sequences, incorporating 55 localities from southern New Zealand, revealed 26 haplotypes, with only one shared among rivers. We detect strong hierarchical genetic differentiation both among and within river systems. Genetic structuring is particularly pronounced across the Taieri River system (63 individuals from 35 sites, 18 haplotypes), with 92% of variation partitioned among locations. Distinctive within-river genetic clusters are invariably associated with major subcatchment units, typically isolated by substantial gorges. The anomalous distribution of a single lineage across a major drainage divide is consistent with local, tectonically driven headwater capture. We conclude that major landscape features such as gorges can strongly partition riverine fish diversity and constrain freshwater biodiversity.

Ancient mitochondrial genomes unveil the origins and evolutionary history of New Zealand's enigmatic takahē and moho

Many avian species endemic to Aotearoa New Zealand were driven to extinction or reduced to relict populations following successive waves of human arrival, due to hunting, habitat destruction and the introduction of mammalian predators. Among the affected species were the large flightless South Island takahē (Porphyrio hochstetteri) and the moho (North Island takahē; P. mantelli), with the latter rendered extinct and the former reduced to a single relictual population. Little is known about the evolutionary history of these species prior to their decline and/or extinction. Here we sequenced mitochondrial genomes from takahē and moho subfossils (12 takahē and 4 moho) and retrieved comparable sequence data from takahē museum skins (n = 5) and contemporary individuals (n = 17) to examine the phylogeny and recent evolutionary history of these species. Our analyses suggest that prehistoric takahē populations lacked deep phylogeographic structure, in contrast to moho, which exhibited significant spatial genetic structure, albeit based on limited sample sizes (n = 4). Temporal genetic comparisons show that takahē have lost much of their mitochondrial genetic diversity, likely due to a sudden demographic decline soon after human arrival (~750 years ago). Time-calibrated phylogenetic analyses strongly support a sister species relationship between takahē and moho, suggesting these flightless taxa diverged around 1.5 million years ago, following a single colonisation of New Zealand by a flighted Porphyrio ancestor approximately 4 million years ago. This study highlights the utility of palaeogenetic approaches for informing the conservation and systematic understanding of endangered species whose ranges have been severely restricted by anthropogenic impacts.

Freshwater eDNA reveals dramatic biological shifts linked to deforestation of New Zealand

Deforestation is considered a major threat to biodiversity across many parts of the globe, but the biological impacts of this dramatic ecosystem disturbance often remain incompletely understood. In New Zealand - the world's last major landmass to be colonised by humans - widespread deforestation over recent centuries has left a highly fragmented suite of relict forest stands, ideal for assessing anthropogenic biological change. We hypothesise that this widespread environmental disturbance has underpinned repeated and predictable ecological shifts across distinct rivers and regions. Here we use freshwater environmental DNA (eDNA) data (113 samples across 38 locations; 89 insect taxa) to test for concordant biological shifts linked to this deforestation. eDNA analyses highlight consistent compositional and functional differentiation between forested versus deforested assemblages, including turnover of 'cryptic' congeneric taxa that are morphologically similar yet ecologically and genetically distinct. These dramatic biological shifts are evident even over fine spatial scales within streams, emphasising the widespread emergence of a novel 'deforested' assemblage. Our results illustrate that environmental change can drive predictable biological shifts across broad geographic regions, and highlight the power of eDNA for assessing anthropogenic ecosystem change over large geographic scales.

Fish biogeography and hybridization: do contemporary distributions predict introgression history?

Freshwater ecosystems frequently house diverse assemblages of closely related fish taxa, which can be particularly prone to hybridization and introgression. While extensive introgression may be expected among biogeographically proximate lineages, recent analyses imply that contemporary distributions do not always accurately predict hybridization history. Here, we use the ABBA-BABA approach to test biogeographic hypotheses regarding the extent of hybridization in the recent evolution of New Zealand's species-rich freshwater Galaxias vulgaris fish complex. Genome-wide comparisons reveal significant increases in introgression associated with increasing geographic overlap of taxa. The estimator DP, which assesses the net proportion of a genome originating from introgression, shows a particularly strong relationship with biogeographic overlap (R2 = .43; p = .005). Our analyses nevertheless reveal surprisingly substantial signatures of introgression among taxa that currently have disjunct distributions within drainages (e.g., separate subcatchments). These "anomalies" imply that current biogeography is not always an accurate predictor of introgression history. Our study suggests that both modern and ancient biogeographic shifts, including recent anthropogenic range fragmentation and tectonically driven riven capture events, have influenced introgression histories in this dynamic freshwater fish radiation.

ebony underpins Batesian mimicry in melanic stoneflies

The evolution of Batesian mimicry - whereby harmless species avoid predation through their resemblance to harmful species - has long intrigued biologists. In rare cases, Batesian mimicry is linked to intraspecific colour variation, in which only some individuals within a population resemble a noxious 'model'. Here, we assess intraspecific colour variation within a widespread New Zealand stonefly, wherein highly melanized individuals of Zelandoperla closely resemble a chemically defended aposematic stonefly, Austroperla cyrene. We assess convergence in the colour pattern of these two species, compare their relative palatability to predators, and use genome-wide association mapping to assess the genetic basis of this resemblance. Our analysis reveals that melanized Zelandoperla overlap significantly with Austroperla in colour space but are significantly more palatable to predators, implying that they are indeed Batesian mimics. Analysis of 194,773 genome-wide SNPs reveals an outlier locus (ebony) strongly differentiating melanic versus non-melanic Zelandoperla. Genotyping of 338 specimens from a single Zelandoperla population indicates that ebony explains nearly 70% of the observed variance in melanism. As ebony has a well-documented role in insect melanin biosynthesis, our findings indicate this locus has a conserved function across deeply divergent hexapod lineages. Distributional records suggest a link between the occurrence of melanic Zelandoperla and the forested ecosystems where the model Austroperla is abundant, suggesting the potential for adaptive shifts in this system underpinned by environmental change.

Integrating kelp genomic analyses and geological data to reveal ancient earthquake impacts

Detached buoyant kelp can disperse thousands of kilometres at sea and can colonize newly available shores in the wake of disturbances that wipe out competitors. Localized earthquake uplift can cause extirpation of intertidal kelp populations followed by recolonization. Sources of recolonizing kelp can be detectable in genomic structure of contemporary populations. Our field observations combined with LiDAR mapping identified a previously unrecognized zone of uplifted rocky coastline in a region that is slowly subsiding. Intertidal kelp (Durvillaea antarctica) on the uplifted section of coast is genetically distinctive from nearby populations, with genomic signatures most similar to that of kelp 300 km to the south. Genetic divergence between these locations suggests reproductive isolation for thousands of years. Combined geological and genetic data suggest that this uplift event occurred during one of four major earthquakes between 6000 and 2000 years ago, with one of the younger events most likely. Extirpation of the pre-existing kelp required sudden uplift of approximately 2 metres, precluding several small incremental uplift events. Our results show the power of integrating biological (genomic) analyses with geological data to understand ancient geological processes and their ecological impacts.

Rapid adaptation in a fast-changing world: Emerging insights from insect genomics

Many researchers have questioned the ability of biota to adapt to rapid anthropogenic environmental shifts. Here, we synthesize emerging genomic evidence for rapid insect evolution in response to human pressure. These new data reveal diverse genomic mechanisms (single locus, polygenic, structural shifts; introgression) underpinning rapid adaptive responses to a variety of anthropogenic selective pressures. While the effects of some human impacts (e.g. pollution; pesticides) have been previously documented, here we highlight startling new evidence for rapid evolutionary responses to additional anthropogenic processes such as deforestation. These recent findings indicate that diverse insect assemblages can indeed respond dynamically to major anthropogenic evolutionary challenges. Our synthesis also emphasizes the critical roles of genomic architecture, standing variation and gene flow in maintaining future adaptive potential. Broadly, it is clear that genomic approaches are essential for predicting, monitoring and responding to ongoing anthropogenic biodiversity shifts in a fast-changing world.

Reduced olfactory acuity in recently flightless insects suggests rapid regressive evolution

Background

Insects have exceptionally fast smelling capabilities, and some can track the temporal structure of odour plumes at rates above 100 Hz. It has been hypothesized that this fast smelling capability is an adaptation for flying. We test this hypothesis by comparing the olfactory acuity of sympatric flighted versus flightless lineages within a wing-polymorphic stonefly species.

Results

Our analyses of olfactory receptor neuron responses reveal that recently-evolved flightless lineages have reduced olfactory acuity. By comparing flighted versus flightless ecotypes with similar genetic backgrounds, we eliminate other confounding factors that might have affected the evolution of their olfactory reception mechanisms. Our detection of different patterns of reduced olfactory response strength and speed in independently wing-reduced lineages suggests parallel evolution of reduced olfactory acuity.

Conclusions

These reductions in olfactory acuity echo the rapid reduction of wings themselves, and represent an olfactory parallel to the convergent phenotypic shifts seen under selective gradients in other sensory systems (e.g. parallel loss of vision in cave fauna). Our study provides evidence for the hypothesis that flight poses a selective pressure on the speed and strength of olfactory receptor neuron responses and emphasizes the energetic costs of rapid olfaction.

Seaweed rafts

Jonathan Waters provides an introduction to seaweed rafts and their role in the dispersal of marine and coastal species.

Concordant phylogeographic responses to large-scale coastal disturbance in intertidal macroalgae and their epibiota

Major ecological disturbance events can provide opportunities to assess multispecies responses to upheaval. In particular, catastrophic disturbances that regionally extirpate habitat-forming species can potentially influence the genetic diversity of large numbers of codistributed taxa. However, due to the rarity of such disturbance events over ecological timeframes, the genetic dynamics of multispecies recolonization processes have remained little understood. Here, we use single nucleotide polymorphism (SNP) data from multiple coastal species to track the dynamics of cocolonization events in response to ancient earthquake disturbance in southern New Zealand. Specifically, we use a comparative phylogeographic approach to understand the extent to which epifauna (with varying ecological associations with their macroalgal hosts) share comparable spatial and temporal recolonization patterns. Our study reveals concordant disturbance-related phylogeographic breaks in two intertidal macroalgal species along with two associated epibiotic species (a chiton and an isopod). By contrast, two codistributed species, one of which is an epibiotic amphipod and the other a subtidal macroalga, show few, if any, genetic effects of palaeoseismic coastal uplift. Phylogeographic model selection reveals similar post-uplift recolonization routes for the epibiotic chiton and isopod and their macroalgal hosts. Additionally, codemographic analyses support synchronous population expansions of these four phylogeographically similar taxa. Our findings indicate that coastal paleoseismic activity has driven concordant impacts on multiple codistributed species, with concerted recolonization events probably facilitated by macroalgal rafting. These results highlight that high-resolution comparative genomic data can help reconstruct concerted multispecies responses to recent ecological disturbance.

Genomic signatures of parallel alpine adaptation in recently evolved flightless insects

Natural selection along elevational gradients has potential to drive predictable adaptations across distinct lineages, but the extent of such repeated evolution remains poorly studied for many widespread alpine taxa. We present parallel genomic analyses of two recently evolved flightless alpine insect lineages to test for molecular signatures of repeated alpine adaptation. Specifically, we compare low-elevation vs. alpine stonefly ecotypes from parallel stream populations in which flightless upland ecotypes have been independently derived. We map 67,922 polymorphic genetic markers, generated across 176 Zelandoperla fenestrata specimens from two independent alpine stream populations in New Zealand's Rock and Pillar Range, to a newly developed plecopteran reference genome. Genome-wide scans revealed 31 regions with outlier single nucleotide polymorphisms (SNPs) differentiating lowland vs. alpine ecotypes in Lug Creek, and 37 regions with outliers differentiating ecotypes in Six Mile Creek. Of these regions, 13% (8/60) yielded outlier SNPs across both within-stream ecotype comparisons, implying comparable genomic shifts contribute to this repeated alpine adaptation. Candidate genes closely linked to repeated outlier regions include several with documented roles in insect wing-development (e.g., dishevelled), suggesting that they may contribute to repeated alpine wing reduction. Additional candidate genes have been shown to influence insect fecundity (e.g., ovo) and lifespan (e.g., Mrp4), implying that they might contribute to life history differentiation between upland and lowland ecotypes. Additional outlier genes have potential roles in the evolution of reproductive isolation among ecotypes (hedgehog and Desaturase 1). These results demonstrate how replicated outlier tests across independent lineages can potentially contribute to the discovery of genes underpinning repeated adaptation.

Northward range extension for Durvillaea poha bull kelp: Response to tectonic disturbance?

Understanding the forces that shape species distributions is increasingly important in a fast-changing world. Although major disturbance events can adversely affect natural populations, they can also present new opportunities, for example by opening up habitat for colonization by other lineages. Following extensive geographic sampling, we use genomic data to infer a range extension following disturbance for an ecologically important intertidal macroalgal species. Specifically, we genotyped 288 southern bull kelp (Durvillaea) plants from 28 localities across central New Zealand. All specimens from the North Island were expected to be D. antarctica, but unexpectedly 10 samples from four sites were identified as D. poha. Extensive sampling from the northern South Island (105 samples at five locations) confirmed the absence of D. poha north of the Kaikōura Peninsula. The North Island specimens of D. poha therefore reveal a biogeographic disjunction, some 150 km northeast of the nearest (South Island) population of this species. Based on strong geographic correspondence between these North Island samples and historic disturbance, we infer that tectonic upheaval, particularly earthquake-generated landslides, likely extirpated local D. antarctica and created an opportunity for a northward range expansion event by D. poha. Close phylogenomic relationships between this new North Island population and South Island samples support a geologically recent northward expansion, rather than a deeper evolutionary origin. These findings indicate the potential of large-scale disturbances to facilitate sudden biogeographic range expansions, and they emphasize the ability of genomic analyses with fine-scale sampling to reveal long-lasting signatures of past disturbance, dispersal, and colonization.

Anthropogenic evolution in an insect wing polymorphism following widespread deforestation

Anthropogenic environmental change can underpin major shifts in natural selective regimes, and can thus alter the evolutionary trajectories of wild populations. However, little is known about the evolutionary impacts of deforestation-one of the most pervasive human-driven changes to terrestrial ecosystems globally. Absence of forest cover (i.e. exposure) has been suggested to play a role in selecting for insect flightlessness in montane ecosystems. Here, we capitalize on human-driven variation in alpine treeline elevation in New Zealand to test whether anthropogenic deforestation has caused shifts in the distributions of flight-capable and flightless phenotypes in a wing-polymorphic lineage of stoneflies from the Zelandoperla fenestrata species complex. Transect sampling revealed sharp transitions from flight-capable to flightless populations with increasing elevation. However, these phenotypic transitions were consistently delineated by the elevation of local treelines, rather than by absolute elevation, providing a novel example of human-driven evolution in response to recent deforestation. The inferred rapid shifts to flightlessness in newly deforested regions have implications for the evolution and conservation of invertebrate biodiversity.

Reinventing the wheel? Reassessing the roles of gene flow, sorting and convergence in repeated evolution

Biologists have long been intrigued by apparently predictable and repetitive evolutionary trajectories inferred across a variety of lineages and systems. In recent years, high-throughput sequencing analyses have started to transform our understanding of such repetitive shifts. While researchers have traditionally categorized such shifts as either "convergent" or "parallel," based on relatedness of the lineages involved, emerging genomic insights provide an opportunity to better describe the actual evolutionary mechanisms at play. A synthesis of recent genomic analyses confirms that convergence is the predominant driver of repetitive evolution among species, whereas repeated sorting of standing variation is the major driver of repeated shifts within species. However, emerging data reveal numerous notable exceptions to these expectations, with recent examples of de novo mutations underpinning convergent shifts among even very closely related lineages, while repetitive sorting processes have occurred among even deeply divergent taxa, sometimes via introgression. A number of very recent analyses have found evidence for both processes occurring on different scales within taxa. We suggest that the relative importance of convergent versus sorting processes depends on the interplay between gene flow among populations, and phylogenetic relatedness of the lineages involved.

Genomics Reveals Widespread Ecological Speciation in Flightless Insects

Recent genomic analyses have highlighted parallel divergence in response to ecological gradients, but the extent to which altitude can underpin such repeated speciation remains unclear. Wing reduction and flight loss have apparently evolved repeatedly in montane insect assemblages and have been suggested as important drivers of hexapod diversification. We test this hypothesis using genomic analyses of a widespread wing-polymorphic stonefly species complex in New Zealand. We identified over 50,000 polymorphic genetic markers generated across almost 200 Zelandoperla fenestrata stonefly specimens using a newly generated plecopteran reference genome, to reveal widespread parallel speciation between sympatric full-winged and wing-reduced ecotypes. Rather than the existence of a single, widespread, flightless taxon (Zelandoperla pennulata), evolutionary genomic data reveal that wing-reduced upland lineages have speciated repeatedly and independently from full-winged Z. fenestrata. This repeated evolution of reproductive isolation between local ecotype pairs that lack mitochondrial DNA differentiation suggests that ecological speciation has evolved recently. A cluster of outlier single-nucleotide polymorphisms detected in independently wing-reduced lineages, tightly linked in an approximately 85 kb genomic region that includes the developmental “supergene” doublesex, suggests that this “island of divergence” may play a key role in rapid ecological speciation. [Ecological speciation; genome assembly; genomic island of differentiation; genotyping-by-sequencing; incipient species; plecoptera; wing reduction.]

The genomic footprint of coastal earthquake uplift

Theory suggests that catastrophic earth-history events can drive rapid biological evolution, but empirical evidence for such processes is scarce. Destructive geological events such as earthquakes can represent large-scale natural experiments for inferring such evolutionary processes. We capitalized on a major prehistoric (800 yr BP) geological uplift event affecting a southern New Zealand coastline to test for the lasting genomic impacts of disturbance. Genome-wide analyses of three co-distributed keystone kelp taxa revealed that post-earthquake recolonization drove the evolution of novel, large-scale intertidal spatial genetic 'sectors' which are tightly linked to geological fault boundaries. Demographic simulations confirmed that, following widespread extirpation, parallel expansions into newly vacant habitats rapidly restructured genome-wide diversity. Interspecific differences in recolonization mode and tempo reflect differing ecological constraints relating to habitat choice and dispersal capacity among taxa. This study highlights the rapid and enduring evolutionary effects of catastrophic ecosystem disturbance and reveals the key role of range expansion in reshaping spatial genetic patterns.

Hitchhiking consequences for genetic and morphological patterns: the influence of kelp-rafting on a brooding chiton

Onithochiton neglectus is a morphologically variable, brooding chiton inhabiting coastal reefs throughout New Zealand and its Sub-Antarctic Islands. Southern O. neglectus populations are typically associated with buoyant kelp (Durvillaea spp.) and are potentially connected via kelp-rafting. Northern O. neglectus populations are less likely to raft, due to lower numbers of Durvillaea in northern New Zealand. To test for the impact of kelp-rafting on the spatial distribution of variation in O. neglectus, we undertook a combined analysis of morphological and genetic variation across the range of the species. Geometric morphometrics were used to assess shell shape. We detected a northern vs. southern split in shell shape, corresponding to the frequency of the O. neglectus/Durvillaea spp. association. To assess O. neglectus genetic patterns across New Zealand, we estimated phylogenetic trees with nuclear (ITS) and mitochondrial (COI and 16S) markers, which revealed distinct northern and southern lineages, and an additional lineage in central New Zealand. Neither the morphological nor genetic groups match existing O. neglectus subspecies, but are concordant with the patterns of association of O. neglectus with Durvillaea. We suggest that shell shape may be linked to O. neglectus’ regionally variable ecological association with kelp holdfasts.

SNP analyses reveal a diverse pool of potential colonists to earthquake-uplifted coastlines

In species that form dense populations, major disturbance events are expected to increase the chance of establishment for immigrant lineages. Real-time tests of the impact of disturbance on patterns of genetic structure are, however, scarce. Central to testing these concepts is determining the pool of potential immigrants dispersing into a disturbed area. In 2016, a 7.8 magnitude earthquake occurred on the South Island of New Zealand. Affecting approximately 100 km of coastline, this quake caused extensive uplift (several metres high), extirpating many intertidal populations, including keystone intertidal kelp species. Following the uplift, we set out to determine the geographic origins of detached kelp specimens which rafted into the disturbed zone. Specifically, we used genotyping-by-sequencing (GBS) approaches to compare beach-cast southern bull-kelp (Durvillaea antarctica and Durvillaea poha) samples to established populations throughout the species' ranges, and thus infer the geographic origins of potential colonists reaching the disturbed coast. Our findings revealed an ongoing supply of diverse lineages dispersing to the newly uplifted coastline, suggesting potential for establishment of "exotic" lineages following disturbance. Furthermore, we found that some drifting individuals of each species came from far-distant regions, some >1,200 km away. These results show that diverse lineages - in many cases from very distant sources - can compete for new space in the wake of an exceptional disturbance event, illustrating the potential of long-distance dispersal as a key mechanism for reassembly of coastal ecosystems. Furthermore, our findings demonstrate that high-resolution genomic baselines can be used to robustly assign the provenance of dispersing individuals.

Phylogenetic divergence of island biotas: Molecular dates, extinction, and "relict" lineages

Island formation is a key driver of biological evolution, and several studies have used geological ages of islands to calibrate rates of DNA change. However, many islands are home to "relict" lineages whose divergence apparently pre-dates island age. The geologically dynamic New Zealand (NZ) archipelago sits upon the ancient, largely submerged continent Zealandia, and the origin and age of its distinctive biota have long been contentious. While some researchers have interpreted NZ's biota as equivalent to that of a post-Oligocene island, a recent review of genetic studies identified a sizeable proportion of pre-Oligocene "relict" lineages, concluding that much of the biota survived an incomplete drowning event. Here, we assemble comparable genetic divergence data sets for two recently formed South Pacific archipelagos (Lord Howe; Chatham Islands) and demonstrate similarly substantial proportions of relict lineages. Similar to the NZ biota, our island reviews provide surprisingly little evidence for major genetic divergence "pulses" associated with island emergence. The dominance of Quaternary divergence estimates in all three biotas may highlight the importance of rapid biological turnover and new arrivals in response to recent climatic and/or geological disturbance and change. We provide a schematic model to help account for discrepancies between expected versus observed divergence-date distributions for island biotas, incorporating the effects of both molecular dating error and lineage extinction. We conclude that oceanic islands can represent both evolutionary "cradles" and "museums" and that the presence of apparently archaic island lineages does not preclude dispersal origins.

Comparative transcriptomic analysis of a wing-dimorphic stonefly reveals candidate wing loss genes

Background

The genetic basis of wing development has been well characterised for model insect species, but remains poorly understood in phylogenetically divergent, non-model taxa. Wing-polymorphic insect species potentially provide ideal systems for unravelling the genetic basis of secondary wing reduction. Stoneflies (Plecoptera) represent an anciently derived insect assemblage for which the genetic basis of wing polymorphism remains unclear. We undertake quantitative RNA-seq of sympatric full-winged versus vestigial-winged nymphs of a widespread wing-dimorphic New Zealand stonefly, Zelandoperla fenestrata, to identify genes potentially involved in wing development and secondary wing loss.

Results

Our analysis reveals substantial differential expression of wing-development genes between full-winged versus vestigial-winged stonefly ecotypes. Specifically, of 23 clusters showing significant similarity to Drosophila wing development-related genes and their pea aphid orthologues, nine were significantly upregulated in full-winged stonefly ecotypes, whereas only one cluster (teashirt) was substantially upregulated in the vestigial-winged ecotype.

Conclusions

These findings suggest remarkable conservation of key wing-development pathways throughout 400 Ma of insect evolution. The finding that two Juvenile Hormone pathway clusters were significantly upregulated in vestigial-winged Zelandoperla supports the hypothesis that Juvenile Hormone may play a key role in modulating insect wing polymorphism, as has previously been suggested for other insect lineages.

More than the eye can see: Genomic insights into the drivers of genetic differentiation in Royal/Macaroni penguins across the Southern Ocean

The study of systematics in wide-ranging seabirds can be challenging due to the vast geographic scales involved, as well as the possible discordance between molecular, morphological and behavioral data. In the Southern Ocean, macaroni penguins (Eudyptes chrysolophus) are distributed over a circumpolar range including populations in Antarctic and sub-Antarctic areas. Macquarie Island, in its relative isolation, is home to a closely related endemic taxon - the royal penguin (Eudyptes schlegeli), which is distinguishable from E. chrysolophus mainly by facial coloration. Although these sister taxa are widely accepted as representing distinct species based on morphological grounds, the extent of their genome-wide differentiation remains uncertain. In this study, we use genome-wide Single Nucleotide Polymorphisms to test genetic differentiation between these geographically isolated taxa and evaluate the main drivers of population structure among breeding colonies of macaroni/royal penguins. Genetic similarity observed between macaroni and royal penguins suggests they constitute a single evolutionary unit. Nevertheless, royal penguins exhibited a tendency to cluster only with macaroni individuals from Kerguelen Island, suggesting that dispersal occurs mainly between these neighboring colonies. A stepping stone model of differentiation of macaroni/royal populations was further supported by a strong pattern of isolation by distance detected across its whole distribution range, possibly driven by large geographic distances between colonies as well as natal philopatry. However, we also detected intraspecific genomic differentiation between Antarctic and sub-Antarctic populations of macaroni penguins, highlighting the role of environmental factors together with geographic distance in the processes of genetic differentiation between Antarctic and sub-Antarctic waters.

Mitogenomes Uncover Extinct Penguin Taxa and Reveal Island Formation as a Key Driver of Speciation

The emergence of islands has been linked to spectacular radiations of diverse organisms. Although penguins spend much of their lives at sea, they rely on land for nesting, and a high proportion of extant species are endemic to geologically young islands. Islands may thus have been crucial to the evolutionary diversification of penguins. We test this hypothesis using a fossil-calibrated phylogeny of mitochondrial genomes (mitogenomes) from all extant and recently extinct penguin taxa. Our temporal analysis demonstrates that numerous recent island-endemic penguin taxa diverged following the formation of their islands during the Plio-Pleistocene, including the Galápagos (Galápagos Islands), northern rockhopper (Gough Island), erect-crested (Antipodes Islands), Snares crested (Snares) and royal (Macquarie Island) penguins. Our analysis also reveals two new recently extinct island-endemic penguin taxa from New Zealand's Chatham Islands: Eudyptes warhami sp. nov. and a dwarf subspecies of the yellow-eyed penguin, Megadyptes antipodes richdalei ssp. nov. Eudyptes warhami diverged from the Antipodes Islands erect-crested penguin between 1.1 and 2.5 Ma, shortly after the emergence of the Chatham Islands (∼3 Ma). This new finding of recently evolved taxa on this young archipelago provides further evidence that the radiation of penguins over the last 5 Ma has been linked to island emergence. Mitogenomic analyses of all penguin species, and the discovery of two new extinct penguin taxa, highlight the importance of island formation in the diversification of penguins, as well as the extent to which anthropogenic extinctions have affected island-endemic taxa across the Southern Hemisphere's isolated archipelagos.

Genotyping-by-sequencing supports a genetic basis for wing reduction in an alpine New Zealand stonefly

Wing polymorphism is a prominent feature of numerous insect groups, but the genomic basis for this diversity remains poorly understood. Wing reduction is a commonly observed trait in many species of stoneflies, particularly in cold or alpine environments. The widespread New Zealand stonefly Zelandoperla fenestrata species group (Z. fenestrata, Z. tillyardi, Z. pennulata) contains populations ranging from fully winged (macropterous) to vestigial-winged (micropterous), with the latter phenotype typically associated with high altitudes. The presence of flightless forms on numerous mountain ranges, separated by lowland fully winged populations, suggests wing reduction has occurred multiple times. We use Genotyping by Sequencing (GBS) to test for genetic differentiation between fully winged (n = 62) and vestigial-winged (n = 34) individuals, sampled from a sympatric population of distinct wing morphotypes, to test for a genetic basis for wing morphology. While we found no population genetic differentiation between these two morphotypes across 6,843 SNP loci, we did detect several outlier loci that strongly differentiated morphotypes across independent tests. These findings indicate that small regions of the genome are likely to be highly differentiated between morphotypes, suggesting a genetic basis for wing reduction. Our results provide a clear basis for ongoing genomic analysis to elucidate critical regulatory pathways for wing development in Pterygota.

Ancient DNA of crested penguins: Testing for temporal genetic shifts in the world's most diverse penguin clade

Human impacts have substantially reduced avian biodiversity in many parts of the world, particularly on isolated islands of the Pacific Ocean. The New Zealand archipelago, including its five subantarctic island groups, holds breeding grounds for a third of the world's penguin species, including several representatives of the diverse crested penguin genus Eudyptes. While this species-rich genus has been little studied genetically, recent population estimates indicate that several Eudyptes taxa are experiencing demographic declines. Although crested penguins are currently limited to southern regions of the New Zealand archipelago, prehistoric fossil and archaeological deposits suggest a wider distribution during prehistoric times, with breeding ranges perhaps extending to the North Island. Here, we analyse ancient, historic and modern DNA sequences to explore two hypotheses regarding the recent history of Eudyptes in New Zealand, testing for (1) human-driven extinction of Eudyptes lineages; and (2) reduced genetic diversity in surviving lineages. From 83 prehistoric bone samples, each tentatively identified as 'Eudyptes spp.', we genetically identified six prehistoric penguin taxa from mainland New Zealand, including one previously undescribed genetic lineage. Moreover, our Bayesian coalescent analyses indicated that, while the range of Fiordland crested penguin (E. pachyrhynchus) may have contracted markedly over the last millennium, genetic DNA diversity within this lineage has remained relatively constant. This result contrasts with human-driven biodiversity reductions previously detected in several New Zealand coastal vertebrate taxa.

Crossing the front: contrasting storm-forced dispersal dynamics revealed by biological, geological and genetic analysis of beach-cast kelp

The subtropical front (STF) generally represents a substantial oceanographic barrier to dispersal between cold-sub-Antarctic and warm-temperate water masses. Recent studies have suggested that storm events can drastically influence marine dispersal and patterns. Here we analyse biological and geological dispersal driven by two major, contrasting storm events in southern New Zealand, 2017. We integrate biological and physical data to show that a severe southerly system in July 2017 disrupted this barrier by promoting movement of substantial numbers of southern sub-Antarctic Durvillaea kelp rafts across the STF, to make landfall in mainland NZ. By contrast, a less intense easterly storm (Cyclone Cook, April 2017) resulted in more moderate dispersal distances, with minimal dispersal between the sub-Antarctic and mainland New Zealand. These quantitative analyses of approximately 200 freshly beach-cast kelp specimens indicate that storm intensity and wind direction can strongly influence marine dispersal and landfall outcomes.

Rafting dispersal in a brooding southern sea star (Asteroidea : Anasterias)

Marine biogeographers have long speculated that macroalgal rafting presents a dispersal mechanism for brooding marine invertebrates of the Southern Ocean, but few direct observations of rafting by echinoderm taxa have been documented. Here we report rafting of the brooding benthic sea star Anasterias suteri, along with two mollusc taxa (Onithochiton neglectus – also a brooder – and Cantharidus roseus), on detached bull-kelp Durvillaea antarctica in Foveaux Strait, southern New Zealand. The rafting journey, intercepted at sea, likely lasted for 2–3 weeks and may have covered several hundred kilometres. We use DNA sequences, together with meteorological and prevailing oceanographic data, to infer the likely Fiordland (mainland) origins of the raft and its epifauna. This rafting dispersal mechanism provides an explanation for the broad (circum-subantarctic) but disjunct distribution of brooding Anasterias populations, and for the genetic connectivity observed between their populations.

Does wing reduction influence the relationship between altitude and insect body size? A case study using New Zealand's diverse stonefly fauna

Researchers have long been intrigued by evolutionary processes that explain biological diversity. Numerous studies have reported strong associations between animal body size and altitude, but insect analyses have often yielded equivocal results. Here, we analyze a collection database of New Zealand's diverse endemic stonefly fauna (106 species across 21 genera) to test for relationships between altitude and plecopteran body size. This insect assemblage includes a variety of wing-reduced (26 spp) and fully winged (80 spp) taxa and covers a broad range of altitudes (0-2,000 m). We detected significant relationships between altitude and body size for wing-reduced, but not fully winged, stonefly taxa. These results suggest that, while the maintenance of flight apparatus might place a constraint on body size in some fully winged species, the loss of flight may free insects from this evolutionary constraint. We suggest that rapid switches in insect dispersal ability may facilitate rapid evolutionary shifts across a number of biological attributes and may explain the inconsistent results from previous macroecological analyses of insect assemblages.

Ancient DNA and morphometric analysis reveal extinction and replacement of New Zealand's unique black swans

Prehistoric human impacts on megafaunal populations have dramatically reshaped ecosystems worldwide. However, the effects of human exploitation on smaller species, such as anatids (ducks, geese, and swans) are less clear. In this study we apply ancient DNA and osteological approaches to reassess the history of Australasia's iconic black swans (Cygnus atratus) including the palaeo-behaviour of prehistoric populations. Our study shows that at the time of human colonization, New Zealand housed a genetically, morphologically, and potentially ecologically distinct swan lineage (C. sumnerensis, Poūwa), divergent from modern (Australian) C. atratus Morphological analyses indicate C. sumnerensis exhibited classic signs of the 'island rule' effect, being larger, and likely flight-reduced compared to C. atratus Our research reveals sudden extinction and replacement events within this anatid species complex, coinciding with recent human colonization of New Zealand. This research highlights the role of anthropogenic processes in rapidly reshaping island ecosystems and raises new questions for avian conservation, ecosystem re-wilding, and de-extinction.

A morphological and phylogenetic investigation into divergence among sympatric Australian southern bull kelps (Durvillaea potatorum and D. amatheiae sp. nov.)

Genetic analyses can reveal a wealth of hitherto undiscovered cryptic biodiversity. For co-occurring and morphologically similar species, the combination of molecular, ecological and morphological analyses provides an excellent opportunity for understanding some of the processes that can lead to divergence and speciation. The Australian endemic brown macroalga Durvillaea potatorum (Phaeophyceae) was examined with a combination of genetic and morphological approaches to confirm the presence of two separate species and to infer the processes that led to their divergence. A total of 331 individuals from 11 sites around coastal Tasmania were collected and measured in situ for a range of morphological and ecological characteristics. Tissue samples were also collected for each individual to allow genetic analyses using mitochondrial (COI) and nuclear (28S) markers. Genetic analyses confirmed the presence of two deeply divergent clades. The significant morphological differentiation, despite high levels of intra-lineage variability, further supported their recognition as distinct species. We describe a new species, D. amatheiae sp. nov., which is characterised by a narrower and proportionately shorter stipe, shorter total length, and higher number of stipitate lateral blades and branches than D. potatorum (sensu stricto). The occurrence of both species in sympatry along Tasmania's eastern and western coasts, as well as their contrasting patterns of haplotype diversity, supports a hypothesis of geographical isolation, allopatric speciation and subsequent secondary contact in response to sea level and ocean current change throughout the Pleistocene glaciation cycles. This research contributes to resolving the phylogenetic relationships, taxonomy and evolution of the ecologically keystone kelp genus Durvillaea.

Invader or resident? Ancient-DNA reveals rapid species turnover in New Zealand little penguins

The expansion of humans into previously unoccupied parts of the globe is thought to have driven the decline and extinction of numerous vertebrate species. In New Zealand, human settlement in the late thirteenth century AD led to the rapid demise of a distinctive vertebrate fauna, and also a number of 'turnover' events where extinct lineages were subsequently replaced by closely related taxa. The recent genetic detection of an Australian little penguin (Eudyptula novaehollandiae) in southeastern New Zealand may potentially represent an additional 'cryptic' invasion. Here we use ancient-DNA (aDNA) analysis and radiocarbon dating of pre-human, archaeological and historical Eudyptula remains to reveal that the arrival of E. novaehollandiae in New Zealand probably occurred between AD 1500 and 1900, following the anthropogenic decline of its sister taxon, the endemic Eudyptula minor. This rapid turnover event, revealed by aDNA, suggests that native species decline can be masked by invasive taxa, and highlights the potential for human-mediated biodiversity shifts.

Managing shifting species: Ancient DNA reveals conservation conundrums in a dynamic world

The spread of exotic species represents a major driver of biological change across the planet. While dispersal and colonization are natural biological processes, we suggest that the failure to recognize increasing rates of human-facilitated self-introductions may represent a threat to native lineages. Notably, recent biogeographic analyses have revealed numerous cases of biological range shifts in response to anthropogenic impacts and climate change. In particular, ancient DNA analyses have revealed several cases in which lineages traditionally thought to be long-established "natives" are in fact recent colonizers. Such range expansion events have apparently occurred in response to human-mediated native biodiversity declines and ecosystem change, particularly in recently colonized, isolated ecosystems such as New Zealand. While such events can potentially boost local biodiversity, the spread of exotic lineages may also hasten the decline of indigenous species, so it is essential that conservation managers recognize these rapid biotic shifts.​.

The importance of replicating genomic analyses to verify phylogenetic signal for recently evolved lineages

Genomewide SNP data generated by nontargeted methods such as RAD and GBS are increasingly being used in phylogenetic and phylogeographic analyses. When these methods are used in the absence of a reference genome, however, little is known about the locations and evolution of the SNPs. In using such data to address phylogenetic questions, researchers risk drawing false conclusions, particularly if a representative number of SNPs is not obtained. Here, we empirically test the robustness of phylogenetic inference based on SNP data for closely related lineages. We conducted a genomewide analysis of 75 712 SNPs, generated via GBS, of southern bull-kelp (Durvillaea). Durvillaea chathamensis co-occurs with D. antarctica on Chatham Island, but the two species have previously been found to be so genetically similar that the status of the former has been questioned. Our results show that D. chathamensis, which differs from D. antarctica ecologically as well as morphologically, is indeed a reproductively isolated species. Furthermore, our replicated analyses show that D. chathamensis cannot be reliably distinguished phylogenetically from closely related D. antarctica using subsets (ranging in size from 400 to 10 000 sites) of the 40 912 parsimony-informative SNPs in our data set and that bootstrap values alone can give misleading impressions of the strength of phylogenetic inferences. These results highlight the importance of independently replicating SNP analyses to verify that phylogenetic inferences based on nontargeted SNP data are robust. Our study also demonstrates that modern genomic approaches can be used to identify cases of recent or incipient speciation that traditional approaches (e.g. Sanger sequencing of a few loci) may be unable to detect or resolve.

Transoceanic dispersal and cryptic diversity in a cosmopolitan rafting nudibranch

The aeolid nudibranch Fiona pinnata (Eschscholtz, 1831) is an obligate rafter that occurs exclusively on macroalgal rafts and other floating substrata, and has a seemingly cosmopolitan marine distribution. Mitochondrial (mtDNA) and nuclear DNA sequence data were generated from specimens collected worldwide to test for global connectivity in this species. Phylogeographic analyses revealed three divergent mtDNA lineages, two of which were abundant and widespread. One of these lineages has an apparent circumequatorial distribution, whereas the other may have an antitropical distribution within the Pacific Ocean. Low genetic divergences within each lineage suggest that rafting can mediate dispersal across transoceanic scales. A third, highly divergent, lineage was detected in a single Indonesian specimen. Broadly concordant phylogeographic relationships were detected for the nuclear ITS1 region, with distinct tropical versus antitropical lineages observed. The substantial genetic divergences and largely allopatric distributions observed among the F. pinnata lineages suggest that they represent a species complex.