Genomic and experimental data provide new insights into luciferin biosynthesis and bioluminescence evolution in fireflies

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.

Cloning and characterization of luciferase from an Asian firefly Pygoluciola qingyu and its comparison with other beetle luciferases

The bioluminescence system of luminescent beetles has extensive applications in biological imaging, protein labeling and drug screening. To explore wild luciferases with excellent catalytic activity and thermal stability, we cloned the luciferase of Pygoluciola qingyu, one species living in areas of high temperature and with strong bioluminescence, by combining transcriptomic sequencing and reverse transcription polymerase chain reaction (RT-PCR). The total length of luciferase gene is 1638 bp and the luciferase consists 544 amino acids. The recombinant P. qingyu luciferase was produced in vitro and its characteristics were compared with those of eight luciferases from China firefly species and two commercial luciferases. Compared with these luciferases, the P. qingyu luciferase shows the highest luminescence activity at room temperature (about 25-28 ℃) with similar KM value for D-luciferin and ATP to the Photinus pyralis luciferase. The P. qingyu luciferase activity was highest at 35 ℃ and can keep high activity at 30-40 ℃, which suggests the potential of P. qingyu luciferase for in vivo and cell application. Our results provide new insights into P. qingyu luciferase and give a new resource for the application of luciferases.

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.

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.

Scalemic natural products

Covering: up to the end of 2022The area of scalemic natural products is often enigmatic from a mechanistic standpoint, since low optical purity is observed in compounds having multiple contiguous stereogenic centers resulting from mechanistically distinct biogenetic steps. A scalemic state is rarely the result of a sloppy enzymatic activity, rather resulting from the expression of antipodal enzymes/directing proteins or from the erosion of optical purity by enzymatic or spontaneous reactions. Evidence for these processes is critically reviewed, identifying the mechanisms most often associated to the enzymatic generation of scalemic natural products and also discussing analytical exploitations of natural products' scalemicity.

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.

Luminescent characteristics and mitochondrial COI barcodes of nine cohabitated Taiwanese fireflies

Background

Over 50 Taiwanese firefly species have been discovered, but scientists lack information regarding most of their genetics, bioluminescent features, and cohabitating phenomena. In this study, we focus on morphological species identification and phylogeny reconstructed by COI barcoding, as well as luminescent characteristics of cohabited Taiwanese firefly species to determine the key factors that influenced how distinct bioluminescent species evolved to coexist and proliferate within the same habitat.

Methods

In this study, 366 specimens from nine species were collected in northern Taiwan from April to August, 2016-2019. First, the species and sex of the specimens were morphologically and genetically identified. Then, their luminescent spectra and intensities were recorded using a spectrometer and a power meter, respectively. The habitat temperature, relative humidity, and environmental light intensity were also measured. The cytochrome oxidase I (COI) gene sequence was used as a DNA barcode to reveal the phylogenetic relationships of cohabitated species.

Results

Nine species-eight adult species (Abscondita chinensis, Abscondita cerata, Aquatica ficta, Luciola curtithorax, Luciola kagiana, Luciola filiformis, Curtos sauteri, and Curtos costipennis) and one larval Pyrocoelia praetexta-were morphologically identified. The nine species could be found in April-August. Six of the eight adult species shared an overlap occurrence period in May. Luminescent spectra analysis revealed that the λ _max of studied species ranged from 552-572 nm (yellow-green to orange-yellow). The average luminescent intensity range of these species was about 1.2-14 lux (182.1-2,048 nW/cm²) for males and 0.8-5.8 lux (122.8-850 nW/cm²) for females, and the maximum luminescent intensity of males was 1.01-7.26-fold higher than that of females. Compared with previous studies, this study demonstrates that different λ _max, species-specific flash patterns, microhabitat choices, nocturnal activity time, and/or an isolated mating season are key factors that may lead to the species-specific courtship of cohabitated fireflies. Moreover, we estimated that the fireflies start flashing or flying when the environmental light intensity decreased to 6.49-28.1 lux. Thus, based on a rough theoretical calculation, the sensing distance between male and female fireflies might be 1.8-2.7 m apart in the dark. In addition, the mitochondrial COI barcode identified species with high resolution and suggested that most of the studied species have been placed correctly with congeners in previous phylogenies. Several cryptic species were revealed by the COI barcode with 3.27%-12.3% variation. This study renews the idea that fireflies' luminescence color originated from the green color of a Lampyridae ancestor, then red-shifted to yellow-green in Luciolinae, and further changed to orange-yellow color in some derived species.

Biosensing firefly luciferin synthesis in bacteria reveals a cysteine-dependent quinone detoxification route in Coleoptera

Luciferin biosynthetic origin and alternative biological functions during the evolution of beetles remain unknown. We have set up a bioluminescent sensing method for luciferin synthesis from cysteine and benzoquinone using E. coli and Pichia pastoris expressing the bright Amydetes vivianii firefly and P. termitilluminans click beetle luciferases. In the presence of d-cysteine and benzoquinone, intense bioluminescence is quickly produced, indicating the expected formation of d-luciferin. Starting with l-cysteine and benzoquinone, the bioluminescence is weaker and delayed, indicating that bacteria produce l-luciferin, and then racemize it to d-luciferin in the presence of endogenous esterases, CoA and luciferase. In bacteria the p-benzoquinone toxicity (I_C50 ~ 25 µM) is considerably reduced in the presence of cysteine, maintaining cell viability at 3.6 mM p-benzoquinone concomitantly with the formation of luciferin. Transcriptional analysis showed the presence of gene products involved with the sclerotization/tanning in the photogenic tissues, suggesting a possible link between these pathways and bioluminescence. The lack of two enzymes involved with the last steps of these pathways, indicate the possible accumulation of toxic quinone intermediates in the lanterns. These results and the abundance of cysteine producing enzymes suggest that luciferin first appeared as a detoxification byproduct of cysteine reaction with accumulated toxic quinone intermediates during the evolution of sclerotization/tanning in Coleoptera.

Integrated Analysis of Transcriptome and Proteome to Reveal Pupal Color Switch in Papilio xuthus Butterflies

Pupal color polyphenism in Papilio butterflies, including green, intermediate, or brown, is an excellent study system for understanding phenotypic plasticity. Previous studies suggested that development of brown pupae may be controlled by a hormone called pupal-cuticle-melanizing-hormone (PCMH) which is synthesized and secreted from brain-suboesophageal ganglion and prothoracic ganglion complexes (Br-SG-TG1) during the pre-pupa stage. However, detailed molecular mechanisms of neuroendocrine regulation in pupal color development remain unknown. In this study, we integrated the expression profiles of transcriptome and proteome at pre-pupa stages [2 h after gut purge (T1) and 3 h after forming the garter around the body (T2)] and pigmentation stages [10 h after ecdysis (T3) and 24 h after ecdysis (T4)] to identify important genes and pathways underlying the development of green and brown pupa in the swallowtail butterfly Papilio xuthus. Combined comparisons of each developmental stage and each tissue under green and brown conditions, a total of 1042 differentially expressed genes (DEGs) and 430 different abundance proteins (DAPs) were identified. Weighted gene co-expression network analysis (WGCNA) and enrichment analysis indicate that these DEGs were mainly related to oxidation-reduction, structural constituent of cuticle, and pigment binding. Soft clustering by Mfuzz and enrichment analysis indicate that these DAPs are mainly involved in tyrosine metabolism, insect hormone biosynthesis, and melanogenesis. By homologous alignment, we further identified those genes encoding neuropeptides (51), GPCRs (116), G-proteins (8), cuticular proteins (226), chitinases (16), and chitin deacetylases (8) in the whole genome of P. xuthus and analyzed their expression profiles. Although we identified no gene satisfying with hypothesized expression profile of PCMH, we found some genes in the neuropeptide cascade showed differentially expressed under two pupal color conditions. We also found that Toll signaling pathway genes, juvenile hormone (JH) related genes, and multiple cuticular proteins play important roles in the formation of selective pupal colors during the prepupal-pupal transition. Our data also suggest that both green and brown pupa include complex pigment system that is regulated by genes involved in black, blue, and yellow pigments. Our results provide important insights into the evolution of pupal protective colors among swallowtail butterflies.

Luciferins Under Construction: A Review of Known Biosynthetic Pathways

Bioluminescence, or the ability of a living organism to generate visible light, occurs as a result of biochemical reaction where enzyme, known as a luciferase, catalyzes the oxidation of a small-molecule substrate, known as luciferin. This advantageous trait has independently evolved dozens of times, with current estimates ranging from the most conservative 40, based on the biochemical diversity found across bioluminescence systems (Haddock et al., 2010) to 100, taking into account the physiological mechanisms involved in the behavioral control of light production across a wide range of taxa (Davis et al., 2016; Verdes and Gruber, 2017; Bessho-Uehara et al., 2020a; Lau and Oakley, 2021). Chemical structures of ten biochemically unrelated luciferins and several luciferase gene families have been described; however, a full biochemical pathway leading to light emission has been elucidated only for two: bacterial and fungal bioluminescence systems. Although the recent years have been marked by extraordinary discoveries and promising breakthroughs in understanding the molecular basis of multiple bioluminescence systems, the mechanisms of luciferin biosynthesis for many organisms remain almost entirely unknown. This article seeks to provide a succinct overview of currently known luciferins’ biosynthetic pathways.

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 afterpotential 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.

Genome-wide identification and gene-editing of pigment transporter genes in the swallowtail butterfly Papilio xuthus

Background

Insect body coloration often functions as camouflage to survive from predators or mate selection. Transportation of pigment precursors or related metabolites from cytoplasm to subcellular pigment granules is one of the key steps in insect pigmentation and usually executed via such transporter proteins as the ATP-binding cassette (ABC) transmembrane transporters and small G-proteins (e.g. Rab protein). However, little is known about the copy numbers of pigment transporter genes in the butterfly genomes and about the roles of pigment transporters in the development of swallowtail butterflies.

Results

Here, we have identified 56 ABC transporters and 58 Rab members in the genome of swallowtail butterfly Papilio xuthus. This is the first case of genome-wide gene copy number identification of ABC transporters in swallowtail butterflies and Rab family in lepidopteran insects. Aiming to investigate the contribution of the five genes which are orthologous to well-studied pigment transporters (ABCG: white, scarlet, brown and ok; Rab: lightoid) of fruit fly or silkworm during the development of swallowtail butterflies, we performed CRISPR/Cas9 gene-editing of these genes using P. xuthus as a model and sequenced the transcriptomes of their morphological mutants. Our results indicate that the disruption of each gene produced mutated phenotypes in the colors of larvae (cuticle, testis) and/or adult eyes in G0 individuals but have no effect on wing color. The transcriptomic data demonstrated that mutations induced by CRISPR/Cas9 can lead to the accumulation of abnormal transcripts and the decrease or dosage compensation of normal transcripts at gene expression level. Comparative transcriptomes revealed 606 ~ 772 differentially expressed genes (DEGs) in the mutants of four ABCG transporters and 1443 DEGs in the mutants of lightoid. GO and KEGG enrichment analysis showed that DEGs in ABCG transporter mutants enriched to the oxidoreductase activity, heme binding, iron ion binding process possibly related to the color display, and DEGs in lightoid mutants are enriched in glycoprotein binding and protein kinases.

Conclusions

Our data indicated these transporter proteins play an important role in body color of P. xuthus. Our study provides new insights into the function of ABC transporters and small G-proteins in the morphological development of butterflies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07400-z.