Host‐associated genomic differentiation in congeneric butterflies: now you see it,now you do not

Ecotypic variation among populations may become associated with widespread genomic differentiation, but theory predicts that this should happen only under particular conditions of gene flow, selection and population size. In closely related species, we might expect the strength of host‐associated genomic differentiation (HAD) to be correlated with the degree of phenotypic differentiation in host‐adaptive traits. Using microsatellite and Amplified Fragment Length Polymorphism (AFLP) markers, and controlling for isolation by distance between populations, we sought HAD in two congeneric species of butterflies with different degrees of host plant specialization. Prior work on Euphydryas editha had shown strong interpopulation differentiation in host‐adapted traits, resulting in incipient reproductive isolation among host‐associated ecotypes. We show here that Euphydryas aurinia had much weaker host‐associated phenotypic differentiation. Contrary to our expectations, we detected HAD in Euphydryas aurinia, but not in E. editha. Even within an E. aurinia population that fed on both hosts, we found weak but significant sympatric HAD that persisted in samples taken 9 years apart. The finding of significantly stronger HAD in the system with less phenotypic differentiation may seem paradoxical. Our findings can be explained by multiple factors, ranging from differences in dispersal or effective population size, to spatial variation in genomic or phenotypic traits and to structure induced by past histories of host‐adapted populations. Other infrequently measured factors, such as differences in recombination rates, may also play a role. Our result adds to recent work as a further caution against assumptions of simple relationships between genomic and adaptive phenotypic differentiation.     

Parallel evolution of site‐specific changes in divergent caribou lineages

The parallel evolution of phenotypes or traits within or between species provides important insight into the basic mechanisms of evolution. Genetic and genomic advances have allowed investigations into the genetic underpinnings of parallel evolution and the independent evolution of similar traits in sympatric species. Parallel evolution may best be exemplified among species where multiple genetic lineages, descended from a common ancestor, colonized analogous environmental niches, and converged on a genotypic or phenotypic trait. Modern North American caribou (Rangifer tarandus) originated from three ancestral sources separated during the Last Glacial Maximum (LGM): the Beringian–Eurasian lineage (BEL), the North American lineage (NAL), and the High Arctic lineage (HAL). Historical introgression between the NAL and the BEL has been found throughout Ontario and eastern Manitoba. In this study, we first characterized the functional differentiation in the cytochrome‐b (cytB) gene by identifying nonsynonymous changes. Second, the caribou lineages were used as a direct means to assess site‐specific parallel changes among lineages. There was greater functional diversity within the NAL despite the BEL having greater neutral diversity. The patterns of amino acid substitutions occurring within different lineages supported the parallel evolution of cytB amino acid substitutions suggesting different selective pressures among lineages. This study highlights the independent evolution of identical amino acid substitutions within a wide‐ranging mammal species that have diversified from different ancestral haplogroups and where ecological niches can invoke parallel evolution.     

High levels of interspecific gene flow in an endemic cichlid fish adaptive radiation from an extreme lake environment

Studying recent adaptive radiations in isolated insular systems avoids complicating causal events and thus may offer clearer insight into mechanisms generating biological diversity. Here, we investigate evolutionary relationships and genomic differentiation within the recent radiation of Alcolapia cichlid fish that exhibit extensive phenotypic diversification, and which are confined to the extreme soda lakes Magadi and Natron in East Africa. We generated an extensive RAD data set of 96 individuals from multiple sampling sites and found evidence for genetic admixture between species within Lake Natron, with the highest levels of admixture between sympatric populations of the most recently diverged species. Despite considerable environmental separation, populations within Lake Natron do not exhibit isolation by distance, indicating panmixia within the lake, although individuals within lineages clustered by population in phylogenomic analysis. Our results indicate exceptionally low genetic differentiation across the radiation despite considerable phenotypic trophic variation, supporting previous findings from smaller data sets; however, with the increased power of densely sampled SNPs, we identify genomic peaks of differentiation (F_ST outliers) between Alcolapia species. While evidence of ongoing gene flow and interspecies hybridization in certain populations suggests that Alcolapia species are incompletely reproductively isolated, the identification of outlier SNPs under diversifying selection indicates the radiation is undergoing adaptive divergence.     

An ecological genetic delineation of local seed‐source provenance for ecological restoration

An increasingly important practical application of the analysis of spatial genetic structure within plant species is to help define the extent of local provenance seed collection zones that minimize negative impacts in ecological restoration programs. Here, we derive seed sourcing guidelines from a novel range‐wide assessment of spatial genetic structure of 24 populations of Banksia menziesii (Proteaceae), a widely distributed Western Australian tree of significance in local ecological restoration programs. An analysis of molecular variance (AMOVA) of 100 amplified fragment length polymorphism (AFLP) markers revealed significant genetic differentiation among populations (Φ_PT = 0.18). Pairwise population genetic dissimilarity was correlated with geographic distance, but not environmental distance derived from 15 climate variables, suggesting overall neutrality of these markers with regard to these climate variables. Nevertheless, Bayesian outlier analysis identified four markers potentially under selection, although these were not correlated with the climate variables. We calculated a global R‐statistic using analysis of similarities (ANOSIM) to test the statistical significance of population differentiation and to infer a threshold seed collection zone distance of ~60 km (all markers) and 100 km (outlier markers) when genetic distance was regressed against geographic distance. Population pairs separated by >60 km were, on average, twice as likely to be significantly genetically differentiated than population pairs separated by <60 km, suggesting that habitat‐matched sites within a 30‐km radius around a restoration site genetically defines a local provenance seed collection zone for B. menziesii. Our approach is a novel probability‐based practical solution for the delineation of a local seed collection zone to minimize negative genetic impacts in ecological restoration.     

Exaggerated male forelegs are not more differentiated than wing morphology in two widespread sister species of black scavenger flies

Sexual selection represents a potent force that can drive rapid population differentiation in traits related to reproductive success. Hence, sexual traits are expected to show greater population divergence than non-sexual traits. We test this prediction by exploring patterns of morphological differentiation of the exaggerated fore femur (a male-specific sexual trait) and the wing (a non-sexual trait) among allopatric and sympatric populations of the widespread sister dung fly species Sepsis neocynipsea and Sepsis cynipsea (Diptera: Sepsidae). While both species occur in Eurasia, S. neocynipsea also abounds in North America, albeit previous studies suggest strong differentiation in morphology, behavior, and mating systems. To evaluate the degree of differentiation expected under neutrality between S. cynipsea, European S. neocynipsea, and North American S. neocynipsea, we genotyped 30 populations at nine microsatellite markers, revealing almost equal differentiation between and minor differentiation among geographic populations within the three lineages. Landmark-based analysis of 18 populations reared at constant 18 and 24°C in a laboratory common garden revealed moderate temperature-dependent phenotypic plasticity and significant heritable differentiation in size and shape of male forelegs and wings among iso-female lines of the three lineages. Following the biological species concept, there was weaker differentiation between cross-continental populations of S. neocynipsea relative to S. cynipsea, and more fore femur differentiation between the two species in sympatry versus allopatry (presumably due to character displacement). Contrary to expectation, wing morphology showed as much shape differentiation between evolutionary independent lineages as fore femora, providing no evidence for faster diversification of traits primarily engaged in mating.     

Biogeography and host‐related factors trump parasite life history: limited congruence among the genetic structures of specific ectoparasitic lice and their rodent hosts

Parasites and hosts interact across both micro‐ and macroevolutionary scales where congruence among their phylogeographic and phylogenetic structures may be observed. Within southern Africa, the four‐striped mouse genus, Rhabdomys, is parasitized by the ectoparasitic sucking louse, Polyplax arvicanthis. Molecular data recently suggested the presence of two cryptic species within P. arvicanthis that are sympatrically distributed across the distributions of four putative Rhabdomys species. We tested the hypotheses of phylogeographic congruence and cophylogeny among the two parasite lineages and the four host taxa, utilizing mitochondrial and nuclear sequence data. Despite the documented host‐specificity of P. arvicanthis, limited phylogeographic correspondence and nonsignificant cophylogeny was observed. Instead, the parasite–host evolutionary history is characterized by limited codivergence and several duplication, sorting and host‐switching events. Despite the elevated mutational rates found for P. arvicanthis, the spatial genetic structure was not more pronounced in the parasite lineages compared with the hosts. These findings may be partly attributed to larger effective population sizes of the parasite lineages, the vagility and social behaviour of Rhabdomys, and the lack of host‐specificity observed in areas of host sympatry. Further, the patterns of genetic divergence within parasite and host lineages may also be largely attributed to historical biogeographic changes (expansion‐contraction cycles). It is thus evident that the association between P. arvicanthis and Rhabdomys has been shaped by the synergistic effects of parasite traits, host‐related factors and biogeography over evolutionary time.     

Are sympatrically speciating Midas cichlid fish special? Patterns of morphological and genetic variation in the closely related species Archocentrus centrarchus

Established empirical cases of sympatric speciation are scarce, although there is an increasing consensus that sympatric speciation might be more common than previously thought. Midas cichlid fish are one of the few substantiated cases of sympatric speciation, and they formed repeated radiations in crater lakes. In contrast, in the same environment, such radiation patterns have not been observed in other species of cichlids and other families of fish. We analyze morphological and genetic variation in a cichlid species (Archocentrus centrarchus) that co‐inhabits several crater lakes with the Midas species complex. In particular, we analyze variation in body and pharyngeal jaw shape (two ecologically important traits in sympatrically divergent Midas cichlids) and relate that to genetic variation in mitochondrial control region and microsatellites. Using these four datasets, we analyze variation between and within two Nicaraguan lakes: a crater lake where multiple Midas cichlids have been described and a lake where the source population lives. We do not observe any within‐lake clustering consistent across morphological traits and genetic markers, suggesting the absence of sympatric divergence in A. centrarchus. Genetic differentiation between lakes was low and morphological divergence absent. Such morphological similarity between lakes is found not only in average morphology, but also when analyzing covariation between traits and degree of morphospace occupation. A combined analysis of the mitochondrial control region in A. centrarchus and Midas cichlids suggests that a difference between lineages in the timing of crater lake colonization cannot be invoked as an explanation for the difference in their levels of diversification. In light of our results, A. centrarchus represents the ideal candidate to study the genomic differences between these two lineages that might explain why some lineages are more likely to speciate and diverge in sympatry than others.     

The interplay of colour and bioacoustic traits in the differentiation of a Southeast Asian songbird complex

Morphological traits have served generations of biologists as a taxonomic indicator, and have been the main basis for defining and classifying species diversity for centuries. A quantitative integration of behavioural characters, such as vocalizations, in studies on biotic differentiation has arisen more recently, and the relative importance of these different traits in the diversification process remains poorly understood. To provide a framework within which to interpret the evolutionary interplay between morphological and behavioural traits, we generated a draft genome of a cryptic Southeast Asian songbird, the limestone wren‐babbler Napothera crispifrons. We resequenced whole genomes of multiple individuals of all three traditional subspecies and of a distinct leucistic population. We demonstrate strong genomic and mitochondrial divergence among all three taxa, pointing to the existence of three species‐level lineages. Despite its great phenotypic distinctness, the leucistic population was characterized by shallow genomic differentiation from its neighbour, with only a few localized regions emerging as highly diverged. Quantitative bioacoustic analysis across multiple traits revealed deep differences especially between the two taxa characterized by limited plumage differentiation. Our study demonstrates that differentiation in these furtive songbirds has resulted in a complex mosaic of colour‐based and bioacoustic differences among populations. Extreme colour differences can be anchored in few genomic loci and may therefore arise and subside rapidly.     

Genetic changes in a novel breeding population of Brassica napus synthesized from hundreds of crosses between B. rapa and B. carinata

Introgression of genomic variation between and within related crop species is a significant evolutionary approach for population differentiation, genome reorganization and trait improvement. Using the Illumina Infinium Brassica 60K SNP array, we investigated genomic changes in a panel of advanced generation new‐type Brassica napus breeding lines developed from hundreds of interspecific crosses between 122 Brassica rapa and 74 Brassica carinata accessions, and compared them with representative accessions of their three parental species. The new‐type B. napus population presented rich genetic diversity and abundant novel genomic alterations, consisting of introgressions from B. rapa and B. carinata, novel allelic combinations, reconstructed linkage disequilibrium patterns and haplotype blocks, and frequent deletions and duplications (nonrandomly distributed), particularly in the C subgenome. After a much shorter, but very intensive, selection history compared to traditional B. napus, a total of 15 genomic regions with strong selective sweeps and 112 genomic regions with putative signals of selective sweeps were identified. Some of these regions were associated with important agronomic traits that were selected for during the breeding process, while others were potentially associated with restoration of genome stability and fertility after interspecific hybridization. Our results demonstrate how a novel method for population‐based crop genetic improvement can lead to rapid adaptation, restoration of genome stability and positive responses to artificial selection.     

Contrasting genomic and phenotypic outcomes of hybridization between pairs of mimetic butterfly taxa across a suture zone

Hybrid zones, whereby divergent lineages come into contact and eventually hybridize, can provide insights on the mechanisms involved in population differentiation and reproductive isolation, and ultimately speciation. Suture zones offer the opportunity to compare these processes across multiple species. In this paper we use reduced‐complexity genomic data to compare the genetic and phenotypic structure and hybridization patterns of two mimetic butterfly species, Ithomia salapia and Oleria onega (Nymphalidae: Ithomiini), each consisting of a pair of lineages differentiated for their wing colour pattern and that come into contact in the Andean foothills of Peru. Despite similarities in their life history, we highlight major differences, both at the genomic and phenotypic level, between the two species. These differences include the presence of hybrids, variations in wing phenotype, and genomic patterns of introgression and differentiation. In I. salapia, the two lineages appear to hybridize only rarely, whereas in O. onega the hybrids are not only more common, but also genetically and phenotypically more variable. We also detected loci statistically associated with wing colour pattern variation, but in both species these loci were not over‐represented among the candidate barrier loci, suggesting that traits other than wing colour pattern may be important for reproductive isolation. Our results contrast with the genomic patterns observed between hybridizing lineages in the mimetic Heliconius butterflies, and call for a broader investigation into the genomics of speciation in Ithomiini ‐ the largest radiation of mimetic butterflies.     

Genetic structure and host–parasite co‐divergence: evidence for trait‐specific local adaptation

Host–parasite co‐evolution can lead to genetic differentiation among isolated host–parasite populations and local adaptation between parasites and their hosts. However, tests of local adaptation rarely consider multiple fitness‐related traits although focus on a single component of fitness can be misleading. Here, we concomitantly examined genetic structure and co‐divergence patterns of the trematode Coitocaecum parvum and its crustacean host Paracalliope fluviatilis among isolated populations using the mitochondrial cytochrome oxidase I gene (COI). We then performed experimental cross‐infections between two genetically divergent host–parasite populations. Both hosts and parasites displayed genetic differentiation among populations, although genetic structure was less pronounced in the parasite. Data also supported a co‐divergence scenario between C. parvum and P. fluviatilis potentially related to local co‐adaptation. Results from cross‐infections indicated that some parasite lineages seemed to be locally adapted to their sympatric (home) hosts in which they achieved higher infection and survival rates than in allopatric (away) amphipods. However, local, intrinsic host and parasite characteristics (host behavioural or immunological resistance to infections, parasite infectivity or growth rate) also influenced patterns of host–parasite interactions. For example, overall host vulnerability to C. parvum varied between populations, regardless of parasite origin (local vs. foreign), potentially swamping apparent local co‐adaptation effects. Furthermore, local adaptation effects seemed trait specific; different components of parasite fitness (infection and survival rates, growth) responded differently to cross‐infections. Overall, data show that genetic differentiation is not inevitably coupled with local adaptation, and that the latter must be interpreted with caution in a multi‐trait context.     

Short‐range phenotypic divergence among genetically distinct parapatric populations of an Australian funnel‐web spider

Speciation involves divergence at genetic and phenotypic levels. Where substantial genetic differentiation exists among populations, examining variation in multiple phenotypic characters may elucidate the mechanisms by which divergence and speciation unfold. Previous work on the Australian funnel‐web spider Atrax sutherlandi Gray (2010; Records of the Australian Museum 62 , 285–392; Mygalomorphae: Hexathelidae: Atracinae) has revealed a marked genetic structure along a 110‐kilometer transect, with six genetically distinct, parapatric populations attributable to past glacial cycles. In the present study, we explore variation in three classes of phenotypic characters (metabolic rate, water loss, and morphological traits) within the context of this phylogeographic structuring. Variation in metabolic and water loss rates shows no detectable association with genetic structure; the little variation observed in these rates may be due to the spiders’ behavioral adaptations (i.e., burrowing), which buffer the effects of climatic gradients across the landscape. However, of 17 morphological traits measured, 10 show significant variation among genetic populations, in a disjunct manner that is clearly not latitudinal. Moreover, patterns of variation observed for morphological traits serving different organismic functions (e.g., prey capture, burrowing, and locomotion) are dissimilar. In contrast, a previous study of an ecologically similar sympatric spider with little genetic structure indicated a strong latitudinal response in 10 traits over the same range. The congruence of morphological variation with deep phylogeographic structure in Tallaganda's A. sutherlandi populations, as well as the inconsistent patterns of variation across separate functional traits, suggest that the spiders are likely in early stages of speciation, with parapatric populations independently responding to local selective forces.     

Habitat fragmentation influences genetic diversity and differentiation: Fine‐scale population structure of Cercis canadensis (eastern redbud)

Forest fragmentation may negatively affect plants through reduced genetic diversity and increased population structure due to habitat isolation, decreased population size, and disturbance of pollen‐seed dispersal mechanisms. However, in the case of tree species, effective pollen‐seed dispersal, mating system, and ecological dynamics may help the species overcome the negative effect of forest fragmentation. A fine‐scale population genetics study can shed light on the postfragmentation genetic diversity and structure of a species. Here, we present the genetic diversity and population structure of Cercis canadensis L. (eastern redbud) wild populations on a fine scale within fragmented areas centered around the borders of Georgia–Tennessee, USA. We hypothesized high genetic diversity among the collections of C. canadensis distributed across smaller geographical ranges. Fifteen microsatellite loci were used to genotype 172 individuals from 18 unmanaged and naturally occurring collection sites. Our results indicated presence of population structure, overall high genetic diversity (H_E = 0.63, H_O = 0.34), and moderate genetic differentiation (F_ST = 0.14) among the collection sites. Two major genetic clusters within the smaller geographical distribution were revealed by STRUCTURE. Our data suggest that native C. canadensis populations in the fragmented area around the Georgia–Tennessee border were able to maintain high levels of genetic diversity, despite the presence of considerable spatial genetic structure. As habitat isolation may negatively affect gene flow of outcrossing species across time, consequences of habitat fragmentation should be regularly monitored for this and other forest species. This study also has important implications for habitat management efforts and future breeding programs.     

A trait‐based approach to predict population genetic structure in bees

Understanding population genetic structure is key to developing predictions about species susceptibility to environmental change, such as habitat fragmentation and climate change. It has been theorized that life‐history traits may constrain some species in their dispersal and lead to greater signatures of population genetic structure. In this study, we use a quantitative comparative approach to assess if patterns of population genetic structure in bees are driven by three key species‐level life‐history traits: body size, sociality, and diet breadth. Specifically, we reviewed the current literature on bee population genetic structure, as measured by the differentiation indices Nei's G_ST, Hedrick's G′_ST, and Jost's D. We then used phylogenetic generalised linear models to estimate the correlation between the evolution of these traits and patterns of genetic differentiation. Our analyses revealed a negative and significant effect of body size on genetic structure, regardless of differentiation index utilized. For Hedrick's G′_ST and Jost's D, we also found a significant impact of sociality, where social species exhibited lower levels of differentiation than solitary species. We did not find an effect of diet specialization on population genetic structure. Overall, our results suggest that physical dispersal or other functions related to body size are among the most critical for mediating population structure for bees. We further highlight the importance of standardizing population genetic measures to more easily compare studies and to identify the most susceptible species to landscape and climatic changes.     

Unexpected cryptic species diversity in the widespread coral Seriatopora hystrix masks spatial‐genetic patterns of connectivity

Mounting evidence of cryptic species in a wide range of taxa highlights the need for careful analyses of population genetic data sets to unravel within‐species diversity from potential interspecies relationships. Here, we use microsatellite loci and hierarchical clustering analysis to investigate cryptic diversity in sympatric and allopatric (separated by 450 km) populations of the widespread coral Seriatopora hystrix on the Great Barrier Reef. Structure analyses delimited unique genetic clusters that were confirmed by phylogenetic and extensive population‐level analyses. Each of four sympatric yet distinct genetic clusters detected within S. hystrix demonstrated greater genetic cohesion across regional scales than between genetic clusters within regions (<10 km). Moreover, the magnitude of genetic differentiation between different clusters (>0.620 G”_ST) was similar to the difference between S. hystrix clusters and the congener S. caliendrum (mean G”_ST 0.720). Multiple lines of evidence, including differences in habitat specificity, mitochondrial identity, Symbiodinium associations and morphology, corroborate the nuclear genetic evidence that these distinct clusters constitute different species. Hierarchical clustering analysis combined with more traditional population genetic methods provides a powerful approach for delimiting species and should be regularly applied to ensure that ecological and evolutionary patterns interpreted for single species are not confounded by the presence of cryptic species.     

Population genomics reveals multiple drivers of population differentiation in a sex‐role‐reversed pipefish

A major goal of molecular ecology is to identify the causes of genetic and phenotypic differentiation among populations. Population genomics is suitably poised to tackle these key questions by diagnosing the evolutionary mechanisms driving divergence in nature. Here, we set out to investigate the evolutionary processes underlying population differentiation in the Gulf pipefish, Syngnathus scovelli. We sampled approximately 50 fish from each of 12 populations distributed from the Gulf coast of Texas to the Atlantic coast of Florida and performed restriction‐site‐associated DNA sequencing to identify SNPs throughout the genome. After imposing quality and stringency filters, we selected a panel of 6348 SNPs present in all 12 populations, 1753 of which were not physically linked. We identified a genome‐wide pattern of isolation by distance, in addition to a more substantial genetic break separating populations in the Gulf of Mexico from those in the Atlantic. We also used several divergence outlier approaches and tests for genotype–environment correlations to identify 400 SNPs putatively involved in local adaptation. Patterns of phenotypic differentiation and variation diverged from the overall genomic pattern, suggesting that selection, phenotypic plasticity or demographic factors may be shaping phenotypes in distinct populations. Overall, our results suggest that population divergence is driven by a variety of factors in S. scovelli, including neutral processes and selection on multiple traits.     

Large‐scale adaptive divergence in Boechera fecunda,an endangered wild relative of Arabidopsis

Many biological species are threatened with extinction because of a number of factors such as climate change and habitat loss, and their preservation depends on an accurate understanding of the extent of their genetic variability within and among populations. In this study, we assessed the genetic divergence of five quantitative traits in 10 populations of an endangered cruciferous species, Boechera fecunda, found in only several populations in each of two geographic regions (WEST and EAST) in southwestern Montana. We analyzed variation in quantitative traits, neutral molecular markers, and environmental factors and provided evidence that despite the restricted geographical distribution of this species, it exhibits a high level of genetic variation and regional adaptation. Conservation efforts therefore should be directed to the preservation of populations in each of these two regions without attempting transplantation between regions. Heritabilities and genetic coefficients of variation estimated from nested ANOVAs were generally high for leaf and rosette traits, although lower (and not significantly different from 0) for water‐use efficiency. Measures of quantitative genetic differentiation, Q_ST, were calculated for each trait from each pair of populations. For three of the five traits, these values were significantly higher between regions compared with those within regions (after adjustment for neutral genetic variation, F_ST). This suggested that natural selection has played an important role in producing regional divergence in this species. Our analysis also revealed that the B. fecunda populations appear to be locally adapted due, at least in part, to differences in environmental conditions in the EAST and WEST regions.     

Prey‐associated head‐size variation in an invasive lizard in the Hawaiian Islands

Biological invasions are recognized as a primary driver of large‐scale changes in global ecosystems. This study addresses ecomorphological variation in head size within and among populations of an ecologically destructive invasive predator, and evaluates the potential roles of environmental components in phenotypic differentiation. We used four size‐corrected measurements of head morphology in Jackson's chameleons, Trioceros jacksonii xantholophus (N = 319), collected from multiple Hawaiian Islands to assess phenotypic variation among and within islands. Results of analysis of variance (ANOVA) comparing chameleon head size (PC1) among islands revealed significant differences (mean difference > 5%) associated with variation in both rainfall and diet composition using Mann–Whitney U‐tests and chi‐squared analyses. These results suggest that morphological differentiation among populations from different islands has occurred over a relatively short ecological timescale, and is likely the result of ecomorphological adaptation to differences in exploited prey hardness. Intra‐island allopatric population variation, however, was also detected in this study. Although we might expect that genetic change is the more likely explanation for differences between islands than within, and that plasticity may be more likely an explanation for the within‐ than the between‐island differences, it is also possible that both within‐ and between‐island patterns are the results of genetic change, or of plasticity.     

Long‐term stasis and short‐term divergence in the phenotypes of microsnails on oceanic islands

Phenotypic divergence is often unrelated to genotypic divergence. An extreme example is rapid phenotypic differentiation despite genetic similarity. Another extreme is morphological stasis despite substantial genetic divergence. These opposite patterns have been viewed as reflecting opposite properties of the lineages. In this study, phenotypic radiation accompanied by both rapid divergence and long‐term conservatism is documented in the inferred molecular phylogeny of the micro land snails Cavernacmella (Assimineidae) on the Ogasawara Islands. The populations of Cavernacmella on the Sekimon limestone outcrop of Hahajima Island showed marked divergence in shell morphology. Within this area, one lineage diversified into types with elongated turret shells, conical shells and flat disc‐like shells without substantial genetic differentiation. Additionally, a co‐occurring species with these types developed a much larger shell size. Moreover, a lineage adapted to live inside caves in this area. In contrast, populations in the other areas exhibited no morphological differences despite high genetic divergence among populations. Accordingly, the phenotypic evolution of Cavernacmella in Ogasawara is characterized by a pattern of long‐term stasis and periodic bursts of change. This pattern suggests that even lineages with phenotypic conservatism could shift to an alternative state allowing rapid phenotypic divergence.     

Chromosomal inversions and ecotypic differentiation in Anopheles gambiae: the perspective from whole‐genome sequencing

The molecular mechanisms and genetic architecture that facilitate adaptive radiation of lineages remain elusive. Polymorphic chromosomal inversions, due to their recombination‐reducing effect, are proposed instruments of ecotypic differentiation. Here, we study an ecologically diversifying lineage of Anopheles gambiae, known as the Bamako chromosomal form based on its unique complement of three chromosomal inversions, to explore the impact of these inversions on ecotypic differentiation. We used pooled and individual genome sequencing of Bamako, typical (non‐Bamako) An. gambiae and the sister species Anopheles coluzzii to investigate evolutionary relationships and genomewide patterns of nucleotide diversity and differentiation among lineages. Despite extensive shared polymorphism and limited differentiation from the other taxa, Bamako clusters apart from the other taxa, and forms a maximally supported clade in neighbour‐joining trees based on whole‐genome data (including inversions) or solely on collinear regions. Nevertheless, F_ST outlier analysis reveals that the majority of differentiated regions between Bamako and typical An. gambiae are located inside chromosomal inversions, consistent with their role in the ecological isolation of Bamako. Exceptionally differentiated genomic regions were enriched for genes implicated in nervous system development and signalling. Candidate genes associated with a selective sweep unique to Bamako contain substitutions not observed in sympatric samples of the other taxa, and several insecticide resistance gene alleles shared between Bamako and other taxa segregate at sharply different frequencies in these samples. Bamako represents a useful window into the initial stages of ecological and genomic differentiation from sympatric populations in this important group of malaria vectors.