A draft genome sequence of Nicotiana benthamiana to enhance molecular plant-microbe biology research

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.

Role of LEAFLESS, an AP2/ERF family transcription factor, in the regulation of trichome specialized metabolism

Acylsugars, specialized metabolites produced by trichomes of many solanaceous species, provide protection against biotic and abiotic stresses. Many acylsugar metabolic enzymes have been identified; however, regulatory factors remain unknown. Our multidisciplinary approaches identified LEAFLESS (APETALA 2/ ETHYLENE RESPONSE FACTOR (AP2/ERF) family member) as a positive regulator of acylsugar biosynthesis. Virus-induced gene silencing (VIGS) of LEAFLESS in Solanum pennellii (SpLFS/Sopen05g008450) revealed its distinct roles in two related but separate processes: acylsugar biosynthesis and trichome development. Most acylsugar (and several flavonoid) metabolic genes were downregulated in SpLFS-silenced plants and showed strong co-expression with SpLFS. Phylogenetic and additional data analyses indicated trichome-enriched expression of SpLFS orthologs in other acylsugar-producing solanaceous species, and VIGS of SpLFS orthologs in Nicotiana benthamiana reduced acylsugar production. Transcriptional reporter showed expression of SpLFS in type I/IV trichome tip cells, the site of acylsugar biosynthesis. Electrophoretic mobility shift assays indicated that SpLFS directly binds to promoters of several acylsugar (and flavonoid) metabolic genes. Additionally, data mining suggested remarkable spatiotemporal functional diversity: from coordinating leaf initiation at incipient primordia (previously reported for the S. lycopersicum ortholog SlLFS/Solyc05g013540) to regulating trichome specialized metabolism (acylsugar and flavonoid). Our work highlights a critical role of LEAFLESS in trichome specialized metabolism, paving the way to disentangle the acylsugar regulatory network.

Evolution of NAC transcription factors from early land plants to domesticated crops

NAC [NO APICAL MERISTEM (NAM), ARABIDOPSIS TRANSCRIPTION ACTIVATOR FACTOR 1/2 (ATAF1/2), and CUP-SHAPED COTYLEDON (CUC2)] transcription factors are key regulators of plant growth, development, and stress responses but were also crucial players during land plant adaptation and crop domestication. Using representative members of green algae, bryophytes, lycophytes, gymnosperms, and angiosperms, we expanded the evolutionary history of NAC transcription factors to unveil the relationships among members of this gene family. We found a massive increase in the number of NAC transcription factors from green algae to lycophytes and an even larger increase in flowering plants. Many of the NAC clades arose later during evolution since we found eudicot- and monocot-specific clades. Cis-elements analysis in NAC promoters showed the presence of abiotic and biotic stress as well as hormonal response elements, which indicate the ancestral function of NAC transcription factor genes in response to environmental stimuli and in plant development. At the transcriptional level, the expression of NAC transcription factors was low or absent in male reproduction, particularly mature pollen, across the plant kingdom. We also identified NAC genes with conserved expression patterns in response to heat stress in Marchantia polymorpha and Oryza sativa. Our study provides further evidence that transcriptional mechanisms associated with stress responses and development emerged early during plant land adaptation and are still conserved in flowering plants and domesticated crops.

A new level of RNA-based plant protection: dsRNAs designed from functionally characterized siRNAs highly effective against Cucumber mosaic virus

RNA-mediated crop protection increasingly becomes a viable alternative to agrochemicals that threaten biodiversity and human health. Pathogen-derived double-stranded RNAs (dsRNAs) are processed into small interfering RNAs (siRNAs), which can then induce silencing of target RNAs, e.g. viral genomes. However, with currently used dsRNAs, which largely consist of undefined regions of the target RNAs, silencing is often ineffective: processing in the plant generates siRNA pools that contain only a few functionally effective siRNAs (esiRNAs). Using an in vitro screen that reliably identifies esiRNAs from siRNA pools, we identified esiRNAs against Cucumber mosaic virus (CMV), a devastating plant pathogen. Topical application of esiRNAs to plants resulted in highly effective protection against massive CMV infection. However, optimal protection was achieved with newly designed multivalent 'effective dsRNAs' (edsRNAs), which contain the sequences of several esiRNAs and are preferentially processed into these esiRNAs. The esiRNA components can attack one or more target RNAs at different sites, be active in different silencing complexes, and provide cross-protection against different viral variants-important properties for combating rapidly mutating pathogens such as CMV. esiRNAs and edsRNAs have thus been established as a new class of 'RNA actives' that significantly increase the efficacy and specificity of RNA-mediated plant protection.

Virus-Induced Gene Silencing to Unravel the Function of Nicotiana benthamiana Genes Linked to Corolla Abscission

Virus-induced gene silencing (VIGS) is a powerful and cost-effective biotechnological tool for studying gene function in plants. It allows for rapid downregulation of gene expression without the need for tissue culture or genetic transformation. While VIGS is transient and may pose challenges for certain developmental genes, it has been successfully used to study organ abscission. Using a VIGS protocol based on Citrus leaf blotch virus (CLBV) viral vectors, we were able to characterize the regulatory role of INFLORESCENCE DEFICIENT IN ABSCISSION (IDA)-like genes and their potential HAESA (HAE)-like receptors in the process of corolla abscission in Nicotiana benthamiana flowers.

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.

Dominant-Negative KAI2d Paralogs Putatively Attenuate Strigolactone Responses in Root Parasitic Plants

Many root parasitic plants in the Orobanchaceae use host-derived strigolactones (SLs) as germination cues. This adaptation facilitates attachment to a host and is particularly important for the success of obligate parasitic weeds that cause substantial crop losses globally. Parasite seeds sense SLs through 'divergent' KARRIKIN INSENSITIVE2 (KAI2d)/HYPOSENSITIVE TO LIGHT α/β-hydrolases that have undergone substantial duplication and diversification in Orobanchaceae genomes. After germination, chemotropic growth of parasite roots toward a SL source also occurs in some species. We investigated which of the seven KAI2d genes found in a facultative hemiparasite, Phtheirospermum japonicum, may enable chemotropic responses to SLs. To do so, we developed a triple mutant Nbd14a,b kai2i line of Nicotiana benthamiana in which SL-induced degradation of SUPPRESSOR OF MORE AXILLARY GROWTH2 (MAX2) 1 (SMAX1), an immediate downstream target of KAI2 signaling, is disrupted. In combination with a transiently expressed, ratiometric reporter of SMAX1 protein abundance, this mutant forms a system for the functional analysis of parasite KAI2d proteins in a plant cellular context. Using this system, we unexpectedly found three PjKAI2d proteins that do not trigger SMAX1 degradation in the presence of SLs. Instead, these PjKAI2d proteins inhibit the perception of low SL concentrations by SL-responsive PjKAI2d in a dominant-negative manner that depends upon an active catalytic triad. Similar dominant-negative KAI2d paralogs were identified in an obligate hemiparasitic weed, Striga hermonthica. These proteins suggest a mechanism for attenuating SL signaling in parasites, which might be used to enhance the perception of shallow SL gradients during root growth toward a host or to restrict germination responses to specific SLs.

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.

NbPTR1 confers resistance against Pseudomonas syringae pv. actinidiae in kiwifruit

Pseudomonas syringae pv. actinidiae biovar 3 (Psa3) causes a devastating canker disease in yellow-fleshed kiwifruit (Actinidia chinensis). The effector HopZ5, which is present in all isolates of Psa3 causing global outbreaks of pandemic kiwifruit canker disease, triggers immunity in Nicotiana benthamiana and is not recognised in susceptible A. chinensis cultivars. In a search for N. benthamiana nonhost resistance genes against HopZ5, we found that the nucleotide-binding leucine-rich repeat receptor NbPTR1 recognised HopZ5. RPM1-interacting protein 4 orthologues from N. benthamiana and A. chinensis formed a complex with NbPTR1 and HopZ5 activity was able to disrupt this interaction. No functional orthologues of NbPTR1 were found in A. chinensis. NbPTR1 transformed into Psa3-susceptible A. chinensis var. chinensis 'Hort16A' plants introduced HopZ5-specific resistance against Psa3. Altogether, this study suggested that expressing NbPTR1 in Psa3-susceptible kiwifruit is a viable approach to acquiring resistance to Psa3 and it provides valuable information for engineering resistance in otherwise susceptible kiwifruit genotypes.

Chlorophyll to zeaxanthin energy transfer in nonphotochemical quenching: An exciton annihilation-free transient absorption study

Zeaxanthin (Zea) is a key component in the energy-dependent, rapidly reversible, nonphotochemical quenching process (qE) that regulates photosynthetic light harvesting. Previous transient absorption (TA) studies suggested that Zea can participate in direct quenching via chlorophyll (Chl) to Zea energy transfer. However, the contamination of intrinsic exciton-exciton annihilation (EEA) makes the assignment of TA signal ambiguous. In this study, we present EEA-free TA data using Nicotiana benthamiana thylakoid membranes, including the wild type and three NPQ mutants (npq1, npq4, and lut2) generated by CRISPR/Cas9 mutagenesis. The results show a strong correlation between excitation energy transfer from excited Chl Qy to Zea S1 and the xanthophyll cycle during qE activation. Notably, a Lut S1 signal is absent in the npq1 thylakoids which lack zeaxanthin. Additionally, the fifth-order response analysis shows a reduction in the exciton diffusion length (LD) from 62 ± 6 nm to 43 ± 3 nm under high light illumination, consistent with the reduced range of exciton motion being a key aspect of plants' response to excess light.

Two plant membrane-shaping reticulon-like proteins play contrasting complex roles in turnip mosaic virus infection

Positive-sense RNA viruses remodel cellular cytoplasmic membranes as the membranous sources for the formation of viral replication organelles (VROs) for viral genome replication. In plants, they traffic through plasmodesmata (PD), plasma membrane-lined pores enabling cytoplasmic connections between cells for intercellular movement and systemic infection. In this study, we employed turnip mosaic virus (TuMV), a plant RNA virus to investigate the involvement of RTNLB3 and RTNLB6, two ER (endoplasmic reticulum) membrane-bending, PD-located reticulon-like (RTNL) non-metazoan group B proteins (RTNLBs) in viral infection. We show that RTNLB3 interacts with TuMV 6K2 integral membrane protein and RTNLB6 binds to TuMV coat protein (CP). Knockdown of RTNLB3 promoted viral infection, whereas downregulation of RTNLB6 restricted viral infection, suggesting that these two RTNLs play contrasting roles in TuMV infection. We further demonstrate that RTNLB3 targets the α-helix motif 42LRKSM46 of 6K2 to interrupt 6K2 self-interactions and compromise 6K2-induced VRO formation. Moreover, overexpression of AtRTNLB3 apparently promoted the selective degradation of the ER and ER-associated protein calnexin, but not 6K2. Intriguingly, mutation of the α-helix motif of 6K2 that is required for induction of VROs severely affected 6K2 stability and abolished TuMV infection. Thus, RTNLB3 attenuates TuMV replication, probably through the suppression of 6K2 function. We also show that RTNLB6 promotes viral intercellular movement but does not affect viral replication. Therefore, the proviral role of RTNLB6 is probably by enhancing viral cell-to-cell trafficking. Taken together, our data demonstrate that RTNL family proteins may play diverse complex, even opposite, roles in viral infection in plants.

Serine proteinase inhibitors from Nicotiana benthamiana, a nonpreferred host plant, inhibit the growth of Myzus persicae (green peach aphid)

BACKGROUND

The green peach aphid (Myzus persicae) is a severe agricultural crop pest that has developed resistance to most current control methods, requiring the urgent development of novel strategies. Plant proteinase inhibitors (PINs) are small proteins that protect plants against pathogens and/or herbivores, likely by preventing efficient protein digestion.

RESULTS

We identified 67 protease genes in the transcriptomes of three M. persicae lineages (USDA-Red, G002 and G006). Comparison of gene expression levels in aphid guts and whole aphids showed that several proteases, including a highly expressed serine protease, are significantly overexpressed in the guts. Furthermore, we identified three genes encoding serine protease inhibitors (SerPIN-II1, 2 and 3) in Nicotiana benthamiana, which is a nonpreferred host for M. persicae. Using virus-induced gene silencing (VIGS) with a tobacco rattle virus (TRV) vector and overexpression with a turnip mosaic virus (TuMV) vector, we demonstrated that N. benthamiana SerPIN-II1 and SerPIN-II2 cause reduced survival and growth, but do not affect aphid protein content. Likewise, SerPIN-II3 overexpression reduced survival and growth, and serpin-II3 knockout mutations, which we generated using CRISPR/Cas9, increased survival and growth. Protein content was significantly increased in aphids fed on SerPIN-II3 overexpressing plants, yet it was decreased in aphids fed on serpin-II3 mutants.

CONCLUSION

Our results show that three PIN-IIs from N. benthamiana, a nonpreferred host plant, effectively inhibit M. persicae survival and growth, thereby representing a new resource for the development of aphid-resistant crop plants. © 2024 Society of Chemical Industry.

Genome-wide identification, characterization and expression analysis of the bZIP transcription factors in garlic (Allium sativum L.)

Introduction

The bZIP genes (bZIPs) are essential in numerous biological processes, including development and stress responses. Despite extensive research on bZIPs in many plants, a comprehensive genome-wide analysis of bZIPs in garlic has yet to be undertaken.

Methods

In this study, we identified and classified 64 AsbZIP genes (AsbZIPs) into 10 subfamilies. A systematic analysis of the evolutionary characteristics of these AsbZIPs, including chromosome location, gene structure, conserved motifs, and gene duplication, was conducted. Furthermore, we also examined the nucleotide diversity, cis-acting elements, and expression profiles of AsbZIPs in various tissues and under different abiotic stresses and hormone treatments.

Results and Discussion

Our findings revealed that gene replication plays a crucial role in the expansion of AsbZIPs, with a minor genetic bottleneck observed during domestication. Moreover, the identification of cis-acting elements suggested potential associations of AsbZIPs with garlic development, hormone, and stress responses. Several AsbZIPs exhibited tissue-preferential and stress/hormone-responsive expression patterns. Additionally, Asa7G01972 and Asa7G01379 were notably differentially expressed under various stresses and hormone treatments. Subsequent yeast two-hybridization and yeast induction experiments validated their interactions with Asa1G01577, a homologue of ABI5, reinforcing their importance in hormone and abiotic stress responses. This study unveiled the characteristics of the AsbZIP superfamily and lays a solid foundation for further functional analysis of AsbZIP in garlic.

A single viral amino acid shapes the root system architecture of a plant host upon virus infection

BACKGROUND

Grapevine fanleaf virus (GFLV) is one of the most detrimental viral pathogens of grapevines worldwide but no information is available on its effect on the root system architecture (RSA) of plant hosts. We used two wildtype GFLV strains and their single amino acid mutants to assess RSA traits in infected Nicotiana benthamiana and evaluate transcriptomic changes in host root gene expression in replicated time course 3'RNA-Seq experiments. Mutations targeted the multi-functional GFLV-encoded protein 1EPol*/Sd, a putative RNA-dependent RNA polymerase and determinant of foliar symptoms in N. benthamiana plants.

RESULTS

Plant infection with wildtype GFLV strain GHu and mutant GFLV strain F13 1EPol G802K, both carrying a lysine in position 802 of protein 1EPol*/Sd, resulted in a significantly lower number of root tips (-30%), and a significantly increased average root diameter (+ 20%) at 17 days post inoculation (dpi) in comparison with roots of mock inoculated plants. In contrast, the RSA of plants infected with wildtype GFLV strain F13 and mutant GFLV strain GHu 1EPol K802G, both carrying a glycine in position 802 of protein 1EPol*/Sd, resembled that of mock inoculated plants. Modifications of RSA traits were not associated with GFLV titer. Root tissue transcriptome analysis at 17 dpi indicated dysregulation of pattern recognition receptors, plant hormones, RNA silencing, and genes related to the production of reactive oxygen species (ROS). For wildtype GFLV strain GHu, RSA modifications were correlated with an abundant accumulation of ROS in the pericycle of primary roots at 7 dpi and the duration of vein clearing symptom expression in apical leaves. Dysegulation of a hypersensitive response was an overarching gene ontology found through enrichment analyses of 3'RNA-Seq data.

CONCLUSIONS

Our findings revealed the causative role of lysine in position 802 of protein 1EPol*/Sd in a novel RSA phenotype during viral infection and documented GFLV-N. benthamiana interactions at the root level based on (i) antiviral response, (ii) receptor mediated production of ROS, and (iii) hormone regulation. A correlation between above and below ground symptoms was reported for the first time in plants infected with wildtype GFLV strain GHu. Further work is warranted to test whether the modified RSA of a plant host might impact GFLV acquisition and transmission by the ectoparasitic dagger nematode Xiphinema index.

Bacterial Spermosphere Inoculants Alter N. benthamiana-Plant Physiology and Host Bacterial Microbiome

In this study, we investigated the interplay between the spermosphere inoculum, host plant physiology, and endophytic compartment (EC) microbial community. Using 16S ribosomal RNA gene sequencing of root, stem, and leaf endophytic compartment communities, we established a baseline microbiome for Nicotiana sp. Phenotypic differences were observed due to the addition of some bacterial inoculants, correlated with endogenous auxin loads using transgenic plants expressing the auxin reporter pB-GFP::P87. When applied as spermosphere inoculants, select bacteria were found to create reproducible variation within the root EC microbiome and, more systematically, the host plant physiology. Our findings support the assertion that the spermosphere of plants is a zone that can influence the EC microbiome when applied in a greenhouse setting.

A novel ilarvirus protein CP-RT is expressed via stop codon readthrough and suppresses RDR6-dependent RNA silencing

Ilarviruses are a relatively understudied but important group of plant RNA viruses that includes a number of crop pathogens. Their genomes comprise three RNA segments encoding two replicase subunits, movement protein, coat protein (CP), and (in some ilarvirus subgroups) a protein that suppresses RNA silencing. Here we report that, in many ilarviruses, RNA3 encodes an additional protein (termed CP-RT) as a result of ribosomal readthrough of the CP stop codon into a short downstream readthrough (RT) ORF. Using asparagus virus 2 as a model, we find that CP-RT is expressed in planta where it functions as a weak suppressor of RNA silencing. CP-RT expression is essential for persistent systemic infection in leaves and shoot apical meristem. CP-RT function is dependent on a putative zinc-finger motif within RT. Replacing the asparagus virus 2 RT with the RT of an ilarvirus from a different subgroup restored the ability to establish persistent infection. These findings open up a new avenue for research on ilarvirus silencing suppression, persistent meristem invasion and vertical transmission.

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.

StMAPKK5 Positively Regulates Response to Drought and Salt Stress in Potato

MAPKKs, as one of the main members of the mitogen-activated protein kinase (MAPK) cascade pathway, are located in the middle of the cascade and are involved in many physiological processes of plant growth and development, as well as stress tolerance. Previous studies have found that StMAPKK5 is responsive to drought and salt stress. To further investigate the function and regulatory mechanism of StMAPKK5 in potato stress response, potato variety 'Atlantic' was subjected to drought and NaCl treatments, and the expression of the StMAPKK5 gene was detected by qRT-PCR. StMAPKK5 overexpression and RNA interference-mediated StMAPKK5 knockdown potato plants were constructed. The relative water content, superoxide dismutase (SOD), catalase (CAT), and peroxidase (POD) activities, as well as proline (Pro) and malondialdehyde (MDA) contents of plant leaves, were also assayed under drought and NaCl stress. The StMAPKK5 interacting proteins were identified and validated by yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC). The results showed that the expression of StMAPKK5 was significantly up-regulated under drought and NaCl stress conditions. The StMAPKK5 protein was localized in the nucleus, cytoplasm, and cell membrane. The expression of StMAPKK5 affected the relative water content, the enzymatic activities of SOD, CAT, and POD, and the proline and MDA contents of potatoes under drought and salt stress conditions. These results suggest that StMAPKK5 plays a significant role in regulating drought and salt tolerance in potato crop. Yeast two-hybrid (Y2H) screening identified four interacting proteins: StMYB19, StZFP8, StPUB-like, and StSKIP19. BiFC confirmed the authenticity of the interactions. These findings suggest that StMAPKK5 is crucial for potato growth, development, and response to adversity.

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.

Guide RNA scaffold variants enabled easy cloning of large gRNA cluster for multiplexed gene editing

Cas9 protein-mediated gene editing has revolutionized genetic manipulation in most organisms. There are many cases where multiplexed gene editing is needed. Cas9 is capable of multiplex gene editing when expressed with multiple guide RNAs. Conventional cloning methods for multiplexed gene editing vector is not efficient due to repeated use of a single-guide RNA scaffold and inefficient ligation. In this study, we conducted structure-guided mutagenesis and random mutagenesis on the original sgRNA scaffold and identified a large number of functional sgRNA scaffold variants. With these scaffold variants and different tRNAs, fusion polymerase chain reaction protocol was developed to rapidly synthesize spacer-scaffold-tRNA-spacer units with up to 9 targets. In conjunction with golden gate cloning, gene editing vectors with up to 24 target sites were efficiently cloned in one-step cloning. One such gene editing vector targeting 12 genes in tomato were tested in stable transformation and 10 out of the 12 genes were found mutated in a single transgenic line. To facilitate the application of multiplexed gene editing using these scaffold variants and tRNAs from different species, a webserver was created to generate primer sets and provide template sequences for the synthesis of large sgRNA expression units based on the user-supplied target sequences and species.

Genetic Cause of Hybrid Lethality Observed in Reciprocal Interspecific Crosses between Nicotiana simulans and N. tabacum

Hybrid lethality, a type of postzygotic reproductive isolation, is an obstacle to wide hybridization breeding. Here, we report the hybrid lethality that was observed in crosses between the cultivated tobacco, Nicotiana tabacum (section Nicotiana), and the wild tobacco species, Nicotiana simulans (section Suaveolentes). Reciprocal hybrid seedlings were inviable at 28 °C, and the lethality was characterized by browning of the hypocotyl and roots, suggesting that hybrid lethality is due to the interaction of nuclear genomes derived from each parental species, and not to a cytoplasmic effect. Hybrid lethality was temperature-sensitive and suppressed at 36 °C. However, when hybrid seedlings cultured at 36 °C were transferred to 28 °C, all of them showed hybrid lethality. After crossing between an N. tabacum monosomic line missing one copy of the Q chromosome and N. simulans, hybrid seedlings with or without the Q chromosome were inviable and viable, respectively. These results indicated that gene(s) on the Q chromosome are responsible for hybrid lethality and also suggested that N. simulans has the same allele at the Hybrid Lethality A1 (HLA1) locus responsible for hybrid lethality as other species in the section Suaveolentes. Haplotype analysis around the HLA1 locus suggested that there are at least six and two haplotypes containing Hla1-1 and hla1-2 alleles, respectively, in the section Suaveolentes.

Identification and characterization of GLYCEROLIPASE A1 for wound-triggered JA biosynthesis in Nicotiana benthamiana leaves

Although many important discoveries have been made regarding the jasmonate signaling pathway, how jasmonate biosynthesis is initiated is still a major unanswered question in the field. Previous evidences suggest that jasmonate biosynthesis is limited by the availability of fatty acid precursor, such as ⍺-linolenic acid (⍺-LA). This indicates that the lipase responsible for releasing α-LA in the chloroplast, where early steps of jasmonate biosynthesis take place, is the key initial step in the jasmonate biosynthetic pathway. Nicotiana benthamiana glycerol lipase A1 (NbGLA1) is homologous to N. attenuata GLA1 (NaGLA1) which has been reported to be a major lipase in leaves for jasmonate biosynthesis. NbGLA1 was studied for its potential usefulness in a species that is more common in laboratories. Virus-induced gene silencing of both NbGLA1 and NbGLA2, another homolog, resulted in more than 80% reduction in jasmonic acid (JA) biosynthesis in wounded leaves. Overexpression of NbGLA1 utilizing an inducible vector system failed to increase JA, indicating that transcriptional induction of NbGLA1 is insufficient to trigger JA biosynthesis. However, co-treatment with wounding in addition to NbGLA1 induction increased JA accumulation several fold higher than the gene expression or wounding alone, indicating an enhancement of the enzyme activity by wounding. Domain-deletion of a 126-bp C-terminal region hypothesized to have regulatory roles increased NbGLA1-induced JA level. Together, the data show NbGLA1 to be a major lipase for wound-induced JA biosynthesis in N. benthamiana leaves and demonstrate the use of inducible promoter-driven construct of NbGLA1 in conjunction with its transient expression in N. benthamiana as a useful system to study its protein function.

Methyltransferase-like (METTL) homologues participate in Nicotiana benthamiana antiviral responses

Methyltransferase (MTase) enzymes catalyze the addition of a methyl group to a variety of biological substrates. MTase-like (METTL) proteins are Class I MTases whose enzymatic activities contribute to the epigenetic and epitranscriptomic regulation of multiple cellular processes. N6-adenosine methylation (m6A) is a common chemical modification of eukaryotic and viral RNA whose abundance is jointly regulated by MTases and METTLs, demethylases, and m6A binding proteins. m6A affects various cellular processes including RNA degradation, post-transcriptional processing, and antiviral immunity. Here, we used Nicotiana benthamiana and plum pox virus (PPV), an RNA virus of the Potyviridae family, to investigated the roles of MTases in plant-virus interaction. RNA sequencing analysis identified MTase transcripts that are differentially expressed during PPV infection; among these, accumulation of a METTL gene was significantly downregulated. Two N. benthamiana METTL transcripts (NbMETTL1 and NbMETTL2) were cloned and further characterized. Sequence and structural analyses of the two encoded proteins identified a conserved S-adenosyl methionine (SAM) binding domain, showing they are SAM-dependent MTases phylogenetically related to human METTL16 and Arabidopsis thaliana FIONA1. Overexpression of NbMETTL1 and NbMETTL2 caused a decrease of PPV accumulation. In sum, our results indicate that METTL homologues participate in plant antiviral responses.

The advent of plant cells in bioreactors

Ever since agriculture started, plants have been bred to obtain better yields, better fruits, or sustainable products under uncertain biotic and abiotic conditions. However, a new way to obtain products from plant cells emerged with the development of recombinant DNA technologies. This led to the possibility of producing exogenous molecules in plants. Furthermore, plant chemodiversity has been the main source of pharmacological molecules, opening a field of plant biotechnology directed to produce high quality plant metabolites. The need for different products by the pharma, cosmetics agriculture and food industry has pushed again to develop new procedures. These include cell production in bioreactors. While plant tissue and cell culture are an established technology, beginning over a hundred years ago, plant cell cultures have shown little impact in biotechnology projects, compared to bacterial, yeasts or animal cells. In this review we address the different types of bioreactors that are currently used for plant cell production and their usage for quality biomolecule production. We make an overview of Nicotiana tabacum, Nicotiana benthamiana, Oryza sativa, Daucus carota, Vitis vinifera and Physcomitrium patens as well-established models for plant cell culture, and some species used to obtain important metabolites, with an insight into the type of bioreactor and production protocols.

A plasma membrane nucleotide-binding leucine-rich repeat receptor mediates the recognition of the Ralstonia pseudosolanacearum effector RipY in Nicotiana benthamiana

Bacterial wilt disease caused by several Ralstonia species is one of the most destructive diseases in Solanaceae crops. Only a few functional resistance genes against bacterial wilt have been cloned to date. Here, we show that the broadly conserved type III secreted effector RipY is recognized by the Nicotiana benthamiana immune system, leading to cell death induction, induction of defense-related gene expression, and restriction of bacterial pathogen growth. Using a multiplexed virus-induced gene-silencing-based N. benthamiana nucleotide-binding and leucine-rich repeat receptor (NbNLR) library, we identified a coiled-coil (CC) nucleotide-binding and leucine-rich repeat receptor (CNL) required for recognition of RipY, which we named RESISTANCE TO RALSTONIA SOLANACEARUM RIPY (RRS-Y). Genetic complementation assays in RRS-Y-silenced plants and stable rrs-y knockout mutants demonstrated that RRS-Y is sufficient to activate RipY-induced cell death and RipY-induced immunity to Ralstonia pseudosolanacearum. RRS-Y function is dependent on the phosphate-binding loop motif of the nucleotide-binding domain but independent of the characterized signaling components ENHANCED DISEASE SUSCEPTIBILITY 1, ACTIVATED DISEASE RESISTANCE 1, and N REQUIREMENT GENE 1 and the NLR helpers NB-LRR REQUIRED FOR HR-ASSOCIATED CELL DEATH-2, -3, and -4 in N. benthamiana. We further show that RRS-Y localization at the plasma membrane is mediated by two cysteine residues in the CC domain and is required for RipY recognition. RRS-Y also broadly recognizes RipY homologs across Ralstonia species. Lastly, we show that the C-terminal region of RipY is indispensable for RRS-Y activation. Together, our findings provide an additional effector/receptor pair system to deepen our understanding of CNL activation in plants.

PSTVd infection in Nicotiana benthamiana plants has a minor yet detectable effect on CG methylation

Viroids are small circular RNAs infecting a wide range of plants. They do not code for any protein or peptide and therefore rely on their structure for their biological cycle. Observed phenotypes of viroid infected plants are thought to occur through changes at the transcriptional/translational level of the host. A mechanism involved in such changes is RNA-directed DNA methylation (RdDM). Till today, there are contradictory works about viroids interference of RdDM. In this study, we investigated the epigenetic effect of viroid infection in Nicotiana benthamiana plants. Using potato spindle tuber viroid (PSTVd) as the triggering pathogen and via bioinformatic analyses, we identified endogenous gene promoters and transposable elements targeted by 24 nt host siRNAs that differentially accumulated in PSTVd-infected and healthy plants. The methylation status of these targets was evaluated following digestion with methylation-sensitive restriction enzymes coupled with PCR amplification, and bisulfite sequencing. In addition, we used Methylation Sensitive Amplification Polymorphism (MSAP) followed by sequencing (MSAP-seq) to study genomic DNA methylation of 5-methylcytosine (5mC) in CG sites upon viroid infection. In this study we identified a limited number of target loci differentially methylated upon PSTVd infection. These results enhance our understanding of the epigenetic host changes as a result of pospiviroid infection.

NLR receptors in plant immunity: making sense of the alphabet soup

Plants coordinately use cell-surface and intracellular immune receptors to perceive pathogens and mount an immune response. Intracellular events of pathogen recognition are largely mediated by immune receptors of the nucleotide binding and leucine rich-repeat (NLR) classes. Upon pathogen perception, NLRs trigger a potent broad-spectrum immune reaction, usually accompanied by a form of programmed cell death termed the hypersensitive response. Some plant NLRs act as multifunctional singleton receptors which combine pathogen detection and immune signaling. However, NLRs can also function in higher order pairs and networks of functionally specialized interconnected receptors. In this article, we cover the basic aspects of plant NLR biology with an emphasis on NLR networks. We highlight some of the recent advances in NLR structure, function, and activation and discuss emerging topics such as modulator NLRs, pathogen suppression of NLRs, and NLR bioengineering. Multi-disciplinary approaches are required to disentangle how these NLR immune receptor pairs and networks function and evolve. Answering these questions holds the potential to deepen our understanding of the plant immune system and unlock a new era of disease resistance breeding.

An improved Nicotiana benthamiana bioproduction chassis provides novel insights into nicotine biosynthesis

The model plant Nicotiana benthamiana is an increasingly attractive organism for the production of high-value, biologically active molecules. However, N. benthamiana accumulates high levels of pyridine alkaloids, in particular nicotine, which complicates the downstream purification processes. Here, we report a new assembly of the N. benthamiana genome as well as the generation of low-nicotine lines by CRISPR/Cas9-based inactivation of berberine bridge enzyme-like proteins (BBLs). Triple as well as quintuple mutants accumulated three to four times less nicotine than the respective control lines. The availability of lines without functional BBLs allowed us to probe their catalytic role in nicotine biosynthesis, which has remained obscure. Notably, chiral analysis revealed that the enantiomeric purity of nicotine was fully lost in the quintuple mutants. In addition, precursor feeding experiments showed that these mutants cannot facilitate the specific loss of C6 hydrogen that characterizes natural nicotine biosynthesis. Our work delivers an improved N. benthamiana chassis for bioproduction and uncovers the crucial role of BBLs in the stereoselectivity of nicotine biosynthesis.

Non-specific LIPID TRANSFER PROTEIN 1 enhances immunity against tobacco mosaic virus in Nicotiana benthamiana

Plant non-specific lipid transfer proteins (nsLTPs) are small, cysteine-rich proteins that play significant roles in biotic and abiotic stress responses; however, the molecular mechanism of their functions against viral infections remains unclear. In this study, we employed virus-induced gene-silencing and transgenic overexpression to functionally analyse a type-I nsLTP in Nicotiana benthamiana, NbLTP1, in the immunity response against tobacco mosaic virus (TMV). NbLTP1 was inducible by TMV infection, and its silencing increased TMV-induced oxidative damage and the production of reactive oxygen species (ROS), compromised local and systemic resistance to TMV, and inactivated the biosynthesis of salicylic acid (SA) and its downstream signaling pathway. The effects of NbLTP1-silencing were partially restored by application of exogenous SA. Overexpressing NbLTP1 activated genes related to ROS scavenging to increase cell membrane stability and maintain redox homeostasis, confirming that an early ROS burst followed by ROS suppression at the later phases of pathogenesis is essential for resistance to TMV infection. The cell-wall localization of NbLTP1 was beneficial to viral resistance. Overall, our results showed that NbLTP1 positively regulates plant immunity against viral infection through up-regulating SA biosynthesis and its downstream signaling component, NONEXPRESSOR OF PATHOGENESIS-RELATED 1 (NPR1), which in turn activates pathogenesis-related genes, and by suppressing ROS accumulation at the later phases of viral pathogenesis.

The Phytophthora nucleolar effector Pi23226 targets host ribosome biogenesis to induce necrotrophic cell death

Pathogen effectors target diverse subcellular organelles to manipulate the plant immune system. Although the nucleolus has emerged as a stress marker and several effectors are localized in the nucleolus, the roles of nucleolar-targeted effectors remain elusive. In this study, we showed that Phytophthora infestans infection of Nicotiana benthamiana results in nucleolar inflation during the transition from the biotrophic to the necrotrophic phase. Multiple P. infestans effectors were localized in the nucleolus: Pi23226 induced cell death in N. benthamiana and nucleolar inflation similar to that observed in the necrotrophic stage of infection, whereas its homolog Pi23015 and a deletion mutant (Pi23226ΔC) did not induce cell death or affect nucleolar size. RNA immunoprecipitation and individual-nucleotide-resolution UV crosslinking and immunoprecipitation sequencing analysis indicated that Pi23226 bound to the 3' end of 25S rRNA precursors, resulting in accumulation of unprocessed 27S pre-rRNAs. The nucleolar stress marker NAC082 was strongly upregulated under Pi23226-expressing conditions. Pi23226 subsequently inhibited global protein translation in host cells by interacting with ribosomes. Pi23226 enhanced P. infestans pathogenicity, indicating that Pi23226-induced ribosome malfunction and cell death were beneficial for pathogenesis in the host. Our results provide evidence for the molecular mechanism underlying RNA-binding effector activity in host ribosome biogenesis and lead to new insights into the nucleolar action of effectors in pathogenesis.

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

Nicotiana benthamiana is an invaluable model plant and biotechnology platform with a ~3 Gb allotetraploid genome. To further improve its usefulness and versatility, we have produced high-quality chromosome-level genome assemblies, coupled with transcriptome, epigenome, microRNA and transposable element datasets, for the ubiquitously used LAB strain and a related wild accession, QLD. In addition, single nucleotide polymorphism maps have been produced for a further two laboratory strains and four wild accessions. Despite the loss of five chromosomes from the ancestral tetraploid, expansion of intergenic regions, widespread segmental allopolyploidy, advanced diploidization and evidence of recent bursts of Copia pseudovirus (Copia) mobility not seen in other Nicotiana genomes, the two subgenomes of N. benthamiana show large regions of synteny across the Solanaceae. LAB and QLD have many genetic, metabolic and phenotypic differences, including disparate RNA interference responses, but are highly interfertile and amenable to genome editing and both transient and stable transformation. The LAB/QLD combination has the potential to be as useful as the Columbia-0/Landsberg errecta partnership, utilized from the early pioneering days of Arabidopsis genomics to today.

Ptr1 and ZAR1 immune receptors confer overlapping and distinct bacterial pathogen effector specificities

Some nucleotide-binding and leucine-rich repeat receptors (NLRs) indirectly detect pathogen effectors by monitoring their host targets. In Arabidopsis thaliana, RIN4 is targeted by multiple sequence-unrelated effectors and activates immune responses mediated by RPM1 and RPS2. These effectors trigger cell death in Nicotiana benthamiana, but the corresponding NLRs have yet not been identified. To identify N. benthamiana NLRs (NbNLRs) that recognize Arabidopsis RIN4-targeting effectors, we conducted a rapid reverse genetic screen using an NbNLR VIGS library. We identified that the N. benthamiana homolog of Ptr1 (Pseudomonas tomato race 1) recognizes the Pseudomonas effectors AvrRpt2, AvrRpm1, and AvrB. We demonstrated that recognition of the Xanthomonas effector AvrBsT and the Pseudomonas effector HopZ5 is conferred independently by the N. benthamiana homolog of Ptr1 and ZAR1. Interestingly, the recognition of HopZ5 and AvrBsT is contributed unequally by Ptr1 and ZAR1 in N. benthamiana and Capsicum annuum. In addition, we showed that the RLCK XII family protein JIM2 is required for the NbZAR1-dependent recognition of AvrBsT and HopZ5. The recognition of sequence-unrelated effectors by NbPtr1 and NbZAR1 provides an additional example of convergently evolved effector recognition. Identification of key components involved in Ptr1 and ZAR1-mediated immunity could reveal unique mechanisms of expanded effector recognition.

NLRs derepress MED10b- and MED7-mediated repression of jasmonate-dependent transcription to activate immunity

Plant intracellular nucleotide-binding domain, leucine-rich repeat-containing receptors (NLRs) activate a robust immune response upon detection of pathogen effectors. How NLRs induce downstream immune defense genes remains poorly understood. The Mediator complex plays a central role in transducing signals from gene-specific transcription factors to the transcription machinery for gene transcription/activation. In this study, we demonstrate that MED10b and MED7 of the Mediator complex mediate jasmonate-dependent transcription repression, and coiled-coil NLRs (CNLs) in Solanaceae modulate MED10b/MED7 to activate immunity. Using the tomato CNL Sw-5b, which confers resistance to tospovirus, as a model, we found that the CC domain of Sw-5b directly interacts with MED10b. Knockout/down of MED10b and other subunits including MED7 of the middle module of Mediator activates plant defense against tospovirus. MED10b was found to directly interact with MED7, and MED7 directly interacts with JAZ proteins, which function as transcriptional repressors of jasmonic acid (JA) signaling. MED10b-MED7-JAZ together can strongly repress the expression of JA-responsive genes. The activated Sw-5b CC interferes with the interaction between MED10b and MED7, leading to the activation of JA-dependent defense signaling against tospovirus. Furthermore, we found that CC domains of various other CNLs including helper NLR NRCs from Solanaceae modulate MED10b/MED7 to activate defense against different pathogens. Together, our findings reveal that MED10b/MED7 serve as a previously unknown repressor of jasmonate-dependent transcription repression and are modulated by diverse CNLs in Solanaceae to activate the JA-specific defense pathways.

Phylogenomic discovery of deleterious mutations facilitates hybrid potato breeding

Hybrid potato breeding will transform the crop from a clonally propagated tetraploid to a seed-reproducing diploid. Historical accumulation of deleterious mutations in potato genomes has hindered the development of elite inbred lines and hybrids. Utilizing a whole-genome phylogeny of 92 Solanaceae and its sister clade species, we employ an evolutionary strategy to identify deleterious mutations. The deep phylogeny reveals the genome-wide landscape of highly constrained sites, comprising ∼2.4% of the genome. Based on a diploid potato diversity panel, we infer 367,499 deleterious variants, of which 50% occur at non-coding and 15% at synonymous sites. Counterintuitively, diploid lines with relatively high homozygous deleterious burden can be better starting material for inbred-line development, despite showing less vigorous growth. Inclusion of inferred deleterious mutations increases genomic-prediction accuracy for yield by 24.7%. Our study generates insights into the genome-wide incidence and properties of deleterious mutations and their far-reaching consequences for breeding.

Horizontally transferred genes as RNA interference targets for aphid and whitefly control

RNA interference (RNAi)-based technologies are starting to be commercialized as a new approach for agricultural pest control. Horizontally transferred genes (HTGs), which have been transferred into insect genomes from viruses, bacteria, fungi or plants, are attractive targets for RNAi-mediated pest control. HTGs are often unique to a specific insect family or even genus, making it unlikely that RNAi constructs targeting such genes will have negative effects on ladybugs, lacewings and other beneficial predatory insect species. In this study, we sequenced the genome of a red, tobacco-adapted isolate of Myzus persicae (green peach aphid) and bioinformatically identified 30 HTGs. We then used plant-mediated virus-induced gene silencing (VIGS) to show that several HTGs of bacterial and plant origin are important for aphid growth and/or survival. Silencing the expression of fungal-origin HTGs did not affect aphid survivorship but decreased aphid reproduction. Importantly, although there was uptake of plant-expressed RNA by Coccinella septempunctata (seven-spotted ladybugs) via the aphids that they consumed, we did not observe negative effects on ladybugs from aphid-targeted VIGS constructs. To demonstrate that this approach is more broadly applicable, we also targeted five Bemisia tabaci (whitefly) HTGs using VIGS and demonstrated that knockdown of some of these genes affected whitefly survival. As functional HTGs have been identified in the genomes of numerous pest species, we propose that these HTGs should be explored further as efficient and safe targets for control of insect pests using plant-mediated RNA interference.

Heritable CRISPR-Cas9 editing of plant genomes using RNA virus vectors

Viral vectors hold enormous potential for genome editing in plants as transient delivery vehicles of CRISPR-Cas components. Here, we describe a protocol to assemble plant viral vectors for single-guide RNA (sgRNA) delivery. The obtained viral constructs are based on compact T-DNA binary vectors of the pLX series and are delivered into Cas9-expressing plants through agroinoculation. This approach allows rapidly assessing sgRNA design for plant genome targeting, as well as the recovery of progeny with heritable mutations at targeted loci. For complete details on the use and execution of this protocol, please refer to Uranga et al. (2021)1 and Aragonés et al. (2022).2.

Genome-wide identification and expression analysis of wall-associated kinase (WAK) and WAK-like kinase gene family in response to tomato yellow leaf curl virus infection in Nicotiana benthamiana

BACKGROUND

Tomato yellow leaf curl virus (TYLCV) is a major monopartite virus in the family Geminiviridae and has caused severe yield losses in tomato and tobacco planting areas worldwide. Wall-associated kinases (WAKs) and WAK-like kinases (WAKLs) are a subfamily of the receptor-like kinase family implicated in cell wall signaling and transmitting extracellular signals to the cytoplasm, thereby regulating plant growth and development and resistance to abiotic and biotic stresses. Recently, many studies on WAK/WAKL family genes have been performed in various plants under different stresses; however, identification and functional survey of the WAK/WAKL gene family of Nicotiana benthamiana have not yet been performed, even though its genome has been sequenced for several years. Therefore, in this study, we aimed to identify the WAK/WAKL gene family in N. benthamiana and explore their possible functions in response to TYLCV infection.

RESULTS

Thirty-eight putative WAK/WAKL genes were identified and named according to their locations in N. benthamiana. Phylogenetic analysis showed that NbWAK/WAKLs are clustered into five groups. The protein motifs and gene structure compositions of NbWAK/WAKLs appear to be highly conserved among the phylogenetic groups. Numerous cis-acting elements involved in phytohormone and/or stress responses were detected in the promoter regions of NbWAK/WAKLs. Moreover, gene expression analysis revealed that most of the NbWAK/WAKLs are expressed in at least one of the examined tissues, suggesting their possible roles in regulating the growth and development of plants. Virus-induced gene silencing and quantitative PCR analyses demonstrated that NbWAK/WAKLs are implicated in regulating the response of N. benthamiana to TYLCV, ten of which were dramatically upregulated in locally or systemically infected leaves of N. benthamiana following TYLCV infection.

CONCLUSIONS

Our study lays an essential base for the further exploration of the potential functions of NbWAK/WAKLs in plant growth and development and response to viral infections in N. benthamiana.

Genome Sequence and Analysis of Nicotiana benthamiana, the Model Plant for Interactions between Organisms

Nicotiana benthamiana is widely used as a model plant for dicotyledonous angiosperms. In fact, the strains used in research are highly susceptible to a wide range of viruses. Accordingly, these strains are subject to plant pathology and plant-microbe interactions. In terms of plant-plant interactions, N. benthamiana is one of the plants that exhibit grafting affinity with plants from different families. Thus, N. benthamiana is a good model for plant biology and has been the subject of genome sequencing analyses for many years. However, N. benthamiana has a complex allopolyploid genome, and its previous reference genome is fragmented into 141,000 scaffolds. As a result, molecular genetic analysis is difficult to perform. To improve this effort, de novo whole-genome assembly was performed in N. benthamiana with Hifi reads, and 1,668 contigs were generated with a total length of 3.1 Gb. The 21 longest scaffolds, regarded as pseudomolecules, contained a 2.8-Gb sequence, occupying 95.6% of the assembled genome. A total of 57,583 high-confidence gene sequences were predicted. Based on a comparison of the genome structures between N. benthamiana and N. tabacum, N. benthamiana was found to have more complex chromosomal rearrangements, reflecting the age of interspecific hybridization. To verify the accuracy of the annotations, the cell wall modification genes involved in grafting were analyzed, which revealed not only the previously indeterminate untranslated region, intron and open reading frame sequences but also the genomic locations of their family genes. Owing to improved genome assembly and annotation, N. benthamiana would increasingly be more widely accessible.

Bringing Plant Immunity to Light: A Genetically Encoded, Bioluminescent Reporter of Pattern-Triggered Immunity in Nicotiana benthamiana

Plants rely on innate immune systems to defend against a wide variety of biotic attackers. Key components of innate immunity include cell-surface pattern-recognition receptors (PRRs), which recognize pest- and pathogen-associated molecular patterns (PAMPs). Unlike other classes of receptors that often have visible cell-death immune outputs upon activation, PRRs generally lack rapid methods for assessing function. Here, we describe a genetically encoded bioluminescent reporter of immune activation by heterologously expressed PRRs in the model organism Nicotiana benthamiana. We characterized N. benthamiana transcriptome changes in response to Agrobacterium tumefaciens and subsequent PAMP treatment to identify pattern-triggered immunity (PTI)-associated marker genes, which were then used to generate promoter-luciferase fusion fungal bioluminescence pathway (FBP) constructs. A reporter construct termed pFBP_2xNbLYS1::LUZ allows for robust detection of PTI activation by heterologously expressed PRRs. Consistent with known PTI signaling pathways, reporter activation by receptor-like protein (RLP) PRRs is dependent on the known adaptor of RLP PRRs, i.e., SOBIR1. The FBP reporter minimizes the amount of labor, reagents, and time needed to assay function of PRRs and displays robust sensitivity at biologically relevant PAMP concentrations, making it ideal for high throughput screens. The tools described in this paper will be powerful for investigations of PRR function and characterization of the structure-function of plant cell-surface receptors. [Formula: see text] The author(s) have dedicated the work to the public domain under the Creative Commons CC0 "No Rights Reserved" license by waiving all of his or her rights to the work worldwide under copyright law, including all related and neighboring rights, to the extent allowed by law, 2023.

The dicot homolog of maize PPR103 carries a C-terminal DYW domain and is required for C-to-U editing of chloroplast RNA transcripts

In plants, cytidine-to-uridine (C-to-U) editing is a crucial step in processing mitochondria and chloroplast-encoded transcripts. This editing requires nuclear-encoded proteins including members of the pentatricopeptide (PPR) family, especially PLS-type proteins carrying the DYW domain. IPI1/emb175/PPR103 is a nuclear gene encoding a PLS-type PPR protein essential for survival in Arabidopsis thaliana and maize. Arabidopsis IPI1 was identified as likely interacting with ISE2, a chloroplast-localized RNA helicase associated with C-to-U RNA editing in Arabidopsis and maize. Notably, while the Arabidopsis and Nicotiana IPI1 homologs possess complete DYW motifs at their C-termini, the maize homolog, ZmPPR103, lacks this triplet of residues which are essential for editing. We examined the function of ISE2 and IPI1 in chloroplast RNA processing in N. benthamiana. A combination of deep sequencing and Sanger sequencing revealed C-to-U editing at 41 sites in 18 transcripts, with 34 sites conserved in the closely related N. tabacum. Virus induced gene silencing of NbISE2 or NbIPI1 led to defective C-to-U revealed that they have overlapping roles at editing a site in the rpoB transcript but have distinct roles in editing other transcripts. This finding contrasts with maize ppr103 mutants that showed no defects in editing. The results indicate that NbISE2 and NbIPI1 are important for C-to-U editing in N. benthamiana chloroplasts, and they may function in a complex to edit specific sites while having antagonistic effects on editing others. That NbIPI1, carrying a DYW domain, is involved in organelle C-to-U RNA editing supports previous work showing that this domain catalyzes RNA editing.

Genome-wide analysis of MYB family in Nicotiana benthamiana and the functional role of the key members in resistance to Bemisia tabaci

MYB transcription factors (TFs) play a key role in plant resistance to abiotic and biotical stresses. However, little is currently known about their involvement in the plant defense to piercing-sucking insects. Here, we studied the MYB TFs that responded to and resisted Bemisia tabaci whitefly in the model plant Nicotiana benthamiana. Firstly, a total of 453 NbMYB TFs in N. benthamiana genome were identified and 182 R2R3-MYB TFs were analyzed for molecular characteristics, phylogenetic analysis, genetic structure, motif composition, and cis-elements. Then, six stress-related NbMYB genes were selected for further study. The expression pattern shows they were highly expressed in mature leaves and intensively induced upon whitefly attack. Combined with bioinformatic analysis, overexpression, β-Glucuronidase (GUS) assay, and virus-induced silencing tests, we determined the transcriptional regulation of these NbMYBs on the genes in lignin biosynthesis and SA-signaling pathways. Meanwhile, we tested the performance of whitefly on plants with increased or silenced NbMYB genes expression and found that NbMYB42, NbMYB107, NbMYB163, and NbMYB423 were resistant to whitefly. Our results contribute to a comprehensive understanding of the MYB TFs in N. benthamiana. Furthermore, our findings will facilitate further studies on the role of MYB TFs in the interaction between plants and piercing-sucking insects.

Analyses of Cullin1 homologs reveal functional redundancy in S-RNase-based self-incompatibility and evolutionary relationships in eudicots

In Petunia (Solanaceae family), self-incompatibility (SI) is regulated by the polymorphic S-locus, which contains the pistil-specific S-RNase and multiple pollen-specific S-Locus F-box (SLF) genes. SLFs assemble into E3 ubiquitin ligase complexes known as Skp1-Cullin1-F-box complexes (SCFSLF). In pollen tubes, these complexes collectively mediate ubiquitination and degradation of all nonself S-RNases, but not self S-RNase, resulting in cross-compatible, but self-incompatible, pollination. Using Petunia inflata, we show that two pollen-expressed Cullin1 (CUL1) proteins, PiCUL1-P and PiCUL1-B, function redundantly in SI. This redundancy is lost in Petunia hybrida, not because of the inability of PhCUL1-B to interact with SSK1, but due to a reduction in the PhCUL1-B transcript level. This is possibly caused by the presence of a DNA transposon in the PhCUL1-B promoter region, which was inherited from Petunia axillaris, one of the parental species of Pe. hybrida. Phylogenetic and syntenic analyses of Cullin genes in various eudicots show that three Solanaceae-specific CUL1 genes share a common origin, with CUL1-P dedicated to S-RNase-related reproductive processes. However, CUL1-B is a dispersed duplicate of CUL1-P present only in Petunia, and not in the other species of the Solanaceae family examined. We suggest that the CUL1s involved (or potentially involved) in the SI response in eudicots share a common origin.

A viral movement protein co-opts endoplasmic reticulum luminal-binding protein and calreticulin to promote intracellular movement

Intracellular movement is an important step for the initial spread of virus in plants during infection. This process requires virus-encoded movement proteins (MPs) and their interaction with host factors. Despite the large number of known host factors involved in the movement of different viruses, little is known about host proteins that interact with one of the MPs encoded by potexviruses, the triple-gene-block protein 3 (TGBp3). The main obstacle lies in the relatively low expression level of potexviral TGBp3 in hosts and the weak or transient nature of interactions. Here, we used TurboID-based proximity labeling to identify the network of proteins directly or indirectly interacting with the TGBp3 of a potexvirus, Bamboo mosaic virus (BaMV). Endoplasmic reticulum (ER) luminal-binding protein 4 and calreticulin 3 of Nicotiana benthamiana (NbBiP4 and NbCRT3, respectively) associated with the functional TGBp3-containing BaMV movement complexes, but not the movement-defective mutant, TGBp3M. Fluorescent microscopy revealed that TGBp3 colocalizes with NbBiP4 or NbCRT3 and the complexes move together along ER networks to cell periphery in N. benthamiana. Loss- and gain-of-function experiments revealed that NbBiP4 or NbCRT3 is required for the efficient spread and accumulation of BaMV in infected leaves. In addition, overexpression of NbBiP4 or NbCRT3 enhanced the targeting of BaMV TGBp1 to plasmodesmata (PD), indicating that NbBiP4 and NbCRT3 interact with TGBp3 to promote the intracellular transport of virion cargo to PD that facilitates virus cell-to-cell movement. Our findings revealed additional roles for NbBiP4 and NbCRT3 in BaMV intracellular movement through ER networks or ER-derived vesicles to PD, which enhances the spread of BaMV in N. benthamiana.

Down, then up: non-parallel genome size changes and a descending chromosome series in a recent radiation of the Australian allotetraploid plant species, Nicotiana section Suaveolentes (Solanaceae)

BACKGROUND AND AIMS

The extent to which genome size and chromosome numbers evolve in concert is little understood, particularly after polyploidy (whole-genome duplication), when a genome returns to a diploid-like condition (diploidization). We study this phenomenon in 46 species of allotetraploid Nicotiana section Suaveolentes (Solanaceae), which formed <6 million years ago and radiated in the arid centre of Australia.

METHODS

We analysed newly assessed genome sizes and chromosome numbers within the context of a restriction site-associated nuclear DNA (RADseq) phylogenetic framework.

KEY RESULTS

RADseq generated a well-supported phylogenetic tree, in which multiple accessions from each species formed unique genetic clusters. Chromosome numbers and genome sizes vary from n = 2x = 15 to 24 and 2.7 to 5.8 pg/1C nucleus, respectively. Decreases in both genome size and chromosome number occur, although neither consistently nor in parallel. Species with the lowest chromosome numbers (n = 15-18) do not possess the smallest genome sizes and, although N. heterantha has retained the ancestral chromosome complement, n = 2x = 24, it nonetheless has the smallest genome size, even smaller than that of the modern representatives of ancestral diploids.

CONCLUSIONS

The results indicate that decreases in genome size and chromosome number occur in parallel down to a chromosome number threshold, n = 20, below which genome size increases, a phenomenon potentially explained by decreasing rates of recombination over fewer chromosomes. We hypothesize that, more generally in plants, major decreases in genome size post-polyploidization take place while chromosome numbers are still high because in these stages elimination of retrotransposons and other repetitive elements is more efficient. Once such major genome size change has been accomplished, then dysploid chromosome reductions take place to reorganize these smaller genomes, producing species with small genomes and low chromosome numbers such as those observed in many annual angiosperms, including Arabidopsis.

Characterization and expression analysis of bHLH transcription factors reveal their putative regulatory effects on nectar spur development in Aquilegia species

Nectar spur is a hollow extension of certain flower parts and shows strikingly diverse size and shape in Aquilegia. Nectar spur development is involved in cell division and expansion processes. The basic helix-loop-helix (bHLH) transcription factors (TFs) control a diversity of organ morphogenesis, including cell division and cell expansion processes. However, the role of bHLH genes in nectar spur development in Aquilegia is mainly unknown. We conducted a genome-wide identification of the bHLH gene family in Aquilegia to determine structural characteristics and phylogenetic relationships, and to analyze expression profiles of these genes during the development of nectar spur in spurless and spurred species. A total of 120 AqbHLH genes were identified from the Aquilegia coerulea genome. The phylogenetic tree showed that AqbHLH proteins were divided into 15 subfamilies, among which S7 and S8 subfamilies occurred marked expansion. The AqbHLH genes in the same clade had similar motif composition and gene structure characteristics. Conserved residue analysis indicated nineteen residues with conservation of more than 50% were found in the four conserved regions. In the upstream sequence of AqbHLH genes, the light-responsive element was the most abundant cis-acting element. Eighteen AqbHLH genes showed syntenic relationships, and eight genes from four syntenic pairs underwent tandem duplications. According to the expression profiling analysis by public RNA-Seq data and qRT-PCR results, five AqbHLH genes, including AqbHLH027, AqbHLH046, AqbHLH082, AqbHLH083 and AqbHLH092, were differentially expressed between different tissues in A. coerulea at early developmental stages, as well as between spurless and spurred Aquilegia species. Of them, AqbHLH046 was not only highly expressed in spur compared with blade, but also showed higher expression levels in spurred species than spurless specie, suggesting it plays an essential role in the development of spur by regulating cell division. This study lays a foundation to investigate the function of AqbHLH genes family in nectar spur development, and has potential implications for speciation and genetic breeding in the genus Aquilegia.

Comparative and expression analyses of AP2/ERF genes reveal copy number expansion and potential functions of ERF genes in Solanaceae

BACKGROUND

The AP2/ERF gene family is a superfamily of transcription factors that are important in the response of plants to abiotic stress and development. However, comprehensive research of the AP2/ERF genes in the Solanaceae family is lacking.

RESULTS

Here, we updated the annotation of AP2/ERF genes in the genomes of eight Solanaceae species, as well as Arabidopsis thaliana and Oryza sativa. We identified 2,195 AP2/ERF genes, of which 368 (17%) were newly identified. Based on phylogenetic analyses, we observed expansion of the copy number of these genes, especially those belonging to specific Ethylene-Responsive Factor (ERF) subgroups of the Solanaceae. From the results of chromosomal location and synteny analyses, we identified that the AP2/ERF genes of the pepper (Capsicum annuum), the tomato (Solanum lycopersicum), and the potato (Solanum tuberosum) belonging to ERF subgroups form a tandem array and most of them are species-specific without orthologs in other species, which has led to differentiation of AP2/ERF gene repertory among Solanaceae. We suggest that these genes mainly emerged through recent gene duplication after the divergence of these species. Transcriptome analyses showed that the genes have a putative function in the response of the pepper and tomato to abiotic stress, especially those in ERF subgroups.

CONCLUSIONS

Our findings will provide comprehensive information on AP2/ERF genes and insights into the structural, evolutionary, and functional understanding of the role of these genes in the Solanaceae.

Targeted Insertion in Nicotiana benthamiana Genomes via Protoplast Regeneration

Insertion of a specific sequence in a targeted region for precise editing is still a major challenge in plants. Current protocols rely on inefficient homology-directed repair or non-homologous end-joining with modified double-stranded oligodeoxyribonucleotides (dsODNs) as donors. We developed a simple protocol that eliminates the need for expensive equipment, chemicals, modifications of donor DNA, and complicated vector construction. The protocol uses polyethylene glycol (PEG)-calcium to deliver low-cost, unmodified single-stranded oligodeoxyribonucleotides (ssODNs) and CRISPR/Cas9 ribonucleoprotein (RNP) complexes into Nicotiana benthamiana protoplasts. Regenerated plants were obtained from edited protoplasts with an editing frequency of up to 50% at the target locus. The inserted sequence was inherited to the next generation; this method thus opens the possibility for the future exploration of genomes by targeted insertion in plants.

Small RNA-seq dataset of wild type and 16C Nicotiana benthamiana leaves sprayed with naked dsRNA using the high-pressure spraying technique

Double-stranded RNA (dsRNA) applications have emerged as promising alternatives to chemical plant pesticides. It has been proposed that the protective effect of dsRNA is mediated by the RNA interference (RNAi) mechanism. Small RNAs (sRNAs) are one of the landmarks of RNAi mechanisms. Two classes of sRNAs appear upon RNAi, triggered by dsRNA: The cleavage products of the dsRNA mapping directly to the dsRNA sequence and the transitive sRNAs mapping to the target transcript outside of the dsRNA sequence. Therefore, the sRNA-seq data obtained from dsRNA-treated plants have been exclusively analysed in the context of the target genes and the outcome has been considered essential to evaluate the underlying mechanism of dsRNA mediated plant protection. Using high-pressure spraying technology (HPST), we have applied a GFP targeting 139bp-long dsRNA on wild type (WT) and GFP expressing (16C) Nicotiana benthamiana plants in biological triplicates. As a control, we applied water with HPST on 16C N. benthamiana. We have acquired sRNA-seq data on the treated and control leaves 5 days post spraying. In this dataset, we have expanded our sRNA-seq analysis from the target GFP transgene sequence to the whole transcriptome of N. benthamiana to provide the community with a resource for the small RNA landscape after high-pressure spraying in 16C and WT samples. Furthermore, we have provided a comparison of sRNA landscape between WT and 16C lines.

Rearranged Endogenized Plant Pararetroviruses as Evidence of Heritable RNA-based Immunity

Eukaryotic genomics frequently revealed historical spontaneous endogenization events of external invading nucleic acids, such as viral elements. In plants, an extensive occurrence of endogenous plant pararetroviruses (EPRVs) is usually believed to endow hosts with an additional layer of internal suppressive weaponry. However, an actual demonstration of this activity remains speculative. We analyzed the EPRV component and accompanying silencing effectors of Solanum lycopersicum, documenting that intronic/intergenic pararetroviral integrations bearing inverted-repeats fuel the plant's RNA-based immune system with suitable transcripts capable of evoking a silencing response. A surprisingly small set of rearrangements explained a substantial fraction of pararetroviral-derived endogenous small-interfering (si)RNAs, enriched in 22-nt forms typically associated with anti-viral post-transcriptional gene silencing. We provide preliminary evidence that such genetic and immunological signals may be found in other species outside the genus Solanum. Based on molecular dating, bioinformatics, and empirical explorations, we propose that homology-dependent silencing emerging from particular immuno-competent rearranged chromosomal areas that constitute an adaptive heritable trans-acting record of past infections, with potential impact against the unlocking of plant latent EPRVs and cognate-free pararetroviruses.

Induction of distinct plant cell death programs by secreted proteins from the wheat pathogen Zymoseptoria tritici

Cell death processes in eukaryotes shape normal development and responses to the environment. For plant-microbe interactions, initiation of host cell death plays an important role in determining disease outcomes. Cell death pathways are frequently initiated following detection of pathogen-derived molecules which can lead to resistance or susceptibility to disease depending on pathogen lifestyle. We previously identified several small secreted proteins (SSPs) from the wheat-infecting fungus Zymoseptoria tritici that induce rapid cell death in Nicotiana benthamiana following Agrobacterium-mediated delivery and expression (agroinfiltration). Here we investigated whether the execution of host cells was mechanistically similar in response to different Z. tritici SSPs. Using RNA sequencing, we found that transient expression of four Z. tritici SSPs led to massive transcriptional reprogramming within 48 h of agroinfiltration. We observed that distinct host gene expression profiles were induced dependent on whether cell death occurs in a cell surface immune receptor-dependent or -independent manner. These gene expression profiles involved differential transcriptional networks mediated by WRKY, NAC and MYB transcription factors. In addition, differential expression of genes belonging to different classes of receptor-like proteins and receptor-like kinases was observed. These data suggest that different Z. tritici SSPs trigger differential transcriptional reprogramming in plant cells.