Identification of serum exosome protein biomarkers of SBA using LC-MS/MS (Label-free)
The proteomics analyses revealed the presence of 397 proteins in the serum exosomes (complete data set in Supplemental Data 2), of which 33 were differentially expressed between the SBA and normal groups. There were 7 proteins extremely upregulated in the SBA group and 12 proteins extremely downregulated in the normal group, while 14 proteins were also at statistically different levels (one upregulated and 13 downregulated; Table 1). GO analysis was conducted to speculate the potential functions of these proteins (except 4 uncharacterized) in three categories (Fig. 1.A). ORA was used to obtain GO enrichment of biological processes with four enrichment processes (Table 2). Accordingly, ACTR2, CORO1A, DNM2, AGT, GLUL, LIMS1, PDCD6, and PKM take part in more than two biological processes.
Table 1
A list of 33 differentially expressed proteins
Protein IDs
|
Protein names
|
Gene names
|
Ratio
|
Upregulated
|
Q6P6T1
|
Complement C1s subcomponent
|
C1s
|
2.99
|
Downregulated
|
F1LZ11
|
Uncharacterized protein
|
N/A
|
0.65
|
Q62636
|
Ras-related protein Rap-1b
|
Rap1b
|
0.65
|
P01015
|
Angiotensinogen
|
Agt
|
0.63
|
D3ZPL2
|
Uncharacterized protein
|
N/A
|
0.59
|
P09606
|
Glutamine synthetase
|
Glul
|
0.55
|
F1M5X4
|
Ig-like domain-containing protein
|
N/A
|
0.55
|
P68136
|
Actin, alpha skeletal muscle
|
Acta1
|
0.54
|
Q8K3U6
|
Coagulation factor VII
|
F7
|
0.53
|
A0A0G2JV65
|
14-3-3 protein zeta/delta
|
Ywhaz
|
0.52
|
B0BNJ1
|
LOC683667 protein
|
Sri
|
0.50
|
A0A0G2K9Z5
|
Uncharacterized protein
|
N/A
|
0.48
|
G3V7W1
|
Programmed cell death protein 6
|
Pdcd6
|
0.45
|
P01883
|
Ig delta chain C region (Fragment)
|
N/A
|
0.09
|
Specific in normal samples
|
A0A0A0MY48
|
Dynamin-2
|
Dnm2
|
|
Q5M7U6
|
Actin-related protein 2
|
Actr2
|
|
Q6P502
|
T-complex protein 1 subunit gamma
|
Cct3
|
|
G3V9N9
|
alpha-1,2-Mannosidase
|
Man1a1
|
|
A0A0G2K393
|
Pleckstrin
|
Plek
|
|
A0A0H2UHM5
|
Protein disulfide-isomerase
|
Pdia3
|
|
C0KUC5
|
LIM and senescent cell antigen-like-containing domain protein
|
Lims1
|
|
Q7M094
|
Destrin-like protein p17a (Fragments)
|
Dstn
|
|
Q5U329
|
Anion exchange protein
|
Slc4a1
|
|
P50115
|
Protein S100-A8
|
S100a8
|
|
Q6LC76
|
Fibronectin (Fragment)
|
Fn1
|
|
Q91ZN1
|
Coronin-1A
|
Coro1a
|
|
Specific in SBA samples
|
A0A0G2K2X4
|
Olfactory receptor
|
LOC100911127
|
|
Q5I0L8
|
Angiopoietin-like 3
|
Angptl3
|
|
M0R7M5
|
Uncharacterized protein
|
LOC100911032
|
|
P01836
|
Ig kappa chain C region, A allele
|
N/A
|
|
P11980
|
Pyruvate kinase PKM
|
Pkm
|
|
Q08420
|
Extracellular superoxide dismutase [Cu-Zn]
|
Sod3
|
|
Q91WX0
|
Complement factor H-related protein
|
RGD1564614
|
|
Table 2
A list of gene ontology analysis of the biological processes associated with the 29 differentially expressed proteins
Gene Set
|
Description
|
Size
|
Expect
|
Ratio
|
P Value
|
FDR
|
GO:0008064
|
regulation of actin polymerization or depolymerization
|
122
|
0.2374
|
21.058
|
3.30E-06
|
0.005217
|
GO:0022603
|
regulation of anatomical structure morphogenesis
|
682
|
1.3273
|
6.7805
|
3.19E-06
|
0.005217
|
GO:1901700
|
response to oxygen-containing compound
|
1542
|
3.0011
|
4.3317
|
1.50E-06
|
0.005217
|
GO:0022607
|
cellular component assembly
|
1997
|
3.8867
|
3.6021
|
4.34E-06
|
0.00563
|
Among these eight proteins, ACTR2,[30, 31] CORO1A,[32] GLUL,[33] and DNM2[34] have been reported to function during neural tube development, but have not been shown to be connected with neural tube defects. Analysis of the four proteins in the String database indicated that ACTR2, CORO1A, and DNM2 may be connected by the same cluster (Fig. 1.B) and require further verification.
Comparative analysis to exclude the effect of atRA on biomarkers in serum exosomes at E18
To clarify whether atRA itself could affect the expression of these biomarkers, we divided into four groups as mentioned above. The expression of ACTR2, CORO1A, and DNM2 in the non-SBA group treated with atRA was slightly decreased compared with that of normal group treated with oil, while their expression in SBA group 1 and SBA group 2 was significantly decreased compared with normal group (Fig. 1.C). The one-way ANOVA analysis with statistical difference intuitively showed this result (Fig. 1.D). The above results indicated that the decreased expression of ACTR2, CORO1A, and DNM2 in serum exosomes were related to the embryogenesis of SBA, and not a result of transcriptional repression by atRA.
Analysis of differential expression of ACTR2, CORO1A, and DNM2 in serum exosomes at E18, E16, E14, and E12
A total of 19 paired E18 serum exosome samples differed from those used in proteomics were examined. The expression of ACTR2, CORO1A, and DNM2 decreased significantly in SBA (Fig. 2.A). Unpaired t-test with statistical difference (ACTR2 P=0.0095, CORO1A P=0.0057, and DNM2 P<0.0001) showed downregulated expression trends in SBA at E18 (Fig. 2.B). To confirm the differential expression of proteins during early gestation, western blotting analysis was performed on serum exosomes at E16 (Fig. 2.C), E14 (Fig. 2.E), and E12 (Fig. 2.G), expressing same trends with E18. Unpaired t-test with statistical difference displayed the same differences at E16 (Fig. 2.D; ACTR2 P=0.0421, CORO1A P=0.0277, and DNM2 P=0.0059), E14 (Fig. 2.F; ACTR2 P=0.0490, CORO1A P=0.0454, and DNM2 P=0.0073) and E12 (Fig. 2.H; ACTR2 P=0.0163, CORO1A P=0.0282, and DNM2 P=0.0157) in expression as those seen at E18. Put together, these results suggest that these proteins maintain a lower expression level in the SBA group from E12 to E18, indicating their diagnostic ability for early gestation.
Analysis of differential expression of ACTR2, CORO1A, and DNM2 in serum without isolation of exosomes
To confirm specific differential expression in the serum exosomes, we quantified their expressions using western blotting in whole serum without exosome isolation. We found no obvious changes between the normal and SBA groups at E18 (Fig. 3.A) and E12 (Fig. 3.C), indicating their downregulation specific to serum exosomes. Coomassie Brilliant Blue staining of total proteins served as the internal reference (calculating the gray value of albumin bands as shown in Fig. 3.A& C by ImageJ) for serum, and staining of the complete gel is given in Supplemental Figure S7& 8. Unpaired t-test showed no statistical difference either at E18 (Fig. 3.B, ACTR2 P=0.5082, CORO1A P=0.7230, DNM2 P=0.1316) or at E12 (Fig. 3.D, ACTR2 P=0.6568, CORO1A P=0.3789, DNM2 P=0.8318).
Analysis of differential expression of ACTR2, CORO1A, and DNM2 in fetal neural tubes and FNEs
To verify whether the downregulated expression of these proteins in the serum exosomes was related to the NTDs, we quantified their expressions using western blotting in SBA neural tubes at E18 and E12. We found that their expression patterns were the same as in the serum exosomes with significantly decreased SBA expression (Fig. 3.E, G). Unpaired t-test with statistical difference showed downregulated expression trends in SBA both at E18 (Fig. 3.F; ACTR2 P=0.0012, CORO1A P=0.0264, and DNM2 P=0.0151) and E12 (Fig. 3.H; ACTR2 P=0.0301, CORO1A P=0.0247, and DNM2 P=0.0012).
Immunohistochemistry analysis was performed in E18 embryos to evaluate the location and level of protein expression in neural tubes. ACTR2 was expressed in all cells within the field of vision with increased expression in neuroepithelium, neural cells, and neural fibers (Fig. 4.A) and downregulated expression in SBA (Fig. 4.D). CORO1A was expressed increasingly in nerve fibers (Fig. 4.B) but downregulated in SBA (Fig. 4.E). DNM2 was expressed increasingly in neuroepithelium and neural cells (Fig. 4.C), but hardly detected in SBA (Fig. 4.F).
To further verify whether the downregulated expression of ACTR2, CORO1A, and DNM2 in total serum exosomes originated from fetal neural source, we tested the proteins in FNEs isolated from total serum exosomes of maternal rats (Fig. 4.G). The change inexpression of CORO1A and DNM2 in FNEs followed the same pattern as that in total serum exosomes, which showed downregulated expression in SBA (Fig. 4.H; CORO1A P=0.0487, and DNM2 P=0.0316), while there was no significant difference in the expression of ACTR2 between normal samples and SBA serum exosome samples (Fig. 4.H; ACTR2 P=0.4987). These suggested that CORO1A and DNM2 expression was specific to neural cells and tissues.
Analysis of differential expression of CORO1A and DNM2 in serum exosomes of pregnant women
We performed ELISA on DNM2 and CORO1A in maternal serum exosomes to validate their expression trends. Samples from pregnant women were paired by gestational week, and paired t-test was used for statistical analysis. DNM2 expression decreased significantly in all 19 pregnant women who were diagnosed as carrying fetuses with NTDs compared with normal controls (Fig. 5.A; P=0.0011) and its ROC curve showed high accuracy (specificity=78.95%, sensitivity=73.68%, AUC=0.806) (Fig. 5.B). Further analysis by the subgroups showed statistical difference of DNM2 expression both in SBA group (Fig. 5.G; P=0.0035), and anencephalus and exencephalus group (Fig. 5.H; P=0.0058). Subsequently, we focused on the expression of DNM2 during different pregnancy periods, which included two groups for statistical analysis: gestational weeks 12–18 and 19–40. There was downregulated expression with significant difference in the group of 12–18 gestational weeks (Fig. 5.C; P=0.0034), which showed extremely high accuracy (specificity=100%, sensitivity=100%, AUC=1.000) (Fig. 5.D). There was also downregulated expression, but with no significant difference, in the group of 19–40 gestational weeks (Fig. 5.E; P=0.0557), which showed lower accuracy (specificity=58.33%, sensitivity=83.33%, AUC=0.736) (Fig. 5.F).
CORO1A expression was also showed decreased significantly in the 19 pregnant women compared with that of normal controls (Fig. 5.I; P=0.0022), and its ROC curve showed high accuracy (specificity=89.47%, sensitivity=68.42%, AUC=0.817) (Fig. 5.J). Further analysis by the subgroups showed statistical difference of CORO1A expression both in SBA group (Fig. 5.O; P<0.0001), and in anencephalus and exencephalus group (Fig. 5.P; P=0.0186). The results showed that there was obvious downregulated expression with significant difference in the 12–18 gestational weeks group (Fig. 5.K; P=0.0108), which showed high accuracy (specificity=85.71%, sensitivity=85.71%, AUC=0.857) (Fig. 5.L). There was also downregulated expression, but with no significant difference in the 19–40 gestational weeks group (Fig. 5.M; P=0.0599), which showed lower accuracy (specificity=83.33%, sensitivity=75%, AUC= 0.785) (Fig. 5.N).
We further analyzed the combined performance of DNM2 and CORO1A in maternal serum exosomes using ELISA. The ROC curve showed better performance in all 19 pair samples and it showed higher accuracy (specificity= 78.95%, sensitivity= 94.74%, AUC=0.889) (Fig. 5.Q). In the 12–18 gestational weeks group, ROC curve also showed better performance and extreme accuracy (specificity= 100%, sensitivity= 100%, AUC=1.000) (Fig. 5.R). In the 19–40 gestational weeks group (n=12), ROC curve also showed better performance (specificity= 66.67%, sensitivity= 100%, AUC=0.854) (Fig. 5.S). In addition, the ROC curve showed good performance in 13 pair SBA samples and high accuracy and specificity (specificity= 92.31%, sensitivity= 76.92%, AUC=0.888) (Fig. 5.T).