- 1Department of Plant Pathology, University of Agriculture Faisalabad, Faisalabad, Pakistan
- 2Centre of Agricultural Biochemistry and Biotechnology, University of Agriculture Faisalabad, Faisalabad, Pakistan
- 3Botany and Microbiology Department, College of Science, King Saud University, Riyadh, Saudi Arabia
- 4Facultad de Ciencias Agrotecnológicas, Universidad Autónoma de Chihuahua, Chihuahua, Chihuahua, Mexico
- 5Plant Production Department, College of Food and Agricultural Sciences, King Saud University, Riyadh, Saudi Arabia
Quantitative real-time PCR is used to quantify gene expression, even to detect low-level transcripts. It detects and quantifies the inoculum level of fungal pathogens in infected hosts. However, reliable expression profiling data require accurate transcript normalization against a stable reference gene. Hence, using stably expressed reference genes under variable conditions is paramount in gene expression analysis. In the current study, reference genes were selected and validated in Colletotrichum gloeosporioides, a guava canker and dieback pathogen. The reference gene selection and validation in C. gloeosporioides were evaluated for germinated conidia and mycelium (in vitro) and in infected guava (Psidium guajava) (interaction with host plant). The CgCAL gene was determined as a highly stable reference gene, followed by the CgTUB2 in C. gloeosporioides for germinating conidia and mycelium. However, the CgTUB2 gene was determined to be a highly stable reference gene, followed by the CgCAL for expression analysis during its interaction with the plant. Expression profiling revealed stable and constant relative expression patterns of selected reference genes for both PR genes by determining their relative transcript level. This study is the first to describe reference gene selection and validation to quantify target gene expression in C. gloeosporioides.
Introduction
The genus Colletotrichum includes a broad array of economically important fungal pathogens that infect various host plant species (Yan et al., 2018; Iqra et al., 2022). Fungal biology, genomics, genetics, colonization, virulence factors, and interaction with host plants are prerequisites to mitigate compromised quality and consequent financial losses in the agricultural sector (Villa-Rivera et al., 2017; Babar et al., 2021). For this purpose, model species have been extensively studied for insight into infection stages and host interactions (Haq and Ijaz, 2020). Different techniques are involved in expression profiling, including northern blotting, semiquantitative real-time PCR, microarray, RNAseq, and quantitative real-time PCR (qRT−PCR) (Ijaz et al., 2020a). Among these, qRT-PCR is the most efficient technique for determining the relative levels of target gene expression in different samples (MaChado et al., 2015; Ijaz et al., 2019). A notable advantage of qRT-PCR over other conventional techniques is its ability to detect low starting material copies of the targeted gene’s mRNA (Huggett et al., 2005; Ijaz and Haq, 2020). Nonetheless, the results obtained through qRT-PCR depend upon the accuracy of target transcript normalization using appropriate reference genes for controlling significant experimental errors (Radonić et al., 2004).
Functional gene expression quantification of the targeted gene is crucial for comprehensive gene transcription and regulation studies through qRT-PCR (Marcial-Quino et al., 2016). A systemic study quantified crucial functional gene expression aspects using qRT-PCR for the targeted gene. This technique for detecting gene expression and profiling analysis is simple, fast, highly sensitive, and reproducible (Hao et al., 2014; Haq et al., 2019). However, data analysis requires suitable reference genes, pivotal in generating reliable qRT-PCR results. In contrast, incompetent reference genes produce misleading expression data (Galli et al., 2015). A reliable reference gene should be constitutively expressed in experimental samples, with unaltered expression under diverse experimental sets (Galli et al., 2015). For Colletotrichum members, several commonly used reference genes have been documented; for example, in Colletotrichum higginasianum, the actin (ACT) gene was used to quantify and normalize fungal growth and transporter (MFS) gene expression (Liu et al., 2017). However, in Colletotrichum fructicola, the Ras-related protein gene (CfRrp) showed better stability (Chen et al., 2022) in quantifying target gene expression. However, to our knowledge, no literature and studies have been found for the reference gene selection to normalize qRT-PCR gene expression analysis in Colletotrichum gloeosporioides.
Studies have illustrated that reference gene expression stability varies considerably between species and could differ across tissue types, stages of development, and experimental conditions (Guenin et al., 2009). Therefore, a constitutively expressed gene is essential for obtaining reliable quantitative real-time PCR results, which should be stable at all growth stages and may serve as internal controls. The most suitable and preferred technique for analyzing gene expression is qRT-PCR. However, in C. gloeosporioides, no suitable, reliable reference genes have been reported for qRT-PCR, but their presence is expected in different experimental conditions, which need to be identified. C. gloeosporioides is a predominant fungal pathogen of guava plants (Lu et al., 2018; Iqra et al., 2022). It is responsible for leaf and other severe guava plant diseases, notably leaf blight, brown blight, and anthracnose (Wang et al., 2016; Haq and Ijaz, 2019; Haq et al., 2021). It is the first study to select and validate suitable reference genes in C. gloeosporioides for transcript normalization in qRT-PCR. Housekeeping genes, including actin (ACT), Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), 18S ribosomal RNA (18S rRNA), calmodulin (CAL), and β-tubulin (TUB2), are considered the most appropriate candidate reference genes for fungal plant pathogens because of their involvement in essential cellular functions (Galli et al., 2015; Haq and Ijaz, 2020). In the present study, four candidate reference genes were appraised for their identification and suitability for transcript normalization in C. gloeosporioides at germinating conidia and mycelial stage and interaction with Psidium guajava. Their stability evaluation was performed with Excel-based software. These data analyses extensively validated the putative reference gene, which can help further gene expression studies associated with this fungus.
Materials and methods
Fungal culture
Germinating conidia and mycelial stage of fungal pathogen Colletotrichum gloeosporioides were considered for selecting suitable reference genes for transcript normalization in qRT−PCR. C. gloeosporioides was cultured with a sterile needle on a Potato Dextrose Agar (PDA) medium to get germinate conidia and mycelium. For this purpose, single spore and hyphae tip methods were employed. The culture plates were incubated at 27°C under a dark regime for 48-72 hrs, and growth was observed. Microscopical features were examined under stereomicroscopy.
Plant inoculation
For reference gene selection, while the interaction of C. gloeosporioides with guava pants, fungal inoculation of plants was performed. A Guava cultivar named “Large Sorahie” (Pyriform) was used in this experiment. Plants were raised in the research area of the Department of Plant Pathology, University of Agriculture Faisalabad, under the greenhouse-controlled condition. Guava plants at the seedling stage (six months old) were inoculated by spraying the plants with a spore suspension (150 ml), which was adjusted to a concentration of 1×106 spores ml−1 (Chen et al., 2021).
RNA isolation
Total RNA was isolated from germinated conidia and mycelium of C. gloeosporioides and inoculated guava leaves sampled at two, three, and four weeks postinoculation. The GeneJET Plant RNA Purification kit (Thermo Scientific USA) was used for RNA extraction. The quantification of extracted RNA was determined using a UV-visible NANODROP 8000 Spectrophotometer (Thermo Scientific USA); however, integrity was analyzed by gel electrophoresis using 1.5% agarose. The isolated RNA from each sample was treated with a RapidOut DNA Removal kit (Thermo Scientific USA) to remove genomic DNA (gDNA).
cDNA synthesis
Treated RNA samples were reverse transcribed with RevertAid M-MuLV RT using a RevertAid First Strand cDNA Synthesis Kit (Thermo Scientific USA) in a final volume, as per the manufacturer’s recommendations, 20-fold diluted and quantified by UV visible NANODROP 8000 Spectrophotometer (Thermo Scientific USA).
Quantitative real-time PCR analysis
Five housekeeping genes were assessed as candidate reference genes in C. gloeosporioides, including calmodulin (CgCAL), Glyceraldehyde-3-phosphate dehydrogenase (CgGAPDH), 18S rRNA, β-tubulin (CgTUB2), and actin (CgACT). The primers employed for qRT-PCR analysis were designed using PrimerQuest. To verify their specificity for C. gloeosporioides, we used general PCR to test amplification in guava. Only primers with no amplification for guava were retained for qRT−PCR analysis (Table 1). This analysis was carried out in a CFX96 Touch Real-Time PCR detection system using a 96 × 0.2 ml plate (Bioplastics Netherlands) with a 5 μl total reaction volume for each sample, containing cDNA (1 μl), Maxima Syber green/ROX/qPCR master mix (2.5 μl) (Thermo Scientific USA), and 1 μM of each primer. For the negative control, water was used in the reaction instead of cDNA (no-template control). The thermal profile for qRT−PCR was denaturation at 95°C (60 sec), 40 amplification cycles at 95°C (20 sec), and an annealing/extension step at 60°C (30 sec). The melting curve analysis was employed at 60°C to 95°C. For each sample, three biological replicates and three technical replicates for each well were made. No-template controls (NTCs) (cDNA zero) were also assessed for each primer pair of candidate reference genes. The data were analyzed using CFX Manager ™ Software v. 3.1 (Bio-Rad Laboratories, Inc. USA).
Stability analysis and validation of selected reference genes
Three mathematical algorithms, geNorm (https://genorm.cmgg.be), NormFinder (https://moma.dk/normfnder-sofware), and BestKeeper (https://www.gene-quantifcation.de/bestkeeper.html), were applied to assess the candidate reference genes’ expression stability for transcript normalization. The best-ranked reference genes in C. gloeosporioides were assessed; the Cq (quantification cycle) conversion value to relative quantitates, except for the BestKeeper algorithm, was made. The overall comprehensive ranking was generated through the RefFinder tool (https://heartcure.com.au). The selected reference genes were validated for normalizing the expression data of qRT−PCR in C. gloeosporioides using the diacylglycerol acyltransferase (DGAT1) gene and 3-dehydroshikimate dehydratase (qutC) gene in qRT−PCR analysis. These pathogenicity-related genes were documented by Liang et al. (2018) in a fungal species of the C. gloeosporioides species complex.
Results
Expression analysis
Target DNA fragments were amplified from the cDNA of C. gloeosporioides by tested primers. Among these, 18S rRNA gave amplification in noninfected guava plants; however, the remaining four tested primers showed no amplification in noninfected guava leaves, which confirmed the absence of C. gloeosporioides in the control. Hence, the 18S rRNA gene was not used as a candidate reference gene for the qRT−PCR study in C. gloeosporioides. Melting curve analysis gave a single peak in all test primers and validated the primer specificity for qRT-PCR analysis. The Cq value was used to determine the gene expression levels, and it has an inverse relationship with gene expression, as a high Cq value reflects low gene expression and vice versa. The gene expression directly in the RNA sample of a fungal pathogen and RNA samples in an infected plant (related to the concentration of that fungal pathogen) cannot be exactly matched. However, it could be similar in these samples. That is why the expression analysis directly in the pathogen and the pathogen in the infected plant could not be compared (Song et al., 2019). The expression levels of the candidate reference genes, CgCAL, CgTUB2, and CgGAPDH, were similar in all samples (Figure 1), except for CgACT, which showed upregulation in infected plant material. Therefore, we discarded this housekeeping gene in further analyses of selecting reference gene(s) for transcript normalization in C. gloeosporioides. For the statistical evaluation of the candidate reference genes, three statistical algorithms were employed to rank the CgCAL, CgTUB2, and CgGAPDH genes based on expression stability. For expression stability analysis, we selected infected leaves of the guava plant after the third and fourth weeks postinoculation. The pathogen’s inoculum level in infected leaves of these periods was higher than that of infected leaves two weeks postinoculation (WPI) based on symptoms that appeared.
Figure 1 qRT−PCR analysis: (A) expression profiling of candidate reference genes in Colletotrichum gloeosporioides using qRT−PCR, (B) expression analysis of candidate genes in pathogen and infected leaves. Boxes denote lower and upper quartiles of the cycle threshold range with medians indicated; whisker caps symbolize minimum and maximum values. White boxes correspond to pathogen samples and gray boxes to infected leaf samples.
Stability analysis of candidate reference genes for normalizing qRT−PCR in Colletotrichum gloeosporioides
A study was carried out to select reference genes (CgGAPDH, CgCAL, and CgTUB2) for normalizing qRT−PCR expression analysis of C. gloeosporioides at four different levels (germinated conidia, mycelium, and infected leaves of guava plants at three and four weeks after inoculation). Three Excel-based statistical software programs were used to select a highly stable reference gene for normalization in expression profiling. The generated Cq values were converted into relative quantities for the analysis. The geNorm analysis revealed that the candidate reference genes CgTUB2 and CgCAL in combination were projected to have the least stability values, “M,” in germinating conidia and mycelium (0.040) and infected leaves (0.025). They ranked as the best reference genes with the highest stability for accurate data normalization in the qRT-PCR analysis of C. gloeosporioides (Table 2).
Table 2 geNorm analysis expression stability of candidate reference genes (a) in germinated conidia and mycelium of Colletotrichum gloeosporioides (b) in infected leaves of guava plants.
The descriptive data of the candidate reference genes based on expression in different samples are given in Table 3. In germinated conidia and mycelium of C. gloeosporioides, the candidate reference gene CgTUB2 with the lowest CV (% CP) ± SD value of 3.43 ± 0.61 was identified as the most stable reference gene for transcript normalization. However, CgGAPDH (5.52 ± 0.67) and CgCAL (4.83 ± 0.71) were ranked as the second and third reference genes, respectively, with stable expression for normalizing qRT−PCR analysis. In the case of infected leaves of guava plants, the ranking was observed; CgCAL ranked as the top best reference gene, with a CV (% CP) ± SD score of 3.86 ± 0.66, followed by CgTUB2, with a value of 4.09 ± 0.78. However, in the BestKeeper ranking, CgGAPDH was proven to be a reference gene with unstable expression by scoring an 8.82 ± 1.21 value because a gene with SD >1 is considered an unstable expression (Table 3).
Table 3 BestKeeper analysis of the expression stability of candidate reference genes: (a) in germinated conidia and mycelium of Colletotrichum gloeosporioides and (b) in infected leaves of guava plants.
NormFinder analysis showed that the CgCAL gene had the most stable expression, with a stability value of 0.020 in germinating conidia and mycelium; this gene’s same stability pattern was observed in infected leaves, with a stability value of 0.013. The ranking of candidate reference genes based on expression stability in germinating conidia and mycelium was CgCAL> CgTUB2>CgGAPDH. The same ranking was found in the NormFinder analysis of infected leaves of the guava plant (Table 4).
Table 4 NormFinder analysis of the expression stability of candidate reference genes: (a) in germinating conidia and mycelium of Colletotrichum gloeosporioides and (b) in infected leaves of guava plants.
In germinating conidia and mycelium, the results obtained using BestKeeper indicated that the reference gene CgTUB2 was stable in C. gloeosporioides. The result from NormFinder revealed the CgCAL gene as a stable reference gene, although geNorm gave the combination of reference genes, CgTUB2 and CgCAL, for normalization of gene expression analysis. However, the comprehensive ranking of candidate reference genes was determined using the RefFinder program, which ranked the CgCAL gene as the most stable reference gene for expression analysis in C. gloeosporioides in all experimental sets (germinating conidia, mycelium, and infected leaves of guava plants after three weeks and four weeks of inoculation). The comprehensive ranking order, CgCAL>CgTUB2>CgGAPDH, for stable reference genes for transcript normalization of qRT−PCR in C. gloeosporioides, was found to be matched with the ranking order given by the NormFinder method (Table 5). Hence, all statistical programs’ overall results determined the CgCAL gene, a highly stable reference gene, followed by the CgTUB2 reference gene, for transcript normalization in the qRT−PCR analysis in C. gloeosporioides for germinating conidia and mycelial stage. However, the CgTUB2 gene was determined to be a highly stable reference, followed by the CgCAL reference gene for expression analysis during its interaction with the plant.
Table 5 Expression stability ranking of candidate reference genes expressed (a) in germinating conidia and mycelium of Colletotrichum gloeosporioides and (b) in infected leaves.
Expression profiling of pathogenicity-related genes for Colletotrichum gloeosporioides for validating selected reference genes
For selected reference gene validation for normalizing the gene expression data in C. gloeosporioides, two pathogenesis-related genes, Diacylglycerol acyltransferase (DGAT1) and qutC, were selected. The expression profiling of both genes was performed using qRT-PCR analysis by adopting the earlier profiles. The expression of these genes was determined relative to selected reference genes (CgCAL and CgTUB2) for normalization. DGAT1 gene expression was higher in germinating conidia, but its conspicuous downregulation was observed in mycelium. However, qutC gene expression was higher in the infected plant sample (4th-week postinoculation), followed by mycelium, which supports the role of this gene in infection (Figure 2). Expression profiling of PR genes for validating the selected reference genes for accurate transcript normalization by determining the relative transcript level of PR genes revealed stable and constant relative expression pattern of selected reference genes for both PR genes.
Figure 2 Expression profiling of DGAT1 and qutC target genes in Colletotrichum gloeosporioides for germinating conidia, mycelium, and interaction with guava plants. The selected reference genes (CgCAL and CgTUB2) were used for data normalization: (A) CgCAL, (B) CgTUB2.
Discussion
Quantitative real-time PCR has developed an imperative technique for studying transcript profiling (Wang et al., 2017). However, the reliability of the results depends upon the reliable reference genes vis-à-vis their expression stability for normalizing expression data (Huang et al., 2018) that should not be altered under varied experimental conditions (Melgar-Rojas et al., 2015; Ijaz et al., 2020b). Under varying experimental conditions, it is ambiguous to use traditional reference genes for normalizing target gene expression without prior investigation of their stability (Galli et al., 2015; Nasir et al., 2023).
There is no report on selecting and validating reliable reference genes for qRT−PCR analysis of target gene expression in C. gloeosporioides and during its interaction with host plants. The specificity of endogenous reference genes is a prerequisite for studying plant-pathogen interactions. Hence, in addition to analyzing candidate reference gene primer specificity for C. gloeosporioides cDNA amplification, the cDNA of the noninfected guava plant was also subjected to PCR analysis using these primers. The 18S rRNA gene was eliminated for qRT−PCR to study C. gloeosporioides gene expression and validate the best reference gene for gene expression analysis on two critical developmental stages, germinating conidia and mycelium, and during its interaction with the guava plant because it gave amplification in the noninfected leaves of guava plants. Statistical software for determining reliable reference genes in fungi and experimental samples influences the results (Song et al., 2019). Therefore, we chose C. gloeosporioides samples and leaf-inoculated guava plants at different stages to detect and validate the reference gene for qRT−PCR.
For qRT-PCR analysis, the candidate genes were evaluated by Cq values. We analyzed the expression stability of each candidate reference gene in germinating conidia and mycelium and three- and four-week-infected guava leaves. The candidate reference gene CgACT was unsuitable for gene expression stability analyses because it showed a higher expression level in infected leaves than in germinating conidia and mycelium samples. In infected plants, fungal pathogen inoculation varies considerably during infection (Vieira et al., 2016). The gene expression of C. gloeosporioides during the germinating conidia and mycelium and in postinoculated leaves of guava plants showed increased expression levels of the CgTUB2 gene.
Several statistical approaches select reference genes with stable expression across biological samples. We used geNorm, NormFinder, and BestKeeper for stability analysis. However, all these algorithms gave varied results, so the RefFinder method was used for comprehensive analysis to avoid the contradiction of individual methods. Based on comprehensive ranking using Ref Finder, CgCAL was selected as the top-ranked, most reliable candidate reference gene for normalizing qRT−PCR of target gene expression in C. gloeosporioides, followed by the CgTUB2 reference gene. However, CgTUB2 was selected as the top-ranked, most reliable candidate reference gene for qRT−PCR of C. gloeosporioides’ target gene expression normalization in infected leaves, followed by the CgCAL reference gene. Our results supported the findings of Vieira et al. (2016). They validated the β-tubulin gene as a stable reference gene for expression studies in C. kahawae. However, Paolinelli-Alfonso et al. (2016) recommended β-tubulin, the best reference gene for normalization in Lasiodiplodia theobromae. Therefore, this study is the first to evaluate the selection and validation of reliable reference genes for normalization in C. gloeosporioides expression analysis.
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 in the article/supplementary material.
Author contributions
Conceptualization, IH and SI. Data curation, IH and SI. Formal analysis, IH and SI. Funding acquisition, EA_A. Writing – original draft, IH and SI. Writing – review & editing, IH, SI, EA_A, AH, and GA-Q. All authors contributed to the article and approved the submitted version.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The authors would like to extend their sincere appreciation to the Researchers Supporting Project Number (RSP2023R134), King Saud University, Riyadh, Saudi Arabia.
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.
References
Babar, M., Ijaz, S., Khan, M. S., Haq, I. U. (2021). Computational genomics based probing of resistance gene analogs (RGAs) in mungbean under cercospora leaf spot disease challenge. Pakistan J. Agric. Sci. 58. doi: 10.3389/fgene.2022.1036029
Chen, A., Liu, T., Wang, Z., Chen, X. (2022). Plant root suberin: A layer of defence against biotic and abiotic stresses. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.1056008
Chen, X. Z., Yang, Q., Wan, Z., Zhou, G. Y., Liu, J. A. (2021). First report of Colletotrichum brevisporum causing anthracnose of Dalbergia odorifera in China. Plant Dis. 105, 2255. doi: 10.1094/PDIS-09-20-1937-PDN
Galli, V., Borowski, J. M., Perin, E. C., da Silva Messias, R., Labonde, J., dos Santos Pereira, I., et al. (2015). Validation of reference genes for accurate normalization of gene expression for real time-quantitative PCR in strawberry fruits using different cultivars and osmotic stresses. Gene. 554, 205–214. doi: 10.1016/j.gene.2014.10.049
Guenin, S., Mauriat, M., Pelloux, J., Van Wuytswinkel, O., Bellini, C., Gutierrez, L. (2009). Normalization of qRT-PCR data: the necessity of adopting a systematic, experimental conditions-specific, validation of references. J. Exp. Bot. 60, 487–493. doi: 10.1093/jxb/ern305
Hao, X., Horvath, D. P., Chao, W. S., Yang, Y., Wang, X., Xiao, B. (2014). Identification and evaluation of reliable reference genes for quantitative real-time PCR analysis in tea plant (Camellia sinensis (L.) O. Kuntze). Int. J. Mol. Sci. 15, 22155–22172. doi: 10.3390/ijms151222155
Haq, I., Ijaz, S. (2019). Assessment of genetic diversity based on ISSR markers in Neopestalotiopsis species collected from guava (Psidium guajava L.) plants affected with canker disease in Pakistan. Appl. Ecol. Environ. Res. 17, 11803–11811. doi: 10.15666/aeer/1705_1180311811
Haq, I., Ijaz, S. (2020). History and recent trends in plant disease control: An overview. Plant Dis. Manage. Strateg. Sustain. Agric. through Tradit. Mod. approaches. 1-13. doi: 10.1007/978-3-030-35955-3_1
Haq, I., Ijaz, S., Khan, N. A. (2021). Genealogical concordance of phylogenetic species recognition-based delimitation of Neopestalotiopsis species associated with leaf spots and fruit canker disease affected guava plants. Pakistan J. Agric. Sci. 58, 1301–1313. doi: 10.21162/PAKJAS/21.1045
Haq, I., Ijaz, S., Latif, M. Z. (2019). Multilocus sequence typing (MLST) based genetic variation analysis of shisham dieback associated strains of Ceratocystis fimbriata sensu lato species complex in Pakistan. Appl. Ecol. Environ. Res. 17, 12573–12582. doi: 10.15666/aeer/1705_1257312582
Huang, N., Ling, H., Liu, F., Su, Y., Su, W., Mao, H., et al. (2018). Identification and evaluation of PCR reference genes for host and pathogen in sugarcane-Sporisorium scitamineum interaction system. BMC Genomics 19, 1–13. doi: 10.1186/s12864-018-4854-z
Huggett, J., Dheda, K., Bustin, S., Zumla, A. (2005). Real-time RT-PCR normalisation; strategies and considerations. Genes Immun. 6, 279–284. doi: 10.1038/sj.gene.6364190
Ijaz, S., Babar, M., Razzaq, H. A., Nasir, B. (2020a). Transgenic approaches in plants: Strategic control for disease management. Plant Dis. Manage. Strateg. Sustain. Agric. through Tradit. Mod. Approaches., 187–215. doi: 10.1007/978-3-030-35955-3_9
Ijaz, S., Haq, I. U. (2020). Genome editing technologies for resistance against phytopathogens: Principles, applications and future prospects. Plant Dis. Manage. Strateg. Sustain. Agric. through Tradit. Mod. Approaches, 237–245. doi: 10.1007/978-3-030-35955-3_11
Ijaz, S., Haq, I. U., Nasir, B. (2020b). In silico identification of expressed sequence tags based simple sequence repeats (EST-SSRs) markers in Trifolium species. ScienceAsia. 46, 6–10. doi: 10.2306/scienceasia1513-1874.2020.001
Ijaz, S., Haq, I. U., Razzaq, H. A., Nasir, B., Babar, M. (2019). ISSR-based population genetics study for tagging a diverse population of shisham (Dalbergia sissoo) in Pakistan. Appl. Ecol. Environ. Res. 17, 5851–5861. doi: 10.15666/aeer/1703_58515861
Iqra, Haq, I. U., Waseem Khan Qadri, R., Amrao, L., Ijaz, S. (2022). Effect of environmental conditions (temperature and precipitation) on severity of guava die-back caused by Colletotrichum spp. under climatic conditions of Pakistan. J. Plant Pathol. 45, 1–12. doi: 10.1007/s42161-021-00968-1
Liang, X., Shang, S., Dong, Q., Wang, B., Zhang, R., Gleason, M. L., et al. (2018). Transcriptomic analysis reveals candidate genes regulating development and host interactions of Colletotrichum fructicola. BMC Genomics 19, 1–21. doi: 10.1007/s42161-021-00968-1
Liu, L., Yan, Y., Huang, J., Hsiang, T., Wei, Y., Li, Y., et al. (2017). A novel MFS transporter gene ChMfs1 is important for hyphal morphology, conidiation, and pathogenicity in Colletotrichum higginsianum. Front. Microbiol. 8. doi: 10.3389/fmicb.2017.01953
Lu, Q., Wang, Y., Li, N., Ni, D., Yang, Y., Wang, X. (2018). Differences in the characteristics and pathogenicity of Colletotrichum camelliae and C. fructicola isolated from the tea plant [Camellia sinensis (L.) O. Kuntze]. Front. Microbiol. 9. doi: 10.3389/fmicb.2018.03060
MaChado, R. D., Christoff, A. P., Loss-Morais, G., Margis-Pinheiro, M., Margis, R., Körbes, A. P. (2015). Comprehensive selection of reference genes for quantitative gene expression analysis during seed development in Brassica napus. Plant Cell Rep. 34, 1139–1149. doi: 10.1007/s00299-015-1773-1
Marcial-Quino, J., Fierro, F., Enríquez-Flores, S., Gómez-Manzo, S., Vanoye-Carlo, A., Garcia-Torres, I., et al. (2016). Validation of housekeeping genes as an internal control for gene expression studies in Giardia lamblia using quantitative real-time PCR. Gene 581, 21–30. doi: 10.1016/j.gene.2016.01.018
Melgar-Rojas, P., Alvarado, J. C., Fuentes-Santamaría, V., Gabaldón-Ull, M. C., Juiz, J. M. (2015). Validation of reference genes for RT-qPCR analysis in noise-induced hearing loss: a study in Wistar rat. PloS One 10, e0138027. doi: 10.1371/journal.pone.0138027
Nasir, B., Ijaz, S., Ul Haq, I., Khan, N. A., Habib, A. (2023). Identification, in silico characterization, and expression analysis of NBS-LRR class of R genes against stem and crown rot disease in. Trifolium alexandrinum L. Sci. Asia. 49. doi: 10.3389/fgene.2022.1036029
Paolinelli-Alfonso, M., Galindo-Sánchez, C. E., Hernandez-Martinez, R. (2016). Quantitative real-time PCR normalization for gene expression studies in the plant pathogenic fungi. Lasiodiplodia theobromae. J. Microbiol. Methods 127, 82–88. doi: 10.1016/j.mimet.2016.05.021
Radonić, A., Thulke, S., Mackay, I. M., Landt, O., Siegert, W., Nitsche, A. (2004). Guideline to reference gene selection for quantitative real-time PCR. Biochem. Biophys. Res. Commun. 313, 856–862. doi: 10.1016/j.bbrc.2003.11.177
Song, Y., Wang, Y., Guo, D., Jing, L. (2019). Selection of reference genes for quantitative real-time PCR normalization in the plant pathogen Puccinia helianthi Schw. BMC Plant Biol. 19, 1–10. doi: 10.1186/s12870-019-1629-x
Vieira, A., Cabral, A., Fino, J., Azinheira, H. G., Loureiro, A., Talhinhas, P., et al. (2016). Comparative validation of conventional and RNA-Seq data-derived reference genes for qPCR expression studies of Colletotrichum kahawae. PloS One 11, e0150651. doi: 10.1371/journal.pone.0150651
Villa-Rivera, M. G., Conejo-Saucedo, U., Lara-Marquez, A., Cano-Camacho, H., Lopez-Romero, E., Zavala-Paramo, M. G. (2017). The role of virulence factors in the pathogenicity of Colletotrichum sp. Curr. Protein Pept. Sci. 18, 1005–1018. doi: 10.2174/1389203717666160813160727
Wang, L., Wang, Y., Cao, H., Hao, X., Zeng, J., Yang, Y., et al. (2016). Transcriptome analysis of an anthracnose-resistant tea plant cultivar reveals genes associated with resistance to Colletotrichum camelliae. PloS One 11, e0148535. doi: 10.1371/journal.pone.0148535
Wang, M.-L., Li, Q.-H., Xin, H.-H., Chen, X., Zhu, X.-J., Li, X.-H. (2017). Reliable reference genes for normalization of gene expression data in tea plants (Camellia sinensis) exposed to metal stresses. PloS One 12, e0175863. doi: 10.1371/journal.pone.0175863
Keywords: qRT−PCR, reference genes, Colletotrichum gloeosporioides, Psidium guajava, pathogenesis-related genes
Citation: Haq IU, Ijaz S, Hashem A, Avila-Quezada GD and Abd_Allah EF (2023) Selection and validation of reference genes for normalizing qRT‒PCR gene expression studies in Colletotrichum gloeosporioides and interaction with the guava plants. Front. Plant Sci. 14:1235848. doi: 10.3389/fpls.2023.1235848
Received: 06 June 2023; Accepted: 07 November 2023;
Published: 23 November 2023.
Edited by:
Heba I. Mohamed, Ain Shams University, EgyptReviewed by:
Rashda Naheed, Government of Punjab, PakistanJun Yang, Chuxiong Normal University, China
Ani Widiastuti, Gadjah Mada University, Indonesia
Copyright © 2023 Haq, Ijaz, Hashem, Avila-Quezada and Abd_Allah. 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: Siddra Ijaz, c2lkZHJhaWphemtoYW5AeWFob28uY29t; ZHJzaWRkcmFpamF6QHVhZi5lZHUucGs=