Gene organisation, sequence variation and isochore structure at the centromeric boundary of the human MHC.

We have mapped and sequenced the region immediately centromeric of the human major histocompatibility complex (MHC). A cluster of 13 genes/pseudogenes was identified in a 175 kb PAC linking the TAPASIN locus with the class II region. It includes two novel human genes (BING4 and SACM2L) and a thus far unnoticed human leucocyte antigen (HLA) class II pseudogene, termed HLA-DPA3. Analysis of the G+C content revealed an isochore boundary which, together with the previously reported telomeric boundary, defines the MHC class II region as one of the first completely sequenced isochores in the human genome. Comparison of the sequence with limited sequence from other cell lines shows that the high sequence variation found within the classical class II region extends beyond the identified isochore boundary leading us to propose the concept of an "extended MHC". By comparative analysis, we have precisely identified the mouse/human synteny breakpoint at the centromeric end of the extended MHC class II region between the genes HSET and PHF1.     

“Beijing Region” (3pter-D3S3397) of the Human Genome: Complete sequence and analysis

The Human Genome Project (HGP), undoubtedly one of the greatest endeavors in the history of natural sciences, aims at determining the DNA sequence of the 3-billion-basepair human genome that holds an extraordinary trove of information about evolution, growth, and development of humans. China officially joined HGP and became one of the six member states of the International Human Genome Sequencing Con-sortium (IHGSC), alongside France, Germany, Japan, United Kingdom, and United State…     

Rhesus macaque class I duplicon structures, organization, and evolution within the alpha block of the major histocompatibility complex

The alpha block of the human and chimpanzee major histocompatibility complex (MHC) class I genomic region contains 10 to 11 duplicated MHC class I genes, including the HLA/Patr-A, -G, and -F genes. In comparison, the alpha block of the rhesus macaque (Macaca mulatta, Mamu) has an additional 20 MHC class I genes within this orthologous region. The present study describes the identification and analysis of the duplicated segmental genomic structures (duplicons) and genomic markers within the alpha block of the rhesus macaque and their use to reconstruct the duplication history of the genes within this region. A variety of MHC class I genes, pseudogenes, transposons, and retrotransposons, such as Alu and ERV16, were used to categorize the 28 duplicons into four distinct structural categories. The phylogenetic relationship of MHC class I genes, Alu, and LTR16B sequences within the duplicons was examined by use of the Neighbor-Joining (NJ) method. Two single-duplicon tandem duplications, two polyduplicon tandem duplications with an accompanying inversion product per duplication, eight polyduplicon tandem duplications steps, 12 deletions, and at least two recombinations were reconstructed to explain the highly complex organization and evolution of the 28 duplicons (nine inversions) within the Mamu alpha block. On the basis of the phylogenetic evidence and the reconstructed tandem duplication history of the 28 duplicons, the Mamu/Patr/HLA-F ortholog was the first MHC class I gene to have been fixed without further duplication within the alpha block of primates. Assuming that the rhesus macaque and the chimpanzee/human lineages had started with the same number of MHC class I duplicons at the time of their divergence approximately 24 to 31 MYA, then the number of genes within the alpha block have been duplicated at an approximately three times greater rate in the rhesus macaque than in either the human or chimpanzee.     

“Beijing Region” (3pter-D3S3397) of the Human Genome: Complete sequence and analysis

The goal of the Human Genome Project (HGP) is to determine a complete and high-quality sequence of the human genome. China, as one of the six member states, takes a region between 3pter and D3S3397 of the human chromosome 3 as its share of this historic project, referred as “Beijing Region”. The complete sequence of this region comprises of 17.4 megabasepairs (Mb) with an average GC content of 42% and an average recombination rate of 2.14 cM/Mb. Within Beijing Region, 122 known and 20 novel genes are identified, as well as 42607 single nucleotide polymorphisms (SNPs). Comprehensive analyses also reveal: (i) gene density and GC-content of Beijing Region are in agreement with human cytogenetic maps, i.e. G-minus bands are GC-rich and of a high gene density, whereas G-plus bands are GC-poor and of a relatively low gene density; (ii) the average recombination rate within Beijing Region is relatively high compared with other regions of chromosome 3, with the highest recombination rate of 6.06 cM/Mb in the subtelomeric area; (iii) it is most likely that a large gene, associated with the mammary gland, may reside in the 1.1 Mb gene-poor area near the telomere; (iv) many disease-related genes are genetically mapped to Beijing Region, including those associated with cancers and metabolic syndromes. All make Beijing Region an important target for in-depth molecular investigations with a purpose of medical applications.     

Physical mapping of the barley stem rust resistance gene rpg4

The barley stem rust resistance gene rpg4 was physically and genetically localized on two overlapping BAC clones covering an estimated 300-kb region of the long arm of barley chromosome 7(5H). Initially, our target was mapped within a 6.0-cM region between the previously described flanking markers MWG740 and ABG391. This region was then saturated by integrating new markers from several existing barley and rice maps and by using BAC libraries of barley cv. Morex and rice cv. Nipponbare. Physical/genetic distances in the vicinity of rpg4 were found to be 1.0 Mb/cM, which is lower than the average for barley (4 Mb/cM) and lower than that determined by translocation breakpoint mapping (1.8 Mb/cM). Synteny at high resolution levels has been established between the region of barley chromosome 7(5H) containing the rpg4 locus and the subtelomeric region of rice chromosome 3 between markers S16474 and E10757. This 1.7-cM segment of the rice genome was covered by two overlapping BAC clones, about 250 kb of total length. In barley the markers S16474 and E10757 genetically delimit rpg4, lying 0.6 cM distal and 0.4 cM proximal to the locus, respectively.     

Evidence for human meiotic recombination interference obtained through construction of a short tandem repeat-polymorphism linkage map of chromosome 19

An improved linkage map for human chromosome 19 containing 35 short tandem repeat polymorphisms (STRPs) and one VNTR (D19S20) was constructed. The map included 12 new (GATA)_n tetranucleotide STRPs. Although total lengths of the male (114 cM) and female (128 cM) maps were similar, at both ends of the chromosome male recombination exceeded female recombination, while in the interior portion of the map female recombination was in excess. Cosmid clones containing the STRP sequences were identified and were positioned along the chromosome by fluorescent in situ hybridization. Four rounds of careful checking and removal of genotyping errors allowed biologically relevant conclusions to be made concerning the numbers and distributions of recombination events on chromosome 19. The average numbers of recombinations per chromosome matched closely the lengths of the genetic maps computed by using the program CRIMAP. Significant numbers of chromosomes with zero, one, two, or three recombinations were detected as products of both female and male meioses. On the basis of the total number of observed pairs of recombination events in which only a single informative marker was situated between the two recombinations, a maximal estimate for the rate of meiotic STRP “gene” conversion without recombination was calculated as 3 × 10⁻⁴/meiosis. For distances up to 30 cM between recombinations, many fewer chromosomes which had undergone exactly two recombinations were observed than were expected on the basis of the assumption of independent recombination locations. This strong new evidence for human meiotic interference will help to improve the accuracy of interpretation of clinical DNA test results involving polymorphisms flanking a genetic abnormality.     

Recombination and mutation shape variations in the major histocompatibility complex

The major histocompatibility complex (MHC) is closely associated with numerous diseases, but its high degree of polymorphism complicates the discovery of disease-associated variants. In principle, recombination and de novo mutations are two critical factors responsible for MHC polymorphisms. However, direct evidence for this hypothesis is lacking. Here, we report the generation of fine-scale MHC recombination and de novo mutation maps of ～5 Mb by deep sequencing (> 100×) of the MHC genome for 17 MHC recombination and 30 non-recombination Han Chinese families (a total of 190 individuals). Recombination hotspots and Han-specific breakpoints are located in close proximity at haplotype block boundaries. The average MHC de novo mutation rate is higher than the genome-wide de novo mutation rate, particularly in MHC recombinant individuals. Notably, mutation and recombination generated polymorphisms are located within and outside linkage disequilibrium regions of the MHC, respectively, and evolution of the MHC locus was mainly controlled by positive selection. These findings provide insights on the evolutionary causes of the MHC diversity and may facilitate the identification of disease-associated genetic variants.     

How strong is the mutagenicity of recombination in mammals?

It is commonly believed that a high recombination rate such as that in a pseudoautosomal region (PAR) greatly increases the mutation rate because a 170-fold increase was estimated for the mouse PAR region. However, sequencing PAR and non-PAR introns of the Fxy gene in four Mus taxa, we found an increase of only twofold to fivefold. Furthermore, analyses of sequence data from human and orangutan PAR and X-linked regions and from autosomal regions showed a weak effect of recombination on mutation rate (a slope of less than 0.2% per cM/Mb), although a much stronger effect on GC content (1% to 2% per cM/Mb). Because typical recombination rates in mammals are much lower than those in PARs, the mutagenicity of recombination is weak or, at best, moderate, although its effect on GC% is much stronger. In addition, contrary to a previous study, we found no Fxy duplicate in Mus spretus.     

An integrated map of human 6q22.3-q24 including a 3-Mb high-resolution BAC/PAC contig encompassing a QTL for fetal hemoglobin

Genetic studies have previously assigned a quantitative trait locus (QTL) for hemoglobin F and F cells to a region of approximately 4 Mb between the markers D6S408 and D6S292 on chromosome 6q23. An initial yeast artificial chromosome contig of 13 clones spanning this region was generated. Further linkage analysis of an extended kindred refined the candidate interval to 1-2 cM, and key recombination events now place the QTL within a region of <800 kb. We describe a high-resolution bacterial clone contig spanning 3 Mb covering this critical region. The map consists of 223 bacterial artificial chromosome (BAC) and 100 P1 artificial chromosome (PAC) clones ordered by sequence-tagged site (STS) content and restriction fragment fingerprinting with a minimum tiling path of 22 BACs and 1 PAC. A total of 194 STSs map to this interval of 3 Mb, giving an average marker resolution of approximately one per 15 kb. About half of the markers were novel and were isolated in the present study, including three CA repeats and 13 single nucleotide polymorphisms. Altogether 24 expressed sequence tags, 6 of which are unique genes, have been mapped to the contig.     

Recombination rates across porcine autosomes inferred from high-density linkage maps

Studies of the variation in recombination rate across the genome provide a better understanding of evolutionary genomics and are also an important step towards mapping and dissecting complex traits in domestic animals. With the recent completion of the porcine genome sequence and the availability of a high‐density porcine single nucleotide polymorphism (SNP) array, it is now possible to construct a high‐density porcine linkage map and estimate recombination rate across the genome. A total of 416 animals were genotyped with the Porcine SNP60BeadChip, and high‐density chromosome linkage maps were constructed using CRI‐MAP, assuming the physical order of the Sscrofa10 assembly. The total linkage map length was 2018.79 cM, using 658 meioses and 14 503 SNPs. The estimated average recombination rate across the porcine autosomes was 0.86 cM/Mb. However, a large variation in recombination rate was observed among chromosomes. The estimated average recombination rates (cM/Mb) per chromosome ranged from 0.48 in SSC1 to 1.48 in SSC10, displaying a significant negative correlation with the chromosome sizes. In addition, the analysis of the variation in the recombination rates taking 1‐Mb sliding windows has allowed us to demonstrate the variation in recombination rates within chromosomes. In general, a larger recombination rate was observed in the extremes than in the centre of the chromosome. Finally, the ratio between female and male recombination rates was also inferred, obtaining a value of 1.38, with the heterogametic sex having the least recombination.     

Genetic dissection of the human leukocyte antigen region by use of haplotypes of Tasmanians with multiple sclerosis

Association of multiple sclerosis (MS) with the human leukocyte antigen (HLA) class II haplotype DRB1*1501-DQB1*0602 is the most consistently replicated finding of genetic studies of the disease. However, the high level of linkage disequilibrium (LD) in the HLA region has hindered the identification of other loci that single-marker tests for association are unlikely to resolve. In order to address this issue, we generated haplotypes spanning 14.754 Mb (5 cM) across the entire HLA region. The haplotypes, which were inferred by genotyping relatives of 152 patients with MS and 105 unaffected control subjects of Tasmanian ancestry, define a genomic segment from D6S276 to D6S291, including 13 microsatellite markers integrated with allele-typing data for DRB1 and DQB1. Association to the DRB1*1501-DQB1*0602 haplotype was replicated. In addition, we found that the class I/extended class I region, defined by a genomic segment of approximately 400 kb between MOGCA and D6S265, harbors genes that independently increase risk of, or provide protection from, MS. Log-linear modeling analysis of constituent haplotypes that represent genomic regions containing class I (MOGCA-D6S265), class III (TNFa-TNFd-D6S273), and class II (DRB1-DQB1) genes indicated that having class I and class II susceptibility variants on the same haplotype provides an additive effect on risk. Moreover, we found no evidence for a disease locus in the class III region defined by a 150-kb genomic segment containing the TNF locus and 14 other genes. A global overview of LD performed using GOLD identified two discrete blocks of LD in the HLA region that correspond well with previous findings. We propose that the analysis of haplotypes, by use of the types of approaches outlined in the present article, should make it possible to more accurately define the contribution of the HLA to MS.     

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.     

Genetic structure of a QTL hotspot on chromosome 2 in sweet cherry indicates positive selection for favorable haplotypes

A genomic region of particular interest for sweet cherry (Prunus avium L.) breeding is a quantitative trait locus (QTL) “hotspot” on chromosome 2. QTLs for fruit size, firmness, sweetness, and flowering time are reported to map to this region. An understanding of genetic diversity, allele sources, linkage relationships, and historical recombinations is critical to enable the combining of favorable alleles at multiple loci. The objectives of this study were to characterize, visualize, and interpret the genetic structure of this previously identified QTL hotspot within North American sweet cherry breeding germplasm, using a pedigree-based haploblocking approach. Across the 29.4 cM (6.3 Mbp) region defined by single nucleotide polymorphism (SNP) information from the RosBREED cherry 6K SNP array v1, a total of 12 recombination events falling into six inter-marker regions were traced within the pedigree of elite and wild germplasm (n = 55). These recombinations defined five haploblocks containing 5–15 markers and exhibiting 7–11 haplotypes each. Over the entire QTL hotspot, 30 extended haplotypes were identified for which parental gametes could be determined. When the haploblocks and their haplotypes were used to explore genetic diversity, ancestry, and recombination patterns, and then integrated with previous QTL results for fruit size, the results indicated that favorable alleles at this QTL hotspot are under positive selection in breeding. The genetic framework provided by a haploblock approach and knowledge of haplotype-level diversity sets the stage for assigning breeding utility to these haplotypes.     

An Integrated Tool to Study MHC Region: Accurate SNV Detection and HLA Genes Typing in Human MHC Region Using Targeted High-Throughput Sequencing

The major histocompatibility complex (MHC) is one of the most variable and gene-dense regions of the human genome. Most studies of the MHC, and associated regions, focus on minor variants and HLA typing, many of which have been demonstrated to be associated with human disease susceptibility and metabolic pathways. However, the detection of variants in the MHC region, and diagnostic HLA typing, still lacks a coherent, standardized, cost effective and high coverage protocol of clinical quality and reliability. In this paper, we presented such a method for the accurate detection of minor variants and HLA types in the human MHC region, using high-throughput, high-coverage sequencing of target regions. A probe set was designed to template upon the 8 annotated human MHC haplotypes, and to encompass the 5 megabases (Mb) of the extended MHC region. We deployed our probes upon three, genetically diverse human samples for probe set evaluation, and sequencing data show that ∼97% of the MHC region, and over 99% of the genes in MHC region, are covered with sufficient depth and good evenness. 98% of genotypes called by this capture sequencing prove consistent with established HapMap genotypes. We have concurrently developed a one-step pipeline for calling any HLA type referenced in the IMGT/HLA database from this target capture sequencing data, which shows over 96% typing accuracy when deployed at 4 digital resolution. This cost-effective and highly accurate approach for variant detection and HLA typing in the MHC region may lend further insight into immune-mediated diseases studies, and may find clinical utility in transplantation medicine research. This one-step pipeline is released for general evaluation and use by the scientific community.     

Analysis of intragenic recombination at wx in rice: correlation between the molecular and genetic maps within the locus.

In plant genomes as well as other eukaryotic genomes, meiotic recombination does not occur uniformly. At the level of the gene, high recombination frequencies are often observed within genetic loci in maize, but this feature of intragenic recombination is not seen at the csr1 locus in Arabidopsis. These observations suggest that meiotic recombination in plant genomes varies considerably among species. In the present study we investigated meiotic recombination at the wx locus in rice. The mutation sites of wx mutants induced by ethyl methanesulfonate (EMS) treatment or gamma-ray irradiation and a spontaneous wx mutant were physically characterized, and the genetic distances between those wx mutation sites were estimated by pollen analysis. Based on these results, the recombination frequency at the wx locus in rice was estimated as 27.3 kb/cM, which was about 10 times higher than the average for the genome, suggesting that there was a radically different rate of meiotic recombination for intra- and intergenic regions in the rice genome.     

Phytophthora sojae avirulence genes Avr4 and Avr6 are located in a 24kb, recombination-rich region of genomic DNA

A cross between two different races (race 7xrace 25) of the soybean root and stem rot pathogen Phytophthora sojae was analyzed to characterize the genomic region flanking two cosegregating avirulence genes, Avr4 and Avr6. Both genes cosegregated in the ratio of 82:17 (avirulent:virulent) in an F(2) population, suggestive of a single locus controlling both phenotypes. A chromosome walk was commenced from RAPD marker OPE7.1C, 2.0cM distant from the Avr4/6 locus. Three overlapping cosmids were isolated which included genetic markers that flank the Avr4/6 locus. The chromosome walk spanned a physical distance of 67kb which represented a genetic map distance of 22.3cM, an average recombination frequency of 3.0kb/cM and 11.7-fold greater than the predicted average recombination frequency of 35.3kb/cM for the entire P. sojae genome. Six genes (cDNA clones) expressed from the Avr4/6 genomic region encompassed by the cosmid contig were identified. Single nucleotide polymorphisms and restriction fragment length polymorphisms showed these six genes were closely linked to the Avr4/6 locus. Physical mapping of the cDNA clones within the cosmid contig made it possible to deduce the precise linkage order of the cDNAs. None of the six cDNA clones appear to be candidates for Avr4/6. We conclude that two of these cDNA clones flank a physical region of approximately 24kb and 4.3cM that appears to include the Avr4/6 locus.     

The 1.5-Mb region spanning the myotonic dystrophy locus shows uniform recombination frequency.

The myotonic dystrophy (DM) mutation has been identified as a heritable unstable CTG trinucleotide repeat sequence. The intergenerational amplification of this sequence is an example of a new class of dynamic mutations responsible for human genetic diseases. To ascertain whether recombination activity influences, or is affected by, the presence of this unique sequence, a comprehensive study of the physical and genetic mapping data for the 1.5-Mb region of human chromosome 19q13.3, which contains the DM locus, was conducted. The recombination rate for this region was examined by correlating genetic distance to physical distance for six selected marker loci. The following markers span the DM region: 19qCEN-p alpha 1.4 (D19S37)-APOC2-CKM-pE0.8 (D19S115)-pGB2.6 (DM)-p134c (D19S51)-19qTER. Initial linear regression analysis of these two parameters failed to reveal a significant linear correlation (coefficient of determination, r2 = .19), suggesting nonuniform rates of recombination. However, the presence of a recombination hot spot was believed to be unlikely, as the marker-to-marker pairs that showed the greatest deviation in recombination frequency were not restricted to a specific region of the 1.5 Mb studied and had relatively broad confidence intervals, as reflected by low LOD values. A second linear regression analysis using only marker intervals with high LOD scores (Zmax > 22) showed linear correlation (r2 = .68) for the entire 1.5-Mb region. This analysis indicated a relatively uniform recombination frequency in the 1.5-Mb region spanning the DM locus. Furthermore, the recombinations observed were neither under- nor overrepresented on DM chromosomes. Consequently, recombination activity is unlikely to influence, or be affected by, the presence of the DM mutation.     

An approach to mapping haplotype-specific recombination sites in human MHC class III

Studies of the major histocompatibility complex (MHC) in mouse indicate that the recombination sites are not randomly distributed and their occurrence is haplotype-dependent. No data concerning haplotype-specific recombination sites in human are available due to the low number of informative families. To investigate haplotype-specific recombination sites in human MHC, we here describe an approach based on identification of recombinant haplotypes derived from one conserved haplotype at the population level. The recombination sites were mapped by comparing polymorphic markers between the recombinant and assumed original haplotypes. We tested this approach on the extended haplotype HLA A3; B47; Bf ^* F; C4A ^* 1; C4B ^* Q0; DR7, which is most suitable for this analysis. First, it carries a number of rare markers, and second, the haplotype, albeit rare in the general population, is frequent in patients with 21-hydroxylase (21OH) defect. We observed recombinants derived from this haplotype in patients with 21OH defect. All these haplotypes had the centromeric part (from Bf to DR) identical to the original haplotype, but they differed in HLA A and B. We therefore assumed that they underwent recombinations in the segment that separates the Bf and HLA B genes. Polymorphic markers indicated that all break points mapped to two segments near the TNF locus. This approach makes possible the mapping of preferential recombination sites in different haplotypes.     

Genetic mapping of RP1 on 8q11-q21 in an Australian family with autosomal dominant retinitis pigmentosa reduces the critical region to 4 cM between D8S601 and D8S285

The locus (RP1) for one form of autosomal dominant retinitis pigmentosa (adRP) was mapped on chromosome 8q11-q22 between D8S589 and D8S285, which are about 8 cM apart, by linkage analysis in an extended family ascertained in the USA. We have studied a multigeneration Australian family with adRP and found close linkage without recombination between the disease locus and D8S591, D8S566, and D8S166 (Z_max = 1.137– 4.650 at θ = 0.00), all mapped in the region known to harbor RP1. Assuming that the mutation of the same gene is responsible for the disease in both families, the analysis of multiply informative meioses in the American and Australian families places the adRP locus between D8S601 and D8S285, which reduces the critical region to about 4 cM, corresponding to approximately 4 Mb, which is completely covered by a yeast artificial chromosome contig assembled recently. Received: 23 April 1996 / Accepted: 3 July 1996     

A high‐density linkage map enables a second‐generation collared flycatcher genome assembly and reveals the patterns of avian recombination rate variation and chromosomal evolution

Detailed linkage and recombination rate maps are necessary to use the full potential of genome sequencing and population genomic analyses. We used a custom collared flycatcher 50 K SNP array to develop a high‐density linkage map with 37 262 markers assigned to 34 linkage groups in 33 autosomes and the Z chromosome. The best‐order map contained 4215 markers, with a total distance of 3132 cM and a mean genetic distance between markers of 0.12 cM . Facilitated by the array being designed to include markers from most scaffolds, we obtained a second‐generation assembly of the flycatcher genome that approaches full chromosome sequences (N50 super‐scaffold size 20.2 Mb and with 1.042 Gb (of 1.116 Gb) anchored to and mostly ordered and oriented along chromosomes). We found that flycatcher and zebra finch chromosomes are entirely syntenic but that inversions at mean rates of 1.5–2.0 event (6.6–7.5 Mb) per My have changed the organization within chromosomes, rates high enough for inversions to potentially have been involved with many speciation events during avian evolution. The mean recombination rate was 3.1 cM /Mb and correlated closely with chromosome size, from 2 cM /Mb for chromosomes >100 Mb to >10 cM /Mb for chromosomes <10 Mb. This size dependence seemed entirely due to an obligate recombination event per chromosome; if 50 cM was subtracted from the genetic lengths of chromosomes, the rate per physical unit DNA was constant across chromosomes. Flycatcher recombination rate showed similar variation along chromosomes as chicken but lacked the large interior recombination deserts characteristic of zebra finch chromosomes.