- Open Access
N7-methylguanosine modification of lncRNAs in a rat model of hypoxic pulmonary hypertension: a comprehensive analysis
BMC Genomics volume 23, Article number: 33 (2022)
Long non-coding RNAs (lncRNAs) play a critical role in the pathogenesis of hypoxic pulmonary hypertension (HPH). The role of N7-methylguanosine (m7G) modification in lncRNAs has received increased attentions in recent years. However, the m7G-methylation of lncRNA in HPH has yet to be determined. We have therefore performed a transcriptome-wide analysis of m7G lncRNAs in HPH.
Differentially-expressed m7Gs were detected in HPH, and m7G lncRNAs were significantly upregulated compared with non-m7G lncRNAs in HPH. Importantly, this was the first time that the upregulated m7G lncXR_591973 and m7G lncXR_592398 were identified in HPH.
This study provides the first m7G transcriptome-wide analysis of HPH. Importantly, two HPH-associated m7G lncRNAs were identified, although their clinical significance requires further validation.
Pulmonary hypertension (PH) is a lethal disease causing increased vascular resistance and eventually leading to right ventricular failure and death . Hypoxic pulmonary hypertension (HPH) is classified as a group III PH based on pathogenesis. It is a progressive disease induced by chronic hypoxia, which normally results from severe interstitial lung diseases and chronic obstructive pulmonary disease (COPD) . The pathophysiology of HPH is characterized by pulmonary arterial remodelling, mediated via unusual proliferation of pulmonary artery endothelial cells (PAECs) and pulmonary artery smooth muscle cells (PASMCs), but also activation of quiescent fibroblasts [3, 4]. However, the underlying molecular mechanisms remain imperfectly understood due to the complexity and malignancy of HPH.
In recent years, long non-coding RNAs (lncRNAs) have received extensive attention because of their critical role in the regulation of biological phenomena such as cell proliferation, apoptosis, migration and invasion [5,6,7,8]. lncRNAs generally refer to a class of transcripts longer than 200 nucleotides (nt) that lack protein-coding ability . They interact with proteins, interfere with microRNA by acting as molecular sponges, modify the epigenome and interfere with gene expression by binding with gene promoters [10,11,12,13]. The crucial role of lncRNAs in HPH pathogenesis in idiopathic PH was highlighted by the up-regulation of 2004 lncRNAs and down-regulation of 507 lncRNAs . Subsequently, multiple lncRNAs have been implicated in the pathogenesis of HPH, such as lncRNA MEG3 and lncRNA Tug1. It is noteworthy that their dysregulated expression has been observed in HPH [15,16,17]. RNA modification plays a critical role in the regulation of gene expression . Nevertheless, post-transcription modification of lncRNAs in HPH remains largely unexplored.
N7-Methylguanosine (m7G) is a ubiquitous post-transcriptional modification of mRNA and lncRNA in eukaryotes, and it is essential for efficient gene expression and cell viability [19, 20]. During transcription, m7G is incorporated at the 5′ end, modulating the various events in the mRNA cycle including RNA splicing, polyadenylation, nuclear export and translation . Although evidence increasingly indicates that m7G modification is closely associated with the initiation and progression of various diseases, the RNA m7G-methylation profile of HPH has yet to be reported. This study provides the first transcriptome-wide analysis of m7G profile in HPH and control group, demonstrating the enormous diversity of m7G modification patterns in these two groups. It is hoped that this study will facilitate further investigation into the potential role of m7G modification in HPH pathogenesis.
General features of m7G methylation of lncRNAs in HPH and control samples
A rat model of hypoxia-induced PH was successfully established . Pulmonary tissues from HPH and control (normoxia, N) rats were used for transcriptome-wide m7G-MeRIP-seq and RNA-seq analyses. The HPH tissues revealed 2685 m7G peaks, while the control rat tissues showed 2644 m7G peaks. The HPH and control tissues shared 255 peaks, which accounted for only 5.0% of all peaks in the two groups (Fig. 1a; Additional file 1: Data S1). The low percentage of shared m7G peaks within the lncRNAs indicated differences in the m7G pattern between the two groups. Additionally, the level of m7G in the total lncRNA derived from HPH rats was lower than that in the control group (Fig. 1b).
Analysis of the distribution of m7G peaks in each lncRNA revealed no significant difference between the two groups. Approximately 80% of the modified lncRNAs showed unique m7G peaks, while almost 10% of the modified lncRNAs revealed two m7G peaks (Fig. 1c).
To further elucidate the m7G methylation pattern of lncRNAs in the HPH and control groups, the samples were divided into the following six categories based on the positional relationship between the m7G methylated lncRNAs and mRNAs: bidirectional, exon sense overlap, intron sense overlap, intron antisense overlap, natural antisense overlap, and intergenic overlap. The results revealed that the majority of lncRNAs contained m7G sites with exon sense overlap. In deed, HPH samples accounted for 53%, while control samples constituted 54% of such lncRNAs. Most lncRNAs carried intergenic m7G sites, with HPH samples accounting for 23% and control samples constituting 25% of such lncRNAs. Only approximately 2% of the methylation sites were located in the intronic antisense group (Fig. 1d-e).
Further analysis of the m7G methylated lncRNAs showed that more than 30% of total lncRNA contained two exons in both HPH and control groups (Fig. 1f-g). Similarly, m7G and non-m7G lncRNAs were mainly encoded by two exons, suggesting that the two exons in the HPH and control groups contributed to most of the m7G methylation in lncRNAs.
Length and distribution of differentially-methylated lncRNAs
As displayed in the heatmap of m7G lncRNA expression (Fig. 2a), 353 differentially methylated m7G sites were identified across 344 lncRNAs. Of these, 55% (193/353) were significantly hypermethylated and 45% (160/353) were significantly hypomethylated (fold change ≥2 and P ≤ 0.00001; Fig. 2b; Additional file 2: Data S2). The 353 differentially methylated m7G sites were located in 344 lncRNAs, including hypermethylated (189) and hypomethylated (155) types (Fig. 2b).
Analysis of the length of the differentially m7G-methylated lncRNAs revealed differences between hypermethylated and hypomethylated types (Fig. 2c). The length of the hypermethylated m7G lncRNAs was primarily 1–1000 bp (24%) and 2000–3000 bp (23%), while the length of the hypomethylated m7G lncRNAs was mostly in the range of 1000 to 2000 bp (29%).
To analyse the distribution of the differentially m7G-methylated lncRNAs across chromosomes, the enrichment level of m7G methylated lncRNAs on each chromosome was further determined. As shown in Fig. 2d, the hypomethylated m7G lncRNAs were primarily located on chromosomes 1 (13.99%), 5 (7.25%) and 4 (6.74%) and the hypomethylated m7G lncRNAs were primarily located on chromosomes 1 (15.63%), 10 (10.0%) and 3 (6.88%).
Functional analysis of genes near differentially methylated lncRNAs
To elucidate the role of differentially methylated lncRNAs in the occurrence and development of HPH, GO and KEGG pathway analyses of genes located near differentially methylated lncRNAs were performed.
In the biological process (BP) category, GO results revealed that genes near hypermethylated m7G lncRNAs primarily participated in complement activation, negative regulation of smooth muscle contraction and the protein activation cascade. In the cellular component (CC) category, genes were primarily associated with the neuronal cell body, intracellular cellular component and postsynaptic membrane. In terms of molecular function (MF), genes were primarily involved in binding, particularly metal and calcium ion (Fig. 3a). Genes located near hypomethylated lncRNAs were primarily involved in the following BPs: protein O-linked glycosylation, muscle cell differentiation and immunoglobulin secretion. In terms of CC, genes were associated with intracellular organelle, the clathrin vesicle coat and cell cellular component. In terms of MF, genes primarily participated in binding, phosphatidylserine binding and acetylgalactosamine transferase activity (Fig. 3b).
KEGG pathway analysis revealed the four most important signalling pathways associated with genes located near hypermethylated lncRNAs were associated with Staphylococcus aureus infection, the cGMP-PKG signalling pathway, complement and coagulation cascades, and Fc gamma R-mediated phagocytosis (Fig. 3c). The four most important signalling pathways for genes located near hypomethylated lncRNAs were associated with graft-versus-host disease, type 1 diabetes mellitus, autoimmune thyroid disease and viral myocarditis (Fig. 3d).
These results suggested that m7G-modified lncRNAs may affect the occurrence and development of HPH via BPs, cellular composition, MF and signalling pathways.
RNA-seq analysis of differentially expressed lncRNAs in HPH
Figure 4a shows a heatmap of RNA-seq data in the HPH and control groups, indicating the differential expression of lncRNAs in the two groups. A total of 111 lncRNAs were differentially expressed in HPH, including 90 up-regulated and 21 down-regulated types (Fig. 4b).
To further analyse the effects of genes located near differentially expressed lncRNAs, GO enrichment and KEGG pathways analysis of genes near lncRNAs were performed. GO analysis revealed that genes located near up-regulated lncRNAs were significantly enriched in terms of positive regulation of immune effector processes, catalytic step 2 spliceosome and unfolded protein binding (Fig. 4c). Genes near down-regulated lncRNAs were significantly associated with blastocyst formation, lysozyme activity and intrinsic components of organelle membranes (Fig. 4d).
KEGG pathway analysis revealed that genes near up-regulated lncRNAs were primarily involved in antigen processing and presentation, endocytosis, and phagocytosis (Fig. 4e). Genes located near down-regulated lncRNAs were mainly involved in the Wnt signalling pathway. Together, these results indicated that genes near differentially expressed lncRNAs may be related to HPH (Fig. 4f).
Association of m7G methylation and the expression of lncRNAs
A total of 972 m7G lncRNAs were identified in the HPH group and 918 m7G lncRNAs were identified in the control group. Among these, 767 m7G-modified lncRNAs were detected in both groups (Fig. 5a). The expression of 767 common m7G lncRNAs was investigated to analyse the effect of hypoxia-regulated m7G methylation. As shown in Fig. 5b, 189 hyper-methylated and 155 hypo-methylated m7G lncRNAs were found in the HPH group compared with the control group (fold change > 2, P < 0.0001).
To further determine whether m7G methylation affected lncRNA expression, all the expressed lncRNAs were divided into m7G and non-m7G categories. The log two-fold change (log2FC) values of these lncRNAs were determined and a cumulative curve was generated. The proportion of m7G-modified lncRNAs was larger than that of non-m7G lncRNAs. These results indicated that m7G significantly up-regulated the expression of lncRNAs under hypoxic conditions (Fig. 5c, P = 0.0267).
The combined results of methylation sequencing and RNA sequencing revealed that only two up-regulated lncRNAs were hypermethylated, which originated from chr4: 7450738–7,450,981 (LOC102554730, XR_591973) and chr13: 95420341–95,420,800 (LOC102555374, XR_592398). However, no down-regulated lncRNAs were hypomethylated.
Construction of lncRNA-miRNA-mRNA competing endogenous RNA (ceRNA) network and identification of differentially expressed m7G lncRNAs
To explore the potential mRNAs, which were regulated by the lncRNAs, lncXR_592398 and lncXR_592398 were used to build a ceRNA network. Among the predicted miRNAs, miR-190a-5p, miR-145-5p, miR-29a-3p, let-7 g-5p, miR-127-3p and miR-221-3p were selected to construct the ceRNA network because they were related to PH . The top 10 genes for each of these miRNAs were predicted using MirTarget and Miranda programs (Table 1). Thus, the network consisted of PH-associated miRNAs combined with the two differentially expressed m7G lncRNAs and the top 10 mRNAs that bound to the miRNAs, including two lncRNAs, six miRNAs and 58 mRNAs. This ceRNA network clearly suggests that lncRNAs regulated miRNAs and mRNAs. For example, LOC102555374 (XR_592398) bound to rno-let-7 g-5p to regulate Tmem198 expression (Fig. 5d).
Up-regulation of m7G XR_591973 and m7G XR_592398 in hypoxic PASMCs
To further corroborate the results, the expression of the two differentially expressed m7G lncRNAs in PASMCs treated with hypoxia (1% O2) and normoxia (21% O2), respectively, was validated. Using cDNA from the PASMCs as templates, lncXR_591973 and lncXR_592398 were amplified by divergent primers. MeRIP-qRT-PCR was used to determine whether lncXR_591973 and lncXR_592398 were modified by m7G. As shown in Fig. 6, the levels of m7G lncXR_591973 and m7G lncXR_592398 were significantly increased in hypoxic PASMCs (P = 0.001 and P = 0.002, respectively), demonstrating the reliability of the transcriptome-wide m7G-MeRIP-seq analysis.
In the present study, transcriptome-wide mapping of m7G lncRNAs in HPH was performed for the first time to elucidate the role of m7G lncRNAs in the pathogenesis of PH. The m7G modification pattern in HPH samples was distinct from that of normal controls, with a tendency towards lower total m7A abundance in the HPH group. In addition, the two groups shared only 5% m7G peaks. A recent study found that internal m7G/G levels were quite low (range from ~ 0.02% to ~ 0.05%) in human and mouse cell lines within mammalian mRNA, which represent roughly 5–10% of their respective m6A/A ratios. Further, m7G was detected internally within mammalian mRNA. Although the m7G MeRIP-seq method has been widely used, the method still has limitations in detecting m7G internally . These limitations might result in the minor overlap in m7G peaks found in this study.
The transcriptome-wide analysis revealed that m7G modifications were associated with a positive correlation with lncRNA expression in HPH samples. This result was further corroborated by m7G MeRIP-qRT-PCR assays. Similarly, the abundance of m6A modification has been positively associated with lncRNA expression . Kouzarideset al. recently reported that m7G modifications promoted human miRNA processing by directly affecting the secondary structure of pri-miRNAs, which suppressed cell migration . Additionally, m7G modification was also demonstrated to be essential for efficient pre-mRNA splicing, translation and protein-RNA interactions . Further, accumulating evidence suggests that m7G modification is significantly related to the development of multiple diseases, including cancers, infections and stem cell-associated diseases .
The KEGG analysis in the present study revealed a significant role of up-regulated m7G lncRNAs in vascular smooth muscle contraction pathway. Pulmonary vasoconstriction is known to play an important role in the pathogenesis of HPH and represents one of the factors contributing to increased PASMC contraction . Clinical studies have also shown that vasodilating agents ameliorate PH and cardiac output of patients diagnosed with HPH, which demonstrates the key role of vasoconstriction in HPH pathophysiology [29, 30]. A previous study further reported that inhibiting PASMC contraction in response to hypoxia significantly decreased the RVP of PH models .
lncRNAs act as ceRNAs to indirectly regulate mRNAs via shared microRNAs, which represents a novel layer of RNA crosstalk and plays a critical role in the development of various diseases . m7G plays a critical role in regulating RNA processing, metabolism and function , suggesting that m7G also regulates the lncRNA-miRNA-mRNA co-expression network. Analysis of the differentially expressed m7G lncRNAs revealed an overlap of 11 pseudogenes between the identified lncRNAs. Among the overlapped pseudogenes, Smad3 is a key downstream element in transforming growth factor (TGF)-β signaling pathways . In recent years, the critical role of TGF-β signaling pathways in PH has been better elucidated. Enhanced TGF-β signaling has been associated with PH accompanied by smooth muscle hypertrophy, perivascular fibrosis, and extracellular matrix remodeling . Another overlapped gene, inositol monophosphatase domain containing 1(Impad1), a novel sulfotransferase that is normally located in golgi bodies, has received increased attention in cancer, as it converts PAP into AMP in human fibroblasts . Increasing evidence indicates that in cancer, Impad1 acted in the ER-Golgi pathway by altering Golgi-mediated secretion of MMPs, resulting in a pro-invasive and metastatic tumour phenotype . Additionally, Impad1 overexpression was shown to inhibit Complex I activity, reducing ROS production in cancer cells, and promoting tumor cell invasion . As a result, the function of Impad1 is diverse and needs further exploration in HPH. Among the non-overlapped pseudogenes, nuclear factor of activated T cells 5 (Nfat5), the Ca2+ / calcineurin sensitive transcription factor, was reported to perpetuate the elevation in cytosolic calcium and promote apoptosis resistance of PASMCs. Increased cytosolic calcium contributes to the contractile, hyperproliferative, and anti-apoptotic phenotype of PASMCs, which leads to PH pathogenesis . However, Wnt5a, a key factor of Wnt signaling pathways, plays a key role in the Wnt signaling pathways in PAH pathobiology. Evidence suggests that alterations in Wnt pathway activation lead to PH by preventing proper small vessel regeneration and allowing excessive PASMC growth . In general, these downstream targets regulated by lncRNA-miRNA of interest were mostly enriched in PH-associated TGF-β, Wnt and Nfat signaling pathways. Our study reported for the first time that m7G influenced lncRNAs, thus regulating the binding of miRNA, resulting in the activation of TGF-β, Wnt and Nfat signaling pathways.
The main limitations of this study were the lack of clinical samples and the inability to determine the precise mechanism of hypoxia-induced m7G modification. Thus, studies focusing on the underlying mechanism of m7G modification, as well as clinical transformation, are required to establish the role of m7G lncRNAs in the pathogenesis of HPH.
This study provides the first m7G transcriptome-wide map of HPH. It demonstrated that the m7G modification pattern in HPH samples is distinct from that of controls, with a tendency towards lower total m7A abundance in the HPH group. Further, two significant m7G lncRNAs were found to be associated with HPH for the first time, and both of them are validated in PASMC cells. The role of m7G lncRNAs in HPH and its underlying mechanism require further exploration. It is hoped that this will be the start of studies investigating m7G functions and modification in HPH. In the future, the integrated analysis of transcription and translation studies might contribute to the elucidation of the mechanism of HPH. Additional evidence derived from large patient populations with HPH will improve the value of this study and identify specific HPH-related m7G methylation sites.
Establishment of HPH rat model
4-week male Sprague-Dawley (SD) rats weighing 180 to 200 g each were used in this study and were supplied by the Animal Experimental Centre, Zhejiang University, China. Rats were randomly allocated to the HPH and the control group (n = 6 in each group). The rat model for this study was established based on our previous study , and the same animal group was used in the two studies. The construction and subsequent verification of a hypoxia-induced PH rat model has been described in our previous study . All rats were sacrificed by exposure to 100% CO2 in a confined and transparent euthanasia device, and heart and lung tissues were removed for the following experiments. All protocols and procedures were approved by the Ethics Committee of Zhejiang University, China and were conducted in accordance with the guidelines of the National Institutes of Health on the care and use of animals.
Isolation and hypoxia treatment of primary PASMCs
The isolation of primary PASMCs from rat pulmonary artery was performed as described previously . PASMCs were cultured in Dulbecco’s modified Eagle medium supplemented with 10% foetal bovine serum (FBS). Before experimental use, the cells were incubated under hypoxic (1% O2) or normoxic (21% O2) conditions at 37 °C for 48 h. PASMCs were grown to no more than 90% confluence and used between passages four and seven.
RNA isolation and RNA-seq analysis
Total RNA from lungs (1 g) of control and HPH rats was extracted with TRIzol reagent (Invitrogen, CA, USA). The extracted RNA was purified using a NEBNext rRNA Depletion Kit (New England Biolabs, MA, USA). Next, RNA libraries were prepared using the NEBNext Ultra Directional RNA Library Prep kit (New England Biolabs, MA, USA). Library sequencing was performed on an IlluminaHiseq instrument with 150 bp paired-end reads.
MeRIP and m7G-seq analysis
The isolated RNA was subjected to immunoprecipitation using the m7G-IPKit (GenSeq, Beijing, China) according to the manufacturer’s instructions. RNA was randomly fragmented to about 200 nt using RNA fragmentation reagents (Millipore Sigma, USA). Protein A/G beads (Thermo Fisher Scientific, USA) were coupled to the m7G antibody via rotation at room temperature for 1 h. The RNA fragments were incubated with the bead-linked antibodies and rotated at 4 °C for 4 h. Then, the bound RNA was eluted from the beads using Proteinase K and purified by phenol-chloroform extraction. Purified RNA was used for RNA-seq library generation with a NEBNext® Ultra II Directional RNA Library Prep Kit (New England Biolabs, MA, USA). Paired-end reads were harvested from an IlluminaHiSeq 4000 sequencer.
Construction of a ceRNA network
The lncRNA-microRNA (miRNA) -mRNA ceRNA network was constructed based on the theory that lncRNA acts as an miRNA sponge to further regulate mRNA. The binding sites of miRNA on differentially expressed m7G lncRNAs as well as the target mRNA of the miRNA were predicted with Target Scan and miRanda. The miRNAs related to PH and the top 10 mRNAs binding to each selected miRNA were retained for further network analysis. Pearson correlation analysis was used to assess the ceRNA network construction. Cytoscape v3.6.0 was used to visualise the ceRNA network.
Methylated RNA immunoprecipitation and MeRIP-qPCR
M7G immunoprecipitation (MeRIP) was performed to measure the specific m7G modification of lncRNA. Total RNA was isolated from hypoxia- and normoxia-treated PASMCs. RNA (18 μg) was fragmented to approximately 200 nt-long fragments in an RNA fragmentation buffer. The fragmented RNA was purified using an RNaseMiniElute Kit according to the manufacturer’s instructions. Protein A/G beads (50 μL) were coupled to m7G antibody by rotating at room temperature for 30 min. The RNA fragments were then incubated with the bead-linked antibodies and rotated at 4 °C for 2 h. The immunopurified RNA was purified and first-strand cDNA synthesis was performed by real-time quantitative PCR (QPCR). All the samples were run with the qPCR SYBR Green Master Mix (CloudSeq). The primers used for the detection of m7G-enriched lncRNA are presented in Table 2.
Sequencing data analysis and statistical analysis
The quality of harvested paired-end reads were controlled by Q30. Cutadapt (v1.9.3) was used for 3′ adaptor trimming and the removal of low-quality reads. Clean reads from all libraries were aligned with the reference genome (UCSC RN5) using HISAT2 (v2.0.4). The m7G MeRIP-sequencing protocol was performed using methylated sites on RNAs (m7G peaks) identified by MACS software. Differentially methylated sites were identified by diffReps and m7G peaks that overlapped with transcription exons were retained for further study. For lncRNA sequencing, the FPKM (fragments per kilobase of exon per million fragments mapped) value was obtained as the expression profiles of lncRNA using Cuffdiff (v2.2.1). Differentially expressed lncRNAs were identified by the fold change and p-value (fold change ≥2 and P ≤ 0.00001). In addition, gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis of the differentially-expressed RNAs and the differentially-methylated protein-coding genes were performed according to previous reports [40, 41]. P values are calculated by DAVID tool for GO and KEGG pathway analysis.
Data were presented as the mean ± standard deviation (SD). Statistical analyses were performed using SPSS 19.0 (Chicago, IL, USA) and GraphPad Prism 5.0 (La Jolla, CA). A two-tailed Student’s t-test or two-sided Wilcoxon-Mann-Whiteney test was used to determine the significant differences between the two groups. The number of samples in each group is shown as n in figure legends. Differences of P < 0.05 were defined as statistically significant. All experiments were independently repeated at least three times.
Availability of data and materials
The raw high-throughput m7G and lncRNAs sequencing generated during the current study can be obtained from the NCBI Gene Expression Omnibus (GEO): accession number GSE182042 and GSE182054, [https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE182042], token:[klsluyqezjkfdir].
Long non-coding RNAs
Hypoxic pulmonary hypertension
Chronic obstructive pulmonary disease
Pulmonary artery endothelial cells
pulmonary artery smooth muscle cells
Foetal bovine serum
Kyoto encyclopedia of genes and genomes
Wnt family member 5a
vomeronasal 2 receptor, 18
- Zfp512b, Zfp282,Zfp608, Zfp26:
zinc finger protein 512b, 282, 608 and 26
SMAD family member 3
Inositol monophosphatase domain containing 1
Adaptor related protein complex 2 subunit beta 1
Core-binding factor subunit beta
Breast cancer metastasis-suppressor 1-like
USP6 N-terminal like
Parathyroid hormone like hormone
Collagen type V alpha 3 chain
Adrenoceptor beta 3
Nuclear factor of activated T cells 5
Nanos C2HC-type zinc finger 2
Trinucleotide repeat containing adaptor 6B
Solute carrier family 30 member 2
MBL associated serine protease 1
Fukutin related protein
Ring finger protein 216
StAR related lipid transfer domain containing 8
HIC ZBTB transcriptional repressor 2
cyclin dependent kinase inhibitor 1B
Retinol dehydrogenase 13
AMME chromosomal region gene 1-like
Pogo transposable element derived with KRAB domain
Zinc finger and BTB domain containing 5
Single-minded family bHLH transcription factor 1
Glial fibrillary acidic protein
Tripartite motif containing 2
BLOC-1 related complex subunit 7
Adenosylhomocysteinase like 2
Platelet-activating factor receptor
Folliculin interacting protein 1
G protein-coupled receptor 155
A disintegrin-like and metallopeptidase (reprolysin type) with thrombospondin type 1 motif, 12
RAB11 family interacting protein 1
Transmembrane protein 198b
CREB/ATF bZIP transcription factor
lipoma HMGIC fusion partner-like 2
Ethanolamine kinase 1
Collagen type XXVI alpha 1 chain
Protein phosphatase 4 regulatory subunit 2
DEAD-box helicase 19B
Neuron navigator 2
Frost A, Badesch D, Gibbs J, et al. Diagnosis of pulmonary hypertension. Eur Respir J. 2019;53(1):1801904.
Simonneau G, Montani D, Celermajer DS, et al. Haemodynamic definitions and updated clinical classification of pulmonary hypertension. Eur Respir J. 2019;53(1):1801913.
Rowan SC, Keane MP, Gaine S, et al. Hypoxic pulmonary hypertension in chronic lung diseases: novel vasoconstrictor pathways. Lancet Respir Med. 2016;4:225–36.
Kylhammar D, Rådegran G. The principal pathways involved in the in vivo modulation of hypoxic pulmonary vasoconstriction, pulmonary arterial remodelling and pulmonary hypertension. Acta Physiol (Oxf). 2017;219:728–56.
Ransohoff JD, Wei Y, Khavari PA. The functions and unique features of long intergenic non-coding RNA. Nat Rev Mol Cell Biol. 2018;19:143–57.
Simion V, Haemmig S, Feinberg MW. LncRNAs in vascular biology and disease. Vasc Pharmacol. 2019;114:145–56.
Gonzalez I, Munita R, Agirre E, et al. AlncRNA regulates alternative splicing via establishment of a splicing-specific chromatin signature. Nat StructMol Biol. 2015;22:370–6.
Liu Y, Sun Z, Zhu J, et al. LncRNA-TCONS_00034812 in cell proliferation and apoptosis of pulmonary artery smooth muscle cells and its mechanism. J Cell Physiol. 2018;233:4801–14.
Derrien T, Johnson R, Bussotti G, et al. The GENCODE v7 catalog of human long noncoding RNAs: analysis of their gene structure, evolution, and expression. Genome Res. 2012;22:1775–89.
Fu XD. Non-coding RNA: a new frontier in regulatory biology. Natl Sci Rev. 2014;1:190–204.
Huang Y. The novel regulatory role of lncRNA-miRNA-mRNA axis in cardiovascular diseases. J Cell Mol Med. 2018;22:5768–75.
Freedman JE, Miano JM. Challenges and opportunities in linking long noncoding RNAs to cardiovascular, lung, and blood diseases. Arterioscler Thromb Vasc Biol. 2017;37:21–5.
Kopp F, Mendell JT. Functional classification and experimental dissection of long noncoding RNAs. Cell. 2018;172:393–407.
Han B, Bu P, Meng X, et al. Microarray profiling of long noncoding RNAs associated with idiopathic pulmonary arterial hypertension. Exp Ther Med. 2019;18:2353.
Zhang CY, Liu M, Wan JM, et al. Role of noncoding RNA in pulmonary arterial hypertension and potential drug therapeutic target. Curr Top Med Chem. 2018;18:975–86.
Wang S, Cao W, Gao S, et al. TUG1 regulates pulmonary arterial smooth muscle cell proliferation in pulmonary arterial hypertension. Can J Cardiol. 2019;35:1534–45.
Su H, Xu X, Yan C, et al. LncRNA H19 promotes the proliferation of pulmonary artery smooth muscle cells through AT(1)R via sponging let-7b in monocrotaline-induced pulmonary arterial hypertension. Respir Res. 2018;19:254.
He C. Grand challenge commentary: RNA epigenetics. Nat Chem Biol. 2010;6:863–5.
Jühling F, Mörl M, Hartmann RK, et al. tRNAdb 2009: compilation of tRNA sequences and tRNA genes. Nucleic Acids Res. 2009;37:D159–62.
Tomikawa C. 7-Methylguanosine modifications in transfer RNA (tRNA). Int J Mol Sci. 2018;19(12):4080.
Furuichi Y. Discovery of m(7)G-cap in eukaryotic mRNAs. Proc Jpn Acad Ser B Phys Biol Sci. 2015;91:394–409.
Su H, Wang G, Wu L, et al. Transcriptome-wide map of m(6)a circRNAs identified in a rat model of hypoxia mediated pulmonary hypertension. BMC Genomics. 2020;21:39.
Caruso P, MacLean MR, Khanin R, et al. Dynamic changes in lung microRNA profiles during the development of pulmonary hypertension due to chronic hypoxia and monocrotaline. Arterioscler Thromb Vasc Biol. 2010;30:716–23.
Zhang LS, Liu C, Ma H, et al. Transcriptome-wide mapping of internal N(7)-MethylguanosineMethylome in mammalian mRNA. Mol Cell. 2019;74:1304–16.e8.
Wu Y, Yang X, Chen Z, et al. m6A-induced lncRNA RP11 triggers the dissemination of colorectal cancer cells via upregulation of Zeb1. Mol Cancer. 2019;18:87.
Pandolfini L, Barbieri I, Bannister AJ, et al. METTL1 promotes let-7 MicroRNA processing via m7G methylation. Mol Cell. 2019;74:1278–90.e9.
Chu JM, Ye TT, Ma CJ, et al. Existence of internal N7-Methylguanosine modification in mRNA determined by differential enzyme treatment coupled with mass spectrometry analysis. ACS Chem Biol. 2018;13:3243–50.
Malbec L, Zhang T, Chen YS, et al. Dynamic methylome of internal mRNA N7-methylguanosine and its regulatory role in translation. Cell Res. 2019;29:927–41.
Blanco I, Santos S, Gea J, et al. Sildenafil to improve respiratory rehabilitation outcomes in COPD: a controlled trial. Eur Respir J. 2013;42:982–92.
Hoeper MM, Halank M, Wilkens H, et al. Riociguat for interstitial lung disease and pulmonary hypertension: a pilot trial. Eur Respir J. 2013;41:853–60.
Lee SJ, Zhang M, Hu K, et al. CCN1 suppresses pulmonary vascular smooth muscle contraction in response to hypoxia. Pulm Circ. 2015;5:716–22.
Braga EA, Fridman MV, Moscovtsev AA, et al. LncRNAs in ovarian Cancer progression, metastasis, and Main pathways: ceRNA and alternative mechanisms. Int J Mol Sci. 2020;21(22):8855.
Sanada TJ, Sun XQ, Happé C, et al. Altered TGFβ/SMAD signaling in human and rat models of pulmonary hypertension: an old target needs attention. Cells. 2021;10(1):84.
Sharmin N, Nganwuchu CC, Nasim MT. Targeting the TGF-β signaling pathway for resolution of pulmonary arterial hypertension. Trends Pharmacol Sci. 2021;42:510–3.
Frederick JP, Tafari AT, Wu SM, et al. A role for a lithium-inhibited Golgi nucleotidase in skeletal development and sulfation. Proc Natl Acad Sci U S A. 2008;105:11605–12.
Bajaj R, Kundu ST, Grzeskowiak CL, et al. IMPAD1 and KDELR2 drive invasion and metastasis by enhancing Golgi-mediated secretion. Oncogene. 2020;39:5979–94.
Yang YF, Wang YY, Hsiao M, et al. IMPAD1 functions as mitochondrial electron transport inhibitor that prevents ROS production and promotes lung cancer metastasis through the AMPK-Notch1-HEY1 pathway. Cancer Lett. 2020;485:27–37.
Kang K, Peng X, Zhang X, et al. MicroRNA-124 suppresses the transactivation of nuclear factor of activated T cells by targeting multiple genes and inhibits the proliferation of pulmonary artery smooth muscle cells. J Biol Chem. 2013;288:25414–27.
de Jesus PV, Yuan K, Alastalo TP, et al. Targeting the Wnt signaling pathways in pulmonary arterial hypertension. Drug Discov Today. 2014;19:1270–6.
Harris MA, Clark J, Ireland A, et al. The gene ontology (GO) database and informatics resource. Nucleic Acids Res. 2004;32(Database issue):D258–61.
Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545–51.
The authors gratefully appreciate Cloud-Seq Biotech Ltd. Co. (Shanghai, China) for the MeRIP-Seq service and the subsequent bioinformatics analysis.
This study was supported by the National Natural Science Foundation of China (Grant Numbers 81570043 and 82170063), the Science Foundation of Ministry of Health of Jiangsu Province in China (Grant Number ZDA2020016) and the 2021 Jiangsu Shuangchuang (Mass Innovation and Entrepreneurship) Talent Program (Grant Number JSSCBS20210147). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Ethics approval and consent to participate
The study protocol has been submitted to and approved by the Institutional Animal Care and Use Committee of Zhejiang University, China (permit number: SYXK 2017–0006). In addition, the study is reported in accordance with ARRIVE guidelines (https://arriveguidelines.org).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Wang, H., Chen, R.B., Zhang, S.N. et al. N7-methylguanosine modification of lncRNAs in a rat model of hypoxic pulmonary hypertension: a comprehensive analysis. BMC Genomics 23, 33 (2022). https://doi.org/10.1186/s12864-021-08188-8
- Hypoxic pulmonary hypertension
- N7-methylguanosine modification