Structural insights into the Smirnoff-Wheeler pathway for vitamin C production in the Amazon fruit camu-camu.

l-Ascorbic acid (AsA, vitamin C) is a pivotal dietary nutrient with multifaceted importance in living organisms. In plants, the Smirnoff-Wheeler pathway is the primary route for AsA biosynthesis, and understanding the mechanistic details behind its component enzymes has implications for plant biology, nutritional science, and biotechnology. As part of an initiative to determine the structures of all six core enzymes of the pathway, the present study focuses on three of them in the model species Myrciaria dubia (camu-camu): GDP-d-mannose 3',5'-epimerase (GME), l-galactose dehydrogenase (l-GalDH), and l-galactono-1,4-lactone dehydrogenase (l-GalLDH). We provide insights into substrate and cofactor binding and the conformational changes they induce. The MdGME structure reveals a distorted substrate in the active site, pertinent to the catalytic mechanism. Mdl-GalDH shows that the way in which NAD+ association affects loop structure over the active site is not conserved when compared with its homologue in spinach. Finally, the structure of Mdl-GalLDH is described for the first time. This allows for the rationalization of previously identified residues which play important roles in the active site or in the formation of the covalent bond with FAD. In conclusion, this study enhances our understanding of AsA biosynthesis in plants, and the information provided should prove useful for biotechnological applications.

The In Planta Gene Expression of Austropuccinia psidii in Resistant and Susceptible Eucalyptus grandis.

Austropuccinia psidii, commonly known as myrtle rust, is an obligate, biotrophic rust pathogen that causes rust disease in a broad host range of Myrtaceae species. Eucalyptus grandis, a widely cultivated hardwood Myrtaceae species, is susceptible to A. psidii infection, with this pathogen threatening both their natural range and various forest plantations across the world. This study aimed to investigate the A. psidii transcriptomic responses in resistant and susceptible E. grandis at four time points. RNA-seq reads were mapped to the A. psidii reference genome to quantify expressed genes at 12 h postinoculation and 1, 2, and 5 days postinoculation (dpi). A total of eight hundred and ninety expressed genes were found, of which 43 were candidate effector protein genes. These included rust transferred protein 1 (RTP1), expressed in susceptible hosts at 5 dpi, and a hydrolase protein gene expressed in both resistant and susceptible hosts over time. Functional categorization of expressed genes revealed processes enriched in susceptible hosts, including malate metabolic and malate dehydrogenase activity, implicating oxalic acid in disease susceptibility. These results highlight putative virulence or pathogenicity mechanisms employed by A. psidii to cause disease, and they provide the first insight into the molecular responses of A. psidii in E. grandis over time.

Genetic associations with resistance to Meloidogyne enterolobii in guava (Psidium sp.) using cross-genera SNPs and comparative genomics to Eucalyptus highlight evolutionary conservation across the Myrtaceae.

Tropical fruit tree species constitute a yet untapped supply of outstanding diversity of taste and nutritional value, barely developed from the genetics standpoint, with scarce or no genomic resources to tackle the challenges arising in modern breeding practice. We generated a de novo genome assembly of the Psidium guajava, the super fruit "apple of the tropics", and successfully transferred 14,268 SNP probesets from Eucalyptus to Psidium at the nucleotide level, to detect genomic loci linked to resistance to the root knot nematode (RKN) Meloidogyne enterolobii derived from the wild relative P. guineense. Significantly associated loci with resistance across alternative analytical frameworks, were detected at two SNPs on chromosome 3 in a pseudo-assembly of Psidium guajava genome built using a syntenic path approach with the Eucalyptus grandis genome to determine the order and orientation of the contigs. The P. guineense-derived resistance response to RKN and disease onset is conceivably triggered by mineral nutrients and phytohormone homeostasis or signaling with the involvement of the miRNA pathway. Hotspots of mapped resistance quantitative trait loci and functional annotation in the same genomic region of Eucalyptus provide further indirect support to our results, highlighting the evolutionary conservation of genomes across genera of Myrtaceae in the adaptation to pathogens. Marker assisted introgression of the resistance loci mapped should accelerate the development of improved guava cultivars and hybrid rootstocks.

Rhodomyrtus tomentosa (Ait.) Hassk fruit phenolic-rich extract mitigates intestinal barrier dysfunction and inflammation in mice.

In this study, the mitigative effect of Rhodomyrtus tomentosa (Ait.) Hassk fruit extract rich in phenolic compounds (RTE) on high fat diet (HFD)-induced intestinal barrier dysfunction of mice and the underlying mechanism were explored. The results revealed that RTE supplementation obviously improved gut microbiota dysbiosis induced by HFD, which was evidenced by elevated alpha diversity, suppressed Firmicutes/Bacteroidetes ratio, enriched short-chain fatty acid-producing bacteria (Odoribacter, Parabacteroides, Blautia and Akkermansia), and depleted harmful bacteria (Helicobacter, norank_f_ Desulfovibrionaceae and Mucispirillum). RTE intervention mitigated intestinal barrier dysfunction and inflammation by elevating tight junction proteins expression levels and decreasing proinflammatory cytokines levels. Furthermore, RTE administration inhibited the HFD-induced trigger of the lipopolysaccharide-toll-like receptor 4-nuclear factor kappa-B (LPS-TLR4-NF-κB) pathway in colonic tissue. Therefore, RTE supplementation may be an effective way to protect the intestinal tract in HFD-induced obese individuals.

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.

Diversity, phylogeny and evolution of the rapidly evolving genus Psidium L. (Myrtaceae, Myrteae).

BACKGROUND AND AIMS: Psidium is the fourthth largest genus of Myrtaceae in the Neotropics. Psidium guajava is widely cultivated in the tropics for its edible fruit. It is commercially under threat due to the disease guava decline. Psidium cattleyanum is one of the 100 most invasive organisms in the world. Knowledge of the phylogenetic relationships within Psidium is poor. We aim to provide a review of the biology, morphology and ecology of Psidium, a phylogenetic tree, an infrageneric classification and a list of species. METHODS: Morphological and geographic data were obtained by studying Psidium in herbaria and in the field between 1988 and 2020. Forty-six herbaria were visited personally. A database of approx. 6000 specimens was constructed, and the literature was reviewed. Thirty species (about a third of the species in the genus) were sampled for molecular phylogenetic inference. Two chloroplast (psbA-trnH and ndhF) and two nuclear (external transcribed spacer and internal transcribed spacer) regions were targeted. Phylogenetic trees were constructed using maximum likelihood (ML; RaxML) and Bayesian inference (BI; MrBayes). KEY RESULTS: Psidium is a monophyletic genus with four major clades recognized as sections. Section Psidium (ten species), to which P. guajava belongs, is sister to the rest of the genus; it is widespread across the Neotropics. Section Obversifolia (six species; restricted to the Brazilian Atlantic Forest), which includes P. cattleyanum, is sister to the innermost clade composed of sister sections Apertiflora (31 species; widespread but most diverse in the Brazilian Atlantic Forest) + Mitranthes (26 species; widespread in dry forests and probably diverse in the Caribbean). Characters associated with diversification within Psidium are discussed. CONCLUSIONS: Research on pre-foliation, colleters, leaf anatomy, leaf physiology, staminal development, placentation and germination associated with the anatomy of the opercular plug is desirable. Studies are biased towards sections Psidium and Obversifolia, with other sections poorly known.

Pests, diseases, and aridity have shaped the genome of Corymbia citriodora.

Corymbia citriodora is a member of the predominantly Southern Hemisphere Myrtaceae family, which includes the eucalypts (Eucalyptus, Corymbia and Angophora; ~800 species). Corymbia is grown for timber, pulp and paper, and essential oils in Australia, South Africa, Asia, and Brazil, maintaining a high-growth rate under marginal conditions due to drought, poor-quality soil, and biotic stresses. To dissect the genetic basis of these desirable traits, we sequenced and assembled the 408 Mb genome of Corymbia citriodora, anchored into eleven chromosomes. Comparative analysis with Eucalyptus grandis reveals high synteny, although the two diverged approximately 60 million years ago and have different genome sizes (408 vs 641 Mb), with few large intra-chromosomal rearrangements. C. citriodora shares an ancient whole-genome duplication event with E. grandis but has undergone tandem gene family expansions related to terpene biosynthesis, innate pathogen resistance, and leaf wax formation, enabling their successful adaptation to biotic/abiotic stresses and arid conditions of the Australian continent.

Development of microsatellite markers for Calyptranthes clusiifolia (Myrtaceae).

The endemic tree species Calyptranthes clusiifolia (Myrtaceae) plays a relevant ecological role in the forest fragments where it has a common occurrence. In this study, we reported the development of microsatellite markers for C. clusiifolia what will allow a better understanding of the relationship between the forest fragmentation process and the genetic structure and diversity of tree populations. Seven microsatellite markers were developed using an enriched genomic library and characterized in 30 individuals (from three populations). These seven loci were polymorphic and resulted in a total of 23 alleles. The expected heterozygosity (HE) varied from 0.14 (Caly 06) to 0.73 (Caly 22). Linkage disequilibrium between the loci (p > 0.0007) pairs was not detected. The parentage exclusion power of the first (Pe-1) and the second (Pe-2) parents were 0.6099 and 0.8548, respectively. The microsatellite markers developed are indicated for future studies of the genetic diversity in natural populations of C. clusiifolia.

Tracing foreign sequences in plant transcriptomes and genomes using OCT4, a POU domain protein.

Contaminations in sequencing data, especially in reference genomes, lead to inevitable errors in downstream analyses. Similarly, presence of contaminants in transcriptomes, misrepresents the molecular basis of various interactions. In this study, we report the presence of a large number of plant transcriptomes contaminated with RNAs encoding POU domain proteins; a family of proteins that has not been reported in plants and fungi. Besides, our findings illustrated that there are four POU domain protein-coding sequences in the reference genome of Rhodamnia argentea. It turned out that the existing foreign fragments are related to arthropods that are considered as plant pests. We also identified two contaminated draft genomes, Humulus lupulus and Cannabis sativa that contained complete rDNA sequences originating from Tetranychus species. As a result, careful screening of sequencing data before releasing them in public databases or checking existing genomes for possible contaminations is recommended.

Phylogeny and biogeography of Myrcia sect. Aguava (Myrtaceae, Myrteae) based on phylogenomic and Sanger data provide evidence for a Cerrado origin and geographically structured clades.

Myrcia is one of the largest exclusively Neotropical angiosperm genera, including ca. 800 species divided into nine sections. Myrcia sect. Aguava is one of most complex sections of Myrcia due to high morphological variation and wide distribution range of some species, including M. guianensis, with distribution throughout South America and a complex taxonomic history. We used complete plastid DNA sequences data generated using next-generation sequencing of 45 terminals, mostly from Myrcia sect. Aguava. These data were combined with five target DNA regions (ITS, psbA-trnH, trnL-trnF, trnQ-rps16, ndhF) of additional terminals to increase taxonomic coverage. Phylogenetic analyses were conducted using a maximum likelihood approach, and divergence times and ancestral range distributions were estimated. Myrcia sect. Aguava is monophyletic and exclusively comprises species with trilocular ovaries but has no relationship with other groups within Myrcia that possess trilocular ovaries. Three main lineages that correspond to geographical distribution are recognized within Myrcia sect. Aguava. Multiple accessions reveal a non-monophyletic Myrcia guianensis and stress the biogeographical structure inside the group. Myrcia sect. Aguava had a probable mid-Miocene origin in the Cerrado, but lineages that persisted there diversified only more recently, when the present-day vegetation started to stabilize. Posterior migrations to Atlantic Forest, Amazon and Caribbean occurred at the end of Miocene, evidencing transitions from open and dry to forested and more humid areas that are less frequent in the Neotropics. Overall, it is observed that related lineages remained in ecologically similar environments. Future perspectives on Myrcia and Myrteae in the phylogenomic era are also discussed.

Monoterpene synthases responsible for the terpene profile of anther glands in Eucalyptus polybractea R.T. Baker (Myrtaceae).

Research on terpene biosynthesis in the genus Eucalyptus (Myrtaceae) is poorly developed, but recently large numbers of terpene synthase (TPS) genes have been identified. Few of these have been characterized or their expression localized to specific tissues. A prime candidate for detailed examination of TPS gene expression is the bisexual eucalypt flower-composed of male and female reproductive organs, and vegetative tissues that may express different TPS genes. We aimed to characterize and compare the terpene profile and TPS genes expressed in anthers and gynoecia in the high oil-yielding Eucalyptus polybractea R.T. Baker. We hypothesized that gynoecia will produce greater amounts of defensive terpenes, whereas anthers will have a terpene profile that is biased towards a role in pollination. Microscopy of isolated anthers showed them to possess a single, prominent oil gland. Chemical analysis of whole floral structures at different stages of development showed total oil per unit dry mass increased as flower buds expanded, with highest concentrations in mature flower buds just prior to flower opening. The oil profile of gynoecia was dominated by the monoterpene 1,8-cineole, whereas that of isolated anthers were enriched with the monoterpene α-pinene. Through transcriptomic analysis and recombinant protein expression, we were able to identify monoterpene synthases responsible for the different profiles. Synthases for α-pinene and 1,8-cineole were expressed in each tissue type, but the relative expression of the former was higher in anthers. Sequence comparison and site-directed mutagenesis of the α-pinene synthase allowed us to identify amino acids that influence the α-pinene to ß-pinene ratio of the product profile. We suggest the terpene constituents of anthers may have multiple roles including attracting pollinators through emission of volatile α-pinene, deterrence of palynivores through emission of volatile 1,8-cineole and adhesion of pollen to pollinators via the release of sticky α-pinene onto the anther surface.

Looks can be deceiving: speciation dynamics of co-distributed Angophora (Myrtaceae) species in a varying landscape.

Understanding the mechanisms underlying species divergence remains a central goal in evolutionary biology. Landscape genetics can be a powerful tool for examining evolutionary processes. We used genome-wide scans to genotype samples from populations of eight Angophora species. Angophora is a small genus within the eucalypts comprising common and rare species in a heterogeneous landscape, making it an appropriate group to study speciation. We found A. hispida was highly differentiated from the other species. Two subspecies of A. costata (subsp. costata and subsp. euryphylla) formed a group, while the third (subsp. leiocarpa, which is only distinguished by its smooth fruits and provenance) was supported as a distinct pseudocryptic species. Other species that are morphologically distinct could not be genetically differentiated (e.g., A. floribunda and A. subvelutina). Distribution and genetic differentiation within Angophora were strongly influenced by temperature and humidity, as well as biogeographic barriers, particularly rivers and higher elevation regions. While extensive introgression was found between many populations of some species (e.g., A. bakeri and A. floribunda), others only hybridized at certain locations. Overall, our findings suggest multiple mechanisms drove evolutionary diversification in Angophora and highlight how genome-wide analyses of related species in a diverse landscape can provide insights into speciation.

Extensive Genetic Connectivity and Historical Persistence Are Features of Two Widespread Tree Species in the Ancient Pilbara Region of Western Australia.

Phylogeographic studies can be used as a tool to understand the evolutionary history of a landscape, including the major drivers of species distributions and diversity. Extensive research has been conducted on phylogeographic patterns of species found in northern hemisphere landscapes that were affected by glaciations, yet the body of literature for older, unaffected landscapes is still underrepresented. The Pilbara region of north-western Australia is an ancient and vast landscape that is topographically complex, consisting of plateaus, gorges, valleys, and ranges, and experiences extreme meteorological phenomena including seasonal cyclonic activity. These features are expected to influence patterns of genetic structuring throughout the landscape either by promoting or restricting the movement of pollen and seed. Whilst a growing body of literature exists for the fauna endemic to this region, less is known about the forces shaping the evolution of plant taxa. In this study we investigate the phylogeography of two iconic Pilbara tree species, the Hamersley Bloodwood (Corymbia hamersleyana) and Western Gidgee (Acacia pruinocarpa), by assessing patterns of variation and structure in several chloroplast DNA regions and nuclear microsatellite loci developed for each species. Gene flow was found to be extensive in both taxa and there was evidence of long-distance seed dispersal across the region (pollen to seed ratios of 6.67 and 2.96 for C. hamersleyana and A. pruinocarpa, respectively), which may result from flooding and strong wind gusts associated with extreme cyclonic activity. Both species possessed high levels of cpDNA genetic diversity in comparison to those from formerly glaciated landscapes (C. hamersleyana = 14 haplotypes, A. pruinocarpa = 37 haplotypes) and showed evidence of deep lineage diversification occurring from the late Miocene, a time of intensifying aridity in this landscape that appears to be a critical driver of evolution in Pilbara taxa. In contrast to another study, we did not find evidence for topographic features acting as refugia for the widely sampled C. hamersleyana.

Metabolomics of Myrcia bella Populations in Brazilian Savanna Reveals Strong Influence of Environmental Factors on Its Specialized Metabolism.

Environmental conditions influence specialized plant metabolism. However, many studies aiming to understand these modulations have been conducted with model plants and/or under controlled conditions, thus not reflecting the complex interaction between plants and environment. To fully grasp these interactions, we investigated the specialized metabolism and genetic diversity of a native plant in its natural environment. We chose Myrcia bella due to its medicinal interest and occurrence in Brazilian savanna regions with diverse climate and soil conditions. An LC-HRMS-based metabolomics approach was applied to analyze 271 samples harvested across seven regions during the dry and rainy season. Genetic diversity was assessed in a subset of 40 samples using amplified fragment length polymorphism. Meteorological factors including rainfall, temperature, radiation, humidity, and soil nutrient and mineral composition were recorded in each region and correlated with chemical variation through multivariate analysis (MVDA). Marker compounds were selected using a statistically informed molecular network and annotated by dereplication against an in silico database of natural products. The integrated results evidenced different chemotypes, with variation in flavonoid and tannin content mainly linked to soil conditions. Different levels of genetic diversity and distance of populations were found to be correlated with the identified chemotypes. These observations and the proposed analytical workflow contribute to the global understanding of the impact of abiotic factors and genotype on the accumulation of given metabolites and, therefore, could be valuable to guide further medicinal exploration of native species.

Sympatric speciation in mountain roses (Metrosideros) on an oceanic island.

Shifts in flowering time have the potential to act as strong prezygotic reproductive barriers in plants. We investigate the role of flowering time divergence in two species of mountain rose (Metrosideros) endemic to Lord Howe Island, Australia, a minute and isolated island in the Tasman Sea. Metrosideros nervulosa and M. sclerocarpa are sister species and have divergent ecological niches on the island but grow sympatrically for much of their range, and likely speciated in situ on the island. We used flowering time and population genomic analyses of population structure and selection, to investigate their evolution, with a particular focus on the role of flowering time in their speciation. Population structure analyses showed the species are highly differentiated and appear to be in the very late stages of speciation. We found flowering times of the species to be significantly displaced, with M. sclerocarpa flowering 53 days later than M. nervulosa. Furthermore, the analyses of selection showed that flowering time genes are under selection between the species. Thus, prezygotic reproductive isolation is mediated by flowering time shifts in the species, and likely evolved under selection, to drive the completion of speciation within a small geographical area. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

The complete chloroplast genome and characteristics analysis of Callistemon rigidus R.Br.

Callistemon rigidus R.Br. one of the traditional Chinese medicinal plants, is acrid-flavored and mild-natured, with the prominent effects reducing swelling, resolving phlegm, and dispelling rheumatism. Clinically, it has been commonly used to treat cold, cough and asthma, pain and swelling from impact injuries, eczema, rheumatic arthralgia. The chloroplast genome study on Callistemon rigidus R.Br. is a few seen. This study demonstrates the data collected from the assembly and annotation of the chloroplast (cp) genome of Callistemon rigidus R.Br., followed by furthers comparative analysis with the cp genomes of closely related species. C. rigidus R.Br. showed a cp genome in the size of 158, 961 bp long with 36.78% GC content, among which a pair of inverted repeats (IRs) of 26, 671 bp separated a large single-copy (LSC) region of 87, 162 bp and a small single-copy (SSC) region of 18, 457 bp. Altogether 131 genes were hosted, including 37 transfer RNAs, 8 ribosomal RNAs, and 86 protein-coding genes. 284 simple sequence repeats (SSRs) were also marked out. A comparative analysis of the genome structure and the sequence data of closely related species unveiled the conserved gene order in the IR and LSC/SSC regions, a quite constructive finding for future phylogenetic research. Overall, this study providing C. rigidus R.Br. genomic resources could positively contribute to the evolutionary study and the phylogenetic reconstruction of Myrtaceae.

Discovery and Biomimetic Synthesis of a Phloroglucinol-Terpene Adduct Collection from Baeckea frutescens and Its Biogenetic Origin Insight.

A phloroglucinol-terpene adduct (PTA) collection consisting of twenty-four molecules featuring three skeletons was discovered from Baeckea frutescens. Inspired by its biosynthetic hypothesis, we synthesized this PTA collection by reductive activation of stable phloroglucinol precursors into highly reactive ortho-quinone methide (o-QM) intermediates and subsequently Diels-Alder cycloaddition. We also demonstrated, for the first time, the generation process of the active o-QM by performing dynamic NMR and HPLC-MS monitoring experiments. Moreover, the PTA collection showed significant antifeedant effect toward the Plutella xylostella larvae.

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.

Updated Genome Assembly and Annotation for Metrosideros polymorpha, an Emerging Model Tree Species of Ecological Divergence.

Accurate feature annotation as well as assembly contiguity are important requisites of a modern genome assembly. They allow large-scale comparison of genomes across and within species and identification of polymorphisms, leading evolutionary and functional studies. We report an updated genome resource for Metrosideros polymorpha, the most dominant tree species in the Hawaiian native forests and a unique example of rapid and remarkable ecological diversification of woody species. Ninety-one percent of the bases in the sequence assembly (304 Mb) were organized into 11 pseudo-molecules, which would represent the chromosome structure of the species assuming the synteny to a close relative Eucalyptus Our complementary approach using manual annotation and automated pipelines identified 11.30% of the assembly to be transposable elements, in contrast to 4.1% in previous automated annotation. By increasing transcript and protein sequence data, we predicted 27,620 gene models with high concordance from the supplied evidence. We believe that this assembly, improved for contiguity, and annotation will be valuable for future evolutionary studies of M. polymorpha and closely related species, facilitating the isolation of specific genes and the investigation of genome-wide polymorphisms associated with ecological divergence.