Tests for a Disease-susceptibility Locus allowing for an Inbreeding Coefficient (F)

We begin by discussing the false positive test results that arise because of cryptic relatedness and population substructure when testing a disease susceptibility locus. We extend and evaluate the Hardy-Weinberg disequilibrium (HWD) method, allowing for an inbreeding coefficient (F) in a similar way that Devlin and Roeder (1999) allowed for inbreeding in a case-control study. Then we compare the HWD measure and the common direct measure of linkage disequilibrium, both when there is no population substructure (F = 0) and when there is population substructure (F not = 0), for a single marker. The HWD test statistic gives rise to false positives caused by population stratification. These false positives can be controlled by adjusting the test statistic for the amount of variance inflation caused by the inbreeding coefficient (F). The power loss for the HWD test that arises when controlling for population structure is much less than that which arises for the common direct measure of linkage disequilibrium. However, in the multiplicative model, the HWD test has virtually no power even when allowing for non-zero F.     

一个基于熵的精细定位数量性状位点的指数

针对数量性状位点的精细定位,本文采用群体的极端样本,利用稠密的标记位点,通过比较标记的熵和条件熵,给出了一个基于熵的指数。该指数是标记基因和性状位点间连锁不平衡系数的函数,它不依赖于标记基因的频率。该指数对应我们之前提出的数量性状位点精细定位的哈迪-温伯格不平衡(HWD)指数,但在精细定位数量性状位点时,本文提出的指数的效能要高于哈迪-温伯格不平衡(HWD)指数。通过计算机模拟,文章调查了不同遗传参数下该指数的性质。模拟结果表明该指数用作精细定位是有效的。     

一个基于熵的对疾病基因精细定位的指数

针对人类疾病基因的精细定位,本文利用稠密的标记位点,通过比较标记的熵和条件熵,给出了一个基于熵的指数。该指数可以度量标记基因和性状位点间连锁不平衡(LD)程度。该指数的特性是它不依赖于标记基因的频率。同时它对应疾病易感位点(DSL)精细定位的哈迪-温伯格不平衡(HWD)指数。通过计算机模拟,文章调查了不同遗传参数下该指数的性质。模拟结果表明该指数用作疾病易感位点精细定位是有效的。     

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.     

Recombination disequilibrium in ideal and natural populations

Following Hardy-Weinberg disequilibrium (HWD) occurring at a single locus and linkage disequilibrium (LD) between two loci in generations, we here proposed the third genetic disequilibrium in a population: recombination disequilibrium (RD). RD is a measurement of crossover interference among multiple loci in a random mating population. In natural populations besides recombination interference, RD may also be due to selection, mutation, gene conversion, drift and/or migration. Therefore, similarly to LD, RD will also reflect the history of natural selection and mutation. In breeding populations, RD purely results from recombination interference and hence can be used to build or evaluate and correct a linkage map. Practical examples from F₂, testcross and human populations indeed demonstrate that RD is useful for measuring recombination interference between two short intervals and evaluating linkage maps. As with LD, RD will be important for studying genetic mapping, association of haplotypes with disease, plant breading and population history.     

High-density genetic map construction and identification of a locus controlling weeping trait in an ornamental woody plant (Prunus mume Sieb. et Zucc)

High-density genetic map is a valuable tool for fine mapping locus controlling a specific trait especially for perennial woody plants. In this study, we firstly constructed a high-density genetic map of mei (Prunus mume) using SLAF markers, developed by specific locus amplified fragment sequencing (SLAF-seq). The linkage map contains 8,007 markers, with a mean marker distance of 0.195 cM, making it the densest genetic map for the genus Prunus. Though weeping trees are used worldwide as landscape plants, little is known about weeping controlling gene(s) (Pl). To test the utility of the high-density genetic map, we did fine-scale mapping of this important ornamental trait. In total, three statistic methods were performed progressively based on the result of inheritance analysis. Quantitative trait loci (QTL) analysis initially revealed that a locus on linkage group 7 was strongly responsible for weeping trait. Mutmap-like strategy and extreme linkage analysis were then applied to fine map this locus within 1.14 cM. Bioinformatics analysis of the locus identified some candidate genes. The successful localization of weeping trait strongly indicates that the high-density map constructed using SLAF markers is a worthy reference for mapping important traits for woody plants.     

Linkage analysis and association analysis in the presence of linkage using age at onset of COGA alcoholism data

Complex disease mapping usually involves a combination of linkage and association techniques. Linkage analysis can scan the entire genome in a few hundred tests. Association tests may involve an even greater number of tests. However, association tests can localize the susceptibility genes more accurately. Using a recently developed combined linkage and association strategy, we analyzed a subset of the Collaborative Study on the Genetics of Alcoholism (COGA) data for the Genetic Analysis Workshop 14 (GAW14). In this analysis, we first employed linkage analysis based on frailty models that take into account age of onset information to establish which regions along the chromosome are likely to harbor disease susceptibility genes for alcohol dependence. Second, we used an association analysis by exploiting linkage disequilibrium to narrow down the peak regions. We also compare the methods with mean identity-by-descent tests and transmission/disequilibrium tests that do not use age of onset information.     

Localization of a small genomic region associated with elevated ACE

Defining the relationship between multiple polymorphisms in a small genomic region and an underlying quantitative trait locus (QTL) represents a major challenge in human genetics. Pedigree analyses have shown that angiotensin I-converting enzyme (ACE) levels are influenced by a QTL located within or close to the ACE gene and most likely resides in the 3' region of this locus. We genotyped seven polymorphisms spanning 13 kb in the 3' end of ACE in 159 Afro-Caribbean subjects to evaluate the linkage disequilibrium between these sites and to narrow the genomic region associated with an elevated ACE level using a cladistic analysis. The linkage disequilibrium measurement D' and a haplotype tree revealed three distinct haplotype segments, presumably because of recombination. The value of the linkage disequilibrium parameter p(excess) was highest for site 22982, which is located in the middle segment. A series of nested, cladistic analyses confirmed that the other two regions are unlikely to be the ACE-linked QTL and that the variant resides in the middle region. Analyses of the same polymorphisms in 98 unrelated Europeans in the Monitoring Trends and Determinants in Cardiovascular Diseases (MONICA) study resulted in fewer haplotypes than were observed among the Afro-Caribbean subjects, suggesting that populations with greater genetic diversity may be especially informative for fine-scale mapping.     

内蒙古地区蒙古族HLA-A、B、DRB1基因座多态性分析

应用序列特异性寡核苷酸探针反向斑点杂交技术对内蒙古地区蒙古族106名无关健康个体的HLA-A、B和DRB1 基因座进行基因分型, 以研究内蒙古地区蒙古族人群HLA-A、B、DRB1基因座的等位基因及其组成的单倍型频率分布特征。 采用最大数学预期值算法计算HLA基因座的等位基因频率和单倍型频率。106 名内蒙古地区蒙古族个体的HLA-A、B、DRB1基因座分别检出13、29、13个等位基因。高频单倍型分别为 HLA-A*02-B*46 (0.0510); HLA-A*02-B*13(0.0495); HLA-A*02-B*51(0.0442); HLA-B*13-DRB1*07 (0.0555); HLA- B*46-DRB1*09(0.0378); HLA-B*35-DRB1*13(0.03300); HLA-A*02-B*13-DRB1*07(0.033019); HLA-A*02-B*46- DRB1*09(0.031985)。研究表明: 内蒙古地区蒙古族人群HLA基因座的等位基因和单倍型具有较高的遗传多态性。HLA- A*24-B*14, HLA-A*32-B*63在该民族具有极强的连锁不平衡。     

Robust indices of Hardy-Weinberg disequilibrium for QTL fine mapping

Hardy-Weinberg disequilibrium (HWD) measures have been proposed using dense markers to fine map a quantitative trait locus (QTL) to regions < approximately 1 cM. Earlier HWD measures may introduce bias in the fine mapping because they are dependent on marker allele frequencies across loci. Hence, HWD indices that do not depend on marker allele frequencies are desired for fine mapping. Based on our earlier work, here we present four new HWD indices that do not depend on marker allele frequencies. Two are for use when marker allele frequencies in a study population are known, and two are for use when marker allele frequencies in a study population are not known and are only known in the extreme samples. The new measures are a function of the genetic distance between the marker locus and a QTL. Through simulations, we investigated and compared the fine mapping performance of the new HWD measures with that of the earlier ones. Our results show that when marker allele frequencies vary across loci, the new measures presented here are more robust and powerful.     

Linkage disequilibrium patterns vary with chromosomal location: a case study from the von Willebrand factor region.

Linkage disequilibrium analysis has been used as a tool for analyzing marker order and locating disease genes. Under appropriate circumstances, disequilibrium patterns reflect recombination events that have occurred throughout a population's history. As a result, disequilibrium mapping may be useful in genomic regions of < 1 cM where the number of informative meioses needed to detect recombinant individuals within pedigrees is exceptionally high. Its utility for refining target areas for candidate disease genes before initiating chromosomal walks and cloning experiments will be enhanced as the relationship between linkage disequilibrium and physical distance is better understood. To address this issue, we have characterized linkage disequilibrium in a 144-kb region of the von Willebrand factor gene on chromosome 12. Sixty CEPH and 12 von Willebrand disease families were genotyped for five PCR-based markers, which include two microsatellite repeats and three single-base-pair substitutions. Linkage disequilibrium and physical distance between polymorphisms are highly correlated (rm = -.76; P < .05) within this region. None of the five markers showed significant disequilibrium with the von Willebrand disease phenotype. The linkage disequilibrium/physical distance relationship was also analyzed as a function of chromosomal location for this and eight previously characterized regions. This analysis revealed a general trend in which linkage disequilibrium dissipates more rapidly with physical distance in telomeric regions than in centromeric regions. This trend is consistent with higher recombination rates near telomeres.     

HLA determinants in an Australian population of hemochromatosis patients and their families.

The frequencies of different HLA-A and -B alleles in 77 Australian patients with hemochromatosis have been compared with frequencies of HLA alleles not associated with hemochromatosis in 63 of their heterozygous relatives and with published population frequencies. As for all other populations reported, an association of HLA-A3 and HLA-B7 with the disease was found. A weak association with HLA-B12 was also detected. No other significant positive or negative associations with HLA alleles were detected. In addition, HLA-A2 and -B12 were in significant linkage disequilibrium in patients but not in controls, which may indicate a new mutation or recent recombination between HLA-A and hemochromatosis either in our patient group or in the founding population. HLA-A1 and -B8 and HLA-A29 and -B12 were in linkage disequilibrium in controls but not in patients, suggesting that this population is not segregating a hemochromatosis allele on either of these haplotypes. Genetic linkage analysis using the program LIPED showed strong linkage in 23/24 families, most of which had additional HLA alleles (other than A3 and B7) associated with hemochromatosis. This provides evidence for a single hemochromatosis locus, possibly with more than one allele.     

Development and characterization of microsatellite markers for the endangered Chinese yew Taxus chinensis var. mairei (Taxaceae)

We isolated and characterized 11 microsatellite loci for the Chinese yew, Taxus chinensis var. mairei, an endangered tree species in China, by constructing a (CA)(12)-enriched library. The number of alleles per locus ranged from 5 to 10. The observed heterozygosities ranged from 0.2500 to 0.8333 and the expected heterozygosities ranged from 0.5196 to 0.8680. No significant linkage disequilibrium was detected at these loci. However, four loci significantly deviated from Hardy-Weinberg equilibrium. The null alleles were found to be present at locus Tach9 and locus Tach11 by the Micro-checker test (P < 0.001). These polymorphic loci could be employed in research of gene flow and spatial genetic patterns of T. chinensis var. mairei.     

How do homozygous parents affect TDT as a test for association?

The traditional transmission disequilibrium test (TDT) (Spielman et al., 1993) is a powerful test for association only in the presence of linkage. Since allele transmissions from homozygous parents do not carry any information on linkage, the TDT statistic uses data only on heterozygous parents. However, homozygous parents carry information on association between alleles at a marker locus and a disease locus. In this article, we explore whether inclusion of homozygous parents increases the power to detect association. The resultant test statistic follows a chi(2) distribution with 2 degrees of freedom. Monte-Carlo simulations are included to compare the performance of this test with the traditional TDT under different disease models.     

The immunoglobulin heavy-chain variable region in insulin-dependent diabetes mellitus: affected-sib-pair analysis and association studies.

We have analyzed immunoglobulin heavy-chain variable-region (VH) polymorphisms and genetic susceptibility to insulin-dependent diabetes mellitus (IDDM) by using a set of polymorphic loci that span approximately 1,000 kb of the VH region on chromosome 14q32. One hundred one Finnish families with at least two children affected with IDDM were studied. Conventional RFLPs determined by hybridization were used, since no microsatellite repeat markers have been available for this gene region. No evidence for linkage between the VH genes and IDDM could be obtained from haplotype-sharing analysis among the 133 diabetic sib pairs. The frequencies of various VH genotypes were also compared between 101 familial IDDM cases and 114 controls derived from the Finnish background population. The distribution of the genotypes at the VH2-B5 locus was significantly different between these groups (P=.004), the 3.4/3.4 genotype being less common in the IDDM cases. In addition, a different genotype distribution at the VH5-B2 locus was observed in the diabetic subjects (P = .022). When the IDDM cases were stratified by presence or absence of the high-risk HLA-DQB1*0302 allele, no differences in VH genotype frequencies were observed between the 0302-positive and 0302-negative cases. In the transmission test for linkage disequilibrium (TDT), no differences were found between the expected and observed frequencies of the transmitted alleles at the VH2-B5 or VH5-B2 locus. In conclusion, significant differences in VH genotype distributions were observed between the familial IDDM cases and the controls, but the observed associations could not be confirmed by the TDT. Haplotype sharing analysis provided no evidence for genetic linkage between the VH gene region and IDDM.     

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.     

Resampling methods to reduce the selection bias in genetic effect estimation in genome-wide scans

Using the simulated data of Problem 2 for Genetic Analysis Workshop 14 (GAW14), we investigated the ability of three bootstrap-based resampling estimators (a shrinkage, an out-of-sample, and a weighted estimator) to reduce the selection bias for genetic effect estimation in genome-wide linkage scans. For the given marker density in the preliminary genome scans (7 cM for microsatellite and 3 cM for SNP), we found that the two sets of markers produce comparable results in terms of power to detect linkage, localization accuracy, and magnitude of test statistic at the peak location. At the locations detected in the scan, application of the three bootstrap-based estimators substantially reduced the upward selection bias in genetic effect estimation for both true and false positives. The relative effectiveness of the estimators depended on the true genetic effect size and the inherent power to detect it. The shrinkage estimator is recommended when the power to detect the disease locus is low. Otherwise, the weighted estimator is recommended.     

人类复杂遗传疾病QTL分析方法的理论探讨

利用连锁不平衡理论,人类遗传学家已能把影响人类疾病的质量基因定位在小至1cM区域内,有些基因已被克隆出来。罗泽伟等进一步发展统计分析方法检测及估算分子标记与QTL之间的连锁不平衡系数,从而提出了人类复杂遗传病高解析度基因定位的理论策略。以此为基础,进一步探讨了供试群体在双亲基因频率存在差异时检测QTL和检测QTL互作的方法,给出了有关的理论结果。     

Models and tests of linkage and association studies of quantitative trait locus for multi-allele marker Loci

In this paper, we explore models and tests for association and linkage studies of a quantitative trait locus (QTL) linked to a multi-allele marker locus. Based on the difference between an offspring's conditional trait means of receiving and not receiving an allele from a parent at marker locus, we propose three statistics T(m), T(m,row) and T(m,col) to test association or linkage disequilibrium between the marker locus and the QTL. These tests are composite tests, and use the offspring marginal sample means including offspring data of both homozygous and heterozygous parents. For the linkage study, we calculate the offspring's conditional trait mean given the allele transmission status of a heterozygous parent at the marker locus. Based on the difference between the conditional means of a transmitted and a nontransmitted allele from a heterozygous parent, we propose statistics T(parsi), T(satur), T(gen) and T(m,het) to perform composite tests of linkage between the marker locus and the quantitative trait locus in the presence of association. These tests only use the offspring data that are related to the heterozygous parents at the marker locus. T(parsi) is a parsimonious or allele-wise statistic, T(satur) and T(gen )are satured or genotype-wise statistics, and T(m,het) compares the row and column sample means for offspring data of heterozygous parents. After comparing the powers and the sample sizes, we conclude that T(parsi) has higher power than those of the bi-allele tests, T(satur), T(gen), and T(m,het). If there is tight linkage between the marker and the trait locus, T(parsi) is powerful in detecting linkage between the marker and the trait locus in the presence of association. By investigating the goodness-of-fit of T(parsi), we find that T(satur) does not gain much power compared to that of T(parsi). Moreover, T(parsi) takes into account the pattern of the data that is consistent with linkage and linkage disequilibrium. As the number of alleles at the marker locus increases, T(parsi) is very conservative, and can be useful even for sparse data. To illustrate the usefulness and the power of the methods proposed in this paper, we analyze the chromosome 6 data of the Oxford asthma data, Genetic Analysis Workshop 12.     

Isolation and characterization of 12 polymorphic microsatellite markers from ladyfish (Elops saurus Linnaeus)

Ladyfish (Elops saurus Linnaeus) is an economically important marine fish species. 76 microsatellite loci were isolated from an enriched genomic library of Elops saurus. Twelve of these markers were polymorphic in a test population with alleles per locus ranging from three to nine. The number of observed, expected heterozygosity and polymorphism information content (PIC) per locus in 20 individuals ranged from 0.2000 to 1.0000, 0.1897–0.8846, 0.1769–0.8476, respectively. One markers significantly deviated from Hardy-Weinberg equilibrium after Bonferroni correction analysis and there was no significant linkage disequilibrium found between pairs of markers. As a result, 12 microsatellite markers probably should provide sufficient level of genetic diversity to investigate the fine-scale population structure, stock management and enhancement, genetic linkage map construction and molecular marker-assisted breeding in Elops saurus Linnaeus.