Intraspecific divergence and evolution of a life-history trade-off along a successional gradient in Hawaii's Metrosideros polymorpha

Foliar functional and genetic variation in a keystone Hawaiian tree species estimated through spectroscopy

Imaging spectroscopy has the potential to map closely related plant taxa at landscape scales. Although spectral investigations at the leaf and canopy levels have revealed relationships between phylogeny and reflectance, understanding how spectra differ across, and are inherited from, genotypes of a single species has received less attention. We used a common-garden population of four varieties of the keystone canopy tree, Metrosideros polymorpha, from Hawaii Island and four F1-hybrid genotypes derived from controlled crosses to determine if reflectance spectra discriminate sympatric, conspecific varieties of this species and their hybrids. With a single exception, pairwise comparisons of leaf reflectance patterns successfully distinguished varieties of M. polymorpha on Hawaii Island as well as populations of the same variety from different islands. Further, spectral variability within a single variety from Hawaii Island and the older island of Oahu was greater than that observed among the four varieties on Hawaii Island. F1 hybrids most frequently displayed leaf spectral patterns intermediate to those of their parent taxa. Spectral reflectance patterns distinguished each of two of the hybrid genotypes from one of their parent varieties, indicating that classifying hybrids may be possible, particularly if sample sizes are increased. This work quantifies a baseline in spectral variability for an endemic Hawaiian tree species and advances the use of imaging spectroscopy in biodiversity studies at the genetic level.

An optimized RNA extraction method for diverse leaves of Hawaiian Metrosideros, a hypervariable tree species complex

Premise

The isolation of RNA from trees is challenging due to the interference of polyphenols and polysaccharides with downstream processes. Furthermore, many RNA extraction protocols are time consuming and involve hazardous chemicals. To address these issues, we aimed to develop a safe protocol for high-quality RNA extraction from diverse Metrosideros taxa representing a broad range of leaf toughness, pubescence, and secondary metabolites.

Methods and Results

We tested popular RNA isolation kits and protocols that were effective on other recalcitrant trees, including a broad range of optimization and purification steps. We optimized a protocol involving two silica-membrane column-based kits that yielded high-quantity RNA with an RNA integrity number >7 and without DNA contamination. All RNA samples were used successfully in a follow-on RNA-Seq experiment.

Conclusions

We present an optimized high-throughput RNA extraction protocol that yielded high-quality and high-quantity RNA from three contrasting leaf phenotypes within a hyperdiverse woody species complex.

Early detection of a tree pathogen using airborne remote sensing

Native forests of Hawai'i Island are experiencing an ecological crisis in the form of Rapid 'Ōhi'a Death (ROD), a recently characterized disease caused by two fungal pathogens in the genus Ceratocystis. Since approximately 2010, this disease has caused extensive mortality of Hawai'i's keystone endemic tree, known as 'ōhi'a (Metrosideros polymorpha). Visible symptoms of ROD include rapid browning of canopy leaves, followed by death of the tree within weeks. This quick progression leading to tree mortality makes early detection critical to understanding where the disease will move at a timescale feasible for controlling the disease. We used repeat laser-guided imaging spectroscopy (LGIS) of forests on Hawai'i Island collected by the Global Airborne Observatory (GAO) in 2018 and 2019 to derive maps of foliar trait indices previously found to be important in distinguishing between ROD-infected and healthy 'ōhi'a canopies. Data from these maps were used to develop a prognostic indicator of tree stress prior to the visible onset of browning. We identified canopies that were green in 2018, but became brown in 2019 (defined as "to become brown"; TBB), and a corresponding set of canopies that remained green. The data mapped in 2018 showed separability of foliar trait indices between TBB and green 'ōhi'a, indicating early detection of canopy stress prior to the onset of ROD. Overall, a combination of linear and non-linear analyses revealed canopy water content (CWC), foliar tannins, leaf mass per area (LMA), phenols, cellulose, and non-structural carbohydrates (NSC) are primary drivers of the prognostic spectral capability which collectively result in strong consistent changes in leaf spectral reflectance in the near-infrared (700-1300 nm) and shortwave-infrared regions (1300-2500 nm). Results provide insight into the underlying foliar traits that are indicative of physiological responses of M. polymorpha trees infected with Ceratocycstis and suggest that imaging spectroscopy is an effective tool for identifying trees likely to succumb to ROD prior to the onset of visible symptoms.

Demography and selection analysis of the incipient adaptive radiation of a Hawaiian woody species

Ecological divergence in a species provides a valuable opportunity to study the early stages of speciation. We focused on Metrosideros polymorpha, a unique example of the incipient radiation of woody species, to examine how an ecological divergence continues in the face of gene flow. We analyzed the whole genomes of 70 plants collected throughout the island of Hawaii, which is the youngest island with the highest altitude in the archipelago and encompasses a wide range of environments. The continuous M. polymorpha forest stands on the island of Hawaii were differentiated into three genetic clusters, each of which grows in a distinctive environment and includes substantial genetic and phenotypic diversity. The three genetic clusters showed signatures of selection in genomic regions encompassing genes relevant to environmental adaptations, including genes associated with light utilization, oxidative stress, and leaf senescence, which are likely associated with the ecological differentiation of the species. Our demographic modeling suggested that the glaberrima cluster in wet environments maintained a relatively large population size and two clusters split: polymorpha in the subalpine zone and incana in dry and hot conditions. This ecological divergence possibly began before the species colonized the island of Hawaii. Interestingly, the three clusters recovered genetic connectivity coincidentally with a recent population bottleneck, in line with the weak reproductive isolation observed in the species. This study highlights that the degree of genetic differentiation between ecologically-diverged populations can vary depending on the strength of natural selection in the very early phases of speciation.

Knowledge about how genetic barriers are formed between populations in distinct environments is valuable to understand the processes of speciation and conserve biodiversity. Metrosideros polymorpha, an endemic woody species in the Hawaiian Islands, is a good system to study developing genetic barriers in a species, because it colonized the diverse environments and diversified the morphology for a relatively short period of time. We analyzed the genomes of 70 M. polymorpha plants from a broad range of environments on the island of Hawaii to infer the current and past genetic barriers among them. Currently, M. polymorpha plants growing in different environments have substantially different genomes, especially at the genomic regions with genes putatively controlling physiology to fit in distinct environment. However, in its history, they had hybridized with one another, possibly because plants formerly growing in different environments came into close contact due to the climate changes. It is suggested that genetic barriers can easily strengthen or weaken depending on environments splitting the ecology of a species before reproductive isolation becomes complete.

Ancestral polymorphisms shape the adaptive radiation of Metrosideros across the Hawaiian Islands

Some of the most spectacular adaptive radiations begin with founder populations on remote islands. How genetically limited founder populations give rise to the striking phenotypic and ecological diversity characteristic of adaptive radiations is a paradox of evolutionary biology. We conducted an evolutionary genomics analysis of genus Metrosideros, a landscape-dominant, incipient adaptive radiation of woody plants that spans a striking range of phenotypes and environments across the Hawaiian Islands. Using nanopore-sequencing, we created a chromosome-level genome assembly for Metrosideros polymorpha var. incana and analyzed whole-genome sequences of 131 individuals from 11 taxa sampled across the islands. Demographic modeling and population genomics analyses suggested that Hawaiian Metrosideros originated from a single colonization event and subsequently spread across the archipelago following the formation of new islands. The evolutionary history of Hawaiian Metrosideros shows evidence of extensive reticulation associated with significant sharing of ancestral variation between taxa and secondarily with admixture. Taking advantage of the highly contiguous genome assembly, we investigated the genomic architecture underlying the adaptive radiation and discovered that divergent selection drove the formation of differentiation outliers in paired taxa representing early stages of speciation/divergence. Analysis of the evolutionary origins of the outlier single nucleotide polymorphisms (SNPs) showed enrichment for ancestral variations under divergent selection. Our findings suggest that Hawaiian Metrosideros possesses an unexpectedly rich pool of ancestral genetic variation, and the reassortment of these variations has fueled the island adaptive radiation.

Isolation of Metrosideros ('Ohi'a) Taxa on O'ahu Increases with Elevation and Extreme Environments

Species radiations should be facilitated by short generation times and limited dispersal among discontinuous populations. Hawaii's hyper-diverse, landscape-dominant tree, Metrosideros, is unique among the islands' radiations for its massive populations that occur continuously over space and time within islands, its exceptional capacity for gene flow by both pollen and seed, and its extended life span (ca. >650 years). Metrosideros shows the greatest phenotypic and microsatellite DNA diversity on O'ahu, where taxa occur in tight sympatry or parapatry in mesic and montane wet forest on 2 volcanoes. We document the nonrandom distributions of 12 taxa (including unnamed morphotypes) along elevation gradients, measure phenotypes of ~6-year-old common-garden plants of 8 taxa to verify heritability of phenotypes, and examine genotypes of 476 wild adults at 9 microsatellite loci to compare the strengths of isolation across taxa, volcanoes, and distance. All 8 taxa retained their diagnostic phenotypes in the common garden. Populations were isolated by taxon to a range of degrees (pairwise FST between taxa: 0.004-0.267), and there was no pattern of isolation by distance or by elevation; however, significant isolation between volcanoes was observed within monotypic species, suggesting limited gene flow between volcanoes. Among the infraspecific taxa of Metrosideros polymorpha, genetic diversity and isolation significantly decreased and increased, respectively, with elevation. Overall, 5 of the 6 most isolated taxa were associated with highest elevations or otherwise extreme environments. These findings suggest a principal role for selection in the origin and maintenance of the exceptional diversity that occurs within continuous Metrosideros stands on O'ahu.

Divergent Selection and Primary Gene Flow Shape Incipient Speciation of a Riparian Tree on Hawaii Island

A long-standing goal of evolutionary biology is to understand the mechanisms underlying the formation of species. Of particular interest is whether or not speciation can occur in the presence of gene flow and without a period of physical isolation. Here, we investigated this process within Hawaiian Metrosideros, a hypervariable and highly dispersible woody species complex that dominates the Hawaiian Islands in continuous stands. Specifically, we investigated the origin of Metrosideros polymorpha var. newellii (newellii), a riparian ecotype endemic to Hawaii Island that is purportedly derived from the archipelago-wide M. polymorpha var. glaberrima (glaberrima). Disruptive selection across a sharp forest-riparian ecotone contributes to the isolation of these varieties and is a likely driver of newellii's origin. We examined genome-wide variation of 42 trees from Hawaii Island and older islands. Results revealed a split between glaberrima and newellii within the past 0.3-1.2 My. Admixture was extensive between lineages within Hawaii Island and between islands, but introgression from populations on older islands (i.e., secondary gene flow) did not appear to contribute to the emergence of newellii. In contrast, recurrent gene flow (i.e., primary gene flow) between glaberrima and newellii contributed to the formation of genomic islands of elevated absolute and relative divergence. These regions were enriched for genes with regulatory functions as well as for signals of positive selection, especially in newellii, consistent with divergent selection underlying their formation. In sum, our results support riparian newellii as a rare case of incipient ecological speciation with primary gene flow in trees.

Comparing Adaptive Radiations Across Space, Time, and Taxa

Adaptive radiation plays a fundamental role in our understanding of the evolutionary process. However, the concept has provoked strong and differing opinions concerning its definition and nature among researchers studying a wide diversity of systems. Here, we take a broad view of what constitutes an adaptive radiation, and seek to find commonalities among disparate examples, ranging from plants to invertebrate and vertebrate animals, and remote islands to lakes and continents, to better understand processes shared across adaptive radiations. We surveyed many groups to evaluate factors considered important in a large variety of species radiations. In each of these studies, ecological opportunity of some form is identified as a prerequisite for adaptive radiation. However, evolvability, which can be enhanced by hybridization between distantly related species, may play a role in seeding entire radiations. Within radiations, the processes that lead to speciation depend largely on (1) whether the primary drivers of ecological shifts are (a) external to the membership of the radiation itself (mostly divergent or disruptive ecological selection) or (b) due to competition within the radiation membership (interactions among members) subsequent to reproductive isolation in similar environments, and (2) the extent and timing of admixture. These differences translate into different patterns of species accumulation and subsequent patterns of diversity across an adaptive radiation. Adaptive radiations occur in an extraordinary diversity of different ways, and continue to provide rich data for a better understanding of the diversification of life.

Varieties of the highly dispersible and hypervariable tree, Metrosideros polymorpha, differ in response to mechanical stress and light across a sharp ecotone

PREMISE

The drivers of isolation between sympatric populations of long-lived and highly dispersible conspecific plants are not well understood. In the Hawaiian Islands, the landscape-dominant tree, Metrosideros polymorpha, displays extraordinary phenotypic differences among sympatric varieties despite high dispersibility of its pollen and seeds, thereby presenting a unique opportunity to investigate how disruptive selection alone can maintain incipient forms. Stenophyllous M. polymorpha var. newellii is a recently evolved tree endemic to the waterways of eastern Hawai'i Island that shows striking neutral genetic differentiation from its ancestor, wet-forest M. polymorpha var. glaberrima, despite sympatry of these forms. We looked for evidence for, and drivers of, differential local adaptation of these varieties across the range of M. polymorpha var. newellii.

METHODS

For paired populations of these varieties, we compared seedling performance under contrasting light conditions and a strong water current characteristic of the riparian zone. We also conducted a reciprocal transplant experiment and contrasted adult leaf anatomy.

RESULTS

Results suggest that the riparian zone is harsh and that selection involving the mechanical stress of rushing water, and secondarily, light, led to significant reciprocal immigrant inviability in adjacent forest and riparian environments. The strongest adaptive divergence between varieties was seen in leaves and seedlings from the site with the sharpest ecotone, coincident with the strongest genetic isolation of M. polymorpha var. newellii observed previously.

CONCLUSIONS

These findings suggest that disruptive selection across a sharp ecotone contributes to the maintenance of an incipient riparian ecotype from within a continuous population of a long-lived and highly dispersible tree species.

Incipient ecological speciation between successional varieties of a dominant tree involves intrinsic postzygotic isolating barriers

Whereas disruptive selection imposed by heterogeneous environments can lead to the evolution of extrinsic isolating barriers between diverging populations, the evolution of intrinsic postzygotic barriers through divergent selection is less certain. Long-lived species such as trees may be especially slow to evolve intrinsic isolating barriers. We examined postpollination reproductive isolating barriers below the species boundary, in an ephemeral hybrid zone between two successional varieties of the landscape-dominant Hawaiian tree, Metrosideros polymorpha, on volcanically active Hawai'i Island. These archipelago-wide sympatric varieties show the weakest neutral genetic divergence of any taxon pair on Hawai'i Island but significant morphological and ecological differentiation consistent with adaptation to new and old lava flows. Cross-fertility between varieties was high and included heterosis of F1 hybrids at the seed germination stage, consistent with a substantial genetic load apparent within varieties through low self-fertility and a lack of self-pollen discrimination. However, a partial, but significant, barrier was observed in the form of reduced female and male fertility of hybrids, especially backcross hybrids, consistent with the accumulation of genetic incompatibilities between varieties. These results suggest that partial intrinsic postzygotic barriers can arise through disruptive selection acting on large, hybridizing populations of a long-lived species.

Genetic analysis of an ephemeral intraspecific hybrid zone in the hypervariable tree, Metrosideros polymorpha, on Hawai'i Island

Intraspecific hybrid zones involving long-lived woody species are rare and can provide insights into the genetic basis of early-diverging traits in speciation. Within the landscape-dominant Hawaiian tree, Metrosideros polymorpha, are morphologically distinct successional varieties, incana and glaberrima, that dominate new and old lava flows, respectively, below 1200 me on volcanically active Hawai'i Island, with var. glaberrima also extending to higher elevations and bogs. Here, we use morphological measurements on 86 adult trees to document the presence of an incana-glaberrima hybrid zone on the 1855 Mauna Loa lava flow on east Hawai'i Island and parent-offspring analysis of 1311 greenhouse seedlings from 71 crosses involving 72 adults to estimate heritabilities and genetic correlations among vegetative traits. Both the variation in adult leaf pubescence at the site and the consistency between adult and offspring phenotypes suggest the presence of two hybrid classes, F1s and var. incana backcrosses, as would be expected on a relatively young lava flow. Nine nuclear microsatellite loci failed to distinguish parental and hybrid genotypes. All four leaf traits examined showed an additive genetic basis with moderate to strong heritabilities, and genetic correlations were stronger for the more range-restricted var. incana. The differences between varieties in trait values, heritabilities and genetic correlations, coupled with high genetic variation within but low genetic variation between varieties, are consistent with a multi-million-year history of alternating periods of disruptive selection in contrasting environments and admixture in ephemeral hybrid zones. Finally, the contrasting genetic architectures suggest different evolutionary trajectories of leaf traits in these forms.

The transcriptional landscape of insect galls: psyllid (Hemiptera) gall formation in Hawaiian Metrosideros polymorpha (Myrtaceae)

BACKGROUND

Recent studies show that galling Hymenoptera and Diptera are able to synthesize the plant hormone indole-3-acetic acid (auxin) from tryptophan and that plant response to insect-produced auxin is implicated in gall formation. We examined the leaf transcriptome of galled and ungalled leaves of individuals of the Hawaiian endemic plant Metrosideros polymorpha (Myrtaceae) subject to infestation by psyllid (Hemiptera) gall-makers in the genus Trioza (Triozidae).

RESULTS

Transcript libraries were sequenced using Illumina technology and the reads assembled de novo into contigs. Functional identification of contigs followed a two-step procedure, first identifying contigs by comparison to the completely sequenced genome of the related Eucalyptus, followed by identifying the equivalent Arabidopsis gene using a pre-computed mapping between Eucalyptus and Arabidopsis genes. This allowed us to use the rich functional annotation of the Arabidopsis genome to assess the transcriptional landscape of galling in Metrosideros. Comparing galled and ungalled leaves, we find a highly significant enrichment of expressed genes with a gene ontology (GO) annotation to auxin response in the former. One gene consistently expressed in all galled trees examined but not detected in any libraries from ungalled leaves was the Metrosideros version of SMALL AUXIN UPREGULATED (SAUR) 67 which appears to be a marker for leaf-galling in Metrosideros.

CONCLUSIONS

We conclude that an auxin response is involved in galling by Metrosideros psyllids. The possibility should therefore be considered that psyllids (like other insects examined) are able to synthesize auxin.

Incipient radiation within the dominant Hawaiian tree Metrosideros polymorpha

Although trees comprise a primary component of terrestrial species richness, the drivers and temporal scale of divergence in trees remain poorly understood. We examined the landscape-dominant tree, Metrosideros polymorpha, for variation at nine microsatellite loci across 23 populations on young Hawai'i Island, sampling each of the island's five varieties throughout its full geographic range. For four varieties, principal coordinate analysis revealed strong clustering of populations by variety across the 10 430 km(2) island, indicating partitioning of the species into multiple evolutionarily significant units. The single island-endemic form, riparian var. newellii, showed especially strong differentiation from other varieties despite occurring in sympatry with other varieties and likely evolved from a bog form on the oldest volcano, Kohala, within the past 500 000 years. Along with comparable riparian forms on other Pacific Islands, var. newellii appears to represent parallel incipient ecological speciation within Metrosideros. Greater genetic distance among the more common varieties on the oldest volcano and an inverse relationship between allelic diversity and substrate age appear consistent with colonization of Hawai'i Island by older, partially diverged varieties followed by increased hybridization among varieties on younger volcanoes. This study demonstrates that broad population-level sampling is required to uncover patterns of diversification within a ubiquitous and long-lived tree species. Hawaiian Metrosideros appears to be a case of incipient radiation in trees and thus should be useful for studies of divergence and the evolution of reproductive isolating barriers at the early stages of speciation.