The common octopus (Octopus vulgaris) is an attractive species for aquaculture, however, several challenges inhibit sustainable commercial production. Little is known about the early paralarval stages in the wild, including diet and intestinal microbiota, which likely play a significant role in development and vitality of this important life stage. High throughput sequencing was used to characterize the gastrointestinal microbiome of wild O. vulgaris paralarvae collected from two different upwelling regions off the coast of North West Spain (n = 41) and Morocco (n = 35). These were compared to that of paralarvae reared with Artemia for up to 25 days in captivity (n = 29). In addition, the gastrointestinal microbiome of zooplankton prey (crabs, copepod and krill) was also analyzed to determine if the microbial communities present in wild paralarvae are derived from their diet. Paralarvae reared in captivity with Artemia showed a depletion of bacterial diversity, particularly after day 5, when almost half the bacterial species present on day 0 were lost and two bacterial families (Mycoplasmataceae and Vibrionaceae) dominated the microbial community. In contrast, bacterial diversity increased in wild paralarvae as they developed in the oceanic realm of both upwelling systems, likely due to the exposure of new bacterial communities via ingestion of a wide diversity of prey. Remarkably, the bacterial diversity of recently hatched paralarvae in captivity was similar to that of wild paralarvae and zooplankton, thus suggesting a marked effect of the diet in both the microbial community species diversity and evenness. This study provides a comprehensive overview of the bacterial communities inhabiting the gastrointestinal tract of O. vulgaris paralarvae, and reveals new research lines to challenge the current bottlenecks preventing sustainable octopus aquaculture.

Common octopus, Octopus vulgaris, is an economically important cephalopod species. However, its rearing under captivity is currently challenged by massive mortalities previous to their juvenile stage due to nutritional and environmental factors. Dissecting the genetic basis and regulatory mechanism behind this mortality requires genomic background knowledge. A transcriptomic sequencing of 10 dph octopus paralarvae from different experimental conditions was constructed via RNA-seq. A total of 613,767,530 raw reads were filtered and de novo assembled into 363,527 contigs of which 82,513 were annotated in UniProt carrying also their GO and KEGG information. Differential gene expression analysis was carried out on paralarvae reared under different diet regimes and temperatures, also including wild paralarvae. Genes related to lipid metabolism exhibited higher transcriptional levels in individuals whose diet includes crustacean zoeas, which had an impact over their development and immune response capability. High temperature induces acclimation processes at the time that increase metabolic demands and oxidative stress. Wild individuals show an expression profile unexpectedly similar to Artemia fed individuals. Proteomic results support the hypothesis revealed by transcriptional analysis. The comparative study of the O. vulgaris transcriptomic profiles allowed the identification of genes that deserve to be further studied as candidates for biomarkers of development and health. The results obtained here on the transcriptional variations of genes caused by diet and temperature will provide new perspectives in understanding the molecular mechanisms behind nutritional and temperature requirements of common octopus that will open new opportunities to deepen in paralarvae rearing requirements.

The fatty acid (FA) profile of wild Octopus vulgaris paralarvae of estimated age was individually analyzed for the very first time in order to establish a reference for comparison in rearing and nutritional studies. Age of each paralarvae was estimated by analysing daily increments on lateral hood surface of beaks. Wild paralarvae age ranged between 6 and 8 days and their FA composition resembled that from hatchlings produced under culture conditions. However, when compared with the FA composition of up to 20 days old cultured paralarvae described in the bibliography, some striking differences were found. Results showed higher levels of docosahexaenoic acid (22:6n − 3, DHA), lower contents of 18:1n − 9, 18:1n − 7 and 18:2n − 6 and negligible levels of 18:3n − 3 in wild paralarvae, when collated to reared one. These results seem to indicate that preys/diets supplied to cultured paralarvae fail to resemble paralarval natural composition and as a result do not fulfil their FA requirement. The individual applied technique developed in this study will allow to refine the study of wild paralarvae along its development, as well as to compare wild and cultured paralarvae of similar age. Statement of relevance Artemia does not fulfil paralarval fatty acid requirements.This study was funded by Spanish Government under Projects OCTOPHYS (Ref. AGL2010-22120-C03), OCTOWELF (Ref. AGL2013-49101-C2-1-R) and LARECO (CTM2011-25929). Most of the authors of this paper participate to the COST network, Action FA1301 (CephsInAction) focussing on Cephalopod welfare. D. Garrido was financed by a Ph.D. grant by Spanish Institute of Oceanography (BOE 3rd November 2011).Peer Reviewe

Octopus vulgaris is a species of great interest in research areas such as neurobiology, ethology and ecology but is also a candidate species for aquaculture as food resource as well as source to alleviating the fishing pressure on its wild populations. This study aimed to characterize the predatory behaviour of O. vulgaris paralarvae and to quantify their digestive activity. Those processes were affordable using the video-recording analysis of 3 dph, mantle-transparent paralarvae feeding on 18 types of live zooplanktonic prey. We show that octopus paralavae attack, immobilize, drill and ingest live cladocerans and copepods with 100% efficiency, which decreases dramatically to 60% on decapod prey (Pisidia longicornis). The majority (85%) of successful attacks targeted the prey cephalothorax while unsuccessful attacks either targeted the dorsal cephalothorax or involved prey defensive strategies (e.g. juvenile crab megalopa) or prey protected by thick carapaces (e.g. Gammaridae). After immobilization, the beak, the buccal mass and the radula were involved in exoskeleton penetration and content ingestion. Ingestion time of prey content was rapid for copepods and cladocerans (73.13 ± 23.34 s) but much slower for decapod zoeae and euphasiids (152.49 ± 29.40 s). Total contact time with prey was always <5min. Contrary to the conventional view of crop filling dynamics observed in adult Octopus vulgaris, food accumulated first in the stomach of paralarvae and the crop filled after the stomach volume plateaued. Peristaltic crop contractions (~18/min) moved food into the stomach (contractions ~30/min) from where it passed to the caecum. Pigmented food tracers were seen to enter the digestive gland, 312 ± 32 s after the crop reached its maximum volume. Digestive tract contents passed into the terminal intestine by peristalsis (contraction frequency ~50/min) and defaecation was accompanied by an increased frequency of mantle contractions. Current results provide novel insights into both, O. vulgaris paralarvae – live prey capture strategies and the physiological mechanisms following ingestion, a key information to develop an effective rearing protocol for O. vulgaris.

The European sardine (Sardina pilchardus, Walbaum 1792) is indisputably a commercially important species. Previous studies using uneven sampling or a limited number of makers have presented sometimes conflicting evidence of the genetic structure of S. pilchardus populations. Here, we show that whole genome data from 108 individuals from 16 sampling areas across 5000 km of the species’ distribution range (from the Eastern Mediterranean to the archipelago of Azores) support at least three genetic clusters. One includes individuals from Azores and Madeira, with evidence of substructure separating these two archipelagos in the Atlantic. Another cluster broadly corresponds to the center of the distribution, including the sampling sites around Iberia, separated by the Almeria–Oran front from the third cluster that includes all of the Mediterranean samples, except those from the Alboran Sea. Individuals from the Canary Islands appear to belong to the Mediterranean cluster. This suggests at least two important geographical barriers to gene flow, even though these do not seem complete, with many individuals from around Iberia and the Mediterranean showing some patterns compatible with admixture with other genetic clusters. Genomic regions corresponding to the top outliers of genetic differentiation are located in areas of low recombination indicative that genetic architecture also has a role in shaping population structure. These regions include genes related to otolith formation, a calcium carbonate structure in the inner ear previously used to distinguish S. pilchardus populations. Our results provide a baseline for further characterization of physical and genetic barriers that divide European sardine populations, and information for transnational stock management of this highly exploited species towards sustainable fisheries.