Front. Plant Sci. Frontiers in Plant Science Front. Plant Sci. 1664-462X Frontiers Media S.A. 10.3389/fpls.2026.1745952 Original Research Transcriptome profiling reveals tissue-wide gene expression in chili pepper (Capsicum annuum L.) under infection by Phytophthora capsici LozadaDennis Nicuh 1 * Conceptualization Formal analysis Funding acquisition Project administration Supervision Writing – original draft Writing – review & editing Cabrales-ArellanoPatricia 2 Conceptualization Formal analysis Funding acquisition Methodology Project administration Writing – review & editing MatresJerlie Mhay 1 Investigation Writing – review & editing AngelesJorge Gil C. 3 Conceptualization Formal analysis Funding acquisition Writing – review & editing Velazquez-MartinezVictor 4 Methodology Writing – review & editing KaurNavdeep 5 Investigation Methodology Writing – review & editing DelgadoEfren 6 Funding acquisition Resources Writing – review & editing SanogoSoum 7 Funding acquisition Resources Writing – review & editing Department of Plant and Environmental Sciences, New Mexico State University, Las Cruces, NM, United States Department of Biology, Eastern New Mexico University, Portales, NM, United States Philippine Genome Center - Program for Agriculture, Livestock, Fisheries and Forestry, University of the Philippines Los Baños, Laguna, Los Baños, Philippines Facultad de Ingenieria Mecanica Electrica, Universidad Veracruzana, Poza Rica, Veracruz, Mexico Department of Horticultural Sciences, Texas A&M University, College Station, TX, United States Department of Family and Consumer Sciences, New Mexico State University, Las Cruces, NM, United States Department of Entomology, Plant Pathology and Weed Science, New Mexico State University, Las Cruces, NM, United States Correspondence: Dennis Nicuh Lozada, dlozada@nmsu.edu 26 02 2026 2026 17 1745952 13 11 2025 22 01 2026 07 01 2026 Copyright © 2026 Lozada, Cabrales-Arellano, Matres, Angeles, Velazquez-Martinez, Kaur, Delgado and Sanogo. 2026 Lozada, Cabrales-Arellano, Matres, Angeles, Velazquez-Martinez, Kaur, Delgado and Sanogo https://creativecommons.org/licenses/by/4.0/ This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Phytophthora blight, manifested by root, stem, and fruit rot, caused by the oomycete Phytophthora capsici, is an important disease affecting Chili pepper production globally. RNA sequencing (RNA-seq) was performed to identify differentially expressed genes (DEGs) and shared genetic resistance mechanisms across different tissues upon infection by the pathogen. RNA-seq revealed the dynamic transcriptome of leaf, stem, and root tissues from resistant (R; CM-334) and susceptible (S; Early Jalapeño) varieties under different times of infection by P. capsici. There were 149,531 differentially expressed genes (DEGs) from 39 different R vs. S, time vs. time, and tissue vs. tissue comparisons. A total of 75,520 DEGs (51%) showed higher expression, whereas 74,011 (49%) demonstrated lower expression across all tissues and times of post-inoculation. The total number of DEGs with higher expression for the different tissue samples decreased across times of post-inoculation, where the 72h post-inoculation showed the least number of genes. The roots generally showed a higher number of DEGs compared to the stems and the leaves. Network analyses of DEGs demonstrated that genes with functions related to defense response to fungal infection were also involved with carbohydrate metabolism and ADP binding. Genes related to immune response to fungal infection and amino acid metabolism (e.g., homoserine kinase activity) showed higher gene expression across all times of infection and tissue samples. Chili pepper transcriptome under P. capsici infection provides evidence of shared gene expression across multiple tissues which can be leveraged for breeding and selection for broad-spectrum resistance in current Capsicum germplasm.

 candidate genes disease resistance functional genomics gene expression phytophthora root rot RNA-sequencing The author(s) declared that financial support was received for this work and/or its publication. This study was funded by the New Mexico Department of Agriculture Specialty Crop Block Grant Program Award No. AM22SCBPNM1121-00, the USDA Specialty Crop Research Initiative Grant. No. 2023-51181-41318, and the New Mexico Chile Association. section-at-acceptance Functional and Applied Plant Genomics Introduction

For more than one hundred years, Phytophthora blight, caused by the oomycete Phytophthora capsici, has been one of the major diseases affecting Chili pepper (Capsicum spp.) production in many growing areas of the world (Quesada-Ocampo et al., 2023; Sanogo et al., 2023). First described by L. H. Leonian in the state of New Mexico, USA (Leonian, 1922), this disease continues to threaten chili pepper production resulting in losses of up to 100% in infected fields. P. capsici has been ranked as one of the top five plant-pathogenic oomycetes with scientific and economic importance, alongside P. infestans, Hyalopernospora arabidopsis, P. ramorum, and P. sojae (Kamoun et al., 2015). P. capsici has a broad and expanding host range including plant species in the Cucurbitaceae, Fabaceae, Malvaceae, and Solanaceae families (Erwin and Ribeiro, 1996; Lamour et al., 2012). Breeding for disease resistance remains the most economical approach to mitigate the effects of P. capsici in pepper and other host plants. Modern omics tools and integrated management strategies have the potential to render a better understanding of the genetic architecture of the disease enabling more informed and robust breeding and selection decisions for improving disease resistance (Lozada et al., 2022; Sanogo et al., 2023).

As a pathogenic oomycete, P. capsici has been known to infect various parts of the chili pepper plant, including the root and stem causing root rot and stem blight (Sy et al., 2005; Kang et al., 2022b). These disease syndromes are influenced by the host species, environmental conditions, and points of infection (Barchenger et al., 2018a). A landrace from Mexico, ‘Criollo de Morellos-334’ (CM-334), has broad-spectrum resistance to P. capsici, and has been regarded as a universally resistant host (Reyes-Tena et al., 2019). CM-334, nonetheless, has little to no market and/or agronomic value. Other sources of resistance including pasilla, piquin, cola de rata, and manzano types have been previously identified (Retes-Manjarrez et al., 2020). More recently, Tipo Ancho, Chilhuacle Orange, and Tipo Pasilla have been identified to possess broad-spectrum resistance to three different isolates of P. capsici (Kaur et al., 2024), providing other resistant sources to examine the Capsicum-P. capsici pathosystem using various genomewide approaches.

Functional genomics is the development and application of genome-wide experimental approaches to evaluate the function and expression of genes (Hieter and Boguski, 1997). One of its goals is to bridge the gap between genomic sequences (as anchor points) and function to render novel insights into the characteristics of biological systems (Hieter and Boguski, 1997; Werner, 2010). As an approach, RNA-seq has become a standard method for studying gene expression, specifically for evaluating transcript abundance and diversity (Griffith et al., 2015). Previous studies have demonstrated the potential of using RNA-seq approaches to explore the genes and gene systems associated with response against plant infection by P. capsici. Different transcription factor and protein families including protein kinases were observed in two contrasting chili pepper landraces (i.e., ‘GojamMecha_9086’ (resistant) and ‘Dabat_80045’ (susceptible)) exposed to P. capsici (Rabuma et al., 2022). In another study, candidate genes related to the modification of cell wall, phytoalexins, and phytohormones were identified to be associated with responses to P. capsici in the pepper line PI 201234 (Wang et al., 2015). The involvement of the phenylpropanoid biosynthesis pathway has been previously implicated in resistance against P. capsici in whole roots of chili pepper (Li et al., 2020) and black pepper (Piper spp.) (Hao et al., 2016). The higher expression of enzymes, including proteases (subtilisin-like protease and xylem cysteine proteinase 1), and various pathways such as Ca2+- and salicylic acid-mediated signaling, and flavonoid biosynthesis pathways were further observed in response to infection by P. capsici (Lei et al., 2023). A total of 14 putative effectors were also previously observed to be differentially expressed between various Capsicum-P. capsici interactions, where 12 genes were classified as RxLR (Arginine-any amino acid-Leucine-Arginine) and two as CRN (Crinkling and Necrosis) effectors (Maillot et al., 2022).

In the current study, we have identified differentially expressed genes in various tissues (roots, stems, and leaves) of resistant and susceptible pepper varieties under infection by P. capsici using an Illumina-based RNA-seq approach. The differential expression of genes and gene systems across different plant tissues contribute to the overall genetic resistance and/or susceptibility to infection of Capsicum by P. capsici. Results from gene expression profiling could be utilized for omics-assisted breeding for disease resistance in current chili pepper germplasm.

 Materials and methods Plant material and inoculation with <italic>P. capsici</italic>

Two varieties, namely, CM-334 (resistant) and Early Jalapeño (susceptible) (C. annuum L.) were used in the current study. These varieties were used for transcriptome analyses as they were previously used as parents to develop the New Mexico recombinant inbred lines (NMRIL) population (Sy et al., 2005). The NMRIL has been extensively used to examine the genetics of P. capsici resistance in C. annuum L. (Monroy-Barbosa and Bosland, 2008; Sy et al., 2008; Jiang et al., 2015; Barchenger et al., 2018b; Lozada et al., 2021a). Seeds were planted in three replications and maintained under standard greenhouse conditions for growing chili pepper at the Fabian Garcia Science Center Greenhouse, New Mexico State University, Las Cruces, NM (Sharma et al., 2017) in March 2022. The greenhouse is comprised of aluminum frames, double layer of polycarbonate sheets for insulation, evaporative cooler, and heaters. The temperature inside the greenhouse is controlled by an automatic control system, with average temperatures of 27.5°C (daytime; mean 12.5 hours) and 21.6°C (nighttime; mean 11.5 hours). At the 4–8 leaf stage, the soil was inoculated with ~10,000 zoospores of 6347, a virulent isolate of P. capsici collected from a commercial farm in New Mexico (Jiang et al., 2015), following the methods previously described (Kaur et al., 2024; Lozada et al., 2021a). At times 0-, 24- and 72- hours post infection (hpi), leaf, stem, and root tissues were sampled and stored in DNA/RNA shield buffer for RNA extraction (Zymo Research, Irvine, CA, USA) in April 2022. In total, there were 18 variety-tissue-timepoint combinations [two varieties (resistant vs. susceptible) x three tissues (leaves, roots, and stems) x three timepoints (0, 24, and 72)] in three biological replicates (total n = 54 samples).

 RNA isolation, library construction, and sequencing

Total RNA was extracted from leaf, roots, and stems of CM-334 and Early Jalapeño at different time points using the Trizol reagent (Invitrogen) following the manufacturer’s instructions. DNA contamination was removed using RNA clean and concentrator (Zymo Research, CA, USA). RNA was quantified using Qubit® RNA Assay Kit (Life Technologies, MA, USA). Total RNA (500–1000 ng) was used for library preparation using poly-A selection method of KAPA mRNA HyperPrep Kits (Roche, CA, USA) following the manufacturer’s protocol. The final concentration of all cDNA libraries was measured using the Qubit® dsDNA HS Assay Kit (Life Technologies). Agilent 2100 Bioanalyzer (Agilent Technologies, CA, USA) was used to determine the average size of the libraries. Paired-end (150-bp) sequencing was performed for all libraries pooled in equimolar ratios of 0.6nM using the Illumina NovaSeq 6000® platform (Illumina, CA, USA). The read files for each sample were analyzed using Cutadapt v4.4 (Martin, 2011) to identify and remove Illumina sequencing adapters. Trimmomatic v0.40 (Bolger et al., 2014) was used to eliminate poor quality reads, or low quality 5’ or 3’ bases or fragments implementing a sliding window approach, where a window of 10 bp was moved from the 5’ to the 3’ end of the read one base at a time. Reads or bases with the following qualities: (a) reads with average quality score < 20; (b) reads < 50 bp in length; and (c) bases with quality scores ≤ 3 at the 5’ or 3’ bases were excluded. Applying these criteria, two samples (viz., S_leaf_24_2 and S_root_72_3) were excluded, resulting in n = 52 biological samples for transcriptome analyses. Supplementary Table 1 summarizes the software tools and packages used for the analyses of RNA-seq data.

 Quality control and quantification of transcript number

Following the quality control filtering steps, the number of reads per transcript was quantified using Kallisto v0.46.1 (Bray et al., 2016), which implements a pseudoalignment approach to rapidly determine the compatibility of reads with transcripts. A list of k-mers (sequence fragments of length k) from the transcriptome were generated using Kallisto. Pseudoalignments were created and quantified through Kallisto by comparing the k-mers of each read to the de Bruijn graph model of the transcripts. The C. annuum L. reference genome UCD10Xv1.1 (RefSeq GCF_002878395.1; GenBank GCA_002878395.3 Annotation release 101; Hulse-Kemp et al., 2018) was used for the estimation of the number of transcripts. The reads per transcript counts reported by Kallisto were then used to generate reads per gene counts with the ‘tximport’ v1.28 package (Soneson et al., 2016) in R.

 Cluster analyses

Principal components analysis (PCA) and Principal coordinates analysis (PCoA) were implemented to visualize the level of similarity of the samples’ transcriptional profiles (i.e., the number of reads assigned to different genes in different samples). Prior to implementing the PCoA and PCA, the reads per gene counts for each sample were normalized using a variance-stabilizing transformation (VST) approach in DeSeq2 v.3.22, which normalizes the data in relation to the library size and removes the dependence of the variance on the mean, where genes with low mean counts show higher variance than genes with high mean counts (Love et al., 2014). The transformation measures and removes the experiment-wide trend of variance over mean. The PCA was performed using a singular value decomposition of the normalized counts per gene matrix (Venables and Ripley, 2013). The normalized counts were used to calculate the Euclidean distances between each pair of samples and the distance matrix was used as input for PCoA (Gower, 1966). The ellipses for groups with more than 3 samples in the PCoA plots show the 95% confidence level for a multivariate t-distribution (Fox and Weisberg, 2011; Friendly et al., 2013).

 Identification and functional annotation of differentially expressed genes

Genes that demonstrated significant expression changes across different groups of samples (pairwise comparisons/contrasts) were identified using DESeq2, where a Wald test with Benjamini-Hochberg correction for sampling was performed (Love et al., 2014). Genes with adjusted P ≤ 0.05 and log2 of fold change (log2FC) ≥ 1 and ≤ -1 were considered to have higher expression and lower expression, respectively. Functional annotation was performed using different information including (a) the protein product obtained from the general feature format (GFF) annotations for the corresponding reference genome and transcriptome in the NCBI RefSeq database; (b) KEGG pathways (Kanehisa and Goto, 2000); and (c) Gene ontology (GO) terms retrieved from the UniProt database using the ‘AnnotationHub’ package (Morgan et al., 2019) in R. Heatmaps showing expression changes (fold changes, FC) for the DEGs based on their corresponding log2FC values were obtained using the ‘pheatmap’ package (Kolde and Kolde, 2015) in R. Enrichment of biological functions among the DEGs were performed using a hypergeometric test and P-values were adjusted for multiple testing using the Benjamini-Hochberg correction.

 Results RNA quality and summary of RNA sequencing

The total RNA concentration for the samples (n = 54) ranged between 33.8 ng/ul and 5,747 ng/ul, with a mean value of 486.30 ng/ul. The average RNA integrity number (RIN) value for a random subset of the samples (n = 21) was 9.03, with a range between 7.1 and 9.8. A total of 592,493,398 reads from n = 54 samples were subjected to transcriptome profiling (Table 1). Out of this number, 447,084,240 (75.5%) were pseudoaligned using Kallisto. Non-duplicate (unique) reads comprised 55.2% of the total number of reads (327,117,203). The average number of reads for the 54 samples was 10,972,100 (Supplementary Table 2). Considering n = 52 samples, there were 562,528,911 total reads, of which 437,616,790 (77.8%) were pseudoaligned. The number of unique reads was 319,695,213 (56.8%), with an average number of 10,817,864 reads, for n = 52 samples.

Summary of RNA sequencing results.

Total number of samples Total number of reads Average number of reads Total number of pseudoaligned reads Pseudoaligned reads (%) Number of unique (non-duplicate) reads Unique reads (%)
54 592,493,398 10,972,100 447,084,240 75.5 327,117,203 55.2
52 562,528,911 10,817,864 437,616,790 77.8 319,695,213 56.8
Differential expression of genes

There were 149,531 differentially expressed genes (DEGs) (P ≤ 0.05 and log2 of fold change (log2FC) ≥ 1 and ≤ -1 for high gene expression and low gene expression, respectively) representing 39 different contrasts/comparisons from 52 biological samples subjected to RNA-seq analyses. Different versions of the same gene (i.e., isoforms) were counted separately. The Top 100 differentially expressed genes across the different tissues and hours post-inoculation are listed on Supplementary Tables 3-10. A total of 75,520 (51%) DEGs showed higher expression, whereas 74,011 (49%) DEGs showed lower expression in the R vs. S contrast. Excluding the transcripts in unplaced (unmapped) scaffolds, there were 134,753 DEGs, with 68,231 (51%) showing higher expression and 66,522 (49%) demonstrating lower expression. Except for the stems, the total number of DEGs with higher expression for the different tissue samples decreased across the different times of post-inoculation, where the 72-h post-inoculation showed the least number of genes. Comparing the resistant (R) vs. susceptible (S) varieties across different tissues, a total of 3,652, 3,037, and 2,069 differentially expressed genes (DEGs) were identified in the roots, stem, and leaves, respectively. In the roots, 1,535 genes showed higher expression, and 2,117 genes were observed to have lower expression; in the stem, 1,446 demonstrated higher expression and 1,591 genes showed lower expression; and in the leaves, 907 genes showed higher expression, and 1,162 genes demonstrated lower expression.

Comparing the R vs. S varieties across different times of post-inoculation, a total of 1,736 DEGs showed higher expression in the roots at 0h (Figure 1). The number of DEGs with higher expression decreased to 1,069 at 24h, and to 587 at 72h post-inoculation with P. capsici. A total of 922, 246, and 115 genes were unique to 0h, 24h, and 72h post-inoculation, respectively, whereas 389 genes were shared in all timepoint comparisons. In comparison to the root tissues, the total number of genes with higher expression was relatively less in the stems and the leaves (2,560 and 1,681 DEGs, respectively). For the stems, the number of distinct genes with higher expression was 725 at 0h, 125 at 24h, and 884 at 72h, whereas in the leaves, there were 950 (0h), 260 (24h), and 91 (72h) genes that showed higher expression. A total of 382 and 144 DEGs with higher expression were shared at all timepoint comparisons for the stem and leaves, respectively.

Number of differentially expressed genes across different times (hours, h) post-infection with P. capsici (0-, 24-, and 72-h for the resistant (R; CM-334) vs. susceptible (S; Early Jalapeno) contrast for root, stem, and leaf tissue samples. .

Six Venn diagrams representing gene expression data at three time points (Time 0, Time 24, Time 72) in root, stem, and leaf samples. The top row shows higher expression: root, stem, and leaf. The bottom row shows lower expression: root, stem, and leaf. Each diagram shows the number and percentage of genes expressed in each section.

At different times of infection in the R vs. S contrast, 2,668 genes showed lower gene expression in the roots, whereas 2,676, and 2,065 genes were observed to have lower expression in the stems and leaves, respectively. The number of genes with lower expression also decreased over time of post-inoculation for the root and leaf tissues, showing a similar trend with the observation for the genes with higher expression. In the root, 880 unique genes showed lower expression at 0h post-inoculation (hpi) whereas, 466 unique genes were observed to have lower expression at 24hpi. The number of unique DEGs with lower expression further decreased to 233 at 72hpi. A total of 465 DEGs in the root tissue were common across all timepoints. The number of DEGs with lower expression was higher in the stem than the root (2,676), with 926, 182, and 753 unique genes for the 0h, 24h, and 72h post-inoculation, respectively. A total of 2,065 genes showed lower expression in the leaf tissues, with 97 genes common across all the times of post-inoculation. Overall, there were 732 genes with lower expression unique for 0h; 647 for 24h; and 307 for 72h in the leaf tissue after inoculation with P. capsici.

 Principal components analysis and multidimensional scaling

Principal components analysis (PCA) was implemented to examine the relationships between the transcriptional profiles of the samples. Samples S_leaf_24_2 and S_root_72_3 (Susceptible_Leaf_ 24hpi_biological replicate 2 and Susceptible_Root_72hpi_biological replicate 3, respectively) and were identified as outliers by visual inspection of the PCA clustering and hence were excluded from further analyses. Clustering based on tissue type was observed when all samples were analyzed using PCA (Figure 2). The first and second principal components (PC1 and PC2) were associated with 35% and 18% of variation, respectively, for the root samples. Grouping based on resistance or susceptibility was most apparent in roots and stems, whereas no clear clustering was observed among the leaf samples. The transcriptional profiles for the varieties were more similar prior to inoculation with P. capsici (0h), relative to after inoculation at 24hpi and 72hpi. The transcriptome at 0hpi tends to form a tight cluster within the resistant and susceptible varieties across the different tissues, whereas there was no clear grouping that could be observed for RNA-seq profiles at 24hpi and 72hpi.

Principal components analysis (PCA) biplots of the transcriptomes of CM-334 (resistant) and Early Jalapeno (susceptible) across various times (t) of P. capsici infection (0, 24, 72 h) for (a) root, (b) stem, and (c) leaf tissue samples.

Three scatter plots showing principal component analysis of plant variants at different times. The first plot shows root data, the second stem data, and the third leaf data. Red dots represent resistant variants and blue dots represent susceptible variants. Time points are indicated by different shapes: circles for time zero, triangles for twenty-four hours, and squares for seventy-two hours. Ellipses encompass clusters of points. Percentages on each axis indicate the variance explained by principal components.

In addition to the PCA, the level of similarity between the transcriptional profiles between the samples was visualized using principal coordinates analysis (PCoA) (Supplementary Figure 1). Clustering based on the tissue source and not based on resistance or susceptibility of the varieties, was apparent when all samples were analyzed altogether using PCoA. When samples were grouped based on resistance or susceptibility, PCoA clustering was more pronounced in the root samples, where the transcripts of the resistant variety grouped together. Conversely, grouping based on resistance/susceptibility was less observable in the stem samples and was not apparent in the leaf samples. More robust clustering between samples at 0hpi was observed, relative to samples at 24hpi and 72hpi when PCoA was performed based on inoculation time.

 Differentially expressed genes with higher expression

The differentially expressed genes with higher expression for the R vs. S contrast belong to different chromosomes including chromosomes 1, 2, 5, 7, 11, and 12 (Table 2). In the roots, the mean log2FC value was highest for Zeatin O-xylosyltransferase with a UDP-glycosyltransferase activity in chromosome 5. This gene was identified on the Ext-Pc5.1 region for P. capsici resistance previously described by Du et al. (2021). Similarly, homoserine kinase also showed higher expression and was mapped within the Ext-Pc5.1 region. Genes with roles related to glutathione metabolic process, glutathione transferase activity, ATP-dependent DNA damage activity, also demonstrated higher expression. In the resistant variety (CM-334) across time vs. time comparisons, genes with biological functions associated with systemic acquired resistance (e.g., putative lipid-transfer protein DIR1), defense response to other organisms (e.g., pathogenesis-related protein R major form; pathogenesis-related protein P2; putative late blight resistance protein homolog R1A-3; disease resistance protein RGA2-like), and defense response to fungus (e.g., acidic endochitinase pcht28, Early nodulin-75, pathogenesis-related protein P2) showed consistent higher expression across multiple tissues (Table 3; Figure 3).

The top 20 differentially expressed genes with higher expression in roots and their functions.

Gene ID Chr. Position (start, end; in bp) Gene Function(s) Average log2FC
107870473 5 26240506, 26244754 Zeatin O-xylosyltransferase UDP-glycosyltransferase activity 10.91
107848341 11 32830222, 32831544 Glutathione S-transferase U17 Glutathione metabolic process; glutathione transferase activity 10.56
107877977 7 161977837, 162066908 DNA mismatch repair protein MSH7 Mismatched DNA binding; ATP binding; mismatch repair; ATP-dependent DNA damage activity 10.54
107870410 5 20998254, 20999886 Homoserine kinase Homoserine kinase activity; ATP binding; threonine metabolic process 10.17
107850467 12 8368711, 8374174 Beta-glucosidase 42 Carbohydrate metabolic process; beta-glucosidase activity 10.08
107877337 7 210507508, 210509477 Protein NODULATION SIGNALING PATHWAY 2 Sequence-specific DNA binding; DNA-binding transcription factor activity; regulation of DNA-templated transcription 9.34
107860513 2 150028093, 150032914 A-kinase anchor protein 7 DNA dealkylation involved in DNA repair; regulation of DNA-templated transcription 9.16
107875488 1 55362182, 55372791 COP9 signalosome complex subunit 2 Protein deneddylation 9.08
107878033 1 96977394, 96986051 Serine/threonine-protein Phosphatase PP2A catalytic subunit Peptidyl-serine dephosphorylation; mitotic cell cycle; protein serine/threonine phosphatase activity 8.82
107868168 4 23444894, 23457437 Vacuolar protein sorting-associated protein 8 homolog Protein binding; endosomal vesicle fusion; late endosome; HOPS complex; protein targeting to vacuole 8.59
107842555 9 65999409, 66002761 Putative glutathione peroxidase 5 Response to oxidative stress; glutathione peroxidase activity 8.31
107859557 2 133594838, 133611861 Pectinesterase Cell wall modification; pectinesterase inhibitor activity; pectinesterase activity 8.04
107872960 6 175798464, 175799168 Protein cornichon homolog 4 Vesicle-mediated transport 7.95
107868113 4 58273595, 58279483 Transcription factor bHLH62 Regulation of DNA-templated transcription; DNA-binding transcription factor activity 7.88
107848161 11 58156110, 58166643 Putative BOI-related E3 ubiquitin-protein ligase 3 Ubiquitin-protein transferase activity 7.58
107840746 9 217422587, 217424606 Putative late blight resistance protein homolog R1B-13 Protein folding; alpha-tubulin binding 7.53
107857159 1 193609050, 193609671 Putative lipid-transfer protein DIR1 Fatty acid binding; systemic acquired resistance 7.30
107871108 5 205442143, 205443599 U4/U6 small nuclear ribonucleoprotein Prp31 homolog U4/U6 x U5 tri-snRNP complex;U4 snRNP; precatalytic spliceosome; mRNA splicing, via spliceosome; spliceosomal tri-snRNP complex assembly 7.04
107863749 3 68261538, 68269405 Glycine cleavage system H protein Glycine decarboxylation via glycine cleavage system; protein lipoylation; glycine cleavage complex; DNA binding 6.73
124897050 3 29468788, 29477960 Cyanidin 3-O-glucoside 7-O-glucosyltransferase (acyl-glucose)-like Beta-glucosidase activity; carbohydrate metabolic process 6.71

Differential expressions based on three contrasts for the resistant (R) vs. susceptible (S) contrast (0h/0h; 24h/24h; 72h/72h post-inoculation). Genes in chromosome 5 (italicized and boldfaced) lie in the extended P. capsici genomic region (Ext-Pc5.1) previously identified by Du et al. (2021).

Significantly enriched (higher expression; P < 0.05) genes and their corresponding biological processes and functions in CM-334 (resistant variety (time vs. time comparison)) across various tissues. Functions in bold and italics are common in multiple tissues.

Biological processes/functions Gene Gene ratio Contrast(s) (time, t in hours post-inoculation) Tissue
Systemic acquired resistance Putative lipid-transfer protein DIR1 0.04 72/0 Root
Putative lipid-transfer protein DIR1 0.05 72/0; 24/0 Stem
Putative lipid-transfer protein Protein DIR1 0.05 72/24 Leaves
Plant-pathogen interaction RPM1-interacting protein 4 0.1 24/0 Root
Putative calcium-binding protein CML19 0.1 24/0 Root
Probable LRR receptor-like serine/threonine-protein kinase At3g47570 0.1 24/0 Root
LRR receptor-like serine/threonine-protein kinase FLS2 0.1 24/0 Root
Response to ethylene Ethylene receptor 2 0.03 24/0 Root
REVERSION-TO-ETHYLENE SENSITIVITY1 0.03 24/0 Root
Protein EIN4 0.03 24/0 Root
Defense response to other organism Protein SRC2 0.05 72/0; 24/0 Stem
MLO-like protein 1 0.05 72/0; 24/0 Stem
Defensin J1-2-like 0.05 72/0; 24/0 Stem
Bon1-association protein 2 0.05 72/0; 24/0 Stem
Pathogenesis-related protein R major form 0.1 72/24 Leaves
Pathogenesis-related protein P2 0.1 72/24 Leaves
Putative late blight resistance protein homolog R1A-3 0.1 72/24 Leaves
Disease resistance protein RGA2-like 0.1 72/24 Leaves
Defense response to fungus Acidic endochitinase pcht28 0.05 72/24 Leaves
Early nodulin-75 0.05 72/24 Leaves
Wound-induced protein WIN2 pathogenesis-related protein P2 0.05 72/24 Leaves
Ethylene-responsive transcription factor 5-like 0.05 72/0; 72/24 Stem
Cell wall biogenesis Cellulose synthase-like protein D3 0.05 24/0 Stem
Xyloglucan endotransglucosylase/hydrolase protein 15-like 0.05 24/0 Stem
Response to reactive oxygen species L-ascorbate peroxidase 2 cytosolic 0.05 24/0 Stem
L-ascorbate peroxidase 3 0.05 24/0 Stem
L-ascorbate peroxidase 1 cytosolic-like 0.05 24/0 Stem
Phenylpropanoid biosynthesis Caffeoyl-CoA O-methyltransferase 0.1 24/0; 72/0 Stem
Plant-pathogen interactionMAPK signaling pathway Respiratory burst oxidase homolog protein B 0.05 72/0 Stem
Mitogen-activated protein kinase 3 0.05 72/0 Stem
WRKY transcription factor 22 0.05 72/0 Stem
Mitogen-activated protein kinase homolog MMK2 0.05 72/0 Stem

Bubble plots showing the significantly enriched KEGG functions for various time (hours, h post-inoculation) contrasts in CM-334 (resistant variety): (a) 24h/0h, (b) 72h/0h, and (c) 72h/24h across all tissue samples.

Figure 3 shows gene enrichment analysis results across all tissue samples for various time contrasts. Plot (a) shows 24h/0h contrast; Plot (b) shows 72h/0h contrast; and Plot (c) focuses on 72h/24h Network analyses of differentially expressed genes

Genes related to immune response against fungal infection were differentially expressed across all time points. Network analyses of DEGs for the R vs. S contrast demonstrated that genes involved in defense response to fungal infection across all timepoints were also involved with carbohydrate metabolism and ADP binding (Figure 4). The DEGs with roles related to defense responses to other organisms were also involved in ADP binding and were differentially expressed in 0hpi and 24hpi with P. capsici. Genes with roles related to transmembrane transporter and oxidoreductase activities were differentially expressed in 24- and 72, and 0-, 24-, and 72hpi, respectively. Genes with functions associated with methylation, lipid transport, mRNA transcription, response to light, water transport, and endopeptidase inhibitor activity were differentially expressed only at the 0hpi. Genes involved in positive regulation of transcription by RNA Pol II were differentially expressed exclusively at 24hpi with P. capsici.

Network of significantly enriched GO functions in resistant (CM-334) vs susceptible (early jalapeño) samples from different tissues. Results from genes with higher expression is shown. The size of a function node is proportional to the number of differentially expressed genes (DEGs) annotated with a given function in all of the comparisons. The DEGs (small nodes) are shown connected to their corresponding functional categories. Function nodes are connected through shared DEGs. The DEG nodes are colored according to their occurrence in different contrasts. A function node is colored by the presence of DEGs annotated with that function in different contrasts, independently of the enrichment level of the function in each contrast.

Network diagram illustrating gene interactions at different times, represented by nodes and circles. Nodes are color-coded: red for Time 0, green for Time 24, and blue for Time 72. Various functions like defense response, protein folding, and mRNA transcription are labeled. The size of the circle represents the number of genes for each biological function.

Genes with functions associated with defense response to other organisms showed higher expression across all tissue samples, and shared roles in binding of ADP and monooxygenase and oxidoreductase activities (Supplementary Figure 2). Similarly, genes with roles in protein folding and alpha-tubulin binding were observed to have higher expression in the roots, stems, and leaves, and shared functions with ADP binding. Genes with roles related to protein dimerization, biosynthetic process of aromatic compounds, methyltransferase activity, and response to light showed higher expression exclusively in the leaves, whereas those involved in secondary shoot formation regulation were unique to the stem. Endopeptidase inhibitor activity, transmembrane transport, and transmembrane transporter activity-related genes were observed to have higher expression in the roots whereas genes with functions related to iron binding and DNA-binding transcription activity showed higher expression on both the leaf and stem tissues.

 Discussion

Phytophthora blight caused by P. capsici remains one of the most destructive diseases affecting chili pepper production in many areas of the world where the crop is cultivated. Understanding the genomic and transcriptomic basis of plant infection by P. capsici will facilitate robust and more informed breeding and selection decisions for improving resistance in current Capsicum germplasm (Barchenger et al., 2018a; Lozada et al., 2022). In the present study, gene expression analyses identified DEGs from different tissues and demonstrated the diversity of candidate genes involved in resistance to P. capsici in chili pepper. Comparative transcriptomics supported the presence of tissue-wide genetic resistance across different organs and tissue-specific mechanisms of the pepper plant under infection by the pathogen.

Root and stem showed similar transcriptome profiles for each biological replicate, where the resistant (R) and susceptible (S) samples clustered into their respective groups in the PCA biplot. In contrast, the leaf tissue demonstrated less consistent, overlapping clusters between the R and S samples. This could demonstrate relatively delayed and/or weaker immune responses in the leaves, compared to the roots and stems. This differential responses to P. capsici infection across the various tissues is consistent with the number of genes with higher expression, where the roots showed the highest number followed by the stem, and the leaf tissues. This could also be a possible consequence of the inoculation method used, where the pathogen was introduced via the soil; therefore, the initial contact occured in the roots.

The decreased number of DEGs in the 72hpi could be due to the temporal nature of the dynamics of the interaction between P. capsici and chili pepper. The initial response in the 0-24hpi, for example, could have resulted in higher transcriptional responses in the pepper plant, relative to the 72hpi, where responses have decreased, resulting in a lesser number of DEGs. In contrast to the roots and leaves where the number of genes with differential expression decreased as time of infection progressed, the stems showed an increased number of unique DEGs at 72hpi. The stems among the tissues evaluated for RNA-seq seem to have distinct mechanisms of resistance to P. capsici infection compared to the roots and leaves, as more genes with diverse functions showed higher differential expression in the stem (Table 3). The genes involved in constitutive or passive defenses (e.g., cell wall formation) were in general differentially expressed in the earlier phase of infection. In the later stage, those genes related to active immunity such as those with roles related to mitogen activated protein kinase signaling pathway, defense response, and systemic acquired resistance were differentially expressed. This indicated a shift from a passive to more active form of immunity as P. capsici infection progressed, albeit the number of differentially expressed genes generally decreased.

The majority of the DEGs (>80%) were exclusive to each tissue, indicating varying responses for each different part of the chili pepper plant upon infection by P. capsici. It was previously noted that infection of stem, leaf (foliar), and root caused by P. capsici demonstrate distinct, multiple disease syndromes, and each tissue will have different modes of action to reduce the effects of the pathogen (Sy et al., 2005). Remarkably, some genes showed higher expression across all tissue samples indicating the presence of shared disease responses among the different tissue types. These include the putative late blight resistance protein homolog R1B-16, alpha-crystallin domain-containing protein 22.3-like, putative disease resistance protein RGA3, and CSC1-like protein ERD4, with functions related to defense response to fungus and other organisms, protein kinase activity, and ADP and ATP binding which showed higher expression across all tissues in the R vs. S contrast.

Gene network analyses supported that genes with functions associated with defense responses to other organisms showed higher expression across all tissue samples in the R vs. S comparison. In other plant-pathogen pathosystems, such genes have been identified in the context of both biotic and abiotic stresses. In soybean (Glycine max), genes associated with defense, replication of DNA, and iron homeostasis demonstrated differential expression in root, stem, and leaf tissues and were regarded as “hallmarks of resistance” to brown stem rot caused by the fungus Phialophora gregata (McCabe et al., 2018). Genes related to phenylpropanoid biosynthesis, phenylalanine ammonia lyase (PAL), and peroxidase enzyme activity observed to have higher expression in the stems, roots, and leaves of black pepper varieties challenged with P. capsici (Hao et al., 2016). Recently, signal transduction pathways of plant hormones including abscisic acid (ABA) biosynthesis and signaling related genes (9-cis-epoxycarotenoid dioxygenase (NCED), phosphatase 2C (PP2C), SNF1-related protein kinase 2 (SnRK2), and ABA-responsive element binding factors (ABF)) were observed to have higher expression in both leaf and root tissue samples in response to drought stress in river tamarind (Leucaena leucocephala) (Zhi et al., 2024).

Chromosome 5 is a known major genomic hotspot for P. capsici-resistant QTLs and metaQTLs in chili pepper (Chunthawodtiporn et al., 2019; Kaur et al., 2024; Liu et al., 2014; Lozada et al., 2021a, 2021b; Mallard et al., 2013). An extended disease resistance QTL region, Ext-Pc5.1, pinpointed at 8.35 Mb–38.13 Mb on chromosome 5, has been previously described (Du et al., 2021). In the current study, transcriptomic analyses across all tissue samples revealed the higher expression of a gene located within this chromosome 5 hotspot between ~20.99 and 21 Mb with a homoserine kinase (HSK) activity function. The C. annuum downy mildew resistance 1 gene (CaDMR1), an HSK, and an orthologue of DMR1 that catalyzes the conversion of homoserine to homoserine-4-phosphate, a step in the aspartic acid-derived threonine biosynthetic pathway in Arabidopsis (van Damme et al., 2009) has been previously identified as a candidate gene for resistance to P. capsici at ~24.36 Mb in chromosome 5 (Rehrig et al., 2014). Due to the colocalization of the HSK gene identified in the present study with the CaDMR1 in chromosome 5, it could be possible that the same resistance gene was identified. The mechanism of homoserine-induced resistance remains unclear and warrants further gene validation in future studies. Genes coding for serine/threonine protein kinases and leucine rich repeat (LRR) receptors showed higher expression in the current study consistent with previous reports (Siddique et al., 2019; Du et al., 2021; Kang et al., 2022a) confirming their involvement in host defense against infection by P. capsici. The observed higher expression of LRR-receptor-like serine/threonine-protein kinase FLS2 in the roots of chili pepper in this study relates to the previous report on the role of flagellin perception FLS2 gene in Arabidopsis (Chinchilla et al., 2006). Overall, results in the present study supported the role of HSK and amino acid metabolism in the regulation of plant immunity (Heinemann and Hildebrandt 2021; Moormann et al., 2022; Zeier, 2013) for resistance to pathogen attack by P. capsici in chili pepper.

Genes with functions related to systemic acquired resistance (SAR) and defense responses to fungi and other organisms were observed to have high gene expression in the resistant cultivar (CM-334) in multiple tissues across various time vs. time comparisons. As a form of inducible resistance, SAR renders a long-lasting systemic immunity against diverse pathogens and is triggered following a previously localized exposure to a pathogen (Choi and Hwang, 2011). The putative lipid-transfer protein DEFECTIVE INDUCED RESISTANCE 1 (DIR1), a member of the lipid transfer proteins (LTPs) and involved in SAR, demonstrated higher expression in the 72h/0h, 24h/0h, and 72h/24h contrasts in all tissue samples for CM-334 in the current study (Table 3). Consistent with these observations, the higher expression of DIR1 has been previously identified in the leaf tissues of two pepper landraces in response to plant infection to P. capsici (Rabuma et al., 2022). Recently, the overexpression of a DIR1 gene in Arabidopsis was observed to enhance response to biotic stress by regulating genes related to defense against diseases and the formation of flavonoids (Duan et al., 2024). LTPs such as DIR1 may function as either a co-signal or translocator for the release of the mobile signals during SAR in plants (Jung et al., 2003). A pepper-specific SAR gene, CaSAR8.2, has been previously identified as a genetic marker for plant infection, and an elicitor of abiotic and environmental stresses in C. annuum (Lee and Hwang, 2003). Kaur et al. (2024) previously identified an SAR8.2-protein coding gene ~470 kb downstream as a candidate gene for SNP S5_14665044 detected using multi-locus association mapping approaches in a diverse population of C. annuum and C. annuum x C. frutescens hybrids. In another study, SAR8.2 was identified in the 27.3–29.2 Mb region of chromosome 5 (within the Ext-Pc5.1 region identified by Du et al. (2021)) as CASAR82A, suggesting its important role as a potential candidate gene for resistance to P. capsici (Siddique et al., 2019). No other SAR gene has been identified in chromosome 5 in the current study; however, various genes within the Ext-Pc5.1 region including a receptor kinase-like protein, Xa21, located < 0.90 Mb downstream of the SNP identified by Kaur et al. (2024) and a kunitz trypsin inhibitor 5 that induces programmed cell death in Arabidopsis (Li et al., 2008) were observed to have high gene expression across the R vs. S contrasts for all tissues infected by P. capsici. These specific genomic regions with known resistance genes can be targeted further for molecular marker-assisted breeding and selection to improve disease resistance in current pepper germplasm.

Genes with roles related to the formation of cell wall, phenylpropanoid biosynthesis, ethylene activated signaling pathway, and mitogen activated protein kinase (MAPK) showed high gene expression exclusively in the stems in the R vs. S contrast across various times of infection. As the primary barrier of defense, the cell wall provides the pepper plants with a preliminary defense and signal perception against pathogen infection through various detection and defense components (Lei et al., 2023). In Arabidopsis mutations in cellulose synthases altered the disease response to the fungus Plectosphaerella cucumerina (Hernandez-Blanco et al., 2007). Cellulose synthase-like D3 (CSLD3) and the xyloglucan endotransglucosylase/hydrolase protein 15-like showed in higher expression the stem at 24h/0h in the current study. While the function of CSLD3 has not been well characterized in Capsicum spp., the overexpression of this gene in cotton (Gossypium hirsutum) has been observed to restore cell elongation and cell wall integrity by enhancing primary cellulose production (Hu et al., 2019). The pathogen damage to the cell wall triggers phenylpropanoid pathway possibly leading to lignin biosynthesis and thereby producing various plant defense metabolites (Yadav et al., 2020). The involvement of phenylpropanoid biosynthesis pathway has been previously demonstrated in the whole root transcriptome of pepper in response to P. capsici (Li et al., 2020) and in the leaf and sheath tissues of rice (Oryza sativa) upon insect herbivory (Li et al., 2021). Another class of enzymes, the MAPKs, are involved in signaling plant defense against pathogen attack through receptors/sensors that transduce extracellular stimuli into intracellular responses (Meng and Zhang, 2013). In the resistant host–pathogenic microbe pathosystem, the resistant hosts directly induce pattern-recognition receptor (PRR)-triggered immunity (PTI), such as calcium ion (Ca2+) influx, reactive oxygen species (ROS) production, and MAPK activation (Ding et al., 2022). Overall, our results also demonstrate a layer of tissue-specific mechanisms to combat pathogen attack in the Capsicum-P. capsici pathosystem.

Current results from RNA-seq highlight various genes showing constant higher expression levels across various tissues under oomycete infection in C. annuum L. (Figure 5). Among these genes, DIR1, homoserine kinase, disease resistance proteins R1B-16 and RGA3, and CSC1-like proteins consistently showed high expression across all tissues, indicating effector-triggered tissue-wide mechanisms of disease resistance against P. capsici. Altogether, resistance to P. capsici infection in chili pepper results from the active expression of genes involved in various biological pathways related to the biosynthesis of amino acids, responses to biotic stresses, plant defense signaling, and adaptation. Overall, our results suggest that P. capsici resistance could result from the active expression of genes involved in various plant innate defenses such as PTI and effector-triggered immunity. These genes could be leveraged for developing broad-spectrum resistance against P. capsici in current Capsicum germplasm.

Model for Phytophthora capsici resistance in chile pepper derived from comparative transcriptome analyses of roots, stems, and leaf tissue samples. RNA-seq revealed a tissue-wide resistance mechanism conferred by genes that consistently showed higher expression across all tissue samples. Figure created using BioRender (https://www.biorender.com/).

Diagram illustrating the gene expression mechanism in a chile pepper plant affected by *P. capsici* infection. The plant is divided into sections: leaves, stems, and roots, each with associated proteins labeled. Biological processes include systemic resistance and defense response. Central starburst indicates “Tissue-wide gene expression mechanism.” Created in BioRender. Conclusions

RNA-seq of resistant and susceptible chili pepper varieties under infection by P. capsici and three different tissue types (leaves, stems, and roots) revealed a multitude of genes that showed higher and lower expression. The genes with higher expression have functions related to homoserine kinase activity, biogenesis of the cell wall, biosynthesis of phenylpropanoid, ethylene-activated signaling pathway, mitogen-activated protein kinase, systemic acquired resistance, plant pathogen interaction, and defense response, among others. Transcriptome analyses supported the role of chromosome 5 as a major genomic region harboring resistance genes, particularly homoserine kinase, which showed higher expression across all tissues. The genes with higher expression across all tissue samples indicate the presence of potential tissue-wide gene expression mechanism that exists in chili pepper. One limitation of the current study is the absence of qRT-PCR validation for genes showing differences in expression levels. The robustness of these results, however, could be supported through using stringent statistical criteria (Benjamini-Hochberg False Discovery Rate) in identifying genes with higher or lower expression. Gene expression profiling will serve as a basis for future studies linking gene expression to resistance to P. capsici. The existence of tissue-wide genetic mechanisms as revealed by transcriptomic profiling can be used for the breeding and selection of P. capsici-resistant chili pepper cultivars.

 Data availability statement

The data presented in this study are publicly available. The raw FASTQ files of the transcriptomes can be found here: https://www.ncbi.nlm.nih.gov, accession PRJNA1256908.

 Author contributions

DL: Conceptualization, Formal analysis, Funding acquisition, Project administration, Supervision, Writing – original draft, Writing – review & editing. PC-A: Conceptualization, Formal analysis, Funding acquisition, Methodology, Project administration, Writing – review & editing. JM: Investigation, Writing – review & editing. JA: Conceptualization, Formal analysis, Funding acquisition, Writing – review & editing. VV-M: Methodology, Writing – review & editing. NK: Investigation, Methodology, Writing – review & editing. ED: Funding acquisition, Resources, Writing – review & editing. SS: Funding acquisition, Resources, Writing – review & editing.

 Acknowledgments

The Authors are grateful to Dr. Joe Song (New Mexico State University, NMSU) for reviewing the manuscript. The technical assistance of Aurora Labastida (OMICS Analysis, Mexico), Dr. Scot Dowd (Molecular Research LP (MR DNA), TX, USA), and the staff and students at the NMSU Chile Pepper Breeding and Genetics Program are also appreciated.

 Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

 Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2026.1745952/full#supplementary-material

Multidimensional scaling (MDS) or Principal coordinates analysis (PCoA) for (A) roots, (B) stems, and (C) leaves under Phytophthora capsici infection demonstrating clustering based on resistance and susceptibility of varieties across various times of infection for the root and stem tissue samples. An apparent overlap of clusters was observed for the leaves.

Network of Significantly Enriched GO Functions in Resistant (CM-334) vs Susceptible (Early Jalapeno) samples from different tissues. Results genes with significant high and lower expression are shown. The size of a function node is proportional to the number of differentially expressed (DE) genes annotated with a given function in all the comparisons. The DE genes (small nodes) are shown connected to their corresponding functional categories. Function nodes are connected through shared DE genes. DE gene nodes are colored according to their occurrence in different contrasts. A function node is colored by the presence of DE genes annotated with that function in different contrasts, independently of the enrichment level of the function in each contrast.

List of software tools and packages used for RNA-sequencing analyses.

Total number of reads, percentage of pseudoaligned reads, number of unique sequences, and percentage of unique sequences identified using Kallisto.

Different comparisons of transcriptome profiles performed for C. annuum under P. capsici infection.

Top 100 genes with higher expression for C. annuum under P. capsici infection (Resistant vs Susceptible) across all tissues.

Top 100 genes with higher expression for C. annuum under P. capsici infection (tissue vs. tissue) at 0 hpi.

Top 100 genes with higher expression for C. annuum under P. capsici infection (tissue vs. tissue) at 24 hpi.

Top 100 genes with higher expression for C. annuum under P. capsici infection (tissue vs. tissue) at 72 hpi.

Top 100 genes with higher expression for C. annuum under P. capsici infection (tissue vs tissue) in CM-334 (Resistant variety).

Top 100 genes with higher expression for C. annuum under P. capsici infection in the roots (time vs. time).

Top 100 genes with higher expression for C. annuum under P. capsici infection (time vs time) in the stems.

 References Barchenger D. W. Lamour K. H. Bosland P. W. (2018a). Challenges and strategies for breeding resistance in Capsicum annuum to the multifarious pathogen, Phytophthora capsici. Front. Plant Sci. 9, 628. doi: 10.3389/fpls.2018.00628, PMID: 29868083 Barchenger D. W. Sheu Z.-M. Kumar S. Lin S.-W. Burlakoti R. R. Bosland P. W. (2018b). Race characterization of Phytophthora root rot on Capsicum in Taiwan as a basis for anticipatory resistance breeding. Phytopathology 108, 964–971. doi: 10.1094/PHYTO-08-17-0289-R, PMID: 29484915 Bolger A. M. Lohse M. Usadel B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170, PMID: 24695404 Bray N. L. Pimentel H. Melsted P. Pachter L. (2016). Near-optimal probabilistic RNA-seq quantification. Nat. Biotechnol. 34, 525–527. doi: 10.1038/nbt.3519, PMID: 27043002 Chinchilla D. Bauer Z. Regenass M. Boller T. Felix G. (2006). The Arabidopsis receptor kinase FLS2 binds flg22 and determines the specificity of flagellin perception. Plant Cell 18, 465–476. doi: 10.1105/tpc.105.036574, PMID: 16377758 Choi H. W. Hwang B. K. (2011). Systemic acquired resistance of pepper to microbial pathogens. J. Phytopathol. 159, 393–400. doi: 10.1111/j.1439-0434.2010.01781.x, PMID: 41711423 Chunthawodtiporn J. Hill T. Stoffel K. Van Deynze A. (2019). Genetic Analysis of Resistance to Multiple Isolates of Phytophthora capsici and Linkage to Horticultural Traits in Bell Pepper. HortScience 54, 1143–1148. doi: 10.21273/HORTSCI13359-18 Ding L.-N. Li Y.-T. Wu Y.-Z. Li T. Geng R. Cao J. . (2022). Plant disease resistance-related signaling pathways: recent progress and future prospects. Int. J. Mol. Sci. 23, 16200. doi: 10.3390/ijms232416200, PMID: 36555841 Du J.-S. Hang L.-F. Hao Q. Yang H.-T. Ali S. Badawy R. S. E. . (2021). The dissection of R genes and locus Pc5. 1 in Phytophthora capsici infection provides a novel view of disease resistance in peppers. BMC Genomics 22, 1–16. doi: 10.1186/s12864-021-07705-z, PMID: 34016054 Duan M. Bao L. Eman M. Han D. Zhang Y. Zheng B. . (2024). The ectopic expression of the mpDIR1 (t) gene enhances the response of plants from arabidopsis thaliana to biotic stress by regulating the defense genes and antioxidant flavonoids. Plants 13, 2692. doi: 10.3390/plants13192692, PMID: 39409562 Erwin D. C. Ribeiro O. K. (1996). Phytophthora diseases worldwide (St. Paul, MN: The American Phytopathological Society). Fox J. Weisberg S. (2011). An R companion to applied regression. (Thousand Oaks, CA, USA: Sage Publishing). Friendly M. Monette G. Fox J. (2013). Elliptical insights: understanding statistical methods through elliptical geometry. Stat. Sci. 28, 1–39. doi: 10.1214/12-STS402, PMID: 32930232 Gower J. C. (1966). Some distance properties of latent root and vector methods used in multivariate analysis. Biometrika 53, 325–338. doi: 10.1093/biomet/53.3-4.325, PMID: 36970824 Griffith M. Walker J. R. Spies N. C. Ainscough B. J. Griffith O. L. (2015). Informatics for RNA sequencing: a web resource for analysis on the cloud. PloS Comput. Biol. 11, e1004393. doi: 10.1371/journal.pcbi.1004393, PMID: 26248053 Hao C. Xia Z. Fan R. Tan L. Hu L. Wu B. . (2016). De novo transcriptome sequencing of black pepper (Piper nigrum L.) and an analysis of genes involved in phenylpropanoid metabolism in response to Phytophthora capsici. BMC Genomics 17, 1–14. doi: 10.1186/s12864-016-3155-7, PMID: 27769171 Heinemann B. Hildebrandt T. M. (2021). The role of amino acid metabolism in signaling and metabolic adaptation to stress-induced energy deficiency in plants. J. Exp. Bot. 72, 4634–4645. doi: 10.1093/jxb/erab182, PMID: 33993299 Hernandez-Blanco C. Feng D. X. Hu J. Sanchez-Vallet A. Deslandes L. Llorente F. . (2007). Impairment of cellulose synthases required for Arabidopsis secondary cell wall formation enhances disease resistance. Plant Cell 19, 890–903. doi: 10.1105/tpc.106.048058, PMID: 17351116 Hieter P. Boguski M. (1997). Functional genomics: it’s all how you read it. Sci. (80-.). 278, 601–602. doi: 10.1126/science.278.5338.601, PMID: 9381168 Hu H. Zhang R. Tang Y. Peng C. Wu L. Feng S. . (2019). Cotton CSLD3 restores cell elongation and cell wall integrity mainly by enhancing primary cellulose production in the Arabidopsis cesa6 mutant. Plant Mol. Biol. 101, 389–401. doi: 10.1007/s11103-019-00910-1, PMID: 31432304 Hulse-Kemp A. M. Maheshwari S. Stoffel K. Hill T. A. Jaffe D. Williams S. R. . (2018). Reference quality assembly of the 3.5-Gb genome of Capsicum annuum from a single linked-read library. Hortic. Res. 5, 1–13. doi: 10.1038/s41438-017-0011-0, PMID: 29423234 Jiang L. Sanogo S. Bosland P. W. (2015). Using recombinant inbred lines to monitor changes in the race structure of Phytophthora capsici in Chile pepper in New Mexico. Plant Heal. Prog. 16, 235–240. doi: 10.1094/PHP-RS-15-0034, PMID: 40211709 Jung H. W. Kim W. Hwang B. K. (2003). Three pathogen-inducible genes encoding lipid transfer protein from pepper are differentially activated by pathogens, abiotic, and environmental stresses. Plant Cell Environ. 26, 915–928. doi: 10.1046/j.1365-3040.2003.01024.x, PMID: 12803619 Kamoun S. Furzer O. Jones J. D. G. Judelson H. S. Ali G. S. Dalio R. J. D. . (2015). The Top 10 oomycete pathogens in molecular plant pathology. Mol. Plant Pathol. 16, 413–434. doi: 10.1111/mpp.12190, PMID: 25178392 Kanehisa M. Goto S. (2000). KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28, 27–30. doi: 10.1093/nar/28.1.27, PMID: 10592173 Kang W.-H. Lee J. Koo N. Kwon J.-S. Park B. Kim Y.-M. . (2022a). Universal gene co-expression network reveals receptor-like protein genes involved in broad-spectrum resistance in pepper (Capsicum annuum L.). Hortic. Res. 9, uhab003. doi: 10.1093/hr/uhab003, PMID: 35043174 Kang W.-H. Nam J.-Y. Kwon J.-S. Yeom S.-I. (2022b). Differential responses of Capsicum accessions to foliar blight and root rot by Phytophthora capsici. Hortic. Sci. Technol. 40, 85–93. doi: 10.7235/HORT.20220009, PMID: 28809037 Kaur N. Lozada D. N. Bhatta M. Barchenger D. W. Khokhar E. S. Nourbakhsh S. S. . (2024). Insights into the genetic architecture of Phytophthora capsici root rot resistance in Chile pepper (Capsicum spp.) from multi-locus genome-wide association study. BMC Plant Biol. 24, 416. doi: 10.1186/s12870-024-05097-2, PMID: 38760676 Kolde R. Kolde M. R. (2015). Package ‘pheatmap.’ R packag, Vol. 1. 790. Lamour K. H. Stam R. Jupe J. Huitema E. (2012). The oomycete broad-host-range pathogen Phytophthora capsici. Mol. Plant Pathol. 13, 329–337. doi: 10.1111/j.1364-3703.2011.00754.x, PMID: 22013895 Lee S. Hwang B. (2003). Identification of the pepper SAR8. 2 gene as a molecular marker for pathogen infection, abiotic elicitors and environmental stresses in Capsicum annuum. Planta 216, 387–396. doi: 10.1007/s00425-002-0875-5, PMID: 12520329 Lei G. Zhou K.-H. Chen X.-J. Huang Y.-Q. Yuan X.-J. Li G.-G. . (2023). Transcriptome and metabolome analyses revealed the response mechanism of pepper roots to Phytophthora capsici infection. BMC Genomics 24, 626. doi: 10.1186/s12864-023-09713-7, PMID: 37864214 Leonian L. H. (1922). Stem and fruit blight of peppers caused by Phytophthora capsici sp. nov. Phytopathology, 12. Li J. Brader G. Palva E. T. (2008). Kunitz trypsin inhibitor: an antagonist of cell death triggered by phytopathogens and fumonisin b1 in Arabidopsis. Mol. Plant 1, 482–495. doi: 10.1093/mp/ssn013, PMID: 19825555 Li Y. Cheah B. H. Fang Y.-F. Kuang Y.-H. Lin S.-C. Liao C.-T. . (2021). Transcriptomics identifies key defense mechanisms in rice resistant to both leaf-feeding and phloem feeding herbivores. BMC Plant Biol. 21, 1–18. doi: 10.1186/s12870-021-03068-5, PMID: 34193042 Li Y. Yu T. Wu T. Wang R. Wang H. Du H. . (2020). The dynamic transcriptome of pepper (Capsicum annuum) whole roots reveals an important role for the phenylpropanoid biosynthesis pathway in root resistance to Phytophthora capsici. Gene 728, 144288. doi: 10.1016/j.gene.2019.144288, PMID: 31846710 Liu W. Y. Kang J. H. Jeong H. S. Choi H. J. Yang H. B. Kim K. T. (2014). Combined use of bulked segregant analysis and microarrays reveals SNP markers pinpointing a major QTL for resistance to Phytophthora capsici in pepper. TAG. Theoretical and Applied Genetics 127, 2503–2513. doi: 10.1007/s00122-014-2394-8, PMID: 25208646 Love M. I. Huber W. Anders S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 1–21. doi: 10.1186/s13059-014-0550-8, PMID: 25516281 Lozada D. N. Bosland P. Barchenger D. W. Haghshenas-Jaryani M. Sanogo S. Walker S. (2022). Chile pepper (Capsicum) breeding and improvement in the “Multi-omics” Era. Front. Plant Sci. 13, 1363. doi: 10.3389/fpls.2022.879182, PMID: 35592583 Lozada D. N. Nunez G. Lujan P. Dura S. Coon D. Barchenger D. W. . (2021a). Genomic regions and candidate genes linked with Phytophthora capsici root rot resistance in Chile pepper (Capsicum annuum L.). BMC Plant Biol. 21, 601. doi: 10.1186/s12870-021-03387-7, PMID: 34922461 Lozada D. N. Whelpley M. Acuña-Galindo A. (2021b). Genetic architecture of Chile pepper (Capsicum spp.) QTLome revealed using meta-QTL analysis. Horticulturae 7, 227. doi: 10.3390/horticulturae7080227, PMID: 41700957 Maillot G. Szadkowski E. Massire A. Brunaud V. Rigaill G. Caromel B. . (2022). Strive or thrive: Trends in Phytophthora capsici gene expression in partially resistant pepper. Front. Plant Sci. 13, 980587. doi: 10.3389/fpls.2022.980587, PMID: 36479518 Mallard S. Cantet M. Massire A. Bachellez A. Ewert S. Lefebvre V. (2013). A key QTL cluster is conserved among accessions and exhibits broad-spectrum resistance to Phytophthora capsici: a valuable locus for pepper breeding. Mol. Breed. 32, 349–364. doi: 10.1007/s11032-013-9875-3, PMID: 41721156 Martin M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 17, 10–12. doi: 10.14806/ej.17.1.200 McCabe C. E. Cianzio S. R. O’Rourke J. A. Graham M. A. (2018). Leveraging RNA-Seq to characterize resistance to Brown stem rot and the Rbs3 locus in soybean. Mol. Plant-Microbe Interact. 31, 1083–1094. doi: 10.1094/MPMI-01-18-0009-R, PMID: 30004290 Meng X. Zhang S. (2013). MAPK cascades in plant disease resistance signaling. Annu. Rev. Phytopathol. 51, 245–266. doi: 10.1146/annurev-phyto-082712-102314, PMID: 23663002 Monroy-Barbosa A. Bosland P. W. (2008). Genetic analysis of Phytophthora root rot race-specific resistance in Chile pepper. J. Am. Soc Hortic. Sci. 133, 825–829. doi: 10.21273/JASHS.133.6.825 Moormann J. Heinemann B. Hildebrandt T. M. (2022). News about amino acid metabolism in plant–microbe interactions. Trends Biochem. Sci. 47, 839–850. doi: 10.1016/j.tibs.2022.07.001, PMID: 35927139 Morgan M. Carlson M. Tenenbaum D. Arora S. (2019). Package ‘AnnotationHub.’. Quesada-Ocampo L. M. Parada-Rojas C. H. Hansen Z. Vogel G. Smart C. Hausbeck M. K. . (2023). Phytophthora capsici: Recent progress on fundamental biology and disease management 100 years after its description. Annu. Rev. Phytopathol. 61, 185–208. doi: 10.1146/annurev-phyto-021622-103801, PMID: 37257056 Rabuma T. Gupta O. P. Yadav M. Chhokar V. (2022). Integrative RNA-Seq analysis of Capsicum annuum L.-Phytophthora capsici L. pathosystem reveals molecular cross-talk and activation of host defence response. Physiol. Mol. Biol. Plants 28, 171–188. doi: 10.1007/s12298-021-01122-y, PMID: 35221578 Rehrig W. Z. Ashrafi H. Hill T. Prince J. Van Deynze A. (2014). CaDMR1 cosegregates with QTL Pc5. 1 for resistance to Phytophthora capsici in pepper (Capsicum annuum). Plant Genome 7, 1–12. doi: 10.3835/plantgenome2014.03.0011 Retes-Manjarrez J. E. Rubio-Aragón W. A. Márques-Zequera I. Cruz-Lachica I. García-Estrada R. S. Sy O. (2020). Novel Sources of Resistance to Phytophthora capsici on Pepper (Capsicum sp.) Landraces from Mexico. Plant Pathol. J. 36, 600. doi: 10.5423/PPJ.OA.07.2020.0131, PMID: 33312095 Reyes-Tena A. Castro-Rocha A. Rodríguez-Alvarado G. Vázquez-Marrufo G. Pedraza-Santos M. E. Lamour K. . (2019). Virulence phenotypes on chili pepper for Phytophthora capsici isolates from Michoacán, Mexico. HortScience 54, 1526–1531. doi: 10.21273/HORTSCI13964-19 Sanogo S. Lamour K. Kousik C. S. Lozada D. N. Parada-Rojas C. H. Quesada-Ocampo L. M. . (2023). Phytophthora capsici, 100 years later: research mile markers from 1922 to 2022. Phytopathology 113, 921–930. doi: 10.1094/PHYTO-08-22-0297-RVW, PMID: 36401843 Sharma H. Shukla M. K. Bosland P. W. Steiner R. (2017). Soil moisture sensor calibration, actual evapotranspiration, and crop coefficients for drip irrigated greenhouse Chile peppers. Agric. Water Manage. 179, 81–91. doi: 10.1016/j.agwat.2016.07.001, PMID: 41723091 Siddique M. I. Lee H.-Y. Ro N.-Y. Han K. Venkatesh J. Solomon A. M. . (2019). Identifying candidate genes for Phytophthora capsici resistance in pepper (Capsicum annuum) via genotyping-by-sequencing-based QTL mapping and genome-wide association study. Sci. Rep. 9, 9962. doi: 10.1038/s41598-019-46342-1, PMID: 31292472 Soneson C. Love M. I. Robinson M. D. (2016). Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences. F1000Research 4, 1521. doi: 10.12688/f1000research.7563.2, PMID: 26925227 Sy O. Bosland P. W. Steiner R. (2005). Inheritance of Phytophthora stem blight resistance as compared to Phytophthora root rot and Phytophthora foliar blight resistance in Capsicum annuum L. J. Am. Soc Hortic. Sci. 130, 75–78. doi: 10.21273/JASHS.130.1.75 Sy O. Steiner R. Bosland P. W. (2008). Recombinant inbred line differential identifies race-specific resistance to Phytophthora root rot in Capsicum annuum. Phytopathology 98, 867–870. doi: 10.1094/PHYTO-98-8-0867, PMID: 18943204 van Damme M. Zeilmaker T. Elberse J. Andel A. de Sain-van der Velden M. van den Ackerveken G. (2009). Downy mildew resistance in Arabidopsis by mutation of HOMOSERINE KINASE. Plant Cell 21, 2179–2189. doi: 10.1105/tpc.109.066811, PMID: 19622802 Venables W. N. Ripley B. D. (2013). Modern applied statistics with S-PLUS. (Berlin, Germany: Springer Science & Business Media). Wang P. Liu X. Guo J. Liu C. Fu N. Shen H. (2015). Identification and expression analysis of candidate genes associated with defense responses to Phytophthora capsici in pepper line “PI 201234. Int. J. Mol. Sci. 16, 11417–11438. doi: 10.3390/ijms160511417, PMID: 25993303 Werner T. (2010). Next generation sequencing in functional genomics. Brief. Bioinform. 11, 499–511. doi: 10.1093/bib/bbq018, PMID: 20501549 Yadav V. Wang Z. Wei C. Amo A. Ahmed B. Yang X. . (2020). Phenylpropanoid pathway engineering: An emerging approach towards plant defense. Pathogens 9, 312. doi: 10.3390/pathogens9040312, PMID: 32340374 Zeier J. (2013). New insights into the regulation of plant immunity by amino acid metabolic pathways. Plant Cell Environ. 36, 2085–2103. doi: 10.1111/pce.12122, PMID: 23611692 Zhi Q.-Q. Chen Y. Hu H. Huang W.-Q. Bao G.-G. Wan X.-R. (2024). Physiological and transcriptome analyses reveal tissue-specific responses of Leucaena plants to drought stress. Plant Physiol. Biochem. 214, 108926. doi: 10.1016/j.plaphy.2024.108926, PMID: 38996715

Edited by: Parimalan Rangan, Indian Council of Agricultural Research (ICAR), India

Reviewed by: Can Baysal, University of Florida, United States

Weier Cui, University of Chile, Chile