ORIGINAL RESEARCH article

Front. Genet., 10 June 2020

Sec. Epigenetics and Genome Architecture

Volume 11 - 2020 | https://doi.org/10.3389/fgene.2020.00609

Involvement of Two Paralogous Methoprene-Tolerant Genes in the Regulation of Vitellogenin and Vitellogenin Receptor Expression in the Rice Stem Borer, Chilo suppressalis

  • LM

    Lijun Miao 1

  • NZ

    Nan Zhang 1

  • HJ

    Heng Jiang 1

  • FD

    Fan Dong 1

  • XY

    Xuemei Yang 1

  • XX

    Xin Xu 1

  • KQ

    Kun Qian 1

  • XM

    Xiangkun Meng 1

  • JW

    Jianjun Wang 1,2*

  • 1. College of Horticulture and Plant Protection, Yangzhou University, Yangzhou, China

  • 2. Joint International Research Laboratory of Agriculture and Agri-Product Safety of the Ministry of Education, Yangzhou University, Yangzhou, China

Abstract

Besides the function of preventing metamorphosis in insects, the juvenile hormone (JH) plays a role in female reproduction; however, the underlying mechanism is largely unknown. The methoprene-tolerant (Met) protein belongs to a family of basic helix-loop-helix–Per-Arnt-Sim (bHLH-PAS) transcription factors and functions as the JH intracellular receptor. In this study, two full length cDNAs encoding Met (CsMet1 and CsMet2) were isolated from the rice stem borer, Chilo suppressalis. Structural analysis revealed that both CsMet1 and CsMet2 exhibited typical bHLH, PAS-A, PAS-B, and PAC (PAS C terminal motif) domains. Comparative analysis of transcript level using reverse transcription-quantitative PCR (RT-qPCR) revealed that CsMet1 was predominant in almost all examined developmental stages and tissues. Treatment with methoprene in vivo induces the transcription of both CsMet1 and CsMet2. Notably, injection of dsCsMet1 and dsCsMet2 suppressed the expression levels of vitellogenin (CsVg) and Vg receptor (CsVgR). These findings revealed the potential JH signaling mechanism regulating C. suppressalis reproduction, and provided evidence that RNAi-mediated knockdown of Met holds great potential as a control strategy of C. suppressalis.

Introduction

The rice stem borer, Chilo suppressalis (Walker) (Lepidoptera: Crambidae), is one of the most serious rice pests in Asia, Middle East, and southern Europe, and causes large crop losses through feeding on the stems of rice. In China, the change of cultivation patterns has led to the frequent outbreaks of C. suppressalis in recent years (). To date, spraying chemical insecticides remains the primary strategy for controlling C. suppressalis. However, intensive use of insecticides has driven C. suppressalis to develop resistance to a wide range of insecticides (; ; ). Hence, the development of RNA interference (RNAi)-mediated disruption of reproduction represents an alternative control strategy of C. suppressalis ().

Insect reproduction is intricately regulated by three hormones, including the juvenile hormones (JHs), ecdysteroids, and insulin (; ; ). JH exerts its function through its intracellular receptor methoprene-tolerant (Met), a member of the family of the basic helix-loop-helix (bHLH)-Per-Arnt-Sim (PAS) transcription factors (; , ). Since the first characterization of Met in Drosophila melanogaster (), Met has been identified from a broad range of insect species, such as Aedes aegypti (), Tribolium castaneum (), Nilaparvata lugens (), and Helicoverpa armigera (). Interestingly, two Met paralogs, Met and Germ cell-expressed (Gce), are present across 12 Drosophila species (, ). While mosquitoes and beetles possess only a single gene, two Met genes have also been identified in several Lepidopteran insects, including Danaus plexippus (), Operophtera brumata (), and Bombyx mori (; ; ). However, the role of two Met genes in the reproduction of Lepidopteran insects remains largely unknown.

The reproductive success of insects depends on vitellogenin (Vg) biosynthesis and uptake of Vg into developing oocytes mediated by vitellogenin receptor (VgR). Recently, we have characterized the C. suppressalis CsVg and CsVgR at the molecular levels (; ). In this study, full-length cDNAs of two paralogous Met genes (named as CsMet1 and CsMet2) were isolated from C. suppressalis. We report the structural features and temporal–spatial expression patterns of CsMet1 and CsMet2. Furthermore, RNAi was employed to reveal the role of CsMet1 and CsMet2 in the regulation of CsVg and CsVgR.

Materials and Methods

Insects Rearing and Sampling

Chilo suppressalis larvae was collected from rice stubbles of Yangzhou (32.39°N, 119.42°E) in 2013, and reared on an artificial diet in an incubator at 28 ± 1°C, 70 ± 5% RH, and a 16-h light/8-h dark photoperiod.

To examine the developmental expression profiles of CsMet, individuals were collected from larvae (third, fourth, fifth, and sixth instar), pupae at intervals of 2 days from pupation, and adults at intervals of 12 h from eclosion. For tissue expression analysis, tissues (including head, epidermis, midgut, hemolymph, and fat body) were dissected from the fifth-instar larva. All the samples were frozen immediately in liquid nitrogen and stored at -80°C until RNA isolation. Each experiment was performed with three biological replicates containing three to 10 individuals.

RNA Isolation, RT-PCR, and RACE

Total RNA was extracted using Trizol reagent (Invitrogen, Carlsbad, CA, United States). First-strand cDNA was synthesized from total RNA using the PrimescriptTM First-Strand cDNA Synthesis Kit (TaKaRa, Dalian, China) according to the manufacturers’ instruction.

The amino acid (aa) sequences of B. mori BmMet1 (GenBank: NP001108458) and BmMet2 (GenBank: BAJ05086) were searched against the transcriptome database of C. suppressalis in Insect Base1, and specific primer pairs were designed based on the sequences of putative transcripts (Table 1). PCR reactions were performed with LA TaqTM DNA polymerase (TaKaRa, Dalian, China). To obtain full length cDNA sequences of CsMet1 and CsMet2, 5′-RACE and 3′-RACE were conducted using the SMART RACE cDNA Amplification Kit (Clontech, Mountain View, CA, United States), following the manufacturers’ instructions. Gene specific primers (GSPs) used for RACE are listed in Table 1.

TABLE 1

Primer nameSequence (5′ to 3′)Description
W988Met1 FATGACATCATTGACTGGAGCCsMet1 RT-PCR
W989Met1 RGGATTACAGGATTTCAGTTTCTGA
W71Met1 RCAACAGGCAGTCAGTCACCAAGTCsMet1 5′-RACE
W72Met1 RTGACATCATTGACTGGAGCCACTG
W73Met1 FTGTTGGTGTAGATTATGGGCGACGCsMet1 3′-RACE
W74Met1 FAAGTGGTGCAAGAAACTGGTGTCC
W990Met2 FAGAGAGATTCGAAACAAAGCGCsMet2 RT-PCR
W991Met2 RTCGTTTGTACCAACACTGTC
W75Met2 RCCAGCATTTCGGCGACCTTGTTCTTCsMet2 5′-RACE
W76Met2 RCCAGTTCACCGATGGATTGGTTCAG
W77Met2 FGTGTTCGTGGGCATTGTCCGCTTGGCsMet2 3′-RACE
W78Met2 FTTTAGTCGGTGAGTCCTGCTATCGT
EF1-α FTGAACCCCCATACAGCGAATCCRT-qPCR
EF1-α RTCTCCGTGCCAACCAGAAATAGG
W465Met1 qFTGGCTTCCTCGAGATTGACART-qPCR
W466Met1 qRTCCTGATGCTACCCCAGATG
W469Met2 qFCAATCCATCGGTGAACTGGCRT-qPCR
W470Met2 qRGTTGAGTATGGACAGCAGCG
W84VgR FAGCCACTTCCCTACCTCCTART-qPCR
W85VgR RTAAGGCATTGGGGACTCGTT
W976Vg FAGCTCAGTCCGCTAAATGGART-qPCR
W977Vg RGCCCAGTTCGTGGTGTCTAT
W463Met1 iFTAATACGACTCACTATAGGGTCTGAC ATAGTGCACGCTCCCsMet1 RNAi
W464Met1 iRTAATACGACTCACTATAGGGTGTGCG TACCCTTCTGACAC
W467Met2 iFTAATACGACTCACTATAGGGCAGGGG GCTCATTGTGGTAGCsMet2 RNAi
W468Met2 iRTAATACGACTCACTATAGGGGTGCCG CGTTCTATACTCCA

Oligonucleotide primers used for RT-PCR, RACE, RT-qPCR, and RNAi.

Molecular Cloning and Sequence Analysis

RT-PCR and RACE products were subcloned into the pMD18−T vector (TaKaRa) and sequenced. The molecular mass and isoelectric point (pI) of the deduced protein sequences were predicted by using the online ExPASy proteomics server2. Conserved domains were predicted by NCBI conserved domain search tool3 or by alignment to other published insect Mets. A phylogenetic tree was constructed with MEGA 7.0 using the neighbor-joining method with a p-distance model and a pairwise deletion of gaps (). The reliability of the NJ tree topology was statistically evaluated by bootstrap analysis with 1000 replicates.

Reverse Transcription-Quantitative PCR (RT-qPCR)

Reverse transcription-quantitative qPCR reactions were performed on the Bio-Rad CFX-96TM Real-time PCR system using TB GreenTM Premix Ex TaqTM (Takara, Dalian, China) following the manufacturers’ instructions. GSPs used are listed in Table 1. The stably expressed gene encoding EF1-α was used as a reference gene (; ). The mRNA levels were normalized to reference gene with the 2–ΔΔCT method (). The means and standard errors for each time point were obtained from the average of three biologically independent samples.

Methoprene Treatment

To study effects of JH on expression of CsMet, newly emerged adult females were injected with 1 μL JH analog methoprene (3 μg/μL, Sigma–Aldrich, St. Louis, MO, United States) or same volume of acetone as control. Insects were collected for RT-qPCR analysis at 1, 6, 12, and 24 h after treatment, respectively, and experiments contained three biological replications.

RNA Interference

Double-strand RNAs (dsRNAs) against CsMet1 or CsMet2 were synthesized using TranscriptAidTM T7 High Yield Transcription Kit (Thermo Fisher Scientific, Waltham, MA, United States). GSPs for dsRNA synthesis are listed in Table 1. Using a Nanoliter 2010 injector system (WPI, Sarasota, FL, United States), 1 μL solution of dsRNA (3 μg/μL) was injected into the abdomen of 6-day-old female pupae under a stereomicroscope. Equal dose of dsRNA for enhanced green fluorescent protein (EGFP) was injected as a control. A total of 40 pupae were injected per treatment, and each treatment was performed in triplicate. The treatment was repeated in newly emerged female adults (within 24 h), and insects were collected at 24 and 48 h after injection, respectively, for expression analysis of CsMet1, CsMet2, CsVg, and CsVgR by RT-qPCR.

Data Analysis

The statistical analysis was done using graphpad.prism.6 by one-way analysis of variance, followed by student’s t-test. All data are presented as the mean ± SE.

Results

Sequence and Structural Analysis of CsMet1 and CsMet2

The full length CsMet1 cDNA (GenBank accession number MN906993) was 2673 bp, with a 460 bp 5’-terminal untranslated region (UTR), a 1587 bp open reading frame (ORF), and a 626 bp 3′-UTR. The deduced CsMet1 protein contained 528 aas with a predicted mass and pI of 60.34 kDa and 8.35, respectively. The 2824 bp CsMet2 cDNA (GenBank accession number MN906994) contained an 85 bp 5′-UTR, a 2598 bp ORF encoding an 865 aa residue protein with a molecular mass of 97.37 kDa and a pI of 6.68, and a 141 bp 3′-UTR.

Structural analysis showed that the deduced protein of CsMet1 and CsMet2 had conserved domains of bHLH-PAS protein family including bHLH (DNA binding and dimerization regions), PAS-A (dimerization region), PAS-B (ligand binding and dimerization region), and PAC (PAS C terminal motif) domains (dimerization region) (Figure 1). Alignment of aa sequences demonstrated that CsMet1 shared identity with other insect Met orthologs including B. mori BmMet1 (60.98%) and BmMet2 (21.88%), D. melanogaster DmMet (27.49%), and T. castaneum TcMet (28.30%). CsMet2 shared 20.55, 44.19, 17.19, and 22.90% identity with BmMet1, BmMet2, DmMet, and TcMet, respectively. The aa identity between CsMet1 and CsMet2 was 21.71% (Figure 1). Phylogenetic analysis revealed that CsMet1 and CsMet2 were clustered into Lepidopteran group 1 and group 2 Met, respectively (Figure 2).

FIGURE 1

FIGURE 2

Temporal and Spatial Expression of CsMet1 and CsMet2

The temporal expression patterns of CsMet1 and CsMet2 were detected from third instar to 3-day-old female adults and determined by RT-qPCR analysis. The results showed that the transcription level of CsMet1 and CsMet2 was detectable during all selected developmental stages. Specifically, the expression of CsMet1 was relatively stable in the larval stage, sharply decreased in the prepupal stage, maintained at a low level during the pupal stage, and reached a peak in 12-h-old female adults (Figure 3A). The expression of CsMet2 gradually increased from fourth instar to 5-day-old pupae. In the adult stage, the highest and lowest expression levels of CsMet2 were observed in 48-h-old and 72-h-old female adults, respectively (Figure 3B).

FIGURE 3

Analysis of the spatial expression profiles of CsMet1 and CsMet2 in fifth-instar larva showed that CsMet1 was highly expressed in head, hemocytes, and midgut, while relatively low expression was observed in epidermis, fat body, and malpighian tube (Figure 3A). The highest expression level of CsMet2 was observed in hemocytes, followed by head and Malpighian tube, whereas low expression was observed in epidermis, fat body, and midgut (Figure 3B).

Effect of JHA Treatment

To examine whether CsMet1 and CsMet2 were regulated by JH, JHA methoprene was injected into the newly emerged female adults. Compared with control, the results showed that CsMet1 expression was upregulated by 2.74, 2.41, and 1.63 times at 1, 6, and 12 h after injection, respectively (Figure 4A). Similarly, transcription level of CsMet2 was increased by 1.95, 1.45, and 3.68 times at 1, 6, and 1 2 h post-treatment, respectively (Figure 4B).

FIGURE 4

Effects of RNAi-Mediated Knockdown of CsMet1 and CsMet2

To explore the role of CsMet1 and CsMet2 in the regulation of CsVg and CsVgR, female adults were collected at 24 and 48 h after the second dsRNA injection, respectively. The transcript level of CsMet1 in dsCsMet1-injected female adults was suppressed by 78.74% (24 h) and 70.32% (48 h) compared with dsEGFP-injected insects, whereas the expression of CsMet2 was significantly reduced by 94.69% (24 h) and 92.79% (48 h) (Figure 5A). Meanwhile, RNAi-mediated suppression of CsMet1 and CsMet2 decreased CsVg expression by 70.85 and 44.37% at 48 h post-treatment, respectively (Figure 5B). CsVgR expression in the dsCsMet1-treated and dsCsMet2-treated group decreased significantly by 52.91 and 62.98% at 24 h post-injection, respectively (Figure 5C). Furthermore, phenotypic observation revealed that silencing of both CsMet1 and CsMet2 resulted in suppressed ovarian development with decreased number of mature follicles (Figure 6).

FIGURE 5

FIGURE 6

Discussion

The important roles played by JH in insect development and reproduction prompted the study on JH signaling, and the breakthrough had been the identification of transcription factor Met as an intracellular receptor for JH (). Binding of JH to Met triggers dimerization of Met with another bHLH-PAS protein Taiman (Tai) to form a functional complex, which interacts with JH response elements (JHREs) of target genes (). While the Met null allele was expected to result in a lethal phenotype, the paralogous Met gene in D. melanogaster, Gce, ensured survival of the met null mutants (, ). Similarly, two Met genes were also recently identified in several Lepidopteran insects (; ; ; ; ). However, only a single Met gene was found in those from other insect orders including A. aegypti (Dipteran, ), T. castaneum (Coleopteran, ), Blattella germanica (Blattarian, ), N. lugens (Hemipteran, ), Sitodiplosis mosellana (Dipteran, ), and so on. In this study, we isolated two paralogous Met genes, CsMet1 and CsMet2, from C. suppressalis. Both proteins showed conserved domain organization of bHLH-PAS protein family (; ; ). Phylogenetic analysis revealed that CsMet1 and CsMet2 were clustered into Lepidopteran group 1 and group 2 Met, respectively, providing further evidence of the occurrence of Met gene duplication in Lepidopteran insects.

The temporal and spatial expression patterns of a gene may provide hints as to its function. In this study, we found a marked reduction in expression levels of CsMet1 in the prepupae stage, supporting its involvement in larval–pupal metamorphosis. Similar results were also observed in HaMet1 of H. armigera () and SmMet of S. mosellana (), and was in agreement with the fact that JH exerts its anti-metamorphic effect through its receptor Met (; ). CsMet2 transcript was hardly detected during the larval stage, and its expression reached the peak in 48-h-old female adults, suggesting that CsMet2 was involved in JH signaling in adults, as had been reported for BmMet2 in B. mori (). For tissue expression, CsMet1 was highly expressed in larval head, midgut, and hemocytes, whereas CsMet2 showed the highest expression in larval hemocytes. Interestingly, the temporal and spatial expression analysis revealed that CsMet1 was dominant at almost all developmental stages and in all tissues examined. This result was somewhat different from the observation in B. mori that BmMet2 transcript was more abundant than BmMet1 transcript in the adult stage (). Whether CsMet1 played a more important role in C. suppressalis needed further study in the future.

Juvenile hormone plays a crucial role in insect reproduction, but its molecular mode of action only became clear recently. It has been reported in H. armigera, Diploptera punctata, N. lugens, and S. mosellana that Met is significantly upregulated by JH III or JHA methoprene (; ; ; ). Our hormone treatment results also indicated that CsMet1 and CsMet2 expression was upregulated after JHA methoprene treatment. Further RNAi analysis revealed that knockdown of CsMet1 and CsMet2 resulted in significant decline of both CsVg and CsVgR expression in female adults, leading to decreased number of mature follicles in ovaries. Similarly, expression of Vg is regulated by JH via Met in A. aegypti, H. armigera, D. punctata, N. lugens, and Bactrocera dorsalis (; ; ; ; ). Downregulation of VgR transcription as well as suppressed ovarian development were also observed in dsMet-treated H. armigera and Colaphellus bowringi (; ). Whether expression of Vg and VgR were directly regulated by Met was not well known. It had been recently reported that knockdown of Krüppel-homolog 1 (Kr-h1), a direct target of JH-Met signaling, also decreased VgR expression in C. bowringi (). On the other hand, Met is also a repressor of 20E pathway gene expression. It is likely that Met regulated the relative genes of 20E pathway (), which contributes to vitellogenis to some extent. Further study needs to be carried out to confirm these presumptions.

Statements

Data availability statement

All datasets generated for this study are included in the article.

Author contributions

JW conceived and designed the study. LM, NZ, HJ, FD, XY, and XX performed the experiments. LM, JW, KQ, and XM wrote the manuscript. All authors reviewed and approved this submission.

Funding

This work was supported by the Special Fund for Agro-scientific Research in the Public Interest under Grant No. 201303017 and National Natural Science Foundation of China (31701807).

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.

References

Summary

Keywords

Chilo suppressalis, juvenile hormone, methoprene-tolerant, vitellogenin, vitellogenin receptor, RNAi

Citation

Miao L, Zhang N, Jiang H, Dong F, Yang X, Xu X, Qian K, Meng X and Wang J (2020) Involvement of Two Paralogous Methoprene-Tolerant Genes in the Regulation of Vitellogenin and Vitellogenin Receptor Expression in the Rice Stem Borer, Chilo suppressalis. Front. Genet. 11:609. doi: 10.3389/fgene.2020.00609

Received

20 February 2020

Accepted

19 May 2020

Published

10 June 2020

Volume

11 - 2020

Edited by

Fei Li, Zhejiang University, China

Reviewed by

Weihua Ma, Huazhong Agricultural University, China; Wei Dou, Southwest University, China

Updates

Copyright

*Correspondence: Jianjun Wang,

This article was submitted to Epigenomics and Epigenetics, a section of the journal Frontiers in Genetics

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics