Genetic Variants Associated With Uncontrolled Blood Pressure on Thiazide Diuretic/β-Blocker Combination Therapy in the PEAR (Pharmacogenomic Evaluation of Antihypertensive Responses) and INVEST (International Verapamil-SR Trandolapril Study) Trials

Oyunbileg Magvanjav, Yan Gong, Caitrin W McDonough, Arlene B Chapman, Stephen T Turner, John G Gums, Kent R Bailey, Eric Boerwinkle, Amber L Beitelshees, Toshihiro Tanaka, Michiaki Kubo, Carl J Pepine, Rhonda M Cooper-DeHoff, Julie A Johnson, Oyunbileg Magvanjav, Yan Gong, Caitrin W McDonough, Arlene B Chapman, Stephen T Turner, John G Gums, Kent R Bailey, Eric Boerwinkle, Amber L Beitelshees, Toshihiro Tanaka, Michiaki Kubo, Carl J Pepine, Rhonda M Cooper-DeHoff, Julie A Johnson

Abstract

Background: The majority of hypertensive individuals require combination antihypertensive therapy to achieve adequate blood pressure (BP) control. This study aimed to identify genetic variants associated with uncontrolled BP on combination therapy with a thiazide diuretic and a β-blocker.

Methods and results: A genome-wide association study of uncontrolled BP on combination therapy was conducted among 314 white participants of the PEAR (Pharmacogenomic Evaluation of Antihypertensive Responses) trial. Multivariable logistic regression analysis was used. Genetic variants meeting a suggestive level of significance (P<1.0E-05) were tested for replication in an external cohort, INVEST (International Verapamil-SR Trandolapril study). We also examined genome-wide variant associations with systolic and diastolic BP response on combination therapy and tested for replication. We discovered a single nucleotide polymorphism, the rs261316 major allele, at chromosome 15 in the gene ALDH1A2 associated with an increased odds of having uncontrolled BP on combination therapy (odds ratio: 2.56, 95% confidence interval, 1.69-3.88, P=8.64E-06). This single nucleotide polymorphism replicated (odds ratio: 1.86, 95% confidence interval, 1.35-2.57, P=0.001) and approached genome-wide significance in the meta-analysis between discovery and replication cohorts (odds ratio: 2.16, 95% confidence interval, 1.63-2.86, P=8.60E-08). Other genes in the region surrounding rs261316 (ALDH1A2) include AQP9 and LIPC.

Conclusions: A single nucleotide polymorphism in the gene ALDH1A2 may be associated with uncontrolled BP following treatment with a thiazide diuretic/β-blocker combination.

Clinical trial registration: URL: https://www.clinicaltrials.gov. Unique identifier: NCT00246519.

Keywords: combination therapy; genomics; high blood pressure; hypertension; pharmacogenomics; thiazide diuretics; β‐blockers.

© 2017 The Authors. Published on behalf of the American Heart Association, Inc., by Wiley.

Figures

Figure 1
Figure 1
Association of rs261316 allele (T) in ALDH1A2 with uncontrolled BP on combination therapy in PEAR with replication analysis in INVEST. P values and 95% confidence intervals for INVEST cohorts are 1‐sided, consistent with a 1‐sided hypothesis for replication; sample sizes: PEAR whites (n=314), INVEST whites (n=221), and INVEST Hispanics (n=193). BP indicates blood pressure; INVEST, International Verapamil‐SR Trandolapril Study; PEAR, Pharmacogenomic Evaluation of Antihypertensive Responses.
Figure 2
Figure 2
Regional plot of top genetic variant associated with uncontrolled BP on combination therapy in PEAR with replication in INVEST. rs261316 associated with uncontrolled BP on combination therapy (PEAR whites: P=8.64E‐06; INVEST white–Hispanic meta‐analysis: 1‐sided P=0.001; PEAR and INVEST meta‐analysis: P=8.60E‐08); BP indicates blood pressure; INVEST, International Verapamil‐SR Trandolapril Study; PEAR, Pharmacogenomic Evaluation of Antihypertensive Responses.

References

    1. Mills KT, Bundy JD, Kelly TN, Reed JE, Kearney PM, Reynolds K, Chen J, He J. Global disparities of hypertension prevalence and control: a systematic analysis of population‐based studies from 90 countries. Circulation. 2016;134:441–450.
    1. Writing Group Members , Mozaffarian D, Benjamin EJ, Go AS, Arnett DK, Blaha MJ, Cushman M, Das SR, de Ferranti S, Després JP, Fullerton HJ, Howard VJ, Huffman MD, Isasi CR, Jiménez MC, Judd SE, Kissela BM, Lichtman JH, Lisabeth LD, Liu S, Mackey RH, Magid DJ, McGuire DK, Mohler ER, Moy CS, Muntner P, Mussolino ME, Nasir K, Neumar RW, Nichol G, Palaniappan L, Pandey DK, Reeves MJ, Rodriguez CJ, Rosamond W, Sorlie PD, Stein J, Towfighi A, Turan TN, Virani SS, Woo D, Yeh RW, Turner MB, American Heart Association Statistics Committee; Stroke Statistics Subcommittee. Heart Disease and Stroke Statistics‐2016 Update: a Report From the American Heart Association. Circulation. 2016;133:e38–e60.
    1. Bakris G, Sarafidis P, Agarwal R, Ruilope L. Review of blood pressure control rates and outcomes. J Am Soc Hypertens. 2014;8:127–141.
    1. James PA, Oparil S, Carter BL, Cushman WC, Dennison‐Himmelfarb C, Handler J, Lackland DT, LeFevre ML, MacKenzie TD, Ogedegbe O, Smith SC, Svetkey LP, Taler SJ, Townsend RR, Wright JT, Narva AS, Ortiz E. 2014 evidence‐based guideline for the management of high blood pressure in adults: report from the panel members appointed to the Eighth Joint National Committee (JNC 8). JAMA. 2014;311:507–520.
    1. Jamerson K, Weber MA, Bakris GL, Dahlöf B, Pitt B, Shi V, Hester A, Gupte J, Gatlin M, Velazquez EJ; ACCOMPLISH Trial Investigators . Benazepril plus amlodipine or hydrochlorothiazide for hypertension in high‐risk patients. N Engl J Med. 2008;359:2417–2428.
    1. Dahlöf B, Sever PS, Poulter NR, Wedel H, Beevers DG, Caulfield M, Collins R, Kjeldsen SE, Kristinsson A, McInnes GT, Mehlsen J, Nieminen M, O'Brien E, Ostergren J; ASCOT Investigators . Prevention of cardiovascular events with an antihypertensive regimen of amlodipine adding perindopril as required versus atenolol adding bendroflumethiazide as required, in the Anglo‐Scandinavian Cardiac Outcomes Trial‐Blood Pressure Lowering Arm (ASCOT‐BPLA): a multicentre randomised controlled trial. Lancet. 2005;366:895–906.
    1. Bakris GL, Cooper‐Dehoff RM, Zhou Q, Kupfer S, Champion A, Pepine CJ; INVEST Investigators . Dual therapy in hypertensive patients with coronary artery disease: the role of calcium channel blockers and beta‐blockers. Am J Cardiovasc Drugs. 2007;7(suppl 1):25–29.
    1. Chowdhury R, Khan H, Heydon E, Shroufi A, Fahimi S, Moore C, Stricker B, Mendis S, Hofman A, Mant J, Franco OH. Adherence to cardiovascular therapy: a meta‐analysis of prevalence and clinical consequences. Eur Heart J. 2013;34:2940–2948.
    1. Cooper‐DeHoff RM, Johnson JA. Hypertension pharmacogenomics: in search of personalized treatment approaches. Nat Rev Nephrol. 2016;12:110–122.
    1. Brunner M, Cooper‐DeHoff RM, Gong Y, Karnes JH, Langaee TY, Pepine CJ, Johnson JA; IINVEST Investigators . Factors influencing blood pressure response to trandolapril add‐on therapy in patients taking verapamil SR (from the International Verapamil SR/Trandolapril [INVEST] Study). Am J Cardiol. 2007;99:1549–1554.
    1. Laffer CL, Elijovich F, Eckert GJ, Tu W, Pratt JH, Brown NJ. Genetic variation in CYP4A11 and blood pressure response to mineralocorticoid receptor antagonism or ENaC inhibition: an exploratory pilot study in African Americans. J Am Soc Hypertens. 2014;8:475–480.
    1. Gupta S, Chattopadhyaya I, Agrawal BK, Sehajpal PK, Goel RK. Correlation of renin angiotensin system (RAS) candidate gene polymorphisms with response to Ramipril in patients with essential hypertension. J Postgrad Med. 2015;61:21–26.
    1. Dumitrescu L, Ritchie MD, Denny JC, El Rouby NM, McDonough CW, Bradford Y, Ramirez AH, Bielinski SJ, Basford MA, Chai HS, Peissig P, Carrell D, Pathak J, Rasmussen LV, Wang X, Pacheco JA, Kho AN, Hayes MG, Matsumoto M, Smith ME, Li R, Cooper‐DeHoff RM, Kullo IJ, Chute CG, Chisholm RL, Jarvik GP, Larson EB, Carey D, McCarty CA, Williams MS, Roden DM, Bottinger E, Johnson JA, de Andrade M, Crawford DC. Genome‐wide study of resistant hypertension identified from electronic health records. PLoS One. 2017;12:e0171745.
    1. Gong Y, McDonough CW, Beitelshees AL, El Rouby N, Hiltunen TP, O'Connell JR, Padmanabhan S, Langaee TY, Hall K, Schmidt SO, Curry RW, Gums JG, Donner KM, Kontula KK, Bailey KR, Boerwinkle E, Takahashi A, Tanaka T, Kubo M, Chapman AB, Turner ST, Pepine CJ, Cooper‐DeHoff RM, Johnson JA. PTPRD gene associated with blood pressure response to atenolol and resistant hypertension. J Hypertens. 2015;33:2278–2285.
    1. Fontana V, McDonough CW, Gong Y, El Rouby NM, Sá AC, Taylor KD, Chen YD, Gums JG, Chapman AB, Turner ST, Pepine CJ, Johnson JA, Cooper‐DeHoff RM. Large‐scale gene‐centric analysis identifies polymorphisms for resistant hypertension. J Am Heart Assoc. 2014;3:e001398 DOI: .
    1. Padmanabhan S, Aman A, Dominiczak AF. Genomics of hypertension. Pharmacol Res. 2017;121:219–229.
    1. Padmanabhan S, Caulfield M, Dominiczak AF. Genetic and molecular aspects of hypertension. Circ Res. 2015;116:937–959.
    1. Johnson JA, Boerwinkle E, Zineh I, Chapman AB, Bailey K, Cooper‐DeHoff RM, Gums J, Curry RW, Gong Y, Beitelshees AL, Schwartz G, Turner ST. Pharmacogenomics of antihypertensive drugs: rationale and design of the Pharmacogenomic Evaluation of Antihypertensive Responses (PEAR) study. Am Heart J. 2009;157:442–449.
    1. Pepine CJ, Handberg‐Thurmond E, Marks RG, Conlon M, Cooper‐DeHoff R, Volkers P, Zellig P. Rationale and design of the International Verapamil SR/Trandolapril Study (INVEST): an Internet‐based randomized trial in coronary artery disease patients with hypertension. J Am Coll Cardiol. 1998;32:1228–1237.
    1. Beitelshees AL, Gong Y, Wang D, Schork NJ, Cooper‐Dehoff RM, Langaee TY, Shriver MD, Sadee W, Knot HJ, Pepine CJ, Johnson JA; INVEST Investigators I . KCNMB1 genotype influences response to verapamil SR and adverse outcomes in the INternational VErapamil SR/Trandolapril STudy (INVEST). Pharmacogenet Genomics. 2007;17:719–729.
    1. Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM, Korbel JO, Marchini JL, McCarthy S, McVean GA, Abecasis GR, Consortium GP. A global reference for human genetic variation. Nature. 2015;526:68–74.
    1. Purcell S, Neale B, Todd‐Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC. PLINK: a tool set for whole‐genome association and population‐based linkage analyses Am J Hum Genet. 2007; 81:559–575.
    1. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D. Principal components analysis corrects for stratification in genome‐wide association studies. Nat Genet. 2006;38:904–909.
    1. Li Y, Willer CJ, Ding J, Scheet P, Abecasis GR. MaCH: using sequence and genotype data to estimate haplotypes and unobserved genotypes. Genet Epidemiol. 2010;34:816–834.
    1. Fuchsberger C, Abecasis GR, Hinds DA. minimac2: faster genotype imputation. Bioinformatics. 2015;31:782–784.
    1. Howie B, Fuchsberger C, Stephens M, Marchini J, Abecasis GR. Fast and accurate genotype imputation in genome‐wide association studies through pre‐phasing. Nat Genet. 2012;44:955–959.
    1. Aulchenko YS, Struchalin MV, van Duijn CM. ProbABEL package for genome‐wide association analysis of imputed data. BMC Bioinformatics. 2010;11:134.
    1. Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, Gliedt TP, Boehnke M, Abecasis GR, Willer CJ. LocusZoom: regional visualization of genome‐wide association scan results. Bioinformatics. 2010;26:2336–2337.
    1. Ward LD, Kellis M. HaploReg: a resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants. Nucleic Acids Res. 2012;40:D930–D934.
    1. Consortium G. Human genomics. The Genotype‐Tissue Expression (GTEx) pilot analysis: multitissue gene regulation in humans. Science. 2015;348:648–660.
    1. Mallipattu SK, He JC. The beneficial role of retinoids in glomerular disease. Front Med (Lausanne). 2015;2:16.
    1. Xavier‐Neto J, Sousa Costa AM, Figueira AC, Caiaffa CD, Amaral FN, Peres LM, da Silva BS, Santos LN, Moise AR, Castillo HA. Signaling through retinoic acid receptors in cardiac development: doing the right things at the right times. Biochem Biophys Acta. 2015;1849:94–111.
    1. El Kares R, Manolescu DC, Lakhal‐Chaieb L, Montpetit A, Zhang Z, Bhat PV, Goodyer P. A human ALDH1A2 gene variant is associated with increased newborn kidney size and serum retinoic acid. Kidney Int. 2010;78:96–102.
    1. Kanzaki G, Tsuboi N, Haruhara K, Koike K, Ogura M, Shimizu A, Yokoo T. Factors associated with a vicious cycle involving a low nephron number, hypertension and chronic kidney disease. Hypertens Res. 2015;38:633–641.
    1. Adeyemo A, Gerry N, Chen G, Herbert A, Doumatey A, Huang H, Zhou J, Lashley K, Chen Y, Christman M, Rotimi C. A genome‐wide association study of hypertension and blood pressure in African Americans. PLoS Genet. 2009;5:e1000564.
    1. Hiltunen TP, Donner KM, Sarin AP, Saarela J, Ripatti S, Chapman AB, Gums JG, Gong Y, Cooper‐DeHoff RM, Frau F, Glorioso V, Zaninello R, Salvi E, Glorioso N, Boerwinkle E, Turner ST, Johnson JA, Kontula KK. Pharmacogenomics of hypertension: a genome‐wide, placebo‐controlled cross‐over study, using four classes of antihypertensive drugs. J Am Heart Assoc. 2015;4:e001521 DOI: .
    1. Zhou TB, Ou C, Rong L, Drummen GP. Effect of all‐trans retinoic acid treatment on prohibitin and renin‐angiotensin‐aldosterone system expression in hypoxia‐induced renal tubular epithelial cell injury. J Renin Angiotensin Aldosterone Syst. 2014;15:243–249.
    1. Zhou TB, Wu WF, Qin YH, Yin SS. Association of all‐trans retinoic acid treatment with the renin‐angiotensin aldosterone system expression in glomerulosclerosis rats. J Renin Angiotensin Aldosterone Syst. 2013;14:299–307.
    1. Zhong JC, Huang DY, Yang YM, Li YF, Liu GF, Song XH, Du K. Upregulation of angiotensin‐converting enzyme 2 by all‐trans retinoic acid in spontaneously hypertensive rats. Hypertension. 2004;44:907–912.
    1. Dechow C, Morath C, Peters J, Lehrke I, Waldherr R, Haxsen V, Ritz E, Wagner J. Effects of all‐trans retinoic acid on renin‐angiotensin system in rats with experimental nephritis. Am J Physiol Renal Physiol. 2001;281:F909–F919.
    1. Todorov VT. PPARgamma‐dependent control of renin expression: molecular mechanisms and pathophysiological relevance. PPAR Res. 2013;2013:451016.
    1. Stump M, Mukohda M, Hu C, Sigmund CD. PPARγ regulation in hypertension and metabolic syndrome. Curr Hypertens Rep. 2015;17:89.
    1. Kumar P, Periyasamy R, Das S, Neerukonda S, Mani I, Pandey KN. All‐trans retinoic acid and sodium butyrate enhance natriuretic peptide receptor a gene transcription: role of histone modification. Mol Pharmacol. 2014;85:946–957.
    1. Wang Y, Han Y, Yang J, Wang Z, Liu L, Wang W, Zhou L, Wang D, Tan X, Fu C, Jose PA, Zeng C. Relaxant effect of all‐trans‐retinoic acid via NO‐sGC‐cGMP pathway and calcium‐activated potassium channels in rat mesenteric artery. Am J Physiol Heart Circ Physiol. 2013;304:H51–H57.
    1. Zhou W, Lin J, Chen H, Wang J, Liu Y, Xia M. Retinoic acid induces macrophage cholesterol efflux and inhibits atherosclerotic plaque formation in apoE‐deficient mice. Br J Nutr. 2015;114:509–518.
    1. Bechor S, Zolberg, Relevy N, Harari A, Almog T, Kamari Y, Ben‐Amotz A, Harats D, Shaish A. 9‐cis β‐carotene increased cholesterol efflux to HDL in macrophages. Nutrients. 2016;8: 435.
    1. Liu Y, Chen H, Mu D, Li D, Zhong Y, Jiang N, Zhang Y, Xia M. Association of Serum Retinoic Acid With Risk of Mortality in Patients With Coronary Artery Disease. Circ Res. 2016;119:557–563.
    1. Kobayashi J, Miyashita K, Nakajima K, Mabuchi H. Hepatic lipase: a comprehensive view of its role on plasma lipid and lipoprotein metabolism. J Atheroscler Thromb. 2015;22:1001–1011.
    1. Brunzell JD, Zambon A, Deeb SS. The effect of hepatic lipase on coronary artery disease in humans is influenced by the underlying lipoprotein phenotype. Biochem Biophys Acta. 2012;1821:365–372.
    1. Aggarwal A, Srivastava S, Velmurugan M. Newer perspectives of coronary artery disease in young. World J Cardiol. 2016;8:728–734.
    1. Ríos‐González BE, Ibarra‐Cortés B, Ramírez‐López G, Sánchez‐Corona J, Magaña‐Torres MT. Association of polymorphisms of genes involved in lipid metabolism with blood pressure and lipid values in mexican hypertensive individuals. Dis Markers. 2014;2014:150358.
    1. Kokubo Y, Tomoike H, Tanaka C, Banno M, Okuda T, Inamoto N, Kamide K, Kawano Y, Miyata T. Association of sixty‐one non‐synonymous polymorphisms in forty‐one hypertension candidate genes with blood pressure variation and hypertension. Hypertens Res. 2006;29:611–619.
    1. Damiano AE. Review: water channel proteins in the human placenta and fetal membranes. Placenta. 2011;32(suppl 2):S207–S211.
    1. Madeira A, Moura TF, Soveral G. Aquaglyceroporins: implications in adipose biology and obesity. Cell Mol Life Sci. 2015;72:759–771.
    1. Procino G, Romano F, Torielli L, Ferrari P, Bianchi G, Svelto M, Valenti G. Altered expression of renal aquaporins and α‐adducin polymorphisms may contribute to the establishment of salt‐sensitive hypertension. Am J Hypertens. 2011;24:822–828.
    1. Montiel V, Leon Gomez E, Bouzin C, Esfahani H, Romero Perez M, Lobysheva I, Devuyst O, Dessy C, Balligand JL. Genetic deletion of aquaporin‐1 results in microcardia and low blood pressure in mouse with intact nitric oxide‐dependent relaxation, but enhanced prostanoids‐dependent relaxation. Pflugers Arch. 2014;466:237–251.
    1. Pluznick JL. Renal and cardiovascular sensory receptors and blood pressure regulation. Am J Physiol Renal Physiol. 2013;305:F439–F444.
    1. Miyamoto J, Kasubuchi M, Nakajima A, Irie J, Itoh H, Kimura I. The role of short‐chain fatty acid on blood pressure regulation. Curr Opin Nephrol Hypertens. 2016;25:379–383.
    1. Paneni F, Costantino S, Battista R, Castello L, Capretti G, Chiandotto S, Scavone G, Villano A, Pitocco D, Lanza G, Volpe M, Lüscher TF, Cosentino F. Adverse epigenetic signatures by histone methyltransferase Set7 contribute to vascular dysfunction in patients with type 2 diabetes mellitus. Circ Cardiovasc Genet. 2015;8:150–158.
    1. Tragante V, Barnes MR, Ganesh SK, Lanktree MB, Guo W, Franceschini N, Smith EN, Johnson T, Holmes MV, Padmanabhan S, Karczewski KJ, Almoguera B, Barnard J, Baumert J, Chang YP, Elbers CC, Farrall M, Fischer ME, Gaunt TR, Gho JM, Gieger C, Goel A, Gong Y, Isaacs A, Kleber ME, Mateo Leach I, McDonough CW, Meijs MF, Melander O, Nelson CP, Nolte IM, Pankratz N, Price TS, Shaffer J, Shah S, Tomaszewski M, van der Most PJ, VanIperen EP, Vonk JM, Witkowska K, Wong CO, Zhang L, Beitelshees AL, Berenson GS, Bhatt DL, Brown M, Burt A, Cooper‐DeHoff RM, Connell JM, Cruickshanks KJ, Curtis SP, Davey‐Smith G, Delles C, Gansevoort RT, Guo X, Haiqing S, Hastie CE, Hofker MH, Hovingh GK, Kim DS, Kirkland SA, Klein BE, Klein R, Li YR, Maiwald S, Newton‐Cheh C, O'Brien ET, Onland‐Moret NC, Palmas W, Parsa A, Penninx BW, Pettinger M, Vasan RS, Ranchalis JE, M Ridker P, Rose LM, Sever P, Shimbo D, Steele L, Stolk RP, Thorand B, Trip MD, van Duijn CM, Verschuren WM, Wijmenga C, Wyatt S, Young JH, Zwinderman AH, Bezzina CR, Boerwinkle E, Casas JP, Caulfield MJ, Chakravarti A, Chasman DI, Davidson KW, Doevendans PA, Dominiczak AF, FitzGerald GA, Gums JG, Fornage M, Hakonarson H, Halder I, Hillege HL, Illig T, Jarvik GP, Johnson JA, Kastelein JJ, Koenig W, Kumari M, März W, Murray SS, O'Connell JR, Oldehinkel AJ, Pankow JS, Rader DJ, Redline S, Reilly MP, Schadt EE, Kottke‐Marchant K, Snieder H, Snyder M, Stanton AV, Tobin MD, Uitterlinden AG, van der Harst P, van der Schouw YT, Samani NJ, Watkins H, Johnson AD, Reiner AP, Zhu X, de Bakker PI, Levy D, Asselbergs FW, Munroe PB, Keating BJ. Gene‐centric meta‐analysis in 87,736 individuals of European ancestry identifies multiple blood‐pressure‐related loci. Am J Hum Genet. 2014;94:349–360.
    1. Cooper‐DeHoff R, Handberg E, Heissenberg C, Johnson K. Electronic prescribing via the internet for a coronary artery disease and hypertension megatrial. Clin Cardiol. 2001;24:V14–V16.
    1. Pepine CJ, Handberg EM, Cooper‐DeHoff RM, Marks RG, Kowey P, Messerli FH, Mancia G, Cangiano JL, Garcia‐Barreto D, Keltai M, Erdine S, Bristol HA, Kolb HR, Bakris GL, Cohen JD, Parmley WW; Investigators I . A calcium antagonist vs a non‐calcium antagonist hypertension treatment strategy for patients with coronary artery disease. The International Verapamil‐Trandolapril Study (INVEST): a randomized controlled trial. JAMA. 2003;290:2805–2816.
    1. Maranville JC, Cox NJ. Pharmacogenomic variants have larger effect sizes than genetic variants associated with other dichotomous complex traits. Pharmacogenomics J. 2016;16:388–392.

Source: PubMed

3
订阅