Bioluminescence of the Largest Luminous Vertebrate, the Kitefin Shark, Dalatias licha: First Insights and Comparative Aspects

Bioluminescence has often been seen as a spectacular yet uncommon event at sea but considering the vastness of the deep sea and the occurrence of luminous organisms in this zone, it is now more and more obvious that producing light at depth must play an important role structuring the biggest ecosystem on our planet. Three species of deepwater sharks (Dalatias licha, Etmopterus lucifer, and Etmopterus granulosus) were collected from the Chatham Rise, off New Zealand, and for the first time, we documented their luminescence. Comparison of glowing shark pictures, combined with histological description of light organs and hormonal control analysis, highlight the evolutive conservation of the bioluminescence process within Dalatiidae and Etmopteridae. A special emphasis is placed on the luminescence of D. licha, the largest known luminous vertebrate. This first experimental study of three luminous shark species from New Zealand provides an insight into the diversity of shark bioluminescence and highlights the need for more research to help understand these unusual deep-sea inhabitants: the glowing sharks.


INTRODUCTION
Bioluminescence, defined as the production of visible light by living organisms, is a widespread phenomenon mainly encountered among various marine taxa (Widder, 1999;Haddock et al., 2010). This living light, also called cold light, occurs through a biochemical reaction; the oxidation of a substrate, a luciferin, by an enzyme, the luciferase, or through a stabilized complex called photoprotein (Shimomura, 2006). Among Squaliformes, bioluminescence is documented for two deep-sea families: Dalatiidae and Etmopteridae (Claes and Mallefet, 2009b;Straube et al., 2015). A third family, Somniosidae was recently suggested to also contain a luminous species, Zameus squamulosus (Günther, 1877), based on density and upper view of putative light organs (i.e., photophores) , new results brought clear evidence Z. squamulosus being a luminous species . The first mentions of shark light emission date back to the nineteenth century (Bennett, 1840;Johann, 1899), but it is only recently that bioluminescence studies, focusing on physiological control, and photophore morphology and function, have been developed. These studies investigated bioluminescence in three etmopterids, Etmopterus spinax (Linnaeus, 1758), Etmopterus molleri (Whitley, 1939), Etmopterus splendidus (Yano, 1988), and one dalatiid, Squaliolus aliae (Teng, 1959) (e.g., Mallefet, 2009b,c, 2015;Claes et al., , 2011bClaes et al., , 2012Renwart et al., 2014Renwart et al., , 2015Duchatelet et al., 2019bDuchatelet et al., , 2020b. Luminous sharks appear to produce blue-green light (between 455 and 486 nm; Claes et al., 2014a) for multiple purposes, such as counterillumination , aposematism (Claes et al., 2013;Duchatelet et al., 2019b), and conspecific recognition (Claes et al., 2014a. Luminescence is achieved via thousands of photophores located within the epidermis. Each photophore is composed of a cupshaped layer of pigmented cells encapsulating one to more than twelve photogenic cells (i.e., photocytes) and topped by one or more lens cells. In E. spinax, a guanine crystal reflector structure is located between the cup-shaped pigmented layer and the photocyte (Renwart et al., 2014(Renwart et al., , 2015. Photophores also display an iris-like structure (ILS), composed mainly of chromatophores, between the photocytes and the lens cells (Renwart et al., 2014;Duchatelet et al., 2020b). Recently, studies of the luminous system of E. spinax failed to identify the reactive compounds underlying the emission of light (i.e., luciferin/luciferase or photoprotein) (Renwart and Mallefet, 2013). Moreover, it has been demonstrated that shark luminescence is not due to symbiotic luminous bacteria (Duchatelet et al., 2019a). Therefore, the nature of the shark luminous system remains enigmatic.
In Metazoans, sharks are the only known bioluminescent organisms to hormonally control light emission. For the studied species, researchers have demonstrated the involvement of several hormones in the control of light emission: melatonin (MT) triggers light production, while alpha-melanocytestimulating (α-MSH) and adrenocorticotropic hormones (ACTH) inhibit it (Claes and Mallefet, 2009c;Duchatelet et al., 2020b). Prolactin triggers brighter and faster light emission than MT in Etmopteridae (Claes and Mallefet, 2009c;Claes et al., 2011b), while this hormone inhibits light production in S. aliae (Claes et al., 2012). More recently, in silico mRNA sequences and expression sites of MT and α-MSH/ACTH receptors were highlighted within the photophores, but neither mRNA sequences nor protein presence was found for the prolactin receptor (Duchatelet et al., 2020a). Other molecules, such as nitric oxide or γ-aminobutyric acid, also exhibited modulatory effects on light emission in some Etmopteridae (Claes et al., , 2011a. Finally, an extraocular opsin (Es-Opn3) has been demonstrated to be involved in a secondary control targeting the ILS and modulating the aperture of this pigmented structure acting as a light organ shutter (Duchatelet et al., 2020c). To establish the conservation of photophore morphology and the control of hormonal light emission in the evolution of luminous Squaliformes, increasing the knowledge on bioluminescent sharks is crucial.
While the majority of Squaliformes never reach more than 60 cm in adulthood, the kitefin shark (also named seal shark or black shark), Dalatias licha (Bonnaterre, 1788), can grow to 180 cm (Compagno, 1984;Roberts et al., 2015). This giant holobenthic dalatiid has a worldwide distribution at depths ranging from 50 to 1800 m but it is usually found in depths below 300 m (Compagno, 1984;Roberts et al., 2015). Recently, through baited-remote video and muscle enzymatic activity analysis, D. licha was suggested to be one of the slowest moving elasmobranch species (Pinte et al., 2020). Reif (1985) assumed that this shark is luminous as it presents pavement-like placoid scales at the ventral side of the body like the related cookie cutter shark, Isistius brasiliensis (Quoy and Gaimard, 1824) (Reif, 1985;Widder, 1998;Delroisse et al., 2021). Nevertheless, no clear evidence has been put forward to confirm its luminescence status.
Here, organization, morphology, density, and physiological control of kitefin shark photophores were investigated. To determine if this species displays the same photophore structure and hormonal control, a comparative analysis was performed on the two most abundant New Zealand lanternshark species, E. lucifer and E. granulosus. Results are compared to previously studied dalatiids and etmopterids. Homogeneity of light emission control among luminous elasmobranch and photophore structures among each shark families are observed, strengthening a conservative evolution of light emission capabilities among sharks. These observations and results raise questions on the luminescence role for the largest luminous vertebrate. The use of counterillumination for this giant luminous shark is here suggested to be co-opted for a camouflage-type approach as a predatory tool.

Specimen Sampling
Shark specimens were captured during the Chatham Rise Trawl survey by the R.V. Tangaroa in January 2020 off the coast of eastern New Zealand. The survey used the same eight-seam hoki bottom trawl and survey methodology that was used on previous surveys (Hurst et al., 1992;Stevens et al., 2018). The net has 100 m sweeps, 50 m bridles, 12 m backstrops, 58.8 m groundrope, 45 m headline, and 60 mm codend mesh. The trawl doors were Super Vee type with an area of 6.1 m 2 .
A total of 37 D. licha [40.9-138.0 cm total length (TL)], 304 E. lucifer (16.2-53.2 cm TL), and 281 E. granulosus (19.3-75.6 cm TL) were captured on the survey, of which 13 D. licha, 7 E. lucifer, and 4 E. granulosus were used for bioluminescence studies. Each specimen was maintained in a tank with fresh cold sea water in a dark cold room until manipulation. Each shark was sexed, measured, weighed (Supplementary Table 1) and photographed in dim daylight and in dark conditions using Sony α7SII camera before having a full incision of the spinal cord at the level of the first vertebrae, according to the European regulation for animal research handling. Ventral skin of a specimen of S. aliae and I. brasiliensis, collected, respectively, as in Delroisse et al. (2021) and Duchatelet et al. (2020b), were used for dalatiid comparative photophore histology.

Photophore Histology and Density
Skin patches of 3 cm 2 were dissected from different locations along the body of D. licha specimens (i.e., rostral, mandibular, pecto-ventral, pectoral, ventral, dorsal, dorsal fin, pelvic, flank, infra-caudal, precaudal, and caudal zones; Figure 1A) to assess photophore presence, size and densities. Skin patches were fixed in 4% formalin at least overnight before being transferred to phosphate buffer saline (PBS). Skin patches were observed and photographed under a transmitted light microscope (Leitz Diaplan, Germany) coupled with a ToupCam camera (UCMOS Series C-mount USB2.0 CMOS camera, ToupTek, Zhejiang, China). Photophore densities (per mm 2 ) and mean diameter (n = 30 or 50 zones per species) were also measured on the two etmopterid species using the same protocol (Supplementary Figure 1).
In parallel, skin patches of D. licha, E. lucifer, E. granulosus, S. aliae, and I. brasiliensis were used to perform histological sections across the photogenic organ. Skin tissues were bathed for 7 days in decalcifying solution (OsteoRAL, Fast decalcifier for Large Anatomical Specimens, RAL Diagnostics, France) with constant agitation and renewal of the solution every 2 days, rinsed in PBS, and placed in PBS with increasing concentrations of sucrose (10% for 1 h, 20% for 1 h, and 30% overnight). Tissues were then embedded in optimal cutting temperature compound (O.C.T. compound, Tissue-Tek, Netherlands) and rapidly frozen at −80 • C. Sections of 10 µm were obtained with a cryostat microtome (CM3050S, Leica, Solms, Germany). Sections were placed on coated Superfrost slides (Thermo Scientific) and left overnight to dry. All sections were observed under a transmitted light microscope (Leitz Diaplan) equipped with a ToupCam camera (ToupTek).
Hormones known to trigger or inhibit light emission in luminous elasmobranchs were applied (Claes and Mallefet, 2009c;Duchatelet et al., 2020b). Here, evaluations of the effect of MT, α-MSH and ACTH were conducted for the first time on the dalatiid species, D. licha, and the etmopterid species, E. lucifer, and E. granulosus.
Experiments were first conducted on 10 D. licha specimens. To obtain a dose response curve for MT application, three different concentrations of MT (i.e., 10 −6 , 10 −7 , 10 −8 mol L −1 ) were used. Skin patches were immersed in 200 µL of MT solution (either 10 −6 , 10 −7 , 10 −8 mol L −1 ). To analyze the effect of α-MSH and ACTH on the light emission of D. licha, another set of skin patches were subjected to an immersion in 100 µL of MT 10 −6 mol L −1 followed after 5 min by an application of 100 µL of either α-MSH 10 −6 mol L −1 or ACTH 10 −5 mol L −1 . Luminescence of ventral skin patches subjected to the various treatments was measured using a FB12 tube-luminometer (Titertek-Berthold, Pforzheim, Germany) calibrated as in Duchatelet et al., 2020b. Lights emissions were recorded through FB12-Sirius, multiple kinetics software (Titertek-Berthold) for at least 30 min with a measurement every 58 s. For comparative purposes, similar treatments were performed on seven specimens of E. lucifer (same experiments) and four E. granulosus specimens (MT dose response and α-MSH treatments). In parallel, for D. licha and E. lucifer, photophore aperture and closure were observed after drug application by taking a time-lapse series of pictures (every 10 min) with a Sony α7SII camera mounted on a binocular microscope.
Luminescence measurements were characterized as follows (Duchatelet et al., 2020b): the maximum intensity of light emission [Lmax, in megaquanta per second (Mq s −1 )], the total amount of light emitted during experimentation [Ltot, in Gigaquanta per hour (Gq h −1 )] and the time to reach maximum light intensity [TLmax, in seconds (s)]. Inhibitory actions of α-MSH and ACTH were measured as the total amount of light emitted after the second drug application [Ltot app , in Gq h −1 ]. All light parameters were standardized according to the surface area of each skin patch (in cm 2 ). A second treatment (α-MSH or ACTH) was added to the first one when the light intensity plateau was reached with the MT application, each timing being speciesspecific. Results of the luminescence decrease were expressed as a percentage of the maximal luminescence value (i.e., plateau MT) measured before the second application. (B) Representation of the dorso-ventral photophore density gradient. Black-dotted photophores (red arrowhead) observed between the placoid scales (delimited areas) at the (C) rostral and (D) ventral areas. Rostral area presents specific leaf-shaped placoid scales, while ventral area harbors typical pavement type placoid scales. (E) Close-up of the black circular-shaped photophores within the integument surrounding the ventral placoid scales. d, placoid scale; e, epidermis; m, melanophore; p, photophore. (F) Photophore density variation across the studied zones. Different lettering indicates statistical differences. All density values are expressed as mean ± SEM.
To evaluate the putative evolutive conservation of the hormonal control of light emission in dalatiids and etmopterids, pharmacological data on shark luminescence were extracted from literature.

Statistical Analyses
All analyses were performed with the software R studio (version 1.1.383, 2009, R Studio Inc., United States). Variance normality and homoscedasticity assumptions were tested by Shapiro-Wilk and Levene's test, respectively, before running ANOVA which reveals significant differences between skin photophore densities or pharmacological treatments. When these parametric assumptions were not met, a non-parametric Kruskal-Wallis ANOVA was used. Post hoc Tukey's tests or Wilcoxon tests allowed pair-wised comparison of means, attributing different letters to significantly different values (P-value < 0.05).

Luminous Pattern and Photophore Morphology
A blue glow was observed on the ventral surface of D. licha, E. lucifer, and E. granulosus specimens kept in a fully dark environment (Figures 2A,D,E). D. licha also emit a faint blue glow from the lateral and dorsal areas and at the two dorsal fins (Figure 3 -Mallefet personal observation). Both etmopterids present a more complex pattern of light emission with flank marks, and lateral, dorsal, and rostral patterns (Figure 4; E. granulosus -Mallefet personal observation). Skin patches observed in toto present black round-shaped photophores distributed between placoid scales for all the observed sharks (Figures 1C-E and Supplementary Figure 1). The mean photophore diameters are 83.9 ± 9.5, 122.4 ± 10.8, and 132.3 ± 14.5 µm, for D. licha, E. lucifer, and E. granulosus, respectively. No statistical differences in photophore diameter were observed between the zones presenting large amount of photophores (ventral, pecto-ventral and infra-caudal) (D. licha Photophores histology highlights a single photocyte within small photophores in dalatiids, while etmopterids harbor bigger and more complex photophores. c, connective tissue; e, epidermis; i, iris-like structure cells; l, lens cell; p, photocyte; s, pigmented sheath. In toto shark picture scale bar: 10 cm; photophore section scale bar: 100 µm. Frontiers in Marine Science | www.frontiersin.org Analyses of photophore density along the D. licha body show an increasing dorso-ventral repartition of photophores reaching up to 20.14 ± 4.01 photophores per mm 2 at the ventral side of the shark (Figures 1A,B,F). The lowest densities were observed for the caudal and dorsal areas with a mean density of 2.85 ± 1.11 and 4.85 ± 1.70 photophores per mm 2 , respectively ( Figure 1A). Statistical differences [ANOVA: F(11,514) = 128.64, P-value < 2.2 × 10 −16 ] in photophore densities are illustrated Figure 1F and Supplementary Table 2A. The scales of the rostral area were leaf-like in shape while the scales of the remaining body parts were pavement-like in shape (Figures 1C,D).
For both studied etmopterids, a high density of photophores was observed at the pectoral zone with 34.00 ± 6.20, and 15.63 ± 2.50 photophores per mm 2 for E. lucifer and E. granulosus, respectively. E. granulosus also have a high density of photophores at the infra-caudal and caudal zones.
Histological sections across photogenic skin highlight the structure of D. licha photophores. Each light organ is embedded in the stratified squamous epidermis and is composed of a cupshaped pigmented sheath containing a unique photocyte, topped by a lens cell with a few diffuse pigmented cells between the photocyte and lens cell (Figure 2A). This structural organization is similar to that found in S. aliae and I. brasiliensis photophores (Figures 2B,C).
Photophore morphologies of E. lucifer and E. granulosus are consistent with those already described for other etmopterids (i.e., E. spinax, Etmopterus molleri, and Etmopterus splendidus) (Figures 2D-F). They are composed of a cup-shaped pigmented sheath embedded with luminous cells and topped a with lens. They are similar to dalatiid photophores, but they harbor a higher number of photocytes, a larger iris-like structure area, and more lens cells (up to 3) (Figures 2D-F).

Light Emission Control
The effect of MT on D. licha, E. lucifer, and E. granulosus was tested through a dose-dependent response. For the studied species, MT 10 −6 mol L −1 triggered a long-lasting light emission, significantly different from the MT 10 −8 mol L −1 application (P-value < 0.05; Figures 5A-C and Supplementary Tables 3A,  4), while MT 10 −7 mol L −1 triggered an intermediate light emission and Ltot value (Figures 5A-C and Supplementary  Tables 3A, 4). All treatments were significantly different from the shark saline control, except for the MT 10 −7 and 10 −8 mol L −1 treatments of E. granulosus (P-value < 0.05; Supplementary  Tables 3A, 4). Although the total amount of light emitted under MT 10 −6 mol L −1 treatment was significantly different [Kruskal-Wallis χ 2 (2) = 10.14, P-value = 0.0063], E. lucifer produced a mean total amount of light during the experiment 2.5 and 5 times higher than D. licha and E. granulosus, respectively. Similar patterns of bioluminescence were observed for the three species (Figures 2, 5).
The effect of α-MSH was evaluated for the three species after reaching the Lmax triggered through MT 10 −6 mol L −1 application. Application of α-MSH 10 −6 mol L −1 induced a rapid decrease of light emission for the studied luminous sharks (Figures 6A-C and Supplementary  Table 3B). After MT-induced bioluminescence, Ltot app values of α-MSH were statistically significant compared with the MT 10 −6 mol L −1 control (Figures 6A-C and  Supplementary Tables 3B, 5).
The effect of ACTH was evaluated on D. licha and E. lucifer bioluminescence. Similar to the results obtained for α-MSH, ACTH 10 −5 mol L −1 applications rapidly induced a decrease in light emission (Figures 6A,B and Supplementary Table 3B).  Each Ltot app value of ACTH 10 −5 mol L −1 was not significantly different from those of α-MSH 10 −6 mol L −1 , respectively, but were statistically different from the MT 10 −6 mol L −1 control (Figures 6A,B and Supplementary Tables 3B, 5). Mean values of Lmax, TLmax, Ltot, Ltot app are presented in Supplementary Table 3.
The time-course of light emission in D. licha under MT 10 −6 mol L −1 stimulation revealed a concomitant opening of the photophore ILS within 15 min of luminescence, in which the ILS stayed open for the next 30 min while the light level remained high (Figure 7A). In the case of E. lucifer MTinduced luminescence, a rapid opening of the photophore ILS was observed within 8 min followed by a slow decrease during which a closure of the ILS was visible (Figure 7B). Aperture and closure of photophores showed pigment movements concomitant to light emission.

DISCUSSION
The three studied shark species inhabit the mesopelagic zone (Roberts et al., 2015), therefore they face an environment with no place to hide, hence the need for glowing camouflage or counterillumination, first proposed by Clarke, 1963. The mesopelagic zone, often called the twilight zone, ranges from 200 to 1000 m depth (maximal depth of solar light penetration) and is the realm of bioluminescence (Martini and Haddock, 2017;Martini et al., 2019). At 200 m the residual solar light is considered too weak to initiate photosynthesis but organisms living there are well adapted to see in low light conditions (Nicol, 1978). Mesopelagic cephalopods, sharks and bony fishes have large eyes with specialized structures such as a large iris, a tapetum, huge rod density, high content of opsins (rhodopsin and chrysopsin), and an elevated integration rate at the optical nerve which allows them to perceive very low light levels down to 800 m depth (Douglas et al., 1998;Warrant, 2004;Warrant and Locket, 2004;Claes et al., 2014a,b).

Luminescent Pattern
The light emission pattern observed in D. licha is similar to that observed in previously studied dalatiids i.e., S. aliae and I. brasiliensis (Claes et al., 2012;Delroisse et al., 2021). The dorso-ventral gradient and the relative homogeneity in ventral photophore densities suggest the luminescence is used for counterillumination. The luminous pelvic zone of D. licha reveals a sexual dimorphism but, contrary to E. spinax and Etmopterus molleri Duchatelet et al., 2020c), it is not brighter than the rest of the ventral body, suggesting it is less important for sexual signaling. The kitefin shark D. licha, like other dalatiids, does not have flank marking or specific dorsal patterns. The lack of these luminescent patterns, previously suggested to be used as conspecific signaling for group aggregation, swimming, or hunting in etmopterids Duchatelet et al., 2019b), rules out this function in D. licha. The aposematic function described for etmopterids (Claes et al., 2013;Duchatelet et al., 2019b) is also ruled out for D. licha luminescence due to the absence of dorsal fin defensive spines. Nevertheless, D. licha is the first shark with fully luminous dorsal fins (Figures 1A, 3), which raises questions about its luminescence function.
The light emission patterns of E. lucifer and E. granulosus, are similar to that of previously studied etmopterids. The dorsal photophores, flank markings, and brighter pectoral fin and claspers are likely to be used for intraspecific communications while the ventrally emitted light is likely to be used for counterillumination. These functions have been documented for E. spinax (Claes and Mallefet, 2009a;, Etmopterus molleri , and Etmopterus splendidus (Claes et al., 2011b). However, a bioluminescence aposematic function through specific spineassociated photophores (Claes et al., 2013;Duchatelet et al., 2019b) was not documented for E. lucifer and E. granulosus. Reif (1985) postulated that a trade-off exists between the space occupied by placoid scales and luminous organs, and that four different types of placoid scales have evolved to allow this trade-off: pavement, cross-, bristle/ needle-, and hook-shaped placoid scales. A new type of squamation with overlapping leafshaped placoid scales is present in the luminous rostral area of D. licha. This new bioluminescent-associated squamation was observed in the somniosid, Zameus squamulosus, which is assumed to be luminous . This new type of bioluminescence-associated placoid scale needs to be highly translucent or possess specific physical characteristics to allow efficient light transmission. The use of Reif placoid scale types to assess the bioluminescent status of a shark species is not a decisive character as shown by a recent study of Ferrón et al. (2018); highlighting the presence of bioluminescent-like squamation in a galeomorph shark, Apristurus ampliceps, a species not known to be luminous.

Luminescence Control Evolutive Conservation
The effect of hormones on light emission in D. licha, E. lucifer and E. granulosus are consistent with increasing literature on light emission control in sharks (Claes and Mallefet, 2009c;Claes et al., 2012;Duchatelet et al., 2020b,d): MT, and α-MSH/ACTH, have been demonstrated as the main triggering and inhibiting agents of shark luminescence, respectively. Similar to observations of E. spinax and Etmopterus molleri photophores Duchatelet et al., 2020b), aperture and closure of D. licha and E. lucifer photophores involved pigment motion within the ILS cells. Simultaneities of curves kinetics and pigment motions highlight the evolutive conservation of hormonally controlled pigment motion regulating luminescence. These data strongly suggest that luminous etmopterids and dalatiids share a common luminescence control mechanism, involving at least MT, and α-MSH/ACTH hormones. This control is assumed to have been successfully and evolutionary co-opted from shark melanophore pigment motion control by a common ancestor of these two squaliform families. For both families, luminescence appears to be dually controlled at the level of (i) the photocyte, site of luminescent reaction, and (ii) the ILS cells, acting as a diaphragm capable of occluding light produced by the photocytes, via melanophore-associated pigment movements (Duchatelet et al., 2020b,d). This was recently demonstrated within ILS cells of the lanternshark, E. spinax (Duchatelet et al., 2020d) i.e., transduction pathways that activate cellular motors such as dynein and kinesin, leading pigment movements within ILS melanophores. The bioluminescence control mechanisms in the two studied etmopterids, as well as in D. licha, might share common features. Moreover, the involvement of extraocular photoreception events in the light emission control of photophores (Duchatelet et al., 2020d), remains to be deciphered for these sharks. Further research are necessary to fully demonstrate the evolutive conservation of luminescence control within etmopterids and dalatiids.

Luminescence of Dalatias licha
The question remains concerning bioluminescence in the largest luminous vertebrate; why does D. licha emit light ventrally to counterilluminate when it has few or no predators? Pinte et al. (2020), analyzed the swimming speed of several New Zealand deep-sea sharks, and found that D. licha possesses one of the slowest cruise swimming speeds ever measured in sharks. Conversely, this species is assumed to possess a high burst capability (Pinte et al., 2020). Stomach content analyses have revealed that this shark species hunts and eats etmopterids, which have a higher cruise swimming speed. Therefore, there are two hypotheses which might explain the ventral luminescence of this holobenthic species: luminescence might be used (i) to illuminate the ocean floor while searching and hunting for prey; or (ii) to stealthily approach toward prey, using counterillumination camouflage, before striking fast when close enough (Zintzen et al., 2011), allowing them to predate etmopterids. In both cases, the principle of counterillumination would have been distorted to serve as a predation tool instead of an avoidance mechanism, a hypothesis already proposed for the cookie cutter shark, I. brasiliensis (Widder, 1998). However, to validate such hypotheses for these dalatiid species, in vivo observations and behavioral studies are essential.

CONCLUSION
Through a histological and pharmacological approach, the bioluminescence of three different shark species was investigated. Our results support evolutive conservation of light organ morphology and luminescence control. For the first time, luminescence was recorded and analyzed for the largest luminous vertebrate, D. licha and two lanternsharks, E. lucifer and E. granulosus. Dalatiid photophores are similar between species and are structurally composed of a single photocyte embedded in a cup-shaped pigmented cell and surmounted by lens cells. The same observation was made for etmopterids, which showed a conservation of photophore structure between species. Etmopterid photophores are slightly more complex than those of dalatiids, with several photocytes and a well-developed ILS between the lens cells and the photocytes. Through this study, the action of MT and α-MSH/ACTH in the bioluminescence control in these two families was shown to be identical and seem to have been co-opted during evolution from the regulation of skin pigment movements. With these data, we can assume that the common luminous ancestor of etmopterids and dalatiids likely had hormonal control of its luminescence and had luminous organs similar to those of the dalatiids (i.e., the simplest structure) for counterillumination.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

ETHICS STATEMENT
Ethical review and approval was not required for this study. The shark specimens were captured as bycatch of a fisheries assessment survey for the New Zealand Ministry for Primary Industries.

AUTHOR CONTRIBUTIONS
JM and DS collected the samples. JM collected the bioluminescence pictures and performed pharmacological studies and fixations on the survey. LD performed the classical histology, pattern, and pharmacological analyses. LD and JM were major contributors to the initial manuscript that was improved by DS revisions. All authors approved the final manuscript.