Vieira Ana, Cabral Ana, Fino Joana, Azinheira Helena G, Loureiro Andreia, Talhinhas Pedro, Pires Ana Sofia, Varzea Vitor, Moncada Pilar, Oliveira Helena, Silva Maria do Céu, Paulo Octávio S, Batista Dora
CIFC-Centro de Investigação das Ferrugens do Cafeeiro, Instituto Superior de Agronomia, Universidade de Lisboa, Oeiras, Portugal.
Computational Biology and Population Genomics group, cE3c - Centre for Ecology, Evolution and Environmental Changes, Faculdade de Ciências, Universidade de Lisboa, Lisboa, Portugal.
PLoS One. 2016 Mar 7;11(3):e0150651. doi: 10.1371/journal.pone.0150651. eCollection 2016.
Colletotrichum kahawae is an emergent fungal pathogen causing severe epidemics of Coffee Berry Disease on Arabica coffee crops in Africa. Currently, the molecular mechanisms underlying the Coffea arabica-C. kahawae interaction are still poorly understood, as well as the differences in pathogen aggressiveness, which makes the development of functional studies for this pathosystem a crucial step. Quantitative real time PCR (qPCR) has been one of the most promising approaches to perform gene expression analyses. However, proper data normalization with suitable reference genes is an absolute requirement. In this study, a set of 8 candidate reference genes were selected based on two different approaches (literature and Illumina RNA-seq datasets) to assess the best normalization factor for qPCR expression analysis of C. kahawae samples. The gene expression stability of candidate reference genes was evaluated for four isolates of C. kahawae bearing different aggressiveness patterns (Ang29, Ang67, Zim12 and Que2), at different stages of fungal development and key time points of the plant-fungus interaction process. Gene expression stability was assessed using the pairwise method incorporated in geNorm and the model-based method used by NormFinder software. For C. arabica-C. kahawae interaction samples, the best normalization factor included the combination of PP1, Act and ck34620 genes, while for C. kahawae samples the combination of PP1, Act and ck20430 revealed to be the most appropriate choice. These results suggest that RNA-seq analyses can provide alternative sources of reference genes in addition to classical reference genes. The analysis of expression profiles of bifunctional catalase-peroxidase (cat2) and trihydroxynaphthalene reductase (thr1) genes further enabled the validation of the selected reference genes. This study provides, for the first time, the tools required to conduct accurate qPCR studies in C. kahawae considering its aggressiveness pattern, developmental stage and host interaction.
咖啡炭疽病菌是一种新出现的真菌病原体,在非洲的阿拉比卡咖啡作物上引发咖啡浆果病的严重流行。目前,阿拉比卡咖啡与咖啡炭疽病菌相互作用的分子机制以及病原体侵袭性的差异仍知之甚少,这使得开展该病理系统的功能研究成为关键一步。定量实时PCR(qPCR)一直是进行基因表达分析最有前景的方法之一。然而,使用合适的内参基因进行恰当的数据标准化是绝对必要的。在本研究中,基于两种不同方法(文献和Illumina RNA测序数据集)选择了一组8个候选内参基因,以评估用于咖啡炭疽病菌样本qPCR表达分析的最佳标准化因子。针对具有不同侵袭性模式的4株咖啡炭疽病菌(Ang29、Ang67、Zim12和Que2),在真菌发育的不同阶段以及植物 - 真菌相互作用过程的关键时间点,评估了候选内参基因的基因表达稳定性。使用geNorm软件中包含的成对法和NormFinder软件使用的基于模型的方法评估基因表达稳定性。对于阿拉比卡咖啡 - 咖啡炭疽病菌相互作用样本,最佳标准化因子包括PP1、Act和ck34620基因的组合,而对于咖啡炭疽病菌样本,PP1、Act和ck20430的组合被证明是最合适的选择。这些结果表明,除了经典内参基因外,RNA测序分析还可以提供内参基因的替代来源。对双功能过氧化氢酶 - 过氧化物酶(cat2)和三羟基萘还原酶(thr1)基因表达谱的分析进一步验证了所选的内参基因。本研究首次提供了在考虑咖啡炭疽病菌的侵袭性模式、发育阶段和宿主相互作用的情况下进行准确qPCR研究所需的工具。