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.

Single-cell landscape of long and short glandular trichomes in Nicotiana tabacum leaves.

Glandular trichomes (GTs) play a crucial role in plant defenses and the synthesis of secondary metabolites. Understanding the developmental trajectory of GTs is essential for unraveling their functional significance and potential applications. Here we established a comprehensive single-cell atlas of Nicotiana tabacum leaves, a model plant for GT studies. The atlas included a total of 40,433 cells and successfully captured both long GTs (LGTs) and short GTs (SGTs) from Nicotiana leaves. The developmental trajectories of these trichomes were delineated, revealing potential disparities in epidermal development. Comparative analysis of Arabidopsis and Nicotiana trichome development indicated limited similarity between Arabidopsis epidermal non-glandular trichomes and Nicotiana LGTs and SGTs, implying the essentiality of studying the genes directly involved in the development of Nicotiana GTs for a proper and comprehensive understanding of GT biology. Overall, our results provide profound insights into the developmental intricacies of the specialized GTs.

Identification of QTLs associated with yield-related traits and superior genotype prediction using recombinant inbred line population in tobacco.

Tobacco is an economically significant industrial crop and model plant for genetic research, yet little is known about its genetic architecture. Quantitative trait loci (QTL) analysis was performed for six agronomic traits on an F_7 population of 341 genotypes, parents, and F1 plants using 1974 SSR markers across two environments. 31 QTLs contributing single-locus additive effects on 13 linkage groups (LGs) and 6 QTL pairs contributing epistatic effects on 6 LGs, were detected by the QTLNetwork 2.0 which was developed for the mixed-linear-model-based composite interval mapping (MCIM). Notably, 5 QTLs and 1 epistatic QTL pair were found to have pleiotropic effects on some genetically related traits. Moreover, the Broad sense heritability of the detected QTLs ranged from 1.05% to 43.33%, while genotype-by-environment interaction heritability spanned from 27.09% to 56.25%. Based on the results of QTL mapping, the potential superior lines for all or specific environments were designed and evaluated. Five major QTLs were finely dissected based on the tobacco reference genome of K326, and 31 candidate genes were predicted. This study offered new insights into the complicated genetic architecture and QTL resources for efficient breeding design for genetic improvement of agronomic traits in tobacco.

Dissection of genetic architecture of nine hazardous component traits of mainstream smoke in tobacco (Nicotiana tabacum L.).

Tobacco (Nicotiana tabacum L.) use is the leading cause of preventable death, due to deleterious chemical components and smoke from tobacco products, and therefore reducing harmful chemical components in tobacco is one of the crucial tobacco breeding targets. However, due to complexity of tobacco smoke and unavailability of high-density genetic maps, the genetic architecture of representative hazardous smoke has not been fully dissected. The present study aimed to explore the genetic architecture of nine hazardous component traits of mainstream smoke through QTL mapping using 271 recombinant inbred lines (RILs) derived from K326 and Y3 in multiple environments. The analysis of genotype and genotype by environment interaction (GE) revealed substantially greater heritability over 95% contributed mostly by GE interaction effects. We also observed strong genetic correlations among most studied hazardous smoke traits, with the highest correlation coefficient of 0.84 between carbon monoxide and crotonaldehyde. Based on a published high-density genetic map, a total of 19 novel QTLs were detected for eight traits using a full QTL model, of which 17 QTLs showed significant additive effects, six showed significant additive-by-environment interaction effects, and one pair showed significant epistasis-by-environment interaction effect. Bioinformatics analysis of sequence in QTL region predicted six genes as candidates for four traits, of which Nt21g04598.1, Nt21g04600.1, and Nt21g04601.1 had pleiotropic effects on PHE and TAR.

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.

CRISPR/Cas9-mediated seamless gene replacement in protoplasts expands the resistance spectrum to TMV-U1 strain in regenerated Nicotiana tabacum.

CRISPR/Cas-based genome editing is now extensively used in plant breeding and continues to evolve. Most CRISPR/Cas current applications in plants focus on gene knock-outs; however, there is a pressing need for new methods to achieve more efficient delivery of CRISPR components and gene knock-ins to improve agronomic traits of crop cultivars. We report here a genome editing system that combines the advantages of protoplast technologies with recent CRISPR/Cas advances to achieve seamless large fragment insertions in the model Solanaceae plant Nicotiana tabacum. With this system, two resistance-related regions of the N' gene were replaced with homologous fragments from the N'alata gene to confer TMV-U1 resistance in the T0 generation of GMO-free plants. Our study establishes a reliable genome-editing tool for efficient gene modifications and provides a detailed description of the optimization process to assist other researchers adapt this system for their needs.

Construction of a high-density genetic map and dissection of genetic architecture of six agronomic traits in tobacco (Nicotiana tabacum L.).

Tobacco (Nicotiana tabacum L.) is an economic crop and a model organism for studies on plant biology and genetics. A population of 271 recombinant inbred lines (RIL) derived from K326 and Y3, two elite flue-cured tobacco parents, has been constructed to investigate the genetic basis of agronomic traits in tobacco. Six agronomic traits including natural plant height (nPH), natural leaf number (nLN), stem girth (SG), inter-node length (IL), length of the largest leaf (LL) and width of the largest leaf (LW) were measured in seven environments, spanning the period between 2018 and 2021. We firstly developed an integrated SNP-indel-SSR linkage map with 43,301 SNPs, 2,086 indels and 937 SSRs, which contained 7,107 bin markers mapped on 24 LGs and covered 3334.88 cM with an average genetic distance of 0.469cM. Based on this high-density genetic map, a total of 70 novel QTLs were detected for six agronomic traits by a full QTL model using the software QTLNetwork, of which 32 QTLs showed significant additive effects, 18 QTLs showed significant additive-by-environment interaction effects, 17 pairs showed significant additive-by-additive epistatic effects and 13 pairs showed significant epistasis-by-environment interaction effects. In addition to additive effect as a major contributor to genetic variation, both epistasis effects and genotype-by-environment interaction effects played an important role in explaining phenotypic variation for each trait. In particular, qnLN6-1 was detected with considerably large main effect and high heritability ( h a 2 =34.80%). Finally, four genes including Nt16g00284.1, Nt16g00767.1, Nt16g00853.1, Nt16g00877.1 were predicted as pleiotropic candidate genes for five traits.

Transcriptomics and iTRAQ-proteomics analyses provide novel insights into the defense mechanism of black shank disease in tobacco.

Black shank disease caused by Phytophthora nicotianae is one of the most important diseases in tobacco worldwide and can result in a devastating loss in tobacco cultivation. Many efforts have been carried out to identify the chromosome segment from Nicotiana plumbaginifolia containing a resistance locus carrying a gene named Php; however, the Php gene has not been cloned, and knowledge of the potential mechanism of the Php gene in the resistant lines is limited. To further characterize the resistance mechanism of the Php gene, we first used the resistant line "RBST" and the susceptible cultivar "Honghuadajinyuan" (HD) to obtain the near-isogenic line RBS89 containing the Php gene from RBST. RBS89 showed high resistance to black shank disease. Transcriptomic and iTRAQ analyses were applied to explore the potential defense mechanisms in RBS89 plants in comparison with HD plants with or without inoculation. Many differentially expressed genes (DEGs) and proteins were identified, and some pathogenesis-related (PR) proteins were extensively abundant in the RBS89 plants when compared with the HD plants in response to black shank disease. Importantly, overexpression of the PR gene NtPR-1B in HD plants improved the resistance of tobacco plants to black shank disease, indicating that NtPR-1B and Php genes might have similar roles in protecting tobacco from black shank disease. However, the relationship between NtPR-1B and Php genes requires further analysis. Therefore, our study provides valuable information for breeding tobacco cultivars with black shank disease resistance and sheds light on the defense mechanism of black shank disease in tobacco for enhancing Phytophthora resistance in other Solanaceae crops.

Screening of Tobacco Genotypes for Phytophthora nicotianae Resistance.

Black shank, caused by the oomycetes Phytophthora nicotianae, is destructive to tobacco, and this pathogen is highly pathogenic to many solanaceous crops. P. nicotianae is well adapted to high temperatures; therefore, research on this pathogen is gaining importance in agriculture worldwide because of global warming. P. nicotianae-resistant varieties of tobacco plants are commonly screened by inoculation with oat grains colonized by P. nicotianae and monitoring for the disease symptoms. However, it is difficult to quantify the inoculation intensity since accurate inoculation is crucial in this case. This study aimed to develop an efficient and reliable method for evaluating the resistance of tobacco to infection with P. nicotianae. This method has been successfully used to identify resistant varieties, and the inoculation efficiency was confirmed by real-time PCR. The resistance evaluation method presented in this study is efficient and practical for precision breeding, as well as molecular mechanism research.

The complete chloroplast genome of Nicotiana plumbaginifolia.

Nicotiana plumbaginifolia Vivianiis 1802 is an annual herb, native to Mexico and South America. It is one of the most widely distributed tobacco species. As a wild tobacco, N. plumbaginifolia has provided several economically important disease-resistance genes to cultivated tobacco. We assembled the complete chloroplast genome of N. plumbaginifolia. The chloroplast genome is 155,945 bp in length, which includes a large single copy region (86,621 bp), a small single copy region (18,528 bp) and two separated inverted repeat regions (25,398 bp). A total of 117 unique genes were annotated, consisting of 84 protein-coding genes, 29 tRNA genes and 4 rRNA genes. Based on chloroplast genomes of 17 Nicotiana species, phylogenetic analyses indicated that N. plumbaginifolia was closely related to N. suaveolens and N. amplexicaulis.

Genome-wide identification and expression analysis of NtbHLH gene family in tobacco (Nicotiana tabacum) and the role of NtbHLH86 in drought adaptation.

The bHLH transcription factors play pivotal roles in plant growth and development, production of secondary metabolites and responses to various environmental stresses. Although the bHLH genes have been well studied in model plant species, a comprehensive investigation of the bHLH genes is required for tobacco with newly obtained high-quality genome. In the present study, a total of 309 NtbHLH genes were identified and can be divided into 23 subfamilies. The conserved amino acids which are essential for their function were predicted for the NtbHLH proteins. Moreover, the NtbHLH genes were conserved during evolution through analyzing the gene structures and conserved motifs. A total of 265 NtbHLH genes were localized in the 24 tobacco chromosomes while the remained 44 NtbHLH genes were mapped to the scaffolds due to the complexity of tobacco genome. Moreover, transcripts of NtbHLH genes were obviously tissue-specific expressed from the gene-chip data from 23 tobacco tissues, and expressions of 20 random selected NtbHLH genes were further confirmed by quantitative real-time PCR, indicating their potential functions in the plant growth and development. Importantly, overexpressed NtbHLH86 gene confers improve drought tolerance in tobacco indicating that it might be involved in the regulation of drought stress. Therefore, our findings here provide a valuable information on the characterization of NtbHLH genes and further investigation of their functions in tobacco.

Genome-wide association study of leaf chemistry traits in tobacco.

Leaf chemistry traits are some of the key factors influencing tobacco quality, which can be significantly reduced by lower chemical components in cured leaf. To improve tobacco quality through breeding, genetic diversity analysis, population structure analysis, and genome-wide association studies were performed in a panel of 347 tobacco germplasms and the markers associated with five leaf chemistry traits, including total sugar (TS), reducing sugar (RS), total nitrogen (TN), nicotine (NIC), and total potassium (TP) contents were identified. Four groups were classified at a genetic distance of 0.316 by genetic diversity analysis based on coefficient parameter NEI72 using a program NTSYS-pc2.10e, whereas four well-differentiated subpopulations were postulated in the 347 tobacco accessions. A total of 47 target trait-associated SNPs was detected in at least three environments as well as the best linear unbiased predictions (BLUPs) across all environments, among which two, two, four, six, and one highly suggestive associated SNPs were repeatedly detected in all environments and BLUPs for TS, RS, TN, NIC, and TP, respectively. On the basis of the phenotypic effects of the alleles corresponding to suggestive associated SNPs, five tobacco accessions harboring favorable alleles with elite phenotypic performance in leaf chemistry traits were identified. The results could facilitate quality tobacco breeding for higher leaf chemistry trait contents through molecular marker-assisted approaches.

Construction of a high-density genetic map with whole genome sequencing in Nicotiana tabacum L.

Tobacco (Nicotiana tabacum L.) is an essential commercial crop and an ideal model plant for biological mechanism studies. As an allopolyploid species, tobacco harbors a massive and complex genome, which makes the application of molecular markers complicated and challenging. In our study, we performed whole-genome sequencing of an intraspecific recombinant inbred line (RIL) population, a F1 generation and their parents. With the Nicotiana tabacum (K326 cultivar) genome as reference, a total of 45,081 markers were characterized to construct the genetic map, which spanned a genetic distance of 3486.78 cM. Evaluation of a two-dimensional heat map proved the high quality of the genetic map. We utilized these markers to anchor scaffolds and analyzed the ancestral genome origin of linkage groups (LGs). Furthermore, such a high-density genetic map will be applied for quantitative trait locus (QTL) detection, gene localization, genome-wide association studies (GWAS), and marker-assisted breeding in tobacco.

A transcriptomic profile of topping responsive non-coding RNAs in tobacco roots (Nicotiana tabacum).

BACKGROUND: Non-coding RNAs (ncRNAs), including microRNAs (miRNAs), long ncRNAs (lncRNAs) and circular RNAs (circRNAs), accomplish remarkable variety of biological functions. However, the composition of ncRNAs and their interactions with coding RNAs in modulating and controlling of cellular process in plants is largely unknown. Using a diverse group of high-throughput sequencing strategies, the mRNA, miRNA, lncRNA and circRNA compositions of tobacco (Nicotiana tabacum) roots determined and their alteration and potential biological functions in response to topping treatment analyzed. RESULTS: A total of 688 miRNAs, 7423 non-redundant lncRNAs and 12,414 circRNAs were identified, among which, some selected differentially expressed RNAs were verified by quantitative real-time PCR. Using the differentially expressed RNAs, a co-expression network was established that included all four types of RNAs. The number of circRNAs identified were higher than that of miRNAs and lncRNAs, but only two circRNAs were present in the co-expression network. LncRNAs appear to be the most active ncRNAs based on their numbers presented in the co-expression network, but none of them seems to be an eTM (endogenous Target Mimicry) of miRNAs. Integrated with analyses of sequence interaction, several mRNA-circRNA-miRNA interaction networks with a potential role in the regulation of nicotine biosynthesis were uncovered, including a QS-circQS-miR6024 interaction network. In this network miR6024 was significantly down-regulated, while the expression levels of its two targets, circQS and its host gene QS, were sharply increased following the topping treatment. CONCLUSIONS: These results illustrated the transcriptomic profiles of tobacco roots, the organ responsible for nicotine biosynthesis. mRNAs always play the most important roles, while ncRNAs are also expressed extensively for topping treatment response, especially circRNAs are the most activated in the ncRNA pool. These studies also provided insights on the coordinated regulation module of coding and non-coding RNAs in a single plant biological sample. The findings reported here indicate that ncRNAs appear to form interaction complex for the regulation of stress response forming regulation networks with transcripts involved in nicotine biosynthesis in tobacco.

Genome-Wide Identification, Gene Structure and Expression Analysis of the MADS-Box Gene Family Indicate Their Function in the Development of Tobacco (Nicotiana tabacum L.).

MADS-box genes play a pivotal role in various processes, including floral and seed development, controlling flowering time, regulation of fruits ripening, and respond to abiotic and biotic stressors in planta. Tobacco (Nicotiana tabacum) has been widely used as a model plant for analyzing the gene function, however, there has been less information on the regulation of flowering, and the associated genes. In the present study, a total of 168 NtMADS-box genes were identified from tobacco, and their phylogenetic relationship, chromosome locations, and gene structures were further analyzed. NtMADS-box genes can be clustered into four sub-families of Mα, Mγ, MIKC*, and MIKCC. A total of 111 NtMADS-box genes were distributed on 20 chromosomes, and 57 NtMADS-box genes were located on the unanchored scaffolds due to the complex and incomplete assembly of the tobacco genome. Expression profiles of NtMADS-box genes by microarray from 23 different tissues indicated that members in different NtMADS-box gene subfamilies might play specific roles in the growth and flower development, and the transcript levels of 24 NtMADS-box genes were confirmed by quantitative real-time PCR. Importantly, overexpressed NtSOC1/NtMADS133 could promote early flowering and dwarfism in transgenic tobacco plants. Therefore, our findings provide insights on the characterization of NtMADS-box genes to further study their functions in plant development.

Characterization and evolution of gene clusters for terpenoid phytoalexin biosynthesis in tobacco.

MAIN CONCLUSION: The study performed genome-wide identification, characterization and evolution analysis of gene clusters for phytoalexin terpenoid biosynthesis in tobacco, and specifically illustrated ones for capsidiol, an efficient defensive specialized metabolite. Terpenoid phytoalexins play an important role in plant self-defense against pest and pathogen attack. Terpenoid biosynthesis involves terpene synthase and cytochrome P450, which always locate and function as cluster(s). In this study, we performed genome-wide investigation of metabolic gene clusters involved in terpenoid production in tobacco (Nicotiana tabacum). Due to the complexity of the tobacco genome, we modified a published prediction pipeline to reduce the influence of the large number of repeats and to improve the annotation of tobacco genes with respect to their metabolic functions. We identified 1181 metabolic gene clusters with 34 of them potentially being involved in terpenoid biosynthesis. Through integration with transcriptome and metabolic pathway annotation analyses, 3 of the 34 terpenoid biosynthesis-related gene clusters were determined to be high-confidence ones, with 2 involved in biosynthesis of capsidiol, a terpenoid recognized as 1 of the effective resistance compounds in the Nicotiana species. The capsidiol-related gene cluster was conserved in N. sylvestris, N. tomentosiformis and N. attenuate. Our findings demonstrate that phytoalexins in tobacco can arise from operon-like gene clusters, a genomic pattern characterized as being beneficial for rapid stress response, gene co-regulation, co-function and co-heredity.

Characterization of Nicotiana tabacum genotypes possessing deletion mutations that affect potyvirus resistance and the production of trichome exudates.

BACKGROUND: Advances in genomics technologies are making it increasingly feasible to characterize breeding lines that carry traits of agronomic interest. Tobacco germplasm lines that carry loci designated VAM and va have been extensively investigated due to their association with potyvirus resistance (both VAM and va) and defects in leaf surface compounds originating from glandular trichomes (VAM only). Molecular studies and classical genetic analyses are consistent with the model that VAM and va represent deletion mutations in the same chromosomal region. In this study, we used RNA-seq analysis, together with emerging tobacco reference genome sequence information to characterize the genomic regions deleted in tobacco lines containing VAM and va. RESULTS: Tobacco genotypes TI 1406 (VAM), K326-va and K326 (wild type) were analyzed using RNA-seq to generate a list of genes differentially expressed in TI 1406 and K326-va, versus the K326 control. Candidate genes were localized onto tobacco genome scaffolds and validated as being absent in only VAM, or missing in both VAM and va, through PCR analysis. These results enabled the construction of a map that predicted the relative extent of the VAM and va mutations on the distal end of chromosome 21. The RNA-seq analyses lead to the discovery that members of the cembratrienol synthase gene family are deleted in TI 1406. Transformation of TI 1406 with a cembratrienol synthase cDNA, however, did not recover the leaf chemistry phenotype. Common to both TI 1406 and K326-va was the absence of a gene encoding a specific isoform of a eukaryotic translation initiation factor (eiF4E1.S). Transformation experiments showed that ectopic expression of eiF4E1.S is sufficient to restore potyvirus susceptibility in plants possessing either the va or VAM mutant loci. CONCLUSIONS: We have demonstrated the feasibility of using RNA-seq and emerging whole genome sequence resources in tobacco to characterize the VAM and va deletion mutants. These results lead to the discovery of genes underlying some of the phenotypic traits associated with these historically important loci. Additionally, initial size estimations were made for the deleted regions, and dominant markers were developed that are very close to one of the deletion junctions that defines va.

CRISPR/Cas9-Mediated Mutagenesis of Carotenoid Cleavage Dioxygenase 8 (CCD8) in Tobacco Affects Shoot and Root Architecture.

Strigolactones (SLs) are a class of phytohormones that regulate plant architecture. Carotenoid cleavage dioxygenase (CCD) genes are involved in the biosynthesis of SLs and are identified and characterized in many plants. However, the function of CCD genes in tobacco remains poorly understood. In this study, two closely related genes NtCCD8A and NtCCD8B were cloned from tobacco (Nicotiana tabacum L.). The two NtCCD8 genes are orthologues of the tomato (Solanum lycopersicum) carotenoid cleavage dioxygenase 8 (SlCCD8) gene. NtCCD8A and NtCCD8B were primarily expressed in tobacco roots, but low expression levels of these genes were detected in all plant tissues, and their transcript levels significantly increased in response to phosphate limitation. NtCCD8A and NtCCD8B mutations were introduced into tobacco using the CRISPR/Cas9 system and transgenic tobacco lines for both ntccd8 mutant alleles were identified. The ntccd8a and ntccd8b mutant alleles were inactivated by a deletion of three nucleotides and insertion of one nucleotide, respectively, both of which led to the production of premature stop codons. The ntccd8 mutants had increased shoot branching, reduced plant height, increased number of leaves and nodes, and reduced total plant biomass compared to wild-type plants; however, the root-to-shoot ratio was unchanged. In addition, mutant lines had shorter primary roots and more of lateral roots than wild type. These results suggest that NtCCD8 genes are important for changes in tobacco plant architecture.

Transcriptomic profile of tobacco in response to Tomato zonate spot orthotospovirus infection.

BACKGROUND: Tomato zonate spot virus (TZSV), a dominant species of thrips-transmitted orthotospoviruses in Yunnan and Guangxi provinces in China, causes significant loss of yield in lots of crops and is a major threat to incomes of rural families. However, the detailed molecular mechanism of crop disease caused by TZSV remains obscure. METHODS: Next-generation sequencing (NGS)-based transcriptome analysis (RNA-seq) was performed to investigate and compare the gene expression changes in systemic leaves of tobacco upon infection with TZSV and mock-inoculated plants as a control. RESULTS: De novo assembly and analysis of tobacco transcriptome data by RNA-Seq identified 135,395 unigenes. 2102 differentially expressed genes (DEGs) were obtained in tobacco with TZSV infection, among which 1518 DEGs were induced and 584 were repressed. Gene Ontology enrichment analysis revealed that these DEGs were associated with multiple biological functions, including metabolic process, oxidation-reduction process, photosynthesis process, protein kinase activity. The KEGG pathway analysis of these DEGs indicated that pathogenesis caused by TZSV may affect multiple processes including primary and secondary metabolism, photosynthesis and plant-pathogen interactions. CONCLUSION: Our global survey of transcriptional changes in TZSV infected tobacco provides crucial information into the precise molecular mechanisms underlying pathogenesis and symptom development. This is the first report on the relationships in the TZSV-plant interaction using transcriptome analysis. Findings of present study will significantly help enhance our understanding of the complicated mechanisms of plant responses to orthotospoviral infection.

Transcriptomic profile of tobacco in response to Phytophthora nicotianae infection.

Black shank, caused by Phytophthora nicotianae (P. nicotianae), is a serious disease of cultivated tobacco (Nicotiana tabacum) worldwide. The interactions between tobacco and P. nicotianae are complex and the outcomes of the interactions depend on the tobacco genotype, P. nicotianae strain, and environmental conditions. In this study, we used RNA-sequencing (RNA-Seq) to investigate and compare transcriptional changes in the stems of tobacco upon inoculation with P. nicotianae strain race 0. We used two tobacco varieties: RBST (named from resistance to black shank and tobacco mosaic virus), which was resistant to the P. nicotianae strain race 0, and Honghuadajinyuan (HD), which was susceptible to P. nicotianae race 0. Samples were collected 12 and 72-hour post inoculation (hpi). Analysis of differentially expressed genes (DEGs) and significantly enriched GO terms indicated that several basic defense mechanisms were suppressed in both varieties, which included response to wounding (GO: 0009611), and defense response to fungus (GO: 0050832). We also found some genes that may especially be related to mechanisms of resistance in RBST, such as the one encoding a chitinase. These results will provide a valuable resource for understanding the interactions between P. nicotianae and tobacco plants.