Metabolic signatures of bisphenol A and genistein in Atlantic salmon liver cells

Probing the bioconcentration and metabolism disruption of bisphenol A and its analogues in adult female zebrafish from integrated AutoQSAR and metabolomics studies

Plenty of emerging bisphenol A (BPA) substitutes rise to wait for assessment of bioconcentration and metabolism disruption. Computational methods are useful to fill the data gap in chemical risk assessment, such as automated quantitative structure-activity relationship (AutoQSAR). It is not clear how AutoQSAR performs in predicting the bioconcentration factor (BCF) in adult zebrafish. Herein, AutoQSAR was used to predict the logBCFs of BPA, bisphenol AF (BPAF), bisphenol B, bisphenol F and bisphenol S (BPS). For the test set, a linear relationship was shown between the observed and predicted logBCFs with a slope of 0.97. The predicted logBCFs of these five bisphenols were quite close to their experimental data with a slope of 0.94, suggesting better performance than directed message passing neural networks and EPI Suite with a slope of 0.69 and 0.61, respectively. Thus, AutoQSAR is powerful in modeling logBCFs in fish with minimal time and expertise. To link bioconcentration with metabolic effects, female zebrafish were exposed to BPA, BPAF and BPS for metabolomics analysis. BPA caused a significant disturbance in amino acid metabolism, while BPAF and BPS significantly altered another three metabolic pathways, showing chemical-specific responses. BPAF with the highest logBCF elicited the strongest metabolomic responses reflected by the metabolic effect level index, followed by BPA and BPS. Thus, BPAF and BPS elicited higher or similar metabolism disruption compared with BPA in female zebrafish, respectively, reflecting consequences of bioconcentration.

Flunitrazepam and its metabolites exposure disturb the zebrafish gut-liver axis: Combined microbiome and metabolomic analysis

Due to clinical treatment and illegal use, psychoactive substances have been widely detected in the aquatic environment. In this study, we investigated the effects of the benzodiazepine drug flunitrazepam (FLZ) and its metabolite 7-aminoflunitrazepam (7-FLZ) on the gut-liver axis of zebrafish. Zebrafish were exposed to two concentrations of FLZ and 7-FLZ (0.05 and 1 μg/L) for 30 days. Results showed that both FLZ and 7-FLZ exposure altered the relative abundance of Proteobacteria at the phylum level, with significant differences observed at the genus level for pathogenic bacteria such as Paracoccus, Shewanella, and Aeromonas. Metabolomics results showed both exposures significantly interfered with nucleotide and amino acid metabolism. The imbalance of gut microbiota and metabolic disorder increased the level of malondialdehyde, which in turn heightened the permeability of the gut mucosal barrier. FLZ and 7-FLZ induced oxidative stress in the liver via the gut-liver axis, leading to decreased levels of glucose, total cholesterol, and triglyceride, as well as the down-regulation of glycolipid metabolism-related genes (PPARα, PPARγ, FABP2, Fabp11, PFKFB3, and LDHA). Metabolomics results revealed that FLZ and 7-FLZ significantly affected the biosynthesis of amino acids and arginine, and other metabolic pathways such as nucleotide, nicotinate and nicotinamide, and purine in the liver. Our results unveiled the mechanisms behind the toxicological effects of psychoactive substances on the gut-liver axis, providing valuable data for ecological and environmental risk assessments.

Hepatotoxic effect of dietary phytoestrogens on juvenile cultured Russian sturgeon (Acipenser gueldenstaedtii)

In the last two decades, much controversy has grown over the use of soybean products in aquafeeds, especially for carnivorous fish like sturgeons. One point of discussion is the effect of soybean phytoestrogens on fish health. There are many aspects of phytoestrogen utilization in aquafeeds, therefore, the aim of this study is to verify if common legume phytoestrogens can affect juvenile cultured sturgeon erythrocyte and hepatocyte genotoxicity and cause liver pathology. Russian sturgeons were fed from 100 till 365 dph¹ with daidzein, genistein, and coumestrol supplemented diets in concentrations: 10, 0.05 and 0.001 g kg^-1 of feed, respectively. The SCGE² method combined with qPCR of three genes involved in DNA repair and genome maintenance, namely cyp1a1, gaad45a and p53 were analyzed. The results were compared with histopathological evaluation of liver tissue. In fish fed with coumestrol supplemented diet, DNA strand damage was the highest in both erythrocytes and hepatocytes, however, simultaneously the lowest level of oxidative DNA damage was found. Additionally, slightly elevated expression of the p53 gene was observed along with a decreased number of apoptotic hepatocytes, which suggests that low concentration of coumestrol may support DNA repair mechanisms in the liver. Although, daidzein showed a preventive effect only against fibrosis. Isoflavones did not show a significant effect on DNA damage in studied cells. Genistein was found to increase macro- and microvesicular steatosis, portal hepatitis and fibrosis, indicating its negative role in the development of liver injuries. Daidzein alleviated some sturgeon liver damage, especially macrovesicular steatosis and interface hepatitis. However, it increased hepatocyte apoptosis, which may suggest daidzein potentially inducing liver injury, though not manifested by other histopathological lesions. Therefore, it can be concluded that at given concentrations, the tested phytoestrogens did not show clearly hepatoprotective effect in sturgeons.

Epigenetic analytical approaches in ecotoxicological aquatic research

Environmental epigenetics has become a key research focus in global climate change studies and environmental pollutant investigations impacting aquatic ecosystems. Specifically, triggered by environmental stress conditions, intergenerational DNA methylation changes contribute to biological adaptive responses and survival of organisms to increase their tolerance towards these conditions. To critically review epigenetic analytical approaches in ecotoxicological aquatic research, we evaluated 78 publications reported over the past five years (2016-2021) that applied these methods to investigate the responses of aquatic organisms to environmental changes and pollution. The results show that DNA methylation appears to be the most robust epigenetic regulatory mark studied in aquatic animals. As such, multiple DNA methylation analysis methods have been developed in aquatic organisms, including enzyme restriction digestion-based and methyl-specific immunoprecipitation methods, and bisulfite (in)dependent sequencing strategies. In contrast, only a handful of aquatic studies, i.e. about 15%, have been focusing on histone variants and post-translational modifications due to the lack of species-specific affinity based immunological reagents, such as specific antibodies for chromatin immunoprecipitation applications. Similarly, ncRNA regulation remains as the least popular method used in the field of environmental epigenetics. Insights into the opportunities and challenges of the DNA methylation and histone variant analysis methods as well as decreasing costs of next generation sequencing approaches suggest that large-scale epigenetic environmental studies in model and non-model organisms will soon become available in the near future. Moreover, antibody-dependent and independent methods, such as mass spectrometry-based methods, can be used as an alternative epigenetic approach to characterize global changes of chromatin histone modifications in future aquatic research. Finally, a systematic guide for DNA methylation and histone variant methods is offered for ecotoxicological aquatic researchers to select the most relevant epigenetic analytical approach in their research.

Bisphenol F induces liver-gut alteration in zebrafish

The unease of consumers with bisphenol A has led to the increased industrial usage of bisphenol F (BPF), which is a new hazard to environmental health. Here, zebrafish were exposed to three BPF concentrations (0.5, 5, and 50 μg/L) from the embryonic stage for 180 days. Results showed that zebrafish body length and weight decreased and hepatosomatic index values increased, even at environmentally relevant concentration. Histological analysis identified the occurrence of hepatic fibrosis and steatosis in 5 and 50 μg/L groups, which indicated the liver injury caused by BPF. Based on the untargeted metabolomics results, a dose-dependent variation in the effects of BPF on liver metabolism was found, and amino acids, purines and one carbon metabolism were the main affected processes in the 0.5, 5, and 50 μg/L treatments, respectively. At the same time, BPF induced a shift in intestinal microbiome composition, including decreased abundance of Erysipelotrichaceae, Rhodobacteraceae and Gemmobacter. In addition, the correlation analysis suggested an association between gut microbiome changes and affected hepatic metabolites after BPF exposure. These findings indicate that a liver-gut alteration is induced by long-term BPF exposure.

Metabolomics coupled with network pharmacology study on the protective effect of Keguan-1 granules in LPS-induced acute lung injury

CONTEXT

Keguan-1 (KG-1) plays a vital role in enhancing the curative effects, improving quality of life, and reducing the development of acute lung injury (ALI).

OBJECTIVE

To unravel the protective effect and underlying mechanism of KG-1 against ALI.

MATERIALS AND METHODS

C57BL/6J mice were intratracheally instilled with lipopolysaccharide to establish the ALI model. Then, mice in the KG-1 group received a dose of 5.04 g/kg for 12 h. The levels of proinflammatory cytokines, chemokines, and pathological characteristics were determined to explore the effects of KG-1. Next, untargeted metabolomics was used to identify the differential metabolites and involved pathways for KG-1 anti-ALI. Network pharmacology was carried out to predict the putative active components and drug targets of KG-1 anti-ALI.

RESULTS

KG-1 significantly improved the levels of TNF-α (from 2295.92 ± 529.87 pg/mL to 1167.64 ± 318.91 pg/mL), IL-6 (from 4688.80 ± 481.68 pg/mL to 3604.43 ± 382.00 pg/mL), CXCL1 (from 4361.76 ± 505.73 pg/mL to 2981.04 ± 526.18 pg/mL), CXCL2 (from 5034.09 ± 809.28 pg/mL to 2980.30 ± 747.63 pg/mL), and impaired lung histological damage. Untargeted metabolomics revealed that KG-1 significantly regulated 12 different metabolites, which mainly related to lipid, amino acid, and vitamin metabolism. Network pharmacology showed that KG-1 exhibited anti-ALI effects through 17 potentially active components acting on seven putative drug targets to regulate four metabolites.

DISCUSSION AND CONCLUSIONS

This work elucidated the therapeutic effect and underlying mechanism by which KG-1 protects against ALI from the view of the metabolome, thus providing a scientific basis for the usage of KG-1.

2'-Fucosyllactose Remits Colitis-Induced Liver Oxygen Stress through the Gut-Liver-Metabolites Axis

Liver oxygen stress is one of the main extraintestinal manifestations of colitis and 5% of cases develop into a further liver injury and metabolic disease. 2′-fucosyllactose (2′-FL), a main member of human milk oligosaccharides (HMOs), has been found to exert efficient impacts on remitting colitis. However, whether 2′-FL exerts the function to alleviate colitis-induced liver injury and how 2′-FL influences the metabolism via regulating gut microbiota remain unknown. Herein, in our study, liver oxygen stress was measured by measuring liver weight and oxygen-stress-related indicators. Then, 16S full-length sequencing analysis and non-target metabolome in feces were performed to evaluate the overall responses of metabolites and intestinal bacteria after being treated with 2′-FL (400 mg/kg b.w.) in colitis mice. The results showed that, compared with the control group, the liver weight of colitis mice was significantly decreased by 18.30% (p < 0.05). After 2′-FL treatment, the liver weight was significantly increased by 12.65% compared with colitis mice (p < 0.05). Meanwhile, they exhibited higher levels of oxidation in liver tissue with decreasing total antioxidant capacity (T-AOC) (decreased by 17.15%) and glutathione (GSH) levels (dropped by 22.68%) and an increasing malondialdehyde (MDA) level (increased by 36.24%), and 2′-FL treatment could reverse those tendencies. Full-length 16S rRNA sequencing revealed that there were 39 species/genera differentially enriched in the control, dextran sulphate sodium (DSS), and DSS + 2′-FL groups. After treatment with 2′-FL, the intestinal metabolic patterns, especially glycometabolism and the lipid-metabolism-related process, in DSS mice were strikingly altered with 33 metabolites significantly down-regulated and 26 metabolites up-regulated. Further analysis found DSS induced a 40.01%, 41.12%, 43.81%, and 39.86% decline in acetic acid, propionic acid, butyric acid, and total short chain fatty acids (SCFAs) in colitis mice (all p < 0.05), respectively, while these were up-regulated to different degrees in the DSS + 2′-FL group. By co-analyzing the data of gut microbiota and metabolites, glycometabolism and lipid-metabolism-associated metabolites exhibited strong positive/negative relationships with Akkermansia_muciniphila (all p < 0.01) and Paraprevotella spp. (all p < 0.01), suggesting that the two species might play crucial roles in the process of 2′-FL alleviating colitis-induced liver oxygen stress. In conclusion, in the gut−liver−microbiotas axis, 2′-FL mediated in glucose and lipid-related metabolism and alleviated liver oxygen stress via regulating gut microbiota in the DSS-induced colitis model. The above results provide a new perspective to understand the probiotic function of 2′-FL.

LC/MS/MS-Based Liver Metabolomics to Identify Chronic Liver Injury Biomarkers Following Exposure to Arsenic in Rats

Arsenic is a widespread natural metalloid element. Long-term chronic exposure to arsenic can lead to different degrees of liver injury. Although the etiology of this disease is well known, to date, the underlying mechanism of arsenic-induced liver injury remains unclear, and no specific treatment exists because of the complexity of arsenic. In the present study, potential biomarkers and metabolic pathways in the livers of Wistar rats treated with arsenic for 24 weeks were investigated using an integrated metabolic approach with an LC-Orbitrap Q Exactive™ HF-X mass spectrometer. Markedly increased liver levels of arsenic, alanine aminotransferase (ALT), alkaline phosphatase (ALP), and total bile acid (TBA) were detected in the arsenic treatment groups (P < 0.05). Furthermore, histopathological examination of liver tissues showed obviously swollen, loose cytoplasm and increased necrosis in the arsenic treatment groups compared with those in the control group (P < 0.05). Metabonomics results showed that 109 metabolites (variable importance in the projection (VIP) > 1; fold change > 2 or < 0.5; P adjusted < 0.05) changed significantly after exposure to arsenic and included 71 upregulated metabolites and 38 downregulated metabolites. Additionally, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that 6 metabolic pathways with statistical significance-phenylalanine metabolism, pyruvate metabolism, glycolysis/gluconeogenesis, citrate cycle (TCA cycle), thiamine metabolism, and vitamin B6 metabolism-were selected, and 13 differential metabolites were detected to be involved in regulating these metabolic pathways. The present study could help identify potential biomarkers and their functions, as well as metabolic pathways, likely providing evidence for the early diagnosis, prevention, and mechanistic study of arsenism.

Metabolomic Analysis Reveals the Mechanisms of Hepatotoxicity Induced by Aflatoxin M1 and Ochratoxin A

Aflatoxin M1 (AFM1) is the only toxin with the maximum residue limit in milk, and ochratoxin A (OTA) represents a common toxin in cereals foods. It is common to find the co-occurrence of these two toxins in the environment. However, the interactive effect of these toxins on hepatoxicity and underlying mechanisms is still unclear. The liver and serum metabolomics in mice exposed to individual AFM1 at 3.5 mg/kg b.w., OTA at 3.5 mg/kg b.w., and their combination for 35 days were conducted based on the UPLC-MS method in the present study. Subsequent metabolome on human hepatocellular liver carcinoma (Hep G2) cells was conducted to narrow down the key metabolites. The phenotypic results on liver weight and serum indicators, such as total bilirubin and glutamyltransferase, showed that the combined toxins had more serious adverse effects than an individual one, indicating that the combined AFM1 and OTA displayed synergistic effects on liver damage. Through the metabolic analysis in liver and serum, we found that (i) a synergistic effect was exerted in the combined toxins, because the number of differentially expressed metabolites on combination treatment was higher than the individual toxins, (ii) OTA played a dominant role in the hepatoxicity induced by the combination of AFM1, and OTA and (iii) lysophosphatidylcholines (LysoPCs), more especially, LysoPC (16:1), were identified as the metabolites most affected by AFM1 and OTA. These findings provided a new insight for identifying the potential biomarkers for the hepatoxicity of AFM1 and OTA.

The effects of bisphenol A on liver proteome and mucus vitellogenin in comparison to plasma as a non-invasive biomarker in immature Siberian sturgeons (Acipenser baerii)

This study was done to evaluate the effects of Bisphenol A (BPA) on Siberian sturgeon (Acipenser baerii). As liver is the main organ in the homeostatic adjustments to stress, we used a proteomics method to address molecular response in this tissue. Also, we compared the levels of vitellogenin in plasma and mucus to propose that the last one be a non-invasive method to analyze this biomarker. The fish received 1, 10, and 100 μg g^-1 week^-1 BPA intraperitoneally for two weeks. The samples were taken on days 0, 7, and 14. Plasma vitellogenin level increased as the highest value was recorded in the group with 100 μg g^-1 week^-1 of BPA. Changes in the mucus and blood vitellogenin showed a similar pattern, suggesting that mucus could be used for evaluating the changes in blood vitellogenin. Comparative proteomics was used to determine the proteome of the liver of A. baerii in the highest dose of BPA in comparison with the control. Sixteen proteins were identified that their expression changed at least twice between the studied groups. The proteomic results showed that BPA increased the expression of proteins involved in the detoxification and metabolism, activated glycolysis, and produced necrosis in the liver.

Bisphenol A and Bisphenol S Induce Endocrine and Chromosomal Alterations in Brown Trout

Bisphenol A is a widely used compound found in large amount of consumer products. As concerns have been raised about its toxicological and public health effect, the use of alternatives to bisphenol A are now increasing. Bisphenol S is one of the analogues being used as a replacement for bisphenol A despite the fact that little is known about the effects of bisphenol S on living organisms. In this study, we investigated the potential endocrine and genotoxic effects of bisphenol A and bisphenol S in juvenile brown trout (Salmo trutta). The fish were exposed to the compounds for either 2 weeks or 8 weeks via sustained-release cholesterol implants containing doses of 2 mg/kg fish or 20 mg/kg fish of the substances. The effects on the thyroid hormone levels and the estrogenic disrupting marker vitellogenin were evaluated, along with the genotoxic markers micronucleated cells and erythrocyte nuclear abnormalities. An increase in plasma vitellogenin was observed in fish exposed to the high dose of bisphenol A for 2 weeks. At this experimental time the level of the thyroid hormone triiodothyronine (T3) in plasma was elevated after bisphenol S exposure at the high concentration, and paralleled by an increase of micronucleated cells. Moreover, bisphenol A induced an increase of micronuclei frequency in fish erythrocytes after the exposure at the lowest dose tested. Taken together the results indicate that both bisphenol A and its alternative bisphenol S cause endocrine disrupting and genotoxic effects in brown trout, although suggesting two different mechanisms of damage underlying bisphenol A and bisphenol S activity.

Genistein Induces Adipogenic and Autophagic Effects in Rainbow Trout (Oncorhynchus mykiss) Adipose Tissue: In Vitro and In Vivo Models

Soybeans are one of the most used alternative dietary ingredients in aquafeeds. However, they contain phytoestrogens like genistein (GE), which can have an impact on fish metabolism and health. This study aimed to investigate the in vitro and in vivo effects of GE on lipid metabolism, apoptosis, and autophagy in rainbow trout (Oncorhynchus mykiss). Primary cultured preadipocytes were incubated with GE at different concentrations, 10 or 100 μM, and 1 μM 17β-estradiol (E2). Furthermore, juveniles received an intraperitoneal injection of GE at 5 or 50 µg/g body weight, or E2 at 5 µg/g. In vitro, GE 100 μM increased lipid accumulation and reduced cell viability, apparently involving an autophagic process, indicated by the higher LC3-II protein levels, and higher lc3b and cathepsin d transcript levels achieved after GE 10 μM. In vivo, GE 50 µg/g upregulated the gene expression of fatty acid synthase (fas) and glyceraldehyde-3-phosphate dehydrogenase in adipose tissue, suggesting enhanced lipogenesis, whereas it increased hormone-sensitive lipase in liver, indicating a lipolytic response. Besides, autophagy-related genes increased in the tissues analyzed mainly after GE 50 µg/g treatment. Overall, these findings suggest that an elevated GE administration could lead to impaired adipocyte viability and lipid metabolism dysregulation in rainbow trout.

Identification of multiple dysregulated metabolic pathways by GC-MS-based profiling of liver tissue in mice with OVA-induced asthma exposed to PM2.5

Particulate matter (PM) exposure increases the risk of asthma. However, the effect of PM_2.5 exposure on liver metabolism in mice with asthma symptoms remains unclear. We established an ovalbumin (OVA)-induced asthma model in mice and divided the animals into four groups: control group (C), PM_2.5 exposure group (P), OVA-induced asthma group (O) and OVA-induced asthma PM_2.5 exposure group (OP). Gas chromatography-mass spectrometry (GC-MS) was used to identify the metabolite markers and related perturbed metabolic pathways in mouse liver tissue after PM_2.5 exposure. Multivariate analysis showed 9 and 12 potential metabolite markers in the P and OP groups, respectively, after PM_2.5 exposure that were significantly correlated with lipid peroxidation indices. PM_2.5 exposure perturbed 5 and 7 metabolic pathways in the P and OP groups, respectively. These metabolic pathways mainly involve the lipid metabolism, amino acid metabolism, carbohydrate metabolism, and nucleotide metabolism. These results highlight the potential to study PM_2.5-triggered alterations via liver tissue in normal and OVA-induced asthmatic mice to gain a more realistic appraisal of the resulting early toxicity events. Additionally, these results revealed potential metabolite markers of early antioxidant defense events triggered by PM_2.5 and indicated that metabolite markers are more sensitive than antioxidant indicators.

Identification of multiple dysregulated metabolic pathways by GC-MS-based profiling of lung tissue in mice with PM2.5-induced asthma

The risk of development of asthma, a multi-faceted chronic disease, increases as a result of exposure to PM_2.5. However, the mechanism underlying asthma-related metabolic changes caused by PM_2.5 exposure is unclear. Here, we investigated the major metabolic changes, metabolic pathways involved, and underlying molecular mechanisms in mice with PM_2.5 exposure-induced asthma. Forty-eight adult female mice were randomly assigned to control (C), low concentration-PM_2.5 exposure: 0.50 mg kg^-1 (L), medium concentration-PM_2.5 exposure: 1.58 mg kg^-1 (M), and high concentration-PM_2.5 exposure: 4.98 mg kg^-1 (H) groups. M and H groups presented significantly higher IL-4, IL-8, IL-1β, IL-5, IL-13, and OVA-specific IgE levels, and significantly lower IFN-γ levels, than the C group, as well as significantly increased eosinophil count and MUC5AC expression in the lung tissue. These findings indicate that exposure to medium and high concentrations of PM_2.5 induced asthma in mice. Statistical analyses identified 13 asthma-related major metabolites, which were analyzed by gas chromatography-mass spectrometry (GC-MS). Meta Mapp Software revealed 4 major metabolic pathways. PM_2.5-induced ATP requirement and oxidative stress may perturb metabolic processes in asthma. The present findings increase our understanding of the toxic effect of PM_2.5 in the development of asthma and identify potentially useful biomarkers.

Prenatal Bisphenol A Exposure in Mice Induces Multitissue Multiomics Disruptions Linking to Cardiometabolic Disorders

The health impacts of endocrine-disrupting chemicals (EDCs) remain debated, and their tissue and molecular targets are poorly understood. In this study, we leveraged systems biology approaches to assess the target tissues, molecular pathways, and gene regulatory networks associated with prenatal exposure to the model EDC bisphenol A (BPA). Prenatal BPA exposure at 5 mg/kg/d, a dose below most reported no-observed-adverse-effect levels, led to tens to thousands of transcriptomic and methylomic alterations in the adipose, hypothalamus, and liver tissues in male offspring in mice, with cross-tissue perturbations in lipid metabolism as well as tissue-specific alterations in histone subunits, glucose metabolism, and extracellular matrix. Network modeling prioritized main molecular targets of BPA, including Pparg, Hnf4a, Esr1, Srebf1, and Fasn as well as numerous less studied targets such as Cyp51 and long noncoding RNAs across tissues, Fa2h in hypothalamus, and Nfya in adipose tissue. Lastly, integrative analyses identified the association of BPA molecular signatures with cardiometabolic phenotypes in mouse and human. Our multitissue, multiomics investigation provides strong evidence that BPA perturbs diverse molecular networks in central and peripheral tissues and offers insights into the molecular targets that link BPA to human cardiometabolic disorders.

Molecular characterisation of oestrogen receptor ERα and the effects of bisphenol A on its expression during sexual development in the Chinese giant salamander (Andrias davidianus)

The aim of this study was to characterise the molecular structure of the oestrogen receptor ERα and to evaluate the effect of bisphenol A (BPA) on ERα expression during sexual development of the Chinese giant salamander (Andrias davidianus). The ERα cDNA of A. davidianus includes an open reading frame of 1755bp (encoding 584 amino acids), a 219-bp 5' untranslated region (UTR) and a 611-bp 3'UTR. A polyadenylation signal was not found in the 3'UTR. Amino acid sequence analysis showed high homology between ERα of A. davidianus and that of other amphibians, such as Andrias japonicas (99.66% identity) and Rana rugose (81.06% identity). In 3-year-old A. davidianus, highest ERα expression was observed in the liver and gonads. During different developmental stages in A. davidianus (from 1 to 3 years of age), ERα expression in the testes increased gradually. ERα was localised in the epithelial cells of seminiferous lobules and in interstitial cells. ERα-positive cells were more abundant in the interstitial tissue during testicular development. ERα was located in the nucleus of oocytes during ovary development. We found that the sex of 6-month-old A. davidianus larvae could not be distinguished anatomically. The sex ratio did not change after larvae were treated with 10μM BPA for 1 month. However, BPA treatment reduced bodyweight and ERα expression in the gonads in male larvae.