Skip to main content

Systemic lupus erythematosus dysregulates the expression of long noncoding RNAs in placentas

Abstract

Background

Systemic lupus erythematosus (SLE) can cause placental dysfunctions, which may result in pregnancy complications. Long noncoding RNAs (lncRNAs) are actively involved in the regulation of immune responses during pregnancy. The present study aimed to determine the lncRNA expression profiles in placentas from women with SLE to gain new insights into the underlying molecular mechanisms in SLE pregnancies.

Methods

RNA sequencing (RNA-seq) analysis was performed to identify SLE-dysregulated lncRNAs and mRNAs in placentas from women with SLE and normal full-term (NT) pregnancies. Bioinformatics analysis was conducted to predict the biological functions of these SLE-dysregulated lncRNAs and mRNAs.

Results

RNA-seq analysis identified 52 dysregulated lncRNAs in SLE placentas, including 37 that were upregulated and 15 downregulated. Additional 130 SLE-dysregulated mRNAs were discovered, including 122 upregulated and 8 downregulated. Bioinformatics analysis revealed that SLE-dysregulated genes were associated with biological functions and gene networks, such as regulation of type I interferon-mediated signaling pathway, response to hypoxia, regulation of MAPK (mitogen-activated protein kinase) cascade, response to steroid hormone, complement and coagulation cascades, and Th1 and Th2 cell differentiation.

Conclusions

This is the first report of the lncRNA profiles in placentas from SLE pregnancies. These results suggest that the aberrant expression and the potential regulatory function of lncRNAs in placentas may play comprehensive roles in the pathogenesis of SLE pregnancies. SLE-dysregulated lncRNAs may potentially serve as biomarkers for SLE.

Background

Systemic lupus erythematosus (SLE) is a systemic autoimmune disease that predominantly affects women of reproductive age [1, 2]. During pregnancy, the loss of immune tolerance to the fetus in SLE may cause numerous maternal and fetal complications, including lupus flare, hypertension, preeclampsia (PE), eclampsia, spontaneous abortion, stillbirth, fetal growth restriction (FGR), neonatal lupus, and neonatal deaths [3,4,5,6]. The pregnancy outcomes in mothers with SLE and the well beings of fetuses born to SLE mothers have been improved tremendously over the last five decades due to the development of disease management [7, 8]. However, the mechanism underlying SLE-induced pregnancy complications remains unclear. Therefore, optimal disease control and multidisciplinary obstetrical care throughout gestation are essential to improve pregnancy outcomes in SLE.

SLE can cause placental dysfunction and insufficiencies such as decreased placental weight, ischemic hypoxic change, decidual vasculopathy and thrombi, fetal thrombi, and chronic villitis of implied unknown etiology during pregnancy [5]. These impaired placental functions may lead to numerous maternal and fetal complications, as mentioned above [9, 10].

Long noncoding RNAs (lncRNAs) are transcripts longer than 200 nucleotides in length without recognizable protein-coding potential [11]. LncRNAs play an essential role in regulating gene expression through multiple mechanisms [11,12,13]. LncRNAs also actively regulate key biological processes, including immune cell differentiation and immune responses [14]. Meanwhile, lncRNAs have been related to autoimmune diseases such as SLE, Sjögren’s syndrome, and rheumatoid arthritis [15,16,17,18,19,20,21,22].

LncRNAs are involved in the initiation and development of SLE via different signaling pathways, e.g., the nuclear factor-κB (NF-κB) signaling pathway [19], mitogen-activated protein kinase (MAPK) signaling pathway [20], and tumor necrosis factor (TNF) signaling pathway [16]. LncRNAs are also associated with SLE disease activity [20,21,22]. However, little is known about lncRNA expression profiles and functions in placentas from SLE pregnancy.

This study aims to investigate the SLE-dysregulated lncRNA expression profiles in placentas. Identifying SLE-dysregulated lncRNAs in placentas may aid us in identifying novel therapeutic targets and disease biomarkers for SLE. RNA sequencing (RNA-seq) was performed on placentas from SLE and normal full-term (NT) pregnancies. Bioinformatics analysis was conducted to reveal the underlying biological functions of dysregulated lncRNAs.

Methods

Subjects

All procedures were conducted in accordance with the Declaration of Helsinki. Tissue collection protocols were 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. Pregnant women with SLE (n = 10, with five female and five male fetuses) and NT (n =10, with five female and five male fetuses) were recruited from Qilu Hospital, Shandong University, and Shanghai First Maternity and Infant Hospital affiliated with Tongji University, respectively. Individuals excluded from the study include smokers and patients with cancer or diabetes mellitus. SLE was defined according to the American College of Rheumatology classification criteria [23]. The SLE disease activity index (SLEDAI) [24] was used to evaluate the disease activity of SLE patients. The questionnaire survey was used to collect demographic data and clinical manifestations of all patients. Laboratory data were obtained from the medical record. Clinical manifestations and laboratory records of all subjects are shown in Table 1 and Additional file 1: Table S1.

Table 1 Clinical and laboratory characteristics of the patients in the study

Sampling and RNA isolation

Placentas were collected immediately after C-section or vaginal delivery. Six random biopsies from the fetal side of each placenta were sampled (using a 1 × 1 cm grid), snap-frozen in liquid nitrogen, and stored at −80 °C for further experiments.

Total RNA was isolated using the RNeasy mini kit (Qiagen, Germany). RNA concentration and quality were determined by the Qubit ® 2.0 Fluorometer (Life Technologies, USA) and the Nanodrop One spectrophotometer (Thermo Fisher Scientific Inc., USA). The integrity of total RNA was assessed using the Agilent 2100 Bioanalyzer (Agilent Technologies Inc., USA), and samples with RNA integrity number (RIN) values above 7.0 were used for sequencing.

RNA-seq and data analysis for gene expression

RNA-seq analysis was performed on RNA samples from SLE and NT placentas (n = 8/group with three female and five male fetuses; Additional file 2: Table S2) using the VAHTS Total RNA-seq (H/M/R) Library Prep Kit (Vazyme, China) and Illumina NovaSeq 6000 platform (Illumina, USA) as described in Methods in Data Supplement. Differential expression gene (DEG) analysis for lncRNA/mRNA was performed using R package edgeR [25]. The P-value significance threshold in multiple tests was set by the false discovery rate (FDR) [26]. The fold-changes were estimated according to the fragments per kilobase of transcript sequence per million base pairs sequenced (FPKM) in each sample [27]. A cutoff value of 1 FPKM was used as the detection limit. Differentially expressed RNAs with fold change > |2| and q-value (FDR adjusted P-value) < 0.05, considered significantly modulated, were retained for further analysis. The RNA-seq data have been deposited in Gene Expression Omnibus (GEO) under the accession number GSE177029.

Transcriptomic differences of placentas have been reported in a large cohort of subjects with 94 PE, 56 FGR, and 155 control placental samples using RNA-seq (European genome-phenome archive, https://www.ebi.ac.uk/ega, accession codes: EGAD00001003457, EGAD00001003507, EGAD00001003508, EGAD00001006304, and EGAD00001004860) [28]. We compared the gene expression profiles of normal control placental samples from the current study and the above report. The correlation of FPKM values of the two groups was analyzed. Meanwhile, transcriptomic changes of SLE placentas were also compared with those from the PE and FGR. Principal component analysis [29] was used to predict disparities in expression patterns among SLE-, PE-, and FGR-DEGs.

Functional genomic analysis of dysregulated genes

To explore biological functions and involved signaling pathways of DEGs, the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of the DEGs were conducted using the Metascape analysis tool [30]. A minimum count of 3 was considered positive expression. Genes with enrichment factor > 1.5 and P < 0.05 were considered statistically significant. The GO and KEGG terms fulfilling this condition were defined as significantly enriched GO and KEGG terms.

CIBERSORTx analysis was conducted to estimate the fraction of twenty-two types of immune cells in placental tissues [31]. The non-parametric test (Mann-Whitney U test) was performed to illuminate the different cell components between SLE and normal term placental samples.

Gene co-expression network analysis

To study relationships between dysregulated lncRNAs and mRNAs in SLE placentas, we calculated co-expression relationships between lncRNAs and mRNAs according to the dynamic change of gene expression signal value and obtained the expression regulation relationship and direction between genes to construct the gene expression regulation network. Using the co-expression network, we can analyze the gene regulation ability and obtain the core regulatory genes. The co-expression network was constructed using Cytoscape [32]. Pearson correlation coefficients were used for the lncRNAs-mRNAs co-expression network when they were above 0.90.

Cis- and trans-target gene prediction of lncRNAs

RNAplex was used to identify potential targets of differentially expressed lncRNAs. The lncRNA cis-action is predicted by searching all coding genes within the 10-kb upstream and downstream of the target lncRNA, and these neighboring genes may be regulated by lncRNAs. The lncRNA trans-action was predicted based on nucleic acid base pairing [33]. The GO and KEGG enrichment analyses of the lncRNA target genes were conducted.

Real-time quantitative PCR (RT-qPCR)

To validate the RNA-seq data, six lncRNAs with different expression patterns were selected for RT-qPCR [34] using NuHi Robustic SYBR Green Mix. A new set of samples (n = 8/group with four female and four male fetuses; Additional file 2: Table S2) was used, in which six samples overlapped with ones used in RNA-seq analysis. Data were normalized to GAPDH. Primers are listed in Additional file 3: Table S3. The normalized data were analyzed using the 2−ΔΔCT method [34].

Statistical analyses

Data were expressed as means ± standard deviation or medians with range. Comparison of continuous data was performed by independent Student’s t test. Microsoft Excel (2016) and SigmaPlot (13.0) for Windows were used for statistical analyses. P-values < 0.05 were considered statistically significant.

Results

Patient characteristics

Demographic and clinical characteristics are shown in Table 1 and Additional file 1: Table S1. All SLE patients were on maintenance medication (prednisone ≤ 15 mg daily and hydroxychloroquine 200–400 mg daily). Maternal age and body mass index (BMI) were similar between SLE and NT. The mean newborn body weight and gestation age in SLE were significantly (P < 0.05) lower than those in NT. Eight patients, each in NT and SLE groups, underwent vaginal delivery.

RNA-seq analysis

Deep sequencing generated the mean number of raw reads, 67,745,864 and 66,918,866, for SLE and NT groups, respectively. Removing reads with non-canonical letters or with low quality and discarding the sequences shorter than 25 nucleotides, the mean numbers of clean reads were 66,385,880 and 65,846,828 for SLE and NT groups, respectively, which were retained for analysis.

SLE-dysregulated lncRNAs in placentas

Fifty-two dysregulated lncRNAs were identified (Fig. 1a, b; Additional file 4: Table S4). Of these lncRNAs, 37 were upregulated, and 15 were downregulated. Among these dysregulated lncRNAs, NONHSAT152368.1 (280-fold) was the most upregulated lncRNA, whereas ENST00000602755 (0.0023-fold) was the most downregulated lncRNA. One dysregulated lncRNA (NONHSAT222664.1) was located on X-chromosome and was upregulated. None of these dysregulated lncRNAs was located on Y-chromosome.

Fig. 1
figure 1

SLE-dysregulated transcriptomic profiles of placentas. a Circos plot illustrating the chromosomal position of dysregulated lncRNAs (gray dots) and mRNAs (green dots) between SLE vs. NT. Each dot represents one gene. The numbers and letters in the outer ring indicate the chromosomal location. For each scatter plot track, dots outside and inside of the centerline are upregulated and downregulated genes, respectively. Circos plot was created using Circa program for Windows (OMGenomics). b lncRNA volcano plot of SLE vs. NT. c mRNA volcano plots of SLE vs. NT. n = 8/group

SLE-dysregulated mRNAs in placentas

We also identified 130 dysregulated mRNAs (Fig. 1a, c; Additional file 5: Table S5). Of these dysregulated mRNAs, 122 were upregulated, and 8 were downregulated. KRT24 (246-fold) was the most upregulated mRNA, and ANGPT2 (0.28-fold) was the most downregulated mRNA. One dysregulated mRNA (MAOB) was located on X-chromosome and was upregulated. None of these dysregulated mRNAs was located on Y-chromosome.

RNA-seq data comparisons

After trimming off genes with the FPKM <1, we observed that 9780 (93.3%) protein-coding genes identified in our normal term cohorts (Additional file 6: Table S6) were identical to those reported in Gong’s study [28]. Analysis of FPKM values showed a significant (r = 0.473, P <0.001, Additional file 7: Fig. S1) correlation between these two groups, suggesting similar expression profiles between these two different cohorts. In addition, 17 DEGs were overlapped in SLE (Additional file 8: Fig. S2) and PE placentas, in which 11 were upregulated in SLE but downregulated in PE placentas and one downregulated in SLE but upregulated in PE placentas. Likewise, four DEGs were overlapped between SLE- and FGR-placentas, in which one was upregulated in SLE but downregulated in FGR placentas. Principal component analysis showed disparities in expression patterns among SLE-, PE-, and FGR-DEGs (Additional file 9: Fig. S3). No overlap was identified among SLE-, PE-, and FGR-dysregulated lncRNAs.

Functional genomics analysis of dysregulated genes

We performed bioinformatics analysis with 130 SLE-dysregulated mRNAs, aiming to explore their biological functions and gene networks in SLE placentas. We found that SLE-dysregulated genes were highly enriched in signaling pathways, including regulation of type I interferon-mediated signaling pathway, response to hypoxia, regulation of MAPK (mitogen-activated protein kinase) cascade, response to steroid hormone, complement and coagulation cascades, and Th1 and Th2 cell differentiation (Fig. 2a). In addition, positive or negative regulation of biological functions, such as regulation of vasculature development, cell adhesion, blood coagulation, hemostasis, and response to wounding, was also enriched (Fig. 2b). These findings suggest that a variety of placental function-associated genes/pathways were involved in the pathogenesis of SLE pregnancy.

Fig. 2
figure 2

Enrichment analysis for SLE-dysregulated mRNAs in placentas. a Selected GOs and KEGGs in SLE. b Selected two-way enrichment

CIBERSORTx algorithm was used to estimate the fraction of twenty-two types of immune cells in placentas (Fig. 3 and Additional file 10: Table S7). Macrophages M2 (alternative activated macrophages), T cell CD4 memory resting, monocytes, and neutrophils were top enriched immune cells (Fig. 3a and Additional file 10: Table S7). The results of the Mann-Whitney U test suggested that monocytes were differentially enriched between SLE and NT placentas (Fig. 3b).

Fig. 3
figure 3

CIBERSORTx algorithm analysis. a Bar chart showing the percentage of 22 kinds of immune cells in each sample by estimating relative subsets of RNA transcripts. b Box plot showing the fraction of 22 immune cells in NT and SLE placentas by Mann-Whitney U test. * P < 0.05, # P = 0.05

Co-expression network of lncRNAs and mRNAs

A co-expression network was established to determine the interactions of dysregulated lncRNAs and mRNAs (Additional file 11: Fig. S4). We found a co-expression network composed of 14 lncRNAs and 39 mRNAs. This network showed that one lncRNA (NONHSAT192274.1) could target at least 16 mRNAs, and one mRNA (FN1) could correlate with at least six lncRNAs. NONHSAT192274.1 and NONHSAT112918.2 were the most connected lncRNAs. VEGFA was one of the most connected mRNAs. These results strongly support a network in which lncRNAs and mRNAs function together in placentas during SLE pregnancy.

LncRNA target prediction analysis

The “cis” and “trans” analyses are shown in Additional file 12: Table S8 and Additional file 13: Table S9, respectively. The prediction results indicated that 49 of 52 SLE-dysregulated lncRNAs may target 1429 genes. The two remaining lncRNAs had no potential target genes.

Bioinformatics analysis was conducted to predict the biological roles of target genes of SLE-dysregulated lncRNAs (Fig. 4). The GO project enrichment of target genes included 165 biological processes, 39 cellular components, and seven molecular functions (P < 0.05), mainly involving intracellular protein transport, protein catabolic process, kinase activity, DNA damage checkpoint signaling, and mitochondrial membrane. Hence, lncRNAs may affect these biological processes, cellular components, and functions, contributing to the pathogenesis of SLE.

Fig. 4
figure 4

Enrichment analysis for lncRNA target genes. SLE-dysregulated lncRNAs may target 2387 genes including both “cis” and “trans” genes. Biological functions and involved signaling pathways of these 2387 genes were explored using Metascape online analysis tool

Confirmation of dysregulated lncRNAs

Consistent with RNA-seq data (Fig. 5a), NONHSAT159677.1, NONHSAT198272.1, NONHSAT244274.1, and NONHSAT244275.1 were upregulated, and NONHSAT246155.1 was downregulated (Fig. 5b). The level of NONHSAT209043.1 exhibited a similar upregulation trend as seen in the RNA-seq analysis, but this upregulation did not reach significance (Fig. 5).

Fig. 5
figure 5

Validation of SLE-dysregulated lncRNAs by RT-qPCR. a Relative lncRNA levels (fold of NT control) in SLE. b The change patterns between RNA-seq and RT-qPCR. *Differ (P < 0.05) from NT. n = 8/group

Discussion

In this study, we defined, for the first time, the lncRNA expression profiles of human placentas from SLE pregnancies, though lncRNA expression profiles of human SLE have been reported in peripheral blood mononuclear cells (PBMCs) [17, 18, 20,21,22, 35, 36], monocyte-derived dendritic cells [37], and whole blood [38]. Our RNA-seq analysis revealed 52 dysregulated lncRNAs in SLE placentas with more upregulated genes than downregulated, which agrees with the pattern in monocyte-derived dendritic cells [37] but differs from the patterns reported in PBMCs and whole blood from SLE patients [17, 38]. Among these dysregulated lncRNAs, NONHSAT192272.1 and NONHSAT192274.1 have been reported in whole blood [38], while NONHSAT022132.2 has been found in PBMCs from SLE patients [36]. Interestingly, all the above three genes in SLE placentas show different expression patterns compared with blood samples. For instance, NONHSAT192272.1 and NONHSAT192274.1 are upregulated in SLE placentas whereas downregulated in SLE whole blood. NONHSAT022132.2 is downregulated in SLE placentas but upregulated in SLE PBMCs. This differential dysregulation of lncRNAs might be due to tissue-specific expression of lncRNAs in human tissues [39]. Specifically, the placenta comprises different cell types, including syncytiotrophoblasts, cytotrophoblasts, mesenchymal cells, mesenchymal-derived macrophages, fibroblasts, vascular smooth muscle cells, perivascular cells, endothelial cells, and blood cells in the intervillous space and fetal vessels [40, 41]. The diverse cell types in the placenta tissue may also contribute to the above different profiles.

CIBERSORTx analysis predicted an enrichment of monocytes in SLE placentas, suggesting that SLE robustly dysregulated RNA transcripts in placental monocytes. However, CIBERSORTx, and other analytical tools, such as xCell, TIMER, and EPIC, may not yet be able to predict gene expression profiles of trophoblasts and other placental cells [42], single-cell RNA-seq may do so. SLE and PE pregnancies have similar patterns in placental protein expressions [43]. In this study, we compared our RNA-seq data with the other dataset from a large cohort of subjects with 94 PE, 56 FGR, and 155 control placentas [28]. However, only a small portion of genes are commonly dysregulated (Additional file 8: Fig. S2). Further principal component analysis indicates huge disparities in expression patterns among SLE-, PE-, and FGR-DEGs. Interestingly, one gene (Follistatin Like 3, FSTL3) was upregulated in placentas from SLE, PE, and FGR. FSTL3 is a regulatory glycoprotein that is upregulated in PE and FGR maternal serum and in PE placentas [28, 44]. Hypoxia, which is a biological function enriched in SLE placentas (Fig. 2a), increased the protein levels of FSTL3 in human trophoblasts [45]. In addition, knockdown of FSTL3 suppressed proliferation and migration of human trophoblasts in vitro [45]. Similar to PE and FGR, which are characterized by maternal systemic inflammation [46, 47], SLE is defined as an autoimmune disease with acute and chronic inflammation [48]. Thus, this commonly upregulated gene in placentas from SLE, PE, and FGR may be important in controlling inflammation responses in placentas.

The distribution of SLE-dysregulated mRNAs and lncRNAs on different chromosomes was also observed. One SLE-dysregulated mRNA and one lncRNA are located on X-chromosome, while none is on Y-chromosome, suggesting that sex chromosome may be minimally involved in sex-dimorphic regulation of gene expression in SLE placentas. The X chromosome contains a large number of genes that are partly responsible for the hyperresponsiveness of the female immune system [49], and the Y chromosome contributes to autoimmunity in mouse models of SLE [50, 51].

Our bioinformatics analysis revealed intriguing biological functions of these dysregulated protein-coding genes in SLE placentas. In particular, both positive and negative regulation of biological processes, such as positive regulation of cell adhesion and negative regulation of cell adhesion were enriched (Fig. 2b). The above two-way regulations indicate that SLE placentas are undergoing comprehensive regulations during pregnancy.

It remains to be elusive what causes these dysregulations of lncRNAs. One major factor might be tissue oxygen levels, as the response to hypoxia is recognized as an SLE-dysregulated network (Fig. 2a). This is supported by the current observation that higher expression of VEGFA, which is hypoxia-driven [52]. Enrichment of response to hypoxia in SLE-induced differential expression lncRNAs suggests that placentas from SLE pregnancies may undergo hypoxia compared with controls. Thus, as hypoxia (~1.5% O2) inhibits extravillous trophoblasts outgrowth and proliferation [53], SLE-associated hypoxia may decrease placental weight, leading to FGR [4,5,6]. This is supported by our current finding that the mean newborn body weight from SLE pregnancies is significantly lower than controls.

SLE-dysregulated genes in placentas were also enriched in type I interferon-mediated signaling pathway, which plays a pivotal role in the pathogenesis of SLE [54]. Genome-wide association studies (GWAS) of SLE patients have shown that key genetic variants are involved in over-activation or regulatory deficits in the innate immune responses that are closely correlated to type I interferons (IFNs) [55]. Additionally, risk alleles that operate in IFN pathway genes have also been implicated in the pathogenesis of lupus in GWAS [56]. The placenta is a transient organ. However, its involvement with the type I interferon signaling pathway suggests that the placenta is indeed a target of SLE during pregnancy, which may cause placental dysfunctions, contributing to maternal and fetal complications commonly seen during SLE pregnancy.

Our findings from the lncRNA-mRNA co-expression and lncRNA target prediction strongly support that lncRNAs and mRNAs interact and function together in placentas during SLE pregnancy. This is in line with the reports that similar lncRNA-mRNA co-expression networks in SLE monocyte-derived dendritic cells and PBMCs [17, 37]. Moreover, VEGFA, one of the potent regulators for endothelial growth and function [57], was identified as one of the most connected mRNAs. Given that dysregulated mRNAs in SLE placentas were also enriched in the function of angiogenesis (Fig. 2a), our data suggest that the lncRNA-mRNA co-expression network may be actively involved in the dysregulation of vascular functions during SLE pregnancies.

Dysregulation of lncRNA has been reported to be associated with the severity of SLE. For example, lncRNA NEAT1 levels are significantly higher in PBMCs of SLE patients compared with the healthy group and are positively correlated with SLE disease activity [20]. GAS5 is downregulated in PBMCs and plasma from SLE patients and is negatively correlated with SLE disease activity [21, 58]. Linc0597 and linc0949 are downregulated in PBMCs from SLE patients and negatively correlated with SLE disease activity [21]. In the current study, these lncRNAs (NEAT1, GAS5, linc0597, and linc0949) were not significantly dysregulated in SLE placentas, indicating different lncRNA expression profiles in different SLE tissues. Searching and validating biomarkers for SLE during pregnancy requires a combination of clinical manifestations and potential biomarkers. In the present study, NONHSAT159677.1, NONHSAT198272.1, NONHSAT244274.1, NONHSAT244275.1, and NONHSAT246155.1 were demonstrated to be dysregulated in SLE placentas. Expression profiles of NONHSAT198272.1 and NONHSAT246155.1 in human tissues were reported in the lncRNAs database (http://www.noncode.org/). NONHSAT198272.1 has been observed to be expressed in kidney and placenta tissues, whereas NONHSAT246155.1 was reported to be expressed in various human tissues, and the expression level of NONHSAT246155.1 is the second-highest among these tissues. This is the first report that NONHSAT159677.1, NONHSAT244274.1, and NONHSAT244275.1 are expressed in human placentas. Based on the above observations, we speculate that the above SLE-dysregulated lncRNAs may potentially serve as biomarkers for SLE. However, the functions of the above SLE-dysregulated lncRNAs require further study.

In the current study, we define, for the first time, the transcriptomic profiles in placentas from SLE pregnancies. Our findings support the concept that lncRNAs play comprehensive roles in the pathogenesis of SLE placentas. However, the limitations of our study should be acknowledged. First, this study is limited by the small patient number, and further validation of these dysregulated mRNAs and lncRNAs in a larger cohort of patients is needed to confirm our results. Second, the narrow SLEDAI score range (0–6) of SLE patients recruited strongly limits the analysis of the correlation between DEGs expression levels and disease activity. In addition, further functional studies are needed to define the role of these dysregulated lncRNAs in SLE placentas.

Conclusions

This study demonstrates a comprehensive expression profile of lncRNAs and mRNAs in SLE placentas. The findings suggest regulatory functions of lncRNAs and mRNAs, which are implicated in the development and pathogenesis of SLE pregnancy. SLE-dysregulated lncRNAs may potentially serve as biomarkers for SLE.

Availability of data and materials

The RNA-seq data have been deposited in Gene Expression Omnibus (GEO) under the accession number GSE177029. The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

SLE:

Systemic lupus erythematosus

FGR:

Fetal growth restriction

lncRNAs:

Long noncoding RNAs

NF-κB:

Nuclear factor-κB

MAPK:

Mitogen-activated protein kinase

TNF:

Tumor necrosis factor

RNA-seq:

RNA sequencing

NT:

Normal full-term

SLEDAI:

SLE disease activity index

RIN:

RNA integrity number

DEG:

Differential expression gene

GO:

Gene Ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

FDR:

False discovery rate

FPKM:

Fragments Per Kilobase of transcript sequence per Millions base pairs sequenced

RT-qPCR:

Real-time quantitative PCR

PBMCs:

Peripheral blood mononuclear cells

GWAS:

Genome-wide association studies

IFNs:

Interferons

References

  1. D'Cruz DP, Khamashta MA, Hughes GR. Systemic lupus erythematosus. Lancet (London, England). 2007;369(9561):587–96.

    Article  Google Scholar 

  2. Kudsi M, Nahas LD, Alsawah R, Hamsho A, Omar A. The prevalence of oral mucosal lesions and related factors in systemic lupus erythematosus patients. Arthr Res Ther. 2021;23(1):229.

    CAS  Article  Google Scholar 

  3. Shimada H, Wakiya R, Kanenishi K, Miyatake N, Nakashima S, Mansour MMF, et al. Preterm birth is strongly affected by the glucocorticoid dose during pregnancy in women complicated by systemic lupus erythematosus. Arthr Res Ther. 2022;24(1):10.

    CAS  Article  Google Scholar 

  4. Baer AN, Witter FR, Petri M. Lupus and pregnancy. Obstet Gynecol Surv. 2011;66(10):639–53.

    PubMed  Article  Google Scholar 

  5. 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(1):226–34.

    CAS  PubMed  Article  Google Scholar 

  6. 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(11):2060–8.

    PubMed  PubMed Central  Article  Google Scholar 

  7. Clark CA, Spitzer KA, Laskin CA. Decrease in pregnancy loss rates in patients with systemic lupus erythematosus over a 40-year period. J Rheumatol. 2005;32(9):1709–12.

    PubMed  Google Scholar 

  8. Yan Yuen S, Krizova A, Ouimet JM, Pope JE. Pregnancy outcome in systemic lupus erythematosus (SLE) is improving: results from a case control study and literature review. Open Rheumatol J. 2008;2:89–98.

    PubMed  PubMed Central  Article  Google Scholar 

  9. Gluhovschi C, Gluhovschi G, Petrica L, Velciov S, Gluhovschi A. Pregnancy associated with systemic lupus erythematosus: immune tolerance in pregnancy and its deficiency in systemic lupus erythematosus--an immunological dilemma. J Immunol Res. 2015;2015:241547.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  10. Hayslett JP. Maternal and fetal complications in pregnant women with systemic lupus erythematosus. Am J Kidney Dis. 1991;17(2):123–6.

    CAS  PubMed  Article  Google Scholar 

  11. Ponting CP, Oliver PL, Reik W. Evolution and functions of long noncoding RNAs. Cell. 2009;136(4):629–41.

    CAS  PubMed  Article  Google Scholar 

  12. Wilusz JE, Freier SM, Spector DL. 3' end processing of a long nuclear-retained noncoding RNA yields a tRNA-like cytoplasmic RNA. Cell. 2008;135(5):919–32.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  13. Xiong XD, Ren X, Cai MY, Yang JW, Liu X, Yang JM. Long noncoding RNAs: an emerging powerhouse in the battle between life and death of tumor cells. Drug Resist Updates. 2016;26:28–42.

    Article  Google Scholar 

  14. Gao Y, Li S, Zhang Z, Yu X, Zheng J. The role of long non-coding RNAs in the pathogenesis of RA, SLE, and SS. Front Med. 2018;5:193.

    Article  Google Scholar 

  15. Sandhya P, Joshi K, Scaria V. Long noncoding RNAs could be potential key players in the pathophysiology of Sjögren's syndrome. Int J Rheum Dis. 2015;18(8):898–905.

    CAS  PubMed  Article  Google Scholar 

  16. Zhao CN, Mao YM, Liu LN, Li XM, Wang DG, Pan HF. Emerging role of lncRNAs in systemic lupus erythematosus. Biomed Pharmacother. 2018;106:584–92.

    CAS  PubMed  Article  Google Scholar 

  17. Ye H, Wang X, Wang L, Chu X, Hu X, Sun L, et al. Full high-throughput sequencing analysis of differences in expression profiles of long noncoding RNAs and their mechanisms of action in systemic lupus erythematosus. Arthr Res Ther. 2019;21(1):70.

    Article  Google Scholar 

  18. Luo Q, Li X, Xu C, Zeng L, Ye J, Guo Y, et al. Integrative analysis of long noncoding RNAs and messenger RNA expression profiles in systemic lupus erythematosus. Mol Med Rep. 2018;17(3):3489–96.

    CAS  PubMed  Google Scholar 

  19. Li S, Li C, Zhang J, Tan X, Deng J, Jiang R, et al. Expression profile of long noncoding RNAs in children with systemic lupus erythematosus: a microarray analysis. Clin Exp Rheumatol. 2019;37(1):156–63.

    PubMed  Google Scholar 

  20. Zhang F, Wu L, Qian J, Qu B, Xia S, La T, et al. Identification of the long noncoding RNA NEAT1 as a novel inflammatory regulator acting through MAPK pathway in human lupus. J Autoimmun. 2016;75:96–104.

    CAS  PubMed  Article  Google Scholar 

  21. Li J, Wu GC, Zhang TP, Yang XK, Chen SS, Li LJ, et al. Association of long noncoding RNAs expression levels and their gene polymorphisms with systemic lupus erythematosus. Sci Rep. 2017;7(1):15119.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  22. Wu Y, Zhang F, Ma J, Zhang X, Wu L, Qu B, et al. Association of large intergenic noncoding RNA expression with disease activity and organ damage in systemic lupus erythematosus. Arthr Res Ther. 2015;17(1):131.

    Article  CAS  Google Scholar 

  23. Hochberg MC. Updating the American College of Rheumatology revised criteria for the classification of systemic lupus erythematosus. Arthritis Rheum. 1997;40(9):1725.

    CAS  PubMed  Article  Google Scholar 

  24. 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(6):630–40.

    CAS  PubMed  Article  Google Scholar 

  25. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics (Oxford, England). 2010;26(1):139–40.

    CAS  Article  Google Scholar 

  26. Iqbal J, Tan ZN, Li MX, Chen HB, Ma B, Zhou X, et al. Estradiol alters hippocampal gene expression during the estrous cycle. Endocr Res. 2020;45(2):84–101.

    CAS  PubMed  Article  Google Scholar 

  27. Li CL, Li KC, Wu D, Chen Y, Luo H, Zhao JR, et al. Somatosensory neuron types identified by high-coverage single-cell RNA-sequencing and functional heterogeneity. Cell Res. 2016;26(1):83–102.

    CAS  PubMed  Article  Google Scholar 

  28. Gong S, Gaccioli F, Dopierala J, Sovio U, Cook E, Volders PJ, et al. The RNA landscape of the human placenta in health and disease. Nat Commun. 2021;12(1):2639.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  29. Ringnér M. What is principal component analysis? Nat Biotechnol. 2008;26(3):303–4.

    PubMed  Article  CAS  Google Scholar 

  30. 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(1):1523.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  31. Newman AM, Steen CB, Liu CL, Gentles AJ, Chaudhuri AA, Scherer F, et al. Determining cell type abundance and expression from bulk tissues with digital cytometry. Nat Biotechnol. 2019;37(7):773–82.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  32. Wu Y, Liu F, Luo S, Yin X, He D, Liu J, et al. Co-expression of key gene modules and pathways of human breast cancer cell lines. Biosci Rep. 2019;39(7):BSR20181925.

  33. Song Y, Chen P, Liu P, Bu C, Zhang D. High-Temperature-Responsive Poplar lncRNAs Modulate Target Gene Expression via RNA Interference and Act as RNA Scaffolds to Enhance Heat Tolerance. Int J Mol Sci. 2020;21(18):6808.

  34. Li Y, Wang K, Zou Q-Y, Jiang Y-Z, 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.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  35. Shi L, Zhang Z, Yu AM, Wang W, Wei Z, Akhter E, et al. The SLE transcriptome exhibits evidence of chronic endotoxin exposure and has widespread dysregulation of noncoding and coding RNAs. PLoS One. 2014;9(5):e93846.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  36. Yang H, Liang N, Wang M, Fei Y, Sun J, Li Z, et al. Long noncoding RNA MALAT-1 is a novel inflammatory regulator in human systemic lupus erythematosus. Oncotarget. 2017;8(44):77400–6.

    PubMed  PubMed Central  Article  Google Scholar 

  37. Wang Y, Chen S, Chen S, Du J, Lin J, Qin H, et al. Long noncoding RNA expression profile and association with SLEDAI score in monocyte-derived dendritic cells from patients with systematic lupus erythematosus. Arthr Res Ther. 2018;20(1):138.

    Article  CAS  Google Scholar 

  38. Zhang Q, Liang Y, Yuan H, Li S, Wang JB, Li XM, et al. Integrated analysis of lncRNA, miRNA and mRNA expression profiling in patients with systemic lupus erythematosus. Arch Med Sci. 2019;15(4):872–9.

    CAS  PubMed  Article  Google Scholar 

  39. Washietl S, Kellis M, Garber M. Evolutionary dynamics and tissue specificity of human long noncoding RNAs in six mammals. Genome Res. 2014;24(4):616–28.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  40. Wang Y, Zhao S. Cell types of the placenta. In: Vascular biology of the placenta. San Rafael: Morgan & Claypool Life Sciences, San Rafael (CA); 2010.

    Chapter  Google Scholar 

  41. Jayonta Bhattacharjee FI, Romagnoli R, Bechi N, Caniggia I. Luana Paulesu: ABC transporters in human placenta and their role in maternal-fetal cholesterol transfer: ABCA1 candidate target. In: Zheng J, editor. Recent advances in research on the human placenta: InTech; 2012. p. 335–54.

    Google Scholar 

  42. Aran D, Hu Z, Butte AJ. xCell: digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. 2017;18(1):220.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  43. Matrai CE, Rand JH, Baergen RN. Absence of distinct immunohistochemical distribution of annexin A5, C3b, C4d, and C5b-9 in placentas from patients with antiphospholipid antibodies, preeclampsia, and systemic lupus erythematosus. Pediatr Dev Pathol. 2019;22(5):431–9.

    PubMed  Article  Google Scholar 

  44. Guo J, Tian T, Lu D, Xia G, Wang H, Dong M. Alterations of maternal serum and placental follistatin-like 3 and myostatin in pre-eclampsia. J Obstet Gynaecol Res. 2012;38(7):988–96.

    CAS  PubMed  Article  Google Scholar 

  45. Xie J, Xu Y, Wan L, Wang P, Wang M, Dong M. Involvement of follistatin-like 3 in preeclampsia. Biochem Biophys Res Commun. 2018;506(3):692–7.

    CAS  PubMed  Article  Google Scholar 

  46. Tang L, He G, Liu X, Xu W. Progress in the understanding of the etiology and predictability of fetal growth restriction. Reproduction. 2017;153(6):R227–r240.

    CAS  PubMed  Article  Google Scholar 

  47. 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(1):154–63.

    CAS  PubMed  Article  Google Scholar 

  48. 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(7):615–35.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  49. Libert C, Dejager L, Pinheiro I. The X chromosome in immune functions: when a chromosome makes the difference. Nat Rev Immunol. 2010;10(8):594–604.

    CAS  PubMed  Article  Google Scholar 

  50. Santiago-Raber ML, Kikuchi S, Borel P, Uematsu S, Akira S, Kotzin BL, et al. Evidence for genes in addition to Tlr7 in the Yaa translocation linked with acceleration of systemic lupus erythematosus. J Immunol. 2008;181(2):1556–62.

    CAS  PubMed  Article  Google Scholar 

  51. Subramanian S, Tus K, Li QZ, Wang A, Tian XH, Zhou J, et al. A Tlr7 translocation accelerates systemic autoimmunity in murine lupus. Proc Natl Acad Sci U S A. 2006;103(26):9970–5.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  52. Zhou C, Zou QY, Jiang YZ, Zheng J. Role of oxygen in fetoplacental endothelial responses: hypoxia, physiological normoxia, or hyperoxia? Am J Phys Cell Phys. 2020;318(5):C943–c953.

    CAS  Google Scholar 

  53. James JL, Stone PR, Chamley LW. The effects of oxygen concentration and gestational age on extravillous trophoblast outgrowth in a human first trimester villous explant model. Hum Reprod. 2006;21(10):2699–705.

    PubMed  Article  Google Scholar 

  54. Brohawn PZ, Streicher K, Higgs BW, Morehouse C, Liu H, Illei G, et al. Type I interferon gene signature test-low and -high patients with systemic lupus erythematosus have distinct gene expression signatures. Lupus. 2019;28(13):1524–33.

    CAS  PubMed  Article  Google Scholar 

  55. Sigurdsson S, Nordmark G, Garnier S, Grundberg E, Kwan T, Nilsson O, et al. A risk haplotype of STAT4 for systemic lupus erythematosus is over-expressed, correlates with anti-dsDNA and shows additive effects with two risk alleles of IRF5. Hum Mol Genet. 2008;17(18):2868–76.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  56. Harley IT, Kaufman KM, Langefeld CD, Harley JB, Kelly JA. Genetic susceptibility to SLE: new insights from fine mapping and genome-wide association studies. Nat Rev Genet. 2009;10(5):285–90.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  57. Saleh A, Stathopoulou MG, Dadé S, Ndiaye NC, Azimi-Nezhad M, Murray H, et al. Angiogenesis related genes NOS3, CD14, MMP3 and IL4R are associated to VEGF gene expression and circulating levels in healthy adults. BMC Med Genet. 2015;16:90.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  58. Wu GC, Li J, Leng RX, Li XP, Li XM, Wang DG, et al. Identification of long noncoding RNAs GAS5, linc0597 and lnc-DC in plasma as novel biomarkers for systemic lupus erythematosus. Oncotarget. 2017;8(14):23650–63.

    PubMed  PubMed Central  Article  Google Scholar 

Download references

Acknowledgements

We thank Yi-xin Yin of Shanghai Biotechnology Corporation for technical assistance.

Funding

This study is supported by Shandong provincial key research and development program 2018GSF118071 (YJZ), 2018GSF118025(HHL) and Shandong University Clinical Research Center Key Special Funds for emergency and critical diseases (2021SDUCRCB010).

Author information

Authors and Affiliations

Authors

Contributions

HHL collected tissues and patient information, and drafted the manuscript. LTS collected the patient information, collected tissues, and participated in manuscript preparation. YL and CF participated in manuscript preparation. KW collected tissues and edited the manuscript critically. CZ wrote and edited the manuscript critically. YJZ, JZ, and QS conceived the concept and wrote the manuscript. YJZ designed the work. All authors read and approved the final version of manuscript.

Corresponding authors

Correspondence to Qiang Shu or Ying-jie Zhao.

Ethics declarations

Ethics approval and consent to participate

The study was 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. Informed consent was obtained from all patients at enrollment.

Consent for publication

Written informed consents were obtained from the patients for publication of their individual details and accompanying images in this manuscript. The consent form is held by the authors and is available for review by the Editor-in-Chief.

Competing interests

The authors declare no that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Clinical manifestations and laboratory records.

Additional file 2: Table S2.

Patient demographics for RNA-seq and RT-qPCR.

Additional file 3: Table S3.

RT-qPCR primers.

Additional file 4: Table S4.

SLE dysregulated lncRNAs.

Additional file 5: Table S5.

SLE dysregulated mRNAs.

Additional file 6: Table S6.

Gene expression levels of controls from two separate cohorts.

Additional file 7: Figure S1.

Correlation analysis of FPKM values of two groups.

Additional file 8: Figure S2.

Venn diagram showing the overlap between SLE-, PE-, and FGR-dysregulated protein-coding genes.

Additional file 9: Figure S3.

Principal component analysis of SLE-, PE-, and FGR- dysregulated protein-coding genes.

Additional file 10: Table S7.

CIBERSORTx analysis.

Additional file 11: Figure S4.

Co-expression network with dysregulated lncRNAs and mRNAs.

Additional file 12: Table S8.

The predicted cis-associated genes of the SLE-dysregulated lncRNAs.

Additional file 13: Table S9.

The predicted trans-associated genes of the SLE-dysregulated lncRNAs.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Li, Hh., Sai, Lt., Liu, Y. et al. Systemic lupus erythematosus dysregulates the expression of long noncoding RNAs in placentas. Arthritis Res Ther 24, 142 (2022). https://doi.org/10.1186/s13075-022-02825-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13075-022-02825-7

Keywords

  • lncRNA
  • RNA sequencing
  • Systemic lupus erythematosus
  • Placenta