Comparison of genome screens for two independent cohorts provides replication of suggestive linkage of bone mineral density to 3p21 and 1p36

Low bone mineral density (BMD) is a major risk factor for osteoporotic fracture. Studies of BMD in families and twins have shown that this trait is under strong genetic control. To identify regions of the genome that contain quantitative trait loci (QTL) for BMD, we performed independent genomewide screens, using two complementary study designs. We analyzed unselected nonidentical twin pairs (1,094 pedigrees) and highly selected, extremely discordant or concordant (EDAC) sib pairs (254 pedigrees). Nonparametric multipoint linkage (NPL) analyses were undertaken for lumbar spine and total-hip BMD in both cohorts and for whole-body BMD in the unselected twin pairs. The maximum evidence of linkage in the unselected twins (spine BMD, LOD 2.7) and the EDAC pedigrees (spine BMD, LOD 2.1) was observed at chromosome 3p21 (76 cM and 69 cM, respectively). These combined data indicate the presence, in this region, of a gene that regulates BMD. Furthermore, evidence of linkage in the twin cohort (whole-body BMD; LOD 2.4) at chromosome 1p36 (17 cM) supports previous findings of suggestive linkage to BMD in the region. Weaker evidence of linkage (LOD 1.0-2.3) in either cohort, but not both, indicates the locality of additional QTLs. These studies validate the use, in linkage analysis, of large cohorts of unselected twins phenotyped for multiple traits, and they highlight the importance of conducting genome scans in replicate populations as a prelude to positional cloning and gene discovery.     

A unified sampling approach for multipoint analysis of qualitative and quantitative traits in sib pairs

Recent advances in molecular biology have enhanced the opportunity to conduct multipoint mapping for complex diseases. Concurrently, one sees a growing interest in the use of quantitative traits in linkage studies. Here, we present a multipoint sib-pair approach to locate the map position (tau) of a trait locus that controls the observed phenotype (qualitative or quantitative), along with a measure of statistical uncertainty. This method builds on a parametric representation for the expected identical-by-descent statistic at an arbitrary locus, conditional on an event reflecting the sampling scheme, such as affected sib pairs, for qualitative traits, or extreme discordant (ED) sib pairs, for quantitative traits. Our results suggest that the variance about tau&d4;, the estimator of tau, can be reduced by as much as 60%-70% by reducing the length of intervals between markers by one half. For quantitative traits, we examine the precision gain (measured by the variance reduction in tau&d4;) by genotyping extremely concordant (EC) sib pairs and including them along with ED sib pairs in the statistical analysis. The precision gain depends heavily on the residual correlation of the quantitative trait for sib pairs but considerably less on the allele frequency and exact genetic mechanism. Since complex traits involve multiple loci and, hence, the residual correlation cannot be ignored, our finding strongly suggests that one should incorporate EC sib pairs along with ED sib pairs, in both design and analysis. Finally, we empirically establish a simple linear relationship between the magnitude of precision gain and the ratio of the number of ED pairs to the number of EC pairs. This relationship allows investigators to address issues of cost effectiveness that are due to the need for phenotyping and genotyping subjects.     

Molecular cytogenetic maps of sorghum linkage groups 2 and 8

To integrate genetic, physical, and cytological perspectives of the Sorghum bicolor genome, we selected 40 landed bacterial artificial chromosome (BAC) clones that contain different linkage map markers, 21 from linkage group 2 (LG-02) and 19 from linkage group 8 (LG-08). Multi-BAC probe cocktails were constructed for each chromosome from the landed BACs, which were also preevaluated for FISH signal quality, relative position, and collective chromosome coverage. Comparison to the corresponding linkage map revealed full concordance of locus order between cytological and prior segregation analyses. The pericentromeric heterochromatin constituted a large quasi-uniform block in each bivalent and was especially large in the bivalent corresponding to LG-08. Centromere positions in LG-02 and LG-08 were progressively delimited using FISH to identify landed BACs for which the FISH signals visibly flanked the centromere. Alignment of linkage and cytological maps revealed that pericentromeric heterochromatin of these sorghum chromosomes is largely devoid of recombination, which is mostly relegated to the more distal regions, which are largely euchromatic. This suggests that the sorghum genome is thus even more amenable to physical mapping of genes and positional cloning than the C-value alone might suggest. As a prelude to positional cloning of the fertility restorer, Rf1, FISH of BAC clones flanking the Rf1 locus was used to delimit the chromosomal position of the gene. FISH of BACs that contain the most proximal linkage markers enabled localization of Rf1 to a approximately 0.4-Mbp euchromatic region of LG-08. Cytogenetic analyses of Rf1 and other trait loci will aid in assessing the feasibility of positional cloning and help formulate strategies required for cloning this and other agriculturally critical genes.     

A robust identity-by-descent procedure using affected sib pairs: multipoint mapping for complex diseases

Multipoint linkage analysis is a powerful tool to localize susceptibility genes for complex diseases. However, the conventional lod score method relies critically on the correct specification of mode of inheritance for accurate estimation of gene position. On the other hand, allele-sharing methods, as currently practiced, are designed to test the null hypothesis of no linkage rather than estimate the location of the susceptibility gene(s). In this paper, we propose an identity-by-descent (IBD)-based procedure to estimate the location of an unobserved susceptibility gene within a chromosomal region framed by multiple markers. Here we deal with the practical situation where some of the markers might not be fully informative. Rather the IBD statistic at an arbitrary within the region is imputed using the multipoint marker information. The method is robust in that no assumption about the genetic mechanism is required other than that the region contains no more than one susceptibility gene. In particular, this approach builds upon a simple representation for the expected IBD at any arbitrary locus within the region using data from affected sib pairs. With this representation, one can carry out a parametric inference procedure to locate an unobserved susceptibility gene. In addition, here we derive a sample size formula for the number of affected sib pairs needed to detect linkage with multiple markers. Throughout, the proposed method is illustrated through simulated data. We have implemented this method including exploratory and formal model-fitting procedures to locate susceptibility genes, plus sample size and power calculations in a program, GENEFINDER, which will be made available shortly.     

A linkage strategy for detection of human quantitative-trait loci. I. Generalized relative risk ratios and power of sib pairs with extreme trait values.

We generalize the concept of the relative risk ratio (lambda) to the case of quantitative traits, to take into account the various trait outcomes of a relative pair. Formulas are derived to express the expected proportions of genes shared identical by descent by a sib pair, in terms of the generalized lambda's for sib pairs (lambda S), parent-offspring pairs (lambda O), and monozygotic twins (lambda M) and in terms of the recombination fraction, with the assumption of no residual correlations. If residual correlations are nonzero among relative pairs, we assume that they are the same among sib pairs, parent-offspring pairs, and monozygotic twins, and we employ a slightly different definition for the generalized lambda so that the same set of formulas still hold. The power (or, the sample size necessary) to detect quantitative-trait loci (QTLs) by use of extreme sib pairs (ESPs) is shown to be a function of the three generalized lambda's. Since lambda M can be derived by use of values of lambda S and lambda O, estimates of the latter two lambda's will suffice for the analysis of power and the necessary sample sizes of ESPs, for a QTL linkage study.     

Assignment of the muscle-eye-brain disease gene to 1p32-p34 by linkage analysis and homozygosity mapping.

Muscle-eye-brain disease (MEB) is an autosomal recessive disease of unknown etiology characterized by severe mental retardation, ocular abnormalities, congenital muscular dystrophy, and a polymicrogyria-pachygyria-type neuronal migration disorder of the brain. A similar combination of muscle and brain involvement is also seen in Walker-Warburg syndrome (WWS) and Fukuyama congenital muscular dystrophy (FCMD). Whereas the gene underlying FCMD has been mapped and cloned, the genetic location of the WWS gene is still unknown. Here we report the assignment of the MEB gene to chromosome 1p32-p34 by linkage analysis and homozygosity mapping in eight families with 12 affected individuals. After a genomewide search for linkage in four affected sib pairs had pinpointed the assignment to 1p, the MEB locus was more precisely assigned to a 9-cM interval flanked by markers D1S200 proximally and D1S211 distally. Multipoint linkage analysis gave a maximum LOD score of 6.17 at locus D1S2677. These findings provide a starting point for the positional cloning of the disease gene, which may play an important role in muscle function and brain development. It also provides an opportunity to test other congenital muscular dystrophy phenotypes, in particular WWS, for linkage to the same locus.     

Mapping the Flv locus controlling resistance to flaviviruses on mouse chromosome 5.

Genetically determined resistance to flaviviruses in mice is a dominant trait conferred by alleles at a single autosomal locus designated Flv, but no gene products have been associated with this locus and the mechanism of resistance is not well understood. To further characterize this model of genetic resistance, we conducted mapping studies to determine the chromosomal location of Flv. Because of evidence suggesting that the Flv locus is on chromosome 5, three-point backcross linkage analyses were used to define the location of Flv relative to previously assigned chromosome 5 markers. The results confirm the chromosome 5 location of Flv and indicate a map position between the anchor loci rd and Gus-s. The chromosomal localization of Flv is the first step in the production of a detailed linkage map of the Flv region, which may open approaches to positional cloning of the resistance gene.     

Selection strategies for linkage studies using twins.

Genetic linkage analysis for complex diseases offers a major challenge to geneticists. In these complex diseases multiple genetic loci are responsible for the disease and they may vary in the size of their contribution; the effect of any single one of them is likely to be small. In many situations, like in extensive twin registries, trait values have been recorded for a large number of individuals, and preliminary studies have revealed summary measures for those traits, like mean, variance and components of variance, including heritability. Given the small effect size, a random sample of twins will require a prohibitively large sample size. It is well known that selective sampling is far more efficient in terms of genotyping effort. In this paper we derive easy expressions for the information contributed by sib pairs for the detection of linkage to a quantitative trait locus (QTL). We consider random samples as well as samples of sib pairs selected on the basis of their trait values. These expressions can be rapidly computed and do not involve simulation. We extend our results for quantitative traits to dichotomous traits using the concept of a liability threshold model. We present tables with required sample sizes for height, insulin levels and migraine, three of the traits studied in the GenomEUtwin project.     

A simple formula useful for positional cloning.

We derive a formula for the distribution of the length T of the recombination interval containing a target gene and using N gametes in a region where R kilobases correspond to 1 cM. The formula can be used to calculate the number of meiotic events required to narrow a target gene down to a specific interval size and hence should be useful for planning positional cloning experiments. The predictions of this formula agree well with the results from a number of published experiments in Arabidopsis.     

Pseudoautosomal Linkage of Hodgkin Disease

Heritable factors appear to account for much of the risk for Hodgkin disease (HD). There is evidence for an HLA-linked gene, but other predisposing loci remain unaccounted for. The observation of a family coinheriting both HD and Leri-Weill dyschondrosteosis (LWD) suggests that a gene conferring risk for HD resides adjacent to the LWD locus. The gene responsible for LWD, SHOX, localizes to the short-arm pseudoautosomal region (PAR) of the X and Y chromosomes. A unique segregation pattern for PAR-linked genes has been predicted-that affected sibs will tend to be same sex. An excess of sex-concordant affected sib pairs with HD has been noted but has been attributed to an environmental etiology. These two observations-sex concordance in sib pairs with HD and cosegregation of HD and LWD-impelled a test of the hypothesis that there is a PAR-localized gene for HD. By first scoring recombinations dissociating sex from phenotype in individuals from pedigrees with LWD, we determined a male maximum recombination frequency (thetamax) of.405. This places SHOX near the short-arm telomeres of the sex chromosome and supports the prediction that PAR recombination is obligatory for spermatogenesis. By inferring recombinations between HD and sexual phenotype in sib pairs, we predict, for the postulated HD gene, a male thetamax as high as .254, which places it in proximity to SHOX. Morton's nonparametric affected-sib-pair "beta" model was used in the evaluation of linkage between HD and phenotypic sex and gave a LOD score of 2.41. Using this approach, we reevaluated evidence for HLA linkage in HD in haplotyped sib pairs and found a LOD score of 2.00. The resulting beta values indicate that the putative PAR- and HLA-linked loci account for 29% and 40%, respectively, of the heritability of HD in an American population.     

Evidence for a new Graves disease susceptibility locus at chromosome 18q21

Graves disease (GD) is a common autoimmune thyroid disorder that is inherited as a complex multigenic trait. By using a single microsatellite marker at each locus, we screened the type 1 diabetes loci IDDM4, IDDM5, IDDM6, IDDM8, and IDDM10 and the fucosyltransferase-2 locus for linkage in sib pairs with GD. This showed a two-point nonparametric linkage (NPL) score of 1.57 (P=.06) at the IDDM6 marker D18S41, but NPL scores were <1.0 at the other five loci. Thus, the investigation of the IDDM6 locus was extended by genotyping 11 microsatellite markers spanning 48 cM across chromosome 18q12-q22 in 81 sib pairs affected with autoimmune thyroid disease (AITD). Multipoint analysis, designating all AITD sib pairs as affected, showed a peak NPL score of 3.46 (P=.0003), at the marker D18S487. Designation of only GD cases as affected (74 sib pairs) showed a peak NPL score of 3.09 (P=.001). Linkage to this region has been demonstrated in type 1 diabetes (IDDM6), rheumatoid arthritis, and systemic lupus erythematosus, which suggests that this locus may have a role in several forms of autoimmunity.     

Targeted isolation of a designated region of the Bacillus subtilis genome by recombinational transfer

A method for positional cloning of the Bacillus subtilis genome was developed. The method requires a set of two small DNA fragments that flank the region to be copied. A 38-kb segment that carries genes ppsABCDE encoding five enzymes for antibiotic plipastatin synthesis and another genome locus as large as 100 kb including one essential gene were examined for positional cloning. The positional cloning vector for ppsABCDE was constructed using a B. subtilis low-copy-number plasmid that faithfully copied the precise length of the 38-kb DNA in vivo via the recombinational transfer system of this bacterium. Structure of the copied DNA was confirmed by restriction enzyme analyses. Furthermore, the unaltered structure of the 38-kb DNA was demonstrated by complementation of a ppsABCDE deletion mutant.     

Adding further power to the Haseman and Elston method for detecting linkage in larger sibships: weighting sums and differences

Haseman and Elston (H-E) proposed a robust test to detect linkage between a quantitative trait and a genetic marker. In their method the squared sib-pair trait difference is regressed on the estimated proportion of alleles at a locus shared identical by descent by sib pairs. This method has recently been improved by changing the dependent variable from the squared difference to the mean-corrected product of the sib-pair trait values, a significantly positive regression indicating linkage. Because situations arise in which the original test is more powerful, a further improvement of the H-E method occurs when the dependent variable is changed to a weighted average of the squared sib-pair trait difference and the squared sib-pair mean-corrected trait sum. Here we propose an optimal method of performing this weighting for larger sibships, allowing for the correlation between pairs within a sibship. The optimal weights are inversely proportional to the residual variances obtained from the two different regressions based on the squared sib-pair trait differences and the squared sib-pair mean-corrected trait sums, respectively, allowing for correlations among sib pairs. The proposed method is compared with the existing extension of the H-E approach for larger sibships. Control of the type I error probabilities for sibships of any size can be improved by using a generalized estimating equation approach and the robust sandwich estimate of the variance, or a Monte-Carlo permutation test.     

Consanguinity and the sib-pair method: an approach using identity by descent between and within individuals.

To test for linkage between a trait and a marker, one can consider identical marker alleles in related individuals, for instance, sibs. For recessive diseases, it has been shown that some information may be gained from the identity by descent (IBD) of the two alleles of an affected inbred individual at the marker locus. The aim of this paper is to extend the sib-pair method of linkage analysis to the situation of sib pairs sampled from consanguineous populations. This extension takes maximum advantage of the information provided by both the IBD pattern between sibs and allelic identity within each sib of the pair. This is possible through the use of the condensed identity coefficients. Here, we propose a new test of linkage based on a chi2. We compare the performance of this test with that of the classical chi2 test based on the distribution of sib pairs sharing 0, 1, or 2 alleles IBD. For sib pairs from first-cousin matings, the proposed test can better detect the role of a disease-susceptibility (DS) locus. Its power is shown to be greater than that of the classical test, especially for models where the DS allele may be common and incompletely penetrant; that is to say for situations that may be encountered in multifactorial diseases. A study of the impact of inbreeding on the expected proportions of sib pairs sharing 0, 1, or 2 alleles IBD is also performed here. Ignoring inbreeding, when in fact inbreeding exists, increases the rate of type I errors in tests of linkage.     

Linkage analysis for diseases with variable age of onset

We present a method for the multivariate linkage analysis of the age of onset of a disease. The approach allows the incorporation of covariates for the study of gene by environment interactions. It is applicable to general pedigrees. The likelihood of the data is expressed as a function of the number of alleles identical by descent at a marker, the censored ages of onset and disease status, and environmental exposures. In a simulation study, we compare the power to detect linkage under different sampling schemes for either a dominant or recessive trait when approximately 10% of individuals are gene carriers. The majority of the linkage information from a sample of randomly selected sib pairs was retained when the analyses were limited to sibships with one sibling having early-onset disease (<59 years old). Incorporating parental phenotypes could improve the power to detect the gene. When the sample consists of affected sib pairs (ASPs) having variable age of onset, the likelihood ratio (LR) test had higher power than the means (t(2)) test for detecting a locus with a large genetic relative risk (R(g) = 20). However, the power of the two tests was similar when ASPs are selected so that the proband has an early onset of disease. Lastly, the LR test had more power than the t(2) test to detect linkage in the presence of gene by environment interactions.     

The relative efficiency of penetrance estimators for sib pairs

Here we present analytical studies to evaluate the relative efficiency of commonly used penetrance estimators using linkage designs. We investigated three different methods of estimating penetrance using sib pairs: Maximum likehood estimation (MLE) with trait information alone, MLE with both trait and marker information and the MOD score approach. Modeling sib pairs with unknown phase, we evaluated the asymptotic relative efficiency between estimators under either random sampling or single ascertainment for an autosomal dominant or recessive disease. We then provide plots of the asymptotic relative efficiency, enabling researchers to easily determine regions where the MOD score or segregation alone performs with comparable efficiency relative to joint segregation and linkage.     

Limits of resolution of genetic linkage studies: implications for the positional cloning of human disease genes.

Positional cloning studies to identify disease genes are being carried out for many human genetic diseases. Such studies often include a genome-scan linkage analysis to identify the rough chromosomal location of a disease gene, fine structure genetic mapping to define and narrow the chromosomal interval in which the disease gene may be located, and physical mapping and gene identification in the genetically defined interval to clone the disease gene. During the planning of a positional cloning study, it is important to know that, if linkage is found, the genetic interval identified is likely to be sufficiently narrow to be dissected efficiently by methods of physical mapping and gene identification. Thus, we wish to know the limits of resolution of a genetic linkage study. In this paper, I determine for Mendelian diseases the distributions and moments of three measures of linkage resolution: (1) in a set of N chromosomes, the distance between the nearest crossovers that flank a disease locus, (2) the distance between the nearest genetic markers that flank the pair of flanking crossovers after a genome scan, and (3) the distance between the nearest flanking markers after additional randomly placed markers are generated and typed in an identified interval. These results provide explicit sample-size guidelines for future positional cloning studies of Mendelian diseases and make possible a more objective evaluation of whether a proposed positional cloning study is likely to be successful. I also briefly discuss the more difficult problem of linkage resolution for complex genetic diseases.     

Optimal strategies for mapping complex diseases in the presence of multiple loci.

Recent advances in genome technology have led to mapping and subsequent isolation, by positional cloning, of a number of genes for common and/or complex human diseases. It therefore will be possible to utilize information about a known locus in the search for additional, perhaps less penetrant, genes for a particular disease. It is also unclear, under these situations, what the optimal sampling strategy should be. To address these questions, we have calculated the expected LOD score for localizing one locus in a variety of two-locus models of disease, for four different pedigree structures, and under three different scenarios regarding knowledge/testing of one of the two loci. These design considerations are evaluated by use of a cost function that incorporates the costs of ascertaining different family structures, the relative costs of genotyping and mutation testing family members, and the amount of information provided by each family structure and testing scenario. The results indicate that, in most cases, affected sib pairs are a particularly poor strategy, especially when linkage or mutation data are available at the known locus. We also demonstrate that prescreening the sample of families for mutations at known susceptibility loci is, in general, a cost-effective strategy.     

70个水稻微卫星标记染色体位置的更正

微卫星标记(SSR)因其操作简单和稳定可靠的特点而成为一种重要的分子标记,被广泛应用于遗传作图和种质鉴定等方面。但其在染色体上位置的正确性将直接影响到基因定位的正确性和后续研究的方向。利用美国国家生物信息技术中心(NCBI)网站的Blast程序,将2740个SSR标记的前后引物序列与水稻粳稻品种日本晴基因组进行比对,共发现70个标记位于另一条染色体,对这70个标记重新锚定的染色体进行了更正。这将有助于今后水稻分子标记遗传连锁图的正确构建。