Genetic bias in repeated evolution of pigment loss in cave populations of the Asellus aquaticus species complex.

Similar phenotypes can evolve repeatedly under the same evolutionary pressures. A compelling example is the evolution of pigment loss and eye loss in cave-dwelling animals. While specific genomic regions or genes associated with these phenotypes have been identified in model species, it remains uncertain whether a bias towards particular genetic mechanisms exists. An isopod crustacean, Asellus aquaticus, is an ideal model organism to investigate this phenomenon. It inhabits surface freshwaters throughout Europe but has colonized groundwater on multiple independent occasions and evolved several cave populations with distinct ecomorphology. Previous studies have demonstrated that three different cave populations utilized common genetic regions, potentially the same genes, in the evolution of pigment and eye loss. Expanding on this, we conducted analysis on two additional cave populations, distinct either phylogenetically or biogeographically from those previously examined. We generated F2 hybrids from cave × surface crosses and tested phenotype-genotype associations, as well as conducted complementation tests by crossing individuals from different cave populations. Our findings revealed that pigment loss and orange eye pigment in additional cave populations were associated with the same genomic regions as observed in the three previously tested cave populations. Moreover, the lack of complementation across all cross combinations suggests that the same gene likely drives pigment loss. These results substantiate a genetic bias in the recurrent evolution of pigment loss in this model system. Future investigations should focus on the cause behind this bias, possibly arising from allele recruitment from ancestral surface populations' genetic variation or advantageous allele effects via pleiotropy.

Phylogenetic relationships and species delimitation in Haemopis (Annelida: Hirudinea: Haemopidae).

The Holarctic leech genus Haemopis currently includes 11 species, all of which are macrophagous, as opposed to their more infamous bloodfeeding counterparts among hirudiniform leeches. In spite of their ecological importance as fish food and predators of freshwater invertebrates, there is a paucity of data regarding morphology and genetic variation that might guide future identification efforts for members of the genus. The lack of detailed descriptions of distinguishing morphological features, coupled with the absence of a robust phylogenetic hypothesis for the genus, have conspired to prevent meaningful inferences on the natural history of the group. In an attempt to remedy this, we present new genetic (using COI, 12S rDNA, 28S rDNA and 18S rDNA) data for the majority of the known species diversity within the genus in order to both infer a phylogenetic hypothesis and to introduce authoritative DNA barcodes for the newly collected species. The potential of these barcodes is increased through rigorous morphological investigations of the specimens, with comparisons to the original literature. Our resulting phylogenetic hypothesis is agnostic as to the geographic origin of the genus, with equal probability afforded to both a Nearctic and Palearctic origin. Beyond this, we show that there is a strong tendency towards a barcoding gap within the genus, but that a distinct gap is lacking due to the relatively high genetic variation found within H. marmorata. Taken together, our results shed light on species delimitation within, and evolutionary history of, this often-neglected group of leeches.

From Cave Dragons to Genomics: Advancements in the Study of Subterranean Tetrapods.

Throughout most of the kingdom Animalia, evolutionary transitions from surface life to a life permanently bound to caves and other subterranean habitats have occurred innumerous times. Not so in tetrapods, where a mere 14 cave-obligate species-all plethodontid and proteid salamanders-are known. We discuss why cave tetrapods are so exceptional and why only salamanders have made the transition. Their evolution follows predictable and convergent, albeit independent pathways. Among the many known changes associated with transitions to subterranean life, eye degeneration, starvation resistance, and longevity are especially relevant to human biomedical research. Recently, sequences of salamander genomes have become available opening up genomic research for cave tetrapods. We discuss new genomic methods that can spur our understanding of the evolutionary mechanisms behind convergent phenotypic change, the relative roles of selective and neutral evolution, cryptic species diversity, and data relevant for conservation such as effective population size and demography.

Genetic basis of eye and pigment loss in the cave crustacean, Asellus aquaticus.

Understanding the process of evolution is one of the great challenges in biology. Cave animals are one group with immense potential to address the mechanisms of evolutionary change. Amazingly, similar morphological alterations, such as enhancement of sensory systems and the loss of eyes and pigmentation, have evolved multiple times in a diverse assemblage of cave animals. Our goal is to develop an invertebrate model to study cave evolution so that, in combination with a previously established vertebrate cave system, we can address genetic questions concerning evolutionary parallelism and convergence. We chose the isopod crustacean, Asellus aquaticus, and generated a genome-wide linkage map for this species. Our map, composed of 117 markers, of which the majority are associated with genes known to be involved in pigmentation, eye, and appendage development, was used to identify loci of large effect responsible for several pigmentation traits and eye loss. Our study provides support for the prediction that significant morphological change can be mediated through one or a few genes. Surprisingly, we found that within population variability in eye size occurs through multiple mechanisms; eye loss has a different genetic basis than reduced eye size. Similarly, again within a population, the phenotype of albinism can be achieved by two different genetic pathways--either by a recessive genotype at one locus or doubly recessive genotypes at two other loci. Our work shows the potential of Asellus for studying the extremes of parallel and convergent evolution-spanning comparisons within populations to comparisons between vertebrate and arthropod systems.

Generation of genome-wide SNP markers from minimally invasive sampling in endangered animals and applications in species ecology and conservation.

High-density genotyping methods have revolutionized the field of population and conservation genetics in the past decade. To exploit the technological and analytical advances in the field, access to high-quality genetic material is a key component. However, access to such samples in endangered and rare animals is often challenging or even impossible. Here, we used a minimally invasive sampling method (MIS) in the endangered cave salamander Proteus anguinus, the olm, to generate thousands of genetic markers using ddRADseq for population and conservation genomic analyses. Using tail clips and MIS skin swabs taken from the same individual, we investigated genotyping data properties of the two different sampling types. We found that sufficient DNA can be extracted from swab samples to generate up to 200,000 polymorphic SNPs in divergent Proteus lineages. Swab and tissue samples were highly reproducible exhibiting low SNP genotyping error rates. We found that SNPs were most frequently (~50%) located within genic regions, while the rest mapped to mostly flanking regions of repetitive DNA. The vast majority of DNA recovered from swabbing was host DNA. However, a fraction of DNA recovered from swabs contained additional ecological information on the species, including eDNA from the surrounding environment and bacterial skin fauna. Most exogenous DNA recovered from swabs were bacteria (~80%), followed by vertebrates (~20%). Our results demonstrate that MIS can be used to (i) generate tens of thousands of ddRADseq markers for conservation and population genomic analyses and (ii) inform on the species health status and ecology from exogenous DNA.

A molecular phylogeny of nephilid spiders: evolutionary history of a model lineage.

The pantropical orb web spider family Nephilidae is known for the most extreme sexual size dimorphism among terrestrial animals. Numerous studies have made Nephilidae, particularly Nephila, a model lineage in evolutionary research. However, a poorly understood phylogeny of this lineage, relying only on morphology, has prevented thorough evolutionary syntheses of nephilid biology. We here use three nuclear and five mitochondrial genes for 28 out of 40 nephilid species to provide a more robust nephilid phylogeny and infer clade ages in a fossil-calibrated Bayesian framework. We complement the molecular analyses with total evidence analysis including morphology. All analyses find strong support for nephilid monophyly and exclusivity and the monophyly of the genera Herennia and Clitaetra. The inferred phylogenetic structure within Nephilidae is novel and conflicts with morphological phylogeny and traditional taxonomy. Nephilengys species fall into two clades, one with Australasian species (true Nephilengys) as sister to Herennia, and another with Afrotropical species (Nephilingis Kuntner new genus) as sister to a clade containing Clitaetra plus most currently described Nephila. Surprisingly, Nephila is also diphyletic, with true Nephila containing N. pilipes+N. constricta, and the second clade with all other species sister to Clitaetra; this "Nephila" clade is further split into an Australasian clade that also contains the South American N. sexpunctata and the Eurasian N. clavata, and an African clade that also contains the Panamerican N. clavipes. An approximately unbiased test constraining the monophyly of Nephilengys, Nephila, and Nephilinae (Nephila, Nephilengys, Herennia), respectively, rejected Nephilengys monophyly, but not that of Nephila and Nephilinae. Further data are therefore necessary to robustly test these two new, but inconclusive findings, and also to further test the precise placement of Nephilidae within the Araneoidea. For divergence date estimation we set the minimum bound for the stems of Nephilidae at 40 Ma and of Nephila at 16 Ma to accommodate Palaeonephila from Baltic amber and Dominican Nephila species, respectively. We also calibrated and dated the phylogeny under three different interpretations of the enigmatic 165 Ma fossil Nephila jurassica, which we suspected based on morphology to be misplaced. We found that by treating N. jurassica as stem Nephila or nephilid the inferred clade ages were vastly older, and the mitochondrial substitution rates much slower than expected from other empirical spider data. This suggests that N. jurassica is not a Nephila nor a nephilid, but possibly a stem orbicularian. The estimated nephilid ancestral age (40-60 Ma) rejects a Gondwanan origin of the family as most of the southern continents were already split at that time. The origin of the family is equally likely to be African, Asian, or Australasian, with a global biogeographic history dominated by dispersal events. A reinterpretation of web architecture evolution suggests that a partially arboricolous, asymmetric orb web with a retreat, as exemplified by both groups of "Nephilengys", is plesiomorphic in Nephilidae, that this architecture was modified into specialized arboricolous webs in Herennia and independently in Clitaetra, and that the web became aerial, gigantic, and golden independently in both "Nephila" groups. The new topology questions previously hypothesized gradual evolution of female size from small to large, and rather suggests a more mosaic evolutionary pattern with independent female size increases from medium to giant in both "Nephila" clades, and two reversals back to medium and small; combined with male size evolution, this pattern will help detect gross evolutionary events leading to extreme sexual size dimorphism, and its morphological and behavioral correlates.

Wolf genetic diversity compared across Europe using the yardstick method.

Integrating data across studies with traditional microsatellite genetic markers requires careful calibration and represents an obstacle for investigation of wide-ranging species where populations require transboundary management. We used the "yardstick" method to compare results published across Europe since 2002 and new wolf (Canis lupus) genetic profiles from the Carpathian Mountains in Central Europe and the Dinaric Mountains in Southeastern Europe, with the latter as our reference population. We compared each population with Dinaric wolves, considering only shared markers (range 4-17). For each population, we calculated standard genetic diversity indices plus calibrated heterozygosity (Hec) and allelic richness (Ac). Hec and Ac in Dinaric (0.704 and 9.394) and Carpathian wolves (0.695 and 7.023) were comparable to those observed in other large and mid-sized European populations, but smaller than those of northeastern Europe. Major discrepancies in marker choices among some studies made comparisons more difficult. However, the yardstick method, including the new measures of Hec and Ac, provided a direct comparison of genetic diversity values among wolf populations and an intuitive interpretation of the results. The yardstick method thus permitted the integration of diverse sources of publicly available microsatellite data for spatiotemporal genetic monitoring of evolutionary potential.

Monitoring the effective population size of a brown bear (Ursus arctos) population using new single-sample approaches.

The effective population size (N(e) ) could be the ideal parameter for monitoring populations of conservation concern as it conveniently summarizes both the evolutionary potential of the population and its sensitivity to genetic stochasticity. However, tracing its change through time is difficult in natural populations. We applied four new methods for estimating N(e) from a single sample of genotypes to trace temporal change in N(e) for bears in the Northern Dinaric Mountains. We genotyped 510 bears using 20 microsatellite loci and determined their age. The samples were organized into cohorts with regard to the year when the animals were born and yearly samples with age categories for every year when they were alive. We used the Estimator by Parentage Assignment (EPA) to directly estimate both N(e) and generation interval for each yearly sample. For cohorts, we estimated the effective number of breeders (N(b) ) using linkage disequilibrium, sibship assignment and approximate Bayesian computation methods and extrapolated these estimates to N(e) using the generation interval. The N(e) estimate by EPA is 276 (183-350 95% CI), meeting the inbreeding-avoidance criterion of N(e) > 50 but short of the long-term minimum viable population goal of N(e) > 500. The results obtained by the other methods are highly consistent with this result, and all indicate a rapid increase in N(e) probably in the late 1990s and early 2000s. The new single-sample approaches to the estimation of N(e) provide efficient means for including N(e) in monitoring frameworks and will be of great importance for future management and conservation.

Phylogeny and phylogeography of medicinal leeches (genus Hirudo): fast dispersal and shallow genetic structure.

Medicinal leeches (Hirudo spp.) are among the best-studied invertebrates in many aspects of their biology. Yet, relatively little is known about their biogeography, ecology and evolution. Previous studies found vast ranges but suggested low genetic diversity for some species. To examine this apparent contradiction, the phylogeny and phylogeography of the widespread Hirudo verbana, Hirudo medicinalis and Hirudo orientalis were investigated in a comparative manner. Populations from across their ranges in Europe, Asia Minor, the Caucasus and Central Asia, were analyzed by various phylogenetic and population genetic approaches using both mitochondrial (COI and 12S) and nuclear DNA sequences (ITS1, 5.8S and ITS2). The populations showed surprisingly little genetic differentiation despite vast ranges. The only clear structure was observed in H. verbana. This species is subdivided into an Eastern (southern Ukraine, North Caucasus, Turkey and Uzbekistan) and a Western phylogroup (Balkans and Italy). The two phylogroups do not overlap, suggesting distinct postglacial colonization from separate refugia. Leeches supplied by commercial facilities belong to the Eastern phylogroup of H. verbana; they originate from Turkey and the Krasnodar Territory in Russia, two leading areas of leech export. H. verbana and H. medicinalis have experienced recent rapid population growth and range expansion, while isolation by distance has shaped the genetic setup of H. orientalis. The habitat of the latter is patchy and scattered about inhospitable arid and alpine areas of Central Asia and Transcaucasia. Centuries of leech collecting and transport across Europe seem not to have affected the natural distribution of genetic diversity, as the observed patterns can be explained by a combination of historical factors and present day climatic influences.

A method for measuring support for synapomorphy using character state distributions on phylogenetic trees.

Synapomorphies are fundamental to phylogenetic systematics as they offer empirical evidence of monophyletic groups. However, no method exists to directly measure synapomorphy. Here, we propose a method that quantifies synapomorphy using the pattern of character state distribution over a cladogram separately for each character and for each clade. We define a fully synapomorphic character state as one shared by all of a clade's terminal taxa and at the same time completely absent from all terminal taxa outside that clade. The extent to which this condition is met corresponds to the support for the character state being synapomorphic or, in short, support for synapomorphy. It is calculated as the probability of randomly selecting, by multi-stage sampling following the topology of the tree, two terminals from inside a clade sharing the same character state and one terminal from outside the clade bearing a different character state. The method is independent of tree inference and free of transformational assumptions, and so can be applied to any tree and used for any type of discrete character. By measuring synapomorphy, the method offers a potential tool for determining diagnostic character states for taxa on different hierarchical levels, for evaluating alternative systems of character coding, and for evaluating clade support. We show how the method differs from ancestral character state reconstruction methods and goodness-of-fit indices. We demonstrate the behaviour of our method with several hypothetical scenarios and its potential use with two real-life examples.

Niche-based mechanisms operating within extreme habitats: a case study of subterranean amphipod communities.

It has been suggested that both niche-based and neutral mechanisms are important for biological communities to evolve and persist. For communities in extreme and isolated environments such as caves, theoretical and empirical considerations (low species turnover, high stress, strong convergence owing to strong directional selection) predict neutral mechanisms and functional equivalence of species. We tested this prediction using subterranean amphipod communities from caves and interstitial groundwater. Contrary to expectations, functional morphological diversity within communities in both habitats turned out to be significantly higher than the null model of randomly assembled communities. This suggests that even the most extreme, energy-poor environments still maintain the potential for diversification via differentiation of niches.

The minnow Phoxinus lumaireul (Leuciscidae) shifts the Adriatic-Black Sea basin divide in the north-western Dinaric Karst region.

Karst landscapes are characterized by intermittent and sinking streams. The most common method used to study underground hydrological connections in karst is tracing tests. However, a more biologically oriented approach has been suggested: analysis of the genetic structure of aquatic organisms. Biological tracers can be sought among trogloxenes, that is, surface species that occasionally enter caves and groundwater. One such example is the fish genus Phoxinus, which exhibits high genetic diversity and complex phylogeography in the Balkan Peninsula. In the north-western Dinaric Karst, the complex hydrological network was digitalized in 2020. Contemporaneously, Phoxinus lumaireul populations in the Slovenian Dinaric Karst were intensively sampled and analysed for fragments of two mitochondrial genes and one nuclear gene. The derived phylogeographic structure and data on hydrological connections were compared to evaluate support for three alternative scenarios: The genetic structure (1) is a consequence of the ongoing geneflow through underground connections, (2) reflects a previous hydrological network or (3) is an outcome of anthropogenic translocations. The results suggest that the first two scenarios seem to have played a major role, while the third has not had profound effects on the genetic composition. Comparison between the genetic structure of Slovenian Dinaric Karst sampling sites and that of hydrologically isolated reference sampling sites indicated a greater genetic connectivity in the former. Moreover, the range of Adriatic (1a) and Black Sea (1c) haplotypes does not correspond to the Adriatic-Black Sea basin divide but is shifted northwards.

Toward the massive genome of Proteus anguinus-illuminating longevity, regeneration, convergent evolution, and metabolic disorders.

Deciphering the genetic code of organisms with unusual phenotypes can help answer fundamental biological questions and provide insight into mechanisms relevant to human biomedical research. The cave salamander Proteus anguinus (Urodela: Proteidae), also known as the olm, is an example of a species with unique morphological and physiological adaptations to its subterranean environment, including regenerative abilities, resistance to prolonged starvation, and a life span of more than 100 years. However, the structure and sequence of the olm genome is still largely unknown owing to its enormous size, estimated at nearly 50 gigabases. An international Proteus Genome Research Consortium has been formed to decipher the olm genome. This perspective provides the scientific and biomedical rationale for exploring the olm genome and outlines potential outcomes, challenges, and methodological approaches required to analyze and annotate the genome of this unique amphibian.

A subterranean adaptive radiation of amphipods in Europe.

Adaptive radiations are bursts of evolutionary species diversification that have contributed to much of the species diversity on Earth. An exception is modern Europe, where descendants of ancient adaptive radiations went extinct, and extant adaptive radiations are small, recent and narrowly confined. However, not all legacy of old radiations has been lost. Subterranean environments, which are dark and food-deprived, yet buffered from climate change, have preserved ancient lineages. Here we provide evidence of an entirely subterranean adaptive radiation of the amphipod genus Niphargus, counting hundreds of species. Our modelling of lineage diversification and evolution of morphological and ecological traits using a time-calibrated multilocus phylogeny suggests a major adaptive radiation, comprised of multiple subordinate adaptive radiations. Their spatio-temporal origin coincides with the uplift of carbonate massifs in South-Eastern Europe 15 million years ago. Emerging subterranean environments likely provided unoccupied, predator-free space, constituting ecological opportunity, a key trigger of adaptive radiation. This discovery sheds new light on the biodiversity of Europe.

Origins of Allium ampeloprasum horticultural groups and a molecular phylogeny of the section Allium (Allium: Alliaceae).

The subgenus Allium section Allium includes economically important species, such as garlic and leek, as well as other polyploid minor crops. Phylogenetic studies within this section, with a focus on horticultural groups within A. ampeloprasum, were performed on 31 accessions of 17 species using the nuclear ribosomal DNA internal transcribed spacer (ITS) region and the chloroplast trnL-F and trnD-T regions. The results confirmed the monophyly of section Allium. Four main clades were identified on all ITS analyses but the relationships among those and the remaining species studied within section Allium remained unresolved. Trees based on cpDNA recovered two major clades and a topology only partly congruent with that of the ITS tree. Intra-individual polymorphism of the ITS region proved useful in tracking putative parent species of polyploid taxa. The allopolyploid origin of great headed garlic (GHG), A. iranicum and A. polyanthum was confirmed. No signs of hybridization in leek or kurrat were detected but possible introgression events were identified in pearl onion and bulbous leek. Although GHG is often used as a garlic substitute, molecular analysis revealed only a distant relationship with garlic. We also clarified the previous incorrect classification of cultivated forms within A. ampeloprasum, by showing that leek, kurrat, pearl onion, and bulbous leek should be considered separately from GHG.

Lack of genetic structure in the jellyfish Pelagia noctiluca (Cnidaria: Scyphozoa: Semaeostomeae) across European seas.

The genetic structure of the holopelagic scyphozoan Pelagia noctiluca was inferred based on the study of 144 adult medusae. The areas of study were five geographic regions in two European seas (Eastern Atlantic and Mediterranean Sea). A 655-bp sequence of mitochondrial cytochrome c oxidase subunit I (COI), and a 645-bp sequence of two nuclear internal transcribed spacers (ITS1 and ITS2) were analyzed. The protein coding COI gene showed a higher level of divergence than the combined nuclear ITS fragment (haplotype diversity 0.962 vs. 0.723, nucleotide diversity 1.16% vs. 0.31%). Phylogeographic analysis on COI gene revealed two clades, the larger consisting of specimens from all sampling sites, and the smaller mostly formed of specimens from the Mediterranean Sea. Haplotype diversity was very high throughout the sampled area, and within sample diversity was higher than diversity among geographical regions. No strongly supported genetically or geographically distinct groups of P. noctiluca were found. The results - long distance dispersal, insignificant F(ST) values, lack of isolation by distance - pointed toward an admixture among Mediterranean and East Atlantic populations.

Draft genome of the European medicinal leech Hirudo medicinalis (Annelida, Clitellata, Hirudiniformes) with emphasis on anticoagulants.

The European medicinal leech has been used for medicinal purposes for millennia, and continues to be used today in modern hospital settings. Its utility is granted by the extremely potent anticoagulation factors that the leech secretes into the incision wound during feeding and, although a handful of studies have targeted certain anticoagulants, the full range of anticoagulation factors expressed by this species remains unknown. Here, we present the first draft genome of the European medicinal leech, Hirudo medicinalis, and estimate that we have sequenced between 79-94% of the full genome. Leveraging these data, we searched for anticoagulation factors across the genome of H. medicinalis. Following orthology determination through a series of BLAST searches, as well as phylogenetic analyses, we estimate that fully 15 different known anticoagulation factors are utilized by the species, and that 17 other proteins that have been linked to antihemostasis are also present in the genome. We underscore the utility of the draft genome for comparative studies of leeches and discuss our results in an evolutionary context.

The limits of cryptic diversity in groundwater: phylogeography of the cave shrimp Troglocaris anophthalmus (Crustacea: Decapoda: Atyidae).

Recent studies have revealed high local diversity and endemism in groundwaters, and showed that species with large ranges are extremely rare. One of such species is the cave shrimp Troglocaris anophthalmus from the Dinaric Karst on the western Balkan Peninsula, apparently uniform across a range of more than 500 kilometers. As such it contradicts the paradigm that subterranean organisms form localized, long-term stable populations that cannot disperse over long distances. We tested it for possible cryptic diversity and/or unexpected evolutionary processes, analysing mitochondrial (COI, 16S rRNA) and nuclear (ITS2) genes of 232 specimens from the entire range. The results of an array of phylogeographical procedures congruently suggested that the picture of a widespread, continuously distributed and homogenous T. anophthalmus was wrong. The taxon is composed of four or possibly five monophyletic, geographically defined phylogroups that meet several species delimitation criteria, two of them showing evidence of biological reproductive isolation in sympatry. COI genetic distances between phylogroups turned out to be a poor predictor, as they were much lower than the sometimes suggested crustacean threshold value of 0.16 substitutions per site. Most results confirmed the nondispersal hypothesis of subterranean fauna, but the southern Adriatic phylogroup displayed a paradoxical pattern of recent dispersal across 300 kilometers of hydrographically fragmented karst terrain. We suggest a model of migration under extreme water-level conditions, when flooded poljes could act as stepping-stones. In the north of the range (Slovenia), the results confirmed the existence of a zone of unique biogeographical conflict, where surface fauna is concordant with the current watershed, and subterranean fauna is not.

Testing the uniqueness of deep terrestrial life.

Terrestrial life typically does not occur at depths greater than a few meters. Notable exceptions are massifs of fissured rock with caves and hollow spaces reaching depths of two kilometres and more. Recent biological discoveries from extremely deep caves have been reported as sensations analogous to wondrous deep sea creatures. However, the existence of unique deep terrestrial communities is questionable when caves are understood as integral parts of a bedrock fissure network (BFN) interconnecting all parts of a massif horizontally and vertically. We tested these two opposing hypotheses - unique deep cave fauna vs. BFN - by sampling subterranean communities within the 3D matrix of a typical karst massif. There was no distinction between deep core and shallow upper zone communities. Beta diversity patterns analysed against null models of random distribution were generally congruent with the BFN hypothesis, but suggested gravity-assisted concentration of fauna in deep caves and temperature-dependent horizontal distribution. We propose that the idea of a unique deep terrestrial fauna akin to deep oceanic life is unsupported by data and unwarranted by ecological considerations. Instead, the BFN hypothesis and local ecological and structural factors sufficiently explain the distribution of subterranean terrestrial life even in the deepest karst massifs.

Composition of the cutaneous bacterial community of a cave amphibian, Proteus anguinus.

The European cave salamander Proteus anguinus is a charismatic amphibian endemic to the concealed and inaccessible subterranean waters of the Dinaric Karst. Despite its exceptional conservation importance not much is known about its ecology and interactions with the groundwater microbiome. The cutaneous microbiota of amphibians is an important driver of metabolic capabilities and immunity, and thus a key factor in their wellbeing and survival. We used high-throughput 16S rRNA gene sequencing based on seven variable regions to examine the bacteriome of the skin of five distinct evolutionary lineages of P. anguinus and in their groundwater environment. The skin bacteriomes turned out to be strongly filtered subsamples of the environmental microbial community. The resident microbiota of the analyzed individuals was dominated by five bacterial taxa. Despite an indicated functional redundancy, the cutaneous bacteriome of P. anguinus presumably provides protection against invading microbes by occupying the niche, and thus could serve as an indicator of health status. Besides conservation implications for P. anguinus, our results provide a baseline for future studies on other endangered neotenic salamanders.