Skip to main content

Mutational and clinical analysis of the ENG gene in patients with pulmonary arterial hypertension

Abstract

Background

Pulmonary arterial hypertension (PAH) is a rare vascular disorder characterized by a capillary wedge pressure ≤ 15 mmHg and a mean pulmonary arterial pressure ≥ 25 mmHg at rest. PAH can be idiopathic, heritable or associated with other conditions. The aim of this study was to analyze the Endoglin (ENG) gene and assess the influence of the c.572G > A (p.G191D) mutation in patients with idiopathic or associated PAH. The correlation between the pathogenic mutations and clinical and functional parameters was further analyzed.

Results

Sixteen different changes in the ENG gene were found in 44 out of 57 patients. After in silico analysis, we classified eight mutations as pathogenic in 16 of patients. The c.572G>A (p.G191D) variation was observed in ten patients, and the analysis for the splicing process using hybrid minigenes, with pSPL3 vector to assess splicing alterations, do not generate a new transcript. Age at diagnosis (p = 0.049) and the 6-min walking test (p = 0.041) exhibited statistically significant differences between carriers and non-carriers of pathogenic mutations. Patients with pathogenic mutations exhibited disease symptoms 8 years before non-carriers. Five patients with pathogenic mutations were carriers of another mutation in the BMPR2 or ACVRL1 genes.

Conclusions

We present a series of PAH patients with mutations in the ENG gene, some of them not previously described, exhibiting clinical and hemodynamic alterations suggesting that the presence of these mutations may be associated with the severity of the disease. Moreover, genetic analysis in patients with PAH may be of clinical relevance and indicates the complexity of the genetic background.

The role of ENG gene in patients with Pulmonary Arterial Hypertension.

Background

Pulmonary arterial hypertension (PAH; OMIM #178600; ORPHA 422) is a severe disease affecting small pulmonary arteries that results in progressive remodeling leading to elevated pulmonary vascular resistance and right ventricular failure [1]. PAH can be idiopathic (IPAH), heritable (HPAH) or associated with other conditions (APAH) [2]. PAH is characterized by a capillary wedge pressure ≤ 15 mmHg and mean pulmonary arterial pressure ≥ 25 mmHg at rest [1, 2]. Symptoms of PAH are mixed but include dyspnea, syncope and chest pain. Eventually, the disease in these patients leads to right-sided heart failure and death [1]. The main pathologic changes associated with increased pulmonary vascular resistance are thrombus development, thickened intima, proliferation of smooth muscles cells, and growth of plexiform lesions in pulmonary vessels [3]. The estimated incidence is approximately 2–5 cases per million per year [3], and the gender ratio is 1.7:1 female vs male [4, 5]. Without treatment, the disease progresses to right ventricular failure and death within 3 years of diagnosis [6].

Heterozygous germline mutations in the bone morphogenetic protein type 2 receptor (BMPR2; MIM #600799) have been identified in approximately 10 to 40 % of IPAH patients without a reported familial history of the disease and in over 80 % of patients with HPAH [4, 79]. PAH patients with BMPR2 mutations are reported to develop more severe disease, are less likely to respond to treatment and are diagnosed approximately 10 years earlier than non-carriers [10]. In a few PAH patients, mutations in other genes participating in the BMPR2 signaling pathway have been reported, including Endoglin, also known as CD105, (ENG; MIM #601284) [11]. ENG gene mutations are less common than BMPR2 gene mutations in patients with PAH. Accordingly, a more complicated genetic background has been proposed for PAH [7].

The ENG gene encodes a type I integral membrane glycoprotein receptor that is a member of the Transforming growth factor beta (TGF-β) signaling superfamily. This receptor is expressed on proliferating vascular endothelial cells and in other cell types associated with cardiovascular system and controls diverse cellular processes, including cell differentiation, proliferation, angiogenesis, inflammation, and wound healings and has been linked to psoriatic skin, inflamed synovial arthritis, vascular injury, tumor vessels and apoptosis in embryonic and mature tissues [1215]. The human ENG gene is located on chromosome 9q33-34 [7, 13, 14], and the encoded protein exhibits an extracellular domain, hydrophobic transmembrane domain and a cytosolic domain. The extracellular domain contains 561 amino acids and is the largest of the domains [13]. This gene is implicated in hereditary hemorrhagic telangiectasia (HHT) type 1, an autosomal dominant syndrome characterized by vascular dysplasia. Mutations found in the ENG gene are an important factor for the development of HHT and may contribute to PAH in some HHT patients due to the gene’s function as a TGF-β receptor [7, 1316]. Mutations in this gene are frequently identified in HHT but are uncommon in PAH patients [4, 15, 17].

The potential role of ENG gene in patients with PAH remains unknown. To analyze its impact in patients with IPAH and APAH, we analyzed the coding region and intronic junctions of this gene and try to associate hemodynamic and clinical characteristics between carriers and non-carriers of ENG mutations. To evaluate the effect of ENG mutations on clinical outcomes of PAH, the phenotypical characteristic of carriers of missense mutations and carriers of mutations that alter the splicing in this gene were compared.

Methods

Patients and samples

As described previously [8], patients with idiopathic or associated PAH (group 1 of the classification of Nice) [18] treated in our clinic were included in this study. All patients are included in the CHUVI DNA Biobank (Biobanco del Complejo Hospitalario Universitario de Vigo). Patients signed an informed consent and the Autonomic Ethics Committee approved the study (Comité Autonómico de Ética da Investigación de Galicia-CAEI de Galicia).

In all cases, cardiac catheterization was performed using the latest consensus diagnostic criteria of the ERS-ESC (European Respiratory Society-European Society of Cardiology) [19]. PAH was considered idiopathic after the exclusion of any of the possible causes associated with the disease. Clinical histories included use of drugs, especially appetite suppressants, and screening for connective tissue diseases and hepatic disease. The study included serology for Human immunodeficiency virus (HIV), autoimmunity, thoracic tomography computerized scan (TC scan) and echocardiography. Patients with PAH that could be related to chronic lung disease were excluded. Fifty-five healthy individuals were used as controls.

Mutational analysis

Venous blood was collected from patients and healthy volunteers to extract genomic DNA using the FlexiGene DNA Kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol.

Amplification of the ENG gene was performed with 50 ng of genomic DNA from each patient and control. We amplified the exon regions and intronic junctions and did not analyzed changes in other regions for this study. The primers use to amplified this region by PCR (Polymerase chain reaction) were described by Gallione et al [20], with minor modifications (Table 1). The PCR mix was GoTaq® Green Master Mix (Promega Corporation, Madison, Wisconsin, USA), which contained Taq DNA polymerase, dNTPs, MgCl2 and reaction buffer. A second independent PCR and sequencing reaction in both the forward and reverse strands was performed to check for the detected mutations. PCR was performed in an MJ MiniTM Gradient Thermal Cycler (Bio-Rad, Hercules, California, USA). Electrophoresis on a 2 % agarose gel containing ethidium bromide was performed to confirmed PCR products in a Sub-Cell GT (Bio-Rad, Hercules, California, USA). HyperLadder IV-V (New England Biolabs, Ipswich, Massachusetts, USA) was used as the molecular weight marker. PCR fragments were purified using the ExoSAP-IT kit (USB Corporation, Cleveland, Ohio, USA) and sequenced with the BigDye Terminator version 3.1 Cycle Sequencing Kit (Applied Biosystems, Carlsbad, California, USA). The sequencing reactions were precipitated and analyzed on an ABI PRISM 3100 genetic analyzer (Applied Biosystems, Carlsbad, California, USA).

Table 1 Primers used to amplify the ENG gene

Sequence data were aligned to the reference Ensembl cDNA sequence ENSG00000106991 for the ENG gene and examined for sequence variations. To align and compare sequences in different organisms we use the Basic Local Alignment Search Tool (BLAST) software. Polyphen-2 (available at http://genetics.bwh.harvard.edu/pph/) characterize an amino acid substitution as “benign”, “possibly damaging” or “probably damaging” [21], Pmut (available at http://mmb2.pcb.ub.es:8080/PMut/) provides a binary prediction of “neutral” or “pathologic” [22], Sort Intolerant from Tolerant (SIFT) (available at http://sift.jcvi.org) predict whether a change is “tolerated” or “damaging” [23] and MutationTaster2 software (available at http://www.mutationtaster.org/) characterize an amino acid substitution as “polymorphism” or “disease causing” [24] computer algorithms were used to predict whether missense variants were pathological. A brief explanation for these software programs is provided in Pousada et al [8]. The mutations were classified as pathogenic if the score were equal or greater than two.

NNSplice (available at http://fruitfly.org:9005/seq_tools/splice.html), NetGene2 (available at http://www.cbs.dtu.dk/services/NetGene2/), Splice View (available at http://zeus2.itb.cnr.it/~webgene/wwwspliceview_ex.html) and HSF Human (available at http://www.umd.be/HSF/) were used to predict whether changes could affect, create or eliminate donor/acceptor splice sites [8]. The mutations were classified as pathogenic if the score were equal or greater than two.

Minigene constructions and expression

For the c.572G>A (p.G191D) change, we amplified the exon and 200 bp of intronic junctions from the control DNA with High Fidelity Phusion polymerase (Finnzymes, Espoo, Finland) to obtain the wild-type (WT). The amplification conditions were as follows: 98 °C for 30 s, 35 cycles of 98 °C for 10 s, 60 °C for 30 s, 72 °C for 30 s and, finally, 72 °C for 7 min. The amplified fragments were digested and cloned into the XhoI/NheI restriction sites (Nzytech, Lisbon, Portugal) using T4 DNA ligase (New England Biolabs, Ipswich, Massachusetts, USA) in the Exon Trapping Expression Vector p.SPL3 (Invitrogen, San Diego, California, USA). The c.572G>A (p.G191D) construct was generated by site-directed mutagenesis. The primers used for mutagenesis were designed using QuikChange Primer Design (Agilent Technologies, Santa Clara, California, USA). The forward and reverse primers were 5'-gccaggacatggaccgcacgctcga-3' and 5'-tcgagcgtgcggtccatgtcctggc-3', respectively. All constructs were confirmed by direct sequencing.

COS-7 cells (from kidney of Cercopithecus aethiops) were transfected in duplicated by the minigene constructs. Lipofectamine 2000 reagents (Invitrogen, San Diego, CA, USA) were used according to the manufacturer’s instructions. RNA extraction was performed using the Nucleic Acid and Protein Purification kit (NucleoSpin RNA II, Macherey-Nagel, Düren, Germany). RNA was subjected to reverse transcription using the GeneAmp Gold RNA PCR Core Kit (Applied Biosystems, Carlsbad, California). cDNA was amplified and PCR products were sequenced in both senses.

Statistical analysis

A non-parametric test (U Mann-Whitney) was used for comparisons between patients and controls and this study is exploratory. The Chi-square test was used to compare genotypes with clinical and hemodynamic variables and variables were categorized according to the best cut off point by ROC curve. Analyses were supported by the statistical package SPSS v19 for Microsoft and we considered differences statistically significant at values <0.05. Values were expressed as the mean ± SD (standard deviation).

Results

Description of the cohort

This cohort has been described previously by our group [8, 25] and included 57 unrelated PAH patients (29 idiopathic, 19 associated with connective tissue disease, four related to HIV and five porto-pulmonary) (Fig. 1). Samples from PAH patients who agreed to participate in the study were collected between 2008 and 2014. At the time of diagnosis, eight patients were functional class (FC) I, 20 patients were FC II, 25 patients were FC III and four were FC IV. The clinical features of the patients are presented in Table 2.

Fig. 1
figure 1

Nature of the patient cohort. This figure describes the patients involved in this analysis separated by PAH type, the proportion of mutation carriers in the study, the female to male proportion and the mean age at diagnosis. PAH: Pulmonary Arterial Hypertension; IPAH: Idiopathic Pulmonary Arterial Hypertension; Associated Pulmonary Arterial Hypertension; CTD: connective tissue disease; HIV: Human Immunodeficiency virus; P-P: Porto-pulmonary hypertension

Table 2 Clinical features and hemodynamic parameters of patients

In the present study, 55 controls from the general population without a familial history of PAH were included to determine the frequency of the mutations detected in the ENG gene. Samples were kindly provided by the Complexo Hospitalario Universitario de Vigo (Vigo, Spain).

Mutational study of the ENG gene

We found 15 variants of the ENG gene in 44 out of 57 patients. We detected eight different variations first described here and seven changes that have been described elsewhere. The vast majority of these changes were detected in amplicon 7 and 11 (Fig. 2), but we detected the exons 6 and 12 as hotspots for pathogenic mutations. The novel variations did not appear in 55 analyzed controls (110 chromosomes). After an exhaustive in silico analysis, we could identify 8 variations as pathogenic mutations.

Fig. 2
figure 2

Mutational frequency of each of the exons of the ENG gene. The pink color indicates the number of different mutations found in each exon, and the purple color indicates the total of mutations found in each exon for the ENG gene

Missense variations were analyzed with different software programs (PolyPhen, Pmut, Sift and Mutation Taster) to predict their pathogenicity and the impact on the disease. We classified the mutation as potentially pathogenic if two or more programs classified it as pathogenic (Table 3). These analyses classified five missense mutations as pathogenic mutations; however, c.572G>A (p.G191D) has been classified as polymorphism by other studies [2628]. Figure 3 presents the amino acid conservation involved in these missense changes. We observed that the wild-type residues in the p.(S432C) and p.(R554C) mutations are not perfectly conserved between Homo sapiens (human) and ten other species, but are conserved amongst some of the species analyzed.

Table 3 Missense changes found in the coding region of the ENG gene and their classification according to computer algorithms (PolyPhen-2, Pmut, SIFT and MutationTaster2)
Fig. 3
figure 3

Representative sequence electropherograms for the missense variations for the ENG gene in PAH patients with their orthologs. 1: Homo sapiens (sp|P17813#1); 2: Homo sapiens mutated (sp|P17813#1); 3: Mus musculus (sp|Q63961#1); 4: Rattus norvegicus (sp|Q6Q3E8#1); 5: Macaca mulatta (sp|F7BB68#1); 6: Sus scrofa (sp|P37176#1); 7: Oryctolagus cuniculus (sp|G1SSF2#1); 8: Canis familiaris (sp|F1P847#1); 9: Bos taurus (sp|Q1RMV1#1); 10: Equus caballus (sp|F6 W046#1); 11: Loxodonta africana (sp|G3SR82#1); 12: Ailuropoda melanoleuca (sp|G1 M9D6#1)

For the six intronic changes detected, only a duplication (c.991 + 21_991 + 26dupCCTCCC) had been described previously as a polymorphism. This duplication was detected in 35 % of patients included in this study but also in 8 % of controls.

We used other algorithms (NNSplice, NetGene2, Splice View and HSF Human) to predict whether these missense, synonymous and intronic changes could affect donor/acceptor splice sites. We classified the mutation as potentially pathogenic if two or more programs classified it as pathogenic (Table 4).

Table 4 Results from the four different bioinformatic programs used to predict the effect of missense, synonymous and intronic changes on the splicing process in the ENG gene (NNSplice, NetGene2, Splice View and HSF Human)

These pathogenic mutations were detected in 16 patients, four mutations were missense (except c.572G > A (p.G191D), as has been classified as polymorphism by other authors), one synonymous and three were located in the intronic region. Of these patients, seven were classified as IPAH and in nine as APAH.

Study of the c.572G>A (p.G191D) change

This change c.572G>A (p.G191D) was found in ten patients included in this study and was more frequent in IPAH than in patients with APAH. This change was not detected in 110 control alleles (p = 0.001). In patients, the G allele frequency was 0.909 (90 %). Allele A was not detected in controls. This change was not in Hardy-Weinberg Equilibrium (H-WE) in patients (p = 0.617), in contrast to the controls (p < 0.001). BLAST software indicated that the G amino acid (glycine) is an evolutionarily conserved residue (Fig. 4). We checked for alterations in the splicing process using hybrid minigenes for this gene in comparison to the wild type sequence. The mutant construct did not generate a new transcript (Fig. 5). All experiments were performed in duplicate.

Fig. 4
figure 4

a Representative sequence electropherograms for the c.572G > A (p.(D191G)) mutation for the ENG gene in PAH patients. b Different orthologs for this mutation. c Mutational frequency for this pathogenic mutation in IPAH and APAH patients. 1: Homo sapiens (sp|P17813#1); 2: Homo sapiens mutated (sp|P17813#1); 3: Mus musculus (sp|Q63961#1); 4: Rattus norvegicus (sp|Q6Q3E8#1); 5: Macaca mulatta (sp|F7BB68#1); 6: Sus scrofa (sp|P37176#1); 7: Oryctolagus cuniculus (sp|G1SSF2#1); 8: Canis familiaris (sp|F1P847#1); 9: Bos taurus (sp|Q1RMV1#1); 10: Equus caballus (sp|F6W046#1); 11: Loxodonta africana (sp|G3SR82#1); 12: Ailuropoda melanoleuca (sp|G1M9D6#1)

Fig. 5
figure 5

In vitro splicing assay for the c.572G > A (p.G191D) change identified in the ENG gene. a Electropherogram of the transcript obtained indicates the molecular characterization of the effect of the studied variant. b Graphical representation of the effect of p.G191D change in mRNA processing. c Electrophoresis of wild-type and mutant construction. SDS and SA2: pSPL3 vector exons, where the inserts to study are cloned

Association with clinical features and hemodynamic parameters

None of the clinical features or hemodynamic parameters exhibited statistically significant differences, except for age at diagnosis (p = 0.040) and the 6-min walking test (p = 0.040). Patients with pathogenic mutations in ENG gene exhibited disease symptoms 8 years earlier and were diagnosed earlier than patients with a negative mutational screening for ENG, BMPR2, ACVRL1 (Activin A Receptor Type II-Like 1) and KCNA5 (Potassium voltage-gated channel, shaker-related subfamily, member 5) genes (Table 2). However, five patients with ENG pathogenic mutations were also carriers for another mutation in the BMPR2 or ACVRL1 genes (Fig. 6) as described by Pousada et al [8]. When removing these patients for statistical analysis, only age at diagnosis was significantly different (mean 9 years early, p = 0.040).

Fig. 6
figure 6

Mutational analysis of patients with multiple pathogenic mutations in analyzed genes

The c.572G>A (p.G191D) change was associated with an early age at diagnosis (mean 10 years earlier, p = 0.035) and lower CI (p = 0.049). Finally, this change was more prevalent in IPAH patients (p = 0.040). Other clinical and hemodynamic parameters exhibited no statistically significant differences. These results should be analyzed carefully because all carriers for c.572G>A (p.G191D) variation but one, were also carriers for mutations in others genes (BMPR2, ACVRL1 and KCNA5).

Discussion

Mutations in the ENG gene have been described in up to 88 % of HHT patients, including some with PAH associated with HHT [29, 30]. In this study we have identified a higher number of pathogenic mutations in comparison with the results showed by other analysis [4, 7, 17, 3133]. All research conducted in ENG gene have been performed in IPAH or HPAH patients, but the study by Pfarr et al [7] described a small number of pathogenic mutations in patients with congenital heart disease associated to PAH. In 29 children with IPAH or HPAH and 11 with APAH due to congenital heart disease without any symptoms or familial history of HHT, Pfarr et al [7] found 2 patients (5 %) carriers of mutations in the ENG gene. However, in our cohort we included patients with IPAH and associated with other pathologies. This is the first mutational analysis of the ENG gene in PAH patients associated to connective tissue disease, human immunodeficiency virus and porto-pulmonary hypertension.

We identified ENG mutations in 16 subjects, a significantly higher percentage. We detected 5 ENG mutations with potential pathogenicity not yet described and three described sequence variants. Furthermore, with the in silico analysis we were able to classify synonymous mutations and mutations located in intronic junctions as pathogenic mutations. However, other studies only focused on the analysis of missense and nonsense mutations [7, 32]. Perhaps this fact can significantly increase the percentage of pathogenic mutations in our patients. For these analyses we used eight bioinformatic softwares that analyzed the pathogenicity of the mutations. We considered PolyPhen, Pmut, Sift and Mutation Taster softwares that analyze the amino acid conservation, the protein function or the protein structure [2124]. However, these softwares show some differences in the criteria used to establish the pathogenicity character of the variation. Some of them included more information as the description of the variants when is possible, the implication in the splicing process or the presence of enhancer sequences. Besides, we used four softwares that analyze the implication of the splicing changes in the mRNA processing. In silico analysis is not totally reliable, and for this reason we believe that analyze this variants with several softwares is necessary to give us a greater approach to catalogue a variant as polymorphism or as pathogenic mutation. Functional studies would be necessary to confirm the pathogenicity of these variants. Aparisi et al [34] described that after exhaustive in silico analysis with splicing softwares, only a few mutations classified as pathogenic resulted really pathogenic in the functional splicing analysis performed. In our study, we have to take into account the fact that none of these variations classified as pathogenic have been found in healthy controls and the c.1633G>A (p.(G545S)) mutation was been classified as pathogenic by other research group [7].

We detected two hot spots for exons 6 and 12 in the ENG gene. These exons are located in the extracellular region (Zona pellucida-like domain) [7], a very important area rich in glycosylation sites and cysteine residues [15, 20]. This region has a characteristic pattern of preserved residues [15, 20, 35]. Furthermore, Ali et al have reported that missense mutations in this region for ENG gene led to a decrease or disappearance of cell surface expression of the protein [36]. Likewise, many missense mutations located in an orphan domain, situated in a Zona pellucida-like domain, resulted in protein misfolding, altering the subcellular localization [35]. It is likely that the mutated protein was retained by the endoplasmic reticulum (ER) quality control machinery [26, 36]. As a result, the protein becomes trapped in the rough ER and is subjected to ER associated protein decay [26]. Thus, disruption of the downstream signaling of the TGF-β pathway might be caused by mutations affecting both the TGF-β/ALK1 and TGF-β/ALK5 balance and the endothelial-cell growth potential [3739]. The number and class of molecules involved in this pathway, which differ among cells, underlie the great complexity and versatility of TGF-β signaling [31]. Moreover, in vitro studies on pulmonary artery smooth muscle cells from IPAH patients have indicated growth abnormalities [40].

Missense changes found in these patients are located in aminoacidc residues highly conserved except p.(S432C) and p.(R554C). These variations could be explained as polymorphic change with evolutionary effects. Serine is a non-essential polar amino acid that is neutrally charged, and arginine is non-polar, essential and neutrally charged. However, cysteine is non-essential and negatively charged. The change in charge could be compensated with another mutation in a region in close three-dimensional proximity. Gallione et al [20] reported that cysteine amino acids are involved in disulfide bridging. These mutations can produce alterations in the protein’s structure that affect its functionality; the mutant allele could have a dominant negative effect over the wild type allele, causing serious consequences for carrier patients as have been described by John et al [41] in the BMPR2 gene in patients with PAH.

The c.572G>A (p.(G191D) change has been previously described as a polymorphism or rare variant [2628] despite being classified as pathogenic with four of the computer programs used. For this reason and because it is found at a very low frequency in the Spanish and European control population, we performed functional studies for this mutation to verify in vitro its pathogenicity. The analysis with the minigenes assay did not detect any change in the splicing process. Förg T et al [26] performed several colocalization experiments with fluorescence microscopy, and the authors also classified it as a polymorphism. Nonetheless, it is possible that this change may act through other mechanisms, as the complete role of ENG is still unknown and requires further functional studies.

Furthermore, we found a pathogenic synonymous change. Synonymous changes could interfere with the splicing accuracy, translation fidelity, mRNA structure and protein folding. Furthermore, these mutations may decrease the half-life of mRNA, leading to downregulation of the protein expression [8, 33]. Synonymous codons are translated at lower levels than standard codons, since specific tRNA levels are decreased [42]. Functional studies for synonymous mutations, intronic changes and intronic duplication would be very interesting, as the role of these changes is unknown, and a functional approach could help us to improve our knowledge of the disease.

In addition, we found that carriers of pathogenic mutations were younger at diagnosis. This fact, together with previous studies, indicates significant heterogeneity in the genetic background of PAH. Mutations in the BMPR2 gene are most common in PAH patients, but other genes may be related, including ACVRL1 or KCNA5 [8]. All patients in this study were analyzed for mutations in these three genes (BMPR2, ACVRL1 and KCNA5) [8]. For the 57 patients analyzed for ENG gene, 11 out of 16 patients exhibited only a mutation in the ENG gene. Mutations in the ENG gene are quite prevalent in our cohort of PAH patients, can influence the development of the pathology and did not appear in 55 control samples.

The ability of ENG to collaborate in the pathogenesis is highly variable, as described by Mallet et al [43]. The mutant protein could act in a haploinsufficient manner, interacting with the wild type protein and interfering in the normal endoglin function; alternatively, reduction or loss of the cell surface expression of the mutant protein has been described. As noted by John et al [41], we cannot exclude other mechanisms, including the ability to interact with other partners or to activate other signaling pathways.

When we compared the hemodynamic and clinical parameters between patients with and without pathogenic mutations, patients with mutations exhibited a significantly earlier age at diagnosis (8 years compared with patients without mutations) and a lower 6MWT. Therefore, we cannot exclude the possibility that these differences may be due to the small number of patients in our series. PAH exhibits highly variable clinical parameters, and clinical diagnosis is complicated by the heterogeneous outcome of disease manifestation; hence, additional diagnostic tools are required to perform early diagnosis in affected individuals.

Considering the patients with mutations only in the ENG gene, we did not find significant differences in clinical or hemodynamic parameters, but patients were diagnosed at an earlier age compare with patients without mutations. Endoglin exhibits two different splice isoforms, short (S) and long (L). Although the most common isoform of endoglin in endothelial cells is L-endoglin, Blanco et al [44] reported that short S-endoglin expression contributes to the cardiovascular pathology associated with age in vivo and in vitro. These results suggest that S-endoglin expression affects the senescent program of endothelial cells when S-endoglin is upregulated instead of being solely responsible for senescence. Furthermore, Liu et al [45] reported that endoglin is also related to critical function in the development of the vascular system in mouse embryonic stem cells, this could explain that patients with pathogenic mutations have an early presentation of the disease.

Previous studies in the BMPR2 gene indicate that PAH patients carrying a mutation have an onset of disease approximately 10 years earlier than non-carriers [4] and Liu et al [46] suggest that the phenotype of PAH patients with BMPR2 mutations are influenced by gender. These male patients have a more penetrant phenotype [46]. The former statement of the BMPR2 gene could be extrapolated to the ENG gene, according to our results, but we did not detect gender differences in this study.

As almost all of our patients with the c.572G>A (p.G191D) change exhibited a pathogenic mutation in other genes (BMPR2, ACVRL1 and KCNA5), we investigated whether the presence of this change could modify the phenotype. Pfarr et al [7] found significant differences for a low PVR value when they compared carriers of mutations in the BMPR2, ACVRL1, ENG and SMAD genes with non-pathogenic mutation carriers. Moreover, we found significant differences in the age at diagnosis, CI and PAH types when comparing hemodynamics and clinical parameters between patients with the c.572G>A (p.G191D) change vs patients without pathogenic mutations in none analyzed genes. Patients harboring this mutation exhibited significantly smaller CI values. We found that this mutation was more prevalent in patients with IPAH than in those APAH. Finally, this mutation appears in patients who are diagnosed 10 years earlier than non-carriers. As the specific mechanism for ENG is not yet characterized and its relation with other PAH genes remain unclear, these data should be cautiously interpreted.

Five patients with pathogenic mutations in the ENG gene also exhibited a mutation in another gene. Two of these patients with p.(Q166Q) mutation in the ENG gene [47] are carriers of p.(Q92L) and p.(V341M) BMPR2 gene mutations, classified as pathogenic [8]. Patient 3, with a c.1272+6A > T mutation, was also a carrier of the p.(W298*) mutation in the BMPR2 gene [8]. The last two patients, with p.(G545S) [7] and c.817+17T>A mutations, also harbored the p.(S232T) and p.(T243N) ACVRL1 gene mutations, respectively. Mallet et al [43] described several patients with pathogenic mutations in different genes, including ENG, in HHT patients. The authors proposed that one of the two mutations classified as pathogenic could be a rare variant [43], unlikely to cause PAH. However, as observed in other human pathologies, oligogenic inheritance of PAH with a major causal gene should not be excluded [48]. Rodríguez-Viales et al [49] proposed that additional variations can produce a more severe phenotype and an early disease. The evaluation of the total mutation load could be the way to understand how mutations in different genes could be responsible for the disease [50, 51]. This fact further supports the importance of the functional analysis of these mutants.

The low number of patients included in this study represents a handicap to our ability to draw stronger conclusions, but the comprehensive investigation and complete follow-up of all cases adds value to our data. The high clinical evaluation of these patients and the expertise in the molecular methodology field make us confident about the elevated presence of putative pathogenic mutations in these patients. A functional study should be necessary in order to corroborate the functional impact of these mutations.

Conclusion

To conclude, we described mutations in the ENG gene in IPAH and APAH patients, some of which have not been previously described. The vast majority of the mutations found in this study are private, making difficult to establish a correlation between the phenotype and a particular mutation. For this reason, the genotype-phenotype correlation was performed according to all mutations found in a group of patients. Mutations in the ENG gene may influence the clinical status of the disease. Moreover, genetic analysis in patients with PAH may be of clinical relevance and demonstrates the complexity of the genetic background. A better understanding of the molecular basis will support the future design of individualized treatments according to the genetic background of each patient.

Abbreviations

6MWT, 6 min walking test; ACVRL1, activin A Receptor Type II-Like 1; APAH, associated pulmonary arterial hypertension; BLAST, basic local alignment search tool; BMPR2, bone morphogenetic protein type 2 receptor; CAEI, Comité Autonómico de Ética da Investigación; CI, cardiac index; CT, tomography computerized; ENG, Endoglin; ER, endoplasmic reticulum; ERS-ESC, European Respiratory Society-European Society of Cardiology; FC, functional class; HHT, hereditary hemorrhagic telangiectasia; HIV, human immunodeficiency virus; HPAH, heritable pulmonary arterial hypertension; H-WE, Hardy-Weinberg Equilibrium; IPAH, idiopathic pulmonary arterial hypertension; KCNA5, potassium voltage-gated channel, shaker-related subfamily, member 5; L: long; mPaP, mean pulmonary artery pressure; mRNA, messenger RNA; PAH, pulmonary arterial hypertension; PCR, polymerase chain reaction; PVR, pulmonary vascular resistance; S, short; SD, standard deviation; sPaP, systolic pulmonary artery pressure; TGF-β, transforming growth factor beta.

References

  1. Galié N, Hoeper MM, Humbert M, Torbicki A, Vachiery JL, Barbera JA, et al. Guidelines for the diagnosis and treatment of pulmonary hypertension. Eur Heart J. 2009;30(20):2493–537.

    Article  PubMed  Google Scholar 

  2. Simonneau G, Gatzoulis MA, Adatia I, Celermajer D, Denton C, Ghofrani A, et al. Updated clinical classification of pulmonary hypertension. J Am Coll Cardiol. 2013;62(25 Suppl):D34–41.

    Article  PubMed  Google Scholar 

  3. Peacock AJ, Murphy NF, McMurray JJV, Caballero L, Stewart S. An epidemiological study of pulmonary arterial hypertension. Eur Respir J. 2007;30:104–9.

    Article  CAS  PubMed  Google Scholar 

  4. Machado RD, Eickelberg O, Elliott G, Geraci MW, Hanaoka M, Loyd JE, et al. Genetics and genomics of pulmonary arterial hypertension. J Am Coll Cardiol. 2009;54 Suppl 1:S32–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Sanchez O, Marié E, Lerolle U, Wermert D, Israël-Biel D, Meyer G. Pulmonary arterial hypertension in women. Rev Mal Respir. 2010;27:e79–87.

    Article  CAS  PubMed  Google Scholar 

  6. Yang X, Long L, Reynolds PN, Morrell NW. Expression of mutant BMPR-II in pulmonary endothelial cells promotes apoptosis and a release of factors that stimulate proliferation of pulmonary arterial smooth muscle cells. Pulmonary Circulation. 2010;1(1):103–11.

    Article  CAS  Google Scholar 

  7. Pfarr N, Fischer C, Ehlken N, Becker-Grünig T, López-González V, Gorenflo M, et al. Hemodynamic and genetic analysis in children with idiopathic, heritable, and congenital heart disease associated pulmonary arterial hypertension. Respir Res. 2013;14:3–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Pousada G, Baloira A, Vilariño 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(6), e100261.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Davies RJ, Morrell NW. Molecular mechanisms of Pulmonary Arterial Hypertension: Morphogenetic protein type II receptor. Chest. 2008;134:1271–7.

    Article  CAS  PubMed  Google Scholar 

  10. Girerd B, Montani D, Eyries M, Yaici A, Sztrymf B, Coulet F, et al. Absence of influence of gender and BMPR2 mutation type on clinical phenotypes of pulmonary arterial hypertension. Respir Res. 2010;1:73.

    Article  Google Scholar 

  11. Pffar N, Szamalek-Hoegel J, Fischer C, Hinderhofer K, Nagel C, Eklken N, et al. Hemodynamic and clinical onset in patients with hereditary pulmonary arterial hypertension and BMPR2 mutations. Respir Res. 2011;12:99–109.

    Article  Google Scholar 

  12. Harrison RE, Flanagan JA, Sankelo M, Abdalla SA, Rowell J, Machado RD, et al. Molecular and functional analysis identifies ALK-1 as the predominant cause of pulmonary hypertension related to hereditary haemorrhagic telangiectasia. J Med Genet. 2003;40:865–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. López-Novoa JM, Bernabeu C. The physiological role of endoglin in the cardiovascular system. Am J Physiol Heart Circ Physiol. 2010;299:H959–74.

    Article  PubMed  Google Scholar 

  14. Sadick H, Hage J, Goessler U, Stern-Straeter J, Riedel F, Hoermann K, et al. Mutation analysis of "Endoglin" and "Activin receptor-like kinase" genes in German patients with hereditary hemorrhagic telangiectasia and the value of rapid genotyping using an allele-specific PCR-technique. BMC Med Genet. 2009;10:53.

    PubMed  PubMed Central  Google Scholar 

  15. Abdalla SA, Letarte M. Hereditary haemorrhagic telangiectasia: current views on genetics and mechanisms of disease. J Med Genet. 2006;43(2):97–110.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Gore B, Izikki M, Mercier O, Dewachter L, Fadel E, Humbert M, et al. Key role of the endothelial TGF-β/ALK1/endoglin signaling pathway in humans and rodents pulmonary hypertension. PLoS One. 2014;9(6), e100310.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Montani D, Günther S, Dorfmüller P, Perros F, Girerd B, García G, et al. Pulmonary arterial hypertension. Orphanet J Rare Dis. 2013;8:97.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Simonneau G, Gatzoulis MA, Adatia I, Celermajer D, Denton C, Ghofrani A, et al. Updated clinical classification of pulmonary hypertension. J Am Coll Cardiol. 2013;62:25s.

    Article  Google Scholar 

  19. Galiè N, Hoeper MM, Humbert M, Torbicki A, Vachiery JL, Barbera JA, et al. Guidelines for the diagnosis and treatment of pulmonary hypertension: the Task Force for the Diagnosis and Treatment of Pulmonary Hypertension of the European Society of Cardiology (ESC) and the European Respiratory Society (ERS), endorsed by the International Society of Heart and Lung Transplantation (ISHLT). Eur Heart J. 2009;30(20):2493–537.

    Article  PubMed  Google Scholar 

  20. Gallione CJ, Klaus DJ, Yeh EY, Stenzel TT, Xue Y, Anthony KB, et al. Mutation and expression analysis of the Endoglin Gene in Hereditary Hemorrhagic Telangiectasia reveals null alleles. Hum Mutat. 1998;11:286–94.

    Article  CAS  PubMed  Google Scholar 

  21. Adzhubei I, Jordan DM, Sunyaev SR. Predicting functional effect of human missense mutations using PolyPhen-2. Curr Protoc Hum Genet. 2013;7(7):20.

    PubMed  Google Scholar 

  22. Ferrer-Costa C, Orozco M, de la Cruz X. Sequence-based prediction of pathological mutations. Proteins. 2004;57(4):811–9.

    Article  CAS  PubMed  Google Scholar 

  23. Kumar P, Henikoff S, Ng PC. Predicting the effects of coding nonsynonymous variants on protein function using the SIFT algorithm. Nat Protoc. 2009;4(7):1073–81.

    Article  CAS  PubMed  Google Scholar 

  24. Schwarz JM, Cooper DN, Schuelke M, Seelow D. MutationTaster2: mutation prediction for the deep-sequencing age. Nat Methods. 2014;11(4):361–2.

    Article  CAS  PubMed  Google Scholar 

  25. Pousada G, Baloira A, Valverde D. Molecular and clinical analysis of TRPC6 and AGTR1 genes in patients with pulmonary arterial hypertension. Orphanet J Rare Dis. 2015;10:1.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Förg T, Hafner M, Lux A. Investigation of endoglin wild-type and missense mutant protein heterodimerisation using fluorescence microscopy based IF, BiFC and FRET analyses. PLoS One. 2014;9(7), e102998.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Lesca G, Plauchu H, Coulet F, Lefebvre S, Plessis G, Odent S, et al. Molecular screening of ALK1/ACVRL1 and ENG genes in hereditary hemorrhagic telangiectasia in France. Hum Mutat. 2004;23(4):289–99.

    Article  CAS  PubMed  Google Scholar 

  28. Abdalla SA, Cymerman U, Rushlow D, Chen N, Stoeber GP, Lemire EG, et al. Novel mutations and polymorphisms in genes causing hereditary hemorrhagic telangiectasia. Hum Mutat. 2005;25(3):320–1.

    Article  PubMed  Google Scholar 

  29. Harrison RE, Berger R, Haworth SG, Tulloh R, Mache CJ, Morrell NW, et al. Transforming growth factor-receptor mutations and pulmonary arterial hypertension in childhood. Circulation. 2005;111(4):435–41.

    Article  CAS  PubMed  Google Scholar 

  30. McDonald J, Wooderchak-Donahue W, VanSant-Webb C, Whitehead K, Stevenson DA, Bayrak-Toydemir P. Hereditary hemorrhagic telangiectasia: genetics and molecular diagnostics in a new era. Front Genet. 2015;6:1.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Wipff J, Kahan A, Hachulla E, Sibilia J, Cabane J, Meyer O, et al. Association between an endoglin gene polymorphism and systemic sclerosis-related pulmonary arterial hypertension. Rheumatology. 2007;46(4):622–5.

    Article  CAS  PubMed  Google Scholar 

  32. Machado RD, Southgate L, Eichstaedt CA, Aldred MA, Austin ED, Best DH, et al. Pulmonary arterial hypertension: a current perspective on established and emerging molecular genetic defects. Hum Mutat. 2015;36(12):1113–27.

    Article  CAS  PubMed  Google Scholar 

  33. Best DH, Austin ED, Chung WK, Elliott CG. Genetics of pulmonary hypertension. Curr Opin Cardiol. 2014;29(6):520–7.

    Article  PubMed  Google Scholar 

  34. Aparisi MJ, García-García G, Aller E, Sequedo MD, Martínez-Fernández de la Cámara C, Rodrigo R, et al. Study of USH1 splicing variants through minigenes and transcript analysis from nasal epithelial cells. PLoS One. 2013;8(2), e57506.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Schulte C, Geisthoff U, Lux A, Kupka S, Zenner HP, Blin N, et al. High frequency of ENG and ALK1/ACVRL1 mutations in German HHT patients. Hum Mutat. 2005;25(6):595.

    Article  PubMed  Google Scholar 

  36. Ali BR, Ben-Rebeh I, John A, Akawi NA, Milhem RM, Al-Shehhi NA, et al. Endoplasmic reticulum quality control is involved in the mechanism of endoglin-mediated hereditary Haemorrhagic Telangiectasia. PLoS One. 2011;6(10), e26206.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Czech A, Fedyunin I, Zhang G, Ignatova Z. Silent mutations in sight: co-variations in tRNA abundance as a key to unravel consequences of silent mutations. Mol BioSyst. 2010;6:1767–72.

    Article  CAS  PubMed  Google Scholar 

  38. Fernández LA, Sanz-Rodriguez F, Blanco FJ, Bernabéu C, Botella LM. Hereditary hemorrhagic telangiectasia, a vascular dysplasia affecting the TGF-beta signaling pathway. Clin Med Res. 2006;4(1):66–78.

    Article  Google Scholar 

  39. Upton PD, Morrell NW. The transforming growth factor-β-bone morphogenetic protein type signalling pathway in pulmonary vascular homeostasis and disease. Exp Physiol. 2013;98(8):1262–6.

    Article  CAS  PubMed  Google Scholar 

  40. Chaouat A, Coulet F, Favre C, Simonneau G, Weitzenblum E, Soubrier F, et al. Endoglin germline mutation in a patient with hereditary haemorrhagic telangiectasia and dexfenfluramine associated pulmonary arterial hypertension. Thorax. 2004;59(5):446–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. John A, Kizhakkedath P, Al-Gazali L, Ali BR. Defective celular trafficking of the bone morphogenetic protein receptor type II by mutations underlying familial pulmonary arterial hypertension. Gene. 2015;561:148–56.

    Article  CAS  PubMed  Google Scholar 

  42. Chaney JL, Clark PL. Roles for synonymous codon usage in protein biogenesis. Annu Rev Biophys. 2015;44:13–48.

    Article  Google Scholar 

  43. Mallet C, Lamribet K, Giraud S, Dupuis-Girod S, Feige JJ, Bailly S, et al. Functional analysis of endoglin mutations from hereditary hemorrhagic telangiectasia type 1 patients reveals different mechanisms for endoglin loss of function. Hum Mol Genet. 2015;144(6):261–4.

    Google Scholar 

  44. Blanco FJ, Grande MT, Langa C, Oujo B, Velasco S, Rodríguez-Barbero A, et al. S-endoglin expression is induced in senescent endothelial cells and contributes to vascular pathology. Circ Res. 2008;103(12):1383–92.

    Article  CAS  PubMed  Google Scholar 

  45. Liu Z, Lebrin F, Maring JA, Van Den Driesche S, van der Brink S, van Dinther M, et al. ENDOGLIN is dispensable for vasculogenesis, but required for vascular endothelial growth factor-induced angiogenesis. PLoS One. 2014;9(1), e86273.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Liu D, Wu WH, Mao YM, Yuan P, Zhang R, Ju FL, et al. BMPR2 mutations influence phenotype more obviously in male patients with pulmonary arterial hypertension. Circ Cardiovasc Genet. 2012;5(5):511–8.

    Article  CAS  PubMed  Google Scholar 

  47. Pousada G, Baloira A, Valverde D. Pulmonary arterial hypertension and portal hypertension in a patient with hereditary hemorrhagic telangiectasia. Med Clin. 2014;144(6):261–4.

    Article  Google Scholar 

  48. Daniels AB, Sandberg MA, Chen J, Weigel-DiFranco C, Fielding Hejtmancic J, Berson EL. Genotype-phenotype correlations in Bardet-Biedl syndrome. Arch Ophtthalmol. 2012;130(7):901–7.

    Google Scholar 

  49. Rodríguez-Viales R, Eichstaedt CA, Ehlken N, Fischer C, Lichtblau M, Grünig E, et al. Mutation in BMPR2 promoter: a ‘second hit’ for manifestation of pulmonary arterial hypertension? PLoS One. 2015;10(7), e0133042.

    Article  Google Scholar 

  50. Forsythe E, Beales PL. Bardet-Biedl syndrome. Eur J Hum Genet. 2013;21:8–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Castro-Sánchez S, Álvarez-Satta M, Cortón M, Guillén E, Ayuso C, Valverde D. Exploring genotype-phenotype relationships in Bardet-Biedl syndrome families. J Med Genet. 2015;52(8):503–13.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We are grateful to the patients who participated in our research. We thank the physicians who participated in the collection of patients and data (Carlos Vilariño Pombo, José Manuel Cifrián, Olalla Cástro-Añón), and we acknowledge the Asociación Española de Hipertensión Pulmonar, who provided DNA samples and their consent for inclusion in this investigation.

Funding

This study was supported by the grants IN-202-05 from SOGAPAR, CO-0085-10 from Actelion Pharmaceuticals and INBIOMED 2009-063 Xunta de Galicia. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Availability of data and materials

Additional data and information are available upon request.

Authors' contributions

GP conceived of the study and its design, performed genetic research, data statistical analysis and interpretation, correlation genotype/phenotype and draft the manuscript. AB conceived the study, realized the collection of patients and data, statistical analysis and drafts the manuscript. DF performed genetic research, data statistical analysis and draft the manuscript. MN realized the collection of patients and data and statistical analysis. DV conceived the study and participated in its design and coordination, helped to experiments and draft the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

All patients are included in the CHUVI DNA Biobank (Biobanco del Complejo Hospitalario Universitario de Vigo). Patients signed an informed consent and the Autonomic Ethics Committee approved the study (Comité Autonómico de Ética da Investigación de Galicia - CAEI de Galicia), following the clinical-ethical guidelines of the Spanish Government and the Helsinki Declaration.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Diana Valverde.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Pousada, G., Baloira, A., Fontán, D. et al. Mutational and clinical analysis of the ENG gene in patients with pulmonary arterial hypertension. BMC Genet 17, 72 (2016). https://doi.org/10.1186/s12863-016-0384-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12863-016-0384-3

Keywords