共查询到20条相似文献，搜索用时 15 毫秒
1.
Tom&#x; Flouri Jun Huang Xiyun Jiao Paschalia Kapli Bruce Rannala Ziheng Yang 《Molecular biology and evolution》2022,39(8)
The multispecies coalescent (MSC) model accommodates both species divergences and within-species coalescent and provides a natural framework for phylogenetic analysis of genomic data when the gene trees vary across the genome. The MSC model implemented in the program bpp assumes a molecular clock and the Jukes–Cantor model, and is suitable for analyzing genomic data from closely related species. Here we extend our implementation to more general substitution models and relaxed clocks to allow the rate to vary among species. The MSC-with-relaxed-clock model allows the estimation of species divergence times and ancestral population sizes using genomic sequences sampled from contemporary species when the strict clock assumption is violated, and provides a simulation framework for evaluating species tree estimation methods. We conducted simulations and analyzed two real datasets to evaluate the utility of the new models. We confirm that the clock-JC model is adequate for inference of shallow trees with closely related species, but it is important to account for clock violation for distant species. Our simulation suggests that there is valuable phylogenetic information in the gene-tree branch lengths even if the molecular clock assumption is seriously violated, and the relaxed-clock models implemented in bpp are able to extract such information. Our Markov chain Monte Carlo algorithms suffer from mixing problems when used for species tree estimation under the relaxed clock and we discuss possible improvements. We conclude that the new models are currently most effective for estimating population parameters such as species divergence times when the species tree is fixed. 相似文献
2.
Joly S 《Molecular ecology resources》2012,12(1):179-184
I introduce the software JML that tests for the presence of hybridization in multispecies sequence data sets by posterior predictive checking following Joly, McLenachan and Lockhart (2009, American Naturalist 174, e54). Although their method could potentially be applied on any data set, the lack of appropriate software made its application difficult. The software JML thus fills a need for an easy application of the method but also includes improvements such as the possibility to incorporate uncertainty in the species tree topology. The JML software uses a posterior distribution of species trees, population sizes and branch lengths to simulate replicate sequence data sets using the coalescent with no migration. A test quantity, defined as the minimum pairwise sequence distance between sequences of two species, is then evaluated on the simulated data sets and compared to the one estimated from the original data. Because the test quantity is a good predictor of hybridization events, departure from the bifurcating species tree model could be interpreted as evidence of hybridization. Software performance in terms of computing time is evaluated for several parameters. I also show an application example of the software for detecting hybridization among native diploid North American roses. 相似文献
3.
Interest in methods that estimate speciation and extinction rates from molecular phylogenies has increased over the last decade. The application of such methods requires reliable estimates of tree topology and node ages, which are frequently obtained using standard phylogenetic inference combining concatenated loci and molecular dating. However, this practice disregards population‐level processes that generate gene tree/species tree discordance. We evaluated the impact of employing concatenation and coalescent‐based phylogeny inference in recovering the correct macroevolutionary regime using simulated data based on the well‐established diversification rate shift of delphinids in Cetacea. We found that under scenarios of strong incomplete lineage sorting, macroevolutionary analysis of phylogenies inferred by concatenating loci failed to recover the delphinid diversification shift, while the coalescent‐based tree consistently retrieved the correct rate regime. We suggest that ignoring microevolutionary processes reduces the power of methods that estimate macroevolutionary regimes from molecular data. 相似文献
4.
Jun Huang Yuttapong Thawornwattana Tom&#x; Flouri James Mallet Ziheng Yang 《Molecular biology and evolution》2022,39(12)
Genomic sequence data provide a rich source of information about the history of species divergence and interspecific hybridization or introgression. Despite recent advances in genomics and statistical methods, it remains challenging to infer gene flow, and as a result, one may have to estimate introgression rates and times under misspecified models. Here we use mathematical analysis and computer simulation to examine estimation bias and issues of interpretation when the model of gene flow is misspecified in analysis of genomic datasets, for example, if introgression is assigned to the wrong lineages. In the case of two species, we establish a correspondence between the migration rate in the continuous migration model and the introgression probability in the introgression model. When gene flow occurs continuously through time but in the analysis is assumed to occur at a fixed time point, common evolutionary parameters such as species divergence times are surprisingly well estimated. However, the time of introgression tends to be estimated towards the recent end of the period of continuous gene flow. When introgression events are assigned incorrectly to the parental or daughter lineages, introgression times tend to collapse onto species divergence times, with introgression probabilities underestimated. Overall, our analyses suggest that the simple introgression model is useful for extracting information concerning between-specific gene flow and divergence even when the model may be misspecified. However, for reliable inference of gene flow it is important to include multiple samples per species, in particular, from hybridizing species. 相似文献
5.
Michael Gruenstaeudl Noah M. Reid Gregory L. Wheeler Bryan C. Carstens 《Molecular ecology resources》2016,16(1):193-205
Bayesian inference operates under the assumption that the empirical data are a good statistical fit to the analytical model, but this assumption can be challenging to evaluate. Here, we introduce a novel r package that utilizes posterior predictive simulation to evaluate the fit of the multispecies coalescent model used to estimate species trees. We conduct a simulation study to evaluate the consistency of different summary statistics in comparing posterior and posterior predictive distributions, the use of simulation replication in reducing error rates and the utility of parallel process invocation towards improving computation times. We also test P2C2M on two empirical data sets in which hybridization and gene flow are suspected of contributing to shared polymorphism, which is in violation with the coalescent model: Tamias chipmunks and Myotis bats. Our results indicate that (i) probability‐based summary statistics display the lowest error rates, (ii) the implementation of simulation replication decreases the rate of type II errors, and (iii) our r package displays improved statistical power compared to previous implementations of this approach. When probabilistic summary statistics are used, P2C2M corroborates the assumption that genealogies collected from Tamias and Myotis are not a good fit to the multispecies coalescent model. Taken as a whole, our findings argue that an assessment of the fit of the multispecies coalescent model should accompany any phylogenetic analysis that estimates a species tree. 相似文献
6.
Oaks JR 《Evolution; international journal of organic evolution》2011,65(11):3285-3297
True crocodiles (Crocodylus) are the most broadly distributed, ecologically diverse, and species-rich crocodylian genus, comprising about half of extant crocodylian diversity and exhibiting a circumtropical distribution. Crocodylus traditionally has been viewed as an ancient group of morphologically conserved species that originated in Africa prior to continental breakup. In this study, these long-held notions about the temporal and geographic origin of Crocodylus are tested using DNA sequence data of 10 loci from 76 individuals representing all 23 crocodylian species. I infer a time-calibrated species tree of all Crocodylia and estimate the spatial pattern of diversification within Crocodylus. For the first time, a fully resolved phylogenetic estimate of all Crocodylia is well-supported. The results overturn traditional views of the evolution of Crocodylus by demonstrating that the true crocodiles are not \"living-fossils\" that originated in Africa. Rather, Crocodylus originated from an ancestor in the tropics of the Late Miocene Indo-Pacific, and rapidly radiated and dispersed around the globe during a period marked by mass extinctions of fellow crocodylians. The findings also reveal more diversity within the genus than is recognized by current taxonomy. 相似文献
8.
Carlos G. Schrago Beatriz Mello Anieli G. Pereira Carolina Furtado Hector N. Seuánez 《Ecology and evolution》2018,8(2):1206-1216
Multispecies coalescent (MSC) theory assumes that gene trees inferred from individual loci are independent trials of the MSC process. As genes might be physically close in syntenic associations spanning along chromosome regions, these assumptions might be flawed in evolutionary lineages with substantial karyotypic shuffling. Neotropical primates (NP) represent an ideal case for assessing the performance of MSC methods in such scenarios because chromosome diploid number varies significantly in this lineage. To this end, we investigated the effect of sequence length on the theoretical expectations of MSC model, as well as the results of coalescent‐based tree inference methods. This was carried out by comparing NP with hominids, a lineage in which chromosome macrostructure has been stable for at least 15 million years. We found that departure from the MSC model in Neotropical primates decreased with smaller sequence fragments, where sites sharing the same evolutionary history were more frequently found than in longer fragments. This scenario probably resulted from extensive karyotypic rearrangement occurring during the radiation of NP, contrary to the comparatively stable chromosome evolution in hominids. 相似文献
9.
10.
11.
The multispecies coalescent model provides a natural framework for species tree estimation accounting for gene-tree conflicts. Although a number of species tree methods under the multispecies coalescent have been suggested and evaluated using simulation, their statistical properties remain poorly understood. Here, we use mathematical analysis aided by computer simulation to examine the identifiability, consistency, and efficiency of different species tree methods in the case of three species and three sequences under the molecular clock. We consider four major species-tree methods including concatenation, two-step, independent-sites maximum likelihood, and maximum likelihood. We develop approximations that predict that the probit transform of the species tree estimation error decreases linearly with the square root of the number of loci. Even in this simplest case, major differences exist among the methods. Full-likelihood methods are considerably more efficient than summary methods such as concatenation and two-step. They also provide estimates of important parameters such as species divergence times and ancestral population sizes,whereas these parameters are not identifiable by summary methods. Our results highlight the need to improve the statistical efficiency of summary methods and the computational efficiency of full likelihood methods of species tree estimation. 相似文献
12.
Bryan C. Carstens Michael Gruenstaeudl Noah M. Reid 《Evolution; international journal of organic evolution》2016,70(5):1080-1093
Groups of codistributed species that responded in a concerted manner to environmental events are expected to share patterns of evolutionary diversification. However, the identification of such groups has largely been based on qualitative, post hoc analyses. We develop here two methods (posterior predictive simulation [PPS], Kuhner–Felsenstein [K–F] analysis of variance [ANOVA]) for the analysis of codistributed species that, given a group of species with a shared pattern of diversification, allow empiricists to identify those taxa that do not codiversify (i.e., “outlier” species). The identification of outlier species makes it possible to jointly estimate the evolutionary history of co‐diversifying taxa. To evaluate the approaches presented here, we collected data from Páramo dipterans, identified outlier species, and estimated a “community tree” from species that are identified as having codiversified. Our results demonstrate that dipteran communities from different Páramo habitats in the same mountain range are more closely related than communities in other ranges. We also conduct simulation testing to evaluate this approach. Results suggest that our approach provides a useful addition to comparative phylogeographic methods, while identifying aspects of the analysis that require careful interpretation. In particular, both the PPS and K–F ANOVA perform acceptably when there are one or two outlier species, but less so as the number of outliers increases. This is likely a function of the corresponding degradation of the signal of community divergence; without a strong signal from a codiversifying community, there is no dominant pattern from which to detect an outlier species. For this reason, both the magnitude of K–F distance distribution and outside knowledge about the phylogeographic history of each putative member of the community should be considered when interpreting the results. 相似文献
13.
Accurate species delimitation is critical for biodiversity studies. However, species complexes characterized by introgression, high levels of population structure and subtle phenotypic differentiation can be challenging to delimit. Here, we report on a molecular systematic investigation of the woodland salamanders Plethodon wehrlei and Plethodon punctatus, which traditionally have been placed in the Plethodon wehrlei species group. To quantify patterns of genetic variation, we collected genetic samples from throughout the range of both species, including 22 individuals from nine populations of P. punctatus, and 60 individuals from 26 populations of P. wehrlei. From these samples, we sequenced three mtDNA loci (5596 base pairs) and five nuclear loci (3377 base pairs). We inferred time‐calibrated gene trees and species trees using BEAST 2.4.6, and we delimited putative species using a Bayesian implementation of the general mixed Yule‐coalescent model (bGMYC) and STRUCTURE. Finally, we validated putative species using the multispecies coalescent as implemented in Bayesian Phylogenetics and Phylogeography (BPP). We found substantial phylogeographic diversity in P. wehrlei, including multiple geographically cohesive clades and an inferred mitochondrial common ancestor at 11.5 myr (95% HPD: 9.6–13.6 myr) that separated populations formerly assigned to P. dixi from all other populations. We also found that P. punctatus is deeply nested within P. wehrlei, rendering the latter paraphyletic. After discussing the challenges faced by modern species delimitation methods, we recommend retaining P. punctatus because it is ecologically and phenotypically distinct. We further recommend that P. dixi be recognized as a valid species. 相似文献
14.
A multispecies-based taxonomic microarray targeting coding sequences of diverged orthologous genes in Saccharomyces cerevisiae, Saccharomyces paradoxus, Saccharomyces mikatae, Saccharomyces bayanus, Saccharomyces kudriavzevii, Naumovia castellii, Lachancea kluyveri and Candida glabrata was designed to allow identification of isolates of these species and their interspecies hybrids. Analysis of isolates of several Saccharomyces species and interspecies hybrids demonstrated the ability of the microarray to differentiate these yeasts on the basis of their specific hybridization patterns. Subsequent analysis of 183 supposed S. cerevisiae isolates of various ecological and geographical backgrounds revealed one misclassified S. bayanus or Saccharomyces uvarum isolate and four aneuploid interspecies hybrids, one between S. cerevisiae and S. bayanus and three between S. cerevisiae and S. kudriavzevii . Furthermore, this microarray design allowed the detection of multiple introgressed S. paradoxus DNA fragments in the genomes of three different S. cerevisiae isolates. These results show the power of multispecies-based microarrays as taxonomic tools for the identification of species and interspecies hybrids, and their ability to provide a more detailed characterization of interspecies hybrids and recombinants. 相似文献
15.
Accurate species delimitation is the key to precise estimation of species diversity and is fundamental to most branches of biology. Unclear species boundaries within species complexes could lead to the underestimation of species diversity. However, species delimitation of species complexes remains challenging due to the continuum of phenotypic variations. To robustly examine species boundaries within a species complex, integrative approaches in phylogeny, ecology, and morphology were applied to the Stewartia sinensis complex (Theaceae) endemic to China. Multispecies coalescent-based species delimitation using 572 nuclear ortholog sequences (anchored enrichment) supported reciprocal phylogenetic monophyly of the northern lineage (NL) and southern lineage (SL), which were not sister clades. Niche equivalency and similarity tests demonstrated significant climatic niche differentiation between NL and SL with observed Warren et al.'s I = 0.0073 and Schoener's D = 0.0021. Species distribution modeling also separated their potential distribution. Morphometric analyses suggested significant interlineage differentiation of multiple traits including the ratio of length and width, leaf width, and pedicel length, although overall similarity did not differ. Based on the integrative species concept, two distinct species were proposed with legitimate names of Stewartia gemmata for SL and S. sinensis for NL. Our empirical study of the S. sinensis complex highlights the importance of applying multiple species criteria, in particular the underappreciated niche differentiation, to species delimitation in species complexes pervasive in plants. 相似文献
16.
İsmaİl K. Sağlam Jason Baumsteiger Matt J. Smith Javier Linares‐Casenave Andrew L. Nichols Sean M. O'Rourke Michael R. Miller 《Molecular ecology》2016,25(16):3962-3973
The Devils Hole pupfish (Cyprinodon diabolis; DHP) is an icon of conservation biology. Isolated in a 50 m2 pool (Devils Hole), DHP is one of the rarest vertebrate species known and an evolutionary anomaly, having survived in complete isolation for thousands of years. However, recent findings suggest DHP might be younger than commonly thought, potentially introduced to Devils Hole by humans in the past thousand years. As a result, the significance of DHP from an evolutionary and conservation perspective has been questioned. Here we present a high‐resolution genomic analysis of DHP and two closely related species, with the goal of thoroughly examining the temporal divergence of DHP. To this end, we inferred the evolutionary history of DHP from multiple random genomic subsets and evaluated four historical scenarios using the multispecies coalescent. Our results provide substantial information regarding the evolutionary history of DHP. Genomic patterns of secondary contact present strong evidence that DHP were isolated in Devils Hole prior to 20–10 ka and the model best supported by geological history and known mutation rates predicts DHP diverged around 60 ka, approximately the same time Devils Hole opened to the surface. We make the novel prediction that DHP colonized and have survived in Devils Hole since the cavern opened, and the two events (colonization and collapse of the cavern's roof) were caused by a common geologic event. Our results emphasize the power of evolutionary theory as a predictive framework and reaffirm DHP as an important evolutionary novelty, worthy of continued conservation and exploration. 相似文献
17.
Ulrike Obertegger Adam Cieplinski Diego Fontaneto Spiros Papakostas 《Zoologica scripta》2018,47(1):122-132
Discordance between mitochondrial and nuclear phylogenies is being increasingly recognized in animals and may confound DNA‐based taxonomy. This is especially relevant for taxa whose microscopic size often challenges any effort to distinguish between cryptic species without the assistance of molecular data. Regarding mitonuclear discordance, two strikingly contrasting scenarios have been recently demonstrated in the monogonont rotifers of the genus Brachionus. While strict mitonuclear concordance was observed in the marine B. plicatilis species complex, widespread hybridization‐driven mitonuclear discordance was revealed in the freshwater B. calyciflorus species complex. Here, we investigated the frequency of occurrence and the potential drivers of mitonuclear discordance in three additional freshwater monogonont rotifer taxa, and assessed its potential impact on the reliability of DNA taxonomy results based on commonly used single markers. We studied the cryptic species complexes of Keratella cochlearis, Polyarthra dolichoptera and Synchaeta pectinata. Phylogenetic reconstructions were based on the mitochondrial barcoding marker cytochrome c oxidase subunit I gene and the nuclear internal transcribed spacer 1 locus, which currently represent the two most typical genetic markers used in rotifer DNA taxonomy. Species were delimited according to each marker separately using a combination of tree‐based coalescent, distance‐based and allele‐sharing‐based approaches. Mitonuclear discordance was observed in all species complexes with incomplete lineage sorting and unresolved phylogenetic reconstructions recognized as the likely drivers. Evidence from additional sources, such as morphology and ecology, is thus advisable for deciding between often contrasting mitochondrial and nuclear species scenarios in these organisms. 相似文献
18.
Secondary contact in close relatives can result in hybridization and the admixture of previously isolated gene pools. However, after an initial period of hybridization, reproductive isolation can evolve through different processes and lead to the interruption of gene flow and the completion of the speciation process. Omocestus minutissimus and O. uhagonii are two closely related grasshoppers with partially overlapping distributions in the Central System mountains of the Iberian Peninsula. To analyse spatial patterns of historical and/or contemporary hybridization between these two taxa and understand how species boundaries are maintained in the region of secondary contact, we sampled sympatric and allopatric populations of the two species and obtained genome‐wide single nucleotide polymorphism data using a restriction site‐associated DNA sequencing approach. We used Bayesian clustering analyses to test the hypothesis of contemporary hybridization in sympatric populations and employed a suite of phylogenomic approaches and a coalescent‐based simulation framework to evaluate alternative hypothetical scenarios of interspecific gene flow. Our analyses rejected the hypothesis of contemporary hybridization but revealed past introgression in the area where the distributions of the two species overlap. Overall, these results point to a scenario of historical gene flow after secondary contact followed by the evolution of reproductive isolation that currently prevents hybridization among sympatric populations. 相似文献
19.
Crispin Y. Jordan Konrad Lohse Frances Turner Marian Thomson Karim Gharbi Richard A. Ennos 《Molecular ecology》2018,27(5):1214-1228
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. 相似文献
20.
Florian Wagner Sabine Härtl Robert Vogt Christoph Oberprieler 《Molecular ecology》2017,26(16):4260-4283