Firefly genomes illuminate parallel origins of bioluminescence in beetles

Genomics investigation of the potentially invasive firefly Photinus signaticollis Blanchard 1845: Complete mitochondrial genome, multigene phylogenies and obtention of the luciferase and luciferin-regenerating genes

A genomic investigation of the potentially invasive firefly Photinus signaticollis Blanchard1845 has been performed and led to the obtention of its complete 16,411 bp long mitochondrial genome. The mitogenome encodes 13 protein-coding genes, 22 tRNA genes and 2 rRNA genes. With other species of the Photinus complex it shares several premature terminations of some protein-coding genes and also an overlap between cox1 and tRNA-Tyr. By data-mining, the complete luciferase and luciferin-regenerating genes were also identified from the contigs file and compared with existing data, in addition to WG and CAD, two genes used in pioneering phylogenetic studies on fireflies. Three maximum likelihood phylogenies were derived from all these data. The multigene phylogeny based on all mitochondrial protein-coding genes strongly associates P. signaticollis with Photinus pyralis Linnaeus, 1758 and the lantern-less daily "winter firefly", Photinus corruscus Linnaeus, 1767. A second phylogeny based on concatenated sequences of the cox1, WG and CAD genes positions P. signaticollis as a sister clade to a large cluster of species containing the 7 sub-groups previously evidenced among the North American species of the Photinus complex. A third phylogeny based on the amino-acid sequence of the luciferase protein associates P. signaticollis to Photinus scintillans. The analysis presented here will most certainly help to come to a better understanding of the very complex inter-relationships in the very large Photinus genus.

The grasshopper genome reveals long-term gene content conservation of the X Chromosome and temporal variation in X Chromosome evolution

We present the first chromosome-level genome assembly of the grasshopper, Locusta migratoria, one of the largest insect genomes. We use coverage differences between females (XX) and males (X0) to identify the X Chromosome gene content, and find that the X Chromosome shows both complete dosage compensation in somatic tissues and an underrepresentation of testis-expressed genes. X-linked gene content from L. migratoria is highly conserved across seven insect orders, namely Orthoptera, Odonata, Phasmatodea, Hemiptera, Neuroptera, Coleoptera, and Diptera, and the 800 Mb grasshopper X Chromosome is homologous to the fly ancestral X Chromosome despite 400 million years of divergence, suggesting either repeated origin of sex chromosomes with highly similar gene content, or long-term conservation of the X Chromosome. We use this broad conservation of the X Chromosome to test for temporal dynamics to Fast-X evolution, and find evidence of a recent burst evolution for new X-linked genes in contrast to slow evolution of X-conserved genes.

Two novel genomes of fireflies with different degrees of sexual dimorphism reveal insights into sex-biased gene expression and dosage compensation

Sexual dimorphism arises because of divergent fitness optima between the sexes. Phenotypic divergence between sexes can range from mild to extreme. Fireflies, bioluminescent beetles, present various degrees of sexual dimorphism, with species showing very mild sexual dimorphism to species presenting female-specific neoteny, posing a unique framework to investigate the evolution of sexually dimorphic traits across species. In this work, we present novel assembled genomes of two firefly species, Lamprohiza splendidula and Luciola italica, species with different degrees of sexual dimorphism. We uncover high synteny conservation of the X-chromosome across ~ 180 Mya and find full X-chromosome dosage compensation in our two fireflies, hinting at common mechanism upregulating the single male X-chromosome. Different degrees of sex-biased expressed genes were found across two body parts showing different proportions of expression conservation between species. Interestingly, we do not find X-chromosome enrichment of sex-biased genes, but retrieve autosomal enrichment of sex-biased genes. We further uncover higher nucleotide diversity in the intronic regions of sex-biased genes, hinting at a maintenance of heterozygosity through sexual selection. We identify different levels of sex-biased gene expression divergence including a set of genes showing conserved sex-biased gene expression between species. Divergent and conserved sex-biased genes are good candidates to test their role in the maintenance of sexually dimorphic traits.

Firefly toxin lucibufagins evolved after the origin of bioluminescence

Fireflies were believed to originally evolve their novel bioluminescence as warning signals to advertise their toxicity to predators, which was later adopted in adult mating. Although the evolution of bioluminescence has been investigated extensively, the warning signal hypothesis of its origin has not been tested. In this study, we test this hypothesis by systematically determining the presence or absence of firefly toxin lucibufagins (LBGs) across firefly species and inferring the time of origin of LBGs. We confirm the presence of LBGs in the subfamily Lampyrinae, but more importantly, we reveal the absence of LBGs in other lineages, including the subfamilies of Luciolinae, Ototretinae, and Psilocladinae, two incertae sedis lineages, and the Rhagophthalmidae family. Ancestral state reconstructions for LBGs based on firefly phylogeny constructed using genomic data suggest that the presence of LBGs in the common ancestor of the Lampyrinae subfamily is highly supported but unsupported in more ancient nodes, including firefly common ancestors. Our results suggest that firefly LBGs probably evolved much later than the evolution of bioluminescence. We thus conclude that firefly bioluminescence did not originally evolve as direct warning signals for toxic LBGs and advise that future studies should focus on other hypotheses. Moreover, LBG toxins are known to directly target and inhibit the α subunit of Na+, K+-ATPase (ATPα). We further examine the effects of amino acid substitutions in firefly ATPα on its interactions with LBGs. We find that ATPα in LBG-containing fireflies is relatively insensitive to LBGs, which suggests that target-site insensitivity contributes to LBG-containing fireflies' ability to deal with their own toxins.

Evolution of bioluminescence in Anthozoa with emphasis on Octocorallia

Bioluminescence is a widespread phenomenon that has evolved multiple times across the tree of life, converging among diverse fauna and habitat types. The ubiquity of bioluminescence, particularly in marine environments where it is commonly used for communication and defense, highlights the adaptive value of this trait, though the evolutionary origins and timing of emergence remain elusive for a majority of luminous organisms. Anthozoan cnidarians are a diverse group of animals with numerous bioluminescent species found throughout the world's oceans, from shallow waters to the light-limited deep sea where bioluminescence is particularly prominent. This study documents the presence of bioluminescent Anthozoa across depth and explores the diversity and evolutionary origins of bioluminescence among Octocorallia-a major anthozoan group of marine luminous organisms. Using a phylogenomic approach and ancestral state reconstruction, we provide evidence for a single origin of bioluminescence in Octocorallia and infer the age of occurrence to around the Cambrian era, approximately 540 Ma-setting a new record for the earliest timing of emergence of bioluminescence in the marine environment. Our results further suggest this trait was largely maintained in descendants of a deep-water ancestor and bioluminescent capabilities may have facilitated anthozoan diversification in the deep sea.

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.

Genome assembly of Genji firefly (Nipponoluciola cruciata) reveals novel luciferase-like luminescent proteins without peroxisome targeting signal

The Genji firefly, Nipponoluciola cruciata, is an aquatic firefly endemic to Japan, inhabiting a wide area of the Japanese archipelago. The luminescence of fireflies is a scientifically interesting phenomenon, and many studies have evaluated this species in Japan. In this study, we sequenced the whole genome of male N. cruciata and constructed a high-quality genome assembly of 662 Mb with a BUSCO completeness of 99.1% in the genome mode. Using the detected set of 15,169 protein-coding genes, the genomic structures and genetic background of luminescence-related genes were also investigated. We found four new firefly luciferase-like genes in the genome. The highest bioluminescent activity was observed for LLa2, which originated from ancestral PDGY, a mitochondrial acyl-CoA synthetase. A thioesterase candidate, NcruACOT1, which is involved in d-luciferin biosynthesis, was expressed in the lantern. Two opsins were also detected and the absorption wavelength of the UV-type opsin candidate shifted from UV to blue. These findings provide an important resource for unravelling the adaptive evolution of fireflies in terms of luminescence and vision.

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.

Key homeobox transcription factors regulate the development of the firefly's adult light organ and bioluminescence

Adult fireflies exhibit unique flashing courtship signals, emitted by specialized light organs, which develop mostly independently from larval light organs during the pupal stage. The mechanisms of adult light organ development have not been thoroughly studied until now. Here we show that key homeobox transcription factors AlABD-B and AlUNC-4 regulate the development of adult light organs and bioluminescence in the firefly Aquatica leii. Interference with the expression of AlAbd-B and AlUnc-4 genes results in undeveloped or non-luminescent adult light organs. AlABD-B regulates AlUnc-4, and they interact with each other. AlABD-B and AlUNC-4 activate the expression of the luciferase gene AlLuc1 and some peroxins. Four peroxins are involved in the import of AlLUC1 into peroxisomes. Our study provides key insights into the development of adult light organs and flash signal control in fireflies.

Genome evolution and divergence in cis-regulatory architecture is associated with condition-responsive development in horned dung beetles

Phenotypic plasticity is thought to be an important driver of diversification and adaptation to environmental variation, yet the genomic mechanisms mediating plastic trait development and evolution remain poorly understood. The Scarabaeinae, or true dung beetles, are a species-rich clade of insects recognized for their highly diversified nutrition-responsive development including that of cephalic horns-evolutionarily novel, secondary sexual weapons that exhibit remarkable intra- and interspecific variation. Here, we investigate the evolutionary basis for horns as well as other key dung beetle traits via comparative genomic and developmental assays. We begin by presenting chromosome-level genome assemblies of three dung beetle species in the tribe Onthophagini (> 2500 extant species) including Onthophagus taurus, O. sagittarius, and Digitonthophagus gazella. Comparing these assemblies to those of seven other species across the order Coleoptera identifies evolutionary changes in coding sequence associated with metabolic regulation of plasticity and metamorphosis. We then contrast chromatin accessibility in developing head horn tissues of high- and low-nutrition O. taurus males and females and identify distinct cis-regulatory architectures underlying nutrition- compared to sex-responsive development, including a large proportion of recently evolved regulatory elements sensitive to horn morph determination. Binding motifs of known and new candidate transcription factors are enriched in these nutrition-responsive open chromatin regions. Our work highlights the importance of chromatin state regulation in mediating the development and evolution of plastic traits, demonstrates gene networks are highly evolvable transducers of environmental and genetic signals, and provides new reference-quality genomes for three species that will bolster future developmental, ecological, and evolutionary studies of this insect group.

A Comprehensive Exploration of Bioluminescence Systems, Mechanisms, and Advanced Assays for Versatile Applications

Bioluminescence has been a fascinating natural phenomenon of light emission from living creatures. It happens when the enzyme luciferase facilitates the oxidation of luciferin, resulting in the creation of an excited-state species that emits light. Although there are many bioluminescent systems, few have been identified. D-luciferin-dependent systems, coelenterazine-dependent systems, Cypridina luciferin-based systems, tetrapyrrole-based luciferins, bacterial bioluminescent systems, and fungal bioluminescent systems are natural bioluminescent systems. Since different bioluminescence systems, such as various combinations of luciferin-luciferase pair reactions, have different light emission wavelengths, they benefit industrial applications such as drug discovery, protein-protein interactions, in vivo imaging in small animals, and controlling neurons. Due to the expression of luciferase and easy permeation of luciferin into most cells and tissues, bioluminescence assays are applied nowadays with modern technologies in most cell and tissue types. It is a versatile technique in a variety of biomedical research. Furthermore, there are some investigated blue-sky research projects, such as bioluminescent plants and lamps. This review article is mainly based on the theory of diverse bioluminescence systems and their past, present, and future applications.

An Improved Chromosome-Level Genome Assembly of the Firefly Pyrocoelia pectoralis

The endemic and endangered Chinese firefly Pyrocoelia pectoralis is a sexually dimorphic, nocturnal species. A previous attempt by this team to assemble a draft genome of P. pectoralis using PacBio and Illumina HiSeq X Ten platforms was limited in its usefulness by high redundancy and contamination. This prompted us to conduct an improved chromosome-level genome assembly of P. pectoralis. Ten chromosomes were further assembled based on Hi-C data to a 532.25 Mb final size with a 52.87 Mb scaffold N50. The total repeat lengths in the genome of P. pectoralis amount to 227.69 Mb; 42.78%. In total, 12,789 genes could be functionally annotated using at least one public database. Phylogenetic inference indicated that P. pectoralis and P. pyralis diverged ~51.41 million years ago. Gene family expansion and contraction analysis of 12 species were performed, and 546 expanded and 2660 contracted gene families were identified in P. pectoralis. We generated a high-quality draft of the P. pectoralis genome. This genome assembly should help promote research on the species' sexual dimorphism and its unique courtship behavior, which involves a combination of pheromonal and bioluminescent signals. It also can serve as a resource for accelerating genome-assisted improvements in the conservation of this species.

Multiple Origins of Bioluminescence in Beetles and Evolution of Luciferase Function

Bioluminescence in beetles has long fascinated biologists, with diverse applications in biotechnology. To date, however, our understanding of its evolutionary origin and functional variation mechanisms remains poor. To address these questions, we obtained high-quality reference genomes of luminous and nonluminous beetles in 6 Elateroidea families. We then reconstructed a robust phylogenetic relationship for all luminous families and related nonluminous families. Comparative genomic analyses and biochemical functional experiments suggested that gene evolution within Elateroidea played a crucial role in the origin of bioluminescence, with multiple parallel origins observed in the luminous beetle families. While most luciferase-like proteins exhibited a conserved nonluminous amino acid pattern (TLA346 to 348) in the luciferin-binding sites, luciferases in the different luminous beetle families showed divergent luminous patterns at these sites (TSA/CCA/CSA/LVA). Comparisons of the structural and enzymatic properties of ancestral, extant, and site-directed mutant luciferases further reinforced the important role of these sites in the trade-off between acyl-CoA synthetase and luciferase activities. Furthermore, the evolution of bioluminescent color demonstrated a tendency toward hypsochromic shifts and variations among the luminous families. Taken together, our results revealed multiple parallel origins of bioluminescence and functional divergence within the beetle bioluminescent system.

Predatory fireflies and their toxic firefly prey have evolved distinct toxin resistance strategies

Toxic cardiotonic steroids (CTSs) act as a defense mechanism in many firefly species (Lampyridae) by inhibiting a crucial enzyme called Na+,K+-ATPase (NKA). Although most fireflies produce these toxins internally, species of the genus Photuris acquire them from a surprising source: predation on other fireflies. The contrasting physiology of toxin exposure and sequestration between Photuris and other firefly genera suggests that distinct strategies may be required to prevent self-intoxication. Our study demonstrates that both Photuris and their firefly prey have evolved highly resistant NKAs. Using an evolutionary analysis of the specific target of CTS (ATPα) in fireflies and gene editing in Drosophila, we find that the initial steps toward resistance were shared among Photuris and other firefly lineages. However, the Photuris lineage subsequently underwent multiple rounds of gene duplication and neofunctionalization, resulting in the development of ATPα paralogs that are differentially expressed and exhibit increasing resistance to CTS. By contrast, other firefly species have maintained a single copy. Our results implicate gene duplication as a facilitator in the transition of Photuris to its distinct ecological role as a predator of toxic firefly prey.

Predatory fireflies and their toxic firefly prey have evolved distinct toxin resistance strategies

Toxic cardiotonic steroids (CTS) act as a defense mechanism in many firefly species (Lampyridae) by inhibiting a crucial enzyme called Na+,K+-ATPase (NKA). While most fireflies produce these toxins internally, species of the genus Photuris acquire them from a surprising source: predation on other fireflies. The contrasting physiology of toxin exposure and sequestration between Photuris and other firefly genera suggests that distinct strategies may be required to prevent self-intoxication. Our study demonstrates that both Photuris and their firefly prey have evolved highly-resistant NKAs. Using an evolutionary analysis of the specific target of CTS (ATPα) in fireflies, and gene-editing in Drosophila, we find that the initial steps towards resistance were shared among Photuris and other firefly lineages. However, the Photuris lineage subsequently underwent multiple rounds of gene duplication and neofunctionalization, resulting in the development of ATPα paralogs that are differentially expressed and exhibit increasing resistance to CTS. In contrast, other firefly species have maintained a single copy. Our results implicate gene duplication as a facilitator in the transition of Photuris to its distinct ecological role as predator of toxic firefly prey.

Firefly luciferin methyl ester illuminates the activity of multiple serine hydrolases

Firefly luciferin methyl ester is hydrolyzed by monoacylglycerol lipase MAGL, amidase FAAH, poorly-characterized hydrolase ABHD11, and hydrolases known for S-depalmitoylation (LYPLA1/2), not just esterase CES1. This enables activity-based bioluminescent assays for serine hydrolases and suggests that the 'esterase activity' responsible for hydrolyzing ester prodrugs is more diverse than previously supposed.

A chromosome-scale genome assembly and evaluation of mtDNA variation in the willow leaf beetle Chrysomela aeneicollis

The leaf beetle Chrysomela aeneicollis has a broad geographic range across Western North America but is restricted to cool habitats at high elevations along the west coast. Central California populations occur only at high altitudes (2,700-3,500 m) where they are limited by reduced oxygen supply and recent drought conditions that are associated with climate change. Here, we report a chromosome-scale genome assembly alongside a complete mitochondrial genome and characterize differences among mitochondrial genomes along a latitudinal gradient over which beetles show substantial population structure and adaptation to fluctuating temperatures. Our scaffolded genome assembly consists of 21 linkage groups; one of which we identified as the X chromosome based on female/male whole genome sequencing coverage and orthology with Tribolium castaneum. We identified repetitive sequences in the genome and found them to be broadly distributed across all linkage groups. Using a reference transcriptome, we annotated a total of 12,586 protein-coding genes. We also describe differences in putative secondary structures of mitochondrial RNA molecules, which may generate functional differences important in adaptation to harsh abiotic conditions. We document substitutions at mitochondrial tRNA molecules and substitutions and insertions in the 16S rRNA region that could affect intermolecular interactions with products from the nuclear genome. This first chromosome-level reference genome will enable genomic research in this important model organism for understanding the biological impacts of climate change on montane insects.

FruitFire: a luciferase based on a fruit fly metabolic enzyme

Firefly luciferase is homologous to fatty acyl-CoA synthetases from insects that are not bioluminescent. Here, we determined the crystal structure of the fruit fly fatty acyl-CoA synthetase CG6178 to 2.5 Å. Based on this structure, we mutated a steric protrusion in the active site to create the artificial luciferase FruitFire, which prefers the synthetic luciferin CycLuc2 to d-luciferin by >1000-fold. FruitFire enabled in vivo bioluminescence imaging in the brains of mice using the pro-luciferin CycLuc2-amide. The conversion of a fruit fly enzyme into a luciferase capable of in vivo imaging underscores the potential for bioluminescence with a range of adenylating enzymes from nonluminescent organisms, and the possibilities for application-focused design of enzyme-substrate pairs.

Pest status, molecular evolution, and epigenetic factors derived from the genome assembly of Frankliniella fusca, a thysanopteran phytovirus vector

BACKGROUND

The tobacco thrips (Frankliniella fusca Hinds; family Thripidae; order Thysanoptera) is an important pest that can transmit viruses such as the tomato spotted wilt orthotospovirus to numerous economically important agricultural row crops and vegetables. The structural and functional genomics within the order Thysanoptera has only begun to be explored. Within the > 7000 known thysanopteran species, the melon thrips (Thrips palmi Karny) and the western flower thrips (Frankliniella occidentalis Pergrande) are the only two thysanopteran species with assembled genomes.

RESULTS

A genome of F. fusca was assembled by long-read sequencing of DNA from an inbred line. The final assembly size was 370 Mb with a single copy ortholog completeness of ~ 99% with respect to Insecta. The annotated genome of F. fusca was compared with the genome of its congener, F. occidentalis. Results revealed many instances of lineage-specific differences in gene content. Analyses of sequence divergence between the two Frankliniella species' genomes revealed substitution patterns consistent with positive selection in ~ 5% of the protein-coding genes with 1:1 orthologs. Further, gene content related to its pest status, such as xenobiotic detoxification and response to an ambisense-tripartite RNA virus (orthotospovirus) infection was compared with F. occidentalis. Several F. fusca genes related to virus infection possessed signatures of positive selection. Estimation of CpG depletion, a mutational consequence of DNA methylation, revealed that F. fusca genes that were downregulated and alternatively spliced in response to virus infection were preferentially targeted by DNA methylation. As in many other insects, DNA methylation was enriched in exons in Frankliniella, but gene copies with homology to DNA methyltransferase 3 were numerous and fragmented. This phenomenon seems to be relatively unique to thrips among other insect groups.

CONCLUSIONS

The F. fusca genome assembly provides an important resource for comparative genomic analyses of thysanopterans. This genomic foundation allows for insights into molecular evolution, gene regulation, and loci important to agricultural pest status.

Flickering flash signals and mate recognition in the Asian firefly, Aquatica lateralis

Nocturnal fireflies sometimes use intricate bioluminescent signal systems for sexual communication. In this study, we examined flash signals and mate recognition in the Asian firefly, Aquatica lateralis, under natural field conditions. We found that the flash pattern of females changes after copulation, from simple short flashes to flashes with longer duration and flickering. To understand the functions of flickering, we video-recorded and analyzed the flashes of sedentary males, receptive females, and mated females. The results showed that the flashes of these three adult phases can be discriminated from each other by two parameters, flash duration and flicker intensity, with little overlap. Male attraction experiments using an artificial LED device termed 'e-firefly' confirmed that flying and sedentary males are attracted to flashes with shorter durations and lower flicker intensities. The range of attraction success was much wider for flying males and narrower for sedentary males, and the latter was close to the range of receptive female's flashes. These findings suggest that in addition to flash duration, flicker intensity is a flash signal parameter of mate recognition in A. lateralis males.

An Efficient Feature Selection Algorithm for Gene Families Using NMF and ReliefF

Gene families, which are parts of a genome's information storage hierarchy, play a significant role in the development and diversity of multicellular organisms. Several studies have focused on the characteristics of gene families, such as function, homology, or phenotype. However, statistical and correlation analyses on the distribution of gene family members in the genome have yet to be conducted. Here, a novel framework incorporating gene family analysis and genome selection based on NMF-ReliefF is reported. Specifically, the proposed method starts by obtaining gene families from the TreeFam database and determining the number of gene families within the feature matrix. Then, NMF-ReliefF is used to select features from the gene feature matrix, which is a new feature selection algorithm that overcomes the inefficiencies of traditional methods. Finally, a support vector machine is utilized to classify the acquired features. The results show that the framework achieved an accuracy of 89.1% and an AUC of 0.919 on the insect genome test set. We also employed four microarray gene data sets to evaluate the performance of the NMF-ReliefF algorithm. The outcomes show that the proposed method may strike a delicate balance between robustness and discrimination. Additionally, the proposed method's categorization is superior to state-of-the-art feature selection approaches.

Contribution of cryptochromes and photolyases for insect life under sunlight

The cryptochrome/photolyase (CRY/PL) family is essential for life under sunlight because photolyases repair UV-damaged DNA and cryptochromes are normally part of the circadian clock that controls the activity-sleep cycle within the 24-h day. In this study, we aim to understand how the lineage and habitat of an insect affects its CRY/PL composition. To this end, we searched the large number of annotated protein sequences of 340 insect species already available in databases for CRY/PLs. Using phylogenetic tree and motif analyses, we identified four frequent CRY/PLs in insects: the photolyases 6-4 PL and CPDII PL, as well as the mammalian-type cryptochrome (MCRY) and Drosophila-type cryptochrome (DCRY). Assignment of CRY/PLs to the corresponding insects confirmed that light-exposed insects tend to have more CRY/PLs than insects with little light exposure. Nevertheless, even insects with greatly reduced CRY/PLs still possess MCRY, which can be regarded as the major insect cryptochrome. Only flies of the genus Schizophora, which includes Drosophila melanogaster, lost MCRY. Moreover, we found that MCRY and CPDII PL as well as DCRY and 6-4 PL occur very frequently together, suggesting an interaction between the two pairs.

Rhagophthalmidae Olivier, 1907 (Coleoptera, Elateroidea): described genera and species, current problems, and prospects for the bioluminescent and paedomorphic beetle lineage

Rhagophthalmidae are a small beetle family known from the eastern Palaearctic and Oriental realms. Rhagophthalmidae are closely related to railroad worms (Phengodidae) and fireflies (Lampyridae) with which they share highly modified paedomorphic females and the ability to emit light. Currently, Rhagophthalmidae include 66 species classified in the following 12 genera: Bicladodrilus Pic, 1921 (two spp.), Bicladum Pic, 1921 (two spp.), Dioptoma Pascoe, 1860 (two spp.), Diplocladon Gorham, 1883 (two spp.), Dodecatoma Westwood, 1849 (eight spp.), Falsophrixothrix Pic, 1937 (six spp.), Haplocladon Gorham, 1883 (two spp.), Menghuoius Kawashima, 2000 (three spp.), Mimoochotyra Pic, 1937 (one sp.), Monodrilus Pic, 1921 (two spp. in two subgenera), Pseudothilmanus Pic, 1918 (two spp.), and Rhagophthalmus Motschulsky, 1854 (34 spp.). The replacement name Haplocladongorhami Kundrata, nom. nov. is proposed for Diplocladonhasseltii Gorham, 1883b (described in subgenus Haplocladon) which is preoccupied by Diplocladonhasseltii Gorham, 1883a. The genus Reductodrilus Pic, 1943 is tentatively placed in Lampyridae: Ototretinae. Lectotypes are designated for Pseudothilmanusalatus Pic, 1918 and P.marginalis Pic, 1918. Interestingly, in the eastern part of their distribution, Rhagophthalmidae have remained within the boundaries of the Sunda Shelf and the Philippines demarcated by the Wallace Line, which separates the Oriental and Australasian realms. This study is intended to be a first step towards a comprehensive revision of the group on both genus and species levels. Additionally, critical problems and prospects for rhagophthalmid research are briefly discussed.

Inferring the demographic history of the North American firefly Photinus pyralis

The firefly Photinus pyralis inhabits a wide range of latitudinal and ecological niches, with populations living from temperate to tropical habitats. Despite its broad distribution, its demographic history is unknown. In this study, we modelled and inferred different demographic scenarios for North American populations of P. pyralis, which were collected from Texas to New Jersey. We used a combination of ABC techniques (for multi-population/colonization analyses) and likelihood inference (dadi, StairwayPlot2, PoMo) for single-population demographic inference, which proved useful with our RAD data. We uncovered that the most ancestral North American population lays in Texas, which further colonized the Central region of the US and more recently the North Eastern coast. Our study confidently rejects a demographic scenario where the North Eastern populations colonized more southern populations until reaching Texas. To estimate the age of divergence between of P. pyralis, which provides deeper insights into the history of the entire species, we assembled a multi-locus phylogenetic data covering the genus Photinus. We uncovered that the phylogenetic node leading to P. pyralis lies at the end of the Miocene. Importantly, modelling the demographic history of North American P. pyralis serves as a null model of nucleotide diversity patterns in a widespread native insect species, which will serve in future studies for the detection of adaptation events in this firefly species, as well as a comparison for future studies of other North American insect taxa.

Genome analysis of Phrixothrix hirtus (Phengodidae) railroad worm shows the expansion of odorant-binding gene families and positive selection on morphogenesis and sex determination genes

Among bioluminescent beetles of the Elateroidea superfamily, Phengodidae is the third largest family, with 244 bioluminescent species distributed only in the Americas, but is still the least studied from the phylogenetic and evolutionary points of view. The railroad worm Phrixothrix hirtus is an essential biological model and symbolic species due to its bicolor bioluminescence, being the only organism that produces true red light among bioluminescent terrestrial species. Here, we performed partial genome assembly of P. hirtus, combining short and long reads generated with Illumina sequencing, providing the first source of genomic information and a framework for comparative analyses of the bioluminescent system in Elateroidea. This is the largest genome described in the Elateroidea superfamily, with an estimated size of ∼3.4 Gb, displaying 32 % GC content, and 67 % transposable elements. Comparative genomic analyses showed a positive selection of genes and gene family expansion events of growths and morphogenesis gene products, which could be associated with the atypical anatomical development and morphogenesis found in paedomorphic females and underdeveloped males. We also observed gene family expansion among distinct odorant-binding receptors, which could be associated with the pheromone communication system typical of these beetles, and retrotransposable elements. Common genes putatively regulating bioluminescence production and control, including two luciferase genes corresponding to lateral lanterns green-emitting and head lanterns red-emitting luciferases with 7 exons and 6 introns, and genes potentially involved in luciferin biosynthesis were found, indicating that there are no clear differences about the presence or absence of gene families associated with bioluminescence in Elateroidea.

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.

Beetle bioluminescence outshines extant aerial predators

We understand very little about the timing and origins of bioluminescence, particularly as a predator avoidance strategy. Understanding the timing of its origins, however, can help elucidate the evolution of this ecologically important signal. Using fireflies, a prevalent bioluminescent group where bioluminescence primarily functions as aposematic and sexual signals, we explore the origins of this signal in the context of their potential predators. Divergence time estimations were performed using genomic-scale datasets providing a robust estimate for the origin of firefly bioluminescence as both a terrestrial and as an aerial signal. Our results recover the origin of terrestrial beetle bioluminescence at 141.17 (122.63-161.17) Ma and firefly aerial bioluminescence at 133.18 (117.86-152.47) Ma using a large dataset focused on Lampyridae; and terrestrial bioluminescence at 148.03 (130.12-166.80) Ma, with the age of aerial bioluminescence at 104.97 (99.00-120.90) Ma using a complementary Elateroidea dataset. These ages pre-date the origins of all known extant aerial predators (i.e. bats and birds) and support much older terrestrial predators (assassin bugs, frogs, ground beetles, lizards, snakes, hunting spiders and harvestmen) as the drivers of terrestrial bioluminescence in beetles. These ages also support the hypothesis that sexual signalling was probably the original function of this signal in aerial fireflies.

A high-quality genome of the dobsonfly Neoneuromus ignobilis reveals molecular convergences in aquatic insects

Neoneuromus ignobilis is an archaic holometabolous aquatic predatory insect. However, a lack of genomic resources hinders the use of whole genome sequencing to explore their genetic basis and molecular mechanisms for adaptive evolution. Here, we provided a high-contiguity, chromosome-level genome assembly of N. ignobilis using high coverage Nanopore and PacBio reads with the Hi-C technique. The final assembly is 480.67 MB in size, containing 12 telomere-ended pseudochromosomes with only 17 gaps. We compared 42 hexapod species genomes including six independent lineages comprising 11 aquatic insects, and found convergent expansions of long wavelength-sensitive and blue-sensitive opsins, thermal stress response TRP channels, and sulfotransferases in aquatic insects, which may be related to their aquatic adaptation. We also detected strong nonrandom signals of convergent amino acid substitutions in aquatic insects. Collectively, our comparative genomic analysis revealed the evidence of molecular convergences in aquatic insects during both gene family evolution and convergent amino acid substitutions.

Firefly genomes illuminate the evolution of beetle bioluminescent systems

Fireflies are one of the best-known bioluminescent organisms, and the reaction mechanism and ecological utility of bioluminescence have been well-studied. Genome assemblies of six species of bioluminescent beetles have recently been published. These studies have focused on the evolution of novelties; luciferase, and the biosynthesis of luciferin and defensive chemicals. For example, clustering of the luciferase gene with acyl-CoA synthetase genes on a chromosome in luminous beetle genomes suggests the involvement of tandem gene duplications and neofunctionalization during the evolution of beetle bioluminescence. Several candidate genes for critical roles in beetle bioluminescence have been identified, but their functional analyses are still ongoing. The establishment of a long-term mass-rearing system and strain will be the key for the post-genome research on bioluminescent beetles. Lastly, the application of contemporary chromosome-scale genome assembly techniques to luminous beetles will help resolve outstanding evolutionary questions, such as how many times bioluminescence evolved in this clade.

Genomic insight into the scale specialization of the biological control agent Novius pumilus (Weise, 1892)

BACKGROUND

Members of the genus Novius Mulsant, 1846 (= Rodolia Mulsant, 1850) (Coleoptera, Coccinellidae), play important roles in the biological control of cotton cushion scale pests, especially those belonging to Icerya. Since the best-known species, the vedalia beetle Novius cardinalis (Mulsant, 1850) was introduced into California from Australia, more than a century of successful use in classical biological control, some species of Novius have begun to exhibit some field adaptations to novel but related prey species. Despite their economic importance, relatively little is known about the underlying genetic adaptations associated with their feeding habits. Knowledge of the genome sequence of Novius is a major step towards further understanding its biology and potential applications in pest control.

RESULTS

We report the first high-quality genome sequence for Novius pumilus (Weise, 1892), a representative specialist of Novius. Computational Analysis of gene Family Evolution (CAFE) analysis showed that several orthogroups encoding chemosensors, digestive, and immunity-related enzymes were significantly expanded (P < 0.05) in N. pumilus compared to the published genomes of other four ladybirds. Furthermore, some of these orthogroups were under significant positive selection pressure (P < 0.05). Notably, transcriptome profiling demonstrated that many genes among the significantly expanded and positively selected orthogroups, as well as genes related to detoxification were differentially expressed, when N. pumilus feeding on the nature prey Icerya compared with the no feeding set. We speculate that these genes are vital in the Icerya adaptation of Novius species.

CONCLUSIONS

We report the first Novius genome thus far. In addition, we provide comprehensive transcriptomic resources for N. pumilus. The results from this study may be helpful for understanding the association of the evolution of genes related to chemosensing, digestion, detoxification and immunity with the prey adaptation of insect predators. This will provide a reference for future research and utilization of Novius in biological control programs. Moreover, understanding the possible molecular mechanisms of prey adaptation also inform mass rearing of N. pumilus and other Novius, which may benefit pest control.

Involvement of CYP2 and mitochondrial clan P450s of Helicoverpa armigera in xenobiotic metabolism

Insect CYP2 and mitochondrial clan P450s are relatively conserved genes encoding enzymes generally thought to be involved in biosynthesis or metabolism of endobiotics. However, emerging evidence argues they have potential roles in chemical defense as well, but their actual detoxification functions remain largely unknown. Here, we focused on the full complement of 8 CYP2 and 10 mitochondrial P450s in the generalist herbivore, Helicoverpa armigera. Their varied spatiotemporal expression profiles were analyzed and reflected their specific functions. For functional study of the mitochondrial clan P450s, the redox partners, adrenodoxin reductase (AdR) and adrenodoxin (Adx), were identified from genomes of eight insects and an efficient in vitro electron transfer system of mitochondrial P450 was established by co-expression with Adx and AdR of H. armigera. All CYP2 clan P450s and 8 mitochondrial P450s were successfully expressed in Sf9 cells and compared functionally. In vitro metabolism assays showed that two CYP2 clan P450s (CYP305B1 and CYP18A1) and CYP333B3 (mito clan) could epoxidize aldrin to dieldrin, while CYP305B1 and CYP339A1 (mito clan) have limited but significant hydroxylation capacities to esfenvalerate. CYP303A1 of the CYP2 clan exhibits high metabolic efficiency to 2-tridecanone. Screening the xenobiotic metabolism competence of CYP2 and mitochondrial clan P450s not only provides new insights on insect chemical defense but also can give indications on their physiological functions in H. armigera and other insects.

Species-Specific Flash Patterns Track the Nocturnal Behavior of Sympatric Taiwanese Fireflies

It is highly challenging to evaluate the species' content and behavior changes in wild fireflies, especially for a sympatric population. Here, the flash interval (FI) and flash duration (FD) of flying males from three sympatric species (Abscondita cerata, Luciola kagiana, and Luciola curtithorax) were investigated for their potentials in assessing species composition and nocturnal behaviors during the A. cerata mating season. Both FI and FD were quantified from the continuous flashes of adult fireflies (lasting 5-30 s) via spatiotemporal analyses of video recorded along the Genliao hiking trail in Taipei, Taiwan. Compared to FD patterns and flash colors, FI patterns exhibited the highest species specificity, making them a suitable reference for differentiating firefly species. Through the case study of a massive occurrence of A. cerata (21 April 2018), the species contents (~85% of the flying population) and active periods of a sympatric population comprising A. cerata and L. kagiana were successfully evaluated by FI pattern matching, as well as field specimen collections. Our study suggests that FI patterns may be a reliable species-specific luminous marker for monitoring the behavioral changes in a sympatric firefly population in the field, and has implication values for firefly conservation.

InsectBase 2.0: a comprehensive gene resource for insects

Insects are the largest group of animals on the planet and have a huge impact on human life by providing resources, transmitting diseases, and damaging agricultural crop production. Recently, a large amount of insect genome and gene data has been generated. A comprehensive database is highly desirable for managing, sharing, and mining these resources. Here, we present an updated database, InsectBase 2.0 (http://v2.insect-genome.com/), covering 815 insect genomes, 25 805 transcriptomes and >16 million genes, including 15 045 111 coding sequences, 3 436 022 3'UTRs, 4 345 664 5'UTRs, 112 162 miRNAs and 1 293 430 lncRNAs. In addition, we used an in-house standard pipeline to annotate 1 434 653 genes belonging to 164 gene families; 215 986 potential horizontally transferred genes; and 419 KEGG pathways. Web services such as BLAST, JBrowse2 and Synteny Viewer are provided for searching and visualization. InsectBase 2.0 serves as a valuable platform for entomologists and researchers in the related communities of animal evolution and invertebrate comparative genomics.

Acquisition of bioluminescent trait by non-luminous organisms from luminous organisms through various origins

Bioluminescence is a natural light emitting phenomenon that occurs due to a chemical reaction between luciferin and luciferase. It is primarily an innate and inherited trait in most terrestrial luminous organisms. However, most luminous organisms produce light in the ocean by acquiring luminous symbionts, luciferin (substrate), and/or luciferase (enzyme) through various transmission pathways. For instance, coelenterazine, a well-known luciferin, is obtained by cnidarians, crustaceans, and deep-sea fish through multi-level dietary linkages from coelenterazine producers such as ctenophores, decapods, and copepods. In contrast, some non-luminous Vibrio bacteria became bioluminescent by obtaining lux genes from luminous Vibrio species by horizontal gene transfer. Various examples detailed in this review show how non-luminescent organisms became luminescent by acquiring symbionts, dietary luciferins and luciferases, and genes. This review highlights three modes (symbiosis, ingestion, and horizontal gene transfer) that allow organisms lacking genes for autonomous bioluminescent systems to obtain the ability to produce light. In addition to bioluminescence, this manuscript discusses the acquisition of other traits such as pigments, fluorescence, toxins, and others, to infer the potential processes of acquisition.

Whole-Genome Sequence of the Trypoxylus dichotomus Japanese rhinoceros beetle

The draft whole-genome sequence of the Japanese rhinoceros beetle, Trypoxylus dichotomus was obtained using long-read PacBio sequence technology. The final assembled genome consisted of 739 Mbp in 2,347 contigs, with 24.5× mean coverage and a G+C content of 35.99%.

A comprehensive review and call for studies on firefly larvae

BACKGROUND

Fireflies (Coleoptera: Lampyridae) are commonly recognized by adult traits, such as a soft exoskeleton, lanterns and associated glow and flash patterns, but their larval stage is far less appreciated. However, fireflies spend most of their lives as larvae, and adults of most species rely solely on resources previously obtained. Therefore, studying the immature stages is imperative towards a comprehensive understanding of fireflies. This paper reviews and indicates key gaps in the biology of firefly larvae based on available literature.

METHODOLOGY

We reviewed the literature on firefly larvae to identify key issues and important taxonomic, geographic, and subject biases and gaps.

RESULTS

We found 376 papers that included information on firefly larvae. Only 139 species in 47 genera across eight of eleven lampyrid subfamilies have been studied during larval stages. These numbers reveal a staggering gap, since 94% of species and over half of the genera of fireflies were never studied in a crucial stage of their life cycle. Most studies on firefly larvae focus on two subfamilies (Luciolinae and Lampyrinae) in four zoogeographic regions (Sino-Japanese, Oriental, Nearctic, and Palearctic), whereas the other subfamilies and regions remain largely unstudied. These studies mainly dealt with morphology and behavior, other subjects remaining greatly understudied by comparison, including habitats, life cycle, physiology and interactions.

CONCLUSIONS

Together, these literature biases and gaps highlight how little is known about firefly larvae, and warmly invite basic and applied research, in the field and in the lab, to overcome these limitations and improve our understanding of firefly biology to better preserve them.

Molecular cloning, characterization, and evolution analysis of the luciferase genes from three sympatric sibling fireflies (Lampyridae: Lampyrinae, Diaphanes)

Firefly adult bioluminescence functions as signal communication between sexes. How sympatric sibling species with similar glow pattern recognize their conspecific mates remains largely unknown. To better understand the role of the luciferases of sympatric fireflies in recognizing mates, we cloned the luciferase genes of three sympatric forest dwelling fireflies (Diaphanes nubilus, Diaphanes pectinealis, and Diaphanes sp2) and evaluated their enzyme characteristics. Our data show that the amino acid (AA) sequences of all three luciferases are highly conserved, including the identities (D. nubilus vs D. pectinealis: 99%; D. nubilus vs Diaphanes sp2: 98.5%; D. pectinealis vs Diaphanes sp2: 99.4%) and the protein structures. Three recombinant luciferases produced in vitro all possess significant luminescence activity at pH 7.8, and similar maximum emission spectrum (D. nubilus: 562 nm; D. pectinealis and Diaphanes sp2: 564 nm). They show the highest activity at 10 °C (D. pectinealis, Diaphanes sp2) and 15 °C (D. nubilus), and completely inactivation at 45 °C. Their KM for D-luciferin and ATP were 2.7 μM and 92 μM (D. nubilus), 3.7 μM and 49 μM (D. pectinealis), 3.5 μM and 46 μM (Diaphanes sp2). Phylogenetic analyses support that D. nubilus is sister to D. pectinealis with Diaphanes sp2 at their base, which further cluster with Pyrocoelia. All combined data indicate that sympatric Diaphanes species have similar luciferase characteristics, suggesting that other strategies (e.g., pheromone, active time, etc.) may be adopted to recognize mates. Our data provide new insights into Diaphanes luciferases and their evolution.

The Easter Egg Weevil (Pachyrhynchus) genome reveals syntenic patterns in Coleoptera across 200 million years of evolution

Patterns of genomic architecture across insects remain largely undocumented or decoupled from a broader phylogenetic context. For instance, it is unknown whether translocation rates differ between insect orders. We address broad scale patterns of genome architecture across Insecta by examining synteny in a phylogenetic framework from open-source insect genomes. To accomplish this, we add a chromosome level genome to a crucial lineage, Coleoptera. Our assembly of the Pachyrhynchus sulphureomaculatus genome is the first chromosome scale genome for the hyperdiverse Phytophaga lineage and currently the largest insect genome assembled to this scale. The genome is significantly larger than those of other weevils, and this increase in size is caused by repetitive elements. Our results also indicate that, among beetles, there are instances of long-lasting (>200 Ma) localization of genes to a particular chromosome with few translocation events. While some chromosomes have a paucity of translocations, intra-chromosomal synteny was almost absent, with gene order thoroughly shuffled along a chromosome. This large amount of reshuffling within chromosomes with few inter-chromosomal events contrasts with patterns seen in mammals in which the chromosomes tend to exchange larger blocks of material more readily. To place our findings in an evolutionary context, we compared syntenic patterns across Insecta in a phylogenetic framework. For the first time, we find that synteny decays at an exponential rate relative to phylogenetic distance. Additionally, there are significant differences in decay rates between insect orders, this pattern was not driven by Lepidoptera alone which has a substantially different rate.

Theoretical Study on Storage and Release of Firefly Luciferin

Among numerous bioluminescent organisms, firefly is the most studied one. Recent experiment proposed that sulfoluciferin (SLH₂ ) may serve as a storage form of luciferin (LH₂ ). In the present article, we employed density functional theory calculation to uncover the mechanism and detailed process of the storage and release reactions. Due to lack of available crystallographic structure of the related enzyme, the calculation was performed on a model system. For the storage reaction, possible amino acid residues were used for imitating the protein environment. For the release reaction, the dielectric constant of 3.0 was employed to simulate the polarity of the protein cavity. The computational results indicated that the reactions from LH₂ to SLH₂ and from SLH₂ to LH₂ are both exergonic, which favor the storage and release processes and coincide with the experimental observation. Basing on experimental and current theoretical study, we supplemented the stages of LH₂ storage and release in the entire bioluminescent cycle of firefly. The current theoretical calculation could inspire the study on LH₂ storage and release of other bioluminescent organisms.

Insect Collections as an Untapped Source of Bioactive Compounds-Fireflies (Coleoptera: Lampyridae) and Cardiotonic Steroids as a Proof of Concept

Natural history collections provide an invaluable basis for systematics, ecology, and conservation. Besides being an important source of DNA, museum specimens may also contain a plethora of natural products. Especially, dried insect collections represent a global repository with billions of inventoried vouchers. Due to their vast diversity, insects possess a great variety of defensive compounds, which they either produce autogenously or derive from the environment. Here, we present a case study on fireflies (Coleoptera: Lampyridae), which produce bufadienolides as a defense against predators. These toxins belong to the cardiotonic steroids, which are used for the treatment of cardiac diseases and specifically inhibit the animal enzyme Na+/K+-ATPase. Bufadienolides have been reported from only seven out of approximately 2000 described firefly species. Using a non-destructive approach, we screened 72 dry coleopteran specimens for bufadienolides using HPLC-DAD and HPLC-MS. We found bufadienolides including five novel compounds in 21 species of the subfamily Lampyrinae. The absence of bufadienolides in the phylogenetically related net-winged beetles (Lycidae) and the lampyrid subfamilies Luciolinae and Lamprohizinae indicates a phylogenetic pattern of bufadienolide synthesis. Our results emphasize the value of natural history collections as an archive of chemical information for ecological and evolutionary basic research and as an untapped source for novel bioactive compounds.

Multiple Cypridina Luciferase Genes in the Genome of Individual Ostracods, Vargula hilgendorfii (Cypridina hilgendorfii)

The genomic structure of the Cypridina luciferase gene in Vargula hilgendorfii (formerly Cypridina hilgendorfii) was determined with three λ phage clones (λ34, λ45, and λ61). The luciferase genes in clones λ34 and λ61 consisted of 13 exons and 12 introns, and clone λ45 only contained exons 1-5. The splicing sites of the luciferase genes in λ34 and λ61 were conserved completely with the consensus sequence. The translated luciferases had 555 amino acid residues, which were over 98.6% identical to those of cDNA clones as previously reported. In contrast, each intron in clones λ34, λ45, and λ61 varied significantly in length. To explain the variation of intron length among the three V. hilgendorfii luciferase genes, genomic DNA was isolated from a single V. hilgendorfii specimen and the regions from exon 1-3 of the luciferase gene were amplified by polymerase chain reaction (PCR). PCR products with various lengths were detected and were confirmed as the luciferase gene fragments by Southern blot analysis. Furthermore, DNA sequence analysis indicated that at least seven luciferase gene groups might be present in the genome of a single specimen. Thus, multiple Cypridina luciferase genes exist in the genome of a single V. hilgendorfii specimen.

Recent advances in chemical ecology: complex interactions mediated by molecules

Chemical ecology is the highly interdisciplinary study of biochemicals that mediate the behavior of organisms and the regulation of physiological changes that alter intraspecific and/or interspecific interactions. Significant advances are often achieved through the collaboration of chemists and biologists working to understand organismal survival strategies with an eye on the development of targeted technologies for controlling agricultural, forestry, medical, and veterinary pests in a sustainable world. We highlight recent advances in chemical ecology from multiple viewpoints and discuss future prospects for applications.

Characterization of the First Complete Mitochondrial Genome of Cyphonocerinae (Coleoptera: Lampyridae) with Implications for Phylogeny and Evolution of Fireflies

Complete mitochondrial genomes are valuable resources for phylogenetics in insects. The Cyphonoceridae represents an important lineage of fireflies. However, no complete mitogenome is available until now. Here, the first complete mitochondrial genome from this subfamily was reported, with Cyphonocerus sanguineus klapperichi as a representative. The mitogenome of C. sanguineus klapperichi was conserved in the structure and comparable to that of others in size and A+T content. Nucleotide composition was A+T-biased, and all genes exhibited a positive AT-skew and negative GC-skew. Two types of tandem repeat sequence units were present in the control region (136 bp × 2; 171 bp × 2 + 9 bp). For reconstruction of Lampyridae's phylogeny, three different datasets were analyzed by both maximum likelihood (ML) and Bayesian inference (BI) methods. As a result, the same topology was produced by both ML analysis of 13 protein-coding genes and 2rRNA and BI analysis of 37 genes. The results indicated that Lampyridae, Lampyrinae, Luciolinae (excluding Emeia) were monophyletic, but Ototretinae was paraphyletic, of which Stenocladius was recovered as the sister taxon to all others, while Drilaster was more closely related to Cyphonocerinae; Phturinae + Emeia were included in a monophyletic clade, which comprised sister groups with Lampyridae. Vesta was deeply rooted in the Luciolinae.

Molecular Evolution of Phototransduction Pathway Genes in Nocturnal and Diurnal Fireflies (Coleoptera: Lampyridae)

Most organisms are dependent on sensory cues from their environment for survival and reproduction. Fireflies (Coleoptera: Lampyridae) represent an ideal system for studying sensory niche adaptation due to many species relying on bioluminescent communication; as well as a diversity of ecologies. Here; using transcriptomics; we examine the phototransduction pathway in this non-model organism; and provide some of the first evidence for positive selection in the phototransduction pathway beyond opsins in beetles. Evidence for gene duplications within Lampyridae are found in inactivation no afterpotential C and inactivation no&nbsp;afterpotential&nbsp;D. We also find strong support for positive selection in arrestin-2; inactivation no afterpotential D; and transient receptor potential-like; with weak support for positive selection in guanine nucleotide-binding protein G(q) subunit alpha and neither inactivation nor afterpotential C. Taken with other recent work in flies; butterflies; and moths; this represents an exciting new avenue of study as we seek to further understand diversification and constraint on the phototransduction pathway in light of organism ecology.

Light Polarization by Biological Nanocoatings

Light plays paramount functions for living beings in nature. In addition to color, the polarization of light is used by many animals for navigation and communication. In this study, we describe the light polarizing role of special nanostructures coating cuticular surfaces of diverse arthropods. These structures are built as parallel nanoscale ridges covering the eyes of the sunlight-navigating spider Drassodes lapidosus and of the water pond-swarming black fly Simulium vittatum, as well as the light-emitting abdominal lantern of the firefly Aquatica lateralis. Exact topography and dimensions of the parallel nanoridges provide different light polarizing efficiencies and wavelength sensitivity. Optical modeling confirms that the nanoscale ridges are responsible for the spectral polarization dependency. Co-opting from our recent work on the self-assembly of Drosophila corneal nanostructures, we engineer arthropod-like parallel nanoridges on artificial surfaces, which recapitulate the light polarization effects. Our work highlights the fundamental importance of nanocoatings in arthropods for the light polarization management and provides a new biomimetic approach to produce ordered nanostructures under mild conditions.

Anchored Phylogenomics, Evolution and Systematics of Elateridae: Are All Bioluminescent Elateroidea Derived Click Beetles?

Simple Summary

In the era of phylogenomics, new molecular sequencing and computational techniques can aid in resolving phylogenetic relationships that were previously intractable by morphological or limited molecular data. In this study, we used anchored hybrid enrichment—designed to recover DNA sequences from hundreds of single-copy orthologous genes—to resolve the phylogeny of the Elateridae (click-beetles) and establish their placement within superfamily Elateroidea. The resulting data were compatible with published transcriptomes, allowing for integrating our dataset with previously published data. Using a wide range of analyses on these molecular data, we tested hypotheses long-debated in the morphological literature and also the robustness of our phylogenetic inferences. Our results placed the bioluminescent lampyroids (fireflies and relatives) within the click-beetles, challenging the current classification of Elateridae, Lampyridae, Phengodidae, and Rhagophthalmidae. However, despite the large amount of molecular data analyzed, a few nodes with conflicting phylogenetic signals could not be unambiguously resolved. Overall, we recovered well-resolved tree topologies that will serve as a framework for further systematic and evolutionary studies of click-beetles. This work further demonstrates that the click-beetle lineage contains not only pest wireworms, but also many species that benefit agriculture.

Abstract

Click-beetles (Coleoptera: Elateridae) are an abundant, diverse, and economically important beetle family that includes bioluminescent species. To date, molecular phylogenies have sampled relatively few taxa and genes, incompletely resolving subfamily level relationships. We present a novel probe set for anchored hybrid enrichment of 2260 single-copy orthologous genes in Elateroidea. Using these probes, we undertook the largest phylogenomic study of Elateroidea to date (99 Elateroidea, including 86 Elateridae, plus 5 non-elateroid outgroups). We sequenced specimens from 88 taxa to test the monophyly of families, subfamilies and tribes. Maximum likelihood and coalescent phylogenetic analyses produced well-resolved topologies. Notably, the included non-elaterid bioluminescent families (Lampyridae + Phengodidae + Rhagophthalmidae) form a clade within the otherwise monophyletic Elateridae, and Sinopyrophoridae may not warrant recognition as a family. All analyses recovered the elaterid subfamilies Elaterinae, Agrypninae, Cardiophorinae, Negastriinae, Pityobiinae, and Tetralobinae as monophyletic. Our results were conflicting on whether the hypnoidines are sister to Dendrometrinae or Cardiophorinae + Negastriinae. Moreover, we show that fossils with the eucnemid-type frons and elongate cylindrical shape may belong to Eucnemidae, Elateridae: Thylacosterninae, ancestral hard-bodied cantharoids or related extinct groups. Proposed taxonomic changes include recognition of Plastocerini as a tribe in Dendrometrinae and Hypnoidinae stat. nov. as a subfamily within Elateridae.

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

The hyperdiverse order Coleoptera comprises a staggering ~25% of known species on Earth. Despite recent breakthroughs in next generation sequencing, there remains a limited representation of beetle diversity in assembled genomes. Most notably, the ground beetle family Carabidae, comprising more than 40,000 described species, has not been studied in a comparative genomics framework using whole genome data. Here we generate a high-quality genome assembly for Nebria riversi, to examine sources of novelty in the genome evolution of beetles, as well as genetic changes associated with specialization to high-elevation alpine habitats. In particular, this genome resource provides a foundation for expanding comparative molecular research into mechanisms of insect cold adaptation. Comparison to other beetles shows a strong signature of genome compaction, with N. riversi possessing a relatively small genome (~147 Mb) compared to other beetles, with associated reductions in repeat element content and intron length. Small genome size is not, however, associated with fewer protein-coding genes, and an analysis of gene family diversity shows significant expansions of genes associated with cellular membranes and membrane transport, as well as protein phosphorylation and muscle filament structure. Finally, our genomic analyses show that these high-elevation beetles have endosymbiotic Spiroplasma, with several metabolic pathways (e.g., propanoate biosynthesis) that might complement N. riversi, although its role as a beneficial symbiont or as a reproductive parasite remains equivocal.

Selection and Validation of Suitable Reference Genes for RT-qPCR Analysis in the Rare Aquatic Firefly Aquatica leii (Coleoptera: Lampyridae)

Simple Summary

The purpose of this paper was to determine which genes in Aquatica leii, an endemic Chinese firefly, are suitable for use as reference genes in future RT-qPCR studies. This is the first example of an effective knockdown in fireflies and an important methodological advance for continuing research in A. leii. The results of this study will help improve accuracy and reliability to normalize RT-qPCR data in A. leii for further molecular analysis. Since luciferase is very abundant in the adult light organ of fireflies generally, this is an encouraging result and of interest not just for A. leii researchers but other researchers who plan to perform RNAi or RT-qPCR in fireflies, or other nonmodel insects, as well.

Abstract

Aquatica leii Fu and Ballantyne is a species of rare aquatic firefly and endemic in China. It is considered good material to study the molecular mechanism of sexual flash communication systems. To improve conservation and behavioral research strategies, large-scale genetic studies involving gene-expression analysis are required and reverse transcription-quantitative polymerase chain reaction (RT-qPCR) is the most commonly used method. However, there have been very few reports on appropriate reference genes in any species of firefly. Here, we evaluated eight widely utilized reference genes including 18S, Actin, Reep5, Odc1, Tub, Gapdh, Ef1a and S27Ae for their expression stabilities in A. leii under three different conditions, i.e., life stage, tissue and dsRNA injection. Based on the gene stability ranking calculated by RefFinder, which integrates four algorithms (geNorm, delta Ct method, NormFinder, and BestKeeper), we recommend S27Ae and Reep5 as the most appropriate reference genes for molecular studies in different life stages; Ef1a and Odc1 for different tissues; Tub and Odc1 for RNAi studies. The most appropriate reference genes in all treatments are S27Ae and Tub. The results of this study will help improve accuracy and reliability to normalize RT-qPCR data in A. leii for further molecular analysis.