Multipoint Genetic Mapping with Trisomy Data

Trisomy is the most common genetic abnormality in humans and is the leading cause of mental retardation. Although molecular studies that use a large number of highly polymorphic markers have been undertaken to understand the recombination patterns for chromosome abnormalities, there is a lack of multilocus approaches to incorporating crossover interference in the analysis of human trisomy data. In the present article, we develop two statistical methods that simultaneously use all genetic information in trisomy data. The first approach relies on a general relationship between multilocus trisomy probabilities and multilocus ordered-tetrad probabilities. Under the assumption that no more than one chiasma exists in each marker interval, we describe how to use the expectation-maximization algorithm to examine the probability distribution of the recombination events underlying meioses that lead to trisomy. One limitation of the first approach is that the amount of computation increases exponentially with the number of markers. The second approach models the crossover process as a chi(2) model. We describe how to use hidden Markov models to evaluate multilocus trisomy probabilities. Our methods are applicable when both parents are available or when only the nondisjoining parent is available. For both methods, genetic distances among a set of markers can be estimated and the pattern of overall chiasma distribution can be inspected for differences in recombination between meioses exhibiting trisomy and normal meioses. We illustrate the proposed approaches through their application to a set of trisomy 21 data.     

A primary genetic map of chromosome 13q.

We have constructed a primary genetic map spanning most of human chromosome 13. A total of 14 polymorphic DNA sequences and one protein polymorphism provided, after construction of haplotypes, seven markers for the long arm of this chromosome. A panel of cell lines from 30 three-generation families with large sibship size served as the sample set. Pairwise cross analysis of the inheritance patterns of the marker loci established that six of the seven loci constituted a single linkage group; the seventh was localized by physical means. Significantly higher recombination rates were found in female than in male meioses in several intervals. The six closely linked loci were arranged, based on the two-point data, in three clusters, and a number of alternate gene orders were excluded by three-point linkage tests. The order and spacing of the individual loci were refined by linkage analyses that considered five loci jointly.     

Pedigree data analysis with crossover interference

We propose a new method for calculating probabilities for pedigree genetic data that incorporates crossover interference using the chi-square models. Applications include relationship inference, genetic map construction, and linkage analysis. The method is based on importance sampling of unobserved inheritance patterns conditional on the observed genotype data and takes advantage of fast algorithms for no-interference models while using reweighting to allow for interference. We show that the method is effective for arbitrarily many markers with small pedigrees.     

An assay for chromosomal and chromatid interference in chromosome V ofSaccharomyces

This survey included 1250 tetrads from 19 families heterozygous for at least three genes on chromosome V. A number of tests were used for the detection of interference. The Papazian tests (table 1) indicated that either positive chromosomal interference or negative chromatid interference occurs in some regions; this effect is too weak to be manifest in individual families. Although chromosomal interference was not evident in the three-point tests of any one family, there was some evidence from lumped data that it occurs. Negative chromatid interference is clearly evident in the three-point tetrad data. This effect varies from family to family, and does not appear to be localized to certain regions. The data indicate that chromatid interference is a property conferred by the specific hybrid rather than a fixed property of any given chromosomal region. In this connection it is interesting to note that family 118, which exhibits chromatid interference, was obtained from a intraascal hybrid from family 108, which also shows chromatid interference, suggesting that this family-specific effect may be hereditary.This work has been supported by the Americal Cancer Society.     

Patterns of postzygotic isolation in Lepidoptera

I present patterns characterizing the evolution of intrinsic postzygotic isolation in Lepidoptera by analyzing data from the literature on genetic distance, strength of hybrid sterility and inviability, biogeography, and natural hybridization. Using genetic distance as a proxy for time, I investigate the time-course of the evolution of postzygotic isolation and the waiting times to particular hybrid fitness problems. The results show that postzygotic isolation increases gradually as species diverge, but that hybrid sterility evolves faster than hybrid inviability. The overwhelming preponderance of female-specific hybrid problems in Lepidoptera shows that Haldane's rule (the preferential sterility or inviability of the heterogametic sex) is well obeyed. Together the rates and patterns characterizing the accumulation of postzygotic isolation allow several tests of the composite theory of Haldane's rule. Interestingly, comparing these data with those from Drosophila reveals that Haldane's rule for sterility evolves as fast (if not faster) in Lepidoptera. Finally, I show that a substantial fraction of sympatric species hybridizes in nature and that the majority of these suffer some level of hybrid sterility or inviability.     

Multipoint genetic mapping with uniparental disomy data

Uniparental disomy (UPD) refers to the presence of two copies of a chromosome from one parent and none from the other parent. In genetic studies of UPDs, many genetic markers are usually used to identify the stage of nondisjunction that leads to UPD and to uncover the associated unusual patterns of recombinations. However, genetic information in such data has not been fully utilized because of the limitations of the existing statistical methods for UPD data. In the present article, we develop a multilocus statistical approach that has the advantages of being able to simultaneously consider all genetic markers for all individuals in the same analysis and to allow general models for the crossover process to incorporate crossover interference. In particular, for a general crossover-process model that assumes only that there exists in each interval at most one crossover, we describe how to use the expectation-maximization algorithm to examine the probability distribution of the recombination events underlying meioses leading to UPD. We can also use this flexible approach to create genetic maps based on UPD data and to inspect recombination differences between meioses exhibiting UPD and normal meioses. The proposed method has been implemented in a computer program, and we illustrate the proposed approach through its application to a set of UPD15 data.     

A mapped set of DNA markers for human chromosome 15

A primary genetic linkage map for human chromosome 15 has been constructed from 16 arbitrary DNA markers genotyped in 59 large reference families. The map spans a genetic distance of 146 cM in males and 187 cM in females. The ratio of female/male genetic distance was approximately 2.1 overall within the region of the chromosome covered by our map, but three segments showed a significant male excess in recombination frequency. A subset of seven of the linked markers would be enough to detect linkage of a genetic defect within the mapped region of chromosome 15, if at least 48 phase-known meioses in affected families were available for analysis.     

A mapped set of genetic markers for human chromosome 9

A genetic map of markers for human chromosome 9, spanning a genetic distance of 147 cM in males and 231 cM in females, has been constructed from linkage studies with 19 loci in a large panel of reference families. The markers included four classical systems previously assigned to chromosome 9, and restriction fragment length polymorphisms of two cloned genes, ABL oncogene and argininosuccinase synthetase pseudogene 3 (ASSP3). The remaining 13 marker loci, with an average heterozygosity of 42%, were defined by arbitrary DNA probes newly ascertained from genomic libraries; seven of them were variable number of tandem repeat (VNTR) loci. A subset of 7 of the 19 linked markers is proposed for a primary map that could detect linkage with a genetic defect within the covered region of chromosome 9, provided that at least 45 phase-known meioses were available for study in an affected family.     

Defining the Proximal Border of the Huntington Disease Candidate Region by Multipoint Recombination Analyses

The candidate region for the Huntington disease (HD) gene has been narrowed down to a 2.2-Mb region between D4S10 and D4S98 on the short arm of chromosome 4. To map the HD gene within this candidate region 65 Dutch HD families were studied. In total 338 informative meioses were analyzed and 11 multiple informative crossovers were detected. Assuming a minimum number of recombinations and no double recombinations, our multiple informative crossovers are consistent with one specific genetic order for 12 loci: D4S10-(D4S81, D4S126)-D4S125-(D4S127, D4S95)-D4S43-(D4S115, D4S96, D4S111, D4S90, D4S141). This is in agreement with the known data derived from similar and other methods. The loci between brackets could not be mapped relative to each other. In our family material, two informative three-point marker recombination events were detected in the proximal HD candidate region, which are also informative for HD. Both recombination events map the HD gene distal to D4S81 and most likely distal to D4S125, narrowing down the HD candidate region to a 1.7-Mb region between D4S125 and D4S98.     

Mapping of recombinants near the Huntington disease locus by using G8 (D4S10) and newly isolated markers in the D4S10 region.

Genetic linkage between the marker G8 (D4S10) and Huntington disease (HD) was studied in six Dutch pedigrees. The informativeness of the D4S10 locus was increased by isolation of a cosmid, C5.5, with a G8 subclone used as probe. We present a restriction map of 70 kb in the D4S10 region. Two subclones of C5.5, H5.52 and F5.53, detect MspI and SinI RFLPs, respectively. These probes increase the informativeness of D4S10 in the Dutch HD population from 55% to 95%. Seven recombinations were found in 124 informative meioses in which multipoint segregation of D4S10 haplotypes and the HD locus was studied. Two of the recombinations occurred within the D4S10 region. The other five recombinations are highly valuable for the mapping of present and future markers relative to each other and to the HD gene. In addition, several recombinations between markers in meioses from unaffected parents were noted, which will also be useful in ordering new markers. On the basis of our three-point recombination data, the orientation of the D4S10 region relative to HD is HD-H5.52-G8-F5.53, which independently confirms the previously derived polarity for D4S10.     

Testing for linkage under robust genetic models.

Robust genetic models are used to assess linkage between a quantitative trait and genetic variation at a specific locus using allele-sharing data. Little is known about the relative performance of different possible significance tests under these models. Under the robust variance components model approach there are several alternatives: standard Wald and likelihood ratio tests, a quasilikelihood Wald test, and a Monte Carlo test. This paper reports on the relative performance (significance level and power) of the robust sibling pair test and the different alternatives under the robust variance components model. Simulations show that (1) for a fixed sample size of nuclear families, the variance components model approach is more powerful than the robust sibling pair approach; (2) when the number of nuclear families is at least approximately 100 and heritability at the trait locus is moderate to high (>0.20) all tests based on the variance components model are equally effective; (3) when the number of nuclear families is less than approximately 100 or heritability at the trait locus is low (<0. 20), on balance, the Monte Carlo test provides the best power and is the most valid. The different testing procedures are applied to determine which are able to detect the known association between low density lipoprotein cholesterol and the common genotypes at the locus encoding apolipoprotein E. Results from this application show that the robust sibling pair method may be more effective in practice than that indicated by simulations.     

Epigenetics and phenotypic variation in mammals

What causes phenotypic variation? By now it is clear that phenotype is a result of the interaction between genotype and environment, in addition to variation not readily attributable to either. Epigenetic phenomena associated with phenotypic variation at the biochemical, cellular, tissue, and organism level are now well recognized and are likely to contribute to the “intangible variation” alluded to. While it is clear that epigenetic modifications are mitotically heritable, the fidelity of this process is not well understood. Inheritance through more than one generation of meioses is even less well studied. So it remains unclear to what extent epigenetic changes contribute to phenotypic variation in natural populations. How might such evidence be obtained? What are the features of phenotypes that might suggest an epigenetic component? How much of the epigenetic component is truly independent of genetic changes? The answers to such questions must come from studies designed specifically to detect subtle, stochastically determined phenotypic variation in suitable animal models.     

A Maximum Likelihood Method for Estimating Genome Length Using Genetic Linkage Data

The genetic length of a genome, in units of Morgans or centimorgans, is a fundamental characteristic of an organism. We propose a maximum likelihood method for estimating this quantity from counts of recombinants and nonrecombinants between marker locus pairs studied from a backcross linkage experiment, assuming no interference and equal chromosome lengths. This method allows the calculation of the standard deviation of the estimate and a confidence interval containing the estimate. Computer simulations have been performed to evaluate and compare the accuracy of the maximum likelihood method and a previously suggested method-of-moments estimator. Specifically, we have investigated the effects of the number of meioses, the number of marker loci, and variation in the genetic lengths of individual chromosomes on the estimate. The effect of missing data, obtained when the results of two separate linkage studies with a fraction of marker loci in common are pooled, is also investigated. The maximum likelihood estimator, in contrast to the method-of-moments estimator, is relatively insensitive to violation of the assumptions made during analysis and is the method of choice. The various methods are compared by application to partial linkage data from Xiphophorus.     

Evidence for Mendelian inheritance of serum IgE levels in Hispanic and non-Hispanic white families.

Considerable evidence is available suggesting a significant genetic component in the pathogenesis of asthma, but the mechanism of inheritance is not well understood. The main objective of this study was to assess if total serum IgE level, a known intermediate phenotype for asthma, is under the control of a major autosomal gene. We studied nuclear families participating in the Tucson Children's Respiratory Study in Tucson and originally selected because they belonged to a health maintenance organization. One hundred twenty-five Hispanic and 673 non-Hispanic White nuclear families were eligible; 50 Hispanic families (with 191 subjects) and 241 non-Hispanic White families (with 886 subjects) were included. Prevalence of asthma, hay fever, and parental smoking was similar among eligible families who were included and those who were not. Segregation analyses using regressive models for continuous traits showed that the best fit to the data was given by a model of Mendelian codominant inheritance of a major autosomal gene associated with higher serum IgE level. Log-likelihood for this model was not significantly different from that of the best-fitting ("unrestricted") model (P = .3) and was significantly better than log-likelihood for a dominant model (P < .0001) and a recessive model (P < .0001). An environmental model showed significant departure (P < .0001) from the unrestricted model. Tests for genetic heterogeneity showed no significant difference between the two ethnic groups. The data strongly suggest that total serum IgE levels are controlled by a major autosomal codominant gene.     

Crossover interference in the mouse

We present an analysis of crossover interference in the mouse genome, on the basis of high-density genotype data from two reciprocal interspecific backcrosses, comprising 188 meioses. Overwhelming evidence was found for strong positive crossover interference with average strength greater than that implied by the Carter-Falconer map function. There was some evidence for interchromosomal variation in the level of interference, with smaller chromosomes exhibiting stronger interference. We further compared the observed numbers of crossovers to previous cytological observations on the numbers of chiasmata and evaluated evidence for the obligate chiasma hypothesis.     

The structure of genetic linkage data: from LIPED to 1M SNPs

There are three assumptions of independence or conditional independence that underlie linkage likelihood computations on sets of related individuals. The first is the independence of meioses, which gives rise to the conditional independence of haplotypes carried by offspring, given those of their parents. The second derives from the assumption of absence of genetic interference, which gives rise to the conditional independence of inheritance vectors, given the inheritance vector at an intermediate location. The third is the assumption of independence of allelic types, at the population level, both among haplotypes of unrelated individuals and also over the loci along a given haplotype. These three assumptions have been integral to likelihood computations since the first lod scores were computed, and remain key components in analysis of modern genetic data. In this paper we trace the role of these assumptions through the history of linkage likelihood computation, through to a new framework of genetic linkage analysis in the era of dense genomic marker data.     

Genetic dissection of hybrid incompatibilities between Drosophila simulans and D. mauritiana. III. Heterogeneous accumulation of hybrid incompatibilities, degree of dominance, and implications for Haldane's rule

The genetic basis of Haldane's rule was investigated through estimating the accumulation of hybrid incompatibilities between Drosophila simulans and D. mauritiana by means of introgression. The accumulation of hybrid male sterility (HMS) is at least 10 times greater than that of hybrid female sterility (HFS) or hybrid lethality (HL). The degree of dominance for HMS and HL in a pure D. simulans background is estimated as 0.23-0.29 and 0.33-0.39, respectively; that for HL in an F1 background is unlikely to be very small. Evidence obtained here was used to test the Turelli-Orr model of Haldane's rule. Composite causes, especially, faster-male evolution and recessive hybrid incompatibilities, underlie Haldane's rule in heterogametic male taxa such as Drosophila (XY male and XX female). However, if faster-male evolution is driven by sexual selection, it contradicts Haldane's rule for sterility in heterogametic-female taxa such as Lepidoptera (ZW female and ZZ male). The hypothesis of a faster-heterogametic-sex evolution seems to fit the current data best. This hypothesis states that gametogenesis in the heterogametic sex, instead of in males per se, evolves much faster than in the homogametic sex, in part because of sex-ratio selection. This hypothesis not only explains Haldane's rule in a simple way, but also suggests that genomic conflicts play a major role in evolution and speciation.     

Genetic and physical mapping of the Chediak-Higashi syndrome on chromosome 1q42-43.

The Chediak-Higashi syndrome (CHS) is a severe autosomal recessive condition, features of which are partial oculocutaneous albinism, increased susceptibility to infections, deficient natural killer cell activity, and the presence of large intracytoplasmic granulations in various cell types. Similar genetic disorders have been described in other species, including the beige mouse. On the basis of the hypothesis that the murine chromosome 13 region containing the beige locus was homologous to human chromosome 1, we have mapped the CHS locus to a 5-cM interval in chromosome segment 1q42.1-q42.2. The highest LOD score was obtained with the marker D1S235 (Zmax = 5.38; theta = 0). Haplo-type analysis enabled us to establish D1S2680 and D1S163, respectively, as the telomeric and the centromeric flanking markers. Multipoint linkage analysis confirms the localization of the CHS locus in this interval. Three YAC clones were found to cover the entire region in a conting established by YAC end-sequence characterization and sequence-tagged site mapping. The YAC contig contains all genetic markers that are nonrecombinant for the disease in the nine CHS families studied. This mapping confirms the previous hypothesis that the same gene defect causes CHS in human and beige pheno-type in mice and provides a genetic framework for the identification of candidate genes.     

A genetic linkage map of 32 loci on human chromosome 10

We have constructed a genetic linkage map of human chromosome 10 based on DNA probes that detect 47 restriction fragment length polymorphisms (RFLPs) at 32 different loci. Segregation data were collected on a set of multigenerational families provided by the Centre d'Etude du Polymorphisme Humain and maps were constructed using recently developed multipoint analysis techniques. The length of the sex-averaged map is 178 cM and the sex-specific map lengths are 131 cM in males and 255 cM in females. Recombination is significantly higher in female meioses. The mean distance between loci is 5.6 cM for the sex-averaged map. The genetic map spans the length of the chromosome as judged by physical localization of probes by in situ hybridization techniques and mapping of the probes on human-hamster hybrid cell lines containing all or part of chromosome 10. The informativeness of two loci near the locus responsible for multiple endocrine neoplasia type 2A (MEN-2A) has been increased by isolation of cosmids that reveal additional RFLPs at these loci.     

The misuse of the three-point cross for estimating chiasma interference

The building of genetic maps in diploid organisms by crosses between different genotypes and estimation of recombination frequencies from the obtained segregation data has been successfully used since a very young step in the birth of genetics. The three-point cross methodology has facilitated this task and has demonstrated at the same time that genetic distances are not additive, as some recombinant products are not recognised in the progeny. Three-point cross also allows to examine if chiasma interference exists and its evaluation. Here I show that the classical method of this estimation is erroneous and inevitably determines the apparition of a spurious, positive interference, which has been claimed to be an almost general phenomenon. Interference can only be estimated from a precise knowledge of the number of crossing over events occurring in meiocytes.