Traditional medicinal use is linked with apparency, not specialized metabolite profiles in the order Caryophyllales.

PREMISE: Better understanding of the relationship between plant specialized metabolism and traditional medicine has the potential to aid in bioprospecting and untangling of cross-cultural use patterns. However, given the limited information available for metabolites in most plant species, understanding medicinal use-metabolite relationships can be difficult. The order Caryophyllales has a unique pattern of lineages of tyrosine- or phenylalanine-dominated specialized metabolism, represented by mutually exclusive anthocyanin and betalain pigments, making Caryophyllales a compelling system to explore the relationship between medicine and metabolites by using pigment as a proxy for dominant metabolism. METHODS: We compiled a list of medicinal species in select tyrosine- or phenylalanine-dominant families of Caryophyllales (Nepenthaceae, Polygonaceae, Simmondsiaceae, Microteaceae, Caryophyllaceae, Amaranthaceae, Limeaceae, Molluginaceae, Portulacaceae, Cactaceae, and Nyctaginaceae) by searching scientific literature until no new uses were recovered. We then tested for phylogenetic clustering of uses using a "hot nodes" approach. To test potential non-metabolite drivers of medicinal use, like how often humans encounter a species (apparency), we repeated the analysis using only North American species across the entire order and performed phylogenetic generalized least squares regression (PGLS) with occurrence data from the Global Biodiversity Information Facility (GBIF). RESULTS: We hypothesized families with tyrosine-enriched metabolism would show clustering of different types of medicinal use compared to phenylalanine-enriched metabolism. Instead, wide-ranging, apparent clades in Polygonaceae and Amaranthaceae are overrepresented across nearly all types of medicinal use. CONCLUSIONS: Our results suggest that apparency is a better predictor of medicinal use than metabolism, although metabolism type may still be a contributing factor.

Evolutionary paths to new phenotypes.

Ecological model systems inform on innovative traits in plants and animals.

Convergence in carnivorous pitcher plants reveals a mechanism for composite trait evolution.

Composite traits involve multiple components that, only when combined, gain a new synergistic function. Thus, how they evolve remains a puzzle. We combined field experiments, microscopy, chemical analyses, and laser Doppler vibrometry with comparative phylogenetic analyses to show that two carnivorous Nepenthes pitcher plant species independently evolved similar adaptations in three distinct traits to acquire a new, composite trapping mechanism. Comparative analyses suggest that this new trait arose convergently through "spontaneous coincidence" of the required trait combination, rather than directional selection in the component traits. Our results indicate a plausible mechanism for composite trait evolution and highlight the importance of stochastic phenotypic variation as a facilitator of evolutionary novelty.

Multiple mechanisms explain loss of anthocyanins from betalain-pigmented Caryophyllales, including repeated wholesale loss of a key anthocyanidin synthesis enzyme.

In this study, we investigate the genetic mechanisms responsible for the loss of anthocyanins in betalain-pigmented Caryophyllales, considering our hypothesis of multiple transitions to betalain pigmentation. Utilizing transcriptomic and genomic datasets across 357 species and 31 families, we scrutinize 18 flavonoid pathway genes and six regulatory genes spanning four transitions to betalain pigmentation. We examined evidence for hypotheses of wholesale gene loss, modified gene function, altered gene expression, and degeneration of the MBW (MYB-bHLH-WD40) trasnscription factor complex, within betalain-pigmented lineages. Our analyses reveal that most flavonoid synthesis genes remain conserved in betalain-pigmented lineages, with the notable exception of TT19 orthologs, essential for the final step in anthocyanidin synthesis, which appear to have been repeatedly and entirely lost. Additional late-stage flavonoid pathway genes upstream of TT19 also manifest strikingly reduced expression in betalain-pigmented species. Additionally, we find repeated loss and alteration in the MBW transcription complex essential for canonical anthocyanin synthesis. Consequently, the loss and exclusion of anthocyanins in betalain-pigmented species appear to be orchestrated through several mechanisms: loss of a key enzyme, downregulation of synthesis genes, and degeneration of regulatory complexes. These changes have occurred iteratively in Caryophyllales, often coinciding with evolutionary transitions to betalain pigmentation.

DNA Barcoding, Phylogenetic Analysis and Secondary Structure Predictions of Nepenthes ampullaria, Nepenthes gracilis and Nepenthes rafflesiana.

Nepentheceae, the most prominent carnivorous family in the Caryophyllales order, comprises the Nepenthes genus, which has modified leaf trap characteristics. Although most Nepenthes species have unique morphologies, their vegetative stages are identical, making identification based on morphology difficult. DNA barcoding is seen as a potential tool for plant identification, with small DNA segments amplified for species identification. In this study, three barcode loci; ribulose-bisphosphate carboxylase (rbcL), intergenic spacer 1 (ITS1) and intergenic spacer 2 (ITS2) and the usefulness of the ITS1 and ITS2 secondary structure for the molecular identification of Nepenthes species were investigated. An analysis of barcodes was conducted using BLASTn, pairwise genetic distance and diversity, followed by secondary structure prediction. The findings reveal that PCR and sequencing were both 100% successful. The present study showed the successful amplification of all targeted DNA barcodes at different sizes. Among the three barcodes, rbcL was the least efficient as a DNA barcode compared to ITS1 and ITS2. The ITS1 nucleotide analysis revealed that the ITS1 barcode had more variations compared to ITS2. The mean genetic distance (K2P) between them was higher for interspecies compared to intraspecies. The results showed that the DNA barcoding gap existed among Nepenthes species, and differences in the secondary structure distinguish the Nepenthes. The secondary structure generated in this study was found to successfully discriminate between the Nepenthes species, leading to enhanced resolutions.

Two independently evolved natural mutations additively deregulate TyrA enzymes and boost tyrosine production in planta.

l-Tyrosine is an essential amino acid for protein synthesis and is also used in plants to synthesize diverse natural products. Plants primarily synthesize tyrosine via TyrA arogenate dehydrogenase (TyrAa or ADH), which are typically strongly feedback inhibited by tyrosine. However, two plant lineages, Fabaceae (legumes) and Caryophyllales, have TyrA enzymes that exhibit relaxed sensitivity to tyrosine inhibition and are associated with elevated production of tyrosine-derived compounds, such as betalain pigments uniquely produced in core Caryophyllales. Although we previously showed that a single D222N substitution is primarily responsible for the deregulation of legume TyrAs, it is unknown when and how the deregulated Caryophyllales TyrA emerged. Here, through phylogeny-guided TyrA structure-function analysis, we found that functionally deregulated TyrAs evolved early in the core Caryophyllales before the origin of betalains, where the E208D amino acid substitution in the active site, which is at a different and opposite location from D222N found in legume TyrAs, played a key role in the TyrA functionalization. Unlike legumes, however, additional substitutions on non-active site residues further contributed to the deregulation of TyrAs in Caryophyllales. The introduction of a mutation analogous to E208D partially deregulated tyrosine-sensitive TyrAs, such as Arabidopsis TyrA2 (AtTyrA2). Moreover, the combined introduction of D222N and E208D additively deregulated AtTyrA2, for which the expression in Nicotiana benthamiana led to highly elevated accumulation of tyrosine in planta. The present study demonstrates that phylogeny-guided characterization of key residues underlying primary metabolic innovations can provide powerful tools to boost the production of essential plant natural products.

The jojoba genome reveals wide divergence of the sex chromosomes in a dioecious plant.

Most flowering plants are hermaphrodites, but around 6% of species are dioecious, having separate male and female plants. Sex chromosomes and some sex-specific genes have been reported in plants, but the genome sequences have not been compared. We now report the genome sequence of male and female jojoba (Simmondsia chinensis) plants, revealing a very large difference in the sex chromosomes. The male genome assembly was 832 Mb and the female 822 Mb. This was explained by the large size differences in the Y chromosome (37.6 Mb) compared with the X chromosome (26.9 Mb). Relative to the X chromosome, the Y chromosome had two large insertions each of more than 5 Mb containing more than 400 genes. Many of the genes in the chromosome-specific regions were novel. These male-specific regions included many flowering-related and stress response genes. Smaller insertions found only in the X chromosome totalled 877 kb. The wide divergence of the sex chromosomes suggests a long period of adaptation to diverging sex-specific roles. Male and female plants may have evolved to accommodate factors such as differing reproductive resource allocation requirements under the stress of the desert environment in which the plants are found. The sex-determining regions accumulate genes beneficial to each sex. This has required the evolution of many more novel sex-specific genes than has been reported for other organisms. This suggest that dioecious plants provide a novel source of genes for manipulation of reproductive performance and environmental adaptation in crops.

Anthocyanin synthesis potential in betalain-producing Caryophyllales plants.

Although anthocyanins are widely distributed in higher plants, betalains have replaced anthocyanins in most species of the order Caryophyllales. The accumulation of flavonols in Caryophyllales plants implies that the late step of anthocyanin biosynthesis from dihydroflavonols to anthocyanins may be blocked in Caryophyllales. The isolation and characterization of functional dihydroflavonol 4-reductase (DFR) and anthocyanidin synthase (ANS) from Caryophyllales plants has indicated a lack of anthocyanins due to suppression of DFR and ANS. In this study, we demonstrated that overexpression of DFR and ANS from Spinacia oleracea (SoDFR and SoANS, respectively) with PhAN9, which encodes glutathione S-transferase (required for anthocyanin sequestration) from Petunia induces ectopic anthocyanin accumulation in yellow tepals of the cactus Astrophytum myriostigma. A promoter assay of SoANS showed that the Arabidopsis MYB transcription factor PRODUCTION OF ANTHOCYANIN PIGMENT1 (PAP1) activated the SoANS promoter in Arabidopsis leaves. The overexpression of Arabidopsis transcription factors with PhAN9 also induced ectopic anthocyanin accumulation in yellow cactus tepals. PAP homologs from betalain-producing Caryophyllales did not activate the promoter of ANS. In-depth characterization of Caryophyllales PAPs and site-directed mutagenesis in the R2R3-MYB domains identified the amino acid residues affecting transactivation of Caryophyllales PAPs. The substitution of amino acid residues recovered the transactivation ability of Caryophyllales PAPs. Therefore, loss of function in MYB transcription factors may suppress expression of genes involved in the late stage of anthocyanin synthesis, resulting in a lack of anthocyanin in betalain-producing Caryophyllales plants.

Pitchers of Nepenthes khasiana express several digestive-enzyme encoding genes, harbor mostly fungi and probably evolved through changes in the expression of leaf polarity genes.

BACKGROUND: A structural phenomenon seen in certain lineages of angiosperms that has captivated many scholars including Charles Darwin is the evolution of plant carnivory. Evidently, these structural features collectively termed carnivorous syndrome, evolved to aid nutritional acquisition from attracted, captured and digested prey. We now understand why plant carnivory evolved but how carnivorous plants acquired these attributes remains a mystery. In an attempt to understand the evolution of Nepenthes pitcher and to shed more light on its role in prey digestion, we analyzed the transcriptome data of the highly specialized Nepenthes khasiana leaf comprising the leaf base lamina, tendril and the different parts/zones of the pitcher tube viz. digestive zone, waxy zone and lid. RESULTS: In total, we generated around 262 million high-quality Illumina reads. Reads were pooled, normalized and de novo assembled to generate a reference transcriptome of about 412,224 transcripts. We then estimated transcript abundance along the N. khasiana leaf by mapping individual reads from each part/zone to the reference transcriptome. Correlation-based hierarchical clustering analysis of 27,208 commonly expressed genes indicated functional relationship and similar cellular processes underlying the development of the leaf base and the pitcher, thereby implying that the Nepenthes pitcher is indeed a modified leaf. From a list of 2386 differentially expressed genes (DEGs), we identified transcripts encoding key enzymes involved in prey digestion and protection against pathogen attack, some of which are expressed at high levels in the digestive zone. Interestingly, many of these enzyme-encoding genes are also expressed in the unopened N. khasiana pitcher. Transcripts showing homology to both bacteria and fungi were also detected; and in the digestive zone, fungi are more predominant as compared to bacteria. Taking cues from histology and scanning electron microscopy (SEM) photomicrographs, we found altered expressions of key regulatory genes involved in leaf development. Of particular interest, the expression of class III HOMEODOMAIN-LEUCINE ZIPPER (HD-ZIPIII) and ARGONAUTE (AGO) genes were upregulated in the tendril. CONCLUSIONS: Our findings suggest that N. khasiana pitchers employ a wide range of enzymes for prey digestion and plant defense, harbor microbes and probably evolved through altered expression of leaf polarity genes.

In vitro plant regeneration and Agrobacterium-mediated genetic transformation of a carnivorous plant, Nepenthes mirabilis.

In nutrient-poor habitats, carnivorous plants have developed novel feeding strategies based on the capture and digestion of prey and the assimilation of prey-derived nutrients by specialized traps. The Nepenthes genus, comprising nearly 160 species, presents a remarkable pitcher-shaped trap, leading to great interest among biologists, but the species of this genus are listed as threatened. In this work, we developed a protocol for reproducing Nepenthes mirabilis through shoot regeneration from calli. The cultivation of stem segments of N. mirabilis on MS medium containing thidiazuron induced organogenic calli after 10 weeks. Subcultured calli exposed to 6-benzylaminopurine showed shoot regeneration in 3 weeks with considerable yields (143 shoots/g of calli). Excised shoots transferred to medium with indole-3-butyric acid allowed rooting in 4 weeks, and rooted plantlets had a 100% survival rate. Based on this method, we also developed an Agrobacterium-mediated genetic transformation protocol using calli as explants and ipt as a positive method of selection. Twelve weeks post infection, regenerated shoots were observed at the surface of calli. Their transgenic status was confirmed by PCR and RT-PCR. In conclusion, this study provides an efficient method for regenerating Nepenthes and the first protocol for its stable genetic transformation, a new tool for studying carnivory.

The genome of jojoba (Simmondsia chinensis): A taxonomically isolated species that directs wax ester accumulation in its seeds.

Seeds of the desert shrub, jojoba (Simmondsia chinensis), are an abundant, renewable source of liquid wax esters, which are valued additives in cosmetic products and industrial lubricants. Jojoba is relegated to its own taxonomic family, and there is little genetic information available to elucidate its phylogeny. Here, we report the high-quality, 887-Mb genome of jojoba assembled into 26 chromosomes with 23,490 protein-coding genes. The jojoba genome has only the whole-genome triplication (γ) shared among eudicots and no recent duplications. These genomic resources coupled with extensive transcriptome, proteome, and lipidome data helped to define heterogeneous pathways and machinery for lipid synthesis and storage, provided missing evolutionary history information for this taxonomically segregated dioecious plant species, and will support efforts to improve the agronomic properties of jojoba.

A phylogenomic analysis of Nepenthes (Nepenthaceae).

Nepenthaceae is one of the largest carnivorous plant families and features ecological and morphological adaptations indicating an impressive adaptive radiation. However, investigation of evolutionary and taxonomic questions is hindered by poor phylogenetic understanding, with previous molecular studies based on limited loci and taxa. We use high-throughput sequencing with a target-capture methodology based on a 353-loci, probe set to recover sequences for 197 samples, representing 151 described or putative Nepenthes species. Phylogenetic analyses were performed using supermatrix and maximum quartet species tree approaches. Our analyses confirm five Western outlier taxa, followed by N. danseri, as successively sister to the remainder of the group. We also find mostly consistent recovery of two major Southeast Asian clades. The first contains common or widespread lowland species plus a Wallacean-New Guinean clade. Within the second clade, sects. Insignes and Tentaculatae are well supported, while geographically defined clades representing Sumatra, Indochina, Peninsular Malaysia, Palawan, Mindanao and Borneo are also consistently recovered. However, we find considerable conflicting signal at the site and locus level, and often unstable backbone relationships. A handful of Bornean taxa are inconsistently placed and require further investigation. We make further suggestions for a modified infra-generic classification of genus Nepenthes.

Protein replenishment in pitcher fluids of Nepenthes × ventrata revealed by quantitative proteomics (SWATH-MS) informed by transcriptomics.

Carnivorous plants capture and digest insects for nutrients, allowing them to survive in soil deprived of nitrogenous nutrients. Plants from the genus Nepenthes produce unique pitchers containing secretory glands, which secrete enzymes into the digestive fluid. We performed RNA-seq analysis on the pitcher tissues and liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis on the pitcher fluids of Nepenthes × ventrata to study protein expression in this carnivory organ during early days of pitcher opening. This transcriptome provides a sequence database for pitcher fluid protein identification. A total of 32 proteins of diverse functions were successfully identified in which 19 proteins can be quantified based on label-free quantitative proteomics (SWATH-MS) analysis while 16 proteins were not reported previously. Our findings show that certain proteins in the pitcher fluid were continuously secreted or replenished after pitcher opening, even without any prey or chitin induction. We also discovered a new aspartic proteinase, Nep6, secreted into pitcher fluid. This is the first SWATH-MS analysis of protein expression in Nepenthes pitcher fluid using a species-specific reference transcriptome. Taken together, our study using a gel-free shotgun proteomics informed by transcriptomics (PIT) approach showed the dynamics of endogenous protein secretion in the digestive organ of N. × ventrata and provides insights on protein regulation during early pitcher opening prior to prey capture.

The evolution of betalain biosynthesis in Caryophyllales.

Within the angiosperm order Caryophyllales, an unusual class of pigments known as betalains can replace the otherwise ubiquitous anthocyanins. In contrast to the phenylalanine-derived anthocyanins, betalains are tyrosine-derived pigments which contain the chromophore betalamic acid. The origin of betalain pigments within Caryophyllales and their mutual exclusion with anthocyanin pigments have been the subject of considerable research. In recent years, numerous discoveries, accelerated by -omic scale data, phylogenetics and synthetic biology, have shed light on the evolution of the betalain biosynthetic pathway in Caryophyllales. These advances include the elucidation of the biosynthetic steps in the betalain pathway, identification of transcriptional regulators of betalain synthesis, resolution of the phylogenetic history of key genes, and insight into a role for modulation of primary metabolism in betalain synthesis. Here we review how molecular genetics have advanced our understanding of the betalain biosynthetic pathway, and discuss the impact of phylogenetics in revealing its evolutionary history. In light of these insights, we explore our new understanding of the origin of betalains, the mutual exclusion of betalains and anthocyanins, and the homoplastic distribution of betalain pigmentation within Caryophyllales. We conclude with a speculative conceptual model for the stepwise emergence of betalain pigmentation.

Identification and Expression Analysis of Chitinase Genes in Pitchers of Nepenthes sp. during Development.

Fifteen chitinases of classes I-V were identified in the transcriptomes of pitchers and adult leaves of the carnivorous plant Nepenthes sp. Ten of these chitinases were identified for the first time, including the chitinases of classes II and V. The expression levels of all found chitinase genes in leaves and at three stages of pitcher development were determined. The maximum level of transcriptional activity in an open pitcher was observed for the genes encoding chitinase NChi4 (class II) and its isoforms. The expression levels of these genes significantly increased as the pitcher developed. In addition, for the first time, transcription of the genes encoding chitinases of all five classes was detected in the leaves of this plant.

ASYMMETRIC LEAVES1 and REVOLUTA are the key regulatory genes associated with pitcher development in Nepenthes khasiana.

Nepenthes develops highly specialized insect-eating organs called pitchers that provide adequate insect-derived nutrients to the plants to offset low nutrient availability in their natural habitat. But so far, the molecular basis of Nepenthes pitcher development remains largely unknown. In an attempt to unravel the underlying mechanisms of pitcher formation, we made morphological observations of the developing N. khasiana leaf and performed RNA-seq to identify genes controlling pitcher development. Histology and scanning electron microscopy photomicrographs show that pitcher formation in N. khasiana occurs early in development and shares anatomical features with the young in-rolled leaf base lamina. Analysis of the RNA-seq data indicated that the modification of the leaf into a pitcher is associated with the altered expressions of leaf polarity genes ASYMMETRIC LEAVES1 (AS1) and REVOLUTA (REV). In fact, both genes displayed exclusive or relatively higher expressions in the tip of the leaf that later developed into a pitcher. We propose that NkAS1 may act to inhibit lamina outgrowth and promote the formation of the tendril. Increased NkREV expression may have been involved in the formation of the N. khasiana pitcher. This dataset will allow further research into this area and serve as the basis for understanding Nepenthes pitcher development.

Plastid phylogenomic insights into the evolution of Caryophyllales.

The Caryophyllales includes 40 families and 12,500 species, representing a large and diverse clade of angiosperms. Collectively, members of the clade grow on all continents and in all terrestrial biomes and often occupy extreme habitats (e.g., xeric, salty). The order is characterized by many taxa with unusual adaptations including carnivory, halophytism, and multiple origins of C4 photosynthesis. However, deep phylogenetic relationships within the order have long been problematic due to putative rapid divergence. To resolve the deep-level relationships of Caryophyllales, we performed phylogenomic analyses of all 40 families of Caryophyllales. We time-calibrated the molecular phylogeny of this clade, and evaluated putative correlations among plastid structural changes and rates of molecular substitution. We recovered a well-resolved and well-supported phylogeny of the Caryophyllales that was largely congruent with previous estimates of this order. Our results provide improved support for the phylogenetic position of several key families within this clade. The crown age of Caryophyllales was estimated at ca. 114.4 million years ago (Ma), with periods of rapid divergence in the mid-Cretaceous. A strong, positive correlation between nucleotide substitution rate and plastid structural changes was detected. Our study highlights the importance of broad taxon sampling in phylogenomic inference and provides a firm basis for future investigations of molecular, morphological, and ecophysiological evolution in Caryophyllales.

Comparative cytogenetics of the ACPT clade (Anacampserotaceae, Cactaceae, Portulacaceae, and Talinaceae): a very diverse group of the suborder Cactineae, Caryophyllales.

The clade ACPT (Anacampserotaceae, Cactaceae, Portulacaceae, and Talinaceae) is the most diverse lineage of the subordem Cactineae. The relationships between these families are still uncertain, with different topologies suggested by phylogenetic analyses with several combinations of markers. Different basic numbers (x) have been suggested for each family and for the subord, often in a contestable way. Comparative cytogenetic has helped to understand the evolutionary relationships of phylogenetically poorly resolved groups, as well as their mechanisms of karyotype evolution. The karyotype evolution in representatives of Cactineae was analyzed, focusing on the ACPT clade, through the analysis of chromosome number in a phylogenetic bias. The phylogeny obtained showed a well-resolved topology with support for the monophyly of the five families. Although a chromosomal number is known for less than 30% of the Cactineae species, the analyses revealed a high karyotype variability, from 2n = 8 to 2n = 110. The analysis of character reconstruction of the ancestral haploid numbers (p) suggested p = 12 for Cactineae, with distinct basic numbers for the clade family ACPT: Cactaceae and Montiaceae (p = 11), Talinaceae (p = 12), and Anacampserotaceae and Portulacaceae (p = 9). Talinaceae, Anacampserotaceae, and Cactaceae were stable, while Portulaca and Montiaceae were karyotypically variable. The chromosome evolution of this group was mainly due to events of descending disploidy and poliploidy. Our data confirm that the low phylogenetic resolution among the families of the ACPT clade is due to a divergence of this clade in a short period of time. However, each of these families can be characterized by basic chromosome numbers and unique karyotype evolution events.

Transcriptomic de novo analysis of pitaya (Hylocereus polyrhizus) canker disease caused by Neoscytalidium dimidiatum.

BACKGROUND: Canker disease caused by Neoscytalidium dimidiatum is the most serious disease that attacks the pitaya industry. One pathogenic fungus, referred to as ND8, was isolated from the wild-type red-fleshed pitaya (Hylocereus polyrhizus) of Hainan Province. In the early stages of this disease, stems show little spots and a loss of green color. These spots then gradually spread until the stems became rotten due to infection by various strains. Canker disease caused by Neoscytalidium dimidiatum poses a significant threat to pitaya commercial plantations with the growth of stems and the yields, quality of pitaya fruits. However, a lack of transcriptomic and genomic information hinders our understanding of the molecular mechanisms underlying the pitaya defense response. RESULTS: We investigated the host responses of red-fleshed pitaya (H. polyrhizus) cultivars against N. dimidiatum using Illumina RNA-Seq technology. Significant expression profiles of 23 defense-related genes were further analyzed by qRT-PCR. The total read length based on RNA-Seq was 25,010,007; mean length was 744, the N50 was 1206, and the guanine-cytosine content was 44.48%. Our investigation evaluated 33,584 unigenes, of which 6209 (18.49%) and 27,375 (81.51%) were contigs and singlets, respectively. These unigenes shared a similarity of 16.62% with Vitis vinifera, 7.48% with Theobroma cacao, 6.6% with Nelumbo nucifera and 5.35% with Jatropha curcas. The assembled unigenes were annotated into non-redundant (NR, 25161 unigenes), Kyoto Encyclopedia of Genes and Genomes (KEGG, 17895 unigenes), Clusters of Orthologous Groups (COG, 10475 unigenes), InterPro (19,045 unigenes), and Swiss-Prot public protein databases (16,458 unigenes). In addition, 24 differentially expressed genes, which were mainly associated with plant pathology pathways, were analyzed in-depth. CONCLUSIONS: This study provides a basis for further in-depth research on the protein function of the annotated unigene assembly with cDNA sequences.

Evolution of Portulacineae Marked by Gene Tree Conflict and Gene Family Expansion Associated with Adaptation to Harsh Environments.

Several plant lineages have evolved adaptations that allow survival in extreme and harsh environments including many families within the plant clade Portulacineae (Caryophyllales) such as the Cactaceae, Didiereaceae, and Montiaceae. Here, using newly generated transcriptomic data, we reconstructed the phylogeny of Portulacineae and examined potential correlates between molecular evolution and adaptation to harsh environments. Our phylogenetic results were largely congruent with previous analyses, but we identified several early diverging nodes characterized by extensive gene tree conflict. For particularly contentious nodes, we present detailed information about the phylogenetic signal for alternative relationships. We also analyzed the frequency of gene duplications, confirmed previously identified whole genome duplications (WGD), and proposed a previously unidentified WGD event within the Didiereaceae. We found that the WGD events were typically associated with shifts in climatic niche but did not find a direct association with WGDs and diversification rate shifts. Diversification shifts occurred within the Portulacaceae, Cactaceae, and Anacampserotaceae, and whereas these did not experience WGDs, the Cactaceae experienced extensive gene duplications. We examined gene family expansion and molecular evolutionary patterns with a focus on genes associated with environmental stress responses and found evidence for significant gene family expansion in genes with stress adaptation and clades found in extreme environments. These results provide important directions for further and deeper examination of the potential links between molecular evolutionary patterns and adaptation to harsh environments.