
冬虫夏草菌实时荧光定量PCR内参基因的筛选
苏强军,夏樱霞,谢放,Uwitugabiye VESTINE,陈照禾,周刚
菌物学报 ›› 2021, Vol. 40 ›› Issue (7) : 1712-1722.
冬虫夏草菌实时荧光定量PCR内参基因的筛选
Screening of the reference genes for qRT-PCR analysis of gene expression in Ophiocordyceps sinensis
以冬虫夏草单子囊孢子分离得到的菌株TZ8-1的3种菌丝形态为实验材料,提取RNA,经反转录获取cDNA,选择了11个持家基因为候选内参基因(18S rRNA、APRTase、β-TUB、RPL2、EF1-α、PGI、PGM、H+-ATPase、ACT1、UBQ和GAPDH),根据该菌基因组注释结果来设计引物,采用实时荧光定量PCR(qRT-PCR)技术进行定量扩增,利用geNorm、NormFinder和BestKeeper算法程序进行表达稳定性评估,并用RefFinder对评估结果进行综合排比,最终筛选得到了最适内参基因。结果表明,所选取的11个候选内参基因均可作为冬虫夏草菌菌丝体时期的内参基因,稳定性最好的3个内参基因分别是UBQ、PGE和ACT1,稳定性较差的3个内参基因分别是GAPDH、RPL2和β-TUB。
Mycelia with three different morphology of Ophiocordyceps sinensis TZ8-1 isolated from monoascospores were used as experimental materials. RNA of mycelium was extracted, and cDNA was obtained by reverse transcription. Eleven housekeeping genes were selected as candidate reference genes (18S rRNA, APRTase, β-TUB, RPL2, EF1-α, PGI, PGM, H+-ATPase, ACT1, UBQ and GAPDH). Primers were designed based on the genome annotation results of Ophiocordyceps sinensis and quantitative amplification was conducted by real-time fluorescence quantitative PCR (qRT-PCR). The expression of stability was evaluated by geNorm, NormFinder and BestKeeper algorithms, and the evaluation of results were combined by RefFinder to screen the optimal internal reference genes. The results showed that all the 11 candidate reference genes could be used as reference genes in the mycelium stage of Ophiocordyceps sinensis. The three most stable ones were UBQ, PGE and ACT1, and the three less stable ones were GAPDH, RPL2 and β-TUB.
内参基因 / 冬虫夏草菌 / qRT-PCR / 菌丝形态 {{custom_keyword}} /
reference genes / Ophiocordyceps sinensis / real-time fluorescence quantitative-PCR / mycelium morphology {{custom_keyword}} /
表1 候选内参基因的引物序列Table 1 Primer sequences of candidate internal reference genes |
基因名称 Gene name | 正向引物 Forward primer (5’-3’) | 反向引物 Reverse primer (5’-3’) | 扩增产物 Amplicon size (bp) | Tm (°C) |
---|---|---|---|---|
18S rRNA | GCAGTGGCATCTCTCAGTC | TCATCGATGCCAGAACC | 128 | 60 |
APRTase | ATGCTGAGCTGTTTAGCG | TGCCCTGTTCGTCGTAGA | 90 | 60 |
β-TUB | TACGCCTCTTCGACGATAG | GCCGTTGTAGACACCATT | 148 | 60 |
RPL2 | ACCTACCGTCTCCATCAT | GTGAACCTGCTGGACAAT | 109 | 60 |
EF1-α | CAAGGGCTCTTTCAAGTATGC | GTGACATAGTACCTGGGAGT | 114 | 60 |
PGI | TTCGACCAGTATCTTCATCGC | GTGTACTTGACCGACGATCC | 83 | 60 |
PGM | AAGCCCTTTCAGGACCAA | GAACGACTCGGTGTAGTG | 87 | 60 |
H+-ATPase | CGCTTCGCGGAAATCTATAC | AGACGTTCTTGTTGACGG | 90 | 60 |
ACT1 | CAATCGGCACAACTGGACA | GACGACCTGAGCGGAATA | 95 | 60 |
UBQ | CGACATCGAGTTGGACTAC | ATACCTGCAATCTGTCCG | 82 | 60 |
GAPDH | GAGGCCGAGAGCCAACTA | TTCATCACGACAGCACCA | 99 | 60 |
表2 RNA质检结果Table 2 Results of RNA quality control |
样本 Sample | 浓度 Concentration (μg/μL) | A260/280 | A260/230 | 体积 Volume (μL) | 总量 Total (μg) | 28S/18S | RIN |
---|---|---|---|---|---|---|---|
A1 | 0.3864 | 2.18 | 1.84 | 20 | 7.73 | 1.7 | 9.5 |
A2 | 0.1979 | 2.19 | 2.28 | 20 | 3.96 | 0.9 | 7.7 |
A3 | 0.1701 | 2.14 | 1.78 | 20 | 3.40 | 1.0 | 7.2 |
B1 | 0.3490 | 2.18 | 1.81 | 20 | 6.98 | 1.7 | 10.0 |
B2 | 0.2323 | 2.17 | 1.77 | 20 | 4.65 | 2.1 | 9.2 |
B3 | 0.1352 | 2.16 | 1.75 | 15 | 2.03 | 1.2 | 7.1 |
C1 | 0.3873 | 2.16 | 2.24 | 20 | 7.75 | 1.5 | 8.8 |
C2 | 0.5354 | 2.23 | 2.00 | 20 | 10.71 | 1.1 | 7.9 |
C3 | 0.2250 | 2.16 | 2.10 | 20 | 4.50 | 1.3 | 7.9 |
表3 NormFinder分析候选内参基因的表达稳定值Table 3 NormFinder analysis of the stable expression values of the candidate reference genes |
基因名称 Gene name | 稳定值 Stability value | 最优基因 Best gene |
---|---|---|
UBQ | 0.096 | UBQ |
PGM | 0.102 | |
ACT1 | 0.122 | |
H+-ATPase | 0.140 | |
EF-1α | 0.155 | |
18S rRNA | 0.223 | |
APRTase | 0.269 | |
PGI | 0.292 | |
β-TUB | 0.350 | |
RPL2 | 0.363 | |
GAPDH | 0.646 |
表4 BestKeeper程序分析候选内参基因的稳定性表达Table 4 BestKeeper program analyzes the stability expression of candidate reference genes |
基因名称 Gene name | 变异系数 CV | 标准差 SD |
---|---|---|
APRTase | 1.23 | 0.33 |
ACT1 | 1.68 | 0.43 |
H+-ATPase | 1.8 | 0.53 |
UBQ | 2.17 | 0.65 |
PGM | 2.97 | 0.69 |
18S rRNA | 4.31 | 0.7 |
GAPDH | 2.37 | 0.7 |
EF-1A | 3.9 | 0.73 |
PGI | 3.65 | 0.86 |
β-TUB | 3.9 | 0.88 |
RPL2 | 3.86 | 0.9 |
表5 RefFinder综合排名结果Table 5 Comprehensive ranking results of RefFinder |
基因名称 Gene name | Delta CT排名 Delta CT ranking | geNorm排名 geNorm ranking | NormFinder排名 NormFinder ranking | BestKeeper排名 BestKeeper ranking | 几何均值 Geomean | 综合排名 Comprehensive ranking |
---|---|---|---|---|---|---|
UBQ | 2 | 2 | 1 | 4 | 1.78 | 1 |
PGM | 1 | 1 | 2 | 5 | 1.86 | 2 |
ACT1 | 4 | 8 | 3 | 2 | 4.12 | 3 |
EF-1α | 3 | 3 | 5 | 8 | 4.49 | 4 |
H+-ATPase | 5 | 9 | 4 | 3 | 5.18 | 5 |
APRTase | 8 | 10 | 7 | 1 | 5.79 | 6 |
18S rRNA | 6 | 7 | 6 | 6 | 6.48 | 7 |
PGI | 7 | 6 | 8 | 9 | 7.16 | 8 |
β-TUB | 9 | 5 | 9 | 10 | 8.17 | 9 |
RPL2 | 10 | 4 | 10 | 11 | 8.32 | 10 |
GAPDH | 11 | 11 | 11 | 7 | 10.84 | 11 |
[1] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[2] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[3] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[4] |
Pleurotus ostreatus is an important edible mushroom and a model lignin degrading organism, whose genome contains nine genes of ligninolytic peroxidases, characteristic of white-rot fungi. These genes encode six manganese peroxidase (MnP) and three versatile peroxidase (VP) isoenzymes. Using liquid chromatography coupled to tandem mass spectrometry, secretion of four of these peroxidase isoenzymes (VP1, VP2, MnP2 and MnP6) was confirmed when P. ostreatus grows in a lignocellulose medium at 25°C (three more isoenzymes were identified by only one unique peptide). Then, the effect of environmental parameters on the expression of the above nine genes was studied by reverse transcription-quantitative PCR by changing the incubation temperature and medium pH of P. ostreatus cultures pre-grown under the above conditions (using specific primers and two reference genes for result normalization). The cultures maintained at 25°C (without pH adjustment) provided the highest levels of peroxidase transcripts and the highest total activity on Mn(2+) (a substrate of both MnP and VP) and Reactive Black 5 (a VP specific substrate). The global analysis of the expression patterns divides peroxidase genes into three main groups according to the level of expression at optimal conditions (vp1/mnp3>vp2/vp3/mnp1/mnp2/mnp6>mnp4/mnp5). Decreasing or increasing the incubation temperature (to 10°C or 37°C) and adjusting the culture pH to acidic or alkaline conditions (pH 3 and 8) generally led to downregulation of most of the peroxidase genes (and decrease of the enzymatic activity), as shown when the transcription levels were referred to those found in the cultures maintained at the initial conditions. Temperature modification produced less dramatic effects than pH modification, with most genes being downregulated during the whole 10°C treatment, while many of them were alternatively upregulated (often 6h after the thermal shock) and downregulated (12h) at 37°C. Interestingly, mnp4 and mnp5 were the only peroxidase genes upregulated under alkaline pH conditions. The differences in the transcription levels of the peroxidase genes when the culture temperature and pH parameters were changed suggest an adaptive expression according to environmental conditions. Finally, the intracellular proteome was analyzed, under the same conditions used in the secretomic analysis, and the protein product of the highly-transcribed gene mnp3 was detected. Therefore, it was concluded that the absence of MnP3 from the secretome of the P. ostreatus lignocellulose cultures was related to impaired secretion. Copyright © 2014 The Authors. Published by Elsevier Inc. All rights reserved.
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[5] |
Hypoxic culturing of chondrocytes is gaining increasing interest in cartilage research. Culturing of chondrocytes under low oxygen tension has shown several advantages, among them increased synthesis of extracellular matrix and increased redifferentiation of dedifferentiated chondrocytes. Quantitative gene expression analyses such as quantitative real-time PCR (qRT-PCR) are powerful tools in the investigation of underlying mechanisms of cell behavior and are used routinely for differentiation and phenotype assays. However, the genes used for normalization in normoxic cell-cultures might not be suitable in the hypoxic environment. The objective of this study was to determine hypoxia-stable housekeeping genes (HKG) for quantitative real-time PCR (qRT-PCR) in human chondrocytes cultured in 21%, 5% and 1% oxygen by geNorm and NormFinder analyses.The chondrocytic response to the hypoxic challange was validated by a significant increase in expression of the hypoxia-inducible gene ankyrin repeat 37 as well as SOX9 in hypoxia. When cultured on the 3-dimentional (3D) scaffold TATA-binding protein (TBP) exhibited the highest expression stability with NormFinder while Ribosomal protein L13a (RPL13A) and beta2-microglobulin (B2M) were the most stable using geNorm analysis. In monolayer RPL13A were the most stable gene using NormFinder, while geNorm assessed RPL13A and human RNA polymerase II (RPII) as most stable. When examining the combination of (3D) culturing and monolayer RPL13A and B2M showed the highest expression stability from geNorm analysis while RPL13A also showed the highest expression stability using NormFinder. Often used HKG such as beta actin (ACTB) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) were the most unstable genes investigated in all comparisons. The pairwise variations for the two most stable HKG in each group were all below the cut-off value of 0.15, suggesting that the two most stable HKG from geNorm analysis would be sufficient for qRT-PCR.All data combined we recommend RPL13A, B2M and RPII as the best choice for qRT-PCR analyses when comparing normoxic and hypoxic cultured human chondrocytes although other genes might also be suitable. However, the matching of HKG to target genes by means of a thorough investigation of the stability in each study would always be preferable.
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[6] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[7] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[8] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[9] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[10] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[11] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[12] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[13] |
We examined eight putative consistently expressed genes-actin (ACT), beta-tubulin, elongation factor 1alpha (EF1), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), phosphoglycerate kinase (PGK), ribosomal protein L2 (RPL2), ubiquitin (UBI), and a catalytic subunit of protein phosphatase 2A (PP2Acs)-for their potential as references for the normalization of gene expression in tomato leaves. Expression stability of candidate reference genes was tested during growth conditions of nitrogen (N) starvation, low temperature, and suboptimal light. The geNorm algorithm, using reciprocal cross-validation among a larger group of candidate references, was applied for this purpose. The widely used reference genes GAPDH and PGK were top ranked during light stress but poorly ranked during N and cold stress. In contrast, EF1 was top ranked during N and cold stress but poorly ranked during light stress. The novel references RPL2 and PP2Acs, as well as the traditional references ACT and UBI, appeared to be stably expressed when looking at the data set as a whole. No gene was identified that exhibited such a constant level of expression as to outperform the other candidates under all experimental conditions. Thus, the results highlight the need for normalizing gene expression in tomato using the geometric average of multiple carefully selected reference genes.
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[14] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[15] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[16] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[17] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[18] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[19] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[20] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[21] |
Lycoris aurea (L' Hér.) Herb, a perennial grass species, produces a unique variety of pharmacologically active Amaryllidaceae alkaloids. However, the key enzymes and their expression pattern involved in the biosynthesis of Amaryllidaceae alkaloids (especially for galanthamine) are far from being fully understood. Quantitative real-time polymerase chain reaction (qRT-PCR), a commonly used method for quantifying gene expression, requires stable reference genes to normalize its data. In this study, to choose the appropriate reference genes under different experimental conditions, 14 genes including YLS8 (mitosis protein YLS8), CYP2 (Cyclophilin 2), CYP 1 (Cyclophilin 1), TIP41 (TIP41-like protein), EXP2 (Expressed protein 2), PTBP1 (Polypyrimidine tract-binding protein 1), EXP1 (Expressed protein 1), PP2A (Serine/threonine-protein phosphatase 2A), β-TUB (β-tubulin), α-TUB (α-tubulin), EF1-α (Elongation factor 1-α), UBC (Ubiquitin-conjugating enzyme), ACT (Actin) and GAPDH (Glyceraldehyde 3-phosphate dehydrogenase) were selected from the transcriptome datasets of L. aurea. And then, expressions of these genes were assessed by qRT-PCR in various tissues and the roots under different treatments. The expression stability of the 14 candidates was analyzed by three commonly used software programs (geNorm, NormFinder, and BestKeeper), and their results were further integrated into a comprehensive ranking based on the geometric mean. The results show the relatively stable genes for each subset as follows: (1) EXP1 and TIP41 for all samples; (2) UBC and EXP1 for NaCl stress; (3) PTBP1 and EXP1 for heat stress, polyethylene glycol (PEG) stress and ABA treatment; (4) UBC and CYP2 for cold stress; (5) PTBP1 and PP2A for sodium nitroprusside (SNP) treatment; (6) CYP1 and TIP41 for methyl jasmonate (MeJA) treatment; and (7) EXP1 and TIP41 for various tissues. The reliability of these results was further enhanced through comparison between part qRT-PCR result and RNA sequencing (RNA-seq) data. In summary, our results identified appropriate reference genes for qRT-PCR in L. aurea, and will facilitate gene expression studies under these conditions.
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[22] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[23] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[24] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[25] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[26] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[27] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[28] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[29] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[30] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[31] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[32] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[33] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[34] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[35] |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[36] |
安丹丹, 李健, 隋飞飞, 曲柏宏, 陈艳秋, 张宇婷, 2016. 桦褐孔菌实时荧光定量PCR内参基因的筛选. 食用菌学报, 23(3):15-24
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[37] |
李增智, 黄勃, 李春如, 樊美珍, 2000. 确证冬虫夏草无性型的分子生物学证据:Ⅰ. 中国被毛孢与冬虫夏草之间的关系. 菌物学报, 19(1):60-64
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[38] |
刘锡琎, 郭英兰, 俞永信, 曾纬, 1989. 冬虫夏草菌无性阶段的分离和鉴定. 真菌学报, 8(1):35-39
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[39] |
刘作易, 梁宗琦, 刘爱英, 2003. 冬虫夏草子囊孢子萌发及其无性型观察. 贵州农业科学, 2003(1):3-5
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[40] |
莫明和, 迟胜起, 张克勤, 2001. 冬虫夏草的微循环产孢及其无性型的分离. 菌物学报, 20(4):482-485
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[41] |
秦晓艺, 王杰, 2015. 杏鲍菇采后木质化相关基因实时定量PCR中内参基因的选择. 西北农林科技大学学报(自然科学版), 43(7):219-227
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[42] |
吴建阳, 何冰, 杜玉洁, 李伟才, 魏永赞, 2017. 利用geNorm、NormFinder和BestKeeper软件进行内参基因稳定性分析的方法. 现代农业科技, 2017(5):278-281
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[43] |
谢放, 朱子雄, 陈京津, 2011. 甘肃省冬虫夏草菌的RAPD和ITS序列分析. 中国微生态学杂志, 23(3):219-22
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[44] |
易乐飞, 穆亮亮, 仲改, 2016. 凡纳滨对虾总RNA完整性的正确评估. 浙江农业学报, 28(12):2007-2013
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[45] |
张倩倩, 刘伟, 边银丙, 2018. 羊肚菌实时荧光定量PCR内参基因的筛选. 中国菌物学会学术年会论文汇编,北京. 78
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[46] |
张越, 姚方杰, 孙文娟, 方明, 武春爽, 2020. 黑木耳实时荧光定量PCR内参基因的筛选. 菌物学报, 39(8):1510-1519
{{custom_citation.content}}
{{custom_citation.annotation}}
|
[47] |
赵小龙, 2016. 茯苓qRT-PCR内参基因筛选的研究. 武汉轻工大学硕士论文,武汉. 32-35
{{custom_citation.content}}
{{custom_citation.annotation}}
|
{{custom_ref.label}} |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
/
〈 |
|
〉 |