From gene trees to a dated allopolyploid network: insights from the angiosperm genus Viola (Violaceae)

Cyclotides derived from Viola dalatensis Gagnep: A novel approach for enrichment and evaluation of antimicrobial activity

Cyclotides, plant-derived cysteine-rich peptides, exhibit a wide range of beneficial biological activities and possess exceptional structural stability. Cyclotides are commonly distributed throughout the Violaceae family. Viola dalatensis Gagnep, a Vietnamese species, has not been well studied, especially for cyclotides. This pioneering research explores cyclotides from V. dalatensis as antimicrobials. This study used a novel approach to enhance cyclotides after extraction. The approach combined 30% ammonium sulfate salt precipitation and RP-HPLC. A comprehensive analysis was performed to ascertain the overall protein content, flavonoids content, polyphenol content, and free radical scavenging capacity of compounds derived from V. dalatensis. Six known cyclotides were sequenced utilizing MS tandem. Semi-purified cyclotide mixtures (M1, M2, and M3) exhibited antibacterial efficacy against Bacillus subtilis (inhibitory diameters: 19.67-23.50 mm), Pseudomonas aeruginosa (22.17-23.50 mm), and Aspergillus flavus (14.67-21.33 mm). The enriched cyclotide precipitate from the stem extract demonstrated a minimum inhibitory concentration (MIC) of 0.08 mg/mL against P. aeruginosa, showcasing significant antibacterial effectiveness compared to the stem extract (MIC: 12.50 mg/mL). Considerable advancements have been achieved in the realm of cyclotides, specifically in their application as antimicrobial agents.

Phylogenomics reveals widespread hybridization and polyploidization in Henckelia (Gesneriaceae)

BACKGROUND AND AIMS

Hybridization has long been recognized as an important process for plant evolution and is often accompanied by polyploidization, another prominent force in generating biodiversity. Despite its pivotal importance in evolution, the actual prevalence and distribution of hybridization across the tree of life remain unclear.

METHODS

We used whole-genome shotgun (WGS) sequencing and cytological data to investigate the evolutionary history of Henckelia, a large genus in the family Gesneriaceae with a high frequency of suspected hybridization and polyploidization events. We generated WGS sequencing data at about 10× coverage for 26 Chinese Henckelia species plus one Sri Lankan species. To untangle the hybridization history, we separately extracted whole plastomes and thousands of single-copy nuclear genes from the sequencing data, and reconstructed phylogenies based on both nuclear and plastid data. We also explored sources of both genealogical and cytonuclear conflicts and identified signals of hybridization and introgression within our phylogenomic dataset using several statistical methods. Additionally, to test the polyploidization history, we evaluated chromosome counts for 45 populations of the 27 Henckelia species studied.

KEY RESULTS

We obtained well-supported phylogenetic relationships using both concatenation- and coalescent-based methods. However, the nuclear phylogenies were highly inconsistent with the plastid phylogeny, and we observed intensive discordance among nuclear gene trees. Further analyses suggested that both incomplete lineage sorting and gene flow contributed to the observed cytonuclear and genealogical discordance. Our analyses of introgression and phylogenetic networks revealed a complex history of hybridization within the genus Henckelia. In addition, based on chromosome counts for 27 Henckelia species, we found independent polyploidization events occurred within Henckelia after different hybridization events.

CONCLUSIONS

Our findings demonstrated that hybridization and polyploidization are common in Henckelia. Furthermore, our results revealed that H. oblongifolia is not a member of the redefined Henckelia and they suggested several other taxonomic treatments in this genus.

Genetic Variation versus Morphological Variability in European Peatland Violets (Viola epipsila—V. palustris Group)

In Europe, the V. epipsila—V. palustris group comprises V. epipsila Ledeb., V. palustris L., V. pubifolia (Kuta) G. H. Loos (=V. palustris subsp. pubifolia Kuta), interspecific hybrids, and putative introgressants. The genetic affinity of V. pubifolia to V. palustris, and their shared origin via hybridization followed by polyploidization, were confirmed using inter simple sequence repeat (ISSR) markers, restriction site-associated DNA sequencing (RAD-Seq), and a low-copy nuclear gene, GPI, which encodes glucose-6-phosphate isomerase. The other taxa of subsect. Stolonosae were not identified as putative parents of V. pubifolia by GPI. Our analyses indicated that V. pubifolia can be included in the morphological and genetic variation of V. palustris. The ISSR, RAD-Seq, and genome size value separated well V. palustris from V. epipsila and hybrids. The results also reopen the discussion on intraspecific variation in the context of taxa ranks and species concepts. The reduced tolerance of V. epipsila in Europe to changing environmental conditions might result from low genetic differentiation and heterozygosity, as well as the increased number of interspecific hybrids (V. epipsila × V. palustris), and eventually can possibly lead to its extinction. The disappearance of populations/individuals of this species may indicate anthropogenic changes occurring in peatlands.

Picks in the Fabric of a Polyploidy Complex: Integrative Species Delimitation in the Tetraploid Leucanthemum Mill. (Compositae, Anthemideae) Representatives

Based on the results of a preceding species-delimitation analysis for the diploid representatives of the genus Leucanthemum (Compositae, Anthemideae), the present study aims at the elaboration of a specific and subspecific taxonomic treatment of the tetraploid members of the genus. Following an integrative taxonomic approach, species-level decisions on eight predefined morphotaxon hypotheses were based on genetic/genealogical, morphological, ecological, and geographical differentiation patterns. ddRADseq fingerprinting and SNP-based clustering revealed genetic integrity for six of the eight morphotaxa, with no clear differentiation patterns observed between the widespread L. ircutianum subsp. ircutianum and the N Spanish (Cordillera Cantábrica) L. cantabricum and the S French L. delarbrei subsp. delabrei (northern Massif Central) and L. meridionale (western Massif Central). The inclusion of differentiation patterns in morphological (leaf dissection and shape), ecological (climatological and edaphic niches), and geographical respects (pair-wise tests of sympatry vs. allopatry) together with the application of a procedural protocol for species-rank decisions (the 'Wettstein tesseract') led to the proposal of an acknowledgement of the eight predefined morphotaxon hypotheses as six species (two of them with two subspecies). Nomenclatural consequences following from these results are drawn and lead to the following new combinations: Leucanthemum delarbrei subsp. meridionale (Legrand) Oberpr., T.Ott & Vogt, comb. nov. and Leucanthemum ruscinonense (Jeanb. & Timb.-Lagr.) Oberpr., T.Ott & Vogt, comb. et stat. nov.

Distinct-Cluster Tree-Child Phylogenetic Networks and Possible Uses to Study Polyploidy

As phylogenetic networks become more widely studied and the networks grow larger, it may be useful to "simplify" such networks into especially tractable networks. Recent results have found methods to simplify networks into normal networks. By definition, normal networks contain no redundant arcs. Nevertheless, there may be redundant arcs in networks where speciation events involving allopolyploidy occur. It is therefore desirable to find a different tractable class of networks that may contain redundant arcs. This paper proposes distinct-cluster tree-child networks as such a class, here abbreviated as DCTC networks. They are shown to have a number of useful properties, such as quadratic growth of the number of vertices with the number of leaves. A DCTC network is shown to be essentially a normal network to which some redundant arcs may have been added without losing the tree-child property. Every phylogenetic network can be simplified into a DCTC network depending only on the structure of the original network. There is always a CSD map from the original network to the resulting DCTC network. As a result, the simplified network can readily be interpreted via a "wired lift" in which the original network is redrawn with each arc represented in one of two ways.

The hybrid number of a ploidy profile

Polyploidization, whereby an organism inherits multiple copies of the genome of their parents, is an important evolutionary event that has been observed in plants and animals. One way to study such events is in terms of the ploidy number of the species that make up a dataset of interest. It is therefore natural to ask: How much information about the evolutionary past of the set of species that form a dataset can be gleaned from the ploidy numbers of the species? To help answer this question, we introduce and study the novel concept of a ploidy profile which allows us to formalize it in terms of a multiplicity vector indexed by the species the dataset is comprised of. Using the framework of a phylogenetic network, we present a closed formula for computing the hybrid number (i.e. the minimal number of polyploidization events required to explain a ploidy profile) of a large class of ploidy profiles. This formula relies on the construction of a certain phylogenetic network from the simplification sequence of a ploidy profile and the hybrid number of the ploidy profile with which this construction is initialized. Both of them can be computed easily in case the ploidy numbers that make up the ploidy profile are not too large. To help illustrate the applicability of our approach, we apply it to a simplified version of a publicly available Viola dataset.

Phylogenomics and Systematics of Overlooked Mesoamerican and South American Polyploid Broad-Leaved Festuca Grasses Differentiate F. sects. Glabricarpae and Ruprechtia and F. subgen. Asperifolia, Erosiflorae, Mallopetalon and Coironhuecu (subgen. nov.)

Allopolyploidy is considered a driver of diversity in subtribe Loliinae. We investigate the evolution and systematics of the poorly studied Mesoamerican and South American polyploid broad-leaved Festuca L. species of uncertain origin and unclear taxonomy. A taxonomic study of seven diagnostic morphological traits was conducted on a representation of 22 species. Phylogenomic analyses were performed on a representation of these supraspecific taxa and all other Loliinae lineages using separate data from the entire plastome, nuclear rDNA 45S and 5S genes, and repetitive DNA elements. F. subgen. Mallopetalon falls within the fine-leaved (FL) Loliinae clade, whereas the remaining taxa are nested within the broad-leaved (BL) Loliinae clade forming two separate Mexico–Central–South American (MCSAI, MCSAII) lineages. MCSAI includes representatives of F. sect. Glabricarpae and F. subgen. Asperifolia plus F. superba, and MCSAII of F. subgen. Erosiflorae and F. sect. Ruprechtia plus F. argentina. MCSAII likely had a BL Leucopoa paternal ancestor, MCSAI and MCSAII a BL Meso-South American maternal ancestor, and Mallopetalon FL, American I–II ancestors. Plastome vs. nuclear topological discordances corroborated the hybrid allopolyploid origins of these taxa, some of which probably originated from Northern Hemisphere ancestors. The observed data indicate rapid reticulate radiations in the Central–South American subcontinent. Our systematic study supports the reclassification of some studied taxa in different supraspecific Festuca ranks.

A Revised Phylogenetic Classification for Viola (Violaceae)

The genus Viola (Violaceae) is among the 40–50 largest genera among angiosperms, yet its taxonomy has not been revised for nearly a century. In the most recent revision, by Wilhelm Becker in 1925, the then-known 400 species were distributed among 14 sections and numerous unranked groups. Here, we provide an updated, comprehensive classification of the genus, based on data from phylogeny, morphology, chromosome counts, and ploidy, and based on modern principles of monophyly. The revision is presented as an annotated global checklist of accepted species of Viola, an updated multigene phylogenetic network and an ITS phylogeny with denser taxon sampling, a brief summary of the taxonomic changes from Becker’s classification and their justification, a morphological binary key to the accepted subgenera, sections and subsections, and an account of each infrageneric subdivision with justifications for delimitation and rank including a description, a list of apomorphies, molecular phylogenies where possible or relevant, a distribution map, and a list of included species. We distribute the 664 species accepted by us into 2 subgenera, 31 sections, and 20 subsections. We erect one new subgenus of Viola (subg. Neoandinium, a replacement name for the illegitimate subg. Andinium), six new sections (sect. Abyssinium, sect. Himalayum, sect. Melvio, sect. Nematocaulon, sect. Spathulidium, sect. Xanthidium), and seven new subsections (subsect. Australasiaticae, subsect. Bulbosae, subsect. Clausenianae, subsect. Cleistogamae, subsect. Dispares, subsect. Formosanae, subsect. Pseudorupestres). Evolution within the genus is discussed in light of biogeography, the fossil record, morphology, and particular traits. Viola is among very few temperate and widespread genera that originated in South America. The biggest identified knowledge gaps for Viola concern the South American taxa, for which basic knowledge from phylogeny, chromosome counts, and fossil data is virtually absent. Viola has also never been subject to comprehensive anatomical study. Studies into seed anatomy and morphology are required to understand the fossil record of the genus.

Application of chloroplast genome in the identification of Traditional Chinese Medicine Viola philippica

BACKGROUND

Viola philippica Cav. is the only source plant of "Zi Hua Di Ding", which is a Traditional Chinese Medicine (TCM) that is utilized as an antifebrile and detoxicant agent for the treatment of acute pyogenic infections. Historically, many Viola species with violet flowers have been misused in "Zi Hua Di Ding". Viola have been recognized as a taxonomically difficult genera due to their highly similar morphological characteristics. Here, all common V. philippica adulterants were sampled. A total of 24 complete chloroplast (cp) genomes were analyzed, among these 5 cp genome sequences were downloaded from GenBank and 19 cp genomes, including 2 "Zi Hua Di Ding" purchased from a local TCM pharmacy, were newly sequenced.

RESULTS

The Viola cp genomes ranged from 156,483 bp to 158,940 bp in length. A total of 110 unique genes were annotated, including 76 protein-coding genes, 30 tRNAs, and four rRNAs. Sequence divergence analysis screening identified 16 highly diverged sequences; these could be used as markers for the identification of Viola species. The morphological, maximum likelihood and Bayesian inference trees of whole cp genome sequences and highly diverged sequences were divided into five monophyletic clades. The species in each of the five clades were identical in their positions within the morphological and cp genome tree. The shared morphological characters belonging to each clade was summarized. Interestingly, unique variable sites were found in ndhF, rpl22, and ycf1 of V. philippica, and these sites can be selected to distinguish V. philippica from samples all other Viola species, including its most closely related species. In addition, important morphological characteristics were proposed to assist the identification of V. philippica. We applied these methods to examine 2 "Zi Hua Di Ding" randomly purchased from the local TCM pharmacy, and this analysis revealed that the morphological and molecular characteristics were valid for the identification of V. philippica.

CONCLUSIONS

This study provides invaluable data for the improvement of species identification and germplasm of V. philippica that may facilitate the application of a super-barcode in TCM identification and enable future studies on phylogenetic evolution and safe medical applications.

Violashiweii, a new species of Viola (Violaceae) from karst forest in Guizhou, China

Violashiweii Xiao C. Li & Z. W. Wang (Violaceae), a new species from Guizhou, China, is described, based on morphological and molecular evidence. The new species is morphologically most similar to V.kwangtungensis Melchior in its glabrous lateral petals and stoloniferous habit, but can be distinguished by its obtuse teeth along the leaf margin, its abaxially greyish-green leaf blade and its broader and entire sepals with a distinct basal appendage.

Complex Scenarios of Reticulation, Polyploidization, and Species Diversity within Annual Pansies of Subsect. Bracteolatae (Viola Sect. Melanium, Violaceae) in Italy: Insights from 5S-IGS High-Throughput Sequencing and Plastid DNA Variation

Viola sect. Melanium, the so-called pansy, is an allopolyploid morphologically well-defined lineage of ca. 110 perennial and annual species in the northern hemisphere, characterized by markedly complex genomic configurations. Five annual pansies occur in Italy, four of which are morphologically very similar and belong to the informal ‘V. tricolor species complex’: V. arvensis (2n = 34), V. hymettia (2n = 16), V. kitaibeliana (2n = 16), and V. tricolor (2n = 26). Their field recognition is difficult and reflects a long-debated taxonomy often resulting in doubtful records in field inventories and across European herbaria. The current lack of comprehensive intra- and interspecific comparative studies and a relative scarcity of appropriate genetic markers coupled with unambiguous cytological descriptions are hindering clear taxa circumscription and phylogenetic inferences within this group. In this work, we tested DNA sequence variation of three highly variable plastid markers and High-Throughput Sequencing (HTS) of the nuclear ribosomal 5S-IGS region in an attempt to decipher species identity within the V. tricolor species complex and to obtain an insight on their genome organization and evolution. Our results document the close relationships within this species group, a reliable molecular resolution for V. tricolor, and the common ancestry of V. arvensis and the poorly differentiated V. kitaibeliana and V. hymettia. Evidence of an important inter-population geographical divergence was recorded in V. tricolor and V. arvensis, pointing at the existence of different eco-cytotypes within these entities. Overall diversity patterns and the occurrence of two to four differently diverging 5S-IGS lineages are discussed in the light of the acknowledged taxonomy and genomic evolutive trajectories of sect. Melanium.

An integrated study of Violae Herba (Viola philippica) and five adulterants by morphology, chemical compositions and chloroplast genomes: insights into its certified plant origin

Background

Viola philippica Cav. is the only original plant for Violae Herba, as described in the Chinese Pharmacopoeia. The quality of this crude drug is affected by several adulterants from congeneric Viola species, and the authentic plant origin of Violae Herba is still controversial. Genome-based identification offers abundant genetic information and potential molecular markers that can be used for the authentication of closely related species. This study aims to investigate the certified origin of Violae Herba and to develop more effective markers for these easily confused species at the genetic level.

Methods

We compared the morphology and chemical composition of 18 batches of commercial samples and six widespread medicinal Viola plants used as Violae Herba or its substitutes by TLC and HPLC-Triple-TOF–MS/MS analyses. The complete chloroplast genomes of these species were sequenced and analyzed, including the general features, repeat sequences, mutational hotspots and phylogeny. The complete chloroplast genomes used as superbarcodes and some specific barcodes screened from mutational hotspots were tested for their ability to distinguish Viola species.

Results

A comparative study showed that Violae Herba is a multi-origin traditional Chinese medicine. Commercial decoction pieces and the standard reference drug were mainly derived from V. prionantha, clashing with the record in the Chinese Pharmacopoeia. Chloroplast genome analyses of V. philippica and five adulterants indicated that sequence divergence was relatively low within Viola species. By tree-based approaches, the complete chloroplast genomes showed a better discrimination ability and phylogenetic resolution for each Viola species. These results indicate that the whole chloroplast genomes can be used as superbarcodes to differentiate Viola medicinal plants. More specific DNA barcodes could be further developed from the Viola chloroplast genomes for more efficient and rapid identification of commercial Violae Herba and its adulterants.

Conclusions

This study has implications for chloroplast genome-based phylogenetic analysis and the authentication of multiple Viola species used as Violae Herba. The legal origin recorded in the Chinese Pharmacopoeia should be further revised to V. prionantha, in line with the commercial Violae Herba in the TCM markets.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13020-022-00585-9.

Classical Morphometrics in V. arvensis and V. kitaibeliana (Viola sect. Melanium) Reveals Intraspecific Variation with Implications for Species Delimitation: Inferences from a Case Study in Central Italy

The high morphological variability of Viola arvensis may hinder the proper identification of the closely related species with an implication for biodiversity surveys. Variation in floral and vegetative morphology was explored in V. arvensis, compared to V. kitaibeliana, based upon 14 wild Italian populations, to provide new insights into their diagnostic features. Species were characterized using 32 morphological descriptors assessed on 272 flowers and as many leaves and scored as quantitative and categorical variables. Statistical methods, including Linear Discriminant Analysis (LDA), were applied to test species delimitation. Data highlighted variations in sepal size, petal size, leaves shape, stylar dark spot, and pollen magazine higher within V. arvensis than between V. arvensis and V. kitaibeliana. LDA partitioned the V. arvensis samples into two distinct clusters; no clear distinction was found between the cluster combining individuals from grasslands and V. kitaibeliana. The separation of V. arvensis and V. kitaibeliana from V. tricolor, included as a reference, was noticeable. Correlations were found in all species between the flower/leaf position on the stem and some floral and vegetative features. The shape and margin of the lower sepal, the stylar flap, and the lamina margin and apex were diagnostic in field identification. The results support the recognition of an undescribed V. arvensis eco-phenotype linked to seminatural dry grasslands, easily distinguishable from the field-grown type of V. arvensis but hardly distinguishable from the dwarf pansy. Data further corroborate the assumption of general deep-rooted confusion in ascribing poorly developed individuals of V. arvensis to the rare and locally threatened V. kitaibeliana, leading to potential implications for its conservation.

Merging Arcs to Produce Acyclic Phylogenetic Networks and Normal Networks

As phylogenetic networks grow increasingly complicated, systematic methods for simplifying them to reveal properties will become more useful. This paper considers how to modify acyclic phylogenetic networks into other acyclic networks by contracting specific arcs that include a set D. The networks need not be binary, so vertices in the networks may have more than two parents and/or more than two children. In general, in order to make the resulting network acyclic, additional arcs not in D must also be contracted. This paper shows how to choose D so that the resulting acyclic network is “pre-normal”. As a result, removal of all redundant arcs yields a normal network. The set D can be selected based only on the geometry of the network, giving a well-defined normal phylogenetic network depending only on the given network. There are CSD maps relating most of the networks. The resulting network can be visualized as a “wired lift” in the original network, which appears as the original network with each arc drawn in one of three ways.

The complete chloroplast genome sequence of Viola verecunda (Violaceae)

Viola verecunda is a perennial plant native to the mountainous areas of Northeast Asia. Here the complete chloroplast genome of V. verecunda and its phylogenetic relationships to other species within the genus Viola are reported. The complete chloroplast genome of V. verecunda is 157,843 bp in length and circular in structure with four regions: a large single-copy region (86,345 bp), a small single-copy region (17,292 bp), and a pair of inverted repeat regions (27,103 bp each). The chloroplast genome contains 111 unique genes comprising 77 protein-coding, 30 unique tRNA, and 4 unique rRNA genes. Based on the protein-coding gene sequences from eight Viola chloroplast genomes, with Balanops balansae designated as the outgroup, maximum likelihood tree analysis indicates that V. verecunda is more closely related to V. raddeana than to other Viola species. The complete chloroplast genome of V. verecunda contributes to a better understanding of the phylogenetic relationships among Viola species.

Seasonal and Simultaneous Cleistogamy in Rostrate Violets (Viola, subsect. Rostratae, Violaceae)

The special mixed reproductive system, i.e., the ability of an individual plant to develop both open, chasmogamous (CH) flowers adapted to cross-pollination and closed, cleistogamous (CL) flowers with obligate self-pollinating, is a common phenomenon in Viola L. In most sections of Northern Hemisphere violets, cleistogamy is seasonal, and CH and CL flowers develop sequentially in the season. Non-seasonal cleistogamy (simultaneous) is a rare phenomenon in rostrate violets. In the current study, we focused on modification of the CH/CL mating system in V. caspia by environmental conditions, resulting in a gradual switch from temporal cleistogamy, occurring in nature, to simultaneous cleistogamy under greenhouse conditions. V. reichenbachiana with seasonal cleistogamy was the control for V. caspia with the labile seasonal/simultaneous cleistogamy system. In simultaneous cleistogamy, the CH and CL flowers, fruits and seeds developed on an individual plant at the same time on the same branch. The typical difference between CH and CL flowers’ pistils is a straight style ending with a head-like stigma in CH and a curved style in CL adapted to self-pollination. This trait persists in the fruit and seed stages, allowing for easy recognition of fruit of CL and CH flowers in simultaneous cleistogamy. Floral meristems of CH flowers of V. reichenbachiana developed on the rhizome at the end of the growing season under short-day conditions and remained dormant until the following season. The CL floral meristems formed under long-day conditions on elongating lateral branches in the upper leaf axils. The daily temperature influenced the variable CH/CL ratio of V. caspia in nature and greenhouse conditions. Regulation of the CL/CH flower ratio by modifying environmental factors is important for basic research on genetic/epigenetic regulation of cleistogamy and for practical use to produce genetically stable lines of economically important species via CL seeds.

Characterization and Dynamics of Intracellular Gene Transfer in Plastid Genomes of Viola (Violaceae) and Order Malpighiales

Functional gene transfer from organelles to the nucleus, known as intracellular gene transfer (IGT), is an ongoing process in flowering plants. The complete plastid genomes (plastomes) of two Ulleung island endemic violets, Viola ulleungdoensis and V. woosanensis, were characterized, revealing a lack of the plastid-encoded infA, rpl32, and rps16 genes. In addition, functional replacement of the three plastid-encoded genes in the nucleus was confirmed within the genus Viola and the order Malpighiales. Three strategies for the acquisition of a novel transit peptide for successful IGT were identified in the genus Viola. Nuclear INFA acquired a novel transit peptide with very low identity between these proteins, whereas the nuclear RPL32 gene acquired an existing transit peptide via fusion with the nuclear-encoded plastid-targeted SOD gene (Cu-Zn superoxide dismutase superfamily) as one exon, and translated both proteins in the cytosol using alternative mRNA splicing. Nuclear RPS16 contains an internal transit peptide without an N-terminal extension. Gene loss or pseudogenization of the plastid-borne rpl32 and rps16 loci was inferred to occur in the common ancestor of the genus Viola based on an infrageneric phylogenetic framework in Korea. Although infA was lost in the common ancestor of the order Malpighiales, the rpl32 and rps16 genes were lost multiple times independently within the order. Our current study sheds additional light on plastid genome composition and IGT mechanisms in the violet genus and in the order Malpighiales.

HybPhaser: A workflow for the detection and phasing of hybrids in target capture data sets

PREMISE

Hybrids contain divergent alleles that can confound phylogenetic analyses but can provide insights into reticulated evolution when identified and phased. We developed a workflow to detect hybrids in target capture data sets and phase reads into parental lineages using a similarity and phylogenetic framework.

METHODS

We used Angiosperms353 target capture data for Nepenthes, including known hybrids to test the novel workflow. Reference mapping was used to assess heterozygous sites across the data set and to detect hybrid accessions and paralogous genes. Hybrid samples were phased by mapping reads to multiple references and sorting reads according to similarity. Phased accessions were included in the phylogenetic framework.

RESULTS

All known Nepenthes hybrids and nine additional samples had high levels of heterozygous sites, had reads associated with multiple divergent clades, and were phased into accessions resembling divergent haplotypes. Phylogenetic analysis including phased accessions increased clade support and confirmed parental lineages of hybrids.

DISCUSSION

HybPhaser provides a novel approach to detect and phase hybrids in target capture data sets, which can provide insights into reticulations by revealing origins of hybrids and reduce conflicting signal, leading to more robust phylogenetic analyses.

Cyclotide host-defense tailored for species and environments in violets from the Canary Islands

Cyclotides are cyclic peptides produced by plants. Due to their insecticidal properties, they are thought to be involved in host defense. Violets produce complex mixtures of cyclotides, that are characteristic for each species and variable in different environments. Herein, we utilized mass spectrometry (LC–MS, MALDI-MS), transcriptomics and biological assays to investigate the diversity, differences in cyclotide expression based on species and different environment, and antimicrobial activity of cyclotides found in violets from the Canary Islands. A wide range of different habitats can be found on these islands, from subtropical forests to dry volcano peaks at high altitudes. The islands are inhabited by the endemic Viola palmensis, V. cheiranthifolia, V. anagae and the common V. odorata. The number of cyclotides produced by a given species varied in plants from different environments. The highest diversity was noted in V. anagae which resides in subtropical forest and the lowest in V. cheiranthifolia from the Teide volcano. Transcriptome sequencing and LC–MS were used to identify 23 cyclotide sequences from V. anagae. Cyclotide extracts exhibited antifungal activities with the lowest minimal inhibitory concentrations noted for V. anagae (15.62 μg/ml against Fusarium culmorum). The analysis of the relative abundance of 30 selected cyclotides revealed patterns characteristic to both species and populations, which can be the result of genetic variability or environmental conditions in different habitats. The current study exemplifies how plants tailor their host defense peptides for various habitats, and the usefulness of cyclotides as markers for chemosystematics.

Normalising phylogenetic networks

Rooted phylogenetic networks provide a way to describe species' relationships when evolution departs from the simple model of a tree. However, networks inferred from genomic data can be highly tangled, making it difficult to discern the main reticulation signals present. In this paper, we describe a natural way to transform any rooted phylogenetic network into a simpler canonical network, which has desirable mathematical and computational properties, and is based only on the 'visible' vertices in the original network. The method has been implemented and we demonstrate its application to some examples.

A new species of Viola (Violaceae) from Guangdong Province, China

Viola huizhouensis (Violaceae), a new species from Xiangtoushan National Nature Reserve of Guangdong Province in China, is described and illustrated. The new species is most similar to V. guangzhouensis, but it can be easily distinguished by its much stouter rhizome, lack of aerial stem, dense pubescence of the basal pedicel and the whole plant. Our phylogenetic analysis, based on ITS sequences, confirms that the new species belongs to V. sect. Diffusae.

Resolving phylogeny and polyploid parentage using genus-wide genome-wide sequence data from birch trees

Numerous plant genera have a history including frequent hybridisation and polyploidisation (allopolyploidisation), which means that their phylogeny is a network of reticulate evolution that cannot be accurately depicted as a bifurcating tree with a single tip per species. The genus Betula, which contains many ecologically important tree species, is a case in point. We generated genome-wide sequence reads for 27 diploid and 36 polyploid Betula species or subspecies using restriction site associated DNA (RAD) sequences. These reads were assembled into contigs with a mean length of 675 bp. We reconstructed the evolutionary relationships among diploid Betula species using both supermatrix (concatenation) and species tree methods. We identified the closest diploid relatives of the polyploids according to the relative rates at which reads from polyploids mapped to contigs from different diploid species within a concatenated reference sequence. By mapping reads from allopolyploids to their different putative diploid relatives we assembled contigs from the putative sub-genomes of allopolyploid taxa. We used these to build new phylogenies that included allopolyploid sub-genomes as separate tips. This approach yielded a highly evidenced phylogenetic hypothesis for the genus Betula, including the complex reticulate origins of the majority of its polyploid taxa. Our phylogeny divides the genus into two well supported clades, which, interestingly, differ in their seed-wing morphology. We therefore propose to split Betula into two subgenera.

Defining phylogenetic networks using ancestral profiles

Rooted phylogenetic networks provide a more complete representation of the ancestral relationship between species than phylogenetic trees when reticulate evolutionary processes are at play. One way to reconstruct a phylogenetic network is to consider its 'ancestral profile' (the number of paths from each ancestral vertex to each leaf). In general, this information does not uniquely determine the underlying phylogenetic network. A recent paper considered a new class of phylogenetic networks called 'orchard networks' where this uniqueness was claimed to hold. Here we show that an additional restriction on the network, that of being 'stack-free', is required in order for the original uniqueness claim to hold. On the other hand, if the additional stack-free restriction is lifted, we establish an alternative result; namely, there is uniqueness within the class of orchard networks up to the resolution of vertices of high in-degree.

Detecting destabilizing species in the phylogenetic backbone of Potentilla (Rosaceae) using low-copy nuclear markers

The genus Potentilla (Rosaceae) has been subjected to several phylogenetic studies, but resolving its evolutionary history has proven challenging. Previous analyses recovered six, informally named, groups: the Argentea, Ivesioid, Fragarioides, Reptans, Alba and Anserina clades, but the relationships among some of these clades differ between data sets. The Reptans clade, which includes the type species of Potentilla, has been noticed to shift position between plastid and nuclear ribosomal data sets. We studied this incongruence by analysing four low-copy nuclear markers, in addition to chloroplast and nuclear ribosomal data, with a set of Bayesian phylogenetic and Multispecies Coalescent (MSC) analyses. A selective taxon removal strategy demonstrated that the included representatives from the Fragarioides clade, P. dickinsii and P. fragarioides, were the main sources of the instability seen in the trees. The Fragarioides species showed different relationships in each gene tree, and were only supported as a monophyletic group in a single marker when the Reptans clade was excluded from the analysis. The incongruences could not be explained by allopolyploidy, but rather by homoploid hybridization, incomplete lineage sorting or taxon sampling effects. When P. dickinsii and P. fragarioides were removed from the data set, a fully resolved, supported backbone phylogeny of Potentilla was obtained in the MSC analysis. Additionally, indications of autopolyploid origins of the Reptans and Ivesioid clades were discovered in the low-copy gene trees.

Implications of plastome evolution in the true lilies (monocot order Liliales)

The families of the monocot order Liliales exhibit highly contrasting characteristic of photosynthetic and mycoheterotrophic life histories. Although previous phylogenetic and morphological studies of Liliales have been conducted, they have not examined molecular evolution associated with this contrasting phenomenon. Here, we conduct the first comparative plastome study of all ten families of Liliales using 29 newly sequenced plastid genomes analyzed together with previously published data. We also present a phylogenetic analysis for Liliales of 78 plastid genes combined with 22 genes from all three genomes (nuclear 18S rDNA and phyC; 17 plastid genes; and mitochondrial matR, atpA, and cob). Within the newly generated phylogenetic tree of Liliales, we evaluate the ancestral state changes of selected morphological traits in the order. There are no significant differences in plastid genome features among species that show divergent characteristics correlated with family circumscriptions. However, the results clearly differentiate between photosynthetic and mycoheterotrophic taxa of Liliales in terms of genome structure, and gene content and order. The newly sequenced plastid genomes and combined three-genome data revealed Smilacaceae as sister to Liliaceae instead of Philesiaceae and Ripogonaceae. Additionally, we propose a revised familial classification system of Liliales that consists of nine families, considering Ripogonaceae a synonym of Philesiaceae. The ancestral state reconstruction indicated synapomorphies for each family of Liliales, except Liliaceae, Melanthiaceae and Colchicaceae. A taxonomic key for all nine families of Liliales is also provided.

Complex patterns of reticulate evolution in opportunistic weeds (Potentilla L., Rosaceae), as revealed by low-copy nuclear markers

Background

Most cinquefoils (Potentilla L., Rosaceae) are polyploids, ranging from tetraploid (4x) to dodecaploid (12x), diploids being a rare exception. Previous studies based on ribosomal and chloroplast data indicated that Norwegian cinquefoil (P. norvegica L.) has genetic material from two separate clades within Potentilla; the Argentea and the Ivesioid clades – and thus a possible history of hybridization and polyploidization (allopolyploidy). In order to trace the putative allopolyploid origin of the species, sequence data from low-copy, biparentally inherited, nuclear markers were used. Specimens covering the circumpolar distribution of P. norvegica and its two subspecies were included, along with the morphologically similar P. intermedia. Potentilla species of low ploidy level known to belong to other relevant clades were also included.

Results

Gene trees based on three low-copy nuclear markers, obtained by Bayesian Inference and Maximum Likelihood analyses, showed slightly different topologies. This is likely due to genomic reorganizations following genome duplication, but the gene trees were not in conflict with a species tree of presumably diploid taxa obtained by Multispecies Coalescent analysis. The results show that both P. norvegica and P. intermedia are allopolyploids with a shared evolutionary history involving at least four parental lineages, three from the Argentea clade and one from the Ivesioid clade.

Conclusions

This is the first time that reticulate evolution has been proven in the genus Potentilla, and shows the importance of continuing working with low-copy markers in order to properly resolve its evolutionary history. Several hybridization events between the Argentea and Ivesioid clades may have given rise to the species of Wolf’s grex Rivales. To better estimate when and where these hybridizations occurred, other Argentea, Ivesioid and Rivales species should be included in future studies.

Hybridization preceded radiation in diploid wheats

Evolutionary relationships among the Aegilops-Triticum relatives of cultivated wheats have been difficult to resolve owing to incomplete lineage sorting and reticulate evolution. Recent studies have suggested that the wheat D-genome lineage (progenitor of Ae. tauschii) originated through homoploid hybridization between the A-genome lineage (progenitor of Triticum s.str.) and the B-genome lineage (progenitor of Ae. speltoides). This scenario of reticulation has been debated, calling for adequate phylogenetic analyses based on comprehensive sampling. To reconstruct the evolution of Aegilops-Triticum diploids, we here combined high-throughput sequencing of 38 nuclear low-copy loci of multiple accessions of all 13 species with inferences of the species phylogeny using the full-parameterized MCMC_SEQ method. Phylogenies recovered a monophyletic Aegilops-Triticum lineage that began diversifying ~6.6 Ma ago and gave rise to four sublineages, i.e. the A- (2 species), B- (1 species), D- (9 species) and T- (Ae. mutica) genome lineage. Full-parameterized phylogenies as well as patterns of tree dilation and tree compression supported a hybrid origin of the D-genome lineage from A and B ~3.0-4.0 Ma ago, and did not indicate additional hybridization events. Conflicting ABBA-BABA tests suggestive of further reticulation were shown here to result from ancestral population structure rather than hybridization. This comprehensive and dated phylogeny of wheat relatives indicates that the origin of the hybrid D-genome was followed by intense diversification into the majority of extant diploid as well as allopolyploid wild wheats.

The complete chloroplast genome sequences of four Viola species (Violaceae) and comparative analyses with its congeneric species

We report the complete chloroplast genomes of four Viola species (V. mirabilis, V. phalacrocarpa, V. raddeana, and V. websteri) and the results of a comparative analysis between these species and the published plastid genome of the congeneric species V. seoulensis. The total genome length of the five Viola species, including the four species analyzed in this study and the species analyzed in the previous study, ranged from 156,507 (V. seoulensis) to 158,162 bp (V. mirabilis). The overall GC contents of the genomes were almost identical (36.2-36.3%). The five Viola plastomes each contained 111 unique genes comprising 77 protein-coding genes, 30 transfer RNA (tRNA) genes, and 4 ribosomal RNA (rRNA) genes. Among the annotated genes, 16 contained one or two introns. Based on the results of a chloroplast genome structure comparison using MAUVE, all five Viola plastomes were almost identical. Additionally, the large single copy (LSC), inverted repeat (IR), and small single copy (SSC) junction regions were conserved among the Viola species. A total of 259 exon, intron, and intergenic spacer (IGS) fragments were compared to verify the divergence hotspot regions. The nucleotide diversity (Pi) values ranged from 0 to 0.7544. The IR region was relatively more conserved than the LSC and SSC regions. The Pi values in ten noncoding regions were relatively high (>0.03). Among these regions, all but rps19-trnH, petG-trnW, rpl16-rps3, and rpl2-rpl23 represent useful molecular markers for phylogenetic studies and will be helpful to resolve the phylogenetic relationships of Viola. The phylogenetic tree, which used 76 protein-coding genes from 21 Malpighiales species and one outgroup species (Averrhoa carambola), revealed that Malpighiales is divided into five clades at the family level: Erythroxylaceae, Chrysobalanaceae, Euphorbiaceae, Salicaceae, and Violaceae. Additionally, Violaceae was monophyletic, with a bootstrap value of 100% and was divided into two subclades.

Extensive allopolyploidy in the neotropical genus Lachemilla (Rosaceae) revealed by PCR-based target enrichment of the nuclear ribosomal DNA cistron and plastid phylogenomics

PREMISE OF THE STUDY

Polyploidy has been long recognized as an important force in plant evolution. Previous studies had suggested widespread occurrence of polyploidy and the allopolyploid origin of several species in the diverse neotropical genus Lachemilla (Rosaceae). Nonetheless, this evidence has relied mostly on patterns of cytonuclear discordance, and direct evidence from nuclear allelic markers is still needed.

METHODS

Here we used PCR target enrichment in combination with high throughput sequencing to obtain multiple copies of the nuclear ribosomal (nr) DNA cistron and 45 regions of the plastid genome (cpDNA) from 219 accessions representing 48 species of Lachemilla and to explore the allopolyploid origin of species in this group.

KEY RESULTS

We were able to identify multiple nrDNA ribotypes and establish clear evidence of allopolyploidy in 33 species of Lachemilla, showing that this condition is common and widespread in the genus. Additionally, we found evidence for three autopolyploid species. We also established multiple, independent origins of several allopolyploid species. Finally, based solely on the cpDNA phylogeny, we identified that the monotypic genus Farinopsis is the sister group of Lachemilla and allied genera within subtribe Fragariinae.

CONCLUSIONS

Our study demonstrates the utility of the nuclear ribosomal DNA cistron to detect allopolyploidy when concerted evolution of this region is not complete. Additionally, with a robust chloroplast phylogeny in place, the direction of hybridization events can be established, and multiple, independent origins of allopolyploid species can be identified.

Transcriptomics Identifies Modules of Differentially Expressed Genes and Novel Cyclotides in Viola pubescens

Viola is a large genus with worldwide distribution and many traits not currently exemplified in model plants including unique breeding systems and the production of cyclotides. Here we report de novo genome assembly and transcriptomic analyses of the non-model species Viola pubescens using short-read DNA sequencing data and RNA-Seq from eight diverse tissues. First, V. pubescens genome size was estimated through flow cytometry, resulting in an approximate haploid genome of 455 Mbp. Next, the draft V. pubescens genome was sequenced and assembled resulting in 264,035,065 read pairs and 161,038 contigs with an N50 length of 3,455 base pairs (bp). RNA-Seq data were then assembled into tissue-specific transcripts. Together, the DNA and transcript data generated 38,081 ab initio gene models which were functionally annotated based on homology to Arabidopsis thaliana genes and Pfam domains. Gene expression was visualized for each tissue via principal component analysis and hierarchical clustering, and gene co-expression analysis identified 20 modules of tissue-specific transcriptional networks. Some of these modules highlight genetic differences between chasmogamous and cleistogamous flowers and may provide insight into V. pubescens' mixed breeding system. Orthologous clustering with the proteomes of A. thaliana and Populus trichocarpa revealed 8,531 sequences unique to V. pubescens, including 81 novel cyclotide precursor sequences. Cyclotides are plant peptides characterized by a stable, cyclic cystine knot motif, making them strong candidates for drug scaffolding and protein engineering. Analysis of the RNA-Seq data for these cyclotide transcripts revealed diverse expression patterns both between transcripts and tissues. The diversity of these cyclotides was also highlighted in a maximum likelihood protein cladogram containing V. pubescens cyclotides and published cyclotide sequences from other Violaceae and Rubiaceae species. Collectively, this work provides the most comprehensive sequence resource for Viola, offers valuable transcriptomic insight into V. pubescens, and will facilitate future functional genomics research in Viola and other diverse plant groups.

Embracing heterogeneity: coalescing the Tree of Life and the future of phylogenomics

Building the Tree of Life (ToL) is a major challenge of modern biology, requiring advances in cyberinfrastructure, data collection, theory, and more. Here, we argue that phylogenomics stands to benefit by embracing the many heterogeneous genomic signals emerging from the first decade of large-scale phylogenetic analysis spawned by high-throughput sequencing (HTS). Such signals include those most commonly encountered in phylogenomic datasets, such as incomplete lineage sorting, but also those reticulate processes emerging with greater frequency, such as recombination and introgression. Here we focus specifically on how phylogenetic methods can accommodate the heterogeneity incurred by such population genetic processes; we do not discuss phylogenetic methods that ignore such processes, such as concatenation or supermatrix approaches or supertrees. We suggest that methods of data acquisition and the types of markers used in phylogenomics will remain restricted until a posteriori methods of marker choice are made possible with routine whole-genome sequencing of taxa of interest. We discuss limitations and potential extensions of a model supporting innovation in phylogenomics today, the multispecies coalescent model (MSC). Macroevolutionary models that use phylogenies, such as character mapping, often ignore the heterogeneity on which building phylogenies increasingly rely and suggest that assimilating such heterogeneity is an important goal moving forward. Finally, we argue that an integrative cyberinfrastructure linking all steps of the process of building the ToL, from specimen acquisition in the field to publication and tracking of phylogenomic data, as well as a culture that values contributors at each step, are essential for progress.

Widespread ancient whole-genome duplications in Malpighiales coincide with Eocene global climatic upheaval

Whole-genome duplications (WGDs) are widespread and prevalent in vascular plants and frequently coincide with major episodes of global and climatic upheaval, including the mass extinction at the Cretaceous-Tertiary boundary (c. 65 Ma) and during more recent periods of global aridification in the Miocene (c. 10-5 Ma). Here, we explore WGDs in the diverse flowering plant clade Malpighiales. Using transcriptomes and complete genomes from 42 species, we applied a multipronged phylogenomic pipeline to identify, locate, and determine the age of WGDs in Malpighiales using three means of inference: distributions of synonymous substitutions per synonymous site (Ks ) among paralogs, phylogenomic (gene tree) reconciliation, and a likelihood-based gene-count method. We conservatively identify 22 ancient WGDs, widely distributed across Malpighiales subclades. Importantly, these events are clustered around the Eocene-Paleocene transition (c. 54 Ma), during which time the planet was warmer and wetter than any period in the Cenozoic. These results establish that the Eocene Climatic Optimum likely represents a previously unrecognized period of prolific WGDs in plants, and lends further support to the hypothesis that polyploidization promotes adaptation and enhances plant survival during episodes of global change, especially for tropical organisms like Malpighiales, which have tight thermal tolerances.

Dating the Species Network: Allopolyploidy and Repetitive DNA Evolution in American Daisies (Melampodium sect. Melampodium, Asteraceae)

Allopolyploidy has played an important role in the evolution of the flowering plants. Genome mergers are often accompanied by significant and rapid alterations of genome size and structure via chromosomal rearrangements and altered dynamics of tandem and dispersed repetitive DNA families. Recent developments in sequencing technologies and bioinformatic methods allow for a comprehensive investigation of the repetitive component of plant genomes. Interpretation of evolutionary dynamics following allopolyploidization requires both the knowledge of parentage and the age of origin of an allopolyploid. Whereas parentage is typically inferred from cytogenetic and phylogenetic data, age inference is hampered by the reticulate nature of the phylogenetic relationships. Treating subgenomes of allopolyploids as if they belonged to different species (i.e., no recombination among subgenomes) and applying cross-bracing (i.e., putting a constraint on the age difference of nodes pertaining to the same event), we can infer the age of allopolyploids within the framework of the multispecies coalescent within BEAST2. Together with a comprehensive characterization of the repetitive DNA fraction using the RepeatExplorer pipeline, we apply the dating approach in a group of closely related allopolyploids and their progenitor species in the plant genus Melampodium (Asteraceae). We dated the origin of both the allotetraploid, Melampodium strigosum, and its two allohexaploid derivatives, Melampodium pringlei and Melampodium sericeum, which share both parentage and the direction of the cross, to the Pleistocene ($<$1.4 Ma). Thus, Pleistocene climatic fluctuations may have triggered formation of allopolyploids possibly in short intervals, contributing to difficulties in inferring the precise temporal order of allopolyploid species divergence of M. sericeum and M. pringlei. The relatively recent origin of the allopolyploids likely played a role in the near-absence of major changes in the repetitive fraction of the polyploids' genomes. The repetitive elements most affected by the postpolyploidization changes represented retrotransposons of the Ty1-copia lineage Maximus and, to a lesser extent, also Athila elements of Ty3-gypsy family.

Reconstructing the origins and the biogeography of species' genomes in the highly reticulate allopolyploid-rich model grass genus Brachypodium using minimum evolution, coalescence and maximum likelihood approaches

The identification of homeologous genomes and the biogeographical analyses of highly reticulate allopolyploid-rich groups face the challenge of incorrectly inferring the genomic origins and the biogeographical patterns of the polyploids from unreliable strictly bifurcating trees. Here we reconstruct a plausible evolutionary scenario of the diverging and merging genomes inherited by the diploid and allopolyploid species and cytotypes of the model grass genus Brachypodium. We have identified the ancestral Brachypodium genomes and inferred the paleogeographical ranges for potential hybridization events that originated its allopolyploid taxa. We also constructed a comprehensive phylogeny of Brachypodium from five nuclear and plastid genes using Species Tree Minimum Evolution allele grafting and Species Network analysis. The divergence ages of the lineages were estimated from a consensus maximum clade credibility tree using fossil calibrations, whereas ages of origin of the diploid and allopolyploid species were inferred from coalescence Bayesian methods. The biogeographical events of the genomes were reconstructed using a stratified Dispersal-Extinction-Colonization model with three temporal windows. Our combined Minimum Evolution-coalescence-Bayesian approach allowed us to infer the origins and the identities of the homeologous genomes of the Brachypodium allopolyploids, matching the expected ploidy levels of the hybrids. To date, the current extant progenitor genomes (species) are only known for B. hybridum. Putative ancestral homeologous genome have been inherited by B. mexicanum, ancestral and recent genomes by B. boissieri, and only recently evolved genomes by B. retusum and the core perennial clade allopolyploids (B. phoenicoides, B. pinnatum 4x, B. rupestre 4x). We dissected the complex spatio-temporal evolution of ancestral and recent genomes and have detected successive splitting, dispersal and merging events for dysploid homeologous genomes in diverse geographical scenarios that have led to the current extant taxa. Our data support Mid-Miocene splits of the Holarctic ancestral genomes that preceded the Late Miocene origins of Brachypodium ancestors of the modern diploid species. Successive divergences of the annual B. stacei and B. distachyon diploid genomes were implied to have occurred in the Mediterranean region during the Late Miocene-Pliocene. By contrast, a profusion of splits, range expansions and different genome mergings were inferred for the perennial diploid genomes in the Mediterranean and Eurasian regions, with sporadic colonizations and further mergings in other continents during the Quaternary. A reliable biogeographical scenario was obtained for the Brachypodium genomes and allopolyploids where homeologous genomes split from their respective diploid counterpart lineages in the same ancestral areas, showing similar or distinct dispersals. By contrast, the allopolyploid taxa remained in the same ancestral ranges after hybridization and genome doubling events. Our approach should have utility in deciphering the genomic composition and the historical biogeography of other allopolyploid-rich organismal groups, which are predominant in eukaryotes.

Lack of correlation between pollen aperture number and environmental factors in pansies (Viola L., sect. Melanium Ging.) - pollen heteromorphism re-examined

The development of different pollen morphs by one specimen - pollen heteromorphism - occurs in ca. one-third species of the genus Viola. Melanium section species (pansies) stand out in producing the widest range of apertures among Viola species. Aperture number decreases with elevation increase, and faster germination of five-aperturate pollen as compared with three-aperturate has previously been postulated. We re-examined pollen heteromorphism in the context of its viability, and made correlations with elevation (>1500 m a.s.l. versus <1500 m a.s.l.), soil type (metalliferous versus non-metalliferous; MET versus NMET) and chromosome number based on selected study criteria of ca. 20% karyologically and morphologically strongly differentiated but genetically closely related pansies. A total of 79% of analysed species were heteromorphic, forming three- to six-colp or ate pollen per individual flower. Mean aperture number and pollen viability were not affected by soil type (MET versus NMET). Mean aperture number was also not influenced by elevation or species chromosome number. Positive correlations were established between aperture number and pollen viability, negative between pollen viability and elevation (increasing altitude of 100 m decreased pollen viability by 0.4%) and lack of correlation between chromosome number and pollen viability. The varied frequencies of different pollen morphs among species are not under the general pressure of ecological conditions, as previously postulated for the species of Melanium section. Rather, this trait in pansies, similar to other floral characters (e.g. long, curved nectar spur, 'landing platform', posterior petals with nectar guides), is adaptive but dependent on the breeding system (inbreeding versus outbreeding) of the individual species.

Cladogenetic and Anagenetic Models of Chromosome Number Evolution: A Bayesian Model Averaging Approach

ABSSTRACT

Chromosome number is a key feature of the higher-order organization of the genome, and changes in chromosome number play a fundamental role in evolution. Dysploid gains and losses in chromosome number, as well as polyploidization events, may drive reproductive isolation and lineage diversification. The recent development of probabilistic models of chromosome number evolution in the groundbreaking work by Mayrose et al. (2010, ChromEvol) have enabled the inference of ancestral chromosome numbers over molecular phylogenies and generated new interest in studying the role of chromosome changes in evolution. However, the ChromEvol approach assumes all changes occur anagenetically (along branches), and does not model events that are specifically cladogenetic. Cladogenetic changes may be expected if chromosome changes result in reproductive isolation. Here we present a new class of models of chromosome number evolution (called ChromoSSE) that incorporate both anagenetic and cladogenetic change. The ChromoSSE models allow us to determine the mode of chromosome number evolution; is chromosome evolution occurring primarily within lineages, primarily at lineage splitting, or in clade-specific combinations of both? Furthermore, we can estimate the location and timing of possible chromosome speciation events over the phylogeny. We implemented ChromoSSE in a Bayesian statistical framework, specifically in the software RevBayes, to accommodate uncertainty in parameter estimates while leveraging the full power of likelihood based methods. We tested ChromoSSE's accuracy with simulations and re-examined chromosomal evolution in Aristolochia, Carex section Spirostachyae, Helianthus, Mimulus sensu lato (s.l.), and Primula section Aleuritia, finding evidence for clade-specific combinations of anagenetic and cladogenetic dysploid and polyploid modes of chromosome evolution. [Anagenetic; Bayes factors; chromosome evolution; chromosome speciation; chromoSSE; cladogenetic; dysploidy; phylogenetic models; polyploidy; reversible-jump Markov chain Monte Carlo; whole genome duplication.].

Gene-Tree Reconciliation with MUL-Trees to Resolve Polyploidy Events

Polyploidy can have a huge impact on the evolution of species, and it is a common occurrence, especially in plants. The two types of polyploids-autopolyploids and allopolyploids-differ in the level of divergence between the genes that are brought together in the new polyploid lineage. Because allopolyploids are formed via hybridization, the homoeologous copies of genes within them are at least as divergent as orthologs in the parental species that came together to form them. This means that common methods for estimating the parental lineages of allopolyploidy events are not accurate, and can lead to incorrect inferences about the number of gene duplications and losses. Here, we have adapted an algorithm for topology-based gene-tree reconciliation to work with multi-labeled trees (MUL-trees). By definition, MUL-trees have some tips with identical labels, which makes them a natural representation of the genomes of polyploids. Using this new reconciliation algorithm we can: accurately place allopolyploidy events on a phylogeny, identify the parental lineages that hybridized to form allopolyploids, distinguish between allo-, auto-, and (in most cases) no polyploidy, and correctly count the number of duplications and losses in a set of gene trees. We validate our method using gene trees simulated with and without polyploidy, and revisit the history of polyploidy in data from the clades including both baker's yeast and bread wheat. Our re-analysis of the yeast data confirms the allopolyploid origin and parental lineages previously identified for this group. The method presented here should find wide use in the growing number of genomes from species with a history of polyploidy. [Polyploidy; reconciliation; whole-genome duplication.].

Displayed Trees Do Not Determine Distinguishability Under the Network Multispecies Coalescent

Recent work in estimating species relationships from gene trees has included inferring networks assuming that past hybridization has occurred between species. Probabilistic models using the multispecies coalescent can be used in this framework for likelihood-based inference of both network topologies and parameters, including branch lengths and hybridization parameters. A difficulty for such methods is that it is not always clear whether, or to what extent, networks are identifiable-that is whether there could be two distinct networks that lead to the same distribution of gene trees. For cases in which incomplete lineage sorting occurs in addition to hybridization, we demonstrate a new representation of the species network likelihood that expresses the probability distribution of the gene tree topologies as a linear combination of gene tree distributions given a set of species trees. This representation makes it clear that in some cases in which two distinct networks give the same distribution of gene trees when sampling one allele per species, the two networks can be distinguished theoretically when multiple individuals are sampled per species. This result means that network identifiability is not only a function of the trees displayed by the networks but also depends on allele sampling within species. We additionally give an example in which two networks that display exactly the same trees can be distinguished from their gene trees even when there is only one lineage sampled per species. [gene tree, hybridization, identifiability, maximum likelihood, species tree, phylogeny.].

Improved transcriptome sampling pinpoints 26 ancient and more recent polyploidy events in Caryophyllales, including two allopolyploidy events

Studies of the macroevolutionary legacy of polyploidy are limited by an incomplete sampling of these events across the tree of life. To better locate and understand these events, we need comprehensive taxonomic sampling as well as homology inference methods that accurately reconstruct the frequency and location of gene duplications. We assembled a data set of transcriptomes and genomes from 168 species in Caryophyllales, of which 43 transcriptomes were newly generated for this study, representing one of the most densely sampled genomic-scale data sets available. We carried out phylogenomic analyses using a modified phylome strategy to reconstruct the species tree. We mapped the phylogenetic distribution of polyploidy events by both tree-based and distance-based methods, and explicitly tested scenarios for allopolyploidy. We identified 26 ancient and more recent polyploidy events distributed throughout Caryophyllales. Two of these events were inferred to be allopolyploidy. Through dense phylogenomic sampling, we show the propensity of polyploidy throughout the evolutionary history of Caryophyllales. We also provide a framework for utilizing transcriptome data to detect allopolyploidy, which is important as it may have different macroevolutionary implications compared with autopolyploidy.

Nonbinary Tree-Based Phylogenetic Networks

Rooted phylogenetic networks are used to describe evolutionary histories that contain non-treelike evolutionary events such as hybridization and horizontal gene transfer. In some cases, such histories can be described by a phylogenetic base-tree with additional linking arcs, which can, for example, represent gene transfer events. Such phylogenetic networks are called tree-based. Here, we consider two possible generalizations of this concept to nonbinary networks, which we call tree-based and strictly-tree-based nonbinary phylogenetic networks. We give simple graph-theoretic characterizations of tree-based and strictly-tree-based nonbinary phylogenetic networks. Moreover, we show for each of these two classes that it can be decided in polynomial time whether a given network is contained in the class. Our approach also provides a new view on tree-based binary phylogenetic networks. Finally, we discuss two examples of nonbinary phylogenetic networks in biology and show how our results can be applied to them.

Evaluating allopolyploid origins in strawberries (Fragaria) using haplotypes generated from target capture sequencing

BACKGROUND

Hybridization is observed in many eukaryotic lineages and can lead to the formation of polyploid species. The study of hybridization and polyploidization faces challenges both in data generation and in accounting for population-level phenomena such as coalescence processes in phylogenetic analysis. Genus Fragaria is one example of a set of plant taxa in which a range of ploidy levels is observed across species, but phylogenetic origins are unknown.

RESULTS

Here, using 20 diploid and polyploid Fragaria species, we combine approaches from NGS data analysis and phylogenetics to infer evolutionary origins of polyploid strawberries, taking into account coalescence processes. We generate haplotype sequences for 257 low-copy nuclear markers assembled from Illumina target capture sequence data. We then identify putative hybridization events by analyzing gene tree topologies, and further test predicted hybridizations in a coalescence framework. This approach confirms the allopolyploid ancestry of F. chiloensis and F. virginiana, and provides new allopolyploid ancestry hypotheses for F. iturupensis, F. moschata, and F. orientalis. Evidence of gene flow between diploids F. bucharica and F. vesca is also detected, suggesting that it might be appropriate to consider these groups as conspecifics.

CONCLUSIONS

This study is one of the first in which target capture sequencing followed by computational deconvolution of individual haplotypes is used for tracing origins of polyploid taxa. The study also provides new perspectives on the evolutionary history of Fragaria.

No evidence of contemporary interploidy gene flow between the closely related European woodland violets Viola reichenbachiana and V. riviniana (sect. Viola, Violaceae)

Viola reichenbachiana (2n = 4x = 20) and V. riviniana (2n = 8x = 40) are closely related species widely distributed in Europe, often sharing the same habitat throughout their overlapping ranges. It has been suggested in numerous studies that their high intraspecific morphological variability and plasticity might have been further increased by interspecific hybridisation in contact zones, given the sympatry of the species and the incomplete sterility of their hybrid. The aims of this study were to: (i) confirm that V. reichenbachiana and V. riviniana have one 4x genome in common, and (ii) determine the impact of hybridisation and introgression on genetic variation of these two species in selected European populations. For our study, we used 31 Viola populations from four European countries, which were analysed using AFLP and sequencing of a variable plastid intergenic spacer, trnH-psbA. Our analysis revealed that V. reichenbachiana exhibited larger haplotype diversity, having three species-specific haplotypes versus one in V. riviniana. The relationships among haplotypes suggest transfer of common haplotypes into V. riviniana from both V. reichenbachiana and hypothetically the other, now extinct, parental species. AFLP analysis showed low overall genetic diversity of both species, with V. riviniana showing higher among-population diversity. None of the morphologically designated hybrid populations had additive AFLP polymorphisms that would have indicated recent hybridisation. Also, kinship coefficients between both species did not indicate gene flow. V. riviniana showed significant population subdivision and significant isolation by distance, in contrast to V. reichenbachiana. The results indicate lack of gene flow between species, high influence of selfing on genetic variability, as well as probably only localised introgression toward V. riviniana.

Next-generation polyploid phylogenetics: rapid resolution of hybrid polyploid complexes using PacBio single-molecule sequencing

Difficulties in generating nuclear data for polyploids have impeded phylogenetic study of these groups. We describe a high-throughput protocol and an associated bioinformatics pipeline (Pipeline for Untangling Reticulate Complexes (Purc)) that is able to generate these data quickly and conveniently, and demonstrate its efficacy on accessions from the fern family Cystopteridaceae. We conclude with a demonstration of the downstream utility of these data by inferring a multi-labeled species tree for a subset of our accessions. We amplified four c. 1-kb-long nuclear loci and sequenced them in a parallel-tagged amplicon sequencing approach using the PacBio platform. Purc infers the final sequences from the raw reads via an iterative approach that corrects PCR and sequencing errors and removes PCR-mediated recombinant sequences (chimeras). We generated data for all gene copies (homeologs, paralogs, and segregating alleles) present in each of three sets of 50 mostly polyploid accessions, for four loci, in three PacBio runs (one run per set). From the raw sequencing reads, Purc was able to accurately infer the underlying sequences. This approach makes it easy and economical to study the phylogenetics of polyploids, and, in conjunction with recent analytical advances, facilitates investigation of broad patterns of polyploid evolution.

Cyclotide Evolution: Insights from the Analyses of Their Precursor Sequences, Structures and Distribution in Violets (Viola)

Cyclotides are a family of plant proteins that are characterized by a cyclic backbone and a knotted disulfide topology. Their cyclic cystine knot (CCK) motif makes them exceptionally resistant to thermal, chemical, and enzymatic degradation. By disrupting cell membranes, the cyclotides function as host defense peptides by exhibiting insecticidal, anthelmintic, antifouling, and molluscicidal activities. In this work, we provide the first insight into the evolution of this family of plant proteins by studying the Violaceae, in particular species of the genus Viola. We discovered 157 novel precursor sequences by the transcriptomic analysis of six Viola species: V. albida var. takahashii, V. mandshurica, V. orientalis, V. verecunda, V. acuminata, and V. canadensis. By combining these precursor sequences with the phylogenetic classification of Viola, we infer the distribution of cyclotides across 63% of the species in the genus (i.e., ~380 species). Using full precursor sequences from transcriptomes, we show an evolutionary link to the structural diversity of the cyclotides, and further classify the cyclotides by sequence signatures from the non-cyclotide domain. Also, transcriptomes were compared to cyclotide expression on a peptide level determined using liquid chromatography-mass spectrometry. Furthermore, the novel cyclotides discovered were associated with the emergence of new biological functions.

Genome origin, historical hybridization and genetic differentiation in Anthosachne australasica (Triticeae; Poaceae), inferred from chloroplast rbcL, trnH-psbA and nuclear Acc1 gene sequences

BACKGROUND AND AIMS

Anthosachne Steudel is a group of allopolyploid species that was derived from hexaploidization between the Asian StY genome Roegneria entity and the Australasia W genome Australopyrum species. Polyploidization and apomixis contribute to taxonomic complexity in Anthosachne Here, a study is presented on the phylogeny and evolutionary history of Anthosachne australasica The aims are to demonstrate the process of polyploidization events and to explore the differentiation patterns of the St genome following geographic isolation.

METHODS

Chloroplast rbcL and trnH-psbA and nuclear Acc1 gene sequences of 60 Anthosachne taxa and nine Roegneria species were analysed with those of 33 diploid taxa representing 20 basic genomes in Triticeae. The phylogenetic relationships were reconstructed. A time-calibrated phylogeny was generated to estimate the evolutionary history of A. australasica Nucleotide diversity patterns were used to assess the divergence within A. australasica and between Anthosachne and its putative progenitors.

KEY RESULTS

Three homoeologous copies of the Acc1 sequences from Anthosachne were grouped with the Acc1 sequences from Roegneria, Pseudoroegneria, Australopyrum, Dasypyrum and Peridictyon The chloroplast sequences of Anthosachne were clustered with those from Roegneria and Pseudoroegneria Divergence time for Anthosachne was dated to 4·66 million years ago (MYA). The level of nucleotide diversity in Australasian Anthosachne was higher than that in continental Roegneria A low level of genetic differentiation within the A. australasica complex was found.

CONCLUSIONS

Anthosachne originated from historical hybridization between Australopyrum species and a Roegneria entity colonized from Asia to Australasia via South-east Asia during the late Miocene. The St lineage served as the maternal donor during the speciation of Anthosachne A contrasting pattern of population genetic structure exists in the A. australasica complex. Greater diversity in island Anthosachne compared with continental Roegneria might be associated with mutation, polyploidization, apomixis and expansion. It is reasonable to consider that A. australasica var. scabra and A. australasica var. plurinervisa should be included in the A. australasica complex.

Immunolocalization of cyclotides in plant cells, tissues and organ supports their role in host defense

The distribution of cyclotides was visualized in plant cells, tissues and organs using immunohistochemistry. Finding of cyclotides in tissues potentially vulnerable to pathogen attacks supports their role as defense molecules. The cyclotide family of plant peptides is characterized by the cyclic cystine knot motif and its diverse biological activities. Given their insecticidal and antimicrobial properties, the role of cyclotides in planta is probably associated with host defense. Our current understanding of the cellular compartmentalization of cyclotides in the vacuole is based on indirect studies on transgenic model plants that do not express cyclotides naturally. Matrix-assisted laser desorption ionization (MALDI) imaging has also been used to study the distribution of cyclotides, but the technique's resolution was insufficient to determine their tissue or cell distribution. To avoid the limitations of these approaches, immunohistochemical visualization methods were used. Antibodies were raised in rabbits using cycloviolacin O2 (cyO2), and their specificity was determined by Western and dot blot experiments. Slides for immunohistochemical analysis were prepared from leaf, petiole and root fragments of Viola odorata and Viola uliginosa, and specimens were visualized using indirect epifluorescence microscopy. The antibodies against cyclotides were specific against selected bracelet cyclotides with high similarity (cyO2, cyO3, cyO8, cyO13) and suitable for immunohistochemistry. The tissue distribution of the cyclotides visualized in this way is consistent with their proposed role in host defense-relatively large quantities were observed in the leaf and petiole epidermis in both Viola species. Cyclotides were also found in vascular tissue in all the assessed plant organs. The vacuole storage of cyclotides was directly shown.

Immunolocalization of cyclotides in plant cells, tissues and organ supports their role in host defense

The distribution of cyclotides was visualized in plant cells, tissues and organs using immunohistochemistry. Finding of cyclotides in tissues potentially vulnerable to pathogen attacks supports their role as defense molecules. The cyclotide family of plant peptides is characterized by the cyclic cystine knot motif and its diverse biological activities. Given their insecticidal and antimicrobial properties, the role of cyclotides in planta is probably associated with host defense. Our current understanding of the cellular compartmentalization of cyclotides in the vacuole is based on indirect studies on transgenic model plants that do not express cyclotides naturally. Matrix-assisted laser desorption ionization (MALDI) imaging has also been used to study the distribution of cyclotides, but the technique's resolution was insufficient to determine their tissue or cell distribution. To avoid the limitations of these approaches, immunohistochemical visualization methods were used. Antibodies were raised in rabbits using cycloviolacin O2 (cyO2), and their specificity was determined by Western and dot blot experiments. Slides for immunohistochemical analysis were prepared from leaf, petiole and root fragments of Viola odorata and Viola uliginosa, and specimens were visualized using indirect epifluorescence microscopy. The antibodies against cyclotides were specific against selected bracelet cyclotides with high similarity (cyO2, cyO3, cyO8, cyO13) and suitable for immunohistochemistry. The tissue distribution of the cyclotides visualized in this way is consistent with their proposed role in host defense-relatively large quantities were observed in the leaf and petiole epidermis in both Viola species. Cyclotides were also found in vascular tissue in all the assessed plant organs. The vacuole storage of cyclotides was directly shown.

Viola woosanensis, a recurrent spontaneous hybrid between V. ulleungdoensis and V. chaerophylloides (Violaceae) endemic to Ulleung Island, Korea

Ulleung Island is an oceanic volcanic island in Korea, which has never been connected to the adjacent continent. Previous studies highlighted Ulleung Island as an excellent system to study the pattern and process of early stages of flowering plant evolutions on oceanic island. The predominant mode of speciation in flowering plants on Ulleung Island appears to be anagenesis. However, the potentially important role of hybrid speciation among incompletely reproductively isolated lineages cannot be ruled out. Viola woosanensis (Violaceae) is of purportedly hybrid origin between V. ulleungdoensis (i.e., formerly recognized as V. selkirkii in Ulleung Island) and V. chaerophylloides, based on morphology. To examine the origin of V. woosanensis, we sampled a total of 80 accessions, including V. woosanensis and its putative parental species and sequenced nrDNA ITS, and four highly variable chloroplast noncoding regions (trnL-trnF, rpl16 intron, atpF-atpH, and psbA-trnH). Representative species of Viola from Korea were also included in the phylogenetic analyses (maximum parsimony, maximum likelihood, and Bayesian inference). Additive polymorphic sites in the nrDNA ITS regions were confirmed by cloning amplicons from representative species. The molecular data strongly supported the hybrid origin of V. woosanensis, and the maternal and paternal parent were determined to be V. ulleungdoensis and V. chaerophylloides, respectively. The presence of two parental ribotypes in V. woosanensis (with the exception in one population) was confirmed by cloning, suggesting V. woosanensis is primarily the F1 generation. No trace of backcrossing and introgression with its parents was detected due to low fertility of hybrid species. We found a multiple and unidirectional hybrid origin of V. woosanensis. Additional studies are required to determine which factors contribute to asymmetric gene flow of Viola species in Ulleung Island.