Molecular differentiation and phylogeny of entomopathogenic nematodes (rhabditida: heterorhabditidae) based on ND4 gene sequences of Mitochondrial DNA.

We determined partial ND4 gene sequences of mitochondrial DNA from 15 heterorhabditid nematode isolates, representing 5 species collected from different regions of the world, by using polymerase chain reaction (PCR) and direct-sequencing of PCR products. Aligned nucleotide as well as amino acid sequences were used to differentiate nematode species by comparing sequence divergence and to infer phylogeny of the nematodes by using maximum parsimony and likelihood methods. Robustness of our phylogenetic trees was checked by bootstrap tests. The 15 nematode isolates can be divided into 7 haplotypes based on DNA sequences. On a larger scale, the sequence divergence revealed 4 distinct groups corresponding to 4 described species. No sequence divergence was detected from 5 isolates of Heterorhabditis bacteriophora or between Heterorhabditis marelatus to Heterorhabditis hepialius. Our sequence data yielded phylogenetic trees with identical topologies when different tree-building methods were used. Most relationships were also confirmed by using amino acid sequences in maximum parsimony analysis. Our molecular phylogeny of Heterorhabditis species support an existing taxonomy that is based largely on morphology and the sequence divergence of the ND4 gene permits species identification.     

The relationship of nitrate reducing bacteria on the basis of narH gene sequences and comparison of narH and 16S rDNA based phylogeny

On the basis of available nitrate reductase gene sequences primer pairs were designed to specifically amplify gene stretches of the beta-subunit of the membrane-bound nitrate reductase (narH). Additional sequences of this gene were amplified and sequenced from pure cultures of reference species and new isolates. The distribution and phylogeny of this gene in denitrifying and nitrate-reducing bacteria was analysed. Comparison of phylogenetic trees based on 16S rDNA sequences with those based on narH sequences revealed highly similar relationships of both genes from most of the bacteria analysed. Since highly conserved functional cysteine clusters within bacterial and archaeal narH sequences support a linear evolution from one common progenitor a long evolutionary history of the respiratory membrane-bound nitrate reductase can be inferred from our phylogenetic data.     

Free from mitochondrial DNA: Nuclear genes and the inference of species trees among closely related darter lineages (Teleostei: Percidae: Etheostomatinae)

Investigations into the phylogenetics of closely related animal species are dominated by the use of mitochondrial DNA (mtDNA) sequence data. However, the near-ubiquitous use of mtDNA to infer phylogeny among closely related animal lineages is tempered by an increasing number of studies that document high rates of transfer of mtDNA genomes among closely related species through hybridization, leading to substantial discordance between phylogenies inferred from mtDNA and nuclear gene sequences. In addition, the recent development of methods that simultaneously infer a species phylogeny and estimate divergence times, while accounting for incongruence among individual gene trees, has ushered in a new era in the investigation of phylogeny among closely related species. In this study we assess if DNA sequence data sampled from a modest number of nuclear genes can resolve relationships of a species-rich clade of North American freshwater teleost fishes, the darters. We articulate and expand on a recently introduced method to infer a time-calibrated multi-species coalescent phylogeny using the computer program ^*BEAST. Our analyses result in well-resolved and strongly supported time-calibrated darter species tree. Contrary to the expectation that mtDNA will provide greater phylogenetic resolution than nuclear gene data; the darter species tree inferred exclusively from nuclear genes exhibits a higher frequency of strongly supported nodes than the mtDNA time-calibrated gene tree.     

Molecular and evolutionary relationships among enteric bacteria

Classification of bacterial species into genera has traditionally relied upon variation in phenotypic characteristics. However, these phenotypes often have a multifactorial genetic basis, making unambiguous taxonomic placement of new species difficult. By designing evolutionarily conserved oligonucleotide primers, it is possible to amplify homologous regions of genes in diverse taxa using the polymerase chain reaction and determine their nucleotide sequences. We have constructed a phylogeny of some enteric bacteria, including five species classified as members of the genus Escherichia, based on nucleotide sequence variation at the loci encoding glyceraldehyde-3-phosphate dehydrogenase and outer membrane protein 3A, and compared this genealogy with the relationships inferred by biotyping. The DNA sequences of these genes defined congruent and robust phylogenetic trees indicating that they are an accurate reflection of the evolutionary history of the bacterial species. The five species of Escherichia were found to be distantly related and, contrary to their placement in the same genus, do not form a monophyletic group. These data provide a framework which allows the relationships of additional species of enteric bacteria to be inferred. These procedures have general applicability for analysis of the classification, evolution, and epidemiology of bacterial taxa.     

Inferring angiosperm phylogeny from EST data with widespread gene duplication

Background

Most studies inferring species phylogenies use sequences from single copy genes or sets of orthologs culled from gene families. For taxa such as plants, with very high levels of gene duplication in their nuclear genomes, this has limited the exploitation of nuclear sequences for phylogenetic studies, such as those available in large EST libraries. One rarely used method of inference, gene tree parsimony, can infer species trees from gene families undergoing duplication and loss, but its performance has not been evaluated at a phylogenomic scale for EST data in plants.

Results

A gene tree parsimony analysis based on EST data was undertaken for six angiosperm model species and Pinus, an outgroup. Although a large fraction of the tentative consensus sequences obtained from the TIGR database of ESTs was assembled into homologous clusters too small to be phylogenetically informative, some 557 clusters contained promising levels of information. Based on maximum likelihood estimates of the gene trees obtained from these clusters, gene tree parsimony correctly inferred the accepted species tree with strong statistical support. A slight variant of this species tree was obtained when maximum parsimony was used to infer the individual gene trees instead.

Conclusion

Despite the complexity of the EST data and the relatively small fraction eventually used in inferring a species tree, the gene tree parsimony method performed well in the face of very high apparent rates of duplication.

     

Phylogenetic diversity of culturable bacteria from Antarctic sandy intertidal sediments

The diversity of culturable bacteria associated with sandy intertidal sediments from the coastal regions of the Chinese Antarctic Zhongshan Station on the Larsemann Hills (Princess Elizabeth Land, East Antarctica) was investigated. A total of 65 aerobic heterotrophic bacterial strains were isolated at 4°C. Microscopy and 16S rRNA gene sequence analysis indicated that the isolates were dominated by Gram-negative bacteria, while only 16 Gram-positive strains were isolated. The bacterial isolates fell in five phylogenetic groups: Alpha- and Gammaproteobacteria, Bacteroidetes, Actinobacteria and Firmicutes. Based on phylogenetic trees, all the 65 isolates were sorted into 29 main clusters, corresponding to at least 29 different genera. Based on sequence analysis (<98% sequence similarity), the Antarctic isolates belonged to at least 37 different bacterial species, and 14 of the 37 bacterial species (37.8%) represented potentially novel taxa. These results indicated a high culturable diversity within the bacterial community of the Antarctic sandy intertidal sediments.     

Diversity and relationships of bradyrhizobia from legumes native to eastern North America

DNA sequencing and polymerase chain reaction (PCR) assays with lineage-specific primers were used to analyze the diversity of 276 isolates of Bradyrhizobium sp. nodule bacteria associated with 13 native legumes species in the northeastern United States, representing eight genera in six legume tribes. A PCR screen with two primer pairs in the rRNA region indicated that seven of the legume species were exclusively associated with strains having markers resembling Bradyrhizobium elkanii, while the remaining six host species harbored strains related to both B. elkanii and Bradyrhizobium japonicum. Sequence analysis of 22 isolates for portions of 16S rRNA and 23S rRNA yielded congruent phylogenetic trees and showed that isolates from different legume genera often shared similar or identical sequences. However, trees inferred from portions of two other genes (alpha-ketoglutarate dioxygenase gene (tfdA), the alpha-subunit of nitrogenase (nifD)) differed significantly from the rRNA phylogeny. Thus, for Bradyrhizobium populations in this region, lateral gene transfer events appear to have altered genealogical relationships of different portions of the genome. These results extend the number of likely cases of gene transfer between divergent taxa of Bradyrhizobium (from members of the B. elkanii lineage to the B. japonicum group) and suggest that transfers have also occurred among separate subgroups of the B. elkanii lineage.     

Evolution of aromatic amino acid biosynthesis and application to the fine-tuned phylogenetic positioning of enteric bacteria.

A comprehensive phylogenetic tree for virtually the entire assemblage of enteric bacteria is presented. Character states of aromatic amino acid biosynthesis are used as criteria, and the results are compared with partial trees based upon sequencing of 16S rRNA, 5S rRNA, and tryptophan leader peptide. Three major clusters are apparent. Enterocluster 1 possesses a gene fusion (trpG-trpD) encoding anthranilate synthase: anthranilate 5-phosphoribosylpyrophosphate phosphoribosyltransferase of tryptophan biosynthesis. This cluster includes the genera Escherichia, Shigella, Citrobacter, Salmonella, Klebsiella, and Enterobacter. The remaining two clusters lack the trpG-trpD gene fusion, but differ in the presence (enterocluster 2) or absence (enterocluster 3) of the three-step overflow pathway to L-phenylalanine. Enterocluster 2 consists of the genera Serratia and Erwinia. Enterocluster 3 includes the genera Cedecea, Kluyvera, Edwardsiella, Hafnia, Yersinia, Proteus, Providencia, and Morganella. Within these three major clusters, a tentative hierarchy of subcluster ordering is formulated on the basis of all data available. This hierarchical framework is proposed as a general working basis for continued refinement of the phylogenetic relationships of enteric bacteria.     

Multigene Phylogenetics Reveals Temporal Diversification of Major African Malaria Vectors

The major vectors of malaria in sub-Saharan Africa belong to subgenus Cellia. Yet, phylogenetic relationships and temporal diversification among African mosquito species have not been unambiguously determined. Knowledge about vector evolutionary history is crucial for correct interpretation of genetic changes identified through comparative genomics analyses. In this study, we estimated a molecular phylogeny using 49 gene sequences for the African malaria vectors An. gambiae, An. funestus, An. nili, the Asian malaria mosquito An. stephensi, and the outgroup species Culex quinquefasciatus and Aedes aegypti. To infer the phylogeny, we identified orthologous sequences uniformly distributed approximately every 5 Mb in the five chromosomal arms. The sequences were aligned and the phylogenetic trees were inferred using maximum likelihood and neighbor-joining methods. Bayesian molecular dating using a relaxed log normal model was used to infer divergence times. Trees from individual genes agreed with each other, placing An. nili as a basal clade that diversified from the studied malaria mosquito species 47.6 million years ago (mya). Other African malaria vectors originated more recently, and independently acquired traits related to vectorial capacity. The lineage leading to An. gambiae diverged 30.4 mya, while the African vector An. funestus and the Asian vector An. stephensi were the most closely related sister taxa that split 20.8 mya. These results were supported by consistently high bootstrap values in concatenated phylogenetic trees generated individually for each chromosomal arm. Genome-wide multigene phylogenetic analysis is a useful approach for discerning historic relationships among malaria vectors, providing a framework for the correct interpretation of genomic changes across species, and comprehending the evolutionary origins of this ubiquitous and deadly insect-borne disease.     

rpoB sequence analysis as a novel basis for bacterial identification

Comparison of the sequences of conserved genes, most commonly those encoding 16S rRNA, is used for bacterial genotypic identification. Among some taxa, such as the Enterobacteriaceae, variation within this gene does not allow confident species identification. We investigated the usefulness of RNA polymerase beta-subunit encoding gene ( rpoB  ) sequences as an alternative tool for universal bacterial genotypic identification. We generated a database of partial rpoB for 14 Enterobacteriaceae species and then assessed the intra- and interspecies divergence between the rpoB and the 16S rRNA genes by pairwise comparisons. We found that levels of divergence between the rpoB sequences of different strains were markedly higher than those between their 16S rRNA genes. This higher discriminatory power was further confirmed by assigning 20 blindly selected clinical isolates to the correct enteric species on the basis of rpoB sequence comparison. Comparison of rpoB sequences from Enterobacteriaceae was also used as the basis for their phylogenetic analysis and demonstrated the genus Klebsiella to be polyphyletic. The trees obtained with rpoB were more compatible with the currently accepted classification of Enterobacteriaceae than those obtained with 16S rRNA. These data indicate that rpoB is a powerful identification tool, which may be useful for universal bacterial identification.     

Numerical taxonomy and ecology of petroleum-degrading bacteria.

A total of 99 strains of petroleum-degrading bacteria isolated from Chesapeake Bay water and sediment were identified by using numerical taxonomy procedures. The isolates, together with 33 reference cultures, were examined for 48 biochemical, cultural, morphological, and physiological characters. The data were analyzed by computer, using both the simple matching and the Jaccard coefficients. Clustering was achieved by the unweighted average linkage method. From the sorted similarity matrix and dendrogram, 14 phenetic groups, comprising 85 of the petroleum-degrading bacteria, were defined at the 80 to 85% similarity level. These groups were identified as actinomycetes (mycelial forms, four clusters), coryneforms, Enterobacteriaceae, Klebsiella aerogenes, Micrococcus spp. (two clusters), Nocardia species (two clusters), Pseudomonas spp. (two clusters), and Sphaerotilus natans. It is concluded that the degradation of petroleum is accomplished by a diverse range of bacterial taxa, some of which were isolated only at given sampling stations and, more specifically, from sediment collected at a given station.     

Brassicaceae phylogeny inferred from phytochrome A and ndhF sequence data: tribes and trichomes revisited

The family Brassicaceae comprises 3710 species in 338 genera, 25 recently delimited tribes, and three major lineages based on phylogenetic results from the chloroplast gene ndhF. To assess the credibility of the lineages and newly delimited tribes, we sequenced an approximately 1.8-kb region of the nuclear phytochrome A (PHYA) gene for taxa previously sampled for the chloroplast gene ndhF. Using parsimony, likelihood, and Bayesian methods, we reconstructed the phylogeny of the gene and used the approximately unbiased (AU) test to compare phylogenetic results from PHYA with findings from ndhF. We also combined ndhF and PHYA data and used a Bayesian mixed model approach to infer phylogeny. PHYA and combined analyses recovered the same three large lineages as those recovered in ndhF trees, increasing confidence in these lineages. The combined tree confirms the monophyly of most of the recently delimited tribes (only Alysseae, Anchonieae, and Descurainieae are not monophyletic), while 13 of the 23 sampled tribes are monophyletic in PHYA trees. In addition to phylogenetic results, we documented the trichome branching morphology of species across the phylogeny and explored the evolution of different trichome morphologies using the AU test. Our results indicate that dendritic, medifixed, and stellate trichomes likely evolved independently several times in the Brassicaceae.     

Sorting through the chaff, nDNA gene trees for phylogenetic inference and hybrid identification of annual sunflowers (Helianthus sect. Helianthus)

The annual sunflowers (Helianthus sect. Helianthus) present a formidable challenge for phylogenetic inference because of ancient hybrid speciation, recent introgression, and suspected issues with deep coalescence. Here we analyze sequence data from 11 nuclear DNA (nDNA) genes for multiple genotypes of species within the section to (1) reconstruct the phylogeny of this group, (2) explore the utility of nDNA gene trees for detecting hybrid speciation and introgression; and (3) test an empirical method of hybrid identification based on the phylogenetic congruence of nDNA gene trees from tightly linked genes. We uncovered considerable topological heterogeneity among gene trees with or without three previously identified hybrid species included in the analyses, as well as a general lack of reciprocal monophyly of species. Nonetheless, partitioned Bayesian analyses provided strong support for the reciprocal monophyly of all species except H. annuus (0.89 PP), the most widespread and abundant annual sunflower. Previous hypotheses of relationships among taxa were generally strongly supported (1.0 PP), except among taxa typically associated with H. annuus, apparently due to the paraphyly of the latter in all gene trees. While the individual nDNA gene trees provided a useful means for detecting recent hybridization, identification of ancient hybridization was problematic for all ancient hybrid species, even when linkage was considered. We discuss biological factors that affect the efficacy of phylogenetic methods for hybrid identification.     

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.     

Phylogenetic analysis of gamma-proteobacteria inferred from nucleotide sequence comparisons of the house-keeping genes adk, aroE and gdh: comparisons with phylogeny inferred from 16S rRNA gene sequences

Nucleotide sequence comparisons of three house-keeping genes, adenylate kinase (adk), shikimate dehydrogenase (aroE), and glucose-6-phosphate dehydrogenase (gdh), were used to infer the phylogeny of 33 gamma-proteobacteria. Phylogenetic trees inferred from each gene, and from the concatenated sequences of all three genes, are, in general, similar to a 16S rRNA gene-inferred tree. Similar grouping of bacteria are revealed at the family, genus, species and strain levels in all five trees. The house-keeping genes, however, show a higher rate of nucleotide sequence substitutions. Consequently, they can possibly probe deeper branches of a phylogenetic tree than the 16S rRNA gene. However, because their nucleotide sequences are not as highly conserved among gamma-proteobacteria, family- or genus-specific primers would need to be designed for the amplification of any of these three house-keeping genes. Since these genes are used in multilocus sequence typing, it is expected that the number of sequences publicly available for many taxa will increase over time proving them very useful either at complementing 16S rRNA-inferred phylogenies or for specific, targeted, phylogenetic analysis.     

A comparative analysis of Y chromosome and mtDNA phylogenies of the Hylobates gibbons

ABSTRACT: BACKGROUND: The evolutionary relationships of closely related species have long been of interest to biologists since these species experienced different evolutionary processes in a relatively short period of time. Comparison of phylogenies inferred from DNA sequences with differing inheritance patterns, such as mitochondrial, autosomal, and X and Y chromosomal loci, can provide more comprehensive inferences of the evolutionary histories of species. Gibbons, especially the genus Hylobates, are particularly intriguing as they consist of multiple closely related species which emerged rapidly and live in close geographic proximity. Our current understanding of relationships among Hylobates species is largely based on data from the maternally-inherited mitochondrial DNAs (mtDNAs). RESULTS: To infer the paternal histories of gibbon taxa, we sequenced multiple Y chromosomal loci from 26 gibbons representing 10 species. As expected, we find levels of sequence variation some five times lower than observed for the mitochondrial genome (mtgenome). Although our Y chromosome phylogenetic tree shows relatively low resolution compared to the mtgenome tree, our results are consistent with the monophyly of gibbon genera suggested by the mtgenome tree. In a comparison of the molecular dating of divergences and on the branching patterns of phylogeny trees between mtgenome and Y chromosome data, we found: 1) the inferred divergence estimates were more recent for the Y chromosome than for the mtgenome, 2) the species H. lar and H. pileatus are reciprocally monophyletic in the mtgenome phylogeny but a H. pileatus individual falls into the H. lar Y chromosome clade. CONCLUSIONS: Based on the ~6.4 kb of Y chromosomal DNA sequence data generated for each of the 26 individuals in this study, we provide molecular inferences on gibbon and particularly on Hylobates evolution complementary to those from mtDNA data. Overall, our results illustrate the utility of comparative studies of loci with different inheritance patterns for investigating potential sex specific processes on the evolutionary histories of closely related taxa, and emphasize the need for further sampling of gibbons of known provenance.     

Trees of trees: an approach to comparing multiple alternative phylogenies

Phylogenetic analysis very commonly produces several alternative trees for a given fixed set of taxa. For example, different sets of orthologous genes may be analyzed, or the analysis may sample from a distribution of probable trees. This article describes an approach to comparing and visualizing multiple alternative phylogenies via the idea of a "tree of trees" or "meta-tree." A meta-tree clusters phylogenies with similar topologies together in the same way that a phylogeny clusters species with similar DNA sequences. Leaf nodes on a meta-tree correspond to the original set of phylogenies given by some analysis, whereas interior nodes correspond to certain consensus topologies. The construction of meta-trees is motivated by analogy with construction of a most parsimonious tree for DNA data, but instead of using DNA letters, in a meta-tree the characters are partitions or splits of the set of taxa. An efficient algorithm for meta-tree construction is described that makes use of a known relationship between the majority consensus and parsimony in terms of gain and loss of splits. To illustrate these ideas meta-trees are constructed for two datasets: a set of gene trees for species of yeast and trees from a bootstrap analysis of a set of gene trees in ray-finned fish. A software tool for constructing meta-trees and comparing alternative phylogenies is available online, and the source code can be obtained from the author.     

Gene trees,species and species trees in the <Emphasis Type="Italic">Ctenosaura palearis</Emphasis> clade

The growing use of molecular systematics in conservation has increased the importance of accurate resolution of taxonomic units and relationships. DNA data relate most directly to genealogies, which need not have perfect relationships with species limits and phylogenies. We used a multilocus gene tree approach to elucidate the relationships between four endangered Central American iguanas. We found support for the proposition that the described species taxa correspond to distinct evolutionary lineages warranting individual protection. We combined gene trees to estimate a phylogeny using Bayesian Estimation of Species Trees (BEST), minimizing deep coalescence, Species Trees from Average Ranks (STAR), and traditional concatenation. The estimate from concatenation conflicted with the other methods, likely owing to the disproportionate effect of mtDNA on concatenated analyses. This illustrates the importance of appropriate treatment of multilocus sequence data in phylogenetics. Our results indicate that these species have gone through recent and rapid speciation, resulting in four closely related narrow-range endemics.     

Reconciling gene and genome duplication events: using multiple nuclear gene families to infer the phylogeny of the aquatic plant family Pontederiaceae

Most plant phylogenetic inference has used DNA sequence data from the plastid genome. This genome represents a single genealogical sample with no recombination among genes, potentially limiting the resolution of evolutionary relationships in some contexts. In contrast, nuclear DNA is inherently more difficult to employ for phylogeny reconstruction because major mutational events in the genome, including polyploidization, gene duplication, and gene extinction can result in homologous gene copies that are difficult to identify as orthologs or paralogs. Gene tree parsimony (GTP) can be used to infer the rooted species tree by fitting gene genealogies to species trees while simultaneously minimizing the estimated number of duplications needed to reconcile conflicts among them. Here, we use GTP for five nuclear gene families and a previously published plastid data set to reconstruct the phylogenetic backbone of the aquatic plant family Pontederiaceae. Plastid-based phylogenetic studies strongly supported extensive paraphyly of Eichhornia (one of the four major genera) but also depicted considerable ambiguity concerning the true root placement for the family. Our results indicate that species trees inferred from the nuclear genes (alone and in combination with the plastid data) are highly congruent with gene trees inferred from plastid data alone. Consideration of optimal and suboptimal gene tree reconciliations place the root of the family at (or near) a branch leading to the rare and locally restricted E. meyeri. We also explore methods to incorporate uncertainty in individual gene trees during reconciliation by considering their individual bootstrap profiles and relate inferred excesses of gene duplication events on individual branches to whole-genome duplication events inferred for the same branches. Our study improves understanding of the phylogenetic history of Pontederiaceae and also demonstrates the utility of GTP for phylogenetic analysis.     

Different filtering strategies of genotyping‐by‐sequencing data provide complementary resolutions of species boundaries and relationships in a clade of sexually deceptive orchids

Ongoing hybridization and retained ancestral polymorphism in rapidly radiating lineages could mask recent cladogenetic events. This presents a challenge for the application of molecular phylogenetic methods to resolve differences between closely related taxa. We reanalyzed published genotyping‐by‐sequencing (GBS) data to infer the phylogeny of four species within the Ophrys sphegodes complex, a recently radiated clade of orchids. We used different data filtering approaches to detect different signals contained in the dataset generated by GBS and estimated their effects on maximum likelihood trees, global F_ST and bootstrap support values. We obtained a maximum likelihood tree with high bootstrap support, separating the species by using a large dataset based on loci shared by at least 30% of accessions. Bootstrap and F_ST values progressively decreased when filtering for loci shared by a higher number of accessions. However, when filtering more stringently to retain homozygous and organellar loci, we identified two main clades. These clades group individuals independently from their a priori species assignment, but were associated with two organellar haplotype clusters. We infer that a less stringent filtering preferentially selects for rapidly evolving lineage‐specific loci, which might better delimit lineages. In contrast, when using homozygous/organellar DNA loci the signature of a putative hybridization event in the lineage prevails over the most recent phylogenetic signal. These results show that using differing filtering strategies on GBS data could dissect the organellar and nuclear DNA phylogenetic signal and yield novel insights into relationships between closely related species.