INTRODUCTION
An increasing number of studies are being conducted into the role of changes in gene expression on the physiology, biology and reproduction of organisms. These gene expression studies often employ real-time quantitative polymerase chain reaction (RT-qPCR). RT-qPCR can also be used to validate other studies, including microarray data and also more frequently RNA-seq data using a measure called RPKM (Mortazavi et al. 2008). While microarray data and RNA-seq cover a large array of genes, RT-qPCR is still the gold standard for expression studies of target genes. RT-qPCR relies on an internal PCR standard for which to compare changes in gene expression of target genes between experimental conditions. These standards, so called ‘housekeeping’ genes or reference genes, are from genes that are assumed to be necessary for regular physiological maintenance of the organism, and are considered to not alter during the experimental manipulations or different life history stages. Studies have shown that equal expression of ‘standard’ housekeeping genes (HKGs) can not be assumed and must be tested for each organism under study (e.g., Ransbotyn and Reusch 2006, Hibbeler et al. 2008, Wu et al. 2013).
Bostrychia moritziana is a red alga and a member of the Florideophyceae. The class Florideophyceae contains approximately 90% of all described red algal species (Brodie and Zuccarello 2007). Bostrychia has been proposed as a model system for evolutionary studies of red algae (Zuccarello and West 2011). Like many Florideophyceae, Bostrychia has an alternation of generations between haploid gametophytes (separate male and female plants) and morphologically identical free-living diploid sporophytes (tetrasporophyte), that produce spores by meiosis. Like nearly all red algae, Bostrychia has an alternate stage between the gametophytes and the free-living sporophyte, producing what has been called a triphasic red alga life cycle (Saunders and Hommersand 2004). This alternate stage is the carposporophyte, the parasitic development of the zygote on the female plant (Kamiya and Kawai 2002). The genetics of sexual differentiation and transcriptomic difference between the various stages is poorly known in red algae and the genes involved are mostly unexplored (Asamizu et al. 2003, Kamiya et al. 2011). At the proteomic level only 3.7-7.4% of proteins are differentially expressed in males versus females (Kim et al. 2008). Without an expressed sequence tag (EST) dataset and understanding of gene homology, it is difficult to hypothesize on the function of these proteins.
Transcriptomics is a powerful way to study an array of genes that could be important in various biological process. To find genes that are over- or under-expressed in particular stages, and therefore could be causative, it is critical that a good internal standard is established. These internal standards (HKGs) ideally are expressed at near similar levels in all experimental conditions, so to that changes in expression of target genes can be calibrated for handling errors (normalization). HKGs have been evaluated for the red alga Pyropia yezoensis (Wu et al. 2013) and Chondrus crispus (Kowalczyk et al. 2014) and these two studies suggested different sets of genes for the appropriate normalization of expression. P. yezoensis, in the class Bangiophyceae, shares a common ancestor with C. crispus (Floridephyceae) over 1.2 billion years ago (Yoon et al. 2004), it is likely that substantial changes in expression stability have occurred in this time. So what are the appropriate HKGs to use? We undertook the screening of eleven commonly used candidate HKGs in different life history stages of Bostrychia and in an asexual isolate to establish a good candidate gene. These data will add to our understanding on which genes are appropriate within the Florideophyceae for RT-qPCR calibration.
MATERIALS AND METHODS
Plant materials
Haploid gametophytes (male and female), females containing carposporophytes and free-living diploid sporophytes of B. moritziana were obtained from the John A. West culture collection (isolate 2746). An asexual isolate that recycles apomictic tetraspores was also grown (isolate 3568-2). Female gametophytes that had several carposporophytes referred to as the carposporophyte was also used. The thalli were maintained in unialgal cultures in IMR medium (West et al. 2001, Klochkova et al. 2006) at 20°C in a 16 h light : 8 h dark cycle with illumination of >20 μmol photons m-2 s-1 provided by cool-white fluorescent lighting. Plants were harvested and cleaned with IMR medium. After removal of moisture from the samples, they were kept at -70°C until use.
RNA extraction and cDNA synthesis
Total RNA were isolated from the different life history stages using the RNeasy Plant RNA extraction kit (Qiagen, Gaithersburgh, MD, USA) according to manufacturer’s protocol. RNA isolation was performed three times for the plants of each life history stages. Genomic DNA contamination was eliminated by on-column digestion with RNase-free DNase I. RNA quality and concentration were determined using gel electrophoresis and spectrophotometrically. RNA samples with 260/280 ratio from 1.9 to 2.1, 260/230 ratio from 2.0-2.5 were used for further analysis. First-strand cDNA was synthesized from 3 μg of total RNA using the iScriptTM cDNA synthesis kit (Bio- Rad Laboratories, Hercules, CA, USA) in a final reaction volume of 20 μL with random primers following the manufacturer’s protocol. All cDNAs were diluted 10 times with autoclaved milli-Q water and stored at -20°C until use.
HKG selection and primers design
Potential HKGs which are popularly used for RT-qPCR analysis were selected from a B. mortiziana transcriptome database (http://genebank.kongju.ac.kr). Eleven genes were selected: elongation factor 1-alpha (EF-1α), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), beta-actin (ACTB), polyubiquitin (PUB), 30S ribosomal protein (30S), 60S ribosomal protein (60S), beta-tubulin (TUBB), alpha-tubulin (TUBA), translation initiation factor (TIF), ubiquitin-conjugating enzyme (UbCE), and isocitrate dehydrogenase (IcdH) (Table 1). Primers for selected genes were designed considering melting temperature and products size using the web-based Primer3 program provided by Bioneer (http://web.bioneer.co.kr/cgi-bin/primer/primer3.cgi).
qPCR and statistical analysis
Real-time PCR was performed using iQ SYBR Green supermix (Bio-Rad Laboratories) in a real-time PCR detection system (Bio-Rad CFX96; Bio-Rad Laboratories). cDNA samples were amplified in triplicate from the same RNA preparation. The final reaction volume was 20 μL and included 5 μL of diluted cDNA, 10 μL of iQ SYBR Green supermix and 250 nM of each primers (Table 1). The reaction protocol was first set at 95°C for 3 min and then 40 cycles of 95°C for 15 s and 60°C for 20 s. For each candidate HKG, the specificity of primers and amplified single product was verified by analyzing the melting curves and gel electrophoresis. The expression of each gene was determined from Cq values (synonym, Ct value) that corresponds to a number of cycles necessary for the PCR amplification to reach a threshold. Cq variability data in HKG comparisons were expressed as whiskers-box plot for all pooled samples. Selected HKGs were tested for the stability of their expression in different life history stages following a normalization strategy using two different algorithms: geNorm using the geNorm plugin of qbase + program (Vandesompele et al. 2002, Hellemans et al. 2007) and Normfinder (Andersen et al. 2004). The expression stability of all candidate HKGs was evaluated by comparing the differences of observed Cq values across life-history stages. Pairwise variation (Vn/n+1) was calculated using geNorm. The pairwise variation value (V) shows the optimal number of internal control genes required for accurate normalization (Andersen et al. 2004). If the value (V) is lower than 0.15 it is regarded as reliable gene set.
RESULTS
The haploid gametophytes and diploid sporophytes of B. moritziana are morphologically identical (isomorphic) in thallus structure except for their reproductive structure (Fig. 1). The basal for selecting HKGs in this study was that they are continuously expressed throughout the entire life cycle of B. moritziana, possess different functions from target genes, and have been used as internal controls in other gene expression studies. We searched our Bostrychia transcriptome library data, and selected 11 HKGs (Table 1). The primers of all the selected HKGs were designed for real-time PCR based on our transcriptome library data (Table 1).
The stability of expression of HKGs in Bostrychia was investigated with Cq values in real-time PCR. Melting curves for each gene were analyzed, all the genes were confirmed to have a single amplification curve. All genes did not generate non-specific products (Appendix 1). The Cq values of each gene showed a value between 19-36 (standard deviation 0.5-2) except for GAPDH, 30S, and IcdH (Fig. 2). We found that PUB showed the highest expression level with a mean Cq value of 20. Conversely, 60S showed the lowest expression level with a mean Cq value of 35 (Fig. 2).
The stability of the gene expression in the different life history stages was analyzed (Fig. 3). Analysis using ge- Norm showed that GAPDH was the most stable HKG with an M value of 0.29, followed by 30S (0.32) and IcdH (0.35) (Fig. 3A). HKGs are regarded as reliable when M values are lower than 0.5 (Vandesompele et al. 2002). Analysis using NormFinder also showed that GAPDH had the most candistable expression levels between the life history stages followed by 30S (Fig. 3B).
In B. moritziana, our V value (V2/3 was 0.124) suggests that a combination of two HKGs is good enough to yield accurate normalization (Fig. 4). These results indicate that the combination of GAPDH and IcdH gives the best gene set for internal control in comparisons between the life history stages of B. moritziana.
DISCUSSION
The temporal, developmental and physiological pattern of gene expression is a topic that is now amiable to scientific study. Changes in gene expression, even fairly subtle changes, could have significant effects on the physiology of the organism of study (Kim et al. 2016). To accurately quantify these gene expression changes a stable comparison is needed, to account for variation in RNA quantification and other experimental error. Accurate normalization is an absolute prerequisite for accurate measurements of gene expression changes, yet the determination of these HKGs (reference genes) could be a circular problem as the expression level of these reference genes also needs to be normalized (Hibbeler et al. 2008). The search therefore for genes that are stable over multiple conditions or life history stages is critical and several methods have been established (geNorm in Vandesompele et al. 2002, NormFinder in Andersen et al. 2004). Often normalization procedures to validate candidate genes are not done and HKGs are used that have been used in other organisms or studies, without any validation of their stability in the specific study (Chan et al. 2014). While validation methods are essential to have a stable platform from which changes can be measured between samples, in red algae, few studies have validated HKGs before measuring differences in target genes (e.g., Asamizu et al. 2003 used 60S L3 protein). This could lead to erroneous conclusions on the changes in expression level in target genes. Our results showed that, at least between life history stages of B. moritziana in constant culture conditions, the combination of GAPDH and 30S was the best gene set for this internal control. More than half of the tested HKGs including ACTB and 60S protein were no good for internal controls.
Tests for HKG stability in life history stages of the red alga has recently been performed in P. yezoensis using seven candidate genes (Wu et al. 2013). Their study also showed significant differences in expression of candidate genes. Their conclusion was that GAPDH was the most stable candidate gene. Our study using 11 genes also showed GAPDH as a useful normalization gene (best with geNorm and NormFinder) in comparisons between the life history stages of B. mortiziana. Kowalczyk et al. (2014) performed similar validation experiments for HKGs, using C. crispus, monitored the expression of 14 genes in conditions using two sets of samples: cultured samples stressed in laboratory conditions and C. crispus collected on the shore and stressed in situ. Their results showed that translation initiation factor 4A-2 and eukaryotic polypeptide chain release factor 3 were the best normalization genes and GAPDH was one of the worst genes for the internal control in stress-related experiment. Both C. crispus and B. moritziana are members of the class Florideophyceae (subclass Rhodymeniophycidae). These results show the importance of testing reference genes according to the nature of the experiments rather than the evolutionary background of the organism. Our results showed that IcdH and 30S are also pertinent candidates for HKGs in life history study of B. moritziana but previous studies showed that these genes are not stable in other red algae, P. yezoensis and C. crispus (Wu et al. 2013, Kowalczyk et al. 2014).
It is clear from our analysis, and several other studies, that housekeeping genes are not equally expressed and the expression of ‘standard’ housekeeping genes must be tested for each organism under study. Although RT-qPCR is expensive and time consuming a survey of potential candidate genes at the start of a study would be essential in reducing uncertainty, increasing confidence and may lessen the need for additional experiments to confirm any results.