Speciation in Passerina buntings: introgression patterns of sex-linked loci identify a candidate gene region for reproductive isolation

Sex-chromosomes are thought to play an important role in speciation, but few studies of non-model organisms have investigated the relative influence of multiple sex-linked markers on reproductive isolation. We collected 222 individuals along a geographical transect spanning the hybrid zone between Passerina amoena and P. cyanea (Aves: Cardinalidae). Using maximum-likelihood cline fitting methods, we estimated locus-specific introgression rates for 10 z-linked markers. Although the cline width estimates ranged from 2.8 to 584 km, eight of 10 loci had cline widths between 224 and 271 km. We also used coalescent-based estimates of locus-specific divergence times between P. amoena and P. cyanea to test a recently proposed hypothesis of an inverse relationship between divergence time and cline width but did not find a significant association. The narrow width (2.8 km) of the cline estimated from the VLDLR9 locus indicates strong selection retarding introgression of alleles at this locus across the hybrid zone. Interestingly, a mutation in the very low density lipoprotein receptor ( VLDLR ) gene, in which VLDLR9 is an intron, is known to reduce the egg-laying ability of some chickens, suggesting a possible link between this gene region and reproductive isolation between P. amoena and P. cyanea . These results underscore the importance of sampling multiple loci to investigate introgression patterns across a chromosome or genome and support previous findings of the importance of sex-linked genes in speciation.     

Spatio-temporal changes in the genetic structure of the Passerina bunting hybrid zone

Although theoretical models predict that the structure of a hybrid zone can change under a variety of scenarios, only a few empirical studies of hybrid zones have unequivocally demonstrated zone movement. These studies are rare because few data sets exist that include repeated, temporally spaced, samples of the same hybrid zone. We analysed mitochondrial DNA haplotype data from samples separated by 40-45 years from across the Passerina amoena (Lazuli Bunting) and Passerina cyanea (Indigo Bunting) hybrid zone to investigate whether the genetic structure of this zone has changed during that interval. Both cline and generalized linear mixed modelling analyses uncovered a significant narrowing and a substantial westward shift of the Passerina bunting hybrid zone, clearly illustrating hybrid zone movement. The cause of the change may be due to a combination of ecological, demographic and behavioural factors. Our results predict that the width of the hybrid zone will continue to narrow over time, a finding consistent with reinforcement theory.     

Haldane's rule in an avian system: using cline theory and divergence population genetics to test for differential introgression of mitochondrial, autosomal, and sex-linked loci across the Passerina bunting hybrid zone

Using cline fitting and divergence population genetics, we tested a prediction of Haldane's rule: autosomal alleles should introgress more than z-linked alleles or mitochondrial haplotypes across the Passerina amoena/Passerina cyanea (Aves: Cardinalidae) hybrid zone. We screened 222 individuals collected along a transect in the Great Plains of North America that spans the contact zone for mitochondrial (two genes), autosomal (four loci) and z-linked (two loci) markers. Maximum-likelihood cline widths estimated from the mitochondrial (223 km) and z-linked (309 km) datasets were significantly narrower on average than the autosomal cline widths (466 km). We also found that mean coalescent-based estimates of introgression were larger for the autosomal loci (0.63 genes/generation, scaled to the mutation rate mu) than for both the mitochondrial (0.27) and z-linked loci (0.59). These patterns are consistent with Haldane's rule, but the among-locus variation also suggests many independently segregating loci are required to investigate introgression patterns across the genome. These results provide the first comprehensive comparison of mitochondrial, sex-linked, and autosomal loci across an avian hybrid zone and add to the body of evidence suggesting that sex chromosomes play an important role in the formation and maintenance of reproductive isolation between closely related species.     

Phylogeography and systematics of the Peromyscus eremicus species group and the historical biogeography of North American warm regional deserts

Phylogeographic relationships among 26 populations from throughout the geographic range of the Peromyscus eremicus species group are described based on sequence data for a 699-bp fragment of the mitochondrial DNA COIII gene. Distance, maximum-likelihood, and maximum-parsimony analyses of phylogenetic trees generated under four separate character-weighting strategies and representing five alternative biogeographic hypotheses revealed the existence of a cryptic species (Peromyscus fraterculus, previously included under P. eremicus) on the Baja California Peninsula and adjacent southwestern California and two distinct forms of P. eremicus, one from the Mojave, Sonoran, and northwestern Chihuahuan regional deserts (West) and one from the remainder of the Chihuahuan Desert (East). Distinctiveness of P. fraterculus is supported by previous morphometric and allozyme analyses, including comparisons with neighboring P. eremicus and parapatric P. eva, with which P. fraterculus shares a sister taxon relationship. Divergence of the eva + fraterculus, West + East eremicus, and P. merriami haplotype lineages likely occurred in the late Neogene (3 Ma), in response to northern extension of the Sea of Cortéz and elevation of the Sierra Madre Occidental; divergence of eva from fraterculus is concordant with the existence of a trans-Peninsular seaway during the Pleistocene (1 Ma); and divergence of West from East eremicus occurred during the Pleistocene pluvial-interpluvial cycles, but well before the Wisconsinan glacial interval. The sequence of divergence within the eremicus species group and causal association of geological events of the Neogene and Holocene provide a working hypothesis against which phylogeographic patterns among other arid-adapted species of the warm regional deserts of North America may be compared.     

Molecular phylogenetics at the population/species interface in cave spiders of the southern Appalachians (Araneae:Nesticidae:Nesticus)

This paper focuses on the relationship between population genetic structure and speciation mechanisms in a monophyletic species group of Appalachian cave spiders (Nesticus). Using mtDNA sequence data gathered from 256 individuals, I analyzed patterns of genetic variation within and between populations for three pairs of closely related sister species. Each sister-pair comparison involves taxa with differing distributional and ecological attributes; if these ecological attributes are reflected in basic demographic differences, then speciation might proceed differently across these sister taxa comparisons. Both frequency-based and gene tree analyses reveal that the genetic structure of the Nesticus species studied is characterized by similar and essentially complete population subdivision, regardless of differences in general ecology. These findings contrast with results of prior genetic studies of cave-dwelling arthropods that have typically revealed variation in population structure corresponding to differences in general ecology. Species fragmentation through both extrinsic and intrinsic evolutionary forces has resulted in discrete, perhaps independent, populations within morphologically defined species. Large sequence divergence values observed between populations suggest that this independence may extend well into the past. These patterns of mtDNA genealogical structure and divergence imply that species as morphological lineages are currently more inclusive than basal evolutionary or phylogenetic units, a suggestion that has important implications for the study of speciation mechanisms.      

HISTORICAL DIVERGENCE AND GENE FLOW: COALESCENT ANALYSES OF MITOCHONDRIAL,AUTOSOMAL AND SEX‐LINKED LOCI IN PASSERINA BUNTINGS

Quantifying the role of gene flow during the divergence of closely related species is crucial to understanding the process of speciation. We collected DNA sequence data from 20 loci (one mitochondrial, 13 autosomal, and six sex‐linked) for population samples of Lazuli Buntings (Passerina amoena) and Indigo Buntings (Passerina cyanea) (Aves: Cardinalidae) to test explicitly between a strict allopatric speciation model and a model in which divergence occurred despite postdivergence gene flow. Likelihood ratio tests of coalescent‐based population genetic parameter estimates indicated a strong signal of postdivergence gene flow and a strict allopatric speciation model was rejected. Analyses of partitioned datasets (mitochondrial, autosomal, and sex‐linked) suggest the overall gene flow patterns are driven primarily by autosomal gene flow, as there is no evidence of mitochondrial gene flow and we were unable to reject an allopatric speciation model for the sex‐linked data. This pattern is consistent with either a parapatric divergence model or repeated periods of allopatry with gene flow occurring via secondary contact. These results are consistent with the low fitness of female avian hybrids under Haldane's rule and demonstrate that sex‐linked loci likely are important in the initial generation of reproductive isolation, not just its maintenance.     

Fine-scale comparative phylogeography of a sympatric sister species triplet of subterranean diving beetles from a single calcrete aquifer in Western Australia

Calcrete aquifers in the arid Yilgarn region of central Western Australia are a biodiversity hotspot for stygofauna. A distinct pattern of interspecific size class variation among subterranean dytiscid beetle species has been observed in 29 of these aquifers where either two or three small, medium and/or large sympatric species are found that are in some cases sister species. We used a 3.5 km² grid of bores to sample dytiscids on a fine-scale and employed a comparative phylogeographical and population genetic approach to investigate the origins of a sympatric sister species triplet of diving beetles from a single aquifer. Mitochondrial DNA sequence data from the Cytochrome oxidase c subunit I gene revealed that all three species have high levels of haplotype diversity with ancient (∼1 million years ago) intra-specific coalescence of haplotypes, but low levels of nucleotide diversity. Population analyses provide evidence for multiple expansion events within each species. There was spatial heterogeneity in the distribution of genetic variation and abundance both within and among the three taxa. Population analyses revealed significant fine-scale differentiation with isolation by distance for Paroster macrosturtensis and P. mesosturtensis , but not the smallest species P. microsturtensis . Haplotype network analyses provided limited or no evidence for past population fragmentation within the large and small species, but substantial historical divergence was observed in P. mesosturtensis that was not spatially structured. A patchy population structure with contemporaneous and historical isolation by distance in the three species is likely to have been a significant isolating and diversifying force, preventing us from ruling out a potential role for allopatric divergence during speciation of this beetle sister triplet.     

Delimiting species without monophyletic gene trees

Genetic data are frequently used to delimit species, where species status is determined on the basis of an exclusivity criterium, such as reciprocal monophyly. Not only are there numerous empirical examples of incongruence between the boundaries inferred from such data compared to other sources like morphology -- especially with recently derived species, but population genetic theory also clearly shows that an inevitable bias in species status results because genetic thresholds do not explicitly take into account how the timing of speciation influences patterns of genetic differentiation. This study represents a fundamental shift in how genetic data might be used to delimit species. Rather than equating gene trees with a species tree or basing species status on some genetic threshold, the relationship between the gene trees and the species history is modeled probabilistically. Here we show that the same theory that is used to calculate the probability of reciprocal monophyly can also be used to delimit species despite widespread incomplete lineage sorting. The results from a preliminary simulation study suggest that very recently derived species can be accurately identified long before the requisite time for reciprocal monophyly to be achieved following speciation. The study also indicates the importance of sampling, both with regards to loci and individuals. Withstanding a thorough investigation into the conditions under which the coalescent-based approach will be effective, namely how the timing of divergence relative to the effective population size of species affects accurate species delimitation, the results are nevertheless consistent with other recent studies (aimed at inferring species relationships), showing that despite the lack of monophyletic gene trees, a signal of species divergence persists and can be extracted. Using an explicit model-based approach also avoids two primary problems with species delimitation that result when genetic thresholds are applied with genetic data -- the inherent biases in species detection arising from when and how speciation occurred, and failure to take into account the high stochastic variance of genetic processes. Both the utility and sensitivities of the coalescent-based approach outlined here are discussed; most notably, a model-based approach is essential for determining whether incompletely sorted gene lineages are (or are not) consistent with separate species lineages, and such inferences require accurate model parameterization (i.e., a range of realistic effective population sizes relative to potential times of divergence for the purported species). It is the goal (and motivation of this study) that genetic data might be used effectively as a source of complementation to other sources of data for diagnosing species, as opposed to the exclusion of other evidence for species delimitation, which will require an explicit consideration of the effects of the temporal dynamic of lineage splitting on genetic data.     

Perspective: gene divergence, population divergence, and the variance in coalescence time in phylogeographic studies

Molecular methods as applied to the biogeography of single species (phylogeography) or multiple codistributed species (comparative phylogeography) have been productively and extensively used to elucidate common historical features in the diversification of the Earth's biota. However, only recently have methods for estimating population divergence times or their confidence limits while taking into account the critical effects of genetic polymorphism in ancestral species become available, and earlier methods for doing so are underutilized. We review models that address the crucial distinction between the gene divergence, the parameter that is typically recovered in molecular phylogeographic studies, and the population divergence, which is in most cases the parameter of interest and will almost always postdate the gene divergence. Assuming that population sizes of ancestral species are distributed similarly to those of extant species, we show that phylogeographic studies in vertebrates suggest that divergence of alleles in ancestral species can comprise from less than 10% to over 50% of the total divergence between sister species, suggesting that the problem of ancestral polymorphism in dating population divergence can be substantial. The variance in the number of substitutions (among loci for a given species or among species for a given gene) resulting from the stochastic nature of DNA change is generally smaller than the variance due to substitutions along allelic lines whose coalescence times vary due to genetic drift in the ancestral population. Whereas the former variance can be reduced by further DNA sequencing at a single locus, the latter cannot. Contrary to phylogeographic intuition, dating population divergence times when allelic lines have achieved reciprocal monophyly is in some ways more challenging than when allelic lines have not achieved monophyly, because in the former case critical data on ancestral population size provided by residual ancestral polymorphism is lost. In the former case differences in coalescence time between species pairs can in principle be explained entirely by differences in ancestral population size without resorting to explanations involving differences in divergence time. Furthermore, the confidence limits on population divergence times are severely underestimated when those for number of substitutions per site in the DNA sequences examined are used as a proxy. This uncertainty highlights the importance of multilocus data in estimating population divergence times; multilocus data can in principle distinguish differences in coalescence time (T) resulting from differences in population divergence time and differences in T due to differences in ancestral population sizes and will reduce the confidence limits on the estimates. We analyze the contribution of ancestral population size (theta) to T and the effect of uncertainty in theta on estimates of population divergence (tau) for single loci under reciprocal monophyly using a simple Bayesian extension of Takahata and Satta's and Yang's recent coalescent methods. The confidence limits on tau decrease when the range over which ancestral population size theta is assumed to be distributed decreases and when tau increases; they generally exclude zero when tau/(4Ne) > 1. We also apply a maximum-likelihood method to several single and multilocus data sets. With multilocus data, the criterion for excluding tau = 0 is roughly that l tau/(4Ne) > 1, where l is the number of loci. Our analyses corroborate recent suggestions that increasing the number of loci is critical to decreasing the uncertainty in estimates of population divergence time.     

Investigating niche and lineage diversification in widely distributed taxa: phylogeography and ecological niche modeling of the Peromyscus maniculatus species group

The relationship between lineage formation and variation in the ecological niche is a fundamental evolutionary question. Two prevailing hypotheses reflect this relationship: niche conservatism and niche divergence. Niche conservatism predicts a pattern where sister taxa will occupy similar niche spaces; whereas niche divergence predicts that sister taxa will occupy different niche spaces. Widely distributed species often show distinct phylogeographic structure, but little research has been conducted on how the environment may be related to these phylogenetic patterns. We investigated the relationship between lineage divergence and environmental space for the closely related species Peromyscus maniculatus and P. polionotus utilizing phylogenetic techniques and ecological niche modeling (ENM). We estimated the phylogenetic relationship among individuals based on complete cytochrome b sequences that represent individuals from a majority of the species ranges. Niche spaces that lineages occupy were estimated by using 12 environmental layers. Differences in niche space were tested using multivariate statistics based on location data, and ENMs were employed using maximum entropy algorithms. Two similarity indices estimated significant divergence in environmental space based on the ENM. Six geographically structured lineages were identified within P. maniculatus. Nested within P. maniculatus we found that P. polionotus recently diverged from a clade occupying central and western United States. We estimated that the majority of the genetic lineages occupy distinct environmental niches, which supports a pattern of niche divergence. Two sister taxa showed niche divergence and represent different ecomorphs, suggesting morphological, genetic and ecological divergence between the two lineages. Two other sister taxa were observed in the same environmental space based on multivariate statistics, suggesting niche conservatism. Overall our results indicate that a widely distributed species may exhibit both niche conservatism and niche divergence, and that most lineages seem to occupy distinct environmental niches.     

从核rDNA序列探讨隙蛛亚科Coelotinae的分类地位

隙蛛亚科Coelotinae主要分布于东亚地区,其中我国的已有种类占到全世界种数的一半以上,因此对于我国隙蛛类蜘蛛的研究已经成为世界暗蛛科研究的重点之一。隙蛛亚科属于无筛器类群,于1893年,由Cambridge以隙蛛属为模式属而建立,归属于无筛器的漏斗蛛科。之后,虽然经历了数次修订     

Tree hybridization and genotypic variation drive cryptic speciation of a specialist mite herbivore

Few studies have investigated the roles that plant hybridization and individual plant genotype play in promoting population divergence within arthropod species. Using nrDNA sequence information and reciprocal transfer experiments, we examined how tree cross type (i.e., pure Populus angustifolia and P. angustifolia x P. fremontii F(1) type hybrids) and individual tree genotype influence host race formation in the bud-galling mite Aceria parapopuli. Three main findings emerged: (1) Strong genetic differentiation of mite populations found on pure P. angustifolia and F(1) type hybrids indicates that these mites represent morphologically cryptic species. (2) Within the F(1) type hybrids, population genetic analyses indicate migration among individual trees; however, (3) transfer experiments show that the mites found on heavily infested F(1) type trees perform best on their natal host genotype, suggesting that genetic interactions between mites and their host trees drive population structure, local adaptation, and host race formation. These findings argue that hybridization and genotypic differences in foundation tree species may drive herbivore population structure, and have evolutionary consequences for dependent arthropod species.     

Molecular systematics of the old world monkey tribe papionini: analysis of the total available genetic sequences

The phylogenetic relationships among the genera of the tribe Papionini are inferred using a taxonomic congruence approach in which gene trees derived for eight unlinked genetic sequence datasets are compared. Population genetics theory predicts that species relationships will be revealed with greater probability when the topology of gene trees from many unlinked loci are found to be congruent. The theory underlying this approach is described. Monophyly of the mangabeys is not supported by any of the gene trees; instead, they are polyphyletic with Cercocebus found to be the sister taxon to Mandrillus in five gene trees (with no conflicting trees), and Lophocebus found to be closely related to Papio and/or Theropithecus in all trees. Theropithecus and Papio are not strongly supported as sister taxa (present in one or two trees only);Lophocebus and Papio are supported as sister taxa in the majority of trees. A close relationship between Mandrillus and Papio is not supported in any of the trees.The relationships among Papio, Lophocebus, and Theropithecus cannot be resolved by congruence, probably due to the short time interval estimated between their divergences. The mtDNA COII sequences are used to estimate divergence dates within the papionins. The internode between the divergences of these species is estimated to be between 290 ka and 370 ka. Lastly, the evolution of morphological features such as long faces, suborbital facial fossae, and terrestrial skeletal adaptations is discussed.     

Phylogeny of the genera Entamoeba and Endolimax as deduced from small-subunit ribosomal RNA sequences

We sequenced small-subunit ribosomal RNA genes (16S-like rDNAs) of 10 species belonging to the genera Entamoeba and Endolimax. This study was undertaken to (1) resolve the relationships among the major lineages of Entamoeba previously identified by riboprinting; (2) examine the validity of grouping the genera Entamoeba and Endolimax in the same family, the Entamoebidae; and (3) examine how different models of nucleotide evolution influence the position of Entamoeba in eukaryotic phylogenetic reconstructions. The results obtained with distance, parsimony, and maximum-likelihood analyses support monophyly of the genus Entamoeba and are largely in accord with riboprinting results. Species of Entamoeba producing cysts with the same number of nuclei from monophyletic groups. The most basal Entamoeba species are those that produce cysts with eight nuclei, while the group producing four-nucleated cysts is most derived. Most phylogenetic reconstructions support monophyly of the Entamoebidae. In maximum-likelihood and parsimony analyses, Endolimax is a sister taxon to Entamoeba, while in some distance analyses, it represents a separate lineage. The secondary loss of mitochondria and other organelles from these genera is confirmed by their relatively late divergence in eukaryotic 16S-like rDNA phylogenies. Finally, we show that the positions of some (fast-evolving) eukaryotic lineages are uncertain in trees constructed with models that make corrections for among-site rate variation.     

Speciation‐by‐depth on coral reefs: Sympatric divergence with gene flow or cryptic transient isolation?

The distributions of many sister species in the sea overlap geographically but are partitioned along depth gradients. The genetic changes leading to depth segregation may evolve in geographic isolation as a prerequisite to coexistence or may emerge during primary divergence leading to new species. These alternatives can now be distinguished via the power endowed by the thousands of scorable loci provided by second‐generation sequence data. Here, we revisit the case of two depth‐segregated, genetically isolated ecotypes of the nominal Caribbean candelabrum coral Eunicea flexuosa. Previous analyses based on a handful of markers could not distinguish between models of genetic exchange after a period of isolation (consistent with secondary contact) and divergence with gene flow (consistent with primary divergence). Analyses of the history of isolation, genetic exchange and population size based on 15,640 new SNP markers derived from RNAseq data best support models where divergence began 800K BP and include epochs of divergence with gene flow, but with an intermediate period of transient isolation. Results also supported the previous conclusion that recent exchange between the ecotypes occurs asymmetrically from the Shallow lineage to the Deep. Parallel analyses of data from two other corals with depth‐segregated populations (Agaricia fragilis and Pocillopora damicornis) suggest divergence leading to depth‐segregated populations may begin with a period of symmetric exchange, but that an epoch of population isolation precedes more complete isolation marked by asymmetric introgression. Thus, while divergence‐with‐gene flow may account for much of the differentiation that separates closely related, depth‐segregated species, it remains to be seen whether any critical steps in the speciation process only occur when populations are isolated.     

Distribution modelling and statistical phylogeography: an integrative framework for generating and testing alternative biogeographical hypotheses

Statistical phylogeographic studies contribute to our understanding of the factors that influence population divergence and speciation, and that ultimately generate biogeographical patterns. The use of coalescent modelling for analyses of genetic data provides a framework for statistically testing alternative hypotheses about the timing and pattern of divergence. However, the extent to which such approaches contribute to our understanding of biogeography depends on how well the alternative hypotheses chosen capture relevant aspects of species histories. New modelling techniques, which explicitly incorporate spatio-geographic data external to the gene trees themselves, provide a means for generating realistic phylogeographic hypotheses, even for taxa without a detailed fossil record. Here we illustrate how two such techniques – species distribution modelling and its historical extension, palaeodistribution modelling – in conjunction with coalescent simulations can be used to generate and test alternative hypotheses. In doing so, we highlight a few key studies that have creatively integrated both historical geographic and genetic data and argue for the wider incorporation of such explicit integrations in biogeographical studies.     

Species trees from gene trees: reconstructing Bayesian posterior distributions of a species phylogeny using estimated gene tree distributions

The desire to infer the evolutionary history of a group of species should be more viable now that a considerable amount of multilocus molecular data is available. However, the current molecular phylogenetic paradigm still reconstructs gene trees to represent the species tree. Further, commonly used methods of combining data, such as the concatenation method, are known to be inconsistent in some circumstances. In this paper, we propose a Bayesian hierarchical model to estimate the phylogeny of a group of species using multiple estimated gene tree distributions, such as those that arise in a Bayesian analysis of DNA sequence data. Our model employs substitution models used in traditional phylogenetics but also uses coalescent theory to explain genealogical signals from species trees to gene trees and from gene trees to sequence data, thereby forming a complete stochastic model to estimate gene trees, species trees, ancestral population sizes, and species divergence times simultaneously. Our model is founded on the assumption that gene trees, even of unlinked loci, are correlated due to being derived from a single species tree and therefore should be estimated jointly. We apply the method to two multilocus data sets of DNA sequences. The estimates of the species tree topology and divergence times appear to be robust to the prior of the population size, whereas the estimates of effective population sizes are sensitive to the prior used in the analysis. These analyses also suggest that the model is superior to the concatenation method in fitting these data sets and thus provides a more realistic assessment of the variability in the distribution of the species tree that may have produced the molecular information at hand. Future improvements of our model and algorithm should include consideration of other factors that can cause discordance of gene trees and species trees, such as horizontal transfer or gene duplication.     

Global population genetic structure and biogeography of the oceanic copepods Eucalanus hyalinus and E. spinifer

Although theory dictates that limited gene flow between populations is a necessary precursor to speciation under allopatric and parapatric models, it is currently unclear how genetic differentiation between conspecific populations can arise in open-ocean plankton species. I examined two recently distinguished sympatric, circumglobal sister species, Eucalanus hyalinus and Eucalanus spinifer, for population genetic structure throughout their global biogeographic ranges. Here I show that oceanic zooplankton species can be highly genetically structured on macrogeographic spatial scales, despite experiencing extensive gene flow within features of the large-scale ocean circulation. Mitochondrial DNA analyses of 450 and 383 individuals of E. hyalinus and E. spinifer, respectively, revealed that habitat discontinuities at the boundaries of subtropical gyres in the North and South Pacific, as well as continental land masses, acted as effective barriers to gene flow for both species. However, the impact of specific barriers on population genetic structure varied between the sister species, despite their close phylogenetic relationship and similar circumglobal biogeogeographic distributions. The sister species differed in their oceanographic distributions, with E. spinifer dominating oligotrophic waters of the subtropical gyres and E. hyalinus more abundant along central water mass boundaries and in frontal zones and upwelling systems. This species-specific difference in the oceanographic habitat is an important factor determining the historical and contemporary patterns of dispersal of the two species. I suggest that species-specific ecological differences are likely to be a primary determinant of population genetic structure of open-ocean plankton.     

啄木鸟科内谱系生物地理学新探

啄木鸟群体共包含254种,但早期对啄木鸟的分类学研究局限于形态学、行为学和鸣唱法等方向,分子研究结果又彼此存在差异,因此啄木鸟分类系统目前还有争议。分子数据的大量累积为啄木鸟系统发育的研究提供了条件。利用从GenBank获取的ND2基因序列,对啄木鸟科179种啄木鸟采用NJ法、MP法以及ML法分别进行了分子进化树构建,同时用NJ法重建了Picumnus属分子进化树。此外又利用cytb基因序列采用NJ法对Dryobates属和Picoides属进行分子进化树重建。我们的结果显示,姬啄木亚科Picumnus属新大陆物种彼此的进化关系仍需进一步研究,而该属唯一一种旧大陆种P.innominatus,种内遗传差异在种间水平依然显著,这说明其种内遗传差异一直被低估了。在啄木鸟亚科内,Hemicircus属和Nesoctites属为姐妹关系,而Geocolaptes olivaceus插入了Campethera属支系中,并非与其构成姐妹群体。此外,一些姐妹群在新旧大陆上的分布,推测这种分布可能与第四纪冰川有关。啄木鸟从旧大陆向新大陆迁移,在第四纪冰川作用下不断发生分布范围的扩张与收缩,最终形成现存的分布格局。     

Speciation processes in putative island endemic sister bat species: false impressions from mitochondrial DNA and microsatellite data

Cases of geographically restricted co‐occurring sister taxa are rare and may point to potential divergence with gene flow. The two bat species Murina gracilis and Murina recondita are both endemic to Taiwan and are putative sister species. To test for nonallopatric divergence and gene flow in these taxa, we generated sequences using Sanger and next‐generation sequencing, and combined these with microsatellite data for coalescent‐based analyses. MtDNA phylogenies supported the reciprocally monophyletic sister relationship between M. gracilis and M. recondita; however, clustering of microsatellite genotypes revealed several cases of species admixture suggesting possible introgression. Sequencing of microsatellite flanking regions revealed that admixture signatures stemmed from microsatellite allele homoplasy rather than recent introgressive hybridization, and also uncovered an unexpected sister relationship between M. recondita and the continental species Murina eleryi, to the exclusion of M. gracilis. To dissect the basis of these conflicts between ncDNA and mtDNA, we analysed sequences from 10 anonymous ncDNA loci with *beast and isolation‐with‐migration and found two distinct clades of M. eleryi, one of which was sister to M. recondita. We conclude that Taiwan was colonized by the ancestor of M. gracilis first, followed by the ancestor of M. recondita after a period of allopatric divergence. After colonization, the mitochondrial genome of M. recondita was replaced by that of the resident M. gracilis. This study illustrates how apparent signatures of sympatric divergence can arise from complex histories of allopatric divergence, colonization and hybridization, thus highlighting the need for rigorous analyses to distinguish between such scenarios.