The Madagascar palm genome provides new insights on the evolution of Apocynaceae specialized metabolism.

Specialized metabolites possess diverse interesting biological activities and some cardenolides- and monoterpene indole alkaloids- (MIAs) derived pharmaceuticals are currently used to treat human diseases such as cancers or hypertension. While these two families of biocompounds are produced by specific subfamilies of Apocynaceae, one member of this medicinal plant family, the succulent tree Pachypodium lamerei Drake (also known as Madagascar palm), does not produce such specialized metabolites. To explore the evolutionary paths that have led to the emergence and loss of cardenolide and MIA biosynthesis in Apocynaceae, we sequenced and assembled the P. lamerei genome by combining Oxford Nanopore Technologies long-reads and Illumina short-reads. Phylogenomics revealed that, among the Apocynaceae whose genomes have been sequenced, the Madagascar palm is so far the species closest to the common ancestor between MIA producers/non-MIA producers. Transposable elements, constituting 72.48% of the genome, emerge as potential key players in shaping genomic architecture and influencing specialized metabolic pathways. The absence of crucial MIA biosynthetic genes such as strictosidine synthase in P. lamerei and non-Rauvolfioideae species hints at a transposon-mediated mechanism behind gene loss. Phylogenetic analysis not only showcases the evolutionary divergence of specialized metabolite biosynthesis within Apocynaceae but also underscores the role of transposable elements in this intricate process. Moreover, we shed light on the low conservation of enzymes involved in the final stages of MIA biosynthesis in the distinct MIA-producing plant families, inferring independent gains of these specialized enzymes along the evolution of these medicinal plant clades. Overall, this study marks a leap forward in understanding the genomic dynamics underpinning the evolution of specialized metabolites biosynthesis in the Apocynaceae family, with transposons emerging as potential architects of genomics restructuring and gene loss.

Coupling High-Throughput and Targeted Screening for Identification of Nonobvious Metabolic Engineering Targets.

Identification of metabolic engineering targets is a fundamental challenge in strain development programs. While high-throughput (HTP) genetic engineering methodologies capable of generating vast diversity are being developed at a rapid rate, a majority of industrially interesting molecules cannot be screened at sufficient throughput to leverage these techniques. We propose a workflow that couples HTP screening of common precursors (e.g., amino acids) that can be screened either directly or by artificial biosensors, with low-throughput targeted validation of the molecule of interest to uncover nonintuitive beneficial metabolic engineering targets and combinations hereof. Using this workflow, we identified several nonobvious novel targets for improving p-coumaric acid (p-CA) and l-DOPA production from two large 4k gRNA libraries each deregulating 1000 metabolic genes in the yeast Saccharomyces cerevisiae. We initially screened yeast cells transformed with gRNA library plasmids for individual regulatory targets improving the production of l-tyrosine-derived betaxanthins, identifying 30 targets that increased intracellular betaxanthin content 3.5-5.7 fold. Hereafter, we screened the targets individually in a high-producing p-CA strain, narrowing down the targets to six that increased the secreted titer by up to 15%. To investigate whether any of the six targets could be additively combined to improve p-CA production further, we created a gRNA multiplexing library and subjected it to our proposed coupled workflow. The combination of regulating PYC1 and NTH2 simultaneously resulted in the highest (threefold) improvement of the betaxanthin content, and an additive trend was also observed in the p-CA strain. Lastly, we tested the initial 30 targets in a l-DOPA producing strain, identifying 10 targets that increased the secreted titer by up to 89%, further validating our screening by proxy workflow. This coupled approach is useful for strain development in the absence of direct HTP screening assays for products of interest.

The Rauvolfia tetraphylla genome suggests multiple distinct biosynthetic routes for yohimbane monoterpene indole alkaloids.

Monoterpene indole alkaloids (MIAs) are a structurally diverse family of specialized metabolites mainly produced in Gentianales to cope with environmental challenges. Due to their pharmacological properties, the biosynthetic modalities of several MIA types have been elucidated but not that of the yohimbanes. Here, we combine metabolomics, proteomics, transcriptomics and genome sequencing of Rauvolfia tetraphylla with machine learning to discover the unexpected multiple actors of this natural product synthesis. We identify a medium chain dehydrogenase/reductase (MDR) that produces a mixture of four diastereomers of yohimbanes including the well-known yohimbine and rauwolscine. In addition to this multifunctional yohimbane synthase (YOS), an MDR synthesizing mainly heteroyohimbanes and the short chain dehydrogenase vitrosamine synthase also display a yohimbane synthase side activity. Lastly, we establish that the combination of geissoschizine synthase with at least three other MDRs also produces a yohimbane mixture thus shedding light on the complex mechanisms evolved for the synthesis of these plant bioactives.

The Vinca minor genome highlights conserved evolutionary traits in monoterpene indole alkaloid synthesis.

Vinca minor, also known as the lesser periwinkle, is a well-known species from the Apocynaceae, native to central and southern Europe. This plant synthesizes monoterpene indole alkaloids, which are a class of specialized metabolites displaying a wide range of bioactive- and pharmacologically important properties. Within the almost 50 monoterpene indole alkaloids it produces, V. minor mainly accumulates vincamine, which is commercially used as a nootropic. Using a combination of Oxford Nanopore Technologies long read- and Illumina short-read sequencing, a 679,098 Mb V. minor genome was assembled into 296 scaffolds with an N50 scaffold length of 6 Mb, and encoding 29,624 genes. These genes were functionally annotated and used in a comparative genomic analysis to establish gene families and to investigate gene family expansion and contraction across the phylogenetic tree. Furthermore, homology-based monoterpene indole alkaloid gene predictions together with a metabolic analysis across 4 different V. minor tissue types guided the identification of candidate monoterpene indole alkaloid genes. These candidates were finally used to identify monoterpene indole alkaloid gene clusters, which combined with synteny analysis allowed for the discovery of a functionally validated vincadifformine-16-hydroxylase, reinforcing the potential of this dataset for monoterpene indole alkaloids gene discovery. It is expected that access to these resources will facilitate the elucidation of unknown monoterpene indole alkaloid biosynthetic routes with the potential of transferring these pathways to heterologous expression systems for large-scale monoterpene indole alkaloid production.

Genome Assembly of the Medicinal Plant Voacanga thouarsii.

The Apocynaceae tree Voacanga thouarsii, native to southern Africa and Madagascar, produces monoterpene indole alkaloids (MIA), which are specialized metabolites with a wide range of bioactive properties. Voacanga species mainly accumulates tabersonine in seeds making these species valuable medicinal plants currently used for industrial MIA production. Despite their importance, the MIA biosynthesis in Voacanga species remains poorly studied. Here, we report the first genome assembly and annotation of a Voacanga species. The combined assembly of Oxford Nanopore Technologies long-reads and Illumina short-reads resulted in 3,406 scaffolds with a total length of 1,354.26 Mb and an N50 of 3.04 Mb. A total of 33,300 protein-coding genes were predicted and functionally annotated. These genes were then used to establish gene families and to investigate gene family expansion and contraction across the phylogenetic tree. A transposable element (TE) analysis showed the highest proportion of TE in Voacanga thouarsii compared with all other MIA-producing plants. In a nutshell, this first reference genome of V. thouarsii will thus contribute to strengthen future comparative and evolutionary studies in MIA-producing plants leading to a better understanding of MIA pathway evolution. This will also allow the potential identification of new MIA biosynthetic genes for metabolic engineering purposes.

An updated version of the Madagascar periwinkle genome.

The Madagascar periwinkle, Catharanthus roseus, belongs to the Apocynaceae family. This medicinal plant, endemic to Madagascar, produces many important drugs including the monoterpene indole alkaloids (MIA) vincristine and vinblastine used to treat cancer worldwide. Here, we provide a new version of the C. roseus genome sequence obtained through the combination of Oxford Nanopore Technologies long-reads and Illumina short-reads. This more contiguous assembly consists of 173 scaffolds with a total length of 581.128 Mb and an N50 of 12.241 Mb. Using publicly available RNAseq data, 21,061 protein coding genes were predicted and functionally annotated. A total of 42.87% of the genome was annotated as transposable elements, most of them being long-terminal repeats. Together with the increasing access to MIA-producing plant genomes, this updated version should ease evolutionary studies leading to a better understanding of MIA biosynthetic pathway evolution.

Burden Imposed by Heterologous Protein Production in Two Major Industrial Yeast Cell Factories: Identifying Sources and Mitigation Strategies.

Production of heterologous proteins, especially biopharmaceuticals and industrial enzymes, in living cell factories consumes cellular resources. Such resources are reallocated from normal cellular processes toward production of the heterologous protein that is often of no benefit to the host cell. This competition for resources is a burden to host cells, has a negative impact on cell fitness, and may consequently trigger stress responses. Importantly, this often causes a reduction in final protein titers. Engineering strategies to generate more burden resilient production strains offer sustainable opportunities to increase production and profitability for this growing billion-dollar global industry. We review recently reported impacts of burden derived from resource competition in two commonly used protein-producing yeast cell factories: Saccharomyces cerevisiae and Komagataella phaffii (syn. Pichia pastoris). We dissect possible sources of burden in these organisms, from aspects related to genetic engineering to protein translation and export of soluble protein. We also summarize advances as well as challenges for cell factory design to mitigate burden and increase overall heterologous protein production from metabolic engineering, systems biology, and synthetic biology perspectives. Lastly, future profiling and engineering strategies are highlighted that may lead to constructing robust burden-resistant cell factories. This includes incorporation of systems-level data into mathematical models for rational design and engineering dynamical regulation circuits in production strains.

Transportome-wide engineering of Saccharomyces cerevisiae.

Synthetic biology enables the production of small molecules by recombinant microbes for pharma, food, and materials applications. The secretion of products reduces the cost of separation and purification, but it is challenging to engineer due to the limited understanding of the transporter proteins' functions. Here we describe a method for genome-wide transporter disruption that, in combination with a metabolite biosensor, enables the identification of transporters impacting the production of a given target metabolite in yeast Saccharomyces cerevisiae. We applied the method to study the transport of xenobiotic compounds, cis,cis-muconic acid (CCM), protocatechuic acid (PCA), and betaxanthins. We found 22 transporters that influenced the production of CCM or PCA. The transporter of the 12-spanner drug:H(+) antiporter (DHA1) family Tpo2p was further confirmed to import CCM and PCA in Xenopus expression assays. We also identified three transporter proteins (Qdr1p, Qdr2p, and Apl1p) involved in betaxanthins transport. In summary, the described method enables high-throughput transporter identification for small molecules in cell factories.

Exploring small-scale chemostats to scale up microbial processes: 3-hydroxypropionic acid production in S. cerevisiae.

BACKGROUND: The physiological characterization of microorganisms provides valuable information for bioprocess development. Chemostat cultivations are a powerful tool for this purpose, as they allow defined changes to one single parameter at a time, which is most commonly the growth rate. The subsequent establishment of a steady state then permits constant variables enabling the acquisition of reproducible data sets for comparing microbial performance under different conditions. We performed physiological characterizations of a 3-hydroxypropionic acid (3-HP) producing Saccharomyces cerevisiae strain in a miniaturized and parallelized chemostat cultivation system. The physiological conditions under investigation were various growth rates controlled by different nutrient limitations (C, N, P). Based on the cultivation parameters obtained subsequent fed-batch cultivations were designed. RESULTS: We report technical advancements of a small-scale chemostat cultivation system and its applicability for reliable strain screening under different physiological conditions, i.e. varying dilution rates and different substrate limitations (C, N, P). Exploring the performance of an engineered 3-HP producing S. cerevisiae strain under carbon-limiting conditions revealed the highest 3-HP yields per substrate and biomass of 16.6 %C-mol and 0.43 g gCDW-1, respectively, at the lowest set dilution rate of 0.04 h-1. 3-HP production was further optimized by applying N- and P-limiting conditions, which resulted in a further increase in 3-HP yields revealing values of 21.1 %C-mol and 0.50 g gCDW-1 under phosphate-limiting conditions. The corresponding parameters favoring an increased 3-HP production, i.e. dilution rate as well as C- and P-limiting conditions, were transferred from the small-scale chemostat cultivation system to 1-L bench-top fermenters operating in fed-batch conditions, revealing 3-HP yields of 15.9 %C-mol and 0.45 g gCDW-1 under C-limiting, as well as 25.6 %C-mol and 0.50 g gCDW-1 under phosphate-limiting conditions. CONCLUSIONS: Small-scale chemostat cultures are well suited for the physiological characterization of microorganisms, particularly for investigating the effect of changing cultivation parameters on microbial performance. In our study, optimal conditions for 3-HP production comprised (i) a low dilution rate of 0.04 h-1 under carbon-limiting conditions and (ii) the use of phosphate-limiting conditions. Similar 3-HP yields were achieved in chemostat and fed-batch cultures under both C- and P-limiting conditions proving the growth rate as robust parameter for process transfer and thus the small-scale chemostat system as powerful tool for process optimization.

Matching NLR Immune Receptors to Autoimmunity in camta3 Mutants Using Antimorphic NLR Alleles.

To establish infection, pathogens deploy effectors to modify or remove host proteins. Plant immune receptors with nucleotide-binding, leucine-rich repeat domains (NLRs) detect these modifications and trigger immunity. Plant NLRs thus guard host "guardees." A corollary is that autoimmunity may result from inappropriate NLR activation because mutations in plant guardees could trigger corresponding NLR guards. To explore these hypotheses, we expressed 108 dominant-negative (DN) Arabidopsis NLRs in various lesion mimic mutants, including camta3, which exhibits autoimmunity. CAMTA3 was previously described as a negative regulator of immunity, and we find that autoimmunity in camta3 is fully suppressed by expressing DNs of two NLRs, DSC1 and DSC2. Additionally, expression of either NLR triggers cell death that can be suppressed by CAMTA3 expression. These findings support a model in which DSC1 and DSC2 guard CAMTA3, and they suggest that other negative regulators of immunity may similarly represent guardees.

CrEdit: CRISPR mediated multi-loci gene integration in Saccharomyces cerevisiae.

BACKGROUND: One of the bottlenecks in production of biochemicals and pharmaceuticals in Saccharomyces cerevisiae is stable and homogeneous expression of pathway genes. Integration of genes into the genome of the production organism is often a preferred option when compared to expression from episomal vectors. Existing approaches for achieving stable simultaneous genome integrations of multiple DNA fragments often result in relatively low integration efficiencies and furthermore rely on the use of selection markers. RESULTS: Here, we have developed a novel method, CrEdit (CRISPR/Cas9 mediated genome Editing), which utilizes targeted double strand breaks caused by CRISPR/Cas9 to significantly increase the efficiency of homologous integration in order to edit and manipulate genomic DNA. Using CrEdit, the efficiency and locus specificity of targeted genome integrations reach close to 100% for single gene integration using short homology arms down to 60 base pairs both with and without selection. This enables direct and cost efficient inclusion of homology arms in PCR primers. As a proof of concept, a non-native ß-carotene pathway was reconstructed in S. cerevisiae by simultaneous integration of three pathway genes into individual intergenic genomic sites. Using longer homology arms, we demonstrate highly efficient and locus-specific genome integration even without selection with up to 84% correct clones for simultaneous integration of three gene expression cassettes. CONCLUSIONS: The CrEdit approach enables fast and cost effective genome integration for engineering of S. cerevisiae. Since the choice of the targeting sites is flexible, CrEdit is a powerful tool for diverse genome engineering applications.

Transcriptome and genome size analysis of the Venus flytrap.

The insectivorous Venus flytrap (Dionaea muscipula) is renowned from Darwin's studies of plant carnivory and the origins of species. To provide tools to analyze the evolution and functional genomics of D. muscipula, we sequenced a normalized cDNA library synthesized from mRNA isolated from D. muscipula flowers and traps. Using the Oases transcriptome assembler 79,165,657 quality trimmed reads were assembled into 80,806 cDNA contigs, with an average length of 679 bp and an N50 length of 1,051 bp. A total of 17,047 unique proteins were identified, and assigned to Gene Ontology (GO) and classified into functional categories. A total of 15,547 full-length cDNA sequences were identified, from which open reading frames were detected in 10,941. Comparative GO analyses revealed that D. muscipula is highly represented in molecular functions related to catalytic, antioxidant, and electron carrier activities. Also, using a single copy sequence PCR-based method, we estimated that the genome size of D. muscipula is approx. 3 Gb. Our genome size estimate and transcriptome analyses will contribute to future research on this fascinating, monotypic species and its heterotrophic adaptations.

ATAF1 transcription factor directly regulates abscisic acid biosynthetic gene NCED3 in Arabidopsis thaliana.

ATAF1, an Arabidopsis thaliana NAC transcription factor, plays important roles in plant adaptation to environmental stress and development. To search for ATAF1 target genes, we used protein binding microarrays and chromatin-immunoprecipitation (ChIP). This identified T[A,C,G]CGT[A,G] and TT[A,C,G]CGT as ATAF1 consensus binding sequences. Co-expression analysis across publicly available microarray experiments identified 25 genes co-expressed with ATAF1. The promoter regions of ATAF1 co-expressors were significantly enriched for ATAF1 binding sites, and TTGCGTA was identified in the promoter of the key abscisic acid (ABA) phytohormone biosynthetic gene NCED3. ChIP-qPCR and expression analysis showed that ATAF1 binding to the NCED3 promoter correlated with increased NCED3 expression and ABA hormone levels. These results indicate that ATAF1 regulates ABA biosynthesis.

Structure, function and networks of transcription factors involved in abiotic stress responses.

Transcription factors (TFs) are master regulators of abiotic stress responses in plants. This review focuses on TFs from seven major TF families, known to play functional roles in response to abiotic stresses, including drought, high salinity, high osmolarity, temperature extremes and the phytohormone ABA. Although ectopic expression of several TFs has improved abiotic stress tolerance in plants, fine-tuning of TF expression and protein levels remains a challenge to avoid crop yield loss. To further our understanding of TFs in abiotic stress responses, emerging gene regulatory networks based on TFs and their direct targets genes are presented. These revealed components shared between ABA-dependent and independent signaling as well as abiotic and biotic stress signaling. Protein structure analysis suggested that TFs hubs of large interactomes have extended regions with protein intrinsic disorder (ID), referring to their lack of fixed tertiary structures. ID is now an emerging topic in plant science. Furthermore, the importance of the ubiquitin-proteasome protein degradation systems and modification by sumoylation is also apparent from the interactomes. Therefore; TF interaction partners such as E3 ubiquitin ligases and TF regions with ID represent future targets for engineering improved abiotic stress tolerance in crops.

Order by disorder in plant signaling.

Protein intrinsic disorder (ID), referring to the lack of a fixed tertiary structure, is an emerging topic in plant science. Proteins with ID challenge our perception of protein interactions because of their malleable behavior. They are abundant in highly regulated processes such as cellular signaling and transcription, where they exploit the flexibility of ID. In this opinion article we highlight trends in the field of protein ID and discuss its implications for interactions between plant transcription factors (TFs) and the cellular signaling hub protein RADICAL-INDUCED CELL DEATH 1 (RCD1). We envision RCD1-TF interactions as models for translating knowledge of ID-based interactions in vitro to the organismal level in vivo, and urge increased focus on ID in basic plant research and agricultural sciences.

NAC genes: time-specific regulators of hormonal signaling in Arabidopsis.

Environmental stresses on both animals and plants impose massive transcriptional perturbations. Successful adaptations to such stresses are being orchestrated by both activating and repressing effects of transcription factors on specific target genes. We have recently published a systematic characterization of members of the large NAC gene transcription factor family in the model weed Arabidopsis thaliana. Our analysis revealed interesting sub-groupings of the Arabidopsis NAC genes, relating structure and function. Here we present a meta-analysis revealing distinct temporal expression profiles of NAC genes upon stimuli with seven phytohormones. Our analysis could be a first indication of NAC-centered transcriptional networks, which coordinate timely hormonal signaling in plants.

The HvNAC6 transcription factor: a positive regulator of penetration resistance in barley and Arabidopsis.

Pathogens induce the expression of many genes encoding plant transcription factors, though specific knowledge of the biological function of individual transcription factors remains scarce. NAC transcription factors are encoded in plants by a gene family with proposed functions in both abiotic and biotic stress adaptation, as well as in developmental processes. In this paper, we provide convincing evidence that a barley NAC transcription factor has a direct role in regulating basal defence. The gene transcript was isolated by differential display from barley leaves infected with the biotrophic powdery mildew fungus, Blumeria graminis f.sp. hordei (Bgh). The full-length cDNA clone was obtained using 5'-RACE and termed HvNAC6, due to its high similarity to the rice homologue, OsNAC6. Gene silencing of HvNAC6 during Bgh inoculation compromises penetration resistance in barley epidermal cells towards virulent Bgh. Complementing the effect of HvNAC6 gene silencing, transient overexpression of HvNAC6 increases the occurrence of penetration resistant cells towards Bgh attack. Quantitative RT-PCR shows the early and transient induction of HvNAC6 in barley epidermis upon Bgh infection. Additionally, our results show that the Arabidopsis HvNAC6 homologue ATAF1 is also induced by Bgh and the ataf1-1 mutant line shows decreased penetration resistance to this non-host pathogen. Collectively, these data suggest a conserved role of HvNAC6 and ATAF1 in the regulation of penetration resistance in monocots and dicots, respectively.