Increased prevalence of the CVD-associated ANRIL allele in the Roma/Gypsy population in comparison with the majority Czech population

Abstract

Cardiovascular disease (CVD) is a major cause of death around the world, with highest prevalence reported in minority Roma/Gypsy populations living in developed countries. Whether these differences are caused by unhealthy lifestyles or genetic factors remain unknown. The aim of our study was to examine the genotype frequencies of the rs10757274 polymorphism in the 9p.21 locus within ANRIL (antisense non-coding RNA in the INK4 locus), a long non-coding RNA located in the vicinity of the CDKN2A/2B inhibitors loci. ANRIL is understood to be the strongest genetic determinant of CVD in Caucasians. Using PCR-RFLP, we analysed the ANRIL rs10757274 polymorphism in 298 non-Roma (50% male) and 302 Roma/Gypsy (50% male) adult (39.5 ± 15.1 years and 39.2 ± 12.8 years, respectively) subjects. We found that frequencies of the ANRIL GG, GA and AA genotypes were 20.1%, 52.4% and 27.5% in the majority population and 32.9%, 47.9% and 19.2% in Roma/Gypsy subjects, respectively. The distribution of genotypes was deemed significantly different at P < 0.001. Within the Roma/Gypsy population, we detected increased prevalence of the CVD-associated GG genotype. Increased prevalence of CVD among Roma/Gypsies subjects may be significantly linked to genetic background.

Keywords:
ANRIL; polymorphism; cardiovascular disease; Roma/Gypsy; ethnic differences

Cardiovascular diseases and myocardial infarction (CVD, MI) are the major causes of death in industrial countries. It is widely acknowledged that, next to environmental and lifestyle risk factors (Bhatnagar, 2017), genetic predisposition plays an important role in CVD determination (Dainis and Ashley, 2018).

Although previous publications on this topic are sparse and typically characterised by low numbers of examined subjects and non-representative selections, they report that prevalence of CVD risk factors in Roma/Gypsy communities is higher than in majority populations (Dobranici et al., 2012). A Czech national survey (Antosova et al., 2014) reported that prevalence of heart disease in the Roma/Gypsy population living in the Czech Republic was more than double that of the majority population. The question remains whether increased prevalence of CVD is caused not just by unhealthy lifestyles but also by genetic factors, as genetic differences between Roma/Gypsy and non-Roma populations have been widely documented (Mydlárová Blaščáková et al., 2017; Hubacek et al., 2017, 2020; Nagy et al., 2017; Dlouhá et al., 2020).

A number of genome-wide association studies (GWAS) have detected a comprehensive list of common polymorphisms associated with increased risk of MI (for a summary, see Erdmann et al., 2018). Among the strongest are the single-nucleotide polymorphisms (rs10757274, rs10965215, rs1333040 and rs4977574 are among the most commonly studied) within ANRIL (antisense non-coding RNA in the INK4 locus; OMIM acc. N. 613149). A long non-coding regulatory RNA located in the 9p.21 locus (Helgadottir et al., 2007), ANRIL is clustered within the cyclin-dependent kinase inhibitors CDKN2A and CDKN2B. ANRIL is understood to regulate - through a complex mechanism - a number of processes implicated in endothelial injury and subsequent development of atherosclerosis. That ANRIL is expressed in almost all human tissue highlights the extraordinary importance of this regulatory RNA.

At-risk ANRIL alleles are associated with an approximate 30-35% increased risk of MI per each allele (Palomaki et al., 2010). Interestingly, despite the homogeneous and strong association between ANRIL SNPs and CVD risk among different ethnicities (Hu et al., 2019), identical variants are not typically associated with traditional CVD risk factors such as smoking, dyslipidaemia, obesity and hypertension (with the exception of diabetes). Nevertheless, even after adjusting for diabetes, the highly significant association between ANRIL polymorphisms and CVD indicates the independent influence of this locus on these two major non-communicable diseases.

Two groups of adult populations (the non-Roma Czech majority and the Roma/Gypsy ethnic minority) inhabiting a Czech region South Bohemia (Table 1; Hubacek et al., 2020) have been included in the study. Snowball sampling (Hughes et al., 1995) has been used to the Czech Roma/Gypsy subpopulation (N = 302) selection and quota sampling (Walter, 1989) for recruiting Czech Caucasians/Slavs (N = 298). All subjects were older than 18 years at the time of data and samples collection and theirs ethnicity was based on self-reported information (Adámková et al., 2015, Šedová et al., 2015, Hubáček et al., 2018).

Ethics committee at South-Bohemia University approved the study protocol and all participants signed the informed consent with participation in the study.

Table 1 -
General characteristics of examined subjects (data identifying cohorts have been published in Hubacek et al., 2020).

	Roma/Gypsies	Majority Slavs	P
N	302	298
Age (years)	39.2 ± 12.8	39.5 ± 15.1	n.s.
% of females	50	50	n.s.
BMI (kg/m²)	29.9 ± 5.6	25.0 ± 6.0	0.01
SBP (mm Hg)	124.3 ± 20.9	124.9 ± 14.4	n.s.
DBP (mm Hg)	76.6 ± 12.7	77.2 ±10.6	n.s.
Total cholesterol (mmol/L)	5.1 ± 1.4	5.1 ± 1.1	n.s.
Glycaemia (>6 mmol/L) (%)	58.6	23.3	0.00001
Total body fat (>25%) (%)	69.4	61.7	0.05

DNA were isolated (“Xtreme DNA Isolation Kit”), from buccal cells collected through “DNA buccal swabs” (both Isohelix, Cell Projects Ltd, UK) according to manufacturer conditions.

The ANRIL rs10757274 genotype was analysed using the polymerase chain reaction-restriction fragment length polymorphism method, as described in detail elsewhere (Hubacek et al., 2016). Restriction fragments were separated using a 10% polyacrylamide gel. Fermentas International Inc. (Burlington, Ontario, Canada) provided all PCR chemicals, with PCRs performed on the MJ Research DYAD Disciple PCR device.

Deviance of genotype frequencies among the groups examined was analysed using the Hardy-Weinberg equilibrium principle (Court, 2005-2008). Differences in allele and genotype frequencies were compared using an online chi-square test (Stangroom). Comparisons were performed for dominant, co-dominant and recessive models.

The call rates achieved for the ANRIL rs10757274 SNP were 96.6% in the Czech majority population and 96.7% in the minority population. No significant gender differences in genotype frequencies were observed in either ethnic group (not presented in detail). Distributions of genotypes were in agreement with the Hardy-Weinberg equilibrium (P = 0.36 for the majority population and P = 0.70 for the minority population).

Within the minority Roma/Gypsy population, there were significantly more carriers of the CVD-associated GG genotype (P < 0.001; for further details and comparisons, see Table 2). Consequently, prevalence of the G allele was also significantly higher (56.8% vs. 46.4%; P < 0.0005).

Table 2 -
Distribution of ANRIL rs10757274 genotypes in the majority Czech non-Roma and minority Roma/Gypsy populations.

ANRIL	Czech majority		Roma/Gypsy minority		P
rs10757274	N	%	N	%
GG	58	20.1	96	32.9	*0.0005
GA	151	52.4	140	47.9	#0.001
AA	79	27.5	56	19.2	§0.02
G	267	46.4	332	56.8	0.0005
A	309	53.6	252	43.2

• P-values are given for *GG vs. +A, ^#GG vs. GA vs. AA, and ^§+G vs. AA comparisons.

Despite the relatively low numbers of subjects enrolled in this study, the genotype frequencies observed for the majority population are in agreement with the frequency (22% of ANRIL rs10757274 GG homozygosity) reported in our previous study (Hubacek et al., 2016), as well as in neighbouring German population (Scheffold et al., 2011). In contrast, detected frequencies of the G allele in the minority Roma/Gypsy population (57%); are among the highest in the world and even higher than in the Indian population (Kumar et al., 2011), which is understood to be the origin of the settled Roma/Gypsy populations in Europe (Mendizabal et al., 2011).

Although examining a sole SNP is evidently insufficient at drawing major conclusions, we speculate that the increased prevalence of at-risk ANRIL genotypes in the Roma/Gypsy population could be an important modifier associated with higher prevalence of CVD in this ethnic group. To extend our findings, a further case-control study involving CVD-focused tagging of ANRIL SNPs within this ethnic group is required.

We conclude that CVD-associated ANRIL genotypes are more common in the Roma/Gypsy population compared to the majority Caucasian population.

Acknowledgments

The study was supported by Ministry of Health of the Czech Republic, grant nr. NV18-01-00046, all rights reserved; projects number LD14114 implemented under the financial support of the Ministry of Education, Youth and Sports within COST (Cooperation On Scientific and Technical Research), entitled “Obesity and overweight in the Romany minority in the Region of South Bohemia”, and by the project supported by Ministry of Health Czech Republic - conceptual development of research organisation (“Institute for Clinical and Experimental Medicine - IKEM, IN 00023001”).

References

• Adámková V, Hubáček JA, Nováková D, Dolák F, Adámek V, Lánská V, Tóthová V and Šedová L (2015) Genetic and biochemical characteristics in the Roma minority in the South Bohemia Region. Neuro Endocrinol Lett 36:29-34.
• Bhatnagar A (2017) Environmental determinants of cardiovascular disease. Circ Res 121:162-180.
• Dainis AM and Ashley EA (2018) Cardiovascular precision medicine in the genomics era. JACC Basic Transl Sci 3:313-326.
• Dlouhá L, Adámková V, Šedová L, Olišarová V, Hubáček JA and Tóthová V (2020) Five genetic polymorphisms of cytochrome P450 enzymes in the Czech non-Roma and Czech Roma population samples. Drug Metab Pers Ther. DOI: 10.1515/dmdi-2020-0103.
  » https://doi.org/10.1515/dmdi-2020-0103
• Dobranici M, Buzea A and Popescu R (2012) The cardiovascular risk factors of the Roma (Gypsies) people in Central-Eastern Europe: a review of the published literature. J Med Life 5:382-389.
• Erdmann J, Kessler T, Munoz Venegas L and Schunkert H (2018) A decade of genome-wide association studies for coronary artery disease: the challenges ahead. Cardiovasc Res 114:1241-1257.
• Helgadottir A, Thorleifsson G, Manolescu A, Gretarsdottir S, Blondal T, Jonasdottir A, Jonasdottir A, Sigurdsson A, Baker A, Palsson A, et al (2007) A common variant on chromosome 9p21 affects the risk of myocardial infarction. Science 316:1491-1493.
• Hu L, Su G and Wang X (2019) The roles of ANRIL polymorphisms in coronary artery disease: a meta-analysis. Biosci Rep 39:BSR20181559.
• Hubacek JA, Vrablik M, Dlouha D, Stanek V, Gebauerova M, Adamkova V, Ceska R, Dostálová G, Linhart A, Vitek L, et al (2016) Gene variants at FTO, 9p21, and 2q36.3 are age-independently associated with myocardial infarction in Czech men. Clin Chim Acta 454:119-123.
• Hubácek JA, Adámková V, Šedová L, Olišarová V, Adámek V and Tóthová V (2017) Frequency of adult type-associated lactase persistence LCT-13910C/T genotypes in the Czech/Slav and Czech Roma/Gypsy populations. Genet Mol Biol 40:450-452.
• Hubáček JA, Šedová L, Olišarová V, Adámková V, Adámek V and Tóthová V (2018) Distribution of ADH1B genotypes predisposed to enhanced alcohol consumption in the Czech Roma/Gypsy population. Cent Eur J Public Health 26:284-288.
• Hubáček JA, Šedová L, Olišarová V, Adámková V and Tóthová V (2020) Different prevalence of T2DM risk alleles in Roma population in comparison with the majority Czech population. Mol Genet Genomic Med 8:e1361.
• Hughes AO, Fenton S, Hine CE, Pilgrim S and Tibbs N (1995) Strategies for sampling black and ethnic minority populations. J Publ Health Med. 17:187-192.
• Kumar J, Yumnam S, Basu T, Ghosh A, Garg G, Karthikeyan G and Sengupta S (2011) Association of polymorphisms in 9p21 region with CAD in North Indian population: replication of SNPs identified through GWAS. Clin Genet 79:588-593.
• Mendizabal I, Valente C, Gusmão A, Alves C, Gomes V, Goios A, Parson W, Calafell F, Alvarez L, Amorim A et al (2011) Reconstructing the Indian origin and dispersal of the European Roma: a maternal genetic perspective. PLoS One 6:e15988.
• Mydlárová Blaščáková M, Blaščáková Ľ, Poráčová J, Mydlár J, Vašková J, Bernasovská J, Boroňová I, Petrejčíková E and Bernasovský I (2017) Relationship between A163G osteoprotegerin gene polymorphism and other osteoporosis parameters in Roma and non-Roma postmenopausal women in eastern Slovakia. J Clin Lab Anal 31:e22093.
• Nagy K, Fiatal S, Sándor J and Ádány R (2017) Distinct penetrance of obesity-associated susceptibility alleles in the hungarian general and Roma populations. Obes Facts 10:444-457.
• Palomaki GE, Melillo S and Bradley LA (2010) Association between 9p21 genomic markers and heart disease: a meta-analysis. JAMA 303:648-656.
• Scheffold T, Kullmann S, Huge A, Binner P, Ochs HR, Schöls W, Thale J, Motz W, Hegge FJ, Stellbrink C, et al (2011) Six sequence variants on chromosome 9p21.3 are associated with a positive family history of myocardial infarction: a multicenter registry. BMC Cardiovasc Disord 11:9.
• Šedová L, Tóthová V, Olišarová V, Adámkova V, Bártlová S, Dolák F, Kajanová A, Mauritzová I, Nováková D and Prokešová R (2015) Evaluation of selected indicators of overweight and obesity of Roma minority in the region of South Bohemia. Neuro Endocrinol Lett 36:35-42.
• Walter SD (1989) The feasibility of matching and quota sampling in epidemiologic studies. Am J Epidemiol 130:379-389.

Internet resources

• Antosova D, Benes C, Csemy L, Castkova J, Fabianova K, Filipova V, Grolmusova L, Jakubu V, Janatova H, Jindrak V, et al (2014). Table. 4.13.1 Chronická neinfekční onemocnění u romské a většinové populace v ČR v roce 2001, 1 Chronická neinfekční onemocnění u romské a většinové populace v ČR v roce 2001, https://www.mzcr.cz/wp-content/uploads/wepub/9420/20954/Zpráva%20o%20zdraví%20obyvatel%20ČR%202014.pdf (accessed 5 February 2021)
  » https://www.mzcr.cz/wp-content/uploads/wepub/9420/20954/Zpráva%20o%20zdraví%20obyvatel%20ČR%202014.pdf
• Stangroom J, Stangroom J, http://www.socscistatistics.com (accessed 5 February 2021).
  » http://www.socscistatistics.com
• Court (2005-2008), Court (2005-2008), https://www.researchgate.net/profile/Khalid-Hassan/post/Does-anyone-know-how-to-calculate-Hardy-Weinberg-Equilibrium-for-14-different-alleles/attachment/59d64ad979197b80779a50fd/AS%3A477303651213312%401490809466020/download/Court+lab+-+HW+calculator-+important.xls (accessed 5 February 2021).
  » http://www.tufts.edu/~mcourt01/documents/court%20lab%20-%20hw%20calculator.xls

• Associate Editor:

  Mara Hutz

Publication Dates

• Publication in this collection
  30 Apr 2021
• Date of issue
  2021

History

• Received
  06 Nov 2020
• Accepted
  21 Feb 2021