Reconciling gene trees of eastern North American minnows

Most eastern North American cyprinid fishes belong to a clade known as the "open posterior myodome" (OPM) minnows, but phylogenetic relationships within this clade have been difficult to ascertain. Previous attempts to resolve relationships among the generally benthic "chubs" and the more pelagic "shiners", that constitute the majority of OPM minnows, have led to highly discordant phylogenetic hypotheses. To further examine relationships among the OPM minnows, we utilized both a concatenated Bayesian approach and a coalescent-based species tree method to analyze data from six protein coding nuclear loci (Enc1, Ptr, Ryr3, Sh3px3, Tbr1, and Zic1), as well as the mitochondrial locus (Cytb). We focused our analyses on the chub-like genus Phenacobius, a group that has drifted topologically between other benthic chubs and the more pelagic shiners, and also included exemplar taxa from 11 other OPM lineages. Individual gene trees were highly discordant regarding relationships within Phenacobius and across the OPM clade. The concatenated Bayesian analysis and coalescent-based species tree reconstruction recovered slightly different phylogenetic topologies. Additionally, the posterior support values for clades using the coalescent-based approach were consistently lower than the concatenated analysis. However, Phenacobius was resolved as monophyletic and as the sister lineage to Erimystax regardless of the combined data approach taken. Furthermore, Phenacobius+Erimystax was recovered as more closely related to the shiners we examined than to other chubs. Relationships within Phenacobius varied depending on the combined phylogenetic method utilized. Our results highlight the importance of multi-locus, coalescent-based approaches for resolving the phylogeny of diverse clades like the eastern North American OPM minnows.     

Species tree estimation and the historical biogeography of heroine cichlids

Heroine cichlids are major components of the fish faunas in both Central America and the Caribbean. To examine the evolutionary patterns of how cichlids colonized both of these regions, we reconstructed the phylogenetic relationships among 23 cichlid lineages. We used three phylogenetically novel nuclear markers (Dystropin b, Myomesin1, and Wnt7b) in combination with sequence data from seven other gene regions (Nd2, Rag1, Enc1, Sreb2, Ptr, Plagl2, and Zic1) to elucidate the species tree of these cichlids. The species examined represent major heroine lineages in South America, Central America, and the Greater Antilles. The individual gene trees of these groups were topologically quite discordant. Therefore, we combined the genetic partitions and inferred the species tree using both concatenation and a coalescent-based Bayesian method. The two resulting phylogenetic topologies were largely concordant but differed in two fundamental ways. First, more nodes in the concatenated tree were supported with substantial or 100% Bayesian posterior support than in the coalescent-based tree. Second, there was a minor, but biogeographically critical, topological difference between the concatenated and coalescent-based trees. Nevertheless, both analyses recovered topologies consistent with the Greater Antillean heroines being phylogenetically nested within the largely Central American heroine radiation. This study suggests that reconstructions of cichlid phylogeny and historical biogeography should account for the vagaries of individual gene histories.     

Phylogeny of the Tachycineta genus of New World swallows: insights from complete mitochondrial genomes

The Tachycineta genus of swallows is comprised of nine species that range from Alaska to southern Chile. We sequenced the entire mitochondrial genome of each member of Tachycineta and generated a completely resolved phylogenetic hypothesis for the corresponding mitochondrial gene tree. Our analyses confirm the presence of two sub-clades within Tachycineta that are associated with geography: a North American/Caribbean clade and a South/Central American clade. We found considerable variation among regions of the mitochondrial genome in both substitution rates and the level of information that each region supplied for phylogenetic reconstruction. We found no evidence of positive directional selection within mitochondrial coding regions, but we identified numerous sites under purifying selection. This finding suggests that, despite differences in life history traits and distributions, mitochondrial genes in Tachycineta are predominantly under purifying selection for conserved function.     

Lizards from the end of the world: phylogenetic relationships of the Liolaemus lineomaculatus section (Squamata: Iguania: Liolaemini)

The Liolaemus lineomaculatus section is a geographically widely distributed group of lizards from the Patagonian region of southern South America, and includes 18 described species representing the most southerly distributed Liolaemus taxa (the genus includes 228 species and extends from Tierra del Fuego north to south-central Peru). Despite high species diversity, the phylogenetic relationships of this section are unknown. In the present work we sampled all described species in the L. lineomaculatus section as well as currently undescribed candidate species to reconstruct the first complete phylogenetic hypothesis for the clade. Our data set included four anonymous nuclear loci, three nuclear protein-coding loci, and two mitochondrial genes. We compared results obtained with three different phylogenetic methods for the concatenated data set (Maximum Parsimony, Maximum Likelihood and Bayesian Inference) with a coalescent-based species tree approach (BEST), and recovered congruent, strongly-supported topological arrangements across all methods. We identified four main clades within the L. lineomaculatus section: the lineomaculatus, magellanicus, somuncurae, and kingii+archeforus groups, for which we estimated divergence times. We discuss the taxonomic implications of these results and how the future integration of phylogeographic, niche modeling and morphological approaches will allow testing biogeographical hypotheses in this clade.     

Cryptic Genetic Diversity Is Paramount in Small-Bodied Amphibians of the Genus Euparkerella (Anura: Craugastoridae) Endemic to the Brazilian Atlantic Forest

Morphological similarity associated to restricted distributions and low dispersal abilities make the direct developing “Terrarana” frogs of the genus Euparkerella a good model for examining diversification processes. We here infer phylogenetic relationships within the genus Euparkerella, using DNA sequence data from one mitochondrial and four nuclear genes coupled with traditional Bayesian phylogenetic reconstruction approaches and more recent coalescent methods of species tree inference. We also used Bayesian clustering analysis and a recent Bayesian coalescent-based approach specifically to infer species delimitation. The analysis of 39 individuals from the four known Euparkerella species uncovered high levels of genetic diversity, especially within the two previously morphologically-defined E. cochranae and E. brasiliensis. Within these species, the gene trees at five independent loci and trees from combined data (concatenated dataset and the species tree) uncovered six deeply diverged and geographically coherent evolutionary units, which may have diverged between the Miocene and the Pleistocene. These six units were also uncovered in the Bayesian clustering analysis, and supported by the Bayesian coalescent-based species delimitation (BPP), and Genealogical Sorting Index (GSI), providing thus strong evidence for underestimation of the current levels of diversity within Euparkerella. The cryptic diversity now uncovered opens new opportunities to examine the origins and maintenance of microendemism in the context of spatial heterogeneity and/or human induced fragmentation of the highly threatened Brazilian Atlantic forest hotspot.     

Multilocus phylogeny and systematics of Iberian endemic Squalius (Actinopterygii,Leuciscidae)

Inferring the evolutionary history of a group of species can be challenging given the many factors involved. In recent years, the increased availability of sequences of multiple genes per species has spurred the development of new methodologies to analyse multilocus data sets. Two approaches that analyse such data are concatenated supermatrix and coalescent-based species-tree analyses. In this study, we used both of these methods to infer the phylogenetic relationships of Iberian species of the genus Squalius from one mitochondrial and six nuclear genes. We found mitonuclear discordance in the phylogenetic relationships of the group. According to the mitochondrial gene analysis, all species were recovered as monophyletic except S. pyrenaicus; besides, in the concatenated supermatrix analysis of the nuclear markers, this species resolved as polyphyletic with three divergent evolutionary lineages. The coalescent-based nuclear species-tree analysis rendered a well-resolved phylogeny compared with the supermatrix analysis, which was unable to discern between S. carolitertii, S. castellanus and one of the evolutionary lineages of S. pyrenaicus. This result is likely due to the better integration of population uncertainty in the coalescent approach. Furthermore, Bayesian multilocus species delimitation analyses based on a BPP approach strongly supported the distinct nuclear lineages as different species. Nevertheless, the supermatrix analysis was able to obtain well-supported relationships in the divergent lineages with low numbers of individuals. Our study highlights the usefulness of different analytical methodologies to obtain a more complete picture of the evolutionary history of taxa, especially when discordant patterns among genes are found.     

Multi-locus species tree of the chub genus <Emphasis Type="Italic">Squalius</Emphasis> (Leuciscinae: Cyprinidae) from western Iberia: new insights into its evolutionary history

The phylogenetic relationships of the genus Squalius are believed to be well established based on the mitochondrial cytochrome b gene. Here, we inferred the phylogenetic relationships of all species inhabiting most of the western Iberia river systems using a nuclear multi-locus approach and different species tree methods: concatenation and coalescent-based methods (BEST and minimize-deep-coalescence). The dataset comprised sequences of seven coding and three non-coding regions belonging to seven nuclear genes, which were chosen to cover multiple biological functions: amh, bmp4, ef1a, egr2, irbp, rh and rpl8. We provide evidence for a conflicting topology between the nuDNA species tree and the widely reported mtDNA gene tree. S. pyrenaicus is rendered paraphyletic in all nuDNA species trees, with populations of the Tagus/Colares clustering with S. carolitertii, while populations from the Guadiana, Sado and Almargem form a separate clade. Although a larger sampling size encompassing the full spectrum of Squalius populations in western Iberia is still needed to fully elucidate the phylogeography and species delimitation of this genus, our results suggest that the two S. pyrenaicus clades may represent different species.     

Species tree inference in a recent radiation of orioles (Genus Icterus): multiple markers and methods reveal cytonuclear discordance in the northern oriole group

Recent computational advances provide novel opportunities to infer species trees based on multiple independent loci. Thus, single gene trees no longer need suffice as proxies for species phylogenies. Several methods have been developed to deal with the challenges posed by incomplete and stochastic lineage sorting. In this study, we employed four Bayesian methods to infer the phylogeny of a clade of 11 recently diverged oriole species within the genus Icterus. We obtained well-resolved and mostly congruent phylogenies using a set of seven unlinked nuclear intron loci and sampling multiple individuals per species. Most notably, Bayesian concordance analysis generally agreed well with concatenation; the two methods agreed fully on eight of nine nodes. The coalescent-based method ^∗BEAST further supported six of these eight nodes. The fourth method used, BEST, failed to converge despite exhaustive efforts to optimize the tree search. Overall, the results obtained by new species tree methods and concatenation generally corroborate our findings from previous analyses and data sets. However, we found striking disagreement between mitochondrial and nuclear DNA involving relationships within the northern oriole group. Our results highlight the danger of reliance on mtDNA alone for phylogenetic inference. We demonstrate that in spite of low variability and incomplete lineage sorting, multiple nuclear loci can produce largely congruent phylogenies based on multiple species tree methods, even for very closely-related species.     

Molecular systematics of the deep-sea hydrothermal vent endemic Brachyuran family Bythograeidae: a comparison of three Bayesian species tree methods

Brachyuran crabs of the family Bythograeidae are endemic to deep-sea hydrothermal vents and represent one of the most successful groups of macroinvertebrates that have colonized this extreme environment. Occurring worldwide, the family includes six genera (Allograea, Austinograea, Bythograea, Cyanagraea, Gandalfus, and Segonzacia) and fourteen formally described species. To investigate their evolutionary relationships, we conducted Maximum Likelihood and Bayesian molecular phylogenetic analyses, based on DNA sequences from fragments of three mitochondrial genes (16S rDNA, Cytochrome oxidase I, and Cytochrome b) and three nuclear genes (28S rDNA, the sodium-potassium ATPase a-subunit 'NaK', and Histone H3A). We employed traditional concatenated (i.e., supermatrix) phylogenetic methods, as well as three recently developed Bayesian multilocus methods aimed at inferring species trees from potentially discordant gene trees. We found strong support for two main clades within Bythograeidae: one comprising the members of the genus Bythograea; and the other comprising the remaining genera. Relationships within each of these two clades were partially resolved. We compare our results with an earlier hypothesis on the phylogenetic relationships among bythograeid genera based on morphology. We also discuss the biogeography of the family in the light of our results. Our species tree analyses reveal differences in how each of the three methods weighs conflicting phylogenetic signal from different gene partitions and how limits on the number of outgroup taxa may affect the results.     

Incongruence of mitochondrial and nuclear gene trees in the Carabid beetles Ohomopterus

We studied the molecular phylogeny of the carabid subgenus Ohomopterus (genus Carabus), using two mitochondrial (mt) DNA regions (16SrRNA and NADH dehydrogenase subunit 5) and three nuclear DNA regions (wingless, phosphoenolpyruvate carboxykinase, and an anonymous locus). We revisited the previously reported incongruence between the distribution of mtDNA markers and morphologically defined species (Su et al., 1996; J. Mol. Evol. 43:662-671), which those authors attributed to "type switching", a concerted change in many morphological characters that results in the repeated evolution of a particular morphological type. Our mtDNA gene tree obtained from 44 individuals representing all 15 currently recognized species of Ohomopterus revealed that haplotypes isolated from individuals of a single "species" were frequently separated into distant clades, confirming the previous report. The three nuclear markers generally conformed better-with the morphologically defined species than did the mitochondrial markers. The phylogenetic signal in mtDNA and nuclear DNA data differed strongly, and these two partitions were significantly incongruent with each other according to the incongruence length difference test of Farris et al. (1994; Cladistics 10:315-320), although the three nuclear partitions were not homogeneous either. Our results did not support the type-switching hypothesis that had been proposed to fit the morphological data to the mitochondrial gene tree: The incongruence of the mtDNA tree with other nuclear markers indicates that the mtDNA-based tree does not reflect species history any better than the morphological data do. Incongruence of gene trees in Ohomopterus may have been promoted by the complex processes of geographic isolation and hybridization in the Japanese Archipelago that have led to occasional gene flow and recombination between separated entities. The occurrence of reticulate patterns in this group is intriguing, because species of Ohomopterus exhibit extremely divergent genitalic structures that represent a highly efficient reproductive isolation mechanism.     

Mitonuclear discordance is caused by rampant mitochondrial introgression in Neodiprion (Hymenoptera: Diprionidae) sawflies

We investigate the pervasiveness of hybridization and mitochondrial introgression in Neodiprion Rohwer (Hymenoptera; Diprionidae), a Holarctic genus of conifer-feeding sawflies. A phylogenetic analysis of the lecontei species group revealed extensive discordance between a contiguous mitochondrial region spanning three genes (COI, tRNA-leucine, and COII) and three nuclear loci (EF1alpha, CAD, and an anonymous nuclear locus). Bayesian tests of monophyly and Shimodaira-Hasegawa (SH) tests of topological congruence were consistent with mitochondrial introgression; however, these patterns could also be explained by lineage sorting (i.e., deep coalescence). Therefore, to explicitly test the mitochondrial introgression hypothesis, we used a novel application of coalescent-based isolation with migration (IM) models to measure interspecific gene flow at each locus. In support of our hypothesis, mitochondrial gene flow was consistently higher than nuclear gene flow across 120 pairwise species comparisons (P < 1 x 10(-12)). We combine phylogenetic and coalescent evidence to identify likely cases of recent and ancient introgression in Neodiprion, and based on these observations, we hypothesize that shared hosts and/or pheromones facilitate hybridization, whereas disparate abundances between hybridizing species promote mitochondrial introgression. Our results carry implications for phylogenetic analysis, and we advocate the separation of high and low gene flow regions to inform analyses of hybridization and speciational history, respectively.     

Phylogenetic relationships of phrynosomatid lizards based on nuclear and mitochondrial data,and a revised phylogeny for Sceloporus

Phrynosomatid lizards are among the most common and diverse groups of reptiles in western North America, Mexico, and Central America. Phrynosomatidae includes 136 species in 10 genera. Phrynosomatids are used as model systems in many research programs in evolution and ecology, and much of this research has been undertaken in a comparative phylogenetic framework. However, relationships among many phrynosomatid genera are poorly supported and in conflict between recent studies. Further, previous studies based on mitochondrial DNA sequences suggested that the most species-rich genus (Sceloporus) is possibly paraphyletic with respect to as many as four other genera (Petrosaurus, Sator, Urosaurus, and Uta). Here, we collect new sequence data from five nuclear genes and combine them with published data from one additional nuclear gene and five mitochondrial gene regions. We compare trees from nuclear and mitochondrial data from 37 phrynosomatid taxa, including a “species tree” (from BEST) for the nuclear data. We also present a phylogeny for 122 phrynosomatid species based on maximum likelihood analysis of the combined data, which provides a strongly-supported hypothesis for relationships among most phrynosomatid genera and includes most phrynosomatid species. Our results strongly support the monophyly of Sceloporus (including Sator) and many of the relationships within it. We present a new classification for phrynosomatid lizards and the genus Sceloporus, and offer a new tree with branch lengths for use in comparative studies.     

Multilocus analyses of an Antarctic fish species flock (Teleostei, Notothenioidei, Trematominae): phylogenetic approach and test of the early-radiation event

Clades that have undergone episodes of rapid cladogenesis are challenging from a phylogenetic point of view. They are generally characterised by short or missing internal branches in phylogenetic trees and by conflicting topologies among individual gene trees. This may be the case of the subfamily Trematominae, a group of marine teleosts of coastal Antarctic waters, which is considered to have passed through a period of rapid diversification. Despite much phylogenetic attention, the relationships among Trematominae species remain unclear. In contrast to previous studies that were mostly based on concatenated datasets of mitochondrial and/or single nuclear loci, we applied various single-locus and multilocus phylogenetic approaches to sequences from 11 loci (eight nuclear) and we also used several methods to assess the hypothesis of a radiation event in Trematominae evolution. Diversification rate analyses support the hypothesis of a period of rapid diversification during Trematominae history and only a few nodes in the hypothetical species tree were consistently resolved with various phylogenetic methods. We detected significant discrepancies among trees from individual genes of these species, most probably resulting from incomplete lineage sorting, suggesting that concatenation of loci is not the most appropriate way to investigate Trematominae species interrelationships. These data also provide information about the possible effects of historic climate changes on the diversification rate of this group of fish.     

Molecular phylogenetic relationships and phenotypic diversity in miniaturized toadlets, genus Brachycephalus (Amphibia: Anura: Brachycephalidae)

Toadlets of the genus Brachycephalus are endemic to the Atlantic rainforests of southeastern and southern Brazil. The 14 species currently described have snout-vent lengths less than 18 mm and are thought to have evolved through miniaturization: an evolutionary process leading to an extremely small adult body size. Here, we present the first comprehensive phylogenetic analysis for Brachycephalus, using a multilocus approach based on two nuclear (Rag-1 and Tyr) and three mitochondrial (Cyt b, 12S, and 16S rRNA) gene regions. Phylogenetic relationships were inferred using a partitioned Bayesian analysis of concatenated sequences and the hierarchical Bayesian method (BEST) that estimates species trees based on the multispecies coalescent model. Individual gene trees showed conflict and also varied in resolution. With the exception of the mitochondrial gene tree, no gene tree was completely resolved. The concatenated gene tree was completely resolved and is identical in topology and degree of statistical support to the individual mtDNA gene tree. On the other hand, the BEST species tree showed reduced significant node support relative to the concatenate tree and recovered a basal trichotomy, although some bipartitions were significantly supported at the tips of the species tree. Comparison of the log likelihoods for the concatenated and BEST trees suggests that the method implemented in BEST explains the multilocus data for Brachycephalus better than the Bayesian analysis of concatenated data. Landmark-based geometric morphometrics revealed marked variation in cranial shape between the species of Brachycephalus. In addition, a statistically significant association was demonstrated between variation in cranial shape and genetic distances estimated from the mtDNA and nuclear loci. Notably, B. ephippium and B. garbeana that are predicted to be sister-species in the individual and concatenated gene trees and the BEST species tree share an evolutionary novelty, the hyperossified dorsal plate.     

Mosaic genealogy of the Mus musculus genome revealed by 21 nuclear genes from its three subspecies

Patterns of genetic variation provide insight into the evolutionary history of a species. Mouse (Mus musculus) is a good model for this purpose. Here we present the analysis of genealogies of the 21 nuclear loci and one mitochondrial DNA region in M. musculus based on our nucleotide sequences of nine inbred strains from three M. musculus subspecies (musculus, domesticus, and castaneus) and one M. spicilegus strain as an outgroup. The mitochondrial DNA gene genealogy of those strains confirmed the introgression pattern of one musculus strain. When all the nuclear DNA data were concatenated to produce a phylogenetic tree of nine strains, musculus and domesticus strains formed monophyletic clusters with each other, while the two castaneus strains were paraphyletic. When each DNA region was treated independently, the phylogenetic networks revealed an unnegligibly high level of subspecies admixture and the mosaic nature of their genome. Estimation of ancestral and derived population sizes and migration rates suggests the effects of ancestral polymorphism and gene flow on the pattern of genetic variation of the current subspecies. Gene genealogies of Fut4 and Dfy loci also suggested existence of the gene flow between M. musculus and M. spicilegus or other distant species.     

Parallel evolution in eight-barbel loaches of the genus Lefua (Balitoridae, Cypriniformes) revealed by mitochondrial and nuclear DNA phylogenies

The evolutionary history of eight-barbel loaches of the genus Lefua contains important phylogenetic information that will aid in resolution of the faunal formations and evolutionary histories of Japanese and East Asian freshwater fishes. Our sequencing of the mitochondrial D-loop region in a large number of samples allowed construction of the most comprehensive phylogeny of these loaches to date; we demonstrated monophyly of five Lefua species and identified populations of Lufua. sp. and Lefua echigonia. Loaches inhabiting the Tokai region in Japan were morphologically and ecologically indistinguishable from Lefua sp. However, they were included in the L. echigonia lineage. We determined a novel phylogeny by sequencing the nuclear ribosomal S7 subunit and showed that nuclear DNA phylogeny essentially matched the mitochondrial DNA phylogeny. Loaches from the Tokai region were part of the L. echigonia lineage, indicating parallel evolution between Tokai loaches and Lefua sp. in western Japan. We presented the most robust phylogeny to date using concatenated mitochondrial and nuclear sequences. The wealth of molecular information allowed us to speculate on evolutionary processes in the genus Lefua.     

Gene tree parsimony of multilocus snake venom protein families reveals species tree conflict as a result of multiple parallel gene loss

The proliferation of gene data from multiple loci of large multigene families has been greatly facilitated by considerable recent advances in sequence generation. The evolution of such gene families, which often undergo complex histories and different rates of change, combined with increases in sequence data, pose complex problems for traditional phylogenetic analyses, and in particular, those that aim to successfully recover species relationships from gene trees. Here, we implement gene tree parsimony analyses on multicopy gene family data sets of snake venom proteins for two separate groups of taxa, incorporating Bayesian posterior distributions as a rigorous strategy to account for the uncertainty present in gene trees. Gene tree parsimony largely failed to infer species trees congruent with each other or with species phylogenies derived from mitochondrial and single-copy nuclear sequences. Analysis of four toxin gene families from a large expressed sequence tag data set from the viper genus Echis failed to produce a consistent topology, and reanalysis of a previously published gene tree parsimony data set, from the family Elapidae, suggested that species tree topologies were predominantly unsupported. We suggest that gene tree parsimony failure in the family Elapidae is likely the result of unequal and/or incomplete sampling of paralogous genes and demonstrate that multiple parallel gene losses are likely responsible for the significant species tree conflict observed in the genus Echis. These results highlight the potential for gene tree parsimony analyses to be undermined by rapidly evolving multilocus gene families under strong natural selection.     

Phylogeny estimation of the radiation of western North American chipmunks (Tamias) in the face of introgression using reproductive protein genes

The causes and consequences of rapid radiations are major unresolved issues in evolutionary biology. This is in part because phylogeny estimation is confounded by processes such as stochastic lineage sorting and hybridization. Because these processes are expected to be heterogeneous across the genome, comparison among marker classes may provide a means of disentangling these elements. Here we use introns from nuclear-encoded reproductive protein genes expected to be resistant to introgression to estimate the phylogeny of the western chipmunks (Tamias: subgenus: Neotamias), a rapid radiation that has experienced introgressive hybridization of mitochondrial DNA (mtDNA). We analyze the nuclear loci using coalescent-based species-tree estimation methods and concatenation to estimate a species tree and we use parametric bootstraps and coalescent simulations to differentiate between phylogenetic error, coalescent stochasticity and introgressive hybridization. Results indicate that the mtDNA gene tree reflects several introgression events that have occurred between taxa of varying levels of divergence and at different time points in the tree. T. panamintinus and T. speciosus appear to be fixed for ancient mitochondrial introgressions from T. minimus. A southern Rocky Mountains clade appears well sorted (i.e., species are largely monophyletic) at multiple nuclear loci, while five of six taxa are nonmonophyletic based on cytochrome b. Our simulations reject phylogenetic error and coalescent stochasticity as causes. The results represent an advance in our understanding of the processes at work during the radiation of Tamias and suggest that sampling reproductive-protein genes may be a viable strategy for phylogeny estimation of rapid radiations in which reproductive isolation is incomplete. However, a genome-scale survey that can statistically compare heterogeneity of genealogical process at many more loci will be necessary to test this conclusion.     

Multilocus sequence analysis of the actinobacterial genus Kribbella

Multilocus sequence analysis (MLSA) was used to refine the phylogenetic analysis of the genus Kribbella, which currently contains 17 species with validly-published names. Sequences were obtained for the 16S rRNA, gyrB, rpoB, recA, relA and atpD genes for 16 of the 17 type strains of the genus plus seven non-type strains. A five-gene concatenated sequence of 4099 nt was used to examine the phylogenetic relationships between the species of the genus Kribbella. Using the concatenated sequence of the gyrB-rpoB-recA-relA and atpD genes, most Kribbella type strains can be distinguished by a genetic distance of >0.04. Each single-gene tree had an overall topology similar to that of the concatenated sequence tree. The single-gene relA tree, used here for the first time in MLSA of actinobacteria, had good bootstrap support, comparable to the rpoB and atpD gene trees, which had topologies closest to that of the concatenated sequence tree. This illustrates that relA is a useful addition in MLSA studies of the genus Kribbella. We propose that concatenated gyrB-rpoB-recA-relA-atpD gene sequences be used for examining the phylogenetic relationships within the genus Kribbella and for determining the closest phylogenetic relatives to be used for taxonomic comparisons.