Molecular phylogenetic analysis of Violaceae (Malpighiales) based on plastid and nuclear DNA sequences

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.

Galactosylated iron oxide nanoparticles for enhancing oral bioavailability of ceftriaxone

The current study focuses on the development, characterization, biocompatibility investigation and oral bioavailability evaluation of ceftriaxone (CFT)-loaded lactobionic acid (LBA)-functionalized iron oxide magnetic nanoparticles (MNP-LBA). Atomic force microscopy and dynamic light scattering showed that the developed CFT-loaded MNP-LBA is spherical, with a measured hydrodynamic size of 147 ± 15.9 nm and negative zeta potential values (-35 ± 0.58 mV). Fourier transformed infrared analysis revealed interactions between the nanocarrier and the drug. Nanoparticles showed high drug entrapment efficiencies of 91.5 ± 2.2%, and the drug was released gradually in vitro and shows prolonged in vitro stability using simulated gastrointestinal (GI) fluids. The formulations were found to be highly biocompatible (up to 100 µg/mL) and hemocompatible (up to 1.0 mg/mL). Using an albino rabbit model, the formulation showed a significant enhancement in drug plasma concentration up to 14.46 ± 2.5 µg/mL in comparison with its control (1.96 ± 0.58 µg/mL). Overall, the developed MNP-LBA formulation was found promising for provision of high-drug entrapment, gradual drug release and was appropriate for enhancing the oral delivery of CFT.

Cyclotides: From Structure to Function

This Review explores the class of plant-derived macrocyclic peptides called cyclotides. We include an account of their discovery, characterization, and distribution in the plant kingdom as well as a detailed analysis of their sequences and structures, biosynthesis and chemical synthesis, biological functions, and applications. These macrocyclic peptides are around 30 amino acids in size and are characterized by their head-to-tail cyclic backbone and cystine knot motif, which render them to be exceptionally stable, with resistance to thermal or enzymatic degradation. Routes to their chemical synthesis have been developed over the past two decades, and this capability has facilitated a wide range of mutagenesis and structure-activity relationship studies. In turn, these studies have both led to an increased understanding of their mechanisms of action as well as facilitated a range of applications in agriculture and medicine, as ecofriendly crop protection agents, and as drug leads or scaffolds for pharmaceutical design. Our overall objective in this Review is to provide readers with a comprehensive overview of cyclotides that we hope will stimulate further work on this fascinating family of peptides.

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.

Distribution of circular proteins in plants: large-scale mapping of cyclotides in the Violaceae

During the last decade there has been increasing interest in small circular proteins found in plants of the violet family (Violaceae). These so-called cyclotides consist of a circular chain of approximately 30 amino acids, including six cysteines forming three disulfide bonds, arranged in a cyclic cystine knot (CCK) motif. In this study we map the occurrence and distribution of cyclotides throughout the Violaceae. Plant material was obtained from herbarium sheets containing samples up to 200 years of age. Even the oldest specimens contained cyclotides in the preserved leaves, with no degradation products observable, confirming their place as one of the most stable proteins in nature. Over 200 samples covering 17 of the 23-31 genera in Violaceae were analyzed, and cyclotides were positively identified in 150 species. Each species contained a unique set of between one and 25 cyclotides, with many exclusive to individual plant species. We estimate the number of different cyclotides in the Violaceae to be 5000-25,000, and propose that cyclotides are ubiquitous among all Violaceae species. Twelve new cyclotides from six phylogenetically dispersed genera were sequenced. Furthermore, the first glycosylated derivatives of cyclotides were identified and characterized, further increasing the diversity and complexity of this unique protein family.

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

Allopolyploidization accounts for a significant fraction of speciation events in many eukaryotic lineages. However, existing phylogenetic and dating methods require tree-like topologies and are unable to handle the network-like phylogenetic relationships of lineages containing allopolyploids. No explicit framework has so far been established for evaluating competing network topologies, and few attempts have been made to date phylogenetic networks. We used a four-step approach to generate a dated polyploid species network for the cosmopolitan angiosperm genus Viola L. (Violaceae Batch.). The genus contains ca 600 species and both recent (neo-) and more ancient (meso-) polyploid lineages distributed over 16 sections. First, we obtained DNA sequences of three low-copy nuclear genes and one chloroplast region, from 42 species representing all 16 sections. Second, we obtained fossil-calibrated chronograms for each nuclear gene marker. Third, we determined the most parsimonious multilabeled genome tree and its corresponding network, resolved at the section (not the species) level. Reconstructing the "correct" network for a set of polyploids depends on recovering all homoeologs, i.e., all subgenomes, in these polyploids. Assuming the presence of Viola subgenome lineages that were not detected by the nuclear gene phylogenies ("ghost subgenome lineages") significantly reduced the number of inferred polyploidization events. We identified the most parsimonious network topology from a set of five competing scenarios differing in the interpretation of homoeolog extinctions and lineage sorting, based on (i) fewest possible ghost subgenome lineages, (ii) fewest possible polyploidization events, and (iii) least possible deviation from expected ploidy as inferred from available chromosome counts of the involved polyploid taxa. Finally, we estimated the homoploid and polyploid speciation times of the most parsimonious network. Homoploid speciation times were estimated by coalescent analysis of gene tree node ages. Polyploid speciation times were estimated by comparing branch lengths and speciation rates of lineages with and without ploidy shifts. Our analyses recognize Viola as an old genus (crown age 31 Ma) whose evolutionary history has been profoundly affected by allopolyploidy. Between 16 and 21 allopolyploidizations are necessary to explain the diversification of the 16 major lineages (sections) of Viola, suggesting that allopolyploidy has accounted for a high percentage-between 67% and 88%-of the speciation events at this level. The theoretical and methodological approaches presented here for (i) constructing networks and (ii) dating speciation events within a network, have general applicability for phylogenetic studies of groups where allopolyploidization has occurred. They make explicit use of a hitherto underexplored source of ploidy information from chromosome counts to help resolve phylogenetic cases where incomplete sequence data hampers network inference. Importantly, the coalescent-based method used herein circumvents the assumption of tree-like evolution required by most techniques for dating speciation events.

Molecular phylogenetic analysis of Passifloraceae sensu lato (Malpighiales) based on plastid and nuclear DNA sequences

A phylogenetic analysis of Passifloraceae sensu lato was performed using rbcL, atpB, matK, and 18S rDNA sequences from 25 genera and 42 species. Parsimony analyses of combined data sets resulted in a single most parsimonious tree, which was very similar to the 50% majority consensus tree from the Bayesian analysis. All nodes except three were supported by more than 50% bootstrap. The monophyly of Passifloraceae s.l. as well as the former families, Malesherbiaceae, Passifloraceae sensu stricto, and Turneraceae were strongly supported. Passifloraceae s.s. and the Turneraceae are sisters, and form a strongly supported clade. Within Passifloraceae s.s., the tribes Passifloreae and Paropsieae are both monophyletic. The intergeneric relationships within Passifloraceae s.s. and Turneraceae are roughly correlated with previous classification systems. The morphological character of an androgynophore/gynophore is better used for characterizing genera grouping within Passifloraceae s.s. Other morphological characters such as the corona and aril are discussed.

Inferring species networks from gene trees in high-polyploid North American and Hawaiian violets (Viola, Violaceae)

The phylogenies of allopolyploids take the shape of networks and cannot be adequately represented as bifurcating trees. Especially for high polyploids (i.e., organisms with more than six sets of nuclear chromosomes), the signatures of gene homoeolog loss, deep coalescence, and polyploidy may become confounded, with the result that gene trees may be congruent with more than one species network. Herein, we obtained the most parsimonious species network by objective comparison of competing scenarios involving polyploidization and homoeolog loss in a high-polyploid lineage of violets (Viola, Violaceae) mostly or entirely restricted to North America, Central America, or Hawaii. We amplified homoeologs of the low-copy nuclear gene, glucose-6-phosphate isomerase (GPI), by single-molecule polymerase chain reaction (PCR) and the chloroplast trnL-F region by conventional PCR for 51 species and subspecies. Topological incongruence among GPI homoeolog subclades, owing to deep coalescence and two instances of putative loss (or lack of detection) of homoeologs, were reconciled by applying the maximum tree topology for each subclade. The most parsimonious species network and the fossil-based calibration of the homoeolog tree favored monophyly of the high polyploids, which has resulted from allodecaploidization 9-14 Ma, involving sympatric ancestors from the extant Viola sections Chamaemelanium (diploid), Plagiostigma (paleotetraploid), and Viola (paleotetraploid). Although two of the high-polyploid lineages (Boreali-Americanae, Pedatae) remained decaploid, recurrent polyploidization with tetraploids of section Plagiostigma within the last 5 Ma has resulted in two 14-ploid lineages (Mexicanae, Nosphinium) and one 18-ploid lineage (Langsdorffianae). This implies a more complex phylogenetic and biogeographic origin of the Hawaiian violets (Nosphinium) than that previously inferred from rDNA data and illustrates the necessity of considering polyploidy in phylogenetic and biogeographic reconstruction.

DNA barcode sequence identification incorporating taxonomic hierarchy and within taxon variability

For DNA barcoding to succeed as a scientific endeavor an accurate and expeditious query sequence identification method is needed. Although a global multiple-sequence alignment can be generated for some barcoding markers (e.g. COI, rbcL), not all barcoding markers are as structurally conserved (e.g. matK). Thus, algorithms that depend on global multiple-sequence alignments are not universally applicable. Some sequence identification methods that use local pairwise alignments (e.g. BLAST) are unable to accurately differentiate between highly similar sequences and are not designed to cope with hierarchic phylogenetic relationships or within taxon variability. Here, I present a novel alignment-free sequence identification algorithm--BRONX--that accounts for observed within taxon variability and hierarchic relationships among taxa. BRONX identifies short variable segments and corresponding invariant flanking regions in reference sequences. These flanking regions are used to score variable regions in the query sequence without the production of a global multiple-sequence alignment. By incorporating observed within taxon variability into the scoring procedure, misidentifications arising from shared alleles/haplotypes are minimized. An explicit treatment of more inclusive terminals allows for separate identifications to be made for each taxonomic level and/or for user-defined terminals. BRONX performs better than all other methods when there is imperfect overlap between query and reference sequences (e.g. mini-barcode queries against a full-length barcode database). BRONX consistently produced better identifications at the genus-level for all query types.

Increased genetic diversity of Viola tricolor L. (Violaceae) in metal-polluted environments

Changes in DNA sequences affecting cryptic intraspecific variability are very important mechanisms of plant microevolutionary processes, initiating species diversification. In polluted environments, intra- and interpopulation changes at the molecular level proceed rapidly and lead to the formation of new ecotypes in a relatively short time. We used ISSR PCR fingerprinting data to analyze the genetic diversity and genetic structure of seven populations of Viola tricolor: four growing on soil contaminated with heavy metals (Zn, Pb, Cd; waste heaps) and three from control soil. The populations from the polluted sites showed higher genetic polymorphism (%(poly)=84%) and gene diversity (H(T)=0.1709) than the control populations (%(poly)=75% and H(T)=0.1448). The number of private markers we detected within metallicolous (MET) populations was more than double that found within non-metallicolous (NON) populations (15 vs. 7). The STRUCTURE and UPGMA analyses showed clear genetic differences between the NON and MET populations. Based on broad analyses of the genetic parameters, we conclude that the effect of these polluted environments on the genetic diversity of the MET populations, separating them from the NON populations, is evidence of microevolutionary processes at species level, leading to species divergence and the emergence of local ecotypes better adapted to their different environments.

Hawaiian angiosperm radiations of North American origin

BACKGROUND

Putative phytogeographical links between America (especially North America) and the Hawaiian Islands have figured prominently in disagreement and debate about the origin of Pacific floras and the efficacy of long-distance (oversea) plant dispersal, given the obstacles to explaining such major disjunctions by vicariance.

SCOPE

Review of past efforts, and of progress over the last 20 years, toward understanding relationships of Hawaiian angiosperms allows for a historically informed re-evaluation of the American (New World) contribution to Hawaiian diversity and evolutionary activity of American lineages in an insular setting.

CONCLUSIONS

Temperate and boreal North America is a much more important source of Hawaiian flora than suggested by most 20th century authorities on Pacific plant life, such as Fosberg and Skottsberg. Early views of evolution as too slow to account for divergence of highly distinctive endemics within the Hawaiian geological time frame evidently impeded biogeographical understanding, as did lack of appreciation for the importance of rare, often biotically mediated dispersal events and ecological opportunity in island ecosystems. Molecular phylogenetic evidence for North American ancestry of Hawaiian plant radiations, such as the silversword alliance, mints, sanicles, violets, schiedeas and spurges, underlines the potential of long-distance dispersal to shape floras, in accordance with hypotheses championed by Carlquist. Characteristics important to colonization of the islands, such as dispersibility by birds and ancestral hybridization or polyploidy, and ecological opportunities associated with 'sky islands' of temperate or boreal climate in the tropical Hawaiian archipelago may have been key to extensive diversification of endemic lineages of North American origin that are among the most species-rich clades of Hawaiian plants. Evident youth of flowering-plant lineages from North America is highly consistent with recent geological evidence for lack of high-elevation settings in the Hawaiian chain immediately prior to formation of the oldest, modern high-elevation island, Kaua'i.

Evolution of plant RNA polymerase IV/V genes: evidence of subneofunctionalization of duplicated NRPD2/NRPE2-like paralogs in Viola (Violaceae)

Background

DNA-dependent RNA polymerase IV and V (Pol IV and V) are multi-subunit enzymes occurring in plants. The origin of Pol V, specific to angiosperms, from Pol IV, which is present in all land plants, is linked to the duplication of the gene encoding the largest subunit and the subsequent subneofunctionalization of the two paralogs (NRPD1 and NRPE1). Additional duplication of the second-largest subunit, NRPD2/NRPE2, has happened independently in at least some eudicot lineages, but its paralogs are often subject to concerted evolution and gene death and little is known about their evolution nor their affinity with Pol IV and Pol V.

Results

We sequenced a ~1500 bp NRPD2/E2-like fragment from 18 Viola species, mostly paleopolyploids, and 6 non-Viola Violaceae species. Incongruence between the NRPD2/E2-like gene phylogeny and species phylogeny indicates a first duplication of NRPD2 relatively basally in Violaceae, with subsequent sorting of paralogs in the descendants, followed by a second duplication in the common ancestor of Viola and Allexis. In Viola, the mutation pattern suggested (sub-) neofunctionalization of the two NRPD2/E2-like paralogs, NRPD2/E2-a and NRPD2/E2-b. The d_N/d_S ratios indicated that a 54 bp region exerted strong positive selection for both paralogs immediately following duplication. This 54 bp region encodes a domain that is involved in the binding of the Nrpd2 subunit with other Pol IV/V subunits, and may be important for correct recognition of subunits specific to Pol IV and Pol V. Across all Viola taxa 73 NRPD2/E2-like sequences were obtained, of which 23 (32%) were putative pseudogenes - all occurring in polyploids. The NRPD2 duplication was conserved in all lineages except the diploid MELVIO clade, in which NRPD2/E2-b was lost, and its allopolyploid derivates from hybridization with the CHAM clade, section Viola and section Melanium, in which NRPD2/E2-a occurred in multiple copies while NRPD2/E2-b paralogs were either absent or pseudogenized.

Conclusions

Following the relatively recent split of Pol IV and Pol V, our data indicate that these two multi-subunit enzymes are still in the process of specialization and each acquiring fully subfunctionalized copies of their subunit genes. Even after specialization, the NRPD2/E2-like paralogs are prone to pseudogenization and gene conversion and NRPD2 and NRPE2 copy number is a highly dynamic process modulated by allopolyploidy and gene death.

Malpighiales phylogenetics: Gaining ground on one of the most recalcitrant clades in the angiosperm tree of life

The eudicot order Malpighiales contains ∼16000 species and is the most poorly resolved large rosid clade. To clarify phylogenetic relationships in the order, we used maximum likelihood, Bayesian, and parsimony analyses of DNA sequence data from 13 gene regions, totaling 15604 bp, and representing all three genomic compartments (i.e., plastid: atpB, matK, ndhF, and rbcL; mitochondrial: ccmB, cob, matR, nad1B-C, nad6, and rps3; and nuclear: 18S rDNA, PHYC, and newly developed low-copy EMB2765). Our sampling of 190 taxa includes representatives from all families of Malpighiales. These data provide greatly increased support for the recent additions of Aneulophus, Bhesa, Centroplacus, Ploiarium, and Rafflesiaceae to Malpighiales; sister relations of Phyllanthaceae + Picrodendraceae, monophyly of Hypericaceae, and polyphyly of Clusiaceae. Oxalidales + Huaceae, followed by Celastrales are successive sisters to Malpighiales. Parasitic Rafflesiaceae, which produce the world's largest flowers, are confirmed as embedded within a paraphyletic Euphorbiaceae. Novel findings show a well-supported placement of Ctenolophonaceae with Erythroxylaceae + Rhizophoraceae, sister-group relationships of Bhesa + Centroplacus, and the exclusion of Medusandra from Malpighiales. New taxonomic circumscriptions include the addition of Bhesa to Centroplacaceae, Medusandra to Peridiscaceae (Saxifragales), Calophyllaceae applied to Clusiaceae subfamily Kielmeyeroideae, Peraceae applied to Euphorbiaceae subfamily Peroideae, and Huaceae included in Oxalidales.

Circling the enemy: cyclic proteins in plant defence

Cyclotides are ultra-stable plant proteins that have a circular peptide backbone crosslinked by a cystine knot of disulfide bonds. They are produced in large quantities by plants of the Violaceae and Rubiaceae families and have a role in plant defence against insect predation. As I discuss here, recent studies have begun to reveal how their unique circular topology evolved. Cyclization is achieved by hijacking existing plant proteolytic enzymes and operating them in 'reverse' to form a peptide bond between the N- and C-termini of a linear precursor. Such studies suggest that circular proteins are more common in the plant kingdom than was previously thought, and their exceptional stability has led to their application as protein-engineering templates in drug design.