High resolution analysis of haplotype diversity and meiotic crossover in the human TAP2 recombination hotspot

Little is known about the nature of recombination hotspots in the human genome and the relationship between crossover activity and patterns of linkage disequilibrium. We have therefore used both haplotype analysis and direct detection of crossovers in sperm to characterize a putative recombination hotspot in the TAP2 gene within the class II region of the MHC. Haplotype diversity provided evidence for a localized hotspot within intron 2 of this gene. Sperm DNA typing using allele-specific PCR primers to selectively amplify recombinant TAP2 molecules revealed a highly localized meiotic crossover hotspot approximately 1.2 kb long, unusually abundant in sequence polymorphisms and flanked by DNA much less active in recombination. Sperm crossover appeared to be fully reciprocal, and almost all crossover products were simple, involving a single exchange between adjacent heterozygous markers. This hotspot appears to be much more active in female than male meiosis. No primary sequence similarities could be found between any of the very few well defined crossover hotspots in the human genome, all of which show recombinationally active domains 1-2 kb long. Direct comparison of recombination frequency and haplotype diversity in TAP2 showed that linkage disequilibrium measures were a poor predictor of crossover frequency in this region, with non-recombining markers sometimes in free association and with examples of pairs of markers spanning the recombination hotspot showing substantial or even absolute linkage disequilibrium.     

Recombination hotspots and block structure of linkage disequilibrium in the human genome exemplified by detailed analysis of PGM1 on 1p31

The distribution of linkage disequilibrium (LD) in the human genome has important consequences for the design of experiments that infer susceptibility genes for complex disease using association studies. Recent studies have shown a non-random distribution of human meiotic recombination associated with intervening tracts of LD. Little is known about the processes, patterns and frequency of reciprocal meiotic recombination in humans. However, this phenomenon can be better understood by the fine structure analysis of several genomic regions by mapping hotspots and characterizing regions with variable LD. Here, we report clustered hotspot activity with intervening blocks of LD within the human PGM1 gene (1p31) using data derived from meiotic and population studies. Earlier work has suggested a high recombination rate in two regions within the PGM1 gene, site A (exons 4-8) and site B (exons 1A-4). Sequencing of eight individuals across 6 kb of targeted regions in site B identified 18 informative SNPs. Individuals from three distinct populations, Caucasian (n=264), Chinese (n=222) and Vietnamese (n=187), were genotyped, and haplotypes were determined using estimate of haplotypes, ldmax and Arlequin. Allelic association and haplotype analysis in these samples revealed variable recombination rates across PGM1, demonstrating the presence of: (i) three hotspots and (ii) three haplotype blocks. The spatial arrangement of haplotype blocks was identical in all populations studied. The pattern of association within PGM1 represents a region decomposed into small blocks of LD, where increased recombination activity has disrupted the ancestral chromosome. Additionally, crossovers in phased data mapped preferentially to regions where LD collapses, which also overlap with sequence motifs.     

Recombination hotspots rather than population history dominate linkage disequilibrium in the MHC class II region

Recombination, demographic history, drift and selection influence the extent of linkage disequilibrium (LD) in the human genome, but their relative contributions remain unclear. To investigate the effect of meiotic recombination versus population history on LD, three populations with different demographic histories (UK north Europeans, Saami and Zimbabweans) were genotyped for high-frequency single-nucleotide polymorphisms (SNPs) across a 75 kb DNA segment of the MHC class II region. This region spans three well-characterized recombination hotspots and a 60 kb long LD block. Despite a high level of underlying haplotype diversity and considerable divergence in haplotype composition between populations, all three populations showed very similar patterns of LD. Surprisingly, the entire 60 kb LD block was present even in Africans, although it was relatively difficult to detect owing to a systematic deficiency of high frequency SNPs. In contrast, DNA within recombination hotspots did not show this low nucleotide diversity in Africans. Thus, while population history has some influence on LD, our findings suggest that recombination hotspots play a major global role in shaping LD patterns as well as helping to maintain localized SNP diversity in this region of the MHC.     

Human haplotype block sizes are negatively correlated with recombination rates

The International Haplotype Map ("HapMap") Project is motivated, in part, by the belief that the organization of the human genome, the mechanics of recombination, and the population-level behavior of alleles at adjacent loci should allow researchers to parse the genome into small segments, or "blocks," that show strong linkage disequilibrium (LD) between alleles at loci within those segments. The discovery and evidence for these blocks is to be based solely on the observed LD strength and patterns between alleles at adjacent loci throughout the genome. Although there are many factors that contribute to LD strength, we assessed the correlation between block structure, in terms of length and percentage of the genome assembled into blocks within a region, and recombination rate obtained from two independent sources. We found evidence of a striking negative correlation between the average recombination rate and average block length, suggesting that recombination rate is a strong contributor to haplotype block structure within the genome. We discuss the potential implications of this negative correlation in the context of the organization, properties, and potential ubiquity of a block-like structure in the human genome.     

Unequal exchange at the Charcot-Marie-Tooth disease type 1A recombination hot-spot is not elevated above the genome average rate

An increasing number of human diseases and syndromes are being found to result from micro-duplications or microdeletions arising from meiotic recombination between homologous repeats on the same chromosome. The first microduplication syndrome delineated, Charcot-Marie-Tooth disease type 1A (CMT1A), results from unequal crossing over between two >98% identical 24 kb repeats (CMT1A-REPs) on chromosome 17. In addition to its medical significance, the CMT1A region has features that make it a unique resource for detailed analysis of human unequal recombination. Previous studies of CMT1A patients showed that the majority of unequal crossovers occurred within a small region (<1 kb) of the REPs suggesting the presence of a recombination hot-spot. We directly measured the frequency of unequal recombination in the hot-spot region using sperm from four normal individuals. Surprisingly, unequal recombination between the REPs occurs at a rate no greater than the average rate for the male genome (approximately 1 cM/Mb) and is the same as that expected for equally aligned REPs. This conclusion extends to humans the findings in yeast that recombination between repeated sequences far apart on the same chromosome may occur at similar frequencies to allelic recombination. Finally, the CMT1A hot-spot stands in sharp contrast to the human MS32 mini-satellite-associated hot-spot that exhibits highly enhanced recombination initiation in addition to positional specificity. One possibility is that the CMT1A hot-spot may consist of a region with genome average recombination potential embedded within a recombination cold-spot.     

Haplotype structure,LD blocks,and uneven recombination within the LRP5 gene

Patterns of linkage disequilibrium (LD) in the human genome are beginning to be characterized, with a paucity of haplotype diversity in "LD blocks," interspersed by apparent "hot spots" of recombination. Previously, we cloned and physically characterized the low-density lipoprotein-receptor-related protein 5 (LRP5) gene. Here, we have extensively analysed both LRP5 and its flanking three genes, spanning 269 kb, for single nucleotide polymorphisms (SNPs), and we present a comprehensive SNP map comprising 95 polymorphisms. Analysis revealed high levels of recombination across LRP5, including a hot-spot region from intron 1 to intron 7 of LRP5, where there are 109 recombinants/Mb (4882 meioses), in contrast to flanking regions of 14.6 recombinants/Mb. This region of high recombination could be delineated into three to four hot spots, one within a 601-bp interval. For LRP5, three haplotype blocks were identified, flanked by the hot spots. Each LD block comprised over 80% common haplotypes, concurring with a previous study of 14 genes that showed that common haplotypes account for at least 80% of all haplotypes. The identification of hot spots in between these LD blocks provides additional evidence that LD blocks are separated by areas of higher recombination.     

Multiple recombination sites at the 5'-end of murine coronavirus RNA

Mouse hepatitis virus (MHV), a murine coronavirus, contains a nonsegmented RNA genome. We have previously shown that MHV could undergo RNA-RNA recombination in crosses between temperature-sensitive mutants and wild-type viruses at a very high frequency (S. Makino, J.G. Keck, S.A. Stohlman, and M.M.C. Lai (1986) J. Virol. 57, 729-737). To better define the mechanism of RNA recombination, we have performed additional crosses involving different sets of MHV strains. Three or possibly four classes of recombinants were isolated. Recombinants in the first class, which are similar to the ones previously reported, contain a single crossover in either gene A or B, which are the 5'-most genes. The second class of recombinants contain double crossovers in gene A. The third class of recombinants have crossovers within the leader sequence located at the 5'-end of the genome. The crossover sites of the third class have been located between 35 and 60 nucleotides from the 5'-end of the leader RNA. One of these recombinants has double crossovers within the short region comprising the leader sequences. Finally, we describe one recombinant which may contain a triple crossover. The presence of so many recombination sites within the 5'-end of the genome of murine coronaviruses confirms that RNA recombination is a frequent event during MHV replication and is consistent with our proposed model of "copy-choice" recombination in which RNA replication occurs in a discontinuous and nonprocessive manner.     

Haplotype-specific linkage disequilibrium patterns define the genetic topography of the human MHC

Detailed knowledge of linkage disequilibrium (LD) is regarded as a prerequisite for population-based disease gene mapping. Variable patterns across the human genome are now recognized, both between regions and populations. Here, we demonstrate that LD may also vary within a genomic region in a haplotype-specific manner. In 864 Caucasian unrelated individuals, we describe haplotype-specific LD patterns across the human MHC by the construction of gene-specific allelic haplotypes at 25 loci between HLA-A and Tapasin. Strong and extensive LD is found across both common and rare haplotypes, suggesting that haplotype structure is influenced by factors other than genetic drift, including both selection and differential haplotype recombination. Knowledge of haplotype-specific LD in the HLA may explain the apparent discrepant data from previous studies of global LD, help delineate key areas in mapping HLA-associated diseases and, together with recombination data, provide valuable information about a population's demographic history and the selective pressures operating on it.     

Factors influencing recombination frequency and distribution in a human meiotic crossover hotspot

Little is known about the factors that influence the frequency and distribution of meiotic recombination events within human crossover hotspots. We now describe the detailed analysis of sperm recombination in the NID1 hotspot. Like the neighbouring MS32 hotspot, the NID1 hotspot is associated with a minisatellite, suggesting that hotspots predispose DNA to tandem repetition. Unlike MS32, crossover resolution breakpoints in NID1 avoid the minisatellite, producing a cold spot within the hotspot. This avoidance may be related to the palindromic nature of the minisatellite interfering with the generation and/or processing of recombination intermediates. The NID1 hotspot also contains a single nucleotide polymorphism (SNP) close to the centre, which appears to directly influence the frequency of crossover initiation. Quantitative gene conversion assays show that this SNP affects the frequency of gene conversion and crossover to a very similar extent, providing evidence that conversions and crossovers are triggered by the same recombination initiating events. The recombination-suppressing allele is over-transmitted to recombinant progeny, and provides the most dramatic example to date of recombination-mediated meiotic drive, of a magnitude sufficient to virtually guarantee that the recombination suppressor will eventually replace the more active allele in human populations.     

High level of unequal meiotic crossovers at the origin of the 22q11. 2 and 7q11.23 deletions

Interstitial chromosomal deletions at 22q11.2 and 7q11.23 are detected in the vast majority of patients affected by CATCH 22 syndromes and the Williams-Beuren syndrome, respectively. In a group of 15 Williams- Beuren patients, we have shown previously that a large number of 7q11.23 deletions occur in association with an interchromosomal rearrangement, indicative of an unequal crossing-over event between the two homologous chromosomes 7. In this study, we show that a similar mechanism also underlies the formation of the 22q11.2 deletions associated with CATCH 22. In eight out of 10 families with a proband affected by CATCH 22, we were able to show that a meiotic recombination had occurred at the critical deleted region based on segregation analysis of grandparental haplotypes. The incidences of crossovers observed between the closest informative markers, proximal and distal to the deletion, were compared with the expected recombination frequencies between the markers. A significant number of recombination events occur at the breakpoint of deletions in CATCH 22 patients (P = 2.99x10(-7)). The segregation analysis of haplotypes in three- generation families was also performed on an extended number of Williams-Beuren cases (22 cases in all). The statistically significant occurrence of meiotic crossovers (P = 4.45x10(-9)) further supports the previous findings. Thus, unequal meiotic crossover events appear to play a relevant role in the formation of the two interstitial deletions. The recurrence risk for healthy parents in cases where such meiotic recombinations can be demonstrated is probably negligible. Such a finding is in agreement with the predominantly sporadic occurrence of the 22q11.2 and 7q11. 23 deletions. No parent-of-origin bias was observed in the two groups of patients with regard to the origin of the deletion and to the occurrence of inter- versus intrachromosomal rearrangements.      

(CA/GT)(n) microsatellites affect homologous recombination during yeast meiosis

One of the most common microsatellites in eukaryotes consists of tandem arrays of the dinucleotide GT. Although the study of the instability of such repetitive DNA has been extremely fruitful over the last decade, no biological function has been demonstrated for these sequences. We investigated the genetic behavior of a region of the yeast Saccharomyces cerevisiae genome containing a 39-CA/GT dinucleotide repeat sequence. When the microsatellite sequence was present at the ARG4 locus on homologous chromosomes, diploid cells undergoing meiosis generated an excess of tetrads containing a conversion of the region restricted to the region of the microsatellite close to the recombination-initiation double-strand break. Moreover, whereas the repetitive sequence had no effect on the frequency of single crossover, its presence strongly stimulated the formation of multiple crossovers. The combined data strongly suggest that numerous recombination events are restricted to the initiation side of the microsatellite as though progression of the strand exchange initiated at the ARG4 promoter locus was impaired by the repetitive sequence. This observation corroborates in vitro experiments that demonstrated that RecA-promoted strand exchange is inhibited by CA/GT dinucleotide tracts. Surprisingly, meiotic instability of the microsatellite was very high (>0.1 alterations per tetrad) in all the spores with parental and recombinant chromosomes.     

Genome-wide variability in recombination activity is associated with meiotic chromatin organization

Recombination enables reciprocal exchange of genomic information between parental chromosomes and successful segregation of homologous chromosomes during meiosis. Errors in this process lead to negative health outcomes, whereas variability in recombination rate affects genome evolution. In mammals, most crossovers occur in hotspots defined by PRDM9 motifs, although PRDM9 binding peaks are not all equally hot. We hypothesize that dynamic patterns of meiotic genome folding are linked to recombination activity. We apply an integrative bioinformatics approach to analyze how three-dimensional (3D) chromosomal organization during meiosis relates to rates of double-strand-break (DSB) and crossover (CO) formation at PRDM9 binding peaks. We show that active, spatially accessible genomic regions during meiotic prophase are associated with DSB-favored loci, which further adopt a transient locally active configuration in early prophase. Conversely, crossover formation is depleted among DSBs in spatially accessible regions during meiotic prophase, particularly within gene bodies. We also find evidence that active chromatin regions have smaller average loop sizes in mammalian meiosis. Collectively, these findings establish that differences in chromatin architecture along chromosomal axes are associated with variable recombination activity. We propose an updated framework describing how 3D organization of brush-loop chromosomes during meiosis may modulate recombination.

The formation of crossovers during meiotic recombination is a highly orchestrated process, enhancing genetic diversity by allowing reciprocal exchange of genomic information to occur between parental chromosomes. Crossover formation also promotes proper segregation of homologous chromosomes (Baker et al. 1976), and errors in this process lead to chromosomal abnormalities such as aneuploidy, which are associated with negative health outcomes (Petronis 1999; Potapova and Gorbsky 2017). In mammals, crossovers are highly enriched (100-fold) in discrete ∼1- to 2-kb stretches along the genome, termed recombination hotspots (Paigen and Petkov 2010). These hotspots are in large part determined by the binding of PRDM9, a meiosis-specific zinc-finger protein that marks loci for potential recombination (Baudat et al. 2010; Myers et al. 2010; Parvanov et al. 2010).Although hotspot initiation is dependent on PRDM9, subsequent DSB and crossover formation are highly stochastic. Although exact numbers vary by species, a mammalian chromosome may harbor hundreds of PRDM9 binding loci, but during a typical meiotic cycle, only 10–20 double-stranded breaks (DSBs) occur (Diagouraga et al. 2018) per chromosome. Out of these DSBs, most are repaired as noncrossover conversion events, and only one or two per chromosome are chosen for crossover formation in mice (Baudat and de Massy 2007; Li et al. 2019). Local chromatin features such as GC content, histone modification, and cofactor binding are known to impact DSB formation at hotspots (Walker et al. 2015; Yamada et al. 2017), whereas nucleosome occupancy, GC content, and chromosomal position are associated with crossover formation (Hinch et al. 2019). Still, a full understanding of why certain hotspots are favored to form DSBs and crossovers remains undetermined.Meiotic chromosomes adopt a brush-loop conformation characterized by chromatin loops attached to a central axis (Møens and Pearlman 1988). Although recombination hotspots are found within loops, DSB machinery, such as DNA-repair proteins, resides on the axis (Blat et al. 2002; Grey et al. 2018; Tock and Henderson 2018; Slotman et al. 2020). This “tethered-loop/axis complex” model of recombination suggests that 3D genome folding could place constraints on the recombination process. Here we apply computational analyses to investigate how 3D chromatin organization relates to PRDM9 binding, DSBs, and crossover formation in male mammalian meiosis. Our analyses aim to integrate observations from multiple recent interphase and meiosis data sets measuring recombination activity and chromatin organization, including Hi-C, leading to an updated framework of how meiotic events related to recombination are associated with brush-loop chromosomal architecture.     

Allelic recombination and de novo deletions in sperm in the human beta-globin gene region

Meiotic recombination is of fundamental importance in creating haplotype diversity in the human genome and has the potential to cause genomic rearrangements by ectopic recombination between repeat sequences and through other changes triggered by recombination-initiating events. However, the relationship between allelic recombination and genome instability in the human germline remains unclear. We have therefore analysed recombination and DNA instability in the delta-, beta-globin gene region and its associated recombination hotspot. Sperm typing has for the first time accurately defined the hotspot and shown it to be the most active autosomal crossover hotspot yet described, although unusually inactive in non-exchange gene conversion. The hotspot just extends into a homology block shared by the delta- and beta-globin genes, within which ectopic exchanges can generate Hb Lepore deletions. We developed a physical selection method for recovering and validating extremely rare de novo deletions in human DNA and used it to characterize the dynamics of these Hb Lepore deletions in sperm as well as other deletions not arising from ectopic exchanges between homologous DNA sequences. Surprisingly, both classes of deletion showed breakpoints that avoided the beta-globin hotspot, establishing that it possesses remarkable fidelity and does not play a significant role in triggering these DNA rearrangements. This study also provides the first direct analysis of de novo deletion in the human germline and points to a possible deletion-controlling element in the beta-globin gene separate from the crossover hotspot.     

A clustering of RNA recombination sites adjacent to a hypervariable region of the peplomer gene of murine coronavirus

Coronaviruses undergo RNA recombination at a very high frequency. To understand the mechanism of recombination in murine coronavirus, we have performed RNA sequencing of viral genomic RNA to determine the precise sites of recombination in a series of recombinants which have crossovers within the gene encoding the peplomer protein. We found that all of the recombination sites are clustered within a region of 278 nucleotides in the 5'-half of the gene. This region in which all of the crossovers occurred represents a small fraction of the distance between the two selection markers used for the isolation of these recombinant viruses. This result suggests that this region may be a preferred site for RNA recombination. The crossover sites are located within and immediately adjacent to a hypervariable area of the gene. This area has undergone deletions of varying sizes in several virus strains which have been passaged either in vivo or in vitro. These results suggest that a similar RNA structure may be involved in the occurrence of both recombination and deletion events.     

Polarity in adenovirus recombination

The distributions of the crossovers necessary to generate ts+ genomes have been examined in a collection of clonally unrelated ts+ recombinants from a set of ts X ts adenovirus crosses. In a cross between two parents that are grossly heterologous between map units 80.2 and 91.5, the distribution of crossovers was significantly skewed toward the left-hand end of the genome, with a declining frequency proceeding rightward. This gradient of recombination was modified by the removal of the right-hand heterology and by the presence of another region of heterology between map units 3.67 and 10.11. In a cross where the ts markers were flanked by both heterologies, no gradient was observed and ts+ recombinants were characterized by a higher rate of supernumerary crossovers. In a cross designed so that one ts marker was internal to two heterologies, crossovers were found disproportionately between the second ts marker and the nearby heterology. In addition, ts+ recombinants formed by crossing over internal to the heterologies again were accompanied by a high frequency of supernumerary crossovers. Finally, ts+ recombinant frequencies in crosses identical except for the presence of either one or two flanking heterologies were markedly lower in the latter case. These data, taken together, suggest that a major pathway of adenovirus recombination initiates at, or near, the molecular termini and is perhaps driven by the displaced single strands produced during DNA replication. Internal initiation, on the other hand, may employ these single strands to form genetic "patches."     

Discontinuities and unsynapsed regions in meiotic chromosomes have a cis effect on meiotic recombination patterns in normal human males

During meiosis, homologous chromosome pairing is essential for subsequent meiotic recombination (crossover). Discontinuous chromosome regions (gaps) or unsynapsed chromosome regions (splits) in the synaptonemal complex (SC) indicate anomalies in chromosome synapsis. Recently developed immunofluorescence techniques (using antibodies against SC proteins and the crossover-associated MLH1 protein) were combined with fluorescence in situ hybridization (using centromere-specific DNA probes) to identify bivalents with gaps/splits and to examine the effect of gaps/splits on meiotic recombination patterns during the pachytene stage of meiotic prophase from three normal human males. Gaps were observed only in the heterochromatic regions of chromosomes 9 and 1, with 9q gaps accounting for 90% of these events. Most splits were also found in chromosomes 9 and 1, with 58% of splits occurring on 9q. Gaps and splits significantly altered the distribution of MLH1 foci on the SC. On gapped SC 9q, the frequency of MLH1 foci was decreased compared with controls, and single 9q crossovers tended toward a more distal distribution. Furthermore, the larger the gap the more distal the location of the MLH1 focus closest to the q arm's telomere. MLH1 foci on split SC 9 had distributions similar to those of gapped SC 9; however, splits did not change the frequencies of MLH1 foci on SC 9. This is the first demonstration that gaps and splits have an effect on meiotic recombination in humans.     

Characterization of susceptible chiasma configurations that increase the risk for maternal nondisjunction of chromosome 21

Recent studies of trisomy 21 have shown that altered levels of recombination are associated with maternal non-disjunction occurring at both meiosis I (MI) and meiosis II (MII). To comprehend better the association of recombination with nondisjunction, an understanding of the pattern of meiotic exchange, i.e. the exchange of genetic material at the four-strand stage during prophase, is required. We examined this underlying exchange pattern to determine if specific meiotic configurations are associated with a higher risk of non-disjunction than others. We examined the crossover frequencies of chromosome 21 for three populations: (i) normal female meiotic events; (ii) meiotic events leading to MI non-disjunction; and (iii) those leading to MII non-disjunction. From these crossover frequencies, we estimated the array of meiotic tetrads that produced the observed crossovers. Using this approach, we found that nearly one-half of MI errors were estimated to be achiasmate. The majority of the remaining MI bivalents had exchanges that clustered at the telomere. In contrast, exchanges occurring among MII cases clustered at the pericentromeric region of the chromosome. Unlike the single exchange distributions, double exchanges from the non-disjoined populations seemed to approximate the distribution in the normal population. These data suggest that the location of certain exchanges makes a tetrad susceptible to non- disjunction. Specifically, this susceptibility is associated with the distance between the centromere and closest exchange. This result challenges the widely held concept that events occurring at MII are largely independent of events occurring at MI, and suggests that all non-disjunction events may be initiated during MI and simply resolved at either of the two meiotic stages.      

Analysis of Arabidopsis genome-wide variations before and after meiosis and meiotic recombination by resequencing Landsberg erecta and all four products of a single meiosis

Meiotic recombination, including crossovers (COs) and gene conversions (GCs), impacts natural variation and is an important evolutionary force. COs increase genetic diversity by redistributing existing variation, whereas GCs can alter allelic frequency. Here, we sequenced Arabidopsis Landsberg erecta (Ler) and two sets of all four meiotic products from a Columbia (Col)/Ler hybrid to investigate genome-wide variation and meiotic recombination at nucleotide resolution. Comparing Ler and Col sequences uncovered 349,171 Single Nucleotide Polymorphisms (SNPs), 58,085 small and 2315 large insertions/deletions (indels), with highly correlated genome-wide distributions of SNPs, and small indels. A total of 443 genes have at least 10 nonsynonymous substitutions in protein-coding regions, with enrichment for disease-resistance genes. Another 316 genes are affected by large indels, including 130 genes with complete deletion of coding regions in Ler. Using the Arabidopsis qrt1 mutant, two sets of four meiotic products were generated and analyzed by sequencing for meiotic recombination, representing the first tetrad analysis with whole-genome sequencing in a nonfungal species. We detected 18 COs, six of which had an associated GC event, and four GCs without COs (NCOs), and revealed that Arabidopsis GCs are likely fewer and with shorter tracts than those in yeast. Meiotic recombination and chromosome assortment events dramatically redistributed genome variation in meiotic products, contributing to population diversity. In particular, meiosis provides a rapid mechanism to generate copy-number variation (CNV) of sequences that have different chromosomal positions in Col and Ler.     

Inference and analysis of haplotypes from combined genotyping studies deposited in dbSNP

In the attempt to understand human variation and the genetic basis of complex disease, a tremendous number of single nucleotide polymorphisms (SNPs) have been discovered and deposited into NCBI's dbSNP public database. More than 2.7 million SNPs in the database have genotype information. This data provides an invaluable resource for understanding the structure of human variation and the design of genetic association studies. The genotypes deposited to dbSNP are unphased, and thus, the haplotype information is unknown. We applied the phasing method HAP to obtain the haplotype information, block partitions, and tag SNPs for all publicly available genotype data and deposited this information into the dbSNP database. We also deposited the orthologous chimpanzee reference sequence for each predicted haplotype block computed using the UCSC BLASTZ alignments of human and chimpanzee. Using dbSNP, researchers can now easily perform analyses using multiple genotype data sets from the same genomic regions. Dense and sparse genotype data sets from the same region were combined to show that the number of common haplotypes is significantly underestimated in whole genome data sets, while the predicted haplotypes over the common SNPs are consistent between studies. To validate the accuracy of the predictions, we bench-marked HAP's running time and phasing accuracy against PHASE. Although HAP is slightly less accurate than PHASE, HAP is over 1000 times faster than PHASE, making it suitable for application to the entire set of genotypes in dbSNP.     

Differential Divergence of Three Human Pseudoautosomal Genes and Their Mouse Homologs: Implications for Sex Chromosome Evolution

The human pseudoautosomal region 1 (PAR1) is essential for meiotic pairing and recombination, and its deletion causes male sterility. Comparative studies of human and mouse pseudoautosomal genes are valuable in charting the evolution of this interesting region, but have been limited by the paucity of genes conserved between the two species. We have cloned a novel human PAR1 gene, DHRSXY, encoding an oxidoreductase of the short-chain dehydrogenase/reductase family, and isolated a mouse ortholog Dhrsxy. We also searched for mouse homologs of recently reported PGPL and TRAMP genes that flank it within PAR1. We recovered a highly conserved mouse ortholog of PGPL by cross-hybridization, but found no mouse homolog of TRAMP. Like Csf2ra and Il3ra, both mouse homologs are autosomal; Pgpl on chromosome 5, and Dhrsxy subtelomeric on chromosome 4. TRAMP, like the human genes within or near PAR1, is probably very divergent or absent in the mouse genome. We interpret the rapid divergence and loss of pseudoautosomal genes in terms of a model of selection for the concentration of repetitive recombinogenic sequences that predispose to high recombination and translocation.