Genome-wide association study of migraine implicates a common susceptibility variant on 8q22.1

Migraine is a common episodic neurological disorder, typically presenting with recurrent attacks of severe headache and autonomic dysfunction. Apart from rare monogenic subtypes, no genetic or molecular markers for migraine have been convincingly established. We identified the minor allele of rs1835740 on chromosome 8q22.1 to be associated with migraine (P = 5.38 x 10(-)(9), odds ratio = 1.23, 95% CI 1.150-1.324) in a genome-wide association study of 2,731 migraine cases ascertained from three European headache clinics and 10,747 population-matched controls. The association was replicated in 3,202 cases and 40,062 controls for an overall meta-analysis P value of 1.69 x 10(-)(1)(1) (odds ratio = 1.18, 95% CI 1.127-1.244). rs1835740 is located between MTDH (astrocyte elevated gene 1, also known as AEG-1) and PGCP (encoding plasma glutamate carboxypeptidase). In an expression quantitative trait study in lymphoblastoid cell lines, transcript levels of the MTDH were found to have a significant correlation to rs1835740 (P = 3.96 x 10(-)(5), permuted threshold for genome-wide significance 7.7 x 10(-)(5). To our knowledge, our data establish rs1835740 as the first genetic risk factor for migraine.

Sequence variants in three loci influence monocyte counts and erythrocyte volume

Blood cells participate in vital physiological processes, and their numbers are tightly regulated so that homeostasis is maintained. Disruption of key regulatory mechanisms underlies many blood-related Mendelian diseases but also contributes to more common disorders, including atherosclerosis. We searched for quantitative trait loci (QTL) for hematology traits through a whole-genome association study, because these could provide new insights into both hemopoeitic and disease mechanisms. We tested 1.8 million variants for association with 13 hematology traits measured in 6015 individuals from the Australian and Dutch populations. These traits included hemoglobin composition, platelet counts, and red blood cell and white blood cell indices. We identified three regions of strong association that, to our knowledge, have not been previously reported in the literature. The first was located in an intergenic region of chromosome 9q31 near LPAR1, explaining 1.5% of the variation in monocyte counts (best SNP rs7023923, p=8.9x10(-14)). The second locus was located on chromosome 6p21 and associated with mean cell erythrocyte volume (rs12661667, p=1.2x10(-9), 0.7% variance explained) in a region that spanned five genes, including CCND3, a member of the D-cyclin gene family that is involved in hematopoietic stem cell expansion. The third region was also associated with erythrocyte volume and was located in an intergenic region on chromosome 6q24 (rs592423, p=5.3x10(-9), 0.6% variance explained). All three loci replicated in an independent panel of 1543 individuals (p values=0.001, 9.9x10(-5), and 7x10(-5), respectively). The identification of these QTL provides new opportunities for furthering our understanding of the mechanisms regulating hemopoietic cell fate.

Replication of the association of common rs9939609 variant of FTO with increased BMI in an Australian adult twin population but no evidence for gene by environment (G x E) interaction

OBJECTIVE: To further investigate a common variant (rs9939609) in the fat mass- and obesity-associated gene (FTO), which recent genome-wide association studies have shown to be associated with body mass index (BMI) and obesity. DESIGN: We examined the effect of this FTO variant on BMI in 3353 Australian adult male and female twins. RESULTS: The minor A allele of rs9939609 was associated with an increased BMI (P=0.0007). Each additional copy of the A allele was associated with a mean BMI increase of approximately 1.04 kg/m(2) (approximately 3.71 kg). Using variance components decomposition, we estimate that this single-nucleotide polymorphism accounts for approximately 3% of the genetic variance in BMI in our sample (approximately 2% of the total variance). By comparing intrapair variances of monozygotic twins of different genotypes we were able to perform a direct test of gene by environment (G x E) interaction in both sexes and gene by parity (G x P) interaction in women, but no evidence was found for either. CONCLUSIONS: In addition to supporting earlier findings that the rs9939609 variant in the FTO gene is associated with an increased BMI, our results indicate that the associated genetic effect does not interact with environment or parity.

Genome-wide association studies of quantitative traits with related individuals: Little (power) lost but much to be gained

For complex disease genetics research in human populations, remarkable progress has been made in recent times with the publication of a number of genome-wide association scans (GWAS) and subsequent statistical replications. These studies have identified new genes and pathways implicated in disease, many of which were not known before. Given these early successes, more GWAS are being conducted and planned, both for disease and quantitative phenotypes. Many researchers and clinicians have DNA samples available on collections of families, including both cases and controls. Twin registries around the world have facilitated the collection of large numbers of families, with DNA and multiple quantitative phenotypes collected on twin pairs and their relatives. In the design of a new GWAS with a fixed budget for the number of chips, the question arises whether to include or exclude related individuals. It is commonly believed to be preferable to use unrelated individuals in the first stage of a GWAS because relatives are 'over-matched' for genotypes. In this study, we quantify that for GWAS of a quantitative phenotype, relative to a sample of unrelated individuals surprisingly little power is lost when using relatives. The advantages of using relatives are manifold, including the ability to perform more quality control, the choice to perform within-family tests of association that are robust to population stratification, and the ability to perform joint linkage and association analysis. Therefore, the advantages of using relatives in GWAS for quantitative traits may well outweigh the small disadvantage in terms of statistical power.

A high-density association screen of 155 ion transport genes for involvement with common migraine

The clinical overlap between monogenic Familial Hemiplegic Migraine (FHM) and common migraine subtypes, and the fact that all three FHM genes are involved in the transport of ions, suggest that ion transport genes may underlie susceptibility to common forms of migraine. To test this leading hypothesis, we examined common variation in 155 ion transport genes using 5257 single nucleotide polymorphisms (SNPs) in a Finnish sample of 841 unrelated migraine with aura cases and 884 unrelated non-migraine controls. The top signals were then tested for replication in four independent migraine case-control samples from the Netherlands, Germany and Australia, totalling 2835 unrelated migraine cases and 2740 unrelated controls. SNPs within 12 genes (KCNB2, KCNQ3, CLIC5, ATP2C2, CACNA1E, CACNB2, KCNE2, KCNK12, KCNK2, KCNS3, SCN5A and SCN9A) with promising nominal association (0.00041 < P < 0.005) in the Finnish sample were selected for replication. Although no variant remained significant after adjusting for multiple testing nor produced consistent evidence for association across all cohorts, a significant epistatic interaction between KCNB2 SNP rs1431656 (chromosome 8q13.3) and CACNB2 SNP rs7076100 (chromosome 10p12.33) (pointwise P = 0.00002; global P = 0.02) was observed in the Finnish case-control sample. We conclude that common variants of moderate effect size in ion transport genes do not play a major role in susceptibility to common migraine within these European populations, although there is some evidence for epistatic interaction between potassium and calcium channel genes, KCNB2 and CACNB2. Multiple rare variants or trans-regulatory elements of these genes are not ruled out.

ssSNPer: identifying statistically similar SNPs to aid interpretation of genetic association studies

ssSNPer is a novel user-friendly web interface that provides easy determination of the number and location of untested HapMap SNPs, in the region surrounding a tested HapMap SNP, which are statistically similar and would thus produce comparable and perhaps more significant association results. Identification of ssSNPs can have crucial implications for the interpretation of the initial association results and the design of follow-up studies. AVAILABILITY: http://fraser.qimr.edu.au/general/daleN/ssSNPer/

Association study of the dystrobrevin-binding gene with schizophrenia in Australian and Indian samples

Numerous studies have reported association between variants in the dystrobrevin binding protein 1 (dysbindin) gene (DTNBP1) and schizophrenia. However, the pattern of results is complex and to date, no specific risk marker or haplotype has been consistently identified. The number of single nucleotide polymorphisms (SNPs) tested in these studies has ranged from 5 to 20. We attempted to replicate previous findings by testing 16 SNPs in samples of 41 Australian pedigrees, 194 Australian cases and 180 controls, and 197 Indian pedigrees. No globally significant evidence for association was observed in any sample, despite power calculations indicating sufficient power to replicate several previous findings. Possible explanations for our results include sample differences in background linkage disequilibrium and/or risk allele effect size, the presence of multiple risk alleles upon different haplotypes, or the presence of a single risk allele upon multiple haplotypes. Some previous associations may also represent false positives. Examination of Caucasian HapMap phase II genotype data spanning the DTNBP1 region indicates upwards of 40 SNPs are required to satisfactorily assess all nonredundant variation within DTNBP1 and its potential regulatory regions for association with schizophrenia. More comprehensive studies in multiple samples will be required to determine whether specific DTNBP1 variants function as risk factors for schizophrenia.

Exclusion of angiotensinogen gene in molecular basis of human hypertension: Sibpair linkage and association analyses in Australian anglo-caucasians

Linkage with essential hypertension has been claimed for a microsatellite marker near the angiotensinogen gene (AGT; chromosome 1q42), as has association for the AGT variants M235T, G(-6)A and A(-20)C. To more rigorously evaluate AGT as a candidate gene for hypertension we performed sibpair analysis with multiple microsatellite markers surrounding this locus and using more sophisticated analysis programs. We also performed an association study of the AGT variants in unrelated subjects with a strong family history (two affected parents). For the linkage study, single and multiplex polymerase chain reaction (PCRs) and automated genescan analysis were conducted on DNA from 175 Australian Anglo-Celtic Caucasian hypertensives for the following markers: D1S2880-(2.1 cM)-D1S213-(2.8 cM)-D1S251-(6.5 cM)-AGT-(2.0 cM) -D1S235. Statistical evaluation of genotype data by nonparametric methods resulted in the following scores: Single-point analysis - SPLINK, P > 0.18; APM method, P > 0.25; ASPEX, MLOD < 0.28; SIB-PAIR, P > 0. 24; Multipoint analysis - MAPMAKER/SIBS, MLOD < 0.24; GENEHUNTER, P > 0.35. Exclusion scores of Lod -4.1 to -5.1 were obtained for these markers using MAPMAKER/SIBS for a lambda(s) of 1.6. The association study of G(-6)A, A(-20)C and M235T variants in 111 hypertensives with strong family history and 190 normotensives with no family history showed significant linkage disequilibrium between particular haplotypes, but we could find no association with hypertension. The present study therefore excludes AGT in the etiology of hypertension, at least in the population of Australian Anglo-Celtic Caucasians studied.

Association of a low density lipoprotein receptor microsatellite variant with obesity

OBJECTIVE: To determine whether a microsatellite polymorphism located towards the 3' end of the low density lipoprotein receptor gene (LDLR) is associated with obesity. DESIGN: A cross-sectional case-control study. SUBJECTS: One hundred and seven obese individuals, defined as a body mass index (BMI) > or = 26 kg/m2, and 163 lean individuals, defined as a BMI < 26 kg/m2. MEASUREMENTS: BMI, blood pressure, serum lipids, alleles of LDLR microsatellite (106 bp, 108 bp and 112 bp). RESULTS: There was a significant association between variants of the LDLR microsatellite and obesity, in the overall tested population, due to a contributing effect in females (chi 2 = 12.3, P = 0.002), but not in males (chi 2 = 0.3, P = 0.87). In females, individuals with the 106 bp allele were more likely to be lean, while individuals with the 112 bp and/or 108 bp alleles tended to be obese. CONCLUSIONS: These results suggest that in females, LDLR may play a role in the development of obesity.

Lens shape and refractive index distribution in type 1 diabetes

Purpose: To compare lens dimensions and refractive index distributions in type 1 diabetes and age-matched control groups. Methods: There were 17 participants with type 1 diabetes, consisting of two subgroups (7 young [23 ± 4 years] and 10 older [54 ± 4 years] participants), with 23 controls (13 young, 24 ± 4 years; 10 older, 55 ± 4 years). For each participant, one eye was tested with relaxed accommodation. A 3T clinical magnetic resonance imaging scanner was used to image the eye, employing a multiple spin echo (MSE) sequence to determine lens dimensions and refractive index profiles along the equatorial and axial directions. Results: The diabetes group had significantly smaller lens equatorial diameters and larger lens axial thicknesses than the control group (diameter mean ± 95% confidence interval [CI]: diabetes group 8.65 ± 0.26 mm, control group 9.42 ± 0.18 mm; axial thickness: diabetes group 4.33 ± 0.30 mm, control group 3.80 ± 0.14 mm). These differences were also significant within each age group. The older group had significantly greater axial thickness than the young group (older group 4.35 ± 0.26 mm, young group 3.70 ± 0.25 mm). Center refractive indices of diabetes and control groups were not significantly different. There were some statistically significant differences between the refractive index fitting parameters of young and older groups, but not between diabetes and control groups of the same age. Conclusions: Smaller lens diameters occurred in the diabetes groups than in the age-matched control groups. Differences in refractive index distribution between persons with and without diabetes are too small to have important effects on instruments measuring axial thickness.

Assessing the combined effect of asbestos exposure and smoking on lung cancer: a Bayesian approach

We review the literature on the combined association between lung cancer and two environmental exposures, asbestos exposure and smoking, and explore a Bayesian approach to assess evidence of interaction between the exposures. The meta-analysis combines separate indices of additive and multiplicative relationships and multivariate relative risk estimates. By making inferences on posterior probabilities we can explore both the form and strength of interaction. This analysis may be more informative than providing evidence to support one relation over another on the basis of statistical significance. Overall, we find evidence for a more than additive and less than multiplicative relation.

Electrophoretic mobility of cell nuclei (EMN index) as a biomarker of the biological aging process: Considering the association between EMN index and age

The present study examined whether a specific property of cell microstructures may be useful as a biomarker of aging. Specifically, the association between age and changes of cellular structures reflected in electrophoretic mobility of cell nuclei index (EMN index) values across the adult lifespan was examined. This report considers findings from cross sections of females (n = 1273) aged 18–98 years, and males (n = 506) aged 19–93 years. A Biotest apparatus was used to perform intracellular microelectrophoresis on buccal epithelial cells collected from each individual. EMN index was calculated on the basis of the number of epithelial cells with mobile nuclei in reference to the cells with immobile nuclei per 100 cells. Regression analyses indicated a significant negative association between EMN index value and age for men (r = −0.71, p < 0.001) and women (r = −0.60, p < 0.001); demonstrating a key requirement that must be met by a biomarker of aging. The strength of association observed between EMN index and age for both men and women was encouraging and supports the potential use of EMN index for determining a biological age of an individual (or a group). In this study, a new attempt of complex explanation of cellular mechanisms contributing to age related changes of the EMN index was made. In this study, a new attempt of complex explanation of cellular mechanisms contributing to age related changes of the EMN index was made. EMN index has demonstrated potential to meet criteria proposed for biomarkers of aging and further investigations are necessary.

The association of chronic kidney disease and dialysis treatment with foot ulceration and major amputation

Objective The objective of this study was to investigate the risk of chronic kidney disease (CKD) stage 4-5 and dialysis treatment on incidence of foot ulceration and major lower extremity amputation in comparison to CKD stage 3. Methods In this retrospective study, all individuals who visited our hospital between 2006 and 2012 because of CKD stages 3 to 5 or dialysis treatment were included. Medical records were reviewed for incidence of foot ulceration and major amputation. The time from CKD 3, CKD 4-5, and dialysis treatment until first foot ulceration and first major lower extremity amputation was calculated and analyzed by Kaplan-Meier curves and multivariate Cox proportional hazards model. Diabetes mellitus, peripheral arterial disease, peripheral neuropathy, and foot deformities were included for potential confounding. Results A total of 669 individuals were included: 539 in CKD 3, 540 in CKD 4-5, and 259 in dialysis treatment (individuals could progress from one group to the next). Unadjusted foot ulcer incidence rates per 1000 patients per year were 12 for CKD 3, 47 for CKD 4-5, and 104 for dialysis (P < .001). In multivariate analyses, the hazard ratio for incidence of foot ulceration was 4.0 (95% confidence interval [CI], 2.6-6.3) in CKD 4-5 and 7.6 (95% CI, 4.8-12.1) in dialysis treatment compared with CKD 3. Hazard ratios for incidence of major amputation were 9.5 (95% CI, 2.1-43.0) and 15 (95% CI, 3.3-71.0), respectively. Conclusions CKD 4-5 and dialysis treatment are independent risk factors for foot ulceration and major amputation compared with CKD 3. Maximum effort is needed in daily clinical practice to prevent foot ulcers and their devastating consequences in all individuals with CKD 4-5 or dialysis treatment.