PRGdb 3.0: a comprehensive platform for prediction and analysis of plant disease resistance genes

Genome architecture of the allotetraploid wild grass Aegilops ventricosa reveals its evolutionary history and contributions to wheat improvement

The allotetraploid wild grass Aegilops ventricosa (2n = 4x = 28, genome DvDvNvNv) has been recognized as an important germplasm resource for wheat improvement owing to its ability to tolerate biotic stresses. In particular, the 2NvS segment from Ae. ventricosa, as a stable and effective resistance source, has contributed greatly to wheat improvement. The 2NvS/2AS translocation is a prevalent chromosomal translocation between common wheat and wild relatives, ranking just behind the 1B/1R translocation in importance for modern wheat breeding. Here, we assembled a high-quality chromosome-level reference genome of Ae. ventricosa RM271 with a total length of 8.67 Gb. Phylogenomic analyses revealed that the progenitor of the Dv subgenome of Ae. ventricosa is Ae. tauschii ssp. tauschii (genome DD); by contrast, the progenitor of the D subgenome of bread wheat (Triticum aestivum L.) is Ae. tauschii ssp. strangulata (genome DD). The oldest polyploidization time of Ae. ventricosa occurred ∼0.7 mya. The Dv subgenome of Ae. ventricosa is less conserved than the D subgenome of bread wheat. Construction of a graph-based pangenome of 2AS/6NvL (originally known as 2NvS) segments from Ae. ventricosa and other genomes in the Triticeae enabled us to identify candidate resistance genes sourced from Ae. ventricosa. We identified 12 nonredundant introgressed segments from the Dv and Nv subgenomes using a large winter wheat collection representing the full diversity of the European wheat genetic pool, and 29.40% of European wheat varieties inherit at least one of these segments. The high-quality RM271 reference genome will provide a basis for cloning key genes, including the Yr17-Lr37-Sr38-Cre5 resistance gene cluster in Ae. ventricosa, and facilitate the full use of elite wild genetic resources to accelerate wheat improvement.

Telomere-to-telomere Citrullus super-pangenome provides direction for watermelon breeding

To decipher the genetic diversity within the cucurbit genus Citrullus, we generated telomere-to-telomere (T2T) assemblies of 27 distinct genotypes, encompassing all seven Citrullus species. This T2T super-pangenome has expanded the previously published reference genome, T2T-G42, by adding 399.2 Mb and 11,225 genes. Comparative analysis has unveiled gene variants and structural variations (SVs), shedding light on watermelon evolution and domestication processes that enhanced attributes such as bitterness and sugar content while compromising disease resistance. Multidisease-resistant loci from Citrullus amarus and Citrullus mucosospermus were successfully introduced into cultivated Citrullus lanatus. The SVs identified in C. lanatus have not only been inherited from cordophanus but also from C. mucosospermus, suggesting additional ancestors beyond cordophanus in the lineage of cultivated watermelon. Our investigation substantially improves the comprehension of watermelon genome diversity, furnishing comprehensive reference genomes for all Citrullus species. This advancement aids in the exploration and genetic enhancement of watermelon using its wild relatives.

Melanaphis sacchari/sorghi complex: current status, challenges and integrated strategies for managing the invasive sap-feeding insect pest of sorghum

Melanaphis sacchari (Zehntner;Hemiptera: Aphididae), sugarcane aphid (SCA), is an invasive phloem-feeder found worldwide with a wide host range of economically important plants including sorghum and sugarcane. Given its high reproductive capacity and ability to rapidly spread over long distances, SCA presents challenges for effective control, leading to substantial economic losses. Recent studies have identified two multiloci SCA genotypes specialized in feeding on sugarcane (MLL-D) and sorghum (MLL-F) in the USA, which raises concerns as the USA is the second largest sorghum-producing country. This has encouraged research towards identifying these two biotypes where some research has stated them as two species; MLL-D clade to be M. sacchari and MLL-F clade to be M. sorghi Theobald (Hemiptera: Aphididae), sorghum aphid (SA). This review aims at compiling research progress that has been made on understanding the SCA/SA species complex. Furthermore, this review also highlights a wide range of management strategies against SCA/SA that includes both biological and chemical methods. In addition, the review emphasizes studies examining host plant resistance to understand and evaluate the role of R-genes and phytohormones such as jasmonic acid, salicylic acid and ethylene against SCA. Beside this, plant volatiles and other secondary metabolites such as flavonoids, terpenes and phytanes are also explored as potential control agents. Being an invasive pest, a single management tactic is inadequate to control SCA population and hence, integrated pest management practices incorporating physical, cultural and biological control methods should be implemented with exclusive chemical control as a last resort, which this review examines in detail. © 2024 The Author(s). Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Transcriptomics and metabolomics analyses provide insights into resistance genes of tree ferns

As ancient organisms, tree ferns play a crucial role as an evolutionary bridge between lower and higher plant species, providing various utilitarian benefits. However, they face challenges such as overexploitation, climate change, adverse environmental conditions, and insect pests, resulting in conservation concerns. In this study, we provide an overview of metabolic and transcriptomic resources of leaves in two typical tree ferns, A. spinulosa and A. metteniana, and explore the resistance genes for the first time. The landscape of metabolome showed that the compound skimmin may hold medicinal significance. A total of 111 differentially accumulated metabolites (DAMs) were detected, with pathway enrichment analysis highlighting 14 significantly enriched pathways, including 2-oxocarboxylic acid metabolism possibly associated with environmental adaptations. A total of 14,639 differentially expressed genes (DEGs) were found, among which 606 were resistance (R) genes. We identified BAM1 as a significantly differentially expressed R gene, which is one of the core genes within the R gene interaction network. Both the maximum-likelihood phylogenetic tree and the PPI network revealed a close relationship between BAM1, FLS2, and TMK. Moreover, BAM1 showed a significant positive correlation with neochlorogenic acid and kaempferol-7-O-glucoside. These metabolites, known for their antioxidant and anti-inflammatory properties, likely play a crucial role in the defense response of tree ferns. This research provides valuable insights into the metabolic and transcriptomic differences between A. spinulosa and A. metteniana, enhancing our understanding of resistance genes in tree ferns.

A haplotype-resolved chromosome-level assembly and annotation of European hazelnut (C. avellana cv. Jefferson) provides insight into mechanisms of eastern filbert blight resistance

European hazelnut (Corylus avellana L.) is an important tree nut crop. Hazelnut production in North America is currently limited in scalability due to Anisogramma anomala, a fungal pathogen that causes Eastern Filbert Blight (EFB) disease in hazelnut. Successful deployment of EFB resistant cultivars has been limited to the state of Oregon, where the breeding program at Oregon State University (OSU) has released cultivars with a dominant allele at a single resistance locus identified by classical breeding, linkage mapping, and molecular markers. C. avellana cultivar "Jefferson" is resistant to the predominant EFB biotype in Oregon and has been selected by the OSU breeding program as a model for hazelnut genetic and genomic research. Here, we present a near complete, haplotype-resolved chromosome-level hazelnut genome assembly for "Jefferson". This new assembly is a significant improvement over a previously published genome draft. Analysis of genomic regions linked to EFB resistance and self-incompatibility confirmed haplotype splitting and identified new gene candidates that are essential for downstream molecular marker development, thereby facilitating breeding efforts.

Transcriptome Analysis Reveals Key Genes Involved in Trichome Formation in Pepper (Capsicum annuum L.)

Trichomes are specialized organs located in the plant epidermis that play important defense roles against biotic and abiotic stresses. However, the mechanisms regulating the development of pepper epidermal trichomes and the related regulatory genes at the molecular level are not clear. Therefore, we performed transcriptome analyses of A114 (less trichome) and A115 (more trichome) to dig deeper into the genes involved in the regulatory mechanisms of epidermal trichome development in peppers. In this study, the epidermal trichome density of A115 was found to be higher by phenotypic observation and was highest in the leaves at the flowering stage. A total of 39,261 genes were quantified by RNA-Seq, including 11,939 genes not annotated in the previous genome analysis and 18,833 differentially expressed genes. Based on KEGG functional enrichment, it was found that DEGs were mainly concentrated in three pathways: plant-pathogen interaction, MAPK signaling pathway-plant, and plant hormone signal transduction. We further screened the DEGs associated with the development of epidermal trichomes in peppers, and the expression of the plant signaling genes GID1B-like (Capana03g003488) and PR-6 (Capana09g001847), the transcription factors MYB108 (Capana05g002225) and ABR1-like (Capana04g001261), and the plant resistance genes PGIP-like (Capana09g002077) and At5g49770 (Capana08g001721) in the DEGs were higher at A115 compared to A114, and were highly expressed in leaves at the flowering stage. In addition, based on the WGCNA results and the establishment of co-expression networks showed that the above genes were highly positively correlated with each other. The transcriptomic data and analysis of this study provide a basis for the study of the regulatory mechanisms of pepper epidermal trichomes.

Genome-wide analysis of the peanut CaM/CML gene family reveals that the AhCML69 gene is associated with resistance to Ralstonia solanacearum

BACKGROUND

Calmodulins (CaMs)/CaM-like proteins (CMLs) are crucial Ca2+-binding sensors that can decode and transduce Ca2+ signals during plant development and in response to various stimuli. The CaM/CML gene family has been characterized in many plant species, but this family has not yet been characterized and analyzed in peanut, especially for its functions in response to Ralstonia solanacearum. In this study, we performed a genome-wide analysis to analyze the CaM/CML genes and their functions in resistance to R. solanacearum.

RESULTS

Here, 67, 72, and 214 CaM/CML genes were identified from Arachis duranensis, Arachis ipaensis, and Arachis hypogaea, respectively. The genes were divided into nine subgroups (Groups I-IX) with relatively conserved exon‒intron structures and motif compositions. Gene duplication, which included whole-genome duplication, tandem repeats, scattered repeats, and unconnected repeats, produced approximately 81 pairs of homologous genes in the AhCaM/CML gene family. Allopolyploidization was the main reason for the greater number of AhCaM/CML members. The nonsynonymous (Ka) versus synonymous (Ks) substitution rates (less than 1.0) suggested that all homologous pairs underwent intensive purifying selection pressure during evolution. AhCML69 was constitutively expressed in different tissues of peanut plants and was involved in the response to R. solanacearum infection. The AhCML69 protein was localized in the cytoplasm and nucleus. Transient overexpression of AhCML69 in tobacco leaves increased resistance to R. solanacearum infection and induced the expression of defense-related genes, suggesting that AhCML69 is a positive regulator of disease resistance.

CONCLUSIONS

This study provides the first comprehensive analysis of the AhCaM/CML gene family and potential genetic resources for the molecular design and breeding of peanut bacterial wilt resistance.

Contrasting plant transcriptome responses between a pierce-sucking and a chewing herbivore go beyond the infestation site

BACKGROUND

Plants have acquired a repertoire of mechanisms to combat biotic stressors, which may vary depending on the feeding strategies of herbivores and the plant species. Hormonal regulation crucially modulates this malleable defense response. Jasmonic acid (JA) and salicylic acid (SA) stand out as pivotal regulators of defense, while other hormones like abscisic acid (ABA), ethylene (ET), gibberellic acid (GA) or auxin also play a role in modulating plant-pest interactions. The plant defense response has been described to elicit effects in distal tissues, whereby aboveground herbivory can influence belowground response, and vice versa. This impact on distal tissues may be contingent upon the feeding guild, even affecting both the recovery of infested tissues and those that have not suffered active infestation.

RESULTS

To study how phytophagous with distinct feeding strategies may differently trigger the plant defense response during and after infestation in both infested and distal tissues, Arabidopsis thaliana L. rosettes were infested separately with the chewing herbivore Pieris brassicae L. and the piercing-sucker Tetranychus urticae Koch. Moderate infestation conditions were selected for both pests, though no quantitative control of damage levels was carried out. Feeding mode did distinctly influence the transcriptomic response of the plant under these conditions. Though overall affected processes were similar under either infestation, their magnitude differed significantly. Plants infested with P. brassicae exhibited a short-term response, involving stress-related genes, JA and ABA regulation and suppressing growth-related genes. In contrast, T. urticae elicited a longer transcriptomic response in plants, albeit with a lower degree of differential expression, in particular influencing SA regulation. These distinct defense responses transcended beyond infestation and through the roots, where hormonal response, flavonoid regulation or cell wall reorganization were differentially affected.

CONCLUSION

These outcomes confirm that the existent divergent transcriptomic responses elicited by herbivores employing distinct feeding strategies possess the capacity to extend beyond infestation and even affect tissues that have not been directly infested. This remarks the importance of considering the entire plant's response to localized biotic stresses.

Transcriptome profile of pecan scab resistant and susceptible trees from a pecan provenance collection

Pecan scab is a devastating disease that causes damage to pecan (Carya illinoinensis (Wangenh.) K. Koch) fruit and leaves. The disease is caused by the fungus Venturia effusa (G. Winter) and the main management practice for controlling the disease is by application of fungicides at 2-to-3-week intervals throughout the growing season. Besides disease-related yield loss, application of fungicides can result in considerable cost and increases the likelihood of fungicide resistance developing in the pathogen. Resistant cultivars are available for pecan growers; although, in several cases resistance has been overcome as the pathogen adapts to infect resistant hosts. Despite the importance of host resistance in scab management, there is little information regarding the molecular basis of genetic resistance to pecan scab.The purpose of this study was to elucidate mechanisms of natural pecan scab resistance by analyzing transcripts that are differentially expressed in pecan leaf samples from scab resistant and susceptible trees. The leaf samples were collected from trees in a provenance collection orchard that represents the natural range of pecan in the US and Mexico. Trees in the orchard have been exposed to natural scab infections since planting in 1989, and scab ratings were collected over three seasons. Based on this data, ten susceptible trees and ten resistant trees were selected for analysis. RNA-seq data was collected and analyzed for diseased and non-diseased parts of susceptible trees as well as for resistant trees. A total of 313 genes were found to be differentially expressed when comparing resistant and susceptible trees without disease. For susceptible samples showing scab symptoms, 1,454 genes were identified as differentially expressed compared to non-diseased susceptible samples. Many genes involved in pathogen recognition, defense responses, and signal transduction were up-regulated in diseased samples of susceptible trees, whereas differentially expressed genes in pecan scab resistant samples were generally down-regulated compared to non-diseased susceptible samples.Our results provide the first account of candidate genes involved in resistance/susceptibility to pecan scab under natural conditions in a pecan orchard. This information can be used to aid pecan breeding programs and development of biotechnology-based approaches for generating pecan cultivars with more durable scab resistance.

Transcriptome analysis in Aegilops tauschii unravels further insights into genetic control of stripe rust resistance

MAIN CONCLUSION

The Aegilops tauschii resistant accession prevented the pathogen colonization by controlling the sugar flow and triggering the hypersensitive reaction. This study suggested that NBS-LRRs probably induce resistance through bHLH by controlling JA- and SA-dependent pathways. Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst) is one of wheat's most destructive fungal diseases that causes a severe yield reduction worldwide. The most effective and economically-friendly strategy to manage this disease is genetic resistance which can be achieved through deploying new and effective resistance genes. Aegilops tauschii, due to its small genome and co-evolution with Pst, can provide detailed information about underlying resistance mechanisms. Hence, we used RNA-sequencing approach to identify the transcriptome variations of two contrasting resistant and susceptible Ae. tauschii accessions in interaction with Pst and differentially expressed genes (DEGs) for resistance to stripe rust. Gene ontology, pathway analysis, and search for functional domains, transcription regulators, resistance genes, and protein-protein interactions were used to interpret the results. The genes encoding NBS-LRR, CC-NBS-kinase, and phenylalanine ammonia-lyase, basic helix-loop-helix (bHLH)-, basic-leucine zipper (bZIP)-, APETALA2 (AP2)-, auxin response factor (ARF)-, GATA-, and LSD-like transcription factors were up-regulated exclusively in the resistant accession. The key genes involved in response to salicylic acid, amino sugar and nucleotide sugar metabolism, and hypersensitive response contributed to plant defense against stripe rust. The activation of jasmonic acid biosynthesis and starch and sucrose metabolism pathways under Pst infection in the susceptible accession explained the colonization of the host. Overall, this study can fill the gaps in the literature on host-pathogen interaction and enrich the Ae. tauschii transcriptome sequence information. It also suggests candidate genes that could guide future breeding programs attempting to develop rust-resistant cultivars.

Mining of disease-resistance genes in Crocus sativus based on transcriptome sequencing

Introduction: Crocus sativus L. has an important medicinal and economic value in traditional perennial Chinese medicine. However, due to its unique growth characteristics, during cultivation it is highly susceptible to disease. The absence of effective resistance genes restricts us to breed new resistant varieties of C. sativus. Methods: In present study, comprehensive transcriptome sequencing was introduced to explore the disease resistance of the candidate gene in healthy and corm rot-infected C. sativus. Results and discussion: Totally, 43.72 Gb of clean data was obtained from the assembly to generate 65,337 unigenes. By comparing the gene expression levels, 7,575 differentially expressed genes (DEGs) were primarily screened. A majority of the DEGs were completely in charge of defense and metabolism, and 152 of them were annotated as pathogen recognition genes (PRGs) based on the PGRdb dataset. The expression of some transcription factors including NAC, MYB, and WRKY members, changed significantly based on the dataset of transcriptome sequencing. Therefore, this study provides us some valuable information for exploring candidate genes involved in the disease resistance in C. sativus.

Genomic evolution and insights into agronomic trait innovations of Sesamum species

Sesame is an ancient oilseed crop with high oil content and quality. However, the evolutionary history and genetic mechanisms of its valuable agronomic traits remain unclear. Here, we report chromosome-scale genomes of cultivated sesame (Sesamum indicum L.) and six wild Sesamum species, representing all three karyotypes within this genus. Karyotyping and genome-based phylogenic analysis revealed the evolutionary route of Sesamum species from n = 13 to n = 16 and revealed that allotetraploidization occurred in the wild species Sesamum radiatum. Early divergence of the Sesamum genus (48.5-19.7 million years ago) during the Tertiary period and its ancient phylogenic position within eudicots were observed. Pan-genome analysis revealed 9164 core gene families in the 7 Sesamum species. These families are significantly enriched in various metabolic pathways, including fatty acid (FA) metabolism and FA biosynthesis. Structural variations in SiPT1 and SiDT1 within the phosphatidyl ethanolamine-binding protein gene family lead to the genomic evolution of plant-architecture and inflorescence-development phenotypes in Sesamum. A genome-wide association study (GWAS) of an interspecific population and genome comparisons revealed a long terminal repeat insertion and a sequence deletion in DIR genes of wild Sesamum angustifolium and cultivated sesame, respectively; both variations independently cause high susceptibility to Fusarium wilt disease. A GWAS of 560 sesame accessions combined with an overexpression study confirmed that the NAC1 and PPO genes play an important role in upregulating oil content of sesame. Our study provides high-quality genomic resources for cultivated and wild Sesamum species and insights that can improve molecular breeding strategies for sesame and other oilseed crops.

A haplotype-resolved chromosome-scale genome for Quercus rubra L. provides insights into the genetics of adaptive traits for red oak species

Northern red oak (Quercus rubra L.) is an ecologically and economically important forest tree native to North America. We present a chromosome-scale genome of Q. rubra generated by the combination of PacBio sequences and chromatin conformation capture (Hi-C) scaffolding. This is the first reference genome from the red oak clade (section Lobatae). The Q. rubra assembly spans 739 Mb with 95.27% of the genome in 12 chromosomes and 33,333 protein-coding genes. Comparisons to the genomes of Quercus lobata and Quercus mongolica revealed high collinearity, with intrachromosomal structural variants present. Orthologous gene family analysis with other tree species revealed that gene families associated with defense response were expanding and contracting simultaneously across the Q. rubra genome. Quercus rubra had the most CC-NBS-LRR and TIR-NBS-LRR resistance genes out of the 9 species analyzed. Terpene synthase gene family comparisons further reveal tandem gene duplications in TPS-b subfamily, similar to Quercus robur. Phylogenetic analysis also identified 4 subfamilies of the IGT/LAZY gene family in Q. rubra important for plant structure. Single major QTL regions were identified for vegetative bud break and marcescence, which contain candidate genes for further research, including a putative ortholog of the circadian clock constituent cryptochrome (CRY2) and 8 tandemly duplicated genes for serine protease inhibitors, respectively. Genome-environment associations across natural populations identified candidate abiotic stress tolerance genes and predicted performance in a common garden. This high-quality red oak genome represents an essential resource to the oak genomic community, which will expedite comparative genomics and biological studies in Quercus species.

Cell-type-specific responses to fungal infection in plants revealed by single-cell transcriptomics

Pathogen infection is a dynamic process. Here, we employ single-cell transcriptomics to investigate plant response heterogeneity. By generating an Arabidopsis thaliana leaf atlas encompassing 95,040 cells during infection by a fungal pathogen, Colletotrichum higginsianum, we unveil cell-type-specific gene expression, notably an enrichment of intracellular immune receptors in vasculature cells. Trajectory inference identifies cells that had different interactions with the invading fungus. This analysis divulges transcriptional reprogramming of abscisic acid signaling specifically occurring in guard cells, which is consistent with a stomatal closure dependent on direct contact with the fungus. Furthermore, we investigate the transcriptional plasticity of genes involved in glucosinolate biosynthesis in cells at the fungal infection sites, emphasizing the contribution of the epidermis-expressed MYB122 to disease resistance. This work underscores spatially dynamic, cell-type-specific plant responses to a fungal pathogen and provides a valuable resource that supports in-depth investigations of plant-pathogen interactions.

Development and Genome-Wide Analysis of a Blast-Resistant japonica Rice Variety

Rice is one of the most important crops in the world, and its production is severely affected by the rice blast disease caused by the fungus Magnaporthe oryzae. Several major blast resistance genes and QTLs associated with blast resistance have been described and mostly identified in indica rice varieties. In this work, we report the obtention of a blast-resistant rice breeding line derived from crosses between the resistant indica variety CT13432 and the japonica elite cultivar JSendra (highly susceptible to blast). The breeding line, named COPSEMAR9, was found to exhibit resistance to leaf blast and panicle blast, as demonstrated by disease assays under controlled and field conditions. Furthermore, a high-quality genome sequence of the blast-resistant breeding line was obtained using a strategy that combines short-read sequencing (Illumina sequencing) and long-read sequencing (Pacbio sequencing). The use of a whole-genome approach allowed the fine mapping of DNA regions of indica and japonica origin present in the COPSEMAR9 genome and the identification of parental gene regions potentially contributing to blast resistance in the breeding line. Rice blast resistance genes (including Pi33 derived from the resistant parent) and defense-related genes in the genome of COPSEMAR9 were identified. Whole-genome analyses also revealed the presence of microRNAs (miRNAs) with a known function in the rice response to M. oryzae infection in COPSEMAR9, which might also contribute to its phenotype of blast resistance. From this study, the genomic information and analysis methods provide valuable knowledge that will be useful in breeding programs for blast resistance in japonica rice cultivars.

Population genomics identifies genetic signatures of carrot domestication and improvement and uncovers the origin of high-carotenoid orange carrots

Here an improved carrot reference genome and resequencing of 630 carrot accessions were used to investigate carrot domestication and improvement. The study demonstrated that carrot was domesticated during the Early Middle Ages in the region spanning western Asia to central Asia, and orange carrot was selected during the Renaissance period, probably in western Europe. A progressive reduction of genetic diversity accompanied this process. Genes controlling circadian clock/flowering and carotenoid accumulation were under selection during domestication and improvement. Three recessive genes, at the REC, Or and Y2 quantitative trait loci, were essential to select for the high α- and β-carotene orange phenotype. All three genes control high α- and β-carotene accumulation through molecular mechanisms that regulate the interactions between the carotenoid biosynthetic pathway, the photosynthetic system and chloroplast biogenesis. Overall, this study elucidated carrot domestication and breeding history and carotenoid genetics at a molecular level.

Understanding the mechanism of pathogenicity through interactome studies between Arachis hypogaea L. and Aspergillus flavus

Aspergillus flavus (A. flavus) infects the peanut seeds during pre-and post-harvest stages, causing seed quality destruction for humans and livestock consumption. Even though many resistant varieties were developed, the molecular mechanism of defense interactions of peanut against A. flavus still needs further investigation. Hence, an interologous host-pathogen protein interaction (HPPI) network was constructed to understand the subcellular level interaction mechanism between peanut and A. flavus. Out of the top 10 hub proteins of both organisms, protein phosphatase 2C and cyclic nucleotide-binding/kinase domain-containing protein and different ribosomal proteins were identified as candidate proteins involved in defense. Functional annotation and subcellular localization based characterization of HPPI identified protein SGT1 homolog, calmodulin and Rac-like GTP-binding proteins to be involved in defense response against fungus. The relevance of HPPI in infectious conditions was assessed using two transcriptome data which identified the interplay of host kinase class R proteins, bHLH TFs and cell wall related proteins to impart resistance against pathogen infection. Further, the pathogenicity analysis identified glycogen phosphorylase and molecular chaperone and allergen Mod-E/Hsp90/Hsp1 as potential pathogen targets to enhance the host defense mechanism. Hence, the computationally predicted host-pathogen PPI network could provide valuable support for molecular biology experiments to understand the host-pathogen interaction. SIGNIFICANCE: Protein-protein interactions execute significant cellular interactions in an organism and are influenced majorly by stress conditions. Here we reported the host-pathogen protein-protein interaction between peanut and A. flavus, and a detailed network analysis based on function, subcellular localization, gene co-expression, and pathogenicity was performed. The network analysis identified key proteins such as host kinase class R proteins, calmodulin, SGT1 homolog, Rac-like GTP-binding proteins bHLH TFs and cell wall related to impart resistance against pathogen infection. We observed the interplay of defense related proteins and cell wall related proteins predominantly, which could be subjected to further studies. The network analysis described in this study could be applied to understand other host-pathogen systems generally.

Long-read, chromosome-scale assembly of Vitis rotundifolia cv. Carlos and its unique resistance to Xylella fastidiosa subsp. fastidiosa

BACKGROUND

Muscadine grape (Vitis rotundifolia) is resistant to many of the pathogens that negatively impact the production of common grape (V. vinifera), including the bacterial pathogen Xylella fastidiosa subsp. fastidiosa (Xfsf), which causes Pierce's Disease (PD). Previous studies in common grape have indicated Xfsf delays host immune response with a complex O-chain antigen produced by the wzy gene. Muscadine cultivars range from tolerant to completely resistant to Xfsf, but the mechanism is unknown.

RESULTS

We assembled and annotated a new, long-read genome assembly for 'Carlos', a cultivar of muscadine that exhibits tolerance, to build upon the existing genetic resources available for muscadine. We used these resources to construct an initial pan-genome for three cultivars of muscadine and one cultivar of common grape. This pan-genome contains a total of 34,970 synteny-constrained entries containing genes of similar structure. Comparison of resistance gene content between the 'Carlos' and common grape genomes indicates an expansion of resistance (R) genes in 'Carlos.' We further identified genes involved in Xfsf response by transcriptome sequencing 'Carlos' plants inoculated with Xfsf. We observed 234 differentially expressed genes with functions related to lipid catabolism, oxidation-reduction signaling, and abscisic acid (ABA) signaling as well as seven R genes. Leveraging public data from previous experiments of common grape inoculated with Xfsf, we determined that most differentially expressed genes in the muscadine response were not found in common grape, and three of the R genes identified as differentially expressed in muscadine do not have an ortholog in the common grape genome.

CONCLUSIONS

Our results support the utility of a pan-genome approach to identify candidate genes for traits of interest, particularly disease resistance to Xfsf, within and between muscadine and common grape.

Genome engineering of disease susceptibility genes for enhancing resistance in plants

Introgression of disease resistance genes (R-genes) to fight against an array of phytopathogens takes several years using conventional breeding approaches. Pathogens develop mechanism(s) to escape plants immune system by evolving new strains/races, thus making them susceptible to disease. Conversely, disruption of host susceptibility factors (or S-genes) provides opportunities for resistance breeding in crops. S-genes are often exploited by phytopathogens to promote their growth and infection. Therefore, identification and targeting of disease susceptibility genes (S-genes) are gaining more attention for the acquisition of resistance in plants. Genome engineering of S-genes results in targeted, transgene-free gene modification through CRISPR-Cas-mediated technology and has been reported in several agriculturally important crops. In this review, we discuss the defense mechanism in plants against phytopathogens, tug of war between R-genes and S-genes, in silico techniques for identification of host-target (S-) genes and pathogen effector molecule(s), CRISPR-Cas-mediated S-gene engineering, its applications, challenges, and future prospects.

Chromosome-scale genome assembly-assisted identification of Mi-9 gene in Solanum arcanum accession LA2157, conferring heat-stable resistance to Meloidogyne incognita

Root-knot nematodes (RKNs) are infamous plant pathogens in tomato production, causing considerable losses in agriculture worldwide. Mi-1 is the only commercially available RKN-resistance gene; however, the resistance is inactivated when the soil temperature is over 28 °C. Mi-9 in wild tomato (Solanum arcanum LA2157) has stable resistance to RKNs under high temperature but has not been cloned and applied. In this study, a chromosome-scale genome assembly of S. arcanum LA2157 was constructed through Nanopore and Hi-C sequencing. Based on molecular markers of Mi-9 and comparative genomic analysis, the localization region and candidate Mi-9 genes cluster consisting of seven nucleotide-binding sites and leucine-rich repeat (NBS-LRR) genes were located. Transcriptional expression profiles confirmed that five of the seven candidate genes were expressed in root tissue. Moreover, virus-induced gene silencing of the Sarc_034200 gene resulted in increased susceptibility of S. arcanum LA2157 to Meloidogyne incognita, and genetic transformation of the Sarc_034200 gene in susceptible Solanum pimpinellifolium conferred significant resistance to M. incognita at 25 °C and 30 °C and showed hypersensitive responses at nematode infection sites. This suggested that Sarc_034200 is the Mi-9 gene. In summary, we cloned, confirmed and applied the heat-stable RKN-resistance gene Mi-9, which is of great significance to tomato breeding for nematode resistance.

Genome wide identification and evolutionary analysis of vat like NBS-LRR genes potentially associated with resistance to aphids in cotton

Nucleotide Binding Site - Leucine Rich Repeat (NBS-LRR) genes play a significant role in plant defense against biotic stresses and are an integral part of signal transduction pathways. Vat gene has been well reported for their role in resistance to Aphis gossypii and viruses transmitted by them. Despite their importance, Vat like NBS-LRR resistance genes have not yet been identified and studied in cotton species. This study report hundreds of orthologous Vat like NBS-LRR genes from the genomes of 18 cotton species through homology searches and the distribution of those identified genes were tend to be clustered on different chromosome. Especially, in a majority of the cases, Vat like genes were located on chromosome number 13 and they all shared two conserved NBS-LRR domains, one disease resistant domain and several repeats of LRR on the investigated cotton Vat like proteins. Gene ontology study on Vat like NBS-LRR genes revealed the molecular functions viz., ADP and protein binding. Phylogenetic analysis also revealed that Vat like sequences of two diploid species, viz., G. arboreum and G. anomalum, were closely related to the sequences of the tetraploids than all other diploids. The Vat like genes of G. aridum and G. schwendimanii were distantly related among diploids and tetraploids species. Various hormones and defense related cis-acting regulatory elements were identified from the 2 kb upstream sequences of the Vat like genes implying their defensive response towards the biotic stresses. Interestingly, G. arboreum and G. trilobum were found to have more regulatory elements than larger genomes of tetraploid cotton species. Thus, the present study provides the evidence for the evolution of Vat like genes in defense mechanisms against aphids infestation in cotton genomes and allows further characterization of candidate genes for developing aphid and aphid transmitted viruses resistant crops through cotton breeding.

Dual Transcriptome Analysis Reveals That ChATG8 Is Required for Fungal Development, Melanization and Pathogenicity during the Interaction between Colletotrichum higginsianum and Arabidopsis thaliana

Anthracnose disease of cruciferous plants caused by Colletotrichum higginsianum is a serious fungal disease that affects cruciferous crops such as Chinese cabbage, Chinese flowering cabbage, broccoli, mustard plant, as well as the model plant Arabidopsis thaliana. Dual transcriptome analysis is commonly used to identify the potential mechanisms of interaction between host and pathogen. In order to identify differentially expressed genes (DEGs) in both the pathogen and host, the conidia of wild-type (ChWT) and Chatg8 mutant (Chatg8Δ) strains were inoculated onto leaves of A. thaliana, and the infected leaves of A. thaliana at 8, 22, 40, and 60 h post-inoculation (hpi) were subjected to dual RNA-seq analysis. The results showed that comparison of gene expression between the 'ChWT' and 'Chatg8Δ' samples detected 900 DEGs (306 upregulated and 594 down-regulated) at 8 hpi, 692 DEGs (283 upregulated and 409 down-regulated) at 22 hpi, 496 DEGs (220 upregulated and 276 down-regulated) at 40 hpi, and 3159 DEGs (1544 upregulated and 1615 down-regulated) at 60 hpi. GO and KEGG analyses found that the DEGs were mainly involved in fungal development, biosynthesis of secondary metabolites, plant-fungal interactions, and phytohormone signaling. The regulatory network of key genes annotated in the Pathogen-Host Interactions database (PHI-base) and Plant Resistance Genes database (PRGdb), as well as a number of key genes highly correlated with the 8, 22, 40, and 60 hpi, were identified during the infection. Among the key genes, the most significant enrichment was in the gene encoding the trihydroxynaphthalene reductase (THR1) in the melanin biosynthesis pathway. Both Chatg8Δ and Chthr1Δ strains showed varying degrees of reduction of melanin in appressoria and colonies. The pathogenicity of the Chthr1Δ strain was lost. In addition, six DEGs from C. higginsianum and six DEGs from A. thaliana were selected for real-time quantitative PCR (RT-qPCR) to confirm the RNA-seq results. The information gathered from this study enriches the resources available for research into the role of the gene ChATG8 during the infection of A. thaliana by C. higginsianum, such as potential links between melanin biosynthesis and autophagy, and the response of A. thaliana to different fungal strains, thereby providing a theoretical basis for the breeding of cruciferous green leaf vegetable cultivars with resistance to anthracnose disease.

Identification of QTLs and critical genes related to sugarcane mosaic disease resistance

Mosaic viral diseases affect sugarcane productivity worldwide. Mining disease resistance-associated molecular markers or genes is a key component of disease resistance breeding programs. In the present study, 285 F₁ progeny were produced from a cross between Yuetang 93-159, a moderately resistant variety, and ROC22, a highly susceptible variety. The mosaic disease symptoms of these progenies, with ROC22 as the control, were surveyed by natural infection under 11 different environmental conditions in the field and by artificial infections with a mixed sugarcane mosaic virus (SCMV) and sorghum mosaic virus (SrMV) inoculum. Analysis of consolidated survey data enabled the identification of 29 immune, 55 highly resistant, 70 moderately resistant, 62 susceptible, and 40 highly susceptible progenies. The disease response data and a high-quality SNP genetic map were used in quantitative trait locus (QTL) mapping. The results showed that the correlation coefficients (0.26~0.91) between mosaic disease resistance and test environments were significant (p< 0.001), and that mosaic disease resistance was a highly heritable quantitative trait (H² = 0.85). Seven mosaic resistance QTLs were located to the SNP genetic map, each QTL accounted for 3.57% ~ 17.10% of the phenotypic variation explained (PVE). Furthermore, 110 pathogen response genes and 69 transcription factors were identified in the QTLs interval. The expression levels of nine genes (Soffic.07G0015370-1P, Soffic.09G0015410-2T, Soffic.09G0016460-1T, Soffic.09G0016460-1P, Soffic.09G0017080-3C, Soffic.09G0018730-3P, Soffic.09G0018730-3C, Soffic.09G0019920-3C and Soffic.03G0019710-2C) were significantly different between resistant and susceptible progenies, indicating their key roles in sugarcane resistance to SCMV and SrMV infection. The seven QTLs and nine genes can provide a certain scientific reference to help sugarcane breeders develop varieties resistant to mosaic diseases.

Autopolyploid inheritance and a heterozygous reciprocal translocation shape chromosome genetic behavior in tetraploid blueberry (Vaccinium corymbosum)

Understanding chromosome recombination behavior in polyploidy species is key to advancing genetic discoveries. In blueberry, a tetraploid species, the line of evidences about its genetic behavior still remain poorly understood, owing to the inter-specific, and inter-ploidy admixture of its genome and lack of in depth genome-wide inheritance and comparative structural studies. Here we describe a new high-quality, phased, chromosome-scale genome of a diploid blueberry, clone W85. The genome was integrated with cytogenetics and high-density, genetic maps representing six tetraploid blueberry cultivars, harboring different levels of wild genome admixture, to uncover recombination behavior and structural genome divergence across tetraploid and wild diploid species. Analysis of chromosome inheritance and pairing demonstrated that tetraploid blueberry behaves as an autotetraploid with tetrasomic inheritance. Comparative analysis demonstrated the presence of a reciprocal, heterozygous, translocation spanning one homolog of chr-6 and one of chr-10 in the cultivar Draper. The translocation affects pairing and recombination of chromosomes 6 and 10. Besides the translocation detected in Draper, no other structural genomic divergences were detected across tetraploid cultivars and highly inter-crossable wild diploid species. These findings and resources will facilitate new genetic and comparative genomic studies in Vaccinium and the development of genomic assisted selection strategy for this crop.

LDRGDb - Legumes disease resistance genes database

Legumes comprise one of the world's largest, most diverse, and economically important plant families, known for their nutritional and medicinal benefits. Legumes are susceptible to a wide range of diseases, similar to other agricultural crops. Diseases have a considerable impact on the production of legume crop species, resulting in large yield losses worldwide. Due to continuous interactions between plants and their pathogens in the environment and the evolution of new pathogens under high selection pressure; disease resistant genes emerge in plant cultivars in the field against those pathogens or disease. Thus, disease resistant genes play critical roles in plant resistance responses, and their discovery and subsequent use in breeding programmes aid in reducing yield loss. The genomic era, with its high-throughput and low-cost genomic tools, has revolutionised our understanding of the complex interactions between legumes and pathogens, resulting in the identification of several critical participants in both the resistant and susceptible relationships. However, a substantial amount of existing information about numerous legume species has been disseminated as text or is preserved across fractions in different databases, posing a challenge for researchers. As a result, the range, scope, and complexity of these resources pose challenges to those who manage and use them. Therefore, there is an urgent need to develop tools and a single conjugate database to manage genetic information for the world's plant genetic resources, allowing for the rapid incorporation of essential resistance genes into breeding strategies. Here, developed the first comprehensive database of disease resistance genes named as LDRGDb - LEGUMES DISEASE RESISTANCE GENES DATABASE comprises 10 legumes [Pigeon pea (Cajanus cajan), Chickpea (Cicer arietinum), Soybean (Glycine max), Lentil (Lens culinaris), Alfalfa (Medicago sativa), Barrelclover (Medicago truncatula), Common bean (Phaseolus vulgaris), Pea (Pisum sativum),Faba bean (Vicia faba), and Cowpea (Vigna unguiculata)]. The LDRGDb is a user-friendly database developed by integrating a variety of tools and software that combine knowledge about resistant genes, QTLs, and their loci, with proteomics, pathway interactions, and genomics (https://ldrgdb.in/).

Genome-Wide Association Studies in Sunflower: Towards Sclerotinia sclerotiorum and Diaporthe/Phomopsis Resistance Breeding

Diseases caused by necrotrophic fungi, such as the cosmopolitan Sclerotinia sclerotiorum and the Diaporthe/Phomopsis complex, are among the most destructive diseases of sunflower worldwide. The lack of complete resistance combined with the inefficiency of chemical control makes assisted breeding the best strategy for disease control. In this work, we present an integrated genome-wide association (GWA) study investigating the response of a diverse panel of sunflower inbred lines to both pathogens. Phenotypic data for Sclerotinia head rot (SHR) consisted of five disease descriptors (disease incidence, DI; disease severity, DS; area under the disease progress curve for DI, AUDPCI, and DS, AUDPCS; and incubation period, IP). Two disease descriptors (DI and DS) were evaluated for two manifestations of Diaporthe/Phomopsis: Phomopsis stem canker (PSC) and Phomopsis head rot (PHR). In addition, a principal component (PC) analysis was used to derive transformed phenotypes as inputs to a univariate GWA (PC-GWA). Genotypic data comprised a panel of 4269 single nucleotide polymorphisms (SNP), generated via genotyping-by-sequencing. The GWA analysis revealed 24 unique marker-trait associations for SHR, 19 unique marker-trait associations for Diaporthe/Phomopsis diseases, and 7 markers associated with PC1 and PC2. No common markers were found for the response to the two pathogens. Nevertheless, epistatic interactions were identified between markers significantly associated with the response to S. sclerotiorum and Diaporthe/Phomopsis. This suggests that, while the main determinants of resistance may differ for the two pathogens, there could be an underlying common genetic basis. The exploration of regions physically close to the associated markers yielded 364 genes, of which 19 were predicted as putative disease resistance genes. This work presents the first simultaneous evaluation of two manifestations of Diaporthe/Phomopsis in sunflower, and undertakes a comprehensive GWA study by integrating PSC, PHR, and SHR data. The multiple regions identified, and their exploration to identify candidate genes, contribute not only to the understanding of the genetic basis of resistance, but also to the development of tools for assisted breeding.

Chromosome-scale Amaranthus tricolor genome provides insights into the evolution of the genus Amaranthus and the mechanism of betalain biosynthesis

Amaranthus tricolor is a vegetable and ornamental amaranth, with high lysine, dietary fibre and squalene content. The red cultivar of A. tricolor possesses a high concentration of betalains, which has been used as natural food colorants. Here, we constructed the genome of A. tricolor, the first reference genome for the subgenus Albersia, combining PacBio HiFi, Nanopore ultra-long and Hi-C data. The contig N50 size was 906 kb, and 99.58% of contig sequence was anchored to the 17 chromosomes, totalling 520 Mb. We annotated 27,813 protein-coding genes with an average 1.3 kb coding sequence and 5.3 exons. We inferred that A. tricolor underwent a whole-genome duplication (WGD) and that the WGD shared by amaranths occurred in the last common ancestor of subfamily Amaranthoideae. Moreover, we comprehensively identified candidate genes in betalain biosynthesis pathway. Among them, DODAα1 and CYP76ADα1, located in one topologically associated domain (TAD) of an active (A) compartment on chromosome 16, were more highly expressed in red leaves than in green leaves, and DODAα1 might be the rate-limiting enzyme gene in betalains biosynthesis. This study presents new genome resources and enriches our understanding of amaranth evolution, betalains production, facilitating molecular breeding improvements and the understanding of C4 plants evolution.

De novo assembly and characterization of the draft genome of the cashew (Anacardium occidentale L.)

Cashew is the second most important tree nut crop in the global market. Cashew is a diploid and heterozygous species closely related to the mango and pistachio. Its improvement by conventional breeding is slow due to the long juvenile phase. Despite the economic importance, very little genomics/transcriptomics information is available for cashew. In this study, the Oxford nanopore reads and Illumina reads were used for de novo assembly of the cashew genome. The hybrid assembly yielded a 356.6 Mb genome corresponding to 85% of the estimated genome size (419 Mb). The BUSCO analysis showed 91.8% of genome completeness. Transcriptome mapping showed 92.75% transcripts aligned with the assembled genome. Gene predictions resulted in the identification of 31,263 genes coding for a total of 35,000 gene isoforms. About 46% (165 Mb) of the cashew genome comprised of repetitive sequences. Phylogenetic analyses of the cashew with nine species showed that it was closely related to Mangifera indica. Analysis of cashew genome revealed 3104 putative R-genes. The first draft assembly of the genome, transcriptome and R gene information generated in this study would be the foundation for understanding the molecular basis of economic traits and genomics-assisted breeding in cashew.

Comparative transcriptome profiling and weighted gene co-expression network analysis to identify core genes in maize (Zea mays L.) silks infected by multiple fungi

Maize (Zea mays L.) is the third most popular Poaceae crop after wheat and rice and used in feed and pharmaceutical sectors. The maize silk contains bioactive components explored by traditional Chinese herbal medicine for various pharmacological activities. However, Fusarium graminearum, Fusarium verticillioides, Trichoderma atroviride, and Ustilago maydis can infect the maize, produce mycotoxins, hamper the quantity and quality of silk production, and further harm the primary consumer's health. However, the defense mechanism is not fully understood in multiple fungal infections in the silk of Z. mays. In this study, we applied bioinformatics approaches to use the publicly available transcriptome data of Z. mays silk affected by multiple fungal flora to identify core genes involved in combatting disease response. Differentially expressed genes (DEGs) were identified among intra- and inter-transcriptome data sets of control versus infected Z. mays silks. Upon further comparison between up- and downregulated genes within the control of datasets, 4,519 upregulated and 5,125 downregulated genes were found. The DEGs have been compared with genes in the modules of weighted gene co-expression network analysis to relevant specific traits towards identifying core genes. The expression pattern of transcription factors, carbohydrate-active enzymes (CAZyme), and resistance genes was analyzed. The present investigation is supportive of our findings that the gene ontology, immunity stimulus, and resistance genes are upregulated, but physical and metabolic processes such as cell wall organizations and pectin synthesis were downregulated respectively. Our results are indicative that terpene synthase TPS6 and TPS11 are involved in the defense mechanism against fungal infections in maize silk.

Predicción bioinformática de proteínas NBS-LRR en el genoma de Coffea arabica

Gracias al acceso al genoma completo de Coffea arabica y el Desarrollo de multiples herramientas de bioinformartica que permite la búsqueda de genes de resistencia de plantas (R-genes), ha sido posible implementar estas estrategias en programas de mejora genética. En las plantas, los R-genes codifican proteínas involucradas en mecanismos de defensa contra patógenos. Los genes con dominios tipo Nucleotide-Binding-Site Leucine-Rich-Repeat (NBS-LRR) forman la familia de R-genes de plantas más grande. El objetivo de este estudio fue identificar genes de proteínas NBS-LRR en el genoma de C. arabica utilizando un enfoque bioinformático. Identificamos motivos conservados de R-genes de C. arabica relacionados con genes similares encontrados en Coffea canephora y Coffea eugenoides, dos especies evolutivas relacionadas con C. arabica. Los resultados de estos análisis revelaron proteínas con origen evolutivo provenientes de dicotiledóneo ancestrales, así como proteínas de resistencia específicas del género Coffea. Además, todas las secuencias de los R-genes de C. arabica mostraron una gran similitud con proteína CNL de Arabidopsis thaliana. Finalmente, la presencia de motivos altamente conservados, la distribución cromosómica y las relaciones filogenéticas de los R-genes de C. arabica muestran procesos de coevolución con patógenos adaptados, demostrando de esta manera la importancia del estudio de estos genes en la inmunidad del café. Palabras clave: Café, NBS-LRR, Proteínas de Resistencia, Bioinformática.

CarrotOmics: a genetics and comparative genomics database for carrot (Daucus carota)

CarrotOmics (https://carrotomics.org/) is a comprehensive database for carrot (Daucus carota L.) breeding and research. CarrotOmics was developed using resources available at the MainLab Bioinformatics core (https://www.bioinfo.wsu.edu/) and is implemented using Tripal with Drupal modules. The database delivers access to download or visualize the carrot reference genome with gene predictions, gene annotations and sequence assembly. Other genomic resources include information for 11 224 genetic markers from 73 linkage maps or genotyping-by-sequencing and descriptions of 371 mapped loci. There are records for 1601 Apiales species (or subspecies) and descriptions of 9408 accessions from 11 germplasm collections representing more than 600 of these species. Additionally, 204 Apiales species have phenotypic information, totaling 28 517 observations from 10 041 biological samples. Resources on CarrotOmics are freely available, search functions are provided to find data of interest and video tutorials are available to describe the search functions and genomic tools. CarrotOmics is a timely resource for the Apiaceae research community and for carrot geneticists developing improved cultivars with novel traits addressing challenges including an expanding acreage in tropical climates, an evolving consumer interested in sustainably grown vegetables and a dynamic environment due to climate change. Data from CarrotOmics can be applied in genomic-assisted selection and genetic research to improve basic research and carrot breeding efficiency.

Database URL

https://carrotomics.org/

Long-read-based draft genome sequence of Indian black gram IPU-94-1 'Uttara': Insights into disease resistance and seed storage protein genes

Black gram [Vigna mungo (L.) Hepper var. mungo] is a warm-season legume highly prized for its protein content along with significant folate and iron proportions. To expedite the genetic enhancement of black gram, a high-quality draft genome from the center of origin of the crop is indispensable. Here, we established a draft genome sequence of an Indian black gram cultivar, 'Uttara' (IPU 94-1), known for its high resistance to mungbean yellow mosaic virus. Pacific Biosciences of California, Inc. (PacBio) single-molecule real-time (SMRT) and Illumina sequencing assembled a draft reference-guided assembly with a cumulative size of ∼454.4 Mb, of which, 444.4 Mb was anchored on 11 pseudomolecules corresponding to 11 chromosomes. Uttara assembly denotes features of a high-quality draft genome illustrated through high N50 value (42.88 Mb), gene completeness (benchmarking universal single-copy ortholog [BUSCO] score 94.17%) and low levels of ambiguous nucleotides (N) percent (0.0005%). Gene discovery using transcript evidence predicted 28,881 protein-coding genes, from which, ∼95% were functionally annotated. A global survey of genes associated with disease resistance revealed 119 nucleotide binding site-leucine rich repeat (NBS-LRR) proteins, while 23 genes encoding seed storage proteins (SSPs) were discovered in black gram. A large set of microsatellite loci were discovered for marker development in the crop. Our draft genome of an Indian black gram provides the foundational genomic resources for the improvement of important agronomic traits and ultimately will help in accelerating black gram breeding programs.

MusaRgeneDB: an online comprehensive database for disease resistance genes in Musa spp

Banana is one of the major food crops and its production is subject to many pests and diseases. Banana breeding exploits wild relatives and progenitor species for the introgression of resistant genes (R) into cultivated varieties to overcome these hurdles. With advances in sequencing technologies, whole-genome sequences are available for many Musa spp. and many of them are potential donors of disease resistance genes. Considering their potential role, R genes from these species were explored to develop an user-friendly open-access database that will be useful for studying and implementing disease resistance in bananas. MusaRgene database is complemented with complete details of 3598 R genes identified from eight Musa spp. and rice, Arabidopsis, sorghum along with its classification and separate modules on its expression under various stresses in resistant and susceptible cultivars and corresponding SSRs are also provided. This database can be regarded as the primary resource of information on R genes from bananas and their relatives. R genes from other allele mining studies are also incorporated which will enable the identification of its homolog in related Musa spp. MusaRgene database will aid in the identification of genes and markers associated, cloning of full-length R genes, and genetic transformation or gene editing of the R genes in susceptible cultivars. Multiple R genes can also be identified for pyramiding the genes to increase the level of resistance and durability. Overall, this database will facilitate the understanding of defense mechanisms in bananas against biotic or abiotic stresses leading to the development of promising disease-resistant varieties.

Comparative transcriptome analysis revealed molecular mechanisms of peanut leaves responding to Ralstonia solanacearum and its type III secretion system mutant

Bacterial wilt caused by Ralstonia solanacearum is a serious soil-borne disease that limits peanut production and quality, but the molecular mechanisms of the peanut response to R. solanacearum remain unclear. In this study, we reported the first work analyzing the transcriptomic changes of the resistant and susceptible peanut leaves infected with R. solanacearum HA4-1 and its type III secretion system mutant strains by the cutting leaf method at different timepoints (0, 24, 36, and 72 h post inoculation). A total of 125,978 differentially expressed genes (DEGs) were identified and subsequently classified into six groups to analyze, including resistance-response genes, susceptibility-response genes, PAMPs induced resistance-response genes, PAMPs induced susceptibility-response genes, T3Es induced resistance-response genes, and T3Es induced susceptibility-response genes. KEGG enrichment analyses of these DEGs showed that plant-pathogen interaction, plant hormone signal transduction, and MAPK signaling pathway were the outstanding pathways. Further analysis revealed that CMLs/CDPKs-WRKY module, MEKK1-MKK2-MPK3 cascade, and auxin signaling played important roles in the peanut response to R. solanacearum. Upon R. solanacearum infection (RSI), three early molecular events were possibly induced in peanuts, including Ca²⁺ activating CMLs/CDPKs-WRKY module to regulate the expression of resistance/susceptibility-related genes, auxin signaling was induced by AUX/IAA-ARF module to activate auxin-responsive genes that contribute to susceptibility, and MEKK1-MKK2-MPK3-WRKYs was activated by phosphorylation to induce the expression of resistance/susceptibility-related genes. Our research provides new ideas and abundant data resources to elucidate the molecular mechanism of the peanut response to R. solanacearum and to further improve the bacterial wilt resistance of peanuts.

Identification of Pseudo-R genes in Vitis vinifera and characterization of their role as immunomodulators in host-pathogen interactions

INTRODUCTION

Duplication events are fundamental to co-evolution in host-pathogen interactions. Pseudogenes (Ψs) are dysfunctional paralogs of functional genes and resistance genes (Rs) in plants are the key to disarming pathogenic invasions. Thus, deciphering the roles of pseudo-R genes in plant defense is momentous.

OBJECTIVES

This study aimed to functionally characterize diverse roles of the resistance Ψs as novel gene footprints and as significant gene regulators in the grapevine genome.

METHODS

PlantPseudo pipeline and HMM-profiling identified whole-genome duplication-derived (WGD) Ψs associated with resistance genes (Ψ-Rs). Further, novel antifungal and antimicrobial peptides were characterized for fungal associations using protein-protein docking with Erysiphe necator proteins. miRNA and tasiRNA target sites and transcription factor (TF) binding sites were predicted in Ψ-Rs. Finally, differential co-expression patterns in Ψ-Rs-lncRNAs-coding genes were identified using the UPGMA method.

RESULTS

2,746 Ψ-Rs were identified from 31,032 WGD Ψs in the genome of grapevine. 69-antimicrobial and 81-antifungal novel peptides were generated from Ψ-Rs. The putative genic potential was predicted for five novel antifungal peptides which were further characterized by docking against E. necator proteins. 395 out of 527 resistance loci-specific Ψ-Rs were acting as parental gene mimics. Further, to explore the diverse roles of Ψ-Rs in plant-defense, we identified 37,026 TF-binding sites, 208 miRNA, and 99 tasiRNA targeting sites on these Ψ-Rs. 194 Ψ-Rs were exhibiting tissue-specific expression patterns. The co-expression network analysis between Ψs-lncRNA-genes revealed six out of 79 pathogen-responsive Ψ-Rs as significant during pathogen invasion.

CONCLUSIONS

Our study provides pathogen responsive Ψ-Rs integral for pathogen invasion, which will offer a useful resource for future experimental validations. In addition, our findings on novel peptide generations from Ψ-Rs offer valuable insights which can serve as a useful resource for predicting novel genes with the futuristic potential of being investigated for their bioactivities in the plant system.

Transcriptome Meta-Analysis Identifies Candidate Hub Genes and Pathways of Pathogen Stress Responses in Arabidopsis thaliana

Following a pathogen attack, plants defend themselves using multiple defense mechanisms to prevent infections. We used a meta-analysis and systems-biology analysis to search for general molecular plant defense responses from transcriptomic data reported from different pathogen attacks in Arabidopsis thaliana. Data from seven studies were subjected to meta-analysis, which revealed a total of 3694 differentially expressed genes (DEGs), where both healthy and infected plants were considered. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis further suggested that the DEGs were involved in several biosynthetic metabolic pathways, including those responsible for the biosynthesis of secondary metabolites and pathways central to photosynthesis and plant-pathogen interactions. Using network analysis, we highlight the importance of WRKY40, WRKY46 and STZ, and suggest that they serve as major points in protein-protein interactions. This is especially true regarding networks of composite-metabolic responses by pathogens. In summary, this research provides a new approach that illuminates how different mechanisms of transcriptome responses can be activated in plants under pathogen infection and indicates that common genes vary in their ability to regulate plant responses to the pathogens studied herein.

Chromosome-level genome assembly of Mentha longifolia L. reveals gene organization underlying disease resistance and essential oil traits

Mentha longifolia (L.) Huds., a wild, diploid mint species, has been developed as a model for mint genetic and genomic research to aid breeding efforts that target Verticillium wilt disease resistance and essential oil monoterpene composition. Here, we present a near-complete, chromosome-scale mint genome assembly for M. longifolia USDA accession CMEN 585. This new assembly is an update of a previously published genome draft, with dramatic improvements. A total of 42,107 protein-coding genes were annotated and placed on 12 chromosomal scaffolds. One hundred fifty-three genes contained conserved sequence domains consistent with nucleotide binding site-leucine-rich-repeat plant disease resistance genes. Homologs of genes implicated in Verticillium wilt resistance in other plant species were also identified. Multiple paralogs of genes putatively involved in p-menthane monoterpenoid biosynthesis were identified and several cases of gene clustering documented. Heterologous expression of candidate genes, purification of recombinant target proteins, and subsequent enzyme assays allowed us to identify the genes underlying the pathway that leads to the most abundant monoterpenoid volatiles. The bioinformatic and functional analyses presented here are laying the groundwork for using marker-assisted selection in improving disease resistance and essential oil traits in mints.

Bioinformatics approaches and applications in plant biotechnology

BACKGROUND

In recent years, major advance in molecular biology and genomic technologies have led to an exponential growth in biological information. As the deluge of genomic information, there is a parallel growth in the demands of tools in the storage and management of data, and the development of software for analysis, visualization, modelling, and prediction of large data set.

MAIN BODY

Particularly in plant biotechnology, the amount of information has multiplied exponentially with a large number of databases available from many individual plant species. Efficient bioinformatics tools and methodologies are also developed to allow rapid genome sequence and the study of plant genome in the 'omics' approach. This review focuses on the various bioinformatic applications in plant biotechnology, and their advantages in improving the outcome in agriculture. The challenges or limitations faced in plant biotechnology in the aspect of bioinformatics approach that explained the low progression in plant genomics than in animal genomics are also reviewed and assessed.

CONCLUSION

There is a critical need for effective bioinformatic tools, which are able to provide longer reads with unbiased coverage in order to overcome the complexity of the plant's genome. The advancement in bioinformatics is not only beneficial to the field of plant biotechnology and agriculture sectors, but will also contribute enormously to the future of humanity.

Genome-Wide Identification, Characterization, and Comparative Analysis of NLR Resistance Genes in Coffea spp

The largest family of disease resistance genes in plants are nucleotide-binding site leucine-rich repeat genes (NLRs). The products of these genes are responsible for recognizing avirulence proteins (Avr) of phytopathogens and triggering specific defense responses. Identifying NLRs in plant genomes with standard gene annotation software is challenging due to their multidomain nature, sequence diversity, and clustered genomic distribution. We present the results of a genome-wide scan and comparative analysis of NLR loci in three coffee species (Coffea canephora, Coffea eugenioides and their interspecific hybrid Coffea arabica). A total of 1311 non-redundant NLR loci were identified in C. arabica, 927 in C. canephora, and 1079 in C. eugenioides, of which 809, 562, and 695 are complete loci, respectively. The NLR-Annotator tool used in this study showed extremely high sensitivities and specificities (over 99%) and increased the detection of putative NLRs in the reference coffee genomes. The NLRs loci in coffee are distributed among all chromosomes and are organized mostly in clusters. The C. arabica genome presented a smaller number of NLR loci when compared to the sum of the parental genomes (C. canephora, and C. eugenioides). There are orthologous NLRs (orthogroups) shared between coffee, tomato, potato, and reference NLRs and those that are shared only among coffee species, which provides clues about the functionality and evolutionary history of these orthogroups. Phylogenetic analysis demonstrated orthologous NLRs shared between C. arabica and the parental genomes and those that were possibly lost. The NLR family members in coffee are subdivided into two main groups: TIR-NLR (TNL) and non-TNL. The non-TNLs seem to represent a repertoire of resistance genes that are important in coffee. These results will support functional studies and contribute to a more precise use of these genes for breeding disease-resistant coffee cultivars.

Underlying mechanism of accelerated cell death and multiple disease resistance in a maize lethal leaf spot 1 allele

Multiple disease resistance (MDR) in maize has attracted increasing attention. However, the interplay between cell death and metabolite changes and their contributions to MDR remains elusive in maize. In this study, we identified a mutant named as lesion mimic 30 (les30) that showed 'suicidal' lesion formation in the absence of disease and had enhanced resistance to the fungal pathogen Curvularia lunata. Using map-based cloning, we identified the causal gene encoding pheophorbide a oxidase (PAO), which is known to be involved in chlorophyll degradation and MDR, and is encoded by LETHAL LEAF SPOT1 (LLS1). LLS1 was found to be induced by both biotic and abiotic stresses. Transcriptomics analysis showed that genes involved in defense responses and secondary metabolite biosynthesis were mildly activated in leaves of the les30 mutant without lesions, whilst they were strongly activated in leaves with lesions. In addition, in les30 leaves with lesions, there was overaccumulation of defense-associated phytohormones including jasmonic acid and salicylic acid, and of phytoalexins including phenylpropanoids, lignin, and flavonoids, suggesting that their biosynthesis was activated in a lesion-dependent manner. Taken together, our study implies the existence of an interactive amplification loop of interrupted chlorophyll degradation, cell death, expression of defense-related genes, and metabolite changes that results in suicidal lesion formation and MDR, and this has the potential to be exploited by genetic manipulation to improve maize disease resistance.

Time-course transcriptome analysis identifies rewiring patterns of transcriptional regulatory networks in rice under Rhizoctonia solani infection

Sheath Blight (SB) disease in rice is caused by the infection from the fungal pathogen Rhizoctonia solani (R. solani). SB is one of the most severe rice diseases that can cause up to 50% yield losses in rice. Naturally occurring rice varieties resistant to SB have not been reported yet. We have performed a Time-Series RNA-Seq analysis on a widely cultivated rice variety BPT-5204 for identifying transcriptome level response signatures during R. solani infection at 1st, 2nd and 5th day post infection (dpi). In total, 428, 3225 and 1225 genes were differentially expressed in the treated rice plants on 1, 2 and 5 dpi, respectively. GO and KEGG enrichment analysis identified significant processes and pathways differentially altered in the rice plants during the fungal infection. Machine learning and network based integrative approach was used to construct rice Transcriptional Regulatory Networks (TRNs) for the three time points. TRN analysis identified SUB1B, MYB30 and CCA1 as important regulatory hub transcription factors in rice during R. solani infection. Jasmonic acid, salicylic acid, ethylene biogenesis and signaling were induced on infection. SAR was up regulated, while photosynthesis and carbon fixation processes were significantly down regulated. Involvement of MAPK, CYPs, peroxidase, PAL, chitinase genes were also observed in response to the fungal infection. The integrative analysis identified seven putative SB resistance genes differentially regulated in rice during R. solani infection.

Genome evolution and diversity of wild and cultivated potatoes

Potato (Solanum tuberosum L.) is the world's most important non-cereal food crop, and the vast majority of commercially grown cultivars are highly heterozygous tetraploids. Advances in diploid hybrid breeding based on true seeds have the potential to revolutionize future potato breeding and production^1-4. So far, relatively few studies have examined the genome evolution and diversity of wild and cultivated landrace potatoes, which limits the application of their diversity in potato breeding. Here we assemble 44 high-quality diploid potato genomes from 24 wild and 20 cultivated accessions that are representative of Solanum section Petota, the tuber-bearing clade, as well as 2 genomes from the neighbouring section, Etuberosum. Extensive discordance of phylogenomic relationships suggests the complexity of potato evolution. We find that the potato genome substantially expanded its repertoire of disease-resistance genes when compared with closely related seed-propagated solanaceous crops, indicative of the effect of tuber-based propagation strategies on the evolution of the potato genome. We discover a transcription factor that determines tuber identity and interacts with the mobile tuberization inductive signal SP6A. We also identify 561,433 high-confidence structural variants and construct a map of large inversions, which provides insights for improving inbred lines and precluding potential linkage drag, as exemplified by a 5.8-Mb inversion that is associated with carotenoid content in tubers. This study will accelerate hybrid potato breeding and enrich our understanding of the evolution and biology of potato as a global staple food crop.

Prediction of Plant Resistance Proteins Based on Pairwise Energy Content and Stacking Framework

Plant resistance proteins (R proteins) recognize effector proteins secreted by pathogenic microorganisms and trigger an immune response against pathogenic microbial infestation. Accurate identification of plant R proteins is an important research topic in plant pathology. Plant R protein prediction has achieved many research results. Recently, some machine learning-based methods have emerged to identify plant R proteins. Still, most of them only rely on protein sequence features, which ignore inter-amino acid features, thus limiting the further improvement of plant R protein prediction performance. In this manuscript, we propose a method called StackRPred to predict plant R proteins. Specifically, the StackRPred first obtains plant R protein feature information from the pairwise energy content of residues; then, the obtained feature information is fed into the stacking framework for training to construct a prediction model for plant R proteins. The results of both the five-fold cross-validation and independent test validation show that our proposed method outperforms other state-of-the-art methods, indicating that StackRPred is an effective tool for predicting plant R proteins. It is expected to bring some favorable contribution to the study of plant R proteins.

Construction and characterization of a de novo draft genome of garden cress (Lepidium sativum L.)

Garden cress (Lepidium sativum L.) is a Brassicaceae crop recognized as a healthy vegetable and a medicinal plant. Lepidium is one of the largest genera in Brassicaceae, yet, the genus has not been a focus of extensive genomic research. In the present work, garden cress genome was sequenced using the long read high-fidelity sequencing technology. A de novo, draft genome assembly that spans 336.5 Mb was produced, corresponding to 88.6% of the estimated genome size and representing 90% of the evolutionarily expected orthologous gene content. Protein coding gene content was structurally predicted and functionally annotated, resulting in the identification of 25,668 putative genes. A total of 599 candidate disease resistance genes were identified by predicting resistance gene domains in gene structures, and 37 genes were detected as orthologs of heavy metal associated protein coding genes. In addition, 4289 genes were assigned as "transcription factor coding." Six different machine learning algorithms were trained and tested for their performance in classifying miRNA coding genomic sequences. Logistic regression proved the best performing trained algorithm, thus utilized for pre-miRNA coding loci identification in the assembly. Repetitive DNA analysis involved the characterization of transposable element and microsatellite contents. L. sativum chloroplast genome was also assembled and functionally annotated. Data produced in the present work is expected to constitute a foundation for genomic research in garden cress and contribute to genomics-assisted crop improvement and genome evolution studies in the Brassicaceae family.

A successful defense of the narrow-leafed lupin against anthracnose involves quick and orchestrated reprogramming of oxidation-reduction, photosynthesis and pathogenesis-related genes

Narrow-leafed lupin (NLL, Lupinus angustifolius L.) is a legume plant cultivated for grain production and soil improvement. Worldwide expansion of NLL as a crop attracted various pathogenic fungi, including Colletotrichum lupini causing a devastating disease, anthracnose. Two alleles conferring improved resistance, Lanr1 and AnMan, were exploited in NLL breeding, however, underlying molecular mechanisms remained unknown. In this study, European NLL germplasm was screened with Lanr1 and AnMan markers. Inoculation tests in controlled environment confirmed effectiveness of both resistance donors. Representative resistant and susceptible lines were subjected to differential gene expression profiling. Resistance to anthracnose was associated with overrepresentation of "GO:0006952 defense response", "GO:0055114 oxidation-reduction process" and "GO:0015979 photosynthesis" gene ontology terms. Moreover, the Lanr1 (83A:476) line revealed massive transcriptomic reprogramming quickly after inoculation, whereas other lines showed such a response delayed by about 42 h. Defense response was associated with upregulation of TIR-NBS, CC-NBS-LRR and NBS-LRR genes, pathogenesis-related 10 proteins, lipid transfer proteins, glucan endo-1,3-beta-glucosidases, glycine-rich cell wall proteins and genes from reactive oxygen species pathway. Early response of 83A:476, including orchestrated downregulation of photosynthesis-related genes, coincided with the successful defense during fungus biotrophic growth phase, indicating effector-triggered immunity. Mandelup response was delayed and resembled general horizontal resistance.

Transcriptomic and genomic analysis provides new insights in molecular and genetic processes involved in zucchini ZYMV tolerance

BACKGROUND

Cucurbita pepo is highly susceptible to Zucchini yellow mosaic virus (ZYMV) and the resistance found in several wild species cannot be considered as complete or broad-spectrum resistance. In this study, a source of tolerance introgressed in C. pepo (381e) from C. moschata, in True French (TF) background, was investigated 12 days post-inoculation (DPI) at transcriptomic and genomic levels.

RESULTS

The comparative RNA-sequencing (RNA-Seq) of TF (susceptible to ZYMV) and 381e (tolerant to ZYMV) allowed the evaluation of about 33,000 expressed transcripts and the identification of 146 differentially expressed genes (DEGs) in 381e, mainly involved in photosynthesis, transcription, cytoskeleton organization and callose synthesis. By contrast, the susceptible cultivar TF triggered oxidative processes related to response to biotic stimulus and activated key regulators of plant virus intercellular movement. In addition, the discovery of variants located in transcripts allowed the identification of two chromosome regions rich in Single Nucleotide Polymorphisms (SNPs), putatively introgressed from C. moschata, containing genes exclusively expressed in 381e.

CONCLUSION

381e transcriptome analysis confirmed a global improvement of plant fitness by reducing the virus titer and movement. Furthermore, genes implicated in ZYMV tolerance in C. moschata introgressed regions were detected. Our work provides new insight into the plant virus recovery process and a better understanding of the molecular basis of 381e tolerance.

Transitions of foliar mycobiota community and transcriptome in response to pathogenic conifer needle interactions

Profiling the host–mycobiota interactions in healthy vs. diseased forest ecosystems helps understand the dynamics of understudied yet increasingly important threats to forest health that are emerging due to climate change. We analyzed the structural and functional changes of the mycobiota and the responses of Pinus contorta in the Lophodermella needle cast pathosystem through metabarcoding and metatranscriptomics. When needles transitioned from asymptomatic to symptomatic, dysbiosis of the mycobiota occurred, but with an enrichment of Lophodermella pathogens. Many pathogenicity-related genes were highly expressed by the mycobiota at the necrotrophic phase, showing an active pathogen response that are absent in asymptomatic needles. This study also revealed that Lophodermella spp. are members of a healthy needle mycobiota that have latent lifestyles suggesting that other pine needle pathogens may have similar biology. Interestingly, Pinus contorta upregulated defense genes in healthy needles, indicating response to fungal recognition, while a variety of biotic and abiotic stresses genes were activated in diseased needles. Further investigation to elucidate the possible antagonistic interplay of other biotic members leading to disease progression and/or suppression is warranted. This study provides insights into microbial interactions in non-model pathosystems and contributes to the development of new forest management strategies against emerging latent pathogens.

RFPDR: a random forest approach for plant disease resistance protein prediction

Background

Plant innate immunity relies on a broad repertoire of receptor proteins that can detect pathogens and trigger an effective defense response. Bioinformatic tools based on conserved domain and sequence similarity are within the most popular strategies for protein identification and characterization. However, the multi-domain nature, high sequence diversity and complex evolutionary history of disease resistance (DR) proteins make their prediction a real challenge. Here we present RFPDR, which pioneers the application of Random Forest (RF) for Plant DR protein prediction.

Methods

A recently published collection of experimentally validated DR proteins was used as a positive dataset, while 10x10 nested datasets, ranging from 400-4,000 non-DR proteins, were used as negative datasets. A total of 9,631 features were extracted from each protein sequence, and included in a full dimension (FD) RFPDR model. Sequence selection was performed, to generate a reduced-dimension (RD) RFPDR model. Model performances were evaluated using an 80/20 (training/testing) partition, with 10-cross fold validation, and compared to baseline, sequence-based and state-of-the-art strategies. To gain some insights into the underlying biology, the most discriminatory sequence-based features in the RF classifier were identified.

Results and Discussion

RD-RFPDR showed to be sensitive (86.4 ± 4.0%) and specific (96.9 ± 1.5%) for identifying DR proteins, while robust to data imbalance. Its high performance and robustness, added to the fact that RD-RFPDR provides valuable information related to DR proteins underlying properties, make RD-RFPDR an interesting approach for DR protein prediction, complementing the state-of-the-art strategies.

Genomic Analysis Based on Chromosome-Level Genome Assembly Reveals an Expansion of Terpene Biosynthesis of Azadirachta indica

Azadirachta indica (neem), an evergreen tree of the Meliaceae family, is a source of the potent biopesticide azadirachtin. The lack of a chromosome-level assembly impedes an in-depth understanding of its genome architecture and the comparative genomic analysis of A. indica. Here, a high-quality genome assembly of A. indica was constructed using a combination of data from Illumina, PacBio, and Hi-C technology, which is the first chromosome-scale genome assembly of A. indica. Based on the length of our assembly, the genome size of A. indica is estimated to be 281 Mb anchored to 14 chromosomes (contig N50 = 6 Mb and scaffold N50 = 19 Mb). The genome assembly contained 115 Mb repetitive elements and 25,767 protein-coding genes. Evolutional analysis revealed that A. indica didn't experience any whole-genome duplication (WGD) event after the core eudicot γ event, but some genes and genome segment might likely experienced recent duplications. The secondary metabolite clusters, TPS genes, and CYP genes were also identified. Comparative genomic analysis revealed that most of the A. indica-specific TPS genes and CYP genes were located on the terpene-related clusters on chromosome 13. It is suggested that chromosome 13 may play an important role in the specific terpene biosynthesis of A. indica. The gene duplication events may be responsible for the terpene biosynthesis expansion in A. indica. The genomic dataset and genomic analysis created for A. indica will shed light on terpene biosynthesis in A. indica and facilitate comparative genomic research of the family Meliaceae.

RNA-Seq of Cyst Nematode Infestation of Potato (Solanum tuberosum L.): A Comparative Transcriptome Analysis of Resistant and Susceptible Cultivars

Potato (Solanum tuberosum L.) is an important food crop worldwide, and potato cyst nematodes (PCNs) are among the most serious pests. The identification of disease resistance genes and molecular markers for PCN infestation can aid in crop improvement research programs against PCN infestation. In the present study, we used high-throughput RNA sequencing to investigate the comprehensive resistance mechanisms induced by PCN infestation in the resistant cultivar Kufri Swarna and the susceptible cultivar Kufri Jyoti. PCN infestation induced 791 differentially expressed genes in resistant cultivar Kufri Swarna, comprising 438 upregulated and 353 downregulated genes. In susceptible cultivar Kufri Jyoti, 2225 differentially expressed genes were induced, comprising 1247 upregulated and 978 downregulated genes. We identified several disease resistance genes (KIN) and transcription factors (WRKY, HMG, and MYB) that were upregulated in resistant Kufri Swarna. The differentially expressed genes from several enriched KEGG pathways, including MAPK signaling, contributed to the disease resistance in Kufri Swarna. Functional network analysis showed that several cell wall biogenesis genes were induced in Kufri Swarna in response to infestation. This is the first study to identify underlying resistance mechanisms against PCN and host interaction in Indian potato varieties.