Evidence for heterogeneity in recombination in the human pseudoautosomal region: high resolution analysis by sperm typing and radiation-hybrid mapping

Genomic and demographic processes differentially influence genetic variation across the human X chromosome

Many forces influence genetic variation across the genome including mutation, recombination, selection, and demography. Increased mutation and recombination both lead to increases in genetic diversity in a region-specific manner, while complex demographic patterns shape patterns of diversity on a more global scale. While these processes act across the entire genome, the X chromosome is particularly interesting because it contains several distinct regions that are subject to different combinations and strengths of these forces: the pseudoautosomal regions (PARs) and the X-transposed region (XTR). The X chromosome thus can serve as a unique model for studying how genetic and demographic forces act in different contexts to shape patterns of observed variation. We therefore sought to explore diversity, divergence, and linkage disequilibrium in each region of the X chromosome using genomic data from 26 human populations. Across populations, we find that both diversity and substitution rate are consistently elevated in PAR1 and the XTR compared to the rest of the X chromosome. In contrast, linkage disequilibrium is lowest in PAR1, consistent with the high recombination rate in this region, and highest in the region of the X chromosome that does not recombine in males. However, linkage disequilibrium in the XTR is intermediate between PAR1 and the autosomes, and much lower than the non-recombining X. Finally, in addition to these global patterns, we also observed variation in ratios of X versus autosomal diversity consistent with population-specific evolutionary history as well. While our results were generally consistent with previous work, two unexpected observations emerged. First, our results suggest that the XTR does not behave like the rest of the recombining X and may need to be evaluated separately in future studies. Second, the different regions of the X chromosome appear to exhibit unique patterns of linked selection across different human populations. Together, our results highlight profound regional differences across the X chromosome, simultaneously making it an ideal system for exploring the action of evolutionary forces as well as necessitating its careful consideration and treatment in genomic analyses.

Deviations from Mendelian Inheritance on Bovine X-Chromosome Revealing Recombination, Sex-of-Offspring Effects and Fertility-Related Candidate Genes

Transmission ratio distortion (TRD), or significant deviations from Mendelian inheritance, is a well-studied phenomenon on autosomal chromosomes, but has not yet received attention on sex chromosomes. TRD was analyzed on 3832 heterosomal single nucleotide polymorphisms (SNPs) and 400 pseudoautosomal SNPs spanning the length of the X-chromosome using 436,651 genotyped Holstein cattle. On the pseudoautosomal region, an opposite sire-TRD pattern between male and female offspring was identified for 149 SNPs. This finding revealed unique SNPs linked to a specific-sex (Y- or X-) chromosome and describes the accumulation of recombination events across the pseudoautosomal region. On the heterosomal region, 13 SNPs and 69 haplotype windows were identified with dam-TRD. Functional analyses for TRD regions highlighted relevant biological functions responsible to regulate spermatogenesis, development of Sertoli cells, homeostasis of endometrium tissue and embryonic development. This study uncovered the prevalence of different TRD patterns across both heterosomal and pseudoautosomal regions of the X-chromosome and revealed functional candidate genes for bovine reproduction.

Recent expansion of the non-recombining sex-linked region on Silene latifolia sex chromosomes

Evolution of a non-recombining sex-specific region on the Y (or W) chromosome (NRY) is a key step in sex chromosome evolution, but how recombination suppression evolves is not well understood. Studies in many different organisms indicated that NRY evolution often involves several expansion steps. Why such NRY expansions occur remains unclear, although it is though that they are likely driven by sexually antagonistic selection. This paper describes a recent NRY expansion due to shift of the pseudoautosomal boundary on the sex chromosomes of a dioecious plant Silene latifolia. The shift resulted in inclusion of at least 16 pseudoautosomal genes into the NRY. This region is pseudoautosomal in closely related Silene dioica and Silene diclinis, indicating that the NRY expansion occurred in S. latifolia after it speciated from the other species ~120 thousand years ago. As S. latifolia and S. dioica actively hybridise across Europe, interspecific gene flow could blur the PAR boundary in these species. The pseudoautosomal genes have significantly elevated genetic diversity (π ~ 3% at synonymous sites), which is consistent with balancing selection maintaining diversity in this region. The recent shift of the PAR boundary in S. latifolia offers an opportunity to study the process of on-going NRY expansion.

Decoding sex: Elucidating sex determination and how high-quality genome assemblies are untangling the evolutionary dynamics of sex chromosomes

Sexual reproduction is a diverse and widespread process. In gonochoristic species, the differentiation of sexes occurs through diverse mechanisms, influenced by environmental and genetic factors. In most vertebrates, a master-switch gene is responsible for triggering a sex determination network. However, only a few genes have acquired master-switch functions, and this process is associated with the evolution of sex-chromosomes, which have a significant influence in evolution. Additionally, their highly repetitive regions impose challenges for high-quality sequencing, even using high-throughput, state-of-the-art techniques. Here, we review the mechanisms involved in sex determination and their role in the evolution of species, particularly vertebrates, focusing on sex chromosomes and the challenges involved in sequencing these genomic elements. We also address the improvements provided by the growth of sequencing projects, by generating a massive number of near-gapless, telomere-to-telomere, chromosome-level, phased assemblies, increasing the number and quality of sex-chromosome sequences available for further studies.

Chromosome-Level Genome Assemblies Expand Capabilities of Genomics for Conservation Biology

Genome assemblies are in the process of becoming an increasingly important tool for understanding genetic diversity in threatened species. Unfortunately, due to limited budgets typical for the area of conservation biology, genome assemblies of threatened species, when available, tend to be highly fragmented, represented by tens of thousands of scaffolds not assigned to chromosomal locations. The recent advent of high-throughput chromosome conformation capture (Hi-C) enables more contiguous assemblies containing scaffolds spanning the length of entire chromosomes for little additional cost. These inexpensive contiguous assemblies can be generated using Hi-C scaffolding of existing short-read draft assemblies, where N50 of the draft contigs is larger than 0.1% of the estimated genome size and can greatly improve analyses and facilitate visualization of genome-wide features including distribution of genetic diversity in markers along chromosomes or chromosome-length scaffolds. We compared distribution of genetic diversity along chromosomes of eight mammalian species, including six listed as threatened by IUCN, where both draft genome assemblies and newer chromosome-level assemblies were available. The chromosome-level assemblies showed marked improvement in localization and visualization of genetic diversity, especially where the distribution of low heterozygosity across the genomes of threatened species was not uniform.

Evolutionary dynamics of the human pseudoautosomal regions

Recombination between the X and Y human sex chromosomes is limited to the two pseudoautosomal regions (PARs) that present quite distinct evolutionary origins. Despite the crucial importance for male meiosis, genetic diversity patterns and evolutionary dynamics of these regions are poorly understood. In the present study, we analyzed and compared the genetic diversity of the PAR regions using publicly available genomic sequences encompassing both PAR1 and PAR2. Comparisons were performed through allele diversities, linkage disequilibrium status and recombination frequencies within and between X and Y chromosomes. In agreement with previous studies, we confirmed the role of PAR1 as a male-specific recombination hotspot, but also observed similar characteristic patterns of diversity in both regions although male recombination occurs at PAR2 to a much lower extent (at least one recombination event at PAR1 and in ≈1% in normal male meioses at PAR2). Furthermore, we demonstrate that both PARs harbor significantly different allele frequencies between X and Y chromosomes, which could support that recombination is not sufficient to homogenize the pseudoautosomal gene pool or is counterbalanced by other evolutionary forces. Nevertheless, the observed patterns of diversity are not entirely explainable by sexually antagonistic selection. A better understanding of such processes requires new data from intergenerational transmission studies of PARs, which would be decisive on the elucidation of PARs evolution and their role in male-driven heterosomal aneuploidies.

The Location of the Pseudoautosomal Boundary in Silene latifolia

Y-chromosomes contain a non-recombining region (NRY), and in many organisms it was shown that the NRY expanded over time. How and why the NRY expands remains unclear. Young sex chromosomes, where NRY expansion occurred recently or is on-going, offer an opportunity to study the causes of this process. Here, we used the plant Silene latifolia, where sex chromosomes evolved ~11 million years ago, to study the location of the boundary between the NRY and the recombining pseudoautosomal region (PAR). The previous work devoted to the NRY/PAR boundary in S. latifolia was based on a handful of genes with locations approximately known from the genetic map. Here, we report the analysis of 86 pseudoautosomal and sex-linked genes adjacent to the S. latifolia NRY/PAR boundary to establish the location of the boundary more precisely. We take advantage of the dense genetic map and polymorphism data from wild populations to identify 20 partially sex-linked genes located in the “fuzzy boundary”, that rarely recombines in male meiosis. Genes proximal to this fuzzy boundary show no evidence of recombination in males, while the genes distal to this partially-sex-linked region are actively recombining in males. Our results provide a more accurate location for the PAR boundary in S. latifolia, which will help to elucidate the causes of PAR boundary shifts leading to NRY expansion over time.

Both male and female gamete generating cells produce processed pseudogenes in the human genome

The human genome contains an unusually large number of processed pseudogenes. The fact that processed pseudogenes are roughly 33% more abundant in our X chromosome than in our autosomes suggests that this overabundance is the result of the fact that human oogenesis is much longer than that of non-mammalian species. Here, we analyze the origins of the processed pseudogenes found on the human Y chromosome to determine whether human spermatogenesis also contribute to this overabundance. Our results show that human processed pseudogenes not only retrotranspose to the Y chromosome, but are also produced by genes on the Y chromosome. Furthermore, the fact that X chromosomes are three times more abundant than Y chromosomes likely explains why the euchromatic density of processed pseudogenes is three times higher in the X chromosome than in the Y chromosome. The large number of processed pseudogenes found in our genome is therefore due to the low substrate specificity of the L1 reverse transcriptase responsible for the reverse transcription of germline mRNA molecules into processed pseudogenes, as well as the life-long production of both male and female gametes.

Identification of a novel 15.5 kb SHOX deletion associated with marked intrafamilial phenotypic variability and analysis of its molecular origin

Haploinsufficiency of the short stature homeobox contaning SHOX gene has been shown to result in a spectrum of phenotypes ranging from Leri-Weill dyschondrosteosis (LWD) at the more severe end to SHOX-related short stature at the milder end of the spectrum. Most alterations are whole gene deletions, point mutations within the coding region, or microdeletions in its flanking sequences. Here, we present the clinical and molecular data as well as the potential molecular mechanism underlying a novel microdeletion, causing a variable SHOX-related haploinsufficiency disorder in a three-generation family. The phenotype resembles that of LWD in females, in males, however, the phenotypic expression is milder. The 15523-bp SHOX intragenic deletion, encompassing exons 3-6, was initially detected by array-CGH, followed by MLPA analysis. Sequencing of the breakpoints indicated an Alu recombination-mediated deletion (ARMD) as the potential causative mechanism.

A Track Record on SHOX: From Basic Research to Complex Models and Therapy

SHOX deficiency is the most frequent genetic growth disorder associated with isolated and syndromic forms of short stature. Caused by mutations in the homeobox gene SHOX, its varied clinical manifestations include isolated short stature, Léri-Weill dyschondrosteosis, and Langer mesomelic dysplasia. In addition, SHOX deficiency contributes to the skeletal features in Turner syndrome. Causative SHOX mutations have allowed downstream pathology to be linked to defined molecular lesions. Expression levels of SHOX are tightly regulated, and almost half of the pathogenic mutations have affected enhancers. Clinical severity of SHOX deficiency varies between genders and ranges from normal stature to profound mesomelic skeletal dysplasia. Treatment options for children with SHOX deficiency are available. Two decades of research support the concept of SHOX as a transcription factor that integrates diverse aspects of bone development, growth plate biology, and apoptosis. Due to its absence in mouse, the animal models of choice have become chicken and zebrafish. These models, therefore, together with micromass cultures and primary cell lines, have been used to address SHOX function. Pathway and network analyses have identified interactors, target genes, and regulators. Here, we summarize recent data and give insight into the critical molecular and cellular functions of SHOX in the etiopathogenesis of short stature and limb development.

Determinants of genetic diversity

Genetic polymorphism varies among species and within genomes, and has important implications for the evolution and conservation of species. The determinants of this variation have been poorly understood, but population genomic data from a wide range of organisms now make it possible to delineate the underlying evolutionary processes, notably how variation in the effective population size (Ne) governs genetic diversity. Comparative population genomics is on its way to providing a solution to 'Lewontin's paradox' - the discrepancy between the many orders of magnitude of variation in population size and the much narrower distribution of diversity levels. It seems that linked selection plays an important part both in the overall genetic diversity of a species and in the variation in diversity within the genome. Genetic diversity also seems to be predictable from the life history of a species.

Hybrid Sterility in Rice (Oryza sativa L.) Involves the Tetratricopeptide Repeat Domain Containing Protein

Intersubspecific hybrid sterility is a common form of reproductive isolation in rice (Oryza sativa L.), which significantly hampers the utilization of heterosis between indica and japonica varieties. Here, we elucidated the mechanism of S7, which specially causes Aus-japonica/indica hybrid female sterility, through cytological and genetic analysis, map-based cloning, and transformation experiments. Abnormal positioning of polar nuclei and smaller embryo sac were observed in F1 compared with male and female parents. Female gametes carrying S7(cp) and S7(i) were aborted in S7(ai)/S7(cp) and S7(ai)/S7(i), respectively, whereas they were normal in both N22 and Dular possessing a neutral allele, S7(n) S7 was fine mapped to a 139-kb region in the centromere region on chromosome 7, where the recombination was remarkably suppressed due to aggregation of retrotransposons. Among 16 putative open reading frames (ORFs) localized in the mapping region, ORF3 encoding a tetratricopeptide repeat domain containing protein was highly expressed in the pistil. Transformation experiments demonstrated that ORF3 is the candidate gene: downregulated expression of ORF3 restored spikelet fertility and eliminated absolutely preferential transmission of S7(ai) in heterozygote S7(ai)/S7(cp); sterility occurred in the transformants Cpslo17-S7(ai) Our results may provide implications for overcoming hybrid embryo sac sterility in intersubspecific hybrid rice and utilization of hybrid heterosis for cultivated rice improvement.

Genetic Diversity on the Human X Chromosome Does Not Support a Strict Pseudoautosomal Boundary

Unlike the autosomes, recombination between the X chromosome and the Y chromosome is often thought to be constrained to two small pseudoautosomal regions (PARs) at the tips of each sex chromosome. PAR1 spans the first 2.7 Mb of the proximal arm of the human sex chromosomes, whereas the much smaller PAR2 encompasses the distal 320 kb of the long arm of each sex chromosome. In addition to PAR1 and PAR2, there is a human-specific X-transposed region that was duplicated from the X to the Y chromosome. The X-transposed region is often not excluded from X-specific analyses, unlike the PARs, because it is not thought to routinely recombine. Genetic diversity is expected to be higher in recombining regions than in nonrecombining regions because recombination reduces the effect of linked selection. In this study, we investigated patterns of genetic diversity in noncoding regions across the entire X chromosome of a global sample of 26 unrelated genetic females. We found that genetic diversity in PAR1 is significantly greater than in the nonrecombining regions (nonPARs). However, rather than an abrupt drop in diversity at the pseudoautosomal boundary, there is a gradual reduction in diversity from the recombining through the nonrecombining regions, suggesting that recombination between the human sex chromosomes spans across the currently defined pseudoautosomal boundary. A consequence of recombination spanning this boundary potentially includes increasing the rate of sex-linked disorders (e.g., de la Chapelle) and sex chromosome aneuploidies. In contrast, diversity in PAR2 is not significantly elevated compared to the nonPARs, suggesting that recombination is not obligatory in PAR2. Finally, diversity in the X-transposed region is higher than in the surrounding nonPARs, providing evidence that recombination may occur with some frequency between the X and Y chromosomes in the X-transposed region.

The Pseudoautosomal Regions of the U/V Sex Chromosomes of the Brown Alga Ectocarpus Exhibit Unusual Features

The recombining regions of sex chromosomes (pseudoautosomal regions, PARs) are predicted to exhibit unusual features due to their being genetically linked to the nonrecombining, sex-determining region. This phenomenon is expected to occur in both diploid (XY, ZW) and haploid (UV) sexual systems, with slightly different consequences for UV sexual systems because of the absence of masking during the haploid phase (when sex is expressed) and because there is no homozygous sex in these systems. Despite a considerable amount of theoretical work on PAR genetics and evolution, these genomic regions have remained poorly characterized empirically. We show here that although the PARs of the U/V sex chromosomes of the brown alga Ectocarpus recombine at a similar rate to autosomal regions of the genome, they exhibit many genomic features typical of nonrecombining regions. The PARs were enriched in clusters of genes that are preferentially, and often exclusively, expressed during the sporophyte generation of the life cycle, and many of these genes appear to have evolved since the Ectocarpales diverged from other brown algal lineages. A modeling-based approach was used to investigate possible evolutionary mechanisms underlying this enrichment in sporophyte-biased genes. Our results are consistent with the evolution of the PAR in haploid systems being influenced by differential selection pressures in males and females acting on alleles that are advantageous during the sporophyte generation of the life cycle.

Recombination in the human Pseudoautosomal region PAR1

The pseudoautosomal region (PAR) is a short region of homology between the mammalian X and Y chromosomes, which has undergone rapid evolution. A crossover in the PAR is essential for the proper disjunction of X and Y chromosomes in male meiosis, and PAR deletion results in male sterility. This leads the human PAR with the obligatory crossover, PAR1, to having an exceptionally high male crossover rate, which is 17-fold higher than the genome-wide average. However, the mechanism by which this obligatory crossover occurs remains unknown, as does the fine-scale positioning of crossovers across this region. Recent research in mice has suggested that crossovers in PAR may be mediated independently of the protein PRDM9, which localises virtually all crossovers in the autosomes. To investigate recombination in this region, we construct the most fine-scale genetic map containing directly observed crossovers to date using African-American pedigrees. We leverage recombination rates inferred from the breakdown of linkage disequilibrium in human populations and investigate the signatures of DNA evolution due to recombination. Further, we identify direct PRDM9 binding sites using ChIP-seq in human cells. Using these independent lines of evidence, we show that, in contrast with mouse, PRDM9 does localise peaks of recombination in the human PAR1. We find that recombination is a far more rapid and intense driver of sequence evolution in PAR1 than it is on the autosomes. We also show that PAR1 hotspot activities differ significantly among human populations. Finally, we find evidence that PAR1 hotspot positions have changed between human and chimpanzee, with no evidence of sharing among the hottest hotspots. We anticipate that the genetic maps built and validated in this work will aid research on this vital and fascinating region of the genome.

Variation in genome-wide levels of meiotic recombination is established at the onset of prophase in mammalian males

Segregation of chromosomes during the first meiotic division relies on crossovers established during prophase. Although crossovers are strictly regulated so that at least one occurs per chromosome, individual variation in crossover levels is not uncommon. In an analysis of different inbred strains of male mice, we identified among-strain variation in the number of foci for the crossover-associated protein MLH1. We report studies of strains with “low” (CAST/EiJ), “medium” (C3H/HeJ), and “high” (C57BL/6J) genome-wide MLH1 values to define factors responsible for this variation. We utilized immunofluorescence to analyze the number and distribution of proteins that function at different stages in the recombination pathway: RAD51 and DMC1, strand invasion proteins acting shortly after double-strand break (DSB) formation, MSH4, part of the complex stabilizing double Holliday junctions, and the Bloom helicase BLM, thought to have anti-crossover activity. For each protein, we identified strain-specific differences that mirrored the results for MLH1; i.e., CAST/EiJ mice had the lowest values, C3H/HeJ mice intermediate values, and C57BL/6J mice the highest values. This indicates that differences in the numbers of DSBs (as identified by RAD51 and DMC1) are translated into differences in the number of crossovers, suggesting that variation in crossover levels is established by the time of DSB formation. However, DSBs per se are unlikely to be the primary determinant, since allelic variation for the DSB-inducing locus Spo11 resulted in differences in the numbers of DSBs but not the number of MLH1 foci. Instead, chromatin conformation appears to be a more important contributor, since analysis of synaptonemal complex length and DNA loop size also identified consistent strain-specific differences; i.e., crossover frequency increased with synaptonemal complex length and was inversely related to chromatin loop size. This indicates a relationship between recombination and chromatin compaction that may develop as DSBs form or earlier during establishment of the meiotic axis.

During prophase of meiosis, homologous chromosomes exchange genetic material, in a process known as crossing-over. Crossovers are thought to be essential for proper separation of chromosomes during meiosis but, surprisingly, most mammalian species exhibit substantial individual variation in the number of crossovers per cell. We investigated the basis for this variation by examining localization patterns of crossover-associated proteins in inbred strains of male mice with differing average numbers of crossovers per spermatocyte. Our results indicate that the strain-specific variation is established early in meiotic prophase, possibly even before the DNA is broken in advance of subsequent exchanges between homologous chromosomes.

Expansion of the pseudo-autosomal region and ongoing recombination suppression in the Silene latifolia sex chromosomes

There are two very interesting aspects to the evolution of sex chromosomes: what happens after recombination between these chromosome pairs stops and why suppressed recombination evolves. The former question has been intensively studied in a diversity of organisms, but the latter has been studied largely theoretically. To obtain empirical data, we used codominant genic markers in genetic mapping of the dioecious plant Silene latifolia, together with comparative mapping of S. latifolia sex-linked genes in S. vulgaris (a related hermaphrodite species without sex chromosomes). We mapped 29 S. latifolia fully sex-linked genes (including 21 newly discovered from transcriptome sequencing), plus 6 genes in a recombining pseudo-autosomal region (PAR) whose genetic map length is ∼25 cM in both male and female meiosis, suggesting that the PAR may contain many genes. Our comparative mapping shows that most fully sex-linked genes in S. latifolia are located on a single S. vulgaris linkage group and were probably inherited from a single autosome of an ancestor. However, unexpectedly, our maps suggest that the S. latifolia PAR region expanded through translocation events. Some genes in these regions still recombine in S. latifolia, but some genes from both addition events are now fully sex-linked. Recombination suppression is therefore still ongoing in S. latifolia, and multiple recombination suppression events have occurred in a timescale of few million years, much shorter than the timescale of formation of the most recent evolutionary strata of mammal and bird sex chromosomes.

Clinical and molecular characterization of Chilean patients with Léri-Weill dyschondrosteosis

AIM

Léri-Weill dyschondrosteosis (LWD) is a mesomelic dysplasia with disproportionate short stature associated with short stature homeobox-containing gene (SHOX) haploinsufficiency. The objective of this study was to improve the diagnosis of patients with suspected LWD through molecular analysis.

METHODS

Twelve patients from 11 families with a clinical diagnosis of LWD were analyzed with multiplex ligation-dependent probe amplification to detect deletions and duplications of SHOX and its enhancer regions. High resolution melting and sequencing was employed to screen for mutations in SHOX coding exons.

RESULTS

The molecular-based screening strategy applied in these patients allowed detection of five SHOX deletions and two previously unreported SHOX missense mutations.

CONCLUSION

Molecular studies confirmed the clinical diagnosis of LWD in seven out of 12 patients, which provided support for therapeutic decisions and improved genetic counseling in their families.

Sex modifies genetic effects on residual variance in urinary calcium excretion in rat (Rattus norvegicus)

Conventional genetics assumes common variance among alleles or genetic groups. However, evidence from vertebrate and invertebrate models suggests that residual genotypic variance may itself be under partial genetic control. Such a phenomenon would have great significance: high-variability alleles might confound the detection of "classically" acting genes or scatter predicted evolutionary outcomes among unpredicted trajectories. Of the few works on this phenomenon, many implicate sex in some aspect of its control. We found that female genetic hypercalciuric stone-forming (GHS) rats (Rattus norvegicus) had higher coefficients of variation (CVs) for urinary calcium (CV = 0.14) than GHS males (CV = 0.06), and the reverse in normocalciuric Wistar-Kyoto rats (WKY) (CV(♂) = 0.14; CV(♀) = 0.09), suggesting sex-by-genotype interaction on residual variance. We therefore investigated the effect of sex on absolute-transformed residuals in urinary calcium in an F(2) GHS × WKY mapping cohort. Absolute residuals were associated with genotype at two microsatellites, D3Rat46 (RNO3, 33.9 Mb) and D4Mgh1 (RNO4, 84.8 MB) at Bonferroni thresholds across the entire cohort, and with the microsatellites D3Rat46, D9Mgh2 (RNO9, 84.4 Mb), and D12Rat25 (RNO12, 40.4 Mb) in females (P < 0.05) but not males. In GHS chromosome 1 congenic lines bred onto a WKY genomic background, we found that congenic males had significantly (P < 0.0001) higher CVs for urinary calcium (CV = 0.25) than females (CV = 0.15), supporting the hypothesis of the inheritance of sex-by-genotype interaction on this effect. Our findings suggest that genetic effects on residual variance are sex linked; heritable, sex-specific residuals might have great potential implications for evolution, adaptation, and genetic analysis.

Fine scale analysis of crossover and non-crossover and detection of recombination sequence motifs in the honeybee (Apis mellifera)

BACKGROUND

Meiotic exchanges are non-uniformly distributed across the genome of most studied organisms. This uneven distribution suggests that recombination is initiated by specific signals and/or regulations. Some of these signals were recently identified in humans and mice. However, it is unclear whether or not sequence signals are also involved in chromosomal recombination of insects.

METHODOLOGY

We analyzed recombination frequencies in the honeybee, in which genome sequencing provided a large amount of SNPs spread over the entire set of chromosomes. As the genome sequences were obtained from a pool of haploid males, which were the progeny of a single queen, an oocyte method (study of recombination on haploid males that develop from unfertilized eggs and hence are the direct reflect of female gametes haplotypes) was developed to detect recombined pairs of SNP sites. Sequences were further compared between recombinant and non-recombinant fragments to detect recombination-specific motifs.

CONCLUSIONS

Recombination events between adjacent SNP sites were detected at an average distance of 92 bp and revealed the existence of high rates of recombination events. This study also shows the presence of conversion without crossover (i. e. non-crossover) events, the number of which largely outnumbers that of crossover events. Furthermore the comparison of sequences that have undergone recombination with sequences that have not, led to the discovery of sequence motifs (CGCA, GCCGC, CCGCA), which may correspond to recombination signals.

The jumping SHOX gene--crossover in the pseudoautosomal region resulting in unusual inheritance of Leri-Weill dyschondrosteosis

CONTEXT

During meiosis I, the recombination frequency in the pseudoautosomal region on Xp and Yp (PAR1) in males is very high. As a result, mutated genes located within the PAR1 region can be transferred from the Y-chromosome to the X-chromosome and vice versa.

PATIENTS

Here we describe three families with SHOX abnormalities resulting in Leri-Weill dyschondrosteosis or Langer mesomelic dysplasia.

RESULTS

In about half of the segregations investigated, a transfer of the SHOX abnormality to the alternate sex chromosome was demonstrated.

CONCLUSIONS

Patients with an abnormality of the SHOX gene should receive genetic counseling as to the likelihood that they may transmit the mutation or deletion to a son as well as to a daughter.

Identification of the first de novo PAR1 deletion downstream of SHOX in an individual diagnosed with Léri-Weill dyschondrosteosis (LWD)

Léri-Weill dyschondrosteosis (LWD, MIM 127300), is a dominantly inherited skeletal dysplasia with disproportionate short stature, mesomelic limb shortening, and the characteristic Madelung deformity. Two regions of the pseudoautosomal region 1 (PAR1) have been shown to be involved in LWD, SHOX (short-stature homeobox-containing gene) and the downstream enhancer region. We report our genetic findings of a young girl clinically diagnosed with LWD. We analyzed the proband and her family using MLPA and microsatellite analysis. We identified a deletion, 726-866 kb in size, of the downstream SHOX enhancer region in the proband. Neither parent carried the deletion. Microsatellite analysis showed that the deleted allele was of paternal origin. The mutation is more likely to have arisen from a de novo event but paternal gonadal mosaicism cannot be excluded. In conclusion, we report the clinical and molecular details of the first case of a de novo deletion of the downstream PAR1 region in an LWD individual. De novo deletions of SHOX and the downstream enhancer region must be therefore considered in cases of isolated LWD.

Variability in sex-determining mechanisms influences genome complexity in reptilia

In this review, we describe the history of amniote sex determination as a classic example of Darwinian evolution. We suggest that evolutionary changes in sex determination provide a foundation for understanding important aspects of chromosome and genome organization that otherwise appear haphazard in their origins and contents. Species with genotypic sex determination often possess heteromorphic sex chromosomes, whereas species with environmental sex determination lack them. Through a series of mutations followed by selection at key genes, sex-determining mechanisms have turned over many times throughout the amniote lineage. As a consequence, amniote genomes have undergone gains or losses of sex chromosomes. We review the genomic and ecological contexts in which either temperature-dependent or genotypic sex determination has evolved. Once genotypic sex determination emerges in a lineage, viviparity and heteromorphic sex chromosomes become more likely to evolve. For example, in extinct marine reptiles, genotypic sex determination apparently led to viviparity, which in turn facilitated their pelagic radiation. Sex chromosomes comprise genome regions that differ from autosomes in recombination rate, mutation rate, levels of polymorphism, and the presence of sex-determining and sexually antagonistic genes. In short, many aspects of amniote genome complexity, life history, and adaptive radiation appear contingent on evolutionary changes in sex-determining mechanisms.

Product length, dye choice, and detection chemistry in the bead-emulsion amplification of millions of single DNA molecules in parallel

The amplification of millions of single molecules in parallel can be performed on microscopic magnetic beads that are contained in aqueous compartments of an oil-buffer emulsion. These bead-emulsion amplification (BEA) reactions result in beads that are covered by almost-identical copies derived from a single template. The post-amplification analysis is performed using different fluorophore-labeled probes. We have identified BEA reaction conditions that efficiently produce longer amplicons of up to 450 base pairs. These conditions include the use of a Titanium Taq amplification system. Second, we explored alternate fluorophores coupled to probes for post-PCR DNA analysis. We demonstrate that four different Alexa fluorophores can be used simultaneously with extremely low crosstalk. Finally, we developed an allele-specific extension chemistry that is based on Alexa dyes to query individual nucleotides of the amplified material that is both highly efficient and specific.

Short stature homeoboxcontaining gene and idiopathic short stature

The term idiopathic short stature (ISS) refers to patients who are short due to various unknown reasons. Although it is clear that multiple factors contribute to final height, genetic factors play a crucial role. Mutations of a human homeobox gene, short stature homeobox-containing (SHOX) gene, have been shown to be associated with the short stature phenotype in patients with Turner syndrome, most patients with Leri-Weill dyschondrosteosis and some cases of ISS. The prevalence of SHOX anomalies in subjects previously recognized as having ISS has been estimated at 2.4% in a large series of ISS individuals. This review focuses on the functional properties of the SHOX gene and its linkage to ISS.

Mixed modeling and multiple imputation for unobservable genotype clusters

Understanding the genetic contributions to complex diseases will require consideration of interaction across multiple genes and environmental factors. At the same time, capturing information on allelic phase, that is, whether alleles within a gene are in cis (on the same chromosome) or in trans (on different chromosomes), is critical when using haplotypic approaches in disease association studies. This paper proposes a combination of mixed modeling and multiple imputation for assessing high-order genotype-phenotype associations while accounting for the uncertainty in phase inherent in population-based association studies. This method provides a flexible statistical framework for controlling for potential confounders and assessing gene-environment and gene-gene interactions in studies of unrelated individuals where the haplotypic phase is generally unobservable. The proposed method is applied to a cohort of 626 subjects with human immunodeficiency virus (HIV) to assess the potential contribution of four genes, apolipoprotein-C-III, apolipoprotein-E, endothelial lipase and hepatic lipase in predicting lipid abnormalities. A simulation study is also presented to describe the method performance.

Characterization of the bovine pseudoautosomal boundary: Documenting the evolutionary history of mammalian sex chromosomes

Here, we report the sequence characterization of the bovine pseudoautosomal boundary (PAB) and its neighborhood. We demonstrate that it maps to the 5' end of the GPR143 gene, which has concomitantly lost upstream noncoding exons on the Y chromosome. We show that the bovine PAB was created approximately 20.7 million years ago by illegitimate intrachromatid recombination between inverted, ruminant-specific Bov-tA repeats. Accordingly, we demonstrate that cattle share their PAB with all other examined ruminants including sheep, but not with cetaceans or more distantly related mammals. We provide evidence that, since its creation, the ancestral ruminant PAB has been displaced by attrition, which occurs at variable rates in different species, and that it is capable of retreat by attrition erasure. We have estimated the ratio of male to female mutation rates in the Bovidae family as approximately 1.7, and we provide evidence that the mutation rate is higher in the recombining pseudoautosomal region than in the adjacent, nonrecombining gonosome-specific sequences.

Whole population, genome-wide mapping of hidden relatedness

We present GERMLINE, a robust algorithm for identifying segmental sharing indicative of recent common ancestry between pairs of individuals. Unlike methods with comparable objectives, GERMLINE scales linearly with the number of samples, enabling analysis of whole-genome data in large cohorts. Our approach is based on a dictionary of haplotypes that is used to efficiently discover short exact matches between individuals. We then expand these matches using dynamic programming to identify long, nearly identical segmental sharing that is indicative of relatedness. We use GERMLINE to comprehensively survey hidden relatedness both in the HapMap as well as in a densely typed island population of 3000 individuals. We verify that GERMLINE is in concordance with other methods when they can process the data, and also facilitates analysis of larger scale studies. We bolster these results by demonstrating novel applications of precise analysis of hidden relatedness for (1) identification and resolution of phasing errors and (2) exposing polymorphic deletions that are otherwise challenging to detect. This finding is supported by concordance of detected deletions with other evidence from independent databases and statistical analyses of fluorescence intensity not used by GERMLINE.

The SWI/SNF protein ATRX co-regulates pseudoautosomal genes that have translocated to autosomes in the mouse genome

BACKGROUND

Pseudoautosomal regions (PAR1 and PAR2) in eutherians retain homologous regions between the X and Y chromosomes that play a critical role in the obligatory X-Y crossover during male meiosis. Genes that reside in the PAR1 are exceptional in that they are rich in repetitive sequences and undergo a very high rate of recombination. Remarkably, murine PAR1 homologs have translocated to various autosomes, reflecting the complex recombination history during the evolution of the mammalian X chromosome.

RESULTS

We now report that the SNF2-type chromatin remodeling protein ATRX controls the expression of eutherian ancestral PAR1 genes that have translocated to autosomes in the mouse. In addition, we have identified two potentially novel mouse PAR1 orthologs.

CONCLUSION

We propose that the ancestral PAR1 genes share a common epigenetic environment that allows ATRX to control their expression.

Recombination and nucleotide diversity in the sex chromosomal pseudoautosomal region of the emu, Dromaius novaehollandiae

Pseudoautosomal regions (PARs) shared by avian Z and W sex chromosomes are typically small homologous regions within which recombination still occurs and are hypothesized to share the properties of autosomes. We capitalized on the unusual structure of the sex chromosomes of emus, Dromaius novaehollandiae, which consist almost entirely of PAR shared by both sex chromosomes, to test this hypothesis. We compared recombination, linkage disequilibrium (LD), GC content, and nucleotide diversity between pseudoautosomal and autosomal loci derived from 11 emu bacterial artificial chromosome (BAC) clones that were mapped to chromosomes by fluorescent in situ hybridization. Nucleotide diversity (pi = 4N(e)mu) was not significantly lower in pseudoautosomal loci (14 loci, 1.9 +/- 2.4 x 10(-3)) than autosomal loci (8 loci, 4.2 +/- 6.1 x 10(-3)). By contrast, recombination per site within BAC-end sequences (rho = 4Nc) (pseudoautosomal, 3.9 +/- 6.9 x 10(-2); autosomal, 2.3 +/- 3.7 x 10(-2)) was higher and average LD (D') (pseudoautosomal, 4.2 +/- 0.2 x 10(-1); autosomal, 4.7 +/- 0.5 x 10(-1)) slightly lower in pseudoautosomal sequences. We also report evidence of deviation from a simple neutral model in the PAR and in autosomal loci, possibly caused by departures from demographic equilibrium, such as population growth. This study provides a snapshot of the population genetics of avian sex chromosomes at an early stage of differentiation.

Mammalian meiotic recombination hot spots

Our understanding of the details of mammalian meiotic recombination has recently advanced significantly. Sperm typing technologies, linkage studies, and computational inferences from population genetic data have together provided information in unprecedented detail about the location and activity of the sites of crossing-over in mice and humans. The results show that the vast majority of meiotic recombination events are localized to narrow DNA regions (hot spots) that constitute only a small fraction of the genome. The data also suggest that the molecular basis of hot spot activity is unlikely to be strictly determined by specific DNA sequence motifs in cis. Further molecular studies are needed to understand how hot spots originate, function and evolve.

Variation in patterns of human meiotic recombination

In the last 30 years it has become evident that patterns of meiotic recombination can be highly variable among individuals. The evidence comes from both low and high resolution analyses of hotspots of recombination in human and other species. In addition, a comparison of the recombination profiles in closely related species such as human and chimpanzee reveals essentially no correlation in the position of hotspots. Although the variation in hotspots of meiotic recombination is clearly documented, the mechanisms responsible for such variation are far from being understood. Here we will review the available evidence of natural variation in meiotic recombination and will discuss potential implications of this variation on the functional mechanisms of crossover formation and control.

Genotypes and phenotypes in children with short stature: clinical indicators of SHOX haploinsufficiency

BACKGROUND

Short stature affects approximately 2% of children, representing one of the more frequent disorders for which clinical attention is sought during childhood. Despite assumed genetic heterogeneity, mutations or deletions of the short stature homeobox-containing gene (SHOX) are found quite frequently in subjects with short stature. Haploinsufficiency of the SHOX gene causes short stature with highly variable clinical severity, ranging from isolated short stature without dysmorphic features to Léri-Weill syndrome, and with no functional copy of the SHOX gene, Langer syndrome.

METHODS

To characterise the clinical and molecular spectrum of SHOX deficiency in childhood we assessed the association between genotype and phenotype in a large cohort of children of short stature from 14 countries.

RESULTS

Screening of 1608 unrelated individuals with sporadic or familial short stature revealed SHOX mutations or deletions in 68 individuals (4.2%): complete deletions in 48 (70.6%), partial deletions in 4 (5.9%) and point mutations in 16 individuals (23.5%). Although mean height standard deviation score (SDS) was not different between participants of short stature with or without identified SHOX gene defects (-2.6 vs -2.6), detailed examination revealed that certain bone deformities and dysmorphic signs, such as short forearm and lower leg, cubitus valgus, Madelung deformity, high-arched palate and muscular hypertrophy, differed markedly between participants with or without SHOX gene defects (p<0.001). Phenotypic data were also compared for 33 children with Turner syndrome in whom haploinsufficiency of SHOX is thought to be responsible for the height deficit.

CONCLUSION

A phenotype scoring system was developed that could assist in identifying the most appropriate subjects for SHOX testing. This study offers a detailed genotype-phenotype analysis in a large cohort of children of short stature, and provides quantitative clinical guidelines for testing of the SHOX gene.

An evolutionary view of human recombination

Recombination has essential functions in mammalian meiosis, which impose several constraints on the recombination process. However, recent studies have shown that, in spite of these roles, recombination rates vary tremendously among humans, and show marked differences between humans and closely related species. These findings provide important insights into the determinants of recombination rates and raise new questions about the selective pressures that affect recombination over different genomic scales, with implications for human genetics and evolutionary biology.

Significant positive correlation between the recombination rate and GC content in the human pseudoautosomal region

This paper establishes that recombination drives the evolution of GC content in a significant way. Because the human P-arm pseudoautosomal region (PAR1) has been shown to have a high recombination rate, at least 20-fold more frequent than the genomic average of approximately 1 cM/Mb, this region provides an ideal system to study the role of recombination in the evolution of base composition. Nine non-coding regions of PAR1 are analyzed in this study. We have observed a highly significant positive correlation between the recombination rate and GC content (rho = 0.837, p < or = 0.005). Five regions that lie in the distal part of PAR1 are shown to be significantly higher than genomic average divergence. By comparing the intra- and inter-specific AT->GC -GC->AT ratios, we have detected no fixation bias toward GC alleles except for L254915, which has excessive AT-->GC changes in the human lineage. Thus, we conclude that the high GC content of the PAR1 genes better fits the biased gene conversion (BGC) model.

The distribution and causes of meiotic recombination in the human genome

Using the statistical analysis of genetic variation, we have developed a high-resolution genetic map of recombination hotspots and recombination rate variation across the human genome. This map, which has a resolution several orders of magnitude greater than previous studies, identifies over 25,000 recombination hotspots and gives new insights into the distribution and determination of recombination. Wavelet-based analysis demonstrates scale-specific influences of base composition, coding context and DNA repeats on recombination rates, though, in contrast with other species, no association with DNase I hypersensitivity. We have also identified specific DNA motifs that are strongly associated with recombination hotspots and whose activity is influenced by local context. Comparative analysis of recombination rates in humans and chimpanzees demonstrates very high rates of evolution of the fine-scale structure of the recombination landscape. In the light of these observations, we suggest possible resolutions of the hotspot paradox.

The pseudoautosomal regions, SHOX and disease

The pseudoautosomal regions represent blocks of sequence identity between the mammalian sex chromosomes. In humans, they reside at the ends of the X and Y chromosomes and encompass roughly 2.7 Mb (PAR1) and 0.33 Mb (PAR2). As a major asset of recently available sequence data, our view of their structural characteristics could be refined considerably. While PAR2 resembles the overall sequence composition of the X chromosome and exhibits only slightly elevated recombination rates, PAR1 is characterized by a significantly higher GC content and a completely different repeat structure. In addition, it exhibits one of the highest recombination frequencies throughout the entire human genome and, probably as a consequence of its structural features, displays a significantly faster rate of evolution. It therefore represents an exceptional model to explore the correlation between meiotic recombination and evolutionary forces such as gene mutation and conversion. At least twenty-nine genes lie within the human pseudoautosomal regions, and these genes exhibit 'autosomal' rather than sex-specific inheritance. All genes within PAR1 escape X inactivation and are therefore candidates for the etiology of haploinsufficiency disorders including Turner syndrome (45,X). However, the only known disease gene within the pseudoautosomal regions is the SHORT STATURE HOMEBOX (SHOX) gene, functional loss of which is causally related to various short stature conditions and disturbed bone development. Recent analyses have furthermore revealed that the phosphorylation-sensitive function of SHOX is directly involved in chondrocyte differentiation and maturation.

Molecular features of meiotic recombination hot spots

Meiotic recombination occurs preferentially at certain regions called hot spots and is important for generating genetic diversity and proper segregation of chromosomes during meiosis. Hot spots have been characterized most extensively in yeast, mice and humans. The development of methods based on sperm typing and population genetics has facilitated rapid and high-resolution mapping of hot spots in mice and humans in recent years. With increasing information becoming available on meiotic recombination in different species, it is now possible to compare several molecular features associated with hot-spot loci. Further, there have been advances in our knowledge of the factors influencing hot-spot activity and the role that they play in structuring the genome into haplotype blocks. We review the molecular features associated with hot spots in terms of their properties and mechanisms underlying their function and distribution. A large number of these features seem to be shared among hot spots from different species suggesting common mechanisms for their formation and function.

High-resolution recombination patterns in a region of human chromosome 21 measured by sperm typing

For decades, classical crossover studies and linkage disequilibrium (LD) analysis of genomic regions suggested that human meiotic crossovers may not be randomly distributed along chromosomes but are focused instead in “hot spots.” Recent sperm typing studies provided data at very high resolution and accuracy that defined the physical limits of a number of hot spots. The data were also used to test whether patterns of LD can predict hot spot locations. These sperm typing studies focused on several small regions of the genome already known or suspected of containing a hot spot based on the presence of LD breakdown or previous experimental evidence of hot spot activity. Comparable data on target regions not specifically chosen using these two criteria is lacking but is needed to make an unbiased test of whether LD data alone can accurately predict active hot spots. We used sperm typing to estimate recombination in 17 almost contiguous ~5 kb intervals spanning 103 kb of human Chromosome 21. We found two intervals that contained new hot spots. The comparison of our data with recombination rates predicted by statistical analyses of LD showed that, overall, the two datasets corresponded well, except for one predicted hot spot that showed little crossing over. This study doubles the experimental data on recombination in men at the highest resolution and accuracy and supports the emerging genome-wide picture that recombination is localized in small regions separated by cold areas. Detailed study of one of the new hot spots revealed a sperm donor with a decrease in recombination intensity at the canonical recombination site but an increase in crossover activity nearby. This unique finding suggests that the position and intensity of hot spots may evolve by means of a concerted mechanism that maintains the overall recombination intensity in the region.

Meiotic crossover events are not randomly distributed across the human genome, but are concentrated in many small regions of a few kb with high recombination rates compared to surrounding regions. How the distribution of recombination events affects the association of different alleles along the chromosome (linkage disequilibrium, or LD) was recently addressed using sperm typing in regions already known or suspected to contain unusually high recombination intensities. In the current paper, the authors used sperm typing to examine recombination in a region not known or suspected of containing recombination hot spots. They first established the crossover distribution pattern within a 103-kb region of human Chromosome 21. Then, they compared their data to predictions of crossover distributions estimated by statistical analyses of polymorphism in the region. They found a good concordance between the two, although it was not perfect. To the authors' knowledge, this work is the first to compare LD-based estimates of recombination to sperm-typing data from regions not previously known or suspected of containing recombination hot spots. In addition, one of the studied hot spots revealed an example of a decrease in recombination intensity with a concurrent increase at a nearby site. This unique observation suggests that the activity of hot spots may evolve in a concerted fashion such that the overall recombination activity of the region is maintained.

High incidence of SHOX anomalies in individuals with short stature

OBJECTIVE

To study the SHOX gene and the PAR1 region in individuals with short stature.

METHODS

The study involved 56 cases of dyschondrosteosis and 84 cases of idiopathic short stature (ISS). The study was designed to determine the following: the prevalence of SHOX anomalies in ISS; the frequency of Madelung deformity in individuals with SHOX anomalies; and the value of a family history of short stature in deciding whether to test for the SHOX gene.

RESULTS

54 SHOX anomalies were observed, including 42 (68%) in the dyschondrosteosis group and 12 (15%) in the ISS group. The high frequency of SHOX anomalies in the ISS group can be explained by the large proportion of boys in this group, reflecting the difficulty in diagnosing dyschondrosteosis in young boys. Clinical evidence of Madelung deformity in six parents of ISS individuals emphasised the importance of family evaluation. Among the 54 SHOX anomalies, 33 PAR1 deletions were identified encompassing the SHOX gene (62%), one partial intragenic deletion (2%), nine deletions located downstream of the SHOX gene (16%), and 11 point mutations (20%).

CONCLUSIONS

These data emphasise the value of using microsatellite markers located within and downstream of the SHOX gene.

Word frequency analysis reveals enrichment of dinucleotide repeats on the human X chromosome and [GATA]n in the X escape region

Most of the human genome encodes neither protein nor known functional RNA, yet available approaches to seek meaningful information in the "noncoding" sequence are limited. The unique biology of the X chromosome, one of which is silenced in mammalian females, can yield clues into sequence motifs involved in chromosome packaging and function. Although autosomal chromatin has some capacity for inactivation, evidence indicates that sequences enriched on the X chromosome render it fully competent for silencing, except in specific regions that escape inactivation. Here we have used a linguistic approach by analyzing the frequency and distribution of nine base-pair genomic "words" throughout the human genome. Results identify previously unknown sequence differences on the human X chromosome. Notably, the dinucleotide repeats [AT]n, [AC]n, and [AG]n are significantly enriched across the X chromosome compared with autosomes. Moreover, a striking enrichment (>10-fold) of [GATA]n is revealed throughout the 10-Mb segment at Xp22 that escapes inactivation, and is confirmed by fluorescence in situ hybridization. A similar enrichment is found in other eutherian genomes. Our findings clearly demonstrate sequence differences relevant to the novel biology and evolution of the X chromosome. Furthermore, they implicate simple sequence repeats, linked to gene regulation and unusual DNA structures, in the regulation and formation of facultative heterochromatin. Results suggest a new paradigm whereby a regional escape from X inactivation is due to the presence of elements that prevent heterochromatinization, rather than the lack of other elements that promote it.

An integrated genetic map for linkage analysis

Here I describe an Internet accessible database containing interpolated genetic map positions for 12917 marker loci. These are estimated via locally weighted linear regression (loess) from the Build 35.1 physical map position and the linkage map of Kong, X., and coworkers (2004) Am. J. Hum. Genet. 75:1143-1148. For the pseudoautosomal region, I have interpolated a male map based on the sperm typing data of Lien, S., and coworkers (2000) Am. J. Hum. Genet. 66:557-566.

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.

Human polymorphism and human-chimpanzee divergence in pseudoautosomal region correlate with local recombination rate

Previous studies have shown widespread correlation between nucleotide polymorphism and recombination rate, but the cause of this correlation is unresolved. One explanation is that recombination is associated with point mutations, potentially through mutagenic effects of meiotic crossover. This hypothesis predicts that regions of frequent recombination should show both elevated nucleotide diversity within a species and increased nucleotide divergence between species. Here we tested this hypothesis by studying the human short-arm pseudoautosomal region (PAR1), which recombines between X and Y chromosomes in men at a rate approximately 20 times the genome average. We sequenced dispersed intronic loci within PAR1 in a panel of humans and in the chimpanzee and directly measured sequence variation and recombination rate from these data. In line with previous reports, we saw a correlation between human polymorphism level and local recombination rate. Moreover, we also found a highly significant correlation between human-chimpanzee divergence and recombination rate. These results are consistent with the hypothesis that recombination is associated with point mutations, possibly because recombination is mutagenic.