Ethics statement
This study was approved by the Ethics Committee of Shijiazhuang Fourth Hospital (Approval No. 20200031). Written informed consent was obtained from all placental tissue donors prior to sample collection. All procedures involving human participants were conducted in accordance with the Declaration of Helsinki. A complete list of reagents, consumables, equipment, and software used in this protocol is provided in the Table of Materials.
Data collection
RNA-seq data were obtained from the GEO database. The datasets included GSE75010, which contains gene expression data from 157 PE placentas and 173 non-PE placentas (N = 330). GSE10588 contains gene expression data from 26 normal placentas and 17 severe PE placentas (N = 43). GSE60438 contains transcriptome profiling data of decidua basalis from pre-eclamptic patients and normotensive pregnancies (N = 125). Single-cell transcriptome data were obtained from the GEO dataset GSE183338. It includes single-nucleus samples from the chorionic villi/maternal-fetal interface of PE and healthy pregnancies.
DEGs analysis
DEGs associated with PE, the GSE75010, GSE10588, and GSE60438 datasets were first preprocessed and normalized. Then, differential analysis was performed using the R package “limma”14 based on the sample grouping information. Genes with p < 0.05 and |log2FC| > 0.5 were selected. Volcano plots of DEGs were generated using the R package ggplot2. Heat maps of the top 20 DEGs were drawn using the R package pheatmap. Subsequently, the intersection of the DEGs selected from the three datasets was taken, and a protein-protein interaction (PPI) network was constructed for the candidate genes using the online platform STRING, with an interaction score ≥0.15. The top 20 hub genes were further identified from this PPI network based on their degree of connectivity, ranked using Cytoscape software. The PPI network results were visualized using Cytoscape software or STRING.
Enrichment analysis
Gene enrichment analysis was performed using the ClusterProfiler and DOSE packages in combination with the Metascape website. The databases were obtained from GO and KEGG. Enrichment analysis was conducted using the “EnrichGO” function. Pathways with p < 0.05 were considered significantly enriched. The enrichment results were visualized using the “ggplot2” and “ggpubr” packages.
Machine learning
To identify robust and biologically meaningful DEGs associated with PE, a multi-model machine learning feature selection analysis was performed based on the publicly available transcriptome dataset GSE60438 (platform: GPL6884). This dataset contains expression profiling of decidua basalis samples collected from pre-eclamptic and normotensive pregnancies at Cesarean section. The pre-filtered DEGs were standardized, and the expression matrix, together with corresponding clinical grouping information, was used as input for four distinct machine learning algorithms in order to reduce model bias and enhance the stability of feature selection.
The four algorithms were applied simultaneously, without a specific order. LASSO was performed using the “glmnet” package to conduct regression analysis and select important feature genes. An L1 regularization term was added to the loss function, which shrinks the coefficients of less important features to zero, thereby achieving feature selection. SVM-RFE was implemented using the “e1071” package to construct a support vector machine with recursive feature elimination. A classifier was first trained using SVM, and the least informative features were iteratively removed based on feature weights, yielding an optimal feature subset. XGBoost was applied using the “xgboost” package to build multiple decision trees. Each tree fitted the residuals of the previous tree, and the weighted outputs were accumulated to obtain the final prediction. Boruta was performed using the “randomForest” package, generating shadow features that competed with real features in training a random forest. Features with importance values significantly higher than random noise were retained.
Single-cell transcriptome data analysis
Single-cell transcriptome data were obtained from the GEO database, and the raw count matrix was retrieved from GSE183338. The count matrix was imported using the “Read10X” function of the Seurat package and converted to a dgCMatrix format. Individual objects were merged into a single aggregate object using the “merge” function, and cell labels were made unique using “RenameCells”. Low-quality cells were filtered based on the following criteria: genes expressed in fewer than three cells were removed, and cells expressing fewer than 200 genes were excluded. Quality-controlled cells were normalized and highly variable genes were identified. Global scaling normalization was applied using “LogNormalize” (scale factor = 10,000), highly variable genes (n = 2,000) were selected using “FindVariableFeatures”, and data were scaled using “ScaleData”. Principal component analysis was performed on highly variable features, and the top 30 principal components were retained. Batch effects between samples were corrected using the Harmony method. Cells were visualized and downscaled using UMAP. Shared nearest neighbor graphs were constructed using “FindNeighbors” and “FindClusters” based on the Louvain algorithm. The resolution parameter in “FindClusters” was optimized between 0.1 and 1. The clustering tree was visualized using the “clustree” function, and a resolution of 0.9 was selected to define cell clusters. Potential doublets were removed using the Scrublet algorithm. Cell clusters were annotated by identifying differentially expressed marker genes using the “FindAllMarkers” function. The non-parametric Wilcoxon rank sum test was applied with Bonferroni correction. Cell identities were assigned based on surface markers, relevant literature, and the Cell Classification Database15.
Cell culture
The trophoblast cell line HTR-8/SVneo cells were cultured in RPMI-1640 medium supplemented with 10% fetal bovine serum and 1% penicillin/streptomycin. Hypoxic conditions were established by culturing cells at 1% O₂, 5% CO₂, and 94% N₂ for 24 h; normoxic controls were maintained at 20% O₂ and 5% CO₂16. All cell culture procedures should be performed in a Class II biosafety cabinet using aseptic technique. Culture media, transfection reagents, and cell waste should be disposed of in accordance with institutional biosafety guidelines.
Cell transfection
Plasmids containing sh-SASH1, sh-NC, OE-SASH1, and OE-NC were synthesized. HTR-8/SVneo cells were seeded at a density of 5 × 105 cells per well in six-well plates. Cells were subsequently transfected with 2 µg sh-SASH1, sh-NC, OE-SASH1, or OE-NC plasmid per well using a transfection reagent according to the manufacturer's instructions. Briefly, plasmid DNA and P3000 Reagent were diluted in Opti-MEM, mixed with Lipofectamine 3000 diluted separately in Opti-MEM, incubated for 15 min at room temperature, and added to cells at 70–80% confluency. Forty-eight hours post-transfection, SASH1 expression was assessed by RT-qPCR and Western blot. The shRNA target sequences used for SASH1 knockdown are listed in Supplementary Table 1.
Real-time quantitative PCR
Total RNA was extracted from HTR-8/SVneo cells and reverse-transcribed into cDNA using a reverse transcription kit at 42 °C for 30 min, followed by 85 °C for 5 min. Real-time quantitative PCR (qPCR) was performed using SYBR Green master mix with the following cycling conditions: 95 °C for 10 min, followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min. The relative mRNA expression was calculated using the ΔΔCt method, with β-actin as the internal reference. Primer sequences used in this experiment are listed in Supplementary Table 2.
Western blot assay
Total protein was extracted from HTR-8/SVneo cells using the lysis buffer. Cell lysates were collected, incubated on ice, and centrifuged at 12,000 × g for 30 min at 4 °C to remove insoluble debris. Protein concentration was determined using a spectrophotometer. Equal amounts of protein (50 µg) were separated by SDS-PAGE and subsequently transferred onto PVDF membranes. The membranes were blocked with 5% non-fat milk and incubated with primary antibodies overnight at 4 °C. After washing, the membranes were incubated with the corresponding secondary antibodies, and protein bands were visualized using an enhanced chemiluminescence detection system.
For protein detection, primary antibodies included anti-SASH1 and β-actin. Appropriate horseradish peroxidase (HRP)-conjugated secondary antibodies—goat anti-rabbit and goat anti-mouse—were employed. β-actin was used as the internal loading control to ensure equal protein loading. The intensity of protein bands was measured and quantified using ImageJ software.
Isolation of P-EXOS
P-EXOS were isolated from placental villous tissue obtained from term placentas of healthy women undergoing elective cesarean section. Placental villous tissue was washed thoroughly with sterile PBS, minced into approximately 1 mm3 fragments, and cultured in RPMI-1640 medium supplemented with 10% exosome-depleted FBS at 37 °C in 5% CO2 for 48 h. The conditioned medium was subjected to differential centrifugation as follows: 300 × g for 10 min to remove cells and tissue debris; 2,000 × g for 20 min to remove cell debris; and 10,000 × g for 30 min to remove microvesicles, all at 4 °C. The resulting supernatant was ultracentrifuged at 120,000 × g for 70 min at 4 °C to pellet exosomes. The pellet was washed once with PBS and re-ultracentrifuged at 120,000 × g for 70 min at 4 °C. The final pellet was resuspended in PBS. The isolated exosomes were characterized by Western blot analysis for exosome markers (PLAP, CD63, and TSG101, with GM130 as a negative control) and further examined by transmission electron microscopy for morphological observation.
P-EXOS Cellular Uptake Experiment
To confirm cellular internalization of P-EXOS, exosomes were fluorescently labeled with the lipophilic membrane dye PKH67 according to the manufacturer's protocol. Briefly, P-EXOS were incubated with PKH67 (4 µM) in Diluent C for 5 min at room temperature, and the reaction was quenched with an equal volume of 1% bovine serum albumin (BSA). Labeled exosomes were re-isolated by ultracentrifugation (120,000 × g, 70 min, 4 °C) to remove unbound dye. PKH67-labeled P-EXOS (50 µg/mL) were then added to HTR-8/SVneo cells and co-incubated for 24 h under normoxic or hypoxic (1% O₂) conditions. Cells were subsequently washed three times with PBS, fixed with 4% paraformaldehyde for 15 min, and the nuclei were counterstained with DAPI (1 µg/mL). Internalization of PKH67-labeled exosomes was visualized by confocal laser scanning microscopy (CLSM; excitation 490 nm, emission 502 nm). For functional co-culture experiments, HTR-8/SVneo cells were treated with P-EXOS at a concentration of 50 µg/mL (protein equivalent) in complete RPMI-1640 medium supplemented with 10% exosome-depleted FBS under hypoxic conditions (1% O₂) for 24 h.
Enzyme-linked immunosorbent assay (ELISA)
Cell culture supernatants were collected, and the levels of IL-6, IL-1β, and TNF-α were measured using IL-6 ELISA kit, IL-1β ELISA kit, and TNF-α ELISA kit, respectively, according to the manufacturers' instructions. Absorbance at 450 nm was measured using a microplate reader, and the actual concentrations were calculated from the standard curves.
TdT-mediated dUTP nick-end labeling (TUNEL)
Apoptotic cells were detected using the TUNEL assay kit according to the manufacturer's instructions. Briefly, cells were fixed with 4% paraformaldehyde for 15 min at room temperature, permeabilized with 0.1% Triton X-100 in PBS for 5 min on ice, and incubated with TUNEL reaction mixture for 60 min at 37 °C in the dark. Nuclei were counterstained with DAPI, and TUNEL-positive cells were visualized using a fluorescence microscope and quantified by counting the percentage of TUNEL-positive cells in at least five randomly selected fields per sample.
Statistical analysis
All data were analyzed using R and GraphPad Prism. Continuous variables are presented as mean ±SD. Two-group comparisons were performed using Student’s t-test, whereas multiple-group comparisons were conducted using one-way ANOVA followed by Tukey’s post hoc test. Statistical significance for categorical variables was assessed by the Chi-square test or Fisher’s exact test. Unless otherwise stated, correlations between molecules were calculated using Spearman correlation analysis. Exosome characterization experiments were performed using P-EXOS isolated from three independent placenta donors. Cell-based experiments were performed in three independent biological replicates, representing independent experiments conducted on separate occasions using HTR-8/SVneo cells of different passages, with each replicate using P-EXOS isolated from a different placental donor. A p < 0.05 was considered statistically significant.