Comparative Analysis of Shikonin and Alkannin Acyltransferases Reveals Their Functional Conservation in Boraginaceae

Shikonin and its enantiomer, alkannin, are bioactive naphthoquinones produced in several plants of the family Boraginaceae. The structures of these acylated derivatives, which have various short-chain acyl moieties, differ among plant species. The acylation of shikonin and alkannin in Lithospermum erythrorhizon was previously reported to be catalyzed by two enantioselective BAHD acyltransferases, shikonin O-acyltransferase (LeSAT1) and alkannin O-acyltransferase (LeAAT1). However, the mechanisms by which various shikonin and alkannin derivatives are produced in Boraginaceae plants remain to be determined. In the present study, evaluation of six Boraginaceae plants identified 23 homologs of LeSAT1 and LeAAT1, with 15 of these enzymes found to catalyze the acylation of shikonin or alkannin, utilizing acetyl-CoA, isobutyryl-CoA or isovaleryl-CoA as an acyl donor. Analyses of substrate specificities of these enzymes for both acyl donors and acyl acceptors and determination of their subcellular localization using Nicotiana benthamiana revealed a distinct functional differentiation of BAHD acyltransferases in Boraginaceae plants. Gene expression of these acyltransferases correlated with the enantiomeric ratio of produced shikonin/alkannin derivatives in L. erythrorhizon and Echium plantagineum. These enzymes showed conserved substrate specificities for acyl donors among plant species, indicating that the diversity in acyl moieties of shikonin/alkannin derivatives involved factors other than the differentiation of acyltransferases. These findings provide insight into the chemical diversification and evolutionary processes of shikonin/alkannin derivatives.

Comparative Analyses of Complete Chloroplast Genomes of Microula sikkimensis and Related Species of Boraginaceae

The present study provides a detailed analysis of the chloroplast genome of Microula sikkimensis. The genome consisted of a total of 149,428 bp and four distinct regions, including a large single-copy region (81,329 bp), a small single-copy region (17,261 bp), and an inverted repeat region (25,419 bp). The genome contained 112 genes, including 78 protein-coding genes, 30 tRNA genes, and 4 rRNA genes, and some exhibited duplication in the inverted repeat region. The chloroplast genome displayed different GC content across regions, with the inverted repeat region exhibiting the highest. Codon usage analysis and the identification of simple sequence repeats (SSRs) offer valuable genetic markers. Comparative analysis with other Boraginaceae species highlighted conservation and diversity in coding and noncoding regions. Phylogenetic analysis placed M. sikkimensis within the Boraginaceae family, revealing its distinct relationship with specific species.

[Effect of acetylalkannin from Arnebia euchroma on proliferation, migration, and invasion of human melanoma A375 cells]

This study aimed to explore the effect and mechanism of acetylalkannin from Arnebia euchroma on the proliferation, migration, and invasion of human melanoma A375 cells. A375 cells were divided into a blank group, and low-, medium-, and high-dose acetylalkannin groups(0.5, 1.0, and 2.0 μmol·L~(-1)). The MTT assay was used to detect cell proliferation. Cell scratch and transwell migration assays were used to detect cell migration ability, and the transwell invasion assay was used to detect cell invasion ability. Western blot was used to detect the protein expression of migration and invasion-related N-cadherin, vimentin, matrix metalloproteina-se-9(MMP-9), and Wnt/β-catenin pathway-related Wnt1, Axin2, glycogen synthase kinase-3β(GSK-3β), phosphorylated GSK-3β(p-GSK-3β), β-catenin, cell cycle protein D_1(cyclin D_1), and p21. Real-time fluorescence-based quantitative polymerase chain reaction(real-time PCR) was used to detect the mRNA expression of E-cadherin, matrix metalloproteinase-2(MMP-2), N-cadherin, vimentin, β-catenin, snail-1, and CD44. MTT results showed that the cell inhibition rates in the acetylalkannin groups significantly increased as compared with that in the blank group(P<0.01). The results of cell scratch and transwell assays showed that compared with the blank group, the acetylalkannin groups showed reduced cell migration and invasion, and migration and invasion rates(P<0.05, P<0.01) and weakened horizontal and vertical migration and invasion abilities. Western blot results showed that compared with the blank group, the high-dose acetylalkannin group showed increased expression of Axin2 protein(P<0.05), and decreased expression of N-cadherin, vimentin, MMP-9, Wnt1, p-GSK-3β, β-catenin, cyclin D_1, and p21 proteins(P<0.05, P<0.01). The expression of GSK-3β protein did not change significantly. PCR results showed that the overall trend of MMP-2, N-cadherin, vimentin, β-catenin, snail-1, and CD44 mRNA expression was down-regulated(P<0.01), and the expression of E-cadherin mRNA increased(P<0.01). Acetylalkannin can inhibit the proliferation, migration, and invasion of human melanoma A375 cells, and its mechanism of action may be related to the regulation of Wnt/β-catenin signaling pathway.

Characterization of Arnebia euchroma PGT homologs involved in the biosynthesis of shikonin

Shikonin is a red naphthoquinone natural product from plants with high economical and medical values. The para-hydroxybenzoic acid geranyltransferase (PGT) catalyzes the key regulatory step of shikonin biosynthesis. PGTs from Lithospermum erythrorhizon have been well-characterized and used in industrial shikonin production. However, its perennial medicinal plant Arnebia euchroma accumulates much more pigment and the underlying mechanism remains obscure. Here, we discovered and characterized the different isoforms of AePGTs. Phylogenetic study and structure modeling suggested that the N-terminal of AePGT6 contributed to its highest activity among 7 AePGTs. Indeed, AePGT2 and AePGT3 fused with 60 amino acids from the N-terminal of AePGT6 showed even higher activity than AePGT6, while native AePGT2 and AePGT3 don't have catalytic activity. Our result not only provided a mechanistic explanation of high shikonin contents in Arnebia euchroma but also engineered a best-performing PGT to achieve the highest-to-date production of 3-geranyl-4-hydroxybenzoate acid, an intermedium of shikonin.

De novo transcriptome analysis provides insights into formation of in vitro adventitious root from leaf explants of Arnebia euchroma

BACKGROUND

Adventitious root formation is considered a major developmental step during the propagation of difficult to root plants, especially in horticultural crops. Recently, adventitious roots induced through plant tissue culture methods have also been used for production of phytochemicals such as flavonoids, anthocyanins and anthraquinones. It is rather well understood which horticultural species will easily form adventitious roots, but the factors affecting this process at molecular level or regulating the induction process in in vitro conditions are far less known. The present study was conducted to identify transcripts involved in in vitro induction and formation of adventitious roots using Arnebia euchroma leaves at different time points (intact leaf (control), 3 h, 12 h, 24 h, 3 d, 7 d, 10 d and 15 d). A. euchroma is an endangered medicinal Himalayan herb whose root contains red naphthoquinone pigments. These phytoconstituents are widely used as an herbal ingredient in Asian traditional medicine as well as natural colouring agent in food and cosmetics.

RESULTS

A total of 137.93 to 293.76 million raw reads were generated and assembled to 54,587 transcripts with average length of 1512.27 bps and N50 of 2193 bps, respectively. In addition, 50,107 differentially expressed genes were identified and found to be involved in plant hormone signal transduction, cell wall modification and wound induced mitogen activated protein kinase signalling. The data exhibited dominance of auxin responsive (AUXIN RESPONSE FACTOR8, IAA13, GRETCHEN HAGEN3.1) and sucrose translocation (BETA-31 FRUCTOFURANOSIDASE and MONOSACCHARIDE-SENSING protein1) genes during induction phase. In the initiation phase, the expression of LATERAL ORGAN BOUNDARIES DOMAIN16, EXPANSIN-B15, ENDOGLUCANASE25 and LEUCINE-rich repeat EXTENSION-like proteins was increased. During the expression phase, the same transcripts, with exception of LATERAL ORGAN BOUNDARIES DOMAIN16 were identified. Overall, the transcriptomic analysis revealed a similar patterns of genes, however, their expression level varied in subsequent phases of in vitro adventitious root formation in A. euchroma.

CONCLUSION

The results presented here will be helpful in understanding key regulators of in vitro adventitious root development in Arnebia species, which may be deployed in the future for phytochemical production at a commercial scale.

CRISPR/Cas9-Mediated Genome Editing in Comfrey (Symphytum officinale) Hairy Roots Results in the Complete Eradication of Pyrrolizidine Alkaloids

Comfrey (Symphytum officinale) is a medicinal plant with anti-inflammatory, analgesic, and proliferative properties. However, its pharmaceutical application is hampered by the co-occurrence of toxic pyrrolizidine alkaloids (PAs) in its tissues. Using a CRISPR/Cas9-based approach, we introduced detrimental mutations into the hss gene encoding homospermidine synthase (HSS), the first pathway-specific enzyme of PA biosynthesis. The resulting hairy root (HR) lines were analyzed for the type of gene-editing effect that they exhibited and for their homospermidine and PA content. Inactivation of only one of the two hss alleles resulted in HRs with significantly reduced levels of homospermidine and PAs, whereas no alkaloids were detectable in HRs with two inactivated hss alleles. PAs were detectable once again after the HSS-deficient HRs were fed homospermidine confirming that the inability of these roots to produce PAs was only attributable to the inactivated HSS and not to any unidentified off-target effect of the CRISPR/Cas9 approach. Further analyses showed that PA-free HRs possessed, at least in traces, detectable amounts of homospermidine, and that the PA patterns of manipulated HRs were different from those of control lines. These observations are discussed with regard to the potential use of such a CRISPR/Cas9-mediated approach for the economical exploitation of in vitro systems in a medicinal plant and for further studies of PA biosynthesis in non-model plants.

[Cloning and functional analysis of caffeic acid and rosmarinic acid glycosyltransferases from Arnebia euchroma]

Caffeic acid and its oligomers are the main water-soluble active constituents of the traditional Chinese medicine(TCM) Arnebiae Radix. These compounds possess multiple biological activities such as antimicrobial, antioxidant, cardiovascular protective, liver protective, anti-liver fibrosis, antiviral and anticancer activities. The phenylpropanoid pathway in plants is responsible for the biosynthesis of caffeic acid and its oligomers. Glycosylation can change phenylpropanoid solubility, stability and toxic potential, as well as influencing compartmentalization and biological activity. In view of the important role played by de-glycosylation in the regulation of phenylpropanoid homeostasis, the biosynthesis of caffeic acid and its oligomers are supposed to be under the control of relative UDP-glycosyltransferases(UGTs). Through the data mining of Arnebia euchroma transcriptome, we cloned 15 full-length putative UGT genes. After recombinant expression using the prokaryotic system, the crude enzyme solution of the putative UGTs was examined for the glycosylation activities towards caffeic acid and rosmarinic acid in vitro. AeUGT_01, AeUGT_02, AeUGT_03, AeUGT_04 and AeUGT_10 were able to glycosylate caffeic acid and/or rosmarinic acid resulting in different mono-and/or di-glycosylated products in the UPLC-MS analyses. The characterized UGTs were distantly related to each other and divided into different clades of the phylogenetic tree. Based on the observation that each characterized UGT exhibited substrate or catalytic similarity with the members in their own clade, we supposed the glycosylation abilities towards caffeic acid and/or rosmarinic acid were evolved independently in different clades. The identification of caffeic acid and rosmarinic acid UGTs from A. euchroma could lead to deeper understanding of the caffeic acid oligomers biosynthesis and its regulation. Furthermore, these UGTs might be used for regiospecific glycosylation of caffeic acid and rosmarinic acid to produce bioactive compounds for potential therapeutic applications.

Genetic diversity, reproductive capacity and alkaloids content in three endemic Alkanna species

The Balkans endemic species Alkanna primuliflora Griseb., A. stribrnyi Velen., and A. graeca Boiss. & Spruner have limited distribution in the Balkan Peninsula and a large variation in the morphological characteristics. The populations of the three Alkanna species in the Bulgarian flora are small and fragmented. There are no previous reports on the chemical profile or on the embryology of these species. The hypothesis was that the limited distribution of A. primuliflora, A. stribrnyi, and A. graeca was due to their reproductive capacity and genetic diversity. Furthermore, we hypothesized that the three species will contain pyrrolizidine alkaloids (PAs), as other species of the genus Alkanna (Boraginaceae), but they would have differential alkaloids composition. The population genetic structure and differentiation showed a clear distinction between species and revealed average levels of genetic diversity among the natural populations of the three Alkanna species. The embryological investigation observed stability of the processes in the male and female generative spheres and high viability of mature pollen and embryo of the three species. The normal formation of male and female gametophytes without deviations or degenerative processes, and observed levels of genetic diversity between Alkanna individuals are important in maintaining the size and resilience of the Alkanna populations. Eight alkaloids were identified by GC-MS in A. primuliflora and A. graeca and six alkaloids in A. stribrnyi. The main pyrrolizidine alkaloids (PAs) in all investigated species was triangularine. A. primuliflora and A. graeca showed similar chemical composition that comprised 9-angeloylretronecine, 7-tigloylretronecine, 9-tigloylretronecine, triangularicine, dihydroxytriangularine, dihydroxytriangularicine, whereas, in A. stribrnyi 9-tigloylretronecine, triangularicine and dihydroxytriangularicine were not found. This is the first report on the presence of PAs in A. primuliflora, A. stribrnyi and A. graeca. Besides, this is the first report on the embryology of these endemic species. The results contribute to the knowledge of the three endemic Alkanna species and will facilitate policy-making and defining new strategies for their conservation.

CYP76B74 Catalyzes the 3''-Hydroxylation of Geranylhydroquinone in Shikonin Biosynthesis

Shikonin and its derivatives are the most abundant naphthoquinone pigments formed in species of the medicinally and economically valuable Boraginaceae. A key step in the shikonin biosynthetic pathway, namely the C-3'' hydroxylation of the prenylated phenolic intermediate geranylhydroquinone, is expected to be catalyzed by a cytochrome P450. To identify cytochrome P450 candidates with transcription profiles similar to those of genes known to be involved in shikonin biosynthesis, we carried out coexpression analysis of transcriptome data sets of shikonin-proficient and shikonin-deficient cell lines and examined the spatial expression of candidate genes in different organs of Arnebia euchroma In biochemical assays using geranylhydroquinone as the substrate, CYP76B74 exhibited geranylhydroquinone 3''-hydroxylase activity and produced 3''-hydroxy-geranylhydroquinone. In CYP76B74 RNA interference A. euchroma hairy roots, shikonin derivative accumulation decreased dramatically, which demonstrated that CYP76B74 is required for shikonin biosynthesis in the plant. Phylogenetic analysis confirmed that CYP76B74 belonged to the CYP76B subfamily and was most likely derived from an ancestral geraniol 10-hydroxylase. In a subcellular localization analysis, a GFP-CYP76B74 fusion localized to endoplasmic reticulum membranes. Our results demonstrate that CYP76B74 catalyzes the key hydroxylation step in shikonin biosynthesis with high efficiency. The characterization of the CYP76B74 described here paves the way for further exploration of the ring closure reactions generating the naphthoquinone skeleton as well as for the alternative metabolism of geranylhydroquinone 3''-hydroxylase to dihydroechinofuran.

Gender specialisation and stigma height dimorphism in Mediterranean Lithodora fruticosa (Boraginaceae)

Dimorphism in style height has evolved repeatedly in flowering plants, with some individuals having short and others long styles; in the case of distylous species, stigma position varies reciprocally with that of the anthers. Distyly can be associated with divergence in the functional gender between long- and short-styled individuals, but gender divergence has rarely been investigated in species with a simple stigma height polymorphism in the absence of reciprocal dimorphism in anther position. To evaluate the relation between stigma height polymorphism and gender, I measured the dimensions of floral morphology and seed production for the two morphs of a large population of the Iberian species Lithodora fruticosa (Boraginaceae). Results confirm the existence of a stigma height polymorphism in L. fruticosa, with long- and short-styled individuals at a 1:1 ratio in the studied population. Long-styled individuals produced substantially more seeds than did short-styled individuals, pointing to strong divergence in functional gender between the two morphs. The results of this study are puzzling in light of recent work that suggests that L. fruticosa has a multi-allelic self-incompatibility system. I discuss the significance of gender divergence in L. fruticosa and evaluate hypotheses that might explain it.

Memoirs of a frequent flier: Phylogenomics reveals 18 long-distance dispersals between North America and South America in the popcorn flowers (Amsinckiinae)

PREMISE OF THE STUDY

American amphitropical disjunction (AAD) is an important but understudied New World biogeographic pattern in which related plants occur in extratropical North America and South America, but are absent in the intervening tropics. Subtribe Amsinckiinae (Boraginaceae) is one of the richest groups of plants displaying the AAD pattern. Here, we infer a time-calibrated molecular phylogeny of the group to evaluate the number, timing, and directionality of AAD events, which yields generalizable insights into the mechanism of AAD.

METHODS

We perform a phylogenomic analysis of 139 samples of subtribe Amsinckiinae and infer divergence times using two calibration schemes: with only fossil calibrations and with fossils plus a secondary calibration from a recent family level analysis. Biogeographic analysis was performed in the R package BioGeoBEARS.

KEY RESULTS

We document 18 examples of AAD in the Amsinckiinae. Inferred divergence times of these AAD examples were strongly asynchronous, ranging from Miocene (17.1 million years ago [Ma]) to Pleistocene (0.33 Ma), with most (12) occurring <5 Ma. Four events occurred 10-5 Ma, during the second rise of the Andes. All AAD examples had a North America to South America directionality.

CONCLUSIONS

Second only to the hyperdiverse Poaceae in number of documented AAD examples, the Amsinckiinae is an ideal system for the study of AAD. Asynchronous divergence times support the hypothesis of long-distance dispersal by birds as the mechanism of AAD in the subtribe and more generally. Further comparative phylogenomic studies may permit biogeographic hypothesis testing and examination of the relationship between AAD and fruit morphology, reproductive biology, and ploidy.

Transcriptome analysis explores genes related to shikonin biosynthesis in Lithospermeae plants and provides insights into Boraginales' evolutionary history

Shikonin and its derivatives extracted from Lithospermeae plants' red roots have current applications in food and pharmaceutical industries. Previous studies have cloned some genes related to shikonin biosynthesis. However, most genes related to shikonin biosynthesis remain unclear, because the lack of the genome/transcriptome of the Lithospermeae plants. Therefore, in order to provide a new understanding of shikonin biosynthesis, we obtained transcriptome data and unigenes expression profiles in three shikonin-producing Lithospermeae plants, i.e., Lithospermum erythrorhizon, Arnebia euchroma and Echium plantagineum. As a result, two unigenes (i.e., G10H and 12OPR) that are involved in "shikonin downstream biosynthesis" and "methyl jasmonate biosynthesis" were deemed to relate to shikonin biosynthesis in this study. Furthermore, we conducted a Lamiids phylogenetic model and identified orthologous unigenes under positive selection in above three Lithospermeae plants. The results indicated Boraginales was more relative to Solanales/Gentianales than to Lamiales.

Chemotaxonomic perspectives of the Paracaryum (Cynoglosseae, Boraginaceae) taxa based on fruit fatty acid composition

Paracaryum is a medium sized genus in Cynoglosseae. This study represents the most comprehensive phytochemical investigation of Paracaryum to date. The fatty acid compositions of the fruits of ten Paracaryum taxa belonging to three different subgenera were investigated for chemotaxonomic allocation using gas chromatography. The fatty acid profiles of seven Paracaryum taxa, five of which are endemic to Turkey, were recorded for the first time. Among the twenty-two analysed fatty acids, oleic, linoleic and α-linolenic acids were the major fatty acids represented. The oleic acid content ranged from 22.1% in P. hirsutum to 51.7% in P. lithospermifolium subsp. cariense var. erectum; linolenic acid content ranged from 8.6% in P. lithospermifolium subsp. cariense var. erectum to 20.7% in P. erysimifolium; α-linolenic acid content ranged from 7.5% in P. lithospermifolium subsp. cariense var. erectum to 13.5% in P. cristatum subsp. cristatum; gamma linolenic acid content ranged from 2.8% in P. erysimifolium to 6.0% in P. hirsutum. Additional fatty acids also displayed varying levels in different species; palmitic acid content accounted for 17.7% in P. erysimifolium, erucic acid content was 8.73% in P. strictum, eicosenoic acid content was 6.0% in P. cristatum subsp. cristatum, eicosadienoic acid content was 4.4% in P. hirsutum, and stearic acid content was 4.3% in P. erysimifolium. The classification of the tribe Cynoglosseae remains controversial despite the many intensive morphological and phylogenetic investigations that have been carried out. Our fatty acid data from Paracaryum were analysed together with previously recorded fatty acid data from Cynoglosseae s.l. taxa to examine the chemotaxonomic contribution to the classification among taxa in Cynoglosseae by multivariate methods, including the unweighted pair group method with arithmetic mean and principal component analysis. An assessment of these chemometrics data supported the chemotaxonomic position of the genus Paracaryum in the tribe Cynoglosseae. While the principal component graphic did not depict clear separation of the three subgenera of Paracaryum, the principal component analysis revealed the chemotaxonomic significance of palmitic, linoleic, capric, and oleic acids.

[Transcriptome-based gene mining and bioinformatics analysis of p-hydroxybenzoate geranyltransferase genes in Arnebia euchroma]

The p-hydroxybenzoate geranyltransferases(PGT) play an important role in the biosynthesis pathways of shikonin derivatives. Six PGTs were obtained from transcriptome datebase of Arnebia euchroma by using bioinformatics methods and the proteins＇physiochemical properties they encoded were predicted. The result of protein domain prediction showed all of the six protein sequences contained the conserved domain of Ubia prenyltransferase family and possessed the motif NDxxDxxxD for prenyldiphosphate binding and a GX(K/Y)STAL sequence for putative aromatic ring binding. The phylogenetic tree showed that PGT and p-hydroxybenzoate polyprenyltransferase(PPT) belonged to two different clades. The results of gene expression analyses showed that the expression levels in the red shikonin-proficient line and the overground part of A. euchroma that could produce shikonin derivatives was much higher than the white shikonin-deficient line and the underground part, which suggested a positive correlation between the expression levels of PGT genes and shikonin production. This study aims to lay a foudation for further understanding of the function and enzymatic properties of PGT and provide a basis for the biosynthesis pathways and metabolic regulation of shikonin derivatives.

Morphological and genetic discrepancies in populations of Oreocarya paradoxa and O. revealii: the impact of edaphic selection on recent diversification in the Colorado Plateau

PREMISE OF THE STUDY

Investigations of recently derived and edaphically (soil) defined plant systems have provided insight into important mechanisms of ecological divergence. We investigated the impact of edaphic adaptation on recent divergence between two Colorado Plateau endemics: the gypsum facultative Oreocarya revealii (Boraginaceae) and its more generalist sister species O. paradoxa. We assessed morphological stability, genetic identity, and soil chemistry to determine whether O. revealii is a distinct lineage edaphically adapted from O. paradoxa, as has been described in the literature.

METHODS

We genotyped 21 populations throughout the ranges of both species using 11 microsatellite markers and three plastid regions (trnL-F, trnT-L, trnQ-rps16) for haplotype analysis. We compared these data with soil chemistry (Ca and S concentrations, indicating gypsum levels), location, and morphological identity of populations.

KEY RESULTS

Soil chemistry failed to explain genetic or morphological identity in either taxon. Haplotype analysis suggests ancestral variation in the more geographically restricted O. revealii, along with regional geographic isolation. A discontinuity was identified between the morphological and genetic identity in several populations, suggesting incomplete lineage sorting and the nonfixation of identifying morphological traits.

CONCLUSIONS

Oreocarya revealii is unlikely to have arisen via edaphic selection, because soil chemistry of population sites, morphology of individuals, and genetic identity are not strongly correlated. The nonfixation of identifying traits is likely a result of recent divergence in this system, and the potentiality of such discrepancies should be considered when investigating recently diversified gypsum-associated groups.

Single cell subtractive transcriptomics for identification of cell-specifically expressed candidate genes of pyrrolizidine alkaloid biosynthesis

Progress has recently been made in the elucidation of pathways of secondary metabolism. However, because of its diversity, genetic information concerning biosynthetic details is still missing for many natural products. This is also the case for the biosynthesis of pyrrolizidine alkaloids. To close this gap, we tested strategies using tissues that express this pathway in comparison to tissues in which this pathway is not expressed. As many pathways of secondary metabolism are known to be induced by jasmonates, the pyrrolizidine alkaloid-producing species Heliotropium indicum, Symphytum officinale, and Cynoglossum officinale of the Boraginales order were treated with methyl jasmonate. An effect on pyrrolizidine alkaloid levels and on transcript levels of homospermidine synthase, the first specific enzyme of pyrrolizidine alkaloid biosynthesis, was not detectable. Therefore, a method was developed by making use of the often observed cell-specific production of secondary compounds. H. indicum produces pyrrolizidine alkaloids exclusively in the shoot. Homospermidine synthase is expressed only in the cells of the lower leaf epidermis and the epidermis of the stem. Suggesting that the whole pathway of pyrrolizidine alkaloid biosynthesis might be localized in these cells, we have isolated single cells of the upper and lower epidermis by laser-capture microdissection. The resulting cDNA preparations have been used in a subtractive transcriptomic approach. Quantitative real-time polymerase chain reaction has shown that the resulting library is significantly enriched for homospermidine-synthase-coding transcripts providing a valuable source for the identification of further genes involved in pyrrolizidine alkaloid biosynthesis.

Cytomictic anomalous male meiosis and 2n pollen grain formation in Mertensia echioides Benth. (Boraginaceae) from Kashmir Himalaya

Presently Mertensia echioides Benth. (Boraginaceae) collected from Kashmir Himalaya, India, is cytologically analyzed for the first time revealing 2n = 2x = 24 (diploid). Interestingly we found 4.3-6.2% syncytic meiocytes/PMCs with 2n = 4x = 48 (tetraploid) in addition to normal meiocytes (2n = 24) during male meiosis. These comparatively larger PMCs (pollen mother cells) lead to the formation of fertile giant 2n pollen grains. A frequency of 6.4-13.3% PMCs shows transfer of chromatin material at prophase-I and, therefore, results in aneuploid meiocytes. Whole chromatin transfer by the process of cytomixis could also have led to the formation of tetraploid cells. Translocation heterozygosity is also evident in the form of multivalents in 12-17% diploid (2x) meiocytes at diakinesis and metaphase-I and is reported for the first time in this species. The syncytes formed depict open chain hexavalent and quadrivalent formation in the three populations with different frequencies. Moreover chromatin stickiness at metaphase-I is observed in 45% of PMCs in population-1 (P-1). Syncyte or unreduced PMC formation leading to unreduced fertile gametes is here speculated to act as a possible way out for infraspecific polyploidization in the species.

[Transcriptome-based bioinformatics analysis of Arnebia euchroma ERF transcription factor family]

Twenty-seven ERF transcription factor family genes were isolated from Arnebia euchroma, with an average size of 1,010 bp, each gene encoded a 212 amino acids on average. The gene structure and expression of physicochemical properties, subcellular localization, signal peptides, senior structural domains and conservative forecasting, and analysis of A. euchroma were studied comparing with ERF gene gi261363612 of Lithospermum erythrorhizon, and phylogenetic analysis of A. euchroma ERF family was carried out. The results showed the existence of three conserved domains in this family, the senior structure based on random coil and it clustered into CBF/DREB and ERF subfamilies.

CsIFIT1, an interferon-induced protein with tetratricopeptide repeat, inhibits viral infection in tongue sole (Cynoglossus semilaevis)

Interferon-induced protein with tetratricopeptide repeats (IFIT) is involved in antiviral immunity in mammalian cells. In fish, the function of this group of proteins is unknown. In this study, we report the identification of an IFIT homologue (named CsIFIT1) from tongue sole (Cynoglossus semilaevis) and examined its antiviral potential. CsIFIT1 is composed of 436 residues and shares 30.9%-58.7% overall sequence identities with the known IFITs of a number of teleost species. In silico analysis identified two tetratricopeptide repeats in CsIFIT1. Quantitative real time RT-PCR analysis showed that CsIFIT1 expression was most abundant in immune relevant organs and upregulated, in a time-dependent manner, by experimental infection with the fish viral pathogen megalocytivirus. When CsIFIT1 was overexpressed in tongue sole before megalocytivirus infection, viral replication in fish tissues was blocked to significant extents. Consistently, when CsIFIT1 expression in tongue sole was knocked down by siRNA, viral replication in fish tissues was significantly enhanced. Taken together, these results indicate that CsIFIT1 is likely a key factor of antiviral immunity and is required for optimal defense against viral infection.

[Overexpression and RNAi vectors built for key secondary metabolic pathway genes PAL, HMGR, PGT of Arnebia euchroma]

Arnebia euchroma is the main source for medicinal herb Zicao. and its most important component shikonin compounds have high medicinal and industrial value. This research is aimed to build overexpression vectors and RNAi vectors for key secondary metabolism genes of A. euchroma, and bulid platform for constructions of related transgenic lines using GATEWAY technology. To build genetic material based genetic research platform is to provide a great convenience for digging and functional verification of the genes on secondary metabolic pathway, and also to fill the gaps in transgenic research of A. euchroma. This study is also important for the cultivation of shikonin high-yielding strains of A. euchroma.

[In silico cloning and bioinformatics analysis of an AP2/EFR family gene from Arnebia euchroma]

A cDNA sequence of Arnebia euchroma AP2/ERF named AeAP2/ERF was cloned by in silico cloning in this study, using ACX71873 sequence from Lithospermum erythrorhizon as the probe sequence. Some characters of the AP2/ERF gene and encoded protein sequences were predicted and analyzed by the bioinformatics methods, including general physical and chemical properties, hydrophobieity, signal peptide, secondary structure, localization sites in cells. Results showed that the 876 bp long gene included a 1 077 bp ORF and encoding 205 amino acid. The AeAP2/ERF protein had no signal peptide, it was a hydrophilic proteins located in nucleus. The function of the AP2/ERF protein was mainly involved with metabolism controlling and signal transduction.

Identification and functional analysis of a Δ6-desaturase gene and the effects of temperature and wounding stresses on its expression in Microula sikkimensis leaves

A Δ6-desaturase gene was isolated from Microula sikkimensis. Sequence analysis indicated that the gene, designated MsD6DES, had an open reading frame of 1,357 bp and encoded 448 amino acids. Heterologous expression in tobacco indicated that MsD6DES can use endogenous substrates to synthesize γ-linolenic acid (GLA, 18:3(Δ 6,9,12)) and octadecatetraenoic acid (OTA, 18:4(Δ 6,9,12,15)). MsD6DES transcripts were distributed in all tested tissues, with high expression levels in seeds and young leaves. The effects of temperature and wounding stresses on MsD6DES expression were analyzed. The results indicated that temperature regulates MsD6DES at the transcriptional level. MsD6DES expression increased first, reaching a maximum 4 h after low-temperature treatment. A slight increase in MsD6DES transcript levels was also observed under high temperature. We found that the response of MsD6DES to temperature stress was different from those of fungi and algae. In addition, MsD6DES was found to be wound-inducible.

Distinct cell-specific expression of homospermidine synthase involved in pyrrolizidine alkaloid biosynthesis in three species of the boraginales

Homospermidine synthase (HSS) is the first specific enzyme in pyrrolizidine alkaloid (PA) biosynthesis, a pathway involved in the plant's chemical defense. HSS has been shown to be recruited repeatedly by duplication of a gene involved in primary metabolism. Within the lineage of the Boraginales, only one gene duplication event gave rise to HSS. Here, we demonstrate that the tissue-specific expression of HSS in three boraginaceous species, Heliotropium indicum, Symphytum officinale, and Cynoglossum officinale, is unique with respect to plant organ, tissue, and cell type. Within H. indicum, HSS is expressed exclusively in nonspecialized cells of the lower epidermis of young leaves and shoots. In S. officinale, HSS expression has been detected in the cells of the root endodermis and in leaves directly underneath developing inflorescences. In young roots of C. officinale, HSS is detected only in cells of the endodermis, but in a later developmental stage, additionally in the pericycle. The individual expression patterns are compared with those within the Senecioneae lineage (Asteraceae), where HSS expression is reproducibly found in specific cells of the endodermis and the adjacent cortex parenchyma of the roots. The individual expression patterns within the Boraginales species are discussed as being a requirement for the successful recruitment of HSS after gene duplication. The diversity of HSS expression within this lineage adds a further facet to the already diverse patterns of expression that have been observed for HSS in other PA-producing plant lineages, making this PA-specific enzyme one of the most diverse expressed proteins described in the literature.

Seasonal change in a pollinator community and the maintenance of style length variation in Mertensia fusiformis (Boraginaceae)

BACKGROUND AND AIMS

In sub-alpine habitats, patchiness in snowpack produces marked, small-scale variation in flowering phenology. Plants in early- and late-melting patches are therefore likely to experience very different conditions during their flowering periods. Mertensia fusiformis is an early-flowering perennial that varies conspicuously in style length within and among populations. The hypothesis that style length represents an adaptation to local flowering time was tested. Specifically, it was hypothesized that lower air temperatures and higher frost risk would favour short-styled plants (with stigmas more shielded by corollas) in early-flowering patches, but that the pollen-collecting behaviour of flower visitors in late-flowering patches would favour long-styled plants.

METHODS

Floral morphology was measured, temperatures were monitored and pollinators were observed in several matched pairs of early and late populations. To evaluate effects of cold temperatures on plants of different style lengths, experimental pollinations were conducted during mornings (warm) and evenings (cool), and on flowers that either had or had not experienced a prior frost. The effectiveness of different pollinators was quantified as seed set following single visits to plants with relatively short or long styles.

KEY RESULTS

Late-flowering populations experienced warmer temperatures than early-flowering populations and a different suite of pollinators. Nectar-foraging bumble-bee queens and male solitary bees predominated in early populations, whereas pollen-collecting female solitary bees were more numerous in later sites. Pollinators differed significantly in their abilities to transfer pollen to stigmas at different heights, in accordance with our prediction. However, temperature and frost sensitivity did not differ between long- and short-styled plants. Although plants in late-flowering patches tended to have longer styles than those in early patches, this difference was not consistent.

CONCLUSIONS

Seasonal change in pollinator-mediated selection on style length may help maintain variation in this trait in M. fusiformis, but adaptation to local flowering time is not apparent. The prevalence of short styles in these populations requires further explanation.

Reciprocal style polymorphisms are not easily categorised: the case of heterostyly in Lithodora and Glandora (Boraginaceae)

Interest in reciprocal floral polymorphisms, such as heterostyly, has increased in recent decades because they can be used as suitable model systems to study mechanisms of outbreeding and disassortative mating in plants. Heterostylous plants are characterised by the presence of discrete morphs that differ in sex organ position and in some other ancillary traits. As regards sex organ deployment, different types of polymorphisms have been described, depending on number and type of discrete classes present in populations and degree of reciprocity between them. However, a clear-cut characterisation of stylar polymorphisms does not appear to be the best approach when there is great variability among populations because of continuous variation of some of traits examined. A recent study in Lithodora sensu lato (recently split into two separate genera, Lithodora and Glandora) showed a wide variation in sex organ position across species in the genus, which warrants precise population analysis of stylar polymorphism and its reciprocity. We provide a detailed morphometric analysis of flower sexual traits and include those considered to be ancillary characters. We report a wide variation in these traits in populations of Lithodora s.l. and highlight the subjectivity of the former characterisation of style polymorphism based on visual inspection. Ancillary traits appear repeatedly in Lithodora and Glandora, particularly in the latter. The appearance of these traits seems to be related to greater reciprocity between sexual whorls in Glandora, with the exception of G. prostrata. These results agree with evolutionary steps proposed in the build-up of heterostyly according to some evolutionary models. We also examined variation in polymorphisms in light of current models for evolution of heterostyly, and, more specifically, we sought to verify the prediction that flower traits as a whole (i.e., flower integration) respond to selective pressure to assure the exact location of pollen on the pollinator body. Most reciprocal populations and species, where between-morph pollen transfer is expected to be higher, would show greater integration. Our results confirm this hypothesis.

Expression of 3-hydroxy-3-methylglutaryl-CoA reductase, p-hydroxybenzoate-m-geranyltransferase and genes of phenylpropanoid pathway exhibits positive correlation with shikonins content in arnebia [Arnebia euchroma (Royle) Johnston]

BACKGROUND

Geranyl pyrophosphate (GPP) and p-hydroxybenzoate (PHB) are the basic precursors involved in shikonins biosynthesis. GPP is derived from mevalonate (MVA) and/or 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway(s), depending upon the metabolite and the plant system under consideration. PHB, however, is synthesized by only phenylpropanoid (PP) pathway. GPP and PHB are central moieties to yield shikonins through the synthesis of m-geranyl-p-hydroxybenzoate (GHB). Enzyme p-hydroxybenzoate-m-geranyltransferase (PGT) catalyses the coupling of GPP and PHB to yield GHB. The present research was carried out in shikonins yielding plant arnebia [Arnebia euchroma (Royle) Johnston], wherein no molecular work has been reported so far. The objective of the work was to identify the preferred GPP synthesizing pathway for shikonins biosynthesis, and to determine the regulatory genes involved in the biosynthesis of GPP, PHB and GHB.

RESULTS

A cell suspension culture-based, low and high shikonins production systems were developed to facilitate pathway identification and finding the regulatory gene. Studies with mevinolin and fosmidomycin, inhibitors of MVA and MEP pathway, respectively suggested MVA as a preferred route of GPP supply for shikonins biosynthesis in arnebia. Accordingly, genes of MVA pathway (eight genes), PP pathway (three genes), and GHB biosynthesis were cloned. Expression studies showed down-regulation of all the genes in response to mevinolin treatment, whereas gene expression was not influenced by fosmidomycin. Expression of all the twelve genes vis-à-vis shikonins content in low and high shikonins production system, over a period of twelve days at frequent intervals, identified critical genes of shikonins biosynthesis in arnebia.

CONCLUSION

A positive correlation between shikonins content and expression of 3-hydroxy-3-methylglutaryl-CoA reductase (AeHMGR) and AePGT suggested critical role played by these genes in shikonins biosynthesis. Higher expression of genes of PP pathway was a general feature for higher shikonins biosynthesis.

Differentiation of Symphytum species using RAPD and seed fatty acid patterns

RAPD and fatty acid profiles, and alpha-tocoferol contents of the seeds of some Symphytum species were analyzed for their differentiation. Discriminative patterns were found in the examined species. Genotypic characteristics of three Symphytum species were experienced using 20 decamer RAPD primers. The majority of band positions varied between species. The total amplified products of 20 RAPD primers was 247 (average of 12.35 bands per primer), of which 189 bands were polymorphic, corresponding to nearly 76.1% genetic diversity. The number of bands for each RAPD primer varied from 8 (UBC320) to 18 (OPA7). The percent of polymorphic bands ranged from 45.4% (OPJ20) to 93.3% (OPA16). Major unsaturated fatty acids in the seeds were linoleic, gamma-linolenic and oleic acids; alpha-linolenic and eicosenoic acids exhibited lower levels. Variations in quantities and total ratios of fatty acid groups, and alpha-tocopherol contents were also observed between the species, and the developmental stages of the seeds. Significant differences were found for the whole series of fatty acids between species, in addition to the calculated ratios (p < 0.05). Differences were also significant for fatty acid traits between two developmental stages of the seeds of endemic S. pseudobulbosum (p < 0.05). In some morphological descriptors of the seeds, difference was found at a significant level (p < 0.05). Obtained data based on genetic and biochemical variations seem to be useful for molecular delimitation of Symphytum, in addition to selection of the genotypes expressing a high amount of GLA.

Nitric oxide regulates shikonin formation in suspension-cultured Onosma paniculatum cells

Endogenously occurring nitric oxide (NO) is involved in the regulation of shikonin formation in Onosma paniculatum cells. NO generated after cells were inoculated into shikonin production medium reached the highest level after 2 d of culture, which was 16 times that at the beginning of the experiment, and maintained a high level for 6 d. A nitric oxide synthase (NOS) inhibitor, N(omega)-nitro-L-arginine (L-NNA), and a nitrate reductase (NR) inhibitor, sodium azide (SoA), consistent with their inhibition of NO biosynthesis, decreased shikonin formation significantly. This reduction could be alleviated or even abolished by exogenous NO supplied by sodium nitroprusside (SNP), suggesting that the inhibition of NO biosynthesis resulted in decreased shikonin formation. However, when endogenous NO biosynthesis was up-regulated by the elicitor from Rhizoctonia cerealis, shikonin production was enhanced further, showing a dependence on the elicitor-induced NO burst. Real-time PCR analysis showed that NO could significantly up-regulate the expression of PAL, PGT and HMGR, which encode key enzymes involved in shikonin biosynthesis. These results demonstrated that NO plays a critical role in shikonin formation in O. paniculatum cells.

Expression analysis of light-regulated genes isolated from a full-length-enriched cDNA library of Onosma paniculatum cell cultures

Shikonin and its derivatives are formed in large amounts in dark-cultured Onosma paniculatum cells. In order to isolate and identify the genes regulating shikonin biosynthesis, we constructed and characterized a full-length-enriched cDNA library of dark-cultured cells by using the SMART (Switching Mechanism At 5'-end of RNA Transcript) cDNA synthesis and LD-PCR (long-distance PCR) strategies. The titer of the primary cDNA library was 1.04 x 10(6)pfu/mL with a recombination rate of 99.60%. Most of the cDNA inserts ranged from 1.0 to 2.5 kb, and 78.33% of the 76 randomly selected clones contained full-length coding regions. Expression analysis of randomly selected genes by small scale microarray revealed that 23 genes were down-regulated, including 17 genes with known functions, 2 genes with putative functions, and 4 novel genes, and that 3 genes were up-regulated (two-fold) in cells cultured under white light as compared with those cultured in the dark. Interestingly, two of the down-regulated genes, encoding aci-reductone dioxygenase (ARD)-like protein and ethylene responsive factor (ERF), are involved in ethylene biosynthesis and signal transduction, implying that ethylene might play an important role as a signal molecule in light-regulated shikonin formation. These data contribute to a better understanding of light-involvement in regulating the formation of plant secondary metabolites.

A new approach to the quantification of degree of reciprocity in distylous (sensu lato) plant populations

BACKGROUND AND AIMS

Although evolution of sexual polymorphism has been traditionally analysed using discrete characters, most of these polymorphisms are continuous. This is the case of heterostyly. Heterostyly is a floral polymorphism successfully used as a model to study the evolution of the sexual systems in plants. It involves the reciprocal positioning of anthers and stigmas in flowers of different plants within the same population. Studies of the functioning of heterostyly require the quantification of the degree of reciprocity between morphs of heterostylous species. Some reciprocity indices have been proposed previously, but they show significant limitations that need to be dealt with. This paper analyses these existing indices, and proposes a new index that aims to avoid their main problems (e.g. takes into account population variability and offers a single value per population).

METHODS

The new index is based on the comparison of the position of every single sexual organ in the population with each and every organ of the opposite sex. To carry out all the calculations, a macro was programmed with MS Visual Basic in MS Excel. The behaviour of the index is tested using hypothetical data to simulate different situations of dimorphic populations; the index is also tested with some actual populations of different species of the genus Lithodora.

RESULTS AND CONCLUSIONS

The index of reciprocity proposed here is a sound alternative to previous indices: it compares stigma-stamen height gaps for all potential crosses in the population, it comprises stigma-stamen distance as well as dispersion, it is not skewed by the more frequent sex, and it can be meaningfully compared between populations and species. It has produced solid results for both hypothetical and natural populations.

Tissue distribution, core biosynthesis and diversification of pyrrolizidine alkaloids of the lycopsamine type in three Boraginaceae species

Three species of the Boraginaceae were studied: greenhouse-grown plants of Heliotropium indicum and Agrobacterium rhizogenes transformed roots cultures (hairy roots) of Cynoglossum officinale and Symphytum officinale. The species-specific pyrrolizidine alkaloid (PA) profiles of the three systems were established by GC-MS. All PAs are genuinely present as N-oxides. In H. indicum the tissue-specific PA distribution revealed the presence of PAs in all tissues with the highest levels in the inflorescences which in a flowering plant may account for more than 70% of total plant alkaloid. The sites of PA biosynthesis vary among species. In H. indicum PAs are synthesized in the shoot but not roots whereas they are only made in shoots for C. officinale and in roots of S. officinale. Classical tracer studies with radioactively labelled precursor amines (e.g., putrescine, spermidine and homospermidine) and various necine bases (trachelanthamidine, supinidine, retronecine, heliotridine) and potential ester alkaloid intermediates (e.g., trachelanthamine, supinine) were performed to evaluate the biosynthetic sequences. It was relevant to perform these comparative studies since the key enzyme of the core pathway, homospermidine synthase, evolved independently in the Boraginaceae and, for instance, in the Asteraceae [Reimann, A., Nurhayati, N., Backenkohler, A., Ober, D., 2004. Repeated evolution of the pyrrolizidine alkaloid-mediated defense system in separate angiosperm lineages. Plant Cell 16, 2772-2784.]. These studies showed that the core pathway for the formation of trachelanthamidine from putrescine and spermidine via homospermidine is common to the pathway in Senecio ssp. (Asteraceae). In both pathways homospermidine is further processed by a beta-hydroxyethylhydrazine sensitive diamine oxidase. Further steps of PA biosynthesis starting with trachelanthamidine as common precursor occur in two successive stages. Firstly, the necine bases are structurally modified and either before or after this modification are converted into their O(9)-esters by esterification with one of the stereoisomers of 2,3-dihydroxy-2-isopropylbutyric acid, the unique necic acid of PAs of the lycopsamine type. Secondly, the necine O(9)-esters may be further diversified by O(7)- and/or O(3')-acylation.

Molecular evidence for the age, origin, and evolutionary history of the American desert plant genus Tiquilia (Boraginaceae)

Although the deserts of North America are of very recent origin, their characteristic arid-adapted endemic plant lineages have been suggested to be much older. Earlier researchers have hypothesized that the ancestors of many of these modern desert lineages first adapted to aridity in highly localized arid or semi-arid sites as early as the late Cretaceous or early Tertiary, and that these lineages subsequently spread and diversified as global climate became increasingly arid during the Cenozoic. No study has explicitly examined these hypotheses for any North American arid-adapted plant group. The current paper tests these hypotheses using the genus Tiquilia (Boraginaceae), a diverse North American desert plant group. A strongly supported phylogeny of the genus is estimated using combined sequence data from three chloroplast markers (matK, ndhF, and rps16) and two nuclear markers (ITS and waxy). Ages of divergence events within the genus are estimated using penalized likelihood and a molecular clock approach on the ndhF tree for Tiquilia and representative outgroups, including most of the major lineages of Boraginales. The dating analysis suggests that the stem lineage of Tiquilia split from its nearest extant relative in the Paleocene or Eocene ( approximately 59-48 Ma). This was followed by a relatively long period before the first divergence in the crown group near the Eocene/Oligocene boundary ( approximately 33-29 Ma), shortly after the greatest Cenozoic episode of rapid aridification. Divergence of seven major lineages of Tiquilia is dated to the early-to-mid Miocene ( approximately 23-13 Ma). Several major lineages show a marked increase in diversification concomitant with the onset of more widespread semi-arid and then arid conditions beginning in the late Miocene ( approximately 7 Ma). This sequence of divergence events in Tiquilia agrees well with earlier researchers' ideas concerning North American desert flora assembly.

Inhibitory effect of the Agrobacterium rhizogenes rolC gene on rabdosiin and rosmarinic acid production in Eritrichium sericeum and Lithospermum erythrorhizon transformed cell cultures

Rabdosiin and related caffeic acid metabolites have been proposed as active pharmacological agents demonstrating potent anti-HIV and antiallergic activities. We transformed Eritrichium sericeum and Lithospermum erythrorhizon seedlings by the rolC gene, which has been recently described as an activator of plant secondary metabolism. Surprisingly, the rolC-transformed cell cultures of both plants yielded two- to threefold less levels of rabdosiin and rosmarinic acid (RA) than respective control cultures. This result establishes an interesting precedent when the secondary metabolites are differently regulated by a single gene. We show that the rolC gene affects production of rabdosiin and RA irrespective of the methyl jasmonate (MeJA)-mediated and the Ca(2+)-dependent NADPH oxidase pathways. Cantharidin, an inhibitor of serine/threonine phosphatases, partly diminishes the rolC-gene inhibitory effect that indicates involvement of the rolC-gene-mediated signal in plant regulatory controls, mediated by protein phosphatases. We also show that the control MeJA-stimulated E. sericeum root culture produces (-)-rabdosiin up to 3.41% dry weight, representing the highest level of this substance for plant cell cultures reported so far.

Biogeographic interpretation of splits graphs: least squares optimization of branch lengths

Although most often used to represent phylogenetic uncertainty, network methods are also potentially useful for describing the phylogenetic complexity expected to characterize recent species radiations. One network method with particular advantages in this context is split decomposition. However, in its standard implementation this approach is limited by a conservative criterion for branch length estimation. Here we extend the utility of split decomposition by introducing a least squares optimization technique for correcting branch lengths that may be underestimated by the standard implementation. This optimization of branch lengths is generally expected to improve divergence time estimates calculated from splits graphs. We illustrate the effect of least squares optimization on such estimates using the Australasian Myosotis and the Hawaiian silversword alliance as examples. We also discuss the biogeographic interpretation and limitations of splits graphs.

Life history traits in selfing versus outcrossing annuals: exploring the 'time-limitation' hypothesis for the fitness benefit of self-pollination

Background

Most self-pollinating plants are annuals. According to the 'time-limitation' hypothesis, this association between selfing and the annual life cycle has evolved as a consequence of strong r-selection, involving severe time-limitation for completing the life cycle. Under this model, selection from frequent density-independent mortality in ephemeral habitats minimizes time to flower maturation, with selfing as a trade-off, and / or selection minimizes the time between flower maturation and ovule fertilization, in which case selfing has a direct fitness benefit. Predictions arising from this hypothesis were evaluated using phylogenetically-independent contrasts of several life history traits in predominantly selfing versus outcrossing annuals from a data base of 118 species distributed across 14 families. Data for life history traits specifically related to maturation and pollination times were obtained by monitoring the start and completion of different stages of reproductive development in a greenhouse study of selfing and outcrossing annuals from an unbiased sample of 25 species involving five pair-wise family comparisons and four pair-wise genus comparisons.

Results

Selfing annuals in general had significantly shorter plant heights, smaller flowers, shorter bud development times, shorter flower longevity and smaller seed sizes compared with their outcrossing annual relatives. Age at first flower did not differ significantly between selfing and outcrossing annuals.

Conclusions

This is the first multi-species study to report these general life-history differences between selfers and outcrossers among annuals exclusively. The results are all explained more parsimoniously by selection associated with time-limitation than by selection associated with pollinator/mate limitation. The shorter bud development time reported here for selfing annuals is predicted explicitly by the time-limitation hypothesis for the fitness benefit of selfing (and not by the alternative 'reproductive assurance' hypothesis associated with pollinator/mate limitation). Support for the time-limitation hypothesis is also evident from published surveys: whereas selfers and outcrossers are about equally represented among annual species as a whole, selfers occur in much higher frequencies among the annual species found in two of the most severely time-limited habitats where flowering plants grow – deserts and cultivated habitats.

Molecular systematics of Boraginaceae tribe Boragineae based on ITS1 and trnL sequences, with special reference to Anchusa s.l

BACKGROUND AND AIMS

Boragineae is one of the main tribes of Boraginaceae, but delimitation and intergeneric classification of this group are unclear and have not yet been studied using DNA sequences. In particular, phylogenetic relationships in Anchusa s.l. still need to be elucidated in order to assess its taxonomic boundaries with respect to the controversial segregate genera Hormuzakia, Gastrocotyle, Phyllocara and Cynoglottis.

METHODS

Phylogenetic relationships among 51 taxa of tribe Boragineae were investigated by comparative sequencing of the trnL(UAA) intron of the plastid genome and of the ITS1 region of the nuclear ribosomal DNA. Exemplar taxa from 16 genera of Boragineae and all subgenera of Anchusa s.l. were included, along with two selected outgroups from tribes Lithospermeae and Cynoglosseae.

KEY RESULTS

Phylogenies generated by maximum parsimony and combined ITS1-trnL sequences support the monophyly of the tribe and a split into two clades, Pentaglottis and the remainder of Boragineae. The latter contains two large monophyletic groups. The first consists of three moderately to well-supported branches, Borago-Symphytum, Pulmonaria-Nonea and Brunnera. In the Pulmonaria-Nonea subclade, the rare endemic Paraskevia cesatiana is sister to Pulmonaria, and Nonea appears to be paraphyletic with respect to Elizaldia. The second main group corresponds to the well-supported clade of Anchusa s.l., with the megaphyllic, polyploid herb Trachystemon orientalis as sister taxon, although with low support. Anchusa s.l. is highly paraphyletic to its segregate genera and falls into four subclades: (1) Phyllocara, Hormuzakia, Anchusa subgenus Buglossum and A. subgenus Buglossoides; (2) Gastrocotyle; (3) A. subgenus Buglossellum and Cynoglottis; and (4) A. subgenus Anchusa, Lycopsis and Anchusella. All species of Anchusa subg. Anchusa, including the South African A. capensis, are included in a single unresolved clade. Anchusa subgenus Limbata is also included here despite marked divergence in floral morphology. The low nucleotide variation of ITS1 suggests a recent partly adaptive radiation within this group.

CONCLUSIONS

Molecular data show that nine of the usually accepted genera of the Boragineae consisting of two or more species are monophyletic: Anchusella, Borago, Brunnera, Cynoglottis, Gastrocotyle, Hormuzakia, Nonea, Pulmonaria and Symphytum. In addition, the tribe includes the four monotypic genera Paraskevia, Pentaglottis, Phyllocara and Trachystemon. The morphologically well-characterized segregate genera in Anchusa s.l. are all confirmed by DNA sequences and should be definitively accepted. Most of the traditionally recognized subgenera of Anchusa are also supported as monophyletic groups by both nuclear and plastid sequence data. In order to bring taxonomy in line with phylogeny, the institution of new, independent generic entities for subgenera Buglossum, Buglossellum and Buglossoides and a narrower but more natural concept of Anchusa are advocated.