Introduction
In recent years, over 350,000 chemicals and chemical mixtures have been registered for production, resulting in the continuous exposure of the general population to multiple contaminants, including airborne particles, industrial chemicals, heavy metals, and persistent organic pollutants (POPs) 1. These toxicants can enter the body via inhalation, ingestion, or dermal absorption and, once present in the circulatory system, they may interact with cells and tissues, interfering with essential biological processes. Such effects can be of particular concern during pregnancy, when maternal exposures can impact foetal development.
In this context, the Developmental Origins of Health and Disease (DOHaD) framework, also known as the Barker hypothesis, postulates that insults that occur during sensitive developmental windows can alter disease susceptibility across the lifespan 2. According to the “two-hit hypothesis”, early-life exposures may induce latent biological changes that manifest later in life when a second environmental or physiological stressor occurs. Indeed, there are epidemiological studies which support this model, showing associations between prenatal exposure to environmental toxicants and adverse birth and childhood outcomes, including impaired growth, neurodevelopmental deficits, and immune dysfunction 3, 4. These findings highlight that foetal development is shaped not only by genetic determinants but also by modifiable environmental factors 5.
A growing body of research indicates that epigenetic mechanisms mediate the long-term health effects of prenatal exposures. DNA methylation, histone modification, and gene silencing by microRNAs (miRNAs) are key processes linking environmental insults to altered gene regulation 6, 7, 8, 9. miRNAs, small non-coding RNAs that regulate gene expression and cellular homeostasis, are highly responsive to environmental stimuli 10. Their stability in biofluids and detectability in blood make them promising, non-invasive early biomarkers of exposure and biological impact 11. Several studies have shown that maternal exposure to chemical mixtures is associated with changes in circulating miRNA expression, with potential consequences for offspring development and later disease risk 12, 13, 14. Since many of these molecular alterations likely originate at the maternal–foetal interface, the placenta emerges as a key organ of interest for the way it mediates this interface by ensuring nutrient and gas exchange and acting as a barrier to xenobiotics. However, evidence indicates that this barrier is not impermeable: POPs, polycyclic aromatic hydrocarbons (PAHs), and even particulate matter can cross into the foetal compartment 15, 16. Such transplacental transfer raises concerns about direct foetal exposure to complex contaminant mixtures during critical windows of intrauterine development.
Despite increasing evidence linking pollutant exposure, miRNA dysregulation, and health outcomes, very few studies have simultaneously investigated both maternal and cord blood to assess maternal-foetal transfer and intergenerational effects. To address this knowledge gap, birth cohorts enrolled in highly contaminated areas are a powerful tool for investigation. The Neonatal Environment and Health Outcomes (NEHO) cohort, established in southern Italy, provides such an opportunity. Conducted in the industrialized Augusta-Priolo area, recognized as a National Priority Contaminated Site (NPCS), NEHO was specifically designed to evaluate how environmental pollutants can affect the health of mothers and their children 17, 18, 19.
A
For this study, we investigated the relationships between maternal exposure to environmental pollutants and circulating miRNA expression in maternal serum during the third trimester of pregnancy and in cord serum at birth. We i) assessed the concentration of selected inorganic and organic contaminants and their placental transfer, ii) examined their associations with miRNA expression in both maternal and foetal sides, and iii) explored the biological pathways potentially linking these molecular alterations to health outcomes. By integrating biomonitoring data, the evaluation of biomarkers, and questionnaire data, this work aims to advance our understanding of maternal-foetal health dynamics in polluted environments, as well as to provide insights which are relevant to public health interventions.
Evaluation of metal and organic compound levels in maternal and cord serum
NEHO participants exhibited widespread exposure to a variety of environmental pollutants. Seven POPs — including HCB, p,p′-DDE, PCB-118, PCB-138, PCB-153, and PCB-180 — and five trace elements, comprising three EEs (Cu, Zn, Se) and two non-essential elements known for their toxicity (Hg, As), were quantified in maternal and cord serum samples as described above. The concentrations of analytes in both maternal and cord serum are presented in Table S1, along with the percentage of values below the LOQ in the 95 samples analysed. Compounds with LOQ frequencies exceeding 15% in either matrix were excluded from further analysis.
Table 2 lists the analytes included in the subsequent analyses. All selected chemicals showed statistically significant differences between the two matrix concentrations (p < 0.001). The placental transfer ratio provides insight into placental permeability and potential foetal exposure to environmental contaminants. A higher ratio indicates greater transfer efficiency from mother to foetus. Among EEs, the placental transfer ranged from 0.13 for Cu to 0.76 for Se, with Zn being the only element exhibiting a median ratio above 1 (1.41 [1.17;1.59]). In contrast, OC compounds displayed consistently low placental transfer ratios, with values ranging from 0.21 to 0.52, indicating limited permeability. Similarly, PCBs showed low transplacental transfer, with values ranging from 0.12 to 0.24 across individual congeners, reflecting their restricted transplacental movement.
Table S1
Maternal vs Cord Blood pollutant concentrations. Mean (standard deviation) and median [interquartile range] were reported for both maternal and cord serum chemicals. Chemicals with less than 15% of subjects having values above the limit of quantification (LOQ) are highlighted in bold. The excluded pollutants for each side are reported in red.
| | | MATERNAL N = 95 | CORD BLOOD N = 95 | p-value |
|---|
CHEMICALS | UNITS | MEAN (SD) | MEDIAN [IQR] | LOQ (%) | MEAN (SD) | MEDIAN [IQR] | LOQ (%) | |
Se | ug/L | 77.8 (12.2) | 78.0 [68.0;86.5] | 0.00 | 52.1 (9.50) | 51.0 [46.5;56.0] | 0.00 | < 0.001 |
Hg | 0.84 (0.80) | 0.57 [0.20;1.09] | 31.25 | 0.64 (0.60) | 0.45 [0.20;0.82] | 45.83 | |
As | 1.28 (2.25) | 0.50 [0.50;1.00] | 71.00 | 1.21 (1.84) | 0.50 [0.50;0.50] | 77.10 | |
Zn | 656 (133) | 639 [572;735] | 0.00 | 889 (128) | 899 [804;970] | 0.00 | < 0.001 |
Cu | 2171 (385) | 2120 [1942;2302] | 0.00 | 417 (242) | 371 [295;476] | 0.00 | < 0.001 |
HCB | ng/g of total lipids | 56.2 (19.6) | 53.0 [42.6;66.3] | 0.00 | 25.5 (21.4) | 18.9 [14.0;27.1] | 0.00 | < 0.001 |
TNC | 7.13 (6.65) | 6.19 [2.50;8.95] | 34.37 | 2.56 (0.58) | 2.50 [2.50;2.50] | 97.92 | |
DDE | 422 (445) | 319 [231;486] | 0.00 | 107 (117) | 86.8 [52.3;123] | 0.00 | < 0.001 |
PCB74 | 10.8 (5.93) | 8.96 [7.23;13.5] | 8.33 | 2.81 (1.40) | 2.50 [2.50;2.50] | 92.71 | < 0.001 |
PCB118 | 22.4 (12.5) | 18.9 [14.6;28.5] | 1.04 | 4.42 (3.17) | 2.50 [2.50;5.87] | 64.58 | < 0.001 |
PCB138 | 64.9 (39.8) | 54.1 [41.7;85.9] | 0.00 | 12.5 (8.46) | 10.4 [7.59;14.8] | 7.29 | < 0.001 |
PCB153 | 115 (75.5) | 94.3 [73.1;150] | 0.00 | 21.3 (14.2) | 17.3 [12.8;26.9] | 1.04 | < 0.001 |
PCB156 | 10.5 (7.60) | 8.55 [6.32;13.1] | 17.71 | 2.64 (0.71) | 2.50 [2.50;2.50] | 94.79 | |
PCB180 | 86.4 (66.2) | 67.7 [50.1;108] | 0.00 | 13.8 (10.3) | 10.3 [7.17;18.3] | 10.41 | < 0.001 |
PCB183 | 8.35 (7.00) | 7.01 [2.50;11.0] | 27.08 | 2.58 (0.56) | 2.50 [2.50;2.50] | 96.87 | |
Table 2
Maternal vs cord pollutant concentrations. Mean (standard deviation) and median [interquartile range] were reported for both maternal and cord serum chemicals. Chemicals with more than 15% of subjects having values above the limit of quantification (LOQ) are highlighted in red. Statistically significant differences between cord and maternal chemical concentrations were computed using the paired Wilcoxon test. The placental transfer ratio of chemicals measured in matched samples is reported either as mean (standard deviation) or median [interquartile range].
CHEMICAL | UNITS | MATERNAL N = 95 | CORD N = 95 | p-value PAIRED (Wilcoxon test) | Placental transfer (CORD/MATERNAL) |
|---|
| | | MEAN (SD) | MEDIAN [IQR] | LOQ (%) | MEAN (SD) | MEDIAN [IQR] | LOQ (%) | | MEAN (SD) | MEDIAN [IQR] |
Se | µg/L | 77.8 (12.2) | 78.0 [68.0;86.5] | 0.00 | 52.1 (9.50) | 51.0 [46.5;56.0] | 0.00 | < 0.001 | 0.68 (0.15) | 0.64 [0.58;0.76] |
Zn | 656 (133) | 639 [572;735] | 0.00 | 889 (128) | 899 [804;970] | 0.00 | < 0.001 | 1.40 (0.32) | 1.41 [1.17;1.59] |
Cu | 2171 (385) | 2120 [1942;2302] | 0.00 | 417 (242) | 371 [295;476] | 0.00 | < 0.001 | 0.20 (0.13) | 0.17 [0.13;0.23] |
HCB | ng/g of total lipids | 56.2 (19.6) | 53.0 [42.6;66.3] | 0.00 | 25.5 (21.4) | 18.9 [14.0;27.1] | 0.00 | < 0.001 | 0.52 (0.78) | 0.35 [0.27;0.52] |
DDE | 422 (445) | 319 [231;486] | 0.00 | 107 (117) | 86.8 [52.3;123] | 0.00 | < 0.001 | 0.26 (0.10) | 0.26 [0.21;0.30] |
PCB-74 | 10.8 (5.93) | 8.96 [7.23;13.5] | 8.33 | 2.81 (1.40) | 2.50 [2.50;2.50] | 92.71 | | | |
PCB-118 | 22.4 (12.5) | 18.9 [14.6;28.5] | 1.04 | 4.42 (3.17) | 2.50 [2.50;5.87] | 64.58 | | | |
PCB-138 | 64.9 (39.8) | 54.1 [41.7;85.9] | 0.00 | 12.5 (8.46) | 10.4 [7.59;14.8] | 7.29 | < 0.001 | 0.20 (0.10) | 0.19 [0.15;0.24] |
PCB-153 | 115 (75.5) | 94.3 [73.1;150] | 0.00 | 21.3 (14.2) | 17.3 [12.8;26.9] | 1.04 | < 0.001 | 0.20 (0.10) | 0.19 [0.15;0.23] |
PCB-180 | 86.4 (66.2) | 67.7 [50.1;108] | 0.00 | 13.8 (10.3) | 10.3 [7.17;18.3] | 10.42 | < 0.001 | 0.17 (0.10) | 0.16 [0.12;0.20] |
Table S2
Maternal vs Cord Blood miRNA deltaCt. Mean (standard deviation) and median [interquartile range] were reported for both maternal and cord serum. miRNAs with less than 30% of subjects having Ct values equal to 40 are highlighted in bold. The excluded miRNAs for each side are reported in red.
| | MATERNAL N = 95 | CORD BLOOD N = 95 | p-value |
|---|
MEAN (SD) | MEDIAN [Q1-Q3] | SUB. VALUE (%) | MEAN (SD) | MEDIAN [Q1-Q3] | SUB. VALUE (%) | |
|---|
hsa-let-7a-5p | 2.38 (0.51) | 2.34 [2.13;2.69] | 0.00 | 2.53 (0.83) | 2.71 [2.34;2.94] | 0.00 | < 0.001 |
hsa-miR-1-3p | 10.0 (1.92) | 9.84 [8.94;11.1] | 4.17 | 12.4 (2.17) | 12.3 [11.3;13.6] | 11.46 | < 0.001 |
hsa-miR-100-5p | 6.51 (1.25) | 6.76 [5.64;7.35] | 0.00 | 7.55 (1.28) | 7.62 [6.90;8.18] | 4.17 | < 0.001 |
hsa-miR-106b-5p | 1.91 (0.54) | 1.91 [1.60;2.21] | 0.00 | 1.25 (0.75) | 1.23 [0.81;1.62] | 0.00 | < 0.001 |
hsa-miR-10b-5p | 6.49 (1.33) | 6.37 [5.44;7.45] | 0.00 | 6.65 (1.43) | 6.53 [5.86;7.42] | 3.13 | 0.435 |
hsa-miR-122-5p | 5.21 (1.85) | 5.29 [3.91;6.40] | 0.00 | 4.25 (2.08) | 4.32 [2.87;5.39] | 4.17 | < 0.001 |
hsa-miR-124-3p | 13.9 (3.72) | 13.8 [11.3;16.9] | 32.29 | 12.3 (2.36) | 12.3 [11.3;13.7] | 13.54 | 0.001 |
hsa-miR-125b-5p | 5.93 (1.81) | 5.96 [5.04;6.54] | 1.04 | 5.59 (1.10) | 5.81 [5.03;6.30] | 2.08 | 0.244 |
hsa-miR-126-3p | 0.23 (0.61) | 0.16 [-0.18;0.58] | 0.00 | 0.90 (1.11) | 1.05 [0.62;1.46] | 0.00 | < 0.001 |
hsa-miR-133a-3p | 8.36 (2.38) | 8.01 [6.96;8.94] | 7.29 | 10.4 (2.15) | 10.1 [9.23;11.2] | 13.54 | < 0.001 |
hsa-miR-133b | 8.69 (1.57) | 8.63 [7.79;9.53] | 1.04 | 10.9 (2.08) | 10.8 [9.81;11.8] | 10.42 | < 0.001 |
hsa-miR-134-5p | 7.98 (2.02) | 7.55 [6.69;8.84] | 3.13 | 6.36 (1.03) | 6.21 [5.67;6.85] | 4.17 | < 0.001 |
hsa-miR-141-3p | 6.73 (1.22) | 6.87 [6.02;7.57] | 0.00 | 9.21 (1.76) | 9.34 [8.85;9.80] | 8.33 | < 0.001 |
hsa-miR-143-3p | 4.94 (0.86) | 4.99 [4.50;5.46] | 0.00 | 4.36 (1.31) | 4.32 [3.64;5.14] | 3.13 | < 0.001 |
hsa-miR-146a-5p | 2.23 (0.88) | 2.24 [1.84;2.73] | 0.00 | 2.75 (1.00) | 2.65 [2.28;3.07] | 0.00 | < 0.001 |
hsa-miR-150-5p | 2.44 (1.08) | 2.37 [1.56;3.15] | 0.00 | 2.39 (1.22) | 2.56 [1.82;3.13] | 0.00 | 0.656 |
hsa-miR-155-5p | 8.76 (1.30) | 8.57 [7.96;9.13] | 2.08 | 8.52 (1.17) | 8.67 [8.27;9.08] | 3.13 | 0.626 |
hsa-miR-17-5p | 0.91 (0.51) | 0.88 [0.63;1.11] | 0.00 | 0.89 (1.05) | 0.73 [0.36;1.23] | 0.00 | 0.052 |
hsa-miR-17-3p | 9.04 (1.22) | 8.87 [8.26;9.49] | 1.04 | 7.50 (1.47) | 7.15 [6.71;8.13] | 4.17 | < 0.001 |
hsa-miR-18a-5p | 4.08 (0.55) | 4.06 [3.76;4.49] | 0.00 | 3.73 (0.70) | 3.74 [3.37;3.99] | 2.08 | < 0.001 |
hsa-miR-192-5p | 4.83 (0.95) | 4.89 [4.18;5.51] | 0.00 | 4.32 (0.76) | 4.41 [3.96;4.85] | 2.08 | < 0.001 |
hsa-miR-195-5p | 8.49 (1.00) | 8.34 [7.92;8.83] | 0.00 | 9.29 (1.63) | 9.40 [8.42;10.2] | 6.25 | < 0.001 |
hsa-miR-196a-5p | 15.0 (2.02) | 15.2 [13.6;16.2] | 33.33 | 12.5 (2.09) | 12.7 [11.5;13.6] | 10.42 | < 0.001 |
hsa-miR-19a-3p | -1.17 (0.60) | -1.14 [-1.49;-0.90] | 0.00 | -1.45 (0.77) | -1.56 [-1.87;-1.20] | 0.00 | < 0.001 |
hsa-miR-19b-3p | -0.99 (0.59) | -0.96 [-1.34;-0.69] | 0.00 | -1.06 (0.64) | -1.08 [-1.41;-0.74] | 0.00 | 0.299 |
hsa-miR-200a-3p | 12.7 (2.49) | 12.1 [11.0;14.3] | 25.00 | 11.5 (1.83) | 11.7 [10.7;12.4] | 11.46 | 0.004 |
hsa-miR-200b-3p | 11.6 (1.74) | 11.4 [10.6;12.2] | 8.33 | 11.9 (1.80) | 12.2 [11.3;12.8] | 7.29 | 0.001 |
hsa-miR-200c-3p | 7.97 (1.50) | 7.88 [7.43;8.47] | 2.08 | 8.98 (1.39) | 9.30 [8.70;9.66] | 5.21 | < 0.001 |
hsa-miR-203a-3p | 12.8 (2.85) | 12.3 [10.9;14.6] | 25.00 | 12.0 (2.65) | 12.2 [10.7;13.2] | 18.75 | 0.210 |
hsa-miR-205-5p | 7.67 (1.61) | 7.65 [6.46;8.71] | 0.00 | 8.76 (2.45) | 8.45 [7.37;9.82] | 9.38 | 0.001 |
hsa-miR-208a-3p | 17.6 (1.98) | 17.9 [16.5;18.8] | 87.50 | 16.7 (3.07) | 17.2 [15.7;18.3] | 58.33 | 0.038 |
hsa-miR-20a-5p | 0.04 (0.43) | 0.04 [-0.20;0.31] | 0.00 | -0.13 (0.83) | -0.19 [-0.59;0.10] | 0.00 | 0.002 |
hsa-miR-21-5p | -1.10 (0.42) | -1.11 [-1.28;-0.83] | 0.00 | -0.48 (0.75) | -0.44 [-0.73;-0.24] | 1.04 | < 0.001 |
hsa-miR-210-3p | 6.19 (0.81) | 6.19 [5.75;6.55] | 0.00 | 4.09 (1.27) | 3.78 [3.48;4.13] | 2.08 | < 0.001 |
hsa-miR-214-3p | 9.32 (1.47) | 9.12 [8.38;10.2] | 0.00 | 7.85 (1.79) | 8.08 [7.14;8.94] | 1.04 | < 0.001 |
hsa-miR-215-5p | 5.46 (0.95) | 5.49 [4.83;6.07] | 0.00 | 4.86 (0.81) | 4.93 [4.47;5.43] | 2.08 | < 0.001 |
hsa-miR-221-3p | 2.10 (0.97) | 1.97 [1.40;2.63] | 0.00 | 3.15 (1.26) | 2.76 [2.37;3.66] | 2.08 | < 0.001 |
hsa-miR-222-3p | 3.41 (0.47) | 3.39 [3.13;3.65] | 0.00 | 4.63 (1.12) | 4.45 [4.19;4.74] | 3.13 | < 0.001 |
hsa-miR-223-3p | -2.57 (0.70) | -2.61 [-3.02;-2.17] | 0.00 | -1.84 (1.15) | -1.69 [-2.36;-1.21] | 0.00 | < 0.001 |
hsa-miR-224-5p | 7.31 (1.48) | 7.14 [6.33;7.97] | 1.04 | 6.22 (1.45) | 6.07 [5.31;6.86] | 4.17 | < 0.001 |
hsa-miR-23a-3p | -0.66 (0.60) | -0.70 [-1.08;-0.38] | 0.00 | 0.73 (0.95) | 0.69 [0.39;1.05] | 0.00 | < 0.001 |
hsa-miR-25-3p | 0.61 (0.73) | 0.71 [0.31;1.14] | 0.00 | 0.13 (0.78) | 0.23 [-0.22;0.63] | 0.00 | < 0.001 |
hsa-miR-27a-3p | 0.22 (0.73) | 0.15 [-0.29;0.59] | 0.00 | 1.29 (0.78) | 1.22 [0.81;1.63] | 1.04 | < 0.001 |
hsa-miR-296-5p | 9.66 (1.45) | 9.41 [8.84;10.3] | 2.08 | 8.64 (2.03) | 8.33 [7.69;9.18] | 6.25 | < 0.001 |
hsa-miR-29a-3p | 3.85 (0.77) | 3.79 [3.31;4.26] | 0.00 | 4.60 (0.90) | 4.56 [3.98;5.07] | 2.08 | < 0.001 |
hsa-miR-30d-5p | 1.71 (0.56) | 1.75 [1.33;2.00] | 0.00 | 2.18 (0.49) | 2.25 [2.01;2.45] | 0.00 | < 0.001 |
hsa-miR-34a-5p | 11.2 (1.93) | 10.8 [9.91;11.9] | 7.29 | 10.3 (2.13) | 10.2 [9.05;11.0] | 12.50 | < 0.001 |
hsa-miR-375-3p | 8.97 (1.97) | 8.94 [7.62;9.94] | 2.08 | 10.6 (2.10) | 10.2 [9.39;11.8] | 9.38 | < 0.001 |
hsa-miR-423-5p | 3.86 (0.69) | 3.77 [3.42;4.30] | 0.00 | 3.48 (0.79) | 3.43 [3.00;4.10] | 2.08 | < 0.001 |
hsa-miR-499a-5p | 14.6 (2.92) | 14.3 [12.4;17.1] | 39.58 | 12.5 (2.21) | 12.2 [11.5;13.6] | 19.79 | < 0.001 |
hsa-miR-574-3p | 7.53 (0.85) | 7.46 [7.01;8.04] | 0.00 | 8.17 (1.52) | 8.40 [7.73;8.93] | 1.04 | < 0.001 |
hsa-miR-885-5p | 8.79 (2.11) | 8.59 [7.36;9.85] | 2.08 | 9.16 (2.03) | 8.82 [7.70;10.4] | 5.21 | 0.145 |
hsa-miR-9-5p | 12.7 (2.46) | 11.6 [11.0;15.2] | 22.92 | 11.3 (1.83) | 11.1 [10.4;12.3] | 12.50 | 0.001 |
hsa-miR-92a-3p | -1.18 (0.64) | -1.07 [-1.70;-0.74] | 0.00 | -1.10 (0.59) | -1.09 [-1.38;-0.83] | 0.00 | 0.559 |
hsa-miR-93-5p | 1.36 (0.41) | 1.40 [1.16;1.64] | 0.00 | 0.96 (1.26) | 0.78 [0.33;1.20] | 3.13 | < 0.001 |
hsa-let-7c-5p | 3.79 (0.58) | 3.77 [3.40;4.08] | 0.00 | 3.54 (0.71) | 3.62 [3.30;3.89] | 0.00 | 0.028 |
hsa-miR-107 | 2.13 (0.63) | 2.08 [1.81;2.32] | 0.00 | 1.60 (0.98) | 1.41 [1.24;1.67] | 2.08 | < 0.001 |
hsa-miR-10a-5p | 9.71 (1.00) | 9.59 [9.01;10.4] | 0.00 | 10.1 (1.40) | 10.2 [9.59;10.8] | 5.21 | 0.001 |
hsa-miR-128-3p | 7.21 (0.71) | 7.14 [6.68;7.57] | 0.00 | 7.05 (0.71) | 7.03 [6.73;7.28] | 3.13 | 0.390 |
hsa-miR-130b-3p | 5.39 (0.55) | 5.43 [5.03;5.79] | 0.00 | 4.91 (0.72) | 4.91 [4.54;5.23] | 2.08 | < 0.001 |
hsa-miR-145-5p | 3.87 (0.84) | 3.88 [3.49;4.38] | 0.00 | 4.62 (1.07) | 4.54 [3.96;5.03] | 4.17 | < 0.001 |
hsa-miR-148a-3p | 3.27 (0.61) | 3.23 [2.93;3.49] | 0.00 | 2.48 (0.86) | 2.40 [2.09;2.71] | 2.08 | < 0.001 |
hsa-miR-15a-5p | 0.93 (0.78) | 0.91 [0.48;1.47] | 0.00 | 0.04 (0.82) | -0.09 [-0.39;0.41] | 1.04 | < 0.001 |
hsa-miR-184 | 16.6 (2.79) | 17.2 [15.3;18.5] | 75.00 | 16.2 (3.03) | 16.9 [14.8;18.3] | 61.46 | 0.361 |
hsa-miR-193a-5p | 7.44 (1.08) | 7.45 [6.60;8.11] | 0.00 | 7.35 (1.16) | 7.31 [6.77;7.92] | 4.17 | 0.539 |
hsa-miR-204-5p | 10.1 (1.48) | 9.97 [9.35;10.8] | 2.08 | 11.6 (2.25) | 11.6 [10.4;12.9] | 11.46 | < 0.001 |
hsa-miR-206 | 14.4 (3.69) | 14.3 [11.2;17.7] | 36.46 | 12.8 (2.72) | 12.9 [11.2;14.3] | 18.75 | 0.005 |
hsa-miR-211-5p | 15.4 (2.19) | 15.2 [13.9;17.1] | 51.04 | 16.0 (2.92) | 16.1 [14.9;18.0] | 59.38 | 0.009 |
hsa-miR-26b-5p | 2.63 (0.45) | 2.54 [2.34;3.00] | 0.00 | 2.02 (0.67) | 2.12 [1.84;2.44] | 0.00 | < 0.001 |
hsa-miR-30e-5p | 3.52 (0.89) | 3.26 [2.90;3.96] | 0.00 | 3.91 (1.53) | 3.73 [3.39;4.05] | 3.13 | < 0.001 |
hsa-miR-372-3p | 15.6 (3.18) | 15.9 [13.1;18.2] | 62.50 | 15.1 (3.50) | 15.4 [13.1;17.6] | 55.21 | 0.480 |
hsa-miR-373-3p | 15.2 (2.94) | 15.2 [12.8;17.7] | 54.17 | 13.6 (3.00) | 13.4 [12.0;15.4] | 31.25 | 0.001 |
hsa-miR-374a-5p | 3.27 (0.48) | 3.33 [3.00;3.59] | 0.00 | 3.04 (0.88) | 3.20 [2.76;3.51] | 2.08 | 0.074 |
hsa-miR-376c-3p | 4.14 (1.41) | 3.98 [3.06;4.88] | 0.00 | 2.95 (0.92) | 2.69 [2.38;3.43] | 2.08 | < 0.001 |
hsa-miR-7-5p | 8.46 (1.55) | 8.47 [7.70;9.15] | 2.08 | 6.72 (0.94) | 6.82 [6.31;7.22] | 3.13 | < 0.001 |
hsa-miR-96-5p | 9.95 (1.59) | 9.84 [8.97;10.7] | 1.04 | 9.89 (1.75) | 9.85 [8.90;10.6] | 8.33 | 0.872 |
hsa-miR-103a-3p | 1.95 (0.67) | 1.93 [1.58;2.22] | 0.00 | 1.22 (0.70) | 1.08 [0.94;1.39] | 0.00 | < 0.001 |
hsa-miR-15b-5p | 0.84 (0.36) | 0.88 [0.57;1.10] | 0.00 | 0.77 (0.52) | 0.88 [0.55;1.08] | 0.00 | 0.765 |
hsa-miR-16-5p | -2.59 (0.66) | -2.52 [-2.96;-2.20] | 0.00 | -2.65 (0.82) | -2.58 [-2.92;-2.14] | 0.00 | 0.745 |
hsa-miR-191-5p | 1.51 (0.57) | 1.52 [1.20;1.86] | 0.00 | 2.05 (0.69) | 2.10 [1.78;2.38] | 1.04 | < 0.001 |
hsa-miR-22-3p | 1.57 (0.66) | 1.48 [1.15;1.93] | 0.00 | 1.21 (0.92) | 1.05 [0.69;1.54] | 1.04 | < 0.001 |
hsa-miR-24-3p | 0.29 (0.64) | 0.25 [-0.10;0.65] | 0.00 | 1.80 (0.81) | 1.64 [1.38;2.04] | 2.08 | < 0.001 |
hsa-miR-26a-5p | 1.89 (0.67) | 1.82 [1.51;2.26] | 0.00 | 1.84 (0.82) | 1.82 [1.40;2.30] | 1.04 | 0.805 |
hsa-miR-31-5p | 10.3 (1.74) | 10.3 [9.14;11.4] | 1.04 | 10.9 (2.30) | 11.1 [10.1;12.3] | 1.04 | 0.002 |
Table S3
Regression analysis results for each analyte vs miRNAs in maternal serum.
Analyte | miRNA | Beta (95% CI) | P-value |
|---|
Cu | hsa-miR-7-5p | 0.29 (0.09; 0.50) | 0.043 |
Se | hsa-miR-214-3p | 0.30 (0.10; 0.51) | 0.037 |
Zn | hsa-miR-7-5p | 0.25 (0.08; 0.44) | 0.043 |
hsa-miR-214-3p | 0.33 (0.14; 0.53) | 0.02 |
hsa-miR-10b-5p | 0.38 (0.19; 0.57) | 0.003 |
hsa-miR-143-3p | 0.33 (0.14; 0.53) | 0.02 |
hsa-miR-148a-3p | 0.29 (0.10; 0.49) | 0.037 |
hsa-miR-205-5p | 0.33 (0.14; 0.53) | 0.02 |
hsa-miR-206 | 0.36 (0.17; 0.56) | 0.006 |
hsa-miR-30d-5p | -0.33 (-0.53; -0.13) | 0.02 |
hsa-miR-374a-5p | -0.32 (-0.52; -0.13) | 0.03 |
p,p′-DDE | hsa-miR-1-3p | 0.28 (0.07; 0.48) | 0.04 |
hsa-miR-133a-3p | 0.27 (0.07; 0.47) | 0.03 |
HCB | hsa-miR-1-3p | 0.40 (0.20; 0.60) | 0.002 |
hsa-miR-133a-3p | 0.36 (0.16; 0.55) | 0.007 |
PCB74 | hsa-miR-1-3p | 0.28 (0.08; 0.49) | 0.04 |
hsa-miR-133p | 0.28 (0.07; 0.48) | 0.034 |
PCB118 | hsa-miR-1-3p | 0.36 (0.16; 0.57) | 0.005 |
hsa-miR-133b | 0.35 (0.15; 0.55) | 0.007 |
hsa-miR-195-5p | 0.28 (0.07; 0.49) | 0.02 |
PCB138 | hsa-miR-1-3p | 0.28 (0.06; 0.50) | 0.047 |
hsa-miR-195-5p | 0.36 (0.14; 0.58) | 0.003 |
hsa-miR-30e-5p | 0.31 (0.09; 0.53) | 0.02 |
PCB153 | hsa-miR-195-5p | 0.42 (0.20; 0.63) | 0.0009 |
hsa-miR-30e-5p | 0.34 (0.12; 0.56) | 0.02 |
hsa-miR-30e-5p | 0.34 (0.11; 0.57) | 0.02 |
PCB180 | hsa-miR-195-5p | 0.43 (0.21; 0.66) | 0.0009 |
hsa-miR-30e-5p | 0.34 (0.11; 0.57) | 0.02 |
hsa-miR-30e-5p | 0.32 (0.10; 0.53) | 0.02 |
ΣPCBs | hsa-miR-195-5p | 0.42 (0.20; 0.64) | 0.0009 |
hsa-miR-30e-5p | 0.34 (0.11; 0.57) | 0.01 |
Table S4
Regression analysis results for each analyte vs miRNAs in cord serum.
Analyte | miRNA | Beta (95% CI) | P-value |
|---|
Cu | hsa-miR-30d-5p | 0.35 (0.15;0.54) | 0.002 |
hsa-miR-7-5p | 0.29 (0.08;0.49) | 0.03 |
hsa-miR-195-5p | 0.32 (0.13;0.52) | 0.013 |
hsa-miR-122-5p | 0.49 (0.31;0.68) | < 0.001 |
hsa-miR-134-5p | 0.35 (0.15;0.54) | 0.006 |
hsa-miR-200c-3p | 0.28 (0.08;0.48) | 0.032 |
hsa-miR-574-3p | 0.29 (0.09;0.49) | 0.025 |
hsa-miR-10a-5p | 0.28 (0.08;0.48) | 0.028 |
hsa-miR-376-3p | 0.38 (0.19;0.57) | 0.002 |
Se | hsa-miR-30d-5p | 0.37 (0.17;0.57) | 0.002 |
hsa-miR-7-5p | 0.27 (0.06;0.48) | 0.031 |
hsa-miR-29a-3p | -0.41 (-0.61;-0.21) | 0.001 |
hsa-miR-22-3p | -0.31 (-0.52;-0.10) | 0.04) |
Zn | hsa-miR-25-3p | -0.24 (-0.44;-0.04) | 0.047 |
hsa-miR-16-5p | -0.32 (-0.52;-0.12) | 0.02 |
HCB | hsa-let-7c-5p | 0.41 (0.20;0.61) | 0.001 |
hsa-miR-374a-5p | 0.32 (0.11;0.53) | 0.04 |
hsa-miR-191a-5p | 0.33 (0.12;0.54) | 0.04 |
PCB138 | hsa-miR-25-3p | 0.26 (0.05;0.48) | 0.047 |
hsa-miR-29a-3p | -0.36 (-0.57;-0.15) | 0.004 |
hsa-miR-30d-5p | 0.28 (0.07;0.49) | 0.02 |
hsa-miR-7-5p | 0.29 (0.07;0.51) | 0.03 |
hsa-miR-133a-3p | -0.37 (-0.58;-0.16) | 0.007 |
hsa-miR-133b-3p | -0.34 (-0.55;-0.13) | 0.02 |
PCB153 | hsa-miR-25-3p | 0.27 (0.05;0.48) | 0.047 |
hsa-miR-29a-3p | -0.34 (-0.56;-0.12) | 0.006 |
hsa-miR-30d-5p | 0.26 (0.04;0.47) | 0.03 |
hsa-miR-7-5p | 0.26 (0.04;0.48) | 0.04 |
hsa-miR-133a-3p | -0.35 (-0.56;-0.13) | 0.009 |
hsa-miR-133b-3p | -0.30 (-0.52;-0.10) | 0.03 |
PCB180 | hsa-miR-25-3p | 0.27 (0.04;0.50) | 0.047 |
hsa-miR-29a-3p | -0.31 (-0.55;-0.10) | 0.02 |
∑PCBs | hsa-miR-25-3p | 0.27 (0.01;0.53) | 0.047 |
hsa-miR-29a-3p | -0.35 (-0.57;-0.12) | 0.006 |
hsa-miR-30d-5p | 0.27 (0.04;0.49) | 0.03 |
hsa-miR-7-5p | 0.27 (0.04;0.50) | 0.04 |
hsa-miR-133a-3p | -0.34 (-0.56;-0.12) | 0.009 |
hsa-miR-133b | -0.30 (-0.53;-0.10) | 0.03 |
Table S5
bWQS results from maternal serum.
miRNA | EE mixture β(95%CI) | POP mixture β(95%CI) | Total mixture β(95%CI) | |
|---|
hsa-miR-1-3p | | 0.41 (0.18; 0.64)$ | 0.60 (0.28; 0.87) | |
hsa-miR-7-5p | 0.39 (0.17; 0.60)$ | | 0.44 (0.05; 0.80) | |
hsa-miR-10b-5p | 0.47 (0.21; 0.68)$ | | | |
hsa-miR-15a-5p | | -0.23 (–0.47; − 0.02) | -0.25 (–0.54; − 0.02) | |
hsa-miR-17-5p | | 0.29 (0.07; 0.47) | | |
hsa-miR-18a-3p | -0.35 (–0.64; − 0.11) | | | |
hsa-miR-20a-5p | | 0.27 (0.0006; 0.46) | | |
hsa-miR-26a-5p | | 0.23 (0.01; 0.45) | | |
hsa-miR-27a-3p | -0.31 (–0.57; − 0.07) | | | |
hsa-miR-30d-5p | -0.34 (–0.64; − 0.12)$ | | | |
hsa-miR-30e-5p | | 0.30 (0.07; 0.52)$ | | |
hsa-miR-31-5p | | -0.25 (–0.48; − 0.02) | | |
hsa-miR-122-5p | 0.34 (0.09; 0.56) | | | |
hsa-miR-126-3p | -0.30 (–0.52; − 0.05) | | | |
hsa-miR-128-3p | | -0.25 (–0.48; − 0.06) | -0.30 (–0.57; − 0.02) | |
hsa-miR-130b-3p | -0.35 (–0.62; − 0.10) | | | |
hsa-miR-133a-3p | | 0.47 (0.21; 0.72)$ | 0.64 (0.32; 0.90) | |
hsa-miR-133b | -0.39 (0.09; 0.65) | | 0.54 (0.10; 0.88) | |
hsa-miR-134-5p | -0.34 (–0.63; − 0.12) | 0.28 (0.07; 0.50) | | |
hsa-miR-143-3p | 0.40 (0.20; 0.64)$ | | | |
hsa-miR-145-5p | 0.46 (0.15; 0.71) | | 0.52 (0.12; 0.87) | |
hsa-miR-195-5p | | 0.39 (0.15; 0.58)$ | 0.52 (0.17; 0.81) | |
hsa-miR-193a-3p | 0.37 (0.07; 0.58) | | | |
hsa-miR-200a-3p | 0.33 (0.03; 0.53) | | 0.41 (0.01; 0.74) | |
hsa-miR-200b-3p | | | | |
hsa-miR-203a-3p | | 0.27 (0.04; 0.47) | | |
hsa-miR-205-5p | 0.39 (0.08; 0.63)$ | 0.27 (0.06; 0.46) | | |
hsa-miR-214-3p | 0.65 (0.25; 1.00)$ | | 0.65 (0.25; 1.00) | |
hsa-miR-296-5p | | | -0.28 (–0.67; − 0.10) | |
hsa-miR-374a-5p | -0.40 (–0.67; − 0.21)$ | | | |
hsa-miR-376c-3p | -0.33 (–0.60; − 0.11) | | | |
$ indicates if the miRNA was previously identified in single-chemical regression of analytes included in the same family. | |
Table S6
bWQS results from cord serum.
miRNA | EE mixture β(95%CI) | POP mixture β(95%CI) | Total mixture β(95%CI) |
|---|
hsa-let-7c-5p | | 0.42 (0.10;0.64)$ | |
hsa-miR-7-5p | 0.32 (0.08;0.53)$ | | 0.46 (0.15;0.71) |
Hsa-miR-17-3p | | -0.30 (-0.61;-0.04) | -0.35 (-0.82;-0.002) |
hsa-miR-19b-3p | | | -0.32 (0.75;-0.01) |
hsa-miR-22-3p | -0.39 (-0.64;-0.17)$ | | -0.54 (-0.91;-0.23) |
hsa-miR-27a-5p | | -0.28 (-0.58;-0.06) | |
hsa-miR-29a-3p | -0.42 (-0.66;-0.22)$ | -0.39 (-0.65;-0.15)$ | -0.68 (-1.00;-0.41) |
hsa-miR-30d-5p | 0.41 (0.18;0.60)$ | | |
hsa-miR-34-5p | | -0.28 (-0.58;-0.05) | |
hsa-miR-122-5p | 0.62 (0.38;0.86)$ | | 0.73 (0.41;1.03) |
hsa-miR-124-3p | | -0.28 (-0.57;-0.01) | |
hsa-miR-125b-5p | 0.35 (0.02;0.60) | | 0.35 (0.01; 0.67) |
hsa-miR-128-3p | | -0.23 (-0.51;-0.01) | |
hsa-miR-133a-3p | | -0.43 (-0.69;-0.20)$ | -0.48 (0.88;-0.21) |
hsa-miR-133b | | -0.39 (-0.67;-0.19)$ | -0.42 (-0.72;-0.17) |
hsa-miR-134-5p | 0.37 (0.12;0.59)$ | | |
hsa-miR-141-3p | -0.27 (-0.54;-0.03) | | -0.37 (-0.74;-0.10) |
hsa-miR-150-5p | -0.23 (-0.48;-0.01) | | |
hsa-miR-195-5p | 0.39 (0.18;0.59)$ | | 0.40 (0.07;0.70) |
hsa-miR-200a-3p | | -0.23 (-0.51;-0.01) | |
hsa-miR-200c-3p | | | 0.36 (0.001;0.64) |
hsa-miR-205-5p | | -0.28 (-0.53;-0.01) | |
hsa-miR-206 | 0.30 (0.02;0.57) | | |
hsa-miR-214-3p | -0.29 (-0.61;-0.02) | | |
hsa-miR-221-3p | 0.27 (0.01;0.55) | | |
hsa-miR-296-5p | | -0.22 (-0.51;-0.02) | |
hsa-miR-375-3p | 0.30 (0.07;0.55) | | |
hsa-miR-376c-3p | 0.43 (0.23;0.62)$ | | 0.43 (0.02; 0.75) |
hsa-miR-499a-5p | | -0.32 (-0.59;-0.10) | -0.34 (-0.69;-0.02) |
$ indicates if the miRNA was previously identified in single-chemical regression of analytes included in the same family. |
Table S7
bWQS chemical weight threshold.
| | EEs | POPs | TOTAL |
|---|
tEEm | tPOPm | tTOTm |
|---|
MATERNAL | (1/3) = 0.33 | (1/7) = 0.14 | (1/10) = 0.10 |
| | tEEc | tPOPc | tTOTc |
CORD | 1/3 = 0.33 | (1/5) = 0.20 | (1/8) = 0.125 |
Table S8
Comparison of maternal and cord serum bWQS results.
| | MATERNAL | CORD |
|---|
miRNA | MIXTURE | DIRECTION | MIXTURE | DIRECTION |
|---|
hsa-miR-7-5p | EEs | POS | EEs | POS |
hsa-miR-27a-5p | EEs | NEG | POPs | NEG |
hsa-miR-30d-5p | EEs | NEG | EEs | POS |
| | | TOT | POS |
hsa-miR-122-5p | EEs | POS | EEs | POS |
hsa-miR-128-3p | POPs | NEG | POPs | NEG |
Hsa-miR-133a-3p | POPs | POS | POPs | NEG |
hsa-miR-133b | EEs | POS | POPs | NEG |
hsa-miR-134-5p | EEs | NEG | EEs | POS |
POPs | POS | | |
hsa-miR-195-5p | POPs | POS | EEs | POS |
hsa-miR-200a-3p | EEs | POS | POPs | NEG |
hsa-miR-205-5p | EEs | POS | POPs | NEG |
hsa-miR-214-3p | EEs | POS | EEs | NEG |
hsa-miR-296-5p | TOT | NEG | POPs | NEG |
hsa-miR-376c-3p | EEs | NEG | EES | POS |
To further investigate the distribution of analytes at the maternal–foetal interface, we examined the statistically significant correlations among analytes in maternal serum (Fig. 1A), in cord serum (Fig. 1B), and between the two matrices for each analyte (Fig. 1C). Strong intra-group correlations were observed among substances with similar chemical characteristics in maternal serum (Fig. 1A). EEs were weakly correlated with each other. Zn and Se did not show any significant inter-element correlations. Among OCs, HCB and p,p’-DDE were highly correlated (r = 0.79). HCB also showed significant correlations with all PCB congeners; however, the strength of association showed an evident reduction in the Pearson correlation index with increasing chlorine substitution. Indeed, the correlation was strongest with PCB-118 (r = 0.68, five chlorine atoms) and lowest with PCB-180 (r = 0.33, seven chlorine atoms). A similar trend was observed for p,p’-DDE, but only for the less chlorinated PCBs (PCB-74, PCB-118, PCB-138, and PCB-153). Finally, high inter-correlations were also found among PCB congeners, especially those with similar chlorination patterns.
In cord serum (Fig. 1B), a strong correlation between Se and Cu (r = 0.66) was observed. Both elements also exhibited weaker correlations with compounds having different chemical characteristics. Specifically, Se showed a slight positive correlation with PCB-138 (r = 0.22) and a negative correlation with HCB (r = -0.20), whereas Cu was positively correlated with PCB-138 (r = 0.36) and PCB-153 (r = 0.27). OCs were not correlated with each other in cord serum but showed weak correlations with PCB congeners. As shown for the maternal side, the PCB congeners in cord serum remained strongly inter-correlated. Finally, several analytes showed moderate correlations between both the maternal and cord compartments (Fig. 1C), reinforcing the idea that, although placental transfer is generally limited, foetal exposure is more likely to be influenced by maternal levels. Although the placental transfer values of OCs and PCBs were very low, p,p’-DDE exhibited the highest correlated index between maternal and cord serum (r = 0.96), followed by various PCBs. Among the EEs, although Zn was the only analyte showing placental transfer above 1, Se was the only one to show a slight maternal-cord correlation (r = 0.22).
Mixture component contributions to miRNA expression
To further explore the drivers of the observed miRNA alterations, we examined the contribution of the individual components within pollutant mixtures to overall miRNA expression patterns on both the maternal and foetal side. The bWQS chemical weight thresholds are listed in Table S7.
As shown in Fig. 6 (left side, EE section), Zn exceeded the contribution threshold (tEEm=0.33) in 94% of the significant associations, followed by Cu (22%), and Se (17%) in maternal serum.
In the POP section, PCB-118 emerged as the most frequent contributor (61% of the significant associations), followed by HCB and PCB-74 (both 54%). PCB-138 and PCB-180 exceeded the weight threshold in 31% of the associations, while p,p’-DDE and PCB-153 exceeded the 0.14 threshold (tPOPm=0.14) in 23% of the associations.
Overall, in maternal serum, as highlighted in Fig. 6 (left side, TOTAL section), Se and Zn emerged as the predominant elements across the total mixture, exceeding the weight threshold in 81% and 72% of the significant associations, respectively. In contrast, Cu overcame the threshold in 36% of associations. Among the POPs, PCB-118 and HCB were consistently above the threshold (tTOTm=0.10) in 54% and 45% of the associations, respectively. PCB-74 and PCB-138 showed weights above the threshold in 27% of the cases, while p,p’-DDE and PCB-153 exceeded the 0.10 threshold in 18% of the associations. PCB-180 was the least influential as 10% of associations reported PCB-180 as having weight above the threshold.
Regarding the cord serum, among the elements included in the EE mixture, Cu most frequently emerged as a key driver of the observed associations (62%), while Se and Zn showed a lower contribution frequency (both 37%) (Fig. 6 right side, EE section).
For the POP section, as illustrated in Fig. 6 (right side), PCB-138 emerged as the most recurrent contributor across the significant associations, exceeding the relevant threshold (tPOPc=0.20) in 85% of the cases. PCB-153 and HCB also demonstrated substantial contributions, with weights over the threshold in 61% and 54% of the significant associations, respectively. In contrast, PCB-180 exceeded the threshold in only 23% of the associations, while p,p’-DDE was the least frequent (15% of the associations).
Overall, in cord serum, Se emerged as the most frequent component above the threshold, in 53% of significant associations (Fig. 6 right side, TOTAL section), followed by Cu and PCB-153, each contributing in 47% of the cases. HCB was involved in 40% of associations, while Zn and p,p’-DDE exceeded the threshold (tTOTc=0.125) in 33% of associations. Finally, PCB-138 and PCB-180 were the least represented, contributing above the threshold in 27% and only 0.7% of associations, respectively.
Comparative miRNA-chemical associations in maternal and cord compartments
A
A heatmap was generated to summarize the statistically significant associations between individual chemicals, chemical mixtures, and circulating miRNAs (Fig. 7). The results obtained from maternal serum are represented in pink, whereas those from cord serum are shown in light blue. In the left section of the heatmap, the associations identified through linear regression models (one miRNA vs one chemical) are reported. P-values were corrected for multiple testing using the FDR, and only associations with adjusted p < 0.05 were included. The second section of the heatmap illustrates the results of the bWQS analyses, aimed at evaluating mixture effects on miRNA regulation. Three mixture models were tested, as previously described: EEs, including only micronutrients; POPs, and the TOTAL mixture, including all exposures. In this analysis, the colour scale of the cells reflects the direction of the estimates, with positive associations shown in red and negative ones in blue.
These data summarize the differences highlighted by the two approaches (single versus mixture analysis).
Figure 7: Heatmap showing statistically significant associations between chemicals and miRNAs. Pink results (left up) were obtained from the maternal side, whereas blue results (left bottom) were from the cord serum. The first part of the heatmap reports the results of linear regression models (unique miRNA vs unique chemical), correcting for maternal age, maternal BMI, gestational age (days), gender, smoking habits. The p-value correction for multiple comparison was performed using the FDR correction. Associations having adjusted p-values < 0.05 were included in the heatmap. In the second part, bWQS results were included in order to evaluate the mixture effect in miRNA regulation. Three different analyses were conducted: EEs, in which the mixture included only essential elements, POPs, including CB, DDE, and PCBs, and TOTAL, using the total mixture. The cell colours indicate the estimates (positive in red and negative in blue).
In contrast, a few miRNAs displayed divergent patterns: miR-30d-5p and miR-376-3p were negatively associated with EEs in maternal serum but positively associated in cord serum. In addition, hsa-miR-30d-5p showed a positive association with the total mixture in cord serum. MiR-133a-3p displayed an opposite trend for POPs of the two sides: positive direction in the maternal serum, and negative in the cord serum. A similar trend was observed for miR-134-5p which exhibited an additional positive association with the POP mixture in maternal serum. Conversely, miR-214-3p was positively associated with EEs in maternal serum and negatively in cord serum.
Some miRNAs were linked to different mixtures across matrices. For example, miR-133b, miR-200a-3p, and miR-205-5p were positively associated with EEs in maternal serum but negatively associated with POPs in cord serum. MiR-27a-5p showed negative associations in both matrices but versus different mixtures: EEs in maternal serum and POPs in cord serum. miR-195-5p maintained a consistent positive direction of association in both matrices, although in relation to different mixtures—POPs in maternal serum and EEs in cord serum. Finally, miR-296-5p displayed negative associations with both the total mixture in maternal serum and the POP mixture in cord serum. An additional schematic version of Fig. 8 is reported in Table S8.
Pathway enrichment analysis
Pathway enrichment analysis of miRNA target genes revealed distinct functional signatures associated with the EE and POP mixtures. miRNAs were grouped according to exposure type (EEs, POPs, or TOT mixtures), direction of the association (positive or negative), and biological origin (maternal or cord side) to identify clusters of enriched pathways. These clusters were organized into communities based on functional similarity, with each community centred around a hub pathway representing the most interconnected term. The heatmap in Fig. 9 displays the hub pathways with at least 4 associated miRNAs. Red indicates the presence of one or more miRNAs positively associated in the bWQS analysis, while blue is used for negative associations. Results related to the EE mixture are indicated by the purple bar, followed by those for the POP and TOT mixtures with the orange and green bars, respectively. Maternal-side associations are displayed on the left, and cord-side associations are on the right.
miRNAs detected in maternal serum and positively associated with the EE mixture were predominantly enriched in three key hub pathways with the lowest adjusted p-value: “Signaling by Receptor Tyrosine kinases”, “Signaling by ERBB2”, and “RAF/MAP kinase cascade”. The latter also includes miRNAs that were negatively associated with the EE mixture, along with pathways such as “Disease of signal transduction by growth factor receptors and second messengers” and “Signal Transduction”.
For the POP mixture, the top enriched pathways by maternal serum miRNAs were “Cytokine Signaling in immune system” and “Cell Cycle, Mitotic”. In addition, several hub pathways, “Apoptosis”, “Metabolism”, and “Disease”, were enriched exclusively by miRNAs positively associated with the POP mixture. On the other hand, pathways such as “Developmental Biology”, “Signaling by NOTCH”, “Generic Transcription Pathway”, and “Cellular response to stress” included both positively and negatively associated miRNAs.
Regarding the total mixture, miRNAs with statistically significant associations were primarily included in “Signaling by Nuclear Receptors”, which shows a consistent trend among the three mixture analyses. Another recurrent hub pathway is “Cytokine Signaling in Immune system”. The “Hemostasis”, “Generic Transcription Pathway”, and “Cellular response to Stress” were enriched exclusively by miRNAs negatively associated with the total mixture. The “Metabolism” hub pathway reflected trends observed with the POP mixture, along with “Developmental Biology” and “Signaling by NOTCH”.
A
In cord serum, most of the pathways were enriched by miRNAs positively associated with the EE mixture, with the exception of “Hemostasis” and “Signaling by RHO GTPases” which were enriched by negatively associated miRNAs. The most significantly enriched pathways in this context included “Diseases of signal transduction by growth factor receptors and second messengers” and “Signaling by NTRK1”. Unlike the maternal profile, “Developmental Biology”, “Metabolism”, and “Signaling by NOTCH” were not enriched in cord serum, while the “Toll like receptor 4 cascade” emerged as a pathway specifically enriched on the foetal side. The remaining pathways show a similar trend to those observed in maternal serum. For the POP mixture, we observed that only negatively associated miRNAs contributed to enriching the pathways. The top enriched pathways were “Disease”, “Signaling by Receptor Tyrosine kinase”, and “Signaling by Rho GTPase”. Several hub pathways appeared to be exclusively enriched by cord miRNAs associated with the POP mixture, including “Signaling by Rho GTPases”, “Innate Immune system”, “Signaling by NOTCH1 in Cancer”, “Diseases of glycosylation”, “Adherens junction interactions”, and “Fatty acid metabolism”. Finally, miRNAs associated with the total mixture in cord serum primarily enriched the “Disease” and “Signaling by NTRK1” pathways.
The detailed contribution of each miRNA to the statistically significant enriched pathways is reported
A
in Supplementary Figure S1.
Discussion
This study provides a comprehensive understanding of the epigenetic consequences of prenatal exposure to essential elements and toxic environmental contaminants in a birth cohort enrolled in a highly industrialized area, integrating pollutant and miRNA expression profiling. Our work aimed to investigate i) placental function in pollutant transfer, ii) differential circulating miRNA expression between maternal and foetal compartments, and iii) the influence of chemical mixtures on miRNA expression on both sides.
By simultaneously analysing EEs and POPs in matched maternal and cord serum samples of the NEHO birth cohort, we provide evidence that the placental maternal–foetal interface is not only a site of selective chemical permeability but also a determinant of compartment-specific epigenetic regulation. Placental transfer behaviour varied considerably across compounds. Among the EEs analysed, Zn showed a median transfer ratio greater than 1, consistent with active and regulated transport mechanisms supporting foetal development 20, 21. In contrast, Se and Cu displayed moderate to low transfer efficiencies, while lipophilic POPs such as OCs and PCBs exhibited limited placental passage, with transfer ratios well below 0.5.
These findings align with prior reports which observed that the physicochemical properties of POPs, particularly their hydrophobicity and protein-binding affinity, limit their mobility across the placenta 22. Strong intra-group correlations were found among chemicals with similar structures, while cross-group correlations were generally weak; notable exceptions included Se with Cu and PCB-156 in maternal serum, and Se with Cu in cord serum. Indeed, HCB and p,p’-DDE were highly correlated, as were PCB congeners, with correlation strength decreasing as chlorine substitution increased. These correlations highlight the placenta’s role as a selective barrier that modulates, rather than simplistically prevents, pollutant transfer to the foetus.
Furthermore, the literature data indicates that stressors such as environmental pollutants can induce compartment-specific molecular alterations, including differential miRNA expression, reflecting maternal adaptation and foetal vulnerability 23.
These data support our initial approach in which we analysed the chemical-specific transfer at the placenta maternal–foetal interface shaping differences in miRNA regulation, with profiles serving as a readout of compartment-specific biology during critical developmental windows. For these reasons, we evaluated the association between environmental pollutants and circulating miRNAs in the maternal and cord blood sides of the placenta in the NEHO birth cohort samples. In maternal serum, we observed that higher concentrations of Cu, Se, and Zn were significantly associated with changes in several miRNAs 24. Specifically, Cu and Se concentrations were positively associated with miR-7-5p and miR-214-3p, respectively. In this context, miR-214-3p is known to be as a putative regulator of GPx4, a Se-dependent selenoprotein, whereas miR-7 is recognized as a key player in cancer biology and in the modulation of drug resistance. By contrast, Zn exhibited both positive and negative associations across different miRNAs. Notably, miR-30d-5p and miR-374a-5p were inversely associated with maternal Zn levels, suggesting a possible specific regulatory mechanism of Zn on miRNA expression participating in maternal–foetal communication and pregnancy complications 13.
In cord serum, Cu levels were positively associated with a broader set of miRNAs, including miR-7-5p, miR-10a-5p, miR-30d-5p, miR-134-5p, miR-122-5p, miR-195-5p, miR-200c-3p, miR-376c-3p, and miR-574-3p. Conversely, Se showed opposite trends depending on the selected miRNAs, with inverse associations observed for miR-22-3p and miR-29a-3p showing that many cord plasma miRNAs respond to maternal and neonatal metabolic traits, reinforcing the idea that cord miRNA profiles are broadly exposure-sensitive in pregnancy 25.
Exposure to OCs, including p,p’-DDE, HCB, and various PCB congeners, was also associated with distinct miRNA expression patterns. In maternal serum, p,p’-DDE, HCB, and lower-chlorinated PCBs (e.g., PCB-74, PCB-118) were positively associated with miR-1-3p and miR-133a-3p, which are known to be involved in muscle development and cardiac function (myomiRs) 26. Interestingly, PCB-118 and PCB-138 showed broader effects, being associated with multiple miRNAs. The most chlorinated PCBs (PCB-153 and PCB-180) were primarily linked to miR-195-5p and miR-30e-5p, associations that were also confirmed using the sum of PCB concentrations. Increasing data demonstrates that both miR-195-5p and miR-30e-5p are key regulatory miRNAs in pregnancy. MiR-195-5p supports trophoblast proliferation, migration, and invasion and interfaces with angiogenic and activin/Nodal pathways. Its circulating levels also track the anti-angiogenic imbalance characteristic of preeclampsia 27. MiR-30e-5p has been linked to vascular and immune/inflammatory adaptations since it associates with preeclampsia trajectories and microvascular dysfunction and varies with late-pregnancy blood-pressure phenotypes 28.
In cord serum, HCB modulated a specific subset of miRNAs with increased ΔCt values, let-7c-5p, miR-191a-5, and miR-374a-5p. Importantly, miR-25-3p and miR-29a-3p 29 showed consistent associations with all measured PCB congeners, though miR-25-3p was positively associated whereas miR-29a-3p was inversely associated, suggesting divergent regulatory responses to these chemicals. Data from the literature shows that miR-25-3p regulates myometrial contractility by being exported in trophoblast extracellular vesicles, and its downregulation enhances Ca²⁺ signaling, suggesting a role in both uterine activity and foetal development 30. Furthermore miR-29a-3p, enriched in villus-derived exosomes, promotes immune tolerance by suppressing decidual NK-cell IFN-γ production; its reduction has been reported in recurrent pregnancy loss 31. Finally, miR-7-5p 32and miR-30d-5p 33 were positively associated with PCB-138 and PCB-153, whereas miR-133a-3p and miR-133b-3p were negatively associated with the same congeners. These trends suggest a potential disruption of miRNAs involved in muscle and cardiovascular development and reinforce the sensitivity of the foetal miRNome to environmental exposures.
As a first approach, these single-pollutant models were very informative for identifying specific associations between exposures (chemicals) and biological responses (miRNAs), also showing qualitative and quantitative responses in both the maternal and foetal compartments. However, their use would oversimplify the real-world scenario in which individuals are simultaneously exposed to a complex mixture of contaminants.
Therefore, moving toward multipollutant approaches provides a more realistic framework for analysing the effects of mixtures, highlighting potential key drivers, and can better capture the overall impact of environmental exposures on maternal and foetal health, particularly in highly industrialized areas characterized by complex pollutant profiles. By extending the analysis to include mixture models through bWQS regression, we identified several additional miRNAs associated with combined exposures to EEs and POPs, many of which were not detectable in single-chemical models.
Our multipollutant analysis revealed both consistent and divergent associations between environmental mixtures and circulating miRNAs in maternal and cord serum, suggesting shared as well as compartment-specific regulatory responses. For instance, miR-7-5p and miR-122-5p were positively associated with essential element mixtures in both maternal and cord samples, indicating stable modulation by micronutrient-related exposures and supporting their role as common biomarkers of exposure across compartments. In contrast, several miRNAs demonstrated opposite directions of association between mother and child, including miR-133a-3p, miR-133b, miR-200a-3p, miR-205-5p, and miR-214-3p, which may reflect differential maternal–foetal adaptive responses or compartment-specific regulation of gene expression in response to environmental stressors (consistent with the known roles of these miRNAs in trophoblast invasion, Epithelial–mesenchymal transition, endothelial signaling, and placental dysfunction) 34.
Such divergent associations suggest that maternal and foetal systems develop distinct molecular programs to cope with environmental stress, reflecting compartment-specific priorities in maintaining pregnancy homeostasis. Notably, miR-195-5p showed consistent positive associations with both POP and EE mixtures in maternal and cord serum, reinforcing its central role in placental and vascular biology 35. By contrast, miR-30d-5p, miR-134-5p, and miR-376c-3p shifted direction depending on the exposure mixture considered, underlining the complexity of multipollutant effects and the necessity of integrated analytical approaches to capture them. Together, these findings suggest that, while some miRNAs act as robust biomarkers of shared maternal–foetal exposure, others capture regulatory programmes, underscoring the value of mixture-aware models in elucidating the molecular impact of environmental exposures during pregnancy.
Considering the placental transfer dynamics, our results illustrate how complex chemical exposures shape compartment-specific miRNA responses. Overall, the bWQS regression demonstrated that, in the maternal serum, Zn emerged as the most frequent contributor to mixture-associated miRNA expression, whereas Cu dominated on the foetal side, emphasizing not only the chemical-specific but also the compartment-specific drivers of epigenetic regulation. Among POPs, PCB-118 and HCB were the most influential components in maternal serum, whereas PCB-138 and HCB were predominant in the cord compartment. These patterns suggest that the most biologically active or bioavailable pollutant components may differ across maternal and foetal systems, probably due to differences in toxicokinetics, metabolism, or placental transfer efficiency. Indeed, few miRNA–pollutant associations overlapped across compartments, despite strong correlations in pollutant levels. This divergence points to independent molecular responses shaped by tissue-specific regulation, chemical concentration, developmental context, and possibly differences in cell-type composition of circulating shuttles such as extracellular vesicles.
To further derive the biological significance of these epigenetic signals, we performed pathway enrichment analysis of the predicted miRNA target genes. In order to analyse the structures from these multi-dimensional signals in detail, miRNAs were grouped by exposure class [EEs, POPs, or total mixture (TOT)], direction of association in bWQS (positive vs negative), and biological origin (maternal vs cord). The general structure of the enriched pathways was organized into general hubs. Only pathways supported by at least four miRNAs are displayed, increasing confidence that the observed communities are not driven by single markers.
The analysis indicates that miRNAs associated with pollutant mixtures converge on biological processes of central relevance to pregnancy. Among the most prominent are apoptosis, cell cycle regulation, and cellular stress responses, all of which are essential for placental development and trophoblast turnover 36. Enrichment of immune-related pathways, including cytokine signaling, innate immune responses, and the TLR4 cascade, is particularly noteworthy as these processes underpin the delicate balance between maternal immune tolerance and controlled inflammation at the maternal–foetal interface 37. In addition, pathways such as NOTCH signaling, Rho GTPase signaling, and adherens junction interactions suggest potential effects on trophoblast migration, vascular remodelling, and intercellular communication—mechanisms frequently disrupted in pregnancy complications 13. The enrichment of metabolic pathways (e.g., fatty acid metabolism, glycosylation, nuclear receptor signaling) further points to a possible impact of pollutant exposures on maternal–foetal metabolic homeostasis 38. Although maternal and cord compartments showed overlapping pathway signatures, notable differences were observed: maternal profiles were more strongly enriched in stress- and apoptosis-related pathways, whereas cord profiles were characterized by immune and developmental signaling, suggesting compartment-specific adaptive responses. Collectively, these findings support the view that environmental mixtures influence pregnancy through interconnected networks regulating placental growth, immune balance, and metabolic function, with potential implications for adverse outcomes. These pathway structures demonstrate that mixture composition and biological compartment jointly determine miRNA-mediated regulation.
The maternal programme emphasizes growth factor signaling, nuclear receptor integration, and immuno-metabolic coupling, whereas the foetal programme highlights neurotrophic, innate immune, and cytoskeletal/adhesion axes.
Because miRNAs typically repress their targets, apparent increases or decreases in pathway enrichment must be interpreted in light of target directionality and network feedback, although validation at the mRNA/protein level is warranted.
The strength of this study lies in its integrative design, combining high-resolution chemical analysis with miRNA profiling in matched maternal and foetal serum, and applying both single-pollutant and advanced mixture modelling approaches. To the best of our knowledge, this combined approach represents a significant and original improvement in this field. However, some limitations must be acknowledged. The sample size, while sufficient for detecting robust associations, may have limited our ability to explore more specific interactions, and the cross-sectional design does not permit inference of causality. In addition, while circulating miRNAs are promising biomarkers, they may not fully capture the tissue-specific regulatory changes occurring at the maternal–foetal interface or within target organs. Nonetheless, the consistency of key findings across models and the alignment with established biological pathways support their interpretive validity.
In summary, our data demonstrate that environmental exposures during pregnancy induce distinct and compartment-specific miRNA responses in maternal and foetal circulation, shaped by individual chemical properties, complex mixture effects, and the placenta’s filtering role. The identification of exposure-sensitive miRNAs involved in developmental, oxidative stress, and immune pathways provides mechanistic insight into how early-life exposures may contribute to long-term health trajectories. These results emphasize the importance of considering both the maternal and foetal compartments, accounting for chemical mixtures and integrating molecular biomarkers to better understand the developmental origins of health and disease. Furthermore, our findings also highlight the pivotal role of the placenta as a selective barrier and an active regulator in mediating foetal exposure and shaping molecular responses to environmental stressors. However, whether the placenta actively sends distinct regulatory signals to the maternal and foetal sides or rather acts as a selective barrier that separates or integrates molecular cues, remains to be clarified and warrants further investigation.
Altogether, these results advance our understanding of how the interplay between environmental factors and placental function contributes to the regulation of gene expression in early development, setting the stage for a detailed discussion of the observed pollutant-specific transfer dynamics and miRNA expression profiles.
Methods
Study population
This study was conducted within the framework of the NEHO birth cohort, an ongoing longitudinal study designed to investigate the effects of environmental exposures on maternal and child health. A total of 95 mother-child pairs were selected from the NEHO cohort for the present analysis.
A
Participants were recruited at the General Hospital of Lentini (Syracuse, Italy) and the Umberto I Hospital in Syracuse (Italy). The study involved mothers residing across a broad geographic area that included both municipalities located within the NPCS of Augusta-Priolo and its surrounding areas
39.
The design and recruitment procedures of the NEHO cohort have been described in detail previously 19. In brief, between January 2018 and January 2020, 561 pregnant women were enrolled based on their continuous residence within the Augusta-Priolo area for at least 12 months prior to enrolment. Participants were recruited during prenatal visits at local healthcare facilities, and eligibility was confirmed based on standardized inclusion criteria, including residency, age ≥ 18 years, and the absence of pre-existing major chronic illnesses.
Data collection
Upon enrolment, detailed socio-demographic, lifestyle, medical, and environmental exposure information was obtained through structured web-based questionnaires. Follow-up assessments were conducted at 6, 12, and 24 months postpartum to gather longitudinal data on maternal health, infant development, and potential ongoing exposures. Clinical and delivery-related information—including gestational age at delivery, mode of delivery, and newborn anthropometric parameters (birthweight, length, and head circumference)—was extracted from medical records by trained healthcare personnel at the time of birth. For this analysis, randomly selected mother–child pairs were included to ensure a sample representative of the general NEHO cohort inclusion criteria 18.
Blood sampling and handling
A
A
The research project was approved by the relevant Ethics Committee (Comitato Etico Catania 2, 11 July 2017, n. 38/2017/CECT2) and strictly followed the principles outlined in the Declaration of Helsinki (64th WMA General Assembly, Fortaleza, Brazil, October 2013). All adopted procedures were compliant with the European Union’s General Data Protection Regulation (EU 2016/679) and Italian laws concerning data protection. To this end, all questionnaire data and biological samples were pseudonymized using ID tracking numbers. Blood samples were collected from all participants during the last trimester of pregnancy. Serum samples were prepared as previously described
18. After serum was separated by centrifugation, samples were temporarily stored at − 20°C in each maternity unit and then transported on dry ice to the NEHO biobank for long-term storage at − 80°C.
Analytical methods for trace metals and organic compounds
The analysis of POPs was performed at the Chemical Exposure Unit of the National Institute for Health and Welfare (Kuopio, Finland) using an Agilent 7000B gas chromatograph coupled to a triple quadrupole mass spectrometer (GC–MS/MS), as previously described 14. Trace element analyses were conducted at the micropollutant unit of LERES (Laboratoire d'Étude et de Recherche en Environnement et Santé), at the French School of Public Health (EHESP, Rennes, France), as described in Davies et al. (2021)40. A total of 95 maternal and cord serum samples were analysed by inductively coupled plasma tandem mass spectrometry (ICP-MS/MS) after acid mineralization.
Certified reference materials and internal quality control samples (spiked blood and serum) were included in each analytical batch to ensure accuracy and reproducibility.
Maternal and cord serum concentrations of 15 chemicals were evaluated in this study: three Essential Elements—EEs (Se, Zn, and Cu), arsenic (As), mercury (Hg), three different organochlorine pesticides—OCs (HCB, Trans-Nonachlor TNC, and p’-p’-DDE) and seven polychlorinated biphenyls—PCBs (PCB-74, PCB-118, PCB-138, PCB-153, PCB156, PCB-180, PCB-183). Concentrations below the limit of quantification (LOQ) were imputed using a value equal to LOQ/2, following the recommendations of the statistical analysis plan developed within the HBM4EU project (https://www.hbm4eu.eu/wp-content/uploads/2018/09/Deliverable-10.2-Statistical-Analysis-Plan.pdf). Raw concentrations of POPs were normalized to total lipid levels by dividing both maternal and cord serum chemical concentrations by total lipid levels 40. Both EEs and lipid-normalized concentrations of POPs were then log-transformed and further standardized using the standard normal variate (SNV) transformation in order to minimize data skewness and allow for meaningful comparison across analytes with differing chemical properties. Only chemicals for which at least 15% of samples had concentrations above the LOQ were included in the further analyses. Placental transfer was computed for each analyte by dividing cord serum concentrations by maternal serum concentration between matched samples. The ratios were calculated only when data for both matrices were available (LOQ < 15%).
Serum miRNA extraction and expression analysis
A
Total RNAs enriched with miRNAs were extracted from 200 µl of maternal or cord serum samples according to the miRNeasy Serum/Plasma Advanced Kit manufacturer’s protocol (Qiagen, Milan, Italy). Reverse transcription reactions were performed using the miRCURY LNA RT kit (Qiagen, Milan, Italy), according to the manufacturer’s protocol. Real-time analyses were performed using the Applied Biosystems StepOnePlus™ Real-Time PCR System for maternal samples, and the Applied Biosystems QuantStudio™ 3 Real-Time PCR System for cord samples; levels of 84 miRNAs (listed in Table S2) were detected using the Human Serum/Plasma Small Focus (miScript) array panel (Qiagen, Milan, Italy, cat. YAHS-206ZC, formerly MIHS-106) and the miRCURY LNA Sybr Green PCR kit (Qiagen, Milan, Italy).
Raw cycle threshold (Ct) values were calibrated on internal spike-in controls using the Gene Globe platform (
https://geneglobe.qiagen.com/it/analyze). Undetermined values of Ct were estimated as 40 Ct, which is the last cycle of the reactions. deltaCt (ΔCt) was computed, correcting for endogenous controls, using the following equation:
in which endogenous controls are indicated as end = hsa-miR-30c-5p, hsa-miR-103a-3p, hsa-miR-451a, hsa-miR-23a-3p.
MiRNAs with more than 30% undetermined values were excluded from the analysis.
Covariates
Covariates included self-reported information on maternal and child characteristics and behaviours; these were collected at enrolment using the “Baseline” questionnaire or at delivery through clinical birth records. Specifically, the following variables were collected at enrolment: maternal age (continuous, years), maternal pre-pregnancy BMI (continuous, kg/m2), smoking during pregnancy (binary, yes/no). Additional data were extracted from hospital delivery records, including child sex (binary, male/female) and gestational age (continuous, days). Alcohol consumption during pregnancy and lack of folate supplementation were excluded as covariates due to the low prevalence of self-reported cases (N = 2;2.11% for both).
Statistical analysis
Descriptive statistics were used to compare the main characteristics of the study sample (N = 95) with those of the total sample (N = 561), as well as to compare analyte concentrations in maternal and cord serum within the study sample. Categorical values were reported as absolute counts with percentages (N[%]). Continuous variables were summarized as median and interquartile range [Q1-Q3] for non-normally distributed data, and as means with standard deviation (土SD) for normally distributed variables. Normality was assessed using the Shapiro-Wilk test. Differences between groups were evaluated using the Chi-square test for categorical variables and the Mann-Whitney U-test for non-normally distributed continuous variables. The paired Wilcoxon test was used to assess differences in analyte concentrations between maternal and cord serum. Pearson correlation coefficients were computed to assess linear relationships between analytes or miRNAs in the two groups.
Univariable regression analyses were performed to investigate associations between each miRNA-ΔCt value and single maternal and cord serum exposure concentrations. Multiple comparisons were corrected using the Benjamini-Hochberg false discovery rate (FDR) procedure 41 at an alpha level of 5%. Associations between explanatory variables and chemical exposures were examined using adjusted models, incorporating the covariates described in the previous section. Two separate analyses were conducted: in “Analysis 1”, each chemical was considered individually as the independent variable; in “Analysis 2”, a mixture analysis was performed. For Analysis 1, independent variables included each maternal or cord analyte, as well as the sum of selected PCBs in maternal serum (ΣPCBmat(PCB−74,PCB−118,PCB−138,PCB−153,PCB−180)) and cord serum (ΣPCBcord(PCB−138,PCB−153,PCB−180)), where the included PCBs had less than 10% of substituted values in each side. To assess the mixture effects of EEs (Se, Cu, Zn), POPs (p,p’-DDE, HCB, PCBs), and the total chemical mixture (TOTAL) on miRNA expression, Bayesian Weighted Quantile Sum (BWQS) regression was performed using the R package “bwqs” 42. The Hamiltonian Monte Carlo chain length was set to 1000, with no specified direction for prior distribution; all other parameters were kept at their default settings. Adjusted p-values of < 0.05 were considered significant. We set the chemical of concern weight threshold by dividing 1 by the number of chemicals considered for each analysis. The thresholds are reported in Table S7 of the supplementary materials.
The significant β coefficients estimated from both analyses were reported using a heatmap generated by the “pheatmap” function in the R package. Common associations between bWQS analyses and linear regressions were reported as Venn diagrams 43.
MiRNAs that showed significant associations with chemical exposures were selected for Reactome pathway enrichment analysis to identify key biological pathways which were differentially perturbed between the maternal and cord serum. Validated target genes for each selected miRNA were retrieved from three publicly available miRNA-target interaction databases: miRecords, miRTarBase, and TarBase. This was performed using the multiMir R package 44. To ensure robustness, only target genes that were consistently reported across all the three databases were retained for further analysis. The final list of target genes was grouped based on the side of origin (maternal or cord) and exposure mixture (POPs or EEs), as well as the direction of the associations. These groups were defined by selecting miRNAs that showed statistically significant associations with either EE or POP mixtures in the bWQS analysis, with the aim of identifying enriched functional profiles of each gene cluster using the clusterProfiler R package 45. To reduce pathway redundancy, we first assessed similarities among statistically significant enriched pathways (adjusted p-value < 0.05) using the “pairwise_termsim” function. A similarity threshold of 0.2 was applied to construct a network in which each node represented a pathway, and the edges connected pathways with similarity values exceeding this threshold, with edge weights corresponding to the degree of similarity.
Community detection was then performed using the Louvain algorithm, and each resulting cluster was labelled based on the pathway with the highest degree (i.e., the most connected node within the community). The final results were visualized as a heatmap, highlighting the hub pathway along with its associated subpathways within each community. Only pathways comprising genes from at least four distinct miRNAs were included. For each pathway, we indicated whether the corresponding miRNAs were associated with EE or POP mixtures and specified the direction of the association (positive = red; negative = blue). The heatmap was organized by biological sample origin, with maternal-side results presented in the upper section and cord-side results in the lower section.