A genome-wide association study uncovers novel genomic regions and candidate genes of yield-related traits in upland cotton

Key message

A total of 62 SNPs associated with yield-related traits were identified by a GWAS. Based on significant SNPs, two candidate genes pleiotropically increase lint yield.

Abstract

Improved fibre yield is considered a constant goal of upland cotton (Gossypium hirsutum) breeding worldwide, but the understanding of the genetic basis controlling yield-related traits remains limited. To better decipher the molecular mechanism underlying these traits, we conducted a genome-wide association study to determine candidate loci associated with six yield-related traits in a population of 719 upland cotton germplasm accessions; to accomplish this, we used 10,511 single-nucleotide polymorphisms (SNPs) genotyped by an Illumina CottonSNP63K array. Six traits, including the boll number, boll weight, lint percentage, fruit branch number, seed index and lint index, were assessed in multiple environments; large variation in all phenotypes was detected across accessions. We identified 62 SNP loci that were significantly associated with different traits on chromosomes A07, D03, D05, D09, D10 and D12. A total of 689 candidate genes were screened, and 27 of them contained at least one significant SNP. Furthermore, two genes (Gh_D03G1064 and Gh_D12G2354) that pleiotropically increase lint yield were identified. These identified SNPs and candidate genes provide important insights into the genetic control underlying high yields in G. hirsutum, ultimately facilitating breeding programmes of high-yielding cotton.

     

A GWA study reveals genetic loci for body conformation traits in Chinese Laiwu pigs and its implications for human BMI

Pigs share numerous physiological and phenotypic similarities with human and thus have been considered as a good model in nonrodent mammals for the study of genetic basis of human obesity. Researches on candidate genes for obesity traits have successfully identified some common genes between humans and pigs. However, few studies have assessed how many similarities exist between the genetic architecture of obesity in pigs and humans by large-scale comparative genomics. Here, we performed a genome-wide association study (GWAS) using the porcine 60 K SNP Beadchip for BMI and other four conformation traits at three different ages in a Chinese Laiwu pig population, which shows a large variability in fat deposition. In total, 35 SNPs were found to be significant at Bonferroni-corrected 5 % chromosome-wise level (P = 2.13 × 10^?5) and 88 SNPs had suggestive (P < 10^?4) association with the conformation traits. Some SNPs showed age-dependent association. Intriguingly, out of 32 regions associated with BMI in pigs, 18 were homologous with the loci for BMI in humans. Furthermore, five closest genes to GWAS peaks including HIF1AN, SMYD3, COX10, SLMAP, and GBE1 have been already associated with BMI in humans, which makes them very promising candidates for these QTLs. The result of GO analysis provided strong support to the fact that mitochondria and synapse play important roles in obesity susceptibility, which is consistent with previous findings on human obesity, and it also implicated new gene sets related to chromatin modification and Ig-like C2-type 5 domain. Therefore, these results not only provide new insights into the genetic architecture of BMI in pigs but also highlight that humans and pigs share the significant overlap of obesity-related genes.     

Genome-wide association study on chicken carcass traits using sequence data imputed from SNP array

Chicken carcass traits are economically important for the chicken industry. Detecting which genes affect chicken carcass traits is of great benefit to the genetic improvement of this important agricultural species. To investigate the genetic mechanism of carcass traits in chickens, we carried out a genome-wide association study (GWAS). A total of 435 Chinese indigenous chickens were phenotyped for carcass weight (CW), eviscerated weight with giblets (EWG), and eviscerated weight (EW) after slaughter at 91 days and were genotyped using a 600-K single nucleotide polymorphism (SNP) genotyping array. Twenty-four birds were selected for sequencing, and the 600 K SNP panel data were imputed to sequence data with the 24 birds as the reference. Univariate GWASs were performed with GEMMA software using the whole genome sequence data imputed from SNP chip data. Finally, 3, 25, and 63 suggestively significant SNPs were identified to be associated with carcass weight (CW), eviscerated weight with giblets (EWG), and eviscerated weight (EW), respectively. Six candidate genes, RNF219, SCEL, MYCBP2, ETS1, APLP2, and PRDM10 were detected. SCEL and MYCBP2 were potentially associated with these three traits, RNF219 and APLP2 were potentially associated with EWG and EW, and ETS1 and PRDM10 were only potentially associated with EWG and EW, respectively. Compared with forefathers’ research, 10 reported QTLs associated with CW were located within a 5-Mb distance near the SNPs with P value lower than 1×10^?5. This study enriched the knowledge of the genetic mechanisms of chicken carcass traits.     

Genomewide association study for economic traits in the large yellow croaker with different numbers of extreme phenotypes

A traditional genomewide association study (GWAS) detects genotype–phenotype associations by the vast number of genotyped individuals. This method requires large-scale samples and considerable sequencing costs. Extreme phenotypic sampling proposes make GWAS more cost-efficient and are applied more widely. With extreme phenotypic sampling, we performed a GWAS for n-3 highly unsaturated fatty acids (HUFA) and eviscerated weight (EW) traits in the large yellow croaker population. Of the 32,249 and 29,748 detected SNPs for the two traits, three candidate regions were found in each trait. Three candidate regions associated with HUFA were known near genes on chromosomes 4 and 11, and three candidate regions were on chromosome 6, and 15 for the EW trait. By combing through our GWAS results and the biological functional analysis of the genes, we suggest that the FABP, DGAT, ATP8B1, FAF2 and CERS2 genes,  as well as the IGF2, BORA, CYP1A1, GRTP1 and HOX genes are promising candidate genes for n-3 HUFA and EW, respectively, in the large yellow croaker. Moreover, compared with the different numbers of the extreme phenotypic sampling, we conclude that 60% of the extreme phenotypic subsample can obtain a similar result as GWAS with whole phenotypes. Thus, extreme phenotypic sampling could save 40% of the cost for genotyping and DNA extraction without loss of the candidate regions and functional genes. Our study may provide a basis for further genomic breeding and a reference for others who want to perform GWAS with extreme phenotypes.     

SNP markers associated with body size and pelt length in American mink (<Emphasis Type="Italic">Neovison vison</Emphasis>)

Background

Identification of genes underlying production traits is a key aim of the mink research community. Recent availability of genomic tools have opened the possibility for faster genetic progress in mink breeding. Availability of mink genome assembly allows genome-wide association studies in mink.

Results

In this study, we used genotyping-by-sequencing to obtain single nucleotide polymorphism (SNP) genotypes of 2496 mink. After multiple rounds of filtering, we retained 28,336 high quality SNPs and 2352 individuals for a genome-wide association study (GWAS). We performed the first GWAS for body weight, behavior, along with 10 traits related to fur quality in mink.

Conclusions

Combining association results with existing functional information of genes and mammalian phenotype databases, we proposed WWC3, MAP2K4, SLC7A1 and USP22 as candidate genes for body weight and pelt length in mink.

     

Identification of loci controlling forage yield and nutritive value in diploid alfalfa using GBS-GWAS

Key message

We attempted to identify genomic regions controlling forage yield and nutritive value in alfalfa. Several candidate genes and associated genetic markers were identified that could potentially be useful for alfalfa breeding to more efficiently develop improved cultivars.

Abstract

Alfalfa is one of the most widely cultivated forage legumes worldwide and improving alfalfa forage yield and nutritive value is a major global breeding goal. Genotyping-by-sequencing (GBS) provides cost-effective molecular marker genotyping for genome-wide association studies (GWAS). Using more than 15,000 genome-wide single nucleotide polymorphisms (SNP) identified from GBS, we conducted a GWAS to investigate forage yield and nutritive value-related traits. We have detected a number of associations for all the traits evaluated and a number of associations detected were located on the Medicago truncatula genome. The SNP in a coding region of a cell wall biosynthesis gene was associated with several cell wall-related traits, and we suggest that it may be the causative polymorphism. Two other SNPs residing in meristematic development and early growth genes were found to associate with the total biomass yield. None of the SNPs associated with regrowth after harvest or with spring regrowth were mapped to the M. truncatula genome, possibly reflecting the fact that M. truncatula is an annual species related to alfalfa that typically has limited ability to regrow. The alleles we identify with the major impact on forage yield and nutritive value can be rapidly incorporated into our breeding program.

     

Genome-wide association study of <Emphasis Type="Italic">Mycobacterium avium</Emphasis> subspecies <Emphasis Type="Italic">Paratuberculosis</Emphasis> infection in Chinese Holstein

Background

Paratuberculosis is a contagious, chronic and enteric disease in ruminants, which is caused by Mycobacterium avium subspecies paratuberculosis (MAP) infection, resulting in enormous economic losses worldwide. There is currently no effective cure for MAP infection or a vaccine, it is thus important to explore the genetic variants that contribute to host susceptibility to infection by MAP, which may provide a better understanding of the mechanisms of paratuberculosis and benefit animal genetic improvement. Herein we performed a genome-wide association study (GWAS) to identify genomic regions and candidate genes associated with susceptibility to MAP infection in dairy cattle.

Results

Using Illumina Bovine 50?K (54,609 SNPs) and GeneSeek HD (138,893 SNPs) chips, two analytical approaches were performed, GRAMMAR-GC and ROADTRIPS in 937 Chinese Holstein cows, among which individuals genotyped by the 50?K chip were imputed to HD SNPs with Beagle software. Consequently, 15 and 11 significant SNPs (P?<?5?×?10^??5) were identified with GRAMMAR-GC and ROADTDRIPS, respectively. A total of 10 functional genes were in proximity to (i.e., within 1?Mb) these SNPs, including IL4, IL5, IL13, IRF1, MyD88, PACSIN1, DEF6, TDP2, ZAP70 and CSF2. Functional enrichment analysis showed that these genes were involved in immune related pathways, such as interleukin, T cell receptor signaling pathways and inflammatory bowel disease (IBD), implying their potential associations with susceptibility to MAP infection. In addition, by examining the publicly available cattle QTLdb, a previous QTL for MAP was found to be overlapped with one of regions detected currently at 32.5?Mb on BTA23, where the TDP2 gene was anchored.

Conclusions

In conclusion, we identified 26 SNPs located on 15 chromosomes in the Chinese Holstein population using two GWAS strategies with high density SNPs. Integrated analysis of GWAS, biological functions and the reported QTL information helps to detect positional candidate genes and the identification of regions associated with susceptibility to MAP traits in dairy cattle.

     

Natural variation and genetic make-up of leaf blade area in spring barley

Key message

GWAS analysis for leaf blade area (LA) revealed intriguing genomic regions associated with putatively novel QTL and known plant stature-related phytohormone and sugar-related genes.

Abstract

Despite long-standing studies in the morpho-physiological characters of leaf blade area (LA) in cereal crops, advanced genetic studies to explore its natural variation are lacking. The importance of modifying LA in improving cereal grain yield and the genes controlling leaf traits have been well studied in rice but not in temperate cereals. To better understand the natural genetic variation of LA at four developmental stages, main culm LA was measured from 215 worldwide spring barleys including 92 photoperiod-sensitive accessions [PHOTOPERIOD RESPONSE LOCUS 1 (Ppd-H1)] and 123 accessions with reduced photoperiod sensitivity (ppd-H1) locus under controlled greenhouse conditions (long-day; 16/8 h; ~ 20/~ 16 °C day/night). The LA of Ppd-H1-carrying accessions was always smaller than in ppd-H1-carrying accessions. We found that nine SNPs from the Ppd-H1 gene were present in the collection of which marker 9 (M9; G/T in the CCT-domain) showed the most significant and consistent effect on LA at all studied developmental stages. Genome-wide association scans (GWAS) showed that the accessions carrying the ppd-H1 allele T/M9 (late heading) possessed more genetic variation in LA than the Ppd-H1 group carrying G/M9 (early heading). Several QTL with major effects on LA variation were found close to plant stature-related heading time, phytohormone- and sugar-related genes. The results provide evidence that natural variation of LA is an important source for improving grain yield, adaptation and canopy architecture of temperate cereals.

     

Genome-wide association study of body weight in Wenshang Barred chicken based on the SLAF-seq technology

Chicken body weight (BW) is an economically important trait, and many studies have been conducted on genetic selection for BW. However, previous studies have detected functional chromosome mutations or regions using gene chips. The present study used the specific-locus amplified fragment sequencing (SLAF-seq) technology to perform a genome-wide association study (GWAS) on purebred Wengshang Barred chicken. A total of 1,286,715 single-nucleotide polymorphisms (SNPs) were detected, and 175,211 SNPs were selected as candidate SNPs for genome-wide association analysis using TASSEL general linear models. Six SNP markers reached genome-wide significance. Of these, rs732048524, rs735522839, rs738991545, and rs15837818 were significantly associated with body weight at 28 days (BW28), while rs314086457 and rs315694878 were significantly associated with BW120. These SNPs are close to seven genes (PRSS23, ME3, FAM181B, NABP1, SDPR, TSSK6L2, and RBBP8). Moreover, 24 BW-associated SNPs reached “suggestive” genome-wide significance. Of these, 6, 13, 1, and 4 SNPs were associated with BW28, BW56, BW80, and BW120, respectively. These results would enrich the studies on BW and promote the use of Chinese chicken, especially the Wenshang Barred chicken.     

A genome-wide association study on growth traits in orange-spotted grouper (<Emphasis Type="Italic">Epinephelus coioides</Emphasis>) with RAD-seq genotyping

The orange-spotted grouper, Epinephelus coioides, is one of the most popular fish in China and Southeast Asian countries because of its important economic value. However, molecular mechanism underlying the growth of orange-spotted grouper has never been fully understood. Herein, we performed a genome-wide association study (GWAS) on a natural population of 198 individuals aiming to screen the whole genome of orange-spotted grouper for identification of growth-related loci by restrictionsite associated DNA sequencing. In this research, 261,366 single nucleotide polymorphisms (SNPs) were developed, in which 110 SNPs were identified to be correlated with growth and 20 SNPs were further confirmed to be associated with both body weight and total length. From these identified SNPs, we annotated a total of 34 genes, including adgrb2, csnkza1, cers5, col22a1, creb5, dnd1, dzank1, dnai1, npy2r, fat3, lrrk2, lrp5, map3k9, and so on. Among these candidate genes, npy2r (neuropeptide Y receptor Y2) was reported to play a critical role in growth of the orange-spotted grouper. In addition, population structure, principal component analysis, kinship matrix and linkage disequilibrium were examined to verify the accuracy and reliability of our GWAS results. Our data will also provide a valuable genetic resource for further marker-assisted selection program to improve growth quality in groupers.     

Identification and validation of genomic regions influencing kernel zinc and iron in maize

Key message

Genome-wide association study (GWAS) on 923 maize lines and validation in bi-parental populations identified significant genomic regions for kernel-Zinc and-Iron in maize.

Abstract

Bio-fortification of maize with elevated Zinc (Zn) and Iron (Fe) holds considerable promise for alleviating under-nutrition among the world’s poor. Bio-fortification through molecular breeding could be an economical strategy for developing nutritious maize, and hence in this study, we adopted GWAS to identify markers associated with high kernel-Zn and Fe in maize and subsequently validated marker-trait associations in independent bi-parental populations. For GWAS, we evaluated a diverse maize association mapping panel of 923 inbred lines across three environments and detected trait associations using high-density Single nucleotide polymorphism (SNPs) obtained through genotyping-by-sequencing. Phenotyping trials of the GWAS panel showed high heritability and moderate correlation between kernel-Zn and Fe concentrations. GWAS revealed a total of 46 SNPs (Zn-20 and Fe-26) significantly associated (P?≤?5.03?×?10^?05) with kernel-Zn and Fe concentrations with some of these associated SNPs located within previously reported QTL intervals for these traits. Three double-haploid (DH) populations were developed using lines identified from the panel that were contrasting for these micronutrients. The DH populations were phenotyped at two environments and were used for validating significant SNPs (P?≤?1?×?10^?03) based on single marker QTL analysis. Based on this analysis, 11 (Zn) and 11 (Fe) SNPs were found to have significant effect on the trait variance (P?≤?0.01, R²?≥?0.05) in at least one bi-parental population. These findings are being pursued in the kernel-Zn and Fe breeding program, and could hold great value in functional analysis and possible cloning of high-value genes for these traits in maize.

     

SNP-based pool genotyping and haplotype analysis accelerate fine-mapping of the wheat genomic region containing stripe rust resistance gene <Emphasis Type="Italic">Yr26</Emphasis>

Key message

NGS-assisted super pooling emerging as powerful tool to accelerate gene mapping and haplotype association analysis within target region uncovering specific linkage SNPs or alleles for marker-assisted gene pyramiding.

Abstract

Conventional gene mapping methods to identify genes associated with important agronomic traits require significant amounts of financial support and time. Here, a single nucleotide polymorphism (SNP)-based mapping approach, RNA-Seq and SNP array assisted super pooling analysis, was used for rapid mining of a candidate genomic region for stripe rust resistance gene Yr26 that has been widely used in wheat breeding programs in China. Large DNA and RNA super-pools were genotyped by Wheat SNP Array and sequenced by Illumina HiSeq, respectively. Hundreds of thousands of SNPs were identified and then filtered by multiple filtering criteria. Among selected SNPs, over 900 were found within an overlapping interval of less than 30 Mb as the Yr26 candidate genomic region in the centromeric region of chromosome arm 1BL. The 235 chromosome-specific SNPs were converted into KASP assays to validate the Yr26 interval in different genetic populations. Using a high-resolution mapping population (>?30,000 gametes), we confined Yr26 to a 0.003-cM interval. The Yr26 target region was anchored to the common wheat IWGSC RefSeq v1.0 and wild emmer WEWSeq v.1.0 sequences, from which 488 and 454 kb fragments were obtained. Several candidate genes were identified in the target genomic region, but there was no typical resistance gene in either genome region. Haplotype analysis identified specific SNPs linked to Yr26 and developed robust and breeder-friendly KASP markers. This integration strategy can be applied to accelerate generating many markers closely linked to target genes/QTL for a trait of interest in wheat and other polyploid species.

     

Targeted next-generation sequencing for identifying genes related to horse temperament

It is a fundamental challenge to discover the association of genetic traits with phenotypic traits. In this study, we aimed to identify possible genetic traits related to horse temperament. Based on previous findings, we selected 71 candidate genes related to temperamental trait and examined them in the human and horse reference genomes (hg38 and equCab2, respectively). We found 16 orthologous genes and, by comparing with the human reference genome, 17 homologous genes in the horse reference genome. We designed probes specific for the 33 horse genes. Using the probes, we built sequencing libraries of the genomic DNA samples from eight aggressive and eight docile horses, and sequenced the constructed libraries using the Illumina Hiseq2500 platform. Through the analysis of the targeted exome sequences, we identified single nucleotide polymorphisms (SNPs) in the genes. SNPs could be served as genetic markers to evaluate aggressive or docile levels of horses. To examine whether any genetic variants are associated with horse temperament, we performed genome-wide association study (GWAS) using the SNP data. GWAS analysis identified ten variants (p-value?<0.05) which could be related to horse temperament. We validated the variants using Sanger sequencing. The most significant variants were found in MAOA (c.1164+41T>C) and AR (c.1047+27G>T) genes with 8.09?×?10^?4 p-value. We suggest that the variants might be used to assess horse temperament and to determine superior horses for riding or racing.     

The effects of polymorphisms in <Emphasis Type="Italic">growth hormone</Emphasis> and <Emphasis Type="Italic">growth hormone receptor</Emphasis> genes on production and reproduction traits in Aberdeen-Angus cattle (<Emphasis Type="Italic">Bos taurus</Emphasis> L., 1758)

The study was aimed to analyze the relation between individual genotypes and allelic variants of SNPs g.2141C>G of growth hormone gene, g.914T>A and g.257A>G of growth hormone receptor gene with growth and reproduction traits and to evaluate the populationgenetic structure in Aberdeen-Angus cattle (Bos taurus L., 1758) sample of Eastern Ukraine according SNPs studied. Allele C of SNP g.2141C>G has a positive correlation with birth weight, body stature, bigger rump, udder and total exterior evaluation score, shorter calving interval and better calve birth weight and negative correlation with calve average daily gain. Allele T of SNP g.914T>A has positive correlation with the muscle and udder size; live weight in each age, average daily gain, weight and average daily gain of calves born conform to the principle AA>TT≥TA. SNP g.257A>G showed a positive correlation for G allele with muscle size. The population is in equilibrium for SNPs g.2141C>G and g.257A>G, and in disequilibrium for SNP g.914T>A. The analysis showed no linkage disequilibrium between SNPs g.914T>A and g.257A>G. Inbreeding coefficient F_ST in Aberdeen-Angus group studied was 16.1%.     

Single nucleotide polymorphism markers for rapid detection of the <Emphasis Type="Italic">Rsv4</Emphasis> locus for soybean mosaic virus resistance in diverse germplasm

Soybean mosaic virus (SMV) causes a substantial decrease in soybean yield and reduction of seed quality. The most effective management strategy to control the virus is the deployment of host resistance. Seven SMV strains and three independent multi-allelic loci for SMV resistance have been identified previously. The goal of this research was to detect single nucleotide polymorphisms (SNPs) associated with SMV resistance at the Rsv4 locus. Ten soybean accessions, with confirmed resistance genes, were used for sequencing the candidate gene Glyma.02g121400. Alignment of these sequences revealed three SNPs displaying 100% consistency for genotypes carrying the Rsv4 gene. These SNPs were applied for a rapid screen of diverse soybean germplasm using the Sequenom iPLEX Gold platform, phenotyped with SMV-G1 and G7 strains to determine phenotype and classified into several groups carrying the proposed R-gene. The population of V94-5152 (Rsv4) × Lee 68 (rsv) was screened using novel SNPs to create a genetic map with improved resolution to determine the location of the Rsv4. To observe the recombination frequencies within the population, three additional SNPs on both sides of the Glyma.02g121400 gene were added. A linkage map revealed a distance of 3.6 cM between the Rsv4 locus and the closest SNP, thus shifting the putative Rsv4 region downstream on chromosome 2. With this region, five candidate genes have been proposed. The genomic position of the discovered SNPs, linked to the Rsv4, could increase screening precision and accelerate breeding efforts to develop multi-strain-resistant crops.     

Single nucleotide polymorphisms in <Emphasis Type="Italic">TaER</Emphasis> genes and their association with carbon isotope discrimination in wheat genotypes under drought

Candidate gene association studies implicate the detection of contributing single nucleotide polymorphism (SNP) for the target traits and have been recommended as a promising technique to anatomize the complex characters in plants. The ERECTA gene in plants controls different physiological functions. In this study, we identified SNPs in 1.1 kb partial sequences of TaER-1 and TaER-2 of wheat (Triticum aestivum L.). Thirty-nine SNPs were identified in the coding regions of TaER-1 gene in 33 wheat genotypes, of which 20 SNPs caused non-synonymous mutations while 19 SNPs produced synonymous mutations; 31 SNPs were located in the coding regions of TaER-2 gene in 26 genotypes, of which 18 SNPs caused non-synonymous mutations and 13 SNPs caused synonymous mutations. In addition, 32 SNPs in TaER-1 and 9 SNPs in TaER-2 were also identified in the non-coding regions. Moreover, the significant genetic associations of SNPs of TaER-1 and TaER-2 genes with carbon isotope discrimination, stomatal conductance, photosynthetic rate, transpiration rate, intrinsic water use efficiency (iWUE), leaf length, leaf width, stomatal density, epidermal cell density, and stomatal index were noted in wheat genotypes. This study confirms the importance of TaER-1 and TaER-2 genes which could improve iWUE of wheat by regulating leaf gas exchange and leaf structural traits. These identified SNPs may play a critical role in molecular breeding by means of marker-assisted selection.     

Association of the <Emphasis Type="Italic">rs1870634</Emphasis> Variant in Long Intergenic Non-protein Coding RNA 841 with Coronary Artery Disease: A GWAS-Replication Study in an Iranian Population

Recent genome-wide association studies (GWAS) identified a list of single-nucleotide polymorphisms (SNPs) associated with coronary artery disease (CAD). Replication of GWAS findings in different population corroborated the observed association in the parent GWAS. In this study, we aimed to replicate the association of rs1870634, a GWAS identified SNP, to CAD in an Iranian population. The study population consisted of 267 subjects undergoing coronary angiography coronary angiography including 155 CAD patients and 112 non-CAD age- and gender-matched controls. The genotype determination of rs1870634 SNP performed using high-resolution melting analysis (HRM) technique. Our results revealed that the GG genotype frequency was significantly higher in CAD patients compared with controls (P?=?0.03). The results of binary logistic regression suggested that this genotype was significantly associated with CAD risk adjustment for age, BMI, sex, TC, and LDL-C lipid levels (OR of 2.78, 95% CI (1.10–7.01), P?=?0.03). Moreover, our results showed that the GG+TG genotypes were 2.52 times more likely to develop CAD (95% CI 1.05–6.03) than TT genotype carriers after adjusting for age, sex, and lipid profiles (P?=?0.037). These data showed that the GG genotype could be associated with increased risk of CAD in a sample of Iranian population.