Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 28 April 2022
Sec. Plant Proteomics and Protein Structural Biology
This article is part of the Research Topic Application of Proteomics in Agricultural Biotechnology, Volume II View all 7 articles

Proteomics and Co-expression Network Analysis Reveal the Importance of Hub Proteins and Metabolic Pathways in Nicotine Synthesis and Accumulation in Tobacco (Nicotiana tabacum L.)

\r\nZejun Mo,Zejun Mo1,2Lili Duan,Lili Duan1,2Yuanyuan Pu,Yuanyuan Pu1,2Zonglin Tian,Zonglin Tian2,3Yuzhou Ke,Yuzhou Ke2,3Wen Luo,Wen Luo2,3Kai Pi,Kai Pi2,3Ying Huang,Ying Huang2,3Qiong Nie,Qiong Nie2,3Renxiang Liu,*\r\nRenxiang Liu2,3*
  • 1College of Agriculture, Guizhou University, Guiyang, China
  • 2Key Laboratory of Tobacco Quality in Guizhou Province, Guiyang, China
  • 3College of Tobacco, Guizhou University, Guiyang, China

Nicotine is a unique alkaloid present in tobacco that is widely used in cigarettes and in the agricultural, chemical, and pharmaceutical industries. However, the research on nicotine is mostly limited to its synthesis pathways, and only a few studies have explored the effects of other metabolic pathways on nicotine precursors. Regulating the nicotine content in tobacco can greatly promoting the application of nicotine in other fields. In this study, we performed global data-independent acquisition proteomics analysis of four tobacco varieties. Of the four varieties, one had high nicotine content and three had a low nicotine content. A total of 31,259 distinct peptides and 6,018 proteins across two samples were identified. A total of 45 differentially expressed proteins (DEPs) co-existed in the three comparison groups and were mainly involved in the transport and metallic processes of the substances. Most DEPs were enriched in the biosynthesis of secondary metals, glutathione metabolism, carbon metabolism, and glycolysis/gluconeogenesis. In addition, the weighted gene co-expression network analysis identified an expression module closely related to the nicotine content (Brown, r = 0.74, P = 0.006). Gene Ontology annotation and Kyoto Encyclopaedia of Genes and Genomes enrichment analysis showed that the module proteins were mainly involved in the synthesis and metabolism of nicotine precursors such as arginine, ornithine aspartate, proline, and glutathione. The increased levels of these precursors lead to the synthesis and accumulation of nicotine in plants. More importantly, these proteins regulate nicotine synthesis by affecting the formation of putrescine, which is the core intermediate product in nicotine anabolism. Our results provide a reference for tobacco variety selection with a suitable nicotine content and regulation of the nicotine content. Additionally, the results highlight the importance of other precursor metabolism in nicotine synthesis.

Introduction

Plants produce various secondary metabolites that are vital for growth and environmental adaptation (Sanchez and Demain, 2019). Alkaloids are one of the most important secondary metabolites that play an important role in defence against herbivores and insects (Baldwin, 1998; Steppuhn et al., 2004). Nicotine is a major alkaloid that is produced specifically in tobacco, accounting for approximately 90–95% of total alkaloids and 0.6–3.0% of leaf dry weight (Hibi et al., 1994; Sanchez and Demain, 2019). Nicotine has an extremely high utilisation value in various industries (Mukhopadhyay et al., 2004; Tang et al., 2007), such as medicine (Cloëz-Tayarani and Changeux, 2007; Hao et al., 2017), and agriculture (Thurston et al., 1966; Cloëz-Tayarani and Changeux, 2007; Hao et al., 2017; Lewis et al., 2020). Nicotine is used for the commercial production of tobacco, and the nicotine content directly affects the inherent quality, safety, and availability of cigarette products (Nölke et al., 2018). The World Health Organisation has recommended reducing the nicotine content in cigarettes to reduce overall addiction to such products (World Health Organization, 2009). Therefore, monitoring the nicotine content of tobacco is essential for the tobacco industry and multi-purpose utilisation of nicotine.

Nicotine is synthesised only in the roots of tobacco and is then translocated to the aerial parts of the plant via the xylem and finally accumulates in the leaves (Dewey and Xie, 2013; Shoji and Hashimoto, 2013). These processes are affected by various external factors and internal genetic factors such as plant variety (Tian et al., 2018), climate (Jung et al., 2009; Zhang et al., 2019), soil environment, and cultivation measures (Baldwin, 1989; Cane et al., 2005). Studies have suggested that the tobacco nicotine content is regulated by many factors; however, the ability of nicotine accumulation is mainly determined by the genetic characteristics of varieties (Legg et al., 1969; Dewey and Xie, 2013). Thus, the most effective and preferred method for obtaining new materials with an appropriate nicotine content is the application of modified plant genetics, which involves studying nicotine-related genes and their transcription and translation products (Lewis, 2018).

The biosynthesis of nicotine involves the formation of a pyridine ring and a pyrrole ring, as well as their combination (Baldwin, 1999). This process is mediated by many key enzymes, transcription factors, and associated genes. Several studies have shown that ornithine decarboxylase and arginine decarboxylase provide putrescine for the synthesis of nicotine (Marton and Pegg, 1995; Imanishi et al., 1998; Tiburcio and Galston, 2016). The putrescine N-methyltransferase (PMT) is a rate-limiting enzyme involved in the formation of the pyrrole ring in nicotine biosynthesis, which affects the flow rate of putrescine (Chou and Kutchan, 1998). PMT overexpression can increase the nicotine content in leaves (Sato et al., 2011). The transcription level of the N-methylputrescine oxidase (MPO) gene is also the key for restricting nicotine synthesis (Wagner et al., 1986). The gene catalyses the conversion of N-methylputrescine into 4-methylaminobutyl ether (Chou and Kutchan, 1998). Wagner et al. (1986) showed that quinolinic acid phosphoribosyl transferase (QPT) is another rate-limiting enzyme involved in nicotine synthesis; the stronger the QPT activity in root tissues, the higher is the nicotine content of the tobacco plant. NtMYB305a regulates the expression of multiple genes associated with the nicotine metabolism pathway including the PMT gene; regulates nicotine synthesis; and cooperates with the G-box binding protein NtMYC2a to regulate nicotine metabolism in tobacco (Bian et al., 2022). In addition, the inhibition of MYCla and MYClb expressions was reported to decrease the nicotine content in tobacco B-Y-2 cells, indicating that MYCla and MYClb are the positive regulators of nicotine synthesis (Zhang et al., 2012). The discovery of these key factors is of great significance for the directional regulation of nicotine synthesis and accumulation.

To date, many studies have explored the molecular mechanisms by using “omics” methods such as transcriptomics (Lou et al., 2017; Yuan et al., 2021), proteomics (Hu et al., 2017; Hochholdinger et al., 2018), and metabolomics (Kumar et al., 2017; Sharma et al., 2021). By comparing the omics expression characteristics between different samples, the genes controlling various phenotypes can be identified (Jung et al., 2021). Previous studies have mostly investigated the effect of a certain gene or enzyme in the nicotine synthesis pathway; however, none of the studies have investigated the key regulatory factors affecting nicotine synthesis and metabolism from the perspective of proteomics (Wang et al., 2011; Yin et al., 2018). In this study, we performed a comprehensive proteomics analysis of root tips with different nicotine contents. The results of this study provide a basis and reference for obtaining new tobacco varieties with a suitable nicotine content.

Materials and Methods

Plant Material and Analysis of Nicotine Traits

Through many years of repeated experiments, four tobacco varieties were screened out to elucidate the mechanism of nicotine synthesis and accumulation. All varieties were provided by the Institute of Tobacco, the Chinese Academy of Agricultural Sciences. One variety was sun-cured tobacco (Qinggeng) with a high nicotine content, and the other three varieties (Va116, Basma, and K326) were flue-cured tobacco with a low nicotine content (Supplementary Table 1 and Figure 1).

FIGURE 1
www.frontiersin.org

Figure 1. The difference in the nicotine content among four tobacco varieties from different regions (Pingba County, Xixiu District, and Weining County) and years (2018, 2019, and 2020). The average of three biological replicates per material was used for mapping. Error bars represent significant difference of three replicates. Significant differences among the nicotine contents at P < 0.05 and P < 0.01 were determined using the Duncan’s new multiple range test. The lowercase alphabets represent a significant difference (P < 0.05); uppercase alphabets represent an extremely significant difference (P < 0.01).

Using a randomised block design, a field experiment was performed at the tobacco research experimental base of Guizhou University in 2018 (Pingba County), 2019 (Xixiu District), and 2020 (Weining County), respectively. The experiment was performed in triplicate. The intra- and inter-row spacing of plants was 110 and 55 cm. All plants were topped on the same day when more than 50% of the plants flowered (68th day after transplanting). Samples were collected on day 75 after transplantation. Refer to Mo et al. (2021b) for the collection method of experimental samples. The root tip samples were evenly mixed in equal quantities, washed with clear water, and rinsed with phosphate-buffered saline. Then, the samples were kept in different sterilised centrifugation tubes frozen with liquid nitrogen and immediately stored in a refrigerator at −80°C. These root samples were used for subsequent proteomics analysis, enzyme activity measurement. The enzymes in the leaf samples were denatured at 105°C for 30 min, after which the samples were dried at 75°C, ground into powder, bagged, and sealed for storage. According to the method of Shoji et al. (2009), nicotine was separated from the extract of dry leaf samples and analysed through gas chromatography.

Protein Extraction and Mass Spectrometry

The protein samples were prepared using the trichloroacetic acid-acetone precipitation method. The protein concentration was determined using a BCA protein concentration determination kit according to the manufacturer’s instructions (Source: Bio-rad, Specifications and models: 5000201).

Protein Enzymolysis and Mass Spectrometry

Protein enzymolysis was performed according to the method reported by Li et al. (2019) and Xu (2019). Other parameters and operation methods were performed according to a study by Mo et al. (2021b). Mass spectrometry was performed according to the methods described by Li et al. (2017), Cui (2018), and Dong et al. (2019).

Processing and Analysis of Proteomics Data

Data Dependent Acquisition (DDA) data were searched using the programme Maxquant (Maxquant_1.5.3.17). The database was downloaded from tobacco_uniprot, and the iRT peptide fragment sequence was added to the database. The parameters were set according to the protocols described by Fu et al. (2019) and Dong et al. (2019). The original raw files and the search results were exported to the software Spectronaut (Spectronaut Pulsar Xerox 12.0.20491.4) to construct a spectral library. The following software parameters were set: retention time prediction type was set to dynamic iRT; interference on MS2 level correction was enabled; and cross run normalisation was enabled. All the results were filtered by setting a parameter Q value cut-off of 0.01 (equivalent to the false discovery rate <1%).

Bioinformatics Analyses

In order to monitor and evaluate the stability of the system and the reliability of the experimental data, one quality control (QC) sample (generally a mixed sample mix of all samples) was inserted in the sample cohort, and the data consistency of the QC samples inserted throughout the experiment was evaluated. Within-group Pearson correlation analysis was mainly used to assess the QC. The closer the correlation coefficient is to 1, the more stable the experimental system is quality control. A hierarchical clustering algorithm was used to analyse differentially expressed proteins (DEPs) among the compared groups. The package Complex Heatmap R (R Version 3.4) was used to classify the two dimensions of samples and protein expression simultaneously (distance algorithm: Euclid, connection mode: Average linkage) and to generate a hierarchical clustering heat map (Mo et al., 2021b). The Gene Ontology (GO) function of the identified protein was annotated using the software Blast2Go1 (Conesa et al., 2005). The Cluster of Orthologous Groups of proteins (COG) analysis was performed using the OmicShare tools.2 The Fisher’s exact test was used for the GO functional enrichment analysis of DEPs (Liu et al., 2021). The Kyoto Encyclopaedia of Genes and Genomes (KEGG)3 pathway of the target protein set was annotated using the KEGG Automatic Annotation Server (KAAS) (Moriya et al., 2007; Kanehisa et al., 2012).

Weighted Gene Co-expression Network Analysis

The weighted gene co-expression network analysis (WGCNA) algorithm is commonly used for constructing a gene co-expression network. In the gene co-expression networks, genes with common expression in different samples are present in the same gene network, and the co-expression relationship between the genes is generally determined by the expression correlation coefficient between those genes (Duan et al., 2020). A commonly expressed gene module was obtained and then connected with the concerned phenotypic information to determine the relationship between the gene network and the phenotype, as well as the core genes in the network. WGCNA was performed using the R language package (Langfelder and Horvath, 2008).

Determination of Enzyme Activity

Enzyme activity was determined using a double-antibody one-step sandwich enzyme-linked immunosorbent assay, according to the specific operation methods described by Pu (2020).

Statistical Analyses

Duncan’s new multiple range test was used to analyse variations in nicotine contents among the groups (P < 0.05) by using SPSS version 16.0.

Results

Analysis of the Nicotinic Content in Different Genotypes

To understand the differences in nicotine synthesis and accumulation capacity among the four genotypic materials, we determined the nicotine content of tobacco leaves through gas chromatography. The nicotine content of Qinggeng was found to be significantly higher than those of K326, Basma, and Va116; however, no significant difference was observed in the nicotine content between the three varieties K326, Basma, and Va116 (Figure 1). Simultaneously, the difference between the nicotine content of different varieties was consistent at different time points and regions. These results indicated that the nicotine content of the four varieties exhibit real and stable differences at the genetic level.

Inventory of Root Proteins Identified by Data-Independent Acquisition

To elucidate the mechanisms of different nicotine contents, we performed a data-independent acquisition (DIA) proteomics analysis of roots in the maturation stage. By using the protein analysis systems, we identified 31,259 distinct peptides and 6,018 proteins across all twelve samples. The detailed information of the identified peptides and proteins is presented Supplementary Files 1, 2. The coefficient of variation (CV) and Pearson’s correlation coefficient analysis showed that the whole experimental system was stable and that the data were reliable (Figures 2A,B). A total of 3,850 proteins were found to co-exist in all four varieties. We identified 42 proteins only in Qinggeng and 570 proteins co-existing in the other three varieties (Figure 2C).

FIGURE 2
www.frontiersin.org

Figure 2. The protein expression abundance test of four tobacco varieties. (A) The distribution of coefficient of variation (CV) of quality control samples. (B) The correlation analysis of quality control samples. (C) Display of differentially expressed proteins in the three comparison groups.

Proteins with the threshold criteria (P-values < 0.05 and up fold change >1.5 or down fold change <0.67) were considered as DEPs and identified through pairwise comparisons. We identified a total of 382 DEPs, including 304 upregulated and 78 downregulated proteins, by comparing Qinggeng with Basma. Furthermore, by comparing Qinggeng with Va116, we obtained 334 DEPs, including 237 upregulated and 97 downregulated proteins. Moreover, we observed 207 DEPs, including 143 upregulated and 64 downregulated proteins, by comparing Qinggeng with K326 (Figure 3A). The Venn diagram analysis revealed that 45 DEPs (44 increased and 1 decreased in abundance) were common in the three comparisons (Table 1 and Figure 3B). Subsequently, we performed a hierarchical clustering analysis of these 45 DEPs, and the results indicated that the protein expression patterns in Qinggeng were contrary to those in K326, Basma, and Va116 (Figure 3C). We speculate that these specifically expressed proteins might be involved in nicotine synthesis and accumulation processes.

FIGURE 3
www.frontiersin.org

Figure 3. The differentially expressed protein (DEP) profiles. (A) Venn diagram of DEPs between the four tobacco varieties. (B) Venn diagram of DEPs in the comparison groups Qinggeng vs. Basma, Qinggeng vs. Va116 and Qinggeng vs. K326. (C) Hierarchical cluster analysis of DEPs among the varieties.

TABLE 1
www.frontiersin.org

Table 1. Verification of 45 candidate proteins relating to nicotinic content by data-independent acquisition (DIA) proteomics.

Classification of Differentially Expressed Proteins Identified Through Data-Independent Acquisition

To investigate the biological function of these DEPs, we classified the identified DEPs into 13 categories by using the Clusters of Orthologous Groups (COG) database. The largest category was carbohydrate transport and metabolism (5 DEPs), followed by intracellular trafficking, secretion, and vesicular transport (3 DEPs), post-translational modification, protein turnover, and chaperones (3 DEPs); energy production and conversion (3 DEPs); amino acid transport and metabolism (3 DEPs); and general function prediction only (3 DEPs) (Figure 4A). The classification showed that many co-expressed differential proteins were involved in the transport and metabolic processes of the substances. Proteins with these functions may play an important role in the transport of nicotine, so as to promote the transfer and accumulation of nicotine synthesised by roots to leaves.

FIGURE 4
www.frontiersin.org

Figure 4. Clusters of orthologous groups (COG) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of differentially expressed protein (DEPs). (A) COG function classification of 45 co-expressed DEPs. (B) KEGG function classification of 45 co-expressed DEPs.

To understand the functional consequences of the DEPs associated with the nicotine content, we performed KEGG pathway mapping. The annotation results were divided into 11 categories, mainly involving global and overview maps categories, followed by carbohydrate metabolism and amino acid metabolism (Figure 4B). The KEGG terms such as biosynthesis of secondary metabolites, glutathione metabolism, carbon metabolism, and glycolysis/gluconeogenesis were highly enriched in the co-expressed DEPs. Some of the terms were found to be related to the metabolic processes of nicotine such as amino sugar and nucleotide sugar metabolism and phenylpropanoid biosynthesis. The results suggested that these proteins from different physiological and metabolic pathways jointly regulate nicotine accumulation.

Visualisation of the Nicotine Content Hub Protein

Weighted gene co-expression network analysis (WGCNA) is a common systematic network analysis method. Using WGCNA to analyze transcriptomic and genomic data, we can identify trait related gene modules and hub genes and find biomarkers (Lu et al., 2019). To further analyse changes in proteins during nicotine synthesis and identify the hub proteins that affect the nicotine content, we finally selected 4,594 proteins from root tip samples to construct a co-expression module by using the WGCNA package tool. The power Estimate function in the WGCNA software package was used to estimate the optimal power value (Supplementary File 3). When the power value was 12, the independence degree was up to 0.9 and the average connectivity tended to be 0 (Figure 5A). The results showed that the power value provided a scale-free network that met the requirements and contained sufficient information. Therefore, the cluster analysis was performed on the network topology overlap matrix calculated with a power value of 12. Finally, 32 co-expression modules were identified, with different colours representing different modules (Figure 5B). The results of correlation analysis between each protein co-expression module and nicotine content are shown in Figure 5C. The brown module exhibited a significant positive correlation with the nicotine content (r = 0.74, P = 0.006), which could be used as a key module for subsequent analysis.

FIGURE 5
www.frontiersin.org

Figure 5. Weighted gene co-expression network analysis of differentially expressed protein (DEPs) in the four tobacco varieties. (A) Determination of soft threshold (power). Left: Scale-free topology fit index as a function of different soft threshold (power), the line represents that R2 = 0.9; Right: Mean connectivity as a function of different soft threshold (power). (B) Clustering dendrograms of protein and module division, with dissimilarity based on the topological overlap, together with assigned module colours. Overall, 32 co-expression modules were constructed and are shown in different colours. These modules ranged from large to small according to the number of genes included. (C) Module-sample group association analysis. Each row corresponds to a module, labelled with colour as in panel (B), and each column corresponds to a sample group. The colour of each cell at the row-column intersection indicates the correlation coefficient between the module and the sample group.

The GO annotation of the brown module indicated that the biological processes involved in proteins mainly include two aspects, namely transport and metabolism (Figure 6A). The transport processes included protein transport, intracellular protein transport, and vesicle transport, whereas the metabolic processes included the ornithine metabolic process, aspartate family amino acid metabolic process, and secondary metabolic process. In terms of cellular components, proteins were mainly associated with the plastid, vacuole, and vesicle composition. Moreover, some proteins were associated with the molecular functions such as enzyme activity. The KEGG analysis for the critical modules indicated that the major mapped pathways included those related to the biosynthesis of other secondary metabolites and amino acid metabolism, such as alanine, aspartate, and glutamate metabolism; arginine and proline metabolism; and glutathione metabolism (Figure 6B). According to these important GO terms and metabolic pathway, some hub proteins were extracted (Table 2). More importantly, most of these pathways directly or indirectly affected the biosynthesis of putrescine and then regulated the synthesis and transport of nicotine in tobacco roots. The results showed that putrescine was the core intermediate product in nicotine anabolism.

FIGURE 6
www.frontiersin.org

Figure 6. Gene Ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of brown module proteins. (A) Gene Ontology function classification of the brown module. (B) Kyoto Encyclopaedia of Genes and Genomes function classification of the brown module.

TABLE 2
www.frontiersin.org

Table 2. Description information of the key proteins of nicotine synthesis.

Enzyme Activity Validation

To verify the quality of proteomics and differential expression level data, four proteins were selected for enzyme activity determination. In the nicotine synthesis pathway, PMT catalyses the conversion of putrescine into N-methyl putrescine, MPO catalyses the conversion of N-methyl putrescine into 4-methylaminobutyl ether, and A622 and BBLs combine nicotinic acid with N-methyl dihydropyrrole salt to synthesise nicotine. The changes in the expression of the selected proteins according to enzyme activity showed a similar expression tendency to the proteomics data (Figure 7), indicating that the proteomics profiling data were reliable.

FIGURE 7
www.frontiersin.org

Figure 7. Validation of proteomics data through enzyme activity determination. A0A1S3ZFC4, A0A1S4BDC4, F1T160, and B6VRE6 are the protein accessions numbers of genes putrescine N-methyltransferase (PMT), N-methylputrescine oxidase (MPO), berberine bridge enzyme-like (BBL), and isoflavone reductase-like protein (A622), respectively.

Discussion

Nicotine has a multi-channel utilisation value, and obtaining a high concentration of nicotine can promote the development of agricultural, pharmaceutical, and chemical industries (Yin, 2018). However, to avoid addiction to cigarette products, the nicotine content in cigarettes has been recommended to be limited to a specific range (2.5–3.5%) (Chen et al., 2004). Therefore, regulating the nicotine content in tobacco can be useful in promoting the application of nicotine in other fields. In recent years, the molecular genetics research on nicotine metabolism has been deepened, and some important genes related to nicotine synthesis, transport and transformation have been cloned successively (Hildreth et al., 2011; Bian et al., 2022), which has played an important role in promoting the research on the mechanism of nicotine synthesis and metabolism and tobacco genetics and breeding. Nicotinoids is a long-distance transport metabolite, and the whole process is also involved in synthases and transporters (Dewey and Xie, 2013). Moreover, the related intermediate metabolites are jointly controlled by adjacent metabolic pathways. Therefore, research on the mechanism of nicotine synthesis should not be limited to the nicotine synthesis pathway.

Proteomics is considered the main pillar of functional genomics, which holds great importance in the analysis of the mechanisms of trait differences between materials (Hochholdinger et al., 2018). Currently available mass spectrometers can scan in the form of DIA to achieve efficient and ready screening of biomarkers (Tsou et al., 2015; Demichev et al., 2020). Our early-stage experiment showed that DIA LC-MS is suitable for tobacco root tips for the detection and analysis of proteomics samples (Mo et al., 2021a). Therefore, in this study, we used the DIA technology to perform a comprehensive proteomics analysis of 12 samples of four tobacco varieties, and a total of 31,259 distinct peptides and 6,018 proteins were identified, which provided a material premise to further understand differences in the nicotine content and synthesis mechanisms. Further analysis showed that the differentially expressed proteins in different materials were mainly involved in the transport and metabolism of substances, and could affect the accumulation of intermediate metabolites of nicotine synthesis from different metabolic pathways.

After nicotinic synthesis, it is first necessary to cross the plasma membrane in the root to be transported to extracellular and thus transferred to the leaves for accumulation (Hashimoto and Yamada, 2003). It was found that the vesicle transport regulator, NtGEF, may indirectly affect the accumulation of nicotine content in leaves by mediating the localisation of nicotine transporter in vacuoles (Yin et al., 2018), and the results showed that vesicle transport was also involved in nicotine transport. Vesicle transport plays an indispensable role in the process of substances transport. Firstly, vesicles must fuse with target cells, and this process is promoted by a protein called soluble N-ethylmaleimide sensitive factor adaptor protein receptors (SNAREs) (Gao et al., 2017). Associated Sec1/Munc18 (SM) proteins in Arabidopsis were found to involve intimal transport between the ER and the Golgi apparatus (Blatt and Thiel, 2018). Vehicle associated membrane proteins 721 and 722 were involved in the transport of substances necessary for plant growth and defence response (Yi et al., 2013; Yun et al., 2013). In this study, we found that SEC1 family transport protein SLY1-like isoform X1 (A0A1S4BG77) and vesicle-associated membrane protein 711-like (A0A1S4AX25) were highly expressed in the high nicotine material Qinggeng. And Go annotation showed that they all bear the function of vacuole mediated transporter in the biological process. Therefore, we speculated that the increased nicotinic content of the Qinggeng material may be driven by these transmembrane transporters, so that the nicotine synthesised in roots can be efficiently transferred to the leaves for accumulation.

Nicotine is a nitrogen-containing compound (Yamada and Sato, 2021), and the intake and metabolism of nitrogen in plants can affect the synthesis and accumulation of nicotine (Shang et al., 2017). Previous studies have shown that over-expression of aspartate aminotransferase (AAT) genes in rice resulted in altered nitrogen metabolism (Zhou et al., 2009). In this study, we found that AAT was significantly upregulated, so we speculated that the upregulation of AAT in tobacco also improved nitrogen metabolism in plants and thus was able to synthesise more nicotinoids. The synthesis of nicotinoids involves many important metabolites, among which arginine and ornithine are the two most important precursor substances (Rubén Alcázar, 1999). Our results showed that these proteins are concentrated in tropane, piperidine and pyridine alkaloid biosynthesis, ornithine metabolism, arginine metabolism, alanine, aspartate and glutamate metabolism, glutathione metabolism, and arginine and proline metabolism. And these pathways can directly or indirectly affect the synthesis of putrescine, which provides a rich material source for the synthesis of nicotine. Putrescine N-methyltransferase (PMT) transferase is the largest rate limiting enzyme, which catalyses putrescine to form N-methylputrescine (Biastoff et al., 2009). In this study, we found that the expression of probable methyltransferase PMT26 (A0A1S3ZFC4) in Qinggeng was significantly higher than that in the other three low nicotine materials, indicating that the increase of catalytic enzyme content also played an important role in the synthesis of nicotine.

A majority of the DEPs identified in plant proteomics are usually enzymes (Lee et al., 2008). They play a vital role in many physiological, biochemical, and signalling pathways and can connect proteome with metabolome (Niu et al., 2018). Proteomics data can be explained and verified by enzyme activity determination to elucidate the mechanisms of protein (enzyme) function and enzyme-related biochemical or metabolic pathways (Boersema et al., 2015; Dhakarey et al., 2017). Xia et al. (2018) used the enzymatic reaction in combination with proteomics to elucidate the seed metabolic mechanism affected by temperature during seed dormancy and germination. In our study, the proteomics data were complemented with the enzyme activity data of four enzymes associated with nicotine synthesis, and the results indicated changes in the expression of the selected proteins. Enzyme activity showed a similar expression tendency to the proteomics data in all four varieties, which indicated that the proteomics profiling data were reliable.

In conclusion, our research explored the hub regulatory proteins of nicotine synthesis from the perspective of proteomics, and revealed that various metabolic pathways affect nicotine synthesis by regulating the formation of putrescine. The results are of great significance for regulating the nicotine content in tobacco and breeding new tobacco varieties with appropriate nicotine content. At the same time, nicotine is a secondary metabolite transported over long distances. Its synthesis and accumulation are not only controlled by its own genotype, but also affected by the growth state of roots, stems, and leaves. Therefore, in the follow-up research, we should also pay more attention to the synergy of source, sink and flow in plants to broaden the network of nicotine metabolism at the overall level.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: The data access link in ProteomeXchange is: http://proteomecentral.proteomexchange.org/cgi/GetDataset?ID=PXD028624; The data access link in iProX is: https://www.iprox.cn/page/project.html?id=IPX00035 10000.

Author Contributions

ZM and RL conceived and designed the experiments. ZM, LD, ZT, KP, WL, and YK performed the experiments. ZM, YP, QN, YH, and RL analysed the data and contributed to the preparation of reagents and materials, use of analytical tools, and manuscript writing. All authors contributed to the article and approved the submitted version.

Funding

This research was supported by the National Natural Science Foundation of China (No. 32060510), the Talent Cultivation Project for High-Level Renovation Plan of Guizhou Province [No. (2016)5663], and Key Projects from the Guizhou Provincial Science and Technology Department [No. (2019)1405].

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.860455/full#supplementary-material

Footnotes

  1. ^ https://www.blast2go.com/
  2. ^ https://www.omicshare.com/tools
  3. ^ http://www.kegg.jp/

References

Baldwin, I. T. (1989). Mechanism of damage-induced alkaloid production in wild tobacco. J. Chem. Ecol. 15, 1661–1680. doi: 10.1007/BF01012392

PubMed Abstract | CrossRef Full Text | Google Scholar

Baldwin, I. T. (1998). Jasmonate-induced responses are costly but benefit plants under attack in native populations. Proc. Natl. Acad. Sci. U. S. A. 95, 8113–8118. doi: 10.1073/pnas.95.14.8113

PubMed Abstract | CrossRef Full Text | Google Scholar

Baldwin, I. T. (1999). Inducible nicotine production in native Nicotiana as an example of adaptive phenotypic plasticity. J. Chem. Ecol. 25, 3–30. doi: 10.1023/A:1020880931488

CrossRef Full Text | Google Scholar

Bian, S., Sui, X., Wang, J., Tian, T., Wang, C., Zhao, X., et al. (2022). NtMYB305a binds to the jasmonate-responsive GAG region of NtPMT1a promoter to regulate nicotine biosynthesis. Plant Physiol. 188, 151–166. doi: 10.1093/plphys/kiab458

PubMed Abstract | CrossRef Full Text | Google Scholar

Biastoff, S., Brandt, W., and Dräger, B. (2009). Putrescine N-methyltransferase - The start for alkaloids. Phytochemistry 70, 1708–1718. doi: 10.1016/j.phytochem.2009.06.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Blatt, M. R., and Thiel, G. (2018). SNARE Components and Mechanisms of Exocytosis in Plants. Annu. Plant Rev. 9, 230–262. doi: 10.1002/9781119312994.apr0089

CrossRef Full Text | Google Scholar

Boersema, P. J., Kahraman, A., and Picotti, P. (2015). Proteomics beyond large-scale protein expression analysis. Curr. Opin. Biotechnol. 34, 162–170. doi: 10.1016/j.copbio.2015.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Cane, K. A., Mayer, M., Lidgett, A. J., Michael, A. J., and Hamill, J. D. (2005). Molecular analysis of alkaloid metabolism in AABB v. aabb genotype Nicotiana tabacum in response to wounding of aerial tissues and methyl jasmonate treatment of cultured roots. Funct. Plant Biol. 32, 305–320. doi: 10.1071/FP04008

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, J., Liu, J., and Long, H. (2004). The distribution characteristics of nutrition elements and main chemical composition in China’s tobacco leaves. Acta Tabacaria Sin. 10, 1–7.

Google Scholar

Chou, W. M., and Kutchan, T. M. (1998). Enzymatic oxidations in the biosynthesis of complex alkaloids. Plant J. 15, 289–300. doi: 10.1046/j.1365-313X.1998.00220.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Cloëz-Tayarani, I., and Changeux, J. (2007). Nicotine and serotonin in immune regulation and inflammatory processes: a perspective. J. Leukoc. Biol. 81, 599–606. doi: 10.1189/jlb.0906544

PubMed Abstract | CrossRef Full Text | Google Scholar

Conesa, A., Götz, S., García-Gómez, J. M., Terol, J., Talón, M., and Robles, M. (2005). Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21, 3674–3676. doi: 10.1093/bioinformatics/bti610

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, H. (2018). Studies on the Function of the PiGK5 Gene and the Phosphorylated Proteomics Analysis of Sexual Reproduction Induced by 1 Sex Hormone in Phytophthora Species. Hohhot: Inner Mongolia Agricultural University.

Google Scholar

Demichev, V., Messner, C. B., Vernardis, S. I., Lilley, K. S., and Ralser, M. (2020). DIA-NN: neural networks and interference correction enable deep proteome coverage in high throughput. Nat. Meth. 17, 41–44. doi: 10.1038/s41592-019-0638-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Dewey, R. E., and Xie, J. (2013). Molecular genetics of alkaloid biosynthesis in Nicotiana tabacum. Phytochemistry 94, 10–27. doi: 10.1016/j.phytochem.2013.06.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Dhakarey, R., Raorane, M. L., Treumann, A., Peethambaran, P. K., Schendel, R. R., Sahi, V. P., et al. (2017). Physiological and proteomic analysis of the rice mutant cpm2 suggests a negative regulatory role of jasmonic acid in drought tolerance. Front. Plant Sci. 8:1903. doi: 10.3389/fpls.2017.01903

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, H., Li, Y., Fan, H., Zhou, D., and Li, H. (2019). Quantitative proteomics analysis reveals resistance differences of banana cultivar ‘Brazilian’ to Fusarium oxysporum f. sp. cubense races 1 and 4. J. Proteomics 203:103376. doi: 10.1016/j.jprot.2019.05.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Duan, H. R., Wang, L. R., Cui, G. X., Zhou, X. H., Duan, X. R., and Yang, H. S. (2020). Identification of the regulatory networks and hub genes controlling alfalfa floral pigmentation variation using RNA-sequencing analysis. BMC Plant Biol. 20:110. doi: 10.1186/s12870-020-2322-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, L., An, Q., Zhang, K., Liu, Y., Tong, Y., Xu, J., et al. (2019). Quantitative proteomic characterization of human sperm cryopreservation: using data-independent acquisition mass spectrometry. BMC Urol. 19:133. doi: 10.1186/s12894-019-0565-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, Y. Q., Chen, J. G., Chen, Z. R., An, D., Lv, Q. Y., Han, M. L., et al. (2017). A new vesicle trafficking regulator CTL1 plays a crucial role in ion homeostasis. PLoS Biol. 15:e2002978. doi: 10.1371/journal.pbio.2002978

PubMed Abstract | CrossRef Full Text | Google Scholar

Hao, S., Hu, J., and Chen, Y. (2017). Effects and mechanisms of nicotine on dopamine release in Parkinson’s disease mouse model Bio - Rad Gel Doc / Chemi Doc. J. Nanjing Med. Univ. Sci. 37, 686–690. doi: 10.7655/NYDXBNS20170607

CrossRef Full Text | Google Scholar

Hashimoto, T., and Yamada, Y. (2003). New genes in alkaloid metabolism and transport. Curr. Opin. Biotechnol. 14, 163–168. doi: 10.1016/S0958-1669(03)00027-2

CrossRef Full Text | Google Scholar

Hibi, N., Higashiguchi, S., Hashimoto, T., Yamada, Y., The, S., Cell, P., et al. (1994). Gene Expression in Tobacco Low-Nicotine Mutants Linked references are available on JSTOR for this article: gene Expression in Tobacco Low-Nicotine Mutants. Plant Cell 6, 723–735. doi: 10.2307/3869875

CrossRef Full Text | Google Scholar

Hildreth, S. B., Gehman, E. A., Yang, H., Lu, R. H., Ritesh, K. C., Harich, K. C., et al. (2011). Tobacco nicotine uptake permease (NUP1) affects alkaloid metabolism. Proc. Natl. Acad. Sci. U. S. A. 108, 18179–18184. doi: 10.1073/pnas.1108620108

PubMed Abstract | CrossRef Full Text | Google Scholar

Hochholdinger, F., Marcon, C., Baldauf, J. A., Yu, P., and Frey, F. P. (2018). Proteomics of maize root development. Front. Plant Sci. 9:143. doi: 10.3389/fpls.2018.00143

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, X., Wang, H., Li, K., Wu, Y., Liu, Z., and Huang, C. (2017). Genome-wide proteomic profiling reveals the role of dominance protein expression in heterosis in immature maize ears. Sci. Rep. 7:16130. doi: 10.1038/s41598-017-15985-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Imanishi, S., Hashizume, K., Nakakita, M., Kojima, H., Matsubayashi, Y., Hashimoto, T., et al. (1998). Differential induction by methyl jasmonate of genes encoding ornithine decarboxylase and other enzymes involved in nicotine biosynthesis in tobacco cell cultures. Plant Mol. Biol. 38, 1101–1111. doi: 10.1023/A:1006058700949

CrossRef Full Text | Google Scholar

Jung, J. W., Park, S. Y., Oh, S. D., Jang, Y., Suh, S. J., Park, S. K., et al. (2021). Metabolomic Variability of Different Soybean Genotypes: β-Carotene-Enhanced (Glycine max), Wild (Glycine soja), and Hybrid (Glycine max × Glycine soja) Soybeans. Foods 10:2421. doi: 10.3390/foods10102421

PubMed Abstract | CrossRef Full Text | Google Scholar

Jung, J. Y., Shin, R., and Schachtman, D. P. (2009). Ethylene mediates response and tolerance to potassium deprivation in arabidopsis. Plant Cell 21, 607–621. doi: 10.1105/tpc.108.063099

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanehisa, M., Goto, S., Sato, Y., Furumichi, M., and Tanabe, M. (2012). KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 40, 109–114. doi: 10.1093/nar/gkr988

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, R., Bohra, A., Pandey, A. K., Pandey, M. K., and Kumar, A. (2017). Metabolomics for plant improvement: status and prospects. Front. Plant Sci. 8:1302. doi: 10.3389/fpls.2017.01302

PubMed Abstract | CrossRef Full Text | Google Scholar

Langfelder, P., and Horvath, S. (2008). WGCNA: an R package for weighted correlation network analysis. BMC Bioinform. 9:559. doi: 10.1186/1471-2105-9-559

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, C. P., Eubel, H., O’Toole, N., and Millar, A. H. (2008). Heterogeneity of the mitochondrial proteome for photosynthetic and non-photosynthetic Arabidopsis metabolism. Mol. Cell. Proteom. 7, 1297–1316. doi: 10.1074/mcp.M700535-MCP200

PubMed Abstract | CrossRef Full Text | Google Scholar

Legg, P. D., Chaplin, J. F., and Collins, G. B. (1969). Inheritance of Percent Total Alkaloids in Nicotiana tabacum L. J. Hered. 60, 213–217. doi: 10.1093/oxfordjournals.jhered.a107974

CrossRef Full Text | Google Scholar

Lewis, R. S. (2018). Potential Mandated Lowering of Nicotine Levels in Cigarettes: A Plant Perspective. Nicotine Tob. Res. 21, 991–995. doi: 10.1093/ntr/nty022/4833192

PubMed Abstract | CrossRef Full Text | Google Scholar

Lewis, R. S., Drake-Stowe, K. E., Heim, C., Steede, T., Smith, W., and Dewey, R. E. (2020). Genetic and Agronomic Analysis of Tobacco Genotypes Exhibiting Reduced Nicotine Accumulation Due to Induced Mutations in Berberine Bridge Like (BBL) Genes. Front. Plant Sci. 11:368. doi: 10.3389/fpls.2020.00368

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C., Liu, J. C., Xiao, X., Chen, X., Yue, S., Yu, H., et al. (2017). Psychological distress and type 2 diabetes mellitus: A 4-year policemen cohort study in China. BMJ Open 7:e014235. doi: 10.1136/bmjopen-2016-014235

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, S., Yu, Q., Han, L., Zhang, Y., Tian, X., and Zhao, S. (2019). Effects of proteome changes on the tenderness of yak rumen smooth muscle during postmortem storage based on the label-free mass spectrometry. Food Res. Int. 116, 1336–1343. doi: 10.1016/j.foodres.2018.10.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, P., Zheng, W., Qiao, J., Li, Z., Deng, Z., Yuan, Y., et al. (2021). Early embryonic transcriptomes of Zeugodacus tau provide insight into sex determination and differentiation genes. Insect Sci. 10, 2962–2968. doi: 10.1111/1744-7917.12974

PubMed Abstract | CrossRef Full Text | Google Scholar

Lou, Q., Chen, L., Mei, H., Xu, K., Wei, H., Feng, F., et al. (2017). Root transcriptomic analysis revealing the importance of energy metabolism to the development of deep roots in rice (Oryza sativa L.). Front. Plant Sci. 8:1314. doi: 10.3389/fpls.2017.01314

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, C., Pu, Y., Liu, Y., Li, Y., Qu, J., Huang, H., et al. (2019). Comparative transcriptomics and weighted gene co-expression correlation network analysis (WGCNA) reveal potential regulation mechanism of carotenoid accumulation in Chrysanthemum × morifolium. Plant Physiol. Biochem. 142, 415–428. doi: 10.1016/j.plaphy.2019.07.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Marton, L. J., and Pegg, A. E. (1995). Polyamines as targets for therapeutic intervention. Annu. Rev. Pharmacol. Toxicol. 35, 55–91. doi: 10.1146/annurev.pa.35.040195.000415

PubMed Abstract | CrossRef Full Text | Google Scholar

Mo, Z., Pu, Y., Zhou, J., Tian, Z., Teng, J., Chen, Q., et al. (2021b). Effect of the over-dominant expression of proteins on nicotine heterosis via proteomic analysis. Sci. Rep. 11:21063. doi: 10.1038/s41598-021-00614-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Mo, Z., Chen, X., Pu, Y., Duan, L., Tian, Z., Quan, Y., et al. (2021a). Detection and analysis of protein samples of tobacco root tip. J. South. Agric. 52, 2962–2968. doi: 10.3969/j.issn.2095-1191.2021.11.006

CrossRef Full Text | Google Scholar

Moriya, Y., Itoh, M., Okuda, S., Yoshizawa, A. C., and Kanehisa, M. (2007). KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 35, 182–185. doi: 10.1093/nar/gkm321

PubMed Abstract | CrossRef Full Text | Google Scholar

Mukhopadhyay, I., Chowdhuri, D. K., Bajpayee, M., and Dhawan, A. (2004). Evaluation of in vivo genotoxicity of cypermethrin in Drosophila melanogaster using the alkaline Comet assay. Mutagenesis 19, 85–90. doi: 10.1093/mutage/geh007

PubMed Abstract | CrossRef Full Text | Google Scholar

Niu, L., Xu, X., Liu, H., Wu, X., and Wang, W. (2018). On the Promising Role of Enzyme Activity Assay in Interpreting Comparative Proteomic Data in Plants. Proteomics 18:e1800234. doi: 10.1002/pmic.201800234

PubMed Abstract | CrossRef Full Text | Google Scholar

Nölke, G., Volke, D., Chudobová, I., Houdelet, M., Lusso, M., Frederick, J., et al. (2018). Polyamines delay leaf maturation in low-alkaloid tobacco varieties. Plant Direct 2:e00077. doi: 10.1002/pld3.77

PubMed Abstract | CrossRef Full Text | Google Scholar

Pu, Y. (2020). The effect of key enzyme activity on the formation of heterosis in nicotine synthesis pathway. Guizhou University 4, 1–60. doi: 10.27047/d.cnki.ggudu.2020.001108

CrossRef Full Text | Google Scholar

Rubén Alcázar, A. F. T. (1999). Determination of Arginine and Ornithine Decarboxylase Activities in Plants. Orthod. Treat. Cl. III Malocclusion 36, 117–122. doi: 10.1007/978-1-4939-7398-9

CrossRef Full Text | Google Scholar

Sanchez, S., and Demain, A. L. (2019). Secondary metabolites. Compr. Biotechnol. 10, 131–143. doi: 10.1016/B978-0-444-64046-8.00012-4

CrossRef Full Text | Google Scholar

Sato, F., Hashimoto, T., Hachiya, A., Tamura, K., Choi, K., Morishige, T., et al. (2011). Metabolic engineering of plant alkaloid biosynthesis. Proc. Natl. Acad. Sci. U.S.A. 22, 926–927. doi: 10.1111/nep.12989

PubMed Abstract | CrossRef Full Text | Google Scholar

Shang, C., Chen, A., Chen, G., Li, H., Guan, S., and He, J. (2017). Microbial Biofertilizer Decreases Nicotine Content by Improving Soil Nitrogen Supply. Appl. Biochem. Biotechnol. 181, 1–14. doi: 10.1007/s12010-016-2195-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, V., Gupta, P., Priscilla, K., Sharankumar Hangargi, B., Veershetty, A., et al. (2021). Metabolomics intervention towards better understanding of plant traits. Cells 10:346. doi: 10.3390/cells10020346

PubMed Abstract | CrossRef Full Text | Google Scholar

Shoji, T., and Hashimoto, T. (2013). Smoking out the masters: transcriptional regulators for nicotine biosynthesis in tobacco. Plant Biotechnol. 30, 217–224. doi: 10.5511/plantbiotechnology.13.0221a

CrossRef Full Text | Google Scholar

Shoji, T., Inai, K., Yazaki, Y., Sato, Y., Takase, H., Shitan, N., et al. (2009). Multidrug and toxic compound extrusion-type transporters implicated in vacuolar sequestration of nicotine in tobacco roots. Plant Physiol. 149, 708–718. doi: 10.1104/pp.108.132811

PubMed Abstract | CrossRef Full Text | Google Scholar

Steppuhn, A., Gase, K., Krock, B., Halitschke, R., and Baldwin, I. T. (2004). Nicotine’s defensive function in nature. PLoS Biol. 2:1074–1080. doi: 10.1371/journal.pbio.0020217

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, D. S., Zhang, L., Chen, H. L., Liang, Y. R., Lu, J. L., Liang, H. L., et al. (2007). Extraction and purification of solanesol from tobacco. (I). Extraction and silica gel column chromatography separation of solanesol. Sep. Purif. Technol. 56, 291–295. doi: 10.1016/j.seppur.2007.01.040

CrossRef Full Text | Google Scholar

Thurston, R., Smith, W. T., and Cooper, B. P. (1966). Alkaloid secretion by trichomes of Nicotiana species and resistance to aphids. Entomol. Exp. Appl. 9, 428–432. doi: 10.1111/j.1570-7458.1966.tb01003.x

CrossRef Full Text | Google Scholar

Tian, M., Nie, Q., Li, Z., Zhang, J., Liu, Y., Long, Y., et al. (2018). Transcriptomic analysis reveals overdominance playing a critical role in nicotine heterosis in Nicotiana tabacum L. BMC Plant Biol. 18:48. doi: 10.1186/s12870-018-1257-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Tiburcio, A. F., and Galston, A. W. (2016). Alkaloid Biosynthesis: the Basis for Metabolic Engineering of Medicinal Plants. Plant Cell 7, 1059–1070. doi: 10.2307/3870057

CrossRef Full Text | Google Scholar

Tsou, C. C., Avtonomov, D., Larsen, B., Tucholska, M., Choi, H., Gingras, A. C., et al. (2015). DIA-Umpire: comprehensive computational framework for data-independent acquisition proteomics. Nat. Meth. 12, 258–264. doi: 10.1038/nmeth.3255

PubMed Abstract | CrossRef Full Text | Google Scholar

Wagner, R., Feth, F., and Wagner, K. G. (1986). Regulation in tobacco callus of enzyme activities of the nicotine pathway: II. The pyridine-nucleotide cycle. Enzyme 168, 408–413. doi: 10.1007/BF00392369

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Yang, S., An, B., Wang, S., Yin, Y., Lu, Y., et al. (2011). Molecular dynamics analysis reveals structural insights into mechanism of nicotine N-demethylation catalyzed by tobacco cytochrome P450 mono-oxygenase. PLoS One 6:e23342. doi: 10.1371/journal.pone.0023342

PubMed Abstract | CrossRef Full Text | Google Scholar

World Health Organization (2009). WHO Study Group on Tobacco Product Regulation. Available online at: http://pesquisa.bvsalud.org/bvsms/resource/pt/mis-33477 (accessed May 4, 2021).

Google Scholar

Xia, Q., Ponnaiah, M., Cueff, G., Rajjou, L., Prodhomme, D., Gibon, Y., et al. (2018). Corrigendum to “Integrating proteomics and enzymatic profiling to decipher seed metabolism affected by temperature in seed dormancy and germination”. Plant Sci. 269, 118–125. doi: 10.1016/j.plantsci.2018.05.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, H. (2019). Study on the mechanism of protein-related changes of wall polysaccharides in different growth stages of Lycium barbarum l. Ningxia University 2, 1–60. doi: 10.27257/d.cnki.gnxhc.2019.000924

CrossRef Full Text | Google Scholar

Yamada, Y., and Sato, F. (2021). Transcription factors in alkaloid engineering. Biomolecules 11:1719. doi: 10.3390/biom11111719

PubMed Abstract | CrossRef Full Text | Google Scholar

Yi, C., Park, S., Yun, H. S., and Kwon, C. (2013). Vesicle-associated membrane proteins 721 and 722 are required for unimpeded growth of Arabidopsis under ABA application. J. Plant Physiol. 170, 529–533. doi: 10.1016/j.jplph.2012.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Yin, P. (2018). The gene expression of NUP1, GEF and relationship with transport efficiency of nicotine in tobacco. Chinese Academy of Agricultural Sciences 12, 1–58.

Google Scholar

Yin, P., Ren, M., Chen, A., Sun, X., Zhou, S., and Zhang, Y. (2018). Study on the Correlation between Nicotine Accumulation and Nicotine. CHIN. Tob. Sci. 39, 25–32. doi: 10.13496/j.issn.1007-5119.2018.05.004

CrossRef Full Text | Google Scholar

Yuan, Y., Qin, L., Su, H., Yang, S., Wei, X., Wang, Z., et al. (2021). Transcriptome and Coexpression Network Analyses Reveal Hub Genes in Chinese Cabbage (Brassica rapa L. ssp. pekinensis) During Different Stages of Plasmodiophora brassicae Infection. Front. Plant Sci. 12:650252. doi: 10.3389/fpls.2021.650252

PubMed Abstract | CrossRef Full Text | Google Scholar

Yun, H. S., Kwaaitaal, M., Kato, N., Yi, C., Park, S., Sato, M. H., et al. (2013). Requirement of vesicle-associated membrane protein 721 and 722 for sustained growth during immune responses in Arabidopsis. Mol. Cells 35, 481–488. doi: 10.1007/s10059-013-2130-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H., Wang, F., Luo, J., Chen, Y., Ye, W., and Chen, J. (2019). Effects of soil and climate on photosynthetic characteristics and chemical constituents of strong-flavor flue-cured tobacco in different ecological subregions. J. South CHIN. Agric. Univ. 40, 23–31. doi: 10.7671/j.issn.1001-411X.201804023

CrossRef Full Text | Google Scholar

Zhang, H. B., Bokowiec, M. T., Rushton, P. J., Han, S. C., and Timko, M. P. (2012). Tobacco transcription factors NtMYC2a and NtMYC2b form nuclear complexes with the NtJAZ1 repressor and regulate multiple jasmonate-inducible steps in nicotine biosynthesis. Mol. Plant 5, 73–84. doi: 10.1093/mp/ssr056

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Y., Cai, H., Xiao, J., Li, X., Zhang, Q., and Lian, X. (2009). Over-expression of aspartate aminotransferase genes in rice resulted in altered nitrogen metabolism and increased amino acid content in seeds. Theor. Appl. Genet. 118, 1381–1390. doi: 10.1007/s00122-009-0988-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: nicotine, proteomics, weighted gene co-expression network analysis, metabolites, enzyme activity, Nicotiana tabacum

Citation: Mo Z, Duan L, Pu Y, Tian Z, Ke Y, Luo W, Pi K, Huang Y, Nie Q and Liu R (2022) Proteomics and Co-expression Network Analysis Reveal the Importance of Hub Proteins and Metabolic Pathways in Nicotine Synthesis and Accumulation in Tobacco (Nicotiana tabacum L.). Front. Plant Sci. 13:860455. doi: 10.3389/fpls.2022.860455

Received: 23 January 2022; Accepted: 24 March 2022;
Published: 28 April 2022.

Edited by:

Mohamed Suhail Rafudeen, University of Cape Town, South Africa

Reviewed by:

Atsushi Fukushima, Kyoto Prefectural University, Japan
Haiying Li, Heilongjiang University, China

Copyright © 2022 Mo, Duan, Pu, Tian, Ke, Luo, Pi, Huang, Nie and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Renxiang Liu, rxliu@gzu.edu.cn

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.