Fire and ice: volcanic and glacial impacts on the phylogeography of the New Zealand forest fernAsplenium hookerianum

Genetic Diversity and Phylogeography of the Relict Tree Fern Culcita macrocarpa: Influence of Clonality and Breeding System on Genetic Variation

The tree fern Culcita macrocarpa, a threatened Iberian-Macaronesian endemism, represents the sole European species of the order Cyatheales. Considered a Tertiary relict of European Palaeotropical flora, its evolutionary history and genetic diversity, potentially influenced by presumed high clonal propagation, remain largely unknown. This study elucidates the phylogeographic history of C. macrocarpa, assessing the impact of vegetative reproduction on population dynamics and genetic variability. We provide genetic data from eight newly identified nuclear microsatellite loci and one plastid DNA region for 17 populations spanning the species' range, together with species distribution modeling data. Microsatellites reveal pervasive clonality in C. macrocarpa, which has varied among populations. We assess the impact of clonality on genetic diversity and evaluate how estimates of intra-population genetic diversity indices and genetic structuring are affected by the chosen definition of "individual" (focusing exclusively on genetically distinct individuals, genets, as opposed to considering all independent clonal replicates, ramets). We identify two main population groups, one in the northern Iberian Peninsula and the other in the Macaronesian archipelagos and southern Iberian Peninsula. Within each group, we found relict populations (in the Azores and the Cantabrian Cornice) as well as recent originated populations. This population structure suggests colonization dynamics in which recent populations originated from one or a few genets of relict populations and became established through intra-gametophytic self-fertilization and vegetative expansion. DAPC analysis facilitated the identification of alleles that most significantly contributed to the observed population structure. The current Andalusian populations appear to have originated from colonization events from the Azores and the Cantabrian Cornice. Our findings suggest that C. macrocarpa persisted through the Last Glacial Maximum in two refugia: the Azores and the Cantabrian Cornice. Colonization into new areas occurred presumably from these refuges, generating two large population groups with structured genetic diversity. This study underscores the significance of clonality in establishing new populations and shaping genetic structure.

Legacy of supervolcanic eruptions on population genetic structure of brown kiwi

Supervolcanoes are volcanoes capable of mega-colossal eruptions that emit more than 1,000 km³ of ash and other particles.¹ The earth's most recent mega-colossal eruption was the Oruanui eruption of the Taupo supervolcano 25,580 years before present (YBP) on the central North Island of New Zealand.² This eruption blanketed major swaths of the North Island in thick layers of ash and igneous rock,^2,3 devastating habitats and likely causing widespread population extinctions.^4-7 An additional devastating super-colossal eruption (>100 km³) of the Taupo supervolcano occurred approximately 1,690 YBP.⁸ The impacts of such massive but ephemeral natural disasters on contemporary population genetic structure remain underexplored. Here, we combined data for 4,951 SNPs with spatially explicit demographic and coalescent models within an approximate Bayesian computation framework to test the drivers of genetic structure in brown kiwi (Apteryx mantelli). Our results strongly support the importance of eruptions of the Taupo supervolcano in restructuring pre-existing geographic patterns of population differentiation and genetic diversity. Range shifts due to climatic oscillations-a frequent explanation for genetic structure⁹-are insufficient to fully explain the empirical data. Meanwhile, recent range contraction and fragmentation due to historically documented anthropogenic habitat alteration adds no explanatory power to our models. Our results support a major role for cycles of destruction and post-volcanic recolonization in restructuring the population genomic landscape of brown kiwi and highlight how ancient and ephemeral mega-disasters may leave a lasting legacy on patterns of intraspecific genetic variation.

Phylogeographical Analyses of a Relict Fern of Palaeotropical Flora (Vandenboschia speciosa): Distribution and Diversity Model in Relation to the Geological and Climate Events of the Late Miocene and Early Pliocene

Fern phylogeographic studies have mostly focused on the influence of the Pleistocene climate on fern distributions and the prevalence of long-distance dispersal. The effect of pre-Pleistocene events on the distributions of fern species is largely unexplored. Here, we elucidate a hypothetical scenario for the evolutionary history of Vandenboschia speciosa, hypothesised to be of Tertiary palaeotropical flora with a peculiar perennial gametophyte. We sequenced 40 populations across the species range in one plastid region and two variants of the nuclear gapCp gene and conducted time-calibrated phylogenetic, phylogeographical, and species distribution modelling analyses. Vandenboschia speciosa is an allopolyploid and had a Tertiary origin. Late Miocene aridification possibly caused the long persistence in independent refugia on the Eurosiberian Atlantic and Mediterranean coasts, with the independent evolution of gene pools resulting in two evolutionary units. The Cantabrian Cornice, a major refugium, could also be a secondary contact zone during Quaternary glacial cycles. Central European populations resulted from multiple post-glacial, long-distance dispersals. Vandenboschia speciosa reached Macaronesia during the Pliocene–Pleistocene, with a phylogeographical link between the Canary Islands, Madeira, and southern Iberia, and between the Azores and northwestern Europe. Our results support the idea that the geological and climate events of the Late Miocene/Early Pliocene shifted Tertiary fern distribution patterns in Europe.

Explaining large mitochondrial sequence differences within a population sample

Mitochondrial DNA sequence is frequently used to infer species' boundaries, as divergence is relatively rapid when populations are reproductively isolated. However, the shared history of a non-recombining gene naturally leads to correlation of pairwise differences, resulting in mtDNA clusters that might be mistaken for evidence of multiple species. There are four distinct processes that can explain high levels of mtDNA sequence difference within a single sample. Here, we examine one case in detail as an exemplar to distinguish among competing hypotheses. Within our sample of tree wētā (Hemideina crassidens; Orthoptera), we found multiple mtDNA haplotypes for a protein-coding region (cytb/ND1) that differed by a maximum of 7.9%. From sequencing the whole mitochondrial genome of two representative individuals, we found evidence of constraining selection. Heterozygotes were as common as expected under random mating at five nuclear loci. Morphological traits and nuclear markers did not resolve the mtDNA groupings of individuals. We concluded that the large differences found among our sample of mtDNA sequences were simply owing to a large population size over an extended period of time allowing an equilibrium between mutation and drift to retain a great deal of genetic diversity within a single species.

Phylogeographical pattern and evolutionary history of an important Peninsular Malaysian timber species, Neobalanocarpus heimii (Dipterocarpaceae)

Tectonic movements, climatic oscillations, and marine transgressions during the Cenozoic have had a dramatic effect on the biota of the tropical rain forest. This study aims to reveal the phylogeography and evolutionary history of a Peninsular Malaysian endemic tropical timber species, Neobalanocarpus heimii (Dipterocarpaceae). A total of 32 natural populations of N. heimii, with 8 samples from each population were investigated. Fifteen haplotypes were identified from five noncoding chloroplast DNA (cpDNA) regions. Overall, two major genealogical cpDNA lineages of N. heimii were elucidated: a widespread southern and a northern region. The species is predicted to have survived in multiple refugia during climatic oscillations: the northwestern region (R1), the northeastern region (R2), and the southern region (R3). These putative glacial refugia exhibited higher levels of genetic diversity, population differentiation, and the presence of unique haplotypes. Recolonization of refugia R1 and R2 could have first expanded into the northern region and migrated both northeastwards and northwestwards. Meanwhile, recolonization of N. heimii throughout the southern region could have commenced from refugia R3 and migrated toward the northeast and northwest, respectively. The populations of Tersang, Pasir Raja, and Rotan Tunggal exhibited remarkably high haplotype diversity, which could have been the contact zones that have received an admixture of gene pools from the northerly and also southerly regions. As a whole, the populations of N. heimii derived from glacial refugia and contact zones should be considered in the conservation strategies in order to safeguard the long-term survival of the species.

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.

Ancient DNA analyses reveal contrasting phylogeographic patterns amongst kiwi (Apteryx spp.) and a recently extinct lineage of spotted kiwi

The little spotted kiwi (Apteryx owenii) is a flightless ratite formerly found throughout New Zealand but now greatly reduced in distribution. Previous phylogeographic studies of the related brown kiwi (A. mantelli, A. rowi and A. australis), with which little spotted kiwi was once sympatric, revealed extremely high levels of genetic structuring, with mitochondrial DNA haplotypes often restricted to populations. We surveyed genetic variation throughout the present and pre-human range of little spotted kiwi by obtaining mitochondrial DNA sequences from contemporary and ancient samples. Little spotted kiwi and great spotted kiwi (A. haastii) formed a monophyletic clade sister to brown kiwi. Ancient samples of little spotted kiwi from the northern North Island, where it is now extinct, formed a lineage that was distinct from remaining little spotted kiwi and great spotted kiwi lineages, potentially indicating unrecognized taxonomic diversity. Overall, little spotted kiwi exhibited much lower levels of genetic diversity and structuring than brown kiwi, particularly through the South Island. Our results also indicate that little spotted kiwi (or at least hybrids involving this species) survived on the South Island mainland until more recently than previously thought.

Efficient long-distance gene flow into an isolated relict oak stand

Geographically isolated and small populations outside a species' central distribution range are likely to be of major importance to a species' ability to quickly adjust its distribution range to global change dynamics. Gene flow from the outside plays a pivotal role in the fate of these marginal populations. It has been proposed that spatial fragmentation and perceived geographic isolation do not necessarily reflect a loss of genetic connectivity in tree species. However, the spatial limits of long-distance gene flow, as well as its magnitude and impact, are still generally unknown. In the present study, we analyzed long-distance pollen-mediated gene flow into an isolated relict stand consisting of 7 individuals of Quercus robur L. based on a total sample of 177 trees and 9 microsatellite loci. We show that pollen-mediated gene flow across more than 80 km in this wind-pollinated tree species contributed at least 35% of all successful pollinations in the investigated isolated and small oak stand at the eastern limit of the species' distribution. The observed pollen immigration shaped the genetic diversity of acorn progenies in the stand and might explain the comparably high genetic diversity in the persisting adult population.

Two micro-scale protocols for the isolation of DNA from polysaccharide-rich plant tissue

The high polysaccharide content of some plant species hinders the successful isolation of their DNA. As an alternative to the macro-extraction methods previously published for polysaccharide-rich plants, we present two techniques (STE/CTAB and HEPES/CTAB), which are performed in microcentrifuge tubes. These protocols are suitable for small amounts of silica gel-preserved plant tissue such as are commonly available from endangered plants. The critical step to remove polysaccharides was performing initial washes in either STE (0.25 M sucrose, 0.03 M Tris, 0.05 M EDTA) or HEPES (2% β-mercaptoethanol, 0.2% PVP, 0.1 M HEPES, pH 8.0) buffer. Precipitating the DNA at room temperature with isopropanol also aided in decreasing polysaccharide co-precipitation. Of the two protocols we present the STE/CTAB method has the advantages of being more cost-effective and avoiding the use of the hazardous chemical β-mercaptoethanol.

Microsatellite DNA analyses of a highly disjunct New Zealand tree reveal strong differentiation and imply a formerly more continuous distribution

Although New Zealand is a biodiversity hotspot, there has been little genetic investigation of why so many of its threatened and uncommon plants have naturally disjunct distributions. We investigated the small tree Pseudopanax ferox (Araliaceae), which has a widespread but highly disjunct lowland distribution within New Zealand. Genotyping of nuclear microsatellites and a chloroplast locus revealed pronounced genetic differentiation and four principal genetic clusters. Our results indicate that the disjunct distribution is a product of vicariance rather than long-distance dispersal. This highlights the need to preserve multiple populations when disjunct distributions are the result of vicariance, rather than focusing conservation efforts on a core area, in order to retain as much as possible of a species' evolutionary legacy and potential. Additionally, based on our genetic findings and the ecology of P. ferox, we hypothesize that it was more continuously distributed during the drier (but not maximally colder) interstadials of glacial periods and/or on the fertile soils available immediately postglacial. We further hypothesize that P. ferox belongs to a suite of species of drought-prone and/or fertile habitats whose distributions are actually restricted during warmer and wetter interglacial periods, despite being principally of the lowlands. Our genetic data for P. ferox are also the first consistent with the survival during the Last Glacial Maxima of a lowland tree at high latitudes in the south-eastern South Island.

Parallel polyploid speciation: distinct sympatric gene-pools of recurrently derived allo-octoploid Asplenium ferns

Although polyploidy is widespread, its significance to the generation of biodiversity remains unclear. Many polyploids have been derived recurrently. For a particular polyploid, gene-flow between the products of independent origin is typical where they come into contact. Here, we use AFLP DNA-fingerprinting and chloroplast DNA sequences to demonstrate parallel polyploid speciation within both of the ferns Asplenium cimmeriorum and A. gracillimum. Both of these taxa comprise at least two allopolyploids, recurrently derived from the same progenitor pair. Each of these allopolyploids remain genetically distinguishable even with extensive sympatry, and could therefore be considered distinct species. To our knowledge, parallel speciation on this scale amongst recurrent polyploids has not been previously reported. With their parallel origins, these 'evolutionary replicates' provide an unrivalled opportunity to investigate how the reproductive barriers and ecological differentiation necessary for speciation arise following polyploidy.

Tree endurance on the Tibetan Plateau marks the world's highest known tree line of the Last Glacial Maximum

Because of heterogeneous topographies, high-mountain areas could harbor a significant pool of cryptic forest refugia (glacial microrefugia unrecognized by palaeodata), which, as a result of poor accessibility, have been largely overlooked. The juniper forests of the southern Tibetan Plateau, with one of the highest tree lines worldwide, are ideal for assessing the potential of high-mountain areas to harbor glacial refugia. Genetic evidence for Last Glacial Maximum (LGM) endurance of these microrefugia is presented using paternally inherited chloroplast markers. Five-hundred and ninety individuals from 102 populations of the Juniperus tibetica complex were sequenced at three polymorphic chloroplast regions. Significant interpopulation differentiation and phylogeographic structure were detected (G(ST) = 0.49, N(ST) = 0.72, N(ST) > G(ST), P < 0.01), indicating limited among-population gene flow. Of 62 haplotypes recovered, 40 were restricted to single populations. These private haplotypes and overall degrees of diversity were evenly spread among plateau and edge populations, strongly supporting the existence of LGM microrefugia throughout the present distribution range, partly well above 3500 m. These results mark the highest LGM tree lines known, illustrating the potential significance of high-mountain areas for glacial refugia. Furthermore, as the close vicinity of orographic rear-edge and leading-edge populations potentially allows gene flow, surviving populations could preserve the complete spectrum of rear-edge and leading-edge adaptations.

Reconstructing the species phylogeny of Pseudopanax (Araliaceae), a genus of hybridising trees

Pseudopanax (Araliaceae) comprises 12 tree species of diverse morphology and ecology, and is endemic to New Zealand. It is notable for the hybridisation that occurs between P. crassifolius and P. lessonii, which have very different leaves and habits. To provide context for the study of this hybridisation and other investigations, we examined the phylogeny of Pseudopanax using chloroplast DNA sequences (c.5900 base-pairs) and AFLP DNA-fingerprinting. Both approaches resolve two principal groups within Pseudopanax--the Arboreus group and the Crassifolius+Lessonii union--but how they are related to other genera remains unclear. There is, nevertheless, little compelling evidence against the monophyly of Pseudopanax, making unnecessary the recent re-segregation of the Arboreus group as Neopanax. The chloroplast data provided minimal additional resolution, although the position of P. linearis was discordant compared to other data. Analyses of the AFLP data strongly recovered each species, aside from the morphologically heterogeneous P. colensoi, and the two mainland species (P. arboreus and P. crassifolius) that each contained a nested island-endemic (P. kermadecensis and P. chathamicus, respectively). However, relationships amongst species within the two principal groups were poorly resolved. An example was the uncertainty of whether P. crassifolius grouped with P. lessonii and its allies, or the morphologically similar species it had been previously placed with. This in turn raises the issue of how hybridisation might affect phylogenies and the ability to reconstruct them, even when using multiple, independent markers.

Conflict amongst chloroplast DNA sequences obscures the phylogeny of a group of Asplenium ferns

A previous study of the relationships amongst three subgroups of the Austral Asplenium ferns found conflicting signal between the two chloroplast loci investigated. Because organelle genomes like those of chloroplasts and mitochondria are thought to be non-recombining, with a single evolutionary history, we sequenced four additional chloroplast loci with the expectation that this would resolve these relationships. Instead, the conflict was only magnified. Although tree-building analyses favoured one of the three possible trees, one of the alternative trees actually had one more supporting site (six versus five) and received greater support in spectral and neighbor-net analyses. Simulations suggested that chance alone was unlikely to produce strong support for two of the possible trees and none for the third. Likelihood permutation tests indicated that the concatenated chloroplast sequence data appeared to have experienced recombination. However, recombination between the chloroplast genomes of different species would be highly atypical, and corollary supporting observations, like chloroplast heteroplasmy, are lacking. Wider taxon sampling clarified the composition of the Austral group, but the conflicting signal meant analyses (e.g., morphological evolution, biogeographic) conditional on a well-supported phylogeny could not be performed.