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.

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.

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.

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.

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.

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.
蝴蝶的多样性体现在其个体发育的各个方面，如形态、生活史、行为等。然而，蝴蝶的许多重要表型性状的发育调控机制尚不清楚。该研究运用转座酶研究染色质可及性 (Assay of transposase accessible chromatin sequencing, ATAC-seq) 技术对两种蝴蝶 (柑橘凤蝶 Papilio xuthus和枯叶蛱蝶 Kallima inachus) 的三个发育时期 (4龄、5龄和蛹期) 的染色质可及性区域进行了研究。在柑橘凤蝶三个发育时期分别鉴定到20 367 (4龄)、9 648 (5龄) 和12 001 (蛹期) 个染色质可及性区域；在枯叶蛱蝶的三个发育时期分别鉴定到28 751 (4龄)、85 433 (5龄) 和105 067 (蛹期) 个染色质可及性区域；有15%–47%的染色质可及性区域位于相关基因转录起始位点 (TSS) 上游3kb (即启动子区域)。另外，对在两个蝴蝶物种三个发育时期中鉴定到的染色质可及性区域进行了Motif预测，再根据已知转录因子 (Transcription Factors, TFs) 数据库搜索匹配，鉴定到了一些与这些Motif结合的TFs。通过对物种内不同发育阶段的比较分析表明，染色质可及性区域是一个动态过程，具有差异可及性区域的相关基因对其表型性状具有相应的功能。我们发现柑橘凤蝶4龄幼虫表皮的黑色模式可能与黑色素生成途径中相关基因启动子区域的染色质可及性区域有关；枯叶蛱蝶5龄幼虫和蛹期的很多长寿基因启动子区域有染色质可及性区域，这些基因可能与其成虫长寿命滞育越冬密切相关。结合RNA-seq分析，这些富集在黑色素生成途径和长寿机制途径的基因均表现出较高的表达水平，进一步证实了这些基因的表达可能与其相应的表型性状密切相关。该研究为蝴蝶幼虫表皮颜色和成虫寿命机制的研究提供了新的思路，也为研究蝴蝶个体发育的调控机制提供了新的数据来源。.

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.

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.

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.

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.