PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of bmcpmBioMed Centralsearchsubmit a manuscriptregisterthis articleBMC Pulmonary Medicine
 
BMC Pulm Med. 2017; 17: 57.
Published online 2017 April 7. doi:  10.1186/s12890-017-0400-z
PMCID: PMC5383973

A burden of rare variants in BMPR2 and KCNK3 contributes to a risk of familial pulmonary arterial hypertension

Abstract

Background

Pulmonary arterial hypertension (PAH) is a severe lung disease with only few effective treatments available. Familial cases of PAH are usually recognized as an autosomal dominant disease, but incomplete penetrance of the disease makes it difficult to identify pathogenic variants in accordance with a Mendelian pattern of inheritance.

Methods

To elucidate the complex genetic basis of PAH, we obtained whole exome- or genome-sequencing data of 17 subjects from 9 families with heritable PAH and applied gene-based association analysis with 9 index patients and 300 PAH-free controls.

Results

A burden of rare variants in BMPR2 significantly contributed to the risk of the disease (p = 6.0 × 10−8). Eight of nine families carried four previously reported single nucleotide variants and four novel insertion/deletion variants in the gene. One of the novel variants was a large 6.5 kilobase-deletion. In the remaining one family, the patient carried a pathogenic variant in a member of potassium channels, KCNK3, which was the first replicative finding of channelopathy in an Asian population.

Conclusions

The variety of rare pathogenic variants suggests that gene-based association analysis using genome-wide sequencing data from increased number of samples is essential to tracing the genetic heterogeneity and developing an appropriate panel for genetic testing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12890-017-0400-z) contains supplementary material, which is available to authorized users.

Keywords: Clinical genetic testing, Gene-based association study, Next generation sequencing, Pulmonary arterial hypertension

Background

Pulmonary arterial hypertension (PAH, MIM #178600) is a rare vascular lung disease presenting increased pulmonary vascular resistance and elevation of mean pulmonary arterial pressure, leading to a grave prognosis of right heart failure without treatment. While survival rates are increasing with a number of recently developed treatments such as epoprostenol, the optimal care of patients with these therapies is unclear due to phenotypic variations and genetic backgrounds.

Despite the complex etiology of PAH due to incomplete penetrance and genetic heterogeneity, multiple genes that cause PAH have been discovered during the last few decades. BMPR2 variants have been identified in not only >70% of heritable PAH (HPAH), but also 10–40% of idiopathic PAH (IPAH) [16]. Variants of other genes including KCNK3, ACVRL1, ENG, CAV1, and the SMAD family are also rare causes of PAH [711]. Furthermore, common variants in CBLN2 and KCNA5 were also reported to be associated with the risk of PAH in European Caucasians [12, 13].

Although identification of individuals who carry genetic variants that increase the risk of developing PAH offers an opportunity for earlier diagnosis and finding a therapeutic strategy, the majority of previous studies was only focused on the protein coding regions of the most frequently mutated gene BMPR2 with the use of conventional methods such as Sanger sequencing [6, 14]. Thus, for the patients who have no mutation in BMPR2, i.e., around 30% of HPAH and 60–90% of IPAH, another approach, which is practical for multiple genes, is necessary. Considering that current state-of-the-art sequencing technologies allow us to access exome- and genome-wide variants with reasonable costs, unbiased screening of pathogenic variants is beneficial to continue expanding the genetic diagnosis catalogue. However, conventional segregation-based approaches, e.g., linkage analysis, do not have enough power to pinpoint the true pathogenic variants from these genome-wide candidate variants without functional evaluations, especially for diseases with phenotypic and genetic heterogeneity. In this study, applying a gene-based association test, we statistically evaluate the significance of rare variant enrichment in genes responsible for PAH. The strategy we employed here would be useful to unbiasedly elucidate the pathogenicity of multiple rare variants arising from independent founder events in Mendelian diseases as well as common diseases.

Methods

Subjects

We consecutively enrolled five families with PAH and four individual cases with a family history of PAH to this study (Fig. 1). The patients have been diagnosed in the National Hospital Organization Okayama Medical Center between 1996 and 2014. All subjects who participated in our study were approved by the Institutional Review Board of our institutes in which donors gave written informed consent in accordance with institutional and national guidelines.

Fig. 1
Segregation of the pathogenic variants identified in 9 familial PAH. a Eight families carried BMPR2 variants and one family carried a KCNK3 variant. Nucleotide and amino acid changes for BMPR2 and KCNK3 are described on NM_001204.6 and NM_002246.2, respectively. ...

Next-generation sequencing and data analysis

To understand comprehensive genetic background of these 9 PAH families, we applied whole exome- or genome-sequencing to 12 patients and 5 healthy family members whose DNA were available (Fig. 1 and Additional file 1: Table S1). For the exome sequencing, DNA fragments were enriched by SureSelectXT Human All Exon v4 + UTR (Agilent Technologies, Santa Clara, CA, USA) and then applied to SOLiD™ 5500XL sequencer (Thermo Fisher Scientific inc., Waltham, MA, USA). The whole genome sequencing was conducted with the Illumina HiSeq X sequencer (Illumina Inc., San Diego, CA, USA). After aligning the sequence reads onto the reference genome (NCBI Build 37) using the Burrows-Wheeler Aligner [15], downstream processes including the duplication removal, the recalibration of base quality values, the local realignment, the variant call, and the variant quality score recalibration were analyzed using GATK [16]. The variants were called with an exome sequencing data set of 300 control samples obtained from the Human Genome Variation Database (accession ID: HGV0000004) [17]. The resulting VCF file has been deposited on the same database under accession HGV0000005.

Quality control for association study

After removing variants that were assigned as low quality by the GATK VariantRecalibrator, additional filters were applied to extract high quality variants such as low call rate (<0.9), excessive strand bias (FS > 50), haplotype score (> 5), deviation from Hardy-Weinberg equilibrium (InbreedingCoeff > 0.3), mapping quality of the reads (MQ < 35), excess of zero mapping quality (MQ0 > 100), bias of mapping quality between reference and alternative alleles (MQRankSum < 13), coverage over sample (DP/sample < 10), positional bias of the reads (ReadPosRankSum > 5), quality over depth (QD < 8), and low LOD score (VQSLOD < 0). For the gene-based association analysis, we selected likely protein damaging variants (premature termination, splice site, missense, and indels on exons) to perform the Variable Threshold (VT) test [18] implemented in Variant Association Tools [19].

Annotation and screening of pathogenic variants

All identified variants were annotated using ANNOVAR [20]. Candidate pathogenic variants were screened according to the registrations and frequencies of the variants in the public databases: dbSNP (Build 147) [21], The 1000 Genomes (November 2010 data release) [22], The 10Gen Data Set (version 1.04) [23], NHLBI GO Exome Sequencing Project (ESP6500SI) [24], the Human Genetic Variation Database [17] or ClinVar [25]. For missense variants, PolyPhen-2 [26] and Mutation Taster [27], LRT [28] and PhyloP [29] score were obtained from the dbNSFP database [30]. Damaging effects of splice site variants were evaluated with MaxEntScan [31] and Human Splicing Finder [32].

Results

Gene-based association study

To identify genes responsible for the pathogenesis of PAH, we applied a gene-based association test (Variable Threshold test [18]) to the 60,367 damaging variants extracted from the nine PAH patients and the 300 control samples. These variants were located within 10,744 gene regions. Despite the small sample size, the burden of association between BMPR2 and PAH was highly significant (p = 6.0 × 10−8) compared to the genome-wide significance threshold (p < 2.4 × 10−6) after Bonferroni correction for approximately 21,000 genes (Fig. 2). No other gene was found beyond the threshold.

Fig. 2
Gene-based association analysis. Quantile-quantile plot for VT test of the nine cases and 300 controls are shown. BMPR2 was the only gene that surpassed the genome-wide significance threshold for 21,000 genes

Identifying the pathogenic variants

The spectrum of rare variants found in BMPR2 is summarized in Table 1. Of the nine families, four carried previously reported single nucleotide pathogenic variants (2 missenses, 1 nonsense, and 1 splice site) [1, 4, 6] and four carried novel insertions/deletions (indels) in this gene (88.9%). One of the novel indels was a large deletion of 6.5 kilobases in length by which one allele lacks the entire region of exon 3 (Additional file 2: Figure S1). Two variants are suspected to be pathogenic although showing incomplete penetrance, since clinically unaffected subjects in the families harbored the same variants found in the patients (Table 1). There were no BMPR2 variants observed in the remaining one of nine families, but we identified one heterozygous missense variant (p.Gly203Asp) in KCNK3 by screening the previously reported pathogenic variants [25]. This variant was shown to disrupt the ion-channel function by patch-clump analysis [11]. None of the pathogenic variants we identified was observed in the 300 control samples or in the public database for the Japanese population [17]. Of these, all three missense variants were occurred at highly conserved nucleotides among vertebrates and were assumed to be damaging to the protein function by at least three in silico prediction programs [2629] (Table 1).

Table 1
Summary of pathogenic variants and clinical features

Discussion

To our knowledge, this is the first report of gene-based genome-wide association analysis of HPAH. A burden of rare variants in BMPR2 significantly contributes to risk of the disease (p = 6.0 × 10−8). The approach robustly detected the gene having a large effect on the pathogenesis of PAH, despite the genetic heterogeneity. Eight probands in the nine families harbored possible pathogenic variants in BMPR2. Half of these variants were novel indels. One of the novel indels was a large 6.5 kilobase deletion spanning the entire region of exon 3. Another novel indel was a three base insertion (NM_001204.6:c.1277-10_1277-9insGGG) in intron 9 (Additional file 3: Figure S2). Although we could not dismiss that this insertion has no responsibility to the disease pathogenicity, a potential creation of a new splice acceptor site by this insertion was strongly suggested from the multiple splice site prediction tools (Table 1, Additional file 4: Figure S3 and Additional file 5: Figure S4) [31, 32]. The remaining one patient harbored a missense variant in KCNK3, which was the first replicative finding of channelopathy in Japanese population. Among the nine families, all variants identified here were mutually exclusive, suggesting that the variants have originated from independent genetic founder events.

Patients who suffered from chronic lung diseases such as chronic obstructive pulmonary disease (COPD) and pulmonary fibrosis are prone to pulmonary hypertension (PH) development. They are categorized as Group 3 in the latest guidelines [33]. Most patients with COPD develop mild PH but 3–5% of them show a further rise in mean pulmonary arterial pressure >35 mmHg. It is unknown how “severe PH-COPD”, formerly known as “out-of-proportion PH” is induced. Furthermore, the PAH-approved drugs are yet to be approved for the patients with Group 3 PH. In this study, a male patient (OM0195 in HPAH005) had been treated at another hospital for COPD. He was later diagnosed with PH and referred to our hospital. This patient could be categorized as “severe PH-COPD”, if none of his family members developed PAH. Since we had treated his daughter for IPAH at our hospital, we clinically diagnosed them as HPAH. Genetic testing revealed an in-frame-deletion (c.1443_1445delGAA) in BMPR2 in both patients. Underlying genetic predisposition might be one of the reasons for developing “severe PH-COPD”. Given our finding and a similarity of morphological appearance of vascular lesions between Group 1 and “severe PH-COPD” patients [33], PAH-approved drug treatments tailored to genetic diagnosis could well be a therapeutic strategy for such patients.

Conclusions

According to the genetic testing registry at the National Institutes of Health, the available panels for clinical genetic testing for PAH do not include KCNK3 and the detection methods are limited. Considering that pathogenic variants could occur within or spanning non-coding regions with a variety of sizes, the sequencing of the entire region of candidate genes is recommended to further understand the genetic factors relevant to PAH. This strategy will be essential for improving genetic diagnosis and counseling for PAH.

Additional files

Additional file 1: Table S1.(23K, pdf)

Summary of mapping statistics. (PDF 22 kb)

Additional file 2: Figure S1.(529K, pdf)

Large deletion found in the patient in family 8. (a) The large deletion was supported by depth of coverage and discordant paired-end reads. (b) The length of the deletion was approximately 6.5 kilobases and spanned the entire region of exon 3. Breakpoints of the large heterozygous deletion were located in the Alu repeat sequences, which were supported by read depths and paired-end reads spanning the long region. (PDF 529 kb)

Additional file 3: Figure S2.(412K, pdf)

Possible pathogenic variants found in the intron of BMPR2 in family 4. A rare three base insertion (GGG) at 10 bases upstream from exon 10 creates identical sequences around the canonical splicing site (ACAGGG). By the insertion, an out-of-frame protein could be translated due to the new splice acceptor site. (PDF 411 kb)

Additional file 4: Figure S3.(219K, pdf)

Creation of the new splice acceptor site by the three base insertion (NM_001204.6:c.1277-10_1277-9insGGG) predicted from in silico programs. (PDF 218 kb)

Additional file 5: Figure S4.(184K, pdf)

Disruption of wild type acceptor site by the single nucleotide variants (c.853-2A > G) predicted from in silico programs. (PDF 184 kb)

Acknowledgements

We wish to acknowledge the Nagahama City Office and the nonprofit organization “Zeroji Club” for their help in collecting control subjects. We thank Dr. Meiko Takahashi for critical reading of the manuscript.

Funding

This work was supported by the Practical Research Project for Rare/Intractable Diseases from the Ministry of Health, Labour and Welfare of Japan and Japan Agency for Medical Research and Development, AMED (15ek0109070h0002).

Availability of data and materials

The datasets generated and analyzed during the current study are available in the Human Genetic Variation Database repository; http://www.hgvd.genome.med.kyoto-u.ac.jp/repository/HGV0000004.html. http://www.hgvd.genome.med.kyoto-u.ac.jp/repository/HGV0000005.html.

Authors’ contributions

RY, HD, HM and FM designed the study. AO, CT, SK and HM collected the DNA and clinical information of pulmonary arterial hypertension patients. MS performed sequencing analysis. KH analyzed the data. KH, AO and FM wrote the manuscript. All authors critically reviewed the final version of the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Written informed consent for publication with data was obtained from patients and their family members.

Ethics approval and consent to participate

The study protocol and all subjects who participated in our study were approved by the Institutional Review Board of Kyoto University (approval number G755) and the National Hospital Organization Okayama Medical Center (approval number H24-RINKEN-24/H26-RINKEN-06) in which donors gave written informed consent in accordance with institutional and national guidelines.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Abbreviations

COPD
Chronic obstructive pulmonary disease
HPAH
Heritable PAH
IPAH
Idiopathic PAH
PAH
Pulmonary arterial hypertension

References

1. Deng Z, Morse JH, Slager SL, Cuervo N, Moore KJ, Venetos G, et al. Familial primary pulmonary hypertension (gene PPH1) is caused by mutations in the bone morphogenetic protein receptor-II gene. Am J Hum Genet. 2000;67:737–744. doi: 10.1086/303059. [PubMed] [Cross Ref]
2. Girerd B, Montani D, Jais X, Eyries M, Yaici A, Sztrymf B, et al. Genetic counselling in a national referral centre for pulmonary hypertension. Eur Respir J. 2016;47:541–552. doi: 10.1183/13993003.00717-2015. [PubMed] [Cross Ref]
3. Ma L, Chung WK. The genetic basis of pulmonary arterial hypertension. Hum Genet. 2014;133:471–479. doi: 10.1007/s00439-014-1419-3. [PubMed] [Cross Ref]
4. Machado RD, Aldred MA, James V, Harrison RE, Patel B, Schwalbe EC, et al. Mutations of the TGF-beta type II receptor BMPR2 in pulmonary arterial hypertension. Hum Mutat. 2006;27:121–132. doi: 10.1002/humu.20285. [PubMed] [Cross Ref]
5. 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:1113–1127. doi: 10.1002/humu.22904. [PMC free article] [PubMed] [Cross Ref]
6. Momose Y, Aimi Y, Hirayama T, Kataoka M, Ono M, Yoshino H, et al. De novo mutations in the BMPR2 gene in patients with heritable pulmonary arterial hypertension. Ann Hum Genet. 2015;79:85–91. doi: 10.1111/ahg.12096. [PubMed] [Cross Ref]
7. Trembath RC. Mutations in the TGF-beta type 1 receptor, ALK1, in combined primary pulmonary hypertension and hereditary haemorrhagic telangiectasia, implies pathway specificity. J Heart Lung Transplant. 2001;20:175. doi: 10.1016/S1053-2498(00)00352-1. [PubMed] [Cross Ref]
8. 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:446–448. doi: 10.1136/thx.2003.11890. [PMC free article] [PubMed] [Cross Ref]
9. Shintani M, Yagi H, Nakayama T, Saji T, Matsuoka R. A new nonsense mutation of SMAD8 associated with pulmonary arterial hypertension. J Med Genet. 2009;46:331–337. doi: 10.1136/jmg.2008.062703. [PubMed] [Cross Ref]
10. Austin ED, Ma L, LeDuc C, Berman Rosenzweig E, Borczuk A, Phillips JA, et al. Whole exome sequencing to identify a novel gene (caveolin-1) associated with human pulmonary arterial hypertension. Circ Cardiovasc Genet. 2012;5:336–343. doi: 10.1161/CIRCGENETICS.111.961888. [PMC free article] [PubMed] [Cross Ref]
11. Ma L, Roman-Campos D, Austin ED, Eyries M, Sampson KS, Soubrier F, et al. A novel channelopathy in pulmonary arterial hypertension. N Engl J Med. 2013;369:351–361. doi: 10.1056/NEJMoa1211097. [PMC free article] [PubMed] [Cross Ref]
12. Germain M, Eyries M, Montani D, Poirier O, Girerd B, Dorfmuller P, et al. Genome-wide association analysis identifies a susceptibility locus for pulmonary arterial hypertension. Nat Genet. 2013;45:518–521. doi: 10.1038/ng.2581. [PMC free article] [PubMed] [Cross Ref]
13. Wipff J, Dieude P, Guedj M, Ruiz B, Riemekasten G, Cracowski JL, et al. Association of a KCNA5 gene polymorphism with systemic sclerosis-associated pulmonary arterial hypertension in the European Caucasian population. Arthritis Rheum. 2010;62:3093–3100. doi: 10.1002/art.27607. [PubMed] [Cross Ref]
14. Kabata H, Satoh T, Kataoka M, Tamura Y, Ono T, Yamamoto M, et al. Bone morphogenetic protein receptor type 2 mutations, clinical phenotypes and outcomes of Japanese patients with sporadic or familial pulmonary hypertension. Respirology. 2013;18:1076–1082. [PubMed]
15. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754–1760. doi: 10.1093/bioinformatics/btp324. [PMC free article] [PubMed] [Cross Ref]
16. DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011;43:491–498. doi: 10.1038/ng.806. [PMC free article] [PubMed] [Cross Ref]
17. Higasa K, Miyake N, Yoshimura J, Okamura K, Niihori T, Saitsu H, et al. Human genetic variation database, a reference database of genetic variations in the Japanese population. J Hum Genet. 2016;61:547–553. doi: 10.1038/jhg.2016.12. [PMC free article] [PubMed] [Cross Ref]
18. Price AL, Kryukov GV, de Bakker PI, Purcell SM, Staples J, Wei LJ, et al. Pooled association tests for rare variants in exon-resequencing studies. Am J Hum Genet. 2010;86:832–838. doi: 10.1016/j.ajhg.2010.04.005. [PubMed] [Cross Ref]
19. Wang GT, Peng B, Leal SM. Variant association tools for quality control and analysis of large-scale sequence and genotyping array data. Am J Hum Genet. 2014;94:770–783. doi: 10.1016/j.ajhg.2014.04.004. [PubMed] [Cross Ref]
20. Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38:e164. doi: 10.1093/nar/gkq603. [PMC free article] [PubMed] [Cross Ref]
21. Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, et al. dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 2001;29:308–311. doi: 10.1093/nar/29.1.308. [PMC free article] [PubMed] [Cross Ref]
22. Clarke L, Zheng-Bradley X, Smith R, Kulesha E, Xiao C, Toneva I, et al. The 1000 genomes project: data management and community access. Nat Methods. 2012;9:459–462. doi: 10.1038/nmeth.1974. [PMC free article] [PubMed] [Cross Ref]
23. Reese MG, Moore B, Batchelor C, Salas F, Cunningham F, Marth GT, et al. A standard variation file format for human genome sequences. Genome Biol. 2010;11:R88. doi: 10.1186/gb-2010-11-8-r88. [PMC free article] [PubMed] [Cross Ref]
24. Tennessen JA, Bigham AW, O’Connor TD, Fu W, Kenny EE, Gravel S, et al. Evolution and functional impact of rare coding variation from deep sequencing of human exomes. Science. 2012;337:64–69. doi: 10.1126/science.1219240. [PMC free article] [PubMed] [Cross Ref]
25. Landrum MJ, Lee JM, Benson M, Brown G, Chao C, Chitipiralla S, et al. ClinVar: public archive of interpretations of clinically relevant variants. Nucleic Acids Res. 2016;44:D862–868. doi: 10.1093/nar/gkv1222. [PMC free article] [PubMed] [Cross Ref]
26. Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7:248–249. doi: 10.1038/nmeth0410-248. [PMC free article] [PubMed] [Cross Ref]
27. Schwarz JM, Cooper DN, Schuelke M, Seelow D. MutationTaster2: mutation prediction for the deep-sequencing age. Nat Methods. 2014;11:361–362. doi: 10.1038/nmeth.2890. [PubMed] [Cross Ref]
28. Chun S, Fay JC. Identification of deleterious mutations within three human genomes. Genome Res. 2009;19:1553–1561. doi: 10.1101/gr.092619.109. [PubMed] [Cross Ref]
29. Siepel A, Pollard K, Haussler D, Apostolico A, Guerra C, Istrail S, et al. New methods for detecting lineage-specific selection. Res Comput Mol Biol Proc. 2006;3909:190–205. doi: 10.1007/11732990_17. [Cross Ref]
30. Liu X, Jian X, Boerwinkle E. dbNSFP v2.0: a database of human non-synonymous SNVs and their functional predictions and annotations. Hum Mutat. 2013;34:E2393–2402. doi: 10.1002/humu.22376. [PMC free article] [PubMed] [Cross Ref]
31. Yeo G, Burge CB. Maximum entropy modeling of short sequence motifs with applications to RNA splicing signals. J Comput Biol. 2004;11:377–394. doi: 10.1089/1066527041410418. [PubMed] [Cross Ref]
32. Desmet FO, Hamroun D, Lalande M, Collod-Beroud G, Claustres M, Beroud C. Human splicing finder: an online bioinformatics tool to predict splicing signals. Nucleic Acids Res. 2009;37:e67. doi: 10.1093/nar/gkp215. [PMC free article] [PubMed] [Cross Ref]
33. Seeger W, Adir Y, Barbera JA, Champion H, Coghlan JG, Cottin V, et al. Pulmonary hypertension in chronic lung diseases. J Am Coll Cardiol. 2013;62:D109–116. doi: 10.1016/j.jacc.2013.10.036. [PubMed] [Cross Ref]

Articles from BMC Pulmonary Medicine are provided here courtesy of BioMed Central