Multilocus Data Research Articles

Species Tree Estimation from Gene Trees by Minimizing Deep Coalescence and Maximizing Quartet Consistency: A Comparative Study and the Presence of Pseudo Species Tree Terraces.

Species tree estimation from multilocus data sets is extremely challenging, especially in the presence of gene tree heterogeneity across the genome due to incomplete lineage sorting (ILS). Summary methods have been developed which estimate gene trees and then combine the gene trees to estimate a species tree by optimizing various optimization scores. In this study, we have extended and adapted the concept of phylogenetic terraces to species tree estimation by "summarizing" a set of gene trees, where multiple species trees with distinct topologies may have exactly the same optimality score (i.e., quartet score, extra lineage score, etc.). We particularly investigated the presence and impacts of equally optimal trees in species tree estimation from multilocus data using summary methods by taking ILS into account. We analyzed two of the most popular ILS-aware optimization criteria: maximize quartet consistency (MQC) and minimize deep coalescence (MDC). Methods based on MQC are provably statistically consistent, whereas MDC is not a consistent criterion for species tree estimation. We present a comprehensive comparative study of these two optimality criteria. Our experiments, on a collection of data sets simulated under ILS, indicate that MDC may result in competitive or identical quartet consistency score as MQC, but could be significantly worse than MQC in terms of tree accuracy-demonstrating the presence and impacts of equally optimal species trees. This is the first known study that provides the conditions for the data sets to have equally optimal trees in the context of phylogenomic inference using summary methods. [Gene tree; incomplete lineage sorting; phylogenomic analysis, species tree; summary method.].

Fossil-calibrated time tree of Podarcis wall lizards provides limited support for biogeographic calibration models

Podarcis wall lizards are endemic to the Mediterranean Basin where they represent the predominant reptile group. Despite being extensively used as model organisms in evolutionary and ecological studies their phylogeny and historical biogeography are still incompletely resolved. Moreover, molecular clock calibrations used in wall lizard phylogeography are based on the assumption of vicariant speciation triggered by the abrupt Mediterranean Sea level rise at the end of the Messinian salinity crisis (MSC). However, the validity of this biogeographic calibration remains untested. In this study we inferred a robust time tree based on multilocus data and fossil calibrations using both gene concatenation and species-tree approaches and including models with gene-flow. We found five deeply divergent, geographically coherent, and well-supported clades comprising species from i) Iberian Peninsula and North Africa; ii) Western Mediterranean islands, iii) Sicilian and Maltese islands; and iv-v) Balkan region and Aegean islands. The mitochondrial tree shows some inconsistencies with the species tree that warrant future investigation. Diversification of main clades is estimated in a short time frame during the Middle Miocene and might have been associated with a period of global climate cooling with the establishment of a marked climatic zonation in Europe. Cladogenetic events within the main clades are scattered throughout the time tree, from the Late Miocene to the Early Pleistocene, suggesting that speciation events in wall lizards reflect a complex interplay between regional topography, climate and geological history rather than a shared major climatic or paleogeographic event. Our absolute time estimates, as well as a relative dating approach, demonstrate that the assumption of a causal link between sea-level rise at the end of the MSC and the diversification of many island endemics is not justified. This study reinforces the notion that multiple dispersal and vicariant events, at different time frames, are required to explain current allopatric distributions and to account for the historical assembly of Mediterranean biota, and cautions against the use of biogeographic calibrations based on the assumption of vicariance.

Intraspecific variation and taxonomy of northern pikaOchotona hyperborea(Mammalia, Lagomorpha)

AbstractRock-dwelling pikas of the genus Ochotona have fragmented montane geographical distribution presumably associated with geological and climatic events of Late Quaternary. Within this genus, the northern pika, O. hyperborea has the widest distribution. Taxonomic composition and subspecific structure of this species has been controversial during the past century. Sixteen available names were combined in different manners into a variable number of weakly overlapping subspecies by previous authors. A comparison of data sets from the entire distribution range of northern pikas is presented here for the first time, including: craniometric, bioacoustic, and multilocus genetic data. We examined 301 skulls, and the mitochondrial COI gene and two nuclear introns in 79 specimens, as well as the structure of alarm calls from 107 individuals. Our results show that the six subspecies within O. hyperborea correspond to the six genetic parapatric lineages. Five of these are hypothesized to be involved in hybridization at the edges of their distribution. Three acoustic races completely correspond with the six genetic lineages: each acoustic race consists of two lineages. Morphometric data do not display any phylogenetic signal in our study. The subspecies from the Mountains of Khabarovsk Territory, north of Amur River, is described de novo as O. h. fedoseevi ssp. n. The type locality of O. hyperborea is defined more accurately. A neotype is designated for O. h. cinereofusca to stabilize use of names of pikas from the Amur region.

One, two or three? Integrative species delimitation of short-range endemic Hemicycla species (Gastropoda: Helicidae) from the Canary Islands based on morphology, barcoding, AFLP and ddRADseq data

Hemicycla mascaensis and H. diegoi are short-range endemics that occur allopatrically in small areas in the Teno Mountains in the western part of Tenerife (Canary Islands). Both taxa have been recognised as distinct species based on differences in shell morphology and genital anatomy. Preliminary molecular analyses using mitochondrial markers suggested a potential paraphyly of H. diegoi with regard to H. mascaensis. We here use multilocus AFLP data and ddRADseq data as well as distribution data, data on shell morphology and genital anatomy to assess the status of these taxa using phylogenetic analyses, species tree reconstruction and molecular species delimitation based on the multispecies coalescent as implemented in BFD* and BPP in an integrative approach. Our analyses show that, based on the analysis of multilocus data, the two taxa are reciprocally monophyletic. Species delimitation methods, however, tend to recognise all investigated populations as distinct species, albeit neither lending unambiguous support to any of the species hypotheses. The comparison of the anatomy of distal genital organs further suggests differentiation within H. mascaensis. This highlights the need for a balanced weighting of arguments from different lines of evidence to determine species status and calls for cautious interpretations of the results of molecular species delimitation analyses, especially in organisms with low active dispersal capacities and expected distinct population structuring such as land snails. Taking all available evidence into account, we favour to recognise H. mascaensis and H. diegoi as distinct species, acknowledging, though, that the recognition of both taxa as subspecies (with possibly a third yet undescribed) would also be an option as morphological differentiation is within the limits of other land snail species that are traditionally subdivided into subspecies.

Inocybe brijunica sp. nov., a New Ectomycorrhizal Fungus from Mediterranean Croatia Revealed by Morphology and Multilocus Phylogenetic Analysis.

A new ectomycorrhizal species was discovered during the first survey of fungal diversity at Brijuni National Park (Croatia), which consists of 14 islands and islets. The National Park is located in the Mediterranean Biogeographical Region, a prominent climate change hot-spot. Inocybe brijunica sp. nov., from sect. Hysterices (Agaricales, Inocybaceae), is described based on morphology and multilocus phylogenetic data. The holotype collection was found at the edge between grassland and Quercus ilex forest with a few planted Pinus pinea trees, on Veli Brijun Island, the largest island of the archipelago. It is easily recognized by a conspicuous orange to orange–red–brown membranaceous surface layer located at or just above the basal part of the stipe. Other distinctive features of I. brijunica are the medium brown, radially fibrillose to rimose pileus; pale to medium brown stipe with fugacious cortina; relatively small, amygdaliform to phaseoliform, and smooth basidiospores, measuring ca. 6.5–9 × 4–5.5 µm; thick-walled, utriform, lageniform or fusiform pleurocystidia (lamprocystidia) with crystals and mostly not yellowing in alkaline solutions; cheilocystidia of two types (lamprocystidia and leptocystidia); and the presence of abundant caulocystidia only in the upper 2–3 mm of the stipe. Phylogenetic reconstruction of a concatenated dataset of the internal transcribed spacer region (ITS), the nuclear 28S rRNA gene (nrLSU), and the second largest subunit of RNA polymerase II (rpb2) resolved I. brijunica and I. glabripes as sister species.

Taxonomic and pathogenic characterization of Phytopythium species from West Azarbaijan, Iran, and description of two new species

ABSTRACT Phytopythium is a globally distributed genus, species of which occur in various ecological niches. Despite their importance as plant pathogens, information on species diversity and their distribution in West Azarbaijan Province of Iran is limited. This study aimed to investigate the taxonomy and plant pathogenicity of Phytopythium isolates recovered from soils in the province. A polyphasic taxonomy based on morphological, cultural, and multilocus sequence data revealed two new species, described here as P. babaiaharii and P. longitubum. In addition, a putatively new species and five known species, P. boreale, P. carbonicum, P. mercuriale, P. ostracodes, and P. vexans, were found in the studied region. Phytopythium ostracodes, P. mercuriale, and P. boreale were the three most frequent species isolated from soil, although P. mercuriale was only found in one field. Among the studied species, P. mercuriale, P. ostracodes, and P. vexans exhibited various levels of pathogenicity on sugar beet, sunflower, and tomato, and P. boreale was documented for first time as a plant pathogen. Phytopythium ostracodes was highly aggressive on sugar beet and sunflower but lowly aggressive on tomato, whereas P. vexans exhibited high aggressiveness toward the three crop plants. Both P. boreale and P. mercuriale were found to be lowly aggressive on the host plants. The results indicated that members of Phytopythium, particularly plant-pathogenic species, are common in arable soils of West Azarbaijan Province and can pose a threat to agricultural crops in the region.

Sinuicella denisonii, a new genus and species in thePeltigeraceaefrom western North America

AbstractThe new genusSinuicella, an early successional lichen, was found on bare soil in Oregon, USA. The thallus is minute fruticose, grey to nearly black, branching isotomic dichotomous, branches round, 20–90 μm wide in water mount. The cortex is composed of interlocking cells shaped like jigsaw puzzle pieces. Spores are hyaline, 1-septate, 25–40(–50) × 6.5–9(–11) μm. Maximum likelihood phylogenetic analyses on multilocus data sets, first spanning the entire orderPeltigeralesand then restricted toPeltigeraceaewith extended sampling fromSolorinaandPeltigera, revealed the placement ofSinuicellaoutside of currently recognized genera, sister toPeltigera, with high support. Based on the phylogenetic, morphological and ecological distinctness ofSinuicella, we formally introduce a new genus represented by the single speciesS. denisonii. The cyanobiont ofS. denisoniiisNostocfrom phylogroup XL, Clade 2, Subclade 3 based on therbcLXmarker.

Performance comparison of two reduced-representation based genome-wide marker-discovery strategies in a multi-taxon phylogeographic framework

Multi-locus genetic data are pivotal in phylogenetics. Today, high-throughput sequencing (HTS) allows scientists to generate an unprecedented amount of such data from any organism. However, HTS is resource intense and may not be accessible to wide parts of the scientific community. In phylogeography, the use of HTS has concentrated on a few taxonomic groups, and the amount of data used to resolve a phylogeographic pattern often seems arbitrary. We explore the performance of two genetic marker sampling strategies and the effect of marker quantity in a comparative phylogeographic framework focusing on six species (arthropods and plants). The same analyses were applied to data inferred from amplified fragment length polymorphism fingerprinting (AFLP), a cheap, non-HTS based technique that is able to straightforwardly produce several hundred markers, and from restriction site associated DNA sequencing (RADseq), a more expensive, HTS-based technique that produces thousands of single nucleotide polymorphisms. We show that in four of six study species, AFLP leads to results comparable with those of RADseq. While we do not aim to contest the advantages of HTS techniques, we also show that AFLP is a robust technique to delimit evolutionary entities in both plants and animals. The demonstrated similarity of results from the two techniques also strengthens biological conclusions that were based on AFLP data in the past, an important finding given the wide utilization of AFLP over the last decades. We emphasize that whenever the delimitation of evolutionary entities is the central goal, as it is in many fields of biodiversity research, AFLP is still an adequate technique.

Phylogeography and demographic expansion in the widely distributed horned passalus beetle, Odontotaenius disjunctus (coleoptera: Passalidae)

Dynamic climatic oscillations during the Pleistocene had profound effects on the distribution of species across North America. Although the role of historical climate change on speciation remains controversial, the impact on genetic variation within species has been well documented. Analyses of mtDNA sequences from the cytochrome oxidase I gene (911 bp) for 115 individuals of Odontotaenius disjunctus was combined with ecological niche modelling (ENM) to infer the demographic and population differentiation scenarios under present and past conditions. We inferred three lineages that diverged during the Pleistocene and replace each other geographically across the eastern United States. One of these lineages traverses previously identified genetic barriers for terrestrial animals including the Mississippi and Apalachicola Rivers and the Appalachian Mountains. We observed overlapping ranges between two haplotype groups as well as a region of secondary contact associated with ecological transition zone in northern Florida. The two continental lineages depict a genetic signature of a recent population increase associated with expanding niche envelope, whereas the clade restricted to peninsular Florida shows stable populations in a shrinking niche envelope. Given the lack of ecological separation, overlapping distribution of haplogroups and the presence of secondary contact zones, the taxonomic status of these lineages must await robust testing using multilocus DNA data to assess species boundaries.

Gene Tree Discord, Simplex Plots, and Statistical Tests under the Coalescent.

A simple graphical device, the simplex plot of quartet concordance factors, is introduced to aid in the exploration of a collection of gene trees on a common set of taxa. A single plot summarizes all gene tree discord and allows for visual comparison to the expected discord from the multispecies coalescent model (MSC) of incomplete lineage sorting on a species tree. A formal statistical procedure is described that can quantify the deviation from expectation for each subset of four taxa, suggesting when the data are not in accord with the MSC, and thus that either gene tree inference error is substantial or a more complex model such as that on a network may be required. If the collection of gene trees is in accord with the MSC, the plots reveal when substantial incomplete lineage sorting is present. Applications to both simulated and empirical multilocus data sets illustrate the insights provided. [Gene tree discordance; hypothesis test; multispecies coalescent model; quartet concordance factor; simplex plot; species tree].

MixIndependR: a R package for statistical independence testing of loci in database of multi-locus genotypes

BackgroundMulti-locus genotype data are widely used in population genetics and disease studies. In evaluating the utility of multi-locus data, the independence of markers is commonly considered in many genomic assessments. Generally, pairwise non-random associations are tested by linkage disequilibrium; however, the dependence of one panel might be triplet, quartet, or other. Therefore, a compatible and user-friendly software is necessary for testing and assessing the global linkage disequilibrium among mixed genetic data.ResultsThis study describes a software package for testing the mutual independence of mixed genetic datasets. Mutual independence is defined as no non-random associations among all subsets of the tested panel. The new R package “mixIndependR” calculates basic genetic parameters like allele frequency, genotype frequency, heterozygosity, Hardy–Weinberg equilibrium, and linkage disequilibrium (LD) by mutual independence from population data, regardless of the type of markers, such as simple nucleotide polymorphisms, short tandem repeats, insertions and deletions, and any other genetic markers. A novel method of assessing the dependence of mixed genetic panels is developed in this study and functionally analyzed in the software package. By comparing the observed distribution of two common summary statistics (the number of heterozygous loci [K] and the number of share alleles [X]) with their expected distributions under the assumption of mutual independence, the overall independence is tested.ConclusionThe package “mixIndependR” is compatible to all categories of genetic markers and detects the overall non-random associations. Compared to pairwise disequilibrium, the approach described herein tends to have higher power, especially when number of markers is large. With this package, more multi-functional or stronger genetic panels can be developed, like mixed panels with different kinds of markers. In population genetics, the package “mixIndependR” makes it possible to discover more about admixture of populations, natural selection, genetic drift, and population demographics, as a more powerful method of detecting LD. Moreover, this new approach can optimize variants selection in disease studies and contribute to panel combination for treatments in multimorbidity. Application of this approach in real data is expected in the future, and this might bring a leap in the field of genetic technology.AvailabilityThe R package mixIndependR, is available on the Comprehensive R Archive Network (CRAN) at: https://cran.r-project.org/web/packages/mixIndependR/index.html.

Phylogenomic re-evaluation of Triaenonychoidea (Opiliones : Laniatores), and systematics of Triaenonychidae, including new families, genera and species

The Opiliones superfamily Triaenonychoidea currently includes two families, the monogeneric New Zealand–endemic Synthetonychiidae Forster, 1954 and Triaenonychidae Sørensen, 1886, a diverse family distributed mostly throughout the temperate Gondwanan terranes, with ~110 genera and ~500 species and subspecies currently described. Traditionally, Triaenonychidae has been divided into subfamilies diagnosed by very few morphological characters largely derived from the troublesome ‘Roewerian system’ of morphology, and classifications based on this system led to many complications. Recent research within Triaenonychoidea using morphology and traditional multilocus data has shown multiple deeply divergent lineages, non-monophyly of Triaenonychidae, and non-monophyly of subfamilies, necessitating a revision based on phylogenomic data. We used sequence capture of ultraconserved elements across 164 samples to create a 50% taxon occupancy matrix with 704 loci. Using phylogenomic and morphological examinations, we explored family-level relationships within Triaenonychoidea, including describing two new families: (1) Lomanellidae Mendes & Derkarabetian, fam. nov., consisting of Lomanella Pocock, 1903, and a newly described genus Abaddon Derkarabetian & Baker, gen. nov. with one species, A. despoliator Derkarabetian, sp. nov.; and (2) the elevation to family of Buemarinoidae Karaman, 2019, consisting of Buemarinoa Roewer, 1956, Fumontana Shear, 1977, Flavonuncia Lawrence, 1959, and a newly described genus Turonychus Derkarabetian, Prieto & Giribet, gen. nov., with one species, T. fadriquei Derkarabetian, Prieto & Giribet, sp. nov. With our dataset we also explored phylogenomic relationships within Triaenonychidae with an extensive taxon set including samples representing ~80% of the genus-level diversity. Based on our results we (1) discuss systematics of this family including the historical use of subfamilies, (2) reassess morphology in the context of our phylogeny, (3) hypothesise placement for all unsampled genera, (4) highlight lineages most in need of taxonomic revision, and (5) provide an updated species-level checklist. Aside from describing new taxa, our study provides the phylogenomic context necessary for future evolutionary and systematic research across this diverse lineage. ZooBank Registration: urn:lsid:zoobank.org:pub:81683834-98AB-43AA-B25A-C28C6A404F41

Integrative Description of Pseudechiniscus lalitae sp. nov. (Tardigrada: Heterotardigrada: Echiniscidae) from the Azores Archipelago (Portugal)

We applied an integrative taxonomy approach to describe a new tardigrade species from Pico Island (the Azores Archipelago). Analysis was based on morphological and morphometrical data (using both light and scanning electron microscopy) and combined with multilocus molecular data (nuclear barcode sequence ITS-2 and mitochondrial COI sequence). Molecular data allowed us to determine that genetically (based on COI sequences) Pseudechiniscus lalitaesp. nov. is most similar to Pse. aff. facettalis and the least similar to Pse. indistinctus. The new species differs morphologically from other Pseudechiniscus species, mostly, by the ventral sculpture. Moreover, it differs from other species from suillus-facettalis complex also by some other more subtle morphological and morphometrical characters.

Polynomial-Time Statistical Estimation of Species Trees Under Gene Duplication and Loss.

Phylogenomics-the estimation of species trees from multilocus data sets-is a common step in many biological studies. However, this estimation is challenged by the fact that genes can evolve under processes, including incomplete lineage sorting (ILS) and gene duplication and loss (GDL), that make their trees different from the species tree. In this article, we address the challenge of estimating the species tree under GDL. We show that species trees are identifiable under a standard stochastic model for GDL, and that the polynomial-time algorithm ASTRAL-multi, a recent development in the ASTRAL suite of methods, is statistically consistent under this GDL model. We also provide a simulation study evaluating ASTRAL-multi for species tree estimation under GDL.

Assessing the fit of the multi-species network coalescent to multi-locus data.

With growing genome-wide molecular datasets from next-generation sequencing, phylogenetic networks can be estimated using a variety of approaches. These phylogenetic networks include events like hybridization, gene flow or horizontal gene transfer explicitly. However, the most accurate network inference methods are computationally heavy. Methods that scale to larger datasets do not calculate a full likelihood, such that traditional likelihood-based tools for model selection are not applicable to decide how many past hybridization events best fit the data. We propose here a goodness-of-fit test to quantify the fit between data observed from genome-wide multi-locus data, and patterns expected under the multi-species coalescent model on a candidate phylogenetic network. We identified weaknesses in the previously proposed TICR test, and proposed corrections. The performance of our new test was validated by simulations on real-world phylogenetic networks. Our test provides one of the first rigorous tools for model selection, to select the adequate network complexity for the data at hand. The test can also work for identifying poorly inferred areas on a network. Software for the goodness-of-fit test is available as a Julia package at https://github.com/cecileane/QuartetNetworkGoodnessFit.jl. Supplementary data are available at Bioinformatics online.

Source attributed case-control study of campylobacteriosis in New Zealand

BackgroundFollowing an initial reduction in human campylobacteriosis in New Zealand after the implementation of poultry food chain-focused interventions during 2006–2008, further decline has been relatively small.We report a year-long study of notified campylobacteriosis cases, incorporating a case control study combined with a source attribution study. The purpose was to generate up-to-date evidence on the relative contributions of different sources of campylobacteriosis in New Zealand. MethodsThe study approach included:• A case-control study of notified cases (aged six months or more) sampled in a major urban centre (Auckland, every second case) and a mixed urban/rural area (Manawatū/Whanganui, every case), between 12 March 2018 and 11 March 2019.• Source attribution of human campylobacteriosis cases sampled from these two regions over the study period by modelling of multilocus sequence typing data of Campylobacter jejuni and C. coli isolates from faecal samples of notified human cases and relevant sources (poultry, cattle, sheep). ResultsMost cases (84%) were infected with strains attributed to a poultry source, while 14% were attributed to a cattle source. Approximately 90% of urban campylobacteriosis cases were attributed to poultry sources, compared to almost 75% of rural cases.Poultry consumption per se was not identified as a significant risk factor. However specific risk factors related to poultry meat preparation and consumption did result in statistically significantly elevated odds ratios. ConclusionsThe overall findings combining source attribution and analysis of specific risk factors indicate that poultry meat remains a dominant pathway for exposure and infection.

Genetic Divergence Across Glacial Refugia Despite Interglacial Gene Flow in a Crested Newt

MtDNA-based phylogeography has illuminated the impact of the Pleistocene Ice Age on species distribution dynamics and the build-up of genetic divergence. The well-known shortcomings of mtDNA in biogeographical inference can be compensated by integrating multilocus data and species distribution modelling into phylogeography. We re-visit the phylogeography of the Italian crested newt (Triturus carnifex), a species distributed in two of Europe’s main glacial refugia, the Balkan and Italian Peninsulas. While a new 51 nuclear DNA marker dataset supports the existence of three lineages previously suggested by mtDNA (Balkan, northern Italy and southern Italy), the nuclear DNA dataset also provides improved resolution where these lineages have obtained secondary contact. We observe geographically restricted admixture at the contact between the Balkan and northern Italy gene pools and identify a potential mtDNA ghost lineage here. At the contact between the northern and southern Italy gene pools we find admixture over a broader area, as well as asymmetric mtDNA introgression. Our species distribution model is in agreement with a distribution restricted to distinct refugia during Pleistocene glacial cycles and postglacial expansion with secondary contact. Our study supports: (1) the relevance of the north-western Balkan Peninsula as a discrete glacial refugium; (2) the importance of north-eastern Italy and the northern Apennine as suture zones; and (3) the applicability of a refugia-within-refugia scenario within the Italian Peninsula.

A systematic revision of Draculoides (Schizomida: Hubbardiidae) of the Pilbara, Western Australia, Part I: the Western Pilbara.

The schizomid fauna of mainland Australia currently comprises 60 species within seven named genera, of which five are endemic to the continent: Attenuizomus Harvey, 2000, Brignolizomus Harvey, 2000, Draculoides Harvey, 1992, Julattenius Harvey, 1992, Notozomus Harvey, 2000. Most Australian schizomids have been described from eastern and northern Australia, but there is also a significant subterranean fauna that has been found in hypogean habitats in the semi-arid Pilbara region of Western Australia. The vast majority of these species can be assigned to the genus Draculoides and this study is the first in a proposed series to revise this highly diverse genus. We treat the species found in the western Pilbara region, which includes 13 new species and 13 previously named species, using morphological characters and multi-locus sequence data. We also incorporate a molecular "mini-barcode" approach for COI, 12S and ITS2 to diagnose the new species. The new species are named: Draculoides akashae Abrams and Harvey, n. sp., D. belalugosii Abrams and Harvey, n. sp., D. carmillae Abrams and Harvey, n. sp., D. christopherleei Abrams and Harvey, n. sp., D. claudiae Abrams and Harvey, n. sp., D. immortalis Abrams and Harvey, n. sp., D. karenbassettae Abrams and Harvey, n. sp., D. mckechnieorum Abrams and Harvey, n. sp., D. minae Abrams and Harvey, n. sp., D. noctigrassator Abrams and Harvey, n. sp., D. nosferatu Abrams and Harvey, n. sp., D. piscivultus Abrams and Harvey, n. sp. and D. warramboo Abrams and Harvey, n. sp. We also provide the first descriptions of males of D. anachoretus (Harvey, Berry, Edward and Humphreys, 2008) and D. gnophicola (Harvey, Berry, Edward and Humphreys, 2008). All of the new species are subterranean-dwelling, short-range endemic species that occur in regions subject to mining activities, rendering them of high conservation significance.

A globally threatened shark, Carcharias taurus, shows no population decline in South Africa

Knowledge about the demographic histories of natural populations helps to evaluate their conservation status, and potential impacts of natural and anthropogenic pressures. In particular, estimates of effective population size obtained through molecular data can provide useful information to guide management decisions for vulnerable populations. The spotted ragged-tooth shark, Carcharias taurus (also known as the sandtiger or grey nurse shark), is widely distributed in warm-temperate and subtropical waters, but has suffered severe population declines across much of its range as a result of overexploitation. Here, we used multilocus genotype data to investigate the demographic history of the South African C. taurus population. Using approximate Bayesian computation and likelihood-based importance sampling, we found that the population underwent a historical range expansion that may have been linked to climatic changes during the late Pleistocene. There was no evidence for a recent anthropogenic decline. Together with census data suggesting a stable population, these results support the idea that fishing pressure and other threats have so far not been detrimental to the local C. taurus population. The results reported here indicate that South Africa could possibly harbour the last remaining, relatively pristine population of this widespread but vulnerable top predator.

AllCoPol: inferring allele co-ancestry in polyploids

BackgroundInferring phylogenetic relationships of polyploid species and their diploid ancestors (leading to reticulate phylogenies in the case of an allopolyploid origin) based on multi-locus sequence data is complicated by the unknown assignment of alleles found in polyploids to diploid subgenomes. A parsimony-based approach to this problem has been proposed by Oberprieler et al. (Methods Ecol Evol 8:835–849, 2017), however, its implementation is of limited practical value. In addition to previously identified shortcomings, it has been found that in some cases, the obtained results barely satisfy the applied criterion. To be of better use to other researchers, a reimplementation with methodological refinement appears to be indispensable.ResultsWe present the AllCoPol package, which provides a heuristic method for assigning alleles from polyploids to diploid subgenomes based on the Minimizing Deep Coalescences (MDC) criterion in multi-locus sequence datasets. An additional consensus approach further allows to assess the confidence of phylogenetic reconstructions. Simulations of tetra- and hexaploids show that under simplifying assumptions such as completely disomic inheritance, the topological errors of reconstructed phylogenies are similar to those of MDC species trees based on the true allele partition.ConclusionsAllCoPol is a Python package for phylogenetic reconstructions of polyploids offering enhanced functionality as well as improved usability. The included methods are supplied as command line tools without the need for prior programming knowledge.