Strong divergent selection at multiple loci in two closely related species of ragworts adapted to high and low elevations on Mount Etna

Recently diverged species present particularly informative systems for studying speciation and maintenance of genetic divergence in the face of gene flow. We investigated speciation in two closely related Senecio species, S. aethnensis and S. chrysanthemifolius, which grow at high and low elevations, respectively, on Mount Etna, Sicily and form a hybrid zone at intermediate elevations. We used a newly generated genome‐wide single nucleotide polymorphism (SNP) dataset from 192 individuals collected over 18 localities along an elevational gradient to reconstruct the likely history of speciation, identify highly differentiated SNPs, and estimate the strength of divergent selection. We found that speciation in this system involved heterogeneous and bidirectional gene flow along the genome, and species experienced marked population size changes in the past. Furthermore, we identified highly‐differentiated SNPs between the species, some of which are located in genes potentially involved in ecological differences between species (such as photosynthesis and UV response). We analysed the shape of these SNPs’ allele frequency clines along the elevational gradient. These clines show significantly variable coincidence and concordance, indicative of the presence of multifarious selective forces. Selection against hybrids is estimated to be very strong (0.16–0.78) and one of the highest reported in literature. The combination of strong cumulative selection across the genome and previously identified intrinsic incompatibilities probably work together to maintain the genetic and phenotypic differentiation between these species – pointing to the importance of considering both intrinsic and extrinsic factors when studying divergence and speciation.     

Chromosomal variation segregates within incipient species and correlates with reproductive isolation

Reproductive isolation is a critical step in the process of speciation. Among the most important factors driving reproductive isolation are genetic incompatibilities. Whether these incompatibilities are already present before extrinsic factors prevent gene flow between incipient species remains largely unresolved in natural systems. This question is particularly challenging because it requires that we catch speciating populations in the act before they reach the full‐fledged species status. We measured the extent of intrinsic postzygotic isolation within and between phenotypically and genetically divergent lineages of the wild yeast Saccharomyces paradoxus that have partially overlapping geographical distributions. We find that hybrid viability between lineages progressively decreases with genetic divergence. A large proportion of postzygotic inviability within lineages is associated with chromosomal rearrangements, suggesting that chromosomal differences substantially contribute to the early steps of reproductive isolation within lineages before reaching fixation. Our observations show that polymorphic intrinsic factors may segregate within incipient species before they contribute to their full reproductive isolation and highlight the role of chromosomal rearrangements in speciation. We propose different hypotheses based on adaptation, biogeographical events and life history evolution that could explain these observations.     

Along the speciation continuum: Quantifying intrinsic and extrinsic isolating barriers across five million years of evolutionary divergence in California jewelflowers

Understanding the relative roles of intrinsic and extrinsic reproductive barriers, and their interplay within the geographic context of diverging taxa, remains an outstanding challenge in the study of speciation. We conducted a comparative analysis of reproductive isolation in California Jewelflowers (Streptanthus, s.l., Brassicaceae) by quantifying potential barriers to gene flow at multiple life history stages in 39 species pairs spanning five million years of evolutionary divergence. We quantified nine potential pre‐ and postzygotic barriers and explored patterns of reproductive isolation in relation to genetic distance. Intrinsic postzygotic isolation was initially weak, increased at intermediate genetic distances, and reached a threshold characterized by complete genetic incompatibility. Climatic niche differences were strong at shallow genetic distances, and species pairs with overlapping ranges showed slight but appreciable phenological isolation, highlighting the potential for ecological barriers to contribute to speciation. Geographic analyses suggest that speciation is not regionally allopatric in the California Jewelflowers, as recently diverged taxa occur in relatively close proximity and display substantial range overlap. Young pairs are characterized by incomplete intrinsic postzygotic isolation, suggesting that extrinsic barriers or fine‐scale spatial segregation are more important early in the divergence process than genetic incompatibilities.     

The genomic bases of morphological divergence and reproductive isolation driven by ecological speciation in Senecio (Asteraceae)

Ecological speciation, driven by adaptation to contrasting environments, provides an attractive opportunity to study the formation of distinct species, and the role of selection and genomic divergence in this process. Here, we focus on a particularly clear‐cut case of ecological speciation to reveal the genomic bases of reproductive isolation and morphological differences between closely related Senecio species, whose recent divergence within the last ~200 000 years was likely driven by the uplift of Mt. Etna (Sicily). These species form a hybrid zone, yet remain morphologically and ecologically distinct, despite active gene exchange. Here, we report a high‐density genetic map of the Senecio genome and map hybrid breakdown to one large and several small quantitative trait loci (QTL). Loci under diversifying selection cluster in three 5 cM regions which are characterized by a significant increase in relative (F_ST), but not absolute (d_XY), interspecific differentiation. They also correspond to some of the regions of greatest marker density, possibly corresponding to ‘cold‐spots’ of recombination, such as centromeres or chromosomal inversions. Morphological QTL for leaf and floral traits overlap these clusters. We also detected three genomic regions with significant transmission ratio distortion (TRD), possibly indicating accumulation of intrinsic genetic incompatibilities between these recently diverged species. One of the TRD regions overlapped with a cluster of high species differentiation, and another overlaps the large QTL for hybrid breakdown, indicating that divergence of these species may have occurred due to a complex interplay of ecological divergence and accumulation of intrinsic genetic incompatibilities.     

Climatic adaptation and ecological divergence between two closely related pine species in Southeast China

Climate is one of the most important drivers for adaptive evolution in forest trees. Climatic selection contributes greatly to local adaptation and intraspecific differentiation, but this kind of selection could also have promoted interspecific divergence through ecological speciation. To test this hypothesis, we examined intra‐ and interspecific genetic variation at 25 climate‐related candidate genes and 12 reference loci in two closely related pine species, Pinus massoniana Lamb. and Pinus hwangshanensis Hisa, using population genetic and landscape genetic approaches. These two species occur in Southeast China but have contrasting ecological preferences in terms of several environmental variables, notably altitude, although hybrids form where their distributions overlap. One or more robust tests detected signals of recent and/or ancient selection at two‐thirds (17) of the 25 candidate genes, at varying evolutionary timescales, but only three of the 12 reference loci. The signals of recent selection were species specific, but signals of ancient selection were mostly shared by the two species likely because of the shared evolutionary history. F_ST outlier analysis identified six SNPs in five climate‐related candidate genes under divergent selection between the two species. In addition, a total of 24 candidate SNPs representing nine candidate genes showed significant correlation with altitudinal divergence in the two species based on the covariance matrix of population history derived from reference SNPs. Genetic differentiation between these two species was higher at the candidate genes than at the reference loci. Moreover, analysis using the isolation‐with‐migration model indicated that gene flow between the species has been more restricted for climate‐related candidate genes than the reference loci, in both directions. Taken together, our results suggest that species‐specific and divergent climatic selection at the candidate genes might have counteracted interspecific gene flow and played a key role in the ecological divergence of these two closely related pine species.     

STRUCTURE OF AN ASYMMETRIC HYBRID ZONE BETWEEN TWO WATER STRIDER SPECIES (HEMIPTERA: GERRIDAE: LIMNOPORUS)

The water stricter species Limnoporus dissortis and L. notabilis hybridize across a broad zone in western Canada. Body length and alleles at four allozyme loci show a steep cline along the east slope of the Rocky Mountains in western Alberta, while in central British Columbia the parental phenotypes coexist without merging fully. One sex-linked locus shows little introgression, while there is apparently considerable gene flow at three autosomal loci. Although the hybrid zone has characteristics of a broad tension zone, the spatial distribution of introgression suggests that habitat patchiness and differential habitat associations of the two species also contribute to the pattern of hybridization. Asymmetry in interspecific mating success and incompatibilities of sex chromosomes with each other or with cytoplasmic factors appear to account for the occurrence of L. dissortis genotypes within the range of L. notabilis, and the lack of L. notabilis genotypes within the range of L. dissortis. The genetic structure of this hybrid zone supports the importance of sex-linked traits in maintaining the integrity of species, while its spatial structure suggests that extrinsic habitat features can combine with intrinsic genetic incompatibilities to produce complex hybrid interactions.     

COMPONENTS OF REPRODUCTIVE ISOLATION BETWEEN ORCHIS MASCULA AND ORCHIS PAUCIFLORA

Studies of the strength and nature of reproductive isolation (RI) between species can greatly contribute to our understanding of speciation. Although the role of RI in speciation is well recognized, there is a dearth of information on the contributions of different barriers between related plant species. Here, we estimated multiple components of RI between two Mediterranean orchid sister species (Orchis mascula and Orchis pauciflora), disentangling the strength and absolute contributions of seven different isolating mechanisms. Our survey includes one prepollination, two postpollination prezygotic (pollen–stigma incompatibility, conspecific pollen precedence), two intrinsic postzygotic (embryo mortality and hybrid sterility) and two extrinsic postzygotic (hybrid habitat differentiation and hybrid pollination) isolating mechanisms. We found strong RI between the investigated species, although none of the barriers were able to completely impede gene flow. Five isolating mechanisms contributed positively to the maintenance of species boundaries. Contrary to most surveys of isolating mechanisms, our data speak against a clear predominance of prepollination or of prezygotic barriers but confirm the emerging pattern of multiple barriers contributing to the maintenance of species integrity. These findings suggest an allopatric condition during early phases of species divergence. We discuss our data in the wider context of previous studies carried out in this orchid group by using a comparative approach.     

Evolutionary history of the scorpionfly Dicerapanorpa magna (Mecoptera,Panorpidae)

Climate changes can have fundamental impacts on the distributional patterns of montane species, and range shifts frequently lead to allopatric divergence followed by the establishment of secondary contact zones. Many European and North American organisms have retreated to southern refugia during glacial periods and colonized northward during postglacial periods, but little is known about the evolutionary response of cold‐adapted insects to Pleistocene climate changes in eastern Asia. The scorpionfly Dicerapanorpa magna (Chou), with cold temperate habitat preference and weak dispersal ability, provides a good model system to explore how climate changes have influenced the distribution and divergence of cold‐adapted insects in eastern Asia. This study reconstructed the demographic dynamics and evolutionary history of D. magna with phylogeographic approaches, and predicted the species’ suitable areas under the Last Glacial Maximum (LGM) and current scenarios with the ecological niche modelling analysis. The mitochondrial cytochrome c oxidase subunit I resolved three phylogenetic lineages in D. magna dating back to Pleistocene, corresponding well with the geographically isolated Qinling, Bashan and Minshan Mountains. The ecological niche modelling recovered the suitable habitats for D. magna were the Qinling and Bashan Mountains under LGM and current conditions. The three lineages of D. magna might be in a process of incipient speciation, and likely derived their current distribution from separate glacial origins, followed by vicariance and divergence.     

The genetic basis of intrinsic and extrinsic post-zygotic reproductive isolation jointly promoting speciation in the lake whitefish species complex (Coregonus clupeaformis)

Understanding the genetic architecture of reproductive barriers and the evolutionary forces that drove their divergence represents a considerable challenge towards understanding speciation. The objective of this study was to determine the genetic basis of intrinsic and extrinsic post-zygotic isolation in diverging populations of dwarf and normal lake whitefish with allopatric glacial origins. We found that the rate of embryonic mortality was 5.3-6.5 times higher in dwarf-normal hybrid backcrosses during development than in F1 dwarf and normal crosses. When comparing embryos that died during development against larvae that successfully hatched, patterns of Mendelian segregation at 101 loci whose linkage is known identified 13 loci distributed over seven linkage groups that exhibited significant shifts in segregation ratios leading to significant segregation distortion at these loci in the surviving progeny. Controlled crosses and quantitative trait loci analysis revealed a significant genetic basis for developmental time until emergence, a trait critical to fish larval survival in nature. Hatching backcross progeny exhibited asynchronous emergence and transgressive segregation, suggesting that extrinsic post-zygotic isolation may select against hybridization in specific environmental contexts. Evidence of a genetic basis for increased embryonic mortality followed by asynchronous emergence indicated that intrinsic and extrinsic mechanisms are not mutually exclusive in the formation and maintenance of reproductive isolation, but may be jointly promoting population divergence and ultimately speciation.     

Speciation history of three closely related pines Pinus mugo (T.), P. uliginosa (N.) and P. sylvestris (L.)

Nucleotide polymorphisms at genomic regions including 17 nuclear loci, two chloroplast and one mitochondrial DNA fragments were used to study the speciation history of three pine species: dwarf mountain pine (Pinus mugo), peat‐bog pine (P. uliginosa) and Scots pine (P. sylvestris). We set out to investigate three specific speciation scenarios: (I) P. uliginosa is a homoploid hybrid between the other two, (II) the species have evolved without gene flow after divergence and (III) there has been substantial gene flow between the species since their divergence. Overall, the genetic data suggest that P. mugo and P. uliginosa share the same gene pool (average net divergence of 0.0001) and that the phenotypic differences (e.g. growth form) are most likely due to very limited areas of the genome. P. mugo and P. uliginosa are more diverged from P. sylvestris than from each other (average net divergence of 0.0027 and 0.0026, respectively). The nucleotide patterns can best be explained by the divergence with migration speciation scenario, although the hybrid speciation scenario with small genomic contribution from P. sylvestris cannot be completely ruled out. We suggest that the large amount of shared polymorphisms between the pine taxa and the lack of monophyly at all loci studied between P. sylvestris and P. mugo–P. uliginosa can largely be explained by relatively recent speciation history and large effective population sizes but also by interspecific gene flow. These closely related pine taxa form an excellent system for searching for loci involved in adaptive variation as they are differentiated in phenotype and ecology but have very similar genetic background.     

Ecological divergence and speciation between lemur (Eulemur) sister species in Madagascar

Understanding ecological niche evolution over evolutionary timescales is crucial to elucidating the biogeographic history of organisms. Here, we used, for the first time, climate‐based ecological niche models (ENMs) to test hypotheses about ecological divergence and speciation processes between sister species pairs of lemurs (genus Eulemur) in Madagascar. We produced ENMs for eight species, all of which had significant validation support. Among the four sister species pairs, we found nonequivalent niches between sisters, varying degrees of niche overlap in ecological and geographic space, and support for multiple divergence processes. Specifically, three sister‐pair comparisons supported the null model that niches are no more divergent than the available background region. These findings are consistent with an allopatric speciation model, and for two sister pairs (E. collaris–E. cinereiceps and E. rufus–E. rufifrons), a riverine barrier has been previously proposed for driving allopatric speciation. However, for the fourth sister pair E. flavifrons–E. macaco, we found support for significant niche divergence, and consistent with their parapatric distribution on an ecotone and the lack of obvious geographic barriers, these findings most strongly support a parapatric model of speciation. These analyses thus suggest that various speciation processes have led to diversification among closely related Eulemur species.     

Genetic divergence, range expansion and possible homoploid hybrid speciation among pine species in Northeast China

Although homoploid hybrid speciation in plants is probably more common than previously realized, there are few well-documented cases of homoploid hybrid origin in conifers. We examined genetic divergence between two currently widespread pines in Northeast China, Pinus sylvestris var. mongolica and Pinus densiflora, and also whether two narrowly distributed pines in the same region, Pinus funebris and Pinus takahasii, might have originated from the two widespread species by homoploid hybrid speciation. Our results, based on population genetic analysis of chloroplast (cp), mitochondrial (mt) DNA, and nuclear gene sequence variation, showed that the two widespread species were divergent for both cp- and mtDNA variation, and also for haplotype variation at two of eight nuclear gene loci surveyed. Our analysis further indicated that P. sylvestris var. mongolica and P. densiflora remained allopatric during the most severe Quaternary glacial period that occurred in Northeast China, but subsequently exhibited rapid range expansions. P. funebris and P. takahasii, were found to contain a mixture of chlorotypes and nuclear haplotypes that distinguish P. sylvestris var. mongolica and P. densiflora, in support of the hypothesis that they possibly originated via homoploid hybrid speciation following secondary contact and hybridization between P. sylvestris var. mongolica and P. densiflora.     

Demography and speciation history of the homoploid hybrid pine Pinus densata on the Tibetan Plateau

Pinus densata is an ecologically successful homoploid hybrid that inhabits vast areas of heterogeneous terrain on the south‐eastern Tibetan Plateau as a result of multiple waves of colonization. Its region of origin, route of colonization onto the plateau and the directions of introgression with its parental species have previously been defined, but little is known about the isolation and divergence history of its populations. In this study, we surveyed nucleotide polymorphism over eight nuclear loci in 19 representative populations of P. densata and its parental species. Using this information and coalescence simulations, we assessed the historical changes in its population size, gene flow and divergence in time and space. The results indicate a late Miocene origin for P. densata associated with the recent uplift of south‐eastern Tibet. The subsequent differentiation between geographical regions of this species began in the late Pliocene and was induced by regional topographical changes and Pleistocene glaciations. The ancestral P. densata population had a large effective population size but the central and western populations were established by limited founders, suggesting that there were severe bottlenecks during the westward migration out of the ancestral hybrid zone. After separating from their ancestral populations, population expansion occurred in all geographical regions especially in the western range. Gene flow in P. densata was restricted to geographically neighbouring populations, resulting in significant differentiation between regional groups. The new information on the divergence and demographic history of P. densata reported herein enhances our understanding of its speciation process on the Tibetan Plateau.     

Testing the role of climate in speciation: New methods and applications to squamate reptiles (lizards and snakes)

Climate may play important roles in speciation, such as causing the range fragmentation that underlies allopatric speciation (through niche conservatism) or driving divergence of parapatric populations along climatic gradients (through niche divergence). Here, we developed new methods to test the frequency of climate niche conservatism and divergence in speciation, and applied it to species pairs of squamate reptiles (lizards and snakes). We used a large‐scale phylogeny to identify 242 sister species pairs for analysis. From these, we selected all terrestrial allopatric pairs with sufficient occurrence records (n = 49 pairs) and inferred whether each originated via climatic niche conservatism or climatic niche divergence. Among the 242 pairs, allopatric pairs were most common (41.3%), rather than parapatric (19.4%), partially sympatric (17.7%), or fully sympatric species pairs (21.5%). Among the 49 selected allopatric pairs, most appeared to have originated via climatic niche divergence (61–76%, depending on the details of the methods). Surprisingly, we found greater climatic niche divergence between allopatric sister species than between parapatric pairs, even after correcting for geographic distance. We also found that niche divergence did not increase with time, further implicating niche divergence in driving lineage splitting. Overall, our results suggest that climatic niche divergence may often play an important role in allopatric speciation, and the methodology developed here can be used to address the generality of these findings in other organisms.     

The genetic architecture of ecological speciation and the association with signatures of selection in natural lake whitefish (Coregonus sp. Salmonidae) species pairs

Adaptive evolutionary change is contingent on variation and selection; thus, understanding adaptive divergence and ultimately speciation requires information on both the genetic basis of adaptive traits as well as an understanding of the role of divergent natural selection on those traits. The lake whitefish (Coregonus clupeaformis) consists of several sympatric "dwarf" (limnetic) and normal (benthic) species pairs that co-inhabit northern postglacial lakes. These young species pairs have evolved independently and display parallelism in life history, behavioral, and morphological divergence associated with the use of distinct trophic resources. We identified phenotype-environment associations and determined the genetic architecture and the role of selection modulating population genetic divergence in sympatric dwarf and normal lake whitefish. The genetic architecture of 9 adaptive traits was analyzed in 2 hybrid backcrosses individually phenotyped throughout their life history. Significant quantitative trait loci (QTL) were associated with swimming behavior (habitat selection and predator avoidance), growth rate, morphology (condition factor and gill rakers), and life history (onset of maturity and fecundity). Genome scans among 4 natural sympatric pairs, using loci segregating in the map, revealed a signature of selection for 24 loci. Loci exhibiting a signature of selection were associated with QTL relative to other regions of the genome more often than expected by chance alone. Two parallel QTL outliers for growth and condition factor exhibited segregation distortion in both mapping families, supporting the hypothesis that adaptive divergence contributing to parallel reductions of gene flow among natural populations may cause genetic incompatibilities. Overall, these findings offer evidence that the genetic architecture of ecological speciation is associated with signatures of selection in nature, providing strong support for the hypothesis that divergent natural selection is currently maintaining adaptive differentiation and promoting ecological speciation in lake whitefish species pairs.     

Differential barrier strength and allele frequencies in hybrid zones maintained by sex-biased hybrid incompatibilities

In animals, hybrid sterility and inviability between closely related species often affect only the heterogametic sex (XY). This widespread phenomenon, known as Haldane's rule, is an early speciation event found across broad taxa, but the role of heterogametic hybrid incompatibilities, as opposed to homogametic ones, as a barrier in a speciation process remains obscure. It has been hypothesized that heterogametic incompatibility may be a more efficient mechanism in driving speciation. The population dynamics after (rather than before) the occurrence of sex-biased incompatibilities may account for Haldane's rule. In this study, a recursion model of hybrid zones was developed to investigate the differences between heterogametic and homogametic incompatibilities. The selection strengths and selection patterns of sex chromosome-linked, two-locus Bateson-Dobzhansky-Muller (BDM) genetic incompatibilities were examined. It is noted that a sex-biased hybrid incompatibility in a hybrid zone confers asymmetric and uneven impedance to gene flow. The clines of different loci in such a hybrid zone displayed diverse differentiation in their width, steepness and asymmetry. Alleles involved in the incompatibility face much stronger resistance to cross a hybrid zone. Different sex-biased BDM incompatibilities also affect the flow of neutral alleles differently. Compared to a homogametic one, heterogametic incompatibility is a weaker but more asymmetric barrier. These unique patterns of gene flow may explain uneven divergence among different genomic regions during speciation between some closely related species.     

Isolation and gene flow: inferring the speciation history of European house mice

Inferring the history of isolation and gene flow during species differentiation can inform us on the processes underlying their formation. Following their recent expansion in Europe, two subspecies of the house mouse (Mus musculus domesticus and Mus musculus musculus) have formed a hybrid zone maintained by hybrid incompatibilities and possibly behavioural reinforcement, offering a good model of incipient speciation. We reconstruct the history of their divergence using an approximate Bayesian computation framework and sequence variation at 57 autosomal loci. We find support for a long isolation period preceding the advent of gene flow around 200,000 generations ago, much before the formation of the European hybrid zone a few thousand years ago. The duration of the allopatric episode appears long enough (74% of divergence time) to explain the accumulation of many post-zygotic incompatibilities expressed in the present hybrid zone. The ancient contact inferred could have played a role in mating behaviour divergence and laid the ground for further reinforcement. We suggest that both subspecies originally colonized the Middle East from the northern Indian subcontinent, domesticus settling on the shores of the Persian Gulf and musculus on those of the Caspian Sea. Range expansions during interglacials would have induced secondary contacts, presumably in Iran, where they must have also interacted with Mus musculus castaneus. Future studies should incorporate this possibility, and we point to Iran and its surroundings as a hot spot for house mouse diversity and speciation studies.     

Metabolic Rate and Climatic Fluctuations Shape Continental Wide Pattern of Genetic Divergence and Biodiversity in Fishes

Taxonomically exhaustive and continent wide patterns of genetic divergence within and between species have rarely been described and the underlying evolutionary causes shaping biodiversity distribution remain contentious. Here, we show that geographic patterns of intraspecific and interspecific genetic divergence among nearly all of the North American freshwater fish species (>750 species) support a dual role involving both the late Pliocene-Pleistocene climatic fluctuations and metabolic rate in determining latitudinal gradients of genetic divergence and very likely influencing speciation rates. Results indicate that the recurrent glacial cycles caused global reduction in intraspecific diversity, interspecific genetic divergence, and species richness at higher latitudes. At the opposite, longer geographic isolation, higher metabolic rate increasing substitution rate and possibly the rapid accumulation of genetic incompatibilities, led to an increasing biodiversity towards lower latitudes. This indicates that both intrinsic and extrinsic factors similarly affect micro and macro evolutionary processes shaping global patterns of biodiversity distribution. These results also indicate that factors favouring allopatric speciation are the main drivers underlying the diversification of North American freshwater fishes.     

The influence of life‐history strategy on genetic differentiation and lineage divergence in darters (Percidae: Etheostomatinae)

Recent studies determined that darters with specialized breeding strategies can exhibit deep lineage divergence over fine geographic scales without apparent physical barriers to gene flow. However, the extent to which intrinsic characteristics interact with extrinsic factors to influence population divergence and lineage diversification in darters is not well understood. This study employed comparative phylogeographic and population genetic methods to investigate the influence of life history on gene flow, dispersal ability, and lineage divergence in two sympatric sister darters with differing breeding strategies. Our results revealed highly disparate phylogeographic histories, patterns of genetic structure, and dispersal abilities between the two species suggesting that life history may contribute to lineage diversification in darters, especially by limiting dispersal among large river courses. Both species also showed striking differences in demographic history, indicating that extrinsic factors differentially affected each species during the Pleistocene. Collectively, our results indicate that intrinsic and extrinsic factors have influenced levels of gene flow among populations within both species examined. However, we suggest that life‐history strategy may play a more important role in lineage diversification in darters than previously appreciated, a finding that has potentially important implications for understanding diversification of the rich North American freshwater fish fauna.