Speciation in the White-breasted Nuthatch (Sitta carolinensis): a multilocus perspective

Inferring the evolutionary and ecological processes that have shaped contemporary species distributions using the geographic distribution of gene lineages is the principal goal of phylogeographic research. Researchers in the field have recognized that inferences made from a single gene, often mitochondrial, can be informative regarding the pattern of diversification but lack conclusive information regarding the evolutionary mechanisms that led to the observed patterns. Here, we use a multilocus (20 loci) data set to explore the evolutionary history of the White‐breasted Nuthatch (Sitta carolinensis). A previous single‐locus study found S. carolinensis is comprised of four reciprocally monophyletic clades geographically restricted to the pine and oak forests of: (i) eastern North America, (ii) southern Rocky Mountain and Mexican Mountain ranges, (iii) Eastern Sierra Nevada and Northern Rocky Mountains and (iv) Pacific slope of North America. The diversification of the clades was attributed to the fragmentation of North American pine and oak woodlands in the Pliocene with subsequent divergences owing to the Pleistocene glacial cycles. Principal component, clustering and species tree analyses of the multilocus data resolved the same four groups or lineages found in the single‐locus study. Coalescent analyses and hypothesis testing of nested isolation and migration models indicate that isolation and not gene flow has been the major evolutionary mechanism responsible for shaping genetic variation, and all the divergence events within S. carolinensis have occurred in response to the Pleistocene glacial cycles.     

The anatomy of a ‘suture zone’ in Amazonian butterflies: a coalescent‐based test for vicariant geographic divergence and speciation

Attempts by biogeographers to understand biotic diversification in the Amazon have often employed contemporary species distribution patterns to support particular theories, such as Pleistocene rainforest refugia, rather than to test among alternative hypotheses. Suture zones, narrow regions where multiple contact zones and hybrid zones between taxa cluster, have been seen as evidence for past expansion of whole biotas that have undergone allopatric divergence in vicariant refuges. We used coalescent analysis of mutilocus sequence data to examine population split times in 22 pairs of geminate taxa in ithomiine and heliconiine butterflies. We test a hypothesis of simultaneous divergence across a suture zone in NE Peru. Our results reveal a scattered time course of diversification in this suture zone, rather than a tight cluster of split times. Additionally, we find rapid diversification within some lineages such as Melinaea contrasting with older divergence within lineages such as the Oleriina (Hyposcada and Oleria). These results strongly reject simple vicariance as a cause of the suture zone. At the same time, observed lineage effects are incompatible with a series of geographically coincident vicariant events which should affect all lineages similarly. Our results suggest that Pleistocene climatic forcing cannot readily explain this Peruvian suture zone. Lineage‐specific biological traits, such as characteristic distances of gene flow or varying rates of parapatric divergence, may be of greater importance.     

Cryptic diversity in a widespread North American songbird: phylogeography of the Brown Creeper (Certhia americana)

The identification of species via morphological characteristics has traditionally left cryptic species undescribed in taxa under selection for morphological conservation (or a lack of selection for morphological change). Treecreepers (Genus: Certhia) have a conserved morphological appearance, making it difficult to ascertain relationships in the genus based on morphology alone. Recent genetic and song structure studies of Eurasian Treecreepers identified cryptic species within Old World Certhia that were previously undescribed using morphological characteristics. Here, we use mtDNA to investigate cryptic diversity and patterns of diversification in the Brown Creeper (Certhia americana), the single described Certhia species in the Americas. Phylogenetic analyses identified six well-supported geographically-structured clades; the basal divergence separates a northern and a southern lineage in the Brown Creeper, likely cryptic species previously characterized as many subspecies. Sympatry is prevalent between clades in western North America, where possible contact zones warrant further investigation. Allopatry appears to be the primary driver of deep phylogeographic structure within the Brown Creeper; however, within clade diversity is highly correlated with the life history traits of the populations that comprise the geographically structured phylogroups.     

Host‐driven divergence in the parasitic plant Orobanche minor Sm. (Orobanchaceae)

Many parasitic angiosperms have a broad host range and are therefore considered to be host generalists. Orobanche minor is a nonphotosynthetic root parasite that attacks a range of hosts from taxonomically disparate families. In the present study, we show that O. minor sensu lato may comprise distinct, genetically divergent races isolated by the different ecologies of their hosts. Using a three‐pronged approach, we tested the hypothesis that intraspecific taxa O. minor var. minor and O. minor ssp. maritima parasitizing either clover (Trifolium pratense) or sea carrot (Daucus carota ssp. gummifer), respectively, are in allopatric isolation. Morphometric analysis revealed evidence of divergence but this was insufficient to define discrete, host‐specific taxa. Intersimple sequence repeat (ISSR) marker‐based data provided stronger evidence of divergence, suggesting that populations were isolated from gene flow. Phylogenetic analysis, using sequence‐characterized amplified region (SCAR) markers derived from ISSR loci, provided strong evidence for divergence by clearly differentiating sea carrot‐specific clades and mixed‐host clades. Low levels of intrapopulation SCAR marker sequence variation and floral morphology suggest that populations on different hosts are probably selfing and inbreeding. Morphologically cryptic Orobanche taxa may therefore be isolated from gene flow by host ecology. Together, these data suggest that host specificity may be an important driver of allopatric speciation in parasitic plants.     

Low genetic diversity and local adaptive divergence of Dracaena cambodiana (Liliaceae) populations associated with historical population bottlenecks and natural selection: an endangered long‐lived tree endemic to Hainan Island,China

Historical population bottlenecks and natural selection have important effects on the current genetic diversity and structure of long‐lived trees. Dracaena cambodiana is an endangered, long‐lived tree endemic to Hainan Island, China. Our field investigations showed that only 10 populations remain on Hainan Island and that almost all have been seriously isolated and grow in distinct habitats. A considerable amount of genetic variation at the species level, but little variation at the population level, and a high level of genetic differentiation among the populations with limited gene flow in D. cambodiana were detected using inter‐simple sequence repeat (ISSR) and random amplified polymorphic DNA (RAPD) analyses. No significant correlation was found between genetic diversity and actual population size, as the genetic diversities were similar regardless of population size. The Mantel test revealed that there was no correlation between genetic and geographic distances among the 10 populations. The UPGMA, PCoA and Bayesian analyses showed that local adaptive divergence has occurred among the D. cambodiana populations, which was further supported by habitat‐private fragments. We suggest that the current genetic diversity and population differentiation of D. cambodiana resulted from historical population bottlenecks and natural selection followed by historical isolation. However, the lack of natural regeneration of D. cambodiana indicates that former local adaptations with low genetic diversity may have been genetically weak and are unable to adapt to the current ecological environments.     

Molecular and morphological divergence in the butterfly genus Lycaeides (Lepidoptera: Lycaenidae) in North America: evidence of recent speciation

Male genital morphology, allozyme allele frequencies and mtDNA sequence variation were surveyed in the butterfly species Lycaeides idas and L. melissa from across much of their range in North America. Despite clear differences in male genital morphology, wing colour patterns and habitat characteristics, genetic variation was not taxonomically or geographically structured and the species were not identifiable by either genetic data set. Genetic distances (Nei's D=0.002–0.078, calculated from allozyme data) between all populations of both species were within the range commonly observed for conspecific populations of other butterflies. The most frequent mtDNA haplotype was present in individuals of both species in populations from southern California to Wisconsin. We conclude that speciation has probably happened recently and the lack of genetic differentiation between the species is the product of either (1) recent or ongoing gene flow at neutral loci, and/or (2) an insufficiency of time for lineage sorting. The evolution of male genital morphology, wing colour patterns and ecological characteristics has proceeded more rapidly than allozyme or mtDNA evolution.     

Generic species richness and body mass in North American mammals: Support for the inverse relationship of body size and speciation rate

Generic species richness, the number of species per genus, is examined as a function of mean generic body mass for extant North American mammals. Species richness decreases as an inverse power function with increased mass, and the Spearman rank correlation coefficient of the logio transformed data is significant (r_s= ‐0.37). When the data are partitioned by trophic level, the relationship is not statistically significant for carnivores but strengthens for herbivores (r_s= ‐0.46). This interesting but incidental effect is due to the negligible number of diminutive and excessively large carnivores, which is in turn determined by foraging strategies. Alternate hypotheses for the “right‐skewed”; size distribution of modern North American mammals, such as disproportionate extinction of large species, differential species longevity, and a geographical scaling function, are rejected in favor of the proposition that elevated levels of speciation are restricted to animals of small body mass, as originally proposed by Gould and Eldredge (1977). This phenomenon is explained as a function of habitat restriction and particularly in herbivores, limited home range size. Aquatic mammals, regardless of body size, speciate rarely. Cope's Rule, the tendency of many animal groups to evolve towards large size, is understood as a probabilistic statement reflecting the phylogenetic tendencies of a disproportionately high number of small species alive at any given point in time.     

Detecting cryptic speciation in the widespread and morphologically conservative carpet chameleon (Furcifer lateralis) of Madagascar

Species delimitation within recently evolved groups can be challenging because species may be difficult to distinguish morphologically. Following the General Lineage Concept, we apply a multiple evidence approach to assess species limits within the carpet chameleon Furcifer lateralis, which is endemic to Madagascar and exported in large numbers for the pet trade. Cryptic speciation within F. lateralis was considered likely because this species (1) has a vast distribution, (2) occupies exceptionally diverse habitats and (3) exhibits subtle regional differences in morphology. Phylogenetic trees reconstructed using nuclear and mitochondrial genes recovered three well-supported clades corresponding with geography. Morphological results based on canonical variates analysis show that these clades exhibit subtle differences in head casque morphology. Ecological niche modelling results found that these phylogenetic groups also occupy unique environmental space and exhibit patterns of regional endemism typical of other endemic reptiles. Combined, our findings provide diverse yet consistent evidence for the existence of three species. Consequently, we elevate the subspecies F. lateralis major to species rank and name a new species distributed in northern and western Madagascar. Initial ecological divergence, associated with speciation of F. lateralis in humid eastern habitat, fits the Ecographic Constraint model for species diversification in Madagascar. By contrast, the second speciation event provides some support for the Riverine Barrier model, with the Mangoky River possibly causing initial isolation between species. These findings thus support two contrasting models of speciation within closely related species and demonstrate the utility of applying a combined-evidence approach for detecting cryptic speciation.     

Isolation and characterization of polymorphic microsatellite loci in the scarlet ibis (Eudocimus ruber— Threskiornithidae‐Aves)

Here we presented 10 polymorphic microsatellite loci obtained from scarlet ibis through an enriched genomic library. The analysis of 45 individuals from three Brazilian natural populations showed allelic diversity ranged from three to 17 alleles, observed heterozygosity ranged from 0.03 to 0.92, and expected heterozygosity ranged from 0.06 to 0.92. These highly variable microsatellite loci can provide means for assessing overall genetic variation in its remnant natural populations, which may help the development of effective conservation programs.     

Molecular genetic and quantitative trait divergence associated with recent homoploid hybrid speciation: a study of Senecio squalidus (Asteraceae)

Hybridization is increasingly seen as a trigger for rapid evolution and speciation. To quantify and qualify divergence associated with recent homoploid hybrid speciation, we compared quantitative trait (QT) and molecular genetic variation between the homoploid hybrid species Senecio squalidus and its parental species, S. aethnensis and S. chrysanthemifolius, and also their naturally occurring Sicilian hybrids. S. squalidus originated and became invasive in the United Kingdom following the introduction of hybrid plants from Mount Etna, Sicily, about 300 years ago. We recorded considerable molecular genetic differentiation between S. squalidus and its parents and their Sicilian hybrids in terms of both reduced genetic diversity and altered allele frequencies, potentially due to the genetic bottleneck associated with introduction to the United Kingdom. S. squalidus is also distinct from its parents and Sicilian hybrids for QTs, but less so than for molecular genetic markers. We suggest that this is due to resilience of polygenic QTs to changes in allele frequency or lack of selection for hybrid niche divergence in geographic isolation. While S. squalidus is intermediate or parental-like for most QTs, some trangressively distinct traits were observed, which might indicate emerging local adaptation in its invasive range. This study emphasizes the important contribution of founder events and geographic isolation to successful homoploid hybrid speciation.     

Evidence of recent and continuous speciation in a biodiversity hotspot: a population genetic approach in southern African gladioli (Gladiolus; Iridaceae)

There has been much debate over the origin of species diversity in biodiversity hotspots, particularly the rate of speciation over extinction and the geographic mode of speciation. Here, we looked at speciation with varying degrees of sympatry in a biodiversity hotspot, focusing on a distinct morphological clade in the Cape Floristic Region in southern Africa, the Gladiolus carinatus species complex (Iridaceae). We investigate the mechanisms involved in population and species differentiation through a combination of ecological and genomic approaches. We estimated spatial and phenological overlap, differences in floral morphology, genetic isolation and genomic selection. A genetic coalescent analysis estimated that the time of divergence between lineages followed the establishment of available habitat in the Cape littoral plain where these species currently overlap geographically. Marked shifts in flowering time and morphology, which act as barriers to gene flow, have developed to varying degrees over the last 0.3–1.4 million years. An amplified fragment length polymorphism genome scan revealed signatures of divergent and balancing selection, although half of the loci consistently behaved neutrally. Divergent species outliers (1%) and floral morph outliers (3%) represent a small proportion of the genome, but these loci produced clear genetic clusters of species and significant associations with floral traits. These results indicate that the G. carinatus complex represents a continuum of recent speciation. We provide further evidence for ecological adaptation in the face of gene flow.     

Host-associated speciation in the coral barnacle Wanella milleporae (Cirripedia: Pyrgomatidae) inhabiting the Millepora coral

Speciation by host shift is a common phenomenon observed in many symbiotic animals. The symbiont–host interaction is highly dynamic, but it is poorly documented in the marine realm. In the present study, we examined the genetic and morphological differentiation of the coral barnacle Wanella milleporae (obligate to fire corals) collected from four different Millepora host species in Taiwan to investigate the host specificity of this barnacle. Phylogenetic analysis of mitochondrial COI gene for 241 individuals of Wanella revealed five distinct clades, whose sequence divergences are comparable to values between other cogeneric barnacle species. The five clades also differ in shell and opercular plate morphology and colour. Genetic and morphological differentiations together strongly suggest the presence of cryptic species. Although the five clades do not display species-level host specificity, they showed a significant difference in preference on host growth form. Clades 1 and 2 were predominantly found on encrusting Millepora exaesa and Millepora platyphylla , while clades 3, 4 and 5 live exclusively on branching-form fire corals Millepora dichotoma and Millepora tenella . Phylogeny inferred from the combined mitochondrial COI, 16S and 12S (2182 bp) analysis suggests the division of the five clades into two major lineages congruent with the morphology of the host coral. Multiple independent invasions to the same form of host and subsequent speciation are evident in the Red Sea and Taiwan. Our results indicate that ecological/sympatric speciation could occur in marine symbiotic invertebrates through host shift and specialization. It appears that, as in their terrestrial counterparts, host–symbiont radiations in the marine realm are more prevalent than we expected and thus warrant further investigation.     

Quantitative trait and allozyme divergence in the Greenfinch (Carduelis chloris, Aves: Fringillidae)

Quantitative trait divergence and variability among 12 greenfinch populations across continental Europe was examined and compared to divergence in neutral genetic markers (allozymes). The added among locality variance component for 16 skeletal traits was large (mean 28%, range 4–48%) equalling a divergence of up to three SD units. The divergence in quantitative traits (Qst = 0.04-0.48) greatly exceeded that in alloymes (F_ST= 0.01-0.07), indicating the differentiation in quantitative traits to be larger than expected by mutation and drift alone. This conclusion was consistent also with results from the multivariate approach of Rogers & Harpending. However, genetic and morphometric distances between populations were positively correlated, even when controlling for the geographic distance separating pairs of populations. In concordance with Bergmann's rule, most traits were strongly and positively correlated with latitude, indicating latitudinally ordered genetic or/and environmental effects. However, the correlation between lower mandible width and latitude was strongly negative, demonstrating an inverse relationship between beak size and body size across the populations. These results are interpreted to reflect the re-colonization of history of northern Europe (genetic and geographic distances correlated) which has been paralleled by selection acting on quantitative traits (Q_ST>F_ST)- In particular, the counter-gradient variation in beak width, a functionally important trait, is suggestive of an adaptive basis for quantitative trait divergence.     

Genetic divergence and migration patterns in a North American passerine bird: implications for evolution and conservation

Like many other migratory birds, the black-throated blue warbler (Dendroica caerulescens) shows pronounced differences in migratory behaviour and other traits between populations: birds in the southern part of the breeding range have darker plumage and migrate to the eastern Caribbean during the winter, whereas those in the north have lighter plumage and migrate to the western Caribbean. We examined the phylogeography of this species, using samples collected from northern and southern populations, to determine whether differentiation between these populations dates to the Pleistocene or earlier, or whether differences in plumage and migratory behaviour have arisen more recently. We analysed variation at 369 bp of the mitochondrial control region domain I and also at seven nuclear microsatellites. Analyses revealed considerable genetic variation, but the vast majority of this variation was found within rather than between populations, and there was little differentiation between northern and southern populations. Phylogeographic analyses revealed a very shallow phylogenetic tree, a star-like haplotype network, and a unimodal mismatch distribution, all indicative of a recent range expansion from a single refugium. Coalescent modelling approaches also indicated a recent common ancestor for the entire group of birds analysed, no split between northern and southern populations, and high levels of gene flow. These results show that Pleistocene or earlier events have played little role in creating differences between northern and southern populations, suggesting that migratory and other differences between populations have arisen very recently. The implications of these results for the evolution of migration and defining taxonomic groups for conservation efforts are discussed.     

Isolation and characterization of microsatellite markers in the southern emu‐wren (Stipiturus malachurus: Aves)

We isolated and characterized eight novel microsatellite loci in the southern emu‐wren (Stipiturus malachurus). We used nonradioactive polymerase chain reaction (PCR)‐based techniques to screen an enriched genomic DNA library. Based on genotypes from a single population, six loci showed no evidence of null alleles and were polymorphic (allele range = 2–9, mean heterozygosity = 0.57), and one locus was sex‐linked (N_A = 4). These loci were variable and had different allele size ranges in three other populations of southern emu‐wrens, and are therefore useful for determining levels of genetic diversity within and between populations of the species.     

Patterns of genetic variation in the adaptive radiation of New World crossbills (Aves: Loxia)

Incipient species groups or young adaptive radiations such as crossbills (Aves: Loxia) present the opportunity to investigate directly the processes occurring during speciation. New World crossbills include white-winged crossbills (Loxia leucoptera), Hispaniolan crossbills (Loxia megaplaga), and red crossbills (Loxia curvirostra complex), the last of which is comprised of at least nine morphologically and vocally differentiated forms ('call types') where divergent natural selection for specialization on different conifer resources has been strongly implicated as driving diversification. Here we use amplified fragment length polymorphism (AFLP) markers to investigate patterns of genetic variation across populations, call types, and species of New World crossbills. Tree-based analyses using 440 AFLP loci reveal strongly supported clustering of the formally recognized species, but did not separate individuals from the eight call types in the red crossbill complex, consistent with recent divergence and ongoing gene flow. Analyses of genetic differentiation based on inferred allele frequency variation however, reveal subtle but significant levels of genetic differentiation among the different call types of the complex and indicate that between call-type differentiation is greater than that found among different geographic locations within call types. Interpreted in light of evidence of divergent natural selection and strong premating reproductive isolation, the observed genetic differentiation suggests restricted gene flow among sympatric call types consistent with the early stages of ecological speciation.     

Reconciling molecular signatures across markers: mitochondrial DNA confirms founder effect in invasive North American house finches (<Emphasis Type="Italic">Carpodacus mexicanus</Emphasis>)

Well-characterized species introductions provide opportunities to compare the genetic signatures of known founder effects across classes of molecular markers. The release of small numbers of house finches (Carpodacus mexicanus) into the eastern United States in the 1940s led to substantial interest in the effects of this introduction on genetic diversity in this now abundant species, an issue that has been highlighted by a recent Mycoplasma disease epidemic that most intensively affects the introduced and potentially genetically depauperate house finch populations. Previous studies comparing genetic diversity levels in native and introduced house finch populations produced seemingly disparate results: comparisons based on amplified fragment length polymorphism, RFLP mtDNA, and allozyme markers found essentially equivalent levels of diversity in eastern and western populations, whereas microsatellite markers showed clear reductions in diversity in the introduced populations. Here we employ sequence variation at the ND2 mtDNA locus to further compare levels of diversity between the four native and five introduced house finch populations that were previously examined in the microsatellite study. We found substantially lower ND2 haplotype richness and diversity across all introduced populations of house finches. The majority of sequence variation (78%) was detected within subpopulations, with the remainder (22%) explained by the historical status of each population (native or introduced). Our results are consistent with previous microsatellite evidence for a founder effect during the introduction of eastern house finches, and suggest that the mtDNA founder effect was particularly severe, likely owing to a male-biased sex ratio at the time of introduction coupled with the lower effective population size of clonally inherited markers. We discuss how the inconsistencies between past studies of house finch diversity can inform the usefulness of distinct marker sets for detecting molecular signatures of founder events.     

Phylogeography and colonization history of Lesser Black‐backed Gulls (Larus fuscus) as revealed by mtDNA sequences

Abstract Because of the differential amplitude of climatic oscillations, species living at northern latitudes are subject to more frequent and more severe range oscillations than species at southern latitudes. As a consequence, northern populations should, on average, be phylogenetically younger and possess less phylogeographical structure than closely related taxa further south. To test these predictions, we studied the mitochondrial‐genetic population structure of NW Palearctic Lesser Black‐backed Gulls (Larus fuscus group [=LBBG], five taxa) breeding at temperate to boreal latitudes from Iceland to the Taimyr Peninsula. Results were compared with those previously obtained (Liebers et al. 2001. Mol. Ecol. 10 : 2447) for more southerly breeding Yellow‐legged Gulls (Larus cachinnans group, six taxa from the Atlantic Islands to Mongolia). Sequences of the hypervariable region I (HVR‐I) of the mitochondrial control region revealed low within‐ and between‐taxon sequence divergence, little genetic variation, a shallow haplotype phylogeny and poor phylogeographical structure in LBBGs compared with Yellow‐legged Gulls. Haplotype frequencies among the five northern taxa formed a stepped cline with significant gene flow restriction between the forms heuglini and fuscus, probably indicating a secondary contact with (partial?) reproductive isolation. Western forms of LBBG, among which mitochondrial gene flow appears unrestricted, show genetic signs of postglacial range expansion and population growth. The Larus fuscus group is derived from a cachinnans‐like ancestral population, probably in the Aralo‐Caspian basin, and spread from east (NW Siberia) to west within the Palearctic.     

Re-examination of a proposed case of stasipatric speciation: phylogeography of the Australian morabine grasshoppers (Vandiemenella viatica species group)

Karyotypic differences have been used for delimiting populations or species, although whether these mutations provide strong barriers to gene flow between populations and promote speciation remains contentious. In this study, we assessed whether 11 chromosomal races of Australian morabine grasshoppers ( Vandiemenella viatica species group) represent genetically distinct populations by analyses of cytological and allozyme (35 loci) data and DNA sequences of the elongation factor-1 alpha ( EF -1α), anonymous Mvia11 , and mitochondrial cytochrome c oxidase subunit I ( COI ) loci. While the Vandiemenella chromosomal taxa generally represent genetically distinct units, a substantial portion of the total genetic variation in our samples was not explained by the chromosomal variation. Mantel tests indicated that Vandiemenella populations were spatially structured and have maintained gene flow at a local scale within each of the taxa. The group was subdivided into 13 genetic clusters; four chromosomal taxa comprised single exclusive clusters, while others comprised more than one cluster or clusters shared with other taxa. Boundaries of these cryptic population subdivisions correspond with several biogeographical barriers, such as straits, gulfs, the Murray River, and an ancient mega-lake, Lake Bungunnia. The viatica species group was previously proposed to have diversified without major geographical separation based on the stasipatric speciation model; however, the present study suggests the involvement of allopatric fragmentation. Given extensive nonmonophyly of chromosomal taxa and incomplete barriers to gene flow among taxa, all Vandiemenella chromosomal taxa and genetically distinct populations within chromosomal taxa, except Vandiemenella pichirichi , should be regarded as populations of one species: Vandiemenella viatica .