A genetic switch for male UV iridescence in an incipient species pair of sulphur butterflies.

Mating cues evolve rapidly and can contribute to species formation and maintenance. However, little is known about how sexual signals diverge and how this variation integrates with other barrier loci to shape the genomic landscape of reproductive isolation. Here, we elucidate the genetic basis of ultraviolet (UV) iridescence, a courtship signal that differentiates the males of Colias eurytheme butterflies from a sister species, allowing females to avoid costly heterospecific matings. Anthropogenic range expansion of the two incipient species established a large zone of secondary contact across the eastern United States with strong signatures of genomic admixtures spanning all autosomes. In contrast, Z chromosomes are highly differentiated between the two species, supporting a disproportionate role of sex chromosomes in speciation known as the large-X (or large-Z) effect. Within this chromosome-wide reproductive barrier, linkage mapping indicates that cis-regulatory variation of bric a brac (bab) underlies the male UV-iridescence polymorphism between the two species. Bab is expressed in all non-UV scales, and butterflies of either species or sex acquire widespread ectopic iridescence following its CRISPR knockout, demonstrating that Bab functions as a suppressor of UV-scale differentiation that potentiates mating cue divergence. These results highlight how a genetic switch can regulate a premating signal and integrate with other reproductive barriers during intermediate phases of speciation.

Evidence for a single, ancient origin of a genus-wide alternative life history strategy.

Understanding the evolutionary origins and factors maintaining alternative life history strategies (ALHS) within species is a major goal of evolutionary research. While alternative alleles causing discrete ALHS are expected to purge or fix over time, one-third of the ~90 species of Colias butterflies are polymorphic for a female-limited ALHS called Alba. Whether Alba arose once, evolved in parallel, or has been exchanged among taxa is currently unknown. Using comparative genome-wide association study (GWAS) and population genomic analyses, we placed the genetic basis of Alba in time-calibrated phylogenomic framework, revealing that Alba evolved once near the base of the genus and has been subsequently maintained via introgression and balancing selection. CRISPR-Cas9 mutagenesis was then used to verify a putative cis-regulatory region of Alba, which we identified using phylogenetic foot printing. We hypothesize that this cis-regulatory region acts as a modular enhancer for the induction of the Alba ALHS, which has likely facilitated its long evolutionary persistence.

Genetics of yellow-orange color variation in a pair of sympatric sulphur butterflies.

Continuous color polymorphisms can serve as a tractable model for the genetic and developmental architecture of traits. Here we investigated continuous color variation in Colias eurytheme and Colias philodice, two species of sulphur butterflies that hybridize in sympatry. Using quantitative trait locus (QTL) analysis and high-throughput color quantification, we found two interacting large-effect loci affecting orange-to-yellow chromaticity. Knockouts of red Malpighian tubules (red), likely involved in endosomal maturation, result in depigmented wing scales. Additionally, the transcription factor bric-a-brac can act as a modulator of orange pigmentation. We also describe the QTL architecture of other continuously varying traits, together supporting a large-X effect model where the genetic control of species-defining traits is enriched on sex chromosomes. This study sheds light on the range of possible genetic architectures that can underpin a continuously varying trait and illustrates the power of using automated measurement to score phenotypes that are not always conspicuous to the human eye.

AFLP linkage map of hybridizing swallowtail butterflies, Papilio glaucus and Papilio canadensis.

High-density linkage maps provide powerful tools for studying the genetic basis of ecologically relevant adaptations and the genomic scope of introgression. We backcrossed an F(1) hybrid male Papilio glaucus/Papilio canadensis tiger swallowtail butterfly to a pure P. glaucus female and constructed amplified fragment length polymorphism linkage maps from the progeny. The paternal map contains 309 markers distributed among 29 linkage groups, with a corrected map distance of 1167 cM (logarithm of the odds [LOD] = 4.0). The average linkage group contained 10.65 +/- 4.85 markers separated by 32.7 +/- 3.8 cM, with statistically significant clustering. The paternal hybrid map had 18.65% more markers than the maternal P. glaucus map, which provides a rough estimate of the extent of genetic differentiation between the species. The maternal map contains 253 markers among 28 linkage groups, without the X and Y chromosomes. Segregation distortion from expected Mendelian ratios was observed for 94/1096 scored loci (8.6%, P < 0.05). The X chromosome map includes 7 markers spanning 29.3 cM (LOD = 3.0). These naturally hybridizing, female heterogametic species are used to study important questions in the maintenance of species boundaries, sex chromosome introgression, sex-limited mimicry, and host plant use. The map will facilitate research into the physiological, ecological, and evolutionary genetics of these phenomena.

Lateral eyes direct principal eyes as jumping spiders track objects.

One way of circumventing the functional tradeoffs on eye design [1,2] is to have different eyes for different tasks. For example, jumping spiders (Salticidae), known for elaborate, visually guided courtship and predatory behavior [3], view the same object simultaneously with two of their four pairs of eyes: the antero-lateral eyes (ALEs) and the principal eyes (reviewed in [2]; Figure 1A). The ALEs, with immobile lenses and retinas, wide fields of view, and hyperacute sensitivity to moving stimuli [4], are structurally distinct from the principal eyes, which have the best spatial acuity known for terrestrial invertebrates and can discern fine details of stationary objects [5]. Behind the immobile corneal lenses of the principal eyes are miniature, boomerang-shaped retinas with correspondingly small fields of view (Figure 1B). The principal-eye visual fields are greatly expanded and overlap because of eye movements: these retinas are at the proximal ends of long, moveable tubes within the spider's cephalothorax [6]. By designing and using a specialized eyetracker, we tested whether principal-eye gaze direction is influenced by what the ALEs see. The principal eyes scanned stationary objects regardless of whether the ALEs were masked, but only when the ALEs were unmasked did the principal eyes smoothly track moving disks. The principal eyes, with high acuity but a narrow field of view, can thus precisely target moving stimuli, but only with the guidance of the secondary eyes.

Persistence and inheritance of costs of resistance to imidacloprid in Colorado potato beetle.

Reduced fitness among resistant versus susceptible individuals slows resistance evolution and makes it easier to manage. A loss of resistance costs could indicate novel adaptations or mutations contributing to resistance. We measured costs of resistance to imidacloprid in a Massachusetts resistant population compared with a Massachusetts susceptible population in 1999 in terms of fecundity, hatching success, egg development time, and sprint speed. Resistance was additive and seemed to be polygenic with high heritability. The fecundity cost appeared overdominant in 1999, and the hatch rate cost was partly recessive in 1999, but neither was significantly different from dominant or recessive. In 2004, we repeated our measures of resistance costs in Massachusetts in terms of fecundity and hatching success, and we added a new resistant population from Maine. In 2005, we compared development time of Maine resistant and the laboratory susceptible colony eggs. Significant fecundity costs of resistance were found in both population in both 1999 and 2004, and significant egg developmental time costs were found in 1999 and 2005. However, the hatching success costs of resistance were significant in 1999 and not apparent in 2004, suggesting some modification or replacement of the resistance genes in the intervening time.

A New Mechanism for Mendelian Dominance in Regulatory Genetic Pathways: Competitive Binding by Transcription Factors.

We report a new mechanism for allelic dominance in regulatory genetic interactions that we call binding dominance. We investigated a biophysical model of gene regulation, where the fractional occupancy of a transcription factor (TF) on the cis-regulated promoter site it binds to is determined by binding energy (-ΔG) and TF dosage. Transcription and gene expression proceed when the TF is bound to the promoter. In diploids, individuals may be heterozygous at the cis-site, at the TF's coding region, or at the TF's own promoter, which determines allele-specific dosage. We find that when the TF's coding region is heterozygous, TF alleles compete for occupancy at the cis-sites and the tighter-binding TF is dominant in proportion to the difference in binding strength. When the TF's own promoter is heterozygous, the TF produced at the higher dosage is also dominant. Cis-site heterozygotes have additive expression and therefore codominant phenotypes. Binding dominance propagates to affect the expression of downstream loci and it is sensitive in both magnitude and direction to genetic background, but its detectability often attenuates. While binding dominance is inevitable at the molecular level, it is difficult to detect in the phenotype under some biophysical conditions, more so when TF dosage is high and allele-specific binding affinities are similar. A body of empirical research on the biophysics of TF binding demonstrates the plausibility of this mechanism of dominance, but studies of gene expression under competitive binding in heterozygotes in a diversity of genetic backgrounds are needed.

An AFLP-based interspecific linkage map of sympatric, hybridizing Colias butterflies.

Colias eurytheme and C. philodice are sister species with broad sympatry in North America. They hybridize frequently and likely share a significant portion of their genomes through introgression. Both taxa have been ecologically well characterized and exploited to address a broad spectrum of evolutionary issues. Using AFLP markers, we constructed the first linkage map of Colias butterflies. The map is composed of 452 markers spanning 2541.7 cM distributed over 51 linkage groups (40 major groups and 11 small groups with 2-4 markers). Statistical tests indicate that these AFLP markers tend to cluster over the map, with the coefficient of variation of interval sizes being 1.236 (95% C.I. is 1.234-1.240). This nonrandom marker distribution can account for the nonequivalence between the number of linkage groups and the actual haploid chromosome number (N = 31). This study presents the initial step for further marker-assisted research on Colias butterflies, including QTL and introgression analyses. Further investigation of the genomes will help us understand better the roles of introgression and natural selection in the evolution of hybridizing species and devise more appropriate strategies to control these pests.

Speciation despite gene flow when developmental pathways evolve.

Evolutionary biologists assume that species formation requires a drastic reduction in gene exchange between populations, but the rate sufficient to prevent speciation is unknown. To study speciation, we use a new class of population genetic models that incorporate simple developmental genetic rules, likely present in all organisms, to construct the phenotype. When we allow replicate populations to evolve in parallel to a new, shared optimal phenotype, often their hybrids acquire poorly regulated phenotypes: Dobzhansky-Muller incompatibilities arise and postzygotic reproductive isolation evolves. Here we show that, although gene exchange does inhibit this process, it is the proportion of migrants exchanged (m) rather than the number of migrants (Nm) that is critical, and rates as high as 16 individuals exchanged per generation still permit the evolution of postzygotic isolation. Stronger directional selection counters the inhibitory effect of gene flow, increasing the speciation probability. We see similar results when populations in a standard two-locus, two-allele Dobzhansky-Muller model are subject to simultaneous directional selection and gene flow. However, in developmental pathway models with more than two loci, gene flow is more able to impede speciation. Genetic incompatibilities arise as frequent by-products of adaptive evolution of traits determined by regulatory pathways, something that does not occur when phenotypes are modeled using the standard, additive genetic framework. Development therefore not only constrains the microevolutionary process, it also facilitates the interactions among genes and gene products that make speciation more likely-even in the face of strong gene flow.

Hybrid incompatibility despite pleiotropic constraint in a sequence-based bioenergetic model of transcription factor binding.

Hybrid incompatibility can result from gene misregulation produced by divergence in trans-acting regulatory factors and their cis-regulatory targets. However, change in trans-acting factors may be constrained by pleiotropy, which would in turn limit the evolution of incompatibility. We employed a mechanistically explicit bioenergetic model of gene expression wherein parameter combinations (number of transcription factor molecules, energetic properties of binding to the regulatory site, and genomic background size) determine the shape of the genotype-phenotype (G-P) map, and interacting allelic variants of mutable cis and trans sites determine the phenotype along that map. Misregulation occurs when the phenotype differs from its optimal value. We simulated a pleiotropic regulatory pathway involving a positively selected and a conserved trait regulated by a shared transcription factor (TF), with two populations evolving in parallel. Pleiotropic constraints shifted evolution in the positively selected trait to its cis-regulatory locus. We nevertheless found that the TF genotypes often evolved, accompanied by compensatory evolution in the conserved trait, and both traits contributed to hybrid misregulation. Compensatory evolution resulted in "developmental system drift," whereby the regulatory basis of the conserved phenotype changed although the phenotype itself did not. Pleiotropic constraints became stronger and in some cases prohibitive when the bioenergetic properties of the molecular interaction produced a G-P map that was too steep. Likewise, compensatory evolution slowed and hybrid misregulation was not evident when the G-P map was too shallow. A broad pleiotropic "sweet spot" nevertheless existed where evolutionary constraints were moderate to weak, permitting substantial hybrid misregulation in both traits. None of these pleiotropic constraints manifested when the TF contained nonrecombining domains independently regulating the respective traits.

Hybrid incompatibility arises in a sequence-based bioenergetic model of transcription factor binding.

Postzygotic isolation between incipient species results from the accumulation of incompatibilities that arise as a consequence of genetic divergence. When phenotypes are determined by regulatory interactions, hybrid incompatibility can evolve even as a consequence of parallel adaptation in parental populations because interacting genes can produce the same phenotype through incompatible allelic combinations. We explore the evolutionary conditions that promote and constrain hybrid incompatibility in regulatory networks using a bioenergetic model (combining thermodynamics and kinetics) of transcriptional regulation, considering the bioenergetic basis of molecular interactions between transcription factors (TFs) and their binding sites. The bioenergetic parameters consider the free energy of formation of the bond between the TF and its binding site and the availability of TFs in the intracellular environment. Together these determine fractional occupancy of the TF on the promoter site, the degree of subsequent gene expression and in diploids, and the degree of dominance among allelic interactions. This results in a sigmoid genotype-phenotype map and fitness landscape, with the details of the shape determining the degree of bioenergetic evolutionary constraint on hybrid incompatibility. Using individual-based simulations, we subjected two allopatric populations to parallel directional or stabilizing selection. Misregulation of hybrid gene expression occurred under either type of selection, although it evolved faster under directional selection. Under directional selection, the extent of hybrid incompatibility increased with the slope of the genotype-phenotype map near the derived parental expression level. Under stabilizing selection, hybrid incompatibility arose from compensatory mutations and was greater when the bioenergetic properties of the interaction caused the space of nearly neutral genotypes around the stable expression level to be wide. F2's showed higher hybrid incompatibility than F1's to the extent that the bioenergetic properties favored dominant regulatory interactions. The present model is a mechanistically explicit case of the Bateson-Dobzhansky-Muller model, connecting environmental selective pressure to hybrid incompatibility through the molecular mechanism of regulatory divergence. The bioenergetic parameters that determine expression represent measurable properties of transcriptional regulation, providing a predictive framework for empirical studies of how phenotypic evolution results in epistatic incompatibility at the molecular level in hybrids.

Delimiting species boundaries and the conservation genetics of the endangered maritime ringlet butterfly (Coenonympha nipisiquit McDunnough).

Species delimitation is a difficult problem that has implications across organismal biology, yet no single method has proved wholly satisfactory. We tested the utility of combining species-delimitation methods based on phylogeny and gene flow statistics using two parapatric members of the Coenonympha tullia group as an example: the endangered maritime ringlet butterfly (Coenonympha nipisiquit McDunnough) and the common inornate ringlet butterfly (Coenonympha inornata Edwards). We reconstructed the phylogeny of the nearctic C. tullia-group taxa from mitochondrial DNA (mtDNA) sequences (cytochrome oxidase I and mitochondrial control region) to explore the ancestry of the C. nipisiquit lineage within the group. We investigated the extent of gene flow between the two taxa with F-statistics using 587 nuclear amplified fragment length polymorphism markers, accounting for the effect of potential scoring 'collisions' where a marker may represent more than one DNA fragment. Combining species-delimitation methods was especially effective because it uncovered both historical and recent evolutionary patterns. Phylogenetic analysis of mtDNA revealed the early divergence of C. nipisiquit from other C. tullia-group taxa, including the morphologically similar C. inornata. F-statistics and gene-by-gene introgression profiles demonstrated clear isolation between the two taxa and revealed strong population structure within C. nipisiquit. C. nipisiquit is the first taxon in the nearctic C. tullia group showing strong evidence of genetic isolation. The methods we used are relatively inexpensive and can be widely used to delimit taxonomic boundaries near the species level, both generally and in particular for taxa that may be targets of conservation efforts.

Evolution of branched regulatory genetic pathways: directional selection on pleiotropic loci accelerates developmental system drift.

Developmental systems are regulated by a web of interacting loci. One common and useful approach in studying the evolution of development is to focus on classes of interacting elements within these systems. Here, we use individual-based simulations to study the evolution of traits controlled by branched developmental pathways involving three loci, where one locus regulates two different traits. We examined the system under a variety of selective regimes. In the case where one branch was under stabilizing selection and the other under directional selection, we observed "developmental system drift": the trait under stabilizing selection showed little phenotypic change even though the loci underlying that trait showed considerable evolutionary divergence. This occurs because the pleiotropic locus responds to directional selection and compensatory mutants are then favored in the pathway under stabilizing selection. Though developmental system drift may be caused by other mechanisms, it seems likely that it is accelerated by the same underlying genetic mechanism as that producing the Dobzhansky-Muller incompatibilities that lead to speciation in both linear and branched pathways. We also discuss predictions of our model for developmental system drift and how different selective regimes affect probabilities of speciation in the branched pathway system.

A test for deviation from island-model population structure.

The neutral island model forms the basis for several estimation models that relate patterns of genetic structure to microevolutionary processes. Estimates of gene flow are often based on this model and may be biased when the model's assumptions are violated. An appropriate test for violations is to compare FST scores for individual loci to a null distribution based on the average FST taken over multiple loci. A parametric bootstrap method is described here based on Wright's beta-distribution to generate null distributions of FST for each locus. These null distributions account for error introduced by sampling populations, individuals and loci, and also biological sources of error, including variable alleles/locus and inbreeding. Confidence limits can be obtained directly from these distributions. Significant deviations from the island model may be the result of selection, deviations from the island model's migration pattern, nonequilibrium conditions, or other deviations from island-model assumptions. Only strong biases are likely to be detected because of the inherently large sampling variation of FST. Nevertheless, a coefficient, Nb, describing bias in the spread of the beta-distribution in units comparable to the gene flow parameter, Nm, can be obtained for each locus. In samples from populations of the butterfly Coenonympha tullia, the loci Idh-1, Mdh-1, Pgi and Pgm showed significantly lower FST than expected.

THE PONTIA DAPLIDICE-ED USA HYBRID ZONE IN NORTHWESTERN ITALY.

The pierid butterflies Pontia daplidice and P. edusa, parapatrically distributed in southern Europe, have very similar morphologies and life histories, but show fixed differences at four allozyme markers. We sampled these allozymes in a 28-population transect north of Genoa in Italy, through the hybrid zone where these taxa meet. We used the numerical techniques developed for hybrid zone analysis to study the patterns of genetic differentiation and their underlying evolutionary causes. The hybrid zone is characterized by a very short and steep central region, flanked by broad tails of introgression extended up to 100 km in either direction. From mean two-locus disequilibium of D = 0.148 (maximum-likelihood two-unit support limits 0.139-0.153), and after accounting for minor differences in the center locations of the single-locus clines, which act to bias the dispersal estimate, we estimated a dispersal rate of σ = 4.4 (3.7-5.5) km/gen1/2 . The effective selection needed to maintain the steep central portion is strong, 0.47 < s∗ < 0.64, when combined over potential intrinsic (genetic background) and extrinsic (ecological) sources of selection. The clines in allozyme loci showed variation that was significantly different between the most divergent shapes, and the differences are attributable to different degrees of introgression on the edusa side of the zone. The average selection acting on individual allozyme loci was high at s∗∗∗e  1.5%, but because of the narrowness of the central region of the cline, we suspect that this estimate is somewhat biased by selection on loci closely linked to the allozyme markers. A common question for taxa that show fixed allozyme differences in parapatry is whether or not they are genetically isolated. A fairly general measure of genetic isolation across hybrid zones is the time, T, that it takes a neutral allele to cross the hybrid zone and recombine into the opposite genetic background, given by T = (ß/σ)2 , where ß is the barrier strength of the hybrid zone. Genetic isolation in the Pontia zone is weak, with T  25 generations for most allozyme markers. By this measure, populations of daplidice and edusa on opposite sides of the hybrid zone share more identical-by-descent alleles than do populations of phenotypically pure daplidice in, say, France and Morocco. Accordingly, we think it best for systematists to consider edusa as a well-marked subspecies of P. daplidice.