Int J Med Sci 2022; 19(1):1-12. doi:10.7150/ijms.65002 This issue

Research Paper

Expression profiling of lncRNAs and mRNAs in placental site trophoblastic tumor (PSTT) by microarray

Jianfeng Gan1,2, Zhixian Chen1,2, Xuan Feng1,2, Zhi Wei1,2, Sai Zhang1,2, Yan Du1,2, Congjian Xu1,2 Corresponding address, Hongbo Zhao1,2 Corresponding address

1. Shanghai Key Laboratory of Female Reproductive Endocrine Related Diseases, Obstetrics and Gynecology Hospital, Fudan University, Shanghai 200011, People's Republic of China.
2. Department of Obstetrics and Gynecology of Shanghai Medical School, Fudan University, Shanghai 200032, People's Republic of China.

This is an open access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/). See http://ivyspring.com/terms for full terms and conditions.
Citation:
Gan J, Chen Z, Feng X, Wei Z, Zhang S, Du Y, Xu C, Zhao H. Expression profiling of lncRNAs and mRNAs in placental site trophoblastic tumor (PSTT) by microarray. Int J Med Sci 2022; 19(1):1-12. doi:10.7150/ijms.65002. Available from https://www.medsci.org/v19p0001.htm

File import instruction

Abstract

Graphic abstract

As a rare type of gestational trophoblastic disease, placental site trophoblastic tumor (PSTT) is originated from intermediate trophoblast cells. Long noncoding RNAs (lncRNAs) regulate numerous biological process. However, the role of lncRNAs in PSTT remains poorly understood. In the present study, expression levels of lncRNAs and mRNAs in four human PSTT tissues and four normal placental villi were investigated. The results of microarray were validated by the reverse transcription and quantitative real-time polymerase reaction (RT-qPCR) and immunohistochemistry analyses. Furthermore, GO and KEGG pathway analyses were performed to identify the underlying biological processes and signaling pathways of aberrantly expressed lncRNAs and mRNAs. We also conducted the coding-non-coding gene co-expression (CNC) network to explore the interaction of altered lncRNAs and mRNAs. In total, we identified 1247 up-regulated lncRNAs and 1013 down-regulated lncRNAs as well as 828 up-regulated mRNAs and 1393 down-regulated mRNAs in PSTT tissues compared to normal villi (fold change ≥ 2.0, p < 0.05). GO analysis showed that mitochondrion was the most significantly down-regulated GO term, and immune response was the most significantly up-regulated term. A CNC network profile based on six confirmed lncRNAs (NONHSAT114519, NR_103711, NONHSAT003875, NONHSAT136587, NONHSAT134431, NONHSAT102500) as well as 354 mRNAs was composed of 497 edges. GO and KEGG analyses indicated that interacted mRNAs were enriched in the signal-recognition particle (SRP)-dependent cotranslational protein targeting to membrane and Ribosome pathway. It contributes to expand the understanding of the aberrant lncRNAs and mRNAs profiles of PSTT, which may be helpful for the exploration of new diagnosis and treatment of PSTT.

Keywords: placental site trophoblastic tumor (PSTT), microarray analysis, lncRNA

Introduction

Placental site trophoblastic tumor (PSTT) is a rare type of gestational trophoblastic neoplasms (GTN), which was first described in 1976[1]. Compared with other types of GTN, PSTT tends to have a more unpredictable biological behavior, chemotherapy resistance and poor prognosis [2]. Therefore, hysterectomy is the first-consideration for PSTT patients who may face the fertility problems after the operations [3]. The proportion of PSTT in GTN is about 0.2-3% with an estimated incidence of 1/100000 pregnancies [4-10]. The early diagnosis of PSTT remains unclear, and distinguishing the benign and malignant forms of PSTT in the early stage is important to determine next clinical therapy [11]. Therefore, it is in urgent need for a better comprehension of PSTT to explore new diagnostic and therapeutic targets.

Long non-coding RNAs (lncRNAs) are non-protein-coding RNAs that longer than 200 nucleotides, which have been confirmed to function as regulators of many cellular processes, such as development, differentiation, cell fate and disease pathogenesis as well as tumorigenesis [12-14]. It has been reported that lncRNAs may play important roles in chemoresistance in many malignant tumors including ovarian cancer [15], breast cancer [16], gastric cancer [17], glioblastoma [18]. For example, LINC00261 is abnormally expressed in a number of tumors including pancreatic cancer, gastric cancer, choriocarcinoma and so on [19]. It has been confirmed to alleviate cisplatin resistance in colon-cancer via Wnt/β-catenin pathway [20]. Nevertheless, little is known about the lncRNA and mRNA expression profiles of PSTT and their involvement in the progression of PSTT. The present study performed the first microarray analysis of lncRNA and mRNA expression profiles of PSTT, which may contribute to expand our understanding of PSTT as well as explore more effective treatment strategy.

Materials and methods

Tissues and ethics statement

In total, 4 human first-trimester placental villi, 18 PSTT tissues were obtained from the tissue bank of Obstetrics and Gynecology Hospital of Fudan University. Written informed consent was obtained. The pathological diagnosis of choriocarcinoma and PSTT was confirmed in the department of pathology in Obstetrics and Gynecology Hospital. Human first-trimester placental villi (6-7 weeks of gestation) were collected from 25-30 years old women without underlying health conditions that had normal pregnancies and terminated for non-medical reasons. PSTT tissues were collected from surgeries for primary tumors. Patients were 25.5-38 years old with at least 1 gravidity and 1 parity. And the stage of PSTT is FIGO Ⅰ. All the PSTT tissues were affirmed by frozen sections to contain more than 75% lesions. This study was approved by the Institution Ethics Committee of Obstetrics and Gynecology Hospital of Fudan University.

Microarray analysis

Total RNAs were isolated by using Trizol RNA extraction kit (Invitrogen, Carlsbad,CA, USA) and was quantified by the NanoDrop ND-2000 (Thermo Scientific), and the RNAs integrity and 28S/18S were assessed using Agilent Bioanalyzer 2100 (Agilent Technologies). The sample labeling, microarray hybridization and washing were performed based on the manufacturer's standard protocols. Briefly, total RNAs were transcribed to double strand cDNAs and then synthesized cRNAs. Next, 2nd cycle cDNAs were synthesized from cRNAs. Followed fragmentation and biotin labeling, the 2nd cycle cDNAs were hybridized onto the microarray. After washing and staining, the arrays were scanned by the Affymetrix Scanner 3000 (Affymetrix). The sample with RIN (RNA integrity number) values of RIN ≥ 7.0, 28S/18S ≥ 0.7 and the total amount can meet 2 or more experiments were processed for subsequent experiments. In this case, 4 human first-trimester placental villi and 4 PSTT tissues were selected and subjected to the Affymetrix Human OElncRNA Array.

Reverse transcription and quantitative real-time polymerase chain reaction (RT-qPCR)

Quantification was performed with a two-step reaction process: reverse transcription (RT) and quantitative PCR. Firstly, mRNA was transcribed to cDNA in a GeneAmp® PCR System 9700 (Applied Biosystems, USA). The second step was performed in a GeneAmp® PCR System 9700 (Applied Biosystems, USA) with 5 × HiScript II Q RT SuperMix IIa. Quantitative Real-time PCR was performed using LightCycler® 480 Ⅱ Real-time PCR Instrument (Roche, Swiss). The primer sequences were designed and synthesized by Generay Biotech (Generay, PRC) based on the mRNA sequences obtained from the NCBI database as shown in Table 1.

Immunohistochemistry

The first-trimester placental villi, choriocarcinoma tissues, and PSTT tissues sections were immunohistochemically stained by using primary antibodies of PLAC8, EGR1 and ADAMTS6 at 4°C overnight followed by the secondary antibody for 1 h at 37°C. The band was then visualized using the EnVision Detection Systems (Dako, Glostrup, Denmark). Paraffin-embedded tissue sections were stained with hematoxylin-eosin (HE) on Leica automated Staining/Coverslipping Workstations (Leica Biosystems, Nussloch, Germany). PLAC8 antibody (12284- 1-AP) for immunohistochemistry (IHC) was purchased from Proteintech (Chicago, IL, USA). EGR1 antibody (4154S) was purchased from Cell Signaling Technology (San Diego, CA, USA). ADAMTS6 antibody (PA5-60365) was obtained from Thermo Fisher Scientific (Waltham, MA, USA). HRP- conjugated secondary antibodies were purchased from Jackson ImmunoResearch Laboratories (West Grove, PA, USA).

Gene Oncology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and coding-non-coding gene co-expression (CNC) network

GO analysis and KEGG analysis were applied to determine the potential biological functions and significant pathways of altered mRNAs or lncRNA-interacted associated mRNAs. Afterwards, we established a CNC network between validated 6 lncRNAs and interacted mRNAs to explore their relationships. Pearson correlation coefficients (PCCs) ≥ 0.95 were selected as the baseline of correlation analysis. Cytoscape V3.8.0 (The Cytoscape Consortium, San Diego, CA, USA) was used to represent the interaction in pictures.

Gene set enrichment analysis (GSEA)

GSEA was conducted using the Molecular Signatures Database hallmark gene set collection (v7.4) to identify differences in pathways among choriocarcinoma (CC), PSTT, epithelioid trophoblastic tumor (ETT) and normal villi. Data set GSE135727 [21] was downloaded from Gene Expression Omnibus [22].

Statistical analysis

Affymetrix GeneChip Command Console (version 4.0, Affymetrix) software was used to extract raw data. Next, Expression Console (version1.3.1, Affymetrix) software offered RMA normalization for both gene and exon level analysis. Genesrping software (version 13.1; Agilent Technologies) was employed to finish the basic analysis. Differentially expressed genes were then identified through fold change as well as P value calculated with t-test. P values were also adjusted for multiple comparison. The threshold set for up- and down-regulated genes was fold change ≥ 2.0 and P value ≤ 0.05. Hierarchical Clustering was performed to display the distinguishable genes' expression pattern among samples.

Results

Expression profiles of lncRNAs and mRNAs in PSTT

In order to identify the expression of lncRNAs in PSTT compared to normal villi, microarray was conducted as described. Principal component analysis (PCA) was performed to explore the expression patterns associated with differentially expressed lncRNAs and mRNAs in our datasets. Both mRNA and lncRNA expression levels could be distinguished between the PSTT and control samples (Figure S1A and S1B). In brief, 1247 up-regulated lncRNAs and 1013 down-regulated lncRNAs were detected (FC ≥ 2.0 and P ≤ 0.05). The top 20 up-regulated and down-regulated lncRNAs are shown in Table 2. Among our data, the most significantly up-regulated lncRNA is NONHSAT114519 with FC of 801.27435, whereas NONHSAT102523 is the most significantly down-regulated lncRNA. Clustering analysis of top 20 up- and down-regulated lncRNAs and volcano plot of lncRNAs expression were shown in Figure 1D, F and B.

 Table 1 

The sequence of the primers for lncRNAs and mRNAs

Gene nameForward and reverse primerTm (℃)Product length (bp)
DPP4F: 5'CTAGGGCAGGGACAGGATAA3'60126
R: 5'TGTGAACAGCTCTTCTCCG3'
PARP14F: 5'ACATTGTGTGCCAGGTAG3'60116
R: 5'GCTTCTTGCACTCTGAGC3'
ADAMTS9F: 5'AATGCTTTGAGTCTTTCCGA3'60112
R: 5'GCTTCCCTTCATCAGCTTG3'
F: 5'ATGAATGCTGTCTGTGTGGAA3'60106
R:5'CAGCAAAGAGTTGCCATATAGT3'
AQPEP//LVRNF: 5'CCATCAGCACATCTCCATTC3'60122
R: 5'TTTACTCACAGCTTGCCAG3'
EGR1F: 5'CAAACCAATGGTGATCCTCTAT3'60101
R: 5'CTGACACATGCTCTGAGAAT3'
PLK2F: 5'ATTAGTCAAGTGACGGTGC3'60115
R: 5'GAAGGAGGTAGAGCCGAG3'
ADAMTS6F: 5'GATCTAATGCAGATGACTAGGC3'60107
R: 5'ATTCCATGCTGATTGTCCAC3'
HTRA4F: 5'CCAATGCCCATGTTGTCAG3'60110
R: 5'CACCGCAAGATCCAATTTAAG3'
CDH11F: 5'CCTGGGTCATTGTGACATA3'60101
R: 5'CCTCTTCTGCTCAGAGACT3'
GAPDHF: 5'TGTTGCCATCAATGACCCCTT3'60201
R: 5'CTCCACGACGTACTCAGCG3'
NONHSAT003875F: 5'GGTGGGGCAGAGAACATAGAAAAAGA3'60246
R: 5'TCAAGGAAGAGTTGGGAAGGAAGAGA3'
NONHSAT102500F: 5'GGGAGCCTTTTCGTTTTGTGCTTTTT3'58202
R: 5'ATTTCGTGCCCTTTGCCTCACTTTTC3'
NONHSAT114519F: 5'AAAGAGGTAGGAGCAAGAAAGAGGAG3'58214
R: 5'TCTATGTGCATATTTGGGATGAGATT3'
NONHSAT134431F: 5'CATTATATAGATGGAAACATCGAGGG3'60264
R: 5'TGCTTATTAGAATTTTTGTGGTGAGA3'
NR_103711F: 5'ACTGGACTGTGCAGTGTGGTTCTGAG3'58190
R: 5'TGGGCATTTTGTTTATTTGTTTGGTG3'
NONHSAT136587F: 5'AAGAGATTTTGTCTAAAGGGCAGCAT3'60326
R: 5'CAGAGAAAGCCAGTCGGTAAGTTCTG3'

Tm: temperature. bp: base pair

 Table 2 

Top 20 significantly up-regulated and down-regulated lncRNAs

Probe Set IDpFC (abs)RegulationPSTT1PSTT2PSTT3PSTT4Villi1Villi2Villi3Villi4GeneSymbolNonCodeID
TC0600002926.oe.12.876E-05801.27435up14.5726617.6193514.2119115.350685.1619945.3357385.8208376.851423lnc-LAMA4-4NONHSAT114519
TC0500001030.oe.12.075E-05734.45276up14.6295116.5072316.4653913.99975.0546275.0781946.5402466.846658lnc-COMMD10-8NONHSAT103347
TC0300003446.oe.10.0001658650.9045up18.8609518.1665818.0566718.815327.4988717.04242410.3655211.6075lnc-TNK2-9NONHSAT094251
TC0300003445.oe.18.53E-05460.6303up16.5580415.3681815.3998216.552255.7956435.4170858.3819658.893738lnc-TNK2-10NONHSAT094250
TC0500001029.oe.16.498E-05330.9917up12.9102615.3704215.5642712.906016.140064.6754546.6972595.755586lnc-COMMD10-9NONHSAT103345
TC0400000688.oe.15.345E-05202.31433up14.1391216.1197114.5787314.293856.0611056.2818387.4723018.674349lnc-COPS4-1NONHSAT097204
TC0500002484.oe.10.0034374105.16381up10.3977311.6398614.56637.6462064.4061094.0587874.5082394.410984RP11-78C3.1NONHSAT102500
TC0900000988.oe.10.0005695100.53162up14.1092116.456717.1416516.37128.7286367.74515311.295179.703777lnc-TRIM32-11NONHSAT134431
TC0900002221.oe.10.000536899.72805up12.9607414.5268314.3840114.131856.125725.5338239.4888538.295328PAPPA-AS1NR_103711
TC0400002102.oe.10.000486889.59299up11.6012715.6326614.0649613.148616.6168276.7947096.6601078.434599lnc-LIN54-3NONHSAT097205
TC0500002714.oe.10.00194886.06065up10.9378913.9987712.59518.3666444.8255024.8579475.6684184.837413lnc-ATG12-2NONHSAT103341
TC1300001424.oe.10.005582771.04717up11.4767310.9982610.8108816.382935.0189815.456837.2696857.320484RP11-318G21.4NONHSAT034434
TC0300003444.oe.19.85E-0556.742313up12.165710.8038411.832312.611625.3557865.1374346.2253557.389472lnc-TNK2-11---
TC1400000340.oe.10.00257654.023994up7.7303511.6516712.146598.1608993.8311984.0729214.6639834.099294lnc-LRR1-6NONHSAT036699
TC0400000156.oe.10.00082653.37624up9.91871512.8425212.9464213.540165.9806945.6785417.0296277.606449lnc-CPZ-1---
TC1500001204.oe.10.00071452.510838up12.266398.41967612.0019211.291064.8423885.1583455.5423425.577796RP11-73C9.1NONHSAT040839
TC0100001151.oe.10.012575751.62064up12.9244913.5507913.0301211.931196.58435911.716534.4203715.955827lnc-IFI44-5NONHSAT004079
TSUnmapped00002029.oe.10.004154951.021782up14.787569.15075813.0463111.741176.1858867.7754356.2551815.817128lnc-RTDR1-1---
TC1200003068.oe.10.000113348.14622up9.0920612.0200910.9762211.147915.5613614.9524835.2112535.153783lnc-HCAR3-1NONHSAT031531
TC1600001762.oe.10.00692743.289265up14.5136112.5437914.50811.075358.79909310.19186.0026275.903477lnc-KIAA0430-2NONHSAT140763
TC0500002493.oe.10.0370893288.59473down5.6389685.0177364.9477585.2094217.5998888.59555818.26619.04404lnc-EDIL3-6NONHSAT102523
TC1100001750.oe.10.001197695.452down5.7167516.6663996.1268715.76496310.4544410.9711214.4962214.66002NONHSAG007410NONHSAT017479
TC0200004232.oe.14.448E-0575.3497down8.2180219.4050649.3793679.75737216.3894316.1581514.5875214.56685lnc-GCG-3NONHSAT075191
TC0500002278.oe.15.254E-0573.73616down7.368565.656326.7686176.70167111.4604212.7385713.5009713.61241lnc-CENPK-2NONHSAT101743
TC0200004231.oe.16.523E-0566.59564down8.5895239.653439.1790549.52641116.2733216.2594414.0654414.57964lnc-GCG-4NONHSAT075190
TC0Y00000312.oe.13.009E-0548.222668down7.4167267.0513437.6810147.15298813.7042113.7582711.8866212.31953lnc-HSFY2-10NONHSAT139611
TC0500000786.oe.10.000359845.157536down8.9373918.6563069.23492211.0287716.1497815.5981113.8834614.21362lnc-XRCC4-6NONHSAT102522
TC1200000179.oe.10.01356935.71296down8.3670528.45798410.008089.47909417.0908516.3340911.9049611.61581lnc-RIMKLB-1NONHSAT026363
TC0900001993.oe.10.042837534.091114down15.4773510.672613.606766.98858116.3627414.9202218.0839117.74372lnc-OMD-1NONHSAT133188
TC14000659.hg.40.005272730.134325down6.7586596.7551268.2070186.67537610.2187710.0337813.8184413.97853MEG8NR_024149
TC1200000178.oe.10.020378325.076464down6.3393857.0469548.7300948.73053114.9169214.4946410.0195210.00893lnc-RIMKLB-3NONHSAT026360
TC1400000032.oe.14.427E-0523.937801down5.9084656.8093186.1802696.86119110.7788510.0851211.6023711.61778lnc-APEX1-3NONHSAT035586
TC0800000155.oe.10.021928521.791851down7.0522046.4858747.2141297.36055313.4618214.357958.4996999.576158lnc-PDGFRL-3NONHSAT125277
TSUnmapped00000268.oe.10.001886321.46142down6.5924076.483017.715037.0815410.669549.54328212.3083913.04547lnc-DLK1-5---
TC1400000966.oe.10.004921620.716341down6.7278866.5157256.8676926.9531299.0942639.74558712.6448213.07055OTTHUMG00000029060NONHSAT039835
TC0500001872.oe.10.004010219.861301down9.4338716.5349477.7194276.94948213.2309813.1002910.404711.14931lnc-TAS2R1-24NONHSAT100322
TC0300002485.oe.10.001471519.597805down7.213317.917227.5397167.54604613.6319712.5888210.7367610.42922lnc-PRICKLE2-6NONHSAT090260
TC0800000049.oe.11.843E-0518.799055down7.1307745.7139986.7990216.78656411.2065610.6759410.4929310.98528lnc-MCPH1-3NONHSAT124787
TC0600001704.oe.10.000224517.906603down6.2368464.6434256.3695435.588889.4457369.7242210.911749.406675lnc-SLC22A1-6NONHSAT115901
TSUnmapped00000213.oe.10.00630917.698973down6.0898417.5336547.4966036.49641413.1971112.048529.2609249.692332lnc-COLEC10-1---

P values were calculated using unpaired t-test and adjusted for multiple comparison. FC (abs): the absolute ratio (no log scale) of average normalized intensities between PSST and normal villi. PSTT 1 to 4 and Villi 1 to 4: each sample's normalized intensity (log2 scale). The same goes for the follow-up table.

We also demonstrated 828 up-regulated mRNAs and 1393 down-regulated mRNAs in PSTT tissues compared to normal villi (fold change ≥ 2.0, p ≤ 0.05). The top 20 up-regulated and down-regulated mRNAs are shown in Table 3. Among those upregulated mRNAs, AQPEP also known as LVRN ranges as top with an FC of 591.0085, while expression of CXCL14 tops among downregulated genes with FC of 1304.5027. The heatmap plot showed the clustering analysis among mRNA expression profiles by revealing the top 20 up- and downregulated mRNAs (Figure 1 C, E). A volcano plot (Figure 1A) showed the variation of mRNA expression between PSTT tissues and normal villi. The heatmaps and hierarchical clustering were drawn by the TBtools [23].

Validation of the aberrantly expressed lncRNAs and mRNAs by RT-qPCR

To confirm the reliability of the microarray profiling results, six lncRNAs (NONHSAT114519, NR_103711, NONHSAT003875, NONHSAT136587, NONHSAT134431, NONHSAT102500) as well as ten mRNAs (DPP4, PARP14, ADAMTS9, PLAC8, AQPEP, EGR1, PLK2, ADAMTS6, HTRA4, CDH11) were chosen randomly for RT-qPCR. The microarray samples including four PSTT tissues and normal villi counterparts were further used for RT-qPCR. The results of RT-qPCR were in accordance with microarray observations. In PSTT tissues, the lncRNA expression of NONHSAT114519, NR_103711, NONHSAT003875, NONHSAT136587, NONHSAT134431, NONHSAT102500 was up-regulated when compared with those in normal placental villi (Figure 2C, D). Further, the mRNA expression of PARP14, PLAC8, AQPEP, EGR1, PLK2, and HTRA4 was significantly up-regulated whereas the expression of DPP4, ADAMTS9, ADAMTS6 and CDH11 in PSTT tissues was significantly down-regulated compared with those in normal placental villi (Figure 2A, B).

Validation of the PLAC8, EGR1 and ADAMTS6 by immunohistochemistry

To further confirm the microarray data, the protein products of several genes were measured by immunohistochemistry analysis. We randomly chose two up-regulated mRNAs (PLAC8 and EGR1) as well as one down-regulated mRNAs (ADAMTS6) for immunohistochemistry. As shown in Figure 2F, PLAC8 and EGR1 expression was dramatically elevated in PSTT tissues compared with those in normal villi. Comparatively, ADAMTS6 was diffusely positive in normal villi but only partly positive in PSTT. Immunohistochemistry data were consistent with microarray analysis. Furthermore, we also found that PLAC8 was highly expressed in PSTT but not choriocarcinoma, which indicates that PLAC8 may act as a potential marker to distinguish PSTT from choriocarcinoma Figure S2.

 Table 3 

Top 20 significantly up-regulated and down-regulated mRNAs

Probe Set IDpFC (abs)RegulationPSTT1PSTT2PSTT3PSTT4Villi1Villi2Villi3Villi4GeneSymbolstrand
TC05000549.hg.43.19235E-05591.0085up14.2750616.4270115.9755513.535964.7186535.5112046.8632366.292351AQPEP+
TC08000297.hg.40.004722227458.7278up17.7326218.9531317.8910517.817025.7015776.16989411.1787413.97764HTRA4+
TC04001304.hg.40.001240178298.35522up17.3165711.5782613.0454616.068375.2851857.7904335.5999196.449578CXCL9-
TC05001550.hg.40.001369429220.83633up10.2493111.6082514.992168.4929873.1933233.2666553.9474323.787961RP11-78C3.1-
TC02004950.hg.40.000845117215.45534up16.6813712.5343415.6168214.568167.1220769.6269415.7311625.915533IGKV3-15-
TC04001753.hg.40.000546863160.13506up13.5139217.1688415.1103715.331657.0050786.2161959.0349139.57601HPGD-
TC21000287.hg.40.000300168154.46967up10.5937213.8992714.5251711.974994.8895614.7693786.1801886.069306AF127936.7-
TC06000935.hg.40.002012451119.142746up15.0382512.5938314.7477911.787787.7603298.9652665.6003944.255473PLN+
TC04001344.hg.40.000890978111.79973up12.0553515.7325814.565313.716136.1785755.5593758.1196978.99262PLAC8-
TC01002774.hg.40.001239457107.338455up14.5702312.0693714.9040911.532057.0590748.579035.2663125.187234PTGER3-
TC09001525.hg.40.000788934105.381905up12.4972613.8394813.6767613.498055.1655124.8039379.221767.44241PAPPA-AS1-
TC04001270.hg.40.00418030893.43535up14.726287.91144311.891210.378144.4476765.5327824.0420934.700925IGJ-
TC06001773.hg.40.00068225490.94472up17.5207117.2110717.0734317.8725512.4259511.8365311.409557.978059PLA2G7-
TC09000038.hg.41.45655E-0589.93031up12.5437613.8334811.7746213.45326.0690046.9978766.1030816.472157PDCD1LG2+
TC08000150.hg.40.00150732877.45381up9.45816413.665419.14114912.476324.5753694.4673565.2685245.328737AC100802.3+
TC16000898.hg.40.00488627972.98015up15.871713.4805815.8759912.171579.26975310.689776.4778796.204714MYH11-
TC02004978.hg.40.00753926668.028145up12.409445.7040459.0748511.440053.0221044.6920993.4583963.103549IGKV2D-28+
TC09000593.hg.40.00325279766.812675up13.5935416.2480415.8784115.894787.7472417.05989711.4590411.10039PAPPA+
TC04001305.hg.40.00330423756.946953up13.4206110.2169.27474913.002675.7174867.5890984.44574.835558CXCL10-
TC07001008.hg.40.00671853356.503063up8.28817711.4720510.5050814.019445.6175467.3704143.7476924.268066AOC1+
TC05001804.hg.41.13295E-051304.5027down4.7896944.7813644.2353594.6192616.1616116.3994513.4728113.78895CXCL14-
TC09000799.hg.40.018049004212.94751down4.279575.4846414.44045.48187416.7512416.768828.0939139.009928PAEP+
TC02002478.hg.46.23932E-0597.702down9.55595211.2769310.5813811.2376318.2853518.1081216.4398716.25981DPP4-
TC07000722.hg.48.31021E-0584.22738down5.3715434.9376155.3336885.77031613.1277312.6002910.3906710.87934MET+
TC05001406.hg.41.12984E-0581.3294down6.726225.4532246.3828955.67525111.4857112.1207812.8900113.12391ADAMTS6-
TC04001335.hg.42.90216E-0567.747826down6.1261984.3394624.1209815.20458511.4027111.5587710.8573410.30082TMEM150C-
TC12000624.hg.40.00499546561.509308down4.4026995.3014053.5621414.0609877.5291988.50219412.187512.87927LGR5+
TC17001112.hg.44.01987E-0650.456596down5.3169125.6620825.718195.52193111.6467911.9393910.4881510.77267SLC25A35-
TC11001274.hg.40.00217468649.247604down6.9862329.4915817.0713698.25341411.6362212.3872214.771715.49538INS-
TC14000407.hg.40.00122055848.467384down5.2330554.433025.9840236.45733810.242699.10938512.4546812.69645PLEKHH1+
TC11001505.hg.40.00362468547.256374down5.2703978.9293166.7620515.76006912.4211110.8254311.0314814.69356MUC15-
TC17001474.hg.40.0179039947.178997down6.294736.8186488.2455326.97100511.134938.63837615.3930915.40381KRT23-
TC09001334.hg.40.04252794846.31269down15.7860311.1422414.074096.86964116.9930615.1875219.4082118.41655OGN-
TC05001151.hg.40.00090672845.62948down9.5843216.1915128.450698.7039814.7029414.6260312.4779313.17118SEMA5A-
TC13000335.hg.40.03522612545.062977down11.858077.95248211.398163.79209415.2373615.6691712.1462313.92353CLDN10+
TC16000223.hg.40.00738961343.19478down4.7414534.3902894.6575234.42030712.7883611.864167.8206787.467516ACSM3+
TC08000022.hg.47.42939E-0642.766193down8.3762947.9542879.1394689.32145614.2242213.5247514.4393914.27674AGPAT5+
TC02000293.hg.40.04033672839.249073down7.8927645.3278167.7822516.06566415.5164315.246318.2151179.268984EPCAM+
TC01002886.hg.40.00053754437.86066down6.8179555.4227426.1180526.6531559.51977411.3545312.5775812.53053F3-
TC14000056.hg.40.00010286337.460266down6.5685048.0841416.3706897.93447212.0322711.6753112.9264413.23294PNP+
 Figure 1 

Differential mRNA and lncRNA expression of PSTT and villi by microarray. The volcano plot represents total identified mRNAs (A) and lncRNAs (B) expression between the PSTT and normal villi. The horizontal black lines show the default 2.0-fold change and the vertical black lines represent a P-value of 0.05. The red and blue plots represent up- and downregulated RNAs with FC ≥ 2.0 and P < 0.05. Heat map and hierarchical clustering analysis of the top 20 up- and downregulated mRNAs (C, E) and lncRNAs (D, F) between PSTT and normal villi.

Int J Med Sci Image

(View in new window)

 Figure 2 

Validation of the microarray results of mRNAs and lncRNAs by RT-qPCR and immunohistochemistry. RT-qPCR was performed to test the differentially expressed mRNAs (A) and lncRNAs (C) (n=4 for PSTT group and control group, respectively); the fold change of each mRNA (B) and lncRNA (D) between PSTT tissues and control group was determined microarray and RT-qPCR. The expression of PLAC8 and EGR1 was obviously higher than normal villi, while ADAMTS6 expression is much lower in PSTT compared to normal villi (F). ns, No significance; *, P < 0.05; **, P < 0.01, P < 0.001, Student t-test. Scale bar = 100μm.

Int J Med Sci Image

(View in new window)

GO and KEGG analyses of differential genes

GO analysis of the altered 2221 mRNAs showed that up-regulated genes were enriched in immune response (ontology: biological process, GO:0006955), actin cytoskeleton (ontology: cellular component, GO:0015629) and actin binding (ontology: molecular function, GO:0003779) (Figure 3B); On the other hand, the top 10 enriched GO terms on down-regulated genes included small molecule metabolic process (ontology: biological process, GO:0044281), mitochondrion (ontology: cellular component, GO:0005739 and structural constituent of ribosome (ontology: molecular function, GO:0003735) (Figure 3A).

 Figure 3 

GO and KEGG analysis of altered mRNAs. The GO analysis of down-regulated mRNAs (A) and up-regulated mRNAs (B). KEGG pathways showed TOP 10 significantly enriched pathways of down-regulated mRNAs(C) and up-regulated mRNAs(D).

Int J Med Sci Image

(View in new window)

KEGG analysis was of these differentially expressed mRNAs was also conducted. Up-regulated mRNAs were involved in Measles, Herpes simplex infection, Chemokine signaling pathway (Figure 3D), while down-regulated mRNAs were most enriched in Ribosome, Metabolic pathways, Biosynthesis of antibiotics (Figure 3C).

CNC network with GO and KEGG analyses

We then constructed lncRNA-mRNA co-expression networks on the basis of 6 validated lncRNAs as well as 354 mRNAs. In brief, the CNC networks were composed of 497 edges and 360 nodes with 207 positive interactions (continuous lines) and 290 negative interactions (dotted lines) (Figure 4). The lncRNA NONHSAT114519 is connected with 151 mRNAs, NR_103711 is connected with 97 mRNAs, NONHSAT003875 is connected with 77 mRNAs, NONHSAT136587, NONHSAT134431, and NONHSAT102500 are connected with 71, 66 and 35 mRNAs respectively.

In order to explore the functions of the lncRNAs, GO and KEGG analysis were also performed based on the interacted mRNAs from CNC networks. SRP-dependent cotranslational protein targeting to membrane (ontology: biological process, GO:0006614), focal adhesion (ontology: cellular component, GO:000592) and structural constituent of ribosome (ontology: molecular function, GO:0003735) were the most significantly enriched GO terms of the interacted genes (Figure 5A). KEGG data indicated that those targeted mRNAs were enriched in Ribosome, Proteoglycans in cancer, cGMP-PKG signaling pathway (Figure 5B).

Discussion

A number of studies have demonstrated that lncRNAs play important roles in the differentiation [24, 25] and biological behaviors [26-32] of trophoblast cells. The impaired trophoblast migration and invasion abilities may further lead to preeclampsia and miscarriage. PSTT is originated from abnormally differentiated intermediate trophoblasts [2]. However, aberrant lncRNAs expression profiles remain unknown in GTN, especially in PSTT. Presumably, the aberrantly expressed lncRNAs may contribute to the development and progression of PSTT. In the present study, expression levels of lncRNAs and mRNAs in four human PSTT tissues and four normal villi tissues were investigated by integrated microarray analysis. 2221 significantly altered mRNAs were identified (including 828 up-regulated; 1393 down-regulated) and 2260 differentially expressed lncRNAs (1247 upregulated; 1013 downregulated) in the PSTT compared to normal villi.

 Figure 4 

CNC networks by validated lncRNAs. The CNC networks were performed by 6 validated lncRNAs and 354 interacted mRNAs. The CNC networks were composed of 497 edges and 360 nodes with 207 positive interactions (continuous lines) and 290 negative interactions (dotted lines). The green nodes represented down-regulated mRNAs and blue ones denoted up-regulated mRNAs.

Int J Med Sci Image

(View in new window)

 Figure 5 

GO and KEGG analyses of interacted mRNAs in CNC networks. GO (A) and KEGG (B) analyses of CNC interacted mRNAs

Int J Med Sci Image

(View in new window)

Our GO analysis showed that the most up-regulated mRNAs was involved in immune system process, while the most down-regulated mRNAs participate in mitochondrion.

Furthermore, the KEGG pathway analysis demonstrated that the aberrantly expressed mRNAs are mainly enriched in ribosome, metabolic pathways and biosynthesis of antibiotics (down-regulated mRNAs); measles, herpes simplex infection and chemokine signaling pathway (up-regulated mRNAs). Similarly, the most enriched GO terms are SRP-dependent cotranslational protein targeting to membrane, focal adhesion and structural constituent of ribosome of the lncRNA-interacted mRNAs. Furthermore, the top one enriched KEGG pathway is also ribosome. Such analysis inferred that PSTT might be mainly immune-regulated by ribosome. Besides, we also found that HMOX1, EPHB2, CLIC4 and CCL2 were up-regulated in PSTT, which were identified to be angiogenesis-related genes in previous studies [33-36]. Protein-protein interaction network based on ten validated mRNAs were shown in Figure S3.

The gene expressions of our research were correlated with previous reported biomarkers including GATA3, which was known expressed in 71% of PSTT [37]. Our results showed that compared to normal villi MAPK14, MAPK13, MAPK7 and MAPK1 were significantly down-regulated in PSTT, but Köbel et al. have reported that MAPK was highly expressed [38]. PLAC8, also named onzin, is a small protein (∼16 kDa) that was originally described to be highly expressed in mouse placenta [39] and its deficiency in mice may result in innate immunity deficiency [40]. We also found that PLAC8 is highly expressed in PSTT but not choriocarcinoma (Figure S2). Our previous study has indicated that PLAC8 may promote autophagic activity and improves the growth priority of human trophoblast cells [41]. We have also demonstrated that the overexpression of AQPEP//LVRN could influence the invasion of trophoblast [42]. In addition, NONHSAT003875/miR-363/EGR1 regulatory network in the carcinoma -associated fibroblasts was confirmed to control the angiogenesis of PSTT [43]. Further studies focusing on these proteins may be promising methods to a better understanding of the pathogenesis process in PSTT.

To further study the changed signaling pathways between different types of GTN. GSEA analysis were performed on the basis of data set GSE135727. Significant differences in signaling pathways were identified by nominal p < 0.05 and FDR q-value < 0.25 (Table S1). Four representative GSEA-enrichment plots were indicated in Figure S4. Compared to PSTT, choriocarcinoma tended to alter in Coagulation, KRAS, P53 and MYC targets pathways. KRAS and P53 signaling pathways were significantly enriched in ETT and CC. GSEA analysis also suggest that alteration in the MYC targets pathway might be one of the major pathways altered in PSTT when compared to normal villi.

In conclusion, we revealed that lncRNAs and mRNAs are differentially expressed in PSTT and normal villi by microarray analysis. GO and KEGG showed that immunotherapy may be effective in PSTT. Furthermore, GO and KEGG analyses of lncRNA-interacted mRNAs based on 6 validated lncRNAs were also performed, which indicated that these specific lncRNAs may be involved in the biological processes of ribosome that might contribute to PSTT pathogenesis. GSEA analysis based on GSE135727 and microarray data set among different types GTN were performed to study the altered signaling pathways. The function of these lncRNAs in PSTT required further research.

Supplementary Material

Supplementary figures and table.

Attachment

Acknowledgements

This work was funded by the National Natural Science Foundation of China [grant number: 81971394, 81571457], the Shanghai Pujiang Program [grant number: 15PJ1400900].

Author Contributions

Conceptualization, Jianfeng Gan and Hongbo Zhao; Methodology, Sai Zhang, Xuan Feng and Zhi Wei; Investigation, Jianfeng Gan, Zhixian Chen and Yan Du; Writing - Original Draft, Jianfeng Gan; Writing -Review & Editing, Jianfeng Gan and Hongbo Zhao; Funding Acquisition, Hongbo Zhao and Congjian Xu; Resources, Hongbo Zhao and Congjian Xu.; Supervision, Hongbo Zhao and Congjian Xu; Literature search and critical revision, Jianfeng Gan; Conception of the work and final version approval, Jianfeng Gan and Hongbo Zhao. All authors read and approved the final manuscript.

Abbreviations

PSTT: placental site trophoblastic tumor; RT-qPCR: reverse transcription and quantitative real-time polymerase reaction; GO: Gene Ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes; CNC: coding-non-coding gene co-expression; PCCs: Pearson correlation coefficients; SRP: signal-recognition particle; GTN: gestational trophoblastic neoplasms; lncRNAs: long non-coding RNAs; GSEA: gene set enrichment analysis; ETT: epithelioid trophoblastic tumor; CC: choriocarcinoma; PCA: principal component analysis.

Competing Interests

The authors have declared that no competing interest exists.

References

1. Kurman RJ, Scully RE, Norris HJ. Trophoblastic pseudotumor of the uterus: an exaggerated form of "syncytial endometritis" simulating a malignant tumor. Cancer. 1976;38:1214-26

2. Feng X, Wei Z, Zhang S, Du Y, Zhao H. A Review on the Pathogenesis and Clinical Management of Placental Site Trophoblastic. Frontiers in oncology. 2019;9:937

3. Feng X, Wei Z, Zhang S, Du Y, Zhao H. A Review on the Pathogenesis and Clinical Management of Placental Site Trophoblastic. Frontiers in oncology. 2019;9:937

4. Seckl MJ, Sebire NJ, Fisher RA, Golfier F, Massuger L, Sessa C. Gestational trophoblastic disease: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up. Annals of oncology: official journal of the European Society for Medical Oncology. 2013;24(Suppl 6):vi39-50

5. Schmid P, Nagai Y, Agarwal R, Hancock B, Savage PM, Sebire NJ, Lindsay I, Wells M, Fisher RA, Short D, Newlands ES, Wischnewsky MB, Seckl MJ. Prognostic markers and long-term outcome of placental-site trophoblastic tumours: a retrospective observational study. Lancet (London, England). 2009;374:48-55

6. Lan C, Li Y, He J, Liu J. Placental site trophoblastic tumor: lymphatic spread and possible target markers. In. 2010 p. 430-7

7. Papadopoulos AJ, Foskett M, Seckl MJ, McNeish I, Paradinas FJ, Rees H, Newlands ES. Twenty-five years' clinical experience with placental site trophoblastic tumors. The Journal of reproductive medicine. 2002;47:460-4

8. Hassadia A, Gillespie A, Tidy J, Everard R G N J, Wells M, Coleman R, Hancock B. Placental site trophoblastic tumour: clinical features and management. Gynecologic oncology. 2005;99:603-7

9. Zhao J, Lv WG, Feng FZ, Wan XR, Liu JH, Yi XF, Qu PP, Xue FX, Wu YM, Zhao X, Ren T, Yang JJ, Xie X, Xiang Y. Placental site trophoblastic tumor: A review of 108 cases and their implications for prognosis and treatment. Gynecologic oncology. 2016;142:102-108

10. Lee H, Shin W, Jang YJ, Choi CH, Lee J, Bae D, Kim B. Clinical characteristics and outcomes of placental site trophoblastic tumor: experience of single institution in Korea. Obstetrics & gynecology science. 2018;61:319-327

11. Braga A, Mora P, de Melo AC, Nogueira-Rodrigues A, Amim-Junior J, Rezende-Filho J, Seckl MJ. Challenges in the diagnosis and treatment of gestational trophoblastic neoplasia worldwide. In. 2019 p. 28-37

12. Wang C, Zhu C, Xu J, Wang M, Zhao W, Liu Q, Zhao G, Zhang Z. The lncRNA UCA1 promotes proliferation, migration, immune escape and inhibits apoptosis in gastric cancer by sponging anti-tumor miRNAs. Molecular Cancer. 2019 18

13. Yang L, Lin C, Jin C, Yang JC, Tanasa B, Li W, Merkurjev D, Ohgi KA, Meng D, Zhang J, Evans CP, Rosenfeld MG. lncRNA-dependent mechanisms of androgen-receptor-regulated gene activation programs. Nature. 2013;500:598-602

14. Sun X, Yuan Y, Xiao Y, Lu Q, Yang L, Chen C, Guo Q. Long non-coding RNA, Bmcob, regulates osteoblastic differentiation of bone marrow mesenchymal stem cells. Biochemical and biophysical research communications. 2018;506:536-542

15. Wambecke A, Ahmad M, Morice P, Lambert B, Weiswald L, Vernon M, Vigneron N, Abeilard E, Brotin E, Figeac M, Gauduchon P, Poulain L, Denoyelle C, Meryet-Figuiere M. The lncRNA 'UCA1' modulates the response to chemotherapy of ovarian cancer through. Molecular oncology. 2021

16. Thakur KK, Kumar A, Banik K, Verma E, Khatoon E, Harsha C, Sethi G, Gupta SC, Kunnumakkara AB. Long noncoding RNAs in triple-negative breast cancer: A new frontier in the. Journal of cellular physiology. 2021

17. Xu W, Ding M, Wang B, Cai Y, Guo C, Yuan C. Molecular mechanism of the canonical oncogenic lncRNA MALAT1 in gastric cancer. Current medicinal chemistry. 2021

18. Wang X, Li X, Zhou Y, Huang X, Jiang X. Long non-coding RNA OIP5-AS1 inhibition upregulates microRNA-129-5p to repress. Cell biology and toxicology. 2021

19. Zhang M, Gao F, Yu X, Zhang Q, Sun Z, He Y, Guo W. LINC00261: a burgeoning long noncoding RNA related to cancer. Cancer Cell International. 2021 21

20. Wang ZK, Yang L, Wu LL, Mao H, Zhou YH, Zhang PF, Dai GH. Long non-coding RNA LINC00261 sensitizes human colon cancer cells to cisplatin therapy. Brazilian journal of medical and biological research = Revista brasileira de pesquisas medicas e biologicas. 2017;51:e6793

21. Cho EJ, Chun S, Park H, Sung CO, Kim K. Whole transcriptome analysis of gestational trophoblastic neoplasms reveals altered PI3K signaling pathway in epithelioid trophoblastic tumor. Gynecologic Oncology. 2020;157:151-160

22. Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic acids research. 2002;30:207-10

23. Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, Xia R. TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Molecular plant. 2020;13:1194-1202

24. Chu T, Mouillet J, Cao Z, Barak O, Ouyang Y, Sadovsky Y. RNA Network Interactions During Differentiation of Human Trophoblasts. Frontiers in Cell and Developmental Biology. 2021 9

25. You J, Wang W, Chang H, Yi Y, Zhao H, Zhu H, Sun Y, Tang M, Wang C, Sang Y, Feng G, Cheng S, Leung PCK, Zhu Y. The BMP2 Signaling Axis Promotes Invasive Differentiation of Human Trophoblasts. Frontiers in cell and developmental biology. 2021;9:607332

26. Wang R, Zou L, Yang X. microRNA-210/ Long non-coding RNA MEG3 axis inhibits trophoblast cell migration and invasion by suppressing EMT process. Placenta. 2021;109:64-71

27. Chen X, Guo D, Yin T, Yang J. Non-Coding RNAs Regulate Placental Trophoblast Function and Participate in Recurrent Abortion. Frontiers in Pharmacology. 2021 12

28. Chi Z, Gao Q, Sun Y, Zhou F, Wang H, Shu X, Zhang M. LINC00473 downregulation facilitates trophoblast cell migration and invasion via the. Environmental toxicology. 2021

29. Fan X, Lou J, Zheng X, Wang Y, Wang J, Luo M, Hu M. Interference with lncRNA NEAT1 promotes the proliferation, migration, and invasion. Immunopharmacology and immunotoxicology. 2021;43:334-342

30. Xie J, Liang T, Zhao J, Xu Z, Tian P, Wang R, Mi C, Huang W, Chen W, Zhang H. Lnc-HZ08 regulates BPDE-induced trophoblast cell dysfunctions by promoting PI3K. Cell biology and toxicology. 2021

31. Xu Z, Tian P, Guo J, Mi C, Liang T, Xie J, Huang W, Dai M, Chen W, Zhang H. Lnc-HZ01 with m6A RNA methylation inhibits human trophoblast cell proliferation and. The Science of the total environment. 2021;776:145950

32. Chen Y, Li Z, Chen X, Zhang S. Long non-coding RNAs: From disease code to drug role. Acta pharmaceutica Sinica. B. 2021;11:340-354

33. Luu Hoang KN, Anstee JE, Arnold JN. The Diverse Roles of Heme Oxygenase-1 in Tumor Progression. Frontiers in immunology. 2021;12:658315

34. Sato S, Vasaikar S, Eskaros A, Kim Y, Lewis JS, Zhang B, Zijlstra A, Weaver AM. EPHB2 carried on small extracellular vesicles induces tumor angiogenesis via. JCI insight. 2019 4

35. Edwards JC, Bruno J, Key P, Cheng Y. Absence of chloride intracellular channel 4 (CLIC4) predisposes to acute kidney. BMC nephrology. 2014;15:54

36. Bonapace L, Coissieux M, Wyckoff J, Mertz KD, Varga Z, Junt T, Bentires-Alj M. Cessation of CCL2 inhibition accelerates breast cancer metastasis by promoting angiogenesis. Nature. 2014;515:130-133

37. Banet N, Gown AM, Shih I, Kay Li Q, Roden RBS, Nucci MR, Cheng L, Przybycin CG, Nasseri-Nik N, Wu L, Netto GJ, Ronnett BM, Vang R. GATA-3 expression in trophoblastic tissues: an immunohistochemical study of 445 cases, including diagnostic utility. In. 2015 p. 101-8

38. Köbel M, Pohl G, Schmitt WD, Hauptmann S, Wang T, Shih I. Activation of Mitogen-Activated Protein Kinase Is Required for Migration and Invasion of Placental Site Trophoblastic Tumor. The American Journal of Pathology. 2005;167:879-885

39. Tatura M, Schmidt H, Haijat M, Stark M, Rinke A, Diels R, Lawlor RT, Scarpa A, Schrader J, Hackert T, Schimmack S, Gress TM, Buchholz M. Placenta-Specific 8 Is Overexpressed and Regulates Cell Proliferation in Low-Grade Human Pancreatic Neuroendocrine Tumors. Neuroendocrinology. 2020;110:23-34

40. Ledford JG, Kovarova M, Koller BH. Impaired host defense in mice lacking ONZIN. Journal of immunology (Baltimore, Md.: 1950). 2007;178:5132-43

41. Feng X, Wei Z, Tao X, Du Y, Wu J, Yu Y, Yu H, Zhao H. PLAC8 promotes the autophagic activity and improves the growth priority of human trophoblast cells. The FASEB Journal. 2021 35

42. Feng X, Wei Z, Zhang S, Zhou J, Wu J, Luan B, Du Y, Zhao H. Overexpression of LVRN impedes the invasion of trophoblasts by inhibiting epithelial-mesenchymal transition. Acta Biochimica et Biophysica Sinica. 2021;53:249-257

43. Zhang S, Tao X, Cao Q, Feng X, Wu J, Yu H, Yu Y, Xu C, Zhao H. lnc003875/miR-363/EGR1 regulatory network in the carcinoma -associated fibroblasts controls the angiogenesis of human placental site trophoblastic tumor (PSTT). Experimental Cell Research. 2020;387:111783

Author contact

Corresponding address Corresponding authors: Hongbo Zhao, 413 Zhaozhou Rd, Shanghai 200011, China. Congjian Xu, 128 Shenyang Rd, Shanghai 200090, China. E-mail addresses: zhaohongbo01com (H. Zhao); xucongjiancom (C. Xu).


Received 2021-8-10
Accepted 2021-10-19
Published 2022-1-1