Epidemiology and molecular typing of multidrug-resistant bacteria in tertiary hospitals and nursing homes in Flanders, Belgium.

This study aimed to map MDRO carriage and potential transmission within and between three Flemish tertiary care hospitals and their neighbouring nursing homes. A cross-sectional MDRO prevalence survey was organized between October 2017 and February 2019. Perianal swabs were cultured for detection of MDRO. Determination of clonal relatedness based on wgMLST allelic profiles was performed. The prevalence of MDRO in Belgian hospitals and NHs is on the rise, compared to previous studies, and transmission in and between institutions is observed. These results re-emphasize the need for a healthcare network-wide infection prevention strategy in which WGS of MDRO strains can be supportive.

Epidemiology and molecular typing of multidrug-resistant bacteria in day care centres in Flanders, Belgium.

The global prevalence and spread of multidrug-resistant organisms (MDROs) represent an emerging public health threat. Day care centre (DCC) attendance is a risk factor for MDRO carriage in children and their environment. This study aimed to map the epidemiology of carriage and potential transmission of these organisms within 18 Flemish DDCs (Belgium). An MDRO prevalence survey was organised between November 2018 and February 2019 among children attending the centres. Selective chromogenic culture media were used for the detection of extended-spectrum beta-lactamase-producing Enterobacterales (ESBL-E), carbapenemase-producing Enterobacterales (CPE), and vancomycin-resistant Enterococci (VRE) in faecal swabs obtained from diapers or jars (n = 448). All isolated MDROs were subjected to resistance gene sequencing. A total of 71 of 448 samples (15.8%) yielded isolates of ESBL-E with a predominance of Escherichia coli (92.2% of ESBL-E) and ESBL resistance gene blaCTX-M-15 (50.7% of ESBL coding genes in E. coli). ESBL-E prevalence varied between DCCs, ranging from 0 to 50%. Transmission, based on the clonal relatedness of ESBL-E strains, was observed. CPE was identified in only one child carrying an E. coli with an OXA-244 gene. VRE was absent from all samples. The observed prevalence of ESBL-E in Flemish DCCs is high compared with previous studies, and our findings re-emphasise the need for rigorous hygiene measures within such centres to control the further spread of MDROs in the community.

Genetic characterization of ESBL-producing and ciprofloxacin-resistant Escherichia coli from Belgian broilers and pigs.

Background: The increasing number of infections caused by Escherichia coli resistant to clinically important antibiotics is a global concern for human and animal health. High overall levels of extended-spectrum beta-lactamase (ESBL)-producing and ciprofloxacin-resistant (ciproR) Escherichia coli in livestock are reported in Belgium. This cross-sectional study aimed to genotypically characterize and trace ESBL-and ciproR-E. coli of Belgian food-producing animals. Methods: A total of 798 fecal samples were collected in a stratified-random sampling design from Belgian broilers and sows. Consequently, 77 ESBL-E. coli and 84 ciproR-E. coli were sequenced using Illumina MiSeq. Minimum inhibitory concentration (MIC) for fluoroquinolones and cephalosporins were determined. Molecular in silico typing, resistance and virulence gene determination, and plasmid identification was performed. Scaffolds harboring ESBL or plasmid-mediated quinolone resistance (PMQR) genes were analyzed to detect mobile genetic elements (MGEs) and plasmid origins. Core genome allelic distances were used to determine genetic relationships among isolates. Results: A variety of E. coli sequence types (ST) (n = 63), resistance genes and virulence profiles was detected. ST10 was the most frequently encountered ST (8.1%, n = 13). The pandemic multidrug-resistant clone ST131 was not detected. Most farms harbored more than one ESBL type, with bla CTX-M-1 (41.6% of ESBL-E. coli) being the most prevalent and bla CTX M-15 (n = 3) being the least prevalent. PMQR genes (15.5%, n = 13) played a limited role in the occurrence of ciproR-E. coli. More importantly, sequential acquisition of mutations in quinolone resistance-determining regions (QRDR) of gyrA and parC led to increasing MICs for fluoroquinolones. GyrA S83L, D87N and ParC S80I mutations were strongly associated with high-level fluoroquinolone resistance. Genetically related isolates identified within the farms or among different farms highlight transmission of resistant E. coli or the presence of a common reservoir. IncI1-I(alpha) replicon type plasmids carried different ESBL genes (bla CTX-M-1, bla CTX-M-32 and bla TEM-52C). In addition, the detection of plasmid replicons with associated insertion sequence (IS) elements and ESBL/PMQR genes in different farms and among several STs (e.g., IncI1-I(alpha)/IncX3) underline that plasmid transmission could be another important contributor to transmission of resistance in these farms. Conclusion: Our findings reveal a multifaceted narrative of transmission pathways. These findings could be relevant in understanding and battling the problem of antibiotic resistance in farms.

Harmonization of whole-genome sequencing for outbreak surveillance of Enterobacteriaceae and Enterococci.

Whole-genome sequencing (WGS) is becoming the de facto standard for bacterial typing and outbreak surveillance of resistant bacterial pathogens. However, interoperability for WGS of bacterial outbreaks is poorly understood. We hypothesized that harmonization of WGS for outbreak surveillance is achievable through the use of identical protocols for both data generation and data analysis. A set of 30 bacterial isolates, comprising of various species belonging to the Enterobacteriaceae family and Enterococcus genera, were selected and sequenced using the same protocol on the Illumina MiSeq platform in each individual centre. All generated sequencing data were analysed by one centre using BioNumerics (6.7.3) for (i) genotyping origin of replications and antimicrobial resistance genes, (ii) core-genome multi-locus sequence typing (cgMLST) for Escherichia coli and Klebsiella pneumoniae and whole-genome multi-locus sequencing typing (wgMLST) for all species. Additionally, a split k-mer analysis was performed to determine the number of SNPs between samples. A precision of 99.0% and an accuracy of 99.2% was achieved for genotyping. Based on cgMLST, a discrepant allele was called only in 2/27 and 3/15 comparisons between two genomes, for E. coli and K. pneumoniae, respectively. Based on wgMLST, the number of discrepant alleles ranged from 0 to 7 (average 1.6). For SNPs, this ranged from 0 to 11 SNPs (average 3.4). Furthermore, we demonstrate that using different de novo assemblers to analyse the same dataset introduces up to 150 SNPs, which surpasses most thresholds for bacterial outbreaks. This shows the importance of harmonization of data-processing surveillance of bacterial outbreaks. In summary, multi-centre WGS for bacterial surveillance is achievable, but only if protocols are harmonized.

Whole-genome sequencing to explore nosocomial transmission and virulence in neonatal methicillin-susceptible Staphylococcus aureus bacteremia.

BACKGROUND: Neonatal Staphylococcus aureus (S. aureus) bacteremia is an important cause of morbidity and mortality. In this study, we examined whether methicillin-susceptible S. aureus (MSSA) transmission and genetic makeup contribute to the occurrence of neonatal S. aureus bacteremia. METHODS: A retrospective, single-centre study was performed. All patients were included who suffered from S. aureus bacteremia in the neonatal intensive care unit (NICU), Erasmus MC-Sophia, Rotterdam, the Netherlands, between January 2011 and November 2017. Whole-genome sequencing (WGS) was used to characterize the S. aureus isolates, as was also done in comparison to reference genomes. Transmission was considered likely in case of genetically indistinguishable S. aureus isolates. RESULTS: Excluding coagulase-negative staphylococci (CoNS), S. aureus was the most common cause of neonatal bacteremia. Twelve percent (n = 112) of all 926 positive blood cultures from neonates grew S. aureus. Based on core genome multilocus sequence typing (cgMLST), 12 clusters of genetically indistinguishable MSSA isolates were found, containing 33 isolates in total (2-4 isolates per cluster). In seven of these clusters, at least two of the identified MSSA isolates were collected within a time period of one month. Six virulence genes were present in 98-100% of all MSSA isolates. In comparison to S. aureus reference genomes, toxin genes encoding staphylococcal enterotoxin A (sea) and toxic shock syndrome toxin 1 (tsst-1) were present more often in the genomes of bacteremia isolates. CONCLUSION: Transmission of MSSA is a contributing factor to the occurrence of S. aureus bacteremia in neonates. Sea and tsst-1 might play a role in neonatal S. aureus bacteremia.

Massively parallel sequencing of micro-manipulated cells targeting a comprehensive panel of disease-causing genes: A comparative evaluation of upstream whole-genome amplification methods.

Single Gene Disorders (SGD) are still routinely diagnosed using PCR-based assays that need to be developed and validated for each individual disease-specific gene fragment. The TruSight One sequencing panel currently covers 12 Mb of genomic content, including 4813 genes associated with a clinical phenotype. When only a limited number of cells are available, whole genome amplification (WGA) is required prior to DNA target capture techniques such as the TruSight One panel. In this study, we compared 4 different WGA methods in combination with the TruSight One sequencing panel to perform single nucleotide polymorphism (SNP) genotyping starting from 3 micro-manipulated cells. This setting simulates clinical settings such as day-5 blastocyst biopsy for Preimplantation Genetic Testing (PGT), liquid biopsy of circulating tumor cells (CTCs) and cancer-cell profiling. Bulk cell samples were processed alongside these WGA samples to serve as a performance reference. Target coverage, coverage uniformity and SNP calling accuracy obtained using any of the WGA, is inferior to the results obtained on bulk cell samples. However, results after REPLI-g come close. Compared to the other WGA methods, the method using REPLI-g WGA results in a better coverage of the targeted genomic regions with a more uniform read depth. Consequently, this method also results in a more accurate SNP calling and could be considered for clinical genotyping of a limited number of cells.

Genome-Wide Copy Number Alteration Detection in Preimplantation Genetic Diagnosis.

Shallow whole genome sequencing has recently been introduced for genome-wide detection of chromosomal copy number alterations (CNAs) in preimplantation genetic diagnosis (PGD), using only 4-7 trophectoderm cells biopsied from day-5 embryos. This chapter describes the complete method, starting from whole genome amplification (WGA) on isolated blastomere(s), up to data analysis for CNA detection. The process is described generically and can also be used to perform CNA analysis on a limited number of cells (down to a single cell) in other applications. This unique description also includes some tips and tricks to increase the chance of success.

Microbial community profiling of fresh basil and pitfalls in taxonomic assignment of enterobacterial pathogenic species based upon 16S rRNA amplicon sequencing.

Application of 16S rRNA (gene) amplicon sequencing on food samples is increasingly applied for assessing microbial diversity but may as unintended advantage also enable simultaneous detection of any human pathogens without a priori definition. In the present study high-throughput next-generation sequencing (NGS) of the V1-V2-V3 regions of the 16S rRNA gene was applied to identify the bacteria present on fresh basil leaves. However, results were strongly impacted by variations in the bioinformatics analysis pipelines (MEGAN, SILVAngs, QIIME and MG-RAST), including the database choice (Greengenes, RDP and M5RNA) and the annotation algorithm (best hit, representative hit and lowest common ancestor). The use of pipelines with default parameters will lead to discrepancies. The estimate of microbial diversity of fresh basil using 16S rRNA (gene) amplicon sequencing is thus indicative but subject to biases. Salmonella enterica was detected at low frequencies, between 0.1% and 0.4% of bacterial sequences, corresponding with 37 to 166 reads. However, this result was dependent upon the pipeline used: Salmonella was detected by MEGAN, SILVAngs and MG-RAST, but not by QIIME. Confirmation of Salmonella sequences by real-time PCR was unsuccessful. It was shown that taxonomic resolution obtained from the short (500bp) sequence reads of the 16S rRNA gene containing the hypervariable regions V1-V3 cannot allow distinction of Salmonella with closely related enterobacterial species. In conclusion 16S amplicon sequencing, getting the status of standard method in microbial ecology studies of foods, needs expertise on both bioinformatics and microbiology for analysis of results. It is a powerful tool to estimate bacterial diversity but amenable to biases. Limitations concerning taxonomic resolution for some bacterial species or its inability to detect sub-dominant (pathogenic) species should be acknowledged in order to avoid overinterpretation of results.

Infections of virulent and avirulent viruses differentially influenced the expression of dicer-1, ago-1, and microRNAs in Bombus terrestris.

The microRNA (miRNA) pathway is well established to be involved in host-pathogen interactions. As key insect pollinators, bees are suffering from widely spreading viruses, especially honeybees and bumblebees. In order to better understand bee-virus interaction, we comparatively analyzed the involvement of the bumblebee miRNA pathway upon infection by two different viruses. In our setup, an avirulent infection is induced by slow bee paralysis virus (SBPV) and a virulent infection is induced by Israeli acute paralysis virus (IAPV). Our results showed the increased expressions of dicer-1 and ago-1 upon SBPV infection. There were 17 and 12 bumblebee miRNAs differentially expressed upon SBPV and IAPV infections, respectively. These results may indicate the involvement of the host miRNA pathway in bumblebee-virus interaction. However, silencing of dicer-1 did not influence the genome copy number of SBPV. Target prediction for these differentially expressed miRNAs showed their possible involvement in targeting viral genomic RNA and in the regulation of networks in bumblebee. Our study opens a new insight into bee-virus interaction meditated by host miRNAs.

Bisulfite Sequencing with Daphnia Highlights a Role for Epigenetics in Regulating Stress Response to Microcystis through Preferential Differential Methylation of Serine and Threonine Amino Acids.

Little is known about the influence that environmental stressors may have on genome-wide methylation patterns, and to what extent epigenetics may be involved in environmental stress response. Yet, studies of methylation patterns under stress could provide crucial insights on stress response and toxicity pathways. Here, we focus on genome-wide methylation patterns in the microcrustacean Daphnia magna, a model organism in ecotoxicology and risk assessment, exposed to the toxic cyanobacterium Microcystis aeruginosa. Bisulfite sequencing of exposed and control animals highlighted differential methylation patterns in Daphnia upon exposure to Microcystis primarily in exonic regions. These patterns are enriched for serine/threonine amino acid codons and genes related to protein synthesis, transport and degradation. Furthermore, we observed that genes with differential methylation corresponded well with genes susceptible to alternative splicing in response to Microcystis stress. Overall, our results suggest a complex mechanistic response in Daphnia characterized by interactions between DNA methylation and gene regulation mechanisms. These results underscore that DNA methylation is modulated by environmental stress and can also be an integral part of the toxicity response in our study species.

Performance of a TthPrimPol-based whole genome amplification kit for copy number alteration detection using massively parallel sequencing.

Starting from only a few cells, current whole genome amplification (WGA) methods provide enough DNA to perform massively parallel sequencing (MPS). Unfortunately, all current WGA methods introduce representation bias which limits detection of copy number aberrations (CNAs) smaller than 3 Mb. A recent WGA method, called TruePrime single cell WGA, uses a recently discovered DNA primase, TthPrimPol, instead of artificial primers to initiate DNA amplification. This method could lead to a lower representation bias, and consequently to a better detection of CNAs. The enzyme requires no complementarity and thus should generate random primers, equally distributed across the genome. The performance of TruePrime WGA was assessed for aneuploidy screening and CNA analysis after MPS, starting from 1, 3 or 5 cells. Although the method looks promising, the single cell TruePrime WGA kit v1 is not suited for high resolution CNA detection after MPS because too much representation bias is introduced.

Gene Body Methylation Patterns in Daphnia Are Associated with Gene Family Size.

The relation between gene body methylation and gene function remains elusive. Yet, our understanding of this relationship can contribute significant knowledge on how and why organisms target specific gene bodies for methylation. Here, we studied gene body methylation patterns in two Daphnia species. We observed both highly methylated genes and genes devoid of methylation in a background of low global methylation levels. A small but highly significant number of genes was highly methylated in both species. Remarkably, functional analyses indicate that variation in methylation within and between Daphnia species is primarily targeted to small gene families whereas large gene families tend to lack variation. The degree of sequence similarity could not explain the observed pattern. Furthermore, a significant negative correlation between gene family size and the degree of methylation suggests that gene body methylation may help regulate gene family expansion and functional diversification of gene families leading to phenotypic variation.

The Equine Embryo Influences Immune-Related Gene Expression in the Oviduct.

Although the equine oviduct clearly affects early embryo development and the selective transport of equine embryos through the oviduct indicates a reciprocal interaction, the influence of the embryo on gene expression in the oviduct remains to be determined in the horse. The aim of this study was to examine this by means of RNA sequencing. Four days after ovulation, epithelial cells ipsilateral and contralateral to the ovulation side from five cyclic and five pregnant mares were collected from the oviduct. RNA was extracted, samples were sequenced, and data analysis was performed to determine differentially expressed genes (DEGs) (P value ≤0.05 and absolute fold change ≥2) and to provide functional interpretation. A total of 10 743 transcripts were identified and 253 genes were found to be upregulated and 108 to be downregulated in the pregnant ipsilateral oviduct when compared to the cyclic ipsilateral oviduct. Comparison of the ipsilateral and the contralateral oviduct indicated 164 DEGs in pregnant mares and 77 DEGs in cyclic mares. Enriched functional categories were detected only in the comparison of pregnant and cyclic ipsilateral oviducts and showed that the equine embryo affects the expression of immune response-related genes in the oviduct, with marked upregulation of interferon-associated genes. This research represents the foundation for further assessment of the role of specific genes in the early embryo-maternal dialogue of the horse.

Suboptimal culture conditions induce more deviations in gene expression in male than female bovine blastocysts.

BACKGROUND: Since the development of in vitro embryo production in cattle, different supplements have been added to culture media to support embryo development, with serum being the most popular. However, the addition of serum during embryo culture can induce high birthweights and low viability in calves (Large Offspring Syndrome). Analysis of global gene expression in bovine embryos produced under different conditions can provide valuable information to optimize culture media for in vitro embryo production. RESULTS: We used RNA sequencing to examine the effect of in vitro embryo production, in either serum-containing or serum-free media, on the global gene expression pattern of individual bovine blastocysts. Compared to in vivo derived embryos, embryos produced in serum-containing medium had five times more differentially expressed genes than embryos produced in serum-free conditions (1109 vs. 207). Importantly, in vitro production in the presence of serum appeared to have a different impact on the embryos according to their sex, with male embryos having three times more genes differentially expressed than their female counterparts (1283 vs. 456). On the contrary, male and female embryos produced in serum-free conditions showed the same number (191 vs. 192) of genes expressed differentially; however, only 44 of those genes were common in both comparisons. The pathways affected by in vitro production differed depending on the type of supplementation. For example, embryos produced in serum-containing conditions had a lower expression of genes related to metabolism while embryos produced in serum-free conditions showed aberrations in genes involved in lipid metabolism. CONCLUSIONS: Serum supplementation had a major impact on the gene expression pattern of embryos, with male embryos being the most affected. The transcriptome of embryos produced in serum-free conditions showed a greater resemblance to that of in vivo derived embryos, although genes involved in lipid metabolism were altered. Male embryos appeared to be most affected by suboptimal in vitro culture, i.e. in the presence of serum.

In vivo study of Dicer-2-mediated immune response of the small interfering RNA pathway upon systemic infections of virulent and avirulent viruses in Bombus terrestris.

Recent studies suggest a potent role of the small interfering RNA (siRNA) pathway in the control of bee viruses and its usefulness to tackle these viral diseases. However, the involvement of the siRNA pathway in the defense against different bee viruses is still poorly understood. Therefore, in this report, we comprehensively analyzed the response of the siRNA pathway in bumblebees of Bombus terrestris to systemic infections of the virulent Israeli acute paralysis virus (IAPV) and the avirulent slow bee paralysis virus (SBPV). Our results showed that IAPV and SBPV infections induced the expression of Dicer-2. IAPV infections also triggered the production of predominantly 22 nt-long virus-derived siRNAs (vsiRNAs). Intriguingly, these 22 nt-long vsiRNAs showed a high proportion of antigenomic IAPV sequences. Conversely, these predominantly 22 nt-long vsiRNAs of SBPV were not detected in SBPV infected bees. Furthermore, an "RNAi-of-RNAi" experiment on Dicer-2 did not result in altered genome copy numbers of IAPV (n = 17-18) and also not of SBPV (n = 11-12). Based on these results, we can speculate about the importance of the siRNA pathway in bumblebees for the antiviral response. During infection of IAPV, this pathway is probably recruited but it might be insufficient to control viral infection in our experimental setup. The host can control SBPV infection, but aside from the induction of Dicer-2 expression, no further evidence of the antiviral activity of the siRNA pathway was observed. This report may also enhance the current understanding of the siRNA pathway in the innate immunity of non-model insects upon different viral infections.

Shallow whole genome sequencing is well suited for the detection of chromosomal aberrations in human blastocysts.

OBJECTIVE: To add evidence that massive parallel sequencing (MPS) is a valuable substitute for array comparative genomic hybridization (arrayCGH) with a resolution that is more appropriate for preimplantation genetic diagnosis (PGD) in translocation carriers. DESIGN: Study of diagnostic accuracy. SETTING: University hospital. PATIENT(S): Fifteen patients with a balanced structural rearrangement were included in the study: eight reciprocal translocations, four Robertsonian translocations, two inversions, and one insertional translocation. INTERVENTION(S): Trophectoderm biopsy was performed on 47 blastocysts. MAIN OUTCOME MEASURE(S): In the current study, shallow whole genome MPS on a NextSeq500 (Illumina) and Ion Proton (Life Technologies) instrument was performed in parallel on 47 whole genome amplified trophectoderm samples. Data analyses were performed using the QDNAseq algorithm implemented in Vivar. RESULT(S): In total, 5 normal and 42 abnormal embryos were analyzed. All aberrations previously detected with arrayCGH could be readily detected in the MPS data using both technologies and were correctly identified. The smallest detected abnormality was a ∼ 4.5 Mb deletion/duplication. CONCLUSION(S): This study demonstrates that shallow whole genome sequencing can be applied efficiently for the detection of numerical and structural chromosomal aberrations in embryos, equaling or even exceeding the resolution of the routinely used microarrays.

Characterization of the Bacterial Community Naturally Present on Commercially Grown Basil Leaves: Evaluation of Sample Preparation Prior to Culture-Independent Techniques.

Fresh herbs such as basil constitute an important food commodity worldwide. Basil provides considerable culinary and health benefits, but has also been implicated in foodborne illnesses. The naturally occurring bacterial community on basil leaves is currently unknown, so the epiphytic bacterial community was investigated using the culture-independent techniques denaturing gradient gel electrophoresis (DGGE) and next-generation sequencing (NGS). Sample preparation had a major influence on the results from DGGE and NGS: Novosphingobium was the dominant genus for three different basil batches obtained by maceration of basil leaves, while washing of the leaves yielded lower numbers but more variable dominant bacterial genera including Klebsiella, Pantoea, Flavobacterium, Sphingobacterium and Pseudomonas. During storage of basil, bacterial growth and shifts in the bacterial community were observed with DGGE and NGS. Spoilage was not associated with specific bacterial groups and presumably caused by physiological tissue deterioration and visual defects, rather than by bacterial growth.

Whole genome amplification with SurePlex results in better copy number alteration detection using sequencing data compared to the MALBAC method.

Current whole genome amplification (WGA) methods lead to amplification bias resulting in over- and under-represented regions in the genome. Nevertheless, certain WGA methods, such as SurePlex and subsequent arrayCGH analysis, make it possible to detect copy number alterations (CNAs) at a 10 Mb resolution. A more uniform WGA combined with massive parallel sequencing (MPS), however, could allow detection at higher resolution and lower cost. Recently, MALBAC, a new WGA method, claims unparalleled performance. Here, we compared the well-established SurePlex and MALBAC WGA for their ability to detect CNAs in MPS generated data and, in addition, compared PCR-free MPS library preparation with the standard enrichment PCR library preparation. Results showed that SurePlex amplification led to more uniformity across the genome, allowing for a better CNA detection with less false positives compared to MALBAC amplified samples. An even more uniform coverage was observed in samples following a PCR-free library preparation. In general, the combination of SurePlex and MPS led to the same chromosomal profile compared to a reference arrayCGH from unamplified genomic DNA, underlining the large potential of MPS techniques in CNA detection from a limited number of DNA material.

Conserved transcriptional responses to cyanobacterial stressors are mediated by alternate regulation of paralogous genes in Daphnia.

Despite a significant increase in genomic data, our knowledge of gene functions and their transcriptional responses to environmental stimuli remains limited. Here, we use the model keystone species Daphnia pulex to study environmental responses of genes in the context of their gene family history to better understand the relationship between genome structure and gene function in response to environmental stimuli. Daphnia were exposed to five different treatments, each consisting of a diet supplemented with one of five cyanobacterial species, and a control treatment consisting of a diet of only green algae. Differential gene expression profiles of Daphnia exposed to each of these five cyanobacterial species showed that genes with known functions are more likely to be shared by different expression profiles, whereas genes specific to the lineage of Daphnia are more likely to be unique to a given expression profile. Furthermore, while only a small number of nonlineage-specific genes were conserved across treatment type, there was a high degree of overlap in expression profiles at the functional level. The conservation of functional responses across the different cyanobacterial treatments can be attributed to the treatment-specific expression of different paralogous genes within the same gene family. Comparison with available gene expression data in the literature suggests differences in nutritional composition in diets with cyanobacterial species compared to diets of green algae as a primary driver for cyanobacterial effects on Daphnia. We conclude that conserved functional responses in Daphnia across different cyanobacterial treatments are mediated through alternate regulation of paralogous gene families.