Study population and sample collection
The participants were selected from the Guangxi Zhuang Birth Cohort (GZBC) between June 2015 and May 2017. Details on the ongoing cohort have been described in our previous study (Liang et al. 2020). In the present study, the mothers met the following enrollment criteria: (1) singleton pregnancy and gestational age < 13 weeks at enrollment; (2) intend to deliver at the study hospitals; (3) agreed to provide blood samples and birth tissues at delivery. A total of 801 mothers with singleton pregnancy who met the enrollment criteria were included in the final analysis. Ethics approval was approved by Guangxi Medical University (No.20140305-001). All participants were required to sign informed consent.
The demographic characteristics and lifestyles of pregnant women (e.g., age, pre-pregnancy weight, drinking, smoking, etc.) were collected via a face-to-face questionnaire. The maternal information (e.g., parity, pre-existing conditions, and pregnancy complications) and birth characteristics (e.g., gender, gestational age, and birth outcomes) were extracted from their medical records. Maternal pre-pregnancy body mass index (BMI, kg/m2) was calculated by their self-reported pre-pregnancy weight and height determined at their first examination. The last menstrual cycle of the mother was used to calculate gestational age (weeks). Ultrasound was used to estimate their gestational age if the menstrual date was uncertain.
Maternal blood samples were collected in a non-anticoagulated vacutainer during the first trimester. Cord blood was taken immediately from the umbilical vein into one ethylenediaminetetraacetic acid (EDTA) vacutainer after deliver. The blood samples were centrifuged at 3500 rpm for 10 min, and then separated into the plasma, serum, and hemocytes. After which, the samples were frozen at -80°C until assayed.
Maternal serum bisphenols measurement
Bisphenol levels in maternal serum were measured via ultra-high performance liquid chromatography-tandem mass spectrometry (UPLC-MS; Waters Xevo TQD) as previously described (Liang et al. 2020). Sample extraction and purification was used by liquid-liquid partition method. Briefly, 500 μl of serum and 2 μl of 1 μg/ml isotope internal standard containing BPA-d16 (98.0%, Sigma, St. Louis), BPS-12C13 (98.0%, Cambridge Osotope Lab, Inc.), and TBBPA-12C13 (99.0%, Cambridge Osotope Lab, Inc.) were buffered with 500 μl of 1.0 M sodium dihydrogen phosphate dihydrate buffer (pH=5.4) and hydrolysed enzymatically with 50 μl of β-glucuronidase/sulfatase (2,00 units/ml) at 37°C overnight. After incubation, 2 ml of n-hexane/acetone (7:3, v/v) was added, and then the mixtures were vortexed for 2 min, sonicated for 15 min, and centrifuged at 4,000 rpm for 10 min. After centrifugation, the upper clear liquid was collected into a new and clean glass tube. Subsequently, 2 ml of methyl-tert-butyl ether was added to the underlying mixtures instead of n-hexane/acetone (7:3, v/v). The above extraction steps (n-hexane/acetone) were repeated, and the supernatants were collected. The extracts were evaporated to dryness at 40°C. The residue was dissolved with 100 μl of methanol: 0.1% ammonia solution (50:50, v/v), filtered using a microporous membrane filter, and 10 μl of these resultants were injected for instrumental analysis.
The mass spectrometer was achieved with the negative-ion electrospray ionization mass spectrometer and multiple reaction monitoring (MRM) mode. The capillary voltage and desolvation temperature were 2.9 kV and 550°C, respectively. The desolvation gas and cone gas were set at flow rates of 50 L/h and 1000 L/h, respectively. The chromatographic separation was operated with an Acquity UPLC BEH C18 (1.7μm 2.1x100 mm, Waters, USA) analytical column. The column temperature was set at 45 °C. The mobile phase was 0.1% ammonia solution (A) and acetonitrile (B), delivered at a flow velocity of 0.3 ml/min; and linear gradient program, 30% B (0–1 min), 30%–100% B (1–4 min), 100% B (4–8 min), 100%–30% (8–8.5 min), and 100% B (8.5–10 min). BPA and its analogues levels were measured by comparing the chromatographic peak area of analytes with the area of the corresponding internal standards. BPA-d16 was used as internal standard for BPA, BPB and BPF, and BPS-12C13 and TBBPA-12C13 for BPS and TBBPA, respectively. The standards substance (BPA (≥ 98.0%); BPB (≥ 98.0%); BPF (≥ 98.0%); BPS (≥ 98.0%), and TBBPA (≥ 97%)) were obtained from Sigma, St. Louis. The calibration curves with ten points displayed satisfactory linearity (R2 > 0.990) over the range of bisphenol concentrations from 0.0 to 50.0 ng/ml. The limit of detection (LOD) was evaluated by the signal-to-noise ratio (S/N = 3). The LODs of BPA, BPB, BPF, BPS, and TBBPA were 0.193, 0.232, 0.507, 0.046, and 0.454 ng/ml, respectively. The average recoveries of the bisphenols ranged from 86.8% to 101.8%, and the relative standard deviations (RSD) were between 8.0% to 12.9%. Bisphenol was undetectable in the procedure blank and solvent blank. Glassware was used to avoid bisphenol contamination by plastic goods. The glassware was washed by methanol and ultra-pure water and then was baked.
Newborn relative leukocyte TL examination
The genomic DNA extraction of cord blood leukocytes was performed using a whole-genome DNA extraction kit and following the manufacturer's instructions (Aidlab Biotechnologies Co, Ltd, China). The concentration and purity of DNA were detected by the ultraviolet spectrophotometer (BioTek, USA). DNA purity was considered eligible if the OD260/OD280 rate was between 1.8 and 2.0.
Relative TL in cord blood was measured by StepOne Plus real-time polymerase chain reaction (RT-PCR; Applied Biosystems, USA). The relative TL was calculated as T/S ratio (telomere/single copy). The single-copy gene, 36B4, served as a reference gene in the telomere quantitative PCR signal. The primer sequence information of telomere and 36B4 gene is as follows: TeloF: 5′-CGGTTTGTTTGGGTTTGGGTTTGGGTTTGGGTTTGGGTT-3′; TeloR: 5′-GGCTTGCCTTACCCTTACCCTTACCCTTACCCTTACCCT-3′; 36B4F: 5′-CAGCAAGTGGGAAGGTGTAATCC-3′; 36B4R: 5′-CCCATTCTATCATCAACGGGTACAA-3′ (Cawthon 2002). Primers were purchased from Invitrogen, USA. The thermal cycling profile for the telomere PCR was: 50°C for 2 min, 95°C for 2 min, followed by 40 cycles of 95°C for 15 s and 62°C for 1 min. For the 36B4 PCR, the thermal cycling profile was: 50°C for 2 min, 95°C for 2 min, followed by 40 cycles of 95°C for 15 s, 62°C for 20 s, and 72°C for 1 min. The reference standard curve was created by using seven serial dilutions of human DNA (200, 50, 25, 125.5, 6.25, 3.125 ng/μl) and validated in each reaction batch. The 35 ng/ml standard concentration was selected as the positive control, and the DNA of the remaining samples were diluted with Tris-EDTA buffer to 25-50 ng/μl. All samples were run in triplicates. Each 96-well plate contained negative and positive controls to create a standard curve. Good linearity was observed across this range (R2 > 0.980). The Ct values generated were used to calculate the T/S ratio using the equation: Relative T/S = 2–(∆Ct Sample − ∆Ct control sample), where ∆Ct Sample=Sample Ct(TL) - Sample Ct(36B4) and ∆Ct control = Control Ct(TL) - Control Ct(36B4) (Cawthon 2002).
Statistical analysis
Bisphenol levels below the LOD were replaced by LOD/2. The Kolmogorov-Smirnov normality test was used to examine the distribution of data. We found that bisphenol concentrations and relative TL showed a skewed distribution. Thus, the relative TL and bisphenol concentrations were converted to the natural log (Ln) –transformation to enhance normality. Demographic characteristics were expressed as mean and standard deviation (Mean ± SD) values for normal distribution data, median (quartile range, IQR) values for skew distributional data, and number (percentages) for categorized data. The correlations of bisphenols were determined by using Spearman's correlation test.
Generalized linear models were used to estimate the regression coefficient (β) and standard error (SE) of the associations between bisphenol levels (Ln-transformed) and relative TL (Ln-transformed). The percentage change in relative TL for per doubling of bisphenol levels was calculated as (e(ln2 × β) − 1) × 100%, and the 95% confidence interval (95% CI) was calculated as (e[ ln 2 × (β ± 1.96 × SE] − 1) × 100% (Zhang et al. 2019). The covariates in the regression models included infant characteristics (gestational age, and infant gender) and maternal characteristics (maternal age, pre-pregnancy BMI, occupation, drinking before pregnancy, passive smoking during pregnancy, parity, hypertensive disorders in pregnancy (HDP), gestational diabetes mellitus (GDM)). We constructed different regression models, adjusted for infant characteristics only, maternal characteristics only and both of them. Considering the potential effect of maternal age and infant gender on the association between bisphenols and relative TL in offspring (Liang et al. 2020, Zhang et al. 2019), we conducted sub-group analysis stratified by maternal age (< 28 or ≥ 28 years, the median value) and infant gender. To ensure the robustness of the results, we also conducted sensitivity analyses by excluding women with GDM (n = 37) or HDP (n = 37).
To further investigate the potential non-linear relationship between bisphenol concentrations (Ln-transformed) and relative TL (Ln-transformed), restricted cubic spline (RCS) function was used with three knots located at the 10th, 50th, and 90th percentiles of the Ln-transformed bisphenol levels. The median Ln-transformed bisphenol levels were selected as reference. Wald chi-square test was used to test for the overall and non-linear associations (Desquilbet &Mariotti 2010). RCS analysis was also adjusted for the above potential confounders. Data analysis was performed using the R 4.0. All P values were two-tailed, and P < 0.05 was defined as statistical significance.