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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

iTRAQ-based quantitative proteomic analysis reveals proteomic changes in leaves of cultivated tobacco (Nicotiana tabacum) in response to drought stress.

Drought is one of the most severe forms of abiotic stresses that threaten the survival of plants, including crops. In turn, plants dramatically change their physiology to increase drought tolerance, including reconfiguration of proteomes. Here, we studied drought-induced proteomic changes in leaves of cultivated tobacco (Nicotiana tabacum), a solanaceous plant, using the isobaric tags for relative and absolute quantitation (iTRAQ)-based protein labeling technology. Of identified 5570 proteins totally, drought treatment increased and decreased abundance of 260 and 206 proteins, respectively, compared with control condition. Most of these differentially regulated proteins are involved in photosynthesis, metabolism, and stress and defense. Although abscisic acid (ABA) levels greatly increased in drought-treated tobacco leaves, abundance of detected ABA biosynthetic enzymes showed no obvious changes. In contrast, heat shock proteins (HSPs), thioredoxins, ascorbate-, glutathione-, and hydrogen peroxide (H2O2)-related proteins were up- or down-regulated in drought-treated tobacco leaves, suggesting that chaperones and redox signaling are important for tobacco tolerance to drought, and it is likely that redox-induced posttranslational modifications play an important role in modulating protein activity. This study not only provides a comprehensive dataset on overall protein changes in drought-treated tobacco leaves, but also shed light on the mechanism by which solanaceous plants adapt to drought stress.

Regulation of Nicotine Biosynthesis by an Endogenous Target Mimicry of MicroRNA in Tobacco.

The interaction between noncoding endogenous target mimicry (eTM) and its corresponding microRNA (miRNA) is a newly discovered regulatory mechanism and plays pivotal roles in various biological processes in plants. Tobacco (Nicotiana tabacum) is a model plant for studying secondary metabolite alkaloids, of which nicotine accounts for approximately 90%. In this work, we identified four unique tobacco-specific miRNAs that were predicted to target key genes of the nicotine biosynthesis and catabolism pathways and an eTM, novel tobacco miRNA (nta)-eTMX27, for nta-miRX27 that targets QUINOLINATE PHOSPHORIBOSYLTRANSFERASE2 (QPT2) encoding a quinolinate phosphoribosyltransferase. The expression level of nta-miRX27 was significantly down-regulated, while that of QPT2 and nta-eTMX27 was significantly up-regulated after topping, and consequently, nicotine content increased in the topping-treated plants. The topping-induced down-regulation of nta-miRX27 and up-regulation of QPT2 were only observed in plants with a functional nta-eTMX27 but not in transgenic plants containing an RNA interference construct targeting nta-eTMX27. Our results demonstrated that enhanced nicotine biosynthesis in the topping-treated tobacco plants is achieved by nta-eTMX27-mediated inhibition of the expression and functions of nta-miRX27. To our knowledge, this is the first report about regulation of secondary metabolite biosynthesis by an miRNA-eTM regulatory module in plants.

Large-scale development of SSR markers in tobacco and construction of a linkage map in flue-cured tobacco.

Tobacco (Nicotiana tabacum L.), particularly flue-cured tobacco, is one of the most economically important nonfood crops and is also an important model system in plant biotechnology. Despite its importance, only limited molecular marker resources are available for genome analysis, genetic mapping, and breeding. Simple sequence repeats (SSR) are one of the most widely-used molecular markers, having significant advantages including that they are generally co-dominant, easy to use, abundant in eukaryotic organisms, and produce highly reproducible results. In this study, based on the genome sequence data of flue-cured tobacco (K326), we developed a total of 13,645 mostly novel SSR markers, which were working in a set of eighteen tobacco varieties of four different types. A mapping population of 213 backcross (BC1) individuals, which were derived from an intra-type cross between two flue-cured tobacco varieties, Y3 and K326, was selected for mapping. Based on the newly developed SSR markers as well as published SSR markers, we constructed a genetic map consisting of 626 SSR loci distributed across 24 linkage groups and covering a total length of 1120.45 cM with an average distance of 1.79 cM between adjacent markers, which is the highest density map of flue-cured tobacco till date.

A diverse set of miRNAs responsive to begomovirus-associated betasatellite in Nicotiana benthamiana.

BACKGROUND: Roles of microRNAs (miRNAs) and short interfering RNAs (siRNAs) in biotic stress responses, e.g., viral infection, have been demonstrated in plants by many studies. Tomato yellow leaf curl China virus (TYLCCNV) is a monopartite begomovirus that can systemically infect Solanaceae plants, and induces leaf curling, yellowing and enation symptoms when co-inoculated with a betasatellite (TYLCCNB). The released genome sequence of Nicotiana benthamiana provides an opportunity to identify miRNAs and siRNAs responsive to begomovirus-associated betasatellite in N. benthamiana. RESULTS: miRNAs were identified in three small RNA libraries generated using RNA isolated from N. benthamiana plants systemically infected with TYLCCNV (Y10A) alone, co-infected with Y10A and its betasatellite TYLCCNB (Y10ß) or a TYLCCNB mutant (Y10mß) that contains a mutated ßC1, the sole betasatellite-encoded protein. A total of 196 conserved miRNAs from 38 families and 197 novel miRNAs from 160 families were identified. Northern blot analysis confirmed that expression of species-specific miRNAs was much lower than that of conserved miRNAs. Several conserved and novel miRNAs were found to be responsive to co-infection of Y10A and Y10ß but not to co-infection of Y10A and Y10mß, suggesting that these miRNAs might play a role unique to interaction between Y10ß and N. benthamiana. Additionally, we identified miRNAs that can trigger the production of phased secondary siRNAs (phasiRNAs). CONCLUSIONS: Identification of miRNAs with differential expression profiles in N. benthamiana co-infected with Y10A and Y10ß and co-infected with Y10A and Y10mß indicates that these miRNAs are betasatellite-responsive. Our result also suggested a potential role of miRNA-mediated production of phasiRNAs in interaction between begomovirus and N. benthamiana.

Alteration of the alkaloid profile in genetically modified tobacco reveals a role of methylenetetrahydrofolate reductase in nicotine N-demethylation.

Methylenetetrahydrofolate reductase (MTHFR) is a key enzyme of the tetrahydrofolate (THF)-mediated one-carbon (C1) metabolic network. This enzyme catalyzes the reduction of 5,10-methylene-THF to 5-methyl-THF. The latter donates its methyl group to homocysteine, forming methionine, which is then used for the synthesis of S-adenosyl-methionine, a universal methyl donor for numerous methylation reactions, to produce primary and secondary metabolites. Here, we demonstrate that manipulating tobacco (Nicotiana tabacum) MTHFR gene (NtMTHFR1) expression dramatically alters the alkaloid profile in transgenic tobacco plants by negatively regulating the expression of a secondary metabolic pathway nicotine N-demethylase gene, CYP82E4. Quantitative real-time polymerase chain reaction and alkaloid analyses revealed that reducing NtMTHFR expression by RNA interference dramatically induced CYP82E4 expression, resulting in higher nicotine-to-nornicotine conversion rates. Conversely, overexpressing NtMTHFR1 suppressed CYP82E4 expression, leading to lower nicotine-to-nornicotine conversion rates. However, the reduced expression of NtMTHFR did not affect the methionine and S-adenosyl-methionine levels in the knockdown lines. Our finding reveals a new regulatory role of NtMTHFR1 in nicotine N-demethylation and suggests that the negative regulation of CYP82E4 expression may serve to recruit methyl groups from nicotine into the C1 pool under C1-deficient conditions.

Establishment and application of Agrobacterium-delivered CRISPR/Cas9 system for wild tobacco (Nicotiana alata) genome editing.

Clustered regularly interspaced short palindromic repeats (CRISPR)-associated protein 9 (CRISPR-Cas9) system has been widely applied in cultivated crops, but limited in their wild relatives. Nicotiana alata is a typical wild species of genus Nicotiana that is globally distributed as a horticultural plant and well-studied as a self-incompatibility model. It also has valuable genes for disease resistance and ornamental traits. However, it lacks an efficient genetic transformation and genome editing system, which hampers its gene function and breeding research. In this study, we developed an optimized hypocotyl-mediated transformation method for CRISPR-Cas9 delivery. The genetic transformation efficiency was significantly improved from approximately 1% to over 80%. We also applied the CRISPR-Cas9 system to target the phytoene desaturase (NalaPDS) gene in N. alata and obtained edited plants with PDS mutations with over 50% editing efficiency. To generate self-compatible N. alata lines, a polycistronic tRNA-gRNA (PTG) strategy was used to target exonic regions of allelic S-RNase genes and generate targeted knockouts simultaneously. We demonstrated that our system is feasible, stable, and high-efficiency for N. alata genome editing. Our study provides a powerful tool for basic research and genetic improvement of N. alata and an example for other wild tobacco species.

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

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

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

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

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

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

Identification of wounding and topping responsive small RNAs in tobacco (Nicotiana tabacum).

BACKGROUND: MicroRNAs (miRNAs) and short interfering RNAs (siRNAs) are two major classes of small RNAs. They play important regulatory roles in plants and animals by regulating transcription, stability and/or translation of target genes in a sequence-complementary dependent manner. Over 4,000 miRNAs and several classes of siRNAs have been identified in plants, but in tobacco only computational prediction has been performed and no tobacco-specific miRNA has been experimentally identified. Wounding is believed to induce defensive response in tobacco, but the mechanism responsible for this response is yet to be uncovered. RESULTS: To get insight into the role of small RNAs in damage-induced responses, we sequenced and analysed small RNA populations in roots and leaves from wounding or topping treated tobacco plants. In addition to confirmation of expression of 27 known miRNA families, we identified 59 novel tobacco-specific miRNA members of 38 families and a large number of loci generating phased 21- or 24-nt small RNAs (including ta-siRNAs). A number of miRNAs and phased small RNAs were found to be responsive to wounding or topping treatment. Targets of small RNAs were further surveyed by degradome sequencing. CONCLUSIONS: The expression changes of miRNAs and phased small RNAs responsive to wounding or topping and identification of defense related targets for these small RNAs suggest that the inducible defense response in tobacco might be controlled by pathways involving small RNAs.

Screening of Tobacco Genotypes for Phytophthora nicotianae Resistance.

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