Characterization of a hybrid zone between two chromosomal races of the weta Hemideina thoracica following a geologically recent volcanic eruption

Spatial Variation of Acanthophlebia cruentata (Ephemeroptera), a Mayfly Endemic to Te Ika-a-Māui—North Island of Aotearoa, New Zealand

Simple Summary

Patterns of variation within a widespread species can provide evidence of population history. Adjacent stable populations with gene flow show clinal genetic divergence forming a pattern of isolation by distance. Populations that have grown due to an increase in potential habitat such as forest cover since the last glacial maximum will have low genetic variation showing patterns of range expansion. The mayfly Acanthophlebia cruentata of Aotearoa, New Zealand, is widespread in streams of North Island but absent from the cooler South Island. Mayfly nymphs are restricted to streams but adults fly, facilitating gene flow among catchments. We detected higher genetic diversity at lower latitudes of this mayfly’s range compared to most of its distribution, concordant with predictions of limited forest cover in New Zealand during Pleistocene glacial periods. A signature of recent range expansion was observed in the higher latitudes. Despite initial observation suggesting mayfly size correlated with latitude, we found sex, elevation and sampling date were significant predictors of size, and some size variation is also explained by three regional groups based on haplotype distribution.

Abstract

The mayfly Acanthophlebia cruentata of Aotearoa, New Zealand, is widespread in Te Ika-a-Māui North Island streams, but has never been collected from South Island despite land connection during the last glacial maximum. Population structure of this mayfly might reflect re-colonisation after volcanic eruptions in North Island c1800 years ago, climate cycling or conceal older, cryptic diversity. We collected population samples from 33 locations to estimate levels of population genetic diversity and to document phenotypic variation. Relatively low intraspecific haplotype divergence was recorded among mitochondrial cytb sequences from 492 individuals, but these resolved three geographic-haplotype regions (north, west, east). We detected a signature of isolation by distance at low latitudes (north) but evidence of recent population growth in the west and east. We did not detect an effect of volcanic eruptions but infer range expansion into higher latitudes from a common ancestor during the last glacial period. As judged from wing length, both sexes of adult mayflies were larger at higher elevation and we found that haplotype region was also a significant predictor of Acanthophlebia cruentata size. This suggests that our mitochondrial marker is concordant with nuclear genetic differences that might be explained by founder effect during range expansion.

Dispersal Predicts Hybrid Zone Widths across Animal Diversity: Implications for Species Borders under Incomplete Reproductive Isolation

Hybrid zones occur as range boundaries for many animal taxa. One model for how hybrid zones form and stabilize is the tension zone model, a version of which predicts that hybrid zone widths are determined by a balance between random dispersal into hybrid zones and selection against hybrids. Here, we examine whether random dispersal and proxies for selection against hybrids (genetic distances between hybridizing pairs) can explain variation in hybrid zone widths across 131 hybridizing pairs of animals. We show that these factors alone can explain ∼40% of the variation in zone width among animal hybrid zones, with dispersal explaining far more of the variation than genetic distances. Patterns within clades were idiosyncratic. Genetic distances predicted hybrid zone widths particularly well for reptiles, while this relationship was opposite tension zone predictions in birds. Last, the data suggest that dispersal and molecular divergence set lower bounds on hybrid zone widths in animals, indicating that there are geographic restrictions on hybrid zone formation. Overall, our analyses reinforce the fundamental importance of dispersal in hybrid zone formation and more generally in the ecology of range boundaries.

Closing the gap: Avian lineage splits at a young, narrow seaway imply a protracted history of mixed population response

The evolutionary significance of spatial habitat gaps has been well recognized since Alfred Russel Wallace compared the faunas of Bali and Lombok. Gaps between islands influence population structuring of some species, and flightless birds are expected to show strong partitioning even where habitat gaps are narrow. We examined the population structure of the most numerous living flightless land bird in New Zealand, Weka (Gallirallus australis). We surveyed Weka and their feather lice in native and introduced populations using genetic data gathered from DNA sequences of mitochondrial genes and nuclear β-fibrinogen and five microsatellite loci. We found low genetic diversity among extant Weka population samples. Two genetic clusters were evident in the mtDNA from Weka and their lice, but partitioning at nuclear loci was less abrupt. Many formerly recognized subspecies/species were not supported; instead, we infer one subspecies for each of the two main New Zealand islands. Although currently range restricted, North Island Weka have higher mtDNA diversity than the more wide-ranging southern Weka. Mismatch and neutrality statistics indicate North Island Weka experienced rapid and recent population reduction, while South Island Weka display the signature of recent expansion. Similar haplotype data from a widespread flying relative of Weka and other New Zealand birds revealed instances of North Island-South Island partitioning associated with a narrow habitat gap (Cook Strait). However, contrasting patterns indicate priority effects and other ecological factors have a strong influence on spatial exchange at this scale.

Little or no gene flow despite F1 hybrids at two interspecific contact zones

Hybridization can create the selective force that promotes assortative mating but hybridization can also select for increased hybrid fitness. Gene flow resulting from hybridization can increase genetic diversity but also reduce distinctiveness. Thus the formation of hybrids has important implications for long‐term species coexistence. This study compares the interaction between the tree wētā Hemideina thoracica and its two neighboring species; H. crassidens and H. trewicki. We examined the ratio of parent and hybrid forms in natural areas of sympatry. Individuals with intermediate phenotype were confirmed as first generation hybrids using nine independent genetic markers. Evidence of gene flow from successful hybridization was sought from the distribution of morphological and genetic characters. Both species pairs appear to be largely retaining their own identity where they live in sympatry, each with a distinct karyotype. Hemideina thoracica and H. trewicki are probably reproductively isolated, with sterile F₁ hybrids. This species pair shows evidence of niche differences with adult size and timing of maturity differing where Hemideina thoracica is sympatric with H. trewicki. In contrast, evidence of a low level of introgression was detected in phenotypes and genotypes where H. thoracica and H. crassidens are sympatric. We found no evidence of size divergence although color traits in combination with hind tibia spines reliably distinguish the two species. This species pair show a bimodal hybrid zone in the absence of assortative mating and possible sexual exclusion by H. thoracica males in the formation of F₁ hybrids.

Little and large: body size and genetic clines in a New Zealand gecko (Woodworthia maculata) along a coastal transect

Clinal variation can result from primary differentiation or secondary contact and determining which of these two processes is responsible for the existence of a cline is not a trivial problem. Samples from a coastal transect of New Zealand geckos (Woodworthia maculata) identified for the first time a body size cline 7–10 km wide. The larger geckos are almost twice the mass of the small adult geckos. Clines in allele and haplotype frequency were found at two of the four genetic loci examined. Estimated width of the morphological cline was concordant with neither the narrower mtDNA cline (3–7 m) nor the wider nuclear cline (RAG-2; 34–42 km), and cline centers were not coincident. Although the body size cline is narrow compared to the entire range of the species, it is 2–3 orders of magnitude greater than estimates of dispersal distance per generation for these geckos. No evidence of assortative mating, nor of hybrid disadvantage was identified, thus there is little evidence to infer that endogenous selection is maintaining a hybrid zone. We cannot distinguish secondary contact from primary origin of this body size cline but conclude that secondary contact is likely due to the occurrence of mtDNA haplotypes from three distinct clades within the coastal transect and the presence of two frequency clines within this region.

The Invertebrate Life of New Zealand: A Phylogeographic Approach

Phylogeography contributes to our knowledge of regional biotas by integrating spatial and genetic information. In New Zealand, comprising two main islands and hundreds of smaller ones, phylogeography has transformed the way we view our biology and allowed comparison with other parts of the world. Here we review studies on New Zealand terrestrial and freshwater invertebrates. We find little evidence of congruence among studies of different taxa; instead there are signatures of partitioning in many different regions and expansion in different directions. A number of studies have revealed unusually high genetic distances within putative species, and in those where other data confirm this taxonomy, the revealed phylogeographic structure contrasts with northern hemisphere continental systems. Some taxa show a signature indicative of Pliocene tectonic events encompassing land extension and mountain building, whereas others are consistent with range expansion following the last glacial maximum (LGM) of the Pleistocene. There is some indication that montane taxa are more partitioned than lowland ones, but this observation is obscured by a broad range of patterns within the sample of lowland/forest taxa. We note that several geophysical processes make similar phylogeographic predictions for the same landscape, rendering confirmation of the drivers of partitioning difficult. Future multi-gene analyses where applied to testable alternative hypotheses may help resolve further the rich evolutionary history of New Zealand's invertebrates.

A comparison of five hybrid zones of the weta Hemideina thoracica (Orthoptera: Anostostomatidae): degree of cytogenetic differentiation fails to predict zone width

Tension zones are maintained by the interaction between selection against hybrids and dispersal of individuals. Investigating multiple hybrid zones within a single species provides the opportunity to examine differences in zone structure on a background of differences in extrinsic factors (e.g., age of the zone, ecology) or intrinsic factors (e.g., chromosomes). The New Zealand tree weta Hemideina thoracica comprises at least eight distinct chromosomal races with diploid numbers ranging from 2n = 11 (XO) to 2n = 23 (XO). Five independent hybrid zones were located that involve races differing from one another by a variety of chromosomal rearrangements. The predicted negative correlation between extent of karyotypic differentiation (measured in terms of both percent of genome and number of rearrangements) and zone width was not found. Conversely, the widest zones were those characterized by two chromosome rearrangements involving up to 35% of the genome. The narrowest zone occurred where the two races differ by a single chromosome rearrangement involving approximately 2% of the genome. The five estimates of chromosomal cline width ranged from 0.5 km to 47 km. A comparative investigation of cline width for both chromosomal and mitochondrial markers revealed a complex pattern of zone characteristics. Three of the five zones in this study showed cline concordance for the nuclear and cytoplasmic markers, and at two of the zones the clines were also coincident. Zones with the widest chromosomal clines had the widest mitochondrial DNA clines. It appears that, even within a single species, the extent of karyotypic differentiation between pairs of races is not a good predictor of the level of disadvantage suffered by hybrids.