Linkage disequilibrium and haplotype diversity in the genes of the renin-angiotensin system: findings from the family blood pressure program

Association studies of candidate genes with complex traits have generally used one or a few single nucleotide polymorphisms (SNPs), although variation in the extent of linkage disequilibrium (LD) within genes markedly influences the sensitivity and precision of association studies. The extent of LD and the underlying haplotype structure for most candidate genes are still unavailable. We sampled 193 blacks (African-Americans) and 160 whites (European-Americans) and estimated the intragenic LD and the haplotype structure in four genes of the renin-angiotensin system. We genotyped 25 SNPs, with all but one of the pairs spaced between 1 and 20 kb, thus providing resolution at small scale. The pattern of LD within a gene was very heterogeneous. Using a robust method to define haplotype blocks, blocks of limited haplotype diversity were identified at each locus; between these blocks, LD was lost owing to the history of recombination events. As anticipated, there was less LD among blacks, the number of haplotypes was substantially larger, and shorter haplotype segments were found, compared with whites. These findings have implications for candidate-gene association studies and indicate that variation between populations of European and African origin in haplotype diversity is characteristic of most genes.     

Haplotype architecture of the norepinephrine transporter gene <Emphasis Type="Italic">SLC6A2</Emphasis> in four populations

The norepinephrine transporter (NET) regulates levels of monoamine neurotransmitters integral to a variety of behaviors and autonomic functions. Two SLC6A2 polymorphisms have been used in genetic association studies, generating intriguing but nondefinitive results on traits such as hypertension and mood. One of these SLC6A2 variants is functional but rare. The other is common but not informative over the entire 48 kb SLC6A2 region and is insufficient to capture the functional diversity potentially contained within any SLC6A2 region. To elucidate SLC6A2 haplotype structure and define markers sufficient to capture haplotype diversity within detected haplotype blocks, 26 single-nucleotide polymorphisms (SNPs) were genotyped in 384 individuals evenly divided across Finnish Caucasian, US Caucasian, Plains American Indian, and African American populations. Three conserved blocks, 13.6, 12.5, and 25 kb in size and showing little evidence for historical recombination were observed in all populations. Haplotype diversity in block 1 and numbers of common haplotypes were highest in African Americans, among whom 5–6 optimal markers were sufficient to maximize diversity of each block. For other populations, 2–3 markers/block sufficed, but the optimal markers differed across populations. The SLC6A2 haplotype map and 25-marker panel (excluding the monomorphic one) is a comprehensive tool for genetic linkage studies on phenotypes related to NET function.Parts of this work were presented at the 53rd Annual Meeting of American Society of Human Genetics, Los Angeles, CA, USA, November 2003.     

Absence of strong linkage disequilibrium between odorant receptor alleles and the major histocompatibility complex

The odorant receptor (OR) genes constitute the largest gene family among vertebrates. While over 800 loci are present in the human genome, their allele diversity is still poorly characterized. It has been hypothesized that the products of OR genes can be relevant in the reproductive context, thereby interacting with products of genes of the major histocompatibility complex (MHC). Here we investigated the genetic diversity of the OR2W6P, OR2B8P, OR1F12 and OR12D2 genes, in order to define haplotypes and haplotype frequencies. We measured levels of linkage disequilibrium (LD) between these OR genes and the MHC genes HLA-A, HLA-B and HLA-DRB1. This was accomplished through the assessment of 30 single nucleotide polymorphisms (SNPs) in samples from 61 family trios. We characterized 26 alleles among the four OR genes and identified three SNPs that had not yet been reported. Based on our haplotype analysis, LD spanning the OR-HLA region is not very strong, and therefore not enough to enable selection regarding specific HLA-OR haplotypes.     

基于标签单核苷酸多态性单倍型和单倍域的构建及其在关联研究中的应用

人类基因组中单倍型(haplotype)和单倍域(haplotype block)的结构提供了人类进化的宝贵信息,并成为发现人类复杂疾病易感基因的有效策略。一个单倍域可分割成多个具有有限单倍型多样性的离散的区域,代表每个区域结构特征的少量标签单核苷酸多态性(tag single nucleotide polymorphism,tSNP)可使绝大部分单倍型相互区分开来。因此,标签SNP在单倍型和单倍域的构建和关联研究中具有重要地位。构建单倍型和单倍域的方法分为两类,分别是基于大家系中基因分型数据和基于统计学的算法。通过系统回顾几种单倍型和单倍域的构建方法,了解它们在不同的疾病模型或根据不同的分割标准,进行关联研究的检验效能,客观评价每种方法的优、缺点、应用前景及其在关联研究中的应用。随着国际人类基因组单倍型图的完成和单倍型构建统计学运算规则的完善,融合数学、物理学、计算机科学等学科的单倍型构建方法将对人类遗传学、复杂疾病易感基因的定位和克隆鉴定等生命科学的相关领域产生深远的影响。     

Genetic architecture of regulatory variation in Arabidopsis thaliana

Studying the genetic regulation of expression variation is a key method to dissect complex phenotypic traits. To examine the genetic architecture of regulatory variation in Arabidopsis thaliana, we performed genome-wide association (GWA) mapping of gene expression in an F(1) hybrid diversity panel. At a genome-wide false discovery rate (FDR) of 0.2, an associated single nucleotide polymorphism (SNP) explains >38% of trait variation. In comparison with SNPs that are distant from the genes to which they were associated, locally associated SNPs are preferentially found in regions with extended linkage disequilibrium (LD) and have distinct population frequencies of the derived alleles (where Arabidopsis lyrata has the ancestral allele), suggesting that different selective forces are acting. Locally associated SNPs tend to have additive inheritance, whereas distantly associated SNPs are primarily dominant. In contrast to results from mapping of expression quantitative trait loci (eQTL) in linkage studies, we observe extensive allelic heterogeneity for local regulatory loci in our diversity panel. By association mapping of allele-specific expression (ASE), we detect a significant enrichment for cis-acting variation in local regulatory variation. In addition to gene expression variation, association mapping of splicing variation reveals both local and distant genetic regulation for intron and exon level traits. Finally, we identify candidate genes for 59 diverse phenotypic traits that were mapped to eQTL.     

Polymorphisms in multidrug resistance 1 (MDR1) gene are associated with refractory Crohn disease and ulcerative colitis

We used coding and noncoding polymorphisms evenly spaced across the ABCB1/MDR1 gene to perform association analysis in Slovenian patients with inflammatory bowel diseases and to obtain haplotype structure and patterns of linkage disequilibrium (LD) in the MDR1 gene. A disease association study was performed in 307 IBD patients, including 144 patients with ulcerative colitis (UC) and 163 patients with Crohn's disease (CD), and 355 healthy controls. Here we report an association between MDR1 alleles, polymorphisms and haplotypes and refractory CD patients, who do not respond to standard therapy, including patients who develop fistulas. We also report an association with UC and MDR1 polymorphisms in a Slovenian population. Haplotypes significantly associated with diseases were defined by single-nucleotide polymorphisms (SNPs) in exons 12 (1236 C>A), 21(A893S), and 26 (3435 C>T). In addition, two intronic SNPs in LD with the disease haplotype, one in intron 13 (rs2235035) and another in intron 16 (rs1922242), were significantly associated with refractory Crohn (P=0.026, odds ratio (OR) 2.7 and P=0.025, OR 2.8, respectively), as well as with UC (P=0.006, OR 1.8 and P=0.026, OR 1.9, respectively). Our results suggest that MDR1 is a potential target for therapy in refractory CD patients and in patients with UC.     

Haplotype analysis at the alcohol dehydrogenase gene region in New Zealand Māori

Alcohol response is a genetically influenced trait, and there is significant variation in the patterns of alcohol consumption between Māori and Caucasians in New Zealand. Previous studies have found that a variant of the alcohol dehydrogenase (ADH) gene (ADH1B*47His) is associated with protection against alcohol dependence in Māori. Here we extend our investigation of the ADH genes, hypothesising a different haplotype signature in Maori compared to Caucasians. We analysed nine single nucleotide polymorphisms (SNPs) spanning a 500-kb region on chromosome 4q surrounding the ADH1B variant and several other alcohol-metabolising genes (ADH 4, 5, 6, 7). Genotyping was carried out on 47 unrelated Māori individuals, and allele frequencies were compared to the Caucasian population. Large differences in minor allele frequencies were observed between Māori and Caucasian populations for six SNPs (P < 0.01). There was also strong linkage disequilibrium (LD) observed among SNP alleles in Maori indicating the presence of extended ancestral haplotype blocks (P < 0.01). Our results suggest that the Māori population has a different haplotype signature at the ADH gene region compared to Caucasians. These findings probably reflect the unique gene flow history of this genomic region in Maori and should be beneficial for designing future genetic association studies of alcohol-response traits and associated disorders in Polynesians.     

MICA, HLA-B haplotypic variation in five population groups of sub-Saharan African ancestry

The human major histocompatibility complex (MHC) class I chain-related gene A (MICA), located 46 kb centromeric to HLA-B, encodes a stress-inducible protein, which is a ligand for the NKG2D receptor. In addition to its primary role in immune surveillance, data suggest that MICA is involved in the immune response to transplants and in susceptibility to some diseases. In this study, 152 subjects from the Yoruba (n=74), Efik (n=32), and Igbo (n=46) tribes of southern Nigeria, 39 nationwide African-American stem cell donors, and 60 African-American individuals residing in the metropolitan Boston area were studied for MICA, HLA-B allelic variation, haplotypic diversity, and linkage disequilibrium (LD). MICA and HLA-B exhibited a high degree of genetic diversity among the populations studied. In particular, MICA allele and HLA-B-MICA haplotype frequencies and LD in the Efik and Igbo tribes were significantly different from the other study groups. HLA-B and MICA loci demonstrated significant global LD in all five populations (P-values &<0.00001). LD also varied in a haplotype-specific manner. A novel MICA allele was detected in the Boston population. These findings are important from an anthropologic perspective, and will inform future HLA-linked disease association studies in related ethnic groups of African-derived ancestry.     

Heterogeneity of linkage disequilibrium in human genes has implications for association studies of common diseases

Linkage disequilibrium (LD) is the central concept of genetic association studies. Although LD has been shown not to be uniformly distributed across the genome, limited information is available about the characteristics of LD within candidate genes at large. We screened coding and regulatory regions of 50 candidate genes for cardiovascular diseases and identified 228 polymorphisms. The overall sequence diversity was 3.81 +/- 0.31 x 10(-4). Intragenic LD was generally very strong, with 40% of the 464 pairs of polymorphisms exhibiting a complete LD. However, if we consider /D'/ = 0.7 as an arbitrary limit for useful LD in association studies, 26% of the pairs fell below this threshold, half of which being in negative LD, a situation where LD is even more difficult to detect. Non-synonymous coding polymorphisms, which are more likely to have a functional role, were more represented among low-frequency alleles and were more often in complete negative LD with other polymorphisms. This implies that in many situations the power to detect the effect of a non-synonymous polymorphism by measuring a nearby marker might be low. Although intragenic LD was partly a function of physical distance, gene-specific patterns of LD were observed, making it difficult to provide general guidelines for selecting the most useful polymorphisms in association studies. For all these reasons, association studies should concentrate on the overall sequence variation of functionally important regions of candidate genes and not only on a few polymorphisms. The variability of important intergenic regions identified by different approaches including comparative genomics will also have to be assessed.     

Haplotypes vs single marker linkage disequilibrium tests: what do we gain?

The genetic dissection of complex diseases represents a formidable challenge for modern human genetics. Recently, it has been suggested that linkage disequilibrium (LD) based methods will be a powerful approach for delineating complex disease genes. Most proposed LD test statistics search for association between a single marker and a putative trait locus. However, the power of a single marker association test may suffer because LD information contained in flanking markers is ignored. Intuitively, haplotypes (which can be regarded as a collection of ordered markers) may be more powerful than individual, unorganised markers. In this study, we derive the analytical tools based on standard chi-square statistics to directly investigate and compare the power between multilocus haplotypes and single marker LD tests. More specifically, novel formulas are obtained in order to calculate expected haplotype frequencies of unlimited size. This study demonstrates that the use of haplotypes can significantly improve the power and robustness of mapping disease genes. Additionally, we detail how the power of haplotype based association tests are affected by important population genetic parameters such as the genetic distance between markers and disease locus, mode of disease inheritance, age of trait causing mutation, frequency of associated marker allele, and level of initial LD. Finally, published data from the Hereditary Hemochromatosis disease region is used to illustrate the utility of haplotypes.     

Association between ADAM33 polymorphisms and adult asthma in the Japanese population

Background ADAM33, a member of the ADAM (a disintegrin and metalloprotease) family, is a putative asthma susceptibility gene recently identified by positional cloning. It is important to know whether the association exists in ethnically diverse populations. Objective To assess whether genetic functional variants of ADAM33 relate to the susceptibility or some phenotypes in adult patients with bronchial asthma in a Japanese population. Methods We searched for single nucleotide polymorphisms (SNPs) in ADAM33 by PCR‐directed sequencing and identified 48 SNPs. Fourteen SNPs were selected with regard to the LD pattern, and genotyped by Taq‐Man and PCR–RFLP methods. We conducted an association study of ADAM33 with 504 adult asthmatic patients and 651 controls, and haplotype analyses of related variants were performed. Results Significant associations with asthma were found for the SNPs T1 (Met764Thr), T2 (Pro774Ser), S2 and V?3 (with the lowest P‐value for T1, P=0.0015; OR 0.63). We analysed the haplotype using these four polymorphisms, and found a positive association with haplotype CCTG (P=0.0024). Conclusion Our results replicate associations reported recently in other ethnic populations, and suggest that the ADAM33 gene is involved in the development of asthma through genetic polymorphisms.     

High Resolution T2 Association Tests of Complex Diseases Based on Family Data

This paper proposes family based Hotelling's T² tests for high resolution linkage disequilibrium (LD) mapping or association studies of complex diseases. Assume that genotype data of multiple markers or haplotype blocks are available for a sample of nuclear families, in which some offspring are affected. Paired Hotelling's T² test statistics are proposed for a high resolution association study using parents as controls for affected offspring, based on two coding methods: haplotype/allele coding and genotype coding. The paired Hotelling's T² tests take not only the correlation between the haplotype blocks or markers into account, but also take the correlation within each parent‐offspring pair into account. The method extends two sample Hotelling's T² test statistics for population case control association studies, which are not valid for family data due to correlation of genetic data among family members. The validity of the proposed method is justified by rigorous mathematical and statistical proof under the large sample theory. The non‐centrality parameter approximations of the test statistics are calculated for power and sample size calculations. From power comparison and type I error calculations, it is shown that the test statistic based on haplotype/allele coding is advantageous over the test statistic of genotype coding. Analysis using multiple markers may provide higher power than single marker analysis. If only one marker is utilized the power of the test statistic based on haplotype/allele coding is nearly identical to that of 1‐TDT. Moreover, a permutation procedure is provided for data analysis. The method is applied to data from a German asthma family study. The results based on the paired Hotelling's T² statistic tests confirm the previous findings. However, the paired Hotelling's T² tests produce much smaller P‐values than those of the previous study. The permutation tests produce similar results to those of the previous study; moreover, additional marker combinations are shown to be significant by permutation tests. The proposed paired Hotelling's T² statistic tests are potentially powerful in mapping complex diseases. A SAS Macro, Hotel_fam.sas, has been written to implement the method for data analysis.     

Sequence variation, linkage disequilibrium and association with Crohn's disease on chromosome 5q31

Chromosome 5q31 contains a cluster of genes involved in immune response, including a 250 kb risk haplotype associated with Crohn's disease (CD) susceptibility. Recently, two functional variants in SLC22A4 and SLC22A5 (L503F and G-207C), encoding the cation transporters OCTN1 and OCTN2, were proposed as causal variants for CD, but with conflicting genetic evidence regarding their contribution. We investigated this locus by resequencing the coding regions of 10 genes in 24 CD cases and deriving a linkage disequilibrium (LD) map of the 27 single nucleotide polymorphisms (SNPs) detected. Ten SNPs representative of the LD groups observed, were tested for CD association. L503F in SLC22A4 was the only nonsynonymous SNP significantly associated with CD (P=0.003), but was not associated with disease in the absence of other markers of the 250 kb risk haplotype. Two other SNPs, rs11242115 in IRF1 and rs17166050 in RAD50, lying outside the 250 kb risk haplotype, also showed CD association (P=0.019 and P=0.0080, respectively). The RAD50 gene contains a locus control region regulating expression of the Th2 cytokine genes at this locus. Other as yet undiscovered SNPs in this region may therefore modulate gene expression and contribute to the risk of CD, and perhaps of other inflammatory phenotypes.     

Clinical and genetic data on Lafora disease patients of Serbian/Montenegrin origin

Lafora disease (LD) is an autosomal recessive, progressive disorder characterized by myoclonus and seizures, inexorable neurologic deterioration, cognitive decline and poor prognosis. LD is caused by mutations either in the EPM2A or in NHLRC1 genes. Here we report clinical and genetic findings on 14 LD patients from 10 families of Serbian/Montenegrin origin. Molecular diagnostics was performed by sequencing the coding regions of the EPM2A and NHLRC1 genes. In addition, haplotype analysis of the chromosomes carrying the two most frequent mutations (c.1048‐1049delGA and deletion of the whole NHLRC1 gene) using eight different markers flanking the NHLRC1 gene was conducted. We identified one new mutation (c.1028T>C) along with the 3 previously reported mutations (c.1048‐1049delGA, c.990delG, deletion of the whole NHLRC1 gene), all of which were located on the NHLRC1 gene. The two predominant mutations (c.1048‐1049delGA and complete NHLRC1 gene deletion) appear to be founder mutations. In addition to documenting the genetic heterogeneity observed for LD, our study suggests that mutations in the NHLRC1 gene may be a common cause of LD in the Serbian/Montenegrin population, primarily because of a founder effect.     

Patterns of Genetic Variation in the Hypertension Candidate Gene GRK4: Ethnic Variation and Haplotype Structure

Association studies using single nucleotide polymorphisms (SNPs) have the potential to help unravel the genetic basis of hypertension. Nevertheless, to date, association studies of hypertension have yielded ambiguous results. It is becoming clear that such association studies must be interpreted within the context of the genetic structure of the populations being studied, and patterns of variation within specific genomic regions. With this in mind we analyzed genetic variation in the G protein‐coupled receptor kinase 4 (GRK4) gene, a gene whose product has recently been shown to inhibit the dopamine receptor D1 (DRD1) from increasing sodium excretion. We genotyped three previously identified GRK4 SNPs, as well as ten additional SNPs, over 71.6 kb of the GRK4 locus in four populations: African Americans, Asians, Hispanics and Caucasians. Haplotype structure varied among populations, with Hispanics and Caucasians having the most linkage disequilibrium (LD) among SNPs. African Americans had three shorter haplotype blocks, while patterns of markers in the Asian populations demonstrated less LD among markers, a pattern inconsistent with block structure. We observed limited haplotype diversity in each of the four populations, with differing haplotype frequencies among the ethnic groups. We also found substantial evidence for population differentiation, with the largest differences between the African‐American and Asian samples with F_ST values in the upper 90^th percentile when compared to a genome‐wide distribution. However, for all population comparisons, F_ST values decreased sharply in the 3′ region of the gene. This pattern of differentiation among populations is consistent with selection in this part of the gene maintaining similar patterns of variation among otherwise divergent populations. Our results document not only different allele frequencies between populations, but differences in haplotype structure that may be important in evaluating association studies between hypertension and GRK4.     

Genetic variants of the HLA-A, HLA-B and AIF1 loci show independent associations with type 1 diabetes in Norwegian families

The main genetic predisposition to type 1 diabetes (T1D) is known to be conferred by the HLA-DRB1, -DQA1 and -DQB1 genes in the major histocompatibility complex (MHC). Other genetic factors within this complex are known to contribute, but their identity has often been controversial. This picture is shared with several other autoimmune diseases (AIDs). Moreover, as common genetic factors are known to exist between AIDs, associations reported with other AIDs may also be involved in T1D. In this study, we have used these observations in a candidate gene approach to look for additional MHC risk factors in T1D. Using complementary conditional methods (involving conditional logistic regression and family-based haplotype tests) and analyses of linkage disequilibrium (LD) patterns, we confirmed association for alleles of the HLA-A and HLA-B genes and found preliminary evidence for a novel association of a single-nucleotide polymorphism (rs2259571) in the AIF1 gene, independent of the DRB1-DQA1-DQB1 genes and of each other. However, no evidence of independent associations for a number of previously suggested candidate polymorphisms was detected. Our results illustrate the importance of a comprehensive adjustment for LD effects when performing association studies in this complex.     

Natural genetic variation caused by small insertions and deletions in the human genome

Human genetic variation is expected to play a central role in personalized medicine. Yet only a fraction of the natural genetic variation that is harbored by humans has been discovered to date. Here we report almost 2 million small insertions and deletions (INDELs) that range from 1 bp to 10,000 bp in length in the genomes of 79 diverse humans. These variants include 819,363 small INDELs that map to human genes. Small INDELs frequently were found in the coding exons of these genes, and several lines of evidence indicate that such variation is a major determinant of human biological diversity. Microarray-based genotyping experiments revealed several interesting observations regarding the population genetics of small INDEL variation. For example, we found that many of our INDELs had high levels of linkage disequilibrium (LD) with both HapMap SNPs and with high-scoring SNPs from genome-wide association studies. Overall, our study indicates that small INDEL variation is likely to be a key factor underlying inherited traits and diseases in humans.     

小家鼠群体基因组连锁不平衡研究

模式动物小鼠是复杂性状相关疾病遗传研究的重要资源.连锁不平衡(linkage disequilib-rium,LD)是群体基因组遗传的重要信息.如果群体LD程度高,相关连锁区段大,有助于用少量遗传位标来对目标基因进行初步定位;反之,如果群体LD程度低,相关连锁区段小,则利用高密度的遗传位标可以对目标基因进行精细定位.本文介绍了实验小鼠群体和野生小鼠群体相关的LD、与LD相关的部分遗传参数以及利用LD进行相关基因定位研究.并提示实验小鼠和野生小鼠各具优势,都是复杂性状基因定位的重要遗传资源.     

Divergent patterns of linkage disequilibrium and haplotype structure across global populations at the interleukin-13 (IL13) locus

Interleukin-13 (IL-13) is a cytokine involved in Th2 immune response, which plays a role in susceptibility to infection by extracellular parasites as well as complex diseases of the immune system such as asthma and allergies. To determine the pattern of genetic diversity at the IL13 gene, we sequenced 3950 bp encompassing the IL13 gene and its promoter in 264 chromosomes from individuals originating from East and West Africa, Europe, China and South America. Thirty-one single-nucleotide polymorphisms (SNPs) arranged in 88 haplotypes were indentified, including the nonsynonymous substitution Arg130Gln in exon 4, which differs in frequency across ethnic groups. We show that genetic diversity and linkage disequilibrium (LD) are not evenly distributed across the gene and that sites in the 5' and 3' regions of the gene show strong differentiation among continental groups. We observe a divergent pattern of haplotype variation and LD across geographic regions and we identify a set of htSNPs that will be useful for functional genetic association studies of complex disease. We use several statistical tests to distinguish the effects of natural selection and demographic history on patterns of genetic diversity at the IL13 locus.