Elysia timida was collected on Giglio, Italy (42°22′ N, 10°52′ E and 42°21′ N, 10°52′ E) between 3 and 6 m depth.
Elysia cornigera was collected at the Florida Keys (24°38′ N, 81°18′ W) at about 1 m depth. Both species were kept at a 12 L : 12 D rhythm at 25 µmol quanta m
−2 s
−1 in 3.7% salt water (Tropic Marine) and at 21°C.
Acetabularia acetabulum DI1 (isolate of Professor Diedrik Menzel, Bonn University, Germany) was grown at the same conditions as the animals, but in algal f/2 medium. For starvation experiments, animals were moved to Petri dishes (∅︀200 mm) and 50% of the water was changed every other day. Monolinuron (JBL GmbH) was added from a stock solution of 4000 mg l
−1 to a final concentration of 2 µg ml
−1. Bleaching was performed by illuminating the slugs with 1000 µmol quanta m
−2 s
−1 for 1 h every day at the same time. For each pigment extraction, pooled slugs were transferred to 90% acetone for each time point, homogenized and kept at −20°C for 1 day. Extractions for every time point were carried out as biological triplicates, each replicate consisting of at least three
E. cornigera and two
E. timida. Crude extracts were filtered through a 200 nM polytetrafluoroethylene membrane and then analysed by reversed-phases high pressure liquid chromatography with ultraviolet/visible spectroscopy detection (Hitachi/Merk) as described earlier [32 (
link)]. Pigment concentrations were determined using external pigment standards isolated from spinach thylakoids. Dry weight was determined on dried slug homogenate after extraction.
Photosystem II maximum quantum yield (
Fv/
Fm) was determined using a WALZ
MINI-PAM as described previously [33 (
link)], on the same group of individuals that were used for
14CO
2 incorporation measurements. For every time point for which
14CO
2 incorporation was determined, an individual group of slugs (at least 12
E. cornigera or nine
E. timida) was kept in isolated dishes. For each of these groups, separate
Fv/
Fm measurements were performed, and each individual slug of a group was measured at least three times. For each data point, the mean of all individual triplicate measurements was calculated. The light-driven incorporation of
14CO
2 was determined as described earlier [5 (
link)]. Briefly, for each triplicate measurement four individuals of
E. cornigera and three of
E. timida were used. Slugs were incubated in 1.2 ml artificial seawater supplemented with 0.40 mM [
14C]-NaHCO
2 (25 µCi per incubation, NEN-radiochemicals, MA, USA) for 2 h at room temperature and 72 µmol quanta m
−2 s
−1 illumination. After rinsing, homogenization and acidification of the material with 150 µl 1 M HCl, all the substrate was expelled from the homogenate. Incorporation of
14C by the slugs was determined in a scintillation counter after the addition of 12 ml LUMA-Gel scintillation cocktail (LUMAC, The Netherlands).
For ROS-imaging, slugs were incubated for 45 min with 100 µM DHE (DHE; excitation/emission HeNe 543/610 nm; Sigma) or 100 µM DCF (2′,7′-dichlorofluorescin; excitation/emission Argon 488/529 nm) plus 2 µM MitoTracker Red CMXRos (excitation/emission HeNe 543/599 nm; LifeTechnologies) in artificial seawater. Slugs were rinsed twice with seawater and decapitated before mounting. Confocal laser scanning microscopy was carried out using a Zeiss LSM 710. Imaging was always performed with the same settings and at a similar position along the parapodial rim at the same depth relative to the epidermis. Images were processed using Fiji/I
mageJ 1.48f [34 (
link)] and statistics were performed using R [35 ]. Normality was tested via a Shapiro–Wilk test [36 (
link)] and significance using a Mann–Whitney
U-test [37 (
link)].
Total RNA was extracted three times from seven
E. cornigera (a total of 21) and three times from five
E. timida (a total of 15) feeding on
A. acetabulum (
t0, feeding control group), 28 individuals of
E. cornigera that were starved for 4 days, of which nine individuals had only starved (S), nine had starved and were treated with monolinuron (S + M), and 10 of which had starved and were treated with 1000 µmol quanta m
−2 s
−1 for 1 h each day (S + B). The RNA of 20
E. timida's (7 S, 7 S + M, 6 S + B) was isolated after four days of starvation (
t1), of 30
E. cornigera (10 S, 9 S + M, 11 S + B) and of 20
E. timida (6 S, 8 S + M, 6 S + B) after 7 days of starvation (
t2), and finally of 32
E. timida (11 S, 10 S + M, 11 S + B) after 30 days of starvation (
t3). RNA was isolated using TRIzol (Life Technologies) according to the manufacturer's protocol and an additional DNAse treatment (Thermo Fisher Scientific). Poly(A) mRNA enrichment, library preparation using the TruSeq kit (Illumina) and 100 bp paired end sequencing using the Illumina HiSeq2000 system was performed by the Beijing Genome Institute (Hong Kong). Reads with remaining adapter sequences, more than 5% of unknown nucleotides or more than 20% of nucleotides with quality scores less than 10 (Illumina 1.5) were removed as part of the raw data extraction.
A total of 1 186 405 486 reads (a minimum of 52 million reads/library) were obtained (electronic supplementary material, table S2). Reads were inspected by F
astQC v. 0.10.1 [38 ] and filtered and trimmed using T
rimmomatic 0.32 [39 (
link)] (parameters: ILLUMINACLIP:TruSeq3-PE.fa:2 : 30 : 10; HEADCROP:10; TRAILING:3; SLIDINGWINDOW:4 : 20; MINLEN:36). Reads of all samples were assembled using T
rinity r20131110 [40 (
link)], which yielded 249 855 and 150 314 transcripts for
E. cornigera and
E. timida, respectively. Expression values for transcript clusters (or unigenes) were extracted and analysed using the T
rinity pipeline [41 (
link)] via RSEM 1.2.11 [42 (
link)], B
owtie 1.0.1 [43 (
link)] and
edgeR 3.4.2 [44 (
link)]. Only unigenes with raw read counts of greater than or equal to 100 over at least two samples were analysed, resulting in 36 380 and 32 897 unigenes, respectively. Differential expression was determined by log
2fc of at least ±1 compared to the control set of unigenes (with raw read counts of greater than or equal to 100). Expression change significance was assessed based upon the dispersion of available replicate information in
edgeR. The longest isoform of each cluster was defined as the representative sequence for a unigene.
To plot the overall taxonomic distribution of the assembled transcripts, unigenes were subjected to a BLASTx-based [45 (
link)] search (
e-value cut-off 10
−10) against a database consisting of protein sequences of R
efS
eq version 64 [46 (
link)] plus those of the genomes of
Crassostrea gigas [47 (
link)] and
Lottia gigantean [48 (
link)] (electronic supplementary material, figure S2). For all downstream analyses, only those genes were included for which top hits to the mentioned organisms or other metazoans were retrieved. Excluded were those with best blast hits to plants, protists, prokaryotes and viruses. Protein annotations for the 14 848
E. cornigera and 13 875
E. timida unigenes were extracted from mentioned BLAST hits and from a second BLASTx search to the UniProtKB/Swiss-Prot database [49 (
link)]. KOG/NOG (Eu
Karyotic/
Non-supervised
Orthologous
Groups [50 (
link)–52 (
link)]) categories were determined based on the best BLASTx hit to protein sequences of
eggNOG v. 4 [53 (
link)]. KEGG [54 (
link)] accessions were obtained using the KAAS 1.6a webservice [55 (
link)] against metazoa.
For qRT-PCR, 200 ng RNA was reverse transcribed using random hexamers and the iScript Select cDNA Synthesis Kit (Bio-Rad). Two-step qRT-PCR was performed on biological triplicates for each time point and treatment (each containing pooled RNA from greater than or equal to five slugs) using the SsoAdvanced Universal SYBR Green Supermix (Bio-Rad) according to the manufacturer's instructions on a StepOne Plus (Applied Biosystems) real-time PCR system. Primers were designed using Primer-BLAST [56 (
link)] (for primer sequences, see the electronic supplementary material, table S4) and data were analysed according to Pfaffl [57 (
link)], using the
EcRPL38 (Eco000149),
EtRPL19 (Eti000121) and
EtSETMAR (Eti000317) as reference genes.
de Vries J., Woehle C., Christa G., Wägele H., Tielens A.G., Jahns P, & Gould S.B. (2015). Comparison of sister species identifies factors underpinning plastid compatibility in green sea slugs. Proceedings of the Royal Society B: Biological Sciences, 282(1802), 20142519.