Genome-wide survey of open chromatin regions in two swallowtail butterflies Papilio machaon and P. bianor.

Papilio machaon was assigned as the type species for all butterflies by Linnaeus and P. bianor is a congener but exhibits a great difference in morphology (especially larva and adult color pattern) and larval host plants from P. machaon. Thus, they are the ideal models to investigate genetic mechanisms underlying morphology and plasticity between congeners. The reference genomes of both species were dissected in our previous studies, but little is known about their regulatory genome and the epigenetic regulation of gene expression throughout developmental stages. Here, we profiled the chromatin accessibility and gene expression of three developmental stages (the 4th instar larva [L4], the 5th instar larva [L5], and pupa [P]) using transposase accessible chromatin sequencing (ATAC-seq) and RNA-seq. Results showed that many accessible chromatin peaks were identified at three developmental stages (peak number, P. machaon: 44,977 [L4], 36,919 [L5], 47,147 [P]; P. bianor: 20,341 [L4], 44,668 [L5], 62,249 [P]). Moreover, the number of differentially accessible peaks and differentially expressed genes between larval stages of each butterfly species are significantly fewer than that between larval and pupal stages, suggesting a higher similarity within larvae and a significant difference between larvae and pupae. This study added the annotated information of chromatin accessibility genome-wide of the two papilionid species and will promote the investigation of gene regulation in butterfly evolution.

Genome-wide SNP discovery and genetic diversity evaluation of Liangshan cattle in China.

Liangshan cattle are a very small indigenous breed with adult weight of less than 300 Kg and have been mainly distributed in the Liangshan Yi Autonomous Prefecture of Southwestern Sichuan, China. Due to its long-term adaptation to local environments, Liangshan cattle is a valuable genetic resource and should be paid with more attentions. However, the genetic diversity of Liangshan cattle have not been specifically investigated yet, which would be required when designing the appropriate conservation and utilization programs. In this study, we successfully employed the restriction-site-associated DNA sequencing (RAD-seq) approach to explore a total of 84,854 genome-wide and high-confidence SNPs of Liangshan cattle. All these SNPs were evenly distributed through all chromosomes with an average of 98 SNPs per 1-Mb region. The nucleotide diversity, expected heterozygosity, polymorphism information content of Liangshan cattle were 0.227, 0.223 and 0.183, respectively. Furthermore, there was no obvious difference on the genetic diversity among the three studied geographical populations. In conclusion, we provided a list of SNPs that could be used in the follow-up studies for Liangshan cattle and revealed a relatively high genetic variation in this gene pool.

Chromosome-Level Genome Assembly of Apoderus dimidiatus Voss (Coleoptera: Attelabidae): Insights into Evolution and Behavior.

Attelabidae insects have attracted much attention due to their unique leaf rolling behavior before oviposition. However, the lack of genomic data makes it difficult to understand the molecular mechanism behind their behavior and their evolutionary relationship with other species. To address this gap, we utilized Illumina and Nanopore sequencing platforms along with Hi-C technology to establish a highly accurate whole genome of A. dimidiatus at the chromosome level. The resulting genome size was determined to be 619.26 Mb, with a contig N50 of 50.89 Mb and GC content of 33.89%. Moreover, a total of 12,572 genes were identified, with 82.59% being functionally annotated, and 64.78% designated as repeat sequences. Our subsequent phylogenetic tree analysis revealed that Attelabidae's divergence from Curculionidae occurred approximately 161.52 million years ago. Furthermore, the genome of A. dimidiatus contained 334 expanded gene families and 1718 contracted gene families. In addition, using Phylogenetic Analysis by Maximum Likelihood (PAML), we identified 106 rapidly evolved genes exhibiting significant signals and 540 positively selected genes. Our research endeavors to serve as an invaluable genomic data resource for the study of Attelabidae, offering fresh perspectives for the exploration of its leaf rolling behavior.

The first chromosome-level genome of the stag beetle Dorcus hopei Saunders, 1854 (Coleoptera: Lucanidae).

Stag beetles (Coleoptera: Lucanidae) represent a significant saproxylic assemblage in forest ecosystems and are noted for their enlarged mandibles and male polymorphism. Despite their relevance as ideal models for the study of exaggerated mandibles that aid in attracting mates, the regulatory mechanisms associated with these traits remain understudied, and restricted by the lack of high-quality reference genomes for stag beetles. To address this limitation, we successfully assembled the first chromosome-level genome of a representative species Dorcus hopei. The genome was 496.58 Mb in length, with a scaffold N50 size of 54.61 Mb, BUSCO values of 99.8%, and 96.8% of scaffolds anchored to nine pairs of chromosomes. We identified 285.27 Mb (57.45%) of repeat sequences and annotated 11,231 protein-coding genes. This genome will be a valuable resource for further understanding the evolution and ecology of stag beetles, and provides a basis for studying the mechanisms of exaggerated mandibles through comparative analysis.

Xenosyangi sp. nov.: A new twisted-wing parasite species (Strepsiptera, Xenidae) from Gaoligong Mountains, Southwest China.

Here we report a new twisted-wing parasite species of the family Xenidae based on both morphological and molecular evidence. By using nearly complete mitogenomes, we confirmed the twisted-wing parasites on two wasps (Vespavelutina and Vespabicolor) (China: Yunnan) as the same species, and associated its neotenic females and alate males. Combining the mitogenomic data (COI) and morphological traits, this species was identified to be a new species of the genus Xenos, namely Xenosyangi Dong, Liu & Li, sp. nov. Detailed descriptions and illustrations are provided for the new species.

Complete mitogenome and phylogenetic significance of Metoecus javanus (Pic, 1913) (Coleoptera: Ripiphoridae) from Southwest China, with notes on morphological traits of adult and immature stages.

Wedge-shaped beetles (Ripiphoridae) not only exhibit enigmatic morphological and biological traits but also disputable phylogenetic positions. At present, however, genetic information regarding this family remains limited. In this study, we report on the complete mitogenome of one ripiphorid beetle, Metoecus javanus (Pic, 1913), from Southwest China, as well as its different developmental stages, populations, and morphological variability. The complete mitogenome of M. javanus was 16 109 bp in length, containing 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and a noncoding control region. Of the 37 genes, 23 were located on the majority strand (J-strand) and 14 were located on the minority strand (N-strand). All PCGs started with "ATN" (N represents A, T, G, and C), and terminated with "TAA", except for NAD1 with "TAG" and COX2 with a single "T". The five most used codons in the PCGs were UUA(L), UCU(S2), CCU(P), UCA(S2), and GGA(G), indicating a strong bias toward A + T-rich codons. All 22 tRNAs showed typical cloverleaf structures, except trnS1, which lacked a dihydrouridine (DHU) stem. The control region contained five types of tandem repeats, with the repeat units ranging from 17 to 24 bp. Phylogenetic analysis of the concatenated set of 13 PCGs and two rRNAs (rrnL and rrnS) of M. javanus and 17 other Tenebrionoidea species indicated that M. javanus did not cluster with Pelecotoma fennica (Ripiphoridae: Pelecotominae), another wedge-shaped beetle, but was located at the base of the Mordellidae + P. fennica clade. This reconstruction supported the paraphyly of Ripiphoridae with respect to Mordellidae. Using the mitogenome COX1 data, wedge-shaped beetles from different stages (male adult, female adult, and pupa), different geographical populations (Nujiang and Lincang), and different wasp hosts (Vespidae: Vespa velutina Lepeletier, 1836 and Vespa bicolor Fabricius, 1787) were identified as a same species (i.e., M. javanus). Based on morphological observations of all specimens, we identified and described variability in the adult pronotum, male genitalia, and pupa of M. javanus. The present results provide important genetic and morphological information for further investigations on the phylogenetic position of Ripiphoridae and its evolutionary diversity.

High-quality reference genomes of swallowtail butterflies provide insights into their coloration evolution.

Swallowtail butterflies (Papilionidae) are a historically significant butterfly group due to their colorful wing patterns, extensive morphological diversity, and phylogenetically important position as a sister group to all other butterflies and have been widely studied regarding ecological adaption, phylogeny, genetics, and evolution. Notably, they contain a unique class of pigments, i.e., papiliochromes, which contribute to their color diversity and various biological functions such as predator avoidance and mate preference. To date, however, the genomic and genetic basis of their color diversity and papiliochrome origin in a phylogenetic and evolutionary context remain largely unknown. Here, we obtained high-quality reference genomes of 11 swallowtail butterfly species covering all tribes of Papilioninae and Parnassiinae using long-read sequencing technology. Combined with previously published butterfly genomes, we obtained robust phylogenetic relationships among tribes, overcoming the challenges of incomplete lineage sorting (ILS) and gene flow. Comprehensive genomic analyses indicated that the evolution of Papilionidae-specific conserved non-exonic elements (PSCNEs) and transcription factor binding sites (TFBSs) of patterning and transporter/cofactor genes, together with the rapid evolution of transporters/cofactors, likely promoted the origin and evolution of papiliochromes. These findings not only provide novel insights into the genomic basis of color diversity, especially papiliochrome origin in swallowtail butterflies, but also provide important data resources for exploring the evolution, ecology, and conservation of butterflies.

The mitochondrial genome of a leaf insect Phyllium westwoodii (Phasmatodea: Phylliidae) in Southeast Asia.

The nearly complete mitochondrial genome (mitogenome) of Phyllium westwoodii, a typical leaf mimic insect in Phasmatodea, was obtained in this study. This mitogenome is 17,222 bp in length and contains 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and almost complete control regions. All PCGs initiate with 'ATN' except for NAD4L that uses 'TTG' as the start codon, and terminate with 'TAA' except for COX2 that uses a single 'T' residue as the stop codon. The phylogenetic analysis based on the concatenated sequences of 13 PCGs and two rRNAs shows that P. westwoodii is closer to Phyllium tibetense than Phyllium giganteum.

The mitochondrial genome of one 'twisted-wing parasite' Xenos cf. moutoni (Insecta, Strepsiptera, Xenidae) from Gaoligong Mountains, Southwest of China.

The nearly complete mitochondrial genome (mitogenome) of Xenos cf. moutoni, one twisted-wing parasite on wasp Vespa velutina from Southwest of China, is described in this study. The total length of this mitogenome is 16,717 bp, containing 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and an incomplete A + T-rich control region . All of the 13 PCGs are initiated with canonical ATN (N represents A, T, G, C) as start codons; 8 PCGs are terminated with a complete typical stop codon TAA, and the remaining five PCGs (cox2, cox3, nad3, nad4 and nad5) have an incomplete stop codon with just a T. The phylogenetic analysis based on the nucleotide sequences of PCGs and rRNAs indicates that Xenos cf. moutoni has a close relationship with Xenos vesparum, confirming its placement in the family Xenidae.

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.