Benchmarking sample pooling for epigenomics of natural populations.

DNA methylation (DNAm) is a mechanism for rapid acclimation to environmental conditions. In natural systems, small effect sizes relative to noise necessitates large sampling efforts to detect differences. Large numbers of individually sequenced libraries are costly. Pooling DNA prior to library preparation may be an efficient way to reduce costs and increase sample size, yet there are to date no recommendations in ecological epigenetics research. We test whether pooled and individual libraries yield comparable DNAm signals in a natural system exposed to different pollution levels by generating whole-epigenome data from two invasive molluscs (Corbicula fluminea, Dreissena polymorpha) collected from polluted and unpolluted localities (Italy). DNA of the same individuals were used for pooled and individual epigenomic libraries and sequenced with equivalent resources per individual. We found that pooling effectively captures similar genome-wide and global methylation signals as individual libraries, highlighting that pooled libraries are representative of the global population signal. However, pooled libraries yielded orders of magnitude more data than individual libraries, which was a consequence of higher coverage. We would therefore recommend aiming for a high initial coverage of individual libraries (15×) in future studies. Consequently, we detected many more differentially methylated regions (DMRs) with the pooled libraries and a significantly lower statistical power for regions from individual libraries. Computationally pooled data from the individual libraries produced fewer DMRs and the overlap with wet-lab pooled DMRs was relatively low. We discuss possible causes for discrepancies, list benefits and drawbacks of pooling, and provide recommendations for future epigenomic studies.

The genome of the rayed Mediterranean limpet Patella caerulea (Linnaeus, 1758).

Patella caerulea (Linnaeus, 1758) is a mollusc limpet species of the class Gastropoda. Endemic to the Mediterranean Sea, it is considered a keystone species due to its primary role in structuring and regulating the ecological balance of tidal and subtidal habitats. It is currently being used as a bioindicator to assess the environmental quality of coastal marine waters and as a model species to understand adaptation to ocean acidification. Here, we provide a high-quality reference genome assembly and annotation for P. caerulea. We generated ∼30 Gb of Pacific Biosciences high-fidelity data from a single individual and provide a final 749.8 Mb assembly containing 62 contigs, including the mitochondrial genome (14,938 bp). With an N50 of 48.8 Mb and 98% of the assembly contained in the 18 largest contigs, this assembly is near chromosome-scale. Benchmarking Universal Single-Copy Orthologs scores were high (Mollusca, 87.8% complete; Metazoa, 97.2% complete) and similar to metrics observed for other chromosome-level Patella genomes, highlighting a possible bias in the Mollusca database for Patellids. We generated transcriptomic Illumina data from a second individual collected at the same locality and used it together with protein evidence to annotate the genome. A total of 23,938 protein-coding gene models were found. By comparing this annotation with other published Patella annotations, we found that the distribution and median values of exon and gene lengths was comparable with other Patella species despite different annotation approaches. The present high-quality P. caerulea reference genome, available on GenBank (BioProject: PRJNA1045377; assembly: GCA_036850965.1), is an important resource for future ecological and evolutionary studies.

Effects of anthropogenic stress on hosts and their microbiomes: Treated wastewater alters performance and gut microbiome of a key detritivore (Asellus aquaticus).

Human activity is a major driver of ecological and evolutionary change in wild populations and can have diverse effects on eukaryotic organisms as well as on environmental and host-associated microbial communities. Although host-microbiome interactions can be a major determinant of host fitness, few studies consider the joint responses of hosts and their microbiomes to anthropogenic changes. In freshwater ecosystems, wastewater is a widespread anthropogenic stressor that represents a multifarious environmental perturbation. Here, we experimentally tested the impact of treated wastewater on a keystone host (the freshwater isopod Asellus aquaticus) and its gut microbiome. We used a semi-natural flume experiment, in combination with 16S rRNA amplicon sequencing, to assess how different concentrations (0%, 30%, and 80%) of nonfiltered wastewater (i.e. with chemical toxicants, nutrients, organic particles, and microbes) versus ultrafiltered wastewater (i.e. only dissolved pollutants and nutrients) affected host survival, growth, and food consumption as well as mid- and hindgut bacterial community composition and diversity. Our results show that while host survival was not affected by the treatments, host growth increased and host feeding rate decreased with nonfiltered wastewater - potentially indicating that A. aquaticus fed on organic matter and microbes available in nonfiltered wastewater. Furthermore, even though the midgut microbiome (diversity and composition) was not affected by any of our treatments, nonfiltered wastewater influenced bacterial composition (but not diversity) in the hindgut. Ultrafiltered wastewater, on the other hand, affected both community composition and bacterial diversity in the hindgut, an effect that in our system differed between sexes. While the functional consequences of microbiome changes and their sex specificity are yet to be tested, our results indicate that different components of multifactorial stressors (i.e. different constituents of wastewater) can affect hosts and their microbiome in distinct (even opposing) manners and have a substantial impact on eco-evolutionary responses to anthropogenic stressors.