Association of a functional variant downstream of TNFAIP3 with systemic lupus erythematosus

Systemic lupus erythematosus (SLE, MIM152700) is an autoimmune disease characterized by self-reactive antibodies resulting in systemic inflammation and organ failure. TNFAIP3, encoding the ubiquitin-modifying enzyme A20, is an established susceptibility locus for SLE. By fine mapping and genomic re-sequencing in ethnically diverse populations, we fully characterized the TNFAIP3 risk haplotype and identified a TT>A polymorphic dinucleotide (deletion T followed by a T to A transversion) associated with SLE in subjects of European (P = 1.58 × 10(-8), odds ratio = 1.70) and Korean (P = 8.33 × 10(-10), odds ratio = 2.54) ancestry. This variant, located in a region of high conservation and regulatory potential, bound a nuclear protein complex composed of NF-κB subunits with reduced avidity. Further, compared with the non-risk haplotype, the haplotype carrying this variant resulted in reduced TNFAIP3 mRNA and A20 protein expression. These results establish this TT>A variant as the most likely functional polymorphism responsible for the association between TNFAIP3 and SLE.     

Genetic variation in the 5q31 cytokine gene cluster confers susceptibility to Crohn disease

Linkage disequilibrium (LD) mapping provides a powerful method for fine-structure localization of rare disease genes, but has not yet been widely applied to common disease. We sought to design a systematic approach for LD mapping and apply it to the localization of a gene (IBD5) conferring susceptibility to Crohn disease. The key issues are: (i) to detect a significant LD signal (ii) to rigorously bound the critical region and (iii) to identify the causal genetic variant within this region. We previously mapped the IBD5 locus to a large region spanning 18 cM of chromosome 5q31 (P<10(-4)). Using dense genetic maps of microsatellite markers and single-nucleotide polymorphisms (SNPs) across the entire region, we found strong evidence of LD. We bound the region to a common haplotype spanning 250 kb that shows strong association with the disease (P< 2 x 10(-7)) and contains the cytokine gene cluster. This finding provides overwhelming evidence that a specific common haplotype of the cytokine region in 5q31 confers susceptibility to Crohn disease. However, genetic evidence alone is not sufficient to identify the causal mutation within this region, as strong LD across the region results in multiple SNPs having equivalent genetic evidence-each consistent with the expected properties of the IBD5 locus. These results have important implications for Crohn disease in particular and LD mapping in general.     

The PHYTOCHROME C photoreceptor gene mediates natural variation in flowering and growth responses of Arabidopsis thaliana

Light has an important role in modulating seedling growth and flowering time. We show that allelic variation at the PHYTOCHROME C (PHYC) photoreceptor locus affects both traits in natural populations of A. thaliana. Two functionally distinct PHYC haplotype groups are distributed in a latitudinal cline dependent on FRIGIDA, a locus that together with FLOWERING LOCUS C explains a large portion of the variation in A. thaliana flowering time. In a genome-wide scan for association of 65 loci with latitude, there was an excess of significant P values, indicative of population structure. Nevertheless, PHYC was the most strongly associated locus across 163 strains, suggesting that PHYC alleles are under diversifying selection in A. thaliana. Our work, together with previous findings, suggests that photoreceptor genes are major agents of natural variation in plant flowering and growth response.     

Evidence for two independent prostate cancer risk-associated loci in the HNF1B gene at 17q12

We carried out a fine-mapping study in the HNF1B gene at 17q12 in two study populations and identified a second locus associated with prostate cancer risk, approximately 26 kb centromeric to the first known locus (rs4430796); these loci are separated by a recombination hot spot. We confirmed the association with a SNP in the second locus (rs11649743) in five additional populations, with P = 1.7 x 10(-9) for an allelic test of the seven studies combined. The association at each SNP remained significant after adjustment for the other SNP.     

A worldwide survey of haplotype variation and linkage disequilibrium in the human genome

Recent genomic surveys have produced high-resolution haplotype information, but only in a small number of human populations. We report haplotype structure across 12 Mb of DNA sequence in 927 individuals representing 52 populations. The geographic distribution of haplotypes reflects human history, with a loss of haplotype diversity as distance increases from Africa. Although the extent of linkage disequilibrium (LD) varies markedly across populations, considerable sharing of haplotype structure exists, and inferred recombination hotspot locations generally match across groups. The four samples in the International HapMap Project contain the majority of common haplotypes found in most populations: averaging across populations, 83% of common 20-kb haplotypes in a population are also common in the most similar HapMap sample. Consequently, although the portability of tag SNPs based on the HapMap is reduced in low-LD Africans, the HapMap will be helpful for the design of genome-wide association mapping studies in nearly all human populations.     

A new gene, encoding an anion transporter, is mutated in sialic acid storage diseases

Sialic acid storage diseases (SASD, MIM 269920) are autosomal recessive neurodegenerative disorders that may present as a severe infantile form (ISSD) or a slowly progressive adult form, which is prevalent in Finland (Salla disease). The main symptoms are hypotonia, cerebellar ataxia and mental retardation; visceromegaly and coarse features are also present in infantile cases. Progressive cerebellar atrophy and dysmyelination have been documented by magnetic resonance imaging (ref. 4). Enlarged lysosomes are seen on electron microscopic studies and patients excrete large amounts of free sialic acid in urine. A H+/anionic sugar symporter mechanism for sialic acid and glucuronic acid is impaired in lysosomal membranes from Salla and ISSD patients. The locus for Salla disease was assigned to a region of approximately 200 kb on chromosome 6q14-q15 in a linkage study using Finnish families. Salla disease and ISSD were further shown to be allelic disorders. A physical map with P1 and PAC clones was constructed to cover the 200-kb area flanked by the loci D6S280 and D6S1622, providing the basis for precise physical positioning of the gene. Here we describe a new gene, SLC17A5 (also known as AST), encoding a protein (sialin) with a predicted transport function that belongs to a family of anion/cation symporters (ACS). We found a homozygous SLC17A5 mutation (R39C) in five Finnish patients with Salla disease and six different SLC17A5 mutations in six ISSD patients of different ethnic origins. Our observations suggest that mutations in SLC17A5 are the primary cause of lysosomal sialic acid storage diseases.     

Chromosome 4q DNA rearrangements associated with facioscapulohumeral muscular dystrophy.

Facioscapulohumeral muscular dystrophy (FSHD) is an autosomal dominant neuromuscular disorder which maps to chromosome 4qter, distal to the D4S139 locus. The cosmid clone 13E, isolated in a search for homeobox genes, was subsequently mapped to 4q35, also distal to D4S139. A subclone, p13E-11, detects in normal individuals a polymorphic EcoRI fragment usually larger than 28 kilobases (kb). Surprisingly, using the same probe we detected de novo DNA rearrangements, characterized by shorter EcoRI fragments (14-28 kb), in 5 out of 6 new FSHD cases. In 10 Dutch families analysed, a specific shorter fragment between 14-28 kb cosegregates with FSHD. Both observations indicate that FSHD is caused by independent de novo DNA rearrangements in the EcoRI fragment detected by p13E-11.     

A predominant role for the HLA class II region in the association of the MHC region with multiple sclerosis

Genetic susceptibility to multiple sclerosis is associated with genes of the major histocompatibility complex (MHC), particularly HLA-DRB1 and HLA-DQB1 (ref. 1). Both locus and allelic heterogeneity have been reported in this genomic region. To clarify whether HLA-DRB1 itself, nearby genes in the region encoding the MHC or combinations of these loci underlie susceptibility to multiple sclerosis, we genotyped 1,185 Canadian and Finnish families with multiple sclerosis (n = 4,203 individuals) with a high-density SNP panel spanning the genes encoding the MHC and flanking genomic regions. Strong associations in Canadian and Finnish samples were observed with blocks in the HLA class II genomic region (P < 4.9 x 10(-13) and P < 2.0 x 10(-16), respectively), but the strongest association was with HLA-DRB1 (P < 4.4 x 10(-17)). Conditioning on either HLA-DRB1 or the most significant HLA class II haplotype block found no additional block or SNP association independent of the HLA class II genomic region. This study therefore indicates that MHC-associated susceptibility to multiple sclerosis is determined by HLA class II alleles, their interactions and closely neighboring variants.     

Positional cloning of the combined hyperlipidemia gene Hyplip1.

Familial combined hyperlipidemia (FCHL, MIM-144250) is a common, multifactorial and heterogeneous dyslipidemia predisposing to premature coronary artery disease and characterized by elevated plasma triglycerides, cholesterol, or both. We identified a mutant mouse strain, HcB-19/Dem (HcB-19), that shares features with FCHL, including hypertriglyceridemia, hypercholesterolemia, elevated plasma apolipoprotein B and increased secretion of triglyceride-rich lipoproteins. The hyperlipidemia results from spontaneous mutation at a locus, Hyplip1, on distal mouse chromosome 3 in a region syntenic with a 1q21-q23 FCHL locus identified in Finnish, German, Chinese and US families. We fine-mapped Hyplip1 to roughly 160 kb, constructed a BAC contig and sequenced overlapping BACs to identify 13 candidate genes. We found substantially decreased mRNA expression for thioredoxin interacting protein (Txnip). Sequencing of the critical region revealed a Txnip nonsense mutation in HcB-19 that is absent in its normolipidemic parental strains. Txnip encodes a cytoplasmic protein that binds and inhibits thioredoxin, a major regulator of cellular redox state. The mutant mice have decreased CO2 production but increased ketone body synthesis, suggesting that altered redox status down-regulates the citric-acid cycle, sparing fatty acids for triglyceride and ketone body production. These results reveal a new pathway of potential clinical significance that contributes to plasma lipid metabolism.     

Delineation of a 50 kilobase DNA segment containing the recombination site in a sporadic case of Huntington's disease.

No detectable rearrangements involving chromosome 4p16.3 have been observed in patients with Huntington's disease (HD). New mutations for HD could involve structural alterations which might aid the localization of the defective gene. We have reinvestigated a well documented sporadic case of HD. DNA haplotyping with markers between D4S10 and the telomeric locus D4S141 reveals a recombination event in one chromosome of the sporadic HD patient. The site of recombination maps within a 50 kilobase (kb) region, about 700 kb from the 4p telomere. Based on the extremely low HD mutation rate and significantly decreased recombination in the distal region of 4p, we hypothesize a direct link between the site of the recombination and HD in this patient.     

The gene encoding 5-lipoxygenase activating protein confers risk of myocardial infarction and stroke

We mapped a gene predisposing to myocardial infarction to a locus on chromosome 13q12-13. A four-marker single-nucleotide polymorphism (SNP) haplotype in this locus spanning the gene ALOX5AP encoding 5-lipoxygenase activating protein (FLAP) is associated with a two times greater risk of myocardial infarction in Iceland. This haplotype also confers almost two times greater risk of stroke. Another ALOX5AP haplotype is associated with myocardial infarction in individuals from the UK. Stimulated neutrophils from individuals with myocardial infarction produce more leukotriene B4, a key product in the 5-lipoxygenase pathway, than do neutrophils from controls, and this difference is largely attributed to cells from males who carry the at-risk haplotype. We conclude that variants of ALOX5AP are involved in the pathogenesis of both myocardial infarction and stroke by increasing leukotriene production and inflammation in the arterial wall.     

Patterns of single-nucleotide polymorphisms in candidate genes for blood-pressure homeostasis.

Sequence variation in human genes is largely confined to single-nucleotide polymorphisms (SNPs) and is valuable in tests of association with common diseases and pharmacogenetic traits. We performed a systematic and comprehensive survey of molecular variation to assess the nature, pattern and frequency of SNPs in 75 candidate human genes for blood-pressure homeostasis and hypertension. We assayed 28 Mb (190 kb in 148 alleles) of genomic sequence, comprising the 5' and 3' untranslated regions (UTRs), introns and coding sequence of these genes, for sequence differences in individuals of African and Northern European descent using high-density variant detection arrays (VDAs). We identified 874 candidate human SNPs, of which 22% were confirmed by DNA sequencing to reveal a discordancy rate of 21% for VDA detection. The SNPs detected have an average minor allele frequency of 11%, and 387 are within the coding sequence (cSNPs). Of all cSNPs, 54% lead to a predicted change in the protein sequence, implying a high level of human protein diversity. These protein-altering SNPs are 38% of the total number of such SNPs expected, are more likely to be population-specific and are rarer in the human population, directly demonstrating the effects of natural selection on human genes. Overall, the degree of nucleotide polymorphism across these human genes, and orthologous great ape sequences, is highly variable and is correlated with the effects of functional conservation on gene sequences.     

Variation in factor B (BF) and complement component 2 (C2) genes is associated with age-related macular degeneration

Age-related macular degeneration (AMD) is the most common form of irreversible blindness in developed countries. Variants in the factor H gene (CFH, also known as HF1), which encodes a major inhibitor of the alternative complement pathway, are associated with the risk for developing AMD. Here we test the hypothesis that variation in genes encoding other regulatory proteins of the same pathway is associated with AMD. We screened factor B (BF) and complement component 2 (C2) genes, located in the major histocompatibility complex class III region, for genetic variation in two independent cohorts comprising approximately 900 individuals with AMD and approximately 400 matched controls. Haplotype analyses identify a statistically significant common risk haplotype (H1) and two protective haplotypes. The L9H variant of BF and the E318D variant of C2 (H10), as well as a variant in intron 10 of C2 and the R32Q variant of BF (H7), confer a significantly reduced risk of AMD (odds ratio = 0.45 and 0.36, respectively). Combined analysis of the C2 and BF haplotypes and CFH variants shows that variation in the two loci can predict the clinical outcome in 74% of the affected individuals and 56% of the controls. These data expand and refine our understanding of the genetic risk for AMD.     

The extent of linkage disequilibrium in Arabidopsis thaliana.

Linkage disequilibrium (LD), the nonrandom occurrence of alleles in haplotypes, has long been of interest to population geneticists. Recently, the rapidly increasing availability of genomic polymorphism data has fueled interest in LD as a tool for fine-scale mapping, in particular for human disease loci. The chromosomal extent of LD is crucial in this context, because it determines how dense a map must be for associations to be detected and, conversely, limits how finely loci may be mapped. Arabidopsis thaliana is expected to harbor unusually extensive LD because of its high degree of selfing. Several polymorphism studies have found very strong LD within individual loci, but also evidence of some recombination. Here we investigate the pattern of LD on a genomic scale and show that in global samples, LD decays within approximately 1 cM, or 250 kb. We also show that LD in local populations may be much stronger than that of global populations, presumably as a result of founder events. The combination of a relatively high level of polymorphism and extensive haplotype structure bodes well for developing a genome-wide LD map in A. thaliana.     

Variants modulating the expression of a chromosome domain encompassing PLAG1 influence bovine stature

We report mapping of a quantitative trait locus (QTL) with a major effect on bovine stature to a ～780-kb interval using a Hidden Markov Model-based approach that simultaneously exploits linkage and linkage disequilibrium. We re-sequenced the interval in six sires with known QTL genotype and identified 13 clustered candidate quantitative trait nucleotides (QTNs) out of >9,572 discovered variants. We eliminated five candidate QTNs by studying the phenotypic effect of a recombinant haplotype identified in a breed diversity panel. We show that the QTL influences fetal expression of seven of the nine genes mapping to the ～780-kb interval. We further show that two of the eight candidate QTNs, mapping to the PLAG1-CHCHD7 intergenic region, influence bidirectional promoter strength and affect binding of nuclear factors. By performing expression QTL analyses, we identified a splice site variant in CHCHD7 and exploited this naturally occurring null allele to exclude CHCHD7 as single causative gene.     

Prospects for whole-genome linkage disequilibrium mapping of common disease genes.

Recently, attention has focused on the use of whole-genome linkage disequilibrium (LD) studies to map common disease genes. Such studies would employ a dense map of single nucleotide polymorphisms (SNPs) to detect association between a marker and disease. Construction of SNP maps is currently underway. An essential issue yet to be settled is the required marker density of such maps. Here, I use population simulations to estimate the extent of LD surrounding common gene variants in the general human population as well as in isolated populations. Two main conclusions emerge from these investigations. First, a useful level of LD is unlikely to extend beyond an average distance of roughly 3 kb in the general population, which implies that approximately 500,000 SNPs will be required for whole-genome studies. Second, the extent of LD is similar in isolated populations unless the founding bottleneck is very narrow or the frequency of the variant is low (<5%).     

A common variant at the TERT-CLPTM1L locus is associated with estrogen receptor-negative breast cancer

Estrogen receptor (ER)-negative breast cancer shows a higher incidence in women of African ancestry compared to women of European ancestry. In search of common risk alleles for ER-negative breast cancer, we combined genome-wide association study (GWAS) data from women of African ancestry (1,004 ER-negative cases and 2,745 controls) and European ancestry (1,718 ER-negative cases and 3,670 controls), with replication testing conducted in an additional 2,292 ER-negative cases and 16,901 controls of European ancestry. We identified a common risk variant for ER-negative breast cancer at the TERT-CLPTM1L locus on chromosome 5p15 (rs10069690: per-allele odds ratio (OR) = 1.18 per allele, P = 1.0 × 10(-10)). The variant was also significantly associated with triple-negative (ER-negative, progesterone receptor (PR)-negative and human epidermal growth factor-2 (HER2)-negative) breast cancer (OR = 1.25, P = 1.1 × 10(-9)), particularly in younger women (<50 years of age) (OR = 1.48, P = 1.9 × 10(-9)). Our results identify a genetic locus associated with estrogen receptor negative breast cancer subtypes in multiple populations.