- Open Access
Lack of association of the HMGA1 IVS5-13insC variant with type 2 diabetes in an ethnically diverse hypertensive case control cohort
Journal of Translational Medicine volume 11, Article number: 12 (2013)
Recently, the high-mobility group A1 gene (HMGA1) variant IVS5-13insC has been associated with type 2 diabetes, but reported associations are inconsistent and data are lacking in Hispanic and African American populations. We sought to investigate the HMGA1-diabetes association and to characterize IVS5-13insC allele frequencies and linkage disequilibrium (LD) in 3,070 Caucasian, Hispanic, and African American patients from the INternational VErapamil SR-Trandolapril STudy (INVEST).
INVEST was a randomized, multicenter trial comparing two antihypertensive treatment strategies in an ethnically diverse cohort of hypertensive, coronary artery disease patients. Controls, who were diabetes-free throughout the study, and type 2 diabetes cases, either prevalent or incident, were genotyped for IVS5-13insC using Taqman®, confirmed with Pyrosequencing and Sanger sequencing. For LD analysis, genotyping for eight additional HMGA1 single nucleotide polymorphisms (SNPs) was performed using the Illumina® HumanCVD BeadChip. We used logistic regression to test association of the HMGA1 IVS5-13insC and diabetes, adjusted for age, gender, body mass index, and percentage European, African, and Native American ancestry.
We observed IVS5-13insC minor allele frequencies consistent with previous literature in Caucasians and African Americans (0.03 in cases and 0.04 in controls for both race/ethnic groups), and higher frequencies in Hispanics (0.07 in cases and 0.07 in controls). The IVS5-13insC was not associated with type 2 diabetes overall (odds ratio 0.98 [0.76-1.26], p=0.88) or in any race/ethnic group. Pairwise LD (r2) of IVS5-13insC and rs9394200, a SNP previously used as a tag SNP for IVS5-13insC, was low (r2=0.47 in Caucasians, r2=0.25 in Hispanics, and r2=0.06 in African Americans). Furthermore, in silico analysis suggested a lack of functional consequences for the IVS5-13insC variant.
Our results suggest that IVS5-13insC is not a functional variant and not associated with type 2 diabetes in an ethnically diverse, hypertensive, coronary artery disease population. Larger, more adequately powered studies need to be performed to confirm our findings.
Type 2 diabetes constitutes a major and growing health problem worldwide and is predicted to afflict 490 million by 2030 . Type 2 diabetes has strong genetic influences and many polymorphisms have now been reproducibly associated with type 2 diabetes [2, 3]. However, genome wide association studies (GWAS) explain only 10-15% of heritability and have not consistently improved diabetes risk prediction . Low frequency variation may account for much of the missing heritability in type 2 diabetes risk and may help translate genetic association study results into clinical type 2 diabetes risk prediction.
Recently, the low frequency insertion polymorphism IVS5-13insC (c.136-14_136-13insC) in the high-mobility group A1 gene (HMGA1), a transcriptional regulator of the insulin receptor gene (INSR), was identified and associated with type 2 diabetes . Whereas two studies have observed a significant positive association between IVS5-13insC and type 2 diabetes in Caucasian and Chinese populations [5, 6], another study in Caucasians observed no association . Data are lacking for IVS5-13insC in populations with African and Hispanic descent, which have disproportionately high type 2 diabetes prevalence . Conflicting results for an HMGA1 association with type 2 diabetes and the lack of data in diverse race/ethnic groups make clinical translation of the HMGA1 IVS5-13insC genotyping especially difficult.
Evidence for the functional impact of the HMGA1 IVS5-13insC variant is also conflicting. One study observed that HMGA1 and INSR expression was decreased in diabetic carriers of IVS5-13insC versus wild type diabetic and non-diabetic patients . Additionally, INSR protein expression and insulin-binding capacity was restored in lymphoblasts obtained from diabetic IVS5-13insC carriers by HMGA1 DNA transfection. Another study observed no effect of IVS5-13insC on HMGA1 or INSR expression in adipose tissue of normoglycemic patients . IVS5-13insC occurs at position −13 of HMGA1 exon 6, but the direct mechanism of the variant’s effects on mRNA expression or amino acid sequence remains unclear.
We tested the association of HMGA1 IVS5-13insC with type 2 diabetes in an ethnically diverse population from the INternational VErapamil SR-Trandolapril STudy (INVEST). INVEST compared CV outcomes and NOD in hypertensive coronary artery disease patients treated with two antihypertensive treatment strategies. We also determined minor allele frequencies (MAF) and linkage disequilibrium (LD) for HMGA1 variants and tested the functional impact of IVS5-13insC in silico.
Study design and participants
INVEST compared CV outcomes and incident diabetes in hypertensive, coronary artery disease patients at least 50 years of age during randomized treatment with either an atenolol-based or a verapamil sustained release (SR)-based antihypertensive treatment strategy. The design, primary outcome, and NOD results have been previously published in detail [9–12]. Briefly, the verapamil SR strategy consisted of stepped therapy with verapamil SR, trandolapril add-on, dose titration, then HCTZ add-on treatment for BP control and end organ protection as necessary. The atenolol-based strategy consisted of atenolol, HCTZ add-on, dose titration, then trandolapril add-on treatment as necessary. The INVEST GENEtic Substudy (INVEST-GENES) collected DNA samples from 5,979 INVEST patients at 187 sites in the United States and Puerto Rico.
We conducted a nested case control study including cases with type 2 diabetes at baseline (prevalent diabetes) or type 2 diabetes that developed during a mean 2.8 years follow-up (incident diabetes). Type 2 diabetes was determined by patient report and by site investigators from a review of all available patient data, including use of diabetic medication and lab measures . We identified age, gender, and race/ethnicity-matched controls who remained diabetes-free over a mean 2.8 years follow-up. Age matching was performed after stratification by decade and we attempted to match cases and controls in a 1:1 ratio. The institutional review boards of participating study centers approved the study protocol and all patients provided written informed consent for participation in INVEST and additional written informed consent for genetic studies. INVEST is registered at clinicaltrials.gov (NCT00133692).
Genotyping for HMGA1 IVS5-13insC was performed using TaqMan® (Applied Biosystems, Foster City, CA, USA) with PCR primers and probe for IVS5-13insC (PN4331349) purchased from Applied Biosystems. For IVS5-13insC genotyping quality control, 5% of samples were genotyped in duplicate on Taqman®. A total of 612 Taqman genotypes were confirmed using pyrosequencing (Biotage AB, Uppsala, Sweden), using the following PCR and sequencing primers respectively: forward-biotinylated-5′-GGGGTGGAAACAGGTGATG-3′, reverse-5′-CACTTCGCTGGGCTCCTT-3′, and reverse-5′-TTCTGTAAAGACAGAGG-3′. Sanger sequencing was used to genotype 58 samples that showed discrepancies between Taqman® and Pyrosequencing platforms.
Genotyping for eight additional HMGA1 single nucleotide polymorphisms (SNPs) was performed using the HumanCVD BeadChip and Infinium II Assay (Illumina, San Diego, CA) on 1,489 INVEST patients to perform LD analyses. The HumanCVD BeadChip contains approximately 50,000 cosmopolitan tag SNPs for 2,100 CV and metabolic-related genes . HumanCVD BeadChip data quality was ensured in PLINK using genotype and sample call rates, concordance rates for blind duplicates, gender confirmation, cryptic relatedness using pairwise identity-by-descent, and estimation of heterozygosity using the inbreeding coefficient F . Individuals were excluded if call rates were below 90 percent and SNPs were excluded if call rates were below 95 percent. In addition, 87 ancestry informative markers were genotyped in 2,860 INVEST patients to estimate Caucasian, African, and Native American ancestry using STRUCTURE . Race/ethnic groups were determined by patient self-report with interaction by the study investigator  and confirmed using principal components analysis generated from LD-pruned HumanCVD BeadChip data and ancestry informative markers.
All statistical analyses were performed using SAS version 9.2 (SAS, Cary, NC). Differences in patient characteristics comparing cases and controls at baseline were determined using t-tests and chi square tests, as appropriate. Deviations from Hardy Weinberg Equilibrium were assessed using a chi square test. Multi-variable logistic regressions were performed overall and by race/ethnic group to determine odds ratios (ORs) and 95% confidence intervals (95%CIs) for type 2 diabetes in IVS5-13insC variant carriers versus non-variant carriers. Variables for adjustment included age, gender, and body mass index (BMI) in order to maintain consistency with previously published HMGA1 analyses [5, 7]. Percentage of Caucasian, African, and Native American ancestry, as estimated by ancestry informative markers, was included as a variable for adjustment in analyses where all race/ethnic groups were combined.
The combined race/ethnic group analysis was considered primary with alpha=0.05, since IVS5-13insC was presumed to be a functional SNP with similar consequences across race/ethnic groups, based on previously published reports [5, 7]. Assuming a MAF of 0.05 and OR of 1.40, we had 84% power to detect an association between the IVS5-13insC variant and type 2 diabetes in a dominant model in the overall population. Assuming a MAF of 0.05 and OR of 1.40, we had 61% power in Hispanics, 44% power in whites, and 18% power in blacks to detect an association by race/ethnicity. Pairwise LD values (r2) and LD plots were generated by race/ethnic group using Haploview . We predicted in silico functional consequences of the IVS5-13insC variant using SNPNexus  and ESEfinder 3.0 .
We identified 446 incident type 2 diabetes cases over a mean 2.8 years of follow-up and genotyped an additional 1329 prevalent type 2 diabetes cases in INVEST-GENES. At baseline, patients with incident or prevalent type 2 diabetes had higher BMI and lower diastolic BP versus age, race/ethnicity, and gender-matched controls (Table 1). Cases also had a higher prevalence of hypercholesterolemia, left ventricular hypertrophy, and congestive heart failure.
The IVS5-13insC variant did not deviate from HWE in any race/ethnic group (Table 2). All HMGA1 SNPs had call rates above 95 percent and no SNPs failed quality control based on call rate, Hardy Weinberg equilibrium testing, concordance with blind duplicates, or estimation of heterozygosity using the inbreeding coefficient F. The genomic inflation factor lambda was 1.03 for INVEST individuals genotyped for the HMGA1 IVS5-13insC variant, suggesting minimal population stratification in genotyped individuals. For the IVS5-13insC variant, concordance of duplicates on Taqman® was 99% and concordance between Pyrosequencing and Taqman® was 95%. In Caucasians, the MAF of IVS5-13insC was 0.03 in diabetic cases and 0.04 in controls. The frequency of IVS5-13insC was similar in African American individuals (0.03 in cases and 0.04 in controls) and highest in Hispanics (0.07 in cases and 0.07 in controls). No significant associations were observed between the IVS5-13insC variant and diabetes overall (OR 0.98 95%CI 0.76-1.26, p=0.88). The IVS5-13insC variant was not associated with diabetes in any race/ethnic group (Table 2).
Pairwise LD of HMGA1 SNPs by race/ethnicity is presented in Figure 1. Pairwise LD (r2) of IVS5-13insC and rs9394200, the SNP previously used as a tag SNP for IVS5-13insC , was 0.47 in Caucasians, 0.25 in Hispanics, 0.06 in African Americans. In silico functional analysis of IVS5-13insC revealed no predicted changes to amino acid sequence, conserved transcription factor binding sites, CpG islands, or miRNA regulatory sites. IVS5-13insC was not predicted to cause the creation or disruption of splice sites. However, in silico evaluation of exonic splice enhancer (ESE) sites revealed that the C insertion of IVS5-13insC creates an additional ESE motif for SF2/ASF (IgM-BRCA1).
In the present study, we observed that the HMGA1 IVS5-13insC variant was not associated with diabetes in the INVEST-GENES population overall or in any race/ethnic group. We observed a frequency of IVS5-13insC consistent with previously published studies in Caucasians [5, 7] and observed an increased frequency of this variant in Hispanics. Furthermore, pairwise LD in our population suggests that IVS5-13insC is not effectively tagged by rs9394200 in any race/ethnic group and in silico analysis revealed minimal evidence of putative functional consequences of the IVS5-13insC variant.
Our observations extend existing knowledge of genetic determinants of type 2 diabetes and are consistent with the findings of Marquez et al., which suggest a lack of association between IVS5-13insC and type 2 diabetes . Marquez et al. also found that the IVS5-13insC variant is not functional, which is supported by our lack of observed association in multiple race/ethnic groups and our in silico functional results which predicted no important consequences from this variability. Our findings are inconsistent with two studies that observed a significant association between IVS5-13insC and type 2 diabetes [5, 6]. This discrepancy may be due to differences in race/ethnicity of our population compared with other studies. However, a functional IVS5-13insC variant would be expected to have similar effects across multiple race/ethnicity groups.
Although our association analysis is consistent with the findings of Marquez et al. , our observation of low LD between IVS5-13insC and rs9394200 in all race/ethnic groups suggests that rs9394200 is not an appropriate tag SNP for IVS5-13insC. In Hapmap, the MAF of the T allele for rs9394200 is similar to the MAF for IVS5-13insC in Caucasians (0.03), but is 0.45 in Yorubans, further suggesting that rs9394200 does not adequately tag IVS5-13insC, especially in non-Caucasian populations. Furthermore, rs9394200 is 5000 base pairs downstream from the of HMGA1 3′ end and is represented on arrays utilized by type 2 diabetes GWAS [3, 20]. Therefore, if rs9394200 were a significant contributor to risk for diabetes, it likely would have been identified in the GWAS analyses.
Our observation of a lack of association in Hispanics and African Americans suggests that IVS5-13insC may be associated with diabetes only in Caucasian individuals. While we did not observe even a trend towards association in our Caucasian population, we acknowledge that our power to observe an association in Caucasians is lower than in the previously published studies. Although we had 84 percent power to detect an association in our overall population, we had inadequate power to definitively conclude an association in by race/ethnicity analyses.
Our in silico analysis did not reveal any direct mechanism of the IVS5-13insC variant’s effects on mRNA expression or amino acid sequence, suggesting that the variant is not likely to be functional. Although the results of Chiefari et al. suggested a functional effect of the IVS5-13insC variant on HMGA1 and INSR mRNA and protein expression in monocytes , Marquez et al. found no effect of the variant on mRNA expression in adipose tissue. The differences in observed effect on mRNA expression may potentially be explained by differential effects of the variant on transcription in monocytes and adipose tissue. The apparent functional effect of IVS5-13insC observed by Chiefari et al. may also be confounded by the lack of evaluation of IVS5-13insC variant carriers versus non-carriers among non-diabetic controls  or potential treatment with glucose-lowering medications in diabetic patients from whom monocytes were collected . Finally, IVS5-13insC may have an unknown functional mechanism not identified by in silico tools used in this study. Our study has several limitations worthy of mention. We recognize the potential for false negative results in our analyses, especially in African Americans, considering the low frequency of the variant and the limited power to detect associations within each race/ethic group. Our observations require replication in independent populations of similar race/ethnic makeup. In INVEST, incident and prevalent diabetes diagnosis was based on investigator reports, but the diabetes phenotype is well described in a previous publication , the accuracy of such reporting has been verified by others , and has been used in other trials [23–25]. INVEST investigator-reported diabetes phenotypes have also been used in a large scale gene-centric meta-analysis with HumanCVD BeadChip data , suggesting validity with regard to the diabetes phenotype and genetic association analysis. Although our HMGA1 association analysis may be confounded if control patients eventually develop diabetes after study follow-up, the high mean age and lower mean BMI suggests that our control patients are less likely to develop type 2 diabetes. In addition, although concordance within the Taqman® genotyping platform was high, our concordance between Taqman® and Pyrosequencing was 95%, suggesting some disagreement between platforms. However, Hardy Weinberg tests did not indicate genotype error and genotype error was minimized by utilization of Sanger sequencing for discrepancy confirmation.
Our results suggest that the HMGA1 IVS5-13insC is not associated with type 2 diabetes and may not have an important functional role in diabetes pathogenesis. We also provide frequency and LD data for Hispanic and African American populations, which have higher prevalence of type 2 diabetes. Our results also suggest that rs9394200 is not an effective tag SNP for HMGA1 IVS5-13insC, especially in non-Caucasian populations. Functional studies and replication of these associations are needed to better define the potential role of HMGA1 variants in predicting type 2 diabetes development. Although the current study suggests the lack of a functional role for IVS5-13insC, further study of HMGA1 is warranted to clarify the role of this gene in diabetes pathogenesis. Larger, more adequately powered studies need to be performed to confirm our findings.
Zhang P, Zhang X, Brown J, Vistisen D, Sicree R, Shaw J, Nichols G: Global healthcare expenditure on diabetes for 2010 and 2030. Diabetes Res Clin Pract. 2010, 87: 293-301. 10.1016/j.diabres.2010.01.026.
Herder C, Roden M: Genetics of type 2 diabetes: pathophysiologic and clinical relevance. Eur J Clin Invest. 2011, 41: 679-692. 10.1111/j.1365-2362.2010.02454.x.
Voight BF, Scott LJ, Steinthorsdottir V, Morris AP, Dina C, Welch RP, Zeggini E, Huth C, Aulchenko YS, Thorleifsson G: Twelve type 2 diabetes susceptibility loci identified through large-scale association analysis. Nat Genet. 2010, 42: 579-589. 10.1038/ng.609.
Talmud PJ, Hingorani AD, Cooper JA, Marmot MG, Brunner EJ, Kumari M, Kivimaki M, Humphries SE: Utility of genetic and non-genetic risk factors in prediction of type 2 diabetes: Whitehall II prospective cohort study. BMJ. 2010, 340: b4838-10.1136/bmj.b4838.
Chiefari E, Tanyolac S, Paonessa F, Pullinger CR, Capula C, Iiritano S, Mazza T, Forlin M, Fusco A, Durlach V: Functional variants of the HMGA1 gene and type 2 diabetes mellitus. JAMA. 2011, 305: 903-912. 10.1001/jama.2011.207.
Liu L, Ding H, Wang HR, Xu YJ, Cui GL, Wang PH, Yuan G, Yu XF, Wang DW: Polymorphism of HMGA1 is associated with increased risk of type 2 diabetes among Chinese individuals. Diabetologia. 2012, 55: 1685-1688. 10.1007/s00125-012-2518-0.
Marquez M, Huyvaert M, Perry JR, Pearson RD, Falchi M, Morris AP, Vivequin S, Lobbens S, Yengo L, Gaget S: Low-frequency variants in HMGA1 are not associated with type 2 diabetes risk. Diabetes. 2012, 61: 524-530. 10.2337/db11-0728.
Roger VL, Go AS, Lloyd-Jones DM, Benjamin EJ, Berry JD, Borden WB, Bravata DM, Dai S, Ford ES, Fox CS: Executive Summary: Heart Disease and Stroke Statistics–2012 Update: A Report From the American Heart Association. Circulation. 2012, 125: 188-197.
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. 10.1016/S0735-1097(98)00423-9.
Pepine CJ, Handberg EM, Cooper-DeHoff RM, Marks RG, Kowey P, Messerli FH, Mancia G, Cangiano JL, Garcia-Barreto D, Keltai M: 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. 10.1001/jama.290.21.2805.
Cooper-Dehoff R, Cohen JD, Bakris GL, Messerli FH, Erdine S, Hewkin AC, Kupfer S, Pepine CJ: Predictors of development of diabetes mellitus in patients with coronary artery disease taking antihypertensive medications (findings from the INternational VErapamil SR-Trandolapril STudy [INVEST]). Am J Cardiol. 2006, 98: 890-894. 10.1016/j.amjcard.2006.04.030.
Cooper-DeHoff RM, Handberg EM, Mancia G, Zhou Q, Champion A, Legler UF, Pepine CJ: INVEST revisited: review of findings from the International Verapamil SR-Trandolapril Study. Expert Rev Cardiovasc Ther. 2009, 7: 1329-1340. 10.1586/erc.09.102.
Keating BJ, Tischfield S, Murray SS, Bhangale T, Price TS, Glessner JT, Galver L, Barrett JC, Grant SF, Farlow DN: Concept, design and implementation of a cardiovascular gene-centric 50 k SNP array for large-scale genomic association studies. PLoS One. 2008, 3: e3583-10.1371/journal.pone.0003583.
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. 10.1086/519795.
Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.
Cooper-DeHoff RM, Aranda JM, Gaxiola E, Cangiano JL, Garcia-Barreto D, Conti CR, Hewkin A, Pepine CJ: Blood pressure control and cardiovascular outcomes in high-risk Hispanic patients–findings from the International Verapamil SR/Trandolapril Study (INVEST). Am Heart J. 2006, 151: 1072-1079. 10.1016/j.ahj.2005.05.024.
Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21: 263-265. 10.1093/bioinformatics/bth457.
Chelala C, Khan A, Lemoine NR: SNPnexus: a web database for functional annotation of newly discovered and public domain single nucleotide polymorphisms. Bioinformatics. 2009, 25: 655-661. 10.1093/bioinformatics/btn653.
Cartegni L, Wang J, Zhu Z, Zhang MQ, Krainer AR: ESEfinder: A web resource to identify exonic splicing enhancers. Nucleic Acids Res. 2003, 31: 3568-3571. 10.1093/nar/gkg616.
Burton PR, Clayton DG, Cardon LR, Craddock N, Deloukas P, Duncanson A, Kwiatkowski DP, McCarthy MI, Ouwehand WH, Samani NJ: Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 2007, 447: 661-678. 10.1038/nature05911.
Froguel P, Marquez M, Cauchi S: Response to comment on: Marquez et al. Low-frequency variants in HMGA1 are not associated with type 2 diabetes risk. Diabetes 2012;61:524–530. Diabetes. 2012, 61: e15-10.2337/db12-0800.
Goldman N, Lin IF, Weinstein M, Lin YH: Evaluating the quality of self-reports of hypertension and diabetes. J Clin Epidemiol. 2003, 56: 148-154. 10.1016/S0895-4356(02)00580-2.
Aguilar D, Solomon SD, Kober L, Rouleau JL, Skali H, McMurray JJ, Francis GS, Henis M, O’Connor CM, Diaz R: Newly diagnosed and previously known diabetes mellitus and 1-year outcomes of acute myocardial infarction: the VALsartan In Acute myocardial iNfarcTion (VALIANT) trial. Circulation. 2004, 110: 1572-1578. 10.1161/01.CIR.0000142047.28024.F2.
Bosch J, Lonn E, Pogue J, Arnold JM, Dagenais GR, Yusuf S: Long-term effects of ramipril on cardiovascular events and on diabetes: results of the HOPE study extension. Circulation. 2005, 112: 1339-1346.
Yusuf S, Gerstein H, Hoogwerf B, Pogue J, Bosch J, Wolffenbuttel BH, Zinman B: Ramipril and the development of diabetes. JAMA. 2001, 286: 1882-1885. 10.1001/jama.286.15.1882.
Saxena R, Elbers CC, Guo Y, Peter I, Gaunt TR, Mega JL, Lanktree MB, Tare A, Castillo BA, Li YR: Large-scale gene-centric meta-analysis across 39 studies identifies type 2 diabetes loci. Am J Hum Genet. 2012, 90: 410-425. 10.1016/j.ajhg.2011.12.022.
This study was funded by NIH grant U01 GM074492, R01 HL74730, HL086558 (RCD), and TL1RR029888 (JHK), and grants from Abbott Laboratories and the University of Florida Opportunity Fund. We thank Ben Burkley, Lynda Stauffer, and Cheryl Galloway for processing and genotyping samples.
The authors declare that they have no competing interests.
JHK drafted the manuscript and performed statistical analysis. TYL, JRC, MR, OEC, and SWC performed genotyping and helped draft the manuscript. CWM and YG performed HumanCVD BeadChip quality control procedures, contributed to statistical analysis and helped draft the manuscript. CJP, JAJ, and RCD conceived of the study, participated in its design and coordination, and helped draft the manuscript. All authors read and approved the final manuscript.