A multi-omic Nicotiana benthamiana resource for fundamental research and biotechnology

Broadening the Nicotiana benthamiana research toolbox through the generation of dicer-like mutants using CRISPR/Cas9 approaches

RNA silencing in plants plays a pivotal role in various biological processes, including development, epigenetic modifications and stress response. Key components of this network are Dicer-like (DCL) proteins. Nicotiana benthamiana encodes four DCLs, each responsible for the generation of distinct small RNA (sRNA) populations, which regulate different functions. However, elucidating the precise role of each DCL has been proven challenging, as overlapping functions exist within DCLs. In our present study, we have successfully generated dcl2, dcl3 and dcl4 homozygous mutants, employing two different CRISPR/Cas9 approaches. The first approach is based on a transgene-mediated delivery of the single-guide RNA (sgRNA), while the second approach employs a viral vector for sgRNA delivery. By utilizing a suite of screening techniques, including polymerase chain reaction (PCR), T7 endonuclease I (T7E1) assay, high-resolution melt analysis (HRMA) and DNA sequencing, we successfully generated dcl2, dcl3 and dcl4 homozygous mutants harboring identical mutations in every allele. To evaluate these dcl mutants, we examined their sRNA profiles and phenotypes. We further have indications that homozygous mutations of a gene do not always lead to the desired loss-of-function, highlighting the importance of mutant evaluation. dcl mutants represent invaluable tools to explore how overlapping silencing pathways are connected to essential plant functions, including development, stress responses and pathogen defense. Additionally, they hold potential for biotechnological applications, such as crop improvement and gene silencing tools. We anticipate that our study will make significant contributions to enhance understanding of the role of DCLs in plants.

Nicotiana benthamiana VASCULAR-RELATED NAC-DOMAIN7-2 (NbVND7-2) has a role in xylem formation during interfamily grafting

VASCULAR-RELATED NAC-DOMAIN7 (VND7) is a transcription factor gene that plays a critical role in xylem differentiation. The ectopic expression of VND7 induces the formation of secondary cell walls with spiral patterns in multiple plant cell types. In the present study, we identified four homologs of VND7 in Nicotiana benthamiana and assigned them the names NbVND7-1 to NbVND7-4. In particular, NbVND7-1 and NbVND7-2 were highly expressed during N. benthamiana and Arabidopsis thaliana (Nb/At) interfamily grafting. Analysis of the promoter β-glucuronidase (GUS) reporter lines of NbVND7 genes elucidated the expression of NbVND7-1 and NbVND7-2 in xylem tissues of intact and grafted plants, and those of NbVND7-3 and NbVND7-4 in internal phloem tissues. Gene network analysis revealed the downstream genes of each NbVND7 homolog and highlighted the association of NbVND7-1 and NbVND7-2 with xylem formation. A β-estradiol-inducible system for NbVND7-2 demonstrated that NbVND7-2 promotes ectopic xylem vessel differentiation in N. benthamiana seedlings and in the stem tissues at graft junctions. Induction of NbVND7-2 at graft junctions enhanced ectopic xylem formation in the callus tissues proliferated at the graft boundary, accelerated the initiation of water transport from stock to scion, and enhanced scion stem growth after grafting. This study reveals a role for NbVND7 genes in xylem formation that can enhance Nicotiana interfamily grafting.

In defense of funding foundational plant science

Plants are essential for life as we know it on Earth. They oxygenate the atmosphere, regulate the climate, and comprise much of the primary producers underpinning complex food systems. In the 1980s, a multinational group of plant scientists chose the small angiosperm-Arabidopsis thaliana-to serve as the model flowering plant for genetic and molecular studies that would be leveraged to produce vast new datasets, resources, and tools. The rationale they used to persuade funding agencies to make significant investments and focus intense effort on this single plant species was to produce a deep fundamental knowledge of the biology of plants and to apply this knowledge to valuable, but typically less tractable, plant species. Over the past 40 yr, Arabidopsis has emerged as the most powerful and versatile plant model to uncover core biological principles and served as a prototyping system to test advanced molecular and genetic concepts. We argue that the emerging challenges of accelerating climate instability and a rapidly growing global population call for renewed and robust investments in fundamental plant biology research. Leveraging the power of Arabidopsis research, resources, datasets, and global collaborative community is more important than ever. This commentary lays out a vigorous defense of foundational, i.e. "basic," plant science research; describes that often, Arabidopsis is preferable to working directly in crops; highlights several transformative applications generated from basic plant research; and makes the argument that plant science is vital to the survival of humanity.

Identification of β-galactosidases along the secretory pathway of Nicotiana benthamiana that collectively hamper engineering of galactose-extended glycans on recombinant glycoproteins

Glycosylation is an important aspect for many biopharmaceuticals, including vaccines against parasitic helminths. Plants, especially Nicotiana benthamiana, have proven to be excellent production hosts for biopharmaceuticals with tailor-made glycosylation. If desired, galactosylation can be introduced on biopharmaceuticals through co-expression of the appropriate glycosyltransferase. However, achieving homogenous glycoforms with terminal galactose residues remains difficult as native N. benthamiana β-galactosidases (NbBGALs) truncate these glycans. Recently, the first NbBGAL has been identified, but a knockout line was insufficient to achieve near complete galactosylation, suggesting that other enzymes could have similar activity. In this study, we selected 10 NbBGALs for further investigation into subcellular localization, in vitro and in vivo activity against β1,4-linked galactose on N-glycans and β1,3-linked galactose on O-glycans. We show that NbBGAL3B is localized in the apoplast and has similar specificity for β1,4-linked galactose on N-glycans as the previously identified NbBGAL1. In contrast, none of the selected NbBGALs cleaved β1,3-linked galactose from O-glycans besides BGAL1. In addition, we provide a novel strategy to achieve near complete galactosylation on galactosidase-prone glycoproteins by using the protective capacity of the Lewis X motif and subsequent removal of the antennary fucose residues. Taken together, our results provide a broad view of the ability of NbBGALs to cleave galactoses and have identified NbBGAL3B as the second major contributor of undesired β-galactosidase activity while engineering N-glycans. This work lays the foundation for generating knockout lines that are devoid of undesired NbBGALs and therefore do not hamper the production of recombinant glycoproteins with galactose-extended glycans.

Modulation of key sterol-related genes of Nicotiana benthamiana by phosphite treatment during infection with Phytophthora cinnamomi

Phytophthora cinnamomi is a globally destructive pathogen causing disease in over 5000 plant species. As sterol auxotrophs, Phytophthora species rely on host-derived phytosterols for reproduction, yet the effects of pathogen infection on plant sterol biosynthesis remains unclear. We utilised a soil-free plant growth system to analyze the impacts of P. cinnamomi on Nicotiana benthamiana roots, a new model for studying P. cinnamomi -plant root interactions. Our results show that P. cinnamomi successfully infected all ecotypes tested, but infection was inhibited by the systemic chemical, phosphite. While phosphite is traditionally associated with the activation of plant defence mechanisms, we show that phosphite also modulates plant immune receptors and phytosterol biosynthesis. qPCR analyses revealed a two-fold upregulation of the N. benthamiana elicitin receptor, Responsive to Elicitins (REL ), and its co-receptor, suppressor of BIR1-1 (SOBIR ) during P. cinnamomi infection when compared with infected, phosphite-treated plants. Furthermore, key genes related to plant sterol biosynthesis were upregulated in their expression during pathogen infection but were suppressed in phosphite-treated and infected plants. Notably, the cytochrome P450 family 710 (CYP710A ) gene encoding a C22-sterol desaturase, involved in stigmasterol production, a phytosterol known to be linked to plant susceptibility to pathogens, was downregulated in phosphite-treated plants, independent of infection status. These findings reveal novel insights into the role of phosphite in modulating plant immune responses and sterol metabolism, with potential in managing diseases caused by P. cinnamomi .

Toward an integrated omics approach for plant biosynthetic pathway discovery in the age of AI

Elucidating plant biosynthetic pathways is key to advancing a sustainable bioeconomy by enabling access to complex natural products through synthetic biology. Despite progress from genomic, transcriptomic, and metabolomic approaches, much multiomics data remain underutilized. This review highlights state-of-the-art multiomics strategies for discovering plant biosynthetic pathways, addressing challenges in data acquisition and interpretation with emerging computational tools. We propose an integrated workflow combining molecular networking, reaction pair analysis, and gene expression patterns to enhance data utilization. Additionally, artificial intelligence (AI)-driven approaches promise to revolutionize pathway discovery by streamlining data analysis and validation. Integrating multiomics data, chemical insights, and advanced algorithms can accelerate understanding of plant metabolism and bioengineering valuable natural products efficiently.

Steroidal compounds in Paris polyphylla:structure, biological activities, and biosynthesis

Steroidal compounds are chemical constituents found in the traditional medicinal plant Paris polyphylla, known for their significant pharmacological activities. Due to their complex structures, the biosynthetic pathways of these compounds have garnered considerable attention. In recent years, substantial progress has been made in elucidating the biosynthetic pathways of steroidal compounds from P. polyphylla, with several complete biosynthetic routes being fully characterized. The de novo synthesis of diosgenin has been successfully achieved in both Saccharomyces cerevisiae and Nicotiana benthamiana using various metabolic engineering techniques. Herein, we summarize the latest research progress regarding the structural classification, biological activities, and biosynthesis studies of steroidal compounds from P. polyphylla.

Genome-wide identification of SWEET family genes and functional analysis of NtSWEET12i under drought and saline-alkali stresses in tobacco

BACKGROUND

SWEET (Sugar Will Eventually be Exported Transporter) proteins play vital roles in the transport of sugars, contributing to the regulation of plant development, hormone signaling, and responses to abiotic stress.

RESULTS

In this study, we identified 57 NtSWEETs in tobacco (Nicotiana tabacum L.), and then their physicochemical properties, chromosomal localization, synteny, phylogenetic relationships, genomic structure, promoter cis-elements, protein interaction network, tissue specificity, and expression pattern were systematically analyzed. In addition, NtSWEET12i improves drought and saline-alkali tolerance in tobacco by enhancing soluble sugars transport, ABA signaling, proline biosynthesis, and ROS scavenging. These findings illuminate the role of NtSWEETs, particularly NtSWEET12i, in regulating plant tolerance to drought and saline-alkali stresses.

CONCLUSIONS

This study offers new insights to enhance our understanding of the roles of NtSWEETs and identify potential genes involved in drought and saline-alkali tolerance of plants.

CLINICAL TRIAL NUMBER

Not applicable.

The Deubiquitinating Enzyme AMSH1 Contributes to Plant Immunity Through Regulating the Stability of BDA1

Plants utilize plasma membrane localized receptors like kinases (RLKs) or receptor-like proteins (RLPs) to recognize pathogens and activate pattern-triggered immunity (PTI) responses. A gain-of-function mutation in the Arabidopsis RLP SNC2 (SUPPRESSOR OF NPR1-1, CONSTITUTIVE 2) leads to constitutive activation of defense responses in snc2-1D mutant plants. Transcription factors, SYSTEMIC ACQUIRED RESISTANCE DEFICIENT 1 (SARD1) and CALMODULIN-BINDING PROTEIN 60g (CBP60g), define two parallel pathways downstream of SNC2. The autoimmunity of snc2-1D was partially affected by single mutations in SARD1 or CBP60g but completely suppressed by the sard1 cbp60g double mutant. From a suppressor screen using sard1-1 snc2-1D, we identified a deubiquitinating enzyme ASSOCIATED MOLECULE WITH THE SH3 DOMAIN OF STAM 1 (AMSH1) as a key component in SNC2-mediated plant immunity. A loss-of-function mutation in AMSH1 can suppress the autoimmune responses of sard1-1 snc2-1D. In eukaryotes, selective protein degradation often occurs through the ubiquitination/deubiquitination system. The deubiquitinating enzymes that remove ubiquitin from target proteins play essential roles in controlling the level of target protein ubiquitination and degradation. As loss of AMSH1 results in decreased BDA1 abundance and BDA1 is a transmembrane protein required for SNC2-mediated immunity, AMSH1 likely contributes to immunity regulation through controlling BDA1 stability.

Expression of kiwifruit-derived actinidin in Nicotiana benthamiana leaves

Kiwifruit (Actinidia deliciosa)-derived actinidin, a cysteine protease, is renowned for its meat-tenderizing and milk-clotting activities. Despite its potential in various biotechnological applications, an efficient expression platform for actinidin production has not yet been developed. Instead, actinidin has traditionally been purified directly from the fruits of various plants. This study aimed to produce kiwifruit-derived actinidin in the leaves of Nicotiana benthamiana. The expressed actinidin was directed to the lumen of the endoplasmic reticulum (ER) using the binding immunoglobulin protein (BiP) signal sequence and an ER retention signal. To facilitate cost-effective purification, the family 3 cellulose-binding module (CBM3) was employed as an affinity tag, along with microcrystalline cellulose beads that bind efficiently to CBM3. A significant portion of the expressed actinidin was recovered in the soluble fraction without proteolytic degradation. The purified actinidin exhibited β-casein-degrading activity, with optimal efficiency observed at 55°C and pH 7.0. These results establish a promising plant-based platform for the efficient production and functional application of kiwifruit-derived actinidin in diverse biotechnological processes.

A necrotizing toxin enables Pseudomonas syringae infection across evolutionarily divergent plants

The Pseudomonas syringae species complex harbors a diverse range of pathogenic bacteria that can infect hosts across the plant kingdom. However, much of our current understanding of P. syringae is centered on its infection of flowering plants. We took a comparative approach to understand how P. syringae infects evolutionarily divergent plants. We identified P. syringae isolates causing disease in the liverwort Marchantia polymorpha, the fern Ceratopteris richardii, and the flowering plant Nicotiana benthamiana, which last shared a common ancestor >500 million years ago. Phytotoxin-enriched phylogroup (PG) 2 isolates of P. syringae are virulent in non-flowering plants, relying on type-3 effectors and the lipopeptide phytotoxin syringomycin. Ectopic syringomycin promotes tissue necrosis, activates conserved stress-related genes, and enhances in planta bacterial growth of toxin-deficient PGs in Marchantia. Collectively, our research reveals a key role for syringomycin in promoting Pseudomonas colonization, which works alongside effectors to antagonize an exceptionally wide spectrum of land plants.

Establishing a comprehensive web-based analysis platform for Nicotiana benthamiana genome and transcriptome

Nicotiana benthamiana has long served as a crucial plant material extensively used in plant physiology research, particularly in the field of plant pathology, because of its high susceptibility to plant viruses. Additionally, it serves as a production platform to test vaccines and other valuable substances. Among its approximately 3.1 Gb genome, 57 583 genes have been annotated within a 61 Mb region. We created a comprehensive and easy-to-use platform to use transcriptomes for modern annotation. These tools allow to visualize gene expression profiles, draw molecular evolutionary phylogenetic trees of gene families, perform functional enrichment analyses, and facilitate output downloads. To demonstrate their utility, we analyzed the gene expression profiles of enzymes within the nicotine biosynthesis pathway, a secondary metabolic pathway characteristic of the Nicotiana genus. Using the developed tool, expression profiles of the nicotine biosynthesis pathway genes were generated. The expression patterns of eight gene groups in the pathway were strongly expressed in the roots and weakly expressed in leaves and flowers of N. benthamiana. The results were consistent with the established gene expression profiles in Nicotiana tabacum and provided insights into gene family composition and expression trends. The compilation of this database tool can facilitate genetic analysis of N. benthamiana in the future.

Incorporation of nitrogen in antinutritional Solanum alkaloid biosynthesis

Steroidal glycoalkaloids (SGAs) are specialized metabolites produced by hundreds of Solanum species including food crops, such as tomato, potato and eggplant. Unlike true alkaloids, nitrogen is introduced at a late stage of SGA biosynthesis through an unknown transamination reaction. Here, we reveal the mechanism by which GLYCOALKALOID METABOLISM12 (GAME12) directs the biosynthesis of nitrogen-containing steroidal alkaloid aglycone in Solanum. We report that GAME12, a neofunctionalized γ-aminobutyric acid (GABA) transaminase, undergoes changes in both active site specificity and subcellular localization to switch from its renown and generic activity in core metabolism to function in a specialized metabolic pathway. Moreover, overexpression of GAME12 alone in engineered S. nigrum leaves is sufficient for de novo production of nitrogen-containing SGAs. Our results highlight how hijacking a core metabolism GABA shunt enzyme is crucial in numerous Solanum species for incorporating a nitrogen to a steroidal-specialized metabolite backbone and form defensive alkaloids.

Characterization of the complete mitochondrial genome of the rice bean (Vigna umbellata)

BACKGROUND

Rice bean (Vigna umbellata), an underrated legume crop, demonstrates strong adaptability to poor soil fertility and has significant potential to enhance global food security. It is valuable both as a vegetable and fodder crop due to its high protein content, essential fatty acids, and micronutrients. Despite the sequencing of a high-quality genome of rice bean, its mitochondrial genome (mitogenome) sequence has not yet been reported.

RESULTS

For the first time, the rice bean mitogenome was assembled and annotated using PacBio HiFi sequencing and Geseq software. The mitogenome is a circular molecule with a length of 404,493 bp, containing 32 protein-coding genes, 17 tRNAs, and 3 rRNAs. Codon usage and sequence repeats were also determined. Six gene migration events from the chloroplast to the mitogenome were detected in rice bean. A phylogenetic analysis, including the rice bean mitogenome and 25 other taxa (23 of which are Fabales species), clarified the evolutionary and taxonomic status of rice bean. Additionally, a collinearity analysis of seven Fabales mitogenomes revealed high structural variability. In total, 473 RNA editing sites in protein-coding genes were identified.

CONCLUSIONS

This study presents the first sequencing, assembly, annotation, and analysis of the rice bean mitogenome, providing valuable background information for understanding the evolution of this species. These findings lay the groundwork for future genetic studies and molecular breeding efforts aimed at improving rice bean.

Pairing omics to decode the diversity of plant specialized metabolism

Plants have evolved complex bouquets of specialized natural products that are utilized in medicine, agriculture, and industry. Untargeted natural product discovery has benefitted from growing plant omics data resources. Yet, plant genome complexity limits the identification and curation of biosynthetic pathways via single omics. Pairing multi-omics types within experiments provides multiple layers of evidence for biosynthetic pathway mining. The extraction of paired biological information facilitates connecting genes to transcripts and metabolites, especially when captured across time points, conditions and chemotypes. Experimental design requires specific adaptations to enable effective paired-omics analysis. Ultimately, metadata standards are required to support the integration of paired and unpaired public datasets and to accelerate collaborative efforts for natural product discovery in the plant research community.

The complete genome assembly of Nicotiana benthamiana reveals the genetic and epigenetic landscape of centromeres

Nicotiana benthamiana is a model organism widely adopted in plant biology. Its complete assembly remains unavailable despite several recent improvements. To further improve its usefulness, we generate and phase the complete 2.85 Gb genome assembly of allotetraploid N. benthamiana. We find that although Solanaceae centromeres are widely dominated by Ty3/Gypsy retrotransposons, satellite-based centromeres are surprisingly common in N. benthamiana, with 11 of 19 centromeres featured by megabase-scale satellite arrays. Interestingly, the satellite-enriched and satellite-free centromeres are extensively invaded by distinct Gypsy retrotransposons which CENH3 protein more preferentially occupies, suggestive of their crucial roles in centromere function. We demonstrate that ribosomal DNA is a major origin of centromeric satellites, and mitochondrial DNA could be employed as a core component of the centromere. Subgenome analysis indicates that the emergence of satellite arrays probably drives new centromere formation. Altogether, we propose that N. benthamiana centromeres evolved via neocentromere formation, satellite expansion, retrotransposon enrichment and mtDNA integration.

Differential prolyl hydroxylation by six Physcomitrella prolyl-4 hydroxylases

Hydroxylation of prolines to 4-trans-hydroxyproline (Hyp) is mediated by prolyl-4 hydroxylases (P4Hs). In plants, Hyps occur in Hydroxyproline-rich glycoproteins (HRGPs), and are frequently O-glycosylated. While both modifications are important, e.g. for cell wall stability, they are undesired in plant-made pharmaceuticals. Sequence motifs for prolyl-hydroxylation were proposed but did not include data from mosses, such as Physcomitrella. We identified six moss P4Hs by phylogenetic reconstruction. Our analysis of 73 Hyps in 24 secretory proteins from multiple mass spectrometry datasets revealed that prolines near other prolines, alanine, serine, threonine and valine were preferentially hydroxylated. About 95 % of Hyps were predictable with combined established methods. In our data, AOV was the most frequent pattern. A combination of 443 AlphaFold models and MS data with 3000 prolines found Hyps mainly on protein surfaces in disordered regions. Moss-produced human erythropoietin (EPO) exhibited O-glycosylation with arabinose chains on two Hyps. This modification was significantly reduced in a p4h1 knock-out (KO) Physcomitrella mutant. Quantitative proteomics with different p4h mutants revealed specific changes in protein amounts, and a modified prolyl-hydroxylation pattern, suggesting a differential function of the Physcomitrella P4Hs. Quantitative RT-PCR revealed a differential effect of single p4h KOs on the expression of the other five p4h genes, suggesting a partial compensation of the mutation. AlphaFold-Multimer models for Physcomitrella P4H1 and its target EPO peptide superposed with the crystal structure of Chlamydomonas P4H1 suggested significant amino acids in the active centre of the enzyme and revealed differences between P4H1 and the other Physcomitrella P4Hs.

Viral proteins resolve the virus-vector conundrum during hemipteran-mediated transmission by subverting salicylic acid signaling pathway

Hemipteran insects transmit viruses when infesting plants, during which vectors activate salicylic acid (SA)-regulated antiviral defenses. How vector-borne plant viruses circumvent these antiviral defenses is largely unexplored. During co-infections of begomoviruses and betasatellites in plants, betasatellite-encoded βC1 proteins interfere with SA signaling and reduce the activation of antiviral resistance. βC1 inhibits SA-induced degradation of NbNPR3 (Nicotiana benthamiana nonexpressor of pathogenesis-related genes 3), a negative regulator of SA signaling. βC1 does not bind directly to NbNPR3, but regulates NbNPR3 degradation via heat shock protein 90s (NbHSP90s). NbHSP90s bind to both NbNPR3 and βC1 and suppress SA signaling. This viral success strategy appears to be conserved as it is also documented for viral proteins encoded by two aphid-borne viruses. Our findings reveal an exquisite mechanism that facilitates the persistence of vector-borne plant viruses and provide important insights into the intricacies of the virus life cycle.

Extracellular plant subtilases dampen cold-shock peptide elicitor levels

Recognizing pathogen-associated molecular patterns on the cell surface is crucial for plant immunity. The proteinaceous nature of many of these patterns suggests that secreted proteases play important roles in their formation and stability. Here we demonstrate that the apoplastic subtilase SBT5.2a inactivates the immunogenicity of cold-shock proteins (CSPs) of the bacterial plant pathogen Pseudomonas syringae by cleaving within the immunogenic csp22 epitope. Consequently, mutant plants lacking SBT5.2a activity retain higher levels of csp22, leading to enhanced immune responses and reduced pathogen growth. SBT5.2 sensitivity is influenced by sequence variation surrounding the cleavage site and probably extends to CSPs from other bacterial species. These findings suggest that variations in csp22 stability among bacterial pathogens are a crucial factor in plant-bacteria interactions and that pathogens exploit plant proteases to avoid pattern recognition.

Engineering Nicotiana benthamiana as a platform for natural product biosynthesis

Nicotiana benthamiana is a model plant, widely used for research. The susceptibility of young plants to Agrobacterium tumefaciens has been utilised for transient gene expression, enabling the production of recombinant proteins at laboratory and commercial scales. More recently, this technique has been used for the rapid prototyping of synthetic genetic circuits and for the elucidation and reconstruction of metabolic pathways. In the last few years, many complex metabolic pathways have been successfully reconstructed in this species. In addition, the availability of improved genomic resources and efficient gene editing tools have enabled the application of sophisticated metabolic engineering approaches to increase the purity and yield of target compounds. In this review, we discuss recent advances in the use of N. benthamiana for understanding and engineering plant metabolism, as well as efforts to improve the utility of this species as a production chassis for natural products.

Evidence for within-species transition between drought response strategies in Nicotiana benthamiana

Nicotiana benthamiana is predominantly distributed in arid habitats across northern Australia. However, none of six geographically isolated accessions shows obvious xerophytic morphological features. To investigate how these tender-looking plants withstand drought, we examined their responses to water deprivation, assessed phenotypic, physiological, and cellular responses, and analysed cuticular wax composition and wax biosynthesis gene expression profiles. Results showed that the Central Australia (CA) accession, globally known as a research tool, has evolved a drought escape strategy with early vigour, short life cycle, and weak, water loss-limiting responses. By contrast, a northern Queensland (NQ) accession responded to drought by slowing growth, inhibiting flowering, increasing leaf cuticle thickness, and altering cuticular wax composition. Under water stress, NQ increased the heat stability and water impermeability of its cuticle by extending the carbon backbone of cuticular long-chain alkanes from c. 25 to 33. This correlated with rapid upregulation of at least five wax biosynthesis genes. In CA, the alkane chain lengths (c. 25) and gene expression profiles remained largely unaltered. This study highlights complex genetic and environmental control over cuticle composition and provides evidence for divergence into at least two fundamentally different drought response strategies within the N. benthamiana species in < 1 million years.

Development of a tightly regulated copper-inducible transient gene expression system in Nicotiana benthamiana incorporating a suicide exon and Cre recombinase

Chemical-inducible gene expression systems are commonly used to regulate gene expression for functional genomics in various plant species. However, a convenient system that can tightly regulate transgene expression in Nicotiana benthamiana is still lacking. In this study, we developed a tightly regulated copper-inducible system that can control transgene expression and conduct cell death assays in N. benthamiana. We tested several chemical-inducible systems using Agrobacterium-mediated transient expression and found that the copper-inducible system exhibited the least concerns regarding leakiness in N. benthamiana. Although the copper-inducible system can control the expression of some tested reporters, it is not sufficiently tight to regulate certain tested hypersensitive cell death responses. Using the MoClo-based synthetic biology approach, we incorporated the suicide exon HyP5SM/OsL5 and Cre/LoxP as additional regulatory elements to enhance the tightness of the regulation. This new design allowed us to tightly control the hypersensitive cell death induced by several tested leucine-rich repeat-containing proteins and their matching avirulence factors, and it can be easily applied to regulate the expression of other transgenes in transient expression assays. Our findings offer new approaches for both fundamental and translational studies in plant functional genomics.

Retrospect and prospect of Nicotiana tabacum genome sequencing

Investigating plant genomes offers crucial foundational resources for exploring various aspects of plant biology and applications, such as functional genomics and breeding practices. With the development in sequencing and assembly technology, several Nicotiana tabacum genomes have been published. In this paper, we reviewed the progress on N. tabacum genome assembly and quality, from the initial draft genomes to the recent high-quality chromosome-level assemblies. The application of long-read sequencing, optical mapping, and Hi-C technologies has significantly improved the contiguity and completeness of N. tabacum genome assemblies, with the latest assemblies having a contig N50 size over 50 Mb. Despite these advancements, further improvements are still required and possible, particularly on the development of pan-genome and telomere-to-telomere (T2T) genomes. These new genomes will capture the genomic diversity and variations among different N. tabacum cultivars and species, and provide a comprehensive view of the N. tabacum genome structure and gene content, so to deepen our understanding of the N. tabacum genome and facilitate precise breeding and functional genomics.

Basal ABA signaling balances transpiration and photosynthesis

The balance between the CO2 entry for photosynthesis and transpiration water loss is crucial for plant growth, and ABA signaling can affect this equilibrium. To test how ABA balances plant growth and environmental adaptation, we performed molecular genetics studies in the biotech crop Nicotiana benthamiana under well-watered or drought conditions. Studies on ABA signaling in crops are complicated by the multigenic nature of the PYR/PYL/RCAR ABA receptor family and its functional redundancy, which is particularly challenging in polyploid plants. We have generated a pentuple pyl mutant in the allotetraploid Nicotiana benthamiana through CRISPR/Cas9 gene editing. The pentuple mutant is impaired in 2 NbPYL1-like and 3 NbPYL8-like receptors, affecting the regulation of transpiration and several ABA-dependent transcriptional processes. RNA-seq and metabolite analysis revealed that the synthesis of galactinol, an essential precursor for the osmoprotective raffinose family of oligosaccharides, is ABA-dependent and impaired in the mutant under osmotic stress. In contrast, our results show that, under well-watered conditions, partial inactivation of ABA signaling leads to higher CO2 entry and photosynthesis in the mutant than in WT. Photosynthesis analyses revealed an increased CO2 diffusion capacity mediated by higher stomatal and mesophyll conductances, and higher substomatal CO2 concentration in the pentuple mutant. RNA-seq analyses revealed that genes associated with cell wall loosening (e.g., expansins) and porosity were strongly downregulated by ABA in WT. In summary, a partial relief of the ABA control on transpiration mediated by ABA receptors positively affects photosynthesis when water is not limited, at the expense of reduced water use efficiency.

The proteome of Nicotiana benthamiana is shaped by extensive protein processing

Processing by proteases irreversibly regulates the fate of plant proteins and hampers the production of recombinant proteins in plants, yet only few processing events have been described in agroinfiltrated Nicotiana benthamiana, which has emerged as the main transient protein expression platform in plant science and molecular pharming. Here, we used in-gel digests and mass spectrometry to monitor the migration and topography of 5040 plant proteins within a protein gel. By plotting the peptides over the gel slices, we generated peptographs that reveal where which part of each protein was detected within the protein gel. These data uncovered that 60% of the detected proteins have proteoforms that migrate at lower than predicted molecular weights, implicating extensive proteolytic processing. This analysis confirms the proteolytic removal and degradation of autoinhibitory prodomains of most but not all proteases, and revealed differential processing within pectinemethylesterase and lipase families. This analysis also uncovered intricate processing of glycosidases and uncovered that ectodomain shedding might be common for a diverse range of receptor-like kinases. Transient expression of double-tagged candidate proteins confirmed processing events in vivo. This large proteomic dataset implicates an elaborate proteolytic machinery shaping the proteome of N. benthamiana.

Beyond co-expression: pathway discovery for plant pharmaceuticals

Plant natural products have been an important source of medicinal molecules since ancient times. To gain access to the whole diversity of these molecules for pharmaceutical applications, it is important to understand their biosynthetic origins. Whilst co-expression is a reliable tool for identifying gene candidates, a variety of complementary methods can aid in screening or refining candidate selection. Here, we review recently employed plant biosynthetic pathway discovery approaches, and highlight future directions in the field.

Genetically encoded Boolean logic operators to sense and integrate phenylpropanoid metabolite levels in plants

Synthetic biology has the potential to revolutionize biotechnology, public health, and agriculture. Recent studies have shown the enormous potential of plants as chassis for synthetic biology applications. However, tools to precisely manipulate metabolic pathways for bioproduction in plants are still needed. We used bacterial allosteric transcription factors (aTFs) that control gene expression in a ligand-specific manner and tested their ability to repress semi-synthetic promoters in plants. We also tested the modulation of their repression activity in response to specific plant metabolites, especially phenylpropanoid-related molecules. Using these aTFs, we also designed synthetic genetic circuits capable of computing Boolean logic operations. Three aTFs, CouR, FapR, and TtgR, achieved c. 95% repression of their respective target promoters. For TtgR, a sixfold de-repression could be triggered by inducing its ligand accumulation, showing its use as biosensor. Moreover, we designed synthetic genetic circuits that use AND, NAND, IMPLY, and NIMPLY Boolean logic operations and integrate metabolite levels as input to the circuit. We showed that biosensors can be implemented in plants to detect phenylpropanoid-related metabolites and activate a genetic circuit that follows a predefined logic, demonstrating their potential as tools for exerting control over plant metabolic pathways and facilitating the bioproduction of natural products.

A vicinal oxygen chelate protein facilitates viral infection by triggering the unfolded protein response in Nicotiana benthamiana

Vicinal oxygen chelate (VOC) proteins are members of an enzyme superfamily with dioxygenase or non-dioxygenase activities. However, the biological functions of VOC proteins in plants are poorly understood. Here, we show that a VOC in Nicotiana benthamiana (NbVOC1) facilitates viral infection. NbVOC1 was significantly induced by infection by beet necrotic yellow vein virus (BNYVV). Transient overexpression of NbVOC1 or its homolog from Beta vulgaris (BvVOC1) enhanced BNYVV infection in N. benthamiana, which required the nuclear localization of VOC1. Consistent with this result, overexpressing NbVOC1 facilitated BNYVV infection, whereas, knockdown and knockout of NbVOC1 inhibited BNYVV infection in transgenic N. benthamiana plants. NbVOC1 interacts with the basic leucine zipper transcription factors bZIP17/28, which enhances their self-interaction and DNA binding to the promoters of unfolded protein response (UPR)-related genes. We propose that bZIP17/28 directly binds to the NbVOC1 promoter and induces its transcription, forming a positive feedback loop to induce the UPR and facilitating BNYVV infection. Collectively, our results demonstrate that NbVOC1 positively regulates the UPR that enhances viral infection in plants.

Using Synthetic Biology to Understand the Function of Plant Specialized Metabolites

Plant specialized metabolites (PSMs) are variably distributed across taxa, tissues, and ecological contexts; this variability has inspired many theories about PSM function, which, to date, remain poorly tested because predictions have outpaced the available data. Advances in mass spectrometry-based metabolomics have enabled unbiased PSM profiling, and molecular biology techniques have produced PSM-free plants; the combination of these methods has accelerated our understanding of the complex ecological roles that PSMs play in plants. Synthetic biology techniques and workflows are producing high-value, structurally complex PSMs in quantities and purities sufficient for both medicinal and functional studies. These workflows enable the reengineering of PSM transport, externalization, structural diversity, and production in novel taxa, facilitating rigorous tests of long-standing theoretical predictions about why plants produce so many different PSMs in particular tissues and ecological contexts. Plants use their chemical prowess to solve ecological challenges, and synthetic biology workflows are accelerating our understanding of these evolved functions.

Deficiency of multiple RNA silencing-associated genes may contribute to the increased susceptibility of Nicotiana benthamiana to viruses

KEY MESSAGE

Recently published high-quality reference genome assemblies indicate that, in addition to RDR1-deficiency, the loss of several key RNA silencing-associated genes may contribute to the hypersusceptibility of Nicotiana benthamiana to viruses.

Exploring the source of TYLCV resistance in Nicotiana benthamiana

Tomato Yellow Leaf Curl Virus (TYLCV) is one of the most devastating pathogens of tomato, worldwide. It is vectored by the globally prevalent whitefly, Bemisia tabaci, and is asymptomatic in a wide range of plant species that act as a virus reservoir. The most successful crop protection for tomato in the field has been from resistance genes, of which five loci have been introgressed fromwild relatives. Of these, the Ty-1/Ty-3 locus, which encodes an RNA-dependent RNA polymerase 3 (RDR3), has been the most effective. Nevertheless, several TYLCV strains that break this resistance are beginning to emerge, increasing the need for new sources of resistance. Here we use segregation analysis and CRISPR-mediated gene dysfunctionalisation to dissect the differential response of two isolates of Nicotiana benthamiana to TYLCV infection. Our study indicates the presence of a novel non-RDR3, but yet to be identified, TYLCV resistance gene in a wild accession of N. benthamiana. This gene has the potential to be incorporated into tomatoes.

Heterologous expression of influenza haemagglutinin leads to early and transient activation of the unfolded protein response in Nicotiana benthamiana

The unfolded protein response (UPR) allows cells to cope with endoplasmic reticulum (ER) stress induced by accumulation of misfolded proteins in the ER. Due to its sensitivity to Agrobacterium tumefaciens, the model plant Nicotiana benthamiana is widely employed for transient expression of recombinant proteins of biopharmaceutical interest, including antibodies and virus surface proteins used for vaccine production. As such, study of the plant UPR is of practical significance, since enforced expression of complex secreted proteins often results in ER stress. After 6 days of expression, we recently reported that influenza haemagglutinin H5 induces accumulation of UPR proteins. Since up-regulation of corresponding UPR genes was not detected at this time, accumulation of UPR proteins was hypothesized to be independent of transcriptional induction, or associated with early but transient UPR gene up-regulation. Using time course sampling, we here show that H5 expression does result in early and transient activation of the UPR, as inferred from unconventional splicing of NbbZIP60 transcripts and induction of UPR genes with varied functions. Transient nature of H5-induced UPR suggests that this response was sufficient to cope with ER stress provoked by expression of the secreted protein, as opposed to an antibody that triggered stronger and more sustained UPR activation. As up-regulation of defence genes responding to H5 expression was detected after the peak of UPR activation and correlated with high increase in H5 protein accumulation, we hypothesize that these immune responses, rather than the UPR, were responsible for onset of the necrotic symptoms on H5-expressing leaves.

A transient expression tool box for anthocyanin biosynthesis in Nicotiana benthamiana

Transient expression in Nicotiana benthamiana offers a robust platform for the rapid production of complex secondary metabolites. It has proven highly effective in helping identify genes associated with pathways responsible for synthesizing various valuable natural compounds. While this approach has seen considerable success, it has yet to be applied to uncovering genes involved in anthocyanin biosynthetic pathways. This is because only a single anthocyanin, delphinidin 3-O-rutinoside, can be produced in N. benthamiana by activation of anthocyanin biosynthesis using transcription factors. The production of other anthocyanins would necessitate the suppression of certain endogenous flavonoid biosynthesis genes while transiently expressing others. In this work, we present a series of tools for the reconstitution of anthocyanin biosynthetic pathways in N. benthamiana leaves. These tools include constructs for the expression or silencing of anthocyanin biosynthetic genes and a mutant N. benthamiana line generated using CRISPR. By infiltration of defined sets of constructs, the basic anthocyanins pelargonidin 3-O-glucoside, cyanidin 3-O-glucoside and delphinidin 3-O-glucoside could be obtained in high amounts in a few days. Additionally, co-infiltration of supplementary pathway genes enabled the synthesis of more complex anthocyanins. These tools should be useful to identify genes involved in the biosynthesis of complex anthocyanins. They also make it possible to produce novel anthocyanins not found in nature. As an example, we reconstituted the pathway for biosynthesis of Arabidopsis anthocyanin A5, a cyanidin derivative and achieved the biosynthesis of the pelargonidin and delphinidin variants of A5, pelargonidin A5 and delphinidin A5.

High-quality chromosome-level genome assembly of Nicotiana benthamiana

Nicotiana benthamiana is a fundamental model organism in plant research. Recent advancements in genomic sequencing have revealed significant intraspecific genetic variations. This study addresses the pressing need for a precise genome sequence specific to its geographic origin by presenting a comprehensive genome assembly of the N. benthamiana LAB strain from the Republic of Korea (NbKLAB). We compare this assembly with the widely used NbLAB360 strain, shedding light on essential genomic differences between them. The outcome is a high-quality, chromosome-level genome assembly comprising 19 chromosomes, spanning 2,762 Mb, with an N50 of 142.6 Mb. Comparative analyses revealed notable variations, including 46,215 protein-coding genes, with an impressive 99.5% BUSCO completeness score. Furthermore, the NbKLAB assembly substantially improved the QV from 33% for NbLAB360 to 49%. This refined chromosomal genome assembly for N. benthamiana, in conjunction with comparative insights, provides a valuable resource for genomics research and molecular biology. This accomplishment forms a strong foundation for in-depth exploration into the intricacies of plant genetics and genomics, improved precision, and a comparative framework.

High-quality assembled and annotated genomes of Nicotiana tabacum and Nicotiana benthamiana reveal chromosome evolution and changes in defense arsenals

Nicotiana tabacum and Nicotiana benthamiana are widely used models in plant biology research. However, genomic studies of these species have lagged. Here we report the chromosome-level reference genome assemblies for N. benthamiana and N. tabacum with an estimated 99.5% and 99.8% completeness, respectively. Sensitive transcription start and termination site sequencing methods were developed and used for accurate gene annotation in N. tabacum. Comparative analyses revealed evidence for the parental origins and chromosome structural changes, leading to hybrid genome formation of each species. Interestingly, the antiviral silencing genes RDR1, RDR6, DCL2, DCL3, and AGO2 were lost from one or both subgenomes in N. benthamiana, while both homeologs were kept in N. tabacum. Furthermore, the N. benthamiana genome encodes fewer immune receptors and signaling components than that of N. tabacum. These findings uncover possible reasons underlying the hypersusceptible nature of N. benthamiana. We developed the user-friendly Nicomics (http://lifenglab.hzau.edu.cn/Nicomics/) web server to facilitate better use of Nicotiana genomic resources as well as gene structure and expression analyses.

Establishment and application of Agrobacterium-delivered CRISPR/Cas9 system for wild tobacco (Nicotiana alata) genome editing

Clustered regularly interspaced short palindromic repeats (CRISPR)-associated protein 9 (CRISPR-Cas9) system has been widely applied in cultivated crops, but limited in their wild relatives. Nicotiana alata is a typical wild species of genus Nicotiana that is globally distributed as a horticultural plant and well-studied as a self-incompatibility model. It also has valuable genes for disease resistance and ornamental traits. However, it lacks an efficient genetic transformation and genome editing system, which hampers its gene function and breeding research. In this study, we developed an optimized hypocotyl-mediated transformation method for CRISPR-Cas9 delivery. The genetic transformation efficiency was significantly improved from approximately 1% to over 80%. We also applied the CRISPR-Cas9 system to target the phytoene desaturase (NalaPDS) gene in N. alata and obtained edited plants with PDS mutations with over 50% editing efficiency. To generate self-compatible N. alata lines, a polycistronic tRNA-gRNA (PTG) strategy was used to target exonic regions of allelic S-RNase genes and generate targeted knockouts simultaneously. We demonstrated that our system is feasible, stable, and high-efficiency for N. alata genome editing. Our study provides a powerful tool for basic research and genetic improvement of N. alata and an example for other wild tobacco species.

Generation and physiological characterization of genome-edited Nicotiana benthamiana plants containing zeaxanthin as the only leaf xanthophyll

MAIN CONCLUSION

Simultaneous genome editing of the two homeologous LCYe and ZEP genes of Nicotiana benthamiana results in plants in which all xanthophylls are replaced by zeaxanthin. Plant carotenoids act both as photoreceptors and photoprotectants in photosynthesis and as precursors of apocarotenoids, which include signaling molecules such as abscisic acid (ABA). As dietary components, the xanthophylls lutein and zeaxanthin have photoprotective functions in the human macula. We developed transient and stable combinatorial genome editing methods, followed by direct LC-MS screening for zeaxanthin accumulation, for the simultaneous genome editing of the two homeologous Lycopene Epsilon Cyclase (LCYe) and the two Zeaxanthin Epoxidase (ZEP) genes present in the allopolyploid Nicotiana benthamiana genome. Editing of the four genes resulted in plants in which all leaf xanthophylls were substituted by zeaxanthin, but with different ABA levels and growth habits, depending on the severity of the ZEP1 mutation. In high-zeaxanthin lines, the abundance of the major photosystem II antenna LHCII was reduced with respect to wild-type plants and the LHCII trimeric state became unstable upon thylakoid solubilization. Consistent with the depletion in LHCII, edited plants underwent a compensatory increase in PSII/PSI ratios and a loss of the large-size PSII supercomplexes, while the level of PSI-LHCI supercomplex was unaffected. Reduced activity of the photoprotective mechanism NPQ was shown in high-zeaxanthin plants, while PSII photoinhibition was similar for all genotypes upon exposure to excess light, consistent with the antioxidant and photoprotective role of zeaxanthin in vivo.