GOLD--graphical overview of linkage disequilibrium

Linkage disequilibrium mapping of schizophrenia susceptibility to the CAPON region of chromosome 1q22

Previously, we have reported linkage of markers from chromosome 1q22 to schizophrenia, a finding supported by several independent studies. We have now examined the region of strongest linkage for evidence of linkage disequilibrium (LD) in a sample of 24 Canadian familial-schizophrenia pedigrees. Analysis of 14 microsatellites and 15 single-nucleotide polymorphisms (SNPs) from the 5.4-Mb region between D1S1653 and D1S1677 produced significant evidence (nominal P<.05) of LD between schizophrenia and 2 microsatellites and 6 SNPs. All of the markers exhibiting significant LD to schizophrenia fall within the genomic extent of the gene for carboxyl-terminal PDZ ligand of neuronal nitric oxide synthase (CAPON), making it a prime positional candidate for the schizophrenia-susceptibility locus on 1q22, although initial mutation analysis of this gene has not identified any schizophrenia-associated changes within exons. Consistent with several recently identified candidate genes for schizophrenia, CAPON is involved in signal transduction in the NMDA receptor system, highlighting the potential importance of this pathway in the etiology of schizophrenia.

Association of neural cell adhesion molecule 1 gene polymorphisms with bipolar affective disorder in Japanese individuals

BACKGROUND

Although the pathogenesis of mood disorders remains unclear, heritable factors have been shown to be involved. Neural cell adhesion molecule 1 (NCAM1) is known to play important roles in cell migration, neurite growth, axonal guidance, and synaptic plasticity. Disturbance of these neurodevelopmental processes is proposed as one etiology for mood disorder. We therefore undertook genetic analysis of NCAM1 in mood disorders.

METHODS

We determined the complete genomic organization of human NCAM1 gene by comparing complementary deoxyribonucleic acid and genomic sequences; mutation screening detected 11 polymorphisms. The genotypic, allelic, and haplotype distributions of these variants were analyzed in unrelated control individuals (n = 357) and patients with bipolar disorder (n = 151) and unipolar disorder (n = 78), all from central Japan.

RESULTS

Three single nucleotide polymorphisms, IVS6+32T>C, IVS7+11G>C and IVS12+21C>A, displayed significant associations with bipolar disorder (for allelic associations, nominal p =.04, p =.02, and p =.004, respectively, all p >.05 after Bonferroni corrections). Furthermore, the haplotype located in a linkage disequilibrium block was strongly associated with bipolar disorder (the p value of the most significant three-marker haplotype is .005).

CONCLUSIONS

Our results suggest that genetic variations in NCAM1 or nearby genes could confer risks associated with bipolar affective disorder in Japanese individuals.

Haplotype block partitioning and tag SNP selection using genotype data and their applications to association studies

Recent studies have revealed that linkage disequilibrium (LD) patterns vary across the human genome with some regions of high LD interspersed by regions of low LD. A small fraction of SNPs (tag SNPs) is sufficient to capture most of the haplotype structure of the human genome. In this paper, we develop a method to partition haplotypes into blocks and to identify tag SNPs based on genotype data by combining a dynamic programming algorithm for haplotype block partitioning and tag SNP selection based on haplotype data with a variation of the expectation maximization (EM) algorithm for haplotype inference. We assess the effects of using either haplotype or genotype data in haplotype block identification and tag SNP selection as a function of several factors, including sample size, density or number of SNPs studied, allele frequencies, fraction of missing data, and genotyping error rate, using extensive simulations. We find that a modest number of haplotype or genotype samples will result in consistent block partitions and tag SNP selection. The power of association studies based on tag SNPs using genotype data is similar to that using haplotype data.

A simple correction for multiple testing for single-nucleotide polymorphisms in linkage disequilibrium with each other

In this report, we describe a simple correction for multiple testing of single-nucleotide polymorphisms (SNPs) in linkage disequilibrium (LD) with each other, on the basis of the spectral decomposition (SpD) of matrices of pairwise LD between SNPs. This method provides a useful alternative to more computationally intensive permutation tests. A user-friendly interface (SNPSpD) for performing this correction is available online (http://genepi.qimr.edu.au/general/daleN/SNPSpD/). Additionally, output from SNPSpD includes eigenvalues, principal-component coefficients, and factor "loadings" after varimax rotation, enabling the selection of a subset of SNPs that optimize the information in a genomic region.

Investigation of polymorphisms in the CREM gene in panic disorder

Clinical and animal studies suggest a role for pathways regulated by cyclic-AMP in anxiety. Mouse gene deletion studies, our own linkage findings on chromosome 10, and a recently published genetic association study by Domschke et al. [2003: Am J Med Genet 117B:70-78] suggest that the cAMP responsive element modulator (CREM) may be involved in panic disorder. We have employed a family-based design to investigate the role of DNA sequence variations in the gene for CREM in panic disorder. We have genotyped 613 individuals in 70 panic disorder pedigrees, as well as 42 parent/offspring triads. Subjects were genotyped at two informative single nucleotide polymorphisms (SNPs) and three polymorphic microsatellites in the CREM genomic region; and the data were analyzed for genetic association and linkage. Linkage analysis employing several diagnostic/genetic models produced a maximum lod score of 0.63 for SNP-1, located in the 5' UTR of CREM, under a dominant model with a broad diagnostic definition of panic disorder. Non-parametric analysis, using the NPL statistic or FBAT, also did not support any linkage or association between the markers and panic disorder. All five markers (spanning 77 kb) used in the study showed modest, but significant linkage disequilibrium. Analysis of 2-, 3-, 4-, or 5-marker haplotypes using TRANSMIT failed to find any globally significant results; however, individual haplotypes containing a single allele of MS-3 were nominally associated with panic disorder. These findings provide little additional evidence for a susceptibility locus for panic disorder either within the CREM gene or in a nearby region of chromosome 10p11 in our sample.

Genetic variation near the hepatocyte nuclear factor-4 alpha gene predicts susceptibility to type 2 diabetes

The Finland-United States Investigation Of NIDDM Genetics (FUSION) study aims to identify genetic variants that predispose to type 2 diabetes by studying affected sibling pair families from Finland. Chromosome 20 showed our strongest initial evidence for linkage. It currently has a maximum logarithm of odds (LOD) score of 2.48 at 70 cM in a set of 495 families. In this study, we searched for diabetes susceptibility variant(s) at 20q13 by genotyping single nucleotide polymorphism (SNP) markers in case and control DNA pools. Of 291 SNPs successfully typed in a 7.5-Mb interval, the strongest association confirmed by individual genotyping was with SNP rs2144908, located 1.3 kb downstream of the primary beta-cell promoter P2 of hepatocyte nuclear factor-4 alpha (HNF4A). This SNP showed association with diabetes disease status (odds ratio [OR] 1.33, 95% CI 1.06-1.65, P = 0.011) and with several diabetes-related traits. Most of the evidence for linkage at 20q13 could be attributed to the families carrying the risk allele. We subsequently found nine additional associated SNPs spanning a 64-kb region, including the P2 and P1 promoters and exons 1-3. Our results and the independent observation of association of SNPs near the P2 promoter with diabetes in a separate study population of Ashkenazi Jewish origin suggests that variant(s) located near or within HNF4A increases susceptibility to type 2 diabetes.

Variations in GABRA2, encoding the alpha 2 subunit of the GABA(A) receptor, are associated with alcohol dependence and with brain oscillations

Alcoholism is a complex disease with both genetic and environmental risk factors. To identify genes that affect the risk for alcoholism, we systematically ascertained and carefully assessed individuals in families with multiple alcoholics. Linkage and association analyses suggested that a region of chromosome 4p contained genes affecting a quantitative endophenotype, brain oscillations in the beta frequency range (13-28 Hz), and the risk for alcoholism. To identify the individual genes that affect these phenotypes, we performed linkage disequilibrium analyses of 69 single-nucleotide polymorphism (SNPs) within a cluster of four GABA(A) receptor genes, GABRG1, GABRA2, GABRA4, and GABRB1, at the center of the linked region. GABA(A) receptors mediate important effects of alcohol and also modulate beta frequencies. Thirty-one SNPs in GABRA2, but only 1 of the 20 SNPs in the flanking genes, showed significant association with alcoholism. Twenty-five of the GABRA2 SNPs, but only one of the SNPs in the flanking genes, were associated with the brain oscillations in the beta frequency. The region of strongest association with alcohol dependence extended from intron 3 past the 3' end of GABRA2; all 43 of the consecutive three-SNP haplotypes in this region of GABRA2 were highly significant. A three-SNP haplotype was associated with alcoholism, with P=.000000022. No coding differences were found between the high-risk and low-risk haplotypes, suggesting that the effect is mediated through gene regulation. The very strong association of GABRA2 with both alcohol dependence and the beta frequency of the electroencephalogram, combined with biological evidence for a role of this gene in both phenotypes, suggest that GABRA2 might influence susceptibility to alcohol dependence by modulating the level of neural excitation.

Extended haplotype analysis in the HLA complex reveals an increased frequency of the HFE-C282Y mutation in individuals with multiple sclerosis

In order to resolve a multiple sclerosis (MS) susceptibility locus that we had identified in earlier work at the telomeric end of the HLA complex, we genotyped another 34 microsatellite markers (47 in total) across the class I/extended class I region in 166 Tasmanian MS case and 104 control families (D6S299-D6S265). Extended MS susceptibility haplotypes, up to 9 Mb in length, were observed in 11% of MS cases and 4% of controls. Direct comparison of the telomerically extended portion of the MS susceptibility haplotype in HFE-Cys282Tyr (C282Y)-homozygous haemochromatosis patients identified a common ancestry for this genomic segment, which translated into an increased frequency of the C282Y allele in 489 MS cases from Tasmania and Victoria (10.2%) compared with controls (6.7%). Six C282Y homozygotes (1.2%), a three-fold increased rate over the general population, and 88 heterozygotes (18%) were identified. One C282Y-homozygous female was identified who had MS and was being treated for symptoms of iron overload. Interestingly, for 71 Victorian MS cases not of north western European (NWE) ancestry, a DR15-independent reduction in the frequency of the C282Y allele was observed, supporting the theory of a NWE origin for the C282Y-variant of the DR15 ancestral haplotype (C282Y-HLA-A*0301-B*0702-DRB1*1501-DQB1*0602). The results of linkage disequilibrium (LD) and log linear modelling analyses suggest that C282Y is increased in MS cases of NWE ancestry because it is in LD with the ancestral DR15 susceptibility haplotype (7.1) and that it does not play an independent role in predisposition to MS. However, our findings provide the impetus for further investigations into the role of iron metabolism in the severity of MS.

Genetic variation in COL17A1 and the development of bullous pemphigoid

BACKGROUND

Bullous pemphigoid (BP) is an autoimmune blistering disease of the skin characterized by autoantibody attack on collagen XVII.

OBJECTIVES

To characterize the genetic complexity of COL17A1, the gene which encodes for the autoantigen collagen XVII. The data will be used to determine whether there is an association between polymorphisms and haplotypes of COL17A1 and genetic susceptibility to development of BP.

METHODS

The genetic complexity in COL17A1 was deduced by screening and then sequencing the gene. Haplotypes were constructed from the resulting polymorphisms using the statistical programme PHASE. The linkage disequilibrium (D') between the polymorphisms was deduced from haplotypic data using the statistical programme GOLD. Association of the polymorphisms and haplotypes was tested for, in a cohort of BP patients and controls.

RESULTS

Screening of COL17A1 for genetic variation was carried out in 29 individuals of North European caucasoid origin, and it revealed 19 single-nucleotide polymorphisms in approximately 14.7 kb of sequence. These variants resulted in 60 different haplotypes in 191 individuals, of which 13 occurred above 1% in the population. D' between the variants was found to be extensive, have a low correlation with physical distance and to extend over 33.8 kb. No association was found with any of the polymorphisms or haplotypes and development of BP, when tested for, in a cohort of patients and controls.

CONCLUSION

This study provides an extensive description of the genetic variation in COL17A1 and shows no association of the genetic variants with susceptibility to BP.

Confirmation and fine mapping of the chromosome 1 alcohol dependence risk locus

Two previous large genetic linkage studies in the US population have implicated an area in chromosome 1p to contain a susceptibility gene for alcohol dependence. The 1-LOD support interval of the linkage signal spans about 30 cM and contains >30000000 DNA base pairs (bp) and 700 predicted genes. In order to reduce the size of the candidate area and potentially identify novel candidate genes within this region, we fine-mapped this area using closely spaced short tandem repeat (STR) markers and the transmission disequilibrium test (TDT) in small nuclear families. The subjects were 87 European-American families including one or more alcohol-dependent offspring (93 children and 174 parents). The initial marker set consisted of 30 STR markers, spanning the Marshfield map interval between 101.48 and 130.73 cM. Using the TDTPHASE program, we identified three markers in the distal part of this region (125-126 cM), which showed evidence of transmission disequilibrium. On the basis of this result, an additional 12 STR markers were genotyped in this region; some of these markers provided additional evidence for linkage disequilibrium. The strongest evidence for transmission disequilibrium was obtained at the marker D1S406 (P=0.005, 126.16 cM), with supporting evidence from three neighboring STR markers D1S424 (126.16 cM, P=0.01), D1S2804 (126.16 cM, P=0.04), and D1S2776 (126.16 cM, P=0.02), which are all located within a <350000 bp interval. These findings suggest that a gene (or genes) causing susceptibility to alcohol dependence resides near location 126.16 cM on chromosome 1. In addition, these results provide independent confirmation of the linkage finding regarding the identification of at least one gene in this region increasing the risk for alcohol dependence.

Evidence for genetic linkage between a polymorphism in the adenosine 2A receptor and panic disorder

Data from clinical and behavioral pharmacological studies have implicated adenosine in anxiety behaviors, while genetic studies have suggested that adenosine receptors may be associated with panic disorder. We have undertaken an analysis of several DNA sequence variations in the adenosine 2A receptor (ADORA2A) in a large sample of panic disorder pedigrees. Individuals from 70 panic disorder pedigrees, and 83 child-parent 'trios', were genotyped at five single-nucleotide polymorphisms (SNPs) in and near the ADORA2A gene and were analyzed for genetic linkage and association. Linkage analysis revealed elevated LOD scores for a silent substitution (1083C/T, SNP-4) in the second coding exon. This SNP has been previously reported to be associated with panic disorder. We observed a maximal heterogeneity LOD score of 2.98 (theta=0) under a recessive genetic model and narrow diagnostic model. Other SNPs showed no evidence for linkage. Association tests were not significant for any of the five ADORA2A SNPs. When SNP haplotypes were assessed in the triads with TRANSMIT, one 3-marker haplotype (SNPs 1, 4, 5) was nominally significantly associated with panic disorder (p=0.029). Pairwise estimations of linkage disequilibrium between the SNPs showed strong patterns of linkage disequilibrium across the ADORA2A locus. Analyses carried out by broadening the panic disorder phenotype to include agoraphobia continued to support linkage to ADORA2A. Our findings provide evidence for a susceptibility locus for panic disorder, and possibly including agoraphobia, either within the ADORA2A gene or in a nearby region of chromosome 22, and serves as the first successful candidate gene replication study in panic disorder.

No association between the APOE gene and autism

Autism is a neurodevelopmental disorder characterized by stereotypic and repetitive behavior and interests, together with social and communicative deficiencies. The results of several genomic screens suggest the presence of an autism susceptibility locus on chromosome 19p13.2-q13.4. The apolipoprotein E (APOE) gene on chromosome 19 encodes for a protein, apoE, whose different isoforms (E2, E3, E4) influence neuronal growth. APOE participates in lipid transport and metabolism, repair, growth, and maintenance of axons and myelin during neuronal development. The APOE protein competes with the Reelin protein for VLDL/APOER2 receptor binding. Several studies have reported evidence for an association between autism and the Reelin gene. Based on these data we tested for association between APOE and autism using family-based association methods in a data set of 322 autism families. Three promoter, one intronic, and one 3' UTR single nucleotide polymorphisms (SNPs) in the APOE gene (-491a/t, -427c/t, -219g/t, 113c/g, and 5361c/t) as well as the APOE functional polymorphism (E2, E3, E4) were examined and failed to reveal significant evidence that autism is associated with APOE.

Association of tryptophan 2,3 dioxygenase gene polymorphism with autism

Although elevation of blood and platelet serotonin has been documented in autism, genetic analyses of serotonin transporter gene have given conflicting results. Tryptophan 2,3 dioxygenase (TDO2) is the rate-limiting enzyme in the catabolism of tryptophan, the precursor of serotonin. A mutation that results in decreased activity of the TDO2 can decrease catabolism of tryptophan and increase the level of whole body serotonin. As such it is a potential candidate gene for autism. We have investigated five single nucleotide polymorphisms in the TDO2 gene for association with autistic disorder. One hundred and ninety six multiplex autistic disorder families were tested using transmission disequilibrium test. There was a significant difference in the transmission of a promoter variant to autistic subjects (P = 0.0006). Haplotype analysis also demonstrated significant difference in the transmission of TDO2 haplotypes to autistic subjects (P = 0.0027). Our results suggest the presence of a susceptibility mutation in the TDO2 or a nearby gene, but may also represent a chance finding.

Haplotype mapping of the bronchiolitis susceptibility locus near IL8

Susceptibility to viral bronchiolitis, the commonest cause of infant admissions to hospital in the industrialised world, is associated with polymorphism at the IL8 locus. Here we map the genomic boundaries of the disease association by case-control analysis and TDT in 580 affected UK infants. Markers for association mapping were chosen after determining patterns of linkage disequilibrium across the surrounding region of chromosome 4q, a 550-kb segment containing nine genes, extending from AFP to PPBP. The region has three major clusters of high linkage disequilibrium and is notable for its low haplotypic diversity. We exclude adjacent chemokine genes as the cause of the association, and identify a disease-associated haplotype that spans a 250-kb region from AFM to IL8. In between these two genes there is only one structural feature of interest, a novel gene RASSF6, which is predicted to encode a Ras effector protein.

Both rare and common polymorphisms contribute functional variation at CHGA, a regulator of catecholamine physiology

The chromogranin/secretogranin proteins are costored and coreleased with catecholamines from secretory vesicles in chromaffin cells and noradrenergic neurons. Chromogranin A (CHGA) regulates catecholamine storage and release through intracellular (vesiculogenic) and extracellular (catecholamine release-inhibitory) mechanisms. CHGA is a candidate gene for autonomic dysfunction syndromes, including intermediate phenotypes that contribute to human hypertension. Here, we show a surprising pattern of CHGA variants that alter the expression and function of this gene, both in vivo and in vitro. Functional variants include both common alleles that quantitatively alter gene expression and rare alleles that qualitatively change the encoded product to alter the signaling potency of CHGA-derived catecholamine release-inhibitory catestatin peptides.

Susceptibility to leprosy is associated with PARK2 and PACRG

Leprosy is caused by Mycobacterium leprae and affects about 700,000 individuals each year. It has long been thought that leprosy has a strong genetic component, and recently we mapped a leprosy susceptibility locus to chromosome 6 region q25-q26 (ref. 3). Here we investigate this region further by using a systematic association scan of the chromosomal interval most likely to harbour this leprosy susceptibility locus. In 197 Vietnamese families we found a significant association between leprosy and 17 markers located in a block of approx. 80 kilobases overlapping the 5' regulatory region shared by the Parkinson's disease gene PARK2 and the co-regulated gene PACRG. Possession of as few as two of the 17 risk alleles was highly predictive of leprosy. This was confirmed in a sample of 975 unrelated leprosy cases and controls from Brazil in whom the same alleles were strongly associated with leprosy. Variants in the regulatory region shared by PARK2 and PACRG therefore act as common risk factors for leprosy.

Evidence for a cluster of genes on chromosome 17q11–q21 controlling susceptibility to tuberculosis and leprosy in Brazilians

The region of conserved synteny on mouse chromosome 11/human 17q11-q21 is known to carry a susceptibility gene(s) for intramacrophage pathogens. The region is rich in candidates including NOS2A, CCL2/MCP-1, CCL3/MIP-1alpha, CCL4/MIP-1beta, CCL5/RANTES, CCR7, STAT3 and STAT5A/5B. To examine the region in man, we studied 92 multicase tuberculosis (627 individuals) and 72 multicase leprosy (372 individuals) families from Brazil. Multipoint nonparametric analysis (ALLEGRO) using 16 microsatellites shows two peaks of linkage for leprosy at D17S250 (Z(lr) score 2.34; P=0.01) and D17S1795 (Z(lr) 2.67; P=0.004) and a single peak for tuberculosis at D17S250 (Z(lr) 2.04; P=0.02). Combined analysis shows significant linkage (peak Z(lr) 3.38) at D17S250, equivalent to an allele sharing LOD score 2.48 (P=0.0004). To determine whether one or multiple genes contribute, 49 informative single nucleotide polymorphisms were typed in candidate genes. Family-based allelic association testing that was robust to family clustering demonstrated significant associations with tuberculosis susceptibility at four loci separated by intervals (NOS2A-8.4 Mb-CCL18-32.3 kb-CCL4-6.04 Mb-STAT5B) up to several Mb. Stepwise conditional logistic regression analysis using a case/pseudo-control data set showed that the four genes contributed separate main effects, consistent with a cluster of susceptibility genes across 17q11.2.

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.

Positional identification of hypertension susceptibility genes on chromosome 2

Chromosome 2 has been consistently identified as a genomic region with genetic linkage evidence suggesting that one or more loci contributes to blood pressure and hypertension status. As with all complex disease traits, following-up linkage evidence to identify the underlying susceptibility gene(s) is an arduous yet biologically and clinically important task. Using combined positional candidate gene methods, the Family Blood Pressure Program (FBPP) has concentrated efforts in narrowing a large region of chromosome 2, demonstrating evidence for linkage in several populations, and identifying underlying candidate hypertension susceptibility gene(s). Initial informatics efforts identified the boundaries of the region and the known genes within it. A total of 82 polymorphic sites in 8 genes were genotyped in a large hypothesis-generating sample consisting of 1640 African Americans, 1339 whites, and 1616 Mexican Americans. After resampling-based false discovery adjustment, SLC4A5, a sodium bicarbonate transporter, was identified as a primary candidate gene for hypertension. Polymorphisms in SLC4A5 were subsequently genotyped and analyzed for validation in two other subcomponents of the FBPP, each contributing African Americans (N=461; N=778) and whites (N=550; N=967). Again, single nucleotide polymorphisms within this gene were significantly associated with blood pressure levels and hypertension status. Although not identifying a single causal gene variant that is significantly associated with blood pressure levels and hypertension status across all samples, the results further implicate SLC4A5 as a candidate hypertension susceptibility gene. Moreover, the present study validates previous evidence for one or more genes on chromosome 2 that influence hypertension-related phenotypes in the population-at-large.

Endogenous Opioid Receptor Genes and Alcohol Dependence Among Taiwanese Han

BACKGROUND

Nonselective opioid antagonists reduce alcohol consumption under various experimental situations, and several association studies have examined possible roles of opioid receptor mu (OPRM), delta (OPRD), and kappa (OPRK) genes in the development of alcohol dependence.

METHODS

We examined 20 single nucleotide polymorphisms (SNPs) across the OPRM, OPRD, and OPRK genes in 158 alcohol-dependent subjects and 149 controls. Differences in allele frequency and genotype distribution between case subjects and controls, as well as the deviation from Hardy-Weinberg equilibrium, were examined using Fisher's exact tests.

RESULTS

No significant difference in either allele or genotype frequency was found between case subjects and controls for each of the SNPs.

CONCLUSIONS

Our findings do not support a possible role of the opioid receptor genes for the proclivity to alcohol dependence in the Taiwanese Han.

Examination of candidate genes in language disorder: A model of genetic association for treatment studies

The purpose of this review is to provide a model for studying genetic association of response to intervention in child language disorders. In addition to a theoretical overview and review of different approaches to studying candidate genes, a specific methodology for completing this type of analysis is presented. The goal of the analysis is to provide detail beyond simple broad phenotyping for affected and nonaffected individuals and to take advantage of data yielded from concise behavioral phenotyping often available in treatment studies.

Findings in an independent sample support an association between bipolar affective disorder and the G72/G30 locus on chromosome 13q33

Markers near the nested genes G72 and G30 on chromosome 13q33 have been implicated in the etiology of schizophrenia and, recently, bipolar affective disorder (BPAD). Hattori et al (2003) reported that single-nucleotide polymorphisms (SNPs) near the G72/G30 locus were associated with BPAD in a sample of 22 pedigrees, and that SNP haplotypes were associated in a second, larger sample of triads. The present study attempts to replicate this finding in an independent case-control sample. Six SNPs near the G72/G30 locus, including the most strongly associated markers in the previous study, were tested in 139 cases and 113 ethnically matched controls. Significant association was detected between BPAD and two adjacent SNPs (smallest P=0.007; global P=0.024). Haplotype analysis produced additional support for association (smallest P=0.004; global P=0.004). Analysis of 31 unlinked microsatellite markers detected no population stratification in the cases or controls studied. Although the associated alleles and haplotypes differ from those previously reported, these new results provide further evidence, in an independent sample, for an association between BPAD and genetic variation in the vicinity of the genes G72 and G30.

Association of GABRG3 With Alcohol Dependence

BACKGROUND

Evidence from human, animal, and in vitro cell models suggests that gamma-aminobutyric acid (GABA), the major inhibitory neurotransmitter in the human central nervous system, is involved in many of the neurochemical pathways that affect alcohol use, abuse, and dependence. Both linkage and association to the region on chromosome 15q that contains a cluster of GABAA receptor genes have previously been reported, but the role of these genes in alcoholism remains inconclusive.

METHODS

We conducted family-based association analyses by using a large sample of multiplex alcoholic families collected as part of the Collaborative Study on the Genetics of Alcoholism, to test for an association between alcohol dependence and the GABAA receptor genes clustered on chromosome 15q. Multiple single-nucleotide polymorphisms were tested in each of the three chromosome 15q GABAA receptor genes: GABRA5, GABRB3, and GABRG3.

RESULTS

Using both classic trio-based analyses and extended-family analyses, we found consistent evidence of association between alcohol dependence and GABRG3. Nearly all single-nucleotide polymorphisms across the gene yielded evidence of association, and haplotype analyses were highly significant. No consistent evidence of association was observed with either GABRA5 or GABRB3, nor was there evidence for parent-of-origin effects with any of the genes.

CONCLUSIONS

These analyses suggest that GABRG3 may be involved in the risk for alcohol dependence. These findings support the theory that the predisposition to alcoholism may be inherited as a general state of central nervous system disinhibition/hyperexcitability that results from an altered responsiveness to GABA.

Association of T-786C eNOS gene polymorphism with increased susceptibility to acute chest syndrome in females with sickle cell disease

Acute chest syndrome (ACS) is a life-threatening complication of sickle cell disease (SCD). A retrospective study was performed to evaluate the role of endothelial nitric oxide synthase (eNOS) gene polymorphisms (E298D and T-786C) in African-American SCD patients. The D298 allele showed no association; the C-786 allele showed a statistically significant association (P = 0.0061) in female ACS cases. Multiple logistic regression analysis showed that relative risk of ACS was 8.695 (P = 0.0076, 95% confidence interval 1.761-42.920) for female carriers of C-786. eNOS T-786C is a gender-specific genetic modifier that is associated with increased susceptibility to ACS in female SCD patients.

Defining haplotype blocks and tag single-nucleotide polymorphisms in the human genome

Recent studies suggest that the genome is organized into blocks of haplotypes, and efforts to create a genome-wide haplotype map of single-nucleotide polymorphisms (SNPs) are already underway. Haplotype blocks are defined algorithmically and to date several algorithms have been proposed. However, little is known about their relative performance in real data or about the impact of allele frequencies and parameter choices on the detection of haplotype blocks and the markers that tag them. Here we present a formal comparison of two major algorithms, a linkage disequilibrium (LD)-based method and a dynamic programming algorithm (DPA), in three chromosomal regions differing in gene content and recombination rate. The two methods produced strikingly different results. DPA identified fewer and larger haplotype blocks as well as a smaller set of tag SNPs than the LD method. For both methods, the results were strongly dependent on the allele frequency. Decreasing the minor allele frequency led to an up to 3.7-fold increase in the number of haplotype blocks and tag SNPs. Definition of haploytpe blocks and tag SNPs was also sensitive to parameter changes, but the results could not be reconciled simply by parameter adjustment. These results show that two major methods for detecting haplotype blocks and tag SNPs can produce different results in the same data and that these results are sensitive to marker allele frequencies and parameter choices. More information is needed to guide the choice of method, marker allele frequencies, and parameters in the development of a haplotype map.

Screening of nine candidate genes for autism on chromosome 2q reveals rare nonsynonymous variants in the cAMP-GEFII gene

The results from several genome scans indicate that chromosome 2q21-q33 is likely to contain an autism susceptibility locus. We studied the potential contribution of nine positional and functional candidate genes: TBR-1; GAD1; DLX1; DLX2; cAMP-GEFII; CHN1; ATF2; HOXD1 and NEUROD1. Screening these genes for DNA variants and association analysis using intragenic single nucleotide polymorphisms did not provide evidence for a major role in the aetiology of autism. Four rare nonsynonymous variants were identified, however, in the cAMP-GEFII gene. These variants were present in five families, where they segregate with the autistic phenotype, and were not observed in control individuals. The significance of these variants is unclear, as their low frequency in IMGSAC families does not account for the relatively strong linkage signal at the 2q locus. Further studies are needed to clarify the contribution of cAMP-GEFII gene variants to autism susceptibility.

Polymorphisms in the 5'-untranslated region of the human serotonin receptor 1B (HTR1B) gene affect gene expression

We present evidence of complex balancing regulation of HTR1B transcription by common polymorphisms in its promoter. Computational analysis of the HTR1B gene predicted that a 5' segment, spanning common DNA sequence variations, T-261G, A-161T, and -182INS/DEL-181, contained a putative functional promoter. Using a secreted alkaline phosphatase (SEAP) reporter gene system, we found that the haplotype -261G_-182INS-181_A-161 enhanced transcriptional activity 2.3-fold compared with the haplotype T-261_-182INS-181_A-161. Conversely, -161T reversed this, and the net effect when -261G and -161T were in the same haplotype (-261G_-182INS-181_-161T) was equivalent to the major haplotype (T-261_-182INS-181_A-161). Electrophoretic mobility shift experiments showed that -261G and -161T modify the binding of transcription factors (TFs): -261G generates a new AP2 binding site, while alleles A-161 and -161T exhibit different binding characteristics to AP1. T-261G and A-161T were found to be in linkage disequilibrium (LD) with G861C in a European ancestry population. Interestingly, G861C has been reported to be associated with several psychiatric disorders. Our results indicate that HTR1B is the target of substantial transcriptional genetic regulation by common haplotypes, which are in LD with the HTR1B single-nucleotide polymorphism (SNP) most commonly used in association studies.

Positional cloning of a novel gene influencing asthma from chromosome 2q14

Asthma is a common disease in children and young adults. Four separate reports have linked asthma and related phenotypes to an ill-defined interval between 2q14 and 2q32 (refs. 1-4), and two mouse genome screens have linked bronchial hyper-responsiveness to the region homologous to 2q14 (refs. 5,6). We found and replicated association between asthma and the D2S308 microsatellite, 800 kb distal to the IL1 cluster on 2q14. We sequenced the surrounding region and constructed a comprehensive, high-density, single-nucleotide polymorphism (SNP) linkage disequilibrium (LD) map. SNP association was limited to the initial exons of a solitary gene of 3.6 kb (DPP10), which extends over 1 Mb of genomic DNA. DPP10 encodes a homolog of dipeptidyl peptidases (DPPs) that cleave terminal dipeptides from cytokines and chemokines, and it presents a potential new target for asthma therapy.

LD mapping of maternally and non-maternally derived alleles and atopy in FcepsilonRI-beta

Polymorphisms in the beta chain of the high affinity receptor for IgE (Fc epsilon RI-beta, MS4A2) are consistently associated with traits underlying asthma and atopy (immunoglobulin E-mediated allergy). However, the causal variants and haplotypes underlying disease have not yet been identified. Maternal effects, with association confined to maternally derived alleles, have been shown in some studies but not in others. We have therefore extended the known sequence and systematically detected polymorphisms across an 18.1 Kb genomic region that includes Fc epsilon RI-beta. Association testing in two panels of subjects showed the presence of single-nucleotide polymorphisms (SNPs) affecting prick skin tests and specific IgE responses in several clusters. Stepwise analyses indicated that the clusters represent independent effects. Interferon regulatory factor 2 (IRF-2) sites were altered by significantly associated SNPs in two regions. Strong association to maternally derived alleles was seen in one panel of subjects and not in the other. Maternal and non-maternally derived associations tended to share the same SNP clusters, but associations were stronger in the presence of maternal effects. Two regions of increased CpG concentration were identified in Fc epsilon RI-beta. One of these approximated a SNP cluster that showed strong association and maternal effects, providing a potential substrate for epigenetic effects.

Fine mapping of the alpha-T catenin gene to a quantitative trait locus on chromosome 10 in late-onset Alzheimer's disease pedigrees

Using plasma amyloid beta protein (Abeta42) levels as an intermediate, quantitative phenotype for late onset Alzheimer's disease (LOAD), we previously obtained significant linkage at approximately 80 cM on chromosome 10. Linkage to the same region was obtained independently in a study of affected LOAD sib-pairs. Together, these two studies provide strong evidence for a novel LOAD locus on chromosome 10 that acts to increase Abeta42. VR22 is a large (1.7 Mb) gene located at 80 cM that encodes alpha-T catenin, which is a binding partner of beta catenin. This makes VR22 an attractive candidate gene because beta catenin interacts with presenilin 1, which has many mutations that elevate Abeta42 and cause early onset familial AD. We identified two intronic VR22 SNPs (4360 and 4783) in strong linkage disequilibrium (LD) that showed highly significant association (P=0.0001 and 0.0006) with plasma Abeta42 in 10 extended LOAD families. This association clearly contributed to the linkage at approximately 80 cM because the lod scores decreased when linkage analysis was performed conditional upon the VR22 association. This association replicated in another independent set of 12 LOAD families (P=0.04 for 4783 and P=0.08 for 4360). Bounding of the association region using multiple SNPs showed VR22 to be the only confirmed gene within the region of association. These findings indicate that VR22 has variant(s) which influence Abeta42 and contribute to the previously reported linkage for plasma Abeta42 in LOAD families.

No evidence of type 1 diabetes susceptibility genes in the region centromeric of the HLA complex

There is strong evidence that DQB1, DQA1, and DRB1 alleles are not the only contributors to the human leukocyte antigen (HLA) linked type 1 diabetes (T1D) predisposition. Although the HLA complex is much studied for disease association, little is known about the neighboring centromeric region. We have previously found suggestive association on DQ2-DR3 haplotypes for marker D6S291, located 3.6-Mb centromeric of HLA-DQB1. This region on human chromosome 6 is syntenic to a part of the region adjacent to the mouse major histocompatibility complex (MHC) on chromosome 17, which has been suggested to harbor a susceptibility gene in mouse (Idd16). To evaluate a possible role of the region centromeric of HLA-DQB1 in human T1D, we have scanned the region with nine microsatellite markers in 267 T1D families from five different populations. Our results indicate that the characteristic strong linkage disequilibrium in the HLA complex does not extend into this region. Furthermore, we did not detect any consistent T1D association for the markers analyzed in the study. In conclusion, our data argue against the presence of any strong genetic susceptibility factors for T1D in the region centromeric of the HLA complex.

Variants of the CD40 ligand gene are not associated with increased susceptibility to tuberculosis in West Africa

Evidence for linkage between tuberculosis and human chromosomal region Xq26 has previously been described. The costimulatory molecule CD40 ligand, encoded by TNFSF5 and located at Xq26.3, is a promising positional candidate. Interactions between CD40 ligand and CD40 are involved in the development of humoral- and cell-mediated immunity, as well as the activation of macrophages, which are the primary host and effector cells for Mycobacterium tuberculosis. We hypothesised that common variation within TNFSF5 might affect susceptibility to tuberculosis disease and, thus, might be responsible for the observed linkage to Xq26. Sequencing 32 chromosomes from a Gambian population identified nine common polymorphisms within the coding, 3' and 5' regulatory sequences of the gene. Six single nucleotide polymorphisms (SNPs) and a 3' microsatellite were genotyped in 121 tuberculosis patients and their available parents. No association with tuberculosis was detected for these variants using a transmission disequilibrium test, although one SNP at -726 showed some evidence of association in males. This finding, however, did not replicate in a separate case control study of over 1,200 West African individuals. We conclude that common genetic variation in TNFSF5 is not likely to affect tuberculosis susceptibility in West Africa and the linkage observed in this region is not due to variation in TNFSF5.

Analysis of reelin as a candidate gene for autism

Genetic studies indicate that chromosome 7q is likely to contain an autism susceptibility locus (AUTS1). We have followed a positional candidate gene approach to identify relevant gene(s) and report here the analysis of reelin (RELN), a gene located under our peak of linkage. Screening RELN for DNA changes identified novel missense variants absent in a large control group; however, the low frequency of these mutations does not explain the relatively strong linkage results on 7q. Furthermore, analysis of a previously reported triplet repeat polymorphism and intragenic single nucleotide polymorphisms, using the transmission disequilibrium test, provided no evidence for association with autism in IMGSAC and German singleton families. The analysis of RELN suggests that it probably does not play a major role in autism aetiology, although further analysis of several missense mutations is warranted in additional affected individuals.

Titin is a candidate gene for stroke volume response to endurance training: the HERITAGE Family Study

A genome-wide linkage scan for endurance training-induced changes in submaximal exercise stroke volume (DeltaSV50) in the HERITAGE Family Study revealed two chromosomal regions (2q31-q32 and 10p11.2) with at least suggestive evidence of linkage among white families. Here we report a further characterization of the quantitative trait locus (QTL) in chromosome 2q31 and provide evidence that titin (TTN) is likely a candidate gene involved. The original linkage was detected with two markers (D2S335 and D2S1391), and the QTL covered approximately 25 million base pairs (Mb). We added 12 microsatellite markers resulting in an average marker density of one marker per 2.3 Mb. The evidence of linkage increased from P = 0.006 to P = 0.0002 and 0.00002 in the multi- and single-point analyses, respectively. The strongest evidence of linkage was seen with two markers in and near the TTN gene. Transmission/disequilibrium test (TDT) with the same marker set provided evidence for association with one of the TTN markers (D2S385; P = 0.004). TTN is a major contributor to the elasticity of cardiomyocytes and a key regulator of the Frank-Starling mechanism. Since TTN is the largest gene in the human genome, the challenge is to identify the DNA sequence variants contributing to the interindividual differences in cardiac adaptation to endurance training.

Haplotype information and linkage disequilibrium mapping for single nucleotide polymorphisms

Single nucleotide polymorphisms in the human genome have become an increasingly popular topic in that their analyses promise to be a key step toward personalized medicine. We investigate two related questions, how much the haplotype information contributes to linkage disequilibrium (LD) mapping and whether an in silico haplotype construction preceding the LD analysis can help. For disease gene mapping, using both simulated and real data sets on cystic fibrosis and the Alzheimer disease, we reached the following conclusions: (1) for simple Mendelian diseases, in which case a tractable full statistical model can be developed, the loss of haplotype information for either control or disease data do not have a great impact on LD fine mapping, and haplotype inference should be carried out jointly with LD mapping; (2) for complex diseases, inferring haplotype phases for individuals prior to LD mapping helps achieve a better accuracy. An improved version of the linkage disequilibrium mapping program, BLADE v2, is available at http://www.fas.harvard.edu/junliu/TechRept/03folder/bladev2.tgz.

Genotype frequencies and linkage disequilibrium in the CEPH human diversity panel for variants in folate pathway genesMTHFR,MTHFD,MTRR,RFC1, andGCP2

BACKGROUND

Genetic variation in enzymes involved in vitamin metabolism is a candidate for analysis in studies of how nutritional covariates may impact a disease state. The role of folate pathway genes in birth defects and cardiovascular disease in humans has been widely studied. Since incidence rates for these disorders vary by geographic origins, it is useful to know which variants are the best candidates for studies based on genotype and allele frequency, as well as linkage disequilibrium (LD) in founder populations.

METHODS

Six polymorphisms in five folate metabolism-related genes (MTHFR, MTHFD, MTRR, GCP2, and RFC1) were genotyped on a collection of 1064 DNA samples from populations around the world, which were made available by the Centre d'Etude du Polymorphisme Humain (CEPH) consortium for analysis.

RESULTS

In this study we report the genotype frequencies for variants in the MTHFR, MTHFD, MTRR, GCP2, and RFC1 genes, and the LD for two variants (C677T and A1298C) in MTHFR.

CONCLUSIONS

The rare allele frequency for each of the five genes studied varied widely. LD is strongest in Pakistani and Brazilian populations (D' = 1.0) and weakest in Mexican populations (D' = 0.45). These findings will allow the selection of variants that will provide the most power in studies of folate pathway genes involving different ancestral populations, and contribute to our knowledge of the population distribution of selected nutritional gene variants.

Linkage and association between serotonin 2A receptor gene polymorphisms and bipolar I disorder

Several inconsistent associations between bipolar I disorder (BD1) and polymorphisms of the genes encoding the serotonin 2A receptor (HTR2A) have been published. We conducted the Transmission Disequilibrium Test (TDT) and case-control comparisons involving nine single nucleotide polymorphisms at the serotonin 2A receptor gene (four SNPs of HTR2A exons and five flanking SNPs). Comparison of BD1 cases (n = 93) with a group of unrelated population based controls (n = 92) revealed associations with SNPs on exons 2 and 3 (516C/T and 1354C/T, respectively), consistent with haplotype-based differences. Analysis of the cases and their available parents using the TDT suggested significant linkage and associations with 1354C/T, as well as haplotypes bearing this SNP. Our results support an etiological role for HTR2A in BD1. In view of the relatively small sample, replicate studies using large samples are needed.

Phenotypic heterogeneity and associations of two aldose reductase gene polymorphisms with nephropathy and retinopathy in type 2 diabetes

OBJECTIVE

We investigated the phenotypic features of diabetic microvascular complications and their association with a (CA)(n) microsatellite and a C/T polymorphism at the 5' region of the aldose reductase gene (ALR2) in a consecutive cohort of 738 Chinese type 2 diabetic patients.

RESEARCH DESIGN AND METHODS

Of the entire patient cohort, 392 were free of diabetes complications, or uncomplicated, 159 had diabetic nephropathy, 66 had diabetic retinopathy, and 121 had both diabetic nephropathy and retinopathy. Nephropathy was defined as urinary albumin excretion rate (AER) >or=20 micro g/min and albumin-to-creatinine ratio >or=3.5 mg/mmol in two urine collections. Retinopathy was defined by the presence of hemorrhages, exudates, laser marks, and fibrous proliferation or by a history of vitrectomy. (CA)(n) and C/T polymorphisms were examined by PCR followed by capillary electrophoresis and digestion with BfaI, respectively.

RESULTS

In the whole cohort, patients with diabetic retinopathy (n = 187) had higher blood pressure and lower BMI, while those with diabetic nephropathy (n = 280) had higher blood pressure, waist-to-hip ratio, and lipid profile than those without the respective complications. The z+6 carriers of the (CA)(n) polymorphism were less common in patients with diabetic retinopathy than those without diabetic retinopathy (n = 551) (4.3 vs. 9.3%, P = 0.04). The CT/TT carriers had a higher AER than the CC carriers (30.2 x/divided by 7.2 vs. 21.9 x/divided by 6.9 micro g/min, P = 0.03). Further subgroup analysis was performed after excluding uncomplicated patients with <5 years disease duration. The group with both diabetic nephropathy and retinopathy had higher frequencies of the z-2 allele (25.7 vs. 16.9%, P = 0.03) and T allele (26.4 vs. 18.5%, P = 0.04) and a lower frequency of the z+6 allele (1.7 vs. 5.5%, P = 0.054) than the uncomplicated group. Multiple logistic regression analysis confirmed that z-2 carrying (odds ratio 2.6, 95% CI 1.20-5.83, P = 0.02) and CT/TT genotypes (OR 2.5, 95% CI 1.19-5.19, P = 0.02) were independent predictors for both diabetic nephropathy and retinopathy.

CONCLUSIONS

Chinese type 2 diabetic patients exhibited phenotypic differences in terms of risk factors for both diabetic nephropathy and diabetic retinopathy. Both the z-2 allele of (CA)(n) polymorphism and T allele of ALR2 were independently associated with severe diabetic microvascular complications.

Insulin degrading enzyme (IDE) genetic variants and risk of Alzheimer's disease: evidence of effect modification by apolipoprotein E (APOE)

Insulin degrading enzyme (IDE) is a protease that degrades insulin and the beta-amyloid peptide implicated in Alzheimer's disease (AD). We reexamined data on five previously reported IDE polymorphisms stratifying the analysis by the presence or absence of an apolipoprotein E (APOE) epsilon4 allele. Three IDE variants were associated with AD within epsilon4-negative subjects (genotype exact test P-values < or =0.02). A haplotype containing the minor variant at each of these sites represented an estimated 4.2% of case haplotypes versus 12.3% of control haplotypes among epsilon4-negative subjects. Lack of this minor haplotype may be predictive of AD, potentially explaining some fraction of disease within subjects without the APOE epsilon4 risk allele. Confirmation of this finding with a larger sample of cases and controls is warranted.

Linkage disequilibrium analysis of the human adenosine deaminase (ada) gene provides evidence for a lack of correlation between hot spots of equal and unequal homologous recombination

The linkage disequilibrium (LD) pattern within the adenosine deaminase (ADA) gene was analyzed by studying 13 polymorphic loci in 137 families from two European and three African populations. Evidence for the presence of a 12-kb meiotic crossover hot spot, spanning part of the first and the second intron and flanked by regions of reduced recombination activity, was obtained. Moreover, segregation analysis of 113 informative meioses revealed two recombination events that are internal or overlap the 12-kb region, thus suggesting a recombination rate for the hot-spot region about 50-fold higher than the mean rate across the human genome. Within the hot spot, a 144-bp palindromic sequence was also identified and its possible involvement in the recombination process is discussed. The 12-kb region characterized by the low degree of LD does not include the 3.2-kb region that is deleted, as a result of recurrent unequal homologous recombination between two Alu elements, in patients affected by autosomal severe combined immunodeficiency. This observation provides the first evidence for an absence of correlation between hot spots of equal and unequal homologous recombination.

Interferon-alpha receptor-1 (IFNAR1) variants are associated with protection against cerebral malaria in the Gambia

The chromosome 21q22.11 cytokine receptor cluster contains four genes that encode subunits of the receptors for the cytokines interleukin-10 and interferon-alpha, -beta and -gamma that may have a role in malaria pathogenesis. A total of 15 polymorphic markers located within these genes were initially genotyped in 190 controls and 190 severe malaria cases from The Gambia. Two interferon-alpha receptor-1 (IFNAR1) gene SNPs (17470 and L168 V) showed evidence for an association with severe malaria phenotypes and were typed in a larger series of samples comprising 538 severe malaria cases, 338 mild malaria cases and 562 controls. Both the 17470-G/G and L168V-G/G genotypes were associated with protection against severe malaria, in general, and cerebral malaria, in particular (P=0.004 and 0.003, respectively). IFNAR1 diplotypes were then constructed for these two markers using the PHASE software package. The (17470-G L168V-G/17470-G L168V-G) diplotype was found to be associated with a reduced risk of cerebral malaria and the (17470-C L168V-C/17470-G L168V-G) diplotype with an increased risk of cerebral malaria (overall 3 x 2 chi(2)=12.8, d.f.=2, P=0.002 and 3 x 2 chi(2)=15.2, d.f.=2, P=0.0005, respectively). These data suggest a role for the type I interferon pathway in resistance to cerebral malaria.

Association of erosive hand osteoarthritis with a single nucleotide polymorphism on the gene encoding interleukin-1 beta

OBJECTIVE

Certain forms of primary osteoarthritis (OA), particularly those affecting hand joints, have a genetic component. Recent studies have shown suggestive evidence that hand and knee OA are linked with the interleukin-1 (IL-1) region on human chromosome 2q. This study was undertaken to assess the association of primary OA of the hand (hand OA) with IL-1 region markers.

METHODS

Sixty-eight US Caucasoid cases and 51 US Caucasoid controls aged 60 years or older were recruited from the Mid-Atlantic region of the United States. Hand OA was classified by American College of Rheumatology (ACR) Clinical Criteria, and cases were subjected to radiographic examination for subgrouping. Genotyping was done for seven previously described single nucleotide polymorphisms (SNPs) of genes for IL-1alpha (encoded by IL1A), IL-1beta (IL1B), and the IL-1 receptor antagonist (IL1RN), as well as an IL1RN variable number of tandem repeat (VNTR) marker. Six microsatellite markers on other chromosomes (null loci) were also typed.

RESULTS

The IL1B 5810 G>A SNP genotypes marker were not in Hardy-Weinberg equilibrium (p<0.05 in both non-erosive and erosive hand OA subgroups). Statistically significant association with the IL1B 5810 AA genotype was found in the erosive hand OA subgroup (relative risk 3.8, p=0.007). This IL1B 5810 AA genotype association was also significant between erosive and non-erosive hand OA subjects (relative risk 4.01, p=0.008). As expected, significant linkage disequilibrium was present between IL1B 5810 SNP and IL1A (-)889 SNP, other IL1B SNPs, and the nearest IL1RN SNP examined. The IL1B 5810A allele occurs most frequently on haplotypes with the SNP alleles IL1B 1423C, IL1B 1903T, IL1B 5887C, and IL1A (-)889C. Genotypes at null loci failed to show evidence suggesting population stratification that might account for spurious association.

CONCLUSION

Statistical evidence shows association between erosive hand OA and a genomic region containing the IL1B 5810 SNP in a US Caucasoid population. This supports a potential role for IL-1 in the pathogenesis of a severe phenotype of hand OA.

Dense linkage disequilibrium mapping in the 15q11-q13 maternal expression domain yields evidence for association in autism

Autism [MIM 209850] is a neurodevelopmental disorder exhibiting a complex genetic etiology with clinical and locus heterogeneity. Chromosome 15q11-q13 has been proposed to harbor a gene for autism susceptibility based on (1) maternal-specific chromosomal duplications seen in autism and (2) positive evidence for linkage disequilibrium (LD) at 15q markers in chromosomally normal autism families. To investigate and localize a potential susceptibility variant, we developed a dense single nucleotide polymorphism (SNP) map of the maternal expression domain in proximal 15q. We analyzed 29 SNPs spanning the two known imprinted, maternally expressed genes in the interval (UBE3A and ATP10C) and putative imprinting control regions. With a marker coverage of 1/10 kb in coding regions and 1/15 kb in large 5' introns, this map was employed to thoroughly dissect LD in autism families. Two SNPs within ATP10C demonstrated evidence for preferential allelic transmission to affected offspring. The signal detected at these SNPs was stronger in singleton families, and an adjacent SNP demonstrated transmission distortion in this subset. All SNPs showing allelic association lie within islands of sequence homology between human and mouse genomes that may be part of an ancestral haplotype containing a functional susceptibility allele. The region was further explored for recombination hot spots and haplotype blocks to evaluate haplotype transmission. Five haplotype blocks were defined within this region. One haplotype within ATP10C displayed suggestive evidence for preferential transmission. Interpretation of these data will require replication across data sets, evaluation of potential functional effects of associated alleles, and a thorough assessment of haplotype transmission within ATP10C and neighboring genes. Nevertheless, these findings are consistent with the presence of an autism susceptibility locus in 15q11-q13.

Complete sequencing shows a role for MSX1 in non-syndromic cleft lip and palate

MSX1 has been proposed as a gene in which mutations may contribute to non-syndromic forms of cleft lip and/or cleft palate. Support for this comes from human linkage and linkage disequilibrium studies, chromosomal deletions resulting in haploinsufficiency, a large family with a stop codon mutation that includes clefting as a phenotype, and the Msx1 phenotype in a knockout mouse. This report describes a population based scan for mutations encompassing the sense and antisense transcribed sequence of MSX1 (two exons, one intron). We compare the completed genomic sequence of MSX1 to the mouse Msx1 sequence to identify non-coding homology regions, and sequence highly conserved elements. The samples studied were drawn from a panethnic collection including people of European, Asian, and native South American ancestry. The gene was sequenced in 917 people and potentially aetiological mutations were identified in 16. These included missense mutations in conserved amino acids and point mutations in conserved regions not identified in any of 500 controls sequenced. Five different missense mutations in seven unrelated subjects with clefting are described. Evolutionary sequence comparisons of all known Msx1 orthologues placed the amino acid substitutions in context. Four rare mutations were found in non-coding regions that are highly conserved and disrupt probable regulatory regions. In addition, a panel of 18 population specific polymorphic variants were identified that will be useful in future haplotype analyses of MSX1. MSX1 mutations are found in 2% of cases of clefting and should be considered for genetic counselling implications, particularly in those families in which autosomal dominant inheritance patterns or dental anomalies appear to be cosegregating with the clefting phenotype.

Positional cloning of a quantitative trait locus on chromosome 13q14 that influences immunoglobulin E levels and asthma

Atopic or immunoglobulin E (IgE)-mediated diseases include the common disorders of asthma, atopic dermatitis and allergic rhinitis. Chromosome 13q14 shows consistent linkage to atopy and the total serum IgE concentration. We previously identified association between total serum IgE levels and a novel 13q14 microsatellite (USAT24G1; ref. 7) and have now localized the underlying quantitative-trait locus (QTL) in a comprehensive single-nucleotide polymorphism (SNP) map. We found replicated association to IgE levels that was attributed to several alleles in a single gene, PHF11. We also found association with these variants to severe clinical asthma. The gene product (PHF11) contains two PHD zinc fingers and probably regulates transcription. Distinctive splice variants were expressed in immune tissues and cells.

Novel association suggests multiple independent QTLs within chromosome 5q21-33 region control variation in total humans IgE levels

Asthma is a common, heterogeneous, complex disease accompanied by raised total and specific immunoglobulin-E (IgE) antibody levels. Despite numerous previous reports of linkage and association of asthma, atopy and serum IgE levels to genes within the 5q21-33 region, definitive, replicable results are still not available. We used the classical twin design to (i) estimate the relative contributions of genes and environment to variation in total IgE levels, (ii) assess genetic linkage, and (iii) examine allelic association of 11 microsatellite markers spanning the 5q21-33 region to total IgE. Variation in total IgE level was shown to be highly heritable (65%). Although evidence for linkage of the 11 microsatellites to IgE was not observed, the omnibus test of association, not confounded by population substructure, showed positive association of D5S393 and D5S673 to IgE. Genes in the vicinity of D5S673 include hepatitis A virus receptor (HAVCR-1) and IL-12B. Recently, the mouse orthologue of HAVCR-1, the T-cell membrane family of proteins, have been shown to be in strong association with expression of airway hyperactivity in a mouse model of human asthma and atopy. IL-12B subserves many proinflammatory functions and also induces B cells proliferation.

Association of specific language impairment (SLI) to the region of 7q31

FOXP2 (forkhead box P2) was the first gene characterized in which a mutation affects human speech and language abilities. A common developmental language disorder, specific language impairment (SLI), affects 6%-7% of children with normal nonverbal intelligence and has evidence of a genetic basis in familial and twin studies. FOXP2 is located on chromosome 7q31, and studies of other disorders with speech and language impairment, including autism, have found linkage to this region. In the present study, samples from children with SLI and their family members were used to study linkage and association of SLI to markers within and around FOXP2, and samples from 96 probands with SLI were directly sequenced for the mutation in exon 14 of FOXP2. No mutations were found in exon 14 of FOXP2, but strong association was found to a marker within the CFTR gene and another marker on 7q31, D7S3052, both adjacent to FOXP2, suggesting that genetic factors for regulation of common language impairment reside in the vicinity of FOXP2.

Endothelial nitric oxide synthase polymorphisms are associated with type 2 diabetes and the insulin resistance syndrome

Endothelial nitric oxide synthase (eNOS) variants were previously demonstrated in cardiovascular disease. To evaluate whether eNOS gene variants are associated with insulin resistance and type 2 diabetes, we evaluated polymorphisms in Exon7 (E298D), intron 18 (IVS18 + 27A-->C), and intron 23 (IVS23 + 10G-->T) in 159 type 2 diabetic patients without macrovascular complications and in 207 healthy control subjects. Samples for all hormonal and metabolic variables were obtained after an overnight fast. The D298 and IVS18 + 27C alleles, but not the IVS23 + 10G-->T variant, were significantly more frequent in type 2 diabetic patients than in control subjects. The two- and three-loci haplotype analysis showed that there is a statistically significant association between the eNOS variants and type 2 diabetes. No significant differences were observed in the clinical characteristics of type 2 diabetic patients according to genotypes (except for visceral obesity [waist-to-hip ratio], which was significantly more present in D298 homozygotes). Healthy control subjects homozygous for both D298 and IVS18 + 27C presented higher insulin, C-peptide, and nitric oxide levels, as well as higher HOMA (homeostasis model assessment) values than the double wild-type homozygotes, with values superimposable on those found in type 2 diabetic patients. In conclusion, we described a significant association between eNOS gene polymorphisms and type 2 diabetes, suggesting a new genetic susceptibility factor for hyperinsulinemia, insulin resistance, and type 2 diabetes.

Patterns of linkage disequilibrium in mitochondrial DNA of 16 ruminant populations

Mitochondrial DNA (mtDNA) is a widely employed molecular tool in phylogeography, in the inference of human evolutionary history, in dating the domestication of livestock and in forensic science. In humans and other vertebrates the popularity of mtDNA can be partially attributed to an assumption of strict maternal inheritance, such that there is no recombination between mitochondrial lineages. The recent demonstration that linkage disequilibrium (LD) declines as a function of distance between polymorphic sites in hominid mitochondrial genomes has been interpreted as evidence of recombination between mtDNA haplotypes, and hence nonclonal inheritance. However, critics of mtDNA recombination have suggested that this association is an artefact of an inappropriate measure of LD or of sequencing error, and subsequent studies of other populations have failed to replicate the initial finding. Here we report the analysis of 16 ruminant populations and present evidence that LD significantly declines with distance in five of them. A meta-analysis of the data indicates a nonsignificant trend of LD declining with distance. Most of the earlier criticisms of patterns between LD and distance in hominid mtDNA are not applicable to this data set. Our results suggest that either ruminant mtDNA is not strictly clonal or that compensatory selection has influenced patterns of variation at closely linked sites within the mitochondrial control region. The potential impact of these processes should be considered when using mtDNA as a tool in vertebrate population genetic, phylogenetic and forensic studies.

Estimation of linkage disequilibrium in a sample of the United Kingdom dairy cattle population using unphased genotypes

The association between genetic marker alleles was estimated for two regions of the bovine genome from a random sample of 50 young dairy bulls born in the United Kingdom between 1988 and 1995. Microsatellite marker genotypes were obtained for six markers on chromosome 2 and seven markers on chromosome 6, spanning 38 and 20 cM, respectively. Two different methods, which do not require family information, were used to estimate population haplotype frequencies. Haplotype frequencies were estimated for pairs of loci using the expectation-maximization algorithm and for all linked loci using a Bayesian approach via a Markov chain-Monte Carlo algorithm. Significant (P = 0.0007) linkage disequilibrium was detected between pairs of loci in syntenic groups (that is, loci in the same linkage group), extending to about 10 cM. No significant linkage disequilibrium was detected between markers in nonsyntenic regions. Given the observed level of linkage disequilibrium, mapping methods based on population-wide association might provide a better resolution than traditional quantitative trait loci mapping methods in the U.K. dairy cattle population and may reduce the required sample sizes of the experiments.