Genomic and Phylogenetic Analysis of the Human CD1 and HLA Class I Multicopy Genes

The human CD1 proteins belong to a lipid-glycolipid antigen-presenting gene family and are related in structure and function to the MHC class I molecules. Previous mapping and DNA hybridization studies have shown that five linked genes located within a cluster on human chromosome 1q22-23 encode the CD1 protein family. We have analyzed the complete genomic sequence of the human CD1 gene cluster and found that the five active genes are distributed over 175,600 nucleotides and separated by four expanded intervening genomic regions (IGRs) ranging in length between 20 and 68 kb. The IGRs are composed mostly of retroelements including five full-length L1 PA sequences and various pseudogenes. Some L1 sequences have acted as receptors for other subtypes or families of retroelements. Alu molecular clocks that have evolved during primate history are found distributed within the HLA class I duplicated segments (duplicons) but not within the duplicons of CD1. Phylogeny of the alpha3 domain of the class I-like superfamily of proteins shows that the CD1 cluster is well separated from HLA class I by a number of superfamily members including MIC (PERB11), HFE, Zn-alpha2-GP, FcRn, and MR1. Phylogenetically, the human CD1 sequences are interspersed by CD1 sequences from other mammalian species, whereas the human HLA class I sequences cluster together and are separated from the other mammalian sequences. Genomic and phylogenetic analyses support the view that the human CD1 gene copies were duplicated prior to the evolution of primates and the bulk of the HLA class I genes found in humans. In contrast to the HLA class I genomic structure, the human CD1 duplicons are smaller in size, they lack Alu clocks, and they are interrupted by IGRs at least 4 to 14 times longer than the CD1 genes themselves. The IGRs seem to have been created as "buffer zones" to protect the CD1 genes from disruption by transposable elements.     

Interchromosomal duplication of major histocompatibility complex class I regions in rainbow trout (Oncorhynchus mykiss), a species with a presumably recent tetraploid ancestry

Salmonid fishes are among the few animal taxa with a probable recent tetraploid ancestor. The present study is the first to compare large (>100 kb) duplicated genomic sequence fragments in such species. Two contiguous stretches with major histocompatibility complex (MHC) class I genes were detected in a rainbow trout BAC library, mapped and sequenced. The MHC class I duplicated regions, mapped by fluorescence in situ hybridization (FISH), were shown to be located on different metaphase chromosomes, Chr 14 and 18. Gene organization in both duplications is similar to that in other fishes, in that the class I loci are tightly linked with the PSMB8, PSMB9, PSMB10 and ABCB3 genes. Whereas one region, Onmy-IA, has a classical MHC class I locus (UBA), Onmy-IB encodes only non-classical class Ib proteins. The nucleotide diversity between the Onmy-IA and Onmy-IB noncoding regions is about 14%. This suggests that the MHC class I duplication event has occurred about 60 mya close to the time of an hypothesized ancestral tetraploid event. The present article is the first convincing report on the co-existence of two closely related MHC class I core regions on two different chromosomes. The interchromosomal duplication and the homology levels are supportive of the tetraploid model.Nucleotide sequence data reported are available in the DDBJ/EMBL/GenBank database under the accession numbers AB162342, AB162343 and from AY525774 to AY525776.     

Induction of genetic duplications and frameshift mutations in Salmonella typhimurium by acridines and acridine mustards: dependence on covalent binding of the mutagen to DNA

Summary The aroC321 allele permits positive selection for the detection of a large genetic duplication that arises in the Salmonella typhimurium chromosome by homologous recombination. Strains that contain both aroC321 and the hisC3076 allele were constructed so that the induction of genetic duplications and frameshift mutations in a run of GC base pairs could be studied simultaneously by selecting for tryptophan and histidine prototrophy, respectively. Using these strains, we examined the ability of 9-aminoacridine, quinacrine, four acridine mustards (ICR-170, ICR-191, ICR-372, and quinacrine mustard) and the nitroacridine Entozon to induce genetic duplications and frameshift mutations. Although all these compounds induce reversion of hisC3076, only the four mustards and Entozon are effective as inducers of genetic duplications under identical treatment conditions. The induction of genetic duplications by acridine mustards, like the toxic and mutagenic effects of these compounds, is enhanced by a deficiency for excision repair caused by a deletion through the uvrB gene. The ineffectiveness of 9-aminoacridine and quinacrine in the test for genetic duplications indicates that simple intercalation is sufficient for the mutagenic effect measured with the hisC3076 allele but that the induction of duplications by the acridine mustards and Entozon requires covalent binding of the chemical to DNA.     

Multiplex real-time PCR for detection of deletions and duplications in dystrophin gene

Genetic testing of Duchenne and Becker muscular dystrophies (DMD/BMD) is a difficult task due to the occurrence of deletions or duplications within dystrophin (DMD) gene that requires dose sensitive tests. We developed three multiplex quantitative real-time PCR assays for dystrophin exon 5, 45, and 51 within two major hotspots of deletion/duplication. Each exon was co-amplified with a reference X-linked gene and the copy number of the target fragment was calculated by comparative threshold cycle method (delta deltaC(t)). We compared the performance of this method with previously described end-point PCR fluorescent analysis (EPFA) by studying 24 subjects carrying DMD deletions or duplications. We showed that Q-PCR is an accurate and sensitive technique for the identification of deletions and duplications in DMD/BMD. Q-PCR is a valuable tool for independent confirmation of EPFA screening, particularly when deletions/duplications of single exons occur or for rapid identification of known mutations in at risk carriers.     

Comparison Between Two Human Endogenous Retrovirus (HERV)-Rich Regions Within the Major Histocompatibility Complex

Sixteen human endogenous retrovirus (HERV) sequences were detected within 656 kb of genomic sequence obtained from the alpha- and beta-block of the class I region of the major histocompatibility complex (MHC). The HERVs were identified and characterized as family members of HERV-16 (11 copies), HERV-L (1 copy), HERV-I (2 copies), HERV-K91 (1 copy), and HARLEQUIN (1 copy) by sequence comparison using CENSOR or Repeat Masker, BLAST searches, and dot plots. The 11 copies of HERV-16 arose as products of duplication of genomic segments containing HLA class I (HLAcI) and PERB11 (MIC) genes inter alia, whereas the other five HERVs arose after duplication probably as a consequence of single insertion events or translocations. HERV-L and HERV-I are located between the duplicated genes PERB11.2 (MICB) and PERB11.1 (MICA), and HLA-B and HLA-C, respectively, whereas HERV-K91 and HARLEQUIN are located telomeric of HLA-C. A highly fragmented copy of HERV-I was also found telomeric of PERB11.4. Structural analysis of open reading frames (ORFs) revealed the absence of intact coding sequence within the putative gag, pol, and env gene regions of all the HERVs with the exception of HERV-K91, which had two large ORFs within the region of the putative protease and pol genes. In addition, the 5′-LTR of HERV-L contained a 2.5-kb element that was AT-rich and large ORFs with putative amino acid sequences rich in tyrosines and isoleucines. HERV-I, HARLEQUIN, and at least four copies of HERV-16 appear to have been receptors for the insertion of other retrotransposons including Alu elements and fragments of L1 and THE1. Examination of flanking sequences suggests that HERV-I and HERV-L had occurred by insertion into ancient L1 fragments. This study has revealed that the alpha- and beta-block region within the MHC is rich in HERV sequences occurring at a much higher ratio (10 to 1) than normally observed in the human genome. These HERV sequences will therefore enhance further studies on disease associations and differences between human haplotypes and primates and their role in the evolution of class I genes in the MHC. Received: 17 September 1998 / Accepted: 8 January 1999     

Phylogenomic Analysis of the α Proteasome Gene Family from Early-Diverging Eukaryotes

We employed a phylogenomic approach to study the evolution of α subunits of the proteasome gene family from early diverging eukaryotes. BLAST similarity searches of the Giardia lamblia genome identified all seven α proteasome genes characteristic of eukaryotes from the crown group. In addition, a PCR strategy for the amplification of multiple α subunit sequences generated single α proteasome products for representatives of the Kinetoplastida (Leishmania major), the Parabasalia (Trichomonas vaginalis), and the Microsporidia (Vairimorpha sp., Nosema sp., Endoreticulata sp., and Spraguea lophii). The kinetoplastid Trypanosoma cruzi and the eukaryote crown group Acanthamoeba castellanii yielded two distinct α proteasome genes each. The presence of seven distinct α proteasome genes in G. lamblia, one of the earliest-diverging eukaryotes, indicates that the α proteasome gene family evolved rapidly from a minimum of one gene in Archaea to seven or more in Eukarya. Results from the phylogenomic analysis are consistent with the idea that the Diplomonida (as represented by G. lamblia), the Kinetoplastida, the Parabasalia, and the Microsporidia diverged after the duplication events that originated the α proteasome gene family. A model for the early origin and evolution of the proteasome gene family is presented. Received: 14 February 2000 / Accepted: 14 August 2000     

Sequence evolution in and around the mitochondrial control region in birds

By cloning and sequencing 3.4 kilobases of snow goose mtDNA we found that the ND5 gene is followed by the genes for cytochrome b, tRNA^Thr, tRNA^Pro, ND6, tRNA^Glu, the control region, tRNA^Phe, and srRNA. This order is identical to that of chicken, quail, and duck mtDNA but differs from that of mammals and a frog (Xenopus). The mean extent of difference due to base substitution between goose and chicken is generally closer to the same comparison between rat and mouse but less than that between human and cow. For one of the nine regions compared (tRNA^Glu), the bird differences appear to be anomalous, possibly implicating altered functional constraints. Within the control region, several short sequences common to mammals are also conserved in the birds. Comparison of the goose control region with that of quail and chicken suggests that a sequence element with similarity to CSB-1 duplicated once prior to the divergence of goose and chicken and again on the lineage leading to chicken. Between goose (or duck) and chicken there are four times more transversions at the third positions of fourfold-degenerate codons in mitochondrial than in nuclear genes.Abbreviations CSB conserved sequence block - cytb cytochrome b - ND NADH dehydrogenase - srRNA small-subunit ribosomal RNA Deceased July 21, 1991 Correspondence to: T.W. Quinn at the University of Denver     

Genomic organization of the extracellular coding region of the human FGFR4 and FLT4 genes: Evolution of the genes encoding receptor tyrosine kinases with immunoglobulin-like domains

Receptor tyrosine kinases with five, seven, and three Ig-like domains in their extracellular region are grouped in subclasses IIIa, IIIb, and IIIc, respectively. Here, we describe the genomic organization of the extracellular coding region of the human FGFR4 (IIIc) and FLT4 (IIIb) genes and compare it to that of the human FGFR1(IIIc), KIT, and FMS (IIIa). The results show that while genes belonging to the same subclass have an identical exon/intron structure in their extracellular coding region—as they do in their intracellular coding region—genes of related subclasses only have a similar exon/intron structure. These results strongly support the hypothesis that the genes of the three subclasses evolved from a common ancestor by duplications involving entire genes, already in pieces. Hypotheses on the origin of introns and on the difference in the number of extracellular Ig-like domains in the three gene subclasses are discussed. Received: 19 August 1996 / Accepted: 2 January 1997     

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