Characterization of microsatellite loci in neotropical Heliconius butterflies

The Heliconius butterflies offer exceptional opportunities for the study of the ecology and evolution of mimicry. Despite previous reports of difficulties in the development of microsatellite loci in Lepidoptera, we characterize 15 polymorphic loci in H. erato that show promise for genetic mapping and population studies in this and other species. Levels of variation were high, in both numbers and size ranges of alleles. The loci showed broad amplification success across the genus and in two other genera. All loci that amplified in a population of H. melpomene were polymorphic.     

A major locus controls a biologically active pheromone component in Heliconius melpomene

Understanding the production, response, and genetics of signals used in mate choice can inform our understanding of the evolution of both intraspecific mate choice and reproductive isolation. Sex pheromones are important for courtship and mate choice in many insects, but we know relatively little of their role in butterflies. The butterfly Heliconius melpomene uses a complex blend of wing androconial compounds during courtship. Electroantennography in H. melpomene and its close relative Heliconius cydno showed that responses to androconial extracts were not species specific. Females of both species responded equally strongly to extracts of both species, suggesting conservation of peripheral nervous system elements across the two species. Individual blend components provoked little to no response, with the exception of octadecanal, a major component of the H. melpomene blend. Supplementing octadecanal on the wings of octadecanal-rich H. melpomene males led to an increase in the time until mating, demonstrating the bioactivity of octadecanal in Heliconius. Using quantitative trait locus (QTL) mapping, we identified a single locus on chromosome 20 responsible for 41% of the parental species’ difference in octadecanal production. This QTL does not overlap with any of the major wing color or mate choice loci, nor does it overlap with known regions of elevated or reduced F_ST. A set of 16 candidate fatty acid biosynthesis genes lies underneath the QTL. Pheromones in Heliconius carry information relevant for mate choice and are under simple genetic control, suggesting they could be important during speciation.     

Beyond magic traits: Multimodal mating cues in Heliconius butterflies

Species coexistence involves the evolution of reproductive barriers opposing gene flow. Heliconius butterflies display colorful patterns affecting mate choice and survival through warning signaling and mimicry. These patterns are called “magic traits” for speciation because divergent natural selection may promote mimicry shifts in pattern whose role as mating cue facilitates reproductive isolation. By contrast, between comimetic species, natural selection promotes pattern convergence. We addressed whether visual convergence interferes with reproductive isolation by testing for sexual isolation between two closely related species with similar patterns, H. timareta thelxinoe and H. melpomene amaryllis. Experiments with models confirmed visual attraction based on wing phenotype, leading to indiscriminate approach. Nevertheless, mate choice experiments showed assortative mating. Monitoring male behavior toward live females revealed asymmetry in male preference, H. melpomene males courting both species equally while H. timareta males strongly preferred conspecifics. Experiments with hybrid males suggested an important genetic component for such asymmetry. Behavioral observations support a key role for short‐distance cues in determining male choice in H. timareta. Scents extracts from wings and genitalia revealed interspecific divergence in chemical signatures, and hybrid female scent composition was significantly associated with courtship intensity by H. timareta males, providing candidate chemical mating cues involved in sexual isolation.     

Morphological traits defining species differences in wild relatives of maize are controlled by multiple quantitative trait loci

We analyzed the genetic basis of morphological differences between two wild species of teosinte (Zea diploperennis and Z. mays ssp. parviglumis), which are relatives of maize. These two species differ in a number of taxonomically important traits including the structure of the tassel (male inflorescence), which is the focus of this report. To investigate the genetic inheritance of six tassel traits, quantitative trait locus (QTL) mapping with 95 RFLP markers was employed on a population of 425 F2 plants. Each trait was analyzed by interval mapping (IM) and composite interval mapping (CIM) to identify and characterize the QTL controlling the differences in tassel morphology. We detected two to eight QTL for each trait. In total, 30 QTL with IM and 33 QTL with CIM were found for tassel morphology. QTL for several of the traits mapped near each other, suggesting pleiotropy and/or linkage of QTL. The QTL showed small to moderate magnitudes of effect. No QTL of exceptionally large effect were found as seen under domestication and in the case of some other natural species. Thus, the model involving major QTL of large effect seems not to apply to the traits and species analyzed. A mixture of QTL with positive and negative allelic effects was found for most tassel traits and may suggest a history of periodic changes in the direction of selection during the divergence of Z. diploperennis and Z. mays ssp. parviglumis or fixation of QTL alleles by random genetic drift rather than selection.     

Sequestration and biosynthesis of cyanogenic glucosides in passion vine butterflies and consequences for the diversification of their host plants

The colorful heliconiine butterflies are distasteful to predators due to their content of defense compounds called cyanogenic glucosides (CNglcs), which they biosynthesize from aliphatic amino acids. Heliconiine larvae feed exclusively on Passiflora plants where ~30 kinds of CNglcs have been reported. Among them, some CNglcs derived from cyclopentenyl glycine can be sequestered by some Heliconius species. In order to understand the evolution of biosynthesis and sequestration of CNglcs in these butterflies and its consequences for their arms race with Passiflora plants, we analyzed the CNglc distribution in selected heliconiine and Passiflora species. Sequestration of cyclopentenyl CNglcs is not an exclusive trait of Heliconius, since these compounds were present in other heliconiines such as Philaethria, Dryas and Agraulis, and in more distantly related genera Cethosia and Euptoieta. Thus, it is likely that the ability to sequester cyclopentenyl CNglcs arose in an ancestor of the Heliconiinae subfamily. Biosynthesis of aliphatic CNglcs is widespread in these butterflies, although some species from the sara‐sapho group seem to have lost this ability. The CNglc distribution within Passiflora suggests that they might have diversified their cyanogenic profile to escape heliconiine herbivory. This systematic analysis improves our understanding on the evolution of cyanogenesis in the heliconiine–Passiflora system.     

Ecological and genetic factors influencing the transition between host‐use strategies in sympatric Heliconius butterflies

Shifts in host‐plant use by phytophagous insects have played a central role in their diversification. Evolving host‐use strategies will reflect a trade‐off between selection pressures. The ecological niche of herbivorous insects is partitioned along several dimensions, and if populations remain in contact, recombination will break down associations between relevant loci. As such, genetic architecture can profoundly affect the coordinated divergence of traits and subsequently the ability to exploit novel habitats. The closely related species Heliconius cydno and H. melpomene differ in mimetic colour pattern, habitat and host‐plant use. We investigate the selection pressures and genetic basis underlying host‐use differences in these two species. Host‐plant surveys reveal that H. melpomene specializes on a single species of Passiflora. This is also true for the majority of other Heliconius species in secondary growth forest at our study site, as expected under a model of interspecific competition. In contrast, H. cydno, which uses closed‐forest habitats where both Heliconius and Passiflora are less common, appears not to be restricted by competition and uses a broad selection of the available Passiflora. However, other selection pressures are likely involved, and field experiments reveal that early larval survival of both butterfly species is highest on Passiflora menispermifolia, but most markedly so for H. melpomene, the specialist on that host. Finally, we demonstrate an association between host‐plant acceptance and colour pattern amongst interspecific hybrids, suggesting that major loci underlying these important ecological traits are physically linked in the genome. Together, our results reveal ecological and genetic associations between shifts in habitat, host use and mimetic colour pattern that have likely facilitated both speciation and coexistence.     

Species specificity and intraspecific variation in the chemical profiles of Heliconius butterflies across a large geographic range

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 mimetic Heliconius butterflies. 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.     

Pollen feeding,resource allocation and the evolution of chemical defence in passion vine butterflies

Evolution of pollen feeding in Heliconius has allowed exploitation of rich amino acid sources and dramatically reorganized life‐history traits. In Heliconius, eggs are produced mainly from adult‐acquired resources, leaving somatic development and maintenance to larva effort. This innovation may also have spurred evolution of chemical defence via amino acid‐derived cyanogenic glycosides. In contrast, nonpollen‐feeding heliconiines must rely almost exclusively on larval‐acquired resources for both reproduction and defence. We tested whether adult amino acid intake has an immediate influence on cyanogenesis in Heliconius. Because Heliconius are more distasteful to bird predators than close relatives that do not utilize pollen, we also compared cyanogenesis due to larval input across Heliconius species and nonpollen‐feeding relatives. Except for one species, we found that varying the amino acid diet of an adult Heliconius has negligible effect on its cyanide concentration. Adults denied amino acids showed no decrease in cyanide and no adults showed cyanide increase when fed amino acids. Yet, pollen‐feeding butterflies were capable of producing more defence than nonpollen‐feeding relatives and differences were detectable in freshly emerged adults, before input of adult resources. Our data points to a larger role of larval input in adult chemical defence. This coupled with the compartmentalization of adult nutrition to reproduction and longevity suggests that one evolutionary consequence of pollen feeding, shifting the burden of reproduction to adults, is to allow the evolution of greater allocation of host plant amino acids to defensive compounds by larvae.     

Sharp genetic discontinuity across a unimodal Heliconius hybrid zone

Hybrid zones are powerful natural systems to study evolutionary processes to gain an understanding of adaptation and speciation. In the Cauca Valley (Colombia), two butterfly races, Heliconius cydno cydnides and Heliconius cydno weymeri, meet and hybridize. We characterized this hybrid zone using a combination of mitochondrial DNA (mtDNA) sequences, amplified fragment length polymorphisms (AFLPs), microsatellites and sequences for nuclear loci within and outside of the genomic regions that cause differences in wing colour pattern. The hybrid zone is largely composed of individuals of mixed ancestry. However, there is strong genetic discontinuity between the hybridizing races in mtDNA and, to a lesser extent, in all nuclear markers surveyed. The mtDNA clustering of H. c. cydnides with the H. cydno race from the Magdalena Valley and H. c. weymeri with the H. cydno race from the pacific coast suggests that H. c. cydnides colonized the Cauca Valley from the north, whereas H. c. weymeri did so by crossing the Andes in the southern part, implying a secondary contact origin. Colonization of the valley by H. cydno was accompanied by mimicry shift. Strong ecological isolation, driven by locally adaptive differences in mimetic wing patterns, is playing an important role in maintaining the hybrid zone. However, selection on wing pattern alone is not sufficient to explain the genetic discontinuity observed. There is evidence for differences in male mating preference, but the contribution of additional barriers needs further investigation. Overall, our results support the idea that speciation is a cumulative process, where the combination of multiple isolation barriers, combined with major phenotypic differences, facilitates population divergence in face of gene flow.     

Movement of a Heliconius hybrid zone over 30 years: A Bayesian approach

Hybrid zones have long been of interest to biologists as natural laboratories where we can gain insight into the processes of adaptation and speciation. Repeated sampling of individual hybrid zones has been particularly useful in elucidating the dynamic balance between selection and dispersal that maintains most hybrid zones. Here, we revisit a hybrid zone between Heliconius erato butterflies in Panamá for a third time over more than 30 years. We combine a novel Bayesian extension of stepped‐cline hybrid zone models with environmental data to understand the genetic and environmental causes of cline dynamics in this species. The cline has continued to move west, likely due to dominance drive, but has slowed and broadened. Environmental analyses suggest that widespread deforestation in Panamá could be leading to decreased avian predation and relaxed selection, causing the observed changes in cline dynamics.     

Evolution of novel mimicry rings facilitated by adaptive introgression in tropical butterflies

Understanding the genetic basis of phenotypic variation and the mechanisms involved in the evolution of adaptive novelty, especially in adaptive radiations, is a major goal in evolutionary biology. Here, we used whole‐genome sequence data to investigate the origin of the yellow hindwing bar in the Heliconius cydno radiation. We found modular variation associated with hindwing phenotype in two narrow noncoding regions upstream and downstream of the cortex gene, which was recently identified as a pigmentation pattern controller in multiple species of Heliconius. Genetic variation at each of these modules suggests an independent control of the dorsal and ventral hindwing patterning, with the upstream module associated with the ventral phenotype and the downstream module with the dorsal one. Furthermore, we detected introgression between H. cydno and its closely related species Heliconius melpomene in these modules, likely allowing both species to participate in novel mimicry rings. In sum, our findings support the role of regulatory modularity coupled with adaptive introgression as an elegant mechanism by which novel phenotypic combinations can evolve and fuel an adaptive radiation.     

Linkage disequilibrium fine mapping of quantitative trait loci: A simulation study

Recently, the use of linkage disequilibrium (LD) to locate genes which affect quantitative traits (QTL) has received an increasing interest, but the plausibility of fine mapping using linkage disequilibrium techniques for QTL has not been well studied. The main objectives of this work were to (1) measure the extent and pattern of LD between a putative QTL and nearby markers in finite populations and (2) investigate the usefulness of LD in fine mapping QTL in simulated populations using a dense map of multiallelic or biallelic marker loci. The test of association between a marker and QTL and the power of the test were calculated based on single-marker regression analysis. The results show the presence of substantial linkage disequilibrium with closely linked marker loci after 100 to 200 generations of random mating. Although the power to test the association with a frequent QTL of large effect was satisfactory, the power was low for the QTL with a small effect and/or low frequency. More powerful, multi-locus methods may be required to map low frequent QTL with small genetic effects, as well as combining both linkage and linkage disequilibrium information. The results also showed that multiallelic markers are more useful than biallelic markers to detect linkage disequilibrium and association at an equal distance.     

Variation of chemical compounds in wild Heliconiini reveals ecological factors involved in the evolution of chemical defenses in mimetic butterflies

Evolutionary convergence of color pattern in mimetic species is tightly linked with the evolution of chemical defenses. Yet, the evolutionary forces involved in natural variations of chemical defenses in aposematic species are still understudied. Herein, we focus on the evolution of chemical defenses in the butterfly tribe Heliconiini. These neotropical butterflies contain large concentrations of cyanogenic glucosides, cyanide‐releasing compounds acting as predator deterrent. These compounds are either de novo synthesized or sequestered from their Passiflora host plant, so that their concentrations may depend on host plant specialization and host plant availability. We sampled 375 wild Heliconiini butterflies across Central and South America, covering 43% species of this clade, and quantify individual variations in the different CGs using liquid chromatography coupled with tandem mass spectrometry. We detected new compounds and important variations in chemical defenses both within and among species. Based on the most recent and well‐studied phylogeny of Heliconiini, we show that ecological factors such as mimetic interactions and host plant specialization have a significant association with chemical profiles, but these effects are largely explained by phylogenetic relationships. Our results therefore suggest that shared ancestries largely contribute to chemical defense variation, pointing out at the interaction between historical and ecological factors in the evolution of Müllerian mimicry.     

Detection of quantitative trait loci for carcass composition traits in pigs

A quantitative trait locus (QTL) analysis of carcass composition data from a three-generation experimental cross between Meishan (MS) and Large White (LW) pig breeds is presented. A total of 488 F2 males issued from six F1 boars and 23 F1 sows, the progeny of six LW boars and six MS sows, were slaughtered at approximately 80 kg live weight and were submitted to a standardised cutting of the carcass. Fifteen traits, i.e. dressing percentage, loin, ham, shoulder, belly, backfat, leaf fat, feet and head weights, two backfat thickness and one muscle depth measurements, ham + loin and back + leaf fat percentages and estimated carcass lean content were analysed. Animals were typed for a total of 137 markers covering the entire porcine genome. Analyses were performed using a line-cross (LC) regression method where founder lines were assumed to be fixed for different QTL alleles and a half/full sib (HFS) maximum likelihood method where allele substitution effects were estimated within each half-/full-sib family. Additional analyses were performed to search for multiple linked QTL and imprinting effects. Significant gene effects were evidenced for both leanness and fatness traits in the telomeric regions of SSC 1q and SSC 2p, on SSC 4, SSC 7 and SSC X. Additional significant QTL were identified for ham weight on SSC 5, for head weight on SSC 1 and SSC 7, for feet weight on SSC 7 and for dressing percentage on SSC X. LW alleles were associated with a higher lean content and a lower fat content of the carcass, except for the fatness trait on SSC 7. Suggestive evidence of linked QTL on SSC 7 and of imprinting effects on SSC 6, SSC 7, SSC 9 and SSC 17 were also obtained.     

Adenylate cyclase 5 is required for melanophore and male pattern development in the guppy (Poecilia reticulata)

Guppies (Poecilia reticulata) are colorful fish that have attracted the attention of pigmentation researchers for almost a century. Here, we report that the blond phenotype of the guppy is caused by a spontaneous mutation in the guppy ortholog of adenylate cyclase 5 (adcy5). Using double digest restriction site‐associated DNA sequencing (ddRADseq) and quantitative trait locus (QTL) mapping, we linked the blond phenotype to a candidate region of 118 kb, in which we subsequently identified a 2‐bp deletion in adcy5 that alters splicing and leads to a premature stop codon. We show that adcy5, which affects life span and melanoma growth in mouse, is required for melanophore development and formation of male orange pigmentation traits in the guppy. We find that some components of the male orange pattern are particularly sensitive to loss of Adcy5 function. Our work thus reveals a function for Adcy5 in patterning of fish color ornaments.     

Mitogenomics reveals high synteny and long evolutionary histories of sympatric cryptic nematode species

Species with seemingly identical morphology but with distinct genetic differences are abundant in the marine environment and frequently co‐occur in the same habitat. Such cryptic species are typically delineated using a limited number of mitochondrial and/or nuclear marker genes, which do not yield information on gene order and gene content of the genomes under consideration. We used next‐generation sequencing to study the composition of the mitochondrial genomes of four sympatrically distributed cryptic species of the Litoditis marina species complex (PmI, PmII, PmIII, and PmIV). The ecology, biology, and natural occurrence of these four species are well known, but the evolutionary processes behind this cryptic speciation remain largely unknown. The gene order of the mitochondrial genomes of the four species was conserved, but differences in genome length, gene length, and codon usage were observed. The atp8 gene was lacking in all four species. Phylogenetic analyses confirm that PmI and PmIV are sister species and that PmIII diverged earliest. The most recent common ancestor of the four cryptic species was estimated to have diverged 16 MYA. Synonymous mutations outnumbered nonsynonymous changes in all protein‐encoding genes, with the Complex IV genes (coxI‐III) experiencing the strongest purifying selection. Our mitogenomic results show that morphologically similar species can have long evolutionary histories and that PmIII has several differences in genetic makeup compared to the three other species, which may explain why it is better adapted to higher temperatures than the other species.     

Facultative pupal mating in Heliconius erato: Implications for mate choice,female preference,and speciation

Mating systems have broad impacts on how sexual selection and mate choice operate within a species, but studies of mating behavior in the laboratory may not reflect how these processes occur in the wild. Here, we examined the mating behavior of the neotropical butterfly Heliconius erato in the field by releasing larvae and virgin females and observing how they mated. H. erato is considered a pupal‐mating species (i.e., males mate with females as they emerge from the pupal case). However, we observed only two teneral mating events, and experimentally released virgins were almost all mated upon recapture. Our study confirms the presence of some pupal‐mating behavior in H. erato, but suggests that adult mating is likely the prevalent mating strategy in this species. These findings have important implications for the role of color pattern and female mate choice in the generation of reproductive isolation in this diverse genus.     

Detection of quantitative trait loci for growth and fatness in pigs

A quantitative trait locus (QTL) analysis of growth and fatness data from a three-generation experimental cross between Meishan (MS) and Large White (LW) pig breeds is presented. Six boars and 23 F1 sows, the progeny of six LW boars and six MS sows, produced 530 F2 males and 573 F2 females. Nine growth traits, i.e. body weight at birth and at 3, 10, 13, 17 and 22 weeks of age, average daily gain from birth to 3 weeks, from 3 to 10 weeks and from 10 to 22 weeks of age, as well as backfat thickness at 13, 17 and 22 weeks of age and at 40 and 60 kg live weight were analysed. Animals were typed for a total of 137 markers covering the entire porcine genome. Analyses were performed using two interval mapping methods: a line-cross (LC) regression method where founder lines were assumed to be fixed for different QTL alleles and a half-/full-sib (HFS) maximum likelihood method where allele substitution effects were estimated within each half-/full-sib family. Both methods revealed highly significant gene effects for growth on chromosomes 1, 4 and 7 and for backfat thickness on chromosomes 1, 4, 5, 7 and X, and significant gene effects on chromosome 6 for growth and backfat thickness. Suggestive QTLs were also revealed by both methods on chromosomes 2 and 3 for growth and 2 for backfat thickness. Significant gene effects were detected for growth on chromosomes 11, 13, 14, 16 and 18 and for backfat thickness on chromosome 8, 10, 13 and 14. LW alleles were associated with high growth rate and low backfat thickness, except for those of chromosome 7 and to a lesser extent early-growth alleles on chromosomes 1 and 2 and backfat thickness alleles on chromosome 6.     

Evidence for nonallopatric speciation among closely related sympatric Heliotropium species in the Atacama Desert

The genetic structure of populations of closely related, sympatric species may hold the signature of the geographical mode of the speciation process. In fully allopatric speciation, it is expected that genetic differentiation between species is homogeneously distributed across the genome. In nonallopatric speciation, the genomes may remain undifferentiated to a large extent. In this article, we analyzed the genetic structure of five sympatric species from the plant genus Heliotropium in the Atacama Desert. We used amplified fragment length polymorphisms (AFLPs) to characterize the genetic structure of these species and evaluate their genetic differentiation as well as the number of loci subject to positive selection using divergence outlier analysis (DOA). The five species form distinguishable groups in the genetic space, with zones of overlap, indicating that they are possibly not completely isolated. Among‐species differentiation accounts for 35% of the total genetic differentiation (F_ST = 0.35), and F_ST between species pairs is positively correlated with phylogenetic distance. DOA suggests that few loci are subject to positive selection, which is in line with a scenario of nonallopatric speciation. These results support the idea that sympatric species of Heliotropium sect. Cochranea are under an ongoing speciation process, characterized by a fluctuation of population ranges in response to pulses of arid and humid periods during Quaternary times.     

Model selection for quantitative trait loci mapping in a full-sib family

Statistical methods for mapping quantitative trait loci (QTLs) in full-sib forest trees, in which the number of alleles and linkage phase can vary from locus to locus, are still not well established. Previous studies assumed that the QTL segregation pattern was fixed throughout the genome in a full-sib family, despite the fact that this pattern can vary among regions of the genome. In this paper, we propose a method for selecting the appropriate model for QTL mapping based on the segregation of different types of markers and QTLs in a full-sib family. The QTL segregation patterns were classified into three types: test cross (1:1 segregation), F₂ cross (1:2:1 segregation) and full cross (1:1:1:1 segregation). Akaike’s information criterion (AIC), the Bayesian information criterion (BIC) and the Laplace-empirical criterion (LEC) were used to select the most likely QTL segregation pattern. Simulations were used to evaluate the power of these criteria and the precision of parameter estimates. A Windows-based software was developed to run the selected QTL mapping method. A real example is presented to illustrate QTL mapping in forest trees based on an integrated linkage map with various segregation markers. The implications of this method for accurate QTL mapping in outbred species are discussed.