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.)

Abstract

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

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 Blast2Go¹ (Conesa et al., 2005). The Cluster of Orthologous Groups of proteins (COG) analysis was performed using the OmicShare tools.² 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)³ 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

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

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

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

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

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

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.

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.

References

• 1

  BaldwinI. T. (1989). Mechanism of damage-induced alkaloid production in wild tobacco.J. Chem. Ecol.151661–1680. 10.1007/BF01012392

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  BaldwinI. T. (1998). Jasmonate-induced responses are costly but benefit plants under attack in native populations.Proc. Natl. Acad. Sci. U. S. A.958113–8118. 10.1073/pnas.95.14.8113

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  BaldwinI. T. (1999). Inducible nicotine production in native Nicotiana as an example of adaptive phenotypic plasticity.J. Chem. Ecol.253–30. 10.1023/A:1020880931488

  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  BianS.SuiX.WangJ.TianT.WangC.ZhaoX.et al (2022). NtMYB305a binds to the jasmonate-responsive GAG region of NtPMT1a promoter to regulate nicotine biosynthesis.Plant Physiol.188151–166. 10.1093/plphys/kiab458

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  BiastoffS.BrandtW.DrägerB. (2009). Putrescine N-methyltransferase - The start for alkaloids.Phytochemistry701708–1718. 10.1016/j.phytochem.2009.06.012

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  BlattM. R.ThielG. (2018). SNARE Components and Mechanisms of Exocytosis in Plants.Annu. Plant Rev.9230–262. 10.1002/9781119312994.apr0089

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  BoersemaP. J.KahramanA.PicottiP. (2015). Proteomics beyond large-scale protein expression analysis.Curr. Opin. Biotechnol.34162–170. 10.1016/j.copbio.2015.01.005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  CaneK. A.MayerM.LidgettA. J.MichaelA. J.HamillJ. 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.32305–320. 10.1071/FP04008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  ChenJ.LiuJ.LongH. (2004). The distribution characteristics of nutrition elements and main chemical composition in China’s tobacco leaves.Acta Tabacaria Sin.101–7.

  • Google Scholar
  • View reference in article

• 10

  ChouW. M.KutchanT. M. (1998). Enzymatic oxidations in the biosynthesis of complex alkaloids.Plant J.15289–300. 10.1046/j.1365-313X.1998.00220.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  Cloëz-TayaraniI.ChangeuxJ. (2007). Nicotine and serotonin in immune regulation and inflammatory processes: a perspective.J. Leukoc. Biol.81599–606. 10.1189/jlb.0906544

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  ConesaA.GötzS.García-GómezJ. M.TerolJ.TalónM.RoblesM. (2005). Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research.Bioinformatics213674–3676. 10.1093/bioinformatics/bti610

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  CuiH. (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
  • View reference in article

• 14

  DemichevV.MessnerC. B.VernardisS. I.LilleyK. S.RalserM. (2020). DIA-NN: neural networks and interference correction enable deep proteome coverage in high throughput.Nat. Meth.1741–44. 10.1038/s41592-019-0638-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  DeweyR. E.XieJ. (2013). Molecular genetics of alkaloid biosynthesis in Nicotiana tabacum.Phytochemistry9410–27. 10.1016/j.phytochem.2013.06.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  DhakareyR.RaoraneM. L.TreumannA.PeethambaranP. K.SchendelR. R.SahiV. 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. 10.3389/fpls.2017.01903

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  DongH.LiY.FanH.ZhouD.LiH. (2019). Quantitative proteomics analysis reveals resistance differences of banana cultivar ‘Brazilian’ to Fusarium oxysporum f. sp. cubense races 1 and 4.J. Proteomics203:103376. 10.1016/j.jprot.2019.05.004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  DuanH. R.WangL. R.CuiG. X.ZhouX. H.DuanX. R.YangH. S. (2020). Identification of the regulatory networks and hub genes controlling alfalfa floral pigmentation variation using RNA-sequencing analysis.BMC Plant Biol.20:110. 10.1186/s12870-020-2322-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  FuL.AnQ.ZhangK.LiuY.TongY.XuJ.et al (2019). Quantitative proteomic characterization of human sperm cryopreservation: using data-independent acquisition mass spectrometry.BMC Urol.19:133. 10.1186/s12894-019-0565-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  GaoY. Q.ChenJ. G.ChenZ. R.AnD.LvQ. Y.HanM. L.et al (2017). A new vesicle trafficking regulator CTL1 plays a crucial role in ion homeostasis.PLoS Biol.15:e2002978. 10.1371/journal.pbio.2002978

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  HaoS.HuJ.ChenY. (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.37686–690. 10.7655/NYDXBNS20170607

  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  HashimotoT.YamadaY. (2003). New genes in alkaloid metabolism and transport.Curr. Opin. Biotechnol.14163–168. 10.1016/S0958-1669(03)00027-2

  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  HibiN.HigashiguchiS.HashimotoT.YamadaY.TheS.CellP.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 Cell6723–735. 10.2307/3869875

  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  HildrethS. B.GehmanE. A.YangH.LuR. H.RiteshK. C.HarichK. C.et al (2011). Tobacco nicotine uptake permease (NUP1) affects alkaloid metabolism.Proc. Natl. Acad. Sci. U. S. A.10818179–18184. 10.1073/pnas.1108620108

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  HochholdingerF.MarconC.BaldaufJ. A.YuP.FreyF. P. (2018). Proteomics of maize root development.Front. Plant Sci.9:143. 10.3389/fpls.2018.00143

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  HuX.WangH.LiK.WuY.LiuZ.HuangC. (2017). Genome-wide proteomic profiling reveals the role of dominance protein expression in heterosis in immature maize ears.Sci. Rep.7:16130. 10.1038/s41598-017-15985-3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  ImanishiS.HashizumeK.NakakitaM.KojimaH.MatsubayashiY.HashimotoT.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.381101–1111. 10.1023/A:1006058700949

  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  JungJ. W.ParkS. Y.OhS. D.JangY.SuhS. J.ParkS. K.et al (2021). Metabolomic Variability of Different Soybean Genotypes: β-Carotene-Enhanced (Glycine max), Wild (Glycine soja), and Hybrid (Glycine max × Glycine soja) Soybeans.Foods10:2421. 10.3390/foods10102421

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  JungJ. Y.ShinR.SchachtmanD. P. (2009). Ethylene mediates response and tolerance to potassium deprivation in arabidopsis.Plant Cell21607–621. 10.1105/tpc.108.063099

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  KanehisaM.GotoS.SatoY.FurumichiM.TanabeM. (2012). KEGG for integration and interpretation of large-scale molecular data sets.Nucleic Acids Res.40109–114. 10.1093/nar/gkr988

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  KumarR.BohraA.PandeyA. K.PandeyM. K.KumarA. (2017). Metabolomics for plant improvement: status and prospects.Front. Plant Sci.8:1302. 10.3389/fpls.2017.01302

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  LangfelderP.HorvathS. (2008). WGCNA: an R package for weighted correlation network analysis.BMC Bioinform.9:559. 10.1186/1471-2105-9-559

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  LeeC. P.EubelH.O’TooleN.MillarA. H. (2008). Heterogeneity of the mitochondrial proteome for photosynthetic and non-photosynthetic Arabidopsis metabolism.Mol. Cell. Proteom.71297–1316. 10.1074/mcp.M700535-MCP200

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  LeggP. D.ChaplinJ. F.CollinsG. B. (1969). Inheritance of Percent Total Alkaloids in Nicotiana tabacum L.J. Hered.60213–217. 10.1093/oxfordjournals.jhered.a107974

  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  LewisR. S. (2018). Potential Mandated Lowering of Nicotine Levels in Cigarettes: A Plant Perspective.Nicotine Tob. Res.21991–995. 10.1093/ntr/nty022/4833192

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  LewisR. S.Drake-StoweK. E.HeimC.SteedeT.SmithW.DeweyR. 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. 10.3389/fpls.2020.00368

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  LiC.LiuJ. C.XiaoX.ChenX.YueS.YuH.et al (2017). Psychological distress and type 2 diabetes mellitus: A 4-year policemen cohort study in China.BMJ Open7:e014235. 10.1136/bmjopen-2016-014235

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  LiS.YuQ.HanL.ZhangY.TianX.ZhaoS. (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.1161336–1343. 10.1016/j.foodres.2018.10.023

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  LiuP.ZhengW.QiaoJ.LiZ.DengZ.YuanY.et al (2021). Early embryonic transcriptomes of Zeugodacus tau provide insight into sex determination and differentiation genes.Insect Sci.102962–2968. 10.1111/1744-7917.12974

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 40

  LouQ.ChenL.MeiH.XuK.WeiH.FengF.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. 10.3389/fpls.2017.01314

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 41

  LuC.PuY.LiuY.LiY.QuJ.HuangH.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.142415–428. 10.1016/j.plaphy.2019.07.023

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 42

  MartonL. J.PeggA. E. (1995). Polyamines as targets for therapeutic intervention.Annu. Rev. Pharmacol. Toxicol.3555–91. 10.1146/annurev.pa.35.040195.000415

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 43

  MoZ.PuY.ZhouJ.TianZ.TengJ.ChenQ.et al (2021b). Effect of the over-dominant expression of proteins on nicotine heterosis via proteomic analysis.Sci. Rep.11:21063. 10.1038/s41598-021-00614-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 44

  MoZ.ChenX.PuY.DuanL.TianZ.QuanY.et al (2021a). Detection and analysis of protein samples of tobacco root tip.J. South. Agric.522962–2968. 10.3969/j.issn.2095-1191.2021.11.006

  • CrossRef
  • Google Scholar
  • View reference in article

• 45

  MoriyaY.ItohM.OkudaS.YoshizawaA. C.KanehisaM. (2007). KAAS: an automatic genome annotation and pathway reconstruction server.Nucleic Acids Res.35182–185. 10.1093/nar/gkm321

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 46

  MukhopadhyayI.ChowdhuriD. K.BajpayeeM.DhawanA. (2004). Evaluation of in vivo genotoxicity of cypermethrin in Drosophila melanogaster using the alkaline Comet assay.Mutagenesis1985–90. 10.1093/mutage/geh007

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 47

  NiuL.XuX.LiuH.WuX.WangW. (2018). On the Promising Role of Enzyme Activity Assay in Interpreting Comparative Proteomic Data in Plants.Proteomics18:e1800234. 10.1002/pmic.201800234

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 48

  NölkeG.VolkeD.ChudobováI.HoudeletM.LussoM.FrederickJ.et al (2018). Polyamines delay leaf maturation in low-alkaloid tobacco varieties.Plant Direct2:e00077. 10.1002/pld3.77

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 49

  PuY. (2020). The effect of key enzyme activity on the formation of heterosis in nicotine synthesis pathway. Guizhou University4, 1–60. 10.27047/d.cnki.ggudu.2020.001108

  • CrossRef
  • Google Scholar
  • View reference in article

• 50

  Rubén AlcázarA. F. T. (1999). Determination of Arginine and Ornithine Decarboxylase Activities in Plants.Orthod. Treat. Cl. III Malocclusion36117–122. 10.1007/978-1-4939-7398-9

  • CrossRef
  • Google Scholar
  • View reference in article

• 51

  SanchezS.DemainA. L. (2019). Secondary metabolites.Compr. Biotechnol.10131–143. 10.1016/B978-0-444-64046-8.00012-4

  • CrossRef
  • Google Scholar
  • View reference in article

• 52

  SatoF.HashimotoT.HachiyaA.TamuraK.ChoiK.MorishigeT.et al (2011). Metabolic engineering of plant alkaloid biosynthesis.Proc. Natl. Acad. Sci. U.S.A.22926–927. 10.1111/nep.12989

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 53

  ShangC.ChenA.ChenG.LiH.GuanS.HeJ. (2017). Microbial Biofertilizer Decreases Nicotine Content by Improving Soil Nitrogen Supply.Appl. Biochem. Biotechnol.1811–14. 10.1007/s12010-016-2195-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 54

  SharmaV.GuptaP.PriscillaK.SharankumarHangargiB.VeershettyA.et al (2021). Metabolomics intervention towards better understanding of plant traits.Cells10:346. 10.3390/cells10020346

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 55

  ShojiT.HashimotoT. (2013). Smoking out the masters: transcriptional regulators for nicotine biosynthesis in tobacco.Plant Biotechnol.30217–224. 10.5511/plantbiotechnology.13.0221a

  • CrossRef
  • Google Scholar
  • View reference in article

• 56

  ShojiT.InaiK.YazakiY.SatoY.TakaseH.ShitanN.et al (2009). Multidrug and toxic compound extrusion-type transporters implicated in vacuolar sequestration of nicotine in tobacco roots.Plant Physiol.149708–718. 10.1104/pp.108.132811

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 57

  SteppuhnA.GaseK.KrockB.HalitschkeR.BaldwinI. T. (2004). Nicotine’s defensive function in nature.PLoS Biol.2:1074–1080. 10.1371/journal.pbio.0020217

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 58

  TangD. S.ZhangL.ChenH. L.LiangY. R.LuJ. L.LiangH. L.et al (2007). Extraction and purification of solanesol from tobacco. (I). Extraction and silica gel column chromatography separation of solanesol.Sep. Purif. Technol.56291–295. 10.1016/j.seppur.2007.01.040

  • CrossRef
  • Google Scholar
  • View reference in article

• 59

  ThurstonR.SmithW. T.CooperB. P. (1966). Alkaloid secretion by trichomes of Nicotiana species and resistance to aphids.Entomol. Exp. Appl.9428–432. 10.1111/j.1570-7458.1966.tb01003.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 60

  TianM.NieQ.LiZ.ZhangJ.LiuY.LongY.et al (2018). Transcriptomic analysis reveals overdominance playing a critical role in nicotine heterosis in Nicotiana tabacum L.BMC Plant Biol.18:48. 10.1186/s12870-018-1257-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 61

  TiburcioA. F.GalstonA. W. (2016). Alkaloid Biosynthesis: the Basis for Metabolic Engineering of Medicinal Plants.Plant Cell71059–1070. 10.2307/3870057

  • CrossRef
  • Google Scholar
  • View reference in article

• 62

  TsouC. C.AvtonomovD.LarsenB.TucholskaM.ChoiH.GingrasA. C.et al (2015). DIA-Umpire: comprehensive computational framework for data-independent acquisition proteomics.Nat. Meth.12258–264. 10.1038/nmeth.3255

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 63

  WagnerR.FethF.WagnerK. G. (1986). Regulation in tobacco callus of enzyme activities of the nicotine pathway: II. The pyridine-nucleotide cycle.Enzyme168408–413. 10.1007/BF00392369

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 64

  WangS.YangS.AnB.WangS.YinY.LuY.et al (2011). Molecular dynamics analysis reveals structural insights into mechanism of nicotine N-demethylation catalyzed by tobacco cytochrome P450 mono-oxygenase.PLoS One6:e23342. 10.1371/journal.pone.0023342

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 65

  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
  • View reference in article

• 66

  XiaQ.PonnaiahM.CueffG.RajjouL.ProdhommeD.GibonY.et al (2018). Corrigendum to “Integrating proteomics and enzymatic profiling to decipher seed metabolism affected by temperature in seed dormancy and germination”.Plant Sci.269118–125. 10.1016/j.plantsci.2018.05.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 67

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

  • CrossRef
  • Google Scholar
  • View reference in article

• 68

  YamadaY.SatoF. (2021). Transcription factors in alkaloid engineering.Biomolecules11:1719. 10.3390/biom11111719

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 69

  YiC.ParkS.YunH. S.KwonC. (2013). Vesicle-associated membrane proteins 721 and 722 are required for unimpeded growth of Arabidopsis under ABA application.J. Plant Physiol.170529–533. 10.1016/j.jplph.2012.11.001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 70

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

  • Google Scholar

• 71

  YinP.RenM.ChenA.SunX.ZhouS.ZhangY. (2018). Study on the Correlation between Nicotine Accumulation and Nicotine.CHIN. Tob. Sci.3925–32. 10.13496/j.issn.1007-5119.2018.05.004

  • CrossRef
  • Google Scholar

• 72

  YuanY.QinL.SuH.YangS.WeiX.WangZ.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. 10.3389/fpls.2021.650252

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 73

  YunH. S.KwaaitaalM.KatoN.YiC.ParkS.SatoM. H.et al (2013). Requirement of vesicle-associated membrane protein 721 and 722 for sustained growth during immune responses in Arabidopsis.Mol. Cells35481–488. 10.1007/s10059-013-2130-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 74

  ZhangH.WangF.LuoJ.ChenY.YeW.ChenJ. (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.4023–31. 10.7671/j.issn.1001-411X.201804023

  • CrossRef
  • Google Scholar

• 75

  ZhangH. B.BokowiecM. T.RushtonP. J.HanS. C.TimkoM. 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. Plant573–84. 10.1093/mp/ssr056

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 76

  ZhouY.CaiH.XiaoJ.LiX.ZhangQ.LianX. (2009). Over-expression of aspartate aminotransferase genes in rice resulted in altered nitrogen metabolism and increased amino acid content in seeds.Theor. Appl. Genet.1181381–1390. 10.1007/s00122-009-0988-3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar