skip to main content


This content will become publicly available on May 24, 2024

Title: Genome sequence of Ophryocystis elektroscirrha, an apicomplexan parasite of monarch butterflies: cryptic diversity and response to host-sequestered plant chemicals
Abstract

Apicomplexa are ancient and diverse organisms which have been poorly characterized by modern genomics. To better understand the evolution and diversity of these single-celled eukaryotes, we sequenced the genome ofOphryocystis elektroscirrha, a parasite of monarch butterflies,Danaus plexippus. We contextualize our newly generated resources within apicomplexan genomics before answering longstanding questions specific to this host-parasite system. To start, the genome is miniscule, totaling only 9 million bases and containing fewer than 3,000 genes, half the gene content of two other sequenced invertebrate-infecting apicomplexans,Porospora giganteaandGregarina niphandrodes. We found thatO. elektroscirrhashares different orthologs with each sequenced relative, suggesting the true set of universally conserved apicomplexan genes is very small indeed. Next, we show that sequencing data from other potential host butterflies can be used to diagnose infection status as well as to study diversity of parasite sequences. We recovered a similarly sized parasite genome from another butterfly,Danaus chrysippus, that was highly diverged from theO. elektroscirrhareference, possibly representing a distinct species. Using these two new genomes, we investigated potential evolutionary response by parasites to toxic phytochemicals their hosts ingest and sequester. Monarch butterflies are well-known to tolerate toxic cardenolides thanks to changes in the sequence of their Type II ATPase sodium pumps. We show thatOphryocystiscompletely lacks Type II or Type 4 sodium pumps, and related proteins PMCA calcium pumps show extreme sequence divergence compared to other Apicomplexa, demonstrating new avenues of research opened by genome sequencing of non-model Apicomplexa.

 
more » « less
Award ID(s):
1661454
NSF-PAR ID:
10481062
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
BioMed Central
Date Published:
Journal Name:
BMC Genomics
Volume:
24
Issue:
1
ISSN:
1471-2164
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Animals rely on a balance of endogenous and exogenous sources of immunity to mitigate parasite attack. Understanding how environmental context affects that balance is increasingly urgent under rapid environmental change. In herbivores, immunity is determined, in part, by phytochemistry which is plastic in response to environmental conditions. Monarch butterfliesDanaus plexippus, consistently experience infection by a virulent parasiteOphryocystis elektroscirrha, and some medicinal milkweed (Asclepias) species, with high concentrations of toxic steroids (cardenolides), provide a potent source of exogenous immunity.

    We investigated plant‐mediated influences of elevated CO2(eCO2) on endogenous immune responses of monarch larvae to infection byO. elektroscirrha. Recently, transcriptomics have revealed that infection byO. elektroscirrhadoes not alter monarch immune gene regulation in larvae, corroborating that monarchs rely more on exogenous than endogenous immunity. However, monarchs feeding on medicinal milkweed grown under eCO2lose tolerance to the parasite, associated with changes in phytochemistry. Whether changes in milkweed phytochemistry induced by eCO2alter the balance between exogenous and endogenous sources of immunity remains unknown.

    We fed monarchs two species of milkweed;A. curassavica(medicinal) andA. incarnata(non‐medicinal) grown under ambient CO2(aCO2) or eCO2. We then measured endogenous immune responses (phenoloxidase activity, haemocyte concentration and melanization strength), along with foliar chemistry, to assess mechanisms of monarch immunity under future atmospheric conditions.

    The melanization response of late‐instar larvae was reduced on medicinal milkweed in comparison to non‐medicinal milkweed. Moreover, the endogenous immune responses of early‐instar larvae to infection byO. elektroscirrhawere generally lower in larvae reared on foliage from aCO2plants and higher in larvae reared on foliage from eCO2plants. When grown under eCO2, milkweed plants exhibited lower cardenolide concentrations, lower phytochemical diversity and lower nutritional quality (higher C:N ratios). Together, these results suggest that the loss of exogenous immunity from foliage under eCO2results in increased endogenous immune function.

    Animal populations face multiple threats induced by anthropogenic environmental change. Our results suggest that shifts in the balance between exogenous and endogenous sources of immunity to parasite attack may represent an underappreciated consequence of environmental change.

     
    more » « less
  2. Abstract

    Herbivorous insects have evolved many mechanisms to overcome plant chemical defences, including detoxification and sequestration. Herbivores may also use toxic plants to reduce parasite infection. Plant toxins could directly interfere with parasites or could enhance endogenous immunity. Alternatively, plant toxins could favour down‐regulation of endogenous immunity by providing an alternative (exogenous) defence against parasitism. However, studies on genomewide transcriptomic responses to plant defences and the interplay between plant toxicity and parasite infection remain rare. Monarch butterflies (Danaus plexippus) are specialist herbivores of milkweeds (Asclepiasspp.), which contain toxic cardenolides. Monarchs have adapted to cardenolides through multiple resistance mechanisms and can sequester cardenolides to defend against bird predators. In addition, high‐cardenolide milkweeds confer monarch resistance to a specialist protozoan parasite (Ophryocystis elektroscirrha). We used this system to study the interplay between the effects of plant toxicity and parasite infection on global gene expression. We compared transcriptional profiles between parasite‐infected and uninfected monarch larvae reared on two milkweed species. Our results demonstrate that monarch differentially express several hundred genes when feeding onA. curassavicaandA. incarnata, two species that differ substantially in cardenolide concentrations. These differentially expressed genes include genes within multiple families of canonical insect detoxification genes, suggesting that they play a role in monarch toxin resistance and sequestration. Interestingly, we found little transcriptional response to infection. However, parasite growth was reduced in monarchs reared onA. curassavica, and in these monarchs, several immune genes were down‐regulated, consistent with the hypothesis that medicinal plants can reduce reliance on endogenous immunity.

     
    more » « less
  3. Abstract

    Insect–pathogen dynamics can show seasonal and inter‐annual variations that covary with fluctuations in insect abundance and climate. Long‐term analyses are especially needed to track parasite dynamics in migratory insects, in part because their vast habitat ranges and high mobility might dampen local effects of density and climate on infection prevalence.

    Monarch butterfliesDanaus plexippusare commonly infected with the protozoanOphryocystis elektroscirrha(OE). Because this parasite lowers monarch survival and flight performance, and because migratory monarchs have experienced declines in recent decades, it is important to understand the patterns and drivers of infection.

    Here we compiled data onOEinfection spanning 50 years, from wild monarchs sampled in the United States, Canada and Mexico during summer breeding, fall migrating and overwintering periods. We examined eastern versus western North American monarchs separately, to ask how abundance estimates, resource availability, climate and breeding season length impact infection trends. We further assessed the intensity of migratory culling, which occurs when infected individuals are removed from the population during migration.

    Average infection prevalence was four times higher in western compared to eastern subpopulations. In eastern North America, the proportion of infected monarchs increased threefold since the mid‐2000s. In the western region, the proportion of infected monarchs declined sharply from 2000 to 2015, and increased thereafter. For both eastern and western subpopulations, years with greater summer adult abundance predicted greater infection prevalence, indicating that transmission increases with host breeding density. Environmental variables (temperature and NDVI) were not associated with changes in the proportion of infected adults. We found evidence for migratory culling of infected butterflies, based on declines in parasitism during fall migration. We estimated that tens of millions fewer monarchs reach overwintering sites in Mexico as a result ofOE, highlighting the need to consider the parasite as a potential threat to the monarch population.

    Increases in infection among eastern North American monarchs post‐2002 suggest that changes to the host’s ecology or environment have intensified parasite transmission. Further work is needed to examine the degree to which human practices, such as mass caterpillar rearing and the widespread planting of exotic milkweed, have contributed to this trend.

     
    more » « less
  4. Abstract

    Humoral and cellular immune responses provide animals with major defences against harmful pathogens. While it is often assumed that immune genes undergo rapid diversifying selection, this assumption has not been tested in many species. Moreover, it is likely that different classes of immune genes experience different levels of evolutionary constraint, resulting in varying selection patterns. We examined the evolutionary patterns for a set of 91 canonical immune genes of North American monarch butterflies (Danaus plexippus), using as an outgroup the closely related soldier butterfly (Danaus eresimus). As a comparison to these immune genes, we selected a set of control genes that were paired with each immune for approximate size and genomic location. As a whole, these immune genes had a significant but modest reduction in Tajima'sDrelative to paired‐control genes, but otherwise did not show distinct patterns of population genetic variation or evolutionary rates. When further partitioning these immune genes into four functional classes (recognition, signalling, modulation, and effector), we found distinct differences among these groups. Relative to control genes, recognition genes exhibit increased nonsynonymous diversity and divergence, suggesting reduced constraints on evolution, and supporting the notion that coevolution with pathogens results in diversifying selection. In contrast, signalling genes showed an opposite pattern of reduced diversity and divergence, suggesting evolutionary constraints and conservation. Modulator and effector genes showed no statistical differences from controls. These results are consistent with patterns found in immune genes in fruit flies andPierisbutterflies, suggesting that consistent selective pressures on different classes of immune genes broadly govern the evolution of innate immunity among insects.

     
    more » « less
  5. Infections by maternally inherited bacterial endosymbionts, especially Wolbachia, are common in insects and other invertebrates but infec- tion dynamics across species ranges are largely under studied. Specifically, we lack a broad understanding of the origin of Wolbachia infec- tions in novel hosts, and the historical and geographical dynamics of infections that are critical for identifying the factors governing their spread. We used Genotype-by-Sequencing data from previous population genomics studies for range-wide surveys of Wolbachia pres- ence and genetic diversity in North American butterflies of the genus Lycaeides. As few as one sequence read identified by assembly to a Wolbachia reference genome provided high accuracy in detecting infections in host butterflies as determined by confirmatory PCR tests, and maximum accuracy was achieved with a threshold of only 5 sequence reads per host individual. Using this threshold, we detected Wolbachia in all but 2 of the 107 sampling localities spanning the continent, with infection frequencies within populations ranging from 0% to 100% of individuals, but with most localities having high infection frequencies (mean 1⁄4 91% infection rate). Three major lineages of Wolbachia were identified as separate strains that appear to represent 3 separate invasions of Lycaeides butterflies by Wolbachia. Overall, we found extensive evidence for acquisition of Wolbachia through interspecific transfer between host lineages. Strain wLycC was confined to a single butterfly taxon, hybrid lineages derived from it, and closely adjacent populations in other taxa. While the other 2 strains were detected throughout the rest of the continent, strain wLycB almost always co-occurred with wLycA. Our demographic modeling suggests wLycB is a recent invasion. Within strain wLycA, the 2 most frequent haplotypes are confined almost exclusively to separate butterfly taxa with haplotype A1 observed largely in Lycaeides melissa and haplotype A2 observed most often in Lycaeides idas localities, consistent with either cladogenic mode of infection acquisition from a common ancestor or by hybridization and accompanying mutation. More than 1 major Wolbachia strain was observed in 15 localities. These results demonstrate the utility of using resequencing data from hosts to quan- tify Wolbachia genetic variation and infection frequency and provide evidence of multiple colonizations of novel hosts through hybridiza- tion between butterfly lineages and complex dynamics between Wolbachia strains. 
    more » « less