The association between HLA-A alleles and an Alu dimorphism near HLA-G

The AluYb8 sequences are a subfamily of short interspersed Alu retroelements that have been amplified within the human genome during recent evolutionary time and are useful polymorphic markers for studies on the origin of human populations. We have identified a new member of the Yb8 subfamily, AluyHG, located between the HLA-H and -G genes and 88-kb telomeric of the highly polymorphic HLA-A gene within the alpha block of the major histocompatibility complex (MHC). The AluyHG element was characterised with a view to examining the association between AluyHG and HLA-A polymorphism and reconstructing the history of the MHC alpha block. A specific primer pair was designed for a simple PCR assay to detect the absence or presence (dimorphism) of the AluyHG element within the DNA samples prepared from a panel of 46 homozygous cell-lines containing complete or recombinant ancestral haplotypes (AH) of diverse ethnic origin and 92 Caucasoid and Asian subjects on which HLA-A typing was available. The AluyHG insertion was most strongly associated with HLA-A2 and, to a lesser degree with HLA-A1, -A3, -A11, and A-19. The gene frequency of the AluyHG insertion for 146 Caucasians and 94 Chinese-Han was 0.30 and 0.32 and there was no significant difference between the observed and expected frequencies. The results of the association studies and the phylogenetic analysis of HLA-A alleles suggest that the AluyHG sequence was integrated within the progenitor of HLA-A2, but has been transferred by recombination to other human ancestral populations. In this regard, the dimorphic AluyHG element is an important diagnostic marker for HLA association studies and could help in elucidating the evolution and functions of the MHC alpha block and polymorphism within and between ancestral haplotypes. Received: 7 December 2000 / Accepted: 28 February 2001     

Lineage Specificity of the Evolutionary Dynamics of the mtDNA D-Loop Region in Rodents

This paper reports an intraorder study on the D-loop-containing region of the mitochondrial DNA in rodents. A complete multialignment of this region is not feasible with the exception of some conserved regions. The comparative analysis of 25 complete rodent sequences from 23 species plus one lagomorph has revealed that only the central domain (CD), a conserved region of about 80 bp in the extended termination-associated sequences (ETAS) domain, adjacent to the CD, the ETAS1, and conserved sequence block (CSB) 1 blocks are present in all rodent species, whereas the presence of CSB2 and CSB3 is erratic within the order. We have also found a conserved region of 90 bp located between tRNAPro and ETAS1 present in fat dormouse, squirrel, guinea pig, and rabbit. Repeated sequences are present in both the ETAS and the CSB domain, but the repeats differ in length, copy number, and base composition in different species. The potential use of the D-loop for evolutionary studies has been investigated; the presence/absence of conserved blocks and/or repeated sequences cannot be used as a reliable phylogenetic marker, since in some cases they may be shared by distantly related organisms but not by close ones, while in other ones a relationship between tree topology and presence/absence of such motifs is observed. Better results can be obtained by the use of the CD, which, however, due to its reduced size, when used for tracing a phylogenetic tree, shows some nodes with low statistical support. Received: 26 February 2001 / Accepted: 6 June 2001     

Coevolution of HLA-B and PERB11.1 (MICA): Significance of Independent Triplet Expansion Within the Transmembrane Region of PERB11.1 (MICA)

Several highly polymorphic sequences are present in the beta block of the MHC, especially HLA-B, HLA-C, PERB11.1 (MICA), and PERB11.2 (MICB). It is now apparent that the polymorphism of PERB11.1 is of the same order as that of HLA-A, -B, and -C and it has been suggested that PERB11 could explain some of the disease associations previously attributed to HLA-B. Phylogenetic analysis of PERB11 α-domain sequences demonstrates relationships with HLA-B cross-reactive serogroups. In contrast, the transmembrane polymorphisms do not appear to be associated with either PERB11 or HLA-B. These data indicate that PERB11 and HLA-B have evolved in concert from their common ancestors and that the transmembrane polymorphisms have arisen independently and more recently. MHC disease associations will need to be reviewed in the light of mechanisms such as receptor binding and signaling. Received: 22 September 1999 / Accepted: 11 January 2000 (Revised)     

Genetic analysis of completely sequenced disease-associated MHC haplotypes identifies shuffling of segments in recent human history

The major histocompatibility complex (MHC) is recognised as one of the most important genetic regions in relation to common human disease. Advancement in identification of MHC genes that confer susceptibility to disease requires greater knowledge of sequence variation across the complex. Highly duplicated and polymorphic regions of the human genome such as the MHC are, however, somewhat refractory to some whole-genome analysis methods. To address this issue, we are employing a bacterial artificial chromosome (BAC) cloning strategy to sequence entire MHC haplotypes from consanguineous cell lines as part of the MHC Haplotype Project. Here we present 4.25 Mb of the human haplotype QBL (HLA-A26-B18-Cw5-DR3-DQ2) and compare it with the MHC reference haplotype and with a second haplotype, COX (HLA-A1-B8-Cw7-DR3-DQ2), that shares the same HLA-DRB1, -DQA1, and -DQB1 alleles. We have defined the complete gene, splice variant, and sequence variation contents of all three haplotypes, comprising over 259 annotated loci and over 20,000 single nucleotide polymorphisms (SNPs). Certain coding sequences vary significantly between different haplotypes, making them candidates for functional and disease-association studies. Analysis of the two DR3 haplotypes allowed delineation of the shared sequence between two HLA class II-related haplotypes differing in disease associations and the identification of at least one of the sites that mediated the original recombination event. The levels of variation across the MHC were similar to those seen for other HLA-disparate haplotypes, except for a 158-kb segment that contained the HLA-DRB1, -DQA1, and -DQB1 genes and showed very limited polymorphism compatible with identity-by-descent and relatively recent common ancestry (<3,400 generations). These results indicate that the differential disease associations of these two DR3 haplotypes are due to sequence variation outside this central 158-kb segment, and that shuffling of ancestral blocks via recombination is a potential mechanism whereby certain DR-DQ allelic combinations, which presumably have favoured immunological functions, can spread across haplotypes and populations.     

Locus-specific conservation of the HLA class I introns by intra-locus homogenization

 The protein-coding sequences of the major histocompatibility complex (MHC) genes are characterized by extraordinarily high polymorphism, apparently maintained by balancing selection, which favors diversity in the peptide-binding domains of the MHC glycoproteins. Here we report that the introns flanking the polymorphic exons of the human MHC class I loci HLA-A, -B, and -C genes have been relatively conserved and have become locus-specific apparently as a result of recombination and subsequent genetic drift, leading to homogenization within loci over evolutionary time. Thus, HLA class I genes have been shaped by contrasting evolutionary forces maintaining polymorphism in the exons and leading to conservation in the introns. This study provides the first extensive analysis of the introns of a highly polymorphic gene family. Received: 10 April 1997 / Revised: 15 July 1997     

Duplication and Diversification of the Apolipoprotein CI (APOCI) Genomic Segment in Association with Retroelements

We have previously shown that several multicopy gene families within the major histocompatibility complex (MHC) arose from a process of segmental duplication. It has also been observed that retroelements play a role in generating diversity within these duplicated segments. The objective of this study was to compare the genomic organization of a gene duplication within another multicopy gene family outside the MHC. Using new continuous genomic sequence encompassing the APOE-CII gene cluster, we show that APOCI and its pseudogene, APOCI′, are contained within large duplicated segments which include sequences from the hepatic control region (HCR). Flanking Alu sequences are observed at both ends of the duplicated unit, suggesting a possible role in the integration of these segments. As observed previously within the MHC, the major differences between the segments are the insertion of sequences (approximately 200–1000 bp in length), consisting predominantly of Alu sequences. Ancestral retroelements also contribute to the generation of sequence diversity between the segments, especially within the 3′ poly(A) tract of Alu sequences. The exonic and regulatory sequences of the APOCI and HCR loci show limited sequence diversity, with exon 3 being an exception. Finally, the typing of pre- and postduplication Alus from both segments indicates an estimated time of duplication of approximately 37 million years ago (mya), some time prior to the separation of Old and New World monkeys. Received: 17 July 1999 / Accepted: 6 November 1999     

Using Alu J Elements as Molecular Clocks to Trace the Evolutionary Relationships Between Duplicated HLA Class I Genomic Segments

The class I region of the major histocompatibility complex contains two subgenomic blocks (250–350 kb each), known as the alpha and beta blocks. These blocks contain members of multicopy gene families including HLA class I, HERV-16 (previously called P5 sequences), and PERB11 (MIC). We have previously shown that each block consists of imperfect duplicated segments (duplicons) containing linked members of different gene families, retroelements and transposons that have coevolved as part of two separate evolutionary events. Another region provisionally designated here as the kappa block is located between the alpha and the beta blocks and contains HLA-E, -30, and -92, HERV-16 (P5.3), and PERB11.3 (MICC) within about 250 kb of sequence. Using Alu elements to trace the evolutionary relationships between different class I duplicons, we have found that (a) the kappa block contains paralogous (duplicated) Alu J sequences and other retroelement patterns more in common with the beta than the alpha block; (b) the retroelement pattern associated with the HLA-E duplicon is different from all other HLA class I duplicons, indicating a more complex evolution; (c) the HLA-92 duplicon, although substantially shorter, is closely related in sequence to the HLA-B and -C duplicons; (d) two of the six paralogous Alu J elements within the HLA-B and -C duplicons are associated with the HLA-X duplicon, confirming their evolutionary relationships within the beta block; and (e) the paralogous Alu J elements within the alpha block are distinctly different from those identified within the beta and kappa blocks. The sequence conservation and location of duplicated (paralogous) Alu J elements in the MHC class I region show that the beta and kappa blocks have evolved separately from the alpha block beginning at a time before or during the evolution of Alu J elements in primates. Received: 22 September 1999 / Accepted: 24 January 2000     

mtDNA Sequence Diversity of Orangutans from the Islands of Borneo and Sumatra

A comparison of mitochondrial DNA sequences was undertaken for two genes among orangutans from Borneo and Sumatra. The distribution of haplotypes among 42 individuals for NADH dehydrogenease subunit 3 and 39 individuals for cytochrome B was used to infer population structure. The haplotypes among all Bornean orangutans form a cluster of closely related individuals for both genes, with two distinct haplotypes occupying different regions of the island. Sumatran haplotypes fall into three distinct, and highly diverged, groups. Strikingly, one of the Sumatran haplotypes shares sequence identity with the most widespread Bornean haplotype. This haplotype distribution is considered in the context of the highly complex geological history for the area around the Malay Archipelago. Alternating periods of geographic isolation and reunion, resulting from glacially induced land bridge formation, presented substantial opportunity for population dispersal between periodically isolated demes. We present a paleodispersal model that is consistent with genetic, geological, paleoecological, and fossil data. The disparity of sequences present in orangutan populations argues against a simple Sumatra–Borneo dichotomy. Our evidence, and that of others, suggests that the Sumatran population alone contains the remnants of at least three separate lineages. Received: 3 February 1999 / Accepted: 2 August 2000     

Diversity and locus specificity of chicken MHC B class I sequences

The major histocompatibility complex B (MHC B) region in a standard haplotype of Leghorn chickens contains two closely linked class I loci, B-FI and B-FIV. Few sequences of B-FI alleles are available, and therefore alleles of the two loci have not been compared with regard to sequence diversity or locus specificity. Here, we report eight new B-F alpha 1/alpha 2-coding sequences from broiler chicken MHC B haplotypes, and a unique recombinant between the two B-F loci. The new sequences were combined with existing B-F sequences from Leghorn and broiler haplotypes for analysis. On the basis of phylogenetic analysis and conserved sequence motifs, B-F sequences separated into two groups (Groups A and B), corresponding to B-FIV and B-FI locus, respectively. Every broiler haplotype had one B-F sequence in Group A and the second B-F sequence, if it existed, clustered in Group B. Group B (presumptive B-FI locus) sequences identified in broiler haplotypes resembled the human MHC class I HLA-C locus in their distinctive pattern of allelic polymorphism. Compared with B-FIV, B-FI alleles were less polymorphic and possessed a conserved locus-specific motif in the alpha1 helix, but nevertheless demonstrated evidence of diversifying selection. One B-FI alpha 1/alpha 2-coding nucleotide sequence was completely conserved in four different broiler haplotypes, but each allele differed in the exon encoding the alpha 3 domain.     

Ancestral and recombinant 16-locus HLA haplotypes in the Hutterites

 Prior studies in the Schmiedeleut Hutterites of South Dakota have demonstrated associations between human leukocyte antigen (HLA) haplotype matching and fetal loss (Ober et al. 1992) and mate preferences (Ober et al. 1997), as well as deficiencies of homozygotes for HLA haplotypes (Kostyu et al. 1993). These studies were based on the serologically-defined five-locus HLA-A, -C, -B, -DR, -DQ haplotype. To further elucidate the effects of specific major histocompatibility (MHC) loci or regions on fetal loss and mate choice, we genotyped a sample of Hutterites for 14 MHC loci by DNA or biochemical methods. Typing for additional loci in the HLA-A to HLA-DPB1 region increased the number of recognized Hutterite MHC haplotypes to 67, and further localized the site of cross-over in 9 of 15 recombinant haplotypes. Hutterite MHC haplotype sequences are similar to those observed in outbred Caucasians, suggesting that the influence of HLA haplotypes on fetal loss and mating structure may be general. Received: 1 May 1998 / Revised: 2 December 1998     

Dominant Sequences of Human Major Histocompatibility Complex Conserved Extended Haplotypes from HLA-DQA2 to DAXX

We resequenced and phased 27 kb of DNA within 580 kb of the MHC class II region in 158 population chromosomes, most of which were conserved extended haplotypes (CEHs) of European descent or contained their centromeric fragments. We determined the single nucleotide polymorphism and deletion-insertion polymorphism alleles of the dominant sequences from HLA-DQA2 to DAXX for these CEHs. Nine of 13 CEHs remained sufficiently intact to possess a dominant sequence extending at least to DAXX, 230 kb centromeric to HLA-DPB1. We identified the regions centromeric to HLA-DQB1 within which single instances of eight “common” European MHC haplotypes previously sequenced by the MHC Haplotype Project (MHP) were representative of those dominant CEH sequences. Only two MHP haplotypes had a dominant CEH sequence throughout the centromeric and extended class II region and one MHP haplotype did not represent a known European CEH anywhere in the region. We identified the centromeric recombination transition points of other MHP sequences from CEH representation to non-representation. Several CEH pairs or groups shared sequence identity in small blocks but had significantly different (although still conserved for each separate CEH) sequences in surrounding regions. These patterns partly explain strong calculated linkage disequilibrium over only short (tens to hundreds of kilobases) distances in the context of a finite number of observed megabase-length CEHs comprising half a population''s haplotypes. Our results provide a clearer picture of European CEH class II allelic structure and population haplotype architecture, improved regional CEH markers, and raise questions concerning regional recombination hotspots.     

Population Genetics of Y-Chromosome Short Tandem Repeats in Humans

Eight human short tandem repeat polymorphisms (STRs) also known as microsatellites—DYS19, DYS388, DYS390, DYS391, DYS392, DYS393, DYS389I, and DYS389II, mapping in the Y chromosome—were analyzed in two Iberian samples (Basques and Catalans). Allele frequency distributions showed significant differences only for DYS392. Fst and gene diversity index (D) were estimated for the Y STRs. The values obtained are comparable to those of autosomal STR if corrections for the smaller effective population size on the Y chromosome are taken into account. This suggests that Y-chromosome microsatellites might be as useful as their autosomal counterparts to both human population genetics and forensics. Our results also reinforce the hypothesis that selective sweeps in the Y chromosome in recent times are unlikely. Haplotypes combining five of the loci were constructed for 71 individuals, showing 29 different haplotypes. A haplotype tree was constructed, from which an estimate of 7,000 to 60,000 years for the age of the Y-chromosome variation in Iberia was derived, in accordance with previous estimates obtained with mtDNA sequences and nuclear markers. Received: 3 January 1997 / Accepted: 25 April 1997     

Distribution of recombination crossovers and the origin of haplotype blocks: the interplay of population history,recombination, and mutation

Recent studies suggest that haplotypes are arranged into discrete blocklike structures throughout the human genome. Here, we present an alternative haplotype block definition that assumes no recombination within each block but allows for recombination between blocks, and we use it to study the combined effects of demographic history and various population genetic parameters on haplotype block characteristics. Through extensive coalescent simulations and analysis of published haplotype data on chromosome 21, we find that (1) the combined effects of population demographic history, recombination, and mutation dictate haplotype block characteristics and (2) haplotype blocks can arise in the absence of recombination hot spots. Finally, we provide practical guidelines for designing and interpreting studies investigating haplotype block structure.     

Conserved Gene Clusters in Bacterial Genomes Provide Further Support for the Primacy of RNA

Five complete bacterial genome sequences have been released to the scientific community. These include four (eu)Bacteria, Haemophilus influenzae, Mycoplasma genitalium, M. pneumoniae, and Synechocystis PCC 6803, as well as one Archaeon, Methanococcus jannaschii. Features of organization shared by these genomes are likely to have arisen very early in the history of the bacteria and thus can be expected to provide further insight into the nature of early ancestors. Results of a genome comparison of these five organisms confirm earlier observations that gene order is remarkably unpreserved. There are, nevertheless, at least 16 clusters of two or more genes whose order remains the same among the four (eu)Bacteria and these are presumed to reflect conserved elements of coordinated gene expression that require gene proximity. Eight of these gene orders are essentially conserved in the Archaea as well. Many of these clusters are known to be regulated by RNA-level mechanisms in Escherichia coli, which supports the earlier suggestion that this type of regulation of gene expression may have arisen very early. We conclude that although the last common ancestor may have had a DNA genome, it likely was preceded by progenotes with an RNA genome. Received: 10 March 1996 / Accepted: 20 May 1997     

Unusual and Strongly Structured Sequence Variation in a Complex Satellite DNA Family from the Nematode Meloidogyne chitwoodi

An AluI satellite DNA family has been isolated in the genome of the root-knot nematode Meloidogyne chitwoodi. This repeated sequence was shown to be present at approximately 11,400 copies per haploid genome, and represents about 3.5% of the total genomic DNA. Nineteen monomers were cloned and sequenced. Their length ranged from 142 to 180 bp, and their A + T content was high (from 65.7 to 79.1%), with frequent runs of As and Ts. An unexpected heterogeneity in primary structure was observed between monomers, and multiple alignment analysis showed that the 19 repeats could be unambiguously clustered in six subfamilies. A consensus sequence has been deduced for each subfamily, within which the number of positions conserved is very high, ranging from 86.7% to 98.6%. Even though blocks of conserved regions could be observed, multiple alignment of the six consensus sequences did not enable the establishment of a general unambiguous consensus sequence. Screening of the six consensus sequences for evidence of internal repeated subunits revealed a 6-bp motif (AAATTT), present in both direct and inverted orientation. This motif was found up to nine times in the consensus sequences, also with the occurrence of degenerated subrepeats. Along with the meiotic parthenogenetic mode of reproduction of this nematode, such structural features may argue for the evolution of this satellite DNA family either (1) from a common ancestral sequence by amplification followed by mechanisms of sequence divergence, or (2) through independent mutations of the ancestral sequence in isolated amphimictic nematode populations and subsequent hybridization events. Overall, our results suggest the ancient origin of this satellite DNA family, and may reflect for M. chitwoodi a phylogenetic position close to the ancestral amphimictic forms of root-knot nematodes. Received: 23 April 1997 / Accepted: 9 July 1997