Concatenated Gene Tree Research Articles

Divergence of the freshwater sleeper, Neodontobutishainanensis (Chen, 1985) (Teleostei, Odontobutidae), in the Pearl River basin and on Hainan Island of southern China.

Study of divergence of freshwater fish populations between island and adjacent mainland areas can shed light on the phylogeographical relationships of these regions. Neodontobutishainanensis is a freshwater fish species restricted to Hainan Island and in Guangdong and Guangxi provinces in the southern mainland China. We examine the phylogenetic relationship and population structure of N.hainanensis based on 3,176 nuclear loci using a gene-capture method. STRUCTURE analysis and principal coordinate analyses (PCA) indicate that populations from Guangdong, Guangxi and Hainan are each distinct, except that some individuals of the Guangdong population share minor genetic components with individuals of the Guangxi population. In the concatenated gene tree, the Hainan population is grouped with the Guangdong population, but the coalescent tree groups the Hainan population as the sister to the Guangxi population. Finally, coalescent simulations confirmed the divergence pattern supported by the coalescent tree and revealed a one-way introgression from the Guangxi population to the Guangdong population, which can explain the discordant results supported by the concatenated and coalescent phylogenetic analyses. Due to recent decline of N.hainanensis populations and the genetic patterns in this species, as revealed in this study, the populations in the three areas should be treated as separate conservation units.

Molecular systematics of Pampus (Perciformes: Stromateidae) based on thousands of nuclear loci using target-gene enrichment

Pomfret fishes of the genus Pampus are commercially important in the Indo-Pacific region, yet the phylogenetic relationships and taxonomy of Pampus remain contentious. Here, we sampled 151 specimens, representing all known species of the genus, as well as two outgroup species (two families). We collected sequences from 17,292 single-copy nuclear coding loci using target-gene enrichment and Illumina sequencing for a subset samples of P. echinogaster, P. argenteus, P. cinereus, P. liuorum, P. chinensis, P. minor, and P. punctatissimus, which were carefully examined according to their species descriptions. Concatenated gene tree and species tree analyses resulted in identical and highly supported phylogenies, in which P. argenteus was sister to P. minor in one clade and P. cinereus sister to P. chinensis and P. punctatissimus in the other clade. Phylogenetic reconstruction using sequences of cytochrome c oxidase subunit I (COI) collected by us and those retrieved from NCBI suggests extensive misidentification of Pampus species in the NCBI database. We also measured morphological characters of each species as well as observed their osteological structure using micro-CT. Both molecular and morphological results suggest that P. echinogaster is a synonym of P. argenteus, and P. liuorum is a synonym of P. cinereus. Pampus cinereus from China, Bangladesh and an uncertain origin were grouped into three clades according to their sampling localities, but we could not find decisive morphological characters to describe the “cryptic species” of P. cinereus. Finally, based on the results of the molecular analyses and morphological reexamination, we created an identification key for the genus of Pampus.

Trypanosoma janseni n. sp. (Trypanosomatida: Trypanosomatidae) isolated from Didelphis aurita (Mammalia: Didelphidae) in the Atlantic Rainforest of Rio de Janeiro, Brazil: integrative taxonomy and phylogeography within the Trypanosoma cruzi clade.

BACKGROUND Didelphis spp. are a South American marsupial species that are among the most ancient hosts for the Trypanosoma spp.OBJECTIVESWe characterise a new species (Trypanosoma janseni n. sp.) isolated from the spleen and liver tissues of Didelphis aurita in the Atlantic Rainforest of Rio de Janeiro, Brazil.METHODSThe parasites were isolated and a growth curve was performed in NNN and Schneider's media containing 10% foetal bovine serum. Parasite morphology was evaluated via light microscopy on Giemsa-stained culture smears, as well as scanning and transmission electron microscopy. Molecular taxonomy was based on a partial region (737-bp) of the small subunit (18S) ribosomal RNA gene and 708 bp of the nuclear marker, glycosomal glyceraldehyde-3-phosphate dehydrogenase (gGAPDH) genes. Maximum likelihood and Bayesian inference methods were used to perform a species coalescent analysis and to generate individual and concatenated gene trees. Divergence times among species that belong to the T. cruzi clade were also inferred.FINDINGS In vitro growth curves demonstrated a very short log phase, achieving a maximum growth rate at day 3 followed by a sharp decline. Only epimastigote forms were observed under light and scanning microscopy. Transmission electron microscopy analysis showed structures typical to Trypanosoma spp., except one structure that presented as single-membraned, usually grouped in stacks of three or four. Phylogeography analyses confirmed the distinct species status of T. janseni n. sp. within the T. cruzi clade. Trypanosoma janseni n. sp. clusters with T. wauwau in a well-supported clade, which is exclusive and monophyletic. The separation of the South American T. wauwau + T. janseni coincides with the separation of the Southern Super Continent.CONCLUSIONSThis clade is a sister group of the trypanosomes found in Australian marsupials and its discovery sheds light on the initial diversification process based on what we currently know about the T. cruzi clade.

Evolutionary analysis of Chironius snakes unveils cryptic diversity and provides clues to diversification in the Neotropics

Recent hypotheses to explain tropical diversity involves the Neogene and Quaternary geoclimatic dynamics, but the absence of unambiguous data permitting the choice between alternative hypotheses makes a general theory for the origin of tropical biodiversity far to be achieved. The occurrence of Chironius snakes in well-defined biogeographical regions led us to adopt Chironius as a model to unveil patterns of vertebrate diversification in the Neotropics. Here, we used molecular markers and records on geographic distribution to investigate Chironius evolution and, subsequently, providing hints on diversification in the Neotropics. To avoid analyzing nominal species that do not constitute exclusive evolutionary lineages, we firstly conducted a species delimitation study prior to carrying out the species distribution modeling analysis. We generated 161 sequences of 12S, 16S, c-mos and rag2 for 15 species and 50 specimens, and included additional data from GenBank yielding a matrix of 137 terminals, and performed the following evolutionary analyses: inference of a concatenated gene tree, estimation of gene divergence times, inference of the coalescent-based phylogeny of Chironius, estimation of effective population sizes and modeling potential distribution of species across the last millennia. We tested for species boundaries within Chironius by implementing a coalescent-based Bayesian species delimitation approach. Our analyses supported the monophyly of Chironius, although our findings underscored cryptic candidate species in C. flavolineatus and C. exoletus. The inferred timetree suggested that Chironius snakes have evolved in the early Miocene (ca. 20.2Mya) and began to diversify from the late Miocene to the early Pliocene, values that are much older than previously reported. Following genetic divergence of virtually all extant Chironius species investigated, the effective sizes of the populations have expanded when compared to their MRCAs. The evolutionary and SDM data from C. brazili and C. diamantina provided additional evidence of the origin of species in the Neotropics. We argue that temperature, instead of altitude, has been the major driving factor in the evolution of both species, and thus we present a case for the consequences of global warming.

Multilocus Sequence Analysis of Root Nodule Bacteria Associated with <i>Lupinus</i> spp. and <i>Glycine max</i>

Lupinus is known to form endophytic associations with both nodulating and non-nodulating bacteria. In this study, multilocus sequence analysis (MLSA) was used to analyze phylogenetic relationships among root nodule bacteria associated with Lupinus and soybean. Out of 17 bacterial strains analyzed, 13 strains isolated from root nodules of Lupinus spp. were obtained from the National Rhizobium Germplasm Resource Collection, USDA. Additionally, two strains of root-nodule bacteria isolated each from native lupinus and domestic soybean were examined. Sequences of the 16S rRNA gene and three house-keeping genes (atpD, dnaK and glnII) were used. All the reference genes were retrieved from the existing complete genome sequences only. The clustering of 12 of the strains was consistent among single and concatenated gene trees, but not USDA strains 3044, 3048, 3504, 3715, and 3060. According to the concatenated phylogeny, we suggest that USDA 3040, 3042, 3044, 3048, 3051, 3060, 3504, 3709 and 3715 are Bradyrhizobium, USDA 3063 and 3717 are Mesorhizobium, USDA 3043 is Burkholderia and USDA 3057a is Microvirga. The two strains isolated from native lupines in this study are Burkholderia and Rhizobium, whereas the two from domestic soybean are Bradyrhizobium. This study emphasizes the robustness of MLSA, the diversity of bacterial species that are capable of nodulating lupine and the substantial capability of Burkholderia spp. to colonize lupine root nodules.

Picking holes in traditional species delimitations: an integrative taxonomic reassessment of theParmotrema perforatumgroup (Parmeliaceae, Ascomycota)

Accurate species delimitation is important as species are a fundamental unit in ecological, evolutionary and conservation biology research. In lichenized fungi, species delimitation has been difficult due to a lack of taxonomically important characteristics and due to the limits of traditional, morphology-based species concepts. In this study we reassess the current taxonomy of the Parmotrema perforatum group, which recognizes six closely related species divided into three species pairs, each pair comprising one apotheciate (sexual) and one sorediate (asexual) species. Each pair is further characterized by a distinct combination of secondary metabolites. It was hypothesized that the three apotheciate species are reproductively isolated sibling species and that each sorediate species evolved once from the chemically identical apotheciate species. In this study, species boundaries were re-examined using an integrative approach incorporating morphological, chemical and molecular sequence data to delimit species boundaries. Phylogenetic trees were inferred from a seven-locus DNA sequence dataset using concatenated gene tree and coalescent-based species-tree inference methods. Furthermore, we employed a multi-species coalescent method to validate candidate species. Micromorphological measurements of conidia were found to be congruent with phylogenetic clusters. Each approach that we applied to the P. perforatum group consistently recovered four of the currently circumscribed species (P. perforatum, P. hypotropum, P. subrigidum and P. louisianae), whereas P. preperforatum and P. hypoleucinum were consistently combined and are thus interpreted as conspecific.

Evolution and taxonomy of native mesorhizobia nodulating medicinal Glycyrrhiza species in China

Previously, 159 bacterial strains were isolated from the root nodules of wild perennial Glycyrrhiza legume species grown on 40 sites in central and north-western China, in which 57 strains were classified as “true symbionts” belonging to the genus Mesorhizobium based on amplified fragment length polymorphism (AFLP) genomic fingerprinting and partial sequences of the 16S rRNA gene [20]. In the present work, the phylogeny of Glycyrrhiza nodulating mesorhizobia was further examined by multilocus sequence analysis (MLSA). The concatenated gene tree of three housekeeping genes (16S rRNA, recA, and rpoB) of 59 strains including the 29 mesorhizobial test strains and 30 type mesorhizobial species, was constructed applying the maximum likelihood method and Bayesian inference. In the concatenated gene tree, the 29 test strains were distributed in seven separate clades. Seventeen test strains clustered with Mesorhizobium tianshanense, Mesorhizobium temperatum, Mesorhizobium muleiense, and Mesorhizobium alhagi with high bootstrap support (BS>85%). Eight test strains did not cluster with any of the described Mesorhizobium species. Based on the results, we proposed these eight test strains might belong to a putative new species of the genus Mesorhizobium. The sequences of three accessory genes (nodA, nodC, and nifH) of the test strains were also analyzed and were compared with those of representatives of the 30 described mesorhizobial species. The results showed that mesorhizobia involved in symbiosis with Glycyrrhiza plants probably have acquired some genetic material from other rhizobia in co-evolution with Glycyrrhiza and other legume species.

Resolving evolutionary relationships in lichen-forming fungi using diverse phylogenomic datasets and analytical approaches.

Evolutionary histories are now being inferred from unprecedented, genome-scale datasets for a broad range of organismal groups. While phylogenomic data has helped in resolving a number of difficult, long-standing questions, constructing appropriate datasets from genomes is not straightforward, particularly in non-model groups. Here we explore the utility of phylogenomic data to infer robust phylogenies for a lineage of closely related lichen-forming fungal species. We assembled multiple, distinct nuclear phylogenomic datasets, ranging from ca. 25 Kb to 16.8 Mb and inferred topologies using both concatenated gene tree approaches and species tree methods based on the multispecies coalescent model. In spite of evidence for rampant incongruence among individual loci, these genome-scale datasets provide a consistent, well-supported phylogenetic hypothesis using both concatenation and multispecies coalescent approaches (ASTRAL-II and SVDquartets). However, the popular full hierarchical coalescent approach implemented in *BEAST provided inconsistent inferences, both in terms of nodal support and topology, with smaller subsets of the phylogenomic data. While comparable, well-supported topologies can be accurately inferred with only a small fraction of the overall genome, consistent results across a variety of datasets and methodological approaches provide reassurance that phylogenomic data can effectively be used to provide robust phylogenies for closely related lichen-forming fungal lineages.

Gene and whole genome analyses reveal that the mycobacterial strain JS623 is not a member of the species Mycobacterium smegmatis.

Unexpected differences were found between the genome of strain JS623, used in bioremediation studies, and the genome of strain mc(2) 155, a model organism for investigating basic biology of mycobacteria. Both strains are currently assigned in the databases to the species Mycobacterium smegmatis and, consequently, the environmental isolate JS623 is increasingly included as a representative of that species in comparative genome-based approaches aiming at identifying distinctive traits of the different members of the genus Mycobacterium. We applied traditional molecular taxonomic procedures--inference of single and concatenated gene trees--to re-evaluate the membership of both strains to the same species, adopting the latest accepted cut-off values for species delimitation. Additionally, modern whole genome-based in silico methods where performed in a comprehensive molecular phylogenetic analysis of JS623 and other members of the genus Mycobacterium. These analyses showed that all relevant genome parameters of JS623 clearly separate this strain from M. smegmatis. The strain JS623 should be corrected as Mycobacterium sp. not only in the literature but, even more importantly, in the database entries, as inclusion of the genome wrongly attributed to the M. smegmatis species in comparative studies will result in misleading conclusions.

Multilocus phylogeny and taxonomic revision of the Hemidactylus robustus species group (Reptilia, Gekkonidae) with descriptions of three new species from Yemen and Ethiopia

The gecko genus Hemidactylus, with its 132 currently recognized species, ranks among the most species-rich reptile genera. Recent phylogenetic studies disclosed unexpectedly high genetic variability and complex biogeographic history within its arid clade distributed in the Mediterranean, Northeast Africa, Middle East and the Indian subcontinent. Particularly, the species from the Arabian Peninsula have been lately the subject of many taxonomic revisions that have resulted in the descriptions of 16 new taxa. Yet not all detected cryptic lineages have been treated taxonomically and thoroughly investigated morphologically. Based on phylogenetic analyses of two mtDNA (12S, cytb) and four nDNA (cmos, mc1r, rag1, rag2) gene fragments of a total length of 4015 bp in combination with analysis of morphological characters, we reinvestigate the systematics of the H. robustus species group consisting of the widespread H. robustus and three undescribed species, two of which occur in Southwest Yemen and one in central Ethiopia. By comparing two phylogenetic inference methods, concatenated gene trees and species-tree estimation, we reconstruct the phylogeny of the H. robustus species group. The coalescent-based species-tree estimation resulted in different tree topology than the concatenation approach, being probably a result of incomplete lineage sorting of ancestral polymorphism, suggesting that the H. robustus species group is in a stage of incipient speciation. The degree of differentiation of the characters examined within the H. robustus species group allowed us to provide a redescription of H. robustus and formally describe three new species of Hemidactylus – H. adensis sp. nov. and H. mandebensis sp. nov. from Yemen and H. awashensis sp. nov. from Ethiopia.

Species delimitation in the Aleochara fucicola species complex (Coleoptera: Staphylinidae: Aleocharinae) and its phylogenetic relationships

To assess the accuracy of species delimitation and phylogenetic relationships of the Aleochara fucicola species complex, we performed molecular phylogenetic analyses. Detailed micromorphological characters were also examined using scanning electron microscopy (SEM). The molecular data set included two mitochondrial (COI and COII) and three nuclear protein‐coding genes (CAD, EF1‐α and wg) for 19 populations of five species. Significant discordance was found between mitochondrial and nuclear gene trees. Two species, A. puetzi (Assing) and A. segregata Yamamoto and Maruyama, were not separated in the mitochondrial gene trees, but clearly were differentiated in the nuclear and concatenated gene trees. The SEM data also supported the separation of these two species, but an analysis of genetic divergence data did not. Adaptation to extremely harsh environments might reduce morphological variation in the A. fucicola species complex during the colonization of seashores. We present a multilocus phylogeny of the species complex. It suggests that the ancestor of the A. fucicola species complex occurred along the southern coasts of Northeast Asia, followed by dispersals to northern coasts.

Coalescent Species Delimitation in Milksnakes (Genus Lampropeltis) and Impacts on Phylogenetic Comparative Analyses

Both gene-tree discordance and unrecognized diversity are sources of error for accurate estimation of species trees, and can affect downstream diversification analyses by obscuring the correct number of nodes, their density, and the lengths of the branches subtending them. Although the theoretical impact of gene-tree discordance on evolutionary analyses has been examined previously, the effect of unsampled and cryptic diversity has not. Here, we examine how delimitation of previously unrecognized diversity in the milksnake (Lampropeltis triangulum) and use of a species-tree approach affects both estimation of the Lampropeltis phylogeny and comparative analyses with respect to the timing of diversification. Coalescent species delimitation indicates that L. triangulum is not monophyletic and that there are multiple species of milksnake, which increases the known species diversity in the genus Lampropeltis by 40%. Both genealogical and temporal discordance occurs between gene trees and the species tree, with evidence that mitochondrial DNA (mtDNA) introgression is a main factor. This discordance is further manifested in the preferred models of diversification, where the concatenated gene tree strongly supports an early burst of speciation during the Miocene, in contrast to species-tree estimates where diversification follows a birth-death model and speciation occurs mostly in the Pliocene and Pleistocene. This study highlights the crucial interaction among coalescent-based phylogeography and species delimitation, systematics, and species diversification analyses.

Evidence for inter-specific recombination among the mitochondrial genomes of Fusarium species in the Gibberella fujikuroi complex

BackgroundThe availability of mitochondrial genomes has allowed for the resolution of numerous questions regarding the evolutionary history of fungi and other eukaryotes. In the Gibberella fujikuroi species complex, the exact relationships among the so-called “African”, “Asian” and “American” Clades remain largely unresolved, irrespective of the markers employed. In this study, we considered the feasibility of using mitochondrial genes to infer the phylogenetic relationships among Fusarium species in this complex. The mitochondrial genomes of representatives of the three Clades (Fusarium circinatum, F. verticillioides and F. fujikuroi) were characterized and we determined whether or not the mitochondrial genomes of these fungi have value in resolving the higher level evolutionary relationships in the complex.ResultsOverall, the mitochondrial genomes of the three species displayed a high degree of synteny, with all the genes (protein coding genes, unique ORFs, ribosomal RNA and tRNA genes) in identical order and orientation, as well as introns that share similar positions within genes. The intergenic regions and introns generally contributed significantly to the size differences and diversity observed among these genomes. Phylogenetic analysis of the concatenated protein-coding dataset separated members of the Gibberella fujikuroi complex from other Fusarium species and suggested that F. fujikuroi (“Asian” Clade) is basal in the complex. However, individual mitochondrial gene trees were largely incongruent with one another and with the concatenated gene tree, because six distinct phylogenetic trees were recovered from the various single gene datasets.ConclusionThe mitochondrial genomes of Fusarium species in the Gibberella fujikuroi complex are remarkably similar to those of the previously characterized Fusarium species and Sordariomycetes. Despite apparently representing a single replicative unit, all of the genes encoded on the mitochondrial genomes of these fungi do not share the same evolutionary history. This incongruence could be due to biased selection on some genes or recombination among mitochondrial genomes. The results thus suggest that the use of individual mitochondrial genes for phylogenetic inference could mask the true relationships between species in this complex.

Neogene‐dominated diversification in neotropical montane lichens: Dating divergence events in the lichen‐forming fungal genus Oropogon (Parmeliaceae)

Diversification in neotropical regions has been attributed to both Tertiary geological events and Pleistocene climatic fluctuations. However, the timing and processes driving speciation in these regions remain unexplored in many important groups. Here, we address the timing of diversification in the neotropical lichenized fungal genus Oropogon (Ascomycota) and assess traditional species boundaries. We analyzed sequence data from three loci to assess phenotypically circumscribed Oropogon species from the Oaxacan Highlands, Mexico. We provide a comparison of dated divergence estimates between concatenated gene trees and a calibrated multilocus species-tree using substitution rates for two DNA regions. We also compare estimates from a data set excluding ambiguously aligned regions and a data set including the hyper-variable regions in two ribosomal markers. Phylogenetic reconstructions were characterized by well-supported monophyletic clades corresponding to traditionally circumscribed species, with the exception of a single taxon. Divergence estimates indicate that most diversification of the sampled Oropogon species occurred throughout the Oligocene and Miocene, although diversification of a single closely related clade appears to have occurred during the late Pliocene and into the Pleistocene. Divergence estimates calculated from a data set with ambiguously aligned regions removed were much more recent than those from the full data set. Overall, our analyses place the majority of divergence events of Oropogon species from the Oaxacan Highlands within the Neogene and provide strong evidence that climatic changes during the Pleistocene were not a major factor driving speciation in the lichenized genus Oropogon in neotropical highlands.

Multilocus phylogeny of the lichen-forming fungal genus Melanohalea (Parmeliaceae, Ascomycota): Insights on diversity, distributions, and a comparison of species tree and concatenated topologies

Accurate species circumscriptions are central for many biological disciplines and have critical implications for ecological and conservation studies. An increasing body of evidence suggests that in some cases traditional morphology-based taxonomy have underestimated diversity in lichen-forming fungi. Therefore, genetic data play an increasing role for recognizing distinct lineages of lichenized fungi that it would otherwise be improbable to recognize using classical phenotypic characters. Melanohalea (Parmeliaceae, Ascomycota) is one of the most widespread and common lichen-forming genera in the northern Hemisphere. In this study, we assess traditional phenotype-based species boundaries, identify previously unrecognized species-level lineages and discuss biogeographic patterns in Melanohalea. We sampled 487 individuals worldwide, representing 18 of the 22 described Melanohalea species, and generated DNA sequence data from mitochondrial, nuclear ribosomal, and protein-coding markers. Diversity previously hidden within traditional species was identified using a genealogical concordance approach. We inferred relationships among sampled species-level lineages within Melanohalea using both concatenated phylogenetic methods and a coalescent-based multilocus species tree approach. Although lineages identified from genetic data are largely congruent with traditional taxonomy, we found strong evidence supporting the presence of previously unrecognized species in six of the 18 sampled taxa. Strong nodal support and overall congruence among independent loci suggest long-term reproductive isolation among most species-level lineages. While some Melanohalea taxa are truly widespread, a limited number of clades appear to have much more restricted distributional ranges. In most instances the concatenated gene tree and multilocus species tree approaches provided similar estimates of relationships. However, nodal support was generally higher in the phylogeny estimated from concatenated data, and relationships among taxa within one major clade were largely unresolved in the species tree. This study contributes to our understanding of diversity and evolution in common lichen-forming fungi by incorporating multiple locus sequence data to circumscribe morphologicallly cryptic lineages and infer relationships within a coalescent-based species tree approach.

Multilocus Sequence Analysis of Isolates of Achromobacter from Patients with Cystic Fibrosis Reveals Infecting Species Other than Achromobacter xylosoxidans

A multilocus sequence analysis (MLSA) scheme was developed for characterization of strains and species from the genus Achromobacter, which are increasingly recovered from patients with cystic fibrosis (CF). Five conserved housekeeping genes were selected for the MLSA, which was applied to a diverse collection of 77 strains originating from Europe, Asia, and South America and including type strains of the seven recognized Achromobacter species, six environmental strains, eight non-CF clinical strains, and 56 CF clinical strains. The discriminatory power of MLSA, based on 2,098 nucleotides (nt), was much superior to a 16S rRNA gene comparison based on 1,309 nt. Congruence was observed between single-gene trees and a concatenated gene tree. MLSA differentiated all seven current Achromobacter species and also demonstrated the presence of at least four novel potential species within the genus. CF isolates were predominantly Achromobacter xylosoxidans (64%), an undescribed Achromobacter species (18%), and Achromobacter ruhlandii (7%). A clone of Achromobacter, which has spread among patients from Danish CF centers in Aarhus and Copenhagen, was identified as Achromobacter ruhlandii. MLSA facilitates the specific identification of isolates of Achromobacter necessary for describing their role in clinical infections.

Miocene and Pliocene dominated diversification of the lichen-forming fungal genus Melanohalea (Parmeliaceae, Ascomycota) and Pleistocene population expansions

BackgroundFactors promoting diversification in lichen symbioses remain largely unexplored. While Pleistocene events have been important for driving diversification and affecting distributions in many groups, recent estimates suggest that major radiations within some genera in the largest clade of macrolichens (Parmeliaceae, Ascomycota) vastly predate the Pleistocene. To better understand the temporal placement and sequence of diversification events in lichens, we estimated divergence times in a common lichen-forming fungal genus, Melanohalea, in the Northern Hemisphere. Divergence times were estimated using both concatenated gene tree and coalescent-based multilocus species tree approaches to assess the temporal context of major radiation events within Melanohalea. In order to complement our understanding of processes impacting genetic differentiation, we also evaluated the effects of Pleistocene glacial cycles on population demographics of distinct Melanohalea lineages, differing in reproductive strategies.ResultsWe found that divergence estimates, from both concatenated gene tree and coalescent-based multilocus species tree approaches, suggest that diversification within Melanohalea occurred predominantly during the Miocene and Pliocene, although estimated of divergence times differed by up to 8.3 million years between the two methods. These results indicate that, in some cases, taxonomically diagnostic characters may be maintained among divergent lineages for millions of years. In other cases, similar phenotypic characters among non-sister taxa, including reproductive strategies, suggest the potential for convergent evolution due to similar selective pressures among distinct lineages. Our analyses provide evidence of population expansions predating the last glacial maximum in the sampled lineages. These results suggest that Pleistocene glaciations were not inherently unfavorable or restrictive for some Melanohalea species, albeit with apparently different demographic histories between sexually and vegetatively reproducing lineages.ConclusionsOur results contribute to the understanding of how major changes during the Miocene and Pliocene have been important in promoting diversification within common lichen-forming fungi in the northern Hemisphere. Additionally, we provide evidence that glacial oscillations have influenced current population structure of broadly distributed lichenized fungal species throughout the Holarctic.

Barcoding against a paradox? Combined molecular species delineations reveal multiple cryptic lineages in elusive meiofaunal sea slugs.

BackgroundMany marine meiofaunal species are reported to have wide distributions, which creates a paradox considering their hypothesized low dispersal abilities. Correlated with this paradox is an especially high taxonomic deficit for meiofauna, partly related to a lower taxonomic effort and partly to a high number of putative cryptic species. Molecular-based species delineation and barcoding approaches have been advocated for meiofaunal biodiversity assessments to speed up description processes and uncover cryptic lineages. However, these approaches show sensitivity to sampling coverage (taxonomic and geographic) and the success rate has never been explored on mesopsammic Mollusca.ResultsWe collected the meiofaunal sea-slug Pontohedyle (Acochlidia, Heterobranchia) from 28 localities worldwide. With a traditional morphological approach, all specimens fall into two morphospecies. However, with a multi-marker genetic approach, we reveal multiple lineages that are reciprocally monophyletic on single and concatenated gene trees in phylogenetic analyses. These lineages are largely concordant with geographical and oceanographic parameters, leading to our primary species hypothesis (PSH). In parallel, we apply four independent methods of molecular based species delineation: General Mixed Yule Coalescent model (GMYC), statistical parsimony, Bayesian Species Delineation (BPP) and Automatic Barcode Gap Discovery (ABGD). The secondary species hypothesis (SSH) is gained by relying only on uncontradicted results of the different approaches (‘minimum consensus approach’), resulting in the discovery of a radiation of (at least) 12 mainly cryptic species, 9 of them new to science, some sympatric and some allopatric with respect to ocean boundaries. However, the meiofaunal paradox still persists in some Pontohedyle species identified here with wide coastal and trans-archipelago distributions.ConclusionsOur study confirms extensive, morphologically cryptic diversity among meiofauna and accentuates the taxonomic deficit that characterizes meiofauna research. We observe for Pontohedyle slugs a high degree of morphological simplicity and uniformity, which we expect might be a general rule for meiofauna. To tackle cryptic diversity in little explored and hard-to-sample invertebrate taxa, at present, a combined approach seems most promising, such as multi-marker-barcoding (i.e., molecular systematics using mitochondrial and nuclear markers and the criterion of reciprocal monophyly) combined with a minimum consensus approach across independent methods of molecular species delineation to define candidate species.

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 18mm 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.

Inferring Species-Level Phylogenies and Taxonomic Distinctiveness Using Multilocus Data in Sistrurus Rattlesnakes

Phylogenetic relationships and taxonomic distinctiveness of closely related species and subspecies are most accurately inferred from data derived from multiple independent loci. Here, we apply several approaches for understanding species-level relationships using data from 18 nuclear DNA loci and 1 mitochondrial DNA locus within currently described species and subspecies of Sistrurus rattlesnakes. Collectively, these methods provide evidence that a currently described species, the massasauga rattlesnake (Sistrurus catenatus), consists of two well-supported clades, one composed of the two western subspecies (S. c. tergeminus and S. c. edwardsii) and the other the eastern subspecies (S. c. catenatus). Within pigmy rattlesnakes (S. miliarius), however, there is not strong support across methods for any particular grouping at the subspecific level. Monophyly based tests for taxonomic distinctiveness show evidence for distinctiveness of all subspecies but this support is strongest by far for the S. c. catenatus clade. Because support for the distinctiveness of S. c. catenatus is both strong and consistent across methods, and due to its morphological distinctiveness and allopatric distribution, we suggest that this subspecies be elevated to full species status, which has significant conservation implications. Finally, most divergence time estimates based upon a fossil-calibrated species tree are > 50% younger than those from a concatenated gene tree analysis and suggest that an active period of speciation within Sistrurus occurred within the late Pliocene/Pleistocene eras.