RNAi efficiency is enhanced through knockdown of double-stranded RNA-degrading enzymes in butterfly Papilio xuthus.

The efficiency of RNA interference (RNAi) has always limited the research on the phenotype innovation of Lepidoptera insects. Previous studies have found that double-stranded RNA-degrading enzyme (dsRNase) is an important factor in RNAi efficiency, but there have been no relevant reports in butterflies (Papilionoidea). Papilio xuthus is one of the important models in butterflies with an extensive experimental application value. To explore the effect of dsRNase in the RNAi efficiency on butterflies, six dsRNase genes (PxdsRNase 1-6) were identified in P. xuthus genome, and their dsRNA-degrading activities were subsequently detected by ex vivo assays. The result shows that the dsRNA-degrading ability of gut content (<1 h) was higher than hemolymph content (>12 h). We then investigated the expression patterns of these PxdsRNase genes during different tissues and developmental stages, and related RNAi experiments were carried out. Our results show that different PxdsRNase genes had different expression levels at different developmental stages and tissues. The expression of PxdsRNase2, PxdsRNase3, and PxdsRNase6 were upregulated significantly through dsGFP injection, and PxdsRNase genes can be silenced effectively by injecting their corresponding dsRNA. RNAi-of-RNAi studies with PxEbony, which acts as a reporter gene, observed that silencing PxdsRNase genes can increase RNAi efficiency significantly. These results confirm that silencing dsRNase genes can improve RNAi efficiency in P. xuthus significantly, providing a reference for the functional study of insects such as butterflies with low RNAi efficiency.

Genome-wide identification and gene-editing of pigment transporter genes in the swallowtail butterfly Papilio xuthus.

BACKGROUND: Insect body coloration often functions as camouflage to survive from predators or mate selection. Transportation of pigment precursors or related metabolites from cytoplasm to subcellular pigment granules is one of the key steps in insect pigmentation and usually executed via such transporter proteins as the ATP-binding cassette (ABC) transmembrane transporters and small G-proteins (e.g. Rab protein). However, little is known about the copy numbers of pigment transporter genes in the butterfly genomes and about the roles of pigment transporters in the development of swallowtail butterflies. RESULTS: Here, we have identified 56 ABC transporters and 58 Rab members in the genome of swallowtail butterfly Papilio xuthus. This is the first case of genome-wide gene copy number identification of ABC transporters in swallowtail butterflies and Rab family in lepidopteran insects. Aiming to investigate the contribution of the five genes which are orthologous to well-studied pigment transporters (ABCG: white, scarlet, brown and ok; Rab: lightoid) of fruit fly or silkworm during the development of swallowtail butterflies, we performed CRISPR/Cas9 gene-editing of these genes using P. xuthus as a model and sequenced the transcriptomes of their morphological mutants. Our results indicate that the disruption of each gene produced mutated phenotypes in the colors of larvae (cuticle, testis) and/or adult eyes in G0 individuals but have no effect on wing color. The transcriptomic data demonstrated that mutations induced by CRISPR/Cas9 can lead to the accumulation of abnormal transcripts and the decrease or dosage compensation of normal transcripts at gene expression level. Comparative transcriptomes revealed 606 ~ 772 differentially expressed genes (DEGs) in the mutants of four ABCG transporters and 1443 DEGs in the mutants of lightoid. GO and KEGG enrichment analysis showed that DEGs in ABCG transporter mutants enriched to the oxidoreductase activity, heme binding, iron ion binding process possibly related to the color display, and DEGs in lightoid mutants are enriched in glycoprotein binding and protein kinases. CONCLUSIONS: Our data indicated these transporter proteins play an important role in body color of P. xuthus. Our study provides new insights into the function of ABC transporters and small G-proteins in the morphological development of butterflies.

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.

Identification and functional analysis of an iron-binding protein, ferritin heavy chain subunit, from the swallowtail butterfly, Papilio xuthus.

Ferritin, which is ubiquitous among all living organisms, plays a crucial role in maintaining iron homeostasis, immune response, and detoxification. In the present research, we identified an iron-binding protein, ferritin heavy chain subunit, from Papilio xuthus and named PxFerHCH. The complete complementary DNA of PxFerHCH was 1,252 bp encoding a sequence of 211 amino acids, which includes an iron-responsive element. Phylogenetic analysis showed that PxFerHCH is clustered with Manduca sexta and Galleria mellonella ferritin heavy chain subunits. Expression levels of PxFerHCH in various tissues were analyzed by reverse transcription quantitative polymerase chain reaction, and the results exhibited that PxFerHCH was expressed in all tissues with the highest expression in the fat body. The relative expression level of PxFerHCH in response to bacterial (Escherichia coli and Staphylococcus aureus) challenges sharply increased by about 12 hr postinfection (hpi) and then decreased at 24 hpi. In addition, the iron-binding capacity and antioxidation activity of recombinant PxFerHCH protein were also investigated. These results reveal that PxFerHCH might play an important role in defense against bacterial infection.

Study of the heterologous gene expression characteristics of a new clone of a cell line derived from Papilio xuthus and its serum-free adaptation.

The baculovirus expression vector system using insect cells as a bioreactor has been used for in vitro expression of recombinant proteins and plays an important role in the fields of biology, agronomy, and medicine. Screening suitable host cell lines is an important part of the construction of insect cell baculovirus expression systems. In previous research, we used a single-cell cloning process with the Papilio xuthus cell line RIRI-PX1 and obtained the monoclonal cell line RIRI-PX1-C31. In this study, we compared the basic biological and recombinant protein expression characteristics of RIRI-PX1-C31 and its parent cell line RIRI-PX1 and found that the expression of recombinant ß-galactosidase in RIRI-PX1-C31 was significantly higher than that in the parental cell line. Further serum-free adaptation studies confirmed that RIRI-PX1-C31 can adapt to the growth environment of Express Five Serum-free medium and that its expression level of recombinant ß-galactosidase was significantly higher than that before adaptation.

Connectome of the lamina reveals the circuit for early color processing in the visual pathway of a butterfly.

Connectomics has become a standard neuroscience methodology in a few model animals,1 with the visual system being a popular target of study.2-5 Combining connectomics with circuit and behavioral physiology, recent studies on the color vision of the fruit fly Drosophila melanogaster have focused on the mechanisms underlying early wavelength processing in the optic ganglia.6-8 However, the color vision capabilities of D. melanogaster are limited,9 compared with many flower-visiting insects.10,11 For example, a butterfly Papilio xuthus has six spectral classes of photoreceptors. Each ommatidium contains nine photoreceptors in one of three fixed combinations, making the eye an array of three spectrally distinct ommatidia types.12 Behaviorally, P. xuthus can detect 1 nm differences in light wavelength across the spectrum from ultraviolet to red, outperforming humans.13 What is the neuronal basis of such precise color vision? How does such a system evolve? Addressing these questions requires comparative studies at the circuit level. Here, we performed a connectome analysis in the first optic ganglion, the lamina, of P. xuthus. The lamina comprises cartridges, each typically containing nine photoreceptor axons from a single ommatidium and four second-order neurons. We found abundant inter-photoreceptor connections, which are absent in the lamina of D. melanogaster. We also identified connections between neighboring cartridges, particularly those receiving inputs from spectrally distinct ommatidia. The linear summation of synaptic connections well explains the spectral sensitivity of photoreceptors and second-order neurons in the lamina.

Chromatin accessibility profiling provides insights into larval cuticle color and adult longevity in butterflies.

Butterflies are diverse in virtually all aspects of their ontogeny, including morphology, life history, and behavior. However, the developmental regulatory mechanisms underlying the important phenotypic traits of butterflies at different developmental stages remain unknown. Here, we investigated the developmental regulatory profiles of butterflies based on transposase accessible chromatin sequencing (ATAC-seq) at three developmental stages in two representative species ( Papilio xuthus and Kallima inachus). Results indicated that 15%-47% of open chromatin peaks appeared in associated genes located 3 kb upstream (i.e., promoter region) of their transcription start site (TSS). Comparative analysis of the different developmental stages indicated that chromatin accessibility is a dynamic process and associated genes with differentially accessible (DA) peaks show functions corresponding to their phenotypic traits. Interestingly, the black color pattern in P. xuthus 4th instar larvae may be attributed to promoter peak-related genes involved in the melanogenesis pathway. Furthermore, many longevity genes in 5th instar larvae and pupae showed open peaks 3 kb upstream of their TSS, which may contribute to the overwintering diapause observed in K. inachus adults. Combined with RNA-seq analysis, our data demonstrated that several genes enriched in the melanogenesis and longevity pathways also exhibit higher expression, confirming that the expression of genes may be closely related to their phenotypic traits. This study offers new insights into larval cuticle color and adult longevity in butterflies and provides a resource for investigating the developmental regulatory mechanisms underlying butterfly ontogeny.

Coexpression of Three Odorant-Binding Protein Genes in the Foreleg Gustatory Sensilla of Swallowtail Butterfly Visualized by Multicolor FISH Analysis.

Lepidopteran insects are mostly monophagous or oligophagous. Female butterflies distinguish their host plants by detecting a combination of specific phytochemicals through the gustatory sensilla densely distributed on their foreleg tarsi, thereby ensuring oviposition on appropriate host plants. In this study, to gain insight into the molecular mechanism underlying host plant recognition by the gustatory sensilla, using Asian swallowtail, Papilio xuthus, we focused on a family of small soluble ligand-binding molecules, odorant-binding proteins (OBPs), and found that three OBP genes showed enriched expression in the foreleg tarsus. Multicolor fluorescence in situ hybridization analyses demonstrated the coexpression of these three OBP genes at the bases of the foreleg gustatory sensilla. Further analyses on other appendages revealed that PxutOBP3 was exclusively expressed in the tissues which could have direct contact with the leaf surface, suggesting that this OBP gene specifically plays an important role in phytochemicals perception.

Integrated Analysis of Transcriptome and Proteome to Reveal Pupal Color Switch in Papilio xuthus Butterflies.

Pupal color polyphenism in Papilio butterflies, including green, intermediate, or brown, is an excellent study system for understanding phenotypic plasticity. Previous studies suggested that development of brown pupae may be controlled by a hormone called pupal-cuticle-melanizing-hormone (PCMH) which is synthesized and secreted from brain-suboesophageal ganglion and prothoracic ganglion complexes (Br-SG-TG1) during the pre-pupa stage. However, detailed molecular mechanisms of neuroendocrine regulation in pupal color development remain unknown. In this study, we integrated the expression profiles of transcriptome and proteome at pre-pupa stages [2 h after gut purge (T1) and 3 h after forming the garter around the body (T2)] and pigmentation stages [10 h after ecdysis (T3) and 24 h after ecdysis (T4)] to identify important genes and pathways underlying the development of green and brown pupa in the swallowtail butterfly Papilio xuthus. Combined comparisons of each developmental stage and each tissue under green and brown conditions, a total of 1042 differentially expressed genes (DEGs) and 430 different abundance proteins (DAPs) were identified. Weighted gene co-expression network analysis (WGCNA) and enrichment analysis indicate that these DEGs were mainly related to oxidation-reduction, structural constituent of cuticle, and pigment binding. Soft clustering by Mfuzz and enrichment analysis indicate that these DAPs are mainly involved in tyrosine metabolism, insect hormone biosynthesis, and melanogenesis. By homologous alignment, we further identified those genes encoding neuropeptides (51), GPCRs (116), G-proteins (8), cuticular proteins (226), chitinases (16), and chitin deacetylases (8) in the whole genome of P. xuthus and analyzed their expression profiles. Although we identified no gene satisfying with hypothesized expression profile of PCMH, we found some genes in the neuropeptide cascade showed differentially expressed under two pupal color conditions. We also found that Toll signaling pathway genes, juvenile hormone (JH) related genes, and multiple cuticular proteins play important roles in the formation of selective pupal colors during the prepupal-pupal transition. Our data also suggest that both green and brown pupa include complex pigment system that is regulated by genes involved in black, blue, and yellow pigments. Our results provide important insights into the evolution of pupal protective colors among swallowtail butterflies.

Deciphering Novel Antimicrobial Peptides from the Transcriptome of Papilio xuthus.

An insect's innate immune system is the front line of defense against many invading microorganisms. One of the important components of this defense system is antimicrobial peptides (AMPs). Papiliocin is a well-studied antimicrobial peptide (AMP) isolated from the swallowtail butterfly, Papilio xuthus, and it was previously reported to be effective against Gram-positive bacteria, Gram-negative bacteria, and fungi, particularly in drug resistant Gram-negative bacteria. Hence, we aimed to identify novel AMPs from Papilio xuthus using its transcriptome. We immunized the swallowtail butterfly with Escherichia coli, Staphylococcus aureus, Candida albicans, and the total RNA was isolated. De novo transcriptome assembly and functional annotations were conducted, and AMPs were predicted using an in-silico pipeline. The obtained 344,804,442 raw reads were then pre-processed to retrieve 312,509,806 (90.6%) total clean reads. A total of 38,272 unigenes were assembled with the average length of 1010 bp. Differential gene expression analysis identified 584 and 1409 upregulated and downregulated genes, respectively. The physicochemical, aggregation, and allergen propensity were used as filtration criteria. A total of 248 peptides were predicted using our in-house pipeline and the known AMPs were removed, resulting in 193 novel peptides. Finally, seven peptides were tested in vitro and three peptides (Px 5, 6, and 7) showed stronger antimicrobial activity against Gram-negative bacteria and yeast. All the tested peptides were non-allergens. The identified novel AMPs may serve as potential candidates for future antimicrobial studies.

Complete mitochondrial genome and phylogenetic analysis of the Asian swallowtail, Papilio xuthus.

The Asian swallowtail, Papilio xuthus, a member of the Papilioniniis family, is one of the most popular swallowtail. Here, we first assembled the complete mitochondrial genome of P. xuthus. The mitogenome is 15 359 bp in length, including 13 protein-coding genes (PCGs), 22 transfer RNA genes, and 2 ribosomal RNA genes. The structure of non-coding control region was also analyzed. The gene content and the order are in accord with the common Invertebrate form. Phylogenetic tree was constructed based on the complete mitogenomes of P. xuthus and closely related nine swallowtail species to assess their phylogenic relationship and evolution. The complete mitochondrial genome data would be useful for further studies of P. xuthus.

The complete mitochondrial genome of Papilio xuthus (Lepidoptera: Papilionidae) and implications for Papilioninae taxonomy.

The complete mitochondrial genome of Papilio xuthus (Lepidoptera: Papilionidae) was determined in this study. The mitochondrial genome is a circular molecule of 15 359 and contains 37 genes including 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes and one control region. The nucleotide composition of the A. chinensis mitogenome is strongly biased toward A + T nucleotides (80.45%). Nine protein-coding genes and 14 tRNA genes are encoded on the H strand, and the other four protein-coding genes and eight tRNA genes are encoded on the L strand. The gene order and the orientation of their mitogenomes were similar to all know Papilionidae species. Finally, the phylogenetic relationships of 11 Papilionidae species were reconstructed based on complete mitochondrial genome using the Bayesian inference (BI) and the maximum-likelihood (ML) method. These molecular-based phylogenies support the traditional morphologically based view of relationships within the Papilionidae.