# Dysregulation of the endothelin pathway in lymphangioleiomyomatosis with no direct effect on cell proliferation and migration

Nader Chebib, Fabienne Archer, Alexandra Bobet Erny, Caroline Leroux, Vincent Cottin

► To cite this version:

Nader Chebib, Fabienne Archer, Alexandra Bobet Erny, Caroline Leroux, Vincent Cottin. Dysregulation of the endothelin pathway in lymphangioleiomyomatosis with no direct effect on cell proliferation and migration. Scientific Reports, 2018, 8, 10.1038/s41598-018-32795-3. hal-01959655

HAL Id: hal-01959655

https://hal.science/hal-01959655

Submitted on 18 Dec 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# Dysregulation of the endothelin pathway in lymphangioleiomyomatosis with no direct effect on cell proliferation and migration

Nader Chebib<fnref n="1" />,², Fabienne Archer<fnref n="1" />, Alexandra Bobet-Erny<fnref n="1" />, Caroline Leroux ID <fnref n="1" /> & Vincent Cottin<fnref n="1" />,²

Received: 18 April 2018

Accepted: 12 September 2018

Published online: 02 October 2018

LAM is a rare low-grade metastasizing lung neoplasm. Inhibitors of mTOR improve clinical outcome of LAM patients by preventing loss of lung function. Nevertheless, other cell targets may be of interest for drug development. Therefore, we explored the potential role of EDN1 (endothelin) in LAM. We report an increased endothelin blood level in LAM patients as well as EDN1 overexpression and EDN1 receptor downregulation in LAM-derived primary cells and in TSC2<sup>NEG</sup> cells mutated in TSC2. We evidenced EDN pathway dysregulation based on EDN1, EDNRA, EDNRB and ARRB1 mRNA expression in LAM-derived primary cells. We showed overexpression of EDN1 and ARRB1 mRNAs in TSC2<sup>NEG</sup> cells; these cells lost their ability to respond to stimulation by endothelin. We analyzed the effects of endothelin receptor antagonists alone or in combination with rapamycin, an mTOR inhibitor, on proliferation and migration of LAM cells. Rapamycin treatment of TSC2<sup>NEG</sup> cells significantly reduced cell proliferation or migration, while none of the tested inhibitors of EDN receptors impaired these functions. We showed that TSC2<sup>NEG</sup> cells have acquired a transformed phenotype as showed by their ability to grow as spheroids in semi-solid medium and that unlike endothelin receptors antagonists, rapamycin reduced anchorage-independent cell growth and prevented expansion of TSC2<sup>NEG</sup> spheroids.

Lymphangioleiomyomatosis (LAM) is a rare pulmonary disease mainly affecting young women<fnref n="1" />. LAM can occur as an isolated disorder, defined as sporadic LAM or in patients with tuberous sclerosis complex, a genetic disease characterized by mutations of the *TSC1* and *TSC2* (*Tuberous Sclerosis Complex*) tumor suppressor genes respectively encoding hamartin and tuberin². LAM cells carry mutations of the *TSC2* gene³, inducing constitutive activation of the PI3K/Akt/mTOR pathway and LAM cell proliferation. LAM causes cystic destruction of the lungs and development of benign renal tumors or angiomyolipomas<fnref n="1" />. Two cell populations are present in LAM lesions: the myofibroblastic-like cells that express markers of smooth muscle cells and fibroblasts, such as α-smooth muscle actin (α-SMA), vimentin and desmin⁴,⁵ and the epithelioid-like cells that express melanocytic markers such as MLANA (Melan A) and proteins evidenced with HMB45 and PNL2 antibodies⁵,⁶. In LAM patients, circulating VEGF-D (Vascular Endothelial Growth Factor D) is increased in the blood and is associated with lymphangio-genesis, a major pathogenic mechanism in LAM progression⁷,⁸.

LAM is considered as a low-grade, destructive, metastasizing neoplasm⁹. Circulating LAM cells have been found in the blood, urine and chylous effusions¹⁰,¹¹ of LAM patients. LAM cells invade organs through degradation of the extracellular matrix by metalloproteinases, similarly to metastatic cancer cells¹²,¹³.

Although mTOR inhibitors (everolimus, sirolimus) have been shown to improve clinical outcomes in preventing loss of lung function¹⁴,¹⁵ and have been approved to treat LAM, other pathways must be explored to improve

<footnote n="1">¹Université de Lyon, Université Claude Bernard Lyon 1, INRA, EPHE, IVPC, Viral Infections and Comparative Pathology, UMR754, F69007, Lyon, France. ²Hospices Civils de Lyon, Groupement Hospitalier Est, Department of Respiratory Diseases, National Reference Center for Rare Pulmonary Diseases, Lyon, France. Nader Chebib and Fabienne Archer contributed equally. Caroline Leroux and Vincent Cottin jointly supervised this work. Correspondence and requests for materials should be addressed to C.L. (email: caroline.leroux@univ-lyon1.fr)</footnote>
patient treatment. In human cancer cells, high expression levels of EDN1 (Endothelin 1) and of endothelin receptors A and B (EDNRA and EDNRB) are associated with the increase of circulating VEGF and of microvessel density¹⁶⁻¹⁹. The EDN1/EDNR/ARRB1 (β Arrestin 1) pathway is implicated in cell proliferation, migration, invasion, survival and angiogenesis in several diseases, among them lung, ovary, prostate and breast cancers²⁰,²¹. The development of endothelin receptor antagonists (ERAs) such as bosentan, a dual EDNRA and EDNRB receptor antagonist, or BQ123 targeting EDNRA, provided targeted treatments for pulmonary arterial hypertension and cancer²²⁻²⁶.

In this study, we explored the role of EDN1 and of its receptors in LAM-derived primary cells and in angiomyolipoma-derived cells lines. We report an increased blood level of endothelin in LAM patients as compared to controls, and the overexpression of EDN1 and downregulation of its receptors in LAM-derived primary cells as well as in TSC2<sup>NEG</sup> cell lines. We analyzed the effects of ERAs, alone or in combination with rapamycin, on LAM cell proliferation and migration.

## Materials and Methods

**Cell lines.** The 621-101 TSC2<sup>NEG</sup> and 621-103 TSC2<sup>POS</sup> cell lines (respectively named “TSC2<sup>NEG</sup>” and “TSC2<sup>POS</sup>” along our study) were generously provided by Pr E.P. Henske (Boston, United States)²⁷. The TSC2<sup>NEG</sup> cell line was derived from a renal angiomyolipoma of a LAM patient. They carry a missense mutation in exon 16 of the *TSC2* gene (G1832A) leading to a loss of heterozygosity. The TSC2<sup>POS</sup> cell line has been developed by re-expression of normal TSC2 gene in the 621-101 TSC2<sup>NEG</sup> cells. These cell lines were cultured in DMEM medium (Sigma) supplemented with 10% inactivated fetal calf serum (Gibco), 100 U/mL penicillin (Sigma), 100 µg/mL streptomycin and with 50 µg/ml zeocin (Thermo Fisher) for the TSC2<sup>POS</sup> cells to maintain the selective pressure for TSC2 expression.

Human primary PASMC (Pulmonary Artery Smooth Muscle Cells) (Lonza) were used as controls and maintained for a short time in culture as recommended.

**Lung-derived primary LAM cells.** LAM pulmonary tissues and associated data from five patients (1300, 1444, 1720, 2634, 2749) were obtained from the Cardiobiotec biobank (CRB-HCL Hospices Civils de Lyon BB-0033-00046), a center for biological resources authorized by the French Ministry of Social Affairs and Health. All samples were collected and used in accordance with the ethical rules of the Biobank and in agreement with the French legislation. All patients signed a written informed consent. After surgical removal, tissue samples were immediately incubated in DMEM medium containing 100 U/mL penicillin, 100 µg/mL streptomycin, and 2.5 µg/mL amphotericin B (PAA). Tissues were dissociated in DMEM medium with 4.5 g/L glucose, 100 U/mL penicillin, 100 µg/mL streptomycin, 2.5 µg/mL amphotericin B, 2 mg/mL type Ia collagenase, 0.1 mg/mL soybean and 3 mg/mL elastase for 60 minutes at 37 °C under gentle agitation. After filtration on gauze, the suspension was centrifuged at 340 g for 10 minutes at 4 °C then cell pellets were cultured in DMEM/HAM'S F12 (50:50, V/V) medium (Gibco) supplemented with 10% inactivated fetal calf serum, 100 U/mL penicillin, 100 µg/mL streptomycin and 2.5 µg/mL amphotericin B.

**Phenotypic characterization of the LAM-derived cells.** Cells fixed with ice-cold acetone were labeled overnight at 4 °C with mouse PNL2 (1/25, Diagnostic Biosystem, Mob421-05), mouse HMB45 (1/50, Diagnostic Biosystem, Mob079-05), mouse anti-desmin (1/100, Diagnostic Biosystem, Mob060-05), mouse anti-αSMA (1/1000, Sigma, A2527) and mouse anti-vimentin (1/200, Sigma, V2258) antibodies then incubated for 1 hour with donkey anti-Goat IgG-DyLight® 488, (1/500, Abcam, ab96931) or mouse anti-IgG-Alexa488 (1/500, Invitrogen A1101) antibodies. Cells were stained with 1 µg/mL DAPI (4'-6' diamidino-2-phenylindole), observed with an “AxioImager Z1 epifluorescence” microscope (Zeiss) and analyzed with the “Zen software” (Zeiss).

**Circulating EDN1 and VEGF-D in LAM patients.** Sera from 28 LAM and 12 control patients with unrelated disorders (6 multiple sclerosis, 6 myocardial infarction) were obtained from the Cardiobiotec biobank. EDN1 and VEGF-D were measured in the serum of LAM patients and controls using “solid-phase ELISA Quantikine Endothelin-1” and “VEGF-D immunoassay” (R&D systems) as recommended.

**Expression of EDNRA and EDNRB proteins.** Total cellular proteins were extracted from cell pellets with 20 mM Tris pH 8; 10% glycerol; 150 mM NaCl; 1% Triton X-100; 5 mM EDTA; 1 mM Na₂VO₄; 1 mM PMSF; 10 µg/mL aprotinin; 10 µg/mL leupeptin, then homogenized with four sonication cycles (20 seconds at 300 W). EDNRA and ENRB receptors were detected on nitrocellulose membrane (Biorad), pre-incubated in TBST-milk (Tris pH 7.6, NaCl, 0.05% Tween 20, 5% non-fat dry milk), by overnight incubation at 4 °C with rabbit anti-EDNRA (SantaCruz, sc-33535) or rabbit anti-EDNRB (SantaCruz, sc-33537) antibodies followed by incubation with anti-rabbit peroxidase conjugated IgG antibodies (Sigma, A0545). The β-actin was detected with mouse anti-β-actin-peroxidase antibody (Sigma, A3854, 1/50 000). Immunoreactive bands were revealed using “West Pico chemiluminescent substrate” (Thermoscientific).

**Cell proliferation assay.** TSC2<sup>NEG</sup> and TSC2<sup>POS</sup> cells were plated onto 96 wells plate at 20,000 cells/well and treated with 10 µM bosentan (Selleckchem), 10 µM BQ123 (Sigma-Aldrich), and/or 10 nM rapamycin (Calbiochem). Cellular proliferation was measured 1 and 3 days post treatment with the “CellTiter-Glo kit” (Promega) as recommended. Each experiment was done in triplicate.

**Cell migration assay.** Forty-eight hours serum-deprived TSC2<sup>NEG</sup> and TSC2<sup>POS</sup> cells were plated in the top chamber of 8 µm pore diameter Fluoroblock cell culture inserts (Corning) in serum-free medium, while the bottom chamber was filled with DMEM medium supplemented with 10% fetal calf serum. Cells were cultivated with 10 nM rapamycin, 10 µM BQ123 or 10 µM bosentan. Fifteen hours after plating, cells were stained
<table><caption>Table 1. PCR primers used in the study. *Reported in<sup>23</sup>.</caption><thead><tr><th>Gene</th><th>Approved name</th><th>5'-3' sequences (Forward Reverse)</th><th>Amplicon size (pb)</th><th>GenBank accession #</th></tr></thead><tbody><tr><td>YWHAZ</td><td>Tyrosine 3-monooxygenase/ tryptophan 5-monooxygenase activation protein zeta</td><td>F AACTCCCCTGAGAAAGCCTG<br/>R CCGATGTCCACAATGTCAAGTT</td><td>145</td><td>NC_019466</td></tr><tr><td>PPIA</td><td>Peptidylprolyl isomerase A</td><td>F GTGCCAGGGTGGTGACTTC<br/>R CCACATGCTTGCCATCCAAC-5'</td><td>200</td><td>XM_004013990</td></tr><tr><td>RPL32</td><td>Ribosomal protein L32</td><td>F GCACATGCTGCCAGTGGCT<br/>R CGCAGCCTGGCATTGGGGTT</td><td>189</td><td>NM_000994.3</td></tr><tr><td>ARRB1</td><td>Arrestin beta 1</td><td>F CGTGTGACTGGTGGGATGAA'<br/>R GCACAGAATTCCGCTTTGTGG</td><td>183</td><td>NM_004041.4</td></tr><tr><td>EDN1</td><td>Endothelin 1</td><td>F GCAGGAAAAGAACTCAGGGCTGAA<br/>R TGCCTTTCAGCTTGGGATCATGAAA</td><td>220</td><td>NM_001955.4</td></tr><tr><td>EDNRA</td><td>Endothelin receptor A</td><td>F ACGAGATGGACAAGAACCGATGTGA<br/>R TCCGTTCATGGGGACCGAGGT</td><td>206</td><td>NM_001957</td></tr><tr><td>EDNRB</td><td>Endothelin receptor B</td><td>F ACCTCAGCAGGATTCTGAAGCTCA<br/>R TCAAATGACTGGCACCAGCAGCA</td><td>209</td><td>NM_000115</td></tr><tr><td>MLANA</td><td>Melan-A</td><td>F TGTGCCTTAACAAGAAGATGCCAC<br/>R CTGCAGAACAGTCACCACCACCTT</td><td>145</td><td>NM_005511</td></tr><tr><td>PMEL</td><td>Premelanosome protein</td><td>F CGCTGATCGTGGGCATCTTGC<br/>R CGGGGTAGACGCAGCCAGTG</td><td>136</td><td>NM_001200054</td></tr><tr><td>TYRP1</td><td>Tyrosinase related protein 1</td><td>F ATGGCCAAGTCGGGAGTTTAGTGT<br/>R TCCATACTGCGTCTGGCACGAA</td><td>123</td><td>NM_000550</td></tr><tr><td>ACTA2</td><td>Actin, alpha 2, smooth muscle, aorta</td><td>F TCCGCTGCCCAGAGACCCTG<br/>R GGTGCCCCCTGATAGGACATTGTT</td><td>146</td><td>NM_001141945</td></tr><tr><td>*TSC2<br/>Exon 40/41</td><td rowspan="2">TSC complex subunit 2</td><td>F ATGGAGGGCCTTGTGGACAC<br/>R CGGAGCCGCTTGATGTG</td><td>249</td><td>NC_000016.10</td></tr><tr><td>*TSC2<br/>Exon 6</td><td>F GGAGATGTAGATTCGGCGTC<br/>R CTGCGGAGCTGAACTTAGG</td><td>222</td><td></td></tr></tbody></table>

with calcein-AM fluorescent dye (Corning) and the fluorescence measured with a “Victor II” spectrophotometer (Wallac) using excitation filter at 488 nm and emission filter at 510 nm.

**Expression of EDN1, EDNRA, EDNRB and ARRB1 mRNAs.** Total cellular RNAs were extracted with the “PureLink RNA mini kit” (Ambion), treated with the “Turbo DNA-free kit” (Ambion), and reverse transcribed with the “iScript cDNA synthesis” kit (Biorad). The mRNA expression was quantified in triplicate with the “PCR KAPA SYBR FAST qPCR” kit (Clinisciences) using 20 ng cDNA. Primer efficiency was determined by amplification of 10-fold dilutions of cDNAs. PPIA (Peptidylprolyl Isomerase A), YWHAZ (Tyrosine 3 monooxygenase) and RPL32 (Ribosomal Protein L32) (Table 1) have been selected as reference genes based on their stability in LAM samples evaluated with the geNorm software. Normalized relative expression was calculated with the 2<sup>-ΔΔCq</sup> method<sup>28</sup>.

**TSC2 mutations.** Mutations in exons 40–41 region (nt 13535–13756, GenBank # NC_000016.10) and exon 6 splicing site (nt 46232–46500, GenBank # NC_000016.10) of the TSC2 gene were screened in TSC2<sup>NEG</sup>, TSC2<sup>POS</sup> cell lines and in LAM-derived primary cells as previously reported<sup>29</sup>. Genomic DNA was extracted with Quick DNA Universal Kit (Zymo Research) and amplified with the “KAPA HIFI Hotstart kit” (Kapa Biosystem). PCR products were purified with the “Zymoclean Gel DNA recovery” kit (Zymo Research), sequenced (GATC, France) and analyzed with the Vector NTI software (Invitrogen).

**Statistical analysis.** Normalized gene expression was done with REST 2009 software (Qiagen). Statistical analysis was done with SPSS Statistics (Version 21, International Business Machines, USA) and PRISM software (GraphPad). Quantitative variables were compared with Mann-Whitney tests. All tests were done with a significance threshold of α = 0.05.

## Results

**Elevated levels of circulating EDN1 and VEGF-D in LAM patients.** EDN1 and VEGF-D levels were measured in sera obtained from 28 LAM patients and 12 female controls with unrelated disease. The circulating level of VEGF-D was significantly higher in LAM patients compared to controls (Fig. 1) with a VEGF-D level above 800 pg/ml in 82% (23/28) of LAM patients (and over 4000 pg/ml VEGF-D for four patients). The level of circulating EDN1 was significantly higher in LAM patients when compared to controls, suggesting a stimulation of the EDN1 production in LAM patients.
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Circulating EDN1 and VEGF-D in LAM patients. Concentrations (pg/ml) of EDN1 and VEGF-D were measured by ELISA in control (n = 12) and LAM (n = 28) sera. P values calculated with the Mann-Whitney test; **p** < 0.01; ****p** < 0.0001.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Characterization of lung-derived primary LAM cells, TSC2<sup>NEG</sup> and TSC2<sup>POS</sup> cells. (A) Morphology by phase contrast microscopy and expression of α-SMA, vimentin and desmin proteins (in green) in LAM cells derived from patient #2634, and in TSC2<sup>NEG</sup> and TSC2<sup>POS</sup> cells. Blue: nucleus stained with DAPI, magnification x400. (B) LAM primary cells and TSC2 cell lines expressed mRNAs of melanocytic (PMEL, TYRP1, MLANA) and smooth-muscle (ACTA2) genes. PPIA was used as reference gene. The figure was constructed from individual gels of PMEL, TYRP1, MLANA, ACTA2 and PPIA PCR products (see supplements). (C) TSC2 point mutation in lung-derived primary LAM from patient #2634.</figcaption></figure>

**Characterization of LAM-derived cells.** We generated lung-derived primary cells from five LAM patients. We observed mixed cell populations in the cultures, with mostly spindle-shaped cells and rare epithelioid-like cells (Fig. 2A), their respective proportion being patient-dependent. Primary LAM cells expressed α-SMA, vimentin and rarely desmin (Fig. 2A) and they were hardly positive for the melanocytic markers HMB45 and PNL2 (not shown). LAM cells expressed PMEL (pre-melanosome protein) but not MLANA (Melan-A) mRNAs (Fig. 2B). All but the cells derived from patient #2634 expressed TYRP1 (Tyrosine related protein 1)
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Dysregulation of expression of genes implicated in the endothelin pathway. (A) Expression of *EDN1*, *EDNRA*, *EDNRB* and *ARRB1* mRNAs has been semi-quantified by RT-qPCR in LAM-derived cells from patients 1300, 1444, 1720, 2634 and 2749 and PASMC (human primary Pulmonary Artery Smooth Muscle Cells) at early passages and in TSC2<sup>NEG</sup> and TSC2<sup>POS</sup> cell lines. The relative mRNA level is reported as mean +/− sem for each cell type. (B) Protein expression of EDNRA and EDNRB in TSC2<sup>NEG</sup> and TSC2<sup>POS</sup> cell lines. Figure was obtained from two individual membranes used to detect ENDRA or EDNRB. Each membrane was incubated with anti β-actin, to estimate the protein load (see supplementals). (C) EDN1, EDNRA and EDNRB gene expression was analyzed after treatment of TSC2<sup>NEG</sup> and TSC2<sup>POS</sup> cells with 100 nM endothelin for 24 h. P values were calculated with Mann-Whitney test.</figcaption></figure>

mRNAs (Fig. 2B). Consistent with the protein expression, ACTA2 (Actin alpha 2) mRNA was expressed in all cells. To summarize, the primary cells derived from LAM patients were consistent with myofibroblast-like cells expressing some of the melanocytic markers.

We morphologically characterized the TSC2<sup>NEG</sup> and TSC2<sup>POS</sup> cell lines; while the TSC2<sup>NEG</sup> cells were spindle-shaped, the TSC2<sup>POS</sup> cells were polygonal (Fig. 2A). TSC2<sup>NEG</sup> and TSC2<sup>POS</sup> cells expressed α-SMA, vimentin and barely desmin (Fig. 2A), but did not express melanocytic proteins as demonstrated by the absence of staining with HMB45 and PNL2 antibodies (not shown). However, both cell lines expressed mRNAs from melanocytic (PMEL and TYRP1) and myofibroblastic (ACTA2) genes (Fig. 2B). TSC2<sup>NEG</sup> and TSC2<sup>POS</sup> cell lines were therefore phenotypically compatible with LAM cells.

Because LAM is often associated with *TSC* mutations, pulmonary LAM cells were screened for mutations in DNA regions corresponding to exons 4041 and intron 5-exon 6 splicing region. We found two silent mutations (as compared to *TSC2* reference gene # NC_000016.10) in LAM patient #2634 with an A/C (nt 46399) transversion in the intronic region and a T/C (nt 46450) transversion in exon 41 (Fig. 2C).

**Expression of EDN1, EDNRA, EDNRB and β-arrestin (ARRB1) is dysregulated in LAM-derived cells.**
We analyzed mRNA expression of EDN1, of its receptors (EDNRA and EDNRB), and of ARRB1 in lung-derived LAM cells and in TSC2<sup>NEG</sup> (Fig. 3A). EDN1, EDNRA, EDNRB and ARRB1 mRNAs were differentially expressed in the primary LAM cells derived from the five patients of our study. LAM cells derived from patient #2634 expressed high levels of EDN1, EDNRA, EDNRB and ARRB1 mRNAs when compared to primary PASMC (Pulmonary Artery Smooth Muscle Cells) (Fig. 3A). Overexpression of EDN1 mRNAs and of at least one of two EDNRs was found in three (#1444, #2634, #2749) out of the five lung-derived LAM cultures (Fig. 3A). Despite some heterogeneity of EDN1, EDNRA, EDNRB and ARRB1 mRNA levels, we evidenced a transcriptional dysregulation of the EDN1 pathway in primary lung-derived LAM cells.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Proliferation of TSC2<sup>NEG</sup> and TSC2<sup>POS</sup> cells. (A) TSC2<sup>NEG</sup> cells have a proliferative advantage when compared to TSC2<sup>POS</sup> cells (***p < 0.001) (n = 3). Proliferation rate has been measured in triplicates and expressed in µM ATP. (B) TSC2<sup>NEG</sup> and TSC2<sup>POS</sup> cell lines were treated with 10µM bosentan, or 10µM BQ123 or 10 nM rapamycin, alone or in combination. Proliferation rate was measured 24h and 72h after inhibitor treatment and expressed as percentage of the control condition (DMSO or mock). Statistical analysis with Mann-Whitney test, *p* < 0.1, ***p* < 0.001.</figcaption></figure>

To investigate TSC2 involvement in EDN pathway dysregulation, we compared mRNA expression of *EDN1*, *EDNR* and *ARRB1* genes in TSC2<sup>NEG</sup> and TSC2<sup>POS</sup> cell lines (Fig. 3A). The mRNA expression of *EDN1* and *ARRB1* was significantly higher (respectively 3.5 and 6.5-fold) while *EDNRB* mRNA expression was significantly lower (10 fold) in TSC<sup>NEG</sup> as compared to TSC<sup>POS</sup> cells (Fig. 3A); ENDRB protein expression was lower in TSC<sup>NEG</sup> cells as compared to TSC<sup>POS</sup> cells (Fig. 3A). Moreover, while EDNRA mRNAs were not differentially expressed in the two cell lines (Fig. 3A), higher levels of EDNRA proteins were detectable in TSC2<sup>NEG</sup> cells as compared to TSC2<sup>POS</sup> cells (Fig. 3B). Conversely, EDNRB proteins were expressed at a weaker level in TSC<sup>NEG</sup> compared to TSC<sup>POS</sup> cells (Fig. 3B).

Next, we tested the effect of endothelin stimulation by incubating both cell lines with 100 nM EDN1 for 24 hours. Following endothelin stimulation, expression of *EDN1* and *EDNRB* mRNAs significantly increased in TSC2<sup>POS</sup> cells but stayed unchanged in TSC2<sup>NEG</sup> cells (Fig. 3B) suggesting that TSC2<sup>NEG</sup> cells lost their ability to respond to endothelin stimulation.

**Proliferation of TSC2<sup>NEG</sup> cells upon EDNR and mTOR inhibitor treatment.** We first showed that the TSC2<sup>NEG</sup> cells had a significant proliferative advantage over the TSC2<sup>POS</sup> cells (Fig. 4A). To investigate EDN pathway involvement in LAM, we tested the effect of ERA such as bosentan (EDNRA and EDNRB inhibitors) and BQ123 (EDNRA inhibitor), with or without rapamycin (mTOR inhibitor) on cell proliferation. Twenty-four hours of treatment with 10 µM bosentan, 10 µM BQ123 or 10 nM rapamycin alone or in combination did not significantly reduce TSC2<sup>POS</sup> or TSC2<sup>NEG</sup> cell proliferation (Fig. 4B). However, at 72 h post treatment, rapamycin was able to reduce the proliferation of both cell line, in a significant way for TSC2<sup>NEG</sup> cells (Fig. 4B). Cell proliferation was also reduced when both cell lines were treated with a combination of BQ123 and rapamycin (Fig. 4B). Surprisingly, bosentan alone significantly increased TSC2<sup>POS</sup> and TSC2<sup>NEG</sup> proliferation (Fig. 4B). These results demonstrated the efficient inhibition of cell proliferation upon rapamycin treatment and suggested that inhibitor of EDNRA alone such as BQ123 do not alter proliferative rate of TSC2<sup>NEG</sup> cells, while bosentan (targeting both EDNRA and EDNRB) activates cell proliferation.

**Migration of TSC2<sup>NEG</sup> cells treated with inhibitors of EDNR or mTOR.** We first established that TSC2<sup>NEG</sup> cells had a significant migration advantage over TSC2<sup>POS</sup> cells (Fig. 5A). Using fluorescent labeling to follow migration, we showed that while TSC2<sup>POS</sup> cells poorly migrated from the upper to the lower compartment of the insert, TSC2<sup>NEG</sup> cells efficiently migrated as measured by the level of fluorescence signal (Fig. 5A). We then analyzed EDN involvement in migration of TSC2<sup>POS</sup> and TSC2<sup>NEG</sup> cells treated with endothelin receptor antagonists, alone or in combination with rapamycin. The addition of 10 µM bosentan or 10 µM BQ123 did not affect cellular migration of TSC2<sup>NEG</sup> or TSC2<sup>POS</sup> cells (Fig. 5B), suggesting that the EDN1 pathway was not implicated in migration. However, rapamycin alone or in combination with BQ123 significantly reduced migration of TSC2<sup>NEG</sup> cells but not of TSC2<sup>POS</sup> cells (Fig. 5B). These results suggest that unlike the EDN pathway, the mTOR pathway is implicated in LAM cell migration *in vitro*.

**TSC2<sup>NEG</sup> cell growth is anchorage-independent and is reduced by rapamycin in combination to ERA.** When maintained in semi-solid medium, TSC2<sup>NEG</sup> cells grew as spheroids, while TSC2<sup>POS</sup> cells were unable to proliferate in 3D conditions (Fig. 6A), suggesting that TSC2<sup>NEG</sup> cells had acquired a transformed
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

**Figure 5.** Migration properties of TSC2<sup>NEG</sup>. Cell migration has been determined by plating cells in Fluoroblock cell culture inserts, 8 µm pore diameter (Corning) and staining with fluorescent dye Calcein-AM 15 hours post treatment. (A) TSC2<sup>NEG</sup> cells have higher migratory potential than TSC2<sup>POS</sup> cells (***p < 0.0001, N = 12), as measured by the fluorescent cells on the lower side of the insert. (B) Cell migration was analyzed in presence of 10 nM rapamycin, 10 µM bosentan or 10 µM BQ123 alone or in combination (rapamycin + bosentan or rapamycin + BQ123). Rapamycin alone (**p < 0.01) or combined to BQ123 (*p < 0.05) significantly inhibited TSC2<sup>NEG</sup> cell migration. Each condition was tested in four points, and experiments replicated 3 times. Statistical analyses were performed using the Mann Whitney test.

or neoplastic phenotype. Interestingly, treatment with rapamycin alone or in combination with BQ123 or bosentan significantly reduced spheroids size when compared to untreated TSC<sup>NEG</sup> cells (Fig. 6B). Endothelin receptor inhibitors alone (BQ123 or bosentan) were unable to decrease spheroid size (Fig. 6B) while bosentan alone had a reverse effect, favoring the growth of TSC2<sup>NEG</sup> spheroids. In conclusion, while ERAs failed to alter anchorage-independent cell growth, rapamycin prevented expansion of TSC2<sup>NEG</sup> spheroids.

## Discussion

LAM cells are characterized by increased proliferation associated to TSC1 and TSC2 mutations, leading to the activation of the mTOR pathway<sup>30</sup>. We hypothesized that dysregulation of the endothelin pathway could be implicated in LAM. The endothelin/β-arrestin/β-catenin pathway is activated in many cancers such as lung, ovary, prostate or breast cancers. It activates proliferation, migration, invasion and survival of cancer cells and stimulates angiogenesis<sup>23</sup>. To explore our hypothesis, we have successfully derived and characterized primary LAM cells that were compatible with the myofibroblast-like subpopulation of LAM, expressing α-smooth actin, desmin and vimentin. Lung-derived LAM primary cells hardly expressed melanocytic markers recognized by HMB45 and PNL2 antibodies but expressed PMEL and TYRP1 mRNAs. The derivation of pulmonary LAM cells revealed to be challenging, as these cells proliferated slowly and were not easily maintained in culture. Despite their limited lifespan, they may be of interest as an *ex vivo* model of LAM. We studied the expression of key genes of the endothelin pathway as an insight of what occurs in LAM patients. We have analyzed the expression of EDN1 and of its receptors EDNRA and EDNRB in primary LAM cells as well as in a TSC2 mutated cell line derived from an angiomyolipoma. Interestingly, mRNA expression of EDN1, EDNRA and EDNRB genes was dysregulated. The TSC2<sup>NEG</sup> cell line and most primary LAM cells expressed high levels of EDN1 and ARRB1 mRNAs and a low level of EDNRB mRNAs when compared to control cells. Unlike TSC2<sup>POS</sup> cells, endothelin stimulation of TSC2<sup>NEG</sup> cells did not affect EDN1, EDNRA or EDNRB mRNA expression, suggesting the loss of positive feedback, possibly due to an independent activation of the EDN1 pathway in these cells.

To investigate the role of these dysregulations in LAM, we studied the effect of EDN receptor inhibitors, combined or not with rapamycin in the LAM cell lines. Proliferation, migration and anchorage-independent growth of TSC2<sup>NEG</sup> were not modified in presence of BQ123, an inhibitor of EDNRA; surprisingly bosentan, which targets EDNRA and EDNRB, activated cell proliferation. However, we showed that rapamycin was an inhibitor of TSC2<sup>NEG</sup> cell migration and anchorage-independent growth.

Despite the expression of EDN in LAM cells, we showed that the EDN1 pathway was not implicated in cellular proliferation or migration or anchorage-independent growth, as demonstrated by the lack of inhibitory effects of ERAs such as BQ123. These results are in contrast with the beneficial effects of ERAs with anti-proliferative action on arterial smooth muscle cells in other diseases such as pulmonary arterial hypertension, where EDN1 is expressed in the pulmonary vasculature<sup>31–33</sup>.
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Anchorage-independent growth of TSC2<sup>NEG</sup> cells. TSC2<sup>POS</sup> and TSC2<sup>NEG</sup> cells were cultured in semi-solid medium for 6 days and monitored for the formation of spheroids. (A) Microscopic observation of cells growing in semi-solid medium. TSC2<sup>NEG</sup> formed multiple spheroids while TSC2<sup>POS</sup> were not able to grow without anchorage. (B) Cells were treated for 24 hours with 10 nM rapamycin, 10 µM BQ123, 10 µM bosentan alone or in dual combination and transferred in semi-solid medium. Spheroids formation was monitored at day 6. Spheroid size has been determined and compared between cells maintained in presence of inhibitors or in non-treated semi-solid medium (BQ123, rapamycin or rapamycin + BQ123) or semi-solid medium containing DMSO (bosentan or rapamycin + bosentan). Statistical analyses were performed using the Mann Whitney test (***p < 0.001).</figcaption></figure>

We have shown that TSC2<sup>NEG</sup> were highly proliferating and migrating cells. It has been reported that TSC2 controls cell migration through its N-terminus and that RhoA activity is necessary<sup>34</sup>. TSC2<sup>POS</sup> control cells used in this study have been generated by re-establishment of normal TSC2 gene expression in the TSC2<sup>NEG</sup> cell line *via* a viral vector. The TSC2<sup>POS</sup> cells that have regained TSC2 expression are characterized by lower proliferation rate, no migration and no anchorage-independent growth; the restoration of TSC2 expression in patients carrying mutations could open to new therapeutic approaches.

Randomized controlled trials have shown that mTOR inhibitors such as sirolimus and everolimus stabilize lung function and reduce angiomyolipoma size in LAM patients<sup>15,35</sup>. Our study is the first report on the inhibitory effect of rapamycin on cellular migration of TSC2-deficient cells. Apart from its effect on cell proliferation, mTOR inhibition may control cell migration in LAM patients; this could participate to the beneficial effect of mTOR inhibitors in LAM patients.

To summarize our results, we have established that the EDN1 pathway is dysregulated in TSC2-deficient cell lines as well as lung-derived primary LAM cells associated with high levels of *EDN* and *ARRB1* mRNA. While targeting EDN receptors with existing antagonists did not prevent proliferation, migration or anchorage-independent growth of LAM cells, β-arrestin 1 could be considered as an alternative target in LAM patients. β-arrestin 1 is a scaffold protein interacting with G protein-coupled receptors and is a major effector of the EDN1 pathway. The binding of EDN1 to EDNRA allows the activation and recruitment of β-arrestin 1 and the formation of β-arrestin 1 and β-catenin complexes<sup>36</sup> that translocate to the nucleus and activate *EDN1* transcription<sup>21,37-39</sup>. Barbadin, a β-arrestin inhibitor, has been recently reported<sup>40</sup> and could help to decipher the role of ARRB1 in LAM as a potential target for new therapeutic tools.

Interestingly, we have shown that rapamycin has an inhibitory effect on LAM cell migration and anchorage-independent growth; this may be of importance in understanding the beneficial effect of mTOR inhibitors in LAM patients.
References

1. Ryu, J. H. et al. The NHLBI lymphangioleiomyomatosis registry: characteristics of 230 patients at enrollment. *Am J Respir Crit Care Med* **173**, 105–111, https://doi.org/10.1164/rccm.200409-1298OC (2006).

2. Moss, J. et al. Prevalence and clinical characteristics of lymphangioleiomyomatosis (LAM) in patients with tuberous sclerosis complex. *Am J Respir Crit Care Med* **164**, 669–671, https://doi.org/10.1164/ajrcmm.164.4.2101154 (2001).

3. Carsillo, T., Astrinidis, A. & Henske, E. P. Mutations in the tuberous sclerosis complex gene TSC2 are a cause of sporadic pulmonary lymphangioleiomyomatosis. *Proc Natl Acad Sci USA* **97**, 6085–6090 (2000).

4. Goncharova, E. A. et al. Abnormal growth of smooth muscle-like cells in lymphangioleiomyomatosis: Role for tumor suppressor TSC2. *Am J Respir Cell Mol Biol* **34**, 561–572, https://doi.org/10.1165/rccm.2005-0300OC (2006).

5. Zhe, X. & Schuger, L. Combined smooth muscle and melanocytic differentiation in lymphangioleiomyomatosis. *J Histochem Cytochem* **52**, 1537–1542, https://doi.org/10.1369/jhc.4A6438.2004 (2004).

6. Matsumoto, Y. et al. Markers of cell proliferation and expression of melanosomal antigen in lymphangioleiomyomatosis. *Am J Respir Cell Mol Biol* **21**, 327–336, https://doi.org/10.1165/ajrcmb.21.3.3693 (1999).

7. Glasgow, C. G., Avila, N. A., Lin, J. P., Stylianou, M. P. & Moss, J. Serum vascular endothelial growth factor-D levels in patients with lymphangioleiomyomatosis reflect lymphatic involvement. *Chest* **135**, 1293–1300, https://doi.org/10.1378/chest.08-1160 (2009).

8. Young, L. R. et al. Serum vascular endothelial growth factor-D prospectively distinguishes lymphangioleiomyomatosis from other diseases. *Chest* **138**, 674–681, https://doi.org/10.1378/chest.10-0573 (2010).

9. McCormack, F. X., Travis, W. D., Colby, T. V., Henske, E. P. & Moss, J. Lymphangioleiomyomatosis: calling it what it is: a low-grade, destructive, metastasizing neoplasm. *Am J Respir Crit Care Med* **186**, 1210–1212, https://doi.org/10.1164/rccm.201205-0848OE (2012).

10. Crooks, D. M. et al. Molecular and genetic analysis of disseminated neoplastic cells in lymphangioleiomyomatosis. *Proc Natl Acad Sci USA* **101**, 17462–17467, https://doi.org/10.1073/pnas.0407971101 (2004).

11. Kumasaka, T. et al. Lymphangiogenesis-mediated shedding of LAM cell clusters as a mechanism for dissemination in lymphangioleiomyomatosis. *Am J Surg Pathol* **29**, 1356–1366 (2005).

12. Krymskaya, V. P. & Shipley, J. M. Lymphangioleiomyomatosis: a complex tale of serum response factor-mediated tissue inhibitor of metalloproteinase-3 regulation. *Am J Respir Cell Mol Biol* **28**, 546–550, https://doi.org/10.1165/rccm.F267 (2003).

13. Zhe, X., Yang, Y., Jakkaraju, S. & Schuger, L. Tissue inhibitor of metalloproteinase-3 downregulation in lymphangioleiomyomatosis: potential consequence of abnormal serum response factor expression. *Am J Respir Cell Mol Biol* **28**, 504–511, https://doi.org/10.1165/rccmb.2002-0124OC (2003).

14. Taveira-DaSilva, A. M., Jones, A. M., Julien-Williams, P., Stylianou, M. & Moss, J. Long-Term Effect of Sirolimus on Serum Vascular Endothelial Growth Factor D Levels in Patients With Lymphangioleiomyomatosis. *Chest*, https://doi.org/10.1016/j.chest.2017.05.012 (2017).

15. McCormack, F. X. et al. Efficacy and safety of sirolimus in lymphangioleiomyomatosis. *N Engl J Med* **364**, 1595–1606, https://doi.org/10.1056/NEJMoa1100391 (2011).

16. Spinella, F. et al. The interplay between hypoxia, endothelial and melanoma cells regulates vascularization and cell motility through endothelin-1 and vascular endothelial growth factor. *Carcinogenesis* **35**, 840–848, https://doi.org/10.1093/carcin/bgu018 (2014).

17. Wulfing, P. et al. Endothelin-1-, endothelin-A-, and endothelin-B-receptor expression is correlated with vascular endothelial growth factor expression and angiogenesis in breast cancer. *Clin Cancer Res* **10**, 2393–2400 (2004).

18. Nelson, J., Bagnato, A., Battistini, B. & Nisen, P. The endothelin axis: emerging role in cancer. *Nat Rev Cancer* **3**, 110–116, https://doi.org/10.1038/nrc990 (2003).

19. Bagnato, A., Spinella, F. & Rosano, L. The endothelin axis in cancer: the promise and the challenges of molecularly targeted therapy. *Can J Physiol Pharmacol* **86**, 473–484, https://doi.org/10.1139/Y08-058 (2008).

20. Rosano, L. & Bagnato, A. Endothelin therapeutics in cancer: Where are we? *Am J Physiol Regul Integr Comp Physiol* **310**, R469–R475, https://doi.org/10.1152/ajpregu.00532.2015 (2016).

21. Rosano, L., Spinella, F. & Bagnato, A. Endothelin 1 in cancer: biological implications and therapeutic opportunities. *Nat Rev Cancer* **13**, 637–651, https://doi.org/10.1038/nrc3546 (2013).

22. Kuntz, M., Leiva-Juarez, M. M. & Luthra, S. Systematic Review of Randomized Controlled Trials of Endothelin Receptor Antagonists for Pulmonary Arterial Hypertension. *Lung* **194**, 723–732, https://doi.org/10.1007/s00408-016-9928-6 (2016).

23. Rosano, L. & Bagnato, A. beta-arrestin1 at the cross-road of endothelin-1 signaling in cancer. *J Exp Clin Cancer Res* **35**, 121, https://doi.org/10.1186/s13046-016-0401-4 (2016).

24. Kim, S. J. et al. Macitentan, a Dual Endothelin Receptor Antagonist, in Combination with Temozolomide Leads to Glioblastoma Regression and Long-term Survival in Mice. *Clin Cancer Res* **21**, 4630–4641, https://doi.org/10.1158/1078-0432.CCR-14-3195 (2015).

25. Smith, M. P. et al. Targeting endothelin receptor signalling overcomes heterogeneity driven therapy failure. *EMBO molecular medicine* **9**, 1011–1029, https://doi.org/10.15252/emmm.201607156 (2017).

26. Cianfroco, R. et al. Blocking endothelin-1-receptor/beta-catenin circuit sensitizes to chemotherapy in colorectal cancer. *Cell Death Differ* **24**, 1811–1820, https://doi.org/10.1038/cdd.2017.121 (2017).

27. Yu, J., Astrinidis, A., Howard, S. & Henske, E. P. Estradiol and tamoxifen stimulate LAM-associated angiomyolipoma cell growth and activate both genomic and nongenomic signaling pathways. *Am J Physiol Lung Cell Mol Physiol* **286**, L694–700, https://doi.org/10.1152/ajplung.00204.2003 (2004).

28. Pfaffl, M. W. A new mathematical model for relative quantification in real-time RT-PCR. *Nucleic Acids Res* **29**, e45 (2001).

29. Cai, X. et al. Phenotypic characterization of disseminated cells with TSC2 loss of heterozygosity in patients with lymphangioleiomyomatosis. *Am J Respir Crit Care Med* **182**, 1410–1418, https://doi.org/10.1164/rccm.201003-0489OC (2010).

30. Goncharova, E. A. et al. Tuberin regulates p70 S6 kinase activation and ribosomal protein S6 phosphorylation. A role for the TSC2 tumor suppressor gene in pulmonary lymphangioleiomyomatosis (LAM). *J Biol Chem* **277**, 30958–30967, https://doi.org/10.1074/jbc.M202678200 (2002).

31. Shao, D., Park, J. E. & Wort, S. J. The role of endothelin-1 in the pathogenesis of pulmonary arterial hypertension. *Pharmacol Res* **63**, 504–511, https://doi.org/10.1016/j.phrs.2011.03.003 (2011).

32. Rubin, L. J. et al. Bosentan therapy for pulmonary arterial hypertension. *N Engl J Med* **346**, 896–903, https://doi.org/10.1056/NEJMoa012212 (2002).

33. Liu, C., Chen, J., Gao, Y., Deng, B. & Liu, K. Endothelin receptor antagonists for pulmonary arterial hypertension. *Cochrane Database Syst Rev*, CD004434, https://doi.org/10.1002/14651858.CD004434.pub5 (2013).

34. Goncharova, E. A., Goncharov, D. A., Lim, P. N., Noonan, D. & Krymskaya, V. P. Modulation of cell migration and invasiveness by tumor suppressor TSC2 in lymphangioleiomyomatosis. *Am J Respir Cell Mol Biol* **34**, 473–480 (2006).

35. Bissler, J. J. et al. Everolimus for angiomyolipoma associated with tuberous sclerosis complex or sporadic lymphangioleiomyomatosis (EXIST-2): a multicentre, randomised, double-blind, placebo-controlled trial. *Lancet* **381**, 817–824, https://doi.org/10.1016/S0140-6736(12)61767-X (2013).

36. Rosano, L. et al. Beta-arrestin links endothelin A receptor to beta-catenin signaling to induce ovarian cancer cell invasion and metastasis. *Proc Natl Acad Sci USA* **106**, 2806–2811, https://doi.org/10.1073/pnas.0807158106 (2009).

37. Spinella, F. et al. Endothelin-1 and endothelin-3 promote invasive behavior via hypoxia-inducible factor-1alpha in human melanoma cells. *Cancer Res* **67**, 1725–1734, https://doi.org/10.1158/0008-5472.CAN-06-2606 (2007).
38. Spinella, F., Rosano, L., Di Castro, V., Natali, P. G. & Bagnato, A. Endothelin-1 induces vascular endothelial growth factor by increasing hypoxia-inducible factor-1alpha in ovarian carcinoma cells. *J Biol Chem* **277**, 27850–27855, https://doi.org/10.1074/jbc. M202421200 (2002).

39. Spinella, F. et al. Endothelin-1 stimulates lymphatic endothelial cells and lymphatic vessels to grow and invade. *Cancer Res* **69**, 2669–2676, https://doi.org/10.1158/0008-5472.CAN-08-1879 (2009).

40. Beautrait, A. et al. A new inhibitor of the beta-arrestin/AP2 endocytic complex reveals interplay between GPCR internalization and signalling. *Nat Commun* **8**, 15054, https://doi.org/10.1038/ncomms15054 (2017).

**Acknowledgements**

NC was the recipient of a Master fellowship from Hospices Civils de Lyon. This project was financially supported by the Hospices Civils de Lyon and the "Fondation du souffle". We thank Pr. EP Henske (Division of pulmonary and critical care medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, USA) for providing us with the 621-101 and 621-103 cell lines and the PLATIM core facility of the SFR Biosciences (UMS 3444/US8), M Gomes for her technical help, F Gonnot and M Chanteperdrix for preliminary results. We thank Pr J Bienvenu and L Garnier, Immunology department, Hospices civils de Lyon for the titration of VEGF and END1 in patient sera.

**Author Contributions**

N.C., F.A., A.B.E., C.L., V.C. had substantial contributions to the conception or design of the work; or the acquisition, analysis, or interpretation of data for the work; drafted the work (N.C., F.A., C.L., V.C.) or revising it critically for important intellectual content (N.C., F.A., C.L., V.C.); All the authors approved the final version and agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.

**Additional Information**

**Supplementary information** accompanies this paper at https://doi.org/10.1038/s41598-018-32795-3.

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_10.png" /></figure>

**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 license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018