Global transcriptome analyses reveal the molecular signatures in the early response of potato (Solanum tuberosum L.) to Phytophthora infestans, Ralstonia solanacearum, and Potato virus Y infection

StCoExpNet: a global co-expression network analysis facilitates identifying genes underlying agronomic traits in potatoes

KEY MESSAGE

We constructed a gene expression atlas and co-expression network for potatoes and identified several novel genes associated with various agronomic traits. This resource will accelerate potato genetics and genomics research. Potato (Solanum tuberosum L.) is the world's most crucial non-cereal food crop and ranks third in food production after wheat and rice. Despite the availability of several potato transcriptome datasets at public databases like NCBI SRA, an effort has yet to be put into developing a global transcriptome atlas and a co-expression network for potatoes. The objectives of our study were to construct a global expression atlas for potatoes using publicly available transcriptome datasets, identify housekeeping and tissue-specific genes, construct a global co-expression network and identify co-expression clusters, investigate the transcriptional complexity of genes involved in various essential biological processes related to agronomic traits, and provide a web server (StCoExpNet) to easily access the newly constructed expression atlas and co-expression network to investigate the expression and co-expression of genes of interest. In this study, we used data from 2299 publicly available potato transcriptome samples obtained from 15 different tissues to construct a global transcriptome atlas. We found that roughly 87% of the annotated genes exhibited detectable expression in at least one sample. Among these, we identified 281 genes with consistent and stable expression levels, indicating their role as housekeeping genes. Conversely, 308 genes exhibited marked tissue-specific expression patterns. We exemplarily linked some co-expression clusters to important agronomic traits of potatoes, such as self-incompatibility, anthocyanin biosynthesis, tuberization, and defense responses against multiple pathogens. The dataset compiled here constitutes a new resource (StCoExpNet), which can be accessed at https://stcoexpnet.julius-kuehn.de . This transcriptome atlas and the co-expression network will accelerate potato genetics and genomics research.

Abscisic-acid-responsive StlncRNA13558 induces StPRL expression to increase potato resistance to Phytophthora infestans infection

Late blight, caused by Phytophthora infestans, is one of the most serious diseases affecting potatoes (Solanum tuberosum L.). Long non-coding RNAs (lncRNAs) are transcripts with a length of more than 200 nucleotides that have no protein-coding potential. Few studies have been conducted on lncRNAs related to plant immune regulation in plants, and the molecular mechanisms involved in this regulation require further investigation. We identified and screened an lncRNA that specifically responds to P. infestans infection, namely, StlncRNA13558. P. infestans infection activates the abscisic acid (ABA) pathway, and ABA induces StlncRNA13558 to enhance potato resistance to P. infestans. StlncRNA13558 positively regulates the expression of its co-expressed PR-related gene StPRL. StPRL promotes the accumulation of reactive oxygen species and transmits a resistance response by affecting the salicylic acid hormone pathway, thereby enhancing potato resistance to P. infestans. In summary, we identified the potato late blight resistance lncRNA StlncRNA13558 and revealed its upstream and downstream regulatory relationship of StlncRNA13558. These results improve our understanding of plant-pathogen interactions' immune mechanism and elucidate the response mechanism of lncRNA-target genes regulating potato resistance to P. infestans infection.

Screening of wild species and transcriptome profiling to identify differentially regulated genes in response to late blight resistance in potato

Late blight (Phytophthora infestans) is a serious disease of potatoes. The aim of this study was to screen wild potato species and identify differentially expressed genes (DEGs) associated with late blight resistance. Wild potato species such as PIN45 (Solanum pinnatisectum), CPH62 (Solanum cardiophyllum), JAM07 (Solanum jamesii), MCD24 (Solanum microdontum), PLD47 (Solanum polyadenium), and cv. Kufri Bahar (control) were tested by artificial inoculation of P. infestans under controlled conditions. Transcriptomes of the leaf tissues (96 h post-inoculation) were sequenced using the Illumina platform. Statistically significant (p < 0.05) DEGs were analyzed in wild species by comparison with the control, and upregulated (>2 log₂ fold change, FC) and downregulated (<-2 log₂ FC) genes were identified. DEGs were functionally characterized with Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Selected genes were validated by real-time PCR analysis to confirm RNA-seq results. We identified some upregulated genes associated with late blight resistance in wild species such as cytochrome P450, proline-rich protein, MYB transcription factor MYB139, ankyrin repeat-containing protein, and LRR receptor-like serine/threonine-protein kinase in PIN45; glucosyltransferase, fructose-bisphosphate aldolase, and phytophthora-inhibited protease 1 in CPH62; steroid binding protein and cysteine proteinase 3 in JAM07; glycine-rich cell wall structural protein 1 and RING finger protein in MCD24; and cysteine proteinase 3 and major latex protein in PLD47. On the other hand, downregulated genes in these species were snakin-2 and WRKY transcription factor 3 in PIN45; lichenase and phenylalanine ammonia-lyase 1 in CPH62; metallothionein and LRR receptor-like serine/threonine-protein kinase in JAM07; UDP-glucoronosyl/UDP-glucosyl transferase family protein and steroid binding protein in MCD24; and cytoplasmic small heat shock protein class I and phosphatase PLD47. Our study identified highly resistant wild potato species and underlying genes such as disease resistance, stress response, phytohormones, and transcription factors (e.g., MYB, WRKY, AP2/ERF, and AN1) associated with late blight resistance in wild potato species.

Identification of a Comprehensive Gene Co-Expression Network Associated with Autotetraploid Potato (Solanum tuberosum L.) Development Using WGCNA Analysis

The formation and development of potato tissues and organs is a complex process regulated by a variety of genes and environmental factors. The regulatory mechanisms underlying the growth and development are still unclear. In this work, we aimed to explore the changes in gene expression patterns and genetic characteristics of potato tissues throughout different developmental stages. To achieve this, we used autotetraploid potato JC14 as an experimental subject to analyze the transcriptome of the root, stem, and leaf at the seedling, tuber formation, and tuber expansion stages. The results revealed thousands of differentially expressed genes, predominantly involved in defense response and carbohydrate metabolism according to KEGG pathway enrichment analysis. Weighted gene co-expression network analysis (WGCNA) revealed a total of 12 co-expressed gene modules, with 4 modules showing the highest correlation with potato stem development. By calculating the connectivity of genes within the module, hub genes were identified, and functional annotations were subsequently performed. A total of 40 hub genes from the four modules were identified, and their functions were found to be related to carbohydrate metabolism, defense response, and transcription factors. These findings provide important insights for further understanding of the molecular regulation and genetic mechanisms involved in potato tissue development.

Global transcriptome and targeted metabolite analyses of roots reveal different defence mechanisms against Ralstonia solanacearum infection in two resistant potato cultivars

Ralstonia solanacearum (Rs), the causal agent of bacterial wilt disease in an unusually wide range of host plants, including potato (Solanum tuberosum), is one of the most destructive phytopathogens that seriously reduces crop yields worldwide. Identification of defence mechanisms underlying bacterial wilt resistance is a prerequisite for biotechnological approaches to resistance breeding. Resistance to Rs has been reported only in a few potato landraces and cultivars. Our in vitro inoculation bioassays confirmed that the cultivars 'Calalo Gaspar' (CG) and 'Cruza 148' (CR) are resistant to Rs infection. Comparative transcriptome analyses of CG and CR roots, as well as of the roots of an Rs-susceptible cultivar, 'Désirée' (DES), were carried out two days after Rs infection, in parallel with their respective noninfected controls. In CR and DES, the upregulation of chitin interactions and cell wall-related genes was detected. The phenylpropanoid biosynthesis and glutathione metabolism pathways were induced only in CR, as confirmed by high levels of lignification over the whole stele in CR roots six days after Rs infection. At the same time, Rs infection greatly increased the concentrations of chlorogenic acid and quercetin derivatives in CG roots as it was detected using ultra-performance liquid chromatography - tandem mass spectrometry. Characteristic increases in the expression of MAP kinase signalling pathway genes and in the concentrations of jasmonic, salicylic, abscisic and indoleacetic acid were measured in DES roots. These results indicate different Rs defence mechanisms in the two resistant potato cultivars and a different response to Rs infection in the susceptible cultivar.

Dual RNA-seq Reveals the Global Transcriptome Dynamics of Ralstonia solanacearum and Pepper (Capsicum annuum) Hypocotyls During Bacterial Wilt Pathogenesis

Bacterial wilt, caused by Ralstonia solanacearum, is a serious disease in pepper. However, the interaction between the pathogen and pepper remains largely unknown. This study aimed to gain insights into determinants of pepper susceptibility and R. solanacearum pathogenesis. We assembled the complete genome of R. solanacearum strain Rs-SY1 and identified 5,106 predicted genes, including 84 type III effectors (T3E). RNA-seq was used to identify differentially expressed genes (DEGs) in susceptible pepper CM334 at 1 and 5 days postinoculation (dpi) with R. solanacearum. Dual RNA-seq was used to simultaneously capture transcriptome changes in the host and pathogen at 3 and 7 dpi. A total of 1,400, 3,335, 2,878, and 4,484 DEGs of pepper (PDEGs) were identified in the CM334 hypocotyls at 1, 3, 5, and 7 dpi, respectively. Functional enrichment of the PDEGs suggests that inducing ethylene production, suppression of photosynthesis, downregulation of polysaccharide metabolism, and weakening of cell wall defenses may contribute to successful infection by R. solanacearum. When comparing in planta and nutrient agar growth of the R. solanacearum, 218 and 1,042 DEGs of R. solanacearum (RDEGs) were detected at 3 and 7 dpi, respectively. Additional analysis of the RDEGs suggested that enhanced starch and sucrose metabolism, and upregulation of virulence factors may promote R. solanacearum colonization. Strikingly, 26 R. solanacearum genes were found to have similar DEG patterns during a variety of host-R. solanacearum interactions. This study provides a foundation for a better understanding of the transcriptional changes during pepper-R. solanacearum interactions and will aid in the discovery of potential susceptibility and virulence factors.

Germplasm, Breeding, and Genomics in Potato Improvement of Biotic and Abiotic Stresses Tolerance

Potato is one of the most important food crops in the world. Late blight, viruses, soil and tuber-borne diseases, insect-pests mainly aphids, whiteflies, and potato tuber moths are the major biotic stresses affecting potato production. Potato is an irrigated and highly fertilizer-responsive crop, and therefore, heat, drought, and nutrient stresses are the key abiotic stresses. The genus Solanum is a reservoir of genetic diversity, however, a little fraction of total diversity has been utilized in potato breeding. The conventional breeding has contributed significantly to the development of potato varieties. In recent years, a tremendous progress has been achieved in the sequencing technologies from short-reads to long-reads sequence data, genomes of Solanum species (i.e., pan-genomics), bioinformatics and multi-omics platforms such as genomics, transcriptomics, proteomics, metabolomics, ionomics, and phenomics. As such, genome editing has been extensively explored as a next-generation breeding tool. With the available high-throughput genotyping facilities and tetraploid allele calling softwares, genomic selection would be a reality in potato in the near future. This mini-review covers an update on germplasm, breeding, and genomics in potato improvement for biotic and abiotic stress tolerance.

Weighted Gene Co-Expression Analysis Network-Based Analysis on the Candidate Pathways and Hub Genes in Eggplant Bacterial Wilt-Resistance: A Plant Research Study

Solanum melongena L. (eggplant) bacterial wilt is a severe soil borne disease. Here, this study aimed to explore the regulation mechanism of eggplant bacterial wilt-resistance by transcriptomics with weighted gene co-expression analysis network (WGCNA). The different expression genes (DEGs) of roots and stems were divided into 21 modules. The module of interest (root: indianred4, stem: coral3) with the highest correlation with the target traits was selected to elucidate resistance genes and pathways. The selected module of roots and stems co-enriched the pathways of MAPK signalling pathway, plant pathogen interaction, and glutathione metabolism. Each top 30 hub genes of the roots and stems co-enriched a large number of receptor kinase genes. A total of 14 interesting resistance-related genes were selected and verified with quantitative polymerase chain reaction (qPCR). The qPCR results were consistent with those of WGCNA. The hub gene of EGP00814 (namely SmRPP13L4) was further functionally verified; SmRPP13L4 positively regulated the resistance of eggplant to bacterial wilt by qPCR and virus-induced gene silencing (VIGS). Our study provides a reference for the interaction between eggplants and bacterial wilt and the breeding of broad-spectrum and specific eggplant varieties that are bacterial wilt-resistant.

Transcriptional Changes in Potato Sprouts upon Interaction with Rhizoctonia solani Indicate Pathogen-Induced Interference in the Defence Pathways of Potato

Rhizoctonia solani is the causer of black scurf disease on potatoes and is responsible for high economical losses in global agriculture. In order to increase the limited knowledge of the plants' molecular response to this pathogen, we inoculated potatoes with R. solani AG3-PT isolate Ben3 and carried out RNA sequencing with total RNA extracted from potato sprouts at three and eight days post inoculation (dpi). In this dual RNA-sequencing experiment, the necrotrophic lifestyle of R. solani AG3-PT during early phases of interaction with its host has already been characterised. Here the potato plants' comprehensive transcriptional response to inoculation with R. solani AG3 was evaluated for the first time based on significantly different expressed plant genes extracted with DESeq analysis. Overall, 1640 genes were differentially expressed, comparing control (-Rs) and with R. solani AG3-PT isolate Ben3 inoculated plants (+Rs). Genes involved in the production of anti-fungal proteins and secondary metabolites with antifungal properties were significantly up regulated upon inoculation with R. solani. Gene ontology (GO) terms involved in the regulation of hormone levels (i.e., ethylene (ET) and jasmonic acid (JA) at 3 dpi and salicylic acid (SA) and JA response pathways at 8 dpi) were significantly enriched. Contrastingly, the GO term "response to abiotic stimulus" was down regulated at both time points analysed. These results may support future breeding efforts toward the development of cultivars with higher resistance level to black scurf disease or the development of new control strategies.

Transcriptome Analysis of Resistance to Fusarium Wilt in Mung Bean (Vigna radiata L.)

Fusarium wilt is a destructive soil-borne disease that threatens the production of mung bean. Mung bean lines Zheng8-4 and Zheng8-20 show high resistance and high susceptibility to Fusarium wilt, respectively. Transcriptome analysis was carried out to identify candidate genes involved in Fusarium wilt resistance using Zheng8-4 and Zheng8-20 at 0, 0.5, 1, 2, and 4 days post inoculation (dpi). Differential expression analysis showed that 3,254 genes responded to pathogen infection and were differentially expressed in the resistant and susceptible lines. Weighted gene co-expression network analysis (WGCNA) was also performed to identify five modules highly correlated with Fusarium wilt resistance, in which 453 differentially expressed genes (DEGs) were considered likely to be involved in Fusarium wilt resistance. Among these DEGs, we found 24 genes encoding resistance (R) proteins, 22 encoding protein kinases, 20 belonging to transcription factor families, 34 encoding proteins with oxidoreductase activity, 17 involved in stimulation/stress responses, and 54 annotated to pathogen resistance-related pathways. Finally, 27 annotated genes were further selected as candidate genes of Fusarium wilt resistance in mung bean. This study identifies novel potential resistance-related genes against Fusarium wilt and provides a theoretical basis for further investigation of Fusarium wilt resistance in mung bean breeding.