The cobia Rachycentron canadum, mainly distributed in the warm waters of tropical and subtropical regions around the world, remains a fish of considerable economic importance. Detailed diversity and the number of microsatellite sequences in the cobia genome are still unintelligible. The primary aim of this work was to identify and quantify the miscellaneous SSR sequences in the cobia genome. More than 280,000 sequences were sequenced and screened using next-generation sequencing technology and microsatellite identification. Perfect mononucleotide repeats, dinucleotide microsatellites, and trinucleotide microsatellites contain (A)10 /(T)10 , (AC)6 /(TG)6 , and (AAT)5-32 as the largest number of motifs in each type of microsatellite, respectively. The tetranucleotide and pentanucleotide microsatellites (TTM and PTM) consist of the largest number of motifs of both (ATCT)5-32 and (TCAT)5-31 in TTMs, and (CTCTC)5-9 in PTMs, whereas the hexanucleotide microsatellites are rarely observed in the cobia genome. All c. 38000 sequences of composite microsatellites are extremely diverse, including compound (11.71%), interrupted compound (71.77%), complex (0.45%), and interrupted complex (16.07%). In this study, we developed a convenient and useful recording system for writing down and categorizing diverse composite microsatellite types. This system will provide great support for exploring repeat origins, evolutionary mechanisms, and the application of polymorphic microsatellites.

Mutations in mismatch repair genes leading to mismatch repair (MMR) deficiency (dMMR) and microsatellite instability (MSI) have been implicated in multiple types of gynecologic malignancies. Endometrial carcinoma represents the largest group, with approximately 30% of these cancers caused by dMMR/MSI. Thus, testing for dMMR is now routine for endometrial cancer. Somatic mutations leading to dMMR account for approximately 90% of these cancers. However, in 5-10% of cases, MMR protein deficiency is due to a germline mutation in the mismatch repair genes MLH1, MSH2, MSH6, PMS2, or EPCAM. These germline mutations, known as Lynch syndrome, are associated with an increased risk of both endometrial and ovarian cancer, in addition to colorectal, gastric, urinary tract, and brain malignancies. So far, gynecological cancers with dMMR/MSI are not well characterized and markers for detection of MSI in gynecological cancers are not well defined. In addition, currently advanced endometrial cancers have a poor prognosis and are treated without regard to MSI status. Elucidation of the mechanism causing dMMR/MSI gynecological cancers would aid in diagnosis and therapeutic intervention. Recently, a new immunotherapy was approved for the treatment of solid tumors with MSI that have recurred or progressed after failing traditional treatment strategies. In this review, we summarize the MMR defects and MSI observed in gynecological cancers, their prognostic value, and advances in therapeutic strategies to treat these cancers.

The human leukocyte antigen (HLA) complex locus has shaped a framework for evolutionary processes because of the dense clustering and strong linkage disequilibrium (LD) of polymorphic genes. Although the landscape of LD among conventional single-nucleotide polymorphisms (SNPs) has been described, the data on the lineage of major histocompatibility complex (MHC) haplotype are limited to pairwise comparisons of several haplotypes in Caucasoid populations. Multi-allelic markers, including microsatellite markers, may provide us with a larger power to analyze the MHC haplotype lineage because the mutation rate of microsatellite exceeds that of SNPs by several orders of magnitude. In this study, we investigated the complex structure of repeat motifs in a microsatellite to figure out the structural lineage of HLA-Cw/-B segments in Japanese. It was found that the genetic differences of HLA-Cw/-B haplotype lineage were reflected by repeat motif patterns at C1_2_5 locus, suggesting that unique mutational dynamics of microsatellites may be a useful marker to chase the haplotype lineage.

High intra-population genetic diversity and multiple measures of genetic variability at STR loci are useful in inferring past evolutionary history. However, STRs, categorized by their repeat motif size, differ in a number of aspects, requiring separate analyses. We analyzed 783 STRs in 36 worldwide populations to examine marker suitability as well as correlations between various measurements, to evaluate the extent of genomic diversity present in modern human populations. The loci were grouped by type and analyzed separately for each population group. Genetic variation defined by gene diversity and allele size variance, shows different trends of variation across four types of STRs. Additionally, there is little variation of genetic diversity, but there is decreased allelic size variance with increasing repeat motifs. A poor correlation between genetic diversity and allelic size variance across loci in all groups for Di-STRs is probably caused by the presence of allelic size gaps. In contrast, allelic size variance, genetic diversity, and number of alleles are strongly correlated with both tri- and tetra-STRs. The positive correlation of allelic size variance and presence of gaps within the range of allelic sizes in Di-STRs alone explains these observations. An unexpected high imbalance index (beta) at Di-STRs due to high allelic size variance also supports this assertion.

Although microsatellite mutation rates generally increase with increasing length of the repeat tract, interruptions in a microsatellite may stabilize it. We have performed a direct analysis of the effect of microsatellite interruptions on mutation rate and spectrum in cultured mammalian cells. Two mononucleotide sequences (G(17) and A(17)) and a dinucleotide [(CA)(17)] were compared with interrupted repeats of the same size and with sequences of 8 repeat units. MMR-deficient (MMR(-)) cells were used for these studies to eliminate effects of this repair process. Mutation rates were determined by fluctuation analysis on cells containing a microsatellite sequence at the 5' end of an antibiotic-resistance gene; the vector carrying this sequence was integrated in the genome of the cells. In general, interrupted sequences had lower mutation rates than perfect ones of the same size, but the magnitude of the difference was dependent upon the sequence of the interrupting base(s). Some interrupted repeats had mutation rates that were lower than those of perfect sequences of the same length but similar to those of half the length. This suggests that interrupting bases effectively divide microsatellites into smaller repeat runs with mutational characteristics different from those of the corresponding full-length microsatellite. We conclude that interruptions decrease microsatellite mutation rate and influence the spectrum of frameshift mutations. The sequence of the interrupting base(s) determines the magnitude of the effect on mutation rate.

Mammalian genome contains a wide variety of repetitive DNA sequences of relatively unknown function. We report a novel 227 bp simple repeat DNA (3.3 DNA) with a d{(GA)7A(AG)7} dinucleotide mirror repeat from the rat (Rattus norvegicus) genome. 3.3 DNA showed 75-85% homology with several eukaryotic mRNAs due to (GA/CU)n dinucleotide repeats by nBlast search and a dispersed distribution in the rat genome by Southern blot hybridization with [32P]3.3 DNA. The d{(GA)7A(AG)7} mirror repeat formed a triplex (H-DNA)-like structure in vitro. Two large RNAs of 9.1 and 7.5 kb were detected by [32P]3.3 DNA in rat brain by Northern blot hybridization indicating expression of such simple sequence repeats at RNA level in vivo. Further, several cDNAs were isolated from a rat cDNA library by [32P]3.3 DNA probe. Three such cDNAs showed tissue-specific RNA expression in rat. pRT 4.1 cDNA showed strong expression of a 2.39 kb RNA in brain and spleen, pRT 5.5 cDNA showed strong expression of a 2.8 kb RNA in brain and a 3.9 kb RNA in lungs, and pRT 11.4 cDNA showed weak expression of a 2.4 kb RNA in lungs. Thus, genomic simple sequence repeats containing d(GA/CT)n dinucleotides are transcriptionally expressed and regulated in rat tissues. Such d(GA/CT)n dinucleotide repeats may form structural elements (e.g., triplex) which may be sites for functional regulation of genomic coding sequences as well as RNAs. This may be a general function of such transcriptionally active simple sequence repeats widely dispersed in mammalian genome.

The patterns of mutation and evolution at 13 microsatellite loci were studied in the filamentous fungal genus Neurospora. First, a detailed investigation was performed on five microsatellite loci by sequencing each microsatellite, together with its nonrepetitive flanking regions, from a set of 147 individuals from eight species of Neurospora. To elucidate the genealogical relationships among microsatellite alleles, repeat number was mapped onto trees constructed from flanking-sequence data. This approach allowed the potentially convergent microsatellite mutations to be placed in the evolutionary context of the less rapidly evolving flanking regions, revealing the complexities of the mutational processes that have generated the allelic diversity conventionally assessed in population genetic studies. In addition to changes in repeat number, frequent substitution mutations within the microsatellites were detected, as were substitutions and insertion/deletions within the flanking regions. By comparing microsatellite and flanking-sequence divergence, clear evidence of interspecific allele length homoplasy and microsatellite mutational saturation was observed, suggesting that these loci are not appropriate for inferring phylogenetic relationships among species. In contrast, little evidence of intraspecific mutational saturation was observed, confirming the utility of these loci for population-level analyses. Frequency distributions of alleles within species were generally consistent with the stepwise mutational model. By comparing variation within species at the microsatellites and the flanking-sequence, estimated microsatellite mutation rates were approximately 2500 times greater than mutation rates of flanking DNA and were consistent with estimates from yeast and fruit flies. A positive relationship between repeat number and variance in repeat number was significant across three genealogical depths, suggesting that longer microsatellite alleles are more mutable than shorter alleles. To test if the observed patterns of microsatellite variation and mutation could be generalized, an additional eight microsatellite loci were characterized and sequenced from a subset of the same Neurospora individuals.

Precise estimates of mutation rates at Y-chromosomal microsatellite STR (short tandem repeat) loci make an important basis for paternity diagnostics and dating of Y chromosome lineage origins. There are indications of considerable locus mutation rate variability between (inter-) and within (intra-) loci. We have studied nine Y-STR loci-DYS19, DYS389I/II, DYS390, DYS391, DYS392, DYS393, DYS385, and DYS388-in 1,766 father-son pairs of confirmed paternity (a total of 15,894 meioses). Five biallelic markers were also analyzed in the fathers-Tat, YAP, 12f2, SRY1532, and 92R7-defining haplogroups 1, 2, 3, 4, 9, and 16, respectively. A total of 36 fragment length mutations were observed: 24 gains (22 single-step, two double-step) and 12 single-step losses. Thus, there was a significant surplus of gains (p=0.045). Overall, the mutation rate was positively correlated to STR repeat length and there was a significant relative excess of losses in long alleles and gains in short alleles (p=0.043). In contrast to the situation in autosomal STR loci and in MSY-1, no noteworthy correlation between mutation rate and the father's age at the child's birth was observed. We observed significant interlocus differences in Y-STR mutation rates (p<0.01). The number of observed mutations ranged from zero in DYS392 to eight in DYS391 and DYS390. We have also demonstrated obvious differences in mutation rates between the haplogroups studied (p=0.024), a phenomenon that is a reflection of the dependence of mutation rate on allele size. Our study has thus demonstrated the necessity of not only locus-specific, but even allele-specific, mutation rate estimates for forensic and population genetic purposes, and provides a considerable basis for such estimates.

It has been argued that the level of genetic diversity in the modern durum wheat ( Triticum turgidum L. var. durum) elite germplasm may have declined due to the high selection pressure applied in breeding programs. In this study, 58 accessions covering a wide spectrum of genetic diversity of the cultivated durum wheat gene pool were characterized with 70 microsatellite loci (or simple sequence repeats, SSRs). On average, SSRs detected 5.6 different allelic variants per locus, with a mean diversity index (DI) equal to 0.56, thus revealing a diversity content comparable to those previously observed with SSRs in other small-grain cereal gene pools. The mean genetic similarity value was equal to 0.44. A highly diagnostic SSR set has been identified. A high variation in allele size was detected among SSR loci, suggesting a different suitability of these loci for estimating genetic diversity. The B genome was characterized by an overall polymorphism significantly higher than that of the A genome. Genetic diversity is organised in well-distinct sub-groups identified by the corresponding foundation-genotypes. A large portion (92.7%) of the molecular variation detected within the group of 45 modern cvs was accounted for by SSR alleles tracing back to ten foundation-genotypes; among those, the most recent CIMMYT-derived founders were genetically distant from the old Mediterranean ones. On the other hand, rare alleles were abundant, suggesting that a large number of genetic introgressions contributed to the foundation of the well-diversified germplasm herein considered. The profiles of recently released varieties indicate that the level of genetic diversity present in the modern durum wheat germplasm has actually increased over time.

Investigators have intensively evaluated the major histocompatibility (MHC) complex for sarcoidosis susceptibility genes with the majority of reports implicating the human leukocyte antigen (HLA)-DRB1 gene. Because most studies have been performed in white and Asian populations, we sought to determine which MHC genes might be risk factors for sarcoidosis in African Americans. We genotyped six microsatellite markers spanning 11.6 megabases that overlapped the MHC region on chromosome 6p21-22 in 225 nuclear families ascertained by African American probands with a history of sarcoidosis. Using a family-based association methods approach, we performed multiallelic tests of association between each marker and sarcoidosis. A statistically significant association was detected between sarcoidosis and the DQCAR marker (p = 0.002) less than two kilobases telomeric from the HLA-DQB1 gene. Typing two additional markers in this region revealed that DQCAR-G51152 haplotypes, spanning a 38-kilobase region across the HLA-DQB1 gene, were associated with sarcoidosis on a global level (p = 0.022). Analysis of individual DQCAR and G51152 alleles showed that the DQCAR 178 (expected = 21.0; observed = 10; p = 0.0005) and G51152 217 (expected = 25.6; observed = 14; p = 0.0009) alleles were transmitted to affected offspring less often than expected; whereas the DQCAR 182 allele was transmitted more often than expected (expected = 52.6; observed = 66; p = 0.002). Our results indicate that HLA-DQB1 and not HLA-DRB1 plays an important role in sarcoidosis susceptibility in African Americans. Identification of the specific HLA-DQB1 alleles that influence sarcoidosis susceptibility in African Americans and the study of their antigenic-binding properties may reveal why African Americans suffer disproportionately from this disease.

Each independently evolving segment of the genomes of a sexually reproducing organism has a separate history reflecting part of the evolutionary history of that organism. Uniparentally or clonally inherited DNA segments such as the mitochondrial and chloroplast genomes and the nonrecombining portion of the Y chromosome have provided, to date, most of the known data regarding compound haplotypic variation within and among populations. These comparatively small segments include numerous polymorphic sites and undergo little or no recombination. Recombining autosomes, however, comprise the major repository of genetic variation. Technical challenges and recombination have limited large-scale application of autosomal haplotypes. We have overcome this barrier through development of a general approach to the assessment of short autosomal DNA segments. Each such segment includes one or more single nucleotide polymorphisms (SNPs) and exactly one short tandem repeat (STR) locus. With dramatically different mutation rates, these two types of genetic markers provide complementary evolutionary information. We call the combination of a SNP and a STR polymorphism a SNPSTR, and have developed a simple, rapid method for empirically determining gametic phase for double and triple heterozygotes. Here, we illustrate the approach with two SNPSTR systems. Although even one system provides insight into population history, the power of the approach lies in combining results from multiple SNPSTR systems.

The present study is a contribution to the definition of the linkage disequilibrium relationship of MICA and MICB with adjacent loci and to the characterization of extended HLA haplotypes. These issues are of importance for the identification of disease associations and for a better definition of donor-recipient compatibility in bone-marrow grafts through the typing of haplospecific markers. The distribution of the five alleles of MICA and the 13 alleles of MICB microsatellites, located, respectively, in MICA transmembrane exon 5 and in MICB intron 1, was examined in 133 healthy Italian individuals previously typed for HLA class I, class II and complement loci and for the TNFa microsatellite. The MICB microsatellite was also analysed in 49 HTCLs for which MICA typing was already available. Very strong linkage disequilibria with HLA-B and TNFa were detected in the Italian population for both MICA and MICB microsatellite alleles, in spite of the high mutability rate of the larger MICB alleles. Some strong associations were also detected between MICB and DRB1. The strongest associations (P < 0.001, D' > 0.7) were those of MICA-A4 with HLA-B18, B27 and TNFa1, MICA-A5 with HLA-B35, B61 and B62, MICA-A5.1 with HLA-B7, B8, B13, B63 and MICB-CA24, MICA-A6 with HLA-B51, MICA-A9 with HLA-B39, B57 and TNFa2, MICB-CA14 with HLA-B14, B27 and TNFa1, MICB-CA15 with HLA-B52, TNFa4 and TNFa13, MICB-CA17 with HLA-B7 and TNFa11, MICB-CA18 with HLA-B13 and TNFa7, MICB-CA22 with HLA-B57, and MICB-CA24 with HLA-B8 and TNFa2. From pairwise associations in the random panel and results for the homozygous cell lines it was possible to deduce the MICA and MICB microsatellite alleles present in many of the well-known Caucasoid extended haplotypes.

Advances in sequencing and genotyping technologies over the last decade have enabled geneticists to easily characterize genetic variation at the nucleotide level. Hundreds of genes harboring mutations associated with genetic disease have now been identified by positional cloning. Using variation at closely linked genetic markers, it is possible to predict the times in the past at which particular mutations arose. Such studies suggest that many of the rare mutations underlying human genetic disorders are relatively young. Studies of variation at genetic markers linked to particular mutations can provide insights into human geographic history, and historical patterns of natural selection and disease, that are not available from other sources. We review two approaches for estimating allele age using variation at linked genetic markers. A phylogenetic approach aims to reconstruct the gene tree underlying a sample of chromosomes carrying a particular mutation, obtaining a "direct" estimate of allele age from the age of the root of this tree. A population genetic approach relies on models of demography, mutation, and/or recombination to estimate allele age without explicitly reconstructing the gene tree. Phylogenetic methods are best suited for studies of ancient mutations, while population genetic methods are better suited for studies of recent mutations. Methods that rely on recombination to infer the ages of alleles can be fine-tuned by choosing linked markers at optimal map distances to maximize the information available about allele age. A limitation of methods that rely on recombination is the frequent lack of a fine-scale linkage map. Maximum likelihood and Bayesian methods for estimating allele age that rely on intensive numerical computation are described, as well as "composite" likelihood and moment-based methods that lead to simple estimators. The former provide more accurate estimates (particularly for large samples of chromosomes) and should be employed if computationally practical.

Most evidence about nonrandom association of alleles at different loci, or gametic disequilibrium, across extensive anonymous regions of the human genome is based on the analysis of overall disequilibrium between pairs of microsatellites. However, analysis of interallelic associations is also necessary for a more complete description of disequilibrium. Here, we report a study characterizing the frequency and strength of both overall and interallelic disequilibrium between pairs of 12 microsatellite loci (CA repeats) spanning 19 cM (14 Mb) on human chromosome 11p15, in a large sample (810 haplotypes deduced from 405 individuals) drawn from a single population. Characterization of disequilibrium was carried out, taking into account the sign of the observed disequilibria. This strategy facilitates detection of associations and gives more accurate estimates of their intensities. Our results demonstrate that the incidence of disequilibrium over an extensive human chromosomal region is much greater than is commonly considered for populations that have expanded in size. In total, 44% of the pairs of microsatellite loci and 18% of the pairs of alleles showed significant nonrandom association. All the loci were involved in disequilibrium, although both the frequency and strength of interallelic disequilibrium were distributed nonuniformly along 11p15. These findings are especially relevant since significant associations were detected between loci separated by as much as 17-19 cM (7 cM on average). It was also found that the overall disequilibrium masks complicated patterns of association between pairs of alleles, dependent on their frequency and size. We suggest that the complex mutational dynamics at microsatellite loci could explain the allele-dependent disequilibrium patterns. These observations are also relevant to evaluation of the usefulness of microsatellite markers for fine-scale localization of disease genes.

The extent of microsatellite size homoplasy, as well as its effect on several population genetics statistics, was investigated in natural populations using the single-strand conformation polymorphism (SSCP) method. The analysis was conducted using 240 individuals from 13 populations of the freshwater snail Bulinus truncatus at a GT(n)CT(m) compound microsatellite locus. We showed that SSCP can be used to uncover, at least partly, size homoplasy in the core sequence of this category of loci. Eight conformers (SSCP variants) were detected among the three size variants (electromorphs). Sequencing revealed that each conformer corresponded to a different combination of repeats in the GT(n) and CT(m) arrays. Part of this additional variability was detected within populations, resulting in a substantial increase in gene diversity in four populations. Additional variability also changed the values of parameters used to analyze population differentiation among populations: pairwise tests of differentiation were significant much more often with conformers than with electromorphs. On the other hand, pairwise estimates of F(st) were either smaller or larger with conformers than with electromorphs, depending on whether or not electromorphs were shared among populations. However, estimates of F(st) (or analogs) over all populations were very similar, ranging between 0.66 and 0.75. Our results were consistent with the theoretical prediction that homoplasy should not always lead to stronger population structure. Finally, conformer sequences and electromorph size distribution suggested that single-point and/or stepwise mutations occurring simultaneously in the different repeated arrays of compound microsatellites produce sequence variation without size variation and hence generate more size homoplasy than expected under a simple stepwise mutation model.

A novel microsatellite homologous to DYS391, a (GATA)(n) short tandem repeat on the human Y chromosome, was identified and characterized in the present work. Employing somatic cell hybrid and deletion panels in a PCR-based approach, we found out that the new microsatellite is located in Xp21.2-22.3, while its Y counterpart mapped to Yq11.21. This X-linked locus (provisionally called DXYS391) and its Y homolog constitute one more example of similarity outside the pseudoautosomal regions between the two human sex chromosomes. Sequencing data showed high levels of homology in the flanking regions of DXYS391 and DYS391 that differ primarily by the presence of a (GACA)(3) motif in the Y locus. Both loci were detected in chimpanzee DNA, suggesting that a putative transposition from the X to the Y occurred before the human/chimpanzee split. The allele frequencies of DYS391 and DXYS391 were investigated, respectively, in 271 Y and 337 X chromosomes from distinct human populations worldwide. DYS391 consistently displayed greater among-population component of the variance of the allele frequencies than DXYS391, as expected due to the three-times lower effective population size of Y chromosomes relative to the X. The intra-population diversity of DYS391, measured by Nei's locus diversity as well as by allele size variance, was lowest in Amerindians, while very low diversity of DXYS391 was seen in Africans. Since our African data are based on a small sample, further studies will be necessary to evaluate better this observation.

DNA sequence polymorphism and divergence was examined in the vicinity of the human beta-globin gene cluster origin of replication initiation region (IR), a 1.3-kb genomic region located immediately 5' of the adult-expressed beta-globin gene. DNA sequence variation in the replication origin IR and 5 kb of flanking DNA was surveyed in samples drawn from two populations, one African (from the Gambia, West Africa) and the other European (from Oxford, England). In these samples, levels of nucleotide and length polymorphism in the IR were found to be more than two times as high as adjacent non-IR-associated regions (estimates of per-nucleotide heterozygosity were 0.30% and 0.12%, respectively). Most polymorphic positions identified in the origin IR fall within or just adjacent to a 52-bp alternating purine-pyrimidine ((RY)n) sequence repeat. Within- and between-populations divergence is highest in this portion of the IR, and interspecific divergence in the same region, determined by comparison with an orthologous sequence from the chimpanzee, is also pronounced. Higher levels of diversity in this subregion are not, however, primarily attributable to slippage-mediated repeat unit changes, as nucleotide substitution contributes disproportionately to allelic heterogeneity. An estimate of helical stability in the sequenced region suggests that the hypervariable (RY)n constitutes the major DNA unwinding element (DUE) of the replication origin IR, the location at which the DNA duplex first unwinds and new strand synthesis begins. These findings suggest that the beta-globin IR experiences a higher underlying rate of neutral mutation than do adjacent genomic regions and that enzyme fidelity associated with the initiation of DNA replication at this origin may be compromised. The significance of these findings for our understanding of eukaryotic replication origin biology is discussed.

We analyzed the origin of allelic diversity at the class II HLA-DRB1 locus, using a complex microsatellite located in intron 2, close to the polymorphic second exon. A phylogenetic analysis of human, gorilla, and chimpanzee DRB1 sequences indicated that the structure of the microsatellite has evolved, primarily by point mutations, from a putative ancestral (GT)x(GA)y-complex-dinucleotide repeat. In all contemporary DRB1 allelic lineages, with the exception of the human *04 and the gorilla *08 lineages, the (GA)y repeat is interrupted, often by a G-->C substitution. In general, the length of the 3' (GA)y repeat correlates with the allelic lineage and thus evolves more slowly than a middle (GA)z repeat, whose length correlates with specific alleles within the lineage. Comparison of the microsatellite sequence from 30 human DRB1 alleles showed the longer 5' (GT)x to be more variable than the shorter middle (GA)z and 3' (GA)y repeats. Analysis of multiple samples with the same exon sequence, derived from different continents, showed that the 5' (GT)x repeat evolves more rapidly than the middle (GA)z and the 3' (GA)y repeats, which is consistent with findings of a higher mutation rate for longer tracts. The microsatellite-repeat-length variation was used to trace the origin of new DRB1 alleles, such as the new *08 alleles found in the Cayapa people of Ecuador and the Ticuna people of Brazil.

The genomic region encompassing the Major Histocompatibility Complex (MHC) contains polymorphic frozen blocks which have developed by local imperfect sequential duplication associated with insertion and deletion (indels). In the alpha block surrounding HLA-A, there are ten duplication units or beads on the 62.1 ancestral haplotype. Each bead contains or contained sequences representing Class I, PERB11 (MHC Class I chain related (MIC) and human endogenous retrovirus (HERV) 16. Here we consider explanations for co-occurrence of genomic polymorphism, duplication and HERVs and we ask how these features encode susceptibility to numerous and very diverse diseases. Ancestral haplotypes differ in their copy number and indels in addition to their coding regions. Disease susceptibility could be a function of all of these differences. We propose a model of the evolution of the human MHC. Population-specific integration of retroviral sequences could explain rapid diversification through duplication and differential disease susceptibility. If HERV sequences can be protective, there are exciting prospects for manipulation. In the meanwhile, it will be necessary to understand the function of MHC genes such as PERB11 (MIC) and many others discovered by genomic sequencing.

Sequencing studies were performed in three previously described microsatellite and minisatellite markers located within the HLA-DQ region, DQCAR, DQCARII and G51152. Multiple nucleotide substitutions that did not change size polymorphisms were observed in all three markers. In all loci, the number of core repeats did not correlate with neighboring DQ allele sequence motifs while single nucleotide changes within or flanking the microsatellite sequence did. This result indicates higher mutation rates for microsatellite expansions/contractions than for nucleotide substitutions in these loci. Further analysis indicated an almost complete phylogenetic correspondence between DQCAR single nucleotide polymorphisms (SNPs) and DQB1 sequences on one side (1.0-1.5 kb apart) and a complete relationship between DQCARII and DQA1 sequences on the other (4.5 kb apart). In contrast, G51152 sequences did not correspond perfectly with DQB1 allelic sequences, thus suggesting the existence of several ancestral crossovers between this marker and DQB1 (20-25 kb). Sequencing microsatellites might be useful in disease mapping studies by increasing marker informativeness and by helping in the interpretation of association study results. It is also proposed that SNPs within the flanking region of CA repeats could be used to develop biallelic markers from already available mapped microsatellite markers.