BioMed Research International
Volume 2020 (2020), Article ID 8926120, 11 pages
https://doi.org/10.1155/2020/8926120

Research Article

Selection of Suitable Reference Genes for qPCR Gene Expression Analysis of HepG2 and L02 in Four Different Liver Cell Injured Models

Jiyu Chen,1,2 Zhenzhen Bao,¹ Yanli Huang,² Zhenglong Wang,² and Yucheng Zhao³

¹School of pharmacy, Jiangsu Health Vocational College, Nanjing, Jiangsu, China
²School of Life Science and Technology, China Pharmaceutical University, Nanjing, China
³Jiangsu Key Laboratory of Bioactive Natural Product Research and State Key Laboratory of Natural Medicines, School of Traditional Chinese Pharmacy, China Pharmaceutical University, Nanjing, Jiangsu, China

Correspondence should be addressed to Yucheng Zhao; zhaoyucheng1986@126.com

Received 1 December 2019; Accepted 2 June 2020; Published 14 July 2020

Academic Editor: Sercan Ergün

Copyright © 2020 Jiyu Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Quantitative real-time PCR (qPCR) has become a widely used approach to analyze the expression level of selected genes. However, owing to variations in cell types and drug treatments, a suitable reference gene should be selected according to special experimental design. In this study, we investigated the expression level of ten candidate reference genes in hepatoma carcinoma cell (HepG2) and human hepatocyte cell line (L02) treated with ethanol (EtOH), hydrogen peroxide (H₂O₂), acetaminophen (APAP), and carbon tetrachloride (CCl₄), respectively. To analyze raw cycle threshold values (Cp values) from qPCR run, three reference gene validation programs, including Bestkeeper, geNorm, and NormFinder, were used to evaluate the stability of ten candidate reference genes. The results showed that TATA-box binding protein (TBP) and tubulin beta 2a (TUBB2a) presented the highest stability for normalization under different treatments and were regarded as the most suitable reference genes of HepG2 and L02. In addition, this study not only identified the most stable reference genes of each treatment, but also suggested that β-actin (ACTB), glyceraldehade-3-phosphate dehydrogenase (GAPDH), tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein zeta (YWHAZ), and beta-2 microglobulin (B2M) were the least stable reference genes in HepG2 and L02. This work was the first report to systematically explore the stability of reference genes in injured models of HepG2 and L02.

1. Introduction

Quantitative real-time PCR (qPCR) is commonly used in analyzing gene expression levels owing to its credible precision and high-throughput competence [1, 2]. However, quantitative analysis of gene expressions is unavoidably affected by several factors such as sample amount, cell activity, RNA integrity, and cDNA quality [3–5]. Hence, in order to avoid quantitative errors and obtain a reliable experimental result, one or several reference genes should be applied as a suitable endogenous control for quantitative measurement of gene expression. Some literature indicated that at least three reference genes were needed to normalize the analysis of qPCR [6, 7]. In addition, numerous reports affirmed that the stability of reference genes might change based on various experimental designs and samples [8, 9]. Hence, a stable reference gene, which ensures the stability in various experimental conditions, should be identified.

Traditionally, GAPDH and ACTB are most frequently used for normalization; however, they have been demonstrated unsuitable for internal control because their stability varies in different experiments and samples [10–12]. For instance, in Li et al.’s [13] study, the mean Cp value of GAPDH was 23.88 in H2O2 treated human umbilical vein endothelial cells (HUVEC), while in cytokines treated HUVEC, the mean Cp values of GAPDH was distinctly below 20 [14]. Undoubtedly, the varied expressions of reference genes lead to the inaccuracy of results. Fortunately, an increasing number of researches have focused on selecting and identifying suitable reference genes of humans [15], plant [16], cell line [17], algae [18], animal [19], and bacteria [20]. However, a systematic research about the validation of suitable reference genes for liver cell (HepG2 and L02) injured models has not been reported.

HepG2 is an immortalized human hepatoma cell line, and L02 is an immortalized hepatocyte cell line [21, 22]. Additionally, HepG2 and L02 are widely accepted model systems for investigating hepatotoxicity, intracellular trafficking, and drug targeting in vitro [23–25]. Owing to the stability of reference genes varied with drug-treatments and differed in different cell lines [26]. Hence, in this study, we chose four liver cell injured models commonly used in pharmacology and toxicology: ethanol (EtOH) [27], hydrogen peroxide (H₂O₂) [28], acetaminophen (APAP) [29], and carbon tetrachloride (CCl₄) [30], which represented alcoholic liver injury (EtOH), hepatic oxidative stress (H₂O₂), drug liver injury (APAP), and acute liver damage (CCl₄), respectively, to find the most appropriate reference genes in different cell injured models of HepG2 and L02.

In this study, ten candidate reference genes, ACTB, B2M, GAPDH, TUBB2a, hypoxanthine phosphoribosyltransferase 1 (HPRT1), succinate dehydrogenase complex flavoprotein subunit A (SDHA), TBP, YWHAZ, cytochrome c isoform 1 (CYC1), and glucuronidase beta (GUSB), were selected to investigate the most stable reference genes for normalization in liver cell injured models. To evaluate the stability of candidate reference genes comprehensively, four types of experimental treatments (EtOH, H₂O₂, APAP, and CCl₄) were investigated in two cell types (HepG2 and L02) in vitro. In addition, in order to analyze the correlation between different concentrations of drug treatment and expression levels of reference genes, we selected three groups of different concentrations (low dose group, middle dose group, and high dose group) for each treatment. All concentrations were chosen based on previous studies [31–36] which had performed a cell viability assay proving varying degrees cytotoxicity. To analyze the original data, three statistical algorithms named, geNorm [37], NormFinder [38], and Bestkeeper [39] were used based on the manufacturers’ procedures. The calculation results of three kinds of software showed that TBP and TUBB2a were the most stable ones among all treatments. Moreover, geNorm was also used to calculate the optimal number of reference genes needed for normalization, and the results showed that it was sufficient for accuracy normalization to choose two reference genes in most groups. To our knowledge, this is the first study about the selection of the best reference genes in liver cell injured models, which would provide a proper choice of reference genes and guarantee a dependable result in liver cell injured model research.

2. Materials and Methods

2.1. Reagents

The ethanol (EtOH, 99.5% pure), hydrogen peroxide (H2O2, 30.0% pure), acetaminophen (APAP, 99.5% pure), and carbon tetrachloride (CCl₄, 99.5% pure) were purchased from Aladdin Biochemical Technology Co., Ltd (Shanghai, China); Penicillin and streptomycin were obtained from Beyotime Institute of Biotechnology (Shanghai, China); Trypsin-EDTA Solution was purchased from Sangon Biotech (Shanghai, China).

2.2. Cell Culture and Treatment

The hepatoma carcinoma cells (HepG2) were obtained from the American Type Culture Collection (HB-8065), and the human hepatocyte cells were purchased from the Cell Bank of Type Culture Collection of the Chinese Academy of Sciences. Cells were grown in Dulbecco’s Modified Eagle Medium (DMEM, Gibco), containing 100 U/ml penicillin-streptomycin and 10% fetal bovine serum (FBS, Bioind) under standard conditions (37°C and 5% CO₂). The cells were grown to 80% confluence and then passaged using Trypsin-EDTA Solution (0.25% Trypsin, 0.02% EDTA). All cells were divided into four groups for treatments: (a) control group; (b) low dose group; (c) middle dose group; (d) high dose group. For HepG2, cells were treated with four different treatments, including ethanol (100 mM, 200 mM, 400 mM), H₂O₂ (200 μM, 400 μM, 800 μM), APAP (2.5 mM, 5 mM, 10 mM), and CCl₄ (0.1%, 0.2%, 0.4%). For L02, cells were treated with four different treatments, including ethanol (100 mM, 200 mM, 400 mM), H₂O₂ (100 μM, 200 μM, 400 μM), APAP (2.5 mM, 5 mM, 10 mM), and CCl₄ (0.05%, 0.1%, 0.2%). The CCl₄ were dissolved into 0.25% DMSO and then were added to the serum-free DMEM; the ethanol, H₂O₂, and APAP were dissolved into serum-free DMEM directly. Cells were seeded in six-well plates before being subjected to treatments. For all groups, cells were incubated in the presence or absence of various treatments and different concentrations for 24 h.

2.3. Screening of Candidate Reference Genes and Primer Design

According to previous studies [9, 40], a total of ten candidate reference genes (ACTB, B2M, GAPDH, TUBB2a, HPRT1, SDHA, TBP, YWHAZ, CYC1, and GUSB) were selected to ascertain the best reference genes of HepG2 and L02 in liver cell injured conditions. The nucleotide sequences were downloaded, using Primer 5 to design primers. Full gene names and accession numbers, as well as primer length and intron-spanning primers, were listed in Table 1. The data of qPCR were repeated three times of biological and technical replicates.

Table 1: Details of the ten candidate reference genes and primers used in the qPCR.

Gene
Description
Primer: forward/reverse(5-3)
Length (bp)
Accession number

ACTB
β-Actin
F: AAGGCCAACCGCGAGAAGAT
R: GCCAGAGGCGTACAGGGATA
102
NM_001101

B2M
Beta-2 microglobulin
F: GTTTACTCACGTCATCCAGC
R:AGACAAGTCTGAATGCTCCA
141
NM_004048

GAPDH
Glyceraldehade-3-phosphate dehydrogenase
F: GCCTCCTGCACCACCAACTG
R: CCATCACGCCACAGTTTCCC
149
NM_002046

TUBB2a
Tubulin beta 2a
F: AACGCCACCCTCTCTGTCCA
R: GCCGACACCAGGTGGTTGAG
143
NM_001069

HPRT1
Hypoxanthine phosphoribosyltransferase 1
F: ACTGAACGTCTTGCTCGAGA
R: TGATGTAATCCAGCAGGTCA
112
NM_000194

SDHA
Succinate dehydrogenase complex flavoprotein subunit A
F: AAAGATCACGTCTACCTGCA
R: CATGTTATAATGCACGGTGG
150
NM_004168

TBP
TATA-box binding protein
F: GTTCAGCAGTCAACGTCCCA
R: TCATGGGGGAGGGATACAGT
127
NM_003194

YWHAZ
Tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein zeta
F: CAGGCTGAGCGATATGATGA
R: CCTACGGGCTCCTACAACAT
126
NM_003406

CYC1
Cytochrome c isoform 1
F: CCAAAACCATACCCCAACAG
R: AGTCCTCACCACCATGCCTA
103
NM_001916

GUSB
Glucuronidase beta
F: GTTCCTTTTGCGAGAGAGAT
R: ACACGCAGGTGGTATCAGTC
124
NM_000181

2.4. Total RNA, DNA Extraction and cDNA Synthesis

Total RNA was extracted from HepG2 and L02 and purified using the RNAiso Plus total RNA kit (TransGen Biotech, Dalian, China) according to the manufacturer’s instruction. And then, DNase I (Takara, Dalian, China) was added to the sample to eliminate DNA contamination for RNA purity. The purity of the total RNA was assessed by measuring the absorbance ration at 260/280 nm of the samples. In addition, the quality of the RNA was confirmed by agarose gel electrophoresis. Purified RNA was reverse transcribed immediately after extraction. For qPCR experiments, HiScript® Q RT SuperMix for qPCR Kit (Vazyme, Nanjing, China) and a quantity of 1 μg total RNA were added into a 20 μl reaction volume to synthesize cDNA.

2.5. Quantitative Real-Time PCR

The sample reaction was run in 96-well plate. Real-time quantitative PCR with AceQ qPCR SYBR Green Master Mix (Vazyme, Nangjing, China) was performed at LightCycler 480 (Roche Molecular Biochemicals, Mannheim, Germany). Each reaction system was 20 μl, respectively. AceQ qPCR SYBR Green Master Mix 10 μl, forward and reverse primers were 0.4 μM each, template cDNA 2 μl, and added ddH₂O to the final volume of 20 μl. Each sample was repeated 3 times. The optimizing reaction conditions of real-time quantitative PCR as follows: 1 cycle of 95°C for 5 min, 40 cycles of 95°C for 10 sec, and then 60°C for 30 sec.

2.6. Analysis of Reference Genes Stability

In order to evaluate the stability of ten selected reference genes, three reference gene validation programs (geNorm, NormFinder, and BestKeeper) were used under the manufacture’s instruction. NormFinder was applied to calculate the stability value () for finding the steadiest candidate genes. For geNorm, the calculation could determine the optimal number of reference genes and, similar to geNorm, evaluate the stability of candidate genes. BestKeeper was based on the coefficient of variance (CV) and the standard deviation (SD) of the Cp values to assess the steadiness of reference genes. Three biological and technical repeats were used for different experimental conditions.

3. Results

3.1. Verification of the Primers Specificity

We used PCR to identify the specificity of the designed primers by agarose gel electrophoresis, as S2 and S3 Figs shows, the single band and peak of a melting curve indicated primers possessed the good specificity.

3.2. Evaluating Expression of the Reference Gene

The most suitable reference genes would have stable expression levels in various treatments and concentrations. And the Cp value of ten candidate reference genes underwent diver treatments were listed in Figure 1, ranging from 14.7 to 34.49 (HepG2 14.7 to 34.49, L02 14.85 to 34.48), suggesting that they have a noticeable variance in expression level. Particularly, most of the Cp values were in a range of 20 to 27. ACTB, B2M, GAPDH, HPRT1, and YWHAZ expressed lower Cp value around 20, while the rest of the genes showed that higher Cp value was greater than 25, especially the GUSB, which had the highest mean Cp values (HepG2 , L02 ). Notably, ACTB showed the minimal change of Cp values from 18.44 to 24.90 under different treatments in HepG2, and meanwhile, the Cp values of HPRT1 from 21.46 to 27.02 showed the low variation in L02, suggesting that the two genes might have a constant expression under various treatments and could be a suitable reference gene. In short, Cp values, combined with box-plot, presented the expression of the reference genes, and as well provided us a general understanding of gene stability.

(a)

(b)

Figure 1: Expression levels of ten reference genes (ACTB, B2M, GAPDH, TUBB2a, HPRT1, SDHA, TBP, YWHAZ, CYC1, GUSB) in HepG2 (a) and L02 (b). Squares in the middle of the box represent the mean values, horizontal lines in the box represent the median, and whiskers represent the highest and lowest values.

3.3. Expression Stability of Candidate Reference Genes

The data obtained from different treatments (wild-type APAP, CCl₄, ethanol, and H₂O₂) and each reference gene were analyzed with three Excel-based programs (geNorm, NormFinder, and BestKeeper) for further evaluation on the stability of putative reference genes.

3.4. geNorm Analysis

To ascertain the stability of candidate reference genes, geNorm was applied to evaluate the expression stability measurement () value by Cp values of each gene in groups. According to the analysis of geNorm, genes with the highest values were considered as the least stable ones and the lowest the most. As shown in Figure 2 and S1 Figure, different reference genes had different values in different treatments. For instance, in the L02 groups, TBP with the value of 0.51 in 400 μM H₂O₂ treatment would be the steadiest reference genes, while the GUSB was more than twice TBP, with value was 1.37 in the same treatment. More interestingly, even the same reference gene had different expressions in different treatments. In the HepG2 groups, the gene with the lowest values in 10 mM APAP treatment was HPRT1, which owned the highest values in 200 μM H₂O₂ treatment, meaning that HPRT1 was the most stable genes in 10 mM APAP treatment and the least stable ones in 200 μM H₂O₂ treatment.

Figure 2: Expression stability of the reference genes in HepG2 evaluated by geNorm values represents the average expression stability. From left to right, the value of decreased in turn, indicating the stability gradually increased. Smaller value means higher stability. The control group, ethanol, hydrogen peroxide, acetaminophen, and carbon tetrachloride were abbreviated to WT, EtOH, H₂O₂, APAP, and CCl₄, respectively.

3.5. NormFinder Analysis

NormFinder was used to evaluate the optimal gene for normalization in each experiment. The raw Cp values obtained from qPCR were firstly log-transformed and used as the input value for the NormFinder, and then used to analyze the expression stability according to the similarity of the expression profiles of candidate genes. Genes with lower values that were close to zero were regarded as the best candidate ones. As shown in Table 2 and S1 Table, the rank of values was increasing from top to bottom of the table, whereas genes on the top of the table were the most stable reference genes. Therefore, the most stable candidate genes could be easily found from the table. In the HepG2 group, CYC1 (7 times to be the top 3 candidate genes in 13 treatments) and HPRT1 (8 times to be the top 3 in 13 treatments) were considered as the most stable reference genes; the results were similar to that of geNorm analysis. Nevertheless, in the L02 group, the results of geNorm and NormFinder analysis seemed to be different. For instance, in the NormFinder analysis, TUBB2a (9 times to be the top 3 in all treatments) was the steadiest ones, while GAPDH, the most stable genes in the geNorm analysis, appeared only once to be the top 3 of all treatments. Hence, the third analysis method should be used.

Table 2: Expression stability values of ten candidate reference genes in HepG2 analyzed by NormFinder.

Rank
WT
EtOH
100 mM
EtOH
200 mM
EtOH
400 mM
H₂O₂
200 μM
H₂O₂
400 μM
H₂O₂
800 μM
APAP
2.5 mM
APAP
5 mM
APAP
10 mM
CCl₄
0.1%
CCl₄
0.2%
CCl₄
0.4%

1
HPRT1
0.011
HPRT1
0.004
CYC1
0.007
HPRT1
0.008
B2M
0.018
CYC1
0.013
TBP
0.011
HPRT1
0.013
GUSB
0.004
HPRT1
0.012
TBP
0.015
TBP
0.011
HPRT1
0.005

2
TUBB2a
0.012
TBP
0.005
TBP
0.009
TUBB2a
0.013
TUBB2a
0.019
TUBB2a
0.017
GAPDH
0.011
GUSB
0.015
CYC1
0.010
B2M
0.018
CYC1
0.017
ACTB
0.013
CYC1
0.008

3
CYC1
0.013
TUBB2a
0.005
SDHA
0.011
B2M
0.014
CYC1
0.020
TBP
0.019
SDHA
0.012
SDHA
0.018
TUBB2a
0.014
SDHA
0.019
SDHA
0.020
B2M
0.021
GUSB
0.009

4
B2M
0.017
YWHAZ
0.012
B2M
0.011
TBP
0.016
TBP
0.027
HPRT1
0.024
HPRT1
0.015
YWHAZ
0.022
ACTB
0.020
TUBB2a
0.020
HPRT1
0.021
TUBB2a
0.021
TUBB2a
0.016

5
TBP
0.018
SDHA
0.013
HPRT1
0.013
SDHA
0.018
GUSB
0.030
GUSB
0.025
B2M
0.017
CYC1
0.023
HPRT1
0.021
CYC1
0.022
TUBB2a
0.024
CYC1
0.021
SDHA
0.019

6
SDHA
0.019
ACTB
0.014
GAPDH
0.015
CYC1
0.020
SDHA
0.030
YWHAZ
0.027
TUBB2a
0.017
TBP
0.024
TBP
0.022
ACTB
0.024
GUSB
0.024
GUSB
0.022
TBP
0.024

7
ACTB
0.020
CYC1
0.022
GUSB
0.015
ACTB
0.025
YWHAZ
0.037
ACTB
0.031
CYC1
0.017
ACTB
0.028
YWHAZ
0.025
GUSB
0.025
YWHAZ
0.028
YWHAZ
0.022
B2M
0.029

8
GUSB
0.021
B2M
0.032
TUBB2a
0.016
GUSB
0.027
ACTB
0.039
B2M
0.037
ACTB
0.019
GAPDH
0.029
B2M
0.029
TBP
0.029
GAPDH
0.032
SDHA
0.024
GAPDH
0.031

9
YWHAZ
0.023
GAPDH
0.041
ACTB
0.020
YWHAZ
0.028
HPRT1
0.046
SDHA
0.043
GUSB
0.022
TUBB2a
0.030
SDHA
0.031
GAPDH
0.035
B2M
0.034
HPRT1
0.029
ACTB
0.037

10
GAPDH
0.029
GUSB
0.042
YWHAZ
0.023
GAPDH
0.041
GAPDH
0.061
GAPDH
0.048
YWHAZ
0.024
B2M
0.032
GAPDH
0.039
YWHAZ
0.037
ACTB
0.037
GAPDH
0.044
YWHAZ
0.039

3.6. BestKeeper Analysis

BestKeeper was an Excel-based tool used to analyze the expression stability of the candidate reference gene. The standard deviation (SD) and coefficient of variation (CV) were calculated by BestKeeper to assess the stability of candidate reference genes in each group. Genes with the lowest SD and CV would be the most stable reference ones. As shown in Table 3 and S2 Table, the () values of ten candidate reference genes progressively increased from top to bottom of tables, showing their decreasingly stability. As an example, TBP was listed on the top of Table 3, with a lower () value of (), representing the most stable genes in 800 μM H₂O₂ induced oxidative stress in HepG2, and meanwhile, GUSB, having a () value of (), was listed at the bottom of the table. In HepG2 groups, some reference genes, namely, TBP, CYC1, and TUBB2a, might be the best suitable genes for the reason that they were listed on top 3 of the rank in majority treatments. Similarly, CYC1 and TBP occupied most of the top 3 in the table, suggesting that the two candidate genes would be the steadiest genes in L02 treatments.

Table 3: Expression stability values of the reference genes calculated by BestKeeper in HepG2.

Rank
WT
EtOH
100 mM
EtOH
200 mM
EtOH
400 mM
H₂O₂
200 μM
H₂O₂
400 μM
H₂O₂
800 μM
APAP
2.5 mM
APAP
5 mM
APAP
10 mM
CCl₄
0.1%
CCl₄
0.2%
CCl₄
0.4%

1
GUSB

HPRT1

SDHA

B2M

CYC1

CYC1

TBP

SDHA

CYC1

GUSB

GUSB

HPRT1

TBP

2
ACTB

CYC1

TBP

TBP

TUBB2a

YWHAZ

CYC1

TBP

GUSB

CYC1

CYC1

GUSB

SDHA

3
TUBB2a

TBP

CYC1

SDHA

SDHA

GUSB

B2M

TUBB2a

TUBB2a

TUBB2a

HPRT1

TBP

CYC1

4
CYC1

SDHA

GUSB

YWHAZ

TBP

HPRT1

SDHA

GUSB

HPRT1

ACTB

B2M

CYC1

TUBB2a

5
B2M

TUBB2a

ACTB

HPRT1

ACTB

TUBB2a

HPRT1

YWHAZ

B2M

TBP

TBP

YWHAZ

HPRT1

6
TBP

ACTB

B2M

CYC1

B2M

TBP

GAPDH

HPRT1

YWHAZ

B2M

SDHA

TUBB2a

GUSB

7
HPRT1

YWHAZ

TUBB2a

TUBB2a

GUSB

ACTB

GUSB

CYC1

ACTB

HPRT1

YWHAZ

ACTB

YWHAZ

8
SDHA

GUSB

GAPDH

GUSB

YWHAZ

B2M

TUBB2a

ACTB

TBP

SDHA

TUBB2a

B2M

ACTB

9
YWHAZ

B2M

HPRT1

ACTB

HPRT1

SDHA

ACTB

B2M

SDHA

GAPDH

ACTB

SDHA

B2M

10
GAPDH

GAPDH

YWHAZ

GAPDH

GAPDH

GAPDH

YWHAZ

GAPDH

GAPDH

YWHAZ

GAPDH

GAPDH

GAPDH

3.7. Optimal Numbers of Reference Genes for Normalization

The minimal numbers of reference genes for accurate normalization could also be determined by geNorm, according to the calculation of pairwise variation (variation coefficient, V) between the normalization factors (NF) in various treatment sets using as a criterion cut-off value [37]. Based on this rule, the calculation was listed in Figure 3. As we can see, there were enough to choose two or three reference genes in most treatments of HepG2 and L02 for normalization. Moreover, 10 mM APAP treatment in HepG2, 200 mM EtOH, and 400 mM EtOH in L02, respectively, required four, five, and nine reference genes for normalization.

(a)

(b)

Figure 3: Calculation of the optimal number of reference genes for quantitative analysis using geNorm. Pairwise variation () of reference genes analyzed in different treatments was listed in (a, b) set the cut-off threshold value 0.15, calculating the optimal number of reference genes for precise quantitative in this work. WT, E100 (E200, E400), H100 (H200, H400, H800), A2.5 (A5, A10), and C0.05% (C0.1%, C0.2%, C0.4%), respectively, were the abbreviation for the control group, EtOH-100 mM (200 mM, 400 mM), H₂O₂-100 μM (200 μM, 400 μM, 800 μM), APAP-2.5 mM (5 mM, 10 mM, 20 mM), CCl₄-0.05% (0.1%, 0.2%, 0.4%).

4. Discussion

Quantitative real-time PCR is one of the most accurate and commonly used techniques for analysis of gene transcript levels. Selection of suitable reference gene is indispensable to guarantee the accuracy and consistency of the data and minimize the experimental errors. To confirm the precise expression analysis of putative genes, numerous steady reference genes have been verified in different experimental designs [41, 42]. Hence, the purpose of this study was to investigate the expression stability of ten candidate reference genes (ACTB, B2M, GAPDH, TUBB2a, HPRT1, SDHA, TBP, YWHAZ, CYC1, GUSB) in two in-vitro cell types, namely, HepG2 cells and L02 cells, and determine the optimal candidate genes under the treatment of alcoholic liver injury (EtOH), hepatic oxidative stress (H₂O₂), drug liver injury (APAP), and acute liver damage (CCl₄). After that, raw data was input and calculated in three Excel-based programs: geNorm, NormFinder, and BestKeeper.

The data from qPCR run of ten candidate genes were listed in Figure 1, where the expression level and mean Cp values of the candidate genes ranging from 14.7 to 34.49 (HepG2 14.7 to 34.49, L02 14.85 to 34.48) could be easily seen. However, the scope of the Cp values of some selected genes inconsistent with the previous study [43, 44] might be a result of liver damage treatments. Based on the fact that the genes with the highest expression levels owned minimal Cp values, GAPDH with the lowest mean Cp values of 17.47 in HepG2 and 17.25 in L02 means that the GAPDH was abundantly distributed in the two cell types. Considering that a wide distribution range tends to be low stability and moreover, Cp values with low variation would be more suitable for reference gene selection. The variation of Cp values suggested that ACTB and HPRT1 were the best reference genes in HepG2 and L02, while HPRT1 and YWHAZ were the least ones. The verification above was a little different from the calculation of the three Excel-based programs. For instance, the Cp values of ACTB and HPRT1 might not fluctuate significantly, but in the calculations of the three kinds of software, the two genes appeared at the bottom of ranking frequently, which reflected more instability. Accordingly, more calculation results needed to be combined.

Furthermore, some literature reported that the expression level of reference genes would change under different concentration treatments [26, 45]. Hence, we investigated the stability of the reference genes in the same treatment of different concentrations with the results of three software calculations. Since the results of three software varied based on different algorithms, we selected the top five (six in some groups) reference genes for each calculation to evaluate them comprehensively. We identified the top five genes under three concentrations, and then found that genes were common in geNorm, NormFinder, and Bestkeeper. For instance, in CCl₄ of three concentrations treated L02 cells, GUSB, TUBB2a, and TBP commonly appeared to be the top five of Bestkeeper calculation results, while SDHA, GAPDH, and TUBB2a of geNorm and SDHA, TBP, and TUBB2a of NormFinder, respectively. Apparently, TUBB2a appeared to be one of the top five of each calculation. Hence, we recommended TUBB2a as the most stable reference genes in CCl₄ treated L02 cells. Likewise, in L02 cells, TBP was considered as the most stable genes in EtOH, H₂O₂, and APAP treatments, and TUBB2a was the steadiest in CCl₄ treatments. In HepG2 cells, we suggested TBP being the most stable reference genes in EtOH and H₂O₂ treatments, and GUSB and CYC1 in APAP and CCl₄ treatments, respectively. In the same way, we evaluated the least stable reference genes in each group. On the whole, ACTB, GAPDH, YWHAZ, and B2M always ranked the last, meaning that they were considered as the least stable. Hence, we did not recommend ACTB, GAPDH, YWHAZ, or B2M as internal control for normalization. However, in Bridget’s study [44], they identified GAPDH as the most stable reference genes in APAP treated HepG2, opposite to our research, which was mainly because they evaluated the stability only using geNorm, which might lead to inaccurate results.

Based on the analysis above, we found the best suitable reference genes in different treatments. Nevertheless, how many reference genes were required to be chosen for optimal data normalization required further investigation. Hence, we chose the geNorm software, which could calculate the optimal number of reference genes in a qPCR experiment to solve these problems. According to the handbook [37], the V score of 0.15 as a criterion value was recommended, and an additional gene was included until () was lower than 0.15. In this study, the results showed that the majority of pairwise values were lower than 0.15 after a total of 26 treated groups. and the calculation was shown in Figure 3, among which 23 out of 26 groups with a low , signifying the inclusion of additional reference genes was unnecessary in those 23 groups. Hence, two or three reference genes would suffice for reliable normalization in these groups above. Unfortunately, not all groups had a suitable value of lower than 0.15. For instance, in the 400 mM EtOH treated L02 group, all of pairwise values were greater than 0.15. Vandesompele et al. recommended that it was a waste of resources to quantify more genes than necessary, and hence, the value (V6/7 and V9/10 values were close to 0.15) indicated that six reference genes should be a good choice for normalization in 400 mM treated L02 group.

5. Conclusion

In this study, 10 candidate genes were selected and evaluated in two types of liver cells (HepG2 and L02) for four types of liver cell injured treatments using the three different algorithms, namely, BestKeeper, geNorm, and NormFinder. To the best of our knowledge, this was the first systematic selection of reference genes in the liver cell injured model and laid the basis for further research in HepG2 and L02. Based on the analysis, we identified the best reference genes of HepG2 and L02 under the treatments of EtOH, H₂O₂, APAP, and CCl₄. The results of gene expression revealed that TBP and TUBB2a were the most stable reference genes for normalization in different treatments. On one hand, in the HepG2, the most stable reference genes of EtOH and H₂O₂ treatments were TBP, while GUSB and CYC1 were, respectively, the most suitable reference genes of APAP and CCl₄ treatments. In the L02, TBP was identified as the most stable reference genes of EtOH, H₂O₂, and APAP treatments, while TUBB2a was the steadiest reference genes of CCl₄ treatment. On the other hand, ACTB, GAPDH, YWHAZ, and B2M were the least stable reference genes in EtOH, H₂O₂, APAP, and CCl₄ treated HepG2 and L02. In short, our study provided a credible selection of reference gene in HepG2 and L02 injured models.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors did not report any conflict of interest.

Authors’ Contributions

Jiyu Chen, Zhenzhen Bao, and Yanli Huang contributed equally to this work.

Acknowledgments

This work was funded by the Practical Training Program for Young Teachers in Higher Vocational Colleges in Jiangsu Province (2019QYSJPX171) and Natural Science Foundation of the Jiangsu Higher Education Institutions of China (19KJD430005). This project was also funded by the Qinglan project of excellent teaching team in Jiangsu and teaching and research project of Jiangsu Health Vocational College (JKKYTD201701, JKA201812).

Supplementary Materials

S1 Figure: expression stability of the housekeeping genes in L02 evaluated by geNorm values represents the average expression stability. From left to right, the value of decreased in turn, indicating the stability gradually increased. Smaller value means higher stability. The control group, ethanol, hydrogen peroxide, acetaminophen, and carbon tetrachloride were abbreviated to WT, EtOH, H2O2, APAP, and CCl4, respectively. S2 Figure: agarose gel (1%) electrophoresis of the ten candidate housekeeping genes. 1-10 represent ACTB, B2M, GAPDH, TUBB2a, HPRT1, SDHA, TBP, YWHAZ, CYC1, and GUSB, respectively. S3 Figure: melt curves of the ten candidate housekeeping genes. S1 Table: expression stability values of ten candidate housekeeping genes in L02 analyzed by NormFinder. S2 Table: expression stability values of the housekeeping genes calculated by BestKeeper in L02. (Supplementary Materials)

References

S. A. Bustin, V. Benes, J. A. Garson et al., “The MIQE Guidelines: Minimum Information for Publication of Quantitative Real-Time PCR Experiments,” Clinical Chemistry, vol. 55, no. 4, pp. 611–622, 2009.

Y. Zhang, D. Zhang, W. Li, J. Chen, Y. Peng, and W. Cao, “A novel real-time quantitative PCR method using attached universal template probe,” Nucleic Acids Research, vol. 31, no. 20, pp. 123e–1123, 2003.

S. Alemayehu, K. C. Feghali, J. Cowden, J. Komisar, C. F. Ockenhouse, and E. Kamau, “Comparative evaluation of published real-time PCR assays for the detection of malaria following MIQE guidelines,” Malaria Journal, vol. 12, no. 1, 2013.

E. Forsgren, B. Locke, E. Semberg, A. T. Laugen, and J. R. . Miranda, “Sample preservation, transport and processing strategies for honeybee RNA extraction: Influence on RNA yield, quality, target quantification and data normalization,” Journal of Virological Methods, vol. 246, pp. 81–89, 2017.

X. Huang, M. Baumann, L. Nikitina et al., “RNA degradation differentially affects quantitative mRNA measurements of endogenous reference genes in human placenta,” Placenta, vol. 34, no. 7, pp. 544–547, 2013.

Q. Kong, J. Yuan, L. Gao et al., “Evaluation of appropriate reference genes for gene expression normalization during watermelon fruit development,” Plos One, vol. 10, no. 6, article e0130865, 2015.

N. Nicot, J. F. Hausman, L. Hoffmann, and D. Evers, “Housekeeping gene selection for real-time RT-PCR normalization in potato during biotic and abiotic stress,” Journal of Experimental Botany, vol. 56, no. 421, pp. 2907–2914, 2005.

L. H. Heckmann, R. Connon, T. H. Hutchinson, S. J. Maund, R. M. Sibly, and A. Callaghan, “Expression of target and reference genes in Daphnia magna exposed to ibuprofen,” BMC Genomics, vol. 7, no. 1, 2006.

X. Zhang, L. Ding, and A. J. Sandford, “Selection of reference genes for gene expression studies in human neutrophils by real-time PCR,” BMC Molecular Biology, vol. 6, no. 1, 2005.

R. Mori, Q. Wang, K. D. Danenberg, J. K. Pinski, and P. V. Danenberg, “Both beta-actin and GAPDH are useful reference genes for normalization of quantitative RT-PCR in human FFPE tissue samples of prostate cancer,” Prostate, vol. 68, no. 14, pp. 1555–1560, 2008.

P. Murthi, E. Fitzpatrick, A. J. Borg, S. Donath, S. P. Brennecke, and B. Kalionis, “GAPDH, 18S rRNA and YWHAZ are Suitable Endogenous Reference Genes for Relative Gene Expression Studies in Placental Tissues from Human Idiopathic Fetal Growth Restriction,” Placenta, vol. 29, no. 9, pp. 798–801, 2008.

A. Pérez-Rico, F. Crespo, M. L. Sanmartín, A. De Santiago, and J. L. Vega-Pla, “Determining ACTB, ATP5B and RPL32 as optimal reference genes for quantitative RT-PCR studies of cryopreserved stallion semen,” Animal Reproduction Science, vol. 149, no. 3-4, pp. 204–211, 2014.

T. Li, H. Diao, L. Zhao et al., “Identification of suitable reference genes for real-time quantitative PCR analysis of hydrogen peroxide-treated human umbilical vein endothelial cells,” BMC Molecular Biology, vol. 18, no. 1, 2017.

J. J. Garcı́a-Vallejo, B. Van het Hof, J. Robben et al., “Approach for defining endogenous reference genes in gene expression experiments,” Analytical Biochemistry, vol. 329, no. 2, pp. 293–299, 2004.

F. Ohl, M. Jung, A. Radonić, M. Sachs, S. A. Loening, and K. Jung, “Identification and validation of suitable endogenous reference genes for gene expression studies of human bladder cancer,” Journal of Urology, vol. 175, no. 5, pp. 1915–1920, 2006.

Y. Zhao, J. Luo, S. Xu et al., “Selection of Reference Genes for Gene Expression Normalization in Peucedanum praeruptorum Dunn under Abiotic Stresses, Hormone Treatments and Different Tissues,” Plos One, vol. 11, no. 3, article e0152356, 2016.

L. Ceelen, W. De Spiegelaere, M. David et al., “Critical selection of reliable reference genes for gene expression study in the HepaRG cell line,” Biochemical Pharmacology, vol. 81, no. 10, pp. 1255–1261, 2011.

M. Dong, X. Zhang, X. Chi et al., “The validity of a reference gene is highly dependent on the experimental conditions in green alga Ulva linza,” Current Genetics, vol. 58, no. 1, pp. 13–20, 2012.

S. Yuanyuan, S. Qin, X. Rongrong et al., “Reference gene selection for real-time quantitative PCR analysis on ovarian cryopreservation by vitrification in mice,” Journal of Assisted Reproduction and Genetics, vol. 32, no. 8, pp. 1277–1284, 2015.

L. Reiter, A.-B. Kolstø, and A. P. Piehler, “Reference genes for quantitative, reverse-transcription PCR in Bacillus cereus group strains throughout the bacterial life cycle,” Journal of Microbiological Methods, vol. 86, no. 2, pp. 210–217, 2011.

A. L. Roe, J. E. Snawder, R. W. Benson, D. W. Roberts, and D. A. Casciano, “HepG2 Cells: An in Vitro Model for P450-Dependent Metabolism of Acetaminophen,” Biochemical and Biophysical Research Communications, vol. 190, no. 1, pp. 15–19, 1993.

P. de Diesbach, F. N'Kuli, C. Berens et al., “Receptor-mediated endocytosis of phosphodiester oligonucleotides in the HepG2 cell line: evidence for non-conventional intracellular trafficking,” Nucleic Acids Research, vol. 30, no. 7, pp. 1512–1521, 2002.

H. Jaeschke, M. R. McGill, and A. Ramachandran, “Pathophysiological Relevance of Proteomics Investigations of Drug-Induced Hepatotoxicity in HepG2 Cells,” Toxicological Sciences, vol. 121, no. 2, pp. 428–430, 2011.

Y. Lihua, Z. Shi, H. Liu, X. Yang, and K. Wang, “Gadolinium induced apoptosis of human embryo liver L02 cell line by ROS- mediated AIF pathway,” Journal of Rare Earths, vol. 29, no. 2, pp. 178–184, 2011.

D. Faibish, M. Suzuki, and J. D. Bartlett, “Appropriate real-time PCR reference genes for fluoride treatment studies performed in vitro or in vivo,” Archives of Oral Biology, vol. 62, pp. 33–42, 2016.

A. Szuster-Ciesielska and M. Kandefer-Szerszen, “Protective effects of betulin and betulinic acid against ethanol-induced cytotoxicity in HepG2 cells,” Pharmacological Reports, vol. 57, no. 5, pp. 588–595, 2005.

S. H. Sohn, S. K. Kim, Y. O. Kim et al., “A comparison of antioxidant activity of Korean White and Red Ginsengs on H₂O₂-induced oxidative stress in HepG2 hepatoma cells,” Journal of Ginseng Research, vol. 37, no. 4, pp. 442–450, 2013.

Y. He, C. Liu, X. Xia, and L. Liu, “Conformal microcapsules encapsulating microcarrier-L02 cell complexes for treatment of acetaminophen-induced liver injury in rats,” Journal of Materials Chemistry B, vol. 5, no. 10, pp. 1962–1970, 2017.

Y. Ou, S. Zheng, L. Lin, Q. Jiang, and X. Yang, “Protective effect of C-phycocyanin against carbon tetrachloride-induced hepatocyte damage in vitro and in vivo,” Chemico-Biological Interactions, vol. 185, no. 2, pp. 94–100, 2010.

S. Yarnpakdee, S. Benjakul, H. G. Kristinsson, and H. E. Bakken, “Preventive effect of Nile tilapia hydrolysate against oxidative damage of HepG2 cells and DNA mediated by H₂O₂ and AAPH,” Journal of Food Science and Technology, vol. 52, no. 10, pp. 6194–6205, 2015.

D. Ozerkan, N. Ozsoy, and E. Yilmaz, “Vitamin D and melatonin protect the cell’s viability and ameliorate the CCl₄ induced cytotoxicity in HepG2 and Hep₃B hepatoma cell lines,” Cytotechnology, vol. 67, no. 6, pp. 995–1002, 2015.

Q. Ma, L. H. Wang, and J. G. Jiang, “Hepatoprotective effect of flavonoids from Cirsium japonicum DC on hepatotoxicity in comparison with silymarin,” Food & Function, vol. 7, no. 5, pp. 2179–2184, 2016.

Y. Liu, H. Huang, J. Lin, Q. Zhang, J. Tan, and J. Ren, “Suppressive effects of genomic imprinted gene PEG10 on hydrogen peroxide-induced apoptosis in L0₂ cells,” Journal of Huazhong University of Science and Technology [Medical Sciences], vol. 29, no. 6, pp. 705–709, 2009.

C. F. Hu, L. P. Sun, Q. H. Yang, D. X. Lu, and S. Luo, “Ginsenosides from stems and leaves of ginseng prevent ethanol-induced lipid accumulation in human L₀₂ hepatocytes,” Chinese Journal of Integrative Medicine, vol. 23, no. 6, pp. 438–444, 2017.

K. J. S. Kumar, J.-W. Liao, J.-H. Xiao, M. G. Vani, and S.-Y. Wang, “Hepatoprotective effect of lucidone against alcohol-induced oxidative stress in human hepatic HepG2 cells through the up-regulation of HO-1/Nrf-2 antioxidant genes,” Toxicology in Vitro, vol. 26, no. 5, pp. 700–708, 2012.

J. Vandesompele, K. De Preter, F. Pattyn et al., “Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes,” Genome Biology, vol. 3, no. 7, 2002.

C. L. Andersen, J. L. Jensen, and T. F. Orntoft, “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 Research, vol. 64, no. 15, pp. 5245–5250, 2004.

M. W. Pfaffl, A. Tichopad, C. Prgomet, and T. P. Neuvians, “Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper - excel-based tool using pair-wise correlations,” Biotechnology Letters, vol. 26, no. 6, pp. 509–515, 2004.

L. Chen, H. Y. Zhong, J. F. Kuang, J. G. Li, W. J. Lu, and J. Y. Chen, “Validation of reference genes for RT-qPCR studies of gene expression in banana fruit under different experimental conditions,” Planta, vol. 234, no. 2, pp. 377–390, 2011.

M. Jung, A. Ramankulov, J. Roigas et al., “In search of suitable reference genes for gene expression studies of human renal cell carcinoma by real-time PCR,” BMC Molecular Biology, vol. 8, no. 1, 2007.

W. Sun, Y. Jin, L. He, W. C. Lu, and M. Li, “Suitable reference gene selection for different strains and developmental stages of the carmine spider mite, Tetranychus cinnabarinus, using quantitative real-time PCR,” Journal Of Insect Science, vol. 10, no. 208, pp. 1–12, 2010.

S. L. Chua, W. C. S. Too, B. Y. Khoo, and L. L. Few, “UBC and YWHAZ as suitable reference genes for accurate normalisation of gene expression using MCF₇, HCT116 and HepG2 cell lines,” Cytotechnology, vol. 63, no. 6, pp. 645–654, 2011.

V. R. Cicinnati, Q. Shen, G. C. Sotiropoulos, A. Radtke, G. Gerken, and S. Beckebaum, “Validation of putative reference genes for gene expression studies in human hepatocellular carcinoma using real-time quantitative RT-PCR,” BMC Cancer, vol. 8, no. 1, 2008.

B. C. Fox, A. S. Devonshire, M. E. Schutte et al., “Validation of reference gene stability for APAP hepatotoxicity studies in different in vitro systems and identification of novel potential toxicity biomarkers,” Toxicology in Vitro, vol. 24, no. 7, pp. 1962–1970, 2010.

A. K. Gebeh, E. L. Marczylo, A. A. Amoako, J. M. Willets, and J. C. Konje, “Variation in stability of endogenous reference genes in fallopian tubes and endometrium from healthy and ectopic pregnant women,” International Journal of Molecular Sciences, vol. 13, no. 3, pp. 2810–2826, 2012.