Abstract
Endocrine disruption is suspected in cnidarians, but questions remain how occurs. Steroid sex hormones are detected in corals and sea anemones even though these animals do not have estrogen receptors and their repertoire of steroidogenic enzymes appears to be incomplete. Pathways associated with sex hormone biosynthesis and sterol signaling are an understudied area in cnidarian biology. The objective of this study was to identify a suite of genes that can be linked to exposure of endocrine disruptors. Exaiptasia diaphana were exposed to nominal 20ppb concentrations of estradiol (E2), testosterone (T), cholesterol, oxybenzone (BP-3), or benzyl butyl phthalate (BBP) for 4 h. Eleven genes of interest (GOIs) were chosen from a previously generated EST library. The GOIs are 17β-hydroxysteroid dehydrogenases type 14 (17β HSD14) and type 12 (17β HSD12), Niemann-Pick C type 2 (NPC2), Equistatin (EI), Complement component C3 (C3), Cathepsin L (CTSL), Patched domain-containing protein 3 (PTCH3), Smoothened (SMO), Desert Hedgehog (DHH), Zinc finger protein GLI2 (GLI2), and Vitellogenin (VTG). These GOIs were selected because of functional associations with steroid hormone biosynthesis; cholesterol binding/transport; immunity; phagocytosis; or Hedgehog signaling. Quantitative Real-Time PCR quantified expression of GOIs. In silico modelling utilized protein structures from Protein Data Bank as well as creating protein structures with SWISS-MODEL. Results show transcription of steroidogenic enzymes, and cholesterol binding/transport proteins have similar transcription profiles for E2, T, and cholesterol treatments, but different profiles when BP-3 or BBP is present. C3 expression can differentiate between exposures to BP-3 versus BBP as well as exposure to cholesterol versus sex hormones. In silico modelling revealed all ligands (E2, T, cholesterol, BBP, and BP-3) have favorable binding affinities with 17β HSD14, 17β HSD12, NPC2, SMO, and PTCH proteins. VTG expression was down-regulated in the sterol treatments but up-regulated in BP-3 and BBP treatments. In summary, these eleven GOIs collectively generate unique transcriptional profiles capable of discriminating between the five chemical exposures used in this investigation. This suite of GOIs are candidate biomarkers for detecting transcriptional changes in steroidogenesis, gametogenesis, sterol transport, and Hedgehog signaling. Detection of disruptions in these pathways offers new insight into endocrine disruption in cnidarians.
Introduction
In the past quarter century, public awareness of the dramatic declines in coral reefs has become increasingly evident. Extensive research has characterized how coral reefs are shifting in structure and biodiversity due to climate change, over-fishing, coastal development, and pollution (Gardner et al., 2003; Hughes et al., 2003; Pandolfi et al., 2003; Jones et al., 2004; Bruno and Selig, 2007). Understanding the impact of land-based pollution onto reefs is further complicated by pulses of anthropogenic activity, variations in tides, variations in seasonal precipitation, and sediments which can act as carriers for a variety of organic compounds (Dachs et al., 1999; Gavio et al., 2010; Singh et al., 2010; Burns and Brinkman, 2011; Edge et al., 2013). Both land-based pollution and sewage discharge produce terrestrial runoff which culminates in diminished water quality in coastal marine environments (McKenna et al., 2001; Fabricius, 2005; Gavio et al., 2010; Singh et al., 2010; Vidal-Dorsch et al., 2012; Bahr et al., 2015; French et al., 2015; Wear and Thurber, 2015; Yoshioka et al., 2016; Al-Jandal et al., 2018). The composition of sewage is a complex mixture of freshwater, heavy metals, inorganic nutrients, microplastics, pathogens, pesticides, pharmaceuticals, plasticizers, sediments, suspended solids, and toxins (Jones et al., 2003; Kerswell and Jones, 2003; Flood et al., 2005; Jones, 2010; Singh et al., 2010; Sutherland et al., 2010; Tijani et al., 2013; French et al., 2015; Mahon et al., 2017). Many components of sewage are recognized as endocrine disrupting chemicals (EDCs) whose impact on marine invertebrates needs to be characterized (Tan et al., 2007; Singh et al., 2010; Wear and Thurber, 2015).
One form of pollution in the marine environment that is ubiquitous even in sewage is plastic (Law, 2017; Mahon et al., 2017; Lebreton et al., 2018). Plastics are known to accumulate and persist longer in sediments than on land (Worm et al., 2017). With the global explosion of plastic pollution, there is evidence cnidarians are capable of ingesting microplastics which are known to release phthalates (Hall et al., 2015; Axworthy and Padilla-Gamiño, 2019; Rotjan et al., 2019; Deng et al., 2020). As a leachable compound from plastics, phthalates are one of the most frequently detected persistent organic pollutants in the environment (Gao and Wen, 2016). In addition to plastics, oxybenzone is becoming an increasingly significant environmental concern for coral reefs because it has been detected in coastal surface waters and sediments, and has been shown to influence larval development (Downs et al., 2016; Mitchelmore et al., 2019). Steroid sex hormones, phthalates, and oxybenzone are all recognized as EDCs and all have been detected in sewage effluent and coastal marine environments where anthropogenic activity is prevalent (Atkinson et al., 2003; Thibaut and Porte, 2004; Singh et al., 2010; Bargar et al., 2013; Tsui et al., 2014; French et al., 2015; Thorn et al., 2015; Al-Jandal et al., 2018; LaPlante et al., 2018; Matouskova et al., 2020). Sterols represent the second highest proportion of lipids (14–17%) found in cnidarians (Revel et al., 2016). Because steroid hormones, phthalates, and oxybenzone all have lipophilic characteristics, the lipid-rich tissues of cnidarians are likely targets for EDCs. It has been proposed that lipophilic compounds will easily diffuse into cnidarian tissues and bioaccumulate (Peters et al., 1997; Imbs, 2013; Ko et al., 2014). Using a variety of mass-spectrometry assays, steroid hormones (and/or their conjugates), phthalates, and oxybenzone have been detected in coral tissues (from <10 ng up to 650 ng g−1 dry weight) and/or the environment (from <10 ng L−1 to >100 μg L−1), (Solbakken et al., 1985; Tarrant et al., 2003; Twan et al., 2003; Twan et al., 2006; Downs et al., 2016; Mitchelmore et al., 2019). Cnidarians have demonstrated the ability to uptake estrogens in the water column at concentrations as low as 300 pg/L (Tarrant et al., 2001).
Steroid sex hormones are known to be important in cnidarian development and reproduction (Atkinson and Atkinson, 1992; Gassman and Kennedy, 1992; Slattery et al., 1999; Tarrant et al., 1999; Twan et al., 2003; Twan et al., 2006; Armoza-Zvuloni et al., 2014). Exogenous estrogen is known to impact asexual reproduction (Thorn et al., 2015). A conundrum exists because anthozoans do not have nuclear estrogen receptors even though the hormone is detectable in tissues and functions in their reproductive biology (Reitzel et al., 2008; Reitzel and Tarrant, 2009). Nuclear receptors capable of binding estrogen have been identified in other cnidarians (medusozoans) (Khalturin et al., 2018), and although not yet well-characterized, a membrane-bound G-protein coupled estrogen receptor 1-like transcript (GPER1) genomic locus (LOC114576065) has been identified for E. diaphana (NCBI Nucleotide). GPER1 has recently been identified in humans in addition to the two traditional nuclear estrogen receptors ERα and ERβ (Fuentes and Silveyra, 2019). Perhaps GPER1 is also involved in physiologic estrogen-signaling in anthozoans and may be a target of EDCs, leading to downstream transcriptional changes.
Another conundrum exists because cnidarians are known to synthesize a diverse set of sterols (Tarrant et al., 2009; Revel et al., 2016), but their repertoire of steroidogenic enzymes appears to be incomplete. There is no genomic evidence that cnidarians produce aromatase, and yet aromatase activity has been detected (Tarrant et al., 2003; Twan et al., 2003; Twan et al., 2006). It is also possible that steroid sex hormones may play other significant roles (e.g., possibly chemical communication) which have yet to be fully elucidated (Tarrant et al., 2003; Blomquist et al., 2006; Lafont and Mathieu, 2007; Shikina et al., 2016). Since cnidarians use steroid hormones as signaling molecules but they do not have an endocrine system, Tarrant (2007) asked a fundamental question: Can endocrine disruption occur in cnidarians? Previous investigations have identified steroid hormone-like signaling molecules but there are still significant gaps in understanding which metabolic pathways are impacted.
One approach of ecotoxicological investigations is to identify changes in gene transcription in aquatic organisms exposed to an exogenous substance (Snell et al., 2003; Snape et al., 2004). Calls have been made to develop a comprehensive assessment of steroid hormones, phthalates, and personal care products impacting marine environments (Singh et al., 2010; Tijani et al., 2013; Archer et al., 2017; DiNardo and Downs 2018). More detailed information regarding the impact of EDCs on specific pathways involved in both synthesis and metabolism of sex hormones would be highly desirable. Developing biomarkers for anthropogenic stressors can be useful tools for monitoring health on coral reefs (Kenkel et al., 2014). Validated mechanistic biomarkers for aquatic invertebrates which can identify responses to a particular class of chemical stressor a have been theorized (Snell et al., 2003; Hutchinson, 2007), but difficult to produce. Presently, biomarkers for cnidarians exposed to EDCs have not been developed and validated. This study seeks to identify a suite of genes that characterize exogenous EDC exposures. Such a suite of genes could help elucidate the cellular pathways leading to cnidarian endocrine disruption.
Materials and Methods
Toxicant Exposure
Sea anemones (Exaiptasia diaphana previously known as Aiptasia pallida) were purchased from a supplier (Wards Natural Science, Rochester, NY, United States) and acclimated to laboratory conditions (45 L recirculating artificial seawater at 22°C and 36 ppt salinity) for 1 week. Anemones (100 individuals) were randomly subdivided into a control group or 20 ppb treatments of 17β-estradiol (E2), testosterone [T], cholesterol, benzyl butyl phthalate (BBP), or oxybenzone (BP-3). Toxicant purity was pharmaceutical secondary standard grade for E2, T, cholesterol, and BP-3 while BBP was analytical standard grade (Sigma-Aldrich, St. Louis, MO, United States). Acetone (150 μg/L) was used to solubilized toxicants were prior to dilution in seawater. Acetone has previously been used as a wetting agent (Morgan et al., 2001; Morgan et al., 2012). The 20 ppb concentration was chosen to induce responses and does not necessarily reflect environmental conditions even though BP-3 has been detected at concentrations >20 ppb (Downs et al., 2016) and cholesterol at 11ppb (French et al., 2015). All treatments were nominal concentrations for 4 h in 1L artificial seawater under ambient laboratory lighting during the early spring season. The 4-h exposure was chosen to detect responses at an early stage of exposure.
Selecting Candidate Genes of Interest
Complementary DNA sequences previously generated by Representational Difference Analysis (RDA) (Edman et al., 1997; Hubank and Schatz, 1999; Pastorian et al., 2000) were screened for candidate ESTs to be used in this investigation. Candidate ESTs were selected based on sequence homology (BLASTX analysis) to genes of interest (GOI) with known functions involving sex hormone synthesis, sterol binding, lipid transport, and/or indirect sterol interactions. Subsets of these ESTs have previously been used to characterize cnidarian stress responses (Morgan et al., 2012; Morgan et al., 2015).
Reverse Transcription Reactions
Anemones from each treatment were pooled and total RNA was isolated using Trizol (Invitrogen, United States). Two milliliter phase-lock gels (5′Prime, United States) aided in the recovery of the aqueous phase. Total RNA was DNase I digested and purified (New England BioLabs, United States). Messenger RNA (mRNA) was isolated (Oligotex, Qiagen, United States) from 100 μg of DNase I treated total RNA. First stand synthesis used SuperScript IV (Invitrogen, United States) along with random hexamers and oligo-dT primers to reverse transcribe 1 μg of poly-A enriched mRNA. The RT reaction conditions were 1 h at 37°C, followed by 1 min at each temperature between 42 and 50°C.
Quantitative Real-Time PCR
A QuantStudio 7 Flex Real-Time PCR system (Applied Biosystems Waltham, MA, United States) used a SYBR Green-based assay to perform qPCR. A 1/100 dilution of first-strand synthesis reactions were used as templates for all qPCR reactions. Primers for each GOI were created using Primer3 (https://primer3.org) (Table 1). Components for each 20 µl qPCR reaction included: 10 µl Luna® Universal qPCR Mix (New England BioLabs, United States), 2.5 µl Forward primer (10 µM), 2.5 µl Reverse primer (10 µM), 2.5 µl dH2O, 2.5 µl sample. Thermocycling conditions were 1 cycle at 95°C for 1 min; 40 cycles of 95°C for 15 s and then 60°C for 30 s; and concluding with 1 cycle of melt curve analysis. Four replicate reactions were used for analyzing the relative expression of each GOI. Ribosomal protein L11 (RPL11) was used as qPCR reference gene (Kenkel et al., 2011). Melt-curve analysis, primer efficiencies, and gel electrophoresis confirmed specificity of priming. Replicate Cq values were averaged to determine ∆Cq and ∆∆Cq for each treatment and GOI. All ∆Cq and ∆∆Cq values are based on the consistent expression of RPL11 across all treatments. The ∆∆Cq method was used to determine the differences between targeted GOIs and a single reference gene (Bustin et al., 2009). Univariate ANOVA was performed on ∆∆Cq data to identify significant expression of individual GOIs across treatments. Similarities in variance between treatments were determined by Levene’s Test of Equality of Error Variances. If the variance between treatments was similar, then the Student-Neuman-Keuls (SNK) posthoc test was performed to determine which treatment(s) were significantly different from the rest. If variance between treatments was different, then Tamhane’s T2 posthoc test was applied since Univariate ANOVA is generally insensitive to heteroscedasticity.
TABLE 1
| Gene ID | Putative Gene Homolog | BLASTX | Organism ID | Homolog Accession # |
|---|---|---|---|---|
| E-value | ||||
| 17βHSD14 | 17β-hydroxysteroid dehydrogenase type 14 | 5e-155 | Exaiptasia diaphana | KXJ20962 |
| 17βHSD12 | 17β-hydroxysteroid dehydrogenase type 12 | 1e-61 | Exaiptasia diaphana | KXJ12187 |
| C3 | Complement component C3 | 3e-57 | Exaiptasia diaphana | KXJ11955 |
| CTSL | Cathepsin L | 6e-95 | Exaiptasia diaphana | KXJ27439 |
| NPC2 | Niemann-Pick C type 2 | 1e-84 | Exaiptasia diaphana | KXJ29862 |
| EI | Equistatin | 7e-57 | Exaiptasia diaphana | KXJ18222 |
| PTCH3 | Patched domain-containing protein 3 | 3e-79 | Exaiptasia diaphana | KXJ20037 |
| DHH | Desert hedgehog protein A | 3e-47 | Exaiptasia diaphana | KXJ11164 |
| SMO | Smoothened-like | 1e-87 | Exaiptasia diaphana | KXJ11374 |
| GLI2 | zinc finger protein GLi2 | 4e-146 | Exaiptasia diaphana | XP_020898684 |
| VTG | Vitellogenin 2 | 7e-38 | Exaiptasia diaphana | KXJ14544 |
BLAST results for RDA probes. Searches performed at NCBI using BLASTX algorithm and the non-redundant database (nr) with default search parameters.
In silico Molecular Modelling
Human protein homologs have previously been used as substitutes for cnidarian proteins in molecular modelling (Khalturin et al., 2018). Some cnidarian proteins have such significant homology to human proteins that they are even capable of stimulating human responses (Dunn et al., 2006; Wood-Charlson and Weis, 2009; Duffy and Frank, 2011; Dani et al., 2014; Brennan et al., 2017; Mansfield et al., 2017). The most similar vertebrate structures (according to BLASTX E-value) were used for selecting which Protein Databank (PDB) files should be used for initial docking experiments. Sequence alignments of cnidarian and relevant homologs determined which relevant protein sequences (i.e., PDB files) were to be used as representatives of corresponding cnidarian proteins. Protein coordinates (PDB files) for homologous constructs were downloaded from the PDB, and any existing ligands or small molecules (e.g., water or salt ions) were removed. Cnidarian protein structures were created with the SWISS-MODEL web server (https://swissmodel.expasy.org) using the homologs as templates. Ligands were created in Avogadro (Hanwell et al., 2012) and energy-minimized prior to docking. Docking was performed using AutoDock Vina (Trott and Olson, 2010); the search space was restricted to the known ligand binding site. After docking, the results were analyzed in UCSF Chimera (Pettersen et al., 2004) with the Find Clashes/Contacts function which was set to detect any atoms within 0.4 Å of the ligand with the best docking score.
Results
Genes of Interest
BLASTX analysis of ESTs identified eleven candidate GOIs with of functions associated with steroidogenesis, gametogenesis, cholesterol transport, immunity, phagocytosis, or the Hedgehog signaling pathway (Table 2). These GOIs were: 17β-hydroxysteroid dehydrogenase type 14 (17βHSD14), 17β-hydroxysteroid dehydrogenase type 12 (17βHSD12), Niemann-Pick C type 2 (NPC2), Equistatin (EI), Complement component C3 (C3), Cathepsin L (CTSL), Patched domain-containing protein 3 (PTCH3), Smoothened (SMO), Desert Hedgehog protein A (DHH), zinc finger protein Gli2 (GLI2), and Vitellogenin (VTG). Two GOIs (17βHSD14; 17βHSD12) are associated with steroidogenesis. 17βHSD14 converts E2 into E1 and T into androstenedione (A4), while 17βHSD12 performs the reverse reactions, E1 into E2, and A4 into T (Tarrant et al., 2003; Luu-The et al., 2006; Shikina et al., 2016; Salah et al., 2019). Three GOIs (CTSL, EI, NPC2) are associated with lysosome/endosomes and/or symbiosome/phagosome (Dani et al., 2014; Mohamed et al., 2016; Dani et al., 2017; Hambleton et al., 2019). One GOI (C3) is associated with innate immunity (Miller et al., 2007; Kvennefors et al., 2010; Palmer and Traylor-Knowles, 2012; Ocampo et al., 2015). Four GOIs (PTCH3, SMO, DHH, GLI2) are representative of Hedgehog (HH) signaling which is highly conserved within vertebrate mammals (Sinkovics, 2015). One GOI (VTG) is the precursor of the egg yolk protein vitellin (Matozzo et al., 2008; Shikina et al., 2013; Lotan et al., 2014).
TABLE 2
| Gene ID | Ligand | Binding energy (kcal/mol) | Binding pocket | ||
|---|---|---|---|---|---|
| Total # of aa/total # contacts | |||||
| Cnid | CS | Cnidarian (accession #) | Crystal structure (PDB file) | ||
| 17βHSD14 | Estradiol | −8.8 | −8.4 | 16aa/238 (KXJ20962) | 17aa/243 (5HS6) |
| Testosterone | −8 | −8.5 | |||
| Cholesterol | −7.5 | −6.4 | |||
| BBP | −7.8 | −7.2 | |||
| BP-3 | −7.4 | −7.4 | |||
| 17βHSD12 | Cholesterol | −6.9 | −7.4 | 17aa/217 (KXJ12187) | 22aa/202 (2ET6) |
| Estradiol | −6.4 | −7.5 | |||
| Testosterone | −6.4 | −7.3 | |||
| BBP | −5.2 | −6.1 | |||
| BP-3 | −5.8 | −5.7 | |||
| PTCH3* | Cholesterol | −9.4 | −7.3 | 25aa/281KXJ20037 | 14aa/2136OEU |
| Estradiol | −7.9 | −7.8 | |||
| Testosterone | −8.5 | −8.1 | |||
| BBP | −6.9 | −6.2 | |||
| BP-3 | −6.6 | −7.0 | |||
| NPC2 | Cholesterol | −9.5 | −12.1 | 18aa/177KXJ29862 | 19aa/2375KWY |
| Estradiol | −8.5 | −10.6 | |||
| Testosterone | −8.7 | −10.3 | |||
| BBP | −6.9 | −8.5 | |||
| BP-3 | −6.6 | −8.0 | |||
| SMO | Cholesterol | −8.4 | −9.9 | 12aa/190KXJ11374 | 20aa/2106XBM |
| Estradiol | −7.6 | −9.0 | |||
| Testosterone | −8.2 | −9.2 | |||
| BBP | −6.6 | −7.9 | |||
| BP-3 | −6.5 | −7.5 | |||
Binding affinities (in kcal/mol) of ligands to proteins, estimated by docking simulations. Values represent the lowest energy binding mode for each docking experiment. Cnid: cnidarian; CS: crystal structure. Total number of amino acids and the total number of physical contacts represent binding pocket interactions. * indicates no direct homolog for PTCH3 so PTCH1 PDB file (6OEU) was used.
Transcriptional Responses of Genes of Interests
For the GOIs associated with steroidogenesis, 17βHSD14 expression was significantly up-regulated (p < 0.01) in each exogenous chemical treatment compared to the control, whereas 17βHSD12 expression was significantly up-regulated (p < 0.01) only in the BP-3 and BBP treatments (Figure 1). VTG exhibited significant down-regulated expression (p < 0.01) in E2, T, and cholesterol treatments (see Figure 1).
FIGURE 1

Expression profile for transcripts associated with sex hormone steroidogenesis or gametogenesis. The ∆∆Cq values represent transformed expression of a GOI relative to RPL11 expression. An * represents a treatment that was significantly different (p < 0.01) relative to the control condition. Error bars represent ± SE. E2: estradiol; T: testosterone; BP-3: oxybenzone; BBP: benzyl butyl phthalate; l7§HSD14: 17a-hydroxysteroid dehydrogenase type 14; 17bHSD12: 17b-hydroxysteroid dehydrogenase type 12: VTG: Vitellogenin 2
GOIs associated with phagocytosis had significant differences in expression in various treatments when compared to the control condition (Figure 2). Cathepsin (CTSL) expression was significantly up-regulated in the BP-3 and BBP treatments (p < 0.01). Equistatin (EI) was significantly up-regulated in E2, BP-3, and BBP treatments (p < 0.01) and significantly down-regulated in the T and cholesterol treatments (p < 0.01) (see Figure 2). NPC2 was significantly down-regulated in cholesterol (p < 0.01) and significantly up-regulated (p < 0.01) in BP-3 and BBP treatments (see Figure 2). Lastly, C3 was significantly up-regulated in E2, T, and BP-3 treatments (see Figure 2).
FIGURE 2

Expression profile of transcripts associated with phagocytosis and cellular structures such as symbiosome, lysosome, or endosome. The ∆∆Cq values represent transformed expression of a GOI relative to RPL11 expression. An * represents a treatment that was significantly different (p < 0.01) relative to the control condition. Error bars represent ± SE. E2: estradiol; T: testosterone; BP-3: oxybenzone; BBP: benzyl butyl phthalate. CSTL: cathepsin L; EI: equistatin; NPC2: Niemann-Pick C type 2; C3: complement component C3.
For the GOIs associated with HH signaling, PTCH3 was significantly down-regulated (p < 0.01) in response to all three sterols, but significantly up-regulated (p < 0.01) in BP-3 and BBP treatments (Figure 3). SMO was significantly down-regulated (p < 0.01) in all treatments except E2 (see Figure 3). DHH was significantly down-regulated (p < 0.01) in T and cholesterol treatments, but significantly up-regulated (p < 0.01) in BP-3 and up-regulated in BBP treatments (p < 0.05) (see Figure 3). GLI2 was significantly up-regulated (p < 0.01) in all treatments except E2 (see Figure 3).
FIGURE 3

Expression profiles of transcripts associated with Hedgehog signaling pathway. The ∆∆Cq values represent transformed expression of GOI relative to RPL11 expression. An * represents a treatment that was significantly different (p < 0.01) relative to the control condition. Error bars represent ± SE. E2: estradiol; T: testosterone; BP-3: oxybenzone; BBP: benzyl butyl phthalate. PTCH3: patched domain-containing protein 3; SMO: smoothened-like; DHH: Desert hedgehog protein A; GLI2: zinc finger protein GLi2.
Univariate ANOVA reveals all sterol treatments were significantly different from each other as well as the xenobiotics (p < 0.01). While the BP-3 and BBP treatments were significantly different from all the sterols (p < 0.01), they were not significantly different from each other (Figure 4).
FIGURE 4

Unique expression profiles for each treatment. (A) control, (B) cholesterol, (C) estradiol, (D) testosterone, (E) oxybenzone, (F) benzyl butyl phthalate. The ∆∆Cq values represent transformed expression of GOI relative to RPL11 expression. Error bars represent ± SE.
Docking Simulations
In order to better understand how sterols and xenobiotics are able to interact with targeted proteins, in silico modelling was used to predict the most energetically favorable ligand/protein binding interactions. The binding affinities of each ligand to the proteins of interest were estimated using docking simulations. Comparative genomic analyses indicate that cnidarian genomes show important similarities to vertebrates in gene content, genomic structure, and organization (Putnam et al., 2007). However, the cnidarian proteome is not as well characterized as the human proteome. Therefore, macromolecular structures needed for the in silico modelling were generated in one of two ways. The first model was the structure in the Protein Databank with the highest degree of homology to the GOI. These PDB files were: 6oeu (PTCHD3), 5kwy (NPC2), 6xbm (SMO), 5hs6 (17βHSD14), and 2et6 (17βHSD12). The best PDB structures corresponded to the human sequences in all cases except for 17βHSD12, where the protein with the closest homology was from Candida tropicalis. The second model was generated using SWISS-MODEL, with the Exaiptasia sequence as input and the closest PDB structure as a template. Each ligand was docked to both protein models. The docking algorithm used flexible fitting with an energy-minimized ligand to determine possible binding conformations and relative binding affinities. All ligands bind to each protein, although the strength of the interactions varied. The binding free energy of the best docking pose for each ligand/receptor combination is shown in Table 3.
TABLE 3
| Gene ID | Primers | Amplicon length |
|---|---|---|
| 17βHSD14 | F: TGCACCCTTTGTTGTGACAT | 209bp |
| R: GATGGCATCCTCCAGAAAGA | ||
| 17βHSD12 | F: AGTCCAGATTTTCTTGCAACCA | 226bp |
| R: TAGACTTCAGTGGTGGGCAG | ||
| VTG | F: GCTGTAGTGGTTTTGGTCGG | 198bp |
| R: TGGTGCTTCTTGGCTTGTTC | ||
| CTSL | F: CATTGCCATTGCATTGATTC | 215bp |
| R: CTGCAAATGCCTACAAGCAA | ||
| EI | F: AGTTGCCCTGGTTTCAAAGA | 200bp |
| R: CCGTCGTCTGTACATTGTGG | ||
| NPC2 | F: TCTTGCAGTTGCCACTTGAC | 204bp |
| R: AATGTTACCGATGCCGAGTC | ||
| C3 | F: TTATCATGGTCCTGGGTGCT | 208bp |
| R: GCGTCAAACTCGAACGTTTT | ||
| PTCH3 | F: TGGATGATTGAGGCTGTGGT | 180bp |
| R: CCTACGCAGCCATTTCCATC | ||
| SMO | F: GAACAGGGTTGGTTGCTCAG | 174bp |
| R: ATTGAAGGCGCTGCTGTTAG | ||
| DHH | F: CGCGTCCTCTCCCTAAACTA | 161bp |
| R: CCCACTCCAACATTCTCCCT | ||
| GLI2 | F: GTGTGTGAAATGCAGCCTCA | 191bp |
| R: GCATCACCTGTCAAGTCCAC |
Genes of Interest and their corresponding primers used in qPCR reaction. Annealing temperature for all primers was 60°C.
Since a crystal structure of 17βHSD14 with a similar ligand (estrone) is available on the PDB (5HS6), it is possible to compare the binding site of an experimentally-determined protein-ligand structure with our docked models (Laskowski and Swindells, 2011). Figure 5 shows a comparison between human 17βHSD14 with bound estrone, human 17βHSD14 with docked estradiol, and Exaiptasia 17βHSD14 with docked estradiol. The search space was confined to the binding pocket for each protein-ligand interaction, and the orientation of the ligand within the pocket was fairly well-conserved (see Figure 5). Additionally, the residues that interact with each ligand were well-conserved between cnidarian and homologous models. For example, estrone and estradiol formed hydrophobic contacts with His93, Gln148, Trp192, and Leu195 in both crystal structure and the docked model of human 17βHSD14. In the cnidarian protein, these residues correspond to His102, Gln157, Trp201, and Leu204, and are predicted to interact the estradiol ligand in the docked model as well (see Figure 5).
FIGURE 5

Comparison of the docking poses of 17§HSD14 with similar ligands generated using LigPlot+ v1.4 software (Laskowski and Swindells 2011). Circled residues indicate side chains that are involved in docking for at least two of the three models. (A) crystal structure of estrone bound to human 17§HSD14 (5HS6). (B) crystal structure of human 17bHSD14 docked to estradiol. (C) crystal structure of Exaiptasia l7§HSD14 docked to estradiol.
Discussion
Attempts at detecting endocrine disruption in cnidarians has been enigmatic. Genomic studies indicate that cnidarians lack CYP19a1a (aka: aromatase) that vertebrates use for sex hormone steroidogenesis (Goldstone, 2008), and yet some cnidarians have demonstrated aromatase activity (Twan et al., 2003). This study sought to identify a suite of GOIs that exhibit differential transcriptional expression in anemones exposed to chemicals known to be endocrine disruptors. The EDCs used in this investigation are lipophilic and detected in the marine environment including sediments and sewage effluent in ng/L to µg/L concentrations (Atkinson et al., 2003; Singh et al., 2010; French et al., 2015; Wear and Thurber, 2015; Al-Jandal et al., 2018; Mitchelmore et al., 2019). Furthermore, all GOIs in this investigation have direct interactions or indirect associations with sterols. Using cholesterol in this investigation was critical as a positive control for evaluating GOIs binding of, and/or interactions with, cholesterol. Transcriptional activities and binding affinities of 17βHSDs in the presence of E2 and T are positive controls since these steroidal hormones are the natural ligands for 17βHSDs.
Links to Steroidogenesis
Endocrine disruption involves chemicals that alter steroidogenesis which interfere with the hormonal signaling leading to down-stream physiological abnormalities in reproduction, development, increased risk of disease, and/or interference with immune and nervous system functions (Heindel et al., 2012). The 17β-hydroxysteroid dehydrogenases (17βHSDs) are involved with steroidogenesis of both estrogens and androgens. 17βHSDs have been detected in cnidarians (Tarrant et al., 2003; Blomquist et al., 2006; Shikina et al., 2016). 17βHSD12 and 17βHSD14 represent the last step in the steroidogenesis pathway prior to the aromatization of T into E2. 17βHSD12 reduces estrone (E1) into E2 as well as androstenedione (A4) into T, while 17βHSD14 performs the opposite reaction by oxidizing E2 into E1, and T into A4 (Luu-The et al., 2006; Caldwell et al., 2012; Sivik et al., 2012; Thomas and Potter, 2013; Shikina et al., 2016; Salah et al., 2019). Concentrations of competing sterols may drive which sterol gets converted by a specific 17βHSD (Lathe and Kotelevtsev, 2014). E2 concentration is known to regulate transcription of 17βHSDs (Rotinen et al., 2009). Transcriptional profiles herein demonstrate the organism (regardless of sex) responds to the environmental “overload” of E2 or T by modulating the transcription of the enzyme in order to transform both sex hormones into less potent sterols. These profiles also reveal that cholesterol and xenobiotics (BBP and BP-3) interfere with steroidogenesis by altering normal transcriptional demands for 17βHSDs (see Figure 1).
Links to Sterol Binding and Transport
Cnidarians get cholesterol and other sterols by either dietary intake and/or from their endosymbiont (Revel et al., 2016; Hambleton et al., 2019). Cholesterol transport and metabolism is tightly regulated (Maxfield et al., 2016). Due to lipid-rich tissues, it has been proposed that lipophilic compounds will easily diffuse into cnidarian tissues (Peters et al., 1997). Sterols are most concentrated in cnidarian epidermis and some sterols are known to increase in concentration after exposure to environmental pollutants (Revel et al., 2016; Stien et al., 2020). The lack of significant expression of C3 in the cholesterol treatment is consistent with the fact that it is a natural constituent of animal membranes.
Five GOIs (C3, PTCH3, CTSL, EI, NPC2) have functional associations with the endocytic processes whereas a sixth GOI (DHH) is a secreted morphogen associated with exocytosis/secretion. C3 is associated with phagocytosis of foreign particles into cnidarians (Palmer and Traylor-Knowles, 2012), while PTCH regulates functions of: receptor-mediated endocytosis of sterol-protein ligands; down-stream morphogens; and the removal of oxysterols (Incardona et al., 2000; Strutt et al., 2001; Corcoran and Scott, 2006; Zhong et al., 2014). Three GOIs (CTSL, EI, NPC2) are associated with lysosomes and phagocytosis. Cathepsins are lysosomal proteases and EI is the cnidarian homolog of an inhibitor of cathepsins (Lenarčič et al., 1997; Štrukelj et al., 2000). Lysosomes contribute to the termination of hormone effects by degrading proteins interacting with the hormone (Qualmann et al., 2000; Totta et al., 2014). NPC2 functions as a key regulator of sterol transport between Symbiodinium and cnidarian hosts which is representative of stable symbioses in a structure known as the symbiosome (Infante et al., 2008; Vance, 2010; Ganot et al., 2011; Dani et al., 2014; Bucher et al., 2016; Revel et al., 2016; Hambleton et al., 2019; Lu et al., 2020). Seven GOIs (NPC2, 17βHSD14, C3, PTCH3, SMO, DHH, VTG) have known linkages to endoderm. (Kvennefors et al., 2010; Shikina et al., 2013; Sanders et al., 2014; Shikina et al., 2016).
Five GOIs (17βHSD14, 17βHSD12, NPC2, PTCH3, SMO) that have sterol binding pockets were used in silico modelling. In silico modelling has previously demonstrated phthalates and steroids without an aromatic A-ring are capable of binding to sterol-binding proteins (Sheikh and Beg, 2017; Khalturin et al., 2018). Patched genes have multiple functionally redundant homologs that cannot discriminate between ligands (Carpenter et al., 1998; Zhulyn et al., 2015; Hu and Song, 2019). Justification for model docking is based on the fact that cnidarian genomes have representative vertebrate homologs of relevant proteins (Miller et al., 2007; Putnam et al., 2007; Kvennefors et al., 2010). In silico modelling demonstrates the lowest energy docking poses for the crystal structures of the homologous proteins and the models of the Exaiptasia proteins correlate well with each other (see Table 3). The binding scores for the Exaiptasia protein models generally follow the same trends as the homologs from the PDB. In both cases, the steroid ligands (cholesterol, E2, and T) bound more strongly than BP-3 or BBP, although the differences in binding energy seem to be smaller for the 17βHSD12 and 17βHSD14 proteins. This result was expected, since BBP and BP-3 are considerably smaller and as such, would make less contact with binding pocket residues (see Table 3). Although there are a few minor differences between homologs, such as hydrogen bonding patterns, the orientation of the ligand in the 17βHSD14 binding pocket and the side chains involved in binding are strikingly similar (see Figure 5). All sterols have similarities in orientations, binding affinities, and expression profiles (see Table 3 and Figure 5). 17βHSDs are known to exhibit variability in substrate binding (Blomquist et al., 2006; Lathe and Kotelevtsev, 2014).
Data herein provides new information about the diversity of ligands, including xenobiotics that bind to NPC2, PTCH3, SMO, and the 17βHSDs. The differential expression of these sterol-binding proteins suggests that these proteins handle natural sterols and xenobiotics differently. Laboratory-induced exposure to cholesterol, E2, and T requires modulating transcript copy numbers to handle the environmental flux of natural ligands. However, laboratory-induced exposures of BBP and BP-3 appear to require transcriptional replacement for proteins after binding xenobiotics. All GOIs except SMO exhibit up-regulated expression in the BP-3 and BBP treatments (see Figures 4E,F).
Links to Development
Hedgehog (HH) signaling is a highly conserved developmental pathway that relies on cholesterol interacting with specific sterol-binding proteins. Endocrine disruption is known to interfere with HH signaling (Johansson and Svingen, 2020). Developmentally-regulated signaling pathways involving apoptosis, embryogenesis, cell proliferation, and development of diseases present in higher animals are also found in cnidarians including Exaiptasia (Technau et al., 2005; Matus et al., 2008; Sinkovics, 2015; Bucher et al., 2016). Four GOIs (PTCH3, SMO, DHH, GLI2) can be linked to HH signaling. E2 acts as a direct HH pathway antagonist (Lipinski and Bushman, 2010) and the E2 expression profile in Figure 3 reaffirms this inhibitory effect. Recognized inhibitors of HH signaling also include the well-characterized steroidal alkaloid cyclopamine (Chen 2002). The down-regulated expression profiles of PTCH3 in the E2, T, and cholesterol treatments and the down-regulated expression of SMO in the T and cholesterol treatments suggests all three sterols exhibit similar inhibitory effects (see Figure 3). The HH pathway is also responsive to xenobiotics as well as sterols (Johansson and Svingen, 2020). Phthalates are known to impact DHH and PTCH expression, and BP-3 is considered estrogenic (Thibaut and Porte, 2004; Borch et al., 2006; Downs et al., 2016; Hong et al., 2016; LaPlante et al., 2018; Matouskova et al., 2020). Although PCTH3 and DHH were down-regulated in response to sterols, PTCH3, DHH, and GLI2 were up-regulated in response to BP-3 and BBP (see Figure 3), suggesting opposing effects of sterols and xenobiotics on HH signaling, presumably with sterols exerting antagonists effects and xenobiotics exerting agonist effects. Regardless, the expression profiles in Figure 3 highlight the inherent vulnerability of the HH pathway to diverse chemical signals.
Links to Gametogenesis
Gametogenesis is associated with HH signaling. HH signaling molecules are considered stem-cell factors (Zhang and Kalderon, 2001). Evidence exists in other metazoans that secreted molecules in the HH signaling pathway are one of the primary regulators of stem cells which are necessary precursors of oogenesis (Zhang and Kalderon, 2001). Shikina and Chang (2016) identifies cnidarian germline stem cells as originating in endodermal mesentery. HH signaling GOIs (PTCH3, SMO, DHH) are known to be restricted to the endoderm of male and female gonophores (Sanders et al., 2014). Germline stem cells are the essential cells responsible for this process in various metazoans. In Drosophila, somatic cells in the ovary of cannot proliferate normally in the absence of Smoothened activity (Zhang and Kalderon, 2000), and excessive HH signaling can result from aberrant expansion of stem cell pools (Zhang and Kalderon, 2001). Results show all exogenous treatments (except E2) caused down-regulation of SMO (see Figure 4). With the exception of PTCH3, results also show that HH signaling GOIs in the E2 treatment are not significantly expressed (see Figure 3).
Exaiptasia reproductive strategies include both sexual and asexual reproduction (Clayton, 1985; Grawunder et al., 2015). Subsets of anemones used for each treatment were taken from the same population as the control anemones. Although not confirmed, it is highly likely that all anemones used in this investigation were clonal. Detecting the presence of sex-specific tissues enhances the ability to interpret anemone responses to the EDCs. Although the sex of the colony was not known, expression profiles of the GOIs could provide clues about the sex of these animals. Vitellogenin is a female-specific protein commonly found across the animal kingdom and is the most abundant protein associated with oocyte development in cnidarians (Shikina et al., 2013; Levitan et al., 2015). Hormonal and/or environmental cues are known to stimulate release of oocytes (Levitan et al., 2015). Exogenous E2 is known to influence pedal laceration and egg bundle formation in cnidarians (Tarrant et al., 2004; Thorn et al., 2015). Blue-light cues have also been used to induce Aiptasia sp gametogenesis in a laboratory setting (Grawunder et al., 2015).
This investigation also did not attempt to identify the sex of the anemones nor their reproductive status, per se. Hermaphroditism and even Trioecy has been observed in closely related species (Schlesinger et al., 2010; Armoza-Zvuloni et al., 2014). A small amount of data from a previous study suggests male tissues from hermaphroditic corals also express VTG (Hayakawa et al., 2005). Results reveal VTG expression in the control treatment suggesting anemones are female and already transcribing VTG (see Figure 1). Decreases in vitellogenin in females can be useful for identifying EDCs (Ankley et al., 2002). If these anemones were male, then up-regulation of vitellogenin would have been expected. Up-regulation of vitellogenin in males is widely recognized as the prominent estrogenic effect in many aquatic animals exposed to exogenous E2 (Denslow et al., 1999; Matozzo et al., 2008). These results cannot confirm nor reject the possibility that male reproductive tissues were also present.
A second piece of evidence suggest that these anemones were female was the expression of 17βHSD14 in E2 which was significantly up-regulated (p < 0.0001) compared to the T, cholesterol, BBP, and BP-3 treatments. Female anemones perform oogenesis which is stimulated by the presence of E2 so when exogenous E2 is introduced, additional 17βHSD14 is required to reduce E2 potency by transformation into E1. In this investigation, VTG expression was crucial for differentiating E2 exposure from other EDC exposures. The expression profiles of 17βHSD14 and 17βHSD12 can be indirect measures of cellular demands for sex hormones. When exposed to exogenous sex hormones, transcription of 17βHSD14 is up-regulated while 17βHSD12 is not significantly different from the control. These results suggest these anemones were females responding to the EDCs in this investigation. Future investigations using these GOIs, coupled with lab-controlled methodologies to induce gametogenesis in recognized clonal lines of anemones should provide greater resolution in how these organisms respond to EDCs.
Conclusion
This data demonstrates that endocrine disrupting exogenous sterols (E2, T, cholesterol) and xenobiotics (BP-3, BBP) alter transcription of genes associated with steroidogenesis, sterol transport, oogenesis, and the Hedgehog signaling pathway (Figures 1–4). In silico modelling demonstrates that EDC’s bind favorably into pockets of proteins involved with steroidogenesis, cholesterol transport, and HH signaling. Because these GOIs were expressed at basal levels in the control treatment (i.e. unstressed conditions) and significantly differentially expressed when exposed to EDCs (Figure 4), they may serve as transcriptional biomarkers. Biomarkers have the potential to link field and laboratory studies by acting as ‘‘mechanistic signposts’’ of exposure to anthropogenic chemicals (Hutchinson, 2007; Traylor-Knowles and Palumbi, 2014). The most validated biomarker of estrogenic exposure is VTG, an estrogen-dependent glycolipophosphoprotein naturally expressed only in oviparous females during the reproductive season (Verderame and Scudiero, 2017). The biomarkers used in this investigation offer new insight into BBP and BP-3 exposures. Understanding mechanisms of action for xenobiotics can be advanced by utilizing known biochemical pathways involving cholesterol and sex hormones.
Collectively, the expression profiles of this suite of GOIs are capable of discriminating seven distinct outcomes. 1) It is possible to discern differences between control vs exposure to exogenous chemicals (i.e. unstressed vs stressed) (see Figures 1–4). 2) Expression profiles of GOIs involved with HH signaling can differentiate E2 exposure from other sterols (see Figure 3). 3) Cholesterol exposure can be discriminated from sex hormones by expression of C3 (see Figures 4B–D). 4) C3 expression can also differentiate BP-3 versus BBP exposure (see Figure 2). 5) It is possible to differentiate natural sterols from xenobiotic exposure (see Figure 4). 6) Expression profiles of VTG and 17βHSD 12/14 provide a mechanism for discerning the responses of female reproductive tissues in E2 exposure compared to T or cholesterol exposures (see Figure 1). 7) Expression of 17βHSD14 and 17βHSD12 were critical for discerning BBP and BP-3 treatments compared to exogenous sterols (see Figure 1). Coupled with expression profiling, in silico results reveal proteins that bind cholesterol or modify sex hormones are vulnerable to endocrine disruption in the presence of either xenobiotics or exogenous sterols.
Future ecotoxicological investigations using more sensitive analytical measurements of chemical concentrations should enable identifying minimal concentrations necessary to alter expression of these GOIs. A primary goal of this investigation was to identify GOIs that are responsive to these EDCs. Moving forward, future work can focus on measuring GOI responses over recognized acute toxicity timeframes such as 24–96 hrs. Lastly, it will also be important to establish the baseline temporal variance in GOI expression before, during, and after the peak reproductive season.
In summary, this investigation provides new insight into anemone responsiveness to a small set of chemicals recognized as endocrine disruptors. Individually, each GOI provides a clue about a specific biochemical pathway, cellular process, and/or developmental pathway. Collectively, these GOIs offer greater resolution into understanding how anemones respond to different EDCs. More than a decade ago, Tarrant (2007) questioned if sex hormones are disrupting or even impacting cnidarians and the enigma of aromatase remains. However, the data herein offers new information towards deciphering endocrine disruption in cnidarians.
Statements
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.
Author contributions
MM, JR, and JW conceived the investigation. JR, HY, JW, JE, RP, MM, DQ provided resources and financial support. JR, HY, JW, JE, MM, RP, DQ performed data collection and computational analysis. JW, HY, MM, DQ performed bioinformatics. MM wrote the manuscript with contributions from DQ, JR, JW, and HY. All authors approved the final manuscript.
Funding
Financial support was provided by Berry College for three Richards Scholarships (JR; JW; HY), and two Richards Undergraduate Support Grants (RP; JE).
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
Al-JandalN.SaeedT.AzadI.Al-SubiaiS.Al-ZekriW.HussainS.et al (2018). Impact of Endocrine Disrupting Compounds in Sewage Impacted Coastal Area on Seabream. Ecotoxicology Environ. Saf.150, 280–288. 10.1016/j.ecoenv.2017.12.035
2
AnkleyG. T.KahlM. D.JensenK. M.HornungM. W.KorteJ. J.MakynenE. A.et al (2002). Evaluation of the Aromatase Inhibitor Fadrozole in a Short-Term Reproduction Assay with the Fathead Minnow (Pimephales promelas). Toxicol. Sci.67 (1), 121–130. 10.1093/toxsci/67.1.121
3
ArcherE.PetrieB.Kasprzyk-HordernB.WolfaardtG. M. (2017). The Fate of Pharmaceuticals and Personal Care Products (PPCPs), Endocrine Disrupting Contaminants (EDCs), Metabolites and Illicit Drugs in a WWTW and Environmental Waters. Chemosphere174, 437–446. 10.1016/j.chemosphere.2017.01.101
4
Armoza-ZvuloniR.Kramarsky-WinterE.LoyaY.SchlesingerA.RosenfeldH. (2014). Trioecy, a Unique Breeding Strategy in the Sea Anemone Aiptasia Diaphana and its Association with Sex Steroids1. Biol. Repro90 (6), 1–8. 10.1095/biolreprod.113.114116
5
AtkinsonS.AtkinsonM. J. (1992). Detection of Estradiol-17? during a Mass Coral Spawn. Coral Reefs11 (1), 33–35. 10.1007/bf00291932
6
AtkinsonS.AtkinsonM. J.TarrantA. M. (2003). Estrogens from Sewage in Coastal Marine Environments. Environ. Health Perspect.111 (4), 531–535. 10.1289/ehp.5233
7
AxworthyJ. B.Padilla-GamiñoJ. L. (2019). Microplastics Ingestion and Heterotrophy in Thermally Stressed Corals. Sci. Rep.9 (1), 18193. 10.1038/s41598-019-54698-7
8
BahrK. D.JokielP. L.ToonenR. J. (2015). The Unnatural History of Kāne'ohe Bay: Coral Reef Resilience in the Face of Centuries of Anthropogenic Impacts. PeerJ3, e950. 10.7717/peerj.950
9
BargarT. A.GarrisonV. H.AlvarezD. A.EcholsK. R. (2013). Contaminants Assessment in the Coral Reefs of Virgin Islands National Park and Virgin Islands Coral Reef National Monument. Mar. Pollut. Bull.70 (1-2), 281–288. 10.1016/j.marpolbul.2013.03.001
10
BlomquistC. H.LimaP. H.TarrantA. M.AtkinsonM. J.AtkinsonS. (2006). 17β-Hydroxysteroid Dehydrogenase (17β-HSD) in Scleractinian Corals and Zooxanthellae. Comp. Biochem. Physiol. B: Biochem. Mol. Biol.143 (4), 397–403. 10.1016/j.cbpb.2005.12.017
11
BorchJ.MetzdorffS. B.VinggaardA. M.BrokkenL.DalgaardM. (2006). Mechanisms Underlying the Anti-androgenic Effects of Diethylhexyl Phthalate in Fetal Rat Testis. Toxicology223 (1-2), 144–155. 10.1016/j.tox.2006.03.015
12
BrennanJ. J.MesserschmidtJ. L.WilliamsL. M.MatthewsB. J.ReynosoM.GilmoreT. D. (2017). Sea Anemone Model Has a Single Toll-like Receptor that Can Function in Pathogen Detection, NF-Κb Signal Transduction, and Development. Proc. Natl. Acad. Sci. USA114 (47), E10122–E10131. 10.1073/pnas.1711530114
13
BrunoJ. F.SeligE. R. (2007). Regional Decline of Coral Cover in the Indo-Pacific: Timing, Extent, and Subregional Comparisons. PLoS One2 (8), e711. 10.1371/journal.pone.0000711
14
BucherM.WolfowiczI.VossP. A.HambletonE. A.GuseA. (2016). Development and Symbiosis Establishment in the Cnidarian Endosymbiosis Model Aiptasia Sp. Sci. Rep.6 (1), 19867. 10.1038/srep19867
15
BurnsK.BrinkmanD. (2011). Organic Biomarkers to Describe the Major Carbon Inputs and Cycling of Organic Matter in the central Great Barrier Reef Region. Estuarine, Coastal Shelf Sci.93 (2), 132–141. 10.1016/j.ecss.2011.04.001
16
BustinS. A.BenesV.GarsonJ. A.HellemansJ.HuggettJ.KubistaM.et al (2009). The MIQE Guidelines: Minimum Information for Publication of Quantitative Real-Time PCR Experiments. Clin. Chem.55 (4), 611–622. 10.1373/clinchem.2008.112797
17
CaldwellD. J.MastroccoF.AndersonP. D.LängeR.SumpterJ. P. (2012). Predicted-no-effect Concentrations for the Steroid Estrogens Estrone, 17β-Estradiol, Estriol, and 17α-Ethinylestradiol. Environ. Toxicol. Chem.31 (6), 1396–1406. 10.1002/etc.1825
18
CarpenterD.StoneD. M.BrushJ.RyanA.ArmaniniM.FrantzG.et al (1998). Characterization of Two Patched Receptors for the Vertebrate Hedgehog Protein Family. Proc. Natl. Acad. Sci.95 (23), 13630–13634. 10.1073/pnas.95.23.13630
19
ChenJ. K.TaipaleJ.CooperM. K.BeachyP. A. (2002). Inhibition of Hedgehog Signaling by Direct Binding of Cyclopamine to Smoothened. Genes Dev.16 (21), 2743–2748. 10.1101/gad.1025302
20
ClaytonW. (1985). Pedal Laceration by the Anemone Aiptasia Pallida. Mar. Ecol. Prog. Ser.21 (1-2), 75–80. 10.3354/meps021075
21
CorcoranR. B.ScottM. P. (2006). Oxysterols Stimulate Sonic Hedgehog Signal Transduction and Proliferation of Medulloblastoma Cells. Proc. Natl. Acad. Sci.103 (22), 8408–8413. 10.1073/pnas.0602852103
22
DachsJ.BayonaJ. M.IttekkotV.AlbaigésJ. (1999). Monsoon-Driven Vertical Fluxes of Organic Pollutants in the Western Arabian Sea. Environ. Sci. Technol.33 (22), 3949–3956. 10.1021/es990200e
23
DaniV.GanotP.PriouzeauF.FurlaP.SabouraultC. (2014). Are Niemann-Pick Type C Proteins Key Players in Cnidarian-Dinoflagellate Endosymbioses?Mol. Ecol.23 (18), 4527–4540. 10.1111/mec.12876
24
DaniV.PriouzeauF.MertzM.MondinM.PagnottaS.Lacas-GervaisS.et al (2017). Expression Patterns of Sterol Transporters NPC1 and NPC2 in the Cnidarian-Dinoflagellate Symbiosis. Cell Microbiol.19 (10), e12753. 10.1111/cmi.12753
25
DengY.YanZ.ShenR.WangM.HuangY.RenH.et al (2020). Microplastics Release Phthalate Esters and Cause Aggravated Adverse Effects in the Mouse Gut. Environ. Int.143, 105916. 10.1016/j.envint.2020.105916
26
DenslowN. D.BowmanC. J.RobinsonG.LeeH. S.FergusonR. J.HemmerM. J.et al (1999). “Biomarkers of Endocrine Disruption at the mRNA Level,” in Environmental Toxicology and Risk Assessment: Standardization of Biomarkers of Endocrine Disruption and Environmental Assessment. Editors HenshelD. S.BlackM. C.HarrassM. C. (West Conshohocken: American Society for Testing Materials), 8, 24–35.
27
DiNardoJ. C.DownsC. A. (2018). Dermatological and Environmental Toxicological Impact of the Sunscreen Ingredient Oxybenzone/benzophenone-3. J. Cosmet. Dermatol.17 (1), 15–19. 10.1111/jocd.12449
28
DownsC. A.Kramarsky-WinterE.SegalR.FauthJ.KnutsonS.BronsteinO.et al (2016). Toxicopathological Effects of the Sunscreen UV Filter, Oxybenzone (Benzophenone-3), on Coral Planulae and Cultured Primary Cells and its Environmental Contamination in Hawaii and the U.S. Virgin Islands. Arch. Environ. Contam. Toxicol.70 (2), 265–288. 10.1007/s00244-015-0227-7
29
DuffyD. J.FrankU. (2011). Modulation of COUP-TF Expression in a Cnidarian by Ectopic Wnt Signalling and Allorecognition. PLoS One6 (4), e19443. 10.1371/journal.pone.0019443
30
DunnS. R.PhillipsW. S.SpataforaJ. W.GreenD. R.WeisV. M. (2006). Highly Conserved Caspase and Bcl-2 Homologues from the Sea Anemone Aiptasia Pallida: Lower Metazoans as Models for the Study of Apoptosis Evolution. J. Mol. Evol.63 (1), 95–107. 10.1007/s00239-005-0236-7
31
EdgeS. E.ShearerT. L.MorganM. B.SnellT. W. (2013). Sub-lethal Coral Stress: Detecting Molecular Responses of Coral Populations to Environmental Conditions over Space and Time. Aquat. Toxicol.128-129, 135–146. 10.1016/j.aquatox.2012.11.014
32
EdmanC. F.PrigentS. A.SchipperA.FeramiscoJ. R. (1997). Identification of ErbB3-Stimulated Genes Using Modified Representational Difference Analysis. Biochem. J323, 113–118. 10.1042/bj3230113
33
FabriciusK. E. (2005). Effects of Terrestrial Runoff on the Ecology of Corals and Coral Reefs: Review and Synthesis. Mar. Pollut. Bull.50 (2), 125–146. 10.1016/j.marpolbul.2004.11.028
34
FloodV. S.PittJ. M.SmithS. R. (2005). Coral Communities on Patch Reefs in Castle Harbour Bermuda after More Than a century of Perturbations. Mar. Poll Bull.51 (5-7), 545–557. 10.1016/j.marpolbul.2005.07.010
35
FrenchV. A.Codi KingS.KumarA.NorthcottG.McGuinnessK.ParryD. (2015). Characterisation of Microcontaminants in Darwin Harbour, a Tropical Estuary of Northern Australia Undergoing Rapid Development. Sci. Total Environ.536, 639–647. 10.1016/j.scitotenv.2015.07.114
36
FuentesN.SilveyraP. (2019). Estrogen Receptor Signaling Mechanisms. Adv. Protein Chem. Str116, 135–170. 10.1016/bs.apcsb.2019.01.001
37
GanotP.MoyaA.MagnoneV.AllemandD.FurlaP.SabouraultC. (2011). Adaptations to Endosymbiosis in a Cnidarian-Dinoflagellate Association: Differential Gene Expression and Specific Gene Duplications. Plos Genet.7 (7), e1002187. 10.1371/journal.pgen.1002187
38
GaoD.-W.WenZ.-D. (2016). Phthalate Esters in the Environment: A Critical Review of Their Occurrence, Biodegradation, and Removal during Wastewater Treatment Processes. Sci. Total Environ.541, 986–1001. 10.1016/j.scitotenv.2015.09.148
39
GardnerT. A.CôtéI. M.GillJ. A.GrantA.WatkinsonA. R. (2003). Long-Term Region-wide Declines in Caribbean Corals. Science301 (5635), 958–960. 10.1126/science.1086050
40
GassmanN.KennedyC. (1992). Cytochrome P-450 Content and Xenobiotic Metabolizing Enzyme Activities in the Scleractinian Coral, Favia fragum (Esper). Bull. Mar. Sci.50 (2), 320–330.
41
GavioB.Palmer-CantilloS.ManceraJ. E. (2010). Historical Analysis (2000-2005) of the Coastal Water Quality in San Andrés Island, SeaFlower Biosphere Reserve, Caribbean Colombia. Mar. Pollut. Bull.60 (7), 1018–1030. 10.1016/j.marpolbul.2010.01.025
42
GoldstoneJ. V. (2008). Environmental Sensing and Response Genes in Cnidaria: the Chemical Defensome in the Sea Anemone Nematostella vectensis. Cell Biol Toxicol24 (6), 483–502. 10.1007/s10565-008-9107-5
43
GrawunderD.HambletonE. A.BucherM.WolfowiczI.BechtoldtN.GuseA. (2015). Induction of Gametogenesis in the Cnidarian Endosymbiosis Model Aiptasia Sp. Sci. Rep.5 (1), 15677. 10.1038/srep15677
44
HallN. M.BerryK. L. E.RintoulL.HoogenboomM. O. (2015). Microplastic Ingestion by Scleractinian Corals. Mar. Biol.162 (3), 725–732. 10.1007/s00227-015-2619-7
45
HambletonE. A.JonesV. A. S.MaegeleI.KvaskoffD.SachsenheimerT.GuseA. (2019). Sterol Transfer by Atypical Cholesterol-Binding NPC2 Proteins in Coral-Algal Symbiosis. eLife8, e43923. 10.7554/eLife.43923
46
HanwellM. D.CurtisD. E.LonieD. C.VandermeerschT.ZurekE.HutchisonG. R. (2012). Avogadro: an Advanced Semantic Chemical Editor, Visualization, and Analysis Platform. J. Cheminform4 (1), 17. 10.1186/1758-2946-4-17
47
HayakawaH.NakanoY.AndohT.WatanabeT. (2005). Sex-dependent Expression of mRNA Encoding a Major Egg Protein in the Gonochoric Coral Galaxea Fascicularis. Coral Reefs24 (3), 488–494. 10.1007/s00338-005-0485-3
48
HeindelJ. J.ZoellerR. T.JoblingS.IguchiT.VandenbergL.WoodruffT. J. (2012). “What Is Endocrine Disruption All about,” in State of the Science of Endocrine Disrupting Chemicals - 2012. Editors BergmanÅ.HeindelJ. J.JoblingS.KiddK. A.ZoellerR. T. (Geneva, Switzerland: World Health Organization), 289.
49
HongH.RuaD.SakkiahS.SelvarajC.GeW.TongW. (2016). Consensus Modeling for Prediction of Estrogenic Activity of Ingredients Commonly Used in Sunscreen Products. Ijerph13 (10), 958. 10.3390/ijerph13100958
50
HuA.SongB.-L. (2019). The Interplay of Patched, Smoothened and Cholesterol in Hedgehog Signaling. Curr. Opin. Cel Biol.61, 31–38. 10.1016/j.ceb.2019.06.008
51
HubankM.SchatzD. G. (1999). “cDNA Representational Difference Analysis: A Sensitive and Flexible Method for Identification of Differentially Expressed Genes,” in Methods in Enzymology; cDNA Preparation and Characterization. Editor WeissmanS. M. (San Diego, California, USA: Academic Press), 303, 325–349. 10.1016/s0076-6879(99)03021-9
52
HughesT. P.BairdA. H.BellwoodD. R.CardM.ConnollyS. R.FolkeC.et al (2003). Climate Change, Human Impacts, and the Resilience of Coral Reefs. Science301 (5635), 929–933. 10.1126/science.1085046
53
HutchinsonT. H. (2007). Small Is Useful in Endocrine Disrupter Assessment-Four Key Recommendations for Aquatic Invertebrate Research. Ecotoxicology16 (1), 231–238. 10.1007/s10646-006-0107-z
54
ImbsA. B. (2013). Fatty Acids and Other Lipids of Corals: Composition, Distribution, and Biosynthesis. Russ. J. Mar. Biol.39 (3), 153–168. 10.1134/s1063074013030061
55
IncardonaJ. P.LeeJ. H.RobertsonC. P.EngaK.KapurR. P.RoelinkH. (2000). Receptor-mediated Endocytosis of Soluble and Membrane-Tethered Sonic Hedgehog by Patched-1. Proc. Natl. Acad. Sci. U S A.97 (22), 12044–12049. 10.1073/pnas.220251997
56
InfanteR. E.WangM. L.RadhakrishnanA.KwonH. J.BrownM. S.GoldsteinJ. L. (2008). NPC2 Facilitates Bidirectional Transfer of Cholesterol between NPC1 and Lipid Bilayers, a Step in Cholesterol Egress from Lysosomes. Proc. Natl. Acad. Sci. U S A.105 (40), 15287–15292. 10.1073/pnas.0807328105
57
JohanssonH. K. L.SvingenT. (2020). Hedgehog Signal Disruption, Gonadal Dysgenesis and Reproductive Disorders: Is There a Link to Endocrine Disrupting Chemicals?Curr. Res. Toxicol.1, 116–123. 10.1016/j.crtox.2020.10.001
58
JonesG. P.McCormickM. I.SrinivasanM.EagleJ. V. (2004). Coral Decline Threatens Fish Biodiversity in marine Reserves. Proc. Natl. Acad. Sci.101 (21), 8251–8253. 10.1073/pnas.0401277101
59
JonesR. (2010). Environmental Contamination Associated with a marine Landfill ('seafill') beside a Coral Reef. Mar. Pollut. Bull.60 (11), 1993–2006. 10.1016/j.marpolbul.2010.07.028
60
JonesR.MullerJ.HaynesD.SchreiberU. (2003). Effects of Herbicides Diuron and Atrazine on Corals of the Great Barrier Reef, Australia. Mar. Ecol. Prog. Ser.251, 153–167. 10.3354/meps251153
61
KenkelC. D.AglyamovaG.AlamaruA.BhagooliR.CapperR.CunningR.et al (2011). Development of Gene Expression Markers of Acute Heat-Light Stress in Reef-Building Corals of the Genus Porites. PLoS One6 (10), e26914. 10.1371/journal.pone.0026914
62
KenkelC. D.SheridanC.LealM. C.BhagooliR.CastilloK. D.KurataN.et al (2014). Diagnostic Gene Expression Biomarkers of Coral thermal Stress. Mol. Ecol. Resour.14 (4), 667–678. 10.1111/1755-0998.12218
63
KerswellA.JonesR. (2003). Effects of Hypo-Osmosis on the Coral Stylophora Pistillata: Nature and Cause of 'low-Salinity Bleaching'. Mar. Ecol. Prog. Ser.253, 145–154. 10.3354/meps253145
64
KhalturinK.BillasI. M. L.ChebaroY.ReitzelA. M.TarrantA. M.LaudetV.et al (2018). NR3E Receptors in Cnidarians: A New Family of Steroid Receptor Relatives Extends the Possible Mechanisms for Ligand Binding. J. Steroid Biochem. Mol. Biol.184, 11–19. 10.1016/j.jsbmb.2018.06.014
65
KoF.-C.ChangC.-W.ChengJ.-O. (2014). Comparative Study of Polycyclic Aromatic Hydrocarbons in Coral Tissues and the Ambient Sediments from Kenting National Park, Taiwan. Environ. Pollut.185, 35–43. 10.1016/j.envpol.2013.10.025
66
KvenneforsE. C. E.LeggatW.KerrC. C.AinsworthT. D.Hoegh-GuldbergO.BarnesA. C. (2010). Analysis of Evolutionarily Conserved Innate Immune Components in Coral Links Immunity and Symbiosis. Dev. Comp. Immunol.34 (11), 1219–1229. 10.1016/j.dci.2010.06.016
67
LafontR.MathieuM. (2007). Steroids in Aquatic Invertebrates. Ecotoxicology16 (1), 109–130. 10.1007/s10646-006-0113-1
68
LaPlanteC. D.BansalR.DunphyK. A.JerryD. J.VandenbergL. N. (2018). Oxybenzone Alters Mammary Gland Morphology in Mice Exposed during Pregnancy and Lactation. J. Endocr. Soc.2 (8), 903–921. 10.1210/js.2018-00024
69
LaskowskiR. A.SwindellsM. B. (2011). LigPlot+: Multiple Ligand-Protein Interaction Diagrams for Drug Discovery. J. Chem. Inf. Model.51 (10), 2778–2786. 10.1021/ci200227u
70
LatheR.KotelevtsevY. (2014). Steroid Signaling: Ligand-Binding Promiscuity, Molecular Symmetry, and the Need for Gating. Steroids82, 14–22. 10.1016/j.steroids.2014.01.002
71
LawK. L. (2017). Plastics in the Marine Environment. Annu. Rev. Mar. Sci.9 (1), 205–229. 10.1146/annurev-marine-010816-060409
72
LebretonL.SlatB.FerrariF.Sainte-RoseB.AitkenJ.MarthouseR.et al (2018). Evidence that the Great Pacific Garbage Patch Is Rapidly Accumulating Plastic. Sci. Rep.8 (1), 4666. 10.1038/s41598-018-22939-w
73
LenarčičB.RitonjaA.ŠtrukeljB.TurkV. (1997). Equistatin, a New Inhibitor of Cysteine Proteinases fromActinia Equina, Is Structurally Related to Thyroglobulin Type-1 Domain. J. Biol. Chem.272 (21), 13899–13903. 10.1074/jbc.272.21.13899
74
LevitanS.SherN.BrekhmanV.ZivT.LubzensE.LotanT. (2015). The Making of an Embryo in a Basal Metazoan: Proteomic Analysis in the Sea anemoneNematostella Vectensis. Proteomics15 (23-24), 4096–4104. 10.1002/pmic.201500255
75
LipinskiR. J.BushmanW. (2010). Identification of Hedgehog Signaling Inhibitors with Relevant Human Exposure by Small Molecule Screening. Toxicol. Vitro24 (5), 1404–1409. 10.1016/j.tiv.2010.04.011
76
LotanT.Chalifa-CaspiV.ZivT.BrekhmanV.GordonM. M.AdmonA.et al (2014). Evolutionary Conservation of the Mature Oocyte Proteome. EuPA Open Proteomics3, 27–36. 10.1016/j.euprot.2014.01.003
77
LuY.JiangJ.ZhaoH.HanX.XiangY.ZhouW. (2020). Clade-Specific Sterol Metabolites in Dinoflagellate Endosymbionts Are Associated with Coral Bleaching in Response to Environmental Cues. mSystems5 (5), e00765. 10.1128/msystems.00765-20
78
Luu-TheV.TremblayP.LabrieF. (2006). Characterization of Type 12 17β-Hydroxysteroid Dehydrogenase, an Isoform of Type 3 17β-Hydroxysteroid Dehydrogenase Responsible for Estradiol Formation in Women. Mol. Endocrinol.20 (2), 437–443. 10.1210/me.2005-0058
79
MahonA. M.O’ConnellB.HealyM. G.O’ConnorI.OfficerR.NashR.et al (2017). Microplastics in Sewage Sludge: Effects of Treatment. Environ. Sci. Technol.51 (2), 810–818. 10.1021/acs.est.6b04048
80
MansfieldK. M.CarterN. M.NguyenL.ClevesP. A.AlshanbayevaA.WilliamsL. M.et al (2017). Transcription Factor NF-Κb Is Modulated by Symbiotic Status in a Sea Anemone Model of Cnidarian Bleaching. Sci. Rep.7 (1), 16025. 10.1038/s41598-017-16168-w
81
MatouskovaK.JerryD. J.VandenbergL. N. (2020). Exposure to Low Doses of Oxybenzone during Perinatal Development Alters Mammary Gland Morphology in Male and Female Mice. Reprod. Toxicol.92, 66–77. 10.1016/j.reprotox.2019.08.002
82
MatozzoV.GagnéF.MarinM. G.RicciardiF.BlaiseC. (2008). Vitellogenin as a Biomarker of Exposure to Estrogenic Compounds in Aquatic Invertebrates: A Review. Environ. Int.34 (4), 531–545. 10.1016/j.envint.2007.09.008
83
MatusD. Q.MagieC. R.PangK.MartindaleM. Q.ThomsenG. H. (2008). The Hedgehog Gene Family of the Cnidarian, Nematostella vectensis, and Implications for Understanding Metazoan Hedgehog Pathway Evolution. Dev. Biol.313 (2), 501–518. 10.1016/j.ydbio.2007.09.032
84
MaxfieldF. R.IaeaD. B.PipaliaN. H. (2016). Role of STARD4 and NPC1 in Intracellular Sterol Transport. Biochem. Cel Biol.94 (6), 499–506. 10.1139/bcb-2015-0154
85
McKennaS. A.RichmondR. H.RoosG. (2001). Assessing the Effects of Sewage on Coral Reefs: Developing Techniques to Detect Stress before Coral Mortality. Bull. Mar. Sci.69 (2), 517–523.
86
MillerD. J.HemmrichG.BallE. E.HaywardD. C.KhalturinK.FunayamaN.et al (2007). The Innate Immune Repertoire in Cnidaria - Ancestral Complexity and Stochastic Gene Loss. Genome Biol.8 (4), R59. 10.1186/gb-2007-8-4-r59
87
MitchelmoreC. L.HeK.GonsiorM.HainE.HeyesA.ClarkC.et al (2019). Occurrence and Distribution of UV-Filters and Other Anthropogenic Contaminants in Coastal Surface Water, Sediment, and Coral Tissue from Hawaii. Sci. Total Environ.670, 398–410. 10.1016/j.scitotenv.2019.03.034
88
MohamedA. R.CumboV.HariiS.ShinzatoC.ChanC. X.RaganM. A.et al (2016). The Transcriptomic Response of the Coral Acropora Digitifera to a Competent Symbiodinium Strain: the Symbiosome as an Arrested Early Phagosome. Mol. Ecol.25 (13), 3127–3141. 10.1111/mec.13659
89
MorganM. B.ParkerC. C.RobinsonJ. W.PierceE. M. (2012). Using Representational Difference Analysis to Detect Changes in Transcript Expression of Aiptasia Genes after Laboratory Exposure to Lindane. Aquat. Toxicol.110-111, 66–73. 10.1016/j.aquatox.2012.01.001
90
MorganM. B.VogelienD. L.SnellT. W. (2001). Assessing Coral Stress Responses Using Molecular Biomarkers of Gene Transcription. Environ. Toxicol. Chem.20 (3), 5372–5543. 10.1897/1551-5028(2001)020<0537:acsrum>2.0.co10.1002/etc.5620200312
91
MorganM.GoodnerK.RossJ.PooleA. Z.SteppE.StuartC. H.et al (2015). Development and Application of Molecular Biomarkers for Characterizing Caribbean Yellow Band Disease inOrbicella Faveolata. PeerJ3, e1371. 10.7717/peerj.1371
92
OcampoI. D.Zárate-PotesA.PizarroV.RojasC. A.VeraN. E.CadavidL. F. (2015). The Immunotranscriptome of the Caribbean Reef-Building Coral Pseudodiploria Strigosa. Immunogenetics67 (9), 515–530. 10.1007/s00251-015-0854-1
93
PalmerC. V.Traylor-KnowlesN. (2012). Towards an Integrated Network of Coral Immune Mechanisms. Proc. R. Soc. B.279 (1745), 4106–4114. 10.1098/rspb.2012.1477
94
PandolfiJ. M.BradburyR. H.SalaE.HughesT. P.BjorndalK. A.CookeR. G.et al (2003). Global Trajectories of the Long-Term Decline of Coral Reef Ecosystems. Science301 (5635), 955–958. 10.1126/science.1085706
95
PastorianK.HawelL.ByusC. V. (2000). Optimization of cDNA Representational Difference Analysis for the Identification of Differentially Expressed mRNAs. Anal. Biochem.283 (1), 89–98. 10.1006/abio.2000.4622
96
PetersE. C.GassmanN. J.FirmanJ. C.RichmondR. H.PowerE. A. (1997). Ecotoxicology of Tropical Marine Ecosystems. Environ. Toxicol. Chem.16 (1), 12–40. 10.1002/etc.5620160103
97
PettersenE. F.GoddardT. D.HuangC. C.CouchG. S.GreenblattD. M.MengE. C.et al (2004). UCSF Chimera?A Visualization System for Exploratory Research and Analysis. J. Comput. Chem.25 (13), 1605–1612. 10.1002/jcc.20084
98
PutnamN. H.SrivastavaM.HellstenU.DirksB.ChapmanJ.SalamovA.et al (2007). Sea Anemone Genome Reveals Ancestral Eumetazoan Gene Repertoire and Genomic Organization. Science317 (5834), 86–94. 10.1126/science.1139158
99
QualmannB.KesselsM. M.TholeH. H.SierraltaW. D. (2000). A Hormone Pulse Induces Transient Changes in the Subcellular Distribution and Leads to a Lysosomal Accumulation of the Estradiol Receptor α in Target Tissues. Eur. J. Cel Biol.79 (6), 383–393. 10.1078/0171-9335-00061
100
ReitzelA. M.SullivanJ. C.Traylor-KnowlesN.FinnertyJ. R. (2008). Genomic Survey of Candidate Stress-Response Genes in the Estuarine AnemoneNematostella Vectensis. Biol. Bull.214 (3), 233–254. 10.2307/25470666
101
ReitzelA. M.TarrantA. M. (2009). Nuclear Receptor Complement of the Cnidarian Nematostella vectensis: Phylogenetic Relationships and Developmental Expression Patterns. BMC Evol. Biol.9, 230. 10.1186/1471-2148-9-230
102
RevelJ.MassiL.MehiriM.BoutouteM.MayzaudP.CapronL.et al (2016). Differential Distribution of Lipids in Epidermis, Gastrodermis and Hosted Symbiodinium in the Sea Anemone Anemonia Viridis. Comp. Biochem. Physiol. A: Mol. Integr. Physiol.191, 140–151. 10.1016/j.cbpa.2015.10.017
103
RotinenM.CelayJ.AlonsoM. M.ArrazolaA.EncioI.VillarJ. (2009). Estradiol Induces Type 8 17β-Hydroxysteroid Dehydrogenase Expression: Crosstalk between Estrogen Receptor α and C/EBPβ. J. Endrocrinol200 (1), 85–92. 10.1677/JOE-08-0134
104
RotjanR. D.SharpK. H.GauthierA. E.YeltonR.LopezE. M. B.CarilliJ.et al (2019). Patterns, Dynamics and Consequences of Microplastic Ingestion by the Temperate Coral, Astrangia Poculata. Proc. R. Soc. B.286, 20190726. 10.1098/rspb.2019.0726
105
SalahM.AbdelsamieA. S.FrotscherM. (2019). Inhibitors of 17β-Hydroxysteroid Dehydrogenase Type 1, 2 and 14: Structures, Biological Activities and Future Challenges. Mol. Cell Endocrinol.489, 66–81. 10.1016/j.mce.2018.10.001
106
SandersS. M.ShcheglovitovaM.CartwrightP. (2014). Differential Gene Expression between Functionally Specialized Polyps of the Colonial Hydrozoan Hydractinia Symbiolongicarpus (Phylum Cnidaria). BMC Genomics15 (1), 406. 10.1186/1471-2164-15-406
107
SchlesingerA.Kramarsky-WinterE.RosenfeldH.Armoza-ZvoloniR.LoyaY. (2010). Sexual Plasticity and Self-Fertilization in the Sea Anemone Aiptasia Diaphana. PLoS One5 (7), e11874. 10.1371/journal.pone.0011874
108
SheikhI. A.BegM. A. (2017). Endocrine Disruption: In Silico Interactions between Phthalate Plasticizers and Corticosteroid Binding Globulin. J. Appl. Toxicol.37 (12), 1471–1480. 10.1002/jat.3497
109
ShikinaS.ChangC.-F. (2016). “Sexual Reproduction in Stony Corals and Insight into the Evolution of Oogenesis in Cnidaria,” in The Cnidaria, Past, Present and Future (Switzerland: Springer), 249–268. 10.1007/978-3-319-31305-4_16
110
ShikinaS.ChenC.-J.ChungY.-J.ShaoZ.-F.LiouJ.-Y.TsengH.-P.et al (2013). Yolk Formation in a Stony Coral Euphyllia Ancora (Cnidaria, Anthozoa): Insight into the Evolution of Vitellogenesis in Nonbilaterian Animals. Endocrinology154 (9), 3447–3459. 10.1210/en.2013-1086
111
ShikinaS.ChungY.-J.ChiuY.-L.HuangY.-J.LeeY.-H.ChangC.-F. (2016). Molecular Cloning and Characterization of a Steroidogenic Enzyme, 17β-Hydroxysteroid Dehydrogenase Type 14, from the Stony Coral Euphyllia Ancora (Cnidaria, Anthozoa). Gen. Comp. Endocrinol.228, 95–104. 10.1016/j.ygcen.2016.02.006
112
SinghS. P.AzuaA.ChaudharyA.KhanS.WillettK. L.GardinaliP. R. (2010). Occurrence and Distribution of Steroids, Hormones and Selected Pharmaceuticals in South Florida Coastal Environments. Ecotoxicology19 (2), 338–350. 10.1007/s10646-009-0416-0
113
SinkovicsJ. (2015). The Cnidarian Origin of the Proto-Oncogenes NF-Κb/STAT and WNT-like Oncogenic Pathway Drives the Ctenophores (Review). Int. J. Oncol.47 (4), 1211–1229. 10.3892/ijo.2015.3102
114
SivikT.VikingssonS.GréenH.JanssonA. (2012). Expression Patterns of 17β-Hydroxysteroid Dehydrogenase 14 in Human Tissues. Horm. Metab. Res.44 (13), 949–956. 10.1055/s-0032-1321815
115
SlatteryM.HinesG. A.StarmerJ.PaulV. J. (1999). Chemical Signals in Gametogenesis, Spawning, and Larval Settlement and Defense of the Soft Coral Sinularia Polydactyla. Coral Reefs18 (1), 75–84. 10.1007/s003380050158
116
SnapeJ. R.MaundS. J.PickfordD. B.HutchinsonT. H. (2004). Ecotoxicogenomics: the challenge of Integrating Genomics into Aquatic and Terrestrial Ecotoxicology. Aquat. Toxicol.67 (2), 143–154. 10.1016/j.aquatox.2003.11.011
117
SnellT. W.BrogdonS. E.MorganM. B. (2003). Gene Expression Profiling in Ecotoxicology. Ecotoxicology12 (6), 475–483. 10.1023/B:ECTX.0000003033.09923.a8
118
SolbakkenJ. E.KnapA. H.OrrP. L. (1985). Uptake and Elimination of Lindane and a Phthalate Ester in Tropical Corals and Mussels. Mar. Environ. Res.16, 103–113. 10.1016/0141-1136(85)90012-1
119
StienD.SuzukiM.RodriguesA. M. S.YvinM.ClergeaudF.ThorelE.et al (2020). A Unique Approach to Monitor Stress in Coral Exposed to Emerging Pollutants. Sci. Rep.10 (1), 9601. 10.1038/s41598-020-66117-3
120
ŠtrukeljB.LenarčičB.GrudenK.PungerčarJ.RogeljB.TurkV.et al (2000). Equistatin, a Protease Inhibitor from the Sea Anemone Actinia equina, Is Composed of Three Structural and Functional Domains. Biochem. Biophysical Res. Commun.269 (3), 732–736. 10.1006/bbrc.2000.2356
121
StruttH.ThomasC.NakanoY.StarkD.NeaveB.TaylorA. M.et al (2001). Mutations in the Sterol-Sensing Domain of Patched Suggest a Role for Vesicular Trafficking in Smoothened Regulation. Curr. Biol.11 (8), 608–613. 10.1016/S0960-9822(01)00179-8
122
SutherlandK. P.PorterJ. W.TurnerJ. W.ThomasB. J.LooneyE. E.LunaT. P.et al (2010). Human Sewage Identified as Likely Source of white Pox Disease of the Threatened Caribbean elkhorn Coral, Acropora Palmata. Environ. Microbiol.12 (5), 1122–1131. 10.1111/j.1462-2920.2010.02152.x
123
TanB. L. L.HawkerD. W.MüllerJ. F.LeuschF. D. L.TremblayL. A.ChapmanH. F. (2007). Comprehensive Study of Endocrine Disrupting Compounds Using Grab and Passive Sampling at Selected Wastewater Treatment Plants in South East Queensland, Australia. Environ. Int.33 (5), 654–669. 10.1016/j.envint.2007.01.008
124
TarrantA.AtkinsonM.AtkinsonS. (2001). Uptake of Estrone from the Water Column by a Coral Community. Mar. Biol.139 (2), 321–325. 10.1007/s002270100556
125
TarrantA.AtkinsonM.AtkinsonS. (2004). Effects of Steroidal Estrogens on Coral Growth and Reproduction. Mar. Ecol. Prog. Ser.269, 121–129. 10.3354/meps269121
126
TarrantA. M.AtkinsonS.AtkinsonM. J. (1999). Estrone and Estradiol-17β Concentration in Tissue of the Scleractinian Coral, Montipora Verrucosa. Comp. Biochem. Physiol. Part A: Mol. Integr. Physiol.122 (1), 85–92. 10.1016/s1095-6433(98)10155-1
127
TarrantA. M.BlomquistC. H.LimaP. H.AtkinsonM. J.AtkinsonS. (2003). Metabolism of Estrogens and Androgens by Scleractinian Corals. Comp. Biochem. Physiol. Part B: Biochem. Mol. Biol.136 (3), 473–485. 10.1016/s1096-4959(03)00253-7
128
TarrantA. M. (2007). Hormonal Signaling in Cnidarians: Do We Understand the Pathways Well Enough to Know whether They Are Being Disrupted?Ecotoxicology16 (1), 5–13. 10.1007/s10646-006-0121-1
129
TarrantA. M.ReitzelA. M.BlomquistC. H.HallerF.TokarzJ.AdamskiJ. (2009). Steroid Metabolism in Cnidarians: Insights from Nematostella vectensis. Mol. Cell Endocrinol.301 (1-2), 27–36. 10.1016/j.mce.2008.09.037
130
TechnauU.RuddS.MaxwellP.GordonP. M. K.SainaM.GrassoL. C.et al (2005). Maintenance of Ancestral Complexity and Non-metazoan Genes in Two Basal Cnidarians. Trends Genet.21 (12), 633–639. 10.1016/j.tig.2005.09.007
131
ThibautR.PorteC. (2004). Effects of Endocrine Disrupters on Sex Steroid Synthesis and Metabolism Pathways in Fish. J. Steroid Biochem. Mol. Biol.92 (5), 485–494. 10.1016/j.jsbmb.2004.10.008
132
ThomasM. P.PotterB. V. L. (2013). The Structural Biology of Oestrogen Metabolism. J. Steroid Biochem. Mol. Biol.137, 27–49. 10.1016/j.jsbmb.2012.12.014
133
ThornH. A.QuinnJ. E.RoarkA. M. (2015). Reproductive and Developmental Effects of Tributyltin, Bisphenol A, and 17 β-estradiol in Pale Anemones (Aiptasia Pallida). Endocr. Disruptors3 (1), e1030062. 10.1080/23273747.2015.1030062
134
TijaniJ. O.FatobaO. O.PetrikL. F. (2013). A Review of Pharmaceuticals and Endocrine-Disrupting Compounds: Sources, Effects, Removal, and Detections. Water Air Soil Pollut.224 (11), 1–29. 10.1007/s11270-013-1770-3
135
TottaP.PesiriV.MarinoM.AcconciaF.MigliaccioA. (2014). Lysosomal Function Is Involved in 17β-Estradiol-Induced Estrogen Receptor α Degradation and Cell Proliferation. PLoS One9 (4), e94880. 10.1371/journal.pone.0094880
136
Traylor-KnowlesN.PalumbiS. R. (2014). Translational Environmental Biology: Cell Biology Informing Conservation. Trends Cel Biol.24 (5), 265–267. 10.1016/j.tcb.2014.03.001
137
TrottO.OlsonA. J. (2010). AutoDock Vina: Improving the Speed and Accuracy of Docking with a New Scoring Function, Efficient Optimization, and Multithreading. J. Comput. Chem.31 (2), 455–461. 10.1002/jcc.21334
138
TsuiM. M. P.LeungH. W.LamP. K. S.MurphyM. B. (2014). Seasonal Occurrence, Removal Efficiencies and Preliminary Risk Assessment of Multiple Classes of Organic UV Filters in Wastewater Treatment Plants. Water Res.53, 58–67. 10.1016/j.watres.2014.01.014
139
TwanW.-H.HwangJ.-S.ChangC.-F. (2003). Sex Steroids in Scleractinian Coral, Euphyllia Ancora: Implication in Mass Spawning1. Biol. Repro68 (6), 2255–2260. 10.1095/biolreprod.102.012450
140
TwanW.-H.HwangJ.-S.LeeY.-H.WuH.-F.TungY.-H.ChangC.-F. (2006). Hormones and Reproduction in Scleractinian Corals. Comp. Biochem. Physiol. Part A: Mol. Integr. Physiol.144 (3), 247–253. 10.1016/j.cbpa.2006.01.011
141
VanceJ. E. (2010). Transfer of Cholesterol by the NPC Team. Cell Metab12 (2), 105–106. 10.1016/j.cmet.2010.07.004
142
VerderameM.ScudieroR. (2017). Estrogen-dependent, Extrahepatic Synthesis of Vitellogenin in Male Vertebrates: A Mini-Review. C R. Biol.340 (3), 139–144. 10.1016/j.crvi.2017.01.005
143
Vidal‐DorschD. E.BayS. M.MaruyaK.SnyderS. A.TrenholmR. A.VanderfordB. J. (2012). Contaminants of Emerging Concern in Municipal Wastewater Effluents and marine Receiving Water. Environ. Tox. Chem.31 (12), 2674–2682. 10.1002/etc.2004
144
WearS. L.ThurberR. V. (2015). Sewage Pollution: Mitigation Is Key for Coral Reef Stewardship. Ann. N. Y Acad. Sci.1355 (1), 15–30. 10.1111/nyas.12785
145
Wood-CharlsonE. M.WeisV. M. (2009). The Diversity of C-type Lectins in the Genome of a Basal Metazoan, Nematostella vectensis. Dev. Comp. Immunol.33 (8), 881–889. 10.1016/j.dci.2009.01.008
146
WormB.LotzeH. K.JubinvilleI.WilcoxC.JeambeckJ. (2017). Plastic as a Persistent marine Pollutant. Ann. Rev. Environ. Resour.42, 1–26. 10.1146/annurev-environ-102016-060700
147
YoshiokaR. M.KimC. J.TracyA. M.MostR.HarvellC. D. (2016). Linking Sewage Pollution and Water Quality to Spatial Patterns of Porites Lobata Growth Anomalies in Puako, Hawaii. Mar. Pollut. Bull.104 (1-2), 313–321. 10.1016/j.marpolbul.2016.01.002
148
ZhangY.KalderonD. (2000). Regulation of Cell Proliferation and Patterning in Drosophila Oogenesis by Hedgehog Signaling. Development127 (10), 2165–2176. 10.1242/dev.127.10.2165
149
ZhangY.KalderonD. (2001). Hedgehog Acts as a Somatic Stem Cell Factor in the Drosophila Ovary. Nature410 (6828), 599–604. 10.1038/35069099
150
ZhongY.GuL. J.SunX. G.YangS. H.ZhangX. H. (2014). Comprehensive Analysis of Patched Domain-Containing Genes Reveals a Unique Evolutionary Pattern. Genet. Mol. Res.13 (3), 7318–7331. 10.4238/2014.February.13.11
151
ZhulynO.NieuwenhuisE.LiuY. C.AngersS.HuiC.-C. (2015). Ptch2 Shares Overlapping Functions with Ptch1 in Smo Regulation and Limb Development. Dev. Biol.397 (2), 191–202. 10.1016/j.ydbio.2014.10.023
Summary
Keywords
endocrine disruption chemicals, cnidaria, in silico modelling and docking, biomarkers, sex hormones, xenobiotics, transcriptional profiling, hedgehog signaling
Citation
Morgan MB, Ross J, Ellwanger J, Phrommala RM, Youngblood H, Qualley D and Williams J (2022) Sea Anemones Responding to Sex Hormones, Oxybenzone, and Benzyl Butyl Phthalate: Transcriptional Profiling and in Silico Modelling Provide Clues to Decipher Endocrine Disruption in Cnidarians. Front. Genet. 12:793306. doi: 10.3389/fgene.2021.793306
Received
25 October 2021
Accepted
24 December 2021
Published
11 January 2022
Volume
12 - 2021
Edited by
Douglas Mark Ruden, Wayne State University, United States
Reviewed by
Qing-Ping Xie, Zhejiang Academy of Agricultural Sciences, China
Bohu Pan, National Center for Toxicological Research (FDA), United States
Updates
Copyright
© 2022 Morgan, Ross, Ellwanger, Phrommala, Youngblood, Qualley and Williams.
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.
*Correspondence: Michael B. Morgan, mbmorgan@berry.edu
This article was submitted to Toxicogenomics, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.