bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
1 Feeding intensity and molecular prey identification of
2 the common long-armed octopus, Octopus minor
3 (Mollusca: Octopodidae) in the wild
4
1, 2 1, 3* 1 5 Qi-Kang Bo , Xiao-Dong Zheng , Zhi-Wei Chen
6
7 1 Institute of Evolution and Marine Biodiversity, Ocean University of China
8 2 Breeding station of Bohai Sea aquatic resources, Tianjin Fisheries Research Institute
9 3 Key Laboratory of Mariculture, Ministry of Education, Ocean University of China
10
11 *Corresponding Author
12 E-mail: [email protected] (X-D Z)
1 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
14
15 Abstract
16 The common long-armed octopus, Octopus minor, is an important component
17 of systems and supports the local fisheries in the coastal areas of northern China.
18 And because of the overfishing, the national germplasm reserve was established as
19 conservation area at the Moon Lake for the genetic resource of O.minor in 2012. For
20 the fishery management and artificial breeding, especially for the management of
21 exclusive conservation reserves, its role in the ecosystem requires assessment.
22 Therefore, the feeding intensity of O. minor was studied from April to July 2014
23 when females reaching maturation, and prey composition was identified from
24 stomach contents using a DNA barcoding method. Of the 172 sampled octopuses, 66
25 had stomach contents that were nearly digested into pulp. Maximum feeding
26 intensity occurred during the month of April and the feeding intensity of the females
27 was greater than that of the males in April and May. A considerable overall
28 reduction of feeding intensity in both sexes occurred from April to July. A total of 8
29 species were identified as the prey of O. minor. Based on homology searches and
30 phylogenetic analysis, of 60 sequences, 30 matched with fish (50.00%, by number),
31 13 with crustaceans (21.66%), one with annelid (1.66%), one with nematode (1.66%)
32 and 15 with itself (25.00%). These results confirm that O. minor has habitual nature
33 of strong dietary preference, with Gobiidae families (62.79%, by number) being an
34 important prey during the time when females reach sexual maturation. From April to
35 July, the observed cannibalism showed an increasing trend.
2 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
36 Introduction
37 The common long-armed octopus, Octopus minor (Sasaki 1920) is a benthic and
38 neritic octopod that has become a commercially important species in the north of
39 China and in South Korea [1, 2] and is identified to be furthermore ecologically
40 important as a generalist predator. Because of overfishing, a national germplasm
41 reserve was established at the Moon Lake (Fig 1) in 2012 to serve as conservation for
42 the genetic resource of O. minor. Moon Lake is a shallow (depth <3 m) lagoon with a
43 muddy, bivalve-covered bottom, and patches of sea grass growing in the subtidal zone,
44 which provides a good habitat for O. minor. Additionally, a program of artificial
45 propagation and release of O. minor has been launched in the north of China [2].
46
47 Fig 1. The location of the sampling site in Moon Lake and the dots show the
48 sampling site.
49 The five methods used to assess the natural feeding habits of marine animals
50 include the examination of prey remains from middens [3], direct observation [4],
51 stomach content analysis [5, 6], isotopic assessment [7, 8], trophic tracing [9] and
52 molecular prey identification [10-12]. In some studies, several different methods have
53 been used in combination [11, 13, 14].
54 A range of factors has made it difficult to determine the natural diet of O. minor;
55 the primary challenge that the oesophageal diameter of the octopus is physically
56 limited as it passes through the brain; therefore, the octopus’ beak bites small pieces
57 of tissue to swallow, avoiding the ingestion of hard skeletal material [11], and it is
3 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
58 difficult to make an accurate identification of pray with lack of hard skeletal material.
59 Moreover, it is difficult to be collected on a large scale. O. minor burrow deep
60 (0.3-0.6 meter) interconnecting tunnels as nest in muddy marine bottom, leaving only
61 digging-holes and breathing-holes on the surface, into which it hides itself [15]; their
62 prey remains are left underground. Even if the prey remains were pushed out of the
63 nests, the light material would be easily removed by biotic and abiotic factors, while
64 the heavier material would be buried under the mud. Rapid digestion rates [16] and
65 external ingestion [17] make the stomach contents visually unidentifiable. These
66 specialized feeding strategies tend to bias data on prey species when morphological
67 analysis is used.
68 Food is an important factor to the extent that it governs growth, fecundity and
69 migratory movements. An understanding of the relationship between octopus species
70 and their favorite food items helps to locate potential feeding grounds, which may, in
71 turn, be helpful for the exploitation of these resources. The extent of the variation in
72 feeding must be taken into account when studying the diet of a generalist predator.
73 Octopods are opportunistic and their diets are often affected by the abundance of
74 available prey [3, 18]. The species, size and type of their diet in the wild are often
75 dominated by many factors such as maturity stage, seasonal variation [19], the body
76 size which reflects the prey ability [3, 14] and benthic assemblages which determine
77 prey availability [20]. In addition, what should be taken into account is that octopods
78 exhibit strong dietary preferences when given the same opportunity to different diets
79 items [18, 21-24].
4 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
80 Few studies have focused on the ecology of O. minor, whose diet is roughly
81 known to include crustaceans, molluscs, polychaetes and fish [25, 26]. Investigation
82 into the natural diet of this species during sexual maturation, are of more research
83 value, and can provide trophic relationship information for fishery management,
84 especially for the management of exclusive conservation reserves. It also provides a
85 pray reference for aquaculture and artificial breeding of this commercially important
86 specie. In this study, we applied the method of molecular prey identification, DNA
87 barcoding method, to identify natural prey of O. minor based on stomach contents.
88 Materials and methods
89 Specimen collection
90 All the analyses have been carried out using freshly dead specimens collected
91 from fishermen of the Mashan Group Co. Ltd. No use of live animals has been
92 required for this study and no specific permissions were needed for the sampling
93 activities in all of the investigated areas because our species of interest is
94 commercially harvested (not endangered nor protected) and it was caught in areas
95 where fishing is allowed. One hundred and seventy-two O. minor specimens were
96 captured using traps placed during the evenings and collected at dawn from Moon
97 Lake (122º35’E, 37º21’N, Fig 1) from April to July 2014. The same number of males
98 and females were collected each month. All specimens were brought to the laboratory
99 where their stomach contents were removed. Each stomach was opened and the
100 contents were flushed into cryogenic vials. The stomach contents, potential residue of
101 prey, were weighed and then preserved in 70% ethanol at -20℃ for later DNA
5 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
102 analysis [27]. Stomach content samples were numbered monthly, AP01 to AP19,
103 MA01 to MA13, JN01 to JN19 and JL01 to JL15 from April to July, respectively. To
104 avoid potential contaminants (e.g., blood and tissue attached to the stomach from the
105 predator), the exterior surface of each stomach was washed with sterile, distilled water
106 before removing the stomach contents [28]. The dorsal mantle length (ML), the
107 distance between the posterior midpoint of the mantle and the midpoint of the eyes,
108 and wet body weight (W) were also measured. This basic information and data were
109 showed in Table 1. Feeding intensity during the study months was determined based
110 on the degree of fullness of the stomach. The condition of feeding activity was
111 determined from observations of the degree of stomach distension as described by
112 Pillay (1952) [29]. Octopuses with stomachs that were gorged, full, ¾full, ½full were
113 considered to have been actively feeding (AF), while stomach ¼full, little and empty
114 were considered to denote poor feeding activity (PAF). The percentage of octopus in
115 the AF and PAF condition in both sexes for each month was calculated.
116
117 Table 1. The sample number, stomach content specimens number (Num.),
118 mantle length (ML, ±SD) and total weight (W, ±SD) of O. minor specimens
119 collected each month from Moon Lake
Sample Sample Female Male Num. date number ML (cm) W (g) ML (cm) W(g)
20/4/2014 52 8.68±0.95 120.47±24.99 8.60±0.89 135.65±28.26 19
14/5/2014 40 8.95±1.86 117.29±37.42 8.87±0.71 128.68±25.71 13
6 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
20/6/2014 40 10.55±1.14 192.46±40.28 10.30±1.09 202.88±42.66 19
16/7/2014 40 9.99±1.23 172.95±37.96 9.57±0.89 200.35±36.80 15
120
121
122 DNA extraction and sequence acquisition
123 Stomach contents were evenly ground in a homogenizer. Duplicate 20-40 mg
124 samples of mixed tissue items were subsampled from each stomach content and
125 genomic DNA was then isolated by standard phenol-chloroform purification
126 procedures. A fragment of the mitochondrial cytochrome oxidase Ⅰ (COⅠ) gene
127 was amplified using the universal primers [30].
128 Each polymerase chain reaction (PCR) was carried out in 50-µL volumes
129 containing 2U Taq DNA polymerase (Takara Co.), about 60 ng template DNA, 0.2
130 mM dNTPs, 0.25 µM of each primer, 2 mM MgCl2 and 1×PCR buffer. The PCR
131 amplification was performed on a GeneAmp® 9700 PCR System (Applied
132 Biosystems). Cycling conditions consisted of an initial denaturation at 94℃ for 3 min,
133 followed by 35 cycles of: denaturation at 94℃ for 1 min, annealing at 50℃ for 1 min,
134 extension at 72℃ for 1 min, and a final step of 5 min at 72℃.
135 Amplification products were confirmed by 1.5% TBE agarose gel
136 electrophoresis stained with ethidium bromide. The cleaned product was prepared for
137 sequencing using the BigDye Terminator Cycle Sequencing Kit (ver.3.1, Applied
138 Biosystems) and sequenced bidirectionally using an ABI PRISM 3730 (Applied
139 Biosystems) automatic sequencer. PCR products producing multiple bands indicated
7 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
140 that more than one prey species were present. Those PCR products were cloned using
141 the TOPO TA Cloning Kit (Invitrogen). Eight colonies per sample were selected for
142 colony PCR amplification and sequencing using the primers M13 (forward):
143 GTAAAACGACGGCCAG, andM13 (reverse): CAGGAAACAGCTATGAC.
144 DNA analyses and statistical analysis
145 Sequences of COⅠ were assembled and edited separately using DNASTAR
146 software (DNASTAR, Inc.), and then aligned with CLUSTAL_X 1.81 using the
147 default settings [31]. Sequences were considered to be part of the same ‘‘operational
148 taxonomic unit’’ (OTU), if there was less than a 1% sequence divergence, allowing
149 for the intraspecific variation and Taq polymerase errors [11, 13].
150 All of the obtained sequences were identified using the Identification System
151 (IDS) in the Barcode of Life Database (BOLD, www.boldsystems.org) and the Basic
152 Local Alignment Search Tool (BLAST) query algorithm in GenBank to establish
153 whenever possible the identification of the ingested material. The maximum
154 likelihood (ML) method was chosen to infer evolutionary history. Bootstrap
155 probabilities with 1,000 replications were calculated to assess reliability on each node
156 of the ML tree. Sequence divergence calculation and evolutionary analyses were
157 conducted in MEGA 6 software based on the Kimura 2-parameter model (K2P) [32].
158 The ML tree contained all of the sequences obtained from the stomach contents,
159 together with the closest matches to each sequence that were downloaded from BOLD
160 databases and GenBank. The criteria to assign identification to at the species level
161 required the sequence similarity display >98% in the BOLD database and BLAST [11]
8 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
162 and, if not, identification was restricted to the highest taxonomic lineage supported by
163 bootstrap probabilities higher than 70% in the consensus tree [14, 15]. The AF and
164 PAF differences the between sexes were verified by the Chi-squared fit test (2).
165 Statistical analyses were performed using IBM SPSS statistics version 20 (IBM,
166 Chicago, IL, USA). To corroborate that the number of analyzed stomachs was
167 adequate for diet description, a cumulative prey curve was generated using the
168 Estimate S Version 8.2 based on the prey identified [33]. The number of samples was
169 assumed to be sufficient to describe the diet when the curve approaches the
170 asymptote.
171
172 Results
173 Feeding intensity
174 Of the 172 octopus individuals, 66 had contents in their stomachs (Table 1) and
175 feeding intensity of O. minor varied to an extent in respect to the months, and more
176 than half of the stomachs were found to be empty (Table 2). There was significant
177 difference in the feeding activity between July and other months (p<0.05) and no
178 difference among the first three months (p>0.05). However, the AF percentage of
179 octopus from April is greater than that of other months, meaning the maximum
180 feeding intensity occurred during the month of April. A considerable overall reduction
181 in the intensity of feeding for both sexes occurred from April to July (Table 2). The
182 feeding intensity of the females is greater than that of the males based on the
183 percentage of AF and PAF between the sexes in each month, and in April and May
9 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
184 the feeding intensity of the females was significantly greater than that of the males
185 (Table 3). That is, compared to the males, more female individuals were actively
186 feeding.
187
188 Table 2. The distribution of stomach with different degree of distension and the
189 percentage of octopus with active feeding activity (AF) from April to July. Data
190 in the same column having different superscripted letters are significantly
191 different (p<0.05).
Stomach status (%) Feeding activity (%)
Month Gorge Full ¾Full ½Full ¼Full Little Empty AF d
Apr. 5.77 5.77 5.77 13.46 1.92 3.85 63.46 30.81a
May. 6.00 16.00 5.00 2.50 3.00 0.00 67.50 29.50a
Jun. 5.00 7.50 7.50 10.00 17.50 0.00 52.50 30.00a
Jul. 0.00 0.00 5.00 5.00 17.50 10.00 62.50 10.00b
192
193 Table 3. The percentage of octopus with active feeding activity (AF) and poor
194 feeding activity (PAF) in both sexes from April to July, including the difference
195 of feeding intensity between females and males in each month by Chi square fit
196 test.
Feeding Female Male Level of
activity (%) AF PAF AF PAF significance
10 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
Apr. 50.00 50.00 11.54 88.46 **
May. 45.00 55.00 15.00 85.00 *
Jun. 35.00 65.00 25.00 75.00 NS
Jul. 20.00 80.00 5.00 95.00 NS
197 *, P<0.05; **, P<0.01; NS, not significant
198
199 Molecular prey identification
200 All stomach contents were nearly digested into pulp and as a result, the prey
201 items were impossible to visually identify. A total of 59 stomach contents yielded
202 amplifiable DNA and 60 sequences were obtained, ranging from 500bp to 708 bp
203 (Table 4). All sequences were submitted to GenBank (Accession numbers
204 MK688462-MK688521). Six OTUs were established, with a maximum sequence
205 divergence of 0.1%. Of the 60 sequences, 59 clones showed similarities higher than
206 98% to reference sequences, allowing identification at species level. Only one clone
207 displayed 97.56% similarity to reference sequences, and it was assigned to the
208 Diopatra genus level (Fig 2). One stomach contained two kinds of prey, and the
209 remaining 58 samples contained only one kind. The duplicate parallel PCR products
210 from the contents of each stomach yielded the same sequences.
211
212 Table 4. Pray DNA detected in wild O. minor by cloning the COⅠ fragment gene,
213 including number per month and min-max length of sequences, GenBank
214 Accession numbers and Sequence ID of closest matches, percentages of similarity
11 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
215 obtained from BLAST and BOLD.
216
Ac. Min-max Similarity Pray Order Family Species Number or . length % DNA Seq. ID of Seq.
Acanthogobius JX679022, Perciformes Gobiidae 99.85,100 OUT1 500-700 flavimanus GU440207
Acanthogobius JQ738534, Perciformes Gobiidae 99.85,100 OUT2 640-697 hasta HQ536245
Tridentiger JX679061, MA07, Perciformes Gobiidae 100 684 bifasciatus HQ536534 MA11
Pholis FOJS101-1 Perciformes Pholidae 100 OUT3 556-694 crassispina 3
Charybdis HM237602, Decapoda Protunidae 100 OUT4 633-687 japonica HM237597
Alpheus Decapoda Alpheilae 100 HM180433 JU19 620 brevicristatus
Stomatopod Oratosquilla Squillidae 100 HM180739 OUT5 572-587 a oratoria
Eunicida Onuphidae Diopatra sp. 97.56 FJ428842 JU18 682
Hysterothylaciu Ascaridida Anisakidae 98.91 FJ907319 JU04 679 m aduncum
12 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
Octopodida HQ638215 Octopoda Octopus minor 99.66,100 OUT6 562 e NC015896
217
218 Fig 2. Maximum likelihood tree (ML tree) for all sequences obtained from
219 stomach contents and the closest matches of each sequence that were
220 downloaded from BOLD databases and GenBank, using MEGA 6 software
221 based on K2p model
222
223 Of the 60 sequences, 30 matched with fish (50.00%, by number), 13 with crustaceans
224 (21.66%), one with annelid (1.66%), one with nematode (1.66%) and 15 with its own
225 species (25.00%). A total of 8 species were identified as prey of O. minor (Table 2).
226 When considering the importance of prey to O. minor by number (%N), and ignoring
227 the nematode and its own species as prey, the families with the highest %N were
228 Gobiidae (61.36%), comprising Acanthogobius flavimanus (40.91%), Acanthogobius
229 hasta(15.91%) and Tridentiger bifasciatus (4.55%); Protunidae (18.18%), comprising
230 the one species Charybdis japonica; Pholidae (6.82%), comprising one species Pholis
231 crassispina; and Squillidae (9.09%), comprising one species, Oratosquilla oratoria.
232 Alpheilae and Onuphidae, accounted for same percentage, 2.27%, respectively
233 comprising one species each, Alpheus brevicristatus and Diopatra sp. The occurrence
234 frequency of Family of prey identified in the stomachs of O. minor caught in the
235 Moon Lake was showed in Fig 3.
236
13 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
237 Fig 3. Diet of O. minor in Swan lake using frequency of occurrence of prey items
238 for stomach content analysis by DNA Barcoding from April to July.
239
240 The number of prey species consumed by O. minor in a month was 2-5 from April to
241 July, and in June, the prey diversity was the most abundant, with 5 species. A.
242 flavimanus and C. japonica appeared with the highest occurring frequency during 3
243 months, followed by A. hasta. Cannibalism was found in 15 individuals, and the
244 cannibalism occurred 2, 3, 4 and 6 cases respectively from April to July, showing an
245 increasing trend. The cumulative prey curve approached the asymptote showing that
246 59 stomachs were adequate to describe the diet of this species (Fig 4).
247
248 Fig 4. The species-accumulation cures. The slopes curve rapidly approached
249 asymptotes, indicating that there were enough stomach content specimens
250 collected to detect most prey species.
251
252 Discussion
253 Based on our recent studies, the male octopus reaches maturity in April,
254 meanwhile, the females are at maturing stage when the female octopus begins egg
255 production actively, which mainly involves the synthesis of vitellus. At this stage the
256 female octopuses use energy directly from food, with no storage reserves being
257 transferred from organs to the gonads [34, 35]. Hence, it is unsurprising then that the
258 feeding intensity of the female was greater than that of the male. However, the
14 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
259 difference in feeding intensity between females and males was subsequently
260 indistinguishable as a considerable overall feeding intensity reduction occured in
261 both sexes caused by the listlessness of the male after mating and of the female after
262 spawn-movement [2, 36].
263 One advantage of molecular methods is that when morphological methods are
264 ineffective, sufficient DNA can be recovered by successful DNA amplifications [11,
265 12]. This can be achieved in poor quality samples as well because it requires only a
266 small amount of tissue for DNA extraction. Moreover, Meusnier et al. (2008) have
267 demonstrated that DNA barcoding could identify species with fragments as short as
268 100 bp with at least 90% efficiency [37]. The read lengths of the DNA fragments
269 obtained in this study was 500-700 bp, and the relative large read lengths strengthen
270 the accuracy of the identification. Only one clone was assigned to the genius level on
271 the basis of their supported topographical position on the bootstrap consensus tree
272 (Figure 2).
273 One intriguing discovery we found was that almost all of the stomach contents
274 (58 out of 59) contained only one kind of prey. In consideration of stochastic
275 sampling errors when subsampling stomach contents for DNA extraction and lack of
276 detection of some prey species arising from low-concentration DNA to PCRs [13],
277 duplicate subsamples with 20-40mg material each were used in this study. However,
278 the same sequences were obtained for duplicate subsamples showed that it was
279 unlikely that the lack of detection of diet composition was due to the procedural
280 errors. In addition, the universal COⅠ primers used in this study are able to amplify
15 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
281 COⅠ gene fragments from 11 invertebrate phyla and fish [30], as shown in previous
282 studies [38, 39]. A strong dietary preference probably played a decisive role in this
283 result. When presented with the same opportunity to different diets, O. minor makes
284 a strong preference, which has been observed in other similar octopods, such as O.
285 maorum [21], O. vulgaris [14] and Enteroctopus dofleini [18]. This situation has also
286 been observed in our previous studies, where when presented with the opportunity to
287 manila clams, the bivalve Potamocorbula laevis and the Asian shore crab, O. minor
288 frequently preyed on the Asian shore crab and rarely on the clams. Biandolino et al.
289 (2010) observed that only after octopus had voraciously fed on all available fish and
290 crabs, they then fed on mussels, as the last option [24]. Ambrose (1984) found that
291 the diet of O. bimaculatus in the field was influenced by food preferences, where
292 highly preferred prey (crabs, bivalves, sedentary grazers) were rare in the field,
293 however the most abundant species available, the snail Tegula eiseni Jordan, was the
294 least preferred [20]. The results showed that although there distributed a large of
295 mollusks, 15 molluscan species belonging to 14 families, in Moon Lake [40],
296 octopuses did not feed on mollusks. The diet of O. minor in the wild was also
297 influenced by food preferences. False negatives cannot be ruled out in this study.
298 PCR dropout is a common phenomenon in genotyping studies because PCRs are
299 known to amplify preferentially the DNA of a higher quality, and variable rates of
300 DNA degradation between different prey items may have biased PCR amplification
301 success [13]. The promise of high throughput methods in the future will only
302 improve this approach.
16 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
303 The prey consisted of 8 species in this study, a relatively narrow dietary for
304 octopus species. However, previous studies suggested that at least 12 prey species
305 were consumed by similar octopus though examination of their middens or via
306 morphological analyses of stomach contents [5, 6, 41]. From our analysis, the slopes
307 of the saturation curves rapidly approached asymptotes, which indicate that
308 sufficient stomach content specimens were collected to capture the major prey items
309 (Fig 4). In this study we were interested in the months when the females were
310 maturing and the dietary range wound likely be greater than reported here if the O.
311 minor stomach contents had been sporadically collected from Moon Lake over a
312 longer period [6]. Fish, including Gobiidae and Pholidae, accounted for the vast
313 majority proportion (68.19%) of the prey. The benthic fish, Gobiidae, in the lake
314 where the octopuses were collected, appear to be an especially important food
315 supply for females during sexual maturation and they provided a substantial and
316 favorite food supply for the octopods. When the preferred diet does meet the feeding
317 needs, O. minor would not feed the other candidate items [17, 18, 22, 24, 32], which
318 is a possible reason for the reduced number of types of consumed prey.
319 The genetic data indicate that cannibalism may be a significant feeding strategy
320 in long-armed octopus, with long-armed octopus DNA being detected in 15 of 59
321 stomachs. We do not believe that the host DNA contaminated the stomach contents
322 in these cases, as suckers and tops of arms were present in the stomachs and the host
323 was sound with no wounds. Cannibalism in octopus has been observed in artificial
324 culture [26, 42] and wild research [3, 43]. In this study, the number of cannibalism
17 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
325 cases increased from April to July.
326 Conclusions
327 The results confirm that although O. minor is a generalist and opportunistic
328 predator, it has strong dietary preferences in the wild. The discovery of the feeding
329 voracity of females during maturation provides a reference for aquaculture and
330 artificial breeding of this commercially important species, as well as accumulating
331 data for the rational resource exploitation of this kind animal. The DNA barcoding
332 method was shown to be successful in unraveling the feeding habits of wild O.
333 minor, which can simplify the field operation as stomach contents can be fixed by
334 the infiltration of ethanol in the field. This method shows higher taxonomic
335 resolution of the determination of prey items compared to traditional descriptions of
336 stomach contents. The results provide trophic relationship information for fishery
337 management, especially for the management of exclusive conservation reserves of O.
338 minor.
339 Acknowledgments
340 This study was supported by research grants from National Natural Science
341 Foundation of China (Nos. 31672257; 31172058) and Fundamental Research Funds
342 for the Central Universities(No.201822022); We are indebted to Mr. Cai B of
343 Mashan Group Co. Ltd (Shan-Dong Province, China) and Mr. Song MP for their
344 kind assistance in specimen collection.
345
346 Reference
18 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
347 1. Kim DH, An HC, Lee KH, Hwang JW. Optimal economic fishing efforts in Korean
348 common octopus Octopus minor trap fishery. Fisheries science, 2008; 74:
349 1215-1221. https://doi.org/10.1111/j.1444-2906.2008.01645.x
350 2. Bo QK, Zheng XD, Gao XL, Li Q. Multiple paternity in the common long‐armed
351 octopus Octopus minor (Sasaki, 1920) (Cephalopoda: Octopoda) as revealed by
352 microsatellite DNA analysis. Marine ecology. 2016; 37: 1073-1078.
353 https://doi.org/10.1111/maec.12364
354 3. Smale MJ, Buchan PR. Biology of Octopus vulgaris off the east coast of South
355 Africa. Marine Biology. 1981; 65: 1-12. https://doi.org/10.1007/BF00397061
356 4. Mather JA. Foraging, feeding and prey remains in middens of juvenile Octopus
357 vulgaris (Mollusca: Cephalopoda). Journal of Zoology. 1991; 224: 27-39.
358 https://doi.org/10.1111/j.1469-7998.1991.tb04786.x
359 5. Forsythe JW, Hanlon RT. Foraging and associated behavior by Octopus cyanea
360 Gray, 1849 on a coral atoll, French Polynesia. Journal of Experimental Marine
361 Biology and Ecology. 1997; 209: 15-31.
362 https://doi.org/10.1016/S0022-0981(96)00057-3
363 6. Ambrose RF, Nelson BV. Predation by Octopus vulgaris in the Mediterranean.
364 Marine Ecology. 1983; 4: 251-261.
365 https://doi.org/10.1111/j.1439-0485.1983.tb00299.x
366 7. Quetglas A, Alemany F, Carbonell A, Merella P, Sánchez P. Biology and fishery of
367 Octopus vulgaris Cuvier, 1797, caught by trawlers in Mallorca (Balearic Sea,
368 Western Mediterranean). Fisheries Research. 1998; 36: 237-249.
19 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
369 https://doi.org/10.1016/S0165-7836(98)00093-9
370 8. Grubert MA, Wadley VA, White RW. Diet and feeding strategy of Octopus
371 maorum in southeast Tasmania. Bulletin of Marine Science. 1999; 65: 441-451.
372 9. Riccialdelli L, Newsome SD, Fogel ML, Goodall R. NP. Isotopic assessment of
373 prey and habitat preferences of a cetacean community in the southwestern South
374 Atlantic Ocean. Marine Ecology Progress Series. 2010; 418: 235-248.
375 https://doi.org/10.3354/meps08826
376 10. Valdez-Moreno M, Quintal-Lizama C, Gómez-Lozano R, & del Carmen
377 García-Rivas M. Monitoring an alien invasion: DNA barcoding and the
378 identification of lionfish and their prey on coral reefs of the Mexican Caribbean.
379 PloS one. 2012: 7, e36636. https://doi.org/10.1371/journal.pone.0036636
380 11. Roura Á, González ÁF, Redd K, Guerra Á. Molecular prey identification in wild
381 Octopus vulgaris paralarvae. Marine Biology. 2012; 159: 1335-1345.
382 https://doi.org/10.1007/s00227-012-1914-9
383 12. Paquin MM, Buckley TW, Hibpshman RE, Canino MF. DNA-based identification
384 methods of prey fish from stomach contents of 12 species of eastern North
385 Pacific groundfish. Deep Sea Research Part I: Oceanographic Research Papers.
386 2014; 85: 110-117. https://doi.org/10.1016/j.dsr.2013.12.002
387 13. Braley M, Goldsworthy SD, Page B, Steer M, Austin JJ. Assessing morphological
388 and DNA-based diet analysis techniques in a generalist predator, the arrow squid
389 Nototodarus gouldi. Molecular Ecology Resources. 2010; 10: 466-474.
390 https://doi.org/10.1111/j.1755-0998.2009.02767.x
20 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
391 14. Smith CD. Diet of Octopus vulgaris in false bay, South Africa. Marine Biology,
392 2003; 143: 1127-1133. https://doi.org/10.1007/s00227-003-1144-2
393 15. Yamamoto T. On the ecology of Octopus variabilis typicus (Sasaki), with special
394 reference to its breading habits. The Malacological Society of Japan. 1942; 12:
395 9-20.
396 16. Andrews PLR, Tansey EM. The digestive tract of Octopus vulgaris: the anatomy,
397 physiology and pharmacology of the upper tract. Journal of the Marine
398 Biological Association of the United Kingdom, 1983; 63: 109-134.
399 https://doi.org/10.1017/S0025315400049845
400 17. Nixon M. Capture of prey, diet and feeding of Sepia officinalis and Octopus
401 vulgaris (Mollusca: Cephalopoda) from hatching to adult. Vie et Milieu. 1985.
402 18. Scheel D, Lauster A, Vincent T. Habitat ecology of Enteroctopus dofleini from
403 middens and live prey surveys in Prince William Sound, Alaska. Cephalopods
404 Present and Past: New insights and fresh perspectives: Springer; 2007. p. 434-58.
405 19. Cortez T, Castro BG, Guerra A. Feeding dynamics of Octopus mimus (Mollusca:
406 Cephalopoda) in northern Chile waters. Marine Biology. 1995; 123: 497-503.
407 https://doi.org/10.1007/BF00349228
408 20. Ambrose RF. Food preferences, prey availability, and the diet of Octopus
409 bimaculatus Verrill. Journal of Experimental Marine Biology and Ecology. 1984;
410 77: 29-44. https://doi.org/10.1016/0022-0981(84)90049-2
411 21. Ritchie L. Octopus Predation on Pot-caught Rock Lobster, Hokianga Area, NZ,
412 September-October 1970. Marine Department; 1972.
21 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
413 22. Castro BG, Guerra Á. The diet of Sepia officinalis (Linnaeus, 1758) and Sepia
414 elegans (D'Orbigny, 1835) (Cephalopoda, Sepioidea) from the Ria de Vigo(NW
415 Spain). Sci Mar. 1990; 54: 375-388.
416 23. Domingues P, Sykes A, Sommerfield A, Almansa E, Lorenzo A, Andrade JP.
417 Growth and survival of cuttlefish (Sepia officinalis) of different ages fed
418 crustaceans and fish. Effects of frozen and live prey. Aquaculture. 2004; 229:
419 239-254. https://doi.org/10.1016/S0044-8486(03)00351-X
420 24. Biandolino F, Portacci G, Prato E. Influence of natural diet on growth and
421 biochemical composition of Octopus vulgaris Cuvier, 1797. Aquaculture
422 international. 2010; 18: 1163-1175. https://doi.org/10.1007/s10499-010-9331-x
423 25. Dong ZZ. Fauna Sinica. Beijing: Science Press; 1988.
424 26. Bo QK, Zheng XD, Wang PL, Bi, KZ. Basic growth relations in experimental
425 rearing of newly hatchlings of Octopus minor (Sasaki, 1920). Oceanol ET
426 LimnolSinica. 2014; 45: 583-588. https://doi.org/10.11693/hyhz20140200056
427 27. Villanueva R, Nozais C, Boletzky SV. Swimming behaviour and food searching
428 in planktonic Octopus vulgaris Cuvier from hatching to settlement. Journal of
429 Experimental Marine Biology and Ecology. 1997; 208: 169-184.
430 https://doi.org/10.1016/S0022-0981(96)02670-6
431 28. Carreon‐Martinez L, Heath DD. Revolution in food web analysis and trophic
432 ecology: diet analysis by DNA and stable isotope analysis. Molecular Ecology.
433 2010; 19: 25-27. https://doi.org/10.1111/j.1365-294X.2009.04412.x
434 29. Pillay T. A critique of the methods of study of food of fishes. J zool Soc. 1952;
22 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
435 4:185-200.
436 30. Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for
437 amplification of mitochondrial cytochrome c oxidase subunit I from diverse
438 metazoan invertebrates. Molecular Marine Biology and Biotechnology. 1994; 3:
439 294-299.
440 31. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The
441 CLUSTAL_X windows interface: flexible strategies for multiple sequence
442 alignment aided by quality analysis tools. Nucleic acids research. 1997; 25:
443 4876-4882. https://doi.org/10.1093/nar/25.24.4876
444 32. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular
445 evolutionary genetics analysis version 6.0. Molecular biology and evolution.
446 2013; 30: 2725-2729. https://doi.org/10.1093/molbev/mst197
447 33. Colwell RK. Estimates: statistical estimation of species richness and shared
448 species from samples. Version 8.2. 2009. http://purl.oclc.org/ estimates.
449 34. Guerra A, Castro BG. Reproductive-somatic relationships in Loligo gahi
450 (Cephalopoda: Loliginidae) from the Falkland Islands. Antarctic Science. 1994;
451 6: 175-178. https://doi.org/10.1017/S0954102094000271
452 35. Rosa R, Costa PR, Nunes ML. Effect of sexual maturation on the tissue
453 biochemical composition of Octopus vulgaris and O. defilippi (Mollusca:
454 Cephalopoda). Marine Biology. 2004; 145: 563-574.
455 https://doi.org/10.1007/s00227-004-1340-8
456 36. Voight JR, Drazen JC. Hatchlings of the deep-sea octopus Graneledone
23 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
457 boreopacifica are the largest and most advanced known. Journal of Molluscan
458 Studies. 2004; 70: 400-402. https://doi.org/10.1093/mollus/70.4.400
459 37. Meusnier I, Singer GA, Landry JF, Hickey DA, Hebert PD, Hajibabaei M. A
460 universal DNA mini-barcode for biodiversity analysis. BMC genomics. 2008; 9:
461 214. https://doi.org/10.1186/1471-2164-9-214
462 38. Blankenship LE, Yayanos AA. Universal primers and PCR of gut contents to
463 study marine invertebrate diets. Molecular Ecology. 2005; 14: 891-899.
464 https://doi.org/10.1111/j.1365-294X.2005.02448.x
465 39. Méheust E, Alfonsi E, Le Ménec P, Hassani S, Jung JL. DNA barcoding for the
466 identification of soft remains of prey in the stomach contents of grey seals
467 (Halichoerus grypus) and harbour porpoises (Phocoena phocoena). Marine
468 Biology Research. 2015; 11: 385-395.
469 https://doi.org/10.1080/17451000.2014.943240
470 40. Liu JY, Li WT, Qin LZ, Zhang XM. Spatio-temporal variations in benthic and
471 macrobenthic molluscs in Swan Lake, Shandong, China. Marine Sciences. 2017;
472 41: 113. https://doi.org/10.11759/hykx20151022001
473 41. Rosas-Luis R, Jiménez Badillo MDL, Montoliu-Elena L, Morillo-Velarde PS. ().
474 Food and feeding habits of Octopus insularis in the Veracruz Reef System
475 National Park and confirmation of its presence in the southwest Gulf of Mexico.
476 Marine Ecology, 2019: e12535. https://doi.org/10.1111/maec.12535
477 42. Promboon P, Nabhitabhata J, Duengdee T. Life cycle of the marbled octopus,
478 Amphioctopus aegina (Gray) (Cephalopoda: Octopodidae) reared in the
24 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
479 laboratory. Scientia Marina. 2011; 75: 811-821.
480 https://doi.org/10.3989/scimar.2011.75n4811
481 43. Vincent TLS, Scheel D, Hough KR. Some Aspects of Diet and Foraging Behavior
482 of Octopus dofleini Wülker, 1910 in its Northernmost Range. Marine Ecology,
483 1998; 19: 13-29. https://doi.org/10.1111/j.1439-0485.1998.tb00450.x
25 bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. bioRxiv preprint doi: https://doi.org/10.1101/707349; this version posted July 18, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.