Chromosome-level genome assembly of an oligophagous leaf beetle Ophraella communa (Coleoptera: Chrysomelidae).

The leaf beetle Ophraella communa LeSage (Coleoptera: Chrysomelidae) is an effective biological control agent of the common ragweed. Here, we assembled a chromosome-level genome of the O. communa by combining Illumina, Nanopore, and Hi-C sequencing technologies. The genome size of the final genome assembly is 733.1 Mb, encompassing 17 chromosomes, with an improved contig N50 of 7.05 Mb compared to the original version. Genome annotation reveals 25,873 protein-coding genes, with functional annotations available for 22,084 genes (85.35%). Non-coding sequence annotation identified 204 rRNAs, 626 tRNAs, and 1791 small RNAs. Repetitive elements occupy 414.41 Mb, constituting 57.76% of the genome. This high-quality genome is fundamental for advancing biological control strategies employing O. communa.

Is systemic lupus erythematosus linked to Immunoglobulin G4 Autoantibodies?

Systemic lupus erythematosus (SLE) is a systemic autoimmune disorder characterized by a hyperactive immune system with multiple abnormalities in B-cell proliferation, antibody production, T-cell regulation, and immune complex (IC) formation. In humans, Immunoglobulin (Ig) G is found in four subclasses. IgG1-IgG4, which are distinguished by both structural and biological differences. Fab-arm Exchange (FAE), specific biases in the IgG4 response repertoire, and a decreased capacity to induce effector functions mediated by interactions in the fragment crystallizable (Fc) region are just a few of the distinctive characteristics of IgG4. The recent finding of the presence of double-stranded DNA (dsDNA) and antinuclear antibody (ANA)-IgG4 has raised attention to this IgG subclass and its possible role in SLE. IgG4 was previously believed to just have anti-inflammatory effects by inhibiting immune responses, but recent studies have shown that these antibodies can also play a role in the onset and development of some clinical disorders. To consider the clinical effects of IgG4 presence, it is necessary to discuss its characteristics, which could underlie the potential role it can play in SLE. Therefore, this study aimed to comprehensively review the role of IgG4 in SLE to elucidate the collective incidence of high IgG4 levels reported in some SLE patients.

Nuclear and mitochondrial genomes of the plum fruit moth Grapholita funebrana.

The plum fruit moth Grapholita funebrana (Tortricidae, Lepidoptera) is an important pest of many wild and cultivated stone fruits and other plants in the family Rosaceae. Here, we assembled its nuclear and mitochondrial genomes using Illumina, Nanopore, and Hi-C sequencing technologies. The nuclear genome size is 570.9 Mb, with a repeat rate of 51.28%, and a BUCSO completeness of 97.7%. The karyotype for males is 2n = 56. We identified 17,979 protein-coding genes, 5,643 tRNAs, and 94 rRNAs. We also determined the mitochondrial genome of this species and annotated 13 protein-coding genes, 22 tRNAs, and 2 rRNA. These genomes provide resources to understand the genetics, ecology, and genome evolution of the tortricid moths.

Chromosome-level genome assembly of the western flower thrips Frankliniella occidentalis.

The western flower thrips Frankliniella occidentalis (Thysanoptera: Thripidae) is a global invasive species that causes increasing damage by direct feeding on crops and transmission of plant viruses. Here, we assemble a previously published scaffold-level genome into a chromosomal level using Hi-C sequencing technology. The assembled genome has a size of 302.58 Mb, with a contig N50 of 1533 bp, scaffold N50 of 19.071 Mb, and BUSCO completeness of 97.8%. All contigs are anchored on 15 chromosomes. A total of 16,312 protein-coding genes are annotated in the genome with a BUSCO completeness of 95.2%. The genome contains 492 non-coding RNA, and 0.41% of interspersed repeats. In conclusion, this high-quality genome provides a convenient and high-quality resource for understanding the ecology, genetics, and evolution of thrips.

Mitochondrial genome rearrangements and phylogenomics of the Hymenoptera (Insecta) using an expanded taxon sample.

The order Hymenoptera is one of the most species-rich insect orders, with more than 150,000 described extant species. Many hymenopteran insects have very different mitochondrial genome (mitogenome) organizations compared to the putative ancestral organization of insects. In this study, we sequenced 18 mitogenomes of representatives in the order Hymenoptera to increase taxonomic sampling. A total of 475 species were used in phylogenetic analyses, including 18 new mitogenomes and 457 existing mitogenomes. Using a site-heterogeneous model, Bayesian's inference from amino acid data yielded more resolved relationships among Hymenoptera than maximum-likelihood analysis and coalescent-based species analyses. The monophyly of Symphyta was not supported. The Xyeloidea was the earliest branching clade in the Hymenoptera. The Orussoidea was closely related to Apocrita. Within Apocrita, the Parasitoida was non-monophyletic. The monophyly of most Parasitoida superfamilies received strong support. The Proctotrupomorpha clade was supported in Bayesian's analysis. The Apoidea was monophyletic when excluding Ampulex compressa from consideration. The superfamilies Vespoidea and Chrysidoidea were found to be non-monophyletic. Comparisons of mitochondrial gene order revealed a higher frequency of gene rearrangement among lineages with a parasitoid lifestyle, particularly prominent in Chalcidoidea. The degree of gene rearrangement ranked second in specific taxa of Cynipoidea and Ichneumonoidea.

Nuclear and mitochondrial genomes of Polylopha cassiicola: the first assembly in Chlidanotinae (Tortricidae).

Tortricidae is one of the largest families in Lepidoptera, including subfamilies of Tortricinae, Olethreutinae, and Chlidanotinae. Here, we assembled the gap-free genome for the subfamily Chlidanotinae using Illumina, Nanopore, and Hi-C sequencing from Polylopha cassiicola, a pest of camphor trees in southern China. The nuclear genome is 302.03 Mb in size, with 36.82% of repeats and 98.4% of BUCSO completeness. The karyotype is 2n = 44 for males. We identified 15412 protein-coding genes, 1052 tRNAs, and 67 rRNAs. We also determined the mitochondrial genome of this species and annotated 13 protein-coding genes, 22 tRNAs, and one rRNA. These high-quality genomes provide valuable information for studying phylogeny, karyotypic evolution, and adaptive evolution of tortricid moths.

A chromosome-level genome for the flower thrips Frankliniella intonsa.

The flower thrips Frankliniella intonsa (Thysanoptera: Thripidae) is a common insect found in flowers of many plants. Sometimes, F. intonsa causes damage to crops through direct feeding and transmission of plant viruses. Here, we assembled a chromosomal level genome of F. intonsa using the Illumina, Oxford Nanopore (ONT), and Hi-C technologies. The assembled genome had a size of 209.09 Mb, with a contig N50 of 997 bp, scaffold N50 of 13.415 Mb, and BUSCO completeness of 92.5%. The assembled contigs were anchored on 15 chromosomes. A set of 14,109 protein-coding genes were annotated in the genome with a BUSCO completeness of 95.0%. The genome contained 491 non-coding RNA and 0.57% of interspersed repeats. This high-quality genome provides a valuable resource for understanding the ecology, genetics, and evolution of F. intonsa, as well as for controlling thrips pests.

Plasticity of cold and heat stress tolerance induced by hardening and acclimation in the melon thrips.

Extreme temperatures threaten species under climate change and can limit range expansions. Many species cope with changing environments through plastic changes. This study tested phenotypic changes in heat and cold tolerance under hardening and acclimation in the melon thrips, Thrips palmi Karny (Thysanoptera: Thripidae), an agricultural pest of many vegetables. We first measured the critical thermal maximum (CTmax) of the species by the knockdown time under static temperatures and found support for an injury accumulation model of heat stress. The inferred knockdown time at 39 °C was 82.22 min. Rapid heat hardening for 1 h at 35 °C slightly increased CTmax by 1.04 min but decreased it following exposure to 31 °C by 3.46 min and 39 °C by 6.78 min. Heat acclimation for 2 and 4 days significantly increased CTmax at 35 °C by 1.83, and 6.83 min, respectively. Rapid cold hardening at 0 °C and 4 °C for 2 h, and cold acclimation at 10 °C for 3 days also significantly increased cold tolerance by 6.09, 5.82, and 2.00 min, respectively, while cold hardening at 8 °C for 2 h and acclimation at 4 °C and 10 °C for 5 days did not change cold stress tolerance. Mortality at 4 °C for 3 and 5 days reached 24.07 % and 43.22 % respectively. Our study showed plasticity for heat and cold stress tolerance in T. palmi, but the thermal and temporal space for heat stress induction is narrower than for cold stress induction.

Rapid and Repeated Climate Adaptation Involving Chromosome Inversions following Invasion of an Insect.

Following invasion, insects can become adapted to conditions experienced in their invasive range, but there are few studies on the speed of adaptation and its genomic basis. Here, we examine a small insect pest, Thrips palmi, following its contemporary range expansion across a sharp climate gradient from the subtropics to temperate areas. We first found a geographically associated population genetic structure and inferred a stepping-stone dispersal pattern in this pest from the open fields of southern China to greenhouse environments of northern regions, with limited gene flow after colonization. In common garden experiments, both the field and greenhouse groups exhibited clinal patterns in thermal tolerance as measured by critical thermal maximum (CTmax) closely linked with latitude and temperature variables. A selection experiment reinforced the evolutionary potential of CTmax with an estimated h2 of 6.8% for the trait. We identified 3 inversions in the genome that were closely associated with CTmax, accounting for 49.9%, 19.6%, and 8.6% of the variance in CTmax among populations. Other genomic variations in CTmax outside the inversion region were specific to certain populations but functionally conserved. These findings highlight rapid adaptation to CTmax in both open field and greenhouse populations and reiterate the importance of inversions behaving as large-effect alleles in climate adaptation.

Chromosome-level genome assembly of the Japanese sawyer beetle Monochamus alternatus.

The Japanese sawyer beetle Monochamus alternatus (Coleoptera: Cerambycidae) is a pest in pine forests and acts as a vector for the pine wood nematode Bursaphelenchus xylophilus, which causes the pine wilt disease. We assembled a high-quality genome of M. alternatus at the chromosomal level using Illumina, Nanopore, and Hi-C sequencing technologies. The assembled genome is 767.12 Mb, with a scaffold N50 of 82.0 Mb. All contigs were assembled into ten pseudo-chromosomes. The genome contains 63.95% repeat sequences. We identify 16, 284 protein-coding genes in the genome, of which 11,244 were functionally annotated. The high-quality genome of M. alternatus provides an invaluable resource for the biological, ecological, and genetic study of this beetle and opens new avenues for understanding the transmission of pine wood nematode by insect vectors.

The complete mitochondrial genome of Aphidius colemani (Hymenoptera: Braconidae: Aphidiinae).

The genome-level features are crucial genetic resources for species identification and phylogenetic analysis. Here, the complete mitochondrial genome of Aphidius colemani Viereck 1912 (Hymenoptera: Braconidae: Aphidiinae) was sequenced, determined and analyzed. The circular genome is 16,372 bp in length with an overall base composition of 38.9% for A, 46.2% for T, 6.7% for C, and 8.2% for G. The mitochondrial genome of A. colemani contained 13 protein-coding genes that initiated by the ATN codon, 22 transfer RNA genes, two ribosomal RNA genes (rRNAs), and a control region (CR). It shared the same gene arrangement patterns that occurred in two tRNA clusters of trnI-trnQ-trnM and trnW-trnC-trnY with Aphidius gifuensis. Phylogenetic analyses using Bayesian inference and Maximum-likelihood methods supported that the two species of Aphidiinae formed a clade and sister to other subfamilies of Braconidae.

Stable composition of gut microbiome in the Asian ladybeetle Coccinella septempunctata reared on natural and artificial diets.

The Asian ladybeetle, Coccinella septempunctata, is an important insect of predatory natural enemy, which has a strong control effect and application prospects for aphids, whiteflies, mealybugs, and other small-sized pests of agriculture and forestry crops. Gut microbiota composition posed impacts on development of insects. In order to clarify the effect of artificial feed feeding on the intestinal microbial species and structure of C. septempunctata, we compared the intestinal microbial flora of C. septempunctata reared on bean aphids and artificial food for 15 days. Results show that Proteobacteria was the dominant component in all groups at phylum level, Rhodobacter, Methylovigula, Burkholderia, and Bradyrhizobium were the dominant bacteria among all groups at genus level. As to the differences in bacterial community structure and diversity, there is no significant difference between Shannon index and Simpson index, the principal components analysis of the bacterial communities, and the samples were roughly distributed in different regions. After 15 days of feeding, artificial diet did not significantly reduce the microbial diversity of the gut of C. septempunctata compared to the aphid group, and there was no significant effect on the abundance of dominant flora in the gut of C. septempunctata, C. septempunctata gut has a similar core microbiota. This study clarifies the effects in intestinal microbial diversity and composition structure of the C. septempunctata with artificial diet, and provides a theoretical basis for understanding the intestinal microorganisms and optimizating the artificial diet of C. septempunctata.

Chromosome-level genome assembly of the giant ladybug Megalocaria dilatata.

The giant ladybug Megalocaria dilatata (Fabricius) is a potential biocontrol agent and a valuable model for coccinellid genomics and evolutionary biology. However, the lack of a reference genome for M. dilatata has impeded further explorations into its evolution and constrained its use in pest management. Here, we assembled and annotated a high-quality, chromosome-level genome of M. dilatata. The resulting assembly spans 772.3 Mb, with a scaffold N50 of 72.48 Mb and a GC content of 34.23%. The Hi-C data aided in anchoring the assembly onto 10 chromosomes ranging from 43.35 to 108.16 Mb. We identified 493.33 Mb of repeat sequences, accounting for 63.88% of the assembled genome. Our gene prediction identified 25,346 genes, with 81.89% annotated in public protein databases. The genome data will provide a valuable resource for studying the biology and evolution of Coccinellidae, aiding in pest control strategies and advancing research in the field.

Chromosome-level genome assembly of an agricultural pest Zeugodacus tau (Diptera: Tephritidae).

The fruit fly Zeugodacus tau (Diptera: Tephritidae) is a major pest of melons and other cucurbits in Southeast Asia. In this study, we used Illumina, Nanopore, and Hi-C sequencing technologies to assemble a reference genome of Z. tau at the chromosomal level. The assembled genome was 421.79 Mb and consisted of six chromosomes (one X-chromosome + five autosomes). The contig N50 was 4.23 Mb. We identified 20,922 protein-coding genes, of which 17,251 (82.45%) were functionally annotated. Additionally, we found 247 rRNAs, 435 tRNAs, 67 small nuclear RNAs, and 829 small RNAs in the genome. Repetitive elements accounted for 55.30 Mb (13.15%) of the genome. This high-quality genome assembly is valuable for evolutionary and genetic studies of Z. tau and its relative species.

Effects of fungicides on fitness and Buchnera endosymbiont density in Aphis gossypii.

BACKGROUND: Several agricultural fungicides are known to affect insect pests directly and these effects may be transgenerational and mediated through impacts on endosymbionts, providing opportunities for pest control. The cotton aphid Aphis gossypii is a polyphagous pest that can cause large crop yield losses. Here, we tested the effects of three fungicides, pyraclostrobin, trifloxystrobin and chlorothalonil, on the fitness and Buchnera endosymbiont of A. gossypii. RESULTS: The formulations of trifloxystrobin and pyraclostrobin, and the active ingredient of pyraclostrobin produced dose-dependent mortality in A. gossypii, whereas there was no dose-dependent mortality for chlorothalonil. The formulations of trifloxystrobin and pyraclostrobin significantly reduced the lifespan and fecundity of A. gossypii, and increased the density of Buchnera in the parental generation but not the (unexposed) F1 . When the active ingredient of pyraclostrobin was tested, the lifespan of the F0 generation was also reduced, but the density of Buchnera was not, indicating that non-insecticidal chemicals in the fungicide formulation may affect the density of the endosymbiont of A. gossypii. There was no transgenerational effect of the active ingredient of pyraclostrobin on the lifespan and Buchnera of (unexposed) F1 . CONCLUSIONS: Our results suggest that formulations of two strobilurin fungicides have immediate impacts on the fitness of A. gossypii, and chemicals in the formulation impact the density of the primary Buchnera endosymbiont. Our study highlights the potential effects of non-insecticidal chemicals of fungicides on aphid pests and their primary endosymbionts but direct connections between fitness and Buchnera densities remain unclear. © 2023 Society of Chemical Industry.

Testing for adaptive changes linked to range expansion following a single introduction of the fall webworm.

Adaptive evolution following colonization can affect the impact of invasive species. The fall webworm (FWW) invaded China 40 years ago through a single introduction event involving a severe bottleneck and subsequently diverged into two genetic groups. The well-recorded invasion history of FWW, coupled with a clear pattern of genetic divergence, provides an opportunity to investigate whether there is any sign of adaptive evolution following the invasion. Based on genome-wide SNPs, we identified genetically separated western and eastern groups of FWW and correlated spatial variation in SNPs with geographical and climatic factors. Geographical factors explained a similar proportion of the genetic variation across all populations compared with climatic factors. However, when the two population groups were analysed separately, environmental factors explained more variation than geographical factors. SNP outliers in populations of the western group had relatively stronger response to precipitation than temperature-related variables. Functional annotation of SNP outliers identified genes associated with insect cuticle protein potentially related to desiccation adaptation in the western group and genes associated with lipase biosynthesis potentially related to temperature adaptation in the eastern group. Our study suggests that invasive species may maintain the evolutionary potential to adapt to heterogeneous environments despite a single invasion event. The molecular data suggest that quantitative trait comparisons across environments would be worthwhile.

Interspecific and intraspecific variation in susceptibility of two co-occurring pest thrips, Frankliniella occidentalis and Thrips palmi, to nine insecticides.

BACKGROUND: Field control of pest thrips mainly relies on insecticides, but the toxicity of insecticides can vary among thrips species and populations. In this study, we examined the susceptibility of multiple field populations of two thrips pests, Frankliniella occidentalis, and Thrips palmi, that often co-occur on vegetables, to nine insecticides belonging to seven subgroups. RESULTS: The highest level of variation in susceptibility among F. occidentalis populations was for spinetoram (73.92 fold difference between most resistant and most susceptible population), followed by three neonicotinoids (8.06-15.99 fold), while among T. palmi populations, it was also for spinetoram (257.19 fold), followed by emamectin benzoate, sulfoxaflor, and acetamiprid (23.64-45.50 fold). These findings suggest evolved resistance to these insecticides in some populations of the two thrips. One population of F. occidentalis had a particularly high level of resistance overall, being the most resistant for five of the nine insecticides tested. Likewise, a population of T. palmi had high resistance to all nine insecticides, again suggesting the evolution of resistance to multiple chemicals. For F. occidentalis, the LC95 values of most populations were higher than the field-recommended dosage for all insecticides except chlorfenapyr and emamectin benzoate. For several T. palmi populations, the LC95 values also tended to be higher than recommended dosages, except in the case of emamectin benzoate and spinetoram. CONCLUSIONS: Our study found interspecific and intraspecific variations in the susceptibility of two thrips to nine insecticides and multiple resistance in some populations, highlighting the need for ongoing monitoring and resistance management. © 2023 Society of Chemical Industry.

The Complete Mitochondrial Genome of the Chinese White Wax Scale Insect, Ericerus pela Chavannes (Hemiptera: Coccidae), with Novel Gene Arrangement and Truncated tRNA Genes.

The Chinese white wax scale insect, Ericerus pela Chavannes (Hemiptera: Coccidae), is one of the scale insects with great economic value and has been dispersed and reared in China for over one thousand years. Its mitochondrial genome provides essential information for the molecular identification and genetic study of this species. We assembled the complete mitochondrial genome of E. pela based on PacBio sequencing and analyzed its genomic features. The genome was 17,766 bp in length with 13 protein-coding genes, 22 tRNAs, and two rRNA genes. The analysis results showed E. pela had significant gene rearrangements involving tRNAs compared with other Coccoidea species. Furthermore, E. pela's nine tRNAs were identified to have obvious truncated structures. The phylogenetic tree compiled of the species showed a long branch of the Coccoidea lineage, which indicated the high evolutionary rate in this group. Our study revealed the mitochondrial characteristics of E. pela and enriched the mitochondrial genetic information on Coccoidea species. It also determined the occurrence of gene rearrangement for the species in this superfamily.

A high-throughput KASP assay provides insights into the evolution of multiple resistant mutations in populations of the two-spotted spider mite Tetranychus urticae across China.

BACKGROUND: The two-spotted spider mite (TSSM), Tetranychus urticae (Acari: Tetranychidae), is a cosmopolitan phytophagous pest in agriculture and horticulture. It has developed resistance to many acaricides by target-site mutations. Understanding the status and evolution of resistant mutations in the field is essential for resistance management. Here, we applied a high-throughput Kompetitive allele-specific polymerase chain reaction (KASP) method for detecting six mutations conferring resistance to four acaricides of the TSSM. We genotyped 3274 female adults of TSSM from 43 populations collected across China in 2017, 2020, and 2021. RESULTS: The KASP genotyping of 24 testing individuals showed 99% agreement with Sanger sequencing results. KASP assays showed that most populations had a high frequency of mutations conferring avermectin (G314D and G326E) and pyridaben (H92R) resistance. The frequency of mutation conferring bifenazate (A269V and G126S) and etoxazole (I1017F) resistance was relatively low. Multiple mutations were common in the TSSM, with 70.2% and 24.6% of individuals having 2-6 and 7-10 of 10 possible resistant alleles, respectively. No loci were linked in most populations among the six mutations, indicating the development of multiple resistance is mainly by independent selection. However, G314D and I1017F on the nuclear genome deviated from Hardy-Weinberg equilibrium in most populations, indicating significant selective pressure on TSSM populations by acaricides or fitness cost of the mutations in the absence of acaricide selection. CONCLUSION: Our study revealed that the high frequency of TSSMs evolved multiple resistant mutations in population and individual levels by independent selection across China, alarming for managing multiple-acaricides resistance. © 2023 Society of Chemical Industry.

KASP genotyping and semi-quantitation of G275E mutation in the α6 subunit of Thrips palmi nAChR gene conferring spinetoram resistance.

BACKGROUND: Pesticide resistance is a long-standing and growing problem in the chemical control of invertebrate pests. Molecular diagnostic methods can facilitate pesticide resistance management by accurately and efficiently detecting resistant mutations and their frequency. In this study, the kompetitive allele specific PCR (KASP) approach, a technology for high-throughput single nucleotide polymorphism (SNP) genotyping, is validated as a useful method for characterizing genotypes at a pesticide-resistance locus for the first time. We focus on the spinetoram resistance mutation of G275E in the nicotinic acetylcholine receptor alpha 6 (nAChR α6) subunit gene of Thrips palmi. RESULTS: Of the 341 individuals of Thrips palmi tested, 98.24% were successfully genotyped, with 100% concordance with Sanger sequencing results. We then quantitatively mixed genomic DNA of known genotypes to establish 21 DNA mixtures with a resistant allele frequency ranging from 0 to 100% at steps of 5%. The linear discriminant analysis (LDA) showed that 75.8% of original grouped cases were correctly classified; six groups had no overlap in membership (resistant allele frequency: 0%, 5%, 10-75%, 80-85%, 90-95%, and 100%). When we chose 11 pooled samples with 10% steps for LDA, 84.4% of original grouped cases were correctly classified; seven groups had no overlap in membership (0%, 10%, 20-30%, 40-70%, 80%, 90%, 100%). The results indicated that KASP applied to pooled samples may provide a semi-quantitative estimate of resistance. CONCLUSIONS: Our study points to the suitability of KASP for high-throughput genotyping of genotypes affecting pesticide resistance and semi-quantitative assessments of resistance allele frequencies in populations. © 2023 Society of Chemical Industry.