Sexual Dimorphisms of Protein-Coding Gene Profiles in Placentas From Women With Systemic Lupus Erythematosus

Background: Systemic lupus erythematosus (SLE) may cause pathogenic changes in the placentas during human pregnancy, such as decreased placental weight, intraplacental hematoma, ischemic hypoxic change, placental infarction, and decidual vasculopathy, which contribute to high maternal and fetal mortality and morbidity. Sex-specific adaptations of the fetus are associated with SLE pregnancies. The present study aimed to determine the transcriptomic profiles of female and male placentas from women with SLE.

Methods: RNA sequencing (RNA-seq) was performed to identify differentially expressed protein-coding genes (DEGs) in placentas from women with SLE vs. normal term (NT) pregnancies with female and male fetuses (n = 3-5/sex/group). Real-time-quantitative PCR was performed (n = 4 /sex/group) to validate the RNA-seq results. Bioinformatics functional analysis was performed to predict the biological functions and pathways of SLE-dysregulated protein-coding genes.

Results: Compared with NT-female (NT-F) placentas, 119 DEGs were identified in SLE-female (SLE-F) placentas. Among these 119 DEGs, five and zero are located on X- and Y-chromosomes, respectively, and four are located on the mitochondrial genome. Compared with NT-male (NT-M) placentas, 458 DEGs were identified in SLE-male (SLE-M) placentas, among which 16 are located on the X-chromosome and zero on the Y-chromosome and mitochondrial genome. Twenty-four DEGs were commonly dysregulated in SLE-F and -M placentas. Functional analysis showed that SLE-dysregulated protein-coding genes were associated with diverse biological functions and pathways, including angiogenesis, cellular response to growth factor stimulus, heparin-binding, HIF (hypoxia-inducible factor)-1 signaling pathway, and Interleukin-17 (IL-17) signaling pathway in both SLE-F and -M placentas. Biological regulations were differentially enriched between SLE-F and -M placentas. Regulation of blood circulation, response to glucocorticoid, and rhythmic process were all enriched in SLE-F, but not SLE-M placentas. In contrast, tumor necrosis factor production, Th17 cell differentiation, and MDA (melanoma differentiation-associated gene)-5 signaling pathway were enriched in SLE-M but not SLE-F placentas.

Conclusion: This report investigated the protein-coding gene profiles of placenta tissues from SLE patients using RNA-seq. The results suggest that the SLE-dysregulated protein-coding genes in placentas may contribute to the pathophysiological progress of SLE pregnancies in a fetal sex-specific manner, leading to adverse pregnancy outcomes.

Introduction

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease that predominantly affects women of reproductive age, typically causing damages to multiple organ systems (1, 2). SLE pregnancies are associated with maternal complications (e.g., lupus flare, hypertension, preeclampsia, and eclampsia) and fetal complications (e.g., stillbirth, spontaneous abortion, fetal growth restriction, neonatal lupus, and neonatal deaths), and increase the induced abortion rate (3–5).

SLE pregnancies are associated with many placental dysfunctions, (e.g., decreased placental weight, intraplacental hematoma, chronic villitis, thickening of the trophoblast basement membrane, ischemic hypoxic change, placental infarction, decidual vasculopathy, and fetal thrombi) (3, 6). Many mechanisms underlying SLE-induced placental and fetal dysfunction have been proposed. The fetal-placental immune system could in turn, interact with the maternal immune system and mediate maternal immune response (7). For example, anti-DNA antibodies in SLE pregnancies may inhibit trophoblast attachment and migration via cross-reacting with laminin (8). Antiphospholipid antibodies may also alter the placental phospholipid membrane and cause infarctions and edema/swelling (9, 10). In addition, anti-SSA (Sjogren's-syndrome-related antigen A) and anti-SSB (Sjogren's-syndrome-related antigen B) antibodies can cause neonatal lupus and induce fetal injury after crossing the placenta (11). Thus, defects in the placental and fetal responses to autoimmune processes and inflammation are closely related to SLE pregnancy.

Sex-specific adaptations of the fetus have been reported in many complicated pregnancies such as SLE, asthma, and preeclampsia (12–14). Specifically, a lower proportion of male fetuses were born to women with SLE, which may be partially attributed to different chronic inflammation responses between male and female fetuses in early gestation (13). During asthma pregnancy, significantly reduced birth weights were observed in women with female, but not male fetuses (12), suggesting that the inflammatory activities of asthma impacted female, but not male fetuses. Recently, Zhou et al. demonstrated that preeclampsia impaired fetal endothelial function in a fetal sex-dependent manner (14). Given that the placenta is a key organ that closely regulates fetal growth and function, a sex-specific dysregulation of placental growth and function might contribute to the adverse pregnancy outcomes induced by SLE.

Different mechanisms may govern the fetal sex-specific adaptations in various complicated pregnancies. For example, Murphy et al. (12) have reported that female but not male fetuses born to mothers with asthma are associated with a significant increase of maternal circulating monocytes and decreases in placental steroid hydroxylase activity and fetal estriol. Zhou et al. (14) have also reported that fetal-sex specific expression of genes accompanied by preeclampsia-impaired fetal endothelial function. However, mechanisms controlling sex-specific adaptations of SLE pregnancies remain unknown.

In this study, we hypothesized that SLE alters gene expression of placentas, disturbing placental biological functions in a sex-specific manner. We determined the expression profiles of protein-coding genes of placentas from SLE and normal term (NT) pregnancies using RNA sequencing (RNA-seq). Real-time quantitative PCR (RT-qPCR) was conducted to verify RNA-seq results. Functional analysis was conducted to describe the underlying biological functions of differentially expressed protein-coding genes (DEGs) in female and male placentas from SLE pregnancies.

Materials and Methods

Ethical Approval

All procedures were conducted in accordance with the Declaration of Helsinki. Two sets of placental tissues were collected from two hospitals. The first set (SLE, n = 10 with five female and five male fetuses) was collected in Qilu Hospital, Shandong University. The tissue collection protocol was approved and carried out in accordance with the regulation of the Institutional Review Board of Qilu Hospital, Shandong University. SLE was defined according to the American College of Rheumatology classification criteria (15). The SLE disease activity index (SLEDAI) (16) was used to assess the disease activity of SLE patients. SLEDAI scores were assessed within 1 week before delivery. The second set (NT, n = 10 with five female and five male fetuses) was collected in Shanghai First Maternity and Infant Hospital affiliated with Tongji University. The tissue collection protocol was approved and carried out in accordance with the regulation of the Ethical Committee of Shanghai First Maternity and Infant Hospital affiliated with Tongji University. All individuals included in this study were Han Chinese without information on their ancestries. Smokers and patients with cancer or diabetes mellitus were excluded.

RNA Isolation and Quality Control

Placental tissues were obtained within 30 min after vaginal delivery or cesarean section delivery. Placental villi were dissected beneath the chorionic and basal plates (~1 × 1 cm). Placental tissues were snap-frozen in liquid nitrogen and stored at −80°C. Total RNA was isolated from placental tissues using the RNeasy mini kit (Qiagen, Germany). The concentration and quality of RNA samples were assessed using a NanoDrop One spectrophotometer (Thermo Fisher Scientific, USA) and Agilent 2100 Bioanalyzer (Agilent Technologies, USA). Samples with RNA integrity number (RIN) values ≥ 7.0 were used for sequencing.

RNA-Seq and Bioinformatics Analysis of Data

We performed RNA-seq analysis of total RNA samples from placentas (n = 3-5/sex/group; Supplementary Table 1) as described in supplementary methods. RNA-seq strand-specific libraries were constructed using the VAHTS Total RNA-seq (H/M/R) Library Prep Kit (Vazyme, China) following the manufacturer's instructions. Purified libraries were quantified and validated by Qubit 2.0 Fluorometer and Agilent 2100 bioanalyzer to confirm the insert size and calculate the mole concentration. The library construction and sequencing were performed by Sinotech Genomics Co., Ltd Shanghai, China). Overall, more than 66 million reads per sample were generated. Reads with non-canonical letters or low quality and sequences shorter than 25 nucleotides were removed. Reads were trimmed off using FASTQ software (17). Trimmed reads were mapped to the GRCH38 genome using the HISAT2 software (18). DEGs were identified and analyzed using Cuffdiff (19) and R package edgeR (20), respectively. The fold change (FC) was estimated according to the Fragments Per Kilobase of transcript sequence per Millions base pairs sequenced (FPKM) values (21, 22). The P-value significance threshold was set according to the Benjamini-Hochberg false discovery rate (FDR). DEGs were selected using the following criteria: FC > |2| and FDR adjusted P-value (q-value) <0.05 (14), were considered as significantly modulated and recognized as SLE-dysregulated genes. The RNA-seq data have been uploaded and deposited in NCBI (National Center for Biotechnology Information) Gene Expression Omnibus (GEO) (https://www.ncbi.nlm.nih.gov/geo/, accession number: GSE177029).

Functional Analysis of SLE-Dysregulated Genes

The Metascape online analysis tool (https://metascape.org/) was used to predict the biological functions and signaling pathways of the SLE-dysregulated genes (23). Biological terms were selected using the following criteria: P-value <0.01, enrichment factor >1.5, and term count > 3. The most statistically significant (lowest p-value) biological terms within each cluster were chosen to represent the cluster in bar graphs. A subset of representative terms from the cluster was selected and converted as a network plot. Terms with a similarity score (23) > 0.3 were connected by an edge, and the thickness of the edge represents the similarity score. Cytoscape (v3.1.2) was used to visualize the network (24).

RT-PCR

Total RNA (400 ng/sample) was transcribed into complementary DNA (cDNA) using the HiScript II Q RT SuperMix for qPCR (Vazyme Biotech, cat#: R222-01). Diluted cDNA corresponding to 4 ng of original total RNA was utilized as the template in each RT-qPCR reaction. To validate the RNA-seq results, 10 candidate genes were selected for RT-qPCR analysis (25) (n = 4 /sex/group; Supplementary Table 1) using NuHi Robustic SYBR Green Mix. We chose these candidate genes based on the fold changes in SLE vs. NT, relevance to placental function (i.e., angiogenesis and immune responses), expression abundances, and different expression patterns in SLE-F and -M placentas according to RNA-seq data. Primers are listed in Supplementary Table 2. Data were normalized to GAPDH and then analyzed using the 2^−ΔΔCT method (25, 26).

Correlation Analysis of SLE-Dysregulated Genes and SLEDAI Scores

After RT-qPCR verification, relative expression levels of SLE-dysregulated genes were used to analyze the correlation with SLEDAI scores.

Statistical Analyses

Microsoft Excel (2016) for Windows and SigmaPlot (13.0) for Windows were used for statistical analyses. Data were represented as the medians ± standard deviation (SD) or medians with range. Data were analyzed using student's t-test or Mann-Whitney Rank Sum Test when applicable. The relationship between the relative expression levels of genes and SLEDAI scores was analyzed by Pearson's correlation coefficient. P-values <0.05 were considered statistically significant.

Results

Patient Characteristics

All SLE patients received maintenance corticosteroids (prednisone ≤ 15 mg daily) and hydroxychloroquine ( ≤ 400 mg daily) during pregnancy. Three patients in the SLE group were antiphospholipid positive and were on aspirin (100 mg daily) during pregnancy. One SLE patient was positive for anti-SSA antibody. None of the newborns developed neonatal lupus. Patient ages and body mass index (BMI) were similar between NT and SLE pregnancies. However, the newborn body weight for SLE-F was significantly lower than that for NT-F (P = 0.026) (Table 1 and Supplementary Table 1). There was one sample from the SLE group with a growth-restricted male fetus (27). Two patients each in SLE and NT group underwent scheduled Cesarean section delivery. Demographic and clinical characteristics are shown in Table 1 and Supplementary Table 1.

TABLE 1

Table 1. Clinical characteristics.

Distinct Transcriptional Profile in Placentas From SLE Pregnancies

Compared with NT, SLE dysregulated 119 protein-coding genes in female placentas (Figure 1A; Supplementary Table 3), among which five and zero are located on the X- and Y-chromosomes, respectively, and four are located on the mitochondrial genome. Compared with NT, SLE dysregulated 458 protein-coding genes in male placentas (Figure 1A; Supplementary Table 3), among which 16 are located on the X-chromosomes and zero are located on the Y-chromosome or mitochondrial genome (Figure 1A; Supplementary Tables 3, 4).

FIGURE 1

Figure 1. SLE differentially dysregulates transcriptomic profiles of female and male placentas. (A) Circos plot illustrating the chromosome location of differentially expressed protein-coding genes between SLE-F vs. NT-F (pink dots) and SLE-M vs. NT-M (blue dots). Each dot represents one differentially expressed gene. The letters and numbers in the outer layer represent the chromosome location. For the scatter plot tracks, dots outside and inside of the centerline represent upregulated and downregulated genes, respectively. (B,C) Volcano plots showing differentially expressed genes between SLE-F vs. NT-F, and SLE-M vs. NT-M in RNA-seq. Gray dots: no significant difference; pink and green dots: > two-fold upregulation and downregulation, respectively (FDR-adjusted P-value < 0.05) in SLE vs. NT; n = 3-5/group. (D,E) RT-qPCR validation of SLE-dysregulated genes in female and male placentas. ^*P < 0.05 vs. NT, n = 4/group. MT, mitochondrial DNA; SLE, systemic lupus erythematosus; APOLD1, apolipoprotein L domain containing 1; C3, complement C3; CTGF, connective tissue growth factor; DUSP1, dual specificity phosphatase 1; GH2, growth hormone 2; PLAC1, placenta specific 1; FN1, fibronectin 1; PSG4, pregnancy-specific beta-1-glycoprotein 4; VEGFA, vascular endothelial growth factor A; MMP1, matrix metallopeptidase 1. RNA-seq, RNA sequencing; FDR, false discovery rate. Circos plot showing the location DE-genes was generated using circa software for Windows.

Compared to NT-F, 77 and 42 genes were upregulated and downregulated in SLE-F placentas, respectively (Figure 1B; Supplementary Table 3). Compared to NT-M, 438 and 20 genes were upregulated and downregulated in SLE-M placentas, respectively (Figure 1C; Supplementary Table 3). Twenty-four genes were commonly dysregulated in SLE-F and -M placentas (Supplementary Table 3), among which 21were commonly upregulated in SLE-F and -M placentas. One (MATR3) was upregulated in SLE-F placentas but downregulated in SLE-M placentas, while two (CP and LGALS3BP) were upregulated in SLE-M placentas but downregulated in SLE-F placentas (Supplementary Table 3). None of these commonly dysregulated genes is on the X- chromosome, Y-chromosome or mitochondrial genome.

The correlation RT-qPCR and RNA-seq results was performed using the fold change (SLE/NT) of 10 selected genes in RT-qPCR and RNA-seq analyses. RT-qPCR data were correlated significantly with RNA-seq analysis (r = 0.637, P < 0.05) (Figures 1D,E). Specifically, CTGF, VEGFA, and GH2 were upregulated in SLE-F and SLE-M placentas. MMP1 was downregulated in SLE-F but upregulated in SLE-M placentas. APOLD1 was upregulated in SLE-F but not in SLE-M placentas. FN1, PLAC1, and PSG4 were upregulated in SLE-M but not in SLE-F placentas. SLE did not alter C3 and DUSP1 expression in female and male placentas.

Functional and Pathway Analyses of SLE-Dysregulated Genes in Placentas

SLE-dysregulated protein-coding genes were associated with diverse biological functions and pathways (284 for SLE-F placentas and 422 for SLE-M placentas) (Supplementary Tables 5, 6). Sixty-two biological functions and pathways were commonly enriched in both SLE-F and SLE-M placentas (Supplementary Tables 5, 6). These biological functions and pathways included angiogenesis, cellular response to growth factor stimulus, heparin-binding, kidney development, mononuclear cell differentiation, pathways in cancer, response to calcium ion, HIF (hypoxia-inducible factor)-1 signaling pathway, IL-17 signaling pathway, and PI3K-Akt signaling pathway (Figures 2A,C; Table 2, Supplementary Tables 5, 6).

FIGURE 2

Figure 2. Enrichment analysis of dysregulated protein-coding genes in SLE-F and -M placentas. (A) Commonly enriched GOs terms from both SLE-F and -M placentas. (B) Un-commonly enriched GO terms from SLE-F and -M placentas, respectively. (C) Commonly enriched KEGG terms from both SLE-F and -M placentas. (D) Un-commonly enriched KEGG terms from SLE-F and -M placentas, respectively. (E) Network of GO and KEGG enriched terms from SLE-F placentas colored by P-value, where terms containing more dysregulated genes indicate a more significant P-value. (F) Network of GO and KEGG enriched terms from SLE-M placentas colored by P-value, where terms containing more genes tend to have a more significant P-value. SLE, systemic lupus erythematosus; F, female; M, male; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes.

TABLE 2

Table 2. SLE dysregulated biological functions in female and male placentas.

Differential pathway regulations were observed between SLE-F and -M placentas. Regulation of blood circulation, regulation of cellular response to stress, response to glucocorticoid, response to temperature stimulus, response to toxic substance, response to testosterone, rhythmic process, smooth muscle cell proliferation, GnRH (gonadotropin-releasing hormone) signaling pathway, and estrogen signaling pathway were enriched in SLE-F, but not SLE-M placentas (Figures 2B,D; Table 2, Supplementary Tables 5, 6). In contrast, G protein-coupled chemoattractant receptor activity, tumor necrosis factor production, T cell mediated immunity, Th17 cell differentiation, IL-27-mediated signaling pathway, MDA (melanoma differentiation-associated gene)-5 signaling pathway, phospholipase D signaling pathway, transforming growth factor beta receptor signaling pathway, type I interferon signaling pathway, and Wnt signaling pathway were enriched in SLE-M, but not SLE-F placentas (Figures 2B,D; Table 2, Supplementary Tables 5, 6). Figures 2E,F were networks that exhibited the interactions among cluster of genes enriched in biological processes and pathways mentioned above.

Distinct transcriptional profiles in placentas from SLE pregnancies also showed that no significant correlation was found between the relative expression levels of SLE-dysregulated genes and SLEDAI scores (P > 0.05, Supplementary Table 7).

Discussion

To our knowledge, this is the first report that profiles the protein-coding gene expression of human placental tissues from SLE pregnancies using RNA-seq analysis. Overall, more upregulated DEGs were identified than those downregulated. We have further demonstrated that the expression profiles are differentially dysregulated between SLE-F and -M placentas and are associated with differently regulated biological functions and pathways. These data provide clear evidence that SLE differentially regulates the expression of placental protein-coding genes in a fetal sex-dependent manner which may lead to dysregulated placental biological functions.

The mechanisms underlying the fetal sexual dimorphisms of SLE-dysregulated protein-coding gene profiles remain elusive. Expression of DEGs on X-chromosomes is likely to be a major factor that governs these fetal sexual dimorphisms since 4% of DEGs (5 from SLE-F placentas and 16 from SLE-M placentas) are located on the X-chromosome, but no DEGs are on the Y-chromosome. Given that SLE affects women more frequently than men, our current data suggest that the important contribution of X-chromosome-linked genes expression may be associated with the female sex bias in SLE (28). This is in line with a previous study that has shown that a large number of genes that may contribute to the hyperresponsiveness of the female immune system are located on the X chromosome (29). Sex hormones could be another factor that mediates the fetal sexual dimorphisms, as sexual hormones levels are different in the umbilical vein blood of female and male fetuses (30). However, no significant differential expression of androgen or estrogen receptors in SLE-F or -M placentas were detected in this study.

Differential regulation of the mitochondrial genome may also contribute to sexual dimorphisms since our data showed that 4 (MT-ND2, MT-ND3, MT-CYB, and MT-ATP8; Supplementary Tables 3, 4) SLE-dysregulated genes are located on the mitochondrial genome of SLE-F, but not SLE-M placentas, which comprise 30.8% of mitochondrial protein-coding genes. The primary function of mitochondria is to generate the chemical energy to power cellular responses, which is achieved through the electron-transport chain and oxidative phosphorylation (31). The above four SLE-dysregulated mitochondrial genes participate in encoding polypeptides of the oxidative phosphorylation system (31), suggesting that mitochondrial genes might be actively involved in the pathogenesis of SLE placentas. Mitochondrial dysfunction has been demonstrated to be associated with SLE pathogenesis (32–34). The mitochondrial dysfunction may increase not only oxidative stress but also cell apoptosis in SLE patients and defective bioenergetics (33). Oxidative stress induced by mitochondrial dysfunction is considered an essential downstream contributor for SLE pathogenesis (34). A study has shown that CD4+ T cells from an SLE mouse model have higher basal and activated mitochondrial oxidative metabolism, while inhibition of mitochondrial respiratory-chain complex 1 by treating SLE mice with metformin can prevent autoimmune activation (32). Our findings are in line with the above reports and support the notion that the mitochondrial genome of SLE may contribute to fetal sexual dimorphisms, and mitochondrial dysfunction may be present in SLE placentas.

We have also demonstrated that SLE dysregulates biological processes within placentas in a fetal sex-dependent manner. For instance, regulation of blood circulation was enriched in SLE-F but not in SLE-M placentas. In addition, published studies reported that age-standardized cardiovascular disease incidence, prevalence, and mortality rates are lower in women than men (35, 36). Similar results were reported in SLE cohorts that male lupus patients had more cardiovascular damage than their female counterparts (37). Therefore, the enrichment of regulation of blood circulation suggests a possible protective adaptation of cardiovascular diseases in SLE-F placentas.

Our findings revealed that MDA-5 was upregulated in SLE-M, but not SLE-F placentas in association with enrichment of MDA-5 [melanoma differentiation-associated gene-5, gene name: IFIH1 (interferon induced with helicase C domain 1)] signaling pathway only in SLE-M, but not SLE-F placentas. The relationship between MDA-5 and SLE has been proposed by several groups (38, 39). For example, the single nucleotide polymorphisms (SNP) in MDA-5 has been implicated in the pathogenesis of SLE in Africa-Americans, American, Asian, Brazilian, European, and European-Americans (39). MDA-5 has also been demonstrated to be associated with increased sensitivity to serum IFN-α and anti-dsDNA antibodies among SLE patients (38). Our results showed, for the first time, the aberrant expression of MDA-5 in human SLE placentas. Therefore, it is necessary to further investigate the relationship between MDA-5 and SLE pregnancy, especially the function of MDA-5 in SLE placentas.

SLE-F and -M placentas share some common SLE-dysregulated biological processes, e.g., angiogenesis, innate immune responses, and inflammation. This is consistent with the previous observation since these biological processes are essential in placental development and function (40). These data demonstrate that aberrant functions of angiogenesis, innate immune responses, and inflammation may be present in SLE placentas. Specifically, our observations agree with the study that observed a significant increase in angiogenic activity in the serum samples of SLE patients, which was positively associated with SLE disease activity (41, 42). Guilherme et al. reported that serum Vascular Endothelial Growth Factor A (VEGFA) levels were significantly higher in pregnant women with active SLE nephritis than patients with inactive SLE or preeclampsia (42). VEGFA has been demonstrated to be the most potent pro-angiogenic factor and therefore is actively involved in the development of inflammation (43). The gene encoding VEGFA is located on chromosome 6 at 6p21.1, one of the major SLE susceptibility loci (44). Another pro-angiogenic counterpart of VEGFA, which is higher in SLE placentas in this study, is CTGF. CTGF, also known as Cellular Communication Network (CCN) Factor 2, belongs to the CCN family. CTGF gene is located on chromosome 6 at 6q23.2, which is also closely associated with SLE (45). CTGF is a key regulatory and signaling molecule associated with numerous biological processes, including angiogenesis, wound healing, cell proliferation, tissue regeneration, and fibrosis through interaction with many factors, such as VEGFA and TNF-α (46, 47). Based on our RT-qPCR results, VEGFA and CTGF mRNA levels were significantly higher in SLE placentas, suggesting that placental angiogenic activity is disrupted in SLE, just like in preeclampsia (48).

Although defined as an autoimmune disease, SLE is characterized by chronic and acute inflammation conditions in multiple organs with the generation of autoantibodies abnormally produced by one's immune system (49). Torricelli et al. reported that IL-17 was increased in serum from pregnant women with SLE (50). Our RNA-seq analysis also showed that the levels of IL17D, but not other IL17 family members (e.g., IL17A, IL17B, IL17C, and IL17F) was elevated in SLE-M placentas, while the IL-17 signaling pathway was enriched in both male and female placentas, supporting that IL-17 signaling pathway is actively involved in the dysregulation of placental immune response in SLE placentas.

Correlation analysis in this study failed to show any correlation between DEGs and SLEDAI scores. However, a small case number (n = 8) and narrow SLEDAI scores (ranging from 0 to 6) of SLE patients recruited in this study may have prevented us from elucidating any correlations.

In conclusion, this is the first report of protein-coding gene profiles of placentas tissues from SLE pregnancies with female and male fetuses using RNA-seq analysis. The results indicate that the differential expression of protein-coding genes in the female and male placenta may contribute to the different pathogenesis of SLE pregnancies. There are several limitations in this study, including a relatively small sample size which might not address all sex differences in SLE placental gene expression. In addition, SLEDAI scores of the recruited SLE patients were relatively low, ranging from 0 to 6. Thus, while these low SLEDAI scores may provide a valuable resource, recruiting a large cohort of SLE patients with a broader range of SLEDAI scores is needed to establish a more reliable correlation between DEGs and SLEDAI scores.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Ethics Statement

The studies involving human participants were reviewed and approved by the Institutional Review Board of Qilu Hospital, Shandong University and the Scientific and Ethical Committee of Shanghai First Maternity and Infant Hospital affiliated with Tongji University. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

HHL and LTS collected tissues samples and drafted the manuscript. ST designed the tables and figures. YL participated in manuscript preparation. CF designed the figures and read the manuscript critically. KW collected tissues and read the manuscript critically. CZ analyzed the RNA-seq data and read the manuscript critically. JZ and QS conceived the concept and wrote the manuscript. YJZ designed the work and read the manuscript critically. All authors approved the final version of the manuscript.

Funding

This study was supported by ECCM Program of Clinical Research Center of Shandong University (2021SDUCRCB010) and Key Technology Research and Development Program of Shandong (CN) (2018GSF118071 to YJZ, and 2018GSF118025 to HHL).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2022.798907/full#supplementary-material

References

1. D'Cruz DP, Khamashta MA, Hughes GR. Systemic lupus erythematosus. Lancet. (2007) 369:587–96. doi: 10.1016/S0140-6736(07)60279-7

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Liossis SN, Staveri C. What's new in the treatment of systemic lupus erythematosus. Front Med. (2021) 8:655100. doi: 10.3389/fmed.2021.655100

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Magid MS, Kaplan C, Sammaritano LR, Peterson M, Druzin ML, Lockshin MD. Placental pathology in systemic lupus erythematosus: a prospective study. Am J Obstet Gynecol. (1998) 179:226–34. doi: 10.1016/S0002-9378(98)70277-7

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Smyth A, Oliveira GH, Lahr BD, Bailey KR, Norby SM, Garovic VD. A systematic review and meta-analysis of pregnancy outcomes in patients with systemic lupus erythematosus and lupus nephritis. Clin J Am Soc Nephrol. (2010) 5:2060–8. doi: 10.2215/CJN.00240110

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Matucci-Cerinic C, Viglizzo G, Ravelli A, Occella C. Neonatal lupus erythematosus in dizygotic twins with anti-RNP antibodies. Clin Exp Rheumatol. (2021) 39:1446.

PubMed Abstract | Google Scholar

6. Hanly JG, Gladman DD, Rose TH, Laskin CA, Urowitz MB. Lupus pregnancy. A prospective study of placental changes. Arthritis Rheum. (1988) 31:358–66. doi: 10.1002/art.1780310307

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Mor G, Cardenas I. The immune system in pregnancy: a unique complexity. Am J Reprod Immunol. (2010) 63:425–33. doi: 10.1111/j.1600-0897.2010.00836.x

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Qureshi F, Yang Y, Jaques SM, Johnson MP, Naparstek Y, Ulmansky R, et al. Anti-DNA antibodies cross-reacting with laminin inhibit trophoblast attachment and migration: implications for recurrent pregnancy loss in SLE patients. Am J Reprod Immunol. (2000) 44:136–42. doi: 10.1111/j.8755-8920.2000.440302.x

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Nayar R, Lage JM. Placental changes in a first trimester missed abortion in maternal systemic lupus erythematosus with antiphospholipid syndrome; a case report and review of the literature. Hum Pathol. (1996) 27:201–6. doi: 10.1016/S0046-8177(96)90377-9

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Lockshin MD, Sammaritano LR. Lupus pregnancy. Autoimmunity. (2003) 36:33–40. doi: 10.1080/0891693031000067313

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Maidhof W, Hilas O. Lupus: an overview of the disease and management options. P T. (2012) 37:240–9.

PubMed Abstract | Google Scholar

12. Murphy VE, Gibson PG, Giles WB, Zakar T, Smith R, Bisits AM, et al. Maternal asthma is associated with reduced female fetal growth. Am J Respir Crit Care Med. (2003) 168:1317–23. doi: 10.1164/rccm.200303-374OC

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Arkema EV, Askling J, Salmon JE, Simard JF. Brief report: sex ratio of offspring born to women with systemic lupus erythematosus or rheumatoid arthritis. Arthritis Rheumatol. (2017) 69:143–7. doi: 10.1002/art.39843

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Zhou C, Yan Q, Zou QY, Zhong XQ, Tyler CT, Magness RR, et al. Sexual dimorphisms of preeclampsia-dysregulated transcriptomic profiles and cell function in fetal endothelial cells. Hypertension. (2019) 74:154–63. doi: 10.1161/HYPERTENSIONAHA.118.12569

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Hochberg MC. Updating the American College of Rheumatology revised criteria for the classification of systemic lupus erythematosus. Arthritis Rheum. (1997) 40:1725. doi: 10.1002/art.1780400928

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Bombardier C, Gladman DD, Urowitz MB, Caron D, Chang CH. Derivation of the SLEDAI. A disease activity index for lupus patients The Committee on Prognosis Studies in SLE. Arthritis Rheum. (1992) 35:630–40. doi: 10.1002/art.1780350606

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Chen S, Zhou Y, Chen Y, Gu J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. (2018) 34:i884–90. doi: 10.1093/bioinformatics/bty560

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Kim D, Paggi JM, Park C, Bennett C, Salzberg SL. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol. (2019) 37:907–15. doi: 10.1038/s41587-019-0201-4

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. (2012) 7:562–78. doi: 10.1038/nprot.2012.016

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. (2010) 26:139–40. doi: 10.1093/bioinformatics/btp616

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Shao LL, Jiang YH, Jiang LY, Yang CH, Qi YZ. Long non-coding RNA and mRNA analysis of Ang II-induced neuronal dysfunction. Mol Biol Rep. (2019) 46:3233–46. doi: 10.1007/s11033-019-04783-x

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Shboul ZA, Diawara N, Vossough A, Chen JY, Iftekharuddin KM. Joint modeling of RNAseq and radiomics data for glioma molecular characterization and prediction. Front Med. (2021) 8:705071. doi: 10.3389/fmed.2021.705071

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. (2019) 10:1523. doi: 10.1038/s41467-019-09234-6

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Zhou Q, Zhang F, He Z, Zuo MZ. E2F2/5/8 serve as potential prognostic biomarkers and targets for human ovarian cancer. Front Oncol. (2019) 9:161. doi: 10.3389/fonc.2019.00161

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Li Y, Wang K, Zou QY, Jiang YZ, Zhou C, Zheng J. ITE suppresses angiogenic responses in human artery and vein endothelial cells: differential roles of AhR. Reprod Toxicol. (2017) 74:181–8. doi: 10.1016/j.reprotox.2017.09.010

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Chen H, Mu M, Liu Q, Hu H, Tian C, Zhang G, et al. Hepatocyte endoplasmic reticulum stress inhibits hepatitis B virus secretion and delays intracellular hepatitis B virus clearance after entecavir treatment. Front Med. (2020) 7:589040. doi: 10.3389/fmed.2020.589040

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Zhu L, Zhang R, Zhang S, Shi W, Yan W, Wang X, et al. [Chinese neonatal birth weight curve for different gestational age]. Zhonghua Er Ke Za Zhi. (2015) 53:97–103. doi: 10.1353/bkb.2015.0019

CrossRef Full Text | Google Scholar

28. Pyfrom S, Paneru B, Knox JJ, Cancro MP, Posso S, Buckner JH, et al. The dynamic epigenetic regulation of the inactive X chromosome in healthy human B cells is dysregulated in lupus patients. Proc Natl Acad Sci USA. (2021) 118:e2024624118. doi: 10.1073/pnas.2024624118

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Libert C, Dejager L, Pinheiro I. The X chromosome in immune functions: when a chromosome makes the difference. Nat Rev Immunol. (2010) 10:594–604. doi: 10.1038/nri2815

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Furuhashi N, Suzuki M, Fukaya T, Kono H, Shinkawa O, Tachibana Y, et al. Concentrations of luteinizing hormone–human chorionic gonadotropin, beta subunit of human chorionic gonadotropin, follicle-stimulating hormone, estradiol, cortisol, and testosterone in cord sera and their correlations. Am J Obstet Gynecol. (1982) 143:918–21. doi: 10.1016/0002-9378(82)90474-4

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Leonard JV, Schapira AH. Mitochondrial respiratory chain disorders I: mitochondrial DNA defects. Lancet. (2000) 355:299–304. doi: 10.1016/S0140-6736(99)05225-3

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Yin Y, Choi SC, Xu Z, Zeumer L, Kanda N, Croker BP, et al. Glucose oxidation is critical for CD4+ T cell activation in a mouse model of systemic lupus erythematosus. J Immunol. (2016) 196:80–90. doi: 10.4049/jimmunol.1501537

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Lightfoot YL, Blanco LP, Kaplan MJ. Metabolic abnormalities and oxidative stress in lupus. Curr Opin Rheumatol. (2017) 29:442–9. doi: 10.1097/BOR.0000000000000413

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Tsai CY, Hsieh SC, Lu CS, Wu TH, Liao HT, Wu CH, et al. Cross-talk between mitochondrial dysfunction-provoked oxidative stress and aberrant noncoding RNA expression in the pathogenesis and pathophysiology of SLE. Int J Mol Sci. (2019) 20:5183. doi: 10.3390/ijms20205183

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Bots SH, Peters SAE, Woodward M. Sex differences in coronary heart disease and stroke mortality: a global assessment of the effect of ageing between 1980 and 2010. BMJ Glob Health. (2017) 2:e000298. doi: 10.1136/bmjgh-2017-000298

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Roth GA, Johnson C, Abajobir A, Abd-Allah F, Abera SF, Abyu G, et al. Global, regional, and national burden of cardiovascular diseases for 10 causes, 1990 to 2015. J Am Coll Cardiol. (2017) 70:1–25. doi: 10.1016/j.jacc.2017.04.052

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Tan TC, Fang H, Magder LS, Petri MA. Differences between male and female systemic lupus erythematosus in a multiethnic population. J Rheumatol. (2012) 39:759–69. doi: 10.3899/jrheum.111061

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Robinson T, Kariuki SN, Franek BS, Kumabe M, Kumar AA, Badaracco M, et al. Autoimmune disease risk variant of IFIH1 is associated with increased sensitivity to IFN-α and serologic autoimmunity in lupus patients. J Immunol. (2011) 187:1298–303. doi: 10.4049/jimmunol.1100857

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Yin Q, Wu LC, Zheng L, Han MY, Hu LY, Zhao PP, et al. Comprehensive assessment of the association between genes on JAK-STAT pathway (IFIH1, TYK2, IL-10) and systemic lupus erythematosus: a meta-analysis. Arch Dermatol Res. (2018) 310:711–28. doi: 10.1007/s00403-018-1858-0

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Wang K, Zheng J. Signaling regulation of fetoplacental angiogenesis. J Endocrinol. (2012) 212:243–55. doi: 10.1530/JOE-11-0296

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Liu J, Wang X, Yang X, Yan Q, Wang S, Han W. Investigating the role of angiogenesis in systemic lupus erythematosus. Lupus. (2015) 24:621–7. doi: 10.1177/0961203314556293

PubMed Abstract | CrossRef Full Text | Google Scholar

42. de Jesús GR, Lacerda MI, Rodrigues BC, Dos Santos FC, do Nascimento AP, Cristóvão Porto L, et al. Soluble Flt-1, placental growth factor, and vascular endothelial growth factor serum levels to differentiate between active lupus nephritis during pregnancy and preeclampsia. Arthritis Care Res (Hoboken). (2021) 73:717–21. doi: 10.1002/acr.24360

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Kirk SL, Karlik SJ. VEGF and vascular changes in chronic neuroinflammation. J Autoimmun. (2003) 21:353–63. doi: 10.1016/S0896-8411(03)00139-2

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Wakeland EK, Liu K, Graham RR, Behrens TW. Delineating the genetic basis of systemic lupus erythematosus. Immunity. (2001) 15:397–408. doi: 10.1016/S1074-7613(01)00201-1

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Maxwell JR, Gowers IR, Wilson AG. Complex genetic association of 6q23 with autoimmune rheumatic conditions. Arthritis Res Ther. (2009) 11:107. doi: 10.1186/ar2663

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Ramazani Y, Knops N, Elmonem MA, Nguyen TQ, Arcolino FO, van den Heuvel L, et al. Connective tissue growth factor (CTGF) from basics to clinics. Matrix Biol. (2018) 68-69:44–66. doi: 10.1016/j.matbio.2018.03.007

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Chen Z, Zhang N, Chu HY, Yu Y, Zhang ZK, Zhang G, et al. Connective tissue growth factor: from molecular understandings to drug discovery. Front Cell Dev Biol. (2020) 8:593269. doi: 10.3389/fcell.2020.593269

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Chung JY, Song Y, Wang Y, Magness RR, Zheng J. Differential expression of vascular endothelial growth factor (VEGF), endocrine gland derived-VEGF, and VEGF receptors in human placentas from normal and preeclamptic pregnancies. J Clin Endocrinol Metab. (2004) 89:2484–90. doi: 10.1210/jc.2003-031580

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Moulton VR, Suarez-Fueyo A, Meidan E, Li H, Mizui M, Tsokos GC. Pathogenesis of human systemic lupus erythematosus: a cellular perspective. Trends Mol Med. (2017) 23:615–35. doi: 10.1016/j.molmed.2017.05.006

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Torricelli M, Bellisai F, Novembri R, Galeazzi LR, Iuliano A, Voltolini C, et al. High levels of maternal serum IL-17 and activin A in pregnant women affected by systemic lupus erythematosus. Am J Reprod Immunol. (2011) 66:84–9. doi: 10.1111/j.1600-0897.2011.00978.x

PubMed Abstract | CrossRef Full Text | Google Scholar