Selection of housekeeping genes and demonstration of RNAi in cotton leafhopper, Amrasca biguttula biguttula (Ishida)

Amrasca biguttula biguttula (Ishida) commonly known as cotton leafhopper is a severe pest of cotton and okra. Not much is known on this insect at molecular level due to lack of genomic and transcriptomic data. To prepare for functional genomic studies in this insect, we evaluated 15 common housekeeping genes (Tub, B-Tub, EF alpha, GADPH, UbiCF, RP13, Ubiq, G3PD, VATPase, Actin, 18s, 28s, TATA, ETF, SOD and Cytolytic actin) during different developmental stages and under starvation stress. We selected early (1st and 2nd), late (3rd and 4th) stage nymphs and adults for identification of stable housekeeping genes using geNorm, NormFinder, BestKeeper and RefFinder software. Based on the different algorithms, RP13 and VATPase are identified as the most suitable reference genes for quantification of gene expression by reverse transcriptase quantitative PCR (RT-qPCR). Based on RefFinder which comprehended the results of three algorithms, RP13 in adults, Tubulin (Tub) in late nymphs, 28S in early nymph and UbiCF under starvation stress were identified as the most stable genes. We also developed methods for feeding double-stranded RNA (dsRNA) incorporated in the diet. Feeding dsRNA targeting Snf7, IAP, AQP1, and VATPase caused 56.17–77.12% knockdown of targeted genes compared to control and 16 to 48% mortality of treated insects when compared to control.


Introduction
Gene expression studies are indispensable for molecular biology research. The knowledge on gene expression helps to better understand its regulation and functions. The availability of huge sequence data in the form of transcriptomes and genomes of several organisms can be utilized to understand the transcription of the gene(s) [1]. However, these studies have become more accurate and robust after the development of gene expression quantitation method, reverse transcriptase quantitative PCR (RT-qPCR). This method is highly sensitive, reproducible and accurate to a level that it can identify even minute variations, which are frequently undetected. The RT-qPCR data are influenced by many factors, which include quality and a1111111111 a1111111111 a1111111111 a1111111111 a1111111111 quantity of the starting material, RNA extraction, cDNA synthesis, and other laboratory procedures. Even pipetting errors and reverse transcription efficiency can influence the Ct values significantly [2,3]. Thus normalization is a prerequisite in gene expression studies as it limits variability by comparing target gene expression with housekeeping genes (HKGs). Normalization is based on the assumption that the expression of the HKGs is stable across various biotic and abiotic stresses and treatments. Recent research indicates that a condition-specific reference gene needs to be identified for accurate measurements of gene expression [4]. Likewise, it is also evident that a single reference gene is not desirable for the wider experimental regime [5,6]. The use of a single reference gene can generate up to a 20-fold error in the expression data [7]. In most of the expression studies, actin is taken as a universal HKG [8] or the HKGs validated for certain systems are being directly applied to other without an appropriate validation of their stability in that particular system. To address this, different software programs have been developed to choose HKGs that are most suited for normalization [7][8][9][10]. Thus, it is necessary to choose the most suitable genes for normalization from a panel of candidate genes in a dedicated set of biological samples from a particular organism. Although RT-qPCR has been widely used for detecting gene expression in insects but there is yet no suitable HKG and stable gene quantification system for the cotton leafhopper. In addition, there is no report of RNAi (RNA interference) in this insect. We demonstrated feeding RNAi in this insect for few genes such as Aquaporins (AQPs), inhibitor-of-apoptosis (IAP), VATPases and Snf7. Aquaporins (AQPs) belong to the family of the major intrinsic proteins (MIP), which are integral membrane channel proteins in most living organisms and facilitate mass transfer of water and sometimes other substrates across cell membranes [11][12][13][14]. In insects, the inhibitor-of-apoptosis (IAP) proteins are a family of conserved survival factors that determine cell fate during development, stress, tumorigenesis, and infection by regulating the phenomenon of apoptosis [15][16][17][18][19]. VATPases are highly conserved and ubiquitous proton pumps which acidify specific organelles such as lysosomes, endosomes or secretory vesicles in every eukaryotic cell. These are present in the plasma membrane of different animal cell types where they are involved either in membrane energization or pH homeostasis [20][21][22]. Snf7 functions as a component of the ESCRT (Endosomal Sorting Complex Required for Transport) pathway which plays an essential part in cellular housekeeping by internalization, transport, sorting and lysosomal degradation of transmembrane proteins [23]. It has been shown to be involved in sorting of transmembrane proteins either through recycling to plasma membrane or routed to lysosomal degradation through the endosomal-autophagic pathway in many organisms [24][25][26][27][28].
Amrasca biguttula biguttula (Ishida) (Hemiptera: Cicadellidae), commonly known as cotton leafhopper/jassid, is among the most economically important sucking pests of cotton. Leafhopper infested tender leaves become small, crinkled, yellow and the margin of the leaves starts curling downwards. In the case of severe infestation, leaves get a bronze or brick red color which is typical "hopperburn" symptom. The leaves dry up and are shed, and the growth of the plant is retarded. Apart from drying of leaves, the punctures of leafhoppers induce the shortening of the internodes, which contribute to reduced plant vigor and yield [29]. The pest has the potential of causing 25-45% loss in seed cotton yield [30,31]. Besides cotton, it also infests many crops in malvaceous and solanaceous families. The mainstay of leafhopper control shifted to neonicotinoids after the development of resistance against most of the pyrethroids and organophosphates [32]. Bt-cotton seeds treated with neonicotinoids protect against cotton leafhopper till 40-50 days after sowing. However, there have been reports of the leafhoppers showing a high level of resistance up to 5,000-fold to imidacloprid and other neonicotinoid insecticides introduced barely decade ago in central India [33]. Besides, the usage of synthetic insecticides is also associated with a number of environmental issues such as insecticide residues in soil and water, and effects on non-targets, etc. The ill effects of pesticides can only be addressed through eco-friendly methods of pest control, but on contrary Bt-toxins (sole candidates of insect-resistant transgenics) are either not available or less effective against sap-sucking pests such as aphids, leafhoppers and whitefly. The present-day RNAi technology has the potential of developing novel insect control strategies. The method relies on feeding/injecting sequence-specific double-stranded RNA (dsRNA) targeted towards downregulation or knockdown of essential genes for causing mortality. The dsRNA targeting such essential genes can help in developing new generation insect resistant transgenic plants. Ironically, in spite of attaining serious pest status on cotton, the molecular genetics of this insect has remained unexplored with information available only on a single gene (mtCOI-taxonomic importance) in the NCBI database. To initiate molecular studies related to gene expression, RNAi, and other future functional genomics studies, selection of right HKGs is a prerequisite. Hence, the studies first time report the suitable HKGs for different developmental stages and under starvation stress in cotton leafhopper. In addition downregulation of targeted genes associated with different physiological functions through dsRNA feeding has been successfully demonstrated in this insect.

Selection of genes and primer design
The selection of genes for identifying stable reference genes in Amrasca biguttula biguttula (Cotton leafhopper) was based on previous studies in other insects [34][35][36]. As per information available, we chose 15 commonly used reference genes such as Tub (Tubulin), β-Tub (Beta Tubulin), EF alpha (elongation factor alpha), GADPH (Glyceraldehyde 3-phosphate dehydrogenase), UbiCF (ubiquitin-conjugating factor), RP13 (ribosomal protein 13), Ubiq (ubiquitin), VATPase (V-type adenosine triphosphatase), Actin, 18s (ribosomal protein 18), 28s (ribosomal protein 28), TATA (TATA-binding protein), ETF (electron transfer flavoprotein), SOD (superoxide dismutase) and Cytolytic actin as shown in Table 1. The sequences of these genes were retrieved from our transcriptome sequence data (unpublished) followed by a BLASTX search of each sequence in NCBI database for re-confirmation of their annotation. The primers were designed using Primer3 [37] software to amplify 100-150 bp region of the respective gene (Table 1). The sequences of all the selected genes have been submitted to NCBI database and are available with Accessions MF101761-MF101776.

Rearing of insect
A. biguttula insects were reared on Gossypium hirsutum variety Ganganagar ageti in the walkin environmental chamber at 65-70% RH, 14:10 h light and dark photoperiod and 27± 2˚C. The plants of the variety were changed from time to time for continuous fresh food supply to the insect and availability of its culture throughout the study period.

dsRNA synthesis and feeding
To test RNAi in A. biguttula biguttula, we chose AQP (aquaporin), IAP, (inhibitor of apoptosis) VATPase and Snf7 (vacuolar sorting protein) genes shown in Table 2. Gene specific primers were designed from the transcriptome sequences using Primer3 open source software using RNA sequences [37]. The dsRNA primers were also designed using Primer3.However, an additional T7 promoter sequence was added to 5' end of both reverse and forward primers as given in Table 2. These primers were used to amplify the template for dsRNA using cDNA which was synthesized using PrimeScript™ 1st strand cDNA Synthesis Kit (Clontech Takara) as per manufacturer's protocol. Amplified fragments of genes were purified using nucleospin PCR cleanup (Macherey-Nagel Nucleospin Gel and PCR Cleanup) and used for the synthesis of dsRNA with T7 RiboMAX™ Express RNAi System (Promega) following manufacturer's protocol. The dsRNA was quantified using BioSpectrometer 1 basic (Eppendorf). The 500ng/μl of dsAQP, dsIAP, dsVATPase, dsSnf7, and dsGFP was mixed with diet [38] as described in Table 3. Three biological replicates of five insects each were released in a 1.5 ml tube (1/4 th cut from the bottom side and covered with muslin cloth) with diet incorporated dsRNA between two layers of stretched parafilm. The insects were observed 48 h post release and the mortality  was compared to dsGFP treated insects. Live insects were collected in Tri-Reagent 1 (Sigma-Aldrich), frozen in liquid nitrogen and processed for RNA extraction.

RNA extraction and cDNA synthesis
For validation of housekeeping genes, insects were categorized into different developmental stages, i.e., early nymph (I and II), late nymph (III and IV) and adult. Three biological replicates of each developmental stage were separated and pooled, i.e., early nymphs 25-30 individuals per pool, late nymphs 15-20 per pool and adult hopper five insects per pool. For expression under starvation stress, late stage nymphs (III and IV) were starved for four h in empty ventilated boxes, and four individuals each were pooled in three different tubes. The total RNA was isolated from each developmental stage, starved leafhopper and dsRNA fed insects using Tri Reagent 1 (Sigma-Aldrich) as per manufacturer's protocol. Isolated RNA was given DNase treatment to remove DNA contamination. The RNA was quantified and checked for its quality and quantity on Eppendorf BioSpectrometer 1 basic. The total RNA (2μg) was used for first-strand cDNA preparation using PrimeScript™ 1st strand cDNA Synthesis Kit (Clontech Takara) as per manufacturer's protocol.

RT-qPCR analysis
qPCR analysis used SYBR1 Premix Ex Taq™ II (TliRNase H Plus) (Clontech) with the Light-Cycler196 (Roche Molecular Systems, Inc.). PCR reactions were performed in triplicates in 10μl volumes using 1μl of 1:10 diluted cDNA, 0.2 ul of 10mM of gene-specific primers per reaction. For absolute quantification, three 10-fold serial dilutions were performed to ensure that the cDNA synthesis reagents did not impair PCR efficiency. Thermal cycling conditions constituted initial denaturation of 95˚C for the 30s followed by 40 cycles of 95˚C of 5s and 60˚C for 10s. Melting curve analysis was done to ensure the specificity and consistency of the amplified product. RT-qPCR efficiency was measured for each gene with the slope of linear regression model and standard curves. Amplification efficiency and correlation coefficients for each primer pair were calculated as described in LightCycler196 users guide. In RNAi studies, the relative expression of genes in biological samples was estimated using the ΔΔCt method [39] normalized with RP13 as HKG and compared with dsGFP.

Stability and statistical analysis for reference genes
We chose four algorithms to determine the stability of genes, i.e., geNorm [7], NormFinder [9], Bestkeeper [10] and RefFinder [40]. The qPCR obtained Ct values from the Lightcycler software (Roche) were transferred to MS-excel for the calculation of linear relative values by comparative Ct method (keeping lowest relative quantity for each gene as 1). These linear relative quantities were used as input data for further analysis of gene stability with geNorm and NormFinder. geNorm calculates the expression stability score (M) by averaging the mean pairwise variation of each HKG. Lower 'M' value indicates the stability of the gene, so the HKGs showing M value > 0.5 were not considered for further normalization studies. NormFinder calculates expression stability of the genes within a group and between the groups. It also determines standard deviation (SD) through advanced analysis. BestKeeper is also a freely available algorithm that directly utilizes the Ct value obtained from the software. It evaluates standard deviation, p-values, index and correlation coefficient of each gene to elucidate the most stable gene. The gene with a lower SD value could serve as a better reference gene. The final assessment was made on the basis of geometric mean calculated for each gene using RefFinder (http://leonxie.esy.es/RefFinder/?type=reference), a web-based tool which uses all three algorithms along with ΔΔCt method for stability analysis of genes. This web-based tool is user-friendly and allows to evaluate and screen the candidate reference genes directly on the basis of Ct values across the samples efficiently. It makes a comprehensive analysis of data obtained from all the algorithms and ranks the candidate genes in decreasing order of their stability. On the grounds of rankings from each program, an appropriate weight is assigned to each candidate gene, and the geometric mean of their weights is calculated for the net final ranking.

Verification of expression of selected genes
Fifteen candidate genes based on their significance in the biological processes were tested for their expression. These included the structure-related genes, Tub, B-Tub, Cytolytic Actin and Actin; ribosomal and protein function genes, EF alpha, 18s, 28s, RP13, TATA, ETF; metabolism-related genes, UbiCF, Ubiq, GAPDH, VATPase and SOD. The initial screening of all targeted genes was done by PCR and the amplified products were checked by agarose gel electrophoresis. These analyses verified that all the targeted genes were expressed in A. biguttula biguttula (data not shown). In order to calculate the amplification efficiency for all fifteen candidate genes, the three-point standard curve was obtained using the Lightcycler software following the PCR amplification with the known concentration of cDNA template. All genes except Ubiq and CytoActin (78.04% and 62.32% respectively) showed significant amplification efficiency (Table 4). So these two genes were eliminated and not taken into account for further analysis. The correlation coefficient (R 2 ) for all genes ranged between 0.92~1.00 (Table 4).
Melt-curve analysis was also performed for amplification specificity, and all genes displayed a single peak. Purified qPCR products were sequenced, and the sequences matched 100% with the target sequences.  Fig 1).

BestKeeper analysis
BestKeeper reveals the best genes on the basis of SD value for each candidate reference gene from raw Ct values. An SD > 1 manifests that the diversification in the expression of a gene within a sample of the same origin is high, thus indicating instability in the expression of that gene. Our data signified that expression was not so stable for many genes, i.e., UbiCF, GAPDH, 28s, and SOD in adult samples as it has shown variations across the samples along with high SD values. Similarly, in late nymph, instability was observed in the genes such as UbiCF, GAPDH, RP13, VATPase and EF based on high SD value (Table 5). In adult hopper, 18s was the most stable gene with lowest SD value of 0.34. Similarly, in late nymph 28s (SD = 0.1), early nymph Tub (SD = 0.08) and under starvation stress TATA (SD = 0.06) were the most stable genes. Ranking of all genes in decreasing order of their stability across all the stages is presented in S1  3).

NormFinder analysis
Based on this program 18S was identified as best gene among all the samples showing stability value~0.102. NormFinder also determines the best combination of two HKGs which can be used simultaneously in a single expression studies. Based on results, the combination of 18s and TATA was highly stable (stability value~0.079). Under starvation stress, UbiCF appeared to be the most stable gene with a stability value of 0.034. The rankings of genes on the basis of their stability values are presented in Table 5 and graphically in

RefFinder analysis
RefFinder assembles the output of all three algorithms described earlier along with ΔΔCt to assign final rankings for all the genes. It interprets the stability of genes by ranking of their geometric mean (Table 6), as well as allocates overall values to the genes across different stages (Table 7). These analyses showed that RP13, TUB, 28s and UbiCF were the most stable genes in the adult, late stage, early stage nymph and starvation stress, respectively. However, RefFinder also revealed RP13, VATPase, 18s in the adult; 28s, RP13, TUB in early nymph; TUB, TATA, Actin in late nymph and UbiCF, B-Tub, G3PD in starvation stress as the most stable genes. Based on overall analysis of all developmental stages, RP13 was identified as the best  stable gene across all stages (Table 7). Thus, this gene was used as HKG for relative expression analysis in dsRNA feeding assays.

Bioassay with gene-specific dsRNA:
The target gene knockdown was confirmed by feeding the late stage nymphs with different dsRNAs using membrane feeding assay (Fig 5). Initial experiments using different doses of dsRNA showed that 500 ng/μl of dsRNA in diet caused knockdown of respective genes when compared to that in dsGFP fed control insects at 48 h post feeding. All the targeted genes showed a significant reduction (p 0.05) in their mRNA levels compared to that in dsGFP fed  Selection of housekeeping genes and RNAi in the cotton leafhopper control insects. The dsRNA feeding caused 4.37, 2.99, 2.49 and 2.28 fold reduction in the levels of Snf7, IAP, AQP and VATPase mRNA, respectively compared to that in control insects ( Fig  6). The percent corrected mortality calculated using Abbott's formula [41] varied between 16-48%; highest in dsSnf7 (48%) fed insects followed by dsAQP (27%), dsVATPase (20%) and dsIAP (16.7%) ( Table 8); (S2 Table).

Discussion
In recent years, it has become clear that selection of an appropriate reference gene or genes is the basic requirement for the success of RT-qPCR [42]. However, in the most of the previously reported expression studies with insects, actin has been used as a universal HKG across species as well as across stages of insects [8,43]. Simultaneously, its use as an internal control has been contradictory in many cases [44,45]. Various studies have already been accomplished by identifying HKGs in different species and tissues. Recently, numerous reference genes have been testified under different conditions in stinkbug (Bansal et al, 2016), western flower thrips [46], honey bee [47], fly [48], silkworm [49], moths [50] and beetles [51]. The analysis of reference genes in cotton leafhopper across different developmental stages showed variation in expression as reported in several other insects [52]. These genes have significant involvement in ubiquitous cellular and biological processes, therefore cannot be used as single normalizer in expression studies. In order to improve the efficiency and accuracy of qRT-PCR, we need to validate various internal genes for their expression. RT-qPCR studies revealed variability in the expression of reference genes across samples inferring that a universal reference gene cannot be used for all species or all experimental conditions [53]. Thus validation of reference HKGs is obligatory to gene expression studies wherein a candidate reference gene should have amplification efficiency similar to the target gene.
In the current study, we evaluated fifteen candidate HKGs of cotton leafhopper using multiple statistical models (Bestkeeper, geNorm, NormFinder, and RefFinder) across different developmental stages and starvation stress for identification of reference genes as suggested by [54]. Interestingly, results obtained from different algorithms were highly variable. As each program is using the unique algorithm, incongruities are to be anticipated. For instance, the BestKeeper takes into account the InVar (intrinsic variance), SD (standard deviation) and Pvalues, all of which contribute to the BestKeeper vs. Pearson correlation coefficient value. These compounding factors result in the obvious differences in outcomes. To integrate the results of all three algorithms, we further used RefFinder for comprehensive ranking of the genes. RefFinder abdicated the ambiguity in the results and helped in the compilation of three algorithms as shown in Table 7.
Our results suggested that the RP13 gene was the most stable gene across all the developmental stages tested. The RP13 gene encodes a ribosomal protein that is a component of the 60S subunit, which is involved in the translation initiation [55]. Various ribosomal proteins genes have been evaluated as reference genes for RT-qPCR in many insects, and these genes have been reported to show the most stable expression in Tetranychus cinnabarinus (RPS18: Selection of housekeeping genes and RNAi in the cotton leafhopper [56]), Apis mellifera (RPS18: [47]), Rhodnius prolixus (RPS18: [57]), Cimex lectularius (RPL18: [58]), and Schistocerca gregaria (RP49: [59]). In late and early stage nymphs, Tub and 28s were the most stable genes, respectively. Previous studies also demonstrated 28s rRNA to be the stable and suitable gene for internal control in various organisms ( [47,[60][61][62][63]. Tubulin has also been most commonly used as a reference gene in various expression studies [64][65][66]. These belong to the family of eukaryotic structural genes that form microtubules, fundamental components of the cytoskeleton that mediate cell division, shape, motility, and intracellular trafficking [67]. Nevertheless, tubulin is reported to be the most variable gene in plant hoppers [68]. Under starvation stress, UbiCF has come out to be the most stable gene in a comprehensive analysis. Ubiquitin is the founding member of a family of structurally conserved proteins that regulate a host of processes (protein degradation, DNA repair, signal transduction and transcription regulation by endocytosis) in insect cells [69][70][71]. Various genes coding for ubiquitin proteins such as ubiquitin conjugation factor (Bansal et al, 2016), polyubiquitin [68] and ubiquitin [64] evaluated in different insect species have shown stable expression as HKGs. Earlier reports concluded that it was difficult to identify most stable reference gene across various developmental stages in an organism [72][73]. We set our parameters under three algorithms (i.e. Bestkeeper, geNorm, Normfinder) to find the stable HKGs across developmental stages and starvation stress. However, the results obtained from individual algorithm were variable. Thus, we used RefFinder to comprehend the results of three algorithms and assign a ranking to the genes based on their stability under each and across (overall ranking) developmental stages and starvation stress. The overall ranking concludes RP13, VATPase and TUB as the most stable genes which could be helpful in the expression studies involving starvation stress and mixture of individuals from different life stages of leafhopper. Finally it may be concluded that top two-three genes could be selected and validated to elucidate the appropriate reference gene(s) under a particular set of experimental conditions for future gene expression studies in this insect. Feeding of dsRNA for causing RNAi has been demonstrated in a few hemipteran insects including Halyomorpha halys [36]; Bemisia tabaci [74]; Acyrthosiphon pisum [75]; Bactericera cockerelli [76]; Nilaparvata lugens [77]; Sogatella furcifera [78]; Laodelphax striatellus [79]. Many insects in this order have also shown good RNAi response through injection of dsRNA [80][81]. In planta expressed dsRNA too has been successful in harnessing RNAi response in few of the hemipteran insects [82][83][84][85][86][87]. Gene silencing has been successfully demonstrated for many hemipteran insects through the oral delivery of synthetic dsRNAs dissolved in sucrose via artificial membrane [84,76,88]. The slightly similar and modified procedure has been used for delivering dsRNA in cotton leafhopper in the current study. Robust RNAi efficiency has been demonstrated in cotton leafhopper through the feeding of dsSnf7. Snf7 belongs to the ESCRT (Endosomal Sorting Complex Required for Transport)-III complex, the ESCRT pathway is a key regulator of biological processes important for eukaryotic cell growth and survival [89]. Feeding of dsRNA targeting western corn rootworm, Snf7 homolog in larvae caused severe stunting after five days of exposure followed by the death of the larvae [23]. The highest mortality was observed with the feeding of dsSnf7 to nymphs of cotton leafhopper. This may be correlated to 77.12% (4.3 fold) decline in the expression of the targeted gene compared to the GFP control. Aquaporin water channels have been implicated in mediating the mass transfer of water in a various physiological processes. Hemipteran water-specific aquaporins have been reported in the gut of phloem-feeding leafhopper Cicadella viridis [90] and B. tabaci [91]. The expression of the Aquaporin 1 gene was reduced to 70% at six days post-feeding and caused 84% mortality in B. tabaci adults fed on 20 μg/ml AQP1dsRNA in 20% sucrose (Vyas et al, 2017). The reports showed no mortality in A. pisum fed on one μg/μl dsApAQP in the diet [92]. However, there was a two-fold reduction in the ApAQP expression compared to that in control. Thus, it is evident that the concentration of dsRNA might be a key factor for causing significant RNAi effect. In our experiments, dsRNA at 500 ng /μl diet caused 66.59% (2.99-fold) reduction in AbbAQP mRNA levels compared to that in dsGFP-fed control insects. However the percent corrected mortality caused by dsAbbAQP feeding was low compared with dsAbbSnf7 feeding. The mortality observed with the feeding of dsAbbIAP and dsAbb-VATPase was also low. Nevertheless, there was a significant reduction of 59.87 and 56.17% in the mRNA levels of both the genes in the dsRNA fed insects. In some of the hemipteran insects, the downregulation of VATPase has resulted in a significant mortality as compared to the control [76,74,86,93]. RNAi mediated screening of 290 genes of western corn rootworm (D. virgifera virgifera) revealed that only 125 genes were able to cause larval mortality [94]. In lepidopteran insects, out of 130 genes, only 49 genes (38%) were silenced more than 50% whereas 18 genes (14%) showed lesser silencing while 62 (48%) genes did not show silencing at all [95]. In insects, Hyalophora cecropia, Antheraea pernyi and Manduca sexta even ten ng dsRNA per mg of insect biomass was sufficient to elicit high levels of target gene silencing [96][97][98]. Similarly, silencing of target genes in coleopteran beetles could be achieved with small quantities of dsRNA [99]. Contrastingly, very high concentration (100 μg) of dsRNA per mg of insect biomass was required to silence target gene of Antheraea mylitta [100].The length of dsRNA is also one of the critical parameters that determine dsRNA uptake and RNAi efficiency [101][102]. To conclude, the RNAi efficiency in cotton leafhopper with different genes showed a differential response. Thus exhaustive studies in future may be helpful to better understand the RNAi and its core machinery in this insect. Overall, the current study has identified stable reference genes across various developmental stages and starvation stress in A. biguttula biguttula. Considered in concert, using different software algorithms and comprehensive analysis results, we suggest that RP13 and VATPase are the most suitable HKGs for gene expression studies across developmental stages. RefFinder concluded RP13, VATPase, 18s in adult, 28s, RP13, TUB in early stage nymph, TUB, TATA, Actin in late-stage nymph and UbiCF, B-Tub, G3PD in starvation stress as the top three stables genes. However, considering the individual developmental stage RP13, Tub and 28S are best suited for adult, late, and early nymphal stages, respectively. In addition, UBiCF is most stable HKG under starvation stress. This study also reports for the first time successful gene silencing through RNAi in A. biguttula biguttula.
Supporting information S1