Chinese hamster ovary cells are the predominant cell lines used for bio-therapeutic production. Real-time quantitative PCR (RT-qPCR) and transcriptomics are powerful tools to understand and optimize the Chinese hamster ovary cells for higher productivity or better control of product qualities. Reliable reference genes, which were proved to be experiment-specific, are critical yardsticks. In this study, we compared expression stability of 20 candidate reference genes at mRNA level, including commonly used housekeeping genes, previous literature reported genes in Chinese hamster ovary cells producing an intact antibody, and new candidates suggested by our RNA-seq transcriptomic database, in RT-qPCR reactions in Fc-fusion protein-producing Chinese hamster ovary cells with various productivity during long-term cultivation and fed-batch cultures at 26 different conditions. geNorm, NormFinder, BestKeeper, and ΔCt programs and methods were utilized to analyze the gene expression stability and gave an overall ranking. Akr1a1, Gpx1, and Aprt in long-term cultivation and Akr1a1, Rps16 in fed-batch culture, which have not been reported previously, exhibited the highest stability of gene expression, while Pabpn1, Hirip3, and Actb in both sets of experiments together with Atp5f1 in long-term passage process showed the weakest stability. The results were then validated using GLP1-Fc (Glucagon-like peptide-1 Fc fusion protein) gene as the target with determined expression level which were doubly confirmed by both absolute RT-qPCR and confocal microscopy. These new references should be considered for the investigations on Chinese hamster ovary cells in related research.
Impact statement
In order to reveal potential genotype-phenotype relationship, RT-qPCR reactions are frequently applied which require validated and reliable reference genes. With the investigation on long-term passage and fed-batch cultivation of CHO cells producing an Fc-fusion protein, four new reference genes–Akr1a1, Gpx1, Aprt, and Rps16, were identified from 20 candidates with the aid of geNorm, NormFinder, BestKeeper, and ΔCt programs and methods. This article provided more verified options in reference gene selection in related research on CHO cells.
KaneyoshiKKurodaKUchiyamaKOnitsukaMYamano-AdachiNKogaYOmasaT.Secretion analysis of intracellular “difficult-to-express” immunoglobulin G (IgG) in Chinese hamster ovary (CHO) cells. Cytotechnology2019;
71:305–16
2.
FischerSHandrickROtteK.The art of CHO cell engineering: a comprehensive retrospect and future perspectives. Biotechnol Adv2015;
33:1878–96
3.
StolfaGSmonskeyMTBonifaceRHachmannA-BGuldePJoshiADPierceAP.CHO-Omics review: the impact of current and emerging technologies on Chinese hamster ovary based bioproduction. Biotechnol J2018;
13:17002–27
4.
DattaPLinhardtRJSharfsteinST.An ‘omics approach towards CHO cell engineering. Biotechnol Bioeng2013;
110:1255–71
5.
KimJYKimYGLeeGM.CHO cells in biotechnology for production of recombinant proteins: current state and further potential. Appl Microbiol Biotechnol2012;
93:917–30
6.
KramerOKlausingSNollT.Methods in mammalian cell line engineering: from random mutagenesis to sequence-specific approaches. Appl Microbiol Biotechnol2010;
88:425–36
7.
ShridharSKlanertGAuerNHernandez-LopezIKandulaMMHacklMBorthN.Transcriptomic changes in CHO cells after adaptation to suspension growth in protein-free medium analysed by a species-specific microarray. J Biotechnol2017;
257:13–21
8.
Fomina-YadlinDMujacicMMaggioraKQuesnellGSaleemRMcGrewJT.Transcriptome analysis of a CHO cell line expressing a recombinant therapeutic protein treated with inducers of protein expression. J Biotechnol2015;
212:106–15
9.
PybusLPDeanGWestNRSmithADaramolaOFieldRJamesDC.Model-directed engineering of “difficult-to-express” monoclonal antibody production by Chinese hamster ovary cells. Biotechnol Bioeng2014;
111:372–85
10.
WurmFMHackerD.First CHO genome. Nat Biotechnol2011;
29:718–20
11.
LewisNELiuXLiYNagarajanHYerganianGO’BrienEPalssonBO.Genomic landscapes of Chinese hamster ovary cell lines as revealed by the cricetulus griseus draft genome. Nat Biotechnol2013;
31:759–65
12.
VishwanathanNLeHJacobNMTsaoYSNgSWLooBHuWS.Transcriptome dynamics of transgene amplification in Chinese hamster ovary cells. Biotechnol Bioeng2014;
111:518–28
13.
PaninaYGermondAMasuiSWatanabeTM.Validation of common housekeeping genes as reference for qPCR gene expression analysis during iPS reprogramming process. Sci Rep2018;
8:8716
14.
SuzukiTHigginsPJCrawfordDR.Control selection for RNA quantitation. Biotechniques2000;
29:332–7
15.
DerksNMMullerMGasznerBTilburg-OuwensDTRoubosEWKoziczLT.Housekeeping genes revisited: different expressions depending on gender, brain area and stressor. Neuroscience2008;
156:305–9
16.
GueninSMauriatMPellouxJVan WuytswinkelOBelliniCGutierrezL.Normalization of qRT-PCR data: the necessity of adopting a systematic, experimental conditions-specific, validation of references. J Exp Bot2009;
60:487–93
17.
LeePDSladekRGreenwoodCMHudsonTJ.Control genes and variability: absence of ubiquitous reference transcripts in diverse mammalian expression studies. Genome Res2002;
12:292–7
18.
BrownAJGibsonSHattonDJamesDC.Transcriptome-Based identification of the optimal reference CHO genes for normalisation of qPCR data. Biotechnol J2018;
13:1700259–65
19.
DijkstraJRvan KempenLCNagtegaalIDBustinSA.Critical appraisal of quantitative PCR results in colorectal cancer research: can we rely on published qPCR results?Mol Oncol2014;
8:813–8
20.
BahrSMBorgschulteTKayserKJLinN.Using microarray technology to select housekeeping genes in Chinese hamster ovary cells. Biotechnol Bioeng2009;
104:1041–6
21.
ParedesVParkJSJeongYYoonJBaekK.Unstable expression of recombinant antibody during long-term culture of CHO cells is accompanied by histone H3 hypoacetylation. Biotechnol Lett2013;
35:987–93
22.
KimMO’CallaghanPMDromsKAJamesDC.A mechanistic understanding of production instability in CHO cell lines expressing recombinant monoclonal antibodies. Biotechnol Bioeng2011;
108:2434–46
23.
JiaYLGuoXLuJTWangXYQiuLLWangTY.CRISPR/Cas9-mediated gene knockout for DNA methyltransferase Dnmt3a in CHO cells displays enhanced transgenic expression and long-term stability. J Cell Mol Med2018;
22:4106–16
24.
VandesompeleJDe PreterKPattynFPoppeBVan RoyNDe PaepeASpelemanF.Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol2002;
3:Research0034.1–11
25.
AndersenCLJensenJLOrntoftTF.Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res2004;
64:5245–50
26.
PfafflMWTichopadAPrgometCNeuviansTP.Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper–excel-based tool using pair-wise correlations. Biotechnol Lett2004;
26:509–15
27.
SilverNBestSJiangJTheinS.Selection of housekeeping genes for gene expression studies in human reticulocytes using real-time PCR. BMC Mol Biol2006;
7:33–41
28.
XieFXiaoPChenDXuLZhangB.miRDeepFinder: a miRNA analysis tool for deep sequencing of plant small RNAs. Plant Mol Biol2012;
80:75–84
29.
LangmeadBSalzbergSL.Fast gapped-read alignment with bowtie 2. Nat Methods2012;
9:357–9
30.
LiBDeweyCN.RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics2011;
12:323–38
31.
KhimaniAHMhashilkarAMMikulskisAO’MalleyMLiaoJGolenkoEELottST.Housekeeping genes in cancer: normalization of array data. Biotechniques2005;
38:739–45
32.
MurphyCLPolakJM.Differentiating embryonic stem cells: GAPDH, but neither HPRT nor beta-tubulin is suitable as an internal standard for measuring RNA levels. Tissue Eng2002;
8:551–9
33.
Serazin-LeroyVDenis-HenriotDMorotMde MazancourtPGiudicelliY.Semi-quantitative RT-PCR for comparison of mRNAs in cells with different amounts of housekeeping gene transcripts. Mol Cell Probes1998;
12:283–91
34.
SavonetVMaenhautCMiotFPirsonI.Pitfalls in the use of several “housekeeping” genes as standards for quantitation of mRNA: the example of thyroid cells. Anal Biochem1997;
247:165–7
35.
GabrielssonBGOlofssonLESjogrenAJernasMElanderALonnMCarlssonLM.Evaluation of reference genes for studies of gene expression in human adipose tissue. Obes Res2005;
13:649–52
36.
GorzelniakKJankeJEngeliSSharmaAM.Validation of endogenous controls for gene expression studies in human adipocytes and preadipocytes. Horm Metab Res2001;
33:625–7
37.
HamalainenHKTubmanJCVikmanSKyrolaTYlikoskiEWarringtonJALahesmaaR.Identification and validation of endogenous reference genes for expression profiling of T helper cell differentiation by quantitative real-time RT-PCR. Anal Biochem2001;
299:63–70
38.
BustinSABenesVGarsonJAHellemansJHuggettJKubistaMWittwerCT.The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem2009;
55:611–22
39.
YeJCoulourisGZaretskayaICutcutacheIRozenSMaddenTL.Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics2012;
13:134–44
40.
LeeCKimJShinSGHwangS.Absolute and relative QPCR quantification of plasmid copy number in Escherichia coli. J Biotechnol2006;
123:273–80
41.
ChenDPanXXiaoPFarwellMAZhangB.Evaluation and identification of reliable reference genes for pharmacogenomics, toxicogenomics, and small RNA expression analysis. J Cell Physiol2011;
226:2469–77
42.
ChongWPReddySGYusufiFNLeeDYWongNSHengCK.Metabolomics-driven approach for the improvement of Chinese hamster ovary cell growth: overexpression of malate dehydrogenase II. J Biotechnol2010;
147:116–21
43.
TabuchiHSugiyamaTTanakaSTainakaS.Overexpression of taurine transporter in Chinese hamster ovary cells can enhance cell viability and product yield, while promoting glutamine consumption. Biotechnol Bioeng2010;
107:998–1003
44.
BohrenKMBullockBWermuthBGabbayKH.The aldo-keto reductase superfamily. cDNAs and deduced amino acid sequences of human aldehyde and aldose reductases. J Biol Chem1989;
264:9547–51
45.
KurahashiTKwonMHommaTSaitoYLeeJTakahashiM.Reductive detoxification of acrolein as a potential role for aldehyde reductase (AKR1A) in mammals. Biochem Biophys Res Commun2014;
452:136–41
46.
GabbayKHBohrenKMMorelloRBertinTLiuJVogelP.Ascorbate synthesis pathway: dual role of ascorbate in bone homeostasis. J Biol Chem2010;
285:19510–20
47.
AkiharaRHommaTLeeJYamadaKMiyataSFujiiJ.Ablation of aldehyde reductase aggravates carbon tetrachloride-induced acute hepatic injury involving oxidative stress and endoplasmic reticulum stress. Biochem Biophys Res Commun2016;
478:765–71
LongCXiaoYLiSTangXYuanZBaiY.Identification of optimal endogenous reference RNAs for RT-qPCR normalization in hindgut of rat models with anorectal malformations. Peer J2019;
7:e6829
50.
ZhangYHanXChenSZhengLHeXLiuM.Selection of suitable reference genes for quantitative real-time PCR gene expression analysis in salix matsudana under different abiotic stresses. Sci Rep2017;
7:40290
51.
FurrowEPfeiferRJOsborneCALulichJP.An APRT mutation is strongly associated with and likely causative for 2,8-dihydroxyadenine urolithiasis in dogs. Mol Genet Metab2014;
111:399–403
52.
Garcia-HuertasPMejia-JaramilloAMGonzalezLTriana-ChavezO.Transcriptome and functional genomics reveal the participation of adenine phosphoribosyltransferase in Trypanosoma cruzi resistance to benznidazole. J Cell Biochem2017;
118:1936–45
53.
LingHWuQGuoJXuLQueY.Comprehensive selection of reference genes for gene expression normalization in sugarcane by real time quantitative rt-PCR. PLoS One2014;
9:e97469
54.
XueBGuoJQueYFuZWuLXuL.Selection of suitable endogenous reference genes for relative copy number detection in sugarcane. Int J Mol Sci2014;
15:8846–62
55.
Ian’shina Dd MalyginAAKarpovaGG.[Binding of human ribosomal protein S16 with the 18S rRNA fragment 1203-1236/1521]. Molekul Biol2007;
41:1023–30
56.
LivramentoKGDFreitasNC.Gene expression profile analysis is directly affected by the selected reference gene: the case of Leaf-Cutting atta sexdens. Insects2018;
9:18–33
57.
de Melo NetoOPda Costa LimaTDCMerloKCRomaoTPRochaPOAssisLA.Phosphorylation and interactions associated with the control of the leishmania Poly-A binding protein 1 (PABP1) function during translation initiation. RNA Biol2018;
15:739–55
58.
SmithRWBleeTKGrayNK.Poly(a)-binding proteins are required for diverse biological processes in metazoans. Biochem Soc Trans2014;
42:1229–37
59.
AssrirNFilholOGalissonFLipinskiM.HIRIP3 is a nuclear phosphoprotein interacting with and phosphorylated by the serine-threonine kinase CK2. Biol Chem2007;
388:391–8
60.
LorainSQuivyJPMonier-GavelleFScampsCLecluseYAlmouzniGLipinskiM.Core histones and HIRIP3, a novel histone-binding protein, directly interact with WD repeat protein HIRA. Mol Cell Biol1998;
18:5546–56
