A genome-wide association study of sleep habits and insomnia

Several aspects of sleep behavior such as timing, duration and quality have been demonstrated to be heritable. To identify common variants that influence sleep traits in the population, we conducted a genome-wide association study of six sleep phenotypes assessed by questionnaire in a sample of 2,323 individuals from the Australian Twin Registry. Genotyping was performed on the Illumina 317, 370, and 610K arrays and the SNPs in common between platforms were used to impute non-genotyped SNPs. We tested for association with more than 2,000,000 common polymorphisms across the genome. While no SNPs reached the genome-wide significance threshold, we identified a number of associations in plausible candidate genes. Most notably, a group of SNPs in the third intron of the CACNA1C gene ranked as most significant in the analysis of sleep latency (P = 1.3 x 10(-)(6)). We attempted to replicate this association in an independent sample from the Chronogen Consortium (n = 2,034), but found no evidence of association (P = 0.73). We have identified several other suggestive associations that await replication in an independent sample. We did not replicate the results from previous genome-wide analyses of self-reported sleep phenotypes after correction for multiple testing.

Genome-wide association study of major depressive disorder: New results, meta-analysis, and lessons learned

Major depressive disorder (MDD) is a common complex disorder with a partly genetic etiology. We conducted a genome-wide association study of the MDD2000+ sample (2431 cases, 3673 screened controls and >1 M imputed single-nucleotide polymorphisms (SNPs)). No SNPs achieved genome-wide significance either in the MDD2000+ study, or in meta-analysis with two other studies totaling 5763 cases and 6901 controls. These results imply that common variants of intermediate or large effect do not have main effects in the genetic architecture of MDD. Suggestive but notable results were: (a) gene-based tests suggesting roles for adenylate cyclase 3 (ADCY3, 2p23.3) and galanin (GAL, 11q13.3); published functional evidence relates both of these to MDD and serotonergic signaling; (b) support for the bipolar disorder risk variant SNP rs1006737 in CACNA1C (P=0.020, odds ratio=1.10), and; (c) lack of support for rs2251219, a SNP identified in a meta-analysis of affective disorder studies (P=0.51). We estimate that sample sizes 1.8- to 2.4-fold greater are needed for association studies of MDD compared with those for schizophrenia to detect variants that explain the same proportion of total variance in liability. Larger study cohorts characterized for genetic and environmental risk factors accumulated prospectively are likely to be needed to dissect more fully the etiology of MDD.

Pharmacogenetic Variation at CYP2D6, CYP2C9, and CYP2C19: Population Genetic and Forensic Aspects

Pharmacogenetics deals with genetically determined variation in drug response. In this context, three phase I drug-metabolizing enzymes, CYP2D6, CYP2C9, and CYP2C19, have a central role, affecting the metabolism of about 20-30% of clinically used drugs. Since genes coding for these enzymes in human populations exhibit high genetic polymorphism, they are of major pharmacogenetic importance. The aims of this study were to develop new genotyping methods for CYP2D6, CYP2C9, and CYP2C19 that would cover the most important genetic variants altering the enzyme activity, and, for the first time, to describe the distribution of genetic variation at these loci on global and microgeographic scales. In addition, pharmacogenetics was applied to a postmortem forensic setting to elucidate the role of genetic variation in drug intoxications, focusing mainly on cases related to tricyclic antidepressants, which are commonly involved in fatal drug poisonings in Finland. Genetic variability data were obtained by genotyping new population samples by the methods developed based on PCR and multiplex single-nucleotide primer extension reaction, as well as by collecting data from the literature. Data consisted of 138, 129, and 146 population samples for CYP2D6, CYP2C9, and CYP2C19, respectively. In addition, over 200 postmortem forensic cases were examined with respect to drug and metabolite concentrations and genotypic variation at CYP2D6 and CYP2C19. The distribution of genetic variation within and among human populations was analyzed by descriptive statistics and variance analysis and by correlating the genetic and geographic distances using Mantel tests and spatial autocorrelation. The correlation between phenotypic and genotypic variation in drug metabolism observed in postmortem cases was also analyzed statistically. The genotyping methods developed proved to be informative, technically feasible, and cost-effective. Detailed molecular analysis of CYP2D6 genetic variation in a global survey of human populations revealed that the pattern of variation was similar to those of neutral genomic markers. Most of the CYP2D6 diversity was observed within populations, and the spatial pattern of variation was best described as clinal. On the other hand, genetic variants of CYP2D6, CYP2C9, and CYP2C19 associated with altered enzymatic activity could reach extremely high frequencies in certain geographic regions. Pharmacogenetic variation may also be significantly affected by population-specific demographic histories, as seen within the Finnish population. When pharmacogenetics was applied to a postmortem forensic setting, a correlation between amitriptyline metabolic ratios and genetic variation at CYP2D6 and CYP2C19 was observed in the sample material, even in the presence of confounding factors typical for these cases. In addition, a case of doxepin-related fatal poisoning was shown to be associated with a genetic defect at CYP2D6. Each of the genes studied showed a distinct variation pattern in human populations and high frequencies of altered activity variants, which may reflect the neutral evolution and/or selective pressures caused by dietary or environmental exposure. The results are relevant also from the clinical point of view since the genetic variation at CYP2D6, CYP2C9, and CYP2C19 already has a range of clinical applications, e.g. in cancer treatment and oral anticoagulation therapy. This study revealed that pharmacogenetics may also contribute valuable information to the medicolegal investigation of sudden, unexpected deaths.

Genome-wide association analysis identifies susceptibility loci for migraine without aura

Migraine without aura is the most common form of migraine, characterized by recurrent disabling headache and associated autonomic symptoms. To identify common genetic variants associated with this migraine type, we analyzed genome-wide association data of 2,326 clinic-based German and Dutch individuals with migraine without aura and 4,580 population-matched controls. We selected SNPs from 12 loci with 2 or more SNPs associated with P values of <1 x 10(-5) for replication testing in 2,508 individuals with migraine without aura and 2,652 controls. SNPs at two of these loci showed convincing replication: at 1q22 (in MEF2D; replication P = 4.9 x 10(-4); combined P = 7.06 x 10(-11)) and at 3p24 (near TGFBR2; replication P = 1.0 x 10(-4); combined P = 1.17 x 10(-9)). In addition, SNPs at the PHACTR1 and ASTN2 loci showed suggestive evidence of replication (P = 0.01; combined P = 3.20 x 10(-8) and P = 0.02; combined P = 3.86 x 10(-8), respectively). We also replicated associations at two previously reported migraine loci in or near TRPM8 and LRP1. This study identifies the first susceptibility loci for migraine without aura, thereby expanding our knowledge of this debilitating neurological disorder.

A genome-wide association study of caffeine-related sleep disturbance: Confirmation of a role for a common variant in the adenosine receptor

OBJECTIVES To identify common genetic variants that predispose to caffeine-induced insomnia and to test whether genes whose expression changes in the presence of caffeine are enriched for association with caffeine-induced insomnia. DESIGN A hypothesis-free, genome-wide association study. SETTING Community-based sample of Australian twins from the Australian Twin Registry. PARTICIPANTS After removal of individuals who said that they do not drink coffee, a total of 2,402 individuals from 1,470 families in the Australian Twin Registry provided both phenotype and genotype information. MEASUREMENTS AND RESULTS A dichotomized scale based on whether participants reported ever or never experiencing caffeine-induced insomnia. A factor score based on responses to a number of questions regarding normal sleep habits was included as a covariate in the analysis. More than 2 million common single nucleotide polymorphisms (SNPs) were tested for association with caffeine-induced insomnia. No SNPs reached the genome-wide significance threshold. In the analysis that did not include the insomnia factor score as a covariate, the most significant SNP identified was an intronic SNP in the PRIMA1 gene (P = 1.4 x 10(-)(6), odds ratio = 0.68 [0.53 - 0.89]). An intergenic SNP near the GBP4 gene on chromosome 1 was the most significant upon inclusion of the insomnia factor score into the model (P = 1.9 x 10(-)(6), odds ratio = 0.70 [0.62 - 0.78]). A previously identified association with a polymorphism in the ADORA2A gene was replicated. CONCLUSIONS Several genes have been identified in the study as potentially influencing caffeine-induced insomnia. They will require replication in another sample. The results may have implications for understanding the biologic mechanisms underlying insomnia.

Genome-wide meta-analysis of common variant differences between men and women

The male-to-female sex ratio at birth is constant across world populations with an average of 1.06 (106 male to 100 female live births) for populations of European descent. The sex ratio is considered to be affected by numerous biological and environmental factors and to have a heritable component. The aim of this study was to investigate the presence of common allele modest effects at autosomal and chromosome X variants that could explain the observed sex ratio at birth. We conducted a large-scale genome-wide association scan (GWAS) meta-analysis across 51 studies, comprising overall 114 863 individuals (61 094 women and 53 769 men) of European ancestry and 2 623 828 common (minor allele frequency >0.05) single-nucleotide polymorphisms (SNPs). Allele frequencies were compared between men and women for directly-typed and imputed variants within each study. Forward-time simulations for unlinked, neutral, autosomal, common loci were performed under the demographic model for European populations with a fixed sex ratio and a random mating scheme to assess the probability of detecting significant allele frequency differences. We do not detect any genome-wide significant (P < 5 x 10(-8)) common SNP differences between men and women in this well-powered meta-analysis. The simulated data provided results entirely consistent with these findings. This large-scale investigation across ~115 000 individuals shows no detectable contribution from common genetic variants to the observed skew in the sex ratio. The absence of sex-specific differences is useful in guiding genetic association study design, for example when using mixed controls for sex-biased traits.

The Functions and Regulation of Ornithine Decarboxylase

Ornithine decarboxylase (ODC) regulates the synthesis of polyamines which are involved in many cellular functions e.g. proliferation and differentiation. Due to its critical role, ODC is a tightly regulated enzyme by antizymes and antizyme inhibitors. If the regulation fails, the activity of ODC increases and may lead to malignant transformation of a cell. Increased ODC activity is found in many common cancers, including colon, prostate, and breast cancer. In a transformed cell, dynamics of the actin cytoskeleton is disturbed. A small G-protein, RhoA regulates organization of the cytoskeleton, and its overactivity increases malignant potential of the cell. The present results indicate that covalent attachment of polyamines by transglutaminase is a physiological means of regulating the activity of RhoA. The translocation of RhoA to the plasma membrane, where it exerts its activity is dependent on the presence of catalytically active ODC. As the overactivity of ODC and RhoA are implicated in cell transformation, the results provide a mechanistic explanation of the interrelationship between the polyamine metabolism and the reorganization of the actin cytoskeleton occurring in cancer cells. ODC and polyamines have also an important role in the function of central nervous system. They participate in the regulation of brain morphogenesis in embryos. In adult nervous tissue, polyamines regulate K+ and glutamate channels. K+ inward rectifying channels control membrane potentials and NMDA-type glutamate receptors (NMDAR) regulate synaptic plasticity. High ODC activity and polyamine levels are considered important in the development of ischemic brain damage and they are implicated in the pathogenesis of Alzheimer s disease (AD). A homolog of ODC was cloned from a human brain cDNA library, and several alternatively spliced variants were detected in human brain and testis. The novel protein was nevertheless devoid of ODC catalytic activity. It was subsequently found to be a novel inductor of ODC activity and polyamine synthesis, called antizyme inhibitor 2 (AZIN2). The accumulation of AZIN2 in vesicle-like formations along the axons and beneath the plasma membrane of neurons as well as in steroid hormone producing Leydig cells and luteal cells of the gonads implies that AZIN2 plays a role in secretion and vesicle trafficking. An accumulation of AZIN2 was detected also in specimens of AD brains. This increased expression of AZIN2 was specific for AD and was not found in brains with other neurodegenerative diseases including CADASIL or dementia with Lewy bodies.

Palladin, a novel microfilament protein

Palladin is a novel actin microfilament associated protein, which together with myotilin and myopalladin forms a novel cytoskeletal IgC2 domain protein family. Whereas the expression of myotilin and myopalladin is limited mainly to striated muscle, palladin is widely expressed in both epithelial and mesenchymal tissues, including heart and the nervous system. Palladin has a complex genetic structure and it is expressed as several different sized and structured splice variants, which also display differences in their expression pattern and interactions. In muscle cells, all the family members localize to the sarcomeric Z-disc, and in non-muscle cells palladin also localizes to the stress-fiber-dense regions, lamellipodia, podosomes and focal adhesions. A common feature of this protein family is the binding to α-actinin, but other interactions are mostly unique to each member. Palladin has been shown to interact with several proteins, including VASP, profilin, Eps8, LASP-1 and LPP. Its domain structure, lack of enzymatic activity and multiple interactions define it as a molecular scaffolding protein, which links together proteins with different functional modalities into large complexes. Palladin has an important role in cytoskeletal regulation, particularly in stress fiber formation and stabilization. This assumption is supported by several experimental results. First, over-expression of palladin in non-muscle cells results in rapid reorganization of the actin cytoskeleton and formation of thick actin bundles. Second, the knock-down of palladin with anti-sense and siRNA techniques or knock-out by genetic methods leads to defective stress fiber formation. Furthermore, palladin is usually up-regulated in situations requiring a highly organized cytoskeleton, such as differentiation of dendritic cells, trophoblasts and myofibroblasts, and activation of astrocytes during glial scar formation. The protein family members have also direct disease linkages; myotilin missense mutations are the cause of LGMD1A and myofibrillar myopathy. Palladin mutations and polymorphisms, on the other hand, have been linked to hereditary pancreatic cancer and myocardial infarction, respectively. In this study we set out to characterize human palladin. We identified several palladin isoforms, studied their tissue distribution and sub-cellular localization. Four novel interaction partners were identified; ezrin, ArgBP2, SPIN90 and Src-kinase.The previously identified interaction between palladin and α-actinin was also characterized in detail. All the identified new binding partners are actin cytoskeleton associated proteins; ezrin links the plasma membrane to the cytoskeleton, ArgBP2 and SPIN90 localize, among other structures, to the lamellipodia and in cardiomyocytes to the Z-disc. Src is a transforming tyrosine kinase, which besides its role in oncogenesis has also important cytoskeletal associations. We also studied palladin in myofibroblasts, which are specialized cells involved in diverse physiological and pathological processes, such as wound healing and tissue fibrosis. We demonstrated that palladin is up-regulated during the differentiation of myofibroblasts in an isoform specific manner, and that this up-regulation is induced by TGF-β via activation of both the SMAD and MAPK signalling cascades. In summary, the results presented here describe the initial characterization of human palladin and offer a basis for further studies.

FTO genotype is associated with phenotypic variability of body mass index

There is evidence across several species for genetic control of phenotypic variation of complex traits1, 2, 3, 4, such that the variance among phenotypes is genotype dependent. Understanding genetic control of variability is important in evolutionary biology, agricultural selection programmes and human medicine, yet for complex traits, no individual genetic variants associated with variance, as opposed to the mean, have been identified. Here we perform a meta-analysis of genome-wide association studies of phenotypic variation using ~170,000 samples on height and body mass index (BMI) in human populations. We report evidence that the single nucleotide polymorphism (SNP) rs7202116 at the FTO gene locus, which is known to be associated with obesity (as measured by mean BMI for each rs7202116 genotype)5, 6, 7, is also associated with phenotypic variability. We show that the results are not due to scale effects or other artefacts, and find no other experiment-wise significant evidence for effects on variability, either at loci other than FTO for BMI or at any locus for height. The difference in variance for BMI among individuals with opposite homozygous genotypes at the FTO locus is approximately 7%, corresponding to a difference of ~0.5 kilograms in the standard deviation of weight. Our results indicate that genetic variants can be discovered that are associated with variability, and that between-person variability in obesity can partly be explained by the genotype at the FTO locus. The results are consistent with reported FTO by environment interactions for BMI8, possibly mediated by DNA methylation9, 10. Our BMI results for other SNPs and our height results for all SNPs suggest that most genetic variants, including those that influence mean height or mean BMI, are not associated with phenotypic variance, or that their effects on variability are too small to detect even with samples sizes greater than 100,000.

Molecular Background of Common Dyslipidemias

The leading cause of death in the Western world continues to be coronary heart disease (CHD). At the root of the disease process is dyslipidemia an aberration in the relevant amounts of circulating blood lipids. Cholesterol builds up in the arterial wall and following rupture of these plaques, myocardial infarction or stroke can occur. Heart disease runs in families and a number of hereditary forms are known. The leading cause of adult dyslipidemia presently however is overweight and obesity. This thesis work presents an investigation of the molecular genetics of common, hereditary dyslipidemia and the tightly related condition of obesity. Familial combined hyperlipidemia (FCHL) is the most common hereditary dyslipidemia in man with an estimated population prevalence of 1-6%. This complex disease is characterized by elevated levels of serum total cholesterol, triglycerides or both and is observed in about 20% of individuals with premature CHD. Our group identified the disease to be associated with genetic variation in the USF1 transcription factor gene. USF1 has a key role in regulating other genes that control lipid and glucose metabolism as well as the inflammatory response all central processes in the progression of atherosclerosis and CHD. The first two works of this thesis aimed at understanding how these USF1 variants result in increased disease risk. Among the many, non-coding single-nucleotide polymorphisms (SNPs) that associated with the disease, one was found to have a functional effect. The risk-enhancing allele of this SNP seems to eradicate the ability of the important hormone insulin to induce the expression of USF1 in peripheral tissues. The resultant changes in the expression of numerous USF1 target genes over time probably enhance and accelerate the atherogenic processes. Dyslipidemias often represent an outcome of obesity and in the final work of this thesis we wanted to address the metabolic pathways related to acquired obesity. It is recognized that active processes in adipose tissue play an important role in the development of dyslipidemia, insulin resistance and other pathological conditions associated with obesity. To minimize the confounding effects of genetic differences present in most human studies, we investigated a rare collection of identical twins that differed significantly in the amount of body fat. In the obese, but otherwise healthy young adults, several notable changes were observed. In addition to chronic inflammation, the adipose tissue of the obese co-twins was characterized by a marked (47%) decrease in amount of mitochondrial DNA (mtDNA) a change associated with mitochondrial dysfunction. The catabolism of branched chain amino acids (BCAAs) was identified as the most down-regulated process in the obese co-twins. A concordant increase in the serum level of these insulin secretagogues was identified. This hyperaminoacidemia may provide the feed-back signal from insulin resistant adipose tissue to the pancreas to ensure an appropriately augmented secretory response. The down regulation of BCAA catabolism correlated closely with liver fat accumulation and insulin. The single most up-regulated gene (5.9 fold) in the obese co-twins was osteopontin (SPP1) a cytokine involved in macrophage recruitment to adipose tissue. SPP1 is here implicated as an important player in the development of insulin resistance. These studies of exceptional study samples provide better understanding of the underlying pathology in common dyslipidemias and other obesity associated diseases important for future improvement of intervention strategies and treatments to combat atherosclerosis and coronary heart disease.

Evaluation of polymorphisms in predicted target sites for micro RNAs differentially expressed in endometriosis

Previous microarray analyses identified 22 microRNAs (miRNAs) differentially expressed in paired ectopic and eutopic endometrium of women with and without endometriosis. To investigate further the role of these miRNAs in women with endometriosis, we conducted an association study aiming to explore the relationship between endometriosis risk and single-nucleotide polymorphisms (SNPs) in miRNA target sites for these differentially expressed miRNAs. A panel of 102 SNPs in the predicted miRNA binding sites were evaluated for an endometriosis association study and an ingenuity pathway analysis was performed. Fourteen rare variants were identified in this study. We found SNP rs14647 in the Wolf-Hirschhorn syndrome candidate gene1 (WHSC1) 3'UTR (untranslated region) was associated with endometriosis-related infertility presenting an odds ratio of 12.2 (95% confidence interval = 2.4-60.7, P = 9.03 x 10(-5)). SNP haplotype AGG in the solute carrier family 22, member 23 (SLC22A23) 3'UTR was associated with endometriosis-related infertility and more severe disease. With the individual genotyping data, ingenuity pathways analysis identified the tumour necrosis factor and cyclin-dependant kinase inhibitor as major factors in the molecular pathways. Significant associations between WHSC1 alleles and endometriosis-related infertility and SLC22A23 haplotypes and the disease severe stage were identified. These findings may help focus future research on subphenotypes of this disease. Replication studies in independent large sample sets to confirm and characterize the involvement of the gene variation in the pathogenesis of endometriosis are needed.

Genomic inflation factors under polygenic inheritance

Population structure, including population stratification and cryptic relatedness, can cause spurious associations in genome-wide association studies (GWAS). Usually, the scaled median or mean test statistic for association calculated from multiple single-nucleotide-polymorphisms across the genome is used to assess such effects, and 'genomic control' can be applied subsequently to adjust test statistics at individual loci by a genomic inflation factor. Published GWAS have clearly shown that there are many loci underlying genetic variation for a wide range of complex diseases and traits, implying that a substantial proportion of the genome should show inflation of the test statistic. Here, we show by theory, simulation and analysis of data that in the absence of population structure and other technical artefacts, but in the presence of polygenic inheritance, substantial genomic inflation is expected. Its magnitude depends on sample size, heritability, linkage disequilibrium structure and the number of causal variants. Our predictions are consistent with empirical observations on height in independent samples of ~4000 and ~133,000 individuals.

Methicillin Resistance In Staphylococci : Horizontal Transfer of Mobile Genetic Element (SCCmec) Between Staphylococcal Species

Regardless of the existence of antibiotics, infectious diseases are the leading causes of death in the world. Staphylococci cause many infections of varying severity, although they can also exist peacefully in many parts of the human body. Most often Staphylococcus aureus colonises the nose, and that colonisation is considered to be a risk factor for spread of this bacterium. S. aureus is considered to be the most important Staphylococcus species. It poses a challenge to the field of medicine, and one of the most problematic aspects is the drastic increase of the methicillin-resistant S. aureus (MRSA) strains in hospitals and community world-wide, including Finland. In addition, most of the clinical coagulase-negative staphylococcus (CNS) isolates express resistance to methicillin. Methicillin-resistance in S. aureus is caused by the mecA gene that encodes an extra penicillin-binding protein (PBP) 2a. The mecA gene is found in a mobile genomic island called staphylococcal chromosome cassette mec (SCCmec). The SCCmec consists of the mec gene and cassette chromosome recombinase (ccr)gene complexes. The areas of the SCCmec element outside the ccr and mec complex are known as the junkyard J regions. So far, eight types of SCCmec(SCCmec I- SCCmec VIII) and a number of variants have been described. The SCCmec island is an acquired element in S. aureus. Lately, it appears that CNS might be the storage place of the SCCmec that aid the S. aureus by providing it with the resistant elements. The SCCmec is known to exist only in the staphylococci. The aim of the present study was to investigate the horizontal transfer of SCCmec between the S. aureus and CNS. One specific aim was to study whether or not some methicillin-sensitive S. aureus (MSSA) strains are more inclined to receive the SCCmec than others. This was done by comparing the genetic background of clinical MSSA isolates in the health care facilities of the Helsinki and Uusimaa Hospital District in 2001 to the representatives of the epidemic MRSA (EMRSA) genotypes, which have been encountered in Finland during 1992-2004. Majority of the clinical MSSA strains were related to the EMRSA strains. This finding suggests that horizontal transfer of SCCmec from unknown donor(s) to several MSSA background genotypes has occurred in Finland. The molecular characteristics of representative clinical methicillin-resistant S. epidermidis (MRSE) isolates recovered in Finnish hospitals between 1990 and 1998 were also studied, examining their genetic relation to each other and to the internationally recognised MRSE clones as well, so as to ascertain the common traits between the SCCmec elements in MRSE and MRSA. The clinical MRSE strains were genetically related to each other; eleven PFGE types were associated with sequence type ST2 that has been identified world-wide. A single MRSE strain may possess two SCCmec types III and IV, which were recognised among the MRSA strains. Moreover, six months after the onset of an outbreak of MRSA possessing a SCCmec type V in a long-term care facility in Northern Finland (LTCF) in 2003, the SCCmec element of nasally carried methicillin-resistant staphylococci was studied. Among the residents of a LTCF, nasal carriage of MR-CNS was common with extreme diversity of SCCmec types. MRSE was the most prevalent CNS species. Horizontal transfer of SCCmec elements is speculated to be based on the sharing of SCCmec type V between MRSA and MRSE in the same person. Additionally, the SCCmec element of the clinical human S. sciuri isolates was studied. Some of the SCCmec regions were present in S. sciuri and the pls gene was common in it. This finding supports the hypothesis of genetic exchange happening between staphylococcal species. Evaluation of the epidemiology of methicillin-resistant staphylococcal colonisation is necessary in order to understand the apparent emergence of these strains and to develop appropriate control strategies. SCCmec typing is essential for understanding the emergence of MRSA strains from CNS, considering that the MR-CNS may represent the gene pool for the continuous creation of new SCCmec types from which MRSA might originate.

Association mapping

Association mapping seeks to identify marker alleles present at significantly different frequencies in cases carrying a particular disease or trait compared with controls. Genome-wide association studies are increasingly replacing candidate gene-based association studies for complex diseases, where a number of loci are likely to contribute to disease risk and the effect size of each particular risk allele is typically modest or low. Good study design is essential to the success of an association study, and factors such as the heritability of the disease under investigation, the choice of controls, statistical power, multiple testing and whether the association can be replicated need to be considered before beginning. Likewise, thorough quality control of the genotype data needs to be undertaken prior to running any association analyses. Finally, it should be kept in mind that a significant genetic association is not proof positive that a particular genetic locus causes a disease, but rather an important first step in discovering the genetic variants underlying a complex disease.

Genome-wide association study identifies a new melanoma susceptibility locus at 1q21.3

We performed a genome-wide association study of melanoma in a discovery cohort of 2,168 Australian individuals with melanoma and 4,387 control individuals. In this discovery phase, we confirm several previously characterized melanoma-associated loci at MC1R, ASIP and MTAP-CDKN2A. We selected variants at nine loci for replication in three independent case-control studies (Europe: 2,804 subjects with melanoma, 7,618 control subjects; United States 1: 1,804 subjects with melanoma, 1,026 control subjects; United States 2: 585 subjects with melanoma, 6,500 control subjects). The combined meta-analysis of all case-control studies identified a new susceptibility locus at 1q21.3 (rs7412746, P = 9.0 x 10(-11), OR in combined replication cohorts of 0.89 (95% CI 0.85-0.95)). We also show evidence suggesting that melanoma associates with 1q42.12 (rs3219090, P = 9.3 x 10(-8)). The associated variants at the 1q21.3 locus span a region with ten genes, and plausible candidate genes for melanoma susceptibility include ARNT and SETDB1. Variants at the 1q21.3 locus do not seem to be associated with human pigmentation or measures of nevus density.