NtMYB27 acts downstream of NtBES1 to modulate flavonoids accumulation in response to UV-B radiation in tobacco.

UV-B radiation can induce the accumulation of many secondary metabolites, including flavonoids, in plants to protect them from oxidative damage. BRI1-EMS-SUPPRESSOR1 (BES1) has been shown to mediate the biosynthesis of flavonoids in response to UV-B. However, the detailed mechanism by which it acts still needs to be further elucidated. Here, we revealed that UV-B significantly inhibited the transcription of multiple transcription factor genes in tobacco, including NtMYB27, which was subsequently shown to be a repressor of flavonoids synthesis in tobacco. We further demonstrated that NtBES1 directly binds to the E-box motifs present in the promoter of NtMYB27 to mediate its transcriptional repression upon UV-B exposure. The UV-B-repressed NtMYB27 could bind to the ACCT-containing element (ACE) in the promoters of Nt4CL and NtCHS and served as a modulator that promoted the biosynthesis of lignin and chlorogenic acid (CGA) but inhibited the accumulation of flavonoids in tobacco. The expression of NtMYB27 was also significantly repressed by heat stress, suggesting its putative roles in regulating heat-induced flavonoids accumulation. Taken together, our results revealed the role of NtBES1 and NtMYB27 in regulating the synthesis of flavonoids during the plant response to UV-B radiation in tobacco.

Transcriptomics and virus-induced gene silencing identify defence-related genes during Ralstonia solanacearum infection in resistant and susceptible tobacco.

Bacterial wilt (BW) caused by Ralstonia solanacearum is a globally prevalent bacterial soil-borne disease. In this study, transcriptome sequencing were subjected to roots after infection with the R. solanacearum in the resistant and susceptible tobacco variety. DEGs that responded to R. solanacearum infection in both resistant and susceptible tobacco contributed to pectinase and peroxidase development and were enriched in plant hormone signal transduction, signal transduction and MAPK signalling pathway KEGG terms. Core DEGs in the resistant tobacco response to R. solanacearum infection were enriched in cell wall, membrane, abscisic acid and ethylene terms. qRT-PCR indicated that Nitab4.5_0004899g0110, Nitab4.5_0004234g0080 and Nitab4.5_0001439g0050 contributed to the response to R. solanacearum infection in different resistant and susceptible tobacco. Silencing the p450 gene Nitab4.5_0001439g0050 reduced tobacco resistance to bacterial wilt. These results improve our understanding of the molecular mechanism of BW resistance in tobacco and solanaceous plants.

[Identification and expression analysis of TCP family members in tobacco (Nicotiana tabacum L.)].

TCP family as plant specific transcription factor, plays an important role in different aspects of plant development. In order to screen TCP family members in tobacco, the homologous sequences of tobacco and Arabidopsis TCP family were identified by genome-wide homologous alignment. The physicochemical properties, phylogenetic relationships and cis-acting elements were analyzed by bioinformatics. The homologous genes of AtTCP3/AtTCP4 were screened, and RT-qPCR was used to detect the changes of gene expression upon 20% PEG6000 treatment. The results show that tobacco contains 63 TCP family members. Their amino acid sequence length ranged from 89 aa to 596 aa, and their protein hydropathicity grand average of hydropathicity (GRAVY) ranged from -1.147 to 0.125. The isoelectric point (pI) ranges from 4.42 to 9.94, the number of introns is 0 to 3, and the subcellular location is all located in the nucleus. The results of conserved domain and phylogenetic relationship analysis showed that the tobacco TCP family can be divided into PCF, CIN and CYC/TB1 subfamilies, and each subfamily has a stable sequence. The results of cis-acting elements in gene promoter region showed that TCP family genes contain low docile acting elements (LTR) and a variety of stress and metabolic regulation related elements (MYB, MYC). Analysis of gene expression patterns showed that AtTCP3/AtTCP4 homologous genes (NtTCP6, NtTCP28, NtTCP30, NtTCP33, NtTCP42, NtTCP57, NtTCP63) accounted for 20% PEG6000 treatment significantly up-regulated/down-regulated expression, and NtTCP30 and NtTCP57 genes were selected as candidate genes in response to drought. The results of this study analyzed the TCP family in the tobacco genome and provided candidate genes for the study of drought-resistance gene function and variety breeding in tobacco.

A dimensionality-reduction genomic prediction method without direct inverse of the genomic relationship matrix for large genomic data.

KEY MESSAGE: A new genomic prediction method (RHPP) was developed via combining randomized Haseman-Elston regression (RHE-reg), PCR based on genomic information of core population, and preconditioned conjugate gradient (PCG) algorithm. Computational efficiency is becoming a hot issue in the practical application of genomic prediction due to the large number of data generated by the high-throughput genotyping technology. In this study, we developed a fast genomic prediction method RHPP via combining randomized Haseman-Elston regression (RHE-reg), PCR based on genomic information of core population, and preconditioned conjugate gradient (PCG) algorithm. The simulation results demonstrated similar prediction accuracy between RHPP and GBLUP, and significantly higher computational efficiency of the former with the increase of individuals. The results of real datasets of both bread wheat and loblolly pine demonstrated that RHPP had a similar or better predictive accuracy in most cases compared with GBLUP. In the future, RHPP may be an attractive choice for analyzing large-scale and high-dimensional data.

Deciphering plant cell-cell communications using single-cell omics data.

Plants have various cell types that respond to different environmental factors, and cell-cell communication is the fundamental process that controls these plant responses. The emergence of single-cell techniques provides opportunities to explore features unique to each cell type and construct a comprehensive cell-cell communication (CCC) network. Although the most current successes of CCC inference were achieved in animal research, computational methods can also be directly applied to plants. This review describes the current major models for cell-cell communication inference and summarizes the computational tools based on single-cell omics datasets. In addition, we discuss the limitations of plant cell-cell communication research and propose new directions to expand the field in meaningful ways.

Analysis of herbivore-responsive long noncoding ribonucleic acids reveals a subset of small peptide-coding transcripts in Nicotiana tabacum.

Long non-coding RNAs (lncRNAs) regulate many biological processes in plants, including defense against pathogens and herbivores. Recently, many small ORFs embedded in lncRNAs have been identified to encode biologically functional peptides (small ORF-encoded peptides [SEPs]) in many species. However, it is unknown whether lncRNAs mediate defense against herbivore attack and whether there are novel functional SEPs for these lncRNAs. By sequencing Spodoptera litura-treated leaves at six time-points in Nicotiana tabacum, 22,436 lncRNAs were identified, of which 787 were differentially expressed. Using a comprehensive mass spectrometry (MS) pipeline, 302 novel SEPs derived from 115 tobacco lncRNAs were identified. Moreover, 61 SEPs showed differential expression after S. litura attack. Importantly, several of these peptides were characterized through 3D structure prediction, subcellular localization validation by laser confocal microscopy, and western blotting. Subsequent bioinformatic analysis revealed some specific chemical and physical properties of these novel SEPs, which probably represent the largest number of SEPs identified in plants to date. Our study not only identifies potential lncRNA regulators of plant response to herbivore attack but also serves as a valuable resource for the functional characterization of SEP-encoding lncRNAs.

The wilt pathogen induces different variations of root-associated microbiomes of plant.

Root-associated compartments, including the rhizosphere, rhizoplane, and endosphere, live with diverse microbial communities which profoundly affect plant growth and health. However, a systematic understanding of the microbiome assembly across the rhizosphere, rhizoplane, and endosphere under pathogen invasion remains elusive. Using 16S high-throughput sequencing, we studied how bacterial wilt disease affected the variation and assembly of the three continuous root-associated microbiomes of tobacco. The results indicated that microorganisms were gradually filtered from the rhizosphere to the endosphere. With the pathogen invasion, the rhizosphere, rhizoplane and endosphere microbiomes selected and recruited different beneficial bacterial taxa. Some recruited bacteria were also identified as keystone members in networks (i.e., Bosea in the endosphere). The microbiomes of endosphere and rhizoplane were more sensitive to plant disease than the rhizosphere microbiome. Still, response strategies of the rhizoplane and endosphere to disease were obviously different. Microbial networks of the rhizoplane became complex in diseased samples and genes involved in sporulation formation and cell cycle were enriched. However, microbial networks of the diseased endosphere were disrupted, and functional genes related to nitrogen utilization and chemotaxis were significantly increased, indicating the importance of nitrogen resources supply of plants for the endosphere microbiome under pathogen invasion. Our results provide novel insights for understanding the different responses of the root-associated microbiomes to plant disease.

Potassium accumulation characteristics and expression of related genes involved in potassium metabolism in a high-potassium variety: tobacco (Nicotiana tabacum) as a model.

We investigated potassium (K) accumulation characteristics and expression of K metabolism related genes in one high-K variety (ND202) and a common variety (NC89) of tobacco (Nicotiana tabacum L.). Results showed that K accumulation and leaf K content in ND202 were higher than those in NC89. The distribution rate and K accumulation in the leaves of ND202 increased significantly, while the distribution rate in the roots and stems had lower values. In addition, the maximum K accumulation rate and high-speed K accumulation duration in ND202 were found to be better than those in NC89. The expression of NKT1 in the upper and middle leaves of ND202 had an advantage, and the relative expression of NtKC1 and NtTPK1 in both the upper and middle leaves, as well as the roots, was also significantly upregulated. Conversely, the expression of NTRK1 in the lower leaves and roots of ND202 was weaker. ND202 had significantly greater expression levels of NtHAK1 than NC89 in the upper and middle leaves and roots; moreover, the expression of NtKT12 in the upper leaves and roots of ND202 was also higher. In comparison with common varieties, high-K varieties had a stronger ability to absorb and accumulate K. They also possessed higher expression of K+ channel- and transporter-related genes and showed a superior K accumulation rate and longer duration of high-speed K accumulation. Furthermore, K accumulation rate at 40-60days can be suggested as an important reference for the selection of high-K tobacco varieties.

Genome-wide association study revealed TaHXK3-2A as a candidate gene controlling stomatal index in wheat seedlings.

Stomata are important channels for the control of gas exchange between plants and the atmosphere. To examine the genetic architecture of wheat stomatal index, we performed a genome-wide association study (GWAS) using a panel of 539 wheat accessions and 450 678 polymorphic single nucleotide polymorphisms (SNPs) that were detected using wheat-specific 660K SNP array. A total of 130 SNPs were detected to be significantly associated with stomatal index in both leaf surfaces of wheat seedlings. These significant SNPs were distributed across 16 chromosomes and involved 2625 candidate genes which participate in stress response, metabolism and cell/organ development. Subsequent bulk segregant analysis (BSA), combined with GWAS identified one major haplotype on chromosome 2A, that is responsible for stomatal index on the abaxial leaf surface. Candidate gene association analysis revealed that genetic variation in the promoter region of the hexokinase gene TaHXK3-2A was significantly associated with the stomatal index. Moreover, transgenic analysis confirmed that TaHXK3-2A overexpression in wheat decreased the size of leaf pavement cells but increased stomatal density through the glucose metabolic pathway, resulting in drought sensitivity among TaHXK3-2A transgenic lines due to an increased transpiration rate. Taken together, these results provide valuable insights into the genetic control of the stomatal index in wheat seedlings.

Tobacco Root Microbial Community Composition Significantly Associated With Root-Knot Nematode Infections: Dynamic Changes in Microbiota and Growth Stage.

The root-knot nematode (RKN) is an important pathogen that affects the growth of many crops. Exploring the interaction of biocontrol bacteria-pathogens-host root microbes is the theoretical basis for improving colonization and controlling the effect of biocontrol bacteria in the rhizosphere. Therefore, 16S and 18S rRNA sequencing technology was used to explore the microbial composition and diversity of tobacco roots (rhizosphere and endophytic) at different growth stages in typical tobacco RKN-infected areas for 2 consecutive years. We observed that RKN infection changed the α-diversity and microbial composition of root microorganisms and drove the transformation of microorganisms from bacteria to fungi. The abundance of Sphingomonas decreased significantly from 18% to less than 3%, while the abundance of Rhizobiaceae increased from 4 to 15% at the early growth stage during the first planting year, and it promoted the proliferation of Chryseobacterium at the late growth stage in rhizosphere microorganisms with the highest abundance of 17%. The overall trend of rhizosphere microorganisms changed in the early growth stage with increasing growth time. The specific results were as follows: (1) Rhizobiaceae and Chryseobacterium increased rapidly after 75 days, became the main abundant bacteria in the rhizosphere microorganisms. (2) The dominant flora in fungi were Fusarium and Setophoma. (3) Comparing the root microbes in 2017 and 2018, RKN infection significantly promoted the proliferation of Pseudomonas and Setophoma in both the rhizosphere and endophytes during the second year of continuous tobacco planting, increasing the relative abundance of Pseudomonas from 2 to 25%. Pseudomonas was determined to play an important role in plant pest control. Finally, a total of 32 strains of growth-promoting bacteria were screened from tobacco rhizosphere bacteria infected with RKN through a combination of 16S rRNA sequencing and life-promoting tests. The results of this research are helpful for analyzing the relationship between RKNs and bacteria in plants, providing reference data for elucidating the pathogenesis of RKNs and new ideas for the biological control of RKNs. GRAPHICAL ABSTRACT.

Nitrogen application and differences in leaf number retained after topping affect the tobacco (Nicotiana tabacum) transcriptome and metabolome.

BACKGROUND: Agronomic treatments such as the application of nitrogen fertilizer and topping (removal of the inflorescence and top leaves) cause substantial changes in plant metabolism. To explore these changes, we conducted comparative transcriptomic and metabolomic analyses of leaves collected from four positions along the stem on plants exposed to two nitrogen doses and with different numbers of leaves retained after topping in tobacco (Nicotiana tabacum). RESULTS: We identified 13,330 unique differentially expressed genes and 32 differentially abundant metabolites. Through RNA-seq and WGCNA analyze, we constructed 2 co-expression networks (green and blue) highly correlation to N application and leaf number retained, predicted a hub gene NtGER3 may play an important role in N metabolism related to amino acid (cysteine) through CK pathway in tobacco leaves, NtARFs may participated in modulating the auxin signal and N in bottom leaves and NtRAP2.12 as key gene involved in N regulation by ethylene pathway. What's more, our data prove C/N transformation and balance affect the "source - flow - sink" redistribution and remobilization in tobacco during growth and development process. CONCLUSIONS: Overall, this comparative transcriptomics study provides novel insight into the complex molecular mechanisms underlying plant responses to different levels of nitrogen application and the number of leaves remaining after topping in plants.

PCMDB: a curated and comprehensive resource of plant cell markers.

The advent of single-cell sequencing opened a new era in transcriptomic and genomic research. To understand cell composition using single-cell studies, a variety of cell markers have been widely used to label individual cell types. However, the specific database of cell markers for use by the plant research community remains very limited. To overcome this problem, we developed the Plant Cell Marker DataBase (PCMDB, http://www.tobaccodb.org/pcmdb/), which is based on a uniform annotation pipeline. By manually curating over 130 000 research publications, we collected a total of 81 117 cell marker genes of 263 cell types in 22 tissues across six plant species. Tissue- and cell-specific expression patterns can be visualized using multiple tools: eFP Browser, Bar, and UMAP/TSNE graph. The PCMDB also supports several analysis tools, including SCSA and SingleR, which allows for user annotation of cell types. To provide information about plant species currently unsupported in PCMDB, potential marker genes for other plant species can be searched based on homology with the supported species. PCMDB is a user-friendly hierarchical platform that contains five built-in search engines. We believe PCMDB will constitute a useful resource for researchers working on cell type annotation and the prediction of the biological function of individual cells.

Comparative transcriptomic analysis of seed coats with high and low lignin contents reveals lignin and flavonoid biosynthesis in Brassica napus.

BACKGROUND: Brassica napus L. (2n = 38, AACC) is one of the most important oil crops and sources of protein for animal feed worldwide. Lignin is a large molecule aromatic polymer and a major cell wall component. However, lignin in the seed coat reduces the availability and restricts the development of rapeseed cake. Therefore, it is critical to reduce the lignin content of the seed coat. Here, high-lignin (H-lignin) and low-lignin (L-lignin) content recombinant inbred lines (RILs) were selected from an RIL population for analysis. RESULTS: The cross-section results indicated that the seed coat of the H-lignin lines was thicker than that of the L-lignin lines, especially the palisade layer. The seed coats and embryos at 35, 40 and 46 days after flowering (DAF) were subjected to RNA sequencing (RNA-Seq), and the expression of the BnPAL and BnC4H gene families in the lignin pathway was significantly higher in the H-lignin seed coat than in the L-lignin seed coat. The Bn4CL gene family also showed this trend. In addition, among the genes related to plant hormone synthesis, BnaC02g01710D was upregulated and BnaA07g11700D and BnaC09g00190D were downregulated in H-lignin lines. Some transcription factors were upregulated, such as BnNAC080, BnNAC083, BnMYB9, BnMYB9-1, BnMYB60 and BnMYB60-1, while BnMYB91 was downregulated in H-lignin lines. Moreover, most genes of the flavonoid pathway, such as BnCHS and BnDFR, were strongly expressed in H-lignin seed coat. CONCLUSIONS: In Our study, some key genes such as hormone synthesis genes, transcription factors and miRNAs related to lignin and flavonoid biosynthesis were identified. A regulatory model of B. napus seed coat lignin was proposed. These results provide new insight into lignin and flavonoid biosynthesis in B. napus.

PLncDB V2.0: a comprehensive encyclopedia of plant long noncoding RNAs.

Long noncoding RNAs (lncRNAs) are transcripts longer than 200 nucleotides with little or no protein coding potential. The expanding list of lncRNAs and accumulating evidence of their functions in plants have necessitated the creation of a comprehensive database for lncRNA research. However, currently available plant lncRNA databases have some deficiencies, including the lack of lncRNA data from some model plants, uneven annotation standards, a lack of visualization for expression patterns, and the absence of epigenetic information. To overcome these problems, we upgraded our Plant Long noncoding RNA Database (PLncDB, http://plncdb.tobaccodb.org/), which was based on a uniform annotation pipeline. PLncDB V2.0 currently contains 1 246 372 lncRNAs for 80 plant species based on 13 834 RNA-Seq datasets, integrating lncRNA information from four other resources including EVLncRNAs, RNAcentral and etc. Expression patterns and epigenetic signals can be visualized using multiple tools (JBrowse, eFP Browser and EPexplorer). Targets and regulatory networks for lncRNAs are also provided for function exploration. In addition, PLncDB V2.0 is hierarchical and user-friendly and has five built-in search engines. We believe PLncDB V2.0 is useful for the plant lncRNA community and data mining studies and provides a comprehensive resource for data-driven lncRNA research in plants.

A large-scale genomic association analysis identifies the candidate causal genes conferring stripe rust resistance under multiple field environments.

The incorporation of resistance genes into wheat commercial varieties is the ideal strategy to combat stripe or yellow rust (YR). In a search for novel resistance genes, we performed a large-scale genomic association analysis with high-density 660K single nucleotide polymorphism (SNP) arrays to determine the genetic components of YR resistance in 411 spring wheat lines. Following quality control, 371 972 SNPs were screened, covering over 50% of the high-confidence annotated gene space. Nineteen stable genomic regions harbouring 292 significant SNPs were associated with adult-plant YR resistance across nine environments. Of these, 14 SNPs were localized in the proximity of known loci widely used in breeding. Obvious candidate SNP variants were identified in certain confidence intervals, such as the cloned gene Yr18 and the major locus on chromosome 2BL, despite a large extent of linkage disequilibrium. The number of causal SNP variants was refined using an independent validation panel and consideration of the estimated functional importance of each nucleotide polymorphism. Interestingly, four natural polymorphisms causing amino acid changes in the gene TraesCS2B01G513100 that encodes a serine/threonine protein kinase (STPK) were significantly involved in YR responses. Gene expression and mutation analysis confirmed that STPK played an important role in YR resistance. PCR markers were developed to identify the favourable TraesCS2B01G513100 haplotype for marker-assisted breeding. These results demonstrate that high-resolution SNP-based GWAS enables the rapid identification of putative resistance genes and can be used to improve the efficiency of marker-assisted selection in wheat disease resistance breeding.

Comparative transcriptome and metabolomic profiling reveal the complex mechanisms underlying the developmental dynamics of tobacco leaves.

Although the leaf is the most important photosynthetic organ in most plants, many of the molecular mechanisms underlying leaf developmental dynamics remain to be explored. To better understand the transcriptional regulatory mechanisms involved in leaf development, we conducted comparative transcriptomic and metabolomic analysis of leaves from seven positions on tobacco (Nicotiana tabacum) plants. A total of 35,622 unique differentially expressed genes and 79 metabolites were identified. A time-series expression analysis detected two interesting transcriptional profiles, one comprising 10,197 genes that displayed continual up-regulation during leaf development and another comprising 4696 genes that displayed continual down-regulation. Combining these data with co-expression network results identified four important regulatory networks involved in photorespiration and the tricarboxylic acid cycle; these networks may regulate carbon/nitrogen balance during leaf development. We also found that the transcription factor NtGATA5 acts as a hub associated with C and N metabolism and chloroplast development during leaf development through regulation of phytohormones. Furthermore, we investigated the transcriptional dynamics of genes involved in the auxin, cytokinin, and jasmonic acid biosynthesis and signaling pathways during tobacco leaf development. Overall, our study greatly expands the understanding of the regulatory network controlling developmental dynamics in plant leaves.

Association Analysis Identifies New Loci for Resistance to Chinese Yr26-Virulent Races of the Stripe Rust Pathogen in a Diverse Panel of Wheat Germplasm.

Stripe rust caused by Puccinia striiformis f. sp. tritici (Pst) is one of the most destructive fungal diseases of wheat worldwide. The expanding Yr26-virulent Pst race (V26) group overcomes almost all currently deployed resistance genes in China and has continued to accumulate new virulence. Investigating the genetic architecture of stripe rust resistance in common wheat is an important basis for a successful utilization of resistance in breeding programs. A panel of 410 exotic wheat germplasms was used for characterizing new stripe rust resistance loci. This panel was genotyped using high-density wheat 660K single-nucleotide polymorphism (SNP) array, and phenotypic evaluation of seedlings for stripe rust resistance was performed using multiple Pst races. Thirty-five loci conferring resistance were identified through genome-wide association mapping, and explained phenotypic variances ranged from 53 to 75%. Of these, 14 were colocated in the proximity of the known loci, including cataloged Yr genes Yr9, Yr10, Yr26, Yr33, Yr47, Yr56, Yr57, Yr64, Yr67, Yr72, and Yr81 and three temporarily designated as YrCen, YrNP63, and YrRC detected in our quantitative trait locus (QTL) mapping studies. Seven of them (Yr9, Yr10, Yr24/26, Yr81, YrCEN, YrNP63, and YrRC) were confirmed by molecular detection or genetic analysis. New loci that were identified to be different from reported Yr genes need further confirmation. Nine QTL with significantly large phenotypic effect on resistance to all tested races were considered as major loci for effective resistance. The identified loci enrich our stripe rust resistance gene pool, and the linked SNPs should be useful for marker-assisted selection in breeding programs.

Genome-Scale Analysis of Homologous Genes among Subgenomes of Bread Wheat (Triticum aestivum L.).

Determining the distribution and correspondence of genome-scale homologous genes in wheat are effective ways to uncover chromosome rearrangement that has occurred during crop evolution and domestication, which can contribute to improvements in crop breeding. High-resolution and comprehensive analysis of the wheat genome by the International Wheat Genome Sequencing Consortium (IWGSC) revealed a total of 88,733 high-confidence homologous genes of four major types (1:1:1, 1:1:0, 0:1:1 and 1:0:1) among the A, B and D subgenomes of wheat. This data was used to compare homologous gene densities among chromosomes, clarify their distribution and correspondence relationship, and compare their functional enrichment. The average density of 1:1:1 homologous genes was about 10 times more than the density of the other three types of homologous genes, although the homologous gene densities of the various chromosomes were similar within each homologous type. Three regions of exceptional density were detected in 1:1:1 homologous genes, the isolate peak on the tail of chromosome 4A, and the desert regions at the start of chromosome 7A and 7D. The correspondence between homologous genes of the wheat subgenomes demonstrated translocation between the tail segments of chromosome 4A and 5A, and the inversion of the segment of original 5A and 7B into the tail of 4A. The homologous genes on the inserting segments of 5A and 7B to 4A were highly enriched in nitrogen, primary metabolite and small molecular metabolism processes, compared with genes on other regions of the original 4A chromosome. This study provides a refined genome-scale reference of homologous genes for wheat molecular research and breeding, which will help to broaden the application of the wheat genome and can be used as a template for research on other polyploid plants.

Activation of potassium released from soil by root-secreted organic acids in different varieties of tobacco (Nicotiana tabacum).

Organic acids secreted from the roots of plants play important roles in nutrient acquisition and metal detoxification; however, the precise underlying mechanisms of these processes remain poorly understood. In the present study we examined the content of organic acids exuded from roots and the effects of these organic acids on the activation of slowly available potassium (K) at different K levels, including normal K supply and K-deficient conditions. In addition, the study system also comprised a high-K tobacco variety (ND202) and two common ones (K326 and NC89). Our results showed that high-K varieties exhibited significantly higher contents of organic acids in its root exudates and available K in both rhizosphere and non-rhizosphere soils than the other varieties. This research also suggested that a cyclic process in which soil was acidified after being complexed by organic acids was involved in the release of slowly available K, and that this process primarily depended on the soil pH at high organic acids concentrations, but the complexation of organic ligands became dominant at low concentrations. In conclusion, tobacco roots secrete organic acids to increase available K content and improve the utilisation rate of soil K. High-K varieties probably enhance slowly available K activation by secreting relatively high amounts of organic acids, thus leading to more available K in soil for absorption by plants.

Systematic study of the stress-responsive Rboh gene family in Nicotiana tabacum: Genome-wide identification, evolution and role in disease resistance.

Plant respiratory burst oxidase homolog (Rboh) gene family encodes the key enzymatic subunits of reactive oxygen species (ROS) production pathways, and play crucial role in plant signaling, development and stress responses. In present work, twenty genes were identified in Nicotiana tabacum Rboh family (NtabRboh) and classified into four phylogenetic groups (I-IV). Fourteen NtabRboh genes were positioned on ten chromosomes (i.e., Ch1, 2, 4, 7-11, 14 and 21), and six scaffolds. Synteny and evolutionary analysis showed that most of the NtabRboh genes have evolved from the genomes of the ancestor species (N. tomentosiformis and N. sylvestris), which afterwards expanded through duplication events. The promoter regions of the NtabRboh genes contained numerous cis-acting regulatory elements for hormones, plant growth, and different biotic and abiotic factors. The NtabRbohF gene transcript comprised target sites for wounding and stress responsive microRNAs: nta-miR166a-d, g and h. The transcript abundance of NtabRboh genes in different tissues reflected their important for plant growth and organ development in tobacco. RT-qPCR-assays demonstrated that the expression of NtabRboh genes are regulated by viral and bacterial pathogens, drought, cold and cadmium stress. The expression levels NtabRbohA, B and C were significantly up-regulated in "black shank and tobacco mosaic virus-inoculated susceptible and transgenic tobacco cultivars, showing that these genes play important roles in disease resistance.