Estimating species phylogeny from gene-tree probabilities despite incomplete lineage sorting: an example from Melanoplus grasshoppers

Estimating phylogenetic relationships among closely related species can be extremely difficult when there is incongruence among gene trees and between the gene trees and the species tree. Here we show that incorporating a model of the stochastic loss of gene lineages by genetic drift into the phylogenetic estimation procedure can provide a robust estimate of species relationships, despite widespread incomplete sorting of ancestral polymorphism. This approach is applied to a group of montane Melanoplus grasshoppers for which genealogical discordance among loci and incomplete lineage sorting obscures any obvious phylogenetic relationships among species. Unlike traditional treatments where gene trees estimated using standard phylogenetic methods are implicitly equated with the species tree, with the coalescent-based approach the species tree is modeled probabilistically from the estimated gene trees. The estimated species phylogeny (the ESP) is calculated for the grasshoppers from multiple gene trees reconstructed for nuclear loci and a mitochondrial gene. This empirical application is coupled with a simulation study to explore the performance of the coalescent-based approach. Specifically, we test the accuracy of the ESP given the data based on analyses of simulated data matching the multilocus data collected in Melanoplus (i.e., data were simulated for each locus with the same number of base pairs and locus-specific mutational models). The results of the study show that ESPs can be computed using the coalescent-based approach long before reciprocal monophyly has been achieved, and that these statistical estimates are accurate. This contrasts with analyses of the empirical data collected in Melanoplus and simulated data based on concatenation of multiple loci, for which the incomplete lineage sorting of recently diverged species posed significant problems. The strengths and potential challenges associated with incorporating an explicit model of gene-lineage coalescence into the phylogenetic procedure to obtain an ESP, as illustrated by application to Melanoplus, versus concatenation and consensus approaches are discussed. This study represents a fundamental shift in how species relationships are estimated - the relationship between the gene trees and the species phylogeny is modeled probabilistically rather than equating gene trees with a species tree.     

Discordance between phylogenetics and coalescent-based divergence modelling: exploring phylogeographic patterns of speciation in the Carex macrocephala species complex

We fit a molecular data set, consisting of the rpL16 cpDNA marker and eight microsatellite loci, to the isolation-with-migration model as implemented in IM a to test a well-supported phylogenetic hypothesis of relationships within the Carex macrocephala species complex (Cyperaceae). The phylogenetic hypothesis suggests C. macrocephala from North America is reciprocally monophyletic and is sister to a reciprocally monophyletic clade of C. kobomugi . The North American C. macrocephala and C. kobomugi clade form a sister clade with a lineage of Asian C. macrocephala , thereby forming a paraphyletic C. macrocephala species. Not only does the phylogenetic hypothesis suggest C. macrocephala is paraphyletic, but it also suggests that the two lineages which share a partially overlapping distribution, Asian C. macrocephala and C. kobomugi , are not the most closely related. To test these relationships, we used coalescent-based population genetic models to infer divergence time for each lineage pair within the species complex. The coalescent-based models account for the stochastic forces which drive population divergence, and can account for the lineage sorting that occurs prior to lineage divergence. A drawback to phylogenetic-based phylogeographical analyses is that they do not account for stochastic lineage sorting that occurs between gene divergence and lineage divergence. By comparing the relative divergence time of the three main lineages within this group, Asian C. macrocephala , North American C. macrocephala , and C. kobomugi , we concluded that the phylogenetic hypothesis is incorrect, and the divergence between these lineages occurred during the Late Pleistocene epoch.     

Species tree of a recent radiation: the subfamily Delphininae (Cetacea, Mammalia)

Lineages undergoing rapid radiations provide exceptional opportunities for studying speciation and adaptation, but also represent a challenge for molecular systematics because retention of ancestral polymorphisms and the occurrence of hybridization can obscure relationships among lineages. Dolphins in the subfamily Delphininae are one such case. Non-monophyly, rapid speciation events, and discordance between morphological and molecular characters have made the inference of phylogenetic relationships within this subfamily very difficult. Here we approach this problem by applying multiple methods intended to estimate species trees using a multi-gene dataset for the Delphininae (Sousa, Sotalia, Stenella, Tursiops, Delphinus and Lagenodelphis). Incongruent gene trees obtained indicate that incomplete lineage sorting and possibly hybridization are confounding the inference of species history in this group. Nonetheless, using coalescent-based methods, we have been able to extract an underlying species-tree signal from divergent histories of independent genes. This is the first time a molecular study provides support for such relationships. This study further illustrates how methods of species-tree inference can be very sensitive both to the characteristics of the dataset and the evolutionary processes affecting the evolution of the group under study.     

Inferring phylogeny despite incomplete lineage sorting

It is now well known that incomplete lineage sorting can cause serious difficulties for phylogenetic inference, but little attention has been paid to methods that attempt to overcome these difficulties by explicitly considering the processes that produce them. Here we explore approaches to phylogenetic inference designed to consider retention and sorting of ancestral polymorphism. We examine how the reconstructability of a species (or population) phylogeny is affected by (a) the number of loci used to estimate the phylogeny and (b) the number of individuals sampled per species. Even in difficult cases with considerable incomplete lineage sorting (times between divergences less than 1 N(e) generations), we found the reconstructed species trees matched the "true" species trees in at least three out of five partitions, as long as a reasonable number of individuals per species were sampled. We also studied the tradeoff between sampling more loci versus more individuals. Although increasing the number of loci gives more accurate trees for a given sampling effort with deeper species trees (e.g., total depth of 10 N(e) generations), sampling more individuals often gives better results than sampling more loci with shallower species trees (e.g., depth = 1 N(e)). Taken together, these results demonstrate that gene sequences retain enough signal to achieve an accurate estimate of phylogeny despite widespread incomplete lineage sorting. Continued improvement in our methods to reconstruct phylogeny near the species level will require a shift to a compound model that considers not only nucleotide or character state substitutions, but also the population genetics processes of lineage sorting. [Coalescence; divergence; population; speciation.].     

Phylogeny of Campanuloideae (Campanulaceae) with Emphasis on the Utility of Nuclear Pentatricopeptide Repeat (PPR) Genes

Background

The Campanuloideae (Campanulaceae) are a highly diverse clade of angiosperms found mostly in the Northern Hemisphere, with the highest diversity in temperate areas of the Old World. Chloroplast markers have greatly improved our understanding of this clade but many relationships remain unclear primarily due to low levels of molecular evolution and recent and rapid divergence. Furthermore, focusing solely on maternally inherited markers such as those from the chloroplast genome may obscure processes such as hybridization. In this study we explore the phylogenetic utility of two low-copy nuclear loci from the pentatricopeptide repeat gene family (PPR). Rapidly evolving nuclear loci may provide increased phylogenetic resolution in clades containing recently diverged or closely related taxa. We present results based on both chloroplast and low-copy nuclear loci and discuss the utility of such markers to resolve evolutionary relationships and infer hybridization events within the Campanuloideae clade.

Results

The inclusion of low-copy nuclear genes into the analyses provides increased phylogenetic resolution in two species-rich clades containing recently diverged taxa. We also obtain support for the placement of two early diverging lineages (Jasione and Musschia-Gadellia clades) that have previously been unresolved. Furthermore, phylogenetic analyses of PPR loci revealed potential hybridization events for a number of taxa (e.g., Campanula pelviformis and Legousia species). These loci offer greater overall topological support than obtained with plastid DNA alone.

Conclusion

This study represents the first inclusion of low-copy nuclear genes for phylogenetic reconstruction in Campanuloideae. The two PPR loci were easy to sequence, required no cloning, and the sequence alignments were straightforward across the entire Campanuloideae clade. Although potentially complicated by incomplete lineage sorting, these markers proved useful for understanding the processes of reticulate evolution and resolving relationships at a wide range of phylogenetic levels. Our results stress the importance of including multiple, independent loci in phylogenetic analyses.     

Molecular phylogeny of extant equids and effects of ancestral polymorphism in resolving species-level phylogenies

Short divergence times and processes such as incomplete lineage sorting and species hybridization are known to hinder the inference of species-level phylogenies due to the lack of sufficient informative genetic variation or the presence of shared but incongruent polymorphism among taxa. Extant equids (horses, zebras, and asses) are an example of a recently evolved group of mammals with an unresolved phylogeny, despite a large number of molecular studies. Previous surveys have proposed trees with rather poorly supported nodes, and the bias caused by genetic introgression or ancestral polymorphism has not been assessed. Here we studied the phylogenetic relationships of all extant species of Equidae by analyzing 22 partial mitochondrial and nuclear genes using maximum likelihood and Bayesian inferences that account for heterogeneous gene histories. We also examined genetic signatures of lineage sorting and/or genetic introgression in zebras by evaluating patterns of intraspecific genetic variation. Our study improved the resolution and support of the Equus phylogeny and in particular the controversial positions of the African wild ass (E. asinus) and mountain zebra (E. zebra): the African wild ass is placed as a sister species of the Asiatic asses and the mountain zebra as the sister taxon of Grevy's and Burchell's zebras. A shared polymorphism (indel) detected among zebra species in the Estrogen receptor 1 gene was likely due to incomplete lineage sorting and not genetic introgression as also indicated by other mitochondrial (Cytochrome b) and nuclear (Y chromosome and microsatellites) markers. Ancestral polymorphism in equids might have contributed to the long-standing lack of clarity in the phylogeny of this highly threatened group of mammals.     

A phylogenomic perspective on gene tree conflict and character evolution in Caprifoliaceae using target enrichment data,with Zabelioideae recognized as a new subfamily

The use of diverse data sets in phylogenetic studies aiming for understanding evolutionary histories of species can yield conflicting inference. Phylogenetic conflicts observed in animal and plant systems have often been explained by hybridization, incomplete lineage sorting (ILS), or horizontal gene transfer. Here, we used target enrichment data, species tree, and species network approaches to infer the backbone phylogeny of the family Caprifoliaceae, while distinguishing among sources of incongruence. We used 713 nuclear loci and 46 complete plastome sequence data from 43 samples representing 38 species from all major clades to reconstruct the phylogeny of the family using concatenation and coalescence approaches. We found significant nuclear gene tree conflict as well as cytonuclear discordance. Additionally, coalescent simulations and phylogenetic species network analyses suggested putative ancient hybridization among subfamilies of Caprifoliaceae, which seems to be the main source of phylogenetic discordance. Ancestral state reconstruction of six morphological characters revealed some homoplasy for each character examined. By dating the branching events, we inferred the origin of Caprifoliaceae at approximately 66.65 Ma in the late Cretaceous. By integrating evidence from molecular phylogeny, divergence times, and morphology, we here recognize Zabelioideae as a new subfamily in Caprifoliaceae. This work shows the necessity of using a combination of multiple approaches to identify the sources of gene tree discordance. Our study also highlights the importance of using data from both nuclear and plastid genomes to reconstruct deep and shallow phylogenies of plants.     

A multilocus phylogeny of the Sulidae (Aves: Pelecaniformes)

Gene trees will often differ from the true species history, the species tree, as a result of processes such as incomplete lineage sorting. New methods such as Bayesian Estimation of the Species Tree (BEST) use the multispecies coalescent to model lineage sorting, and directly infer the species tree from multilocus DNA sequence data. The Sulidae (Aves: Pelecaniformes) is a family of ten booby and gannet species with a global distribution. We sequenced five nuclear intron loci and one mitochondrial locus to estimate a species tree for the Sulidae using both BEST and by concatenating nuclear loci. We also used fossil calibrated strict and relaxed molecular clocks in BEAST to estimate divergence times for major nodes in the sulid phylogeny. Individual gene trees showed little phylogenetic conflict but varied in resolution. With the exception of the mitochondrial gene tree, no gene tree was completely resolved. On the other hand, both the BEST and concatenated species trees were highly resolved, strongly supported, and topologically consistent with each other. The three sulid genera (Morus, Sula, Papasula) were monophyletic and the relationships within genera were mostly consistent with both a previously estimated mtDNA gene tree and the mtDNA gene tree estimated here. However, our species trees conflicted with the mtDNA gene trees in the relationships among the three genera. Most notably, we find that the endemic and endangered Abbott's booby (Papasula abbotti) is likely basal to all other members of the Sulidae and diverged from them approximately 22 million years ago.     

PHYLOGEOGRAPHY OF THE JUMPING SPIDER HABRONATTUS PUGILLIS (ARANEAE: SALTICIDAE): RECENT VICARIANCE OF SKY ISLAND POPULATIONS?

In island systems with diverging populations, the history of island formation and genealogical estimates of divergence dates can be mutually informative. In the "sky islands" of southeastern Arizona, climate-induced contraction of woodlands appears to have fragmented populations of woodland-dwelling species onto disjunct mountain ranges. Montane populations of the jumping spider, Habronattus pugillis, display striking amounts of phenotypic divergence among ranges. Paleoclimatic estimates date woodland fragmentation at approximately 10,000 years ago, suggesting that phenotypic divergence has been extraordinarily rapid in these spiders. This phylogeographic study of populations of H. pugillis attempts to clarify the species' history of isolation and divergence and to address the suitability of available paleoclimatic data for dating divergences among populations of the region's woodland-dwelling organisms. Mitochondrial sequence data of spiders from 13 mountain ranges was used to reconstruct genealogical relationships. Gene trees show that small mountain ranges tend to have populations whose sequences form monophyletic groups, whereas larger ranges do not. Paraphyly among genes from larger ranges could result from either recent migration or incomplete lineage sorting. I use phylogenetic and geographic information to test these alternatives, and conclude that incomplete lineage sorting best explains the observed paraphyly. Gene trees are concordant with some of the predictions of vegetation history generated by examination of topography. Dates estimated for divergence of populations vary from 30,000 years to more than 2 million years ago, suggesting multiple vicariance events that are older than would be inferred from paleoclimatic studies. These findings illustrate that use of any single paleontological dataset to calibrate molecular clocks can potentially greatly underestimate actual divergence times.     

Phylogenetic relationships and ancient incomplete lineage sorting among cichlid fishes in Lake Tanganyika as revealed by analysis of the insertion of retroposons

Lake Tanganyika harbors numerous endemic species of extremely diverse cichlid fish that have been classified into 12 major taxonomic groups known as tribes. Analysis of short interspersed element (SINE) insertion data has been acknowledged to be a powerful tool for the elucidation of phylogenetic relationships, and we applied this method in an attempt to clarify such relationships among these cichlids. We studied insertion patterns of 38 SINEs in total, 24 of which supported the monophyly of three clades. The other 14 loci revealed extensive incongruence in terms of the patterns of SINE insertions. These incongruencies most likely stem from a period of adaptive radiation. One possible explanation for this phenomenon is the extensive incomplete lineage sorting of alleles for the presence or absence of a SINE during successive speciation events which took place about 5-10 MYA. The present study is the first to report the successful application of the SINE method in demonstrating the existence of such possible "ancient" incomplete lineage sorting. We discuss the possibility that it might potentially be very difficult to resolve the species phylogeny of a group that radiated explosively, even by resolving the genealogies of more than 10 nuclear loci, as a consequence of incomplete lineage sorting during speciation.     

Is RAD‐seq suitable for phylogenetic inference? An in silico assessment and optimization

Inferring phylogenetic relationships between closely related taxa can be hindered by three factors: (1) the lack of informative molecular variation at short evolutionary timescale; (2) the lack of established markers in poorly studied taxa; and (3) the potential phylogenetic conflicts among different genomic regions due to incomplete lineage sorting or introgression. In this context, Restriction site Associated DNA sequencing (RAD‐seq) seems promising as this technique can generate sequence data from numerous DNA fragments scattered throughout the genome, from a large number of samples, and without preliminary knowledge on the taxa under study. However, divergence beyond the within‐species level will necessarily reduce the number of conserved and non‐duplicated restriction sites, and therefore the number of loci usable for phylogenetic inference. Here, we assess the suitability of RAD‐seq for phylogeny using a simulated experiment on the 12 Drosophila genomes, with divergence times ranging from 5 to 63 million years. These simulations show that RAD‐seq allows the recovery of the known Drosophila phylogeny with strong statistical support, even for relatively ancient nodes. Notably, this conclusion is robust to the potentially confounding effects of sequencing errors, heterozygosity, and low coverage. We further show that clustering RAD‐seq data using the BLASTN and SiLiX programs significantly improves the recovery of orthologous RAD loci compared with previously proposed approaches, especially for distantly related species. This study therefore validates the view that RAD sequencing is a powerful tool for phylogenetic inference.     

Nuclear loci and coalescent methods support ancient hybridization as cause of mitochondrial paraphyly between gadwall and falcated duck (Anas spp.)

Many species have mitochondrial DNA lineages that are phylogenetically intermixed with other species, but studies have rarely tested the cause of such paraphyly. In this study, we tested two hypotheses that could explain mitochondrial paraphyly of Holarctic gadwalls (Anas strepera) with respect to Asian falcated ducks (A. falcata). First, hybridization could have resulted in falcated duck mitochondrial DNA (mtDNA) introgressing into the gadwall gene pool. Second, gadwalls and falcated ducks could have diverged so recently that mtDNA lineages have not sorted to reciprocal monophyly. We used coalescent analyses of three independent loci to distinguish between these two hypotheses. Two lines of evidence support introgression. First, analyses of the three loci combined show that some introgression is necessary to explain current genetic diversity in gadwalls. Second, we generated alternative predictions regarding time since divergence estimated from mtDNA: falcated ducks and gadwalls would have diverged between 65,000 and 700,000 years before present (ybp) under the introgression hypothesis and between 11,000 and 76,000 ybp under the incomplete lineage sorting hypothesis. The two independent nuclear introns indicated that these species diverged between 210,000 and 5,200,000 ybp, which did not overlap the predicted time for incomplete lineage sorting. These analyses also suggested that ancient introgression ( approximately 14,000 ybp) has resulted in the widespread distribution and high frequency of falcated-like mtDNA (5.5% of haplotypes) in North America. This is the first study to use a rigorous quantitative framework to reject incomplete lineage sorting as the cause of mitochondrial paraphyly.     

Concatenation and concordance in the reconstruction of mouse lemur phylogeny: an empirical demonstration of the effect of allele sampling in phylogenetics

The systematics and speciation literature is rich with discussion relating to the potential for gene tree/species tree discordance. Numerous mechanisms have been proposed to generate discordance, including differential selection, long-branch attraction, gene duplication, genetic introgression, and/or incomplete lineage sorting. For speciose clades in which divergence has occurred recently and rapidly, recovering the true species tree can be particularly problematic due to incomplete lineage sorting. Unfortunately, the availability of multilocus or "phylogenomic" data sets does not simply solve the problem, particularly when the data are analyzed with standard concatenation techniques. In our study, we conduct a phylogenetic study for a nearly complete species sample of the dwarf and mouse lemur clade, Cheirogaleidae. Mouse lemurs (genus, Microcebus) have been intensively studied over the past decade for reasons relating to their high level of cryptic species diversity, and although there has been emerging consensus regarding the evolutionary diversity contained within the genus, there is no agreement as to the inter-specific relationships within the group. We attempt to resolve cheirogaleid phylogeny, focusing especially on the mouse lemurs, by employing a large multilocus data set. We compare the results of Bayesian concordance methods with those of standard gene concatenation, finding that though concatenation yields the strongest results as measured by statistical support, these results are found to be highly misleading. By employing an approach where individual alleles are treated as operational taxonomic units, we show that phylogenetic results are substantially influenced by the selection of alleles in the concatenation process.     

Isolation‐driven divergence: speciation in a widespread North American songbird (Aves: Certhiidae)

Lineage, or true ‘species’, trees may differ from gene trees because of stochastic processes in molecular evolution leading to gene‐tree heterogeneity. Problems with inferring species trees because of excessive incomplete lineage sorting may be exacerbated in lineages with rapid diversification or recent divergences necessitating the use of multiple loci and individuals. Many recent multilocus studies that investigate divergence times identify lineage splitting to be more recent than single‐locus studies, forcing the revision of biogeographic scenarios driving divergence. Here, we use 21 nuclear loci from regional populations to re‐evaluate hypotheses identified in an mtDNA phylogeographic study of the Brown Creeper (Certhia americana), as well as identify processes driving divergence. Nuclear phylogeographic analyses identified hierarchical genetic structure, supporting a basal split at approximately 32°N latitude, splitting northern and southern populations, with mixed patterns of genealogical concordance and discordance between data sets within the major lineages. Coalescent‐based analyses identify isolation, with little to no gene flow, as the primary driver of divergence between lineages. Recent isolation appears to have caused genetic bottlenecks in populations in the Sierra Madre Oriental and coastal mountain ranges of California, which may be targets for conservation concerns.     

Reconstructing the phylogeny of "Buarremon" brush-finches and near relatives (Aves, Emberizidae) from individual gene trees

Gene trees are often assumed to be equivalent to species trees, but processes such as incomplete lineage sorting can generate incongruence among gene topologies and analyzing multilocus data in concatenated matrices can be prone to systematic errors. Accordingly, a variety of new methods have been developed to estimate species trees using multilocus data sets. Here, we apply some of these methods to reconstruct the phylogeny of Buarremon and near relatives, a group in which phylogenetic analyses of mitochondrial DNA sequences produced results that were inconsistent with relationships implied by a taxonomy based on variation in external phenotype. Gene genealogies obtained for seven loci (one mitochondrial, six nuclear) were varied, with some supporting and some rejecting the monophyly of Buarremon. Overall, our species-tree analyses tended to support a monophyletic Buarremon, but due to lack of congruence between methodologies, resolution of the phylogeny of this group remains uncertain. More generally, our study indicates that the number of individuals sampled can have an important effect on phylogenetic reconstruction, that the use of seven markers does not guarantee obtaining a strongly-supported species tree, and that methods for species-tree reconstruction can produce different results using the same data; these are important considerations for researchers using these new phylogenetic approaches in other systems.     

Multilocus phylogenetics of a rapid radiation in the genus Thomomys (Rodentia: Geomyidae)

Species complexes undergoing rapid radiation present a challenge in molecular systematics because of the possibility that ancestral polymorphism is retained in component gene trees. Coalescent theory has demonstrated that gene trees often fail to match lineage trees when taxon divergence times are less than the ancestral effective population sizes. Suggestions to increase the number of loci and the number of individuals per taxon have been proposed; however, phylogenetic methods to adequately analyze these data in a coalescent framework are scarce. We compare two approaches to estimating lineage (species) trees using multiple individuals and multiple loci: the commonly used partitioned Bayesian analysis of concatenated sequences and a modification of a newly developed hierarchical Bayesian method (BEST) that simultaneously estimates gene trees and species trees from multilocus data. We test these approaches on a phylogeny of rapidly radiating species wherein divergence times are likely to be smaller than effective population sizes, and incomplete lineage sorting is known, in the rodent genus, Thomomys. We use seven independent noncoding nuclear sequence loci (total approximately 4300 bp) and between 1 and 12 individuals per taxon to construct a phylogenetic hypothesis for eight Thomomys species. The majority-rule consensus tree from the partitioned concatenated analysis included 14 strongly supported bipartitions, corroborating monophyletic species status of five of the eight named species. The BEST tree strongly supported only the split between the two subgenera and showed very low support for any other clade. Comparison of both lineage trees to individual gene trees revealed that the concatenation method appears to ignore conflicting signals among gene trees, whereas the BEST tree considers conflicting signals and downweights support for those nodes. Bayes factor analysis of posterior tree distributions from both analyses strongly favor the model underlying the BEST analysis. This comparison underscores the risks of overreliance on results from concatenation, and ignoring the properties of coalescence, especially in cases of recent, rapid radiations.     

Ecological speciation in the East Maui-endemic Dubautia (Asteraceae) species

Abstract Adaptive diversification can be viewed as a balance between the conservative force of interpopulation gene flow and selection for differential environments. In this paper, we examine ecological, morphological, and genetic differentiation in a small clade consisting of four East Maui-endemic species of Dubautia : D. menziesii , D. platyphylla , D. reticulata , and D. waianapanapaensis , in the Hawaiian silversword alliance (Asteraceae). The East Maui clade is apparently recently derived (less than 1 million years ago) and is geographically restricted yet displays significant ecological and morphological differences. We used geographic data from historical herbarium specimens, measurements of plant architecture and leaf morphometrics, and measures of genetic differentiation in both microsatellite and nuclear coding loci to examine the correlation of different forms of divergence in this small species flock. We found overlap in large-scale geographic distributions, significant differentiation in most habitat factors, significant micro-satellite differentiation, and many shared alleles at nuclear coding loci suggesting on-going lineage sorting. Despite the presence of apparent hybrids in some populations, microsatellite variation is consistent with isolation among species. Using Mantel tests, we compared the direction and extent of diversification among different datasets, to determine whether ecological/morphological divergence was correlated with genetic divergence. Correlations among different datasets showed that habitat was strongly correlated with plant architecture but not leaf morphology. Taken together, these results indicate that ecological and morphological diversification has driven genetic divergence at rapidly evolving microsatellite loci, whereas there is continuing lineage sorting at neutral sites in nuclear coding loci.     

Distinguishing noise from signal in patterns of genomic divergence in a highly polymorphic avian radiation

Recently diverged taxa provide the opportunity to search for the genetic basis of the phenotypes that distinguish them. Genomic scans aim to identify loci that are diverged with respect to an otherwise weakly differentiated genetic background. These loci are candidates for being past targets of selection because they behave differently from the rest of the genome that has either not yet differentiated or that may cross species barriers through introgressive hybridization. Here we use a reduced‐representation genomic approach to explore divergence among six species of southern capuchino seedeaters, a group of recently radiated sympatric passerine birds in the genus Sporophila. For the first time in these taxa, we discovered a small proportion of markers that appeared differentiated among species. However, when assessing the significance of these signatures of divergence, we found that similar patterns can also be recovered from random grouping of individuals representing different species. A detailed demographic inference indicates that genetic differences among Sporophila species could be the consequence of neutral processes, which include a very large ancestral effective population size that accentuates the effects of incomplete lineage sorting. As these neutral phenomena can generate genomic scan patterns that mimic those of markers involved in speciation and phenotypic differentiation, they highlight the need for caution when ascertaining and interpreting differentiated markers between species, especially when large numbers of markers are surveyed. Our study provides new insights into the demography of the southern capuchino radiation and proposes controls to distinguish signal from noise in similar genomic scans.     

Using patterns of genetic structure based on microsatellite loci to test hypotheses of current hybridization,ancient hybridization and incomplete lineage sorting in Conradina (Lamiaceae)

Hybridization and/or incomplete sorting of ancestral polymorphism are commonly implicated to explain discordant phylogenetic analyses of closely related species complexes. One genus in which these phenomena have been suggested to have played major roles based on phylogenetic data is Conradina, a genus of mints (Lamiaceae) endemic to the southeastern USA containing several endangered species. The goals of this study were to use microsatellite data to better understand patterns of genetic structure in Conradina, to test hypotheses of recent or ancient hybridization and incomplete lineage sorting, and to clarify species boundaries. Individuals from 55 populations representing all Conradina species were genotyped using 10 microsatellite loci. Analyses of the patterns of genetic structure in Conradina revealed a clear differentiation of populations following recognized species boundaries, indicating that species have diverged from one another genetically and interspecific hybridization has not occurred recently. Neither ancient hybridization nor incomplete lineage sorting is supported as the sole cause of species nonmonophyly, suggesting that both may have contributed to patterns found in phylogenetic trees; however, analyses of other types of data may be more appropriate to distinguish between these two hypotheses. Because all described species appear to be valid entities, the current listing status of most endangered species of Conradina is appropriate; however, populations of Conradina canescens are genetically differentiated into three groups, each of which may merit species status, and several recently discovered populations of Conradina in Dunn's Creek State Park in Florida are highly differentiated genetically and also appear to represent a new species.