Phylogeography and natural selection in the Tenerife gecko Tarentola delalandii: testing historical and adaptive hypotheses

Molecular analyses of the Kalotermes dispar-complex (Blattodea: Kalotermitidae) from the Canary Islands reveal cryptic intraspecific divergence and a connection to a lone Nearctic congener

The Canary Islands is a Macaronesian volcanic archipelago with a depauperate community of three species of Kalotermitidae, including Kalotermes dispar. A total of 54 Kalotermes colonies were collected from Gran Canaria, Tenerife, La Gomera, La Palma, and El Hierro islands. Soldiers and imagos were morphologically examined and sequenced for four mitochondrial markers. Although morphological differences could not be detected, phylogenetic analysis of both cox1/tRNA/cox2 and rrnL markers revealed two distinct clades of K. dispar, suggesting cryptic diversity. The diversification within the Canary Kalotermes lineage most likely occurred around 7.5 Mya, while the divergence within the two clades was reconstructed at about 3.6 Mya and 1.9 Mya. Kalotermes approximatus from the southeastern Nearctic constitutes a sister to the Canary Kalotermes, while the Palearctic K. flavicollis, K. italicus, and K. phoenicae form a separate clade. It is hypothesized that a faunal exchange of Kalotermes from the Nearctic to the Canary Islands occurred via transoceanic rafting during the mid-Miocene.

Habitat-associated Genomic Variation in a Wall Lizard from an Oceanic Island

The lizard Teira dugesii exhibits morphological divergence between beach and inland habitats in the face of gene flow, within the volcanic island of Madeira, Portugal. Here, we analyzed genomic data obtained by genotyping-by-sequencing, which provided 16,378 single nucleotide polymorphisms (SNPs) from 94 individuals sampled from 15 sites across Madeira. Ancient within-island divergence in allopatry appears to have mediated divergence in similar species within other Atlantic islands, but this hypothesis was not supported for T. dugesii. Across all samples, a total of 168 SNPs were classified as statistical outliers using pcadapt and OutFLANK. Redundancy analysis (RDA) revealed that 17 of these outliers were associated with beach/inland habitats. The SNPs were located within 16 sequence tags and 15 of these were homologous with sequences in a 31 Mb region on chromosome 3 of a reference wall lizard genome (the remaining tag could not be associated with any chromosome). We further investigated outliers through contingency analyses of allele frequencies at each of four pairs of adjacent beach-inland sites. The majority of the outliers detected by the RDA were confirmed at two pairs of these matched sites. These analyses also suggested some parallel divergence at different localities. Six other outliers were associated with site elevation, four of which were located on chromosome 5 of the reference genome. Our study lends support to a previous hypothesis that divergent selection between gray shingle beaches and inland regions overcomes gene flow and leads to the observed morphological divergence between populations in these adjacent habitats.

Genetic Diversities and Historical Dynamics of Native Ethiopian Horse Populations (Equus caballus) Inferred from Mitochondrial DNA Polymorphisms

Matrilineal genetic diversity and relationship were investigated among eight morphologically identified native Ethiopian horse populations using polymorphisms in 46 mtDNA D-loop sequences (454 base pairs). The horse populations identified were Abyssinian, Bale, Borana, Horro, Kafa, Kundido feral horses, Ogaden and Selale. Mitochondrial DNA D-loop sequences were characterized by 15 variable sites that defined five different haplotypes. All genetic diversity estimates, including Reynolds' linearized genetic distance, genetic differentiation (F_ST) and nucleotide sequence divergence (D_A), revealed a low genetic differentiation in native Ethiopian horse populations. However, Kundido feral and Borana domestic horses were slightly diverged from the rest of the Ethiopian horse populations. We also tried to shed some light on the matrilineal genetic root of native Ethiopian horses from a network constructed by combining newly generated haplotypes and reference haplotypes deposited in the GenBank for Eurasian type Turkish Anatolian horses that were used as a genetic conduit between Eurasian and African horse populations. Ninety-two haplotypes were generated from the combined Ethio-Eurasian mtDNA D-loop sequences. A network reconstructed from the combined haplotypes using Median-Joining algorithm showed that haplotypes generated from native Ethiopian horses formed separate clusters. The present result encourages further investigation of the genetic origin of native African horses by retrieving additional mtDNA sequences deposited in the GenBank for African and Eurasian type horses.

Ecological, behavioral, and phylogenetic influences on the evolution of dorsal color pattern in geckos

The dorsal surfaces of many taxonomic groups often feature repetitive pattern elements consisting of stripes, spots, or bands. Here, we investigate how distinct categories of camouflage pattern work by relating them to ecological and behavioral traits in 439 species of gecko. We use phylogenetic comparative methods to test outstanding hypotheses based on camouflage theory and research in other taxa. We found that bands are associated with nocturnal activity, suggesting bands provide effective camouflage for motionless geckos resting in refugia during the day. A predicted association between stripes and diurnal activity was not supported, suggesting that stripes do not work via dazzle camouflage mechanisms in geckos. This, along with a lack of support for our prediction that plain patterning should be associated with open habitats, suggests that similar camouflage patterns do not work in consistent ways across taxa. We also found that plain and striped lineages frequently switched between using open or closed habitats, whereas spotted lineages rarely transitioned. This suggests that pattern categories differ in how specialized or generalized their camouflage is. This result has ramifications for theory on how camouflage compromises to background heterogeneity and how camouflage pattern might influence evolutionary trajectories.

Along for the ride or missing it altogether: exploring the host specificity and diversity of haemogregarines in the Canary Islands

Background

Host-parasite relationships are expected to be strongly shaped by host specificity, a crucial factor in parasite adaptability and diversification. Because whole host communities have to be considered to assess host specificity, oceanic islands are ideal study systems given their simplified biotic assemblages. Previous studies on insular parasites suggest host range broadening during colonization. Here, we investigate the association between one parasite group (haemogregarines) and multiple sympatric hosts (of three lizard genera: Gallotia, Chalcides and Tarentola) in the Canary Islands. Given haemogregarine characteristics and insular conditions, we hypothesized low host specificity and/or occurrence of host-switching events.

Methods

A total of 825 samples were collected from the three host taxa inhabiting the seven main islands of the Canarian Archipelago, including locations where the different lizards occurred in sympatry. Blood slides were screened to assess prevalence and parasitaemia, while parasite genetic diversity and phylogenetic relationships were inferred from 18S rRNA gene sequences.

Results

Infection levels and diversity of haplotypes varied geographically and across host groups. Infections were found in all species of Gallotia across the seven islands, in Tarentola from Tenerife, La Gomera and La Palma, and in Chalcides from Tenerife, La Gomera and El Hierro. Gallotia lizards presented the highest parasite prevalence, parasitaemia and diversity (seven haplotypes), while the other two host groups (Chalcides and Tarentola) harbored one haplotype each, with low prevalence and parasitaemia levels, and very restricted geographical ranges. Host-sharing of the same haemogregarine haplotype was only detected twice, but these rare instances likely represent occasional cross-infections.

Conclusions

Our results suggest that: (i) Canarian haemogregarine haplotypes are highly host-specific, which might have restricted parasite host expansion; (ii) haemogregarines most probably reached the Canary Islands in three colonization events with each host genus; and (iii) the high number of parasite haplotypes infecting Gallotia hosts and their restricted geographical distribution suggest co-diversification. These findings contrast with our expectations derived from results on other insular parasites, highlighting how host specificity depends on parasite characteristics and evolutionary history.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-2760-5) contains supplementary material, which is available to authorized users.

Predictability in evolution: Adaptation of the Bonaire anole (Anolis bonairensis) to an extreme environment

The extent to which evolution is deterministic (predictable), or random, is a fundamental question in evolution. This case study attempts to determine the extent to which interspecific divergence can be predicted from intraspecific trends related species. The mountainous Lesser Antilles are occupied by one or two anole species with very substantial intraspecific differences in the quantitative traits between xeric and rainforest habitats. These ecologically determined differences tend to be in parallel in each island species. A related species (Anolis bonairensis) lives on the far more xeric island of Bonaire, and this study tests the extent to which its interspecific divergence in hue and pattern traits can be predicted from the parallel intraspecific variation exhibited in Lesser Antillean anoles. Regression against a multivariate climate variable suggests that the hue and pattern of the Bonaire anole are consistently predicted from the ecologically determined intraspecific variation of its Lesser Antillean relatives. However, this predictability may be less consistent with other character systems, for example, scalation.

Genetic structure of Micromeria (Lamiaceae) in Tenerife, the imprint of geological history and hybridization on within-island diversification

Geological history of oceanic islands can have a profound effect on the evolutionary history of insular flora, especially in complex islands such as Tenerife in the Canary Islands. Tenerife results from the secondary connection of three paleo‐islands by a central volcano, and other geological events that further shaped it. This geological history has been shown to influence the phylogenetic history of several taxa, including genus Micromeria (Lamiaceae). Screening 15 microsatellite markers in 289 individuals representing the eight species of Micromeria present in Tenerife, this study aims to assess the genetic diversity and structure of these species and its relation with the geological events on the island. In addition, we evaluate the extent of hybridization among species and discuss its influence on the speciation process. We found that the species restricted to the paleo‐islands present lower levels of genetic diversity but the highest levels of genetic differentiation suggesting that their ranges might have contracted over time. The two most widespread species in the island, M. hyssopifolia and M. varia, present the highest genetic diversity levels and a genetic structure that seems correlated with the geological composition of the island. Samples from M. hyssopifolia from the oldest paleo‐island, Adeje, appear as distinct while samples from M. varia segregate into two main clusters corresponding to the paleo‐islands of Anaga and Teno. Evidence of hybridization and intraspecific migration between species was found. We argue that species boundaries would be retained despite hybridization in response to the habitat's specific conditions causing postzygotic isolation and preserving morphological differentiation.

Predation drives interpopulation differences in parental care expression

Expressing parental care after oviposition or parturition is usually an obligate (evolved) trait within a species, despite evolutionary theory predicting that widespread species should vary in whether or not they express parental care according to local selection pressures. The lizard Eutropis longicaudata expresses maternal care only in a single population throughout its large geographical range, but why this pattern occurs is unknown. We used reciprocal translocation and predator exclusion experiments to test whether this intraspecific variation is a fixed trait within populations and whether predator abundance explains this perplexing pattern. Wild-caught female lizards that were reciprocally translocated consistently guarded or abandoned eggs in line with their population of origin. By contrast, most lizards raised in a common garden environment and subsequently released as adults adopted the maternal care strategy of the recipient population, even when the parents originated from a population lacking maternal care. Egg predation represents a significant fitness cost in the populations where females display egg-guarding behaviour, but guarding eggs outweighs this potential cost by increasing hatching success. These results imply that predators can be a driving force in the expression of parental care in instances where it is normally absent and that local selection pressure is sufficient to cause behavioural divergence in whether or not parental care is expressed.

Little and large: body size and genetic clines in a New Zealand gecko (Woodworthia maculata) along a coastal transect

Clinal variation can result from primary differentiation or secondary contact and determining which of these two processes is responsible for the existence of a cline is not a trivial problem. Samples from a coastal transect of New Zealand geckos (Woodworthia maculata) identified for the first time a body size cline 7–10 km wide. The larger geckos are almost twice the mass of the small adult geckos. Clines in allele and haplotype frequency were found at two of the four genetic loci examined. Estimated width of the morphological cline was concordant with neither the narrower mtDNA cline (3–7 m) nor the wider nuclear cline (RAG-2; 34–42 km), and cline centers were not coincident. Although the body size cline is narrow compared to the entire range of the species, it is 2–3 orders of magnitude greater than estimates of dispersal distance per generation for these geckos. No evidence of assortative mating, nor of hybrid disadvantage was identified, thus there is little evidence to infer that endogenous selection is maintaining a hybrid zone. We cannot distinguish secondary contact from primary origin of this body size cline but conclude that secondary contact is likely due to the occurrence of mtDNA haplotypes from three distinct clades within the coastal transect and the presence of two frequency clines within this region.

Factors determining the direction of ecological specialization in snail-feeding carabid beetles

A stout-slender dimorphism in body shape is observed among carabid beetles of the subtribe Carabina, which feed on land snails. We hypothesized that this dimorphism has resulted from divergent ecological specialization for feeding on different-sized land snails. Therefore, we examined whether the geographic variation in the body shape of Damaster blaptoides, a representative snail-feeding species in Japan, is correlated with the size of Euhadra, a genus of land snails frequently consumed by D. blaptoides. An analysis of beetle specimens from the whole distribution area of D. blaptoides determined that more slender beetle populations occurred in localities harboring larger snails, whereas more stout beetles inhabited localities harboring smaller snails. This pattern could be adaptive because slender beetles exhibit high feeding performance for large snails by inserting their heads into the shells, whereas stout beetles do so for small snails by crushing the shells. The D. blaptoides populations showed a clear genetic isolation-by-distance pattern, which could be effective in promoting such local adaptation. Thus, food resources as well as geographic isolation may have promoted adaptive divergence of external morphology in the snail-feeding carabid beetles.

The divergence of two independent lineages of an endemic Chinese gecko, Gekko swinhonis, launched by the Qinling orogenic belt

The genetic structure and demographic history of an endemic Chinese gecko, Gekko swinhonis, were investigated by analysing the mitochondrial cytochrome b gene and 10 microsatellite loci for samples collected from 27 localities. Mitochondrial DNA data provided a detailed distribution of two highly divergent evolutionary lineages, between which the average pairwise distance achieved was 0.14. The geographic division of the two lineages coincided with a plate boundary consisting of the Qinling and Taihang Mts, suggesting a historical vicariant pattern. The orogeny of the Qinling Mts, a dispersal and major climatic barrier of the region, may have launched the independent lineage divergence. Both lineages have experienced recent expansion, and the current sympatric localities comprised the region of contact between the lineages. Individual-based phylogenetic analyses of nucDNA and Bayesian-clustering approaches revealed a deep genetic structure analogous to mtDNA. Incongruence between nucDNA and mtDNA at the individual level at localities outside of the contact region can be explained by the different inheritance patterns and male-biased dispersal in this species. High genetic divergence, long-term isolation and ecological adaptation, as well as the morphological differences, suggest the presence of a cryptic species.

Nuclear and mitochondrial sequences confirm complex colonization patterns and clear species boundaries for flightless weevils in the Galápagos archipelago

Nuclear sequence data were collected from endemic Galápagos species and an introduced close relative, and contrasted with mitochondrial DNA sequences, continuing investigation into the colonization history and modes of diversification in the weevil genus Galapaganus. The current combined phylogeny together with previously published penalized likelihood age estimates builds a complex picture of the archipelago's colonization history. The present reconstruction relies on submerged platforms to explain the early divergence of the young southern Isabela endemics or the Española or San Cristobal populations. Diversity is later built through inter-island divergence starting on older islands and continuing on two simultaneous tracks towards younger islands. The amount of diversity generated through intra-island processes is skewed towards older islands, suggesting that island age significantly influences diversity. Phylogenetic concordance between nuclear and mitochondrial datasets and well-supported monophyletic species in mitochondrial derived topologies appear to reject the possibility of inter-species hybridization. These clear species boundaries might be related to the tight host associations of adult weevils in discrete ecological zones. If shared hosts facilitate hybridization, then host- or habitat-promoted divergences could prevent it, even in the case of species that share islands, since the altitudinal partitioning of habitats minimizes range overlap.

Colonization history, ecological shifts and diversification in the evolution of endemic Galápagos weevils

Mitochondrial DNA sequence data were obtained for eight species of flightless Galapaganus endemic weevils and one winged close relative in order to study their colonization history and modes of diversification in the Galápagos Archipelago. Contrary to most other insular radiations, the phylogeny estimates we recovered for Galapaganus do not follow the progression rule of island biogeography. The penalized likelihood age estimates of colonization of the archipelago exceed the age of the emerged islands and underscore the potential role of now sunken seamounts for the early evolution of Galapaganus. The phylogeny proposes one intra-island origin for Galapaganus endemics, but monophyly tests suggest a larger contribution of in-situ speciation on older islands. Generalist habitat preferences were reconstructed as ancestral while shifts to highland habitats were reconstructed as having evolved independently on different islands. Magnitudes and patterns of diversification rate were found to differ between older and younger islands. Our analyses reveal that the colonization sequence of islands and timing of colonization of Galapaganus could be linked with the geological and volcanic history of the islands in a rather complex scenario. Even though most islands appear to have been colonized soon after their emergence, there are notable deviations from the pattern of sequential colonization expected under the progression rule when considering only the extant emerged islands. Patterns of diversification rate variation on older and younger islands correspond to the volcanic activity or remnants of such activity, while the pattern of independent evolution of restricted habitat preferences in different islands suggests that habitat shifts could also have contributed to species diversity in Galapaganus.

Body size evolution simultaneously creates and collapses species boundaries in a clade of scincid lizards

Speciation is generally viewed as an irreversible process, although habitat alterations can erase reproductive barriers if divergence between ecologically differentiated species is recent. Reversed speciation might also occur if geographical contact is established between species that have evolved the same reproductive isolating barrier in parallel. Here, we demonstrate a loss of intrinsic reproductive isolation in a clade of scincid lizards as a result of parallel body size evolution, which has allowed for gene flow where large-bodied lineages are in secondary contact. An mtDNA phylogeny confirms the monophyly of the Plestiodon skiltonianus species complex, but rejects that of two size-differentiated ecomorphs. Mate compatibility experiments show that the high degree of body size divergence imposes a strong reproductive barrier between the two morphs; however, the strength of the barrier is greatly diminished between parallel-evolved forms. Since two large-bodied lineages are in geographical contact in the Sierra Nevada Mountains of California, we were also able to test for postzygotic isolation under natural conditions. Analyses of amplified fragment length polymorphisms show that extensive gene exchange is occurring across the contact zone, resulting in an overall pattern consistent with isolation by distance. These results provide evidence of reversed speciation between clades that diverged from a common ancestor more than 12Myr ago.

Adaptive radiation in Lesser Antillean lizards: molecular phylogenetics and species recognition in the Lesser Antillean dwarf gecko complex, Sphaerodactylus fantasticus

The time associated with speciation varies dramatically among lower vertebrates. The nature and timing of divergence is investigated in the fantastic dwarf gecko Sphaerodactylus fantasticus complex, a nominal species that occurs on the central Lesser Antillean island of Guadeloupe and adjacent islands and islets. This is compared to the divergence in the sympatric anole clade from the Anolis bimaculatus group. A molecular phylogenetic analysis of numerous gecko populations from across these islands, based on three mitochondrial DNA genes, reveals several monophyletic groups occupying distinct geographical areas, these being Les Saintes, western Basse Terre plus Dominica, eastern Basse Terre, Grand Terre, and the northern and eastern islands (Montserrat, Marie Galante, Petite Terre, Desirade). Although part of the same nominal species, the molecular divergence within this species complex is extraordinarily high (27% patristic distance between the most divergent lineages) and is compatible with this group occupying the region long before the origin of the younger island arc. Tests show that several quantitative morphological traits are correlated with the phylogeny, but in general the lineages are not uniquely defined by these traits. The dwarf geckos show notably less nominal species-level adaptive radiation than that found in the sympatric southern clade of Anolis bimculatus, although both appear to have occupied the region for a broadly similar period of time. Nevertheless, the dwarf gecko populations on Les Saintes islets are the most morphologically distinct and are recognized as a full species (Sphaerodactylus phyzacinus), as are anoles on Les Saintes (Anolis terraealtae).

Population structure, genetic variation and morphological diversity in indigenous sheep of Ethiopia

We investigated genetic and morphological diversity and population structure of 14 traditional sheep populations originating from four ecological zones in Ethiopia (sub-alpine, wet highland, sub-humid lowland and arid lowland). All animals (n = 672) were genotyped for 17 microsatellite markers and scored for 12 morphological characters. The sheep were initially classified as fat-tailed (11 populations), thin-tailed (one population) and fat-rumped sheep (two populations). These classifications are thought to correspond to three consecutive introduction events of sheep from the Near-East into East Africa. For the 14 populations, allelic richness ranged from 5.87 to 7.51 and expected heterozygosity (H(E)) from 0.66 to 0.75. Genetic differentiations (F(ST) values) between all pairs of populations, except between sub-alpine populations, were significantly different from zero (P < 0.001). Cluster analysis of morphological characters and a dendrogram constructed from genetic distances were broadly consistent with the classification into fat-tailed, thin-tailed and fat-rumped sheep. Bayesian cluster analysis using microsatellite markers indicated that there has been further genetic differentiation after the initial introduction of sheep into Ethiopia. Investigation of factors associated with genetic variation showed that an isolation-by-distance model, independently of other factors, explained most of the observed genetic variation. We also obtained a strong indication of adaptive divergence in morphological characters, patterns of morphological variation being highly associated with ecology even when the effect of neutral genetic divergence (F(ST)) was parcelled out in partial Mantel tests. Using a combination of F(ST) values, Bayesian clustering analysis and morphological divergence, we propose a classification of Ethiopian sheep into six breed groups and nine breeds.

B chromosome polymorphism in maize landraces: adaptive vs. demographic hypothesis of clinal variation

Cytogenetic analysis of maize landraces from northwestern Argentina has revealed an altitudinal cline in the mean number of B chromosomes (B's) per plant, with cultivars growing at higher altitudes exhibiting a higher number of B's. Altitudinal and longitudinal clines are frequently interpreted as evidence of selection, however, they can also be produced by the interplay between drift and spatially restricted gene flow or by admixture between previously isolated populations that have come into secondary contact. Here, we test the adaptive significance of the observed altitudinal gradient by comparing the levels of differentiation in the mean number of B's to those obtained from 18 selectively neutral loci [simple sequence repeats (SSRs)] among seven populations of the cline. The adequacy of alternative genetic-differentiation measures was determined, and associations between cytogenetic, genetic, and altitudinal distances were assessed by means of matrix- correspondence tests. No evidence for association between pairwise F(ST) and altitudinal distance or B-chromosome differentiation was found. The contrasting pattern of altitudinal divergence between the mean number of B's per plant and the genetic differentiation at SSR loci indicates that demographic processes cannot account for the observed levels of divergence in the mean number of B's.

Population structure and genetic diversity in tristylous Narcissus triandrus: insights from microsatellite and chloroplast DNA variation

We investigated cpDNA sequence and nuclear microsatellite variation among populations of the wild daffodil Narcissus triandrus to examine the role of historical vs. contemporary forces in shaping population structure, morphological differentiation and sexual-system evolution. This wide-ranging heterostylous species of the Iberian Peninsula is largely composed of two allopatric varieties (vars. cernuus and triandrus), and populations with either stylar trimorphism or dimorphism. Dimorphic populations only occur in var. triandrus, are mainly restricted to the northwestern portion of the species range, and uniformly lack the mid-styled morph (M-morph). Chloroplast DNA (cpDNA) sequence variation revealed strong geographical structuring and evidence for a fragmentation event associated with differentiation of the two varieties. In var. triandrus, population fragmentation, restricted gene flow and isolation-by-distance were also inferred. Significant differences in genetic diversity and population structure between the two varieties likely reflect historical and contemporary differences in demography and gene flow among populations. Discordance between cpDNA markers and both microsatellite and morphological variation indicate that hybridization has occurred between the two varieties at contact zones. There were no differences in genetic diversity or population structure between dimorphic and trimorphic populations, and chloroplast haplotypes were not associated with either sexual system, indicating transitions in morph structure within each maternal lineage. M-morph frequencies were positively correlated with differentiation at microsatellite loci, indicating that the evolutionary processes influencing these neutral markers also influence alleles controlling the style morphs.

Geological Dates and Molecular Rates: Rapid Divergence of Rivers and Their Biotas

We highlight a novel molecular clock calibration system based on geologically dated river reversal and river capture events. Changes in drainage pattern may effect vicariant isolation of freshwater taxa, and thus provide a predictive framework for associated phylogeographic study. As a case in point, New Zealand's Pelorus and Kaituna rivers became geologically isolated from the larger Wairau River system 70 to 130 kyr BP. We conducted mitochondrial DNA phylogeographic analyses of two unrelated freshwater-limited fish taxa native to these river systems (Gobiomorphus breviceps, n = 63; Galaxias divergens, n = 95). Phylogenetic analysis of combined control region and cytochrome b sequences yielded reciprocally monophyletic clades of Pelorus-Kaituna and Wairau haplotypes for each species. Calibrated rates of molecular change based on this freshwater vicariant event are substantially faster than traditionally accepted rates for fishes but consistent with other recent inferences based on geologically young calibration points. A survey of freshwater phylogeographic literature reveals numerous examples in which the ages of recent evolutionary events may have been substantially overestimated through the use of "accepted" calibrations. We recommend that--wherever possible--biologists should start to reassess the conclusions of such studies by using more appropriate molecular calibrations derived from recent geological events.

Geological history and within-island diversity: a debris avalanche and the Tenerife lizard Gallotia galloti

Several processes have been described that could explain geographical variation and speciation within small islands, including fragmentation of populations through volcanic eruptions. Massive landslides, or debris avalanches, could cause similar effects. Here we analyse the potential impact of the 0.8 million-year-ago (Ma) Güimar valley debris avalanche on the phylogeography of the lizard Gallotia galloti on the Canary Island of Tenerife. Distributions of mitochondrial DNA lineages (based on cytochrome b sequences) were analysed on a 60-km southeastern coast transect centred on this area. Three main clades were detected, which can be divided into northern (one clade) and southern (two clades) groups that introgress across the valley. Maximum-likelihood estimates of migration rates (scaled for mutation rate) revealed highly asymmetric patterns, indicating that long-term gene flow into this region from both the northern and the southern populations greatly exceeded that in the opposite directions, consistent with recolonization of the area. The ancestral Tenerife node on the G. galloti tree is estimated at 0.80 Ma, matching closely with the geological estimate for the debris avalanche. Morphological variation (body dimensions and scalation) was also analysed and indicated a stepped cline in female scalation across the valley, although the patterns for male scalation and male and female body dimensions were not as clear. Together these findings provide support for the hypothesis that the debris avalanche has shaped the phylogeography of G. galloti and may even have been a primary cause of the within-island cladogenesis through population fragmentation and isolation. Current estimates of timing of island unification mean that the original hypothesis that within-island diversity is explained by the secondary contact of populations from the two ancient precursor islands of Teno and Anaga is less plausible for this and some other Tenerife species. Large-scale landslides have occurred on many volcanic islands, and so may have been instrumental in shaping within-island diversities.

Phylogeography and the geographic cline in the armament of a seed-predatory weevil: effects of historical events vs. natural selection from the host plant

Japanese camellia (Camellia japonica) and its seed predator, the camellia weevil (Curculio camelliae), provide a notable example of a geographic mosaic of coevolution. In the species interaction, the offensive trait of the weevil (rostrum length) and the defensive trait of the plant (pericarp thickness) are involved in a geographically-structured arms race, and these traits and selective pressures acting on the plant defence vary greatly across a geographical landscape. To further explore the geographical structure of this interspecific interaction, we tested whether the geographical variation in the weevil rostrum over an 800-km range along latitude is attributed to local natural selection or constrained by historical (phylogeographical) events of local populations. Phylogeographical analyses of the mitochondrial DNA sequences of the camellia weevil revealed that this species has experienced differentiation into two regions, with a population bottleneck and subsequent range and/or population expansion within each region. Although these phylogeographical factors have affected the variation in rostrum length, analyses of competing factors for the geographical variation revealed that this pattern is primarily determined by the defensive trait of the host plant rather than by the effects of historical events of populations and a climatic factor (annual mean temperature). Thus, our study suggests the overwhelming strength of coevolutionary selection against the effect of historical events, which may have limited local adaptation.

Phenotypic variation in Galerida larks in Morocco: the role of history and natural selection

The relative influences of history, natural selection and hybridization in shaping phenotypic variation in closely related taxa is a crucial issue in current evolutionary biology. In this study, we used as a model two sibling but paradoxically highly variable species of larks (Galerida theklae and Galerida cristata) of Morocco to separate the impacts of these evolutionary forces. In the former species, variation is manifested mainly in colouration, while in the latter, variation also encompasses bill size and shape. Mitochondrial and nuclear DNA sequencing were used to identify the historical relationships among the subspecies and species. According to our analyses, G. cristata and G. theklae diverged about 3.7 million years ago (Ma), and we found no evidence for a role of hybridization in maintaining their similarity. In G. theklae, there was no further subdivision, while in G. cristata two major mtDNA groups were identified (divergence approximately 1.1 Ma). These two lineages are parapatric and regroup, respectively, the three short-billed subspecies [G. (cristata) cristata] and the two long-billed subspecies [G. (cristata) randonii]. Patterns of morphological variation were then contrasted to this pattern of neutral relationships: we found that G. (c.) cristata was morphologically more similar to G. theklae than to G. (c.) randonii. Overall, these results point towards the prominent role of (i) natural selection and/or phenotypic plasticity in adapting the plumage to local conditions and (ii) natural selection in combination with historical isolation in driving the divergence in size and bill morphology in the crested larks.

Molecular clocks and geological dates: cytochrome b of Anolis extremus substantially contradicts dating of Barbados emergence

Even though molecular clocks vary in rate to some extent, they are widely used and very important in a range of evolutionary studies, not least in interpreting cause and colonization in phylogeography. Evolutionists may use island age and emergence to give the earliest possible date for colonization by a species and hence give the lower limit in a molecular clock calibration. The geology of the Lesser Antilles is well studied and Barbados, although composed of some ancient rocks, is thought to have emerged only about 1 million years ago (Ma). The cytochrome b mitochondrial gene is the most widely used gene in vertebrate phylogeography, and generally evolves at a rate of 1-2% per million years (Myr) for poikilothermic vertebrates. Divergence measured across almost all of this gene in the endemic anole (Anolis extremus) reveals a mean patristic distance of approximately 8.3% between this clade and its sister, together with distinct divergence and phylogeographical structure within Barbados. The divergence time, estimated by a range of procedures using four calibration points, is not in the least compatible with the proposed geological time of emergence of Barbados. Hence, either the molecular clock rate does not apply to the Barbadian anole population, or the geological dating of the emergence of Barbados is erroneous. The compatibility of geological times and molecular divergence of this complex on Martinique, together with relative rates tests comparing the rates on Barbados and Martinique, do not suggest atypical clock rates. The question of whether Barbados emerged much earlier than is currently thought, or whether the molecular clock assumptions are inappropriate, remains open.

The dynamics of genetic and morphological variation on volcanic islands

Oceanic archipelagos of volcanic origin have been important in the study of evolution because they provide repeated natural experiments allowing rigorous tests of evolutionary hypotheses. Ongoing volcanism on these islands may, however, affect the evolutionary diversification of species. Analysis of population structure and phylogeographic patterns in island populations can provide insight into evolutionary dynamics on volcanic islands. We analysed genetic and morphological variation in the gecko Tarentola boettgeri on the island of Gran Canaria and compared it with Tarentola delalandii on Tenerife, a neighbouring volcanic island of similar age but distinctly different geological past. Intraspecific divergence of mitochondrial haplotypes indicates long-term persistence of Tarentola on each island, with a phylogeographic signal left by older volcanic events. More recent volcanic eruptions (approximately 0.2 million years ago on Tenerife, approximately 2.2 million years ago on Gran Canaria) have left a signature of population expansion in the population genetic structure, the strength of which depends on the time since the last major volcanic eruption on each island. While these stochastic events have left traces in morphological variation in Tenerife, in Gran Canaria geographical variation was solely associated with environmental variables. This suggests that historically caused patterns in morphology may be overwritten by natural selection within 2 million years.