Haplotype tagging for the identification of common disease genes.

Genome-wide linkage disequilibrium (LD) mapping of common disease genes could be more powerful than linkage analysis if the appropriate density of polymorphic markers were known and if the genotyping effort and cost of producing such an LD map could be reduced. Although different metrics that measure the extent of LD have been evaluated, even the most recent studies have not placed significant emphasis on the most informative and cost-effective method of LD mapping-that based on haplotypes. We have scanned 135 kb of DNA from nine genes, genotyped 122 single-nucleotide polymorphisms (SNPs; approximately 184,000 genotypes) and determined the common haplotypes in a minimum of 384 European individuals for each gene. Here we show how knowledge of the common haplotypes and the SNPs that tag them can be used to (i) explain the often complex patterns of LD between adjacent markers, (ii) reduce genotyping significantly (in this case from 122 to 34 SNPs), (iii) scan the common variation of a gene sensitively and comprehensively and (iv) provide key fine-mapping data within regions of strong LD. Our results also indicate that, at least for the genes studied here, the current version of dbSNP would have been of limited utility for LD mapping because many common haplotypes could not be defined. A directed re-sequencing effort of the approximately 10% of the genome in or near genes in the major ethnic groups would aid the systematic evaluation of the common variant model of common disease.

Cloning of a novel member of the low-density lipoprotein receptor family.

A gene encoding a novel transmembrane protein was identified by DNA sequence analysis within the insulin-dependent diabetes mellitus (IDDM) locus IDDM4 on chromosome 11q13. Based on its chromosomal position, this gene is a candidate for conferring susceptibility to diabetes. The gene, termed low-density lipoprotein receptor related protein 5 (LRP5), encodes a protein of 1615 amino acids that contains conserved modules which are characteristic of the low-density lipoprotein (LDL) receptor family. These modules include a putative signal peptide for protein export, four epidermal growth factor (EGF) repeats with associated spacer domains, three LDL-receptor (LDLR) repeats, a single transmembrane spanning domain, and a cytoplasmic domain. The encoded protein has a unique organization of EGF and LDLR repeats; therefore, LRP5 likely represents a new category of the LDLR family. Both human and mouse LRP5 cDNAs have been isolated and the encoded mature proteins are 95% identical, indicating a high degree of evolutionary conservation.

Expression of the type I diabetes-associated gene LRP5 in macrophages, vitamin A system cells, and the Islets of Langerhans suggests multiple potential roles in diabetes.

LRP5 is a novel member of the low-density lipoprotein receptor family that is genetically associated with Type 1 diabetes. As a start to defining the normal function of LRP5 and to generate testable hypotheses of its potential role in Type 1 diabetes pathogenesis, we carried out an extensive expression analysis of this gene at the mRNA and protein levels in normal human, monkey, and mouse, as well as in non-obese diabetic (NOD) mice at several stages of diabetes development. In all species, expression of LRP5 was found in four functionally important cell types: the distributed mononuclear phagocyte system, the islets of Langerhans, vitamin A-metabolizing cells, and CNS neurons. Given the critical role of macrophages in the onset and progression of islet cell destruction in Type 1 diabetes and the hypothesized role of retinoids as modifiers of diabetes progression, these findings suggest that LRP5 may confer Type 1 diabetes risk by altering the normal functioning of one or more of these regulatory systems. Specifically, given that the LRP5 polymorphisms associated with diabetes are in the promoter region of the gene, alterations in LRP5 expression may be responsible for diabetes susceptibility and therefore may be potential targets for therapeutic intervention. (J Histochem Cytochem 48:1357-1368, 2000)

The sequence and gene characterization of a 400-kb candidate region for IDDM4 on chromosome 11q13.

Type 1 diabetes is a complex disorder with interaction of both genetic and environmental factors. One of the loci, IDDM4, has been mapped to chromosome 11q13, with evidence of association to two markers, D11S1917 and H0570polyA. To identify putative candidate genes for IDDM4, we have constructed a 400-kb clone contig in this region and sequenced the clones. We have also sequenced the orthologous DNA from mouse. Previously, we identified a cDNA for the low-density lipoprotein receptor-related protein 5 gene (LRP5) 3 kb distal to H0570polyA. We have now determined the exon-intron structure of this gene. Detailed sequence analysis has identified a further three genes in this region: the CGI-85 gene (previously identified by W.-C. Lin) and two novel genes, C11orf24 and C11orf23. The C11orf24 gene has no known similarity to other genes, and its function is unknown. C11orf23 has similarity to the SIT4 (sporulation-induced transcript 4)-associated protein (SAP) family of yeast proteins, which are involved in regulation of the cell cycle. The full-length C11orf23 cDNA is the first mammalian orthologue of the yeast SAP family to be identified. Identification of these four genes in a 400-kb region of the IDDM4 region underpins our strategy to identify the IDDM4 locus.

Isolation and characterization of LRP6, a novel member of the low density lipoprotein receptor gene family.

A novel member of the low density lipoprotein receptor (LDLR) gene family has been identified and characterized. This gene, termed LDL receptor-related protein 6 (LRP6), encodes a transmembrane protein which has 71% identity and is structurally similar to the protein encoded by LRP5, a proposed candidate gene for type 1 diabetes located on human chromosome 11q13. LRP6 maps to human chromosome 12p11-p13. Mouse Lrp6 encodes a protein that has 98% identity to human LRP6 and maps to chromosome 6. Unlike other members of the LDLR family, LRP6 and LRP5 display a unique pattern of four epidermal growth factor (EGF) and three LDLR repeats in the extracellular domain. The cytoplasmic domain of LRP6 is not similar to other members of the LDLR family, while comparison with LRP5 reveals proline-rich motifs that may mediate protein-protein interactions. Thus, it is likely that LRP6 and LRP5 comprise a new class of the LDLR family.

Cost-effective analysis of candidate genes using htSNPs: a staged approach.

We have previously shown that the selection of haplotype tag single nucleotide polymorphisms (htSNPs) and their statistical analysis in a multi-locus transmission/disequilibrium test (TDT) results in a more cost-effective genotyping strategy in disease association studies of genes by minimising redundancy due to linkage disequilibrium between SNPs. Further savings can be achieved by the use of a two-stage genotyping strategy. This approach is illustrated here in conjunction with the multi-locus TDT in determining whether common alleles of the immune regulatory genes RANK and its ligand TRANCE (RANKL) are associated with type 1 diabetes (T1D). A saving of approximately 75% of potential genotyping reactions could be made with minimal loss of power. There was little evidence from our analysis for association between the TRANCE and RANK genes and T1D in the populations tested.

Transmission of haplotypes of microsatellite markers rather than single marker alleles in the mapping of a putative type 1 diabetes susceptibility gene (IDDM6).

Allelic association methods based on increased transmission of marker alleles will have to be employed for the mapping of complex disease susceptibility genes. However, because the extent of association of single marker alleles with disease is a function of the relative frequency of the allele on disease-associated chromosomes versus non disease-predisposing chromosomes, the most associated marker allele in a region will not necessarily be closest to the disease locus. To overcome this problem we describe a haplotype-based approach developed for mapping of the putative type 1 diabetes susceptibility gene IDDM6. Ten microsatellite markers spanning a 550 kb segment of chromosome 18q21 in the putative IDDM6 region were genotyped in 1708 type 1 diabetic Caucasian families from seven countries. The most likely ancestral diabetogenic chromosome was reconstructed in a stepwise fashion by analysing linkage disequilibrium between a previously defined haplotype of three adjacent markers and the next marker along the chromosome. A plot of transmission from heterozygous parents to affected offspring of single marker alleles present on the ancestral chromosome versus the physical distance between them, was compared with a plot of transmission of haplotypes of groups of three adjacent markers. Analysing transmission of haplotypes largely negated apparent decreases in transmission of single marker alleles. Peak support for association of the D18S487 region with IDDM6 is P = 0.0002 (corrected P = 0.01). The results also demonstrate the utility of polymorphic microsatellite markers to trace and delineate extended and presumably ancient haplotypes in the analysis of common disease and in the search for identical-by-descent chromosome regions that carry an aetiological variant.

Testing the possible negative association of type 1 diabetes and atopic disease by analysis of the interleukin 4 receptor gene.

Variations in the interleukin 4 receptor A (IL4RA) gene have been reported to be associated with atopy, asthma, and allergy, which may occur less frequently in subjects with type 1 diabetes (T1D). Since atopy shows a humoral immune reactivity pattern, and T1D results from a cellular (T lymphocyte) response, we hypothesised that alleles predisposing to atopy could be protective for T1D and transmitted less often than the expected 50% from heterozygous parents to offspring with T1D. We genotyped seven exonic single nucleotide polymorphisms (SNPs) and the -3223 C>T SNP in the putative promoter region of IL4RA in up to 3475 T1D families, including 1244 Finnish T1D families. Only the -3223 C>T SNP showed evidence of negative association (P=0.014). There was some evidence for an interaction between -3233 C>T and the T1D locus IDDM2 in the insulin gene region (P=0.001 in the combined and P=0.02 in the Finnish data set). We, therefore, cannot rule out a genetic effect of IL4RA in T1D, but it is not a major one.

Fine mapping of the diabetes-susceptibility locus, IDDM4, on chromosome 11q13.

Genomewide linkage studies of type 1 diabetes (or insulin-dependent diabetes mellitus [IDDM]) indicate that several unlinked susceptibility loci can explain the clustering of the disease in families. One such locus has been mapped to chromosome 11q13 (IDDM4). In the present report we have analyzed 707 affected sib pairs, obtaining a peak multipoint maximum LOD score (MLS) of 2.7 (lambda(s)=1.09) with linkage (MLS>=0.7) extending over a 15-cM region. The problem is, therefore, to fine map the locus to permit structural analysis of positional candidate genes. In a two-stage approach, we first scanned the 15-cM linked region for increased or decreased transmission, from heterozygous parents to affected siblings in 340 families, of the three most common alleles of each of 12 microsatellite loci. One of the 36 alleles showed decreased transmission (50% expected, 45.1% observed [P=.02, corrected P=.72]) at marker D11S1917. Analysis of an additional 1,702 families provided further support for negative transmission (48%) of D11S1917 allele 3 to affected offspring and positive transmission (55%) to unaffected siblings (test of heterogeneity P=3x10-4, corrected P=. 01]). A second polymorphic marker, H0570polyA, was isolated from a cosmid clone containing D11S1917, and genotyping of 2,042 families revealed strong linkage disequilibrium between the two markers (15 kb apart), with a specific haplotype, D11S1917*03-H0570polyA*02, showing decreased transmission (46.4%) to affected offspring and increased transmission (56.6%) to unaffected siblings (test of heterogeneity P=1.5x10-6, corrected P=4.3x10-4). These results not only provide sufficient justification for analysis of the gene content of the D11S1917 region for positional candidates but also show that, in the mapping of genes for common multifactorial diseases, analysis of both affected and unaffected siblings is of value and that both predisposing and nonpredisposing alleles should be anticipated.