Matthew L. Aardema

Parallel Molecular Evolution in an Herbivore Community

Making the Change The Apocynaceae plant family produces toxic cardenolids. However, many insects have managed to escape the deleterious effects of these chemicals and even, in some cases, use them in their own defense. Zhen et al. (p. 1634) investigated a broad range of taxa and found several examples of parallel changes, as well as duplications, in members of the ATPα family that likely explain the shift to allow insects to avoid the toxic effects of these plants. Thus, natural selection can harness a combination of gene duplication, protein evolution, and regulatory evolution to allow distantly related species to adapt to specific niches. Parallel mutations in the alpha subunit of the sodium pump allow insects to specialize on host plants that produce ouabain. Numerous insects have independently evolved the ability to feed on plants that produce toxic secondary compounds called cardenolides and can sequester these compounds for use in their defense. We surveyed the protein target for cardenolides, the alpha subunit of the sodium pump, Na+,K+-ATPase (ATPα), in 14 species that feed on cardenolide-producing plants and 15 outgroups spanning three insect orders. Despite the large number of potential targets for modulating cardenolide sensitivity, amino acid substitutions associated with host-plant specialization are highly clustered, with many parallel substitutions. Additionally, we document four independent duplications of ATPα with convergent tissue-specific expression patterns. We find that unique substitutions are disproportionately associated with recent duplications relative to parallel substitutions. Together, these findings support the hypothesis that adaptation tends to take evolutionary paths that minimize negative pleiotropy.

Analysis of the Population Structure of Anaplasma phagocytophilum Using Multilocus Sequence Typing

Anaplasma phagocytophilum is a Gram-negative obligate intracellular bacterium that replicates in neutrophils. It is transmitted via tick-bite and causes febrile disease in humans and animals. Human granulocytic anaplasmosis is regarded as an emerging infectious disease in North America, Europe and Asia. However, although increasingly detected, it is still rare in Europe. Clinically apparent A. phagocytophilum infections in animals are mainly found in horses, dogs, cats, sheep and cattle. Evidence from cross-infection experiments that A. phagocytophilum isolates of distinct host origin are not uniformly infectious for heterologous hosts has led to several approaches of molecular strain characterization. Unfortunately, the results of these studies are not always easily comparable, because different gene regions and fragment lengths were investigated. Multilocus sequence typing is a widely accepted method for molecular characterization of bacteria. We here provide for the first time a universal typing method that is easily transferable between different laboratories. We validated our approach on an unprecedented large data set of almost 400 A. phagocytophilum strains from humans and animals mostly from Europe. The typability was 74% (284/383). One major clonal complex containing 177 strains was detected. However, 54% (49/90) of the sequence types were not part of a clonal complex indicating that the population structure of A. phagocytophilum is probably semiclonal. All strains from humans, dogs and horses from Europe belonged to the same clonal complex. As canine and equine granulocytic anaplasmosis occurs frequently in Europe, human granulocytic anaplasmosis is likely to be underdiagnosed in Europe. Further, wild boars and hedgehogs may serve as reservoir hosts of the disease in humans and domestic animals in Europe, because their strains belonged to the same clonal complex. In contrast, as they were only distantly related, roe deer, voles and shrews are unlikely to harbor A. phagocytophilum strains infectious for humans, domestic or farm animals.

Allochronic isolation and incipient hybrid speciation in tiger swallowtail butterflies

Hybridization leading to reproductively isolated, novel genotypes is poorly understood as a means of speciation and few empirical examples have been studied. In 1999, a previously non-existent delayed flight of what appeared to be the Canadian tiger swallowtail butterfly, Papilio canadensis, was observed in the Battenkill River Valley, USA. Allozyme frequencies and morphology suggest that this delayed flight was the product of hybridization between Papilio canadensis and its sibling species Papilio glaucus. The mitochondrial DNA (mtDNA) restriction fragment length polymorphisms presented here indicate that only P. canadensis-like mtDNA occurs in this population, suggesting that introgression likely occurred from hybrid males mating with P. canadensis females. Preliminary studies of this population indicated that delayed post-diapause pupal emergence in this hybrid genotype was the root cause behind the observed delayed flight, which suggests a potential empirical example of a mechanism leading to reproductive isolation. Here we provide further evidence of the role of adult pupal emergence as a reproductive barrier likely leading to reproductive isolation. In particular, we present results from pupal emergence studies using four different spring and two different winter temperature treatments. The results indicate a clear separation of adult emergences between the hybrid population and both parental species. However, our results indicate that exceptionally hot springs are likely to lead to greater potential for overlap between the local parental species, P. canadensis, and this delayed population with hybrid origins. Conversely, our results also show that warmer winters are likely to increase the temporal separation of the hybrid population and the parental species. Finally, we report recently collected evidence that this hybrid population remains morphologically distinct.

Hybridization leads to host-use divergence in a polyphagous butterfly sibling species pair

Climate warming has lead to increased genetic introgression across a narrow hybrid zone separating the eastern and Canadian tiger swallowtails (Papilio glaucus and Papilio canadensis). This situation has led to the formation of an allochronically separated hybrid population with a delayed emerging phenotype or “late flight”. Here, we assess how the recombination of the parental genomes that lead to this phenotype may have facilitated another major ecological shift, host-use divergence. We first contrast the ovipositional profiles of the late flight population to that of the parental species P. glaucus and P. canadensis. Subsequently we contrast the larval survival and growth of the late flight, a P. canadensis and a P. glaucus population, and a population from the northern edge of the hybrid zone on five hosts. Our results indicate that the ovipositional preference of this hybrid swarm is identical to that of the introgressing parental species, P. glaucus. Due to the absence of the preferred hosts of P. glaucus (Liriodendron tulipifera L. and Ptelea trifoliata L.) where the late flight occurs, this ovipositional pattern implies a functional specialization onto a secondary host of both parental species, Fraxinus americana L. In contrast, the larval host-use abilities represent a mixture of P. glaucus and P. canadensis, indicating divergence in larval host-use abilities has not taken place. However, high genetic variability (genetic coefficient of variation) is present for growth on F. americana in the late flight hybrid swarm and tradeoffs for larval performance on the preferred hosts of the parental species are evident; indicating a strong potential for future specialization in larval host-use abilities. This current scenario represents an instance where a shift in a major ecological trait, host use, is likely occurring as a byproduct of a shift in an unrelated trait (delayed emergence) leading to partial reproductive isolation.

The evolution of cardenolide‐resistant forms of Na+,K+‐ATPase in Danainae butterflies

Cardenolides are a class of plant secondary compounds that inhibit the proper functioning of the Na+, K+‐ATPase enzyme in susceptible animals. Nonetheless, many insect species are able to sequester cardenolides for their own defence. These include butterflies in the subfamily Danainae (Family: Nymphalidae) such as the monarch (Danaus plexippus). Previous studies demonstrated that monarchs harbour an asparagine (N) to histidine (H) substitution (N122H) in the α subunit of Na+, K+‐ATPase (ATPα) that reduces this enzyme’s sensitivity to cardenolides. More recently, it has been suggested that at ATPα position 111, monarchs may also harbour a leucine (L)/glutamine (Q) polymorphism. This later amino acid could also contribute to cardenolide insensitivity. However, here we find that incorrect annotation of the initially reported DNA sequence for ATPα has led to several erroneous conclusions. Using a population genetic and phylogenetic analysis of monarchs and their close relatives, we show that an ancient Q111L substitution occurred prior to the radiation of all Danainae, followed by a second substitution at the same site to valine (V), which arose before the diversification of the Danaus genus. In contrast, N122H appears to be a recent substitution specific to monarchs. Surprisingly, examination of a broader insect phylogeny reveals that the same progression of amino acid substitutions (Q111L → L111V + N122H) has also occurred in Chyrsochus beetles (Family: Chrysomelidae, Subfamily: Eumolpinae) that feed on cardenolide‐containing host plants. The parallel pattern of amino acid substitution in these two distantly related lineages is consistent with an adaptive role for these substitutions in reducing cardenolide sensitivity and suggests that their temporal order may be limited by epistatic interactions.

Considering Local Adaptation in Issues of Lepidopteran Conservation—a Review and Recommendations

Abstract Butterflies are among the most charismatic and alluring of invertebrate organisms. Their complex life histories and diverse habitat affiliations often produce regional variation and local adaptation among different populations. Unfortunately, numerous lepidopteran species are facing extinction. In this paper we address the importance of considering local adaptation in butterfly conservation and encourage investigations of regional specialization that may enhance the effectiveness of conservation strategies such as captive rearing, habitat restoration and the introduction of populations into new localities. We also give suggestions for dealing with problems arising from local adaptation in particular instances.

Differential effects of short term winter thermal stress on diapausing tiger swallowtail butterflies (Papilio spp.)

Abstract  It is generally thought that insects inhabiting lower latitudes are more severely impacted by changes in their thermal environment than are high latitude species. This is attributed to the wider range of temperatures to which high‐latitude species are exposed. By contrast, low‐latitude species have typically evolved in more thermally stable environments with a narrower range of temperature variation. However, deviation from this pattern can occur and here we report that under variable winter conditions a higher latitude species may be more sensitive to thermal variation than its lower latitude sister species. Using split broods, we examined the survival and adult emergence success of diapausing pupae of Papilio canadensis and P. glaucus, as well as a unique, recombinant hybrid population (“late‐flight”) to short periods of mid‐winter cold and heat stress. Our results indicate that the higher latitude, univoltine populations (P. canadensis and late‐flights) exhibit lower pupal survival than the lower latitude, facultative diapauser (P. glaucus) for all mid‐winter thermal stress treatments, both high and low. Size differences alone do not appear to account for the observed differences in survival or metabolic costs in these three phenotypes, as late‐flight individuals are similar in size to P. glaucus. We attribute the observed differences in survival and weight loss to potential metabolic differences and variation in the intensity of diapause, in addition to divergent adaptation to winter precipitation levels (e.g. snow cover) and the influences this may have on microhabitat temperature moderation.

Sexual Dimorphism and Retinal Mosaic Diversification Following the Evolution of a Violet Receptor in Butterflies.

Numerous animal lineages have expanded and diversified the opsin-based photoreceptors in their eyes underlying color vision behavior. However, the selective pressures giving rise to new photoreceptors and their spectral tuning remain mostly obscure. Previously, we identified a violet receptor (UV2) that is the result of a UV opsin gene duplication specific to Heliconius butterflies. At the same time the violet receptor evolved, Heliconius evolved UV-yellow coloration on their wings, due to the pigment 3-hydroxykynurenine (3-OHK) and the nanostructure architecture of the scale cells. In order to better understand the selective pressures giving rise to the violet receptor, we characterized opsin expression patterns using immunostaining (14 species) and RNA-Seq (18 species), and reconstructed evolutionary histories of visual traits in five major lineages within Heliconius and one species from the genus Eueides. Opsin expression patterns are hyperdiverse within Heliconius. We identified six unique retinal mosaics and three distinct forms of sexual dimorphism based on ommatidial types within the genus Heliconius. Additionally, phylogenetic analysis revealed independent losses of opsin expression, pseudogenization events, and relaxation of selection on UVRh2 in one lineage. Despite this diversity, the newly evolved violet receptor is retained across most species and sexes surveyed. Discriminability modeling of behaviorally preferred 3-OHK yellow wing coloration suggests that the violet receptor may facilitate Heliconius color vision in the context of conspecific recognition. Our observations give insights into the selective pressures underlying the origins of new visual receptors.

Phylogenetic incongruence and the evolutionary origins of cardenolide‐resistant forms of Na+,K+‐ATPase in Danaus butterflies

Many distantly related insect species are specialized feeders of cardenolide‐containing host plants such as milkweed (Asclepias spp.). Previous studies have revealed frequent, parallel substitution of a functionally important amino acid substitution (N122H) in the alpha subunit of Na+,K+‐ATPase in a number of these species. This substitution facilitates the ability of these insects to feed on their toxic hosts and sequester cardenolides for their own use in defense. Among milkweed butterflies of the genus Danaus, the previously established phylogeny for this group suggests that N122H arose independently and fixed in two distinct lineages. We reevaluate this conclusion by examining Danaus phylogenetic relationships using >400 orthologous gene sequences assembled from transcriptome data. Our results indicate that the three Danaus species known to harbor the N122H substitution are more closely related than previously thought, consistent with a single, common origin for N122H. However, we also find evidence of both incomplete lineage sorting and post‐speciation genetic exchange among these butterfly species, raising the possibility of collateral evolution of cardenolide‐insensitivity in this species group.

The role of latitudinal, genetic and temperature variation in the induction of diapause of Papilio glaucus (Lepidoptera: Papilionidae)

A key adaptation in insects for dealing with variable environmental conditions is the ability to diapause. The tiger swallowtail butterflies, Papilio glaucus and P. canadensis are ideal species to explore the genetic causes and population genetic consequences of diapause because divergence in this trait is believed to be a salient factor in maintaining a hybrid zone between these species. Yet little is known about the factors that influence diapause induction in this system. Here we explored how spatial (latitudinal), environmental (temperature) and genetic (hybridization) factors affect diapause induction in this system. Specifically, a series of growth chamber experiments using wild caught individuals from across the eastern United States were performed to: (1) evaluate how critical photoperiod varies with latitude, (2) isolate the stage in which induction occurs, (3) test whether changes in temperature affected rates of diapause induction, and (4) explore how the incidence of diapause is affected in hybrid offspring. We find that induction occurs in the larval stage, is not sensitive to a relatively broad range of temperatures, appears to have a complex genetic basis (i.e., is not simply a dominant trait following a Mendelian inheritance pattern) and that the critical photoperiod increases by 0.4 h with each increasing degree in latitude. This work deepens our understanding of how spatial, environmental and genetic variation influences a key seasonal adaptation (diapause induction) in a well‐developed ecological model system and will make possible future studies that explore how climatic variation affects the population dynamics and genetics of this system.