Examining the Developmental Trajectory of an in Vitro Model of Mouse Primordial Germ Cells following Exposure to Environmentally Relevant Bisphenol A Levels

Background: Animal-based studies indicate that bisphenol A (BPA) exposure is detrimental to reproductive health, but its impact on the earliest stages of germ cell development remains poorly defined. Objectives: Using a murine in vitro model of early germ cell specification and differentiation, we sought to assess whether exposure to low levels of BPA prior to formation of primordial germ cells (PGCs) alters their differentiation trajectory and unique molecular program. Methods: We used an established method of in vitro differentiation of mouse embryonic stem cells (ESCs) into epiblast-like cells (EpiLCs) followed by PGC-like cells (PGCLCs), which together recapitulate defined stages of early germ cell development. Cellular consequences were determined using hemocytometer-based cell counting, fixation, and intracellular staining, followed by flow cytometry/fluorescence-activated cell sorting (FACS) of cells exposed to increasing concentrations (range: 1 nM–10 μM) of BPA. To interrogate and characterize gene expression differences resulting from BPA exposure, we also generated RNA-seq libraries from RNA extracted from FACS-purified PGCLCs and performed transcriptome analysis using bioinformatics-based approaches. Results: Exposure of EpiLCs to BPA resulted in higher numbers of cells that were associated with a higher proportion of cells in S-phase as well as a lower proportion undergoing apoptosis; this difference occurred in a concentration-dependent manner. Exposure also resulted in a greater fraction of EpiLCs showing signs of DNA damage. Remarkably, EpiLC exposure did not negatively affect PGC specification and resulted in a concentration-dependent effect on PGCLC proliferation in XX but not XY cells. PGCLC transcriptome analysis revealed an aberrant program with significant deregulation of X-linked genes and retrotransposon expression. Differential gene expression analysis also revealed the deregulation of genes associated with lipid metabolism as well as deregulated expression of genes associated with later stages of gametogenesis. Conclusions: To the best of our knowledge our findings represent the first characterization of the consequences of early BPA exposure on a model of mammalian PGC development, highlighting altered cell behavior, altered underlying pathways, and altered molecular processes. https://doi.org/10.1289/EHP8196


Introduction
Germ cells represent the bridge between generations, essential for the formation of new individuals in sexually reproducing organisms. It is therefore crucial that the information they contain is both correctly established, maintained, and executed over the course of development. Errors in these crucial developmental steps manifest as developmental defects, disease, and an overall reduction in reproductive fitness (Larose 2019).
In mammals, gametes are specified initially as primordial germ cells (PGCs), which in mice arise in embryonic day 6.5 (E6.5) postimplantation embryos, from cells derived in the epiblast layer (Gardner and Rossant 1979;Snow 1981). Over the course of their development, PGCs undergo a unique and complex reprogramming in both transcriptional and epigenetic signatures, necessary for the correct formation of germ cells and resulting gametes in the adult animal. Errors in reprogramming and proper germ cell function lie at the root of a number of adverse health outcomes, including chromosome anomalies, birth defects (Hassold et al. 1993;Hassold and Hunt 2001;Hunt and Hassold 2008), infertility (Matzuk and Lamb 2008), and cancer (Oosterhuis and Looijenga 2005) Perturbations in normal germ cell development are therefore something that must be minimized or, ideally, altogether avoided.
Germ cells at late stages of their development are extremely sensitive to environmental exposures, as highlighted by the breadth of chemicals exhibiting germ cell toxicity in both humans and rodents. These environmental agents include compounds as diverse as components of plastics (Hunt et al. 2003), pesticides (Harkonen 2005), cigarette smoke (Mailhes et al. 2000), antianxiety medication (Tanrikut et al. 2010), chemotherapeutic agents (Hales et al. 2005), and many more (Pacchierotti and Ranaldi 2006). However, there is a remarkable paucity of information concerning the impact of environmental agents on the earlier stages of germ cell development described above. Existing studies on the impact of toxicants on the window of germ cell development when they are potentially most vulnerable (i.e., when PGCs are present) are obfuscated by the difficulties intrinsic to their study: PGCs exist for a short developmental window between E6.5 and E11.5 in the developing embryo (Ginsburg 1990;Molyneaux 2001), and toxicant exposure typically covers multiple stages of germ cell development and differentiation; consequences of exposure can only be analyzed well after PGCs have differentiated into later developmental stages (e.g., mature oocytes or spermatozoa) or in resultant offspring derived from such gametes (reviewed in Matuszczak 2019).
In this context, the development of in vitro methods to differentiate embryonic stem cells (ESCs) into PGC-like cells (PGCLCs), such as those pioneered by Saitou et al., provide an invaluable system to study environmental impacts on early germ cells (Hayashi et al. 2011). Transcriptome and methylome analyses indicate that PGCLCs closely mirror the expression program and epigenetic reprogramming observed in in vivo PGCs and have the ability to give rise to live offspring (Hayashi et al. 2011;Shirane et al. 2016), indicating the recapitulation of all essential steps for gamete generation. Furthermore, because differentiation is ordered and controlled by the provision of specified growth factors in the culture media, the system is highly tractable, allowing spatial and temporal control of environmental exposure.
At moderate to high exposures, the plastic-manufacturing chemical bisphenol A (BPA) is a well-described reproductive toxicant in a wide variety of species, including nematodes (Allard and Colaiacovo 2010), fish (Chen et al. 2017), mouse (Hunt et al. 2003;Susiarjo et al. 2007), and humans (Brieño-Enríquez et al. 2011). Studies focused on the mammalian germline indicate detrimental consequences of BPA exposure on both male and female gametogenesis (Prins et al. 2019). However, owing to the complexity of gametogenesis, the delineation of critical phases for exposure, especially at the earlier stages of germ cell development, the capture of more subtle effects resulting from exposure to low environmentally relevant BPA levels, and the potential mechanisms involved all remain poorly defined.
Here, we leveraged the stepwise generation of PGCLCs to perform developmental stage-specific exposure. Specifically, we aimed to test whether exposure to low, environmentally relevant levels of BPA prior to initiation of the PGCLC program could alter their induction and/or developmental trajectory.

Differentiation Model and Exposure Paradigm
ES cells were maintained on 2i+LIF conditioned media for 2-3 d before being induced into EpiLCs by switching to EpiLC media ( Figures 1A and 2A). Because the transition from ESCs to EpiLCs itself involves profound molecular reprogramming in DNA methylation (reviewed in Wu and Zhang 2014), exposures were performed after that first reprogramming phase during a 24-h window prior to PGCLC induction and washed out when cells were switched to PGCLC culture medium. Because even low quantities of dimethylsulfoxide (DMSO) can alter ESC differentiation (Adler et al. 2006) and BPA is soluble in water at the concentration range tested (Shareef et al. 2006;Plahuta et al. 2015), BPA (239658; Sigma-Aldrich) was dissolved in autoclaved, double-distilled water (ddH 2 O) at 55°C for several hours at a concentration of 1 mM (equivalent to 0:228 mg=mL), before filtering and aliquoting, followed by storage at −20 C. The amount of BPA available in frozen stocks was verified by ELISA (BPA ELISA Kit; Detroit R & D, Inc.; see "BPA ELISA" section, Figure S1, and Excel Table S1). For each exposure experiment, aliquots of 1 mM BPA stocks were thawed and diluted in appropriate media (2i+LIF, EpiLC media, or PGCLC induction media). For water controls, the same volume of autoclaved ddH 2 O as the highest concentration of BPA dissolved in water being used was added to the culture media. Of the 4-order magnitude range of BPA concentrations tested, the 1-100 nM range was chosen to encompass the reported human BPA environmentally relevant range detected in serum, urine, and reproductive fluids (Vandenberg 2007), whereas the 1 lM and 10 lM concentrations represent moderate exposure levels comparable to those found in rodent in vivo models (Allard 2014).

BPA ELISA
ELISA was performed according to the manufacturer's protocol. Briefly, a six-point standard curve was set up, with BPA ranging in concentration between 1 × 10 6 pg=mL and 1 × 10 1 pg=mL, generated by 10-fold serial dilution using kit-provided BPA or dissolved BPA stock and sample dilution buffer. Samples were mixed with an equal volume (100 lL) of diluted BPA-HRP conjugate and incubated at room temperature for 2 h. Plates were then washed three times in 400 lL 1 × wash buffer before incubation with 200 lL TMB substrate at room temperature for 30 min. Development was stopped by addition of 50 lL of 2N sulfuric acid, and absorbance at 450 nM was read on a Tecan Infinite ® M1000 plate reader (Tecan). Corrected absorbance readings were then used to generate standard curves, according to the manufacturer's instructions.
"Intracellular Staining" section) and DAPI-incorporation in cells pulsed with BrdU (10 lM) 30-60 min prior to cell harvesting; and by proliferation tracing using cells loaded with CellTrace™ Yellow (5 lM, loaded for 20 min at 37°C, protected from light) (ThermoFisher Scientific). This dye passively diffuses through the cell membrane, where it remains covalently bound to intracellular amines, resulting in a stable, well-retained fluorescent signal with no/low cellular toxicity (Tempany 2018). Cell division results in a progressive diminution of signal intensity, which can be used to estimate the number of cell divisions a cell has undergone. To calculate numbers of cell division, the Proliferation Platform on FlowJo™ software was used (version 10; FlowJo, LLC).
in cell sort buffer [1 × Dulbecco's phosphate buffered saline (DPBS), 1% BSA, 1 mM EDTA, 25 mM HEPES], passed through a cell strainer (70 lm), and sorted on a BD FACSAria III (BD Biosciences), gating and collecting cell populations of interest in Eppendorf tubes containing GK15 media. Purified cells were used for the extraction of total RNA and subsequent library construction (see below). For analysis, cells were resuspended in fluorescence-activated cell sorting (FACS) buffer (1 × DPBS, 2% BSA, 1 mM EDTA, 25 mM HEPES, 100 U=mL penicillin, 0:1 mg=mL streptomycin, 2 mM L-glutamine), and analyzed on an LSR II (BD Biosciences). Cells were initially identified by forward-and side-scatter gating, with back-gating used to verify the  Table 2 and Excel Table S8, respectively. (C) FACS contour plots showing cell populations indicated from d5 aggregates. EpiLCs were treated with the conditions indicated for 24 h prior to aggregate formation. (D) Representative histograms showing CellTrace™ Yellow staining profile of the different cell populations indicated. The Y-axis represents the number of cells, whereas the X-axis represents the fluorescence intensity. Recurring cell divisions create the appearance of secondary peaks, which widens the original peak of undivided cells. (E-G) Scatter-bar plots showing proportion of cells in the different populations indicated that remain undivided or have undergone one, two, or three cell divisions in day 5 aggregates. Number of divisions calculated using FlowJo Proliferation tool (version 10; FlowJo, LLC). Scatter plots show mean ± standard deviation. n = 3-4. For calculated one-way ANOVA with Šidák adjusted p-values for comparison and cell proportions, see Table 3 and Excel Table S11, respectively. Note: ANOVA, analysis of variance; bFGF, basic fibroblast growth factor; BMP4, bone morphogenetic protein 4; BMP8b, bone morphogenetic protein 8b; BPA, bisphenol A; BVSC, B lymphocyte-induced maturation protein 1 Blimp1-mVenus and Stella-ECFP reporter transgenes; DAPI, 4 0 ,6-diamidino-2-phenylindole; DN, double negative; DP, double positive; EGF, epidermal growth factor; ESC, embryonic stem cell; EpiLC, epiblast-like cell; FACS, fluorescence-activated cell sorting; LIF, leukemia inhibitory factor; KSR, knockout serum replacement; PGC, primordial germ cell; PGCLC, PGC-like cell; SP, single positive. accuracy by which target cell populations were identified. Cell populations of interest were identified by 2-D plots displaying the parameter of interest. Manually defined gates as well as quadrants were used, as indicated. For positive controls to show the presence of DNA damage and apoptosis, v6.5 ES cells were treated with etoposide (10 lM; E1383; Sigma-Aldrich) and doxorubicin (500 nM; 44583; Sigma-Aldrich), both dissolved in DMSO, for 12 h prior to harvesting, fixing, and staining.

Intracellular Staining
Cells were permeabilized, fixed, and stained using the Apoptosis, DNA Damage and Cell Proliferation Kit (BD Biosciences) according to the manufacturer's protocol. In summary, cells were harvested and pelleted by centrifugation [5 min, 300 relative centrifugal force (rcf)]. Typically, between 2 × 10 5 and 5 × 10 5 cells were collected for each staining experiment. Pelleted cells were fixed in 100 lL BD Cytofix™/Cytoperm™ fixation/permeabilization solution before incubating on ice for 30 min. Cells were then washed with 100 lL 1 × BD Perm/Wash™ buffer (BD Biosciences) and pelleted by centrifugation. Following aspiration of supernatant, cells were resuspended in 100 lL BD Cytofix™/ Cytoperm™ Plus permeabilization buffer and incubated on ice for 10 min. Cells were then washed with 100 lL 1 × BD Perm/ Wash™ buffer (BD Biosciences) and pelleted by centrifugation. Following aspiration of supernatant, cells were resuspended in 100 lL BD Cytofix™/Cytoperm™ fixation/permeabilization solution before incubating on ice for 5 min. Cells were then washed with 100 lL 1 × BD Perm/Wash™ buffer and pelleted by centrifugation. Following aspiration of supernatant, cells were resuspended in 50 lL 1 × BD Perm/Wash™ buffer supplemented with 150 lg=mL DNase I and incubated at 37°C for 1 h. Following incubation, cells were washed in 150 lL 1 × BD Perm/Wash™ buffer and pelleted by centrifugation. Following aspiration of supernatant, cells were resuspended in 30 lL of 1 × BD Perm/Wash™ buffer containing antibodies [2 lL of each of anti-BrdU PerCP-Cy5.5 (clone 3D4); anti-!-H2AX AF 647 (clone N1-431); and anti-cleaved poly(ADP-ribose) polymerase (cPARP) (Asp214) PE (clone F21-852), for a final antibody dilution of 1:15]. All antibodies were derived from BD Biosciences as part of the Apoptosis, DNA Damage and Cell Proliferation Kit. Cells were incubated in the dark for 20 min, before washing with 180 lL 1 × BD Perm/Wash™ buffer and pelleting by centrifugation. Following aspiration of supernatant, cells were resuspended in 300 lL of 1 × BD Perm/Wash™ buffer with DAPI (final concentration 1 lg=mL DAPI). Cells were kept at 4°C and shielded from light until ready for analysis.

Total RNA Extraction
RNA was extracted from cell pellets (fresh or thawed snapfrozen) by using either an AllPrep DNA/RNA Mini or Micro Kit (Qiagen), according to the manufacturer's protocol. RNA concentration was measured using a NanoDrop™ 2000 UV spectrophotometer (ThermoFisher Scientific).

RNA-Seq Library Construction
Strand-specific RNA-seq libraries were prepared with Universal Plus mRNA-Seq kit (Nugen), according to the manufacturer's protocol. Briefly, this process consisted of poly(A) RNA selection, RNA fragmentation, and double-stranded cDNA generation using a mixture of random and oligo(dT) priming, followed by end repair to generate blunt ends, adaptor ligation, strand selection, and polymerase chain reaction (PCR) amplification to produce the final library. Different index adaptors were used for multiplexing samples in one sequencing lane. Sequencing was performed on Illumina HiSeq 3000 and NovaSeq 6000 sequencers for paired end (PE), 2 × 150 base pair (bp) runs. Data quality check was performed using Illumina Sequencing Analysis Viewer (SAV) software. Demultiplexing was performed with Illumina Bcl2fastq2 program (version 2.19.1.403; Illumina Inc.).

Expression Analysis of RNA-Seq Data
Fastq reads were checked for overall quality by FastQC (https:// www.bioinformatics.babraham.ac.uk/projects/fastqc/). For differential gene expression analysis, reads were aligned using STAR (Dobin et al. 2013) with the following commands: STAR -runThreadN 20 -readFilesCommand zcat -runMode alignReads -outFilterMultimapNmax 1 -outSAMtype BAM Unsorted -quantMode GeneCounts -twopassMode Basic. After sorting, indexing and converting into SAM file format using SAMtools (Li et al. 2009), read counts were obtained using HTSeq (https://htseq.readthedocs.io/en/master/) with the following commands: htseq-count − − mode = union − − stranded = no − − idattr = gene id -r pos -f sam. Output files were filtered to remove genes with nine or fewer read counts. These were then used for differential gene expression analysis using the edgeR Bioconductor package (https://bioconductor.org/packages/release/ bioc/html/edgeR.html). Gene counts were normalized using the trimmed mean of M-values normalization (TMM) method (Robinson and Oshlack, 2010) before determining counts per million (cpm) values. For a gene to be classified as showing differential expression between treated and untreated germ cells, a threshold fold-change of ≥2 and Benjamini-Hochberg adjusted p ≤ 0:05 had to be met.
For analysis of X-linked gene expression levels, genes were binned based on chromosomal location, before determining median expression levels for each chromosome. Median expression for all autosomal genes (all chromosomes except X) was determined and then used to normalize the median expression level for all chromosomes.
For analysis of repeat expression, fastq files were aligned using HiSat2 (Kim et al. 2019) with the -k parameter set to allow up to 100 alignments. Resultant sam files were then used to generate read counts for repeats using TETools (Lerat et al. 2017) and a mouse-specific repeat reference file generated using RepeatMasker. For each biological repeat between conditions (untreated and exposed), total numbers of reads following alignment using TETools was determined before down sampling of aligned sam files to ensure that an equivalent number of reads was analyzed across all samples. For an element to be considered differentially expressed between untreated and exposed samples, there needed to be 10 or more reads aligning to the repeat in either or both conditions, as well as a two-fold or greater difference in read number. Fastq files and Excel table showing

Gene Ontology (GO) Analysis and CirGO Plot Generation
Lists of differentially expressed genes as well as a list of all genes analyzed were generated from read counts using edgeR Bioconductor package . Enrichment of GO terms in lists of up-and down-regulated genes was determined using Gene Ontology enRIchment anaLysis and visuaLizAtion tool (GOrilla, http://cbl-gorilla.cs.technion.ac.il) (Eden, 2009). Redundant GO terms were removed using reduce + visualize gene ontology (REVIGO, http://revigo.irb.hr) (Supek et al. 2011) to generate .csv files. These were then used to generate Circular Gene Ontology (CirGO) plots using the open source CirGO software (version 1.0; https://github.com/IrinaVKuznetsova/CirGO.git) (Kuznetsova 2019). Terms were included if the Benjamini-Hochberg-adjusted p-value was less than 0.05.

Statistical Methods
For analysis of proliferation and cell staining data, calculated oneway analysis of variance (ANOVA) was initially performed, followed by Holm-Šidák testing for pairwise comparison between exposed and water/untreated control samples. For differential gene expression analysis, adjusted p-values were calculated using the Benjamini-Hochberg method. For statistical analysis of gene expression data derived from differential expression of RNA-seq data, Wilcoxon-signed rank test was used. In all cases, significance was determined by p-values less than or equal to 0.05. For analysis of repeat representation, chi-square analysis was performed. For each experiment, unless otherwise noted, technical n = 3.

Cell Proliferation and Viability Analysis of BPA-Exposed Murine EpiLCs
Because our exposure window corresponds to the stage of EpiLCs prior to PGCLC differentiation, we first examined the consequences of BPA exposure on EpiLCs. Exposure of these cells for 24 h resulted in a significantly higher numbers of cells recovered, which was in a concentration-dependent manner ( Figure 1B; p = 0:0086, one-way ANOVA, and Excel Table S2). By contrast, exposure of BVSC ESCs to BPA for up to 48 h had no significant effect on the number of cells recovered ( Figure S2; p = 0:62, one-way ANOVA, and Excel Table  S3). To determine the generality of BPA's effect on cell number in EpiLCs, we also examined two additional male (XY) ESC clones, v6.5 and BVSC-R8 ( Figure S3 and Excel Table S4). An interesting finding was that we did not detect any significant difference in cell number recovered when comparing unexposed to BPA-exposed cells over a range of different concentrations. This extended to the lack of perturbation in cell number in day 5 (d5) aggregates generated from BVSC-R8, ( Figure S4A and Excel Table S5). Transcriptome analysis revealed an absence of differentially expressed genes between BVSC-R8 PGCLCs derived from untreated and 100-nM BPAexposed EpiLCs ( Figure S4B and Excel Table S6). We found it interesting that d5 PGCLCs derived from BVSC-H18 and BVSC-R8 displayed a marked difference in their transcriptional identity, despite showing no significant difference in their expression of germ cell markers associated with the PGCLC state ( Figure S4B and S4C).
The higher number of BVSC-H18 cells recovered when exposed at the EpiLC stage could be accounted for by a higher rate of proliferation, lower apoptosis levels, or a combination of both. To investigate the mechanisms involved, we performed intracellular staining and analysis by flow cytometry. Because BPA exposure has previously been reported to result in DNA damage in germ cells (Barbonetti et al. 2016;Sahu et al. 2020;Yin et al. 2020), we also examined the proportion of cells staining positive for cH2AX, a marker of DNA damage (Rogakou et al. 1998). To control for the ability to detect cH2AX and cPARP signals, separate samples of EpiLCs were also exposed to the DNA-damaging agent etoposide (10 lM), the cytotoxin doxorubicin (500 nM), or both ( Figure S5). Analysis of BrdU and DAPI incorporation revealed that BPA exposure resulted in the presence of a significantly higher proportion of cells in S-phase/M, which occurred in a concentration-responsive manner ( Figure 1C and D; p = 0:0181, one-way ANOVA; Table 1 and Excel Table S7). Analysis of either BrdU with DAPI incorporation or detection of cleaved PARP (cPARP, a marker of early apoptosis), revealed a significantly higher proportion of cells undergoing apoptosis ( Figure 1C, E, and I; p = 0:0095, and 0.0002, respectively, one-way ANOVA; Table 1 and Excel Table S7). Examination of live events that were positive for cH2AX revealed the presence of a significantly higher proportion of cells with DNA damage in response to BPA exposure ( Figure 1C and F; p = 0:029, one-way ANOVA; Table 1 and Excel Table S7). We were surprised to observe a significantly higher proportion of both viable and proliferating cells with DNA damage ( Figure 1C, G, and H; p = 0:0015 and p = 0:02, respectively; Table 1 and Excel Table S7).

Proliferation Analysis of PGCLCs Derived from BPA-Exposed EpiLCs
We next investigated whether the brief exposure of EpiLCs to BPA had consequences on induction of PGCLCs in day 5 (d5) aggregates ( Figure 2A). During PGCLC induction, expression of the Blimp1-mVenus (BV) reporter is initially up-regulated before the Stella-ECFP (SC) reporter, allowing discrimination of three different cell populations: Blimp1-; Stella-(double negative, DN), Blimp1+; Stella-(single positive, SP), and finally PGCLCs that are Blimp1+; Stella+ (double positive, DP). In all induction experiments, equal numbers of control or exposed cells were    seeded at the start of each induction phase (see "Methods" section). Cell populations analyzed were initially identified using a "live" cell gate, based on forward-and side-scatter profiles. When measuring absolute numbers of events, BPA exposure of EpiLCs resulted in the presence of significantly higher numbers of SP and DP cells in d5 aggregates, ( Figure 2B and Table 2; p = 0:0084; p = 0:0071 for SP and DP, respectively, one-way ANOVA and Excel Table S8). As a proportion of all live events, significantly higher levels were observed for DNs and SPs ( Figure S6 and Excel Table S9 and S10; p = 0:0045 and 0.013, respectively, one-way ANOVA). For PGCLCs, higher numbers of cells were observed when derived from EpiLCs exposed to 10 nM BPA or above ( Figure 2B,C; Excel Table S8). In addition to the aforementioned methods of measuring proliferation, to further investigate the proliferative kinetics of the three different cell populations, we also used proliferation tracing of EpiLCs preloaded with CellTrace™ Yellow dye on day 0 of aggregate formation. This assay revealed that a greater proportion of DN cells had undergone three mitotic divisions, with concomitant reductions in the proportion of cells that had undergone either only one or two divisions ( Figure 2E; p = 0:0059, 0.0001, and 0.0003, calculated one-way ANOVA with Šidák-adjusted p-value, control vs. 100 nM, one, two, and three divisions, respectively; Table 3; Excel Table S11). A greater proportion of SP cells derived from EpiLCs exposed to BPA underwent three mitotic divisions, compared to control ( Figure 2C, D, and F; p = 0:013, one-way ANOVA; Table 3; Excel Table S11).
Although not significant, there was a similar trend for a higher proportion of DP/PGCLCs that had undergone three mitotic divisions ( Figure 2C, D, and G; Table 3; Excel Table S11).

Assessment of Gene Expression Differences in PGCLCs
Derived from EpiLCs Exposed to BPA Because differences in EpiLC and, subsequently, PGCLC behavior were observed at a concentration of 100 nM, we examined the nature of any underlying transcriptional alterations using RNAseq at this concentration ( Figure 3A). When compared to EpiLCs exposed to vehicle (water) control, unsupervised differential gene expression analysis of the transcriptome of EpiLCs exposed to 100 nM BPA revealed an absence of differentially expressed genes (Excel Table S12). However, as we observed that BPA exposure resulted in a higher proportion of EpiLCs labeling positive for cH2AX and a lower proportion of cells showing signs of apoptosis, we also performed a targeted analysis of the expression of genes in GO Biological Process-terms associated with either DNA damage or apoptosis to determine whether our transcriptome data could provide some insight into any potential mechanisms involved ( Figure S7 and Excel Table S13). For apoptosisassociated GO terms, nine GO terms showed significant differences between control and exposed samples ( Figure S7A and Excel Table S13). As well as higher median expression levels of GO terms associated with positive regulation of apoptosis (e.g., GO:0,043,065, positive regulation of apoptotic process), there was also higher expression of genes in one GO term associated with negative regulation of apoptosis (GO:2,001,237, negative regulation of extrinsic apoptotic signaling pathway). For GO terms associated with DNA damage, there was higher expression in 10 GO terms associated with checkpoint control of DNA damage ( Figure S7B and Excel Table S13). We next examined the transcriptome of PGCLCs derived from EpiLCs that were either untreated or exposed to 100 nM BPA. Expression profiling of known markers of PGCLC-identity revealed no significant differences between PGCLCs derived from control or 100 nM BPA-exposed EpiLCs ( Figure 3B; p = 0:1324, Wilcoxon matched pairs signed rank test). However, transcriptome analysis globally using principal component analysis (PCA) indicated a clear divergence between control and exposed cells ( Figure  3C). Analysis of differentially expressed genes using stringent conditions (fold-change ≥ 2; False Discovery Rate ðFDRÞ ≤ 0:05) revealed that although 161 genes were down-regulated ( Figure 3D; Excel Table S14), the majority of genes (501) were up-regulated ( Figure 3D; Excel Table S15). Examples of down-and upregulated genes are shown in Figure 3E.
To investigate biological functions/pathways affected in control and exposed germ cells, gene ontology analysis of enriched Biological Process terms was performed. Up-regulated genes were enriched in biological processes, such as ion and lipid transport [e.g., GO:0006811 (ion transport)]; GO:0006869 (lipid transport), lipoprotein metabolism [e.g., GO:0042157 (lipoprotein metabolic process)] and steroid metabolism [e.g., GO:0008202 (steroid metabolic process)] (Figure 4; Excel Table S16 and S17). We also noted that enriched terms in the GO Cellular Component category featured proteins associated with the extracellular space and cell membrane ( Figure S8).

BPA-Induced Perturbations in Gene Expression in PGCLCs
We determined whether genes associated with retinoic acid biosynthesis and/or signaling showed altered expression in PGCLCs derived from exposed EpiLCs. There was higher expression of genes in two gene ontology terms associated with retinoic acid biosynthesis: GO: 0042573 (retinoic acid metabolic process) and GO: 0042572 (retinol metabolic processes) ( Figure 5A,B; Figure S9A and 9B). Expression levels of genes containing a Retinoic Acid Response Element (RARE, Lalevée et al. 2011) and therefore inferred to be targets of retinoic acid signaling were found to be higher ( Figure 5C; Figure S9C, Wilcoxon-signed rank test p = 0:0022). An interesting finding was that although expression of Stra8 was not significantly different, there was a trend toward lower levels in exposed cells ( Figure S9D). Examination of whether BPA-exposure of EpiLCs had any impact on nuclear hormone receptor (NHR) expression levels in PGCLCs had revealed that neither retinoic acid nor estrogen nuclear receptor family members were significantly altered ( Figure S9E). Despite the lack of any differences in the expression of RA and estrogen nuclear receptor family members, differential gene expression analysis revealed two NHRs that were up-regulated: these were Hepatocyte nuclear factor 4 alpha (Hnf4a) and Nuclear receptor subfamily 2 group f member 2 (Nr2f2) (Excel Table S18). Table 3. Calculated p-values for differences in cell populations in PGCLCs derived from BPA-exposed EpiLCs.  Table 4; Excel Table S19 and Figure S10). However, other meiosis-specific GO terms revealed no significant differences resulting from BPA exposure ( Figure S11 and Excel Table S20).
Remarkably, among down-regulated genes, enrichment for transcripts associated with gametogenesis was observed (Excel Table S21). Several genes are members of the X-linked, lymphocyte regulated (Xlr) family, which generate coding and noncoding transcripts that are highly expressed in meiotic germ cells (Garchon et al. 1989;Shi et al. 2013). Of the 65 annotated Xlr family members listed in the latest build of the Mus musculus genome (GRCm38/mm10), 7 were found to be expressed and significantly lower in d5 PGCLCs derived from BPA-exposed EpiLCs ( Figure 5H; Wilcoxon-signed rank test p = 0:0005). In addition, we noted that a substantial proportion, 25% (40 of 161) of downregulated genes were X-linked, including the long noncoding RNA (lncRNA), X-inactive specific transcript (Xist) ( Figure 5I). We [Illustration in part created with ©BioRender -(biorender.com), per the Biorender terms and conditions.] (B) Log 10 -expression level-ranked heatmap of germ cell markers indicated in RNA-seq data generated from d5 aggregate FACSpurified DP/PGCLCs, derived from EpiLCs exposed to control (untreated) or 100 nM BPA. Candidate germ cell marker list taken from Hayashi et al. (Hayashi et al. 2011). (C) PCA of transcriptome of d5 DP/BVSC cells, derived from untreated or 100 nM BPA-exposed EpiLCs. (D) Volcano plot of RNA-seq data displaying gene expression pattern of d5 DP/PGCLCs, derived from untreated or 100 nM BPA-exposed EpiLCs. Significant differentially expressed genes (log 2-fold change ≥2, Benjamini-Hochberg-adjusted p ≤ 0:05) are highlighted in red. (E) Genome browser view of expression data for representative genes, up-(Ttr and Afp) and down-regulated (Lars2 and 1700013H16Rik) in d5 DP/BVSC cells, derived from 100 nM BPA-exposed EpiLCs. Y-axis values are normalized cpm/10 3 . For each locus, top panel represents "Untreated" and bottom panel represents " 100nM BPA exposed." For full list of differentially up-and down-regulated genes, see Excel Tables S14 and S15, respectively. Note: bFGF, basic fibroblast growth factor; BMP4, bone morphogenetic protein 4; BMP8b, bone morphogenetic protein 8b; BPA, bisphenol A; BVSC, B lymphocyte-induced maturation protein 1 Blimp1-mVenus and stella-ECFP reporter transgenes; cpm, counts per million; DN, double negative; DP, double positive; EGF, epidermal growth factor; embryonic stem cell; EpiLC, epiblast-like cell; FACS, fluorescence-activated cell sorting; KSR, knockout serum replacement; LIF, leukemia inhibitory factor; PGC, primordial germ cell; PGCLC, PGC-like cell; SCF, stem cell factor. also examined whether global levels of expression of X-linked genes relative to those on autosomes were affected. Comparison of expression values by chromosome, normalized by global median expression of all autosomal genes, in each sample revealed that BPA-exposed germ cells had significantly lower expression of X-linked genes ( Figure 5J; Benjamini-Hochberg-adjusted p < 0:05). In addition, this analysis revealed that BPA-exposure also resulted in lower expression of genes located on chromosomes 12 and higher expression of genes on chromosomes 1, 8, and 19 ( Figure 5J; Benjamini-Hochberg-adjusted p < 0:05). vi ta m in tr an sp or t tr a n sm e m b ra n e tr a n sp o rt o r g a n ic h y d r o x y c o m p o u n d t r a n s p o r t li p id t r a n s p o r t io n t r a n s p o r t an ion tra ns po rt GO1 transmembrane transport, 26.2% GO2 regulaƟon of lipid metabolic process, 18.1% GO3 system process, 14.4% GO4 regulaƟon of hormone levels, 13.5% GO5 regulaƟon of anatomical structure morphogenesis, 10.6% GO6 signal transducƟon, 4.3% GO7 organic hydroxy compound metabolic process, 3.8% plasma lipoprotein parƟcle remodeling, 3.0% mulƟcellular organismal process, 2.3% regulaƟon of microvillus organizaƟon, 2.3% regulaƟon of growth, 1.6% Figure 4. GO biological processes categories enriched in genes up-regulated in day 5 (d5) DP/ PGCLCs, derived from 100 nM BPA-exposed EpiLCs. CirGO plots showing parent representative records (inner ring) and child, description records (outer ring). For calculated p-values, see Excel Table S17. Note: BPA, bisphenol A; DP, double positive; EpiLC, epiblast-like cell; GO, gene ontology; PGC, primordial germ cell; PGCLC, PGC-like cell.  Table 4. GO biological process terms listed were obtained from http://www.informatics.jax.org/vocab/gene_ontology/GO:0,008,150. Genes within 10 kb of RAREs, and conserved between Mus musculus genome and at least six other different species, including humans were obtained from Lalevée et al. (Lalevée et al. 2011). For full lists of genes analyzed expression values and p-values, see Excel Tables S18-S21. (I) Chromosome ideogram indicating chromosome location of down-regulated genes in d5 PGCLCs derived from BPA-exposed EpiLCs. Down-regulated genes defined as those with a fold change of >2 and Benjamini-Hochberg-adjusted p < 0:05. Genes indicated by red triangles. (J) Scatter-bar plot showing autosome-normalized mean expression values for genes binned by chromosome. Median gene expression values for each chromosome was determined before normalizing by total autosomal gene expression value. Asterisks indicate chromosomes with differential expression values significantly different (Benjamini-Hochberg-adjusted p ≤ 0:05) between untreated and 100 nM BPA exposed samples. For all plots, the first in the series represents "Untreated" followed by "100 nM BPA exposed." Note: BPA, bisphenol A; EpiLCs, epiblast-like cell; GO, gene ontology; RARE, retinoic acid response element; PGC, primordial germ cell; PGCLC, PGC-like cell.
The continued enhanced proliferation observed in PGCLCs derived from BPA-exposed EpiLCs prompted us to examine whether there were any differences in expression levels of different factors listed in GO biological process terms associated with epigenetic processes. This examination revealed upregulation of genes in two terms: GO: 000012 (negative regulation of transcription by RNA polymerase II) ( Figure S12A) and GO: 1903319 (positive regulation of protein maturation) ( Figure S12B). Interestingly, genes in GO terms associated with proteins involved in DNA methylation (GO: 0003606), histone methylation (GO: 0016571), and piRNA metabolic process (GO: 0034587) (Figure S12C-E) did not show significant differences in expression between control and treated germ cells.

Analysis of Repeat Expression in PGCLCs Derived from BPA-Exposed EpiLCs
We examined whether there were any perturbations in repeat expression in untreated and exposed germ cells. We analyzed read counts mapping to all annotated mouse repeats listed on RepeatMasker (http://repeatmasker.org), focusing on those repeats that showed differential expression in untreated in comparison with exposed germ cells. This analysis revealed that Figure 6. Analysis of repeat DNA expression in PGCLCs derived from EpiLCs exposed to 100 nM BPA. (A) Pie chart showing percentage distribution of different repeat families up-regulated in exposed germ cells. For full list of repeat families, see Excel Table S22. (B) Violin plots showing aggregate expression levels of different repeat families indicated, up-regulated in exposed germ cells. Dashed lines represent median and quartiles. Wilcoxon matched-pairs signed rank test p-values: <0:0001 (ERVK, ERV1, ERVL and L1); <0:001 for ERVL-MaR. First plot represents "Untreated" followed by "100 nM exposed." (C) Bar graph comparing observed with expected frequency of different repeat families indicated, up-regulated in exposed germ cells. Expected frequencies of each family calculated based on total number of upregulated repeat families and proportion occurrence in RepeatMasker database. First plot represents "Observed," followed by "Expected" frequency. For the full list of up-regulated elements and repeat superfamilies, see Excel Tables S23 and S24, respectively. Note: BPA, bisphenol A; EpiLCs, epiblast-like cell; PGC, primordial germ cell; PGCLC, PGC-like cell. Table 4. Calculated Wilcoxon paired sign rank test p-values of expression levels for different gene sets in PGCLCs derived from BPA-exposed EpiLCs.
although 10 repeat elements were down-regulated in PGCLCs derived from BPA-exposed EpiLCs (Excel Table S22), the vast majority (180 repeat elements) were up-regulated ( Figure 6A,B; Excel Table S23, S24). More than three-quarters of the upregulated elements belonged to five groups: ERVK, ERV1, ERVL, L1, and ERVL-MaLR. Using the proportion of all different repeat groups annotated in RepeatMasker, chi-square analysis was performed to determine if any groups displayed over-or underrepresentation ( Figure 6C). This analysis revealed that aside from ERVL-MaLR, there was overrepresentation of ERVK, ERV1, ERVL, and L1 group elements (v 2 = 0:025). Because these groups contain young and presumably active retrotransposons, this finding supports the notion that these elements might therefore be coordinately regulated and therefore transcribed. However, our analysis of the RNA-seq data did not reveal significant differences in transcript levels of well-known restriction factors ( Figure S13).

Discussion
Although a large body of animal-based studies exists on the reproductive impacts of BPA, to the best of our knowledge the data presented in this study are the first to reveal the specific impact of BPA-exposure in a model of the earliest stages of mammalian germ cell development. We found that EpiLCs displayed greater sensitivity toward BPA-exposure compared with ESCs, with exposed cells showing higher proliferation, lower apoptosis, and a higher proportion of cells staining positive for cH2AX, a well-established marker of DNA damage. This higher proliferation had consequences on PGCLC formation, resulting in a higher number of PGCLCs formed in day 5 germ cell aggregates. Although expression analysis of known germ cell markers clearly showed that these cells retained a germ cell identity, they harbored altered expression of genes involved in ion and lipid transport, lipoprotein and steroid metabolism, and elevated retrotransposon expression. We were surprised to find that this behavior was not observed in a different cell line (BVSC-R8). Because BVSC-H18 is genetically female (XX), a possible explanation for this difference is the existence of sex-specific differences in the cell lines analyzed and their responsiveness to BPA. Female PGCLCs showed lower expression of X-linked genes, indicating that BPA exposure affected normal dosage compensation. Our results therefore provide evidence for the sensitivity of the PGC differentiation program and homeostasis to environmentally relevant levels of BPA. BPA exposure of EpiLCs for a period of 24 h was sufficient to result in a higher number of cells recovered, in a concentrationdependent manner ( Figure 1B). Intracellular staining followed by FACS indicated that this was correlated with a higher proportion of cells in S-phase as well as a lower proportion of apoptotic cells ( Figure 1C-I). Collectively, these results show that the higher number of cells recovered following BPA exposure was in part a result of higher cell proliferation and lower apoptosis, despite signs of DNA damage in a higher proportion of cells. The effect on cell numbers in vertebrates is not unprecedented: A study in zebrafish involving BPA exposure of male embryos for 24 h post fertilization, during which time PGCs migrate and colonize the gonadal ridge, revealed that although exposure resulted in a reduction in the number of PGCs, there was no consequence on final testes development or fertility in the adult (Lombo et al. 2019), suggesting that BPA exposure results in "hyperproliferation" of the fewer founding PGCs and/or derivative cells as a compensatory mechanism in the future gonad.
Unsupervised differential gene expression analysis did not uncover significant transcriptional differences between control and BPA-exposed EpiLCs; this lack of significant difference might be attributed to the fact that cells were exposed for only 24 h, resulting in insufficient time to detect widespread differences in gene expression. However, our targeted analysis of genes associated with DNA damage indicated higher expression of genes involved in DNA damage checkpoint control, specifically expression of factors involved in detecting the presence of DNA damage. Although cH2AX is a molecular marker of DNA damage (Rogakou et al. 1998), it is also prevalent in highly proliferative cells, such as ES cells. In addition to this, it has been proposed that the presence of cH2AX in ES cells is a consequence of global chromatin decondensation (Banáth et al. 2009). Although we cannot exclude the possibility that BPA exposure results in a further increase in chromatin de-condensation, the fact that our data clearly showed higher proliferation as well higher expression of genes associated with DNA damage checkpoint control supports the notion that BPA is responsible for an elevation in DNA damage; the exact mechanisms, whether direct, indirect or both, remain to be determined.
The molecular underpinning responsible for enhanced proliferation remains unclear. Transcriptome analysis of untreated or 100 nM BPA-exposed EpiLCs did not reveal significant differences in gene expression (Excel Table S12). Despite this finding, there was evidence suggestive of epigenetic consequences on day 5 PGCLCs (DP). Although no longer exposed, these cells also showed higher proliferation ( Figure 2B,G). We note that because nonPGCLCs (DN) and presumed transitioning cells (SP) also showed higher proliferation, this proliferation is likely a general effect on EpiLCs and resultant cell types, rather than a consequence specific to PGCLCs. An intriguing observation was the higher proportion of cells showing signs of DNA damage ( Figure  1F-H). Previous data from a testicular coculture model indicated that BPA had no effect on the proportion of cells with signs of DNA damage when cultured for up to 72 h and at BPA concentrations up to 100 lM (Yin et al. 2020). Consistent with this, 1-wk exposure of spermatogonial stem cell (SSC) cultures to 10 lM BPA followed by testicular transplantation did not result in significant alterations in the presence of double-strand breaks in pachytene spermatocytes (Karmakar et al. 2017). Our results highlight differences in sensitivity based on developmental stage when exposure arises and suggests this aspect deserves greater attention when considering chemical toxicity on basic cell biological processes. Of particular concern, the presence of more cells with DNA damage has implications for the transmission of genetic lesions/mutations into later stages of germ cell development, which could have consequences for disease or the overall fitness of any resultant offspring.
The observation that BPA exposure had no discernible effect on PGCLCs derived from BVSC-R8 EpiLCs was unexpected but, in retrospect, perhaps not as surprising as it initially seems. Although we cannot rule out clonal differences between the BVSC-R8 and BVSC-H18 cell lines, we note that a key difference between these two cell lines is sex. Sex-specific gene expression differences in early embryonic cells, including ESCs and EpiLCs, have previously been reported (Lowe et al. 2015;Shirane et al. 2016), which presumably could affect responsiveness to BPA exposure. This hypothesis is supported by our finding that BPA exposure led to a significant deregulation of X-linked transcripts. In mammals, dosage of X-linked genes is regulated by X inactivation (Brockdorff et al. 1991), whereby expression of Xist coats and inactivates one of the two X chromosomes. Both X chromosomes are reportedly active in female EpiLCs (Hayashi et al. 2012) and in day 6 PGCLCs (Shirane et al. 2016). In light of this, our observation that PGCLCs derived from BPA-exposed BVSC-H18/female EpiLCs showed a disproportionate impact on X-linked genes suggests a disturbance in correct dosage regulation of genes located on the X chromosome and warrants further investigation. Taken together with our results, these observations suggest the existence of sex-specific differences in the response of early germ cells to environmental exposures, prior to commitment to sexual fate.
Gene ontology analysis of transcriptome data from d5 PGCLCs derived from BPA-exposed EpiLCs indicated the enrichment of biological processes involved in lipid transport, lipoprotein metabolism and steroid metabolism (Figure 4; Excel Tables S16, S17). Our previous studies in C. elegans revealed a novel mechanism by which BPA interferes with steroidogenic acute regulatory (StAR) protein mediated cholesterol transport, resulting in reduced fertility ). Up-regulation of genes involved in cholesterol biosynthesis supports the notion of a cholesterol imbalance in developing PGCLCs that might affect their normal function. In addition, the fact that pathways involved in cholesterol metabolism are also perturbed in this murine system underscores the high conservation of the pathways/processes affected by BPA in nematodes and mammals.
PGCs represent the immediate precursor stage before sexual differentiation. In the case of male gametogenesis, cells enter mitotic arrest and differentiate into prospermatogonia/gonocytes, whereas in females, cells enter meiosis I. Because PGCLCs derived from BPA-exposed EpiLCs exhibited transcriptional differences, we used a candidate approach to examine whether there were any effects on the expression levels of genes involved in later stages of gametogenesis ( Figure 5). The mechanism by which meiosis is initiated depends on retinoic acid signaling, which induces, among other factors, Stra8 expression (Baltus et al. 2006;Bowles et al. 2006). Our analyses suggest that although we did not find evidence of precocious meiotic entry, BPA exposure nonetheless resulted in transcriptional perturbations in GO categories implicated in RA signaling and meiotic regulation. Our inability to corroborate previous work suggesting a specific effect on RA signaling could be explained by methodological reasons. In contrast to our use of a defined protocol to specifically derive murine germ cells in the form of PGCLCs, Aoki and Takada investigated the effect of BPA exposure on embryoid bodies, which results in the formation of cells derived from all three germ layers (Aoki and Takada 2012).
Our transcriptomics analysis in PGCLCs does not unambiguously support a mechanism by which BPA may be acting on EpiLCs (i.e., binding to/activation of specific NHRs/pathways), and the effects on hormone signaling could be indirect and involve multiple pathways. Interestingly, our finding that Hnf4a transcript levels were higher recapitulates a previous study that showed that this gene is similarly perturbed by BPA in several tissues (Shu et al. 2019), suggesting this gene may be a general target deregulated by BPA, irrespective of cell type. Owing to the perturbations we detected in gene expression, our results suggest that BPA exposure may result in the formation of germ cells that are compromised in their function. For example, expression of alpha-fetoprotein (Afp) is substantially elevated in PGCLCs derived from BPA-exposed EpiLCs ( Figure 3E; Excel Table S15). Through its capacity to bind estrogens at high affinity and thereby compete with cognate receptors (Uriel et al. 1976;Savu et al. 1981;Nishi et al. 1991) Afp is proposed to have antiestrogenic effects. Assuming this translates to protein expression, BPA-induced higher Afp expression level is therefore likely to have a negative rather than positive impact on estrogen signaling in PGCLCs. In utero exposure of mice to BPA during the formation of primary follicles has demonstrated an impact on germ cell nest breakdown, crucial in determining the size of the primary follicle and future fertility (Wang et al. 2014). Because this process is triggered by a drop in estrogen-mediated signaling, considering our data, it will be interesting to determine whether estrogen signaling is affected in PGCLCs derived from BPA-exposed EpiLCs that are subsequently used to generate primary oocytes. A modification of the PGCLC protocol has been successfully used to generate functional oocytes in culture (Hayashi et al. 2017), and the use of this expanded protocol in future experiments could test the functional consequences of BPA-induced impacts on aspects of fertility such as primary follicle size.
A vital aspect of PGC differentiation is the need to silence repetitive elements, principally retrotransposons. The absence of mechanisms correctly regulating their expression has consequences on later fertility. For example, in E13.5 PGCs, endogenous retroviruses (ERVs) are repressed by the histone methyltransferase SET domain bifurcated histone lysine methyltransferase 1 (Setdb1), responsible for trimethylation of lysine position 9 of histone H3 (H3K9Me3) (Liu et al. 2014). Although germline conditional Setdb1 mutant females do not show reduced numbers of PGCs, at postnatal day 10 (P10) these animals have smaller ovaries, with gonadal hypotrophy persisting through puberty and into adulthood (Liu et al. 2014). Our analysis of d5 PGCLC transcriptome data indicated that BPA exposure has an impact on retrotransposon regulation. We note that the libraries generated were enriched for polyadenylated mRNA. Given that the vast majority of repeats are not full-length but truncated and therefore lack a poly A tail (Lanciano and Cristofari 2020), it is likely that a large proportion of the repeat compartment, was underrepresented in our data set. Correct repression of repeat elements, primarily RNA-based retrotransposons, is a crucial feature of normal germ cell development. Following global reductions in genome-wide levels of 5mC in CG contexts, and alterations in posttranslational histone modifications (including H3K27me3 and H3K9me2/3) retrotransposons are activated; future studies examining the epigenome of PGCLCs derived from BPA-exposed EpiLCs should provide insight into whether repeat compartments harbor correct patterns of modifications and levels. Ordinarily, L1 elements are thought to be active between E13.5 to E16.5 in male and female germ cells (Seisenberger et al. 2012;Kobayashi et al. 2013). The higher expression of repeats induced by BPA is consistent with disruption in the normal mechanisms involved in their repression. Aside from histone methylation, the mammalian genome has evolved numerous host restriction mechanisms to control repeat expression at both the transcriptional and posttranscriptional levels (Goodier 2016). The nature of the mechanism(s) responsible for the elevation observed in repeat transcripts remains unclear. Although expression analysis of epigenetic factors and known host-encoded restriction factors are unchanged at the RNA level ( Figure S13), it remains possible that prior exposure to BPA could alter protein levels, subcellular localization, and/or catalytic activity of one or more factors involved in regulating repeat expression. Recently in the context of epigenetic mechanisms, interest has centered on the importance of phase-separated condensates with liquidlike properties in the cell nucleus (reviewed by Bergeron-Sandoval et al. 2016). Germ cells, including PGCs, have such a structure termed "nuage," which consists of (among other features) the chromatoid body (CB) (Kotaja and Sassone-Corsi 2007), itself enriched in repetitive RNA and RNA-binding proteins involved in repeat restriction such as Ddx4/Mvh (Meikar et al. 2014) and Maelstrom (Soper et al. 2008). Future experiments involving high-resolution imaging and using the PGCLC system could provide insight into whether the structure and/or composition of CB components is altered in cells derived from BPA-exposed EpiLCs. This approach could therefore provide a mechanism for the observed elevation of repeat transcripts.
In summary, to the best of our knowledge this study represents the first description of the consequences of toxicant exposure at a low, environmentally relevant concentration on the earliest stages of germ cell development. Together, our observations suggest that despite BPA-exposure increasing PGCLC formation, these cells are abnormal, and we posit that both their genetic and epigenetic content carry perturbations that are likely to affect the function of derivative gametes. Future studies will examine the nature of these consequences.