MCMC multilocus lod scores: application of a new approach

On extended pedigrees with extensive missing data, the calculation of multilocus likelihoods for linkage analysis is often beyond the computational bounds of exact methods. Growing interest therefore surrounds the implementation of Monte Carlo estimation methods. In this paper, we demonstrate the speed and accuracy of a new Markov chain Monte Carlo method for the estimation of linkage likelihoods through an analysis of real data from a study of early-onset Alzheimer's disease. For those data sets where comparison with exact analysis is possible, we achieved up to a 100-fold increase in speed. Our approach is implemented in the program lm_bayes within the framework of the freely available MORGAN 2.6 package for Monte Carlo genetic analysis (http://www.stat.washington.edu/thompson/Genepi/MORGAN/Morgan.shtml).     

Re-localization of Actsk-1 to mouse Chromosome 8, a new region of homology with human Chromosome 1

We present here the genetic mapping of the -skeletal actin locus (Actsk-1) on mouse Chromosome (Chr) 8, on the basis of the PCR analysis of a microsatellite in an interspecific backcross. Linkage and genetic distances were established for four loci by analysis of 192 (or 222) meiotic events and indicated the following gene order: (centromere)-Es-1-11.7 cM-Tat-8.3 cM-Actsk-1-0.5 cM-Aprt. Mapping of ACTSK to human Chr 1 and of TAT and APRT to human Chr 16 demonstrates the existence of a new short region of homology between mouse Chr 8 and human Chr 1. Intermingling on this scale between human and mouse chromosomal homologies that occurred during evolution creates disorders in comparative linkage studies.     

Sensitivity of lod scores to changes in diagnostic status.

This paper investigates effects on lod scores when one individual in a data set changes diagnostic or recombinant status. First we examine the situation in which a single offspring in a nuclear family changes status. The nuclear-family situation, in addition to being of interest in its own right, also has general theoretical importance, since nuclear families are "transparent"; that is, one can track genetic events more precisely in nuclear families than in complex pedigrees. We demonstrate that in nuclear families log10 [(1-theta)/theta] gives an upper limit on the impact that a single offspring's change in status can have on the lod score at that recombination fraction (theta). These limits hold for a fully penetrant dominant condition and fully informative marker, in either phase-known or phase-unknown matings. Moreover, log10 [(1-theta)/theta] (where theta denotes the value of theta at which Zmax occurs) gives an upper limit on the impact of a single offspring's status change on the maximum lod score (Zmax). In extended pedigrees, in contrast to nuclear families, no comparable limit can be set on the impact of a single individual on the lod score. Complex pedigrees are subject to both stabilizing and destabilizing influences, and these are described. Finally, we describe a "sensitivity analysis," in which, after all linkage analysis is completed, every informative individual in the data set is changed, one at a time, to see the effect which each separate change has on the lod scores. The procedure includes identifying "critical individuals," i.e., those who would have the greatest impact on the lod scores, should their diagnostic status in fact change. To illustrate use of the sensitivity analysis, we apply it to the large bipolar pedigree reported by Egeland et al. and Kelsoe et al. We show that the changes in lod scores observed there, on the order of 1.1-1.2 per person, are not unusual. We recommend that investigators include a sensitivity analysis as a standard part of reporting the results of a linkage analysis.     

Chromosome walking to the AVR1-CO39 avirulence gene of Magnaporthe grisea: discrepancy between the physical and genetic maps.

The avrCO39 gene conferring avirulence toward rice cultivar CO39 was previously mapped to chromosome 1 of Magnaporthe grisea between cosegregating markers CH5-120H and 1.2H and marker 5-10-F. In the present study, this region of the chromosome was physically mapped using RecA-mediated Achilles'' cleavage. Cleavage of genomic DNA sequences within CH5-120H and 5-10-F liberated a 610-kb restriction fragment, representing the physical distance between these markers. Chromosome walking was initiated from both markers but was curtailed due to the presence of repetitive DNA sequences and the absence of overlapping clones in cosmid libraries representing several genome equivalents. These obstacles were overcome by directly subcloning the target region after release by Achilles'' cleavage and a contig spanning avrCO39 was thus assembled. Transformation of two cosmids into a virulent recipient strain conferred a cultivar-specific avirulence phenotype thus confirming the cloning of avrCO39. Meiotic crossover points were unevenly distributed across this chromosomal region and were clustered around the avrCO39 locus. A 14-fold variation in the relationship between genetic and physical distance was measured over the avrCO39 chromosomal region. Thus the poor correlation of physical to genetic distance previously observed in M. grisea appears to be manifested over relatively short distances.     

Chromosome assignments of the human TNF p55 and p75 receptor genes.

At least two different receptor molecules have been described that are capable of binding tumor necrosis factor alpha, a cytokine that plays an important role in inflammation and antitumor activity. Comparative analyses at the nucleotide sequence level suggest that these receptors are members of a newly defined protein family that also includes human and rat nerve growth factor receptors. In this study, we determine the chromosome assignments of the human TNF alpha receptor genes, one of which may have evolved as part of a conserved Hox locus-containing chromosome segment.     

NMR assignments of the binary hvDHFR1:folate complex

A better understanding of how salt affects enzyme activity can be gained via NMR studies of binary hvDHFR1:folate complex. Chemical shift assignments of the 17.9 kDa enzyme with bound substrate prepare the way for ongoing research of the effects of salt on enzyme flexibility through relaxation studies.     

Application of the lod method to the detection of linkage between a quantitative trait and a qualitative marker: a simulation experiment.

Simulation experiments were performed to assess the power of the lod method for detecting linkage between a quantitative trait and a qualitative marker. Using family data on the Amy2-Duffy linkage, it was found that linkage detection is feasible in certain limited circumstances. The same qualitative data yielded no evidence for a significant difference in male and female recombination fractions.     

Reduction of the surface-volume ratio: a physical mechanism contributing to the loss of red cell deformability in malaria

Plasmodia and other intraerythrocytic parasites reduce the deformability of the red cells they infect. One mechanism potentially responsible for this reduction in deformability is the decrease in the surface:volume (S/V) ratio of the red cell which occurs with parasite growth. To examine this hypothesis, normal red cells were allowed to phagocytize polylysine-coated latex spheres 1.0 to 2.9 microns in diameter. Deformability decreased progressively with spheres of increasing size, consistent with the decreasing S/V ratios of those cells (from an initial length:width [L/W] ratio of 2.398 +/- 0.549 for normal red cells to 1.559 +/- 0.249 for red cells containing 2.92 microns latex spheres at 40 dynes per cm2, p less than 0.001). Nevertheless, red cells containing latex spheres 2.0-2.9 microns in diameter remained deformable and continued to tank tread, in contrast to red cells containing Plasmodium falciparum parasites of that size, which are not deformable and do not tank tread. The progressive decrease in S/V produced by the latex spheres is consistent with their effect on the L/W ratio. However, the total loss of deformability observed with red cells containing parasites of similar or smaller size cannot be explained on these grounds alone. It suggests an additional mechanism, such as calcium-induced crosslinking of the red cell cytoskeleton.     

Two-locus maximum lod score analysis of a multifactorial trait: joint consideration of IDDM2 and IDDM4 with IDDM1 in type 1 diabetes.

To investigate the genetic component of multifactorial diseases such as type 1 (insulin-dependent) diabetes mellitus (IDDM), models involving the joint action of several disease loci are important. These models can give increased power to detect an effect and a greater understanding of etiological mechanisms. Here, we present an extension of the maximum lod score method of N. Risch, which allows the simultaneous detection and modeling of two unlinked disease loci. Genetic constraints on the identical-by-descent sharing probabilities, analogous to the "triangle" restrictions in the single-locus method, are derived, and the size and power of the test statistics are investigated. The method is applied to affected-sib-pair data, and the joint effects of IDDM1 (HLA) and IDDM2 (the INS VNTR) and of IDDM1 and IDDM4 (FGF3-linked) are assessed with relation to the development of IDDM. In the presence of genetic heterogeneity, there is seen to be a significant advantage in analyzing more than one locus simultaneously. Analysis of these families indicates that the effects at IDDM1 and IDDM2 are well described by a multiplicative genetic model, while those at IDDM1 and IDDM4 follow a heterogeneity model.     

Linkage of Agt and Actsk-1 to distal mouse Chromosome 8 loci: a new conserved linkage

Angiotensinogen is an ₂ involved in the maintenance of blood pressure and electrolyte balance. We have refined the position of the mouse angiotensinogen locus (Agt) on Chromosome (Chr) 8 and have also confirmed the assignment of the human angiotensinogen locus (AGT) to Chr 1. The segregation of several restriction fragment length variants (RFLVs) was followed in two interspecific backcross sets and in four recombinant inbred (RI) mouse sets. Analysis of the segregation patterns closely linked Agt to Aprt and Emv-2, which places the angiotensionogen locus on the distal end of mouse Chr 8. Additionally, a literature search has revealed that the strain distribution pattern (SDP) for the mouse skeletal -actin locus 1 (Actsk-1, previously Actal, Acta, or Acts) is nearly identical to the SDP for Agt in two RI sets. On the basis of this information we were able to reassign Actsk-1 to mouse Chr 8. By screening a panel of human-mouse somatic cell hybrids, we confirmed that the human angiotensioogen locus lies on Chr 1. This information describes a new region of conserved linkage homology between mouse Chr 8 and human Chr 1. It also defines the end of a large region of conserved linkage homology between mouse Chr 8 and human Chr 16.     

A genetic and physical map of bovine Chromosome 11

A genetic map of bovine Chromosome (Chr) 11 (BTA11, synteny group U16) has been constructed from 330 animals belonging to 21 families, which constitute the international bovine reference panel (IBRP). This map is based on 13 polymorphic microsatellite markers, two of which were chosen in previously published maps. Three markers have been isolated from cosmids. Two of the three cosmids have been physically localized by fluorescence in situ hybridization (FISH), to anchor the genetic map on the chromosome. In addition, a biallelic polymorphism in the -lactoglobulin gene (LGB) has been genetically positioned relative to the microsatellite markers. The most probable order of the markers is: cen-INRA044-BM716-INRA177-(TGLA 327, INRA198, INRA131)-INRA111-INRABERN169-(INRA115, INRA032)-INRA108-INRABERN162-INRA195-LGB. T The total linkage group spans 126 cM, which probably corresponds to most of the chromosome length. The average intermarker distance is about 10.5 cM, allowing the potential detection of a genetic linkage with any Economic Trait Loci (ETL) of this chromosome.