A high-quality carabid genome assembly provides insights into beetle genome evolution and cold adaptation

Evidence for Admixture and Rapid Evolution During Glacial Climate Change in an Alpine Specialist

The pace of current climate change is expected to be problematic for alpine flora and fauna, as their adaptive capacity may be limited by small population size. Yet, despite substantial genetic drift following post-glacial recolonization of alpine habitats, alpine species are notable for their success surviving in highly heterogeneous environments. Population genomic analyses demonstrating how alpine species have adapted to novel environments with limited genetic diversity remain rare, yet are important in understanding the potential for species to respond to contemporary climate change. In this study, we explored the evolutionary history of alpine ground beetles in the Nebria ingens complex, including the demographic and adaptive changes that followed the last glacier retreat. We first tested alternative models of evolutionary divergence in the species complex. Using millions of genome-wide SNP markers from hundreds of beetles, we found evidence that the N. ingens complex has been formed by past admixture of lineages responding to glacial cycles. Recolonization of alpine sites involved a distributional range shift to higher elevation, which was accompanied by a reduction in suitable habitat and the emergence of complex spatial genetic structure. We tested several possible genetic pathways involved in adaptation to heterogeneous local environments using genome scan and genotype-environment association approaches. From the identified genes, we found enriched functions associated with abiotic stress responses, with strong evidence for adaptation to hypoxia-related pathways. The results demonstrate that despite rapid demographic change, alpine beetles in the N. ingens complex underwent rapid physiological evolution.

The chromosome-level genome assembly of the giant dobsonfly Acanthacorydalis orientalis (McLachlan, 1899)

Acanthacorydalis orientalis (McLachlan, 1899) (Megaloptera: Corydalidae) is an important freshwater-benthic invertebrate species that serves as an indicator for water-quality biomonitoring and is valuable for conservation from East Asia. Here, a high-quality reference genome for A. orientalis was constructed using Oxford Nanopore sequencing and High throughput Chromosome Conformation Capture (Hi-C) technology. The final genome size is 547.98 Mb, with the N50 values of contig and scaffold being 7.77 Mb and 50.53 Mb, respectively. The longest contig and scaffold are 20.57 Mb and 62.26 Mb in length, respectively. There are 99.75% contigs anchored onto 13 pseudo-chromosomes. Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis showed that the completeness of the genome assembly is 99.01%. There are 10,977 protein-coding genes identified, of which 84.00% are functionally annotated. The genome contains 44.86% repeat sequences. This high-quality genome provides substantial data for future studies on population genetics, aquatic adaptation, and evolution of Megaloptera and other related insect groups.

Six complete mitochondrial genomes of ground beetles from the Harpalinae and Carabinae (Coleoptera, Carabidae) with phylogenetic analysis based on mitogenomic data

In this study, we employed high-throughput sequencing technology to determine the complete mitochondrial genomes of six ground beetles, encompassing five Harpalinae species and one Carabinae species. The sizes of mitochondrial genomes ranged from 15,334 to 16,972 bp, encompassing 37 genes, including 13 protein-coding genes, 22 transfer RNA genes, and 2 ribosomal RNA genes. Furthermore, each species was found to possess a putative control region. Combining with 65 published mitochondrial genome sequences of Carabidae as ingroups and four species from Trachypachidae, Gyrinidae and Dytiscidae as outgroups, we conducted phylogenetic analyses utilizing Maximum likelihood and Bayesian inference methods. Moreover, we reconstructed a species tree of Carabidae based on mitochondrial genome data using the coalescent-based species tree method (ASTRAL). The results revealed that the family Carabidae was not a monophyletic group. The subfamily Harpalinae was supported to be a monophyletic group in Maximum likelihood analysis. Although the subfamily Carabinae was found to be nonmonophyletic in the concatenation analyses under both Maximum likelihood and Bayesian inference criteria, it was identified as a monophyletic group in the species tree analysis.

A near-complete genome assembly of Monochamus alternatus a major vector beetle of pinewood nematode

The Japanese sawyer beetle, Monochamus alternatus, is not only one of the most important wood boring pest itself, but also a major vector of the invasive pinewood nematode (PWN), which is the causal agent of the devastative pine wilt disease (PWD) and threats the global pine forest. Here, we present a near-complete genome of M. alternatus at the chromosome level. The assembled genome was 792.05 Mb with contig N50 length of 55.99 Mb, which is the largest N50 size among the sequenced Coleoptera insects currently. 99.57% of sequence was anchored onto ten pseudochromosomes (one X-chromosome and nine autosomes), and the final genome harbored only 13 gaps. BUSCO evaluation revealed the presence of 99.0% of complete core genes. Thus, our genome assembly represented the highest-contiguity genome assembly as well as high completeness in insects so far. We identified 20,471 protein-coding genes, of which 20,070 (98.04%) were functionally annotated. The genome assembly of M. alternatus provides a valuable resource for exploring the evolution of the symbiosis between PWN and the vector insects.

The toxins of vertically transmitted Spiroplasma

Vertically transmitted (VT) microbial symbionts play a vital role in the evolution of their insect hosts. A longstanding question in symbiont research is what genes help promote long-term stability of vertically transmitted lifestyles. Symbiont success in insect hosts is due in part to expression of beneficial or manipulative phenotypes that favor symbiont persistence in host populations. In Spiroplasma, these phenotypes have been linked to toxin and virulence domains among a few related strains. However, these domains also appear frequently in phylogenetically distant Spiroplasma, and little is known about their distribution across the Spiroplasma genus. In this study, we present the complete genome sequence of the Spiroplasma symbiont of Drosophila atripex, a non-manipulating member of the Ixodetis clade of Spiroplasma, for which genomic data are still limited. We perform a genus-wide comparative analysis of toxin domains implicated in defensive and reproductive phenotypes. From 12 VT and 31 non-VT Spiroplasma genomes, ribosome-inactivating proteins (RIPs), OTU-like cysteine proteases (OTUs), ankyrins, and ETX/MTX2 domains show high propensity for VT Spiroplasma compared to non-VT Spiroplasma. Specifically, OTU and ankyrin domains can be found only in VT-Spiroplasma, and RIP domains are found in all VT Spiroplasma and three non-VT Spiroplasma. These domains are frequently associated with Spiroplasma plasmids, suggesting a possible mechanism for dispersal and maintenance among heritable strains. Searching insect genome assemblies available on public databases uncovered uncharacterized Spiroplasma genomes from which we identified several spaid-like genes encoding RIP, OTU, and ankyrin domains, suggesting functional interactions among those domain types. Our results suggest a conserved core of symbiont domains play an important role in the evolution and persistence of VT Spiroplasma in insects.

Comparative genomics analysis and genome assembly integration with the recombination landscape contribute to Takifugu bimaculatus assembly refinement

Takifugu genus has been brought to the fore in scientific and practical research due to its compact genome, explosive speciation progress and economic value. Here we updated the chromosome-level genome of Takifugu bimaculatus by an ultra-high-density linkage map, a classic and accurate way of chromosome assembly. The map constituted a robust assembly frame, with 92.2% (372.77 Mb) of the draft genome cumulatively placed. With intraspecies and interspecies comparative genomic analysis, we developed a criterion to quantify the differences between assemblies and established a novel way to integrate information from multiple assemblies. The integrated assembly rectified potential mis-assemblies, greatly improving the genome contiguity and correctness. Our results rendered profound information on the genetic recombination of T. bimaculatus and provided new insights into effective genome assembly. The consolidated assembly will be a contributory tool of T. bimaculatus and broadly across the Takifugu by providing a convincing reference for genomic research.

Chromosome-level genome assembly of Monochamus saltuarius reveals its adaptation and interaction mechanism with pine wood nematode

Monochamus saltuarius (Coleoptera: Cerambycidae) was reported as the vector beetle of the pine wood nematode (PWN, Bursaphelenchus xylophilus) in Japan and Europe. It was first reported to transmitted the PWN to native Pinus species in 2018 in Liaoning Province, China. However, the lack of genomic resources has limited the in-depth understanding of its interspecific relationship with PWN. Here, we obtained a chromosome-level reference genome of M. saltuarius combining Illumina, Nanopore and Hi-C sequencing technologies. We assembled the scaffolds into ten chromosomes (including an X chromosome) and obtained a 682.23 Mb chromosome-level genome with a N50 of 73.69 Mb. In total, 427.67 Mb (62.69 %) repeat sequences were identified and 14, 492 protein-coding genes were predicted, of which 93.06 % were annotated. We described the mth/mthl, P450, OBP and OR gene families associated with the vector beetle's development and resistance, as well as the host selection and adaptation, which serve as a valuable resource for understanding the host adaptation in insects during evolution. This high quality reference genome of M. saltuarius also provide new avenues for researching the mechanism of this synergistic damage between vector beetles and PWN.

Chromosome-level genome assembly and population genomic analyses provide insights into adaptive evolution of the red turpentine beetle, Dendroctonus valens

BACKGROUND

Biological invasions are responsible for substantial environmental and economic losses. The red turpentine beetle (RTB), Dendroctonus valens LeConte, is an important invasive bark beetle from North America that has caused substantial tree mortality in China. The lack of a high-quality reference genome seriously limits deciphering the extent to which genetic adaptions resulted in a secondary pest becoming so destructive in its invaded area.

RESULTS

Here, we present a 322.41 Mb chromosome-scale reference genome of RTB, of which 98% of assembled sequences are anchored onto fourteen linkage groups including the X chromosome with a N50 size of 24.36 Mb, which is significantly greater than other Coleoptera species. Repetitive sequences make up 45.22% of the genome, which is higher than four other Coleoptera species, i.e., Mountain pine beetle Dendroctonus ponderosae, red flour beetle Tribolium castaneum, blister beetle Hycleus cichorii, and Colorado potato beetle Leptinotarsa decemlineata. We identify rapidly expanded gene families and positively selected genes in RTB, which may be responsible for its rapid environmental adaptation. Population genetic structure of RTB was revealed by genome resequencing of geographic populations in native and invaded regions, suggesting substantial divergence of the North American population and illustrates the possible invasion and spread route in China. Selective sweep analysis highlighted the enhanced ability of Chinese populations in environmental adaptation.

CONCLUSIONS

Overall, our high-quality reference genome represents an important resource for genomics study of invasive bark beetles, which will facilitate the functional study and decipher mechanism underlying invasion success of RTB by integrating the Pinus tabuliformis genome.

Illumina Short-Read Sequencing of the Mitogenomes of Novel Scarites subterraneus Isolates Allows for Taxonomic Refinement of the Genus Scarites Fabricius 1775, within the Carabidae Family

Simple Summary

Ground beetles (Carabidae) have important ecological functions and serve as food, biological pest control, and models in biological research. Although there are over 40,000 ground beetle species worldwide, only a small fraction of those has genomic information currently available, which limits their classification and understanding of diversity. Since next-generation genome sequencing has become more mainstream, we used Illumina short-read sequencing to obtain complete mitogenomes from two Scarites beetles that we collected ourselves in Nebraska and Arkansas. Scarites are large ground beetles that resemble tropical beetles with a big head and large mandibles, and their role as predator and prey helps maintain sustainability in local ecosystems. This straightforward sequencing and analysis was found to be accurate and sufficient to help classify these isolates to the subspecies level. This is the first report of mitogenomes for Scarites subterraneus and only the second in that genus. This method is easily applicable to more beetle species and can be used to increase our understanding of beetles worldwide.

Abstract

We sequenced the complete mitogenomes, 18S and 28S rRNA of two new Scarites isolates, collected in Eastern Nebraska and Northern Arkansas (US). Based on molecular sequence data comparison and morphological characteristics, the new isolates were identified as a subspecies of Scarites subterraneus Fabricius 1775, for which we propose the subspecies names ‘nebraskensis’ and ‘arkansensis’. The new 18S and 28S rRNA sequences were found to be 99% and 98% identical to Scarites subterraneus. There are no other Scarites 18S or 28S rRNA sequences in the Genbank database, however, phylogenetic analysis of the Cox1 genes showed S. vicinus Chaudoir, 1843, and S. aterrimus Morawitz, 1863, as the closest relatives. This is the first report of a mitogenome for S. subterraneus, and only the second mitogenome for that genus. The nucleotide sequence identity between the mitogenomes of the two isolates is 98.8%, while the earlier sequenced S. buparius Forster 1771 mitogenome is more distantly related, with only 90% (to ssp. nebraskensis) and 89% (to ssp. arkansensis) overall nucleotide sequence identity. These new mitogenomes, and their phylogenetic analysis, firmly establish the position of Scarites on the Carabidae family tree and further refine the genus. In addition to the molecular data provided for the Scarites species, this approach also allowed us to identify bacterial and viral signatures for Providencia, Myroides, Spiroplasma, and a giant Nucleocytoviricota virus, associated with the Scarites species. We hereby present a simple and efficient protocol for identification and phylogenetic analysis of Scarites, that is applicable to other Coleoptera, based on total DNA extraction and Illumina short-read Next-Gen sequencing.