Prenatal Environmental Exposure Drives Distinct MicroRNA Patterns in Maternal and Cord Serum: the NEHO cohort
1Institute for Biomedical Research and Innovation, National Research Council of Italy (IRIB-CNR), Palermo, Italy.
Ilaria Cosentini1*, Alessia Li Vigni1*, Noemi Aloi1,2, Alessia Maria Sampino1, Fabio Massimo Contarino3, Paolo Colombo1, Silvia Ruggieri1#, Gaspare Drago1# and Valeria Longo1#
2Department of Biological, Chemical and Pharmaceutical Sciences and Technologies (STEBICEF), University of Palermo, Palermo, Italy.
3Department of Public Health, Siracusa Local Health Authority, Siracusa, Italy.
*Equally contributed to this paper
#Corresponding Authors:
Valeria Longo e-mail: valeria.longo@irib.cnr.it
Gaspare Drago e-mail: gaspare.drago@irib.cnr.it
Silvia Ruggieri e-mail: silvia.ruggieri@irib.cnr.it
ABSTRACT
Environmental exposures during pregnancy can influence cross-talk between the maternal and foetal compartments, eliciting distinct yet coordinated adaptive responses at the placental interface. miRNAs act as key mediators, translating environmental stress into post-transcriptional regulatory signals. This study investigates the effects of prenatal exposure, examining associations between essential elements (EEs), persistent organic pollutants (POPs), and circulating miRNA profiles in 95 maternal–cord pairs from the NEHO birth cohort. Single-chemical regression models identified distinct sets of miRNAs altered by EEs or
POPs in the maternal and cord compartments, while mixture analysis (bWQS) confirmed these associations and revealed both shared and compartment-specific miRNA signatures. Pathway enrichment analysis indicated stress- and apoptosis-related responses in mothers and immune and developmental signaling in foetuses. These findings demonstrate that prenatal exposure to chemical mixtures shapes distinct maternal and foetal miRNA responses, highlighting the placenta’s dual role as a selective barrier and active regulator
of foetal programming.
Key words:
birth cohorts
pollutants exposure
circulating miRNAs
multipollutant analysis
maternal-foetus interface
Abbrevations
bWQS
bayesian Weighted Qantile Sum regression
BMI
Body Mass Index
Ct
Cycle treshold
DDE
Dichlorodiphenyldichloroethylene
DOHaD
Developmental Origins of Health and Disease
EE
Essential Element
FDR
False Discovery Rate
HCB
Hexachlorobenzene
miRNA
microRNA
LOQ
Limit of Quantification
NEHO
Neonatal Environment and Health Outcomes
NPCS
National Priority Contaminated Site
PCB
Poly-Chlorinated Biphenyl
PAH
Polycyclic Aromatic Hydrocarbon
POP
Persistent Organic Pollutant
OC
OrganoChlorine
SNV
Standard normal variate
TNC
Trans-Nonachlor
A
A
A
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.
Results
Characteristics of the study population
A
The NEHO cohort methods are described in detail elsewhere 18. Here, 95 mother/child pairs were randomly selected from within the larger context of the NEHO birth cohort. Table 1 shows a descriptive analysis comparing mothers belonging to the NEHO cohort living in the Priolo NPCS (561 subjects) and the subset of these mothers analysed in this study (N = 95; 16.9%). The analysis of demographic and gestational information (Table 1) shows that the two groups of mothers are homogeneous for age, pre-pregnancy BMI, marital status, educational levels, parity (Primiparous or Multiparous) and gestational length (days). No significant differences were observed between the two groups regarding pregnancy-related habits, such as folate supplementation, smoking, or alcohol consumption.
Focusing on the subpopulation being studied, these mothers had a median age of 31.0 years (IQR 28.0–34.0) and a median pre-pregnancy BMI of 23.1 kg/m² (IQR 20.4–25.1). Smoking during pregnancy was reported by 6 out of 95 mothers (6.32%), while 2 mothers (2.11%) reported alcohol consumption during pregnancy, and an equal proportion (2.11%) declared that they had not taken folic acid. All pregnancies reached full term, and all infants were born alive. The sex distribution of the children was balanced, with 47.4% females and 52.6% males.
Table 1
Descriptive table comparing mothers belonging to the NEHO cohort living in the Priolo NPCS (N = 561) and the subset of these mothers analysed in this study (N = 95).
 
TOTAL
N = 561
SUBSET
N = 95
p-value
Mother’s age (years) - median[Q1-Q3]
31.0
[27.0;34.0]
31.0
[28.0;34.0]
0.467
Mother’s BMI (Kg/m2) - median[Q1-Q3]
22.5
[20.3;25.5]
23.1
[20.4;25.1]
0.410
Marital Status - N(%):
  
0.164
Married
341 (61.0%)
65 (69.1%)
 
Unmarried
218 (39.0%)
29 (30.9%)
 
Educational level - N(%):
  
0.244
Low
158 (28.3%)
19 (20.0%)
 
Medium
277 (49.6%)
53 (55.8%)
 
High
124 (22.2%)
23 (24.2%)
 
Parity - N(%):
  
0.948
Primiparous
246 (52.7%)
47 (51.6%)
 
Multiparous
221 (47.3%)
44 (48.4%)
 
Gestational age (days) - median[Q1-Q3]
275
[270;282]
276
[270;282]
0.970
Folate supplementation - N(%):
  
0.269
no
5 (0.89%)
2 (2.11%)
 
yes
555 (99.1%)
93 (97.9%)
 
Smoking during pregnancy - N(%):
  
0.203
no
496 (88.7%)
89 (93.7%)
 
yes
63 (11.3%)
6 (6.32%)
 
Alcohol consumption during pregnancy – N (%):
  
0.625
no
553 (98.8%)
93 (97.9%)
 
yes
7 (1.25%)
2 (2.11%)
 
Sex – N (%):
  
0.916
Female
259 (46.2%)
45 (47.4%)
 
Male
302 (53.8%)
50 (52.6%)
 
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).
Fig. 1
Statistically significant Pearson correlation indices are reported for analyte concentrations in maternal (A) and cord (B) serum. In (C) the results combining both maternal and cord samples are reported. Positive correlations are shown in red, and negative correlations in blue, with intensity proportional to the strength of the correlation.
Click here to Correct
Analysis of miRNAs in maternal and cord serum
A total of 84 miRNAs were analysed across maternal and cord serum samples. Of these, 75 from maternal serum and 79 from cord serum met quality thresholds (Ct < 40) in more than 70% of samples. The full list of miRNAs, along with their ΔCt values and the percentage of imputed values in each compartment, is reported in Table S2. Five miRNAs were excluded from further analysis in both maternal and cord serum due to having more than 30% imputed values. Additionally, four miRNAs were removed from maternal serum because their threshold Ct values were below 40 in fewer than 70% of subjects. We compared the median ΔCt values between maternal and cord serum to identify statistically significant differences in their expression levels. 31 miRNAs showed significantly lower ΔCt in maternal serum, indicating higher relative expression compared to the cord group. On the other hand, 35 miRNAs showed a higher median ΔCt value in maternal serum, reflecting greater expression in cord serum. The remaining 18 did not show significant differences between compartments.
These results suggest that specific circulating miRNAs are differentially expressed in maternal and cord blood, potentially reflecting foetal-specific gene regulation, maternal-foetal exchange, and exposure-related molecular signaling.
Maternal serum exposures associated with miRNA expression
Associations between individual maternal miRNAs and analytes, correcting for covariates indicated in section 2.6 of “Materials and Methods”, are reported in Fig. 2. Only statistically significant associations after FDR correction (FDR-adjusted p < 0.05) are listed.
Fig. 2
Heatmap showing statistically significant associations (FDR-adjusted p value < 0.05) between the concentrations of each analyte (y-axis) and the ΔCt values of miRNAs (x-axis) in maternal serum. The colour of each box represents the magnitude and direction of the beta coefficient (reported inside the box): red indicates positive associations, while blue indicates negative associations.
Click here to Correct
Maternal concentrations of Cu and Se were positively associated with a higher ΔCt of miR-7-5p and miR-214-3p, respectively. Both miRNAs also showed statistically significant associations with maternal Zn levels. In addition, increasing Zn concentrations were positively associated with the higher ΔCt values of five additional miRNAs. In contrast, miR-30d-5p and miR-374a-5p were inversely associated with Zn concentrations in maternal serum. Regarding POPs, maternal p,p’-DDE and HCB concentrations were associated with increased ΔCt values of miR-1-3p and miR-133a-3p. Although PCB-74 and PCB-118 are grouped in a different chemical class, we identified a similar trend with respect to the OCs. This included positive associations with miR-1-3p and miR-133p. Furthermore, PCB-118 was also associated with miR-195-5p. Higher concentrations of PCB-138 in maternal serum were positively associated with miR-1-3p, miR-195-5p, and miR-30e-5p. The most chlorinated PCBs (PCB-153 and PCB-180) were primarily linked to the latter two of these miRNAs.
A
These patterns were consistent with results obtained using the sum of maternal PCB concentrations: high concentrations of PCBs were positively associated with miR-195-5p and miR-30e-5p. Detailed information about beta estimates are reported in Table S3.
Cord serum exposures associated with miRNA Expression
We investigated associations between miRNA ΔCt values and analyte concentrations in cord serum, correcting for covariates indicated in section 2.6 of “Materials and Methods”. Figure 3 reports the results, focusing on statistically significant associations after FDR correction (FDR-adjusted p < 0.05).
Fig. 3
Heatmap illustrating the associations passing the FDR-adjusted significant threshold (p < 0.05) between analyte concentrations and miRNA ΔCt values in cord serum. The colour gradient reflects the direction and strength of the beta coefficient, with red shades indicating positive associations and blue shades indicating negative associations.
Click here to Correct
Increasing concentrations of Cu and Se in cord serum were positively associated with increasing ΔCt values of miR-30d-5p and miR-7-5p. Additionally, cord Cu concentrations were positively associated with seven miRNAs: miR-195-5p, miR-122-5p, miR-134-5p, miR-200c-3p, miR-574-3p, miR-10a-5p, miR-376-3p. Higher Se concentrations in cord serum were associated with lower ΔCt values of miR-29a-3p and miR-22-3p. Moreover, miR-25-3p and miR-16-5p were inversely associated with Zn concentrations.
Across the POPs, HCB and PCB congeners were associated with distinct miRNAs. Specifically, cord HCB concentrations were associated with increased ΔCt values of let-7c-5p, miR-374a-5p, and miR-191a-5p. Two miRNAs showed consistent associations with all measured PCB congeners in cord serum (PCB-138, PCB-153, and PCB-180) and confirmed their behaviour when analysing the sum of PCB congeners: miR-25-3p was positively associated with each congener, whereas miR-29a-3p was negatively associated with the same congeners. In addition, higher levels of cord PCB-138 and PCB-153 were positively associated with increased ΔCt values of miR-30d-5p and miR-7-5p, they were also negatively associated with miR-133a-3p and miR-133b. These results were confirmed by analyses using the sum of PCBs, while showing no statistically significant results with PCB-180. The β coefficients and confidence intervals for each significant association are reported in Table S4.
Mixture effects of maternal serum chemicals on miRNA Expression
In order to identify mixture effects, we performed bWQS regressions to assess the impact of three different mixtures: EEs (Cu, Se, Zn), POPs (including p,p’-DDE, HCB, PCB-74, PCB-118, PCB-138, PCB-153, PCB-180), and the total mixture comprising both EEs and POPs. The statistically significant associations for each mixture are reported in Fig. 4.
Fig. 4
(A) Heatmap reporting associations between maternal miRNAs and exposure mixtures (EEs, POPs, and TOT) obtained from the bWQS analysis. The colour scale reflects the direction and magnitude of the β coefficients, with red indicating positive and blue indicating negative relationships. Associations also identified in single-analyte regression models are marked with an asterisk “*”. (B-C) Venn diagrams summarizing the overlap between miRNAs associated with each bWQS mixture and those identified in single-exposure analysis. (D) Venn diagram reporting miRNAs commonly associated across the three exposure mixtures.
Click here to Correct
The bWQS analysis of maternal serum EEs identified eighteen significant associations, seven of which confirmed findings from the single-chemical regression models (Fig. 4 panel A, see asterisks, and Fig. 4 panel B). Among these, five emerged with a strong positive association, while two were confirmed as having negative associations. Specifically, five are in common between Se, Cu, and Zn, while two are respectively in common with Cu and Se (Fig. 4 panel B). In addition, eleven new associations were identified (Fig. 4 panels A and B). Five maternal miRNAs were identified as positively associated with the EE mixture, while six additional negative associations were observed.
As for POPs, bWQS analysis confirmed four significant associations previously reported in relation to individual chemicals in maternal serum (see asterisks - Fig. 4 panel A and panel C). Two were previously associated with HCB/DDE, and two were previously associated with at least one PCB. Additionally, nine novel miRNAs which had not previously been associated with single pollutants were identified as significantly associated with the POP mixture (Fig. 4 panel A and C). Of these, six miRNAs exhibited positive associations, while three showed negative associations with the POP mixture. In particular, miR-134-5p was also linked to the EE mixture, displaying a positive association with POPs but a negative association with EEs.
When examining the total mixture, bWQS identified eleven significant associations with maternal miRNA ΔCt values (Fig. 4 panel A and D). Of these, five positive associations overlapped with those observed for the EE mixture. Three positive and two negative associations identified analysing the POP mixture were confirmed in the total mixture. Finally, a novel negative association was observed with miR-296-5p. Detailed information about beta estimates are reported in Table S5.
Mixture effects of cord serum chemicals on miRNA expression
Following the same approach, we applied bWQS regression to examine the impact of three distinct exposure groups: EEs, POPs, and a combined total mixture of both EEs and POPs. Statistically significant associations for each mixture are presented in Fig. 5.
Fig. 5
(A) Associations between cord miRNAs and chemical mixtures (EEs, POPs, and TOT) derived from bWQS analysis. Colours reflect the sign and intensity of the β coefficients (red: positive; blue: negative). Asterisks (“*”) mark miRNAs also identified in single-exposure models. (B–C) Overlaps between bWQS- and single-exposure-associated miRNAs are shown as Venn diagrams. (D) Shared miRNAs across the three mixtures are depicted in a final Venn diagram.
Click here to Correct
Eight out of eleven associations were confirmed through the bWQS analysis performed on the EE mixture in cord serum (Fig. 5 panel A and B). Specifically, the four positively associated miRNAs previously reported for Cu and the two negative associations previously reported for Se were confirmed in the mixture analysis. No associations in common with Zn were observed. In addition, two associations, previously observed in the single-chemical models for both Se and Cu, were confirmed by bWQS analysis (Fig. 5 panel A and B).
Eight cord serum miRNAs were newly identified as being associated with the EE mixture, revealing five positive and three negative associations.
Analysing the POP mixture in cord serum, four associations previously identified in the single-chemical regression were confirmed in the mixture model (Fig. 5 panels and C). Of these, three miRNAs were also negatively associated with at least one PCBs, while the remaining miRNA was positively associated with HCB. In addition, the bWQS analysis of the POP mixture revealed nine novel associations which had not previously been identified in single-chemical regression, with a general negative trend.
The bWQS analysis of the total mixture in cord serum identified fifteen associations (Fig. 5 panel A and D). Among these, one negative association was commonly identified analysing both the EE and POP mixtures. Eight associations overlapped with the EE mixture results, of which six were positive. In addition, two new associations emerged: a positive association with miR-200c-3p and negative association for miR-19b-3p were observed.
The β coefficients and confidence intervals for each significant association are reported in Table S6.
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.
Fig. 6
Mixture component contributions to the outcome. The figure reports maternal information on the left and cord results on the right. In addition, the first line reports details for the EE mixtures, followed by the POP mixture and Total mixture. In each plot, bars correspond to the bottom axis and indicate how many times each chemical in the mixture showed a weight higher than threshold across statistically significant miRNAs. Box plot and circle points correspond to the upper axis (%). Individual data points indicate weights for each chemical. Box plots display the 25th, 50th, and 75th percentiles, where the mean weights are indicated by a rhombus. The dotted black lines indicate the threshold value.
Click here to Correct
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).
Fig. 8
illustrates the main concordances and discrepancies between maternal and cord results obtained through the bWQS analysis. Several miRNAs exhibited consistent associations across the two biological matrices, both in direction and mixture type. For instance, miR-7-5p and miR-122-5p were positively associated with EEs in both maternal and cord serum, whereas miR-128-3p showed a negative association with the POP mixture in both compartments.
Click here to Correct
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.
Fig. 8
Schematic representation of similarity/discrepancy associations between miRNAs and chemical exposure mixtures in maternal and cord compartments. Each miRNA is represented by a pair of adjacent squares (maternal on the left, cord on the right), coloured according to the direction of the association (red = positive; blue = negative). From each square, coloured edges connect to three circular nodes corresponding to the exposure mixtures: EEs, POPs, and TOT (pink-left = maternal mixture, light blue-right = cord mixtures). The colour of each edge matches that of the square, indicating the direction of the association for that specific mixture.
Click here to Correct
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.
Fig. 9
Reactome pathway enrichment analysis. Hub pathways with at least 4 associated miRNAs. Red indicates the presence of miRNAs positively associated in the bWQS analysis, while blue is 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.
Click here to Correct
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.
A
A
A
A
Funding:
This research was funded by the “Progetto Obiettivo Piano Sanitario Nazionale 2019 – Linea 4 – Piano Nazionale Prevenzione della Regione Siciliana (Italy)”.
Acknowledgments: The authors express their gratitude to Mrs. Antonina Azzolina and Mr. Gaetano Spinelli for their assistance
References
1.
Wang Z, Walker GW, Muir DCG, Nagatani-Yoshida K (2020) Toward a Global Understanding of Chemical Pollution: A First Comprehensive Analysis of National and Regional Chemical Inventories. Environ Sci Technol 54:2575–2584
2.
Barker DJ (2007) The origins of the developmental origins theory. J Intern Med 261:412–417
3.
Bouchard MF et al (2011) Prenatal exposure to organophosphate pesticides and IQ in 7-year-old children. Environ Health Perspect 119:1189–1195
4.
Grandjean P et al (2012) Serum vaccine antibody concentrations in children exposed to perfluorinated compounds. JAMA 307:391–397
5.
Ward G, Correia Watts MP, Hansson SR (2025) The unintended consequences of modernity: Pollution and its effect on reproductive, maternal and fetal health. Pregnancy Hypertens 40:101204
6.
Barouki R et al (2018) Epigenetics as a mechanism linking developmental exposures to long-term toxicity. Environ Int 114:77–86
7.
Dieckmann L, Czamara D (2024) Epigenetics of prenatal stress in humans: the current research landscape. Clin Epigenetics 16:20
8.
Goyal D, Limesand SW, Goyal R (2019) Epigenetic responses and the developmental origins of health and disease. J Endocrinol 242:T105–T119
9.
Jiang C et al (2024) Maternal oral exposure to low-dose BPA accelerates the onset of puberty by promoting prepubertal Kiss1 expression in the AVPV nucleus of female offspring. Reprod Toxicol 124:108543
10.
Chaiwangyen W, Khantamat O, Kangwan N, Tipsuwan W, de Sousa FLP (2025) MicroRNA expression in response to environmental hazards: Implications for health. Ecotoxicol Environ Saf 300:118420
11.
Metcalf GAD (2024) MicroRNAs: circulating biomarkers for the early detection of imperceptible cancers via biosensor and machine-learning advances. Oncogene 43:2135–2142
12.
Foley HB et al (2024) EV-miRNA associated with environmental air pollution exposures in the MADRES cohort. Environ Epigenet 10:dvae019
13.
Dietrich B, Haider S, Meinhardt G, Pollheimer J, Knofler M (2022) WNT and NOTCH signaling in human trophoblast development and differentiation. Cell Mol Life Sci 79:292
14.
Longo V et al (2022) A multipollutant low-grade exposure regulates the expression of miR-30b, Let-7a and miR-223 in maternal sera: Evidence from the NEHO cohort. Sci Total Environ 844:157051
15.
Agarwal P, Anand M, Chakraborty P, Singh L, Masih J, Taneja A (2022) Placental levels of polycyclic aromatic hydrocarbons (PAHs) and their association with birth weight of infants. Drug Chem Toxicol 45:868–877
16.
Bove H et al (2019) Ambient black carbon particles reach the fetal side of human placenta. Nat Commun 10:3866
17.
Cosentini I et al (2025) Influence of prenatal hexachlorobenzene, PCB and selenium levels on growth trajectories in the first year of life: Findings from the NEHO birth cohort. Environ Int 195:109225
18.
Ruggieri S et al (2019) Three contaminated sites in southern Italy. The Neonatal Environment and Health Outcomes cohort: protocol for a longitudinal birth cohort study. BMJ Open 9:e029471
19.
Ruggieri S, Maltese S, Drago G, Cibella F, Panunzi S (2021) The Neonatal Environment and Health Outcomes (NEHO) Birth Cohort Study: Behavioral and Socioeconomic Characteristics and Drop-Out Rate from a Longitudinal Birth Cohort in Three Industrially Contaminated Sites in Southern Italy. Int J Environ Res Public Health 18
20.
Brion LP, Heyne R, Lair CS (2021) Role of zinc in neonatal growth and brain growth: review and scoping review. Pediatr Res 89:1627–1640
21.
Willekens J, Runnels LW (2022) Impact of Zinc Transport Mechanisms on Embryonic and Brain Development. Nutrients 14
22.
Jeong Y et al (2018) Placental transfer of persistent organic pollutants and feasibility using the placenta as a non-invasive biomonitoring matrix. Sci Total Environ 612:1498–1505
23.
Addo KA, Palakodety N, Hartwell HJ, Tingare A, Fry RC (2020) Placental microRNAs: Responders to environmental chemicals and mediators of pathophysiology of the human placenta. Toxicol Rep 7:1046–1056
24.
Grzeszczak K, Kwiatkowski S, Kosik-Bogacka D (2020) The Role of Fe, Zn, and Cu in Pregnancy. Biomolecules 10
25.
Zaunschirm-Strutz J, Rieder A, Tocantins C, Diniz MS, Weiss E, Hiden U (2025) MicroRNA profiling in umbilical cord plasma: links to maternal metabolism and neonatal metabolic and inflammatory traits. J Physiol 603:1663–1680
26.
Mitchelson KR, Qin WY (2015) Roles of the canonical myomiRs miR-1, -133 and – 206 in cell development and disease. World J Biol Chem 6:162–208
27.
Ghosh S, Thamotharan S, Fong J, Lei MYY, Janzen C, Devaskar SU (2024) Circulating extracellular vesicular microRNA signatures in early gestation show an association with subsequent clinical features of pre-eclampsia. Sci Rep 14:16770
28.
Timofeeva AV, Fedorov IS, Sukhova YV, Ivanets TY, Sukhikh GT (2023) Prediction of Early- and Late-Onset Pre-Eclampsia in the Preclinical Stage via Placenta-Specific Extracellular miRNA Profiling. Int J Mol Sci 24
29.
Krauskopf J et al (2017) MicroRNA profile for health risk assessment: Environmental exposure to persistent organic pollutants strongly affects the human blood microRNA machinery. Sci Rep 7:9262
30.
Loscalzo G et al (2021) Overexpression of microRNAs miR-25-3p, miR-185-5p and miR-132-3p in Late Onset Fetal Growth Restriction, Validation of Results and Study of the Biochemical Pathways Involved. Int J Mol Sci 23
31.
Fang Z et al (2024) Increased levels of villus-derived exosomal miR-29a-3p in normal pregnancy than uRPL patients suppresses decidual NK cell production of interferon-gamma and exerts a therapeutic effect in abortion-prone mice. Cell Commun Signal 22:230
32.
Yang B et al (2022) Evaluation of Early Biomarkers of Atherosclerosis Associated with Polychlorinated Biphenyl Exposure: An in Vitro and in Vivo Study. Environ Health Perspect 130:37011
33.
Cai JL et al (2016) Polychlorinated biphenyls impair endometrial receptivity in vitro via regulating mir-30d expression and epithelial mesenchymal transition. Toxicology 365:25–34
34.
Jin M, Xu Q, Li J, Xu S, Tang C (2022) Micro-RNAs in Human Placenta: Tiny Molecules, Immense Power. Molecules 27
35.
Zhou D et al (2022) MiR-195-5p facilitates the proliferation, migration, and invasion of human trophoblast cells by targeting FGF2. J Obstet Gynaecol Res 48:2122–2133
36.
Lawless L, Qin Y, Xie L, Zhang K (2023) Trophoblast Differentiation: Mechanisms and Implications for Pregnancy Complications. Nutrients 15
37.
Colamatteo A et al (2023) Immunobiology of pregnancy: from basic science to translational medicine. Trends Mol Med 29:711–725
38.
Duttaroy AK, Basak S (2021) Maternal Fatty Acid Metabolism in Pregnancy and Its Consequences in the Feto-Placental Development. Front Physiol 12:787848
39.
Zona A et al (2019) [SENTIERI: Epidemiological Study of Residents in National Priority Contaminated Sites. Fifth Report]. Epidemiol Prev 43:1–208
40.
Drago G et al (2023) Exposure profiles in pregnant women from a birth cohort in a highly contaminated area of southern Italy. Sci Rep 13:14815
41.
Benjamini Y, Yekutieli D (2001) The control of the false discovery rate in multiple testing under dependency
42.
Colicino E, Pedretti NF, Busgang SA, Gennings C (2020) Per- and poly-fluoroalkyl substances and bone mineral density: Results from the Bayesian weighted quantile sum regression. Environ Epidemiol 4:e092
43.
Wilkinson L (2012) Exact and approximate area-proportional circular Venn and Euler diagrams. IEEE Trans Vis Comput Graph 18:321–331
44.
Ru Y et al (2014) The multiMiR R package and database: integration of microRNA-target interactions along with their disease and drug associations. Nucleic Acids Res 42:e133
45.
Yu G (2024) Thirteen years of clusterProfiler. Innov (Camb) 5:100722
Supplementary Materials
Total words in MS: 11393
Total words in Title: 15
Total words in Abstract: 148
Total Keyword count: 5
Total Images in MS: 9
Total Tables in MS: 11
Total Reference count: 45