Upscaling biodiversity monitoring: Metabarcoding estimates 31,846 insect species from Malaise traps across Germany.

Mitigating ongoing losses of insects and their key functions (e.g. pollination) requires tracking large-scale and long-term community changes. However, doing so has been hindered by the high diversity of insect species that requires prohibitively high investments of time, funding and taxonomic expertise when addressed with conventional tools. Here, we show that these concerns can be addressed through a comprehensive, scalable and cost-efficient DNA metabarcoding workflow. We use 1815 samples from 75 Malaise traps across Germany from 2019 and 2020 to demonstrate how metabarcoding can be incorporated into large-scale insect monitoring networks for less than 50 € per sample, including supplies, labour and maintenance. We validated the detected species using two publicly available databases (GBOL and GBIF) and the judgement of taxonomic experts. With an average of 1.4 M sequence reads per sample we uncovered 10,803 validated insect species, of which 83.9% were represented by a single Operational Taxonomic Unit (OTU). We estimated another 21,043 plausible species, which we argue either lack a reference barcode or are undescribed. The total of 31,846 species is similar to the number of insect species known for Germany (~35,500). Because Malaise traps capture only a subset of insects, our approach identified many species likely unknown from Germany or new to science. Our reproducible workflow (~80% OTU-similarity among years) provides a blueprint for large-scale biodiversity monitoring of insects and other biodiversity components in near real time.

Rapid growth and the evolution of complete metamorphosis in insects.

More than 50% of all animal species are insects that undergo complete metamorphosis. The key innovation of these holometabolous insects is a pupal stage between the larva and adult when most structures are completely rebuilt. Why this extreme lifestyle evolved is unclear. Here, we test the hypothesis that a trade-off between growth and differentiation explains the evolution of this novelty. Using a comparative approach, we find that holometabolous insects grow much faster than hemimetabolous insects. Using a theoretical model, we then show how holometaboly evolves under a growth-differentiation trade-off and identify conditions under which such temporal decoupling of growth and differentiation is favored. Our work supports the notion that the holometabolous life history evolved to remove developmental constraints on fast growth, primarily under high mortality.

The haplotype-resolved Prymnesium parvum (type B) microalga genome reveals the genetic basis of its fish-killing toxins.

The catastrophic loss of aquatic life in the Central European Oder River in 2022, caused by a toxic bloom of the haptophyte microalga Prymnesium parvum (in a wide sense, s.l.), underscores the need to improve our understanding of the genomic basis of the toxin. Previous morphological, phylogenetic, and genomic studies have revealed cryptic diversity within P. parvum s.l. and uncovered three clade-specific (types A, B, and C) prymnesin toxins. Here, we used state-of-the-art long-read sequencing and assembled the first haplotype-resolved diploid genome of a P. parvum type B from the strain responsible for the Oder disaster. Comparative analyses with type A genomes uncovered a genome-size expansion driven by repetitive elements in type B. We also found conserved synteny but divergent evolution in several polyketide synthase (PKS) genes, which are known to underlie toxin production in combination with environmental cues. We identified an approximately 20-kbp deletion in the largest PKS gene of type B that we link to differences in the chemical structure of types A and B prymnesins. Flow cytometry and electron microscopy analyses confirmed diploidy in the Oder River strain and revealed differences to closely related strains in both ploidy and morphology. Our results provide unprecedented resolution of strain diversity in P. parvum s.l. and a better understanding of the genomic basis of toxin variability in haptophytes. The reference-quality genome will enable us to better understand changes in microbial diversity in the face of increasing environmental pressures and provides a basis for strain-level monitoring of invasive Prymnesium in the future.

Unpredicted ecosystem response to compound human impacts in a European river.

Climate change elevates the threat of compound heat and drought events, with their ecological and socioeconomic impacts exacerbated by human ecosystem alterations such as eutrophication, salinization, and river engineering. Here, we study how multiple stressors produced an environmental disaster in a large European river, the Oder River, where a toxic bloom of the brackish-water planktonic haptophyte Prymnesium parvum (the "golden algae") killed approximately 1000 metric tons of fish and most mussels and snails. We uncovered the complexity of this event using hydroclimatic data, remote sensing, cell counts, hydrochemical and toxin analyses, and genetics. After incubation in impounded upstream channels with drastically elevated concentrations of salts and nutrients, only a critical combination of chronic salt and nutrient pollution, acute high water temperatures, and low river discharge during a heatwave enabled the riverine mass proliferation of B-type P. parvum along a 500 km river section. The dramatic losses of large filter feeders and the spreading of vegetative cells and resting stages make the system more susceptible to new harmful algal blooms. Our findings show that global warming, water use intensification, and chronic ecosystem pollution could increase likelihood and severity of such compound ecoclimatic events, necessitating consideration in future impact models.

Phylogenomics including new sequence data of phytoplankton-infecting chytrids reveals multiple independent lifestyle transitions across the phylum.

Parasitism is the most common lifestyle on Earth and has emerged many times independently across the eukaryotic tree of life. It is frequently found among chytrids (Chytridiomycota), which are early-branching unicellular fungi that feed osmotrophically via rhizoids as saprotrophs or parasites. Chytrids are abundant in most aquatic and terrestrial environments and fulfil important ecosystem functions. As parasites, they can have significant impacts on host populations. They cause global amphibian declines and influence the Earth's carbon cycle by terminating algal blooms. To date, the evolution of parasitism within the chytrid phylum remains unclear due to the low phylogenetic resolution of rRNA genes for the early diversification of fungi, and because few parasitic lineages have been cultured and genomic data for parasites is scarce. Here, we combine transcriptomics, culture-independent single-cell genomics and a phylogenomic approach to overcome these limitations. We newly sequenced 29 parasitic taxa and combined these with existing data to provide a robust backbone topology for the diversification of Chytridiomycota. Our analyses reveal multiple independent lifestyle transitions between parasitism and saprotrophy among chytrids and multiple host shifts by parasites. Based on these results and the parasitic lifestyle of other early-branching holomycotan lineages, we hypothesise that the chytrid last common ancestor was a parasite of phytoplankton.

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.

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,6 having diverged from the animals ∼1.35 billion years ago.7 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 lineage10 have brought the number of proposed phyla to 20,11 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.

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 CO2, 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, blaOX58) were found in 85% of sites, and were also more prevalent in urban streams and rivers. Light pollution and CO2 were significantly related to microbial community composition, with CO2 inversely related to microbial phototrophy. Results of synchronous nationwide sampling indicate that pollution-driven alterations to the freshwater microbiome lead to changes in CO2 production in natural waters and highlight the vulnerability of running waters to anthropogenic stressors.

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.

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 FST 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.

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.

Temporal dynamics of freshwater planktonic parasites inferred using a DNA metabarcoding time-series.

Parasites are important components of biodiversity and contributors to ecosystem functioning, but are often neglected in ecological studies. Most studies examine model parasite systems or single taxa, thus our understanding of community composition is lacking. Here, the seasonal and annual dynamics of parasites was quantified using a 5-year metabarcoding time-series of freshwater plankton, collected weekly. We first identified parasites in the dataset using literature searches of the taxonomic match and using sequence metadata from the National Center for Biotechnology Information (NCBI) nucleotide database. In total, 441 amplicon sequence variants (belonging to 18 phyla/clades) were classified as parasites. The four phyla/clades with the highest relative read abundance and richness were Chytridiomycota, Dinoflagellata, Oomycota and Perkinsozoa. Relative read abundance of total parasite taxa, Dinoflagellata and Perkinsozoa significantly varied with season and was highest in summer. Parasite richness varied significantly with season and year, and was generally lowest in spring. Each season had distinct parasite communities, and the difference between summer and winter communities was most pronounced. Combining DNA metabarcoding with searches of the literature and NCBI metadata allowed us to characterize parasite diversity and community dynamics and revealed the extent to which parasites contribute to the diversity of freshwater plankton communities.

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.

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.