Skip to main content

RNA-sequencing profiling analysis of pericyte-derived extracellular vesicle–mimetic nanovesicles-regulated genes in primary cultured fibroblasts from normal and Peyronie’s disease penile tunica albuginea



Peyronie’s disease (PD) is a severe fibrotic disease of the tunica albuginea that causes penis curvature and leads to penile pain, deformity, and erectile dysfunction. The role of pericytes in the pathogenesis of fibrosis has recently been determined. Extracellular vesicle (EV)–mimetic nanovesicles (NVs) have attracted attention regarding intercellular communication between cells in the field of fibrosis. However, the global gene expression of pericyte-derived EV–mimetic NVs (PC–NVs) in regulating fibrosis remains unknown. Here, we used RNA-sequencing technology to investigate the potential target genes regulated by PC–NVs in primary fibroblasts derived from human PD plaque.


Human primary fibroblasts derived from normal and PD patients was cultured and treated with cavernosum pericytes isolated extracellular vesicle (EV)–mimetic nanovesicles (NVs). A global gene expression RNA-sequencing assay was performed on normal fibroblasts, PD fibroblasts, and PD fibroblasts treated with PC–NVs. Reverse transcription polymerase chain reaction (RT-PCR) was used for sequencing data validation.


A total of 4135 genes showed significantly differential expression in the normal fibroblasts, PD fibroblasts, and PD fibroblasts treated with PC–NVs. However, only 91 contra-regulated genes were detected among the three libraries. Furthermore, 20 contra-regulated genes were selected and 11 showed consistent changes in the RNA-sequencing assay, which were validated by RT-PCR.


The gene expression profiling results suggested that these validated genes may be good targets for understanding potential mechanisms and conducting molecular studies into PD.

Peer Review reports


Peyronie’s disease (PD) is caused by excessive fibrosis and scar tissue formation in the tunica albuginea (TA), resulting in penile pain, abnormal curvature, and erectile dysfunction (ED) [1, 2]. Although the existence of PD has been known for a long time, the pathophysiology of PD has not been studied as widely as fibrosis in other organs, such as the kidneys, liver, or lungs. Currently, the most available medical therapy is collagenase and interferon injection and surgical intervention [3, 4]. However, these treatments can cause glandular hypoesthesia and a high risk of new onset ED [5]. Therefore, the identification of novel therapeutic targets related to PD fibrosis is required.

Pericytes play a fundamental role in vascular contractility and stability, regulation of vascular development, and as a storage vault of mesenchymal stem cells [6, 7]. In vitro studies have shown that pericytes exhibit fibrogenic potential [8, 9] and transition to myofibroblasts [10]. Moreover, the inhibition of angiogenesis may be effective in the suppression of fibrosis [11]. However, recently studies have shown that the inhibition of angiogenesis may aggravate fibrosis [12, 13]. These findings suggest that different antiangiogenic and molecular targets produce different results in the treatment of fibrosis. We recently reported in a mouse model of diabetic ED that pericyte-derived angiogenic factor restored erectile function by enhancing cavernous angiogenesis [14, 15].

Extracellular vesicles (EVs) were previously believed to be cell excretions. However, a number of studies have shown that EVs contain proteins, lipids, and RNA, which can affect the physiological and pathological communications between cells [16, 17]. Many studies regarding the potential role of EVs have been conducted for human diseases, including strokes [18], tumor metastasis [19], and kidney disease [20]. Therefore, clarifying the role of EVs in fibrosis would be beneficial to aid in the understanding of fibrosis mechanisms. However, one of the major limitations of EVs is the low production yield [21]. Therefore, to maximize the production of vesicles, we used a mini extruder system and extracted more than 100-fold greater EV–mimetic NVs from pericytes. The cell-derived EV–mimetic NVs showed similar characteristics to the natural EVs [22]. We recently demonstrated that intracavernous injection of PC–NVs induced penile angiogenesis and rescued erectile function in a mouse model of cavernous nerve injury [15]. Along with potential relationship between angiogenesis and suppression of tissue fibrosis [11,12,13], the angiogenic potential of pericytes and PC–NVs led us to investigate the potential role of PC–NVs in the fibrogenic process of PD.

Gene expression profiling analysis in physiological and pathological conditions can provide a foundation for studying the mechanisms of fibrosis in PD. In the present study, we performed an RNA-sequencing assay on normal fibroblasts, PD fibroblasts, and PD fibroblasts treated with PC–NVs.


Ethics statement and study design

All TA tissues and animals used in this study were approved by the Institutional Review Board (IRB No: 2007-730) and the Institutional Animal Care and Use Committee of our University (approval number: 171129-527), respectively. The plaque tissue of a patient with PD (48 years old) and the normal TA tissue from control patients (undergoing penoplasty for congenital curvature, 21 years old) were used for the human fibroblast culture study. In addition, 10 adult male C57BL/6J mice (8 weeks old, Orient Bio, Korea) were used for the mouse cavernous pericytes (MCPs) primary culture.

Primary culture and characterization of human fibroblasts

The TA tissues were used for the primary fibroblast culture as described previously [23, 24]. Briefly, PD plaque and normal TA tissues were maintained in sterile vials with Hank’s balanced salt solution (HBSS, Gibco, Carlsbad, CA, USA) and washed three times with phosphate-buffered saline (PBS). The TA tissues were cut into 1–2 mm sections and incubated in 12.5 mL Dulbecco’s modified Eagle’s medium (DMEM, Gibco) supplemented with 0.06% collagenase A (Sigma-Aldrich, St. Louis, MO, USA) at 37 °C for 1 h in a 5% CO2 atmosphere. The cells and tissue fragments were collected by centrifugation (400 g for 5 min), washed with fresh culture medium, and placed in 100 mm cell culture dishes (Falcon-Becton Dickinson Labware, Franklin Lakes, NJ, USA) with DMEM containing 10% fetal bovine serum (FBS), penicillin (100 U/mL), and streptomycin (100 μg/mL) at 37 °C in a 5% CO2 atmosphere. Media were changed every 2 days and the cells were characterized as previously described [23, 24]. Passages 5 to 8 were used for the experiments.

To determine cell type, the cells were cultured on sterile cover glasses, placed into 12-well plates and grown until nearly confluent. The cells were stained with antibody to CD90 (fibroblast marker, R&D Systems Inc., Minneapolis, MN, USA; 1:100), Vimentin (fibroblast marker, R&D Systems Inc., 1:50), NG2 chondroitin sulfate proteoglycan (NG2, pericyte marker, Millipore, San Francisco, CA, USA; 1:50), CD31 (endothelial cell marker, Chemicon, Temecula, CA, USA; 1:50), or DAPI (nucleus marker; Vector Laboratories Inc., Burlingame, CA, USA), as previously described [25]. Signals were visualized and digital images were obtained using a confocal fluorescence microscope (K1-Fluo, Nanoscope Systems, Inc., Daejeon, Korea).

Primary culture of MCPs

The primary cultures of MCPs were performed as described previously [26, 27]. Shortly, 8 weeks old male C57BL/6J mice were anesthetized with ketamine (100 mg/kg) and xylazine (5 mg/kg) intramuscularly, and sacrificed by cervical dislocation. Then, the penis tissues were harvested and maintained in sterile vials with HBSS (Gibco). After washing three times with PBS, the urethra and dorsal neurovascular bundle were removed, and only the corpus cavernosum tissues were used. The corpus cavernosum tissues were cut into approximately 1–2 mm sections and settled via gravity into collagen I-coated 35 mm cell culture dishes with 300 µL complement DMEM (GIBCO) at 37 °C for 20 min in a 5% CO2 atmosphere. Thereafter, 900 µL of complement medium was added and incubated at 37 °C with 5% CO2. The complement medium contained 20% FBS, 1% penicillin/streptomycin, and 10 nM human pigment epithelium-derived factor (PEDF; Sigma-Aldrich). The medium was changed every 2 days, and after approximately 10 days sprouting cells were sub-cultured into collagen I (Advanced BioMatrix, San Diego, CA, USA)-coated dishes. Cells from passages 2 to 3 were used for the experiments.

Preparation and characterization of MCP-derived EV–mimetic nanovesicles (NVs)

MCP-derived EV–mimetic NVs (PC–NVs) were prepared using a mini extruder system (Avanti Polar Lipids, Birmingham, AL, USA), as described previously [28, 29]. Briefly, MCPs were washed three times with PBS, detached with 0.25% Trypsin–EDTA (Invitrogen, Carlsbad, CA, USA) and re-suspended in 4-(2-Hydroxyethyl)piperazine-1-ethanesulfonic acid (HEPES)buffer solution (Gibco). The cell suspension was sequentially extruded 10 times across 10, 5, and 1 μm pore-sized polycarbonate membranes (Nuclepore, Whatman Inc., Clifon, NJ, USA), respectively. Next, ultracentrifugation was performed at 100,000 g for 2 h at 4 °C with a step gradient, which was formed with 50% iodixanol (1 mL; Axis-Shield PoC AS, Oslo, Norway) overlaid with 10% iodixanol (2 mL) and topside with the extruded samples (7 mL). PC–NVs were filtered with a 0.45 μm filter and stored at − 80 °C until further use. To quantify the PC–NVs, the EXOCET exosome quantitation assay kit (System Biosciences, Palo Alto, CA, USA) was used, and 1 µg/µL of the final concentration of the PC–NVs was prepared for all experiments.

Western blotting

For the immunoblot analyses of PC–NVs, equal protein amounts (10 µg) of purified PC–NVs and whole cells extracted using RIPA lysis buffer (Sigma-Aldrich) were separated by SDS-PAGE (12% gel) and transferred to polyvinylidene fluoride (PVDF) membranes. Each blot was blocked and incubated with antibodies to GM130 (BD Biosciences, San Jose, CA, USA; 1:1000), CD9 (Abcam, Cambridge, UK; 1:1000), CD81 (Novus Biologicals; 1:1000), or TSG101 (Novus Biologicals; 1:500).

RNA-sequencing assay

For the RNA-sequencing study, the normal and PD TA-derived fibroblasts were cultured and treated with PC–NVs (n = 4 per group). The RNA-sequencing assay was performed by E-Biogen Inc. (Korea). Briefly, total RNA was isolated 24 h after exposure to PC–NVs using TRIzol reagent (Invitrogen). RNA quality was assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies, Amstelveen, The Netherlands), and RNA quantification was performed using an ND-2000 Spectrophotometer (Thermo Inc., DE, USA).

Library sequencing and data analysis

Libraries were prepared from total RNA using the SMARTer Stranded RNA-Seq Kit (Clontech Laboratories, Inc., USA). The isolation of mRNA was performed using the Poly(A) RNA Selection Kit (LEXOGEN, Inc., Austria). Indexing was performed using the Illumina indices 1–12. The enrichment step was performed using PCR. Subsequently, libraries were checked using the Agilent 2100 Bioanalyzer (DNA High Sensitivity Kit) to evaluate the mean fragment size. Quantification was performed using the library quantification kit using a StepOne Real-Time PCR System (Life Technologies, Inc., USA). High-throughput sequencing was performed as paired-end 100 sequencing using HiSeq 2500 (Illumina, Inc., USA).

Quality control of the raw sequencing data was performed using FastQC ( Adapter and low-quality reads (< Q20) were removed using FASTX_Trimmer ( and BBMap ( Then, the trimmed reads were mapped to the reference genome using TopHat [30]. Gene expression levels were estimated using RC (read count) and FPKM (fragments per kb per million reads) values by BEDTools [31] and Cufflinks [32]. The expression values were normalized with the Quantile normalization method using EdgeR within R ( Data mining and graphic visualization were performed using ExDEGA (E-Biogen, Inc., Korea). The datasets generated and/or analyzed during the current study are available in the Gene Expression Omnibus repository ( Accession no. GSE146500).

Validation of sequencing data by RT-PCR

Total RNA was extracted from cultured cells using TRIzol (Invitrogen) following the manufacturer’s protocols. Reverse transcription was performed using 1 µg of RNA in 20 µL of reaction buffer with oligo dT primer and AccuPower RT Premix (Bioneer Inc., Korea). The PCR reaction was performed with denaturation at 94 °C for 30 s, annealing at 60 °C for 30 s, and extension at 72 °C for 1 min in a DNA Engine Tetrad Peltier Thermal Cycler. For the analysis of PCR products, 10 µL of each PCR productwas electrophoresed on a 1% agarose gel and detected under ultraviolet light. GAPDH was used as an internal control [14].

All digital image, western blot, and PCR band densitometry analyses were performed using an image analyzer system (National Institutes of Health [NIH] ImageJ 1.34,

Statistical analysis

All data are expressed as means ± standard errors. Statistical analysis was performed using Student t-test. p values less than 0.05 were considered statistically significant.


Identification of human fibroblasts

The fibroblasts were isolated from human normal and PD plaque tissues. Representative images showed high positive staining for CD90 and Vimentin (fibroblast markers) of more than 95%, but not for pericyte (NG2) or endothelial cell (CD31) markers (Fig. 1A, B).

Fig. 1
figure 1

Isolation and characterization of fibroblasts from human tunica albuginea (TA) and extracellular vesicle (EV)–mimetic nanovesicles (NVs) from mouse cavernous pericytes (MCPs). A Representative phage images (screen magnification ×40 and ×100) of primary cultured fibroblasts from the TA tissues of healthy subjects and those with PD at passage 5. B Immunofluorescent staining of fibroblasts with antibodies against CD90 (fibroblast marker), Vimentin (fibroblast marker), NG2 (pericyte marker), and CD31 (endothelial cell marker). Nuclei were labeled with the DNA dye DAPI. Scale bar indicates 100 μm. DAPI = 4,6-diamidino-2-phenylindole. C Representative phage images (screen magnification ×40) of primary cultured mouse cavernous pericytes at passage 2. D Representative western blot for EV positive markers (CD9, CD81, and TSG101) or an EV negative marker (GM130) in MCPs and pericyte-derived EV–mimetic NVs

PC–NV preparation and characterization

PC–NVs were prepared from MCPs according to previous methods [28]. Western blot analysis showed that PC–NVs displayed positive exosomes markers, including CD9, CD81, and TSG101, but not for negative marker GM130 (Fig. 1C, D).

Transcriptional profiling and gene ontology (GO) category analysis

For this study, three gene libraries for the normal fibroblast (NF), PD fibroblast (PF), and PC–NVs-treated PF (PFPC) groups were constructed for an RNA-sequencing assay (n = 4 for each group). In total, 25,737 genes were detected in three libraries. Significant gene selection was performed with three conditions: fold-change > 2.0, log2 > 4, and p value < 0.05. Among all detected genes, 3961 showed significant differential expression in the PF group compared with the NF group, and 174 were significantly differentially expressed in the PFPC group compared with the PF group (Fig. 2A–C). Only 91 contra-regulated genes (Additional file 1: Table S1) were detected between PF/NF and PFPC/PF through Venn diagram analysis (Fig. 2D).

Fig. 2
figure 2

Analysis of differentially expressed genes (DEGs) from three libraries. A Total gene expression of normal fibroblasts (NF) and PD fibroblasts (PF). B Total gene expression of PF and PC–NV-treated PF (PFPC) groups. C DEG analysis for significantly altered gene selection according to the conditions set: fold-change > 2.0, log2 > 4, and p value < 0.05. D Venn diagram analysis between PF/NF and PFPC/PF

To evaluate the GO categories, all significantly differentially expressed genes (DEGs) were classified in 16 GO categories by ExDEGA software (E-Biogen, Inc., Korea). Among these, extracellular matrix (33.61%), angiogenesis (32.33%), and fibrosis (31.07%) made up the largest proportion in the PF group compared with the NF group (Fig. 3A, B). However, fibrosis (3.88%), extracellular matrix (2.92%), and angiogenesis (2.59%) made up the largest proportion in the PFPC group compared with the PF group (Fig. 3C, D).

Fig. 3
figure 3

Significantly differentially expressed genes (DEGs) of the RNA-sequencing data were distributed to gene ontology (GO) categories. A, B The distribution (A total percentage; B detailed percentage and numbers of upregulated and downregulated genes, respectively) of total significantly DEGs in 16 GO categories in the PD fibroblast (PF) group compared with the normal fibroblast (NF) group. C, D The distribution (C total percentage; D detailed percentage and numbers of upregulated and downregulated genes, respectively) of total significantly DEGs in 16 GO categories in the PC–NV-treated PF group compared with the PF group

Validation of RNA-sequencing results by RT-PCR

To validate the RNA-sequencing results, we selected 20 genes (Additional file 1: Table S2) from 91 contra-regulated DEGs, and 11 (primers as shown in Additional file 1: Table S3) showed results consistent with the RNA-sequencing assay by RT-PCR. Among these genes, MMP3, AKR1C1, SMOC1, ANGPTL2, SEMA3A, TRIM15, EGR1, and BMP2 were downregulated in the PF group compared with the NF group, and were significantly recovered in the PFPC group (Fig. 4A, C). Only TFPI2, SFRP4, and SERPINE1 were induced in the PF group and recovered in the PFPC group (Fig. 4B, D).

Fig. 4
figure 4

RT-PCR validation of differentially expressed genes selected from the RNA-sequencing assay. A, B Eleven contra-regulated genes were found to be consistent with the RNA-sequencing result in normal fibroblast (NF), PD fibroblast (PF), and PC–NV-treated PF (PFPC) groups. C, D Each bar depicts the mean value (± standard error) from three separate experiments. *p < 0.001 compared with the NF group. #p < 0.05 compared with the PF group


The accurate physiological and pathological mechanisms of PD remain poorly understood. To date, most gene expression studies have focused on human PD plaque at a tissue level in vivo [33, 34] but not at a cellular level in vitro. Therefore, to investigate the exact mechanisms and potential target genes for PD, we cultured human fibroblasts from human PD plaque and performed RNA-sequencing assays.

EVs display a potential role in kidney fibrosis and other fibrotic diseases [20, 35]; however, little is currently known regarding the detailed mechanisms. Considering the low yield of EVs, we extracted more than 100-fold greater EV–mimetic NVs from MCPs, which were primarily cultured from mouse corpus cavernosum tissues. Many previous studies have found that pericytes display diverse features in relation to fibrosis that are dependent on different molecular targets [8, 9, 11]. We used pericytes isolated from mouse cavernous tissue, but not a human tissue, because it was not easy to harvest human erectile tissue in a large quantity. Although there are animal models of PD by using intratunical injection of fibrin [36] or transforming growth factor-β1 [37, 38], these animal models cannot represent the complex pathophysiologic processes of the human PD. Meanwhile, we can easily obtain fibrotic plaque from the patients with PD during reconstructive surgery. In the present study, therefore, human PD fibroblasts exposed to mouse origin of PC–NVs were compared with human PD fibroblasts to investigate the regulation of gene expression by PC–NVs in PD.

From the RNA-sequencing assay, 3961 DEGs were detected, and the 16 top GO categories were assessed in this study. GO analysis showed that significantly altered genes were enriched in the extracellular matrix, angiogenesis, and fibrosis. The extracellular matrix is a driver of progressive fibrosis [39], and angiogenesis is closely associated with chronic liver fibrosis [20]. These data suggest that our DEG detection is credible. However, the molecular basis of PC–NVs in regulating the extracellular matrix or angiogenesis pathway in PD remains largely unknown. In this study, only 91 contra-regulated genes were identified from the three libraries (NF, PF, and PFPC). After precision screening, 20 genes were selected and validated by RT-PCR in same conditions. However, only 11 genes were validated to be consistent with the RNA-sequencing results as showed in Fig. 4. Through literature analysis, we found that most target genes have multiple roles, such as angiogenesis, extracellular matrix (ECM), and fibrosis. Some studies have shown that the MMP3 attenuated the process of fibrosis by degrading α2-antiplasmin and ECM [40] and these enzymes inhibitors are up-regulated (TFPI2) [41]. Many studies have shown that fibrosis and angiogenesis have common characteristics. Wound healing is a typical example of the synergistic interaction between the molecular mechanisms of fibrosis and angiogenesis [42, 43]. Although some studies have shown that AKR1C1 [44], SMOC1 [45], ANGPTL2 [46], SEMA3A [47], TRIM15 [48], EGR1 [49], BMP2 [50], SFRP4 [51], and SERPINE1 [52] are related to angiogenesis or fibrosis, but the detailed mechanism of action of these genes is still unknown, and it will be interesting to study the network and pathway of these genes in fibrosis and angiogenesis in the future. Therefore, on the one hand, we hypothesized that pericytes might be releasing nanovesicles, which contain these target genes to directly affect nearby fibroblasts, such as increase angiogenesis and reducing excessive pathological ECM deposition during fibrosis. On the other hand, pericyte-derived nanovesicles might indirectly affect the expression of these target genes in fibroblasts through microRNA, and further research is needed to resolve these hypotheses. These genes may be the key to understanding how PC–NVs regulate the extracellular matrix, angiogenesis, and fibrosis mechanisms in PD.

To the best of our knowledge, this is the first study to demonstrate the systematic profiling of gene alterations in NF, PF, and PFPC. However, the present study has some limitations. First, a small number of cultured human fibroblast samples were used in target gene validation and age differences existed among groups, and need to collect more age-matched samples for further studies. It would be more appropriate to use healthy TA from an age-matched PD patient instead of using TA from a young patient with congenital penile curvature. Second, we were unable to demonstrate the network of these validated genes in the extracellular matrix, angiogenesis, and fibrosis pathways. Further studies are needed to document the function of each gene at cellular and molecular levels as well as in animal models of PD. Third, mouse corpus cavernous pericytes were used for EV–mimetic NVs isolation, and further research is needed to cultivate human corpus cavernous pericytes and isolate human PC–NVs and process them in human PD fibroblasts to study the detailed mechanism of these initial selected target genes.


In summary, we profiled the DEGs of human TA cultured fibroblasts in NF, PF, and PFPC groups. We hypothesize that these validated genes are good candidates for the study of the mechanism of PC–NVs in PD. Further studies exploring the effect of these target genes will be beneficial to further our understanding of the detailed mechanisms of the extracellular matrix, angiogenesis, and fibrosis in PD.

Availability of data and materials

The datasets used during the current study available from the corresponding author on request.



Peyronie’s disease


Extracellular vesicle




Pericyte-derived EV–mimetic NVs


Reverse transcription polymerase chain reaction


Tunica albuginea


Rrectile dysfunction


Mouse cavernous pericytes


Hank’s balanced salt solution


Phosphate-buffered saline


Dulbecco’s modified Eagle’s medium


Fetal bovine serum


Pigment epithelium-derived factor


4-(2-Hydroxyethyl)piperazine-1-ethanesulfonic acid


Polyvinylidene fluoride


  1. Gholami SS, Gonzalez-Cadavid NF, Lin CS, Rajfer J, Lue TF. Peyronie’s disease: a review. J Urol. 2003;169:1234–41.

    Article  PubMed  Google Scholar 

  2. Hellstrom WJ, Bivalacqua TJ. Peyronie’s disease: etiology, medical, and surgical therapy. J Androl. 2000;21:347–54.

    CAS  PubMed  Google Scholar 

  3. Duncan MR, Berman B, Nseyo UO. Regulation of the proliferation and biosynthetic activities of cultured human Peyronie’s disease fibroblasts by interferons-alpha, -beta and -gamma. Scand J Urol Nephrol. 1991;25:89–94.

    Article  CAS  PubMed  Google Scholar 

  4. Russo GI, Milenkovic U, Hellstrom W, Levine LA, Ralph D, Albersen M. Clinical efficacy of injection and mechanical therapy for Peyronie’s disease: a systematic review of the literature. Eur Urol. 2018;74:767–81.

    Article  PubMed  Google Scholar 

  5. Levine LA, Greenfield JM, Estrada CR. Erectile dysfunction following surgical correction of Peyronie’s disease and a pilot study of the use of sildenafil citrate rehabilitation for postoperative erectile dysfunction. J Sex Med. 2005;2:241–7.

    Article  PubMed  Google Scholar 

  6. Armulik A, Abramsson A, Betsholtz C. Endothelial/pericyte interactions. Circ Res. 2005;97:512–23.

    Article  CAS  PubMed  Google Scholar 

  7. Armulik A, Genove G, Betsholtz C. Pericytes: developmental, physiological, and pathological perspectives, problems, and promises. Dev Cell. 2011;21:193–215.

    Article  CAS  PubMed  Google Scholar 

  8. Crisan M, Yap S, Casteilla L, et al. A perivascular origin for mesenchymal stem cells in multiple human organs. Cell Stem Cell. 2008;3:301–13.

    Article  CAS  PubMed  Google Scholar 

  9. Cai X, Lin Y, Friedrich CC, et al. Bone marrow derived pluripotent cells are pericytes which contribute to vascularization. Stem Cell Rev Rep. 2009;5:437–45.

    Article  PubMed  Google Scholar 

  10. Lin SL, Kisseleva T, Brenner DA, Duffield JS. Pericytes and perivascular fibroblasts are the primary source of collagen-producing cells in obstructive fibrosis of the kidney. Am J Pathol. 2008;173:1617–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Yoshiji H, Kuriyama S, Yoshii J, et al. Vascular endothelial growth factor and receptor interaction is a prerequisite for murine hepatic fibrogenesis. Gut. 2003;52:1347–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Patsenker E, Popov Y, Stickel F, et al. Pharmacological inhibition of integrin alphavbeta3 aggravates experimental liver fibrosis and suppresses hepatic angiogenesis. Hepatology. 2009;50:1501–11.

    Article  CAS  PubMed  Google Scholar 

  13. Stockmann C, Kerdiles Y, Nomaksteinsky M, et al. Loss of myeloid cell-derived vascular endothelial growth factor accelerates fibrosis. Proc Natl Acad Sci U S A. 2010;107:4329–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Yin GN, Jin HR, Choi MJ, et al. Pericyte-derived Dickkopf2 regenerates damaged penile neurovasculature through an angiopoietin-1-Tie2 pathway. Diabetes. 2018;67:1149–61.

    Article  CAS  PubMed  Google Scholar 

  15. Yin GN, Park SH, Ock J, et al. Pericyte-derived extracellular vesicle–mimetic nanovesicles restore erectile function by enhancing neurovascular regeneration in a mouse model of cavernous nerve injury. J Sex Med. 2020;17:2118–28.

    Article  CAS  PubMed  Google Scholar 

  16. Colombo M, Raposo G, Thery C. Biogenesis, secretion, and intercellular interactions of exosomes and other extracellular vesicles. Annu Rev Cell Dev Biol. 2014;30:255–89.

    Article  CAS  PubMed  Google Scholar 

  17. Record M, Carayon K, Poirot M, Silvente-Poirot S. Exosomes as new vesicular lipid transporters involved in cell–cell communication and various pathophysiologies. Biochim Biophys Acta. 2014;1841:108–20.

    Article  CAS  PubMed  Google Scholar 

  18. Zhang ZG, Chopp M. Exosomes in stroke pathogenesis and therapy. J Clin Investig. 2016;126:1190–7.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Li K, Chen Y, Li A, Tan C, Liu X. Exosomes play roles in sequential processes of tumor metastasis. Int J Cancer. 2019;144:1486–95.

    Article  CAS  PubMed  Google Scholar 

  20. Zhang W, Zhou X, Zhang H, Yao Q, Liu Y, Dong Z. Extracellular vesicles in diagnosis and therapy of kidney diseases. Am J Physiol Renal Physiol. 2016;311:F844–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Kim OY, Lee J, Gho YS. Extracellular vesicle mimetics: novel alternatives to extracellular vesicle-based theranostics, drug delivery, and vaccines. Semin Cell Dev Biol. 2017;67:74–82.

    Article  CAS  PubMed  Google Scholar 

  22. Jang SC, Kim OY, Yoon CM, et al. Bioinspired exosome-mimetic nanovesicles for targeted delivery of chemotherapeutics to malignant tumors. ACS Nano. 2013;7:7698–710.

    Article  CAS  PubMed  Google Scholar 

  23. Piao S, Choi MJ, Tumurbaatar M, et al. Transforming growth factor (TGF)-beta type I receptor kinase (ALK5) inhibitor alleviates profibrotic TGF-beta1 responses in fibroblasts derived from Peyronie’s plaque. J Sex Med. 2010;7:3385–95.

    Article  CAS  PubMed  Google Scholar 

  24. Ryu JK, Kim WJ, Choi MJ, et al. Inhibition of histone deacetylase 2 mitigates profibrotic TGF-beta1 responses in fibroblasts derived from Peyronie’s plaque. Asian J Androl. 2013;15:640–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Yin GN, Ryu JK, Kwon MH, et al. Matrigel-based sprouting endothelial cell culture system from mouse corpus cavernosum is potentially useful for the study of endothelial and erectile dysfunction related to high-glucose exposure. J Sex Med. 2012;9:1760–72.

    Article  CAS  PubMed  Google Scholar 

  26. Neng L, Zhang W, Hassan A, et al. Isolation and culture of endothelial cells, pericytes and perivascular resident macrophage-like melanocytes from the young mouse ear. Nat Protoc. 2013;8:709–20.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Yin GN, Das ND, Choi MJ, et al. The pericyte as a cellular regulator of penile erection and a novel therapeutic target for erectile dysfunction. Sci Rep. 2015;5:10891.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Kwon MH, Song KM, Limanjaya A, et al. Embryonic stem cell-derived extracellular vesicle–mimetic nanovesicles rescue erectile function by enhancing penile neurovascular regeneration in the streptozotocin-induced diabetic mouse. Sci Rep. 2019;9:20072.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Oh K, Kim SR, Kim DK, et al. In vivo differentiation of therapeutic insulin-producing cells from bone marrow cells via extracellular vesicle–mimetic nanovesicles. ACS Nano. 2015;9:11718–27.

    Article  CAS  PubMed  Google Scholar 

  30. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25:1105–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Roberts A, Trapnell C, Donaghey J, Rinn JL, Pachter L. Improving RNA-Seq expression estimates by correcting for fragment bias. Genome Biol. 2011;12:R22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Wang Z, Liu JF, Zhou ZH, et al. Gene expression profiles and effects of transforming growth factor-beta1 intervention in Peyronie’s disease. Zhonghua Wai Ke Za Zhi. 2004;42:182–6.

    PubMed  Google Scholar 

  34. Sampaio FJ. Comparison of gene expression profiles between Peyronie’s disease and Dupuytren’s contracture. Int Braz J Urol. 2004;30:349–50.

    Article  PubMed  Google Scholar 

  35. Asef A, Mortaz E, Jamaati H, Velayati A. Immunologic role of extracellular vesicles and exosomes in the pathogenesis of cystic fibrosis. Tanaffos. 2018;17:66–72.

    PubMed  PubMed Central  Google Scholar 

  36. Davila HH, Ferrini MG, Rajfer J, Gonzalez-Cadavid NF. Fibrin as an inducer of fibrosis in the tunica albuginea of the rat: a new animal model of Peyronie’s disease. BJU Int. 2003;91:830–8.

    Article  CAS  PubMed  Google Scholar 

  37. El-Sakka AI, Hassoba HM, Chui RM, Bhatnagar RS, Dahiya R, Lue TF. An animal model of Peyronie’s-like condition associated with an increase of transforming growth factor beta mRNA and protein expression. J Urol. 1997;158:2284–90.

    Article  CAS  PubMed  Google Scholar 

  38. Piao S, Ryu JK, Shin HY, et al. Repeated intratunical injection of adenovirus expressing transforming growth factor-beta1 in a rat induces penile curvature with tunical fibrotic plaque: a useful model for the study of Peyronie’s disease. Int J Androl. 2008;31:346–53.

    Article  CAS  PubMed  Google Scholar 

  39. Bruno S, Porta S, Bussolati B. Extracellular vesicles in renal tissue damage and regeneration. Eur J Pharmacol. 2016;790:83–91.

    Article  CAS  PubMed  Google Scholar 

  40. Niwa H, Kanno Y, Shu E, Seishima M. Decrease in matrix metalloproteinase3 activity in systemic sclerosis fibroblasts causes alpha2antiplasmin and extracellular matrix deposition, and contributes to fibrosis development. Mol Med Rep. 2020;22:3001–7.

    CAS  PubMed  Google Scholar 

  41. Hanahan D, Folkman J. Patterns and emerging mechanisms of the angiogenic switch during tumorigenesis. Cell. 1996;86:353–64.

    Article  CAS  PubMed  Google Scholar 

  42. Folkman J, Shing Y. Angiogenesis. J Biol Chem. 1992;267:10931–4.

    Article  CAS  PubMed  Google Scholar 

  43. Kalluri R, Sukhatme VP. Fibrosis and angiogenesis. Curr Opin Nephrol Hypertens. 2000;9:413–8.

    Article  CAS  PubMed  Google Scholar 

  44. Zhu H, Chang LL, Yan FJ, et al. AKR1C1 activates STAT3 to promote the metastasis of non-small cell lung cancer. Theranostics. 2018;8:676–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Awwad K, Hu J, Shi L, et al. Role of secreted modular calcium-binding protein 1 (SMOC1) in transforming growth factor beta signalling and angiogenesis. Cardiovasc Res. 2015;106:284–94.

    Article  CAS  PubMed  Google Scholar 

  46. Carbone C, Piro G, Merz V, et al. Angiopoietin-like proteins in angiogenesis, inflammation and cancer. Int J Mol Sci. 2018;19:431.

    Article  PubMed Central  CAS  Google Scholar 

  47. Chen G, Sima J, Jin M, et al. Semaphorin-3A guides radial migration of cortical neurons during development. Nat Neurosci. 2008;11:36–44.

    Article  PubMed  CAS  Google Scholar 

  48. Kornberg LJ, Shaw LC, Spoerri PE, Caballero S, Grant MB. Focal adhesion kinase overexpression induces enhanced pathological retinal angiogenesis. Investig Ophthalmol Vis Sci. 2004;45:4463–9.

    Article  Google Scholar 

  49. Sheng J, Liu D, Kang X, Chen Y, Jiang K, Zheng W. Egr-1 increases angiogenesis in cartilage via binding Netrin-1 receptor DCC promoter. J Orthop Surg Res. 2018;13:125.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Yang YL, Ju HZ, Liu SF, et al. BMP-2 suppresses renal interstitial fibrosis by regulating epithelial–mesenchymal transition. J Cell Biochem. 2011;112:2558–65.

    Article  CAS  PubMed  Google Scholar 

  51. Muley A, Majumder S, Kolluru GK, et al. Secreted frizzled-related protein 4: an angiogenesis inhibitor. Am J Pathol. 2010;176:1505–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Wu J, Strawn TL, Luo M, et al. Plasminogen activator inhibitor-1 inhibits angiogenic signaling by uncoupling vascular endothelial growth factor receptor-2-alphaVbeta3 integrin cross talk. Arterioscler Thromb Vasc Biol. 2015;35:111–20.

    Article  CAS  PubMed  Google Scholar 

Download references


The manuscript is edited by a professional, native English-speaking editor at Editage (; Vikas Narang, Chief Operating Officer, Editage).


The design, data collection, device and PC–NVs isolation system of this study were supported by a National Research Foundation of Korea (NRF) grant (Guo Nan Yin, 2018R1C1B6003829) and the interpretation of data, writing the manuscript were supported by a Medical Research Center Grant (Ji-Kan Ryu, 2014R1A5A2009392) funded by the Korean government (Ministry of Science, ICT and Future Planning), and the purchase of experimental animals and regents of this study was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (Mi-Hye Kwon, 2018R1A6A3A11050852).

Author information

Authors and Affiliations



GNY, SGP, JKS and JKR drafted and designed the research. GNY, SGP, JYO, MHK performed the experiments. GNY, SGP, ZYL, LW, JYO, DKK, analyzed the data. GNY, SGP prepared the figures and drafted the manuscript. YSG, JKS and JKR edited and revised the manuscript. GNY, SGP, ZYL, LW, JYO, MHK, DKK, YSG, JKS and JKR have read and approved the final version of the manuscript. All authors are accountable for all aspects of work. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Jun-Kyu Suh or Ji-Kan Ryu.

Ethics declarations

Ethics approval and consent to participate

All TA tissues donors provided written informed consent and the procedures were approved by the Institutional Review Board (IRB No: 2007-730). And all male C57BL/6J mice (8 weeks old, Orient Bio, Korea) used in this study were approved by the Institutional Animal Care and Use Committee of our University (approval number: 171129-527), respectively. All authors of the manuscript have read and agreed to its content of the abstract.

Consent for publication

All participating authors agree to publish in “BMC urology”.

Competing interests

The authors declare 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.

Supplementary tables for contra-regulated genes, primers for RT-PCR and un-cutted membrane or gels images.

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 The Creative Commons Public Domain Dedication waiver ( 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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yin, G.N., Piao, S., Liu, Z. et al. RNA-sequencing profiling analysis of pericyte-derived extracellular vesicle–mimetic nanovesicles-regulated genes in primary cultured fibroblasts from normal and Peyronie’s disease penile tunica albuginea. BMC Urol 21, 103 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Peyronie’s disease
  • Pericytes
  • Nanovesicles
  • Gene expression
  • RNA-sequencing