High-quality haploid genomes corroborate 29 chromosomes and highly conserved synteny of genes in Hyles hawkmoths (Lepidoptera: Sphingidae)

BACKGROUND

Morphological and traditional genetic studies of the young Pliocene genus Hyles have led to the understanding that despite its importance for taxonomy, phenotypic similarity of wing patterns does not correlate with phylogenetic relationship. To gain insights into various aspects of speciation in the Spurge Hawkmoth (Hyles euphorbiae), we assembled a chromosome-level genome and investigated some of its characteristics.

RESULTS

The genome of a male H. euphorbiae was sequenced using PacBio and Hi-C data, yielding a 504 Mb assembly (scaffold N50 of 18.2 Mb) with 99.9% of data represented by the 29 largest scaffolds forming the haploid chromosome set. Consistent with this, FISH analysis of the karyotype revealed n = 29 chromosomes and a WZ/ZZ (female/male) sex chromosome system. Estimates of chromosome length based on the karyotype image provided an additional quality metric of assembled chromosome size. Rescaffolding the published male H. vespertilio genome resulted in a high-quality assembly (651 Mb, scaffold N50 of 22 Mb) with 98% of sequence data in the 29 chromosomes. The larger genome size of H. vespertilio (average 1C DNA value of 562 Mb) was accompanied by a proportional increase in repeats from 45% in H. euphorbiae (measured as 472 Mb) to almost 55% in H. vespertilio. Several wing pattern genes were found on the same chromosomes in the two species, with varying amounts and positions of repetitive elements and inversions possibly corrupting their function.

CONCLUSIONS

Our two-fold comparative genomics approach revealed high gene synteny of the Hyles genomes to other Sphingidae and high correspondence to intact Merian elements, the ancestral linkage groups of Lepidoptera, with the exception of three simple fusion events. We propose a standardized approach for genome taxonomy using nucleotide homology via scaffold chaining as the primary tool combined with Oxford plots based on Merian elements to infer and visualize directionality of chromosomal rearrangements. The identification of wing pattern genes promises future understanding of the evolution of forewing patterns in the genus Hyles, although further sequencing data from more individuals are needed. The genomic data obtained provide additional reliable references for further comparative studies in hawkmoths (Sphingidae).

The Hox gene Abdominal-B regulates the appendage development during the embryogenesis of scorpionflies

The Homeotic Complex (Hox) genes encode conserved homeodomain transcription factors that specify segment identity and appendage morphology along the antero-posterior axis in bilaterian animals. The Hox gene Abdominal-B (Abd-B) is mainly expressed in the posterior segments of the abdomen and plays an important role in insect organogenesis. In Mecoptera, the potential function of this gene remains unclear yet. Here, we performed a de novo transcriptome assembly and identified an Abd-B ortholog in the scorpionfly Panorpa liui. Quantitative real-time reverse transcription PCR showed that Abd-B expression increased gradually in embryos 76 h post oviposition, and was mainly present in the more posterior abdominal segments. Embryonic RNA interference of Abd-B resulted in a set of abnormalities, including developmental arrest, malformed suckers and misspecification of posterior segment identity. These results suggest that Abd-B is required for the proper development of the posterior abdomen. Furthermore, in Abd-B RNAi embryos, the expression of the appendage marker Distal-less (Dll) was up-regulated and was additionally present on abdominal segments IX and X compared with wild embryos, suggesting that scorpionfly Abd-B may act to suppress proleg development and has gained the ability to repress Dll expression on the more posterior abdominal segments. This study provides additional information on both the functional and evolutionary roles of Abd-B across different insects.

Distal-less and spalt are distal organisers of pierid wing patterns

Two genes, Distal-less (Dll) and spalt (sal), are known to be involved in establishing nymphalid butterfly wing patterns. They function in several ways: in the differentiation of the eyespot’s central signalling cells, or foci; in the differentiation of the surrounding black disc; in overall scale melanisation (Dll); and in elaborating marginal patterns, such as parafocal elements. However, little is known about the functions of these genes in the development of wing patterns in other butterfly families. Here, we study the expression and function of Dll and sal in the development of spots and other melanic wing patterns of the Indian cabbage white, Pieris canidia, a pierid butterfly. In P. canidia, both Dll and Sal proteins are expressed in the scale-building cells at the wing tips, in chevron patterns along the pupal wing margins, and in areas of future scale melanisation. Additionally, Sal alone is expressed in the future black spots. CRISPR knockouts of Dll and sal showed that each gene is required for the development of melanic wing pattern elements, and repressing pteridine granule formation, in the areas where they are expressed. We conclude that both genes likely play ancestral roles in organising distal butterfly wing patterns, across pierid and nymphalid butterflies, but are unlikely to be differentiating signalling centres in pierids black spots. The genetic and developmental mechanisms that set up the location of spots and eyespots are likely distinct in each lineage.

Transcription factors underlying wing margin color patterns and pupal cuticle markings in butterflies

Background

The diversity of butterfly color patterns can be attributed to a relatively small number of pattern elements that are homologous across Lepidoptera. Although genes involved in patterning some of these elements have been identified, the development of several major elements remains poorly understood. To identify genes underlying wing pupal cuticle markings and wing margin color patterns, we examined expression of the candidate transcription factors Engrailed/Invected (En/Inv), Distal-less (Dll), Cubitus interruptus (Ci), and Spalt in two nymphalids: Junonia coenia and Bicyclus anynana.

Results

We found that En/Inv, Dll, and Ci mark domains on the J. coenia last-instar forewing disc that closely correspond to the position and shape of pupal cuticle markings. We also found that Spalt demarcates wing margin color patterns in both J. coenia and B. anynana, and that CRISPR/Cas9 deletions in the spalt gene result in reduction and loss of wing margin color patterns in J. coenia. These data demonstrate a role for spalt in promoting wing margin color patterning, in addition to its previously described role in eyespot patterning.

Conclusion

Our observations support the model that a core set of regulatory genes are redeployed multiple times, and in multiple roles, during butterfly wing pattern development. Of these genes, spalt is of special interest as it plays a dual role in both eyespot and margin color pattern development.

Activation of butterfly eyespots by Distal-less is consistent with a reaction-diffusion process

Eyespots on the wings of nymphalid butterflies represent colorful examples of pattern formation, yet the developmental origins and mechanisms underlying eyespot center differentiation are still poorly understood. Using CRISPR-Cas9 we re-examine the function of Distal-less (Dll) as an activator or repressor of eyespots, a topic that remains controversial. We show that the phenotypic outcome of CRISPR mutations depends upon which specific exon is targeted. In Bicyclus anynana, exon 2 mutations are associated with both missing and ectopic eyespots, and also exon skipping. Exon 3 mutations, which do not lead to exon skipping, produce only null phenotypes, including missing eyespots, lighter wing coloration and loss of scales. Reaction-diffusion modeling of Dll function, using Wnt and Dpp as candidate morphogens, accurately replicates these complex crispant phenotypes. These results provide new insight into the function of Dll as a potential activator of eyespot development, scale growth and melanization, and suggest that the tuning of Dll expression levels can generate a diversity of eyespot phenotypes, including their appearance on the wing.This article has an associated 'The people behind the papers' interview.

UTR-specific knockdown of Distal-less and Sp8 leads to new phenotypic variants in the flour beetle Tribolium

RNA interference (RNAi)-mediated knockdown serves as an effective technique for the functional analysis of developmental genes that is well established in many organisms. In the beetle Tribolium castaneum, double-stranded RNA is applied by simple injection and distributes systemically within the tissue. Thus, systematic testing for RNAi specificity and efficiency is easily possible in this organism. Generally, the use of non-overlapping dsRNA fragments yielding qualitatively identical phenotypes is the method of choice to verify target-specific knockdown effects. Here, we show that UTR-specific RNAi results in different effects regarding quality, severity and penetrance when compared to RNAi fragments directed at the coding region. Furthermore, when using 3'UTR-specific dsRNA, we first describe the Distal-less^RNAi antenna-to-leg transformation phenotype in the Tribolium larva, which has only been observed in the adult beetle and Drosophila so far. In addition, we unexpectedly observed sterility effects caused by 3'UTR-specific knockdown of the Tribolium-Sp8 orthologue that is not seen when dsRNA targeted a sequence within the coding-region or the 5'UTR that itself led to early embryonic lethality. We conclude that targeting UTR sequences by region-specific RNAi can reveal unexpected new aspects of gene function applicable in basic research and crop protection.

Distal-less homeobox genes of insects and spiders: genomic organization, function, regulation and evolution

The Distal-less (Dll) genes are homeodomain transcription factors that are present in most Metazoa and in representatives of all investigated arthropod groups. In Drosophila, the best studied insect, Dll plays an essential role in forming the proximodistal axis of the legs, antennae and analia, and in specifying antennal identity. The initiation of Dll expression in clusters of cells in mid-lateral regions of the Drosophila embryo represents the earliest genetic marker of limbs. Dll genes are involved in the development of the peripheral nervous system and sensitive organs, and they also function as master regulators of black pigmentation in some insect lineages. Here we analyze the complete genomes of six insects, the nematode Caenorhabditis elegans and Homo sapiens, as well as multiple Dll sequences available in databases in order to examine the structure and protein features of these genes. We also review the function, expression, regulation and evolution of arthropod Dll genes with emphasis on insects and spiders.

Formation of body appendages during caudal regeneration in Platynereis dumerilii: adaptation of conserved molecular toolsets

Background

Platynereis and other polychaete annelids with homonomous segmentation are regarded to closely resemble ancestral forms of bilateria. The head region comprises the prostomium, the peristomium, a variable number of cephalized body segments and several appendages, like cirri, antennae and palps. The trunk of such polychaetes shows numerous, nearly identical segments. Each segment bears a parapodium with species-specific morphology on either side. The posterior end of the trunk features a segment proliferation zone and a terminal pygidium with the anus and anal cirri. The removal of a substantial part of the posterior trunk is by no means lethal. Cells at the site of injury dedifferentiate and proliferate forming a blastema to regenerate both the pygidium and the proliferation zone. The pygidium forms new anal cirri, and the proliferation zone generates new segments at a rapid pace. The formation of body appendages like the cirri and the segmental parapodia can thus be studied in the caudal regenerate of Platynereis within only a few days.

Results

The development of body appendages in Platynereis is regulated by a network of genes common to polychaetes but also shared by distant taxa. We isolated DNA sequences from P. dumerilii of five genes known to be involved in appendage formation within other groups: Meis/homothorax, Pbx1/extradenticle, Dlx/Distal-less, decapentaplegic and specificprotein1/buttonhead. Analyses of expression patterns during caudal regeneration by in situ hybridization reveal striking similarities related to expression in arthropods and vertebrates. All genes exhibit transient expression during differentiation and growth of segments. As was shown previously in other phyla Pdu-Meis/hth and Pdu-Pbx1/exd are co-expressed, although the expression is not limited to the proximal part of the parapodia. Pdu-Dll is prominent in parapodia but upregulated in the anal cirri. No direct dependence concerning Pdu-Dll and Pdu-sp/btd expression is observed in Platynereis. Pdu-dpp shows an expression pattern not comparable to its expression in other taxa.

Conclusions

The expression patterns observed suggest conserved roles of these genes during appendage formation across different clades, but the underlying mechanisms utilizing this toolset might not be identical. Some genes show broad expression along the proximodistal axis indicating a possible role in proximodistal patterning of body appendages. Other genes exhibit expression patterns limited to specific parts and tissues of the growing parapodia, thus presumably being involved in formation of taxon-specific morphological differences.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-016-0046-6) contains supplementary material, which is available to authorized users.

A Single-Wing Removal Method to Assess Correspondence Between Gene Expression and Phenotype in Butterflies: The Case of Distal-less

It is often desirable but difficult to retrieve information on the mature phenotype of an immature tissue sample that has been subjected to gene expression analysis. This problem cannot be ignored when individual variation within a species is large. To circumvent this problem in the butterfly wing system, we developed a new surgical method for removing a single forewing from a pupa using Junonia orithya; the operated pupa was left to develop to an adult without eclosion. The removed right forewing was subjected to gene expression analysis, whereas the non-removed left forewing was examined for color patterns. As a test case, we focused on Distal-less (Dll), which likely plays an active role in inducing elemental patterns, including eyespots. The Dll expression level in forewings was paired with eyespot size data from the same individual. One third of the operated pupae survived and developed wing color patterns. Dll expression levels were significantly higher in males than in females, although male eyespots were smaller in size than female eyespots. Eyespot size data showed weak but significant correlations with the Dll expression level in females. These results demonstrate that a single-wing removal method was successfully applied to the butterfly wing system and suggest the weak and non-exclusive contribution of Dll to eyespot size determination in this butterfly. Our novel methodology for establishing correspondence between gene expression and phenotype can be applied to other candidate genes for color pattern development in butterflies. Conceptually similar methods may also be applicable in other developmental systems.

Distal-less induces elemental color patterns in Junonia butterfly wings

BACKGROUND

The border ocellus, or eyespot, is a conspicuous color pattern element in butterfly wings. For two decades, it has been hypothesized that transcription factors such as Distal-less (Dll) are responsible for eyespot pattern development in butterfly wings, based on their expression in the prospective eyespots. In particular, it has been suggested that Dll is a determinant for eyespot size. However, functional evidence for this hypothesis has remained incomplete, due to technical difficulties.

RESULTS

Here, we show that ectopically expressed Dll induces ectopic elemental color patterns in the adult wings of the blue pansy butterfly, Junonia orithya (Lepidoptera, Nymphalidae). Using baculovirus-mediated gene transfer, we misexpressed Dll protein fused with green fluorescent protein (GFP) in pupal wings, resulting in ectopic color patterns, but not the formation of intact eyespots. Induced changes included clusters of black and orange scales (a basic feature of eyespot patterns), black and gray scales, and inhibition of cover scale development. In contrast, ectopic expression of GFP alone did not induce any color pattern changes using the same baculovirus-mediated gene transfer system.

CONCLUSIONS

These results suggest that Dll plays an instructive role in the development of color pattern elements in butterfly wings, although Dll alone may not be sufficient to induce a complete eyespot. This study thus experimentally supports the hypothesis of Dll function in eyespot development.

Molecular phylogeny of four homeobox genes from the purple sea star Pisaster ochraceus

Homeobox genes cloned from the purple sea star Pisaster ochraceus (Phylum Echinodermata/Class Asteroidea) were used along with related sequences available from members of other representative animal phyla to generate molecular phylogenies for Distal-less/Dlx, Hox5, Hox7, and Hox9/10 homeobox genes. Phylogenetic relationships were inferred based on the predicted 60 amino acid homeodomain, using amino acid (AA) and nucleotide (NT) models as well as the recently developed codon substitution models of sequence evolution. The resulting phylogenetic trees were mostly congruent with the consensus species-tree, grouping these newly identified genes with those isolated from other Asteroidea. This analysis also allowed a preliminary comparison of the performance of codon models with that of NT and AA evolutionary models in the inference of homeobox phylogeny. We found that, overall, the NT models displayed low reliability in recovering major clades at the Superphylum/Phylum level, and that codon models were slightly more dependable than AA models. Remarkably, in the majority of cases, codon substitution models seemed to outperform both AA and NT models at both the Class level and homeobox paralogy-group level of classification.