Normative genetic profiles of RAAS pathway gene polymorphisms in North Indian and South Indian populations

Population-based genetic association studies, popularly known as case-control studies, have continued to be the most preferred method for deciphering the genetic basis of various complex diseases, even in the post-human genome sequencing era. However, interpopulation differences in allele, genotype, and haplotype frequencies and linkage disequilibrium patterns lead to inconsistent results in candidate gene association studies. Therefore, for any meaningful disease association study, knowledge of the normative genetic background of the baseline population is a prerequisite. In addition, such genetic variation data also provide a ready-made menu of allele frequencies and linkage disequilibrium patterns of various polymorphisms in specific candidate genes in a particular population, which is a useful reference for further genetic association studies. Such genetic variation data are lacking for the Indian population, which represents about one-sixth of the world's population. In the present study we have reported the allele, genotype, and haplotype frequencies, Hardy-Weinberg equilibrium status, and linkage disequilibrium patterns of 12 polymorphisms in six candidate genes from the renin-angiotensin-aldosterone system among Indians. Because of their different history of origin, the Indian population is broadly divided into two subpopulations: North Indians (Caucasian Europeans) and South Indians (Dravidians). Considering this well-documented difference in gene pools, we have presented a comparative account of the normative genetic data of North Indian and South Indian populations with at least four individuals of urban and suburban origin from each of the representative states of northern and southern India.     

Linkage Disequilibrium in Subdivided Populations

The linkage disequilibrium in a subdivided populaton is shown to be equal to the sum of the average linkage disequilibrium for all subpopulations and the covariance between gene frequencies of the loci concerned. Thus, in a subdivided population the linkage disequilibrium may not be 0 even if the linkage disequilibrium in each subpopulation is 0. If a population is divided into two subpopulations between which migration occurs, the asymptotic rate of approach to linkage equilibrium is equal to either r or 2(m(1) + m(2)) - (m(1) + m(2))(2), whichever is smaller, where r is the recombination value and m(1) and m(2) are the proportions of immigrants in subpopulations 1 and 2, respectively. Thus, if migration rate is high compared with recombination value, the change of linkage disequilibrium in subdivided populations is similar to that of a single random mating population. On the other hand, if migration rate is low, the approach to lnkage equilibrium may be retarded in subdivided populations. If isolated populations begin to exchange genes by migration, linkage disequilibrium may increase temporarily even for neutral loci. If overdominant selection operates and the equilibrium gene frequencies are different in the two subpopulations, a permanent linkage disequilibrium may be produced without epistasis in each subpopulation.     

HL-A antigens in North American black families

The polymorphic HL-A histocompatibility system has been studied in North American black families. The family studies show that haplotype frequencies differ between black and white populations. Seven haplotypes (W28,W5; W28,W17; W28, undefined four; W23,W5; W19,W5; undefined LA,W5; and undefined LA, undefined four) were significantly more frequent in blacks than whites, while haplotypes 1,8 and 3,7 were significantly less frequent. Some of these differences may be accounted for by differences in gene frequencies between the two groups; other differences may be explained by linkage disequilibrium in the white population. No significant linkage disequilibrium between the LA and FOUR loci was found in the black population.     

An efficient haplotyping method with DNA pools

Determination of haplotype frequencies (the joint distribution of genetic markers) in large population samples is a powerful tool for association studies. This is due to their greater extent of polymorphism since any two bi-allelic single nucleotide polymorphisms (SNPs) generate a potential four-allele genetic marker. Therefore, a haplotype may capture a given functional polymorphism with higher statistical power than its SNP components. The statistical estimation of haplotype frequencies, usually employed in linkage disequilibrium studies, requires individual genotyping for each SNP in the haplotype, thus making it an expensive process. In this study, we describe a new method for direct measurement of haplotype frequencies in DNA pools by allele-specific, long-range haplotype amplification. The proposed method allows the efficient determination of haplotypes composed of two SNPs in close vicinity (up to 20 kb).     

Linkage disequilibrium in the North American Holstein population

Linkage disequilibrium was estimated using 7119 single nucleotide polymorphism markers across the genome and 200 animals from the North American Holstein cattle population. The analysis of maternally inherited haplotypes revealed strong linkage disequilibrium ( r ²   >   0.8) in genomic regions of ∼50 kb or less. While linkage disequilibrium decays as a function of genomic distance, genomic regions within genes showed greater linkage disequilibrium and greater variation in linkage disequilibrium compared with intergenic regions. Identification of haplotype blocks could characterize the most common haplotypes. Although maximum haplotype block size was over 1 Mb, mean block size was 26–113 kb by various definitions, which was larger than that observed in humans (∼10 kb). Effective population size of the dairy cattle population was estimated from linkage disequilibrium between single nucleotide polymorphism marker pairs in various haplotype ranges. Rapid reduction of effective population size of dairy cattle was inferred from linkage disequilibrium in recent generations. This result implies a loss of genetic diversity because of the high rate of inbreeding and high selection intensity in dairy cattle. The pattern observed in this study indicated linkage disequilibrium in the current dairy cattle population could be exploited to refine mapping resolution. Changes in effective population size during past generations imply a necessity of plans to maintain polymorphism in the Holstein population.     

Population genetics of HLA polymorphism in blood donors of Bremen and Hannover

501 blood donors from Bremen have been typed for HLA-ABC and -DR. The results are compared with HLA data obtained on 474 blood donors from Hannover. The gene frequencies do not differ significantly between these two population samples. Comparisons with population samples from Kiel, Hamburg, Essen, Frankfurt/M., Mainz, Mannheim, Freiburg/Br., Munich and Vienna did also not reveal any remarkable differences concerning the gene frequencies. Analysis of linkage disequilibrium of two-factor and three-factor haplotypes could show that the typical Caucasoid allele combination A1/B8 is not a constituent part of three-factor haplotype combinations. Between the population samples from Bremen and Hannover no marked differences in the distribution of two-factor and three-factor haplotype frequencies could be found.     

An analysis of linkage disequilibrium in the interleukin-1 gene cluster, using a novel grouping method for multiallelic markers.

In population- and family-based association studies, it is useful to have some knowledge of the patterns of linkage disequilibrium that exist between markers in candidate regions. When such studies are carried out with multiallelic markers, it is often convenient to group the alleles into a biallelic system, for analysis. In this study, we specifically examined the interleukin-1 (IL-1) gene cluster on chromosome 2, a region containing candidates for many inflammatory and autoimmune disorders. Data were collected on eight markers, four of which were multiallelic. Using these data, we investigated the effect of three allele-grouping strategies, including a novel method, on the detection of linkage disequilibrium. The novel approach, termed the "delta method," measures the deviation from the expected haplotype frequencies under linkage equilibrium, for each allelic combination. This information is then used to group the alleles, in an attempt to avoid the grouping together of alleles at one locus that are in opposite disequilibrium with the same allele at the second locus. The estimate haplotype frequencies (EH) program was used to estimate haplotype frequencies and the disequilibrium measure. In our data it was found that the delta method compared well with the other two strategies. Using this method, we found that there was a reasonable correlation between disequilibrium and physical distance in the region (r=-.540, P=.001, one-tailed). We also identified a common, eight-locus haplotype of the IL-1 gene cluster.     

Disequilibrium Pattern Analysis. I. Theory

We have developed a method, disequilibrium pattern analysis, for examining the disequilibrium distribution of the entire array of two locus multiallelic haplotypes in a population. It is shown that a selected haplotype will produce a distinct pattern of linkage disequilibrium values for all generations while the selection is acting. This pattern will also presumably be maintained for many generations after the selection event, until the disequilibrium pattern is eventually broken down by genetic drift and recombination. Related haplotypes, sharing an allele with a selected haplotype, assume a value of linkage disequilibrium proportional to the frequency of the unshared allele and have a single negative value of the normalized linkage disequilibrium. The analysis assumes zero linkage disequilibrium for all allelic combinations initially. The same basic results continue to apply if the selection involves a new mutant, the occurrence of which creates linkage disequilibrium for some haplotypes. The disequilibrium pattern predicted under selection is robust with respect to the influence of migration and random genetic drift. This method is applicable to population data having linked polymorphic loci including that determined from protein or DNA sequencing.     

Accuracy of haplotype estimation in a region of low linkage disequilibrium

We compared the accuracy of haplotype inferences at a 6 Mb region on chromosome 7 where significant linkage between a brain oscillation phenotype and a cholinergic muscarinic receptor gene was previously reported. Individual haplotype assignments and haplotype frequencies were estimated using 5, 10, and 14 consecutive Illumina single-nucleotide polymorphisms (SNPs) within the 1-LOD unit support interval of the chromosome 7 linkage peak. Initially, haplotypes were constructed incorporating phase information provided by relatives using the pedigree analysis package MERLIN. Population-based haplotypes were inferred using the haplotype estimation software HAPLO.STATS and PHASE, using unrelated individuals. The 14 SNPs within this region exhibited markedly low linkage disequilibrium, and the average D' estimate between SNPs was 0.18 (range: 0.01-0.97). In comparison to the family-based haplotypes calculated in MERLIN, the computational inferences of individual haplotype assignments were most accurate when considering 5 consecutive SNPs, but decayed dramatically when considering 10 or 14 SNPs in both PHASE and HAPLO.STATS. When comparing the two haplotype inference methods, both PHASE and HAPLO.STATS performed poorly. These analyses underscore the difficulties of haplotype estimation in the presence of low linkage disequilibrium and stress the importance of careful consideration of confidence measures when using estimated haplotype frequencies and individual assignments in biomedical research.     

Association of neuregulin 1 with schizophrenia confirmed in a Scottish population

Recently, we identified neuregulin 1 (NRG1) as a susceptibility gene for schizophrenia in the Icelandic population, by a combined linkage and association approach. Here, we report the first study evaluating the relevance of NRG1 to schizophrenia in a population outside Iceland. Markers representing a core at-risk haplotype found in Icelanders at the 5' end of the NRG1 gene were genotyped in 609 unrelated Scottish patients and 618 unrelated Scottish control individuals. This haplotype consisted of five SNP markers and two microsatellites, which all appear to be in strong linkage disequilibrium. For the Scottish patients and control subjects, haplotype frequencies were estimated by maximum likelihood, using the expectation-maximization algorithm. The frequency of the seven-marker haplotype among the Scottish patients was significantly greater than that among the control subjects (10.2% vs. 5.9%, P=.00031). The estimated risk ratio was 1.8, which is in keeping with our report of unrelated Icelandic patients (2.1). Three of the seven markers in the haplotype gave single-point P values ranging from .000064 to .0021 for the allele contributing to the at-risk haplotype. This direct replication of haplotype association in a second population further implicates NRG1 as a factor that contributes to the etiology of schizophrenia.     

中国汉族人群MASP2基因单核苷酸多态性的连锁不平衡和单体型分析

目的：探讨中国汉族人群甘露糖结合凝集素相关的丝氨酸蛋白酶（MASP2）基因单核苷酸多态性（SNP）的连锁不平衡和单体型,为相关研究提供更为详实的数据。方法：选取30例广州汉族无关个体,用重新测序的方法进行MASP2基因SNP的发掘,构建连锁不平衡（LD）模式,与HapMap公布的其他4个人群数据相比,并选择4个标签SNP（tagSNP）在232例北京汉族个体和291例广州汉族个体中分析MASP2的多态性和单体型。结果和结论：对MASP2的重测序共检测出16个SNP,LD分析显示来自中国不同地区人群的LD模式基本相同,与来自美国和日本的人群存在差异;北京和广州地区汉族人群tagSNP的等位基因和基因型频率分布不存在显著差异。     

A haplotype-based algorithm for multilocus linkage disequilibrium mapping of quantitative trait loci with epistasis

For tightly linked loci, cosegregation may lead to nonrandom associations between alleles in a population. Because of its evolutionary relationship with linkage, this phenomenon is called linkage disequilibrium. Today, linkage disequilibrium-based mapping has become a major focus of recent genome research into mapping complex traits. In this article, we present a new statistical method for mapping quantitative trait loci (QTL) of additive, dominant, and epistatic effects in equilibrium natural populations. Our method is based on haplotype analysis of multilocus linkage disequilibrium and exhibits two significant advantages over current disequilibrium mapping methods. First, we have derived closed-form solutions for estimating the marker-QTL haplotype frequencies within the maximum-likelihood framework implemented by the EM algorithm. The allele frequencies of putative QTL and their linkage disequilibria with the markers are estimated by solving a system of regular equations. This procedure has significantly improved the computational efficiency and the precision of parameter estimation. Second, our method can detect marker-QTL disequilibria of different orders and QTL epistatic interactions of various kinds on the basis of a multilocus analysis. This can not only enhance the precision of parameter estimation, but also make it possible to perform whole-genome association studies. We carried out extensive simulation studies to examine the robustness and statistical performance of our method. The application of the new method was validated using a case study from humans, in which we successfully detected significant QTL affecting human body heights. Finally, we discuss the implications of our method for genome projects and its extension to a broader circumstance. The computer program for the method proposed in this article is available at the webpage http://www.ifasstat.ufl.edu/genome/~LD.     

SimPed: a simulation program to generate haplotype and genotype data for pedigree structures

With the widespread availability of SNP genotype data, there is great interest in analyzing pedigree haplotype data. Intermarker linkage disequilibrium for microsatellite markers is usually low due to their physical distance; however, for dense maps of SNP markers, there can be strong linkage disequilibrium between marker loci. Linkage analysis (parametric and nonparametric) and family-based association studies are currently being carried out using dense maps of SNP marker loci. Monte Carlo methods are often used for both linkage and association studies; however, to date there are no programs available which can generate haplotype and/or genotype data consisting of a large number of loci for pedigree structures. SimPed is a program that quickly generates haplotype and/or genotype data for pedigrees of virtually any size and complexity. Marker data either in linkage disequilibrium or equilibrium can be generated for greater than 20,000 diallelic or multiallelic marker loci. Haplotypes and/or genotypes are generated for pedigree structures using specified genetic map distances and haplotype and/or allele frequencies. The simulated data generated by SimPed is useful for a variety of purposes, including evaluating methods that estimate haplotype frequencies for pedigree data, evaluating type I error due to intermarker linkage disequilibrium and estimating empirical p values for linkage and family-based association studies.     

HLA phenotype and haplotype frequencies in the Cantabria (middle north Spain) population.

The gene and haplotype frequencies of the HLA-A and -B locus antigens were determined in 502 unrelated individuals from Cantabria (middle North Spain). Our results were compared with those reported for other European and Spanish populations. The haplotypes with significant linkage disequilibrium were also analyzed in various Spanish population samples in order to establish possible relationships with geographic situation and historical events.     

Interleukin–10 promoter haplotypes are differently distributed in the Brazilian versus the Dutch population

The frequency of five different single nucleotide polymorphisms of the promoter interleukin-10 (IL-10) gene (-3575, -2849, 2763, -1082, -819) was compared between two healthy populations, one originating from the Netherlands and one from Rio de Janeiro, Brazil. A total of 321 Caucasian Dutch individuals and 293 Brazilians, grouped as Afro-Brazilians and Euro-Brazilians, were genotyped using PCR-RFLP. The frequencies of the genotypes in the Brazilian population were different (P<0.05) from the frequencies in the Dutch population in all but one (-2763) genotype. The comparison of genotype frequencies between Afro- and Euro-Brazilians did not demonstrate any differences. The haplotype combination of the most-distant three polymorphisms showed strong linkage disequilibrium. All eight possible combinations were observed in Brazilians, but only seven in Dutch Caucasians. The haplotype frequencies were also significantly different between Brazilians when compared with Dutch and also between Euro-Brazilians and Dutch. No differences were observed in haplotype frequencies between Afro-Brazilians and Euro-Brazilians. The -3575T/-2849G/-2763C is more frequent, while the AAA haplotype was much less represented in the Brazilian than in the Dutch population. The haplotype TAC, which was described in African-Americans, was observed only in Brazilians, almost exclusively among those of European origin. The results corroborate the data indicating that the Brazilian population exhibits a genetic admixture of Africans, Europeans, and Amerindians, and the data may serve as a background for clinical and immunological studies involving the IL-10 locus.     

Amino Acid Covariation in a Functionally Important Human Immunodeficiency Virus Type 1 Protein Region Is Associated With Population Subdivision

The frequently reported amino acid covariation of the highly polymorphic human immunodeficiency virus type 1 (HIV-1) exterior envelope glycoprotein V3 region has been assumed to reflect fitness epistasis between residues. However, nonrandom association of amino acids, or linkage disequilibrium, has many possible causes, including population subdivision. If the amino acids at a set of sequence sites differ in frequencies between subpopulations, then analysis of the whole population may reveal linkage disequilibrium even if it does not exist in any subpopulation. HIV-1 has a complex population structure, and the effects of this structure on linkage disequilibrium were investigated by estimating within- and among-subpopulation components of variance in linkage disequilibrium. The amino acid covariation previously reported is explained by differences in amino acid frequencies among virus subpopulations in different patients and by nonsystematic disequilibrium among patients. Disequilibrium within patients appears to be entirely due to differences in amino acid frequencies among sampling time points and among chemokine coreceptor usage phenotypes of virus particles, but not source tissues. Positive selection explains differences in allele frequencies among time points and phenotypes, indicating that these differences are adaptive rather than due to genetic drift. However, the absence of a correlation between linkage disequilibrium and phenotype suggests that fitness epistasis is an unlikely cause of disequilibrium. Indeed, when population structure is removed by analyzing sequences from a single time point and phenotype, no disequilibrium is detectable within patients. These results caution against interpreting amino acid covariation and coevolution as evidence for fitness epistasis.     

A global perspective on genetic variation at the ADH genes reveals unusual patterns of linkage disequilibrium and diversity

Variants of different Class I alcohol dehydrogenase (ADH) genes have been shown to be associated with an effect that is protective against alcoholism. Previous work from our laboratory has shown that the two sites showing the association are in linkage disequilibrium and has identified the ADH1B Arg47His site as causative, with the ADH1C Ile349Val site showing association only because of the disequilibrium. Here, we describe an initial study of the nature of linkage disequilibrium and genetic variation, in population samples from different regions of the world, in a larger segment of the ADH cluster (including the three Class I ADH genes and ADH7). Linkage disequilibrium across approximately 40 kb of the Class I ADH cluster is moderate to strong in all population samples that we studied. We observed nominally significant pairwise linkage disequilibrium, in some populations, between the ADH7 site and some Class I ADH sites, at moderate values and at a molecular distance as great as 100 kb. Our data indicate (1) that most ADH-alcoholism association studies have failed to consider many sites in the ADH cluster that may harbor etiologically significant alleles and (2) that the relevance of the various ADH sites will be population dependent. Some individual sites in the Class I ADH cluster show Fst values that are among the highest seen among several dozen unlinked sites that were studied in the same subset of populations. The high Fst values can be attributed to the discrepant frequencies of specific alleles in eastern Asia relative to those in other regions of the world. These alleles are part of a single haplotype that exists at high (>65%) frequency only in the eastern-Asian samples. It seems unlikely that this haplotype, which is rare or unobserved in other populations, reached such high frequency because of random genetic drift alone.     

Effectiveness of computational methods in haplotype prediction

Haplotype analysis has been used for narrowing down the location of disease-susceptibility genes and for investigating many population processes. Computational algorithms have been developed to estimate haplotype frequencies and to predict haplotype phases from genotype data for unrelated individuals. However, the accuracy of such computational methods needs to be evaluated before their applications can be advocated. We have experimentally determined the haplotypes at two loci, the N-acetyltransferase 2 gene ( NAT2, 850 bp, n=81) and a 140-kb region on chromosome X ( n=77), each consisting of five single nucleotide polymorphisms (SNPs). We empirically evaluated and compared the accuracy of the subtraction method, the expectation-maximization (EM) method, and the PHASE method in haplotype frequency estimation and in haplotype phase prediction. Where there was near complete linkage disequilibrium (LD) between SNPs (the NAT2 gene), all three methods provided effective and accurate estimates for haplotype frequencies and individual haplotype phases. For a genomic region in which marked LD was not maintained (the chromosome X locus), the computational methods were adequate in estimating overall haplotype frequencies. However, none of the methods was accurate in predicting individual haplotype phases. The EM and the PHASE methods provided better estimates for overall haplotype frequencies than the subtraction method for both genomic regions.     

Simultaneous detection of linkage disequilibrium and genetic differentiation of subdivided populations

We propose a new method for simultaneously detecting linkage disequilibrium and genetic structure in subdivided populations. Taking subpopulation structure into account with a hierarchical model, we estimate the magnitude of genetic differentiation and linkage disequilibrium in a metapopulation on the basis of geographical samples, rather than decompose a population into a finite number of random-mating subpopulations. We assume that Hardy-Weinberg equilibrium is satisfied in each locality, but do not assume independence between marker loci. Linkage states remain unknown. Genetic differentiation and linkage disequilibrium are expressed as hyperparameters describing the prior distribution of genotypes or haplotypes. We estimate related parameters by maximizing marginal-likelihood functions and detect linkage equilibrium or disequilibrium by the Akaike information criterion. Our empirical Bayesian model analyzes genotype and haplotype frequencies regardless of haploid or diploid data, so it can be applied to most commonly used genetic markers. The performance of our procedure is examined via numerical simulations in comparison with classical procedures. Finally, we analyze isozyme data of ayu, a severely exploited fish species, and single-nucleotide polymorphisms in human ALDH2.