Characterization, comparison, and phylogenetic analyses of chloroplast genomes of Euphorbia species.

The genus Euphorbia (Euphorbiaceae) has near-cosmopolitan distribution and serves as a significant resource for both ornamental and medicinal purposes. Despite its economic importance, Euphorbia's taxonomy has long been challenged by the intricate nature of morphological traits exhibiting high levels of convergence. While molecular markers are essential for phylogenetic studies, their availability for Euphorbia has been limited. To address this gap, we conducted comparative analyses focusing on the chloroplast (CP) genomes of nine Euphorbia species, incorporating three newly sequenced and annotated accessions. In addition, phylogenetic informativeness and nucleotide diversity were computed to identify candidate markers for phylogenetic analyses among closely related taxa in the genus. Our investigation revealed relatively conserved sizes and structures of CP genomes across the studied species, with notable interspecific variations observed primarily in non-coding regions and IR/SC borders. By leveraging phylogenetic informativeness and nucleotide diversity, we identified rpoB gene as the optimal candidate for species delimitation and shallow-level phylogenetic inference within the genus. Through this comprehensive analysis of CP genomes across multiple taxa, our study sheds light on the evolutionary dynamics and taxonomic intricacies of Euphorbia, offering valuable insights into its CP genome evolution and taxonomy.

Chromosome-level genome assembly of Euphorbia tirucalli (Euphorbiaceae), a highly stress-tolerant oil plant.

Euphorbia, one of the largest genera of flowering plants, is well-known for containing many biofuel crops. Euphorbia tirucalli, an evergreen succulent mainly native to the Africa continent but cultivated worldwide, is a promising petroleum plant with high tolerance to drought and salt stress. However, the exploration of such an important plant resource is severely hampered by the lack of a reference genome. Here, we present the chromosome-level genome assembly of E. tirucalli using PacBio HiFi sequencing and Hi-C technology. Its genome size was approximately 745.62 Mb, with a contig N50 of 74.16 Mb. A total of 743.63 Mb (99.73%) of the assembled sequences were anchored to 10 chromosomes with a complete BUSCO score of 97.80%. Genome annotation revealed 26,304 protein-coding genes, and 76.37% of the genome was identified as repeat elements. The high-quality genome provides valuable genetic resources that would be useful for unraveling the genetic mechanisms of biofuel synthesis and evolutionary adaptation of E. tirucalli.

Intraspecific divergence in a coastal plant, Euphorbia jolkinii, at a major biogeographic boundary in East Asia.

PREMISE: Quaternary climatic fluctuations and long-distance seed dispersal across the sea are critical factors affecting the distribution of coastal plants, but the spatiotemporal nature of population expansion and distribution change of East Asian coastal plants during this period are rarely examined. To explore this process, we investigated the genome-wide phylogenetic patterns of Euphorbia jolkinii Boiss. (Euphorbiaceae), which grows widely on littoral areas of Japan, Korea, and Taiwan. METHODS: We used plastome sequences and genome-wide single nucleotide polymorphisms in samples across the species range to reveal phylogeographic patterns and spatiotemporal distributional changes. We conducted ecological niche modeling for the present and the last glacial maximum (LGM). RESULTS: Genetic differentiation was observed between the northern and southern populations of E. jolkinii, separated by the major biogeographic boundary, the Tokara Gap. These two groups of populations differentiated during the glacial period and subsequently intermingled in the intermorainic areas of the central Ryukyu Islands after the LGM. Ecological niche models suggested that the potential range of E. jolkinii was restricted to southern Kyushu; however, it was widespread in the southern Ryukyu Islands and Taiwan during the LGM. CONCLUSIONS: This study provides evidence of genetic differentiation among coastal plant populations separated by the prominent biogeographical boundary. Although coastal plants are typically expected to maintain population connectivity through sea-drifted seed dispersal, our findings suggest that genetic differences may arise because of a combination of limited gene flow and changes in climate during the glacial period.

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.

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.

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.

Expression Level of a Phenylalanine Ammonia-Lyase Gene in Poinsettia Is Negatively Correlated with Poinsettia Branch-Inducing Phytoplasma Titer.

Poinsettia is an important ornamental cultivated worldwide. Commercial poinsettias are almost universally infected with a pathogen known as the poinsettia branch-inducing phytoplasma (PoiBI), which can increase the level of branching in host plants and make the plants more desirable to consumers. Despite PoiBI's crucial role in poinsettia production, little is known about PoiBI-poinsettia interactions in regard to the pathogen's in planta population dynamics. The expression profiles of a phenylalanine ammonia-lyase gene (Euphorbia pulcherrima PAL [EpPAL]) and the PoiBI titers in poinsettia tissues were investigated. Differential gene expression analyses using quantitative PCR (qPCR) showed that EpPAL expression levels differed significantly across tissue types. The highest expression levels were detected in stems, followed by root. Lower EpPAL expression levels were detected in leaf tissues, particularly in source leaves closer to the base; the average expression level in these leaves was only one-seventh of that detected in stems. Phytoplasma concentrations in source leaves close to the base were significantly greater than the other tissue types; the average value was 7.6-fold of that detected in stem tissues, which had the lowest phytoplasma titers. A negative correlation between EpPAL expression level and PoiBI load was detected, suggesting that the products of EpPAL-associated pathways or other genes indirectly associated with EpPAL may interfere with PoiBI's growth. While additional studies are needed to validate these interpretations, the results from this work provide new insights into PoiBI-poinsettia interaction and showed that correlations between pathogen load and defense-related genes could be detected in phytoplasma-associated pathosystems. IMPORTANCE Phytoplasma-plant interactions are interesting subjects for fundamental and applicative research. Although many studies have characterized molecular interplays between these pathogens and hosts, knowledge on relationships between phytoplasmas' in planta population dynamics and host gene expression remains scarce. By using the poinsettia branch-inducing phytoplasma (PoiBI) and poinsettia as a model system, a negative correlation was observed between the expression level of a plant defense-related gene and the pathogen's titer. The findings provide potential explanations to PoiBI's distribution patterns in the plant and highlight the importance of studying phytoplasma-plant interactions in regard to the pathogen's population dynamics in other pathosystems.

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.

Gene discovery and virus-induced gene silencing reveal branched pathways to major classes of bioactive diterpenoids in Euphorbia peplus.

Most macro- and polycyclic Euphorbiaceae diterpenoids derive from the common C20 precursor casbene. While the biosynthetic pathway from casbene to the lathyrane jolkinol C is characterized, pathways to other more complex classes of bioactive diterpenoids remain to be elucidated. A metabolomics-guided transcriptomic approach and a genomics approach that led to the discovery of two casbene-derived diterpenoid gene clusters yielded a total of 68 candidate genes that were transiently expressed in Nicotiana benthamiana for activity toward jolkinol C and other lathyranes. We report two short-chain dehydrogenases/reductases (SDRs), identified by RNA sequencing to be highly expressed in Euphorbia peplus latex. One of these, EpSDR-5, is a C3-ketoreductase, converting jolkinol C to the lathyrane jolkinol E. Gene function of EpSDR-5 was further confirmed by heterologous expression in Saccharomyces cerevisiae. To investigate the in vivo role of EpSDR-5, we established virus-induced gene silencing (VIGS) in E. peplus, resulting in a significant reduction in jatrophanes and a corresponding increase in ingenanes. VIGS of Casbene Synthase results in a major reduction in both jatrophanes and ingenanes, the two most abundant classes of E. peplus diterpenoids. VIGS of CYP71D365 had a similar effect, consistent with the previously determined role of this gene in the pathway to jolkinol C. These results point to jolkinol C being a branch point intermediate in the pathways to ingenanes and jatrophanes with EpSDR-5 responsible for the first step from jolkinol C to jatrophane production.

High-quality genome assembly of an important biodiesel plant, Euphorbia lathyris L.

Caper spurge, Euphorbia lathyris L., is an important energy crop and medicinal crop. Here, we generated a high-quality, chromosome-level genome assembly of caper spurge using Oxford Nanopore sequencing, Illumina sequencing, and Hi-C technology. The final genome assembly was ∼988.9 Mb in size, 99.8% of which could be grouped into 10 pseudochromosomes, with contig and scaffold N50 values of 32.6 and 95.7 Mb, respectively. A total of 651.4 Mb repetitive sequences and 36,342 protein-coding genes were predicted in the genome assembly. Comparative genomic analysis showed that caper spurge and castor bean clustered together. We found that no independent whole-genome duplication event had occurred in caper spurge after its split from the castor bean, and recent substantial amplification of long terminal repeat retrotransposons has contributed significantly to its genome expansion. Furthermore, based on gene homology searching, we identified a number of candidate genes involved in the biosynthesis of fatty acids and triacylglycerols. The reference genome presented here will be highly useful for the further study of the genetics, genomics, and breeding of this high-value crop, as well as for evolutionary studies of spurge family and angiosperms.

The evolutionary fate of rpl32 and rps16 losses in the Euphorbia schimperi (Euphorbiaceae) plastome.

Gene transfers from mitochondria and plastids to the nucleus are an important process in the evolution of the eukaryotic cell. Plastid (pt) gene losses have been documented in multiple angiosperm lineages and are often associated with functional transfers to the nucleus or substitutions by duplicated nuclear genes targeted to both the plastid and mitochondrion. The plastid genome sequence of Euphorbia schimperi was assembled and three major genomic changes were detected, the complete loss of rpl32 and pseudogenization of rps16 and infA. The nuclear transcriptome of E. schimperi was sequenced to investigate the transfer/substitution of the rpl32 and rps16 genes to the nucleus. Transfer of plastid-encoded rpl32 to the nucleus was identified previously in three families of Malpighiales, Rhizophoraceae, Salicaceae and Passifloraceae. An E. schimperi transcript of pt SOD-1-RPL32 confirmed that the transfer in Euphorbiaceae is similar to other Malpighiales indicating that it occurred early in the divergence of the order. Ribosomal protein S16 (rps16) is encoded in the plastome in most angiosperms but not in Salicaceae and Passifloraceae. Substitution of the E. schimperi pt rps16 was likely due to a duplication of nuclear-encoded mitochondrial-targeted rps16 resulting in copies dually targeted to the mitochondrion and plastid. Sequences of RPS16-1 and RPS16-2 in the three families of Malpighiales (Salicaceae, Passifloraceae and Euphorbiaceae) have high sequence identity suggesting that the substitution event dates to the early divergence within Malpighiales.

A highly mutable GST is essential for bract colouration in Euphorbia pulcherrima Willd. Ex Klotsch.

BACKGROUND: Mutation breeding is an extraordinary tool in plant breeding to increase the genetic variability, where mutations in anthocyanin biosynthesis are targets to generate distinctive phenotypes in ornamental species. In poinsettia, ionizing radiation is routinely applied in breeding programs to obtaining a range of colours, with nearly all pink and white varieties being obtained after γ- or X-ray mutagenesis of red varieties. In the present study we performed a thorough characterization of a potential mutagenesis target gene as the main responsible for the 'white paradox' in poinsettia. RESULTS: We identified a GST gene in poinsettia (Bract1) as an essential factor for the expression of anthocyanin-based red colouration of bracts, which presents a high phylogenetic similarity to known anthocyanin-related GSTs. Red poinsettia varieties and white mutants generated from these varieties by X-ray were analysed for polymorphisms related to the 'white paradox' in the species. A 4 bp mutation in a short repeat within the coding region of Bract1 is most likely responsible for the appearance of white phenotypes upon irradiation treatment. The polymorphism between wild-type and mutant alleles co-segregates with the phenotype in progeny from heterozygous red and white parents. Moreover, overexpression of Bract1 wild-type allele in Arabidopsis tt19 mutants restored the anthocyanin phenotype, while the Bract1 mutated allele showed to be non-functional. CONCLUSIONS: The identified repeat seems to be highly unstable, since mutated plants can be easily detected among fewer than 200 shoots derived from 10 mutated plants. Our data indicate that particular short repeat sequences, similar to microsatellite sequences or so-called dynamic mutations, might be hot spots for genetic variability. Moreover, the identification of the Bract1 mutation fills a gap on the understanding on the molecular mechanism of colour formation in poinsettia.

Large-scale advances in SSR markers with high-throughput sequencing in Euphorbia fischeriana Steud

BACKGROUND: Euphorbia fischeriana Steud is a very important medicinal herb and has significant medical value for healing cancer, edema and tuberculosis in China. The lack of molecular markers for Euphorbia fischeriana Steud is a dominant barrier to genetic research. For the purpose of developing many simple sequence repeat (SSR) molecular markers, we completed transcriptome analysis with the Illumina HiSeq 2000 platform. RESULTS: Approximately 9.1 million clean reads were acquired and then assembled into approximately 186.3 thousand nonredundant unigenes, 53,146 of which were SSR-containing unigenes. A total of 76,193 SSR loci were identified. Of these SSR loci, 28,491 were detected at the terminal position of ESTs, which made it difficult to design SSR primers for these SSR-containing sequences, and the residual SSRs were thus used to design primer pairs. Analyzing the results of these markers revealed that the mononucleotide motif A/T (44,067, 57.83% of all SSRs) was the most abundant, followed by the dinucleotide type AG/CT (9430, 12.38%). Using 100 randomly selected primer pairs, 77 primers were successfully amplified in Euphorbia fischeriana Steud, and 79 were successfully amplified in three other related species. The markers developed displayed relatively high quality and cross-species transferability. CONCLUSIONS: The large number of EST-SSRs exploited successfully in Euphorbia fischeriana Steud for the first time could provide genetic information for research on linkage maps, variety identification, genetic diversity analysis, and molecular marker-assisted breeding.

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.

De novo assembly and characterization of transcriptome in the medicinal plant Euphorbia jolkini.

BACKGROUND: Euphorbia jolkini, a medicinal herb that grows on the warm beaches in Japan and South Korea, is known to be used for traditional medicines to treat a variety of ailments, including bruises, stiffness, indigestion, toothache, and diabetes. OBJECTIVE: It is to analyze the whole transcriptome and identify the genes related to the phenylpropanoid biosynthesis in the medicinally important herb E jolkini. METHODS: Paired-end Illumina HiSeq™ 2500 sequencing technology was employed for cDNA library construction and Illumina sequencing. Public databases like TAIR (The Arabidopsis Information Resource), Swissprot and KEGG (Kyoto Encyclopedia of Genes and Genomes) were used for annotations of unigenes obtained. RESULTS: The transcriptome of E. jolkini generated 139,215 assembled transcripts with an average length of 868 bp and an N50 value of 1460 bp that were further clustered using CD-HIT into 93,801 unigenes with an average length of 847 bp (N50-1410 bp). Sixty-three percent of the coding sequences (CDS) were annotated from the longest open reading frame (ORF). A remarkable percentage of unigenes were annotated against various databases. The differentially expressed gene analysis revealed that the expression of genes related to the terpenoid backbone biosynthesis pathway was higher in the flowers, whereas that of genes related to the phenylpropanoid biosynthesis pathway was both up- and downregulated in flowers and leaves. A search of against the transcription factor domain found 1023 transcription factors (TFs) that were from 54 TF families. CONCLUSION: Assembled sequences of the E. jolkini transcriptome are made available for the first time in this study E. jolkini and lay a foundation for the investigation of secondary metabolite biosynthesis.

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.

Hybrid de novo transcriptome assembly of poinsettia (Euphorbia pulcherrima Willd. Ex Klotsch) bracts.

BACKGROUND: Poinsettia is a popular and important ornamental crop, mostly during the Christmas season. Its bract coloration ranges from pink/red to creamy/white shades. Despite its ornamental value, there is a lack of knowledge about the genetics and molecular biology of poinsettia, especially on the mechanisms of color formation. We performed an RNA-Seq analysis in order to shed light on the transcriptome of poinsettia bracts. Moreover, we analyzed the transcriptome differences of red- and white-bracted poinsettia varieties during bract development and coloration. For the assembly of a bract transcriptome, two paired-end cDNA libraries from a red and white poinsettia pair were sequenced with the Illumina technology, and one library from a red-bracted variety was used for PacBio sequencing. Both short and long reads were assembled using a hybrid de novo strategy. Samples of red- and white-bracted poinsettias were sequenced and comparatively analyzed in three color developmental stages in order to understand the mechanisms of color formation and accumulation in the species. RESULTS: The final transcriptome contains 288,524 contigs, with 33% showing confident protein annotation against the TAIR10 database. The BUSCO pipeline, which is based on near-universal orthologous gene groups, was applied to assess the transcriptome completeness. From a total of 1440 BUSCO groups searched, 77% were categorized as complete (41% as single-copy and 36% as duplicated), 10% as fragmented and 13% as missing BUSCOs. The gene expression comparison between red and white varieties of poinsettia showed a differential regulation of the flavonoid biosynthesis pathway only at particular stages of bract development. An initial impairment of the flavonoid pathway early in the color accumulation process for the white poinsettia variety was observed, but these differences were no longer present in the subsequent stages of bract development. Nonetheless, GSTF11 and UGT79B10 showed a lower expression in the last stage of bract development for the white variety and, therefore, are potential candidates for further studies on poinsettia coloration. CONCLUSIONS: In summary, this transcriptome analysis provides a valuable foundation for further studies on poinsettia, such as plant breeding and genetics, and highlights crucial information on the molecular mechanism of color formation.

Target site as the main mechanism of resistance to imazamox in a Euphorbia heterophylla biotype.

Euphorbia heterophylla is a weed species that invades extensive crop areas in subtropical regions of Brazil. This species was previously controlled by imazamox, but the continuous use of this herbicide has selected for resistant biotypes. Two biotypes of E. heterophylla from southern Brazil, one resistant (R) and one susceptible (S) to imazamox, were compared. The resistance of the R biotype was confirmed by dose-response assays since it required 1250.2 g ai ha-1 to reduce the fresh weight by 50% versus 7.4 g ai ha-1 for the S biotype. The acetolactate synthase (ALS) enzyme activity was studied using ALS-inhibiting herbicides from five different chemical families. The R biotype required the highest concentrations to reduce this enzyme activity by 50%. A Ser653Asn mutation was found in the ALS gene of the R biotype. The experiments carried out showed that imazamox absorption and metabolism were not involved in resistance. However, greater 14C-imazamox root exudation was found in the R biotype (~70% of the total absorbed imazamox). Target site mutation in the ALS gene is the principal mechanism that explains the imazamox resistance of the R biotype, but root exudation seems to also contribute to the resistance of this biotype.