# Patterns of Element Incorporation in Calcium Carbonate Biominerals Recapitulate Phylogeny for a Diverse Range of Marine Calcifiers

Robert Ulrich, Maxence Guillermic, Julia Campbell, Abbas Hakim, Rachel

Han, Shayleen Singh, Justin Stewart, Cristian Román-Palacios, Hannah

Carroll, Ilian de Corte, et al.

► To cite this version:

Robert Ulrich, Maxence Guillermic, Julia Campbell, Abbas Hakim, Rachel Han, et al. Patterns of Element Incorporation in Calcium Carbonate Biominerals Recapitulate Phylogeny for a Diverse Range of Marine Calcifiers. Frontiers in Earth Science, 2021, 9, 10.3389/feart.2021.641760. http://hal-03263910

HAL Id: hal-03263910

https://hal.univ-brest.fr/hal-03263910

Submitted on 18 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN ACCESS

**Edited by:**
Claire Rollion-Bard,
UMR 7154 Institut de Physique du
Globe de Paris (IPGP), France

**Reviewed by:**
Toshihiro Yoshimura,
Japan Agency for Marine-Earth
Science and Technology (JAMSTEC),
Japan
Dong Feng,
Shanghai Ocean University, China

*Correspondence:
Robert N. Ulrich
robu@g.ucla.edu
Robert A. Eagle
robeagle@ucla.edu

†ORCID:
Justin D. Stewart
orcid.org/0000-0002-7812-5095
Robert N. Ulrich
orcid.org/0000-0002-9541-2467
Robert A. Eagle
orcid.org/0000-0003-0304-6111
Hannah M. Carroll
orcid.org/0000-0003-3343-3358
Aradhna Tripati
orcid.org/0000-0002-1695-1754

**Specialty section:**
This article was submitted to
Geochemistry,
a section of the journal
Frontiers in Earth Science

Received: 14 December 2020

Accepted: 24 March 2021

Published: 04 May 2021

**Citation:**
Ulrich RN, Guillermic M,
Campbell J, Hakim A, Han R, Singh S,
Stewart JD, Román-Palacios C,
Carroll HM, De Corte I, Gilmore RE,
Doss W, Tripati A, Ries JB and
Eagle RA (2021) Patterns of Element
Incorporation in Calcium Carbonate
Biominerals Recapitulate Phylogeny
for a Diverse Range of Marine
Calcifiers. Front. Earth Sci. 9:641760.
doi: 10.3389/feart.2021.641760

# Patterns of Element Incorporation in Calcium Carbonate Biominerals Recapitulate Phylogeny for a Diverse Range of Marine Calcifiers

Robert N. Ulrich<sup>1,2,3*†</sup>, Maxence Guillermic<sup>1,2,3,4</sup>, Julia Campbell<sup>2,4</sup>, Abbas Hakim<sup>3,5,6</sup>, Rachel Han<sup>2,3,4</sup>, Shayleen Singh<sup>3,7</sup>, Justin D. Stewart<sup>8†</sup>, Cristian Román-Palacios<sup>3,9</sup>, Hannah M. Carroll<sup>1,3†</sup>, Ilian De Corte<sup>2,3,4,10</sup>, Rosaleen E. Gilmore<sup>4</sup>, Whitney Doss<sup>1,10</sup>, Aradhna Tripati<sup>1,2,3,4,10,11†</sup>, Justin B. Ries<sup>12</sup> and Robert A. Eagle<sup>2,3,4,6,10*†</sup>

Elemental ratios in biogenic marine calcium carbonates are widely used in geobiology, environmental science, and paleoenvironmental reconstructions. It is generally accepted that the elemental abundance of biogenic marine carbonates reflects a combination of the abundance of that ion in seawater, the physical properties of seawater, the mineralogy of the biomineral, and the pathways and mechanisms of biomineralization. Here we report measurements of a suite of nine elemental ratios (Li/Ca, B/Ca, Na/Ca, Mg/Ca, Zn/Ca, Sr/Ca, Cd/Ca, Ba/Ca, and U/Ca) in 18 species of benthic marine invertebrates spanning a range of biogenic carbonate polymorph mineralogies (low-Mg calcite, high-Mg calcite, aragonite, mixed mineralogy) and of phyla (including Mollusca, Echinodermata, Arthropoda, Annelida, Cnidaria, Chlorophyta, and Rhodophyta) cultured at a single temperature (25°C) and a range of pCO₂ treatments (ca. 409, 606, 903, and 2856 ppm). This dataset was used to explore various controls over elemental partitioning in biogenic marine carbonates, including species-level and biomineralization-pathway-level controls, the influence of internal pH regulation compared to external pH changes, and biocalcification responses to changes in seawater carbonate chemistry. The dataset also enables exploration of broad scale phylogenetic patterns of elemental partitioning across calcifying species, exhibiting high phylogenetic signals estimated from both uni- and multivariate analyses of the elemental ratio data (univariate: λ = 0–0.889; multivariate: λ = 0.895–0.99). Comparing partial R<fnref n="2" /> values returned from non-phylogenetic and phylogenetic regression analyses echo the importance of and show

<footnote n="1"><sup>1</sup> Department of Earth, Planetary, and Space Sciences, University of California, Los Angeles, Los Angeles, CA, United States,</footnote>

<footnote n="2"><sup>2</sup> Institute of the Environment and Sustainability, University of California, Los Angeles, Los Angeles, CA, United States,</footnote>

<footnote n="3"><sup>3</sup> Center for Diverse Leadership in Science, University of California, Los Angeles, Los Angeles, CA, United States,</footnote>

<footnote n="4"><sup>4</sup> Department of Atmospheric and Oceanic Sciences, University of California, Los Angeles, Los Angeles, CA, United States,</footnote>

<footnote n="5"><sup>5</sup> Department of Microbiology, Immunology, and Molecular Genetics, University of California, Los Angeles, Los Angeles, CA, United States, <sup>6</sup> Department of Ecology and Evolutionary Biology, University of California, Los Angeles, Los Angeles, CA, United States, <sup>7</sup> Department of Integrative Biology and Physiology, University of California, Los Angeles, Los Angeles, CA, United States, <sup>8</sup> Department of Ecological Science, Vrije Universiteit Amsterdam, Amsterdam, Netherlands, <sup>9</sup> Department of Ecology and Evolutionary Biology, The University of Arizona, Tucson, AZ, United States, <sup>10</sup> Institut Universitaire Européen de la Mer, Plouzané, France, <sup>11</sup> American Indian Studies Center, University of California, Los Angeles, Los Angeles, CA, United States, <sup>12</sup> Department of Marine and Environmental Sciences, Marine Science Center, Northeastern University, Boston, MA, United States</footnote>
that phylogeny explains the elemental ratio data 1.4–59 times better than mineralogy in five out of nine of the elements analyzed. Therefore, the strong associations between biomineral elemental chemistry and species relatedness suggests mechanistic controls over element incorporation rooted in the evolution of biomineralization mechanisms.

**Keywords:** marine calcification, calcite, aragonite, trace elements, ocean acidification, biomineralization, phylogeny

# INTRODUCTION

Calcareous biominerals can possess chemical compositions and morphologies as diverse as the organisms that form them (Lowenstam, 1981; Addadi et al., 2003; Weiner and Dove, 2003). Biominerals form within localized regions, isolated to varying degrees from the external environment. In the marine realm, where seawater is generally the source of the ions for biomineralization, the parent fluid for mineralization is chemically modified from seawater. It is thought that passive leakage of seawater and/or ions to the site of calcification, or active engulfment and vacuolization of seawater and/or transport of ions, may play a role in some organisms, such as corals and foraminifera (Erez, 2003; Gagnon et al., 2012; Venn et al., 2020). To control the composition of this fluid, biomineral constituents (typically cations) are transported into and out of the calcifying fluid by active (ATP-consuming) transport processes or by passive diffusion across membranes and chemical gradients (Simkiss et al., 1986). Organisms also vary widely with respect to the location of biomineral formation, with some such as coccolithophores and echinoderms beginning the process intracellularly, while others, such as molluscs, mineralizing extracellularly onto organic templates (Simkiss et al., 1986; Addadi et al., 2003; Weiner and Dove, 2003). The complex nature of biomineral chemistry, such as the inclusion of organic molecules and metastable mineral phases, has also been suggested to contribute to the variation in the dissolution kinetics of biogenic marine carbonates (Ries et al., 2016).

Marine calcifying organisms typically form either calcite or aragonite, or a mixture of the two calcium carbonate polymorphs. Elemental substitutions for Ca²⁺ in the CaCO₃ crystal lattice can be predicted by ionic radius and thermodynamic equilibrium (Brand and Veizer, 1980; Menadakis et al., 2009). Elements larger than Ca²⁺ are predicted to be discriminated against in the tight rhombohedral lattice of calcite relative to the orthorhombic lattice of aragonite (e.g., Shen et al., 1991; Gabitov et al., 2011, 2019; DeCarlo et al., 2015; Mavromatis et al., 2018). For calcite, there is also considerable variation in the abundance of Mg²⁺ cations substituting for Ca²⁺ in the crystal lattice. Magnesium is the most abundant secondary element incorporated into marine biogenic CaCO₃, and secular variation in the Mg/Ca ratio of seawater is hypothesized to be linked to shifts in the dominant CaCO₃ polymorph found in the geologic record of marine calcifiers (Chavé, 1954; Stanley and Hardie, 1999; Ries, 2010; Drake et al., 2020). Studies on the inorganic precipitation of CaCO₃ have shown that solution Mg²⁺/Ca²⁺ ratio influences CaCO₃ polymorph mineralogy, with molar Mg/Ca ratios > 2 favoring precipitation of high-Mg calcite and aragonite and

Mg/Ca < 2 favoring low-Mg calcite (Mucci and Morse, 1983; Park et al., 2008; De Choudens-Sánchez and González, 2009; Konrad et al., 2018; Purgstaller et al., 2019; Mergelsberg et al., 2020). Physiological regulation of [Mg²⁺] in the parent fluids for calcification has also been observed within some types of marine calcifiers. For example, *Mytilus edulis* have been observed to actively exclude Mg from their extrapallial fluid relative to seawater (e.g., Lorens and Bender, 1977). Furthermore, the extent of Mg incorporation into calcareous shells appears to influence the incorporation of other elements that substitute for Ca²⁺, such as Sr (Mucci and Morse, 1983). Additionally, some dissolved anions in seawater may substitute for CO₃²⁻ in the CaCO₃ lattice (e.g., SO₄²⁻; Burdett et al., 1989) or may be incorporated by occlusion in the crystal lattice (e.g., HCO₃⁻; Feng et al., 2006). Cation substitution for Ca²⁺ may also promote these latter two mechanisms of anion incorporation by distorting the crystal lattice (Mucci and Morse, 1983).

In addition to the paleobiological study of the evolution of marine biomineralization, a principal driver of the study of the chemistry of biogenic marine carbonates is the development of tracers for past changes in the ocean. The extensive literature on this topic notably includes the study of the incorporation of Mg and Sr as proxies for ocean temperature, the use of B and Li as proxies for the oceanic carbonate system, the use of Cd and Zn as links to nutrient availability, Ba as a tracer of freshwater input, Na as a tracer of salinity, and the use of Mo, Mn, and U as seawater redox proxies (e.g., Lorens and Bender, 1980; Boyle et al., 1981; Broecker and Peng, 1983; Delaney and Boyle, 1985; Delaney et al., 1989; Rosenthal et al., 1997; Elderfield and Rickaby, 2000; Marchitto, et al., 2002; Russell et al., 2004; Gaetani et al., 2011; Lea, 2013). In addition to these studies exploring relationships between elemental ratios and seawater physicochemical conditions of interest, some also explore drivers underlying so-called “vital effects,” or deviations in biogenic samples from apparent abiogenic equilibrium values (e.g., Weiner and Dove, 2003; Kunioka et al., 2006; Sinclair and Risk, 2006). Each element considered here has its own proposed mechanism of incorporation into CaCO₃, and we do not attempt to thoroughly review this large and evolving body of work.

Two notable examples are B and U incorporation, as both elements in marine carbonates can be used as tracers for the ocean pH and/or carbonate system (Yu et al., 2007, 2010; Tripati et al., 2009, 2011; Dawber and Tripati, 2012; DeCarlo et al., 2015; Haynes et al., 2019; Guillermic et al., 2020). U and B speciation in seawater are pH dependent; however, for the pH range of the present experiment, the relative abundances of uranium forms does not change (Endrizzi and Rao, 2014), and borate ion (thought to be the main form incorporated within the carbonate)
ranges from 24% (highest-pH treatment) to 6% of the total boron (lowest-pH treatment). According to the canonical systematics of B and U, both are incorporated into CaCO₃ by substitution of CO₃²⁻ within the crystal lattice, allowing a relatively simple relationship to be derived that relates B/Ca, U/Ca, and δ¹¹B in carbonates to ocean pH or DIC species (Hemming and Hanson, 1992; DeCarlo et al., 2015; McCulloch et al., 2017). Although the δ¹¹B-based proxy of seawater pH appears to work well in aragonite because aragonite incorporates most B as BO₄⁻, some recent studies—mainly on experimentally precipitated inorganic CaCO₃—suggest that the δ¹¹B pH proxy may be more complicated to interpret in calcite due to the apparent incorporation of some B as BO₃²⁻ in that polymorph (e.g., Branson, 2018; Farmer et al., 2019). As the δ¹¹B composition of BO₄⁻ and BO₃²⁻ in solution is markedly different, this could complicate the use of boron proxies in calcite shells in paleoceanography. In contrast, a recent study found that a diverse range of calcite-precipitating organisms, including high-Mg calcite-producing urchins, yielded δ¹¹B values incompatible with significant BO₄⁻ incorporation, with the possible exception of a coralline red alga species (Sutton et al., 2018; Liu et al., 2020). Furthermore, as many of the investigated species appeared to regulate the pH of the parent fluid for calcification, as indicated by independent measurements (e.g., pH microelectrodes, pH-sensitive dyes), deviations from ideal relationships between B shell chemistry (δ¹¹B, B/Ca) and seawater pH are not surprising (Sutton et al., 2018; Liu et al., 2020). These examples highlight that the diverse mechanisms used by marine calcifiers for biomineralization can cause deviations from relationships observed in inorganic precipitation experiments.

Attempts to model the mechanisms driving trace element
fractionation between CaCO₃ and solution have focused on
crystal growth rate and processes at the surface of the growing
crystal (e.g., Watson, 2004; DePaolo, 2011; Gabitov et al., 2014).
Several studies have highlighted the possible importance of
Rayleigh fractionation (e.g., Elderfield et al., 1996; Gaetani et al.,
2011). However, additional complexity is emerging with the
appreciation that many types of organisms, such as molluscs
and urchins, produce amorphous precursor phases from which
the biomineral ultimately forms via transient and/or metastable
phases (Addadi et al., 2003; De Yoreo et al., 2015). The
transformation from the amorphous to the crystalline phase
may also involve intermediate metastable CaCO₃ minerals,
such as vaterite (e.g., Jacob et al., 2017). Indeed, evidence
has emerged that [Mg²⁺], as well as other trace and minor
elements, may modulate the transformation of these precursor
or and/or metastable phases to crystalline phases (Loste et al.,
2003; Purgstaller et al., 2019; Mergelsberg et al., 2020), but
with different sensitivities for different trace/minor elements
and for different crystalline CaCO₃ end-products (Evans et al.,
2020). Furthermore, it is often assumed that trace elements
are hosted in the CaCO₃ crystal structure, but recent studies
of some organisms indicate that this is not always the
case. For example, Mg²⁺ in a coral skeletal aragonite was
reported to be present in a disordered non-crystalline phase of
potential organic affinity (Finch and Allison, 2008). Likewise,
a study on ostracods found that Mg²⁺ identified in the

shell was not incorporated into the calcite lattice of the shell
(Branson et al., 2018).

There is much yet unknown about the controls of trace and minor element incorporation in biominerals, and the literature is strongly focused on exploring this process in taxa that are typically used in palaeoceanographic reconstruction, such as foraminifera, corals, and molluscs. Additionally, previous studies tend to focus on only one to a few species and elemental ratios at a time and use species collected from different geographic locations, impeding direct comparisons between species and their respective geochemistry (e.g., Rosenheim et al., 2005; Gillikin et al., 2006; Dellinger et al., 2018; Piwoni-Piórewicz et al., 2021), emphasizing the need for systematic studies surveying geochemical data across diverse species. Here, we use univariate, multivariate, and phylogenetic analyses to explore a diverse array of 18 species of marine calcifying organisms, including Mollusca, Echinodermata, Arthropoda, Annelida, Chlorophyta, Rhodophyta, and Cnidaria, that were cultured in a controlled laboratory experiment using a single seawater source, a single temperature, and four controlled pCO₂ conditions (Ries, 2009). The calcification responses (Ries et al., 2009), polymorph mineralogies (Ries, 2011), and δ¹¹B compositions (Liu et al., 2020) of these cultured specimens has been previously characterized and reported. This sample set represents a unique opportunity to examine broad scale patterns of trace and minor element incorporation in a diverse array of ecologically and economically important marine calcifying organisms.

MATERIALS AND METHODS

Culturing Experiments

A detailed description of the culturing experiment and accompanying seawater carbonate chemistry measurements are described elsewhere (Ries et al., 2009), but we briefly summarize key details here. Eighteen species of calcifying marine organisms were cultured for 60 days at 25°C and at four controlled pCO₂ conditions (409, 606, 903, and 2856 ppm) (Table 1). Carbonate system parameters, including salinity, pH<sub>sw</sub>, dissolved inorganic carbon (DIC), total alkalinity, and calcite/aragonite saturation state were monitored and the net calcification rates of each specimen was calculated as the percent-change in buoyant weight of the specimen (Supplementary Tables 1, 2). Polymorph mineralogy of the shells (Table 1) was quantified via powder X-ray diffraction as described in Ries (2011).

Sample Selection

For this study, Li/Ca, B/Ca, Mg/Ca, Zn/Ca, Sr/Ca, Cd/Ca, Ba/Ca, U/Ca, and Na/Ca was quantified in the shells of the 18 calcifying species cultured under four pCO₂ conditions (409, 606, 903, and 2856 ppm), as described Ries et al. (2009). The organisms sampled include the American lobster (*Homarus americanus*), blue crab (*Callinectes sapidus*), gulf shrimp (*Penaeus plebejus*), conch (*Strombus alatus*), limpet (*Crepidula fornicata*), periwinkle (*Littorina littorea*), whelk (*Urosalpinx cinerea*), coralline red alga (*Neogoniolithon* sp.), Halimeda green alga (Halimeda incrassata), temperate purple urchin (*Arbacia punctulata*),
<table><caption>**TABLE 1** | Species examined in this study with supporting data.</caption><thead><tr><th>Organism</th><th>Latin</th><th>Phylum</th><th>Mineralogy<sup>#.21</sup></th><th>Net calcification rate response<sup>22</sup></th><th>pH<sub>cf</sub> (±SD)<sup>23</sup></th><th>ΔpH<sup>23</sup></th><th>ΔpH Response<sup>23</sup></th><th>Amorphous precursor?</th><th>Biomineralization mode</th></tr></thead><tbody><tr><td>American Lobster</td><td>Homarus americanus</td><td>Arthropoda</td><td>HMC</td><td>Threshold positive</td><td>—</td><td>—</td><td>—</td><td>Yes<sup>25,26,27</sup></td><td>Extracellular</td></tr><tr><td>Blue Crab</td><td>Callinectes sapidus</td><td>Arthropoda</td><td>HMC</td><td>Positive</td><td>7.965 (0.06)</td><td>-0.065</td><td>Neutral</td><td>Yes<sup>24,26,27</sup></td><td>Extracellular</td></tr><tr><td>Gulf Shrimp</td><td>Penaeus plebejus</td><td>Arthropoda</td><td>HMC</td><td>Positive</td><td>8.6 (0.12)</td><td>0.57</td><td>Positive</td><td>Yes<sup>26,27</sup></td><td>Extracellular</td></tr><tr><td>Blue Mussel</td><td>Mytilus edulis</td><td>Mollusca</td><td>97LMC:3A</td><td>Neutral</td><td>7.62 (0.36)</td><td>-0.53</td><td>Non-linear positive</td><td>Yes</td><td>Extracellular</td></tr><tr><td>Oyster</td><td>Crassostrea virginica</td><td>Mollusca</td><td>LMC</td><td>Negative</td><td>7.76 (0.38)</td><td>-0.39</td><td>Non-linear positive</td><td>Yes<sup>1,2,3</sup></td><td>Extracellular</td></tr><tr><td>Hard Clam</td><td>Mercenaria mercenaria</td><td>Mollusca</td><td>97A: 3HMC</td><td>Threshold negative</td><td>7.89 (0.30)</td><td>-0.22</td><td>Positive</td><td>Yes<sup>1,2,3</sup></td><td>Extracellular</td></tr><tr><td>Soft Clam</td><td>Mya arenaria</td><td>Mollusca</td><td>97A:3HMC</td><td>Negative</td><td>—</td><td>—</td><td>—</td><td>Yes<sup>4,5,6</sup></td><td>Extracellular</td></tr><tr><td>Bay Scallop</td><td>Argopecten irradians</td><td>Mollusca</td><td>97LMC:3A</td><td>Negative</td><td>—</td><td>—</td><td>—</td><td>Yes<sup>7</sup></td><td>Extracellular</td></tr><tr><td>Conch</td><td>Strombus alatus</td><td>Mollusca</td><td>97A:3LMC</td><td>Threshold negative</td><td>—</td><td>—</td><td>—</td><td>Yes<sup>2,8</sup></td><td>Extracellular</td></tr><tr><td>Periwinkle</td><td>Littorina littorea</td><td>Mollusca</td><td>97LMC:3A</td><td>Negative</td><td>—</td><td>—</td><td>—</td><td>Yes<sup>2,8</sup></td><td>Extracellular</td></tr><tr><td>Whelk</td><td>Urosalpinx cinerea</td><td>Mollusca</td><td>75A:25LMC</td><td>Negative</td><td>—</td><td>—</td><td>—</td><td>Yes<sup>2,8</sup></td><td>Extracellular</td></tr><tr><td>Limpet</td><td>Crepidula fornicata</td><td>Mollusca</td><td>97A:3LMC</td><td>Parabolic</td><td>—</td><td>—</td><td>—</td><td>Yes<sup>9</sup></td><td>Extracellular</td></tr><tr><td>Purple Urchin</td><td>Arbacia punctulata</td><td>Echinodermata</td><td>HMC</td><td>Parabolic</td><td>7.81 (0.15)</td><td>-0.23</td><td>Positive</td><td>Yes</td><td>Intercellular<sup>15</sup></td></tr><tr><td>Pencil Urchin</td><td>Eucidar is tribuloides</td><td>Echinodermata</td><td>HMC</td><td>Threshold negative</td><td>8.13 (0.02)</td><td>0.09</td><td>Non-linear positive</td><td>Yes</td><td>Intercellular<sup>15</sup></td></tr><tr><td>Red Coralline Algae</td><td>Neogoniolithon spectabile</td><td>Rhodophyta</td><td>HMC</td><td>Parabolic</td><td>9.22 (0.16)</td><td>1.03</td><td>Neutral</td><td>Yes<sup>10,11</sup></td><td>Intercellular<sup>16,17</sup></td></tr><tr><td>Green Halimeda Algae</td><td>Halimeda incrassata</td><td>Chlorophyta</td><td>A</td><td>Parabolic</td><td>—</td><td>—</td><td>—</td><td>—</td><td>Intercellular<sup>16,18</sup></td></tr><tr><td>Serpulid worm</td><td>Hydroides crucigera</td><td>Annelida</td><td>~50A:50HMC</td><td>Negative</td><td>8.17 (0.02)</td><td>0.06</td><td>Non-linear positive</td><td>Yes<sup>12,13,14</sup></td><td>Intercellular<sup>13,14</sup></td></tr><tr><td>Temperate Coral</td><td>Oculina arbuscula</td><td>Cnidaria</td><td>A</td><td>Threshold negative</td><td>8.52 (0.01)</td><td>0.41</td><td>Positive</td><td>Yes (debated)</td><td>Extracellular<sup>19,20</sup></td></tr></tbody></table>

tropical pencil urchin (*Eucidar is tribuloides*), serpulid worm (*Hydroides crucigera*), temperate coral (*Oculina arbuscula*), American oyster (*Crassostrea virginica*), bay scallop (*Argopecten irradians*), blue mussel (*M. edulis*), hard clam (*Mercenaria mercenaria*), and soft clam (*Mya arenaria*; see **Table 1** for details). Calcium carbonate powders were collected by scraping the growing edge of shells or skeletons (e.g., molluscs, corals, algae, and serpulid worms), by sampling the tip of individual spines (e.g., urchins), or by pulverizing and homogenizing the entire carapace if it was formed entirely under the experimental conditions (e.g., crabs, lobster, and shrimp). Shell or skeleton formed exclusively under the experimental conditions was identified relative to a ¹³⁷Ba spike that was emplaced in the shells/skeletons at the start of the experiment (Ries, 2011).

## Measurement of Element/Calcium Ratios

Element/calcium ratios were determined for all taxa (**Supplementary Table 3**) via inductively coupled plasma atomic emission spectrophotometry (ICP-AES). Samples were first cleaned to remove organic matter using an oxidative procedure described in Barker et al. (2003) for the analysis of foraminiferal samples. The sample set contains a number of biomaterials from taxa that have seldom or never been analyzed for X/Ca ratios, including the crustacean carapace, which is a complex biomineral containing chitin and with variable carbonate content. Given the compositional complexity and diversity of the sample set, it was not practical to optimize the cleaning and sample preparation protocol for each species, including assessing potential matrix effects for every species, as

<footnote># High Mg-calcite (HMC), low Mg-calcite (LMC), aragonite (A); mineralogy is considered mixed if the proportions exceed a 90:10 split; otherwise, if less, the biomineral is considered to be called the dominant mineral [e.g., 97A:3LMC is considered to be A in the following analyses as the detection limit for the XRD used in Ries (2011) as 3%]. More specific quantification regarding relative proportions of each mineralogy can be found in Ries, 2011;¹Weiss et al., 2002;²Addadi et al., 2006;³Immenhauser et al., 2016;⁴Harper et al., 1997;⁵Gazeau et al., 2013;⁶Meseck et al., 2018;⁷Brecevic and Nielsen, 1989;⁸Nassif et al., 2005;⁹Heinemann et al., 2011;¹⁰Leptophytum foecundum; Rahman et al., 2019;¹¹Corallina elongate; Raz et al., 2000;¹²De Yoreo et al., 2015;¹³Chan et al., 2015a;¹⁴Chan et al., 2015b;¹⁵Witt, 2002;¹⁶Halimeda sp.; Cuif et al., 2013;¹⁷Borowitzka and Larkum, 1987;¹⁸Halimeda; Wilbur et al., 1969;¹⁹Constantz, 1986;²⁰Constantz and Meike, 1989;²¹Ries, 2011;²²Ries et al., 2009;²³Liu et al., 2020;²⁴Dillaman et al., 2005;²⁵Mergelsberg et al., 2019;²⁶Luquet, 2012;²⁷Roer and Dillaman, 1984.</footnote>
this would have been prohibitively costly in terms of time and financial resources. Instead, we used the established protocol of Barker et al. (2003) that is widely used for trace element determinations within biogenic carbonates (e.g., Yu et al., 2005).

Calcium concentrations and X/Ca ratios were determined using a Varian Vista ICP-AES at the University of Cambridge, following the methods of de Villiers et al. (2002) and Barker et al. (2003). Pursuant to the procedures outlined in these studies, samples were diluted to be matrix-matched with standards at a Ca²⁺ concentration of 60 ppm and analyzed in duplicate for consistency. The relative standard deviation (%RSD) of repeated analyses of standards and samples was better than 0.2% for Mg/Ca, consistent with Barker et al. (2003).

All samples were then analyzed on a PerkinElmer Elan DRC II quadrupole inductively coupled plasma mass spectrometer at the University of Cambridge to determine Li/Ca, B/Ca, Zn/Ca, Sr/Ca, Cd/Ca, Ba/Ca, U/Ca, and Na/Ca ratios, with the exception of the crustacea (crabs, lobsters, and shrimp) and Halimeda algae due to issues with nebulizer blockage or because elemental ratios were outside of the working range established by the standards. An intensity-ratio calibration was used for determinations, following the procedure of Yu et al. (2005). The relative standard deviation (%RSD) of repeat analyses of standards and samples is better than 1.2% for all elemental ratios, consistent with Yu et al. (2005).

## Apparent Partition Coefficient and Inorganic Partition Coefficient Selection

Apparent partition coefficients (Table 2) of the carbonate biominerals were calculated by dividing the measured element-to-calcium ratios by established corresponding element-to-calcium ratios of seawater (Bruland, 1983), such that:

$$K_X = \frac{(X/Ca)_{\text{mineral}}}{(X/Ca)_{\text{seawater}}} \quad (1)$$

Inorganic elemental partition coefficients for comparison to biominerals were selected from inorganic calcite or aragonite precipitation studies (See Table 2 and Supplementary Table 6) that employed experimental conditions similar to the culturing experiment of Ries et al. (2009). The inorganic precipitation experiments included in this comparison generally used filtered or artificial seawater and were conducted at or near 25°C and within a pH₃₄ range of 7.2–8.5. A summary of the selected partition coefficients for synthetic aragonite, calcite, and amorphous calcium carbonate are presented in Table 2 (see more detailed descriptions of experimental conditions in Supplementary Table 6 and a more detailed review of the literature on element partitioning between inorganic CaCO₃ and fluid in Supplementary Material section “Description of Inorganic Partition Coefficient Selection”).

## Statistics

All statistical analyses were conducted in R version 3.6.2 (R Core Team, 2019). Outlier testing was conducted using Rosner's Test as implemented in the EnviroStat R package (v0.4-2, Le et al., 2015). This test was used because it was designed to avoid masking effects, which are losses in test power due to larger

numbers of outliers being present than being tested for. Outliers were not considered in the analyses and are not illustrated in plots. A complete spreadsheet of the raw data is provided as Supplementary Data Sheet 2.

## Models and Akaike Information Criterion (AIC)

Models relating X/Ca ratios to seawater carbonate parameters (pH₃₄, [CO₃²⁻]₃₄), net calcification rates (Ries et al., 2009), and boron isotope derived pH₄₅ (Liu et al., 2020) were conducted by fitting both linear and quadratic models to the data. Both the linear and quadratic models only include sums of terms; thus, both types of models can be considered linear regressions, evaluated, and compared. The regressions were conducted using the glm function in base R, which returned R² values and F-test statistics. Since both linear and quadratic models were considered linear regressions, R² values associated with the regression analyses were comparable and used to model the relationship between the tested variables. The F-test within the glm function is a test of overall significance that compares the specified model to a model with no predictors and subsequently generates a p-value. Regression models that returned a p-value < 0.05 were considered to indicate a statistically significant relationship between the variables of interest. Scatterplots were generated using the ggplot2, ggtext, and Rmisc packages in R and are displayed in Supplementary Figures 10–45 (v3.32, Wickham, 2016; v0.1.1, Wilke, 2020; v1.5, Hope, 2013). Modeled relationships with a p-value < 0.05 were outlined in red (Supplementary Figures 10–45). The code used to generate the plots is available on Github with links provided in Supplementary Material section “Phylogenetic Analysis.”

Akaike information criterion (AIC), which ranks models based on their fit and parsimony (Anderson and Burnham, 2004), was used to select between linear or quadratic models for the data. The AIC model that returned the smallest AIC score was selected as the optimal model for describing the relationship between the X/Ca ratio and carbonate chemistry or other measured parameter. The results of the AIC tests are summarized in Tables 3–5, which show the adjusted R² and p-values for all species exhibiting statistically significant relationships.

## Generalized Additive Models (GAMs)

We modeled the relationships between eight X/Ca ratios and Mg/Ca ratios using generalized additive models (GAMs) from the mgcv package (v1.8–31, Marra and Wood, 2011; Wood, 2017). GAMs are an extension of the generalized linear model family, with the benefit that they have no assumption of a linear relationship between the predictor and response variables, and no assumption of a specific underlying distribution in the data (Hastie and Tibshirani, 1990). Smoothing parameters (penalized regression splines) are introduced and automatically calculated via generalized cross validation (GCV) in the GAM implementation in mgcv (Wood, 2011). GAMs are frequently employed in ecological and environmental contexts where the nature of the relationship between predictors and responses is not known a priori, and/or a flexible nonparametric model is required.
<table><caption>**TABLE 2** | Summary of partition coefficients, including calculated apparent partition coefficients of the species cultured in this study; imputed values; synthetic calcite, aragonite, and amorphous calcium carbonate (ACC); and a number of published biogenic apparent partition coefficients.</caption><thead><tr><th>Organism</th><th></th><th>K<sub>Li</sub> (SD)</th><th>K<sub>B</sub> (SD)</th><th>K<sub>Na</sub> (SD)</th><th>K<sub>Mg</sub> (SD)</th><th>K<sub>Zn</sub> (SD)</th><th>K<sub>Sr</sub> (SD)</th><th>K<sub>Cd</sub> (SD)</th><th>K<sub>Ba</sub> (SD)</th><th>K<sub>U</sub> (SD)</th></tr></thead><tbody><tr><td>Temperate Coral</td><td>This study</td><td>0.0028 (0.0004)</td><td>0.014 (0.0009)</td><td>0.506*</td><td>0.0008 (0.00007)</td><td>7.853 (4.144)</td><td>0.0011 (0.00004)</td><td>0.647 (1.044)</td><td>1.439 (0.195)</td><td>0.676 (0.0948)</td></tr><tr><td>Halmeda Algae</td><td>This study</td><td>0.0081*</td><td>0.0024*</td><td>0.259 (0.029)</td><td>0.0009 (0.0007)</td><td>0.252 (0.192)</td><td>0.001 (0.00008)</td><td>5.736*</td><td>2.525 (0.697)</td><td>0.0255*</td></tr><tr><td>Soft Clam</td><td>This study</td><td>0.00215 (0.00023)</td><td>0.00064 (0.00011)</td><td>0.490*</td><td>0.00015 (0.00003)</td><td>3.966 (2.022)</td><td>0.00029 (0.00003)</td><td>0.906 (1.220)</td><td>0.7146 (0.4464)</td><td>0.0311 (0.0153)</td></tr><tr><td>Conch</td><td>This study</td><td>0.002 (0.0003)</td><td>0.0003 (0.00005)</td><td>0.567 (0.038)</td><td>0.0002 (0.00006)</td><td>2.84 (1.67)</td><td>0.0002 (0.00003)</td><td>0.187 (0.113)</td><td>0.418 (0.242)</td><td>0.014 (0.007)</td></tr><tr><td>Hard Clam (Quahog)</td><td>This study</td><td>0.00185 (0.00131)</td><td>0.00066 (0.00050)</td><td>0.482*</td><td>0.0031 (0.0044)</td><td>185.46 (184.71)</td><td>0.00026 (0.00008)</td><td>0.10166 (0.92659)</td><td>0.81877 (0.48562)</td><td>0.12254 (0.13280)</td></tr><tr><td>Limpet</td><td>This study</td><td>0.0019 (0.00013)</td><td>0.0004 (0.00009)</td><td>0.492 (0.050)</td><td>0.0002 (0.00008)</td><td>4.28 (4.39)</td><td>0.00018 (0.000021)</td><td>0.182 (0.133)</td><td>0.652 (0.793)</td><td>0.045 (0.015)</td></tr><tr><td>Whelk</td><td>This study</td><td>0.002 (0.0003)</td><td>0.0004 (0.00007)</td><td>0.508 (0.033)</td><td>0.0002 (0.00009)</td><td>12.31 (20.43)</td><td>0.0002 (0.00003)</td><td>17.84 (48.82)</td><td>0.492 (0.555)</td><td>0.049 (0.022)</td></tr><tr><td>Serpulid Worm</td><td>This study</td><td>0.0195 (0.0039)</td><td>0.0117 (0.0014)</td><td>0.433*</td><td>0.0316 (0.0044)</td><td>11.75 (6.06)</td><td>0.0007 (0.00009)</td><td>9.775 (8.302)</td><td>8.177 (6.382)</td><td>0.388 (0.098)</td></tr><tr><td>Blue Mussel</td><td>This study</td><td>0.0028 (0.0005)</td><td>0.00102 (0.00025)</td><td>0.481*</td><td>0.0052 (0.0058)</td><td>13.323 (5.825)</td><td>0.00025 (0.00017)</td><td>0.5478 (0.3158)</td><td>3.1824 (6.9398)</td><td>0.0048 (0.0010)</td></tr><tr><td>Periwinkle</td><td>This study</td><td>0.006 (0.002)</td><td>0.00007 (0.00002)</td><td>0.366 (0.065)</td><td>0.0008 (0.0003)</td><td>12.63 (10.09)</td><td>0.0001 (0.00001)</td><td>0.064 (0.061)</td><td>0.144 (0.051)</td><td>0.007 (0.004)</td></tr><tr><td>American Oyster</td><td>This study</td><td>0.0095 (0.0017)</td><td>0.0025 (0.00025)</td><td>0.459*</td><td>0.0021 (0.0003)</td><td>52.67 (29.13)</td><td>0.00012 (0.000013)</td><td>0.365 (0.121)</td><td>0.0958 (0.0144)</td><td>0.00697 (0.0194)</td></tr><tr><td>Bay Scallop</td><td>This study</td><td>0.0047 (0.001)</td><td>0.0008 (0.0002)</td><td>0.4787*</td><td>0.0038 (0.0008)</td><td>18.098 (11.507)</td><td>0.00016 (0.000008)</td><td>0.0396 (0.0124)</td><td>0.5489 (0.1901)</td><td>0.00287 (0.00156)</td></tr><tr><td>American Lobster</td><td>This study</td><td>0.0122*</td><td>0.0058*</td><td>0.688 (0.353)</td><td>0.015 (0.004)</td><td>1378.00 (698.64)</td><td>0.0006 (0.0001)</td><td>6.65*</td><td>11.69 (4.64)</td><td>0.0184*</td></tr><tr><td>Blue Crab</td><td>This study</td><td>0.0132*</td><td>0.0052*</td><td>0.438 (0.083)</td><td>0.011 (0.001)</td><td>398.35 (162.77)</td><td>0.0006 (0.00004)</td><td>6.88*</td><td>4.17 (2.62)</td><td>0.4380*</td></tr><tr><td>Gulf Shrimp</td><td>This study</td><td>0.0125*</td><td>0.0049*</td><td>0.381 (0.101)</td><td>0.009 (0.004)</td><td>2009.32 (1101.27)</td><td>0.0006 (0.0002)</td><td>6.918*</td><td>6.66 (4.69)</td><td>0.0258*</td></tr><tr><td>Coralline Red Alga</td><td>This study</td><td>0.018 (0.001)</td><td>0.012 (0.002)</td><td>0.337 (0.027)</td><td>0.052 (0.002)</td><td>142.16 (64.98)</td><td>0.0004 (0.00002)</td><td>22.87 (8.86)</td><td>2.31 (0.52)</td><td>0.108 (0.012)</td></tr><tr><td>Pencil Urchin</td><td>This study</td><td>0.02 (0.0005)</td><td>0.008 (0.001)</td><td>0.450 (0.031)</td><td>0.013 (0.0014)</td><td>9.945 (5.854)</td><td>0.0002 (0.000008)</td><td>7.118 (1.893)</td><td>0.9967 (0.2938)</td><td>0.0102 (0.0043)</td></tr><tr><td>Purple Urchin</td><td>This study</td><td>0.023 (0.002)</td><td>0.010 (0.0004)</td><td>0.517 (0.030)</td><td>0.015 (0.001)</td><td>3.99 (1.99)</td><td>0.0003 (0.00001)</td><td>4.0000 (0.4585)</td><td>0.6995 (0.2267)</td><td>0.005 (0.0008)</td></tr><tr><td>Published inorganic aragonite</td><td>See Supplementary Table 5</td><td>0.0321</td><td>0.02–2.48</td><td>–</td><td>4.72E-04</td><td>5.7</td><td>1.13</td><td>–</td><td>0.22–2.11</td><td>0.248</td></tr><tr><td>Published inorganic calcite</td><td>See Supplementary Table 5</td><td>1.4E-4–4.9E-3</td><td>1.4E-6–4.14E-3</td><td>–</td><td>2.66 × 10<sup>-2</sup></td><td>54</td><td>2.59 × 10<sup>-1</sup></td><td>18.5</td><td>4.0E-3–96.3</td><td>2.0–6.0 × 10<sup>-2</sup></td></tr><tr><td>Published ACC</td><td>See Supplementary Table 5</td><td>6.5–14.5E-3</td><td>0.01–2.34</td><td>0.6–4.06E-3</td><td>0.014–0.221</td><td>–</td><td>0.18–1.02</td><td>–</td><td>0.35–5.97</td><td>0.01–0.94</td></tr><tr><td>Coral</td><td>Published</td><td>1.22E-3–2.12E-3<sup>1</sup></td><td>2.63E-05<sup>2</sup></td><td>–</td><td>3.41E-4<sup>1,2</sup>–2.37<sup>3</sup></td><td>–</td><td>1.035<sup>2</sup>–2051<sup>3</sup></td><td>–</td><td>3770<sup>3</sup></td><td>0.0647<sup>2</sup></td></tr><tr><td>Mollusc</td><td>Published</td><td>0.003–0.01<sup>4</sup></td><td>–</td><td>–</td><td>0.0002–0.002<sup>4</sup></td><td>–</td><td>0.116–0.23<sup>4,6</sup></td><td>–</td><td>0.1–1.29<sup>5,6</sup></td><td>0.4125<sup>7</sup></td></tr><tr><td>Echinoderms</td><td>Published</td><td>0.0111–0.0958<sup>9</sup></td><td>0.0203–0.0376<sup>9</sup></td><td>–</td><td>0.0218–0.0338<sup>8,9</sup></td><td>–</td><td>0.222–0.269<sup>8</sup></td><td>–</td><td>–</td><td>–</td></tr><tr><td>Arthropods</td><td>Published</td><td>–</td><td>–</td><td>–</td><td>0.00358–0.022<sup>10,11</sup></td><td>–</td><td>36.65–58.61<sup>10</sup></td><td>–</td><td>–</td><td>–</td></tr><tr><td>Algae</td><td>Published</td><td>–</td><td>0.013<sup>13</sup></td><td>–</td><td>0.03–0.055<sup>12</sup></td><td>–</td><td>0.34<sup>12</sup></td><td>–</td><td>–</td><td>–</td></tr><tr><td>Brachiopods</td><td>Published</td><td>0.0091–0.0188<sup>14</sup></td><td>0.0297–1.2351<sup>14</sup></td><td>0.0001–0.0003<sup>14</sup></td><td>0.0008–0.0019<sup>14</sup></td><td>–</td><td>0.1332–0.2807<sup>14</sup></td><td>–</td><td>0–0.4239<sup>14</sup></td><td>–</td></tr><tr><td>Foraminifera</td><td>Published</td><td>5.85E-03<sup>19</sup></td><td>1.19E-03<sup>15,19</sup></td><td>1.35E-04<sup>19</sup></td><td>0.00031–0.00278<sup>18,19</sup></td><td>–</td><td>0.159–0.33<sup>18,19</sup></td><td>0.8–4<sup>15,17,19</sup></td><td>0.20–0.7<sup>17,19</sup></td><td>31<sup>16,19</sup></td></tr></tbody></table>

<footnote n="1">¹Montagna et al., 2014; ²Sinclair and Risk, 2006 (model); ³Gaetani et al., 2011 (model); ⁴Calculated from Dellinger et al., 2018; ⁵Gillikin et al., 2006; ⁶Zacherl et al., 2003; ⁷Frieder et al., 2014; ⁸Lavigne et al., 2013; ⁹Nguyen, 2013; ¹⁰Wickins, 1984; ¹¹Kondo et al., 2005; ¹²Ries, 2006; ¹³Donald et al., 2017; ¹⁴Junikova et al., 2020; ¹⁵Dang et al., 2019; ¹⁶Keul et al., 2013; ¹⁷Havach et al., 2001; ¹⁸Raitzsch et al., 2010; ¹⁹Allen et al., 2016. *Refers to the apparent partition coefficients that were calculated using imputed values.</footnote>
We tested models with and without the presence of the variable “Phylum” as both an additive and interactive effect to better explore the influence of taxonomy on mineralogy. It should be noted that these analyses are distinct from the phylogenetic analyses described below in section “Phylogenetic Analyses,” as they are not informed by the phylogenetic relationships between species. Models were fit using maximum likelihood, and a likelihood ratio test from package *lmtest* (v0.9–38, Zeileis and Hothorn, 2002) was used to select from among candidate models. In cases where the likelihood ratio test indicated that two models were equally likely, we selected the less complex model. Each final model fit was bootstrapped using 1,000 iterations in package *boot* (v1.3–25, Davison and Hinkley, 1997; Canty and Ripley, 2020) to produce a more robust estimate of the R² and is accompanied by a 95% confidence interval around the bootstrapped R² (**Supplementary Tables 4, 5; see Table 6** for a summary of the best-performing models).

## Imputation, Hierarchical Clustering, and Ordination

A K-Nearest-Neighbors (KNN) machine learning algorithm was employed to estimate missing values (with a pseudocount of +1 added to place the values on a log-10 scale) in the X/Ca ratio dataset (14.21% of the dataset) using the *impute* R package (v1.64.0, Beretta and Santaniello, 2016; Hastie et al., 2020). For this model, KNN selects an *a priori* neighborhood (3 in this study to get a range of values and not overfit value estimation) and imputes them using a weighted mean of the neighbors where those with the smallest Euclidean distance have more impact on the final estimation. Each “object” represents a species, and its location in the space is determined by each of the log-10 transformed X/Ca ratios. The distance between objects is their straight-line distance in the space. The analysis begins by treating each organism as a separate cluster. Then, the algorithm iteratively merges the two closest (i.e., most similar) clusters until all are merged. The main output of the hierarchical clustering analysis is the dendrogram on the left of **Figure 4**, which shows the results of the iterative clustering and reveals relationships in X/Ca ratios separated by mineralogy (**Figure 4**). The resulting dataset was used for ordination and cluster analysis.

Clustering and ordination of elemental ratios by taxon and annotation of phylum was conducted using the *vegan* and *ggplot2* packages (v3.32, Oksanen et al., 2008; v2.5–7, Wickham, 2009). Bray-Curtis dissimilarity (a statistic used to quantify the compositional dissimilarity between two different sites) was calculated for each pairwise imputed elemental ratio by species and subject to complete linkage hierarchical clustering, where the maximum distance between two clusters is identified and the new clusters are merged. This step is repeated until the dataset has agglomerated into a complete dendrogram. The same values were used to model species dissimilarity in ordination space using Non-metric Multidimensional Scaling (NMDS), which was selected because it is rank-based and detects non-linear relationships well. The model collapses information from multiple dimensions and projects this into a lower state (2 dimensions) and iteratively (999 iterations) places objects in ordination space while minimizing distance in projected dimensions as compared to the original dissimilarity matrix.

## Phylogenetic Analyses

Initial phylogenetic trees reflecting the taxonomy for the invertebrate species studied here was constructed with the phyloT application<fnref n="1" /> based on the NCBI taxonomic backbone (**Figures 4B–D**). NCBI taxonomy IDs used for the studied species were 6706 (*H. americanus*), 6763 (*C. sapidus*), 161926 (*P. plebejus*), 6550 (*M. edulis*), 6565 (*C. virginica*), 6596 (*M. mercenaria*), 6604 (*M. arenaria*), 31199 (*A. irradians*), 387448 (*S. alatus*), 31216 (*L. littorea*), 399971 (*U. cinerea*), 176853 (*C. fornicata*), 7641 (*A. punctulata*), 7632 (*E. tribuloides*), 2071691 (*Neogoniolithon* sp.), 170419 (*H. incrassata*), 1196087 (*H. crucigera*), and 1282862 (*O. arbuscula*). This tree constructed using the NCBI taxonomy was used to initially examine the correlation of species’ geochemistry and mineralogy against taxonomy.

Next, two mitochondrial (COI, cytb) gene regions and a single nuclear (18S) gene region were used to infer the phylogenetic relationships and times of divergence for all the 18 focal species. **Supplementary Material** section “Methods and Results of Phylogenetic Trees” includes details on tree assembly in BEAST 1.10 (Suchard et al., 2018). The resulting phylogeny is provided in **Supplementary Material** section “Link to Nexus File for Time-Calibrated Phylogenetic Tree.” The newly assembled time-calibrated phylogeny was first used to test for phylogenetic signal in geochemistry. Phylogenetic signal, defined as the tendency for closely related species to resemble each other more than species selected at random from the tree, is here used as a test for phylogenetic conservatism or inertia in geochemistry (Pagel, 1999). However, the coupling between evolutionary history and patterns of X/Ca divergence between species can also be confounded by alternative variables such as mineralogy (**Table 7**). Thus, regression models incorporating the newly inferred phylogeny were used to estimate the relative contributions of mineralogy and evolutionary history to differences in geochemistry among species (**Figure 5**).

Phylogenetic signals in the elemental geochemistry were estimated using both uni- and multivariate approaches. The univariate phylogenetic signal (i.e., for each elemental ratio) was computed using the lambda statistics in the phylosig function in the *phytools* R package (v0.7–70, Revell, 2012). Next, multivariate phylogenetic signals (i.e., across all elemental ratios) were estimated with Pagel’s lambda (Pagel, 1999) using the *phyl.pca* in *phytools* (Revell, 2009, 2012). For both approaches, lambda estimates range from 0 to 1, with values closer to 1 indicating that differences in elemental ratio(s) among species reflect evolutionary history (phylogenetic conservatism). Finally, the initial dataset was partitioned into five different datasets with variable number of imputed values to determine whether they affect the resulting phylogenetic signal (see details in **Supplementary Table 9**). First, the phylogenetic signal of an elemental dataset, including all of the elements except Cd/Ca and Na/Ca Second, the original dataset was subsampled for only Zn/Ca, Mg/Ca, Sr/Ca, and Li/Ca data. This dataset did not include any imputed values. Third, since almost half of elemental data for the arthropods were imputed, the shrimp,

<footnote n="1">¹https://phylot.biobyte.de</footnote>
crabs, and lobsters were excluded from the analyses. Fourth, a dataset with all species was generated that included only Mg/Ca and Sr/Ca data. Fifth, a dataset of all species was generated that included only Mg/Ca, Sr/Ca, and Li/Ca. This set of five additional datasets was used to examine the effects of data imputation on the multivariate phylogenetic signal estimate. Lastly, the relative contributions of the phylogeny and mineralogy to geochemistry patterns among species was estimated using phylogenetic and non-phylogenetic regression models. It should be noted that patterns of divergence between X/Ca ratios and phylogenetic signal can still be explained by differences in polymorph mineralogy. Thus, a set of three regression models were fit for each of the elemental ratios. Each set of regression models included an intercept-only linear regression model (model 1), a linear regression model between the target elemental ratios and mineralogy (model 2), and a phylogenetic regression model between the target elemental ratios and mineralogy (model 3). Non-phylogenetic models were fit using the lm function in the stats R package. Phylogenetic regression models, with lambda values estimated from the data, were fit using the phylolm function in the *phylolm* R package (v2.6.2; Ho et al., 2016). The relative contribution of differences in biogenic carbonate mineralogies among species to elemental ratios among species was estimated by comparing partial *R*²-values between non-phylogenetic models (models 1 and 2). Similarly, the relative contribution of evolutionary history to explaining differences in elemental ratios among species was estimated by comparing partial *R*²-values between the phylogenetic regression and the linear regression models of elemental ratios versus mineralogy across species and *pCO₂* conditions (models 2 and 3). Comparison of partial *R*²-values, following Ives (2018), was conducted using the *rr2* R package (v1.0.2, Ives and Li, 2018).

## RESULTS

The reported X/Ca ratios for each species represent means of the measured data across the *pCO₂* culture conditions (Supplementary Table 3). These mean values were used for the multivariate and phylogenetic analyses described below because the range of X/Ca values for a given species across *pCO₂* treatments is relatively small compared to the differences in mean X/Ca of different species. For example, the pencil urchin, *E. tribuloides*, exhibits Li/Ca and Mg/Ca values of 50–53 µmol/mol and 60–78 mmol/mol across *pCO₂* treatments, respectively, whereas the oyster, *C. virginica*, exhibits values of 15–30 µmol/mol and 9.5–13 mmol/mol, respectively (Figure 1).

## Element/Calcium Ratios Amongst Diverse Marine Calcifiers

As expected, Mg, Sr, and Na were the most abundant elemental substitutions across all mineralogies (Figure 1 and Supplementary Table 3; Reeder, 1983; Sugawara and Kato, 2000; Allison et al., 2001; Iglikowska et al., 2016). The calcite in crustaceans’ carapaces and urchins’ spicules exhibits Mg/Ca in the range of 50–100 mmol/mol and the mixed mineralogy serpulid

worm and calcitic coralline red algae exhibit relatively high Mg/Ca in the range of 100–300 mmol/mol, while the aragonite-producing species were low in Mg/Ca. Molluscs exhibited low Mg/Ca, regardless of the *CaCO₃* polymorph produced. Sr/Ca was highest in the range of 9–12 mmol/mol in the aragonite-forming coral and green algae. Strikingly, mollusc *CaCO₃* of all polymorphs was relatively low in all elemental ratios compared to other species, except for Li/Ca compared to the temperate coral and Na/Ca compared to most other species. Of the other elements, Li/Ca and B/Ca was notably higher in organisms that formed high-Mg calcite, including the serpulid worm, which produces a worm tube of aragonite and high-Mg calcite (see more analyses below in section “Relationships Between Mg/Ca and Other Elemental Ratios”). B/Ca was also high in the aragonite coral skeleton. The coralline red algae and serpulid worm generally had comparatively high levels of all measured elements compared to the other species. Amongst the aragonite producing species, the coral was noticeably different in that it exhibited much higher B/Ca, Sr/Ca, and U/Ca. The *Halimeda* algae exhibited elemental ratios similar to those of the coral.

## Relationships Between Mg/Ca and Other Elemental Ratios

Generally, higher incorporation of Mg in the biomineral lattice appears to be linked to higher incorporation of other elements (Figure 2). For example, for Li/Ca, the molluscs and coral group together with lower Mg content (Li/Ca: 6.5–25 µmol/mol; Mg/Ca: 0.5–6 mmol/mol), while the echinoderms, serpulid worms, and coralline red algae are relatively enriched in both Mg and Li (Li/Ca: 39–60 µmol/mol; Mg/Ca: 60–275 mmol/mol). Potential increases in Zn/Ca with Mg incorporation are less clearly resolved.

In this X/Ca-Mg/Ca space, the data appear to group by phylum, so generalized additive models (GAMs) were used to explore the influences of phylum and mineralogy on elemental clustering. The final GAMs explain between 26.4 and 99% of the variation in the elemental ratio of interest, with model selection indicating a preferred final model with phylum, or a combination of phylum and mineralogy, for seven of the eight elements (see Table 6 for final models; see Supplementary Table 4 for complete list of candidate models; see in Supplementary Table 5 for likelihood ratio test results; see Supplementary Figure 9 for final model components). Models that only consider Mg/Ca explain up to 71.9% of the variation in the element of interest but were in no case preferred by model selection. In seven of eight cases, including mineralogy with Mg/Ca substantially improves the models and explains between 16.6 and 95.6% of the variation in the element of interest. Cd/Ca is the only case where the model including only mineralogy was preferred. Six of the eight models were improved by including phylum, explaining up to 98.7% of the variation. All final models produced normally distributed residuals (Supplementary Figures 1–8). Overall, we find evidence that both mineralogy and phylum are key predictors of elemental ratios in the study organisms. Please note that these models consider the phylum that were used to classify the samples (i.e., Rhodophyta, Chlorophyta,
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Elemental ratio data across all pCO2 culture conditions combined for the invertebrates studied. Color and shape indicate phylum and mineralogy, respectively. The species are organized from left to right according to mineralogy: Aragonite, Mixed Mineralogy, and Calcite. Mixed Mineralogy refers to carbonates with less than 90% of either aragonite or calcite. Differences between species and phyla are typically much greater than the variation within species data that could be attributed to varying seawater conditions. Note some elemental ratios are not displayed to avoid scale compression or due to low confidence in their analytical accuracy (e.g., arthropod Zn/Ca; see section "Materials and Methods").</figcaption></figure>

Cnidaria, etc.); the models are not informed by the phylogenetic relationships between species. The results of the phylogenetically informed analyses are described below in section "Phylogenetic Signal and Regression Analyses."

## Comparison of Inorganic Carbonate Element Partitioning to Marine Biominerals

Patterns of elemental partition coefficients ($K_x$) observed here are of course similar to those observed for the X/Ca data. However, calculating the apparent partition coefficients allows for comparison to other cultured biogenic samples and published

partition coefficients derived from inorganic synthesis of calcite, aragonite, and amorphous calcium carbonate (ACC; **Table 2** and **Figure 3**). Though some apparent partition coefficients align with those of inorganically precipitated CaCO$_3$ minerals (e.g., $K_{Li}$ for aragonitic organisms and ACC; $K_{Li}$ for calcitic organisms and calcite), there are many that are relatively enriched or depleted with respect to those observed for inorganic precipitates. Positioning relative to the inorganic partition coefficients appears to be consistent amongst organisms of the same phylum regardless of mineralogy, as observed, for example, in the mollusca.

The molluscs are relatively depleted with respect to most of the measured elements compared the other organisms,
<figure><img src="image_4.png" /></figure>

**FIGURE 2** | Elemental ratio data plotted against Mg/Ca (log10-transformed values are shown. Color and shape indicate phylum and mineralogy, respectively. Generalized additive models (GAMs) of this data show a number of significant relationships between elemental ratios and Mg/Ca, explaining 0–71% of the variance in the data (**Supplementary Table 4**). Grouping the organisms by phylum or mineralogy substantially improves model performance, and the best-performing models consider either or both grouping variables (38.70–99.00% of variance explained; **Table 6**).

while the coralline red algae and echinoderms are relatively enriched with respect to these elements. For $K_{Li}$, the coral and the molluscs, regardless of mineralogy, exhibited values aligned with the published inorganic partition coefficients of ACC. The apparent $K_B$ values of the calcitic organisms all lie closest to the published partition coefficient for calcite (Figure 3 and Table 2). The apparent $K_B$ of aragonitic coral (0.014 ± 0.0009) is depleted relative to inorganic aragonite (0.02–2.48). The apparent $K_B$ of the serpulid worm tube (0.0117 ± 0.0014) falls between the published values for inorganic calcite (1.4E-6-4.14E-3) and inorganic aragonite (0.02–2.48), consistent with its mixed mineralogy.

$K_{Zn}$ values spanned the greatest range across the organisms, with the crustaceans exhibiting values >200. Only the coralline red algae exhibited $K_{Zn}$ that was comparable to that of the crustacea, ranging from approximately 100 to 200. With the exception of the calcitic Eastern Oyster, which has an apparent partition coefficient similar to that of inorganic calcite, the remainder of the organisms sampled align with the $K_{Zn}$ published for inorganic aragonite.

Surprisingly, apparent $K_{Sr}$ for all species studied were lower relative to published values for inorganic precipitates of corresponding carbonate polymorph (Schöne et al., 2010). The highest apparent $K_{Sr}$ for organisms were calculated
<figure><img src="image_5.png" /></figure>

**FIGURE 3** | Heatmaps displaying the log10-transformed differences between apparent partition coefficients and published partition coefficients from inorganic mineral precipitation experiments (See **Table 2** for values; **Supplementary Table 5** for experimental conditions and citations for inorganic values). The apparent elemental partition coefficients calculated for the cultured species studied typically do not align with values derived from synthetic minerals, and species belonging to the same phyla appear to exhibit similar patterns. The warm color gradient (i.e., the left, yellow-black gradient) represents when biogenic carbonates are enriched in an element relative to the published synthetic values (KSynthetic < KBiogenic). Alternatively, a green-yellow color gradient (i.e., the right gradient) represents when biogenic carbonates are depleted relative to the published synthetic values (KBiogenic < KSynthetic). White pixels indicate that no published, synthetic values for that element were found. Circles within pixels denote imputed values (See section "Imputation, Hierarchical Clustering, and Ordination" and **Figure 4** for details). Apparent partition coefficients were calculated by dividing the measured and imputed elemental ratios of the biogenic carbonates by established seawater minor and trace element concentrations (Bruland, 1983). The top row of panels presents the aragonitic organisms, while the bottom panels represent the calcitic organisms. The label above each panel denotes the inorganic, synthetic partition coefficients being compared to (i.e., aragonite, calcite, and amorphous calcium carbonate; ACC). The color of the pixels was determined by log10([ KX, Biogenic - KX, Synthetic]). The absolute value is used since y cannot be less than 0 in log10(y).

from the aragonitic species of coral and *Halimeda* algae sampled, which range from 0.001 to 0.00125. The next highest values were calculated for the serpulid worm and the arthropods—approximately half of that of the aragonite species—which possess shells of mixed mineralogy and high-Mg calcite, respectively. The $K_{Sr}$ of the coralline red algae (4E-4 ± 0.2), molluscs (1.8-2.9E-4), and echinoderms (2.0-3.0E-4) are even lower.

All of the molluscs, as well as the aragonitic coral and *Halimeda* algae, have apparent $K_{Mg}$ similar to that of synthetic aragonite (4.72E-4). The high-Mg calcite arthropods and echinoderms are enriched relative to the values for

published inorganic calcite (2.66E-2), while the serpulid worm is comparable to value for inorganic calcite. Coralline red algae high-Mg calcite exhibits the highest $K_{Mg}$ at approximately 0.06.

There were no published partition coefficients for Cd in synthetic aragonite. The published $K_{Cd}$ for synthetic calcite (18.5) overlaps with the highest and lowest ends of the respective ranges of apparent $K_{Cd}$ calculated for the serpulid worm (9.78) and the coralline red algae (22.9). The coral, molluscs and echinoderms are all extremely depleted with respect to Cd relative to synthetic calcite.

The apparent $K_{Ba}$ values calculated for the crustacea and the serpulid worm are higher than the other species and published
partition coefficients for synthetic CaCO₃ minerals, though their values range widely from 3 to 20. The lowest limit of this range as well as both the *Halimeda* and the coralline red algae appear to align with the *K<sub>Ba</sub>* measured for synthetic ACC. The echinoderms, which produce high-Mg calcite biominerals, provide an apparent *K<sub>Ba</sub>* consistent with that derived from synthetic calcite experiments. The apparent *K<sub>Ba</sub>* of the mollusc species appear to be consistent, regardless of mineralogy, and lie between the published *K<sub>Ba</sub>* values for synthetic calcite and aragonite. The values calculated for coral fall between the published values for synthetic calcite and ACC.

The highest apparent *K*<sub>U</sub> are calculated for the coral and the serpulid worm, ranging from 0.5 to 0.75 and 0.25 to 0.40, respectively, and are both enriched relative to the published *K*<sub>U</sub> values for synthetic ACC, calcite, and aragonite. The molluscs, regardless of mineralogy, and the echinoderms partition U similarly to synthetic calcite, except for the mixed mineralogy hard clam, which has an average calculated *K*<sub>U</sub> similar to values published for synthetic aragonite. Coralline red algae appear to partition U similarly to synthetic ACC.

## Imputed Values and Hierarchical Clustering

Some of the analyzed specimens were missing elemental ratios (14.2% of the dataset). These missing data were imputed using the data of similar points in the space (see Methods 2.5.3 for details; see imputed results in **Supplementary Table 3**) in order to facilitate multivariate and phylogenetic analysis of the dataset. These mean and imputed values for elemental ratios are used for the hierarchical clustering (**Figure 4A**) and phylogenetic signal analyses (**Table 7** and **Supplementary Table 8**), as well as for the non-metric multidimensional scaling ordinations (NMDS; **Figures 5, 6A,B**) described below.

When compared to the taxonomic trees (**Figures 4B–D**), the hierarchical clustering results of the elemental ratio data (**Figure 4A**) appear to recapitulate the mineralogically specific taxonomic trees (**Figures 4C,D**). When including species of differing polymorph mineralogy, the similarities between the clustering analysis and the taxonomic tree (**Figure 4B**) seem to be confounded by polymorph mineralogy. Parsing out the confounding effects of mineralogy from the effects of phylogeny is explored below.

## Phylogenetic Signal and Regression Analyses

Uni- and multivariate analyses of phylogenetic signal based on the Pagel's lambda statistics, reflecting the extent to which closely related species show similar X/Ca ratios, generally indicate a strong coupling between the phylogenetic tree and patterns of elemental ratios between species. While a Pagel's lambda value of 0 for X/Ca ratios would indicate that closely related species are more likely to differ in X/Ca ratios than distantly related species (Pagel, 1999), the multivariate analyses of all nine X/Ca ratios recovered an estimated lambda of 0.89–0.99 (**Table 7**; range based on alternative datasets), indicating that closely related species tend to have similar X/Ca ratios. However, patterns of congruent

evolutionary history and X/Ca ratios are not generalizable across all X/Ca ratios. For instance, three out of the nine X/Ca ratios examined in this study had phylogenetically decoupled patterns of change in ratios between species (Li/Ca, Zn/Ca, and Ba/Ca; **Table 7**).

We also explored the potential impact of the imputed values on the phylogenetic signal results by determining Pagel's lambda with subsets of the data that minimize the number of imputed values (e.g., 0.938 for all X/Ca ratios except Cd/Ca and Na/Ca; 0.968 for all species except Arthropods) and found no substantial effects of missing data on estimate of phylogenetic signal (**Supplementary Table 9**). Furthermore, the impact of using a lower number of X/Ca was assessed through determining Pagel's lambda using subsets of two, three, and four elemental ratios (**Table 7**). Pagel's lambda remains high (0.925 for multivariate phylogenetic signal across Mg/Ca, Sr/Ca; 0.803 for Mg/Ca; 0.835 for Sr/Ca; **Table 7**). However, when an elemental ratio with no apparent phylogenetic signal is included in the analysis, lambda slightly decreases (0.898; **Table 7**). Thus, the analyses suggest that for some, but not all, of the examined X/Ca ratios, phylogenetic closeness is potentially a major driver of differences in X/Ca values amongst species.

Further, regression analysis was used to examine whether interspecific differences in X/Ca ratios are better explained by phylogeny or differences in mineralogy. As shown above, phylogenetic signal varies between X/Ca ratios (**Table 7**), so the effects of phylogenetic history relative to mineralogy are expected to vary amongst elemental ratios. Phylogenetic and non-phylogenetic models were used to resolve the relative contributions of evolutionary history and mineralogy to explaining patterns of X/Ca ratios between species. When classifying mineralogy as calcite, aragonite, or mixed, phylogeny was 1.4–59 times better at explaining variance in X/Ca ratios than mineralogy in five out of the nine X/Ca ratios (i.e., B/Ca, Mg/Ca, Sr/Ca, Cd/Ca, and Na/Ca; **Figure 7** and **Supplementary Table 9**). However, when mineralogical classification is broken down further (high-Mg calcite, low-Mg calcite, aragonite, and mixed), the amount of variance explained by phylogeny decreases to 1.3–15.7 times better in four out of the nine X/Ca ratios, demonstrating the importance of distinguishing between low-Mg and high-Mg calcites. Also notable is that in all but one of the X/Ca ratios (U/Ca), high phylogenetic signal was coupled with higher explanatory effects of phylogeny relative to mineralogy. Overall, these phylogenetic analyses suggest that elevated phylogenetic signals generally indicate that differences in X/Ca ratios are better explained by phylogenetic closeness, not mineralogy. Thus, the elevated phylogenetic signal recovered across all X/Ca ratios using multivariate analyses reflects the major contributions of phylogeny to generating extant patterns of X/Ca ratios amongst species.

Finally, mollusc elemental data were used to examine the relative contributions of mineralogy and phylogeny to explaining differences in X/Ca ratios between species within the same phylum. Regression analyses were conducted comparing the relative contributions of mineralogy and phylogenetic history using phylogenetic and non-phylogenetic models. The mollusk dataset includes four species assigned to the calcite category,
A **Hierarchical Clustering & Heatmap**

<figure><img src="image_6.png" /></figure>

B **Taxonomic Tree: All Organisms**

<figure><img src="image_7.png" /></figure>

C **Taxonomic Tree: Calcite Only**

<figure><img src="image_8.png" /></figure>

D **Taxonomic Tree: Aragonite Only**

<figure><img src="image_9.png" /></figure>

**FIGURE 4** | (A) Hierarchically clustered (Bray-Curtis dissimilarity) heatmap of log10 values of element to calcium ratios separated by species. Letters in ellipses correspond to organism biomineralogy: A, aragonite; HMC, high Mg calcite; LMC, low Mg calcite; and M, mixed. Lighter colors indicate higher values from log10-transformations of elemental ratios, and an “O” on top of a pixel indicates that the value was imputed using a K-Nearest Neighbors model estimation (see section “Materials and Methods” for details). (B–D) Taxonomic trees of the studied organisms based on genetic similarity/difference for all, calcitic, and aragonitic species, respectively. Taxonomic trees were built using the phyloT application based on the NCBI taxonomic backbone (see section “Phylogenetic Analyses”; Letunic and Bork, 2019). Qualitative similarities between the resulting dendrogram of the clustering analysis and the mineral-specific taxonomic trees suggest that taxonomic patterns arise in the elemental ratio measurements but are potentially convoluted by mineralogical effects.

four assigned to the aragonite category, and a single species to the mixed mineralogy category. It should be noted that the limited sample size may compromise the calculations assessing phylogenetic signal at the class level presented below. Regression models run for X/Ca ratios with high phylogenetic signals (B/Ca, Mg/Ca, and Sr/Ca; Table 7) found that for two of the three X/Ca ratios examined within molluscs, phylogenetic history better explains patterns of divergence in X/Ca ratios than mineralogy. Phylogeny is about two times better than mineralogy (partial R²phylogeny = 0.239 vs. mineralogy = 0.097) in predicting B/Ca. Similarly, phylogeny explains nearly two times more variance in Mg/Ca than mineralogy (partial R², phylogeny = 0.60 vs. mineralogy = 0.33). For Sr/Ca, both phylogeny and mineralogy explain a similar fraction of differences in ratios between species (partial R², phylogeny = 0.632 vs. mineralogy = 0.65). Overall, results within molluscs for X/Ca ratios with high phylogenetic signal indicate that phylogeny is generally a better predictor of differences in X/Ca ratios between species than mineralogy. Nevertheless, as indicated above, these results were potentially affected by the limited sample size and should be re-examined in future studies including more mollusc species. This analysis of

a mollusc-only dataset reaffirms that, for specific X/Ca with high phylogenetic signal, differences in X/Ca between species largely reflect phylogenetic closeness, not mineralogy.

**Non-metric Multidimensional Scaling
(NMDS) Ordination**

Non-metric Multidimensional Scaling ordinations project the all of X/Ca ratios, including the imputed values, project in two dimensions (Figures 5, 6). As opposed to Euclidean distance, which is the distance computed for the hierarchical clustering analysis, NMDS uses rank-based ordering and is better suited for detecting non-linear relationships. The model iteratively collapses the nine dimensions of X/Ca ratios onto two dimensions, seeking to minimize the distances between species.

Not surprisingly, there is a split with the dataset between the calcitic and aragonitic organisms, with the 50HMC:50A serpulid worm occupying an intermediate position between the groups (Figure 5). Within mineralogy, however, species cluster according to phylogenetic relationships (Figure 5) much more strongly than in the hierarchical clustering analysis. NMDS
<figure><img src="image_10.png" /></figure>

**FIGURE 5** | Non-metric Multidimensional Scaling (NMDS) ordination of the elemental ratio data projected to two-dimensions and colored by taxon phylum. Shape denotes mineralogy. PERMANOVA determined whether clustering by phyla is significant (*p* = 0.001). Multivariate Levene's test of homogeneity (*p* = 0.183) and the stress value for this analysis (0.086) are reported. Visualizing the elemental ratio data in this way further illustrates the relative influences of mineralogy and phylogeny and potentially provides a framework for interpreting diagenetic effects of fossil carbonates.

ordination of all organisms shows statistically significant clusters of X/Ca ratios by phylum (i.e., for species of same mineralogy), where the ratios within phylum show equal variances across organisms (PERMANOVA: *p* = 0.001; permutational multivariate Levene's test: *p* = 0.183). For example, the bivalve molluscs fall close together, as do the arthropods and echinoderms. The aragonitic *Halimeda* green algae and temperature coral aragonitic species lie proximal, as do the high-Mg calcite producing coralline red algae and the serpulid worms of mixed mineralogy.

Separating the calcite-forming (Figure 6A) and the aragonite-forming (Figure 6B) organisms removes the role of mineralogy while revealing rank-based positioning of the X/Ca ratios, themselves. In the calcite-only NMDS, significant clustering of organisms by phylum is also observed (Figure 6A; PERMANOVA: *p* = 0.003), along with variation within phyla. The X/Ca ratios position themselves in the ordination space with a similar multivariate structure as the organisms; they lie closest to organisms that incorporate them the most. Of the arthropods, Ba/Ca has the most similar data structure to the American lobster and Gulf shrimp. Li/Ca, Na/Ca, and Sr/Ca in the calcitic organisms have data structures similar to those elements in the calcitic molluscs, although there is some variation amongst them. Also, Li/Ca in addition to B/Ca in the calcite ordination to have data structures similar to those elements in the echinoderms. The coralline red alga exhibits the most distinct data structure,

exhibiting high ratios of Mg/Ca, Cd/Ca, and U/Ca. Zn/Ca in the calcite ordination does not appear to have a data structure similar to that of any of the other calcitic organisms. In the aragonite-only NMDS, significant clustering by phyla is again observed (Figure 6B; PERMANOVA: *p* = 0.136). The temperate coral and the *Halimeda* green algae have data structures more similar to B/Ca and Cd/Ca, respectively. The aragonitic molluscs, with the exception of the hard clam, exhibit relatively high Na/Ca, Li/Ca, and Ba/Ca. Also of note is the clustering of Sr/Ca, Mg/Ca, Ba/Ca, and Li/Ca for the aragonitic species, potentially indicating similar mechanisms of partitioning.

## Relationships Between Elemental Ratios, Carbonate Chemistry, and Other Measured Parameters

Relationships between the elemental data and a number of the measured experimental parameters were explored via linear and quadratic regres sional analysis, including net calcification rate (Ries et al., 2009), culture water carbonate ion concentration ([CO₃²⁻]), culture seawater pH (pHsw), and δ¹¹B-derived internal calcifying fluid pH (pHcf; Liu et al., 2020; Supplementary Figures 10–45). If both linear and quadratic fits were significant, Akaike Information Criterion (AIC) was used to compare and identify the best fit (Tables 3–5).
<figure><img src="image_11.png" /></figure>

**FIGURE 6** | NMDS ordination of elemental ratios for (A) calcitic and (B) aragonitic organisms only, which are separated by species. The elemental ratios themselves are also represented in the ordination space. The positioning of the elemental ratios depends on the apparent level of incorporation in the respective biogenic carbonates; the closer to an organism the trace element ratio, the higher the level of incorporation of that element relative to the other organisms represented in the space. Considering only one mineralogy at a time excludes the influence of mineralogical differences. The color of a point denotes phyla. (A) Similar to **Figure 4B**, PERMANOVA determines significant clustering of points by phylum (*p* = 0.003). The stress value for this analysis is 0.035. (B) Similar to **Figure 4B**, PERMANOVA determines significant clustering of points by phylum (*p* = 0.136). The stress value for this analysis is 8.84E-5.

A total of 16 species displayed at least one significant relationship between an X/Ca ratio and seawater [CO₃²⁻]. Nine species exhibited more than four significant linear or quadratic relationships, including the coral, coralline red algae, hard clam, limpet, oyster, pencil urchin, serpulid worm, soft clam, and whelk. The pencil urchin exhibited the most significant relationships between the measured elemental ratios versus the carbonate chemistry and other measured parameters. Of the elements, Sr/Ca had the most significant relationships, while only one significant relationship was identified for Na/Ca between both linear and quadratic fits for all species versus the carbonate chemistry and other measured parameters.

Net calcification rate, as determined by the buoyant weight method described in Ries et al. (2009), has 24 significant relationships with X/Ca ratios across all species. The species with the most significant relationships with changing net calcification rate was the hard clam, and the elements most frequently correlated with calcification rate across species are Sr and Ba. Sr/Ca significantly increases with increasing net calcification rate in a few molluscs, including the soft clam, blue mussel, and limpet. The aragonitic molluscs, the soft clam and the limpet, exhibit Sr/Ca that increases linearly with higher net calcification rates (R² = 0.74, *p* = 0.0002; R² = 0.342, *p* = 0.02), while the calcitic blue mussel Sr/Ca exhibits a parabolic variation with net calcification rate (R² = 0.452, *p* = 0.01). The pencil urchin Sr/Ca exhibits an inverse parabolic pattern with respect to net calcification rate (R² = 0.461, *p* = 0.01). Measured Ba/Ca significantly varies a parabolic pattern with increasing net calcification rate in the coralline red alga, lobster, and oyster (R² = 0.733, *p* = 0.02; R² = 0.99, *p* = 0.05; R² = 0.661, *p* = 0.03), while increasing linearly for the serpulid worm tube and hard clam (R² = 0.571, *p* = 0.003; R² = 0.211, *p* = 0.02).

The [CO₃²⁻] of the culture seawater exhibited 28 significant relationships with X/Ca ratios across all of the studied species (Table 4). Organisms that exhibited the most significant relationships with [CO₃²⁻] are the coral, limpet, pencil urchin, and serpulid worm. The elements that most frequently varied with [CO₃²⁻] were Sr/Ca, Zn/Ca, and Ba/Ca. Sr/Ca significantly decreases with decreasing [CO₃²⁻] in the serpulid worm and a number of molluscs, including the oyster, soft clam, conch, whelk, and limpet. The mixed mineralogy serpulid worm and the oyster both decrease parabolically with decreasing [CO₃²⁻] (R² = 0.309, *p* = 0.008; R² = 0.541, *p* = 0.004), while the aragonitic soft clam and conch decrease following inverse parabolic patterns (R² = 0.66, *p* = 0.001; R² = 0.354, *p* = 0.004). The aragonitic gastropod molluscs, the whelk and limpet, both linearly decrease with decreasing [CO₃²⁻] (R² = 0.156, *p* = 0.03; R² = 0.201, *p* = 0.03). Zn/Ca also systematically varies with decreasing seawater [CO₃²⁻] in the coralline red algae, serpulid worm, coral, blue mussel, and gulf shrimp. Yet, despite sharing the same high-Mg calcite mineralogy, Zn/Ca levels in shrimp are an order of magnitude higher than those measured in the coralline red algae. For the coral, Zn/Ca linearly increases with decreasing seawater [CO₃²⁻] (R² = 0.0982, *p* = 0.0005), while the Zn/Ca ratio in the serpulid worms and blue mussels and span comparable ranges and exhibit similar inverse parabolic patterns with decreasing [CO₃²⁻] (R² = 0.258, *p* = 0.02; R² = 0.539, *p* = 0.04). Significant changes in biomineral Ba/Ca occur with decreasing seawater [CO₃²⁻] were observed for the lobsters, echinoderms, and a suite of molluscs, including the hard clam, oyster, conch, and whelk. In the calcitic oyster, the decrease is slight and linear (R² = 0.632; *p* = 0.004), while in the aragonitic hard clam, conch, and whelk, Ba/Ca exhibits an inverse parabolic pattern with decreasing [CO₃²⁻] (R² = 0.784, *p* = 0.0004; R² = 0.257, *p* = 0.004;
<figure><img src="image_12.png" /></figure>

**FIGURE 7** | Estimated relative contributions of phylogeny and mineralogy on X/Ca ratios using phylogenetic and non-phylogenetic models on the full species-level dataset. (A) Bar graphs illustrating the phylogenetic signal, λ, (B,C) relative contributions of evolutionary history (phylogeny) and mineralogy on differences in X/Ca ratios among 18 focal species. Partial R²s for each predictor within each given element-to-calcium ratio follow Ives (2018) and Ives and Li (2018). Full results are provided in **Supplementary Table 9**. "*" denotes elemental ratios that have three or less values imputed. Mg, Sr, Li, and Zn datasets have zero imputed values. Striped bars denote elemental ratios (Cd and Na) that consider five or more imputed values, potentially influencing the resulting λ. and R². The difference between graphs (B,C) are how mineralogy was defined for the analysis: (B) Mineralogy is categorized as whether an organism produces calcite or aragonite, or a mixture of the two; (C) Mineralogy is broken down further, categorizing mineralogy as HMC, high-Mg calcite; LMC, low-Mg calcite; A, aragonite; or M, mixture (Table 1). The high phylogenetic signals and high partial-R² Phylogeny values deconvolve phylogenetic and mineralogical impacts on the elemental ratios, highlighting the importance of considering phylogeny when assessing elemental ratio data.

$R^2 = 0.193$, $p = 0.04$). A similar inverse parabolic pattern in Ba/Ca with decreasing seawater [CO₃²⁻] emerges for the high-Mg calcite lobster and pencil urchin ($R^2 = 0.763$, $p = 0.006$; $R^2 = 0.752$, $p < 0.0001$), and the purple urchin exhibits a linear increase in Ba/Ca with decreasing seawater [CO₃²⁻] ($R^2 = 0.335$, $p = 0.01$).

The pH<sub>sw</sub> of the culture tanks was measured, and the pH<sub>cf</sub> has been calculated for 10 of the species based on δ<sup>11</sup>B measurements of the same specimens evaluated in the present study (Liu et al., 2020). Twenty-five significant linear or quadratic relationships were identified between pH<sub>sw</sub> and elemental ratios across all species. Of the 10 species with pH<sub>cf</sub> data, there were 13 significant relationships between pH<sub>cf</sub> and elemental ratios. The species that exhibit the most significant relationships between elemental ratios and

pH<sub>sw</sub> are the coral (3; **Supplementary Figures 22, 23**), coralline red alga (3; **Supplementary Figures 24, 25**), limpet (3; **Supplementary Figures 32, 33**) and pencil urchin (3; **Supplementary Figures 34, 35**). The coralline red alga exhibits the most significant relationships between elemental ratios and pH<sub>cf</sub> (4; **Supplementary Figures 24, 25**). The elements that most frequently vary with pH<sub>sw</sub> are B/Ca, Sr/Ca, and Ba/Ca. Likewise, B/Ca and Sr/Ca most often correlate with pH<sub>cf</sub>. B/Ca significantly decreases with decreasing pH<sub>sw</sub> in the coralline red alga, pencil urchin, and limpet, while significantly increasing with pH<sub>cf</sub> in the coralline red algae, serpulid worm, and hard clam. For both pH<sub>sw</sub> and pH<sub>cf</sub>, the high-Mg calcite coralline red alga appears to maintain its B/Ca under moderately lower pH<sub>sw</sub> conditions before exhibiting a decrease under even more acidic conditions ($R^2 = 0.695$, $p = 0.002$; $R^2 = 0.681$, $p = 0.004$).
<table><caption>**TABLE 3** | Summary of AIC test results for significant relationships (*p* < 0.05) between element/Ca ratios and net calcification rates across all organisms.</caption><thead><tr><th rowspan="2">Organism</th><th rowspan="2">Phylum</th><th rowspan="2">X/Ca</th><th colspan="4">Best-fit regression</th></tr><tr><th>TYPE</th><th>p</th><th>R<sup>2</sup></th><th>Relationship</th></tr></thead><tbody><tr><td>American lobster</td><td>Arthropoda</td><td>Ba</td><td>Para</td><td>0.00498</td><td>0.99</td><td>+</td></tr><tr><td rowspan="2">Blue mussel</td><td rowspan="2">Mollusca</td><td>Sr</td><td>Para</td><td>0.0146</td><td>0.452</td><td>+</td></tr><tr><td>Cd</td><td>Para</td><td>0.0192</td><td>0.999</td><td>-</td></tr><tr><td rowspan="2">American oyster</td><td rowspan="2">Mollusca</td><td>Li</td><td>Para</td><td>0.00317</td><td>0.916</td><td>+</td></tr><tr><td>Ba</td><td>Para</td><td>0.0288</td><td>0.661</td><td>+</td></tr><tr><td rowspan="4">Hard clam</td><td rowspan="4">Mollusca</td><td>Mg</td><td>Linear</td><td>0.0276</td><td>0.251</td><td>+</td></tr><tr><td>Sr</td><td>Linear</td><td>0.0413</td><td>0.213</td><td>+</td></tr><tr><td>Ba</td><td>Linear</td><td>0.0026</td><td>0.571</td><td>+</td></tr><tr><td>U</td><td>Linear</td><td>0.0284</td><td>0.508</td><td>+</td></tr><tr><td rowspan="3">Soft clam</td><td rowspan="3">Mollusca</td><td>Zn</td><td>Linear</td><td>0.0446</td><td>0.595</td><td>-</td></tr><tr><td>Sr</td><td>Linear</td><td>0.000196</td><td>0.742</td><td>+</td></tr><tr><td>Cd</td><td>Linear</td><td>0.0149</td><td>0.86</td><td>-</td></tr><tr><td>Bay scallop</td><td>Mollusca</td><td>Mg</td><td>Para</td><td>0.00587</td><td>0.61</td><td>-</td></tr><tr><td rowspan="3">Limpet</td><td rowspan="3">Mollusca</td><td>Li</td><td>Para</td><td>0.0127</td><td>0.831</td><td>+</td></tr><tr><td>Mg</td><td>Para</td><td>0.00632</td><td>0.471</td><td>+</td></tr><tr><td>Sr</td><td>Linear</td><td>0.0196</td><td>0.325</td><td>+</td></tr><tr><td rowspan="3">Pencil urchin</td><td rowspan="3">Echinodermata</td><td>B</td><td>Linear</td><td>0.00694</td><td>0.623</td><td>+</td></tr><tr><td>Sr</td><td>Para</td><td>0.00707</td><td>0.461</td><td>-</td></tr><tr><td>U</td><td>Para</td><td>0.0151</td><td>0.816</td><td>+</td></tr><tr><td>Coralline red alga</td><td>Rhodophyta</td><td>Ba</td><td>Para</td><td>0.0159</td><td>0.733</td><td>+</td></tr><tr><td>Serpulid worm</td><td>Annelida</td><td>Ba</td><td>Linear</td><td>0.0182</td><td>0.211</td><td>+</td></tr><tr><td rowspan="3">Temperate coral</td><td rowspan="3">Cnidaria</td><td>Li</td><td>Linear</td><td>0.00222</td><td>0.165</td><td>-</td></tr><tr><td>Mg</td><td>Para</td><td>0.0318</td><td>0.0462</td><td>+</td></tr><tr><td>Zn</td><td>Linear</td><td>0.00321</td><td>0.0798</td><td>-</td></tr></tbody></table>

Measured Sr/Ca decreases significantly with decreasing pH<sub>sw</sub> in the serpulid worm and a suite of molluscs, including the oyster, soft clam, conch, whelk, and limpet. Sr/Ca significantly increases in the hard clam and pencil urchin, and significantly decreases in the oyster with increasing pH<sub>cf</sub>. The Sr/Ca of the serpulid worm and the calcitic oyster exhibit parabolic trends with decreasing pH<sub>sw</sub>604781529604781529 (R<sup>2</sup> = 0.285, *p* = 0.01; R<sup>2</sup> = 0.356, *p* = 0.03). The Sr/Ca of the soft clam, whelk and limpet decrease with decreasing pH<sub>sw</sub>; however, they differ in that Sr/Ca of the soft clam follows an inverse parabolic pattern, while Sr/Ca of the whelk and limpet decrease linearly (R<sup>2</sup> = 0.723, *p* = 0.0003; R<sup>2</sup> = 0.171, *p* = 0.0254). For pH<sub>cf</sub>, similar to pH<sub>sw</sub>, the molluscan Sr/Ca varies in different ways: Sr/Ca of the aragonitic hard clam decreases parabolically with pH<sub>cf</sub> (R<sup>2</sup> = 0.639; *p* = 0.001), while Sr/Ca of the low-Mg calcite oyster decreases linearly with pH<sub>cf</sub> (R<sup>2</sup> = 0.602; *p* = 0.001).

DISCUSSION

The data presented survey nine elemental ratios in 18 marine calcifying species and thus represent a unique opportunity to examine carbonate elemental geochemistry from the perspective of comparative organismal biology. One of the main strengths of this analysis is its revelation of large-scale patterning of element incorporation amongst different CaCO<sub>3</sub> mineralogies, species,

<table><caption>**TABLE 4** | Summary of AIC test results for significant relationships (*p* < 0.05) between element/Ca ratios and seawater carbonate ion concentrations across all organisms.</caption><thead><tr><th rowspan="2">Organism</th><th rowspan="2">Phylum</th><th rowspan="2">X/Ca</th><th colspan="4">Best-fit regression</th></tr><tr><th>TYPE</th><th>p</th><th>R<sup>2</sup></th><th>Relationship</th></tr></thead><tbody><tr><td>American lobster</td><td>Arthropoda</td><td>Ba</td><td>Para</td><td>0.00562</td><td>0.763</td><td>-</td></tr><tr><td>Gulf shrimp</td><td>Arthropoda</td><td>Zn</td><td>Para</td><td>0.0302</td><td>0.479</td><td>+</td></tr><tr><td>Blue mussel</td><td>Mollusca</td><td>Zn</td><td>Para</td><td>0.0414</td><td>0.539</td><td>-</td></tr><tr><td>American oyster</td><td>Mollusca</td><td>Sr</td><td>Para</td><td>0.00373</td><td>0.541</td><td>+</td></tr><tr><td rowspan="2">Hard clam</td><td rowspan="2">Mollusca</td><td>Ba</td><td>Linear</td><td>0.00365</td><td>0.632</td><td>+</td></tr><tr><td>Ba</td><td>Para</td><td>0.000406</td><td>0.784</td><td>-</td></tr><tr><td>Soft clam</td><td>Mollusca</td><td>Sr</td><td>Para</td><td>0.00106</td><td>0.66</td><td>-</td></tr><tr><td rowspan="2">Conch</td><td rowspan="2">Mollusca</td><td>Sr</td><td>Para</td><td>0.00391</td><td>0.354</td><td>-</td></tr><tr><td>Ba</td><td>Para</td><td>0.0171</td><td>0.257</td><td>-</td></tr><tr><td rowspan="3">Whelk</td><td rowspan="3">Mollusca</td><td>Mg</td><td>Linear</td><td>0.00817</td><td>0.245</td><td>+</td></tr><tr><td>Sr</td><td>Linear</td><td>0.0319</td><td>0.156</td><td>+</td></tr><tr><td>Ba</td><td>Para</td><td>0.0403</td><td>0.193</td><td>-</td></tr><tr><td rowspan="3">Limpet</td><td rowspan="3">Mollusca</td><td>B</td><td>Linear</td><td>0.00406</td><td>0.675</td><td>+</td></tr><tr><td>Sr</td><td>Linear</td><td>0.0271</td><td>0.201</td><td>+</td></tr><tr><td>U</td><td>Linear</td><td>0.00535</td><td>0.597</td><td>+</td></tr><tr><td rowspan="3">Pencil urchin</td><td rowspan="3">Echinodermata</td><td>B</td><td>Linear</td><td>&lt;0.0001</td><td>0.88</td><td>+</td></tr><tr><td>Ba</td><td>Para</td><td>&lt;0.00001</td><td>0.752</td><td>-</td></tr><tr><td>U</td><td>Para</td><td>&lt;0.0001</td><td>0.982</td><td>+</td></tr><tr><td>Purple urchin</td><td>Echinodermata</td><td>Ba</td><td>Linear</td><td>0.014</td><td>0.335</td><td>-</td></tr><tr><td>Coralline red alga</td><td>Rhodophyta</td><td>B</td><td>Para</td><td>0.00357</td><td>0.651</td><td>-</td></tr><tr><td rowspan="5">Serpulid worm</td><td rowspan="5">Annelida</td><td>Zn</td><td>Para</td><td>&lt;0.00001</td><td>0.648</td><td>+</td></tr><tr><td>Li</td><td>Linear</td><td>0.0422</td><td>0.287</td><td>-</td></tr><tr><td>Mg</td><td>Para</td><td>0.0107</td><td>0.554</td><td>-</td></tr><tr><td>Zn</td><td>Para</td><td>0.0227</td><td>0.258</td><td>-</td></tr><tr><td>Sr</td><td>Para</td><td>0.00799</td><td>0.309</td><td>+</td></tr><tr><td rowspan="3">Temperate coral</td><td rowspan="3">Cnidaria</td><td>Li</td><td>Linear</td><td>0.00208</td><td>0.147</td><td>-</td></tr><tr><td>Zn</td><td>Linear</td><td>0.000476</td><td>0.0982</td><td>-</td></tr><tr><td>U</td><td>Linear</td><td>0.0435</td><td>0.0524</td><td>+</td></tr></tbody></table>

and seawater carbonate chemistries. It also allows examination of the sensitivity of elemental partitioning in biogenic carbonates to changes in external (pH<sub>sw</sub>) and internal pH (pH<sub>cf</sub>; boron-isotope-based; Liu et al., 2020) conditions. For many of the investigated species, elemental data were only available for 2–3 individuals per pCO<sub>2</sub> treatment. Although some relationships may not be apparent due to small sample size, interesting and statistically significant relationships were still evident (Tables 3–5; red outlines, Supplementary Figures 10–45). Nevertheless, these relationships warrant further investigation due to the small sample sizes involved.

Precipitation rates are known to affect elemental partition coefficients in inorganically precipitated calcium carbonates to varying degrees (See Supplementary Table 5 and references therein). Ries et al. (2009) illustrated that the cultured species exhibit diverse responses in net calcification rate to increasing pCO<sub>2</sub> (Table 1). Though net calcification rates of whole
<table><caption>**TABLE 5** | Summary of AIC test results for significant relationships (*p* < 0.05) between element/Ca ratios and seawater pH and ¹¹B-derived calcifying fluid pH (indicated with bold italic text) across all organisms.</caption><thead><tr><th rowspan="2">Organism</th><th rowspan="2">Phylum</th><th rowspan="2">X/Ca</th><th rowspan="2">pH<sub>sw</sub> or pH<sub>cf</sub></th><th colspan="4">Best-fit regression</th></tr><tr><th>TYPE</th><th>p</th><th>R<sup>2</sup></th><th>Relationship</th></tr></thead><tbody><tr><td rowspan="3">American lobster</td><td rowspan="3">Arthropoda</td><td>Ba</td><td>pH<sub>sw</sub></td><td>Para</td><td>0.00562</td><td>0.763</td><td>-</td></tr><tr><td>Sr</td><td>pH<sub>sw</sub></td><td>Para</td><td>0.0284</td><td>0.356</td><td>+</td></tr><tr><td></td><td><b>pH<sub>cf</sub></b></td><td><b>Linear</b></td><td><b>0.000675</b></td><td><b>0.602</b></td><td>-</td></tr><tr><td rowspan="3">American oyster</td><td rowspan="3">Mollusca</td><td>Ba</td><td>pH<sub>sw</sub></td><td>Linear</td><td>0.00136</td><td>0.71</td><td>+</td></tr><tr><td>Sr</td><td>pH<sub>sw</sub></td><td>Para</td><td>0.0221</td><td>0.695</td><td>+</td></tr><tr><td></td><td><b>pH<sub>cf</sub></b></td><td>Linear</td><td>0.0308</td><td>0.326</td><td>+</td></tr><tr><td rowspan="3">Hard clam</td><td rowspan="3">Mollusca</td><td>Ba</td><td>pH<sub>sw</sub></td><td><b>Para</b></td><td><b>0.000625</b></td><td><b>0.963</b></td><td>+</td></tr><tr><td><b>B</b></td><td><b>pH<sub>cf</sub></b></td><td><b>Para</b></td><td><b>0.00146</b></td><td><b>0.639</b></td><td>+</td></tr><tr><td><b>Sr</b></td><td><b>pH<sub>cf</sub></b></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Soft clam</td><td rowspan="2">Mollusca</td><td>Sr</td><td>pH<sub>sw</sub></td><td>Para</td><td>0.000343</td><td>0.723</td><td>-</td></tr><tr><td>Cd</td><td>pH<sub>sw</sub></td><td>Linear</td><td>0.0491</td><td>0.577</td><td>-</td></tr><tr><td>Conch</td><td>Mollusca</td><td>Sr</td><td>pH<sub>sw</sub></td><td>Para</td><td>0.0156</td><td>0.263</td><td>-</td></tr><tr><td>Periwinkle</td><td>Mollusca</td><td>Ba</td><td>pH<sub>sw</sub></td><td>Linear</td><td>0.0482</td><td>0.155</td><td>-</td></tr><tr><td rowspan="2">Whelk</td><td rowspan="2">Mollusca</td><td>Mg</td><td>pH<sub>sw</sub></td><td>Linear</td><td>0.00441</td><td>0.283</td><td>+</td></tr><tr><td>Sr</td><td>pH<sub>sw</sub></td><td>Linear</td><td>0.0254</td><td>0.171</td><td>+</td></tr><tr><td rowspan="3">Limpet</td><td rowspan="3">Mollusca</td><td>B</td><td>pH<sub>sw</sub></td><td>Linear</td><td>0.00404</td><td>0.675</td><td>+</td></tr><tr><td>Sr</td><td>pH<sub>sw</sub></td><td>Linear</td><td>0.0261</td><td>0.204</td><td>+</td></tr><tr><td>U</td><td>pH<sub>sw</sub></td><td>Para</td><td>0.0082</td><td>0.674</td><td>+</td></tr><tr><td rowspan="3">Tropical urchin</td><td rowspan="3">Echinodermata</td><td>B</td><td>pH<sub>sw</sub></td><td>Linear</td><td>&lt;0.00001</td><td>0.905</td><td>+</td></tr><tr><td>Ba</td><td>pH<sub>sw</sub></td><td>Para</td><td>&lt;0.00001</td><td>0.833</td><td>-</td></tr><tr><td></td><td><b>pH<sub>cf</sub></b></td><td><b>Para</b></td><td><b>0.0086</b></td><td><b>0.502</b></td><td>-</td></tr><tr><td rowspan="2">Temperate urchin</td><td rowspan="2">Echinodermata</td><td>Sr</td><td><b>pH<sub>cf</sub></b></td><td><b>Para</b></td><td><b>0.000357</b></td><td><b>0.791</b></td><td>-</td></tr><tr><td>U</td><td>pH<sub>sw</sub></td><td>Para</td><td>&lt;0.0001</td><td>0.982</td><td>+</td></tr><tr><td>Coralline red alga</td><td>Rhodophyta</td><td>Ba</td><td>pH<sub>sw</sub></td><td>Linear</td><td>0.147</td><td>0.33</td><td>-</td></tr><tr><td></td><td></td><td>B</td><td>pH<sub>sw</sub></td><td>Para</td><td>0.00194</td><td>0.695</td><td>-</td></tr><tr><td></td><td></td><td></td><td><b>pH<sub>cf</sub></b></td><td><b>Linear</b></td><td><b>0.0038</b></td><td><b>0.681</b></td><td>+</td></tr><tr><td></td><td></td><td><b>Ba</b></td><td><b>pH<sub>cf</sub></b></td><td><b>Para</b></td><td><b>0.0119</b></td><td><b>0.472</b></td><td>-</td></tr><tr><td rowspan="4">Serpulid worm</td><td rowspan="4">Annelida</td><td>Mg</td><td>pH<sub>sw</sub></td><td>Linear</td><td>0.0306</td><td>0.44</td><td>+</td></tr><tr><td>Zn</td><td>pH<sub>sw</sub></td><td>Linear</td><td>0.0379</td><td>0.3</td><td>+</td></tr><tr><td>Sr</td><td>pH<sub>sw</sub></td><td>Linear*</td><td>&lt;0.00001</td><td>0.663</td><td>-</td></tr><tr><td></td><td><b>pH<sub>cf</sub></b></td><td><b>Para</b></td><td><b>0.000377</b></td><td><b>0.529</b></td><td>-</td></tr><tr><td rowspan="3">Temperate coral</td><td rowspan="3">Cnidaria</td><td><b>B</b></td><td><b>pH<sub>cf</sub></b></td><td><b>Linear</b></td><td><b>0.0161</b></td><td><b>0.528</b></td><td>-</td></tr><tr><td>Li</td><td>pH<sub>sw</sub></td><td>Para</td><td>0.0151</td><td>0.519</td><td>-</td></tr><tr><td>Zn</td><td>pH<sub>sw</sub></td><td>Para</td><td>0.0129</td><td>0.301</td><td>-</td></tr><tr><td></td><td></td><td></td><td><b>pH<sub>cf</sub></b></td><td><b>Linear</b></td><td><b>0.0147</b></td><td><b>0.447</b></td><td>+</td></tr><tr><td rowspan="3">Temperate coral</td><td rowspan="3">Cnidaria</td><td><b>B</b></td><td><b>pH<sub>sw</sub></b></td><td>Linear</td><td>0.00129</td><td>0.161</td><td>-</td></tr><tr><td>Zn</td><td>pH<sub>sw</sub></td><td>Linear</td><td>0.0000893</td><td>0.0884</td><td>-</td></tr><tr><td><b>Cd</b></td><td><b>pH<sub>cf</sub></b></td><td><b>Para</b></td><td><b>0.0125</b></td><td><b>0.5</b></td><td>+</td></tr><tr><td></td><td></td><td>U</td><td>pH<sub>sw</sub></td><td>Para</td><td>0.0356</td><td>0.0792</td><td>+</td></tr></tbody></table>

organisms cannot easily be related to the growth rates of calcium carbonate crystals, there are a number of significant relationships between X/Ca and the determined net calcification rate of the species across the pCO₂ conditions for 10 of the 18 species (Table 3). Mg/Ca, Sr/Ca, and Ba/Ca were the most common elemental ratios that varied significantly with net calcification rates across the 10 species, while Li/Ca, Zn/Ca, Cd/Ca, and U/Ca were the least common. B/Ca only significantly varied with net calcification rate in the pencil urchin, *E. tribuloides*. However, despite there being a number of significant relationships between X/Ca and net calcification rate in the cultured species, the

signs and shapes of many of the relationships differ from those published for inorganic minerals (e.g., Gabitov et al., 2014; Holcomb et al., 2016; Mavromatis et al., 2018).

Boron was among the most responsive elements to seawater carbonate chemistry, showing statistically significant trends in five of the 18 species. For example, the B/Ca in the pencil urchin was highly sensitive (over 200 µmol/mol, 62% of the mean, across treatments; Supplementary Figures 34, 35) to manipulation of seawater carbonate chemistry, while B/Ca in the purple urchin was not (Supplementary Figures 38, 39). B/Ca in the coralline red algae was also highly sensitive to
<table><caption>**TABLE 6** | Summary of the simplest (X/Ca vs. Mg/Ca) as well as the final GAMs (bolded) as determined by likelihood ratio tests results for the X/Ca vs. Mg/Ca scatterplots (Figure 2).</caption><thead><tr><th>Elemental ratio</th><th>Formula</th><th>Bootstrapped R<sup>2</sup></th><th>95% confidence interval<br>(lower, higher)</th><th>% deviance explained</th><th>Model likelihood</th><th>n</th></tr></thead><tbody><tr><td rowspan="2">Li/Ca</td><td>Li.Ca ~ Mg.Ca</td><td>0.65</td><td>0.56, 0.72</td><td>64.60%</td><td>678.43</td><td>181</td></tr><tr><td>Li.Ca ~ Mg.Ca *<br>Phylum +<br>Carbonate.Material</td><td>0.96</td><td>0.93, 0.97</td><td>96%</td><td>482.14</td><td>181</td></tr><tr><td rowspan="2">B/Ca</td><td>B.Ca ~ Mg.Ca</td><td>0.06</td><td>0.02, 0.10</td><td>5.54%</td><td>1302</td><td>188</td></tr><tr><td>B.Ca ~ Mg.Ca *<br>Phylum</td><td>0.98</td><td>0.98, 0.99</td><td>98.20%</td><td>927.98</td><td>188</td></tr><tr><td rowspan="2">Zn/Ca</td><td>Zn.Ca ~ Mg.Ca</td><td>0.03</td><td>0.01, 0.05</td><td>2.89%</td><td>1928.7</td><td>279</td></tr><tr><td>Zn.Ca ~ Mg.Ca +<br>Phylum</td><td>0.62</td><td>0.49, 0.70</td><td>61.70%</td><td>1799</td><td>279</td></tr><tr><td rowspan="2">Sr/Ca</td><td>Sr.Ca ~ Mg.Ca</td><td>n.s.</td><td></td><td></td><td></td><td></td></tr><tr><td>Sr.Ca ~ Mg.Ca *<br>Phylum +<br>Carbonate.Material</td><td>0.99</td><td>0.99, 0.99</td><td>99%</td><td>85.993</td><td>239</td></tr><tr><td rowspan="2">Cd/Ca</td><td>Cd.Ca ~ Mg.Ca</td><td>0.72</td><td>0.51, 0.83</td><td>71.90%</td><td>-12.355</td><td>148</td></tr><tr><td>Cd.Ca ~<br>Carbonate.Material</td><td>0.51</td><td>0.39, 0.60</td><td>50.90%</td><td>29.034</td><td>148</td></tr><tr><td rowspan="2">Ba/Ca</td><td>Ba.Ca ~ Mg.Ca</td><td>0.17</td><td>0.10, 0.27</td><td>17.00%</td><td>1281.5</td><td>270</td></tr><tr><td>Ba.Ca ~ Mg.Ca *<br>Phylum</td><td>0.69</td><td>0.54, 0.79</td><td>68.70%</td><td>1149.9</td><td>270</td></tr><tr><td rowspan="2">U/Ca</td><td>U.Ca ~ Mg.Ca</td><td>n.s.</td><td></td><td></td><td></td><td></td></tr><tr><td>U.Ca ~ Mg.Ca *<br>Phylum +<br>Carbonate.Material</td><td>0.93</td><td>0.87, 0.95</td><td>93%</td><td>1063.6</td><td>175</td></tr><tr><td rowspan="2">Na/Ca</td><td>Na.Ca ~ Mg.Ca</td><td>n.s.</td><td></td><td></td><td></td><td></td></tr><tr><td>Na.Ca ~ Mg.Ca *<br>Phylum +<br>Carbonate. Material</td><td>0.39</td><td>0.15, 0.56</td><td>38.70%</td><td>383.96</td><td>121</td></tr></tbody></table>

In the case of Cd/Ca, the simplest model was also determined to be the most likely. Model fit statistics have been bootstrapped over 1,000 iterations. The models tested explored the influence of considering grouping the data by phylum and/or mineralogy. Interactive model components are denoted by an asterisk (*) while additive model components are denoted by a plus sign. Note that "Phylum" here does not account for the phylogenetic relatedness of the species as we do in **Figure 7** and **Table 2** and **Supplementary Table 9** (see methods 2.5.2). For the full list of candidate models, see **Supplementary Table 4**, and model comparisons in **Supplementary Table 5**. Models for Cd/Ca and Na/Ca may not be accurate due to the large amount of data missing for five and seven organisms, respectively.

seawater [CO₃²⁻] (Supplementary Figures 24, 25). U/Ca was also found to be sensitive to seawater [CO₃²⁻] in the temperate coral (Supplementary Figure 23), the limpet (Supplementary Figure 32) and the pencil urchin (Supplementary Figure 34). Significant trends were also observed for many other elements, including Mg, Ba, Sr, Zn, and U. In some cases, significant relationships between element/Ca ratios and pH<sub>cf</sub> derived from boron isotopes (Liu et al., 2020) were apparent when relationships with external pH<sub>sw</sub> were not (e.g., Hard clam; Table 5 and Supplementary Figures 30, 31). These data suggest that seawater carbonate chemistry influences X/Ca ratios in a range of calcifying species, provide launching points for exploring the viability of these elements as palaeoceanographic proxies of seawater carbonate chemistry.

Regarding large-scale patterns in the study, most elements showed significant variability amongst species, and apparent biomineral partition coefficients were often quite different from empirically derived elemental partition coefficients for calcite, aragonite, or ACC precipitated inorganically from seawater (Table 2 and Figure 3). Deviations from inorganic mineral

relationships reveal the influence that different biomineralization mechanisms have on the pathways that ions take to the site of calcification. Some clear patterns emerged in this regard. For example, Mg/Ca in mollusc shells was usually low regardless of shell polymorph mineralogy, and Sr/Ca was low in all the biominerals compared to ratios in inorganic precipitates, with coral and *Halimeda* aragonite exhibiting higher Sr/Ca compared to other species, and aragonite molluscs exhibiting especially low Sr/Ca. Strikingly, mollusc CaCO₃ (of all polymorphs) was relatively low in all X/Ca ratios compared to other species, except for Li/Ca (compared to the temperate coral) and Na/Ca.

The partition coefficients for some elements in mollusc shells, including for Li and B, were notably similar to recently published elemental data for inorganically precipitated ACC (e.g., K<sub>Li</sub> in ACC: 6.5E-3–14.5E-3; K<sub>Li</sub> in the oyster: 9.5E-3; Table 2; Evans et al., 2020), regardless of mollusc polymorph mineralogy (Figure 3). This is noteworthy as molluscs are known to form their shells from ACC precursors (Addadi et al., 2003), whereas the evidence and importance of ACC in some of the other species examined here are unclear. However, it should be noted
<table><caption>**TABLE 7** | Estimates of phylogenetic signal for X/Ca ratios.</caption><thead><tr><th>Elemental Ratio</th><th>Lambda</th></tr></thead><tbody><tr><td><strong>Multivariate</strong></td><td></td></tr><tr><td>All</td><td>0.99</td></tr><tr><td>Li/Ca, Mg/Ca, Zn/Ca, and Sr/Ca</td><td>0.895</td></tr><tr><td>Li/Ca, Mg/Ca, Sr/Ca</td><td>0.898</td></tr><tr><td>Mg/Ca and Sr/Ca</td><td>0.925</td></tr><tr><td><strong>Univariate</strong></td><td></td></tr><tr><td>Li/Ca</td><td>0</td></tr><tr><td>B/Ca</td><td>0.488</td></tr><tr><td>Mg/Ca</td><td>0.803</td></tr><tr><td>Zn/Ca</td><td>0</td></tr><tr><td>Sr/Ca</td><td>0.835</td></tr><tr><td>Cd/Ca</td><td>0.665</td></tr><tr><td>Ba/Ca</td><td>0</td></tr><tr><td>U/Ca</td><td>0.889</td></tr><tr><td>Na/Ca</td><td>0.516</td></tr></tbody></table>

Multivariate analyses were conducted using all and combinations of the X/Ca ratios together, as based on a phylogenetic principal component analysis in phytools (Revell, 2012; row = "All"). Univariate analyses estimated phylogenetic signals for individual X/Ca ratios using the phylosig function also in phytools.

that ACC precursor phases may dissolve and reprecipitate in their transformation to crystalline calcium carbonate, potentially involving the exchange of elements with the surrounding fluid (e.g., Giuffre et al., 2015). Since most of the apparent partition coefficients do not resemble those reported for ACC, it is not known how pertinent the elemental data for ACC (e.g., Evans et al., 2020) are to crystalline biominerals. Nevertheless, the data in the present study provide an interesting point of comparison.

Varying levels of biological control on biomineralization are apparent across the phyla studied when assessing elements individually (Figures 1, 3; Dove, 2010). However, clear groupings of species arise when evaluating the elemental ratios together (Figures 4–6). A possible explanation is that species most enriched in Mg tend to also exhibit high ratios of other elements (Figure 2). This is most apparent for the organisms that produce high-Mg calcite: the arthropods, echinoderms, coralline red algae, and serpulid worms. This scenario has been demonstrated in synthetic calcites, where higher incorporation of Mg²⁺ led to higher incorporation of Sr²⁺ due to distortion of the crystal lattice (Mucci and Morse, 1983). Although this scenario may explain the behavior of some of the measured elements (e.g., Li/Ca; Figure 2), many cannot be explained in this manner. To explore potential explanations for the variability in these patterns, the GAMs demonstrate that accounting for phylum and/or mineralogy substantially improves the models' power in predicting X/Ca in the studied species (Table 6). Furthermore, when analyzed in an evolutionary context (i.e., time-calibrated phylogenetic tree), X/Ca ratios exhibit strong phylogenetic signals, with the differences in elemental ratios between species largely reflecting phylogenetic distance (Table 2). Clearly, Mg/Ca alone cannot explain the abundance of other elements via lattice distortion, and our GAMs provide evidence that mineralogy and phylum are key predictors of elemental ratios in the studies species. Thus, explanations for differences in

elemental ratios rooted in differences in the species' mechanisms of biomineralization must be explored.

Calcitic bivalve molluscs are known to strongly and actively exclude Mg from their biominerals to form low-Mg calcite, as shown specifically for the blue mussel (M. edulis; Lorens and Bender, 1977). The exclusion of trace and minor elements from molluscan biominerals is likely associated with regulation of the ionic composition of the parent fluid for calcification, as well as potentially other processes, such as detoxification (L. littorea; Mason and Nott, 1981). The data presented here suggest that this or similar processes are also active in these bivalves and gastropods, and that other elements are also excluded during the mollusc biomineralization process, given the low elemental ratios and low apparent partition coefficients compared to those of inorganically precipitated calcite (Figures 1, 3).

Amongst the high-Mg calcitic species investigated, comparisons between the arthropods and echinoderms are particularly revealing. Both exhibit similar Mg/Ca, but the crustaceans are relatively enriched in the other elemental ratios. Echinoderms and arthropods form their calcareous exoskeletons by transporting amorphous calcium carbonate to be co-precipitated with an organic matrix that appears to serve as a mineralogical template, but has different B, Ba, and Sr compositions (Roer and Dillaman, 1984; Märkel et al., 1986; Beniash et al., 1997; Dillaman et al., 2005). The high-Mg content can be explained by way of an amorphous precursor, as ACC has been shown to incorporate up to 20 mol% Mg during its transformation into Mg-calcite (Purgstaller et al., 2016); however, the source of discrepancy between the other elements is less clear (see reviews by Long et al. (2014) and Evans et al. (2020)). Some crustacea, including the American lobster evaluated here, maintain a stabilized ACC phase as a predominant component of their fully calcified cuticle in addition to the layer of high-Mg calcite in the exocuticle (Levi-Kalisman et al., 2002; Mergelsberg et al., 2019; Evans et al., 2020). Arthropods are also unique in their ability to bi-directionally transport Ca²⁺ and CO₃²⁻ when resorbing the old cuticle during pre-molting, which could lead to element accumulation and possibly explain the enrichment in X/Ca ratios relative to the echinoderms.

Echinoderms and arthropods also differ in the extent to which they isolate their site of calcification from seawater and utilize different types of cells to transport ACC to the site of calcification—primary mesenchyme cells (PMC) in echinoderms and epithelial cells in the arthropods (Roer and Dillaman, 1984; Märkel et al., 1986). Echinoderm PMCs transport ACC to the site of calcification and form a syncytium of other cells that maintain favorable calcification conditions. The pH within PMCs is maintained at approximately 6.9; however, the pH<sub>cf</sub> within the syncytium is typically equivalent to that of seawater (Stumpp et al., 2012). Despite this, [HCO₃⁻] in the syncytium is maintained at higher levels relative to seawater (Holtmann, 2013). Therefore, these fundamentally different modes of biomineralization in arthropods and echinoderms may both support the formation of high-Mg calcite, but with differing compositions of other elements.

The species that exhibit the highest Mg/Ca ratios, the serpulid worm and the coralline red alga, also exhibit the highest levels
of the other X/Ca ratios. Here, too, there are no significant relationships of any X/Ca ratio with Mg/Ca (Figure 2). Given that these two species are phylogenetically distant, one an animal and the other an alga, it is interesting that they cluster together in the multivariate analyses (Figures 4, 5). Figure 3 demonstrates that the apparent partition coefficients calculated for the serpulid worm typically either agree with the published values for inorganic calcite or aragonite or lie in the middle of the two (Table 2). This may indicate that, although generally enriched in X/Ca ratios, the differences between the serpulid worm and the coralline red algae potentially arise from the bimineralic composition of the serpulid worm tube sampled here. Serpulid worms and coralline red algae, including the species investigated here, have been shown to elevate pH at the site of calcification (Chan et al., 2015a; Short et al., 2015; Cornwall et al., 2017; Donald et al., 2017; Sutton et al., 2018). Although these processes may support calcification, they do not explain the high X/Ca ratios measured in the coralline red algae. Coralline red algae biomineralize intercellularly, outside of the cell membrane but within and between cell walls (Ries, 2006; Long et al., 2014). Evidence of ACC has been observed at the growing edge of a coralline red algal biomineral (Raz et al., 2000; Rahman et al., 2019). Stabilized ACC does not dominate the biomineral phase produced by coralline red algae, as it does for the lobster, but the potential formation of coralline red algal high-Mg calcite via ACC could explain the very high Mg component (Purgstaller et al., 2016; Mergelsberg et al., 2019; Evans et al., 2020). It has also been shown that various biomolecules that compose the organic matrices and templates, in conjunction with Mg²⁺ ions in solution, can lead to higher levels of Mg²⁺-incorporation in calcite (Long et al., 2014).

The green alga *H. incrassata* and the temperature coral *O. arbuscula* group together in the hierarchical clustering and multivariate analyses of the measured elements, but this may simply be driven by their shared aragonitic mineralogy. Evidence for pH elevation at the site of calcification has also been shown for this coral species and this and other species of *Halimeda* green algae (Liu et al., 2020; see review by Borowitzka and Larkum, 1987). Both of the species photosynthesize and precipitate aragonite biominerals (Ries et al., 2009; Ries, 2011), but that is generally where their similarities in biomineralization end. For the *Halimeda* alga, precipitation occurs extracellularly within invaginated sections of a syncytial cell wall (Borowitzka and Larkum, 1987; Ries, 2009; Cuif et al., 2013). Conversely, the *Oculina* coral apparently mineralizes in fluid pockets beneath the calicoblastic epithelium (Tambutté et al., 2011). Though there is currently debate surrounding the precise mechanisms of ion transport and the potential role of organic templates for mineral nucleation in scleractinian corals, there is general agreement that the conditions at the site of calcification and the mineralization process are highly controlled (Ogilvie, 1896; Cohen and McConnaughey, 2003; Mass et al., 2017; Von Euw et al., 2017).

Qualitatively, it does not appear that elemental chemistry perfectly recapitulates taxonomy due to the confounding influence of calcium carbonate polymorph (Figures 4B-D). However, the dendrogram produced by the hierarchical

clustering analysis of the elemental chemistry of the studied species (Figure 4A) visually recapitulates phylogenetic relationships amongst the aragonitic (Figure 4B) and calcitic species (Figure 4C). However, mixed mineralogy whelk and serpulid worm do not exhibit this taxonomic clustering. Even within the molluscs, the bivalve and gastropod molluscs groups are distinct in Figure 4A, potentially revealing clustering at the class level. A similar class-level clustering for P/Ca and Mg/Ca has also been identified in the Malacostraca (Mergelsberg et al., 2019). The whelk is further removed from the other aragonitic gastropods in the hierarchical cluster analysis, most likely due to its bimineralic nature (~75% aragonite: 25% calcite; Ries, 2011). The dislocation of the serpulid worm tube in the hierarchical clustering analysis, which is relatively closely related to the molluscs, probably arised from its bimineralic composition (ca. 50% aragonite: 50% calcite; Ries, 2011). The aragonitic coral and green algae are the least related to everything else in the aragonitic tree and are also closest in the clustering analysis of the elemental data. Collectively, these results suggest that the divergence in X/Ca ratios between species is related to both phylogeny and mineralogy.

The wide array of biomineralization pathways amongst the studied species makes it difficult to deconvolve the effects of biological processes from physicochemical processes. Although the strong phylogenetic signal detected across all X/Ca ratios using phylogenetic analyses is partly driven by patterns of elemental congruence within aragonitic and calcitic groups (Figures 4A-D), using both phylogenetic and non-phylogenetic regression models to examine interspecific differences in X/Ca ratios show that phylogeny explains more of the variance in the X/Ca ratios than mineralogy (Figure 7 and Supplementary Table 9).604781538604781538

NMDS ordination visualizes this primary clustering of organisms by phyla as well as the secondary clustering by mineralogy (Figures 5, 6). Clustering by mineralogy is not unexpected, as differences in X/Ca ratios between calcite and aragonite are well documented via inorganic precipitation experiments (e.g., Gabitov et al., 2011, 2019; Mavromatis et al., 2015, 2018). Furthermore, the phylogenetic analyses only including the Molluscan data show that clustering of elemental ratios by phylogenetic proximity is relevant even at the class level (see section “Phylogenetic Signal and Regression Analyses”). Future studies should examine whether the apparent accord between phylogeny and elemental clustering at the class level within Mollusca holds true when including more species, including ones with mixed mineralogies. Regardless, an interesting and new result of this work is the strong phylogenetic signal within the elemental clustering of the studied species, even after controlling for the effects of mineralogy (Figures 3-7).

The present study has implications for interpreting carbonate preservation and diagenesis in the sedimentary record. By systematically defining the plausible ranges and patterns of elemental partitioning in a diverse range of modern calcifiers, for example by multivariate ordination analysis (e.g., Figures 5, 6), measurements of paleo-carbonates can be placed within this framework to assess whether they fall within the range of compositions of evolutionarily-related and mineralogically
similar organisms. If they do not, it may indicate diagenetic alteration of the elemental composition of the fossil—rendering it unsuitable for paleo-reconstructions.

In sum, it is challenging to unequivocally attribute specific aspects of species' biomineralization pathways to differing patterns of elemental partitioning in biogenic carbonates. The present study shows that the broad patterning of elemental partitioning within marine calcifiers reflects not only environmental conditions, but also phylogeny, potentially via the molecular and structural mechanisms that they evolved to support the formation of their protective shells and skeletons. With more work, mechanistic explanations of elemental partitioning in biogenic carbonates should be achievable.

## 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.

## AUTHOR CONTRIBUTIONS

RE designed the project, carried out the trace element analyses, and advised students and the postdoctoral researchers who contributed to the project on data analysis and manuscript drafting. RU wrote the manuscript with input from RE, as well as from MG, JS, CR-P, HC, AT, and JR. JR carried out the culturing experiment and characterization of polymorph mineralogy, calcification rates, and carbonate seawater chemistry. RH and SS surveyed the literature and curated the list of synthetic carbonate partition coefficients with input from RU and MG. AH constructed the SI plots of element-to-calcium ratios vs. net calcification rate and measured carbonate chemistry parameters and conducted the AIC analyses. JS conducted the hierarchical clustering analysis and constructed the NMDS ordination with all

organisms. JC constructed the calcite- and aragonite-only NMDS ordinations with input from RU. HC consulted on much of the statistics employed and constructed the figures showing the final GAM components. CR-P constructed the time-calibrated phylogenetic tree and conducted the phylogenetic analyses with input from RU. All authors contributed to the article and approved the submitted version.

## FUNDING

RU was funded by an NSF GRFP DGE-1650604 and a fellowship from UCLA's Center for Diverse Leadership in Science. MG was supported by DOE BES grant (grant no. DE-FG02-13ER16402). WD was supported an award from the "Laboratoire d'Excellence" LabexMER (ANR-10-LABX-19) and co-funded by a grant from the French government under the program "Investissements d'Avenir" as well as NSF grant OCE-1437166 to RE. SS was supported by the UCLA CARE SEM program. JR acknowledges support from NSF grant OCE-1437371. RE acknowledges support from the Pritzker Endowment to UCLA IoES. HC was supported through a postdoctoral fellowship by the Institutional Research and Academic Career Development Awards (IRACDA) program at UCLA (Award #K12 GM106996).

## ACKNOWLEDGMENTS

We also thank Mervyn Greaves (Cambridge) for assistance with the trace element analyses.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart. 2021.641760/full#supplementary-material

## REFERENCES

* Addadi, L., Joester, D., Nudelman, F., and Weiner, S. (2006). Mollusc shell formation: a source of new concepts for understanding biomineralization processes. *Chem. Eur. J.* 12, 980–987. doi: 10.1002/chem.200500980

* Addadi, L., Raz, S., and Weiner, S. (2003). Taking advantage of disorder: amorphous calcium carbonate and its roles in biomineralization. *Adv. Mater.* 15, 959–970. doi: 10.1002/adma.200300381

* Allen, K. A., Hönisch, B., Eggins, S. M., Haynes, L. L., Rosenthal, Y., and Yu, J. (2016). Trace element proxies for surface ocean conditions: a synthesis of culture calibrations with planktic foraminifera. *Geochem. Cosmochim. Acta* 193, 197–221. doi: 10.1016/j.gca.2016.08.015

* Allison, N., Finch, A. A., Sutton, S. R., and Newville, M. (2001). Strontium heterogeneity and speciation in coral aragonite: implications for the strontium paleothermometer. *Geochem. Cosmochim. Acta* 65, 2669–2676. doi: 10.1016/S0016-7037(01)00628-7

* Anderson, D., and Burnham, K. (2004). *Model Selection and Multi-Model Inference*, 2nd Edn, Vol. 63. New York, NY: Springer-Verlag, 10.

* Barker, S., Greaves, M., and Elderfield, H. (2003). A study of cleaning procedures used for foraminiferal Mg/Ca paleothermometry. *Geochem. Geophys. Geosyst.* 4:8407. doi: 10.1029/2003GC000559

* Beniash, E., Aizenberg, J., Addadi, L., and Weiner, S. (1997). Amorphous calcium carbonate transforms into calcite during sea urchin larval spicule growth. *Proc. R. S. B Biol. Sci.* 264, 461–465.

* Beretta, L., and Santaniello, A. (2016). Nearest neighbor imputation algorithms: a critical evaluation. *BMC Med. Inform. Decis. Mak.* 16(Suppl 3):74. doi: 10.1186/s12911-016-0318-z

* Borowitzka, M. A., and Larkum, A. W. D. (1987). Calcification in algae: mechanisms and the role of metabolism. *Crit. Rev. Plant Sci.* 6, 1–45. doi: 10.1080/07352688709382246

* Boyle, E. A., Huested, S. S., and Jones, S. P. (1981). On the distribution of copper, nickel, and cadmium in the surface waters of the North Atlantic and North Pacific Ocean. *J. Geophys. Res.* 86, 8048–8066. doi: 10.1029/jc086ic09p08048

* Brand, U., and Veizer, J. (1980). Chemical diagenesis of a multicomponent carbonate system – 1: trace Elements. *J. Sediment. Petrol.* 50, 1219–1236. doi: 10.1306/212f7df6-2b24-11d7-8648000102c1865d

* Branson, O. (2018). "Boron incorporation into marine CaCO₃," in *Advances in Isotope Geochemistry*, eds H. Marschall and G. Foster (Cham: Springer), 71–105. doi: 10.1007/978-3-319-64666-4_4

* Branson, O., Redfern, S. A. T., Elmore, A. C., Read, E., Valencia, S., and Elderfield, H. (2018). The distribution and coordination of trace elements in Krithe
ostracods and their implications for paleothermometry. *Geochim. Cosmochim. Acta* 236, 230–239. doi: 10.1016/j.gca.2017.12.005

Brecevic, L., and Nielsen, A. E. (1989). Solubility of amorphous calcium carbonate. *J. Cryst. Growth* 98, 504–510.

Broecker, W. S., and Peng, T. H. (1983). *Tracers in the SEA*. Palisades, NY, Eldigio Press, 23–24.

Bruland, K. (1983). *Trace Elements in Seawater. Chemical Oceanography*, 2nd Edn, Vol. 8, eds J. P. Riley and R. Chester (London: Academic).

Burdett, J. W., Arthur, M. A., and Richardson, M. (1989). A Neogene seawater sulfur isotope age curve from calcareous pelagic microfossils. *Earth Planet. Sci. Lett.* 94, 189–198. doi: 10.1016/0012-821X(89)90138-6

Canty, A., and Ripley, B. (2020). *boot: Bootstrap R (S-Plus) Functions*. R package version 1.3-25.

Chan, V. B. S., Toyofuku, T., Wetzel, G., Saraf, L., Thiyagarajan, V., and Mount, A. S. (2015a). Direct deposition of crystalline aragonite in the controlled biomineralization of the calcareous tubeworm. *Front. Mar. Sci.* 2:97. doi: 10.3389/fmars.2015.00097

Chan, V. B. S., Vinn, O., Li, C., Lu, X., Kudryavtsev, A. B., Schopf, J. W., and Thiyagarajan, V. (2015b). Evidence of compositional and ultrastructural shifts during the development of calcareous tubes in the biofouling tubeworm, Hydroides elegans. *Journal of Structural Biology*, 189(3), 230–237. doi: 10.1016/j.jsb.2015.01.004

Chavé, K. E. (1954). Aspects of the biogeochemistry of magnesium 1. Calcareous marine organisms. *J. Geol.* 62, 266–283. doi: 10.1086/626162

Cohen, A. L., and McConnaughey, T. A. (2003). Geochemical perspectives on coral mineralization. *Rev. Mineral. Geochem.* 54, 151–187. doi: 10.2113/0540151

Constantz, B. R. (1986). Coral skeleton construction: a physicochemically dominated process. *Palaios* 1, 152–157. doi: 10.2307/3514508

Constantz, B. R., and Meike, A. (1989). "Origin, evolution, and modern aspects of biomineralization in plants and animals," in *Origin, Evolution, and Modern Aspects of Biomineralization in Plants and Animals*, ed. R. E. Crick (New York, NY: Springer Science+Business Media), 201–208. doi: 10.1007/978-1-4757-6114-6

Cornwall, C. E., Comeau, S., and McCulloch, M. T. (2017). Coralline algae elevate pH at the site of calcification under ocean acidification. *Glob. Change Biol.* 23, 4245–4256. doi: 10.1111/gcb.13673

Cuif, J.-P., Dauphin, Y., and Sorauf, J. E. (2013). *Biominerals and Fossils Through Time. Journal of Chemical Information and Modeling*, Vol. 53. Cambridge: Cambridge University Press, doi: 10.1017/CBO9781107415324.004

Dang, H., Wang, T., Qiao, P., Bassinot, F., and Jian, Z. (2019). The B/Ca and Cd/Ca of a subsurface-dwelling foraminifera Pulleniatina obliquiluculata in the tropical Indo-Pacific Ocean: implications for the subsurface carbonate chemistry estimation. *Acta Oceanol. Sin.* 38, 138–150. doi: 10.1007/s13131-019-1406-6

Davison, A. C., and Hinkley, D. V. (1997). *Bootstrap Methods and Their Applications*. Cambridge: Cambridge University Press.

Dawber, C. F., and Tripati, A. (2012). Relationships between bottom water carbonate saturation and element/Ca ratios in coretop samples of the benthic foraminifera Oridorsalis umbonatus. *Biogeosciences* 9, 3029–3045. doi: 10.5194/bg-9-3029-2012

De Choudens-Sánchez, V., and González, L. A. (2009). Calcite and aragonite precipitation under controlled instantaneous supersaturation: elucidating the role of CACO 3 saturation state and Mg/Ca ratio on calcium carbonate polymorphism. *J. Sediment. Res.* 79, 363–376. doi: 10.2110/jsr.2009.043

de Villiers, S., Greaves, M., and Elderfield, H. (2002). An intensity ratio calibration method for the accurate determination of Mg/Ca and Sr/Ca of marine carbonates by ICP-AES. *Geochem. Geophys. Geosyst.* 3, 623–634. doi: 10.1029/2001gc000169

De Yoreo, J. J., Gilbert, P. U. P. A., Sommerdijk, N. A. J. M., Penn, R. L., Whitelam, S., Joester, D., et al. (2015). Crystallization by particle attachment in synthetic, biogenic, and geologic environments. *Science* 349:aaa6760. doi: 10.1126/science.aaa6760

DeCarlo, T. M., Gaetani, G. A., Holcomb, M., and Cohen, A. L. (2015). Experimental determination of factors controlling U/Ca of aragonite precipitated from seawater: implications for interpreting coral skeleton. *Geochim. Cosmochim. Acta* 162, 151–165. doi: 10.1016/j.gca.2015.04.016

Delaney, M. L., and Boyle, E. A. (1985). Li, Sr, Mg, and Na in foraminiferal calcite shells from laboratory culture, sediment traps, and sediment cores. *Deep Sea Res. Part B Oceanogr. Lit. Rev.* 32:1025. doi: 10.1016/0198-0254(85)93853-1

Delaney, M. L., Popp, B. N., Lepzelter, C. G., and Anderson, T. F. (1989). Lithium-to-calcium ratios in modern Cenozoic, and Paleosoic articulate brachiopod shells. *Paleoceanography* 4, 681–691.

Dellinger, M., West, A. J., Paris, G., Adkins, J. F., Pogge von Strandmann, P. A. E., Ullmann, C. V., et al. (2018). The Li isotope composition of marine biogenic carbonates: patterns and mechanisms. *Geochim. Cosmochim. Acta* 236, 315–335. doi: 10.1016/j.gca.2018.03.014

DePaolo, D. J. (2011). Surface kinetic model for isotopic and trace element fractionation during precipitation of calcite from aqueous solutions. *Geochim. Cosmochim. Acta* 64, 67–72. doi: 10.1016/j.gca.2010.11.020

Dillaman, R., Hequembourg, S., and Gay, M. (2005). Early pattern of calcification in the dorsal carapace of the blue crab, Callinectes sapidus. *J. Morphol.* 263, 356–374. doi: 10.1002/jmor.10311

Donald, H. K., Ries, J. B., Stewart, J. A., Fewell, S. E., and Foster, G. L. (2017). Boron isotope sensitivity to seawater pH change in a species of Neogoniolithon coralline red alga. *Geochim. Cosmochim. Acta* 217, 240–253. doi: 10.1016/j.gca.2017.08.021

Dove, P. M. (2010). The rise of skeletal biominerals. *Elements* 6, 37–42. doi: 10.2113/gselements.6.1.37

Drake, J. L., Mass, T., Stolarski, J., Von Euw, S., van de Schootbrugge, B., and Falkowski, P. G. (2020). How corals made rocks through the ages. *Glob. Change Biol.* 26, 31–53. doi: 10.1111/gcb.14912

Elderfield, H., and Rickaby, R. E. M. (2000). Oceanic Cd/P ratio and nutrient utilization in the glacial Southern Ocean. *Nature* 405, 305–310. doi: 10.1038/35012507

Elderfield, H., Bertram, C. J., and Erez, J. (1996). A biomineralization model for the incorporation of trace elements into foraminiferal calcium carbonate. *Earth Planet. Sci. Lett.* 142, 409–423. doi: 10.1016/0012-821X(96)00105-7

Endrizzi, F., and Rao, L. (2014). Chemical speciation of uranium(VI) in marine environments: complexation of calcium and magne- sium ions with [(UO2)(CO3)3]4? and the effect on the extraction of uranium from seawater. *Chem. Eur. J.* 20, 14499–14506.

Erez, J. (2003). The source of ions for biomineralization in foraminifera and their implications for paleoceanographic proxies. *Rev. Mineral. Geochem.* 54, 115–149. doi: 10.2113/0540115

Evans, D., Gray, W. R., Rae, J. W. B., Greenop, R., Webb, P. B., Penkman, K., et al. (2020). Trace and major element incorporation into amorphous calcium carbonate (ACC) precipitated from seawater. *Geochim. Cosmochim. Acta* 290, 293–311. doi: 10.1016/j.gca.2020.08.034

Farmer, J. R., Branson, O., Uchikawa, J., Penman, D. E., Hönisch, B., and Zeebe, R. E. (2019). Boric acid and borate incorporation in inorganic calcite inferred from B/Ca, boron isotopes and surface kinetic modeling. *Geochim. Cosmochim. Acta* 244, 229–247. doi: 10.1016/j.gca.2018.10.008

Feng, J., Lee, J. Y., Reeder, R. J., and Phillips, B. L. (2006). Observation of bicarbonate in calcite by NMR spectroscopy. *Am. Mineral.* 91, 957–960. doi: 10.2138/am.2006.2206

Finch, A. A., and Allison, N. (2008). Mg structural state in coral aragonite and implications for the paleoenvironmental proxy. *Geophys. Res. Lett.* 35:L08704. doi: 10.1029/2008GL033543

Frieder, C. A., Gonzalez, J. P., and Levin, L. A. (2014). Uranium in larval shells as a barometer of molluscan ocean acidification exposure. *Environ. Sci. Technol.* 48, 6401–6408. doi: 10.1021/es500514j

Gabitov, R. I., Sadekov, A., and Leinweber, A. (2014). Crystal growth rate effect on Mg/Ca and Sr/Ca partitioning between calcite and fluid: an in situ approach. *Chem. Geol.* 367, 70–82. doi: 10.1016/j.chemgeo.2013.12.019

Gabitov, R. I., Schmitt, A. K., Rosner, M., McKeegan, K. D., Gaetani, G. A., Cohen, A. L., et al. (2011). In situ $^{87}$Li, Li/Ca, and Mg/Ca analyses of synthetic aragonites. *Geochem. Geophys. Geosystems* 12, 1–16. doi: 10.1029/2010GC003322

Gabitov, R., Sadekov, A., Yapaskurt, V., Borrelli, C., Bychkov, A., Sabourin, K., et al. (2019). Elemental uptake by calcite slowly grown from seawater solution: an In-situ study via depth profiling. *Front. Earth Sci.* 7:51. doi: 10.3389/feart.2019.00051
Gaetani, G. A., Cohen, A. L., Wang, Z., and Crusius, J. (2011). Rayleigh-based, multi-element coral thermometry: a biomineralization approach to developing climate proxies. *Geochim. Cosmochim. Acta* 75, 1920-1932. doi: 10.1016/j.gca. 2011.01.010

Gagnon, A. C., Adkins, J. F., and Erez, J. (2012). Seawater transport during coral biomineralization. *Earth Planet. Sci. Lett.* 32, 150-161. doi: 10.1016/j.epsl.2012. 03.005

Gazeau, F., Parker, L. M., Comeau, S., Gattuso, J.-P. J.-P., O'Connor, W. A., Martin, S., et al. (2013). Impacts of ocean acidification on marine shelled molluscs. *Mar. Biol.* 160, 2207-2245. doi: 10.1007/s00227-013-2219-3

Gillikin, D. P., Dehairs, F., Lorrain, A., Steenmans, D., Baeyens, W., and André, L. (2006). Barium uptake into the shells of the common mussel (*Mytilus edulis*) and the potential for estuarine paleo-chemistry reconstruction. *Geochim. Cosmochim. Acta* 70, 395-407. doi: 10.1016/j.gca.2005.09.015

Giuffre, A. J., Gagnon, A. C., De Yoreo, J. J., and Dove, P. M. (2015). Isotopic tracer evidence for the amorphous calcium carbonate to calcite transformation by dissolution-reprecipitation. *Geochim. Cosmochim. Acta* 165, 407-417. doi: 10.1016/j.gca.2015.06.002

Guillermic, M., Misra, S., Eagle, R., Villa, A., Chang, F., and Tripati, A. (2020). Seawater pH reconstruction using boron isotopes in multiple planktonic foraminifera species with different depth habitats and their potential to constrain pH and pCO2 gradients. *Biogeosciences* 17, 3487-3510.

Harper, E. M., Palmer, T. J., and Alphey, J. R. (1997). Evolutionary response by bivalves to changing Phanerozoic sea-water chemistry. *Geol. Mag.* 134, 403-407. doi: 10.1017/s0016756897007061

Hastie, T. J., and Tibshirani, R. J. (1990). *Generalized Additive Models*, Vol. 43. Boca Raton, FL: CRC press.

Hastie, T., Tibshirani, R., Narasimhan, B., and Gilbert, C. (2020). *impute: Imputation for Microarray data*. R Package. doi: 10.18219/B9.bioc.impute

Havach, S. M., Chandler, G. T., Wilson-Finelli, A., and Shaw, T. J. (2001). Experimental determination of trace element partition coefficients in cultured benthic foraminifera. *Geochim. Cosmochim. Acta* 65, 1277-1283. doi: 10.1016/S0016-7037(00)00563-9

Haynes, L. L., Hönisch, B., Holland, K., Rosenthal, Y., and Eggins, S. M. (2019). Evaluating the planktic foraminiferal B/Ca proxy for application to deep time paleoceanography. *Earth Planet. Sci. Lett.* 528, 115824. doi: 10.1016/j.epsl.2019. 115824

Heinemann, F., Launspach, M., Gries, K., and Fritz, M. (2011). Gastropod nacre: structure, properties and growth - Biological, chemical and physical basics. *Biophys. Chem.* 153, 126-153. doi: 10.1016/j.bpc.2010.11.003

Hemming, N. G., and Hanson, G. N. (1992). Boron isotopic composition and concentration in modern marine carbonates. *Geochim. Cosmochim. Acta* 56, 537-543. doi: 10.1016/0016-7037(92)90151-8

Ho, L. S. T., Ane, C., Lachlan, R., Tarpinian, K., Feldman, R., Yu, Q., et al. (2016). *Package 'Phylolm'. Available online at: http://cran.r-project. org/web/packages/phylolm/index.html* (accessed February, 2018).

Holcomb, M., DeCarlo, T. M., Gaetani, G. A., and McCulloch, M. (2016). Factors affecting B/Ca ratios in synthetic aragonite. *Chem. Geol.* 437, 67-76. doi: 10. 1016/j.chemgeo.2016.05.007

Holtmann, W. C. (2013). *Sea Urchin Membrane Transport Mechanisms for Calcification and pH Homeostasis*. Kiel: Christian Albrechts University.

Hope, R. M. (2013). *Rmisc: Rmisc: Ryan Miscellaneous*. R package version 1.5. Available online at: https://CRAN.R-project.org/package=Rmisc (accessed October 22, 2013).

Iglikowska, A., Beldowski, J., Chelchowski, M., Chierici, M., Kędra, M., Przytarska, J., et al. (2016). Chemical composition of two mineralogically contrasting Arctic bivalves' shells and their relationships to environmental variables. *Mar. Pollut. Bull.* 114, 903-916. doi: 10.1016/j.marpolbul.2016.10.071

Immenhauser, A., Schöne, B. R., Hoffmann, R., and Niedermayr, A. (2016). Mollusc and brachiopod skeletal hard parts: intricate archives of their marine environment. *Sedimentology* 63, 1-59. doi: 10.1111/sed.12231

Ives, A. R. (2018). R²s for correlated data: phylogenetic models, LMMs, and GLMMs. *Syst. Biol.* 68, 234-251. doi: 10.1093/sysbio/syy060

Ives, A. R., and Li, D. (2018). rr2: an R package to calculate R² s for regression models. *J. Open Source Softw.* 3:1028.

Jacob, D. E., Wirth, R., Agbaje, O. B. A., Branson, O., and Eggins, S. M. (2017). Planktic foraminifera form their shells via metastable carbonate phases. *Nat. Commun.* 8:1265. doi: 10.1038/s41467-017-00955-0

Jurikova, H., Ippach, M., Liebetrau, V., Gutjahr, M., Krause, S., Büsse, S., et al. (2020). Incorporation of minor and trace elements into cultured brachiopods: implications for proxy application with new insights from a biomineralisation model. *Geochim. Cosmochim. Acta* 286, 418-440. doi: 10.1016/j.gca.2020.07. 026

Keul, N., Langer, G., de Nooijer, L. J., Nehrke, G., Reichart, G.-J., and Bijma, J. (2013). Incorporation of uranium in benthic foraminiferal calcite reflects seawater carbonate ion concentration. *Geochem. Geophys. Geosyst.* 14, 102-111. doi: 10.1029/2012GC004330

Kondo, H., Toyofuku, T., and Ikeya, N. (2005). Mg/Ca ratios in the shells of cultured specimens and natural populations of the marine ostracode Xestoleberis hanaii (Crustacea). *Palaeogeogr. Palaeoclimatol. Palaeoecol.* 225, 3-13. doi: 10.1016/j.palaeo.2004.05.026

Konrad, F., Purgstaller, B., Gallien, G., Mavromatis, V., Gane, P., and Dietzel, M. (2018). Influence of aqueous Mg concentration on the transformation of amorphous calcium carbonate. *J. Cryst. Growth* 498, 381-390. doi: 10.1016/j. jcrystgro.2018.07.018

Kunioka, D., Shirai, K., Takahata, N., Sano, Y., Toyofuku, T., and Ujiie, Y. (2006). Microdistribution of Mg/Ca, Sr/Ca, and Ba/Ca ratios in Pulleniatina obliquiloculata test by using a NanoSIMS: implication for the vital effect mechanism. *Geochem. Geophys. Geosyst.* 7:Q1220. doi: 10.1029/2006GC00 1280

Lavigne, M., Hill, T. M., Sanford, E., Gaylord, B., Russell, A. D., Lenz, E. A., et al. (2013). The elemental composition of purple sea urchin (*Strongylocentrotus purpuratus*) calcite and potential effects of pCO2 during early life stages. *Biogeosciences* 10, 3465-3477. doi: 10.5194/bg-10-3465-2013

Le, N., Zidek, J., White, R., Cubranic, D., Sampson, P. D., Guttorp, P., et al. (2015). *EnviroStat*: Statistical Analysis of Environmental Space-Time Processes. R package version 0.4-2. Available online at: https://CRAN.R-project.org/ package=EnviroStat (accessed June 03, 2015).

Lea, D. W. (2013). *Elemental and Isotopic Proxies of Past Ocean Temperatures. Treatise on Geochemistry*, 2nd Edn, Vol. 8. Oxford: Elsevier Ltd, doi: 10.1016/ B978-0-08-095975-7.00614-8

Letunic, I., and Bork, P. (2019). Interactive tree Of Life (iTOL) v4: recent updates and new development. *Nucleic Acids Res.* 47, 256-259.

Levi-Kalisman, Y., Raz, S., Weiner, S., Addadi, L., and Sagì, I. (2002). Structural differences between biogenic amorphous calcium carbonate phases using X-ray absorption spectroscopy. *Adv. Funct. Mater.* 12, 43-48. doi: 10.1002/1616- 3028(20020101)12:1<43::AID-ADFM43<3.0.CO;2-C

Liu, Y.-W. W., Sutton, J. N., Ries, J. B., and Eagle, R. A. (2020). Regulation of calcification site pH is a polyphyletic but not always governing response to ocean acidification. *Sci. Adv.* 6, 7-8. doi: 10.1126/sciadv.aax1314

Long, X., Ma, Y., and Qi, L. (2014). Biogenic and synthetic high magnesium calcite – A review. *Journal of Structural Biology* 185, 1-14. doi: 10.1016/j.jsb. 2013.11.004

Lorens, R. B., and Bender, M. L. (1977). Physiological exclusion of magnesium from *Mytilus edulis* calcite. *Nature* 269, 793-794.

Lorens, R. B., and Bender, M. L. (1980). The impact of solution chemistry on *Mytilus edulis* calcite and aragonite. *Geochim. Cosmochim. Acta* 44, 1265-1278. doi: 10.1016/0016-7037

Loste, E., Wilson, R. M., Seshadri, R., and Meldrum, F. C. (2003). The role of magnesium in stabilising amorphous calcium carbonate and controlling calcite morphologies. *J. Cryst. Growth* 254, 206-218. doi: 10.1016/S0022-0248(03) 01153-9

Lowenstam, H. A. (1981). Minerals formed by weathering. *Science* 211, 1126-1131.

Luquet, G. (2012). Biomineralizations: insights and prospects from crustaceans. *ZooKeys* 176, 103-121. doi: 10.3897/zookeys.176.2318

Marchitto, T. M., Oppo, D. W., and Curry, W. B. (2002). Paired benthic foraminiferal Cd/Ca and Zn/Ca evidence for a greatly increased presence of Southern Ocean Water in the glacial North Atlantic. *Paleoceanography* 17:1038. doi: 10.1029/2000PA000598

Märkel, K., Röser, U., Mackenstedt, U., and Klostermann, M. (1986). Ultrastructural investigation of matrix-mediated biomineralization in echinoids (Echinodermata, Echinoida). *Zoomorphology* 106, 232-243. doi: 10.1007/BF00312044

Marra, G., and Wood, S. N. (2011). Practical variable selection for generalized additive models. *Comput. Stat. Data Anal.* 55, 2372-2387. doi: 10.1016/j.csda. 2011.02.004
Mason, A. Z., and Nott, J. A. (1981). The role of intracellular biomineralized granules in the regulation and detoxification of metals in gastropods with special reference to the marine prosobranch *Littorina littorea*. *Aquat. Toxicol.* 1, 239–256. doi: 10.1016/0166-445X(81)90018-7

Mass, T., Giuffre, A. J., Sun, C.-Y. Y., Stifler, C. A., Frazier, M. J., Neder, M., et al. (2017). Amorphous calcium carbonate particles form coral skeletons. *Proc. Natl. Acad. Sci. U.S.A.* 114, E7670–E7678. doi: 10.1073/pnas.1707890114

Mavromatis, V., Goetschl, K. E., Grengg, C., Konrad, F., Purgstaller, B., and Dietzel, M. (2018). Barium partitioning in calcite and aragonite as a function of growth rate. *Geochim. Cosmochim. Acta* 237, 65–78. doi: 10.1016/j.gca.2018.06.018

Mavromatis, V., Montouillout, V., Noireaux, J., Gaillardet, J., and Schott, J. (2015). Characterization of boron incorporation and speciation in calcite and aragonite from co-precipitation experiments under controlled pH, temperature and precipitation rate. *Geochim. Cosmochim. Acta* 150, 299–313. doi: 10.1016/j.gca.2014.10.024

McCulloch, M. T., D'Olive, J. P., Falter, J., Holcomb, M., and Trotter, J. A. (2017). Coral calcification in a changing World and the interactive dynamics of pH and DIC upregulation. *Nat. Commun.* 8, 1–8. doi: 10.1038/ncomms15686

Menadakis, M., Maroulis, G., and Koutskos, P. G. (2009). Incorporation of Mg²⁺, Sr²⁺, Ba²⁺ and Zn²⁺ into aragonite and comparison with calcite. *J. Math. Chem.* 46, 484–491. doi: 10.1007/s10910-008-9490-4

Mergelsberg, S. T., De Yoreo, J. J., Miller, Q. R. S., Marc Michel, F., Ulrich, R. N., and Dove, P. M. (2020). Metastable solubility and local structure of amorphous calcium carbonate (ACC). *Geochim. Cosmochim. Acta* 289, 196–206. doi: 10.1016/j.gca.2020.06.030

Mergelsberg, S. T., Ulrich, R. N., Xiao, S., and Dove, P. M. (2019). Composition systematics in the exoskeleton of the american lobster, *Homarus americanus* and implications for malacostraca. *Front. Earth Sci.* 7:69. doi: 10.3389/feart. 2019.00069

Meseck, S. L., Mercaldo-Allen, R., Kuropat, C., Clark, P., and Goldberg, R. (2018). Variability in sediment-water carbonate chemistry and bivalve abundance after bivalve settlement in Long Island Sound, Milford, Connecticut. *Mar. Pollut. Bull.* 135, 165–175. doi: 10.1016/j.marpolbul.2018.07.025

Montagna, P., McCulloch, M., Douville, E., López Correa, M., Trotter, J., Rodolfo-Metalpa, R., et al. (2014). Li/Mg systematics in scleractinian corals: calibration of the thermometer. *Geochim. Cosmochim. Acta* 132, 288–310. doi: 10.1016/j.gca.2014.02.005

Mucci, A., and Morse, J. W. (1983). The incorporation of Mg²⁺ and Sr²⁺ into calcite overgrowths: influences of growth rate and solution composition. *Geochim. Cosmochim. Acta* 47, 217–233. doi: 10.1016/0016-7037(83)90135-7

Nassif, N., Pinna, N., Gehrke, N., Antonietti, M., Jäger, C., and Cölfen, H. (2005). Amorphous layer around aragonite platelets in nacre. *Proc. Natl. Acad. Sci. U.S.A.* 102, 12653–12655. doi: 10.1073/pnas.0502577102

Nguyen, T. T. (2013). *Li/Ca, B/Ca, and Mg/Ca Composition of Cultured Sea Urchin Spines and Paleo-Echinoderms Measured Using a Secondary Ion Mass Spectrometer*. Los Angeles, CA: University of California.

Ogilvie, M. (1896). Corals: microscopic and systematic study of Madreporian types of corals. *Phil. Trans. R Soc. Lond.* 187, 83–345.

Oksanen, J., Kindt, R., Legendre, P., O'Hara, B., Simpson, G. L., Solymos, P. M., et al. (2008). *The Vegan Package. Community Ecology Package*, (January). 190. Available online at: https://bcrc.bio.umass.edu/biometry/images/8/85/Vegan. pdf (accessed January, 2008).

Pagel, M. (1999). Inferring the historical patterns of biological evolution. *Nature* 401, 877–884.

Park, W. K., Ko, S. J., Lee, S. W., Cho, K. H., Ahn, J. W., and Han, C. (2008). Effects of magnesium chloride and organic additives on the synthesis of aragonite precipitated calcium carbonate. *J. Cryst. Growth* 310, 2593–2601. doi: 10.1016/j.jcrysgro.2008.01.023

Piwoni-Piórewicz, A., Strekopytov, S., Humphreys-Williams, E., and Kukliński, P. (2021). The patterns of elemental concentration (Ca, Na, Sr, Mg, Mn, Ba, Cu, Pb, V, Y, U and Cd) in shells of invertebrates representing different CaCO₃ polymorphs: a case study from the brackish Gulf of Gdańsk (the Baltic Sea). *Biogeosciences* 18, 707–728. doi: 10.5194/bg-18-707-2021

Purgstaller, B., Goetschl, K. E., Mavromatis, V., and Dietzel, M. (2019). Solubility investigations in the amorphous calcium magnesium carbonate system. *Cryst. Eng. Commun.* 21, 155–164. doi: 10.1039/c8ce01596a

Purgstaller, B., Mavromatis, V., Immenhauser, A., and Dietzel, M. (2016). Transformation of Mg-bearing amorphous calcium carbonate to Mg-calcite – In situ monitoring. *Geochim. Cosmochim. Acta* 174, 180–195. doi: 10.1016/j.gca. 2015.10.030

R Core Team (2019). *R: A Language and Environment for Statistical Computing*. Vienna: R Foundation for Statistical Computing.

Rahman, M. A., Halfar, J., Adey, W. H., Nash, M., Paulo, C., and Dittrich, M. (2019). The role of chitin-rich skeletal organic matrix on the crystallization of calcium carbonate in the crustose coralline alga *Leptophytum foecundum*. *Sci. Rep.* 9:11869. doi: 10.1038/s41598-019-47785-2

Raitzsch, M., Dueñas-Bohórquez, A., Reichart, G.-J., de Nooijer, L. J., and Bickert, T. (2010). Incorporation of Mg and Sr in calcite of cultured benthic foraminifera: impact of calcium concentration and associated calcite saturation state. *Biogeosciences* 7, 869–881. doi: 10.5194/bg-7-869- 2010

Raz, S., Weiner, S., and Addadi, L. (2000). Formation of high-magnesian calcites via an amorphous precursor phase: possible biological implications. *Adv. Mater.* 12, 38–42. doi: 10.1002/(SICI)1521-4095(200001)12:1<38::AID-ADMA38<3.0.CO;2-1

Reeder, R. J. (1983). Crystal chemistry of the rhombohedral carbonates. *Rev. Mineral.* 11, 1–47.

Revell, L. J. (2009). Size-correction and principal components for interspecific comparative studies. *Evolution* 63, 3258–3268.

Revell, L. J. (2012). phytools: an R package for phylogenetic comparative biology (and other things). *Methods Ecol. Evol.* 3, 217–223.

Ries, J. B. (2006). Mg fractionation in crustose coralline algae: geochemical, biological, and sedimentological implications of secular variation in the Mg/Ca ratio of seawater. *Geochim. Cosmochim. Acta* 70, 891–900. doi: 10.1016/j.gca. 2005.10.025

Ries, J. B. (2009). Effects of secular variation in seawater Mg/Ca ratio (calcite- aragonite seas) on CaCO₃ sediment production by the calcareous algae Halimeda, Penicillus and Udotea – Evidence from recent experiments and the geological record. *Terra Nova* 21, 323–339. doi: 10.1111/j.1365-3121.2009. 00899.x

Ries, J. B. (2010). Review: geological and experimental evidence for secular variation in seawater Mg/Ca (calcite-aragonite seas) and its effects on marine biological calcification. *Biogeosciences* 7, 2795–2849. doi: 10.5194/bg-7-2795- 2010

Ries, J. B. (2011). Skeletal mineralogy in a high-CO₂ world. *J. Exp. Mar. Biol. Ecol.* 403, 54–64. doi: 10.1016/j.jembe.2011.04.006

Ries, J. B., Cohen, A. L., and McCorkle, D. C. (2009). Marine calcifiers exhibit mixed responses to CO₂-induced ocean acidification. *Geology* 37, 1131–1134. doi: 10.1130/G30210A.1

Ries, J. B., Ghazaleh, M. N., Connolly, B., Westfield, I., and Castillo, K. D. (2016). Impacts of seawater saturation state (ΩA=0.4–4.6) and temperature (10, 25 °C) on the dissolution kinetics of whole-shell biogenic carbonates. *Geochim. Cosmochim. Acta* 192, 318–337. doi: 10.1016/j.gca.2016. 07.001

Roer, R., and Dillaman, R. (1984). The structure and calcification of the crustacean cuticle. *Integr. Comp. Biol.* 24, 893–909. doi: 10.1093/icb/24.4.893

Rosenheim, B. E., Swart, P. K., and Thorrold, S. R. (2005). Minor and trace elements in sclerosponge *Ceratoporella nicholsoni*: biogenic aragonite near the inorganic endmember? *Palaeogeogr. Palaeoclimatol. Palaeoecol.* 228, 109–129. doi: 10.1016/j.palaeo.2005.03.055

Rosenthal, Y., Boyle, E. A., and Slowey, N. (1997). Temperature control on the incorporation of magnesium, strontium, fluorine, and cadmium into benthic foraminiferal shells from Little Bahama Bank: prospects for thermocline paleoceanography. *Geochim. Cosmochim. Acta* 61, 3633–3643.

Russell, A. D., Hönisch, B., Spero, H. J., and Lea, D. W. (2004.). Effects of seawater carbonate ion concentration and temperature on shell U, Mg, and Sr in cultured planktonic foraminifera. *Geochim. Cosmochim. Acta* 68, 4347–4361. doi: 10.1016/j.gca.2004.03.013

Schöne, B. R., Zhang, Z., Jacob, D., Gillikin, D. P., Tütken, T., Garbe-Schönberg, D., et al. (2010). Effect of organic matrices on the determination of the trace element chemistry (Mg, Sr, Mg/Ca, Sr/Ca) of aragonitic bivalve shells (Arctica islandica) – comparison of ICP-OES and LA-ICP-MS data. *Geochem. J.* 44, 23–37. doi: 10.2343/geochemj.1.0045
Shen, G. T., Campbell, T. M., Dunbar, R. B., Wellington, G. M., Colgan, M. W., and Glynn, P. W. (1991). Paleochemistry of manganese in corals from the Galapagos Islands. *Coral Reefs* 10, 91–100.

Short, J., Foster, T., Falter, J., Kendrick, G. A., and McCulloch, M. T. (2015). Crustose coralline algal growth, calcification and mortality following a marine heatwave in Western Australia. *Cont. Shelf Res.* 106, 38–44. doi: 10.1016/j.csr.2015.07.003

Simkiss, K., Taylor, M., and Greaves, G. N. (1986). Amorphous structure of intracellular mineral granules. *Biochem. Soc. Trans.* 14, 549–552. doi: 10.1042/bst0140549

Sinclair, D. J., and Risk, M. J. (2006). A numerical model of trace-element coprecipitation in a physicochemical calcification system: application to coral biomineralization and trace-element "vital effects". *Geochim. Cosmochim. Acta* 70, 3855–3868. doi: 10.1016/j.gca.2006.05.019

Stanley, S. M., and Hardie, L. A. (1999). Hypercalcification: paleontology links plate tectonics. *GSA Today* 9, 1–7.

Stumpp, M., Hu, M. Y., Melzner, F., Gutowska, M. A., Dorey, N., Himmerkus, N., et al. (2012). Acidified seawater impacts sea urchin larvae pH regulatory systems relevant for calcification. *Proc. Natl. Acad. Sci. U.S.A.* 109, 18192–18197. doi: 10.1073/pnas.1209174109

Suchard, M. A., Lemey, P., Baele, G., Ayres, D. L., Drummond, A. J., and Rambaut, A. (2018). Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. *Virus Evol.* 4:vey016.

Sugawara, A., and Kato, T. (2000). Aragonite CaCO₃ thin-film formation by cooperation of Mg2+ and organic polymer matrices. *Chem. Commun.* 6, 487–488. doi: 10.1039/A909566G

Sutton, J. N., Liu, Y. W., Ries, J. B., Guillermic, M., Ponzevera, E., and Eagle, R. A. (2018). δ11B as monitor of calcification site pH in divergent marine calcifying organisms. *Biogeosciences* 15, 1447–1467. doi: 10.5194/bg-15-1447-2018

Tambutté, E., Tambutté, S., Segonds, N., Zoccola, D., Venn, A., Erez, J., et al. (2011). Calcein labelling and electrophysiology: insights on coral tissue permeability and calcification. *Proc. R. Soc. B Biol. Sci.* 279, 19–27. doi: 10.1098/rspb.2011.0733

Tripati, A. K., Roberts, C. D., and Eagle, R. A. (2009). Coupling of CO₂ and Ice sheet stability over major climate transitions of the last 20 million years. *Science* 326, 1394–1397. doi: 10.1126/science.1178296

Tripati, A. K., Roberts, C. D., Eagle, R. A., and Li, G. (2011). A 20 million year record of planktic foraminiferal B/Ca ratios: systematics and uncertainties in pCO₂ reconstructions. *Geochim. Cosmochim. Acta* 75, 2582–2610. doi: 10.1016/j.gca.2011.01.018

Venn, A. A., Bernardet, C., Chabenat, A., Tambutté, E., and Tambutté, S. (2020). Paracellular transport to the coral calcifying medium: effects of environmental parameters. *J. Exp. Bio.* 223:jeb227074. doi: 10.1242/jeb.227074

Von Euw, S., Zhang, Q., Manichev, V., Murali, N., Gross, J., Feldman, L. C., et al. (2017). Biological control of aragonite formation in stony corals. *Science* 356, 933–938. doi: 10.1126/science.aam6371

Watson, E. B. (2004). A conceptual model for near-surface kinetic controls on the trace-element and stable isotope composition of abiogenic calcite crystals. *Geochim. Cosmochim. Acta* 68, 1473–1488. doi: 10.1016/j.gca.2003.10.003

Weiner, S., and Dove, P. M. (2003). An overview of biomineralization processes and the problem of the vital effect. *Rev. Mineral. Geochem.* 54, 1–29. doi: 10.2113/0540001

Weiss, I. M., Tuross, N., Addadi, L., and Weiner, S. (2002). Mollusc larval shell formation: amorphous calcium carbonate is a precursor phase for aragonite. *J. Exp. Zool.* 293, 478–491. doi: 10.1002/jez.90004

Wickham, H. (2009). *ggplot2: Elegant Graphics for Data Analysis – Hadley Wickham – Google Books. Springer Science & Business Media*. Houston, TX: Springer Nature. doi: 10.1007/978-3-319-24277-4

Wickham, H. (2016). *ggplot2: Elegant Graphics for Data Analysis*. New York, NY: Springer-Verlag.

Wickins, J. F. (1984). The effect of reduced pH on carapace calcium, strontium and magnesium levels in rapidly growing prawns (*Penaeus monodon fabricius*). *Aquaculture* 41, 49–60. doi: 10.1016/0044-8486(84)90389-2

Wilbur, K. M., Colinvaux, L. H., and Watabe, N. (1969). Electron microscope study of calcification in the alga Halimeda (order Siphonales). *Phycologia* 8, 27–35. doi: 10.2216/i0031-8884-8-1-27.1

Wilke, C. O. (2020). *ggtext: Improved Text Rendering Support for 'ggplot2'. R package version 0.1.0*. Available online at: https://CRAN.R-project.org/package=ggtext (accessed Dec 17, 2020).

Wilt, F. H. (2002). Biomineralization of the spicules of sea urchin embryos. *Zool. Sci.* 19, 253–261. doi: 10.2108/zzj.19.253

Wood, S.N. (2011). Fast stable restricted maximum likelihood and marginal likelihood estimation of semiparametric generalized linear models. *J. R. Stat. Soc. B* 73, 3–36.

Wood, S. N. (2017). *Generalized Additive Models: An Introduction with R*, 2nd Edn. Boca Raton, FL: Chapman and Hall, doi: 10.1201/9781315370279

Yu, J., Day, J., Greaves, M., and Elderfield, H. (2005). Determination of multiple element/calcium ratios in foraminiferal calcite by quadrupole ICP-MS. *Geochem. Geophys. Geosys.* 6:Q0801. doi: 10.1029/2005GC000964

Yu, J., Elderfield, H., and Hönisch, B. (2007). B/Ca in planktonic foraminifera as a proxy for surface seawater pH. *Paleoceanography* 22:A2202. doi: 10.1029/2006PA001347

Yu, J., Foster, G. L., Elderfield, H., Broecker, W. S., and Clark, E. (2010). An evaluation of benthic foraminiferal B/Ca and δ11B for deep ocean carbonate ion and pH reconstructions. *Earth Planet. Sci. Lett.* 293, 114–120. doi: 10.1016/j.epsl.2010.02.029

Zacherl, D. C., Paradis, G., and Lea, D. W. (2003). Barium and strontium uptake into larval protoconchs and statoliths of the marine neogastropod Kelletia kelledi. *Geochim. Cosmochim. Acta* 67, 4091–4099.

Zeileis, A., and Hothorn, T. (2002). Diagnostic checking in regression relationships. *R News* 2, 7–10.

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2021 Ulrich, Guillermic, Campbell, Hakim, Han, Singh, Stewart, Román-Palacios, Carroll, De Corte, Gilmore, Doss, Tripati, Ries and Eagle. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.