Next Article in Journal
Analysis of Various Machine Learning Algorithms for Using Drone Images in Livestock Farms
Previous Article in Journal
Bridging the Gap: Scaling Up the Sustainable Production of the Yellow Mealworm with Agricultural By-Products—Insights into Larval Growth and Body Composition
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Agriculture
Volume 14
Issue 4
10.3390/agriculture14040521
agriculture-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcript Analysis Reveals Positive Regulation of CA12g04950 on Carotenoids of Pigment Pepper Fruit under Nitrogen Reduction

by
Lingfeng Shen
1,2,
Chenfei Zhang
1,2,
Yahui Xia
1,2,
Shasha Yang
1,2,
Tian Chang
1,2,
Saleem Ullah
1,2 and
Xuehua Ji
1,2,*
1
Department of Horticulture, Agricultural College, Shihezi University, Shihezi 832000, China
2
Xinjiang Production and Construction Corps Key Laboratory of Special Fruits and Vegetables Cultivation Physiology and Germplasm Resources Utilization, College of Agriculture, Shihezi University, Shihezi 832003, China
*
Author to whom correspondence should be addressed.
Agriculture 2024, 14(4), 521; https://doi.org/10.3390/agriculture14040521
Submission received: 22 February 2024 / Revised: 20 March 2024 / Accepted: 21 March 2024 / Published: 25 March 2024
(This article belongs to the Section Crop Production)
Browse Figures
Versions Notes

Abstract

:
This study investigates the relationship between nitrogen fertilization and pepper fruit color by employing five different nitrogen treatments (N1: 750 kg/hm2, N2: 562.5 kg/hm2, N3: 375 kg/hm2, N4: 187.5, and N0: 0 kg/hm2). Fruits were harvested at 30 (S1: green ripening stage), 45 (S2: color transition stage), and 60 days (S3: red ripening stage) after flowering. Subsequently, pigment content, carotenoid component content, carotenoid enzyme activity, and transcriptome sequence were analyzed, and CA12g04950 function was validated through virus-induced gene silencing (VIGS). The results indicate that a reduction in nitrogen application led to an earlier onset of fruit color breakdown, and increased the contents of total carotenoid, capsanthin, phytoene and PSY (phytoene synthase) activity, LCYB (lycopene β-cyclase) activity and CCS (capsanthin/capsorubin synthase) activity. The analysis of different expression genes indicated that the most differently expressed genes were enriched in the N1 vs. N4 comparison, with 18 genes involved in carotenoid metabolism and 16 genes involved in nitrogen metabolism. Most DE genes were enriched in the pathways of photosynthesis, porphyrin, carotenoid biosynthesis, seleno-compounds, and nitrogen metabolism. There were numerous differential transcription factor families, including ERF, bHLH, MYB, C2H2, and NAC. Pearson correlation analysis revealed a significant positive correlation between CA12g04950 expression and 11 carotenoid genes in the N4 treatment. Subsequent silencing of CA12g04950 using VIGS resulted in delayed color ripening while a significant decrease in total carotenoid content and the expression levels of carotenoid genes. In conclusion, nitrogen reduction led to an increase in carotenoid content in pigment pepper fruits. Furthermore, under nitrogen reduction, CA12g04950 positively influenced the redness of the fruits.

1. Introduction

Pepper, one of the most extensively cultivated and highly demanded vegetables worldwide, accounted for approximately 40 million tons of global production and consumption in 2019 [1]. Xinjiang Province, China, serves as the foremost production hub for pigment peppers in China, with the current cultivation area of pigment peppers exceeding 33,000 hectares. Pepper fruits are widely utilized as natural colorants in both the food and cosmetics industries due to their richness in flavonoids and carotenoids [2,3].
Carotenoids are a class of tetraterpene pigments that exhibit yellow, orange-red, or red hues and contain eight double bonds [4]. The distinctive capsanthin found in chili peppers is a significant contributor to the red coloration of chili fruits, constituting approximately 46% of the total carotenoid content. In orange pepper fruits, α-carotene, β-carotene, and zeaxanthin play crucial roles in imparting the orange-yellow color, with zeaxanthin constituting 16% of the total carotenoid content [2,5]. Benefiting from the potent antioxidant properties of carotenoids, they are now extensively utilized in pharmaceuticals, food, and cosmetics industries [2], with capsanthin exhibiting the highest market demand [1].
Previous investigations have focused on elucidating the biochemical properties and physiological significance of carotenoids and flavonoids, as well as singular constituents within diverse cultivars of chili fruits. Several studies have undertaken research on the phenotypic manifestations, quantitative compositions, and chromosomal allocations of carotenoids and flavonoids within botanical specimens [6,7,8,9,10]. Feng [11] used three wild-type red peppers and their corresponding mutants (yellow and orange) as research materials to clone the genes involved in carotenoid biosynthesis, including GGPS, PSY, LCYB, CRTZ, ZEP, and CCS. It was discovered that in two mutant types, the complete CCS gene could not be cloned. In another mutant type (orange), although the complete CCS gene could be cloned, its gene expression level was significantly lower than that of the wild type. Therefore, it is believed that the color variation in this mutant fruit is related to the regulatory mechanism of the CCS gene. Through the integration of transcriptomics and metabolomics approaches, it was discerned that the genes CCS, NCED2, AAO4, VDE1, and CYP97C1 serve as pivotal regulators in the synthesis of carotenoids [12].
Modulating the expression levels of key genes involved in the carotenoid biosynthesis pathway significantly influences the biosynthesis of carotenoids. Using gene editing techniques to induce specific overexpression of the PSY gene has been demonstrated as an effective strategy to enhance the content of carotenoids in plant tissues [13,14]. Concomitant overexpression of the DXS, PSY, and CRTI genes leads to a substantial augmentation in carotenoid content in rice endosperm [15]. After silencing the PSY, CCS, LCYB, and CRTZ genes in peppers, there are changes in the levels of precursor compounds essential for capsanthin synthesis. Consequently, this disrupts the normal coloration of chili fruit, causing a transition from red to yellow or orange [16,17]. Transcription factors, as a crucial class of signaling molecules, can directly or indirectly regulate the synthesis and accumulation of carotenoids in plants [18]. In citrus, CsMADS5 positively regulates carotenoid synthesis in fruits by directly binding to the promoters of PSY, PDS, and CYB1 and by initiating their transcription [19]. In tomato, a model plant of the Solanaceae family, numerous studies have identified the involvement of transcription factors in the regulation of carotenoid biosynthesis [20]. Overexpression of SlWRKY35 can increase the content of lutein in tomato fruit, consequently enhancing the accumulation of carotenoids in the fruit [21]. Overexpression of SlNAC1 reduces ethylene synthesis in tomato fruit, leading to a decrease in carotenoid content and resulting in fruit appearing yellow or orange [22]. Recent studies in peppers have revealed that a specific R-R-type MYB transcription factor can promote carotenoid biosynthesis, thereby facilitating the accumulation of capsanthin in peppers [23].
In addition to the influence of genes on the composition and content of carotenoids, research has found that cultivation measures such as fertilization and pruning can change the content and type of plant pigments [24,25]. Nitrogen (N) is one of the significant influencing factors affecting plant growth and development, as well as fruit quality and yield. Studies have demonstrated that within a specific range of nitrogen application, an increase in the amount of urea applied leads to an increase in capsanthin content in pepper fruits [26,27]. Research indicates that when 750 kg/ha of nitrogen is applied to Lycium chinense Miller, the content of flavonoids and carotenoids in the fruits reaches its highest level [28]. There are also studies indicating that the content of carotenoids in plants does not continuously increase with nitrogen fertilizer concentration. Research indicates that increasing nitrogen fertilization decreases the content of lycopene and total carotenoids in tomato fruit [29]. Currently, research on carotenoids in peppers primarily focuses on the variations in carotenoid components and total carotenoid content across different varieties. However, there is relatively scant investigation into the impacts of nitrogen fertilizer application on metabolites related to chili pigment and gene regulation. The relationship between cultivation factors like nitrogen fertilizer and fruit color, as well as carotenoid formation, remains ambiguous [30].
To elucidate the effects of nitrogen fertilizer on pepper fruit color, this study analyzed changes in carotenoid components’ content, enzyme activities, and gene expression of pigment pepper fruits under different urea concentrations. The transcription mechanism of nitrogen reduction on fruit carotenoids is explored through the analysis of DE gene (DEG) and metabolic pathways. Then, virus-induced gene silencing (VIGS) is employed to confirm the function of the carotenoid DE gene that is significantly expressed under the N4 treatment. This study aims to provide a theoretical guide for optimizing fertilization to enhance the color and commercial quality of pepper fruits.

2. Materials and Methods

2.1. Materials and Design

The experiment was conducted in Shihezi City, Xinjiang Province, China. The selected variety was the pigment pepper named ‘Honglong 23’, developed by Xinjiang Tianjiaohongan Company (Shihezi, China). Pepper seedlings, aged 60 days according to the calendar, were transplanted into plots measuring 12 m in length and 1 m in width. Each plot comprised two rows, with a spacing of 50 cm between the rows, and 30 cm between the plants. Drip irrigation under mulch was utilized for water and fertilizer supply. Nitrogen fertilizer, in the form of urea (N ≥ 46%), was applied 10 days after transplanting.
The urea dosages was set as follows: 750 kg/hm2 (N1, widely used in the fields and serving as the control for this study), 562.5 kg/hm2 (N2, a 25% reduction in urea compared to N1), 375 kg/hm2 (N3, a 50% reduction in urea compared to N1), 187.5 kg/hm2 (N4, a 75% reduction in urea compared to N1), and N0 without urea application. Urea was applied once every 10 days for a total of 8 times after 10 days of planting (Table 1). The soil’s organic matter, alkaline-dissolved nitrogen, available phosphorus, and available potassium contents were 11.8, 43.7, 18.9, and 155 mg/kg, respectively. Three replicate plots were established for each urea dosage. Fruits were sampled at 30 (S1: green ripening period), 45 (S2: color transition period), and 60 (S3: red ripening period) days after blooming. Each parameter was measured with three biological repeats.

2.2. Determination of Fruit Pigment Content and Enzyme Activity

Three pepper plants from each treatment were selected, and three fruits were sampled, totaling 9 fruits to serve as experimental samples. Approximately 0.1 g of mixed fruit was frozen and ground into powder using liquid nitrogen. The resulting sample underwent extraction with a 10 mL mixture (2:1, v/v) of 100% ethanol and acetone using an ultrasonic instrument for 30 min. Subsequently, the mixture was centrifuged at 12,000 rpm for 20 min, and the supernatant was collected. This extraction process was repeated 2 to 3 times, and the collected supernatant was adjusted to a constant volume of 25 mL. The solution was then stored in a brown glass bottle for the measurement of carotenoid contents using a spectrophotometer (UV-2600i) (Shimadzu, Shanghai, China) as described by Burgos et al. [31].
Enzyme activities of phytoene synthase (PSY), capsanthin/capsorubin synthase (CCS), and lycopene β-cyclase (LCYB) were analyzed using ELISA kits sourced from Shanghai Enzyme-linked Biotechnology Co., Ltd., Shanghai, China [32]. The samples were sent to Wuhan Maiwei Metabolism Company (Wuhan, China) for carotenoid friction determination using HPLC-MS, as described in [30,33].

2.3. RNA Extraction

Total RNA was isolated from N1, N3, and N4 fruits at stage S1 using Trizol (Invitrgen, Beijing, China) according to the manufacturer’s protocol. The extracted RNAs were cleaned using an RNeasy MinElute Cleanup Kit (Qiagen, Valencia, CA, USA).

2.4. RNA Sequencing

RNA was examined by Nanodrop spectrophotometry (Nanodrop Technologies, Wilmington, DE, USA) and a Nano 6000 Assay Kit from the Agilent Bioanalyzer 2100 system (Agilent Technologies, Santa Clara, CA, USA). A total of 1 μg RNA per sample was used to generate a cDNA library by NEBNext Ultra II RNA Library Prep Kit for Illumina (New England Biolabs Inc.; Ipswich, MA, USA). Messenger RNA (mRNA) was enriched by oligo beads, and synthesized cDNA. Used AMPure XP (Beckman Coulter, Brea, CA, USA) beads to select cDNA of approximately 400–500 bp in size, followed by PCR amplification, The PCR products were then purified again using AMPure XP beads to obtain 9 libraries, and the Sequel™ Sequencing Kit 2.0 was used to sequence 9 cDNA libraries. Then, 9 RNA-seq libraries were sequenced by the Illumina HiSeq 6000 (San Diego, CA, USA). cDNA library construction and RNA-seq were performed by Personal Bioinformatics Technology Co. (Shanghai, China). Low-quality reads were removed using Cutadapt 1.9.1 software, and the clean reads were mapped to the Pepper_1.55 reference genome (from the Sol Genomics Network) with the help of HISAT2 [34]. Based on the established reference genome, a mapping-based assembly approach was employed using the software StringTie 2.1.4 for assembling and splicing the mapped reads. The fragments per kilobase of transcript per million mapped reads (FPKM) for each gene was computed by considering the gene length and the count of reads mapped to it [35]. The raw transcriptome data underwent principal components analysis (PCA) and Pearson correlation based on gene expression profiles obtained from fragments per kilobase of exon per million fragments mapped (FPKM) among the nine samples using the package R 4.1.1. The identified transcripts underwent functional annotation against seven public databases, namely NR [36], KEGG [37], Swissprot [38], GO [39], and eggNOG [40], utilizing predefined parameters. GO enrichment analysis was carried out using GO-seq, and KEGG enrichment was performed with the hypergeometric test [41,42]. Differential expression analysis was performed using the DESeq R package [43]. Differentially expressed genes were identified with |log2FoldChange| > 1 and a significance p-value < 0.05. Transcription factor families were analyzed using the PlanTFDB database. The raw transcriptome data were deposited into the NCBI SRA (BioProject ID: PRJNA1087067).

2.5. Quantitative Real-Time PCR

To validate the RNA-seq data, 15 genes from carotenoid pathways detected by KEGG functional analysis were selected for quantitative real-time PCR (qRT-PCR). RNA extraction was performed using the RNAprep Pure Plant Kit (DP432; Tiangen Technologies, Beijing, China), and cDNA was synthesized with the HyperScript III RT Super Mix for qPCR with gDNA Remover kit (EnzyArtisan Biotech Co., Ltd., Eveleigh, NSW, Australia). The qRT-PCR was conducted on a real-time PCR detection system (CFX96) using 2 × S6 Universal SYBR qPCR Mix (EnzyArtisan Biotech Co., Ltd., Eveleigh, NSW, Australia), and each sample had three replications. The 2−ΔΔCT method was used to calculate the relative expressions of target genes [44]. To standardize the relative expression, the actin gene (Actin Actin-depolymerizing factor 7 [T. aestivum]) was employed as the housekeeping gene [12]. The qRT-PCR primers were designed using Primer 5.0 software [45], and all primer sequences are provided in Table 2.

2.6. VIGS-Induced Transcription Factor CA12g04950 Silencing

VIGS silencing was conducted following the method described by Shi et al. [46]. CA12g04950 was amplified using P505 high-fidelity enzyme (Vazyme Biotech Co., Ltd., Nanjing, China), and the target bands were purified using gel recovery kit (DP204, Tiangen, Beijing, China). The pTRV vector was digested with EcoR1 and BamH1, and then transformed into E. coli DH5α (Beijing TransGen Biotech Co., Ltd., Beijing, China) by ligating the TRV2 vector with c112 ligase (Vazyme). The recombinant bacterial solution was spread on LB solid medium (KanR, Yuanye, Shanghai, China), and a single colony was selected after 24 h of incubation for PCR identification and verification.
After successful identification, the viral plasmid was extracted and transformed into the Agrobacterium strain GV3101 to obtain the recombinant vector. Monoclonal Agrobacterium was confirmed by PCR, and subsequently, the infection solution was prepared (10 mM MgCl2, 10 mM MES, 200 μM AS).
The pTRV1 was mixed with TRV2 in a 1:1 volume ratio, with an empty TRV2 vector serving as the negative control (TRV), and no TRV as the blank control (CK). Four points at the base of green maturing pepper (S1 stage) fruits were selectively marked, and 150 μL of bacterial solution was injected into each point, with three fruits serving as replicates. Fruits treated with TRV were placed in a box with agar half filled. The box was then transferred to an artificial climate chamber set at 23 °C, 70% relative humidity, and a light cycle of 16/8 h (200 μmol/m2/s).
Fresh samples were collected in a circular pattern around the infected area at 3, 6, 9, 12, and 15 days after silencing. The fresh samples were immediately frozen in liquid nitrogen and stored at −80 °C for the determination of carotenoid content, friction, enzyme activity, and gene expression.

2.7. Statistical Analysis

Dates were analyzed with SPSS 26.0 (SPSS Inc., Chicago, IL, USA) and Graphpad Prims 9.5 through one-way analysis of variance (ANOVA) followed by Dunnett’s test. Data are shown as the mean ± standard deviation (SD), and p values < 0.05 were considered statistically significant.

3. Results

3.1. Effect of Nitrogen Reduction on Pigment Content of Pepper Fruit

Figure 1A illustrates a noticeable change in fruit color 30 days after urea application (S1), reaching a vibrant red hue at 45 days (S2). Fruits subjected to nitrogen reduction exhibited a consistently darker red color across stages S1 to S3, with the order of redness being N0 > N4 > N3 > N2 > N1 (Figure 1A). Analysis of total carotenoids (Figure 1B) revealed that at stages S1 and S2, compared to N1, the total carotenoids of N0 increased by 4813.71% and 2446.23%, respectively. In contrast, at the S3 stage, the total carotenoids of N3 increased by 40.12%.
Examining the fractions of carotenoids in fruits (Figure 1C–H) unveiled that the contents of phytoene (Figure 1C), α-carotene (Figure 1D), β-carotene (Figure 1E), zeaxanthin (Figure 1G), and capsanthin (Figure 1H) were notably higher at the S3 stage, except for lutein (Figure 1F). Overall, except for N0, the contents of α-carotene and β-carotene in fruits with nitrogen reduction decreased at stages S2 and S3 compared to N1. Simultaneously, the contents of zeaxanthin and capsanthin increased. At the S2 stage, the contents of zeaxanthin and capsanthin in N4 were 86.17 μg/g-FW and 116.90 μg/g-FW, respectively, representing increases of 124.94% and 85.08% compared to N1.

3.2. Effect of Nitrogen Reduction on Carotenoid Enzyme Activity and Gene Expression in Pigment Pepper Fruit

In general, enzyme activities of PSY, LCYB, and CCS increased gradually with fruit development. Figure 2A–C illustrates that urea reduction led to an increase in PSY activity at the S3 stage and CCS activity at the S2 stage. qRT-PCR analysis (Figure 2D–I) revealed that the expressions of LCYB, CRTZ, and CCS genes increased with fruit development. There were slight decreases in GGPS and lutein gene expression with fruit development. At the S2 stage, fruits treated with urea reduction exhibited higher expression levels of GGPS, PSY, LCYB, CRTZ, and CCS genes compared to N1, and the expression of PSY, LCYB, CRTZ, and CCS in N4 increased by 5 to 25 times.

3.3. Transcription and Metabolic Analysis of Pigment Pepper Fruits Treated with Nitrogen Reduction

3.3.1. RNA Sequence Quality

Table 3 shows that a minimum of 84.01% of the reads aligned with the pepper reference genome, with Q30 scores ranging from 93.16% to 94.05%. These results indicate that the sequences are of high quality.
FPKM analysis indicated high quality and repeatability in the nine samples (Figure 3A). Principal component analysis (PCA) was used to analyze the gene expression levels of multiple samples from the same period. It was found that the differences in gene expression between N1 and N4 treatments were most significant in pepper fruits (Figure 3B). In the PCA analysis, PCA1 accounted for 79.65% and PCA2 accounted for 10.2%. The grouping in PCA suggests that the differences between N1 and N4 treatment fruits were most pronounced during the same developmental stages.
Differential genes of transcription metabolism (DEGs) were calculated based on log2 fold change (FC) ≥ 1 and false discovery rate (FDR)-corrected p-value ≤ 0.05. In the comparison of N1 vs. N4 and N3 vs. N4, there were 5321 and 4590 differentially expressed genes, respectively, which were significantly higher than that of N1 vs. N3 (Figure 4A). The Venn plot (Figure 4B) showed that there were 1522 co-expressed genes, while 778 and 3799 genes were specifically expressed in the comparison of N1 vs. N3 and N1 vs. N4. Among DEGs in the two pairwise groups, 8134–17,192 genes were successfully annotated to one of the five public databases (Figure S1).

3.3.2. Analysis of GO and KEGG Enrichment in Pigment Pepper Fruits under Nitrogen Reduction

Differentially expressed genes (DEGs) can be categorized into various metabolic pathways based on biological processes, molecular functions, and cellular components. Figure 5A presents the top 10 enriched metabolic pathways of DEGs, including porphyrin-containing compound metabolic processes, tetrapyrrole metabolic processes, pigment metabolic processes, chloroplast envelope, cell wall, external encapsulating structure, hydrolyzing O-glycosyl compounds, catalytic activity, and acting on glycosyl bonds. Notably, the most differently expressed genes were enriched in the N1 vs. N4 comparison (Table S1; Figure 5A, right column). In the pigment metabolic processes and pigment biosynthesis processes, −log10 (p-value) of DEGs were 7.61 and 6.96 in N1 vs. N4, which were much higher than those of N1 vs. N3 (0.96, 0.94) and N3 vs. N4 (4.61, 4.29). There was a similar trend in the chloroplast envelope, where carotenoid was produced (Figure 5A, right column).
In the KEGG analysis (Figure 5B), the DEGs from N1 vs. N3, N1 vs. N4, and N3 vs. N4 were annotated to 124, 124, and 126 metabolic pathways, respectively (Table S2), and the top five enriched pathways were photosynthesis-antenna proteins, porphyrin metabolism, carotenoid biosynthesis, seleno-compound metabolism, and nitrogen metabolism. log10(p-value) of porphyrin metabolism and carotenoid biosynthesis were 8.27 and 4.10 in N1 vs. N4 which were higher than that of N1 vs. N3 (2.69, 2.02) (Figure 5B, right column). Specifically, the most differently expressed genes were enriched in N1 vs. N4, with 18 genes involved in carotenoid metabolism and 16 genes involved in nitrogen metabolism. Furthermore, carotenoid content genes Z-ISO and AO 4 were upregulated, while the SDRs were downregulated under both the N3 and N4 treatments. Nitrogen metabolism genes NRT 2.5 and NRT 2.7 were upregulated in N3 and N4 treatments. Z-ISO participated in regulating the biosynthesis of strigolactone (SL) [47], and SL is an important mediator involved in nitrogen regulation within the plant [48]. The NRTs were involved in the absorption of NO3− under low nitrogen conditions [49].

3.3.3. Analysis of Differentially Expressed Genes (DEGs) in Carotenoid Metabolism of Pigment Pepper Fruits under Nitrogen Reduction

To validate the precision of the RNA-seq data, the expression of nine carotenoid genes was measured using qRT–PCR (Figure 6). The results demonstrated that the transcription levels of these nine carotenoid genes were consistent with the RNA-seq expression, indicating the reliability of the transcription data in this study.

3.3.4. Analysis of Transcription Factors in Pepper Fruits under Nitrogen Reduction

The analysis revealed the presence of numerous differential transcription factor families, including ERF, bHLH, MYB, C2H2, and NAC (Figure 7), under different urea dosages. The NAC family, crucially related to pigment metabolism [50,51], was selected as the focus of subsequent research. The analysis of NAC transcription factors in N1 vs. N4 showed that 12 NAC members, with the names CA02g22730, CA03g21000, CA04g14740, CA05g04410, CA06g11310, CA06g21590, CA07g21320, and CA08g15460, exhibited significant changes in expression levels. Correlation analysis (Figure 7B,C) demonstrated a high correlation between the expression of carotenoid genes such as ABA 8′, BCH, PSY, SDRS, CCS, and VDE and NAC transcription factors. Among these genes, the expression of CA12g04950 in N4 was significantly higher than that in N1.

3.4. Regulation of CA12g04950 on Pepper Fruits’ Carotenoids

CA12g04950 was silenced using the VIGS technique (Figure 8). As shown in Figure 8A, the fruit color of the control group (CK) and the empty vector (TRV) began to turn red at 9 days after silencing, while the fruits of TRV-CA12g04950 began to change color at 12 days. By the 15th day of silencing, the fruits of CK and TRV were fully red, whereas only part of the fruit surface of TRV-CA12g04950 became red. At day 6, the total chlorophyll content of TRV-CA12g04950 fruits (Figure 8C,D) increased by 13.03% and 24.75% compared to CK and TRV. At day 12, compared to CK and TRV, the chlorophyll content of TRV−CA12g04950 increased by 74.67% and 42.33%, respectively, while the total carotenoid content decreased by 87.31% and 69.15%, respectively. qRT−PCR indicated that the expression of the CA12g04950 gene in CK and TRV increased initially and then decreased (Figure 8E). In contrast, the expression level of CA12g04950 in silenced fruit was significantly lower than that in CK and TRV from day 3 to day 12. Specifically, at day 9, the CA12g04950 expression in silenced fruits decreased by 58.17% and 42.39%, respectively, compared to CK and TRV, but by day 15, the CA12g04950 expression increased by 51.3% and 39.2%, respectively. These results suggest that silencing CA12g04950 can inhibit the carotenoid content within a certain time.
Determining the content of carotenoid components in the fruit after silencing (Figure 9A–D) revealed significant differences. The content of capsaicinoids in the sixth-day silenced fruit was 2.02 mg/kg·FW, which marked a decrease of 96.6% and 83.3% compared to the CK and TRV groups (Figure 9D), respectively. Similar trends were observed in the content of other carotenoids compared to the control and empty vector groups. At 12D, the content of capsaicinoids was 93.08% and 85.44% lower than that of the CK and TRV groups, respectively (Figure 9D), consistent with the phenotypic results. These results were consistent with the fruit phenotype. After 6 days of silencing, the enzyme activities of PSY, LCYB, and CCS in TRV-CA12g04950 were lower than that in CK and TRV; however, the CCS enzyme activity increased rapidly after 12 days.
qRT-PCR analysis indicated that within 9 days of silencing, the transcription levels of PSY, LCYB, CRTZ, ZEP, and CCS genes were significantly suppressed (Figure 10A–F). After 6 days of silencing, PSY expression in TRV-CA12g04950 fruits decreased by 93.74% and 85.67% compared to CK and TRV, and by 41.38% and 62.49% after 12 days. Within 6 days, the transcription levels of ZEP decreased by 94.29% and 93.51% compared to CK and TRV. Throughout the silencing period, the expression level of CCS was significantly reduced. From days 6 to 15, the CCS expression in TRV-CA12g04950 was reduced by 97.29% to 24.18% compared to CK and by 97.68% to 25.57% compared to TRV.

4. Discussion

4.1. Nitrogen Fertilizer and Pigment

Color serves as a critical quality indicator for crops, with plant color being primarily determined by flavonoids, chlorophyll, and carotenoids [52]. Nitrogen emerges as a crucial factor influencing plant color, as its concentration affects metabolic processes and chlorophyll content [53]. Notably, the contents of chlorophyll, lutein, β-carotene, neoxanthin, purple xanthin, and anther xanthin in sugar beets increased significantly as the ratio of NH4+ to NO3− decreased [54,55].
Chang [56] observed that when nitrogen fertilizer application decreased from 51 kg/hm2 to 38 kg/hm2, the carotenoid content in pepper fruits increased. However, when the nitrogen application rate decreased to 25 kg/hm2, the carotenoid content decreased. In this experiment, it was observed that during stages S1 and S2, as the urea application rate decreased from 750 kg/hm2 to 0 kg/hm2, the carotenoid content in pepper fruits gradually increased. Remarkably, during the red ripening period (S3), N3 exhibited the highest carotenoid content (Figure 1B). This aligns with the findings of Kopsell [54] in kale, where carotenoid content increased with reduced ammonium nitrogen application. However, some discrepancies exist with the study by Francisco M. del Amor and T. Casey Barickman [55,57], who found that carotenoid content increased when ammonium nitrogen application was reduced early in development but increased at the peak of reproductive growth [58].
Barickman’s [55] results demonstrated that as the ammonium nitrogen usage decreased from 100% to 50%, the carotenoid content increased in Swiss chard leaf tissue. Furthermore, carotenoid content rose with decreasing nitrogen application during the full reproductive stage. Enzyme activities of PSY, LCYB, and CCS in this study generally exhibited a trend of gradual increase with decreasing nitrogen application. Specifically, the enzyme activity of PSY increased significantly during the S3 period, and the expression of PSY also displayed significant upregulation during the same period. Zhang et al.’s [59] study demonstrates that N250 (250 kg N/hm2) is beneficial for capsanthin accumulation and secondary metabolites in pepper fruits actively respond to nitrogen supply, leading to changes in the composition, content, and enzyme activity of pigment-related substances in the fruits. This study’s conclusion is similar to Zhang’s [59]. According to GO and KEGG (Figure 5), the differential expression of genes related to secondary metabolites at the same developmental stage, N1 vs. N4, is significantly higher compared to N1 vs. N3h and N3 vs. N4. Except for GGPSS, other key genes of carotenoid content exhibited an increase in expression levels with the reduction in nitrogen application. GGPP, a precursor for the synthesis of numerous plant olefins and belonging to the upstream carotenoid metabolic pathway [12], did not show a significant effect of nitrogen. This further confirms that plants respond positively to environmental changes by regulating transcript levels to adapt to varying nutrient levels [60].

4.2. Transcription Factors and Carotenoid

Transcription factors have garnered significant attention in recent years. In addition to NAC and ERF, a notable change in urea dosage led to significant differences in the expression of transcription factor families such as bHLH, MYB, and C2H2 (Figure 7A). Recent research has emphasized the regulatory role of transcription factors in determining the color of plant fruits. Jiali Song [61] conducted an analysis of the ERF family in peppers and identified CaERF82, CaERF97, CaERF66, CaERF107, and CaERF101 as participants in the carotenoid synthesis of pepper peel. Kang et al. [62] indicated that the early synthesis of anthocyanins in tomato fruit peels is regulated by the complex of MYB, bHLH, and WD40. Further findings revealed that bHLH and MYB regulated the expression of flavonoid biosynthesis genes, subsequently influencing fruit color [63].
Pearson correlation analysis between the expressions of upregulated transcription factors and carotenoid genes in N1 and N4 treatments revealed that NAC transcription factors significantly correlated with carotenoid gene expression in N4. Among them, CA12g04950, belonging to the NAC TF family, exhibited a significant positive correlation with 11 carotenoid genes. In this study, the expressions of CA12g04950 increased significantly in N4 (control) and decreased in the silenced fruits of N4. On the sixth day, the expression level of PSY decreased by 93.74% in silenced fruits compared to the control, while the expression level of CCS decreased by 97.29% compared to the control. It can be inferred that CA12g04950 is involved in the regulation of pepper carotenoid synthesis under low nitrogen conditions.
The NAC family, as the most abundant and functionally comprehensive transcription factors in plants, plays a crucial role in regulating immune signal transmission. NAC proteins have been shown to modulate plant hormones, including GA, ABA, and IAA, to respond to various stresses [51,64,65,66]. Zhang et al. [67] discovered positive regulations of abiotic stresses such as salt, cold, and heat on the expression level of CaNAC035.
Gao et al. [50] silenced Solyc07g063420, a member of the NAC family, and found that the maturation of tomato fruit was delayed, with fruit softening and chlorophyll degradation being delayed, leading to a decrease in lycopene accumulation. Knocking out NOR-like1 by CRISPR/Cas9 led to a decrease in ethylene content in tomato fruit, subsequently changing the expressions of carotenoid genes such as DXS, GGPP2, PSY1, and PSY, thereby affecting fruit color [68,69]. Research by Feng et al. [70] demonstrated that silencing the SINACs of tomato led to a decrease in ethylene content and lycopene content, resulting in yellow fruit. The expressions of carotenoid genes (PSY1, DXS2, SGR1, CrtR-b2) in the SINAC9 mutant changed significantly, affecting the transformation process of chloroplasts and chromosomes. Many NAC family members are thus involved in the regulation of carotenoid genes.
PSY, as the upstream gene of the carotenoid metabolic pathway, is extensively researched. Liu et al. [30] found that the most important function of PSY1 in peppers was to provide precursors for downstream carotenoid synthesis by regulating phytoene accumulation, rather than directly participating in pigment synthesis. Yeast one-hybrid experiments demonstrated that SINAC9 interacted with PSY1 at the protein level, affecting the chlorophyll degradation and carotenoid synthesis of tomatoes [21,70]. This study similarly found that silencing CA12g04950 had markedly different effects on PSY and CCS. Silencing 78 resulted in a significant decrease in the expression level of CCS compared to PSY, similar to the results of the study mentioned above. This suggests that PSY may indirectly regulate carotenoid synthesis in peppers.

5. Conclusions

This study utilized transcriptomic analysis to investigate the differential expression of carotenoid-related genes and transcription factors in pepper fruits under varying nitrogen conditions. Pearson correlation analysis revealed a positive correlation between CA12g04950 and differentially expressed carotenoid genes under the N4 treatment. Then, silencing CA12g04950 under N4 treatment significantly decreased the contents of phytoene, β-carotenoids, zeaxanthin, and capsanthin, as well as the activity of carotenoid enzymes. These findings suggest that under nitrogen deficiency conditions, CA12g04950 may directly activate the biosynthetic genes of carotenoids in peppers, thereby increasing the levels of carotenoids. In this process, CA12g04950 simultaneously acts on multiple carotenoid biosynthetic genes, with the downstream CCS showing a more pronounced effect (Figure 11). Future research will focus on verifying the function of CA12g04950 in the coloration of pepper fruits through protein interaction, overexpression, and silencing techniques.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/agriculture14040521/s1, Figure S1: Functional annotations of transcripts to public databases; Table S1: Analysis of GO enrichment for DEGs; Table S2: KEGG pathways of differentially expressed genes.

Author Contributions

All authors contributed to the discussion. L.S. and X.J. conceived and designed the experiment; L.S. performed experiments and data analyses and wrote the manuscript; X.J. and S.U. were involved in revising the manuscript; L.S. and X.J. collected the literature and revised the manuscript; X.J. and S.U. checked the language. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant numbers 32060676 and 31860548.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Materials; further inquiries can be directed to the corresponding authors. The raw data supporting the conclusions of this article will be made available by the authors upon request. This study not involving humans or animals.

Acknowledgments

We thank Yanjun Li and Xiaohu Ma (Shihezi University) for helping with the TRV1 and TRV2 vectors.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Barik, S.; Ponnam, N.; Reddy, A.C.; DC, L.R.; Saha, K.; Acharya, G.C.; Reddy, M. Breeding peppers for industrial uses: Progress and prospects. Ind. Crops Prod. 2022, 178, 114626. [Google Scholar] [CrossRef]
  2. Baenas, N.; Belović, M.; Ilic, N.; Moreno, D.A.; García-Viguera, C. Industrial use of pepper (Capsicum annum L.) derived products: Technological benefits and biological advantages. Food Chem. 2019, 274, 872–885. [Google Scholar] [CrossRef]
  3. Gómez-García, M.; Ochoa-Alejo, N. Biochemistry and Molecular Biology of Carotenoid Biosynthesis in Chili Peppers (Capsicum spp.). Int. J. Mol. Sci. 2013, 14, 19025–19053. [Google Scholar] [CrossRef] [PubMed]
  4. Lee, S.Y.; Jang, S.J.; Jeong, H.B.; Lee, S.Y.; Venkatesh, J.; Lee, J.H.; Kwon, J.K.; Kang, B.C. A mutation in Zeaxanthin epoxidase contributes to orange coloration and alters carotenoid contents in pepper fruit (Capsicum annuum). Plant J. 2021, 106, 1692–1707. [Google Scholar] [CrossRef] [PubMed]
  5. Villa-Rivera, M.G.; Ochoa-Alejo, N. Chili Pepper Carotenoids: Nutraceutical Properties and Mechanisms of Action. Molecules 2020, 25, 5573. [Google Scholar] [CrossRef] [PubMed]
  6. Giuffrida, D.; Dugo, P.; Torre, G.; Bignardi, C.; Cavazza, A.; Corradini, C.; Dugo, G. Characterization of 12 Capsicum varieties by evaluation of their carotenoid profile and pungency determination. Food Chem. 2013, 140, 794–802. [Google Scholar] [CrossRef] [PubMed]
  7. Wahyuni, Y.; Ballester, A.R.; Tikunov, Y.; de Vos, R.C.; Pelgrom, K.T.; Maharijaya, A.; Sudarmonowati, E.; Bino, R.J.; Bovy, A.G. Metabolomics and molecular marker analysis to explore pepper (Capsicum sp.) biodiversity. Metabolomics 2013, 9, 130–144. [Google Scholar] [CrossRef] [PubMed]
  8. Mi, S.; Yu, W.; Li, J.; Liu, M.; Sang, Y.; Wang, X. Characterization and discrimination of chilli peppers based on multi-element and non-targeted metabolomics analysis. LWT 2020, 131, 109742. [Google Scholar] [CrossRef]
  9. Zoccali, M.; Giuffrida, D.; Salafia, F.; Rigano, F.; Dugo, P.; Casale, M.; Mondello, L. Apocarotenoids profiling in different Capsicum species. Food Chem. 2021, 334, 127595. [Google Scholar] [CrossRef]
  10. Mi, S.; Zhang, X.; Wang, Y.; Zheng, M.; Zhao, J.; Gong, H.; Wang, X. Effect of different genotypes on the fruit volatile profiles, flavonoid composition and antioxidant activities of chilli peppers. Food Chem. 2022, 374, 131751. [Google Scholar] [CrossRef]
  11. Feng, P.; Wang, Y.; Wen, J.; Ren, Y.; Zhong, Q.; Li, Q. Cloning and Analysis of Expression of Genes Related to Carotenoid Metabolism in Different Fruit Color Mutants of Pepper (Capsicum annuum L.). Genes 2024, 15, 315. [Google Scholar] [CrossRef]
  12. Lv, J.; Zhang, R.; Mo, Y.; Zhou, H.; Li, M.; Wu, R.; Cheng, H.; Zhang, M.; Wang, H.; Hua, W.; et al. Integrative Metabolome and Transcriptome Analyses Provide Insights into Carotenoid Variation in Different-Colored Peppers. Int. J. Mol. Sci. 2023, 24, 16563. [Google Scholar] [CrossRef]
  13. Sun, T.; Yuan, H.; Cao, H.; Yazdani, M.; Tadmor, Y.; Li, L. Carotenoid Metabolism in Plants: The Role of Plastids. Mol. Plant 2018, 11, 58–74. [Google Scholar] [CrossRef]
  14. Giuliano, G. Provitamin A biofortification of crop plants: A gold rush with many miners. Curr. Opin. Biotechnol. 2017, 44, 169–180. [Google Scholar] [CrossRef]
  15. Bai, C.; Capell, T.; Berman, J.; Medina, V.; Sandmann, G.; Christou, P.; Zhu, C. Bottlenecks in carotenoid biosynthesis and accumulation in rice endosperm are influenced by the precursor-product balance. Plant Biotechnol. J. 2016, 14, 195–205. [Google Scholar] [CrossRef] [PubMed]
  16. Tian, S.L.; Li, L.; Chai, W.G.; Shah, S.N.; Gong, Z.H. Effects of silencing key genes in the capsanthin biosynthetic pathway on fruit color of detached pepper fruits. BMC Plant Biol. 2014, 14, 314. [Google Scholar] [CrossRef] [PubMed]
  17. Tian, S.L.; Li, L.; Shah, S.N.M.; Gong, Z.H. The relationship between red fruit colour formation and key genes of capsanthin biosynthesis pathway in Capsicum annuum. Biol. Plant. 2015, 59, 507–513. [Google Scholar] [CrossRef]
  18. Welsch, R.; Medina, J.; Giuliano, G.; Beyer, P.; Von Lintig, J. Structural and functional characterization of the phytoene synthase promoter from Arabidopsis thaliana. Planta 2003, 216, 523–534. [Google Scholar] [CrossRef]
  19. Lu, S.; Ye, J.; Zhu, K.; Zhang, Y.; Zhang, M.; Xu, Q.; Deng, X. A fruit ripening-associated transcription factor CsMADS5 positively regulates carotenoid biosynthesis in citrus. J. Exp. Bot. 2021, 72, 3028–3043. [Google Scholar] [CrossRef] [PubMed]
  20. Sun, Q.; Zhang, N.; Wang, J.; Zhang, H.; Li, D.; Shi, J.; Li, R.; Weeda, S.; Zhao, B.; Ren, S.; et al. Melatonin promotes ripening and improves quality of tomato fruit during postharvest life. J. Exp. Bot. 2015, 66, 657–668. [Google Scholar] [CrossRef]
  21. Yuan, Y.; Ren, S.; Liu, X.; Su, L.; Wu, Y.; Zhang, W.; Li, Y.; Jiang, Y.; Wang, H.; Fu, R.; et al. SlWRKY35 positively regulates carotenoid biosynthesis by activating the MEP pathway in tomato fruit. New Phytol. 2022, 234, 164–178. [Google Scholar] [CrossRef]
  22. Ma, N.; Feng, H.; Meng, X.; Li, D.; Yang, D.; Wu, C.; Meng, Q. Overexpression of tomato SlNAC1transcription factor alters fruit pigmentation and softening. BMC Plant Biol. 2014, 14, 351. [Google Scholar] [CrossRef] [PubMed]
  23. Song, J.; Sun, B.; Chen, C.; Ning, Z.; Zhang, S.; Cai, Y.; Zheng, X.; Cao, B.; Chen, G.; Jin, D.; et al. An R-R-type MYB transcription factor promotes non-climacteric pepper fruit carotenoid pigment biosynthesis. Plant J. Cell Mol. Biol. 2023, 115, 724–741. [Google Scholar] [CrossRef] [PubMed]
  24. Zhang, J.; Lv, J.; Xie, J.; Gan, Y.; Coulter, J.A.; Yu, J.; Li, J.; Wang, J.; Zhang, X. Nitrogen Source Affects the Composition of Metabolites in Pepper (Capsicum annuum L.) and Regulates the Synthesis of Capsaicinoids through the GOGAT-GS Pathway. Foods 2020, 9, 150. [Google Scholar] [CrossRef] [PubMed]
  25. Lamrani, Z.; Belakbir, A.; Ruiz, J.M.; Ragala, L.; López-Cantarero, I.; Romero, L. Influence of nitrogen, phosphorus, and potassium on pigment concentration in cucumber leaves. Commun. Soil Sci. Plant Anal. 2008, 27, 1001–1012. [Google Scholar] [CrossRef]
  26. Berova, M.; Karanatsidis, G. Influence of bio-fertilizer, produced by Lumbricus rubellus on growth, leaf gas-exchange and photosynthetic pigment content of pepper plants (Capsicum annuum L.). Acta Hortic. 2009, 830, 447–452. [Google Scholar] [CrossRef]
  27. Johnson, C.D.; Decoteau, D.R. Nitrogen and Potassium Fertility Affects Jalapeño Pepper Plant Growth, Pod Yield, and Pungency. HortScience 1996, 31, 1119–1123. [Google Scholar] [CrossRef]
  28. Liu, M.; Yang, S.Q.; Fu, X. Effects of Water Nitrogen Interaction on Growth and Yield of Lycium Barbarum under Brackish Water Drip Irrigation. Water Sav. Irrig. 2018, 33–37. [Google Scholar] [CrossRef]
  29. Kuscu, H.; Turhan, A.; Ozmen, N.; Aydinol, P.; Demir, A.O. Optimizing levels of water and nitrogen applied through drip irrigation for yield, quality, and water productivity of processing tomato (Lycopersicon esculentum Mill.). Hortic. Environ. Biotechnol. 2014, 55, 103–114. [Google Scholar] [CrossRef]
  30. Liu, Y.; Lv, J.; Liu, Z.; Wang, J.; Yang, B.; Chen, W.; Ou, L.; Dai, X.; Zhang, Z.; Zou, X. Integrative analysis of metabolome and transcriptome reveals the mechanism of color formation in pepper fruit (Capsicum annuum L.). Food Chem. 2020, 306, 125629. [Google Scholar] [CrossRef]
  31. Burgos, G.; Salas, E.; Amoros, W.; Auqui, M.; Muñoa, L.; Kimura, M.; Bonierbale, M. Total and individual carotenoid profiles in Solanum phureja of cultivated potatoes: I. Concentrations and relationships as determined by spectrophotometry and HPLC. J. Food Compos. Anal. 2009, 22, 503–508. [Google Scholar] [CrossRef]
  32. Dou, J.-L.; Yuan, P.-L.; Zhao, S.-J.; He, N.; Zhu, H.-J.; Gao, L.; Ji, W.-L.; Lu, X.-Q.; Liu, W.-G. Effect of ploidy level on expression of lycopene biosynthesis genes and accumulation of phytohormones during watermelon (Citrullus lanatus) fruit development and ripening. J. Integr. Agric. 2017, 16, 1956–1967. [Google Scholar] [CrossRef]
  33. Rodríguez-Rodríguez, E.; Sánchez-Prieto, M.; Olmedilla-Alonso, B. Assessment of carotenoid concentrations in red peppers (Capsicum annuum) under domestic refrigeration for three weeks as determined by HPLC-DAD. Food Chem. X 2020, 6, 100092. [Google Scholar] [CrossRef] [PubMed]
  34. Kim, D.; Paggi, J.M.; Park, C.; Bennett, C.; Salzberg, S.L. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 2019, 37, 907–915. [Google Scholar] [CrossRef] [PubMed]
  35. Trapnell, C.; Williams, B.A.; Pertea, G.; Mortazavi, A.; Kwan, G.; van Baren, M.J.; Salzberg, S.L.; Wold, B.J.; Pachter, L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 2010, 28, 511–515. [Google Scholar] [CrossRef]
  36. Deng, Y.Y.; Li, J.Q.; Wu, S.F.; Zhu, Y.; Chen, Y.W.; He, F.C. Integrated nr Database in Protein Annotation System and Its Localization. Comput. Eng. 2006, 32, 71–72. [Google Scholar]
  37. Kanehisa, M.; Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  38. Boeckmann, B.; Bairoch, A.; Apweiler, R.; Blatter, M.C.; Estreicher, A.; Gasteiger, E.; Martin, M.J.; Michoud, K.; O’Donovan, C.; Phan, I.; et al. The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res. 2003, 31, 365–370. [Google Scholar] [CrossRef]
  39. Ashburner, M.; Ball, C.A.; Blake, J.A.; Botstein, D.; Butler, H.; Cherry, J.M.; Davis, A.P.; Dolinski, K.; Dwight, S.S.; Eppig, J.T.; et al. Gene ontology: Tool for the unification of biology. The Gene Ontology Consortium. Nat. Genet. 2000, 25, 25–29. [Google Scholar] [CrossRef]
  40. Huerta-Cepas, J.; Forslund, K.; Coelho, L.P.; Szklarczyk, D.; Jensen, L.J.; von Mering, C.; Bork, P. Fast Genome-Wide Functional Annotation through Orthology Assignment by eggNOG-Mapper. Mol. Biol. Evol. 2017, 34, 2115–2122. [Google Scholar] [CrossRef] [PubMed]
  41. Young, M.D.; Wakefield, M.J.; Smyth, G.K.; Oshlack, A. Gene ontology analysis for RNA-seq: Accounting for selection bias. Genome Biol. 2010, 11, R14. [Google Scholar] [CrossRef]
  42. Kanehisa, M.; Sato, Y.; Kawashima, M. KEGG mapping tools for uncovering hidden features in biological data. Protein Sci. 2022, 31, 47–53. [Google Scholar] [CrossRef] [PubMed]
  43. Wang, L.; Feng, Z.; Wang, X.; Wang, X.; Zhang, X. DEGseq: An R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics 2010, 26, 136–138. [Google Scholar] [CrossRef] [PubMed]
  44. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  45. Ren, L.; Zhu, B.Q.; Zhang, Y.B.; Wang, H.Y.; Li, C.Y.; Su, Y.; Ba, C.F. The research of applying primer premier 5.0 to design PCR primer. J. Jinzhou Med. Coll. 2004, 25, 43–46. [Google Scholar]
  46. Shi, G.; Hao, M.; Tian, B.; Cao, G.; Wei, F.; Xie, Z. A Methodological Advance of Tobacco Rattle Virus-Induced Gene Silencing for Functional Genomics in Plants. Front. Plant Sci. 2021, 12, 671091. [Google Scholar] [CrossRef]
  47. Liu, L.; Ren, M.; Peng, P.; Chun, Y.; Li, L.; Zhao, J.; Fang, J.; Peng, L.; Yan, J.; Chu, J.; et al. MIT1, encoding a 15-cis-zeta-carotene isomerase, regulates tiller number and stature in rice. J. Genet. Genom. 2021, 48, 88–91. [Google Scholar] [CrossRef]
  48. Ito, S.; Ito, K.; Abeta, N.; Takahashi, R.; Sasaki, Y.; Yajima, S. Effects of strigolactone signaling on Arabidopsis growth under nitrogen deficient stress condition. Plant Signal. Behav. 2016, 11, e1126031. [Google Scholar] [CrossRef] [PubMed]
  49. Li, Z.; Wu, X.N.; Jaquot, A.; Chaput, V.; Adamo, M.; Neuhauser, B.; Straub, T.; Lejay, L.; Schulze, W.X. Phosphoregulation in the N-terminus of NRT2.1 affects nitrate uptake by controlling the interaction of NRT2.1 with NAR2.1 and kinase HPCAL1. J. Exp. Bot. 2023. [Google Scholar] [CrossRef]
  50. Gao, Y.; Wei, W.; Zhao, X.; Tan, X.; Fan, Z.; Zhang, Y.; Jing, Y.; Meng, L.; Zhu, B.; Zhu, H.; et al. A NAC transcription factor, NOR-like1, is a new positive regulator of tomato fruit ripening. Hortic. Res. 2018, 5, 75. [Google Scholar] [CrossRef]
  51. Martin-Pizarro, C.; Vallarino, J.G.; Osorio, S.; Meco, V.; Urrutia, M.; Pillet, J.; Casanal, A.; Merchante, C.; Amaya, I.; Willmitzer, L.; et al. The NAC transcription factor FaRIF controls fruit ripening in strawberry. Plant Cell 2021, 33, 1574–1593. [Google Scholar] [CrossRef]
  52. Stournaras, K.E.; Lo, E.; Böhning-Gaese, K.; Cazetta, E.; Matthias Dehling, D.; Schleuning, M.; Stoddard, M.C.; Donoghue, M.J.; Prum, R.O.; Martin Schaefer, H. How colorful are fruits? Limited color diversity in fleshy fruits on local and global scales. New Phytol. 2013, 198, 617–629. [Google Scholar] [CrossRef]
  53. Wang, S.; Guan, K.; Wang, Z.; Ainsworth, E.A.; Zheng, T.; Townsend, P.A.; Li, K.; Moller, C.; Wu, G.; Jiang, C. Unique contributions of chlorophyll and nitrogen to predict crop photosynthetic capacity from leaf spectroscopy. J. Exp. Bot. 2021, 72, 341–354. [Google Scholar] [CrossRef] [PubMed]
  54. Kopsell, D.A.; Kopsell, D.E.; Curran-Celentano, J. Carotenoid pigments in kale are influenced by nitrogen concentration and form. J. Sci. Food Agric. 2007, 87, 900–907. [Google Scholar] [CrossRef]
  55. Barickman, T.C.; Kopsell, D.A. Nitrogen form and ratio impact Swiss chard (Beta vulgaris subsp. cicla) shoot tissue carotenoid and chlorophyll concentrations. Sci. Hortic. 2016, 204, 99–105. [Google Scholar] [CrossRef]
  56. Chang, X.; Wei, Y.; Li, A.; Chen, X.; Wei, X.; Duan, X. Evaluation on Several Germplasm Resources of Dried Hot Pepper and Correlative Analysis. Shandong Agric. Sci. 2014, 46, 28–30. [Google Scholar] [CrossRef]
  57. del Amor, F.M.; Cuadra-Crespo, P.; Varó, P.; Gómez, M.C. Influence of foliar urea on the antioxidant response and fruit color of sweet pepper under limited N supply. J. Sci. Food Agric. 2008, 89, 504–510. [Google Scholar] [CrossRef]
  58. Davarpanah, S.; Tehranifar, A.; Davarynejad, G.; Aran, M.; Abadía, J.; Khorassani, R. Effects of Foliar Nano-nitrogen and Urea Fertilizers on the Physical and Chemical Properties of Pomegranate (Punica granatum cv. Ardestani) Fruits. HortScience 2017, 52, 288–294. [Google Scholar] [CrossRef]
  59. Zhang, L.; Zhang, F.; He, X.; Dong, Y.; Sun, K.; Liu, S.; Wang, X.; Yang, H.; Zhang, W.; Lakshmanan, P.; et al. Comparative metabolomics reveals complex metabolic shifts associated with nitrogen-induced color development in mature pepper fruit. Front. Plant Sci. 2024, 15, 1319680. [Google Scholar] [CrossRef]
  60. Brumbarova, T.; Ivanov, R. The Nutrient Response Transcriptional Regulome of Arabidopsis. iScience 2019, 19, 358–368. [Google Scholar] [CrossRef]
  61. Song, J.L.; Chen, C.M.; Zhang, S.L.; Wang, J.T.; Huang, Z.B.; Chen, M.X.; Cao, B.H.; Zhu, Z.S.; Lei, J.J. Systematic analysis of the ERF transcription factor family: Identification of regulatory factors involved in the regulation of species-specific metabolites. BMC Genom. 2020, 21, 573. [Google Scholar] [CrossRef]
  62. Kang, S.-I.; Rahim, M.A.; Afrin, K.S.; Jung, H.-J.; Kim, H.-T.; Park, J.-I.; Nou, I.-S. Expression of anthocyanin biosynthesis-related genes reflects the peel color in purple tomato. Hortic. Environ. Biotechnol. 2018, 59, 435–445. [Google Scholar] [CrossRef]
  63. Duan, L.; Zhang, B.; Dai, G.; He, X.; Zhou, X.; Huang, T.; Liang, X.; Zhao, J.; Qin, K. Integrated Analysis of Transcriptome and Metabolome Reveals New Insights into the Molecular Mechanism Underlying the Color Differences in Wolfberry (Lycium barbarum). Agronomy 2023, 13, 1926. [Google Scholar] [CrossRef]
  64. Peng, Y.; van Wersch, R.; Zhang, Y. Convergent and Divergent Signaling in PAMP-Triggered Immunity and Effector-Triggered Immunity. Mol. Plant-Microbe Interact. MPMI 2018, 31, 403–409. [Google Scholar] [CrossRef]
  65. Wang, J.; Zheng, C.; Shao, X.; Hu, Z.; Li, J.; Wang, P.; Wang, A.; Yu, J.; Shi, K. Transcriptomic and genetic approaches reveal an essential role of the NAC transcription factor SlNAP1 in the growth and defense response of tomato. Hortic. Res. 2020, 7, 209. [Google Scholar] [CrossRef]
  66. Moyano, E.; Martínez-Rivas, F.J.; Blanco-Portales, R.; Molina-Hidalgo, F.J.; Ric-Varas, P.; Matas-Arroyo, A.J.; Caballero, J.L.; Muñoz-Blanco, J.; Rodríguez-Franco, A. Genome-wide analysis of the NAC transcription factor family and their expression during the development and ripening of the Fragaria × ananassa fruits. PLoS ONE 2018, 13, e0196953. [Google Scholar] [CrossRef] [PubMed]
  67. Zhang, H.; Ma, F.; Wang, X.; Liu, S.; Saeed, U.H.; Hou, X.; Zhang, Y.; Luo, D.; Meng, Y.; Zhang, W.; et al. Molecular and Functional Characterization of CaNAC035, an NAC Transcription Factor from Pepper (Capsicum annuum L.). Front. Plant Sci. 2020, 11, 14. [Google Scholar] [CrossRef]
  68. Zhu, M.; Chen, G.; Zhou, S.; Tu, Y.; Wang, Y.; Dong, T.; Hu, Z. A New Tomato NAC (NAM/ATAF1/2/CUC2) Transcription Factor, SlNAC4, Functions as a Positive Regulator of Fruit Ripening and Carotenoid Accumulation. Plant Cell Physiol. 2013, 55, 119–135. [Google Scholar] [CrossRef] [PubMed]
  69. Kou, X.; Zhao, Y.; Wu, C.; Jiang, B.; Zhang, Z.; Rathbun, J.R.; He, Y.; Xue, Z. SNAC4 and SNAC9 transcription factors show contrasting effects on tomato carotenoids biosynthesis and softening. Postharvest Biol. Technol. 2018, 144, 9–19. [Google Scholar] [CrossRef]
  70. Feng, Y.; Kou, X.; Yuan, S.; Wu, C.; Zhao, X.; Xue, Z.; Li, Q.; Huang, Z.; Sun, Y. CRISPR/Cas9-mediated SNAC9 mutants reveal the positive regulation of tomato ripening by SNAC9 and the mechanism of carotenoid metabolism regulation. Hortic. Res. 2023, 10, uhad019. [Google Scholar] [CrossRef]
Agriculture 14 00521 g001
Figure 1. Phenotypes and pigment content of pepper fruits at different stages with nitrogen reduction. (A) Phenotypes of pigment pepper fruits. (B) Total carotenoids of pigment pepper fruits. (C–H) Contents of six carotenoids fractions of pigment pepper fruits. Data are presented as mean ± SD. Different letters (a, b, c, d, and e) indicate the significance of differences between nitrogen treatments (p ≤ 0.05).
Figure 1. Phenotypes and pigment content of pepper fruits at different stages with nitrogen reduction. (A) Phenotypes of pigment pepper fruits. (B) Total carotenoids of pigment pepper fruits. (C–H) Contents of six carotenoids fractions of pigment pepper fruits. Data are presented as mean ± SD. Different letters (a, b, c, d, and e) indicate the significance of differences between nitrogen treatments (p ≤ 0.05).
Agriculture 14 00521 g001
Agriculture 14 00521 g002
Figure 2. Enzyme activities of PSY (A), LCYB (B), CCS (C), and carotenoid gene expression (D–I) of pigment pepper fruits at different developmental stages under varied urea dosage. Different letters (a, b, c, d, and e) indicate the significance of differences between nitrogen treatments (p ≤ 0.05).
Figure 2. Enzyme activities of PSY (A), LCYB (B), CCS (C), and carotenoid gene expression (D–I) of pigment pepper fruits at different developmental stages under varied urea dosage. Different letters (a, b, c, d, and e) indicate the significance of differences between nitrogen treatments (p ≤ 0.05).
Agriculture 14 00521 g002
Agriculture 14 00521 g003
Figure 3. Distribution of kilobase per million exon kilobase fragment mapping (FPKM) in rNA-seq of pigment pepper fruits under urea reduction. (A) Violin plot; (B) FPKM-based principal component analysis.
Figure 3. Distribution of kilobase per million exon kilobase fragment mapping (FPKM) in rNA-seq of pigment pepper fruits under urea reduction. (A) Violin plot; (B) FPKM-based principal component analysis.
Agriculture 14 00521 g003
Agriculture 14 00521 g004
Figure 4. Analysis of differentially expressed genes (DEGs) in RNA-Seq of pigment pepper fruits under urea reduction. (A) Number of upregulated and downregulated DEGs; (B) Venn diagrams representing the intersection of DEGs between different nitrogen treatments.
Figure 4. Analysis of differentially expressed genes (DEGs) in RNA-Seq of pigment pepper fruits under urea reduction. (A) Number of upregulated and downregulated DEGs; (B) Venn diagrams representing the intersection of DEGs between different nitrogen treatments.
Agriculture 14 00521 g004
Agriculture 14 00521 g005
Figure 5. Enrichment analysis of DEGs metabolic pathways of pigment pepper fruits under urea reduction. (A) GO enrichment analysis; (B) KEGG enrichment analysis. The numbers in the figure represent the −log(p-value). Different colors represent the abundance of enriched genes, with orange indicating a higher number of genes, green indicating a lower number of genes, and white representing an intermediate level between the two.
Figure 5. Enrichment analysis of DEGs metabolic pathways of pigment pepper fruits under urea reduction. (A) GO enrichment analysis; (B) KEGG enrichment analysis. The numbers in the figure represent the −log(p-value). Different colors represent the abundance of enriched genes, with orange indicating a higher number of genes, green indicating a lower number of genes, and white representing an intermediate level between the two.
Agriculture 14 00521 g005
Agriculture 14 00521 g006
Figure 6. Correlation analysis between RNA-seq and qRT−PCR for carotenoid genes in pepper fruits under urea reduction.
Figure 6. Correlation analysis between RNA-seq and qRT−PCR for carotenoid genes in pepper fruits under urea reduction.
Agriculture 14 00521 g006
Agriculture 14 00521 g007
Figure 7. Analysis of differential transcription factors in pigment pepper fruits under urea reduction. (A) Top 10 transcription factor families; (B,C) heatmap of Pearson correlation analysis between NAC family and carotenoid genes of N1 (B) and N4 (C).
Figure 7. Analysis of differential transcription factors in pigment pepper fruits under urea reduction. (A) Top 10 transcription factor families; (B,C) heatmap of Pearson correlation analysis between NAC family and carotenoid genes of N1 (B) and N4 (C).
Agriculture 14 00521 g007
Agriculture 14 00521 g008
Figure 8. Silencing of the CA12g04950 gene in pigment pepper fruits using VIGS. (A) Fruits’ colors after silencing; (B) electropherogram of silenced CA12g04950 gene (left: DL2000 marker; middle and right: purified product); (C) total chlorophyll content of silenced fruits; (D) total carotenoid content of silenced fruits; (E) expression level of CA12g04950 in silenced fruits; * indicates significant differences between different treatments (* p ≤ 0.05, ** p ≤ 0.01, *** p ≤ 0.001, **** p ≤ 0.0001), ns indicates no significant differences between different treatments.
Figure 8. Silencing of the CA12g04950 gene in pigment pepper fruits using VIGS. (A) Fruits’ colors after silencing; (B) electropherogram of silenced CA12g04950 gene (left: DL2000 marker; middle and right: purified product); (C) total chlorophyll content of silenced fruits; (D) total carotenoid content of silenced fruits; (E) expression level of CA12g04950 in silenced fruits; * indicates significant differences between different treatments (* p ≤ 0.05, ** p ≤ 0.01, *** p ≤ 0.001, **** p ≤ 0.0001), ns indicates no significant differences between different treatments.
Agriculture 14 00521 g008
Agriculture 14 00521 g009
Figure 9. (A–D) Carotenoid Component Contents and (E–G) Metabolic Enzyme Activities of Pigment Pepper Fruits after Silencing the CA12g04950 Gene. (** p ≤ 0.01, *** p ≤ 0.001).
Figure 9. (A–D) Carotenoid Component Contents and (E–G) Metabolic Enzyme Activities of Pigment Pepper Fruits after Silencing the CA12g04950 Gene. (** p ≤ 0.01, *** p ≤ 0.001).
Agriculture 14 00521 g009
Agriculture 14 00521 g010
Figure 10. (A–F) Expression level of carotenoids’ genes of pigment pepper fruits after silencing the CA12g0495 (* p ≤ 0.05, ** p ≤ 0.01, *** p ≤ 0.001), ns indicates no significant differences between different treatments.
Figure 10. (A–F) Expression level of carotenoids’ genes of pigment pepper fruits after silencing the CA12g0495 (* p ≤ 0.05, ** p ≤ 0.01, *** p ≤ 0.001), ns indicates no significant differences between different treatments.
Agriculture 14 00521 g010
Agriculture 14 00521 g011
Figure 11. Working model of CA12g049501 in the synthesis of carotenoids in pigment pepper fruits under different nitrogen treatments and the expression of related genes. Red arrows represent up−regulation, green arrows represent down−regulation.
Figure 11. Working model of CA12g049501 in the synthesis of carotenoids in pigment pepper fruits under different nitrogen treatments and the expression of related genes. Red arrows represent up−regulation, green arrows represent down−regulation.
Agriculture 14 00521 g011
Table 1. Fertilizer dosage was applied on the pigment pepper.
Table 1. Fertilizer dosage was applied on the pigment pepper.
TreatmentUrea (kg/hm2)KH2PO4
(kg/hm2)		Calcium, Magnesium, and Sulfur Nutrients (kg/hm2)Micronutrient Fertilizer on Leaf (%)
N175020022.50.5%
N2562.520022.50.5%
N337520022.50.5%
N4187.520022.50.5%
N0020022.50.5%
Table 2. qRT-PCR and VIGS primers.
Table 2. qRT-PCR and VIGS primers.
GeneForward PrimerReverse Primer
qRT-PCR Primers
CA01g13640
(Indole-3-acetaldehyde oxidase-like)		CAGATGGATTGGTTGTTTCATGTCCGCTGTTTAGCACTT
CA04g18710
(Abscisic acid 8′-hydroxylase)		AAGGGAAAGGAAAGAGAAGAAGCAAGAAACAACTCAACAACCC
CA11g15430
(Aldehyde oxidase 4)		ACGAACCTTACCAGCAAATCATGACAGAGTCCACCTCCA
CA03g25820 (BCH)CCGCCACTTCTCTTCTTCTCGCTTTGTTTTCCACTTTG
CA01g09070
(Abscisic acid 8′-hydroxylase CYP707A1)		CCAACTTTACTCACAAAACCCTTCCCCAACATTCTCTCTTT
CA12g20770 (Z-ISO)TTTACCATTTCCCACCTCTCTTCACCCACCAAGATTTCA
CA08g10750 (ZDS)CATCAGTTGTGAAGATTGGGTGAGCCAGCAAGAAAGAAA
CA10g15950 (P450)GAAGTTGACAGAGTTTTGGGAGACAGGTGGATGTGGATAAAG
CA12g11380 (VDE)CTAATCCCAGTGTTCTTGTCCCAGGTGTCCGTATTCTCCA
CA01g33030
(Abscisic acid 8′-hydroxylase CYP707A2)		TGTCATTGGTGTCATATTTGCTATTGCCTCTTGTTCCTCTGT
CA03g36860 (PDS)AGATGGTTGCTCGCAAAGGACGCGGAGAAGATCGGAATGA
CA02g10990 (ZEP)TGCCAAACAAGCCAGGAGAAACACCTCATCCGTCACCCTA
CA06g22860 (CCS)AGTGGCCTGTGAGTTGTGTTGAGCCACCATGTACCCAGAC
CA04g04080 (PSY)CAAAGGCAACAACGGAGAGCCGCTCAATTCGGTCACTCCT
CA05g00080 (Lcyb)ACGTGGAGCTCAAGGAGAGAACGAGCCACCATTCGTTCTT
CA02g23040 (GGPS)ATTGAAGCAGCACAGACGGAGAGTATTGCGCCGAGTACGA
CA03g25820 (CrtZ)ACGAGTCACACCACAAACCATTTCCAACTCTTCCAGCCCC
CA12g04950GCAATGGTGATGAAGCGCAAAACAAACGTCAGAGGTCGGG
VIGS Primers (the underlined part represent the enzyme-cleaved homologous arm)
TRV-CA12g04950AAGGTTACCGAATT
TGCCGATGAAGAATTGGACAG		CTCGGTACCGGATC
CTCTTGCAAACAAACGTCAGAG
Table 3. RNA sequence qualities of pigment pepper fruits at S2 stage under nitrogen reduction.
Table 3. RNA sequence qualities of pigment pepper fruits at S2 stage under nitrogen reduction.
SampleRaw Read NumberClean
Reads		Total Mapped% of
Clean Reads		Raw Q30 Rate
N1_139,442,00036,904,79231,243,116 (84.66%)93.5793.61
N1_237,259,70234,983,96629,631,058 (84.70%)93.8993.54
N1_345,200,93042,473,75036,031,689 (84.83%)93.9793.9
N3_138,890,42236,556,26631,022,737 (84.86%)94.0093.32
N3_243,611,22640,987,50834,803,261 (84.91%)93.9893.82
N3_348,362,08045,251,53238,094,608 (84.18%)93.5794.05
N4_146,759,74043,958,72437,229,386 (84.69%)94.0193.7
N4_248,074,91844,964,77837,776,074 (84.01%)93.5393.6
N4_346,333,09843,253,84236,702,783 (84.85%)93.3593.16
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shen, L.; Zhang, C.; Xia, Y.; Yang, S.; Chang, T.; Ullah, S.; Ji, X. Transcript Analysis Reveals Positive Regulation of CA12g04950 on Carotenoids of Pigment Pepper Fruit under Nitrogen Reduction. Agriculture 2024, 14, 521. https://doi.org/10.3390/agriculture14040521

AMA Style

Shen L, Zhang C, Xia Y, Yang S, Chang T, Ullah S, Ji X. Transcript Analysis Reveals Positive Regulation of CA12g04950 on Carotenoids of Pigment Pepper Fruit under Nitrogen Reduction. Agriculture. 2024; 14(4):521. https://doi.org/10.3390/agriculture14040521

Chicago/Turabian Style

Shen, Lingfeng, Chenfei Zhang, Yahui Xia, Shasha Yang, Tian Chang, Saleem Ullah, and Xuehua Ji. 2024. "Transcript Analysis Reveals Positive Regulation of CA12g04950 on Carotenoids of Pigment Pepper Fruit under Nitrogen Reduction" Agriculture 14, no. 4: 521. https://doi.org/10.3390/agriculture14040521

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop