Abstract
G protein-coupled receptors (GPCRs) are the largest family of membrane receptors in animals and humans, which transmit various signals from the extracellular environment into cells. Studies have reported that several GPCRs transmit the same signal; however, the mechanism is unclear. In the present study, we identified all 122 classical GPCRs from the genome of Helicoverpa armigera, a lepidopteran pest species. Twenty-four GPCRs were identified as upregulated at the metamorphic stage by comparing the transcriptomes of the midgut at the metamorphic and feeding stages. Nine of them were confirmed to be upregulated at the metamorphic stage. RNA interference in larvae revealed the prolactin-releasing peptide receptor (PRRPR), smoothened (SMO), adipokinetic hormone receptor (AKHR), and 5-hydroxytryptamine receptor (HTR) are involved in steroid hormone 20-hydroxyecdysone (20E)-promoted pupation. Frizzled 7 (FZD7) is involved in growth, while tachykinin-like peptides receptor 86C (TKR86C) had no effect on growth and pupation. Via these GPCRs, 20E regulated the expression of different genes, respectively, including Pten (encoding phosphatidylinositol-3,4,5-trisphosphate 3-phosphatase), FoxO (encoding forkhead box O), BrZ7 (encoding broad isoform Z7), Kr-h1 (encoding Krüppel homolog 1), Wnt (encoding Wingless/Integrated) and cMyc, with hormone receptor 3 (HHR3) as their common regulating target. PRRPR was identified as a new 20E cell membrane receptor using a binding assay. These data suggested that 20E, via different GPCRs, regulates different gene expression to integrate growth and development.
Introduction
G protein-coupled receptors (GPCRs) are present widely in animals and humans (Hanlon and Andrew, 2015). GPCRs sense and transmit external stimuli into cells to regulate a variety of physiological processes, including cognition, metabolism, inflammation, immunity, and cell proliferation (Rasmussen et al., 2011). There are more than 800 GPCRs encoded in the human genome (Fredriksson and Schioth, 2005), over 1,300 GPCRs in mice, 116 classical GPCRs, which can act as guanine nucleotide exchange factors (GEFs), in Drosophila (Hanlon and Andrew, 2015), and 276 in Anopheles gambiae (Hill et al., 2002). The importance of GPCRs in cellular signaling has resulted in ∼34% of human drugs acting at GPCRs (Hauser et al., 2017). GPCRs are also suggested as targets for next generation pesticides (Audsley and Down, 2015). An interesting phenomenon in GPCRs-mediated signaling is that several GPCRs transmit the same signal of a ligand. For example, nine GPCRs function for adrenaline and five for dopamine (Hauser et al., 2017); however, the mechanism is unclear.
G protein-coupled receptors also transmit animal steroid hormone signals in the cell membrane. For example, GPCR 30 (GPR30/GPER) is an estrogen cell membrane receptor and transmits estrogen signals in mammals (Maggiolini and Picard, 2010). The dopamine/ecdysteroid receptor (DopEcR) transmits the non-genomic signal of insect molting hormone 20-hydroxyecdysone (20E) in Drosophila (Srivastava et al., 2005) and in Helicoverpa armigera (Kang et al., 2019). To date, several GPCRs have been proven to transmit 20E signals in H. armigera (Zhao, 2020), including ecdysone-responsible GPCR 1 (ErGPCR-1), ecdysone-responsible GPCR 2 (ErGPCR-2), and ecdysone-responsible GPCR 3 (ErGPCR-3) (Cai et al., 2014a; Wang et al., 2015; Kang et al., 2021). These data suggest that several GPCRs function as steroid hormone receptors; however, whether any other GPCRs transmit 20E signals, and the mechanism by which several GPCRs function in 20E signaling, are unclear.
Helicoverpa armigera is a widespread lepidopteran agricultural pest (Wu et al., 2008). We used H. armigera as the research model to identify the new GPCRs involved in 20E signaling and addressed the mechanism of their function in the 20E pathway. In all, 122 GPCRs were identified from the genome of H. armigera. Six GPCRs transmit 20E signal for hormone receptor 3 (HHR3) expression, a 20E-induced delayed early gene (Palli et al., 1997). 20E, via different GPCRs, regulates the expression of various genes, including Pten (encoding phosphatidylinositol-3,4,5-trisphosphate 3-phosphatase), FoxO (encoding transcription factor forkhead box O), which are known playing roles in 20E signaling (Cai et al., 2016), BrZ7 (encoding broad isoform Z7), a transcription factor that promotes metamorphosis (Cai et al., 2014b), Kr-h1 (encoding Krüppel homolog 1), the antimetamorphic effector induced by juvenile hormone (JH) (Belles, 2020), and Wnt (encoding Wingless/Integrated) and cMyc, which play significant roles in insect growth and development (Clevers, 2006; Kayukawa et al., 2012; Cai et al., 2014b), to integrate growth and metamorphosis. One GPCR, prolactin-releasing peptide receptor (PRRPR), was determined to bind 20E. Our study presents an example to explain the mechanism by which several GPCRs transmit the same signal.
Results
Identification of Helicoverpa armigera G Protein-Coupled Receptors From the Genome
We searched for all GPCRs from the genome of H. armigera to identify classification of those GPCRs that are involved in 20E signaling. We found 122 genes encoding classical GPCRs in the H. armigera genome1 using BLAST (Basic Local Alignment Search Tool) with Drosophila and Bombyx mori GPCRs, respectively. Having removed four GPCRs with large sequence differences, the sequences of 118 presumed GPCRs, named as they are in the genome, were used to create a phylogenetic tree. These GPCRs could be divided into three clades according to four major categories of GPCRs (Sadowski and Parish, 2003): Class A (89 sequences), class B (15 sequences) and class C or F (14 sequences) (Figure 1 and Supplementary Table 1). Some GPCRs from D. melanogaster and Homo sapiens were used as landmarks of the classes, respectively. Fifteen GPCRs annotated to class A, B, C, or F in the genome were reclassified in different classes according to the sequences, which are marked with the related colors in each class in Figure 1. Twenty-five GPCRs that were not classified in the genome were gathered to different classes according to our phylogenetic analysis, which are marked in black in Figure 1. Four GPCRs known to transmit 20E signals in H. armigera, were classified as Class A (DopEcR) and class B (ErGPCR-1, ErGPCR-2, and ErGPCR-3).
FIGURE 1
Screening of the G Protein-Coupled Receptors in 20E Signaling
To screen for GPCRs involved in 20E signaling pathway, we compared the transcriptomes of the midgut at the feeding stage (6th–24 h) and the metamorphic molting stage (6th–72 h). Twenty-four GPCRs were found to be upregulated and seven were downregulated in the metamorphic stage (Supplementary Figure 1), suggesting that these twenty-four GPCRs are involved in metamorphosis. To examine the transcriptome analysis, 13 of the GPCRs (11 upregulated and 2 downregulated) from different classes were selected and examined for their developmental expression profiles in tissues using quantitative real-time reverse transcription PCR (qRT-PCR) to validate the result of the transcriptome analysis. Three GPCR genes, PrRPR, Akhr (encoding adipokinetic hormone receptor), and Fzd7 (encoding frizzled 7), showed increased expression during metamorphosis (MM to P) in four detected tissues (Figure 2). Six GPCR genes, Smo (encoding smoothened), Htr (encoding 5-hydroxytryptamine receptor), TkR86C (encoding tachykinin-like peptides receptor 86C), Fshr (encoding follicle-stimulating hormone receptor), Rya-R (encoding the RYamide receptor), and Npfr (encoding neuropeptide F receptor) showed increased expression during metamorphosis in some tissues (Figure 3). Four GPCR genes, Galr3 (encoding galanin receptor type 3), GPCR4 (encoding the uncharacterized protein LOC110374861), Opsin (encoding red-sensitive opsin), and Fzd4 (encoding frizzled 4) did not show increased expression during metamorphosis (Supplementary Figure 2). These results confirmed that the expression levels of nine GPCRs genes increased during metamorphosis, with or without tissue specificity, and might play roles in 20E-promoted metamorphosis.
FIGURE 2
FIGURE 3
G Protein-Coupled Receptors Have Different Functions in 20E Pathway
Among the nine GPCRs that showed increased expression during metamorphosis, three GPCRs-Fshr, Rya-R, and Npfr were not successfully knocked down in larvae using RNA interference (RNAi). Other six GPCRs, PrRPR, Smo, Akhr, Htr, Fzd7, and TkR86C, which showed increased expression during metamorphosis in all tissues or in some tissues, were knocked down in larvae using RNAi to examine their roles in 20E-promoted earlier pupation. In the dsGFP plus 20E treatment group, 91.1% of the larvae pupated in 111 h (timed from the 6th instar 6 h to pupae). However, knockdown of PrRPR caused 63.3% of the larvae delayed pupation for 36 h, and increased death, compared with dsGFP plus 20E (Figures 4A–C). In addition, the midgut did not show a red color, a sign of programmed cell death (Wang et al., 2007; Hakim et al., 2010), in the dsPrRPR plus 20E treatment group compared with that in the dsGFP plus 20E control (Figure 4D). Hematoxylin and eosin (HE) staining showed that the imaginal midgut formed after dsGFP control injection, indicating the occurrence of midgut remodeling. In contrast, the imaginal midgut did not form after dsPrRPR injection for 60 h (Figure 4E). Similar results were obtained after knockdown of Smo, Akhr, and Htr. Compared with the dsGFP + 20E group, 53–65% of larvae delayed pupation for 24–43 h, and the midgut did not change to red or remodel on time (Supplementary Figures 3–5). These results suggested that these four GPCRs play roles in 20E-promoted pupation.
FIGURE 4
However, knockdown of Fzd7 caused 57.8% of the larvae to form small pupae (Figures 5A,B). The pupal weight decreased to an average of 0.29 g compared with 0.43 g of the dsGFP injection control, with no significant difference in pupation time compared with the control group (Figures 5C,D). These results suggested that FZD7 is involved in larval growth. However, knockdown of TkR86C resulted in no abnormal phenotype (Supplementary Figure 6). The efficacy of RNAi was confirmed for these GPCRs, and except for Smo and Fzd7, which decreased after knockdown of Htr, no off target effects were detected for the other GPCRs (Supplementary Figure 7). These results showed that different GPCRs in the 20E signaling pathway play different roles in regulating growth and metamorphosis.
FIGURE 5
20E, via Different G Protein-Coupled Receptors, Regulates Gene Expression
The mechanism by which knockdown of the six GPCRs caused different outcomes was addressed by examining gene expression, including HHR3, Pten, FoxO, and BrZ7, which play roles in 20E-induced metamorphosis, Kr-h1, which plays role in JH signaling, Wnt and cMyc, which play roles in growth. qRT-PCR analysis showed that the expression levels of the six GPCR genes were upregulated by 20E in the midgut, confirming their responses to 20E induction. After knockdown of the six GPCR genes by RNAi in larvae, the mRNA levels of HHR3 decreased, suggesting that these six GPCRs play roles in 20E signaling. However, the expression of Pten and FoxO decreased only after PrRPR and Smo knockdown (Figures 6A,B), but not after Akhr, Htr, Fzd7, and TkR86C knockdown (Figures 6C–F). Furthermore, BrZ7 expression decreased after knockdown of PrRPR and Smo, Kr-h1 expression increased after Akhr and Htr knockdown, and Wnt and cMyc decreased after knockdown of Fzd7 (Figure 7). These results revealed that 20E, via different GPCRs, integrates insect pupation and growth by regulating the expression of various genes.
FIGURE 6
FIGURE 7
To support this conclusion, the previous reported GPCRs that transmit the 20E signal, ErGPCR-1, ErGPCR-2, ErGPCR-3, and DopEcR, were examined for their regulation of gene expression. The results showed that the expression of HHR3 decreased after knockdown of ErGPCR-1, ErGPCR-2, ErGPCR-3, and DopEcR; however, Pten and FoxO expression decreased after ErGPCR-1 knockdown, but not after ErGPCR-2, ErGPCR-3, or DopEcR knockdown (Supplementary Figure 8A), which confirmed that different GPCRs regulate the expression of different genes in the 20E pathway. The RNA interference efficiency of these four GPCRs was confirmed (Supplementary Figures 8B–E).
Prolactin-Releasing Peptide Receptor Binds 20E
To identify new GPCR functioning as cell membrane receptor of 20E, PRRPR and SMO were further analyzed for their binding 20E to determine their receptor roles in 20E signaling. Surflex-Dock (SFXC) in SYBYL X2.0 software (Certara, Princeton, NJ, United States) was used to dock 20E to PRRPR and SMO to predict the possibility of PRRPR and SMO binding 20E (Figures 8A,B). 20E forms hydrogen bonds with Ala-61, Gly-64, and Pro-316, of PRRPR, and Gln-314 and Glu-219 of SMO (Figures 8C,D). The scores for PRRPR and SMO binding to 20E were 2.96, and –0.78, respectively. These data predicted that PRRPR has a higher binding ability to 20E than SMO.
FIGURE 8
PRRPR-GFP and SMO-GFP were overexpressed in an H. armigera epidermal cell line (HaEpi) to address their binding to 20E. Green fluorescent protein (GFP) was overexpressed as a tag control. The overexpressed PRRPR-GFP and SMO-GFP were confirmed to be located in the cell membrane (Figure 9A). A binding assay using a 20-hydroxyecdysone enzyme immunoassay (20E-EIA) showed that the amount of 20E bound by the cell membrane from the PRRPR-GFP-overexpressing cells increased significantly compared with that bound by the GFP-overexpressing cells. However, the amount of 20E bound by cell membranes from SMO-GFP overexpressing cells did not increase compared with that of the GFP-overexpressing cells (Figure 9B). These results suggested that PRRPR could bind 20E in the cell membrane.
FIGURE 9
A saturation-binding curve was constructed using 20E-EIA to further examine the affinity of GPCRs to 20E by calculating their dissociation constants (Kds). The saturable specific binding of cell membranes from cells overexpressing PRRPR-GFP to 20E had a Bmax of 2.096 ± 0.1037 nmol/mg protein and a Kd of 12.76 ± 2.192 nM. In comparison, the saturation binding of cell membranes from cells overexpressing GFP to 20E had a Bmax of 1.195 ± 0.1007 nmol/mg protein and a Kd of 30.2 ± 6.452 nM (cells overexpressing GFP still have other GPCRs on their cell membranes) (Figure 9C). The 20E-EIA assay is based on competition between the unlabeled 20E (20E bound to GPCR) and acetyl choline esterase (AChE)-labeled 20E (Tracer) for the limited-specific rabbit anti-20E antiserum; therefore, an inhibition or competitive curve was not detected. These data confirmed that the PRRPR-GFP could bind 20E.
Discussion
This research identified and classified all classical GPCRs in the H. armigera genome. The GPCRs that function as 20E receptors were classified in classes A, B and class C or F. Further study revealed that different GPCRs showed different expression profiles and mediated the expression of different genes in 20E signaling. PRRPR was determined as a new GPCR cell membrane receptor. These data explained the mechanism by which several GPCRs are involved in the signaling of the same ligand.
Identification and Classification of G Protein-Coupled Receptors in Helicoverpa armigera Genome
We identified 122 genes encoding classical GPCRs in the H. armigera genome. The GPCRs were classified into categories A, B, C, or F, which was relatively consistent with the classification of GPCRs in Drosophila (Hanlon and Andrew, 2015). Most of the GPCRs were classified consistently with their classification in the genome; however, some GPCRs were mixed in different classes in our study when using the full open reading frames (ORFs). We found 19 Mth-like GPCRs in the H. armigera genome, which is close to the 16 Mth-like GPCRs in D. melanogaster (Patel et al., 2012), but more than the 7 found in Anopheles, and the 2 found in B. mori (Fan et al., 2010). The Mth-like GPCRs play various roles in regulating the metabolism, aging, and self-balance to high temperature, hunger, dryness, and oxidative damage (Friedrich and Jones, 2016). In insects, Mth-like GPCRs are known to play roles in the setting of the endogenous circadian clocks (Mertens et al., 2007), regulation of fluid and ion secretion (Reagan, 1994), as well as the stress response and longevity (Lin et al., 1998).
20E, via Different G Protein-Coupled Receptors, Regulates Gene Expression
The involvement of several GPCRs in a same signal, such as 20E signaling, is an intriguing phenomenon. The differences among the GPCRs in 20E signaling have been explained by their induced downstream effects, including ErGPCR-1 inducing the Ca2+-PKC signaling, while ErGPCR-2 inducing the GPCR-cAMP-PKA and GPCR-Ca2+-PKC signaling, increasing 20E entry, and being internalized by 20E induction. DopEcR of H. armigera directly interacts with Gαs and Gαq under the induction of 20E to increase the levels of cAMP and Ca2+ (Zhao, 2020). ErGPCR-3 has very similar characteristics to ErGPCR-2 (Kang et al., 2021). Different GPCRs can cross react with different G proteins (Flock et al., 2017). Here, we further revealed that GPCRs have quite different expression profiles in tissues and at different developmental stages. Importantly, 20E, via different GPCRs, regulates the expression of various genes, including via PRRPR and SMO, which upregulate the expression of Pten, FoxO, and BrZ7 to promote pupation. 20E, via AKHR and HTR, represses the expression of Kr-h1 to promote pupation. 20E via FZD7 upregulates the expression of Wnt and cMyc to promote growth. 20E signaling also promotes wing disk development (Mirth et al., 2009). By the integration and competition of different signals induced by different ligands in vivo, 20E regulates pupation. These findings in 20E signaling first revealed the mechanism by which several GPCRs transmit the same signal to regulate the expression of different genes in the network of the cells. In this work, we performed the screen based on the expression levels of GPCRs. There might be GPCRs transmit external signals in an expression-independent manner, which needs further study.
Our results suggested PRRPR, SMO, AKHR, HTR, FZD7, and TKR86C are involved in 20E-inducing HHR3 expression, suggesting that they transmit 20E signals. In humans, PRRPR is a neuropeptide prolactin receptor (Dodd and Luckman, 2013). Human SMO participates in hedgehog signaling to guide cell differentiation, proliferation, and survival (Wu et al., 2017). FZD7 is the most important WNT receptor involved in cancer development and progression in mammals (King et al., 2012). In insects, AKHR binds AKH (adipokinetic hormone) to increase lipolysis, glycogenolysis, and trehalose production (Van der Horst et al., 2001; Bednarova et al., 2013; Baumbach et al., 2014). HTR (5-HT receptor) plays a key role in morphogenesis in the insect nervous system (Blenau and Thamm, 2011). TkR86C is the neurokinin K receptor in D. melanogaster that plays a role in neuromodulation in the central nervous system, participating in the processing of sensory information and the control of motor activities (Vanden Broeck et al., 1999). The role of TkR86C in insect needs further study. Here, we revealed a new function of these GPCRs in 20E signaling.
G Protein-Coupled Receptors Can Transmit 20E Signals Whether They Bind 20E or Not
It has been suggested that cells or cell membranes that overexpress GPCRs can bind steroid hormones in Drosophila (Srivastava et al., 2005) and mammals (Maggiolini and Picard, 2010). We found that PRRPR could bind 20E with the saturable specific binding Kd of 12.76 ± 2.192 nM. However, SMO could not bind 20E, although SMO transmits the 20E signal and is involved in 20E-induced pupation. In our previous study, we found that ErGPCR-1, ErGPCR-2, ErGPCR-3, and DopEcR can transmit 20E signals in H. armigera. ErGPCR-2, ErGPCR-3, and DopEcR can bind 20E, but ErGPCR-1 cannot (Kang et al., 2019, 2021). These data suggested that GPCRs transmit 20E signal with or without binding 20E. This might be because GPCRs loosely or dynamically bind their ligands (Nygaard et al., 2013; Strasser et al., 2017). Another possibility is that 20E competes with another ligand, such as dopamine, in H. armigera (Kang et al., 2019). GPCRs might also play roles in other pathways after upregulation by 20E, which requires further study.
G protein-coupled receptors share a seven transmembrane domain structural architecture (Latorraca et al., 2017). Except ErGPCR-1, ErGPCR-2, ErGPCR-3, and DopEcR have been reported to bind 20E. In addition to DopEcR in class A, ErGPCR-1, ErGPCR-2, and ErGPCR-3 belong to the Mth-like GPCRs in class B. An important feature of Mth-like GPCRs is the presence of 10 cysteine residues that form five disulfide bonds (West et al., 2001). The long N-terminal domains tend to recognize peptide ligands, such as hormones and neuropeptides (Cardoso et al., 2005, 2010), such as secretin; therefore, these GPCRs are also called secretin receptors (Cvejic et al., 2004; Ja et al., 2009). The relationship between the structure of GPCR and its binding 20E is unclear now. Up to date, several GPCRs can bind 20E, therefore, the upregulation of GPCR expression by 20E is likely a key factor to perform their functions. The mechanism that 20E upregulates GPCR expression differentially needs further studied. Mth-like GPCRs are not present in vertebrates, but are more abundant in arthropods (Patel et al., 2012; Araujo et al., 2013; de Mendoza et al., 2016), and thus represent targets for insecticides.
Conclusion
There are 122 classical GPCRs in the genome of H. armigera. The GPCRs that transmit 20E signal were classified in classes A, B and class C or F. Various GPCRs transmit 20E signals according to their different expression patterns in tissues and their increased expression during metamorphosis. 20E, via different GPCRs, regulates the expression of various genes, thus promoting pupation by integrating different signals in vivo. PRRPR binds 20E and is a newly identified 20E cell membrane receptor (Figure 10).
FIGURE 10
Materials and Methods
Identification of G Protein-Coupled Receptors
Putative H. armigera GPCRs were identified in four steps: First, we downloaded all the GPCR protein sequences of Drosophila (Hanlon and Andrew, 2015) and B. mori (Fan et al., 2010). B. mori protein sequences were retrieved from the NCBI sequence repository2. We obtained 90 classical GPCRs from B. mori. Drosophila protein sequences were retrieved from FlyBase3 and the NCBI database. We obtained 116 classical Drosophila GPCRs. Second, preliminary screening to obtain putative GPCRs of H. armigera was performed using NCBI BLAST based on downloaded Drosophila and B. mori GPCRs. We queried the H. armigera proteome4 using each GPCR sequence from Drosophila and B. mori separately and selected the protein sequences with the highest scores. Third, the protein sequences with highest score were then used as query sequences one by one in a BLAST search against the proteome from H. armigera to obtain other GPCRs sequences that were not found in the previous step. Finally, we removed the repetitive sequences in the protein sequence obtained in the above steps. Then, NCBI conserved domain search service (CD search)5 and SMART online software6 were used to predict the structure of these protein sequences, and the GPCRs were seven-transmembrane domain proteins (7TMPs) were obtained.
GPCRs were identified from the transcriptomes of 6th–24 h larvae and 6th–72 h larvae. The transcriptomes were analyzed once without technique replicates. However, the samples were collected from several larvae to normalize the individual differences. The mRNA levels of GPCRs were examined after injection of 20E or Double-Stranded RNA (dsRNA) using qRT-PCR, with an equal amount of diluted DMSO as a solvent control.
Phylogenetics Analyses
The classification of Drosophila and Homo sapiens proteins in each GPCR family is clear and detailed. We classified the potential GPCRs in the H. armigera genome into various categories based on sequence homology. Drosophila and Homo sapiens GPCRs were used as guides, and the MEGA 6.0 software was used to construct phylogenetic trees using the Neighbor-Joining method with 1000 bootstraps (Tamura et al., 2013).
Insects
Helicoverpa armigera larvae were raised in the insect culture room at 25–27°C under a photoperiod of 14 h light/10 h dark. The larvae are reared on a previously described artificial diet (Zhao et al., 1998).
Quantitative Real-Time Reverse Transcription PCR
The total RNA was extracted using the Trizol reagent (TransGen Biotech, Beijing, China). cDNA was synthesized from the total RNA using a FastQuant RT Kit (Tiangen Biotech, Beijing, China). qRT-PCR was then carried out using the cDNA as the template in a CFX96 real-time system (Bio-Rad, Hercules, CA, United States) with 2 × SYBR qRT-PCR pre-mixture (TransGen Biotech). All the primers used are listed in Supplementary Table 2. Actb encodes a type of actin, which is a structural component of the cytoskeleton microfilaments. The Actb gene is highly conserved and highly expressed at the mRNA level (Butet et al., 2014). In many studies, including studies on different developmental stages and different tissues in H. armigera (Di et al., 2020), Actb is considered a suitable internal reference gene (Lu et al., 2013; Gao et al., 2017). H. armigera Actb (encoding beta actin; GenBank accession no. EU52707) was used as the internal standard. All data were from at least three biological replicates and were analyzed using the 2−ΔΔCT method (ΔΔCT = ΔCTtreated sample-ΔCTcontrol, ΔCT = CTgene-CTActb) (Livak and Schmittgen, 2001).
20E Induction in Larvae
The 20E powder (10 mg) was dissolved in 1 mL DMSO as a storage solution and diluted with phosphate-buffered saline (PBS; 140 mM NaCl, 10 mM sodium phosphate, pH 7.4) for experiment. 20E was injected into the hemocoel from the side of the larval abdomen. The control groups were treated with the equal amount of diluted DMSO.
Double-Stranded RNA Synthesis
RNA interference (RNAi) has been used for gene knockdown in many moths (Xu et al., 2016). The long dsRNA is broken down into smaller fragments in vivo (Zamore et al., 2000) and specifically and successfully inhibits the expression of target genes in worms (Fire et al., 1998). DNA fragments– 5′-583 bp-1037 bp-3′ of Fzd7, 5′-586 bp-1109 bp-3′ of Htr, 5′-43 bp-586 bp-3′ of PrRPR, 5′-1195 bp-1858 bp-3′ of Smo, 5′-641 bp-1045 bp-3′ of Akhr, 5′-65 bp-740 bp-3′ of TkR86C– were amplified as the template for dsRNA synthesis using the primers RNAiF and RNAiR. A T7 promoter sequence was added to the RNAi primers (Supplementary Table 2). The cDNA of the target gene was amplified using a single PCR reaction and was used as the template to synthesize dsRNA. The dsRNA was synthesized using MEGAscript RNAi kit (Ambion, Austin, TX, United States) according to the instruction manual. Next, the product was purified using the phenol-chloroform method. The quality of the synthesized dsRNA was quantified using a micro-spectrophotometer and detected using 1% agarose gel electrophoresis.
RNA Interference of Genes in Larvae
The dsRNA was diluted with PBS. The sixth instar 6 h larvae were placed on ice for 15 min until they did not move. A sterile micro syringe was used to inject 1 μg of dsRNA into the hemocoel from the side of the larval abdomen (taking care not to touch the midgut). dsRNAs were injected three times at 24 h intervals. The control groups were treated with the same amount of dsGFP. Each experimental group and control group contained 30 larvae and three independent biological replicates were performed. Total RNA was extracted using Trizol reagent (TransGen Biotech) and qRT-PCR was performed to detect the effects of RNAi at 24 h after the last injection.
Hematoxylin-Eosin Staining
The midgut dissected from the larva was washed with PBS, and then fixed in 4% paraformaldehyde at 4°C overnight. The fixed tissue was submitted to a professional company (Servicebio, Wuhan, China) for processing into glass slides and for Hematoxylin-Eosin (HE) staining.
Overexpression of Prolactin-Releasing Peptide Receptor and Smoothened in HaEpi Cells
The pIEx-4-GFP-His vector that was fused with a sequence encoding the green fluorescent protein (GFP) and used for experiments in the insect cell line. The open reading frames (ORFs) of PrRPR (GenBank accession no. XP_021184170.1) and Smo (GenBank accession no. XP_021189185.1) were amplified using primers (Supplementary Table 1) and inserted into the vector. Then, 5 μg of the recombinant plasmids were transfected into HaEpi cells using the QuickShuttle-enhanced transfection reagent (Biodragon Immunotech, Beijing, China). After 48 h of transfection, further experiments were conducted.
Immunocytochemistry
After PRRPR-GFP and SMO-GFP were overexpressed for 48 h, HaEpi cells were washed three times with 500 μL of Dulbecco’s phosphate-buffered saline (DPBS; 137 mM NaCl, 2.7 mM KCl, 1.5 mM KH2PO4 and 8 mM Na2HPO4, pH 7.4), and fixed with PBS containing 4% paraformaldehyde for 10 min in the dark at room temperature. The fixed cells were washed three times for 3 min each. The plasma membrane was stained using wheat germ agglutinin (WGA; Sigma-Aldrich, St. Louis, MO, United States; 1 μg/mL in PBS) in the dark for 4 min and then washed with PBS six times. Nuclei were stained with 4′-6-diamidino-2-phenylindole dihydrochloride (DAPI; Sangon Biotech, Shanghai, China; 1 μg/mL in PBS) in the dark at room temperature for 10 min and then washed with PBS six times. Fluorescence was detected using an Olympus BX51 fluorescence microscope (Olympus, Tokyo, Japan). The negative control (GFP expression) was treated following the same method.
20-Hydroxyecdysone Enzyme Immunoassay
The 20-hydroxyecdysone enzyme immunoassay (20E-EIA) is based on the competition between unlabeled 20E (free 20E) and acetylcholinesterase (AChE)-labeled 20E (Tracer) for limited specific rabbit anti-20E antiserum. The rabbit anti-20E antiserum was combined with the mouse anti-rabbit monoclonal antibody coated-plate. Then, the plate was washed using the wash buffer included with the 20-Hydroxyecdysone Enzyme Immunoassay kit (20E-EIA kit) (Bertin Pharma, Paris, France) (2 mL of concentrated Wash Buffer #A17000 was diluted by 800 mL of UltraPure water then added 400 μL of Tween20 #A12000) to remove all unbound reagents. Then, tracer and free 20E in samples were added into the wells and the plates were incubated at 4°C overnight. After washing the plate five times with wash buffer, 200 μL Ellman’s reagent (an enzymatic substrate for AChE and a chromogen) was added to the wells, and the plate was then incubated with an orbital shaker at 400 rpm in the dark at room temperature. AChE-labeled 20E acts on the substrate in Ellman’s Reagent to form a yellow compound, which can strongly absorb light at 414 nm. The intensity of the color was detected using a spectrophotometer (Infinite M200PRO NanoQuant, Tecan, Grödig, Austria) at 414 nm. The optical density was proportional to the amount of tracer bound to the well and inversely proportional to the amount of 20E in the sample. The 20E standard curve generated by this method was used to determine the quantity of 20E bound to cell membrane proteins.
Detection of the 20E Quantity Bound by the Cell Membranes of HaEpi Cells
PRRPR-GFP and SMO-GFP were overexpressed in HaEpi cells in a 25 cm2 cell culture bottle, respectively. After washing with DPBS twice, the cells were incubated in Grace’s medium containing 1 μM 20E for 5 min at 27°C to allow 20E to bind to the cell membrane. The cells were then collected by centrifugation at 1,700 × g at 4°C for 5 min and the pellet was resuspended in 500 μL enzyme immunoassay (EIA) buffer (Bertin Pharma, Paris, France) and sonicated for 5 min. The pelleted cell membrane debris was resuspended in 100 μL EIA buffer after centrifugation at 4°C at 48,000 × g for 1 h. Then, 50 μg of cell membrane proteins with fixed 20E in 50 μL EIA buffer was added with 450 μL EIA buffer and used to quantify 20E. The 20E-EIA kit was used to detect cell membrane bound-20E according to the manufacturer’s instructions.
Statistical Analysis
Two-group datasets were analyzed using Student’s t-test and in the figures, an asterisk represents a significant difference (p < 0.05) and two asterisks represent an extremely significant difference (p < 0.01). Analysis of variance (ANOVA) was used for multiple comparisons and in the figures, different lowercase letters indicate significant differences (p < 0.05), and the bars indicate the mean ± standard deviation (SD) of three biological replicates. The details are provided in the figure legends.
Publisher’s Note
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.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
Y-LL and X-FZ designed the experiments and wrote the manuscript. D-JD, J-XW, and X-FZ conceived the idea. Y-LL, Y-XL, X-PW, X-LK, and K-QG performed the experiments. Y-LL and Y-XL performed the data analyses. All the authors contributed to the article and approved the submitted version.
Funding
This project was supported by funding from the National Natural Science Foundation of China (http://www.nsfc.gov.cn/publish/portal0/default.htm) (Grant No. 31730083) to X-FZ. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Acknowledgments
We thank Caiyun Sun from State Key Laboratory of Microbial Technology, Shandong University for the assistance of using Infinite M200PRO NanoQuant.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.753787/full#supplementary-material
Footnotes
1.^https://www.ncbi.nlm.nih.gov/Taxonomy/Browser/wwwtax.cgi?mode=Info&id=29058
2.^http://www.ncbi.nlm.nih.gov/
4.^https://www.ncbi.nlm.nih.gov/genome/?term=Helicoverpa+armigera
5.^https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi
6.^http://smart.embl-heidelberg.de/smart/set_mode.cgi?NORMAL=1
References
1
AraujoA. R.ReisM.RochaH.AguiarB.Morales-HojasR.Macedo-RibeiroS.et al (2013). The Drosophila melanogaster methuselah gene: a novel gene with ancient functions.PLoS One8:e63747. 10.1371/journal.pone.0063747
2
AudsleyN.DownR. E. (2015). G protein-coupled receptors as targets for next generation pesticides.Insect Biochem. Mol. Biol.6727–37. 10.1016/j.ibmb.2015.07.014
3
BaumbachJ.XuY.HehlertP.KuhnleinR. P. (2014). Gαq, Gγ1 and Plc21C control Drosophila body fat storage.J. Genet. Genomics41283–292. 10.1016/j.jgg.2014.03.005
4
BednarovaA.KodrikD.KrishnanN. (2013). Unique roles of glucagon and glucagon-like peptides: parallels in understanding the functions of adipokinetic hormones in stress responses in insects.Comp. Biochem. Physiol. A Mol. Integr. Physiol.16491–100. 10.1016/j.cbpa.2012.10.012
5
BellesX. (2020). Kruppel homolog 1 and E93: the doorkeeper and the key to insect metamorphosis.Arch. Insect Biochem. Physiol.103:e21609. 10.1002/arch.21609
6
BlenauW.ThammM. (2011). Distribution of serotonin (5-HT) and its receptors in the insect brain with focus on the mushroom bodies: lessons from Drosophila melanogaster and Apis mellifera.Arthropod. Struct. Dev.40381–394. 10.1016/j.asd.2011.01.004
7
ButetN.SolihinD.SoewardiK.SaefuddinA. (2014). Actin gene from blood cockle Anadara granosa as a potential housekeeping gene for gene expression analysis.Emir. J. Food Agr.26730–736. 10.9755/ejfa.v26i8.15765
8
CaiM. J.DongD. J.WangY.LiuP. C.LiuW.WangJ. X.et al (2014a). G protein-coupled receptor participates in 20-hydroxyecdysone signaling on the plasma membrane.Cell Commun. Signal.12:9. 10.1186/1478-811X-12-9
9
CaiM. J.LiuW.PeiX. Y.LiX. R.HeH. J.WangJ. X.et al (2014b). Juvenile hormone prevents 20-hydroxyecdysone-induced metamorphosis by regulating the phosphorylation of a newly identified broad protein.J. Biol. Chem.28926630–26641. 10.1074/jbc.M114.581876
10
CaiM. J.ZhaoW. L.JingY. P.SongQ.ZhangX. Q.WangJ. X.et al (2016). 20-Hydroxyecdysone activates Forkhead box O to promote proteolysis during Helicoverpa armigera molting.Development1431005–1015. 10.1242/dev.128694
11
CardosoJ. C. R.ClarkM. S.VieraF. A.BridgeP. D.GillesA.PowerD. M. (2005). The secretin G protein-coupled receptor family: teleost receptors.J. Mol. Endocrinol.34753–765.
12
CardosoJ. C. R.VieiraF. A.GomesA. S.PowerD. M. (2010). The serendipitous origin of chordate secretin peptide family members.BMC Evol. Biol.10:135. 10.1186/1471-2148-10-135
13
CleversH. (2006). Wnt/beta-catenin signaling in development and disease.Cell127469–480. 10.1016/j.cell.2006.10.018
14
CvejicS.ZhuZ.FeliceS. J.BermanY.HuangX. Y. (2004). The endogenous ligand stunted of the GPCR Methuselah extends lifespan in Drosophila.Nat. Cell Biol.6540–546. 10.1038/ncb1133
15
de MendozaA.JonesJ. W.FriedrichM. (2016). Methuselah/Methuselah-like G protein-coupled receptors constitute an ancient metazoan gene family.Sci Rep.6:21801. 10.1038/srep21801
16
DiY.-Q.HanX.-L.KangX.-L.WangD.ChenC.-H.WangJ.-X.et al (2020). Autophagy triggers CTSD (cathepsin D) maturation and localization inside cells to promote apoptosis.Autophagy171170–1192. 10.1080/15548627.2020.1752497
17
DoddG. T.LuckmanS. M. (2013). Physiological Roles of GPR10 and PrRP signaling.Front. Endocrinol. (Lausanne)4:20. 10.3389/fendo.2013.00020
18
FanY.SunP.WangY.HeX.DengX.ChenX.et al (2010). The G protein-coupled receptors in the silkworm, Bombyx mori.Insect Biochem. Mol. Biol.40581–591. 10.1016/j.ibmb.2010.05.005
19
FireA.XuS.MontgomeryM. K.KostasS. A.DriverS. E.MelloC. C. (1998). Potent and specific genetic interference by double-stranded RNA in Caenorhabditis elegans.Nature391806–811. 10.1038/35888
20
FlockT.HauserA. S.LundN.GloriamD. E.BalajiS.BabuM. M. (2017). Selectivity determinants of GPCR-G-protein binding.Nature545317–322. 10.1038/nature22070
21
FredrikssonR.SchiothH. B. (2005). The repertoire of G protein-coupled receptors in fully sequenced genomes.Mol. Pharmacol.671414–1425. 10.1124/mol.104.009001
22
FriedrichM.JonesJ. W. (2016). Gene ages, nomenclatures, and functional diversification of the methuselah/methuselah-like GPCR family in Drosophila and Tribolium.J. Exp. Zool. B Mol. Dev. Evol.326453–463. 10.1002/jez.b.22721
23
GaoX. K.ZhangS.LuoJ. Y.WangC. Y.LuL. M.ZhangL. J.et al (2017). Identification and validation of reference genes for gene expression analysis in Aphidius gifuensis (Hymenoptera: Aphidiidae).PLoS One12:e0188477. 10.1371/journal.pone.0188477
24
HakimR. S.BaldwinK.SmaggheG. (2010). Regulation of midgut growth, development, and metamorphosis.Annu. Rev. Entomol.55593–608. 10.1146/annurev-ento-112408-085450
25
HanlonC. D.AndrewD. J. (2015). Outside-in signaling–a brief review of GPCR signaling with a focus on the Drosophila GPCR family.J. Cell Sci.1283533–3542. 10.1242/jcs.175158
26
HauserA. S.AttwoodM. M.Rask-AndersenM.SchiothH. B.GloriamD. E. (2017). Trends in GPCR drug discovery: new agents, targets and indications.Nat. Rev. Drug Discov.16829–842. 10.1038/nrd.2017.178
27
HillC. A.FoxA. N.PittsR. J.KentL. B.TanP. L.ChrystalM. A.et al (2002). G protein-coupled receptors in Anopheles gambiae.Science298176–178. 10.1126/science.1076196
28
JaW. W.CarvalhoG. B.MadrigalM.RobertsR. W.BenzerS. (2009). The Drosophila G protein-coupled receptor, Methuselah, exhibits a promiscuous response to peptides.Protein Sci.182203–2208. 10.1002/pro.221
29
KangX. L.LiY. X.LiY. L.WangJ. X.ZhaoX. F. (2021). The homotetramerization of a GPCR transmits the 20-hydroxyecdysone signal and increases its entry into cells for insect metamorphosis.Development148:dev196667. 10.1242/dev.196667
30
KangX. L.ZhangJ. Y.WangD.ZhaoY. M.HanX. L.WangJ. X.et al (2019). The steroid hormone 20-hydroxyecdysone binds to dopamine receptor to repress lepidopteran insect feeding and promote pupation.PLoS Genet.15:e1008331. 10.1371/journal.pgen.1008331
31
KayukawaT.MinakuchiC.NamikiT.TogawaT.YoshiyamaM.KamimuraM.et al (2012). Transcriptional regulation of juvenile hormone-mediated induction of Kruppel homolog 1, a repressor of insect metamorphosis.Proc. Natl. Acad. Sci. U.S.A.10911729–11734. 10.1073/pnas.1204951109
32
KingT. D.ZhangW.SutoM. J.LiY. H. (2012). Frizzled7 as an emerging target for cancer therapy.Cell Signal.24846–851.
33
LatorracaN. R.VenkatakrishnanA. J.DrorR. O. (2017). GPCR dynamics: structures in motion.Chem. Rev.117139–155. 10.1021/acs.chemrev.6b00177
34
LinY. J.SeroudeL.BenzerS. (1998). Extended life-span and stress resistance in the Drosophila mutant methuselah.Science282943–946. 10.1126/science.282.5390.943
35
LivakK. J.SchmittgenT. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method.Methods25402–408. 10.1006/meth.2001.1262
36
LuY.YuanM.GaoX.KangT.ZhanS.WanH.et al (2013). Identification and validation of reference genes for gene expression analysis using quantitative PCR in Spodoptera litura (Lepidoptera: Noctuidae).PLoS One8:e68059. 10.1371/journal.pone.0068059
37
MaggioliniM.PicardD. (2010). The unfolding stories of GPR30, a new membrane-bound estrogen receptor.J. Endocrinol.204105–114. 10.1677/JOE-09-0242
38
MertensI.HussonS. J.JanssenT.LindemansM.SchoofsL. (2007). PACAP and PDF signaling in the regulation of mammalian and insect circadian rhythms.Peptides281775–1783. 10.1016/j.peptides.2007.05.005
39
MirthC. K.TrumanJ. W.RiddifordL. M. (2009). The ecdysone receptor controls the post-critical weight switch to nutrition-independent differentiation in Drosophila wing imaginal discs.Development1362345–2353. 10.1242/dev.032672
40
NygaardR.ZouY.DrorR. O.MildorfT. J.ArlowD. H.ManglikA.et al (2013). The dynamic process of beta(2)-adrenergic receptor activation.Cell152532–542. 10.1016/j.cell.2013.01.008
41
PalliS. R.LaddT. R.RetnakaranA. (1997). Cloning and characterization of a new isoform of Choristoneura hormone receptor 3 from the spruce budworm.Arch. Insect Biochem. Physiol.3533–44. 10.1002/(SICI)1520-6327199735:1/2<33::AID-ARCH4>3.0.CO;2-Y
42
PatelM. V.HallalD. A.JonesJ. W.BronnerD. N.ZeinR.CaravasJ.et al (2012). Dramatic expansion and developmental expression diversification of the methuselah gene family during recent Drosophila evolution.J. Exp. Zool. B Mol. Dev. Evol.318368–387. 10.1002/jez.b.22453
43
RasmussenS. G.DeVreeB. T.ZouY.KruseA. C.ChungK. Y.KobilkaT. S.et al (2011). Crystal structure of the beta2 adrenergic receptor-Gs protein complex.Nature477549–555. 10.1038/nature10361
44
ReaganJ. D. (1994). Expression cloning of an insect diuretic hormone receptor. A member of the calcitonin/secretin receptor family.J. Biol. Chem.2699–12.
45
SadowskiM. I.ParishJ. H. (2003). Automated generation and refinement of protein signatures: case study with G protein-coupled receptors.Bioinformatics19727–734. 10.1093/bioinformatics/btg075
46
SrivastavaD. P.YuE. J.KennedyK.ChatwinH.RealeV.HamonM.et al (2005). Rapid, nongenomic responses to ecdysteroids and catecholamines mediated by a novel Drosophila G protein-coupled receptor.J. Neurosci.256145–6155. 10.1523/JNEUROSCI.1005-05.2005
47
StrasserA.WittmannH. J.SeifertR. (2017). Binding kinetics and pathways of ligands to GPCRs.Trends Pharmacol. Sci.38717–732. 10.1016/j.tips.2017.05.005
48
TamuraK.StecherG.PetersonD.FilipskiA.KumarS. (2013). MEGA6: molecular evolutionary genetics analysis version 6.0.Mol. Biol. Evol.302725–2729. 10.1093/molbev/mst197
49
Van der HorstD. J.Van MarrewijkW. J.DiederenJ. H. (2001). Adipokinetic hormones of insect: release, signal transduction, and responses.Int. Rev. Cytol.211179–240. 10.1016/s0074-7696(01)11019-3
50
Vanden BroeckJ.TorfsH.PoelsJ.Van PoyerW.SwinnenE.FerketK.et al (1999). Tachykinin-like peptides and their receptors – A review.Ann. N. Y. Acad. Sci.897374–387. 10.1111/j.1749-6632.1999.tb07907.x
51
WangD.ZhaoW. L.CaiM. J.WangJ. X.ZhaoX. F. (2015). G protein-coupled receptor controls steroid hormone signaling in cell membrane.Sci. Rep.5:8675. 10.1038/srep08675
52
WangJ. L.JiangX. J.WangQ.HouL. J.XuD. W.WangJ. X.et al (2007). Identification and expression profile of a putative basement membrane protein gene in the midgut of Helicoverpa armigera.BMC Dev. Biol.7:76. 10.1186/1471-213X-7-76
53
WestA. P.Jr.LlamasL. L.SnowP. M.BenzerS.BjorkmanP. J. (2001). Crystal structure of the ectodomain of Methuselah, a Drosophila G protein-coupled receptor associated with extended lifespan.Proc. Natl. Acad. Sci. U.S.A.983744–3749. 10.1073/pnas.051625298
54
WuF.ZhangY.SunB.McMahonA. P.WangY. (2017). Hedgehog signaling: from basic biology to cancer therapy.Cell Chem. Biol.24252–280. 10.1016/j.chembiol.2017.02.010
55
WuK. M.LuY. H.FengH. Q.JiangY. Y.ZhaoJ. Z. (2008). Suppression of cotton bollworm in multiple crops in china in areas with Bt toxin-containing cotton.Science3211676–1678. 10.1126/science.1160550
56
XuJ.WangX. F.ChenP.LiuF. T.ZhengS. C.YeH.et al (2016). RNA interference in moths: mechanisms, applications, and progress.Genes (Basel)7:88. 10.3390/genes7100088
57
ZamoreP.TuschlT.SharpP.BartelD. (2000). RNAi: double-stranded RNA directs the ATP-dependent cleavage of mRNA at 21 to 23 nucleotide intervals.Cell10125–33. 10.1016/s0092-8674(00)80620-0
58
ZhaoX. F. (2020). G protein-coupled receptors function as cell membrane receptors for the steroid hormone 20-hydroxyecdysone.Cell Commun. Signal.18:146.
59
ZhaoX. F.WangJ. X.WangY. C. (1998). Purification and characterization of a cysteine proteinase from eggs of the cotton boll worm, Helicoverpa armigera.Insect Biochem. Mol. Biol.28259–264. 10.1016/S0965-1748(98)00015-0
Summary
Keywords
genome, GPCR, 20-hydroxyecdysone (20E), forkhead box O, Pten
Citation
Li Y-L, Li Y-X, Wang X-P, Kang X-L, Guo K-Q, Dong D-J, Wang J-X and Zhao X-F (2021) Identification and Functional Analysis of G Protein-Coupled Receptors in 20-Hydroxyecdysone Signaling From the Helicoverpa armigera Genome. Front. Cell Dev. Biol. 9:753787. doi: 10.3389/fcell.2021.753787
Received
05 August 2021
Accepted
30 September 2021
Published
26 October 2021
Volume
9 - 2021
Edited by
Sameer Mohammad, King Abdullah International Medical Research Center (KAIMRC), Saudi Arabia
Reviewed by
Wen Liu, Huazhong Agricultural University, China; Zhen Li, China Agricultural University, China
Updates
Copyright
© 2021 Li, Li, Wang, Kang, Guo, Dong, Wang and Zhao.
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: Xiao-Fan Zhao, xfzhao@sdu.edu.cn
This article was submitted to Signaling, a section of the journal Frontiers in Cell and Developmental Biology
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.