Root Transcriptome Analysis of Wild Peanut Reveals Candidate Genes for Nematode Resistance

Overexpression of NtEXPA7 promotes seedling growth and resistance to root-knot nematode in tobacco

Root-knot nematode (RKN) is one of the most damaging plant pathogen in the world. They exhibit a wide host range and cause serious crop losses. The cell wall, encasing every plant cell, plays a crucial role in defending of RKN invasion. Expansins are a group of cell wall proteins inducing cell wall loosening and extensibility. They are widely involved in the regulation of plant growth and the response to biotic and abiotic stresses. In this study, we have characterized the biological function of tobacco (Nicotiana tabacum) NtEXPA7, the homologue of Solyc08g080060.2 (SlEXPA18), of which the transcription level was significantly reduced in susceptible tomato upon RKN infection. The expression of NtEXPA7 was up-regulated after inoculation of RKNs. The NtEXPA7 protein resided in the cell wall. Overexpression of NtEXPA7 promoted the seedling growth of transgenic tobacco. Meanwhile the increased expression of NtEXPA7 was beneficial to enhance the resistance against RKNs. This study expands the understanding of biological role of expansin in coordinate plant growth and disease resistance.

Advances in the evolution research and genetic breeding of peanut

Peanut is an important cash crop used in oil, food and feed in our country. The rapid development of sequencing technology has promoted the research on the related aspects of peanut genetic breeding. This paper reviews the research progress of peanut origin and evolution, genetic breeding, molecular markers and their applications, genomics, QTL mapping and genome selection techniques. The main problems of molecular genetic breeding in peanut research worldwide include: the narrow genetic resources of cultivated species, unstable genetic transformation and unclear molecular mechanism of important agronomic traits. Considering the severe challenges regarding the supply of edible oil, and the main problems in peanut production, the urgent research directions of peanut are put forward: The de novo domestication and the exploitation of excellent genes from wild resources to improve modern cultivars; Integration of multi-omics data to enhance the importance of big data in peanut genetics and breeding; Cloning the important genes related to peanut agronomic traits and analyzing their fine regulation mechanisms; Precision molecular design breeding and using gene editing technology to accurately improve the key traits of peanut.

The Stilbene Synthase Family in Arachis: A Genome-Wide Study and Functional Characterization in Response to Stress

Peanut (Arachis hypogaea) and its wild relatives are among the few species that naturally synthesize resveratrol, a well-known stilbenoid phytoalexin that plays a crucial role in plant defense against biotic and abiotic stresses. Resveratrol has received considerable attention due to its health benefits, such as preventing and treating various human diseases and disorders. Chalcone (CHS) and Stilbene (STS) Synthases are plant-specific type III Polyketide Synthases (PKSs) that share the same substrates and are key branch enzymes in the biosynthesis of flavonoids and stilbenoids, respectively. Although resveratrol accumulation in response to external stimulus has been described in peanut, there are no comprehensive studies of the CHS and STS gene families in the genus Arachis. In the present study, we identified and characterized 6 CHS and 46 STS genes in the tetraploid peanut and an average of 4 CHS and 22 STS genes in three diploid wild species (Arachis duranensis, Arachis ipaënsis and Arachis stenosperma). The CHS and STS gene and protein structures, chromosomal distributions, phylogenetic relationships, conserved amino acid domains, and cis-acting elements in the promoter regions were described for all Arachis species studied. Based on gene expression patterns of wild A. stenosperma STS genes in response to different biotic and abiotic stresses, we selected the candidate AsSTS4 gene, which is strongly induced by ultraviolet (UV) light exposure, for further functional investigation. The AsSTS4 overexpression in peanut hairy roots significantly reduced (47%) root-knot nematode infection, confirming that stilbene synthesis activation in transgenic plants can increase resistance to pathogens. These findings contribute to understanding the role of resveratrol in stress responses in Arachis species and provide the basis for genetic engineering for improved production of valuable secondary metabolites in plants.

Genome-Wide Identification and Characterization of the Phytochrome Gene Family in Peanut

To investigate the potential role of phytochrome (PHY) in peanut growth and its response to environmental fluctuations, eight candidate AhPHY genes were identified via genome-wide analysis of cultivated peanut. These AhPHY polypeptides were determined to possess acidic and hydrophilic physiochemical properties and exhibit subcellular localization patterns consistent with residence in the nucleus and cytoplasm. Phylogenetic analysis revealed that the AhPHY gene family members were classified into three subgroups homologous to the PHYA/B/E progenitors of Arabidopsis. AhPHY genes within the same clade largely displayed analogous gene structure, conserved motifs, and phosphorylation sites. AhPHY exhibited symmetrical distribution across peanut chromosomes, with 7 intraspecific syntenic gene pairs in peanut, as well as 4 and 20 interspecific PHY syntenic gene pairs in Arabidopsis and soybean, respectively. A total of 42 cis-elements were predicted in AhPHY promoters, including elements implicated in phytohormone regulation, stress induction, physiology, and photoresponse, suggesting putative fundamental roles across diverse biological processes. Moreover, spatiotemporal transcript profiling of AhPHY genes in various peanut tissues revealed distinct expression patterns for each member, alluding to putative functional specialization. This study contributes novel insights into the classification, structure, molecular evolution, and expression profiles of the peanut phytochrome gene family, and also provides phototransduction gene resources for further mechanistic characterization.

TETRASPANIN 8-1 from Phaseolus vulgaris plays a key role during mutualistic interactions

Arbuscular mycorrhizal (AM) fungi and rhizobia form two of the most important plant-microbe associations for the assimilation of phosphorus (P) and nitrogen (N). Symbiont-derived signals are able to coordinate the infection process by triggering multiple responses in the plant root, such as calcium influxes and oscillations, increased reactive oxygen species (ROS), cytoskeletal rearrangements and altered gene expression. An examination was made of the role of tetraspanins, which are transmembrane proteins that self-organize into tetraspanin web regions, where they recruit specific proteins into platforms required for signal transduction, membrane fusion, cell trafficking, and ROS generation. In plant cells, tetraspanins are scaffolding proteins associated with root radial patterning, biotic and abiotic stress responses, cell fate determination, plasmodesmata and hormonal regulation. Some plant tetraspanins, such as Arabidopsis thaliana TETRASPANIN 8 and TETRASPANIN 9 (AtTET8 and AtTET9) are associated with exosomes during inter-kingdom communication. In this study, a homolog of AtTET8, PvTET8-1, in common bean (Phaseolus vulgaris L. var. Negro Jamapa) was examined in roots during interactions with Rhizobium tropici and Rhizophagus irregularis. The promoter of PvTET8-1 contained several cis-acting regulatory DNA elements potentially related to mutualistic interactions, and PvTET8-1 was transcriptionally activated during AM fungal and rhizobial associations. Silencing it decreased the size and number of nodules, nitrogen fixation, and mycorrhizal arbuscule formation, whereas overexpressing it increased the size and number of nodules, and mycorrhizal arbuscule formation but decreased nitrogen fixation. PvTET8-1 appears to be an important element in both of these mutualistic interactions, perhaps through its interaction with NADPH oxidase and the generation of ROS during the infection processes.

Biotechnological Tools to Elucidate the Mechanism of Plant and Nematode Interactions

Plant-parasitic nematodes (PPNs) pose a threat to global food security in both the developed and developing worlds. PPNs cause crop losses worth a total of more than USD 150 billion worldwide. The sedentary root-knot nematodes (RKNs) also cause severe damage to various agricultural crops and establish compatible relationships with a broad range of host plants. This review aims to provide a broad overview of the strategies used to identify the morpho-physiological and molecular events that occur during RKN parasitism. It describes the most current developments in the transcriptomic, proteomic, and metabolomic strategies of nematodes, which are important for understanding compatible interactions of plants and nematodes, and several strategies for enhancing plant resistance against RKNs. We will highlight recent rapid advances in molecular strategies, such as gene-silencing technologies, RNA interference (RNAi), and small interfering RNA (siRNA) effector proteins, that are leading to considerable progress in understanding the mechanism of plant-nematode interactions. We also take into account genetic engineering strategies, such as targeted genome editing techniques, the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR associated protein 9 (Cas9) (CRISPR/Cas-9) system, and quantitative trait loci (QTL), to enhance the resistance of plants against nematodes.

Advances in omics research on peanut response to biotic stresses

Peanut growth, development, and eventual production are constrained by biotic and abiotic stresses resulting in serious economic losses. To understand the response and tolerance mechanism of peanut to biotic and abiotic stresses, high-throughput Omics approaches have been applied in peanut research. Integrated Omics approaches are essential for elucidating the temporal and spatial changes that occur in peanut facing different stresses. The integration of functional genomics with other Omics highlights the relationships between peanut genomes and phenotypes under specific stress conditions. In this review, we focus on research on peanut biotic stresses. Here we review the primary types of biotic stresses that threaten sustainable peanut production, the multi-Omics technologies for peanut research and breeding, and the recent advances in various peanut Omics under biotic stresses, including genomics, transcriptomics, proteomics, metabolomics, miRNAomics, epigenomics and phenomics, for identification of biotic stress-related genes, proteins, metabolites and their networks as well as the development of potential traits. We also discuss the challenges, opportunities, and future directions for peanut Omics under biotic stresses, aiming sustainable food production. The Omics knowledge is instrumental for improving peanut tolerance to cope with various biotic stresses and for meeting the food demands of the exponentially growing global population.

Engineering Resistance against Sclerotinia sclerotiorum Using a Truncated NLR (TNx) and a Defense-Priming Gene

The association of both cell-surface PRRs (Pattern Recognition Receptors) and intracellular receptor NLRs (Nucleotide-Binding Leucine-Rich Repeat) in engineered plants have the potential to activate strong defenses against a broad range of pathogens. Here, we describe the identification, characterization, and in planta functional analysis of a novel truncated NLR (TNx) gene from the wild species Arachis stenosperma (AsTIR19), with a protein structure lacking the C-terminal LRR (Leucine Rich Repeat) domain involved in pathogen perception. Overexpression of AsTIR19 in tobacco plants led to a significant reduction in infection caused by Sclerotinia sclerotiorum, with a further reduction in pyramid lines containing an expansin-like B gene (AdEXLB8) potentially involved in defense priming. Transcription analysis of tobacco transgenic lines revealed induction of hormone defense pathways (SA; JA-ET) and PRs (Pathogenesis-Related proteins) production. The strong upregulation of the respiratory burst oxidase homolog D (RbohD) gene in the pyramid lines suggests its central role in mediating immune responses in plants co-expressing the two transgenes, with reactive oxygen species (ROS) production enhanced by AdEXLB8 cues leading to stronger defense response. Here, we demonstrate that the association of potential priming elicitors and truncated NLRs can produce a synergistic effect on fungal resistance, constituting a promising strategy for improved, non-specific resistance to plant pathogens.

Transcriptome Responses of Wild Arachis to UV-C Exposure Reveal Genes Involved in General Plant Defense and Priming

Stress priming is an important strategy for enhancing plant defense capacity to deal with environmental challenges and involves reprogrammed transcriptional responses. Although ultraviolet (UV) light exposure is a widely adopted approach to elicit stress memory and tolerance in plants, the molecular mechanisms underlying UV-mediated plant priming tolerance are not fully understood. Here, we investigated the changes in the global transcriptome profile of wild Arachis stenosperma leaves in response to UV-C exposure. A total of 5751 differentially expressed genes (DEGs) were identified, with the majority associated with cell signaling, protein dynamics, hormonal and transcriptional regulation, and secondary metabolic pathways. The expression profiles of DEGs known as indicators of priming state, such as transcription factors, transcriptional regulators and protein kinases, were further characterized. A meta-analysis, followed by qRT-PCR validation, identified 18 metaDEGs as being commonly regulated in response to UV and other primary stresses. These genes are involved in secondary metabolism, basal immunity, cell wall structure and integrity, and may constitute important players in the general defense processes and establishment of a priming state in A. stenosperma. Our findings contribute to a better understanding of transcriptional dynamics involved in wild Arachis adaptation to stressful conditions of their natural habitats.

Ectopic expression of an expansin-like B gene from wild Arachis enhances tolerance to both abiotic and biotic stresses

Plant expansins are structural cell wall-loosening proteins implicated in several developmental processes and responses to environmental constraints and pathogen infection. To date, there is limited information about the biological function of expansins-like B (EXLBs), one of the smallest and less-studied subfamilies of plant expansins. In the present study, we conducted a functional analysis of the wild Arachis AdEXLB8 gene in transgenic tobacco (Nicotiana tabacum) plants to clarify its putative role in mediating defense responses to abiotic and biotic stresses. First, its cell wall localization was confirmed in plants expressing an AdEXLB8:eGFP fusion protein, while nanomechanical assays indicated cell wall reorganization and reassembly due to AdEXLB8 overexpression without compromising the phenotype. We further demonstrated that AdEXLB8 increased tolerance not only to isolated abiotic (drought) and biotic (Sclerotinia sclerotiorum and Meloidogyne incognita) stresses but also to their combination. The jasmonate and abscisic acid signaling pathways were clearly favored in transgenic plants, showing an activated antioxidative defense system. In addition to modifications in the biomechanical properties of the cell wall, we propose that AdEXLB8 overexpression interferes with phytohormone dynamics leading to a defense primed state, which culminates in plant defense responses against isolated and combined abiotic and biotic stresses.

Defining the combined stress response in wild Arachis

Nematodes and drought are major constraints in tropical agriculture and often occur simultaneously. Plant responses to these stresses are complex and require crosstalk between biotic and abiotic signaling pathways. In this study, we explored the transcriptome data of wild Arachis species subjected to drought (A-metaDEG) and the root-knot nematode Meloidogyne arenaria (B-metaDEG) via meta-analysis, to identify core-stress responsive genes to each individual and concurrent stresses in these species. Transcriptome analysis of a nematode/drought bioassay (cross-stress) showed that the set of stress responsive DEGs to concurrent stress is distinct from those resulting from overlapping A- and B-metaDEGs, indicating a specialized and unique response to combined stresses in wild Arachis. Whilst individual biotic and abiotic stresses elicit hormone-responsive genes, most notably in the jasmonic and abscisic acid pathways, combined stresses seem to trigger mainly the ethylene hormone pathway. The overexpression of a cross-stress tolerance candidate gene identified here, an endochitinase-encoding gene (AsECHI) from Arachis stenosperma, reduced up to 30% of M. incognita infection and increased post-drought recovery in Arabidopsis plants submitted to both stresses. The elucidation of the network of cross-stress responsive genes in Arachis contributes to better understanding the complex regulation of biotic and abiotic responses in plants facilitating more adequate crop breeding for combined stress tolerance.

The Defense Response Involved in Sweetpotato Resistance to Root-Knot Nematode Meloidogyne incognita: Comparison of Root Transcriptomes of Resistant and Susceptible Sweetpotato Cultivars With Respect to Induced and Constitutive Defense Responses

Sweetpotato (Ipomoea batatas [L.] Lam) is an economically important, nutrient- and pigment-rich root vegetable used as both food and feed. Root-knot nematode (RKN), Meloidogyne incognita, causes major yield losses in sweetpotato and other crops worldwide. The identification of genes and mechanisms responsible for resistance to RKN will facilitate the development of RKN resistant cultivars not only in sweetpotato but also in other crops. In this study, we performed RNA-seq analysis of RKN resistant cultivars (RCs; Danjami, Pungwonmi and Juhwangmi) and susceptible cultivars (SCs; Dahomi, Shinhwangmi and Yulmi) of sweetpotato infected with M. incognita to examine the induced and constitutive defense response-related transcriptional changes. During induced defense, genes related to defense and secondary metabolites were induced in SCs, whereas those related to receptor protein kinase signaling and protein phosphorylation were induced in RCs. In the uninfected control, genes involved in proteolysis and biotic stimuli showed differential expression levels between RCs and SCs during constitutive defense. Additionally, genes related to redox regulation, lipid and cell wall metabolism, protease inhibitor and proteases were putatively identified as RKN defense-related genes. The root transcriptome of SCs was also analyzed under uninfected conditions, and several potential candidate genes were identified. Overall, our data provide key insights into the transcriptional changes in sweetpotato genes that occur during induced and constitutive defense responses against RKN infection.

Comprehensive analysis of coding sequence architecture features and gene expression in Arachis duranensis

Coding sequence (CDS) architecture affects gene expression levels in organisms. Codon optimization can increase the gene expression level. Therefore, understanding codon usage patterns has important implications for research on genetic engineering and exogenous gene expression. To date, the codon usage patterns of many model plants have been analyzed. However, the relationship between CDS architecture and gene expression in Arachis duranensis remains poorly understood. According to the results of genome sequencing, A. duranensis has many resistant genes that can be used to improve the cultivated peanut. In this study, bioinformatic approaches were used to estimate A. duranensis CDS architectures, including frequency of the optimal codon (Fop), polypeptide length and GC contents at the first (GC1), second (GC2) and third (GC3) codon positions. In addition, Arachis RNA-seq datasets were downloaded from PeanutBase. The relationships between gene expression and CDS architecture were assessed both under normal growth as well as nematode and drought stress conditions. A total of 26 codons with high frequency were identified, which preferentially ended with A or T in A. duranensis CDSs under the above-mentioned three conditions. A similar CDS architecture was found in differentially expressed genes (DEGs) under nematode and drought stresses. The GC1 content differed between DEGs and non-differentially expressed genes (NDEGs) under both drought and nematode stresses. The expression levels of DEGs were affected by different CDS architectures compared with NDEGs under drought stress. In addition, no correlation was found between differential gene expression and CDS architecture neither under nematode nor under drought stress. These results aid the understanding of gene expression in A. duranensis.

Molecular insights into the compatible and incompatible interactions between sugar beet and the beet cyst nematode

BACKGROUND

Sugar beet (Beta vulgaris subsp. vulgaris) is an economically important crop that provides nearly one third of the global sugar production. The beet cyst nematode (BCN), Heterodera schachtii, causes major yield losses in sugar beet and other crops worldwide. The most effective and economic approach to control this nematode is growing tolerant or resistant cultivars. To identify candidate genes involved in susceptibility and resistance, the transcriptome of sugar beet and BCN in compatible and incompatible interactions at two time points was studied using mRNA-seq.

RESULTS

In the susceptible cultivar, most defense-related genes were induced at 4 dai while suppressed at 10 dai but in the resistant cultivar Nemakill, induction of genes involved in the plant defense response was observed at both time points. In the compatible interaction, alterations in phytohormone-related genes were detected. The effect of exogenous application of Methyl Jasmonate and ET-generator ethephon on susceptible plants was therefore investigated and the results revealed significant reduction in plant susceptibility. Genes putatively involved in the resistance of Nemakill were identified, such as genes involved in phenylpropanoid pathway and genes encoding CYSTM domain-containing proteins, F-box proteins, chitinase, galactono-1,4-lactone dehydrogenase and CASP-like protein. Also, the transcriptome of the BCN was analyzed in infected root samples and several novel potential nematode effector genes were found.

CONCLUSIONS

Our data provides detailed insights into the plant and nematode transcriptional changes occurring during compatible and incompatible interactions between sugar beet and BCN. Many important genes playing potential roles in susceptibility or resistance of sugar beet against BCN, as well as some BCN effectors with a potential role as avr proteins were identified. In addition, our findings indicate the effective role of jasmonate and ethylene in enhancing sugar beet defense response against BCN. This research provides new molecular insights into the plant-nematode interactions that can be used to design novel management strategies against BCN.

Genome-wide identification and characterization of nonspecific lipid transfer protein (nsLTP) genes in Arachis duranensis

Nonspecific lipid transfer proteins (nsLTPs) play vital roles in lipid metabolism, cell apoptosis and biotic and abiotic stresses in plants. However, the distribution of nsLTPs in Arachis duranensis has not been fully characterized. In this study, we identified 64 nsLTP genes in A. duranensis (designated AdLTPs), which were classified into six subfamilies and randomly distributed along nine chromosomes. Tandem and segmental duplication events were detected in the evolution of AdLTPs. The K_s and ω values differed significantly between Types 1 and D subfamilies, and eight AdLTPs were under positive selection. The expression levels of AdLTPs were changed after salinity, PEG, low-temperature and ABA treatments. Three AdLTPs were associated with resistance to nematode infection, and DOF and WRI1 transcription factors may regulate the AdLTP response to nematode infection. Our results may provide valuable genomic information for the breeding of peanut cultivars that are resistant to biotic and abiotic stresses.

Evolutionarily conserved plant genes responsive to root-knot nematodes identified by comparative genomics

Root-knot nematodes (RKNs, genus Meloidogyne) affect a large number of crops causing severe yield losses worldwide, more specifically in tropical and sub-tropical regions. Several plant species display high resistance levels to Meloidogyne, but a general view of the plant immune molecular responses underlying resistance to RKNs is still lacking. Combining comparative genomics with differential gene expression analysis may allow the identification of widely conserved plant genes involved in RKN resistance. To identify genes that are evolutionary conserved across plant species, we used OrthoFinder to compared the predicted proteome of 22 plant species, including important crops, spanning 214 Myr of plant evolution. Overall, we identified 35,238 protein orthogroups, of which 6,132 were evolutionarily conserved and universal to all the 22 plant species (PLAnts Common Orthogroups-PLACO). To identify host genes responsive to RKN infection, we analyzed the RNA-seq transcriptome data from RKN-resistant genotypes of a peanut wild relative (Arachis stenosperma), coffee (Coffea arabica L.), soybean (Glycine max L.), and African rice (Oryza glaberrima Steud.) challenged by Meloidogyne spp. using EdgeR and DESeq tools, and we found 2,597 (O. glaberrima), 743 (C. arabica), 665 (A. stenosperma), and 653 (G. max) differentially expressed genes (DEGs) during the resistance response to the nematode. DEGs' classification into the previously characterized 35,238 protein orthogroups allowed identifying 17 orthogroups containing at least one DEG of each resistant Arachis, coffee, soybean, and rice genotype analyzed. Orthogroups contain 364 DEGs related to signaling, secondary metabolite production, cell wall-related functions, peptide transport, transcription regulation, and plant defense, thus revealing evolutionarily conserved RKN-responsive genes. Interestingly, the 17 DEGs-containing orthogroups (belonging to the PLACO) were also universal to the 22 plant species studied, suggesting that these core genes may be involved in ancestrally conserved immune responses triggered by RKN infection. The comparative genomic approach that we used here represents a promising predictive tool for the identification of other core plant defense-related genes of broad interest that are involved in different plant-pathogen interactions.

Genomics of Plant Disease Resistance in Legumes

The constant interactions between plants and pathogens in the environment and the resulting outcomes are of significant importance for agriculture and agricultural scientists. Disease resistance genes in plant cultivars can break down in the field due to the evolution of pathogens under high selection pressure. Thus, the protection of crop plants against pathogens is a continuous arms race. Like any other type of crop plant, legumes are susceptible to many pathogens. The dawn of the genomic era, in which high-throughput and cost-effective genomic tools have become available, has revolutionized our understanding of the complex interactions between legumes and pathogens. Genomic tools have enabled a global view of transcriptome changes during these interactions, from which several key players in both the resistant and susceptible interactions have been identified. This review summarizes some of the large-scale genomic studies that have clarified the host transcriptional changes during interactions between legumes and their plant pathogens while highlighting some of the molecular breeding tools that are available to introgress the traits into breeding programs. These studies provide valuable insights into the molecular basis of different levels of host defenses in resistant and susceptible interactions.

Comparative transcriptome analyses provide novel insights into the differential response of Pigeonpea (Cajanus cajan L.) and its wild relative (Cajanus platycarpus (Benth.) Maesen) to herbivory by Helicoverpa armigera (Hübner)

Deeper insights into the resistance response of Cajanus platycarpus were obtained based on comparative transcriptomics under Helicoverpa armigera infestation. Devastation by pod borer, Helicoverpa armigera is one of the major factors for stagnated productivity in Pigeonpea. Despite possessing a multitude of desirable traits including pod borer resistance, wild relatives of Cajanus spp. have remained under-utilized due to linkage drag and cross-incompatibility. Discovery and deployment of genes from them can provide means to tackle key pests like H. armigera. Transcriptomic differences between Cajanus platycarpus and Cajanus cajan during different time points (0, 18, 38, 96 h) of pod borer infestation were elucidated in this study. For the first ever time, we demonstrated captivating variations in their response; C. platycarpus apparently being reasonably agile with effectual transcriptomic reprogramming to deter the insect. Deeper insights into the differential response were obtained by identification of significant GO-terms related to herbivory followed by combined KEGG and ontology analyses. C. platycarpus portrayed a multilevel response with cardinal involvement of SAR, redox homeostasis and reconfiguration of primary metabolites leading to a comprehensive defense response. The credibility of RNA-seq analyses was ascertained by transient expression of selected putative insect resistance genes from C. platycarpus viz., chitinase (CHI4), Alpha-amylase/subtilisin inhibitor (IAAS) and Flavonoid 3_5 hydroxylase (C75A1) in Nicotiana benthamiana followed by efficacy analysis against H. armigera. qPCR validated results of the study provided innovative insights and useful leads for development of durable pod borer resistance.

Differential tetraspanin genes expression and subcellular localization during mutualistic interactions in Phaseolus vulgaris

Arbuscular mycorrhizal fungi and rhizobia association with plants are two of the most successful plant-microbe associations that allow the assimilation of P and N by plants, respectively. These mutualistic interactions require a molecular dialogue, i.e., legume roots exude flavonoids or strigolactones which induce the Nod factors or Myc factors synthesis and secretion from the rhizobia or fungi, respectively. These Nod or Myc factors trigger several responses in the plant root, including calcium oscillations, and reactive oxygen species (ROS). Furthermore, superoxide and H₂O₂ have emerged as key components that regulate the transitions from proliferation to differentiation in the plant meristems. Similar to the root meristem, the nodule meristem accumulates superoxide and H₂O₂. Tetraspanins are transmembrane proteins that organize into tetraspanin web regions, where they recruit specific proteins into platforms required for signal transduction, membrane fusion, cell trafficking and ROS generation. Plant tetraspanins are scaffolding proteins associated with root radial patterning, biotic and abiotic stress responses, cell fate determination, and hormonal regulation and recently have been reported as a specific marker of exosomes in animal and plant cells and key players at the site of plant fungal infection. In this study, we conducted transcriptional profiling of the tetraspanin family in common bean (Phaseolus vulgaris L. var. Negro Jamapa) to determine the specific expression patterns and subcellular localization of tetraspanins during nodulation or under mycorrhizal association. Our results demonstrate that the tetraspanins are transcriptionally modulated during the mycorrhizal association, but are also expressed in the infection thread and nodule meristem development. Subcellular localization indicates that tetraspanins have a key role in vesicular trafficking, cell division, and root hair polar growth.

The characteristic of Arachis duranensis-specific genes and their potential function

Arachis species produce flowers aerially, and then grow into the ground, where they develop into fruits; a feature that is unique to Arachis species. We hypothesized that Arachis species evolved genes specifically involved in the control of aerial flowers and the formation of underground fruits. Arachis duranensis is more resistant to biotic and abiotic stressors. Here, we compared different legume species and identified Arachis duranensis-specific genes. We analyzed gene expression patterns, base substitution patterns and sequence features between genes that are conserved across legume plants and A. duranensis-specific genes. Furthermore, we tested the role of A. duranensis-specific genes during seed development, response to nematode Meloidogyne arenaria infection and drought stress. We found that A. duranensis-specific genes had characteristics of young genes. The gene expression level and breadth were lower in the A. duranensis-specific genes compared to conserved genes. The A. duranensis-specific genes had higher codon usage bias than conserved genes, and the polypeptide length and GC content at the three codon sites were lower compared to conserved genes. Of the A. duranensis-specific genes, single-copy and duplicated genes had different features. The RNA-seq result showed A. duranensis-specific genes were involved in seed development, as well as response to nematode infection and drought stress. In addition, we detected asymmetric functions in A. duranensis-specific duplicated genes in response to nematode infection and drought stress.

Comparative transcriptome analysis of pigeonpea, Cajanus cajan (L.) and one of its wild relatives Cajanus platycarpus (Benth.) Maesen

Pigeonpea is a major source of dietary protein to the vegetarian population of the Indian sub-continent. Crop improvement to mitigate biotic and abiotic stresses for realization of its potential yield and bridging yield gap is the need of the hour. Availability of limited genomic resources in the cultivated germplasm, however, is a serious bottleneck towards successful molecular breeding for the development of superior genotypes in pigeonpea. In view of this, improvement of pigeonpea can be attempted through transgenesis or by exploiting genetic resources from its wild relatives. Pigeonpea wild relatives are known to be bestowed with agronomic traits of importance; discovery and deployment of genes from them can provide a lucrative option for crop improvement. Understanding molecular signatures of wild relatives would not only provide information about the mechanism behind desired traits but also enable us to extrapolate the information to cultivated pigeonpea. The present study deals with the characterization of leaf transcriptomes of Cajanus cajan and one of its wild relatives, Cajanus platycarpus. Illumina sequencing revealed 0.11 million transcripts in both the species with an annotation of 0.09 million (82%) transcripts using BLASTX. Comparative transcriptome analyses on the whole, divulged cues about the wild relative being vigilant and agile. Gene ontology and Mapman analysis depicted higher number of transcripts in the wild relative pertaining to signaling, transcription factors and stress responsive genes. Further, networking between the differentially expressed MapMan bins demonstrated conspicuous interactions between different bins through 535 nodes (512 Genes and 23 Pathways) and 1857 edges. The authenticity of RNA-seq analysis was confirmed by qRT-PCR. The information emanating from this study can provide valuable information and resource for future translational research including genome editing to alleviate varied stresses. Further, this learning can be a platform for in-depth investigations to decipher molecular mechanisms for mitigation of various stresses in the wild relative.

Transcriptomic changes in sweetpotato peroxidases in response to infection with the root-knot nematode Meloidogyne incognita

A previous transcriptomic analysis of the roots of susceptible and resistant cultivars of sweetpotato (Ipomoea batatas) identified genes that were likely to contribute to protection against infection with the root-knot nematode Meloidogyne incognita. The current study examined the roles of peroxidase genes in sweetpotato defense responses during root-knot nematode infection, using the susceptible (cv. Yulmi) and resistant (cv. Juhwangmi) cultivars. Differentially expressed genes were assigned to gene ontology categories to predict their functional roles and associated biological processes. Comparison with Arabidopsis peroxidases identified a group of genes orthologous to Arabidopsis PEROXIDASE 52 (AtPrx52). An analysis of sweetpotato peroxidase genes determined their roles in protecting plants against root-knot nematode infection and enabled identification of important peroxidases. The interactions involved in sweetpotato resistance to nematode infection are discussed.

Old and young duplicate genes reveal different responses to environmental changes in Arachis duranensis

Old and young duplicate genes have been reported in some organisms. However, little is known about the properties of old and young duplicate genes in Arachis. Here, we have identified old and young duplicate genes in Arachis duranensis, and analyzed the evolution, gene complexity, gene expression pattern, and functional divergence between old and young duplicate genes. Our results showed different evolutionary, gene complexity and gene expression patterns, as well as differing correlations between old and young duplicate genes. Gene ontology results showed that old duplicate genes play a crucial role in lipid and amino acid biosynthesis and the oxidation-reduction process and that young duplicate genes are preferentially involved in photosynthesis and response to biotic stimulus. Transcriptome data sets revealed that most old and young duplicate genes had asymmetric function, and only a few duplicate genes exhibited symmetric function under drought and nematode stress. We found that old duplicate genes are preferentially involved in lipid and amino acid metabolism and response to abiotic stress, while young duplicate genes are likely to participate in photosynthesis and response to biotic stress. This work provides a better understanding of the evolution and functional divergence of old and young duplicate genes in A. duranensis.

Contrasting Effects of Wild Arachis Dehydrin Under Abiotic and Biotic Stresses

Plant dehydrins (DNHs) belong to the LEA (Late Embryogenesis Abundant) protein family and are involved in responses to multiple abiotic stresses. DHNs are classified into five subclasses according to the organization of three conserved motifs (K-; Y-; and S-segments). In the present study, the DHN protein family was characterized by molecular phylogeny, exon/intron organization, protein structure, and tissue-specificity expression in eight Fabaceae species. We identified 20 DHN genes, encompassing three (YnSKn, SKn, and Kn) subclasses sharing similar gene organization and protein structure. Two additional low conserved DHN Φ-segments specific to the legume SKn-type of proteins were also found. The in silico expression patterns of DHN genes in four legume species (Arachis duranensis, A. ipaënsis, Glycine max, and Medicago truncatula) revealed that their tissue-specific regulation is associated with the presence or absence of the Y-segment. Indeed, DHN genes containing a Y-segment are mainly expressed in seeds, whereas those without the Y-segment are ubiquitously expressed. Further qRT-PCR analysis revealed that, amongst stress responsive dehydrins, a SKn-type DHN gene from A. duranensis (AdDHN1) showed opposite response to biotic and abiotic stress with a positive regulation under water deficit and negative regulation upon nematode infection. Furthermore, transgenic Arabidopsis lines overexpressing (OE) AdDHN1 displayed improved tolerance to multiple abiotic stresses (freezing and drought) but increased susceptibility to the biotrophic root-knot nematode (RKN) Meloidogyne incognita. This contradictory role of AdDHN1 in responses to abiotic and biotic stresses was further investigated by qRT-PCR analysis of transgenic plants using a set of stress-responsive genes involved in the abscisic acid (ABA) and jasmonic acid (JA) signaling pathways and suggested an involvement of DHN overexpression in these stress-signaling pathways.

Genome-wide identification of LRR-containing sequences and the response of these sequences to nematode infection in Arachis duranensis

BACKGROUND

Leucine-rich repeat (LRR)-containing genes are involved in responses to various diseases. Recently, RNA-seq data from A. duranensis after nematode (Meloidogyne arenaria) infection were released. However, the number of LRR-containing genes present in A. duranensis and the response of LRR-containing genes to nematode infection are poorly understood.

RESULTS

In this study, we found 509 amino acid sequences containing nine types of LRR domains in A. duranensis. The inferred phylogenetic relationships revealed that the nine types of LRR domains had two originations. The inferred selective pressure was mainly consistent with LRR domains undergoing purifying selection. Twenty-one LRR-containing genes were associated with possible resistance to nematode infection after 3, 6, and 9 days. Among them, Aradu.T5WNW, Aradu.JM17V, and Aradu.MKP1A were up-regulate at these three time points, while Aradu.QD5DS and Aradu.M0ENQ were up-regulated 6 and 9 days after nematode infection. The expression of the above mentioned five genes was significantly and negatively correlated with the number of LRR8 domain, indicating that fewer LRR8 domains are associated with the promotion of LRR-containing genes that resist nematode infection. Patterns of co-expression and cis-acting elements indicated that WRKY possibly regulate the responses of LRR-containing genes to nematode infection and that expansin genes may work together with LRR-containing genes in response to nematode infection.

CONCLUSIONS

We identified the number and type of LRR-containing genes in A. duranensis. The LRR-containing genes that were found appear to be involved in responses to nematode infection. The number of LRR8 domains was negatively correlated with expression after nematode infection. The WRKY transcription factor may regulate resistance to nematode infection based on LRR-containing genes. Our results could improve the understanding of resistance to nematodes and molecular breeding in peanuts.

Comparative analysis of selection mode reveals different evolutionary rate and expression pattern in Arachis duranensis and Arachis ipaënsis duplicated genes

Our results reveal that K_s is a determining factor affecting selective pressure and different evolution and expression patterns are detected between PSGs and NSGs in wild Arachis duplicates. Selective pressure, including purifying (negative) and positive selection, can be detected in organisms. However, studies on comparative evolutionary rates, gene expression patterns and gene features between negatively selected genes (NSGs) and positively selected genes (PSGs) are lagging in paralogs of plants. Arachis duranensis and Arachis ipaënsis are ancestors of the cultivated peanut, an important oil and protein crop. Here, we carried out a series of systematic analyses, comparing NSG and PSG in paralogs, using genome sequences and transcriptome datasets in A. duranensis and A. ipaënsis. We found that synonymous substitution rate (K_s) is a determining factor affecting selective pressure in A. duranensis and A. ipaënsis duplicated genes. Lower expression level, lower gene expression breadth, higher codon bias and shorter polypeptide length were found in PSGs and not in NSGs. The correlation analyses showed that gene expression breadth was positively correlated with polypeptide length and GC content at the first codon site (GC1) in PSGs and NSGs, respectively. There was a negative correlation between expression level and polypeptide length in PSGs. In NSGs, the K_s was positively correlated with expression level, gene expression breadth, GC1, and GC content at the third codon site (GC3), but selective pressure was negatively correlated with expression level, gene expression breadth, polypeptide length, GC1, and GC3 content. The function of most duplicated gene pairs was divergent under drought and nematode stress. Taken together, our results show that different evolution and expression patterns occur between PSGs and NSGs in paralogs of two wild Arachis species.

Comparative Transcriptome Profiling Reveals Defense-Related Genes against Meloidogyne incognita Invasion in Tobacco

Root-knot nematodes Meloidogyne incognita are one of the most destructive pathogens, causing severe losses to tobacco productivity and quality. However, the underlying resistance mechanism of tobacco to M. incognita is not clear. In this study, two tobacco genotypes, K326 and Changbohuang, which are resistant and susceptible to M. incognita, respectively, were used for RNA-sequencing analysis. An average of 35 million clean reads were obtained. Compared with their expression levels in non-infected plants of the same genotype, 4354 and 545 differentially expressed genes (DEGs) were detected in the resistant and susceptible genotype, respectively, after M. incognita invasion. Overall, 291 DEGs, involved in diverse biological processes, were common between the two genotypes. Genes encoding toxic compound synthesis, cell wall modification, reactive oxygen species and the oxidative burst, salicylic acid signal transduction, and production of some other metabolites were putatively associated with tobacco resistance to M. incognita. In particular, the complex resistance response needed to overcome M. incognita invasion may be regulated by several transcription factors, such as the ethylene response factor, MYB, basic helix⁻loop⁻helix transcription factor, and indole acetic acid⁻leucine-resistant transcription factor. These results may aid in the identification of potential genes of resistance to M. incognita for tobacco cultivar improvement.

Comparative root transcriptome of wild Arachis reveals NBS-LRR genes related to nematode resistance

BACKGROUND

The Root-Knot Nematode (RKN), Meloidogyne arenaria, significantly reduces peanut grain quality and yield worldwide. Whilst the cultivated species has low levels of resistance to RKN and other pests and diseases, peanut wild relatives (Arachis spp.) show rich genetic diversity and harbor high levels of resistance to many pathogens and environmental constraints. Comparative transcriptome analysis can be applied to identify candidate resistance genes.

RESULTS

Transcriptome analysis during the early stages of RKN infection of two peanut wild relatives, the highly RKN resistant Arachis stenosperma and the moderately susceptible A. duranensis, revealed genes related to plant immunity with contrasting expression profiles. These included genes involved in hormone signaling and secondary metabolites production and also members of the NBS-LRR class of plant disease resistance (R) genes. From 345 NBS-LRRs identified in A.duranensis reference genome, 52 were differentially expressed between inoculated and control samples, with the majority occurring in physical clusters unevenly distributed on eight chromosomes with preferential tandem duplication. The majority of these NBS-LRR genes showed contrasting expression behaviour between A. duranensis and A. stenosperma, particularly at 6 days after nematode inoculation, coinciding with the onset of the Hypersensitive Response in the resistant species. The physical clustering of some of these NBS-LRR genes correlated with their expression patterns in the contrasting genotypes. Four NBS-LRR genes exclusively expressed in A. stenosperma are located within clusters on chromosome Aradu. A09, which harbors a QTL for RKN resistance, suggesting a functional role for their physical arrangement and their potential involvement in this defense response.

CONCLUSION

The identification of functional novel R genes in wild Arachis species responsible for triggering effective defense cascades can contribute to the crop genetic improvement and enhance peanut resilience to RKN.

Early responses to dehydration in contrasting wild Arachis species

Wild peanut relatives (Arachis spp.) are genetically diverse and were selected throughout evolution to a range of environments constituting, therefore, an important source of allelic diversity for abiotic stress tolerance. In particular, A. duranensis and A. stenosperma, the parents of the reference Arachis A-genome genetic map, show contrasting transpiration behavior under limited water conditions. This study aimed to build a comprehensive gene expression profile of these two wild species under dehydration stress caused by the withdrawal of hydroponic nutrient solution. For this purpose, roots of both genotypes were collected at seven time-points during the early stages of dehydration and used to construct cDNA paired-end libraries. Physiological analyses indicated initial differences in gas exchange parameters between the drought-tolerant genotype of A. duranensis and the drought-sensitive genotype of A. stenosperma. High-quality Illumina reads were mapped against the A. duranensis reference genome and resulted in the identification of 1,235 and 799 Differentially Expressed Genes (DEGs) that responded to the stress treatment in roots of A. duranensis and A. stenosperma, respectively. Further analysis, including functional annotation and identification of biological pathways represented by these DEGs confirmed the distinct gene expression behavior of the two contrasting Arachis species genotypes under dehydration stress. Some species-exclusive and common DEGs were then selected for qRT-PCR analysis, which corroborated the in silico expression profiling. These included genes coding for regulators and effectors involved in drought tolerance responses, such as activation of osmosensing molecular cascades, control of hormone and osmolyte content, and protection of macromolecules. This dataset of transcripts induced during the dehydration process in two wild Arachis genotypes constitute new tools for the understanding of the distinct gene regulation processes in these closely related species but with contrasting drought responsiveness. In addition, our findings provide insights into the nature of drought tolerance in wild germoplasm, which might be explored as novel sources of diversity and useful wild alleles to develop climate-resilient crop varieties.

Genome-wide analysis of expansin superfamily in wild Arachis discloses a stress-responsive expansin-like B gene

Expansins are plant cell wall-loosening proteins involved in adaptive responses to environmental stimuli and various developmental processes. The first genome-wide analysis of the expansin superfamily in the Arachis genus identified 40 members in A. duranensis and 44 in A. ipaënsis, the wild progenitors of cultivated peanut (A. hypogaea). These expansins were further characterized regarding their subfamily classification, distribution along the genomes, duplication events, molecular structure, and phylogeny. A RNA-seq expression analysis in different Arachis species showed that the majority of these expansins are modulated in response to diverse stresses such as water deficit, root-knot nematode (RKN) infection, and UV exposure, with an expansin-like B gene (AraEXLB8) displaying a highly distinct stress-responsive expression profile. Further analysis of the AraEXLB8 coding sequences showed high conservation across the Arachis genotypes, with eight haplotypes identified. The modulation of AraEXLB8 expression in response to the aforementioned stresses was confirmed by qRT-PCR analysis in distinct Arachis genotypes, whilst in situ hybridization revealed transcripts in different root tissues according to the stress imposed. The overexpression of AraEXLB8 in soybean (Glycine max) composite plants remarkably decreased the number of galls in transformed hairy roots inoculated with RKN. This study improves the current understanding of the molecular evolution, divergence, and gene expression of expansins in Arachis, and provides molecular and functional insights into the role of expansin-like B, the less-studied plant expansin subfamily.

Ex vitro hairy root induction in detached peanut leaves for plant-nematode interaction studies

BACKGROUND

Peanut (Arachis hypogaea) production is largely affected by a variety of abiotic and biotic stresses, including the root-knot nematode (RKN) Meloidogyne arenaria that causes yield losses worldwide. Transcriptome studies of wild Arachis species, which harbor resistance to a number of pests and diseases, disclosed several candidate genes for M. arenaria resistance. Peanut is recalcitrant to genetic transformation, so the use of Agrobacterium rhizogenes-derived hairy roots emerged as an alternative for in-root functional characterization of these candidate genes.

RESULTS

The present report describes an ex vitro methodology for hairy root induction in detached leaves based on the well-known ability of peanut to produce roots spontaneously from its petiole, which can be maintained for extended periods under high-humidity conditions. Thirty days after infection with the A. rhizogenes 'K599' strain, 90% of the detached leaves developed transgenic hairy roots with 5 cm of length in average, which were then inoculated with M. arenaria. For improved results, plant transformation, and nematode inoculation parameters were adjusted, such as bacterial cell density and growth stage; moist chamber conditions and nematode inoculum concentration. Using this methodology, a candidate gene for nematode resistance, AdEXLB8, was successfully overexpressed in hairy roots of the nematode-susceptible peanut cultivar 'Runner', resulting in 98% reduction in the number of galls and egg masses compared to the control, 60 days after M. arenaria infection.

CONCLUSIONS

This methodology proved to be more practical and cost-effective for functional validation of peanut candidate genes than in vitro and composite plant approaches, as it requires less space, reduces analysis costs and displays high transformation efficiency. The reduction in the number of RKN galls and egg masses in peanut hairy roots overexpressing AdEXLB8 corroborated the use of this strategy for functional characterization of root expressing candidate genes. This approach could be applicable not only for peanut-nematode interaction studies but also to other peanut root diseases, such as those caused by fungi and bacteria, being also potentially extended to other crop species displaying similar petiole-rooting competence.

Transcriptomic and histological responses of African rice (Oryza glaberrima) to Meloidogyne graminicola provide new insights into root-knot nematode resistance in monocots

BACKGROUND AND AIMS

The root-knot nematode Meloidogyne graminicola is responsible for production losses in rice ( Oryza sativa ) in Asia and Latin America. The accession TOG5681 of African rice, O. glaberrima , presents improved resistance to several biotic and abiotic factors, including nematodes. The aim of this study was to assess the cytological and molecular mechanisms underlying nematode resistance in this accession.

METHODS

Penetration and development in M. graminicola in TOG5681 and the susceptible O. sativa genotype 'Nipponbare' were compared by microscopic observation of infected roots and histological analysis of galls. In parallel, host molecular responses to M. graminicola were assessed by root transcriptome profiling at 2, 4 and 8 d post-infection (dpi). Specific treatments with hormone inhibitors were conducted in TOG5681 to assess the impact of the jasmonic acid and salicylic acid pathways on nematode penetration and reproduction.

KEY RESULTS

Penetration and development of M. graminicola juveniles were reduced in the resistant TOG5681 in comparison with the susceptible accession, with degeneration of giant cells observed in the resistant genotype from 15 dpi onwards. Transcriptome changes were observed as early as 2 dpi, with genes predicted to be involved in defence responses, phenylpropanoid and hormone pathways strongly induced in TOG5681, in contrast to 'Nipponbare'. No specific hormonal pathway could be identified as the major determinant of resistance in the rice-nematode incompatible interaction. Candidate genes proposed as involved in resistance to M. graminicola in TOG5681 were identified based on their expression pattern and quantitative trait locus (QTL) position, including chalcone synthase, isoflavone reductase, phenylalanine ammonia lyase, WRKY62 transcription factor, thionin, stripe rust resistance protein, thaumatins and ATPase3.

CONCLUSIONS

This study provides a novel set of candidate genes for O. glaberrima resistance to nematodes and highlights the rice- M. graminicola pathosystem as a model to study plant-nematode incompatible interactions.

Gene expression analysis in Musa acuminata during compatible interactions with Meloidogyne incognita

BACKGROUND AND AIMS

Endoparasitic root-knot nematodes (RKNs) ( Meloidogyne spp.) cause considerable losses in banana ( Musa spp.), with Meloidogyne incognita a predominant species in Cavendish sub-group bananas. This study investigates the root transcriptome in Musa acuminata genotypes 4297-06 (AA) and Cavendish Grande Naine (CAV; AAA) during early compatible interactions with M. incognita .

METHODS

Roots were analysed by brightfield light microscopy over a 35 d period to examine nematode penetration and morphological cell transformation. RNA samples were extracted 3, 7 and 10 days after inoculation (DAI) with nematode J2 juveniles, and cDNA libraries were sequenced using lllumina HiSeq technology. Sequences were mapped to the M. acuminata ssp. malaccensis var. Pahang genome sequence, differentially expressed genes (DEGs) identified and transcript representation determined by gene set enrichment and pathway mapping.

KEY RESULTS

Microscopic analysis revealed a life cycle of M. incognita completing in 24 d in CAV and 27 d in 4279-06. Comparable numbers of DEGs were up- and downregulated in each genotype, with potential involvement of many in early host defence responses involving reactive oxygen species and jasmonate/ethylene signalling. DEGs revealed concomitant auxin metabolism and cell wall modification processes likely to be involved in giant cell formation. Notable transcripts related to host defence included those coding for leucine-rich repeat receptor-like serine/threonine-protein kinases, peroxidases, thaumatin-like pathogenesis-related proteins, and DREB, ERF, MYB, NAC and WRKY transcription factors. Transcripts related to giant cell development included indole acetic acid-amido synthetase GH3.8 genes, involved in auxin metabolism, as well as genes encoding expansins and hydrolases, involved in cell wall modification.

CONCLUSIONS

Expression analysis in M. acuminata during compatible interactions with RKNs provides insights into genes modulated during infection and giant cell formation. Increased understanding of both defence responses to limit parasitism during compatible interactions and effector-targeted host genes in this complex interaction will facilitate the development of genetic improvement measures for RKNs.

Transcriptome analysis of resistant and susceptible tobacco (Nicotiana tabacum) in response to root-knot nematode Meloidogyne incognita infection

The root-knot nematode (RKN) Meloidogyne incognita reproduces on the roots of tobacco (Nicotiana tabacum), damaging crops, reducing crop yield, and causing economic losses annually. The development of resistant genotypes is an alternative strategy to effectively control these losses. However, the molecular mechanism responsible for host pathogenesis and defense responses in tobacco specifically against RKNs remain poorly understood. Here, root transcriptome analysis of resistant (Yuyan12) and susceptible (Changbohuang) tobacco varieties infected with RKNs was performed. Moreover, 2623 and 545 differentially expressed genes (DEGs) in RKN-infected roots were observed in Yuyan12 and Changbohuang, respectively, compared to those in non-infected roots, including 289 DEGs commonly expressed in the two genotypes. Among these DEGs, genes encoding cell wall modifying proteins, auxin-related proteins, the ROS scavenging system, and transcription factors involved in various biological and physiochemical processes were significantly expressed in both the resistant and susceptible genotypes. This work is thus the first report on the relationships in the RKN-tobacco interaction using transcriptome analysis, and the results provide important information on the mechanism of RKN resistance in tobacco.