Competition and evolution along environmental gradients: patterns, boundaries and sympatric divergence

The present study examined how competitive interactions and environmental conditions generate species boundaries and determine species distributions. A spatially explicit, quantitative genetic, two-species competition model was used to manipulate the strengths of competition, gene flow and local adaptation along environmental gradients. This allowed us to assess the long-term persistence of each species and whether the ranges they inhabited had boundaries in space or were unlimited. We found that a species boundary arises along less steep environmental gradients when the strength of stabilizing selection and diversifying selection are similar. We also found that a species boundary may arise along shallow environmental gradients if interspecific competition is more intense than intraspecific, which relaxes previous requirements for steep gradients for generating range limits. We determined an analytical form for the critical environmental gradient as a function of ecological and genetic parameters at which a species boundary is expected to arise by competition. Results suggest an alternative to resource competition as an explanation for phenotypic divergence between sympatric competitors. Competitors sharing a trait that is under stabilizing selection along an environmental gradient may segregate spatially and evolve in different regions, with phenotypic sympatric divergence reflecting the resulting clines. Along various types of environmental gradients, variation in stabilizing selection intensities could lead to contrasting patterns in the distribution of species. For stabilizing selection strengths in accord with field data estimates, this study predicts that the level of sympatric character divergence would be limited along environmental gradients.     

An intraspecific comparative analysis of character divergence between sympatric species

Although sympatric character divergence between closely related species has been described in a wide variety of taxa, the evolutionary processes responsible for generating these patterns are difficult to identify. One hypothesis that can explain sympatric differences is ecological character displacement: the sympatric origin of morphologically divergent phenotypes in response to selection caused by interspecific competition. Alternatively, populations may adapt to different conditions in allopatry, with sympatric distributions evolving through selective colonization and proliferation of ecologically compatible phenotypes. In this study, I characterize geographic variation within two sibling species of rocky-shore gastropods that have partially overlapping distributions in central California. In sympatry, both Nucella emarginata and N. ostrina show significant differences in shell shape and shell ornamentation that together suggest that where the two species co-exist, divergent phenotypes arose as an evolutionary consequence of competition. To examine the evolutionary origins of divergent characters in sympatry, I used a comparative method based on spatial autocorrelation to remove the portion of the phenotypic variance among populations that is explained by genetic distance (using mitochondrial DNA sequences and allozyme frequency data). Because the remaining portion of the phenotypic variance represents the independent divergence of individual populations, a significant sympatric difference in the corrected dataset provides evidence of true character displacement: significant sympatric character evolution that is independent of population history. After removal of genetic distance effects in Nucella, shell shape differences remain statistically significant in N. emarginata, providing evidence of significant sympatric character divergence. However, for external shell ornamentation in both species and shell shape in N. ostrina, the significance of sympatric differences is lost in the corrected dataset, indicating that colonization events and gene flow have played important roles in the evolutionary history of character divergence in sympatry. Although the absence of a widely dispersing planktonic larva in the life cycle of Nucella will promote local adaptation, the results here indicate that once advantageous traits arise, demographic processes, such as recurrent gene flow between established populations and extinction and recolonization, are important factors contributing to the geographic pattern of sympatric character divergence.     

Location-specific sympatric morphological divergence as a possible response to species interactions in West Virginia Plethodon salamander communities

1. The competitive interactions of closely related species have long been considered important determinants of community composition and a major cause of phenotypic diversification. However, while patterns such as character displacement are well documented, less is known about how local adaptation influences diversifying selection from interspecific competition. 2. We examined body size and head shape variation among allopatric and sympatric populations of two salamander species, the widespread Plethodon cinereus and the geographically restricted P. nettingi. We quantified morphology from 724 individuals from 20 geographical localities throughout the range of P. nettingi. 3. Plethodon nettingi was more robust in cranial morphology relative to P. cinereus, and sympatric localities were more robust relative to allopatric localities. Additionally, there was significantly greater sympatric head shape divergence between species relative to allopatric communities, and sympatric localities of P. cinereus exhibited greater morphological variation than sympatric P. nettingi. 4. The sympatric morphological divergence and increase in cranial robustness of one species (P. nettingi) were similar to observations in other Plethodon communities, and were consistent with the hypothesis of interspecific competition. These findings suggest that interspecific competition in Plethodon may play an important role in phenotypic diversification in this group. 5. The increase in among-population variance in sympatric P. cinereus suggests a species-specific response to divergent natural selection that is influenced in part by other factors. We hypothesize that enhanced morphological flexibility and ecological tolerance allow P. cinereus to more rapidly adapt to local environmental conditions, and initial differences among populations have allowed the evolutionary response of P. cinereus to vary across replicate sympatric locations, resulting in distinct evolutionary trajectories of morphological change.     

Parallel evolution of character displacement driven by competitive selection in terrestrial salamanders

Background  

Parallel evolution can occur when common environmental factors exert similar selective forces on morphological variation in populations in different geographic localities. Competition can also generate morphological shifts, and if competing species co-occur in multiple geographic regions, then repeated instances of competitively-driven morphological divergence (character displacement) can occur. Despite the importance of character displacement for inferring the role of selection in morphological evolution however, replicated instances of sympatric morphological divergence are understudied.     

The role of selection and historical factors in driving population differentiation along an elevational gradient in an island bird

Adaptation to local environmental conditions and the range dynamics of populations can influence evolutionary divergence along environmental gradients. Thus, it is important to investigate patterns of both phenotypic and genetic variations among populations to reveal the respective roles of these two types of factors in driving population differentiation. Here, we test for evidence of phenotypic and genetic structure across populations of a passerine bird (Zosterops borbonicus) distributed along a steep elevational gradient on the island of Réunion. Using 11 microsatellite loci screened in 401 individuals from 18 localities distributed along the gradient, we found that genetic differentiation occurred at two spatial levels: (i) between two main population groups corresponding to highland and lowland areas, respectively, and (ii) within each of these two groups. In contrast, several morphological traits varied gradually along the gradient. Comparison of neutral genetic differentiation (F_ST) and phenotypic differentiation (P_ST) showed that P_ST largely exceeds F_ST at several morphological traits, which is consistent with a role for local adaptation in driving morphological divergence along the gradient. Overall, our results revealed an area of secondary contact midway up the gradient between two major, cryptic, population groups likely diverged in allopatry. Remarkably, local adaptation has shaped phenotypic differentiation irrespective of population history, resulting in different patterns of variation along the elevational gradient. Our findings underscore the importance of understanding both historical and selective factors when trying to explain variation along environmental gradients.     

Character shifts in the defensive armor of sympatric sticklebacks

Natural enemies may contribute to the morphological divergence of sympatric species, yet their role has received little attention to date. We tested for character shifts in defensive armor of sympatric threespine sticklebacks (Gasterosteus aculeatus complex) previously shown to exhibit ecological character displacement in traits related to resource use. We scored five defensive armor traits in sympatric benthic and limnetic stickleback species from southwestern British Columbia and compared them with the same traits in nearby allopatric populations in the presence of the same predatory fish (Oncorhynchus sp.). This approach is analogous to tests of ecological character displacement that compare trophic traits of sympatric and allopatric species in the presence of the same community of resource types. Three patterns consistent with character displacement in defensive armor were found. First, limnetics in different lakes had consistently more armor than sympatric benthics. Second, the average amount of armor, averaged over both species, was reduced in sympatry compared to allopatric populations. This reduction was almost entirely the result of shifts by benthic species, whereas armor in limnetics was more similar to that in allopatric populations. Third, differences between sympatric benthics and limnetics in total armor were greater than expected from comparisons with allopatric populations. We interpret these patterns as the result of differences in habitat-specific predation regimes accompanying ecological character displacement and indirect interactions between sympatric stickleback species mediated by their top predators. These results suggest that predation may facilitate, rather than hinder, the process of divergence in sympatry.     

Do different rates of gene flow underlie variation in phenotypic and phenological clines in a montane grasshopper community?

Species responses to environmental change are likely to depend on existing genetic and phenotypic variation, as well as evolutionary potential. A key challenge is to determine whether gene flow might facilitate or impede genomic divergence among populations responding to environmental change, and if emergent phenotypic variation is dependent on gene flow rates. A general expectation is that patterns of genetic differentiation in a set of codistributed species reflect differences in dispersal ability. In less dispersive species, we predict greater genetic divergence and reduced gene flow. This could lead to covariation in life‐history traits due to local adaptation, although plasticity or drift could mirror these patterns. We compare genome‐wide patterns of genetic structure in four phenotypically variable grasshopper species along a steep elevation gradient near Boulder, Colorado, and test the hypothesis that genomic differentiation is greater in short‐winged grasshopper species, and statistically associated with variation in growth, reproductive, and physiological traits along this gradient. In addition, we estimate rates of gene flow under competing demographic models, as well as potential gene flow through surveys of phenological overlap among populations within a species. All species exhibit genetic structure along the elevation gradient and limited gene flow. The most pronounced genetic divergence appears in short‐winged (less dispersive) species, which also exhibit less phenological overlap among populations. A high‐elevation population of the most widespread species, Melanoplus sanguinipes, appears to be a sink population derived from low elevation populations. While dispersal ability has a clear connection to the genetic structure in different species, genetic distance does not predict growth, reproductive, or physiological trait variation in any species, requiring further investigation to clearly link phenotypic divergence to local adaptation.     

Differentiation in phenotypic plasticity among populations of Arabis serrata Fr. & Sav. (Brassicaceae)

Studies on divergence of phenotypic plasticity in closely related species have suggested that character means and plasticity of these characters may evolve independently. Similar patterns of divergence between populations within a species have been reported although few plant species have been studied. Thus, in this paper, the patterns of differentiation between character means and phenotypic plasticity among eight populations of Arabis serrata are documented. Mean response and magnitude and pattern of phenotypic plasticity were measured and compared in plants growing under an environmental gradient of nutrients. Differences in means and coefficients of variation (CV as indicators of plasticity) among populations were compared using the Canberra metric and generating unrooted Wagner trees. Populations showed significant differences in character means in nine morphological traits. Magnitude and patterns of phenotypic plasticity showed a complex pattern of differentiation for each trait and population. Biomass traits were more plastic, in general, than characters associated with linear size. Comparisons between pairs of populations for nine morphological traits showed that in 28.6% of 252 possible cases, populations differed in means, magnitude and patterns of phenotypic plasticity. In almost 90% of the cases, populations differed in magnitude and/or pattern of plasticity. Considering all characters together, populations from similar habitats and with common life history features tended to respond in similar ways. The patterns of divergence, however, suggest that character means and character plasticities among populations are able to evolve independently.     

Reproductive character displacement and environmental filtering shape floral variation between sympatric sister taxa

Divergence in reproductive traits between closely related species that co‐occur contributes to speciation by reducing interspecific gene flow. In flowering plants, greater floral divergence in sympatry than allopatry may reflect reproductive character displacement (RCD) by means of divergent pollinator‐mediated selection or mating system evolution. However, environmental filtering (EF) would prevail for floral traits under stronger selection by abiotic factors than pollination, and lead to sympatric taxa being more phenotypically similar. We determine whether floral UV pigmentation and size show signatures of RCD or EF using a biogeographically informed sister taxa comparison. We determine whether 35 sister pairs in the Potentilleae tribe (Rosaceae) are allopatric or sympatric and confirm that sympatric sisters experience more similar bioclimatic conditions, an assumption of the EF hypothesis. We test whether interspecific differences are greater in allopatry or sympatry while accounting for divergence time. For UV pigmentation, sympatric sisters are more phenotypically similar than allopatric ones. For flower size, sympatric sisters show increased divergence with time since speciation but allopatric ones do not. We conclude that floral UV pigmentation shows a signature of EF, whereas flower size shows a signature of RCD. Discordant results between the traits suggest that the dominant selective agent differs between them.     

Moose body mass variation revisited: disentangling effects of environmental conditions and genetics

Large-scale geographical variation in phenotypic traits within species is often correlated to local environmental conditions and population density. Such phenotypic variation has recently been shown to also be influenced by genetic structuring of populations. In ungulates, large-scale geographical variation in phenotypic traits, such as body mass, has been related to environmental conditions and population density, but little is known about the genetic influences. Research on the genetic structure of moose suggests two distinct genetic lineages in Norway, structured along a north-south gradient. This corresponds with many environmental gradients, thus genetic structuring provides an additional factor affecting geographical phenotypic variation in Norwegian moose. We investigated if genetic structure explained geographical variation in body mass in Norwegian moose while accounting for environmental conditions, age and sex, and if it captured some of the variance in body mass that previously was attributed to environmental factors. Genetic structuring of moose was the most important variable in explaining the geographic variation in body mass within age and sex classes. Several environmental variables also had strong explanatory power, related to habitat diversity, environmental seasonality and winter harshness. The results suggest that environmental conditions, landscape characteristics, and genetic structure should be evaluated together when explaining large-scale patterns in phenotypic characters or life history traits. However, to better understand the role of genetic and environmental effects on phenotypic traits in moose, an extended individual-based study of variation in fitness-related characters is needed, preferably in an area of convergence between different genetic lineages.     

Sympatric divergence and clinal variation in multiple coloration traits of Ficedula flycatchers

Geographic variation in phenotypes plays a key role in fundamental evolutionary processes such as local adaptation, population differentiation and speciation, but the selective forces behind it are rarely known. We found support for the hypothesis that geographic variation in plumage traits of the pied flycatcher Ficedula hypoleuca is explained by character displacement with the collared flycatcher Ficedula albicollis in the contact zone. The plumage traits of the pied flycatcher differed strongly from the more conspicuous collared flycatcher in a sympatric area but increased in conspicuousness with increasing distance to there. Phenotypic differentiation (P_ST) was higher than that in neutral genetic markers (F_ST), and the effect of geographic distance remained when statistically controlling for neutral genetic differentiation. This suggests that a cline created by character displacement and gene flow explains phenotypic variation across the distribution of this species. The different plumage traits of the pied flycatcher are strongly to moderately correlated, indicating that they evolve non‐independently from each other. The flycatchers provide an example of plumage patterns diverging in two species that differ in several aspects of appearance. The divergence in sympatry and convergence in allopatry in these birds provide a possibility to study the evolutionary mechanisms behind the highly divergent avian plumage patterns.     

Post‐glacial colonization routes coincide with a life‐history breakpoint along a latitudinal gradient

Although adaptive divergence along environmental gradients has repeatedly been demonstrated, the role of post‐glacial colonization routes in determining phenotypic variation along gradients has received little attention. Here, we used a hierarchical Q_ST–F_ST approach to separate the roles of adaptive and neutral processes in shaping phenotypic variation in moor frog (Rana arvalis) larval life histories along a 1,700 km latitudinal gradient across northern Europe. This species has colonized Scandinavia via two routes with a contact zone in northern Sweden. By using neutral SNP and common garden phenotypic data from 13 populations at two temperatures, we showed that most of the variation along the gradient occurred between the two colonizing lineages. We found little phenotypic divergence within the lineages; however, all phenotypic traits were strongly diverged between the southern and northern colonization routes, with higher growth and development rates and larger body size in the north. The Q_ST estimates between the colonization routes were four times higher than F_ST, indicating a prominent role for natural selection. Q_ST within the colonization routes did not generally differ from F_ST, but we found temperature‐dependent adaptive divergence close to the contact zone. These results indicate that lineage‐specific variation can account for much of the adaptive divergence along a latitudinal gradient.     

Reproductive character displacement and the genetics of gamete recognition in tropical sea urchins

Reproductive character displacement occurs when sympatric and allopatric populations of a species differ in traits crucial to reproduction, and it is commonly thought of as a signal of selection acting to limit hybridization. Most documented cases of reproductive character displacement involve characters that are poorly understood at the genetic level, and rejecting alternative hypotheses for biogeographic shifts in reproductive traits is often very difficult. In sea urchins, the gamete recognition protein bindin evolves under positive selection when species are broadly sympatric, suggesting character displacement may be operating in this system. We sampled sympatric and allopatric populations of two species in the sea urchin genus Echinometra for variation in bindin and for the mitochondrial cytochrome oxidase I to examine patterns of population differentiation and molecular evolution at a reproductive gene. We found a major shift in bindin alleles between central Pacific (allopatric) and western Pacific (sympatric) populations of E. oblonga. Allopatric populations of E. oblonga are polyphyletic with E. sp. C at bindin, whereas sympatric populations of the two species are reciprocally monophyletic. There is a strong signal of positive selection (P(N)/P(S) = 4.5) in the variable region of the first exon of bindin, which is associated with alleles found in sympatric populations of E. oblonga. These results indicate that there is a strong pattern of reproductive character displacement between E. oblonga and E. sp. C and that the divergence is driven by selection. There is much higher population structure in sympatric populations at the bindin locus than at the neutral mitochondrial locus, but this difference is not seen in allopatric populations. These data suggest a pattern of speciation driven by selection for local gamete coevolution as a result of interactions between sympatric species. Although this pattern is highly suggestive of speciation by reinforcement, further research into hybrid fitness and egg-sperm interactions is required to address this potential mechanism for character displacement.     

Genetic constraints on the evolution of mate recognition under natural selection

Field populations of Drosophila serrata display reproductive character displacement in cuticular hydrocarbons (CHCs) when sympatric with Drosophila birchii. We have previously shown that the naturally occurring pattern of reproductive character displacement can be experimentally replicated by exposing field allopatric populations of D. serrata to experimental sympatry with D. birchii. Here, we tested whether the repeated evolution of reproductive character displacement in natural and experimental populations was a consequence of genetic constraints on the evolution of CHCs. The genetic variance-covariance (G) matrices for CHCs were determined for populations of D. serrata that had evolved in either the presence or absence of D. birchii under field and experimental conditions. Natural selection on mate recognition under both field and experimental sympatric conditions increased the genetic variance in CHCs consistent with a response to selection based on rare alleles. A close association between G eigenstructure and the eigenstructure of the phenotypic divergence (D) matrix in natural and experimental populations suggested that G matrix eigenstructure may have determined the direction in which reproductive character displacement evolved during the reinforcement of mate recognition.     

Experimental evidence that competition between stickleback species favours adaptive character divergence

The ecological character displacement hypothesis assumes that the effects of interspecific resource competition cause divergent selection to favour phenotypes that exploit non-shared resources. This model predicts that interspecific competition declines with increased divergence. Direct tests of this decline are rare despite much comparative evidence for character displacement. We tested this prediction using a pair of divergent brook sticklebacks. Brook sticklebacks sympatric with ninespine sticklebacks have diverged from local allopatric brook populations, and so the two types of brook sticklebacks potentially represent pre- and post-displacement forms. We used enclosures placed in a lake to compare short-term fitness (growth) of sympatric (post-displacement) and allopatric (pre-displacement) brook forms in the presence and absence of ninespine sticklebacks. Brook sticklebacks grew less in the presence vs. absence of ninespine sticklebacks, indicating that interspecific competition occurred. As expected, allopatric brook forms had lower growth than sympatric forms when ninespine sticklebacks were present. This result suggests that ecological character displacement has occurred.     

Frequency dependent natural selection during character displacement in sticklebacks

Abstract We know little about how natural selection on a species is altered when a closely related species consuming similar resources appears in its environment. In a pond experiment with threespine sticklebacks I tested the prediction that divergent natural selection between competitors is frequency-dependent, changing with the distribution of phe-notypes in the environment. Differential growth and survival of phenotypes in a target stickleback population were contrasted between two treatments. In one treatment an offshore zooplankton feeder (the limnetic stickleback species) was added to the same pond as the target. In the other treatment I added the benthic stickleback instead, a species adapted to feeding on invertebrates from sediments and inshore vegetation. The target population was ecologically and morphologically intermediate with phenotypic variance artificially inflated by hybridization. Growth rates of phenotypes within the target population differed between treatments as predicted by character displacement. The impact of adding a second species always fell most heavily on those phenotypes in the target population resembling the added species most closely. However, those individuals in the target population that most resembled the added species did not experience reduced survival. Instead, consistent survival differences between populations suggested the presence of an inshore –offshore gradient in mortality risk. These results provide further support for the hypothesis of character displacement in sympatric sticklebacks. They suggest that displacement along the resource gradient also led to divergence in vulnerability to agents of mortality, probably including predation.     

Contemporary and historical evolutionary processes interact to shape patterns of within-lake phenotypic divergences in polyphenic pumpkinseed sunfish, Lepomis gibbosus

Historical and contemporary evolutionary processes can both contribute to patterns of phenotypic variation among populations of a species. Recent studies are revealing how interactions between historical and contemporary processes better explain observed patterns of phenotypic divergence than either process alone. Here, we investigate the roles of evolutionary history and adaptation to current environmental conditions in structuring phenotypic variation among polyphenic populations of sunfish inhabiting 12 postglacial lakes in eastern North America. The pumpkinseed sunfish polyphenism includes sympatric ecomorphs specialized for littoral or pelagic lake habitats. First, we use population genetic methods to test the evolutionary independence of within-lake phenotypic divergences of ecomorphs and to describe patterns of genetic structure among lake populations that clustered into three geographical groupings. We then used multivariate analysis of covariance (MANCOVA) to partition body shape variation (quantified with geometric morphometrics) among the effects of evolutionary history (reflecting phenotypic variation among genetic clusters), the shared phenotypic response of all populations to alternate habitats within lakes (reflecting adaptation to contemporary conditions), and unique phenotypic responses to habitats within lakes nested within genetic clusters. All effects had a significant influence on body form, but the effects of history and the interaction between history and contemporary habitat were larger than contemporary processes in structuring phenotypic variation. This highlights how divergence can be better understood against a known backdrop of evolutionary history.     

Adaptive divergence in moor frog (Rana arvalis) populations along an acidification gradient: inferences from Q(st) -F(st) correlations

Microevolutionary responses to spatial variation in the environment seem ubiquitous, but the relative role of selection and neutral processes in driving phenotypic diversification remain often unknown. The moor frog (Rana arvalis) shows strong phenotypic divergence along an acidification gradient in Sweden. We here used correlations among population pairwise estimates of quantitative trait (P(ST) or Q(ST) from common garden estimates of embryonic acid tolerance and larval life-history traits) and neutral genetic divergence (F(ST) from neutral microsatellite markers), as well as environmental differences (pond pH, predator density, and latitude), to test whether this phenotypic divergence is more likely due to divergent selection or neutral processes. We found that trait divergence was more strongly correlated with environmental differences than the neutral marker divergence, suggesting that divergent natural selection has driven phenotypic divergence along the acidification gradient. Moreover, pairwise P(ST) s of embryonic acid tolerance and Q(ST) s of metamorphic size were strongly correlated with breeding pond pH, whereas pairwise Q(ST) s of larval period and growth rate were more strongly correlated with geographic distance/latitude and predator density, respectively. We suggest that incorporating measurements of environmental variation into Q(ST) -F(ST) studies can improve our inferential power about the agents of natural selection in natural populations.     

Adaptation to local ultraviolet radiation conditions among neighbouring Daphnia populations

Understanding the historical processes that generated current patterns of phenotypic diversity in nature is particularly challenging in subdivided populations. Populations often exhibit heritable genetic differences that correlate with environmental variables, but the non-independence among neighbouring populations complicates statistical inference of adaptation. To understand the relative influence of adaptive and non-adaptive processes in generating phenotypes requires joint evaluation of genetic and phenotypic divergence in an integrated and statistically appropriate analysis. We investigated phenotypic divergence, population-genetic structure and potential fitness trade-offs in populations of Daphnia melanica inhabiting neighbouring subalpine ponds of widely differing transparency to ultraviolet radiation (UVR). Using a combination of experimental, population-genetic and statistical techniques, we separated the effects of shared population ancestry and environmental variables in predicting phenotypic divergence among populations. We found that native water transparency significantly predicted divergence in phenotypes among populations even after accounting for significant population structure. This result demonstrates that environmental factors such as UVR can at least partially account for phenotypic divergence. However, a lack of evidence for a hypothesized trade-off between UVR tolerance and growth rates in the absence of UVR prevents us from ruling out the possibility that non-adaptive processes are partially responsible for phenotypic differentiation in this system.     

Speciation in ninespine stickleback: reproductive isolation and phenotypic divergence among cryptic species of Japanese ninespine stickleback

Although similar patterns of phenotypic diversification are often observed in phylogenetically independent lineages, differences in the magnitude and direction of phenotypic divergence have been also observed among independent lineages, even when exposed to the same ecological gradients. The stickleback family is a good model with which to explore the ecological and genetic basis of parallel and nonparallel patterns of phenotypic evolution, because there are a variety of populations and species that are locally adapted to divergent environments. Although the patterns of phenotypic divergence as well as the genetic and ecological mechanisms have been well characterized in threespine sticklebacks, Gasterosteus aculeatus, we know little about the patterns of phenotypic diversification in other stickleback lineages. In eastern Hokkaido, Japan, there are three species of ninespine sticklebacks, Pungitius tymensis and the freshwater type and the brackish‐water type of the P. pungitius–P. sinensis species complex. They utilize divergent habitats along coast–stream gradients of rivers. Here, we investigated genetic, ecological and phenotypic divergence among three species of Japanese ninespine sticklebacks. Divergence in trophic morphology and salinity tolerance occurred in the direction predicted by the patterns observed in threespine sticklebacks. However, the patterns of divergence in armour plate were different from those previously found in threespine sticklebacks. Furthermore, the genetic basis of plate variation may differ from that in threespine sticklebacks. Because threespine sticklebacks are well‐established model for evolutionary research, the sympatric trio of ninespine sticklebacks will be an invaluable resource for ecological and genetic studies on both common and lineage‐specific patterns of phenotypic diversification.