A chromosome-level genome assembly of rice leaffolder, Cnaphalocrocis medinalis

Deep Insight into the genome of a major rice pest, Scripophaga incertulas

The yellow stem borer, Scripophaga incertulas (Crambidae, Lepidoptera), is a highly destructive pest damaging rice crops. Here, we combined Illumina, PacBio sequencing and Hi-C scaffolding to generate a high-quality chromosome-level genome assembly of S. incertulas. We initially obtained an 839.78 Mb assembly with scoffold N50 of 1.13 Mb. Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis demonstrated that this genome assembly has a high-level completeness of 97.5% gene coverage. The assembly contained 45.67% repetitive sequences (383.50 Mb), and has been annotated with 19,874 protein-coding genes. Next, 3,369 scaffolds were anchored to 26 chromosomes, resulting in an increase of the N50 to 32.16 Mb. Our research provides insights into the evolution and ecology of S. incertulas, offering valuable resources for future pest management.

Chromosome-level genome assembly of a specialist walnut pest Atrijuglans aristata

Gelechioidea represents the most diverse superfamily of tiny boring pests in Lepidoptera that pose a serious threat to agricultural and forestry economic crops. However, the lack of high-quality genome of highly specialized species makes it difficult to draw general conclusions about the mechanism of the close binding relationship between pests and crops. In this study, based on second- and third-generation sequencing reads, we constructed a chromosome-level genome for the Atrijuglans aristata, a specialized boring pest, that specifically harms the green husk of cultivated walnuts. The genome is 480.99 Mb and spans 31 pseudo-chromosomes, including Z and a portion of W chromosome. Contig N50 and scaffold N50 of the genome are 2.68 Mb and 16.01 Mb respectively. The BUSCO completeness achieves 95.6% with a total of 22,542 protein-coding genes are annotated. As the first sequenced genome of Stathmopodidae family, this high-quality genome provides a genetic basis for the mining of genes for important functional traits in A. aristata, as well as an important reference for the study of host adaptation of the (insect) specialists.

Chromosome-level genome assembly of Z strain European corn borer Ostrinia nubilalis (Lepidoptera: Crambidae)

European corn borer Ostrinia nubilalis (Hübner) is the most important pest of maize around world, and an ideal model for the sympatric host races, evolutionary and speciation research. In this study, we assembled a chromosome-level genome of Z strain O. nubilalis by the integrated Illumina short reads, PacBio Revio long reads, and Hi-C sequencing data. The chromosome-level genome was 480.04 Mb in total length with a contig N50 length of 16.51 Mb, which 98.59% genome anchored into 31 chromosomes. For the annotation, 1,046,695 repeat sequences in length of 212.07 Mb, 1,550 non-coding RNAs (including 1,208 tRNAs, 179 rRNAs, 62 miRNAs, 81 snRNAs, and 20 snoRNAs), and 17,145 protein-coding genes were identified. And 100% genes were functional annotated by SwissProt, NR, eggNOG, Go, and KEGG database. This genome provides a valuable genomics resource to elucidate the host plant adaptation, thermal adaptation, diapause induction, Bacillus thuringiensis toxin resistance, sexual communication, sympatric host races, and speciation process of O. nubilalis.

Knockdown of the glucosamine-6-phosphate N-acetyltransferase gene by RNA interference enhances the virulence of entomopathogenic fungi against rice leaffolder Cnaphalocrocis medinalis

Insect cuticle acts as a first line of defense and a physical protective barrier against entomopathogens. Chitin biosynthesis pathway plays a crucial role in chitin formation in the cuticle of insects. Glucosamine-6-phosphate N-acetyltransferase (GNA) is a key enzyme in insect chitin biosynthesis that regulate the chitin formation. However, how GNA-mediated cuticle metabolism influences virulence of entomopathogenic fungi is still unknown. In this study, CmGNA gene was cloned and characterized from the rice leaffolder Cnaphalocrocis medinalis. The CmGNA contains an open read frame (ORF) 600 nucleotides, encoding 199 amino acids with an isoelectric point of 8.65 and a molecular weight of 22.30 kDa. The expression profile showed that CmGNA was highly expressed in 4th instar larvae and in the cuticle. Here, we also reported the impact of CmGNA gene and entomopathogenic fungi, Metarhizium anisopliae and Beauveria bassiana, on expression pattern of chitin biosynthesis genes, feeding behavior, survival rate and average body weight of infected larvae, phenotypic deformities, rate of pupation, and adult emergence. Our results showed that knockdown of CmGNA and application of M. anisopliae and B. bassiana three days after RNA interference (RNAi) significantly decreased the expression of CmGNA and other associated genes, reduced feeding efficiency and survival rate, and caused loss of average body weight, less rate of pupation and adult emergence of infected larvae. Knockdown of CmGNA gene also increased the lethality of larvae caused by M. anisopliae and B. bassiana and resulted in significantly phenotypic deformities of infected larvae. Our findings illustrated that RNAi-mediated CmGNA knockdown disturbed the chitin synthesis genes that led to enhancing the virulence of M. anisopliae and B. bassiana, which can provide us new insights to develop novel biocontrol strategies against C. medinalis.

High-quality chromosome-level genome assembly of female Artemia franciscana reveals sex chromosome and Hox gene organization

Artemia is a crustacean genus belonging to the order Anostraca in the class Branchiopoda and lives in inland hypersaline lakes. Among the genus, A. franciscana is a valuable species as a fish food in the aquaculture industry or as an aquatic model organism for toxicity tests. However, genomic data for A. franciscana remains incomplete. In this study, high-quality genome assembly at the chromosome level of female A. franciscana was conducted by combining various sequencing and assembly technologies. The final A. franciscana assembled genome was 1.27 Gb in length, containing 21 chromosomal scaffolds (>10 Mb). The scaffold N50 was 45.3 Mb, with a complete BUSCO value of 91.0 %, thereby confirming that a high-quality genome was assembled. Gene annotation shows that the A. franciscana genome contained 67.26 % of repetitive sequences, and a total of 26,923 protein-coding genes were predicted. Among the 21 chromosome-scale scaffolds, chromosome 1 was identified as a sex chromosome Z. Additionally, five contigs of putative W chromosome fragments and the candidate sex-determining genes were suggested. Ten homeobox (Hox) genes were identified in A. franciscana on the chromosome 14, which were in two subclusters with a large gap. Hox gene organizations within 13 arthropods showed that four anostracans had conserved synteny. This study provides a new female Artemia genome with sex chromosome and the first complete genomic arrangement of the Hox cluster in Anostraca. This study will be a useful genomic and genetic reference for understanding the evolution and development of A. franciscana.

Chromosome-level genome assembly from a single planthopper Nilaparvata muiri (Hemiptera: Delphacidae)

The planthopper Nilaparvata muiri is a sister species to N. lugens (Hemiptera: Delphacidae), a notorious insect pest in Asian rice fields. N. muiri and N. lugens have a different host preference despite the similarities in many biological features. To better understand the adaptive evolution of planthoppers, comprehensive genomic information on N. muiri and N. lugens are urgently needed. In this study, we used ultra-low input PacBio HiFi libraries and Hi-C sequencing technologies to assemble a reference genome of a single N. muiri at the chromosomal level. The genome size was determined to be 531.62 Mb with a contig N50 size of 2.47 Mb and scaffold N50 size of 38.37 Mb. Totally, 96.61% assembled sequences were anchored to the 15 pseudo-chromosomes. BUSCO analysis yielded an Insecta completeness score of 98.6%. A total of 22,057 protein-coding genes were annotated, and 168.16 Mb repetitive sequences occupying 31.63% of genome were pinpointed. The assembled genome is valuable for evolutionary and genetic studies of planthoppers, and may provide sights to pest control.

Chromosome-level genome assembly and sex chromosome identification of the pink stem borer, Sesamia inferens (Lepidoptera: Noctuidae)

The pink stem borer, Sesamia inferens Walker (Lepidoptera: Noctuidae), is one of the most notorious pest insects of rice and maize crops in the world. Here, we generated a high-quality chromosome-level genome assembly of S. inferens, using a combination of Illumina, PacBio HiFi and Hi-C technologies. The total assembly size was 973.18 Mb with a contig N50 of 33.39 Mb, anchored to 31 chromosomes, revealing a karyotype of 30 + Z. The BUSCO analysis indicated a high completeness of 98.90% (n = 5286), including 5172 (97.8%) single-copy BUSCOs and 58 (1.1%) duplicated BUSCOs. The genome contains 58.59% (564.58 Mb) repeat elements and 26628 predicted protein-coding genes. The chromosome-level genome assembly of S. inferens provides in-depth knowledge and will be a helpful resource for the Lepidoptera and pest control research communities.

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.

Gene gain and loss from the Asian corn borer W chromosome

BACKGROUND

Sex-limited chromosomes Y and W share some characteristics, including the degeneration of protein-coding genes, enrichment of repetitive elements, and heterochromatin. However, although many studies have suggested that Y chromosomes retain genes related to male function, far less is known about W chromosomes and whether they retain genes related to female-specific function.

RESULTS

Here, we built a chromosome-level genome assembly of the Asian corn borer, Ostrinia furnacalis Guenée (Lepidoptera: Crambidae, Pyraloidea), an economically important pest in corn, from a female, including both the Z and W chromosome. Despite deep conservation of the Z chromosome across Lepidoptera, our chromosome-level W assembly reveals little conservation with available W chromosome sequence in related species or with the Z chromosome, consistent with a non-canonical origin of the W chromosome. The W chromosome has accumulated significant repetitive elements and experienced rapid gene gain from the remainder of the genome, with most genes exhibiting pseudogenization after duplication to the W. The genes that retain significant expression are largely enriched for functions in DNA recombination, the nucleosome, chromatin, and DNA binding, likely related to meiotic and mitotic processes within the female gonad.

CONCLUSIONS

Overall, our chromosome-level genome assembly supports the non-canonical origin of the W chromosome in O. furnacalis, which experienced rapid gene gain and loss, with the retention of genes related to female-specific function.

Unraveling the complex evolutionary history of lepidopteran chromosomes through ancestral chromosome reconstruction and novel chromosome nomenclature

BACKGROUND

Lepidoptera is one of the most species-rich animal groups, with substantial karyotype variations among species due to chromosomal rearrangements. Knowledge of the evolutionary patterns of lepidopteran chromosomes still needs to be improved.

RESULTS

Here, we used chromosome-level genome assemblies of 185 lepidopteran insects to reconstruct an ancestral reference genome and proposed a new chromosome nomenclature. Thus, we renamed over 5000 extant chromosomes with this system, revealing the historical events of chromosomal rearrangements and their features. Additionally, our findings indicate that, compared with autosomes, the Z chromosome in Lepidoptera underwent a fast loss of conserved genes, rapid acquisition of lineage-specific genes, and a low rate of gene duplication. Moreover, we presented evidence that all available 67 W chromosomes originated from a common ancestor chromosome, with four neo-W chromosomes identified, including one generated by fusion with an autosome and three derived through horizontal gene transfer. We also detected nearly 4000 inter-chromosomal gene movement events. Notably, Geminin is transferred from the autosome to the Z chromosome. When located on the autosome, Geminin shows female-biased expression, but on the Z chromosome, it exhibits male-biased expression. This contributes to the sexual dimorphism of body size in silkworms.

CONCLUSIONS

Our study sheds light on the complex evolutionary history of lepidopteran chromosomes based on ancestral chromosome reconstruction and novel chromosome nomenclature.

Glucose Dehydrogenases-Mediated Acclimation of an Important Rice Pest to Global Warming

Global warming is posing a threat to animals. As a large group of widely distributed poikilothermal animals, insects are liable to heat stress. How insects deal with heat stress is worth highlighting. Acclimation may improve the heat tolerance of insects, but the underlying mechanism remains vague. In this study, the high temperature of 39 °C was used to select the third instar larvae of the rice leaf folder Cnaphalocrocis medinalis, an important insect pest of rice, for successive generations to establish the heat-acclimated strain (HA39). The molecular mechanism of heat acclimation was explored using this strain. The HA39 larvae showed stronger tolerance to 43 °C than the unacclimated strain (HA27) persistently reared at 27 °C. The HA39 larvae upregulated a glucose dehydrogenase gene, CmGMC10, to decrease the reactive oxygen species (ROS) level and increase the survival rate under heat stress. The HA39 larvae maintained a higher activity of antioxidases than the HA27 when confronted with an exogenous oxidant. Heat acclimation decreased the H₂O₂ level in larvae under heat stress which was associated with the upregulation of CmGMC10. The rice leaf folder larvae may acclimate to global warming via upregulating CmGMC10 to increase the activity of antioxidases and alleviate the oxidative damage of heat stress.

Identification and characterization of ABC proteins in an important rice insect pest, Cnaphalocrocis medinalis unveil their response to Cry1C toxin

Rice leaffolder (Cnaphalocrocis medinalis) is an important insect pest in paddy fields. Due to their essential role in the physiology and insecticidal resistance, ATP-binding cassette (ABC) proteins were studied in many insects. In this study, we identified the ABC proteins in C. medinalis through genomic data and analyzed their molecular characteristics. A total of 37 sequences with nucleotide-binding domain (NBD) were identified as ABC proteins and belonged to eight families (ABCA-ABCH). Four structure styles of ABC proteins were found in C. medinalis, including full structure, half structure, single structure, and ABC2 structure. In addition to these structures, TMD-NBD-TMD, NBD-TMD-NBD, and NBD-TMD-NBD-NBD were found in C. medinalis ABC proteins. Docking studies suggested that in addition to the soluble ABC proteins, other ABC proteins including ABCC4, ABCH1, ABCG3, ABCB5, ABCG1, ABCC7, ABCB3, ABCA3, and ABCC5 binding with Cry1C had higher weighted scores. The upregulation of ABCB1 and downregulation of ABCB3, ABCC1, ABCC7, ABCG1, ABCG3, and ABCG6 were associated with the C. medinalis response to Cry1C toxin. Collectively, these results help elucidate the molecular characteristics of C. medinalis ABC proteins, pave the way for further functional studies of C. medinalis ABC proteins, including their interaction with Cry1C toxin, and provide potential insecticide targets.

Reference Genes for Expression Analysis Using RT-qPCR in Cnaphalocrocis medinalis (Lepidoptera: Pyralidae)

Cnaphalocrocis medinalis is a destructive migratory rice pest. Although many studies have investigated its behavioral and physiological responses to environmental changes and migration-inducing factors, little is known about its molecular mechanisms. This study was conducted to select suitable RT-qPCR reference genes to facilitate future gene expression studies. Here, thirteen candidate housekeeping genes (EF1α, AK, EF1β, GAPDH, PGK, RPL13, RPL18, RPS3, 18S rRNA, TBP1, TBP2, ACT, and UCCR) were selected to evaluate their stabilities under different conditions using the ∆CT method; the geNorm, NormFinder, BestKeeper algorithms; and the online tool RefFinder. The results showed that the most stable reference genes were EF1β, PGK, and RPL18, related to developmental stages; RPS3 and RPL18 in larval tissues; EF1β and PGK in larvae feeding on different rice varieties; EF1α, EF1β, and PGK in larvae temperature treatments; PGK and RPL13, related to different adult ages; PGK, EF1α, and ACT, related to adult nutritional conditions; RPL18 and PGK, related to adult mating status; and, RPS3 and PGK, related to different adult take-off characteristics. Our results reveal reference genes that apply to various experimental conditions and will greatly improve the reliability of RT-qPCR analysis for the further study of gene function in this pest.

Cuticular protein genes involve heat acclimation of insect larvae under global warming

Cuticular proteins (CPs) play important roles in insect growth and development. However, it is unknown whether CPs are related to heat tolerance. Cnaphalocrocis medinalis, a serious pest of rice, occurs in summer and exhibits strong adaptability to high temperature, but the underlying mechanism is unclear. Here, the role of CP genes in heat acclimation was studied. Heat tolerance of the heat-acclimated larvae was significantly stronger than the unacclimated larvae. The cuticular protein content in the heat-acclimated larvae was higher than that of the unacclimated larvae. 191 presumed CP genes of C. medinalis (CmCPs) were identified. Expression patterns of 14 CmCPs were different between the heat acclimated (S39) and unacclimated (S27) larvae under heat stress. CmCPs were specifically expressed in epidermis and the head except CmCPR20 mainly expressed in Malpighian tubules. CmCPR20 was upregulated in S39 while downregulated in S27, but CmTweedle1 and CmCPG1 were upregulated in S27 and downregulated in S39. RNAi CmTweedle1 or CmCPG1 remarkably decreased heat tolerance and cuticular protein content of the heat-acclimated larvae but not the unacclimated larvae. RNAi CmCPR20 decreased heat tolerance and cuticular protein content of the unacclimated larvae but not the heat-acclimated larvae. CmTweedle1 and CmCPG1 genes involve heat acclimation of C. medinalis.

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.

Novel genetic basis of resistance to Bt toxin Cry1Ac in Helicoverpa zea

Crops genetically engineered to produce insecticidal proteins from the bacterium Bacillus thuringiensis have advanced pest management, but their benefits are diminished when pests evolve resistance. Elucidating the genetic basis of pest resistance to Bacillus thuringiensis toxins can improve resistance monitoring, resistance management, and the design of new insecticides. Here, we investigated the genetic basis of resistance to Bacillus thuringiensis toxin Cry1Ac in the lepidopteran Helicoverpa zea, one of the most damaging crop pests in the United States. To facilitate this research, we built the first chromosome-level genome assembly for this species, which has 31 chromosomes containing 375 Mb and 15,482 predicted proteins. Using a genome-wide association study, fine-scale mapping, and RNA-seq, we identified a 250-kb quantitative trait locus on chromosome 13 that was strongly associated with resistance in a strain of Helicoverpa zea that had been selected for resistance in the field and lab. The mutation in this quantitative trait locus contributed to but was not sufficient for resistance, which implies alleles in more than one gene contributed to resistance. This quantitative trait locus contains no genes with a previously reported role in resistance or susceptibility to Bacillus thuringiensis toxins. However, in resistant insects, this quantitative trait locus has a premature stop codon in a kinesin gene, which is a primary candidate as a mutation contributing to resistance. We found no changes in gene sequence or expression consistently associated with resistance for 11 genes previously implicated in lepidopteran resistance to Cry1Ac. Thus, the results reveal a novel and polygenic basis of resistance.

Chromosome-level genome assembly of the mirid predator Cyrtorhinus lividipennis Reuter (Hemiptera: Miridae), an important natural enemy in the rice ecosystem

Though the genomes of many rice herbivorous pests have recently been well characterized, little is known about the genome of their natural enemies. Here, by using the Illumina and PacBio platforms, we sequenced and assembled the whole genome of the mirid species Cyrtorhinus lividipennis Reuter (Hemiptera: Miridae), which is an economically and ecologically important natural enemy in the rice ecosystem acting as a dominant predator for planthoppers and leafhoppers in the field. Through Hi-C scaffolding, 1615 scaffolds with a total size of 338.08 Mb were successfully anchored onto 13 chromosomes. The assembled genome size was 345.75 Mb with a final scaffold N50 of 27.58 Mb. Approximately 107.51 Mb of sequences accounting for 31.10% of the genome were identified as repeat elements, and 14,644 protein-coding genes were annotated. Phylogenetic analysis showed that C. lividipennis clustered with other Hemipteran species and diverged from Apolygus lucorum about 66.7 million years ago. Gene families related to detoxification, environmental adaptation and digestion were analysed comparatively with other Hemipteran species, but no significant expansion or contraction was found in C. lividipennis. We also observed male meiosis in C. lividipennis, which showed a typical post-reduction of sex chromosomes and a karyotype of 2n = 22 + XY. As the first natural-enemy genome in the rice ecosystem, the genomic resource of C. lividipennis not only expands our understanding of the multitrophic interactions (host plant-prey-predator), but also provides a genomic basis for better understanding this dominant predator and therefore promotes sustainable rice pest management and food grain production.

Chromosome-level genome assembly of an agricultural pest, the rice leaffolder Cnaphalocrocis exigua (Crambidae, Lepidoptera)

The rice leaffolder Cnaphalocrocis exigua (Crambidae, Lepidoptera) is an important agricultural pest that damages rice crops and other members of related grass families. C. exigua exhibits a very similar morphological phenotype and feeding behaviour to C. medinalis, another species of rice leaffolder whose genome was recently reported. However, genomic information for C. exigua remains extremely limited. Here, we used a hybrid strategy combining different sequencing technologies, including Illumina, PacBio, 10× Genomics, and Hi-C scaffolding, to generate a high-quality chromosome-level genome assembly of C. exigua. We initially obtained a 798.8 Mb assembly with a contig N50 size of 2.9 Mb, and the N50 size was subsequently increased to 25.7 Mb using Hi-C technology to anchor 1413 scaffolds to 32 chromosomes. We detected a total of 97.7% Benchmarking Universal Single-Copy Orthologues (BUSCO) in the genome assembly, which was comprised of ~52% repetitive sequence and annotated 14,922 protein-coding genes. Of note, the Z and W sex chromosomes were assembled and identified. A comparative genomic analysis demonstrated that despite the high synteny observed between the two rice leaffolders, the species have distinct genomic features associated with expansion and contraction of gene families and selection pressure. In summary, our chromosome-level genome assembly and comparative genomic analysis of C. exigua provide novel insights into the evolution and ecology of this rice insect pests and offer useful information for pest control.