Repeated inversions within a pannier intron drive diversification of intraspecific colour patterns of ladybird beetles

Dissecting roles of pannier splice variants during pupal and adult morphogenesis in Henosepilachna vigintioctopunctata

The GATA transcription factor gene, pannier (pnr), has been extensively studied in Drosophila, revealing its crucial role in dorsal closure, heart development and the regulation of cuticular bristle patterns in adults. However, studies on the functions of pnr in the development of coleopteran insects are still scarce. Herein, we identified the pnr gene in Henosepilachna vigintioctopunctata and discovered two splicing variants named Hvpnr-α and Hvpnr-β respectively. Temporal expression analysis revealed that Hvpnr, Hvpnr-α and Hvpnr-β were expressed at various stages including egg, larval, pupal and adult stages. To investigate the developmental role of Hvpnr in H. vigintioctopunctata, RNA interference (RNAi) assays were conducted on third-instar larvae. Injection of dsHvpnr, dsHvpnr-α and dsHvpnr-β and co-injection of dsHvpnr-α and dsHvpnr-β (dsRNAs mix) all resulted in significant downregulation of the target transcripts. In pupae developed from dsHvpnr-treated larvae, the symmetric black spots on both sides of the mesothorax, metathorax and tergites approached and connected. Pupal morphometric analysis revealed that dsHvpnr, dsHvpnr-α and dsRNAs mix injections significantly narrowed the spacing of dorsal symmetric spots, contracted spiracle distances on tergite sides, diminished pronotum width and markedly reduced inter-compound eye spacing compared to controls. In addition, injections of dsHvpnr and dsRNAs mix significantly reduced the oviposition in female adults. Silencing of Hvpnr led to the disappearance of the scutellum in adults, preventing the elytra from closing and properly attaching to the dorsal side of the abdomen. It is noteworthy that dsHvpnr-α or dsRNA mix induced scutellum formation defects in adults, while knockdown of Hvpnr-β had no impact. Furthermore, in stark contrast to previous studies on ladybird species such as Harmonia axyridis and Coccinella septempunctata, silencing Hvpnr did not affect melanin synthesis in pupae and adults in H. vigintioctopunctata. These findings demonstrate that among the splice variants of Hvpnr, Hvpnr-α plays a dominant regulatory role in the post-embryonic morphogenesis of H. vigintioctopunctata. This study also shows that Hvpnr is not involved in melanin synthesis, indicating significant functional differentiation of pnr during the evolution of ladybirds.

Genomics of insect natural enemies in agroecosystems

Currently, a wealth of genomic data are now accessible for numerous insect natural enemies, serving as valuable resources that deepen our understanding of the genetic basis of biocontrol traits in these organisms. We summarize the current state of genome sequencing and highlight candidate genes related to biocontrol traits that hold promise for genetic improvement. We also review the recent population genomic studies in biological control and the discovery of potential insecticidal genes in parasitoid wasps. Collectively, current genomic works have shown the powerful ability to identify candidate genes responsible for desirable traits or promising effectors. However, further functional study is necessary to gain a mechanistic understanding of these genes, and future efforts are also needed to develop suitable approaches to translate genomic insights into field applications.

Pleiotropic effects of Ebony on pigmentation and development in the Asian multi-coloured ladybird beetle, Harmonia axyridis (Coleoptera: Coccinellidae)

Melanin plays a pivotal role in insect body pigmentation, significantly contributing to their adaptation to diverse biotic and abiotic environmental challenges. Several genes involved in insect melanin synthesis showed pleiotropic effects on insect development and reproduction. Among these, the N-β-alanyl dopamine synthetase gene (Ebony) is integral to the pigmentation process. However, the full spectrum of its pleiotropic impacts is not yet thoroughly understood. In this study, we identified and characterised the HaEbony gene in the Asian multi-coloured ladybird beetle (Harmonia axyridis) and found that HaEbony gene is a conserved gene within the Coleoptera order. We aimed to further explore the multiple roles of HaEbony in the physiology and behaviour in H. axyridis. The CRISPR/Cas9 system was applied to generate multiple HaEbony knockout allele (HaEbony+/-), showing nucleotide deletion in the G0 and G1 generations. Remarkably, the resultant HaEbony+/- mutants consistently displayed darker pigmentation than their wild-type counterparts across larval, pupal and adult stages. Furthermore, these HaEbony+/- individuals (G0) demonstrated an enhanced predatory efficiency, evidenced by a higher number of aphids consumed compared to the wild type. A significant finding was the reduced egg hatchability in both G0 and G1 generations of the HaEbony+/- group, highlighting a potential reproductive fitness cost associated with HaEbony deficiency. In conclusion, our study not only sheds light on the multifaceted roles of HaEbony in H. axyridis but also highlights the potential of employing CRISPR/Cas9-targeted modifications of the Ebony gene. Such genetic interventions could enhance the environmental adaptability and predatory efficacy of ladybirds, presenting a novel strategy in biological control application.

Molecular evolution of dietary shifts in ladybird beetles (Coleoptera: Coccinellidae): from fungivory to carnivory and herbivory

BACKGROUND

Dietary shifts are major evolutionary steps that shape ecological niches and biodiversity. The beetle family Coccinellidae, commonly known as ladybirds, first transitioned from a fungivorous to an insectivorous and subsequently a plant diet. However, the molecular basis of this dietary diversification remained unexplored.

RESULTS

We investigated the molecular evolution of dietary shifts in ladybirds, focusing on the transitions from fungivory to carnivory (Coccinellidae) and from carnivory to herbivory (Epilachnini), by comparing 25 genomes and 62 transcriptomes of beetles. Our analysis shows that chemosensory gene families have undergone significant expansions at both nodes of diet change and were differentially expressed in feeding experiments, suggesting that they may be related to foraging. We found expansions of digestive and detoxifying gene families and losses of chitin-related digestive genes in the herbivorous ladybirds, and absence of most plant cell wall-degrading enzymes in the ladybirds dating from the transition to carnivory, likely indicating the effect of different digestion requirements on the gene repertoire. Immunity effector genes tend to emerge or have specific amino acid sequence compositions in carnivorous ladybirds and are downregulated under suboptimal dietary treatments, suggesting a potential function of these genes related to microbial symbionts in the sternorrhynchan prey.

CONCLUSIONS

Our study provides a comprehensive comparative genomic analysis to address evolution of chemosensory, digestive, detoxifying, and immune genes associated with dietary shifts in ladybirds. Ladybirds can be considered a ubiquitous example of dietary shifts in insects, and thus a promising model system for evolutionary and applied biology.

Molecular and Biochemical Mechanisms of Scutellum Color Variation in Bactrocera dorsalis Adults (Diptera: Tephritidae)

Bactrocera dorsalis (Hendel) is an invasive fruit and vegetable pest, infesting citrus, mango, carambola, etc. We observed that the posterior thoracic scutella of some B. dorsalis adults are yellow, some light yellow, and some white in China. Compared with the B. dorsalis races with a yellow scutellum (YS) and white scutellum (WS), the race with a light-yellow scutellum (LYS) is dominant in citrus and carambola orchards. To reveal genetic correlates among the three races, the genomes of 22 samples (8 with YS, 7 with LYS, and 7 with WS) were sequenced by high-throughput sequencing technology. Single-nucleotide polymorphism (SNP) annotation showed that there were 17,580 non-synonymous mutation sites located in the exonic region. Principal component analysis based on independent SNP data revealed that the SNPs with LYS were more similar to that with YS when compared with WS. Most genes associated with scutellum color variation were involved in three pathways: oxidative phosphorylation, porphyrin and chlorophyll metabolism, and terpenoid backbone biosynthesis. By comparing the sequences among the three races, we screened out 276 differential genes (DGs) in YS vs. WS, 185 DGs in LYS vs. WS, and 104 DGs in YS vs. LYS. Most genes determining color variation in B. dorsalis scutella were located on chromosomes 2-5. Biochemical analysis showed that β-carotene content in YS and LYS was significantly higher than that in WS at any stage of adult days 1, 10, and 20. No significant differences were observed in cytochrome P450 or melanin content in YS, LYS, or WS. Our study provides results on aspects of scutellum color variation in B. dorsalis adults, providing molecular and physiological information for revealing the adaptation and evolution of the B. dorsalis population.

The Genomics Revolution Drives a New Era in Entomology

Thanks to the fast development of sequencing techniques and bioinformatics tools, sequencing the genome of an insect species for specific research purposes has become an increasingly popular practice. Insect genomes not only provide sets of gene sequences but also represent a change in focus from reductionism to systemic biology in the field of entomology. Using insect genomes, researchers are able to identify and study the functions of all members of a gene family, pathway, or gene network associated with a trait of interest. Comparative genomics studies provide new insights into insect evolution, addressing long-lasting controversies in taxonomy. It is also now feasible to uncover the genetic basis of important traits by identifying variants using genome resequencing data of individual insects, followed by genome-wide association analysis. Here, we review the current progress in insect genome sequencing projects and the application of insect genomes in uncovering the phylogenetic relationships between insects and unraveling the mechanisms of important life-history traits. We also summarize the challenges in genome data sharing and possible solutions. Finally, we provide guidance for fully and deeply mining insect genome data.

Egg Microinjection for the Ladybird Beetle Harmonia axyridis

In this paper, we present a detailed protocol for microinjecting DNA, RNA, or protein solutions into fertilized eggs of the multicolored Asian ladybird beetle, Harmonia axyridis, under a stereomicroscope equipped with an injection apparatus. H. axyridis is an emerging model organism for studying various biological fields, showing intraspecific polymorphisms exhibiting highly diverse color patterns on the elytra. Here, we describe how to rear ladybird beetles in a laboratory and obtain fertilized eggs for microinjection experiments. We also provide a constant fluid flow injection method, which enhances the efficiency of microinjection and improves throughput. Our step-by-step protocol is applicable to generating transgenic or genome-edited ladybird beetles, facilitating functional genetics in H. axyridis; the microinjection method should be applicable to other insect eggs. Key features • Egg microinjection for the multicolored Asian ladybird beetle Harmonia axyridis under a stereomicroscope with injection equipment. • Detailed procedures for post-injection care, including rearing Harmonia axyridis. • Step-by-step guide for the efficient collection of Harmonia axyridis eggs.

The genome sequence of the harlequin ladybird, Harmonia axyridis (Pallas, 1773)

We present a genome assembly from an individual female Harmonia axyridis (the harlequin ladybird; Arthropoda; Insecta; Coleoptera; Coccinellidae). The genome sequence is 426 megabases in span. The majority (99.98%) of the assembly is scaffolded into 8 chromosomal pseudomolecules, with the X sex chromosome assembled.

Competing adaptations maintain nonadaptive variation in a wild cricket population

How emerging adaptive variants interact is an important factor in the evolution of wild populations, but the opportunity to empirically study this interaction is rare. We recently documented the emergence of an adaptive phenotype "curly-wing" in Hawaiian populations of field crickets (Teleogryllus oceanicus). Curly-wing inhibits males' ability to sing, protecting them from eavesdropping parasitoid flies (Ormia ochracea). Surprisingly, curly-wing co-occurs with similarly protective silent "flatwing" phenotypes in multiple populations, in which neither phenotype has spread to fixation. These two phenotypes are frequently coexpressed, but since either sufficiently reduces song amplitude to evade the fly, their coexpression confers no additional fitness benefit. Numerous "off-target" phenotypic changes are known to accompany flatwing, and we find that curly-wing, too, negatively impacts male courtship ability and affects mass and survival of females under lab conditions. We show through crosses and genomic and mRNA sequencing that curly-wing expression is associated with variation on a single autosome. In parallel analyses of flatwing, our results reinforce previous findings of X-linked single-locus inheritance. By combining insights into the genetic architecture of these alternative phenotypes with simulations and field observations, we show that the co-occurrence of these two adaptations impedes either from fixing, despite extreme fitness benefits, due to fitness epistasis. This co-occurrence of similar adaptive forms in the same populations might be more common than is generally considered and could be an important force inhibiting adaptive evolution in wild populations of sexually reproducing organisms.

PhyloAln: A Convenient Reference-Based Tool to Align Sequences and High-Throughput Reads for Phylogeny and Evolution in the Omic Era

The current trend in phylogenetic and evolutionary analyses predominantly relies on omic data. However, prior to core analyses, traditional methods typically involve intricate and time-consuming procedures, including assembly from high-throughput reads, decontamination, gene prediction, homology search, orthology assignment, multiple sequence alignment, and matrix trimming. Such processes significantly impede the efficiency of research when dealing with extensive data sets. In this study, we develop PhyloAln, a convenient reference-based tool capable of directly aligning high-throughput reads or complete sequences with existing alignments as a reference for phylogenetic and evolutionary analyses. Through testing with simulated data sets of species spanning the tree of life, PhyloAln demonstrates consistently robust performance compared with other reference-based tools across different data types, sequencing technologies, coverages, and species, with percent completeness and identity at least 50 percentage points higher in the alignments. Additionally, we validate the efficacy of PhyloAln in removing a minimum of 90% foreign and 70% cross-contamination issues, which are prevalent in sequencing data but often overlooked by other tools. Moreover, we showcase the broad applicability of PhyloAln by generating alignments (completeness mostly larger than 80%, identity larger than 90%) and reconstructing robust phylogenies using real data sets of transcriptomes of ladybird beetles, plastid genes of peppers, or ultraconserved elements of turtles. With these advantages, PhyloAln is expected to facilitate phylogenetic and evolutionary analyses in the omic era. The tool is accessible at https://github.com/huangyh45/PhyloAln.

Canola bee pollen is an effective artificial diet additive for improving larval development of predatory coccinellids: a lesson from Harmonia axyridis

BACKGROUND

Pollen is a common plant-derived food source for predatory ladybird beetles under field conditions, yet the potential for pollen to improve the quality of artificial diets remains largely unexplored. In this study, we developed three pollen diets by incorporating varying proportions of canola bee pollen (7.5%, 15.0% and 22.5% with 2.5%, 5.0%, and 7.5% of water, respectively) into a conventional diet. The feeding efficiency of Harmonia axyridis, an omnivorous predator, was evaluated and compared on three pollen diets, a conventional nonpollen diet and pea aphids.

RESULTS

The larvae fed a medium or high pollen diet exhibited significantly higher survival in the 4th instar, pupa and adult stages than those fed a nonpollen diet. These larvae also developed into significantly heavier adults, and their survival rates in adulthood were comparable to those fed pea aphids. Specifically, we revealed the underlying mechanisms through which a high pollen diet enhances pupal development. Consumption of high pollen diet versus nonpollen diet resulted not only in a significant decrease in pupal glycogen content, but also an increase in adult lipid content. Both diet treatments induced similar changes in carbohydrate and glycogen content compared to the aphid diet while exhibiting different alterations in pupal protein content and adult lipid content. Furthermore, the transcriptome analysis revealed that the nutrient metabolism, immune response, and cuticle development pathways were predominantly enriched among the differentially expressed genes (DEGs).

CONCLUSION

Canola bee pollen offers diverse advantages in terms of rearing H. axyridis larvae with an artificial diet, which will advance the development of effective diets for predaceous coccinellids. © 2024 Society of Chemical Industry.

Genetic structure and origin of emu populations in Japanese farms inferred from large-scale SNP genotyping based on double-digest RAD-seq

The emu is a novel poultry species in Japan. However, Japanese farmed emu populations have reduced genetic diversity owing to inbreeding. We have previously suggested that there are genetic resources in the Tohoku Safari Park (TSP) and Fuji/Kakegawa Kachoen Garden Park (FGP/KGP) to extend the genetic diversity of commercial emu farms based on microsatellite (SSR) and mitochondrial DNA. However, those markers provide relatively poor information. Thus, we investigated the genetic structure of farmed Japanese populations based on a large-scale genotyping system using RAD-seq and verified the usefulness of TSP and FGP/KGP as genetic resources for expanding genetic diversity. Admixture, phylogenetic, and principal component analyses based on 28,676 SNPs showed that TSP individuals were ancestors in the Okhotsk Emu Farm (OEF). FGP/KGP individuals showed a unique genetic component that differed from that of the others. We have previously reported that the mitochondrial haplotypes of FGP/KGP were shared with an isolated wild population in eastern Australia. These results suggest that FGP/KGP individuals originated from an eastern Australia isolated population different from other populations including ancestral of OEF/TSP. Our results would provide information for the development of Japanese emu farms and industry and for the conservation of genetic resources in the Australian wild emu.

Transposable Element Insertions Are Associated with Batesian Mimicry in the Pantropical Butterfly Hypolimnas misippus

Hypolimnas misippus is a Batesian mimic of the toxic African Queen butterfly (Danaus chrysippus). Female H. misippus butterflies use two major wing patterning loci (M and A) to imitate three color morphs of D. chrysippus found in different regions of Africa. In this study, we examine the evolution of the M locus and identify it as an example of adaptive atavism. This phenomenon involves a morphological reversion to an ancestral character that results in an adaptive phenotype. We show that H. misippus has re-evolved an ancestral wing pattern present in other Hypolimnas species, repurposing it for Batesian mimicry of a D. chrysippus morph. Using haplotagging, a linked-read sequencing technology, and our new analytical tool, Wrath, we discover two large transposable element insertions located at the M locus and establish that these insertions are present in the dominant allele responsible for producing mimetic phenotype. By conducting a comparative analysis involving additional Hypolimnas species, we demonstrate that the dominant allele is derived. This suggests that, in the derived allele, the transposable elements disrupt a cis-regulatory element, leading to the reversion to an ancestral phenotype that is then utilized for Batesian mimicry of a distinct model, a different morph of D. chrysippus. Our findings present a compelling instance of convergent evolution and adaptive atavism, in which the same pattern element has independently evolved multiple times in Hypolimnas butterflies, repeatedly playing a role in Batesian mimicry of diverse model species.

Chromosome-level genome assembly of the giant ladybug Megalocaria dilatata

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

Structure and mechanical properties of ladybird elytra as biological sandwich panels

The armour of the ladybird, elytra, protect the body from injury and are well-adapted to flight. However, experimental methods to decipher their mechanical performances had been challenging due to the small size, making it unclear how the elytra balance mass and strength. Here, we provide insights to the relationship between the microstructure and multifunctional properties of the elytra by means of structural characterization, mechanical analysis and finite element simulations. Micromorphology analysis on the elytron revealed the thickness ratio of the upper lamination, middle layer and lower lamination is approximately 51:139:7. The upper lamination had multiple cross fibre layers and the thickness of each fibre layer is not the same. In addition, the tensile strength, elastic modulus, fracture strain, bending stiffness and hardness of elytra were obtained through in-situ tensile and nanoindentation-bending under the influence of multiple loading conditions, which also serve as references for finite element models. The finite element model revealed that structural factors such as thickness of each layer, angle of fibre layer and trabeculae are key to affecting the mechanical properties, but the effect is different. When the thickness of upper, middle and lower layers is the same, the tensile strength provided by unit mass of the model is 52.78% lower than that provided by elytra. These findings broaden the relationship between the structural and mechanical properties of the ladybird elytra, and are expected to inspire the development of sandwich structures in biomedical engineering.

The genome sequence of the orange ladybird, Halyzia sedecimguttata (Linnaeus, 1758)

We present a genome assembly from an individual Halyzia sedecimguttata (the orange ladybird, Arthropoda; Insecta; Coleoptera, Coccinellidae). The genome sequence is 919.1 megabases in span. Most of the assembly is scaffolded into 10 chromosomal pseudomolecules, including the X sex chromosome. The mitochondrial genome has also been assembled and is 21.0 kilobases in length. Gene annotation of this assembly on Ensembl identified 27,547 protein coding genes.

A chromosome-level genome assembly of the Henosepilachna vigintioctomaculata provides insights into the evolution of ladybird beetles

The ladybird beetle Henosepilachna vigintioctomaculata is an economically significant oligophagous pest that induces damage to many Solanaceae crops. An increasing number of studies have examined the population and phenotype diversity of ladybird beetles. However, few comparative genome analyses of ladybird beetle species have been conducted. Here, we obtained a high-quality chromosome-level genome assembly of H. vigintioctomaculata using various sequencing technologies, and the chromosome-level genome assembly was ~581.63 Mb, with 11 chromosomes successfully assembled. The phylogenetic analysis showed that H. vigintioctomaculata is a more ancient lineage than the other three sequenced ladybird beetles, Harmonia axyridis, Propylea japonica, and Coccinella septempunctata. We also compared positively selected genes (PSGs), transposable elements (TEs) ratios and insertion times, and key gene families associated with environmental adaptation among these ladybird beetles. The pattern of TEs evolution of H. vigintioctomaculata differs from the other three ladybird beetles. The PSGs were associated with ladybird beetles development. However, the key gene families associated with environmental adaptation in ladybird beetles varied. Overall, the high-quality draft genome sequence of H. vigintioctomaculata provides a useful resource for studies of beetle biology, especially for the invasive biology of ladybird beetles.

Melanin Synthesis Pathway Interruption: CRISPR/Cas9-mediated Knockout of dopa decarboxylase (DDC) in Harmonia axyridis (Coleoptera: Coccinellidae)

CRISPR/Cas9 technology is a very powerful genome editing tool and has been used in many insect species for functional genomics studies through targeted gene mutagenesis. Here, we successfully established CRISPR/Cas9 research platform in Asian multi-colored ladybird beetle, Harmonia axyridis, an important natural enemy in biological control. In this study, one pivotal gene dopa decarboxylase (DDC) in melanin synthesis was targeted by CRISPR/Cas9 to generate mutants in H. axyridis by CRISPR/Cas9 technology. Our results showed that injection of single guide RNA of the DDC and Cas9 protein into preblastoderm eggs induced one insertion and four deletion (indels) mutant H. axyridis. Mutations of HaDDC gene generated 25% mutant rate with melanin missing phenotype in larva, pupa,l and adult stage. The predation ability of the fourth instar larvae has no significant difference between wild (control) and mutant H. axyridis (G0), while these mutant fourth instar larvae had longer developmental period than that of the wild type. Consequently, the total predation of the fourth instar larvae was significantly increased in H. axyridis mutants comparing with the wild type. These results indicated that the success of CRISPR/Cas9 gene editing in H. axyridis. The gene editing platform in H. axyridis would facilitate the gene function research and promote special strain of predatory ladybird beetle generation.

The evolution and diversification of oakleaf butterflies

Oakleaf butterflies in the genus Kallima have a polymorphic wing phenotype, enabling these insects to masquerade as dead leaves. This iconic example of protective resemblance provides an interesting evolutionary paradigm that can be employed to study biodiversity. We integrated multi-omic data analyses and functional validation to infer the evolutionary history of Kallima species and investigate the genetic basis of their variable leaf wing patterns. We find that Kallima butterflies diversified in the eastern Himalayas and dispersed to East and Southeast Asia. Moreover, we find that leaf wing polymorphism is controlled by the wing patterning gene cortex, which has been maintained in Kallima by long-term balancing selection. Our results provide macroevolutionary and microevolutionary insights into a model species originating from a mountain ecosystem.

[Small scale evolution]

Small-scale evolution or microevolution concerns evolution at the intra-specific level or between closely related species. At the intra-specific level, it allows the analysis of the evolutionary forces at work: mutation, genetic drift, migration and selection. Moreover, because of the short evolutionary time, it is easier to identify the genetic basis of observed phenotypic differences. Most studies focus on current populations but more and more analyses are performed on ancient DNA. This provides important information for tracing the history of populations and also allows the reconstruction of phenotypes of individuals that disappeared several thousand years ago. In this short review, I present studies showing how pre-zygotic or post-zygotic barriers involved in species formation are set up using the example of the geographical barrier due to the formation of the Isthmus of Panama and that of the heterochromatin divergence in Drosophilidae. I also describe the different approaches that have been used to identify the genetic basis of well known phenotypic variations: candidate gene approach (about melanism in felines), QTL mapping (variation in the number of lateral bone plates in sticklebacks), association study (pigmentation in the Asian ladybird). Finally, I illustrate the key impact of natural selection with the iconic example of the evolution of the beak of Galapagos finches, and the role of certain developmental genes in its morphological diversification.

Exploring new genomic territories with emerging model insects

Improvements in reference genome generation for insects and across the tree of life are extending the concept and utility of model organisms beyond traditional laboratory-tractable supermodels. Species or groups of species with comprehensive genome resources can be developed into model systems for studying a large variety of biological phenomena. Advances in sequencing and assembly technologies are supporting these emerging genome-enabled model systems by producing resources that are increasingly accurate and complete. Nevertheless, quality controls including assessing gene content completeness are required to ensure that these data can be included in expanding catalogues of high-quality references that will greatly advance understanding of insect biology and evolution.

The Molecular Properties and Roles of Pannier in Harmonia axyridis's Metamorphosis and Melanin Synthesis

The GATA transcription factor Pannier is identified as the major regulatory gene in color pattern formation in the Asian multi-colored ladybird beetle (Harmonia axyridis). however, the mechanisms of Pannier in regulating melanin synthesis and development in H. axyridis remain elusive. In this study, we identified and characterized Pannier in H. axyridis (HaPnr) and showed it to have two alternative spliced variants named HaPnr-α and HaPnr-β. Analyses of developmental stage expression patterns revealed that HaPnr, HaPnr-α and HaPnr-β were constitutively expressed throughout all developmental stages. To examine the role of HaPnr in H. axyridis development, RNA interference was performed in late larvae (the fourth instar) and early pupae (the first day of pupa stage). The transcript levels of HaPnr were effectively suppressed after the injection of double-stranded RNA of HaPnr (dsHaPnr). The fourth instar larvae injected with dsHaPnr reduced the pupation rates to only 61.50%, compared with 88.5% in the dsGFP-injected group. The un-pupated larvae gradually died after 1 week, and visually unaffected pupae emerged into abnormal adults with malformed hind wings and melanin absent from the cuticle. These abnormal adults gradually died 10 days after eclosion. However, when early pupae were injected with dsHaPnr, the normal eclosion rate was achieved at 88.41% on day 6 after the injection. In addition, these successful eclosion adults also showed an absence of melanin in the cuticle, but they could mate normally and have normal fecundity as compared with the control. We further demonstrated that the suppression of HaPnr-α or HaPnr-β individually did not affect the pupation and eclosion process. The suppression of HaPnr-α expression resulted in elytra melanin decreasing in both the conspicua and the succinea subgroup in H. axyridis. Even though the suppression of HaPnr-β expression only affected the melanin synthesis in the succinea subgroup, it significantly prolonged the time taken for melanin synthesis to occur in the conspicua subgroup in H. axyridis. These results indicate that HaPnr plays an essential role in insect development, especially during their metamorphosis, and also support our hypothesis that HaPnr could regulate melanin synthesis in H. axyridis under the combined action with its two splicing variants, HaPnr-α and HaPnr-β.

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.

Recombinant inbred lines and next-generation sequencing enable rapid identification of candidate genes involved in morphological and agronomic traits in foxtail millet

We constructed recombinant inbred lines (RILs) between a Japanese and a Taiwanese landrace of foxtail millet and employed next-generation sequencing, such as flexible ddRAD-seq and Nanopore sequencing to identify the candidate genes involved in the crop evolution of foxtail millet. We successfully constructed a linkage map using flexible ddRAD-seq with parents and RILs and detected major QTLs for each of three traits: leaf sheath colors, spikelet-tipped bristles (stb), and days to heading (DTH). (1) For leaf sheath colors, we identified the C gene on chromosome IV. (2) We identified a homeobox (HOX14) gene for stb on chromosome II, which shows homology with HvVrs1 in barley. (3) Finally, we identified a QTL with a large effect on DTH on chromosome II. A parent of the RILs from Taiwan and Yugu1 had a Harbinger-like TE in intron 3 of this gene. We also investigated the geographical distribution of the TE insertion type of this gene and found that the insertion type is distributed in the northern part of East Asia and intensively in South and Southeast Asia, suggesting that loss/reduction of function of this gene plays an important role in spreading into the northern part of East Asia and subtropical and tropical zones.

The potential pigmentation-related genes in spider mites revealed by comparative transcriptomes of the red form of Tetranychus urticae

Colouration in spider mites is due to the presence of carotenoids with diverse colours, including yellows, oranges and reds. Tetranychus urticae has two main colour forms, red and green. Although a ketolase has been implicated in determining the colour, the underlying genetic basis of body colour divergence between the two forms has remained unclear. Based on a combination of comparative transcriptomes and RNA interference, we found that a gene encoding a cytochrome P450 enzyme of the CYP4 clan (CYP389B1) had remarkably high expression in adult females of the red T. urticae, as well as in hybrids obtained by crossing the red and green forms. Down-regulation of this gene by RNA interference resulted in decreased accumulation of red pigment. Up-regulation of the expressions of a scavenger receptor gene (SCARB1) and a mitochondrial glycine transporter (SLC25A38) also strongly contributed to red colour development in adult females. Suppressing the mRNA levels of these genes also resulted in reduced accumulation of red pigment in the three other spider mites with red body colour. Our results provide evidence that the body colour divergence between the two forms is caused by different expressions of pigmentation-related genes, and point to a possible role of a novel cytochrome P450 gene (CYP389B1) in regulating red-orange body colour. These findings expand the number of candidate cytochrome P450 genes involved in endogenous pigmentation and will help to understand their roles in determining colour patterns in mites and other species.

Aspartate-β-alanine-NBAD pathway regulates pupal melanin pigmentation plasticity of ladybird Harmonia axyridis

Phenotypic plasticity is observed in many animal species and it is effective for them to cope with many types of environmental threats. The multicolored Asian ladybird Harmonia axyridis shows a cuticular pigmentation plasticity that can be rapidly induced by temperature changes, and in the form of changeable melanin spot patterns to adjust heat-absorbing. Here, H. axyridis with thermal stimulation were selected for determining the molecular regulations behind it. First, we confirmed the melanin level changes of H. axyridis pupa could be induced by temperature, and then screened the efficient time window for thermal sensing of H. axyridis pre-pupa; it is suggested that the late stage of pre-pupa (late stage of 4th instar larva) is the critical period to sense thermal signals and adjust its pupal melanin spot area size to adapt to upcoming thermal conditions. The Ha-ADC (aspartate decarboxylase) and Ha-ebony (NBAD synthase) of aspartate-β-alanine-NBAD pathway were then proved in regulation of cuticular melanization for pupa through RNA interference experiments; knockdown of these two genes enlarged the melanin spot size. Finally, we designed a random injection of Ha-ADC at different pre-pupal stages, to further study the regulation window during this process. Combined with all evidence observed, we suggested the spot size determination can be regulated very close to the time point of pupation, and genes of the aspartate-β-alanine-NBAD pathway play an important role at the molecular level. In brief, H. axyridis exhibits a flexible active physiological regulation through transcriptional modification to thermal changes.

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

The draft genome of the specialist flea beetle Altica viridicyanea (Coleoptera: Chrysomelidae)

BACKGROUND

Altica (Coleoptera: Chrysomelidae) is a highly diverse and taxonomically challenging flea beetle genus that has been used to address questions related to host plant specialization, reproductive isolation, and ecological speciation. To further evolutionary studies in this interesting group, here we present a draft genome of a representative specialist, Altica viridicyanea, the first Alticinae genome reported thus far.

RESULTS

The genome is 864.8 Mb and consists of 4490 scaffolds with a N50 size of 557 kb, which covered 98.6% complete and 0.4% partial insect Benchmarking Universal Single-Copy Orthologs. Repetitive sequences accounted for 62.9% of the assembly, and a total of 17,730 protein-coding gene models and 2462 non-coding RNA models were predicted. To provide insight into host plant specialization of this monophagous species, we examined the key gene families involved in chemosensation, detoxification of plant secondary chemistry, and plant cell wall-degradation.

CONCLUSIONS

The genome assembled in this work provides an important resource for further studies on host plant adaptation and functionally affiliated genes. Moreover, this work also opens the way for comparative genomics studies among closely related Altica species, which may provide insight into the molecular evolutionary processes that occur during ecological speciation.

The color pattern inducing gene wingless is expressed in specific cell types of campaniform sensilla of a polka-dotted fruit fly, Drosophila guttifera

A polka-dotted fruit fly, Drosophila guttifera, has a unique pigmentation pattern on its wings and is used as a model for evo-devo studies exploring the mechanism of evolutionary gain of novel traits. In this species, a morphogen-encoding gene, wingless, is expressed in species-specific positions and induces a unique pigmentation pattern. To produce some of the pigmentation spots on wing veins, wingless is thought to be expressed in developing campaniform sensillum cells, but it was unknown which of the four cell types there express(es) wingless. Here we show that two of the cell types, dome cells and socket cells, express wingless, as indicated by in situ hybridization together with immunohistochemistry. This is a unique case in which non-neuronal SOP (sensory organ precursor) progeny cells produce Wingless as an inducer of pigmentation pattern formation. Our finding opens a path to clarifying the mechanism of evolutionary gain of a unique wingless expression pattern by analyzing gene regulation in dome cells and socket cells.

Evo-devo of wing colour patterns in beetles

Insects have evolved tremendously diverse wing colour patterns that fulfil ecologically vital functions, including intraspecific sexual signalling, mimesis, mimicry, and detering predators. Beetles, which form the most species-rich order Coleoptera, have amazingly diverse wing colour patterns; however, the molecular mechanisms that give rise to these patterns remain poorly understood. Recently, the gene pannier (pnr), which encodes a transcription factor of the GATA family, was identified as an essential player in the wing patterning of the multi-coloured Asian ladybird beetle Harmonia axyridis. Here we review recent progress in understanding the molecular underpinnings of wing colour pattern formation in H. axyridis.

A chromosome-level assembly of the harlequin ladybird Harmonia axyridis as a genomic resource to study beetle and invasion biology

The harlequin ladybird, Harmonia axyridis (Pallas), is a well-known model organism for genetic studies and is also a well-studied natural enemy used for pest control. It became an invasive species after being introduced to North America and Europe as a pest control agent. Though two genome assemblies for this insect have been previously reported, a high-quality genome assembly at the chromosome level is still not available. Here, we obtained a new chromosome-level genome assembly of H. axyridis by combining various sequencing technologies, namely Illumina short reads, PacBio long reads, 10X Genomics and Hi-C. The chromosome-level genome assembly is 423 Mb with a scaffold N50 of 45.92 Mb. Using Hi-C data 1,897 scaffolds were anchored to eight chromosomes. A total of 730,068 repeat sequences were identified, making up 51.2% of the assembled genome. After masking these repeat sequences, we annotated 22,810 protein-encoding genes. The X chromosome and Y-linked scaffolds were also identified by resequencing male and female genomes and calculating the male to female coverage ratios. Two gene families associated with environmental adaptation, odorant receptor and cytochrome P450, were analysed and showed no obvious expansion in H. axyridis. We successfully constructed a putative biosynthesis pathway of harmonine, a defence compound in the haemolymph of H. axyridis, which is a key factor for H. axyridis strong immunity. The chromosome-level genome assembly of H. axyridis is a helpful resource for studies of beetle biology and invasive biology.

Horizontally acquired antibacterial genes associated with adaptive radiation of ladybird beetles

BACKGROUND

Horizontal gene transfer (HGT) has been documented in many herbivorous insects, conferring the ability to digest plant material and promoting their remarkable ecological diversification. Previous reports suggest HGT of antibacterial enzymes may have contributed to the insect immune response and limit bacterial growth. Carnivorous insects also display many evolutionary successful lineages, but in contrast to the plant feeders, the potential role of HGTs has been less well-studied.

RESULTS

Using genomic and transcriptomic data from 38 species of ladybird beetles, we identified a set of bacterial cell wall hydrolase (cwh) genes acquired by this group of beetles. Infection with Bacillus subtilis led to upregulated expression of these ladybird cwh genes, and their recombinantly produced proteins limited bacterial proliferation. Moreover, RNAi-mediated cwh knockdown led to downregulation of other antibacterial genes, indicating a role in antibacterial immune defense. cwh genes are rare in eukaryotes, but have been maintained in all tested Coccinellinae species, suggesting that this putative immune-related HGT event played a role in the evolution of this speciose subfamily of predominant predatory ladybirds.

CONCLUSION

Our work demonstrates that, in a manner analogous to HGT-facilitated plant feeding, enhanced immunity through HGT might have played a key role in the prey adaptation and niche expansion that promoted the diversification of carnivorous beetle lineages. We believe that this represents the first example of immune-related HGT in carnivorous insects with an association with a subsequent successful species radiation.

Satellitome Analysis in the Ladybird Beetle Hippodamia variegata (Coleoptera, Coccinellidae)

Hippodamia variegata is one of the most commercialized ladybirds used for the biological control of aphid pest species in many economically important crops. This species is the first Coccinellidae whose satellitome has been studied by applying new sequencing technologies and bioinformatics tools. We found that 47% of the H. variegata genome is composed of repeated sequences. We identified 30 satellite DNA (satDNA) families with a median intragenomic divergence of 5.75% and A+T content between 45.6% and 74.7%. This species shows satDNA families with highly variable sizes although the most common size is 100–200 bp. However, we highlight the existence of a satDNA family with a repeat unit of 2 kb, the largest repeat unit described in Coleoptera. PCR amplifications for fluorescence in situ hybridization (FISH) probe generation were performed for the four most abundant satDNA families. FISH with the most abundant satDNA family as a probe shows its pericentromeric location on all chromosomes. This location is coincident with the heterochromatin revealed by C-banding and DAPI staining, also analyzed in this work. Hybridization signals for other satDNA families were located only on certain bivalents and the X chromosome. These satDNAs could be very useful as chromosomal markers due to their reduced location.

A linked-read approach to museomics: Higher quality de novo genome assemblies from degraded tissues

High-throughput sequencing technologies are a proposed solution for accessing the molecular data in historical specimens. However, degraded DNA combined with the computational demands of short-read assemblies has posed significant laboratory and bioinformatics challenges for de novo genome assembly. Linked-read or "synthetic long-read" sequencing technologies, such as 10× Genomics, may provide a cost-effective alternative solution to assemble higher quality de novo genomes from degraded tissue samples. Here, we compare assembly quality (e.g., genome contiguity and completeness, presence of orthogroups) between four new deer mouse (Peromyscus spp.) genomes assembled using linked-read technology and four published genomes assembled from a single shotgun library. At a similar price-point, these approaches produce vastly different assemblies, with linked-read assemblies having overall higher contiguity and completeness, measured by larger N50 values and greater number of genes assembled, respectively. As a proof-of-concept, we used annotated genes from the four Peromyscus linked-read assemblies and eight additional rodent taxa to generate a phylogeny, which reconstructed the expected relationships among species with 100% support. Although not without caveats, our results suggest that linked-read sequencing approaches are a viable option to build de novo genomes from degraded tissues, which may prove particularly valuable for taxa that are extinct, rare or difficult to collect.

Phylotranscriptomics of the Pentapetalae Reveals Frequent Regulatory Variation in Plant Local Responses to the Fungal Pathogen Sclerotinia sclerotiorum

Quantitative disease resistance (QDR) is a conserved form of plant immunity that limits infections caused by a broad range of pathogens. QDR has a complex genetic determinism. The extent to which molecular components of the QDR response vary across plant species remains elusive. The fungal pathogen Sclerotinia sclerotiorum, causal agent of white mold diseases on hundreds of plant species, triggers QDR in host populations. To document the diversity of local responses to S. sclerotiorum at the molecular level, we analyzed the complete transcriptomes of six species spanning the Pentapetalae (Phaseolus vulgaris, Ricinus communis, Arabidopsis [Arabidopsis thaliana], Helianthus annuus, Solanum lycopersicum, and Beta vulgaris) inoculated with the same strain of S. sclerotiorum About one-third of plant transcriptomes responded locally to S. sclerotiorum, including a high proportion of broadly conserved genes showing frequent regulatory divergence at the interspecific level. Evolutionary inferences suggested a trend toward the acquisition of gene induction relatively recently in several lineages. Focusing on a group of ABCG transporters, we propose that exaptation by regulatory divergence contributed to the evolution of QDR. This evolutionary scenario has implications for understanding the QDR spectrum and durability. Our work provides resources for functional studies of gene regulation and QDR molecular mechanisms across the Pentapetalae.

Mimicry diversification in Papilio dardanus via a genomic inversion in the regulatory region of engrailed-invected

Polymorphic Batesian mimics exhibit multiple protective morphs that each mimic a different noxious model. Here, we study the genomic transitions leading to the evolution of different mimetic wing patterns in the polymorphic Mocker Swallowtail Papilio dardanus. We generated a draft genome (231 Mb over 30 chromosomes) and re-sequenced individuals of three morphs. Genome-wide single nucleotide polymorphism (SNP) analysis revealed elevated linkage disequilibrium and divergence between morphs in the regulatory region of engrailed, a developmental gene previously implicated in the mimicry switch. The diverged region exhibits a discrete chromosomal inversion (of 40 kb) relative to the ancestral orientation that is associated with the cenea morph, but not with the bottom-recessive hippocoonides morph or with non-mimetic allopatric populations. The functional role of this inversion in the expression of the novel phenotype is currently unknown, but by preventing recombination, it allows the stable inheritance of divergent alleles enabling geographic spread and local coexistence of multiple adaptive morphs.

Superficially Similar Adaptation Within One Species Exhibits Similar Morphological Specialization but Different Physiological Regulations and Origins

Animals have developed numerous strategies to contend with environmental pressures. We observed that the same adaptation strategy may be used repeatedly by one species in response to a certain environmental challenge. The ladybird Harmonia axyridis displays thermal phenotypic plasticity at different developmental stages. It is unknown whether these superficially similar temperature-induced specializations share similar physiological mechanisms. We performed various experiments to clarify the differences and similarities between these processes. We examined changes in the numbers and sizes of melanic spots in pupae and adults, and confirmed similar patterns for both. The dopamine pathway controls pigmentation levels at both developmental stages of H. axyridis. However, the aspartate-β-alanine pathway controls spot size and number only in the pupae. An upstream regulation analysis revealed the roles of Hox genes and elytral veins in pupal and adult spot formation. Both the pupae and the adults exhibited similar morphological responses to temperatures. However, they occurred in different body parts and were regulated by different pathways. These phenotypic adaptations are indicative of an effective thermoregulatory system in H. axyridis and explains how insects contend with certain environmental pressure based on various control mechanisms.

The genomics of coloration provides insights into adaptive evolution

Coloration is an easily quantifiable visual trait that has proven to be a highly tractable system for genetic analysis and for studying adaptive evolution. The application of genomic approaches to evolutionary studies of coloration is providing new insight into the genetic architectures underlying colour traits, including the importance of large-effect mutations and supergenes, the role of development in shaping genetic variation and the origins of adaptive variation, which often involves adaptive introgression. Improved knowledge of the genetic basis of traits can facilitate field studies of natural selection and sexual selection, making it possible for strong selection and its influence on the genome to be demonstrated in wild populations.

Evolution of wing pigmentation in Drosophila: Diversity, physiological regulation, and cis-regulatory evolution

Fruit flies (Drosophila and its close relatives, or “drosophilids”) are a group that includes an important model organism, Drosophila melanogaster, and also very diverse species distributed worldwide. Many of these species have black or brown pigmentation patterns on their wings, and have been used as material for evo‐devo research. Pigmentation patterns are thought to have evolved rapidly compared with body plans or body shapes; hence they are advantageous model systems for studying evolutionary gains of traits and parallel evolution. Various groups of drosophilids, including genus Idiomyia (Hawaiian Drosophila), have a variety of pigmentations, ranging from simple black pigmentations around crossveins to a single antero‐distal spot and a more complex mottled pattern. Pigmentation patterns are sometimes obviously used for sexual displays; however, in some cases they may have other functions. The process of wing formation in Drosophila, the general mechanism of pigmentation formation, and the transport of substances necessary for pigmentation, including melanin precursors, through wing veins are summarized here. Lastly, the evolution of the expression of genes regulating pigmentation patterns, the role of cis‐regulatory regions, and the conditions required for the evolutionary emergence of pigmentation patterns are discussed. Future prospects for research on the evolution of wing pigmentation pattern formation in drosophilids are presented, particularly from the point of view of how they compare with other studies of the evolution of new traits.

Molecular and Potential Regulatory Mechanisms of Melanin Synthesis in Harmonia axyridis

Melanization is a common phenomenon in insects, and melanin synthesis is a conserved physiological process that occurs in epidermal cells. Moreover, a comprehensive understanding of the mechanisms of melanin synthesis influencing insect pigmentation are well-suited for investigating phenotype variation. The Asian multi-colored (Harlequin) ladybird beetle, Harmonia axyridis, exhibits intraspecific polymorphism based on relative levels of melanization. However, the specific characteristics of melanin synthesis in H. axyridis remains elusive. In this study, we performed gene-silencing analysis of the pivotal inverting enzyme, tyrosine hydroxylase (TH), and DOPA decarboxylase (DDC) in the tyrosine metabolism pathway to investigate the molecular and regulatory mechanism of melanin synthesis in H. axyridis. Using RNAi of TH and DDC genes in fourth instar larvae, we demonstrated that dopamine melanin was the primary contributor to the overall body melanization of H. axyridis. Furthermore, our study provides the first conclusive evidence that dopamine serves as a melanin precursor for synthesis in the early pupal stage. According to transcription factor Pannier, which is essential for the formation of melanic color on the elytra in H. axyridis, we further demonstrated that suppression of HaPnr can significantly decrease expression levels of HaTH and HaDDC. These results in their entirety lead to the conclusion that transcription factor Pannier can regulate dopamine melanin synthesis in the dorsal elytral epidermis of H. axyridis.

The Role of the Dopamine Melanin Pathway in the Ontogeny of Elytral Melanization in Harmonia axyridis

Polymorphic melanism in insects is a conspicuous phenotype which is derived from specific genotypes, and might be central to speciation processes via assortative sexual selection. At the molecular level, melanism in insects is attributed to the melanin pathway. DOPA decarboxylase (DDC) protein encoded by the DDC gene plays a central role in dopamine-melanin synthesis, the main component of melanin in insects. Although the mechanism of melanism has been elucidated in holometabolous insects, other physiological processes coupled with melanin synthesis are unknown. Herein, we identified DDC from the Asian multi-colored ladybird (Harmonia axyridis), an ideal holometabolous insect for studies of melanization due to highly variable color on their elytra. Analyses revealed that HaDDC (the DDC gene of H. axyridis) was constitutively expressed throughout all developmental stages. We performed RNAi technique to examine the melanin synthesis pathway of elytra in H. axyridis. The transcript levels of HaDDC were significantly suppressed after the injection of double-strand RNA of HaDDC (dsHaDDC) at 300 ng/individual in third instar larvae. Silencing HaDDC in third instar larvae did not result in mortality nor significantly affect pupation and eclosion. We further demonstrated that all adults of H. axyridis (forms succinea, spectabilis, and conspicua) with HaDDC silenced in third larvae showed abnormal phenotype which emerged as decreased elytra melanin. However, melanin was still observed in other parts of the adults such as head or pronotum. These results demonstrate for the first time that dopamine-derived melanin is the main contributor in elytra melanization in H. axyridis. Additionally, we provide evidence for DDC in regulating fecundity by showing that silencing of HaDDC in third instar larvae significantly reduced female egg-laying and egg hatching. As such, DDC is likely pleiotropic in respect of its role in melanin production and fecundity processes. These findings bring novel insights into melanin production in holometabolous insects, and contribute to the framework on which further studies may be conducted on the mechanism of pigment production and patterning in various types of insect coloration.

Coleoptera genome and transcriptome sequences reveal numerous differences in neuropeptide signaling between species

Background

Insect neuropeptides are interesting for the potential their receptors hold as plausible targets for a novel generation of pesticides. Neuropeptide genes have been identified in a number of different species belonging to a variety of insects. Results suggest significant neuropeptide variation between different orders, but much less is known of neuropeptidome variability within an insect order. I therefore compared the neuropeptidomes of a number of Coleoptera.

Methodology

Publicly available genome sequences, transcriptomes and the original sequence data in the form of short sequence read archives were analyzed for the presence or absence of genes coding neuropeptides as well as some neuropeptide receptors in seventeen beetle species.

Results

Significant differences exist between the Coleoptera analyzed here, while many neuropeptides that were previously characterized from Tribolium castaneum appear very similar in all species, some are not and others are lacking in one or more species. On the other hand, leucokinin, which was presumed to be universally absent from Coleoptera, is still present in non-Polyphaga beetles.

Conclusion

The variability in neuropeptidome composition between species from the same insect order may be as large as the one that exists between species from different orders.

Development and evolution of color patterns in ladybird beetles: A case study in Harmonia axyridis

Many organisms show various geometric color patterns on their bodies, and the developmental, evolutionary, genetic, and ecological bases of these patterns have been intensely studied in various organisms. Ladybird beetles display highly diverse patterns of wing (elytral) color and are one of the most attractive model organisms for studying these characteristics. In this study, we reviewed the genetic history of elytral color patterns in the Asian multicolored ladybird beetle Harmonia axyridis from the classical genetic studies led by the pupils of Thomas Hunt Morgan and Theodosius Dobzhansky to recent genomic studies that revealed that a single GATA transcription factor gene, pannier, regulates the highly diverse elytral color patterns in this species. We also reviewed and discussed the developmental and evolutionary mechanisms driven by the pannier locus in H. axyridis. In the development sections, we focused on the following two topics: (a) how the red (carotenoid) and black (melanin) pigmentation of elytra is regulated by the pannier and pigmentation genes and (b) how the diverse color patterns are formed by integrating regulatory inputs from other genes involved in wing development. In the evolution section, we subsequently focused on the highly diversified DNA sequences within the first intron of pannier that are 56-76 kb long and that were generated through recurrent multiple inversions. Furthermore, we discussed how these recurrent inversions have driven the diversification of color patterns throughout evolution.

The Genomic Basis of Color Pattern Polymorphism in the Harlequin Ladybird

Many animal species comprise discrete phenotypic forms. A common example in natural populations of insects is the occurrence of different color patterns, which has motivated a rich body of ecological and genetic research [1-6]. The occurrence of dark, i.e., melanic, forms displaying discrete color patterns is found across multiple taxa, but the underlying genomic basis remains poorly characterized. In numerous ladybird species (Coccinellidae), the spatial arrangement of black and red patches on adult elytra varies wildly within species, forming strikingly different complex color patterns [7, 8]. In the harlequin ladybird, Harmonia axyridis, more than 200 distinct color forms have been described, which classic genetic studies suggest result from allelic variation at a single, unknown, locus [9, 10]. Here, we combined whole-genome sequencing, population-based genome-wide association studies, gene expression, and functional analyses to establish that the transcription factor Pannier controls melanic pattern polymorphism in H. axyridis. We show that pannier is necessary for the formation of melanic elements on the elytra. Allelic variation in pannier leads to protein expression in distinct domains on the elytra and thus determines the distinct color patterns in H. axyridis. Recombination between pannier alleles may be reduced by a highly divergent sequence of ∼170 kb in the cis-regulatory regions of pannier, with a 50 kb inversion between color forms. This most likely helps maintain the distinct alleles found in natural populations. Thus, we propose that highly variable discrete color forms can arise in natural populations through cis-regulatory allelic variation of a single gene.