The evolution and genetics of sexually dimorphic 'dual' mimicry in the butterfly Elymnias hypermnestra.

Sexual dimorphism is a major component of morphological variation across the tree of life, but the mechanisms underlying phenotypic differences between sexes of a single species are poorly understood. We examined the population genomics and biogeography of the common palmfly Elymnias hypermnestra, a dual mimic in which female wing colour patterns are either dark brown (melanic) or bright orange, mimicking toxic Euploea and Danaus species, respectively. As males always have a melanic wing colour pattern, this makes E. hypermnestra a fascinating model organism in which populations vary in sexual dimorphism. Population structure analysis revealed that there were three genetically distinct E. hypermnestra populations, which we further validated by creating a phylogenomic species tree and inferring historical barriers to gene flow. This species tree demonstrated that multiple lineages with orange females do not form a monophyletic group, and the same is true of clades with melanic females. We identified two single nucleotide polymorphisms (SNPs) near the colour patterning gene WntA that were significantly associated with the female colour pattern polymorphism, suggesting that this gene affects sexual dimorphism. Given WntA's role in colour patterning across Nymphalidae, E. hypermnestra females demonstrate the repeatability of the evolution of sexual dimorphism.

The roles of hybridization and habitat fragmentation in the evolution of Brazil's enigmatic longwing butterflies, Heliconius nattereri and H. hermathena.

BACKGROUND: Heliconius butterflies are widely distributed across the Neotropics and have evolved a stunning array of wing color patterns that mediate Müllerian mimicry and mating behavior. Their rapid radiation has been strongly influenced by hybridization, which has created new species and allowed sharing of color patterning alleles between mimetic species pairs. While these processes have frequently been observed in widespread species with contiguous distributions, many Heliconius species inhabit patchy or rare habitats that may strongly influence the origin and spread of species and color patterns. Here, we assess the effects of historical population fragmentation and unique biology on the origins, genetic health, and color pattern evolution of two rare and sparsely distributed Brazilian butterflies, Heliconius hermathena and Heliconius nattereri. RESULTS: We assembled genomes and re-sequenced whole genomes of eight H. nattereri and 71 H. hermathena individuals. These species harbor little genetic diversity, skewed site frequency spectra, and high deleterious mutation loads consistent with recent population bottlenecks. Heliconius hermathena consists of discrete, strongly isolated populations that likely arose from a single population that dispersed after the last glacial maximum. Despite having a unique color pattern combination that suggested a hybrid origin, we found no genome-wide evidence that H. hermathena is a hybrid species. However, H. hermathena mimicry evolved via introgression, from co-mimetic Heliconius erato, of a small genomic region upstream of the color patterning gene cortex. CONCLUSIONS: Heliconius hermathena and H. nattereri population fragmentation, potentially driven by historical climate change and recent deforestation, has significantly reduced the genetic health of these rare species. Our results contribute to a growing body of evidence that introgression of color patterning alleles between co-mimetic species appears to be a general feature of Heliconius evolution.

Genomic evidence for gene flow between monarchs with divergent migratory phenotypes and flight performance.

Monarch butterflies are known for their spectacular annual migration in eastern North America, with millions of monarchs flying up to 4,500 km to overwintering sites in central Mexico. Monarchs also live west of the Rocky Mountains, where they travel shorter distances to overwinter along the Pacific Coast. It is often assumed that eastern and western monarchs form distinct evolutionary units, but genomic studies to support this notion are lacking. We used a tethered flight mill to show that migratory eastern monarchs have greater flight performance than western monarchs, consistent with their greater migratory distances. However, analysing more than 20 million SNPs in 43 monarch genomes, we found no evidence for genomic differentiation between eastern and western monarchs. Genomic analysis also showed identical and low levels of genetic diversity, and demographic analyses indicated similar effective population sizes and ongoing gene flow between eastern and western monarchs. Gene expression analysis of a subset of candidate genes during active flight revealed differential gene expression related to nonmuscular motor activity. Our results demonstrate that eastern and western monarchs maintain migratory differences despite ongoing gene flow, and suggest that migratory differences between eastern and western monarchs are not driven by select major-effects alleles. Instead, variation in migratory distance and destination may be driven by environmentally induced differential gene expression or by many alleles of small effect.

Disentangling Population History and Character Evolution among Hybridizing Lineages.

Understanding the origin and maintenance of adaptive phenotypic novelty is a central goal of evolutionary biology. However, both hybridization and incomplete lineage sorting can lead to genealogical discordance between the regions of the genome underlying adaptive traits and the remainder of the genome, decoupling inferences about character evolution from population history. Here, to disentangle these effects, we investigated the evolutionary origins and maintenance of Batesian mimicry between North American admiral butterflies (Limenitis arthemis) and their chemically defended model (Battus philenor) using a combination of de novo genome sequencing, whole-genome resequencing, and statistical introgression mapping. Our results suggest that balancing selection, arising from geographic variation in the presence or absence of the unpalatable model, has maintained two deeply divergent color patterning haplotypes that have been repeatedly sieved among distinct mimetic and nonmimetic lineages of Limenitis via introgressive hybridization.

Butterfly Mimicry Polymorphisms Highlight Phylogenetic Limits of Gene Reuse in the Evolution of Diverse Adaptations.

Some genes have repeatedly been found to control diverse adaptations in a wide variety of organisms. Such gene reuse reveals not only the diversity of phenotypes these unique genes control but also the composition of developmental gene networks and the genetic routes available to and taken by organisms during adaptation. However, the causes of gene reuse remain unclear. A small number of large-effect Mendelian loci control a huge diversity of mimetic butterfly wing color patterns, but reasons for their reuse are difficult to identify because the genetic basis of mimicry has primarily been studied in two systems with correlated factors: female-limited Batesian mimicry in Papilio swallowtails (Papilionidae) and non-sex-limited Müllerian mimicry in Heliconius longwings (Nymphalidae). Here, we break the correlation between phylogenetic relationship and sex-limited mimicry by identifying loci controlling female-limited mimicry polymorphism Hypolimnas misippus (Nymphalidae) and non-sex-limited mimicry polymorphism in Papilio clytia (Papilionidae). The Papilio clytia polymorphism is controlled by the genome region containing the gene cortex, the classic P supergene in Heliconius numata, and loci controlling color pattern variation across Lepidoptera. In contrast, female-limited mimicry polymorphism in Hypolimnas misippus is associated with a locus not previously implicated in color patterning. Thus, although many species repeatedly converged on cortex and its neighboring genes over 120 My of evolution of diverse color patterns, female-limited mimicry polymorphisms each evolved using a different gene. Our results support conclusions that gene reuse occurs mainly within ∼10 My and highlight the puzzling diversity of genes controlling seemingly complex female-limited mimicry polymorphisms.

The molecular genetic basis of herbivory between butterflies and their host plants.

Interactions between herbivorous insects and their host plants are a central component of terrestrial food webs and a critical topic in agriculture, where a substantial fraction of potential crop yield is lost annually to pests. Important insights into plant-insect interactions have come from research on specific plant defences and insect detoxification mechanisms. Yet, much remains unknown about the molecular mechanisms that mediate plant-insect interactions. Here we use multiple genome-wide approaches to map the molecular basis of herbivory from both plant and insect perspectives, focusing on butterflies and their larval host plants. Parallel genome-wide association studies in the cabbage white butterfly, Pieris rapae, and its host plant, Arabidopsis thaliana, pinpointed a small number of butterfly and plant genes that influenced herbivory. These genes, along with much of the genome, were regulated in a dynamic way over the time course of the feeding interaction. Comparative analyses, including diverse butterfly/plant systems, showed a variety of genome-wide responses to herbivory, as well as a core set of highly conserved genes in butterflies as well as their host plants. These results greatly expand our understanding of the genomic causes and evolutionary consequences of ecological interactions across two of nature's most diverse taxa, butterflies and flowering plants.

The long-term maintenance of a resistance polymorphism through diffuse interactions.

Plant resistance (R) genes are a crucial component in plant defence against pathogens. Although R genes often fail to provide durable resistance in an agricultural context, they frequently persist as long-lived balanced polymorphisms in nature. Standard theory explains the maintenance of such polymorphisms through a balance of the costs and benefits of resistance and virulence in a tightly coevolving host-pathogen pair. However, many plant-pathogen interactions lack such specificity. Whether, and how, balanced polymorphisms are maintained in diffusely interacting species is unknown. Here we identify a naturally interacting R gene and effector pair in Arabidopsis thaliana and its facultative plant pathogen, Pseudomonas syringae. The protein encoded by the R gene RPS5 recognizes an AvrPphB homologue (AvrPphB2) and exhibits a balanced polymorphism that has been maintained for over 2 million years (ref. 3). Consistent with the presence of an ancient balanced polymorphism, the R gene confers a benefit when plants are infected with P. syringae carrying avrPphB2 but also incurs a large cost in the absence of infection. RPS5 alleles are maintained at intermediate frequencies in populations globally, suggesting ubiquitous selection for resistance. However, the presence of P. syringae carrying avrPphB is probably insufficient to explain the RPS5 polymorphism. First, avrPphB homologues occur at very low frequencies in P. syringae populations on A. thaliana. Second, AvrPphB only rarely confers a virulence benefit to P. syringae on A. thaliana. Instead, we find evidence that selection for RPS5 involves multiple non-homologous effectors and multiple pathogen species. These results and an associated model suggest that the R gene polymorphism in A. thaliana may not be maintained through a tightly coupled interaction involving a single coevolved R gene and effector pair. More likely, the stable polymorphism is maintained through complex and diffuse community-wide interactions.

Patterns of divergence of a large family of nodule cysteine-rich peptides in accessions of Medicago truncatula.

The nodule cysteine-rich (NCR) groups of defensin-like (DEFL) genes are one of the largest gene families expressed in the nodules of some legume plants. They have only been observed in the inverted repeat loss clade (IRLC) of legumes, which includes the model legume Medicago truncatula. NCRs are reported to play an important role in plant-microbe interactions. To understand their diversity we analyzed their expression and sequence polymorphisms among four accessions of M. truncatula. A significant expression and nucleotide variation was observed among the genes. We then used 26 accessions to estimate the selection pressures shaping evolution among the accessions by calculating the nucleotide diversity at non-synonymous and synonymous sites in the coding region. The mature peptides of the orthologous NCRs had signatures of both purifying and diversifying selection pressures, unlike the seed DEFLs, which predominantly exhibited purifying selection. The expression, sequence variation and apparent diversifying selection in NCRs within the Medicago species indicates rapid and recent evolution, and suggests that this family of genes is actively evolving to adapt to different environments and is acquiring new functions.

Detecting small plant peptides using SPADA (Small Peptide Alignment Discovery Application).

BACKGROUND: Small peptides encoded as one- or two-exon genes in plants have recently been shown to affect multiple aspects of plant development, reproduction and defense responses. However, popular similarity search tools and gene prediction techniques generally fail to identify most members belonging to this class of genes. This is largely due to the high sequence divergence among family members and the limited availability of experimentally verified small peptides to use as training sets for homology search and ab initio prediction. Consequently, there is an urgent need for both experimental and computational studies in order to further advance the accurate prediction of small peptides. RESULTS: We present here a homology-based gene prediction program to accurately predict small peptides at the genome level. Given a high-quality profile alignment, SPADA identifies and annotates nearly all family members in tested genomes with better performance than all general-purpose gene prediction programs surveyed. We find numerous mis-annotations in the current Arabidopsis thaliana and Medicago truncatula genome databases using SPADA, most of which have RNA-Seq expression support. We also show that SPADA works well on other classes of small secreted peptides in plants (e.g., self-incompatibility protein homologues) as well as non-secreted peptides outside the plant kingdom (e.g., the alpha-amanitin toxin gene family in the mushroom, Amanita bisporigera). CONCLUSIONS: SPADA is a free software tool that accurately identifies and predicts the gene structure for short peptides with one or two exons. SPADA is able to incorporate information from profile alignments into the model prediction process and makes use of it to score different candidate models. SPADA achieves high sensitivity and specificity in predicting small plant peptides such as the cysteine-rich peptide families. A systematic application of SPADA to other classes of small peptides by research communities will greatly improve the genome annotation of different protein families in public genome databases.

Regulatory patterns of a large family of defensin-like genes expressed in nodules of Medicago truncatula.

Root nodules are the symbiotic organ of legumes that house nitrogen-fixing bacteria. Many genes are specifically induced in nodules during the interactions between the host plant and symbiotic rhizobia. Information regarding the regulation of expression for most of these genes is lacking. One of the largest gene families expressed in the nodules of the model legume Medicago truncatula is the nodule cysteine-rich (NCR) group of defensin-like (DEFL) genes. We used a custom Affymetrix microarray to catalog the expression changes of 566 NCRs at different stages of nodule development. Additionally, bacterial mutants were used to understand the importance of the rhizobial partners in induction of NCRs. Expression of early NCRs was detected during the initial infection of rhizobia in nodules and expression continued as nodules became mature. Late NCRs were induced concomitantly with bacteroid development in the nodules. The induction of early and late NCRs was correlated with the number and morphology of rhizobia in the nodule. Conserved 41 to 50 bp motifs identified in the upstream 1,000 bp promoter regions of NCRs were required for promoter activity. These cis-element motifs were found to be unique to the NCR family among all annotated genes in the M. truncatula genome, although they contain sub-regions with clear similarity to known regulatory motifs involved in nodule-specific expression and temporal gene regulation.

Spatio-temporal expression patterns of Arabidopsis thaliana and Medicago truncatula defensin-like genes.

Plant genomes contain several hundred defensin-like (DEFL) genes that encode short cysteine-rich proteins resembling defensins, which are well known antimicrobial polypeptides. Little is known about the expression patterns or functions of many DEFLs because most were discovered recently and hence are not well represented on standard microarrays. We designed a custom Affymetrix chip consisting of probe sets for 317 and 684 DEFLs from Arabidopsis thaliana and Medicago truncatula, respectively for cataloging DEFL expression in a variety of plant organs at different developmental stages and during symbiotic and pathogenic associations. The microarray analysis provided evidence for the transcription of 71% and 90% of the DEFLs identified in Arabidopsis and Medicago, respectively, including many of the recently annotated DEFL genes that previously lacked expression information. Both model plants contain a subset of DEFLs specifically expressed in seeds or fruits. A few DEFLs, including some plant defensins, were significantly up-regulated in Arabidopsis leaves inoculated with Alternaria brassicicola or Pseudomonas syringae pathogens. Among these, some were dependent on jasmonic acid signaling or were associated with specific types of immune responses. There were notable differences in DEFL gene expression patterns between Arabidopsis and Medicago, as the majority of Arabidopsis DEFLs were expressed in inflorescences, while only a few exhibited root-enhanced expression. By contrast, Medicago DEFLs were most prominently expressed in nitrogen-fixing root nodules. Thus, our data document salient differences in DEFL temporal and spatial expression between Arabidopsis and Medicago, suggesting distinct signaling routes and distinct roles for these proteins in the two plant species.