Data access for the 1,000 Plants (1KP) project

Plant peptides - redefining an area of ribosomally synthesized and post-translationally modified peptides

Covering 1965 to February 2024Plants are prolific peptide chemists and are known to make thousands of different peptidic molecules. These peptides vary dramatically in their size, chemistry, and bioactivity. Despite their differences, all plant peptides to date are biosynthesized as ribosomally synthesized and post-translationally modified peptides (RiPPs). Decades of research in plant RiPP biosynthesis have extended the definition and scope of RiPPs from microbial sources, establishing paradigms and discovering new families of biosynthetic enzymes. The discovery and elucidation of plant peptide pathways is challenging due to repurposing and evolution of housekeeping genes as both precursor peptides and biosynthetic enzymes and due to the low rates of gene clustering in plants. In this review, we highlight the chemistry, biosynthesis, and function of the known RiPP classes from plants and recommend a nomenclature for the recent addition of BURP-domain-derived RiPPs termed burpitides. Burpitides are an emerging family of cyclic plant RiPPs characterized by macrocyclic crosslinks between tyrosine or tryptophan side chains and other amino acid side chains or their peptide backbone that are formed by copper-dependent BURP-domain-containing proteins termed burpitide cyclases. Finally, we review the discovery of plant RiPPs through bioactivity-guided, structure-guided, and gene-guided approaches.

Shifts in evolutionary lability underlie independent gains and losses of root-nodule symbiosis in a single clade of plants

Root nodule symbiosis (RNS) is a complex trait that enables plants to access atmospheric nitrogen converted into usable forms through a mutualistic relationship with soil bacteria. Pinpointing the evolutionary origins of RNS is critical for understanding its genetic basis, but building this evolutionary context is complicated by data limitations and the intermittent presence of RNS in a single clade of ca. 30,000 species of flowering plants, i.e., the nitrogen-fixing clade (NFC). We developed the most extensive de novo phylogeny for the NFC and an RNS trait database to reconstruct the evolution of RNS. Our analysis identifies evolutionary rate heterogeneity associated with a two-step process: An ancestral precursor state transitioned to a more labile state from which RNS was rapidly gained at multiple points in the NFC. We illustrate how a two-step process could explain multiple independent gains and losses of RNS, contrary to recent hypotheses suggesting one gain and numerous losses, and suggest a broader phylogenetic and genetic scope may be required for genome-phenome mapping.

Predicting photosynthetic pathway from anatomy using machine learning

Plants with Crassulacean acid metabolism (CAM) have long been associated with a specialized anatomy, including succulence and thick photosynthetic tissues. Firm, quantitative boundaries between non-CAM and CAM plants have yet to be established - if they indeed exist. Using novel computer vision software to measure anatomy, we combined new measurements with published data across flowering plants. We then used machine learning and phylogenetic comparative methods to investigate relationships between CAM and anatomy. We found significant differences in photosynthetic tissue anatomy between plants with differing CAM phenotypes. Machine learning-based classification was over 95% accurate in differentiating CAM from non-CAM anatomy, and had over 70% recall of distinct CAM phenotypes. Phylogenetic least squares regression and threshold analyses revealed that CAM evolution was significantly correlated with increased mesophyll cell size, thicker leaves, and decreased intercellular airspace. Our findings suggest that machine learning may be used to aid the discovery of new CAM species and that the evolutionary trajectory from non-CAM to strong, obligate CAM requires continual anatomical specialization.

Plant genome information facilitates plant functional genomics

MAIN CONCLUSION

In this review, we give an overview of plant sequencing efforts and how this impacts plant functional genomics research. Plant genome sequence information greatly facilitates the studies of plant biology, functional genomics, evolution of genomes and genes, domestication processes, phylogenetic relationships, among many others. More than two decades of sequencing efforts have boosted the number of available sequenced plant genomes. The first plant genome, of Arabidopsis, was published in the year 2000 and currently, 4604 plant genomes from 1482 plant species have been published. Various large sequence initiatives are running, which are planning to produce tens of thousands of sequenced plant genomes in the near future. In this review, we give an overview on the status of sequenced plant genomes and on the use of genome information in different research areas.

Molecular and phylogenetic evidence of parallel expansion of anion channels in plants

Aluminum-activated malate transporters (ALMTs) and slow anion channels (SLACs) are important in various physiological processes in plants, including stomatal regulation, nutrient uptake, and in response to abiotic stress such as aluminum toxicity. To understand their evolutionary history and functional divergence, we conducted phylogenetic and expression analyses of ALMTs and SLACs in green plants. Our findings from phylogenetic studies indicate that ALMTs and SLACs may have originated from green algae and red algae, respectively. The ALMTs of early land plants and charophytes formed a monophyletic clade consisting of three subgroups. A single duplication event of ALMTs was identified in vascular plants and subsequent duplications into six clades occurred in angiosperms, including an identified clade, 1-1. The ALMTs experienced gene number losses in clades 1-1 and 2-1 and expansions in clades 1-2 and 2-2b. Interestingly, the expansion of clade 1-2 was also associated with higher expression levels compared to genes in clades that experienced apparent loss. SLACs first diversified in bryophytes, followed by duplication in vascular plants, giving rise to three distinct clades (I, II, and III), and clade II potentially associated with stomatal control in seed plants. SLACs show losses in clades II and III without substantial expansion in clade I. Additionally, ALMT clade 2-2 and SLAC clade III contain genes specifically expressed in reproductive organs and roots in angiosperms, lycophytes, and mosses, indicating neofunctionalization. In summary, our study demonstrates the evolutionary complexity of ALMTs and SLACs, highlighting their crucial role in the adaptation and diversification of vascular plants.

Adaptive evolution of chloroplast division mechanisms during plant terrestrialization

Despite extensive research, the origin and evolution of the chloroplast division machinery remain unclear. Here, we employ recently sequenced genomes and transcriptomes of Archaeplastida clades to identify the core components of chloroplast division and reconstruct their evolutionary histories, respectively. Our findings show that complete division ring structures emerged in Charophytes. We find that Glaucophytes experienced strong selection pressure, generating diverse variants adapted to the changing terrestrial environments. By integrating the functions of chloroplast division genes (CDGs) annotated in a workflow developed using large-scale multi-omics data, we further show that dispersed duplications acquire more species-specific functions under stronger selection pressures. Notably, PARC6, a dispersed duplicate CDG, regulates leaf color and plant growth in Solanum lycopersicum, demonstrating neofunctionalization. Our findings provide an integrated perspective on the functional evolution of chloroplast division machinery and highlight the potential of dispersed duplicate genes as the primary source of adaptive evolution of chloroplast division.

Comparative analysis of SPL transcription factors from streptophyte algae and embryophytes reveals evolutionary trajectories of SPL family in streptophytes

SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) genes encode plant-specific transcription factors which are important regulators of diverse plant developmental processes. We took advantage of available genome sequences of streptophyte algae representatives to investigate the relationships of SPL genes between freshwater green algae and land plants. Our analysis showed that streptophyte algae, hornwort and liverwort genomes encode from one to four SPL genes which is the smallest set, in comparison to other land plants studied to date. Based on the phylogenetic analysis, four major SPL phylogenetic groups were distinguished with Group 3 and 4 being sister to Group 1 and 2. Comparative motif analysis revealed conserved protein motifs within each phylogenetic group and unique bryophyte-specific motifs within Group 1 which suggests lineage-specific protein speciation processes. Moreover, the gene structure analysis also indicated the specificity of each by identifying differences in exon-intron structures between the phylogenetic groups, suggesting their evolutionary divergence. Since current understanding of SPL genes mostly arises from seed plants, the presented comparative and phylogenetic analyzes from freshwater green algae and land plants provide new insights on the evolutionary trajectories of the SPL gene family in different classes of streptophytes.

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.

A reference genome for the long-term kleptoplast-retaining sea slug Elysia crispata morphotype clarki

Several species of sacoglossan sea slugs possess the incredible ability to sequester chloroplasts from the algae they consume. These "photosynthetic animals" incorporate stolen chloroplasts, called kleptoplasts, into the epithelial cells of tubules that extend from their digestive tracts throughout their bodies. The mechanism by which these slugs maintain functioning kleptoplasts in the absence of an algal nuclear genome is unknown. Here, we report a draft genome of the sacoglossan slug Elysia crispata morphotype clarki, a morphotype native to the Florida Keys that can retain photosynthetically active kleptoplasts for several months without feeding. We used a combination of Oxford Nanopore Technologies long reads and Illumina short reads to produce a 786-Mb assembly (N50 = 0.459 Mb) containing 68,514 predicted protein-coding genes. A phylogenetic analysis found no evidence of horizontal acquisition of genes from algae. We performed gene family and gene expression analyses to identify E. crispata genes unique to kleptoplast-containing slugs that were more highly expressed in fed versus unfed developmental life stages. Consistent with analyses in other kleptoplastic slugs, our investigation suggests that genes encoding lectin carbohydrate-binding proteins and those involved in regulation of reactive oxygen species and immunity may play a role in kleptoplast retention. Lastly, we identified four polyketide synthase genes that could potentially encode proteins producing UV- and oxidation-blocking compounds in slug cell membranes. The genome of E. crispata is a quality resource that provides potential targets for functional analyses and enables further investigation into the evolution and mechanisms of kleptoplasty in animals.

The construction and optimization of engineered yeast chassis for efficient biosynthesis of 8-hydroxygeraniol

Microbial production of monoterpenoid indole alkaloids (MIAs) provides a sustainable and eco-friendly means to obtain compounds with high pharmaceutical values. However, efficient biosynthesis of MIAs in heterologous microorganisms is hindered due to low supply of key precursors such as geraniol and its derivative 8-hydroxygeraniol catalyzed by geraniol 8-hydroxylase (G8H). In this study, we developed a facile evolution platform to screen strains with improved yield of geraniol by using the SCRaMbLE system embedded in the Sc2.0 synthetic yeast and confirmed the causal role of relevant genomic targets. Through genome mining, we identified several G8H enzymes that perform much better than the commonly used CrG8H for 8-hydroxygeraniol production in vivo. We further showed that the N-terminus of these G8H enzymes plays an important role in cellular activity by swapping experiments. Finally, the combination of the engineered chassis, optimized biosynthesis pathway, and utilization of G8H led to the final strain with more than 30-fold improvement in producing 8-hydroxygeraniol compared with the starting strain. Overall, this study will provide insights into the construction and optimization of yeast cells for efficient biosynthesis of 8-hydroxygeraniol and its derivatives.

A phylogeny of the evening primrose family (Onagraceae) using a target enrichment approach with 303 nuclear loci

BACKGROUND

The evening primrose family (Onagraceae) includes 664 species (803 taxa) with a center of diversity in the Americas, especially western North America. Ongoing research in Onagraceae includes exploring striking variation in floral morphology, scent composition, and breeding system, as well as the role of these traits in driving diversity among plants and their interacting pollinators and herbivores. However, these efforts are limited by the lack of a comprehensive, well-resolved phylogeny. Previous phylogenetic studies based on a few loci strongly support the monophyly of the family and the sister relationship of the two largest tribes but fail to resolve several key relationships.

RESULTS

We used a target enrichment approach to reconstruct the phylogeny of Onagraceae using 303 highly conserved, low-copy nuclear loci. We present a phylogeny for Onagraceae with 169 individuals representing 152 taxa sampled across the family, including extensive sampling within the largest tribe, Onagreae. Deep splits within the family are strongly supported, whereas relationships among closely related genera and species are characterized by extensive conflict among individual gene trees.

CONCLUSIONS

This phylogenetic resource will augment current research projects focused throughout the family in genomics, ecology, coevolutionary dynamics, biogeography, and the evolution of characters driving diversification in the family.

Phylogenomics sheds new light on the drivers behind a long-lasting systematic riddle: the figwort family Scrophulariaceae

The figwort family, Scrophulariaceae, comprises c. 2000 species whose evolutionary relationships at the tribal level have proven difficult to resolve, hindering our ability to understand their origin and diversification. We designed a specific probe kit for Scrophulariaceae, targeting 849 nuclear loci and obtaining plastid regions as by-products. We sampled c. 87% of the genera described in the family and use the nuclear dataset to estimate evolutionary relationships, timing of diversification, and biogeographic patterns. Ten tribes, including two new tribes, Androyeae and Camptolomeae, are supported, and the phylogenetic positions of Androya, Camptoloma, and Phygelius are unveiled. Our study reveals a major diversification at c. 60 million yr ago in some Gondwanan landmasses, where two different lineages diversified, one of which gave rise to nearly 81% of extant species. A Southern African origin is estimated for most modern-day tribes, with two exceptions, the American Leucophylleae, and the mainly Australian Myoporeae. The rapid mid-Eocene diversification is aligned with geographic expansion within southern Africa in most tribes, followed by range expansion to tropical Africa and multiple dispersals out of Africa. Our robust phylogeny provides a framework for future studies aimed at understanding the role of macroevolutionary patterns and processes that generated Scrophulariaceae diversity.

A target enrichment probe set for resolving phylogenetic relationships in the coffee family, Rubiaceae

Premise

Rubiaceae is among the most species-rich plant families, as well as one of the most morphologically and geographically diverse. Currently available phylogenies have mostly relied on few genomic and plastid loci, as opposed to large-scale genomic data. Target enrichment provides the ability to generate sequence data for hundreds to thousands of phylogenetically informative, single-copy loci, which often leads to improved phylogenetic resolution at both shallow and deep taxonomic scales; however, a publicly accessible Rubiaceae-specific probe set that allows for comparable phylogenetic inference across clades is lacking.

Methods

Here, we use publicly accessible genomic resources to identify putatively single-copy nuclear loci for target enrichment in two Rubiaceae groups: tribe Hillieae (Cinchonoideae) and tribal complex Palicoureeae+Psychotrieae (Rubioideae). We sequenced 2270 exonic regions corresponding to 1059 loci in our target clades and generated in silico target enrichment sequences for other Rubiaceae taxa using our designed probe set. To test the utility of our probe set for phylogenetic inference across Rubiaceae, we performed a coalescent-aware phylogenetic analysis using a subset of 27 Rubiaceae taxa from 10 different tribes and three subfamilies, and one outgroup in Apocynaceae.

Results

We recovered an average of 75% and 84% of targeted exons and loci, respectively, per Rubiaceae sample. Probes designed using genomic resources from a particular subfamily were most efficient at targeting sequences from taxa in that subfamily. The number of paralogs recovered during assembly varied for each clade. Phylogenetic inference of Rubiaceae with our target regions resolves relationships at various scales. Relationships are largely consistent with previous studies of relationships in the family with high support (≥0.98 local posterior probability) at nearly all nodes and evidence of gene tree discordance.

Discussion

Our probe set, which we call Rubiaceae2270x, was effective for targeting loci in species across and even outside of Rubiaceae. This probe set will facilitate phylogenomic studies in Rubiaceae and advance systematics and macroevolutionary studies in the family.

Recent advances in triterpenoid pathway elucidation and engineering

Triterpenoids are among the most assorted class of specialized metabolites found in all the taxa of living organisms. Triterpenoids are the leading active ingredients sourced from plant species and are utilized in pharmaceutical and cosmetic industries. The triterpenoid precursor 2,3-oxidosqualene, which is biosynthesized via the mevalonate (MVA) pathway is structurally diversified by the oxidosqualene cyclases (OSCs) and other scaffold-decorating enzymes such as cytochrome P450 monooxygenases (P450s), UDP-glycosyltransferases (UGTs) and acyltransferases (ATs). A majority of the bioactive triterpenoids are harvested from the native hosts using the traditional methods of extraction and occasionally semi-synthesized. These methods of supply are time-consuming and do not often align with sustainability goals. Recent advancements in metabolic engineering and synthetic biology have shown prospects for the green routes of triterpenoid pathway reconstruction in heterologous hosts such as Escherichia coli, Saccharomyces cerevisiae and Nicotiana benthamiana, which appear to be quite promising and might lead to the development of alternative source of triterpenoids. The present review describes the biotechnological strategies used to elucidate complex biosynthetic pathways and to understand their regulation and also discusses how the advances in triterpenoid pathway engineering might aid in the scale-up of triterpenoid production in engineered hosts.

Developing the Protocol Infrastructure for DNA Sequencing Natural History Collections

Intentionally preserved biological material in natural history collections represents a vast repository of biodiversity. Advances in laboratory and sequencing technologies have made these specimens increasingly accessible for genomic analyses, offering a window into the genetic past of species and often permitting access to information that can no longer be sampled in the wild. Due to their age, preparation and storage conditions, DNA retrieved from museum and herbarium specimens is often poor in yield, heavily fragmented and biochemically modified. This not only poses methodological challenges in recovering nucleotide sequences, but also makes such investigations susceptible to environmental and laboratory contamination. In this paper, we review the practical challenges associated with making the recovery of DNA sequence data from museum collections more routine. We first review key operational principles and issues to address, to guide the decision-making process and dialogue between researchers and curators about when and how to sample museum specimens for genomic analyses. We then outline the range of steps that can be taken to reduce the likelihood of contamination including laboratory set-ups, workflows and working practices. We finish by presenting a series of case studies, each focusing on protocol practicalities for the application of different mainstream methodologies to museum specimens including: (i) shotgun sequencing of insect mitogenomes, (ii) whole genome sequencing of insects, (iii) genome skimming to recover plant plastid genomes from herbarium specimens, (iv) target capture of multi-locus nuclear sequences from herbarium specimens, (v) RAD-sequencing of bird specimens and (vi) shotgun sequencing of ancient bovid bone samples.

Two disulfide-reducing pathways are required for the maturation of plastid c-type cytochromes in Chlamydomonas reinhardtii

In plastids, conversion of light energy into ATP relies on cytochrome f, a key electron carrier with a heme covalently attached to a CXXCH motif. Covalent heme attachment requires reduction of the disulfide-bonded CXXCH by CCS5 and CCS4. CCS5 receives electrons from the oxidoreductase CCDA, while CCS4 is a protein of unknown function. In Chlamydomonas reinhardtii, loss of CCS4 or CCS5 yields a partial cytochrome f assembly defect. Here, we report that the ccs4ccs5 double mutant displays a synthetic photosynthetic defect characterized by a complete loss of holocytochrome f assembly. This defect is chemically corrected by reducing agents, confirming the placement of CCS4 and CCS5 in a reducing pathway. CCS4-like proteins occur in the green lineage, and we show that HCF153, a distant ortholog from Arabidopsis thaliana, can substitute for Chlamydomonas CCS4. Dominant suppressor mutations mapping to the CCS4 gene were identified in photosynthetic revertants of the ccs4ccs5 mutants. The suppressor mutations yield changes in the stroma-facing domain of CCS4 that restore holocytochrome f assembly above the residual levels detected in ccs5. Because the CCDA protein accumulation is decreased specifically in the ccs4 mutant, we hypothesize the suppressor mutations enhance the supply of reducing power through CCDA in the absence of CCS5. We discuss the operation of a CCS5-dependent and a CCS5-independent pathway controlling the redox status of the heme-binding cysteines of apocytochrome f.

Channelrhodopsins: From Phototaxis to Optogenetics

Channelrhodopsins stand out among other retinal proteins because of their capacity to generate passive ionic currents following photoactivation. Owing to that, channelrhodopsins are widely used in neuroscience and cardiology as instruments for optogenetic manipulation of the activity of excitable cells. Photocurrents generated by channelrhodopsins were first discovered in the cells of green algae in the 1970s. In this review we describe this discovery and discuss the current state of research in the field.

The D14 and KAI2 Orthologs of Gymnosperms Sense Strigolactones and KL Mimics, Respectively, and the Signals Are Transduced to Control Downstream Genes

Strigolactones (SLs), lactone-containing carotenoid derivatives, function as signaling molecules in the rhizosphere, inducing symbiosis with arbuscular mycorrhizal. In addition, as a class of plant hormones, SLs control plant growth and development in flowering plants (angiosperms). Recent studies show that the ancestral function of SLs, which precede terrestrialization of plants, is as rhizosphere signaling molecules. SLs were then recruited as a class of plant hormones through the step-by-step acquisition of signaling components. The D14 gene encoding the SL receptor arose by gene duplication of KARRIKIN INSENSITIVE2 (KAI2), the receptor of karrikins and KAI2 ligand (KL), an unknown ligand, in the common ancestor of seed plants. KL signaling targets SMAX1, a repressor protein. On the other hand, the SL signaling targets SMXL78 subclade repressors, which arose by duplication of SMAX1 in angiosperms. Thus, gymnosperms contain the SL receptor D14 but not SMXL78, the SL signaling-specific repressor proteins. We studied two gymnosperm species, ginkgo (Ginkgo biloba) and Japanese umbrella pine (Sciadopitys verticillata), to clarify whether SLs are perceived and the signals are transduced in gymnosperms. We show that D14 and KAI2 of ginkgo and Japanese umbrella pine specifically perceive an SL analog and KL mimic, respectively. Furthermore, our results suggest that both SL signaling and KL signaling target SMAX1, and the specific localization of the receptor may result in the specificity of the signaling in gymnosperms.

Extreme genome diversity and cryptic speciation in a harmful algal-bloom-forming eukaryote

Harmful algal blooms of the toxic haptophyte Prymnesium parvum are a recurrent problem in many inland and estuarine waters around the world. Strains of P. parvum vary in the toxins they produce and in other physiological traits associated with harmful algal blooms, but the genetic basis for this variation is unknown. To investigate genome diversity in this morphospecies, we generated genome assemblies for 15 phylogenetically and geographically diverse strains of P. parvum, including Hi-C guided, near-chromosome-level assemblies for two strains. Comparative analysis revealed considerable DNA content variation between strains, ranging from 115 to 845 Mbp. Strains included haploids, diploids, and polyploids, but not all differences in DNA content were due to variation in genome copy number. Haploid genome size between strains of different chemotypes differed by as much as 243 Mbp. Syntenic and phylogenetic analyses indicate that UTEX 2797, a common laboratory strain from Texas, is a hybrid that retains two phylogenetically distinct haplotypes. Investigation of gene families variably present across the strains identified several functional categories associated with metabolic and genome size variation in P. parvum, including genes for the biosynthesis of toxic metabolites and proliferation of transposable elements. Together, our results indicate that P. parvum comprises multiple cryptic species. These genomes provide a robust phylogenetic and genomic framework for investigations into the eco-physiological consequences of the intra- and inter-specific genetic variation present in P. parvum and demonstrate the need for similar resources for other harmful algal-bloom-forming morphospecies.

Substrate-binding glycine residues are major determinants for hydrolase and ligase activity of plant legumains

Peptide asparaginyl ligases (PALs) are useful tools for precision modifications of proteins and live-cell surfaces by ligating peptides after Asn/Asp (Asx). They share high sequence and structural similarity to plant legumains that are generally known as asparaginyl endopeptidases (AEPs), thus making it challenging to identify PALs from AEPs. In this study, we investigate 875 plant species from algae to seed plants with available sequence data in public databases to identify new PALs. We conducted evolutionary trace analysis on 1500 plant legumains, including eight known PALs, to identify key residues that could differentiate ligases and proteases, followed by recombinant expression and functional validation of 16 novel legumains. Previously, we showed that the substrate-binding sequences flanking the catalytic site can strongly influence the enzymatic direction of a legumain and which we named as ligase-activity determinants (LADs). Here, we show that two conserved substrate-binding Gly residues of LADs are critical, but negative determinants for ligase activity. Our results suggest that specific glycine residues are molecular determinants to identify PALs and AEPs as two different legumain subfamilies, accounting for c. 1% and 88%, respectively.

A nuclear target sequence capture probe set for phylogeny reconstruction of the charismatic plant family Bignoniaceae

The plant family Bignoniaceae is a conspicuous and charismatic element of the tropical flora. The family has a complex taxonomic history, with substantial changes in the classification of the group during the past two centuries. Recent re-classifications at the tribal and generic levels have been largely possible by the availability of molecular phylogenies reconstructed using Sanger sequencing data. However, our complete understanding of the systematics, evolution, and biogeography of the family remains incomplete, especially due to the low resolution and support of different portions of the Bignoniaceae phylogeny. To overcome these limitations and increase the amount of molecular data available for phylogeny reconstruction within this plant family, we developed a bait kit targeting 762 nuclear genes, including 329 genes selected specifically for the Bignoniaceae; 348 genes obtained from the Angiosperms353 with baits designed specifically for the family; and, 85 low-copy genes of known function. On average, 77.4% of the reads mapped to the targets, and 755 genes were obtained per species. After removing genes with putative paralogs, 677 loci were used for phylogenetic analyses. On-target genes were compared and combined in the Exon-Only dataset, and on-target + off-target regions were combined in the Supercontig dataset. We tested the performance of the bait kit at different taxonomic levels, from family to species-level, using 38 specimens of 36 different species of Bignoniaceae, representing: 1) six (out of eight) tribal level-clades (e.g., Bignonieae, Oroxyleae, Tabebuia Alliance, Paleotropical Clade, Tecomeae, and Jacarandeae), only Tourrettieae and Catalpeae were not sampled; 2) all 20 genera of Bignonieae; 3) seven (out of nine) species of Dolichandra (e.g., D. chodatii, D. cynanchoides, D. dentata, D. hispida, D. quadrivalvis, D. uncata, and D. uniguis-cati), only D. steyermarkii and D. unguiculata were not sampled; and 4) three individuals of Dolichandra unguis-cati. Our data reconstructed a well-supported phylogeny of the Bignoniaceae at different taxonomic scales, opening new perspectives for a comprehensive phylogenetic framework for the family as a whole.

Genomic capacities for Reactive Oxygen Species metabolism across marine phytoplankton

Marine phytoplankton produce and scavenge Reactive Oxygen Species, to support cellular processes, while limiting damaging reactions. Some prokaryotic picophytoplankton have, however, lost all genes encoding scavenging of hydrogen peroxide. Such losses of metabolic function can only apply to Reactive Oxygen Species which potentially traverse the cell membrane outwards, before provoking damaging intracellular reactions. We hypothesized that cell radius influences which elements of Reactive Oxygen Species metabolism are partially or fully dispensable from a cell. We therefore investigated genomes and transcriptomes from diverse marine eukaryotic phytoplankton, ranging from 0.4 to 44 μm radius, to analyze the genomic allocations encoding enzymes metabolizing Reactive Oxygen Species. Superoxide has high reactivity, short lifetimes and limited membrane permeability. Genes encoding superoxide scavenging are ubiquitous across phytoplankton, but the fractional gene allocation decreased with increasing cell radius, consistent with a nearly fixed set of core genes for scavenging superoxide pools. Hydrogen peroxide has lower reactivity, longer intracellular and extracellular lifetimes and readily crosses cell membranes. Genomic allocations to both hydrogen peroxide production and scavenging decrease with increasing cell radius. Nitric Oxide has low reactivity, long intracellular and extracellular lifetimes and readily crosses cell membranes. Neither Nitric Oxide production nor scavenging genomic allocations changed with increasing cell radius. Many taxa, however, lack the genomic capacity for nitric oxide production or scavenging. The probability of presence of capacity to produce nitric oxide decreases with increasing cell size, and is influenced by flagella and colony formation. In contrast, the probability of presence of capacity to scavenge nitric oxide increases with increasing cell size, and is again influenced by flagella and colony formation.

Using landscape genomics to assess local adaptation and genomic vulnerability of a perennial herb Tetrastigma hemsleyanum (Vitaceae) in subtropical China

Understanding adaptive genetic variation of plant populations and their vulnerabilities to climate change are critical to preserve biodiversity and subsequent management interventions. To this end, landscape genomics may represent a cost-efficient approach for investigating molecular signatures underlying local adaptation. Tetrastigma hemsleyanum is, in its native habitat, a widespread perennial herb of warm-temperate evergreen forest in subtropical China. Its ecological and medicinal values constitute a significant revenue for local human populations and ecosystem. Using 30,252 single nucleotide polymorphisms (SNPs) derived from reduced-representation genome sequencing in 156 samples from 24 sites, we conducted a landscape genomics study of the T. hemsleyanum to elucidate its genomic variation across multiple climate gradients and genomic vulnerability to future climate change. Multivariate methods identified that climatic variation explained more genomic variation than that of geographical distance, which implied that local adaptation to heterogeneous environment might represent an important source of genomic variation. Among these climate variables, winter precipitation was the strongest predictor of the contemporary genetic structure. F _ST outlier tests and environment association analysis totally identified 275 candidate adaptive SNPs along the genetic and environmental gradients. SNP annotations of these putatively adaptive loci uncovered gene functions associated with modulating flowering time and regulating plant response to abiotic stresses, which have implications for breeding and other special agricultural aims on the basis of these selection signatures. Critically, modelling revealed that the high genomic vulnerability of our focal species via a mismatch between current and future genotype-environment relationships located in central-northern region of the T. hemsleyanum's range, where populations require proactive management efforts such as assistant adaptation to cope with ongoing climate change. Taken together, our results provide robust evidence of local climate adaption for T. hemsleyanum and further deepen our understanding of adaptation basis of herbs in subtropical China.

Phylotranscriptomic Analyses of Mycoheterotrophic Monocots Show a Continuum of Convergent Evolutionary Changes in Expressed Nuclear Genes From Three Independent Nonphotosynthetic Lineages

Mycoheterotrophy is an alternative nutritional strategy whereby plants obtain sugars and other nutrients from soil fungi. Mycoheterotrophy and associated loss of photosynthesis have evolved repeatedly in plants, particularly in monocots. Although reductive evolution of plastomes in mycoheterotrophs is well documented, the dynamics of nuclear genome evolution remains largely unknown. Transcriptome datasets were generated from four mycoheterotrophs in three families (Orchidaceae, Burmanniaceae, Triuridaceae) and related green plants and used for phylogenomic analyses to resolve relationships among the mycoheterotrophs, their relatives, and representatives across the monocots. Phylogenetic trees based on 602 genes were mostly congruent with plastome phylogenies, except for an Asparagales + Liliales clade inferred in the nuclear trees. Reduction and loss of chlorophyll synthesis and photosynthetic gene expression and relaxation of purifying selection on retained genes were progressive, with greater loss in older nonphotosynthetic lineages. One hundred seventy-four of 1375 plant benchmark universally conserved orthologous genes were undetected in any mycoheterotroph transcriptome or the genome of the mycoheterotrophic orchid Gastrodia but were expressed in green relatives, providing evidence for massively convergent gene loss in nonphotosynthetic lineages. We designate this set of deleted or undetected genes Missing in Mycoheterotrophs (MIM). MIM genes encode not only mainly photosynthetic or plastid membrane proteins but also a diverse set of plastid processes, genes of unknown function, mitochondrial, and cellular processes. Transcription of a photosystem II gene (psb29) in all lineages implies a nonphotosynthetic function for this and other genes retained in mycoheterotrophs. Nonphotosynthetic plants enable novel insights into gene function as well as gene expression shifts, gene loss, and convergence in nuclear genomes.

Alkaloid production and response to natural adverse conditions in Peganum harmala: in silico transcriptome analyses

Peganum harmala is a valuable wild plant that grows and survives under adverse conditions and produces pharmaceutical alkaloid metabolites. Using different assemblers to develop a transcriptome improves the quality of assembled transcriptome. In this study, a concrete and accurate method for detecting stress-responsive transcripts by comparing stress-related gene ontology (GO) terms and public domains was designed. An integrated transcriptome for P. harmala including 42 656 coding sequences was created by merging de novo assembled transcriptomes. Around 35 000 transcripts were annotated with more than 90% resemblance to three closely related species of Citrus, which confirmed the robustness of the assembled transcriptome; 4853 stress-responsive transcripts were identified. CYP82 involved in alkaloid biosynthesis showed a higher number of transcripts in P. harmala than in other plants, indicating its diverse alkaloid biosynthesis attributes. Transcription factors (TFs) and regulatory elements with 3887 transcripts comprised 9% of the transcriptome. Among the TFs of the integrated transcriptome, cystein2/histidine2 (C2H2) and WD40 repeat families were the most abundant. The Kyoto Encyclopedia of Genes and Genomes (KEGG) MAPK (mitogen-activated protein kinase) signaling map and the plant hormone signal transduction map showed the highest assigned genes to these pathways, suggesting their potential stress resistance. The P. harmala whole-transcriptome survey provides important resources and paves the way for functional and comparative genomic studies on this plant to discover stress-tolerance-related markers and response mechanisms in stress physiology, phytochemistry, ecology, biodiversity, and evolution. P. harmala can be a potential model for studying adverse environmental cues and metabolite biosynthesis and a major source for the production of various alkaloids.

JAZ is essential for ligand specificity of the COI1/JAZ co-receptor

Jasmonates are phytohormones that regulate defense and developmental processes in land plants. Despite the chemical diversity of jasmonate ligands in different plant lineages, they are all perceived by COI1/JAZ co-receptor complexes, in which the hormone acts as a molecular glue between the COI1 F-box and a JAZ repressor. It has been shown that COI1 determines ligand specificity based on the receptor crystal structure and the identification of a single COI1 residue, which is responsible for the evolutionary switch in ligand binding. In this work, we show that JAZ proteins contribute to ligand specificity together with COI1. We propose that specific features of JAZ proteins, which are conserved in bryophytes and lycophytes, enable perception of dn-OPDA ligands regardless the size of the COI1 binding pocket. In vascular plant lineages beyond lycophytes, JAZ evolved to limit binding to JA-Ile, thus impeding dn-OPDA recognition by COI1.

Phylogenomic resolution of order- and family-level monocot relationships using 602 single-copy nuclear genes and 1375 BUSCO genes

We assess relationships among 192 species in all 12 monocot orders and 72 of 77 families, using 602 conserved single-copy (CSC) genes and 1375 benchmarking single-copy ortholog (BUSCO) genes extracted from genomic and transcriptomic datasets. Phylogenomic inferences based on these data, using both coalescent-based and supermatrix analyses, are largely congruent with the most comprehensive plastome-based analysis, and nuclear-gene phylogenomic analyses with less comprehensive taxon sampling. The strongest discordance between the plastome and nuclear gene analyses is the monophyly of a clade comprising Asparagales and Liliales in our nuclear gene analyses, versus the placement of Asparagales and Liliales as successive sister clades to the commelinids in the plastome tree. Within orders, around six of 72 families shifted positions relative to the recent plastome analysis, but four of these involve poorly supported inferred relationships in the plastome-based tree. In Poales, the nuclear data place a clade comprising Ecdeiocoleaceae+Joinvilleaceae as sister to the grasses (Poaceae); Typhaceae, (rather than Bromeliaceae) are resolved as sister to all other Poales. In Commelinales, nuclear data place Philydraceae sister to all other families rather than to a clade comprising Haemodoraceae+Pontederiaceae as seen in the plastome tree. In Liliales, nuclear data place Liliaceae sister to Smilacaceae, and Melanthiaceae are placed sister to all other Liliales except Campynemataceae. Finally, in Alismatales, nuclear data strongly place Tofieldiaceae, rather than Araceae, as sister to all the other families, providing an alternative resolution of what has been the most problematic node to resolve using plastid data, outside of those involving achlorophyllous mycoheterotrophs. As seen in numerous prior studies, the placement of orders Acorales and Alismatales as successive sister lineages to all other extant monocots. Only 21.2% of BUSCO genes were demonstrably single-copy, yet phylogenomic inferences based on BUSCO and CSC genes did not differ, and overall functional annotations of the two sets were very similar. Our analyses also reveal significant gene tree-species tree discordance despite high support values, as expected given incomplete lineage sorting (ILS) related to rapid diversification. Our study advances understanding of monocot relationships and the robustness of phylogenetic inferences based on large numbers of nuclear single-copy genes that can be obtained from transcriptomes and genomes.

Integrated network analyses identify MYB4R1 neofunctionalization in the UV-B adaptation of Tartary buckwheat

A hallmark of adaptive evolution is innovation in gene function, which is associated with the development of distinct roles for genes during plant evolution; however, assessing functional innovation over long periods of time is not trivial. Tartary buckwheat (Fagopyrum tataricum) originated in the Himalayan region and has been exposed to intense UV-B radiation for a long time, making it an ideal species for studying novel UV-B response mechanisms in plants. Here, we developed a workflow to obtain a co-functional network of UV-B responses using data from more than 10,000 samples in more than 80 projects with multi-species and multi-omics data. Dissecting the entire network revealed that flavonoid biosynthesis was most significantly related to the UV-B response. Importantly, we found that the regulatory factor MYB4R1, which resides at the core of the network, has undergone neofunctionalization. In vitro and in vivo experiments demonstrated that MYB4R1 regulates flavonoid and anthocyanin accumulation in response to UV-B in buckwheat by binding to L-box motifs in the FtCHS, FtFLS, and FtUFGT promoters. We used deep learning to develop a visual discrimination model of buckwheat flavonoid content based on natural populations exposed to global UV-B radiation. Our study highlights the critical role of gene neofunctionalization in UV-B adaptation.

A cross-species co-functional gene network underlying leaf senescence

The complex leaf senescence process is governed by various levels of transcriptional and translational regulation. Several features of the leaf senescence process are similar across species, yet the extent to which the molecular mechanisms underlying the process of leaf senescence are conserved remains unclear. Currently used experimental approaches permit the identification of individual pathways that regulate various physiological and biochemical processes; however, the large-scale regulatory network underpinning intricate processes like leaf senescence cannot be built using these methods. Here, we discovered a series of conserved genes involved in leaf senescence in a common horticultural crop (Solanum lycopersicum), a monocot plant (Oryza sativa), and a eudicot plant (Arabidopsis thaliana) through analyses of the evolutionary relationships and expression patterns among genes. Our analyses revealed that the genetic basis of leaf senescence is largely conserved across species. We also created a multi-omics workflow using data from more than 10 000 samples from 85 projects and constructed a leaf senescence-associated co-functional gene network with 2769 conserved, high-confidence functions. Furthermore, we found that the mitochondrial unfolded protein response (UPR^mt) is the central biological process underlying leaf senescence. Specifically, UPR^mt responds to leaf senescence by maintaining mitostasis through a few cross-species conserved transcription factors (e.g. NAC13) and metabolites (e.g. ornithine). The co-functional network built in our study indicates that UPR^mt figures prominently in cross-species conserved mechanisms. Generally, the results of our study provide new insights that will aid future studies of leaf senescence.

Head-to-tail polymerization by VEL proteins underpins cold-induced Polycomb silencing in flowering control

Transcriptional silencing through the Polycomb silencing machinery utilizes a "read-write" mechanism involving histone tail modifications. However, nucleation of silencing and long-term stable transmission of the silenced state also requires P-olycomb Repressive Complex 2 (PRC2) accessory proteins, whose molecular role is poorly understood. The Arabidopsis VEL proteins are accessory proteins that interact with PRC2 to nucleate and propagate silencing at the FLOWERING LOCUS C (FLC) locus, enabling early flowering in spring. Here, we report that VEL proteins contain a domain related to an atypical four-helix bundle that engages in spontaneous concentration-dependent head-to-tail polymerization to assemble dynamic biomolecular condensates. Mutations blocking polymerization of this VEL domain prevent Polycomb silencing at FLC. Plant VEL proteins thus facilitate assembly of dynamic multivalent Polycomb complexes required for inheritance of the silenced state.

Divergent evolutionary trajectories of bryophytes and tracheophytes from a complex common ancestor of land plants

The origin of plants and their colonization of land fundamentally transformed the terrestrial environment. Here we elucidate the basis of this formative episode in Earth history through patterns of lineage, gene and genome evolution. We use new fossil calibrations, a relative clade age calibration (informed by horizontal gene transfer) and new phylogenomic methods for mapping gene family origins. Distinct rooting strategies resolve tracheophytes (vascular plants) and bryophytes (non-vascular plants) as monophyletic sister groups that diverged during the Cambrian, 515-494 million years ago. The embryophyte stem is characterized by a burst of gene innovation, while bryophytes subsequently experienced an equally dramatic episode of reductive genome evolution in which they lost genes associated with the elaboration of vasculature and the stomatal complex. Overall, our analyses reveal that extant tracheophytes and bryophytes are both highly derived from a more complex ancestral land plant. Understanding the origin of land plants requires tracing character evolution across a diversity of modern lineages.

Applications of protein engineering in the microbial synthesis of plant triterpenoids

Triterpenoids are a class of natural products widely used in fields related to medicine and health due to their biological activities such as hepatoprotection, anti-inflammation, anti-viral, and anti-tumor. With the advancement in biotechnology, microorganisms have been used as cell factories to produce diverse natural products. Despite the significant progress that has been made in the construction of microbial cell factories for the heterogeneous biosynthesis of triterpenoids, the industrial production of triterpenoids employing microorganisms has been stymied due to the shortage of efficient enzymes as well as the low expression and low catalytic activity of heterologous proteins in microbes. Protein engineering has been demonstrated as an effective way for improving the specificity, catalytic activity, and stability of the enzyme, which can be employed to overcome these challenges. This review summarizes the current progress in the studies of Oxidosqualene cyclases (OSCs), cytochrome P450s (P450s), and UDP-glycosyltransferases (UGTs), the key enzymes in the triterpenoids synthetic pathway. The main obstacles restricting the efficient catalysis of these key enzymes are analyzed, the applications of protein engineering for the three key enzymes in the microbial synthesis of triterpenoids are systematically reviewed, and the challenges and prospects of protein engineering are also discussed.

Genome sequencing and comparative analysis of Ficus benghalensis and Ficus religiosa species reveal evolutionary mechanisms of longevity

Ficus benghalensis and Ficus religiosa are large woody trees well known for their long lifespan, ecological and traditional significance, and medicinal properties. To understand the genomic and evolutionary aspects of these characteristics, the whole genomes of these Ficus species were sequenced using 10x Genomics linked reads and Oxford Nanopore long reads. The draft genomes of F. benghalensis and F. religiosa comprised of 392.89 Mbp and 332.97 Mbp, respectively. We established the genome-wide phylogenetic positions of the two Ficus species with respect to 50 other Angiosperm species. Comparative evolutionary analyses with other phylogenetically closer Eudicot species revealed adaptive evolution in genes involved in key cellular mechanisms associated with prolonged survival including phytohormones signaling, senescence, disease resistance, and abiotic stress tolerance, which provide genomic insights into the mechanisms conferring longevity and suggest that longevity is a multifaceted phenomenon. This study also provides clues on the existence of CAM pathway in these Ficus species.

•

First whole genome assemblies of Ficus benghalensis and Ficus religiosa trees

•

Genome-wide phylogeny with 50 other Angiosperm species

•

Evolution of genes in phytohormone signaling, senescence, and stress tolerance

•

Genomic insights into longevity-regulating mechanisms

Plant-based engineering for production of high-valued natural products

Covering: up to March 2022Plants are a unique source of complex specialized metabolites, many of which play significant roles in human society. In many cases, however, the availability of these metabolites from naturally occurring sources fails to meet current demands. Thus, there is much interest in expanding the production capacity of target plant molecules. Traditionally, plant breeding, chemical synthesis, and microbial fermentation are considered the primary routes towards large scale production of natural products. Here, we explore the advances, challenges, and future of plant engineering as a complementary path. Although plants are an integral part of our food and agricultural systems and sustain an extensive array of chemical constituents, their complex genetics and physiology have prevented the optimal exploitation of plants as a production chassis. We highlight emerging engineering tools and scientific advances developed in recent years that have improved the prospects of using plants as a sustainable and scalable production platform. We also discuss technological limitations and overall economic outlook of plant-based production of natural products.

An overview of bioinformatics, genomics, and transcriptomics resources for bryophytes

Bryophytes are useful models for the study of plant evolution, development, plant-fungal symbiosis, stress responses, and gametogenesis. Additionally, their dominant haploid gametophytic phase makes them great models for functional genomics research, allowing straightforward genome editing and gene knockout via CRISPR or homologous recombination. Until 2016, however, the only bryophyte genome sequence published was that of Physcomitrium patens. Throughout recent years, several other bryophyte genomes and transcriptome datasets became available, enabling better comparative genomics in evolutionary studies. The increase in the number of bryophyte genome and transcriptome resources available has yielded a plethora of annotations, databases, and bioinformatics tools to access the new data, which covers the large diversity of this clade and whose biology comprises features such as association with arbuscular mycorrhiza fungi, sex chromosomes, low gene redundancy, or loss of RNA editing genes for organellar transcripts. Here we provide a guide to resources available for bryophytes with regards to genome and transcriptome databases and bioinformatics tools.

Diversity, phylogeny, and adaptation of bryophytes: insights from genomic and transcriptomic data

Bryophytes including mosses, liverworts, and hornworts are among the earliest land plants, and occupy a crucial phylogenetic position to aid in the understanding of plant terrestrialization. Despite their small size and simple structure, bryophytes are the second largest group of extant land plants. They live ubiquitously in various habitats and are highly diversified, with adaptive strategies to modern ecosystems on Earth. More and more genomes and transcriptomes have been assembled to address fundamental questions in plant biology. Here, we review recent advances in bryophytes associated with diversity, phylogeny, and ecological adaptation. Phylogenomic studies have provided increasing supports for the monophyly of bryophytes, with hornworts sister to the Setaphyta clade including liverworts and mosses. Further comparative genomic analyses revealed that multiple whole-genome duplications might have contributed to the species richness and morphological diversity in mosses. We highlight that the biological changes through gene gain or neofunctionalization that primarily evolved in bryophytes have facilitated the adaptation to early land environments; among the strategies to adapt to modern ecosystems in bryophytes, desiccation tolerance is the most remarkable. More genomic information for bryophytes would shed light on key mechanisms for the ecological success of these 'dwarfs' in the plant kingdom.

Comparative transcriptomics of tropical woody plants supports fast and furious strategy along the leaf economics spectrum in lianas

Lianas, climbing woody plants, influence the structure and function of tropical forests. Climbing traits have evolved multiple times, including ancestral groups such as gymnosperms and pteridophytes, but the genetic basis of the liana strategy is largely unknown. Here, we use a comparative transcriptomic approach for 47 tropical plant species, including ten lianas of diverse taxonomic origins, to identify genes that are consistently expressed or downregulated only in lianas. Our comparative analysis of full-length transcripts enabled the identification of a core interactomic network common to lianas. Sets of transcripts identified from our analysis reveal features related to functional traits pertinent to leaf economics spectrum in lianas, include upregulation of genes controlling epidermal cuticular properties, cell wall remodeling, carbon concentrating mechanism, cell cycle progression, DNA repair and a large suit of downregulated transcription factors and enzymes involved in ABA-mediated stress response as well as lignin and suberin synthesis. All together, these genes are known to be significant in shaping plant morphologies through responses such as gravitropism, phyllotaxy and shade avoidance.

Summary: The full-length fraction of liana transcriptomes mapped on a protein–protein interactome revealed the nature of their convergence through distinct sets of expressed and downregulated genes not observed in free-standing plants.

Transcriptome Analysis of Soursop (Annona muricata L.) Fruit under Postharvest Storage Identifies Genes Families Involved in Ripening

Soursop (Annona muricata L.) is climacteric fruit with a short ripening period and postharvest shelf life, leading to a rapid softening. In this study, transcriptome analysis of soursop fruits was performed to identify key gene families involved in ripening under postharvest storage conditions (Day 0, Day 3 stored at 28 ± 2 °C, Day 6 at 28 ± 2 °C, Day 3 at 15 ± 2 °C, Day 6 at 15 ± 2 °C, Day 9 at 15 ± 2 °C). The transcriptome analysis showed 224,074 transcripts assembled clustering into 95, 832 unigenes, of which 21, 494 had ORF. RNA-seq analysis showed the highest number of differentially expressed genes on Day 9 at 15 ± 2 °C with 9291 genes (4772 up-regulated and 4519 down-regulated), recording the highest logarithmic fold change in pectin-related genes. Enrichment analysis presented significantly represented GO terms and KEGG pathways associated with molecular function, metabolic process, catalytic activity, biological process terms, as well as biosynthesis of secondary metabolites, plant hormone signal, starch, and sucrose metabolism, plant–pathogen interaction, plant–hormone signal transduction, and MAPK-signaling pathways, among others. Network analysis revealed that pectinesterase genes directly regulate the loss of firmness in fruits stored at 15 ± 2 °C.

Using interdisciplinary, phylogeny-guided approaches to understand the evolution of plant metabolism

To cope with relentless environmental pressures, plants produce an arsenal of structurally diverse chemicals, often called specialized metabolites. These lineage-specific compounds are derived from the simple building blocks made by ubiquitous core metabolic pathways. Although the structures of many specialized metabolites are known, the underlying metabolic pathways and the evolutionary events that have shaped the plant chemical diversity landscape are only beginning to be understood. However, with the advent of multi-omics data sets and the relative ease of studying pathways in previously intractable non-model species, plant specialized metabolic pathways are now being systematically identified. These large datasets also provide a foundation for comparative, phylogeny-guided studies of plant metabolism. Comparisons of metabolic traits and features like chemical abundances, enzyme activities, or gene sequences from phylogenetically diverse plants provide insights into how metabolic pathways evolved. This review highlights the power of studying evolution through the lens of comparative biochemistry, particularly how placing metabolism into a phylogenetic context can help a researcher identify the metabolic innovations enabling the evolution of structurally diverse plant metabolites.

Distribution and the evolutionary history of G-protein components in plant and algal lineages

Heterotrimeric G-protein complexes comprising Gα-, Gβ-, and Gγ-subunits and the regulator of G-protein signaling (RGS) are conserved across most eukaryotic lineages. Signaling pathways mediated by these proteins influence overall growth, development, and physiology. In plants, this protein complex has been characterized primarily from angiosperms with the exception of spreading-leaved earth moss (Physcomitrium patens) and Chara braunii (charophytic algae). Even within angiosperms, specific G-protein components are missing in certain species, whereas unique plant-specific variants-the extra-large Gα (XLGα) and the cysteine-rich Gγ proteins-also exist. The distribution and evolutionary history of G-proteins and their function in nonangiosperm lineages remain mostly unknown. We explored this using the wealth of available sequence data spanning algae to angiosperms representing extant species that diverged approximately 1,500 million years ago, using BLAST, synteny analysis, and custom-built Hidden Markov Model profile searches. We show that a minimal set of components forming the XLGαβγ trimer exists in the entire land plant lineage, but their presence is sporadic in algae. Additionally, individual components have distinct evolutionary histories. The XLGα exhibits many lineage-specific gene duplications, whereas Gα and RGS show several instances of gene loss. Similarly, Gβ remained constant in both number and structure, but Gγ diverged before the emergence of land plants and underwent changes in protein domains, which led to three distinct subtypes. These results highlight the evolutionary oddities and summarize the phyletic patterns of this conserved signaling pathway in plants. They also provide a framework to formulate pertinent questions on plant G-protein signaling within an evolutionary context.

Biotechnological Approaches to Optimize the Production of Amaryllidaceae Alkaloids

Amaryllidaceae alkaloids (AAs) are plant specialized metabolites with therapeutic properties exclusively produced by the Amaryllidaceae plant family. The two most studied representatives of the family are galanthamine, an acetylcholinesterase inhibitor used as a treatment of Alzheimer’s disease, and lycorine, displaying potent in vitro and in vivo cytotoxic and antiviral properties. Unfortunately, the variable level of AAs’ production in planta restricts most of the pharmaceutical applications. Several biotechnological alternatives, such as in vitro culture or synthetic biology, are being developed to enhance the production and fulfil the increasing demand for these AAs plant-derived drugs. In this review, current biotechnological approaches to produce different types of bioactive AAs are discussed.

A phylotranscriptome study using silica gel-dried leaf tissues produces an updated robust phylogeny of Ranunculaceae

The utility of transcriptome data in plant phylogenetics has gained popularity in recent years. However, because RNA degrades much more easily than DNA, the logistics of obtaining fresh tissues has become a major limiting factor for widely applying this method. Here, we used Ranunculaceae to test whether silica-dried plant tissues could be used for RNA extraction and subsequent phylogenomic studies. We sequenced 27 transcriptomes, 21 from silica gel-dried (SD-samples) and six from liquid nitrogen-preserved (LN-samples) leaf tissues, and downloaded 27 additional transcriptomes from GenBank. Our results showed that although the LN-samples produced slightly better reads than the SD-samples, there were no significant differences in RNA quality and quantity, assembled contig lengths and numbers, and BUSCO comparisons between two treatments. Using these data, we conducted phylogenomic analyses, including concatenated- and coalescent-based phylogenetic reconstruction, molecular dating, coalescent simulation, phylogenetic network estimation, and whole genome duplication (WGD) inference. The resulting phylogeny was consistent with previous studies with higher resolution and statistical support. The 11 core Ranunculaceae tribes grouped into two chromosome type clades (T- and R-types), with high support. Discordance among gene trees is likely due to hybridization and introgression, ancient genetic polymorphism and incomplete lineage sorting. Our results strongly support one ancient hybridization event within the R-type clade and three WGD events in Ranunculales. Evolution of the three Ranunculaceae chromosome types is likely not directly related to WGD events. By clearly resolving the Ranunculaceae phylogeny, we demonstrated that SD-samples can be used for RNA-seq and phylotranscriptomic studies of angiosperms.

Reticulate evolution in the Pteris fauriei group (Pteridaceae)

The Pteris fauriei group (Pteridaceae) has a wide distribution in Eastern Asia and includes 18 species with similar but varied morphology. We collected more than 300 specimens of the P. fauriei group and determined ploidy by flow cytometry and inferred phylogenies by molecular analyses of chloroplast and nuclear DNA markers. Our results reveal a complicated reticulate evolution, consisting of seven parental taxa and 58 hybrids. The large number of hybrid taxa have added significant morphological complexity to the group leading to difficult taxonomic issues. The hybrids generally had broader ranges and more populations than their parental taxa. Genetic combination of different pairs of parental species created divergent phenotypes of hybrids, exhibited by both morphological characteristics and ecological fidelities. Niche novelty could facilitate hybrid speciation. Apogamy is common in this group and potentially contributes to the sustainability of the whole group. We propose that frequent hybridizations among members of the P. fauriei group generate and maintain genetic diversity, via novel genetic combinations, niche differentiation, and apogamy.

Hololectin Interdomain Linker Determines Asparaginyl Endopeptidase-Mediated Maturation of Antifungal Hevein-Like Peptides in Oats

Heveins and hevein-containing (hev-) lectins play important roles in stress and pathogenic responses in plants but cause health concerns in humans. Hev-hololectins contain multiple modular hev-peptide domains and are abundantly present in cereals and pseudocereals. However, it is unclear why some cereal hev-hololectins are presented as different forms of proteolytically processed proteoforms. Here we show the precursor architectures of hev-hololectins lead to different processing mechanisms to give either hololectins or hevein-like peptides. We used mass spectrometry and datamining to screen hev-peptides from common cereals, and identified from the oat plant Avena sativa nine novel hevein-like peptides, avenatide aV1-aV9. Bioinformatic analysis revealed that asparaginyl endopeptidase (AEP) can be responsible for the maturation of the highly homologous avenatides from five oat hev-hololectin precursors, each containing four tandemly repeating, hev-like avenatide domains connected by AEP-susceptible linkers with 13-16 residues in length. Further analysis of cereal hev-hololectins showed that the linker lengths provide a distinguishing feature between their cleavable and non-cleavable precursors, with the cleavables having considerably longer linkers (>13 amino acids) than the non-cleavables (<6 amino acids). A detailed study of avenatide aV1 revealed that it contains eight cysteine residues which form a structurally compact, metabolic-resistant cystine-knotted framework with a well-defined chitin-binding site. Antimicrobial assays showed that avenatide aV1 is anti-fungal and inhibits the growth of phyto-pathogenic fungi. Together, our findings of cleavable and non-cleavable hololectins found in cereals expand our knowledge to their biosynthesis and provide insights for hololectin-related health concerns in human.

CHD chromatin remodeling protein diversification yields novel clades and domains absent in classic model organisms

Chromatin remodelers play a fundamental role in the assembly of chromatin, regulation of transcription, and DNA repair. Biochemical and functional characterizations of the CHD family of chromatin remodelers from a variety of model organisms have shown that these remodelers participate in a wide range of activities. However, because the evolutionary history of CHD homologs is unclear, it is difficult to predict which of these activities are broadly conserved and which have evolved more recently in individual eukaryotic lineages. Here, we performed a comprehensive phylogenetic analysis of 8,042 CHD homologs from 1,894 species to create a model for the evolution of this family across eukaryotes with a particular focus on the timing of duplications that gave rise to the diverse copies observed in plants, animals, and fungi. Our analysis confirms that the three major subfamilies of CHD remodelers originated in the eukaryotic last common ancestor, and subsequent losses occurred independently in different lineages. Improved taxon sampling identified several subfamilies of CHD remodelers in plants that were absent or highly divergent in the model plant Arabidopsis thaliana. Whereas the timing of CHD subfamily expansions in vertebrates corresponds to whole genome duplication events, the mechanisms underlying CHD diversification in land plants appear more complicated. Analysis of protein domains reveals that CHD remodeler diversification has been accompanied by distinct transitions in domain architecture, contributing to the functional differences observed between these remodelers. This study demonstrates the importance of proper taxon sampling when studying ancient evolutionary events to prevent misinterpretation of subsequent lineage-specific changes and provides an evolutionary framework for functional and comparative analysis of this critical chromatin remodeler family across eukaryotes.

Gene-Guided Discovery and Ribosomal Biosynthesis of Moroidin Peptides

Moroidin is a bicyclic plant octapeptide with tryptophan side-chain cross-links, originally isolated as a pain-causing agent from the Australian stinging tree Dendrocnide moroides. Moroidin and its analog celogentin C, derived from Celosia argentea, are inhibitors of tubulin polymerization and, thus, lead structures for cancer therapy. However, low isolation yields from source plants and challenging organic synthesis hinder moroidin-based drug development. Here, we present biosynthesis as an alternative route to moroidin-type bicyclic peptides and report that they are ribosomally synthesized and posttranslationally modified peptides (RiPPs) derived from BURP-domain peptide cyclases in plants. By mining 793 plant transcriptomes for moroidin core peptide motifs within BURP-domain precursor peptides, we identified a moroidin cyclase in Japanese kerria, which catalyzes the installation of the tryptophan-indole-centered macrocyclic bonds of the moroidin bicyclic motif in the presence of cupric ions. Based on the kerria moroidin cyclase, we demonstrate the feasibility of producing diverse moroidins including celogentin C in transgenic tobacco plants and report specific cytotoxicity of celogentin C against a lung adenocarcinoma cancer cell line. Our study sets the stage for future biosynthetic development of moroidin-based therapeutics and highlights that mining plant transcriptomes can reveal bioactive cyclic peptides and their underlying cyclases from new source plants.

CHAPAO: Likelihood and hierarchical reference-based representation of biomolecular sequences and applications to compressing multiple sequence alignments

Background

High-throughput experimental technologies are generating tremendous amounts of genomic data, offering valuable resources to answer important questions and extract biological insights. Storing this sheer amount of genomic data has become a major concern in bioinformatics. General purpose compression techniques (e.g. gzip, bzip2, 7-zip) are being widely used due to their pervasiveness and relatively good speed. However, they are not customized for genomic data and may fail to leverage special characteristics and redundancy of the biomolecular sequences.

Results

We present a new lossless compression method CHAPAO (COmpressing Alignments using Hierarchical and Probabilistic Approach), which is especially designed for multiple sequence alignments (MSAs) of biomolecular data and offers very good compression gain. We have introduced a novel hierarchical referencing technique to represent biomolecular sequences which combines likelihood based analyses of the sequence similarities and graph theoretic algorithms. We performed an extensive evaluation study using a collection of real biological data from the avian phylogenomics project, 1000 plants project (1KP), and 16S and 23S rRNA datasets. We report the performance of CHAPAO in comparison with general purpose compression techniques as well as with MFCompress and Nucleotide Archival Format (NAF)—two of the best known methods especially designed for FASTA files. Experimental results suggest that CHAPAO offers significant improvements in compression gain over most other alternative methods. CHAPAO is freely available as an open source software at https://github.com/ashiq24/CHAPAO.

Conclusion

CHAPAO advances the state-of-the-art in compression algorithms and represents a potential alternative to the general purpose compression techniques as well as to the existing specialized compression techniques for biomolecular sequences.

Phylotranscriptomics of Theaceae: generic-level relationships, reticulation and whole-genome duplication

BACKGROUND AND AIMS

Theaceae, with three tribes, nine genera and more than 200 species, are of great economic and ecological importance. Recent phylogenetic analyses based on plastomic data resolved the relationships among the three tribes and the intergeneric relationships within two of those tribes. However, generic-level relationships within the largest tribe, Theeae, were not fully resolved. The role of putative whole-genome duplication (WGD) events in the family and possible hybridization events among genera within Theeae also remain to be tested further.

METHODS

Transcriptomes or low-depth whole-genome sequencing of 57 species of Theaceae, as well as additional plastome sequence data, were generated. Using a dataset of low-copy nuclear genes, we reconstructed phylogenetic relationships using concatenated, species tree and phylogenetic network approaches. We further conducted molecular dating analyses and inferred possible WGD events by examining the distribution of the number of synonymous substitutions per synonymous site (Ks) for paralogues in each species. For plastid protein-coding sequences , phylogenies were reconstructed for comparison with the results obtained from analysis of the nuclear dataset.

RESULTS

Based on the 610 low-copy nuclear genes (858 606 bp in length) investigated, Stewartieae was resolved as sister to the other two tribes. Within Theeae, the Apterosperma-Laplacea clade grouped with Pyrenaria, leaving Camellia and Polyspora as sister. The estimated ages within Theaceae were largely consistent with previous studies based mainly on plastome data. Two reticulation events within Camellia and one between the common ancestor of Gordonia and Schima were found. All members of the tea family shared two WGD events, an older At-γ and a recent Ad-β; both events were also shared with the outgroups (Diapensiaceae, Pentaphylacaceae, Styracaceae and Symplocaceae).

CONCLUSIONS

Our analyses using low-copy nuclear genes improved understanding of phylogenetic relationships at the tribal and generic levels previously proposed based on plastome data, but the phylogenetic position of the Apterosperma-Laplacea clade needs more attention. There is no evidence for extensive intergeneric hybridization within Theeae or for a Theaceae-specific WGD event. Land bridges (e.g. the Bering land bridge) during the Late Oligocene may have permitted the intercontinental plant movements that facilitated the putative ancient introgression between the common ancestor of Gordonia and Schima.

Heterologous expression of Bixa orellana cleavage dioxygenase 4-3 drives crocin but not bixin biosynthesis

Annatto (Bixa orellana) is a perennial shrub native to the Americas, and bixin, derived from its seeds, is a methoxylated apocarotenoid used as a food and cosmetic colorant. Two previous reports claimed to have isolated the carotenoid cleavage dioxygenase (CCD) responsible for the production of the putative precursor of bixin, the C24 apocarotenal bixin dialdehyde. We re-assessed the activity of six Bixa CCDs and found that none of them produced substantial amounts of bixin dialdehyde in Escherichia coli. Unexpectedly, BoCCD4-3 cleaved different carotenoids (lycopene, β-carotene, and zeaxanthin) to yield the C20 apocarotenal crocetin dialdehyde, the known precursor of crocins, which are glycosylated apocarotenoids accumulated in saffron stigmas. BoCCD4-3 lacks a recognizable transit peptide but localized to plastids, the main site of carotenoid accumulation in plant cells. Expression of BoCCD4-3 in Nicotiana benthamiana leaves (transient expression), tobacco (Nicotiana tabacum) leaves (chloroplast transformation, under the control of a synthetic riboswitch), and in conjunction with a saffron crocetin glycosyl transferase, in tomato (Solanum lycopersicum) fruits (nuclear transformation) led to high levels of crocin accumulation, reaching the highest levels (>100 µg/g dry weight) in tomato fruits, which also showed a crocin profile similar to that found in saffron, with highly glycosylated crocins as major compounds. Thus, while the bixin biosynthesis pathway remains unresolved, BoCCD4-3 can be used for the metabolic engineering of crocins in a wide range of different plant tissues.