Multilocus phylogeography of Australian teals (Anas spp.): a case study of the relationship between vagility and genetic structure

Biogeographic barriers potentially restrict gene flow but variation in dispersal or vagility can influence the effectiveness of these barriers among different species and produce characteristic patterns of population genetic structure. The objective of this study was to investigate interspecific and intraspecific genetic structure in two closely related species that differ in several life‐history characteristics. The grey teal Anas gracilis is geographically widespread throughout Australia with a distribution that crosses several recognized biogeographic barriers. This species has high vagility as its extensive movements track broad‐scale patterns in rainfall. In contrast, the closely related chestnut teal A. castanea is endemic to the mesic southeastern and southwestern regions of Australia and is more sedentary. We hypothesized that these differences in life‐history characteristics would result in more pronounced population structuring in the chestnut teal. We sequenced five nuclear loci (nuDNA) for 49 grey teal and 23 chestnut teal and compared results to published mitochondrial DNA (mtDNA) sequences. We used analysis of molecular variance to examine population structure, and applied coalescent based approaches to estimate demographic parameters. As predicted, chestnut teal were more strongly structured at both mtDNA and nuDNA (Φ_ST= 0.163 and 0.054, respectively) than were grey teal (Φ_ST < 0.0001 for both sets of loci). Surprisingly, a greater proportion of the total genetic variation was partitioned among populations within species (Φ_SC= 0.014 and 0.047 for nuDNA and mtDNA, respectively) than between the two species (Φ_CT < 0.0001 for both loci). The ‘Isolation with Migration’ coalescent model suggested a late Pleistocene divergence between the taxa, but remarkably, a deeper divergence between the southeastern and southwestern populations of chestnut teal. We conclude that dispersal potential played a prominent role in the structuring of populations within these species and that divergent selection associated with ecology and life history traits likely contributed to rapid and recent speciation in this pair.     

Speciation genomics and a role for the Z chromosome in the early stages of divergence between Mexican ducks and mallards

Speciation is a continuous and dynamic process, and studying organisms during the early stages of this process can aid in identifying speciation mechanisms. The mallard (Anas platyrhynchos) and Mexican duck (A. [p.] diazi) are two recently diverged taxa with a history of hybridization and controversial taxonomy. To understand their evolutionary history, we conducted genomic scans to characterize patterns of genetic diversity and divergence across the mitochondrial DNA (mtDNA) control region, 3523 autosomal loci and 172 Z‐linked sex chromosome loci. Between the two taxa, Z‐linked loci (Φ_ST = 0.088) were 5.2 times more differentiated than autosomal DNA (Φ_ST = 0.017) but comparable to mtDNA (Φ_ST = 0.092). This elevated Z differentiation deviated from neutral expectations inferred from simulated data that incorporated demographic history and differences in effective population sizes between marker types. Furthermore, 3% of Z‐linked loci, compared to <0.1% of autosomal loci, were detected as outlier loci under divergent selection with elevated relative (Φ_ST) and absolute (d_XY) estimates of divergence. In contrast, the ratio of Z‐linked and autosomal differentiation among the seven Mexican duck sampling locations was close to 1:1 (Φ_ST = 0.018 for both markers). We conclude that between mallards and Mexican ducks, divergence at autosomal markers is largely neutral, whereas greater divergence on the Z chromosome (or some portions thereof) is likely the product of selection that has been important in speciation. Our results contribute to a growing body of literature indicating elevated divergence on the Z chromosome and its likely importance in avian speciation.     

Hybrid speciation in sparrows II: a role for sex chromosomes?

Homoploid hybrid speciation in animals is poorly understood, mainly because of the scarcity of well‐documented cases. Here, we present the results of a multilocus sequence analysis on the house sparrow (Passer domesticus), Spanish sparrow (P. hispaniolensis) and their proposed hybrid descendant, the Italian sparrow (P. italiae). The Italian sparrow is shown to be genetically intermediate between the house sparrow and Spanish sparrow, exhibiting genealogical discordance and a mosaic pattern of alleles derived from either of the putative parental species. The average variation on the Z chromosome was significantly reduced compared with autosomal variation in the putative parental species, the house sparrow and Spanish sparrow. Additionally, divergence between the two species was elevated on the Z chromosome relative to the autosomes. This pattern of variation and divergence is consistent with reduced introgression of Z‐linked genes and/or a faster‐Z effect (increased rate of adaptive divergence on the Z). F_ST‐outlier tests were consistent with the faster‐Z hypothesis: two of five Z‐linked loci (CHD1Z and PLAA) were identified as candidates for being subject to positive, divergent selection in the putative parental species. Interestingly, the two latter genes showed a mosaic pattern in the (hybrid) Italian sparrow; that is, the Italian sparrow was found to be fixed for Spanish sparrow alleles at CHD1Z and to mainly have house sparrow alleles at PLAA. Preliminary evidence presented in this study thus suggests that sex chromosomes may play a significant role in this case of homoploid hybrid speciation.     

GENETIC AND PHENOTYPIC DIVERGENCE BETWEEN LOW‐ AND HIGH‐ALTITUDE POPULATIONS OF TWO RECENTLY DIVERGED CINNAMON TEAL SUBSPECIES

Spatial variation in the environment can lead to divergent selection between populations occupying different parts of a species’ range, and ultimately lead to population divergence. The colonization of new areas can thus facilitate divergence in beneficial traits, yet with little differentiation at neutral genetic markers. We investigated genetic and phenotypic patterns of divergence between low‐ and high‐altitude populations of cinnamon teal inhabiting normoxic and hypoxic regions in the Andes and adjacent lowlands of South America. Cinnamon teal showed strong divergence in body size (PC1; P_ST= 0.56) and exhibited significant frequency differences in a single nonsynonymous α‐hemoglobin amino acid polymorphism (Asn/Ser‐α9; F_ST= 0.60) between environmental extremes, despite considerable admixture of mtDNA and intron loci (F_ST= 0.004–0.168). Inferences of strong population segregation were further supported by the observation of few mismatched individuals in either environmental extreme. Coalescent analyses indicated that the highlands were most likely colonized from lowland regions but following divergence, gene flow has been asymmetric from the highlands into the lowlands. Multiple selection pressures associated with high‐altitude habitats, including cold and hypoxia, have likely shaped morphological and genetic divergence within South American cinnamon teal populations.     

Genomic footprints of speciation in Atlantic eels (Anguilla anguilla and A. rostrata)

The importance of speciation‐with‐geneflow scenarios is increasingly appreciated. However, the specific processes and the resulting genomic footprints of selection are subject to much discussion. We studied the genomics of speciation between the two panmictic, sympatrically spawning sister species; European (Anguilla anguilla) and American eel (A. rostrata). Divergence is assumed to have initiated more than 3 Ma, and although low gene flow still occurs, strong postzygotic barriers are present. Restriction‐site‐associated DNA (RAD) sequencing identified 328 300 SNPs for subsequent analysis. However, despite the presence of 3757 strongly differentiated SNPs (F_ST > 0.8), sliding window analyses of F_ST showed no larger genomic regions (i.e. hundreds of thousands to millions of bases) of elevated differentiation. Overall F_ST was 0.041, and linkage disequilibrium was virtually absent for SNPs separated by more than 1000 bp. We suggest this to reflect a case of genomic hitchhiking, where multiple regions are under directional selection between the species. However, low but biologically significant gene flow and high effective population sizes leading to very low genetic drift preclude accumulation of strong background differentiation. Genes containing candidate SNPs for positive selection showed significant enrichment for gene ontology (GO) terms relating to developmental processes and phosphorylation, which seems consistent with assumptions that differences in larval phase duration and migratory distances underlie speciation. Most SNPs under putative selection were found outside coding regions, lending support to emerging views that noncoding regions may be more functionally important than previously assumed. In total, the results demonstrate the necessity of interpreting genomic footprints of selection in the context of demographic parameters and life‐history features of the studied species.     

Low population genetic differentiation in the Orchidaceae: implications for the diversification of the family

A leading hypothesis for the immense diversity of the Orchidaceae is that skewed mating success and small, disjunct populations lead to strong genetic drift and switches between adaptive peaks. This mechanism is only possible under conditions of low gene flow that lead to high genetic differentiation among populations. We tested whether orchids typically exhibit high levels of population genetic differentiation by conducting a meta‐analysis to compare mean levels of population genetic differentiation (F_ST) between orchids and other diverse families and between rare and common orchids. Compared with other families, the Orchidaceae is typically characterized by relatively low genetic differentiation among populations (mean F_ST = 0.146) at allozyme loci. Rare terrestrial orchids showed higher population genetic differentiation than common orchids, although this value was still lower than the mean for most plant families. All lines of evidence suggest that orchids are typically characterized by low levels of population genetic differentiation, even in species with naturally disjunct populations. As such, we found no strong evidence that genetic drift in isolated populations has played a major role in the diversification of the Orchidaceae. Further research into the diversification of the family needs to unravel the relative roles of biotic and environmental selective pressures in the speciation of orchids.     

Genome divergence and diversification within a geographic mosaic of coevolution

Despite substantial interest in coevolution's role in diversification, examples of coevolution contributing to speciation have been elusive. Here, we build upon past studies that have shown both coevolution between South Hills crossbills and lodgepole pine (Pinus contorta), and high levels of reproductive isolation between South Hills crossbills and other ecotypes in the North American red crossbill (Loxia curvirostra) complex. We used genotyping by sequencing to generate population genomic data and applied phylogenetic and population genetic analyses to characterize the genetic structure within and among nine of the ecotypes. Although genome‐wide divergence was slight between ecotypes (F_ST = 0.011–0.035), we found evidence of relative genetic differentiation (as measured by F_ST) between and genetic cohesiveness within many of them. As expected for nomadic and opportunistic breeders, we detected no evidence of isolation by distance. The one sedentary ecotype, the South Hills crossbill, was genetically most distinct because of elevated divergence at a small number of loci rather than pronounced overall genome‐wide divergence. These findings suggest that mechanisms related to recent local coevolution between South Hills crossbills and lodgepole pine (e.g. strong resource‐based density dependence limiting gene flow) have been associated with genome divergence in the face of gene flow. Our results further characterize a striking example of coevolution driving speciation within perhaps as little as 6000 years.     

Population structure of island‐associated dolphins: Evidence from mitochondrial and microsatellite markers for common bottlenose dolphins (Tursiops truncatus) around the main Hawaiian Islands

We used mitochondrial and nuclear genetic markers to investigate population structure of common bottlenose dolphins, Tursiops truncatus, around the main Hawaiian Islands. Though broadly distributed throughout the world's oceans, bottlenose dolphins are known to form small populations in coastal waters. Recent photo‐identification data suggest the same is true in Hawaiian waters. We found genetic differentiation among (mtDNA Φ_ST= 0.014–0.141, microsatellite F’_ST= 0.019–0.050) and low dispersal rates between (0.17–5.77 dispersers per generation) the main Hawaiian Island groups. Our results are consistent with movement rates estimated from photo‐identification data and suggest that each island group supports a demographically independent population. Inclusion in our analyses of samples collected near Palmyra Atoll provided evidence that the Hawaiian Islands are also occasionally visited by members of a genetically distinct, pelagic population. Two of our samples exhibited evidence of partial ancestry from Indo‐Pacific bottlenose dolphins (T. aduncus), a species not known to inhabit the Hawaiian Archipelago. Our findings have important implications for the management of Hawaiian bottlenose dolphins and raise concerns about the vulnerability to human impacts of pelagic species in island ecosystems.     

Historic and recent fragmentation coupled with altitude affect the genetic population structure of one of the world's highest tropical tree line species

Aim To assess the effects of altitude and historic and recent forest fragmentation on the genetic diversity and structure of the wind‐pollinated tropical tree line species Polylepis incana. Location One of the highest mountain forest regions of the world, located in the Eastern Cordillera of the Ecuadorian Andes. Methods We compared genetic diversity and structure of adult trees with those of seedlings (n= 118 in both cases) in nine forest stands spanning an altitudinal gradient from 3500 to 4100 m a.s.l. using amplified fragment length polymorphisms (AFLPs). Genetic diversity was calculated as percentage of polymorphic bands (P) and Nei's expected heterozygosity (He); genetic differentiation was assessed using analysis of molecular variance, Φ_ST statistics and Bayesian cluster analysis. Results Estimates of genetic diversity at the population level were significantly lower in seedlings than in adults. Genetic diversity (He‐value) was, in both cases, negatively correlated to altitude and positively correlated to population size in the seedlings. Genetic differentiation of the seedlings was approximately as high (φ_ST= 0.298) as that of the adults (φ_ST= 0.307), and geographical differentiation was clearly reflected in both AFLP profiles, with mountain ridges acting as barriers to gene flow. Main conclusions Our study provides evidence of a historic upslope migration of P. incana in central Ecuador. In addition, it highlights the detrimental effects of unexpectedly strong genetic isolation, both recent and historical, particularly for our wind‐pollinated species where the distance between forest stands was less than 25 km. We therefore additionally propose that in habitats with pronounced high‐mountain landscape structures, gene flow may be hampered to such an extent that species have a more pronounced sensitivity to habitat fragmentation, even among populations of wind‐pollinated trees.     

Evidence for nonallopatric speciation among closely related sympatric Heliotropium species in the Atacama Desert

The genetic structure of populations of closely related, sympatric species may hold the signature of the geographical mode of the speciation process. In fully allopatric speciation, it is expected that genetic differentiation between species is homogeneously distributed across the genome. In nonallopatric speciation, the genomes may remain undifferentiated to a large extent. In this article, we analyzed the genetic structure of five sympatric species from the plant genus Heliotropium in the Atacama Desert. We used amplified fragment length polymorphisms (AFLPs) to characterize the genetic structure of these species and evaluate their genetic differentiation as well as the number of loci subject to positive selection using divergence outlier analysis (DOA). The five species form distinguishable groups in the genetic space, with zones of overlap, indicating that they are possibly not completely isolated. Among‐species differentiation accounts for 35% of the total genetic differentiation (F_ST = 0.35), and F_ST between species pairs is positively correlated with phylogenetic distance. DOA suggests that few loci are subject to positive selection, which is in line with a scenario of nonallopatric speciation. These results support the idea that sympatric species of Heliotropium sect. Cochranea are under an ongoing speciation process, characterized by a fluctuation of population ranges in response to pulses of arid and humid periods during Quaternary times.     

Global patterns of population genetic differentiation in seed plants

Evaluating the factors that drive patterns of population differentiation in plants is critical for understanding several biological processes such as local adaptation and incipient speciation. Previous studies have given conflicting results regarding the significance of pollination mode, seed dispersal mode, mating system, growth form and latitudinal region in shaping patterns of genetic structure, as estimated by F_ST values, and no study to date has tested their relative importance together across a broad scale. Here, we assembled a 337‐species data set for seed plants from publications with data on F_ST from nuclear markers and species traits, including variables pertaining to the sampling scheme of each study. We used species traits, while accounting for sampling variables, to perform phylogenetic multiple regressions. Results demonstrated that F_ST values were higher for tropical, mixed‐mating, non‐woody species pollinated by small insects, indicating greater population differentiation, and lower for temperate, outcrossing trees pollinated by wind. Among the factors we tested, latitudinal region explained the largest portion of variance, followed by pollination mode, mating system and growth form, while seed dispersal mode did not significantly relate to F_ST. Our analyses provide the most robust and comprehensive evaluation to date of the main ecological factors predicted to drive population differentiation in seed plants, with important implications for understanding the basis of their genetic divergence. Our study supports previous findings showing greater population differentiation in tropical regions and is the first that we are aware of to robustly demonstrate greater population differentiation in species pollinated by small insects.     

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.     

Fine-scale genetic population structure in a mobile marine mammal: inshore bottlenose dolphins in Moreton Bay, Australia

Highly mobile marine species in areas with no obvious geographic barriers are expected to show low levels of genetic differentiation. However, small‐scale variation in habitat may lead to resource polymorphisms and drive local differentiation by adaptive divergence. Using nuclear microsatellite genotyping at 20 loci, and mitochondrial control region sequencing, we investigated fine‐scale population structuring of inshore bottlenose dolphins (Tursiops aduncus) inhabiting a range of habitats in and around Moreton Bay, Australia. Bayesian structure analysis identified two genetic clusters within Moreton Bay, with evidence of admixture between them (F_ST = 0.05, P = 0.001). There was only weak isolation by distance but one cluster of dolphins was more likely to be found in shallow southern areas and the other in the deeper waters of the central northern bay. In further analysis removing admixed individuals, southern dolphins appeared genetically restricted with lower levels of variation (AR = 3.252, π = 0.003) and high mean relatedness (r = 0.239) between individuals. In contrast, northern dolphins were more diverse (AR = 4.850, π = 0.009) and were mixing with a group of dolphins outside the bay (microsatellite‐based STRUCTURE analysis), which appears to have historically been distinct from the bay dolphins (mtDNA Φ_ST = 0.272, P < 0.001). This study demonstrates the ability of genetic techniques to expose fine‐scale patterns of population structure and explore their origins and mechanisms. A complex variety of inter‐related factors including local habitat variation, differential resource use, social behaviour and learning, and anthropogenic disturbances are likely to have played a role in driving fine‐scale population structure among bottlenose dolphins in Moreton Bay.     

ddRAD‐seq data reveal significant genome‐wide population structure and divergent genomic regions that distinguish the mallard and close relatives in North America

Recently evolved species typically share genetic variation across their genomes due to incomplete lineage sorting and/or ongoing gene flow. Given only subtle allele frequency differences at most loci and the expectation that divergent selection may affect only a tiny fraction of the genome, distinguishing closely related species based on multi‐locus data requires substantial genomic coverage. In this study, we used ddRAD‐seq to sample the genomes of five recently diverged, New World “mallards” (Anas spp.), a group of dabbling duck species characterized by diagnosable phenotypic differences but minimal genetic differentiation. With increased genomic sampling, we aimed to characterize population structure within this group and identify genomic regions that may have experienced divergent selection during speciation. We analyzed 3,017 autosomal ddRAD‐seq loci and 177 loci from the Z‐chromosome. In contrast to previous studies, the ddRAD‐seq data were sufficient to assign individuals to their respective species or subspecies and to generate estimates of gene flow in a phylogenetic framework. We find limited evidence of contemporary gene flow between the dichromatic mallard and several monochromatic taxa, but find evidence for historical gene flow between some monochromatic species pairs. We conclude that the overall genetic similarity of these taxa likely reflects retained ancestral polymorphism rather than recent and extensive gene flow. Thus, despite recurring cases of hybridization in this group, our results challenge the current dogma predicting the genetic extinction of the New World monochromatic dabbling ducks via introgressive hybridization with mallards. Moreover, ddRAD‐seq data were sufficient to identify previously unknown outlier regions across the Z‐chromosome and several autosomal chromosomes that may have been involved in the diversification of species in this recent radiation.     

Evidence for adaptive phenotypic differentiation in Baltic Sea sticklebacks

The evidence for adaptive phenotypic differentiation in mobile marine species remains scarce, partly due to the difficulty of obtaining quantitative genetic data to demonstrate the genetic basis of the observed phenotypic differentiation. Using a combination of phenotypic and molecular genetic approaches, we elucidated the relative roles of natural selection and genetic drift in explaining lateral plate number differentiation in threespine sticklebacks (Gasterosteus aculeatus) across the entire Baltic Sea basin (approximately 392 000 km²). We found that phenotypic differentiation (P_ST = 0.213) in plate number exceeded that in neutral markers (F_ST = 0.008), suggesting an adaptive basis for the observed differentiation. Because a close correspondence was found between plate phenotype and genotype at a quantitative trait loci (QTL; STN381) tightly linked to the gene (Ectodysplasin) underlying plate variation, the evidence for adaptive differentiation was confirmed by comparison of F_ST at the QTL (F_STQ = 0.089) with F_ST at neutral marker loci. Hence, the results provide a comprehensive demonstration of adaptive phenotypic differentiation in a high‐gene‐flow marine environment with direct, rather than inferred, verification for the genetic basis of this differentiation. In general, the results illustrate the utility of P_ST–F_ST–F_STQ comparisons in uncovering footprints of natural selection and evolution and add to the growing evidence for adaptive genetic differentiation in high‐gene‐flow marine environments, including that of the relatively young Baltic Sea.     

EXAMINING TWO STANDARD ASSUMPTIONS OF ANCESTRAL RECONSTRUCTIONS: REPEATED LOSS OF DICHROMATISM IN DABBLING DUCKS (ANATINI)

Although phylogenetic reconstruction of ancestral character states is becoming an increasingly common technique for studying evolution, few researchers have assessed the reliability of these reconstructions. Here I test for congruence between a phylogenetic reconstruction and a widely accepted scenario based on independent lines of evidence. I used Livezey's (1991) phylogeny to reconstruct ancestral states of plumage dichromatism in dabbling ducks (Anatini). Character state mapping reconstructs monochromatic ancestors for the genus Anas as well as most of its main clades. This reconstruction differs strongly from the widely accepted scenario of speciation and plumage evolution in the group (e.g., Delacour and Mayr 1945; Sibley 1957). This incongruence may occur because two standard assumptions of character state reconstruction are probably not met in this case. Violating either of these two assumptions would be a source of error sufficient to create misleading reconstructions. The first assumption that probably does not apply to ducks is that terminal taxa, in this case species, are monophyletic. Many of the widespread dichromatic species of ducks may be paraphyletic and ancestral to isolated monochromatic species. Three lines of evidence support this scenario: population-level phylogenies, biogeography, and vestigial plumage patterns. The second assumption that probably does not apply to duck plumage color is that gains and losses of character states are equally likely. Four lines of evidence suggest that dichromatic plumage might be lost more easily than gained: weak female preferences for bright male plumage, biases toward the loss of sexually dichromatic characters, biases toward the loss of complex characters, and repeated loss of dichromatism in other groups of birds. These seven lines of evidence support the accepted scenario that widespread dichromatic species repeatedly budded off isolated monochromatic species. Drift and genetic biases probably caused the easy loss of dichromatism in ducks and other birds during peripatric speciation. In order to recover the accepted scenario using Livezey's tree, losses of dichromatism must be five times more likely than gains. The results of this study caution against the uncritical use of unordered parsimony as the sole criterion for inferring ancestral states. Detailed population-level sampling is needed and altered transformation weighting may be warranted in ducks and in many other groups and character types with similar attributes.     

Adaptive divergence despite strong genetic drift: genomic analysis of the evolutionary mechanisms causing genetic differentiation in the island fox (Urocyon littoralis)

The evolutionary mechanisms generating the tremendous biodiversity of islands have long fascinated evolutionary biologists. Genetic drift and divergent selection are predicted to be strong on islands and both could drive population divergence and speciation. Alternatively, strong genetic drift may preclude adaptation. We conducted a genomic analysis to test the roles of genetic drift and divergent selection in causing genetic differentiation among populations of the island fox (Urocyon littoralis). This species consists of six subspecies, each of which occupies a different California Channel Island. Analysis of 5293 SNP loci generated using Restriction‐site Associated DNA (RAD) sequencing found support for genetic drift as the dominant evolutionary mechanism driving population divergence among island fox populations. In particular, populations had exceptionally low genetic variation, small N_e (range = 2.1–89.7; median = 19.4), and significant genetic signatures of bottlenecks. Moreover, islands with the lowest genetic variation (and, by inference, the strongest historical genetic drift) were most genetically differentiated from mainland grey foxes, and vice versa, indicating genetic drift drives genome‐wide divergence. Nonetheless, outlier tests identified 3.6–6.6% of loci as high F_ST outliers, suggesting that despite strong genetic drift, divergent selection contributes to population divergence. Patterns of similarity among populations based on high F_ST outliers mirrored patterns based on morphology, providing additional evidence that outliers reflect adaptive divergence. Extremely low genetic variation and small N_e in some island fox populations, particularly on San Nicolas Island, suggest that they may be vulnerable to fixation of deleterious alleles, decreased fitness and reduced adaptive potential.     

Refugia within refugia: the case study of a canopy tree (Eurycorymbus cavaleriei) in subtropical China

Aim Eurycorymbus cavaleriei (Lévl.) Rehd. et Hand.‐Mazz. (Sapindaceae) is a Tertiary relict tree endemic to subtropical China. This area is a centre for speciation and evolution within the East Asia biome and one of its most important refugial locations. In this study we aim to elucidate the phylogeographical patterning in E. cavaleriei, in order to identify the locations of the species’ main refugia and the predominant patterns of migration that have led to the contemporary spatial genetic structure of chloroplast variation. Location Subtropical China. Methods We sampled 18 populations of E. cavaleriei throughout its geographical range. Chloroplast DNA (cpDNA) sequence data from two non‐coding regions ((trnS/trnG and pl20/5′‐rps12) were obtained from 170 individuals for phylogeographical analyses. Relationships among cpDNA haplotypes were determined using median‐joining networks. Genetic structure was examined by spatial analysis of molecular variance (SAMOVA). Population differentiation was estimated by G_ST and N_ST statistics. Results Ten distinct haplotypes were identified. The level of differentiation among populations was relatively high (G_ST = 0.817), and N_ST was significantly higher than G_ST (P < 0.05), indicating that strong phylogeographical structure is exhibited by this species. The SAMOVA revealed five diverging groups of related haplotypes, which coincide with major landscape features in this region. Main conclusions The high differentiation among populations of E. cavaleriei may be a combined effect of historical and contemporary processes, such as the low effective population size for the chloroplast genome of a dioecious species, long‐term range fragmentation and limited seed dispersal for the species. Clear‐cut geographical distributions of ancestral haplotypes of the species suggest multiple potential refugia across subtropical China. The identified refugial regions have long been recognized as centres of plant diversity and endemism for China and have also been suggested as glacial refugia for many other plant species. The combination of these factors means that these locations should be considered as the highest priority for inclusion in conservation policies and sustainable forest management strategies for subtropical China.     

Multiple origins of sexual dichromatism and aposematism within large carpenter bees

The evolution of reversed sexual dichromatism and aposematic coloration has long been of interest to both theoreticians and empiricists. Yet despite the potential connections between these phenomena, they have seldom been jointly studied. Large carpenter bees (genus Xylocopa) are a promising group for such comparative investigations as they are a diverse clade in which both aposematism and reversed sexual dichromatism can occur either together or separately. We investigated the evolutionary history of dichromatism and aposematism and a potential correlation of these traits with diversification rates within Xylocopa, using a newly generated phylogeny for 179 Xylocopa species based on ultraconserved elements (UCEs). A monochromatic, inconspicuous ancestor is indicated for the genus, with subsequent convergent evolution of sexual dichromatism and aposematism in multiple lineages. Aposematism is found to covary with reversed sexual dichromatism in many species; however, reversed dichromatism also evolved in non‐aposematic species. Bayesian Analysis of Macroevolutionary Models (BAMM) did not show increased diversification in any specific clade in Xylocopa, whereas support from Hidden State Speciation and Extinction (HiSSE) models remained inconclusive regarding an association of increased diversification rates with dichromatism or aposematism. We discuss the evolution of color patterns and diversification in Xylocopa by considering potential drivers of dichromatism and aposematism.     

Morphological stasis in the Eurytemora affinis species complex (Copepoda: Temoridae)

Morphological stasis has long been regarded as one of the most challenging problems in evolutionary biology. This study focused on the copepod species complex, Eurytemora affinis, as a model system to determine pattern and degree of morphological stasis. This study revealed discordant rates of morphological differentiation, molecular evolution, and reproductive isolation, where speciation was accompanied by lack of morphological differentiation in secondary sex characters. Comparisons were made among phylogenies based on morphometrics, nuclear (allozyme) loci, and mitochondrial DNA (mtDNA) sequences from cytochrome oxidase I, for a total of 43 populations within the complex. These systematic relationships were also compared to patterns of reproductive isolation. In addition, genetic subdivision of nuclear molecular (allozyme) markers (G _ST) and quantitative (morphological) characters (Q _ST) were determined to infer evolutionary forces driving morphological differentiation. The morphometric phylogeny revealed that all clades, excluding the European clade, were morphologically undifferentiated and formed a polytomy (multifurcation). Morphometric distances were not correlated with mtDNA distances, or with patterns of reproductive isolation. In contrast, nuclear and mtDNA phylogenies were mostly congruent. Reproductive isolation proved to be the most sensitive indicator of speciation, given that two genetically and morphologically proximate populations showed evidence of hybrid breakdown. Quantitative genetic (morphological) subdivision (Q _ST = 0.162) was lower than nuclear genetic subdivision (G _ST = 0.617) for four laboratory-reared North American populations, indicating retarded evolution of morphological characters. This result contrasts with most other species, where Q _ST typically exceeds G _ST as a result of directional selection. Thus, in all but the European populations, evolution of the secondary sex characters was marked by morphological stasis, even between reproductively-isolated populations.