Genetic studies in the Pakistani population reveal novel associations with ventricular septal defects (VSDs).

BACKGROUND: With prevalence up to 4%, Ventricular Septal Defect (VSD) is one of the leading causes of neonatal deaths. VSD is a common complex genetic disorder that has been associated with many genetic determinants. Variants from genes for the transcription factors including T-Box TBX5 and NFATc1 (nuclear factor of activated T cells, cytoplasmic 1), Vascular endothelial growth factor (VEGF), ISLET1 (encoded by the ISL1 gene) and enzyme MTHFR, a methylene tetrahydrofolate reductase were selected. Genetic risk score (GRS) is a widely accepted approach used to convert the genetic data into prediction and assessment tool for disease susceptibility. METHODS: A total of 200 participants were recruited for the current study, 100 VSD patients and 100 controls. Genotyping of the ISL1: rs1017, NFATc1: rs7240256, VEGF: rs36208048, TBX5: rs11067075, and MTHFR: rs1801133 variants was performed using tetra primer ARMS PCR and PCR-RFLP. For the statistical analysis, the software SPSS version 23 was used. Genotypic frequencies of cases and controls were calculated using chi-square (χ²) whereas allelic frequencies were calculated by using the SNPStats tool. The association of GRS quartiles with VSD was examined using binary logistic regression. Adjusted p-value 0.01 was used as significance threshold for all analyses. RESULTS: The ISL1 (OD: 0.242, CI: 0.158-0.37, p-value: 2.15 × 10- 4 :), NFATc1 (OD: 2.53, CI: 1.64-3.89, p-value: 2.11 × 10- 5), TBX5 (OD: 2.24, CI: 1.47-3.41, p-value:1.6 × 10- 4) and MTHFR (OD: 10.46, CI: 5.68-19.26, p-value: 2.09 × 10- 9:) variants were found to be in association with VSD. In contrast, the VEGF (OD: 0.952, CI: 0.56-1.62, p-value: 0.8921) variant did not show significance association with the VSD. For cases, the mean GRS score was 3.78 ± 1.285 while in controls it was 2.95 ± 1.290 (p-value: 0.479, CI: 0.474-1.190). Comparison of GRS between cases and control showed that mean GRS of cases was 1.90 ± 0.480 while in controls it was 1.68 ± 0.490 (p-value: 0.001, CI: 0.086-0.354). Higher quartiles were more prevalent in cases whereas lower quartiles were more prevalent in controls. CONCLUSION: GRS of these five loci was strongly associated with VSD. Moreover, genetic risk score can provide better information for the association between variants and disease as compared to a single SNP. We also illustrated that the cumulative power of GRS is greater over the single SNP effect. This is a pilot scale study with a relatively small sample size whose findings should be replicated in a larger sample size for the unique local Pakistani population.

Study of variants associated with ventricular septal defects (VSDs) highlights the unique genetic structure of the Pakistani population.

BACKGROUND: Ventricular septal defects (VSDs) are one of the leading causes of death due to cardiac anomalies during the first months of life. The prevalence of VSD in neonates is reported up to 4%. Despite the remarkable progress in medication, treatment and surgical procedure for VSDs, the genetic etiology of VSDs is still in infancy because of the complex genetic and environmental interactions. METHODS: Three hundred fifty subjects (200 VSD children and 150 healthy controls) were recruited from different pediatric cardiac units. Pediatric clinical and demographic data were collected. A total of six variants, rs1017 (ISL1), rs7240256 (NFATc1), rs36208048 (VEGF), variant of HEY2, rs11067075 (TBX5) and rs1801133 (MTHFR) genes were genotyped by tetra-ARMS PCR and PCR-RFLP methods. RESULTS: The results showed that in cases, the rs1017 (g.16138A > T) variant in the ISL1 gene has an allele frequency of 0.42 and 0.58 respectively for the T and A alleles, and 0.75 and 0.25 respectively in the controls. The frequencies of the AA, TA and TT genotypes were, 52%, 11% and 37% in cases versus 21%, 8% and 71% respectively in the controls. For the NFATc1 variant rs7240256, minor allele frequency (MAF) was 0.43 in cases while 0.23 in controls. For the variant in the VEGF gene, genotype frequencies were 0% (A), 32% (CA) and 68% (CC) in cases and 0.0%, 33% and 67% respectively in controls. The allele frequency of C and A were 0.84 and 0.16 in cases and 0.83 and 0.17 respectively in controls. The TBX5 polymorphism rs11067075 (g.51682G > T) had an allelic frequency of 0.44 and 0.56 respectively for T and G alleles in cases, versus 0.26 and 0.74 in the controls. We did not detect the presence of the HEY2 gene variant (g.126117350A > C) in our pediatric cohort. For the rs1801133 (g.14783C > T) variant in the MTHFR gene, the genotype frequencies were 25% (CC), 62% (CT) and 13% (TT) in cases, versus 88%, 10% and 2% in controls. The ISL1, NFATc1, TBX5 and MTHFR variants were found to be in association with VSD in the Pakistani pediatric cohort whilst the VEGF and HEY2 variants were completely absent in our cohort. CONCLUSION: We propose that a wider programme of genetic screening of the Pakistani population for genetic markers in heart development genes would be helpful in reducing the risk of VSDs.

Association of variable number of tandem repeats (VNTR) and T941G polymorphism of monoamine oxidase (MAO-A) gene with aggression in Pakistani subjects.

BACKGROUND: Human behavioral traits are known to be significantly heritable. Certain individuals have a greater tendency of negative behavioral aspects including aggression. The quest to identify tunderlying genetic causes has led to identification of a number of genetic markers, one of them is the monoamine oxidase-A (MAO-A) gene. OBJECTIVE: We aimed to genotype a variable number of tandem repeats (VNTRs) in the promoter region and a functional SNP within this gene (T941G, dbSNP ID: rs6323) in the recruited cohort of 482 subjects. METHODS: After DNA isolation, genotyping was done by PCR-RFLP and the results were confirmed by sequencing. RESULTS: For VNTRs, the results showed, highest frequency of 3.5 repeats in males and 4 repeats in females in the promoter region. The genotype frequencies for the SNP in cases were GG=16.3%, TG=20.6% and TT=63.1%, while in controls, the frequencies were GG=12.7%, TG=6.3%, and TT=81.0%. The allele frequencies were significantly different between cases and controls (p=0.015; OR=1.51; CI=1.085-2.102). CONCLUSION: The selected VNTR and SNP appeared to be significantly associated with aggression. These VNTRs and SNP have not been studied previously in the Pakistani population, hence they represent a unique ethnic group. These results, however, would have to be replicated in larger cohorts.

First report of polymorphisms in MTRR, GATA4, VEGF, and ISL1 genes in Pakistani children with isolated ventricular septal defects (VSD).

BACKGROUND: Ventricular septal defects (VSDs) are malformations in the septum separating the heart's ventricles. VSDs may present as a single anomaly (isolated/nonsyndromic VSD) or as part of a group of phenotypes (syndromic VSD). The exact location of the defect is crucial in linking the defect to the underlying genetic cause. The number of children visiting cardiac surgery units is constantly increasing. However, there are no representative data available on the genetics of VSDs in Pakistani children. METHODS: Two hundred forty-two subjects (121 VSD children and 121 healthy controls) were recruited from pediatric cardiac units of Lahore. The clinical and demographic data of the subjects were collected. A total of four SNPs, one each from MTRR, GATA4, VEGF, and ISL1 genes were genotyped by PCR-RFLP. RESULTS: The results showed that the minor allele (T) frequency (MAFs) for the MTRR gene variant rs1532268 (c.524C > T) was 0.20 and 0.41 in the controls and the cases, respectively, with the genotype frequencies 3, 35, 62% in the controls and 12, 59 and 29% in the cases for TT, CT, CC genotypes, respectively (allelic OR: 5.73, CI: 3.82-8.61, p-value: 5.11 × 10- 7). For the GATA4 variant rs104894073 (c.886G > A), the MAF for the controls and the cases was 0.16 and 0.37, respectively, the frequencies of AA, GA and GG genotypes were 2, 28, and 70% in the controls and 5, 64 and 31% of the cases (allelic OR: 3.08, CI: 2.00-4.74, p-value: 8.36 × 10- 8). The rs699947 (c.-2578C > A) of VEGF gene showed MAF 0.36 and 0.53 for the controls and cases, respectively, with the genotype frequencies 13, 42, and 45% in the controls and 22, 15, and 63% in the cases for the AA, CA, CC (allelic OR: 2.03, CI: 1.41-2.92, p-value: 0.0001). The ISL1 gene variant rs6867206 (g.51356860 T > C), the MAFs were 0.26 and 0.31 in the controls and cases, respectively. The genotype frequencies were 48, 52, 0% in the controls and 39, 61, 0% in the cases for TT, TC, CC genotypes (allelic OR: 0.27, CI: 0.85-1.89, p-value: 0.227). The MTRR, GATA4 and VEGF variants showed association while ISL1 variant did not appear to be associated with the VSD in the recruited cohort. CONCLUSION: This first report in Pakistani children demonstrates that single nucleotide polymorphisms in genes encoding transcription factors, signaling molecules and structural heart genes involved in fetal heart development are associated with developmental heart defects., however further work is needed to validate the results of the current investigation.

The abnormal lipid profile in obesity and coronary heart disease (CHD) in Pakistani subjects.

BACKGROUND: Obesity has become global epidemic in the last three decades, whereas Coronary Heart Disease (CHD) still remains the most important cause of mortality in the world. The study was aimed at determining the pattern of lipid profile for the obese and CHD population in Pakistan. As obesity is a strong predisposing risk factor for CHD, we aimed to analyze the lipid parameters in both conditions and compare them with the healthy controls of the same ethnicity. METHODS: Blood samples were collected from one thousand individuals (500 with CHD, 250 with obesity, 250 healthy controls). The lipid profile (total Cholesterol, triglycerides, HDL-C, LDL-C and VLDL) was measured using commercially available kits. The pattern of dyslipidemia was then studied by comparing the results in both groups with controls as well as population cutoffs. The quantitative variables were checked for normality and log transformation was done for variables where appropriate. Analysis of variance and logistic regression were done to check the association of lipid parameters with obesity and CHD. RESULTS: The obese and CHD groups showed a dyslipidemic profile than the healthy controls. CHD group had a higher proportion of CHD in any of the first degree blood relatives (36.0% vs. 1.8%), a similar trend was observed in the obese group, where 63.9% cases had positive family history. Among cases, 50.7% had combined lipid abnormalities, i.e., the values of TC, LDL-C, TG and HDL-C, all were deranged. Whereas 49.52% had TC more than normal cut off (> 200 mg/dl), 51.6% had LDL-C > 100 mg/dl. Similarly, 80.4% of patients had TG levels more than upper normal range (> 150 mg/dl) and 64% had HDL values in moderate CHD risk group (< 50 mg/dl). The results show that Pakistani cases are hyperlipidemic for lipid traits except for HDL which is lowered. Patients with comorbidities also had lipid profiles deviated from the normal range. CONCLUSION: The study provides information regarding the aberration of lipid profile in the metabolic disorders that can increase the predisposition to complications.

The combined role of allelic variants of IRS-1 and IRS-2 genes in susceptibility to type2 diabetes in the Punjabi Pakistani subjects.

BACKGROUND: Diabetes mellitus is a multifactorial disorder characterized by a high level of glucose in the blood. Both genetic and environmental factors interact to cause diabetes. Insulin receptor substrate (IRS) proteins have a significant part in insulin signaling pathways. We aimed to investigate the relationship of type 2 diabetes with a Gly972Arg (G972R) variant of the IRS-1 gene and Gly1057Asp (G1057D) polymorphism of IRS-2 gene in the population of Punjab, Pakistan. METHODS: We collected 926 samples, 500 healthy controls (fasting blood sugar < 99 mg/dL, random blood sugar < 126 mg/dL) and 426 cases with diabetes (fasting blood sugar > 99 mg/dL, random blood sugar > 126 mg/dL). Several anthropometric measurements were measured. Statistical analysis was performed by using SPSS to determine the allele group/genotype frequency of the selected variants in the study population. RESULTS: The genotyping results of G972R by RLFP-PCR showed the allelic frequency of G = 0.68 and R = 0.32 in controls while G = 0.71 and R = 0.29 in the cases. The minor R allele had a slightly higher frequency in the cases than the controls (OR = 0.86, CI 0.706-1.052, p = 0.17). The genotyping results of G1057D showed allelic frequency G = 0.74 and D = 0.26 in the controls while G = 0.961 and D = 0.29 in the cases. The minor D allele appeared to be a risk allele for this SNP although the difference in the allele frequencies was not statistically significant (OR = 1.55, CI 0.961-1.41, p = 0.108). The combined genotype analysis showed that the difference in the allele and genotype frequencies reached statistical difference between the cases and the controls as well as the odds ratio substantially increased when the R allele (G972R) was present together with D allele (G1057D) in any combination. When the association of single variants with the lipid traits was observed, only D allele (G1057D) showed significant association with TG, HDL and LDL, however when the analysis was repeated for combined genotypes using general linear model, many more significant associations between the genotype where D allele and R allele are together, were seen with many lipid traits. CONCLUSION: In conclusion, the single nucleotide polymorphisms with low-modest effect size may not affect the phenotype individually but when in combination, the effect becomes stronger and more visible, therefore, for the SNP association studies, the more the number of SNPs included in the analysis, the more meaningful the results.

Genetic risk score (GRS) constructed from polymorphisms in the PON1, IL-6, ITGB3, and ALDH2 genes is associated with the risk of coronary artery disease in Pakistani subjects.

BACKGROUND: Coronary artery disease (CAD) is a major killer in today's world. Pakistan is also affected by this non-communicable disease like other countries. It is a multifactorial disease and is influenced by many gene-gene and gene-environment interactions. METHODS: A total of 623 (219 controls, 404 cases) Pakistani subjects were genotyped for four SNPs, rs662 (PON1), rs5918 (ITGB3), rs671 (ALDH2), rs1800795 (IL-6) by PCR-RFLP. Various anthropometric parameters were noted and serum lipid profile was measured using commercially available kits. Statistical analysis was done by SPSS version 22. A Genetic Risk Score (GRS) was calculated from individual SNPs. The association of the SNPs and the GRS with CAD was checked using logistic regression. RESULTS: The results showed that the risk allele frequencies of all variants were higher in the cases than the controls, however the difference was not statistically significant association (p > 0.0125). The mean GRS in the controls was 3.99 ± 1.42 and in cases, it was 4.29 ± 1.39, the difference between the groups was significant (p = 0.0109). logistic regression of individual SNPs and GRS with the CAD showed that independent SNPs were not significantly associated with the CAD however, the GRS had a strong association (p = 1.4 × 10- 4). The subjects were divided into three groups based on GRS (Gp 1 with GRS 0-2, Gp 2 with GRS 3-5 and Gp 3 with GRS 6-8). The analysis of the effect of the individual SNPs and GRS groups on different lipid profile parameters revealed no significant association of any of the tested SNPs with any lipid parameter, however, the GRS groups showed marginally significant for TC and highly significant association for TG, LDL-c and HDL-c. CONCLUSION: In conclusion, use of a GRS can provide better information than individual SNPs. The larger the number of the SNPs included in the analysis, the better would be the risk prediction.