Gut microbial diversity across a contact zone for California voles: Implications for lineage divergence of hosts and mitonuclear mismatch in the assembly of the mammalian gut microbiome

Gut microbial diversity is thought to reflect the co‐evolution of microbes and their hosts as well as current host‐specific attributes such as genetic background and environmental setting. To explore interactions among these parameters, we characterized variation in gut microbiome composition of California voles (Microtus californicus) across a contact zone between two recently diverged lineages of this species. Because this contact zone contains individuals with mismatched mitochondrial‐nuclear genomes (cybrids), it provides an important opportunity to explore how different components of the genotype contribute to gut microbial diversity. Analyses of bacterial 16S rRNA sequences and joint species distribution modelling revealed that host genotypes and genetic differentiation among host populations together explained more than 50% of microbial community variation across our sampling transect. The ranked importance (most to least) of factors contributing to gut microbial diversity in our study populations were: genome‐wide population differentiation, local environmental conditions, and host genotypes. However, differences in microbial communities among vole populations (β‐diversity) did not follow patterns of lineage divergence (i.e., phylosymbiosis). Instead, among‐population variation was best explained by the spatial distribution of hosts, as expected if the environment is a primary source of gut microbial diversity (i.e., dispersal limitation hypothesis). Across the contact zone, several bacterial taxa differed in relative abundance between the two parental lineages as well as among individuals with mismatched mitochondrial and nuclear genomes. Thus, genetic divergence among host lineages and mitonuclear genomic mismatches may also contribute to microbial diversity by altering interactions between host genomes and gut microbiota (i.e., hologenome speciation hypothesis).     

Divergence and gene flow in the globally distributed blue‐winged ducks

The ability to disperse over long distances can result in a high propensity for colonizing new geographic regions, including uninhabited continents, and lead to lineage diversification via allopatric speciation. However, high vagility can also result in gene flow between otherwise allopatric populations, and in some cases, parapatric or divergence‐with‐gene‐flow models might be more applicable to widely distributed lineages. Here, we use five nuclear introns and the mitochondrial control region along with Bayesian models of isolation with migration to examine divergence, gene flow, and phylogenetic relationships within a cosmopolitan lineage comprising six species, the blue‐winged ducks (genus Anas), which inhabit all continents except Antarctica. We found two primary sub‐lineages, the globally‐distributed shoveler group and the New World blue‐winged/cinnamon teal group. The blue‐winged/cinnamon sub‐lineage is composed of sister taxa from North America and South America, and taxa with parapatric distributions are characterized by low to moderate levels of gene flow. In contrast, our data support strict allopatry for most comparisons within the shovelers. However, we found evidence of gene flow from the migratory, Holarctic northern shoveler A. clypeata and the more sedentary, African Cape shoveler A. smithii into the Australasian shoveler A. rhynchotis, although we could not reject strict allopatry. Given the diverse mechanisms of speciation within this complex, the shovelers and blue‐winged/cinnamon teals can serve as an effective model system for examining how the genome diverges under different evolutionary processes and how genetic variation is partitioned among highly dispersive taxa.     

Cyto‐nuclear discordance suggests complex evolutionary history in the cave‐dwelling salamander,Eurycea lucifuga

Our understanding of the evolutionary history and ecology of cave‐associated species has been driven historically by studies of morphologically adapted cave‐restricted species. Our understanding of the evolutionary history and ecology of nonrestricted cave species, troglophiles, is limited to a few studies, which present differing accounts of troglophiles’ relationship with the cave habitat, and its impact on population dynamics. Here, we used phylogenetics, demographic statistics, and population genetic methods to study lineage divergence, dates of divergence, and population structure in the Cave Salamander, Eurycea lucifuga, across its range. In order to perform these analyses, we sampled 233 individuals from 49 populations, using sequence data from three gene loci as well as genotyping data from 19 newly designed microsatellite markers. We find, as in many other species studied in a phylogeographic context, discordance between patterns inferred from mitochondrial relationships and those inferred by nuclear markers indicating a complicated evolutionary history in this species. Our results suggest Pleistocene‐based divergence among three main lineages within E. lucifuga corresponding to the western, central, and eastern regions of the range, similar to patterns seen in species separated in multiple refugia during climatic shifts. The conflict between mitochondrial and nuclear patterns is consistent with what we would expect from secondary contact between regional populations following expansion from multiple refugia.     

Genetic surfing,not allopatric divergence,explains spatial sorting of mitochondrial haplotypes in venomous coralsnakes

Strong spatial sorting of genetic variation in contiguous populations is often explained by local adaptation or secondary contact following allopatric divergence. A third explanation, spatial sorting by stochastic effects of range expansion, has been considered less often though theoretical models suggest it should be widespread, if ephemeral. In a study designed to delimit species within a clade of venomous coralsnakes, we identified an unusual pattern within the Texas coral snake (Micrurus tener): strong spatial sorting of divergent mitochondrial (mtDNA) lineages over a portion of its range, but weak sorting of these lineages elsewhere. We tested three alternative hypotheses to explain this pattern—local adaptation, secondary contact following allopatric divergence, and range expansion. Collectively, near panmixia of nuclear DNA, the signal of range expansion associated sampling drift, expansion origins in the Gulf Coast of Mexico, and species distribution modeling suggest that the spatial sorting of divergent mtDNA lineages within M. tener has resulted from genetic surfing of standing mtDNA variation—not local adaptation or allopatric divergence. Our findings highlight the potential for the stochastic effects of recent range expansion to mislead estimations of population divergence made from mtDNA, which may be exacerbated in systems with low vagility, ancestral mtDNA polymorphism, and male‐biased dispersal.     

Host specificity in spinturnicid mites: do parasites share a long evolutionary history with their host?

Host specificity in parasites can be explained by spatial isolation from other potential hosts or by specialization and speciation of specific parasite species. The first assertion is based on allopatric speciation, the latter on differential lifetime reproductive success on different available hosts. We investigated the host specificity and cophylogenetic histories of four sympatric European bat species of the genus Myotis and their ectoparasitic wing mites of the genus Spinturnix. We sampled >40 parasite specimens from each bat species and reconstructed their phylogenetic COI trees to assess host specificity. To test for cospeciation, we compared host and parasite trees for congruencies in tree topologies. Corresponding divergence events in host and parasite trees were dated using the molecular clock approach. We found two species of wing mites to be host specific and one species to occur on two unrelated hosts. Host specificity cannot be explained by isolation of host species, because we found individual parasites on other species than their native hosts. Furthermore, we found no evidence for cospeciation, but for one host switch and one sorting event. Host‐specific wing mites were several million years younger than their hosts. Speciation of hosts did not cause speciation in their respective parasites, but we found that diversification of recent host lineages coincided with a lineage split in some parasites.     

DOES NICHE DIVERGENCE ACCOMPANY ALLOPATRIC DIVERGENCE IN APHELOCOMA JAYS AS PREDICTED UNDER ECOLOGICAL SPECIATION?: INSIGHTS FROM TESTS WITH NICHE MODELS

The role of ecology in the origin of species has been the subject of long‐standing interest to evolutionary biologists. New sources of spatially explicit ecological data allow for large‐scale tests of whether speciation is associated with niche divergence or whether closely related species tend to be similar ecologically (niche conservatism). Because of the confounding effects of spatial autocorrelation of environmental variables, we generate null expectations for niche divergence for both an ecological‐niche modeling and a multivariate approach to address the question: do allopatrically distributed taxa occupy similar niches? In a classic system for the study of niche evolution—the Aphelocoma jays—we show that there is little evidence for niche divergence among Mexican Jay (A. ultramarina) lineages in the process of speciation, contrary to previous results. In contrast, Aphelocoma species that exist in partial sympatry in some regions show evidence for niche divergence. Our approach is widely applicable to the many cases of allopatric lineages in the beginning stages of speciation. These results do not support an ecological speciation model for Mexican Jay lineages because, in most cases, the allopatric environments they occupy are not significantly more divergent than expected under a null model.     

Cryptic speciation in the field vole: a multilocus approach confirms three highly divergent lineages in Eurasia

Species are generally described from morphological features, but there is growing recognition of sister forms that show substantial genetic differentiation without obvious morphological variation and may therefore be considered ‘cryptic species’. Here, we investigate the field vole (Microtus agrestis), a Eurasian mammal with little apparent morphological differentiation but which, on the basis of previous sex‐linked nuclear and mitochondrial DNA (mtDNA) analyses, is subdivided into a Northern and a Southern lineage, sufficiently divergent that they may represent two cryptic species. These earlier studies also provided limited evidence for two major mtDNA lineages within Iberia. In our present study, we extend these findings through a multilocus approach. We sampled 163 individuals from 46 localities, mainly in Iberia, and sequenced seven loci, maternally, paternally and biparentally inherited. Our results show that the mtDNA lineage identified in Portugal is indeed a distinct third lineage on the basis of other markers as well. In fact, multilocus coalescent‐based methods clearly support three separate evolutionary units that may represent cryptic species: Northern, Southern and Portuguese. Divergence among these units was inferred to have occurred during the last glacial period; the Portuguese lineage split occurred first (estimated at c. 70 000 bp ), and the Northern and Southern lineages separated at around the last glacial maximum (estimated at c. 18 500 bp ). Such recent formation of evolutionary units that might be considered species has repercussions in terms of understanding evolutionary processes and the diversity of small mammals in a European context.     

Multilocus phylogeography reveals nested endemism in a gecko across the monsoonal tropics of Australia

Multilocus phylogeography can uncover taxonomically unrecognized lineage diversity across complex biomes. The Australian monsoonal tropics include vast, ecologically intact savanna‐woodland plains interspersed with ancient sandstone uplands. Although recognized in general for its high species richness and endemism, the biodiversity of the region remains underexplored due to its remoteness. This is despite a high rate of ongoing species discovery, especially in wetter regions and for rock‐restricted taxa. To provide a baseline for ongoing comparative analyses, we tested for phylogeographic structure in an ecologically generalized and widespread taxon, the gecko Heteronotia binoei. We apply coalescent analyses to multilocus sequence data (mitochondrial DNA and eight nuclear DNA introns) from individuals sampled extensively and at fine scale across the region. The results demonstrate surprisingly deep and geographically nested lineage diversity. Several intra‐specific clades previously shown to be endemic to the region were themselves found to contain multiple, short‐range lineages. To infer landscapes with concentrations of unique phylogeographic diversity, we probabilistically estimate the ranges of lineages from point data and then, combining these estimates with the nDNA species tree, estimate phyloendemism across the region. Highest levels of phyloendemism occur in northern Top End, especially on islands, across the topographically complex Arnhem escarpment, and across the sandstone ranges of the western Gulf region. These results drive home that deep phylogeographic structure is prevalent in tropical low‐dispersal taxa, even ones that are ubiquitous across geography and habitats.     

A genomic assessment of population structure and gene flow in an aquatic salamander identifies the roles of spatial scale,barriers, and river architecture

相似文献   

Cryptic lineages of a common alpine mayfly show strong life‐history divergence

Understanding ecological divergence of morphologically similar but genetically distinct species – previously considered as a single morphospecies – is of key importance in evolutionary ecology and conservation biology. Despite their morphological similarity, cryptic species may have evolved distinct adaptations. If such ecological divergence is unaccounted for, any predictions about their responses to environmental change and biodiversity loss may be biased. We used spatio‐temporally replicated field surveys of larval cohort structure and population genetic analyses (using nuclear microsatellite markers) to test for life‐history divergence between two cryptic lineages of the alpine mayfly Baetis alpinus in the Swiss Alps. We found that the more widespread and abundant cryptic lineage represents a ‘generalist’ with at least two cohorts per year, whereas the less abundant lineage is restricted to higher elevations and represents a ‘specialist’ with a single cohort per year. Importantly, our results indicate partial temporal segregation in reproductive periods between these lineages, potentially facilitating local coexistence and reproductive isolation. Taken together, our findings emphasize the need for a taxonomic revision: widespread and apparently generalist morphospecies can hide cryptic lineages with much narrower ecological niches and distribution ranges.     

New Species of Closely Related Endosymbiotic Dinoflagellates in the Greater Caribbean have Niches Corresponding to Host Coral Phylogeny

Symbiotic dinoflagellates in the genus Breviolum (formerly Symbiodinium Clade B) dominate coral communities in shallow waters across the Greater Caribbean. While some formally described species exist, mounting genetic, and ecological evidence indicate that numerous more comprise this genus, many of which are closely related. To test this, colonies of common reef‐building corals were sampled across a large geographical range. Phylogenetic and population genetic markers then used to examine evolutionary divergence and delineate boundaries of genetic recombination. Three new candidate species were distinguished by fixed differences in nucleotide sequences from nuclear and chloroplast DNA. Population connectivity was evident within each lineage over thousands of kilometers, however, substantial genetic structure persisted between lineages co‐occurring within sampling locations, signifying reproductive isolation. While geographically widespread with overlapping distributions, each species is ecologically distinct, exhibiting specific mutualisms with phylogenetically distinct coral hosts. Moreover, significant differences in mean cell sizes provide some morphological evidence substantiating formal species distinctions. In providing evidence that satisfies the biological, phylogenetic, ecological, and morphological species concepts, we classify and formally name Breviolum faviinorum n. sp., primarily associated with Caribbean corals belonging to the Caribbean subfamily Faviinae; B. meandrinium n. sp., associated with corals belonging to the family Meandrinidae; and B. dendrogyrum n. sp., a symbiont harbored exclusively by the threatened coral Dendrogyra cylindrus. These findings support the primary importance of niche diversification (i.e. host habitat) in the speciation of symbiotic dinoflagellates.     

Range instability leads to cytonuclear discordance in a morphologically cryptic ground squirrel species complex

The processes responsible for cytonuclear discordance frequently remain unclear. Here, we employed an exon capture data set and demographic methods to test hypotheses generated by species distribution models to examine how contrasting histories of range stability vs. fluctuation have caused cytonuclear concordance and discordance in ground squirrel lineages from the Otospermophilus beecheyi species complex. Previous studies in O. beecheyi revealed three morphologically cryptic and highly divergent mitochondrial DNA lineages (named the Northern, Central and Southern lineages based on geography) with only the Northern lineage exhibiting concordant divergence for nuclear genes. Here, we showed that these mtDNA lineages likely formed in allopatry during the Pleistocene, but responded differentially to climatic changes that occurred since the last interglacial (~120,000 years ago). We find that the Northern lineage maintained a stable range throughout this period, correlating with genetic distinctiveness among all genetic markers and low migration rates with the other lineages. In contrast, our results suggested that the Southern lineage expanded from Baja California Sur during the Late Pleistocene to overlap and potentially swamp a contracting Central lineage. High rates of intraspecific gene flow between Southern lineage individuals among expansion origin and expansion edge populations largely eroded Central ancestry from autosomal markers. However, male‐biased dispersal in this system preserved signals of this past hybridization and introgression event in matrilineal‐biased X‐chromosome and mtDNA markers. Our results highlight the importance of range stability in maintaining the persistence of phylogeographic lineages, whereas unstable range dynamics can increase the tendency for lineages to merge upon secondary contact.     

Population and phylogenomic decomposition via genotyping‐by‐sequencing in Australian Pelargonium

Species delimitation has seen a paradigm shift as increasing accessibility of genomic‐scale data enables separation of lineages with convergent morphological traits and the merging of recently diverged ecotypes that have distinguishing characteristics. We inferred the process of lineage formation among Australian species in the widespread and highly variable genus Pelargonium by combining phylogenomic and population genomic analyses along with breeding system studies and character analysis. Phylogenomic analysis and population genetic clustering supported seven of the eight currently described species but provided little evidence for differences in genetic structure within the most widely distributed group that containing P. australe. In contrast, morphometric analysis detected three deep lineages within Australian Pelargonium; with P. australe consisting of five previously unrecognized entities occupying separate geographic ranges. The genomic approach enabled elucidation of parallel evolution in some traits formerly used to delineate species, as well as identification of ecotypic morphological differentiation within recognized species. Highly variable morphology and trait convergence each contribute to the discordance between phylogenomic relationships and morphological taxonomy. Data suggest that genetic divergence among species within the Australian Pelargonium may result from allopatric speciation while morphological differentiation within and among species may be more strongly driven by environmental differences.     

Divergence and biogeography of the recently evolved Macaronesian red Festuca (Gramineae) species inferred from coalescence-based analyses

Studying the biogeography and the phylogeography of the endemic Macaronesian red Festuca species (Loliinae, Poaceae) is of prime interest in understanding the speciation and colonization patterns of recently evolved groups in oceanic archipelagos. Coalescence‐based analyses of plastid trnLF sequences were employed to estimate evolutionary parameters and to test different species‐history scenarios that model the pattern of species divergence. Bayesian IM estimates of species divergence times suggested that ancestral lineages of diploid Macaronesian and Iberian red fescues could have diverged between 1.2 and 1.57 Ma. When empirical data were compared to coalescence‐based simulated distributions of discordance and p‐distance statistics, two species‐history models were chosen in which the first branching lineage derived in Canarian Festuca agustinii. Its sister lineage could have involved a recent polytomy leading to the Madeiran Festuca jubata, the Azorean Festuca francoi + Festuca petraea and the continental Festuca rivularis lineages (Canarian model) or the sequential branching of lineages leading to F. jubata and finally to the sister clades of F. rivularis and F. francoi + F. petraea (Sequential model). Nested clade phylogeographic analysis (NCPA) and a first adapted host–parasite co‐evolutionary ParaFit method were used to detect the phylogeographic signal. NCPA inferred long‐distance colonizations for the entire diploid red Festuca complex, but allopatric‐fragmentation and isolation‐by‐distance (IBD) patterns were inferred within archipelagos. In addition, the ParaFit method suggested a generalized pattern of a stepping‐stone model at all hierarchical levels. Maximum‐likelihood‐based dispersal‐extinction‐cladogenesis (DEC) models were superimposed on the Sequential model species tree. The three‐independent‐colonization (3IC) model was the best supported biogeographic scenario, concurring with previous analysis based on multilocus AFLP data.     

A genomic assessment of species boundaries and hybridization in a group of highly polymorphic anoles (distichus species complex)

Delimiting young species is one of the great challenges of systematic biology, particularly when the species in question exhibit little morphological divergence. Anolis distichus, a trunk anole with more than a dozen subspecies that are defined primarily by dewlap color, may actually represent several independent evolutionary lineages. To test this, we utilized amplified fragment length polymorphisms (AFLP) genome scans and genetic clustering analyses in conjunction with a coalescent‐based species delimitation method. We examined a geographically widespread set of samples and two heavily sampled hybrid zones. We find that genetic divergence is associated with a major biogeographic barrier, the Hispaniolan paleo‐island boundary, but not with dewlap color. Additionally, we find support for hypotheses regarding colonization of two Hispaniolan satellite islands and the Bahamas from mainland Hispaniola. Our results show that A. distichus is composed of seven distinct evolutionary lineages still experiencing a limited degree of gene flow. We suggest that A. distichus merits taxonomic revision, but that dewlap color cannot be relied upon as the primary diagnostic character.     

Range‐wide analysis of genetic structure in a widespread,highly mobile species (Odocoileus hemionus) reveals the importance of historical biogeography

Highly mobile species that thrive in a wide range of habitats are expected to show little genetic differentiation across their range. A limited but growing number of studies have revealed that patterns of broad‐scale genetic differentiation can and do emerge in vagile, continuously distributed species. However, these patterns are complex and often shaped by both historical and ecological factors. Comprehensive surveys of genetic variation at a broad scale and at high resolution are useful for detecting cryptic spatial genetic structure and for investigating the relative roles of historical and ecological processes in structuring widespread, highly mobile species. In this study, we analysed 10 microsatellite loci from over 1900 samples collected across the full range of mule deer (Odocoileus hemionus), one of the most widely distributed and abundant of all large mammal species in North America. Through both individual‐ and population‐based analyses, we found evidence for three main genetic lineages, one corresponding to the ‘mule deer’ morphological type and two to the ‘black‐tailed deer’ type. Historical biogeographic events likely are the primary drivers of genetic divergence in this species; boundaries of the three lineages correspond well with predictions based on Pleistocene glacial cycles, and substructure within each lineage demonstrates island vicariance. However, across large geographic areas, including the entire mule deer lineage, we found that genetic variation fit an isolation‐by‐distance pattern rather than discrete clusters. A lack of genetic structure across wide geographic areas of the continental west indicates that ecological processes have not resulted in restrictions to gene flow sufficient for spatial genetic structure to emerge. Our results have important implications for our understanding of evolutionary mechanisms of divergence, as well as for taxonomy, conservation and management.     

Maintaining their genetic distance: Little evidence for introgression between widely hybridizing species of Geum with contrasting mating systems

Within the plant kingdom, many genera contain sister lineages with contrasting outcrossing and inbreeding mating systems that are known to hybridize. The evolutionary fate of these sister lineages is likely to be influenced by the extent to which they exchange genes. We measured gene flow between outcrossing Geum rivale and selfing Geum urbanum, sister species that hybridize in contemporary populations. We generated and used a draft genome of G. urbanum to develop dd‐RAD data scorable in both species. Coalescent analysis of RAD data from allopatric populations indicated that the species diverged 2–3 Mya, and that historical gene flow between them was extremely low (1 migrant every 25 generations). Comparison of genetic divergence between species in sympatry and allopatry, together with an analysis of allele frequencies in potential parental and hybrid populations, provided no evidence of contemporary introgression in sympatric populations. Cluster‐ and species‐specific marker analyses revealed that, apart from four early‐generation hybrids, individuals in sympatric populations fell into two genetically distinct groups that corresponded exactly to their morphological species classification with maximum individual admixture estimates of only 1–3%. However, we did observe joint segregation of four putatively introgressed SNPs across two scaffolds in the G. urbanum population that was associated with significant morphological variation, interpreted as tentative evidence for rare, recent interspecific gene flow. Overall, our results indicate that despite the presence of hybrids in contemporary populations, genetic exchange between G. rivale and G. urbanum has been extremely limited throughout their evolutionary history.     

Cryptic diversity hides host and habitat specialization in a gorgonian‐algal symbiosis

Shallow water anthozoans, the major builders of modern coral reefs, enhance their metabolic and calcification rates with algal symbionts. Controversy exists over whether these anthozoan–algae associations are flexible over the lifetimes of individual hosts, promoting acclimative plasticity, or are closely linked, such that hosts and symbionts co‐evolve across generations. Given the diversity of algal symbionts and the morphological plasticity of many host species, cryptic variation within either partner could potentially confound studies of anthozoan‐algal associations. Here, we used ribosomal, organelle and nuclear sequences, along with microsatellite variation, to study the relationship between lineages of a common Caribbean gorgonian and its algal symbionts. The gorgonian Eunicea flexuosa is a broadcast spawner, composed of two recently diverged, genetically distinct lineages largely segregated by depth. We sampled colonies of the two lineages across depth gradients at three Caribbean locations. We find that each host lineage is associated with a unique Symbiodinium B1/184 phylotype. This relationship between host and symbiont is maintained when host colonies are reciprocally transplanted, although cases of within phylotype switching were also observed. Even when the phylotypes of both partners are present at intermediate depths, the specificity between host and symbiont lineages remained absolute. Unrecognized cryptic diversity may mask host‐symbiont specificity and change the inference of evolutionary processes in mutualistic associations. Symbiotic specificity thus likely contributes to the ecological divergence of the two partners, generating species diversity within coral reefs.     

Understanding the mechanisms of antitropical divergence in the seabird White‐faced Storm‐petrel (Procellariiformes: Pelagodroma marina) using a multilocus approach

Analytical methods that apply coalescent theory to multilocus data have improved inferences of demographic parameters that are critical to understanding population divergence and speciation. In particular, at the early stages of speciation, it is important to implement models that accommodate conflicting gene trees, and benefit from the presence of shared polymorphisms. Here, we employ eleven nuclear loci and the mitochondrial control region to investigate the phylogeography and historical demography of the pelagic seabird White‐faced Storm‐petrel (Pelagodroma marina) by sampling subspecies across its antitropical distribution. Groups are all highly differentiated: global mitochondrial Φ_ST = 0.89 (P < 0.01) and global nuclear Φ_ST varies between 0.22 and 0.83 (all P < 0.01). The complete lineage sorting of the mitochondrial locus between hemispheres is corroborated by approximately half of the nuclear genealogies, suggesting a long‐term antitropical divergence in isolation. Coalescent‐based estimates of demographic parameters suggest that hemispheric divergence of P. marina occurred approximately 840 000 ya (95% HPD 582 000–1 170 000), in the absence of gene flow, and divergence within the Southern Hemisphere occurred 190 000 ya (95% HPD 96 000–600 000), both probably associated with the profound palaeo‐oceanographic changes of the Pleistocene. A fledgling sampled in St Helena (tropical South Atlantic) suggests recent colonization from the Northern Hemisphere. Despite the great potential for long‐distance dispersal, P. marina antitropical groups have been evolving as independent, allopatric lineages, and divergence is probably maintained by philopatry coupled with asynchronous reproductive phenology and local adaptation.