Next Article in Journal
A Novel Protozoa Parasite-Derived Protein Adjuvant Is Effective in Immunization with Cancer Cells to Activate the Cancer-Specific Protective Immunity and Inhibit the Cancer Growth in a Murine Model of Colorectal Cancer
Previous Article in Journal
CAH3 from Chlamydomonas reinhardtii: Unique Carbonic Anhydrase of the Thylakoid Lumen
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cells
Volume 13
Issue 2
10.3390/cells13020110
cells-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High Stretch Modulates cAMP/ATP Level in Association with Purine Metabolism via miRNA–mRNA Interactions in Cultured Human Airway Smooth Muscle Cells

by
Mingzhi Luo
*,†,
Chunhong Wang
†,
Jia Guo
,
Kang Wen
,
Chongxin Yang
,
Kai Ni
,
Lei Liu
,
Yan Pan
,
Jingjing Li
and
Linhong Deng
*
Changzhou Key Laboratory of Respiratory Medical Engineering, Institute of Biomedical Engineering and Health Sciences, School of Medical and Health Engineering, Changzhou University, Changzhou 213164, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Cells 2024, 13(2), 110; https://doi.org/10.3390/cells13020110
Submission received: 25 October 2023 / Revised: 6 December 2023 / Accepted: 14 December 2023 / Published: 5 January 2024
Browse Figures
Review Reports Versions Notes

Abstract

:
High stretch (>10% strain) of airway smooth muscle cells (ASMCs) due to mechanical ventilation (MV) is postulated to contribute to ventilator-induced lung injury (VILI), but the underlying mechanisms remain largely unknown. We hypothesized that ASMCs may respond to high stretch via regulatory miRNA–mRNA interactions, and thus we aimed to identify high stretch-responsive cellular events and related regulating miRNA–mRNA interactions in cultured human ASMCs with/without high stretch. RNA-Seq analysis of whole genome-wide miRNAs revealed 12 miRNAs differentially expressed (DE) in response to high stretch (7 up and 5 down, fold change >2), which target 283 DE-mRNAs as identified by a parallel mRNA sequencing and bioinformatics analysis. The KEGG and GO analysis further indicated that purine metabolism was the first enriched event in the cells during high stretch, which was linked to miR-370-5p–PDE4D/AK7. Since PDE4D/AK7 have been previously linked to cAMP/ATP metabolism in lung diseases and now to miR-370-5p in ASMCs, we thus evaluated the effect of high stretch on the cAMP/ATP level inside ASMCs. The results demonstrated that high stretch modulated the cAMP/ATP levels inside ASMCs, which could be largely abolished by miR-370-5p mimics. Together, these findings indicate that miR-370-5p–PDE4D/AK7 mediated high stretch-induced modulation of cAMP and ATP synthesis inside ASMCs. Furthermore, such interactive miRNA–mRNA pairs may provide new insights for the discovery of effective biomarkers/therapeutic targets for the diagnosis and treatment of VILI and other MV-associated respiratory diseases.

Graphical Abstract

1. Introduction

Mechanical ventilation (MV) with a large breathing volume or high pressure is clinically critical for providing life support to patients suffering from acute respiratory distress syndrome (ARDS). However, such unphysiological breathing volume or pressure can over-distend the lung and cause a high stretch (>10% strain) of the airways, which has been postulated to aggregate lung injury, a phenomenon known as ventilator-induced lung injury (VILI) [1,2,3,4,5]. Unfortunately, there is currently no effective pharmacological intervention for VILI, and the consequence of uncontrolled VILI can be high mortality (e.g., up to 60% during the COVID-19 pandemic) [6]. Therefore, there is an urgent need to fully elucidate the pathological mechanisms of VILI in order to improve the outcome of MV treatment [7,8,9,10].
Clinical studies of patients under MV for over six days indicated decreased lung compliance and increased airway inflammation [6,11,12,13,14,15], reflecting pathological alterations in respiratory structure and function [16,17], particularly, aberrant behaviors of the airway smooth muscle cells (ASMCs) [18,19,20,21,22]. Since ASMCs are the primary mechanosensitive cells in the airway system and are known to change phenotype in response to stretch [22], alterations of ASMCs subjected to high stretch during MV are likely to play a central role in the underlying pathological mechanisms and thus may provide potential targets for prevention of VILI.
Indeed, several pharmacological targets of ASMCs have been explored previously, including the blockade of transcription factors and neutralization of inflammation, but all failed in the attempt to effectively reduce the mortality due to VILI [7,8,9,23]. More recently, noncanonical molecular pathways, such as microRNAs (miRNAs)-mRNA interaction networks, have emerged as novel potential targets to mediate the mechanoresponse of ASMCs [24]. For example, miRNA let-7a has been demonstrated to not only modulate mechanically associated genes in melanoma cells [25] but also in biomechanical and inflammatory phenotypes of the ASMCs [26]. It has also been identified that TRPV4 and ERK1/2 are involved in the stretch-induced modulation of inflammatory factors in ASMCs [27,28]. However, these conventional hypothesis-driven studies are insufficient to fully reveal the stretch-induced prominent responses from an integral perspective. This limitation can be overcome with emerging high-throughput screening techniques, such as whole genome-wide RNA sequencing (RNA-Seq) [29,30]. Through bioinformatic analysis of genes derived from RNA-Seq, it is possible to obtain tremendous gene expression information and explore disease-associated genes and their biological functions without the presumptions of any particular hypothesis [31].
Inspired by these developments, we thus sought to systematically identify miRNA and mRNA expression profiles and their specific interactions in ASMCs whether or not subjected to high stretch. In doing so, we hope to reveal the full regulatory network of the high stretch-induced behavioral changes of ASMCs so that a potential new target or approach may be discovered for effective prevention or treatment of VILI [32]. Using RNA sequencing and bioinformatics analysis, we identified 12 miRNAs (7 up-regulated and 5 down-regulated) and 283 targeted mRNAs (141 up-regulated and 142 down-regulated) that were differently expressed in ASMCs cultured in static and high-stretch conditions. However, it was unexpectedly identified that purine metabolism was the top enriched event in the ASMCs exposed to high stretch, which was further linked to 7 miRNA–mRNA interactive pairs, including miR-370-5p–PDE4D/AK7. Since PDE4D/AK7 has been previously linked to cAMP/ATP metabolism in lung diseases [33,34] and now to miR-370-5p in ASMCs, we thus evaluated the effect of high stretch on the cAMP/ATP level inside ASMCs. The results show that high stretch modulated cAMP/ATP levels inside ASMCs, which could be largely abolished by miR-370-5p. Additionally, our study revealed purine metabolism and novel miRNA–mRNA interactions for the response of ASMCs to high stretch, which provides a scientific basis for targeting purine metabolism in treatment/prevention of VILI.

2. Materials and Methods

2.1. Materials

Transferrin (#T8158) and insulin (#91077C) were purchased from Sigma-Aldrich (St. Louis, MO, USA). Collagen type I (#08-115) was purchased from Advanced BioMatrix (#5279, Poway, CA, USA). Dulbecco’s modified Eagle’s medium (#11885092, DMEM), fetal bovine serum (FBS, #16000-044), penicillin-streptomycin (#15140122), and trypsin (#25200056) were purchased from Thermo Fisher Scientific (Waltham, MA, USA). Cell culture flasks (#3073) and plates (#CLS3506) were purchased from Corning Incorporated (Corning, NY, USA). All other reagents were purchased from Fisher Scientific unless noted otherwise.

2.2. Human ASMCs Cultured with/without High Stretch

Primary human ASMCs (#BNCC339826) were purchased from the BeNa Culture Collection (Beijing, China) and used in passages 3–10 which maintained normal “hill and valley” morphology and cellular functions, including physiological response to contractile/relaxant agonist and α-smooth muscle actin (SMA) expression (#ab119952, Abcam, Waltham, MA, USA) [26,35,36,37]. These cells were cultured in DMEM supplemented with 10% FBS, 2 mM L-glutamine, 100 units/mL penicillin, and 100 μg/mL streptomycin in an incubator containing 5% CO2 humidified at 37 °C according to the method described previously. Before the experiment, exponential proliferating ASMCs (2 × 104 cells/cm2) were plated on the type I collagen-coated elastic membrane of Bioflex 6-well plates (FlexCell International, Hillsborough, NC, USA) followed by 24 h serum deprivation with transferrin (5 ng/mL) and insulin (5 ng/mL). Subsequently, a 13% cyclic strain in sinusoidal waveform (0.5 Hz, 1 s deformation alternating with 1 s relaxation) was applied to the ASMCs for 72 h using Flexcell 5000 (FlexCell International), which simulated the in vivo high stretch on ASMCs in ARDS during MV as estimated according to the literature reports [3,4,26,37,38,39]. The ASMCs prepared in the same way but maintained in static condition for 72 h were used as controls. After high-stretch treatment, RNAs were collected from the ASMCs and analyzed for the responsive events and signaling molecules, using whole genome-wide RNA-Seq or quantitative PCR (qPCR), bioinformatics analysis, and function identification. This study protocol is schematically shown in Figure 1.

2.3. Whole Genome-Wide RNA Sequencing of Cultured ASMCs

Total RNA was extracted from cultured ASMCs using TRI Reagent RNA Isolation Reagent (#T9424, Sigma-Aldrich, St. Louis, MO, USA). Then the RNA samples were sent to Shanghai Majorbio Bio-Pharm Technology Co. (Shanghai, China) for quality test, library construction, and sequencing of mRNA and miRNA on the Illumina NovaSeq 6000 and HiSeqTM 2000 platform, respectively (n = 3, per group) [40,41,42].

2.4. Analysis of Differentially Expressed mRNAs and miRNAs in ASMCs

More specifically, the RNA-Seq data was analyzed to identify the differentially expressed mRNAs and miRNAs in ASMCs cultured in either static or high-stretch conditions on the Majorbio Cloud Platform (Version 2.0, https://www.majorbio.com/, accessed on 20 February 2022) [42,43]. The mRNA or miRNA abundance was quantified as transcripts per million (TPM), and those of low abundance (TPM < 10) were filtered out. Then the differential expression of mRNAs or miRNAs due to high stretch was determined using DESeq2, with p values calculated using the Benjamini–Hochberg model. A 2-fold change was used as the cut-off value (p < 0.05) to identify the mRNAs and miRNAs that were differentially expressed (labeled as DE-mRNAs or DE-miRNAs, respectively) in ASMCs when cultured with high stretch as compared to those in static condition [44].
The DE-mRNAs that are target mRNAs of the DE-miRNAs in ASMCs treated with high stretch were further predicted by miRanda (Version 3.3a, http://www.miranda.org/, accessed on 25 September 2022), TargetScan (Version 7.0, http://www.targetscan.org/, accessed on 25 September 2022), and RNAhybrid (Version 2.1.2, http://bibiserv.techfak.uni-bielefeld.de/rnahybrid/, accessed on 25 September2022). Those genes overlapping the target mRNAs of down/up-regulated DE-miRNAs with up/down-regulated DE-mRNAs were identified as target DE-mRNAs.

2.5. Bioinformatics Analysis

Once the target DE-mRNAs were identified, the DE-mRNAs-associated gene functions that were enriched in response to high stretch were further evaluated with Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and Gene Ontology (GO) analysis, which included biological processes (BP) and molecular functions (MF) performed by the Database for Annotation, Visualization, and Integrated Discovery database (DAVID, Version 6.8, https://david.ncifcrf.gov/, accessed on 12 October 2022), respectively, and false discovery rate (FDR) < 0.05 was considered significant [45].

2.6. Analysis of Protein–Protein Interaction (PPI)

The protein–protein interaction (PPI) between DE-mRNAs was constructed using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING 11.5, https://string-db.org/, accessed on 25 October 2022) with a confidence score ≥ 0.70 [46]. A visual PPI network was constructed using Cytoscape software (Version 3.9.1, https://cytoscape.org/, accessed on 25 October 2022). The modules in the PPI network were analyzed using the MCODE plugin (Version 1.6.1, https://apps.cytoscape.org/apps/mcode, accessed on 25 October 2022) in Cytoscape software [47].

2.7. Transfection of Cultured ASMCs with miRNA Mimics or Inhibitor

ASMCs were transfected with the mimics or inhibitor of specific DE-miRNA, in this case miR-370-5p, as it was identified as a primary DE-miRNA involved in high stretch-induced response in ASMCs. The ASMCs were first seeded in 24-well plates and cultured in Dulbecco’s modified Eagle’s medium with 10% fetal bovine serum at 37 °C under 5% CO2 atmosphere. When cultured to ~50% confluence, the cells were transfected using Lipofectamine 3000 (#L3000015, Invitrogen, Carlsbad, CA, USA) with either 50 nM of miR-370-5p mimics (small, chemically modified double-stranded RNAs that mimic endogenous miR-370-5p) or scrambled negative control (NC), and either 50 nM of miR-370-5p inhibitor (chemically modified single-stranded RNA molecules designed to specifically bind to and inhibit endogenous miR-370-5p) or a scrambled NC, respectively, for 48 h according to the manufacturer’s instructions.

2.8. Assessment of ATP/cAMP Level in Cultured ASMCs

The ATP level in cultured ASMCs was determined using an ATP determination kit (#S0026, Beyotime Biotechnology, Shanghai, China) [48], following the protocol provided by the manufacturer. In brief, after high-stretch treatment (0.5 Hz, 13%, 72 h) or miR-370-5p mimics/inhibitor transfection (50 nM, 48 h) as described earlier, the cells were first lysed using lysate. The lysis was then centrifuged at 12,000× g for 5 min (4 °C), and the supernatant was collected for measurement of intracellular ATP. Briefly, a 10 μL sample of cell culture medium or cell lysis suspension was added to 90 μL reaction solution in a 96-well white bottom assay plate followed by incubation for 15 min at room temperature. Luminescence was monitored by a microplate reader (Tecan Spark, TM10M, Meilen, Switzerland) at 563 nm at room temperature.
The cAMP level in cultured ASMCs was measured using an enzyme-linked immunosorbent assay (ELISA). In brief, after static and high-stretch treatment (0.5 Hz, 13%, 72 h), or miR-370-5p mimics/inhibitor transfection (50 nM, 48 h) as described earlier, the cells were subjected to ultrasonic lysis in the presence of the phosphodiesterase inhibitor IBMX (10−5 M) in PBS and centrifuged at 14,000× g for 10 min under 4 °C. Then the cAMP enzyme-linked immunosorbent-assay kit (#E-EL-0056c, Elabscience, Houston, TX, USA) was employed following the manufacturer’s protocol. Absorbance at 450 nm was measured using a colorimetric 96-well plate reader (Infinite F50, Tecan, Zürich, Switzerland). Data are expressed as ratios of the optical density of groups treated with DE-miRNA mimics/inhibitor to that of the vehicle control groups.

2.9. Quantitative PCR Analysis of mRNA Expression in ASMCs Transfected with miRNA Mimics/Inhibitor

Quantitative PCR (qPCR) analysis was used to further evaluate the effect of up/down-regulation of specific DE-miRNA (in this case miR-370-5p) on the mRNA expression level of genes related to signaling that was most enriched in response to high stretch in cultured ASMCs (in this case AK7 and PDE4D in purine metabolism). Primers associated with these genes were obtained from General Biosystems (Chuzhou, China), as shown in Table S1. Briefly, total RNA was purified using TRI Reagent RNA Isolation Reagent (#T9424, Sigma-Aldrich, St. Louis, MO, USA), and the extracted RNA was quantified using a Nanodrop 2000 Spectrophotometer (Thermo Scientific, Willmington, DE, USA). And then total RNA (500 ng) was used to generate 1st strand cDNA using the Revert Aid First Strand cDNA Synthesis Kit (#K1622, Thermo Scientific, Willmington, DE, USA). The qPCR experiment was performed with PowerUp SYBR Green Master Mix (#A25742, Applied Biosystems, Foster City, CA, USA) using the StepOne real-time PCR system (Applied Biosystems) at 50 °C for 2 min, 95 °C for 2 min, followed by 40 cycles of 95 °C for 15 s, 55 °C for 15 s, and 72 °C for 60 s. The reaction system (10 µL) contained 1 µL of cDNA in triplicates according to the manufacturer’s instructions. Calibration and normalization were performed using the 2−∆∆CT method, where ∆CT = CT (target gene) − CT (reference gene, GAPDH) and ∆∆CT = ∆CT (experiment groups) − ∆CT (control groups). Fold changes in mRNA expression of different genes were calculated as the ratio of experiment groups to the control groups from the resulting 2−∆∆CT values from three independent experiments.

2.10. Statistical Analysis

Statistical analysis was performed by using GraphPad Prism 8.0 (Graph Pad Software, San Diego, CA, USA). Unless stated otherwise, each sample was tested in triplicate, and each experiment for a given study condition was repeated a number of times depending on group size (n). Data were reported as means ± S.E.M. Comparisons of means between 2 groups were performed using an unpaired Student t-test. Comparisons of means among three or more groups were performed by one-way analysis of variance (ANOVA) followed by a Post Hoc test (Turkey HSD method) or two-way ANOVA followed by the Tukey HSD method. Differences were considered significant at p < 0.05.

3. Results

3.1. High Stretch-Induced Differentially Expressed mRNAs and miRNAs in ASMCs

The mRNA response of ASMCs to high stretch was evaluated with mRNA sequencing (GSE206435). We defined the differentially expressed mRNAs (DE-mRNAs) in experimental groups as mRNAs with count ≥ 10 that had an absolute log2 (fold change) expression ≥1 together with a p < 0.05 when compared with their counterparts in static cells. We found a total of 2858 DE-mRNAs, among which 1466 were downregulated and 1392 were upregulated.
The miRNA response of ASMCs to high stretch was evaluated with miRNA sequencing (GSE206435). Similar to mRNAs, we defined the differentially expressed miRNAs (DE-miRNAs) in experimental groups as miRNAs with a count ≥ 10 that had an absolute log2 (fold change) expression ≥1 together with a p < 0.05 when compared with their counterparts in static cells. We found, in total, only 12 DE-miRNAs, among which 7 were upregulated, and 5 were downregulated, as shown in Table 1.
These 12 miRNAs were predicted to target 4464 genes in total, among which 1942 were targeted by the 7 upregulated DE-miRNAs, and 2522 were targeted by the 5 downregulated DE-miRNAs. We further compared the target genes of down/up-regulated DE-miRNAs with those of up/down-regulated DE-mRNAs and found 283 genes that overlapped these two groups (142 overlapped the target genes of down-regulated DE-miRNAs and the genes of up-regulated mRNAs, and 141 overlapped the target genes of up-regulated miRNAs and the genes of down-regulated mRNAs, as shown in Figure 2A and Table S2). We identified these overlapping genes as target DE-mRNAs, whose expressions were visualized in the volcano plot as shown in Figure 2B. It can be seen that hypoxia up-regulated 1 (HYOU1), activating transcription factor 3 (ATF3), solute carrier family 3 member 2 (SLC3A2), SEL1L adaptor subunit of SYVN1 ubiquitin ligase (SEL1L), fibronectin leucine rich transmembrane protein 2 (FLRT2), KLF transcription factor 15 (KLF15), ATPase plasma membrane Ca2+ transporting 1 (ATP2B1), lysosomal thiol reductase (IFI30), MIA SH3 domain ER export factor 2 (MIA2), and MAPK interacting serine/threonine kinase 2 (MKNK2) were the top 10 upregulated target DE-mRNAs, and LBH regulator of WNT signaling pathway (LBH), neuronal regeneration-related protein (NREP), thrombospondin 1 (THBS1), catalase (CAT), guanylate cyclase 1 soluble subunit beta 1 (GUCY1B1), peripheral myelin protein 22 (PMP22), phosphodiesterase 7B (PDE7B), RAB7B member RAS oncogene family (RAB7B), stathmin 1 (STMN1), and elastin (ELN) were the top 10 downregulated target DE-mRNAs in ASMCs cultured with high stretch.

3.2. High Stretch-Induced Enrichment of Cellular Functions and Signaling Pathways in Cultured ASMCs

The biological function of these 283 target DE-mRNAs was analyzed by KEGG pathway and GO enrichment analysis using DAVID. As shown in Figure 3A and Table S3, KEGG pathway analysis indicated that purine metabolism-related signaling, mechanical stress/strain-related signaling, such as cGMP-PKG signaling pathway, MAPK signaling pathway, fluid shear stress and atherosclerosis, amoebiasis, and Ras signaling were all enriched to various extents in ASMCs when cultured with high stretch. However, the top enriched signaling in ASMCs when cultured with high stretch appeared to be related to purine metabolism, which involved 9 target DE-mRNAs, including phosphoribosyl pyrophosphate synthetase 1 (PRPS1), GUCY1B1, phosphodiesterase 4D (PDE4D), ectonucleotide pyrophosphatase/phosphodiesterase 1 (ENPP1), adenylate kinase 4 (AK4), phosphodiesterase 5A (PDE5A), PDE7B, adenylate kinase 4 (AK7), NME/NM23 nucleoside diphosphate kinase 1 (NME1).
The results of GO analysis (Figure 3B and Table S4) also show that purine metabolism-related biological processes, including regulation of GTPase activity (Rho GTPase activating protein 9 (ARHGAP9), FYVE RhoGEF and PH domain containing 6 (FGD6), RAS p21 protein activator 4B (RASA4B), synaptic Ras GTPase activating protein 1 (SYNGAP1), RAS p21 protein activator 4 (RASA4), chimerin 1 (CHN1), Cbl proto-oncogene B (CBLB)) and ATP metabolic process (ATP binding cassette subfamily C member 6 (ABCC6), ATPase H+ transporting V1 subunit B2 (ATP6V1B2), ENPP1, AK4), as well as purine metabolism-related molecular functions including 3’,5’-cyclic-AMP phosphodiesterase activity (PDE4D, ENPP1, PDE7B), nucleoside diphosphate kinase activity (AK4, AK7, NME1), and 3’,5’-cyclic-nucleotide phosphodiesterase activity (PDE4D, PDE5A, PDE7B), were also enriched. Other events, such as mechanical (basement membrane organization and cell migration) and inflammatory (positive regulation of interleukin-6 (IL-6), tumor necrosis factor-alpha (TNFα), monocyte chemoattractant protein-1 (MCP-1), and bone morphogenetic protein 1 (BMP) production)-related events, were also enriched in ASMCs treated with high stretch.

3.3. High Stretch-Induced Enrichments of PPI in Cultured ASMCs

The protein–protein interaction (PPI) network of DE-mRNAs contained 93 nodes and 115 interactions that clustered into three modules, as shown in Figure 4A. The KEGG analysis (Figure 4B and Table S5) indicates that the 38 target DE-mRNAs in the first cluster of the PPI network were enriched due to high stretch in various functions of the cultured ASMCs related to cellular metabolism, including, in rank-order, purine metabolism (PDE4D, ENPP1, AK4, PDE5A, PDE7B, AK7, NME1), nucleotide metabolism, tryptophan metabolism, metabolism pathways, biosynthesis of cofactors, and thiamine metabolism. The GO analysis (Figure 4C) further indicates that the biological processes and molecular functions enriched in cultured ASMCs due to high stretch include, in rank order, nucleotide diphosphate phosphorylation (AK4, AK7, NME1), ATP metabolic process (ABCC6, ENPP1, AK4), cAMP catabolic process (PDE4D, PDE7B), nucleoside triphosphate biosynthetic process (AK4, AK7), and nucleoside monophosphate phosphorylation (AK4, AK7), 3′,5′-cyclic-cAMP phosphodiesterase activity (PDE4D, ENPP1, PDE7B), nucleoside diphosphate kinase activity (AK4, AK7, NME1), 3′,5′-cyclic-nucleotide phosphodiesterase activity (PDE4D, PDE5A, PDE7B), and adenylate kinase activity (AK4, AK7). All these data indicate that the 7 target DE-mRNAs, including PDE4D, ENPP1, AK4, PDE5A, PDE7B, AK7, and NME1, were the genes most responsive to high stretch in ASMCs. While some of them were also involved in other cell functions, these genes were all involved in regulating purine metabolism. More specifically, PDE5A is involved in guanosine metabolism and the other six DE-mRNAs are involved in adenosine metabolism. Since the second messenger cAMP and ATP are crucial signaling and energy molecules and regulate mechanical and inflammatory phenotypes of ASMCs, the high stretch-induced signaling related to purine metabolism, particularly the cAMP and ATP metabolism in ASMCs, might be key to understanding the regulatory roles as well as the therapeutic potential of miRNAs in VILI [49,50].

3.4. High Stretch-Responsive miRNA–mRNA Interactions in Relation to Regulation of cAMP/ATP in Cultured ASMCs

Figure 5A shows seven pairs of possible interactions between the DE-mRNAs targeted by their corresponding DE-miRNAs, namely, PDE4D–miR-370-5p, PDE4D–miR-543, PDE7B–miR-29b-3p, AK4–miR-148a-3p, AK7–miR-370-5p, ENPP1–miR-194-5p, and NME1–miR-146-5p, with respect to their roles in regulating cAMP and ATP metabolism. Since PDE4D and PDE7B both hydrolyze cAMP (cAMP→AMP), the high stretch-increased mRNA expression of PDE4D in association with the downregulated expression of miR-370-5p or miR-543 may result in decreased cAMP concentration in ASMCs, whereas the decreased PDE7B with upregulated miR-29b-3p may result in increased cAMP concentration in ASMCs.
On the other hand, AK4 and AK7 (in mitochondrial and cytosolic, respectively) either efficiently phosphorylates AMP to generate ATP or catalyzes the reversible transfer of the terminal phosphate group between ATP and AMP (AMP→ATP). Therefore, the high stretch-decreased AK4 with upregulated miR-148a-3p may decrease ATP concentration in ASMCs, while the increased AK7 with downregulated miR-370-5p may increase ATP concentration in ASMCs.
In addition, ENPP1 and NME1 either preferentially hydrolyze ATP with the release of pyrophosphate and adenosine polyphosphates (ATP→AMP) or transfer the ATP gamma phosphate to the NDP beta phosphate via a ping-pong mechanism and thus promote the synthesis of nucleoside triphosphates other than ATP (ATP→ADP). Thus, the high stretch-decreased ENPP1 or NME1 with upregulated miR-194-5p or miR-146-5p may all lead to increased ATP concentration in ASMCs.
Figure 5B and 5C show that the intracellular level of cAMP and ATP as measured by ELISA decreased by 20% (p = 0.0248) and increased by 60% (p = 0.0047), respectively, in ASMCs cultured with high stretch as compared to their counterparts under static condition.

3.5. Effect of miR-370-5p on the High Stretch-Induced Variation of cAMP/ATP Concentration in Cultured ASMCs

Among all the above miRNAs, miR-370-5p can target PDE4D and AK7 to regulate the concentration of cAMP and ATP and enrich the purine metabolisms according to the results of KEGG analysis (Figure S1). Therefore, we speculated that miR-370-5p might be the primary miRNA in regulating high stretch-induced variations of cAMP/ATP levels in ASMCs. To verify this speculation, we transferred the cultured ASMCs with miR-370-5p mimics or inhibitor followed by high-stretch treatment. The cell viability of, and transfection efficiency of, miR-370-5p mimics/inhibitor in ASMCs was evaluated with CCK8 assay and qPCR, respectively. The results indicated that the transfection of miR-370-5p mimics (p = 0.0966) or inhibitor (p = 0.7105) did not significantly affect the viability of the cultured ASMCs (Figure 6A), which was also highly efficient as it increased or decreased the expression level of miR-370-5p by ~280% (p = 0.0001) or ~80% (p = 0.0421), respectively (Figure 6B). In addition, the mRNA expression of PDE4D (Figure 6C, p = 0.0230) and AK7 (Figure 6D, p = 0.0139) was decreased by miR-370-5p, respectively. At the same time, the level of cAMP (Figure 6E, p = 0.0331) and ATP (Figure 6F, p = 0.0036) inside ASMCs was increased and decreased by miR-370-5p, respectively. Consistent with these effects, we found that the miR-370-5p mimics and inhibitor respectively decreased (p = 0.0321) and increased (p = 0.0251) the basal tone, but neither significantly affected the relaxation of ASMCs (Figure S2). The basal tone and relaxation of ASMCs were evaluated by measuring cell stiffness using optical magnetic twisting cytometry (OMTC) before (basal tone) and after (relaxation) exposure to isoproterenol (ISO, a conventional β-agonists) as described elsewhere [51,52].
More importantly, the presence of miR-370-5p mimics in cultured ASMCs abolished the high stretch-induced decrease in cAMP level (p = 0.0212), as shown in Figure 6G, and reversed the increase in the ATP level (p = 0.0044), as shown in Figure 6H. These data suggest that the interactions between miR-370-5p and PDE4D, miR-370-5p, and AK7 may be essential in mediating high stretch-induced pathological responses in ASMCs through modulation of the purine metabolism of the cells.

4. Discussion

In this study, we first identified, using RNA-Seq 12, DE-miRNAs and 283 target DE-mRNAs that were highly up/down-regulated in ASMCs in response to high stretch. Further bioinformatic analysis, together with PPI network clustering of these DE-mRNAs, indicate that the high stretch-induced prominent change in the cultured ASMCs occurred in the cells’ purine metabolism. More specifically, the cAMP level decreased but the ATP level increased inside ASMCs due to high stretch. At the same time, miR-370-5p expression was decreased and AK7/PDE4D expression was enhanced by high stretch, and the high stretch-induced change in the cAMP and ATP level inside ASMCs could be partially abolished by miR-370-5p mimics. These data suggest that purine metabolism and related miRNA–mRNA interactions may serve as candidate biomarkers and therapeutic targets for ASMCs in high stretch-related diseases.
Although ASMCs may respond to high stretch by changing the expression, activation, and/or distribution of signaling molecules, the change of expression is relatively easy to evaluate by RNA-Seq. Therefore, with RNA-Seq data, we identified the genes and miRNAs that are differentially expressed in ASMCs when cultured in either high-stretch or static conditions, as a key to understanding the high stretch-induced events and associated regulating miRNA–mRNA interactions in ASMCs. We determined the differentially expressed mRNA and miRNA by analyzing RNA-Seq data with DESeq2, considering its high specificity and true positive rate [29,30]. In addition, we used DAVID, a highly integrated data-mining environment to determine the functional enrichment of the molecules [53]. Based on miRNA–mRNA integrated bioinformatics, we found that purine metabolism was the most responsive event in ASMCs exposed to high stretch.
Purine metabolism has been widely studied for its essential roles in cell behaviors such as cell energy production, proliferation, and survival. In particular, other than serving as building blocks for DNA, purine nucleotides such as ATP are crucial for providing cellular energy, intracellular signaling, and cofactors to promote cell survival and proliferation [49]. Several diseases have also been linked to purine metabolism [54]. In airways, purine is not only involved in the normal energy flow but can also serve as a VILI marker in bronchoalveolar lavage fluid [55,56,57]. For example, lung tissues dissected from rats after exposure to long-term MV show distinct metabolomic profiles indicating enhanced energy production in the lung [58]. Our results of high stretch-induced dramatical responses in purine metabolism, particularly the enhanced ATP level in ASMCs, are consistent with the in vivo finding. When exposed to cyclic high stretch, cells may experience passive exercise and inevitably enhance ATP consumption. ATP could eventually be consumed during these cellular activities, and, hence, the increased probability of ATP synthesis may compensate for the ATP consumption during cyclic high stretch. Therefore, increased ATP may be a feedback mechanism to enhance cellular response to external mechanical stimuli. In addition to functioning as an essential energy molecule, ATP is secreted from various types of lung cells, including ASMCs, when the cells are subjected to mechanical stimulation such as high stretch. The secreted ATP then can mediate the change of basal tone and/or contractility of ASMCs [59]. ATP can also activate purinergic receptors to potently regulate both the physical and chemical behaviors of ASMCs. We did observe that ASMCs exposed to increased levels of ATP became stiffer, more contractile, and pro-inflammatory but not more proliferative. Therefore, the high stretch-enhanced ATP synthesis in ASMCs may be responsible for increased contractility and thus decreased compliance of the airways of mechanically ventilated patients.
Previous studies have shown that miRNAs may change the function of purine metabolism. For example, miR-151a-5p is reported to participate in the regulation of cellular respiration and target cytochrome B for the production of ATP [60]. In this study, we found that ASMCs responded to high stretch with marked changes in interactions between certain miRNAs and their target genes, including miR-370-5p–PDE4D, miR-543–PDE4D, miR-29b-3p–PDE7B, miR-370-5p–AK7, miR-148a-3p–AK4, miR-194-5p–ENPP1, and miR-146-5p–NME1, which may all impact the mitochondrial respiratory function. Among them, miR-370-5p–AK7 is more likely to be involved in ATP synthesis inside ASMCs during high stretch because AK7 is known to function as phosphotransferase in energy homeostasis, and as adenylate kinase to catalyze the production and breakdown of adenine nucleotide in a reversible transfer method [54]. It has also been shown that depletion of AK7 protein will lead to impaired mitochondrial respiratory function, especially the synthesis of ATP [61].
On the other hand, cAMP is a ubiquitous second messenger and its principal function in airways is to mediate relaxation/contraction of ASMCs [50]. For example, the classical bronchodilator, β2-agonist, activates the Gs-coupled β2 adrenergic receptor (β2AR) on ASMCs to activate adenylyl cyclase and, hence, increase the intracellular level of cAMP to mediate relaxation of ASMCs. Acetylcholine released from postganglionic parasympathetic nerves instead inactivates adenylyl cyclase and, hence, lowers the intracellular level of cAMP to mediate contraction of ASMCs. Thus, the intracellular level of cAMP is essential for balancing the contraction and relaxation of ASMCs. In more detail, cAMP is mainly generated by adenylyl cyclase V and VI from cytoplasmic ATP and broken down by intracellular phosphodiesterases. In this study, we found that high stretch decreased intracellular levels of cAMP in association with modulated mRNA expressions of several cAMP-specific 3’,5’-cyclic phosphodiesterases, including PDE4D. Previous studies have widely used optical magnetic twisting cytometry (OMTC) to measure cell stiffness as a surrogate assay in order to assess contraction/relaxation behaviors, such as basal tone and contractility of ASMCs, especially when the cells are cultured with cyclic stretch, exogenous cAMP, and ATP [51,52]. Generally, cyclic stretch and exogenous ATP can enhance basal tone or contractility, but cAMP does the opposite, which has been commonly attributed to the cytoskeletal reorganization and associated Ca2+ signaling in the cells. In this study, however, we found that high stretch modulated intracellular levels of cAMP/ATP via miRNA–mRNA, particularly miR-370-5p–PDE4D/AK7 interactions. In addition, we found that miR-370-5p did impact cell stiffness, but not the cell stiffness reduction induced by β-agonist, suggesting that high stretch through miR-370-5p–PDE4D interaction is a mechanical modulator of bronchial airway tone [62]
In addition to the modulation of airway tone through contraction/relaxation, ASMCs also modulate airway inflammation by secreting or expressing a wide array of immunomodulatory mediators in response to extracellular stimuli such as high stretch [63]. For example, mechanical stress can induce more pro-inflammatory mediators including IL-1β, IL-8, IL-6, and TNF-α in the distal air spaces of the lung [64]. In this study, we found that high stretch led to enriched positive regulation of production for IL-6, TNF-α, MCP-1, and BMP-1. These inflammatory factors are potential therapeutic targets for VILI since they are upstream targets in response to stretch and initiate many pathophysiological processes linked to VILI, including barrier permeability, alveolar fluid clearance, and apoptosis [65,66]. Although ASMCs share with canonical immune cells the same mechanisms to secrete cytokines and chemokines, the unique miRNA–mRNA interactions revealed in this study for regulating inflammatory factors in ASMCs may provide novel therapeutic targets for VILI. The demonstration of high stretch-modulated cAMP and ATP levels and related purine metabolism inside ASMCs via the miR-370-5p–PDE4D/AK7 pathway also warrants the exploration of these effects in other lung cells that are also important in VILI, such as airway epithelial cells and macrophages. For example, cAMP triggers Na+ absorption, anti-inflammation, mucociliary clearance of the airway epithelium, and has been recognized as a potential therapeutic target of cystic fibrosis [34,67]. Therefore, when other types of lung cells were targeted, the miR-370-5–PDE4D/AK7 pathway may have had other pharmacological functions.
To our knowledge, this is the first report of differential responses of miRNAs and mRNAs in ASMCs when exposed to high stretch using a whole genome-wide approach, which surprisingly linked high stretch to the purine metabolism, especially cAMP and ATP synthesis in ASMCs. However, there are limitations to this work. For instance, the purine metabolism modulated by miR-370-5p–PDE4D/AK7 interactions to regulate cAMP/ATP levels inside ASMCs during high stretch was identified solely based on target prediction, statistical evidence, and mRNA evaluation. The effects of high stretch on PDE4D and AK7 in ASMCs are yet to be experimentally verified at the protein level using methods such as dual-luciferase reporter assay and protein evaluation. It is also known that the effect of cyclic stretch on ASMCs is amplitude-dependent. Although previous studies have already shown that high stretch can cause pathological effects on ASMCs, such as enhanced ATP secretion and disrupted alveolar epithelial barrier integrity [59,68], it is unclear whether cyclic stretch at low tidal volume or gradient amplitude would incrementally change ASMCs as they proceed toward pathological remodeling so that more pathogenic signals/pathways may be identified. Furthermore, it remains unknown whether the high-stretch modulation of intracellular levels of cAMP and ATP depends on specific membrane transporters, such as ATP-binding cassettes, and how cAMP and ATP coordinate to regulate downstream signals such as AMP-activated protein kinase activity. Finally, the in vitro findings need to be investigated with in vivo animal models and ultimately human subjects to further confirm their pharmacological and therapeutic implications in relation to the treatment of VILI.

5. Conclusions

The present study revealed that purine metabolism regulated through miRNA–mRNA interactions was a prominent cellular response of ASMCs to high stretch in association with significant down/up-regulation of intracellular cAMP/ATP levels. Furthermore, miR-370-5p for transcription regulation of PDE4D and AK7 was the major regulatory molecule responsible for high stretch-induced changes related to purine metabolism in ASMCs. Taken together, the link established between the physiological functions, such as purine metabolism and related miRNA–mRNA interactions, in ASMCs during high stretch may help discover potential new biomarkers and targets for diagnosis and treatment of respiratory diseases associated with MV, such as VILI.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cells13020110/s1. Figure S1: Purine metabolism was enriched in the KEGG pathway of the 59 target DE-mRNAs of miR-370-5p. The circle indicates the number of enriched genes in each term. Figure S2: Effect of miR-370-5p mimics, inhibitor, and their controls on the cell stiffness of ASMCs before (A) and after (B) exposure to isoproterenol (ISO, a conventional β-agonist). Cell stiffness was measured by optical magnetic twisting cytometry as described elsewhere [34,50]. Each experiment was repeated three times (n = 3). Table S1: Primer sequence of mRNA. Table S2: 283 target DE-mRNAs by 12 DE-miRNAs. Table S3: Significantly enriched KEGG pathway of 283 target DE-mRNAs. Table S4: Significant enrichment of GO analysis of target DE-mRNAs in ASMCs. Table S5: Significant enrichment of KEGG of the first module of PPI.

Author Contributions

Conceptualization, M.L. and L.D.; methodology, C.W., C.Y., Y.P., J.L., L.L. and J.G.; software, K.W.; validation, K.W.; formal analysis, C.W.; investigation, C.W.; resources, C.W.; data curation, J.G.; writing—original draft preparation, M.L.; writing—review and editing, L.D.; visualization, K.N.; supervision, M.L. and L.D.; project administration, M.L. and L.D.; funding acquisition, M.L. and L.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Natural Science Foundation of China (NSFC), Grants No. 12072048 to M.L., 12272063, and 11532003 to L.D., and partially supported by the Science and Technology Innovation Leading Plan of High Tech Industry in Hunan Province, China, Grant No. 2020SK2018 to L.D.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Acknowledgments

We thank Mingxing Ouyang (Changzhou University) for the technical support and Xiang Wang (Changzhou University) for the design of the experiments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sinclair, S.E.; Molthen, R.C.; Haworth, S.T.; Dawson, C.A.; Waters, C.M. Airway strain during mechanical ventilation in an intact animal model. Am. J Respir. Crit. Care Med. 2007, 176, 786–794. [Google Scholar] [CrossRef] [PubMed]
  2. Ibrahim, I.B.M.; Aghasafari, P.; Pidaparti, R.M. Transient mechanical response of lung airway tissue during mechanical ventilation. Bioengineering 2016, 3, 4. [Google Scholar] [CrossRef] [PubMed]
  3. Bagchi, A.; Vidal Melo, M.F. Follow the voxel-a new method for the analysis of regional strain in lung injury. Crit. Care Med. 2018, 46, 1033–1035. [Google Scholar] [CrossRef]
  4. Retamal, J.; Hurtado, D.; Villarroel, N.; Bruhn, A.; Bugedo, G.; Amato, M.B.P.; Costa, E.L.V.; Hedenstierna, G.; Larsson, A.; Borges, J.B. Does regional lung strain correlate with regional inflammation in acute respiratory distress syndrome during nonprotective ventilation? An experimental porcine study. Crit. Care Med. 2018, 46, e591–e599. [Google Scholar] [CrossRef] [PubMed]
  5. Slutsky, A.S. Lung injury caused by mechanical ventilation. Chest 1999, 116, 9S–15S. [Google Scholar] [CrossRef]
  6. Grieco, D.L.; Bongiovanni, F.; Chen, L.; Menga, L.S.; Cutuli, S.L.; Pintaudi, G.; Carelli, S.; Michi, T.; Torrini, F.; Lombardi, G.; et al. Respiratory physiology of COVID-19-induced respiratory failure compared to ARDS of other etiologies. Crit. Care 2020, 24, 529. [Google Scholar] [CrossRef]
  7. Fan, E.; Del Sorbo, L.; Goligher, E.C.; Hodgson, C.L.; Munshi, L.; Walkey, A.J.; Adhikari, N.K.J.; Amato, M.B.P.; Branson, R.; Brower, R.G.; et al. An official american thoracic society/european society of intensive care medicine/society of critical care medicine clinical practice guideline: Mechanical ventilation in adult patients with acute respiratory distress syndrome. Am. J. Respir. Crit. Care Med. 2017, 195, 1253–1263. [Google Scholar] [CrossRef]
  8. Slutsky, A.S.; Ranieri, V.M. Ventilator-induced lung injury. N. Engl. J. Med. 2013, 369, 2126–2136. [Google Scholar] [CrossRef]
  9. Cabrera-Benitez, N.E.; Laffey, J.G.; Parotto, M.; Spieth, P.M.; Villar, J.; Zhang, H.; Slutsky, A.S. Mechanical ventilation-associated lung fibrosis in acute respiratory distress syndrome: A significant contributor to poor outcome. Anesthesiology 2014, 121, 189–198. [Google Scholar] [CrossRef]
  10. Villar, J.; Ferrando, C.; Tusman, G.; Berra, L.; Rodriguez-Suarez, P.; Suarez-Sipmann, F. Unsuccessful and successful clinical trials in acute respiratory distress syndrome: Addressing physiology-based gaps. Front. Physiol. 2021, 12, 774025. [Google Scholar] [CrossRef]
  11. Tsumura, H.; Harris, E.; Brandon, D.; Pan, W.; Vacchiano, C. Review of the mechanisms of ventilator induced lung injury and the principles of intraoperative lung protective ventilation. AANA J. 2021, 89, 227–233. [Google Scholar] [PubMed]
  12. Cronin, J.N.; Camporota, L.; Formenti, F. Mechanical ventilation in COVID-19: A physiological perspective. Exp. Physiol. 2022, 107, 683–693. [Google Scholar] [CrossRef] [PubMed]
  13. Seeley, E.J.; McAuley, D.F.; Eisner, M.; Miletin, M.; Zhuo, H.; Matthay, M.A.; Kallet, R.H. Decreased respiratory system compliance on the sixth day of mechanical ventilation is a predictor of death in patients with established acute lung injury. Respir. Res. 2011, 12, 52. [Google Scholar] [CrossRef] [PubMed]
  14. Cortes-Puentes, G.A.; Keenan, J.C.; Adams, A.B.; Parker, E.D.; Dries, D.J.; Marini, J.J. Impact of chest wall modifications and lung injury on the correspondence between airway and transpulmonary driving pressures. Crit. Care Med. 2015, 43, e287–e295. [Google Scholar] [CrossRef]
  15. Nair, G.B.; Al-Katib, S.; Podolsky, R.; Quinn, T.; Stevens, C.; Castillo, E. Dynamic lung compliance imaging from 4DCT-derived volume change estimation. Phys. Med. Biol. 2021, 66, 21NT06. [Google Scholar] [CrossRef]
  16. Shapiro, M.B.; Bartlett, R.H. Pulmonary compliance and mechanical ventilation. Arch. Surg. 1992, 127, 485–486. [Google Scholar] [CrossRef] [PubMed]
  17. Katira, B.H. Ventilator-induced lung injury: Classic and novel concepts. Respir. Care 2019, 64, 629–637. [Google Scholar] [CrossRef] [PubMed]
  18. Mauri, T.; Lazzeri, M.; Bellani, G.; Zanella, A.; Grasselli, G. Respiratory mechanics to understand ARDS and guide mechanical ventilation. Physiol. Meas. 2017, 38, R280. [Google Scholar] [CrossRef]
  19. Guérin, C.; Cour, M.; Argaud, L. Airway closure and expiratory flow limitation in acute respiratory distress syndrome. Front. Physiol. 2022, 12, 815601. [Google Scholar] [CrossRef]
  20. Asano, S.; Ito, S.; Morosawa, M.; Furuya, K.; Naruse, K.; Sokabe, M.; Yamaguchi, E.; Hasegawa, Y. Cyclic stretch enhances reorientation and differentiation of 3-D culture model of human airway smooth muscle. Biochem. Biophys. Rep. 2018, 16, 32–38. [Google Scholar] [CrossRef] [PubMed]
  21. Schwingshackl, A. The role of stretch-activated ion channels in acute respiratory distress syndrome: Finally a new target? Am. J. Physiol. Lung Cell. Mol. Physiol. 2016, 311, L639–L652. [Google Scholar] [CrossRef]
  22. Bartolak-Suki, E.; LaPrad, A.S.; Harvey, B.C.; Suki, B.; Lutchen, K.R. Tidal stretches differently regulate the contractile and cytoskeletal elements in intact airways. PLoS ONE 2014, 9, e94828. [Google Scholar] [CrossRef]
  23. Uhlig, S.; Uhlig, U. Pharmacological interventions in ventilator-induced lung injury. Trends Pharmacol. Sci. 2004, 25, 592–600. [Google Scholar] [CrossRef] [PubMed]
  24. Maneechotesuwan, K. Role of microRNA in severe asthma. Respir. Investig. 2019, 57, 9–19. [Google Scholar] [CrossRef]
  25. Müller, D.W.; Bosserhoff, A.K. Integrin β3 expression is regulated by let-7a mirna in malignant melanoma. Oncogene 2008, 27, 6698–6706. [Google Scholar] [CrossRef]
  26. Wen, K.; Ni, K.; Guo, J.; Bu, B.; Liu, L.; Pan, Y.; Li, J.; Luo, M.; Deng, L. MircroRNA let-7a-5p in airway smooth muscle cells is most responsive to high stretch in association with cell mechanics modulation. Front. Physiol. 2022, 13, 830406. [Google Scholar] [CrossRef] [PubMed]
  27. Pairet, N.; Mang, S.; Fois, G.; Keck, M.; Kuhnbach, M.; Gindele, J.; Frick, M.; Dietl, P.; Lamb, D.J. TRPV4 inhibition attenuates stretch-induced inflammatory cellular responses and lung barrier dysfunction during mechanical ventilation. PLoS ONE 2018, 13, e0196055. [Google Scholar] [CrossRef] [PubMed]
  28. Fang, X.; Ni, K.; Guo, J.; Li, Y.; Zhou, Y.; Sheng, H.; Bu, B.; Luo, M.; Ouyang, M.; Deng, L. FRET visualization of cyclic stretch-activated erk via calcium channels mechanosensation while not integrin β1 in airway smooth muscle cells. Front. Cell Dev. Biol. 2022, 10, 847852. [Google Scholar] [CrossRef] [PubMed]
  29. Zhang, Z.H.; Jhaveri, D.J.; Marshall, V.M.; Bauer, D.C.; Edson, J.; Narayanan, R.K.; Robinson, G.J.; Lundberg, A.E.; Bartlett, P.F.; Wray, N.R.; et al. A comparative study of techniques for differential expression analysis on RNA-seq data. PLoS ONE 2014, 9, e103207. [Google Scholar] [CrossRef]
  30. Costa-Silva, J.; Domingues, D.; Lopes, F.M. RNA-seq differential expression analysis: An extended review and a software tool. PLoS ONE 2017, 12, e0190152. [Google Scholar] [CrossRef]
  31. Da Huang, W.; Sherman, B.T.; Lempicki, R.A. Bioinformatics enrichment tools: Paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37, 1–13. [Google Scholar] [CrossRef]
  32. Mukherjee, S.; Banerjee, B.; Karasik, D.; Frenkel-Morgenstern, M. mRNA-lncRNA co-expression network analysis reveals the role of lncrnas in immune dysfunction during severe SARS-CoV-2 infection. Viruses 2021, 13, 402. [Google Scholar] [CrossRef] [PubMed]
  33. Fernandez-Gonzalez, A.; Kourembanas, S.; Wyatt, T.A.; Mitsialis, S.A. Mutation of murine adenylate kinase 7 underlies a primary ciliary dyskinesia phenotype. Am. J. Respir. Cell Mol. Biol. 2009, 40, 305–313. [Google Scholar] [CrossRef] [PubMed]
  34. Turner, M.J.; Abbott-Banner, K.; Thomas, D.Y.; Hanrahan, J.W. Cyclic nucleotide phosphodiesterase inhibitors as therapeutic interventions for cystic fibrosis. Pharmacol. Ther. 2021, 224, 107826. [Google Scholar] [CrossRef]
  35. Wang, Y.; Wang, A.; Zhang, M.; Zeng, H.; Lu, Y.; Liu, L.; Li, J.; Deng, L. Artesunate attenuates airway resistance in vivo and relaxes airway smooth muscle cells in vitro via bitter taste receptor-dependent calcium signalling. Exp. Physiol. 2019, 104, 231–243. [Google Scholar] [CrossRef] [PubMed]
  36. Luo, M.; Ni, K.; Gu, R.; Qin, Y.; Guo, J.; Che, B.; Pan, Y.; Li, J.; Liu, L.; Deng, L. Chemical activation of piezo1 alters biomechanical behaviors toward relaxation of cultured airway smooth muscle cells. Biol. Pharm. Bull. 2023, 46, 1–11. [Google Scholar] [CrossRef] [PubMed]
  37. Yang, C.; Guo, J.; Ni, K.; Wen, K.; Qin, Y.; Gu, R.; Wang, C.; Liu, L.; Pan, Y.; Li, J.; et al. Mechanical ventilation-related high stretch mainly induces endoplasmic reticulum stress and thus mediates inflammation response in cultured human primary airway smooth muscle cells. Int. J. Mol. Sci. 2023, 24, 3811. [Google Scholar] [CrossRef]
  38. Deng, L.; Fairbank, N.J.; Fabry, B.; Smith, P.G.; Maksym, G.N. Localized mechanical stress induces time-dependent actin cytoskeletal remodeling and stiffening in cultured airway smooth muscle cells. Am. J. Physiol. Cell Physiol. 2004, 287, C440–C448. [Google Scholar] [CrossRef]
  39. Deng, L.; Bosse, Y.; Brown, N.; Chin, L.Y.M.; Connolly, S.C.; Fairbank, N.J.; King, G.G.; Maksym, G.N.; Paré, P.D.; Seow, C.Y.; et al. Stress and strain in the contractile and cytoskeletal filaments of airway smooth muscle. Pulm. Pharmacol. Ther. 2009, 22, 407–416. [Google Scholar] [CrossRef]
  40. Liu, G.; Liao, R.; Lv, Y.; Zhu, L.; Lin, Y. Altered expression of mirnas and mrnas reveals the potential regulatory role of mirnas in the developmental process of early weaned goats. PLoS ONE 2019, 14, e0220907. [Google Scholar] [CrossRef]
  41. Yao, Y.; Jiang, C.; Wang, F.; Yan, H.; Long, D.; Zhao, J.; Wang, J.; Zhang, C.; Li, Y.; Tian, X.; et al. Integrative analysis of mirna and mRNA expression profiles associated with human atrial aging. Front. Physiol. 2019, 10, 1226. [Google Scholar] [CrossRef]
  42. Huang, P.; Li, F.; Mo, Z.; Geng, C.; Wen, F.; Zhang, C.; Guo, J.; Wu, S.; Li, L.; Brünner, N.; et al. A comprehensive RNA study to identify circrna and mirna biomarkers for docetaxel resistance in breast cancer. Front. Oncol. 2021, 11, 669270. [Google Scholar] [CrossRef]
  43. Ren, Y.; Yu, G.; Shi, C.; Liu, L.; Guo, Q.; Han, C.; Zhang, D.; Zhang, L.; Liu, B.; Gao, H.; et al. Majorbio cloud: A one-stop, comprehensive bioinformatic platform for multiomics analyses. iMeta 2022, 1, e12. [Google Scholar] [CrossRef]
  44. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. Edger: A bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [PubMed]
  45. Huang, C.; Xiao, X.; Chintagari, N.R.; Breshears, M.; Wang, Y.; Liu, L. MicroRNA and mRNA expression profiling in rat acute respiratory distress syndrome. BMC Med. Genom. 2014, 7, 46. [Google Scholar] [CrossRef] [PubMed]
  46. Szklarczyk, D.; Gable, A.L.; Nastou, K.C.; Lyon, D.; Kirsch, R.; Pyysalo, S.; Doncheva, N.T.; Legeay, M.; Fang, T.; Bork, P.; et al. The STRING database in 2021: Customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 2021, 49, D605–D612. [Google Scholar] [CrossRef] [PubMed]
  47. Gao, X.; Chen, Y.; Chen, M.; Wang, S.; Wen, X.; Zhang, S. Identification of key candidate genes and biological pathways in bladder cancer. PeerJ 2018, 6, e6036. [Google Scholar] [CrossRef] [PubMed]
  48. Ni, K.; Guo, J.; Bu, B.; Pan, Y.; Li, J.; Liu, L.; Luo, M.; Deng, L. Naringin as a plant-derived bitter tastant promotes proliferation of cultured human airway epithelial cells via activation of TAS2R signaling. Phytomedicine 2021, 84, 153491. [Google Scholar] [CrossRef] [PubMed]
  49. Pedley, A.M.; Benkovic, S.J. A new view into the regulation of purine metabolism: The purinosome. Trends Biochem. Sci. 2017, 42, 141–154. [Google Scholar] [CrossRef] [PubMed]
  50. Billington, C.K.; Penn, R.B. Signaling and regulation of g protein-coupled receptors in airway smooth muscle. Respir. Res. 2003, 4, 4. [Google Scholar] [CrossRef]
  51. Smith, P.G.; Deng, L.; Fredberg, J.J.; Maksym, G.N. Mechanical strain increases cell stiffness through cytoskeletal filament reorganization. Am. J. Physiol. Lung Cell Mol. Physiol. 2003, 285, L456–L463. [Google Scholar] [CrossRef] [PubMed]
  52. Cao, G.; Lam, H.; Jude, J.A.; Karmacharya, N.; Kan, M.; Jester, W.; Koziol-White, C.; Himes, B.E.; Chupp, G.L.; An, S.S.; et al. Inhibition of ABCC1 decreases cAMP egress and promotes human airway smooth muscle cell relaxation. Am. J. Respir. Cell Mol. Biol. 2022, 66, 96–106. [Google Scholar] [CrossRef] [PubMed]
  53. Huang, D.W.; Sherman, B.T.; Tan, Q.; Kir, J.; Liu, D.; Bryant, D.; Guo, Y.; Stephens, R.; Baseler, M.W.; Lane, H.C.; et al. DAVID bioinformatics resources: Expanded annotation database and novel algorithms to better extract biology from large gene lists. Nucleic Acids Res. 2007, 35, W169–W175. [Google Scholar] [CrossRef]
  54. Wang, G.; Zou, R.; Liu, L.; Wang, Z.; Zou, Z.; Tan, S.; Xu, W.; Fan, X. A circular network of purine metabolism as coregulators of dilated cardiomyopathy. J. Transl. Med. 2022, 20, 532. [Google Scholar] [CrossRef] [PubMed]
  55. Esther, C.R.; Peden, D.B.; Alexis, N.E.; Hernandez, M.L. Airway purinergic responses in healthy, atopic nonasthmatic, and atopic asthmatic subjects exposed to ozone. Inhal. Toxicol. 2011, 23, 324–330. [Google Scholar] [CrossRef] [PubMed]
  56. Verbrugge, S.J.; de Jong, J.W.; Keijzer, E.; Vazquez de Anda, G.; Lachmann, B. Purine in bronchoalveolar lavage fluid as a marker of ventilation-induced lung injury. Crit. Care Med. 1999, 27, 779–783. [Google Scholar] [CrossRef]
  57. Esther, J.C.R.; Alexis, N.E.; Clas, M.L.; Lazarowski, E.R.; Donaldson, S.H.; Ribeiro, C.M.P.; Moore, C.G.; Davis, S.D.; Boucher, R.C. Extracellular purines are biomarkers of neutrophilic airway inflammation. Eur. Respir. J. 2008, 31, 949–956. [Google Scholar] [CrossRef]
  58. Wen, Y.; Zhang, X.; Larsson, L. Metabolomic profiling of respiratory muscles and lung in response to long-term controlled mechanical ventilation. Front. Cell Dev. Biol. 2022, 10, 849973. [Google Scholar] [CrossRef]
  59. Takahara, N.; Ito, S.; Furuya, K.; Naruse, K.; Aso, H.; Kondo, M.; Sokabe, M.; Hasegawa, Y. Real-time imaging of ATP release induced by mechanical stretch in human airway smooth muscle cells. Am. J. Respir. Cell Mol. Biol. 2014, 51, 772–782. [Google Scholar] [CrossRef]
  60. Zhou, R.; Wang, R.; Qin, Y.; Ji, J.; Xu, M.; Wu, W.; Chen, M.; Wu, D.; Song, L.; Shen, H.; et al. Mitochondria-related mir-151a-5p reduces cellular ATP production by targeting cytb in asthenozoospermia. Sci. Rep. 2015, 5, 17743. [Google Scholar] [CrossRef]
  61. Lorès, P.; Coutton, C.; El Khouri, E.; Stouvenel, L.; Givelet, M.; Thomas, L.; Rode, B.; Schmitt, A.; Louis, B.; Sakheli, Z.; et al. Homozygous missense mutation L673P in adenylate kinase 7 (AK7) leads to primary male infertility and multiple morphological anomalies of the flagella but not to primary ciliary dyskinesia. Hum. Mol. Genet. 2018, 27, 1196–1211. [Google Scholar] [CrossRef] [PubMed]
  62. Alvarez-Santos, M.D.; Alvarez-Gonzalez, M.; Eslava-De-Jesus, E.; Gonzalez-Lopez, A.; Pacheco-Alba, I.; Perez-Del-Valle, Y.; Rojas-Madrid, R.; Bazan-Perkins, B. Role of airway smooth muscle cell phenotypes in airway tone and obstruction in guinea pig asthma model. Allergy Asthma Clin. Immunol. 2022, 18, 3. [Google Scholar] [CrossRef] [PubMed]
  63. O’Sullivan, M.J.; Jang, J.H.; Panariti, A.; Bedrat, A.; Ijpma, G.; Lemos, B.; Park, J.A.; Lauzon, A.M.; Martin, J.G. Airway epithelial cells drive airway smooth muscle cell phenotype switching to the proliferative and pro-inflammatory phenotype. Front. Physiol. 2021, 12, 687654. [Google Scholar] [CrossRef] [PubMed]
  64. Frank, J.A.; Matthay, M.A. Science review: Mechanisms of ventilator-induced injury. Crit. Care 2003, 7, 233. [Google Scholar] [CrossRef] [PubMed]
  65. Wilson, M.R.; Choudhury, S.; Goddard, M.E.; O’Dea, K.P.; Nicholson, A.G.; Takata, M. High tidal volume upregulates intrapulmonary cytokines in an in vivo mouse model of ventilator-induced lung injury. J. Appl. Physiol. 2003, 95, 1385–1393. [Google Scholar] [CrossRef]
  66. Yang, G.; Hamacher, J.; Gorshkov, B.; White, R.; Sridhar, S.; Verin, A.; Chakraborty, T.; Lucas, R. The dual role of TNF in pulmonary edema. J. Cardiovasc. Dis. Res. 2010, 1, 29–36. [Google Scholar] [CrossRef]
  67. Luan, X.; Le, Y.; Jagadeeshan, S.; Murray, B.; Carmalt, J.L.; Duke, T.; Beazley, S.; Fujiyama, M.; Swekla, K.; Gray, B.; et al. cAMP triggers Na+ absorption by distal airway surface epithelium in cystic fibrosis swine. Cell Rep. 2021, 37, 109795. [Google Scholar] [CrossRef]
  68. Cavanaugh, K.J., Jr.; Oswari, J.; Margulies, S.S. Role of stretch on tight junction structure in alveolar epithelial cells. Am. J. Respir. Cell Mol. Biol. 2001, 25, 584–591. [Google Scholar] [CrossRef]
Cells 13 00110 g001
Figure 1. Schematic diagram of the experimental protocol to identify high stretch-responsive events and the underlying regulatory miRNA–mRNA interactions in human airway smooth muscle cells (ASMCs). (A) The control and experimental ASMCs were cultured on an elastic membrane in either static (left) or high-stretch (right, 13% cyclic strain, 0.5 Hz) conditions, respectively. (B) Flow chart of high-stretch treatment, whole genome-wide RNA-Seq, bioinformatics analysis, and functional analysis to identify high stretch-induced enriched events and miRNA–mRNA interactions in cultured ASMCs.
Figure 1. Schematic diagram of the experimental protocol to identify high stretch-responsive events and the underlying regulatory miRNA–mRNA interactions in human airway smooth muscle cells (ASMCs). (A) The control and experimental ASMCs were cultured on an elastic membrane in either static (left) or high-stretch (right, 13% cyclic strain, 0.5 Hz) conditions, respectively. (B) Flow chart of high-stretch treatment, whole genome-wide RNA-Seq, bioinformatics analysis, and functional analysis to identify high stretch-induced enriched events and miRNA–mRNA interactions in cultured ASMCs.
Cells 13 00110 g001
Cells 13 00110 g002
Figure 2. Identification of target differentially expressed mRNAs (DE-mRNAs) in ASMCs cultured under static or high-stretch conditions. (A) Overlapped genes between targets of down-regulated DE-miRNAs and up-regulated DE-mRNAs (up panel) and overlapped genes between targets of up-regulated DE-miRNAs and down-regulated DE-mRNAs (down panel). (B) Volcano plot of 283 target DE-mRNAs. Black dots indicate the genes that were differently expressed in ASMCs cultured in high-stretch groups compared to static groups with p < 0.01 and |log2FC| > 2. Red dots indicate the top 10 upregulated mRNAs and green dots indicate the top 10 down-regulated mRNAs.
Figure 2. Identification of target differentially expressed mRNAs (DE-mRNAs) in ASMCs cultured under static or high-stretch conditions. (A) Overlapped genes between targets of down-regulated DE-miRNAs and up-regulated DE-mRNAs (up panel) and overlapped genes between targets of up-regulated DE-miRNAs and down-regulated DE-mRNAs (down panel). (B) Volcano plot of 283 target DE-mRNAs. Black dots indicate the genes that were differently expressed in ASMCs cultured in high-stretch groups compared to static groups with p < 0.01 and |log2FC| > 2. Red dots indicate the top 10 upregulated mRNAs and green dots indicate the top 10 down-regulated mRNAs.
Cells 13 00110 g002
Cells 13 00110 g003
Figure 3. Purine metabolism was enriched in the KEGG pathway and GO enrichment analysis of the 283 target DE-mRNAs. (A) Significant enrichment of the KEGG pathway of the 283 target DE-mRNAs. (B) Significant enrichment GO of the 283 target DE-mRNAs. Purple bars stand for the BP, and green bars stand for the MF. Box indicates that these events were related to purine metabolism. The circle indicates the number of enriched genes in each term.
Figure 3. Purine metabolism was enriched in the KEGG pathway and GO enrichment analysis of the 283 target DE-mRNAs. (A) Significant enrichment of the KEGG pathway of the 283 target DE-mRNAs. (B) Significant enrichment GO of the 283 target DE-mRNAs. Purple bars stand for the BP, and green bars stand for the MF. Box indicates that these events were related to purine metabolism. The circle indicates the number of enriched genes in each term.
Cells 13 00110 g003
Cells 13 00110 g004
Figure 4. Purine metabolism was enriched in the first module of the PPI network. (A) The PPI network construction and module analysis. The blue nodes represent the up-regulated target DE-mRNAs. The green nodes represent the down-regulated target DE-mRNAs. The size of the nodes indicates the degree of the target DE-mRNAs, and the lines represent the interaction between gene-encoded proteins. (B) Significant enrichment of the KEGG pathway of the 38 DE-mRNAs in the first module of the PPI network. Box indicates that these events were related to purine metabolism. (C) Significant enrichment GO of the 38 DE-mRNAs in the first module of the PPI network. Green bars stand for the BP, and red bars stand for the MF. Box indicates that these events were related to purine metabolism. The circle indicates the number of enriched genes in each term.
Figure 4. Purine metabolism was enriched in the first module of the PPI network. (A) The PPI network construction and module analysis. The blue nodes represent the up-regulated target DE-mRNAs. The green nodes represent the down-regulated target DE-mRNAs. The size of the nodes indicates the degree of the target DE-mRNAs, and the lines represent the interaction between gene-encoded proteins. (B) Significant enrichment of the KEGG pathway of the 38 DE-mRNAs in the first module of the PPI network. Box indicates that these events were related to purine metabolism. (C) Significant enrichment GO of the 38 DE-mRNAs in the first module of the PPI network. Green bars stand for the BP, and red bars stand for the MF. Box indicates that these events were related to purine metabolism. The circle indicates the number of enriched genes in each term.
Cells 13 00110 g004
Cells 13 00110 g005
Figure 5. High stretch modulated purine metabolism. (A) The miRNA–mRNA interactions related to purine metabolism response to high stretch of ASMCs. The red and green colors represent up-regulation and down-regulation, respectively. The bar graph shows the relative cAMP (B) and ATP (C) levels inside ASMCs cultured in static or high-stretch conditions. Each experiment was repeated three times (n = 3). Data are presented as means ± S.E.M. * p < 0.05 and ** p < 0.01 compared to static groups.
Figure 5. High stretch modulated purine metabolism. (A) The miRNA–mRNA interactions related to purine metabolism response to high stretch of ASMCs. The red and green colors represent up-regulation and down-regulation, respectively. The bar graph shows the relative cAMP (B) and ATP (C) levels inside ASMCs cultured in static or high-stretch conditions. Each experiment was repeated three times (n = 3). Data are presented as means ± S.E.M. * p < 0.05 and ** p < 0.01 compared to static groups.
Cells 13 00110 g005
Cells 13 00110 g006
Figure 6. miR-370-5p partially reversed high stretch-induced change in cAMP and ATP levels inside ASMCs cultured in static or high-stretch conditions. (A) The cellular activity of ASMCs treated with miR-370-5p mimics or inhibitor. (B) Relative expression levels of miR-370-5p in ASMCs treated with miR-370-5p mimics or inhibitor. (C,D) The mRNA expression of PDE4D, AK7 of ASMCs treated with miR-370-5p mimics or inhibitor, respectively. (E,F) The relative cAMP, ATP level inside ASMCs cultured with miR-370-5p mimics or inhibitor, respectively. (G,H) The relative cAMP, ATP level inside ASMCs cultured in static or high-stretch conditions with or without miR-370-5p mimics. Each experiment was repeated three times (n = 3). Data are presented as means ± S.E.M. Markers denote not significant (n.s) or significant difference from control conditions on NC mimics or static group (* p < 0.05), NC inhibitor (# p < 0.05), high-stretch group with NC mimics ($ p < 0.05), ** p < 0.01, ## p < 0.01, $$ p < 0.01.
Figure 6. miR-370-5p partially reversed high stretch-induced change in cAMP and ATP levels inside ASMCs cultured in static or high-stretch conditions. (A) The cellular activity of ASMCs treated with miR-370-5p mimics or inhibitor. (B) Relative expression levels of miR-370-5p in ASMCs treated with miR-370-5p mimics or inhibitor. (C,D) The mRNA expression of PDE4D, AK7 of ASMCs treated with miR-370-5p mimics or inhibitor, respectively. (E,F) The relative cAMP, ATP level inside ASMCs cultured with miR-370-5p mimics or inhibitor, respectively. (G,H) The relative cAMP, ATP level inside ASMCs cultured in static or high-stretch conditions with or without miR-370-5p mimics. Each experiment was repeated three times (n = 3). Data are presented as means ± S.E.M. Markers denote not significant (n.s) or significant difference from control conditions on NC mimics or static group (* p < 0.05), NC inhibitor (# p < 0.05), high-stretch group with NC mimics ($ p < 0.05), ** p < 0.01, ## p < 0.01, $$ p < 0.01.
Cells 13 00110 g006
Table 1. High stretch-induced differentially expressed miRNAs (DE-miRNAs) in cultured ASMCs (stretch/static, count ≥ 10; p < 0.05; Fold change > 2).
Table 1. High stretch-induced differentially expressed miRNAs (DE-miRNAs) in cultured ASMCs (stretch/static, count ≥ 10; p < 0.05; Fold change > 2).
No.DE-miRNALog2FC(Stretch/Static)padjustRegulate
1miR-121362.253.07 × 10−62up
2miR-192-5p1.403.22 × 10−15up
3miR-146a-5p1.295.34 × 10−3up
4miR-194-5p1.181.31 × 10−17up
5miR-29b-3p1.162.23 × 10−7up
6miR-148a-3p1.167.06 × 10−20up
7miR-137-3p1.108.46 × 10−9up
8miR-370-5p−1.012.93 × 10−4down
9miR-27b-5p−1.231.63 × 10−11down
10miR-543−1.302.80 × 10−19down
11miR-485-3p−1.479.13 × 10−19down
12miR-335-3p−2.164.85 × 10−7down
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Luo, M.; Wang, C.; Guo, J.; Wen, K.; Yang, C.; Ni, K.; Liu, L.; Pan, Y.; Li, J.; Deng, L. High Stretch Modulates cAMP/ATP Level in Association with Purine Metabolism via miRNA–mRNA Interactions in Cultured Human Airway Smooth Muscle Cells. Cells 2024, 13, 110. https://doi.org/10.3390/cells13020110

AMA Style

Luo M, Wang C, Guo J, Wen K, Yang C, Ni K, Liu L, Pan Y, Li J, Deng L. High Stretch Modulates cAMP/ATP Level in Association with Purine Metabolism via miRNA–mRNA Interactions in Cultured Human Airway Smooth Muscle Cells. Cells. 2024; 13(2):110. https://doi.org/10.3390/cells13020110

Chicago/Turabian Style

Luo, Mingzhi, Chunhong Wang, Jia Guo, Kang Wen, Chongxin Yang, Kai Ni, Lei Liu, Yan Pan, Jingjing Li, and Linhong Deng. 2024. "High Stretch Modulates cAMP/ATP Level in Association with Purine Metabolism via miRNA–mRNA Interactions in Cultured Human Airway Smooth Muscle Cells" Cells 13, no. 2: 110. https://doi.org/10.3390/cells13020110

APA Style

Luo, M., Wang, C., Guo, J., Wen, K., Yang, C., Ni, K., Liu, L., Pan, Y., Li, J., & Deng, L. (2024). High Stretch Modulates cAMP/ATP Level in Association with Purine Metabolism via miRNA–mRNA Interactions in Cultured Human Airway Smooth Muscle Cells. Cells, 13(2), 110. https://doi.org/10.3390/cells13020110

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop