Niche divergence of two closely related Carbula species (Insecta: Hemiptera: Pentatomidae) despite the presence of a hybrid zone

1. The extent to which ecologically divergent selection acts to maintain species boundaries in the presence of hybridisation and gene flow is not well understood. Two parapatric taxa, Carbula humerigera (Uhler) and Carbula putoni (Jakovlev), were used to test the extent to which niche differentiation might sustain divergence of related taxa despite ongoing gene flow. 2. Mitochondrial [cytochrome c oxidase subunit I (COI), cytochrome c oxidase subunit II (COII) and cytochrome b (Cytb)] and nuclear [elongation factor 1‐alpha (EF‐1α)] markers were sequenced from 383 individuals. Clusters of morphological variations were estimated and visualised using principal component analysis (PCA). Haplotype networks were constructed using median‐joining and neighbour‐net algorithm methods. Gene flows were estimated using migrate‐n analyses. Suitable habitats for each species and their sympatric distribution were predicted with ecological niche modelling (ENM). Niche comparisons were conducted using Schoener's D and Warren's I. In addition, PCA, multivariate analysis of variance (manova ) and discriminant function analysis were used to test ecological differentiation. 3. Morphological clusters and network analysis indicated that samples were generally divided into three groups (C. humerigera, C. putoni and hybrids). Ongoing gene flow was detected among the three groups, with abundant magnitudes in the sympatric region. Niche comparison and statistical analysis showed ecological differentiation between C. humerigera and C. putoni. Two potential hybrid zones were predicted in the ENM reconstruction, located along the Yanshan‐Taihang‐Qinling and Taishan mountains. 4. These results reveal a geographically delineated hybrid zone between C. humerigera and C. putoni. These two closely related Carbula species still live in different ecological niches despite hybridisation and ongoing gene flow.     

The genomics of adaptation,divergence and speciation: a congealing theory

In this issue, Flaxman et al. ( 2014 ) report the results of sophisticated whole‐genome simulations of speciation with gene flow, enhancing our understanding of the process by building on previous single‐locus, multilocus and analytical works. Their findings provide us with new insights about how genomes can diverge and the importance of statistical and chromosomal linkage in facilitating reproductive isolation. The authors characterize the conditions under which, even with high gene flow and weak divergent selection, reproductive isolation between populations can occur due to the emergent stochastic process of genomewide congealing, where numerous statistically or physically linked loci of small effect allow selection to limit effective migration rates. The initial congealing event can occur within a broad range conditions, and once initiated, the self‐reinforcing process leads to rapid divergence and ultimately two reproductively isolated populations. Flaxman et al.'s ( 2014 ) work is a valuable contribution to our understanding of speciation with gene flow and in making a more predictive field of evolutionary genomics and speciation.     

Ecological divergence in the presence of gene flow in two closely related Oryza species (Oryza rufipogon and O. nivara)

Ecological divergence plays a prominent role in the process of speciation, but how divergence occurs in the face of gene flow is still less clear, and remains controversial among evolutionists. Here we investigated the nucleotide diversity, divergence and gene flow between Oryza nivara and O. rufipogon using sequences of seven chloroplast and nuclear loci. By analysing samples from 26 wild populations across the geographic ranges of the two species, we showed that both species were highly structured and O. rufipogon maintained a higher level of species‐wide diversity than O. nivara. Notably, phylogenetic, amova and F_ST analyses were unable to detect significant nucleotide differentiation between the two species. We estimated that the two species began to diverge at c. 0.16 million years ago. Our coalescent‐based simulations strongly rejected the simple isolation model of zero migration between species, but rather provided unambiguous evidence of bidirectional gene flow between species, particularly from O. rufipogon to O. nivara. Our simulations also indicated that gene flow was recurrent during the divergence process rather than arising from secondary contact after allopatric divergence. In conjunction with different morphological and life‐history traits and habitat preference in the two species, this study supports the hypothesis that these Oryza species are better treated as ecotypes that diverged quite recently and are still under the process of divergence. Importantly, we demonstrate the ecological divergence between O. rufipogon and O. nivara in the presence of significant gene flow, implying that natural selection plays a primary role in driving the divergence of the two Oryza species.     

Habitat-associated morphological divergence in two Neotropical fish species

We examined intraspecific morphological diversification between river channel and lagoon habitats for two Neotropical fish ( Bryconops caudomaculatus , Characidae; Biotodoma wavrini , Cichlidae). We hypothesized that differences between habitats (e.g. flow regime, foraging opportunities) might create selective pressures resulting in morphological divergence between conspecific populations. We collected fish from four channel-lagoon habitat pairs in the Río Cinaruco, Venezuela, and compared body morphology using geometric morphometrics. There were two aspects of divergence in both species: (1) placement of maximum body depth and (2) orientation of the mouth. For both species, maximum body depth was positioned more anteriorly (i.e. fusiform) in the river channel than in lagoons. Both species exhibited a relatively terminal mouth in lagoons compared to the channel. The mouth of B. caudomaculatus was relatively upturned, whereas the mouth of B. wavrini was relatively subterminal, in channel habitats. Observed morphological patterns are consistent with functional morphological principles suggesting adaptive divergence. We also show that spatial distance between habitats, presumably reflecting rates of population mixing, appears to have constrained diversification. For both species, morphological divergence increased with distance between habitats. Thus morphological divergence between channel and lagoon habitats apparently reflects a balance between diversification driven by natural selection, and homogenization driven by population mixing.  © 2003 The Linnean Society of London, Biological Journal of the Linnean Society , 2003, 80 , 689–698.     

Influence of voltine ecotype and geographic distance on genetic and haplotype variation in the Asian corn borer

Diapause is an adaptive dormancy strategy by which arthropods endure extended periods of adverse climatic conditions. Seasonal variation in larval diapause initiation and duration in Ostrinia furnacalis may influence adult mating generation number (voltinism) across different local environments. The degree to which voltine ecotype, geographic distance, or other ecological factors influence O. furnacalis population genetic structure remains uncertain. Genetic differentiation was estimated between voltine ecotypes collected from 8 locations. Mitochondrial haplotypes were significantly different between historically allopatric univoltine and bivoltine locations, but confounded by a strong correlation with geographic distance. In contrast, single nucleotide polymorphism (SNP) genotypes show low but significant levels of variation and a lack of influence of geographic distance between allopatric voltine locations. Regardless, 11 of 257 SNP loci were predicted to be under selection, suggesting population genetic homogenization except at loci proximal to factors putatively under selection. These findings provide evidence of haplotype divergent voltine ecotypes that may be maintained in allopatric and sympatric areas despite relatively high rates of nuclear gene flow, yet influence of voltinism on maintenance of observed haplotype divergence remains unresolved.     

Population differentiation and restricted gene flow in Spanish crossbills: not isolation-by-distance but isolation-by-ecology

Divergent selection stemming from environmental variation may induce local adaptation and ecological speciation whereas gene flow might have a homogenizing effect. Gene flow among populations using different environments can be reduced by geographical distance (isolation-by-distance) or by divergent selection stemming from resource use (isolation-by-ecology). We tested for and encountered phenotypic and genetic divergence among Spanish crossbills utilizing different species of co-occurring pine trees as their food resource. Morphological, vocal and mtDNA divergence were not correlated with geographical distance, but they were correlated with differences in resource use. Resource diversity has now been found to repeatedly predict crossbill diversity. However, when resource use is not 100% differentiated, additional characters (morphological, vocal, genetic) must be used to uncover and validate hidden population structure. In general, this confirms that ecology drives adaptive divergence and limits neutral gene flow as the first steps towards ecological speciation, unprevented by a high potential for gene flow.     

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.     

Landscape drivers of genomic diversity and divergence in woodland Eucalyptus

Spatial genetic patterns are influenced by numerous factors, and they can vary even among coexisting, closely related species due to differences in dispersal and selection. Eucalyptus (L'Héritier 1789; the “eucalypts”) are foundation tree species that provide essential habitat and modulate ecosystem services throughout Australia. Here we present a study of landscape genomic variation in two woodland eucalypt species, using whole‐genome sequencing of 388 individuals of Eucalyptus albens and Eucalyptus sideroxylon. We found exceptionally high genetic diversity (π ≈ 0.05) and low genome‐wide, interspecific differentiation (F_ST = 0.15) and intraspecific differentiation between localities (F_ST ≈ 0.01–0.02). We found no support for strong, discrete population structure, but found substantial support for isolation by geographic distance (IBD) in both species. Using generalized dissimilarity modelling, we identified additional isolation by environment (IBE). Eucalyptus albens showed moderate IBD, and environmental variables have a small but significant amount of additional predictive power (i.e. IBE). Eucalyptus sideroxylon showed much stronger IBD and moderate IBE. These results highlight the vast adaptive potential of these species and set the stage for testing evolutionary hypotheses of interspecific adaptive differentiation across environments.     

Comparing geographical genetic differentiation between candidate and noncandidate loci for adaptation strengthens support for parallel ecological divergence in the marine snail Littorina saxatilis

The Galician sympatric ecotypes of Littorina saxatilis have been proposed as a model system for studying parallel ecological speciation. Such a model system makes a clear prediction: candidate loci (for divergent adaptation) should present a higher level of geographical differentiation than noncandidate (neutral) loci. We used 2356 amplified fragment length polymorphisms (AFLPs) and four microsatellite loci to identify candidate loci for ecological adaptation using the F _ST outlier method. Three per cent of the studied AFLP loci were identified as candidate loci associated with adaptation, after multitest adjustments, thus contributing to ecotype differentiation (candidate loci were not detected within ecotypes). Candidate and noncandidate loci were analysed separately at four different F _ST partitions: differences between ecotypes (overall and local), differences between localities and micro-geographical differences within ecotypes. The magnitude of F _ST differed between candidate and noncandidate loci for all partitions except in the case of microgeographical differentiation within ecotypes, and the microsatellites (putatively neutral) showed an identical pattern to noncandidate loci. Thus, variation in candidate loci is determined partially independent by divergent natural selection (in addition to stochastic forces) at each locality, while noncandidate loci are exclusively driven by stochastic forces. These results support the evolutionary history described for these particular populations, considered to be a clear example of incomplete sympatric ecological speciation.     

Strong differences in chemical recognition cues between two closely related species of ants from the genus Lasius (Hymenoptera: Formicidae)

Cuticular hydrocarbons (CHCs) are increasingly recognized as important to insects and are used for constructing taxonomies. However, multiple parameters affect the expression of CHCs besides a genetic component. We propose that selection may act differently on the expression of CHCs, depending on the evolutionary context. To explore the influence of selection, the CHCs of two closely related ant species, Lasius niger and Lasius platythorax, were studied in a multidisciplinary approach. We characterized (i) CHCs and (ii) niches (through baiting, activity observations and foraging analysis). The species were distinct in both measures, although to a varying degree. Although they showed moderate niche partitioning along diet and environmental preferences, chemical differences were unexpectedly pronounced. This may be explained by divergent selection on mate recognition cues or by other influences on CHCs. Such striking chemical differences among closely related species may not be the rule and suggest that taxonomies based on CHCs should be interpreted cautiously; though, they remain useful tools for differentiating among cryptic species.     

Polygenic evolution drives species divergence and climate adaptation in corals

Closely related species often show substantial differences in ecological traits that allow them to occupy different environmental niches. For few of these systems is it clear what the genomic basis of adaptation is and whether a few loci of major effect or many genome‐wide differences drive species divergence. Four cryptic species of the tabletop coral Acropora hyacinthus are broadly sympatric in American Samoa; here we show that two common species have differences in key environmental traits such as microhabitat distributions and thermal stress tolerance. We compared gene expression patterns and genetic polymorphism between these two species using RNA‐Seq. The vast majority of polymorphisms are shared between species, but the two species show widespread differences in allele frequencies and gene expression, and tend to host different symbiont types. We find that changes in gene expression are related to changes in the frequencies of many gene regulatory variants, but that many of these differences are consistent with the action of genetic drift. However, we observe greater genetic divergence between species in amino acid replacement polymorphisms compared to synonymous variants. These findings suggest that polygenic evolution plays a major role in driving species differences in ecology and resilience to climate change.     

Contrasting microsatellite variation between subalpine and western larch, two closely related species with different distribution patterns

Subalpine larch (Larix lyallii Parl.) and western larch (Larix occidentalis Nutt.) represent two closely related species with contrasting abundance and distribution patterns in Western North America. Genetic diversity at seven informative microsatellite loci was determined for 19 populations of subalpine larch and nine populations of western larch. Contrasting genetic diversity and patterns of population differentiation were observed between the two species. The overall within-population genetic diversity parameters were lower in subalpine larch (A = 3.2; A(P) = 3.6; H(E) = 0.418) than in western larch (A(P) = 5.51; H(E) = 0.580), a pattern that is likely related to historical or demographic factors. No evidence of interspecific hybridization was observed. Significantly more population differentiation (theta = 0.15; R(ST) = 0.07), consistent with more restricted gene flow, was observed for subalpine larch as compared to western larch (theta = 0.05; R(ST) = 0.04). Under the assumption of an infinite allele mutation model, 12 of the 19 subalpine larch populations showed signs of deviation from the mutation-drift equilibrium, which suggests Holocene population bottlenecks and fluctuations in effective population size for this species. None of the western larch populations deviated significantly from the mutation-drift equilibrium. For both species, Mantel's test revealed a significant positive relationship between geographical and genetic distances indicative of isolation by distance. A similar geographical structure was detected in both species, suggesting at least two genetically distinct glacial populations in each species. The various implications for gene conservation are discussed.     

Interspecific interactions influence contrasting spatial genetic structures in two closely related damselfly species

Spatial genetic structure (SGS) is largely determined by colonization history, landscape and ecological characteristics of the species. Therefore, sympatric and ecologically similar species are expected to exhibit similar SGSs, potentially enabling prediction of the SGS of one species from that of another. On the other hand, due to interspecific interactions, ecologically similar species could have different SGSs. We explored the SGSs of the closely related Calopteryx splendens and Calopteryx virgo within Finland and related the genetic patterns to characteristics of the sampling localities. We observed different SGSs for the two species. Genetic differentiation even within short distances in C. splendens suggests genetic drift as an important driver. However, we also observed indication of previous gene flow (revealed by a negative relationship between genetic differentiation and increasing potential connectivity of the landscape). Interestingly, genetic diversity of C. splendens was negatively related to density of C. virgo, suggesting that interspecific interactions influence the SGS of C. splendens. In contrast, genetic differentiation between C. virgo subpopulations was low and only exhibited relationships with latitude, pointing to high gene flow, colonization history and range margin effects as the drivers of SGS. The different SGSs of the two ecologically similar species caution indirect inferences of SGS based on ecologically similar surrogate species.     

Intrinsic incompatibilities evolving as a by‐product of divergent ecological selection: Considering them in empirical studies on divergence with gene flow

The possibility of intrinsic barriers to gene flow is often neglected in empirical research on local adaptation and speciation with gene flow, for example when interpreting patterns observed in genome scans. However, we draw attention to the fact that, even with gene flow, divergent ecological selection may generate intrinsic barriers involving both ecologically selected and other interacting loci. Mechanistically, the link between the two types of barriers may be generated by genes that have multiple functions (i.e., pleiotropy), and/or by gene interaction networks. Because most genes function in complex networks, and their evolution is not independent of other genes, changes evolving in response to ecological selection can generate intrinsic barriers as a by‐product. A crucial question is to what extent such by‐product barriers contribute to divergence and speciation—that is whether they stably reduce gene flow. We discuss under which conditions by‐product barriers may increase isolation. However, we also highlight that, depending on the conditions (e.g., the amount of gene flow and the strength of selection acting on the intrinsic vs. the ecological barrier component), the intrinsic incompatibility may actually destabilize barriers to gene flow. In practice, intrinsic barriers generated as a by‐product of divergent ecological selection may generate peaks in genome scans that cannot easily be interpreted. We argue that empirical studies on divergence with gene flow should consider the possibility of both ecological and intrinsic barriers. Future progress will likely come from work combining population genomic studies, experiments quantifying fitness and molecular studies on protein function and interactions.     

A parapatric propensity for breeding precludes the completion of speciation in common teal (Anas crecca, sensu lato)

Speciation is a process in which genetic drift and selection cause divergence over time. However, there is no rule dictating the time required for speciation, and even low levels of gene flow hinder divergence, so that taxa may be poised at the threshold of speciation for long periods of evolutionary time. We sequenced mitochondrial DNA (mtDNA) and eight nuclear introns (nuDNA) to estimate genomic levels of differentiation and gene flow between the Eurasian common teal (Anas crecca crecca) and the North American green‐winged teal (Anas crecca carolinensis). These ducks come into contact in Beringia (north‐eastern Asia and north‐western North America) and have probably done so, perhaps cyclically, since the Pliocene–Pleistocene transition, ～2.6 Ma, when they apparently began diverging. They have diagnosable differences in male plumage and are 6.9% divergent in the mtDNA control region, with only 1 of 58 crecca and 2 of 86 carolinensis having haplotypes grouping with the other. Two nuclear loci were likewise strongly structured between these teal (Φ_st ≥ 0.35), but six loci were undifferentiated or only weakly structured (Φ_st = 0.0–0.06). Gene flow between crecca and carolinensis was ～1 individual per generation in both directions in mtDNA, but was asymmetrical in nuDNA, with ～1 and ～20 individuals per generation immigrating into crecca and carolinensis, respectively. This study illustrates that species delimitation using a single marker oversimplifies the complexity of the speciation process, and it suggests that even with divergent selection, moderate levels of gene flow may stall the speciation process short of completion.     

Increased genetic divergence between two closely related fir species in areas of range overlap

Because of introgressive hybridization, closely related species can be more similar to each other in areas of range overlap (parapatry or sympatry) than in areas where they are geographically isolated from each other (allopatry). Here, we report the reverse situation based on nuclear genetic divergence between two fir species, Abies chensiensis and Abies fargesii, in China, at sites where they are parapatric relative to where they are allopatric. We examined genetic divergence across 126 amplified fragment length polymorphism (AFLP) markers in a set of 172 individuals sampled from both allopatric and parapatric populations of the two species. Our analyses demonstrated that AFLP divergence was much greater between the species when comparisons were made between parapatric populations than between allopatric populations. We suggest that selection in parapatry may have largely contributed to this increased divergence.     

Genetic divergence is decoupled from ecological diversification in the hawaiian nesosydne planthoppers

Adaptive radiation involves ecological shifts coupled with isolation of gene pools. However, we know little about what drives the initial stages of divergence. We study a system in which ecological diversification is found within a chronologically well-defined geological matrix to provide insight into this enigmatic phase of radiation. We tested the hypothesis that a period of geographic isolation precedes ecological specialization in an adaptive radiation of host-specialized Hawaiian planthoppers. We examined population structure and history using mitochondrial and multiple independent microsatellite loci in a species whose geographic distribution on the island of Hawaii enabled us to observe the chronology of divergence in its very earliest stages. We found that genetic divergence is associated with geographic features but not different plant hosts and that divergence times are very recent and on the same timescales as the dynamic geology of the island. Our results suggest an important role for geography in the dynamics of the early stages of divergence.