Ancestral range reconstruction of remote oceanic island species of Plantago (Plantaginaceae) reveals differing scales and modes of dispersal

Macroevolutionary decline in mycorrhizal colonization and chemical defense responsiveness to mycorrhization

Arbuscular mycorrhizal fungi (AMF) have evolved associations with roots of 60% plant species, but the net benefit for plants vary broadly from mutualism to parasitism. Yet, we lack a general understanding of the evolutionary and ecological forces driving such variation. To this end, we conducted a comparative phylogenetic experiment with 24 species of Plantago, encompassing worldwide distribution, to address the effect of evolutionary history and environment on plant growth and chemical defenses in response to AMF colonization. We demonstrate that different species within one plant genus vary greatly in their ability to associate with AMF, and that AMF arbuscule colonization intensity decreases monotonically with increasing phylogenetic branch length, but not with concomitant changes in pedological and climatic conditions across species. Moreover, we demonstrate that species with the highest colonization levels are also those that change their defensive chemistry the least. We propose that the costs imposed by high AMF colonization in terms of reduced changes in secondary chemistry might drive the observed macroevolutionary decline in mycorrhization.

A comparison of freezer-stored DNA and herbarium tissue samples for chloroplast assembly and genome skimming

Premise

The use of DNA from herbarium specimens is an increasingly important source for evolutionary studies in plant biology, particularly in cases where species are rare or difficult to obtain. Here we compare the utility of DNA from herbarium tissues to their freezer-stored DNA counterparts via the Hawaiian Plant DNA Library.

Methods

Plants collected for the Hawaiian Plant DNA Library were simultaneously accessioned as herbarium specimens at the time of collection, from 1994-2019. Paired samples were sequenced using short-read sequencing and assessed for chloroplast assembly and nuclear gene recovery.

Results

Herbarium specimen-derived DNA was statistically more fragmented than freezer-stored DNA derived from fresh tissue, leading to poorer chloroplast assembly and overall lower coverage. The number of nuclear targets recovered varied mostly by total sequencing reads per library and age of specimen, but not by storage method (herbarium or long-term freezer). Although there was evidence of DNA damage in the samples, there was no evidence that it was related to the length of time in storage, whether frozen or as herbarium specimens.

Discussion

DNA extracted from herbarium tissues will continue to be invaluable, despite being highly fragmented and degraded. Rare floras would benefit from both traditional herbarium storage methods and extracted DNA freezer banks.

Living on the edge: morphological, karyological and genetic diversity studies of the Hungarian Plantago maxima populations and established ex situ collection

BACKGROUND

The analysis of genetic diversity of protected plant species can greatly support conservation efforts. Plantago maxima Juss. ex Jacq. is a perennial species distributed along the Eurasian steppe. The westernmost range edge of the species' distribution is located in the Pannonian basin, in Hungary where it is represented by a few, fragmented and highly endangered populations. We studied population diversity of all Hungarian range edge, natural populations, and one established ex situ population. One population from the centre of distribution (Kazakhstan) was implemented in the cpDNA haplotype study to compare the peripheral vs. central populations. We performed morphometric trait-based analysis, chromosome studies (morphometric analyses and FISH) and genetic diversity evaluations using inter simple sequence repeats (ISSR) and cpDNA trnL-trnF to evaluate differences between the in situ and ex situ populations as well as central vs. peripheral populations.

RESULTS

Our results showed no obvious morphological differences among the in situ and ex situ populations in the period between 2018 and 2020. One ex situ subpopulation develops flowers three years in a row from 2019, which is a favourable indicator of the introduction success. Hungarian populations are exclusively diploids (2n = 2x = 12). The karyogram consists of 5 metacentric and 1 acrocentric chromosome pair. Plantago maxima has one 35S and two 5S rDNA loci, located on the acrocentric chromosome pair. Eight variable ISSR primers yielded 100 fragments, of which 74.6% were polymorphic (mean H_e = 0.220). A high level of genetic variation within population was observed (92%) while the genetic differentiation among the populations was only 8%. STRUCTURE analysis revealed that the largest Kunpeszér population separated from the rest of the Hungarian populations, indicating a high rate of admixture among the other ones. Based on the trnL-trnF sequence analysis the Hungarian populations represent a single haplotype, which can indicate a reduced diversity due to isolation and recent population decline. By contrast, Kazakh population represents a distinct haplotype compared to the Hungarian samples.

CONCLUSIONS

The present study draws the attention to the high conservation value of the Plantago maxima populations from the westernmost range edge of the species' distribution.

Travel Tales of a Worldwide Weed: Genomic Signatures of Plantago major L. Reveal Distinct Genotypic Groups With Links to Colonial Trade Routes

Retracing pathways of historical species introductions is fundamental to understanding the factors involved in the successful colonization and spread, centuries after a species' establishment in an introduced range. Numerous plants have been introduced to regions outside their native ranges both intentionally and accidentally by European voyagers and early colonists making transoceanic journeys; however, records are scarce to document this. We use genotyping-by-sequencing and genotype-likelihood methods on the selfing, global weed, Plantago major, collected from 50 populations worldwide to investigate how patterns of genomic diversity are distributed among populations of this global weed. Although genomic differentiation among populations is found to be low, we identify six unique genotype groups showing very little sign of admixture and low degree of outcrossing among them. We show that genotype groups are latitudinally restricted, and that more than one successful genotype colonized and spread into the introduced ranges. With the exception of New Zealand, only one genotype group is present in the Southern Hemisphere. Three of the most prevalent genotypes present in the native Eurasian range gave rise to introduced populations in the Americas, Africa, Australia, and New Zealand, which could lend support to the hypothesis that P. major was unknowlingly dispersed by early European colonists. Dispersal of multiple successful genotypes is a likely reason for success. Genomic signatures and phylogeographic methods can provide new perspectives on the drivers behind the historic introductions and the successful colonization of introduced species, contributing to our understanding of the role of genomic variation for successful establishment of introduced taxa.

Predominantly Eastward Long-Distance Dispersal in Pantropical Ochnaceae Inferred From Ancestral Range Estimation and Phylogenomics

Ochnaceae is a pantropical family with multiple transoceanic disjunctions at deep and shallow levels. Earlier attempts to unravel the processes that led to such biogeographic patterns suffered from insufficient phylogenetic resolution and unclear delimitation of some of the genera. In the present study, we estimated divergence time and ancestral ranges based on a phylogenomic framework with a well-resolved phylogenetic backbone to tackle issues of the timing and direction of dispersal that may explain the modern global distribution of Ochnaceae. The nuclear data provided the more robust framework for divergence time estimation compared to the plastome-scale data, although differences in the inferred clade ages were mostly small. While Ochnaceae most likely originated in West Gondwana during the Late Cretaceous, all crown-group disjunctions are inferred as dispersal-based, most of them as transoceanic long-distance dispersal (LDD) during the Cenozoic. All LDDs occurred in an eastward direction except for the SE Asian clade of Sauvagesieae, which was founded by trans-Pacific dispersal from South America. The most species-rich clade by far, Ochninae, originated from either a widespread neotropical-African ancestor or a solely neotropical ancestor which then dispersed to Africa. The ancestors of this clade then diversified in Africa, followed by subsequent dispersal to the Malagasy region and tropical Asia on multiple instances in three genera during the Miocene-Pliocene. In particular, Ochna might have used the South Arabian land corridor to reach South Asia. Thus, the pantropical distribution of Ochnaceae is the result of LDD either transoceanic or via land bridges/corridors, whereas vicariance might have played a role only along the stem of the family.

Molecular and Morphological Data Improve the Classification of Plantagineae (Lamiales)

The tribe Plantagineae (Lamiales) is a group of plants with worldwide distribution, notorious for its complicated taxonomy and still unresolved natural history. We describe the result of a broadly sampled phylogenetic study of tribe. The expanded sampling dataset is based on the trnL-F spacer, rbcL, and ITS2 markers across all three included genera (Aragoa, Littorella and Plantago) and makes this the most comprehensive study to date. The other dataset uses five markers and provides remarkably good resolution throughout the tree, including support for all of the major clades. In addition to the molecular phylogeny, a morphology database of 114 binary characters was assembled to provide comparison with the molecular phylogeny and to develop a means to assign species not sampled in the molecular analysis to their most closely related species that were sampled. Based on the molecular phylogeny and the assignment algorithm to place unsampled species, a key to sections is presented, and a revised classification of the tribe is provided. We also include the description of new species from North America.

Plantago campestris (Plantaginaceae), a rare new species from southern Brazil, supported by phylogenomic and morphological evidence

High-throughput sequencing, when combined with taxonomic expertise, is a powerful tool to refine and advance taxonomic classification, including at the species level. In the present work, a new species, Plantago campestris, is described out of the P. commersoniana species complex, based on phylogenomic and morphological evidence. The main morphological characters that distinguish the new species from P. commersoniana are the glabrous posterior sepals and the slightly broader leaves. The new species is known from only three localities, all in natural high-elevation grasslands in Paraná and Santa Catarina states, southern Brazil. According to the IUCN criteria new species should be assessed as Endangered (EN). We present field photographs of P. campestris and related species, and we provide an identification key to the species previously included within the circumscription of P. commersoniana.

Plastomes from tribe Plantagineae (Plantaginaceae) reveal infrageneric structural synapormorphies and localized hypermutation for Plantago and functional loss of ndh genes from Littorella

Tribe Plantagineae (Plantaginaceae) comprises ~ 270 species in three currently recognized genera (Aragoa, Littorella, Plantago), of which Plantago is most speciose. Plantago plastomes exhibit several atypical features including large inversions, expansions of the inverted repeat, increased repetitiveness, intron losses, and gene-specific increases in substitution rate, but the prevalence of these plastid features among species and subgenera is unknown. To assess phylogenetic relationships and plastomic evolutionary dynamics among Plantagineae genera and Plantago subgenera, we generated 25 complete plastome sequences and compared them with existing plastome sequences from Plantaginaceae. Using whole plastome and partitioned alignments, our phylogenomic analyses provided strong support for relationships among major Plantagineae lineages. General plastid features-including size, GC content, intron content, and indels-provided additional support that reinforced major Plantagineae subdivisions. Plastomes from Plantago subgenera Plantago and Coronopus have synapomorphic expansions and inversions affecting the size and gene order of the inverted repeats, and particular genes near the inversion breakpoints exhibit accelerated nucleotide substitution rates, suggesting localized hypermutation associated with rearrangements. The Littorella plastome lacks functional copies of ndh genes, which may be related to an amphibious lifestyle and partial reliance on CAM photosynthesis.

Plantago spp. as Models for Studying the Ecology and Evolution of Species Interactions across Environmental Gradients

AbstractA central challenge in ecology and evolutionary biology is to understand how variation in abiotic and biotic factors combine to shape the distribution, abundance, and diversity of focal species. Environmental gradients, whether natural (e.g., latitude, elevation, ocean proximity) or anthropogenic (e.g., land-use intensity, urbanization), provide compelling settings for addressing this challenge. However, not all organisms are amenable to the observational and experimental approaches required for untangling the factors that structure species along gradients. Here we highlight herbaceous plants in the genus Plantago as models for studying the ecology and evolution of species interactions along abiotic gradients. Plantago lanceolata and P. major are native to Europe and Asia but distributed globally, and they are established models for studying population ecology and interactions with herbivores, pathogens, and soil microbes. Studying restricted range congeners in comparison with those cosmopolitan species can provide insight into abiotic and biotic determinants of range size and population structure. We highlight one such species, P. rugelii, which is endemic to eastern North America. We give an overview of the literature on these focal Plantago species and explain why they are logical candidates for studies of species interactions across environmental gradients. Finally, we emphasize collaborative and community science approaches that can facilitate such research and note the amenability of Plantago for authentic research projects in science education.

Expanded inverted repeat region with large scale inversion in the first complete plastid genome sequence of Plantago ovata

Plantago ovata (Plantaginaceae) is an economically and medicinally important species, however, least is known about its genomics and evolution. Here, we report the first complete plastome genome of P. ovata and comparison with previously published genomes of related species from Plantaginaceae. The results revealed that P. ovata plastome size was 162,116 bp and that it had typical quadripartite structure containing a large single copy region of 82,084 bp and small single copy region of 5,272 bp. The genome has a markedly higher inverted repeat (IR) size of 37.4 kb, suggesting large-scale inversion of 13.8 kb within the expanded IR regions. In addition, the P. ovata plastome contains 149 different genes, including 43 tRNA, 8 rRNA, and 98 protein-coding genes. The analysis revealed 139 microsatellites, of which 71 were in the non-coding regions. Approximately 32 forward, 34 tandem, and 17 palindromic repeats were detected. The complete genome sequences, 72 shared genes, matK gene, and rbcL gene from related species generated the same phylogenetic signals, and phylogenetic analysis revealed that P. ovata formed a single clade with P. maritima and P. media. The divergence time estimation as employed in BEAST revealed that P. ovata diverged from P. maritima and P. media about 11.0 million years ago (Mya; 95% highest posterior density, 10.06-12.25 Mya). In conclusion, P. ovata had significant variation in the IR region, suggesting a more stable P. ovata plastome genome than that of other Plantaginaceae species.