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.     

Statistical Analysis of Uniparental Disomy Data Using Hidden Markov Models

Genetic studies of uniparental disomy (UPD) employing many markers have helped geneticists to gain a better understanding of the molecular mechanisms underlying nondisjunction. However, most existing methods cannot simultaneously analyze all genetic markers and consistently incorporate crossover interference; they thus fail to make the most use of genetic information in the data. In the present article, we describe a hidden Markov model for multilocus uniparental disomy data. This method is based on the chi-square model for the crossover process and can simultaneously incorporate all marker information including untyped and uninformative markers. We then apply this novel method to analyze a set of UPD15 data.     

Sample size requirement for detecting interference under the chi-square model.

The chi-square model (CHS) of recombination has been studied extensively, in recent years, for its ability of capturing the process and estimating the level of crossover interference. Exploration thus far shows that this model yields much better fits to human genetic data than Haldane's no-interference model, and the explicit level of interference can be easily estimated as well. This paper provides calculations of sample sizes required to detect interference under CHS for a variety of settings. Two data types, fully informative meioses and phase-unknown backcross families, are studied. Under each setting, we calculate the number of meioses/families needed to ensure that the expected log-likelihood difference between the chi(2) interference model and Haldane's no-interference model exceeds a prespecified threshold. It is found that joint consideration of multiple (more than three) markers dramatically reduces the number of meioses/families needed when compared to an analysis based on three-point data, a traditional setting for detecting interference using three-point tests. The results indicate that the numbers of meioses needed to detect interference under CHS are well within the reach of most genetic mapping studies.     

A quantitative assay for crossover and noncrossover molecular events at individual recombination hotspots in both male and female gametes

Meiotic recombination is a fundamental process in all eukaryotes. Among organisms in which recombination initiates prior to synapsis, recombination preferentially occurs in short 1-to 2-kb regions, known as recombination hotspots. Among mammals, genotyping sperm DNA has provided a means of monitoring recombination events at specific hotspots in male meiosis. To complement these current techniques, we developed an assay for amplifying all copies of a hotspot from the DNA of male and female germ cells, cloning the products into Escherichia coli, and SNP genotyping the resulting colonies using fluorescence technology. This approach examines the molecular details of crossover and noncrossover events of individual meioses directly at active hotspots while retaining the simplicity of using pooled DNA. Using this technique, we analyzed recombination events at the Hlx1 hotspot located on mouse chromosome 1, finding that the results agree well with a prior genetic characterization of 3026 male and 3002 female meioses.     

Recombination and maternal age-dependent nondisjunction: molecular studies of trisomy 16.

Trisomy 16 is the most common human trisomy, occurring in > or = 1% of all clinically recognized pregnancies. It is thought to be completely dependent on maternal age and thus provides a useful model for studying the association of increasing maternal age and nondisjunction. We have been conducting a study to determine the parent and meiotic stage of origin of trisomy 16 and the possible association of nondisjunction and aberrant recombination. In the present report, we summarize our observations on 62 spontaneous abortions with trisomy 16. All trisomies were maternally derived, and in virtually all the error occurred at meiosis I. In studies of genetic recombination, we observed a highly significant reduction in recombination in the trisomy-generating meioses by comparison with normal female meioses. However, most cases of trisomy 16 had at least one detectable crossover between the nondisjoined chromosomes, indicating that it is reduced--and not absent--recombination that is the important predisposing factor. Additionally, our data indicate an altered distribution of crossing-over in trisomy 16, as we rarely observed crossovers in the proximal long and short arms. Thus, it may be that, at least for trisomy 16, the association between maternal age and trisomy is due to diminished recombination, particularly in the proximal regions of the chromosome.     

Centromeric genotyping and direct analysis of nondisjunction in humans: Down syndrome

In species with chiasmate meioses, alterations in genetic recombination are an important correlate of nondisjunction. In general, these alterations fall into one of two categories: either homologous chromosomes fail to pair and/or recombine at meiosis I, or they are united by chiasmata that are suboptimally positioned. Recent studies of human nondisjunction suggest that these relationships apply to our species as well. However, methodological limitations in human genetic mapping have made it difficult to determine whether the important determinant(s) in human nondisjunction is absent recombination, altered recombination, or both. In the present report, we describe somatic cell hybrid studies of chromosome 21 nondisjunction aimed at overcoming this limitation. By using hybrids to “capture” individual chromosomes 21 of the proband and parent of origin of trisomy, it is possible to identify complementary recombinant meiotic products, and thereby to uncover crossovers that cannot be detected by conventional mapping methods. In the present report, we summarize studies of 23 cases. Our results indicate that recombination in proximal 21q is infrequent in trisomy-generating meioses and that, in a proportion of the meioses, recombination does not occur anywhere on 21q. Thus, our observations indicate that failure to recombine is responsible for a proportion of trisomy 21 cases. Received: 12 January 1997; in revised form: 16 February 1998 / Accepted: 19 February 1998     

Exploiting spore-autonomous fluorescent protein expression to quantify meiotic chromosome behaviors in Saccharomyces cerevisiae

The budding yeast Saccharomyces cerevisiae has proven to be a rich source of information about the mechanisms and regulation of homologous recombination during meiosis. A common technique for studying this process involves microdissecting the four products (ascospores) of a single meiosis and analyzing the configuration of genetic markers in the spores that are viable. Although this type of analysis is powerful, it can be laborious and time-consuming to characterize the large numbers of meioses needed to generate statistically robust data sets. Moreover, the reliance on viable (euploid) spores has the potential to introduce selection bias, especially when analyzing mutants with elevated frequencies of meiotic chromosome missegregation. To overcome these limitations, we developed a versatile, portable set of reporter constructs that drive fluorescent protein expression specifically in only those spores that inherit the reporter. These spore-autonomous fluorescence constructs allow direct visualization of inheritance patterns in intact tetrads, eliminating the need for microdissection and permitting meiotic segregation patterns to be ascertained even in aneuploid spores. As proof of principle, we demonstrate how different arrangements of reporters can be used to quantify crossover frequency, crossover interference, gene conversion, crossover/noncrossover ratios, and chromosome missegregation.     

Human X chromosome markers and Duchenne muscular dystrophy.

Two DNA markers, a random DNA fragment 754 and the cDNA sequence encoding the gene for ornithine transcarbamylase (OTC) have been studied in kindreds segregating for Duchenne muscular dystrophy. 754 and OTC are located close physically to the mutation in the region Xp21 below the breakpoints in two Duchenne females. The genetic distance was found to be approximately 10cM between 754 and DMD (two crossovers in 26 meioses) and to be approximately 10cM between OTC and DMD (two crossovers in 26 meioses). Physical data suggest the order DMD-754-OTC. The frequency of recombination compared to physical distance between these markers and DMD suggests that there may be a hot spot of recombination. The relevance of these observations for the isolation of the DMD mutation and clinical use of these probes is discussed.     

Trans events associated with crossovers are revealed in the absence of mismatch repair genes in Saccharomyces cerevisiae

Genetic analysis of recombination in Saccharomyces cerevisiae has revealed products with structures not predicted by the double-strand break repair model of meiotic recombination. A particular type of recombinant containing trans heteroduplex DNA has been observed at two loci. Trans events were originally identified only in tetrads in which the non-Mendelian segregations were not associated with a crossover. Because of this, these events were proposed to have arisen from the unwinding of double Holliday junctions. Previous studies used palindromes, refractory to mismatch repair, as genetic markers whereas we have used a complementary approach of deleting mismatch repair proteins to identify heteroduplex DNA. We found that the markers occurred in trans and were associated with crossovers. In both mlh1Delta and msh2Delta strains, the frequency of trans events associated with a crossover exceeded that predicted from the random association of crossovers with noncrossover trans events. We propose two different models to account for trans events associated with crossovers and discuss the relevance to wild-type DSB repair.     

Differences in recombination frequencies during female and male meioses of the sex chromosomes of the medaka, Oryzias latipes

In the medaka, Oryzias latipes, sex is determined chromosomally. The sex chromosomes differ from those of mammals in that the X and Y chromosomes are highly homologous. Using backcross panels for linkage analysis, we mapped 21 sequence tagged site (STS) markers on the sex chromosomes (linkage group 1). The genetic map of the sex chromosome was established using male and female meioses. The genetic length of the sex chromosome was shorter in male than in female meioses. The region where male recombination is suppressed is the region close to the sex-determining gene y, while female recombination was suppressed in both the telomeric regions. The restriction in recombination does not occur uniformly on the sex chromosome, as the genetic map distances of the markers are not proportional in male and female recombination. Thus, this observation seems to support the hypothesis that the heterogeneous sex chromosomes were derived from suppression of recombination between autosomal chromosomes. In two of the markers, Yc-2 and Casp6, which were expressed sequence-tagged (EST) sites, polymorphisms of both X and Y chromosomes were detected. The alleles of the X and Y chromosomes were also detected in O. curvinotus, a species related to the medaka. These markers could be used for genotyping the sex chromosomes in the medaka and other species, and could be used in other studies on sex chromosomes.     

A consensus map for loblolly pine (Pinus taeda L.). I. Construction and integration of individual linkage maps from two outbred three-generation pedigrees.

A consensus map for loblolly pine (Pinus taeda L.) was constructed from the integration of linkage data from two unrelated three-generation outbred pedigrees. The progeny segregation data from restriction fragment length polymorphism, random amplified polymorphic DNA, and isozyme genetic markers from each pedigree were recoded to reflect the two independent populations of parental meioses, and genetic maps were constructed to represent each parent. The rate of meiotic recombination was significantly greater for males than females, as was the average estimate of genome length for males (1983.7 cM [Kosambi mapping function (K)]) and females [1339.5 cM(K)]. The integration of individual maps allows for the synthesis of genetic information from independent sources onto a single consensus map and facilitates the consolidation of linkage groups to represent the chromosomes n = 12 of loblolly pine. The resulting consensus map consists of 357 unique molecular markers and covers approximately 1300 cM(K).     

The consequences of uniparental disomy and copy number neutral loss-of-heterozygosity during human development and cancer

UPD (uniparental disomy) describes the inheritance of a pair of chromosomes from only one parent. Mechanisms that lead to UPD include trisomy rescue, gamete complementation, monosomy rescue and somatic recombination. Most of these mechanisms can involve aberrant chromosomes, particularly isochromosomes and Robertsonian translocations. In the last decade, the number of UPD cases reported in the literature has increased exponentially. This is partly due to the advances in genomic technologies that have allowed for high‐resolution SNP (single nucleotide polymorphism) studies, which have complemented traditional methods relying on polymorphic microsatellite markers. In this review, we discuss aberrant cellular mechanisms leading to UPD and their impact on gene expression. Special emphasis is placed on the unmasking of mutant recessive alleles and the disruption of imprinted gene dosage, which give rise to specific and recurrent imprinting phenotypes. Finally, we discuss how copy number maps determined from SNP array datasets have helped identify not only deletions and duplications but also recurrent copy number neutral regions of loss‐of‐heterozygosity, which have been reported in many cancer types and that may constitute an important driving force in cancer. These tiny regions of UPD also alter imprinted gene dosage, which may have cumulative tumourgenic effects in addition to that of unmasking homozygous cancer‐associated mutations.     

Genetic crosses in the apicomplexan parasite Cryptosporidium parvum define recombination parameters

Recombinant progeny lines of Cryptosporidium parvum were generated by coinfecting immunosuppressed mice with two genetically distinct isolates of C. parvum. Progeny lines were obtained from a cross of parental lines MD x TU114 through targeted propagation in mice of progeny oocysts originating from populations lacking one parental allele at one or more loci. For each infection lineage this process was repeated until only a single allele remained for each marker, indicating that the progeny line was clonal. To study genetic recombination, 16 progeny clones were genotyped at 40 loci located on each of the eight chromosomes. The inheritance of parental alleles was significantly skewed towards the more virulent parent isolate MD. A contiguous 476 kb segment of chromosome V displayed MD allele in all progeny recovered, while MD and TU114 alleles were detected at other loci throughout the genome. The absence of alleles from one parental isolate in this chromosomal region may indicate phenotypic selection for the MD allele during the generation of these lines. A range for the meiotic crossover frequency was determined on the basis of 40 markers and the number of meioses estimated to have taken place during the crossing experiment. C. parvum exhibits a high rate of recombination commensurate with other Apicomplexa.     

Uniparental disomy analysis in trios using genome-wide SNP array and whole-genome sequencing data imply segmental uniparental isodisomy in general populations

Whole chromosomal and segmental uniparental disomy (UPD) is one of the causes of imprinting disorder and other recessive disorders. Most investigations of UPD were performed only using cases with relevant phenotypic features and included few markers. However, the diagnosis of cases with segmental UPD requires a large number of molecular investigations. Currently, the accurate frequency of whole chromosomal and segmental UPD in a normal developing embryo is not well understood. Here, we present whole chromosome and segmental UPD analysis using single nucleotide polymorphism (SNP) microarray data of 173 mother-father-child trios (519 individuals) from six populations (including 170 HapMap trios). For two of these trios, we also investigated the possibility of shorter segmental UPD as a consequence of homologous recombination repair (HR) for DNA double strand breaks (DSBs) during the early developing stage using high-coverage whole-genome sequencing (WGS) data from 1000 Genomes Project. This could be overlooked by SNP microarray. We identified one obvious segmental paternal uniparental isodisomy (iUPD) (8.2 mega bases) in one HapMap sample from 173 trios using Genome-Wide Human SNP Array 6.0 (SNP6.0 array) data. However, we could not identify shorter segmental iUPD in two trios using WGS data. Finally, we estimated the rate of segmental UPD to be one per 173 births (0.578%) based on the UPD screening for 173 trios in general populations. Based on the autosomal chromosome pairs investigated, we estimate the rate of segmental UPD to be one per 3806 chromosome pairs (0.026%). These data imply the possibility of hidden segmental UPD in normal individuals.     

An RFLP linkage map of Lycopersicon peruvianum

In order to map genes determining resistance to bacterial canker in tomato, backcrosses were made between a resistant and a susceptible Lycopersicon peruvianum accession. The linkage study with RFLP markers yielded a genetic map of L. Peruvianum. This map was compared to that derived from a L. esculentum x L. pennellii F₂ population, based on 70 shared RFLP markers. The maps showed a good resemblance in both the order of markers and the length of the chromosomes, with the exception of just one relocated marker on chromosome 9. Because backcrosses were made with the F₁, either as the pollen parent or as the pistil parent, linkage maps from male and female meioses could be estimated. It was concluded that recombination at male meiosis was reduced, and that gametophytic selection for parental genotypes at more than one locus per chromosome might be partly responsible for the reduction of the estimated male map length.     

On Genetic Map Functions

Various genetic map functions have been proposed to infer the unobservable genetic distance between two loci from the observable recombination fraction between them. Some map functions were found to fit data better than others. When there are more than three markers, multilocus recombination probabilities cannot be uniquely determined by the defining property of map functions, and different methods have been proposed to permit the use of map functions to analyze multilocus data. If for a given map function, there is a probability model for recombination that can give rise to it, then joint recombination probabilities can be deduced from this model. This provides another way to use map functions in multilocus analysis. In this paper we show that stationary renewal processes give rise to most of the map functions in the literature. Furthermore, we show that the interevent distributions of these renewal processes can all be approximated quite well by gamma distributions.     

Multilocus genomics of outcrossing plant populations

The structure and organization of natural plant populations can be understood by estimating the genetic parameters related to mating behavior, recombination frequency, and gene associations with DNA-based markers typed throughout the genome. We developed a statistical and computational model for estimating and testing these parameters from multilocus data collected in a natural population. This model, constructed by a maximum likelihood approach and implemented within the EM algorithm, is shown to be robust for simultaneously estimating the outcrossing rate, recombination frequencies and linkage disequilibria. The algorithm built with three or more markers allows the characterization of crossover interference in meiosis and high-order disequilibria among different genes, thus providing a powerful tool for illustrating a detailed picture of genetic diversity and organization in natural populations. Computer simulations demonstrate the statistical properties of the proposed model. This multilocus model will be useful for studying the pattern and amount of genetic variation within and among populations to further infer the evolutionary history of a plant species.     

Refined physical and genetic mapping of the NF1 region on chromosome 17.

A total of 15 polymorphic markers were used to construct a genetic map that encompasses the NF1 locus on chromosome 17. The markers were a subset of a large collection of chromosome 17-specific probes and were selected for marker typing in NF1 families after physical localization to the pericentric region of the chromosome. Multilocus data for a total of 17 informative NF1 families and 39 other families were included in genetic analyses. No recombination was observed between NF1 and four markers, one or more of which was informative in 86% of parents. More-refined physical mapping studies demonstrated that all four of the markers are proximal to the chromosome 17 translocation breakpoints from two NF1 patients bearing balanced translocations. The region flanking the disease locus spans a distance of 1 centimorgan (cM) in males and 9 cM in females. Close flanking markers were informative in 76% of meioses. Sex differences in recombination rates in the pericentric region were highly significant statistically.     

The Theory and Applications of Measuring Broad-Range and Chromosome-Wide Recombination Rate from Allele Frequency Decay around a Selected Locus

Recombination is the exchange of genetic material between homologous chromosomes via physical crossovers. High-throughput sequencing approaches detect crossovers genome wide to produce recombination rate maps but are difficult to scale as they require large numbers of recombinants individually sequenced. We present a simple and scalable pooled-sequencing approach to experimentally infer near chromosome-wide recombination rates by taking advantage of non-Mendelian allele frequency generated from a fitness differential at a locus under selection. As more crossovers decouple the selected locus from distal loci, the distorted allele frequency attenuates distally toward Mendelian and can be used to estimate the genetic distance. Here, we use marker selection to generate distorted allele frequency and theoretically derive the mathematical relationships between allele frequency attenuation, genetic distance, and recombination rate in marker-selected pools. We implemented nonlinear curve-fitting methods that robustly estimate the allele frequency decay from batch sequencing of pooled individuals and derive chromosome-wide genetic distance and recombination rates. Empirically, we show that marker-selected pools closely recapitulate genetic distances inferred from scoring recombinants. Using this method, we generated novel recombination rate maps of three wild-derived strains of Drosophila melanogaster, which strongly correlate with previous measurements. Moreover, we show that this approach can be extended to estimate chromosome-wide crossover interference with reciprocal marker selection and discuss how it can be applied in the absence of visible markers. Altogether, we find that our method is a simple and cost-effective approach to generate chromosome-wide recombination rate maps requiring only one or two libraries.     

A genetic linkage map of human chromosome 5 with 60 RFLP loci.

A genetic map of human chromosome 5 that contains 60 restriction fragment length polymorphism (RFLP) loci in one linkage group has been constructed. Segregation data using these markers and 40 large multigenerational families supplied by the Centre d'Etude du Polymorphisme Humain have been collected. Linkage analyses were performed with the program package CRI-MAP; using odds greater than 1000:1, 30 RFLP loci could be placed on the map. This genetic map spans 289 cM sex-equal, 353 cM in females, and 244 cM in males. While the relative rate of recombination for female meioses is nearly twice that of males over much of the chromosome, several instances of statistically significant excess male recombination were observed. The order of probes on the genetic map has been confirmed by their physical order as determined by somatic cell hybrid lines containing deletions of normal chromosome 5. There is concordance between the physical positions of markers and their genetic positions. Our most distal probes on the genetic map are cytologically localized to the most distal portions of the chromosome. This suggests that our genetic map spans most of chromosome 5.