The dynamics of cyanide defences in the life cycle of an aposematic butterfly: Biosynthesis versus sequestration

Heliconius butterflies are highly specialized in Passiflora plants, laying eggs and feeding as larvae only on them. Interestingly, both Heliconius butterflies and Passiflora plants contain cyanogenic glucosides (CNglcs). While feeding on specific Passiflora species, Heliconius melpomene larvae are able to sequester simple cyclopentenyl CNglcs, the most common CNglcs in this plant genus. Yet, aromatic, aliphatic, and modified CNglcs have been reported in Passiflora species and they were never tested for sequestration by heliconiine larvae. As other cyanogenic lepidopterans, H. melpomene also biosynthesize the aliphatic CNglcs linamarin and lotaustralin, and their toxicity does not rely exclusively on sequestration. Although the genes encoding the enzymes in the CNglc biosynthesis have not yet been biochemically characterized in butterflies, the cytochromes P450 CYP405A4, CYP405A5, CYP405A6 and CYP332A1 have been hypothesized to be involved in this pathway in H. melpomene. In this study, we determine how the CNglc composition and expression of the putative P450s involved in the biosynthesis of these compounds vary at different developmental stages of Heliconius butterflies. We also establish which kind of CNglcs H. melpomene larvae can sequester from Passiflora. By analysing the chemical composition of the haemolymph from larvae fed with different Passiflora diets, we show that H. melpomene is able to sequestered prunasin, an aromatic CNglcs, from P. platyloba. They are also able to sequester amygdalin, gynocardin, [C¹³/C¹⁴]linamarin and [C¹³/C¹⁴]lotaustralin painted on the plant leaves. The CNglc tetraphyllin B-sulphate from P. caerulea is not detected in the larval haemolymph, suggesting that such modified CNglcs cannot be sequestered by Heliconius. Although pupae and virgin adults contain dihydrogynocardin resulting from larval sequestration, this compound was metabolized during adulthood, and not used as nuptial gift or transferred to the offspring. Thus, we speculate that dihydrogynocardin is catabolized to recycle nitrogen and glucose, and/or to produce fitness signals during courtship. Mature adults have a higher concentration of CNglcs than any other developmental stages due to increased de novo biosynthesis of linamarin and lotaustralin. Accordingly, all CYP405As are expressed in adults, whereas larvae mostly express CYP405A4. Our results shed light on the importance of CNglcs for Heliconius biology and their coevolution with Passiflora.

Gene duplication and subsequent functional diversification of sucrose hydrolase in Papilio xuthus

Sucrose is the main product of photosynthesis in plants, providing a rich carbon and energy source for the physiological growth and development of insects. In a previous study, we identified a novel sucrose hydrolase (SUH) in the larval midgut of moths. Intriguingly, there are two copies of Suh, namely Suh1 and Suh2, in several species of butterflies. However, the biochemical characteristics of SUHs in butterflies remain unclear. In this study, we found that this duplication and subsequent diversification produced two Suh genes in Papilio xuthus. These two PxSuh genes were significantly divergent in terms of their expression pattern and enzyme properties. PxSuh messenger RNA expression was highest during the larval stage, reduced in the prepupal and pupal stages and, for PxSuh1, slightly increased again in the adult. The observed levels of PxSuh2 were overall below those of PxSuh1 amongst the development stages examined. Compared with PxSUH2, which has maintained the original gene function of maltose hydrolysis, PxSUH1 exhibits substrate specificity for sucrose with an optimum enzyme activity occurring at an alkaline pH. The data show that PxSuh1 is evolutionarily adapted for effective functioning in an alkaline digestive system. Furthermore, we find that functional diversification of Suh facilitates P. xuthus to digestive carbohydrate of host plants. Thus, our findings offer new insights into the ecological and evolutionary adaptation of digestive enzymes in butterflies.

Biochemical characterization a digestive trypsin in the midgut of large cabbage white butterfly, Pieris brassicae L. (Lepidoptera: Pieridae)

A comprehensive study on digestive trypsin was undertaken in the larval midgut of Pieris brassicae L. Results of enzymatic compartmentalization showed a significantly higher activity of crude trypsin in the anterior larval midgut rather than posterior-midgut. Using Diethylaminoethyl cellulose fast flow column chromatography a purified trypsin was obtained by specific activity of 21 U mg-1 protein, recovery of 22%, purification fold of 28-fold and molecular weight of 25 kDa. This purified enzyme showed the highest activity at pH 8 and the corresponding temperature of 40°C. However, the specific inhibitors used including 4-(2-Aminoethyl) benzenesulfonyl fluroride hydrochloride, N-p-Tosyl-L-lysine methyl ester hydrochloride and Soybean Trypsin Inhibitor significantly lowered the activity of the purified enzyme in vitro. Moreover, the activity of trypsin and likewise the nutritional indices were significantly altered in the larval midgut feeding upon the leaves treated by 1 mM concentration of each inhibitor in comparison with control. Determination of enzymatic characteristics of insect trypsins is crucial in paving the path for controlling pests by potential natural compounds via transgenic plants.

Inhibitory effects of an extract from non-host plants on physiological characteristics of two major cabbage pests

The diamondback moth (Plutella xylostella) and small white cabbage butterfly (Pieris rapae) are the two main serious pests of cruciferous crops (Brassicaceae) that have developed resistance to chemical control methods. In order to avoid such resistance and also the adverse effects of chemical pesticides on the environment, alternative methods have usually been suggested, including the use of plant enzyme inhibitors. Here, the inhibitory effects of proteinaceous inhibitors extracted from wheat, canola, sesame, bean and triticale were evaluated against the digestive α-amylases, larval growth, development and nutritional indecs of the diamondback moth and small white cabbage butterfly. Our results indicated that triticale and wheat extracts inhibited α-amylolytic activity in an alkaline pH, which is in accordance with the moth and butterfly gut α-amylase optimum pH. Dose-dependent inhibition of two crucifer pests by triticale and wheat was observed using spectrophotometry and gel electrophoresis. Implementation of specificity studies showed that wheat and triticale-proteinaceous extract were inactive against Chinese and purple cabbage amylase. Triticale and wheat were resistant against insects' gut proteases. Results of the feeding bioassay indicated that triticale-proteinaceous extract could cause a significant reduction in survival and larval body mass. The results of the nutritional indecs also showed larvae of both species that fed on a Triticale proteinaceous inhibitor-treated diet had the lowest values for the efficiency of conversion of ingested food and relative growth rate. Our observations suggested that triticale shows promise for use in the management of crucifer pests.

Comparative transcriptome analysis provides insights of anti-insect molecular mechanism of Cassia obtusifolia trypsin inhibitor against Pieris rapae

Pieris rapae, a serious Lepidoptera pest of cultivated crucifers, utilizes midgut enzymes to digest food and detoxify secondary metabolites from host plants. A recombinant trypsin inhibitor (COTI) from nonhost plant, Cassia obtusifolia, significantly decreased activities of trypsin-like proteases in the larval midgut on Pieris rapae and could suppress the growth of larvae. In order to know how COTI took effect, transcriptional profiles of P. rapae midgut in response to COTI was studied. A total of 51,544 unigenes were generated and 45.86% of which had homologs in public databases. Most of the regulated genes associated with digestion, detoxification, homeostasis, and resistance were downregulated after ingestion of COTI. Meanwhile, several unigenes in the integrin signaling pathway might be involved in response to COTI. Furthermore, using comparative transcriptome analysis, we detected differently expressing genes and identified a new reference gene, UPF3, by qRT-polymerase chain reaction (PCR). Therefore, it was suggested that not only proteolysis inhibition, but also suppression of expression of genes involved in metabolism, development, signaling, and defense might account for the anti-insect resistance of COTI.

A possible genetic basis for vulnerability in Euphydryas maturna (Lepidoptera: Nymphalidae)

Nearly all of the known populations of Scarce Fritillary, Euphydryas maturna (Linnaeus, 1758), are declining in Western and Central Europe. In order to identify the possible reasons for its vulnerability we surveyed the population genetics of this butterfly species using multi-locus genotype data. Females of our target species lay lots of eggs in one or two batches only and pre-hibernation caterpillars live and feed gregariously in a nest. As a consequence, a random unfavourable event can eliminate most offspring of a particular female resulting in a strong genetic drift effect combined with inbreeding. Thus, our hypothesis regarding the genetic composition of Scarce Fritillary populations suggests that: (1) there will be random fluctuations in allele frequencies from generation to generation; (2) populations should exhibit small effective sizes and a relatively high level of heterozygote deficiency, and; (3) the majority of the individuals in a population will be composed of the offspring of just a few females. In order to test these hypotheses, fine-scale genetic structure was studied in two subpopulations of a Hungarian Scarce Fritillary population for 4 consecutive years (generations) using enzyme polymorphism data. The results supported all of our assumptions. We detected random fluctuation in the frequency of several alleles, small effective population size and the index of heterozygote deficiency (F _IS) indicated a considerable level of inbreeding in most samples. Furthermore, average values of relatedness were also fairly high, and we were able to identify 17 putative sib families in total with the two subpopulations based on estimation of individual gametic phases. Thus, the present study suggests that intrinsic factors (e.g. specific life history) might increase the sensitivity of a species to various threatening factors (e.g. habitat loss or fragmentation) and result in the vulnerability of the given species.

Phylogenetic incongruence and the evolutionary origins of cardenolide-resistant forms of Na(+) ,K(+) -ATPase in Danaus butterflies

Many distantly related insect species are specialized feeders of cardenolide-containing host plants such as milkweed (Asclepias spp.). Previous studies have revealed frequent, parallel substitution of a functionally important amino acid substitution (N122H) in the alpha subunit of Na(+) ,K(+) -ATPase in a number of these species. This substitution facilitates the ability of these insects to feed on their toxic hosts and sequester cardenolides for their own use in defense. Among milkweed butterflies of the genus Danaus, the previously established phylogeny for this group suggests that N122H arose independently and fixed in two distinct lineages. We reevaluate this conclusion by examining Danaus phylogenetic relationships using >400 orthologous gene sequences assembled from transcriptome data. Our results indicate that the three Danaus species known to harbor the N122H substitution are more closely related than previously thought, consistent with a single, common origin for N122H. However, we also find evidence of both incomplete lineage sorting and post-speciation genetic exchange among these butterfly species, raising the possibility of collateral evolution of cardenolide-insensitivity in this species group.

Changes is genes coding for laccases 1 and 2 may contribute to deformation and reduction of wings in apollo butterfly (Parnassius apollo, Lepidoptera: Papilionidae) from the isolated population in Pieniny National Park (Poland)

An isolated population of apollo butterfly (Parnassius apollo, Lepidoptera: Papilionidae) occurs in Pieniny National Park (Poland). Deformations and reductions of wings in a relatively large number of individuals from this population is found, yet the reasons for these defects are unknown. During studies devoted to identify cause(s) of this phenomenon, we found that specific regions of genes coding of enzymes laccases 1 and 2 could not be amplified from DNA samples isolated from large fractions of malformed insects while expected PCR products were detected in almost all (with one exception) normal butterflies. Laccases (p-diphenol:dioxygen oxidoreductases) are oxidases containing several copper atoms. They catalyse single-electron oxidations of phenolic or other compounds with concomitant reduction of oxygen to water. In insects, their enzymatic activities were found previously in epidermis, midgut, Malpighian tubules, salivary glands, and reproductive tissues. Therefore, we suggest that defects in genes coding for laccases might contribute to deformation and reduction of wings in apollo butterflies, though it seems obvious that deficiency in these enzymes could not be the sole cause of these developmental improperties in P. apollo from Pieniny National Park.

Pierisins and CARP-1: ADP-ribosylation of DNA by ARTCs in butterflies and shellfish

The cabbage butterfly, Pieris rapae, and related species possess a previously unknown ADP-ribosylating toxin, guanine specific ADP-ribosyltransferase. This enzyme toxin, known as pierisin, consists of enzymatic N-terminal domain and receptor-binding C-terminal domain, or typical AB-toxin structure. Pierisin efficiently transfers an ADP-ribosyl moiety to the N(2) position of the guanine base of dsDNA. Receptors for pierisin are suggested to be the neutral glycosphingolipids, globotriaosylceramide (Gb3), and globotetraosylceramide (Gb4). This DNA-modifying toxin exhibits strong cytotoxicity and induces apoptosis in various human cell lines, which can be blocked by Bcl-2. Pierisin also produces detrimental effects on the eggs and larvae of the non-habitual parasitoids. In contrast, a natural parasitoid of the cabbage butterfly, Cotesia glomerata, was resistant to this toxin. The physiological role of pierisin in the butterfly is suggested to be a defense factor against parasitization by wasps. Other type of DNA ADP-ribosyltransferase is present in certain kinds of edible clams. For example, the CARP-1 protein found in Meretrix lamarckii consists of an enzymatic domain without a possible receptor-binding domain. Pierisin and CARP-1 are almost fully non-homologous at the amino acid sequence level, but other ADP-ribosyltransferases homologous to pierisin are present in different biological species such as eubacterium Streptomyces. Possible diverse physiological roles of the DNA ADP-ribosyltransferases are discussed.

Cytochrome P450-encoding genes from the Heliconius genome as candidates for cyanogenesis

Cytochrome P450s are important both in the metabolism of xenobiotics and the production of compounds such as cyanogenic glucosides, which insects use in their defence. In the present study, we use transcriptomic and genomic information to isolate and name P450-encoding genes from the butterfly Heliconius melpomene. We classify each of the putative genes into its appropriate superfamily and compare the distribution of P450s across sequenced insects. We also identify homologues of two P450s known to be involved in cyanogenesis in the six-spot Burnet moth, Zygaena filipendulae. Classification of Heliconius P450s should be an important step in the dissection of their role in the exploitation of their host plant, the passion vine Passiflora.

A genetic polymorphism evolving in parallel in two cell compartments and in two clades

BACKGROUND

The enzyme phosphoenolpyruvate carboxykinase, PEPCK, occurs in its guanosine-nucleotide-using form in animals and a few prokaryotes. We study its natural genetic variation in Colias (Lepidoptera, Pieridae). PEPCK offers a route, alternative to pyruvate kinase, for carbon skeletons to move between cytosolic glycolysis and mitochondrial Krebs cycle reactions.

RESULTS

PEPCK is expressed in both cytosol and mitochondrion, but differently in diverse animal clades. In vertebrates and independently in Drosophila, compartment-specific paralogous genes occur. In a contrasting expression strategy, compartment-specific PEPCKs of Colias and of the silkmoth, Bombyx, differ only in their first, 5', exons; these are alternatively spliced onto a common series of following exons. In two Colias species from distinct clades, PEPCK sequence is highly variable at nonsynonymous and synonymous sites, mainly in its common exons. Three major amino acid polymorphisms, Gly 335 ↔ Ser, Asp 503 ↔ Glu, and Ile 629 ↔ Val occur in both species, and in the first two cases are similar in frequency between species. Homology-based structural modelling shows that the variants can alter hydrogen bonding, salt bridging, or van der Waals interactions of amino acid side chains, locally or at one another's sites which are distant in PEPCK's structure, and thus may affect its enzyme function. We ask, using coalescent simulations, if these polymorphisms' cross-species similarities are compatible with neutral evolution by genetic drift, but find the probability of this null hypothesis is 0.001 ≤ P ≤ 0.006 under differing scenarios.

CONCLUSION

Our results make the null hypothesis of neutrality of these PEPCK polymorphisms quite unlikely, but support an alternative hypothesis that they are maintained by natural selection in parallel in the two species. This alternative can now be justifiably tested further via studies of PEPCK genotypes' effects on function, organismal performance, and fitness. This case emphasizes the importance, for evolutionary insight, of studying gene-specific mechanisms affected by natural genetic variation as an essential complement to surveys of such variation.

Effect of milk thistle, Silybium marianum, extract on toxicity, development, nutrition, and enzyme activities of the small white butterfly, Pieris rapae

The methanolic extract of milk thistle, Silybium marianum L. (Asterales: Asteraceae), was investigated for its effects on the mortality, growth, feeding indices, enzymatic activity, and levels of non-enzymatic molecules of the small white butterfly, Pieris rapae L. (Lepidoptera: Pieridae), a pest of cruciferous plants. Feeding indices including approximate digestibility (AD), efficiency of conversion of digested food (ECD), efficiency of conversion of ingested food (ECI), relative growth rate (RGR), and relative consumption rate (RCR) were measured. These indices were variously affected: the RGR, RCR, and AD decreased, but the ECD and ECI increased. The LC50 and LC25 values were estimated as 2.94% and 1.20%, respectively. At the lowest concentration of S. marianum extract (0.625%), the feeding deterrence index was 40.48%. The duration of the pupal stage and the rate of larval growth decreased. These changes may be due to alterations in metabolic activity, such as the increase in alkaline phosphatase activity, which is likely involved in detoxification. Additionally, the activities of alanine aminotransferase and aspartate aminotransferase, which are key components of amino acid catabolism, decreased. The amount of glucose (an energy source) and uric acid (the excreted end product) increased, while total protein (another energy source) and cholesterol decreased. These results indicate that this plant possesses potential secondary metabolites that may be useful for the future study of the control of insect pests.

Efficient targeted mutagenesis in the monarch butterfly using zinc-finger nucleases

The development of reverse-genetic tools in "nonmodel" insect species with distinct biology is critical to establish them as viable model systems. The eastern North American monarch butterfly (Danaus plexippus), whose genome is sequenced, has emerged as a model to study animal clocks, navigational mechanisms, and the genetic basis of long-distance migration. Here, we developed a highly efficient gene-targeting approach in the monarch using zinc-finger nucleases (ZFNs), engineered nucleases that generate mutations at targeted genomic sequences. We focused our ZFN approach on targeting the type 2 vertebrate-like cryptochrome gene of the monarch (designated cry2), which encodes a putative transcriptional repressor of the monarch circadian clockwork. Co-injections of mRNAs encoding ZFNs targeting the second exon of monarch cry2 into "one nucleus" stage embryos led to high-frequency nonhomologous end-joining-mediated, mutagenic lesions in the germline (up to 50%). Heritable ZFN-induced lesions in two independent lines produced truncated, nonfunctional CRY2 proteins, resulting in the in vivo disruption of circadian behavior and the molecular clock mechanism. Our work genetically defines CRY2 as an essential transcriptional repressor of the monarch circadian clock and provides a proof of concept for the use of ZFNs for manipulating genes in the monarch butterfly genome. Importantly, this approach could be used in other lepidopterans and "nonmodel" insects, thus opening new avenues to decipher the molecular underpinnings of a variety of biological processes.

A digestive lipase of Pieris brassicae L. (Lepidoptera: Pieridae): purification, characterization, and host plants effects

The properties of a digestive lipase from the larval midgut of Pieris brassicae were studied by performing biochemical purification, characterization, effect of host plants, and extracted inhibitors. The purification process revealed a lipase with a purification fold of 42, recovery of 18.12%, molecular weight mass of 72.3 kDa, optimal pH at 11, and optimal temperature at 30°C, as well as stability at the optimal temperature for 12 h. The purified enzyme was inhibited by the ions Na(+), Mn(+), Fe(2+), and Cu(2+) and the inhibitors SDS, EDTA, TTHA, and mercaptoethanol. Ca(2+) and Mg(2+) increased activity of the purified lipase, but urea, PMSF, EGTA, and DTC had no effect on enzymatic activity. Feeding of larvae on three host plants, Trepaeolus majus, Brassica olearcea var. alba, and B. olearcea var. rubra revealed the highest lipase activity on T. majus, but the two varieties of B. olearcea significantly decreased lipase activity. Extraction of a crude inhibitor from two varieties of B. olearcea demonstrated that the crude inhibitor inhibited the purified lipase up to 75%. The inhibitor changed the kinetic parameters of the enzyme by elevating the K(m), as in competitive inhibition. The data suggest a possible role for plant lipase inhibitors in host plant resistance.

Phylogenetic relationships of true butterflies (Lepidoptera: Papilionoidea) inferred from COI, 16S rRNA and EF-1α sequences

The molecular phylogenetic relationships among true butterfly families (superfamily Papilionoidea) have been a matter of substantial controversy; this debate has led to several competing hypotheses. Two of the most compelling of those hypotheses involve the relationships of (Nymphalidae + Lycaenidae) + (Pieridae + Papilionidae) and (((Nymphalidae + Lycaenidae) + Pieridae) + Papilionidae). In this study, approximately 3,500 nucleotide sequences from cytochrome oxidase subunit I (COI), 16S ribosomal RNA (16S rRNA), and elongation factor-1 alpha (EF-1α) were sequenced from 83 species belonging to four true butterfly families, along with those of three outgroup species belonging to three lepidopteran superfamilies. These sequences were subjected to phylogenetic reconstruction via Bayesian Inference (BI), Maximum Likelihood (ML), and Maximum Parsimony (MP) algorithms. The monophyletic Pieridae and monophyletic Papilionidae evidenced good recovery in all analyses, but in some analyses, the monophylies of the Lycaenidae and Nymphalidae were hampered by the inclusion of single species of the lycaenid subfamily Miletinae and the nymphalid subfamily Danainae. Excluding those singletons, all phylogenetic analyses among the four true butterfly families clearly identified the Nymphalidae as the sister to the Lycaenidae and identified this group as a sister to the Pieridae, with the Papilionidae identified as the most basal linage to the true butterfly, thus supporting the hypothesis: (Papilionidae + (Pieridae + (Nymphalidae + Lycaenidae))).

Expression of immune-response genes in lepidopteran host is suppressed by venom from an endoparasitoid, Pteromalus puparum

BACKGROUND

The relationships between parasitoids and their insect hosts have attracted attention at two levels. First, the basic biology of host-parasitoid interactions is of fundamental interest. Second, parasitoids are widely used as biological control agents in sustainable agricultural programs. Females of the gregarious endoparasitoid Pteromalus puparum (Hymenoptera: Pteromalidae) inject venom along with eggs into their hosts. P. puparum does not inject polydnaviruses during oviposition. For this reason, P. puparum and its pupal host, the small white butterfly Pieris rapae (Lepidoptera: Pieridae), comprise an excellent model system for studying the influence of an endoparasitoid venom on the biology of the pupal host. P. puparum venom suppresses the immunity of its host, although the suppressive mechanisms are not fully understood. In this study, we tested our hypothesis that P. puparum venom influences host gene expression in the two main immunity-conferring tissues, hemocytes and fat body.

RESULTS

At 1 h post-venom injection, we recorded significant decreases in transcript levels of 217 EST clones (revealing 113 genes identified in silico, including 62 unknown contigs) derived from forward subtractive libraries of host hemocytes and in transcript levels of 288 EST clones (221 genes identified in silico, including 123 unknown contigs) from libraries of host fat body. These genes are related to insect immune response, cytoskeleton, cell cycle and apoptosis, metabolism, transport, stress response and transcriptional and translational regulation. We verified the reliability of the suppression subtractive hybridization (SSH) data with semi-quantitative RT-PCR analysis of a set of randomly selected genes. This analysis showed that most of the selected genes were down-regulated after venom injection.

CONCLUSIONS

Our findings support our hypothesis that P. puparum venom influences gene expression in host hemocytes and fat body. Specifically, the venom treatments led to reductions in expression of a large number of genes. Many of the down-regulated genes act in immunity, although others act in non-immune areas of host biology. We conclude that the actions of venom on host gene expression influence immunity as well as other aspects of host biology in ways that benefit the development and emergence of the next generation of parasitoids.

Ultrafast dynamics and anionic active states of the flavin cofactor in cryptochrome and photolyase

We report here our systematic studies of the dynamics of four redox states of the flavin cofactor in both photolyases and insect type 1 cryptochromes. With femtosecond resolution, we observed ultrafast photoreduction of oxidized state flavin adenine dinucleotide (FAD) in subpicosecond and of neutral radical semiquinone (FADH(*)) in tens of picoseconds through intraprotein electron transfer mainly with a neighboring conserved tryptophan triad. Such ultrafast dynamics make these forms of flavin unlikely to be the functional states of the photolyase/cryptochrome family. In contrast, we find that upon excitation the anionic semiquinone (FAD(*-)) and hydroquinone (FADH(-)) have longer lifetimes that are compatible with high-efficiency intermolecular electron transfer reactions. In photolyases, the excited active state (FADH(-)*) has a long (nanosecond) lifetime optimal for DNA-repair function. In insect type 1 cryptochromes known to be blue-light photoreceptors the excited active form (FAD(*-)*) has complex deactivation dynamics on the time scale from a few to hundreds of picoseconds, which is believed to occur through conical intersection(s) with a flexible bending motion to modulate the functional channel. These unique properties of anionic flavins suggest a universal mechanism of electron transfer for the initial functional steps of the photolyase/cryptochrome blue-light photoreceptor family.

Modifications in the N-terminus of an insect cytochrome P450 enhance production of catalytically active protein in baculovirus-Sf9 cell expression systems

Although baculovirus vectors are powerful tools for the heterologous expression of proteins in insect cell cultures, some insect and plant microsomal P450 proteins are not effectively expressed in this system. Hypothesizing that their expression failures might result from collisions between their N-terminal sequences and adjacent cytosolic sequences, we compared and mutated the N-terminus of Papilio multicaudatus CYP6B33, which is inappropriately folded in Sf9 cells, to sequences present in its Papilio polyxenes CYP6B1 counterpart, which is efficiently expressed and appropriately folded. Molecular modeling of the three differences in the linker separating the signal anchor domain (SAD) and the cytosolic domain identified Val32 in CYP6B33 as a residue potentially important for folding and/or positioning of the cytosolic domain. Mutation of Val32 to Ala32 in the CYP6B33 linker (CYP6B33 V32A mutant) or replacement of the CYP6B33 SAD with that of CYP6B1 (CYP6B1 1-20/CYP6B33 21-500 mutant) allowed for significant P450 expression, indicating that complex interactions involving both the signal anchor and membrane linker affect folding and activity of P450s in this heterologous expression system.

Regulation of 20-hydroxyecdysone on the larval pigmentation and the expression of melanin synthesis enzymes and yellow gene of the swallowtail butterfly, Papilio xuthus

The swallowtail butterfly, Papilio xuthus, changes its larval body pattern dramatically during the fourth ecdysis. Cuticular pigmentation occurs with precise timing just before ecdysis. We previously found that the cuticular pigmentation was regulated by three melanin synthesis genes, tyrosine hydroxylase (TH), dopa decarboxylase (DDC), and ebony. We discovered that yellow is strongly expressed in the presumptive black markings earlier than TH and DDC. Because the ecdysis is triggered by 20-hydroxyecdysone (20E), the effects of 20E on the pigmentation and expression of the melanin synthesis genes were examined. Here, we established a method for the topical application of 20E to molting specimens, so that 20E has only a partial effect, resulting in successful ecdysis. When we applied 20E during the mid-phase of the molting period, when the 20E titer is declining, cuticular pigmentation was completely inhibited. The cessation of hormonal treatments caused delayed pigmentation. yellow expression was promoted by a high titer of 20E, whereas the expression of TH, DDC, and ebony was suppressed, suggesting that a decline in the 20E concentration is necessary for the induction of the expression of the latter three genes. These results indicate that cuticular pigmentation is controlled by the exposure to 20E and its removal.

Cytochrome P450s in Papilio multicaudatus and the transition from oligophagy to polyphagy in the Papilionidae

Although substrate-specific CYP6B1 and CYP6B3 enzymes in Papilio polyxenes contribute to specialization on furanocoumarin-containing host plants, CYP6B4 and CYP6B17 enzymes in the polyphagous Papilio glaucus and Papilio canadensis have a broader range of substrates. Papilio multicaudatus, an oligophage with one furanocoumarin-containing host, is putatively ancestral to polyphagous Papilio species. Furanocoumarin-inducible CYP6B33-CYP6B37 and CYP6AB6 were characterized from this species. Heterologous expression of CYP6B33 revealed furanocoumarin metabolism resembling that of CYP6B4-CYP6B17 enzymes from P. glaucus and P. canadensis. Molecular models of CYP6B33 and CYP6B4 indicate that seven conserved aromatic side chains stabilize their hydrophobic catalytic sites and that a Lys484-Ser484 substitution enlarges the CYP6B4 active site pocket to increase the predicted distance between the substrate and reactive oxygen relative to CYP6B1. Loss of specialization in this lineage may have resulted from relatively few mutational changes, allowing acquisition of broader catalytic activities without loss of ancestral furanocoumarin-metabolizing activities.

Purification and some properties of β-N-acetyl-d-glucosaminidase from the cabbage butterfly (Pieris rapae)

A beta-N-acetyl-D-glucosaminidase (NAGase) from the cabbage butterfly (Pieris rapae) was purified. The purified enzyme was a single band on polyacrylamide gel electrophoresis and the specific activity was determined to be 8715 U/mg. The molecular weight of whole enzyme was determined to be 106 kDa by gel filtration, and the result of SDS-PAGE showed that the enzyme was a heterodimer, which contained two subunits with different mass of 59.5 and 57.2 kDa. The optimum pH and optimum temperature of the enzyme for the hydrolysis of p-nitrophenyl-N-acetyl-beta-D-glucosaminide (pNP-NAG) were investigated to be at pH 6.2 and at 42 degrees C, respectively, and the Michaelis-Menten constant (K(m)) was determined to be 0.285 mM at pH 6.2 and 37 degrees C. The stability of the enzyme was investigated and the results showed that the enzyme was stable at the pH range from 4.0 to 9.0 and at the temperature below 45 degrees C. The activation energy was 83.86 kJ/mol. The reaction of this enzyme with pNP-NAG was judged to be Ordered Bi-Bi mechanism according to the inhibitory behaviors of the products. The ionization constant, pK(e), of ionizing group at the active site of the enzyme was found to be 5.20 at 39.0 degrees C, and the standard dissociation enthalpy (DeltaH(o)) was determined to be 2.18 kcal/mol. These results showed that the ionizing group of the enzyme active center was the carboxyl group. The results of chemical modification also suggested that carboxyl group was essential to the enzyme activity. Moreover, Zn(2+), Hg(2+), Cu(2+) had strongly inhibitory effects on the enzyme activity.

Midgut protease activities in monophagous larvae of Apollo butterfly, Parnassius apollo ssp. frankenbergeri

We assayed the relative activities of midgut proteolytic enzymes in individuals of the fourth (L(4)) and fifth (L(5)) instar of Apollo larvae, inhabiting Pieniny Mts (southern Poland). The comparisons between midgut tissue with glicocalyx (MT) and liquid midgut contents with peritrophic membrane (MC) were made. Optimal media pHs of the assayed proteolytic enzymes in P. apollo midgut samples were similar to those of other lepidopteran species. Endopeptidases, as well as carboxypeptidases, digested effectively in alkaline environment, while aminopeptidases were active in a broad pH range. Trypsin is probably the main endoprotease (correlation with caseinolytic activity in MC of L(5) larvae: r=0.606; p=0.004); however, its activity was low as compared with that in other leaf-eating Lepidoptera. This suggests a minor role of trypsin and chymotrypsin in protein digestion in Apollo larvae, probably due to limited availability of the leaf proteins. Instead, due to very high carboxypeptidase A activity in midgut tissue, the larvae obtain exogenous amino acids either directly or from oligopeptides and glycoproteins. High and significant positive correlations between the enzyme activity and glucosidase as well as galactosidase activities strongly support this opinion.

Evidence of protease in the saliva of the butterfly Heliconius melpomene (L.) (Nymphalidae, Lepidoptera)

Butterflies of the genus Heliconius are well known for their peculiar habits of utilizing pollen as a source of amino acids. Saliva plays a major role in the process of extracting amino acids and proteins from the pollen grains. In this investigation, we obtained samples of saliva from adult Heliconius melpomene by placing pumpkin pollen or fine glass-beads on the proboscis, which stimulates the butterflies to release saliva. Proteolytic activity was determined in the saliva by an insoluble protein-dye that turns blue when cleaved by proteases. Its extinction value was measured with a spectrophotometer at 595 nm. Both the saliva sampled with pollen and the saliva obtained from inert glass-beads exhibit proteolytic activity demonstrating that the saliva contains proteases. The proteolytic activity of the pollen/saliva samples was higher than that of the glass-bead/saliva samples, which we attribute to the stimulating effects of pollen, such as taste, smell, and texture, and not to proteases which might have been liberated from the pollen. This is indicated by the fact that pollen samples without saliva showed only a negligible indication for proteolytic activity. In general, females exhibit higher proteolytic activities than males, presumably due to their greater amino acid investment in reproduction. We present here first evidence for the existence of proteases in the saliva of a butterfly species and suggest that these enzymes are crucial for the use of amino acids and proteins from pollen in Heliconius butterflies.

Genetic diversity in butterflies: Interactive effects of habitat fragmentation and climate-driven range expansion

Some species are expanding their ranges polewards during current climate warming. However, anthropogenic fragmentation of suitable habitat is affecting expansion rates and here we investigate interactions between range expansion, habitat fragmentation and genetic diversity. We examined three closely related Satyrinae butterflies, which differ in their habitat associations, from six sites along a transect in England from distribution core to expanding range margin. There was a significant decline in allozyme variation towards an expanding range margin in Pararge aegeria, which has the most restricted habitat availability, but not in Pyronia tithonus whose habitat is more widely available, or in a non-expanding 'control species' (Maniola jurtina). Moreover, data from another transect in Scotland indicated that declines in genetic diversity in P. aegeria were evident only on the transect in England, which had greater habitat fragmentation. Our results indicate that fragmentation of breeding habitats leads to more severe founder events during colonization, resulting in reduced diversity in marginal populations in more specialist species. The continued widespread loss of suitable habitats in the future may increase the likelihood of loss of genetic diversity in expanding species, which may affect whether or not species can adapt to future environmental change.

Porphobilinogen synthase from the butterfly, Pieris brassicae: purification and comparative characterization

Porphobilinogen represents a key building block of tetrapyrroles serving as functional ligands of many vitally important proteins. Here we report the first purification of porphobilinogen synthase (PBGS) from whole insects by sequentially employing two modes of native electrophoresis on polyacrylamide gels subsequent to more conventional procedures. Using adults of Pieris brassicae L. (Lepidoptera: Pieridae) we achieved approximately 10,000-fold purification with final yields of up to 25% of electrophoretically pure PBGS with a specific activity of approximately 160 micromol PBG h(-1) mg(-1) at 37 degrees C and an affinity of 0.36 mM to its substrate 5-aminolevulinic acid. Enzyme activity was inhibited by the substrate mimics, levulinic acid and succinylacetone, and by chelating agents. PBGS behaved as a relatively heat-stable octameric complex of 292.3 kDa composed of 36.5 kDa subunits. Most general features of this insect PBGS were comparable to those published for other animal PBGS enzymes, while remarkable differences were found to the reported recombinant Drosophila enzyme. Moreover, rabbit antiserum directed against purified Pieris PBGS revealed significant immunological differences among insect PBGS enzymes from a wide range of orders contrasting to the overall evolutionary conserved features of this enzyme.

A candidate locus for variation in dispersal rate in a butterfly metapopulation

Frequent extinctions of local populations in metapopulations create opportunities for migrant females to establish new populations. In a metapopulation of the Glanville fritillary butterfly (Melitaea cinxia), more mobile individuals are more likely to establish new populations, especially in habitat patches that are poorly connected to existing populations. Here we show that flight metabolic rate and the frequency of a specific allele of the metabolic enzyme phosphoglucose isomerase (pgi) were both highest in newly established, isolated populations. Furthermore, genotypes with this pgi allele had elevated flight metabolic rates. These results suggest that genetic variation in pgi or a closely linked locus has a direct effect on flight metabolism, dispersal rate, and thereby on metapopulation dynamics in this species. These results also contribute to an emerging understanding of the mechanisms by which population turnover in heterogeneous landscapes may maintain genetic and phenotypic variation across populations.

Molecular-level variation affects population growth in a butterfly metapopulation

The dynamics of natural populations are thought to be dominated by demographic and environmental processes with little influence of intraspecific genetic variation and natural selection, apart from inbreeding depression possibly reducing population growth in small populations. Here we analyse hundreds of well-characterised local populations in a large metapopulation of the Glanville fritillary butterfly (Melitaea cinxia), which persists in a balance between stochastic local extinctions and recolonisations in a network of 4,000 discrete habitat patches. We show that the allelic composition of the glycolytic enzyme phosphoglucose isomerase (Pgi) has a significant effect on the growth of local populations, consistent with previously reported effects of allelic variation on flight metabolic performance and fecundity in the Glanville fritillary and Colias butterflies. The strength and the sign of the molecular effect on population growth are sensitive to the ecological context (the area and spatial connectivity of the habitat patches), which affects genotype-specific gene flow and the influence of migration on the dynamics of local populations. The biological significance of the results for Pgi is underscored by lack of any association between population growth and allelic variation at six other loci typed in the same material. In demonstrating, to our knowledge for the first time, that molecular variation in a candidate gene affects population growth, this study challenges the perception that differential performance of individual genotypes, leading to differential fitness, is irrelevant to population dynamics. These results also demonstrate that the spatial configuration of habitat and spatial dynamics of populations contribute to maintenance of Pgi polymorphism in this species.

Researchers observe a link between the genetic composition of the Pgi locus (an enzyme linked with flight metabolic performance) and local population growth within a metapopulation of butterflies.

Inhibition of CYP6B1-mediated detoxification of xanthotoxin by plant allelochemicals in the black swallowtail (Papilio polyxenes)

The structural and biosynthetic diversity of allelochemicals in plants is thought to arise from selection for additive toxicity as a consequence of toxin mixture or for enhanced toxicity as a result of synergism. In order to understand how insects cope with this type of plant defense, we tested the effects of some allelochemicals in host plants of the black swallowtail Papilio polyxenes on the xanthotoxin-metabolic activity of CYP6B1, the principal enzyme responsible for the detoxification of furanocoumarins in this caterpillar. Additionally, the effects of some synthetic compounds not normally encountered by P. polyxenes on CYP6B1 were tested. These studies demonstrate that the integrity of furanocoumarin structure is important for competitive binding to the active site of CYP6B1, even though the carbonyl group on the pyranone ring apparently does not affect its inhibitory capacity, as in the case of furanochromones. Angular furanocoumarins are generally less phototoxic to many organisms than linear furanocoumarins due to their reduced capacity for cross-linking DNA strands, yet they are more toxic than linear furanocoumarins to black swallowtail larvae. This enhanced toxicity in vivo may be due to the ability of angular furanocoumarins to bind to the active site of CYP6B1 without being rapidly metabolized. This binding reduces the availability of CYP6B1 to metabolize other linear furanocoumarins. The structure-activity relationships for methylenedioxyphenyl compounds, flavonoids, imidazole, and imidazole derivatives are also discussed in light of their capacity to inhibit the xanthotoxin-metabolic activity of CYP6B1.

From DNA to fitness differences: sequences and structures of adaptive variants of Colias phosphoglucose isomerase (PGI)

Colias eurytheme butterflies display extensive allozyme polymorphism in the enzyme phosphoglucose isomerase (PGI). Earlier studies on biochemical and fitness effects of these genotypes found evidence of strong natural selection maintaining this polymorphism in the wild. Here we analyze the molecular features of this polymorphism by sequencing multiple alleles and modeling their structures. PGI is a dimer with rotational symmetry. Each monomer provides a critical residue to the other monomer's catalytic center. Sequenced alleles differ at multiple amino acid positions, including cryptic charge-neutral variation, but most consistent differences among the electromorph alleles are at the charge-changing amino acid sites. Principal candidate sites of selection, identified by structural and functional analyses and by their variants' population frequencies, occur in interpenetrating loops across the interface between monomers, where they may alter subunit interactions and catalytic center geometry. Comparison to a second (and basal) species, Colias meadii, also polymorphic for PGI under natural selection, reveals one fixed amino acid difference between their PGIs, which is located in the interpenetrating loop and accompanies functional differences among their variants. We also study nucleotide variability among the PGI alleles, comparing these data to similar data from another glycolytic enzyme gene, glyceraldehyde-3-phosphate dehydrogenase. Despite extensive nonsynonymous and synonymous polymorphism at PGI in each species, the only base changes fixed between species are the two causing the amino acid replacement; this absence of synonymous fixation yields a significant McDonald-Kreitman test. Analyses of these data suggest historical population expansion. Positive peaks of Tajima's D statistic, representing regions of neutral "hitchhiking," are found around the principal candidate sites of selection. This study provides novel views of molecular-structural mechanisms, and beginnings of historical evidence, for a long-persistent balanced enzyme polymorphism at PGI in these and perhaps other species.

Disjunct distributions during glacial and interglacial periods in mountain butterflies: Erebia epiphron as an example

Over several decades, the distribution patterns and evolution of alpine disjunct species has become an increasingly discussed subject. Large scale genetic analysis has allowed the resolution of the past range changes and intraspecific evolution of many species, in Europe especially of Mediterranean origin. However, the phylogeographic structures of species with arctic-alpine disjunct distribution patterns are relatively poorly studied. The existing phylogeographic analysis (mostly of alpine plant species) supports disjunct distributions during glacial as well as post-glacial periods for a number of species. However, several questions still remain unresolved and we therefore analysed the Mountain Ringlet Erebia epiphron as a model for such alpine disjunct species. We found strong differentiation into five different lineages supporting five differentiation centres: (i) the eastern Pyrenees, (ii) the mountain ranges between the central Pyrenees and south-western Alps, (iii, iv) two areas along the southern Alps margin and (v) the northern Alps margin. We propose that these patterns evolved due to the humidity requirements of this species, which did not allow survival in the dry glacial steppes, but along the margins of the wetter glaciated high mountain ranges.

Expression of one isoform of GTP cyclohydrolase I coincides with the larval black markings of the swallowtail butterfly, Papilio xuthus

The larva of the swallowtail butterfly Papilio xuthus changes its body markings during the fourth ecdysis. We found that stage-specific cuticular black markings are mainly regulated by co-localization of two melanin synthesis enzymes; tyrosine hydroxylase (TH) and dopa decarboxylase (DDC). TH converts tyrosine to dihydroxyphenylalanine (dopa), and tyrosine itself is converted from phenylalanine by phenylalanine hydroxylase (PAH). Guanosine triphosphate cyclohydrolase I (GTPCHI) is essential for the synthesis of tetrahydrobiopterin (BH4) that is a cofactor of TH and PAH. In this report, we found that a GTPCHI inhibitor prevents pigmentation in cultured integuments, suggesting that the GTPCHI activity is also involved in cuticle pigmentation. We have cloned GTPCHI and PAH cDNAs from P. xuthus and investigated their spatial expression patterns in epidermis by whole-mount in situ hybridization. There are two isoforms of GTPCHI in larval epidermis (GTPCHIa and GTPCHIb). GTPCHIa is expressed at the black markings of the subsequent instar, similar to TH, whereas GTPCHIb is expressed uniformly, similar to PAH. This suggests that the region-specific expression of GTPCHIa supplies sufficient BH(4) reinforcing the TH activity in black marking area. Our results imply that larval markings are regulated by not only melanin synthesis enzymes but also the cofactor supplying enzyme.

Melanin-synthesis enzymes coregulate stage-specific larval cuticular markings in the swallowtail butterfly, Papilio xuthus

Like the adult wing, butterfly larvae are unique in their coloring. However, the molecular mechanisms underlying the formation of insect larval color patterns are largely unknown. The larva of the swallowtail butterfly Papilio xuthus changes its color pattern markedly during the 4th ecdysis. We investigated its cuticular color pattern, which is thought to be composed of melanin and related pigments derived from tyrosine. We cloned three enzymes involved in the melanin-synthesis pathway in P. xuthus: tyrosine hydroxylase (TH), dopa decarboxylase (DDC), and ebony. Whole-mount in situ hybridization showed that the expression of both TH and DDC is strongly correlated with the black markings. ebony is strongly expressed only in the reddish-brown area. The expression pattern of each enzyme coincides with the cuticular color pattern of the subsequent instar. We also investigated the uptake of melanin precursors into cultured integument. Inhibition of either TH or DDC activity prevents in vitro pigmentation completely. Addition of dopamine to integuments in the presence of TH inhibitor causes overall darkening without specific markings. From these results, specific larval cuticular color patterns are regulated by stage-specific colocalization of enzymes in epidermal cells rather than by the differential uptake of melanin precursors into individual epidermal cells. Epidermal cells expressing TH and DDC, but not ebony, produce the black cuticle, and epidermal cells expressing TH, DDC, and ebony produce the reddish-brown cuticle.

Phylogenetic relationships of the New World Troidini swallowtails (Lepidoptera: Papilionidae) based on COI, COII, and EF-1α genes

A phylogeny of the Neotropical members of the Tribe Troidini (Lepidoptera: Papilionidae) was obtained with sequences of three protein-coding genes: two mitochondrial (COI and COII), and one nuclear (EF-1alpha). Parsimony and Bayesian analyses of 33 taxa resulted in very similar trees regardless of method used with the 27 troidines always forming a monophyletic clade. Within Troidini, the genus Battus is sister group to the remaining troidines, followed by a clade formed by the Paleotropical taxa (here represented by three exemplars). The genus Euryades is the next branch, and sister group of Parides. The genus Parides is monophyletic, and is divided into four main groups by Maximum Parsimony analysis, with the most basal group composed of tailed species restricted to SE Brazil. Character optimization of ecological and morphological traits over the phylogeny proposed for troidines indicated that the use of several species of Aristolochia is ancestral over the use of few or a single host-plant. For the other three characters, the ancestral states were the absence of long tails, forest as the primary habitat and oviposition solitary or in loose group of several eggs.

Regulation of an insect cytochrome P450 monooxygenase gene (CYP6B1) by aryl hydrocarbon and xanthotoxin response cascades

Many organisms respond to toxic compounds in their environment by inducing regulatory networks controlling the expression and activity of cytochrome P450 monooxygenase (P450s) detoxificative enzymes. In particular, black swallowtail (Papilio polyxenes) caterpillars respond to xanthotoxin, a toxic phytochemical in their hostplants, by activating transcription of the CYP6B1 promoter via several regions located within 150 nt of the transcription initiation site. One such element is the xenobiotic response element to xanthotoxin (XRE-Xan) that lies upstream of consensus XRE-AhR (xenobiotic response element to the aryl hydrocarbon receptor) and OCT-1 (octamer-1 binding site) element known to be utilized in mammalian aryl hydrocarbon response cascades. Two-plasmid transfections conducted in Sf9 cells have indicated that XRE-Xan, XRE-AhR and a number of other proximal elements, but not OCT-1, are critical for basal as well as xanthotoxin- and benzo[alpha]pyrene-induced transcription of the CYP6B1 promoter. Four-plasmid transfections with vectors co-expressing the spineless (Ss) and tango (Tgo) proteins, the Drosophila melanogaster homologues of mammalian AhR and ARNT, have indicated that these proteins enhance basal expression of the CYP6B1 promoter but not the magnitude of its xanthotoxin and benzo[alpha]pyrene induction. Based on these results, we propose that these Drosophila transcription factors modulate basal expression of this promoter in a ligand-independent manner and attenuate its subsequent responses to planar aryl hydrocarbons (benzo[alpha]pyrene) and allelochemicals (xanthotoxin).

Identification of cytochrome P450 and glutathione-S-transferase genes preferentially expressed in chemosensory organs of the swallowtail butterfly, Papilio xuthus L

The swallowtail butterfly, Papilio xuthus L., feeds exclusively on members of the plant family, Rutaceae. Female butterflies lay eggs in response to specific chemicals contained in their host plants. They perceive a variety of polar compounds as oviposition stimulants through the tarsal chemosensilla of the foreleg by drumming upon the leaf surface. We undertook an expressed sequence tag (EST) analysis to identify the chemosensory-related genes that are expressed in chemosensilla on the tarsus of P. xuthus. Several genes that showed similarity with biotransformation enzymes were identified from the ESTs. Among them, a cytochrome P450 and a glutathione-S-transferase (GST) were preferentially expressed in the chemosensory organs. We have determined the structure of both cDNA and genomic sequences encoding these enzymes and designated the P450 as CYP341A2, a novel member of CYP341A subfamily, and the GST as GST-pxcs1, respectively. We observed a localized expression of CYP341A2 at the base of tarsal chemosensilla by in situ hybridization. These results suggest that these degrading enzymes play a role in the chemosensory reception for host plant recognition.

Ile115Leu mutation in the SRS1 region of an insect cytochrome P450 (CYP6B1) compromises substrate turnover via changes in a predicted product release channel

CYP6B1 represents the principal cytochrome P450 monooxygenase responsible for metabolizing furanocoumarins in Papilio polyxenes, an insect that specializes on host plants containing these toxins. Investigations of the amino acids responsible for the efficient metabolism of these plant toxins has identified Ile115 as one that modulates the rate of furanocoumarin metabolism even though it is predicted to be positioned at the edge of the heme plane and outside substrate contact regions. In contrast to previous expression studies conducted under conditions of limiting P450 reductase showing that the Ile115-to-Leu replacement enhances turnover of xanthotoxin and other furanocoumarins, studies conducted at high P450 reductase indicate that the Ile115-to-Leu replacement reduces turnover of these substrates. Further analysis of substrate binding affinities, heme spin state and NADPH consumption rates indicate that, whereas the I115L replacement mutant displays higher substrate affinity and heme spin state than the wild-type CYP6B1 protein, it utilizes NADPH more slowly than the wild-type CYP6B1 protein at high P450 reductase levels. Molecular models developed for the wild-type CYP6B1 and mutant protein suggest that more constricted channels extending from the catalytic site in the I115L mutant to the P450 surface limit the rate of product release from this mutant catalytic site under conditions not limited by the rate of electron transfer from NADPH.

Enzymatic properties of pierisin-1 and its N-terminal domain, a guanine-specific ADP-ribosyltransferase from the cabbage butterfly

The cabbage butterfly, Pieris rapae, produces an ADP-ribosylating cytotoxic protein, pierisin-1. Unlike other ADP-ribosylating toxins, the acceptor site for ADP-ribosylation by pierisin-1 is the N-2 position of guanine bases in DNA. The present study was designed to characterize this novel guanine-specific ADP-ribosyltransferase, pierisin-1. The N-terminal polypeptide from Met-1 to Arg-233, but not the C-terminal Ser-234-Met-850 polypeptide, was found to exhibit guanine ADP-ribosyltransferase activity. Trypsin-treated pierisin-1, which is considered to be a "nicked" full-length form composed of associated N- and C-terminal fragments, also demonstrated such activity. Optimum conditions for the N-terminal polypeptide of pierisin-1 were pH 8-10, 37-40 degrees C, in the presence of 100-200 mM NaCl or KCl. Other metal ions such as Ca(2+) or Mg(2+) were not required. Kinetic studies demonstrated potent ADP-ribosyltransferase activity with a K(M) value for NAD of 0.17 mM and k(cat) of 55 per second. Under these optimum conditions, the specific activity of trypsin-treated pierisin-1 was about half (k(cat) = 25 per second). When the conditions were changed to pH 5-7 or 10-20 degrees C, some activity (6-55% or 5-20%, respectively, of that under optimal conditions) of the N-terminal polypeptide was still evident; however, almost all of the trypsin-treated enzyme activity disappeared. This implies the inhibition of the N-terminal enzyme domain by the associated C-terminal fragment. Long-term reactions indicated that a single molecule of pierisin-1 has the capacity to generate more than 10(6) ADP-ribosylated DNA adducts, which could cause the death of a mammalian cell.

Conserved regulatory elements in the promoters of two allelochemical-inducible cytochrome P450 genes differentially regulate transcription

CYP6B4, a cytochrome P450 gene from the tiger swallowtail Papilio glaucus, is transcriptionally induced in the midgut by dietary furanocoumarins, plant allelochemicals that can crosslink DNA in their UV-activated form. The CYP6B4 promoter contains an overlapping EcRE/ARE/XRE-xan element similar to that used for basal and xanthotoxin-inducible expression of the CYP6B1 promoter from the black swallowtail Papilio polyxenes. Transfection of the CYP6B4 promoter:CAT reporter construct into Sf9 cells demonstrates that the basal and xanthotoxin-inducible expression levels observed reflect the relative expression levels of this gene in the midguts of tiger swallowtail larvae. Transfections of mutant CYP6B4 promoter constructs into Sf9 cells indicate that the EcRE/ARE/XRE-xan element is necessary for CYP6B4 induction by xanthotoxin but not for its minimal basal expression. In addition to these elements, the CYP6B4 and CYP6B1 promoters also contain putative XRE-AhR elements identical to the aryl hydrocarbon response elements present in mammalian phase I detoxification genes. Transfections of CYP6B4 and CYP6B1 promoters containing EcRE/ARE/XRE-xan and XRE-AhR elements indicate that both are induced significantly by benzo(alpha)pyrene, an aryl hydrocarbon widespread in the environment, as well as by xanthotoxin, an allelochemical encountered in their hostplants.

Purification and properties of isocitrate lyase from pupas of the butterfly Papilio machaon L

Key enzymes of the glyoxylate cycle, isocitrate lyase and malate synthase, were identified in pupas of the butterfly Papilio machaon L. The activities of these enzymes in pupas were 0.056 and 0.108 unit per mg protein, respectively. Isocitrate lyase was purified by a combination of various chromatographic steps including ammonium sulfate fractionation, ion-exchange chromatography on DEAE-Toyopearl, and gel filtration. The specific activity of the purified enzyme was 5.5 units per mg protein, which corresponded to 98-fold purification and 6% yield. The enzyme followed Michaelis-Menten kinetics (Km for isocitrate, 1.4 mM) and was competitively inhibited by succinate (Ki = 1.8 mM) and malate (Ki = 1 mM). The study of physicochemical properties of the enzyme showed that it is a homodimer with a subunit molecular weight of 68 +/- 2 kD and a pH optimum of 7.5 (in Tris-HCl buffer).

Diversification of furanocoumarin-metabolizing cytochrome P450 monooxygenases in two papilionids: Specificity and substrate encounter rate

Diversification of cytochrome P450 monooxygenases (P450s) is thought to result from antagonistic interactions between plants and their herbivorous enemies. However, little direct evidence demonstrates the relationship between selection by plant toxins and adaptive changes in herbivore P450s. Here we show that the furanocoumarin-metabolic activity of CYP6B proteins in two species of swallowtail caterpillars is associated with the probability of encountering host plant furanocoumarins. Catalytic activity was compared in two closely related CYP6B4 and CYP6B17 groups in the polyphagous congeners Papilio glaucus and Papilio canadensis. Generally, P450s from P. glaucus, which feeds occasionally on furanocoumarin-containing host plants, display higher activities against furanocoumarins than those from P. canadensis, which normally does not encounter furanocoumarins. These P450s in turn catalyze a larger range of furanocoumarins at lower efficiency than CYP6B1, a P450 from Papilio polyxenes, which feeds exclusively on furanocoumarin-containing host plants. Reconstruction of the ancestral CYP6B sequences using maximum likelihood predictions and comparisons of the sequence and geometry of their active sites to those of contemporary CYP6B proteins indicate that host plant diversity is directly related to P450 activity and inversely related to substrate specificity. These predictions suggest that, along the lineage leading to Papilio P450s, the ancestral, highly versatile CYP6B protein presumed to exist in a polyphagous species evolved through time into a more efficient and specialized CYP6B1-like protein in Papilio species with continual exposure to furanocoumarins. Further diversification of Papilio CYP6Bs has likely involved interspersed events of positive selection in oligophagous species and relaxation of functional constraints in polyphagous species.

Metabolism of linear and angular furanocoumarins by Papilio polyxenes CYP6B1 co-expressed with NADPH cytochrome P450 reductase

One challenge in the heterologous expression of microsomal cytochrome P450 monooxygenases (P450s) is fulfilling their obligatory requirement for electrons transferred from NADPH P450 reductase. We have established co-expression parameters for Papilio polyxenes CYP6B1 and house fly P450 reductase in baculovirus-infected Sf9 cells that allow for efficient expression of both components and significantly enhance metabolic turnover of this insect P450's substrates. These expression conditions have allowed us to reexamine the turnover capacities of CYP6B1 toward linear and angular furanocoumarins present in the host plants for the specialist caterpillar P. polyxenes. Coexpression of CYP6B1 and P450 reductase at equivalent viral concentrations [MOI (multiplicity of infection) ratio of 1] results in turnover rates for the linear furanocoumarins xanthotoxin and psoralen, which are increased 32-33 fold over the turnover rates obtained with CYP6B1 expressed alone. The turnover rate for the angular furanocoumarin angelicin is also significantly increased to 4.76 nmol/min/nmol P450 compared to its barely detectable level obtained with CYP6B1 expressed alone. Substrate binding analyses indicate that all three of these compounds elicit typical type I binding spectra but with varying magnitudes and affinities that are indicative of each substrate's effectiveness at coordinating with the heme iron. The relative proportions of high spin state generated with these substrates are consistent with CYP6B1 metabolizing these furanocoumarins in the rank order xanthotoxin>psoralen>angelicin. These differential activities for CYP6B1 suggest that it may have been an ancient participant in the coevolutionary arms race between papilionid butterflies and their apiaceous host plants. Due to its ability to handle a range of furanocoumarin structures, CYP6B1 may have contributed to P. polyxenes' early colonization of linear furanocoumarin-containing plants and to its subsequent colonization of angular furanocoumarin-containing plants.

CYP6B cytochrome p450 monooxygenases from Papilio canadensis and Papilio glaucus: potential contributions of sequence divergence to host plant associations

Two groups of furanocoumarin-inducible cytochrome p450 genes, the CYP6B4 group and the CYP6B17 group, characterized in two closely related tiger swallowtails, Papilio glaucus and Papilio canadensis, are induced to different extents, with generally higher levels of CYP6B transcripts in P. glaucus. To investigate the evolutionary history of these CYP6B genes in the context of their association with furanocoumarin detoxification, we isolated thirteen CYP6B genes from these species. Each of these genes contains an intron at a conserved position (1334 nucleotides from the translation start site), which varies in length due to three insertion/deletions. The proximal 5' end flanking sequence from the transcription initiation site is highly conserved (91-98% nt identity). The sequence 5' to -640 is significantly variable due largely to the presence of three insertion/deletions. The sequence at the 3' end of this region contains a putative xenobiotic response element to xanthotoxin (XRE-xan), important for basal and xanthotoxin-inducible transcription of the P. polyxenes CYP6B1v3 gene, and multiple elements known to regulate vertebrate phase I and II promoters, including an XRE-AhR (Xenobiotic Response Element to Aryl hydrocarbon Receptor), an OCT-1 element (octamer protein binding site), an ARE (Antioxidant Response Element), an EcRE (Ecdysone Response Element), and an imperfect PXR (Pregnane X Receptor) responsive element (PRE). Our analyses of CYP6B genes in these two species indicate that these genes are in an early stage of divergence and that differential exposure of these two species to chemically distinct host plants resulting from geographical isolation has had functional impacts not only on the coding regions of these genes but also on their promoter regions. Thus, changes in p450 regulation as well as catalytic activity may play a role in the evolution of host plant associations in herbivorous insects.

Hormonal control of GTP cyclohydrolase I gene expression and enzyme activity during color pattern development in wings of Precis coenia

Color patterns of butterfly wings are composed of single color points represented by each scale. In the case of Precis coenia, at the end of pupal development, different types of pigments are synthesized sequentially in the differently colored scales beginning with white (pterins) followed by red (ommatins) and then black (melanin). In order to explain how formation of these different colors is regulated, we examined the expression of an mRNA-encoding guanosine triphosphate-cyclohydrolase I (GTP-CH I; EC 3.5.4.16), the first key enzyme in the biosynthesis of pteridines, during pigment formation in the wings of P. coenia. The strongest positive signal was recognized around pigment formation one day before butterfly emergence. This GTP-CH I gene expression is paralleled by GTP-CH I enzyme activity measured in wing extracts. We also investigated the effect of 20-hydroxyecdysone on the expression of GTP-CH I mRNA and the enzyme activity during color formation. The results strongly suggest that the onset and duration of the expression of a GTP-CH I mRNA is triggered by a declining ecdysteroid hormone titer during late pupal development.

Herbivory: caterpillar saliva beats plant defences

Blood-feeding arthropods secrete special salivary proteins that suppress the defensive reaction they induce in their hosts. This is in contrast to herbivores, which are thought to be helpless victims of plant defences elicited by their oral secretions. On the basis of the finding that caterpillar regurgitant can reduce the amount of toxic nicotine released by the tobacco plant Nicotiana tabacum, we investigate here whether specific salivary components from the caterpillar Helicoverpa zea might be responsible for this suppression. We find that the enzyme glucose oxidase counteracts the production of nicotine induced by the caterpillar feeding on the plant.

Amino acids in SRS1 and SRS6 are critical for furanocoumarin metabolism by CYP6B1v1, a cytochrome P450 monooxygenase

CYP6B1v1 is the principal cytochrome P450 monooxygenase (P450) that detoxifies dietary furanocoumarins in the guts of Papilio polyxenes, the black swallowtail caterpillar. Sequence alignments and structure comparisons of CYP6B1v1 with the mouse CYP2A5 and bacterial CYP102 proteins, which are also capable of metabolizing the linear furanocoumarin xanthotoxin (8-methoxypsoralen), suggested that Phe116, His117, Val368 and Phe484 might be active site residues. In a homology model developed for CYP6B1v1, the side chains of Phe116 and His117 located in the B'-C loop of SRS1 are predicted to be positioned above the haem plane, while the side chain of Phe484 located in SRS6 is predicted near the entrance of the catalytic pocket. Site-directed mutagenesis of residues Phe116, His117 and Phe484 indicated that these residues represent several of those that determine this protein's stability and substrate specificity. Whereas all aromatic mutants of Phe116 and Phe484 generated CO-difference spectra with maxima at 450 nm indicative of correctly configured monooxygenases, aromatic mutants of Phe116 exhibited reduced reactivities toward some furanocoumarins and aromatic mutants of Phe484 eliminated all reactivities toward furanocoumarins. All single and double aliphatic mutants of Phe116, His117 and Phe484 and aromatic mutants of His117 generated carbon monoxide (CO) difference spectra with maxima at 420 nm (P420) indicative of incorrectly configured monooxygenases. These studies define residues Phe116, His117 and Phe484 as determinants of this insect P450's catalytic site integrity and residues Phe116 and Phe484 as determinants of its substrate specificity. Conservation of Phe116 and His117 in an array of lepidopteran CYP6B proteins implies that these amino acids serve a similar function in other monooxygenases of the insect CYP6B subfamily.

Developmental variation in cytochrome P450 expression in Papilio polyxenes in response to xanthotoxin, a hostplant allelochemical

Although developmental variation in activity and inducibility is typical of cytochrome P450 monooxygenases (P450s) in insects, the adaptive significance of such variation is often unclear, in part because the natural function of insect P450s is rarely known. In this study, we examined developmental variation in expression of CYP6B1 and CYP6B3 in Papilio polyxenes, the black swallowtail. Enzymes encoded by these genes have been implicated in the metabolism of xanthotoxin, a furanocoumarin characteristic of the apiaceous hostplants of P. polyxenes. In each life stage-egg, five larval instars, pupa, and adult-we examined individuals exposed to foliage with and without supplemental xanthotoxin. For each stage, we conducted enzyme assays to estimate xanthotoxin metabolism, Northern analysis to detect constitutive and induced mRNA levels, and RT-PCR amplification and Southern analysis to differentiate among P450 genes expressed. Inducible xanthotoxin metabolism, previously reported in fifth instars, was observed in four of five larval stages but was absent or undetectable in all stages that do not feed on foliage; the highest levels of activity were in early larval instars. The same pattern was observed in both Northern and RT-PCR gel blot analyses. In inducible larval stages, inducibility of CYP6B1 transcripts by xanthotoxin was greater than the inducibility of CYP6B3 transcripts. These findings support earlier suggestions that these two P450s contribute to xanthotoxin metabolism in this species and that expression of these P450 genes is regulated in an adaptive fashion with respect to probability of exposure to hostplant toxins over the course of development.

Molecular analysis of multiple CYP6B genes from polyphagous Papilio species

Papilio glaucus (eastern tiger swallowtail) and Papilio. canadensis (Canadian tiger swallowtail) are two closely related species with broad but overlapping hostplant ranges. P. glaucus encounters toxic furanocoumarins occasionally in its diet in its rutaceous hostplants, whereas P. canadensis rarely if ever encounters these compounds. Analysis of their furanocoumarin-metabolic profiles indicates that these species induce cytochrome P450 monooxygenases (P450s) capable of metabolizing linear and angular furanocoumarins to varying degrees in response to dietary supplementation with xanthotoxin (a linear furanocoumarin). In P. glaucus, metabolism is induced to a significantly higher level than in P. canadensis. Cloning of multiple P450 genes from each species has revealed that both species contain and express two groups of P450s, designated CYP6B4 and CYP6B17, that are related to the P. glaucus CYP6B4v1 enzyme known to metabolize an array of furanocoumarins. Expression patterns of the CYP6B4 and CYP6B17 group transcripts differ in these species in both their basal and furanocoumarin-inducible levels. In P. glaucus, CYP6B4 transcripts, which are not detectable constitutively, are 311-fold induced by xanthotoxin and CYP6B17 transcripts, which are detectable constitutively, are 3-fold induced by xanthotoxin. In P. canadensis, CYP6B4 transcripts are only 8-fold induced and CYP6B17 transcripts are 13-fold induced. These findings are consistent with the postulated evolutionary history of these two species, according to which P. glaucus maintains its association with rutaceous hostplants and P. canadensis has differentiated to utilize hostplants in other families more extensively.

Expression of CYP6B1 and CYP6B3 cytochrome P450 monooxygenases and furanocoumarin metabolism in different tissues of Papilio polyxenes (Lepidoptera: Papilionidae)

The CYP6B1 and CYP6B3 cytochrome P450 monooxygenases in the midgut of the black swallowtail participate in the metabolism of toxic furanocoumarins present in its host plants. In this study, biochemical analyses indicate that the fat body metabolizes significant amounts of the linear furanocoumarins bergapten and xanthotoxin after larvae feed on xanthotoxin. Northern analyses of the combined CYP6B1/3 transcript expression patterns indicate that transcripts in this P450 subfamily are induced in the midgut and fat body by xanthotoxin. Semi-quantitative RT-PCR analyses of individual CYP6B1/CYP6B3 mRNAs indicate that CYP6B1 transcripts are induced by xanthotoxin in all tissues examined and that CYP6B3 transcripts are induced in the fat body only. These results indicate that the fat body participates in the P450-mediated metabolism of excess furanocoumarins unmetabolized by the midgut. Although transcripts of both genes were detected and CYP6B1 transcripts were induced by xanthotoxin in the integument, furanocoumarin metabolism was not detected. Comparison of these P450 promoters with the promoters of alcohol dehydrogenase genes expressed in the fat bodies of several Drosophila species suggest that the xanthotoxin inducibilities of these P450 genes in fat bodies are regulated by elements other than those modulating expression of Adh genes.

Role of proteolysis in determining potency of Bacillus thuringiensis Cry1Ac delta-endotoxin

Bacillus thuringiensis protein delta-endotoxins are toxic to a variety of different insect species. Larvicidal potency depends on the completion of a number of steps in the mode of action of the toxin. Here, we investigated the role of proteolytic processing in determining the potency of the B. thuringiensis Cry1Ac delta-endotoxin towards Pieris brassicae (family: Pieridae) and Mamestra brassicae (family: Noctuidae). In bioassays, Cry1Ac was over 2,000 times more active against P. brassicae than against M. brassicae larvae. Using gut juice purified from both insects, we processed Cry1Ac to soluble forms that had the same N terminus and the same apparent molecular weight. However, extended proteolysis of Cry1Ac in vitro with proteases from both insects resulted in the formation of an insoluble aggregate. With proteases from P. brassicae, the Cry1Ac-susceptible insect, Cry1Ac was processed to an insoluble product with a molecular mass of approximately 56 kDa, whereas proteases from M. brassicae, the non-susceptible insect, generated products with molecular masses of approximately 58, approximately 40, and approximately 20 kDa. N-terminal sequencing of the insoluble products revealed that both insects cleaved Cry1Ac within domain I, but M. brassicae proteases also cleaved the toxin at Arg423 in domain II. A similar pattern of processing was observed in vivo. When Arg423 was replaced with Gln or Ser, the resulting mutant toxins resisted degradation by M. brassicae proteases. However, this mutation had little effect on toxicity to M. brassicae. Differential processing of membrane-bound Cry1Ac was also observed in qualitative binding experiments performed with brush border membrane vesicles from the two insects and in midguts isolated from toxin-treated insects.

Molecular phylogeny of swallowtail butterflies of the tribe Papilionini (Papilionidae, Lepidoptera)

Swallowtail butterflies of the tribe Papilionini number about 225 species and are currently used as model organisms in several research areas, including genetics, chemical ecology and phylogenetics of host plant utilization and mimicry, mechanisms of speciation, and conservation. We have inferred phylogenetic relationships for a sample of 18 species of the genus Papilio (sensu lato) and five outgroup taxa by sequencing two stretches of mitochondrial DNA that correspond to segments 12886-13370 and 12083-12545 of Drosophila melanogaster mitochondrial DNA and consist of sections of the genes for the large ribosomal RNA and subunit 1 of NADH-dehydrogenase. Our data support the monophyly of Papilio and, within it, of several traditionally recognized subgroups. Species belonging to groups that utilize primarily Rutaceae as larval foodplants form two clusters, corresponding to Old World and American taxa, respectively, while two previously recognized clades-of American and South Asian-Austronesian origin-whose members were known to feed mostly on Lauraceae and Magnoliaceae, are observed to form a clade. The sister group of Papilio is found to be the South Asian genus Meandrusa, which also happens to feed on Lauraceae. The latter plant family is therefore the probable larval host of the ancestor Papilio and the shift to Rutaceae (which four-fifths of extant Papilio species use as foodplants) is more likely to have occurred only after the initial diversification of the genus.