Genetic risk and a primary role for cell-mediated immune mechanisms in multiple sclerosis.

Multiple sclerosis is a common disease of the central nervous system in which the interplay between inflammatory and neurodegenerative processes typically results in intermittent neurological disturbance followed by progressive accumulation of disability. Epidemiological studies have shown that genetic factors are primarily responsible for the substantially increased frequency of the disease seen in the relatives of affected individuals, and systematic attempts to identify linkage in multiplex families have confirmed that variation within the major histocompatibility complex (MHC) exerts the greatest individual effect on risk. Modestly powered genome-wide association studies (GWAS) have enabled more than 20 additional risk loci to be identified and have shown that multiple variants exerting modest individual effects have a key role in disease susceptibility. Most of the genetic architecture underlying susceptibility to the disease remains to be defined and is anticipated to require the analysis of sample sizes that are beyond the numbers currently available to individual research groups. In a collaborative GWAS involving 9,772 cases of European descent collected by 23 research groups working in 15 different countries, we have replicated almost all of the previously suggested associations and identified at least a further 29 novel susceptibility loci. Within the MHC we have refined the identity of the HLA-DRB1 risk alleles and confirmed that variation in the HLA-A gene underlies the independent protective effect attributable to the class I region. Immunologically relevant genes are significantly overrepresented among those mapping close to the identified loci and particularly implicate T-helper-cell differentiation in the pathogenesis of multiple sclerosis.

Integrating common and rare genetic variation in diverse human populations.

Despite great progress in identifying genetic variants that influence human disease, most inherited risk remains unexplained. A more complete understanding requires genome-wide studies that fully examine less common alleles in populations with a wide range of ancestry. To inform the design and interpretation of such studies, we genotyped 1.6 million common single nucleotide polymorphisms (SNPs) in 1,184 reference individuals from 11 global populations, and sequenced ten 100-kilobase regions in 692 of these individuals. This integrated data set of common and rare alleles, called 'HapMap 3', includes both SNPs and copy number polymorphisms (CNPs). We characterized population-specific differences among low-frequency variants, measured the improvement in imputation accuracy afforded by the larger reference panel, especially in imputing SNPs with a minor allele frequency of

A high-resolution HLA and SNP haplotype map for disease association studies in the extended human MHC.

The proteins encoded by the classical HLA class I and class II genes in the major histocompatibility complex (MHC) are highly polymorphic and are essential in self versus non-self immune recognition. HLA variation is a crucial determinant of transplant rejection and susceptibility to a large number of infectious and autoimmune diseases. Yet identification of causal variants is problematic owing to linkage disequilibrium that extends across multiple HLA and non-HLA genes in the MHC. We therefore set out to characterize the linkage disequilibrium patterns between the highly polymorphic HLA genes and background variation by typing the classical HLA genes and >7,500 common SNPs and deletion-insertion polymorphisms across four population samples. The analysis provides informative tag SNPs that capture much of the common variation in the MHC region and that could be used in disease association studies, and it provides new insight into the evolutionary dynamics and ancestral origins of the HLA loci and their haplotypes.

Human aging-associated DNA hypermethylation occurs preferentially at bivalent chromatin domains.

There is a growing realization that some aging-associated phenotypes/diseases have an epigenetic basis. Here, we report the first genome-scale study of epigenomic dynamics during normal human aging. We identify aging-associated differentially methylated regions (aDMRs) in whole blood in a discovery cohort, and then replicate these aDMRs in sorted CD4(+) T-cells and CD14(+) monocytes in an independent cohort, suggesting that aDMRs occur in precursor haematopoietic cells. Further replication of the aDMRs in buccal cells, representing a tissue that originates from a different germ layer compared with blood, demonstrates that the aDMR signature is a multitissue phenomenon. Moreover, we demonstrate that aging-associated DNA hypermethylation occurs predominantly at bivalent chromatin domain promoters. This same category of promoters, associated with key developmental genes, is frequently hypermethylated in cancers and in vitro cell culture, pointing to a novel mechanistic link between aberrant hypermethylation in cancer, aging, and cell culture.

A genome-wide survey of the prevalence and evolutionary forces acting on human nonsense SNPs.

Nonsense SNPs introduce premature termination codons into genes and can result in the absence of a gene product or in a truncated and potentially harmful protein, so they are often considered disadvantageous and are associated with disease susceptibility. As such, we might expect the disrupted allele to be rare and, in healthy people, observed only in a heterozygous state. However, some, like those in the CASP12 and ACTN3 genes, are known to be present at high frequencies and to occur often in a homozygous state and seem to have been advantageous in recent human evolution. To evaluate the selective forces acting on nonsense SNPs as a class, we have carried out a large-scale experimental survey of nonsense SNPs in the human genome by genotyping 805 of them (plus control synonymous SNPs) in 1,151 individuals from 56 worldwide populations. We identified 169 genes containing nonsense SNPs that were variable in our samples, of which 99 were found with both copies inactivated in at least one individual. We found that the sampled humans differ on average by 24 genes (out of about 20,000) because of these nonsense SNPs alone. As might be expected, nonsense SNPs as a class were found to be slightly disadvantageous over evolutionary timescales, but a few nevertheless showed signs of being possibly advantageous, as indicated by unusually high levels of population differentiation, long haplotypes, and/or high frequencies of derived alleles. This study underlines the extent of variation in gene content within humans and emphasizes the importance of understanding this type of variation.

A genome-wide association study confirms VKORC1, CYP2C9, and CYP4F2 as principal genetic determinants of warfarin dose.

We report the first genome-wide association study (GWAS) whose sample size (1,053 Swedish subjects) is sufficiently powered to detect genome-wide significance (p<1.5 x 10(-7)) for polymorphisms that modestly alter therapeutic warfarin dose. The anticoagulant drug warfarin is widely prescribed for reducing the risk of stroke, thrombosis, pulmonary embolism, and coronary malfunction. However, Caucasians vary widely (20-fold) in the dose needed for therapeutic anticoagulation, and hence prescribed doses may be too low (risking serious illness) or too high (risking severe bleeding). Prior work established that approximately 30% of the dose variance is explained by single nucleotide polymorphisms (SNPs) in the warfarin drug target VKORC1 and another approximately 12% by two non-synonymous SNPs (*2, *3) in the cytochrome P450 warfarin-metabolizing gene CYP2C9. We initially tested each of 325,997 GWAS SNPs for association with warfarin dose by univariate regression and found the strongest statistical signals (p<10(-78)) at SNPs clustering near VKORC1 and the second lowest p-values (p<10(-31)) emanating from CYP2C9. No other SNPs approached genome-wide significance. To enhance detection of weaker effects, we conducted multiple regression adjusting for known influences on warfarin dose (VKORC1, CYP2C9, age, gender) and identified a single SNP (rs2108622) with genome-wide significance (p = 8.3 x 10(-10)) that alters protein coding of the CYP4F2 gene. We confirmed this result in 588 additional Swedish patients (p<0.0029) and, during our investigation, a second group provided independent confirmation from a scan of warfarin-metabolizing genes. We also thoroughly investigated copy number variations, haplotypes, and imputed SNPs, but found no additional highly significant warfarin associations. We present power analysis of our GWAS that is generalizable to other studies, and conclude we had 80% power to detect genome-wide significance for common causative variants or markers explaining at least 1.5% of dose variance. These GWAS results provide further impetus for conducting large-scale trials assessing patient benefit from genotype-based forecasting of warfarin dose.

Common variants in the region around Osterix are associated with bone mineral density and growth in childhood.

Peak bone mass achieved in adolescence is a determinant of bone mass in later life. In order to identify genetic variants affecting bone mineral density (BMD), we performed a genome-wide association study of BMD and related traits in 1518 children from the Avon Longitudinal Study of Parents and Children (ALSPAC). We compared results with a scan of 134 adults with high or low hip BMD. We identified associations with BMD in an area of chromosome 12 containing the Osterix (SP7) locus, a transcription factor responsible for regulating osteoblast differentiation (ALSPAC: P = 5.8 x 10(-4); Australia: P = 3.7 x 10(-4)). This region has previously shown evidence of association with adult hip and lumbar spine BMD in an Icelandic population, as well as nominal association in a UK population. A meta-analysis of these existing studies revealed strong association between SNPs in the Osterix region and adult lumbar spine BMD (P = 9.9 x 10(-11)). In light of these findings, we genotyped a further 3692 individuals from ALSPAC who had whole body BMD and confirmed the association in children as well (P = 5.4 x 10(-5)). Moreover, all SNPs were related to height in ALSPAC children, but not weight or body mass index, and when height was included as a covariate in the regression equation, the association with total body BMD was attenuated. We conclude that genetic variants in the region of Osterix are associated with BMD in children and adults probably through primary effects on growth.

Variation within the gene encoding the upstream stimulatory factor 1 does not influence susceptibility to type 2 diabetes in samples from populations with replicated evidence of linkage to chromosome 1q.

The gene encoding the transcription factor upstream stimulatory factor (USF)1 influences susceptibility to familial combined hyperlipidemia (FCHL) and triglyceride levels. Phenotypic overlap between FCHL and type 2 diabetes makes USF1 a compelling positional candidate for the widely replicated type 2 diabetes linkage signal on chromosome 1q. We typed 22 variants in the F11R/USF1 region (1 per 3 kb), including those previously implicated in FCHL-susceptibility (or proxies thereof) in 3,726 samples preferentially enriched for 1q linkage. We also examined glucose- and lipid-related continuous traits in an overlapping set of 1,215 subjects of European descent. There was no convincing evidence for association with type 2 diabetes in any of seven case-control comparisons, individually or combined. Family-based association analyses in 832 Pima subjects were similarly negative. At rs3737787 (the variant most strongly associated with FCHL), the combined odds ratio, per copy of the rarer A-allele, was 1.10 (95% CI 0.97-1.24, P = 0.13). In 124 Utah subjects, rs3737787 was significantly associated (P = 0.002) with triglyceride levels, but direction of this association was opposite to previous reports, and there was no corroboration in three other samples. These data exclude USF1 as a major contributor to type 2 diabetes susceptibility and the basis for the chromosome 1q linkage. They reveal only limited evidence for replication of USF1 effects on continuous metabolic traits.

The brain-derived neurotrophic factor rs6265 (Val66Met) polymorphism and depression in Mexican-Americans.

The hypothesis that brain-derived neurotrophic factor (BDNF) is involved in the pathogenesis of major depression is supported by several research findings; however, genetic studies assessing the relationship between BDNF and psychiatric disorders have produced conflicting results. We examined the effect of a BDNF polymorphism on depression susceptibility in Mexican-Americans. The single nucleotide polymorphism (Val66Met), which has been shown to have functional and behavioral effects, was genotyped in 284 depressed participants and 331 controls, showing association with depression (P=0.005). Individuals homozygous for the major allele (GG) had an increased chance of being depressed (OR=1.7 95% CI 1.17-2.47). Our findings support the association of BDNF single nucleotide polymorphism rs6265 and depression, suggesting that this polymorphism may increase susceptibility to major depression in Mexican-Americans.

Australian midwives knowledge, attitude and perceived learning needs around perinatal mental health.

DESIGN AND AIM: a cross sectional survey was undertaken to explore midwives' knowledge of, and attitudes towards, mental health disorders in childbearing women vis-à-vis their perceived mental health learning needs. SETTING AND PARTICIPANTS: a 50.1% response rate included 238 midwives employed in the only public tertiary maternity hospital in Western Australia from March to June 2013. METHOD AND FINDINGS: The survey comprised a mixture of custom-designed questions and vignettes presenting various disorders. Only 37.6% of midwives felt well-equipped to support women, whilst 50.2% reported insufficient access to information. Demand was highest for education on: personality disorders (77.8%); the impact of childbearing on mental health disorders (74.2%); and skills for handling stress and aggression (57.8%). Knowledge scores were variable: on average eight out of a maximum 13 questions were answered correctly, but few (2.7%) answered more than 11 correctly, and 3.7% scored ≤4 correct. Across disorders, recognition from vignettes was highest for depression (93.9%), and lowest for schizophrenia (65.6%). Surprisingly, there were no associations between general knowledge scores and previous mental health experience, recent professional development, or access to information around mental health. The majority endorsed positive beliefs about midwives' role in mental health assessment, and belief in women's recovery (83.5%), however, cluster analysis of warmth and competence ratings revealed negative stereotyping of mental health disorders. KEY CONCLUSIONS: Midwives accept it is their role to assess the mental health status of women but many feel ill-equipped to do so and express a strong desire for further knowledge and skills across a range of perinatal mental health topics. Attitudes to recovery are positive but negative stereotypes exist; therefore awareness of potential bias is important to negate their influence on care. IMPLICATIONS FOR PRACTICE: Learning needs may change due to trends in clinical practice. Strategies are needed to recognise negative beliefs and to ensure education is responsive to local contexts.

Genome-wide association study of ankylosing spondylitis identifies non-MHC susceptibility loci.

To identify susceptibility loci for ankylosing spondylitis, we undertook a genome-wide association study in 2,053 unrelated ankylosing spondylitis cases among people of European descent and 5,140 ethnically matched controls, with replication in an independent cohort of 898 ankylosing spondylitis cases and 1,518 controls. Cases were genotyped with Illumina HumHap370 genotyping chips. In addition to strong association with the major histocompatibility complex (MHC; P < 10(-800)), we found association with SNPs in two gene deserts at 2p15 (rs10865331; combined P = 1.9 x 10(-19)) and 21q22 (rs2242944; P = 8.3 x 10(-20)), as well as in the genes ANTXR2 (rs4333130; P = 9.3 x 10(-8)) and IL1R2 (rs2310173; P = 4.8 x 10(-7)). We also replicated previously reported associations at IL23R (rs11209026; P = 9.1 x 10(-14)) and ERAP1 (rs27434; P = 5.3 x 10(-12)). This study reports four genetic loci associated with ankylosing spondylitis risk and identifies a major role for the interleukin (IL)-23 and IL-1 cytokine pathways in disease susceptibility.

Genome-wide and fine-resolution association analysis of malaria in West Africa.

We report a genome-wide association (GWA) study of severe malaria in The Gambia. The initial GWA scan included 2,500 children genotyped on the Affymetrix 500K GeneChip, and a replication study included 3,400 children. We used this to examine the performance of GWA methods in Africa. We found considerable population stratification, and also that signals of association at known malaria resistance loci were greatly attenuated owing to weak linkage disequilibrium (LD). To investigate possible solutions to the problem of low LD, we focused on the HbS locus, sequencing this region of the genome in 62 Gambian individuals and then using these data to conduct multipoint imputation in the GWA samples. This increased the signal of association, from P = 4 × 10(-7) to P = 4 × 10(-14), with the peak of the signal located precisely at the HbS causal variant. Our findings provide proof of principle that fine-resolution multipoint imputation, based on population-specific sequencing data, can substantially boost authentic GWA signals and enable fine mapping of causal variants in African populations.

Population-specific risk of type 2 diabetes conferred by HNF4A P2 promoter variants: a lesson for replication studies.

OBJECTIVE: Single nucleotide polymorphisms (SNPs) in the P2 promoter region of HNF4A were originally shown to be associated with predisposition for type 2 diabetes in Finnish, Ashkenazi, and, more recently, Scandinavian populations, but they generated conflicting results in additional populations. We aimed to investigate whether data from a large-scale mapping approach would replicate this association in novel Ashkenazi samples and in U.K. populations and whether these data would allow us to refine the association signal. RESEARCH DESIGN AND METHODS: Using a dense linkage disequilibrium map of 20q, we selected SNPs from a 10-Mb interval centered on HNF4A. In a staged approach, we first typed 4,608 SNPs in case-control populations from four U.K. populations and an Ashkenazi population (n = 2,516). In phase 2, a subset of 763 SNPs was genotyped in 2,513 additional samples from the same populations. RESULTS: Combined analysis of both phases demonstrated association between HNF4A P2 SNPs (rs1884613 and rs2144908) and type 2 diabetes in the Ashkenazim (n = 991; P < 1.6 x 10(-6)). Importantly, these associations are significant in a subset of Ashkenazi samples (n = 531) not previously tested for association with P2 SNPs (odds ratio [OR] approximately 1.7; P < 0.002), thus providing replication within the Ashkenazim. In the U.K. populations, this association was not significant (n = 4,022; P > 0.5), and the estimate for the OR was much smaller (OR 1.04; [95%CI 0.91-1.19]). CONCLUSIONS: These data indicate that the risk conferred by HNF4A P2 is significantly different between U.K. and Ashkenazi populations (P < 0.00007), suggesting that the underlying causal variant remains unidentified. Interactions with other genetic or environmental factors may also contribute to this difference in risk between populations.

Phosphodiesterase genes are associated with susceptibility to major depression and antidepressant treatment response.

Cyclic nucleotide phosphodiesterases (PDEs) constitute a family of enzymes that degrade cAMP and cGMP. Intracellular cyclic nucleotide levels increase in response to extracellular stimulation by hormones, neurotransmitters, or growth factors and are down-regulated through hydrolysis catalyzed by PDEs, which are therefore candidate therapeutic targets. cAMP is a second messenger implicated in learning, memory, and mood, and cGMP modulates nervous system processes that are controlled by the nitric oxide (NO)/cGMP pathway. To investigate an association between genes encoding PDEs and susceptibility to major depressive disorder (MDD), we genotyped SNPs in 21 genes of this superfamily in 284 depressed Mexican Americans who participated in a prospective, double-blind, pharmacogenetic study of antidepressant response, and 331 matched controls. Polymorphisms in PDE9A and PDE11A were found to be associated with the diagnosis of MDD. Our data are also suggestive of the association between SNPs in other PDE genes and MDD. Remission on antidepressants was significantly associated with polymorphisms in PDE1A and PDE11A. Thus, we found significant associations with both the diagnosis of MDD and remission in response to antidepressants with SNPs in the PDE11A gene. We show here that PDE11A haplotype GAACC is significantly associated with MDD. We conclude that PDE11A has a role in the pathophysiology of MDD. This study identifies a potential CNS role for the PDE11 family. The hypothesis that drugs affecting PDE function, particularly cGMP-related PDEs, represent a treatment strategy for major depression should therefore be tested.

A high-resolution linkage-disequilibrium map of the human major histocompatibility complex and first generation of tag single-nucleotide polymorphisms.

Autoimmune, inflammatory, and infectious diseases present a major burden to human health and are frequently associated with loci in the human major histocompatibility complex (MHC). Here, we report a high-resolution (1.9 kb) linkage-disequilibrium (LD) map of a 4.46-Mb fragment containing the MHC in U.S. pedigrees with northern and western European ancestry collected by the Centre d'Etude du Polymorphisme Humain (CEPH) and the first generation of haplotype tag single-nucleotide polymorphisms (tagSNPs) that provide up to a fivefold increase in genotyping efficiency for all future MHC-linked disease-association studies. The data confirm previously identified recombination hotspots in the class II region and allow the prediction of numerous novel hotspots in the class I and class III regions. The region of longest LD maps outside the classic MHC to the extended class I region spanning the MHC-linked olfactory-receptor gene cluster. The extended haplotype homozygosity analysis for recent positive selection shows that all 14 outlying haplotype variants map to a single extended haplotype, which most commonly bears HLA-DRB1*1501. The SNP data, haplotype blocks, and tagSNPs analysis reported here have been entered into a multidimensional Web-based database (GLOVAR), where they can be accessed and viewed in the context of relevant genome annotation. This LD map allowed us to give coordinates for the extremely variable LD structure underlying the MHC.

The impact of SNP density on fine-scale patterns of linkage disequilibrium.

Linkage disequilibrium (LD) is a measure of the degree of association between alleles in a population. The detection of disease-causing variants by association with neighbouring single nucleotide polymorphisms (SNPs) depends on the existence of strong LD between them. Previous studies have indicated that the extent of LD is highly variable in different chromosome regions and different populations, demonstrating the importance of genome-wide accurate measurement of LD at high resolution throughout the human genome. A uniform feature of these studies has been the inability to detect LD in regions of low marker density. To investigate the dependence of LD patterns on marker selection we performed a high-resolution study in African-American, Asian and UK Caucasian populations. We selected over 5000 SNPs with an average spacing of approximately 1 SNP per 2 kb after validating ca 12 000 SNPs derived from a dense SNP collection (1 SNP per 0.3 kb on average). Applications of different statistical methods of LD assessment highlight similar areas of high and low LD. However, at high resolution, features such as overall sequence coverage in LD blocks and block boundaries vary substantially with respect to marker density. Model-based linkage disequilibrium unit (LDU) maps appear robust to marker density and consistently influenced by marker allele frequency. The results suggest that very dense marker sets will be required to yield stable views of fine-scale LD in the human genome.