Chromosomal-level reference genome of the moth Heortia vitessoides (Lepidoptera: Crambidae), a major pest of agarwood-producing trees.

The moth Heortia vitessoides Moore (Lepidoptera: Crambidae) is a major pest of ecologically, commercially and culturally important agarwood-producing trees in the genus Aquilaria. In particular, H. vitessoides is one of the most destructive defoliating pests of the incense tree Aquilaria sinesis, which produces a valuable fragrant wood used as incense and in traditional Chinese medicine [33]. Nevertheless, a genomic resource for H. vitessoides is lacking. Here, we present a chromosomal-level assembly for H. vitessoides, consisting of a 517 megabase (Mb) genome assembly with high physical contiguity (scaffold N50 of 18.2 Mb) and high completeness (97.9% complete BUSCO score). To aid gene annotation, 8 messenger RNA transcriptomes from different developmental stages were generated, and a total of 16,421 gene models were predicted. Expansion of gene families involved in xenobiotic metabolism and development were detected, including duplications of cytosolic sulfotransferase (SULT) genes shared among lepidopterans. In addition, small RNA sequencing of 5 developmental stages of H. vitessoides facilitated the identification of 85 lepidopteran conserved microRNAs, 94 lineage-specific microRNAs, as well as several microRNA clusters. A large proportion of the H. vitessoides genome consists of repeats, with a 29.12% total genomic contribution from transposable elements, of which long interspersed nuclear elements (LINEs) are the dominant component (17.41%). A sharp decrease in the genome-wide percentage of LINEs with lower levels of genetic distance to family consensus sequences suggests that LINE activity has peaked in H. vitessoides. In contrast, opposing patterns suggest a substantial recent increase in DNA and LTR element activity. Together with annotations of essential sesquiterpenoid hormonal pathways, neuropeptides, microRNAs and transposable elements, the high-quality genomic and transcriptomic resources we provide for the economically important moth H. vitessoides provide a platform for the development of genomic approaches to pest management, and contribute to addressing fundamental research questions in Lepidoptera.

Biological composition analysis of a natural medicine, Faeces Vespertilionis, with complex sources using DNA metabarcoding.

Faeces Vespertilionis is a commonly used fecal traditional Chinese medicine. Traditionally, it is identified relying only on morphological characters. This poses a serious challenge to the composition analysis accuracy of this complex biological mixture. Thus, for quality control purposes, an accurate and effective method should be provided for taxonomic identification of Faeces Vespertilionis. In this study, 26 samples of Faeces Vespertilionis from ten provinces in China were tested using DNA metabarcoding. Seven operational taxonomic units (OTUs) were detected as belonging to bats. Among them, Hipposideros armiger (Hodgson, 1835) and Rhinolophus ferrumequinum (Schober and Grimmberger, 1997) were the main host sources of Faeces Vespertilionis samples, with average relative abundances of 59.3% and 24.1%, respectively. Biodiversity analysis showed that Diptera and Lepidoptera were the most frequently consumed insects. At the species level, 19 taxa were clearly identified. Overall, our study used DNA metabarcoding to analyze the biological composition of Faeces Vespertilionis, which provides a new idea for the quality control of this special traditional Chinese medicine.

Complete mitochondrial genome of two Thitarodes species (Lepidoptera, Hepialidae), the host moths of Ophiocordyceps sinensis and phylogenetic implications.

Thitarodes (Lepidoptera, Hepialidae) is the only genus that hosts to the Ophiocordyceps sinensis, a traditional Chinese medicine considered as a powerful medicinal supplement. In this study, the complete mitochondrial genomes (mitogenomes) of two species, T. damxungensis and T. pui, have been sequenced, which are 15,928 bp and 15,362 bp in size respectively, and both contain 13 protein-coding genes (PCGs), 2 rRNAs, 22 tRNAs and an AT-rich region. Like other hepialoids, the gene arrangement of the mitogenomes of T. damxungensis and T. pui is identical to the ancestral arrangement but differs from those of other lepidopteran species on account of the different arrangements of trnM, trnI, and trnQ. The size of AT-rich region is 545 bp in T. damxungensis and 1030 bp in T. pui. Tandem repetition in the AT-rich region is responsible for the length difference of the A + T-rich region in both species. In Hepialidae, the phylogenetic study based on the dataset of the sequences that combined the protein-coding genes and RNA genes suggested that the species T. yunnanensis should still belong to the genus Thitarodes rather than Ahamns, which is different from the results based on the traditional phylogeny.

Genome-wide Identification of Tebufenozide Resistant Genes in the smaller tea tortrix, Adoxophyes honmai (Lepidoptera: Tortricidae).

The smaller tea tortrix, Adoxophyes honmai, has developed strong resistance to tebufenozide, a diacylhydrazine-type (DAH) insecticide. Here, we investigated its mechanism by identifying genes responsible for the tebufenozide resistance using various next generation sequencing techniques. First, double-digest restriction site-associated DNA sequencing (ddRAD-seq) identified two candidate loci. Then, synteny analyses using A. honmai draft genome sequences revealed that one locus contained the ecdysone receptor gene (EcR) and the other multiple CYP9A subfamily P450 genes. RNA-seq and direct sequencing of EcR cDNAs found a single nucleotide polymorphism (SNP), which was tightly linked to tebufenozide resistance and generated an amino acid substitution in the ligand-binding domain. The binding affinity to tebufenozide was about 4 times lower in in vitro translated EcR of the resistant strain than in the susceptible strain. RNA-seq analyses identified commonly up-regulated genes in resistant strains, including CYP9A and choline/carboxylesterase (CCE) genes. RT-qPCR analysis and bioassays showed that the expression levels of several CYP9A and CCE genes were moderately correlated with tebufenozide resistance. Collectively, these results suggest that the reduced binding affinity of EcR is the main factor and the enhanced detoxification activity by some CYP9As and CCEs plays a supplementary role in tebufenozide resistance in A. honmai.

Silence of ryanodine receptor gene decreases susceptibility to chlorantraniliprole in the oriental armyworm, Mythimna separata Walker.

The ryanodine receptors of insects are the main target sites of diamide insecticides, which show highly selective insecticidal activity relative to toxicity in mammals and provide a novel option for managing lepidopteran pests. The oriental armyworm, Mythimna separata (Walker), is a destructive pest of agricultural crops, and great efforts have been undertaken to control this pest including repeated insecticide applications. In this study, full-length cDNA of a ryanodine receptor gene from M. separata (MsRyR) was cloned and characterized. The cDNA of MsRyR had a 15,372 bp open reading frame and encoded 5124 amino acids (GenBank ID: MG712298). MsRyR shares 78-97% identity with RyR isoforms of other insects, and <50% identities with Homo sapiens RyRs 1-3. Temporal and spatial expression analysis detected MsRyR at all developmental stages and in all tissues. The highest relative levels of MsRyR were detected in the second instar and head. Exposure to chlorantraniliprole after 24 h significantly increased the expression levels of whole body MsRyR mRNA. In addition, dietary ingestion of dsMsRyR significantly reduced the mRNA level of MsRyR and greatly decreased chlorantraniliprole-induced mortality. Our results revealed that the MsRyR could be the molecular target of chlorantraniliprole, and provided the basis for further understanding the resistance mechanism of chlorantraniliprole.

Molecular characterization and expression profiling of ryanodine receptor gene in the pink stem borer, Sesamia inferens (Walker).

The susceptibilities of three field populations of pink stem borer (PSB), Sesamia inferens (walker) to diamide insecticides, chlorantraniliprole and flubendiamide, were evaluated in this study. The results showed that these PSB field populations were still sensitive to the two diamide insecticides after many years of exposure. To further understand PSB and diamide insecticide, the full-length ryanodine receptor (RyR) cDNA (named as SiRyR), the molecular target of diamide insecticides was cloned from PSB and characterized. The SiRyR gene contains an open reading frame of 15,420 nucleotides, encoding 5140 amino acid residues, which shares 77% to 98% sequence identity with RyR homologous of other insects. All hallmarks of RyR proteins are conserved in the SiRyR protein, including the conserved C-terminal domain with the consensus calcium-biding EF-hands (calcium-binding motif), the six transmembrane domains, as well as mannosyltransferase, IP3R and RyR (pfam02815) (MIR) domains. Real-time qPCR analysis revealed that the highest mRNA expression levels of SiRyR were observed in pupa and adults, especially in males. SiRyR was expressed at the highest level in thorax, and the lowest level in wing. The full genetic characterization of SiRyR could provide useful information for future functional expression studies and for discovery of new insecticides with selective insecticidal activity.

Larval transcriptomic response to host plants in two related phytophagous lepidopteran species: implications for host specialization and species divergence.

BACKGROUND: Most phytophagous insects have morphological, behavioral and physiological adaptations allowing them to specialize on one or a few plant species. Identifying the mechanisms involved in host plant specialization is crucial to understand the role of divergent selection between different environments in species diversification, and to identify sustainable targets for the management of insect pest species. In the present study, we measured larval phenotypic and transcriptomic responses to host plants in two related phytophagous lepidopteran species: the European corn borer (ECB), a worldwide pest of maize, and the adzuki bean borer (ABB), which feeds of various dicotyledons. Our aim was to identify the genes and functions underlying host specialization and/or divergence between ECB and ABB. RESULTS: At the phenotypic level, we observed contrasted patterns of survival, weight gain and developmental time between ECB and ABB, and within ECB and ABB reared on two different host plants. At the transcriptomic level, around 8% of the genes were differentially expressed (DE) between species and/or host plant. 70% of these DE genes displayed a divergent pattern of expression between ECB and ABB, regardless of the host, while the remaining 30% were involved in the plastic response between hosts. We further categorized plastic DE genes according to their parallel or opposite pattern between ECB and ABB to specifically identify candidate genes involved in the species divergence by host specialization. These candidates highlighted a comprehensive response, involving functions related to plant recognition, digestion, detoxification, immunity and development. Last, we detected viral, bacterial, and yeast genes whose incidence contrasted ECB and ABB samples, and maize and mugwort conditions. We suggest that these microorganism communities might influence the survival, metabolism and defense patterns observed in ECB and ABB larvae. CONCLUSIONS: The comprehensive approach developed in the present study allowed to identify phenotypic specialization patterns and underlying candidate molecular mechanisms, and highlighted the putative role of microorganisms in the insect-host plant interaction. These findings offer the opportunity to pinpoint specific and sustainable molecular or physiological targets for the regulation of ECB pest populations.

International Barcode of Life: Focus on big biodiversity in South Africa.

Participants in the 7th International Barcode of Life Conference (Kruger National Park, South Africa, 20-24 November 2017) share the latest findings in DNA barcoding research and its increasingly diversified applications. Here, we review prevailing trends synthesized from among 429 invited and contributed abstracts, which are collated in this open-access special issue of Genome. Hosted for the first time on the African continent, the 7th Conference places special emphasis on the evolutionary origins, biogeography, and conservation of African flora and fauna. Within Africa and elsewhere, DNA barcoding and related techniques are being increasingly used for wildlife forensics and for the validation of commercial products, such as medicinal plants and seafood species. A striking trend of the conference is the dramatic rise of studies on environmental DNA (eDNA) and on diverse uses of high-throughput sequencing techniques. Emerging techniques in these areas are opening new avenues for environmental biomonitoring, managing species-at-risk and invasive species, and revealing species interaction networks in unprecedented detail. Contributors call for the development of validated community standards for high-throughput sequence data generation and analysis, to enable the full potential of these methods to be realized for understanding and managing biodiversity on a global scale.

Ryanodine receptor genes of the rice stem borer, Chilo suppressalis: Molecular cloning, alternative splicing and expression profiling.

The ryanodine receptor (RyR) of the calcium release channel is the main target of anthranilic and phthalic diamide insecticides which have high selective insecticidal activity relative to mammalian toxicity. In this study, the full-length cDNA of Chilo suppressalis RyR (CsRyR) was isolated and characterized. The CsRyR mRNA has an open reading frame (ORF) of 15,387bp nucleotides, which encodes 5128 amino acids with GenBank ID: KR088972. Comparison of protein sequences showed that CsRyR shared high identities with other insects of 77-96% and lower identity to mammals and nematodes with only 42-45%. One alternative splicing site (KENLG) unique to Lepidoptera was found and two exclusive exons of CsRyR (I /II) were revealed. Spatial and temporal expression of CsRyR mRNA was at the highest relative level in 3rd instar larvae and head (including brain and muscle), and at the lowest expression level in egg and fat body. The expression levels of whole body CsRyR mRNA were increased remarkably after injection of 4th instar larvae with chlorantraniliprole at 0.004 to 0.4µg/g. This structural and functional information on CsRyR provides the basis for further understanding the selective action of chlorantraniliprole and possibly other diamide insecticides.

Plant phenolics are detoxified by prophenoloxidase in the insect gut.

Plant phenolics are a group of important secondary metabolites that are toxic to many animals and insects if ingested at high concentrations. Because most insects consume plant phenolics daily, they have likely evolved the capacity to detoxify these compounds. Here, we used Drosophila melanogaster, Bombyx mori and Helicoverpa armigera as models to study the metabolism of plant phenolics by prophenoloxidases. We found that insect foreguts release prophenoloxidases into the lumen, and that the survival of prophenoloxidase-deletion mutants was impaired when fed several plant phenolics and tea extracts. Using l-DOPA as a model substrate, biochemical assays in large Lepidopteran insects demonstrated that low levels of l-DOPA are rapidly metabolized into intermediates by phenoloxidases. Feeding with excess l-DOPA showed that the metabolic intermediate 5,6-dihydroxyindole reached the hindgut either by passing directly through the midgut, or by transport through the hemolymph. In the hindgut, 5,6-dihydroxyindole was further oxidized by prophenoloxidases. Intermediates exerted no toxicity in the hemocoel or midgut. These results show that plant phenolics are not toxic to insects unless prophenoloxidase genes are lost or the levels of phenolics exceed the catalytic activity of the gut prophenoloxidases.

Identification and molecular characterization of a chitin-binding protein from the beet webworm, Loxostege sticticalis L.

As the first crucial barrier in the midgut of insects, the peritrophic membrane (PM) plays an important role in preventing external invasion. PM proteins, as the major components of the PM, determine the structure and function of this membrane. A new PM protein, named LstiCBP, from the PM of Loxostege sticticalis larvae was identified using cDNA library screening. The full cDNA of LstiCBP is 2606 bp in length and contains a 2403 bp ORF that encodes an 808-amino acid preprotein with a 15-amino acid as signal peptide. The deduced protein sequence of the cDNA contains 8 cysteine-rich chitin-binding domains (CBDs). Recombinant LstiCBP was successfully expressed in BL21 cells using recombinant plasmid DNA and showed high chitin-binding activity. LstiCBP expression was detected in the midgut at both the transcriptional and translational levels; however, the biochemical and physiological functions of LstiCBP in L. sticticalis require further investigation.

Detecting deep divergence in seventeen populations of tea geometrid (Ectropis obliqua Prout) in China by COI mtDNA and cross-breeding.

The tea geometrid (Ectropis obliqua Prout, Lepidoptera: Geometridae) is a dominant chewing insect endemic in most tea-growing areas in China. Recently some E. obliqua populations have been found to be resistant to the nucleopolyhedrovirus (EoNPV), a host-specific virus that has so far been found only in E. obliqua. Although the resistant populations are morphologically indistinguishable from susceptible populations, we conducted a nationwide collection and examined the genetic divergence in the COI region of the mtDNA in E. obliqua. Phylogenetic analyses of mtDNA in 17 populations revealed two divergent clades with genetic distance greater than 3.7% between clades and less than 0.7% within clades. Therefore, we suggest that E. obliqua falls into two distinct groups. Further inheritance analyses using reciprocal single-pair mating showed an abnormal F1 generation with an unbalanced sex ratio and the inability to produce fertile eggs (or any eggs) through F1 self-crossing. These data revealed a potential cryptic species complex with deep divergence and reproductive isolation within E. obliqua. Uneven distribution of the groups suggests a possible geographic effect on the divergence. Future investigations will be conducted to examine whether EoNPV selection or other factors prompted the evolution of resistance.

Characterization of four midgut aminopeptidase N isozymes from Ostrinia furnacalis strains with different susceptibilities to Bacillus thuringiensis.

The full-length cDNA of four Ofapn aminopeptidases were cloned and sequenced from susceptible and resistant Ostrinia furnacalis strains. Four sequences were identified as APN because they shared the common structural features with APN from Lepidoptera, including the signal peptide, GPI anchor signal, the zinc binding/gluzincin motif HEX2HX18E and the gluzincin aminopeptidase motif GAMEN. Compared with APN sequences from the susceptible strain, there were 9, 5, 10 and 12 amino acid variations in the deduced protein sequences from the resistant strain. There were also differences in mRNA expression of the four Ofapn genes between resistant and susceptible O. furnacalis strains.

Molecular cloning, characterization and expression profiling of a ryanodine receptor gene in Asian corn borer, Ostrinia furnacalis (Guenée).

Ryanodine receptor (RyR) Ca(2+) release channel is the target of diamide insecticides, which show selective insecticidal activity against lepidopterous insects. To study the molecular mechanisms underlying the species-specific action of diamide insecticides, we have cloned and characterized the entire cDNA sequence of RyR from Ostrinia furnacalis (named as OfRyR). The OfRyR mRNA has an Open Reading Frame of 15324 bp nucleotides and encodes a 5108 amino acid polypeptide that displays 79-97% identity with other insects RyR proteins and shows the greatest identity with Cnaphalocrocis medinalis RyR (97%). Quantitative real-time PCR showed that the OfRyR was expressed at the lowest level in egg and the highest level in adult. The relative expression level of OfRyR in first, third and fifth-instar larva were 1.28, 1.19 and 1.99 times of that in egg. Moreover, two alternative splicing sites were identified in the OfRyR gene. One pair of mutually exclusive exons (a/b) were present in the central part of the predicted SPRY domain, and an optional exon (c) was located between the third and fourth RyR domains. Diagnostic PCR demonstrated that exons a and b existed in all developmental stages of OfRyR cDNA, but exon c was not detected in the egg cDNA. And the usage frequencies of these exons showed a significant difference between different developmental stages. These results provided the crucial basis for the functional expression of OfRyR and for the discovery of compound with potentially selective insect activtity.

Protease inhibitor from insect silk-activities of derivatives expressed in vitro and in transgenic potato.

Several recombinant derivatives of serine protease inhibitor called silk protease inhibitor 2 (SPI2), which is a silk component in Galleria mellonella (Lepidoptera, Insecta), were prepared in the expression vector Pichia pastoris. Both the native and the recombinant protease inhibitors were highly active against subtilisin and proteinase K. The synthetic SPI2 gene with Ala codon in the P1 position was fused with mGFP-5 to facilitate detection of the transgene and its protein product. A construct of the fusion gene with plant regulatory elements (promoter 35S and terminator OCS) was inserted into the binary vector pRD400. The final construct was introduced into Agrobacterium tumefaciens that was then used for genetic transformation of the potato variety Velox. The transgene expression was monitored with the aid of ELISA employing polyclonal antibody against natural SPI2. In vitro tests showed increased resistance to the late blight Phytophthora infestans in several transformed lines. No effect was seen on the growth, mortality, life span or reproduction of Spodoptera littoralis (Lepidoptera, Insecta) caterpillars, while feeding on transformed potato plants expressing the fusion protein, indicating that the transformed potatoes may be harmless to non-target organisms.

Influence of gamma irradiation on productivity indices of the edible Emperor moth caterpillar, Cirina forda (Lepidoptera: Saturniidae).

This study was aimed at generating baseline information for sustainable genetic improvement of Cirana forda larvae for entomophagy, through the use of gamma irradiation. Eggs of C. forda were irradiated with increasing doses of gamma rays from 0 to 200 Gy and raised through larval instal stages under laboratory conditions. The Body Weight (BW) and Head Capsule Width (HCW) of the larval instar stages were monitored as indices of productivity. Successful larval emergence was recorded for all irradiation doses tested and BW of the 1st and 2nd instar larvae were not significantly (p > 0.05) different between the control and treated groups (range = 0.021 +/- 0.003 g/larva in the 200 Gy treatment to 0.028 +/- 0.003 g/larva in the control group and 0.105 +/- 0.003 g/larva in 20 Gy treatment to 0.172 +/- 0.009 g/larva in the control group, respectively). On the other hand, BW during the 3rd and 4th larval instars were significantly (p < 0.05) lower among the irradiated treatments than control. Pattern of distribution of HCW was different from that of BW; as HCW increased with irradiation dose from 10-50 Gy during the 3rd and 4th larval instars. Also, HCW during the 5th instar larvae among the irradiated treatments (range = 5.256 +/- 0.012 to 5.662 +/- 0.026 mm) were not higher than that of the 6th instar in the control group (6.065 +/- 0.010 mm). These results suggest promising potentials of the use of gamma irradiation in sustainably improving the productivity of C. forda larvae for entomophagy.

Molecular characterization of a ryanodine receptor gene in the rice leaffolder, Cnaphalocrocis medinalis (Guenée).

Ryanodine receptors (RyRs) are the targets of two novel classes of synthetic insecticidal chemicals, phthalic acid diamides and anthranilic diamides. Isolation of full-length RyR cDNAs is a critical step towards the structural and functional characterization of insect RyRs and an understanding of the molecular mechanisms underlying the species selective toxicity of diamide insecticides. However, there has been little research on the insect RyR genes due to the high molecular weight of the RyR proteins. In this study, we isolated a full-length RyR cDNA (named as CmRyR) from Cnaphalocrocis medinalis, an important rice pest throughout Southeast Asia. The composite CmRyR gene contains an ORF of 15264 bp encoding a protein of 5087 amino acid residues, which shares 79% overall identity with its Drosophila melanogaster homologue. All hallmarks of the RyR proteins are conserved in the CmRyR protein, suggesting that CmRyR is a structural and functional analogue of known RyRs. A multiple sequence alignment illustrates that the insect RyRs share high levels of amino acid sequence identity at the the COOH-terminal region. However, the amino acid residues analogous to the CmRyR residues N(4922), N(4924), N(4935), L(4950), L(4981), N(5013) and T(5064) are unique to lepidopteran RyRs compared with non-lepidopteran insect RyRs. This finding suggests that these residues may be involved in the differences in channel properties between lepidopteran and non-lepidopteran insect RyRs and in the species selective toxicity of diamide insecticides. Furthermore, two alternative splicing sites were identified in the CmRyR gene, one of which was located in the central part of the predicted second SPRY domain. Diagnostic PCR showed that the inclusion frequencies of two mutually exclusive exons (a/b) and one optional exon (c) differed between developmental stages or adult anatomical regions. Our results imply that alternative splicing may be a major means of generating functional diversity in C. medinalis RyR channel.

Prophenoloxidase from Pieris rapae: gene cloning, activity, and transcription in response to venom/calyx fluid from the endoparasitoid wasp Cotesia glomerata.

Prophenoloxidase (PPO) plays an important role in melanization, necessary for defense against intruding parasitoids. Parasitoids have evolved to inject maternal virulence factors into the host hemocoel to suppress hemolymph melanization for the successful development of their progeny. In this study, the full-length complementary DNA (cDNA) of a Pieris rapae PPO was cloned. Its cDNA contained a 2 076-base pair (bp) open reading frame (ORF) encoding 691 amino acids (aa). Two putative copper-binding sites, a proteolytic activation site, three conserved hemocyanin domains, and a thiol ester motif were found in the deduced amino acid sequence. According to both multiple alignment and phylogenetic analysis, P. rapae PPO gene cloned here is a member of the lepidopteran PPO-2 family. Injection of Cotesia glomerata venom or calyx fluid resulted in reduction of P. rapae hemolymph phenoloxidase activity, demonstrating the ability to inhibit the host's melanization. Real-time reverse transcriptase polymerase chain reaction (RT-PCR) showed that transcripts of P. rapae PPO-2 in the haemocytes from larvae had not significantly changed following venom injection, suggesting that the regulation of PPO messenger RNA (mRNA) expression by venom was not employed by C. glomerata to cause failure of melanization in parasitized host. While decreased P. rapae PPO-2 gene expression was observed in the haemocytes after calyx fluid injection, no detectable transcriptional change was induced by parasitization, indicating that transcriptional down-regulation of PPO by calyx fluid might play a minor role involved in inhibiting the host's melanization.

Mutualism favours higher host specificity than does antagonism in plant-herbivore interaction.

Coevolved mutualisms often exhibit high levels of partner specificity. Obligate pollination mutualisms, such as the fig-fig wasp and yucca-yucca moth systems, represent remarkable examples of such highly species-specific associations; however, the evolutionary processes underlying these patterns are poorly understood. The prevailing hypothesis suggests that the high degree of specificity in pollinating seed parasites is the fortuitous result of specialization in their ancestors because these insects are derived from endophytic herbivores that are themselves highly host-specific. Conversely, we show that in the Glochidion-Epicephala obligate pollination mutualism, pollinators are more host-specific than are closely related endophytic leaf-feeding taxa, which co-occur with Epicephala on the same Glochidion hosts. This difference is probably not because of shifts in larval diet (i.e. from leaf- to seed-feeding), because seed-eating lepidopterans other than Epicephala do not show the same degree of host specificity as Epicephala. Species of a tentative sister group of Epicephala each attack several distantly related plants, suggesting that the evolution of strict host specificity is tied to the evolution of pollinator habit. These results suggest that mutualists can attain higher host specificity than that of their parasitic ancestors and that coevolutionary selection can be a strong promoter of extreme reciprocal specialization in mutualisms.

Molecular cloning of apoptosis-inducing Pierisin-like proteins, from two species of white butterfly, Pieris melete and Aporia crataegi.

Pierisin-1, present in cabbage butterfly, Pieris rapae, induces apoptosis against various kinds of cancer cell lines. Another cabbage butterfly, Pieris brassicae, also has an apoptosis-inducing protein, Pierisin-2. These proteins exhibit DNA ADP-ribosylating activity. Pierisin-like proteins are found to be distributed in subtribes Pierina, Aporiina and Appiadina. In this study, we performed the cDNA cloning of Pierisin-like proteins designated Pierisin-3 from gray-veined white, Pieris melete, and Pierisin-4 from black-veined white, Aporia crataegi. The nucleotide sequences of Pierisin-3 and -4 encode an 850 and an 858 amino acid protein, respectively. The partial peptide sequences of Pierisin-3 and -4 purified from pupae were identical to the deduced amino acid sequence of ORF. The deduced amino acid sequence revealed that Pierisin-3 is 93% similar to Pierisin-1 and Pierisin-4 is 64%. Pierisin-3 and -4 synthesized in vitro with the rabbit reticulocyte lysate exhibited apoptosis-inducing activity against human cervical carcinoma HeLa and human gastric carcinoma TMK-1 cells. Site-directed mutagenesis at a glutamic acid residue comprising the NAD-binding site resulted in a significant decrease in cytotoxicity of both proteins. Moreover, the proteins incubated with calf thymus DNA and beta-NAD resulted in the formation of N(2)-(ADP-ribos-1-yl)-2'-deoxyguanosine, as in the case of Pierisin-1 and -2. These findings could provide useful information for understanding the importance of apoptosis-inducing ability and molecular evolution of Pierisin-like proteins in family Pieridae.