HLA and insulin gene associations with IDDM

Common variants in the TCF7L2 gene help to differentiate autoimmune from non-autoimmune diabetes in young (15-34 years) but not in middle-aged (40-59 years) diabetic patients

AIMS/HYPOTHESIS

Type 1 diabetes in children is characterised by autoimmune destruction of pancreatic beta cells and the presence of certain risk genotypes. In adults the same situation is often referred to as latent autoimmune diabetes in adults (LADA). We tested whether genetic markers associated with type 1 or type 2 diabetes could help to discriminate between autoimmune and non-autoimmune diabetes in young (15-34 years) and middle-aged (40-59 years) diabetic patients.

METHODS

In 1,642 young and 1,619 middle-aged patients we determined: (1) HLA-DQB1 genotypes; (2) PTPN22 and INS variable-number tandem repeat (VNTR) polymorphisms; (3) two single nucleotide polymorphisms (rs7903146 and rs10885406) in the TCF7L2 gene; (4) glutamic acid decarboxylase (GAD) and IA-2-protein tyrosine phosphatase-like protein (IA-2) antibodies; and (5) fasting plasma C-peptide.

RESULTS

Frequency of risk genotypes HLA-DQB1 (60% vs 25%, p = 9.4 x 10(-34); 45% vs 18%, p = 1.4 x 10(-16)), PTPN22 CT/TT (34% vs 26%, p = 0.0023; 31% vs 23%, p = 0.034), INS VNTR class I/I (69% vs 53%, p = 1.3 x 10(-8); 69% vs 51%, p = 8.5 x 10(-5)) and INS VNTR class IIIA/IIIA (75% vs 63%, p = 4.3 x 10(-6); 73% vs 60%, p = 0.008) was increased in young and middle-aged GAD antibodies (GADA)-positive compared with GADA-negative patients. The type 2 diabetes-associated genotypes of TCF7L2 CT/TT of rs7903146 were significantly more common in young GADA-negative than in GADA-positive patients (53% vs 43%; p = 0.0004). No such difference was seen in middle-aged patients, in whom the frequency of the CT/TT genotypes of TCF7L2 was similarly increased in GADA-negative and GADA-positive groups (55% vs 56%).

CONCLUSIONS/INTERPRETATION

Common variants in the TCF7L2 gene help to differentiate young but not middle-aged GADA-positive and GADA-negative diabetic patients, suggesting that young GADA-negative patients have type 2 diabetes and that middle-aged GADA-positive patients are different from their young GADA-positive counterparts and share genetic features with type 2 diabetes.

Association of Single Nucleotide Polymorphisms of the Insulin Gene with Chicken Early Growth and Fat Deposition

Growth rate, body composition, and fat deposition are important traits in chickens. Insulin plays important roles in hepatic cells, muscle cells, and adipose tissue cells. The purpose of the present study was to analyze association of the insulin (INS) gene with chicken growth and body composition traits. Using a F2 design resource population constructed with the crossing of Chinese native Xinghua chickens and White Recessive Rock chickens, the association of 4 single nucleotide polymorphisms (SNP; A+428G, C+1549T, T+3737C, and A+3971G) of INS gene with 13 growth and body composition traits was studied. The T+3737C genotypes were significantly associated with small intestine length (P = 0.0002), and the A+3971G genotypes were significantly associated with early growth (hatch weight and BW at 28 d of age) (P < 0.0001), breast angle (P = 0.0002), and small intestine length (P < 0.0001). None of the 4 SNP was significantly associated with abdominal fat pad weight (P > 0.05). The haplotypes based on the 4 SNP were also significantly associated with early growth (hatch weight and BW at 28 d of age; P < 0.0001) and breast angle (P < 0.0001) but not with small intestine length (P = 0.0505). These results suggested that variation of the insulin gene was significantly associated with chicken early growth but not with fat deposition. In addition, the data from the present study supported the inference that both the one-SNP-at-a-time and the haplotype-based approaches have their own advantages and disadvantages when association analysis of one SNP and haplotypes with chicken complex traits was conducted.

What is the significance of a significant TDT?

This note summarizes the development of the transmission/disequilibrium test (TDT). The initial purpose of the TDT procedure was to test for linkage between a genetic marker and a disease susceptibility locus when an association had been found between the two. An association between disease and marker had sometimes been taken to imply linkage. An association could, however, be due to population stratification even in the absence of linkage. In contrast, the outcome of the TDT is not affected by such stratification. Furthermore, when linkage is not in doubt, the TDT can, in some cases, also provide a test of association between marker and disease. We discuss these various matters in this paper.

Genetic interaction among three genomic regions creates distinct contributions to early- and late-onset type 1 diabetes mellitus

There are two peaks in the distribution of the age of onset of type 1 diabetes mellitus (T1DM)--the first in early childhood and the second at the time of puberty. Although T1DM results from the interaction of genetic and non-genetic factors, it has not been established which factors contribute to the bimodal distribution. The genetic component of T1DM is in large part related to genes from the human leukocyte antigen (HLA) complex (IDDM1); however, loci from the variable nucleotide tandem repeat (VNTR) region of the insulin (INS) gene (IDDM2) and more recently, the cytotoxic T-lymphocyte-associated protein-4 region (CTLA4, IDDM12) have also been implicated. Therefore, we examined the potential interaction between these loci through the influence of the age of onset of T1DM in diabetic and control Caucasian individuals. We discovered that younger individuals with HLA-DRB1*0301/DRB1*04 and INS I/I genotypes exhibited increased susceptibility to T1DM, whereas the interaction of INS I/I and CTLA4 G/G genotypes was more common in older children with T1DM. Combining the age of onset of T1DM with specific genotypes may operate to produce a single disease through different underlying causes.

The stages of type 1A diabetes: 2005

Type 1A diabetes is a chronic autoimmune disease usually preceded by a long prodrome during which autoantibodies to islet autoantigens are present. These antibodies are directed to a variety of antigens, but the best characterized are glutamic acid decarboxylase-65, insulinoma-associated antigen-2, and insulin. We hypothesize that the natural history of type 1A diabetes can be represented by several stages, starting from genetic susceptibility and ending in complete beta-cell destruction and overt diabetes. Type 1A diabetes probably results from a balance between genetic susceptibility and environmental influences. In both humans and animal models, the major determinants of the disease are genes within the major histocompatibility complex. The next best-characterized susceptibility locus is the insulin gene, the variable nucleotide tandem repeat locus. This gene affects the expression of insulin in the thymus and thus may play a role in the modulation of tolerance to this molecule. In a subset of genetically susceptible individuals, the activation of autoimmunity may be triggered by environmental factors such as viruses and/or diet. However, no conclusive association has been established between type 1A diabetes and specific environmental triggers. In this review, we provide evidence that insulin has a fundamental role in anti-islet autoimmunity.

CTLA-4 gene polymorphisms and susceptibility to type 1 diabetes mellitus: a HuGE Review and meta-analysis

The authors performed a meta-analysis of 33 studies examining the association of type 1 diabetes mellitus with polymorphisms in the cytotoxic T-lymphocyte-associated antigen-4 (CTLA-4) gene, including the A49G (29 comparisons), C(-318)T (three comparisons), and (AT)n microsatellite (six comparisons) polymorphisms. The studies included 5,637 cases of type 1 diabetes and 6,759 controls (4,775 and 5,829, respectively, for analysis of the A49G polymorphism). The random-effects odds ratio for the *G (Ala) allele versus the *A (Thr) allele was 1.45 (95% confidence interval (CI): 1.28, 1.65), with significant between-study heterogeneity (p < 0.001). The effect size tended to be higher in type 1 diabetes cases with age of onset <20 years (odds ratio (OR) = 1.61), and there was a significant association between the presence of glutamic acid decarboxylase-65 autoantibodies and the *G allele among type 1 diabetes cases (OR = 1.49). Larger studies showed more conservative results (p = 0.011). After exclusion of studies with fewer than 150 subjects and studies with significant deviation from Hardy-Weinberg equilibrium in the controls, the summary odds ratio was 1.40 (95% CI: 1.28, 1.54). Available data showed no strong association for the 106-base-pair allele of the microsatellite polymorphism (OR = 0.99, 95% CI: 0.64, 1.55) or the *T allele of the C(-318)T polymorphism (OR = 0.92, 95% CI: 0.45, 1.89). This meta-analysis demonstrates that the CTLA-4*G genotype is associated with type 1 diabetes.

Remapping the insulin gene/IDDM2 locus in type 1 diabetes

Type 1 diabetes susceptibility at the IDDM2 locus was previously mapped to a variable number tandem repeat (VNTR) 5' of the insulin gene (INS). However, the observation of associated markers outside a 4.1-kb interval, previously considered to define the limits of IDDM2 association, raised the possibility that the VNTR association might result from linkage disequilibrium (LD) with an unknown polymorphism. We therefore identified a total of 177 polymorphisms and obtained genotypes for 75 of these in up to 434 pedigrees. We found that, whereas disease susceptibility did map to within the 4.1-kb region, there were two equally likely candidates for the causal variant, -23HphI and +1140A/C, in addition to the VNTR. Further analyses in 2,960 pedigrees did not support the difference in association between VNTR lineages that had previously enabled the exclusion of these two polymorphisms. Therefore, we were unable to rule out -23HphI and +1140A/C having an etiological effect. Our mapping results using robust regression methods show how precisely a variant for a common disease can be mapped, even within a region of strong LD, and specifically that IDDM2 maps to one or more of three common variants in a approximately 2-kb region of chromosome 11p15.

Surfactant protein A and B genetic variants in respiratory distress syndrome in singletons and twins

Interactive genetic and environmental factors may influence the differentiation of surfactant and the risk of respiratory distress syndrome (RDS). DNA samples from 441 premature singleton infants and 480 twin or multiple infants were genotyped for surfactant-specific protein (SP)-A1, SP-A2, and SP-B exon 4 polymorphisms and intron 4 size variants in a homogeneous white population. Distributions of the SP-A and SP-B gene variants between RDS and no-RDS infants were determined alone and in combination. SP-A1 allele 6A2 (p = 0.009) and the homozygous genotype 6A2/6A2 (p = 0.003) were overrepresented in RDS of singletons when the SP-B exon 4 genotype was Thr/Thr, and underrepresented in RDS of multiples when the SP-B genotype was Ile/Thr (p = 0.012 for 6A2 and p = 0.03 for 6A2/6A2) or Thr/Thr (p = 0.12 for 6A2 and p = 0.018 for 6A2/6A2, respectively). The SP-A 6A2 allele in the SP-B Thr131 background predisposed the smallest singleton infants to RDS, whereas near-term multiples were protected from RDS. There was a continuous association between fetal mass and risk of RDS, defined by the SP-A and SP-B variants. Labeled lung explants with the Thr/Thr genotype showed proSP-B amino-terminal glycosylation, which was absent in Ile/Ile samples. Genetic and environmental variation may influence intracellular processing of surfactant complex and the susceptibility to RDS.

Nonparametric disequilibrium mapping of functional sites using haplotypes of multiple tightly linked single-nucleotide polymorphism markers

As the speed and efficiency of genotyping single-nucleotide polymorphisms (SNPs) increase, using the SNP map, it becomes possible to evaluate the extent to which a common haplotype contributes to the risk of disease. In this study we propose a new procedure for mapping functional sites or regions of a candidate gene of interest using multiple linked SNPs. Based on a case-parent trio family design, we use expectation-maximization (EM) algorithm-derived haplotype frequency estimates of multiple tightly linked SNPs from both unambiguous and ambiguous families to construct a contingency statistic S for linkage disequilibrium (LD) analysis. In the procedure, a moving-window scan for functional SNP sites or regions can cover an unlimited number of loci except for the limitation of computer storage. Within a window, all possible widths of haplotypes are utilized to find the maximum statistic S* for each site (or locus). Furthermore, this method can be applied to regional or genome-wide scanning for determining linkage disequilibrium using SNPs. The sensitivity of the proposed procedure was examined on the simulated data set from the Genetic Analysis Workshop (GAW) 12. Compared with the conventional and generalized TDT methods, our procedure is more flexible and powerful.

A weak association of HLA-B*2702 with Behçet's disease

This study aimed to analyse the association of HLA-B alleles other than -B51 with Behçet's disease (BD). We also investigated the frequency of HLA-B alleles sharing the same natural killer cell immunoglobulin-like receptor (KIR) binding sequence with HLA-B51. Broad-genotyping of HLA-B locus by PCR-SSOP in 174 Turkish BD patients and 191 healthy controls confirmed the strong association of B*51 with BD (60.9% in BD patients, 24.6% in healthy controls, OR = 4.78). No other HLA-B allele was identified showing an association with BD after adjusting for multiple testing or by using relative predispositional effects (RPE) analysis after the deletion of B*51. HLA-B alleles reacting with the sequence specific oligonucleotide probe 23, which corresponds to the KIR binding site of B*51, were found to be positive in 127 BD patients (73%) and 90 controls (47%) (OR = 3.03, 95% CI 2-4.7). The repeated RPE analysis after separating HLA-B alleles carrying B51-KIR binding sequence as distinct alleles within a broad-type allele group revealed B*2702 allele as the only allele showing an association with BD after the deletion of B*51. Selective increase of B*2702, the only B*27 allele carrying the same KIR binding sequence with B*51, warrants investigation of the possibility of interaction of HLA molecules with KIRs on NK or other T cells in the pathogenesis of BD.

Treasure hunting in a new era: genotyping of single nucleotide polymorphisms and the search for complex disease genes by association scans

Genotyping single nucleotide polymorphisms by technologies like those described in this issue are propelling us towards a new era in the study of human genetics. This new era promises improved prospects for mapping genes responsible for complex human diseases and other complex traits. In this article, I sketch empirical results and statistical theory that together form the scientific rationale for anticipating enhanced success in gene mapping in the new era that is now approaching.

Molecular genetics of schizophrenia: past, present and future

Schizophrenia is a severe neuropsychiatric disorder with a polygenic mode of inheritance which is also governed by non-genetic factors. Candidate genes identified on the basis of biochemical and pharmacological evidence are being tested for linkage and association studies. Neurotransmitters, especially dopamine and serotonin have been widely implicated in its etiology. Genome scan of all human chromosomes with closely spaced polymorphic markers is being used for linkage studies. The completion and availability of the first draft of Human Genome Sequence has provided a treasure-trove that can be utilized to gain insight into the so far inaccessible regions of the human genome. Significant technological advances for identification of single nucleo-tide polymorphisms (SNPs) and use of microarrays have further strengthened research methodologies for genetic analysis of complex traits. In this review, we summarize the evolution of schizophrenia genetics from the past to the present, current trends and future direction of research.

Use of unlinked genetic markers to detect population stratification in association studies

We examine the issue of population stratification in association-mapping studies. In case-control studies of association, population subdivision or recent admixture of populations can lead to spurious associations between a phenotype and unlinked candidate loci. Using a model of sampling from a structured population, we show that if population stratification exists, it can be detected by use of unlinked marker loci. We show that the case-control-study design, using unrelated control individuals, is a valid approach for association mapping, provided that marker loci unlinked to the candidate locus are included in the study, to test for stratification. We suggest guidelines as to the number of unlinked marker loci to use.

Associations of MHC class II alleles with insulin-dependent diabetes mellitus (IDDM) in patients from North India

Thirty-four insulin-dependent diabetes mellitus (IDDM) patients from North India were studied with respect to their HLA class II alleles including those of the DRB1, DQA1, DQB1 and DPB1 loci, using the polymerase chain reaction (PCR) and hybridization with sequence-specific oligonucleotide probes (SSOP). They were compared with the class II alleles of 94 normal adult controls from the same ethnic background. The results show a statistically significant increase of DRB1*03011 (p < 0.00001), DQB1*0201 (p < 0.007), DQA1*0501 (0.0027) and DPB1*2601 (p < 0.0042) in patients compared to controls. DR*04 was not significantly increased. However, homozygosity for DRB1*03011 was increased more than expected. DRB1*1501 and *1502 did not show a significant decrease in the patients. However, DRB1*0701 was decreased significantly, but this difference did not remain significant when the p value was corrected for the number of alleles tested. Similarly, DPB1*2601 was increased significantly in the patients but did not remain significant after p was corrected for the number of alleles tested. However, DPB1*2601 was increased, and remained significant after correction, in patients not having HLA-DR3. We also studied the possible role of aspartic acid at codon 57 of the DQ beta chain in protection against development of diabetes, and arginine at codon 52 of the DQ alpha chain in susceptibility. We observed an increase in non-Asp57 alleles in DQ beta and Arg52 in DQ alpha in the patients, however, this effect seems to be due to the fact that the most prevalent haplotype in diabetic patients: DRB1*03011-DQA1*0501-DQB1*0201, has DQB1 and DQA1 alleles which carry non-Asp57 and Arg52, respectively.

A second-generation screen of the human genome for susceptibility to insulin-dependent diabetes mellitus

During the past decade, the genetics of type 1 (insulin-dependent) diabetes mellitus (IDDM) has been studied extensively and the disorder has become a paradigm for genetically complex diseases. Previous genome screens and studies focused on candidate genes have provided evidence for genetic linkage between polymorphic DNA markers and 15 putative IDDM susceptibility loci, designated IDDM1-IDDM15. We have carried out a second-generation screen of the genome for linkage and analysed the data by multipoint linkage methods. An initial panel of 212 affected sibpairs (ASPs) was genotyped for 438 markers spanning all autosomes, and an additional 467 ASPs were used for follow-up genotyping. Other than the well-established linkage with the HLA region at chromosome 6p21.3, there was only one region, located on chromosome 1q and not previously reported, where the log likelihood ratio (lod) was greater than 3. Lods between 1.0 and 1.8 were found in six other regions, three of which have been reported in other studies. Another reported region, on chromosome 6q and loosely linked to HLA, also had an elevated lod. Little or no support was found for most reported IDDM loci (lods were less than 1), despite larger sample sizes in the present study.

Susceptibility to relapsing-progressive multiple sclerosis is associated with inheritance of genes linked to the variable region of the TcR beta locus: use of affected family-based controls

We tested the hypothesis that susceptibility to relapsing-progressive (RP) (but not to relapsing-remitting [RR]) multiple sclerosis (MS) is associated with a gene linked to the TcR beta-chain variable region delimited by the Vbeta8-BamHI and Vbeta11-BamHI RFLP alleles in DRw15+ MS patients, using a contingency-table test of patient data and affected family-based controls. Control alleles and haplotypes were composed of parental marker alleles and haplotypes not transmitted to the affected child, in 90 simplex and 31 multiplex families from British Columbia. A total of 6,164 alleles at 11 loci were segregated through families of probands with RP MS or RR MS. The Vbeta8-Vbeta11 subhaplotype frequencies in the DRw15+ RP MS (but not RR MS) patients differed from control frequencies, because of an increase of the 2-1 subhaplotype (P=.02). Vbeta8-BamHI and Vbeta11-BamHI allele frequencies (P=.05 and .009, respectively) in the DRw15+ RP MS (but not RR MS) patients differed from control frequencies. The Vbeta1-Vbeta8 subhaplotype frequencies in the DRw15- RP MS (but not RR MS) patients differed from control frequencies (P=.03), with a significantly increased frequency of the 1-1 subhaplotype (P=.01; RR=7.1) in RP MS versus RR MS patients. Susceptibility to RP MS is associated both with a recessive inheritance of a gene linked to the 3' (Vbeta11) end of the 2-1 subhaplotype defined by the Vbeta8-BamHI and Vbeta11-BamHI alleles in DRw15+ patients and with a gene, located on the 1-1 subhaplotype, defined by the Vbeta1-TaqI and Vbeta8-MspI alleles of the TcR beta-chain complex in DRw15- patients.

Genetics of sarcoidosis

Hereditary susceptibility to sarcoidosis is suggested by ethnic preponderance, familial clustering, and multigenerational involvement. The genetics of sarcoidosis cannot be adequately addressed in small samples of patients; a large-scale study with stratification for patient phenotypic differences is necessary. A study that uses both genetic marker and environmental data would be able to control for and examine different causative mechanisms. Until such a well-designed, comprehensive study is carried out, we are left with interesting patterns of disease in families and uncertain allelic associations.

Association mapping of disease loci, by use of a pooled DNA genomic screen

Genomic screening to map disease loci by association requires automation, pooling of DNA samples, and 3,000-6,000 highly polymorphic, evenly spaced microsatellite markers. Case-control samples can be used in an initial screen, followed by family-based data to confirm marker associations. Association mapping is relevant to genetic studies of complex diseases in which linkage analysis may be less effective and to cases in which multigenerational data are difficult to obtain, including rare or late-onset conditions and infectious diseases. The method can also be used effectively to follow up and confirm regions identified in linkage studies or to investigate candidate disease loci. Study designs can incorporate disease heterogeneity and interaction effects by appropriate subdivision of samples before screening. Here we report use of pooled DNA amplifications-the accurate determination of marker-disease associations for both case-control and nuclear family-based data-including application of correction methods for stutter artifact and preferential amplification. These issues, combined with a discussion of both statistical power and experimental design to define the necessary requirements for detecting of disease loci while virtually eliminating false positives, suggest the feasibility and efficiency of association mapping using pooled DNA screening.

Human type 1 diabetes and the insulin gene: principles of mapping polygenes

We review the strategy used to identify a susceptibility locus (IDDM2) for type 1 (insulin dependent) diabetes mellitus. As type 1 diabetes is becoming the paradigm for dissecting multifactorial disease genetics, the approach described provides important general guidelines for positional cloning of human disease polygenes. Main topics include: (a) historical conspectus of the mapping and identification of IDDM2--a critical survey of the work leading up to the conclusion that IDDM2 most likely corresponds to allelic variation at the insulin gene minisatellite (VNTR) locus; (b) the nature of allelic (length and sequence) variation at the VNTR locus; (c) gene interactions and disease pathogenesis; (d) mechanism of action of the INS VNTR in type 1 diabetes--insulin gene expression, parent-of-origin effects (genomic imprinting); and (e) summary and future prospects--alleles of the insulin VNTR that are protective for type 1 diabetes appear to encode susceptibility to type 2 diabetes.

No linkage or association of telomeric and centromeric T-cell receptor beta-chain markers with susceptibility to type 1 insulin-dependent diabetes in HLA-DR4 multiplex families

The T-cell receptor beta locus (TCRB) on chromosome 7q35 was studied as a candidate region for genetic susceptibility to type 1 insulin-dependent diabetes (IDDM). A highly polymorphic microsatellite marker mapping to the TCRBV6.7 gene and a TCRB C-region RFLP were used to genotype the members of a total of 21 multiplex IDDM families from two different geographical areas. There was no evidence to support linkage to either of these markers with IDDM, and conventional two-point analysis excluded linkage to the telomeric end of the TCRB complex, in the region of the highly informative TCRBV6.7 marker. There was significant linkage of IDDM to the class II HLA-D locus with significant lod scores > 3.0 obtained for the HLA-DRB1 and HLA-DQB1 genes. Affected sib-pair (ASP) and transmission disequilibrium (TDT) association tests confirmed these findings.

Panning for gold: genome-wide scanning for linkage in type 1 diabetes

Genome-wide scans for linkage of chromosome regions to type 1 diabetes in affected sib pair families have revealed that the major susceptibility locus resides within the major histocompatibility complex (MHC) on chromosome 6p21 (lambda S = 2.4). It is recognized that the MHC contains multiple susceptibility loci (referred to collectively as IDDM1), including the class II antigen receptor genes, which control the major pathological feature of the disease: T-lymphocyte-mediated autoimmune destruction of the insulin-producing pancreatic beta cells. However, the MHC genes, and a second locus, the insulin gene minisatellite on chromosome 11p15 (IDDM2; lambda S = 1.25), cannot account for all of the observed clustering of disease in families (lambda S = 15), and the scans suggested the presence of other susceptibility loci scattered throughout the genome. There are four additional loci for which there is currently sufficient evidence from linkage and association studies to justify fine mapping experiments: IDDM4 (FGF3/11q13), IDDM5 (ESR/6q22), IDDM8 (D6S281/6q27) and IDDM12 (CTLA-4/2q33). IDDM4, 5 and 8 were detected by genome scanning, and IDDM12 by a candidate gene strategy. Seven other named loci are not discounted but remain to be replicated widely. Multiple susceptibility loci were expected as genome-wide scans of the mouse model of type 1 diabetes had shown that although the MHC is the major mouse locus, at least 13 genes unlinked to the MHC are involved in the development of disease. Genome-wide scans using 1000 affected sibpair families will be required to be confident that all genes with effects on familial clustering equivalent to the insulin gene locus (lambda S = 1.25) have been detected. The identification of aetiological determinants requires exclusion of hitchhiking polymorphisms in regions of linkage disequilibrium, as demonstrated for the MHC and the insulin gene loci, and functional studies implicating the disease-associated variant in pathogenesis. Ultimately, targeting of specific candidate mutations in mice by homologous recombination and replacement will be necessary to prove the primary role of any candidate mutation.

The geneticist's approach to complex disease

Many studies are in progress worldwide to elucidate the genetics of complex diseases. Nevertheless, few articles are available that provide the scientific rationale and give guidelines for such ambitious endeavours. We describe the methodology and background necessary to study the genetics of complex disease and discuss how to analyze the data. We also provide a table of some ongoing studies. In particular, we wish to emphasize the analysis of intermediate, heritable, quantitative traits as a means of dissecting the genetic basis of a complex trait.

Genetic analysis of type 1 diabetes using whole genome approaches

Whole genome linkage analysis of type 1 diabetes using affected sib pair families and semi-automated genotyping and data capture procedures has shown how type 1 diabetes is inherited. A major proportion of clustering of the disease in families can be accounted for by sharing of alleles at susceptibility loci in the major histocompatibility complex on chromosome 6 (IDDM1) and at a minimum of 11 other loci on nine chromosomes. Primary etiological components of IDDM1, the HLA-DQB1 and -DRB1 class II immune response genes, and of IDDM2, the minisatellite repeat sequence in the 5' regulatory region of the insulin gene on chromosome 11p15, have been identified. Identification of the other loci will involve linkage disequilibrium mapping and sequencing of candidate genes in regions of linkage.

Susceptibility to human type 1 diabetes at IDDM2 is determined by tandem repeat variation at the insulin gene minisatellite locus

The IDDM2 locus encoding susceptibility to type 1 diabetes was mapped previously to a 4.1-kb region spanning the insulin gene and a minisatellite or variable number of tandem repeats (VNTR) locus on human chromosome 11p15.5. By 'cross-match' haplotype analysis and linkage disequilibrium mapping, we have mapped the mutation IDDM2 to within the VNTR itself. Other polymorphisms were systematically excluded as primary disease determinants. Transmission of IDDM2 may be influenced by parent-of-origin phenomena. Although we show that the insulin gene is expressed biallelically in the adult pancreas, we present preliminary evidence that the level of transcription in vivo is correlated with allelic variation within the VNTR. Allelic variation at VNTRs may play an important general role in human disease.

Linkage disequilibrium mapping of a type 1 diabetes susceptibility gene (IDDM7) to chromosome 2q31-q33

The role of human chromosome 2 in type 1 diabetes was evaluated by analysing linkage and linkage disequilibrium at 21 microsatellite marker loci, using 348 affected sibpair families and 107 simplex families. The microsatellite D2S152 was linked to, and associated with, disease in families from three different populations. Our evidence localizes a new diabetes susceptibility gene, IDDM7, to within two centiMorgans of D2S152. This places it in a region of chromosome 2q that shows conserved synteny with the region of mouse chromosome 1 containing the murine type 1 diabetes gene, Idd5. These results demonstrate the utility of polymorphic microsatellites for linkage disequilibrium mapping of genes for complex diseases.

HLA disease associations: models for the study of complex human genetic disorders

The genes of the human leukocyte antigen (HLA) region, the major histocompatibility complex (MHC) of humans, control a variety of functions involved in immune response and influence susceptibility to over 40 diseases. Theoretical studies in the development of models to determine the modes of inheritance of the HLA-associated diseases have led to a better understanding of the inheritance patterns in insulin-dependent diabetes mellitus (IDDM), rheumatoid arthritis, multiple sclerosis, ankylosing spondylitis, hemochromatosis, celiac disease, and others. It is now clear that many of the HLA-associated diseases involve heterogeneity in their HLA components, as well as non-HLA genetic factors. This review is presented using HLA-associated diseases, and in particular IDDM, as the example of interest, but the observations and techniques presented have direct relevance to the study of all human diseases with a complex genetic component. Three methods for localizing disease-predisposing genes are presented: (1) association studies, including population, family, and relative predispositional effects, (2) affected sib pair and other affected-relative methods, and (3) lod score analysis. A variety of complementary methods for studying the mode(s) of inheritance of the alleles at the disease-predisposing locus and for identifying the alleles and amino acids directly involved in the disease process also are presented.

Relative-risk regression models using cases and their parents

Relative-risk regression models are presented for studies of the association of genetic markers with disease status when the study design uses affected cases and their parents, with or without unaffected sibs. These models generalize the "haplotype relative-risk" method and allow for censored unaffected sibs; in this sense, these models resemble proportional hazards models that are commonly used in survival analysis. A critical distinction between these models and the usual Cox proportional hazards model is that the frequencies of the genotypes of the cases are compared to controls based on Mendelian expectations, and not simply to the genotypes of the sib controls who are at risk of disease. These models allow modeling of the contribution of specific alleles to the relative risk of disease, as well as interactions of allelic effects with environmental risk factors. To demonstrate the application of these models, we have fit them to the binary affection status of the Problem 2 data set. Four candidate-gene loci were found to have a significant association with affection status, after allowance for relative risks that decrease with age; two of these associations correctly identified two of the major gene loci, and the other two were false-positive associations.

A genome-wide search for human type 1 diabetes susceptibility genes

We have searched the human genome for genes that predispose to type 1 (insulin-dependent) diabetes mellitus using semi-automated fluorescence-based technology and linkage analysis. In addition to IDDM1 (in the major histocompatibility complex on chromosome 6p21) and IDDM2 (in the insulin gene region on chromosome 11p15), eighteen different chromosome regions showed some positive evidence of linkage to disease. Linkages to chromosomes 11q (IDDM4) and 6q (IDDM5) were confirmed by replication, and chromosome 18 may encode a fifth disease locus. There are probably no genes with large effects aside from IDDM1. Therefore polygenic inheritance is indicated, with a major locus at the major histocompatibility complex.

HLA class II and susceptibility and resistance to insulin-dependent diabetes mellitus in a population from the northwest of Spain

The role of HLA class II alleles in genetic predisposition to insulin-dependent diabetes mellitus (IDDM) was examined using Polymerase Chain Reaction/oligonucleotide probe typing (PCR/SSOs) of eight HLA class II loci in 58 IDDM patients and 50 healthy controls from the Northwest of Spain (Asturias). We compared the distribution of HLA class II alleles, haplotypes and genotypes between IDDM patients and controls, and tested three recently proposed HLA-IDDM susceptibility theories. By using the aetiologic fraction (delta) as an almost absolute measure of the strongest linkage of disequilibrium of a HLA marker to the putative Type I susceptibility locus, it has been found that the strength of association of the HLA markers may be quantified as follows: DQA1*03-DQB1*0302 or DQA1*0501-DQB1*0201 > DR3 or DR4; presence of more than one dimer DQ alpha beta of the six proposed by Rønningen > non-Asp57 DQ beta and Arg52 DQ alpha > Arg52 DQ alpha > non-Asp57 DQ beta/non-Asp57 DQ beta > DRB1*0301; DQA1*0501-DQB1*0201 > DQA1*03-DQB1*0302; DQB1*0302. The presence of at least one Asp57 DQ beta allele was the best protection HLA marker to IDDM in our population. Therefore, the above data confirm that IDDM susceptibility to HLA locus is linked more to DQ than DR.

Involvement of human muscle acetylcholine receptor alpha-subunit gene (CHRNA) in susceptibility to myasthenia gravis

The muscle acetylcholine receptor is the major target of the autoimmune response in generalized myasthenia gravis. To investigate the role of the gene encoding the alpha subunit of the receptor (CHRNA), two stable polymorphic d[(GT).(CA)]dinucleotide repeats, designated HB and BB, were characterized within the first intron of CHRNA. The HB*14 allele conferred a relative risk for myasthenia gravis of 2.5 in 81 unrelated patients compared with 100 control subjects. Very significantly, family analysis based on haplotype segregation data indicated that parental haplotypes associated with HB*14 always segregated to the child with myasthenia gravis (P < 0.0002 for the comparison with the transmission of haplotypes not bearing HB*14), whereas their transmission to unaffected siblings was equilibrated. Myasthenia gravis patients also showed a high frequency of microsatellite variants unseen in controls. These findings implicate the CHRNA in susceptibility to myasthenia gravis.

Gene for familial psoriasis susceptibility mapped to the distal end of human chromosome 17q

A gene involved in psoriasis susceptibility was localized to the distal region of human chromosome 17q as a result of a genome-wide linkage analysis with polymorphic microsatellites and eight multiply affected psoriasis kindreds. In the family which showed the strongest evidence for linkage, the recombination fraction between a psoriasis susceptibility locus and D17S784 was 0.04 with a maximum two-point lod score of 5.33. There was also evidence for genetic heterogeneity and although none of the linked families showed any association with HLA-Cw6, two unlinked families showed weak levels of association. This study demonstrates that in some families, psoriasis susceptibility is due to variation at a single major genetic locus other than the human lymphocyte antigen locus.

Additive susceptibility to insulin-dependent diabetes conferred by HLA-DQB1 and insulin genes

Several genomic polymorphisms at the insulin (INS) gene and its flanking regions were analyzed in 197 unrelated Caucasian patients affected by insulin-dependent diabetes (IDDM) and 159 ethnically matched, normal controls ascertained from the South-Eastern United States. We found that the frequency of homozygotes for the common variant at the insulin gene was significantly increased in the diabetic population (RR = 2.0, p < 0.005). However, the polymorphisms in the 5' and 3' regions flanking the INS were not significantly associated with IDDM. These results suggest that the IDDM susceptibility locus on chromosome 11p is located within the region extending from the 5' VNTR to the 3' end of the INS gene. We determined the HLA-DQB1 genotypes by denaturing gradient gel electrophoresis (DGGE) and/or sequence-specific primers (SSP) techniques to assess the possible interactions between INS and HLA. DQB1*0302 had the strongest predisposing effect on IDDM susceptibility (RR = 9.3) and DQB1*0602 the strongest protective effect (RR = 0.02). However, a significant predisposing effect of DQB1*0201 could be demonstrated only after removal of the effects of DQB1*0302 and DQB1*0602. Analyses of the genotypes revealed that all genotypes containing 0602 were protective and that the heterozygous genotype 0201/0302 and homozygous genotype 0302/0302 confer the highest risk (RR = 20.9 and 12.9 respectively). However, heterozygous genotypes 0302/X (X excludes 0201, 0302 and 0602) have a significantly lower predisposing risk. Similarly, there is heterogeneity in risk between predisposing 0201/0201 homozygous individuals and protective 0201/X individuals. When subjects were stratified by HLA genotypes, the relative risks conferred by INS did not vary, thus suggesting that the susceptibility effects conferred by HLA and INS are additive rather than interactive.

Susceptibility to insulin dependent diabetes mellitus maps to a 4.1 kb segment of DNA spanning the insulin gene and associated VNTR

Recent studies have demonstrated that a locus at 11p15.5 confers susceptibility to insulin dependent diabetes mellitus (IDDM). This locus has been shown to lie within a 19 kb region. We present a detailed sequence comparison of the predominant haplotypes found in this region in a population of French Caucasian IDDM patients and controls. Identification of polymorphisms both associated and unassociated with IDDM has allowed us to define further the region of association to 4.1 kb. Ten polymorphisms within this region are in strong linkage disequilibrium with each other and extend across the insulin gene locus and the variable number tandem repeat (VNTR) situated immediately 5' to the insulin gene. These represent a set of candidate disease polymorphisms one or more of which may account for the susceptibility to IDDM.

Susceptibility for multiple sclerosis is determined, in part, by inheritance of a 175-kb region of the TcR V beta chain locus and HLA class II genes

Genetic makeup thought to affect susceptibility to multiple sclerosis (MS) and current evidence suggests that multiple genes may be involved. We have mapped a potential susceptibility gene or genes in the germ-line T cell receptor (TcR) V beta region of multiple sclerosis (MS) patients. Six restriction fragment length polymorphisms (RFLPs) spanning approximately 600 kb of the TcR V beta region were used to define TcR haplotypes in 197 Caucasian controls and 83 Caucasian MS patients in the chronic progressive stage of the disease. The distribution of TcR subhaplotype frequencies was significantly different only in the approx. 175-kb region between RFLPs defined by V beta 8.1 and V beta 11. Stratification of the MS patients into HLA-DR2+ (n = 51) and HLA-DR2- (n = 32) populations demonstrated that the subhaplotype frequencies differed from the control population significantly only in the HLA-DR2+ (corrected P = 0.00007) and not in the HLA-Dr2- (corrected P = 0.46) population. Subhaplotypes which are rare in the normal population are overrepresented in the HLA-DR2+ MS patient population and confer a relative risk of 4.06. These results indicate the existence of an MS susceptibility gene within the TcR V beta region, and provide new evidence for gene complementation between a HLA class II gene and TcR V beta gene(s) in conferring susceptibility to MS.

T cell receptor haplotypes in families of patients with insulin-dependent diabetes mellitus

The frequencies of Bgl 11 and BamH1 restriction fragment length polymorphisms (RFLP) of C beta, V beta 8, V beta 11 and V beta 7.2 have been defined in a healthy Australian population. Linkage disequilibrium between alleles of the T cell receptor (TCR) V beta 8 and V beta 11 gene segments has been confirmed. We have also confirmed the lack of linkage disequilibrium between either of these loci and alleles at C beta or V beta 7.2. Using RFLPs at V beta 11 and V beta 8 loci TCR beta haplotypes have been identified in five families in which the probands have insulin-dependent diabetes mellitus (IDDM). An extremely rare haplotype, marked by the higher molecular weight BamH1 allele (H, H) at each of V beta 11 and V beta 8, was found in the DR4+ DR3- probands of two families (P = 0.004). In three families in which the probands had DR3, the more common TCR haplotype LH (V beta 11, V beta 8) was found. Taken together, these data confirm that linkage disequilibrium does exist in the TCR beta locus, at least in some regions, and suggest that detailed analysis of the relationship between TCR V beta haplotypes and HLA is warranted since these RFLPs may be markers for important allelic V gene sequence variations.

Homozygous parent affected sib pair method for detecting disease predisposing variants: application to insulin dependent diabetes mellitus

For complex genetic diseases involving incomplete penetrance, genetic heterogeneity, and multiple disease genes, it is often difficult to determine the molecular variant(s) responsible for the disease pathogenesis. Linkage and association studies may help identify genetic regions and molecular variants suspected of being directly responsible for disease predisposition or protection, but, especially for complex diseases, they are less useful for determining when a predisposing molecular variant has been identified. In this paper, we expand upon the simple concept that if a genetic factor predisposing to disease has been fully identified, then a parent homozygous for this factor should transmit either of his/her copies at random to any affected children. Closely linked markers are used to determine identity by descent values in affected sib pairs from a parent homozygous for a putative disease predisposing factor. The expected deviation of haplotype sharing from 50%, when not all haplotypes carrying this factor are in fact equally predisposing, has been algebraically determined for a single locus general disease model. Equations to determine expected sharing for multiple disease alleles or multiple disease locus models have been formulated. The recessive case is in practice limiting and therefore can be used to estimate the maximum proportion of putative susceptibility haplotypes which are in fact predisposing to disease when the mode of inheritance of a disease is unknown. This method has been applied to 27 DR3/DR3 parents and 50 DR4/DR4 parents who have at least 2 children affected with insulin dependent diabetes mellitus (IDDM). The transmission of both DR3 and DR4 haplotypes is statistically different from 50% (P < 0.05 and P < 0.001, respectively). An upper estimate for the proportion of DR3 haplotypes associated with a high IDDM susceptibility is 49%, and for DR4 haplotypes 38%. Our results show that the joint presence of non-Asp at DQ beta position 57 and Arg at DQ alpha position 52, which has been proposed as a strong IDDM predisposing factor, is insufficient to explain the HLA component of IDDM predisposition.

The HLA DQB1*0502 allele is neutrally associated with insulin-dependent diabetes mellitus in the Sardinian population

In the Sardinian population a very high incidence of insulin-dependent diabetes mellitus (IDDM) and the lack of HLA-DR2 protective effect due to the high frequency of the A2, Cw7, B17, 3F31, DR2, DQw1 extended haplotype has been reported. This haplotype, carrying a Serine at position 57 of the DQB1*0502 allele, has been previously reported to be underrepresented in patients when compared to controls. In order to provide an explanation for this finding, we defined by RFLP analysis the HLA haplotype of 45 Sardinian IDDM patients and 49 controls. All DR-2DQw1 subjects were molecularly characterized at the HLA DQA and DQB loci. All DR2-positive patients and the vast majority of the DR2-positive controls had the DQB1*0502 allele at the DR2-linked DQB1 locus, with no statistically significant difference between the two groups. All DQA1 genes were the ones expected, with only two exceptions. Nine out of 10 of the DR2-positive patients were compound heterozygotes for DQB1*0201/DQB1*0502 alleles; only this allele combination was significantly increased (p less than 0.0003). Our data suggests that a) the DQB1*0502 allele is neutral for IDDM development and b) the susceptibility to IDDM in our DR2-positive patients is related to the compound heterozygous state between the neutral DQA1*0102/DQB1*0502 and the susceptibility DQA1*0501/DQB1*0201 alleles.

Relationship of HLA to schizophrenia not supported in multiplex families

The role of the human histocompatibility complex (HLA) in the pathogenesis of schizophrenia has been suggested in previous reports. We conducted a genetic study in 33 new families. Our linkage analysis, which used the affected sib-pair method, did not provide evidence for nonrandom assortment. Moreover, the results of an association study using the "haplotype relative risk" method failed to confirm the positive association between HLA A9 and schizophrenia. Taken together, our data did not support any relationship of HLA type to schizophrenia.

HLA-DR and the 5' insulin gene polymorphism in insulin-dependent diabetes

While the human leukocyte antigen (HLA) region provides the major susceptibility for insulin-dependent (type I) diabetes mellitus (IDDM), other (non-HLA) genes must also play a role. Population studies have shown an increased frequency of small insertions (class I alleles) 5' to the insulin gene in individuals with IDDM, suggesting that this region may account for part, if not all, of the non-HLA genetic predisposition. However, no data are available as to whether the relation of the insulin gene polymorphism is to a DR-defined subset of IDDM or with all of IDDM. To test the hypothesis that specific combinations of HLA and insulin gene polymorphism alleles may interact in providing susceptibility for IDDM, HLA-DR and 5' insulin gene insertion size have been determined in 300 individuals with IDDM. The frequency of class 1 insulin gene alleles in the entire sample is 0.79 and the frequency of class 3 alleles (large inserts) is 0.20. The frequencies of class 1 alleles were equal across all DR classes: 0.79 in the DR3/X IDDM subjects, 0.80 in the DR4/X, 0.79 in the DR3/4, and 0.78 in those with DRX/X. Additionally, the frequencies of class 1/1 homozygotes and 1/3 heterozygotes were similar between HLA-DR types. These results suggest that the HLA region and the region 5' to the insulin gene provide independent and nonsynergistic genetic risks for IDDM.

Insulin-IGF2 region on chromosome 11p encodes a gene implicated in HLA-DR4-dependent diabetes susceptibility

A class of alleles at the VNTR (variable number of tandem repeat) locus in the 5' region of the insulin gene (INS) on chromosome 11p is associated with increased risk of insulin-dependent diabetes mellitus (IDDM), but family studies have failed to demonstrate linkage. INS is thought to contribute to IDDM susceptibility but this view has been difficult to reconcile with the lack of linkage evidence. We thus investigated polymorphisms of INS and neighbouring loci in random diabetics, IDDM multiplex families and controls. HLA-DR4-positive diabetics showed an increased risk associated with common variants at polymorphic sites in a 19-kilobase segment spanned by the 5' INS VNTR and the third intron of the gene for insulin-like growth factor II (IGF2). As INS is the major candidate gene from this region, diabetic and control sequence were compared to identify all INS polymorphisms that could contribute to disease susceptibility. In multiplex families the IDDM-associated alleles were transmitted preferentially to HLA-DR4-positive diabetic offspring from heterozygous parents. The effect was strongest in paternal meioses, suggesting a possible role for maternal imprinting. Our results strongly support the existence of a gene or genes affecting HLA-DR4 IDDM susceptibility which is located in a 19-kilobase region of INS-IGF2. Our results also suggest new ways to map susceptibility loci in other common diseases.

Non-HLA region genes in insulin dependent diabetes mellitus

The focus of this chapter is on the contribution of genes outside the HLA region to insulin dependent diabetes mellitus (IDDM) susceptibility. We review laboratory evidence for such genes from published studies and also present unpublished data from our recent research. The existence of genes predisposing to IDDM in the region of the insulin (INS) gene now appears established. Association analysis has demonstrated an increased frequency of class 1 alleles of the 5' INS polymorphism in diabetics compared with controls, and a new method of analysis (AFBAC) has shown that this association is not an artefact of population stratification. Interestingly, the effect of INS region susceptibility on IDDM cannot be detected by linkage analysis, suggesting that if a genetic marker locus is close to a disease susceptibility locus, association analysis may be a more sensitive method than linkage analysis for detecting the susceptibility locus. There is no convincing evidence that genes in the T cell receptor beta chain (TCRB) or alpha chain (TCRA) regions influence predisposition to IDDM, either directly, or indirectly through interaction with HLA region genes. However, we present new evidence for interaction between TCRB and immunoglobulin heavy chain (Gm) region genes in IDDM: diabetics who are positive for the IgG2 allotype G2m(23) have significantly different frequencies of a TCRB restriction fragment length polymorphism (RFLP) than those who are negative for the allotype. Gm region genes also appear to have indirect effects on IDDM susceptibility through interaction with HLA and INS region genes: DR3/4 and non-DR3/4 diabetics have significantly different frequencies of G2m(23), and INS1/1 and non-INS1/1 diabetics also have significantly different frequencies of this allotype. To our knowledge, there are no other studies of Gm-TCRB or Gm-INS interaction in IDDM susceptibility. Evidence for Gm-HLA interaction in IDDM has been published by several other groups of investigators, however the specific phenotypic interaction effects reported have differed. Nevertheless, pooled data from three studies of Gm/HLA haplotype segregation in affected sib pairs shows significantly increased sharing of Gm haplotypes in affected pairs who share both HLA haplotypes. The biological mechanisms underlying the direct (HLA, INS) and indirect (Gm-TCRB, Gm-HLA, Gm-INS) effects of these genetic regions on IDDM susceptibility remain to be elucidated.

Genetic analysis of autoimmune type 1 diabetes mellitus in mice

Two genes, Idd-3 and Idd-4, that influence the onset of autoimmune type 1 diabetes in the nonobese diabetic mouse have been located on chromosomes 3 and 11, outside the chromosome 17 major histocompatibility complex. A genetic map of the mouse genome, analysed using the polymerase chain reaction, has been assembled specifically for the study. On the basis of comparative maps of the mouse and human genomes, the homologue of Idd-3 may reside on human chromosomes 1 or 4 and Idd-4 on chromosome 17.