skip to main content

Attention:

The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Thursday, October 10 until 2:00 AM ET on Friday, October 11 due to maintenance. We apologize for the inconvenience.


Title: Population genomics reveals variable patterns of immune gene evolution in monarch butterflies ( Danaus plexippus )
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
NSF-PAR ID:
10367338
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Molecular Ecology
Volume:
30
Issue:
18
ISSN:
0962-1083
Page Range / eLocation ID:
p. 4381-4391
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Plant pathogens use effector proteins to target host processes involved in pathogen perception, immune signalling, or defence outputs. Unlike foliar pathogens, it is poorly understood how root‐invading pathogens suppress immunity. The Avr2 effector from the tomato root‐ and xylem‐colonizing pathogenFusarium oxysporumsuppresses immune signalling induced by various pathogen‐associated molecular patterns (PAMPs). It is unknown how Avr2 targets the immune system. TransgenicAVR2 Arabidopsis thalianaphenocopies mutants in which the pattern recognition receptor (PRR) co‐receptor BRI1‐ASSOCIATED RECEPTOR KINASE (BAK1) or its downstream signalling kinase BOTRYTIS‐INDUCED KINASE 1 (BIK1) are knocked out. We therefore tested whether these kinases are Avr2 targets. Flg22‐induced complex formation of the PRR FLAGELLIN SENSITIVE 2 and BAK1 occurred in the presence and absence of Avr2, indicating that Avr2 does not affect BAK1 function or PRR complex formation. Bimolecular fluorescence complementation assays showed that Avr2 and BIK1 co‐localize in planta. Although Avr2 did not affect flg22‐induced BIK1 phosphorylation, mono‐ubiquitination was compromised. Furthermore, Avr2 affected BIK1 abundance and shifted its localization from nucleocytoplasmic to the cell periphery/plasma membrane. Together, these data imply that Avr2 may retain BIK1 at the plasma membrane, thereby suppressing its ability to activate immune signalling. Because mono‐ubiquitination of BIK1 is required for its internalization, interference with this process by Avr2 could provide a mechanistic explanation for the compromised BIK1 mobility upon flg22 treatment. The identification of BIK1 as an effector target of a root‐invading vascular pathogen identifies this kinase as a conserved signalling component for both root and shoot immunity.

     
    more » « less
  2. 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
  3. Abstract

    Pathogen pressure on hosts can lead to the evolution of genes regulating the innate immune response. By characterizing naturally occurring polymorphisms in immune receptors, we can better understand the molecular determinants of pathogen recognition. ZAR1 is an ancientArabidopsis thalianaNLR (Nucleotide‐binding [NB] Leucine‐rich‐repeat [LRR] Receptor) that recognizes multiple secreted effector proteins from the pathogenic bacteriaPseudomonas syringaeandXanthomonas campestristhrough its interaction with receptor‐like cytoplasmic kinases (RLCKs). ZAR1 was first identified for its role in recognizingP. syringaeeffector HopZ1a, through its interaction with the RLCK ZED1. To identify additional determinants of HopZ1a recognition, we performed a computational screen for ecotypes from the 1001 Genomes project that were likely to lack HopZ1a recognition, and tested ~300 ecotypes. We identified ecotypes containing polymorphisms in ZAR1 and ZED1. Using our previously establishedNicotiana benthamianatransient assay and Arabidopsis ecotypes, we tested for the effect of naturally occurring polymorphisms on ZAR1 interactions and the immune response. We identified key residues in the NB or LRR domain of ZAR1 that impact the interaction with ZED1. We demonstrate that natural diversity combined with functional assays can help define the molecular determinants and interactions necessary to regulate immune induction in response to pathogens.

     
    more » « less
  4. Abstract

    In many animals, mate choice is important for the maintenance of reproductive isolation between species. Traits important for mate choice and behavioral isolation are predicted to be under strong stabilizing selection within species; however, such traits can also exhibit variation at the population level driven by neutral and adaptive evolutionary processes. Here, we describe patterns of divergence among androconial and genital chemical profiles at inter‐ and intraspecific levels in mimeticHeliconiusbutterflies. Most variation in chemical bouquets was found between species, but there were also quantitative differences at the population level. We found a strong correlation between interspecific chemical and genetic divergence, but this correlation varied in intraspecific comparisons. We identified “indicator” compounds characteristic of particular species that included compounds already known to elicit a behavioral response, suggesting an approach for identification of candidate compounds for future behavioral studies in novel systems. Overall, the strong signal of species identity suggests a role for these compounds in species recognition, but with additional potentially neutral variation at the population level.

     
    more » « less
  5. Abstract Background

    Fungal plant pathogens have dynamic genomes that allow them to rapidly adapt to adverse conditions and overcome host resistance. One way by which this dynamic genome plasticity is expressed is through effector gene loss, which enables plant pathogens to overcome recognition by cognate resistance genes in the host. However, the exact nature of these loses remains elusive in many fungi. This includes the tomato pathogenCladosporium fulvum, which is the first fungal plant pathogen from which avirulence (Avr) genes were ever cloned and in which loss ofAvrgenes is often reported as a means of overcoming recognition by cognate tomatoCfresistance genes. A recent near-complete reference genome assembly ofC. fulvumisolate Race 5 revealed a compartmentalized genome architecture and the presence of an accessory chromosome, thereby creating a basis for studying genome plasticity in fungal plant pathogens and its impact on avirulence genes.

    Results

    Here, we obtained near-complete genome assemblies of four additionalC. fulvumisolates. The genome assemblies had similar sizes (66.96 to 67.78 Mb), number of predicted genes (14,895 to 14,981), and estimated completeness (98.8 to 98.9%). Comparative analysis that included the genome of isolate Race 5 revealed high levels of synteny and colinearity, which extended to the density and distribution of repetitive elements and of repeat-induced point (RIP) mutations across homologous chromosomes. Nonetheless, structural variations, likely mediated by transposable elements and effecting the deletion of the avirulence genesAvr4E,Avr5, andAvr9, were also identified. The isolates further shared a core set of 13 chromosomes, but two accessory chromosomes were identified as well. Accessory chromosomes were significantly smaller in size, and one carried pseudogenized copies of two effector genes. Whole-genome alignments further revealed genomic islands of near-zero nucleotide diversity interspersed with islands of high nucleotide diversity that co-localized with repeat-rich regions. These regions were likely generated by RIP, which generally asymmetrically affected the genome ofC. fulvum.

    Conclusions

    Our results reveal new evolutionary aspects of theC. fulvumgenome and provide new insights on the importance of genomic structural variations in overcoming host resistance in fungal plant pathogens.

     
    more » « less