2.1 Study design
Participants of this study were women of the couples in the pregnancy cohort of the Shanghai Birth Cohort (SBC)[23]. In brief, from 2013 to 2016, pregnant women who visited for booking for prenatal care were approached for enrollment at six SBC participating hospitals in the four administrative districts of Shanghai city. Biological samples (maternal venous blood) were collected during the research process and blood Cd concentrations were measured using the whole blood. The collection and treatment of biological samples used unified standard operating procedures. The subjects’ basic information like social environment, educational levels, family income, smoking status and obstetric records of gestational diabetes mellitus (GDM) status during pregnancy were collected. The data related to newborns, regarding gender, gestational age, and anthropometric parameters, were recorded at birth by trained nurses. From 2018 to 2021, physical and cardiovascular ultrasonography tests were performed in the 4-year follow-up survey. We randomly selected 450 women to test for Cd concentration from this cohort. In our study, women who were miscarriage, stillbirth, multiple pregnancies, lost to follow-up, and without available medical records were excluded. For offspring, children with congenital heart disease, lost to follow-up, and without available records were excluded. Finally, 380 mother-offspring pairs were included in the analysis.
This study was approved by the Institutional Review Boards of Xinhua Hospital affiliated to Shanghai Jiao Tong University School of Medicine (XHEC-C-2013-001-2). Pregnant women were informed of the purposes and procedures of the study and signed an appropriate informed consent form at the beginning of the study.
2.2 Measurements of Cd in maternal blood in the first trimester
Maternal venous blood (gestational week ≤16 weeks) were collected using trace metal-free ethylenediaminetetraacetic acid (EDTA)-coated tubes by trained phlebotomists in the first trimester. Whole blood Cd levels were determined using inductively coupled plasma mass spectrometry (ICP-MS) (7500ce, Agilent, USA) and using validated, previously described laboratory methods[24]. The detection limits (LOD) for Cd using this method were 0.01 ug /L. Quality control for Cd was performed using a Levey-Jennings chart, with the mean±2SD values set as the allowable range. All samples were above these detection limits.
2.3 Measurement of child characteristics at the age of 4
Transthoracic echocardiography examinations were performed for the children according to the American and European guidelines[25] by trained operators using the Philip EPIQ7C (Philips Healthcare, Andover, USA) ultrasound that uses the X5-1 (1-5MHz) or S8-3 (8-3MHz) matrix-array transducers (Philips Healthcare, Andover, USA). Standard techniques were used to obtain M mode, 2D, and Doppler measurements in accordance with the American Society of Echocardiography[25, 26]. The following measurements were assessed: the subxiphoid long-axis view, the subxiphoid short-axis view, the apical 4-chamber view, the apical 3-chamber view, the apical 2-chamber view, the parasternal long-axis view, the parasternal short-axis view, the suprasternal long-axis and the suprasternal short-axis. And the following measurements were assessed: left ventricle diameter in systole (LVDs), left ventricular internal diameter in end-diastole (LVDd), left ventricular posterior wall thickness at systole (LVPWs), left ventricular posterior wall thickness at diastole (LVPWd), interventricular septal thickness at systole (IVSs), interventricular septal thickness at diastole (IVSd), left ventricular mass index (LVMI), end diastolic volume (EDV), end systolic volume (ESV), cardiac output (CO), stroke volume (SV), early diastolic mitral flow/late diastolic mitral flow (E/A ratio), Tei index (measured using the isovolumetric contraction time, isovolumetric relaxation times and ejection time), heart rate (HR), left ventricular ejection fractions (LVEF), and left ventricular fractional shortening (LVFS).
2.3 Statistical Analysis
Continuous variables are expressed as mean (SD) and categorical variables were expressed as numbers and percentages. Generalized Linear Model (GLM) was used to assess the linear association between maternal Cd exposure during early pregnancy and child cardiovascular measurements at the age of 4.
In the association analyses, a set of potential confounders were adjusted. They included maternal age (continuous), BMI before pregnancy (continuous), education status ("below college " and " college and above "), family annual income ("<100,000 CNY" and">100,000 CNY") and child gender. Except for age and BMI before pregnancy, which were continuously modeled in a linear form, the rest of selected covariates were modeled as categorical variables. Race/ ethnicity was not adjusted in all analyses since 99.2% of participants were of Chinese Han.
Causal mediation analysis[27] was conducted estimate the mediation effect of maternal Cd exposure on child's LV structure, which used generalized additive models (GAMs). Maternal Cd concentration was considered as exposure, child's cardiovascular structure as outcome, child's characteristics (gestational age and child BMI at age of 4 years) as potential mediator, maternal perinatal confounders (maternal age, prepregnancy BMI, education status, and annual income) as covariates.
To assess the robustness of our findings, we performed several sensitivity analyses. First, as Cd exposure was reported associated with GDM[28, 29], we excluded GDM women and then again conducted the association of maternal Cd exposure during pregnancy and child cardiovascular structure. Second, we used GAM to test nonlinear relationships between child cardiovascular structure and maternal Cd exposure, which presented by smooth spline curve. Finally, according to the results of GAMs, we further did piecewise linear regression models for LVDd and LVPWs stratified by Cd quartile.
The results were expressed by β (95% CI) and the P-values. The data were analyzed using Empower Stats (version 4, X&Y solutions, MA, USA) and R software (version 4.0.5). P <0.05 (two tail) was considered statistically significant.