Genetics of venous thrombosis: insights from a new genome wide association study

Background

Venous Thrombosis (VT) is a common multifactorial disease associated with a major public health burden. Genetics factors are known to contribute to the susceptibility of the disease but how many genes are involved and their contribution to VT risk still remain obscure. We aimed to identify genetic variants associated with VT risk.

Methodology/Principal Findings

We conducted a genome-wide association study (GWAS) based on 551,141 SNPs genotyped in 1,542 cases and 1,110 controls. Twelve SNPs reached the genome-wide significance level of 2.0×10⁻⁸ and encompassed four known VT-associated loci, ABO, F5, F11 and FGG. By means of haplotype analyses, we also provided novel arguments in favor of a role of HIVEP1, PROCR and STAB2, three loci recently hypothesized to participate in the susceptibility to VT. However, no novel VT-associated loci came out of our GWAS. Using a recently proposed statistical methodology, we also showed that common variants could explain about 35% of the genetic variance underlying VT susceptibility among which 3% could be attributable to the main identified VT loci. This analysis additionally suggested that the common variants left to be identified are not uniformly distributed across the genome and that chromosome 20, itself, could contribute to ∼7% of the total genetic variance.

Conclusions/Significance

This study might also provide a valuable source of information to expand our understanding of biological mechanisms regulating quantitative biomarkers for VT.     

Investigation of transferability of BovineSNP50 BeadChip from cattle to water buffalo for genome wide association study

Cattle and water buffalo belong to the same subfamily Bovinae and share chromosome banding and gene order homology. In this study, we used genome-wide Illumina BovineSNP50 BeadChip to analyze 91 DNA samples from three breeds of water buffalo (Nili-Ravi, Murrah and their crossbred with local GuangXi buffalos in China), to demonstrate the genetic divergence between cattle and water buffalo through a large single nucleotide polymorphism (SNP) transferability study at the whole genome level, and performed association analysis of functional traits in water buffalo as well. A total of 40,766 (75.5 %) bovine SNPs were found in the water buffalo genome, but 49,936 (92.5 %) were with only one allele, and finally 935 were identified to be polymorphic and useful for association analysis in water buffalo. Therefore, the genome sequences of water buffalo and cattle shared a high level of homology but the polymorphic status of the bovine SNPs varied between these two species. The different patterns of mutations between species may associate with their phenotypic divergence due to genome evolution. Among 935 bovine SNPs, we identified a total of 9 and 7 SNPs significantly associated to fertility and milk production traits in water buffalo, respectively. However, more works in larger sample size are needed in future to verify these candidate SNPs for water buffalo.     

Large intraspecific haplotype variability at the Rph7 locus results from rapid and recent divergence in the barley genome

To study genome evolution and diversity in barley (Hordeum vulgare), we have sequenced and compared more than 300 kb of sequence spanning the Rph7 leaf rust disease resistance gene in two barley cultivars. Colinearity was restricted to five genic and two intergenic regions representing <35% of the two sequences. In each interval separating the seven conserved regions, the number and type of repetitive elements were completely different between the two homologous sequences, and a single gene was absent in one cultivar. In both cultivars, the nonconserved regions consisted of approximately 53% repetitive sequences mainly represented by long-terminal repeat retrotransposons that have inserted <1 million years ago. PCR-based analysis of intergenic regions at the Rph7 locus and at three other independent loci in 41 H. vulgare lines indicated large haplotype variability in the cultivated barley gene pool. Together, our data indicate rapid and recent divergence at homologous loci in the genome of H. vulgare, possibly providing the molecular mechanism for the generation of high diversity in the barley gene pool. Finally, comparative analysis of the gene composition in barley, wheat (Triticum aestivum), rice (Oryza sativa), and sorghum (Sorghum bicolor) suggested massive gene movements at the Rph7 locus in the Triticeae lineage.     

A common haplotype of the glucokinase gene alters fasting glucose and birth weight: association in six studies and population-genetics analyses

Fasting glucose is associated with future risk of type 2 diabetes and ischemic heart disease and is tightly regulated despite considerable variation in quantity, type, and timing of food intake. In pregnancy, maternal fasting glucose concentration is an important determinant of offspring birth weight. The key determinant of fasting glucose is the enzyme glucokinase (GCK). Rare mutations of GCK cause fasting hyperglycemia and alter birth weight. The extent to which common variation of GCK explains normal variation of fasting glucose and birth weight is not known. We aimed to comprehensively define the role of variation of GCK in determination of fasting glucose and birth weight, using a tagging SNP (tSNP) approach and studying 19,806 subjects from six population-based studies. Using 22 tSNPs, we showed that the variant rs1799884 is associated with fasting glucose at all ages in the normal population and exceeded genomewide levels of significance (P=10-9). rs3757840 was also highly significantly associated with fasting glucose (P=8x10-7), but haplotype analysis revealed that this is explained by linkage disequilibrium (r2=0.2) with rs1799884. A maternal A allele at rs1799884 was associated with a 32-g (95% confidence interval 11-53 g) increase in offspring birth weight (P=.002). Genetic variation influencing birth weight may have conferred a selective advantage in human populations. We performed extensive population-genetics analyses to look for evidence of recent positive natural selection on patterns of GCK variation. However, we found no strong signature of positive selection. In conclusion, a comprehensive analysis of common variation of the glucokinase gene shows that this is the first gene to be reproducibly associated with fasting glucose and fetal growth.     

SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci

The authors want to changed the web link of the software platform in this Briefing.Page 97,section "INTRODUCTION",the web link of SHEsis is changed from http://www.nhgg.org/analysis tohttp://analysis.bio-x.cn     

SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci

In multiloci-based genetic association studies of complex diseases, a powerful and high efficient tool for analyses of linkage disequilibrium (LD) between markers, haplotype distributions and many chi-square/p values with a large number of samples has been sought for long. In order to achieve the goal of obtaining meaningful results directly from raw data,we developed a robust and user-friendly software platform with a series of tools for analysis in association study with high efficiency. The platform has been well evaluated by several sets of real data.     

The genetic architecture of selection at the human dopamine receptor D4 (DRD4) gene locus

Associations of the seven-repeat (7R) allele of the human dopamine receptor D4 (DRD4) gene with both the personality trait of novelty seeking and attention deficit/hyperactivity disorder have been reported. Recently, on the basis of the unusual DNA sequence organization of the DRD4 7R 48-bp tandem repeat (VNTR), we proposed that the 7R allele originated as a rare mutational event that increased to high frequency by positive selection. We now have resequenced the entire DRD4 locus from 103 individuals homozygous for 2R, 4R, or 7R variants of the VNTR, a method developed to directly estimate haplotype diversity. DNA from individuals of African, European, Asian, North and South American, and Pacific Island ancestry were used. 4R/4R homozygotes exhibit little linkage disequilibrium (LD) over the region examined, with more polymorphisms observed in DNA samples from African individuals. In contrast, the evidence for strong LD surrounding the 7R allele is dramatic, with all 7R/7R individuals (including those from Africa) exhibiting the same alleles at most polymorphic sites. By intra-allelic comparison at 18 high-heterozygosity sites spanning the locus, we estimate that the 7R allele arose prior to the upper Paleolithic era (approximately 40000-50000 years ago). Further, the pattern of recombination at these polymorphic sites is the pattern expected for selection acting at the 7R VNTR itself, rather than at an adjacent site. We propose a model for selection at the DRD4 locus consistent with these observed LD patterns and with the known biochemical and physiological differences between receptor variants.     

Dimerization of the melanocortin 4 receptor: a study using bioluminescence resonance energy transfer

The melanocortin 4 receptor is important in the regulation of satiety. In this study we have investigated the propensity of the MC4 receptor to homodimerize. MC4 receptors with either a modified green fluorescent protein (GFP(2)) or Renilla luciferase (RLuc) at their C-terminus were constructed. These receptors showed equivalent binding and functional properties to the wild-type MC4 receptor. Bioluminescence resonance energy transfer readings indicated that the MC4 receptor exists as a constitutive homodimer, which was not regulated by peptide interaction. The efficiency of MC4 receptor to form homodimers was greatly enhanced compared to its ability to heterodimerize with the kappa opioid receptor.     

Identification of genetic variation and haplotype structure of the canine ABCA4 gene for retinal disease association studies

Over 200 mutations in the retina specific member of the ATP-binding cassette transporter superfamily (ABCA4) have been associated with a diverse group of human retinal diseases. The disease mechanisms, and genotype–phenotype associations, nonetheless, remain elusive in many cases. As orthologous genes are commonly mutated in canine models of human blinding disorders, canine ABCA4 appears to be an ideal candidate gene to identify and study sequence changes in dogs affected by various forms of inherited retinal degeneration. However, the size of the gene and lack of haplotype assignment significantly limit targeted association and/or linkage approaches. This study assessed the naturally observed sequence diversity of ABCA4 in the dog, identifying 80% of novel variations. While none of the observed polymorphisms have been associated with blinding disorders to date, breed and potentially disease specific haplotypes have been identified. Moreover, a tag SNP map of 17 (15) markers has been established that accurately predicts common ABCA4 haplotypes (frequency > 5%) explaining >85% (>80%) of the observed genetic diversity and will considerably advance future studies. Our sequence analysis of the complete canine ABCA4 coding region will clearly provide a baseline and tools for future association studies and comparative genomics to further delineate the role of ABCA4 in canine blinding disorders.     

Marked association of a RFLP for the low density lipoprotein receptor gene with obesity in essential hypertensives.

RFLPs at the low density lipoprotein receptor locus (LDLR) display marked linkage disequilibrium between each other. Cross-sectional analysis of a bi-alleleic ApaLI RFLP of LDLR showed that the 9.4- and 6.6-kb alleles were present in similar frequency between a group of 84 Caucasian essential hypertensive (HT) and a group of 96 normotensive subjects whose parents each had a similar blood pressure status at age > or = 50. After subdividing HTs into lean and obese, however, the frequency of the 6.6-kb allele in the 27 HTs with BMI > or = 26 kg/m2 was 0.63, compared with 0.39 for HTs with BMI < 26 (chi 2 = 8.8; P = 0.004). The difference in genotype frequencies was even more striking (chi 2 = 23; P = 0.00008), with a virtual absence of 9.4-kb homozygotes in the obese HT group (1 vs 22). Genetic variation at LDLR (19p13.2) is thus associated with obesity in HT.     

1H-NMR studies on a potent and selective antagonist at human melanocortin receptor 4 (hMC-4R)

Melanocortin receptor 4 (MC-4R) is involved in the regulation of energy balance and body weight, and recognizes alpha-, beta-, and gamma-melanocyte stimulating hormones (alpha-, beta-, gamma-MSH). In the search for compounds that regulate food intake and body weight, two synthetic lactam-derivative ligands of alpha-MSH were discovered, MTII and SHU9119. MTII is an agonist and reduces food intake in rats, whereas SHU9119 is an antagonist, and increases food intake and body weight in rats. MTII and SHU9119 are nonselective compounds to MC-4R. To enhance the potency and selectivity at the human MC-4R (hMC-4R), we recently synthesized analogs of SHU9119 (M. A. Bednarek, T. MacNeil, R. N. Kalyani, R. Tang, Van der L. H. T. Ploeg, and D. H. Weinberg, Journal of Medicinal Chemistry, 2001, Vol. 44, pp. 401-409), wherein compound 1 was the most selective for hMC-4R. Replacement of D-Nal by L-Nal in compound 1 made compound 2 weakly active. Comparison of the structures by NMR and molecular modeling of compounds 1 and 2 vs SHU9119 and MTII indicated that, even though they existed as an average of several conformations in solution, there were distinctions in their structures. The gamma-methylene protons of Arg in compound 1 were nonequivalent and shielded probably by the aromatic ring of Nal. The NHi-NHi+1 NOE cross peaks and the temperature coefficients of the amide protons around the "essential core" Nal/Phe7-Arg8-Trp9, required for high affinity and high selectivity at hMC-4R, were indicative of a possible turn structure for these compounds but with differences in their NOE strengths and temperature coefficient values. Molecular modeling of these compounds based on their NMR data showed that the essential core appeared as a "V" shape with two different orientations, one for compound 1 and some of the conformers of SHU9119 and MTII, and the other for compound 2 and some other conformers of SHU9119 and MTII. The remaining conformers of SHU9119 and MTII, which did not map to compound 1 or 2, suggested that they were outside of the hMC-4R binding envelop. These observations may lead to conjectures as to why compound 1 is highly active and selective toward hMC-4R.     

Structure-activity relationship and signal transduction of gamma-MSH peptides in GH3 cells: further evidence for a new melanocortin receptor

The structure-activity relationship and signal transduction properties of the pro-opiomelanocortin (POMC)-derived gamma-MSH peptides in the GH3 cell line was compared with that described for the known melanocortin receptors (MCRs). Single alanine replacements showed that, unlike the classical MCRs, the His(5)-Phe(6)-Arg(7)-Trp(8) sequence in gamma2-MSH is not a core sequence for activating the gamma-MSH receptor in GH3 cells, whereas Met(3) is essential. gamma2-MSH increased binding of [35S]GTPgammaS to membrane preparations of GH3 cells. Blockade of protein kinase A abolished the [Ca(2+)](i) responses to gamma3-MSH, and low nanomolar doses of gamma3-MSH increased intracellular cAMP levels, which could be blocked by pertussis toxin (PTX). We conclude that the putative novel gamma-MSH receptor in GH3 cells is a GPCR, but with structure-activity and signal transduction features different from those of the classical MCRs.