- Research
- Open access
- Published:
The prevalence of pathogenic variants in the BMPR2 gene in patients with the idiopathic pulmonary arterial hypertension in the Russian population: sequencing data and meta-analysis
Respiratory Research volume 26, Article number: 146 (2025)
Abstract
Background
Idiopathic pulmonary arterial hypertension (IPAH) is a rare and severe form of pulmonary hypertension, with a genetic basis most commonly associated with mutations in the BMPR2 gene. However, no genetic testing has been reported for IPAH patients in the Russian population, nor have systematic studies been conducted to assess the frequency of pathogenic variants in this group.
Methods
The study cohort included 105 IPAH patients, consisting of 23 males and 82 females, who were managed at the PH care center in Moscow, Russia, from 2014 to 2024. Genetic testing was performed using whole-genome sequencing. Variant identification and annotation were conducted using GATK, DeepVariant, VEP, sv-callers and AnnotSV. A meta-analysis, performed with MOOSE, included 24 studies involving 3124 IPAH patients and 470 P/LP variants. Pathogenicity reassessment was carried out using InterVar, which incorporates ACMG criteria.
Results
Analysis of 105 adult IPAH patients in Russia revealed 11 patients (10.48%) as carriers of pathogenic or likely pathogenetic (P/LP) BMPR2 variants. As the result of reassessment, the number of P/LP BMPR2 variants raised from 394 (59%) to 445 (67%) with 80 pathogenic variants became of uncertain significance, and 152 unclassified variants became P/LP. The meta-analysis of these reevaluated pathogenic variants showed that while the frequency of P/LP variants in our cohort (10.48%) is lower than the overall average of 17.75% from the meta-analysis, the difference is not statistically significant (p = 0.062). Additionally, we report three P/LP BMPR2 variants, not reported in literature, with one being structural, and four P/LP variants in TBX4, ATP13A3 and AQP1 genes from 27 IPAH genes in 3 patients.
Conclusions
For the first time, we present the results of genetic testing in IPAH patients from the Russian population. Despite the considerable heterogeneity in the world-wide data, the prevalence of pathogenic BMPR2 mutations in IPAH patients from the Russian population does not significantly differ from the overall average in the meta-analysis. It is crucial to periodically reassess the pathogenicity of published variants, as half of the pathogenic BMPR2 IPAH variants were reclassified as LP or of uncertain significance.
Introduction
Pulmonary hypertension (PH) is a condition characterized by an increased mean pulmonary artery pressure (mPAP) of more than 20 mmHg, assessed by right heart catheterization at rest. One of the types of PH is idiopathic pulmonary arterial hypertension (IPAH), a rare and severe form with a genetic basis. IPAH is marked by pulmonary artery remodeling, which leads to increased pulmonary vascular resistance, right ventricular hypertrophy, and ultimately right heart failure.
If relatives have the same disease, it is referred to as “hereditary pulmonary arterial hypertension” (HPAH) [1]. The inheritance of the disease is autosomal dominant with low penetrance, which varies between males and females. Therefore, the term “IPAH/HPAH” is commonly used to refer to both forms, as they share the same clinical characteristics. Currently, the disease mechanism is associated with defects in genes encoding components of the bone morphogenetic protein receptor type-2 (BMPR2) signaling pathway. To date, 27 genes have been described as potentially involved in the disease [2]. Based on experimental and clinical evidence, the involvement of 12 of these genes is considered definitive, while 6 others have moderate evidence. The remaining genes have either low or disputed evidence.
The most common IPAH-associated pathogenic variants are located in BMPR2 gene. Its frequency varies from study to study in a wide range from 8% [2] to 52% [3], in the latter case, among lung transplant IPAH/HPAH patients [3].
These differences could arise from several factors. They may be related to variations in the severity of IPAH/HPAH within specific groups, reflecting the distinct characteristics of the study populations. Additionally, the frequency of pathogenic variant carriers may be influenced by ethnic or regional factors. Furthermore, we suggest that the use of different methods and programs to assess the pathogenicity of variants in IPAH/HPAH patients could significantly contribute to these discrepancies.
Thus, the aim of the study was to determine the carriage of pathogenic variants of IPAH/HPAH–associated genes, to re-evaluate all published variants of the BMPR2 gene using a modern method considering the current pathogenicity criteria, and to compare the frequency of carriage of pathogenic variants of the BMPR2 gene in the Russian population with that in other studies using meta-analysis.
Methods
Study cohort
The study cohort consisted of 105 patients with IPAH, including 23 males and 82 females managed in the PH care center from 2014 to 2024. Most of the patients (80%) were from the Moscow region and 87%—ethnic Russians. The diagnosis of idiopathic pulmonary arterial hypertension was verified according to the current guidelines of the European Society of Cardiology for the diagnosis and treatment of pulmonary hypertension (ESC/ERS Guidelines for the Diagnosis and Treatment of Pulmonary Hypertension) [4], by excluding other groups of pulmonary hypertension (associated with left heart disease, lung diseases, chronic thromboembolic pulmonary hypertension) and associated conditions (congenital heart defects, systemic connective tissue diseases, HIV infection, etc.). Hemodynamic inclusion criteria included a mean pulmonary arterial pressure (mPAP) > 25 mmHg, pulmonary vascular resistance (PVR) > 3 Wood units, and a pulmonary artery wedge pressure (PAWP) < 15 mmHg. Following the publication of the ESC/ERS 2022 recommendations, the mPAP threshold was revised to > 20 mmHg, PVR > 2 Wood units, and PAWP < 15 mmHg. In cases of borderline PAWP values, a fluid challenge test was performed to exclude heart failure with preserved ejection fraction (HFpEF). Indeed, there is currently an observed increase in the average age of patients with idiopathic pulmonary arterial hypertension, which is a global trend likely associated with improved disease diagnosis. In the COMPERA registry, the average age of patients with idiopathic/hereditary/drug-induced pulmonary arterial hypertension was 65.9 ± 15.9 years [5].
Baseline characteristics of PAH patient data including age, age of IPAH first symptoms, age at diagnosis, time from IPAH first symptoms to diagnosis, six-minute walking distance, WHO functional class at time of diagnosis, peripheral oedema, mean pulmonary arterial pressure (mPAP), pulmonary vascular resistance (PVR), right atrial pressure (RAP), tricuspid annular plane systolic excursion (TAPSE), NT-proBNP, N-terminal pro-brain natriuretic peptide are presented in Table 1.
Sequencing
Genomic DNA was extracted using the magnetic bead-based sorption method (MGIEasy Magnetic Beads Blood Genomic DNA Extraction Kit, MGI) and subsequently used for the preparation of genomic libraries for sequencing. Library preparation was performed using a PCR-free protocol with enzymatic DNA fragmentation (MGIEasy FS PCR-Free Library Prep Set, 96 reactions (MIX), MGI). Libraries were sequenced using DNBSEQ-T7 (PE150) sequencer technologies following the manufacturer’s recommendations.
Variant calling and annotation
The pipeline for calling germline variants was implemented as follows. Initially, adapter sequences and low-quality nucleotides were trimmed with Trimmomatic [6]. Subsequently, the reads were aligned to the GRCh38 reference genome, including additional contigs, in accordance with GATK Best Practices, utilizing bwa version 0.7.17 [7]. Following alignment, duplicates, which comprise approximately 1% of the reads due to the PCR-free library preparation protocol, were marked using the GATK MarkDuplicatesSpark version 4.3.0.0 [8]. Finally, germline SNVs and InDels calling was performed with DeepVariant version 1.4.0 [9]. Only variants with “pass” filter were included in the analysis.
Annotation of the identified variants was performed using the VEP tool [10]. The classification of the pathogenicity of genetic variants was conducted in accordance with the recommendations of the American College of Medical Genetics and Genomics (ACMG), utilizing the InterVar tool [11]. After that we performed additional search for the discovered variants in the databases HGMD [12] and ClinVar [13].
The aligned bam files of 105 patients were used for structural variants (SVs) calling using sv-callers [14], a union of four callers: (1) Lumpy, (2) GRIDSS, (3) DELLY, and (4) Manta. Each tool detects SVs based on a combination of coverage, split reads, and assembly, and has unique discrete workflows. Briefly, SV detections were performed independently with each tool. SURVIVOR (v1.0.7) was then used to merge SVs longer than 50 bp, of the same type, within 100 bp of each other, and identify consensus calls as well as eligible SR > 5&PE > 3 metrics. Only SVs that overlapped with the gene set of interest were included in the analysis. Further analysis included pathogenicity annotation using AnnotSV, based on the joint consensus recommendation of ACMG and ClinGen.
Selection of genes for analysis
IPAH gene list included 27 genes that are divided into five category: definitive (BMPR2, ACVRL1, ATP13A3, CAV1, EIF2AK4, ENG, GDF2, KCNK3,KDR, SMAD9, SOX17, and TBX4), moderate (ABCC8, GGCX, and TET2), limited (AQP1, BMP10, FBLN2, KLF2, KLK1, and PDGFD), disputed (BMPR1A, BMPR1B, NOTCH3, SMAD1, and SMAD4), and unknown (TOPBP1) [2].
Selection of articles for meta-analysis
For meta-analysis, we selected studies that reported data on the frequency of pathogenic and likely pathogenic variants in the BMPR2 gene among patients with IPAH. A total of 561 articles were initially identified in the PubMed database using the following search query: “(idiopathic pulmonary arterial hypertension OR pulmonary arterial hypertension OR PAH OR IPAH) AND (BMPR2 OR bone morphogenetic protein receptor type 2).” After a preliminary screen based on abstracts, 144 articles were selected for further consideration. Following a comprehensive full-text review, 81 studies were included in the analysis. Of those, 48 articles contained relevant data on the frequency of BMPR2 variants, while detailed variant descriptions were provided in only 24 of the selected studies (Fig. 1).
A meta-analysis of the frequency of pathogenic variants was conducted in accordance with the MOOSE (Meta-Analysis of Observational Studies in Epidemiology) guidelines [15]. The Freeman-Tukey transformation was used to calculate the pooled weighted proportion under both the fixed and random effects models. Heterogeneity of the model was assessed using the Q statistic and the I2 index. In cases of substantial heterogeneity (I2 > 75%), the random effects model was applied, while the fixed effects model was used for lower heterogeneity levels. For the Q-statistic, a p < 0.1 was considered to indicate heterogeneity. The presence of publication bias was assessed using the Egger's test and Begg's test.
Results
Discovered pathogenic and likely pathogenic variants in IPAH/HPAH patients in the Russian population.
In the cohort of 105 patients with IPAH/HPAH nine SNVs and one structural deletion were identified in the BMPR2 gene (NM_001204.7) as pathogenic or likely pathogenic based on ACMG criteria. All variants are heterozygous. Seven of ten variants have been previously reported in dbSNP and literature. A subsequent analysis of these seven variants in HGMD and ClinVar confirmed their pathogenicity. All variants had a “Disease-causing mutation” status associated with “Pulmonary hypertension, primary” in the HGMD database and P or P/LP with 2-star review status for rs863223426, rs1060502581, rs137852751, rs1085307151 and rs137852746, and 1-star for rs946132834 and rs1085307325 status in ClinVar (Table 2 and Fig. 2).
Distribution of the BMPR2 gene variants over genomic regions. The distribution of unique variants from the literature after InterVar reassessment. Only pathogenic (146) and likely pathogenic (299) variants within the BMPR2 gene (NM_001204.7) are presented. Additionally, positions of 10 pathogenic and likely pathogenic variants detected in our cohort are indicated
Four P/LP variants were also identified in other genes from the IPAH gene list (Table 2). All variants are heterozygous. 2 in the gene TBX4 (NM_001321120.2), 1 in the gene ATP13A3 (NM_001367549.1) and 1 in the gene AQP1 (NM_198098.4). Three patients were identified as carriers of these variants.
The clinical data of patients with identified pathogenic or likely pathogenic variants in the 27 IPAH genes are presented in Table 3.
Notably, no significant differences were observed in the overall clinical profiles of these patients. However, individuals with these variants exhibited significantly lower cardiac output (3.08 ± 0.88 L/min, 4.9 ± 1.96 L/min, p = 0.002), cardiac index (1.8 ± 0.41 L/min/sq m vs 2.6 ± 1.04 L/min/sq m, p = 0.005), and higher PVR (21.0 ± 10.67 Wood unit Vs 12.0 ± 8.52 Wood unit, p = 0.002).
Reevaluation BMPR2 gene pathogenicity with ACMG criteria and InterVar software tool
Before conducting the meta-analysis to compare the frequencies of IPAH-associated genes in the Russian population with those in global populations, we first decided to reassess the pathogenicity of the published variants. The studies included in the analysis used heterogeneous methods for determining P/LP variants. To standardize the approach, we applied a single pathogenicity assessment tool, InterVar, to all studies included in the meta-analysis. The meta-analysis focused solely on the BMPR2 gene, as we did not detect a significant number of P/LP variants in other genes, likely due to the limited sample size.
The literature analysis identified 665 unique variants of the BMPR2 gene (Supplemental files 1, 2). For each variant we included pathogenicity evaluation as it was presented in the corresponding study, along with the variant type. The distribution of reassigned P/LP variants of BMPR2 gene over genomic regions is presented in Fig. 2. To overcome heterogeneity in methods for pathogenicity assessment in various studies, for all variants we utilized the same program InterVar, which is the widely used and most cited tool for pathogenicity assessment based on ACMG criteria. Out of 342, initially defined as pathogenic (without mixed classification) variants, only 87 (25.4%) remain pathogenic, 172 (50.3%) were reclassified as likely pathogenic and 80 (23.4%) as variants of uncertain significance. 3 (0.9%) pathogenic variants were even moved to likely benign category. All 230 previously unclassified variants were reclassified to all categories as pathogenic, likely pathogenic, uncertain significance, likely benign and benign. As the result of reassessment, we raised the number of P/LP BMPR2 variants from 394 (59%) to 445 (67%). The distribution of BMPR2 gene variants, identified from the literature, over categories before and after reassessment is given in Fig. 3.
Reassessment of pathogenicity of variants in the BMPR2 gene. The distribution of 665 unique variants in the BMPR2 gene with pathogenicity assessment presented by authors in the literature before and after our reassessment according to the ACMG criteria. Some variants had mixed classifications before re-evaluation, either because they were identified in multiple studies with different assessments or because authors within an individual study were unable to classify the variant conclusively. Some variants reported in the literature do not have a definitive pathogenicity classification and are thus labeled as “Unclassified.”
Meta-analysis of the prevalence of pathogenic variants of BMPR2 gene in IPAH/HPAH patients
The final meta-analysis included 24 studies that evaluated the frequency of pathogenic variants in the BMPR2 gene among patients with IPAH/HPAH (full list with summary is presented in Table 4). Additionally, the analysis included the results from our sequencing data of 105 IPAH/HPAH patients from the Russian population with 11 patients being identified as carriers of pathogenic or likely pathogenic variants.
The final analysis included data from 3124 patients with IPAH/HPAH. Due to significant heterogeneity in the data, as indicated by the Q statistic of 100.4 and I2 value of 76.10%, a random-effects model was applied for the final frequency estimation. After reassessment the number of pathogenic variants’ carriers decreased in 19 studies, increased in one study [29], and not changed in 4 studies. The data heterogeneity index (I2 = 69.02%) allowed to use the fixed effects model. The proportion of patients whose carrier frequency of pathogenic variants was revised comprised 21.9% (18.97–25.15). No systematic error in data selection was recorded.
All analyzed studies fell within the confidence interval of the funnel plot, although asymmetry in the frequency distribution was observed (Fig. 4). In our cohort, the frequency of pathogenic or likely pathogenic variants was [10.48% (5.34–17.97)]. Overall average frequency from the meta-analysis was [17.75%, (14.88–20.72)], though the difference was not statistically significant (p = 0.062). No evidence of selection bias was detected based on Egger’s test (p = 0.661) and Begg’s test (p = 0.851). Here we conclude that the frequency of BMPR2 gene variants in the Russian population is not statistically different from those reported for other populations.
Meta-analysis of the frequency of pathogenic BMPR2 gene variants in patients with IPAH/HPAH. A Funnel plot of the meta-analysis of the frequency of pathogenic BMPR2 gene variants in patients with IPAH/HPAH. B Meta-analysis data of the frequency of pathogenic BMPR2 gene variants in patients with IPAH/HPAH
Discussion
Thus, we provide the first report on the prevalence of P/LP variants in genes associated with IPAH/HPAH in the Russian population. The patients included in the study were managed at Moscow PAH expert center, and the ethnic composition of the cohort was typical of the Russian population. The prevalence of PAH varies widely, ranging from 8 to 50 cases per 1 million people, with approximately one-third of these cases being IPAH patients [38]. The estimated number of IPAH/HPAH patients in Moscow, with a population of about 13 million, ranges from 35 to 217. Therefore, we believe that our study sample of 105 patients with IPAH/HPAH, 80% of whom are local residents, is representative. The clinical characteristics of the genotyped patients were comparable with the data of observational studies—ASPAIRE [39], HOPE [40], USPHSR [41] by the age of the included patients, the proportion of women, and the distribution by functional classes.
The BMPR2 is the primary gene associated with the development of pulmonary arterial hypertension. The presence of pathogenic variants in the BMPR2 among patients with PAH is typically associated with an earlier onset of the disease and more pronounced hemodynamics and respiratory impairments [42]. The data of our study partially coincides with the data of this meta-analysis. Patients with P/LP variants of BMPR2 had higher pulmonary resistance and lower cardiac output and cardiac index. The absence of differences in the age of IPAH diagnosis, pulmonary artery pressure level and some other parameters may be due to the relatively small cohort in our study.
Currently, genetic testing and counseling are indicated for patients with an established diagnosis of IPAH/HPAH [4]. Relatives of the proband with P/LP variants of BMPR2 have a higher risk of the disease. According to the DELPHI-2 study data, if asymptomatic carriers of BMPR2 gene variants are identified, an annual comprehensive examination for at least two years may be recommended. The authors of this study suggest that this approach can help facilitate earlier disease diagnosis and the timely initiation of treatment [43].
However, the significance of the BMPR2 gene is not limited to direct causality in relation to IPAH. The literature suggests that IPAH is oligogenic, meaning it is influenced by multiple genetic factors rather than being caused by a single gene mutation [44, 45]. Carrying multiple genetic variants can influence the progression of the disease, which is why it is crucial to understand the frequency of pathogenic variant carriers within each population. The analysis of sequencing data for the BMPR2 gene revealed that the overall prevalence of pathogenic and likely pathogenic variants in our cohort was 10.48% (11 out of 105 patients). Of these, only 7 out of the 10 variants had been previously described in the literature.
Three variants—frameshift variant c.793_794delGA, splice acceptor variant c.419-1G > A and structural deletion g.202514389_202517603del – have not been previously reported in literature (Table 2). Two nucleotide deletion c.793_794delGA induces a frameshift in the region of 6th exon, affecting key positions essential for protein function. Exon 6 is situated in the N-terminal part of the serine-threonine kinase domain, which is made up of conserved subdomains, including exons 6 through 11. This region plays a key role in ATP-binding and contains specific patterns of conserved amino acids. Mutations within this area can disrupt signaling by interfering with ATP-binding and affect the protein function [46]. There is no available information on frequency of the c.793_794delGA variant, therefore, the InterVar defined this variant as likely pathogenic. The splice acceptor variant c.419-1G > A is located within a cryptic splice site in the third intron with a splice distance of − 1. Splicing prediction tools, including SpliceAI and MaxEntScan, indicate a significant impact of this variant on splicing, suggesting complete disruption of the existing acceptor site. Additional in silico tools, such as CADD, DANN, and FATHMM-MKL, also classify this variant as pathogenic. There is no information on the frequency of the variant. An alternative substitution (T instead of A in our study) at the same position (c.419-1G > T) was reported in the study of Zhu et al. [47] in a patient with PAH. This variant (c.419-1G > T) is recorded in the HGMD as a disease-causing mutation for pulmonary hypertension, primary. According to the combination of factors, the InterVar classified the variant c.419-1G > A as pathogenic. The heterozygous structural variant (deletion) of 3214 bp in length in the BMPR2 gene (NC_000002.12:g.202514389_202517603del) has not been previously reported and, based on AnnotSV predictions, is classified as pathogenic. This deletion completely deletes the 5th exon of the gene and significant regions of surrounding introns. Moreover, this deletion overlaps with the previously described deletion in the dbVar database, identified in a patient with PH: nsv4682448 (NC_000002.12:g.202514864_202514989del).
We also identified three additional carriers of four pathogenic variants in the ATP13A3, AQP1, and TBX4 genes in our group, with one patient carrying two P/LP variants. However, due to the small sample size, it is not possible to draw reliable conclusions about the frequencies of these variants within the Russian population. Future studies will help refine the population frequencies of pathogenic variants in these genes and other IPAH-associated genes.
Reevaluation of the BMPR2 gene pathogenicity using ACMG criteria and InterVar software demonstrated that comparing the frequency of pathogenic variant carriers across different research groups requires a standardized approach to pathogenicity assessment. It is evident that the accumulation of data on the molecular genetic mechanisms of disease development will, in the future, clarify the functional significance of the identified variants. This underscores the importance of regularly reevaluating pathogenicity. To establish an irrefutable association between disease development and the gene, it must be reassessed at least three years after being assigned evidence level as “definitive” [48].
The most extensive data in the literature regarding the frequency of pathogenic variants is available for the BMPR2 gene. We attempted to compare the data from our cohort with findings from other studies; however, the reported frequencies of pathogenic variants in the BMPR2 gene vary significantly across different studies.
Differences in the frequency of detecting P/LP variant carriers may be linked to variations in the clinical characteristics of the cohorts analyzed. A meta-analysis was conducted to aggregate the literature data. Additionally, differences in frequency could be attributed to the methods used to assess the pathogenicity of the identified variants. To evaluate the pathogenicity of the variants, we applied the criteria established by the American College of Medical Genetics and Genomics (ACMG) in collaboration with the Association for Molecular Pathology (AMP) in 2015. These guidelines incorporate 28 distinct factors that facilitate the classification of genetic variants as “pathogenic,” “likely pathogenic,” “variants of uncertain significance (VUS),” “likely benign,” or “benign.” These criteria consider various parameters, including the variant's impact on protein function, allele frequency within the population, functional studies, the variant's position in relation to other known pathogenic variants, and other relevant factors [49].
Currently, there is no universally accepted methodology for the use of specific bioinformatics tools to evaluate the pathogenicity of genetic variants. The challenge of calibrating analytical tools remains highly relevant [50]. One of the most cited and widely used programs for interpretation of variant pathogenicity is InterVar, which automates the classification process based on 18 out of 28 criteria defined by the ACMG. InterVar usage significantly reduces analysis time and minimizes variability in interpretation among different experts [51,52,53,54,55,56,57,58]. Nevertheless, it should be noted that the interpretation of sequencing data using InterVar in patients with pulmonary hypertension has limitations, the main one being the low penetrance of the disease. The authors of this software product suggest that InterVar is most effective for interpreting variants associated with severe congenital or early developmental disorders with nearly 100% penetrance, and its performance may be less reliable for late-onset diseases or recessive conditions [11].
A meta-analysis of the literature, with variants reassessed according to ACMG criteria, indicated a slightly higher average BMPR2 mutation frequency of 17.75% (14.88–20.82). Despite the heterogeneity observed in the global data, our meta-analysis showed that the prevalence of BMPR2 mutations in the Russian population does not significantly differ from the global average.
Conclusion
The prevalence of pathogenic variants in the BMPR2 gene within our group was found to be 10.48%, which does not statistically differ from the 17.75% carrier frequency reported in the meta-analysis. It is essential to consider the revised pathogenicity criteria established since 2015 and to utilize state-of-the-art tools for reevaluating the data.
Study limitations
This study had several limitations. Firstly, InterVar uses only part of the ACMG criteria for pathogenicity assessment. The authors of the studies included in the meta-analysis did not publish the complete set of identified variants, including benign and likely benign variants. As a result, we only reclassified P, LP, and VUS variants, as well as those for which the authors did not provide a classification (N/A). A limitation of this approach is that pathogenic variants may exist among the benign and likely benign variants.
Data availability
The sequencing data of 105 IPAH patients that support the findings of this study are not openly available due to reasons of sensitivity and are available from the corresponding author upon reasonable request. Data are located in controlled access data storage at the Bioinformatics Lab at HSE University.
Abbreviations
- ACMG:
-
American College of Medical Genetics and Genomics
- AMP:
-
Association for Molecular Pathology
- BMPR2:
-
Bone morphogenetic protein receptor type-2
- CTEPH:
-
Chronic thromboembolic pulmonary hypertension
- HPAH:
-
Hereditary pulmonary arterial hypertension
- InDel:
-
Insertion and deletion
- IPAH:
-
Idiopathic pulmonary arterial hypertension
- LP variant:
-
Likely pathogenic variant
- mPAP:
-
Mean pulmonary arterial pressure
- MOOSE:
-
Meta-analysis of observational studies in epidemiology
- NT-proBNP:
-
N-terminal prohormone of Brain Natriuretic Peptide
- P variant:
-
Pathogenic variant
- PH:
-
Pulmonary hypertension
- PAH:
-
Pulmonary arterial hypertension
- Pts:
-
Patients
- PVR:
-
Pulmonary vascular resistance
- PWP:
-
Pulmonary wedge pressure
- RAP:
-
Right atrial pressure
- SNV:
-
Single nucleotide variant
- SV:
-
Structural variant
- TAPSE:
-
Tricuspid annular plane systolic excursion
- VUS:
-
Variants of uncertain significance
- WU:
-
Wood units
References
Aldalaan AM, Ramzan K, Saleemi SA, Weheba I, Alquait L, Abdelsayed A, Alzubi F, Zaytoun H, Alharbi N, Al-Owain M, Imtiaz F. Genetic basis of pulmonary arterial hypertension: a prospective study from a highly inbred population. Pulm Circ. 2021;11:20458940211032056.
Welch CL, Aldred MA, Balachandar S, Dooijes D, Eichstaedt CA, Gräf S, Houweling AC, Machado RD, Pandya D, Prapa M. Defining the clinical validity of genes reported to cause pulmonary arterial hypertension. Genet Med. 2023;25:100925.
Ghigna MR, Guignabert C, Montani D, Girerd B, Jais X, Savale L, Herve P, de Montpreville VT, Mercier O, Sitbon O, et al. BMPR2 mutation status influences bronchial vascular changes in pulmonary arterial hypertension. Eur Respir J. 2016;48:1668–81.
Humbert M, Kovacs G, Hoeper MM, Badagliacca R, Berger RMF, Brida M, Carlsen J, Coats AJS, Escribano-Subias P, Ferrari P, et al. 2022 ESC/ERS Guidelines for the diagnosis and treatment of pulmonary hypertension. Eur Heart J. 2022;43:3618–731.
Hoeper MM, Pausch C, Grünig E, Staehler G, Huscher D, Pittrow D, Olsson KM, Vizza CD, Gall H, Distler O. Temporal trends in pulmonary arterial hypertension: results from the COMPERA registry. Eur Respir J. 2022;59:2102024.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754–60.
Picard toolkit. http://broadinstitute.github.io/picard. Accessed 8 Apr 2025.
Poplin R, Chang PC, Alexander D, Schwartz S, Colthurst T, Ku A, Newburger D, Dijamco J, Nguyen N, Afshar PT, et al. A universal SNP and small-indel variant caller using deep neural networks. Nat Biotechnol. 2018;36:983–7.
McLaren W, Gil L, Hunt SE, Riat HS, Ritchie GR, Thormann A, Flicek P, Cunningham F. The ensembl variant effect predictor. Genome Biol. 2016;17:122.
Li Q, Wang K. InterVar: clinical interpretation of genetic variants by the 2015 ACMG-AMP guidelines. Am J Hum Genet. 2017;100:267–80.
Stenson PD, Ball EV, Mort M, Phillips AD, Shiel JA, Thomas NS, Abeysinghe S, Krawczak M, Cooper DN. Human gene mutation database (HGMD): 2003 update. Hum Mutat. 2003;21:577–81.
Landrum MJ, Lee JM, Riley GR, Jang W, Rubinstein WS, Church DM, Maglott DR. ClinVar: public archive of relationships among sequence variation and human phenotype. Nucleic Acids Res. 2014;42:D980-985.
Kuzniar A, Maassen J, Verhoeven S, Santuari L, Shneider C, Kloosterman WP, de Ridder J. sv-callers: a highly portable parallel workflow for structural variant detection in whole-genome sequence data. PeerJ. 2020;8: e8214.
Brooke BS, Schwartz TA, Pawlik TM. MOOSE reporting guidelines for meta-analyses of observational studies. JAMA Surg. 2021;156:787–8.
Wang MT, Weng KP, Chang SK, Huang WC, Chen LW. Hemodynamic and clinical profiles of pulmonary arterial hypertension patients with GDF2 and BMPR2 variants. Int J Mol Sci. 2024;25:2734.
Liang KW, Chang SK, Chen YW, Lin WW, Tsai WJ, Wang KY. Whole exome sequencing of patients with heritable and idiopathic pulmonary arterial hypertension in Central Taiwan. Front Cardiovasc Med. 2022;9:911649.
Eichstaedt CA, Sassmannshausen Z, Shaukat M, Cao D, Xanthouli P, Gall H, Sommer N, Ghofrani HA, Seyfarth HJ, Lerche M, et al. Gene panel diagnostics reveals new pathogenic variants in pulmonary arterial hypertension. Respir Res. 2022;23:74.
van den Heuvel LM, Jansen SMA, Alsters SIM, Post MC, van der Smagt JJ, Handoko-De Man FS, van Tintelen JP, Gille H, Christiaans I, Vonk Noordegraaf A, et al. Genetic evaluation in a cohort of 126 Dutch pulmonary arterial hypertension Patients. Genes. 2020;11:1191.
Wang XJ, Lian TY, Jiang X, Liu SF, Li SQ, Jiang R, Wu WH, Ye J, Cheng CY, Du Y, et al. Germline BMP9 mutation causes idiopathic pulmonary arterial hypertension. Eur Respir J. 2019;53:1801609.
Yang H, Zeng Q, Ma Y, Liu B, Chen Q, Li W, Xiong C, Zhou Z. Genetic analyses in a cohort of 191 pulmonary arterial hypertension patients. Respir Res. 2018;19:87.
Abou Hassan OK, Haidar W, Nemer G, Skouri H, Haddad F, BouAkl I. Clinical and genetic characteristics of pulmonary arterial hypertension in Lebanon. BMC Med Genet. 2018;19:89.
Navas P, Tenorio J, Quezada CA, Barrios E, Gordo G, Arias P, Lopez Meseguer M, Santos-Lozano A, Palomino Doza J, Lapunzina P, Escribano Subias P. Molecular analysis of BMPR2, TBX4, and KCNK3 and genotype-phenotype correlations in Spanish patients and families with idiopathic and hereditary pulmonary arterial hypertension. Rev Esp Cardiol. 2016;69:1011–9.
Best DH, Sumner KL, Smith BP, Damjanovich-Colmenares K, Nakayama I, Brown LM, Ha Y, Paul E, Morris A, Jama MA, et al. EIF2AK4 mutations in patients diagnosed with pulmonary arterial hypertension. Chest. 2017;151:821–8.
Pousada G, Baloira A, Valverde D. Complex inheritance in pulmonary arterial hypertension patients with several mutations. Sci Rep. 2016;6:33570.
Zhao Q, Zhang R, Shi J, Xie H, Zhang L, Li F, Jiang R, Wu W, Luo C, Qiu H, et al. Imaging features in BMPR2 mutation-associated pulmonary arterial hypertension. Radiology. 2023;307: e222488.
Zhang X, Zhang C, Li Q, Piao C, Zhang H, Gu H. Clinical characteristics and prognosis analysis of idiopathic and hereditary pulmonary hypertension patients with ACVRL1 gene mutations. Pulm Circ. 2021;11:20458940211044576.
Jang AY, Kim BG, Kwon S, Seo J, Kim HK, Chang HJ, Chang SA, Cho GY, Rhee SJ, Jung HO, et al. Prevalence and clinical features of bone morphogenetic protein receptor type 2 mutation in Korean idiopathic pulmonary arterial hypertension patients: the PILGRIM explorative cohort. PLoS ONE. 2020;15: e0238698.
van der Bruggen CE, Happe CM, Dorfmuller P, Trip P, Spruijt OA, Rol N, Hoevenaars FP, Houweling AC, Girerd B, Marcus JT, et al. Bone morphogenetic protein receptor type 2 mutation in pulmonary arterial hypertension: a view on the right ventricle. Circulation. 2016;133:1747–60.
Mutlu Z, Kayikcioglu M, Nalbantgil S, Vuran O, Kemal H, Mogulkoc N, Erturk B, Onay H, Eroglu Z, Kultursay H. Sequencing of mutations in the serine/threonine kinase domain of the bone morphogenetic protein receptor type 2 gene causing pulmonary arterial hypertension. Anatol J Cardiol. 2016;16:491–6.
Zeng Q, Yang H, Liu B, Ma Y, Liu Z, Chen Q, Li W, Luo Q, Zhao Z, Zhou Z, Xiong C. Clinical characteristics and survival of Chinese patients diagnosed with pulmonary arterial hypertension who carry BMPR2 or EIF2KAK4 variants. BMC Pulm Med. 2020;20:150.
Song J, Eichstaedt CA, Viales RR, Benjamin N, Harutyunova S, Fischer C, Grunig E, Hinderhofer K. Identification of genetic defects in pulmonary arterial hypertension by a new gene panel diagnostic tool. Clin Sci. 2016;130:2043–52.
Kataoka M, Aimi Y, Yanagisawa R, Ono M, Oka A, Fukuda K, Yoshino H, Satoh T, Gamou S. Alu-mediated nonallelic homologous and nonhomologous recombination in the BMPR2 gene in heritable pulmonary arterial hypertension. Genet Med. 2013;15:941–7.
Kabata H, Satoh T, Kataoka M, Tamura Y, Ono T, Yamamoto M, Huqun, Hagiwara K, Fukuda K, Betsuyaku T, Asano K. Bone morphogenetic protein receptor type 2 mutations, clinical phenotypes and outcomes of Japanese patients with sporadic or familial pulmonary hypertension. Respirology. 2013;18:1076–82.
Girerd B, Montani D, Jais X, Eyries M, Yaici A, Sztrymf B, Savale L, Parent F, Coulet F, Godinas L, et al. Genetic counselling in a national referral centre for pulmonary hypertension. Eur Respir J. 2016;47:541–52.
Gamou S, Kataoka M, Aimi Y, Chiba T, Momose Y, Isobe S, Hirayama T, Yoshino H, Fukuda K, Satoh T. Genetics in pulmonary arterial hypertension in a large homogeneous Japanese population. Clin Genet. 2018;94:70–80.
Eyries M, Montani D, Nadaud S, Girerd B, Levy M, Bourdin A, Tresorier R, Chaouat A, Cottin V, Sanfiorenzo C, et al. Widening the landscape of heritable pulmonary hypertension mutations in paediatric and adult cases. Eur Respir J. 2019;53:1801371.
Saleh K, Khan N, Dougherty K, Bodi G, Michalickova M, Mohammed S, Kerenidi T, Sadik Z, Mallat J, Farha S, Sabbour H. The first pulmonary hypertension registry in the United Arab Emirates (UAEPH): clinical characteristics, hemodynamic parameters with focus on treatment and outcomes for patients with group 1-PH. J Clin Med. 2023;12:1996.
Hurdman J, Condliffe R, Elliot CA, Davies C, Hill C, Wild JM, Capener D, Sephton P, Hamilton N, Armstrong IJ, et al. ASPIRE registry: assessing the Spectrum of Pulmonary hypertension Identified at a REferral centre. Eur Respir J. 2012;39:945–55.
Arvanitaki A, Vrana E, Boutsikou M, Anthi A, Apostolopoulou S, Avgeropoulou A, Demerouti E, Patrianakos A, Karyofyllis P, Mitrouska I, et al. The impact of cardiovascular comorbidities associated with risk for left heart disease on idiopathic pulmonary arterial hypertension: data from the Hellenic Pulmonary Hypertension Registry (HOPE). Pulm Circ. 2022;12: e12086.
Badlam JB, Badesch DB, Austin ED, Benza RL, Chung WK, Farber HW, Feldkircher K, Frost AE, Poms AD, Lutz KA, et al. United States pulmonary hypertension scientific registry: baseline characteristics. Chest. 2021;159:311–27.
Evans JD, Girerd B, Montani D, Wang XJ, Galie N, Austin ED, Elliott G, Asano K, Grunig E, Yan Y, et al. BMPR2 mutations and survival in pulmonary arterial hypertension: an individual participant data meta-analysis. Lancet Respir Med. 2016;4:129–37.
Montani D, Girerd B, Jais X, Laveneziana P, Lau EMT, Bouchachi A, Hascoet S, Gunther S, Godinas L, Parent F, et al. Screening for pulmonary arterial hypertension in adults carrying a BMPR2 mutation. Eur Respir J. 2021;58:2004229.
Loyd JE. Pulmonary arterial hypertension: insights from genetic studies. Proc Am Thorac Soc. 2011;8:154–7.
Eichstaedt CA, Belge C, Chung WK, Graf S, Grunig E, Montani D, Quarck R, Tenorio-Castano JA, Soubrier F, Trembath RC, et al. Genetic counselling and testing in pulmonary arterial hypertension: a consensus statement on behalf of the International Consortium for Genetic Studies in PAH. Eur Respir J. 2023;61:2201471.
Pousada G, Baloira A, Vilarino C, Cifrian JM, Valverde D. Novel mutations in BMPR2, ACVRL1 and KCNA5 genes and hemodynamic parameters in patients with pulmonary arterial hypertension. PLoS ONE. 2014;9: e100261.
Zhu N, Pauciulo MW, Welch CL, Lutz KA, Coleman AW, Gonzaga-Jauregui C, Wang J, Grimes JM, Martin LJ, He H, et al. Novel risk genes and mechanisms implicated by exome sequencing of 2572 individuals with pulmonary arterial hypertension. Genome Med. 2019;11:69.
Strande NT, Riggs ER, Buchanan AH, Ceyhan-Birsoy O, DiStefano M, Dwight SS, Goldstein J, Ghosh R, Seifert BA, Sneddon TP. Evaluating the clinical validity of gene-disease associations: an evidence-based framework developed by the clinical genome resource. Am J Hum Genet. 2017;100:895–906.
Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, Grody WW, Hegde M, Lyon E, Spector E, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17:405–24.
Bergquist T, Stenton SL, Nadeau EAW, Byrne AB, Greenblatt MS, Harrison SM, Tavtigian SV, O'Donnell-Luria A, Biesecker LG, Radivojac P, et al. Calibration of additional computational tools expands ClinGen recommendation options for variant classification with PP3/BP4 criteria. bioRxiv. 2024.
Lesurf R, Said A, Akinrinade O, Breckpot J, Delfosse K, Liu T, Yao R, Persad G, McKenna F, Noche RR, et al. Whole genome sequencing delineates regulatory, copy number, and cryptic splice variants in early onset cardiomyopathy. NPJ Genom Med. 2022;7:18.
Yuan J, Zhuang YY, Liu X, Zhang Y, Li K, Chen ZJ, Li D, Chen H, Liang J, Yao Y, et al. Exome-wide association study identifies KDELR3 mutations in extreme myopia. Nat Commun. 2024;15:6703.
Zech M, Winkelmann J. Next-generation sequencing and bioinformatics in rare movement disorders. Nat Rev Neurol. 2024;20:114–26.
Gamirova R, Shagimardanova E, Sato T, Kannon T, Gamirova R, Tajima A. Identification of potential disease-associated variants in idiopathic generalized epilepsy using targeted sequencing. J Hum Genet. 2024;69:59–67.
Deyell RJ, Shen Y, Titmuss E, Dixon K, Williamson LM, Pleasance E, Nelson JMT, Abbasi S, Krzywinski M, Armstrong L, et al. Whole genome and transcriptome integrated analyses guide clinical care of pediatric poor prognosis cancers. Nat Commun. 2024;15:4165.
Bai J, Shi J, Li C, Wang S, Zhang T, Hua X, Zhu B, Koka H, Wu HH, Song L, et al. Whole genome sequencing of skull-base chordoma reveals genomic alterations associated with recurrence and chordoma-specific survival. Nat Commun. 2021;12:757.
Jiang Y, Wang J, Sun M, Zuo D, Wang H, Shen J, Jiang W, Mu H, Ma X, Yin F, et al. Multi-omics analysis identifies osteosarcoma subtypes with distinct prognosis indicating stratified treatment. Nat Commun. 2022;13:7207.
Pleasance E, Titmuss E, Williamson L, Kwan H, Culibrk L, Zhao EY, Dixon K, Fan K, Bowlby R, Jones MR, et al. Pan-cancer analysis of advanced patient tumors reveals interactions between therapy and genomic landscapes. Nat Cancer. 2020;1:452–68.
Funding
The whole-genome sequencing data were obtained within the frames of the scientific agreement between the company Biotech Campus Ltd, HSE University and City Clinical Hospital No 29 by N.E. Bauman, Moscow, Russia. GO, VZ and MP were supported by the Basic Research Program at HSE University and the funding of Cardiogenetic group at the Faculty of Computer Science at HSE University.
Author information
Authors and Affiliations
Contributions
ED, EZ, AL, AD, NO performed patients selection and inclusion. GO, IB, VZ, MP performed bioinformatics analysis. IB, DN performed literature analysis. LP performed meta-analysis. DZ and MP supervised the study. GO, IB, VZ, DN, ED, LM, DZ and MP analysed the data and wrote the manusucript. All authors read and approved the manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
This study was conducted in accordance with the Declaration of Helsinki, approved by the Ethics Committee of City Clinical Hospital Nu 51, Moscow, Russia (protocol 19/15, 22 Oct 2014), and confirmed for extension (protocol #03/23, 4 Sept 2023) by the local Ethics Committee of the City Clinical Hospital Nu 29, Moscow, Russia.
Consent for publication
Non-applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
About this article
Cite this article
Okhrimenko, G., Borovikova, I., Dankovtseva, E. et al. The prevalence of pathogenic variants in the BMPR2 gene in patients with the idiopathic pulmonary arterial hypertension in the Russian population: sequencing data and meta-analysis. Respir Res 26, 146 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12931-025-03214-9
Received:
Accepted:
Published:
DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12931-025-03214-9