Genomic evidence for the widespread presence of GH45 cellulases among soil invertebrates

Lignocellulose is a major component of vascular plant biomass. Its decomposition is crucial for the terrestrial carbon cycle. Microorganisms are considered primary decomposers, but evidence increases that some invertebrates may also decompose lignocellulose. We investigated the taxonomic distribution and evolutionary origins of GH45 hydrolases, important enzymes for the decomposition of cellulose and hemicellulose, in a collection of soil invertebrate genomes. We found that these genes are common in springtails and oribatid mites. Phylogenetic analysis revealed that cellulase genes were acquired early in the evolutionary history of these groups. Domain architectures and predicted 3D enzyme structures indicate that these cellulases are functional. Patterns of presence and absence of these genes across different lineages prompt further investigation into their evolutionary and ecological benefits. The ubiquity of cellulase genes suggests that soil invertebrates may play a role in lignocellulose decomposition, independently or in synergy with microorganisms. Understanding the ecological and evolutionary implications might be crucial for understanding soil food webs and the carbon cycle.

The MetaInvert soil invertebrate genome resource provides insights into below-ground biodiversity and evolution

Soil invertebrates are among the least understood metazoans on Earth. Thus far, the lack of taxonomically broad and dense genomic resources has made it hard to thoroughly investigate their evolution and ecology. With MetaInvert we provide draft genome assemblies for 232 soil invertebrate species, representing 14 common groups and 94 families. We show that this data substantially extends the taxonomic scope of DNA- or RNA-based taxonomic identification. Moreover, we confirm that theories of genome evolution cannot be generalised across evolutionarily distinct invertebrate groups. The soil invertebrate genomes presented here will support the management of soil biodiversity through molecular monitoring of community composition and function, and the discovery of evolutionary adaptations to the challenges of soil conditions.

Metatranscriptomics reveals contrasting effects of elevation on the activity of bacteria and bacterial viruses in soil

Soil microbial diversity affects ecosystem functioning and global biogeochemical cycles. Soil bacterial communities catalyse a diversity of biogeochemical reactions and have thus sparked considerable scientific interest. One driver of bacterial community dynamics in natural ecosystems has so far been largely neglected: the predator-prey interactions between bacterial viruses (bacteriophages) and bacteria. To generate ground level knowledge on environmental drivers of these particular predator-prey dynamics, we propose an activity-based ecological framework to simultaneous capture community dynamics of bacteria and bacteriophages in soils. An ecological framework and specifically the analyses of community dynamics across latitudinal and elevational gradients have been widely used in ecology to understand community-wide responses of innumerable taxa to environmental change, in particular to climate. Here, we tested the hypothesis that the activity of bacteria and bacteriophages codeclines across an elevational gradient. We used metatranscriptomics to investigate bacterial and bacteriophage activity patterns at five sites across 400 elevational metres in the Swiss Alps in 2015 and 2017. We found that metabolic activity (transcription levels) of bacteria declined significantly with increasing elevation, but activity of bacteriophages did not. We showed that bacteriophages are consistently active in soil along the entire gradient, making bacteriophage activity patterns divergent from that of their putative bacterial prey. Future efforts will be necessary to link the environment-activity relationship to predator-prey dynamics, and to understand the magnitude of viral contributions to carbon, nitrogen and phosphorus cycling when infection causes bacterial cell death, a process that may represent an overlooked component of soil biogeochemical cycles.

Dissent in the sediment? Lake sediments as archives of short- and long-range impact of anthropogenic activities in northeastern Germany

The suitability of lake sediment cores to reconstruct past inputs, regional pollution, and usage patterns of pesticides has been shown previously. Until now, no such data exist for lakes in eastern Germany. Therefore, 10 sediment cores (length 1 m) of 10 lakes in eastern Germany, the territory of the former German Democratic Republic (GDR), were collected and cut into 5-10-mm layers. In each layer, concentrations of trace elements (TEs) As, Cd, Cr, Cu, Ni, Pb, S, and Zn, as well as of organochlorine pesticides (OCPs), i.e., dichlorodiphenyltrichloroethane (DDT) and hexachlorocyclohexane (HCH), were analyzed. A miniaturized solid-liquid extraction technique in conjunction with headspace solid-phase microextraction (HS-SPME) and gas chromatography-mass spectrometry (GC-MS) was used for the latter. The progression of TE concentrations over time is uniform. It follows a trans-regional pattern and is indicative of activity and policy making in West Germany before 1990 instead of those in the GDR. Of OCPs, only transformation products of DDT were found. Congener ratios indicate a mainly aerial input. In the lakes' profiles, several regional features and responses to national policies and measures are visible. Dichlorodiphenyldichloroethane (DDD) concentrations reflect the history of DDT use in the GDR. Lake sediments proved to be suitable to archive short- and long-range impacts of anthropogenic activity. Our data can be used to complement and validate other forms of environmental pollution long-term monitoring and to check for the efficiency of pollution countermeasures in the past.

RNA interference to combat the Asian tiger mosquito in Europe: A pathway from design of an innovative vector control tool to its application

The Asian tiger mosquito Aedes albopictus is currently spreading across Europe, facilitated by climate change and global transportation. It is a vector of arboviruses causing human diseases such as chikungunya, dengue hemorrhagic fever and Zika fever. For the majority of these diseases, no vaccines or therapeutics are available. Options for the control of Ae. albopictus are limited by European regulations introduced to protect biodiversity by restricting or phasing out the use of pesticides, genetically modified organisms (GMOs) or products of genome editing. Alternative solutions are thus urgently needed to avoid a future scenario in which Europe faces a choice between prioritizing human health or biodiversity when it comes to Aedes-vectored pathogens. To ensure regulatory compliance and public acceptance, these solutions should preferably not be based on chemicals or GMOs and must be cost-efficient and specific. The present review aims to synthesize available evidence on RNAi-based mosquito vector control and its potential for application in the European Union. The recent literature has identified some potential target sites in Ae. albopictus and formulations for delivery. However, we found little information concerning non-target effects on the environment or human health, on social aspects, regulatory frameworks, or on management perspectives. We propose optimal designs for RNAi-based vector control tools against Ae. albopictus (target product profiles), discuss their efficacy and reflect on potential risks to environmental health and the importance of societal aspects. The roadmap from design to application will provide readers with a comprehensive perspective on the application of emerging RNAi-based vector control tools for the suppression of Ae. albopictus populations with special focus on Europe.

How genomics can help biodiversity conservation

The availability of public genomic resources can greatly assist biodiversity assessment, conservation, and restoration efforts by providing evidence for scientifically informed management decisions. Here we survey the main approaches and applications in biodiversity and conservation genomics, considering practical factors, such as cost, time, prerequisite skills, and current shortcomings of applications. Most approaches perform best in combination with reference genomes from the target species or closely related species. We review case studies to illustrate how reference genomes can facilitate biodiversity research and conservation across the tree of life. We conclude that the time is ripe to view reference genomes as fundamental resources and to integrate their use as a best practice in conservation genomics.

Host-pathogen coevolution drives innate immune response to Aphanomyces astaci infection in freshwater crayfish: transcriptomic evidence

BACKGROUND

For over a century, scientists have studied host-pathogen interactions between the crayfish plague disease agent Aphanomyces astaci and freshwater crayfish. It has been hypothesised that North American crayfish hosts are disease-resistant due to the long-lasting coevolution with the pathogen. Similarly, the increasing number of latent infections reported in the historically sensitive European crayfish hosts seems to indicate that similar coevolutionary processes are occurring between European crayfish and A. astaci. Our current understanding of these host-pathogen interactions is largely focused on the innate immunity processes in the crayfish haemolymph and cuticle, but the molecular basis of the observed disease-resistance and susceptibility remain unclear. To understand how coevolution is shaping the host's molecular response to the pathogen, susceptible native European noble crayfish and invasive disease-resistant marbled crayfish were challenged with two A. astaci strains of different origin: a haplogroup A strain (introduced to Europe at least 50 years ago, low virulence) and a haplogroup B strain (signal crayfish in lake Tahoe, USA, high virulence). Here, we compare the gene expression profiles of the hepatopancreas, an integrated organ of crayfish immunity and metabolism.

RESULTS

We characterised several novel innate immune-related gene groups in both crayfish species. Across all challenge groups, we detected 412 differentially expressed genes (DEGs) in the noble crayfish, and 257 DEGs in the marbled crayfish. In the noble crayfish, a clear immune response was detected to the haplogroup B strain, but not to the haplogroup A strain. In contrast, in the marbled crayfish we detected an immune response to the haplogroup A strain, but not to the haplogroup B strain.

CONCLUSIONS

We highlight the hepatopancreas as an important hub for the synthesis of immune molecules in the response to A. astaci. A clear distinction between the innate immune response in the marbled crayfish and the noble crayfish is the capability of the marbled crayfish to mobilise a higher variety of innate immune response effectors. With this study we outline that the type and strength of the host immune response to the pathogen is strongly influenced by the coevolutionary history of the crayfish with specific A. astaci strains.

Dataset of the de novo assembly and annotation of the marbled crayfish and the noble crayfish hepatopancreas transcriptomes

OBJECTIVES

Crayfish plague disease, caused by the oomycete pathogen Aphanomyces astaci represents one of the greatest risks for the biodiversity of the freshwater crayfish. This data article covers the de novo transcriptome assembly and annotation data of the noble crayfish and the marbled crayfish challenged with Ap. astaci. Following the controlled infection experiment (Francesconi et al. in Front Ecol Evol, 2021, https://doi.org/10.3389/fevo.2021.647037 ), we conducted a differential gene expression analysis described in (Boštjančić et al. in BMC Genom, 2022, https://doi.org/10.1186/s12864-022-08571-z ) DATA DESCRIPTION: In total, 25 noble crayfish and 30 marbled crayfish were selected. Hepatopancreas tissue was isolated, followed by RNA sequencing using the Illumina NovaSeq 6000 platform. Raw data was checked for quality with FastQC, adapter and quality trimming were conducted using Trimmomatic followed by de novo assembly with Trinity. Assembly quality was assessed with BUSCO, at 93.30% and 93.98% completeness for the noble crayfish and the marbled crayfish, respectively. Transcripts were annotated using the Dammit! pipeline and assigned to KEGG pathways. Respective transcriptome and raw datasets may be reused as the reference transcriptome assemblies for future expression studies.

Shotgun metagenomics of soil invertebrate communities reflects taxonomy, biomass, and reference genome properties

Metagenomics - shotgun sequencing of all DNA fragments from a community DNA extract - is routinely used to describe the composition, structure, and function of microorganism communities. Advances in DNA sequencing and the availability of genome databases increasingly allow the use of shotgun metagenomics on eukaryotic communities. Metagenomics offers major advances in the recovery of biomass relationships in a sample, in comparison to taxonomic marker gene-based approaches (metabarcoding). However, little is known about the factors which influence metagenomics data from eukaryotic communities, such as differences among organism groups, the properties of reference genomes, and genome assemblies.We evaluated how shotgun metagenomics records composition and biomass in artificial soil invertebrate communities at different sequencing efforts. We generated mock communities of controlled biomass ratios from 28 species from all major soil mesofauna groups: mites, springtails, nematodes, tardigrades, and potworms. We shotgun sequenced these communities and taxonomically assigned them with a database of over 270 soil invertebrate genomes.We recovered over 95% of the species, and observed relatively high false-positive detection rates. We found strong differences in reads assigned to different taxa, with some groups (e.g., springtails) consistently attracting more hits than others (e.g., enchytraeids). Original biomass could be predicted from read counts after considering these taxon-specific differences. Species with larger genomes, and with more complete assemblies, consistently attracted more reads than species with smaller genomes. The GC content of the genome assemblies had no effect on the biomass-read relationships. Results were similar among different sequencing efforts.The results show considerable differences in taxon recovery and taxon specificity of biomass recovery from metagenomic sequence data. The properties of reference genomes and genome assemblies also influence biomass recovery, and they should be considered in metagenomic studies of eukaryotes. We show that low- and high-sequencing efforts yield similar results, suggesting high cost-efficiency of metagenomics for eukaryotic communities. We provide a brief roadmap for investigating factors which influence metagenomics-based eukaryotic community reconstructions. Understanding these factors is timely as accessibility of DNA sequencing and momentum for reference genomes projects show a future where the taxonomic assignment of DNA from any community sample becomes a reality.

The era of reference genomes in conservation genomics

Progress in genome sequencing now enables the large-scale generation of reference genomes. Various international initiatives aim to generate reference genomes representing global biodiversity. These genomes provide unique insights into genomic diversity and architecture, thereby enabling comprehensive analyses of population and functional genomics, and are expected to revolutionize conservation genomics.

Genomic basis for drought resistance in European beech forests threatened by climate change

In the course of global climate change, Central Europe is experiencing more frequent and prolonged periods of drought. The drought years 2018 and 2019 affected European beeches (Fagus sylvatica L.) differently: even in the same stand, drought-damaged trees neighboured healthy trees, suggesting that the genotype rather than the environment was responsible for this conspicuous pattern. We used this natural experiment to study the genomic basis of drought resistance with Pool-GWAS. Contrasting the extreme phenotypes identified 106 significantly associated single-nucleotide polymorphisms (SNPs) throughout the genome. Most annotated genes with associated SNPs (>70%) were previously implicated in the drought reaction of plants. Non-synonymous substitutions led either to a functional amino acid exchange or premature termination. An SNP assay with 70 loci allowed predicting drought phenotype in 98.6% of a validation sample of 92 trees. Drought resistance in European beech is a moderately polygenic trait that should respond well to natural selection, selective management, and breeding.

Climate change is having a serious impact on many ecosystems. In the summer of 2018 and 2019, around two thirds of European beech trees were damaged or killed by extreme drought. It is critical to keep these beech woods healthy, as they are central to the survival of over 6,000 other species of animals and plants.

The level of damage caused by the drought varied between forests. However, not all the trees in each forest responded in the same way, with severely damaged trees often sitting next to fully healthy ones. This suggests that the genetic make-up of each tree determines how well it can adapt to drought rather than its local environment.

To investigate this further, Pfenninger et al. studied the genome of over 400 European beech trees from the Hesse region in Germany. The samples came from pairs of neighbouring trees that had responded differently to the droughts. The analysis found more than 80 parts of the genome that differed between healthy and damaged trees.

Pfenninger et al. then used this information to create a genetic test which can quickly and inexpensively predict how well an individual beech tree might survive in a drought. Applying this test to another 92 trees revealed that it can reliably detect which ones were healthy and which ones were damaged.

Beech forests are typically managed by private owners, agencies or breeders that could use this genetic test to select and reproduce trees that are better adapted to drought. The goal now is to develop the test so that it can be used more widely to manage European beech trees and potentially other species.

Two high-quality de novo genomes from single ethanol-preserved specimens of tiny metazoans (Collembola)

BACKGROUND

Genome sequencing of all known eukaryotes on Earth promises unprecedented advances in biological sciences and in biodiversity-related applied fields such as environmental management and natural product research. Advances in long-read DNA sequencing make it feasible to generate high-quality genomes for many non-genetic model species. However, long-read sequencing today relies on sizable quantities of high-quality, high molecular weight DNA, which is mostly obtained from fresh tissues. This is a challenge for biodiversity genomics of most metazoan species, which are tiny and need to be preserved immediately after collection. Here we present de novo genomes of 2 species of submillimeter Collembola. For each, we prepared the sequencing library from high molecular weight DNA extracted from a single specimen and using a novel ultra-low input protocol from Pacific Biosciences. This protocol requires a DNA input of only 5 ng, permitted by a whole-genome amplification step.

RESULTS

The 2 assembled genomes have N50 values >5.5 and 8.5 Mb, respectively, and both contain ∼96% of BUSCO genes. Thus, they are highly contiguous and complete. The genomes are supported by an integrative taxonomy approach including placement in a genome-based phylogeny of Collembola and designation of a neotype for 1 of the species. Higher heterozygosity values are recorded in the more mobile species. Both species are devoid of the biosynthetic pathway for β-lactam antibiotics known in several Collembola, confirming the tight correlation of antibiotic synthesis with the species way of life.

CONCLUSIONS

It is now possible to generate high-quality genomes from single specimens of minute, field-preserved metazoans, exceeding the minimum contig N50 (1 Mb) required by the Earth BioGenome Project.

Genome-level insights into the operation of an on-site biological wastewater treatment unit reveal the importance of storage time

On-site wastewater treatment systems are gaining popularity in areas where centralized wastewater treatment is not available. In the current case study a domestic activated sludge system was investigated, where treated effluent was stored in a short-term (1 week turn-over time) and a long-term (over 2-3 months) storage tank and was then used for irrigation. This design provided a unique opportunity to assess the chemical and microbial changes of the effluent upon storage. Long-term storage greatly improved both the chemical quality and the degradation efficiency of most organic micropollutants examined, including petroleum hydrocarbons and the pesticide diethyltoluamide. Taxonomic profile of the core microbiome of the effluent was also influenced upon storage. Relative abundance values of the members of Azoarcus and Thauera genera, which are important in degrading polycyclic aromatic hydrocarbons compounds, clearly indicated the biodegrading activity of these microbes across samples. The abundance of xenobiotics degradation functions correlated with the observed organic micropollutant degradation values indicating efficient microbial decomposition of these contaminants. Functions related to infectious diseases also had the highest abundance in the short-term storage tank corresponding well with the relative abundance of indicator organisms and implying to the significance of storage time in the elimination of pathogens. Based on these results, small, on-site wastewater treatment systems could benefit from long-term storage of wastewater effluent.

Larval morphology of the Western Balkans endemic caddisflies Drusus krusniki Malicky 1981, D. vernonensis Malicky 1989, and D. vespertinus Marinković 1976 (Trichoptera, Limnephilidae, Drusinae)

Drusinae (Trichoptera, Limnephilidae) are highland caddisflies inhabiting high-gradient, turbulent running water and spring habitats. They are disjunctly distributed over the Eurasian mountain ranges, and the majority of species is endemic to particular mountain areas. The most diverse of three main groups of the Drusinae, the grazer clade, consists of species in which larvae feed on epiltihic biofilm and algae. In this paper we describe three previously unknown grazer-clade Drusinae larvae: Drusus krusniki Malicky 1981 (endemic to the Dinaric western Balkans), D. vernonensis Malicky 1989 (endemic to the Hellenic western Balkans), and D. vespertinus Marinković 1976 (endemic to the Dinaric western Balkans). The larvae of these species have toothless mandibles typical of the Drusinae grazer clade. Larvae and adults were unambiguously associated using molecular genetic data, i.e., the mitochondrial cytochrome oxidase I gene fragment (mtCOI3-P). Morphological characteristics of the larvae are described and the diagnostic features enabling species-level identification are illustrated. We further discuss the ecology and distribution of three Western Balkan endemic species.

Proper experimental design requires randomization/balancing of molecular ecology experiments

Properly designed (randomized and/or balanced) experiments are standard in ecological research. Molecular methods are increasingly used in ecology, but studies generally do not report the detailed design of sample processing in the laboratory. This may strongly influence the interpretability of results if the laboratory procedures do not account for the confounding effects of unexpected laboratory events. We demonstrate this with a simple experiment where unexpected differences in laboratory processing of samples would have biased results if randomization in DNA extraction and PCR steps do not provide safeguards. We emphasize the need for proper experimental design and reporting of the laboratory phase of molecular ecology research to ensure the reliability and interpretability of results.

Season-Long Experimental Drought Alters Fungal Community Composition but Not Diversity in a Grassland Soil

Using terrestrial model ecosystems (TMEs), we investigated how reduced moisture conditions impact soil fungal communities from a temperate grassland over the course of an entire season. Starting at about 65% of the soil's maximum water holding capacity (WHC_max), TME soils were adjusted to three moisture levels for 15 weeks: 70% WHC_max, approximating starting conditions, 50% WHC_max, and 30% WHC_max, representing reduced moisture conditions. Diversity and abundances of soil fungi at the start and at the end of the experiment were characterized using Illumina meta-barcoding. Community diversity at the end of the experiment did not differ between experimental moisture levels and was comparable to diversity measures from the field. However, fungal communities did change compositionally in both abundances and presence/absence of species. Analyzing class-level and individual contributions of fungi to these changes revealed that only a minor portion reacted significantly, indicating that most compositional change was likely driven by many consistent small-scale shifts in presence/absences or abundances. Together, our results show that prolonged reduction in soil moisture conditions will trigger compositional changes in soil fungal communities but not necessarily change overall diversity. We highlight the cumulative contribution of minor but consistent changes among community members, as opposed to significant responses of individual species. We also detected a strong general experimental effect on soil fungi that are moved from the field to experimental TMEs, suggesting the importance of acclimatization effects in these communities under laboratory conditions.

Integrative taxonomy by molecular species delimitation: multi-locus data corroborate a new species of Balkan Drusinae micro-endemics

BACKGROUND

Taxonomy offers precise species identification and delimitation and thus provides basic information for biological research, e.g. through assessment of species richness. The importance of molecular taxonomy, i.e., the identification and delimitation of taxa based on molecular markers, has increased in the past decade. Recently developed exploratory tools now allow estimating species-level diversity in multi-locus molecular datasets.

RESULTS

Here we use molecular species delimitation tools that either quantify differences in intra- and interspecific variability of loci, or divergence times within and between species, or perform coalescent species tree inference to estimate species-level entities in molecular genetic datasets. We benchmark results from these methods against 14 morphologically readily differentiable species of a well-defined subgroup of the diverse Drusinae subfamily (Trichoptera, Limnephilidae). Using a 3798 bp (6 loci) molecular data set we aim to corroborate a geographically isolated new species by integrating comparative morphological studies and molecular taxonomy.

CONCLUSIONS

Our results indicate that only multi-locus species delimitation provides taxonomically relevant information. The data further corroborate the new species Drusus zivici sp. nov. We provide differential diagnostic characters and describe the male, female and larva of this new species and discuss diversity patterns of Drusinae in the Balkans. We further discuss potential and significance of molecular species delimitation. Finally we argue that enhancing collaborative integrative taxonomy will accelerate assessment of global diversity and completion of reference libraries for applied fields, e.g., conservation and biomonitoring.

Support for the evolutionary speed hypothesis from intraspecific population genetic data in the non-biting midge Chironomus riparius

The evolutionary speed hypothesis (ESH) proposes a causal mechanism for the latitudinal diversity gradient. The central idea of the ESH is that warmer temperatures lead to shorter generation times and increased mutation rates. On an absolute time scale, both should lead to an acceleration of selection and drift. Based on the ESH, we developed predictions regarding the distribution of intraspecific genetic diversity: populations of ectothermic species with more generations per year owing to warmer ambient temperatures should be more differentiated from each other, accumulate more mutations and show evidence for increased mutation rates compared with populations in colder regions. We used the multivoltine insect species Chironomus riparius to test these predictions with cytochrome oxidase I (COI) sequence data and found that populations from warmer regions are indeed significantly more differentiated and have significantly more derived haplotypes than populations from colder regions. We also found a significant correlation of the annual mean temperature with the population mutation parameter θ that serves as a proxy for the per generation mutation rate under certain assumptions. This pattern could be corroborated with two nuclear loci. Overall, our results support the ESH and indicate that the thermal regime experienced may be crucially driving the evolution of ectotherms and may thus ultimately govern their speciation rate.

New records for the Kosovo caddisfly fauna with the description of a new species, Drusus dardanicus sp. nov. (Trichoptera: Limnephilidae)

The Balkan Peninsula is one of the most important European hotspots of freshwater biodiversity. The region is, however, to a large extent insufficiently investigated. Here we present data on distribution of caddisflies in one particularly understudied area, the Republic of Kosovo. Our data include the first records of Adicella altandroconia Botosaneanu & Novak and Halesus tessellatus (Rambur) for the Kosovo caddisfly fauna, and a new locality for the recently described Ecclisopteryx keroveci Previšić, Graf, & Vitecek. Further, we describe the new caddisfly species Drusus dardanicus sp. nov. from the Kopaonik Mountains. The new species belongs to the D. discophorus Species Group and differs morphologically from its most similar congeners (D. discophorus Radovanović, D. balcanicus Kumanski, and D. bureschi Kumanski) mainly in exhibiting (1) subtrianglar superior appendages; (2) a narrow, dorsal spinate area of tergite VIII; and (3) evenly rounded tips of intermediate appendages in caudal view. In phylogenetic analysis, D. dardanicus sp. nov. is well delineated and recovered as a sister taxon to D. osogovicus Kumanski, a species recorded from Bulgaria. The recent discovery of a new species and other rare or microendemic species presents important contributions to the knowledge on the rich freshwater biodiversity in Kosovo. These species face increasing anthropogenic pressure and threats to their conservation.

Description of two new filtering carnivore Drusus species (Limnephilidae, Drusinae) from the Western Balkans

Two new species of the genus Drusus (Trichoptera, Limnephilidae, Drusinae) from the Western Balkans are described. Additionally, observations on the biodiversity and threats to the region's endemic aquatic fauna are discussed. Drususkrpachi sp. n. is a micro-endemic of the Korab Mountains, Macedonia, and Drususmalickyi sp. n. is a micro-endemic of the Prokletije Mountains, Albania. Both new species are most similar to Drususmacedonicus but differ from the latter in the shape of segment IX, the shape of the tips of the intermediate appendages in lateral view, the shape of the inferior appendages, and the form and shape of the parameres. In addition, males of the European species of filtering carnivore Drusinae are diagnosed and illustrated, including Cryptothrixnebulicola McLachlan, Drususchrysotus Rambur, Drususdiscolor Rambur, Drususmacedonicus Schmid, Drususmeridionalis Kumanski, Drususmuelleri McLachlan, Drususromanicus Murgoci and Botosaneanu, and Drusussiveci Malicky. These additions to the Western Balkan fauna demonstrate the significance of this region for European biodiversity and further highlight the importance of faunistic studies in Europe.

The larvae of Drusus franzressli Malicky 1974 and Drusus spelaeus (Ulmer 1920) (Trichoptera: Limnephilidae: Drusinae) with notes on ecology and zoogeography

Water quality monitoring is greatly dependent on identification tools for aquatic and semi-aquatic insects. Species-level identification improves resolution and precision of water quality assessment and requires comprehensive keys. With the aim of increasing the suitability of Drusinae for such applications, this paper gives a description of the hitherto unknown larvae of Drusus franzressli Malicky 1974 and Drusus spelaeus (Ulmer 1920). Information on the morphology of the larvae is given and the most important diagnostic features are illustrated. In the context of already available keys, the larvae of D. franzressli and D. spelaeus key together with Metanoea flavipennis (Pictet 1834), M. rhaetica Schmid 1956, D. improvisus McLachlan 1884, D. nigrescens Meyer-Dür 1875 and Ecclisopteryx malickyi Moretti 1991. These species are easily separated by differences in larval morphology (dorsal outline and sculpturing of pronotum, presence/absence of lateral gills at 2nd and 3rd abdominal segments, start of lateral fringe) and their distribution ranges. Drusus franzressli is endemic to the Hellenic western Balkans whereas D. spelaeus is endemic to the western Alps (Grenoble area). In addition, ecological characteristics are briefly discussed.

Description of a new species of Wormaldia from Sardinia and a new Drusus species from the Western Balkans (Trichoptera, Philopotamidae, Limnephilidae)

New species are described in the genera Wormaldia (Trichoptera, Philopotamidae) and Drusus (Trichoptera, Limnephilidae, Drusinae). Additionally, the larva of the new species Drususcrenophylax sp. n. is described, and a key provided to larval Drusus species of the bosnicus-group, in which the new species belongs. Observations on the threats to regional freshwater biodiversity and caddisfly endemism are discussed. The new species Wormaldiasarda sp. n. is an endemic of the Tyrrhenian island of Sardinia and differs most conspicuously from its congeners in the shape of segment X, which is trilobate in lateral view. The new species Drususcrenophylax sp. n. is a micro-endemic of the Western Balkans, and increases the endemism rate of Balkan Drusinae to 79% of 39 species. Compared to other Western Balkan Drusus, males of the new species are morphologically most similar to Drususdiscophorus Radovanovic and Drususvernonensis Malicky, but differ in the shape of superior and intermediate appendages. The females of Drususcrenophylax sp. n. are most similar to those of Drususvernonensis, but differ distinctly in the outline of segment X. Larvae of Drususcrenophylax sp. n. exhibit toothless mandibles, indicating a scraping grazing-feeding ecology.

Cryptic diversity of caddisflies in the Balkans: the curious case of Ecclisopteryx species (Trichoptera: Limnephilidae)

Adults and larvae of two new cryptic, endemic caddisflies, Ecclisopteryx keroveci sp.n. and Ecclisopteryx ivkae sp.n., are described and illustrated from the Western Balkans. Phylogenetic analysis (Bayesian MCMCMC) and association of different life history stages in both cryptic species were achieved through comparison of morphological characters and mitochondrial (mtCOI and mtLSU) and nuclear (nuWG) gene sequence data. The new species form a sister clade to the widely distributed E. dalecarlica and E. guttulata, with which they were formerly misidentified. Adults differ from each other and other species in the genus by the uniquely shaped inferior appendages in males and segment X in females. The larvae differ from each other and their congeners in the shape of the pronotum, and presence and constitution of additional spines on the parietalia. Larvae of both species are grazers and prefer stony substrate. Ecclisopteryx keroveci sp.n. has a wide distribution in the Western Balkans, while E. ivkae sp.n. is endemic to Dalmatia. Our findings demonstrate the significance of the Western Balkans as a freshwater biodiversity hotspot, and accentuate the importance of research focused on freshwater biodiversity and biogeography in southern Europe.

Misperceived invasion: the Lusitanian slug (Arion lusitanicus auct. non-Mabille or Arion vulgaris Moquin-Tandon 1855) is native to Central Europe

The Lusitanian slug, presumed to be native to south-west Europe, was ranked among the 100 worst invading species in Central Europe. However, from the very beginning of its recognition in the presumed invasion area, there was little evidence that the species was actually anthropogenically introduced. We investigated the invasive status of the species by comparing specific predictions on the population genetic structure in the invasion area with the pattern actually found. In a DNA-taxonomy approach, the species could not be found in its presumed native range. Using statistical phylogeographic techniques on a mitochondrial (COI) and nuclear (ZF) marker and species distribution modelling, we could show that the species is with very high probability not an invasor, but native to Central Europe. The study underlines the value of statistical phylogeography in rigorously testing hypotheses on the dynamics of biological invasions.

An Illumina metabarcoding pipeline for fungi

High-throughput metabarcoding studies on fungi and other eukaryotic microorganisms are rapidly becoming more frequent and more complex, requiring researchers to handle ever increasing amounts of raw sequence data. Here, we provide a flexible pipeline for pruning and analyzing fungal barcode (ITS rDNA) data generated as paired-end reads on Illumina MiSeq sequencers. The pipeline presented includes specific steps fine-tuned for ITS, that are mostly missing from pipelines developed for prokaryotes. It (1) employs state of the art programs and follows best practices in fungal high-throughput metabarcoding; (2) consists of modules and scripts easily modifiable by the user to ensure maximum flexibility with regard to specific needs of a project or future methodological developments; and (3) is straightforward to use, also in classroom settings. We provide detailed descriptions and revision techniques for each step, thus giving the user maximum control over data treatment and avoiding a black-box approach. Employing this pipeline will improve and speed up the tedious and error-prone process of cleaning fungal Illumina metabarcoding data.

Meta-analysis of deep-sequenced fungal communities indicates limited taxon sharing between studies and the presence of biogeographic patterns

High-throughput amplicon sequencing gives new insights into fungal community ecology. Massively generated molecular data lead to the discovery of vast fungal diversity. However, it is unclear to what extent operational taxonomic units (OTUs) overlap among independent studies, because no comparative studies exist. We compared fungal diversity based on the internal transcribed spacer (ITS1) region among 10 published studies. Starting from the raw 454 pyrosequencing data, we used a uniform pipeline to prune the reads. We investigated fungal richness and taxonomic composition among phyllosphere and soil fungal communities, as well as biogeographic signals in the data. We did not find globally distributed OTUs, even when comparing fungal communities from similar habitats (phyllosphere or soil). This suggests that high local fungal diversity scales up to high global diversity. The most OTU-rich classes in the phyllosphere were Dothideomycetes (21%) and Sordariomycetes (14%), and in the soil were Sordariomycetes (13%) and Agaricomycetes (12%). The richness estimates suggest the presence of undiscovered fungal diversity even in deeply sequenced study systems. The small number of OTUs shared among studies indicates that globally distributed taxa and habitat generalists may be rare. Latitudinal diversity decline and distance decay relationships suggest the presence of biogeographic patterns similar to those in plants and animals.

Chaetopteryx bucari sp. n., a new species from the Chaetopteryx rugulosa group from Croatia (Insecta, Trichoptera, Limnephilidae) with molecular, taxonomic and ecological notes on the group

We describe a new autumnal caddisfly species Chaetopteryx bucari sp. n. from 8 localities in the Banovina region of Croatia. We also present molecular, taxonomic and ecological notes (emergence, sex ratio and seasonal dynamics) on the new species and discuss the distribution of Chaetopteryx species in general and the Chaetopteryx rugulosa group in particular. Based on Bayesian phylogenetic analysis Chaetopteryx rugulosa schmidi was separated from the clade containing the other subspecies of Chaetopteryx rugulosa. Thus the subspecies Chaetopteryx rugulosa schmidi is here raised to species level, Chaetopteryx schmidi, as it was described originally. We further present distribution data on rare species in the genus Chaetopteryx in Croatia.

The larva of Drusus vinconi Sipahiler, 1992 (Trichoptera, Limnephilidae, Drusinae)

This paper describes the previously unknown larva of Drusus vinconi Sipahiler, 1992. Information on the morphology of the 5th larval instar is given, and the most important diagnostic features are illustrated. In the context of existing identification keys the larva of Drusus vinconi keys together with Drusus annulatus (Stephens, 1837), Drusus biguttatus (Pictet, 1834), Drusus ingridae Sipahiler, 1993, Hadimina torosensis Sipahiler, 2002 and Leptodrusus budtzi (Ulmer, 1913). These species differ in the contours of the pronotum in lateral view, the presence/absence of the pronotal transverse groove, the shape of the median notch of the pronotum (in anterior view), pronotal sculpturing, presence/absence of the lateral carina of the head capsule, the number of proximo-dorsal setae on the mid-and hind femora, where the lateral fringe starts on the abdomen, and in geographic distribution. With respect to zoogeography, Drusus vinconi is a (micro-)endemic of the Western Pyrenees. The species prefers stony substratum in springs and springbrooks of the montane and subalpine region (Graf et al. 2008; Sipahiler 1992, 1993). As a grazer, the larvae of Drusus vinconi feed on biofilm and epilithic algae.

Methodological framework for projecting the potential loss of intraspecific genetic diversity due to global climate change

BACKGROUND

While research on the impact of global climate change (GCC) on ecosystems and species is flourishing, a fundamental component of biodiversity - molecular variation - has not yet received its due attention in such studies. Here we present a methodological framework for projecting the loss of intraspecific genetic diversity due to GCC.

METHODS

The framework consists of multiple steps that combines 1) hierarchical genetic clustering methods to define comparable units of inference, 2) species accumulation curves (SAC) to infer sampling completeness, and 3) species distribution modelling (SDM) to project the genetic diversity loss under GCC. We suggest procedures for existing data sets as well as specifically designed studies. We illustrate the approach with two worked examples from a land snail (Trochulus villosus) and a caddisfly (Smicridea (S.) mucronata).

RESULTS

Sampling completeness was diagnosed on the third coarsest haplotype clade level for T. villosus and the second coarsest for S. mucronata. For both species, a substantial species range loss was projected under the chosen climate scenario. However, despite substantial differences in data set quality concerning spatial sampling and sampling depth, no loss of haplotype clades due to GCC was predicted for either species.

CONCLUSIONS

The suggested approach presents a feasible method to tap the rich resources of existing phylogeographic data sets and guide the design and analysis of studies explicitly designed to estimate the impact of GCC on a currently still neglected level of biodiversity.

Population dynamics and genetic changes of Picea abies in the South Carpathians revealed by pollen and ancient DNA analyses

Background

Studies on allele length polymorphism designate several glacial refugia for Norway spruce (Picea abies) in the South Carpathian Mountains, but infer only limited expansion from these refugia after the last glaciation. To better understand the genetic dynamics of a South Carpathian spruce lineage, we compared ancient DNA from 10,700 and 11,000-year-old spruce pollen and macrofossils retrieved from Holocene lake sediment in the Retezat Mountains with DNA extracted from extant material from the same site. We used eight primer pairs that amplified short and variable regions of the spruce cpDNA. In addition, from the same lake sediment we obtained a 15,000-years-long pollen accumulation rate (PAR) record for spruce that helped us to infer changes in population size at this site.

Results

We obtained successful amplifications for Norway spruce from 17 out of 462 pollen grains tested, while the macrofossil material provided 22 DNA sequences. Two fossil sequences were found to be unique to the ancient material. Population genetic statistics showed higher genetic diversity in the ancient individuals compared to the extant ones. Similarly, statistically significant Ks and Kst values showed a considerable level of differentiation between extant and ancient populations at the same loci.

Lateglacial and Holocene PAR values suggested that population size of the ancient population was small, in the range of 1/10 or 1/5 of the extant population. PAR analysis also detected two periods of rapid population growths (from ca. 11,100 and 3900 calibrated years before present (cal yr BP)) and three bottlenecks (around 9180, 7200 and 2200 cal yr BP), likely triggered by climatic change and human impact.

Conclusion

Our results suggest that the paternal lineages observed today in the Retezat Mountains persisted at this site at least since the early Holocene. Combination of the results from the genetic and the PAR analyses furthermore suggests that the higher level of genetic variation found in the ancient populations and the loss of ancient allele types detected in the extant individuals were likely due to the repeated bottlenecks during the Holocene; however our limited sample size did not allow us to exclude sampling effect.

This study demonstrates how past population size changes inferred from PAR records can be efficiently used in combination with ancient DNA studies. The joint application of palaeoecological and population genetics analyses proved to be a powerful tool to understand the influence of past population demographic changes on the haplotype diversity and genetic composition of forest tree species.