Transcriptional responses underlying the hormetic and detrimental effects of the plant secondary metabolite gossypol on the generalist herbivore Helicoverpa armigera
© Celorio-Mancera et al; licensee BioMed Central Ltd. 2011
Received: 22 July 2011
Accepted: 23 November 2011
Published: 23 November 2011
Hormesis is a biphasic biological response characterized by the stimulatory effect at relatively low amounts of chemical compounds which are otherwise detrimental at higher concentrations. A hormetic response in larval growth rates has been observed in cotton-feeding insects in response to increasing concentrations of gossypol, a toxic metabolite found in the pigment glands of some plants in the family Malvaceae. We investigated the developmental effect of gossypol in the cotton bollworm, Helicoverpa armigera, an important heliothine pest species, by exposing larvae to different doses of this metabolite in their diet. In addition, we sought to determine the underlying transcriptional responses to different gossypol doses.
Larval weight gain, pupal weight and larval development time were measured in feeding experiments and a hormetic response was seen for the first two characters. On the basis of net larval weight gain responses to gossypol, three concentrations (0%, 0.016% and 0.16%) were selected for transcript profiling in the gut and the rest of the body in a two-color double reference design microarray experiment. Hormesis could be observed at the transcript level, since at the low gossypol dose, genes involved in energy acquisition such as β-fructofuranosidases were up-regulated in the gut, and genes involved in cell adhesion were down-regulated in the body. Genes with products predicted to be integral to the membrane or associated with the proteasome core complex were significantly affected by the detrimental dose treatment in the body. Oxidoreductase activity-related genes were observed to be significantly altered in both tissues at the highest gossypol dose.
This study represents the first transcriptional profiling approach investigating the effects of different concentrations of gossypol in a lepidopteran species. H. armigera's transcriptional response to gossypol feeding is tissue- and dose-dependent and involves diverse detoxifying mechanisms not only to alleviate direct effects of gossypol but also indirect damage such as pH disturbance and oxygen radical formation. Genes discovered through this transcriptional approach may be additional candidates for understanding gossypol detoxification and coping with gossypol-induced stress. In a generalist herbivore that has evolved transcriptionally-regulated responses to a variety of different plant compounds, hormesis may be due to a lower induction threshold of growth-promoting, stress-coping responses and a higher induction threshold of detoxification pathways that are costly and cause collateral damage to the cell.
Gossypol is a polyphenolic binaphthyl dialdehyde found in the pigment glands of plant species belonging to Malvaceae, most notably cotton, Gossypium spp. . The toxicity of gossypol in cottonseed meal has spurred a multitude of studies of its effects on diverse biological systems. The presence of both phenolic and carbonyl groups in the gossypol structure allows the interaction with other molecules through a wide range of weak to strong chemical bonds. This compound can covalently bind to amino acids, particularly lysine, through Schiff's base condensation reactions, and its dimeric structure facilitates cross-linking of proteins. Gossypol can also chelate iron and other metal ions and be both a pro-oxidant and antioxidant [2, 3]. Moreover, the aromatic rings in gossypol render it a fairly hydrophobic compound able to penetrate biological membranes more readily . All this intricate chemical reactivity implies the existence of multiple targets, which is reflected by the different enzyme classes inhibited by gossypol, especially those involved in energy production processes of the cell, and in the variety of biological effects attributed to this compound due to its antitumor, spermicidal, antimalarial, antiparasitic, antiamoebic and antiviral activities [3, 5].
There is only limited knowledge regarding the mode of action of gossypol regardless of its therapeutic or toxic effect on different biological systems. It is clear, however, that the mechanism mostly relies on the disturbance of proteins and membranes. In addition, gossypol occurs in two optically active forms, of which the (-)-enantiomer has been found to be more reactive towards biological systems, particularly non-ruminant animals. In fact, the molecular mechanism behind (-)-gossypol antineoplastic activity appears to rely on its interaction with antiapoptotic proteins in the outer mitochondrial membrane . However, the inhibitory effect of gossypol on insects and fungi is not enantiospecific . Therefore, the ability to extrapolate across systems is limited and the response within a given system is likely to be very complex and dose-dependent.
Both quantity and quality of gossypol in plant tissues vary according to several factors such as cultivar, phenology or plant organ [1, 8–10]. The protective nature of this sesquiterpene dimer to cotton against different insect herbivores became evident after the development of glandless cotton lines . These lines were produced by conventional breeding to reduce the amount of gossypol in cottonseed meal which is fed to livestock. Some members of the "major-pest lineage" within the subfamily Heliothinae (Lepidoptera: Noctuidae) , prefer glandless cotton varieties over glanded ones . Moreover, high concentrations of gossypol inhibit growth and development of tobacco budworm, Heliothis virescens. However, the dose response is not monotonic and in fact gossypol acts as a hormetic agent, since the maximum weight gain was observed when larvae fed on a low dose gossypol-containing diet (0.0125%). Hormesis, first described in the late 1800's, is a widespread biological response towards an environmental stressor which is stimulatory in low amounts (e.g. increased fertility), and detrimental in high amounts (e.g. toxicity) .
Although the implications of the hormetic effect in the area of human toxicology remains controversial , there is substantial research regarding hormesis due to environmental stressors and their impact on aging and longevity in invertebrate species , and the efforts to understand the biology underlying the phenomenon continue . Recent analyses for two animal systems (mammal and annelid) inquired whether gene expression in a global scale differs under the "low-dose effect" of toxic chemicals [18, 19]. Dose-dependent transcriptional responses to gossypol have been observed for a handful of genes including the expression of mitochondrial-related genes in rat liver cells  and P450 monooxygenases in insects [21, 22]. However, there is no unbiased transcriptional profiling approach documented investigating the effects of different doses of gossypol.
The cotton bollworm (CBW), H. armigera, is a generalist herbivore and one of the most injurious insect pests in the world, damaging a large number of plant species  with cotton as one of its most suitable hosts . We investigated whether a hormetic effect of gossypol can be observed in H. armigera by measuring developmental parameters in response to different concentrations in the diet. In addition, we investigated the transcriptional response of H. armigera gut as well as the rest of body towards different dosages of gossypol using microarray expression profile analysis. The results herein provide a general view of the underlying transcriptional response to gossypol with implications for the detoxification of this plant secondary metabolite. They also indicate the role hormesis may play in the adaptation of generalist herbivores to a variety of hostplants presenting different sets of chemical challenges.
H. armigera larvae were collected from Toowoomba, Queensland, Australia, in 2003 and reared on artificial diet under laboratory conditions (26°C, 55% RH, 16:8 hr = L:D) in Jena, Germany. The colony was maintained for 25 generations prior to the start of this study, exclusively by about 50 single-pair crosses per generation, avoiding brother-sister mating to minimize inbreeding. The artificial diet for larval rearing was purchased from BioServ (F9772, Frenchtown, NJ, USA). Gossypol from cotton seeds (G8761, racemic mixture, Sigma-Aldrich, MO, USA) was incorporated into the artificial diet using a carrier material of non-nutritive cellulose, "Alphacel", as reported elsewhere  in various concentrations (CT = 0.0, T1 = 0.0004, T2 = 0.0016, T3 = 0.004, T4 = 0.008, T5 = 0.016, T6 = 0.04, or T7 = 0.16% (w/v)). Cotton seeds contain a racemic mixture of gossypol with generally a higher percentage of the (+)-enantiomer although this may differ depending on the cotton species and cultivar .
Gossypol treatment and larval development
A total of nine hundred and sixty larvae from a pool of newly molted fifth instar larvae were evenly divided into eight groups according to their initial fresh weights (I). Each group of larvae was randomly assigned to each gossypol treatment. Control group of larvae were exposed to the diet supplemented only with the Alphacel carrier (CT). Each group of larvae was composed of four biological replicates of thirty larvae. After those larvae were subjected to each treatment for three days, larvae were individually weighed as a measure of their final fresh weight (F), in order to calculate the net weight gain (= F-I).
Mortality was observed in CT, T4 and T7 treatments (4, 3 and 5 larvae were dead, respectively). Only one individual died in each T1 and T6, while no mortality was observed in the remaining treatments. Ten randomly chosen larvae from each replicate per treatment were separated for microarray analysis (see below), the other eighty larvae (20 larvae × 4 replicates) were returned to the same diet treatments and further observed for other developmental parameters: larval development time to pupation (day) and pupal weight in a day after pupation (mg). The developmental parameter data were analyzed by ANOVA and the statistical differences among treatment means were further tested by Post-hoc analysis.
Expressed Sequence Tag (EST) Project
RNA was extracted from several larval tissues (e.g. midgut, fat body, integument), all larval instars and developmental stages (larvae and adults) with TRIzol Reagent (Invitrogen) according to the manufacturer's protocol. An additional DNAse (Turbo DNAse, Ambion) treatment was included to eliminate any contaminating DNA. The DNAse enzyme was removed and the RNA was further purified by using the RNeasy MinElute Cleanup Kit (Qiagen) following the manufacturer's protocol. RNA integrity and quantity was verified on an Agilent 2100 Bioanalyzer using the RNA Nano chips (Agilent Technologies). RNA quantity was determined on a Nanodrop ND-1000 spectrophotometer. Normalized full length-enriched cDNA libraries were generated using a combination of the SMART cDNA library construction kit and the Trimmer Direct cDNA normalization kit (Evrogen) following the manufacturer's protocol with several modifications . The normalization process facilitates the identification of low abundance transcripts.
Single-pass sequencing of the 5' termini of cDNA library plasmid clones was carried out on an ABI 3730 xl automatic DNA sequencer (PE Applied Biosystems). Vector clipping, quality trimming and sequence assembly was done with the Lasergene software package (DNAStar Inc.). In total, 8 different cDNA libraries were generated from the tissues and developmental stages as described above and ~60,000 clones were sequenced. Additional sequencing was performed with a mixed cDNA pool on a Roche 454 FLX instrument, obtaining 274,607 high quality reads after trimming and quality filtering steps. The H. armigera ESTs generated and all publicly available Genbank sequences for this species were jointly assembled using Seqman NGen (Lasergene) and clustered into 27,381 contigs (putative gene objects) subsequently used for microarray oligo probe design.
Microarray Design, labeling, hybridization and data acquisition
In order to optimize our H. armigera microarray design and maximize the output of subsequent gene expression profiling experiments, a Pre Selection Strategy (PSS, Imagenes) approach was used to select well performing probes based on initial test hybridizations. For the preliminary large array design each gene was tiled by a maximum number of probes. A total of 231,399 oligos for the 27,381 contigs were designed and a 244 K Agilent microarray was hybridized with labeled complex total RNA mixture and genomic DNA. The best performing probes for each gene were selected, for the expressed genes based on the RNA hybridization and for the non-expressed genes based on the DNA hybridization. A final condensed Agilent 4 × 44 K array design based on the eArray platform (Agilent Technologies; https://earray.chem.agilent.com/earray/) contains the few best performing probes of each gene (1-2 for each Gene Object) with a final number of 43803 non-control probe set and 1417 Agilent Technologies built in controls (structural and spike in).
A subset of forty ice-cold anesthetized larvae per treatment (having fed on 0%, 0.016%, or 0.16% gossypol diet for three days) was dissected longitudinally under phosphate-buffered saline and gut tissues (G) were separated from the rest of the body (RB). Four biological replicates were prepared by pooling ten individuals in each replicate of either G or RB tissues, snap-frozen in liquid nitrogen, and stored at -80°C until RNA isolation. Total RNA was purified, quality tested and quantified as mentioned above. Agilent Technologies spike-in RNA was added to 500 ng of total RNA and labelled using the Low RNA Input Linear Amplification kit (Agilent Technologies). Treated RNA and control samples were labelled with cyanine 5-CTP and 3-CTP dyes according to manufacturer instructions following a double reference dye-swap design. Labelled amplified cRNA samples were purified using RNeasy MinElute Cleanup kit (Qiagen) and analyzed on a Nanodrop spectrophotometer using the microarray function. Amplified cRNA samples were used for microarray hybridization only if the yield is > 825 ng and the specific activity is > 8.0 pmol Cy3 or Cy5 per ug cRNA. 825 ng each of cyanine 3 and cyanine5 labeled cRNA were used for each array. Hybridization was carried out at 65°C for 17 hours. Slides were washed in GE Wash Buffer 1 (Agilent Technologies) for 1 min at room temperature and a further minute in GE Wash Buffer 2 pre-warmed overnight to 37°C. Slides were treated in stabilization and drying solution (Agilent Technologies), scanned with the Agilent Microarray Scanner, and data was extracted from the TIFF images with Agilent Feature Extraction software version 9.1. The initial technical validation included visual inspection of images to identify gross abnormalities or background. Prior to normalization the sensitivity of the array and relationship between RNA concentration and fluorescent signal was assessed by calculating the signal intensity generated by reporters complementary to 10 'alien' synthetic RNA spikes introduced at known concentrations (from 1 pmole to 30 nmole prior to labeling). The microarray data reported in this paper have been deposited in the Gene Expression Omnibus (GEO) database, http://www.ncbi.nlm.nih.gov/geo [GEO:GPL14736].
Expression profiling of H. armigera G and RB samples subjected to different gossypol-containing diets was generated by normalizing fluorescence signals to the median intensity and log base 2-transformation of the normalized data. In order to determine the relationship between the samples per tissue, the clustering application (Euclidean distance, average linkage) was applied to normalized to median, log-transformed, statistically significant data after ANOVA (unequal variance, no threshold, Benjamini and Hochberg false discovery rate (B&H FDR) multiple test correction, adjusted P cut off < 0.001) using the Geospiza GeneSifter® genetic analysis software. Data was also filtered by volcano plots comparing each gossypol dosage to its control per tissue treatment by means of an unpaired t-test, unequal variance using Agilent GeneSpring GX11.5.1 software. All 43863 probes passed the data quality filtering based on intensity measurements. Only probes with corrected P values less than 0.001 after B&H FDR were considered statistically differentially expressed. Gene Ontology (GO)  annotations (for the target genes represented in the microarray) obtained through Blast2go (E-value cutoff: 0.001)  software were used to find significant GO terms (P < 0.001) represented in the statistically filtered data using GeneSpring. In order to perform a gene level analysis using GeneSpring program, the probes of the microarray along with the target EST sequences were assembled using Sequencher 4.7, contigs were blasted to public and in-house EST databases to assign correct gene codes to each probe. Nomenclature for important H. armigera detoxification genes such as esterases, cytochrome P450 s (P450 s), UDP-glycosyltransferases (UGTs) and glutathione transferases (GSTs) was included using public databases or by assigning it based on homology. Nomenclature for P450 s and UGTs was approved by the relevant Nomenclature Committees. Based on a list of homologous genes, we also inspected Drosophila melanogaster pathway enrichment in the data set for each t-test comparison using GeneSifter® applying z-score statistics to determine whether a pathway occurs more or less frequently than expected as previously described by others . The pathway source originated from the Kyoto Encyclopaedia of Genes and Genomes (KEGG) (Kanehisa Laboratories). Homology was established by obtaining best BLAST  hits (about 24% of the probes had a hit at e < 5 × 10-2) for H. armigera ESTs to D. melanogaster reviewed reference sequences (RefSeq) (See additional file 1: Table S1 displaying best RefSeq homolog hits to H. armigera). Homology databases containing the protein and nucleotide refseq accession were prepared from NCBI Entrez. In addition, the normalized log-ratios for each gene in all eight biological replicates per treatment comparison (relative to control), were used to apply the rank products (RP) method which has been proven to be a powerful method for identifying biologically relevant gene expression changes .
Quantitative real-time PCR
Single stranded DNA from 500 ng of total RNA was obtained and amplified using Verso SYBR® Green 2-Step QRT-PCR kit (Thermo Fisher Scientific) following manufacturer's instructions. Real-time PCR oligonucleotide primers were designed on the basis of sequences obtained for H. armigera CYP6AE14 and CYP6AE11 and two additional genes used as normalizers (endogenous control genes) i. e. 18 S ribosomal RNA and elongation initiation factor 4. Data was analyzed using the qBase 1.3.1 software (Biogazelle) and graphed as fold expression relative to the lowest expression across treatments.
Effect of gossypol on larval development
Effect of gossypol on gene expression
Since larval net weight gain was affected significantly only at the low T5 and high T7 gossypol doses, the transcriptional profiles of gut tissue (G) and the rest of the body (RB) from larvae collected from these treatments were examined by microarray analysis in a two-color double reference design using the 0% gossypol tissue samples as reference (See additional file 3: Figure S2-A depicting the microarray two-color design). Sample clustering using the expression profiles of the probes which passed the statistical filtering in GeneSifter yielded dendrograms grouping the control (CT) and T5 treatments closer to each other within each tissue type (See additional file 3: Figure S2-B depicting the hierarchical clustering of transcriptional responses across gossypol dose-tissue conditions).
Gene Ontology and KEGG Pathways Analyses
GO categories enriched in microarray comparisons through gossypol-tissue treatments
mitochondrial inner membrane
organelle inner membrane
generation of precursor metabolites and energy
cellular ketone metabolic process
monovalent inorganic cation transmembrane transporter activity
carboxylic acid metabolic process
oxoacid metabolic process
organic acid metabolic process
lipid biosynthetic process
monovalent inorganic cation transport
iron ion binding
T5 - Rest of body/DOWN
homophilic cell adhesion
T7 - Rest of body/DOWN
threonine-type peptidase activity
proteolysis involved in cellular protein catabolic process
negative regulation of ubiquitin-protein ligase activity
positive regulation of ubiquitin-protein ligase activity
regulation of ubiquitin-protein ligase activity involved in mitotic cell cycle
regulation of ubiquitin-protein ligase activity
positive regulation of ubiquitin-protein ligase activity involved in mitotic cell cycle
negative regulation of ubiquitin-protein ligase activity involved in mitotic cell cycle
negative regulation of ligase activity
positive regulation of ligase activity
regulation of ligase activity
cellular macromolecule catabolic process
cellular protein catabolic process
modification-dependent macromolecule catabolic process
proteasomal ubiquitin-dependent protein catabolic process
positive regulation of protein ubiquitination
negative regulation of protein ubiquitination
regulation of protein ubiquitination
anaphase-promoting complex-dependent proteasomal ubiquitin-dependent protein catabolic process
protein catabolic process
proteasome accessory complex
modification-dependent protein catabolic process
proteasome core complex, alpha-subunit complex
proteasomal protein catabolic process
macromolecule catabolic process
ubiquitin-dependent protein catabolic process
proteasome core complex
proteasome regulatory particle
threonine-type endopeptidase activity
T7 - Rest of body/UP
structural constituent of cuticle
FAD or FADH2 binding
alcohol metabolic process
oxidoreductase activity, acting on CH-OH group of donors
At the higher gossypol concentration T7, the oxidoreductase molecular function was significantly enriched in both gut and RB. Among these genes several P450 enzymes are significantly upregulated, notably CYP6AE14 and CYP6AE11, while a few are downregulated (Figure 2). Several UGTs and carboxylesterases, and a fewer number of GSTs are upregulated, predominantly in the midgut (Figure 2). Many genes were down-regulated in G-T7; however, no particular GO enrichment was observed even if a less stringent p-value cut off (0.01) was applied to the analysis. Many ESTs down-regulated in the RB in the gossypol T7 treatment were found to be associated with the proteasome complex. (See additional file 5: Table S3 containing a list of differentially expressed genes across conditions with no GO enrichment.)
Rank Products Analysis
Quantitative real-time PCR
Additionally, we inspected the gene expression of CYP6E14 and CYP6AE11 (with 91% nucleotide identity to CYP6E14) across all gossypol concentrations by qRT-PCR (See additional file 17 Figure S3) CYP6AE14 and CYP6AE11 are highly up-regulated at the highest gossypol concentration. CYP6AE11 expression declines slightly over 0.004 to 0.016% gossypol concentrations compared to the control 0%, and then increases at the highest concentration, roughly mirroring the hormetic growth response.
Effect of gossypol on larval growth
The biphasic response to gossypol concentration is similar to observations for another heliothine species, Heliothis virescens where the highest larval weight was obtained when the insects were exposed to a 0.0125% gossypol-containing diet. Before our study, the effect of gossypol doses lower than 0.0125% on larval development had not been examined. The maximum stimulatory concentrations of gossypol on weight gain in heliothine larvae may lie between 0.008% and 0.04%. Data in Figure 1 are also consistent with the results of Stipanovic et al.  where the developmental effect of 0.16% gossypol and higher towards Helicoverpa zea larvae resulted in a significant extension of the time to reach pupation and decrease of pupal weight regardless of gossypol enantiomer or racemate used in the diet. The results herein indicate that the inhibitory effects of gossypol on H. armigera can be observed even if the insect is exposed later in its development, since we did not rear the insects on the gossypol-containing diet from the neonate stage as done by other researchers [13, 25, 32]. In addition, a hormetic effect on larval weight gain can also be observed on newly molted fifth instar larvae even after three days of exposure to gossypol, although in this regard we have no data on the relative effect of the racemic mixture versus the separate enantiomers. Thus hormesis seems to be a general feature of the response of generalist heliothines to gossypol. The many host-plant specialists within this insect group, which do not normally encounter gossypol in their diet, have not been tested for hormetic effects.
Coping with gossypol
It has been suggested that bacterial infection stimulates epithelial renewal via reactive oxygen species (ROS) production . A similar stimulatory response may be triggered by low gossypol concentrations in the CBW gut, since a consequence of gossypol toxicity may involve the generation of superoxide free radicals damaging the epithelium . However, if present, such a proliferative response does not seem to be mediated by significant transcriptional changes in the midgut at the low gossypol concentration.
Gossypol, as a lipophilic compound, interacts with cellular membranes and forms bonds with the amine groups of proteins via Schiff's base formation; cross-linking of membrane proteins has been suggested to block cell-to-cell communication . Previously, it has been proposed that adaptive changes in lipid content can occur upon gossypol exposure [34, 35]. We speculate that this effect on genes involved in cell adhesion may represent a response to the cytotoxicity of gossypol at this concentration. A down-regulation of cell adhesion genes has also recently been observed as a consequence of gossypol exposure in mammalian cells . A 12-cadherin domain protein is expressed in the midgut of lepidopteran larvae where it is a binding target for Cry1A toxins from Bacillus thuringiensis (Bt) . It has recently been shown that the same cadherin is also expressed in the larval testis of the cotton specialist pink bollworm (Pectinophora gossypiella) suggesting a role in sperm production . Mutations in this protein responsible for Bt toxin resistance reduce fertility somewhat and in turn Bt-resistant pink bollworm is more susceptible to gossypol . We observed down-regulation of the homologous 12-domain cadherin protein, not in the midgut where it abundantly expressed but in the rest of the body which includes the testis. In vertebrates, some cadherins are known to interact with compounds leading to signal transduction events such as the carcinoembryonic antigen-related cell adhesion molecules (CEACAMs); however insect homologs of the latter are unknown.
The mode of action of gossypol has been studied mostly in relation to its toxic or therapeutic effects on mammalian cells. These studies show that certain low doses of gossypol can have antitumor, antiviral and antiparasitic activities mostly due to the fact that gossypol inhibits key proteins belonging to different classes of enzymes such as oxidoreductases, hydrolases and transferases . Gossypol is not considered a bactericidal compound . Therefore, we do not consider that the hormetic effect observed at T5 is a consequence of gossypol acting as a toxin to bacteria possibly present in the artificial diet.
A stimulatory effect on energy metabolism due to gossypol has been observed in cultured mouse cells. Relatively low doses of gossypol were added to the cell medium along with glucose resulting in an increase of lactate production. Additionally, the inhibition of oxygen consumption produced by 5 mM glucose was reversed when gossypol was added to the medium . The phenomenon was explained by the ability of gossypol to uncouple oxidative phosphorylation, leading to a decrease of mitochondrial production of adenosine triphosphate (ATP). Thus, glycolysis and the production of lactate through the Embden-Meyerhof pathway were stimulated as an attempt to maintain the required ATP levels in the cell. The authors concluded that the response towards gossypol may depend on the cell's glycolytic and mitochondrial oxidative phosphorylation capacity and on its ability to maintain acid-base homeostasis . A differential susceptibility and hormetic effect in oxygen production has also been observed between somatic and germ rat cells exposed to different concentrations of gossypol . Therefore, it is possible that the gut epithelium can be maintained at low gossypol doses and even stimulated by the uncoupling of mitochondria. This may explain the fact that we do not see the biological process of glycolysis enriched by differentially regulated genes since the response to gossypol may depend on subtle adjustments in metabolism which take place as acclimatization at the protein level. However, glycolysis-related genes were affected in the RB at both gossypol concentrations examined.
It has previously been described that β-fructofuranosidases of Bombyx mori are resistant to inhibition by mulberry sugar mimic alkaloids which inhibit α-glucosidases, and therefore represent an adaptation to the mulberry host in this specialist moth . If α-glucosidases are also inhibited by gossypol, the up-regulation of β-fructofuranosidase genes by the high dose of gossypol in H. armigera might be explained by a similar compensatory mechanism. Interestingly, we found that in H. armigera, the regulation of at least one β-fructofuranosidase gene (GH32FruA-1), in response to gossypol, is tissue and dose dependent. The same seems to be true for glucose dehydrogenase genes, which are up-regulated at gossypol T7 only in the RB.
Cross-linking amino acids and proteins, phenolic compounds cause oxidative damage to the midgut cells of insects . Thus, the up-regulation of peroxisome-related genes may be part of an antioxidant response to the ROS generated by gossypol activity at such high concentration. Consistently, phenoloxidase inhibitor and oxidase peroxidase genes are respectively down-and up- regulated in both tissues, indicating that homeostasis has been compromised possibly due to major tissue wounding caused by oxidative stress.
The proteasome apparatus, composed of a proteolytic core and two regulatory particles, degrades a variety of cellular proteins involved in many essential functions, such as signal transduction pathways, stress signaling, inflammatory responses, and apoptosis. If increased proteolysis of cellular targets inactivated by binding to gossypol were required, we would expect an up-regulation of this system. The fact that the proteasome is instead down-regulated by gossypol suggests that the response of CBW to this compound may be similar to the response by an annelid towards fluoranthene, a polycyclic aromatic hydrocarbon, known to disrupt biological membranes due to its lipophilicity , characteristics shared with the gossypol molecule. This down-regulation of the ubiquitin-26S proteasome pathway is interpreted as a mechanism to reduce the proteolytic turnover of the aryl hydrocarbon receptor (Ahr) which is considered a mediator in the expression of genes involved in detoxification. In mammals, Ahr regulates cellular responses to certain polycyclic aromatic hydrocarbon toxins similar to fluoranthene and gossypol . Several pathways also mediating gene expression (i.e. mTOR, JAK-STAT, phototransduction and Notch) were affected, representing potential signaling pathways involved in the response towards gossypol which deserve further attention.
Since gossypol is a defensive chemical encountered in some but not all of the hostplants of this generalist herbivore, some transcriptional responses are expected to be directed towards its detoxification, and these can be compared with previous metabolic studies. The metabolic fate of gossypol in Heliothis virescens larvae has been examined by means of 14C-labelling the compound and adding it to an artificial diet. 25% was found to be metabolized by conjugation with six sugar molecules per mole of gossypol and excreted in the frass, whereas glutathione conjugates were not detected . This is consistent with the up-regulation of 10 UDP-glycosyltransferases but only one glutathione transferase in the midgut. If endogenous α-glucosidases were capable of hydrolyzing these glucose conjugates, their down-regulation would be favored in favor of up-regulation of β-fructofuranosidases which can still digest carbohydrates. About 33% of the labelled gossypol was found in the larval tissues, mostly in the fat body and the rest of the compound excreted as free gossypol or bound to components of the frass. About 10% of labelled gossypol was recovered as carbon dioxide but the mechanism has not been elucidated in this species. A similar study done with rats revealed that decarbonylation of gossypol is an important detoxification pathway; labeled 14CO2 appeared in expired air 1 hour after feeding of 14C-gossypol . However, in swine decarbonylation may not be the main detoxification pathway since the main products found were glucuronides, sulphates and unconjugated metabolites . These contrasts between species made the authors consider whether the degree of decarbonylation of gossypol indicates the susceptibility towards the allelochemical (i.e. rats, more tolerant to this phenolic compound, retain less of it in their tissues and detoxify it mostly by decarbonylation) .
A zinc-iron transporter was upregulated in the gut at T7 (Figure 5), which is consistent with gossypol acting as a sink for iron in the midgut lumen. Ferric ions have been shown to precipitate gossypol, while ferrous ions produce a soluble chelate which is precipitated by calcium . There is evidence that at least some of the gossypol decarbonylation observed in the rat digestive tract occurs by an auto-oxidation process catalyzed by ferrous ions through a free-radical chain mechanism similar to that proposed for benzaldehyde decarbonylation .
The oxidoreductase molecular function in the G and RB is one of the most significantly enriched at T7 and within this category P450 enzymes are prominent. In order to metabolize gossypol to gossic acid several oxidation steps are required , for which the P450 s are candidates. Consistent with previous results by others , CYP6AE14 is up-regulated in the G and the RB at T7. However, CYP6AE14 is slightly down-regulated at T5 in both tissues. Moreover, 13 additional P450 s are upregulated and two downregulated in the gut at the highest gossypol concentration. When CYP6AE14 was silenced in H. armigera by feeding dsRNA or transgenic cotton transformed with a construct expressing dsRNA, larval growth was greatly reduced in the presence of gossypol [21, 22], however there are no heterologous expression studies providing information about CYP6AE14 substrate specificity and its direct role in gossypol metabolism. Racemic gossypol promotes the formation of superoxide free-radicals when incubated with rat liver microsomes . There is evidence indicating that the damage caused by these superoxide free radicals is due to the interaction of gossypol with the iron of the P450 enzyme . It might be worth testing whether CYP6AE14 is more susceptible to this effect which may provide an alternative explanation for its up-regulation.
The strong upregulation of several carboxylesterases in the midgut poses somewhat of an enigma, since neither gossypol nor any of its known metabolites are suitable substrates for esterases. In some cases, overexpression of esterases confers resistance to organophosphorus insecticides, due to sequestration rather than metabolism (reviewed in ). If some of these esterases similarly trap gossypol, keeping it from more sensitive cellular targets, this could be an effective but expensive mechanism of tolerance. Alternatively, gossypol may interact with a regulatory molecule such as Ahr which controls a suite of different detoxicative genes, only some of which are directly involved in the detoxification of a given xenobiotic.
Responses to low concentrations
Low concentrations of gossypol could stimulate larval growth by increasing food consumption rate, increasing food conversion efficiency, or reducing the metabolic cost of other physiological activities. Cotton specialists that can efficiently deal with gossypol could be selected to use it as a feeding stimulus, as glucosinolates are by crucifer specialists such as Plutella xylostella. If these responses were caused solely by gossypol's interacting with proteins present in the cell without causing any changes in gene expression, they would go undetected by our approach. Moreover, such responses are unlikely to have been shaped by natural selection on the herbivore specifically to benefit from low concentrations of gossypol; because such low concentrations are rarely encountered in nature.
In evaluating whether transcriptional responses can shed any light on the mechanism of hormesis of gossypol, it is worth remembering that this compound is made by the cotton plant and stored in special glands primarily to deter herbivory, and that H. armigera frequently uses cotton as a host but also consumes many other plants that do not make gossypol. The first plant consumed by a larva is dictated by its mother's oviposition choice, and there is usually limited scope for larval movement to another plant after that. We would therefore expect that natural selection would favor inducible, adaptive responses by the insect to gossypol when it is first consumed, and that many of these responses would be harmful, or at least not beneficial, in the absence of gossypol; otherwise they would be constitutively expressed. This history of coevolutionary interaction is an important difference to most studies of gossypol's effects on mammalian systems; or indeed most studies exhibiting a hormetic effect of any kind, which usually deal with pollutants or drugs that do not have an evolutionary history with the target system.
Beneficial effects of a toxin at low doses are usually attributed to completely different mechanisms than those by which the toxin exerts its ill effects at high doses. These beneficial effects are either completely unpredictable, or fall into a general pattern of activating stress-response mechanisms that are generally beneficial to the organism. However, if they were always beneficial, and incurred no additional costs, natural selection would favor their constitutive expression. Therefore one of the paradoxes of hormesis that is often overlooked is why the organism should require slight stresses to increase its fitness. Rather than having a general explanation, whether or not hormesis is observed would seem to be highly system-dependent and the precise nature of such benefits unpredictable.
The evolution of hormesis
The history of cotton-feeding insects' interaction with cotton, however, does potentially offer an evolutionary explanation for hormetic effects of a specific plant defensive compound on a generalist herbivore. We propose that different types of transcriptional responses to dietary gossypol should have different thresholds, which are shaped by natural selection. Detoxicative mechanisms that are specific to gossypol are likely to be ernegetically costly, and even cause cellular damage such as the generation of reactive oxygen species by unproductive P450 reactions. Additional stress-coping mechanisms that respond to damage directly caused by gossypol and by its detoxification are also expected to be costly. Since gossypol is either completely absent from the hostplant or present in high concentrations; such responses to high gossypol concentrations that are specific to that compound should be selected for. Conversely, highly sensitive induction of detoxicative responses at low gossypol concentrations would not be selected for, and in fact selected against if they trigger inappropriate and costly responses. However, mildly beneficial responses to gossypol that impose modest energetic costs would not be subject to similar selection for a higher threshold. They should even be selected to be induced at effectively low concentrations if initially high gossypol concentrations are reduced by metabolism.
The extremely wide diversity of systems in which hormesis has been observed makes it unlikely that a single general hypothesis could satisfactorily account for them all. However, in H. armigera hormesis emerges as one of two types of specific transcriptional adaptations to gossypol, not merely as a general response to low levels of stress. Gossypol itself exerts no ill effects at unrealistically low concentrations which are still sufficient to induce low-cost growth-promoting adaptations to gossypol. Only at higher concentrations does the damage caused by gossypol and mechanisms that cope with it overcome these benefits. We would expect the same reasoning to apply to any generalist herbivore that has developed adaptations to several different plant toxins, only one of which may be encountered in an individual's lifetime. Hormesis in such systems can thus be seen as one end of a spectrum of plant-insect coevolutionary interactions.
We thank Dr. David Nelson (University of Tennessee) for help with P450 nomenclature, Dr. Peter Mackenzie (Flinders University) for advice on UGT nomenclature, Steffi Gebauer-Jung for providing her excellent skills in Bioinformatics and Regina Seibt for the maintenance of the H. armigera Toowoomba colony in the laboratory. This work was supported by the Max-Planck-Gesellschaft.
- Lukefahr MJ, Fryxell PA: Content of gossypol in plants belonging to Genera related to cotton. Econ Bot. 1967, 21 (2): 128-131. 10.1007/BF02897860.View ArticleGoogle Scholar
- Depeyster A, Wang YY: Genetic toxicity studies of gossypol. Mutation Res. 1993, 297 (3): 293-312.View ArticleGoogle Scholar
- Abou-Donia MB: Physiological effects and metabolism of gossypol. Residue Reviews. 1976, 61: 124-161.Google Scholar
- Laughton MJ, Halliwell B, Evans PJ, Hoult JRS: Antioxidant and pro-oxidant actions of the plant phenolics quercetin, gossypol and myricetin - Effects on lipid-peroxidation, hydroxyl radical generation and bleomycin-dependent damage to DNA. Biochem Pharmacol. 1989, 38 (17): 2859-2865. 10.1016/0006-2952(89)90442-5.PubMedView ArticleGoogle Scholar
- Dodou K: Investigations on gossypol: past and present developments. Expert Opinion Investigational Drugs. 2005, 14 (11): 1419-1434. 10.1517/135437188.8.131.529.View ArticleGoogle Scholar
- Oliver CL, Miranda MB, Shangary S, Land S, Wang SM, Johnson DE: (-)-Gossypol acts directly on the mitochondria to overcome Bcl-2-mediated and Bcl-X-L-mediated apoptosis resistance. Mol Cancer Therapeutics. 2005, 4 (1): 23-31.View ArticleGoogle Scholar
- Gershenzon J, Dudareva N: The function of terpene natural products in the natural world. Nature Chem Biol. 2007, 3 (7): 408-414. 10.1038/nchembio.2007.5.View ArticleGoogle Scholar
- Cai YF, Zhang H, Zeng Y, Mo JC, Bao JK, Miao C, Bai J, Yan F, Chen F: An optimized gossypol high-performance liquid chromatography assay and its application in evaluation of different gland genotypes of cotton. J Biosci. 2004, 29 (1): 67-71. 10.1007/BF02702563.PubMedView ArticleGoogle Scholar
- Cass QB, Oliveira RV, De Pietro AC: Determination of gossypol enantiomer ratio in cotton plants by chiral higher-performance liquid chromatography. J Agric Food Chem. 2004, 52 (19): 5822-5827. 10.1021/jf049626p.PubMedView ArticleGoogle Scholar
- Yang HC, Davis DD: Variations in gossypol concentration of flower buds of cotton. Crop Science. 1976, 16 (4): 485-488. 10.2135/cropsci1976.0011183X001600040011x.View ArticleGoogle Scholar
- Cho S, Mitchell A, Mitter C, Regier J, Matthews M, Robertson R: Molecular phylogenetics of heliothine moths (Lepidoptera: Noctuidae: Heliothinae), with comments on the evolution of host range and pest status. Syst Entomol. 2008, 33 (4): 581-594. 10.1111/j.1365-3113.2008.00427.x.View ArticleGoogle Scholar
- Lusas EW, Jividen GM: A review of the first 25 years of processing and utilization research. J Amer Oil Chemists' Soc. 1987, 64 (6): 839-854. 10.1007/BF02641491.View ArticleGoogle Scholar
- Stipanovic RD, Williams HJ, Smith LA: Cotton terpenoid inhibition of Heliothis virescens development. Amer Chem Soc Symp Ser. 1986, 79-94.Google Scholar
- Kaiser J: Hormesis - Sipping from a poisoned chalice. Science. 2003, 302 (5644): 376-379. 10.1126/science.302.5644.376.PubMedView ArticleGoogle Scholar
- Douglas H: Science, hormesis and regulation. Human Exp Toxicol. 2008, 27 (8): 603-607. 10.1177/0960327108098493.View ArticleGoogle Scholar
- Le Bourg E: Hormesis, aging and longevity. Biochim Biophys Acta. 2009, 1790 (10): 1030-1039. 10.1016/j.bbagen.2009.01.004.PubMedView ArticleGoogle Scholar
- Steinberg CEW, Stürzenbaum SR, Menzel R: Genes and environment -- Striking the fine balance between sophisticated biomonitoring and true functional environmental genomics. Sci Total Environ. 2008, 400 (1-3): 142-161. 10.1016/j.scitotenv.2008.07.023.PubMedView ArticleGoogle Scholar
- Toyoshiba H, Sone H, Yamanaka T, Parham FM, Irwin RD, Boorman GA, Portier CJ: Gene interaction network analysis suggests differences between high and low doses of acetaminophen. Toxicol Appl Pharmacol. 2006, 215 (3): 306-316. 10.1016/j.taap.2006.03.009.PubMedView ArticleGoogle Scholar
- Gong P, Guan X, Inouye LS, Pirooznia M, Indest KJ, Athow RS, Deng YP, Perkins EJ: Toxicogenomic analysis provides new insights into molecular mechanisms of the sublethal toxicity of 2,4,6-trinitrotoluene in Eisenia fetida. Environ Toxicol Chem. 1995, 26: 954-959.View ArticleGoogle Scholar
- Hutchinson RW, Ing NH, Burghardt RC: Induction of c-fos, and cytochrome c oxidase subunits I and II by gossypol acetic acid in rat liver cells. Cell Biol Toxicol. 1998, 14 (6): 391-399. 10.1023/A:1007543510337.PubMedView ArticleGoogle Scholar
- Mao YB, Cai WJ, Wang JW, Hong GJ, Tao XY, Wang LJ, Huang YP, Chen XY: Silencing a cotton bollworm P450 monooxygenase gene by plant-mediated RNAi impairs larval tolerance of gossypol. Nature Biotechnol. 2007, 25 (11): 1307-1313. 10.1038/nbt1352.View ArticleGoogle Scholar
- Zhou XJ, Ma CX, Li M, Sheng CF, Liu HX, Qiu XH: CYP9A12 and CYP9A17 in the cotton bollworm, Helicoverpa armigera: sequence similarity, expression profile and xenobiotic response. Pest Manag Sci. 2010, 66 (1): 65-73. 10.1002/ps.1832.PubMedView ArticleGoogle Scholar
- Fitt GP: The ecology of Heliothis species in relation to agroecosystems. Ann Rev Entomol. 1989, 34: 17-52. 10.1146/annurev.en.34.010189.000313.View ArticleGoogle Scholar
- Liu ZD, Li DM, Gong PY, Wu KJ: Life table studies of the cotton bollworm, Helicoverpa armigera (Hubner) (Lepidoptera: Noctuidae), on different host plants. Environ Entomol. 2004, 33 (6): 1570-1576. 10.1603/0046-225X-33.6.1570.View ArticleGoogle Scholar
- Stipanovic RD, Lopez JD, Dowd MK, Puckhaber LS, Duke SE: Effect of racemic and (+)- and (-)-gossypol on the survival and development of Helicoverpa zea larvae. J Chem Ecol. 2006, 32 (5): 959-968. 10.1007/s10886-006-9052-9.PubMedView ArticleGoogle Scholar
- Fischer HM, Wheat CW, Heckel DG, Vogel H: Evolutionary origins of a novel host plant detoxification gene in butterflies. Mol Biol Evol. 2008, 25 (5): 809-820. 10.1093/molbev/msn014.PubMedView ArticleGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al: Gene Ontology: tool for the unification of biology. Nature Genetics. 2000, 25 (1): 25-29. 10.1038/75556.PubMed CentralPubMedView ArticleGoogle Scholar
- Gotz S, Garcia-Gomez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talon M, Dopazo J, Conesa A: High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008, 36: 3420-3435. 10.1093/nar/gkn176.PubMed CentralPubMedView ArticleGoogle Scholar
- Doniger SW, Salomonis N, Dahlquist KD, Vranizan K, Lawlor SC, Conklin BR: MAPPFinder: using Gene Ontology and GenMAPP to create a global gene-expression profile from microarray data. Genome Biology. 2003, 4 (1): R7-10.1186/gb-2003-4-1-r7.PubMed CentralPubMedView ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang JH, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Research. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralPubMedView ArticleGoogle Scholar
- Breitling R, Armengaud P, Amtmann A, Herzyk P: Rank products: a simple, yet powerful, new method to detect differentially regulated genes in replicated microarray experiments. Febs Letters. 2004, 573 (1-3): 83-92. 10.1016/j.febslet.2004.07.055.PubMedView ArticleGoogle Scholar
- Stipanovic RD, Lopez JD, Dowd MK, Puckhaber LS, Duke SE: Effect of Racemic, (+)- and (-)-Gossypol on Survival and Development of Heliothis virescens Larvae. Environ Entomol. 2008, 37 (5): 1081-1085. 10.1603/0046-225X(2008)37[1081:EORAGO]2.0.CO;2.PubMedView ArticleGoogle Scholar
- Buchon N, Broderick NA, Poidevin M, Pradervand S, Lemaitre B: Drosophila Intestinal Response to Bacterial Infection: Activation of Host Defense and Stem Cell Proliferation. Cell Host Microbe. 2009, 5 (2): 200-211. 10.1016/j.chom.2009.01.003.PubMedView ArticleGoogle Scholar
- Herve JC, Pluciennik F, Bastide B, Cronier L, Verrecchia F, Malassine A, Joffre M, Deleze J: Contraceptive gossypol blocks cell-to-cell communication in human and rat cells. Eur J Pharmacol. 1996, 313 (3): 243-255. 10.1016/0014-2999(96)00476-1.PubMedView ArticleGoogle Scholar
- Depeyster A, Hyslop PA, Kuhn CE, Sauerheber RD: Membrane structural functional perturbations induced by gossypol - effects on membrane order, liposome permeability, and insulin-sensitive hexose-transport. Biochem Pharmacol. 1986, 35 (19): 3293-3300. 10.1016/0006-2952(86)90426-0.View ArticleGoogle Scholar
- Moon DO, Choi YH, Moon SK, Kim WJ, Kim GY: Gossypol decreases tumor necrosis factor-alpha-induced intercellular adhesion molecule-1 expression via suppression of NF-kappa B activity. Food Chem Toxicol. 2011, 49 (4): 999-1005. 10.1016/j.fct.2011.01.006.PubMedView ArticleGoogle Scholar
- Gahan LJ, Gould F, Heckel DG: Identification of a gene associated with Bt resistance in Heliothis virescens. Science. 2001, 293 (5531): 857-860. 10.1126/science.1060949.PubMedView ArticleGoogle Scholar
- Carriere Y, Showalter AM, Fabrick JA, Sollome J, Ellers-Kirk C, Tabashnik BE: Cadherin gene expression and effects of Bt resistance on sperm transfer in pink bollworm. J Insect Physiol. 2009, 55 (11): 1058-1064. 10.1016/j.jinsphys.2009.07.013.PubMedView ArticleGoogle Scholar
- Carriere Y, Ellers-Kirk C, Biggs R, Higginson DM, Dennehy TJ, Tabashnik BE: Effects of gossypol on fitness costs associated with resistance to Bt cotton in pink bollworm. J Econ Entomol. 2004, 97 (5): 1710-1718. 10.1603/0022-0493-97.5.1710.PubMedView ArticleGoogle Scholar
- Przybylski P, Pyta K, Stefanska J, Ratajczak-Sitarz M, Katrusiak A, Huczynski A, Brzezinski B: Synthesis, crystal structures and antibacterial activity studies of aza-derivatives of phytoalexin from cotton plant - gossypol. Eur J Med Chem. 2009, 44 (11): 4393-4403. 10.1016/j.ejmech.2009.05.032.PubMedView ArticleGoogle Scholar
- Reyes J, Borriero L, Tanphaichitr N, Bellve AR, Benos DJ: Energy-metabolism of cultured tm4 cells and the action of gossypol. Biol Reprod. 1986, 34 (5): 809-819. 10.1095/biolreprod34.5.809.PubMedView ArticleGoogle Scholar
- Tvermyr SM, Froysa A, Jutte N, Hansson V: Glucose-metabolism in rat germ-cells - mechanism of action of gossypol. Annals New York Acad Sci. 1984, 438 (DEC): 543-545.View ArticleGoogle Scholar
- Daimon T, Taguchi T, Meng Y, Katsuma S, Mita K, Shimada T: Beta-fructofuranosidase genes of the silkworm, Bombyx mori - Insights into enzymatic adaptation of B. mori to toxic alkaloids in mulberry latex. J Biol Chem. 2008, 283 (22): 15271-15279. 10.1074/jbc.M709350200.PubMed CentralPubMedView ArticleGoogle Scholar
- Krishnan N, Kodrik D, Turanli F, Sehnal F: Stage-specific distribution of oxidative radicals and antioxidant enzymes in the midgut of Leptinotarsa decemlineata. J Insect Physiol. 2007, 53 (1): 67-74. 10.1016/j.jinsphys.2006.10.001.PubMedView ArticleGoogle Scholar
- Owen J, Hedley BA, Svendsen C, Wren J, Jonker MJ, Hankard PK, Lister LJ, Sturzenbaum SR, Morgan AJ, Spurgeon DJ: Transcriptome profiling of developmental and xenobiotic responses in a keystone soil animal, the oligochaete annelid Lumbricus rubellus. BMC Genomics. 2008, 9 (art. 266):
- Rojas MG, Stipanovic RD, Williams HJ, Vinson SB: Metabolism of gossypol by Heliothis virescens (F.) (Lepidoptera, Noctuidae). Environmental Entomology. 1992, 21 (3): 518-526.View ArticleGoogle Scholar
- Abou-Donia MB, Lyman CM, Dieckert JW: Metabolic fate of gossypol - metabolism of 14C-gossypol in rats. Lipids. 1970, 5 (11): 938-946. 10.1007/BF02531126.PubMedView ArticleGoogle Scholar
- Abou-Donia MB, Dieckert JW: Metabolic-fate of gossypol - metabolism of [14C]gossypol in swine. Toxicol Appl Pharmacol. 1975, 31 (1): 32-46. 10.1016/0041-008X(75)90049-6.PubMedView ArticleGoogle Scholar
- Depeyster A, Quintanilha A, Packer L, Smith MT: Oxygen radical formation induced by gossypol in rat liver microsomes and human sperm. Biochem Biophys Res Comm. 1984, 118 (2): 573-579. 10.1016/0006-291X(84)91341-X.View ArticleGoogle Scholar
- Oakeshott J, Claudianos C, Campbell P: Biochemical genetics and genomics of insect esterases. Comprehensive molecular insect science. Edited by: Gilbert L, Iatrou K, Gill S. Oxford: Elsevier, 5: 309-81.
- Thorsteinson AJ: The chemotactic responses that determine host specificity in an oligophagous insect (Plutella maculipennis (Curt.) Lepidoptera). Canadian J Zool. 1953, 31: 52-72. 10.1139/z53-006.View ArticleGoogle Scholar