Skip to main content

METHODS article

Front. Physiol., 07 May 2020
Sec. Invertebrate Physiology

Identification and Validation of Reference Genes for Quantitative Gene Expression Analysis in Ophraella communa

\r\nYan ZhangYan Zhang1Jiqiang ChenJiqiang Chen1Guangmei ChenGuangmei Chen1Chao MaChao Ma1Hongsong Chen,Hongsong Chen1,2Xuyuan Gao,Xuyuan Gao1,2Zhenqi TianZhenqi Tian1Shaowei CuiShaowei Cui1Zhenya TianZhenya Tian1Jianying GuoJianying Guo1Fanghao WanFanghao Wan1Zhongshi Zhou*Zhongshi Zhou1*
  • 1State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing, China
  • 2Guangxi Key Laboratory of Biology for Crop Diseases and Insect Pests, Institute of Plant Protection, Guangxi Academy of Agricultural Sciences, Nanning, China

Ophraella communa is an effective bio-control agent of the invasive common weed. By now, the reference genes in O. communa have not yet been screened and validated. The aim of this study was to screen for the most stable reference genes in different backgrounds, such as different developmental stages, sexes, tissues, and male reproductive system with different body sizes. We selected 12 common housekeeping genes involved in different biological processes, including GAPDH, ACT1, ACT2, ARF1, ARF4, SDH, βTUBC, RPL4, RPL19, RPS18, EF1α, and COX as the candidate reference genes. To analyze the stability of the candidate reference genes, we first used three dedicated algorithms, GeNorm, NormFinder, and BestKeeper, and further comprehensive ranking was provided by ReFinder. The results showed that RPL19 and RPL4 exhibited the least variation in different developmental stages/sexes and in male reproductive systems with different body sizes. COX proved to be most suitable for normalizing the gene expression levels in different tissues, and coincidentally, RPL19 was also found to be second in terms of stability in this study. To the best of our knowledge, this is the first study to identify suitable reference genes for analyzing gene expression in O. communa; thus, this study would lay the foundation for future research on the molecular physiology and biochemistry of O. communa and other insects.

Introduction

Ophraella communa LeSage (Coleoptera: Chrysomelidae) originated in North America (LeSage, 2012). It is a biological control agent of the invasive weed Ambrosia artemisiifolia L. (common ragweed) (Yamazaki et al., 2000; Zhou et al., 2010a, 2014; Gerber et al., 2011; Guo et al., 2011). For a biological control agent, reproductive regulation is a very important issue in order to facilitate its mass rearing. Gratifyingly, we have found that the O. communa males with a large body size can facilitate a female’s fecundity (unpublished data). Therefore, how to regulate the reproductive potential of the male beetle will be focused on in a further study. It is necessary to assess gene function in the beetle; quantification of gene expression is required as well.

In recent years, the ecology and physiology of O. communa has been studied extensively (Zhou et al., 2010a,b; Clark et al., 2014; Kim and Lee, 2018; Zhao et al., 2018). However, the molecular mechanisms of its biology remain largely unknown (Ma et al., 2019a,b). With the aim of investigating the genetic basis of its biology and physiology, the male transcriptome and proteome of O. communa have been sequenced (unpublished data). Remarkably, we found that the male’s body size is an important determinant in the evolutionary process and in the speciation of O. communa, as it affects female mate choice. Then, the large males could further improve the population fitness of their offspring (unpublished data). Currently, the mechanisms for higher fitness in large males are also being studied, so it is necessary to assess gene function in O. communa; quantification of gene expression is required as well.

Reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR), also referred to as qPCR, owing to its excellent repeatability, high sensitivity, accurate quantification, wide applicability, and easy operation, is the most common and sensitive technique for the detection and quantification of gene expression in different experimental samples, such as developmental stages, tissues, and other processing conditions (Bustin, 2002). It is especially useful to detect low-abundance mRNAs in limited samples (Bustin et al., 2005; Valasek and Repa, 2005; Chen et al., 2011; Hellemans and Vandesompele, 2014). However, the results of qPCR are often limited by factors such as the amount of the initial sample, cDNA synthesis efficiency, and primer amplification efficiency, which may result in inaccurate qPCR results (Bustin et al., 2005). Therefore, this approach requires normalization.

To date, internal control genes, also known as housekeeping genes or reference genes, are genes whose expression levels are stable regardless of the cell type and the experimental conditions. It is the most popular tool to normalize variability (Thellin et al., 1999; Vandesompele et al., 2002; Andersen et al., 2004; Huggett et al., 2005; Hellemans and Vandesompele, 2014). Many genes such as Actin, GAPDH, and rRNAs (18s or 28s) that are involved in fundamental biological processes are frequently used as reference genes (Keertan et al., 2004; de Kok et al., 2005). Nevertheless, it has become clear that no single reference gene retains a constant expression in all tissues, developmental stages, or experimental conditions (Silver et al., 2006; Lourenço et al., 2008; Lu et al., 2013, 2018a; Rodrigues et al., 2014; Tan et al., 2015; Dai et al., 2017; Luo et al., 2018). Therefore, the traditional genes are not applicable to all species and conditions (Shu et al., 2018; Zhou et al., 2018; Qu et al., 2019). Blindly choosing “classical” reference genes such as β-actin, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), or ribosomal 18S (18S) is unreliable and inaccurate. The reliability of qPCR largely depends on the stability of the reference genes. Previous studies have shown that it is crucial to select and validate the reference gene for data quantification (Lu et al., 2018b; Shakeel et al., 2018; Shu et al., 2018). Hence, it is extremely important to verify the stability of expression of the intended control gene under different experimental conditions. To our knowledge, research on the reference gene of O. communa is absent, which makes the study of the mechanisms of its physiological and biochemical behavior difficult.

In this study, we aim to identify the solidly expressed reference genes for RT-qPCR investigation under different conditions including developmental stages, sexes, tissues, and male reproductive system with different body sizes of O. communa. Based on prior experimental reports regarding reference genes in coleoptera and other insects, 12 candidate reference genes including glyceraldehyde-3-phosphate dehydrogenase (GAPDH), actin 1 (ACT1), actin 2 (ACT2), succinate dehydrogenase (SDH), ADP-ribosylation factor 1 (ARF1), ADP-ribosylation factor 4 (ARF4), cytochrome oxidase (COX), β-tubulin C (βTUBC), elongation factor-1α (EF1α), ribosomal protein L4 (RPL4), ribosomal protein L19 (RPL19), and ribosomal protein S18 (RPS18), which are involved in various biological processes, were selected from the male transcriptome of O. communa. The results of the qPCR were analyzed using three statistical algorithms, GeNorm (Vandesompele et al., 2002), NormFinder (Andersen et al., 2004), and BestKeeper (Pfaffl et al., 2004). Finally, a comprehensive ranking of the stability of these candidates was given by ReFinder (Yang et al., 2015a). In addition, the expression profile of angiotensin-converting enzyme (ACE) was tested to verify our result. This study provides a valuable platform for gene expression and future functional research and serves as an example for similar studies in other insects.

Materials and Methods

Insect

Ophraella communa adults were collected from the fields in Laibin City in Guangxi Province in the summer of June 2017 and were maintained on common ragweed plants in the laboratory at the Chinese Academy of Agricultural Sciences (BJ, CAAS) for three generations. All the insects used for sampling were reared in an environmental chamber (PRX-450D, Ningbo Haishu Safe Experimental Equipment Co., Ltd., Zhejiang, China) at 27 ± 1°C, with 50 ± 10.0% relative humidity and a 14:10-hour photoperiod provided by artificial lights and with the robust common ragweed plants used as food.

Experimental Sample in This Study

Samples of Developmental Stages and Sexes

The adults deposited their eggs in batches of about 35 on the underside of common ragweed leaves. The number of sampled individuals for each replicate across the different developmental stages and sexes of O. communa included 150 eggs, 30 first instar nymphs, 20 second instar nymphs, eight third instar nymphs, eight pupa (third pupa), five sexually immature (1 day old) male and female adults, and five sexually mature (8 days old) male and female adults. Three independent biological replicates were collected for each sample.

Tissue

Male adults (8 days old) of O. communa were dissected for samples of five different tissues, including head, testis, male accessory glands (MAG), fat body, and gut. Three independent biological replicates were collected with each group of 12 insects.

Male Reproductive System With a Different Body Size

The testis and accessory glands of large and small males, respectively, were dissected for these samples. Three independent biological replicates were collected with each group of 12 insects.

All the above samples were collected in 1.5-ml tubes, immediately frozen in liquid nitrogen, and stored at −80°C until use.

RNA Extraction and cDNA Synthesis

Total RNA from all 48 samples was extracted using the Trizol reagent (Invitrogen, United States) according to the manufacturer’s protocol. The RNA integrity was tested by 1% agarose gel electrophoresis and quantified using NanoPhotometerTM P330 (Implen, Munich, Germany), with A260/A280 ratios ranging from 1.8 to 2.1 for all RNA samples. One microgram of total RNA was used to synthesize cDNAs using the TransScript One-Step gDNA Removal and cDNA Synthesis SuperMix (Transgen Biotech, Beijing, China) according to the manufacturer’s instructions.

Reference Gene Selection and Primer Design

A total of 12 candidate reference genes that are most frequently used in qPCR investigations were selected from the O. communa transcriptome data (Table 1). The primers were designed for gene cloning and qPCR using the Primer five programs. All the candidate reference genes were PCR-amplified from the bettle’s cDNA using the corresponding primers (Table 1) and the pEASY-T3 Simple Cloning Kit (TransGen, Beijing, China). Six independent subclones were sequenced and accurate nucleotide sequences were obtained after the sequence analysis. All the verified reference gene sequences were submitted to the GenBank database.

TABLE 1
www.frontiersin.org

Table 1. qPCR primers for candidate reference genes and characteristics of PCR amplification in O. communa.

Quantitative Real-Time PCR Analysis

Quantitative PCR reactions were performed using an ABI 7500 PCR Detection System (Applied Biosystems, United States). The qPCR reactions consisted of 10 μl of 2 × TransStartTM Green qPCR SuperMix, 1 μl of the cDNA template, 0.4 μl of passive reference dye, and 0.4 μl each of the forward and reverse primers (10 μM) in a 20-μl final volume. The RT-qPCR procedure was as follows: an initial cycle at 94°C (30 s), followed by 40 cycles of 5 s at 94°C and 34 s at 60°C. Following the RT-qPCR procedure, the consistency and the specificity of the PCR-amplified products were subjected to melting curve analysis for all reactions. The amplification efficiencies (E) and correlation coefficients (R2) were determined for each gene using the standard curves with a 5-fold dilution series of the template (1/2, 1/10, 1/50, 1/250, and 1/1,250), where R2 is the slope of the standard curve. E was calculated according to the equation E = (10[−1/slope] −1) × 100 (Michael and Pfaffl, 2001). Each reaction was performed in triplicate.

Determination of Stability of Expression of the Candidate Reference Gene

The stability of the 12 candidate reference genes was evaluated using three different statistical programs, GeNorm (Vandesompele et al., 2002), NormFinder (Andersen et al., 2004), and BestKeeper (Pfaffl et al., 2004) under experimental conditions, and the instructions of all the three programs were consulted for the data arrangement and the subsequent analysis. Finally, the comprehensive ranking of each condition was obtained according to ReFinder. In addition, the number of reference genes for normalizing gene expression was decided by the pairwise variation (Vn/n+1) which was performed using the GeNorm program. Universally, Vn/n+1 less than the threshold value of 0.15 indicates that the most suitable number is n, and there is no need to introduce the n + 1 reference gene for normalization (Vandesompele et al., 2002).

Validation of Reference Genes

Angiotensin-converting enzyme (ACE), a zinc-dependent peptidase (Isaac et al., 1999, 2007; Macours and Biol, 2004), was reported as a seminal fluid protein in Drosophila melanogaster and Tribolium castaneum (Rylett et al., 2007; Xu et al., 2013). The ACE gene (GenBank number: MN641122) cloning primers used were as follows: forward (5′-AATCCATGCCTACGTACGGT AT-3′) and reverse (5′-CTTCCCATCTCCACAGGTCA-3′). The normalization of ACE expression in different tissues was conducted using the most stable gene combinations (COX/RPL19) as determined by geNorm, the most single stable gene (RPL19), and the least stable gene (SDH)as determined by ReFinder. The relative quantification of the ACE gene was calculated using the 2–ΔΔCT method (Schmittgen and Livak, 2008). One-way ANOVA followed by Tukey’s least significant difference test was used to determine the significance of the ACE expression levels between different tissues (SAS 9.4).

Results

Acquisition of Candidate Reference Gene Sequences

The sequences of the 12 candidate gene, including RPL4, RPS18, ACT1, ACT2, ARF1, ARF4, SDH, GAPDH, EF1α, βTUBC, COX, and RPL19, were obtained from the transcriptome database of the testis and the accessory glands of O. communa. The sequences were verified by qPCR. All amplicons were sequenced and confirmed to exhibit 99–100% identity with the corresponding transcriptome sequence. The sequence information is listed in Supplementary Table S1.

Specificity and Efficiency of RT-qPCR Primers

The accuracy of the RT-qPCR primers was confirmed by the presence of a single peak in the melting curve analyses and the absence of primer dimers, both of which also indicate the specificity of each primer pair (Figure 1). The amplification efficiency (E) of each primer (shown in Table 2) ranged from 97.05 to 105.63%, with associated R2 values of 0.996–1.000 (Table 1).

FIGURE 1
www.frontiersin.org

Figure 1. Melting curves of the 12 candidate reference genes in O. communa.

TABLE 2
www.frontiersin.org

Table 2. Stability ranking of the 12 candidate reference genes’ expression in O. communa under different experimental conditions calculated by the three different analytical tools, geNorm, Normfinder, and BestKeeper, respectively.

Expression Profiles of Candidate Reference Genes

The expression profiles of the 12 candidate genes were analyzed by qPCR, and the Ct values exhibited different levels of transcript abundance (Figure 2). For the different developmental stages and sexes, the mean Ct values of the 12 candidate reference genes ranged from 18.01 to 26.98 cycles. In different tissues, the Ct values ranged from 18.13 to 28.83. Even in males of different sizes, it varied from 20.58 to 28.76.

FIGURE 2
www.frontiersin.org

Figure 2. The expression profiles of the 12 candidate reference genes are shown as Ct values under three experimental conditions for O. communa. (A) Different developmental stages and sexes. (B) Different tissues. (C) Male reproductive systems with different body sizes. (D) Total Ct values exhibited different levels of transcript abundance in three biological samples.

Analysis of Stability of Candidate Reference Genes

For different developmental stages and sexes—based on the analysis using GeNorm, NormFinder, and BestKeeper—the expression of RPL19 varied the least, while GAPDH tended to fluctuate the most. According to ReFinder parameters, the comprehensive ranking in these samples was as follows: RPL19 > RPL4 > ACT2 > ACT1 > RPS18 > SDH > EF1α> COX > ARF1 > βTUBC > ARF4 > GAPDH (Figure 3A). For the analysis of different tissue types, the stability of RPL4 and RPL19 were similarly excellent, while SDH and GAPDH were consistently recommended as the least stable by GeNorm, NormFinder, and BestKeeper. The comprehensive ranking was COX > RPL19 > RPL4 = ARF1 = RPS18 > βTUBC > ACT2 > ARF4 > EF1α> ACT1 > GAPDH > SDH (Figure 3B). For the male reproductive system with different body sizes, the overall variability of the reference genes was low except for βTUBC and SDH. The overall ReFinder stability ranking (from most to least stable) was RPL4 > RPL19 > ARF4 > ACT2 > GAPDH > ACT1 > COX > EF1α> RPS18 > ARF1 > SDH > βTUBC (Figure 3C). The stability of each sample’s stability was calculated using GeNorm, Normfinder, and BestKeeper, and all the specific values are shown in Table 2.

FIGURE 3
www.frontiersin.org

Figure 3. Expression stability ranking of the 12 candidate reference genes analyzed using ReFinder under three experimental conditions. (A) Different developmental stages and sexes. (B) Different tissue. (C) Male accessory glands from large and small males. The most stable genes are on the right and the least stable genes are on the left.

Selection of the Optimal Number of Reference Genes to Normalize the Gene Expression

We used GeNorm to compute the variation value to ensure the optimal number of reference genes. As shown in Figure 4, the V2/3 was below 0.15 across the different experimental conditions, which indicated that two reference genes were reliable for an accurate normalization of gene expression.

FIGURE 4
www.frontiersin.org

Figure 4. Pairwise variation in the three groups calculated by geNorm is required for accurate normalization of gene expression. The threshold used was 0.15. The arrow indicates that the value of V2/V3 is less than 0.15.

Validation of Selected Reference Genes in O. communa

The ACE expression level was used to verify the validity of the reference genes between the different tissues. The best paired, best single, and the worst single reference genes—COX/RPL19, RPL19, and SDH, respectively, were used as internal controls. The ACE expression patterns were similar in that the highest expression of ACE was in the male reproductive tissues (testis and MAG) when normalized by the paired COX/RPL19 and single RPL19 (Figure 5). However, when we used SDH to normalize the expression levels, the highest expression level of ACE was in the heads and significantly higher than in the other tissues (Figure 5). This conclusion was completely inconsistent with the results seen through the normalization with RPL19 or COX/RPL19. This indicated that the accuracy of the RT-qPCR results may be influenced by the instability of the reference gene, which might even lead to an inaccurate or wrong conclusion.

FIGURE 5
www.frontiersin.org

Figure 5. Relative quantification of ACE genes in different tissues. The expression levels of ACE genes were normalized using the best reference paired COX/RPL19, the single best stable gene RPL19, and the reference gene SDH that changed the most in tissues. Values are means ± SE. Different letters indicate significant differences in gene expression among the different tissues of O. communa (P < 0.05).

Discussion

Ophraella communa is a very valuable and effective biological control agent of Ambrosia artemisiifolia (Guo et al., 2011). Unfortunately, research on the molecular mechanisms of O. communa is scarce, making the study of its physiological and biochemical behavior difficult. Recently, the male transcriptome of O. communa was sequenced using Illumina high-throughput sequencing technologies (unpublished data), and further research on functional genomics and related gene expression studies will be initiated soon.

The quantitative accuracy of gene expression relies on the reference gene (Lu et al., 2018b). In our study, based on the study of reference genes in Coleoptera, 12 candidate reference genes, including traditional genes and novel genes, were selected. These genes have been evaluated as the most stable genes in different insects: RP18, RP4, AFR1, and ARF4 in Leptinotarsa decemlineata (Shi et al., 2013), β-actin and EF1a in Diabrotica virgifera virgifera (Rodrigues et al., 2014), Actin and RPS18 in T. castaneum (Lord et al., 2010; Toutges et al., 2010; Sang et al., 2015), EF1α, ACT1, and RPL19 in Colaphellus bowringi (Tan et al., 2015), SDHA, ACT and TUB in Galeruca daurica (Tan et al., 2017), GADPH in Hippodamia convergensin (Pan et al., 2015), and COX in Agasicles hygrophila. Besides that, these candidate genes have also been recommended for normalization gene expression in others (Lu et al., 2018b; Shakeel et al., 2018). Three different methods, GeNorm (Vandesompele et al., 2002), NormFinder (Andersen et al., 2004), and BestKeeper (Pfaffl et al., 2004), were used to identify the stability of the 12 candidate genes across the condition of developmental stages and sexes, different tissues, and male reproductive system with different body sizes. However, because each method has its own strengths and appropriate application conditions, the expression stability obtained by different algorithms was different (Table 2). Therefore, another analysis tool, ReFinder (Yang et al., 2015a), was used to integrate the variation in the results between them and to define an overall final ranking (Xie et al., 2012).

Our current research focuses on the molecular mechanisms of the male fitness differences relative to different body sizes. The results of this study displayed that the reference gene transcript levels vary with the above experimental conditions (Figure 2). For the sample of different developmental stages and sexes, contrary to the results of other studies in Coleoptera insects, GAPDH, EF1α, and TUB performed the best stability under this condition (Rodrigues et al., 2014; Pan et al., 2015; Yang et al., 2018). Our results demonstrated that RPL19 and RPL4 showed the most stable gene expression, while GAPDH, EF1α, and βTUBC showed the least stability (Figure 3A). All four algorithms obtained the same two top-ranked genes with few differences in their respective orders. Similar to our results, in L. decemlineata, RPL4 and RPS18 exhibited excellent stability in different developmental stages and sexes. Ribosomal protein (RP) also presented superior stability in Coccinella septempunctata across this condition (Shi et al., 2013; Tan et al., 2015; Yang et al., 2016). For different tissues—consistent with the results in A. hygrophila—a novel reference gene, COX, that encodes for a key enzyme (Liu et al., 2017) and that responds to a variety of metabolic states showed the least variation (Figure 3B). In addition, like the ranking for developmental stages and sexes, the RPs, including RPL19, RPLl4, and RPS18, displayed high stability (Figure 3B), while SDH and GADPH were identified as the least stable genes. For the male reproductive system with different body sizes, the variation of each candidate gene is relatively small compared with the first two conditions (Figure 2C), and a similar ranking trend with previous condition was observed in this sample. RPL19 and RPL4 performed the least variation. To our knowledge, studies on the reference gene for body size are scarce. Our results provide a reference for the gene expression in other insects with different body sizes. Integrating the excellent stable reference genes under three conditions, RPL19 can be used to normalize the gene expression. Studies in C. bowringi similarly concluded that the RPL19 presented excellent stability across developmental stage, sex, and tissue (Tan et al., 2015). In D. virgifera virgifera, EF1a and Actin were respectively identified as the most stable gene among different tissues and developmental stages (Rodrigues et al., 2014). However, our results were almost the opposite. Therefore, the choice of the control gene for O. communa is of utmost importance.

Coincidentally, our results verified that RP gene expression was extremely stable under each experimental condition. Generally, ribosomal protein genes are expressed in all types of cells for the synthesis of new ribosomes, and RPs are among the most highly conserved proteins across all life forms (Hsiao et al., 2001). Recently, researchers found that even in T. castaneum, most RP gene expression was maintained stably following a fungal challenge (Lord et al., 2010; Lu et al., 2018b). In addition, besides Coleoptera insects, ribosomal genes showed high stability in diverse abiotic and biotic conditions in other organisms (Jonge et al., 2007; Popovici et al., 2009; Tong et al., 2009; Zhong et al., 2013; Yang et al., 2015b, 2016; Ma et al., 2016; Yan et al., 2016; Lu et al., 2018a; Luo et al., 2018). Our screening results also demonstrated that the RP gene is transcriptionally conserved under different physiological conditions. However, under some conditions, the expression levels of RP-encoding genes may be unstable. For example, RPS3 was the least stable gene in Spodoptera litura under different populations, exposed to different temperatures, or treated by different insecticides (Lu et al., 2013). Moreover, other classic references are GADPH, TUB, and Actin which have poor stability in O. communa (Figure 3). Therefore, it is very unwise to choose a reference gene blindly regardless of the conditions and species.

Furthermore, to certify the reliability of the reference gene in O. communa, the relative transcriptional level of the ACE gene was determined in different tissues using qPCR. GeNorm recommended that the use of two reference genes was optimal for normalization under each experimental condition (Figure 4). Therefore, ACE expression was normalized by the single stable RPL19, the pairwise stable RPL19 and COX, and the least stable SDH alone. Consistent with a previous study (Xu et al., 2013), when RPL19/COX and RPL19 alone were used as internal controls, ACE expression was found to be abundant in the reproductive organs, including the testis and the MAG. In contrast, the head and the testis showed a high ACE expression when normalized using SDH (Figure 5). This implies that different controls may lead to different conclusions. Therefore, we strongly recommend the validation of the stability of the reference gene expression prior to performing RT-qPCR experiments. Our study found that RPL19 and RPL4 showed excellent stability among the candidate reference genes for all three experimental conditions. This provides a fundamental premise to analyze qPCR data in O. communa and to further clarify the physiological and biochemical mechanisms in this important natural ragweed enemy.

Conclusion

In summary, we identified several stable housekeeping genes that are suitable as reference genes in O. communa. This is the first study to report the reference genes and to establish a standardized RT-qPCR analysis protocol in O. communa, which will contribute to the development of molecular research on this vital insect.

Data Availability Statement

The datasets generated for this study can be found in the NCBI using accession MN641110, MN641111, MN641112, MN641113, MN641114, MN641115, MN641116, MN641117, MN641118, MN641119, MN641121, MN641120, and MN641122.

Ethics Statement

The animal study was reviewed and approved by Ethics of Animal Experiment of Institute of Plant Protection, Chinese Academy of Agricultural Sciences.

Author Contributions

ZZ designed the study and revised the manuscript. YZ performed the research and wrote the manuscript. JC applied software for data analysis. GC searched for the sequence of candidate reference genes. CM, XG, ZQT, SC, and ZYT participated in the sample collection and data sorting. HC, JG, and FW guided in the data analysis.

Funding

This study was supported by the National Natural Science Foundation of China (Nos. 31672089 and 31972340).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We would like to thank Editage for polishing the English.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2020.00355/full#supplementary-material

References

Andersen, C. L., Jensen, J. L., and Ørntoft, T. F. (2004). 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 Res. 64, 5245–5250. doi: 10.1158/0008-5472.Can-04-0496

PubMed Abstract | CrossRef Full Text | Google Scholar

Bustin, S. A. (2002). Quantification of mRNA using real-time reverse transcription PCR (RT-PCR): trends and problems. J. Mol. Endocrinol. 29, 23–39. doi: 10.1677/jme.0.0290023

PubMed Abstract | CrossRef Full Text | Google Scholar

Bustin, S. A., Benes, V., Nolan, T., and Pfaffl, M. W. (2005). Quantitative real-time RT-PCR - a perspective. J. Mol. Endocrinol. 34, 597–601. doi: 10.1677/jme.1.01755

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C., Tan, R., Wong, L., Fekete, R., and Halsey, J. (2011). Quantitation of microRNAs by real-Time RT-qPCR. Methods Mol. Biol. 687, 113–134. doi: 10.1007/978-1-60761-944-4_8

PubMed Abstract | CrossRef Full Text | Google Scholar

Clark, S., Scavarda, G., Chersi, C., Massobrio, V., and Bosio, G. (2014). Spread of the ragweed leaf beetle, Ophraella communa LeSage, 1986 (Coleoptera Chrysomelidae), in piedmont region (northwestern Italy). Boll. Soc. Ent. Ital. 146, 17–30. doi: 10.4081/BollettinoSEI.2014.17

CrossRef Full Text | Google Scholar

Dai, T. M., Lu, Z. C., Liu, W. X., and Wan, F. H. (2017). Selection and validation of reference genes for qRT-PCR analysis during biological invasions: the thermal adaptability of Bemisia tabaci MED. PLoS One 12:e0173821. doi: 10.1371/journal.pone.0173821

PubMed Abstract | CrossRef Full Text | Google Scholar

de Kok, J. B., Roelofs, R. W., Giesendorf, B. A., Pennings, J. L., Waas, E. T., Feuth, T., et al. (2005). Normalization of gene expression measurements in tumor tissues: comparison of 13 endogenous control genes. Lab. Invest. 85, 154–159. doi: 10.1038/labinvest.3700208

PubMed Abstract | CrossRef Full Text | Google Scholar

Gerber, E., Schaffner, U., Gassmann, A., Hinz, H. L., Seier, M., and Müller-Schärer, H. (2011). Prospects for biological control of Ambrosia artemisiifolia in europe: learning from the past. Weed. Res. 51, 559–573. doi: 10.1111/j.1365-3180.2011.00879.x

CrossRef Full Text | Google Scholar

Guo, J. Y., Zhou, Z. S., Zheng, X. W., Chen, H. S., Wan, F. H., and Luo, Y. H. (2011). Control efficiency of leaf beetle, Ophraella communa, on the invasive common ragweed, Ambrosia artemisiifolia, at different growing stages. Biocontrol. Sci. Techn. 21, 1049–1063. doi: 10.1080/09583157.2011.603823

CrossRef Full Text | Google Scholar

Hellemans, J., and Vandesompele, J. (2014). Selection of reliable reference genes for RT-qPCR analysis. Methods Mol. Biol. 1160, 19–26. doi: 10.1007/978-1-4939-0733-5_3

PubMed Abstract | CrossRef Full Text | Google Scholar

Hsiao, L. L., Dangond, F., Yoshida, T., Hong, R., Jensen, R. V., Misra, J., et al. (2001). A compendium of gene expression in normal human tissues. Physiol. Genomics 7, 97–104. doi: 10.1152/physiolgenomics.00040.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

Huggett, J., Dheda, K., Bustin, S., and Zumla, A. (2005). Real-time RT-PCR normalisation; strategies and considerations. Genes Immun. 6, 279–284. doi: 10.1038/sj.gene.6364190

PubMed Abstract | CrossRef Full Text | Google Scholar

Isaac, R. E., Ekbote, U., Coates, D., and Shirras, A. D. (1999). Insect angiotensin-converting enzyme: a processing enzyme with broad substrate specificity and a role in reproduction. Ann. N. Y. Acad. Sci. 897, 342–347. doi: 10.1111/j.1749-6632.1999.tb07904.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Isaac, R. E., Lamango, N. S., Ekbote, U., Taylor, C. A., Hurst, D., Weaver, R. J., et al. (2007). Angiotensin-converting enzyme as a target for the development of novel insect growth regulators. Peptides 28, 153–162. doi: 10.1016/j.peptides.2006.08.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Jonge, H. J. M. D., Fehrmann, R. S. N., Bont, E. S. J. M., Hofstra, R. M. W., Gerbens, F., Kamps, W. A., et al. (2007). Evidence based selection of housekeeping genes. PLoS One 2:e898. doi: 10.1371/journal.pone.0000898

PubMed Abstract | CrossRef Full Text | Google Scholar

Keertan, D., Huggett, J. F., Bustin, S. A., Johnson, M. A., Graham, R., and Alimuddin, Z. (2004). Validation of housekeeping genes for normalizing RNA expression in real-time PCR. Biotechniques 37, 118–119. doi: 10.1016/j.biosystems.2004.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, H. G., and Lee, D. H. (2018). Review of the biology and ecology of a ragweed leaf beetle, Ophraella communa (Coleoptera: Chrysomelidae), which is a biological control agent of an invasive common ragweed, Ambrosia artemisiifolia (Asterales: Asteraceae). Biocontrol. Sci. Techn. 29, 185–200. doi: 10.1080/09583157.2018.1540032

CrossRef Full Text | Google Scholar

LeSage, L. (2012). A taxonomic monograph of the nearctic galerucine genus Ophraella Wilcox (Coleoptera: Chrysomelidae). Mem. Entomol. Soc. Can 118, 3–75. doi: 10.4039/entm118133fv

CrossRef Full Text | Google Scholar

Liu, J., Hiser, C., and Ferguson-Miller, S. (2017). Role of conformational change and K-path ligands in controlling cytochrome coxidase activity. Biochem. Soc. Trans. 45:BST20160138. doi: 10.1042/BST20160138

PubMed Abstract | CrossRef Full Text | Google Scholar

Lord, J. C., Hartzer, K., Toutges, M., and Oppert, B. (2010). Evaluation of quantitative PCR reference genes for gene expression studies in Tribolium castaneum after fungal challenge. J. Microbiol. Methods 80, 219–221. doi: 10.1016/j.mimet.2009.12.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Lourenço, A. P., Mackert, A., dos Santos Cristino, A., and Simões, Z. L. P. (2008). Validation of reference genes for gene expression studies in the honey bee. Apis mellifera, by quantitative real-time RT-PCR. Apidologie 39, 372–385. doi: 10.1051/apido:2008015

CrossRef Full Text | Google Scholar

Lu, J., Chen, S., Guo, M., Ye, C., Qiu, B., Wu, J., et al. (2018a). Selection and validation of reference genes for RT-qPCR analysis of the ladybird beetle Henosepilachna vigintioctomaculata. Front. Physiol. 9:1614. doi: 10.3389/fphys.2018.01614

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, J., Yang, C., Zhang, Y., and Pan, H. (2018b). Selection of reference genes for the normalization of RT-qPCR data in gene expression studies in insects: a systematic review. Front. Physiol. 9:1560. doi: 10.3389/fphys.2018.01560

CrossRef Full Text | Google Scholar

Lu, Y., Yuan, M., Gao, X., Kang, T., Zhan, S., Wan, H., et al. (2013). Identification and validation of reference genes for gene expression analysis using quantitative PCR in Spodoptera litura (Lepidoptera: Noctuidae). PLoS One 8:e68059. doi: 10.1371/journal.pone.0068059

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, J., Ma, C., Li, Z., Zhu, B. Q., Zhang, J., Lei, C. L., et al. (2018). Assessment of suitable reference genes for qRT-PCR analysis in Adelphocoris suturalis. J. Integr. Agr. 17, 2745–2757. doi: 10.1016/s2095-3119(18)61926-4

CrossRef Full Text | Google Scholar

Ma, C., Cui, S., Tian, Z., Zhang, Y., Chen, G., Gao, X., et al. (2019a). OcomCSP12, a chemosensory protein expressed specifically by ovary, mediates reproduction in Ophraella communa (Coleoptera: Chrysomelidae). Front. Physiol. 10:1290. doi: 10.3389/fphys.2019.01290

CrossRef Full Text | Google Scholar

Ma, C., Zhao, C., Cui, S., Zhang, Y., Chen, G., Chen, H., et al. (2019b). Identification of candidate chemosensory genes of Ophraella communa LeSage (Coleoptera: Chrysomelidae) based on antennal transcriptome analysis. Sci. Rep. 9, 1–14. doi: 10.1038/s41598-019-52149-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, K. S., Li, F., Liang, P. Z., Chen, X. W., Liu, Y., and Gao, X. W. (2016). Identification and validation of reference genes for the normalization of gene expression data in qRT-PCR analysis in Aphis gossypii (Hemiptera: Aphididae). J. Insect. Sci. 16:17. doi: 10.1093/jisesa/iew003

PubMed Abstract | CrossRef Full Text | Google Scholar

Macours, N., and Biol, K. H. J. I. B. M. (2004). Zinc-metalloproteases in insects: ACE and ECE. Insect Biochem. Molec. 34, 501–510. doi: 10.1016/j.ibmb.2004.03.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Michael, W., and Pfaffl, A. (2001). A new mathematical model for relative quantification in real-time RT–PCR. Nucleic Acids Res. 29:e45. doi: 10.1093/nar/29.9.e45

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, H., Yang, X., Siegfried, B. D., and Zhou, X. (2015). A comprehensive selection of reference genes for RT-qPCR analysis in a predatory lady beetle, Hippodamia convergens (Coleoptera: Coccinellidae). PLoS One 10:e125868. doi: 10.1371/journal.pone.0125868

PubMed Abstract | CrossRef Full Text | Google Scholar

Pfaffl, M. W., Tichopad, A., Prgomet, C., and Neuvians, T. P. (2004). Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper-Excel-Based tool using pair-wise correlations. Biotechnol. Lett. 26, 509–515. doi: 10.1023/B:BILE.0000019559.84305.47

CrossRef Full Text | Google Scholar

Popovici, V., Goldstein, D. R., Antonov, J., Jaggi, R., Delorenzi, M., and Wirapati, P. (2009). Selecting control genes for RT-qPCR using public microarray data. BMC Bioinformatics 10:42. doi: 10.1186/1471-2105-10-42

PubMed Abstract | CrossRef Full Text | Google Scholar

Qu, R., Miao, Y., Cui, Y., Cao, Y., Zhou, Y., Tang, X., et al. (2019). Selection of reference genes for the quantitative real-time PCR normalization of gene expression in Isatis indigotica fortune. BMC Mol. Biol. 20:9. doi: 10.1186/s12867-019-0126-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodrigues, T. B., Khajuria, C., Wang, H., Matz, N., Cunha Cardoso, D., Valicente, F. H., et al. (2014). Validation of reference housekeeping genes for gene expression studies in western corn rootworm (Diabrotica virgifera virgifera). PLoS One 9:e109825. doi: 10.1371/journal.pone.0109825

PubMed Abstract | CrossRef Full Text | Google Scholar

Rylett, C. M., Walker, M. J., Howell, G. J., Shirras, A. D., and Isaac, R. E. (2007). Male accessory glands of Drosophila melanogaster make a secreted angiotensin I-converting enzyme (ANCE), suggesting a role for the peptide-processing enzyme in seminal fluid. J. Exp. Biol. 210, 3601–3606. doi: 10.1242/jeb.009035

PubMed Abstract | CrossRef Full Text | Google Scholar

Sang, W., He, L., Wang, X. P., Zhu-Salzman, K., and Lei, C. L. (2015). Evaluation of reference genes for RT-qPCR in Tribolium castaneum (Coleoptera: Tenebrionidae) under UVB stress. Environ. Entomol. 44, 418–425. doi: 10.1093/ee/nvv010

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmittgen, T. D., and Livak, K. J. (2008). Analyzing real-time PCR data by the comparative C(T) method. Nat. Protoc. 3, 1101–1108. doi: 10.1038/nprot.2008.73

PubMed Abstract | CrossRef Full Text | Google Scholar

Shakeel, M., Rodriguez, A., Tahir, U. B., and Jin, F. (2018). Gene expression studies of reference genes for quantitative real-time PCR: an overview in insects. Biotechnol. Lett. 40, 227–236. doi: 10.1007/s10529-017-2465-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, X. Q., Guo, W. C., Wan, P. J., Zhou, L. T., Ren, X. L., Ahmat, T., et al. (2013). Validation of reference genes for expression analysis by quantitative real-time PCR in Leptinotarsa decemlineata (Say). BMC Res. Notes 6:93. doi: 10.1186/1756-0500-6-93

PubMed Abstract | CrossRef Full Text | Google Scholar

Shu, B., Zhang, J., Cui, G., Sun, R., Sethuraman, V., Yi, X., et al. (2018). Evaluation of reference genes for real-time quantitative pcr analysis in larvae of Spodoptera litura exposed to azadirachtin stress conditions. Front. Physiol. 9:372. doi: 10.3389/fphys.2018.00372

CrossRef Full Text | Google Scholar

Silver, N., Best, S., Jiang, J., and Thein, S. L. (2006). Selection of housekeeping genes for gene expression studies in human reticulocytes using real-time PCR. BMC Mol. Biol. 7:33. doi: 10.1186/1471-2199-7-33

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, Q. Q., Zhu, L., Li, Y., Liu, W., Ma, W. H., Lei, C. L., et al. (2015). A de novo transcriptome and valid reference genes for quantitative real-time PCR in Colaphellus bowringi. PLoS One 10:e0118693. doi: 10.1371/journal.pone.0118693

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, Y., Zhou, X.-R., and Pang, B.-P. (2017). Reference gene selection and evaluation for expression analysis using qRT-PCR in Galeruca daurica (Joannis). B Entomol. Res. 107, 359–368. doi: 10.1017/s0007485316000948

PubMed Abstract | CrossRef Full Text | Google Scholar

Thellin, O., Zorzi, W., Lakaye, B., De Borman, B., Coumans, B., Hennen, G., et al. (1999). Housekeeping genes as internal standards: use and limits. J. Biotechnol. 75, 291–295. doi: 10.1016/S0168-1656(99)00163-7

CrossRef Full Text | Google Scholar

Tong, Z., Gao, Z., Wang, F., Zhou, J., and Zhang, Z. (2009). Selection of reliable reference genes for gene expression studies in peach using real-time PCR. BMC Mol. Biol. 10:1–13. doi: 10.1186/1471-2199-10-71

PubMed Abstract | CrossRef Full Text | Google Scholar

Toutges, M. J., Hartzer, K., Lord, J., and Oppert, B. (2010). Evaluation of reference genes for quantitative polymerase chain reaction across life cycle stages and tissue types of Tribolium castaneum. J. Agric. Food Chem. 58, 8948–8951. doi: 10.1021/jf101603j

PubMed Abstract | CrossRef Full Text | Google Scholar

Valasek, M. A., and Repa, J. J. (2005). The power of real-time PCR. Adv. Physiol. Educ. 29, 151–159. doi: 10.1152/advan.00019.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Vandesompele, J., De Preter, K., Pattyn, F., Poppe, B., Van Roy, N., De Paepe, A., et al. (2002). Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 3:7. doi: 10.1186/gb-2002-3-7-research0034

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, F., Xiao, P., Chen, D., Xu, L., and Zhang, B. (2012). miRDeepFinder: a miRNA analysis tool for deep sequencing of plant small RNAs. Plant Mol. Biol. 80, 75–84. doi: 10.1007/s11103-012-9885-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, J., Baulding, J., and Palli, S. R. (2013). Proteomics of Tribolium castaneum seminal fluid proteins: identification of an angiotensin-converting enzyme as a key player in regulation of reproduction. J. Proteomics 78, 83–93. doi: 10.1016/j.jprot.2012.11.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamazaki, K., Imai, C., and Natuhara, Y. (2000). Rapid population growth and food-plant exploitation pattern in an exotic leaf beetle, Ophraella communa LeSage (Coleoptera: Chrysomelidae), in western Japan. Appl. Entomol. Zool. 35, 215–223. doi: 10.1303/aez.2000.215

CrossRef Full Text | Google Scholar

Yan, Z., Gao, J., Lv, X., Yang, W., Wen, S., Tong, H., et al. (2016). Quantitative evaluation and selection of reference genes for quantitative RT-PCR in mouse Acute pancreatitis. Biomed. Res. Int. 2016:8367063. doi: 10.1155/2016/8367063

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, C., Pan, H., Liu, Y., and Zhou, X. (2015a). Stably expressed housekeeping genes across developmental stages in the two-spotted spider mite. Tetranychus Urticae 10:1. doi: 10.1371/journal.pone.0120833

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, C., Pan, H., Liu, Y., and Zhou, X. (2015b). Temperature and development impacts on housekeeping gene expression in Cowpea Aphid. Aphis craccivora (Hemiptera: Aphidiae). PLoS One 10:e0130593. doi: 10.1371/journal.pone.0130593

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, C., Preisser, E. L., Zhang, H., Liu, Y., Dai, L., Pan, H., et al. (2016). Selection of reference genes for RT-qPCR analysis in Coccinella septempunctata to assess un-intended effects of RNAi transgenic plants. Front. Plant Sci. 7:1672. doi: 10.3389/fpls.2016.01835

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, X., Pan, H., Yuan, L., and Zhou, X. (2018). Reference gene selection for RT-qPCR analysis in Harmonia axyridis, a global invasive lady beetle. Sci. Rep. 8:2689. doi: 10.1038/s41598-018-20612-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, C., Ma, F., Chen, H., Wan, F., Guo, J., and Zhou, Z. (2018). Heritability and evolutionary potential drive cold hardiness in the overwintering Ophraella communa beetles. Front. Physiol. 9:666. doi: 10.3389/fphys.2018.00666

CrossRef Full Text | Google Scholar

Zhong, M., Wang, X., Wen, J., Cai, J., Wu, C., and Aly, S. M. (2013). Selection of reference genes for quantitative gene expression studies in the house fly (Musca domestica L.) using reverse transcription quantitative real-time PCR. Acta Biochem. Biophys. Sin. 45, 1069–1073. doi: 10.1093/abbs/gmt111

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, L., Chen, F., Ye, J., and Pan, H. (2018). Selection of reliable reference genes for RT-qPCR analysis of Bursaphelenchus mucronatus gene expression from different habitats and developmental stages. Front. Genet. 9:269. doi: 10.3389/fgene.2018.00269

CrossRef Full Text | Google Scholar

Zhou, Z. S., Chen, H. S., Zheng, X. W., Guo, J. Y., Wei, G., Min, L., et al. (2014). Control of the invasive weed Ambrosia artemisiifolia with Ophraella communa and Epiblema strenuana. Biocontrol. Sci. Techn. 24, 950–964. doi: 10.1080/09583157.2014.897305

CrossRef Full Text | Google Scholar

Zhou, Z. S., Guo, J. Y., Chen, H. S., and Wan, F. H. (2010a). Effects of temperature on survival, development, longevity, and fecundity of Ophraella communa (Coleoptera: Chrysomelidae), a potential biological control agent against Ambrosia artemisiifolia (Asterales: Asteraceae). Environ. Entomol. 39, 1021–1027. doi: 10.1603/EN09176

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Z. S., Guo, J. Y., Michaud, J. P., Li, M., and Wan, F.-H. (2010b). Variation in cold hardiness among geographic populations of the ragweed beetle, Ophraella communa LeSage (Coleoptera: Chrysomelidae), a biological control agent of Ambrosia artemisiifolia L. (Asterales: Asteraceae), in China. Biol. Invas. 13, 659–667. doi: 10.1007/s10530-010-9857-x

CrossRef Full Text | Google Scholar

Keywords: Ophraella communa, RT-qPCR, reference gene, normalization, gene expression

Citation: Zhang Y, Chen J, Chen G, Ma C, Chen H, Gao X, Tian Z, Cui S, Tian Z, Guo J, Wan F and Zhou Z (2020) Identification and Validation of Reference Genes for Quantitative Gene Expression Analysis in Ophraella communa. Front. Physiol. 11:355. doi: 10.3389/fphys.2020.00355

Received: 22 January 2020; Accepted: 26 March 2020;
Published: 07 May 2020.

Edited by:

Xanthe Vafopoulou, York University, Canada

Reviewed by:

Yong Liu, Hunan Academy of Agricultural Sciences (CAAS), China
Shaoli Wang, Institute of Vegetables and Flowers (CAAS), China

Copyright © 2020 Zhang, Chen, Chen, Ma, Chen, Gao, Tian, Cui, Tian, Guo, Wan and Zhou. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Zhongshi Zhou, emhvdXpob25nc2hpQGNhYXMuY24=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.