Effect of genetic heterogeneity and assortative mating on linkage analysis: a simulation study.

Linkage studies of complex genetic traits raise questions about the effects of genetic heterogeneity and assortative mating on linkage analysis. To further understand these problems, I have simulated and analyzed family data for a complex genetic disease in which disease phenotype is determined by two unlinked disease loci. Two models were studied, a two-locus threshold model and a two-locus heterogeneity model. Information was generated for a marker locus linked to one of the disease-defining loci. Random-mating and assortative-mating samples were generated. Linkage analysis was then carried out by use of standard methods, under the assumptions of a single-locus disease trait and a random-mating population. Results were compared with those from analysis of a single-locus homogeneous trait in samples with the same levels of assortative mating as those considered for the two-locus traits. The results show that (1) introduction of assortative mating does not, in itself, markedly affect the estimate of the recombination fraction; (2) the power of the analysis, reflected in the LOD scores, is somewhat lower with assortative rather than random mating. Loss of power is greater with increasing levels of assortative mating; and (3) for a heterogeneous genetic disease, regardless of mating type, heterogeneity analysis permits more accurate estimate of the recombination fraction but may be of limited use in distinguishing which families belong to each homogeneous subset. These simulations also confirmed earlier observations that linkage to a disease "locus" can be detected even if the disease is incorrectly defined as a single-locus (homogeneous) trait, although the estimated recombination fraction will be significantly greater than the true recombination fraction between the linked disease-defining locus and the marker locus.     

The detection of linkage and heterogeneity in nuclear families for complex disorders: one versus two marker loci.

Using exact expected likelihoods, we have computed the average number of phase-unknown nuclear families needed to detect linkage and heterogeneity. We have examined the case of both dominant and recessive inheritance with reduced penetrance and phenocopies. Most of our calculations have been carried out under the assumption that 50% of families are linked to a marker locus. We have varied both the number of offspring per family and the sampling scheme. We have also investigated the increased power when the disease locus is midway between two marker loci 10 cM apart. For recessive inheritance, both linkage and heterogeneity can be detected in clinically feasible sample sizes. For dominant inheritance, linkage can be detected but heterogeneity cannot be detected unless larger sibships (four offspring) are sampled or two linked markers are available. As expected, if penetrance is reduced, sampling families with all sibs affected is most efficient. Our results provide a basis for estimating the amount of resources needed to find genes for complex disorders under conditions of heterogeneity.     

Linkage analysis of "necessary" disease loci versus "susceptibility" loci.

The association of some diseases with specific alleles of certain genetic markers has been difficult to explain. Several explanations have been proposed for the phenomenon of association, e.g. the existence of multiple, interacting genes (epistasis) or a disease locus in linkage disequilibrium with the marker locus. One might suppose that when marker data from families with associated diseases are analyzed for linkage, the existence of the association would assure that linkage will be found, and found at a tight recombination fraction. In fact, however, linkage analyses of some diseases associated with HLA, as well as diseases associated with alleles at other loci located throughout the genome, show significant evidence against linkage, and others show loose linkage, to the puzzlement of many researchers. In part, the puzzlement arises because linkage analysis is ideal for looking for loci that are necessary, even if not sufficient, for disease expression but may be much less useful for finding loci that are neither necessary nor sufficient for disease expression (so-called susceptibility loci). This work explores what happens when one looks for linkage to susceptibility loci. A susceptibility locus in this case means that the allele increases risk but is neither necessary nor sufficient for disease expression. It might be either an allele at the marker locus itself that is increasing susceptibility or an allele at a locus in linkage disequilibrium with the marker. This work uses computer simulation to examine how linkage analyses behave when confronted with data from such a model.(ABSTRACT TRUNCATED AT 250 WORDS)     

An AFLP linkage map for Douglas-fir based upon multiple full-sib families

An amplified fragment length polymorphism (AFLP) linkage map for coastal Douglas-fir (Pseudotsuga menziesii) was constructed from eight full-sib families each consisting of 40 progeny. These families were part of the British Columbia Ministry of Forests second-generation progeny test program and represent typical family sizes used in progeny trials. For map construction, ten primer pairs using EcoRI+3 and MseI+4 were employed to identify and assay AFLP loci that segregated in backcross configurations. A new technique was used to obtain a single recombination rate for each pair of marker loci: for each locus pair, a recombination rate and log-odd value were estimated across all segregating families using a joint maximum likelihood function that considered the full dataset of segregating genotypes. The resulting matrix of recombination rates between all pairs of loci was used to construct an integrated linkage map using JoinMap. The final map consisted of 19 linkage groups spanning 938.6 cM at an average distance of 9.3 cM between markers. The simultaneous integration of data from multiple families may provide an effective way to construct a linkage map, using the genetic resources inherent in most tree improvement programs, where progeny tests of small size are conducted. The statistical property of number of families used is briefly discussed. For our data, at least three to four families greatly increased the chance of obtaining an informative locus in at least one family. Families as small as ten are adequate for closely linked loci (<10 cM), while the size used in our study (40) is adequate for loci within 30 cM.     

Sample-size considerations and strategies for linkage analysis in autosomal recessive disorders.

The opportunity raised by recombinant DNA technology to develop a linkage marker panel that spans the human genome requires cost-efficient strategies for its optimal utilization. Questions arise as to whether it is more cost-effective to convert a dimorphic restriction enzyme marker system into a highly polymorphic system or, instead, to increase the number of families studied, simply using the available marker alleles. The choice is highly dependent on the population available for study, and, therefore, an examination of the informational content of the various family structures is important to obtain the most informative data. To guide such decisions, we have developed tables of the average sample number of families required to detect linkage for autosomal recessive disorders under single backcross and under "fully informative" matings. The latter cross consists of a marker locus with highly polymorphic codominant alleles such that the parental marker genotypes can be uniquely distinguished. The sampling scheme considers families with unaffected parents of known mating types ascertained via affected offspring, for sibship sizes ranging from two to four and various numbers of affected individuals. The sample-size tables, calculated for various values of the recombination fractions and lod scores, may serve as a guide to a more efficient application of the restriction fragment length polymorphism technology to sequential linkage analysis.     

Utility and efficiency of linked marker genes for genetic counseling. II. Identification of linkage phase by offspring phenotypes

For a linked marker locus to be useful for genetic counseling, the counselee must be heterozygous for both disease and marker loci and his or her linkage phase must be known. It is shown that when the phenotypes of the counselee's previous children for the disease and marker loci are known, the linkage phase can often be inferred with a high probability, and thus it is possible to conduct genetic counseling. To evaluate the utility of linked marker genes for genetic counseling, the accuracy of prediction of the risk for a prospective child with a given marker gene to develop the genetic disease and the proportion of families in which a particular marker locus can be used for genetic counseling are studied for X-linked recessive, autosomal dominant, and autosomal recessive diseases. In the case of X-linked genetic diseases, information from children is very useful for determining the linkage phase of the counselee and predicting the genetic disease. In the case of autosomal dominant diseases, not all children are informative, but if the number of children is large, the phenotypes of children are often more informative than the information from grandparents. In the case of autosomal recessive diseases, information from grandparents is usually useless, since they show a normal phenotype for the disease locus. If we use information on the phenotypes of children, however, the linkage phase of the counselee and the risk of a prospective child can be inferred with a high probability. The proportion of informative families depends on the dominance relationship and frequencies of marker alleles, and the number of children. In general, codominant markers are more useful than are dominant markers, and a locus with high heterozygosity is more useful than is a locus with low heterozygosity.     

Utility and efficiency of linked marker genes for genetic counseling. III. Proportion of informative families under linkage disequilibrium.

A marker locus closely linked to a disease locus is often useful for genetic counseling provided that a counselee is heterozygous at both disease and marker loci. Furthermore, the linkage phase of these genes in the counselee must be known. When the linkage between the disease and marker loci is very close, one often finds linkage disequilibrium between the loci. To evaluate the effect of such nonrandom associations on the utility of linked marker genes for genetic counseling, the proportion of informative families is studied for X-linked recessive and autosomal dominant diseases. This proportion is higher for X-linked genes than for autosomal genes, if other factors are the same. In general, codominant markers are more useful than dominant markers. Also, under appropriate conditions, the proportion of informative families is higher when linkage disequilibrium is present. The results obtained in this paper are useful for evaluating the utility of polymorphic restriction endonuclease cleavage sites as markers in genetic counseling.     

Detection of linkage for heterogeneous disorders by using multipoint linkage analysis.

We have compared the efficiency of the lod score test which assumes heterogeneity (lod2) to the standard lod score test which assumes homogeneity (lod1) when three-point linkage analysis is used in successive map intervals. If it is assumed that a gene located midway between two linked marker loci is responsible for a proportion of disease cases, then the lod1 test loses power relative to the lod2 test, as the proportion of linked families decreases, as the flanking markers are more closely linked, and as more map intervals are tested. Moreover, when multipoint analysis is used, linkage for a disease gene is more likely to be incorrectly excluded from a complete and dense linkage map if true genetic heterogeneity is ignored. We thus conclude that, in general, the lod2 linkage test is more efficient for detecting a true linkage when a complete genetic marker map is screened for a heterogeneous disorder.     

Efficient strategies for genomic searching using the affected-pedigree-member method of linkage analysis.

The affected-pedigree-member (APM) method of linkage analysis is a nonparametric statistic that tests for nonrandom cosegregation of a disease and marker loci. The APM statistic is based on the observation that if a marker locus is near a disease-susceptibility locus, then affected individuals within a family should be more similar at the marker locus than is expected by chance. The APM statistic measures marker similarity in terms of identity by state (IBS) of marker alleles; that is, two alleles are IBS if they are the same, regardless of their ancestral origin. Since the APM statistic measures increased marker similarity, it makes no assumptions concerning how the disease is inherited; this can be an advantage when dealing with complex diseases for which the mode of inheritance is difficult to determine. We investigate here the power of the APM statistic to detect linkage in the context of a genomewide search. In such a search, the APM statistic is evaluated at a grid of markers. Then regions with high APM statistics are investigated more thoroughly by typing more markers in the region. Using simulated data, we investigate various search strategies and recommend an optimal search strategy that maximizes the power to detect linkage while minimizing the false-positive rate and number of markers. We determine an optimal series of three increasing cut-points and an independent criterion for significance.     

Inter- and intrafamilial heterogeneity: effective sampling strategies and comparison of analysis methods.

Heterogeneity, both inter- and intrafamilial, represents a serious problem in linkage studies of common complex diseases. In this study we simulated different scenarios with families who phenotypically have identical diseases but who genotypically have two different forms of the disease (both forms genetic). We examined the proportion of families displaying intrafamilial heterogeneity, as a function of mode of inheritance, gene frequency, penetrance, and sampling strategies. Furthermore, we compared two different ways of analyzing linkage in these data sets: a two-locus (2L) analysis versus a one-locus (SL) analysis combined with an admixture test. Data were simulated with tight linkage between one disease locus and a marker locus; the other disease locus was not linked to a marker. Our findings are as follows: (1) In contrast to what has been proposed elsewhere to minimize heterogeneity, sampling only "high-density" pedigrees will increase the proportion of families with intrafamilial heterogeneity, especially when the two forms are relatively close in frequency. (2) When one form is dominant and one is recessive, this sampling strategy will greatly decrease the proportions of families with a recessive form and may therefore make it more difficult to detect linkage to the recessive form. (3) An SL analysis combined with an admixture test achieves about the same lod scores and estimate of the recombination fraction as does a 2L analysis. Also, a 2L analysis of a sample of families with intrafamilial heterogeneity does not perform significantly better than an SL analysis. (4) Bilineal pedigrees have little effect on the mean maximum lod score and mean maximum recombination fraction, and therefore there is little danger that including these families will lead to a false exclusion of linkage.     

Genetic linkage between the loci for phosphohexose isomerase (PHI) and a serum protein (Xk) in horses

Genetic linkage between the equine loci for phosphohexose isomerase (PHI) and serum Xk protein was demonstrated by means of segregation data from three sire families. The recombination frequency was estimated from pooled data to be 0.23 +/- 0.02; a significant heterogeneity between sires for estimates of the recombination frequency was observed. No indication of linkage was detected between Xk and 14 other blood marker loci. Linkage between the Xk locus and the locus for soluble malic enzyme (ME1) has recently been reported in horses. An equine linkage group designated LG IV comprising the three loci ME1, PHI, and Xk has thus been established. The possibility that the linkage between PHI and Xk is homologous to the linkage between the loci for PHI and a serum postalbumin (PO-2) in pigs was discussed.     

Multipoint linkage analysis and heterogeneity testing in 20 X-linked retinitis pigmentosa families

Using multipoint linkage analysis in 20 families segregating for X-linked retinitis pigmentosa (XLRP), the lod scores on a map of eight RFLP loci were obtained. Our results indicate that under the hypothesis of homogeneity the maximal multipoint lod score supports one disease locus located slightly distal to OTC at Xp21.1. Heterogeneity testing for two XLRP loci suggested that a second XLRP locus may be located 8.5 cM proximal to DXS28 at Xp21.3. Further heterogeneity testing for three disease loci failed to detect a third XLRP locus proximal to DXS7 in any of our 20 XLRP families.     

Genetic linkage between the loci for phosphohexose isomerase (PHI) and a serum protein (Xk) in horses

Genetic linkage between the equine loci for phosphohexose isomerase (PHI) and serum Xk protein was demonstrated by means of segregation data from three sire families. The recombination frequency was estimated from pooled data to be 0.23 ± 0.02; a significant heterogeneity between sires for estimates of the recombination frequency was observed. No indication of linkage was detected between Xk and 14 other blood marker loci. Linkage between the Xk locus and the locus for soluble malic enzyme ( ME1 ) has recently been reported in horses. An equine linkage group designated LG IV comprising the three loci ME1, PH1 , and Xk has thus been established. The possibility that the linkage between PH1 and Xk is homologous to the linkage between the loci for PHI and a serum postalbumin (PO-2) in pigs was discussed.     

Evaluation of designs for reference families for livestock linkage mapping experiments

Abstract

The development of dense linkage maps consisting of highly polymorphic loci for livestock species is technically feasible. However, linkage mapping experiments are expensive as they involve many animals and marker typings per animal. To minimize costs of developing linkage maps for livestock species, optimizing designs for mapping studies is necessary. This study provides a general framework for evaluating the efficiency of designs for reference families consisting of two‐ or three‐ generation full‐sib or half‐sib families selected from a segregating population. The influence of number of families, number of offspring per family, family structure (either half‐sib or full‐sib) and marker polymorphism is determined. Evaluation is done for two markers with a recombination rate of .20 and for a marker and a dominant single gene with a recombination rate of .20. Two evaluation criteria are used: expected maximum lod score for detection of linkage and accuracy of an estimated recombination rate defined as probability that the true recombination rate is in an interval around the estimated recombination rate. First, for several designs the contribution of reference families to expected maximum lod score and accuracy is given. Second, the required number of families in a design to obtain a certain value for the evaluation criteria is calculated when number of offspring per family, family structure and marker polymorphism are specified. The required numbers increase when designs are optimized not only for expected maximum lod score but also for accuracy. The required number of animals to map a dominant single gene is very large. Therefore, a set of reference families should be designed for strictly mapping marker loci. Examples illustrate how tabulated results can be generalized to determine the values for a wide range of designs containing two‐ or three‐generation full‐sib or half‐sib families.     

Mapping of facioscapulohumeral muscular dystrophy gene to chromosome 4q35-qter by multipoint linkage analysis and in situ hybridization

We have recently assigned the facioscapulohumeral muscular dystrophy (FSHD) gene to chromome 4 by linkage to the microsatellite marker Mfd 22 (locus D4S171). We now report that D4S139, a VNTR locus, is much more closely linked to FSHD. Two-point linkage analysis between FSHD and D4S139 in nine informative families showed a maximum combined lod score (Zmax) of 17.28 at a recombination fraction theta of 0.027. Multipoint linkage analysis between FSHD and the loci D4S139 and D4S171 resulted in a peak lod score of 20.21 at 2.7 cM from D4S139. Due to the small number of recombinants found with D4S139, the position of the FSHD gene relative to that of D4S139 could not be established with certainty. D4S139 was mapped to chromosome 4q35-qter by in situ hybridization, thus firmly establishing the location of the FSHD gene in the subtelomeric region of chromosome 4q. One small family yielded a negative lod score for D4S139. In the other families no significant evidence for genetic heterogeneity was obtained. Studies of additional markers and new families will improve the map of the FSHD region, reveal possible genetic heterogeneity, and allow better diagnostic reliability.     

The evolution of recombination in a heterogeneous environment

Most models describing the evolution of recombination have focused on the case of a single population, implicitly assuming that all individuals are equally likely to mate and that spatial heterogeneity in selection is absent. In these models, the evolution of recombination is driven by linkage disequilibria generated either by epistatic selection or drift. Models based on epistatic selection show that recombination can be favored if epistasis is negative and weak compared to directional selection and if the recombination modifier locus is tightly linked to the selected loci. In this article, we examine the joint effects of spatial heterogeneity in selection and epistasis on the evolution of recombination. In a model with two patches, each subject to different selection regimes, we consider the cases of mutation-selection and migration-selection balance as well as the spread of beneficial alleles. We find that including spatial heterogeneity extends the range of epistasis over which recombination can be favored. Indeed, recombination can be favored without epistasis, with negative and even with positive epistasis depending on environmental circumstances. The selection pressure acting on recombination-modifier loci is often much stronger with spatial heterogeneity, and even loosely linked modifiers and free linkage may evolve. In each case, predicting whether recombination is favored requires knowledge of both the type of environmental heterogeneity and epistasis, as none of these factors alone is sufficient to predict the outcome.     

Bayesian analysis of linkage between genetic markers and quantitative trait loci. I. Prior knowledge

Summary Prior information on gene effects at individual quantitative trait loci (QTL) and on recombination rates between marker loci and QTL is derived. The prior distribution of QTL gene effects is assumed to be exponential with major effects less likely than minor ones. The prior probability of linkage between a marker and another single locus is a function of the number and length of chromosomes, and of the map function relating recombination rate to genetic distance among loci. The prior probability of linkage between a marker locus and a quantitative trait depends additionally on the number of detectable QTL, which may be determined from total additive genetic variance and minimum detectable QTL effect. The use of this prior information should improve linkage tests and estimates of QTL effects.     

Hereditary spastic paraplegia: LOD-score considerations for confirmation of linkage in a heterogeneous trait.

Hereditary spastic paraplegia (HSP) is a degenerative disorder of the motor system, defined by progressive weakness and spasticity of the lower limbs. HSP may be inherited as an autosomal dominant (AD), autosomal recessive, or an X-linked trait. AD HSP is genetically heterogeneous, and three loci have been identified so far: SPG3 maps to chromosome 14q, SPG4 to 2p, and SPG4a to 15q. We have undertaken linkage analysis with 21 uncomplicated AD families to the three AD HSP loci. We report significant linkage for three of our families to the SPG4 locus and exclude several families by multipoint linkage. We used linkage information from several different research teams to evaluate the statistical probability of linkage to the SPG4 locus for uncomplicated AD HSP families and established the critical LOD-score value necessary for confirmation of linkage to the SPG4 locus from Bayesian statistics. In addition, we calculated the empirical P-values for the LOD scores obtained with all families with computer simulation methods. Power to detect significant linkage, as well as type I error probabilities, were evaluated. This combined analytical approach permitted conclusive linkage analyses on small to medium-size families, under the restrictions of genetic heterogeneity.     

Selection in complex genetic systems IV. Multiple alleles and interactions between two loci

Summary A symmetric viability model for two loci with two alleles at one locus and m alleles at the other is suggested and analyzed. The analysis of the equilibria is complete if the two loci are absolutely linked, while if recombination is allowed the analysis is incomplete. The dynamics of the mode! resemble those of the two locus two allele model, namely that for loose linkage there will be no correlation between the loci and for tight linkage there may be strong correlation. The major caveats to this are: 1. The equilibria stable for tight linkage may belong to an array of different structures dependent on the selection and the number of alleles. 2. If both loci are overdominant in viability, the stable equilibria always contain all alleles segregating in the population; otherwise, the stable equilibria may only be two locus two allele high complementarity equilibria for tight linkage. 3. For intermediate linkage values and special selection values the boundary two locus two allele high complementarity equilibria may be stable simultaneously with the totally polymorphic central point at which there is no association between the loci.Dedicated to the memory of Ove Frydenberg.Research supported in part by a grant from the Danish Natural Science Research Council, a grant from National Science Foundation, U.S.A., and by USPHS grant NIH 10452-09-11.     

The deviation from linkage equilibrium with multiple loci varying in a stepping-stone cline

The balance between the creation of associations between alleles at different loci by immigration and the convergence to linkage equilibrium due to the recombination process is studied in a theoretical model. The geographical structure of the model is a stepping-stone chain of populations linking two genetically constant source populations. The model assumes an arbitrary number of autosomal loci and considers genetic variation (two alleles at each locus) that is not subject to natural selection. The gene frequencies at each locus will then show a linear cline through the stepping-stone chain of populations. The deviation from linkage equilibrium through the stepping-stone cline is characterized by an equation for linear measures that provide the linkage disequilibrium measures for a given set of loci in terms of the gene frequencies and the linkage disequilibria in the source populations and in terms of the linkage disequilibrium measures through the cline for lower numbers of loci. Numerical examples of this iterative solution are given, and it is shown that the build-up of the higher order Bennett-disequilibria through the cline is considerably more pronounced than the build-up of two-locus disequilibria.