Endocrine-disrupting potentials of equine estrogens equilin, equilenin, and their metabolites, in the medakaOryzias latipes:in silicoand DNA microarray studies

In vivo and in silico analyses of estrogenic potential of equine estrogens in medaka (Oryzias latipes)

Equine estrogens (EEs) are widely used in hormone replacement therapy pharmaceuticals for postmenopausal women. Previous studies have shown that EEs occur in the aquatic environment; however, the potential estrogenicity and risk of EEs in aquatic organisms, including fish, have yet to be studied in detail. Therefore, we evaluated the estrogenic potential of major EEs, namely equilin (Eq), 17α-dihydroequilin (17α-Eq), 17β-dihydroequilin (17β-Eq), equilenin (Eqn), 17α-dihydroequilenin (17α-Eqn), and 17β-dihydroequilenin (17β-Eqn), on medaka (Oryzias latipes) using in vivo and in silico assays. Quantitative real-time RT-PCR analyses revealed that expression levels of choriogenin L (ChgL) and choriogenin H (ChgH) in medaka embryos responded to various types and concentrations of EEs in a concentration-dependent manner, whereas transcription levels of vitellogenin 1 were not significantly affected by any of the EEs in the concentration range tested. The order of the in vivo estrogenic potencies of EEs was as follows: 17β-Eq > Eq > 17β-Eqn > Eqn > 17α-Eqn > 17α-Eq. Additionally, the 50% effective concentrations (EC₅₀) of 17β-Eq was lower than that of 17β-estradiol. We also investigated the interaction potential of EEs with medaka estrogen receptor (ER) subtypes in silico using a three-dimensional model of the ligand-binding domain (LBD) for each ER and docking simulations. All six EEs were found to interact with the LBDs of ERα, ERβ1, and ERβ2. The order of the in silico interaction potentials of EEs with each ER LBD was as follows: 17β-Eq > 17α-Eq > Eq > 17β-Eqn > 17α-Eqn > Eqn. Furthermore, we identified the key amino acids that interact with EEs in each ER LBD; our findings suggest that amino acids and/or their hydrogen bonding may be responsible for the ligand-specific interactions with each ER. This study is the first to comprehensively analyze the estrogenic potential of EEs in medaka both in vivo and in silico.

Potential mechanisms underlying embryonic developmental toxicity caused by benzo[a]pyrene in Japanese medaka (Oryzias latipes)

Polycyclic aromatic hydrocarbons (PAHs), such as benzo[a]pyrene (BaP), are widely distributed in air, water, and sediments; however, limited data are available regarding their potential adverse effects on the early life stages of fish. In this study, we evaluated the embryonic teratogenicity and developmental toxicity of BaP in Japanese medaka (Oryzias latipes) using a nanosecond pulsed electric field (nsPEF) technique and predicted their molecular mechanisms via transcriptome analysis. The gas chromatography/mass spectrometry analyses revealed that the BaP was efficiently incorporated into the embryos by nsPEF treatment. The embryos incorporating BaP presented typical teratogenic and developmental effects, such as cardiovascular abnormalities, developmental abnormalities, and curvature of backbone. DNA microarray analysis revealed several unique upregulated genes, such as those involved in cardiovascular diseases, various cellular processes, and neural development. Furthermore, the gene set enrichment and network analyses found several genes and hub proteins involved in the developmental effects of BaP on the embryos. These findings suggest a potential mechanism of teratogenicity and developmental toxicity caused by exposure to BaP. The nsPEF and transcriptome analyses in combination can be effective for evaluating the potential effects of chemical substances on medaka embryos.

Occurrence and seasonal variation of equine estrogens, equilin and equilenin, in the river water of Japan: Implication with endocrine-disrupting potentials to Japanese medaka (Oryzias latipes)

In this study, we determined the concentration of equine estrogens, such as equilin (Eq) and equilenin (Eqn), in the river water collected from nine research stations in Hokkaido, Japan. The LC-MS/MS analysis revealed that Eq concentrations were 2.7 ± 6.7, 0.22 ± 0.12, and 1.2 ± 0.64 ng/L in Sep 2015, Feb 2016, and Jul 2016, respectively. Eqn had concentration levels similar to those of Eq. Comparison of the concentrations at nine research stations showed that seasonal variation was observed in the detected Eq and Eqn concentration levels. This study was the first to show the occurrences and seasonal variation of Eq and Eqn in the river water of Japan. We further investigated the reproductive and transgenerational effects of Eq in Japanese medaka (Oryzias latipes) exposed to 10, 100, and 1000 ng/L for 21 days and assessed the transcriptional profiles of the estrogen-responsive genes in the livers of both sexes. The reproduction assay demonstrated that 1000 ng/L of Eq adversely affected the reproduction (i.e. fecundity) in the F₀ generation and that the hatching of F₁ generation fertilized eggs was reduced in the 100 and 1000 ng/L treatment groups. Our qRT-PCR assay revealed that the mRNA expression levels of hepatic vitellogenin 1 and 2, choriogenin L and H, and estrogen receptor α were significantly up-regulated in males exposed to 100 and/or 1000 ng/L of Eq. In contrast, the transcriptional levels of several genes, such as pregnane X receptor and cytochrome P450 3A, were down-regulated in the livers of males after the 21-d exposure. These results suggest that Eq has endocrine-disrupting potential such as reproductive and transgenerational effects by the modulation of hepatic estrogen-responsive genes expression on medaka.

Nanosecond pulsed electric field incorporation technique to predict molecular mechanisms of teratogenicity and developmental toxicity of estradiol-17β on medaka embryos

Herein, we propose using a nanosecond pulsed electric field (nsPEF) technique to assess teratogenicity and embryonic developmental toxicity of estradiol-17β (E₂ ) and predict the molecular mechanisms of teratogenicity and embryonic developmental defects caused by E₂ on medaka (Oryzias latipes). The 5 hour post-fertilization embryos were exposed to co-treatment with 10 μm E₂ and nsPEF for 2 hours and then continuously cultured under non-E₂ and nsPEF conditions until hatching. Results documented that the time to hatching of embryos was significantly delayed in comparison to the control group and that typical abnormal embryo development, such as the delay of blood vessel formation, was observed. For DNA microarray analysis, 6 day post-fertilization embryos that had been continuously cultured under the non-E₂ and nsPEF condition after 2 hour co-treatments were used. DNA microarray analysis identified 542 upregulated genes and one downregulated gene in the 6 day post-fertilization embryos. Furthermore, bioinformatic analyses using differentially expressed genes revealed that E₂ exposure affected various gene ontology terms, such as response to hormone stimulus. The network analysis also documented that the estrogen receptor α in the mitogen-activated protein kinase signaling pathway may be involved in regulating several transcription factors, such as FOX, AKT1 and epidermal growth factor receptor. These results suggest that our nsPEF technique is a powerful tool for assessing teratogenicity and embryonic developmental toxicity of E₂ and predict their molecular mechanisms in medaka embryos.

Neofunctionalization of Androgen Receptor by Gain-of-Function Mutations in Teleost Fish Lineage

Steroid hormone receptor family provides an example of evolution of diverse transcription factors through whole-genome duplication (WGD). However, little is known about how their functions have been evolved after the duplication. Teleosts present a good model to investigate an accurate evolutionary history of protein function after WGD, because a teleost-specific WGD (TSGD) resulted in a variety of duplicated genes in modern fishes. This study focused on the evolution of androgen receptor (AR) gene, as two distinct paralogs, ARα and ARβ, have evolved in teleost lineage after TSGD. ARα showed a unique intracellular localization with a higher transactivation response than that of ARβ. Using site-directed mutagenesis and computational prediction of protein-ligand interactions, we identified two key substitutions generating a new functionality of euteleost ARα. The substitution in the hinge region contributes to the unique intracellular localization of ARα. The substitution on helices 10/11 in the ligand-binding domain possibly modulates hydrogen bonds that stabilize the receptor-ligand complex leading to the higher transactivation response of ARα. These substitutions were conserved in Acanthomorpha (spiny-rayed fish) ARαs, but not in an earlier branching lineage among teleosts, Japanese eel. Insertion of these substitutions into ARs from Japanese eel recapitulates the evolutionary novelty of euteleost ARα. These findings together indicate that the substitutions generating a new functionality of teleost ARα were fixed in teleost genome after the divergence of the Elopomorpha lineage. Our findings provide a molecular explanation for an adaptation process leading to generation of the hyperactive AR subtype after TSGD.

In vivo and in silico analyses of estrogenic potential of bisphenol analogs in medaka (Oryzias latipes) and common carp (Cyprinus carpio)

Various studies have demonstrated the estrogenic effect of bisphenol A (BPA), a member of bisphenol analogs (BPs), in in vitro and in vivo assays. However, limited data are available on the estrogenic potentials and risks of other BPs in aquatic organisms. In addition, the estrogenic effect of chemicals is known to have species-specific responses in teleost fish. The objective of this study was to evaluate the potential estrogenic effects of BPs on the medaka (Oryzias latipes) and common carp (Cyprinus carpio) using in vivo and in silico assays. Our quantitative real-time PCR analyses revealed that the expression levels of several hepatic estrogen-responsive biomarker genes in male medaka responded to various types and concentrations of BPs in a dose-response manner. The order of in vivo estrogenic potencies of BPs was as follows: BPC≈BPAF>BPB>BPA⋙BPP. To further investigate the interaction potential of BPs with medaka estrogen receptor α (ERα) in silico, a three-dimensional model of the ERα ligand-binding domain (LBD) was built and docking simulations were performed. The docking simulation analysis revealed that BPC interaction potential for medaka ERα LBD was the most potent, followed by BPAF and BPA. Comparing this with carp ERα LBD revealed that the interaction potentials of these BPs to medaka ERα LBD were more stable than to carp ERα LBD. Furthermore, we identified key amino acid residues in medaka ERα LBD that interacted with BPC (Glu356, Arg397, and Cys533), BPAF (Thr350 and Glu356), and BPA (Glu356 and Met424), and found some differences in these key amino acid residues between medaka and carp ERα LBDs. These results of in vivo and in silico analyses showed potential estrogenic effects of BPs in teleost fish, and they also indicated that the differences in interaction potentials and key amino acid residues between medaka and carp ERα LBDs may be due to the differences between the species and estrogenic potencies of the selected BPs.