Genomic sequence around butterfly wing development genes: annotation and comparative analysis

apterous A specifies dorsal wing patterns and sexual traits in butterflies

Butterflies have evolved different colour patterns on their dorsal and ventral wing surfaces to serve different signalling functions, yet the developmental mechanisms controlling surface-specific patterning are still unknown. Here, we mutate both copies of the transcription factor apterous in Bicyclus anynana butterflies using CRISPR/Cas9 and show that apterous A, expressed dorsally, functions both as a repressor and modifier of ventral wing colour patterns, as well as a promoter of dorsal sexual ornaments in males. We propose that the surface-specific diversification of wing patterns in butterflies proceeded via the co-option of apterous A or its downstream effectors into various gene regulatory networks involved in the differentiation of discrete wing traits. Further, interactions between apterous and sex-specific factors such as doublesex may have contributed to the origin of sexually dimorphic surface-specific patterns. Finally, we discuss the evolution of eyespot number diversity in the family Nymphalidae within the context of developmental constraints due to apterous regulation.

Multifaceted biological insights from a draft genome sequence of the tobacco hornworm moth, Manduca sexta

Manduca sexta, known as the tobacco hornworm or Carolina sphinx moth, is a lepidopteran insect that is used extensively as a model system for research in insect biochemistry, physiology, neurobiology, development, and immunity. One important benefit of this species as an experimental model is its extremely large size, reaching more than 10 g in the larval stage. M. sexta larvae feed on solanaceous plants and thus must tolerate a substantial challenge from plant allelochemicals, including nicotine. We report the sequence and annotation of the M. sexta genome, and a survey of gene expression in various tissues and developmental stages. The Msex_1.0 genome assembly resulted in a total genome size of 419.4 Mbp. Repetitive sequences accounted for 25.8% of the assembled genome. The official gene set is comprised of 15,451 protein-coding genes, of which 2498 were manually curated. Extensive RNA-seq data from many tissues and developmental stages were used to improve gene models and for insights into gene expression patterns. Genome wide synteny analysis indicated a high level of macrosynteny in the Lepidoptera. Annotation and analyses were carried out for gene families involved in a wide spectrum of biological processes, including apoptosis, vacuole sorting, growth and development, structures of exoskeleton, egg shells, and muscle, vision, chemosensation, ion channels, signal transduction, neuropeptide signaling, neurotransmitter synthesis and transport, nicotine tolerance, lipid metabolism, and immunity. This genome sequence, annotation, and analysis provide an important new resource from a well-studied model insect species and will facilitate further biochemical and mechanistic experimental studies of many biological systems in insects.

Adaptive developmental plasticity: compartmentalized responses to environmental cues and to corresponding internal signals provide phenotypic flexibility

BACKGROUND

The environmental regulation of development can result in the production of distinct phenotypes from the same genotype and provide the means for organisms to cope with environmental heterogeneity. The effect of the environment on developmental outcomes is typically mediated by hormonal signals which convey information about external cues to the developing tissues. While such plasticity is a wide-spread property of development, not all developing tissues are equally plastic. To understand how organisms integrate environmental input into coherent adult phenotypes, we must know how different body parts respond, independently or in concert, to external cues and to the corresponding internal signals.

RESULTS

We quantified the effect of temperature and ecdysone hormone manipulations on post-growth tissue patterning in an experimental model of adaptive developmental plasticity, the butterfly Bicyclus anynana. Following a suite of traits evolving by natural or sexual selection, we found that different groups of cells within the same tissue have sensitivities and patterns of response that are surprisingly distinct for the external environmental cue and for the internal hormonal signal. All but those wing traits presumably involved in mate choice responded to developmental temperature and, of those, all but the wing traits not exposed to predators responded to hormone manipulations. On the other hand, while patterns of significant response to temperature contrasted traits on autonomously-developing wings, significant response to hormone manipulations contrasted neighboring groups of cells with distinct color fates. We also showed that the spatial compartmentalization of these responses cannot be explained by the spatial or temporal compartmentalization of the hormone receptor protein.

CONCLUSIONS

Our results unravel the integration of different aspects of the adult phenotype into developmental and functional units which both reflect and impact evolutionary change. Importantly, our findings underscore the complexity of the interactions between environment and physiology in shaping the development of different body parts.

System-dependent regulations of colour-pattern development: a mutagenesis study of the pale grass blue butterfly

Developmental studies on wing colour patterns have been performed in nymphalid butterflies, but efficient genetic manipulations, including mutagenesis, have not been well established. Here, we have performed mutagenesis experiments in a lycaenid butterfly, the pale grass blue Zizeeria maha, to produce colour-pattern mutants. We fed the P-generation larvae an artificial diet containing the mutagen ethyl methane sulfonate (EMS), and the F₁- and F₂-generation adults showed various aberrant colour patterns: dorsoventral transformation, anterioposterior background colouration gap, weak contrast, disarrangement of spots, reduction of the size of spots, loss of spots, fusion of spots, and ectopic spots. Among them, the disarrangement, reduction, and loss of spots were likely produced by the coordinated changes of many spots of a single wing around the discal spot in a system-dependent manner, demonstrating the existence of the central symmetry system. The present study revealed multiple genetic regulations for system-dependent and wing-wide colour-pattern determination in lycaenid butterflies.

Constraint and contingency in multifunctional gene regulatory circuits

Gene regulatory circuits drive the development, physiology, and behavior of organisms from bacteria to humans. The phenotypes or functions of such circuits are embodied in the gene expression patterns they form. Regulatory circuits are typically multifunctional, forming distinct gene expression patterns in different embryonic stages, tissues, or physiological states. Any one circuit with a single function can be realized by many different regulatory genotypes. Multifunctionality presumably constrains this number, but we do not know to what extent. We here exhaustively characterize a genotype space harboring millions of model regulatory circuits and all their possible functions. As a circuit's number of functions increases, the number of genotypes with a given number of functions decreases exponentially but can remain very large for a modest number of functions. However, the sets of circuits that can form any one set of functions becomes increasingly fragmented. As a result, historical contingency becomes widespread in circuits with many functions. Whether a circuit can acquire an additional function in the course of its evolution becomes increasingly dependent on the function it already has. Circuits with many functions also become increasingly brittle and sensitive to mutation. These observations are generic properties of a broad class of circuits and independent of any one circuit genotype or phenotype.

Taxonomy and distribution pattern of the African rain forest butterfly genus Euphaedra Hübner sensu stricto with the description of three new subspecies of Euphaedra cyparissa (Cramer) and one of E. sarcoptera (Butler) (Lepidoptera, Nymphalidae, Limenitidinae, Adoliadini)

Updated data on the distribution, ecology and taxonomy of Euphaedra cyparissa (Cramer) and Euphaedra sarcoptera (Butler) are presented. Three new subspecies of Euphaedra cyparissa and one of Euphaedra sarcoptera are described and their geographic distribution is presented. The monophyly of the genus Euphaedra sensu Hecq is assessed based on morphological, in particular male and female genitalia, and behavioural traits. Possible evolutionary reasons for the convergence of colour pattern between the sympatric subspecies of Euphaedra cyparissa and Euphaedra sarcoptera are discussed.

Linkage map of the peppered moth, Biston betularia (Lepidoptera, Geometridae): a model of industrial melanism

We have constructed a linkage map for the peppered moth (Biston betularia), the classical ecological genetics model of industrial melanism, aimed both at localizing the network of loci controlling melanism and making inferences about chromosome dynamics. The linkage map, which is based primarily on amplified fragment length polymorphisms (AFLPs) and genes, consists of 31 linkage groups (LGs; consistent with the karyotype). Comparison with the evolutionarily distant Bombyx mori suggests that the gene content of chromosomes is highly conserved. Gene order is conserved on the autosomes, but noticeably less so on the Z chromosome, as confirmed by physical mapping using bacterial artificial chromosome fluorescence in situ hybridization (BAC-FISH). Synteny mapping identified three pairs of B. betularia LGs (11/29, 23/30 and 24/31) as being orthologous to three B. mori chromosomes (11, 23 and 24, respectively). A similar finding in an outgroup moth (Plutella xylostella) indicates that the B. mori karyotype (n=28) is a phylogenetically derived state resulting from three chromosome fusions. As with other Lepidoptera, the B. betularia W chromosome consists largely of repetitive sequence, but exceptionally we found a W homolog of a Z-linked gene (laminin A), possibly resulting from ectopic recombination between the sex chromosomes. The B. betularia linkage map, featuring the network of known melanization genes, serves as a resource for melanism research in Lepidoptera. Moreover, its close resemblance to the ancestral lepidopteran karyotype (n=31) makes it a useful reference point for reconstructing chromosome dynamic events and ancestral genome architectures. Our study highlights the unusual evolutionary stability of lepidopteran autosomes; in contrast, higher rates of intrachromosomal rearrangements support a special role of the Z chromosome in adaptive evolution and speciation.

Identification of microRNAs in Helicoverpa armigera and Spodoptera litura based on deep sequencing and homology analysis

The current identification of microRNAs (miRNAs) in insects is largely dependent on genome sequences. However, the lack of available genome sequences inhibits the identification of miRNAs in various insect species. In this study, we used a miRNA database of the silkworm Bombyx mori as a reference to identify miRNAs in Helicoverpa armigera and Spodoptera litura using deep sequencing and homology analysis. Because all three species belong to the Lepidoptera, the experiment produced reliable results. Our study identified 97 and 91 conserved miRNAs in H. armigera and S. litura, respectively. Using the genome of B. mori and BAC sequences of H. armigera as references, 1 novel miRNA and 8 novel miRNA candidates were identified in H. armigera, and 4 novel miRNA candidates were identified in S. litura. An evolutionary analysis revealed that most of the identified miRNAs were insect-specific, and more than 20 miRNAs were Lepidoptera-specific. The investigation of the expression patterns of miR-2a, miR-34, miR-2796-3p and miR-11 revealed their potential roles in insect development. miRNA target prediction revealed that conserved miRNA target sites exist in various genes in the 3 species. Conserved miRNA target sites for the Hsp90 gene among the 3 species were validated in the mammalian 293T cell line using a dual-luciferase reporter assay. Our study provides a new approach with which to identify miRNAs in insects lacking genome information and contributes to the functional analysis of insect miRNAs.

Comprehensive microarray-based analysis for stage-specific larval camouflage pattern-associated genes in the swallowtail butterfly, Papilio xuthus

Background

Body coloration is an ecologically important trait that is often involved in prey-predator interactions through mimicry and crypsis. Although this subject has attracted the interest of biologists and the general public, our scientific knowledge on the subject remains fragmentary. In the caterpillar of the swallowtail butterfly Papilio xuthus, spectacular changes in the color pattern are observed; the insect mimics bird droppings (mimetic pattern) as a young larva, and switches to a green camouflage coloration (cryptic pattern) in the final instar. Despite the wide variety and significance of larval color patterns, few studies have been conducted at a molecular level compared with the number of studies on adult butterfly wing patterns.

Results

To obtain a catalog of genes involved in larval mimetic and cryptic pattern formation, we constructed expressed sequence tag (EST) libraries of larval epidermis for P. xuthus, and P. polytes that contained 20,736 and 5,376 clones, respectively, representing one of the largest collections available in butterflies. A comparison with silkworm epidermal EST information revealed the high expression of putative blue and yellow pigment-binding proteins in Papilio species. We also designed a microarray from the EST dataset information, analyzed more than five stages each for six markings, and confirmed spatial expression patterns by whole-mount in situ hybridization. Hence, we succeeded in elucidating many novel marking-specific genes for mimetic and cryptic pattern formation, including pigment-binding protein genes, the melanin-associated gene yellow-h3, the ecdysteroid synthesis enzyme gene 3-dehydroecdysone 3b-reductase, and Papilio-specific genes. We also found many cuticular protein genes with marking specificity that may be associated with the unique surface nanostructure of the markings. Furthermore, we identified two transcription factors, spalt and ecdysteroid signal-related E75, as genes expressed in larval eyespot markings. This finding suggests that E75 is a strong candidate mediator of the hormone-dependent coordination of larval pattern formation.

Conclusions

This study is one of the most comprehensive molecular analyses of complicated morphological features, and it will serve as a new resource for studying insect mimetic and cryptic pattern formation in general. The wide variety of marking-associated genes (both regulatory and structural genes) identified by our screening indicates that a similar strategy will be effective for understanding other complex traits.

Genetic basis of stage-specific melanism: a putative role for a cysteine sulfinic acid decarboxylase in insect pigmentation

Melanism, the overall darkening of the body, is a widespread form of animal adaptation to particular environments, and includes bookcase examples of evolution by natural selection, such as industrial melanism in the peppered moth. The major components of the melanin biosynthesis pathway have been characterized in model insects, but little is known about the genetic basis of life-stage specific melanism such as cases described in some lepidopteran species. Here, we investigate two melanic mutations of Bicyclus anynana butterflies, called Chocolate and melanine, that exclusively affect pigmentation of the larval and adult stages, respectively. Our analysis of Mendelian segregation patterns reveals that the larval and adult melanic phenotypes are due to alleles at different, independently segregating loci. Our linkage mapping analysis excludes the pigmentation candidate gene black as the melanine locus, and implicates a gene encoding a putative pyridoxal phosphate-dependant cysteine sulfinic acid decarboxylase as the Chocolate locus. We show variation in coding sequence and in expression levels for this candidate larval melanism locus. This is the first study that suggests a biological function for this gene in insects. Our findings open up exciting opportunities to study the role of this locus in the evolution of adaptive variation in pigmentation, and the uncoupling of regulation of pigment biosynthesis across developmental stages with different ecologies and pressures on body coloration.