The Landscapes of Full-Length Transcripts and Splice Isoforms as Well as Transposons Exonization in the Lepidopteran Model System, Bombyx mori

Multiple independent origins of the female W chromosome in moths and butterflies

Lepidoptera, the most diverse group of insects, exhibit female heterogamy (Z0 or ZW), which is different from most other insects (male heterogamy, XY). Previous studies suggest a single origin of the Z chromosome. However, the origin of the lepidopteran W chromosome remains poorly understood. Here, we assemble the genome from females down to the chromosome level of a model insect (Bombyx mori) and identify a W chromosome of approximately 10.1 megabase using a newly developed tool. In addition, we identify 3593 genes that were not previously annotated in the genomes of B. mori. Comparisons of 21 lepidopteran species (including 17 ZW and four Z0 systems) and three trichopteran species (Z0 system) reveal that the formation of Ditrysia W involves multiple mechanisms, including previously proposed canonical and noncanonical models, as well as a newly proposed mechanism called single-Z turnover. We conclude that there are multiple independent origins of the W chromosome in the Ditrysia (most moths and all butterflies) of Lepidoptera.

SQANTI3: curation of long-read transcriptomes for accurate identification of known and novel isoforms

The emergence of long-read RNA sequencing (lrRNA-seq) has provided an unprecedented opportunity to analyze transcriptomes at isoform resolution. However, the technology is not free from biases, and transcript models inferred from these data require quality control and curation. In this study, we introduce SQANTI3, a tool specifically designed to perform quality analysis on transcriptomes constructed using lrRNA-seq data. SQANTI3 provides an extensive naming framework to describe transcript model diversity in comparison to the reference transcriptome. Additionally, the tool incorporates a wide range of metrics to characterize various structural properties of transcript models, such as transcription start and end sites, splice junctions, and other structural features. These metrics can be utilized to filter out potential artifacts. Moreover, SQANTI3 includes a Rescue module that prevents the loss of known genes and transcripts exhibiting evidence of expression but displaying low-quality features. Lastly, SQANTI3 incorporates IsoAnnotLite, which enables functional annotation at the isoform level and facilitates functional iso-transcriptomics analyses. We demonstrate the versatility of SQANTI3 in analyzing different data types, isoform reconstruction pipelines, and sequencing platforms, and how it provides novel biological insights into isoform biology. The SQANTI3 software is available at https://github.com/ConesaLab/SQANTI3 .

A Blueprint of Microstructures and Stage-Specific Transcriptome Dynamics of Cuticle Formation in Bombyx mori

Insect cuticle is critical for the environmental adaptability and insecticide resistance of insects. However, there is no clear understanding of the structure and protein components of the cuticle during each developmental stage of holometabolous insects, and knowledge about the protein components within each layer is vague. We conducted serial sectioning, cuticular structure analysis, and transcriptome sequencing of the larval, pupal, and adult cuticles of Bombyx mori. The deposition processes of epicuticle, exocuticle, and endocuticle during larval, pupal, and adult cuticle formation were similar. Transcriptome analysis showed that these cuticle formations share 74% of the expressed cuticular protein (CP) genes and 20 other structural protein genes, such as larval serum protein and prisilkin. There are seven, six, and eleven stage-specific expressed CP genes in larval, pupal, and adult cuticles, respectively. The types and levels of CP genes may be the key determinants of the properties of each cuticular layer. For example, the CPs of the RR-2 protein family with high contents of histidine (His) are more essential for the exocuticle. Functional analysis suggested that BmorCPAP1-H is involved in cuticle formation. This study not only offers an in-depth understanding of cuticle morphology and protein components but also facilitates the elucidation of molecular mechanisms underlying cuticle formation in future studies.

The Role of Chitooligosaccharidolytic β-N-Acetylglucosamindase in the Molting and Wing Development of the Silkworm Bombyx mori

The insect glycoside hydrolase family 20 β-N-acetylhexosaminidases (HEXs) are key enzymes involved in chitin degradation. In this study, nine HEX genes in Bombyx mori were identified by genome-wide analysis. Bioinformatic analysis based on the transcriptome database indicated that each gene had a distinct expression pattern. qRT-PCR was performed to detect the expression pattern of the chitooligosaccharidolytic β-N-acetylglucosaminidase (BmChiNAG). BmChiNAG was highly expressed in chitin-rich tissues, such as the epidermis. In the wing disc and epidermis, BmChiNAG has the highest expression level during the wandering stage. CRISPR/Cas9-mediated BmChiNAG deletion was used to study the function. In the BmChiNAG-knockout line, 39.2% of female heterozygotes had small and curly wings. The ultrastructure of a cross-section showed that the lack of BmChiNAG affected the stratification of the wing membrane and the formation of the correct wing vein structure. The molting process of the homozygotes was severely hindered during the larva to pupa transition. Epidermal sections showed that the endocuticle of the pupa was not degraded in the mutant. These results indicate that BmChiNAG is involved in chitin catabolism and plays an important role in the molting and wing development of the silkworm, which highlights the potential of BmChiNAG as a pest control target.