The sweet tabaiba or there and back again: phylogeographical history of the Macaronesian Euphorbia balsamifera

BACKGROUND AND AIMS

Biogeographical relationships between the Canary Islands and north-west Africa are often explained by oceanic dispersal and geographical proximity. Sister-group relationships between Canarian and eastern African/Arabian taxa, the 'Rand Flora' pattern, are rare among plants and have been attributed to the extinction of north-western African populations. Euphorbia balsamifera is the only representative species of this pattern that is distributed in the Canary Islands and north-west Africa; it is also one of few species present in all seven islands. Previous studies placed African populations of E. balsamifera as sister to the Canarian populations, but this relationship was based on herbarium samples with highly degraded DNA. Here, we test the extinction hypothesis by sampling new continental populations; we also expand the Canarian sampling to examine the dynamics of island colonization and diversification.

METHODS

Using target enrichment with genome skimming, we reconstructed phylogenetic relationships within E. balsamifera and between this species and its disjunct relatives. A single nucleotide polymorphism dataset obtained from the target sequences was used to infer population genetic diversity patterns. We used convolutional neural networks to discriminate among alternative Canary Islands colonization scenarios.

KEY RESULTS

The results confirmed the Rand Flora sister-group relationship between western E. balsamifera and Euphorbia adenensis in the Eritreo-Arabian region and recovered an eastern-western geographical structure among E. balsamifera Canarian populations. Convolutional neural networks supported a scenario of east-to-west island colonization, followed by population extinctions in Lanzarote and Fuerteventura and recolonization from Tenerife and Gran Canaria; a signal of admixture between the eastern island and north-west African populations was recovered.

CONCLUSIONS

Our findings support the Surfing Syngameon Hypothesis for the colonization of the Canary Islands by E. balsamifera, but also a recent back-colonization to the continent. Populations of E. balsamifera from northwest Africa are not the remnants of an ancestral continental stock, but originated from migration events from Lanzarote and Fuerteventura. This is further evidence that oceanic archipelagos are not a sink for biodiversity, but may be a source of new genetic variability.

Phylogenetics and phylogeography of Euphorbia canariensis reveal an extreme Canarian-Asian disjunction but limited inter-island colonization

Euphorbia canariensis is an iconic endemic species representative of the lowland xerophytic communities of the Canary Islands. It is widely distributed in the archipelago despite having diasporas unspecialized for long-distance dispersal. Here, we reconstructed the evolutionary history of E. canariensis at two levels: a time-calibrated phylogenetic analysis aimed at clarifying interspecific relationships and large-scale biogeographic patterns; and a phylogeographic study focused on the history of colonization across the Canary Islands. For the phylogenetic study, we sequenced the ITS region for E. canariensis and related species of Euphorbia sect. Euphorbia. For the phylogeographic study, we sequenced two cpDNA regions for 28 populations representing the distribution range of E. canariensis. The number of inter-island colonization events was explored using PAICE, a recently developed method that includes a sample size correction. Additionally, we used species distribution modelling (SDM) to evaluate environmental suitability for E. canariensis through time. Phylogenetic results supported a close relationship between E. canariensis and certain Southeast Asian species (E. epiphylloides, E. lacei, E. sessiliflora). In the Canaries, E. canariensis displayed a west-to-east colonization pattern, not conforming to the "progression rule", i.e. the concordance between phylogeographic patterns and island emergence times. We estimated between 20 and 50 inter-island colonization events, all of them in the Quaternary, and SDM suggested a late Quaternary increase in environmental suitability for E. canariensis. The extreme biogeographic disjunction between Macaronesia and Southeast Asia (ca. 11,000 km) parallels that found in a few other genera (Pinus, Dracaena). We hypothesize that these disjunctions are better explained by extinction across north Africa and southwest Asia rather than long-distance dispersal. The relatively low number of inter-island colonization events across the Canaries is congruent with the low dispersal capabilities of E. canariensis.

Identification of One O-Methyltransferase Gene Involved in Methylated Flavonoid Biosynthesis Related to the UV-B Irradiation Response in Euphorbia lathyris

Flavonoids are ubiquitous polyphenolic compounds that play a vital role in plants' defense response and medicinal efficacy. UV-B radiation is a vital environmental regulator governing flavonoid biosynthesis in plants. Many plants rapidly biosynthesize flavonoids as a response to UV-B stress conditions. Here, we investigated the effects of flavonoid biosynthesis via UV-B irradiation in Euphorbia lathyris. We found that exposure of the E. lathyris callus to UV-B radiation sharply increased the level of one O-methyltransferase (ElOMT1) transcript and led to the biosynthesis of several methylated flavonoids. The methyltransferase ElOMT1 was expressed heterologously in E. coli, and we tested the catalytic activity of recombinant ElOMT1 with possible substrates, including caffeic acid, baicalin, and luteolin, in vitro. ElOMT1 could efficiently methylate when the hydroxyl groups were contained in the core nucleus of the flavonoid. This molecular characterization identifies a methyltransferase responsible for the chemical modification of the core flavonoid structure through methylation and helps reveal the mechanism of methylated flavonoid biosynthesis in Euphorbiaceae. This study identifies the O-methyltransferase that responds to UV-B irradiation and helps shed light on the mechanism of flavonoid biosynthesis in Euphorbia lathyris.

Identification and Functional Characterization of a Novel Sinapyl Alcohol Acyltransferase from Euphorbia lathyris L

Methoxyeugenol is a phenylpropene compound derived from plants and has various bioactivities. The chemical synthesis of methoxyeugenol is accompanied by pollution issues, whereas extraction from plants is associated with problems such as low yield and high cost. The production of methoxyeugenol can be effectively addressed through an enzymatic approach. In this study, the acyltransferase genes of Euphorbia lathyris L. were screened by homologous alignment of the transcriptome data of E. lathyris in the late growth stage and the acyltransferase genes of the closely related plant species. The results showed that ElBAHD10 had the closest relationship with earlier reported ScCFAT and PhCFAT, which were found to catalyze the reaction of coniferyl alcohol to generate coniferyl acetate. The ElBAHD10 gene was successfully cloned from E. lathyris and subsequently expressed in Escherichia coli. The purified protein ElBAHD10 catalyzed the reaction of sinapyl alcohol with acetyl CoA and cinnamoyl CoA to form sinapyl acetate and sinapyl cinnamate, respectively. In contrast, the crude ElBAHD10 protein could catalyze sinapyl alcohol to directly generate methoxyeugenol. The recombinant E. coli strain expressing ElBAHD10 produced methoxyeugenol through whole-cell transformation. This study provides insights and lays the foundation for methoxyeugenol production through biosynthetic approaches.

Evolutionary dynamics of Euphorbia carniolica suggest a complex Plio-Pleistocene history of understorey species of deciduous forest in southeastern Europe

Deciduous forests form the dominant natural vegetation of Europe today, but were restricted to small refugia during Pleistocene cold stages, implying an evolutionary past shaped by recurrent range contractions and expansions. Cold-stage forest refugia were probably widespread in southern and central Europe, with the northwestern Balkan Peninsula being of particular importance. However, the actual number and location of deciduous forest refugia, as well as the connections between them, remain disputed. Here, we address the evolutionary dynamics of the deciduous forest understorey species Euphorbia carniolica as a proxy for past forest dynamics. To do so, we obtained genomic and morphometric data from populations representing the species' entire range, investigated phylogenetic position and intraspecific genetic variation, tested explicit demographic scenarios and applied species distribution models. Our data support two disjoint groups linked to separate refugia on the northwestern and central Balkan Peninsula. We find that genetic differentiation between groups started in the early Pleistocene via vicariance, suggesting a larger distribution in the past. Both refugia acted as sources for founder events to the southeastern Alps and the Carpathians; the latter were likely colonised before the last cold stage. In line with traditional views on the pre-Pleistocene origin of many southeastern European deciduous forest species, the origin of E. carniolica was dated to the late Pliocene. The fact that E. carniolica evolved at a time when a period of continuous forestation was ending in much of Eurasia provides an interesting biogeographical perspective on the past links between Eurasian deciduous forests and their biota.

Identifying Terpenoid Biosynthesis Genes in Euphorbia maculata via Full-Length cDNA Sequencing

The annual herb Euphorbia maculata L. produces anti-inflammatory and biologically active substances such as triterpenoids, tannins, and polyphenols, and it is used in traditional Chinese medicine. Of these bioactive compounds, terpenoids, also called isoprenoids, are major secondary metabolites in E. maculata. Full-length cDNA sequencing was carried out to characterize the transcripts of terpenoid biosynthesis reference genes and determine the copy numbers of their isoforms using PacBio SMRT sequencing technology. The Illumina short-read sequencing platform was also employed to identify differentially expressed genes (DEGs) in the secondary metabolite pathways from leaves, roots, and stems. PacBio generated 62 million polymerase reads, resulting in 81,433 high-quality reads. From these high-quality reads, we reconstructed a genome of 20,722 genes, in which 20,246 genes (97.8%) did not have paralogs. About 33% of the identified genes had two or more isoforms. DEG analysis revealed that the expression level differed among gene paralogs in the leaf, stem, and root. Whole sets of paralogs and isoforms were identified in the mevalonic acid (MVA), methylerythritol phosphate (MEP), and terpenoid biosynthesis pathways in the E. maculata L. The nucleotide information will be useful for identifying orthologous genes in other terpenoid-producing medicinal plants.

Evolution of EPSPS double mutation imparting glyphosate resistance in wild poinsettia (Euphorbia heterophylla L.)

The evolution of glyphosate resistance (GR) in weeds is an increasing problem. Glyphosate has been used intensively on wild poinsettia (Euphorbia heterophylla L.) populations for at least 20 years in GR crops within South America. We investigated the GR mechanisms in a wild poinsettia population from a soybean field in southern Brazil. The GR population required higher glyphosate doses to achieve 50% control (LD50) and 50% dry mass reduction (MR50) compared to a glyphosate susceptible (GS) population. The ratio between the LD50 and MR50 of GR and GS resulted in resistance factors (RF) of 6.9-fold and 6.1-fold, respectively. Shikimate accumulated 6.7 times more in GS than in GR when leaf-discs were incubated with increasing glyphosate concentrations. No differences were found between GR and GS regarding non-target-site mechanisms. Neither population metabolized glyphosate to significant levels following treatment with 850 g ha-1 glyphosate. Similar levels of 14C-glyphosate uptake and translocation were observed between the two populations. No differences in EPSPS expression were found between GS and GR. Two target site mutations were found in all EPSPS alleles of homozygous resistant plants: Thr102Ile + Pro106Thr (TIPT-mutation). Heterozygous individuals harbored both alleles, wild-type and TIPT. Half of GR individuals were heterozygous, suggesting that resistance is still evolving in the population. A genotyping assay was developed based on the Pro106Thr mutation, demonstrating high efficiency to identify homozygous, heterozygous or wild-type EPSPS sequences across different plants. This is the first report of glyphosate-resistant wild-poinsettia harboring an EPSPS double mutation (TIPT) in the same plant.

The Chromosome-Based Rubber Tree Genome Provides New Insights into Spurge Genome Evolution and Rubber Biosynthesis

The rubber tree, Hevea brasiliensis, produces natural rubber that serves as an essential industrial raw material. Here, we present a high-quality reference genome for a rubber tree cultivar GT1 using single-molecule real-time sequencing (SMRT) and Hi-C technologies to anchor the ∼1.47-Gb genome assembly into 18 pseudochromosomes. The chromosome-based genome analysis enabled us to establish a model of spurge chromosome evolution, since the common paleopolyploid event occurred before the split of Hevea and Manihot. We show recent and rapid bursts of the three Hevea-specific LTR-retrotransposon families during the last 10 million years, leading to the massive expansion by ∼65.88% (∼970 Mbp) of the whole rubber tree genome since the divergence from Manihot. We identify large-scale expansion of genes associated with whole rubber biosynthesis processes, such as basal metabolic processes, ethylene biosynthesis, and the activation of polysaccharide and glycoprotein lectin, which are important properties for latex production. A map of genomic variation between the cultivated and wild rubber trees was obtained, which contains ∼15.7 million high-quality single-nucleotide polymorphisms. We identified hundreds of candidate domestication genes with drastically lowered genomic diversity in the cultivated but not wild rubber trees despite a relatively short domestication history of rubber tree, some of which are involved in rubber biosynthesis. This genome assembly represents key resources for future rubber tree research and breeding, providing novel targets for improving plant biotic and abiotic tolerance and rubber production.

iTRAQ-Based Proteomics Analysis of Autophagy-Mediated Responses against MeJA in Laticifers of Euphorbia kansui L

Autophagy is a well-defined catabolic mechanism whereby cytoplasmic materials are engulfed into a structure termed the autophagosome. Methyl jasmonate (MeJA), a plant hormone, mediates diverse developmental process and defense responses which induce a variety of metabolites. In plants, little is known about autophagy-mediated responses against MeJA. In this study, we used high-throughput comparative proteomics to identify proteins of latex in the laticifers. The isobaric tags for relative and absolute quantification (iTRAQ) MS/MS proteomics were performed, and 298 proteins among MeJA treated groups and the control group of Euphorbia kansui were identified. It is interesting to note that 29 significant differentially expressed proteins were identified and their associations with autophagy and ROS pathway were verified for several selected proteins as follows: α-L-fucosidase, β-galactosidase, cysteine proteinase, and Cu/Zn superoxide dismutase. Quantitative real-time PCR analysis of the selected genes confirmed the fact that MeJA might enhance the expression of some genes related to autophagy. The western blotting and immunofluorescence results of ATG8 and ATG18a which are two important proteins for the formation of autophagosomes also demonstrated that MeJA could promote autophagy at the protein level. Using the electron microscope, we observed an increase in autophagosomes after MeJA treatment. These results indicated that MeJA might promote autophagy in E. kansui laticifers; and it was speculated that MeJA mediated autophagy through two possible ways: the increase of ROS induces ATG8 accumulation and then aotophagosome formation, and MeJA promotes ATG18 accumulation and then autophagosome formation. Taken together, our results provide several novel insights for understanding the mechanism between autophagy and MeJA treatment. However, the specific mechanism remains to be further studied in the future.

Convergent evolution in floral morphology in a plant ring species, the Caribbean Euphorbia tithymaloides

PREMISE

Ring species have long fascinated evolutionary biologists for their potential insights into lineage divergence and speciation across space. Few studies have investigated the potential for convergent or parallel evolution along the diverging fronts of ring species. We investigated a potential case of parallel floral variation in the Caribbean spurge Euphorbia tithymaloides, the only plant system with molecular support as a ring species. The terminal populations of each front, despite being the most divergent, exhibit such similar floral traits that they were originally considered each other's closest relative.

METHODS

We evaluated convergence in floral and leaf traits in relation to geography across 95 populations spanning the distribution of E. tithymaloides. We also reanalyzed available genetic data (from previous phylogenetic analyses) in an explicitly spatial framework.

RESULTS

Floral morphology appears to have shifted in a convergent fashion along both geographic fronts of E. tithymaloides, resulting in shorter and more compact inflorescences in Antillean populations compared to the typical elongate "slipper-like" cyathia characteristic of the area of origin. Patterns of spatial genetic variation were more consistent with a two-fronted invasion of the Caribbean than with a simpler model of isolation-by-distance.

CONCLUSIONS

Floral divergence in E. tithymaloides is consistent with convergent evolution along the two fronts of a ring species. We outline several (not mutually exclusive) mechanisms that could be driving patterns in morphology, including shifts toward generalized pollination with reduced reliance on hummingbirds, shifts in floral structure closely matching available hummingbird bill traits, and shifts toward increased selfing.

MicroRNAome Profile of Euphorbia kansui in Response to Methyl Jasmonate

miRNAs play vital regulatory roles in different plant developmental stages and in plant response to biotic and abiotic stresses. However, information is limited on the miRNA regulatory mechanism to methyl jasmonate (MeJA). In this study, we used the microRNAome profile to illustrate the relevant regulatory mechanisms of Euphorbia kansui in response to methyl jasmonate (MeJA) through Illumina RNA-Seq. As a result, we identified 875 miRNAs corresponding to 11,277 target mRNAs, among them, 168 known miRNA families representing 6019 target mRNAs sequences were obtained. 452 miRNA-mRNA pairs presented an anti-correlationship (Cor < −0.50 and p-value of correlation ≤ 0.05). The miRNA with a fold change ≥ 2 and a p (p-Value) < 0.05 in pairwise comparison were identified as significant differentially expressed miRNAs (DEMs). The DEMs in MeJA treatment of 0, 24, 36 and 48 h were compared by using Short Time Expression Miner (STEM) cluster and 4 significant gene profiles (p-value ≤ 0.02) were identified. Through the kyoto encyclopedia of genes and genomes (KEGG) pathway and gene ontology (GO) enrichment analysis on all miRNA targets, we identified 33 mRNAs in terpenoid biosynthesis, which were regulated by miRNAs under MeJA treatment, so the miRNA maybe involved in the response of E. kansui plant to exogenous MeJA and the results would provide very useful information on illustrating the regulatory mechanism of E. kansui and also provide an overall view of the miRNAs response to MeJA stress of a non-model plant.

The rare orange-red colored Euphorbia pulcherrima cultivar 'Harvest Orange' shows a nonsense mutation in a flavonoid 3'-hydroxylase allele expressed in the bracts

BACKGROUND

Commercially available poinsettia (Euphorbia pulcherrima) varieties prevalently accumulate cyanidin derivatives and show intense red coloration. Orange-red bract color is less common. We investigated four cultivars displaying four different red hues with respect to selected enzymes and genes of the anthocyanin pathway, putatively determining the color hue.

RESULTS

Red hues correlated with anthocyanin composition and concentration and showed common dark red coloration in cultivars 'Christmas Beauty' and 'Christmas Feeling' where cyanidin derivatives were prevalent. In contrast, orange-red bract color is based on the prevalent presence of pelargonidin derivatives that comprised 85% of the total anthocyanin content in cv. 'Premium Red' and 96% in cv. 'Harvest Orange' (synonym: 'Orange Spice'). cDNA clones of flavonoid 3'-hydroxylase (F3'H) and dihydroflavonol 4-reductase (DFR) were isolated from the four varieties, and functional activity and substrate specificity of the corresponding recombinant enzymes were studied. Kinetic studies demonstrated that poinsettia DFRs prefer dihydromyricetin and dihydroquercetin over dihydrokaempferol, and thus, favor the formation of cyanidin over pelargonidin. Whereas the F3'H cDNA clones of cultivars 'Christmas Beauty', 'Christmas Feeling', and 'Premium Red' encoded functionally active enzymes, the F3'H cDNA clone of cv. 'Harvest Orange' contained an insertion of 28 bases, which is partly a duplication of 20 bases found close to the insertion site. This causes a frameshift mutation with a premature stop codon after nucleotide 132 and, therefore, a non-functional enzyme. Heterozygosity of the F3'H was demonstrated in this cultivar, but only the mutated allele was expressed in the bracts. No correlation between F3'H-expression and the color hue could be observed in the four species.

CONCLUSIONS

Rare orange-red poinsettia hues caused by pelargonidin based anthocyanins can be achieved by different mechanisms. F3'H is a critical step in the establishment of orange red poinsettia color. Although poinsettia DFR shows a low substrate specificity for dihydrokaempferol, sufficient precursor for pelargonidin formation is available in planta, in the absence of F3'H activity.

LOL2 and LOL5 loci control latex production by laticifer cells in Euphorbia lathyris

Laticifers are specialized plant cells capable of indefinite elongation that ramify extensively and are responsible for latex biosynthesis and accumulation. However, the mechanisms underlying laticifer cell differentiation, growth and production of latex remain largely unknown. In a search for mutants showing enhanced accumulation of latex we identified two LOT OF LATEX (LOL) loci in Euphorbia lathyris. lol2 and lol5 mutants show enhanced production of latex contained within laticifer cells. The recessive lol2 mutant carries increased biosynthesis of the plant hormone jasmonoyl-isoleucine (JA-Ile) and therefore establishes a genetic link between jasmonic acid (JA) signaling and latex production in laticifers. Instead, heightened production of latex in lol5 plants obeys to enhanced proliferation of laticifer cells. Phylogenetic analysis of laticifer-expressed genes in E. lathyris and in two other latex-bearing species, Euphorbia corallioides and Euphorbia palustris, allowed the identification of canonical JA responsive elements present in the gene promoter regions of laticifer marker genes. Moreover, we identified that the hormone JA functions not as a morphogen for laticifer differentiation but as a trigger for the fill out of laticifers with latex and the associated triterpenoids. The identification of LOL loci represents a further step towards the understanding of mechanisms controlling latex production in laticifer cells.

Foliar Glyphosate Treatment Alters Transcript and Hormone Profiles in Crown Buds of Leafy Spurge and Induces Dwarfed and Bushy Phenotypes throughout its Perennial Lifecycle

Leafy spurge ( L.) is an invasive weed of North America and its perennial nature attributed to underground adventitious buds (UABs) that undergo seasonal cycles of para-, endo-, and ecodormancy. Recommended rates of glyphosate (∼1 kg ha) destroy aboveground shoots but plants still regenerate vegetatively; therefore, it is considered glyphosate-tolerant. However, foliar application of glyphosate at higher rates (2.2-6.7 kg ha) causes sublethal effects that induce UABs to produce stunted, bushy phenotypes. We investigated the effects of glyphosate treatment (±2.24 kg ha) on vegetative growth, phytohormone, and transcript profiles in UABs under controlled environments during one simulated seasonal cycle. Because shoots derived from UABs of foliar glyphosate-treated plants produced stunted, bushy phenotypes, we could not directly determine if these UABs transitioned through seasonally induced endo- and ecodormancy. However, transcript abundance for leafy spurge dormancy marker genes and principal component analyses suggested that UABs of foliar glyphosate-treated plants transitioned through endo- and ecodormancy. Glyphosate treatment increased shikimate abundance in UABs 7 d after treatment; however, the abundance of shikimate gradually decreased as UABs transitioned through endo- and ecodormancy. The dissipation of shikimate over time suggests that glyphosate's target site was no longer affected, but these changes did not reverse the altered phenotypes observed from UABs of foliar glyphosate-treated leafy spurge. Transcript profiles further indicated that foliar glyphosate treatment significantly affected phytohormone biosynthesis and signaling, particularly auxin transport; gibberellic acid, abscisic acid and jasmonic acid biosynthesis; ethylene responses; and detoxification and cell cycle processes in UABs. These results correlated well with the available phytohormone profiles and altered phenotypes.

Transcriptome and proteome analyses provide insight into laticifer's defense of Euphorbia tirucalli against pests

The cytoplasm of laticifers, which are plant cells specialized for rubber production and defense against microbes and herbivores, is a latex. Although laticifers share common functions, the protein constituents of latexes are highly variable among plant species and even among organs. In this study, transcriptomic and proteomic analyses of Euphorbia tirucalli's (Euphorbiaceae) latex were conducted to determine the molecular basis of the laticifer's functions in this plant. The hybrid de novo assembly of Illumina mRNA-seq and expressed sequence tags obtained by Sanger's sequencing revealed 26,447 unigenes. A unigene similar to Arabidopsis embryo-specific protein 3 (AT5G62200), which is a PLAT domain-containing protein, and rubber elongation factor showed the highest expression levels. The proteome analysis, studied by liquid chromatography-mass spectrometry with the de novo assembled unigenes as the database, revealed 161 proteins in the latex, 107 of which were not detected in the stem. A gene ontology analysis indicated that the laticifer's proteome was enriched with proteins related to proteolysis, phosphatase, defense against various environmental stresses and lipid metabolisms. D-mannose-binding lectin, ricin (which lacked the N-terminal conserved ribosome-inactivating protein domain), chitinase and peroxidase were highly accumulated, as confirmed by two-dimensional polyacrylamide gel electrophoresis. Thus, the lectins and chitinase may be the major defensive proteins against pests, and the other defense-related proteins and transcripts detected in latex may work in coordination with them. Highly expressing unigenes with unknown functions are candidate novel defense- or rubber production-related genes.

Novel Insights into the Organization of Laticifer Cells: A Cell Comprising a Unified Whole System

Laticifer cells are specialized plant cells that synthesize and accumulate latex. Studies on laticifers have lagged behind in recent years, and data regarding the functional role of laticifers and their fitness benefit still remain elusive. Laticifer differentiation and its impact on plant growth and development also remain to be investigated. Here, cellular, molecular, and genetic tools were developed to examine the distribution, differentiation, ontogeny, and other characteristic features, as well as the potential developmental role of laticifer cells in the latex-bearing plant Euphorbia lathyris. The organization of the laticiferous system within the E. lathyris plant body is reported, emerging as a single elongated and branched coenocytic cell, constituting the largest cell type existing in plants. We also report the ontogeny and organization of laticifer cells in the embryo and the identification of a laticifer-associated gene expression pattern. Moreover, the identification of laticifer- and latex-deficient mutants (pil mutants) allowed for the identification of distinct loci regulating laticifer differentiation, growth, and metabolic activity. Additionally, pil mutants revealed that laticifer cells appear nonessential for plant growth and development, thus pointing toward their importance, instead, for specific ecophysiological adaptations of latex-bearing plants in natural environments.

Effect of endophytic Bacillus cereus ERBP inoculation into non-native host: Potentials and challenges for airborne formaldehyde removal

Phytoremediation could be a cost-effective, environmentally friendly approach for the treatment of indoor air. However, some drawbacks still dispute the expediency of phytotechnology. Our objectives were to investigate the competency of plant growth-promoting (PGP) endophytic Bacillus cereus ERBP (endophyte root blue pea), isolated from the root of Clitoria ternatea, to colonize and stabilize within Zamioculcas zamiifolia and Euphorbia milii as non-native hosts without causing any disease or stress symptoms. Moreover, the impact of B. cereus ERBP on the natural shoot endophytic community and for the airborne formaldehyde removal capability of non-native hosts was assessed. Non-native Z. zamiifolia was effectively inoculated with B. cereus ERBP through soil as the most efficient method of endophyte inoculation. Denaturing gradient gel electrophoresis profiling of the shoot endophytic community verified the colonization and stability of B. cereus ERBP within its non-native host during a 20-d fumigation period without interfering with the natural shoot endophytic diversity of Z. zamiifolia. B. cereus ERBP conferred full protection to its non-native host against formaldehyde phytotoxicity and enhanced airborne formaldehyde removal of Z. zamiifolia whereas non-inoculated plants suffered from formaldehyde phytotoxicity because their natural shoot endophytic community was detrimentally affected by formaldehyde. In contrast, B. cereus ERBP inoculation into non-native E. milii deteriorated airborne formaldehyde removal of the non-native host (compared to a non-inoculated one) as B. cereus ERBP interfered with natural shoot endophytic community of E. milii, which caused stress symptoms and stimulated ethylene biosynthesis. Non-native host inoculation with PGP B. cereus ERBP could bear potentials and challenges for airborne formaldehyde removal.

Oxidation and cyclization of casbene in the biosynthesis of Euphorbia factors from mature seeds of Euphorbia lathyris L

The seed oil of Euphorbia lathyris L. contains a series of macrocyclic diterpenoids known as Euphorbia factors. They are the current industrial source of ingenol mebutate, which is approved for the treatment of actinic keratosis, a precancerous skin condition. Here, we report an alcohol dehydrogenase-mediated cyclization step in the biosynthetic pathway of Euphorbia factors, illustrating the origin of the intramolecular carbon-carbon bonds present in lathyrane and ingenane diterpenoids. This unconventional cyclization describes the ring closure of the macrocyclic diterpene casbene. Through transcriptomic analysis of E. lathyris L. mature seeds and in planta functional characterization, we identified three enzymes involved in the cyclization route from casbene to jolkinol C, a lathyrane diterpene. These enzymes include two cytochromes P450 from the CYP71 clan and an alcohol dehydrogenase (ADH). CYP71D445 and CYP726A27 catalyze regio-specific 9-oxidation and 5-oxidation of casbene, respectively. When coupled with these P450-catalyzed monooxygenations, E. lathyris ADH1 catalyzes dehydrogenation of the hydroxyl groups, leading to the subsequent rearrangement and cyclization. The discovery of this nonconventional cyclization may provide the key link to complete elucidation of the biosynthetic pathways of ingenol mebutate and other bioactive macrocyclic diterpenoids.

Glyphosate's impact on vegetative growth in leafy spurge identifies molecular processes and hormone cross-talk associated with increased branching

BACKGROUND

Leafy spurge (Euphorbia esula) is a perennial weed that is considered glyphosate tolerant, which is partially attributed to escape through establishment of new vegetative shoots from an abundance of underground adventitious buds. Leafy spurge plants treated with sub-lethal concentrations of foliar-applied glyphosate produce new vegetative shoots with reduced main stem elongation and increased branching. Processes associated with the glyphosate-induced phenotype were determined by RNAseq using aerial shoots derived from crown buds of glyphosate-treated and -untreated plants. Comparison between transcript abundance and accumulation of shikimate or phytohormones (abscisic acid, auxin, cytokinins, and gibberellins) from these same samples was also done to reveal correlations.

RESULTS

Transcriptome assembly and analyses confirmed differential abundance among 12,918 transcripts (FDR ≤ 0.05) and highlighted numerous processes associated with shoot apical meristem maintenance and stem growth, which is consistent with the increased number of actively growing meristems in response to glyphosate. Foliar applied glyphosate increased shikimate abundance in crown buds prior to decapitation of aboveground shoots, which induces growth from these buds, indicating that 5-enolpyruvylshikimate 3-phosphate (EPSPS) the target site of glyphosate was inhibited. However, abundance of shikimate was similar in a subsequent generation of aerial shoots derived from crown buds of treated and untreated plants, suggesting EPSPS is no longer inhibited or abundance of shikimate initially observed in crown buds dissipated over time. Overall, auxins, gibberellins (precursors and catabolites of bioactive gibberellins), and cytokinins (precursors and bioactive cytokinins) were more abundant in the aboveground shoots derived from glyphosate-treated plants.

CONCLUSION

Based on the overall data, we propose that the glyphosate-induced phenotype resulted from complex interactions involving shoot apical meristem maintenance, hormone biosynthesis and signaling (auxin, cytokinins, gibberellins, and strigolactones), cellular transport, and detoxification mechanisms.

Coordinated Expression of FLOWERING LOCUS T and DORMANCY ASSOCIATED MADS-BOX-Like Genes in Leafy Spurge

Leafy spurge (Euphorbia esula L.) is a noxious perennial weed that produces underground adventitious buds, which are crucial for generating new vegetative shoots following periods of freezing temperatures or exposure to various control measures. It is also capable of flowering and producing seeds, but requires vernalization in some cases. DORMANCY ASSOCIATED MADS-BOX (DAM) genes have been proposed to play a direct role in the transition to winter-induced dormancy and maintenance through regulation of the FLOWERING LOCUS T (FT) gene, which also is likely involved in the vernalization process. To explore the regulation of FT and DAM during dormancy transitions in leafy spurge, the transcript accumulation of two previously cloned DAM splice variants and two different previously cloned FT genes was characterized. Under long-photoperiods (16 h light), both DAM and FT transcripts accumulate in a diurnal manner. Tissue specific expression patterns indicated the tissues with high DAM expression had low FT expression and vice versa. DAM expression is detected in leaves, stems, shoot tips, and crown buds. FT transcripts were detected mainly in leaves and flowers. Under dormancy inducing conditions, DAM and FT genes had an inverse expression pattern. Additionally, chromatin immunoprecipitation assays were performed using DAM-like protein specific antibodies to demonstrate that DAM or related proteins likely bind to cryptic and/or conserved CArG boxes in the promoter regions of FT genes isolated from endodormant crown buds. These results are consistent with the hypothesis that DAM proteins play a crucial role in leafy spurge dormancy transition and maintenance, potentially by negatively regulating the expression of FT.

Euphorbia characias latex: micromorphology of rubber particles and rubber transferase activity

We have recently characterized a natural rubber in the latex of Euphorbia characias. Following that study, we here investigated the rubber particles and rubber transferase in that Mediterranean shrub. Rubber particles, observed by scanning electron microscopy, are spherical in shape with diameter ranging from 0.02 to 1.2 μm. Washed rubber particles exhibit rubber transferase activity with a rate of radiolabeled [(14)C]IPP incorporation of 4.5 pmol min(-1)mg(-1). Denaturing electrophoresis profile of washed rubber particles reveals a single protein band of 37 kDa that is recognized in western blot analysis by antibodies raised against the synthetic peptide whose sequence, DVVIRTSGETRLSNF, is included in one of the five regions conserved among cis-prenyl chain elongation enzymes. The cDNA nucleotide sequence of E. characias rubber transferase (GenBank JX564541) and the deduced amino acid sequence appear to be highly homologous to the sequence of several plant cis-prenyltransferases.

Dehydration-induced endodormancy in crown buds of leafy spurge highlights involvement of MAF3- and RVE1-like homologs, and hormone signaling cross-talk

Vegetative shoot growth from underground adventitious buds of leafy spurge is critical for survival of this invasive perennial weed after episodes of severe abiotic stress. To determine the impact that dehydration-stress has on molecular mechanisms associated with vegetative reproduction of leafy spurge, greenhouse plants were exposed to mild- (3-day), intermediate- (7-day), severe- (14-day) and extended- (21-day) dehydration treatments. Aerial tissues of treated plants were then decapitated and soil was rehydrated to determine the growth potential of underground adventitious buds. Compared to well-watered plants, mild-dehydration accelerated new vegetative shoot growth, whereas intermediate- through extended-dehydration treatments both delayed and reduced shoot growth. Results of vegetative regrowth further confirmed that 14 days of dehydration induced a full-state of endodormancy in crown buds, which was correlated with a significant (P < 0.05) change in abundance of 2,124 transcripts. Sub-network enrichment analyses of transcriptome data obtained from the various levels of dehydration treatment also identified central hubs of over-represented genes involved in processes such as hormone signaling (i.e., ABA, auxin, ethylene, GA, and JA), response to abiotic stress (DREB1A/2A, RD22) and light (PIF3), phosphorylation (MPK4/6), circadian regulation (CRY2, PHYA), and flowering (AGL20, AP2, FLC). Further, results from this and previous studies highlight homologs most similar to Arabidopsis HY5, MAF3, RVE1 and RD22 as potential molecular markers for endodormancy in crown buds of leafy spurge. Early response to mild dehydration also highlighted involvement of upstream ethylene and JA-signaling, whereas severe dehydration impacted ABA-signaling. The identification of conserved ABRE- and MYC-consensus, cis-acting elements in the promoter of leafy spurge genomic clones similar to Arabidopsis RVE1 (AT5G17300) implicates a potential role for ABA-signaling in its dehydration-induced expression. Response of these molecular mechanisms to dehydration-stress provides insights on the ability of invasive perennial weeds to adapt and survive under harsh environments, which will be beneficial for addressing future management practices.

[Identification of medicinal plants used as Tibetan traditional medicine Chuan-Bu based on nrDNA markers]

OBJECTIVE

To identify the common Tibetan herb Chuan-Bu.

METHOD

Local herbalists were visited to observe which plants were being used as Chuan-Bu. Samples of the indigenous plants were collected at the same time. Leaf materials were collected from field surveys. Total genomic DNA was extracted from silica gel-dried leaf samples. The PCR products were purified and directly sequenced.

RESULT

As the origin of Chuan-Bu in Tibet autonomous region was authenticated, two species were determined, i. e. Euphorbia stracheyiand E. wallichii. Also, based on our earlier research, the origin of Chuan-Bu in Gansu province, is from E. kansuensis. The sequences of ITS1 for E. stracheyi and E. wallichii were 261 bp in size, and 221 bp in ITS2, respectively. The size of the 5.8S coding region was 164 bp for all species examined in the genus. Especially, there was a heterozygous locus in ITS1 (C/G; position 72) for E. stracheyi. The nucleotide divergence between sequences of the 6 species in pairwise comparisons was calculated and the result showed that the variable site could be detected in each pairwise comparison of sequences. Also, there were 8 point mutations in the 5.8S coding region.

CONCLUSION

nrDNA ITS sequences can be used as the molecular markers to identify the Tibetan herb Chuan-Bu and such Traditional Chinese Medicines from the same genus Euphorbia as E. lathyris, E. humifusa and E. pekinensis.

Stem fasciation in cacti and succulent species--tissue anatomy, protein pattern and RAPD polymorphisms

Fasciated and normal stem segments of Opuntia microdasys, Opuntia cylindrica, Huernia primulina and Euphorbia lactea were collected from the same plant and compared for their anatomy, water relations and genetic variations. Anatomical differences in terms of thickness of cuticle, vascular bundle, xylem and phloem were analyzed in both normal and fasciated stems. The mucilage cells were higher in the fasciated form of Opuntia microdasys than that in the normal form. Water status in terms of total water content (TWC), water deficit and relative water content (RWC) was influenced by fasciation. Genetic variations were tested in normal and fasciated stems using randomly amplified polymorphic DNA (RAPD) fingerprints and SDS-PAGE of soluble protein extracts. SDS-PAGE protein and RAPD analysis confirmed that normal and fasciated tissues were genetically different. Polymerase chain reaction (PCR) yielded different polymorphic banding patterns that were unique to each primer and distinguishable over all samples. The PCR results of normal and fasciated samples were significantly different in cases of primers P1, P2 and P3. These results indicate that occurrence of fasciation in Opuntia microdasys, Opuntia cylindrica, Huernia primulina and Euphorbia lactea is an epigenetic mutation of tissues.

Euphorbia tirucalli L.-comprehensive characterization of a drought tolerant plant with a potential as biofuel source

Of late, decrease in mineral oil supplies has stimulated research on use of biomass as an alternative energy source. Climate change has brought problems such as increased drought and erratic rains. This, together with a rise in land degeneration problems with concomitant loss in soil fertility has inspired the scientific world to look for alternative bio-energy species. Euphorbia tirucalli L., a tree with C3/CAM metabolism in leaves/stem, can be cultivated on marginal, arid land and could be a good alternative source of biofuel. We analyzed a broad variety of E. tirucalli plants collected from different countries for their genetic diversity using AFLP. Physiological responses to induced drought stress were determined in a number of genotypes by monitoring growth parameters and influence on photosynthesis. For future breeding of economically interesting genotypes, rubber content and biogas production were quantified. Cluster analysis shows that the studied genotypes are divided into two groups, African and mostly non-African genotypes. Different genotypes respond significantly different to various levels of water. Malate measurement indicates that there is induction of CAM in leaves following drought stress. Rubber content varies strongly between genotypes. An investigation of the biogas production capacities of six E. tirucalli genotypes reveals biogas yields higher than from rapeseed but lower than maize silage.

The transcriptomes of dormant leafy spurge seeds under alternating temperature are differentially affected by a germination-enhancing pretreatment

Seed dormancy is an important stage in the life cycle of many non-domesticated plants, often characterized by the temporary failure to germinate under conditions that normally favor the process. Pre-treating dormant imbibed seeds at a constant temperate accelerated germination of leafy spurge seeds under alternating temperatures. However, dormant seeds will also germinate without a pre-treatment, albeit at a much slower rate, which gives rise to longer periods of imbibition before germination. Transcriptome analyses on seeds exposed to prolonged imbibition highlighted pathways associated with phenylpropanoid biosynthesis and interacting networks of genes involved in plant defense. In addition to the many pathways associated with phenylpropanoid biosynthesis enriched with down-regulated genes upon germination, there were also numerous pathways enriched with up-regulated genes associated with energy metabolism, such as glycolysis. Transcriptome data further suggest that metabolism and signaling by the plant hormones ethylene, gibberellin, and abscisic acid are involved in the developmental transition from dormancy to germination. More specifically, sub-network enrichment analysis identified ABI3 as a central hub of a sub-network at germination including several down-regulated genes such as DELLA (i.e., RGL2), which represses gibberellin signaling processes required for germination. The 595-fold increase in the expression of ACC oxidase (ACO4) at germination also suggests an important role for ethylene biosynthesis in germinating leafy surge seeds. Furthermore, the 10-578-fold difference in expression of many genes such as HY5 and Histone H3 between two populations at germination, which were treated with and without a constant temperature germination-enhancing pretreatment, revealed disparate impacts on various biosynthetic, growth, signaling, and response processes. Overall, our results indicate a constant temperature pretreatment (20°C for 21d) is not required for germination of leafy spurge seeds at an alternating temperature. However, the presence or absence of the pretreatment does affect the rate of germination and the germination transcriptional programs.

Induction of endodormancy in crown buds of leafy spurge (Euphorbia esula L.) implicates a role for ethylene and cross-talk between photoperiod and temperature

Leafy spurge is a model for studying well-defined phases of dormancy in underground adventitious buds (UABs) of herbaceous perennial weeds, which is a primary factor facilitating their escape from conventional control measures. A 12-week ramp down in both temperature (27 → 10 °C) and photoperiod (16 → 8 h light) is required to induce a transition from para- to endo-dormancy in UABs of leafy spurge. To evaluate the effects of photoperiod and temperature on molecular networks of UABs during this transition, we compared global transcriptome data-sets obtained from leafy spurge exposed to a ramp down in both temperature and photoperiod (RDtp) versus a ramp down in temperature (RDt) alone. Analysis of data-sets indicated that transcript abundance for genes associated with circadian clock, photoperiodism, flowering, and hormone responses (CCA1, COP1, HY5, MAF3, MAX2) preferentially increased in endodormant UABs. Gene-set enrichment analyses also highlighted metabolic pathways involved in endodormancy induction that were associated with ethylene, auxin, flavonoids, and carbohydrate metabolism; whereas, sub-network enrichment analyses identified hubs (CCA1, CO, FRI, miR172A, EINs, DREBs) of molecular networks associated with carbohydrate metabolism, circadian clock, flowering, and stress and hormone responses. These results helped refine existing models for the transition to endodormancy in UABs of leafy spurge, which strengthened the roles of circadian clock associated genes, DREBs, COP1-HY5, carbohydrate metabolism, and involvement of hormones (ABA, ethylene, and strigolactones). We further examined the effects of ethylene by application of 1-aminocyclopropane-1-carboxylate (ACC) to paradormant plants without a ramp down treatment. New vegetative growth from UABs of ACC-treated plants resulted in a dwarfed phenotype that mimicked the growth response in RDtp-induced endodormant UABs. The results of this study provide new insights into dormancy regulation suggesting a short-photoperiod treatment provides an additive cross-talk effect with temperature signals that may impact ethylene's effect on AP2/ERF family members.

Selection and validation of endogenous reference genes for qRT-PCR analysis in leafy spurge (Euphorbia esula)

Quantitative real-time polymerase chain reaction (qRT-PCR) is the most important tool in measuring levels of gene expression due to its accuracy, specificity, and sensitivity. However, the accuracy of qRT-PCR analysis strongly depends on transcript normalization using stably expressed reference genes. The aim of this study was to find internal reference genes for qRT-PCR analysis in various experimental conditions for seed, adventitious underground bud, and other organs of leafy spurge. Eleven candidate reference genes (BAM4, PU1, TRP-like, FRO1, ORE9, BAM1, SEU, ARF2, KAPP, ZTL, and MPK4) were selected from among 171 genes based on expression stabilities during seed germination and bud growth. The other ten candidate reference genes were selected from three different sources: (1) 3 stably expressed leafy spurge genes (60S, bZIP21, and MD-100) identified from the analyses of leafy spurge microarray data; (2) 3 orthologs of Arabidopsis “general purpose” traditional reference genes (GAPDH_1, GAPDH_2, and UBC); and (3) 4 orthologs of Arabidopsis stably expressed genes (UBC9, SAND, PTB, and F-box) identified from Affymetrix ATH1 whole-genome GeneChip studies. The expression stabilities of these 21 genes were ranked based on the C_T values of 72 samples using four different computation programs including geNorm, Normfinder, BestKeeper, and the comparative ΔC_T method. Our analyses revealed SAND, PTB, ORE9, and ARF2 to be the most appropriate reference genes for accurate normalization of gene expression data. Since SAND and PTB were obtained from 4 orthologs of Arabidopsis, while ORE9 and ARF2 were selected from 171 leafy spurge genes, it was more efficient to identify good reference genes from the orthologs of other plant species that were known to be stably expressed than that of randomly testing endogenous genes. Nevertheless, the two newly identified leafy spurge genes, ORE9 and ARF2, can serve as orthologous candidates in the search for reference genes from other plant species.

Poinsettia's wild ancestor in the Mexican dry tropics: Historical, genetic, and environmental evidence

PREMISE OF THE STUDY

The poinsettia (Euphorbia pulcherrima) is the world's most economically important potted plant, but despite its preeminence it is not clear which wild populations are ancestral to the varieties cultivated around the world. Tradition holds that the U.S. envoy to Mexico J. R. Poinsett collected the progenitors of the over 300 varieties in global cultivation on an 1828 excursion to northern Guerrero State, Mexico. It is unknown whether the contemporary cultivars are descended from plants from Guerrero or whether germplasm from other parts of poinsettia's 2000 km long distribution entered into cultivation during the nearly 200 yr of subsequent poinsettia horticulture.

METHODS

To identify the wild populations that likely gave rise to the cultivars and test this historical account, we sequenced plastid and nuclear DNA regions and modeled poinsettia's potential distribution.

KEY RESULTS

The combination of nuclear and plastid haplotypes characterizing cultivars was found only in northern Guerrero. Distribution modeling indicated that suitable habitat conditions for wild poinsettias are present in this area, consistent with their likely wild status.

CONCLUSIONS

Our data pinpoint the area of northern Guerrero as the cultivated poinsettia's probable ancestral region, congruent with the traditional account attributing the original collections to Poinsett. Abundant genetic variation likely offers raw material for improving the many shortcomings of cultivars, including vulnerability to cold, stem breakage, and pathogens such as Pythium and Phytophthora. However, genetic differences between populations make conservation of all of poinsettia's diversity difficult.

Phylogenetics of the Chamaesyce clade (Euphorbia, Euphorbiaceae): reticulate evolution and long-distance dispersal in a prominent C4 lineage

PREMISE OF THE STUDY

The Chamaesyce clade of Euphorbia is the largest lineage of C(4) plants among the eudicots, with 350 species including both narrow endemics and cosmopolitan weeds. We sampled this group worldwide to address questions about subclade relationships, the origin of C(4) photosynthesis, the evolution of weeds, and the role of hybridization and long-distance dispersal in the diversification of the group. •

METHODS

Two nuclear (ITS and exon 9 of EMB2765) and three chloroplast markers (matK, rpl16, and trnL-F) were sequenced for 138 ingroup and six outgroup species. Exon 9 of EMB2765 was cloned in accessions with >1% superimposed peaks. •

KEY RESULTS

The Chamaesyce clade is monophyletic and consists of three major subclades [1(2,3)]: (1) the Acuta clade, containing three North American species with C(3) photosynthesis and C(3)-C(4) intermediates; (2) the Peplis clade, mostly North American and entirely C(4); and (3) the Hypericifolia clade, all C(4), with both New World and Old World groups. Incongruence between chloroplast and ITS phylogenies and divergent cloned copies of EMB2765 exon 9 suggest extensive hybridization, especially in the Hawaiian Islands radiation. •

CONCLUSIONS

The Chamaesyce clade originated in warm, arid areas of North America, where it evolved C(4) photosynthesis. From there, it diversified globally with extensive reticulate evolution and frequent long-distance dispersals. Although many species are weedy, there are numerous local adaptations to specific substrates and regional or island radiations, which have contributed to the great diversity of this group.

Linking self-incompatibility, dichogamy, and flowering synchrony in two Euphorbia species: alternative mechanisms for avoiding self-fertilization?

BACKGROUND

Plant species have several mechanisms to avoid selfing such as dichogamy or a self-incompatibility response. Dichogamy in a single flower may reduce autogamy but, to avoid geitonogamy, plants must show flowering synchronization among all their flowers (i.e. synchronous dichogamy). It is hypothesized that one species would not simultaneously show synchronous dichogamy and self-incompatibility because they are redundant mechanisms to reduce selfing; however, this has not been accurately assessed.

METHODOLOGY/PRINCIPAL FINDINGS

This expectation was tested over two years in two natural populations of the closely related Mediterranean spurges Euphorbia boetica and E. nicaeensis, which completely avoid autogamy by protogyny at the cyathia level. Both spurges showed a high population synchrony (Z<79), and their inflorescences flower synchronously. In E. nicaeensis, there was no overlap among the cyathia in anthesis of successive inflorescence levels and the overlap between sexual phases of cyathia of the same inflorescence level was uncommon (4-16%). In contrast, E. boetica showed a high overlap among consecutive inflorescence levels (74-93%) and between sexual phases of cyathia of the same inflorescence level (48-80%). The flowering pattern of both spurges was consistent in the two populations and over the two successive years. A hand-pollination experiment demonstrated that E. nicaeensis was strictly self-compatible whereas E. boetica was partially self-incompatible.

CONCLUSIONS/SIGNIFICANCE

We propose that the complex pattern of synchronized protogyny in E. nicaeensis prevents geitonogamous crosses and, consequently, avoids selfing and inbreeding depression. In E. boetica, a high probability of geitonogamous crosses may occur but, alternatively, this plant escapes selfing through a self-incompatibility response. We posit that synchronous dichogamy and physiological self-incompatibility do not co-occur in the same species because each process is sufficiently effective in avoiding self-fertilization.

Autophosphorylation is crucial for CDK-activating kinase (Ee;CDKF;1) activity and complex formation in leafy spurge

Ee;CDKF;1 protein is a leafy spurge (Euphorbia esula) CDK-activating kinase that is involved in a phosphorylation cascade linked to early stages of cell cycle progression. Yeast two-hybrid screening performed using Ee;CDKF;1 as a bait indicated that one of the interacting proteins was Ee;CDKF;1. Protein-protein interaction of Ee;CDKF;1 was further confirmed by yeast two-hybrid interaction and in vitro pull-down assays. Gel exclusion chromatography and/or native PAGE showed that GST-CDKF;1, MBP-CDKF;1, GST-CDKF;1 devoid of GST, and endogenous Ee;CDKF;1 were capable of forming homo protein complexes which are in dimer, trimer, and/or higher molecular-mass complex in its native state. In addition, Ee;CDKF;1 complexes were autophosphorylated and able to phosphorylate CDK. Moreover, mutant forms of Ee;CDKF;1 (106G/A, 166K/A), which lost autophosphorylation capability completely, were unable to form homo protein complexes in their native state. The result thus demonstrated that autophosphorylation of Ee;CDKF;1 is crucial for both kinase activity and complex formation.

Is LEAFY a useful marker gene for the flower-inflorescence boundary in the Euphorbia cyathium?

The flower-like reproductive structure of Euphorbia s.l. (Euphorbiaceae) is widely believed to have evolved from an inflorescence, and is therefore interpreted as a special type of pseudanthium, termed a cyathium. However, fuzzy morphological boundaries between the inflorescence, individual flowers, and organs have fuelled the suggestion that the cyathium does not merely superficially resemble a flower but could actually share developmental genetic pathways with a conventional flower. To test this hypothesis, immunolocalizations of FLORICAULA/LEAFY (LFY), a protein associated with floral identity in many angiosperm species, were performed in developing cyathia of different species of Euphorbia. Expression of the LFY protein was found not only in individual floral primordia (as predicted from results in the model organisms Arabidopsis and Anthirrhinum), but also in the cyathium primordium and in the primordia of partial male inflorescences. These results provide further evidence that the evolution of floral traits in pseudanthial inflorescences often involves expression of floral development genes in the inflorescence apex. This finding blurs the conventional rigid distinction between flowers and inflorescences.

A reliable and efficient method for total rna isolation from various members of spurge family (Euphorbiaceae)

INTRODUCTION

It is prerequisite and crucial to extract RNA with high quality and integrity in order to carry out molecular biology studies in any plant species of a family. Euphorbiaceae members are known for high levels of their waxes, oils with polysaccharides, polyphenolics and secondary metabolites. These conditions are recognised to interfere unfavourably with various methodologies of RNA isolation.

OBJECTIVE

To develop a simple, rapid and reproducible cetyltrimethylamonium bromide (CTAB)-based protocol, to reduce the time and cost of extraction without reducing quality and yield of RNA extracted from various recalcitrant Euphorbiaceae member plant tissues such as from tree leaves (Hevea brasilensis), woody shrubs leaves (Ricinus communis, Jatropha curcas, Manihot esculenta) and storage root tissue (M. esculenta).

METHODOLOGY

Simple modifications and fast steps were introduced to the original CTAB protocol. All centrifugation steps were carried out at 4°C at 12000 rpm for 10 min, the sample weight was decreased and usage of spermidine and LiCl was omitted, reducing incubation time prior to RNA precipitation. This rapid CTAB protocol was compared with various RNA isolation methods intended for use with plants rich in polysaccharides and secondary metabolites.

RESULTS

The procedure can be completed within 2 h and many samples can be processed at the same time. RNA of high quality could be isolated from all the tissues of species that we tried. The isolated RNA from different species served as a robust template for RT-PCR analysis.

CONCLUSION

The study has shown that the improvement of a CTAB-based protocol allows the rapid isolation of high-quality RNA from various recalcitrant Euphorbiaceae members.

Initial changes in the transcriptome of Euphorbia esula seeds induced to germinate with a combination of constant and diurnal alternating temperatures

We investigated transcriptome changes in Euphorbia esula (leafy spurge) seeds with a focus on the effect of constant and diurnal fluctuating temperature on dormancy and germination. Leafy spurge seeds do not germinate when incubated for 21 days at 20 degrees C constant temperatures, but nearly 30% germinate after 21 days under fluctuating temperatures 20:30 degrees C (16:8 h). Incubation at 20 degrees C for 21 days followed by 20:30 degrees C resulted in approximately 63% germination in about 10 days. A cDNA microarray representing approximately 22,000 unique sequences was used to profile transcriptome changes in the first day after transfer of seeds from constant to alternating temperature conditions. Functional classification based on MIPS and gene ontology revealed active metabolism including up-regulation of energy, protein synthesis, and signal transduction processes. Down-regulated processes included translation elongation, translation, and some biosynthetic processes. Subnetwork analysis identified genes involved in abscisic acid, sugar, and circadian clock signaling as key regulators of physiological activity in seeds soon after the transfer to alternating conditions.

Characterization, expression and function of DORMANCY ASSOCIATED MADS-BOX genes from leafy spurge

DORMANCY ASSOCIATED MADS-BOX (DAM) genes are related to AGAMOUS-LIKE 24 and SHORT VEGETATIVE PHASE genes of arabidopsis and are differentially regulated coordinately with endodormancy induction and release in buds of several perennial plant species. DAM genes were first shown to directly impact endodormancy in peach where a deletion of a series of DAM resulted in loss of endodormancy induction. We have cloned and characterized several MADS box genes from the model perennial weed leafy spurge. Leafy spurge DAM genes are preferentially expressed in shoot tips and buds in response to cold temperatures and day length in a manner that is relative to the level of endodormancy induced by various environmental conditions. Over-expression of one DAM gene in arabidopsis delays flowering. Additionally, we show that at least one DAM gene is differentially regulated by chromatin remodeling. Comparisons of the DAM gene promoters between poplar and leafy spurge have identified several conserved sequences that may be important for their expression patterns in response to dormancy-inducing stimuli.

Low temperatures impact dormancy status, flowering competence, and transcript profiles in crown buds of leafy spurge

Leafy spurge (Euphorbia esula) is an herbaceous perennial weed that produces vegetatively from an abundance of underground adventitious buds. In this study, we report the effects of different environmental conditions on vegetative production and flowering competence, and determine molecular mechanisms associated with dormancy transitions under controlled conditions. Reduction in temperature (27-10 degrees C) and photoperiod (16-8 h) over a 3-month period induced a para- to endo-dormant transition in crown buds. An additional 11 weeks of cold (5-7 degrees C) and short-photoperiod resulted in accelerated shoot growth from crown buds, and 99% floral competence when plants were returned to growth-promoting conditions. Exposure of paradormant plants to short-photoperiod and prolonged cold treatment alone had minimal affect on growth potential and resulted in ~1% flowering. Likewise, endodormant crown buds without prolonged cold treatment displayed delayed shoot growth and ~2% flowering when returned to growth-promoting conditions. Transcriptome analysis revealed that 373 and 260 genes were differentially expressed (P < 0.005) during para- to endo-dormant and endo- to eco-dormant transitions, respectively. Transcripts from flower competent vs. non-flower competent crown buds identified 607 differentially expressed genes. Further, sub-network analysis identified expression targets and binding partners associated with circadian clock, dehydration/cold signaling, phosphorylation cascades, and response to abscisic acid, ethylene, gibberellic acid, and jasmonic acid, suggesting these central regulators affect well-defined phases of dormancy and flowering. Potential genetic pathways associated with these dormancy transitions and flowering were used to develop a proposed conceptual model.

Changes in the expression of carbohydrate metabolism genes during three phases of bud dormancy in leafy spurge

Underground adventitious buds of leafy spurge (Euphorbia esula) undergo three well-defined phases of dormancy, para-, endo-, and ecodormancy. In this study, relationships among genes involved in carbohydrate metabolism and bud dormancy were examined after paradormancy release (growth induction) by decapitation and in response to seasonal signals. Real-time PCR was used to determine the expression levels of carbohydrate metabolism genes at different phases of bud dormancy. Among differentially-regulated genes, expression of a specific Euphorbia esula beta-amylase gene (Ee-BAM1) increased 100-fold after growth induction and 16,000-fold from July (paradormancy) to December (ecodormancy). Sequence data analysis indicated that two genes, Ee-BAM1 and Ee-BAM2, could encode this beta-amylase. However, real-time PCR using gene-specific primer pairs only amplified Ee-BAM1, indicating that Ee-BAM2 is either specific to other organs or not abundant. The deduced amino acid sequences of these two genes are very similar at the N-terminal but differ at the C-terminal. Both contain a nearly identical, predicted 48-amino acid plastid transit peptide. Immunoblot analyses identified a 29 kD (mature Ee-BAM1 after cleavage of the transit peptide) and a 35 kD (unprocessed EeBAM1) protein. Both 35 and 29 kD proteins were constitutively expressed in growth-induced and seasonal samples. Immunolocalization indicated that Ee-BAM1 is in the cytosol of cells at the shoot tip of the bud. Ee-BAM1 also surrounds the amyloplasts in mature cells toward the base of the bud. These observations suggests that Ee-BAM1 may have dual functions; serving as reserve protein in the cytosol and as a degrading enzyme at the surface of amyloplasts.

Cloning and characterization of a squalene synthase gene from a petroleum plant, Euphorbia tirucalli L

Euphorbia tirucalli L., which is also known as a petroleum plant, produces a large amount of phytosterols and triterpenes. During their biosynthesis, squalene synthase converts two molecules of the hydrophilic substrate farnesyl diphosphate into a hydrophobic product, squalene. An E. tirucalli cDNA clone of a putative squalene synthase gene (EtSS) was isolated by RT-PCR followed by 5'- and 3'-RACE. The restriction fragment polymorphisms revealed by Southern blot analysis suggest that EtSS is a single copy gene. The glycine at the 287th residue from the N-terminal end of domain C has replaced alanine, which is conserved among all the other SS sequences deposited in the Genbank database. The N-terminal 380 residues of the hydrophilic sequence was expressed as a peptide-tagged protein in E. coli, and the resultant bacterial crude extract was incubated with farnesyl diphosphate and NADPH. GC-MS analysis showed that squalene was detected in the in vitro reaction mixture. E. tirucalli transgenic callus lines, in which EtSS was overexpressed, accumulated increased amounts of phytosterols as compared with that of wild type callus. RT-PCR analysis of wild type E. tirucalli plants revealed that the EtSS transcript accumulated in almost equal amounts in the stems and the leaves with a stalk, while a lower amount was detected in the roots. In situ hybridization analysis revealed that prominent antisense-probe signal was detected in the cambia within bundle sheathes. These results indicate that EtSS functions prominently in cambia, which are located adjacent to conductive tubes, and that this gene plays important roles in phytosterol accumulation in petroleum plants.

Agrobacterium tumefaciens-mediated transformation of poinsettia, Euphorbia pulcherrima, with virus-derived hairpin RNA constructs confers resistance to Poinsettia mosaic virus

Agrobacterium-mediated transformation for poinsettia (Euphorbia pulcherrima Willd. Ex Klotzsch) is reported here for the first time. Internode stem explants of poinsettia cv. Millenium were transformed by Agrobacterium tumefaciens, strain LBA 4404, harbouring virus-derived hairpin (hp) RNA gene constructs to induce RNA silencing-mediated resistance to Poinsettia mosaic virus (PnMV). Prior to transformation, an efficient somatic embryogenesis system was developed for poinsettia cv. Millenium in which about 75% of the explants produced somatic embryos. In 5 experiments utilizing 868 explants, 18 independent transgenic lines were generated. An average transformation frequency of 2.1% (range 1.2-3.5%) was revealed. Stable integration of transgenes into the poinsettia nuclear genome was confirmed by PCR and Southern blot analysis. Both single- and multiple-copy transgene integration into the poinsettia genome were found among transformants. Transgenic poinsettia plants showing resistance to mechanical inoculation of PnMV were detected by double antibody sandwich enzyme-linked immunosorbent assay (DAS-ELISA). Northern blot analysis of low molecular weight RNA revealed that transgene-derived small interfering (si) RNA molecules were detected among the poinsettia transformants prior to inoculation. The Agrobacterium-mediated transformation methodology developed in the current study should facilitate improvement of this ornamental plant with enhanced disease resistance, quality improvement and desirable colour alteration. Because poinsettia is a non-food, non-feed plant and is not propagated through sexual reproduction, this is likely to be more acceptable even in areas where genetically modified crops are currently not cultivated.

Anthocyanin inhibits propidium iodide DNA fluorescence in Euphorbia pulcherrima: implications for genome size variation and flow cytometry

BACKGROUND

Measuring genome size by flow cytometry assumes direct proportionality between nuclear DNA staining and DNA amount. By 1997 it was recognized that secondary metabolites may affect DNA staining, thereby causing inaccuracy. Here experiments are reported with poinsettia (Euphorbia pulcherrima) with green leaves and red bracts rich in phenolics.

METHODS

DNA content was estimated as fluorescence of propidium iodide (PI)-stained nuclei of poinsettia and/or pea (Pisum sativum) using flow cytometry. Tissue was chopped, or two tissues co-chopped, in Galbraith buffer alone or with six concentrations of cyanidin-3-rutinoside (a cyanidin-3-rhamnoglucoside contributing to red coloration in poinsettia).

KEY RESULTS

There were large differences in PI staining (35-70 %) between 2C nuclei from green leaf and red bract tissue in poinsettia. These largely disappeared when pea leaflets were co-chopped with poinsettia tissue as an internal standard. However, smaller (2.8-6.9 %) differences remained, and red bracts gave significantly lower 1C genome size estimates (1.69-1.76 pg) than green leaves (1.81 pg). Chopping pea or poinsettia tissue in buffer with 0-200 microm cyanidin-3-rutinoside showed that the effects of natural inhibitors in red bracts of poinsettia on PI staining were largely reproduced in a dose-dependent way by this anthocyanin.

CONCLUSIONS

Given their near-ubiquitous distribution, many suspected roles and known affects on DNA staining, anthocyanins are a potent, potential cause of significant error variation in genome size estimations for many plant tissues and taxa. This has important implications of wide practical and theoretical significance. When choosing genome size calibration standards it seems prudent to select materials producing little or no anthocyanin. Reviewing the literature identifies clear examples in which claims of intraspecific variation in genome size are probably artefacts caused by natural variation in anthocyanin levels or correlated with environmental factors known to induce variation in pigmentation.

Differentiation of the traditional Chinese medicinal plants Euphorbia humifusa and E. maculata from adulterants by TaqMan real-time polymerase chain reaction

DNA sequence analysis of the rDNA internal transcribed spacer 1 (ITS1) and TaqMan real-time polymerase chain reaction were exploited for their applications in the differentiation of the traditional chinese medicinal plants euphorbia humifusa and e. maculata from three related adulterants e. hypericifolia, E. atoto and E. prostrata. The data demonstrated that variations in the ITS1 regions were very low at the intra-species level but extremely high at the inter-species level, so that they could be easily distinguished at the DNA level. The sequence difference allowed an effective and reliable differentiation of E. humifusa and E. maculata from the adulterants by TaqMan real-time PCR.

Expression of the gene for sterol-biosynthesis enzyme squalene epoxidase in parenchyma cells of the oil plant, Euphorbia tirucalli

In plants, phytosterols and triterpenes are major secondary metabolites. In an attempt to reveal the mechanism for synthesis and storage of these compounds, we isolated and characterized cDNA clones for squalene epoxidase (SE), from a succulent shrub, Euphorbia tirucalli. Southern-blot analysis of total DNA using cDNA fragment as a probe showed that the E. tirucalli squalene epoxidase gene (EtSE) is single-copy type in terms of restriction fragment length polymorphism (RFLP). Deduced amino-acid sequence of the cDNA showed 83 and 75% identity to those of rice and ginseng, respectively, in an area excluding a less homologous putative transmembrane region in the N-terminal end. Functional characterization with heterologous expression using an erg1-disrupted yeast mutant KLN1 indicated that the EtSE recovered ergosterol auxotrophy of the mutant, and gave rise to an ergosterol accumulation in the EtSE transformant. RT-PCR analysis showed the EtSE transcripts in leaves and stem internodes accumulated in almost equal amounts, which were more abundant than those in roots. In situ hybridization using EtSE antisense probe revealed prominent EtSE expression on a parenchyma cell adjacent to primary laticifers that were located in a rosary orientation in the inner region of cortex. This is the first report of expression of a gene for a rate-limiting enzyme in mevalonate pathway in organs and tissues of a plant.

Catalase and antiquitin from Euphorbia characias: Two proteins involved in plant defense?

Here we report the cDNA nucleotide sequences of a calmodulin-binding catalase and an antiquitin from the latex of the Mediterranean shrub Euphorbia characias. Present findings suggest that catalase and antiquitin might represent additional nodes in the Euphorbia defense systems, and a multi-enzymatic interaction contributing to plant's protection against biotic and abiotic stresses is proposed to occur in E. characias laticifers.

Potential roles for autophosphorylation, kinase activity, and abundance of a CDK-activating kinase (Ee;CDKF;1) during growth in leafy spurge

Leafy spurge (Euphorbia esula L.) is a deep-rooted perennial weed that propagates both by seeds and underground adventitious buds located on the crown and roots. To enhance our understanding of growth and development during seed germination and vegetative propagation, a leafy spurge gene (Accession No. AF230740) encoding a CDK-activating kinase (Ee;CDKF;1) involved in cell-cycle progression was identified, and its function was confirmed based on its ability to rescue a yeast temperature-sensitive CAK mutant (GF2351) and through in vitro kinase assays. Site-directed mutagenesis of Ee;CDKF;1 indicated that two threonine residues (Thr291 and Thr296) were mutually responsible for intra-molecular autophosphorylation and for phosphorylating its substrate protein, cyclin-dependent kinase (CDK). Polyclonal antibodies generated against the Ee;CDKF;1 protein or against a phosphorylated Ee;CDKF;1 peptide [NERYGSL(pT)SC] were used to examine abundance and phosphorylation of CDKF;1 during seed germination and bud growth. The levels of CDKF;1 were lower in dry or imbibed seeds than in germinating seeds or seedlings. Differences in CDKF;1 were also observed during adventitious bud development; small buds appeared to have greater levels of CDKF;1 than large buds. Similar patterns of CDKF;1 expression were detected with either the polyclonal antibody developed using the CDKF;1 protein or the phosphorylated peptide. These results indicated that Thr291 is constitutively phosphorylated in vivo and associated with Ee;CDKF;1 activity. Our results further suggest that a certain level of CDKF;1 activity is maintained in most tissues and may be an important phenomenon for enzymes that regulate early steps in cell-cycle signaling pathways.

Molecular authentication of the traditional Chinese medicinal plant Euphorbia pekinensis

The ITS regions of Euphorbia pekinensis and six other Euphorbia species used as adulterants of E. pekinensis were sequenced to differentiate them. The sequences are identical among the individuals in the seven species studies. Diversity in DNA sequences among various species was found ranging from 8.3% to 43.8% in ITS1 and 7.6% to 36.6% in ITS2 region. Furthermore, based on the divergent ITS regions, species-specific primers, JDJp 1 and JDJp 2, were designed in the polymorphic regions of E. pekinensis to distinguish it from adulterants. These ITS-derived primers amplified a 281-bp-specific DNA fragment from E. pekinensis. No amplified product was observed using DNA of six adulterants.

Characterization of SCP-2 from Euphorbia lagascae reveals that a single Leu/Met exchange enhances sterol transfer activity

Sterol carrier protein-2 (SCP-2) is a small intracellular basic protein domain implicated in peroxisomal beta-oxidation. We extend our knowledge of plant SCP-2 by characterizing SCP-2 from Euphorbia lagascae. This protein consists of 122 amino acids including a PTS1 peroxisomal targeting signal. It has a molecular mass of 13.6 kDa and a pI of 9.5. It shares 67% identity and 84% similarity with SCP-2 from Arabidopsis thaliana. Proteomic analysis revealed that E. lagascae SCP-2 accumulates in the endosperm during seed germination. It showed in vitro transfer activity of BODIPY-phosphatidylcholine (BODIPY-PC). The transfer of BODIPY-PC was almost completely inhibited after addition of phosphatidylinositol, palmitic acid, stearoyl-CoA and vernolic acid, whereas sterols only had a very marginal inhibitory effect. We used protein modelling and site-directed mutagenesis to investigate why the BODIPY-PC transfer mediated by E. lagascae SCP-2 is not sensitive to sterols, whereas the transfer mediated by A. thaliana SCP-2 shows sterol sensitivity. Protein modelling suggested that the ligand-binding cavity of A. thaliana SCP-2 has four methionines (Met12, 14, 15 and 100), which are replaced by leucines (Leu11, 13, 14 and 99) in E. lagascae SCP-2. Changing Leu99 to Met99 was sufficient to convert E. lagascae SCP-2 into a sterol-sensitive BODIPY-PC-transfer protein, and correspondingly, changing Met100 to Leu100 abolished the sterol sensitivity of A. thaliana SCP-2.

Subtractive cDNA libraries identify differentially expressed genes in dormant and growing buds of leafy spurge (Euphorbia esula)

Two subtractive cDNA libraries were developed to study genes associated with bud dormancy (reverse library) and initiation of shoot growth (forward library) in leafy spurge. To identify unique sequences represented in each library, 15744 clones were screened to reduce the level of redundancy within both libraries. A total of 516 unique sequences were obtained from 2304 minimally redundant clones. Radioactive probes developed from RNAs extracted from crown buds of either intact (para-dormant control) or a series of growth-induced (2 h, 2, and 4 d after decapitation) plants were used to identify differentially expressed genes by macroarray analysis. Semi-quantitative RT-PCR was used to confirm results obtained by macroarray analysis and to determine the expression profiles for other transcripts identified within the subtractive libraries. Selected clones were also used to examine gene expression in crown buds after growth induction and/or during normal seasonal growth. In this study, four distinct patterns of gene expression were observed during the transition from para-dormancy to growth-induction. Many of the differentially regulated genes identified have unknown or hypothetical functions while others are known to play important roles in molecular functions. Gene ontology analysis identified a greater proportion of genes involved with catalytic activity in the forward library while the reverse library had a greater proportion of genes involved in DNA/RNA binding.

Cloning and characterization of a cDNA encoding beta-amyrin synthase from petroleum plant Euphorbia tirucalli L

Euphorbia tirucalli L., known as the petroleum plant, produces a large amount of triterpenes, such as beta-amyrin. Degenerate RT-PCR based on the sequences conserved among known beta-amyrin synthases led to cloning of a putative triterpene synthase cDNA, EtAS, from leaves of E. tirucalli. The deduced amino acid sequence of the EtAS cDNA showed the highest identity of 82% to the Panax ginseng beta-amyrin synthase. Heterologous expression of the EtAS ORF in the methylotrophic yeast, Pichia pastoris, resulted in production of beta-amyrin, revealing that the EtAS cDNA codes for a beta-amyrin synthase. This is the first report of a gene involved in the triterpene synthetic pathway from Euphorbiaceae plants.

Identification and expression analysis of the cyclophilin gene in Kandelia candel under stress of salt

Two cDNA fragments, named for SRGKC2 and SRGKC3, encoding cyclophilin in Kandelia candel were isolated by Representational Difference Analysis of cDNA. The two cDNA fragments were 282 bp and 160 bp, respectively. Sequence analysis shows that both of the SRGKC2 and SRGKC3 come from the same gene region, and SRGKC3 is a part of SRGKC2. In addition the SRGKC2 displayed 90% sequence identity over a region of 84 amino acids to the cyclophilin from Euphorbia esula and the SRGKC3 displayed 93% sequence identity over a region of 47 amino acids to the fava bean. The Northern blotting showed that the expression of SRGKC2 was suppressed under stress of salt. Based on the sequence of SRGKC2, a full-length cDNA (KCCYP1) was isolated by RACE reaction (This sequence data has been submitted to the EMBL databases under accession No. AY150052). The full-length cDNA was about 0.9 kb, which contained an open reading frame (ORF) of 516 bp and coded for 172 amino acid residues with isoelectric point of 8.57 and molecular weight of 18.2 kD. The motif A of the ATP/GTP-binding site in KCCYP1 appears at amino acid residues of 41-49, and seven-amino-acids-residue was inserted at 48-54 amino acid residues. The expression patterns of SRGKC2 in various species were also investigated.