Genetic evidence for assortative mating between 13-year cicadas and sympatric "17-year cicadas with 13-year life cycles" provides support for allochronic speciation

Abstract. Thirteen-year cicadas of brood XIX from northern Arkansas, Missouri, and southern Illinois (lineage A) are known to be genetically different at two marker loci (mitochondrial DNA and abdominal color) from 13-year cicadas to the south (lineage B) that emerge in the same year. Because 17-year cicadas from all broods (year classes) are indistinguishable from lineage A at these two marker loci, previous workers suggested that the lineage A cicadas of 13-year brood XIX were derived from 17-year cicadas by life-cycle switching (allochrony). Data presented here show that, over the same northern geographic range, lineage A is also present in 13-year cicadas belonging to brood XXIII (which always emerges four years later than brood XIX). Detailed sampling along the putative life-cycle-switching boundary in 13-year brood XXIII revealed a previously unsuspected broad zone of overlap where populations contained individuals of both lineages A and B. Despite this sympatry, and previous reports of a lack of behavioral barriers to interbreeding, a strong correlation between mitochondrial haplotype and abdominal color suggests that assortative mating has taken place. Lineage A 13-year cicadas from both broods XIX and XXIII are only found within a gap in the spatial distribution of 17-year cicadas. This, in combination with the lack of differentiation between lineage A 13- and 17-year cicadas at the marker loci and new behavioral data for 13-year brood XIX, suggests a recent derivation of all northern 13-year cicadas from the 17-year cicadas via life-cycle switching. We discuss the implications of these allochronic shifts for speciation.     

Reproductive character displacement and speciation in periodical cicadas, with description of new species, 13-year Magicicada neotredecem

Abstract. Acoustic mate-attracting signals of related sympatric, synchronic species are always distinguishable, but those of related allopatric species sometimes are not, thus suggesting that such signals may evolve to "reinforce" premating species isolation when similar species become sympatric. This hypothesis predicts divergences restricted to regions of sympatry in partially overlapping species, but such "reproductive character displacement" has rarely been confirmed. We report such a case in the acoustic signals of a previously unrecognized 13-year periodical cicada species, Magicicada neotredecim , described here as a new species (see Appendix). Where M. neotredecim overlaps M. tredecim in the central United States, the dominant male call pitch (frequency) of M. neotredecim increases from approximately 1.4 kHz to 1.7 kHz, whereas that of M. tredecim remains comparatively stable. The average preferences of female M. neotredecim for call pitch show a similar geographic pattern, changing with the call pitch of conspecific males. Magicicada neotredecim differs from 13-year M. tredecim in abdomen coloration, mitochondrial DNA, and call pitch, but is not consistently distinguishable from 17-year M. septendecim ; thus, like other Magicicada species, M. neotredecim appears most closely related to a geographically adjacent counterpart with the alternative life cycle. Speciation in Magicicada may be facilitated by life-cycle changes that create temporal isolation, and reinforcement could play a role by fostering divergence in premating signals prior to speciation. We present two theories of Magicicada speciation by life-cycle evolution: "nurse-brood facilitation" and "life-cycle canalization."     

Genetic mosaic in a marine species flock

We used molecular approaches to study the status of speciation in coral reef fishes known as hamlets (Serranidae: Hypoplectrus). Several hamlet morphospecies coexist on Caribbean reefs, and mate assortatively with respect to their strikingly distinct colour patterns. We provide evidence that, genetically, the hamlets display characteristics common in species flocks on land and in freshwaters. Substitutions within two mitochondrial DNA (mtDNA) protein-coding genes place hamlets within a monophyletic group relative to members of two related genera (Serranus and Diplectrum), and establish that the hamlet radiation must have been very recent. mtDNA distances separating hamlet morphospecies were slight (0.6 +/- 0.04%), yielding a coalescent estimate for the age of the hamlet flock of approximately 430 000 years. Morphospecies did not sort into distinct mtDNA haplotype phylogroups, and alleles at five hypervariable microsatellite loci were shared broadly across species boundaries. None the less, molecular variation was not distributed at random. Analyses of mtDNA haplotype frequencies and nested clades in haplotype networks revealed significant genetic differences between geographical regions and among colour morphospecies. We also observed significant microsatellite differentiation between geographical regions and in Puerto Rico, among colour morphospecies; the latter providing evidence for reproductive isolation between colour morphospecies at this locale. In our Panama collection, however, colour morphospecies were mostly genetically indistinguishable. This mosaic pattern of DNA differentiation implies a complex interaction between population history, mating behaviour and geography and suggests that porous boundaries separate species in this flock of brilliantly coloured coral reef fishes.     

Reconstructing asymmetrical reproductive character displacement in a periodical cicada contact zone

Selection against costly reproductive interactions can lead to reproductive character displacement (RCD). We use information from patterns of displacement and inferences about predisplacement character states to investigate causes of RCD in periodical cicadas. The 13-year periodical cicada Magicicada neotredecim exhibits RCD and strong reproductive isolation in sympatry with a closely related 13-year species, Magicicada tredecim. Displacement is asymmetrical, because no corresponding pattern of character displacement exists within M. tredecim. Results from playback and hybridization experiments strongly suggest that sexual interactions between members of these species were possible at initial contact. Given these patterns, we evaluate potential sources of selection for displacement. One possible source is 'acoustical interference', or mate-location inefficiencies caused by the presence of heterospecifics. Acoustical interference combined with the species-specificity of song pitch and preference appears to predict the observed asymmetrical pattern of RCD in Magicicada. However, acoustical interference does not appear to be a complete explanation for displacement in Magicicada, because our experiments suggest a significant potential for direct sexual interactions between these species before displacement. Another possible source of selection for displacement is hybrid failure. We evaluate the attractiveness of inferred hybrid mating signals, and we examine the viability of hybrid eggs. Neither of these shows strong evidence of hybrid inferiority. We conclude by presenting a model of hybrid failure related to life cycle differences in Magicicada.     

Incipient allochronic speciation by climatic disruption of the reproductive period

Disruptive selection of life-cycle timing may cause temporal isolation directly and, ultimately, allochronic speciation. Despite the fact that segregation of the reproductive period among related species has been broadly observed across taxa, it remains controversial whether temporal isolation can function as the primary process of speciation. In the Japanese winter geometrid moth Inurois punctigera, allochronic divergence has resulted from climatic disruption of the reproductive period. In habitats with severe midwinter, two sympatric groups of moth reproduce allochronically in early and late winter. These groups are genetically diverging sister lineages and now co-occur allochronically throughout Japan. By contrast, in habitats with milder midwinter these lineages form a continuous adult period and gene flow has been facilitated between the lineages. These results, together with the fact that there is no difference in larval host use, indicate that temporal isolation has been the sole mechanism for allochronic isolation in colder habitats and that allochrony is not a by-product of other adaptations. Thus, the allochronic divergence of sympatric I. punctigera populations represents an incipient speciation process driven by midwinter disruption of the reproductive period.     

Evidence for low‐level hybridization between two allochronic populations of the pine processionary moth,Thaumetopoea pityocampa (Lepidoptera: Notodontidae)

Divergence between populations sharing the same habitat can be initiated by different reproductive times, leading to allochronic differentiation. A spatially localized allochronic summer population (SP) of the pine processionary moth Thaumetopoea pityocampa, recently discovered in Portugal, occurs in sympatry with the local winter population (WP). We examined the level of genetic differentiation between the two populations and estimated the current gene flow within the spatial framework of their co‐occurrence. Mitochondrial data indicated that the two sympatric populations were genetically closer than other WP populations. Conversely, microsatellite genotyping uncovered greater differentiation between the two sympatric populations than between allopatric ones. While male trapping confirmed that reproduction of SP and WP occurred at distinct times, clustering approaches demonstrated the presence of a few LateSP individuals emerging within the WP flight period, although genetically identified as SP. We also identified rare recent hybridization events apparently occurring mainly in the margins of the current SP range. The ongoing gene flow detected between the ancestral and the emerging allochronic populations revealed an incomplete reproductive isolation, which must therefore be taken into account and integrated with studies focussed on ecological drivers, so that a complete understanding of the ongoing speciation process might be achieved.     

Molecular and morphological divergence in the butterfly genus Lycaeides (Lepidoptera: Lycaenidae) in North America: evidence of recent speciation

Male genital morphology, allozyme allele frequencies and mtDNA sequence variation were surveyed in the butterfly species Lycaeides idas and L. melissa from across much of their range in North America. Despite clear differences in male genital morphology, wing colour patterns and habitat characteristics, genetic variation was not taxonomically or geographically structured and the species were not identifiable by either genetic data set. Genetic distances (Nei's D=0.002–0.078, calculated from allozyme data) between all populations of both species were within the range commonly observed for conspecific populations of other butterflies. The most frequent mtDNA haplotype was present in individuals of both species in populations from southern California to Wisconsin. We conclude that speciation has probably happened recently and the lack of genetic differentiation between the species is the product of either (1) recent or ongoing gene flow at neutral loci, and/or (2) an insufficiency of time for lineage sorting. The evolution of male genital morphology, wing colour patterns and ecological characteristics has proceeded more rapidly than allozyme or mtDNA evolution.     

Life cycle replacement by gene introduction under an allee effect in periodical cicadas

Periodical cicadas (Magicicada spp.) in the USA are divided into three species groups (-decim, -cassini, -decula) of similar but distinct morphology and behavior. Each group contains at least one species with a 17-year life cycle and one with a 13-year cycle; each species is most closely related to one with the other cycle. One explanation for the apparent polyphyly of 13- and 17-year life cycles is that populations switch between the two cycles. Using a numerical model, we test the general feasibility of life cycle switching by the introduction of alleles for one cycle into populations of the other cycle. Our results suggest that fitness reductions at low population densities of mating individuals (the Allee effect) could play a role in life cycle switching. In our model, if the 13-year cycle is genetically dominant, a 17-year cycle population will switch to a 13-year cycle given the introduction of a few 13-year cycle alleles under a moderate Allee effect. We also show that under a weak Allee effect, different year-classes ("broods") with 17-year life cycles can be generated. Remarkably, the outcomes of our models depend only on the dominance relationships of the cycle alleles, irrespective of any fitness advantages.     

Incipient allochronic speciation in the pine processionary moth (Thaumetopoea pityocampa, Lepidoptera, Notodontidae)

A plausible case of allochronic differentiation, where barrier to gene flow is primarily attributable to a phenological shift, was recently discovered in Portugal for the pine processionary moth Thaumetopoea pityocampa. Previous results suggested that the observed 'summer population' (SP) originated from the sympatric winter population (WP). Our objectives were to finely analyse these patterns and test their stability in time, through field monitoring and genetic analyses of larvae and adults across different years. Reproductive activity never overlapped between SP and WP. Microsatellites showed a clear differentiation of the SP, consistent with a strong reduction in gene flow owing to the phenological shift. Assignment tests suggested that some individuals shift from the SP to the WP phenology, causing some hybridization. We discuss these patterns and their maintenance over time. This could be a first stage of allochronic speciation, and SP should be considered as a distinct phenological race.     

Chronic speciation in periodical cicadas

Allochronic speciation and reproductive character displacement are two intuitively attractive models of speciation. The first proposes that changes in the timing of life cycles produce new species, whereas the second suggests that speciation is due to the exaggeration of sexual isolation in sympatric populations. Clear examples of either process in nature remain elusive, despite some extensive searches. Two recent studies of mtDNA markers and behaviour of periodical cicadas in North America have identified a new species of cicada that seems to provide good evidence for the involvement of both processes in its origin.     

SIGNIFICANT ROLE FOR HISTORICAL EFFECTS IN THE EVOLUTION OF REPRODUCTIVE ISOLATION: EVIDENCE FROM PATTERNS OF INTROGRESSION BETWEEN THE CYPRINID FISHES,LUXILUS CORNUTUS AND LUXILUS CHRYSOCEPHALUS

Samples of Luxilus cornutus, Luxilus chrysocephalus, and their hybrids were collected along hypothesized routes of dispersal from Pleistocene refugia to examine the significance of geographic variation in patterns of introgression between these species. Patterns of allozyme and mitochondrial DNA (mtDNA) variation were generally consistent with those from previous studies. Tests of Hardy-Weinberg equilibrium revealed significant deficiencies of heterozygotes in all samples, indicating some form of reproductive isolation. Mitochondrial DNAs of each species were not equally represented in F₁ hybrids; however, this bias was eliminated when the two largest samples were excluded from the analysis. Backcross hybrids exhibited biased mtDNA introgression, as samples from Lake Erie (eastern) and Lake Michigan (western) drainages showed significant excesses of mtDNAs from L. chrysocephalus and L. cornutus, respectively, relative to frequencies of diagnostic allozyme markers. The extent and direction of allozyme and mtDNA introgression was quantified by calculating isolation index values from morphologically “pure” individuals of each species from each locality. Analysis of variance of these measures identified limited introgression of allozyme variants with no geographic pattern, but significant differences in direction of mtDNA introgression between drainages (i.e., postglacial dispersal route). Association between patterns of mtDNA introgression and dispersal route across the latitudinal width of the contact zone is best explained by genetic divergence during past isolation of ancestral populations from these drainages. These results identify a significant role for historical effects in the evolution of reproductive isolation and the process of speciation.     

SELECTION FOR PRIME-NUMBER INTERVALS IN A NUMERICAL MODEL OF PERIODICAL CICADA EVOLUTION

Periodical cicadas are known for unusually long and prime-numbered life cycles (13 and 17 years) for insects. To explain the evolution of prime-numbered reproductive intervals (life cycles), the hybridization hypothesis claims that prime numbers greatly reduce the chance of hybridization with other life cycles. We investigate the hybridization hypothesis using a simulation model. This model is a deterministic, discrete population model with three parameters: larval survival per year, clutch size, and emergence success. Reproductive intervals from 10 years to 20 years compete for survival in the simulations. The model makes three key assumptions: a Mendelian genetic system, random mating among broods of different life-cycle lengths, and integer population sizes. Longer life cycles have larger clutch sizes but suffer higher total mortality than shorter life cycles. Our results show that (1) nonprime-numbered reproductive intervals disappear rapidly in comparison to the selection among the various prime-numbered life cycles, (2) the selection of prime-numbered intervals happens only when populations are at the verge of extinction, and (3) the 13- and 17-year prime phenotypes evolve under certain conditions of the model and may coexist. The hybridization hypothesis is discussed in light of other hypotheses for the evolution of periodical cicada life cycles.     

Permanence induced by life-cycle resonances: the periodical cicada problem

Periodical cicadas are known for their unusually long life cycle for insects and their prime periodicity of either 13 or 17 years. One of the explanations for the prime periodicity is that the prime periods are selected to prevent cicadas from resonating with predators with submultiple periods. This paper considers this hypothesis by investigating a population model for periodical predator and prey. The study shows that if the periods of the two periodical species are not coprime, then the predator cannot resist the invasion of the prey. On the other hand, if the periods are coprime, then the predator can resist the invasion of the prey. It is also shown that if the periods are not coprime, then the life-cycle resonance can induce a permanent system, in which no cohorts are missing in both populations. On the other hand, if the periods are coprime, then the system cannot be permanent.     

Genetic differentiation and reproductive isolation in a Cyprinodon fish species flock from Laguna Chichancanab, Mexico

The Cyprinodon species flock from Laguna Chichancanab, aged 8000 years, provides another potential case of sympatric speciation. The flock consists of seven morphologically distinct species, each within partially different trophic niches, and a group of specimens which cannot unequivocally be assigned to one of these species. Genetic analyses, based on mtDNA and five microsatellite loci, revealed significant genetic differentiation of one species, C. maya, from other members of the species flock, providing strong evidence for reproductive isolation. For the remaining members of the flock significant genetic structuring was detected, with some evidence of gene flow with the most abundant species C. beltrani. These analyses suggest that speciation proceeds with ongoing hybridisation, and further suggest that the morphologically unidentifiable specimens found in the lake are probably hybrids. I propose that in the Cyprinodon species flock besides disruptive selection sexual selection plays an important part in achieving and maintaining reproductive isolation.     

Permanence induced by life-cycle resonances: the periodical cicada problem

Periodical cicadas are known for their unusually long life cycle for insects and their prime periodicity of either 13 or 17 years. One of the explanations for the prime periodicity is that the prime periods are selected to prevent cicadas from resonating with predators with submultiple periods. This paper considers this hypothesis by investigating a population model for periodical predator and prey. The study shows that if the periods of the two periodical species are not coprime, then the predator cannot resist the invasion of the prey. On the other hand, if the periods are coprime, then the predator can resist the invasion of the prey. It is also shown that if the periods are not coprime, then the life-cycle resonance can induce a permanent system, in which no cohorts are missing in both populations. On the other hand, if the periods are coprime, then the system cannot be permanent.     

Small body size increases the regional differentiation of populations of tropical mantellid frogs (Anura: Mantellidae)

The processes affecting species diversification may also exert an influence on patterns of genetic variability within species. We evaluated the contributions of five variables potentially influencing clade diversification (body size, reproductive mode, range size, microhabitat and skin texture) on mtDNA divergence and polymorphism among populations of 40 species of frogs (Mantellidae) from two rainforest communities in Madagascar. We report an inverse association between body size and nucleotide divergence between populations but find no influence of other variables on genetic variation. Body size explained ca. 11% of the variation in nucleotide divergence between populations and was coupled with high F_ST levels and an absence of haplotype sharing in small‐bodied and medium‐sized frogs. Low dispersal ability is likely the proximate mechanism producing higher population differentiation in small mantellids. The lack of genetic cohesion among populations establishes regional genetic fragmentation which in turn has the potential to accelerate rates of allopatric speciation in small frogs relative to large species. However, there is little evidence of increased speciation rates in these or other small‐bodied organisms. We reconcile these contradictory observations by suggesting that lower dispersal ability also curbs colonization of new areas, decelerating diversification in weak dispersers. Our results imply that the intermediate dispersal model also applies to amphibians and may explain inconsistent previous results on the correlation of body size and speciation rate.     

Review. Hybrid trait speciation and Heliconius butterflies

Homoploid hybrid speciation (HHS) is the establishment of a novel species through introgressive hybridization without a change in chromosome number. We discuss different routes by which this might occur and propose a novel term, 'hybrid trait speciation', which combines the idea that hybridization can generate adaptive novelty with the 'magic trait' model of ecological speciation. Heliconius butterflies contain many putative examples of hybrid colour patterns, but only recently has the HHS hypothesis been tested explicitly in this group. Molecular data has shown evidence for gene flow between many distinct species. Furthermore, the colour pattern of Heliconius heurippa can be recreated in laboratory crosses between Heliconius melpomene and Heliconius cydno and, crucially, plays a role in assortative mating between the three species. Nonetheless, although the genome of H. heurippa shows evidence for hybridization, it is not a mosaic of the two parental species. Instead, ongoing hybridization has likely blurred any signal of the original speciation event. We argue that where hybridization leads to novel adaptive traits that also cause reproductive isolation, it is likely to trigger speciation.     

Divergent sexual selection via male competition: ecology is key

Sexual selection and ecological differences are important drivers of speciation. Much research has focused on female choice, yet the role of male competition in ecological speciation has been understudied. Here, we test how mating habitats impact sexual selection and speciation through male competition. Using limnetic and benthic species of threespine stickleback fish, we find that different mating habitats select differently on male traits through male competition. In mixed habitat with both vegetated and open areas, selection favours two trait combinations of male body size and nuptial colour: large with little colour and small with lots of colour. This matches what we see in reproductively isolated stickleback species, suggesting male competition could promote trait divergence and reproductive isolation. In contrast, when only open habitat exists, selection favours one trait combination, large with lots of colour, which would hinder trait divergence and reproductive isolation. Other behavioural mechanisms in male competition that might promote divergence, such as avoiding aggression with heterospecifics, are insufficient to maintain separate species. This work highlights the importance of mating habitats in male competition for both sexual selection and speciation.     

The role of allochrony in speciation

The importance of sympatric speciation – the evolution of reproductive isolation between codistributed conspecific individuals – in generating biodiversity is highly controversial. Allochrony, or differences in breeding time (phenology) between conspecific individuals, has the potential to lead to reproductive isolation and therefore speciation. We critically review the literature to test the importance of allochronic speciation over the three timescales over which allochrony can occur – over the day, between seasons or between years – and explore what is known about genomic mechanisms underlying allochrony in the diverse taxa in which it is found. We found that allochrony can be a key contributor to reproductive isolation, especially if populations have little overlap in breeding time and therefore little potential for gene flow, and may sometimes be the initial or key driver of speciation. Shifts in phenology can be caused by several factors, including a new ecological opportunity, environmental change, or reinforcement. The underlying genomic basis of allochrony has been studied mostly in insects, highlighting the need for genomic studies in other taxa; nonetheless, results to date indicate that several cases of allochrony involve changes in circadian genes. This review provides the first comprehensive discussion of the role of allochrony in speciation and demonstrates that allochrony as a contributor to divergence may be more widespread than previously thought. Understanding genomic changes and adaptations allowing organisms to breed at new times may be key in the light of phenological changes required under climate change.     

Parallel allochronic divergence in a winter moth due to disruption of reproductive period by winter harshness

The disruption of reproductive timing by climatic harshness may result in the temporal isolation of conspecific populations and, ultimately, in speciation. However, whether temporal isolation alone can act as the force initiating speciation and how often the same type of climatic disruption results in the divergence of allochronic populations in a lineage are largely unknown. The reproductive period of the winter geometrid moth Inurois punctigera is separated into early and late winter in habitats with severe winters, but not in habitats with mild winters, suggesting that the reproductive season is disrupted by the harshness of the mid-winter period. Here, we show that sympatric pairs of early- and late-winter populations that differ in origin exist in different regions, suggesting a parallel divergence of reproductive timing. In each region, significant genetic differentiation exists between these early- and late-winter populations, suggesting that the temporal reproductive isolation has persisted. Moreover, we demonstrate that the temporal isolation, in comparison with geographic isolation, contributes greatly to the genetic differentiation among geographic and temporal populations by an analysis of molecular variance and by a comparison of genetic differentiations (F(ST) ) between geographic populations with and without difference in reproductive season. Our results suggest that adaptive divergence of allochronically reproducing populations has occurred independently in different regions, implying the generality of the role of temporal isolation in initiating speciation in a winter moth lineage.