We report the finding of two partial specimens of Cryptolepas rhachianecti (Cirripedia, Coronulidae), a coronulid barnacle known only to inhabit the skin of gray whales (Eschrichtius robustus), in Pleistocene-aged sediments from the Canoa Basin, Ecuador. While the historical range of gray whales includes the North Pacific and North Atlantic, to our knowledge this is the first inferred evidence of a gray whale population having resided within the South Pacific. We describe the two Cryptolepas rhachianecti fossils, use isotopic analysis to investigate evidence of migration in their host whales, and discuss their implications for our understanding of gray whale evolutionary history.
Introduction
Coronuloidea consists of a superfamily of commensal barnacles adapted to live on a remarkable variety of organisms, including manatees, crabs, mollusks, sea snakes, sea turtles, and cetaceans (Hayashi, 2013). The whale-living barnacles, or coronulids (family Coronulidae), are the most recently derived lineage, having diverged from the other commensal barnacles within the past five million years (Hayashi et al., 2013). This timing coincides with the onset of gigantism in baleen whales (Slater et al., 2017), and most coronulids make their home on the modern giants, although some species are occasionally seen on delphinids (Seilacher, 2005; Hayashi, 2012).
As a group, the coronulids are relatively abundant throughout the modern oceans. The most abundant coronulid species in the oceans today is Coronula diadema (Linnaeus, 1767), which lives attached to the skin of humpback whales (Megaptera novaeangliae Borowski, 1781). Where whale barnacle fossils are found, they belong nearly exclusively to the genus Coronula (Fleming, 1959; Zullo, 1969; Beu, 1971; Buckeridge, 1983; Bianucci et al., 2006a, b; Dominici et al., 2011; Collareta et al., 2016; Taylor et al., 2019). The genus Cryptolepas is represented by a single extant species, Cryptolepas rhachianecti, which is a host-specific inhabitant of the skin of gray whales (Eschrichtius robustus Lilljeborg, 1861) (Newman and Ross, 1976; Newman and Abbott, 1980; Seilacher, 2005; Bradford et al., 2011; Hayashi, 2012; Swartz, 2018). Like its host, C. rhachianecti previously has been known only from the Northern Hemisphere, both in the present day and within the fossil record. While fossil whale barnacles usually exist in small numbers wherever they are found, an exception to this general rule is the Canoa Basin of Ecuador, where dozens of coronulid fossils have been collected, previously all assigned to the genus Coronula (Bianucci et al., 2006a, b).
Here we report finding two partial shells of the gray whale barnacle Cryptolepas rhachianecti in Pleistocene-aged deposits from the Canoa Basin, Ecuador. To our knowledge, these are the oldest known fossil specimens of Cryptolepas, the first occurrences of the genus in the Southern Hemisphere, and the first evidence of a gray whale population regularly living in or migrating through the southern Pacific Ocean.
Geological setting
Specimens were collected in Pleistocene-aged sediments of the Canoa Basin in the Cabo San Lorenzo area of the Manabi Province, Ecuador (1.1821°S, 80.8622°W). During the Plio-Pleistocene, the Cabo San Lorenzo area was an uplifting island separated from mainland Ecuador by a shallow, narrow strait. Continued uplift eventually connected this island to mainland Ecuador via a land bridge. Glacial-interglacial oscillations caused the shoreline to advance and retreat multiple times, extending west of present-day La Plata Island during glacial periods, when sea level was lower (Di Celma et al., 2002, 2005; Cantalamessa and Di Celma, 2004). During warmer interglacials, rising sea levels resulted in the shoreline retreating eastward to create a southwesterly facing coastal embayment, the Canoa Basin (Fig. 1). Over time, the basin has been progressively filled by ∼120 m of cyclically stacked shallow marine strata (Di Celma et al., 2002, 2005).
Figure 1.
(1) A map showing the location of the Canoa Basin sediments in relation to the range of modern-day gray whales, with inset showing greater detail of the region surrounding the field site. (2) During glacial lowstands, the Canoa coastline extended west past present-day Isla de la Plata. (3) During past interglacial high-stands, the coastline moved east, creating a small coastal embayment, the Canoa Basin.

Stratigraphic, sedimentological, paleoecological, and taphonomic studies of the entire succession indicate it is composed of a series of glacio-eustatically driven parasequences capturing recurring periods of rising and falling sea levels, exposed along ∼10 km of coastline between Punta Canoa and Rio Callejón (Di Celma et al., 2002, 2005). The succession has been divided into the lower Canoa and upper Tablazo formations, with the Canoa Formation being further divided by an angular unconformity (Sheppard, 1930; Pilsbry and Olsson, 1941; Marchant, 1961; Savoyat, 1971; Baldock, 1982; Tsuchi et al., 1988; Whittaker, 1988; Di Celma et al., 2005). Older sediments in the basin have yielded Oligo-Miocene fossils, and the Miocene-aged Tosagua Formation lies just below the lower Canoa (Cadena et al., 2018). The lower Canoa rests unconformably atop the Tosagua, and is composed of fine-grained, bluish-gray, silty to sandy shales (Di Celma et al., 2005). Above, the upper Canoa consists of gray-brown silts to sands, overlain by the brown sands of the Tablazo Formation. Sediments in both formations are rich in fossils (bivalves, gastropods, echinoids, crustaceans, corals, and occasional vertebrate remains) (Bianucci et al., 2006a; Flores et al., 2018).
The lower Canoa consists of four depositional sequences formed under the control of 40 kyr-long sea level fluctuations (Di Celma et al., 2002, 2005; Cantalamessa and Di Celma, 2004). Each sequence is bound below by an erosional surface, above which lies a basal hiatal shell bed followed by several meters of sparsely fossiliferous sediment (Di Celma et al., 2005). While the lowermost shell bed is quite taxa-rich and densely packed, whale barnacles are rare in the lowest two sequences and become more abundant thereafter. The upper Canoa and Tablazo formations are interpreted to be a continuous sedimentary record of two and six depositional sequences created under the control of 100 kyr-long sea level fluctuations (Di Celma et al., 2005). Shell beds are found basally and mid-cycle, separated by sparsely fossiliferous siliciclastics, and exhibit within-habitat time averaging (Bianucci et al., 2006a). The Tablazo Formation was described and is best known from the Santa Elena Province (∼120 km to the south of this site), and has yielded mostly continental fossils (Hoffstetter, 1952; Edmund, 1965; Ficcarelli et al., 2003; Lindsey and Lopez, 2015; Cadena et al., 2017), but also marine remains (Edmund, 1965; Flores, 2018; Flores et al., 2020).
The Canoa and Tablazo formations previously have yielded dozens of whale barnacles belonging to the species Coronula diadema, which lives on the skin of humpback whales, leading to the conclusion that the Canoa Basin was visited by migrating whales in the Pleistocene (Bianucci et al., 2006a, b).
Materials and methods
Fossils were collected in September 2018 from sediments of the upper Canoa Formation in the Canoa Basin, Ecuador, coordinates 1.1821°S, 80.8622°W. For isotopic analysis, a small Dremel handheld drill was used to collect calcite samples of 50–100 micrograms from along the primary (vertical) growth axis of the shell. Samples were analyzed at the Center for Stable Isotope Biogeochemistry at the University of California, Berkeley with a GV IsoPrime mass spectrometer with Dual-Inlet and MultiCarb systems. Several replicates of one international standard NBS19 and two lab standards CaCO3-I and II were measured along with every run of samples. Overall external analytical precision is ±0.07‰ for δ18O and ±0.05‰ for δ13C.
Barnacle calcite δ18O is determined by both the temperature and δ18O of the seawater in which it forms, as described by the balanomorph-barnacle paleotemperature equation of Killingley and Newman (1982):
where δC denotes barnacle calcite δ18O and δW denotes seawater δ18O. Because of latitudinal differences in the whale's feeding and breeding regions, the barnacle experiences the coldest waters in the summer feeding season. Although this will also generally correspond with the lowest seawater δ18O, the temperature-dependent fractionation in barnacles (and other calcifying organisms) generates an enriched shell calcite δ18O in cold temperatures and a depleted δ18O in warm temperatures.
Repositories and institutional abbreviations.—Specimens collected and described in this study are deposited at the University of California Museum of Paleontology (UCMP), Berkeley, California, and in the California Academy of Sciences (CAS), San Francisco, California.
Systematic paleontology
Class Maxillopoda Dahl, 1956
Subclass Cirripedia Burmeister, 1834
Superorder Thoracica Darwin, 1854
Order Sessilia Lamarck, 1818
Suborder Balanomorpha Pilsbry, 1916
Superfamily Coronuloidea Leach, 1817
Family Coronulidae Leach, 1817
Genus Cryptolepas Dall, 1872
Type species.—Cryptolepas rhachianecti Dall, 1872, from beached whale, Monterey, CA, U.S.A., by original designation.
Other species.—Cryptolepas murata Zullo, 1961, from Pleistocene deposits of San Quentin Bay, CA, U.S.A., by original designation.
Cryptolepas rhachianecti Dall, 1872
Figure 2
Holotype.—Shell collected from beached whale, Monterey, CA, U.S.A., (Dall, 1872, p. 300).
Description.—Specimen UCMP 116131 is a single compartment, with the radial lamellae mostly missing. Specimen UCMP 116132 is also represented by a single compartment, but is more complete, with many of the radial lamellae still intact (Fig. 2). Taken together, several characters identify the specimens as Cryptolepas rhachianecti; these include a transversely grooved sheath, the presence of four lamellar folds plus two half (sutural) folds on each compartment, irregularly branching radial lamellae where many branches fail to reach the periphery, and the absence of terminal flanges uniting the radial lamellae into a solid outer wall (Pilsbry, 1916; Davis, 1972). Transverse grooving of the sheath is present in both UCMP 116131 and UCMP 116132. The eroded lamellar folds are clear on the outward-facing walls of UCMP 116131. In UCMP 116132, many of the lamellar folds are well-preserved and display the species' characteristic irregular branching pattern, the failure of several of these branches to reach the periphery, and the absence of the wall-forming terminal flanges found in some other coronulids (Fig. 2). UCMP 116131 bears a strong resemblance to a (sub)fossil Cryptolepas rhachianecti specimen from the Netherlands reported by Bosselaers and Collareta (2016), and closely resembles a worn shell depicted by Pilsbry (1916, pl. 66, fig. 2).
Remarks.—Specimens UCMP 116131 and UCMP 116132 were both collected from the upper Canoa Formation. The features described above distinguish Cryptolepas from Coronula—a coronulid found in much greater abundance in the Canoa deposits. Coronula has an ungrooved sheath, fewer lamellar folds, fewer branches arising from the folds (all of which reach the periphery), and a solid outer wall formed by the uniting of T-shaped terminal flanges of the folds.
While fossils of the genus Cryptolepas are sparse, two species are now recognized in the fossil record. The first is C. rhachianecti, reported here and as a (sub)fossil from the Netherlands (Bosselaers and Collareta, 2016). The second is C. murata, reported by Zullo (1961) from late Pleistocene deposits of California. Cryptolepas traditionally has been suggested to have derived from the lineage that includes Coronula and Cetopirus, with the major differences in shell morphology seen in Cryptolepas interpreted as degenerative changes related to the more embedded, protected station of the shell in the host's skin (Pilsbry, 1916; Monroe, 1981). A primary feature distinguishing Cryptolepas rhachianecti from Coronula and Cetopirus is the absence of a complete outer wall to the shell, which in the latter two genera is formed by the merging of T-shaped flanges at the ends of the lamellar folds and which serves to create coring chambers that envelop prongs of the host's skin. Zullo (1961) interpreted C. murata as being an intermediate form, possessing many of the distinguishing features of C. rhachianecti while retaining an outer wall to the shell (Zullo, 1961, 1969). Bosselaers and Collareta (2016) have questioned whether Cryptolepas murata should be reassigned to the genus Cetopirus, but the rarity of Cryptolepas fossils makes answering that question difficult. While C. murata shares with Cetopirus an outer wall formed by terminal flanges of the lamellar folds, it also shares with C. rhachianecti a grooved sheath, less symmetrical internal branching of the lamellar folds, and the presence of internal branches that fail to reach the periphery of the shell (Fig. 2).
Figure 2.
Fossil Cryptolepas rhachianecti shells UCMP 116131 (1, 2) and UCMP 116132 (3, 4) from Pleistocene sediments of the Canoa Basin, Ecuador, alongside modern C. rhachianecti shell UCMP 34678 (5–7) and fossil C. murata shell UCMP 34677 (8–10). Fragile folds of shell are supported by whale skin in life (dark material seen in 5, 7) and will easily crumble when dislodged; the folds of UCMP 116132 survived by support from fine-grained sediment (3, 4). The grooved sheath (2, 4, 6, 9) can be seen in all specimens, and several display the blind-ended folds, which do not reach the periphery (4, 7, 10). In C. rhachianecti, the terminal ends of the primary folds do not connect to neighboring folds (4, 7), while in C. murata, the terminal ends of the folds fuse to form a rudimentary outer wall (8, 10), similar to that of the genus Coronula.

Results
Oxygen isotope ratios of coronulid shells have been shown to record the movements of their host whales (Killingley, 1980; Collaretta et al., 2018; Taylor et al., 2019). Isotopic analysis of fossil Cryptolepas rhachianecti shells UCMP 116131 and UCMP 116132 revealed δ18O ranges of 1.93‰ and 2.14‰ (Fig. 3), which is less than the 2.61–3.28‰ range measured in a modern-day C. rhachianecti shells of similar size (Fig. 3; Taylor et al., 2019, fig. S2). The fossils also have lower δ18O maxima than their modern counterparts. Fossil shell δ13C is likewise depleted relative to that of modern shells, and δ13C is only loosely correlated with δ18O (Fig. 3; Taylor et al., 2019, fig. S2).
Discussion
To our knowledge, these specimens represent the oldest known occurrence of Cryptolepas rhachianecti, and the first report of the species in the southern hemisphere. The oldest previously reported specimen of C. rhachianecti is a single compartment from late Quaternary deposits in the Netherlands (Bosselaers and Collareta, 2016). The gray whale lineage has a geological and historical presence in much of the North Atlantic, and it is reasonable to expect that C. rhachianecti also once ranged throughout the North Atlantic as well (Bisconti and Varola, 2006; Noakes et al., 2013; Alter et al., 2015; Bosselaers and Collareta, 2016; Hufthammer et al., 2018; Rodrigues et al., 2018). Overharvesting led to the final collapse of this population by sometime in the 18th century, however, and today C. rhachianecti is found only in the North Pacific, where it is a host-specific symbiont of the gray whale (Newman and Ross, 1976; Newman and Abbott, 1980; Scarff, 1986; Bradford et al., 2011; Hayashi, 2012). Modern gray whales live in two distinct populations in the western and eastern North Pacific, although some interchange between these populations does occur (LeDuc et al., 2002; Cooke et al., 2007; IWC, 2011). While the fossil record clearly attests to the presence of the gray whale lineage in the North Pacific, there are no historical or fossil records of gray whale populations living in the equatorial Pacific or southern hemisphere (Tsai et al., 2014; Tsai and Boessenecker, 2015; Kimura et al., 2018).
Besides the Netherlands specimen, no other fossils of C. rhachianecti are known, most likely due to the fragile nature of the shell. Due to this fragility, much of the fine structure of UCMP 116132 was damaged during isotopic sampling (images in Fig. 2 were taken prior to sampling). Whereas other whale barnacles extrude far above the host's skin and consequently are constructed robustly, C. rhachianecti lives deeply buried in the skin of gray whales, with the shell being largely supported by the skin of the host interweaving between the radial lamellae (Fig. 2). This embedded, low-profile lifestyle of C. rhachianecti may be related to the suction-feeding habit of gray whales, which causes the whales to rub their skin on the seafloor. If C. rhachianecti shells protruded much more above the skin, they would risk being dislodged or destroyed via abrasion.
The fragility of C. rhachianecti shells and their resulting scarcity in the fossil record suggests that the barnacles' host whales must have been a common visitor of the ancient Canoa Basin, considering we found two specimens in only three days of fieldwork. Whale barnacles are thought to be shed mostly while host whales are in their winter breeding areas, and accumulations of fossil whale barnacle shells have been interpreted as representing ancient whale breeding areas (Monroe, 1981; Bianucci et al., 2006a, b; Taylor et al., 2019). Our finding of Cryptolepas rhachianecti fossils raises the possibility that a Pleistocene gray whale population used the Canoa Basin as a winter breeding area.
Figure 3.
δ18O (1) and δ13C (2) profiles collected from along the primary growth axis of UCMP 116131, UCMP 116132, and CAS MAM 21149. CAS MAM 21149 is a modern C. rhachianecti shell collected from a gray whale that beached in northern California, likely while migrating south from its summer feeding areas. δ18O of barnacle shells is determined by the temperature and isotopic composition of the seawater in which the barnacle was located during each growth interval. Shell δ13C is affected by several factors and is best interpreted in conjunction with δ18O, where positive correlations may indicate changes in temperature or salinity, while depleted δ13C coupled with enriched δ18O may be a signal of upwelling. Analytical precision is ±0.07‰ for δ18O and ±0.05‰ for δ13C.

Extant gray whales in the eastern North Pacific spend their summer months feeding in cold, poleward waters, primarily in the Bering and Chukchi Seas, before migrating southward to warm, shallow, and sheltered waters along the coast of Baja California Sur to breed and raise their calves (Swartz et al., 2006; Mate et al., 2015). The Pleistocene Canoa Basin would have offered conditions similar to the lagoons and bays where modern gray whales overwinter, and the large number of Coronula diadema fossils found in the region suggests that the region at least served as an ancient breeding area for humpback whales (Bianucci et al., 2006a). While humpback and gray whales do not typically winter together, they do sometimes seek out similar water depths and temperatures (Martins et al., 2000), and wintering populations of both species can be found only a few miles apart along the coast of present-day Baja California.
If a Pleistocene gray whale population once inhabited the South Pacific, it has not survived to the present day. Humpback whales, on the other hand, still breed off the modern Ecuadorian coast. Compared to humpbacks, gray whales are less numerous, less ubiquitous, and more selective about areas where they will feed or breed. These differences may have made prehistoric gray whales comparatively rarer and more sensitive to change. While they are capable of generalist filter feeding like the other baleen whales, gray whales primarily prey on shallow, benthic invertebrate communities via suction feeding. This feeding habitat was greatly reduced during glacial maxima, however, reducing carrying capacity for the species (Pyenson and Lindberg, 2011). It is plausible that one of these disruptions could have winnowed a southern population beyond recovery, and that an intolerance for crossing warm equatorial waters coupled with the fidelity of whale mothers to their natal lagoons has kept the species from recolonizing these areas in the modern day (Lindberg, 1991; Goerlitz et al., 2003; IWC, 2011).
At least in the Canoa region, however, a primary cause of the loss of gray whales must have been the continual uplift of the Ecuador coast, which has eliminated the shallow embayments suitable for gray whale breeding. Gray whales only raise calves in sheltered lagoons and embayments, typically in waters of 10 m depth or less (Gardner and Chavez-Rosales, 2000; Goerlitz et al., 2003; IWC, 2011). While previous periods of high sea level would have created such a suitable breeding habitat within the Canoa Basin, uplift of the Ecuadorian coast has resulted in the modern Canoa coast being dominated by cliffs, while Andean uplift has further eliminated shallow embayments along the South American coast (Lindberg, 1991). Humpback whales, in contrast, will breed in a wider variety of water depths, and so may have been undisturbed by this loss of shallow habitat (Felix and Botero-Acosta, 2011; Pack et al., 2017).
An alternative explanation is that C. rhachianecti once commonly occurred on a different host whale species. However, there is little to support this hypothesis. Settlement of coronulid larvae seems to be initiated by chemical cues from the preferred host's skin (Nogata and Matsumura, 2005), while contact with other hosts initiates an immune response capable of shedding a barnacle (Ridgeway et al., 1997). Consequently, the most parsimonious explanation for the occurrence of C. rhachianecti fossils within the deposits is the prior presence of gray whales in the area.
The δ18O range of the Canoa Basin fossils is smaller than that collected from modern C. rhachianecti specimens, which may reflect the fragmentary nature of the shells, or it may reflect real differences in the migratory extent of the host whales as compared to their modern counterparts. If there was once a population of gray whales visiting the coast of Ecuador, then those whales have no directly comparable modern counterpart, but some insight may come from looking at the behavior of humpback whales that breed off the modern Ecuador coast (Sheidat et al., 2000; Felix and Botero-Acosta, 2011). These modern humpbacks are known to migrate to feeding areas along the coast of Chile, within the Magellan Strait, and along the Antarctic Peninsula (Gibbons et al., 2003; Acevedo et al., 2007; Capella et al., 2008; Felix and Olavarria, 2012; Hucke-Gaete et al., 2013). It is known that small numbers of modern gray whales do not migrate (Pyenson and Lindberg, 2011), but the δ18O ranges of these fossils are difficult to explain via annual variation in the immediate area, suggesting that the whales did migrate to some extent. The δ18O maxima are less enriched than would be expected if the host whales were traveling to Antarctica, but the δ18O ranges are similar to those of some barnacles collected from modern, migrating humpback whales (Taylor et al., 2019). It is plausible these ancient whales may have migrated to the Chilean coast, as some modern humpbacks do (Gibbons et al., 2003; Acevedo et al., 2007; Capella et al., 2008; Felix and Olavarria, 2012; Hucke-Gaete et al., 2013), but they also may have migrated to any number of unknown ancient feeding grounds.
Interpreting shell δ13C is less straight forward because δ13C may be affected by the incorporation of metabolic carbon, kinetic disequilibrium effects, temperature-dependent fractionation, and the effects of upwelling or freshwater input on δ13C of dissolved inorganic carbon. Interpretation of shell δ13C is best done in conjunction with corresponding shell δ18O, where depletion of both δ13C and δ18O may be a signal of freshwater input, and where depletion of δ13C corresponding to enriched δ18O may be a signal of upwelling (Killingley and Berger, 1979; Killingley and Lutcavage, 1983; Bemis and Geary, 1996; Sadler et al., 2012). There is some indication of depletion of δ13C where δ18O is most enriched in modern C. rhachianecti shells (Fig. 3; Taylor et al., 2019, fig. S2), but in both modern and fossil shells δ13C and δ18O are either not correlated or show a slight positive correlation over the majority of the captured signal, indicating that freshwater input and other factors may be exerting a strong influence on shell δ13C.
With only a handful of Cryptolepas fossils known, our finding of C. rhachianecti fossils in the Canoa Basin makes the region a promising area for further study. Collecting and isotopically analyzing more C. rhachianecti fossils from the Canoa Basin may offer more insight into the behavior of this prehistoric whale population.
Acknowledgments
Specimens were collected using excavation permit number M-JA-2018-82. Thank you to F. Flores and X. Guncay, our undergraduate research assistants during field work. We also thank D. Strauss and S. Kahanamoku of the UCMP for their work in photographing the specimens. We thank J. Jijón (Instituto Nacional de Patrimonio Cultural, Ecuador) for his support and help before, during, and after fieldwork. This study was supported by the Generalitat de Catalunya (CERCA Programme, and Beatriu de Pinós contract 2017 BP 00223 from AGAUR to JA), by the Geological Society of America, by the UC Berkeley Department of Integrative Biology, by the University of California Museum of Paleontology, by the American Society of Naturalists, and by the Paleontological Society.