Environmental DNA, hydrochemistry and stable water isotopes as integrative tracers of urban ecohydrology

Urbanization and the persistent environmental changes present a major challenge for urban freshwaters and availability of water for humans and wildlife. In order to increase understanding of urban ecohydrology, we investigated the variability of planktonic bacteria and benthic diatoms - as two key biological indicators - coupled with insights from hydrochemistry and stable water isotopes across four urban streams characterized by different dominant water sources in Berlin, the German capital, over a period of one year (2021-2022). DNA metabarcoding results show that substantial spatio-temporal variability exists across urban streams in terms of microbial diversity and richness, with clear links to abiotic factors and nutrient concentrations. Bacterial communities showed clear distinction between effluent-impacted and non-effluent impacted streams as well as clear seasonal turnover. In-stream benthic diatom assemblages also showed robust seasonal variation as well as high species diversity. Our multiple-tracer approach is relevant for emerging questions regarding the increased use of treated effluent to supplement declining baseflows, the assessment of stream restoration projects and the impact of storm drainage and surface pollution on aquatic ecosystem health. eDNA analysis allows analysis of spatial and temporal patterns not feasibly studied with traditional analyses of macroinvertebrates. This can ultimately be leveraged for future water resource management and restoration planning and monitoring of urban freshwater systems across metropolitan areas.

Single-cell genomics reveals new rozellid lineages and supports their sister relationship to Microsporidia

The phylum Rozellomycota has been proposed for a group of early-branching holomycotan lineages representing obligate parasites and hyperparasites of zoosporic fungi, oomycotes or phytoplankton. Given their predominantly intracellular lifestyle, rozellids are typically known from environmental ribosomal DNA data, except for the well-studied Rozella species. To date, the phylogenetic relationship between rozellids and microsporidians (Microsporidia) is not fully understood and most reliable hypotheses are based on phylogenomic analyses that incorporate the only publicly available rozellid genome of Rozella allomycis. Here, we provide genomic data of three new rozellid lineages obtained by single-cell sequencing from environmental samples and show with a phylogenomic approach that rozellids form a monophyletic group that is sister to microsporidians, corroborating the previously proposed phylum Rozellomycota. Whereas no mitochondrial genes coding for the respiratory Complex I could be found, we discovered a gene coding for a nucleotide phosphate transporter in one of the three draft genomes. The scattered absence of Complex I genes and scattered presence of nucleotide transporter genes across diverse microsporidian and rozellid lineages suggest that these adaptations to a parasitic lifestyle, which reduce the parasite's capability to synthesize ATP but enables it to steal ATP from its host, evolved independently in microsporidians and rozellids.

Persisting roadblocks in arthropod monitoring using non-destructive metabarcoding from collection media of passive traps

Background

Broad-scale monitoring of arthropods is often carried out with passive traps (e.g., Malaise traps) that can collect thousands of specimens per sample. The identification of individual specimens requires time and taxonomic expertise, limiting the geographical and temporal scale of research and monitoring studies. DNA metabarcoding of bulk-sample homogenates has been found to be faster, efficient and reliable, but the destruction of samples prevents a posteriori validation of species occurrences and relative abundances. Non-destructive metabarcoding of DNA extracted from collection medium has been applied in a limited number of studies, but further tests of efficiency are required with different trap types and collection media to assess the consistency of the method.

Methods

We quantified the detection rate of arthropod species when applying non-destructive DNA metabarcoding with a short (127-bp) fragment of mitochondrial COI on two combinations of passive traps and collection media: (1) water with monopropylene glycol (H2O-MPG) used in window-flight traps (WFT, 53 in total); (2) ethanol with monopropylene glycol (EtOH-MPG) used in Malaise traps (MT, 27 in total). We then compared our results with those obtained for the same samples using morphological identification (for WFTs) or destructive metabarcoding of bulk homogenate (for MTs). This comparison was applied as part of a larger study of arthropod species richness in silver fir (Abies alba Mill., 1759) stands across a range of climate-induced tree dieback levels and forest management strategies.

Results

Of the 53 H2O-MPG samples from WFTs, 16 produced no metabarcoding results, while the remaining 37 samples yielded 77 arthropod MOTUs in total, of which none matched any of the 343 beetle species morphologically identified from the same traps. Metabarcoding of 26 EtOH-MPG samples from MTs detected more arthropod MOTUs (233) than destructive metabarcoding of homogenate (146 MOTUs, 8 orders), of which 71 were shared MOTUs, though MOTU richness per trap was similar between treatments. While we acknowledge the failure of metabarcoding from WFT-derived collection medium (H2O-MPG), the treatment of EtOH-based Malaise trapping medium remains promising. We conclude however that DNA metabarcoding from collection medium still requires further methodological developments and cannot replace homogenate metabarcoding as an approach for arthropod monitoring. It can be used nonetheless as a complementary treatment when enhancing the detection of soft-bodied arthropods like spiders and Diptera.

Spatial and phylogenetic structure of Alpine stonefly assemblages across seven habitats using DNA-species

Stream ecosystems are spatially heterogeneous, with many different habitat patches distributed within a small area. The influence of this heterogeneity on the biodiversity of benthic insect communities is well documented; however, studies of the role of habitat heterogeneity in species coexistence and assembly remain limited. Here, we investigated how habitat heterogeneity influences spatial structure (beta biodiversity) and phylogenetic structure (evolutionary processes) of benthic stonefly (Plecoptera, Insecta) communities. We sampled 20 sites along two Alpine rivers, including seven habitats in four different reaches (headwaters, meandering, bar-braided floodplain, and lowland spring-fed). We identified 21 morphological species and delineated 52 DNA-species based on sequences from mitochondrial cox1 and nuclear ITS markers. Using DNA-species, we first analysed the patterns of variation in richness, diversity, and assemblage composition by quantifing the contribution of each reach and habitat to the overall DNA-species diversity using an additive partition analysis and distance-based redundancy analysis. Using gene-tree phylogenies, we assessed whether environmental filtering could lead to the co-occurrence of DNA-species using a two-step analysis to detect a phylogenetic signal. All four reaches significantly contributed to DNA-species richness, with the meandering reach having the highest contribution. Habitats had an effect on DNA-species diversity, where glide, riffle and, pool influenced the spatial structure of stonefly assemblage possibly due to the high habitat heterogeneity. Among the habitats, the pool showed significant phylogenetic clustering, suggesting high levels of evolutionary adaptation and strong habitat filtering. This assemblage structure may be caused by long-term stability of the habitat and the similar requirements for co-occurring species. Our study shows the importance of different habitats for the spatial and phylogenetic structure of stonefly assemblage and sheds light on the habitat-specific diversity that may help improve conservation practices.

Large-scale sampling of the freshwater microbiome suggests pollution-driven ecosystem changes

Freshwater microbes play a crucial role in the global carbon cycle. Anthropogenic stressors that lead to changes in these microbial communities are likely to have profound consequences for freshwater ecosystems. Using field data from the coordinated sampling of 617 lakes, ponds, rivers, and streams by citizen scientists, we observed linkages between microbial community composition, light and chemical pollution, and greenhouse gas concentration. All sampled water bodies were net emitters of CO₂, with higher concentrations in running waters, and increasing concentrations at higher latitudes. Light pollution occurred at 75% of sites, was higher in urban areas and along rivers, and had a measurable effect on the microbial alpha diversity. Genetic elements suggestive of chemical stress and antimicrobial resistances (IntI1, blaOX₅₈) were found in 85% of sites, and were also more prevalent in urban streams and rivers. Light pollution and CO₂ were significantly related to microbial community composition, with CO₂ inversely related to microbial phototrophy. Results of synchronous nationwide sampling indicate that pollution-driven alterations to the freshwater microbiome lead to changes in CO₂ production in natural waters and highlight the vulnerability of running waters to anthropogenic stressors.

Phylogenomic insights into the early diversification of fungi

Phylogenomic analyses have boosted our understanding of the evolutionary trajectories of all living forms by providing continuous improvements to the tree of life.^1-5 Within this tree, fungi represent an ancient eukaryote group,⁶ having diverged from the animals ∼1.35 billion years ago.⁷ Estimates of the number of extant species range between 1.5 and 3.8 million.^8,9 Recent reclassifications and the discovery of the deep-branching Sanchytriomycota lineage¹⁰ have brought the number of proposed phyla to 20,¹¹ 21 if the Microsporidia are included.^12-14 Uncovering how the diverse and globally distributed fungi are related to each other is fundamental for understanding how their lifestyles, morphologies, and metabolic capacities evolved. To date, many of the proposed relationships among the phyla remain controversial and no phylogenomic study has examined the entire fungal tree using a taxonomically comprehensive dataset and suitable models of evolution. We assembled and curated a 299-protein dataset with a taxon sampling broad enough to encompass all recognized fungal diversity with available data, but selective enough to run computationally intensive analyses using best-fitting models. Using a range of reconstruction methods, we were able to resolve many contested nodes, such as a sister relationship of Chytridiomyceta to all other non-Opisthosporidia fungi (with Chytridiomycota being sister to Monoblepharomycota + Neocallimastigomycota), a branching of Blastocladiomycota + Sanchytriomycota after the Chytridiomyceta but before other non-Opisthosporidia fungi, and a branching of Glomeromycota as sister to the Dikarya. Our up-to-date fungal tree of life will serve as a springboard for future investigations on the early evolution of fungi.

Climate-induced forest dieback drives compositional changes in insect communities that are more pronounced for rare species

Species richness, abundance and biomass of insects have recently undergone marked declines in Europe. We metabarcoded 211 Malaise-trap samples to investigate whether drought-induced forest dieback and subsequent salvage logging had an impact on ca. 3000 species of flying insects in silver fir Pyrenean forests. While forest dieback had no measurable impact on species richness, there were significant changes in community composition that were consistent with those observed during natural forest succession. Importantly, most observed changes were driven by rare species. Variation was explained primarily by canopy openness at the local scale, and the tree-related microhabitat diversity and deadwood amount at landscape scales. The levels of salvage logging in our study did not explain compositional changes. We conclude that forest dieback drives changes in species assemblages that mimic natural forest succession, and markedly increases the risk of catastrophic loss of rare species through homogenization of environmental conditions.

A global agenda for advancing freshwater biodiversity research

Global freshwater biodiversity is declining dramatically, and meeting the challenges of this crisis requires bold goals and the mobilisation of substantial resources. While the reasons are varied, investments in both research and conservation of freshwater biodiversity lag far behind those in the terrestrial and marine realms. Inspired by a global consultation, we identify 15 pressing priority needs, grouped into five research areas, in an effort to support informed stewardship of freshwater biodiversity. The proposed agenda aims to advance freshwater biodiversity research globally as a critical step in improving coordinated actions towards its sustainable management and conservation.

Evidence for Lignocellulose-Decomposing Enzymes in the Genome and Transcriptome of the Aquatic Hyphomycete Clavariopsis aquatica

Fungi are ecologically outstanding decomposers of lignocellulose. Fungal lignocellulose degradation is prominent in saprotrophic Ascomycota and Basidiomycota of the subkingdom Dikarya. Despite ascomycetes dominating the Dikarya inventory of aquatic environments, genome and transcriptome data relating to enzymes involved in lignocellulose decay remain limited to terrestrial representatives of these phyla. We sequenced the genome of an exclusively aquatic ascomycete (the aquatic hyphomycete Clavariopsis aquatica), documented the presence of genes for the modification of lignocellulose and its constituents, and compared differential gene expression between C. aquatica cultivated on lignocellulosic and sugar-rich substrates. We identified potential peroxidases, laccases, and cytochrome P450 monooxygenases, several of which were differentially expressed when experimentally grown on different substrates. Additionally, we found indications for the regulation of pathways for cellulose and hemicellulose degradation. Our results suggest that C. aquatica is able to modify lignin to some extent, detoxify aromatic lignin constituents, or both. Such characteristics would be expected to facilitate the use of carbohydrate components of lignocellulose as carbon and energy sources.

Three mitochondrial genomes of early-winged insects (Ephemeroptera: Baetidae and Leptophlebiidae)

Mayflies (Ephemeroptera) are a semi-aquatic insect order with comparatively few genomic data available despite their phylogenetic position at the root of the winged-insects and possession of ancestral traits. Here, we provide three mitochondrial genomes (mtgenomes) from representatives of the two most species-rich families, Baetis rutilocylindratus and Cloeon dipterum (Baetidae), and Habrophlebiodes zijinensis (Leptophlebiidae). All mtgenomes had a complete set of 13 protein-coding genes and a conserved orientation except for two inverted tRNAs in H. zijinensis. Phylogenetic reconstructions using 21 mayfly mtgenomes and representatives of seven additional orders recovered both Baetidae and Leptophlebiidae as well supported monophyletic clades, with Ephemeroptera as the sister-taxon to all other winged insects (i.e. Odonata and Neoptera).

Comparative population genetic structure of two ixodid tick species (Acari:Ixodidae) (Ixodes ovatus and Haemaphysalis flava) in Niigata prefecture, Japan

Ixodid ticks (Acari:Ixodidae) are essential vectors of tick-borne diseases in Japan. In this study, we characterized the population genetic structure and inferred genetic divergence in two widespread and abundant ixodid species, Ixodes ovatus and Haemaphysalis flava. Our hypothesis was that genetic divergence would be high in I. ovatus because of the low mobility of their small rodent hosts of immature I. ovatus would limit their gene flow compared to more mobile avian hosts of immature H. flava. We collected 320 adult I. ovatus from 29 locations and 223 adult H. flava from 17 locations across Niigata Prefecture, Japan, and investigated their genetic structure using DNA sequences from fragments of two mitochondrial gene regions, cox1 and the 16S rRNA gene. For I. ovatus, pairwise F_ST and analysis of molecular variance (AMOVA) analyses of cox1 and 16S sequences indicated significant genetic variation among populations, whereas both markers showed non-significant genetic variation among locations for H. flava. A cox1 gene tree and haplotype network revealed three genetic groups of I. ovatus. One of these groups consisted of haplotypes distributed at lower altitudes (251-471 m.a.s.l.). The cox1 sequences of I. ovatus from Japan clustered separately from I. ovatus sequences reported from China, suggesting the potential for cryptic species in Japan. Our results support our hypothesis and suggest that the host preference of ticks at the immature stage may influence the genetic structure of the ticks. This information may be important for understanding the tick-host interactions in the field to better understand the tick-borne disease transmission and in designing an effective tick control program.

Encapsulation of Anguillicola crassus reduces the abundance of adult parasite stages in the European eel (Anguilla anguilla)

Encapsulation of the parasitic nematode Anguillicola crassus Kuwahara, Niimi & Hagaki is commonly observed in its native host, the Japanese eel (Anguilla japonica Temminck & Schlegel). Encapsulation has also been described in a novel host, the European eel (A. anguilla L.), and there is evidence that encapsulation frequency has increased since the introduction of A. crassus. We examined whether encapsulation of A. crassus provides an advantage to its novel host in Lake Müggelsee, NE Germany. We provide the first evidence that encapsulation was associated with reduced abundance of adult A. crassus. This pattern was consistent in samples taken 3 months apart. There was no influence of infection on the expression of the two metabolic genes studied, but the number of capsules was negatively correlated with the expression of two mhc II genes of the adaptive immune response, suggesting a reduced activation. Interestingly, eels that encapsulated A. crassus had higher abundances of two native parasites compared with non-encapsulating eels. We propose that the response of A. anguilla to infection by A. crassus may interfere with its reaction to other co-occurring parasites.

Trade-offs between reducing complex terminology and producing accurate interpretations from environmental DNA: Comment on "Environmental DNA: What's behind the term?" by Pawlowski et al., (2020)

In a recent paper, “Environmental DNA: What's behind the term? Clarifying the terminology and recommendations for its future use in biomonitoring,” Pawlowski et al. argue that the term eDNA should be used to refer to the pool of DNA isolated from environmental samples, as opposed to only extra‐organismal DNA from macro‐organisms. We agree with this view. However, we are concerned that their proposed two‐level terminology specifying sampling environment and targeted taxa is overly simplistic and might hinder rather than improve clear communication about environmental DNA and its use in biomonitoring. This terminology is based on categories that are often difficult to assign and uninformative, and it overlooks a fundamental distinction within eDNA: the type of DNA (organismal or extra‐organismal) from which ecological interpretations are derived.

Revisiting the phylogenetic position of Caullerya mesnili (Ichthyosporea), a common Daphnia parasite, based on 22 protein-coding genes

Caullerya mesnili is a common and virulent parasite of the water flea, Daphnia. It was classified within the Haplosporidia (Rhizaria) for over a century. However, a recent molecular phylogeny based on the 18S rRNA gene suggested it belonged to the Ichthyosporea, a class of protists closely related to animals within the Opisthokonta clade. The exact phylogenetic position of C. mesnili remained uncertain because it appeared in the 18S rRNA tree with a very long branch and separated from all other taxa, suggesting that its position could be artifactual. A better understanding of its phylogenetic position has been constrained by a lack of molecular markers and the difficulty of obtaining a suitable quantity and quality of DNA from in vitro cultures, as this intracellular parasite cannot be cultured without its host. We isolated and collected spores of C. mesnili and sequenced genomic libraries. Phylogenetic analyses of a newly generated multi-protein data set (22 proteins, 4998 amino acids) and of sequences from the 18S rRNA gene both placed C. mesnili within the Ichthyophonida sub-clade of Ichthyosporea, as sister-taxon to Abeoforma whisleri and Pirum gemmata. Our study highlights the utility of metagenomic approaches for obtaining genomic information from intracellular parasites and for more accurate phylogenetic placement in evolutionary studies.

Gene expression response to a nematode parasite in novel and native eel hosts

Invasive parasites are involved in population declines of new host species worldwide. The high susceptibilities observed in many novel hosts have been attributed to the lack of protective immunity to the parasites which native hosts acquired during their shared evolution. We experimentally infected Japanese eels (Anguilla japonica) and European eels (Anguilla anguilla) with Anguillicola crassus, a nematode parasite that is native to the Japanese eel and invasive in the European eel. We inferred gene expression changes in head kidney tissue from both species, using RNA-seq data to determine the responses at two time points during the early stages of infection (3 and 23 days postinfection). At both time points, the novel host modified the expression of a larger and functionally more diverse set of genes than the native host. Strikingly, the native host regulated immune gene expression only at the earlier time point and to a small extent while the novel host regulated these genes at both time points. A low number of differentially expressed immune genes, especially in the native host, suggest that a systemic immune response was of minor importance during the early stages of infection. Transcript abundance of genes involved in cell respiration was reduced in the novel host which may affect its ability to cope with harsh conditions and energetically demanding activities. The observed gene expression changes in response to a novel parasite that we observed in a fish follow a general pattern observed in amphibians and mammals, and suggest that the disruption of physiological processes, rather than the absence of an immediate immune response, is responsible for the higher susceptibility of the novel host.

Mutual Facilitation Among Invading Nuttall's Waterweed and Quagga Mussels

Nuttall's waterweed (Elodea nuttallii) is the most abundant invasive aquatic plant species in several European countries. Elodea populations often follow a boom-bust cycle, but the causes and consequences of this dynamics are yet unknown. We hypothesize that both boom and bust periods can be affected by dreissenid mussel invasions. While mutual facilitations between these invaders could explain their rapid parallel expansion, subsequent competition for space might occur. To test this hypothesis, we use data on temporal changes in the water quality and the abundance of E. nuttallii and the quagga mussel Dreissena rostriformis bugensis in a temperate shallow lake. Lake Müggelsee (Germany) was turbid and devoid of submerged macrophytes for 20 years (1970-1989), but re-colonization with macrophytes started in 1990 upon reductions in nutrient loading. We mapped macrophyte abundance from 1999 and mussel abundance from 2011 onwards. E. nuttallii was first detected in 2011, spread rapidly, and was the most abundant macrophyte species by 2017. Native macrophyte species were not replaced, but spread more slowly, resulting in an overall increase in macrophyte coverage to 25% of the lake surface. The increased abundance of E. nuttallii was paralleled by increasing water clarity and decreasing total phosphorus concentrations in the water. These changes were attributed to a rapid invasion by quagga mussels in 2012. In 2017, they covered about one-third of the lake area, with mean abundances of 3,600 mussels m^-2, filtering up to twice the lake's volume every day. The increasing light availability in deeper littoral areas supported the rapid spread of waterweed, while in turn waterweed provided surface for mussel colonization. Quantities of dreissenid mussels and E. nuttallii measured at 24 locations were significantly correlated in 2016, and yearly means of E. nuttallii quantities increased with increasing mean dreissenid mussel quantities between 2011 and 2018. In 2018, both E. nuttallii and dreissenid abundances declined. These data imply that invasive waterweed and quagga mussels initially facilitated their establishment, supporting the invasional meltdown hypothesis, while subsequently competition for space may have occurred. Such temporal changes in invasive species interaction might contribute to the boom-bust dynamics that have been observed in Elodea populations.

Daphnia galeata responds to the exposure to an ichthyosporean gut parasite by down-regulation of immunity and lipid metabolism

BACKGROUND

Regulatory circuits of infection in the emerging experimental model system, water flea Daphnia and their microparasites, remain largely unknown. Here we provide the first molecular insights into the response of Daphnia galeata to its highly virulent and common parasite Caullerya mesnili, an ichthyosporean that infects the gut epithelium. We generated a transcriptomic dataset using RNAseq from parasite-exposed (vs. control) Daphnia, at two time points (4 and 48 h) after parasite exposure.

RESULTS

We found a down-regulation of metabolism and immunity-related genes, at 48 h (but not 4 h) after parasite exposure. These genes are involved in lipid metabolism and fatty acid biosynthesis, as well as microbe recognition (e.g. c-type lectins) and pathogen attack (e.g. gut chitin).

CONCLUSIONS

General metabolic suppression implies host energy shift from reproduction to survival, which is in agreement with the known drastic reduction in Daphnia fecundity after Caullerya infection. The down-regulation of gut chitin indicates a possible interaction between the peritrophic matrix and the evading host immune system. Our study provides the first description of host transcriptional responses in this very promising host-parasite experimental system.

A transition to white LED increases ecological impacts of nocturnal illumination on aquatic primary producers in a lowland agricultural drainage ditch

The increasing use of artificial light at night (ALAN) has led to exposure of freshwater ecosystems to light pollution worldwide. Simultaneously, the spectral composition of nocturnal illumination is changing, following the current shift in outdoor lighting technologies from traditional light sources to light emitting diodes (LED). LEDs emit broad-spectrum white light, with a significant amount of photosynthetically active radiation, and typically a high content of blue light that regulates circadian rhythms in many organisms. While effects of the shift to LED have been investigated in nocturnal animals, its impact on primary producers is unknown. We performed three field experiments in a lowland agricultural drainage ditch to assess the impacts of a transition from high-pressure sodium (HPS) to white LED illumination (color temperature 4000 K) on primary producers in periphyton. In all experiments, we compared biomass and pigment composition of periphyton grown under a natural light regime to that of periphyton exposed to nocturnal HPS or, consecutively, LED light of intensities commonly found in urban waters (approximately 20 lux). Periphyton was collected in time series (1-13 weeks). We found no effect of HPS light on periphyton biomass; however, following a shift to LED the biomass decreased up to 62%. Neither light source had a substantial effect on pigment composition. The contrasting effects of the two light sources on biomass may be explained by differences in their spectral composition, and in particular the blue content. Our results suggest that spectral composition of the light source plays a role in determining the impacts of ALAN on periphyton and that the ongoing transition to LED may increase the ecological impacts of artificial lighting on aquatic primary producers. Reduced biomass in the base of the food web can impact ecosystem functions such as productivity and food supply for higher trophic levels in nocturnally-lit ecosystems.

City-Dwellers and Country Folks: Lack of Population Differentiation Along an Urban-Rural Gradient in the Mosquito Culex pipiens (Diptera: Culicidae)

Mosquitoes (Diptera, Culicidae) occur in natural, urban, and peri-urban areas throughout the globe. Although the characteristics of urban and peri-urban habitats differ from those of natural habitats in many ways (e.g., fragmentation, pollution, noise, and light), few studies have examined the population connectivity of mosquitoes in urban areas. To obtain an overview of the species composition, we sampled mosquitoes from 23 sites in and around the city of Berlin, Germany. Of 23 species, five occurred in urban, 10 in peri-urban, and 20 in rural areas. Culex pipiens Linnaeus (Diptera: Culicidae) was the most common species collected (75% of all individuals) and occurred in all habitats. Hence this species was selected to be analysed at 10 microsatellite markers. There was no significant differentiation (FST = 0.016, P = 0.9) or isolation by distance (P = 0.06) among Cx. pipiens populations along an urban-rural gradient. The only significant differences detected were between Cx. pipiens and a laboratory population of Cx. pipiens f. molestus (pairwise FST = 0.114-0.148, P ≤ 0.001 in all comparisons), suggesting that the markers chosen were suitable for the identification of population differentiation. Our results indicate that Cx. pipiens gene flow is widespread within and among urban, peri-urban, and rural areas and that urban habitat does not necessarily impede or enhance gene flow among these populations.

DNA metabarcoding of unfractionated water samples relates phyto-, zoo- and bacterioplankton dynamics and reveals a single-taxon bacterial bloom

Most studies of aquatic plankton focus on either macroscopic or microbial communities, and on either eukaryotes or prokaryotes. This separation is primarily for methodological reasons, but can overlook potential interactions among groups. Here we tested whether DNA metabarcoding of unfractionated water samples with universal primers could be used to qualitatively and quantitatively study the temporal dynamics of the total plankton community in a shallow temperate lake. Significant changes in the relative proportions of normalized sequence reads of eukaryotic and prokaryotic plankton communities over a 3-month period in spring were found. Patterns followed the same trend as plankton estimates measured using traditional microscopic methods. The bloom of a conditionally rare bacterial taxon belonging to Arcicella was characterized, which rapidly came to dominate the whole lake ecosystem and would have remained unnoticed without metabarcoding. The data demonstrate the potential of universal DNA metabarcoding applied to unfractionated samples for providing a more holistic view of plankton communities.

DNA metabarcoding of unfractionated water samples relates phyto-, zoo- and bacterioplankton dynamics and reveals a single-taxon bacterial bloom

Most studies of aquatic plankton focus on either macroscopic or microbial communities, and on either eukaryotes or prokaryotes. This separation is primarily for methodological reasons, but can overlook potential interactions among groups. Here we tested whether DNA metabarcoding of unfractionated water samples with universal primers could be used to qualitatively and quantitatively study the temporal dynamics of the total plankton community in a shallow temperate lake. Significant changes in the relative proportions of normalized sequence reads of eukaryotic and prokaryotic plankton communities over a 3-month period in spring were found. Patterns followed the same trend as plankton estimates measured using traditional microscopic methods. The bloom of a conditionally rare bacterial taxon belonging to Arcicella was characterized, which rapidly came to dominate the whole lake ecosystem and would have remained unnoticed without metabarcoding. The data demonstrate the potential of universal DNA metabarcoding applied to unfractionated samples for providing a more holistic view of plankton communities.

Comparative tests of the species-genetic diversity correlation at neutral and nonneutral loci in four species of stream insect

A fundamental question linking population genetics and community ecology is how adaptive processes (e.g., natural selection) and neutral processes (e.g., drift-migration equilibrium) underpin the species-genetic diversity correlation (SGDC). Here, we combine genome scans and outlier loci detection with community analysis to separately test for neutral and nonneutral SGDCs in four species of stream insect. We sampled 60 localities in Japan and examined the relationships among population AFLP band richness (Br), taxon richness of the total community (S) and of the trophic guild (S_tr ), and 15 habitat parameters that could potentially drive adaptation and influence richness. Neutral Br was positively correlated with S only in the dominant species of these communities, suggesting Br may be constrained when intraspecific competition is pronounced. Nonneutral Br was correlated with S_tr in a species restricted to high elevations where habitat heterogeneity was highest. Community distance and genetic distance (β-SGDC) was correlated in two of the four species at both neutral and nonneutral loci. Distance-based redundancy analysis found geographic isolation and elevation to drive divergence of both communities and populations. This suggests that both neutral and adaptive divergence occurred through the shared influences of geographic isolation and local adaptation at the two levels of diversity.

Shifts among Eukaryota, Bacteria, and Archaea define the vertical organization of a lake sediment

BACKGROUND

Lake sediments harbor diverse microbial communities that cycle carbon and nutrients while being constantly colonized and potentially buried by organic matter sinking from the water column. The interaction of activity and burial remained largely unexplored in aquatic sediments. We aimed to relate taxonomic composition to sediment biogeochemical parameters, test whether community turnover with depth resulted from taxonomic replacement or from richness effects, and to provide a basic model for the vertical community structure in sediments.

METHODS

We analyzed four replicate sediment cores taken from 30-m depth in oligo-mesotrophic Lake Stechlin in northern Germany. Each 30-cm core spanned ca. 170 years of sediment accumulation according to ¹³⁷Cs dating and was sectioned into layers 1-4 cm thick. We examined a full suite of biogeochemical parameters and used DNA metabarcoding to examine community composition of microbial Archaea, Bacteria, and Eukaryota.

RESULTS

Community β-diversity indicated nearly complete turnover within the uppermost 30 cm. We observed a pronounced shift from Eukaryota- and Bacteria-dominated upper layers (<5 cm) to Bacteria-dominated intermediate layers (5-14 cm) and to deep layers (>14 cm) dominated by enigmatic Archaea that typically occur in deep-sea sediments. Taxonomic replacement was the prevalent mechanism in structuring the community composition and was linked to parameters indicative of microbial activity (e.g., CO₂ and CH₄ concentration, bacterial protein production). Richness loss played a lesser role but was linked to conservative parameters (e.g., C, N, P) indicative of past conditions.

CONCLUSIONS

By including all three domains, we were able to directly link the exponential decay of eukaryotes with the active sediment microbial community. The dominance of Archaea in deeper layers confirms earlier findings from marine systems and establishes freshwater sediments as a potential low-energy environment, similar to deep sea sediments. We propose a general model of sediment structure and function based on microbial characteristics and burial processes. An upper "replacement horizon" is dominated by rapid taxonomic turnover with depth, high microbial activity, and biotic interactions. A lower "depauperate horizon" is characterized by low taxonomic richness, more stable "low-energy" conditions, and a dominance of enigmatic Archaea.

Colonization and diversification of aquatic insects on three Macaronesian archipelagos using 59 nuclear loci derived from a draft genome

The study of processes driving diversification requires a fully sampled and well resolved phylogeny, although a lack of phylogenetic markers remains a limitation for many non-model groups. Multilocus approaches to the study of recent diversification provide a powerful means to study the evolutionary process, but their application remains restricted because multiple unlinked loci with suitable variation for phylogenetic or coalescent analysis are not available for most non-model taxa. Here we identify novel, putative single-copy nuclear DNA (nDNA) phylogenetic markers to study the colonization and diversification of an aquatic insect species complex, Cloeon dipterum L. 1761 (Ephemeroptera: Baetidae), in Macaronesia. Whole-genome sequencing data from one member of the species complex were used to identify 59 nDNA loci (32,213 base pairs), followed by Sanger sequencing of 29 individuals sampled from 13 islands of three Macaronesian archipelagos. Multispecies coalescent analyses established six putative species. Three island species formed a monophyletic clade, with one species occurring on the Azores, Europe and North America. Ancestral state reconstruction indicated at least two colonization events from the mainland (to the Canaries, respectively Azores) and one within the archipelago (between Madeira and the Canaries). Random subsets of the 59 loci showed a positive linear relationship between number of loci and node support. In contrast, node support in the multispecies coalescent tree was negatively correlated with mean number of phylogenetically informative sites per locus, suggesting a complex relationship between tree resolution and marker variability. Our approach highlights the value of combining genomics, coalescent-based phylogeography, species delimitation, and phylogenetic reconstruction to resolve recent diversification events in an archipelago species complex.

Molecular phylogeny and timing of diversification in Alpine Rhithrogena (Ephemeroptera: Heptageniidae)

Background

Larvae of the Holarctic mayfly genus Rhithrogena Eaton, 1881 (Ephemeroptera, Heptageniidae) are a diverse and abundant member of stream and river communities and are routinely used as bio-indicators of water quality. Rhithrogena is well diversified in the European Alps, with a number of locally endemic species, and several cryptic species have been recently detected. While several informal species groups are morphologically well defined, a lack of reliable characters for species identification considerably hampers their study. Their relationships, origin, timing of speciation and mechanisms promoting their diversification in the Alps are unknown.

Results

Here we present a species-level phylogeny of Rhithrogena in Europe using two mitochondrial and three nuclear gene regions. To improve sampling in a genus with many cryptic species, individuals were selected for analysis according to a recent DNA-based taxonomy rather than traditional nomenclature. A coalescent-based species tree and a reconstruction based on a supermatrix approach supported five of the species groups as monophyletic. A molecular clock, mapped on the most resolved phylogeny and calibrated using published mitochondrial evolution rates for insects, suggested an origin of Alpine Rhithrogena in the Oligocene/Miocene boundary. A diversification analysis that included simulation of missing species indicated a constant speciation rate over time, rather than any pronounced periods of rapid speciation. Ancestral state reconstructions provided evidence for downstream diversification in at least two species groups.

Conclusions

Our species-level analyses of five gene regions provide clearer definitions of species groups within European Rhithrogena. A constant speciation rate over time suggests that the paleoclimatic fluctuations, including the Pleistocene glaciations, did not significantly influence the tempo of diversification of Alpine species. A downstream diversification trend in the hybrida and alpestris species groups supports a previously proposed headwater origin hypothesis for aquatic insects.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0758-1) contains supplementary material, which is available to authorized users.

discomark: nuclear marker discovery from orthologous sequences using draft genome data

High-throughput sequencing has laid the foundation for fast and cost-effective development of phylogenetic markers. Here we present the program discomark, which streamlines the development of nuclear DNA (nDNA) markers from whole-genome (or whole-transcriptome) sequencing data, combining local alignment, alignment trimming, reference mapping and primer design based on multiple sequence alignments to design primer pairs from input orthologous sequences. To demonstrate the suitability of discomark, we designed markers for two groups of species, one consisting of closely related species and one group of distantly related species. For the closely related members of the species complex of Cloeon dipterum s.l. (Insecta, Ephemeroptera), the program discovered a total of 78 markers. Among these, we selected eight markers for amplification and Sanger sequencing. The exon sequence alignments (2526 base pairs) were used to reconstruct a well-supported phylogeny and to infer clearly structured haplotype networks. For the distantly related species, we designed primers for the insect order Ephemeroptera, using available genomic data from four sequenced species. We developed primer pairs for 23 markers that are designed to amplify across several families. The discomark program will enhance the development of new nDNA markers by providing a streamlined, automated approach to perform genome-scale scans for phylogenetic markers. The program is written in Python, released under a public licence (GNU GPL version 2), and together with a manual and example data set available at: https://github.com/hdetering/discomark.

Microbial diversity and community respiration in freshwater sediments influenced by artificial light at night

An increasing proportion of the Earth's surface is illuminated at night. In aquatic ecosystems, artificial light at night (ALAN) may influence microbial communities living in the sediments. These communities are highly diverse and play an important role in the global carbon cycle. We combined field and laboratory experiments using sediments from an agricultural drainage system to examine how ALAN affects communities and alters carbon mineralization. Two identical light infrastructures were installed parallel to a drainage ditch before the start of the experiment. DNA metabarcoding indicated that both sediment communities were similar. After one was lit for five months (July-December 2012) we observed an increase in photoautotroph abundance (diatoms, Cyanobacteria) in ALAN-exposed sediments. In laboratory incubations mimicking summer and winter (six weeks each), communities in sediments that were exposed to ALAN for 1 year (July 2012-June 2013) showed less overall seasonal change compared with ALAN-naive sediments. Nocturnal community respiration was reduced in ALAN-exposed sediments. In long-term exposed summer-sediments, we observed a shift from negative to positive net ecosystem production. Our results indicate ALAN may alter sediment microbial communities over time, with implications for ecosystem-level functions. It may thus have the potential to transform inland waters to nocturnal carbon sinks.

Sex-specific gene expression in the mosquito Culex pipiens f. molestus in response to artificial light at night

Background

Artificial light at night (ALAN) is a typical feature of urban areas and most organisms living in urban or suburban habitats are exposed to low levels of ALAN. Light is one of the most important environmental cues that organisms use to regulate their activities. Studies have begun to quantify the influence of ALAN on the behavior and ecology of organisms, but research on the effects at the molecular level remains limited. Mosquitoes in the Culex pipiens complex (Diptera, Culicidae) are widespread and abundant in urban areas where they are potential disease vectors. It is thus of particular interest to understand how ALAN may influence biologically and ecologically relevant traits.

Results

We used RNAseq to evaluate the transcriptome response in a Cx. pipiens f. molestus laboratory population that was exposed to near-natural light conditions (light:dark L16:D8 hours, “control”) and ALAN conditions with 3 h of constant low-level light at night (L16 + L_low3:D5 hours, “low-light”). The resulting transcripts were mapped to the reference genome of the closely related Culex quinquefasciatus. Female expression patterns differed significantly between control and treatment conditions at five genes although none showed an absolute fold change greater than two (FC > 2). In contrast, male expression differed at 230 genes (74 with FC > 2). Of these, 216 genes (72 with FC > 2) showed reduced expression in the low-light treatment, most of which were related to gametogenesis, lipid metabolism, and immunity. Of the 14 genes (two with FC > 2) with increased expression, only five had any functional annotation. There was a pronounced sex-bias in gene expression regardless of treatment, with 11,660 genes (51 % of annotated genes; 8694 with FC > 2; 48 % of annotated genes) differentially expressed between males and females, including 14 genes of the circadian clock.

Conclusion

Our data suggest a stronger response to artificial light by males of Cx. pipiens f. molestus than by females, and that a wide range of physiological pathways may be affected by ALAN at the molecular level. The fact that differences in gene expression appear to be sex-specific may have a strong influence at the population level.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2336-0) contains supplementary material, which is available to authorized users.

Occurrence and host preferences of Anopheles maculipennis group mosquitoes in England and Wales

Mosquitoes of the Anopheles maculipennis Meigen (Diptera: Culicidae) group are of public health concern: five of the 11 morphologically indistinct species have been historically considered as vectors of malaria in Europe. Three members of the An. maculipennis group have been reported in the U.K.: Anopheles atroparvus van Thiel; Anopheles messeae Falleroni, and Anopheles daciae Linton, Nicolescu & Harbach. To study the distribution of the three U.K. species, particularly that of An. daciae, we developed a polymerase chain reaction-Restriction fragment length polymorphism (PCR-RFLP) assay using the nuclear ribosomal internal transcribed spacer 2 (ITS-2) gene. Anopheles daciae was found to be widespread, occurring in four of the five counties surveyed in southern England and on the Welsh island of Anglesey, often in sympatry with the closely related species An. messeae. The host preferences of 237 blood-fed females were determined using either direct sequencing or PCR-based fragment analysis of the mitochondrial cytochrome oxidase b gene with DNA from females' abdomens. All three species were found to be opportunistic, having fed on at least three different hosts. Seventeen individuals contained multiple bloodmeals, including two An. daciae that had fed on humans and birds. Our results show that An. daciae is widespread in England and Wales, occurs in sympatry with other members of the An. maculipennis group, and feeds on humans, which suggests it is a potential vector of disease in the U.K.

Adaptive genetic divergence along narrow environmental gradients in four stream insects

A central question linking ecology with evolutionary biology is how environmental heterogeneity can drive adaptive genetic divergence among populations. We examined adaptive divergence of four stream insects from six adjacent catchments in Japan by combining field measures of habitat and resource components with genome scans of non-neutral Amplified Fragment Length Polymorphism (AFLP) loci. Neutral genetic variation was used to measure gene flow and non-neutral genetic variation was used to test for adaptive divergence. We identified the environmental characteristics contributing to divergence by comparing genetic distances at non-neutral loci between sites with Euclidean distances for each of 15 environmental variables. Comparisons were made using partial Mantel tests to control for geographic distance. In all four species, we found strong evidence for non-neutral divergence along environmental gradients at between 6 and 21 loci per species. The relative contribution of these environmental variables to each species' ecological niche was quantified as the specialization index, S, based on ecological data. In each species, the variable most significantly correlated with genetic distance at non-neutral loci was the same variable along which each species was most narrowly distributed (i.e., highest S). These were gradients of elevation (two species), chlorophyll-a, and ammonia-nitrogen. This adaptive divergence occurred in the face of ongoing gene flow (Fst = 0.01-0.04), indicating that selection was strong enough to overcome homogenization at the landscape scale. Our results suggest that adaptive divergence is pronounced, occurs along different environmental gradients for different species, and may consistently occur along the narrowest components of species' niche.

Freshwater biodiversity and aquatic insect diversification

Inland waters cover less than 1% of Earth's surface but harbor more than 6% of all insect species: Nearly 100,000 species from 12 orders spend one or more life stages in freshwater. Little is known about how this remarkable diversity arose, although allopatric speciation and ecological adaptation are thought to be primary mechanisms. Freshwater habitats are highly susceptible to environmental change and exhibit marked ecological gradients. Standing waters appear to harbor more dispersive species than running waters, but there is little understanding of how this fundamental ecological difference has affected diversification. In contrast to the lack of evolutionary studies, the ecology and habitat preferences of aquatic insects have been intensively studied, in part because of their widespread use as bioindicators. The combination of phylogenetics with the extensive ecological data provides a promising avenue for future research, making aquatic insects highly suitable models for the study of ecological diversification.

Effects of predation pressure and resource use on morphological divergence in omnivorous prey fish

Background

Body shape is one of the most variable traits of organisms and responds to a broad array of local selective forces. In freshwater fish, divergent body shapes within single species have been repeatedly observed along the littoral-pelagic axes of lakes, where the structural complexity of near shore habitats provides a more diverse set of resources compared to the open-water zones. It remains poorly understood whether similar resource-driven polymorphism occurs among lakes that vary in structural complexity and predation pressure, and whether this variation is heritable. Here, we analyzed body shape in four populations of omnivorous roach (Rutilus rutilus) inhabiting shallow lakes. We tested the relationship between body shape, gradients of resources, predation pressure, and, in a subset of two lakes, diet composition. We used genome scans of 331 polymorphic AFLP markers to test whether there was a heritable component to the observed morphological diversification.

Results

Body shape differed among lakes and was significantly correlated to differences in predation pressure. Roach from the lake with highest predation pressure were most divergent from the average body shape of all populations, characterized by a more streamlined body and caudally inserted dorsal fins; features that facilitate predator escape. Surprisingly, diet composition was not associated with morphology. AFLP analysis revealed weak genetic differentiation among lakes and no isolation by distance (IBD). Outlier analysis detected three loci under positive selection with differing frequencies in the four populations. General linear models did not support an association of lake-specific genotypes with morphological variation.

Conclusion

Body shape was divergent among lakes, suggesting that processes previously reported from within single lakes may also be operating at the scale of whole lakes. We found no evidence for body shape being heritable, although sample size was small in these natural populations. Rather than habitat structure and diet, we conclude that predation had a stronger effect on the prevalence of local morphotypes. A variable morphotype facilitating the efficient uptake of a variety of spatially and temporarily scattered resources seems to be favored in these small aquatic systems.

Endemism and diversification in freshwater insects of Madagascar revealed by coalescent and phylogenetic analysis of museum and field collections

The biodiversity and endemism of Madagascar are among the most extraordinary and endangered in the world. This includes the island's freshwater biodiversity, although detailed knowledge of the diversity, endemism, and biogeographic origin of freshwater invertebrates is lacking. The aquatic immature stages of mayflies (Ephemeroptera) are widely used as bio-indicators and form an important component of Malagasy freshwater biodiversity. Many species are thought to be microendemics, restricted to single river basins in forested areas, making them particularly sensitive to habitat reduction and degradation. The Heptageniidae are a globally diverse family of mayflies (>500 species) but remain practically unknown in Madagascar except for two species described in 1996. The standard approach to understanding their diversity, endemism, and origin would require extensive field sampling on several continents and years of taxonomic work followed by phylogenetic analysis. Here we circumvent this using museum collections and freshly collected individuals in a combined approach of DNA taxonomy and phylogeny. The coalescent-based GMYC analysis of DNA barcode data (mitochondrial COI) revealed 14 putative species on Madagascar, 70% of which were microendemics. A phylogenetic analysis that included African and Asian species and data from two mitochondrial and four nuclear loci indicated the Malagasy Heptageniidae are monophyletic and sister to African species. The genus Compsoneuria is shown to be paraphyletic and the genus Notonurus is reinstalled for African and Malagasy species previously placed in Compsoneuria. A molecular clock excluded a Gondwanan vicariance origin and instead favoured a more recent overseas colonization of Madagascar. The observed monophyly and high microendemism highlight their conservation importance and suggest the DNA-based approach can rapidly provide information on the diversity, endemism, and origin of freshwater biodiversity. Our results underline the important role that museum collections can play in molecular studies, especially in critically endangered biodiversity hotspots like Madagascar where entire species or populations may go extinct very quickly.

The effect of geographical scale of sampling on DNA barcoding

Eight years after DNA barcoding was formally proposed on a large scale, CO1 sequences are rapidly accumulating from around the world. While studies to date have mostly targeted local or regional species assemblages, the recent launch of the global iBOL project (International Barcode of Life), highlights the need to understand the effects of geographical scale on Barcoding's goals. Sampling has been central in the debate on DNA Barcoding, but the effect of the geographical scale of sampling has not yet been thoroughly and explicitly tested with empirical data. Here, we present a CO1 data set of aquatic predaceous diving beetles of the tribe Agabini, sampled throughout Europe, and use it to investigate how the geographic scale of sampling affects 1) the estimated intraspecific variation of species, 2) the genetic distance to the most closely related heterospecific, 3) the ratio of intraspecific and interspecific variation, 4) the frequency of taxonomically recognized species found to be monophyletic, and 5) query identification performance based on 6 different species assignment methods. Intraspecific variation was significantly correlated with the geographical scale of sampling (R-square = 0.7), and more than half of the species with 10 or more sampled individuals (N = 29) showed higher intraspecific variation than 1% sequence divergence. In contrast, the distance to the closest heterospecific showed a significant decrease with increasing geographical scale of sampling. The average genetic distance dropped from > 7% for samples within 1 km, to < 3.5% for samples up to > 6000 km apart. Over a third of the species were not monophyletic, and the proportion increased through locally, nationally, regionally, and continentally restricted subsets of the data. The success of identifying queries decreased with increasing spatial scale of sampling; liberal methods declined from 100% to around 90%, whereas strict methods dropped to below 50% at continental scales. The proportion of query identifications considered uncertain (more than one species < 1% distance from query) escalated from zero at local, to 50% at continental scale. Finally, by resampling the most widely sampled species we show that even if samples are collected to maximize the geographical coverage, up to 70 individuals are required to sample 95% of intraspecific variation. The results show that the geographical scale of sampling has a critical impact on the global application of DNA barcoding. Scale-effects result from the relative importance of different processes determining the composition of regional species assemblages (dispersal and ecological assembly) and global clades (demography, speciation, and extinction). The incorporation of geographical information, where available, will be required to obtain identification rates at global scales equivalent to those in regional barcoding studies. Our result hence provides an impetus for both smarter barcoding tools and sprouting national barcoding initiatives-smaller geographical scales deliver higher accuracy.

Act together-implications of symbioses in aquatic ciliates

Mutual interactions in the form of symbioses can increase the fitness of organisms and provide them with the capacity to occupy new ecological niches. The formation of obligate symbioses allows for rapid evolution of new life forms including multitrophic consortia. Microbes are important components of many known endosymbioses and their short generation times and strong potential for genetic exchange may be important drivers of speciation. Hosts provide endo- and ectosymbionts with stable, nutrient-rich environments, and protection from grazers. This is of particular importance in aquatic ecosystems, which are often highly variable, harsh, and nutrient-deficient habitats. It is therefore not surprising that symbioses are widespread in both marine and freshwater environments. Symbioses in aquatic ciliates are good model systems for exploring symbiont-host interactions. Many ciliate species are globally distributed and have been intensively studied in the context of plastid evolution. Their relatively large cell size offers an ideal habitat for numerous microorganisms with different functional traits including commensalism and parasitism. Phagocytosis facilitates the formation of symbiotic relationships, particularly since some ingested microorganisms can escape the digestion. For example, photoautotrophic algae and methanogens represent endosymbionts that greatly extend the biogeochemical functions of their hosts. Consequently, symbiotic relationships between protists and prokaryotes are widespread and often result in new ecological functions of the symbiotic communities. This enables ciliates to thrive under a wide range of environmental conditions including ultraoligotrophic or anoxic habitats. We summarize the current understanding of this exciting research topic to identify the many areas in which knowledge is lacking and to stimulate future research by providing an overview on new methodologies and by formulating a number of emerging questions in this field.

A molecular approach to identifying the natural prey of the African creeping water bug Naucoris, a potential reservoir of Mycobacterium ulcerans

The extra-oral digestion of creeping water bugs (Naucoridae: Hemiptera) hinders the study of their diet using the standard method of identifying prey body parts in the gut. Genetic methods are available, but rely on PCR tests or similar diagnostics to confirm suspected prey. Where the potential prey is unknown and a broad search for all possible prey is desirable, methods that can potentially capture any prey item are required. Naucoris sp. is known to harbor Mycobacterium ulcerans (Actinomycetales: Mycobacteriaceae), the causative bacterium of Buruli ulcer. Outbreaks of Buruli ulcer have been associated with disturbed freshwater habitats, but the mode of transmission to humans remains unclear. Here we examine the diet of Naucoris sp., a dominant aquatic predator in water bodies in Ghana where the prevalence of Buruli ulcer is high. We cloned and sequenced 576 PCR products (mtDNA rrnL, cox1) isolated from the gut of 60 Naucoris sp. individuals to determining diet composition as broadly as possible. Using phylogenetic analysis of newly sequenced clones and 6 potential prey taxa collected from the site, sequences isolated from Naucoris sp. guts matched locally collected Coleoptera (Hydrophilidae). Blastn queries to GenBank of other clone sequences produced matches to (Anura) (n = 1), Rotifera (n = 5), and fungi (n = 4) as additional components of the diet. Our results suggest that sp. in this Buruli ulcer-endemic area feeds on a wide range of prey and body sizes, and that the approach could be successfully applied to studies of aquatic food webs where morphological identification of prey is impossible and where little or no a priori knowledge is available.

Toward a DNA taxonomy of Alpine Rhithrogena (Ephemeroptera: Heptageniidae) using a mixed Yule-coalescent analysis of mitochondrial and nuclear DNA

Aquatic larvae of many Rhithrogena mayflies (Ephemeroptera) inhabit sensitive Alpine environments. A number of species are on the IUCN Red List and many recognized species have restricted distributions and are of conservation interest. Despite their ecological and conservation importance, ambiguous morphological differences among closely related species suggest that the current taxonomy may not accurately reflect the evolutionary diversity of the group. Here we examined the species status of nearly 50% of European Rhithrogena diversity using a widespread sampling scheme of Alpine species that included 22 type localities, general mixed Yule-coalescent (GMYC) model analysis of one standard mtDNA marker and one newly developed nDNA marker, and morphological identification where possible. Using sequences from 533 individuals from 144 sampling localities, we observed significant clustering of the mitochondrial (cox1) marker into 31 GMYC species. Twenty-one of these could be identified based on the presence of topotypes (expertly identified specimens from the species' type locality) or unambiguous morphology. These results strongly suggest the presence of both cryptic diversity and taxonomic oversplitting in Rhithrogena. Significant clustering was not detected with protein-coding nuclear PEPCK, although nine GMYC species were congruent with well supported terminal clusters of nDNA. Lack of greater congruence in the two data sets may be the result of incomplete sorting of ancestral polymorphism. Bayesian phylogenetic analyses of both gene regions recovered four of the six recognized Rhithrogena species groups in our samples as monophyletic. Future development of more nuclear markers would facilitate multi-locus analysis of unresolved, closely related species pairs. The DNA taxonomy developed here lays the groundwork for a future revision of the important but cryptic Rhithrogena genus in Europe.

Evolutionary criteria outperform operational approaches in producing ecologically relevant fungal species inventories

Analyses of the structure and function of microbial communities are highly constrained by the diversity of organisms present within most environmental samples. A common approach is to rely almost entirely on DNA sequence data for estimates of microbial diversity, but to date there is no objective method of clustering sequences into groups that is grounded in evolutionary theory of what constitutes a biological lineage. The general mixed Yule-coalescent (GMYC) model uses a likelihood-based approach to distinguish population-level processes within lineages from processes associated with speciation and extinction, thus identifying a distinct point where extant lineages became independent. Using two independent surveys of DNA sequences associated with a group of ubiquitous plant-symbiotic fungi, we compared estimates of species richness derived using the GMYC model to those based on operational taxonomic units (OTUs) defined by fixed levels of sequence similarity. The model predicted lower species richness in these surveys than did traditional methods of sequence similarity. Here, we show for the first time that groups delineated by the GMYC model better explained variation in the distribution of fungi in relation to putative niche-based variables associated with host species identity, edaphic factors, and aspects of how the sampled ecosystems were managed. Our results suggest the coalescent-based GMYC model successfully groups environmental sequences of fungi into clusters that are ecologically more meaningful than more arbitrary approaches for estimating species richness.

Genetic population structure of sympatric and allopatric populations of Baltic ciscoes (Coregonus albula complex, Teleostei, Coregonidae)

Background

Teleost fishes of the Coregonidae are good model systems for studying postglacial evolution, adaptive radiation and ecological speciation. Of particular interest is whether the repeated occurrence of sympatric species pairs results from in-situ divergence from a single lineage or from multiple invasions of one or more different lineages. Here, we analysed the genetic structure of Baltic ciscoes (Coregonus albula complex), examining 271 individuals from 8 lakes in northern Germany using 1244 polymorphic AFLP loci. Six lakes had only one population of C. albula while the remaining two lakes had C. albula as well as a sympatric species (C. lucinensis or C. fontanae).

Results

AFLP demonstrated a significant population structure (Bayesian θ^B = 0.22). Lower differentiation between allopatric (θ^B = 0.028) than sympatric (0.063-0.083) populations contradicts the hypothesis of a sympatric origin of taxa, and there was little evidence for stocking or ongoing hybridization. Genome scans found only three loci that appeared to be under selection in both sympatric population pairs, suggesting a low probability of similar mechanisms of ecological segregation. However, removal of all non-neutral loci decreased the genetic distance between sympatric pairs, suggesting recent adaptive divergence at a few loci. Sympatric pairs in the two lakes were genetically distinct from the six other C. albula populations, suggesting introgression from another lineage may have influenced these two lakes. This was supported by an analysis of isolation-by-distance, where the drift-gene flow equilibrium observed among allopatric populations was disrupted when the sympatric pairs were included.

Conclusions

While the population genetic data alone can not unambiguously uncover the mode of speciation, our data indicate that multiple lineages may be responsible for the complex patterns typically observed in Coregonus. Relative differences within and among lakes raises the possibility that multiple lineages may be present in northern Germany, thus understanding the postglacial evolution and speciation in the C. albula complex requires a large-scale phylogenetic analysis of several potential founder lineages.

Diversity and distribution of Braconidae, a family of parasitoid wasps in the central European peatbogs of South Bohemia, Czech Republic

An ecological overview of seven years investigation of Braconidae, a family of parasitoid wasps (Hymenoptera: Ichneumonoidea) and a tyrpho-classification of parasitoids in peatbog areas of South Bohemia, Czech Republic are given. A total of 350 species were recorded in investigated sites, but only five tyrphobionts (1.4%) are proposed: Microchelonus basalis, Microchelonus koponeni, Coloneura ate, Coloneura danica and Myiocephalus niger. All of these species have a boreal-alpine distribution that, in Central Europe, is associated only with peatbogs. Tyrphophilous behaviour is seen in at least four (1.1%) species: Microchelonus pedator, Microchelonus subpedator, Microchelonus karadagi and Microchelonus gravenhorstii; however, a number of other braconids prefer peatbogs because they were more frequently encountered within, rather than outside, the bog habitat. The rest of the braconids (342 species, 97.5%) are tyrphoneutrals, many of them being eurytopic components of various habitats throughout their current ranges. Lists of tyrphobiontic braconids and a brief commentary on species composition, distributional picture of actual ranges, and parasitoid association to bog landscape are provided. Being true refugial habitats for populations in an ever-changing world, peatbogs play a significant role in harboring insect communities.