Holistic description of new deep sea megafauna (Cephalopoda: Cirrata) using a minimally invasive approach

Background In zoology, species descriptions conventionally rely on invasive morphological techniques, frequently leading to damage of the specimens and thus only a partial understanding of their structural complexity. More recently, non-destructive imaging techniques have successfully been used to describe smaller fauna, but this approach has so far not been applied to identify or describe larger animal species. Here, we present a combination of entirely non-invasive as well as minimally invasive methods that permit taxonomic descriptions of large zoological specimens in a more comprehensive manner. Results Using the single available representative of an allegedly novel species of deep-sea cephalopod (Mollusca: Cephalopoda), digital photography, standardized external measurements, high-field magnetic resonance imaging, micro-computed tomography, and DNA barcoding were combined to gather all morphological and molecular characters relevant for a full species description. The results show that this specimen belongs to the cirrate octopod (Octopoda: Cirrata) genus Grimpoteuthis Robson, 1932. Based on the number of suckers, position of web nodules, cirrus length, presence of a radula, and various shell characters, the specimen is designated as the holotype of a new species of dumbo octopus, G. imperator sp. nov. The digital nature of the acquired data permits a seamless online deposition of raw as well as derived morphological and molecular datasets in publicly accessible repositories. Conclusions Using high-resolution, non-invasive imaging systems intended for the analysis of larger biological objects, all external as well as internal morphological character states relevant for the identification of a new megafaunal species were obtained. Potentially harmful effects on this unique deep-sea cephalopod specimen were avoided by scanning the fixed animal without admixture of a contrast agent. Additional support for the taxonomic placement of the new dumbo octopus species was obtained through DNA barcoding, further underlining the importance of combining morphological and molecular datasets for a holistic description of zoological specimens.


Background
Invasive methods to obtain morphological character states still form the cornerstone of metazoan species descriptions, despite the increased use of molecular techniques [1,2]. For instance, describing a new cephalopod (Mollusca: Cephalopoda) species requires providing information on internal organs obtained through dissection [3]. Inevitably, this approach involves damage to or even partial destruction of the specimen at hand and therefore may preclude analysis of singular, endangered, rare, or otherwise valuable organisms. Furthermore, invasive techniques invariably alter the structural integrity of zoological specimens and thus do not permit studying organ systems in their natural context, ultimately resulting in a significantly restricted representation of the complexity of an organism. In contrast, digital three-dimensional (3D) imaging techniques such as computed tomography (CT) or magnetic resonance imaging (MRI) permit minimally or even entirely non-invasive analyses of whole biological specimens or parts thereof [4][5][6][7][8]. Although the application of such methods-an approach more recently labelled morphomics [9]-has so far been successfully demonstrated for the identification or description of smaller fauna [10][11][12], larger metazoan species have, for various technical or logistical reasons, thus far not been the target of this type of analysis.
Amongst cephalopods, the finned or dumbo octopuses (Octopoda: Cirrata) constitute relatively rare organisms, despite forming a significant part of the megafauna in deep-sea habitats of the World Ocean down to at least 7000 m depth [13]. With several cirrate species known only from a single collected specimen [14], a minimally invasive approach to describing new or redescribing existing taxa would constitute a significant improvement over conventional methods [15]. Based on the application of two complementary 3D imaging techniques supplemented by digital photography, standardized measurements of external features and DNA barcoding, we here present a novel workflow that results in the first minimally invasive, holistic description of a new metazoan species pertaining to a taxon composed entirely of megafaunal representatives.

Specimen collection, fixation, and preservation
The single cirrate specimen was collected in the North Pacific Ocean during scientific cruise SO-249 BERING (R/V SONNE) using a chain bag dredge [16]. The animal arrived dead on deck and was immediately transferred to a bucket filled with cold seawater. Several small tissue samples were collected from a single damaged arm using a scalpel and scissors. These tissue samples were then placed inside small plastic vials each filled with one of the following fixatives: 2 ml of 100% ethanol, 2 ml of 4% paraformaldehyde solution, 2 ml of 1:1 acetone/methanol solution, and 5 ml of RNAlater (Merck KGaA, Darmstadt, Germany). Following photography and measurements, the entire organism was placed inside a large plastic drum filled with 5 l of a 10% formalin solution prepared using distilled water supplemented with pH 7.2 microscopy buffer tablets (Merck KGaA). After several months, the specimen was transferred into 70% ethanol using a graded ethanol series. The specimen is deposited at the Museum für Naturkunde in Berlin, Germany (Zoologisches Museum Berlin Molluskensammlung, ZMB MOLL 240160).

Video, photography, and standardized measurements
A short movie of the adult specimen arriving on deck (Additional file 1) was recorded using a Galaxy S7 smartphone with an integrated digital camera (Samsung Corp., Seoul, South Korea). The movie was edited using the Windows 10 Video-Editor software (Microsoft Corp., Redmond, CA, USA) and saved as an MPEG-4 file. In addition, the specimen was photographed on-board ship using a DSC-HX 400 digital camera with integrated flash (Sony Corp., Tokyo, Japan). Before fixation, external morphological characters were measured to the nearest millimetre using a ruler and calliper following established guidelines [3]. To document gross morphological changes potentially caused by fixation and preservation [14], the specimen was photographed prior to MRI after several months in the 10% formalin solution as well as prior to micro-computed tomography (µCT) after several months in the 70% ethanol solution.

Magnetic resonance imaging
For MRI, the formalin-fixed specimen was placed inside a cylindrical plastic container filled with the original 10% formalin solution. Scanning was performed using a 7-T high-field Magnetom clinical MR (magnetic resonance) system with a 600-mm magnet bore and equipped with an SC72 shielded gradient set with a maximum strength of 70 mT/m (Siemens, Berlin, Germany). The container was placed horizontally inside a 32-channel human head coil (Siemens). For imaging, a 3D magnetizationprepared spiral acquisition gradient echo (MP-SAGE) MR sequence with 3000 ms repetition time, 3.4 ms echo time, 7°flip angle, 37 frames averaged, 150 × 111 × 150 mm field of view, 536 × 396 × 536 px matrix size, 280 μm isotropic voxel resolution, and about 16 h 59 min acquisition time was used. The selected field of view comprised a region of interest extending from the posterior mantle edge to about the middle of the arms and thus covered all the internal structures relevant for a cirrate species description [3]. The acquired data were reconstructed using the software syngo MR B17 (Siemens). For further image processing, the original 16-bit NIfTI file was transformed into an 8-bit TIFF image stack, zero-filled to 140-μm isotropic voxel resolution, and finally rotated to a dorsoventral orientation using the software Fiji 1.52v [17].

Micro-computed tomography
For μCT, the ethanol-preserved specimen was placed inside a cylindrical plastic container filled with the original 70% ethanol solution. Scanning was performed using a Phoenix v|tome|x s 180/240 CT system (GE Sensing & Inspection Technologies, Wunstorf, Germany) equipped with a 180-kV X-ray source and a detector measuring 2024 × 2024 px. Scan parameters were as follows: 100 kV source voltage, 100 μA source current, no filter, 500 ms exposure, 7 averages, 0.24°step size with 1500 projections over 360°, 13.2 × 14.4 × 13.2 mm field of view, 1210 × 1314 × 1204 px matrix size, 10.93 μm isotropic voxel resolution, and about 1 h 28 min acquisition time. The selected field of view focused on the buccal mass area. The resulting 16-bit projection images were reconstructed using the software Phoenix datos|x 2.7 (GE Sensing & Inspection Technologies). For further image processing, the original 16-bit RAW image volume was transformed into an 8-bit TIFF image stack and rotated to a sagittal orientation using the software Fiji 1.52v.

Three-dimensional reconstruction and visualization
Based on the MRI 8-bit TIFF image stack, manual segmentation of selected internal organ systems was carried out using the Segmentation Editor in the software Amira 6.1 (Thermo Fisher Scientific Inc., Waltham, MA, USA). Surface rendering of the labelled organs was performed using the SurfaceGen module with Constrained smoothing activated. The resulting mesh was reduced from about 2,500,000 to 150,000 faces using the Simplifier function. Finally, the 3D model was saved as a WRL file using Amira's VRML-Export module with Render smooth and Render specular activated. A similar approach was chosen to reconstruct the upper and lower beak based on the μCT 8-bit TIFF image stack. In addition, selected elements of the internal anatomy of the specimen were visualized in 3D using the Volren module in Amira. Virtual two-dimensional (2D) sections of the 8-bit MRI and μCT datasets were created using the Volume Viewer 2.0 plugin in the software Fiji 1.52v.

Interactive 3D model creation
Two interactive 3D PDF files (Additional files 2 and 3) were created using the software Adobe 3D Reviewer 9.5.5 and Adobe Acrobat Pro Extended 9.5.5 (Adobe Systems Inc., San José, CA, USA). To this end, the WRL files exported from Amira were loaded into Adobe 3D Reviewer, where lighting, background, orientation, and labelling settings were adjusted. The models were each exported as PDF files, which were then loaded into Adobe Acrobat, where the cover image and pre-saved views were generated in order to finalize the interactive 3D PDF files. Please refer to previously published articles on how to use [18] or create interactive 3D PDF files [19][20][21].

Phylogenetic analysis
General time reversible (GTR) with gamma distribution was selected as the model for phylogenetic reconstruction using MrModeltest2 2.4 [28]. Phylogenetic trees were reconstructed in MEGA 6.06 using maximum likelihood (ML) as optimality criterion [29]. Branch support was estimated using 500 bootstrap replicates. Following an initial comparison, the 16S rRNA gene sequence was favoured over the COI gene sequence due to significantly broader taxon sampling in GenBank for the former and thus a better resolution of the resulting phylogeny. The 16S gene sequence of the new species (MW575539) was compared with one Grimpoteuthis sp. sequence (AF110100) from the North Pacific (note that this specimen is very likely an Opisthoteuthis), other Grimpoteuthis sp. sequences (AF487305-AF487312) from the North Atlantic (in fact all G. discoveryi -M. Collins, personal communication), and all further cirrate (Octopodiformes: Octopoda: Cirrata) 16S sequence data available from GenBank. In addition, further octopodiform specimens were incorporated as putative outgroup taxa, including several incirrate (Octopodiformes: Octopoda: Incirrata) and vampire squid (Octopodiformes: Vampyromorpha) specimens (Table 1).

Results
The combined use of the non-invasive techniques digital photography, external measurements, MRI, and μCT complemented with DNA barcoding based on minimally invasive tissue sampling reveals that the single, wellpreserved cirrate collected in the North Pacific Ocean is the first representative of a previously undescribed species of dumbo octopus. The full species description given below follows previously published guidelines [3] as well as two of the most recent descriptions of other new Grimpoteuthis species [30,31]

Diagnosis
Medium-sized species with moderately long, lateral fins. Cirri short and suckers moderate. Gills compact with eight broad lamellae. Radula present, teeth homodont. Paired anterior and unpaired posterior salivary glands present. Shell U-shaped, smooth with lateral wings parallel, broadly tapering towards distal ends.

Description
Medium-sized species, body semi-gelatinous, bell-shaped ( Fig. 1a-e). Mantle about one third of total length (MLI 32.8), saccular, gelatinous, nearly twice as long as wide, broadly rounded posteriorly. Head wider than mantle (HWI 23.8), no neck region visible (Fig. 1f). Pallial aperture narrow, enveloping base of funnel (Fig. 1e). Funnel long (FuLI 30.4), protruding far beyond mantle margin, distal third of funnel free (Fig. 1b). Funnel broader at the base than at the opening (Fig. 1g). Funnel organ of inverted V shape. Olfactory organs rounded, situated within the pallial aperture on either side of the funnel base (Fig. 1g). Fins moderately long (FSI 68.6, FLI 100.0), about half as wide as long (Fig. 1b). Fins situated laterally, between the eyes and apex of the mantle, but closer to the eyes. Posterior fin margin straight, less gelatinous than the rest of the mantle. Tips of fins broadly rounded. Anterior fin margin slightly convex, tapering posteriorly towards the base. Fin cartilage and fin muscles inserting in the posterior portion of the fins (Fig.  1d). Broad muscular base, attached to the shell. Towards the tips of fins, the fin cartilage becomes narrower, ending in a small tip about ¾ of fin length. Arms semi-gelatinous, subequal, moderately long (ALI 63.4). The fourth arm on ventral side shortest (ca. 54% of total length), first arm on dorsal side longest (ca. 63% of total length), arm formula varying between left and right sides. The arms deeply set into the soft and fleshy primary web, no secondary web present (Fig. 1c). Web formula differing between right and left sides with sector A or B deepest and sector D or E shallowest, respectively. Web nodules large and rounded, located on the ventral side of the arms between suckers 33 and 34, uniting web to the ventral side of the arms (Fig. 1c). Beyond the nodule, the web extends only a few suckers further, becoming very narrow. Nodules absent on the dorsal surface, web extending nearly to the tip of each arm. No nodules present in sector A, so the web extends nearly to the tip of both adjacent arms. In sector E, two nodules present, so the web extends only a short distance beyond the nodules, thus rendering this sector one of the shallowest sectors and sector A one of the deepest. Suckers in a single row extending from the mouth opening to the tips of the arms (Fig. 1c). Arms with 68-74 suckers set into the oral surface, their apertures projecting freely. First five suckers very small, followed by suckers gradually increasing in size (MSDI 4.3). No enlarged suckers around the web nodule, no sign of hectocotylization (Fig. 1c). Suckers cylindrical, with short, narrow apertural rims (Fig. 1h). Simple sucker aperture, the base of suckers embedded into primary web. Cirri short, 1.0× maximum sucker diameter (MCLI 4.3), located on the oral surface of the arms between base of suckers (Fig. 1h). First cirrus between suckers 4 and 5. Cirri extending to the tips of the arms with no change in length along the arms.
Gills compact (GDI 13.1), almost spherical (Fig. 1i), of 'half-orange' type, with eight broad lamellae each (Fig. 1j). Gills partly surrounding ovoid-shaped branchial glands. Teardrop-shaped branchial hearts located dorsoposterior to the gills, in close proximity to the renal appendages. Systemic heart with single large curved ventricle, indistinguishable auricles (Fig. 1k). Small swelling on the dorsal aorta directly dorsal to the ventricle. Dorsal aorta forming single large curve across the dorsal side of the oesophagus. Two large pallial arteries leaving the dorsal aorta as a single vessel each near the posterior end of the oesophagus. Eyes large (EDI 24.6), not protruding, with prominent lens (Fig. 2a, b). Optic lobes almost spherical in shape, with a single bundle of optic nerves passing through the white body (Fig. 2c) before tapering into several smaller nerves closer to the eye. White bodies more than twice the size of the optic lobes, somewhat spherical in shape with two processes directed towards the brain. Statocysts almost as large as the eyes (Fig. 2d). Stellate ganglion ovoid in shape, epistellar body morphologically not distinguishable (Fig. 2e).
The digestive tract in lateral view C-shaped with relatively large buccal mass (Fig. 2f-h). The upper beak without distinct folds, ridges, or thickenings (Fig. 2i, j). Lower beak with rounded hood and broad wings, lateral walls without any folds or ridges (Fig. 2k, l). Paired anterior salivary glands centrally located within the buccal mass, to the left and right side above the radular sac (Fig.  2f). Large median submandibular gland below radular sac (Fig. 2f, h). Radula poorly developed, with homodont dentition arranged in five rows (Fig. 2m, n). Rhachidian tooth slender, smaller than laterals. First lateral tooth broad triangular, second lateral tooth slender triangular, no trace of marginal teeth or plates. Oesophagus long, widening in  (Fig. 2f). Unpaired posterior salivary gland small, located on the ventral side of the oesophagus. Stomach nearly rectangular, tapering towards the ventral side. The caecum slightly smaller than the stomach with two long and slender hepatic ducts leading to the large, unilobular, almost spherical digestive gland (Fig. 2f). Small pancreas present on the posterior surface of the digestive gland. Intestine shorter than the oesophagus, no enlargement. Rectum located close to the funnel. Anal flaps and ink sac absent.
Shell located in the dorsal part of the mantle cavity, robust, U-shaped (Fig. 2a, b). Shell wings not expanding, parallel to each other, broadly tapering towards the distal ends (Fig. 2o). Shoulder blades present, but weak. Outer surface of the shell saddle convex without median ridge or transverse groove. Fin cartilage insertions long, encompassing about ⅔ of the entire length of the shell wing (Fig. 2p).
Male reproductive system with large, almost spherical testis located centrally in the posterior part of the mantle, posterior to the stomach and caecum, ventral to the shell (Fig. 2a). Vas deferens elongated. Spermatophoric glands I-III moderately developed, convoluted (Fig. 2q). Accessory spermatophoric gland complex much larger than the spermatophoric gland complex. Accessory spermatophoric glands 1 and 3 nearly the same size, much larger than accessory spermatophoric gland 2 (Fig. 2r). Several spermatophores located inside the central duct within the accessory spermatophoric gland complex, each measuring ca. 2 × 1 mm. Terminal organ (penis) long and slender, directed towards the rectum (Fig. 2r).
Skin surface smooth. Dorsal surface of the head and mantle white, slightly reddish towards the posterior apex (Fig. 1a). Ventral surface of the head and mantle white with a red pigmentation, funnel darker red in colour, in particular at the distal end (Fig. 1b). Posterior margin of the fins deep red, becoming lighter and almost white towards the anterior margin. Dorsal arms on the surface white, dorsal surface of the web reddish (Fig. 1a). Ventral arms and web dark red, suckers lighter red (Fig. 1c).

Measurements and counts
Most morphometric and meristic characters were obtained from the holotype directly following capture ( Table 2).

Distribution
So far known only from the type locality in the northern part of the Emperor Seamounts, an undersea mountain chain in the northwestern part of the North Pacific [34].

Remarks
Based on the shell form, fin position, optic lobe shape, arm length, web form, and optic nerve arrangement, this specimen is readily identifiable as a cirrate of the genus Grimpoteuthis [14]. Differentiation from previously described species is based on the number of suckers (68-74), absence of enlarged suckers near the web nodules, absence of a thin web between the suckers, position of the web nodules (near suckers 33-34), cirrus length (1.0× maximum sucker diameter), location of proximal cirri (between suckers 4 and 5), presence of a radula, gill shape ('half-orange' type), gill lamellae count (8/8), and various shell characters (Table 3). Apart from the new species described here, seven other Grimpoteuthis species as well as several unidentified specimens ascribed to this genus have been recorded from the Pacific Ocean (Fig. 3). The type localities of the two geographically closest species, i.e. G. bathynectes and G. tuftsi, are also found in the North Pacific, but in areas more than 4000 km east of the type locality of G. imperator sp. nov. [30]. Further type localities of Pacific species pertaining to that genus are found in the tropical and southern parts of the ocean. Whilst only a single individual was analysed here, this is not uncommon in cirrate taxonomy due to the scarcity of suitable material [14] and under certain circumstances (e.g. 'unquestioned distinctiveness of the species' and 'the high probability that no additional material will be forthcoming soon') does not impede designation of a new species ([3]: 49).

DNA barcoding and phylogenetic analysis
In addition to the non-invasively acquired morphological character states, minimally invasive DNA barcoding was used to obtain the first molecular sequence data for an unambiguously identified member of the genus Grimpoteuthis from the Pacific Ocean. Inferences based on this as well as 55 previously deposited 16S rRNA gene sequences result in the first phylogeny with full coverage of all eight extant cirrate genera (Fig. 4).
The phylogenetic analysis provided very good resolution within and between cirrate genera as well as all outgroup taxa. The North Atlantic Grimpoteuthis

Discussion
In recent years, several non-invasive imaging techniques have become available [4][5][6][7][8] that permit acquisition of high quality structural data in 3D from various metazoan taxa, including terrestrial and aquatic megafauna [33,[35][36][37]. In the present study, MRI was performed overnight using a high-field clinical imaging system [33]. Potentially harmful long-term effects on the singular specimen were avoided by scanning the organism in its original fixative without any contrast agent added. At 140 μm isotropic voxel resolution, the signal to noise ratio of the resulting 3D dataset was excellent (Fig. 1f), allowing extraction of all taxonomically relevant internal morphological characters such as shell  Location of proximal cirri between suckers 4-5 Gill diameter (GD) 8.8 Gill lamellae count R/L (GiLC) 8/8 and gill shape or digestive tract morphology, but also of more minute structures such as nervous system and sensory organ morphology (Fig. 2c). In addition, morphological characters so far not used in cirrate species descriptions, but potentially of taxonomic value such as the shape of the systemic heart or the arrangement of the dorsal aorta [38], were here integrated into the description of a new cirrate species for the first time as well (Fig. 1k).
The capacity of MRI to visualize water-rich tissues has in the past been applied to a broad spectrum of zoological specimens using different MR modalities such as post mortem, in vivo, or diffusion tensor imaging [4-8, 15, 20, 33, 35-37, 39-41]. However, due to the tissue properties of cephalopod beaks and radula (i.e. water-poor chitin) as well as a strong susceptibility artefact in the buccal mass area (Fig. 1f) caused by ingested sediment [36,41], a complementary region of interest μCT scan was conducted at 10.93 μm isotropic voxel resolution. The improved hard part contrast of this X-ray-based imaging technique proved more suitable for the visualization of chitinous tissues (Fig. 2) and allowed building the first interactive 3D model of a cephalopod beak (Additional file 3). Although analysis of a new megafaunal species using two complementary imaging modalities with different fields of view was successfully performed here, a single scan of the entire specimen with combined soft and hard tissue contrast would obviously have been more desirable. Such studies are bound to be possible in the near future, for example, on the upcoming BM18 beamline of the European Synchrotron Radiation Facility that will permit phase-contrast tomography at nanometre spatial resolution on specimens of up to 2.5 m length [42]. In addition, recent developments based on the implementation of deep neural networks are poised to significantly increase the efficiency of semi-automated 3D dataset segmentation [43]. Apart from allowing 2D virtual sectioning at arbitrary angles (Figs. 1f, g; 2c, h, m), volume rendering of selected areas (Figs. 1i-k; 2e, n), or the creation of interactive 3D models (Additional files 2, 3), the MRI and μCT datasets underlying the present species description can be employed for further analyses that would not have been possible to this extent using conventional, invasive techniques. For example, volumetric data rapidly obtainable from the organ systems reconstructed here could be compared with similar data previously gathered non-invasively from a Grimpoteuthis hatchling using high-field preclinical MRI [44]. This type of analysis would permit drawing conclusions on the lifestyle and behaviour of these hard-to-observe deepsea organisms [45].
In general, such inferences are undoubtedly facilitated by the deposition of the respective digital raw datasets in online repositories [5,33,46,47], an approach of particular importance for zoological taxonomy [10][11][12]. Apart from allowing improved or ideally full data transparency as well as data mining and modelling, dataset deposition in online repositories could in the future permit designation of so-called cybertypes [48]. Also termed virtual or e-type [49], a cybertype constitutes a digital, graphical representation of the physical specimen on which a description is based [10,11]. Whilst the International Commission on Zoological Nomenclature presently does not consider digital copies of physical specimens as sufficiently adequate type material [50], the provision of digital cybertype infrastructure has nonetheless been identified as an important step for overcoming the so-called taxonomic impediment [51]. The two publicly accessible digital repositories for molecular (NCBI Fig. 3 Geographic distribution of identified and unidentified specimens of Grimpoteuthis in the Pacific Ocean. See Table 3 for a list of the respective type localities GenBank) and protein sequence data (RCSB Protein Data Bank) provide striking evidence for the benefit of a centralized and professionally curated online database [5,48]. In the meantime, the data gathered in the course of this study were deposited on MorphoBank [52].

Conclusions
By extending the morphomics concept to the description of a new species of megafauna, we here show that a minimally invasive approach based on the application of complementary non-invasive 3D imaging techniques supplemented with molecular sequence data can help to advance metazoan taxonomy, in particular, in cases where valuable, larger zoological specimens require a more detailed, holistic analysis.

Supplementary Information
The online version contains supplementary material available at https://doi. org/10.1186/s12915-021-01000-9. interactive 3D PDF model (requires Adobe Reader 9 or higher on all operating systems). Use mouse wheel to zoom in or out. A right-click on the activated model provides access to further options such as a set of pre-saved views, a full-screen mode, the model tree icon, or deactivation of the 3D model. Branchial glands = dark green, branchial hearts = light brown, branchial heart appendages = maroon, central nervous system = medium yellow, circulatory system = dark brown, digestive gland = blue, digestive tract = dark blue, eyeballs = orange, fin cartilages = pink, gills = light green, hepatic ducts = turquoise, lenses = dark grey, optic lobes = dark yellow, pancreas = light blue, Fig. 4 Phylogenetic analysis of the Cirrata and other selected Octopodiformes based on the 16S rRNA gene sequence. Scale bar refers to a phylogenetic distance of 0.02 nucleotide substitutions; new species marked in bold font; numbers on the branches denote bootstrap values after 500 replicates peripheral nervous system = light yellow, posterior salivary gland = bluegreen, renal appendages = magenta, reproductive system = light grey, shell = red, statocysts = purple, ureters = violet, white bodies = white.
Additional file 3. Interactive 3D model of upper and lower beak of Grimpoteuthis imperator sp. nov. (ZMB MOLL 240160). Left-click anywhere on this figure to open the embedded, interactive 3D PDF model (requires Adobe Reader 9 or higher on all operating systems). Use mouse wheel to zoom in or out. A right-click on the activated model provides access to further options such as a set of pre-saved views, a full-screen mode, the model tree icon, or deactivation of the 3D model.