A chromosome-level reference genome of an aromatic medicinal plant Adenosma buchneroides.

Adenosma buchneroides Bonati, belonging to the genus Adenosum (Plantaginaceae), is an aromatic medicinal plant and utilized in traditional Chinese medicine. It has been widely used as plant-based repellents to prevent vector-borne diseases. However, the lack of a reference genome limits the study of conservation management and molecular biology of A. buchneroides. Here, we generated a chromosome-level de novo genome assembly of A. buchneroides which is a high-quality chromosome-scale assembly of aromatic medicinal plant in Plantaginaceae. The genome has a total length of 442.84 Mb with scaffold N50 of 27.98 Mb and 95.55% of the genome assigned to 14 chromosomes. BUSCO assessment yielded a completeness score of 97.2%. Furthermore, we predicted 24,367 protein-coding genes, and 95.79% of them was functionally annotated. The chromosome-scale genome of A. buchneroides will be a significant resource for understanding the genetic basis and evolution of active components biosynthesis, which will facilitate further study and exploit of A. buchneroides.

Genome assembly of the pioneer species Plantago major L. (Plantaginaceae) provides insight into its global distribution and adaptation to metal-contaminated soil.

Plantago is a major genus belonging to the Plantaginaceae family and is used in herbal medicine, functional food, and pastures. Several Plantago species are also characterized by their global distribution, but the mechanism underpinning this is not known. Here, we present a high-quality, chromosome-level genome assembly of Plantago major L., a species of Plantago, by incorporating Oxford Nanopore sequencing and Hi-C technologies. The genome assembly size was approximately 671.27 Mb with a contig N50 length of 31.30 Mb. 31,654 protein-coding genes were identified from the genome. Evolutionary analysis showed that P. major diverged from other Lamiales species at ~62.18 Mya and experienced two rounds of WGD events. Notably, many gene families related to plant acclimation and adaptation expanded. We also found that many polyphenol biosynthesis genes showed high expression patterns in roots. Some amino acid biosynthesis genes, such as those involved in histidine synthesis, were highly induced under metal (Ni) stress that led to the accumulation of corresponding metabolites. These results suggest persuasive arguments for the global distribution of P. major through multiscale analysis. Decoding the P. major genome provides a valuable genomic resource for research on dissecting biological function, molecular evolution, taxonomy, and breeding.

Comprehensive analysis of complete chloroplast genome sequence of Plantago asiatica L. (Plantaginaceae).

Plantago asiatica L. is a representative individual species of Plantaginaceae, whose high reputation is owed to its edible and medicinal values. However, the phylogeny and genes of the P. asiatica chloroplast have not yet been well described. Here we report the findings of a comprehensive analysis of the P. asiatica chloroplast genome. The P. asiatica chloroplast genome is 164,992 bp, circular, and has a GC content of 37.98%. The circular genome contains 141 genes, including 8 rRNAs, 38 tRNAs, and 95 protein-coding genes. Seventy-two simple sequence repeats are detected. Comparative chloroplast genome analysis of six related species suggests that a higher similarity exists in the coding region than the non-coding region, and differences in the degree of preservation is smaller between P. asiatica and Plantago depressa than among others. Our phylogenetic analysis illustrates P. asiatica has a relatively close relationship with P. depressa, which was also divided into different clades with Plantago ovata and Plantago lagopus in the genus Plantago. This analysis of the P. asiatica chloroplast genome contributes to an improved deeply understanding of the evolutionary relationships among Plantaginaceae.

Callitriche as a potential model system for evolutionary studies on the dorsiventral distribution of stomata.

Controlling the distribution of stomata is crucial for the adaptation of plants to new, or changing environments. While many plant species produce stomata predominantly on the abaxial leaf surface (hypostomy), some produce stomata on both surfaces (amphistomy), and the remaining few produce them only on the adaxial surface (hyperstomy). Various selective pressures have driven the evolution of these three modes of stomatal distribution. Despite recent advances in our understanding of stomatal development and dorsiventral leaf polarity, the genetic basis for the evolution of different stomatal distributions is still unclear. Here, we propose the genus Callitriche as a new model system to investigate patterns in the evolution of stomatal distribution. Callitriche comprises species with diverse lifestyles, including terrestrial, amphibious, and obligately aquatic plants. We found that species in this genus cover all three modes of dorsiventral stomatal distribution, making it a desirable model for comparative and evolutionary analyses on distribution modes. We further characterized the genetic basis of the different distribution modes, focusing on the stomatal key transcription factor SPEECHLESS. Future research using the promising model system Callitriche would open a new direction for evolutionary developmental biology studies on stomata.

Identification of the unique molecular framework of heterophylly in the amphibious plant Callitriche palustris L.

Heterophylly is the development of different leaf forms in a single plant depending on the environmental conditions. It is often observed in amphibious aquatic plants that can grow under both aerial and submerged conditions. Although heterophylly is well recognized in aquatic plants, the associated developmental mechanisms and the molecular basis remain unclear. To clarify these underlying developmental and molecular mechanisms, we analyzed heterophyllous leaf formation in an aquatic plant, Callitriche palustris. Morphological analyses revealed extensive cell elongation and the rearrangement of cortical microtubules in the elongated submerged leaves of C. palustris. Our observations also suggested that gibberellin, ethylene, and abscisic acid all regulate the formation of submerged leaves. However, the perturbation of one or more of the hormones was insufficient to induce the formation of submerged leaves under aerial conditions. Finally, we analyzed gene expression changes during aerial and submerged leaf development and narrowed down the candidate genes controlling heterophylly via transcriptomic comparisons, including a comparison with a closely related terrestrial species. We discovered that the molecular mechanism regulating heterophylly in C. palustris is associated with hormonal changes and diverse transcription factor gene expression profiles, suggesting differences from the corresponding mechanisms in previously investigated amphibious plants.

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.

New genomic resources and comparative analyses reveal differences in floral gene expression in selfing and outcrossing Collinsia sister species.

The evolutionary transition from outcross- to self-fertilization is one of the most common in angiosperms and is often associated with a parallel shift in floral morphological and developmental traits, such as reduced flower size and pollen to ovule ratios, known as the "selfing syndrome." How these convergent phenotypes arise, the extent to which they are shaped by selection, and the nature of their underlying genetic basis are unsettled questions in evolutionary biology. The genus Collinsia (Plantaginaceae) includes seven independent transitions from outcrossing or mixed mating to high selfing rates accompanied by selfing syndrome traits. Accordingly, Collinsia represents an ideal system for investigating this parallelism, but requires genomic resource development. We present a high quality de novo genome assembly for the highly selfing species Collinsia rattanii. To begin addressing the basis of selfing syndrome developmental shifts, we evaluate and contrast patterns of gene expression from floral transcriptomes across three stages of bud development for C. rattanii and its outcrossing sister species Collinsia linearis. Relative to C. linearis, total gene expression is less variable among individuals and bud stages in C. rattanii. In addition, there is a common pattern among differentially expressed genes: lower expression levels that are more constant across bud development in C. rattanii relative to C. linearis. Transcriptional regulation of enzymes involved in pollen formation specifically in early bud development may influence floral traits that distinguish selfing and outcrossing Collinsia species through pleiotropic functions. Future work will include additional Collinsia outcrossing-selfing species pairs to identify genomic signatures of parallel evolution.

The diversity of stomatal development regulation in Callitriche is related to the intrageneric diversity in lifestyles.

Stomata, the gas exchange structures of plants, are formed by the division and differentiation of stem cells, or meristemoids. Although diverse patterns of meristemoid behavior have been observed among different lineages of land plants, the ecological significance and diversification processes of these different patterns are not well understood. Here we describe an intrageneric diversity in the patterns of meristemoid division within the ecologically diverse genus Callitriche (Plantaginaceae). Meristemoids underwent a series of divisions before differentiating into stomata in the terrestrial species of Callitriche, but these divisions did not occur in amphibious species, which can grow in both air and water, in which meristemoids differentiated directly into stomata. These findings imply the adaptive significance of diversity in meristemoid division. Molecular genetic analyses showed that the different expression times of the stomatal key transcription factors SPEECHLESS and MUTE, which maintain and terminate the meristemoid division, respectively, underlie the different division patterns of meristemoids. Unlike terrestrial species, amphibious species prematurely expressed MUTE immediately after expressing SPEECHLESS, which corresponded to their early termination of stomatal division. By linking morphological, ecological, and genetic elements of stomatal development, this study provides significant insight that should aid ecological evolutionary developmental biology investigations of stomata.

Direct and indirect selection on mate choice during pollen competition: Effects of male and female sexual traits on offspring performance following two-donor crosses.

Mate choice in plants is poorly understood, in particular its indirect genetic benefits, but also the direct benefits of avoiding harmful matings. In the herb Collinsia heterophylla, delayed stigma receptivity has been suggested to enhance pollen competition, potentially functioning as a female mate choice trait. Previous studies show that this trait can mitigate the cost of early fertilization caused by pollen, thus providing a direct benefit. We performed two-donor pollinations during successive floral stages to assess how this stigma receptivity trait and two pollen traits known to affect siring success influence indirect benefits in terms of offspring performance. We also investigated differential resource allocation by studying the influence of sibling performance in the same capsule. Offspring performance in terms of flower number was mainly affected by parental identities and differential resource allocation. Offspring seed production showed some influence of resource allocation, but was also affected by pollen donor identity and varied positively with late stigma receptivity. However, the effect of late stigma receptivity on offspring seed production was weakened in matings with pollen that advanced stigma receptivity. In conclusion, delayed stigma receptivity may be selected through both direct and indirect fitness effects in C. heterophylla, where pollen-based delay on stigma receptivity might act as a cue for mate choice. However, selection may also be counteracted by antagonistic selection on pollen to advance stigma receptivity. Our results highlight the challenges of studying indirect genetic benefits and other factors that influence mate choice in plants.

Phylogeography of western Mediterranean Cymbalaria (Plantaginaceae) reveals two independent long-distance dispersals and entails new taxonomic circumscriptions.

The Balearic Islands, Corsica and Sardinia (BCS) constitute biodiversity hotspots in the western Mediterranean Basin. Oligocene connections and long distance dispersal events have been suggested to cause presence of BCS shared endemic species. One of them is Cymbalaria aequitriloba, which, together with three additional species, constitute a polyploid clade endemic to BCS. Combining amplified fragment length polymorphism (AFLP) fingerprinting, plastid DNA sequences and morphometrics, we inferred the phylogeography of the group and evaluated the species' current taxonomic circumscriptions. Based on morphometric and AFLP data we propose a new circumscription for C. fragilis to additionally comprise a group of populations with intermediate morphological characters previously included in C. aequitriloba. Consequently, we suggest to change the IUCN category of C. fragilis from critically endangered (CR) to near threatened (NT). Both morphology and AFLP data support the current taxonomy of the single island endemics C. hepaticifolia and C. muelleri. The four species had a common origin in Corsica-Sardinia, and two long-distance dispersal events to the Balearic Islands were inferred. Finally, plastid DNA data suggest that interspecific gene flow took place where two species co-occur.

OneTwoTree: An online tool for phylogeny reconstruction.

Phylogeny reconstruction is a key instrument in numerous biological analyses, ranging from evolutionary and ecology research, to conservation and systems biology. The increasing accumulation of genomic data makes it possible to reconstruct phylogenies with both high accuracy and at increasingly finer resolution. Yet, taking advantage of the enormous amount of sequence data available requires the use of computational tools for efficient data retrieval and processing, or else the process could quickly become an error-prone endeavour. Here, we present OneTwoTree (http://onetwotree.tau.ac.il/), a Web-based tool for tree reconstruction based on the supermatrix paradigm. Given a list of taxa names of interest as the sole input requirement, OneTwoTree retrieves all available sequence data from NCBI GenBank, clusters these into orthology groups, identifies the most informative set of markers, searches for an appropriate outgroup, and assembles a partitioned sequence matrix that is then used for the final phylogeny reconstruction step. OneTwoTree further allows users to control various steps of the process, such as the merging of sequences from similar clusters, or phylogeny reconstruction based on markers from a specific genome type. By comparing the performance of OneTwoTree to a manually reconstructed phylogeny of the Antirrhineae tribe, we show that the use of OneTwoTree resulted in substantially higher data coverage in terms of both taxon sampling and the number of informative markers assembled. OneTwoTree provides a flexible online tool for species-tree reconstruction, aimed to assist researchers ranging in their level of prior expertise in the task of phylogeny reconstruction.

Bayesian inference of phylogenetic networks from bi-allelic genetic markers.

Phylogenetic networks are rooted, directed, acyclic graphs that model reticulate evolutionary histories. Recently, statistical methods were devised for inferring such networks from either gene tree estimates or the sequence alignments of multiple unlinked loci. Bi-allelic markers, most notably single nucleotide polymorphisms (SNPs) and amplified fragment length polymorphisms (AFLPs), provide a powerful source of genome-wide data. In a recent paper, a method called SNAPP was introduced for statistical inference of species trees from unlinked bi-allelic markers. The generative process assumed by the method combined both a model of evolution for the bi-allelic markers, as well as the multispecies coalescent. A novel component of the method was a polynomial-time algorithm for exact computation of the likelihood of a fixed species tree via integration over all possible gene trees for a given marker. Here we report on a method for Bayesian inference of phylogenetic networks from bi-allelic markers. Our method significantly extends the algorithm for exact computation of phylogenetic network likelihood via integration over all possible gene trees. Unlike the case of species trees, the algorithm is no longer polynomial-time on all instances of phylogenetic networks. Furthermore, the method utilizes a reversible-jump MCMC technique to sample the posterior of phylogenetic networks given bi-allelic marker data. Our method has a very good performance in terms of accuracy and robustness as we demonstrate on simulated data, as well as a data set of multiple New Zealand species of the plant genus Ourisia (Plantaginaceae). We implemented the method in the publicly available, open-source PhyloNet software package.

Microsatellite primers in Plantago virginica (Plantaginaceae), an invasive species with both cleistogamous and chasmogamous flowers.

Polymorphic microsatellite markers were developed in Plantago virginica, an invasive species in China with both cleistogamous and chasmogamous flowers, to investigate its genetic structure and mating patterns. Fourteen novel microsatellite primer sets were designed, and the marker loci they amplified were characterized in 96 individuals from four populations. Eleven of these markers showed polymorphism and the number of alleles per locus ranged from two to six. AMOVA and STRUCTURE indicated that there were distinct patterns of genetic differentiation among the one invasive and three native US populations. These markers provide a useful tool for investigating genetic diversity in P. virginica and studying the mechanisms of invasion success.

Sexually antagonistic evolution caused by male-male competition in the pistil.

Although sexual selection and sexual conflict are important evolutionary forces in animals, their significance in plants is uncertain. In hermaphroditic organisms, such as many plants, sexual conflict may occur both between mating partners (interlocus conflict) and between male and female sex roles within an individual (intralocus conflict). We performed experimental evolution, involving lines that were crossed with either one or two pollen donors (monogamous or polyandrous lines), in the hermaphroditic plant (Collinsia heterophylla) where early fertilizations are associated with female fitness costs (reduced seed set). Artificial polyandry for four generations resulted in enhanced pollen performance and increased female fitness costs compared to the monogamous and source (starting material) lines. Female fitness was also reduced in the monogamous line, indicating a possible trade-off between sex roles, resulting from early pollination. We performed a second experiment to investigate a potential harming effect of pollen performance on seed set. We found that high siring success of early arriving pollen competing with later-arriving pollen was associated with high female fitness costs, consistent with an interlocus sexual conflict. Our study provides evidence for the importance of sexual selection in shaping evolution of plant reproductive strategies, but also pinpoints the complexity of sexual conflict in hermaphroditic species.

Travelling to the south: Phylogeographic spatial diffusion model in Monttea aphylla (Plantaginaceae), an endemic plant of the Monte Desert.

Effects of Pleistocene climatic oscillations on plant phylogeographic patterns are relatively well studied in forest, savanna and grassland biomes, but such impacts remain less explored on desert regions of the world, especially in South America. Here, we performed a phylogeographical study of Monttea aphylla, an endemic species of the Monte Desert, to understand the evolutionary history of vegetation communities inhabiting the South American Arid Diagonal. We obtained sequences of three chloroplast (trnS-trnfM, trnH-psbA and trnQ-rps16) and one nuclear (ITS) intergenic spacers from 272 individuals of 34 localities throughout the range of the species. Population genetic and Bayesian coalescent analyses were performed to infer genealogical relationships among haplotypes, population genetic structure, and demographic history of the study species. Timing of demographic events was inferred using Bayesian Skyline Plot and the spatio-temporal patterns of lineage diversification was reconstructed using Bayesian relaxed diffusion models. Palaeo-distribution models (PDM) were performed through three different timescales to validate phylogeographical patterns. Twenty-five and 22 haplotypes were identified in the cpDNA and nDNA data, respectively. that clustered into two main genealogical lineages following a latitudinal pattern, the northern and the southern Monte (south of 35° S). The northern Monte showed two lineages of high genetic structure, and more relative stable demography than the southern Monte that retrieved three groups with little phylogenetic structure and a strong signal of demographic expansion that would have started during the Last Interglacial period (ca. 120 Ka). The PDM and diffusion models analyses agreed in the southeast direction of the range expansion. Differential effect of climatic oscillations across the Monte phytogeographic province was observed in Monttea aphylla lineages. In northern Monte, greater genetic structure and more relative stable demography resulted from a more stable climate than in the southern Monte. Pleistocene glaciations drastically decreased the species area in the southern Monte, which expanded in a southeastern direction to the new available areas during the interglacial periods.

Genetic and geographical structure of boreal plants in their southern range: phylogeography of Hippuris vulgaris in China.

BACKGROUND: Our current understanding of the evolutionary history of boreal and arctic-alpine plants in their southern range in Asia remains relatively poor. Using three cpDNA non-coding regions and nine nuclear microsatellite (nSSR) loci, we examine the phylogeographic pattern in a broad geographic sampling of the boreal plant Hippuris vulgaris to infer its dispersal and diversification in China. In addition, the species distributions at the Last Glacial Maximum (LGM) and at present were predicted using ecological niche modeling (ENM). RESULTS: The cpDNA results revealed two distinct lineages, A and B. A is restricted to Northeast China; B is distributed in Northwest China, the Qinghai-Tibet Plateau (QTP) and North and Northeast (NNE) China; and A and B diverged ca. 1.36 Ma. The nSSR data revealed two genetic clusters corresponding to the two cpDNA lineages and nonreciprocal hybridization with lineage A as the maternal lineage in Northeast China. Cluster B further divided into three subclusters: I, mainly in NNE China and the northeastern border of the QTP; II, in Northwest China and the QTP; and III, on the QTP. ENM predicted a marked range shift on the QTP at the LGM, retreating from the platform to the northeast and southeast edges. CONCLUSIONS: Hippuris vulgaris probably diverged into lineages A and B in high latitudes and then immigrated into Northeast China and Northwest China, respectively. Lineage A persisted and diversified in Northeast China. Lineage B reached the QTP during the mid-Pleistocene, diversified in that region due to the influence of climatic oscillations, migrated into Northeast China and subsequently hybridized with lineage A. Our findings give empirical evidence that boreal plants display complex evolutionary history in their southern range in Asia and provide new insights into the evolution of boreal and arctic-alpine plants.

No evidence for Fabaceae Gametophytic self-incompatibility being determined by Rosaceae, Solanaceae, and Plantaginaceae S-RNase lineage genes.

BACKGROUND: Fabaceae species are important in agronomy and livestock nourishment. They have a long breeding history, and most cultivars have lost self-incompatibility (SI), a genetic barrier to self-fertilization. Nevertheless, to improve legume crop breeding, crosses with wild SI relatives of the cultivated varieties are often performed. Therefore, it is fundamental to characterize Fabaceae SI system(s). We address the hypothesis of Fabaceae gametophytic (G)SI being RNase based, by recruiting the same S-RNase lineage gene of Rosaceae, Solanaceae or Plantaginaceae SI species. RESULTS: We first identify SSK1 like genes (described only in species having RNase based GSI), in the Trifolium pratense, Medicago truncatula, Cicer arietinum, Glycine max, and Lupinus angustifolius genomes. Then, we characterize the S-lineage T2-RNase genes in these genomes. In T. pratense, M. truncatula, and C. arietinum we identify S-RNase lineage genes that in phylogenetic analyses cluster with Pyrinae S-RNases. In M. truncatula and C. arietinum genomes, where large scaffolds are available, these sequences are surrounded by F-box genes that in phylogenetic analyses also cluster with S-pollen genes. In T. pratense the S-RNase lineage genes show, however, expression in tissues not involved in GSI. Moreover, levels of diversity are lower than those observed for other S-RNase genes. The M. truncatula and C. arietinum S-RNase and S-pollen like genes phylogenetically related to Pyrinae S-genes, are also expressed in tissues other than those involved in GSI. To address if other T2-RNases could be determining Fabaceae GSI, here we obtained a style with stigma transcriptome of Cytisus striatus, a species that shows significant difference on the percentage of pollen growth in self and cross-pollinations. Expression and polymorphism analyses of the C. striatus S-RNase like genes revealed that none of these genes, is the S-pistil gene. CONCLUSION: We find no evidence for Fabaceae GSI being determined by Rosaceae, Solanaceae, and Plantaginaceae S-RNase lineage genes. There is no evidence that T2-RNase lineage genes could be determining GSI in C. striatus. Therefore, to characterize the Fabaceae S-pistil gene(s), expression analyses, levels of diversity, and segregation analyses in controlled crosses are needed for those genes showing high expression levels in the tissues where GSI occurs.

Sexual antagonism in the pistil varies among populations of a hermaphroditic mixed-mating plant.

Sexual conflicts and their evolutionary outcomes may be influenced by population-specific features such as mating system and ecological context; however, very few studies have investigated the link between sexual conflict and mating system. The self-compatible, mixed-mating hermaphrodite Collinsia heterophylla (Plantaginaceae) is thought to exhibit a sexual conflict over timing of stigma receptivity. This conflict involves (i) delayed stigma receptivity, which intensifies pollen competition, and (ii) early fertilization forced by pollen, which reduces seed set. We investigated the potential for the conflict to occur under field conditions and performed glasshouse crosses within eight populations to assess its consistency across populations. Flowers were visited, and produced seeds after pollination, at all developmental stages, suggesting that the conflict can be of significance under natural conditions. In the glasshouse, early pollination imposed costs in all populations. Overall, the timing of first seed set was most strongly affected by the maternal parent, denoting stronger female than male ability to influence the onset of stigma receptivity. Crosses also revealed a negative relationship between donor- and recipient-related onset of receptivity within individuals, a novel result hinting at trade-offs in sex allocation or a history of antagonistic selection. Neither timing of stigma receptivity, timing of first seed set, nor pollen competitive ability covaried with population outcrossing rate. In conclusion, these results indicate that sexually antagonistic selection may be present in varying degrees in different populations of C. heterophylla, but this variation does not appear to be directly related to mating system variation.