Dionis Borràs1†Lorenzo Barchi1†Karina Schulz2†Andrea Moglia1†
Alberto Acquadro1†Iman Kamranfar3†
Salma Balazadeh2,4†
Sergio Lanteri1*†- 1Department of Agricultural, Forest and Food Sciences, Plant Genetics and Breeding, University of Torino, Turin, Italy
- 2Max Planck Institute of Molecular Plant Physiology, Potsdam, Germany
- 3Department Molecular Biology, Institute of Biochemistry and Biology, University of Potsdam, Potsdam, Germany
- 4Plant Sciences and Natural Products, Institute of Biology Leiden (IBL), Leiden University, Leiden, Netherlands
Capsicum annuum L. is one of the most cultivated Solanaceae species, and in the open field, water limitation leading to drought stress affects its fruit quality, fruit setting, fruit size and ultimately yield. We identified stage-specific and a common core set of differentially expressed genes, following RNA-seq transcriptome analyses of a breeding line subjected to acute drought stress followed by recovery (rewatering), at three stages of plant development. Among them, two NAC transcription factor (TF) genes, i.e., CaNAC072 and CaNAC104, were always upregulated after drought stress and downregulated after recovery. The two TF proteins were observed to be localized in the nucleus following their transient expression in Nicotiana benthamiana leaves. The expression of the two NACs was also induced by NaCl, polyethylene glycol (PEG) and abscisic acid (ABA) treatments, suggesting that CaNAC072 is an early, while CaNAC104 is a late abiotic stress-responsive gene. Virus-induced gene silencing (VIGS) of CaNAC104 did not affect the pepper plantlet’s tolerance to drought stress, while VIGS of CaNAC072 increased drought tolerance. Heterologous expression of CaNAC072 in Arabidopsis thaliana as well as in plants mutated for its homolog ANAC072 did not increase drought stress tolerance. This highlights a different role of the two NAC homologs in the two species. Here, we discuss the complex role of NACs as transcriptional switches in the response to drought stress in bell pepper.
Introduction
Drought stress is one of the key limiting factors affecting plant growth, development and survival, with a substantial impact on crop yield (Sinclair, 2011; Boyer et al., 2013). The lack of water, due to rising temperatures and changes in precipitation patterns, represents a serious problem in many parts of the world, and it is expected to become a cause of food shortage and malnutrition with the increase in population and food demand (Newton et al., 2011; Fitton et al., 2019). Plants have evolved various molecular mechanisms to adjust their growth to limited water availability, and the elucidation of such mechanisms is essential for implementing breeding strategies aimed at developing crop varieties more resilient and capable to deal with water shortage while maintaining yield (Osakabe et al., 2014; Iovieno et al., 2016). Bell pepper (Capsicum annuum L.), also known as sweet pepper, is one of the most widely cultivated solanaceous vegetables (FAOSTAT, 2018), and particularly at flowering, water deficit causes flower dropping which reduces the setting and size of the fruits, and alters their biochemical composition (Borras et al., 2020).
Comparative transcriptomic analyses have proven to be a valuable approach for identifying key genes controlling the response to low water availability (Seki et al., 2002) in plants such as Arabidopsis (Zhang et al., 2019), rice (Chung et al., 2018), maize (Kakumanu et al., 2012), potato (Sprenger et al., 2018), and tomato (Lee et al., 2018). However, in C. annuum, studies on transcriptional changes occurring in response to water shortage are to date rather limited (Zinselmeier et al., 2002; Ji et al., 2010; Sharoni et al., 2012; Lee and Choi, 2013).
Previous studies in plants have highlighted that transcription factors (TFs) play a crucial role as molecular regulators to activate or inhibit the expression of stress-related genes and promote adaptation to stressful environments. Important abiotic stress-responsive TF families are AP2/ERF, DREB, bZIP, AREB/ABF, MYB, bHLH, WRKY and NAC (Nakashima et al., 2012; Singh and Laxmi, 2015; Gamboa-Tuz et al., 2018; Mittal et al., 2018; Mohanta et al., 2020). With respect to NACs, several members of this family have been identified and functionally characterized in crops, among which is the solanaceous species tomato (Yang et al., 2011; Han et al., 2012; Ma et al., 2013; Thirumalaikumar et al., 2018; Devkar et al., 2020). The typical NAC protein is characterized by a highly conserved DNA-binding NAC domain at its N-terminus and a divergent transcription regulatory region (TRR) at its C-terminus (Mao et al., 2007).
Whole-genome sequences of both peppers and closely related species have recently been made available (Kim et al., 2014; Qin et al., 2014; Kang et al., 2016; Ahn et al., 2018; Acquadro et al., 2020), and a comprehensive analysis of the NAC gene family has been performed. A whole of 104 CaNAC genes, of which 24 are exclusive to the Solanaceae family, were identified and found to be mostly located on chromosomes 1, 2, 3, and 6 (Diao et al., 2018). However, few NAC genes have been functionally characterized in C. annuum with respect to their function in stress responses. Recently, Zhang et al. (2020) isolated and functionally characterized CaNAC035, a positive regulator of abiotic stress tolerance acting through multiple signaling pathways. The involvement in the response to abiotic stress tolerance has also been demonstrated for CaNAC064 (Hou et al., 2020) and CaNAC2 (Guo et al., 2015).
In this study, we employed RNA-sequencing to identify differentially expressed genes (DEGs) in a C. annuum inbred line subjected to acute drought stress, and recovery from it, at different developmental stages. Common core sets and plant developmental stage-specific DEGs were identified. Among them, two members of the NAC family (namely CaNAC072 and CaNAC104) were significantly up-regulated transcriptionally following acute drought stress, and down-regulated after recovery across all developmental stages, and were thus selected for further functional characterization. Transient expression in Nicotiana benthamiana leaves revealed a nuclear localization of CaNAC072 and CaNAC104; however, a different time of transcriptional activation of the two pepper NAC genes was detected following NaCl, PEG or ABA treatments.
Both NACs were also functionally characterized by Virus-Induced Gene Silencing (VIGS) in pepper plantlets, and CaNAC072 was further assessed in three Arabidopsis thaliana genetic backgrounds.
Material and Methods
Drought Stress Treatments, RNA Extraction and Sequencing
Seeds of a yellow-fruited breeding line of Capsicum annuum (Cuneo CCu07), selected at the DISAFA (University of Turin), were placed on filter paper soaked with sterilized tap water. Petri dishes were placed in a growth chamber at 350 μmol photons m−2 s−1, 50% relative humidity, 26–24°C and a 16−8 h light-dark regime. After germination, seedlings were transferred to individual plastic pots filled with peat and grown for about 15 days. Before transfer to the greenhouse, the plantlets were transplanted in 5 L pots filled with a mixture of peat and perlite at a ratio of 1:1 by volume and grown following standard horticultural practices.
Groups of 15 plants were irrigated up to the soil water retention, and grown up to: 1) the production of five true leaves (stage 1); 2) the production of the third flower (stage 2); 3) the setting of the first fruit (stage 3). Ten plants at each developmental stage were subjected to drought stress by stopping irrigation, while irrigation was maintained on five other plants (Controls). The effect of drought stress was assessed by evaluating stomatal conductance (gs) through a Leaf Porometer Model SC-1 (Meter Group) on five fully expanded and randomly chosen leaves of each plant. Each measurement was replicated three times in each plant. After drought stress, a group of five plants was recovered following re-watering up to the soil water retention, until they reached the same gs as control plants.
At each stage of plant development, RNA was extracted using the NucleoSpin RNA kit (Macherey-Nagel) from three leaves of three control plants as well as three plants subjected to drought stress, and three plants subjected to drought stress followed by recovery, chosen randomly. The RNAs were treated with RNase-free DNase I at room temperature to remove contaminant DNA.
Prior to library construction, an Agilent 2,100 Bioanalyzer (Agilent, CA, United States) was used to confirm the quality and quantity of RNA such that the 28S/18 S rRNA ratio was >1.5 and the RNA integrity number >7. For each sample, mRNA was fragmented randomly by adding fragmentation buffer, then the cDNA was synthesized by using mRNA template and random hexamer primers, after which a custom second-strand synthesis buffer (Illumina), dNTPs, RNase H and DNA polymerase I were added to initiate the second-strand synthesis. After terminal-end repair and ligation of sequencing adaptors, the double-stranded cDNA library was completed through size selection and PCR enrichment. Libraries were sequenced on an Illumina NovaSeq 6,000 platform at Novogene, and 150 bp paired-end reads were generated. Raw reads are publicly available in the NCBI Sequence Read Archive under Bioproject PRJNA668245.
Identification of Differentially Expressed Genes and Enrichment Analysis
Raw reads obtained from the Illumina NovaSeq 6,000 platform were processed for standard quality controls (QC) using FastQC (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Then they were cleaned with Scythe (v0.994, https://github.com/vsbuffalo/scythe) for removing contaminant residual adapters, and Sickle (v1.33, https://github.com/najoshi/sickle) to remove reads with poor quality ends (Q < 30). The cleaned reads were mapped to the reference CM334 pepper genome (Kim et al., 2017) with Hisat2 (Kim et al., 2015). Data quality analysis was conducted by visualizing results of principal component analysis (PCA). The counts of mapped reads per gene per sample were obtained with Stringtie (Pertea et al., 2015). Then, the DESeq2 R package was applied to identify differentially expressed genes (DEGs) (Love et al., 2014) at a False Discovery Rate (FDR) threshold of 0.05. The genes with an absolute log2 fold change ≥1 were considered as DEGs.
Functional enrichment analyses of Gene Ontology (GO) annotation was performed to identify which DEGs were significantly enriched in GO terms. GO enrichment analysis was conducted using AGRIGOv2, using the available annotation of the previously published sequence of the inbred line CCu07 (Acquadro et al., 2020).
NAC Gene Prediction and Transcription Factor Identification
The Hidden Markov Model (HMM) profile of the NAC domain (PF02365), downloaded from the Protein families (Pfam) database (El-Gebali et al., 2019), was used for identifying NAC genes in the pepper CM334 genome (Kim et al., 2017) using HMMER (V 3.2.1). Positive matches were manually checked to confirm the presence of NAC domains, using InterProScan (Mitchell et al., 2019). Transcription factors in the reference genome CM334 were identified following a BlastP analysis against the Plant Transcription Factor Database (http://planttfdb.gao-lab.org/). To identify Arabidopsis orthologs of two NACs, which were always upregulated after drought stress and downregulated after recovery (i.e., CA.PGAv.1.6.scaffold213.1 and CA.PGAv.1.6.scaffold1156.46), a multiple alignment with the entire Arabidopsis NAC family was carried out by applying Clustal Omega (1.2.4) (Madeira et al., 2019). A further Pairwise Sequence Alignment with Emboss Needle (Madeira et al., 2019), for both pepper NACs with the closest Arabidopsis gene was used to retrieve identity and similarity percentage.
Subcellular Localization
Subcellular localization of CaNAC072 and CaNAC104 was predicted by identifying Nuclear Localization Signal (NLS) sequences using the Localizer software (http://localizer.csiro.au/).
The coding sequences (without stop codon) of CaNAC072 and CaNAC104 were amplified from CCu07 cDNA by using primers CaNAC072-F, CaNAC072-R, CaNAC104-F and CaNAC104-R (Supplementary Table S1). PCR was performed with Phusion High-Fidelity DNA polymerase (Thermo Fisher Scientific) and the amplicon sequences were validated (Eurofins MWG Operon). Amplicons were cloned into pDONR207 vector using BP clonase (Invitrogen). The sequence-confirmed entry vector was recombined into pK7FWG2, containing kanamycin as selection marker (Karimi et al., 2002), using LR reaction mix II (Invitrogen). The recombined plasmids (35S:CaNAC072-GFP and 35S:CaNAC104-GFP) were transformed into Agrobacterium tumefaciens (strain GV3101) and then used for infiltration of Nicotiana benthamiana leaves (Senthil-Kumar and Mysore, 2014). GFP signals were analyzed 2 days after leaf infiltration using a Leica DM6000B/SP5 confocal laser-scanning microscope (Leica Microsystems, Wetzlar, Germany). The laser excitation wavelength was 488 nm, and the detection was set at 507 nm.
Further confirmation of nuclear signal was assessed with nuclear-specific DAPI (4,6-diamidino-2-phenylindole) staining. Transformed leaves were incubated for 30 min in darkness with a solution containing 10 mg L−1 DAPI. Afterwards, nuclear-specific signal was assessed by confocal laser-scanning microscopy as described by Sampathkumar and Wightman (2015).
Stress and Hormonal Treatments
CCu07 seeds were germinated on full-strength Murashige-Skoog (MS) medium containing 2% (w/v) sucrose, and grown under a 24–26°C and a 16–8 h light-dark regime. Ten-day-old seedlings were transferred to MS liquid media containing 120 mM NaCl (Sigma-Aldrich), 10% (w/v) PEG 6000 (Sigma-Aldrich), or 100 μM ABA (Sigma-Aldrich), and treated for 2, 4, 6, 10 and 24 h. Stressors and hormones were omitted from control seedlings.
Virus-Induced Gene Silencing
For the virus-induced gene silencing (VIGS) experiment, CCu07 seeds were sown directly into soil containing a mixture of peat and vermiculite (2:1, v/v), and plants grown in a phytotron at 450 μmol photons m−2s−1, 50% relative humidity, 21°C day/night under a 14 h day/10 h night regime.
VIGS was performed with vectors pTRV1 (pYL192) (Liu et al., 2002) and pTRV2 (pYL170) following the protocol described by Senthil-Kumar and Mysore, 2014. Briefly, VIGS constructs were designed with the SGN VIGS Tool (Fernandez-Pozo et al., 2015). The selected gene fragments (Supplementary Table S2) were PCR-amplified using primers VIGS_CaNAC072-F, VIGS_CaNAC072-R, and VIGS_CaNAC104-F, VIGS_CaNAC104-R (Supplementary Table S1) from CCu07 cDNA and cloned into pDONR207 using BP clonase (Invitrogen). The sequence-confirmed entry vector was recombined into pK7FWG2 using LR reaction mix II (Life Technologies). The plasmids pTRV1, pTRV2:CaNAC072, pTRV2:CaNAC104 as well as pTRV2:AtPDS (as a positive control for efficient VIGS) and empty vector pTRV2:00 were used to transform A. tumefaciens (strain GV3101), and the latter used for infecting cotyledons of 10 days-old pepper seedlings. After 12 days of infection, drought stress was performed for 7 and 9 days. Recovery treatments were performed by re-watering plants for 7 days. For assessing plant survival rate, plants were subjected to 9 days severe drought stress followed by recovery through re-watering. The percentage of surviving plants was determined after 7 days of recovery. The drought and recovery phenotypic analyses were assessed in 15 biological replicates per treatment.
Identification of Pepper Homologous Genes in Arabidopsis and Plant Transformation
Pepper orthologs of Arabidopsis drought marker genes [i.e., ANAC019 (AT1G52890), ANAC055 (AT3G15500), DREB2A (AT5G05410), DELLA/GAI (AT1G14920), and SnRK2 (AT3G50500)], were identified using PLAZA 4.5 (https://bioinformatics.psb.ugent.be/plaza/versions/plaza_v4_5_dicots/) and annotated using the Sol Genomics Database (https://solgenomics.net/).
The 35S:CaNAC072-GFP and 35S:CaNAC104-GFP constructs were inserted into A. tumefaciens (strain GV3101) and transformed into A. thaliana genotypes Columbia-0 (Col-0), anac072 mutant (rd26-2, SALK_083,756) and anac019/anac055/anac072 triple mutant backgrounds using the floral dip method (Clough and Bent, 1998). NAC T-DNA insertion lines (anac019: SALK_096,295, anac055: SALK_014331, and anac072: SALK_083756) were obtained from NASC (http://arabidopsis.info). Homozygous knockout plants were identified by PCR-based genotyping, qRT-PCR, and end-point PCR. The triple NAC mutant was generated by crossing. Plants were grown at 22°C under long-day (LD) condition, i.e., 16 h light/8 h dark. Seeds of T2 plants were harvested.
From T1 transgenic plants showing a 3:1 segregation ratio when grown on MS plates supplemented with the same antibiotic. Ten-day-old seedlings were transferred to soil [potting and quartz sand; 2:1 (v/v)]. After 14 days, their drought tolerance was assessed by withholding water for 20 days.
Relative Water Content and Water Loss Assay
Relative water content (RWC) was measured as described by Thirumalaikumar et al., 2018. Briefly, second true leaves were harvested and their fresh weights (FW) were determined. Leaves were then immersed in distilled water and incubated overnight at 4°C to obtain the saturated weight (SW). Next, leaves were dried at 70°C for 48 h to measure the dry weight (DW). RWC was calculated as % by the formula (FW-DW)/(SW-DW) x 100. RWC was measured in 15 biological replicates per genotype per treatment.
For water loss assays, second true leaves were detached and immediately weighted (FW). Leaves were then placed on filter paper in a growth chamber at 22°C and the weight measured after 1, 2, 3, 4, 5 and 6 h. Water loss was measured in 15 biological replicates each.
Expression Profiling
Total RNA extraction, cDNA synthesis, and qRT-PCR were performed as described by Balazadeh et al., 2008. The qRT-PCR primers (Supplementary Table S3) were designed using QuantPrime (Arvidsson et al., 2008). ABI-PRISM 7900 HT sequence detection system (Applied Biosystems, Darmstadt, Germany) and SYBR Green (Applied Biosystems, Darmstadt, Germany) were used for qRT-PCR and detection of amplified product.
The transcription pattern of CaNAC072 and CaNAC104 was assessed in ten-day-old wild-type pepper seedlings subjected to salinity induced by 120 mM NaCl, drought induced by 10% (w/v) PEG 6000, as well as after treatments with 100 μM ABA over 24 h. Three biological replicates, each a pool of five seedlings, were used for gene expression profiling. The homologue of the Arabidopsis RD29A gene, known to be induced by various abiotic stressors (Msanne et al., 2011), i.e., CaRD29A (CA.PGAv.1.6.scaffold490.32), was used as a positive control. Data were quantified using the 2− ΔΔCt method based on Ct values of CaNAC072 and CaNAC104, and pepper ACTIN, (CA.PGAv.1.6.scaffold786.11) and GADPH (TC.CA12g18170) as housekeeping genes.
In the VIGS experiment, the expression pattern of CaNAC072 (CA.PGAv.1.6.scaffold213.1), CaNAC104 (CA.PGAv.1.6.scaffold1156.46), CaNAC019 (CA.PGAv.1.6.scaffold643.12), CaNAC055 (TC.CA07g18020), CaDREB2A (TC.CA05g15500), CaDELLA (TC.CA12g02780), and CaSnRK2 (TC.CA05g02700) was assessed. RNA was extracted from pepper plants infiltrated with TRV2:CaNAC072, TRV2:CaNAC104 and TRV2:00 and subjected to 7 days of water deprivation. Three biological replicates were used for gene expression analyses.
Results
Effects of Drought Stress and Recovery on Stomatal Conductance
We studied the effects of drought stress and recovery (re-watering) on a C. annuum inbred line at different plant developmental stages. To this end, plants at the status of five true leaves (stage 1), the third flower (stage 2), and at the setting of the first fruit (stage 3) were exposed to drought stress, by withholding water, and then allowed to recover by re-watering. In the control condition, plants were regularly irrigated. Drought stress was imposed by interruption of irrigation. For recovery experiments, drought-stressed plants were re-watered up to soil water retention (Figure 1). Leaves were sampled after each treatment and used for RNA-seq-based transcriptome profiling.
FIGURE 1. Experimental design of the treatments for control, drought and recovery plants. Grey micro test tubes indicate when the samples were collected. The colour of the bars represents the irrigation status of the plants (light blue = normal irrigation; orange = no irrigation; green = irrigation, 24 h recovery).
The effect of water deprivation was assessed through a repeated evaluation of stomatal conductance (gs). Control plants showed values of 400 ± 60 mmol m−2s−1 on average, while drought-stressed plants displayed a progressive decrease of gs until about 50 mmol m−2s−1 and an evident foliar withering. Afterwards, plants were re-watered, and after 24 h reached a gs value analogous to control plants, and they recovered their leaf turgor.
RNA-Seq Analyses
RNA-sequencing profiles were obtained from pepper leaves of control, drought-stressed, and drought stressed and then recovered plants at the three developmental stages, as described above. About 33 million paired reads were obtained from each sample. After removing low-quality sequences, adapters and possible contaminations, around 30 million reads (PE and singlets) per sample were obtained (Supplementary Table S4). Clean reads were aligned against the CM334 pepper genome (Kim et al., 2017), resulting in approximately 92% mapping reads on average.
Following quality control by PCA analysis, the three replicates for each condition (control, drought and recovery) clustered together. Only one drought stress replicate at fruit stage was found to be separated from the other two replicates. In addition, a clear separation was found between drought stress replicates and control/recovery samples (Supplementary Figure S1).
Using a false discovery rate threshold of 0.05 and an absolute log2 fold change ≥1 (adjusted p ≤ 0.01) as the significance cutoff, the expression of 3,680, 3,976 and 1,282 DEGs was found to be significantly regulated by comparing control to drought-stressed plants at stage 1, 2 and 3, respectively. By comparing drought-stressed and recovered plants a whole of 4,155, 2010, 1,328 DEGs were found to be significantly regulated at stage 1, 2 and 3, respectively. At last, when recovered plants were compared with control plants 1,285, 271 and 56 genes where differentially expressed at stages 1, 2 and 3, respectively (Supplementary Table S5).
Annotation and GO Enrichment of Core Sets of Commonly up- or Downregulated Genes
All DEG sets were compared in order to identify groups of commonly regulated genes. At all stages of plant development, a whole of 197 and 123 genes were found to be commonly up- and downregulated after drought, respectively (Figures 2A,B). On the other hand, 90 and 187 genes were found to be commonly up- and downregulated following plant recovery, respectively (Figures 2C,D).
FIGURE 2. Venn diagram showing intersections of up- or downregulated DEGs, and enriched GO terms, among the DEGs in response to drought and after recovery. (A) Genes upregulated after water stress, compared to control plants. (B) Genes downregulated after water stress. (C) Genes upregulated after recovery. (D) Genes downregulated after recovery. (E) Enriched GO terms are shown for the three stages of plant development (i.e., St.1 = five true leaves, St.2 = production of third flower, and St.3 = setting of first fruit). Red color represents upregulated GO terms whereas blue color represents downregulated GO terms. Colors on the FDR columns represent different significance levels (yellow: FDR <0.05; orange: FDR <0.01; red: FDR <0.001.
In order to compare the function enrichment of commonly up- or downregulated genes at all stages of plant development, a gene ontology (GO) enrichment analysis was performed. The common over-represented GO terms at all stages are summarized in Figure 2E. Nineteen significant GO terms were identified among all developmental stages. When comparing drought vs control treatments, three enriched GO terms were identified among upregulated DEGs (i.e., GO:0009628, “response to abiotic stimulus”; GO:0009415, “response to water”, and GO:0016491, “oxidoreductase activity”), while seven GO terms were identified among the downregulated DEGs (including e.g. GO:0046982, “protein heterodimerization activity”; GO:0032993, “protein-DNA complex”). On the other hand, when comparing recovery vs drought treatments, two enriched GO terms were identified on upregulated DEGs, i.e., GO:0004553, “hydrolase activity, hydrolyzing O-glycosyl compounds”, and GO:0016798, “hydrolase activity, acting on glycosyl bonds”), while seven were identified among the downregulated DEGs (including e.g., GO:0009415, “response to water”; GO:0009628, “response to abiotic stimulus”). No common enriched GO terms were found by comparing recovery vs. control. All GO terms, including stage-specific enriched GO terms, can be inspected in the Supplementary Tables S6−S17.
Among the 197 commonly up-regulated genes after drought stress, we focused our attention on the 30 most expressed genes at all stages of plant development, and 19 were found to be shared (Figure 3A). Among them are two Late Embryogenesis Abundant protein genes, LEA46 and LEA-D29 (CA.PGAv.1.6.scaffold1161.31 and CA.PGAv.1.6.scaffold407.111), several stress-induced genes, such as ABA metabolism-related ASR1, TAS14 and ATHB-7 (CA.PGAv.1.6.scaffold264.2, CA.PGAv.1.6.scaffold358.34, and CA.PGAv.1.6.scaffold567.90) as well as some plant cell wall-related genes, i.e. EXLB1, AAA1 KATANIN (CA.PGAv.1.6.scaffold572.3 and CA.PGAv.1.6.scaffold532.70). Likewise, among the 123 commonly downregulated genes after drought stress, we focused our attention on the 30 most downregulated genes. Six of them were found to be shared and included one gene involved in fatty acid biosynthesis, i.e., 3-KETOACYL-COA SYNTHASE 3 (KCS3; CA.PGAv.1.6.scaffold1289.4), and a gene encoding a serine protease (i.e., SBT1.7; CA.PGAv.1.6.scaffold1394.16) (Figure 3B).
FIGURE 3. Heatmap showing the log2 fold change differences in the expression profiles of the top 30 common DEGs at the three stages of plant development. Genes upregulated after drought stress or recovery (Panel A and D) and genes downregulated after drought stress or recovery (Panel B and C). Data represent the means of three biological replicates. Asterisks represent statistically significant differences (p values) between control plants and stressed (droughted) plants, and between stressed (droughted) plants and recovery plants. The p values were adjusted using the Benjamini and Hochberg method. Corrected p value of 0.005 and |log2 (Fold Change)| of 1 were set as the threshold for significantly different expression. White fields indicate the absence of significant differences between treatments. NA stands for “Not applicable”.
After recovery, 12 upregulated genes were found to be in common among the top 30 most upregulated ones. Among them, MYO-INOSITOL OXYGENASE 1 (MIOX1; CA.PGAv.1.6.scaffold1241.1) is involved in providing nucleotide sugars for cell wall polymers, and 3-KETOACYL-COA SYNTHASE 3 (KCS3). Lastly, 18 genes were downregulated after recovery, 13 of which were observed to be upregulated after stress. Among them are BURP DOMAIN CONTAINING PROTEIN 3 and homeobox-leucine zipper protein gene ATHB-7 (CA.PGAv.1.6.scaffold186.49 and CA.PGAv.1.6.scaffold1045.47) (Figures 3C,D).
Differentially Expressed TFs After Drought Stress and Recovery
After drought, the RNA-seq expression profiling revealed 295, 370 and 108 differentially expressed (DE) TFs at stage 1, 2 and 3 of plant development, respectively. Among the ones at stage 1, 162 were up- and 133 downregulated, while at stage 2, 188 were up- and 183 downregulated, and at stage 3, 80 were up- and 28 downregulated (Supplementary Table S18). On the other hand, after the recovery phase a whole of 325, 176 and 106 DE TFs were identified at stage 1, 2 and 3 of plant development, respectively. Of those, at stage 1, 138 were up- and 187 downregulated; at stage 2, 84 were up- and 92 downregulated, while at stage 3, 33 were up- and 73 downregulated (Supplementary Table S18).
Overall, after drought and recovery significant changes occurred in the transcript levels of the five following TF families: MYB, bHLH, ERF, NAC and HSF. Stage-specific TF families were also identified (Supplementary Figure S2). A core set of 31 TF genes were found to be always up regulated after drought and downregulated after recovery at all stages of plant development, among which where MYB, bHLH, CO-like and NF-YA, as well as two NAC TFs (Supplementary Table S19).
Genome-wide Identification of NAC TFs and Their Expression After Drought and Recovery
We performed a comprehensive analysis to identify NAC genes in the pepper CM334 genome sequence (Kim et al., 2017) on the basis of the NAC domain (PF02365). A whole of 113 NAC TFs were identified (Supplementary Table S20), of which 13, 20 and 9 NACs were differentially expressed after drought stress at stage 1, 2 and 3, respectively (Figures 4A,B). On the other hand, 13, 10, and 8 NAC TFs were found differentially expressed after recovery (Figures 4C,D).
FIGURE 4. Heatmap showing the log2 fold change differences in the expression of the differentially expressed NAC TFs identified. Genes upregulated after drought stress or recovery (Panel A and D) and genes downregulated after drought stress or recovery (Panel B and C). Black arrows highlight the two NACs selected for functional characterization, being always upregulated after drought, and downregulated after recovery. Data represent the means of three biological replicates. Asterisks represent statistically significant differences (p values < 0.05*, Pvalues < 0.01**) between control and stressed plants, and between drought-stressed and recovered plants. The p values were adjusted using the Benjamini and Hochberg method. Corrected p value of 0.005 and |log2 (Fold Change)| of 1 were set as the threshold for significantly different expression. White fields indicate the absence of significant differences between treatments. NA stands for “Not applicable”.
On the basis of our RNA-seq experiment, two members of the NAC TF family, i.e. CA.PGAv.1.6.scaffold213.1 (referred to as CaNAC072 in the following) and CA.PGAv.1.6.scaffold1156.46 (CaNAC104) were found to be always upregulated after drought stress and downregulated after recovery, and thus were selected for functional characterization.
Subcellular Localization of CaNAC072 and CaNAC104
In-silico prediction analyses of the encoding sequences identified one NLS in CaNAC072 (RKNGSSKLDEWVLCRIYKK) and one in CaNAC104 (KKRK), suggesting a nuclear localization of both proteins. This was confirmed by fusing in-frame the coding sequences of CaNAC072 and CaNAC104, without stop codon, to the green fluorescent protein (GFP) coding sequence under the control of the CaMV 35S promoter. The subcellular localization of each protein was analyzed by confocal laser scanning microscopy after transient expression in N. benthamiana leaves. The analysis showed that the GFP signal was localized in the nucleus of leaves transfected with both CaNAC072- and CaNAC104-GFP fusion constructs (Figure 5). This observation is consistent with the role of CaNAC072 and CaNAC104 being transcription factors.
FIGURE 5. Nuclear localization of CaNAC072 and CaNAC104. Confocal microscopy analysis showing the nuclear localization of CaNAC072:GFP and CaNAC104:GFP upon transient expression in Nicotiana benthamiana leaves. Nuclear-specific staining using DAPI confirmed nuclear localization of both TFs. “Merged” indicates the combination of GFP and chlorophyll (Chl) autofluorescence signals. White arrows indicate the nuclear-specific GFP signal. Scale bar, 50 μm.
Effects of Abiotic Stresses on CaNAC072 and CaNAC104 Expression
NAC TFs play key roles in plant stress responses. CaNAC072 is the closest homologue of Arabidopsis ANAC072, with an identity of 56.3% and similarity of 67.6% (E value: 5e-132). Lower similarity percentages were detected for the ANAC072 homologues: ANAC019 (identity 53.0%, similarity 64.4%), and ANAC055 (identity 52.8%, similarity 63.1%).
On the other hand, CaNAC104 was found the closest homolog of Arabidopsis ANAC104, with an identity of 40.5% and similarity of 54.8% (E value: 1e-43), the latter being induced by several abiotic stresses and playing a major role in root xylem vessels (Kilian et al., 2007; Tang et al., 2018). Lower similarity percentages were also detected for the ANAC0104 homologues: ANAC02 (identity 5.7%, similarity 10.1%), and ANAC05 (identity 3.8%, similarity 6.6%).
To validate the response of CaNAC072 and CaNAC104 to abiotic stresses, we assessed the effects of salinity induced by NaCl, and drought induced by PEG, over 24 h in CCu07 pepper seedlings. Since both salinity and drought stimulate the accumulation of ABA (Marusig and Tombesi, 2020), we also validated gene expression induction upon ABA treatment.
Our treatments were effective in inducing the expression of the stress-related gene CaRD29A. The NaCl treatment induced its highest level of expression after 2 h, which progressively decreased in the following hours. Differently, the PEG treatment induced an increase of CaRD29A expression after 4 h and a further increase after 6 h, which was followed by a progressive decrease. At last, treatment with ABA caused an increase of CaRD29A expression after 2 h onwards following the hormonal treatment (Supplementary Figure S3).
As shown in Figure 6, CaNAC072 transcript levels significantly increased after 2 h of NaCl treatment and gradually decreased in the following hours. The PEG treatment induced the expression of CaNAC072 only after 4 and 6 h of treatment. Treatment with ABA resulted in a rapid rise (already after 2 h) in the expression of CaNAC072. Interestingly, this upregulation was maintained for up to 24 h. In contrast, the NaCl treatment induced the upregulation of CaNAC104 after 4 h and its transcript level progressively decreased in the following hours, while the PEG treatment did not induce its expression. On the other hand, ABA treatment induced a progressive increase of CaNAC104 expression from 4 h onwards (Figure 6). These results confirm that CaNAC072 and CaNAC104 are differentially induced by salinity and drought and suggest that they are ABA-dependent, with CaNAC072 being an early-responsive gene and CaNAC104 acting as a late-responsive gene.
FIGURE 6. Transcript levels of CaNAC072 and CaNAC104 after stress treatments. Gene expression was quantified by qRT-PCR after 0, 2, 4, 6, 10 and 24 h of treatments with NaCl, PEG and ABA. Ten-day-old seedlings were transferred to liquid MS media containing 120 mM NaCl, 10% (w/v) PEG 6000, or 100 μM ABA. Data represent the means ± SD (n = 3 biological replicates). Pepper actin and GADPH genes were used as the reference genes. Asterisks (*) indicate statistically significant differences between control condition (0 h of treatment) and stress treatment according to Student’s t-test; *, p < 0.05; **, p < 0.01; n.s. not significant.
Effect of Drought Stress on CaNAC072- and CaNAC104-Silenced Pepper Plants
To investigate the role of CaNAC072 and CaNAC104 for the response to drought, we reduced their expression through Virus-Induced Gene Silencing (VIGS) in Cuneo CCu07 pepper plantlets. The efficiency of VIGS was confirmed by generating TRV2:AtPDS plants, in which the silencing of the PHYTOENE DESATURASE (PDS) is easily displayed by photo-bleaching of leaves. As expected, plants infected with the empty vector TRV2:00 did not show leaf bleaching (Figure 7A).
FIGURE 7. Effect of drought stress on CaNAC072- and CaNAC104-silenced pepper plants. (A) Phenotype of a 25-day-old pepper plant after 12 days of agro-infiltration with TRV2:AtPDS (left) and with the empty vector TRV2:00 (right). (B). Expression of CaNAC072 and CaNAC104 in well-watered control plants (grey bar) and drought-stressed plants (black bar) after 7 days of water withholding (7-days-WW), determined by qRT-PCR. Data represent the means ± SD (n = 3 biological replicates). Values are expressed as the difference between an arbitrary value of 40 and dCt, so that high 40-dCt values indicate high gene expression levels. Pepper actin and GADPH genes were used as the reference genes. (C) Phenotype of TRV2:00 vis-à-vis with TRV2:CaNAC072 and TRV2:CaNAC104 plants in watering conditions (upper panel), and after 7 days of water withholding (lower panel). (D) Relative Water Content (RWC) of TRV2:00, TRV2:CaNAC072 and TRV2:CaNAC104 plants watered (grey bars) and after 7 days of no watering (black bars). Data represent the means ± SD (n = 15 biological replicates). (E) Trend of leaf fresh weights over 6 h (index of water loss) in detached leaves of TRV2:00, TRV2:CaNAC072 and TRV2:CaNAC104 plants. Data represent means ± SD of 15 biological replicates. (F) Heatmap showing the log2 fold change (FC) differences in the expression of five pepper genes (i.e., CaNAC019, CaNAC055, CaDREB2A, CaDELLA and CaSnRK2) which are orthologues of Arabidopsis drought-responsive genes, in TRV2:CaNAC072, TRV2:CaNAC104, and TRV2:00 plants watered and after 7 days of water withholding. Transcript levels were evaluated by qRT-PCR. Data represent the means of three biological replicates per assay. Asterisks indicate statistically significant differences between CaNAC072-non silenced and CaNAC072-silenced plants, and between CaNAC104-non silenced and CaNAC104-silenced plants according to Student’s t-test; *p < 0.05, **p < 0.01.
CaNAC072- and CaNAC104-silenced plantlets (TRV2:CaNAC072 and TRV2:CaNAC104, respectively), as well as the ones infected with the empty vector (TRV2:00), were subjected to drought stress (water withholding for 7 days), followed by recovery (rewatering for 7 days). No phenotype alterations were observed between non-infiltrated plants and TRV2:00-infected plants. Similarly, no phenotype alterations were detected between CaNAC072- and CaNAC104-silenced plants compared to TRV2:00 plants under well-watered conditions. After drought stress, a significant increase in the transcript levels of both CaNAC072 and CaNAC104 was detected in plantlets carrying the empty vector TRV2:00, but not in plantlets in which the two NAC TFs were silenced, further confirming the effectiveness of the VIGS (Figure 7B).
After 7 days of drought (water withholding), the TRV2:00 and TRV2:CaNAC104 plantlets showed severe leaf wilting and a marked decrease in their relative leaf water content (RWC), while TRV2:CaNAC072 plants did not, and their RWC was not significantly different from that of well-watered control plants (Figures 7C,D). In order to assess the rate of water loss in silenced plants, detached leaves of TRV2:00, TRV2:CaNAC072 and TRV2:CaNAC104 plants were analysed over 6 h and, as expected, CaNAC072-silenced plants showed a better water retention than TRV2:00 and CaNAC104-silenced plants (Figure 7E).
We also assessed the survival rate of TRV2:CaNAC072, TRV2:CaNAC104 and TRV2:00 plants after 9 days of water withholding, followed by 7 days of recovery (rewatering) (Supplementary Figure S4). Consistent with their higher RWC and not-withering phenotype, 15.2% of the CaNAC072-silenced plants survived, while none of the TRV2:00 and CaNAC104-silenced plants did. This result suggests that CaNAC072 acts as a negative regulator of the drought stress response.
In pepper infiltrated with TRV2:CaNAC072, TRV2:CaNAC104 and TRV2:00, we assessed the transcript profile of five pepper genes orthologous to drought stress-responsive genes in Arabidopsis, i.e. CaNAC019, CaNAC055, CaDREB2A, CaDELLA and CaSnRK2 (Sakuma et al., 2002; Ye et al., 2017). The qRT-PCR analyses revealed that, after 7 days of water withholding, the expression of all of them significantly increased in TRV2:CaNAC072 leaves, while only a bland activation of CaDELLA was observed in TRV2:CaNAC104 (Figure 7E). This might indicate a more rapid perception and response to drought stress in CaNAC072-silenced plants, which ultimately might explain their improved drought tolerance.
CaNAC072 Does Not Recover the Function of Its Ortholog in Arabidopsis
ANAC072 (also known as RD26) and its homologs ANAC019 and ANAC055 are drought-responsive genes promoting drought stress tolerance in Arabidopsis (Tran et al., 2004). Previous studies showed that the overexpression of either ANAC019, ANAC055 or ANAC072 confers drought tolerance, while the anac019/anac055/anac072 triple mutant shows enhanced drought sensitivity (Tran et al., 2004; Ye et al., 2017; Sukiran et al., 2019). In this study, we observed a drought tolerance increase in CaNAC072-silenced pepper plants, suggesting a different function in comparison to its homolog ANAC072 in Arabidopsis. Furthermore, we generated Arabidopsis plants overexpressing CaNAC072 fused to green fluorescent protein (GFP) (hereafter, CaNAC072:GFP-OX plants) in Col-0 as well as in anac072 single and anac019/anac055/anac072 triple mutants. Neither differences in the growth rate nor in the response to severe drought stress were detected between CaNAC072:GFP/Col-0, CaNAC072:GFP/rd26-2, and CaNAC072:GFP/triple mutant plants compared to Col-0, rd26-2 and the triple mutant (Supplementary Figure S5). Our results thus suggest that CaNAC072 from bell pepper plays a different function than its homolog ANAC072 in Arabidopsis.
Discussion
The response to drought in plants is a complex and coordinated process involving many biological mechanisms. Here, we subjected Capsicum annuum plants of a highly homozygous breeding line to acute water stress until leaf withering, followed by rewatering (recovery phase), at three developmental stages. Following sequencing of RNA extracted from leaves of drought-stressed, recovered and control plants, thousands of DEGs were identified, among which were core sets of commonly and plant development-specific up- and downregulated genes.
Abscisic acid (ABA) is known to act as an endogenous messenger in the regulation of the plant’s water status, and its action can induce stomatal closure, but also systemically signal for adjustment towards severe water shortage (Tuteja, 2007; Todaka et al., 2019). We detected ABA metabolism-related genes downregulated (e.g., ATHB-7) or activated (e.g., ASR1 and TAS14) upon drought stress (Godoy et al., 1994; Söderman et al., 1996; Dominguez and Carrari, 2015). Indeed, TAS14 was found to increase drought and salinity tolerance when overexpressed in transgenic tomato plants, and it was suggested as a promising candidate for enhancing drought tolerance in pepper as well (Muñoz-Mayor et al., 2012).
In plants exposed to drought stress, hydrophilic proteins such as Late Embryogenesis Abundant (LEA) proteins accumulate and cell walls reorganize in order to preserve cell turgor (Shinozaki and Yamaguchi-Shinozaki, 2007; Gong et al., 2010). We detected upregulation of two LEA protein genes: LEA46 and LEA-D29, whose overexpression in Arabidopsis confers tolerance to severe drought (Olvera-Carrillo et al., 2010), as well as of some cell wall-related genes, such as EXLB1 and AAA1 KATANIN, which were upregulated following drought stress and downregulated after recovery. Furthermore, genes such as MIOX1, involved in providing nucleotide sugars for cell wall polymers, and 3-KETOACYL-COA SYNTHASE 3 (KCS3), which is required for cuticular wax and root suberin biosynthesis, showed upregulation just after recovery. All in all, our transcriptomic data confirm that cell wall reorganization is a complex cellular process actively involved in the response of C. annuum to water stress and deserves to be better explored for enhancing drought tolerance in the future.
As expected, we identified significant changes in the expression of members of TF families, which are known to be related to stress responses. Among them are bHLH genes involved in the ABA response, HSF genes acting as regulators that prevent the accumulation of damaged proteins, as well as MYBs, ERFs and NACs (Trujillo et al., 2008; Gong et al., 2010); NACs represent one of the largest plant-specific transcription factor gene families (Nakashima et al., 2012; Franco-Zorrilla et al., 2014; Thirumalaikumar et al., 2018).
Previously in Capsicum annuum a whole of 104 CaNAC genes were identified by using 47 NAC proteins as query (Diao et al., 2018); here, we identified 113 sequences putatively assigned to the NAC family, presumably due to the different bioinformatics approach adopted. However, to date, only three NACs have been functionally characterized in C. annuum: CaNAC035, which acts through multiple signaling pathways and regulates the tolerance to abiotic stress (Zhang et al., 2020); CaNAC064, which positively modulates cold tolerance (Hou et al., 2020); and CaNAC2, which is mainly expressed in seeds and roots, and whose silencing results in increased susceptibility to cold stress (Guo et al., 2015).
On the basis of our RNA sequencing data, only two NACs, i.e., CaNAC072 and CaNAC104, were always upregulated after drought stress and downregulated after recovery, regardless of the stage of plant development; furthermore, following their transient expression in Nicotiana benthamiana leaves, both TFs were localized in the nucleus. Our results thus suggested a possible role of CaNAC072 and CaNAC104 as drought response switchers, and we, therefore, focused our attention on them.
The upregulation of the two NACs was induced at different times following treatments with NaCl, PEG and ABA, whereby CaNAC072 appeared to be a more rapid responder to abiotic stress than CaNAC104. This appears to be confirmed by the not significant upregulation of the slower responsive gene CaNAC104 after PEG treatment, likely due to a less marked drought stress induced by this treatment compared to the one we imposed (i.e., stop of irrigation up to leaf withering) before our transcriptome analyses.
VIGS represents a fast and powerful reverse genetics tool for analyzing the function of genes in many plant species. In the last decade, VIGS has been successfully applied in Capsicum for the functional characterization of the RING Type E3 ligase CaAIR1, involved in ABA signaling and drought stress response (Park et al., 2014), the ethylene-responsive TF CaAIEF1, involved in enhancing ABA sensitivity and drought tolerance (Hong et al., 2017) and, more recently, the CaNAC035 TF (Zhang et al., 2020). Here, a VIGS protocol, whose efficiency was validated by disabling phytoene desaturase (PDS) as a visible reporter/marker, was applied for transient gene silencing of CaNAC072 and CaNAC104. No difference in the tolerance to drought stress was observed in CaNAC104-silenced pepper in comparison to control plants, and a similar relative water content (RWC) and water loss was detected in detached leaves of both plants. This result indicates that the stress-responsive TF CaNAC104 is a component of a complex regulatory network activated by different abiotic stresses, but seemingly it does not play a key role in influencing the response to drought stress. However, considering that CaNAC104 is the ortholog of the Arabidopsis XND1 gene, which plays a major role in root xylem tracheids (vessels) (Tang et al., 2018), it cannot be excluded that its inhibition by VIGS was not effective and thus did not alter the response of the silenced plants in comparison to the control plants.
Unexpectedly, the silencing of CaNAC072 improved the tolerance of bell pepper to drought stress, as silenced plants did not display a wilting phenotype after 7 days of no watering, and showed a higher RWC and lower water loss of detached leaves in comparison to control plants. Notably, in contrast to CaNAC104-silenced plants, the silencing of CaNAC072 induced a marked upregulation of stress-responsive genes, including CaNAC019 and CaNAC055. Furthermore, after severe drought (9 days of water withholding), about 15% of the CaNAC072-silenced plants survived, while none of the TRV2:00 and CaNAC104-silenced plants did.
In Arabidopsis, both ANAC019 and ANAC055 are closely related and functionally redundant homologs of ANAC072, and they act as positive regulators of the drought response (Tran et al., 2004; Ye et al., 2017). On the basis of our RNA-seq experiment, bell pepper plants exposed to acute water stress episodes activated the expression of CaNAC072, but not of CaNAC019 and CaNAC055 (Supplementary Table S21). However, since the latter two were upregulated in drought-stressed CaNAC072-silenced plants, it can be assumed that their upregulation might compensate the disabling of CaNAC072, because of a complex gene regulatory network involved in the response to abiotic stresses (Bouché and Bouchez, 2001; El-Brolosy and Stainier, 2017). Similarly, in the Arabidopsis anac019 mutant, drought stress induced the upregulation of ANAC055 and ANAC072 (Sukiran et al., 2019).
CaNAC072 is the orthologue of the Arabidopsis gene ANAC072, a gene activated following drought stress and conferring drought stress tolerance when overexpressed (Tran et al., 2004). With the goal to shed light on its role, we generated transgenic Arabidopsis plants overexpressing CaNAC072 fused to a green fluorescent protein (GFP) coding sequence (CaNAC072:GFP-OX plants) in Col-0 as well as in anac072 single and anac019/anac055/anac072 triple mutants. After 20 days of water withholding, no difference in water stress tolerance between the Col-0 wild type and CaNAC072 overexpressors was observed. Furthermore, no recovery of function, i.e. increased drought tolerance, was observed when CaNAC072 was overexpressed in the two Arabidopsis mutant backgrounds. This suggests that the interactions between a NAC protein and the cis-acting elements of target genes as well as the gene regulatory pathway and downstream signaling can vary between plant species. Indeed, previous studies highlighted that orthologous genes might not play the same function in different species, like in the case of Arabidopsis CBF homologs identified in tomato (i.e., LeCBF1, LeCBF2 or LeCBF3) and known to be involved in cold acclimation, and of which only LeCBF1 increased freezing tolerance in transgenic Arabidopsis plants (Zhang et al., 2004).
Conclusion
Following RNA-sequencing, we identified two transcription factors, CaNAC072 and CaNAC104, that were always up-regulated by acute drought stress and down-regulated after recovery when water was again supplied, at three stages of bell pepper development. The transcriptional regulation and fine-tuning of NAC gene expression is determined by complex signaling pathways: a single NAC gene often responds to several stress factors and participates in the regulation of different processes, as a positive or negative regulator. In our study, the VIGS-mediated silencing of CaNAC104 did not affect drought tolerance, while silencing of CaNAC072 increased drought stress tolerance in pepper plantlets. An important result of our study is that CaNAC072 did not recover the function of its homolog ANAC072 in Arabidopsis thaliana. This observation suggests that the two TF bind to different cis-acting elements of target genes in the two plant species, or that interactions with other TFs or regulatory proteins that ANAC072 might undergo are not faithfully replicated by CaNAC072.
Future research may take advantage of CRISPR/Cas9-mediated genome editing to stably knockout the NAC genes in C. annuum for further characterization of their mode of action. This, however, requires the availability of a robust in vitro culture and regeneration protocol for this recalcitrant species, which we are currently establishing.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: National Center for Biotechnology Information (NCBI) BioProject database under accession number PRJNA668245.
Author Contributions
DB, SB, and SL planned the experiments. SB, IK, and SL provided plant material. DB carried out the drought stress experiments, sampling, phenotype analyses, plant transformation, and statistical analyses. KS performed the cloning and primer design. LB performed the bioinformatics analyses. DB, AA, SB, AM, LB, and SL contributed to data interpretation and preparation of graphic material. DB and SL wrote the manuscript and DB curated the supplementary materials. DB, SB, AM, AA, LB, and SL discussed the results and revised the manuscript.
Funding
This research was funded by the “Cassa di Risparmio di Cuneo” Foundation, under the research project: Resequencing of sweet pepper ecotypes for enhancing their quality and for their traceability (RISEPP).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
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.
Acknowledgments
The authors thank Anna Maria Milani, Danila Valentino, Dr. Shuchao Dong and Dr. Venkatesh Thirumalaikumar for technical support. This research was funded by the “Cassa di Risparmio di Cuneo” Foundation, under the research project: Resequencing of sweet pepper ecotypes for enhancing their quality and for their traceability (RISEPP). DB thanks the University of Turin for havingfunded his PhD fellowship. SB thanks the Max Planck Institute of Molecular Plant Physiology, Potsdam-Golm, for financial support. Funding of the ERA-CAPS project “AbioSen” by the Deutsche Forschungsgemeinschaft (MU 1199/16-1) is greatly acknowledged.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.743902/full#supplementary-material
References
Acquadro, A., Barchi, L., Portis, E., Nourdine, M., Carli, C., Monge, S., et al. (2020). Whole Genome Resequencing of Four Italian Sweet Pepper Landraces Provides Insights on Sequence Variation in Genes of Agronomic Value. Sci. Rep. 10, 1–16. doi:10.1038/s41598-020-66053-2
Ahn, Y.-K., Manivannan, A., Karna, S., Jun, T.-H., Yang, E.-Y., Choi, S., et al. (2018). Whole Genome Resequencing of Capsicum Baccatum and Capsicum annuum to Discover Single Nucleotide Polymorphism Related to Powdery Mildew Resistance. Sci. Rep. 8, 1–11. doi:10.1038/s41598-018-23279-5
Arvidsson, S., Kwasniewski, M., Riaño-Pachón, D. M., and Mueller-Roeber, B. (2008). QuantPrime - A Flexible Tool for Reliable High-Throughput Primer Design for Quantitative PCR. BMC Bioinformatics 9, 465. doi:10.1186/1471-2105-9-465
Balazadeh, S., Riaño-Pachón, D. M., and Mueller-Roeber, B. (2008). Transcription Factors Regulating Leaf Senescence in Arabidopsis thaliana. Plant Biol. 10, 63–75. doi:10.1111/j.1438-8677.2008.00088.x
Borràs, D., Plazas, M., Moglia, A., and Lanteri, S. (2021). The Influence of Acute Water Stresses on the Biochemical Composition of bell Pepper (Capsicum annuum L.) Berries. J. Sci. Food Agric. 101, 4724–4734. doi:10.1002/jsfa.11118
Bouché, N., and Bouchez, D. (2001). Arabidopsis Gene Knockout: Phenotypes Wanted. Curr. Opin. Plant Biol. 4, 111–117. doi:10.1016/S1369-5266(00)00145-X
Boyer, J. S., Byrne, P., Cassman, K. G., Cooper, M., Delmer, D., Greene, T., et al. (2013). The U.S. Drought of 2012 in Perspective: A Call to Action. Glob. Food Security 2, 139–143. doi:10.1016/j.gfs.2013.08.002
Chung, P. J., Jung, H., Choi, Y. D., and Kim, J.-K. (2018). Genome-wide Analyses of Direct Target Genes of Four rice NAC-Domain Transcription Factors Involved in Drought Tolerance. BMC Genomics 19, 40. doi:10.1186/s12864-017-4367-1
Clough, S. J., and Bent, A. F. (1998). Floral Dip: a Simplified Method for Agrobacterium-Mediated Transformation of Arabidopsis Thaliana. Plant J. 16, 735–743. doi:10.1046/j.1365-313x.1998.00343.x
Devkar, V., Thirumalaikumar, V. P., Xue, G. P., Vallarino, J. G., Turečková, V., Strnad, M., et al. (2020). Multifaceted Regulatory Function of Tomato SlTAF1 in the Response to Salinity Stress. New Phytol. 225, 1681–1698. doi:10.1111/nph.16247
Diao, W., Snyder, J., Wang, S., Liu, J., Pan, B., Guo, G., et al. (2018). Genome-wide Analyses of the NAC Transcription Factor Gene Family in Pepper (Capsicum annuum L.): Chromosome Location, Phylogeny, Structure, Expression Patterns, Cis-Elements in the Promoter, and Interaction Network. Ijms 19, 1028. doi:10.3390/ijms19041028
Dominguez, P. G., and Carrari, F. (2015). ASR1 Transcription Factor and its Role in Metabolism. Plant Signaling Behav. 10, e992751–5. doi:10.4161/15592324.2014.992751
El-Brolosy, M. A., and Stainier, D. Y. R. (2017). Genetic Compensation: A Phenomenon in Search of Mechanisms. PLoS Genet. 13, e1006780. doi:10.1371/journal.pgen.1006780
El-Gebali, S., Mistry, J., Bateman, A., Eddy, S. R., Luciani, A., Potter, S. C., et al. (2019). The Pfam Protein Families Database in 2019. Nucleic Acids Res. 47, D427–D432. doi:10.1093/nar/gky995
Fernandez-Pozo, N., Rosli, H. G., Martin, G. B., and Mueller, L. A. (2015). The SGN VIGS Tool: User-Friendly Software to Design Virus-Induced Gene Silencing (VIGS) Constructs for Functional Genomics. Mol. Plant 8, 486–488. doi:10.1016/j.molp.2014.11.024
Fitton, N., Alexander, P., Arnell, N., Bajzelj, B., Calvin, K., Doelman, J., et al. (2019). The Vulnerabilities of Agricultural Land and Food Production to Future Water Scarcity. Glob. Environ. Change 58, 101944. doi:10.1016/j.gloenvcha.2019.101944
Franco-Zorrilla, J. M., López-Vidriero, I., Carrasco, J. L., Godoy, M., Vera, P., and Solano, R. (2014). DNA-binding Specificities of Plant Transcription Factors and Their Potential to Define Target Genes. Proc. Natl. Acad. Sci. USA 111, 2367–2372. doi:10.1073/pnas.1316278111
Gamboa-Tuz, S. D., Pereira-Santana, A., Zamora-Briseño, J. A., Castano, E., Espadas-Gil, F., Ayala-Sumuano, J. T., et al. (2018). Transcriptomics and Co-expression Networks Reveal Tissue-specific Responses and Regulatory Hubs under Mild and Severe Drought in Papaya (Carica Papaya L.). Sci. Rep. 8, 1–16. doi:10.1038/s41598-018-32904-2
Godoy, J. A., Lunar, R., Torres-Schumann, S., Moreno, J., Rodrigo, R. M., and Pintor-Toro, J. A. (1994). Expression, Tissue Distribution and Subcellular Localization of Dehydrin TAS14 in Salt-Stressed Tomato Plants. Plant Mol. Biol. 26, 1921–1934. doi:10.1007/BF00019503
Gong, P., Zhang, J., Li, H., Yang, C., Zhang, C., Zhang, X., et al. (2010). Transcriptional Profiles of Drought-Responsive Genes in Modulating Transcription Signal Transduction, and Biochemical Pathways in Tomato. J. Exp. Bot. 61, 3563–3575. doi:10.1093/jxb/erq167
Guo, W.-L., Wang, S.-B., Chen, R.-G., Chen, B.-H., Du, X.-H., Yin, Y.-X., et al. (2015). Characterization and Expression Profile of CaNAC2 Pepper Gene. Front. Plant Sci. 6, 755. doi:10.3389/fpls.2015.00755
Han, Q., Zhang, J., Li, H., Luo, Z., Ziaf, K., Ouyang, B., et al. (2012). Identification and Expression Pattern of One Stress-Responsive NAC Gene from Solanum lycopersicum. Mol. Biol. Rep. 39, 1713–1720. doi:10.1007/s11033-011-0911-2
Hong, E., Lim, C. W., Han, S.-W., and Lee, S. C. (2017). Functional Analysis of the Pepper Ethylene-Responsive Transcription Factor, CaAIEF1, in Enhanced ABA Sensitivity and Drought Tolerance. Front. Plant Sci. 8, 1407. doi:10.3389/fpls.2017.01407
Hou, X.-m., Zhang, H.-f., Liu, S.-y., Wang, X.-k., Zhang, Y.-m., Meng, Y.-c., et al. (2020). The NAC Transcription Factor CaNAC064 Is a Regulator of Cold Stress Tolerance in Peppers. Plant Sci. 291, 110346. doi:10.1016/j.plantsci.2019.110346
Iovieno, P., Punzo, P., Guida, G., Mistretta, C., Van Oosten, M. J., Nurcato, R., et al. (2016). Transcriptomic Changes Drive Physiological Responses to Progressive Drought Stress and Rehydration in Tomato. Front. Plant Sci. 7, 371. doi:10.3389/fpls.2016.00371
Ji, X., Shiran, B., Wan, J., Lewis, D. C., Jenkins, C. L. D., Condon, A. G., et al. (2010). Importance of Pre-anthesis Anther Sink Strength for Maintenance of Grain Number during Reproductive Stage Water Stress in Wheat. Plant Cel Environ. 33, 926–942. doi:10.1111/j.1365-3040.2010.02130.x
Kakumanu, A., Ambavaram, M. M. R., Klumas, C., Krishnan, A., Batlang, U., Myers, E., et al. (2012). Effects of Drought on Gene Expression in Maize Reproductive and Leaf Meristem Tissue Revealed by RNA-Seq. Plant Physiol. 160, 846–867. doi:10.1104/pp.112.200444
Kang, Y. J., Ahn, Y.-K., Kim, K.-T., and Jun, T.-H. (2016). Resequencing of Capsicum Annuum Parental Lines (YCM334 and Taean) for the Genetic Analysis of Bacterial Wilt Resistance. BMC Plant Biol. 16, 235. doi:10.1186/s12870-016-0931-0
Karimi, M., Inzé, D., and Depicker, A. (2002). GATEWAY Vectors for Agrobacterium-Mediated Plant Transformation. Trends Plant Sci. 7, 193–195. doi:10.1016/S1360-1385(02)02251-3
Kilian, J., Whitehead, D., Horak, J., Wanke, D., Weinl, S., Batistic, O., et al. (2007). The AtGenExpress Global Stress Expression Data Set: Protocols, Evaluation and Model Data Analysis of UV-B Light, Drought and Cold Stress Responses. Plant J. 50, 347–363. doi:10.1111/j.1365-313X.2007.03052.x
Kim, D., Langmead, B., and Salzberg, S. L. (2015). hisat2. Nat. Methods 12, 357–360. doi:10.1038/nmeth.3317
Kim, S., Park, J., Yeom, S.-I., Kim, Y.-M., Seo, E., Kim, K.-T., et al. (2017). New Reference Genome Sequences of Hot Pepper Reveal the Massive Evolution of Plant Disease-Resistance Genes by Retroduplication. Genome Biol. 18, 1–11. doi:10.1186/s13059-017-1341-9
Kim, S., Park, M., Yeom, S.-I., Kim, Y.-M., Lee, J. M., Lee, H.-A., et al. (2014). Genome Sequence of the Hot Pepper Provides Insights into the Evolution of Pungency in Capsicum Species. Nat. Genet. 46, 270–278. doi:10.1038/ng.2877
Lee, J., Shim, D., Moon, S., Kim, H., Bae, W., Kim, fnm., et al. (2018). Genome-wide Transcriptomic Analysis of BR-Deficient Micro-tom Reveals Correlations between Drought Stress Tolerance and Brassinosteroid Signaling in Tomato. Plant Physiol. Biochem. 127, 553–560. doi:10.1016/j.plaphy.2018.04.031
Lee, S., and Choi, D. (2013). Comparative Transcriptome Analysis of Pepper (Capsicum annuum) Revealed Common Regulons in Multiple Stress Conditions and Hormone Treatments. Plant Cel Rep. 32, 1351–1359. doi:10.1007/s00299-013-1447-9
Liu, Y., Schiff, M., and Dinesh-Kumar, S. P. (2002). Virus-induced Gene Silencing in Tomato. Plant J. 31, 777–786. doi:10.1046/j.1365-313X.2002.01394.x
Love, M. I., Huber, W., and Anders, S. (2014). Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. Genome Biol. 15, 550. doi:10.1186/s13059-014-0550-8
Ma, N.-N., Zuo, Y.-Q., Liang, X.-Q., Yin, B., Wang, G.-D., and Meng, Q.-W. (2013). The Multiple Stress-Responsive Transcription Factor SlNAC1 improves the Chilling Tolerance of Tomato. Physiol. Plantarum 149, 474–486. doi:10.1111/ppl.12049
Madeira, F., Park, Y. M., Lee, J., Buso, N., Gur, T., Madhusoodanan, N., et al. (2019). The EMBL-EBI Search and Sequence Analysis Tools APIs in 2019. Nucleic Acids Res. 47 (W1), W636–W641. doi:10.1093/nar/gkz268
Mao, C., Ding, W., Wu, Y., Yu, J., He, X., Shou, H., et al. (2007). Overexpression of a NAC‐domain Protein Promotes Shoot Branching in Rice. New Phytol. 176, 288–298. doi:10.1111/j.1469-8137.2007.0217710.1111/j.1469-8137.2007.02177.x
Marusig, D., and Tombesi, S. (2020). Abscisic Acid Mediates Drought and Salt Stress Responses in Vitis Vinifera-A Review. Ijms 21, 8648. doi:10.3390/ijms21228648
Mitchell, A. L., Attwood, T. K., Babbitt, P. C., Blum, M., Bork, P., Bridge, A., et al. (2019). InterPro in 2019: Improving Coverage, Classification and Access to Protein Sequence Annotations. Nucleic Acids Res. 47, D351–D360. doi:10.1093/nar/gky1100
Mittal, S., Banduni, P., Mallikarjuna, M. G., Rao, A. R., Jain, P. A., Dash, P. K., et al. (2018). Structural, Functional, and Evolutionary Characterization of Major Drought Transcription Factors Families in Maize. Front. Chem. 6, 177. doi:10.3389/fchem.2018.00177
Mohanta, T. K., Yadav, D., Khan, A., Hashem, A., Tabassum, B., Khan, A. L., et al. (2020). Genomics, Molecular and Evolutionary Perspective of NAC Transcription Factors. PLoS One 15, e0231425. doi:10.1371/journal.pone.0231425
Msanne, J., Lin, J., Stone, J. M., and Awada, T. (2011). Characterization of Abiotic Stress-Responsive Arabidopsis thalianaRD29A and RD29B Genes and Evaluation of Transgenes. Planta 234, 97–107. doi:10.1007/s00425-011-1387-y
Muñoz-Mayor, A., Pineda, B., Garcia-Abellán, J. O., Antón, T., Garcia-Sogo, B., Sanchez-Bel, P., et al. (2012). Overexpression of Dehydrin Tas14 Gene Improves the Osmotic Stress Imposed by Drought and Salinity in Tomato. J. Plant Physiol. 169, 459–468. doi:10.1016/j.jplph.2011.11.018
Nakashima, K., Takasaki, H., Mizoi, J., Shinozaki, K., and Yamaguchi-Shinozaki, K. (2012). NAC Transcription Factors in Plant Abiotic Stress Responses. Biochim. Biophys. Acta (Bba) - Gene Regul. Mech. 1819, 97–103. doi:10.1016/j.bbagrm.2011.10.005
Newton, A. C., Johnson, S. N., Gregory, P. J., Newton, C., Johnson, N., and Gregory, J. (2011). Implications of Climate Change for Diseases, Crop Yields and Food Security. Euphytica 179, 3–18. doi:10.1007/s10681-011-0359-4
Olvera-Carrillo, Y., Campos, F., Reyes, J. L., Garciarrubio, A., and Covarrubias, A. A. (2010). Functional Analysis of the Group 4 Late Embryogenesis Abundant Proteins Reveals Their Relevance in the Adaptive Response during Water Deficit in Arabidopsis. Plant Physiol. 154, 373–390. doi:10.1104/pp.110.158964
Osakabe, Y., Osakabe, K., Shinozaki, K., and Tran, L.-S. P. (2014). Response of Plants to Water Stress. Front. Plant Sci. 5, 86. doi:10.3389/fpls.2014.00086
Park, C., Lim, C. W., Baek, W., and Lee, S. C. (2015). RING Type E3 Ligase CaAIR1 in Pepper Acts in the Regulation of ABA Signaling and Drought Stress Response. Plant Cel Physiol 56, 1808–1819. doi:10.1093/pcp/pcv103
Pertea, M., Pertea, G. M., Antonescu, C. M., Chang, T.-C., Mendell, J. T., and Salzberg, S. L. (2015). StringTie Enables Improved Reconstruction of a Transcriptome from RNA-Seq Reads. Nat. Biotechnol. 33, 290–295. doi:10.1038/nbt.3122
Qin, C., Yu, C., Shen, Y., Fang, X., Chen, L., Min, J., et al. (2014). Whole-genome Sequencing of Cultivated and Wild Peppers Provides Insights into Capsicum Domestication and Specialization. Proc. Natl. Acad. Sci. U. S. A. 111, 5135–5140. doi:10.1073/pnas.1400975111
Sakuma, Y., Liu, Q., Dubouzet, J. G., Abe, H., Shinozaki, K., and Yamaguchi-Shinozaki, K. (2002). DNA-binding Specificity of the ERF/AP2 Domain of Arabidopsis DREBs, Transcription Factors Involved in Dehydration- and Cold-Inducible Gene Expression. Biochem. Biophysical Res. Commun. 290, 998–1009. doi:10.1006/bbrc.2001.6299
Sampathkumar, A., and Wightman, R. (2015). Live Cell Imaging of the Cytoskeleton and Cell wall Enzymes in Plant Cells. Methods Mol. Biol. 1242, 133–141. doi:10.1007/978-1-4939-1902-4_12
Seki, M., Ishida, J., Narusaka, M., Fujita, M., Nanjo, T., Umezawa, T., et al. (2002). Monitoring the Expression Pattern of Around 7,000 Arabidopsis Genes under ABA Treatments Using a Full-Length cDNA Microarray. Funct. Integr. Genomics 2, 282–291. doi:10.1007/s10142-002-0070-6
Senthil-Kumar, M., and Mysore, K. S. (2014). Tobacco Rattle Virus-Based Virus-Induced Gene Silencing in Nicotiana benthamiana. Nat. Protoc. 9, 1549–1562. doi:10.1038/nprot.2014.092
Sharoni, A. M., Nuruzzaman, M., Satoh, K., Moumeni, A., Attia, K., Venuprasad, R., et al. (2012). Comparative Transcriptome Analysis of AP2/EREBP Gene Family under normal and Hormone Treatments, and under Two Drought Stresses in NILs Setup by Aday Selection and IR64. Mol. Genet. Genomics 287, 1–19. doi:10.1007/s00438-011-0659-3
Shinozaki, K., and Yamaguchi-Shinozaki, K. (2007). Gene Networks Involved in Drought Stress Response and Tolerance. J. Exp. Bot. 58, 221–227. doi:10.1093/jxb/erl164
Sinclair, T. R. (2011). Challenges in Breeding for Yield Increase for Drought. Trends Plant Sci. 16, 289–293. doi:10.1016/j.tplants.2011.02.008
Singh, D., and Laxmi, A. (2015). Transcriptional Regulation of Drought Response: A Tortuous Network of Transcriptional Factors. Front. Plant Sci. 6, 895. doi:10.3389/fpls.2015.00895
Söderman, E., Mattsson, J., and Engström, P. (1996). The Arabidopsis Homeobox Gene ATHB-7 Is Induced by Water Deficit and by Abscisic Acid. Plant J. 10, 375–381. doi:10.1046/j.1365-313X.1996.10020375.x
Sprenger, H., Erban, A., Seddig, S., Rudack, K., Thalhammer, A., Le, M. Q., et al. (2018). Metabolite and Transcript Markers for the Prediction of Potato Drought Tolerance. Plant Biotechnol. J. 16, 939–950. doi:10.1111/pbi.12840
Sukiran, N. L., Ma, J. C., Ma, H., and Su, Z. (2019). ANAC019 Is Required for Recovery of Reproductive Development under Drought Stress in Arabidopsis. Plant Mol. Biol. 99, 161–174. doi:10.1007/s11103-018-0810-1
Tang, N., Shahzad, Z., Lonjon, F., Loudet, O., Vailleau, F., and Maurel, C. (2018). Natural Variation at XND1 Impacts Root Hydraulics and Trade-Off for Stress Responses in Arabidopsis. Nat. Commun. 9, 1–12. doi:10.1038/s41467-018-06430-8
Thirumalaikumar, V. P., Devkar, V., Mehterov, N., Ali, S., Ozgur, R., Turkan, I., et al. (2018). NAC Transcription Factor JUNGBRUNNEN1 Enhances Drought Tolerance in Tomato. Plant Biotechnol. J. 16, 354–366. doi:10.1111/pbi.12776
Todaka, D., Takahashi, F., Yamaguchi-Shinozaki, K., and Shinozaki, K. (2019). ABA-responsive Gene Expression in Response to Drought Stress: Cellular Regulation and Long-Distance Signaling. Adv. Bot. Res. 92, 83–113. doi:10.1016/bs.abr.2019.05.001
Tran, L.-S. P., Nakashima, K., Sakuma, Y., Simpson, S. D., Fujita, Y., Maruyama, K., et al. (2004). Isolation and Functional Analysis of Arabidopsis Stress-Inducible NAC Transcription Factors that Bind to a Drought-Responsive Cis-Element in the Early Responsive to Dehydration Stress 1 Promoter[W]. Plant Cell 16, 2481–2498. doi:10.1105/tpc.104.022699
Trujillo, L. E., Sotolongo, M., Menéndez, C., Ochogavía, M. E., Coll, Y., Hernández, I., et al. (2008). SodERF3, a Novel Sugarcane Ethylene Responsive Factor (ERF), Enhances Salt and Drought Tolerance when Overexpressed in Tobacco Plants. Plant Cel Physiol 49, 512–525. doi:10.1093/pcp/pcn025
Tuteja, N. (2007). Abscisic Acid and Abiotic Stress Signaling. Plant Signaling Behav. 2, 135–138. doi:10.4161/psb.2.3.4156
Yang, R., Deng, C., Ouyang, B., and Ye, Z. (2011). Molecular Analysis of Two Salt-Responsive NAC-Family Genes and Their Expression Analysis in Tomato. Mol. Biol. Rep. 38, 857–863. doi:10.1007/s11033-010-0177-0
Ye, H., Liu, S., Tang, B., Chen, J., Xie, Z., Nolan, T. M., et al. (2017). RD26 Mediates Crosstalk between Drought and Brassinosteroid Signalling Pathways. Nat. Commun. 8, 14573. doi:10.1038/ncomms14573
Zhang, H., Ma, F., Wang, X., Liu, S., Saeed, U. H., Hou, X., et al. (2020). Molecular and Functional Characterization of CaNAC035, an NAC Transcription Factor from Pepper (Capsicum annuum L.). Front. Plant Sci. 11, 14. doi:10.3389/fpls.2020.00014
Zhang, P., Fan, Y., Sun, X., Chen, L., Terzaghi, W., Bucher, E., et al. (2019). A Large‐scale Circular RNA Profiling Reveals Universal Molecular Mechanisms Responsive to Drought Stress in maize and Arabidopsis. Plant J. 98, 697–713. doi:10.1111/tpj.14267
Zhang, X., Fowler, S. G., Cheng, H., Lou, Y., Rhee, S. Y., Stockinger, E. J., et al. (2004). Freezing-sensitive Tomato Has a Functional CBF Cold Response Pathway, but a CBF Regulon that Differs from that of Freezing-tolerantArabidopsis. Plant J. 39, 905–919. doi:10.1111/j.1365-313X.2004.02176.x
Keywords: bell pepper, transcriptome, drought tolerance, NAC, VIGS, functional characterization
Citation: Borràs D, Barchi L, Schulz K, Moglia A, Acquadro A, Kamranfar I, Balazadeh S and Lanteri S (2021) Transcriptome-Based Identification and Functional Characterization of NAC Transcription Factors Responsive to Drought Stress in Capsicum annuum L.. Front. Genet. 12:743902. doi: 10.3389/fgene.2021.743902
Received: 19 July 2021; Accepted: 28 September 2021;
Published: 22 October 2021.
Edited by:
Nunzio D’Agostino, University of Naples Federico II, ItalyReviewed by:
Neftali Ochoa-Alejo, Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional, MexicoAllen Van Deynze, University of California, Davis, United States
Copyright © 2021 Borràs, Barchi, Schulz, Moglia, Acquadro, Kamranfar, Balazadeh and Lanteri. 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.
*Correspondence: Sergio Lanteri, c2VyZ2lvLmxhbnRlcmlAdW5pdG8uaXQ=
†ORCID: Dionis Borràs, orcid.org/0000-0003-2793-7832; Lorenzo Barchi, orcid.org/0000-0001-6414-4239; Karina Schulz, orcid.org/0000-0002-6540-2209; Andrea Moglia, orcid.org/0000-0001-7431-844X; Alberto Acquadro, orcid.org/0000-0002-5322-9701; Iman Kamranfar, orcid.org/0000-0001-9632-3717; Salma Balazadeh, orcid.org/0000-0002-5789-4071; Sergio Lanteri, orcid.org/0000-0003-3012-8710