Oncotarget

Clinical Research Papers:

Interleukin-27 polymorphisms are associated with premature coronary artery disease and metabolic parameters in the Mexican population: the genetics of atherosclerotic disease (GEA) Mexican study

PDF |  HTML  |  How to cite

Oncotarget. 2017; 8:64459-64470. https://doi.org/10.18632/oncotarget.16223

Metrics: PDF 1612 views  |   HTML 2969 views  |   ?  

Rosalinda Posadas-Sánchez _, Nonanzit Pérez-Hernández, José Manuel Rodríguez-Pérez, Ramón M. Coral-Vázquez, Bladimir Roque-Ramírez, Luis Llorente, Guadalupe Lima, Carmina Flores-Dominguez, Teresa Villarreal-Molina, Carlos Posadas-Romero and Gilberto Vargas-Alarcón

Abstract

Rosalinda Posadas-Sánchez1, Nonanzit Pérez-Hernández2, José Manuel Rodríguez-Pérez2, Ramón M. Coral-Vázquez3, Bladimir Roque-Ramírez4, Luis Llorente5, Guadalupe Lima5, Carmina Flores-Dominguez7, Teresa Villarreal-Molina6, Carlos Posadas-Romero1 and Gilberto Vargas-Alarcón2

1 Departamento de Endocrinología, Instituto Nacional de Cardiología Ignacio Chávez, Mexico D.F., México

2 Departamento de Biología Molecular, Instituto Nacional de Cardiología Ignacio Chávez, Mexico D.F., México

3 Sección de Estudios de Posgrado e Investigación, Escuela Superior de Medicina, Instituto Politécnico Nacional, Mexico D.F., México

4 Departamento de Farmacobiologia CINVESTAV-Sede Sur, Mexico D.F., México

5 Departamento de Inmunología y Reumatología, Instituto Nacional de Ciencias Médicas y Nutrición Salvador Zubirán, Mexico D.F., México

6 Laboratorio de Genómica Cardiovascular, Instituto Nacional de Medicina Genómica, Mexico D.F., México

7 CICSA, Universidad Anahuac, Estado de México, México

Correspondence to:

Gilberto Vargas-Alarcón, email:

Keywords: association studies, coronary artery disease, inflammation, interleukin 27, polymorphism genetics

Received: December 19, 2016 Accepted: March 03, 2017 Published: March 15, 2017

Abstract

Several studies suggest an important role of Interleukin-27 in the development of atherosclerosis. The aim of this study was to establish whether the IL-27p28 gene polymorphisms are associated with premature coronary artery disease and/or other cardiovascular risk factors. Four IL-27p28 gene polymorphisms were selected and genotyped in 1162 premature coronary artery disease cases and 1107 controls. rs26528 T and rs40837 A alleles were significantly associated with a lower risk of premature coronary artery disease under different inheritance models (Pdominant = 0.046; Pover-dominant = 0.002; Pco-dominant1 = 0.007 for rs26528T; Pover-dominant = 0.008 and Pco-dominant1 = 0.031 for rs40837). The rs40837 A allele was also associated with a lower risk of insulin resistance, in cases (Pover-dominant = 0.037) and controls (Padditive = 0.008; Pdominant = 0.047; Precessive = 0.014; Pco-dominant2 = 0.006), while the rs26528 T allele was associated with a lower risk of insulin resistance only in the control group (Precessive = 0.016; Pco-dominant2 = 0.021). Interleukin-27 plasma levels were measured in 450 controls and 450 cases, and were significantly higher in cases compared to controls (P = 0.004). However, Interleukin-27 plasma levels were not associated with IL-27p28 polymorphisms. Luciferase assays showed that co-transfection of the rs40837 A allele and miR-379-5p significantly decreased luciferase gene expression. Our study shows for the first time, that IL-27p28 gene polymorphisms are associated with premature coronary artery disease and with some metabolic parameters. The rs40837 A allele in presence of miR-379-5p significantly decreased luciferase gene expression.


INTRODUCTION

Cardiovascular disease (CVD) is the main cause of morbidity in developed and emerging countries. Coronary arterial disease (CAD) is the most common CVDs mainly caused by atherosclerosis, a multifactorial disease involving both genetic and environmental factors [1]. Several studies suggest atherosclerosis could be considered an inflammatory disease [2-4]. Macrophage and T cell infiltrates are known to play an important role in atherosclerotic lesions in humans and in animal models [5,6]. Macrophages [7,8] and Th1 cells [9-11] secreted cytokines and chemokines that amplify local immune responses, while Th2 cells are known to have a protective effect on the development of atherosclerosis [12-13].

Interleukin 27 (IL-27), a member of the interleukin 12 family, is a heterodimeric cytokine, conformed by an α p28 and a β EBI3 subunit [14]. IL-27 is an immune/inflammatory response regulator [15] by promoting early Th1 differentiation [16], suppressing Th2 and Th17 differentiation [17-20], and inducing the production of anti-inflammatory cytokines such as IL-10 by activated T cell [21]. IL-27 is expressed in atherosclerotic plaques [22], and its role in atherosclerosis has been studied in cultured cells, animal models and coronary patients, with inconsistent findings. IL-27-deficient (Ldlr−/−Ebi3−/−) and IL-27 receptor-deficient (Ldlr/−WSX-1−/−) Ldlr−/− knockout mice were more susceptible to develop atherosclerosis, and IL-27 administration suppressed macrophage activation and atherosclerosis development [23]. However, coronary patients showed higher IL-27 levels as compared to controls, and IL-27 levels showed a significant correlation with stenosis severity [24]. Moreover, dendritic cells incubated with oxidized LDL (low density lipoprotein) produced IL-27, suggesting these modified lipoproteins could play an important role in dendritic cell activation and IL-27 production [24]. Altogether, these studies suggest IL-27 could play a crucial role in the immunity and inflammation regulatory net in atherosclerosis. The human IL-27p28 gene encodes the IL-27 alpha subunit, located in the 16p11 locus, spans 5 exons and is highly polymorphic [25]. Considering the important role of the IL-27 in the developing of atherosclerosis, the objective of this study was to evaluate whether IL-27p28 polymorphisms are associated with premature CAD (pCAD) and/or cardiovascular risk factors, as well as to evaluate whether the associated polymorphisms have a functional effect.

RESULTS

Clinical characteristics and metabolic parameters

A total of 1107 controls with no tomographic evidence [coronary artery calcification (CAC) score = 0] of subclinical atherosclerosis (SA) and 1162 pCAD cases with complete clinical, demographic, anthropometric and biochemical information belonging to the Genetics of Atherosclerotic Disease (GEA) Mexican Study were selected for the analyses. Clinical characteristics of the pCAD cases and control subjects are shown in Table 1. Age, male percentage, body mass index (BMI), waist circumference, systolic and diastolic blood pressure, visceral abdominal fat (VAF), alanine aminotransferase (ALT), aspartate aminotransferase (AST), hypoalphalipoproteinemia, hypertriglyceridemia, general obesity, abdominal obesity, type 2 diabetes mellitus (T2DM), insulin resistance, metabolic syndrome, hypertension, high VAF, hyperuricemia and hypoadiponectinemia were significantly higher in pCAD cases than in controls (Table 1). On the other hand, hypercholesterolemia [total cholesterol (TC)>200 mg/dL or low density lipoprotein-cholesterol (LDL-C) ≥130 mg/dL], high non-HDL (high density lipoprotein) cholesterol, inflammation [defined as high sensitivity C reactive protein (hsCRP) levels ≥3mg/L] and current smoking habit were significantly more frequent in controls than in pCAD cases most likely due to the effect of statin treatment (Table 1). All differences were statistically significant.

rs40837

Figure 1: rs40837 affects luciferase expression levels in the presence of miR-379-5p. A. Schematic diagram of the IL27 gene indicating the location of rs40837 allelic variants (A/G) in the 3’-UTR region and its sequence complementarity with the miR-379-5p. B. Sequences with each allele were cloned in pMIR-REPORT vector to evaluate the expression of the luciferase reporter gene as described in Materials and Methods. C. The presence of allele A reduced luciferase expression in the cells co-transfected with miR-379-5p as compared to cells co-transfected with miR control (P = 0.0022), and cells with G allele co-transfected with miR-379-5p (P = 0.0087). n = 5 independent duplicate experiments; medians were compared using Mann Whitney’s U test.

Table 1: Clinical characteristics and cardiovascular risk factors in the studied groups.

Controls (n= 1107)

pCAD (n = 1162)

aP

Age (years)

51 ± 9

54 ± 8

<0.001

Sex (% male)

41.2

81.1

<0.001

Body mass index (kg/m2)

27.8 [25.4-30.8]

28.3 [25.9-31.1]

0.003

Waist circumference (cm)

94 ± 11

98 ± 10

<0.001

Systolic blood pressure (mmHg)

102 [104-122]

116 [106-127]

<0.001

Diastolic blood pressure (mmHg)

70 [65-76]

71 [66-78]

0.013

Visceral abdominal fat (cm2)

139 [10.4-180]

168 [129-215]

<0.001

ALT (IU/L)

24 [18-34]

26 [19-36]

0.031

AST (IU/L)

25 [21-30]

26 [22-31]

0.001

Total cholesterol > 200mg/dl (%)

36.5

20.3

<0.001

LDL-C ≥ 130 mg/dL (%)

29.6

16.1

<0.001

Hypoalphalipoproteinemia (%)

52.0

67.2

<0.001

Hypertriglyceridemia (%)

47.3

56.2

<0.001

Non HDL-C>160 mg/dL (%)

27.9

19.5

<0.001

Obesity (%)

30.1

35.0

0.003

Abdominal obesity (%)

81.0

83.6

0.060

Type 2 Diabetes mellitus (%)

10.1

35.4

<0.001

Insulin resistance (%)

54.2

77.0

<0.001

Metabolic syndrome (%)

40.7

71.9

<0.001

Hypertension (%)

6.7

68.1

<0.001

Increased VAF (%)

54.5

64.6

<0.001

Current smoking habit (%)

22.8

11.6

<0.001

Hypoadiponectinemia (%)

41.5

56.5

<0.001

hsCRP ≥ 3 mg/L (%)

26.6

21.3

0.002

Hyperuricemia (%)

19.9

35.9

<0.001

Abbreviations: pCAD, premature coronary artery disease; ALT, alanine aminotransferase; AST, aspartate aminotransferase; LDL-C, low density lipoprotein cholesterol; Non HDL-C, non high density lipoprotein cholesterol; VAF, visceral abdominal fat; hsCRP, high sensitivity C reactive protein. Data are shown as average ± standard deviation, median [interquartile range] or percentage. aStudent’s t test, Mann Whitney’s U test or Chi square test.

Association of IL-27 polymorphisms with pCAD

Genotype distributions of all polymorphisms in pCAD cases and controls are described in Table 2. The polymorphisms were in Hardy-Weinberg equilibrium. While genotype distributions of rs17855750 and rs181206 were similar in cases and controls, rs26528 T and rs40837 A alleles were significantly associated with a lower risk of pCAD. The rs26528 T allele showed significant associations with lower risk of pCAD under dominant (OR = 0.794, P = 0.046), over-dominant (OR = 0.701, P = 0.002) and co-dominant 1 (OR = 0.718, P = 0.007) models, while the rs40837 A allele was significantly associated with a lower risk of pCAD under the over-dominant and co-dominant 1 models (OR = 0.740, P = 0.008 and OR = 0.768, P = 0.031, respectively). All associations were adjusted for age, gender, BMI, smoking habit, total abdominal fat (TAF), homeostasis model assessment of insulin resistance (HOMA-IR), AST, adiponectin and uric acid levels (Table 2).

Table 2: Association of IL-27p28 polymorphisms with premature coronary artery disease.

SNP

Genotype frequency

MAF

Model

OR [95% CI]

P

rs26528

CC

CT

TT

Additive

0.958 [0.817-1.122]

0.594

Control (n = 1107)

0.378

0.481

0.142

0.382

Dominant

0.794 [0.634-0.996]

0.046

Patients (n = 1162)

0.410

0.441

0.149

0.369

Recessive

1.302 [0.958-1.771]

0.092

Over-dominant

0.701 [0.562-0.875]

0.002

Co-dominant 1

0.718 [0.564-0.913]

0.007

Co-dominant 2

1.086 [0.778-1.517]

0.627

rs17855750

AA

AC

CC

Additive

0.967 [0.777-1.205]

0.767

Control (n = 1107)

0.729

0.248

0.023

0.147

Dominant

0.932 [0.724-1.199]

0.582

Patients (n = 1162)

0.731

0.242

0.024

0.146

Recessive

1.240 [0.610-2.522]

0.553

Over-dominant

0.900 [0.694-1.169]

0.431

Co-dominant 1

0.906 [0.697-1.177]

0.460

Co-dominant 2

1.211 [0.594-2.470]

0.599

rs181206

AA

AG

GG

Additive

1.003 [0.855-1.178]

0.967

Control (n = 1107)

0.451

0.426

0.123

0.336

Dominant

0.955 [0.768-1.187]

0.677

Patients (n = 1162)

0.442

0.441

0.117

0.337

Recessive

1.135 [0.810-1.590]

0.463

Over-dominant

0.923 [0.740-1.151]

0.476

Co-dominant 1

0.936 [0.742-1.182]

0.581

Co-dominant 2

1.073 [0.749-1.536]

0.701

rs40837

GG

GA

AA

Additive

0.992 [0.847-1.162]

0.923

Control (n = 1107)

0.378

0.480

0.143

0.382

Dominant

0.847 [0.676-1.062]

0.150

Patients (n = 1162)

0.405

0.445

0.150

0.372

Recessive

1.316 [0.970-1.786]

0.078

Over-dominant

0.740 [0.593-0.923]

0.008

Co-dominant 1

0.768 [0.604-0.977]

0.031

Co-dominant 2

1.140 [0.818-1.589]

0.440

Abbreviations: SNP, single nucleotide polymorphism; MAF, minor allele frequency; OR, odds ratio; 95% CI, 95% confidence interval. All models were adjusted by age, gender, body mass index, smoking habit, total abdominal fat, HOMA-IR, aspartate aminotransferase, adiponectin and uric acid.

Co-dominant 1: comparison of heterozygote vs major allele homozygote genotype frequencies.

Co-dominant 2: comparison of heterozygote vs minor allele homozygote genotype frequencies.

Associations with metabolic parameters

Tables 3 and 4 describe associations of IL-27p28 polymorphisms with metabolic parameters in controls and pCAD cases, respectively. In the control group, rs26528 T and rs40837 A were significantly associated with a lower risk of insulin resistance (OR = 0.623, Precessive = 0.016; OR = 0.610, Pco-dominant2 = 0.021 for rs26528 T; OR = 0.646, Padditive = 0.008; OR = 0.616, Pdominant = 0.047; OR = 0.488, Precessive =0.014 and OR = 0.404, Pco-dominant2 = 0.006 for rs40837 A) and with a lower risk of increased AST activity (OR = 0.777, Padditive = 0.016; OR = 0.654, Precessive = 0.046 and OR = 0.583, Pco-dominant2 = 0.018 for rs26528 T; OR = 0.771, Padditive = 0.013 and OR = 0.603, Pco-dominant2 = 0.025 for rs40837 A); while rs181206 G was associated with a lower risk of hyperuricemia (OR = 0.788, Padditive = 0.044) (Table 3). In pCAD cases, rs26528 T and rs40837 A were significantly associated with a lower risk of elevated AST activity (OR = 0.802, Padditive = 0.014; OR = 0.689, Pdominant = 0.003; OR = 0.737, Pover-dominant= 0.016 and OR = 0.678, Pco-dominant1 = 0.004 for rs26528 T; OR = 0.799, Padditive = 0.013; OR = 0.720, Pover-dominant = 0.009 and OR = 0.662, Pco-dominant1 = 0.002 for rs40837 A). Moreover, rs181206 G was associated with a lower risk of hyperuricemia (OR = 0.546, Precessive = 0.005 and OR = 0.566, Pco-dominant2 = 0.011) and rs40837 A was associated with a lower risk of insulin resistance (OR = 0.702, Pover-dominant = 0.037) (Table 4).

Table 3: IL27p28 gene polymorphism association with metabolic parameters in control subjects.

SNP

Genotype frequency

MAF

Model

OR [95% CI]

P

rs26528

CC

CT

TT

Insulin resistance

No (n=507)

0.358

0.466

0.176

0.408

Recessive

0.623 [0.423-0.916]

0.016

Si (n=600)

0.395

0.492

0.113

0.359

Co-dominant 2

0.610 [0.401-0.930]

0.021

Aspartate aminotransferase >p75

No (n=705)

0.368

0.471

0.161

0.396

Additive

0.777 [0.633-0.954]

0.016

Si (n=402)

0.391

0.500

0.109

0.359

Recessive

0.654 [0.431-0.992]

0.046

Co-dominant 2

0.583 [0.373-0.912]

0.018

rs181206

AA

AG

GG

Hyperuricemia

No (n = 887)

0.438

0.434

0.128

0.345

Additive

0.788 [0.624-0.994]

0.044

Si (n = 220)

0.495

0.405

0.100

0.302

rs40837

GG

GA

AA

Insulin resistance

No (n = 506)

0.354

0.470

0.176

0.411

Additive

0.646 [0.468-0.892]

0.008

Si (n = 601)

0.399

0.486

0.115

0.358

Dominant

0.616 [0.381-0994]

0.047

Recessive

0.488 [0.274-0.867]

0.014

Co-dominant 2

0.404 [0.210-0.775]

0.006

Aspartate aminotransferase >p75

No (n = 705)

0.366

0.475

0.159

0.396

Additive

0.771 [0.628-0.947]

0.013

Si (n = 402)

0.339

0.489

0.112

0.357

Co-dominant 2

0.603 [0.388-0.938]

0.025

Abbreviations: SNP, single nucleotide polymorphism; MAF, minor allele frequency; OR, odds ratio; 95% CI, 95% confidence interval; p75, percentile 75. All models were adjusted by age, gender and body mass index.

Co-dominant 1: comparison of heterozygote vs major allele homozygote genotype frequencies.

Co-dominant 2: comparison of heterozygote vs minor allele homozygote genotype frequencies.

Table 4: IL27p28 gene polymorphisms association with metabolic parameters in pCAD cases.

SNP

Genotype frequency

MAF

Model

OR [95% CI]

P

rs26528

CC

CT

TT

Aspartate aminotransferase>p75

No (n = 710)

0.382

0.459

0.159

0.389

Additive

0.802 [0.673-0.957]

0.014

Si (n = 452)

0.454

0.407

0.139

0.343

Dominant

0.689 [0.539-0.882]

0.003

Over-dominant

0.737 [0.576-0.944]

0.016

Co-dominant 1

0.678 [0.521-0.883]

0.004

rs181206

AA

AG

GG

Hyperuricemia

No (n = 745)

0.439

0.426

0.135

0.348

Recessive

0.546 [0.359-0.830]

0.005

Si (n = 417)

0.455

0.460

0.085

0.314

Co-dominant 2

0.566 [0.365-0.878]

0.011

rs40837

GG

GA

AA

Insulin resistance

No (n = 316)

0.393

0.470

0.137

0.372

Over-dominant

0.702 [0.504-0.978]

0.037

Si (n = 846)

0.414

0.433

0.153

0.369

Aspartate aminotransferase>p75

No (n = 710)

0.375

0.469

0.156

0.391

Additive

0.799 [0.670-0.953]

0.013

Si (n = 452)

0.456

0.407

0.137

0.341

Over-dominant

0.720 [0.563-0.922]

0.009

Co-dominant 1

0.662 [0.509-0.863]

0.002

Abbreviations: SNP, single nucleotide polymorphism; MAF, minor allele frequency; OR, odds ratio; 95% CI, 95% confidence interval; p75, percentile 75. All models were adjusted by age, gender and body mass index.

Co-dominant 1: comparison of heterozygote vs major allele homozygote genotype frequencies.

Co-dominant 2: comparison of heterozygote vs minor allele homozygote genotype frequencies.

Haplotype analysis

Only rs40837 and rs26528 polymorphisms were in high linkage disequilibrium (r2>0.956). Four different haplotypes were observed, but none of them showed a significant association with pCAD (data not shown).

IL-27 plasma levels

In a subsample, pCAD cases showed significantly higher IL-27 levels than control subjects (2.9 pg/mL vs 0.94 pg/mL, respectively; P = 0.004,). Both pCAD cases and healthy controls were non-obese individuals with hsCRP levels < 3 mg/L. IL-27 plasma levels were not significantly associated with any of the IL-27p28 polymorphisms analyzed here.

rs40837 G/A luciferase assays

Because rs40837 is a predicted 3’-UTR region target for miR-379-5p and/or miR-1225-5p, we examined the effect of both G and A alleles co-transfected with miR-379-5p or miR-1225-5p in a luciferase expression system. Cells co-transfected with the A allele construct and miR-379-5p showed a ~5% reduction in luciferase activity (P=0.0022) compared to those co-transfected with miR-Control. In contrast, cells co-transfected with the G allele construct (Figure 1B) and miR-379-5p showed higher mean luciferase activity (p = 0.0476) than those cotransfected with the miR-Control (Figure 1C). Moreover, in the presence of miR-379-5p, the A allele construct showed significantly lower luciferase activity than the G allele construct (p = 0.0087). Co-transfection of both constructs (A or G allele) with the miR-1225-5p did not significantly affect luciferase expression (data not shown).

DISCUSSION

As far as we know, this is the first study reporting the association of IL-27p28 gene polymorphisms with pCAD. We analyzed the distribution of rs26528, rs17855750, rs181206 and rs40837 in pCAD cases and controls in order to determine whether they confer susceptibility to pCAD. Polymorphism selection was based on informatics analyses for informativity (minor allele frequency >5%) and/or predicted functional effects. Both rs26528 T and rs40837 A alleles were significantly associated with a decreased risk of developing pCAD. IL-27p28 gene polymorphisms were also analyzed for associations with cardiometabolic risk factors, independently in cases and controls. rs26528 T and rs40837 A alleles were associated with a lower risk of high AST activity (>p75), rs181206 G was associated with a lower risk of hiperuricemia and rs40837 A with a lower risk of insulin resistance. The associations detected between the polymorphisms and cardiovascular parameters in both groups (cases and controls) were similar, suggesting that these associations are independent of the pathology present in these individuals.

Informatics analyses showed that rs40837 creates DNA binding sites for miR-379-5p and miR-1225-5p. In order to evaluate the functional effect of this polymorphism, luciferase assays were used to test the effect of each allele on gene expression in the presence of both miRNAs. In the presence of miR-379-5p, the rs40837 A allele showed significantly decreased luciferase gene expression. Several studies have provided evidence on the role of several miRNAs in atherosclerosis. They participate in the regulation of lipid metabolism, insulin biosynthesis, adipogenesis, endothelial dysfunction, neoangiogenesis, plaque development and rupture, as well as glucose homeostasis, among others [26-27]. Interestingly rs181206, associated with a lower risk of hyperuricemia, generates a binding site for SF/ASF proteins, which regulates alternative splicing [28]. This polymorphic site could regulate IL-27 isoforms relevant for the development of hiperuricemia.

IL-27 plays an important role in inflammation and atherosclerosis pathogenesis with dual effects, both pro and anti-inflammatory [20, 21]. It promotes early Th1 cell differentiation [16] and suppresses Th2 [18] and Th17 [19] differentiation. The type of effect on inflammation exerted by IL-27 has been reported to differ in various diseases. While it is known to promote inflammation in hepatitis [29] and systemic sclerosis [30], it suppresses inflammation in autoimmune arthritis [31], allergic asthma [32] and autoimmune encephalomyelitis [19]. Moreover, studies on the role of IL-27 in atherosclerosis in animal models and humans have shown inconsistent results. In the murine model, IL-27 administration suppressed macrophage activation and atherosclerosis development [23] and mice deficient for IL-27 or its receptor showed increased atherosclerosis susceptibility [23], suggesting IL-27 has a protective role. In contrast, coronary patients showed significantly higher IL-27 levels, which correlated with the severity of stenosis [24], suggesting a pro-atherogenic role in humans. Similar to the findings of Jin et al. [24], significantly higher IL-27 plasma concentrations were observed in Mexican pCAD cases as compared to controls.

We consider that the main strengths of this study are the following: a) The study included a large cohort of Mexican cases and controls with thorough phenotyping, and with tomographic, clinical and biochemical data, allowing to adjust our analyses for a large number of potential confounders; b) Controls included only individuals without tomographic evidence of SA (CAC score=0); c) Population stratification was ruled out as a potential confounding factor, because the proportions of Caucasian, Native American and African ancestry were similar in cases and controls; and d) A functional effect of rs40837 polymorphism on luciferase expression was observed, which was in accordance with the observed associations. Nevertheless, results should be interpreted with caution, considering the following limitations. First, due to the transversal character of the study, conclusions on causality cannot be made. Second, because the selection of participants was not random, the findings may not be applicable to the general population. However, considering that the participants have no knowledge of their genotypes, the genotype distributions would be expected to be similar in a randomly selected sample. Third, insulin resistance was not evaluated using euglycemic/hyperinsulinemic clamp, nonetheless HOMA-IR index has proven to be a reliable measurement of insulin sensitivity [33].

This study shows for the first time, that IL-27p28 gene polymorphisms are associated with pCAD, AST activity, hyperuricemia and insulin resistance in the Mexican population. Despite all the evidence on the role of IL-27 in atherosclerosis, to best of our knowledge, to date only one cross sectional analysis has evaluated the role of IL-27p28 gene variants in cardiovascular disease. In this study, four IL-27p28 tag SNPs, (rs181206, rs17855750, rs37833 and rs153109) were determined in a large number of CAD cases belonging to the GeneID Chinese Han population [34]. After adjusting for confounder’s variables, the polymorphisms were not associated with CAD, age at disease onset or severity [34]. These results are consistent with the findings of the present study for the rs181206 and rs17855750 IL-27p28 gene variants. Thus, replications of the associations here reported (rs26528, rs40837) should be sought in other cohorts to confirm these results. Our results suggest rs26528T and rs40837A alleles could be considered as potential susceptibility markers for pCAD and insulin resistance in our population. Although these polymorphisms were not significantly associated with IL-27 plasma levels, pCAD cases showed significantly higher IL-27 levels than control subjects as previously reported [24]. The fact that the IL-27 plasma levels were not significantly associated with any of the IL-27p28 polymorphisms analyzed here could be explained considering that like other molecules, the production of IL-27 include a complex mechanism that involve not only changes at DNA level but also epigenetic modifications. Moreover, is important to consider that in our study the levels of IL-27 were measured only in a subsample of pCAD cases and controls with specific characteristics.

Because the Mexican population has particular and different genetic characteristics to other ethnic groups [35-38], the IL-27p28 polymorphism associations observed here should be sought in other populations in order to establish if they are specific for the Mexican population or are shared with other ethnic groups.

MATERIALs AND METHODS

Subjects

The GEA Mexican Study was designed to examine the genetic bases of pCAD and the relationship between traditional and emerging risk factors of SA in an adult Mexican population. This study included 1200 pCAD cases and 1500 healthy individuals as control group aged 30 to 75 years. All participants were unrelated and of self-reported Mexican mestizo ancestry for 3 generations. pCAD was defined as history of myocardial infarction, angioplasty, revascularization surgery or coronary stenosis >50% on angiography, diagnosed before age 55 in men and before age 65 in women. Patients with acute cardiovascular events 3 months prior to the selection were excluded. Controls were apparently healthy asymptomatic individuals without personal or family history of pCAD, recruited from blood bank donors and through brochures posted in Social Services centers. Exclusion criteria for controls included congestive heart failure; liver, renal, thyroid or oncological disease. Standardized questionnaires were applied to all participants to obtain demographic information, family medical history, history of nutritional habits, physical activity, alcohol consumption and pharmacological treatment. The GEA study was approved by Bioethics Committee of the Instituto Nacional de Cardiología Ignacio Chávez (INCICH), and aligned to Helsinki’s Declaration. All participants provided informed consent.

Anthropometric and biochemical measurements

BMI was calculated as weight in kilograms divided by height in meters squared. Waist circumference was measured using a glass fiber measuring tape in the middle point of the distance between the lower side of the waist and the iliac crest. Blood pressure was measured at rest 3 times using a digital 5200 series Welch Allyn sphygmomanometer (Shaneateies Fails, N.Y., USA.) and the last two measurements were averaged. Venous blood samples were obtained after a 12-hour fast, and all biochemical measurements were performed at the Endocrinology Laboratory of the INCICH using standardized procedures as previously described [39-41].

Computed axial tomography study

Computed tomography of the chest and abdomen were performed using a 64-channel multi-detector helical computed tomography system (Somatom Sensation, Siemens) and interpreted by experienced radiologists. Scans were read to assess and quantify the following: 1) CAC score using the Agatston method [42]; 2) TAF, subcutaneous and visceral abdominal fat areas (SAF and VAF) as described by Kvist [43]; and 3) hepatic to splenic attenuation ratio as described by Longo et al [44]. SA was defined as the presence CAC score>0. All pCAD cases and healthy controls underwent computed tomography. Of the 1500 apparently healthy controls, 393 subjects had a CAC score above cero, therefore they were not considered for the present analysis and were thus considered as individuals with SA. The final control group included 1107 individuals (CAC scores = 0).

Definition of risk factors

Metabolic and cardiovascular risk factors were evaluated in both pCAD cases and controls and defined as previously described [39-41].

Genetic analysis and functional prediction

Functional prediction of IL-27p28 single nucleotide polymorphisms (SNPs) was sought using bioinformatics tools including FastSNP, SplicePort: An Interactive Splice Site Analysis Tool (http://spliceport.cbcb.umd.edu/SplicingAnalyser.html), SNPs3D (http://www.snps3d.org/), PESX: Putative Exonic Splicing enharcers/Silencers (http://cubweb.biology.columbia.edu/pesx/), and ESEfinder release 3.0 (http://rulai.cshl.edu/cgi-bin/tools/ESE3/esefinder.cgi).

Four IL-27p28 gene polymorphisms with possible functional consequences and/or minor allele frequencies > 5% were selected for analysis: rs17855750 and rs181206 introduce binding sites for transcriptional factors SF2/ASF2; rs40837 modifies binding sites for miRNAs. rs26528, although not predicted as functional, was informative and in high linkage disequilibrium with rs40837 (r2=0.956) and was thus included in the study. Genomic DNA was extracted from peripheral blood using standard methods. All SNPs were genotyped using TaqMan assays on a real-time PCR Prism 7900HT Fast Real-Time PCR system (Applied Biosystems, Foster City, CA, USA) and analyzed by the allelic exclusion program. Samples previously sequenced of the different genotypes of the polymorphisms studied were included as positive controls.

In order to gain insight into the possible functional implications of the 3’UTR rs40837 polymorphism, the TargetScan (http://www.targetscan.org), Diana-MicroT3.0 (http://diana.cslab.ece.ntua.gr) and miRanda (http://www.microrna.org) software’s were used, revealing that this SNP could be a target for miR-379-5p and miR-1225-5p.

Estimation of ancestry

Because the Mexican-Mestizo population is admixed, in order to assess the possible influence of population stratification, a panel of 265 ancestry informative markers distinguishing mainly Amerindian, European and African ancestry were selected [45] and genotyped on Illumina BeadStation using the GoldenGate assay. Duplicate control samples were genotyped on each chip, which also served as internal controls for quality of clustering and reproducibility. The primary analysis of the genotyping data with the Illumina GenomeStudio software v.2011.1 was followed by visual inspection and assessment of data quality and clustering. Genotyping accuracy was also assessed by genotype clustering using the Illumina GeneTrain score, which is a measure of the clustering confidence of individual SNP alleles. Global Caucasian, Amerindian and African ancestry were determined using the ADMIXTURE software. Mean global ancestry was not significantly different between cases and controls (55.8% vs 54.0 % Amerindian ancestry, 34.3% vs 35.8% Caucasian and 9.8% vs 10.1% African mean ancestry for cases and controls respectively, P>0.05), strongly suggesting that population stratification was not a bias or confounding factor in this study.

Measurement of IL-27 plasma levels

IL-27 plasma concentrations were measured in a carefully selected subsample of 450 pCAD cases and 450 healthy controls (without obesity and hsCRP < 3 mg/L), using a Bioplex system (Bio-Rad, Contra Costa County, State of California, USA) according to manufacturer’s instructions.

Reporter constructs, transfection, and luciferase assays

The 266 bp human IL27 3´ UTR containing the rs40837 A or G allele was amplified with the forward primer 5´-GCGCACGCGTCCCCCACCCTTTAGAACTTT-3´ and the reverse primer 5´-GCGCAAGCTTTGGATGAGAGTGCTTTATTGG-3´ from a homozygous human genomic DNA sample. PCR products were separated on agarose gels, purified and cloned into pMIR-REPORT plasmids (Applied Biosystems, Foster City, CA, USA) with MluI and HindIII digestion (Figure 1A and 1B). HEK293 cells were grown in Dulbecco’s modify Eagle’s medium (DMEM) (InvitrogenTM) supplemented with 10% fetal bovine serum and 1% antibiotics-antimycotics (InvitrogenTM) at 37oC with 5% CO2. A total of 100,000 cells were plated into each well of 12-well plates in Opti-MEM® serum free medium (Invitrogen Life Technologies, Inc., Carlsbad, CA, USA). Forty-eight hours after plating, cells were co-transfected using Lipofectamine 2000 (Invitrogen, Life Technologies, Carlsbad, CA) according to the manufacturer’s instructions. Each co-transfection reaction contained 500 ng of pMIR-REPORT (rs40837) A or G allele vector plus 100 ng pRL/CMV Renilla reniformis luciferase vector plasmids that served as a transfection internal control. Twenty-four hours after co-transfection of the plasmids, 75 nM negative control (scrambled sequence) or miR-379-5p RNA (Dharmacon GE Life Sciences) were transfected with siPORT amine transfection agent (Applied Biosystems, Foster City, CA, USA). Forty-eight hours after miRNA transfection, both firefly and Renilla luciferase activities were quantified by a dual-luciferase reporter assay system (Promega, Madison, WI). The relative luciferase activity was calculated according to the manufacturer’s instructions in a TD-20/20 luminometer (Turner BioSystems, Sunnyvale, CA).

Statistical analysis

Categorical variables were compared between groups using the Chi square test, continuous variables were compared with Student’s t test or Mann Whitney’s U test for parametric and non-parametric variables, respectively. Allele and genotype frequencies were estimated by direct counting. Hardy-Weinberg’s equilibrium was tested using the Chi square test. For all the studied variants, statistical power to detect association with pCAD was greater than 90% as estimated with QUANTO software [http:hydra.usc.edu/GxE/]. Multivariate logistic regression analysis was used to analyze associations with pCAD under different inheritance models: additive (major allele homozygotes vs. heterozygotes vs. minor allele homozygotes), co-dominant (major allele homozygotes vs. heterozygotes and major allele homozygotes vs. minor allele homozygotes), dominant (major allele homozygotes vs. heterozygotes + minor allele homozygotes), over-dominant (heterozygotes vs. major allele homozygotes + minor allele homozygotes), and recessive (major allele homozygotes + heterozygotes vs. minor allele homozygotes). The models were adjusted for age, gender, BMI, smoking habit, TAF, HOMA-IR, AST, adiponectin and uric acid levels as appropriate. Logistic regression analyses were performed to assess associations of IL-27 SNPs with metabolic parameters and cardiovascular risk factors, under different inheritance models and adjusting for age, gender and BMI, as appropriate. Linkage disequilibrium and haplotype analysis were performed with Haploview software (version 4.1, Broad Institute of Massachusetts Institute of Technology and Harvard University, Cambridge, MA, USA). P < 0.05 values were considered statistically significant. All analyses were performed using SPSS software v15.0 (SPSS Chicago, IL).

Abbreviations

ALT: alanine aminotransferase, AST: aspartate aminotransferase, BMI: body mass index, CAC: coronary artery calcification, CAD: coronary artery disease, CI: confidence interval, CVD: cardiovascular disease, hsCRP: high sensitivity C reactive protein, HOMA-IR: homeostasis model assessment of insulin resistance, IL-27: interleukin 27, LDL-C: low density lipoprotein-cholesterol, MAF: minor allele frequency, OR= odds ratio, oxLDL: oxidized LDL, pCAD: premature coronary artery disease, SA: subclinical atherosclerosis, SAF: subcutaneous abdominal fat, SNP: single nucleotide polymorphism, VAF: visceral abdominal fat, TAF: total abdominal fat, TC: total cholesterol, T2DM: type 2 diabetes mellitus.

Author contributions

RPS and GVA designed the study, were responsible of the data analysis and extraction, and wrote the manuscript with final comment and approval by all the authors. NPH, JMRP, RMCV, BRR, LL and GL were responsible for the acquisition, analysis and interpretation of data and advice regarding study design. CFD, TVM and CPR have been involved in the critical revision of the manuscript for important intellectual content. All authors were involved in the drafting of the manuscript. All authors read and approved the final manuscript version to be published.

ACKNOWLEDGMENTS

Rosalinda Posadas-Sánchez is a doctoral student from Programa de Doctorado en Ciencias Biomédicas, Universidad Nacional Autónoma de México (UNAM)

The authors are grateful to the study participants and to Angel Rene López-Uribe for his technical assistance.

CONFLICTs OF INTEREST

No conflicts to disclose.

FUNDING

This work was supported in part by grants from the Consejo Nacional de Ciencia y Tecnología (Project number SALUD-2014-1-233727).

REFERENCES

1. Leeper NJ, Kullo IJ, Cooke JP. Genetics of peripheral artery disease. Circulation. 2012; 125 :3220−3228.

2. Raman K, Chong M, Akhtar-Danesh GG, D’Mello M, Hasso R, Ross S, Xu F, Paré G. Genetic markers of inflammation and their role in cardiovascular disease. Can J Cardiol. 2013; 29: 67−74.

3. McPherson R, Davies RW. Inflammation and coronary artery disease: insights from genetic studies. Can J Cardiol. 2012; 28: 662−666.

4. Hansson GK. Inflammation, atherosclerosis, and coronary artery disease. N Engl J Med. 2005; 352: 1685−1695.

5. Westerterp M, Murphy AJ, Wang M, Pagler TA, Vengrenyuk Y, Kappus MS, Gorman DJ, Nagareddy PR, Zhu X, Abramowicz S, Parks JS, Welch C, Fisher EA, et al. Deficiency of ATP-binding cassette transporters A1 and G1 in Macrophages Increases Inflammation and Accelerates Atherosclerosis in Mice. Circ Res. 2013; 112: 1456−1465.

6. Chinetti-Gbaguidi G, Staels B. Macrophage polarization in metabolic disorders: functions and regulation. Curr Opin Lipidol. 2011; 22: 365−372.

7. Tuttolomondo A, Di Raimondo D, Pecoraro R, Arnao V, Pinto A, Licata G. Atherosclerosis as an inflammatory disease. Curr Pharm Des. 2012; 18: 4266−4288.

8. Salem MK, Vijaynagar B, Sayers RD, West K, Moore D, Robinson TG, Naylor AR, Bown MJ. Histologically unstable asymptomatic carotid plaques have altered expression of genes involved in chemokine signalling leading to localized plaque inflammation and rupture. Eur J Vasc Endovasc Surg. 2013; 45: 121−127.

9. de Boer OJ, van der Wal AC, Verhagen CE, Becker AE. Cytokine secretion profiles of cloned T cells from human aortic atherosclerotic plaques. J Pathol. 1999; 188: 174−179.

10. Davenport P, Tipping PG. The role of interleukin-4 and interleukin-12 in the progression of atherosclerosis in apolipoprotein E-deficient mice. Am J Pathol. 2003; 163: 1117−1125.

11. Whitman SC, Ravisankar P, Daugherty A. Interleukin-18 enhances atherosclerosis in apolipoprotein E(-/-) mice through release of interferon-gamma. Circ Res. 2002; 90: E34−E38.

12. Buono C, Binder CJ, Stavrakis G, Witztum JL, Glimcher LH, Lichtman AH. T-bet deficiency reduces atherosclerosis and alters plaque antigen-specific immune responses. Proc Natl Acad Sci USA. 2005; 102: 1596−1601.

13. Hansson GK, Libby P. The immune response in atherosclerosis: a double-edged sword. Nat Rev Immunol. 2006; 6: 508-519.

14. Pflanz S, Timans JC, Cheung J, Rosales R, Kanzler H, Gilbert J, Hibbert L, Churakova T, Travis M, Vaisberg E, Blumenschein WM, Mattson JD, Wagner JL, et al. IL-27, a heterodimeric cytokine composed of EBI3 and p28 protein, induces proliferation of naive CD4+ T cells. Immunity. 2002; 16: 779−790.

15. Yoshida H, Nakaya M, Miyazaki Y. Interleukin 27: a double-edged sword for offense and defense. J Leukoc Biol. 2009; 86: 1295−1303.

16. Takeda A, Hamano S, Yamanaka A, Hanada T, Ishibashi T, Mak TW, Yoshimura A, Yoshida H. Cutting edge: role of IL-27/WSX-1 signaling for induction of T-bet through activation of STAT1 during initial Th1 commitment. J Immunol. 2003; 170: 4886−4890.

17. Artis D, Villarino A, Silverman M, He W, Thornton EM, Mu S, Summer S, Covey TM, Huang E, Yoshida H, Koretzky G, Goldschmidt M, Wu GD, et al. The IL-27 receptor (WSX-1) is an inhibitor of innate and adaptive elements of type 2 immunity. J Immunol. 2004; 173: 5626−5634.

18. Yoshimoto T, Yoshimoto T, Yasuda K, Mizuguchi J, Nakanishi K. IL-27 suppresses Th2 cell development and Th2 cytokines production from polarized Th2 cells: a novel therapeutic way for Th2-mediated allergic inflammation. J Immunol. 2007; 179: 4415−4423.

19. Batten M, Li J, Yi S, Kljavin NM, Danilenko DM, Lucas S, Lee J, de Sauvage FJ, Ghilardi N. Interleukin 27 limits autoimmune encephalomyelitis by suppressing the development of interleukin 17-producing T cells. Nat Immunol. 2006; 7: 929−936.

20. Stumhofer JS, Laurence A, Wilson EH, Huang E, Tato CM, Johnson LM, Villarino AV, Huang Q, Yoshimura A, Sehy D, Saris CJ, O’Shea JJ, Hennighausen L, et al. Interleukin 27 negatively regulates the development of interleukin 17-producing T helper cells during chronic inflammation of the central nervous system. Nat Immunol. 2006; 7: 937−945.

21. Awasthi A, Carrier Y, Peron JP, Bettelli E, Kamanaka M, Flavell RA, Kuchroo VK, Oukka M, Weiner HL. A dominant function for interleukin 27 in generating interleukin 10-producing anti-inflammatory T cells. Nat Immunol. 2007; 8: 1380−1389.

22. Kempe S, Heinz P, Kokai E, Devergne O, Marx N, Wirth T. Epstein-barr virus-induced gene-3 is expressed in human atheroma plaques. Am J Pathol. 2009; 175: 440−447.

23. Hirase T, Hara H, Miyazaki Y, Ide N, Nishimoto-Hazuku A, Fujimoto H, Saris CJ, Yoshida H, Node K. Interleukin 27 inhibits atherosclerosis via immunoregulation of macrophages in mice. Am J Physiol Heart Circ Physiol. 2013; 305: H420−H429.

24. Jin W, Zhao Y, Yan W, Cao L, Zhang W, Wang M, Zhang T, Fu Q, Li Z. Elevated circulating interleukin-27 in patients with coronary artery disease is associated with dendritic cells, oxidized low-density lipoprotein, and severity of coronary artery stenosis. Mediators Inflamm. 2012; 2012: 506283.

25. Chae SC, Li CS, Kim KM, Yang JY, Zhang Q, Lee YC, Yang YS, Chung HT. Identification of polymorphisms in human interleukin-27 and their association with asthma in a Korean population. J Hum Genet. 2007; 52: 355−361.

26. Fernández-Hernando C, Ramírez CM, Goedeke L, Suárez Y. MicroRNAs in metabolic disease. Arterioscler Thromb Vasc Biol. 2013; 33 :178−185.

27. Madrigal-Matute J, Rotllan N, Aranda JF, Fernández-Hernando C. MicroRNAs and aterosclerosis. Curt Atheroscler Rep. 2013; 15: 322.

28. Sureau A, Gattoni R, Dooghe Y, Stévenin J, Soret J. SC35 autoregulates its expression by promoting splicing events that destabilize its mRNAs. EMBO J. 2001; 20: 1785−1796.

29. Siebler J, Wirtz S, Frenzel C, Schuchmann M, Lohse AW, Galle PR, Neurath MF. Cutting edge: a key pathogenic role of IL-27 in T cell- mediated hepatitis. J Immunol. 2008; 180: 30−33.

30. Yoshizaki A, Yanaba K, Iwata Y, Komura K, Ogawa A, Muroi E, Ogawa F, Takenaka M, Shimizu K, Hasegawa M, Fujimoto M, Sato S. Elevated serum interleukin-27 levels in patients with systemic sclerosis: association with T cell, B cell and fibroblast activation. Ann Rheum Dis. 2011; 70: 194−200.

31. Rajaiah R, Puttabyatappa M, Polumuri SK, Moudgil KD. Interleukin-27 and interferon-gamma are involved in regulation of autoimmune arthritis. J Biol Chem. 2011; 286: 2817−2825.

32. Miyazaki Y, Inoue H, Matsumura M, Matsumoto K, Nakano T, Tsuda M, Hamano S, Yoshimura A, Yoshida H. Exacerbation of experimental allergic asthma by augmented Th2 responses in WSX-1-deficient mice. J Immunol. 2005; 175: 2401−2407.

33. Bonora E, Targher G, Alberiche M, Bonadonna RC, Saggiani F, Zenere MB, Monauni T, Muggeo M. Homeostasis model assessment closely mirrors the glucose clamp technique in the assessment of insulin sensitivity: studies in subjects with various degrees of glucose tolerance and insulin sensitivity. Diabetes Care. 2000; 23: 57−63.

34. Fan Q, Nie S, Li S, Liao Y, Zhang H, Zha L, Wang F, Tang T, Xia N, Xu C, Wang P, Xie T, Xie J, et al. Analysis of the genetic association between IL27 variants and coronary artery disease in a Chinese Han population. Sci Rep. 2016; 6: 25782.

35. Lisker R, Perez-Briceño R, Granados J, Babinsky V, de Rubens J, Armendares S, Buentello L. Gene frequencies and admixture estimates in a Mexico City population. Am J Phys Anthropol. 1986; 71: 203−207.

36. Lisker R, Pérez-Briceno R, Granados J, Babinsky V. Gene frequencies and admixture estimates in the state of Puebla, Mexico. Am J Phys Anthropol. 1988; 76: 331−335.

37. Lisker R, Ramirez E, Briceño RP, Granados J, Babinsky V. Gene frequencies and admixture estimates in four Mexican urban centers. Hum Biol. 1988; 62: 791−801.

38. Juárez-Cedillo T, Zuñiga J, Acuña-Alonzo V, Pérez-Hernández N, Rodríguez-Pérez JM, Barquera R, Gallardo GJ, Sánchez-Arenas R, García-Peña Mdel C, Granados J, Vargas-Alarcón G. Genetic admixture and diversity estimations in the Mexican Mestizo population from Mexico City using 15 STR polymorphic markers. Forensic Sci Int Genet. 2008; 2: e37−39.

39. Posadas-Sánchez R, López-Uribe ÁR, Posadas-Romero C, Pérez-Hernández N, Rodríguez-Pérez JM, Ocampo-Arcos WA, Fragoso JM, Cardoso-Saldaña G, Vargas-Alarcón G. Association of the I148M/PNPLA3 (rs738409) polymorphism with premature coronary artery disease, fatty liver, and insulin resistance in type 2 diabetic patients and healthy controls. The GEA study. Immunobiology. 2016 Sep 3. [Epub ahead of print].

40. Medina-Urrutia A, Posadas-Romero C, Posadas-Sánchez R, Jorge-Galarza E, Villarreal-Molina T, González-Salazar M del C, Cardoso-Saldaña G, Vargas-Alarcón G, Torres-Tamayo M, Juárez-Rojas JG. Role of adiponectin and free fatty acids on the association between abdominal visceral fat and insulin resistance. Cardiovasc Diabetol. 2015; 14: 20.

41. Posadas-Sánchez R, Ocampo-Arcos WA, López-Uribe AR, González-Salazar M-C, Cardoso-Saldaña G, Mendoza-Pérez E, Medina-Urrutia A, Jorge-Galarza E, Posadas-Romero C. Asociación del ácido úrico con factores de riesgo cardiovascular y aterosclerosis subclínica en adultos mexicanos. Rev Mex Endocrinol Metabol Nutr. 2014; 1: 14−21.

42. Mautner GC, Mautner SL, Froehlich J, Feuerstein IM, Proschan MA, Roberts WC, Doppman JL. Coronary artery calcification: assessment with electron beam CT and histomorphometric correlation. Radiology. 1994; 192: 619-623.

43. Kvist H, Chowdhury B, Grangård U, Tylén U, Sjöström L. Total and visceral adipose-tissue volumes derived from measurement with computed tomography in adult men and women: predicted equations. Am J Clin Nutr. 1988; 48: 1351−1361.

44. Longo R1, Ricci C, Masutti F, Vidimari R, Crocé LS, Bercich L, Tiribelli C, Dalla Palma L. Fatty infiltration of the liver. Quantification by 1H localized magnetic resonance spectroscopy and comparison with computed tomography. Invest Radiol. 1993; 28: 297-302.

45. Silva-Zolezzi I, Hidalgo-Miranda A, Estrada-Gil J, Fernandez-Lopez JC, Uribe-Figueroa L, Contreras A, Balam-Ortiz E, del Bosque-Plata L, Velazquez-Fernandez D, Lara C, Goya R, Hernandez-Lemus E, Davila C, et al. Analysis of genomic diversity in Mexican Mestizo populations to develop genomic medicine in Mexico. Proc Natl Acad Sci U S A. 2009; 106: 8611−8616.


Creative Commons License All site content, except where otherwise noted, is licensed under a Creative Commons Attribution 4.0 License.
PII: 16223