Expression, linkage, and polymorphism of MHC-related genes in rainbow trout, Oncorhynchus mykiss.

The architecture of the MHC in teleost fish, which display a lack of linkage between class I and II genes, differs from all other vertebrates. Because rainbow trout have been examined for a variety of immunologically relevant genes, they present a good teleost model for examining both the expression and organization of MHC-related genes. Full-length cDNA and partial gDNA clones for proteasome delta, low molecular mass polypeptide (LMP) 2, TAP1, TAP2A, TAP2B, class Ia, and class IIB were isolated for this study. Aside from the expected polymorphisms associated with class I genes, LMP2 and TAP2 are polygenic. More specifically, we found a unique lineage of LMP2 (LMP2/delta) that shares identity to both LMP2 and delta but is expressed like the standard LMP2. Additionally, two very different TAP2 loci were found, one of which encodes polymorphic alleles. In general, the class I pathway genes are expressed in most tissues, with highest levels in lymphoid tissue. We then analyzed the basic genomic organization of the trout MHC in an isogenic backcross. The main class Ia region does not cosegregate with the class IIB locus, but LMP2, LMP2/delta, TAP1A, and TAP2B are linked to the class Ia locus. Interestingly, TAP2A (second TAP2 locus) is a unique lineage in sequence composition that appears not to be linked to this cluster or to class IIB. These results support and extend the recent findings of nonlinkage between class I and II in a different teleost order (cyprinids), suggesting that this unique arrangement is common to all teleosts.     

Proteasome, transporter associated with antigen processing, and class I genes in the nurse shark Ginglymostoma cirratum: evidence for a stable class I region and MHC haplotype lineages.

Cartilaginous fish (e.g., sharks) are derived from the oldest vertebrate ancestor having an adaptive immune system, and thus are key models for examining MHC evolution. Previously, family studies in two shark species showed that classical class I (UAA) and class II genes are genetically linked. In this study, we show that proteasome genes LMP2 and LMP7, shark-specific LMP7-like, and the TAP1/2 genes are linked to class I/II. Functional LMP7 and LMP7-like genes, as well as multiple LMP2 genes or gene fragments, are found only in some sharks, suggesting that different sets of peptides might be generated depending upon inherited MHC haplotypes. Cosmid clones bearing the MHC-linked classical class I genes were isolated and shown to contain proteasome gene fragments. A non-MHC-linked LMP7 gene also was identified on another cosmid, but only two exons of this gene were detected, closely linked to a class I pseudogene (UAA-NC2); this region probably resulted from a recent duplication and translocation from the functional MHC. Tight linkage of proteasome and class I genes, in comparison with gene organizations of other vertebrates, suggests a primordial MHC organization. Another nonclassical class I gene (UAA-NC1) was detected that is linked neither to MHC nor to UAA-NC2; its high level of sequence similarity to UAA suggests that UAA-NC1 also was recently derived from UAA and translocated from MHC. These data further support the principle of a primordial class I region with few class I genes. Finally, multiple paternities in one family were demonstrated, with potential segregation distortions.     

Sequence organization and matrix attachment regions of the human serine protease inhibitor gene cluster at 14q32.1

The human serine protease inhibitor (serpin) gene cluster at 14q32.1 is a useful model system for studying the regulation of gene activity and chromatin structure. We demonstrated previously that the six known serpin genes in this region were organized into two subclusters of three genes each that occupied ~370 kb of DNA. To more fully understand the genomic organization of this region, we annotated a 1-Mb sequence contig from data from the Genoscope sequencing consortium ( ). We report that 11 different serpin genes reside within the 14q32.1 cluster, including two novel 1-antiproteinase-like gene sequences, a kallistatin-like sequence, and two recently identified serpins that had not been mapped previously to 14q32.1. The genomic regions proximal and distal to the serpin cluster contain a variety of unrelated gene sequences of diverse function. To gain insight into the chromatin organization of the region, sequences with putative nuclear matrix-binding potential were identified by using the MAR-Wiz algorithm, and these MAR-Wiz candidate sequences were tested for nuclear matrix-binding activity in vitro. Several differences between the MAR-Wiz predictions and the results of biochemical tests were observed. The genomic organization of the serpin gene cluster is discussed. These authors (Stephanie J. Namciu and Richard D. Friedman) contributed equally to this work.     

Isolation of a 97-kb minimal essential MHC B locus from a new reverse-4D BAC library of the golden pheasant

The bacterial artificial chromosome (BAC) system is widely used in isolation of large genomic fragments of interest. Construction of a routine BAC library requires several months for picking clones and arraying BACs into superpools in order to employ 4D-PCR to screen positive BACs, which might be time-consuming and laborious. The major histocompatibility complex (MHC) is a cluster of genes involved in the vertebrate immune system, and the classical avian MHC-B locus is a minimal essential one, occupying a 100-kb genomic region. In this study, we constructed a more effective reverse-4D BAC library for the golden pheasant, which first creates sub-libraries and then only picks clones of positive sub-libraries, and identified several MHC clones within thirty days. The full sequencing of a 97-kb reverse-4D BAC demonstrated that the golden pheasant MHC-B locus contained 20 genes and showed good synteny with that of the chicken. The notable differences between these two species were the numbers of class II B loci and NK genes and the inversions of the TAPBP gene and the TAP1-TAP2 region. Furthermore, the inverse TAP2-TAP1 was unique in the golden pheasant in comparison with that of chicken, turkey, and quail. The newly defined genomic structure of the golden pheasant MHC will give an insight into the evolutionary history of the avian MHC.     

Concerted evolution of human amylase genes.

Cosmid clones containing 250 kilobases of genomic DNA from the human amylase gene cluster have been isolated. These clones contain seven distinct amylase genes which appear to comprise the complete multigene family. By sequence comparison with the cDNAs, we have identified two pancreatic amylase genes and three salivary amylase genes. Two truncated pseudogenes were also recovered. Intergenic distances of 17 to 22 kilobases separate the amylase gene copies. Within the past 10 million years, duplications, gene conversions, and unequal crossover events have resulted in a very high level of sequence similarity among human amylase gene copies. To identify sequence elements involved in tissue-specific expression and hormonal regulation, the promoter regions of the human amylase genes were sequenced and compared with those of the corresponding mouse genes. The promoters of the human and mouse pancreatic amylase genes are highly homologous between nucleotide -160 and the cap site. Two sequence elements thought to influence pancreas-specific expression of the rodent genes are present in the human genes. In contrast, similarity in the 5' flanking sequences of the salivary amylase genes is limited to several short sequence elements whose positions and orientations differ in the two species. Some of these sequence elements are also associated with other parotid-specific genes and may be involved in their tissue-specific expression. A glucocorticoid response element and a general enhancer element are closely associated in several of the amylase promoters.     

Analysis of a 26-kb region linked to the Mhc in zebrafish: genomic organization of the proteasome component beta/transporter associated with antigen processing-2 gene cluster and identification of five new proteasome beta subunit genes.

Sequencing of zebrafish (Danio rerio) bacterial artificial chromosome and P1 artificial chromosome genomic clone fragments and of cDNA clones has led to the identification of five new loci coding for beta subunits of proteasomes (PSMB). Together with the four genes identified previously, nine PSMB genes have now been defined in the zebrafish. Six of the nine genes reside in the zebrafish MHC (Mhc) class I region, four of them reside in a single cluster closely associated with TAP2 on a 26-kb long genomic fragment, and two reside at some distance from the fragment. In addition to homologues of the human genes PSMB5 through PSMB9, two new genes, PSMB11 and PSMB12, have been found for which there are no known corresponding genes in humans. The new genes reside in the PSMB cluster in the Mhc. Homology and promoter region analysis suggest that the Mhc-associated genes might be inducible by IFN-gamma. The zebrafish class I region contains representatives of three phylogenetically distinguishable groups of PSMB genes, X, Y, and Z. It is proposed that these genes were present in the ancestral PSMB region before Mhc class I genes became associated with it.     

Sequence, structure, and evolution of a complete human olfactory receptor gene cluster

The olfactory receptor (OR) gene cluster on human chromosome 17p13.3 was subjected to mixed shotgun automated DNA sequencing. The resulting 412 kb of genomic sequence include 17 OR coding regions, 6 of which are pseudogenes. Six of the coding regions were discovered only upon genomic sequencing, while the others were previously reported as partial sequences. A comparison of DNA sequences in the vicinity of the OR coding regions revealed a common gene structure with an intronless coding region and at least one upstream noncoding exon. Potential gene control regions including specific pyrimidine:purine tracts and Olf-1 sites have been identified. One of the pseudogenes apparently has evolved into a CpG island. Four extensive CpG islands can be discerned within the cluster, not coupled to specific OR genes. The cluster is flanked at its telomeric end by an unidentified open reading frame (C17orf2) with no significant similarity to any known protein. A high proportion of the cluster sequence (about 60%) belongs to various families of interspersed repetitive elements, with a clear predominance of LINE repeats. The OR genes in the cluster belong to two families and seven subfamilies, which show a relatively high degree of intermixing along the cluster, in seemingly random orientations. This genomic organization may be best accounted for by a complex series of evolutionary events.     

Sequences, annotation and single nucleotide polymorphism of the major histocompatibility complex in the domestic cat

Two sequences of major histocompatibility complex (MHC) regions in the domestic cat, 2.976 and 0.362 Mbps, which were separated by an ancient chromosome break (55-80 MYA) and followed by a chromosomal inversion were annotated in detail. Gene annotation of this MHC was completed and identified 183 possible coding regions, 147 human homologues, possible functional genes and 36 pseudo/unidentified genes) by GENSCAN and BLASTN, BLASTP RepeatMasker programs. The first region spans 2.976 Mbp sequence, which encodes six classical class II antigens (three DRA and three DRB antigens) lacking the functional DP, DQ regions, nine antigen processing molecules (DOA/DOB, DMA/DMB, TAPASIN, and LMP2/LMP7,TAP1/TAP2), 52 class III genes, nineteen class I genes/gene fragments (FLAI-A to FLAI-S). Three class I genes (FLAI-H, I-K, I-E) may encode functional classical class I antigens based on deduced amino acid sequence and promoter structure. The second region spans 0.362 Mbp sequence encoding no class I genes and 18 cross-species conserved genes, excluding class I, II and their functionally related/associated genes, namely framework genes, including three olfactory receptor genes. One previously identified feline endogenous retrovirus, a baboon retrovirus derived sequence (ECE1) and two new endogenous retrovirus sequences, similar to brown bat endogenous retrovirus (FERVmlu1, FERVmlu2) were found within a 140 Kbp interval in the middle of class I region. MHC SNPs were examined based on comparisons of this BAC sequence and MHC homozygous 1.9x WGS sequences and found that 11,654 SNPs in 2.84 Mbp (0.00411 SNP per bp), which is 2.4 times higher rate than average heterozygous region in the WGS (0.0017 SNP per bp genome), and slightly higher than the SNP rate observed in human MHC (0.00337 SNP per bp).     

Evolution of human immunoglobulin kappa J region genes

Immunoglobulin kappa chain variable region genes are assembled from two discontinuous DNA segments, a V and a J gene. The J region genes, in addition to encoding amino acid positions 96-108 of the kappa polypeptide chain, also provide sequences required for both DNA and RNA splicing reactions. For purposes of evolutionary comparison and to establish the complexity of the kappa J region locus in man, we have determined an approximately 3000 basepair nucleotide sequence in a cloned human DNA fragment that encodes the germline distinct J region segments. Significant blocks of homology have been tightly maintained between this region and an analogous segment of the mouse genome. In particular, the short sequences, GGTTTTTGT and CACTGTG, thought to be involved in V-J recombination, are the most highly conserved regions (97% homology). In addition, from heteroduplex data and computer analysis of the nucleotide sequences, it is clear that the mouse J3 sequence, a pseudogene, is not present in the human cluster. This can be explained by a duplication event in the mouse J region gene cluster that may have been the result of unequal crossing over between homologous chromosomes.     

Molecular characterization, expression and mapping of porcine LMP2 and MECL-1 genes.

Low molecular mass polypeptides 2 (LMP2) and Multicatalytic endopeptidase complex subunit (MECL-1) are two of three catalytic beta-type subunits of the 20s proteasome and upon interferon gamma-induction. LMP2 is critical for the production of major histocompatibility complex (MHC) class I ligand and T-lymphocytes. The LMP2 gene is located in the MHC region, but MECL-1 is located outside the MHC region. They are involved in the antigen presentation and are important candidate genes for an initial exploration of relationships between the antigen processing genes and disease resistance. In this report, the porcine LMP2 and MECL-1 cDNA were cloned and A 5099 bp LMP2 genomic DNA structure was identified, then two single nucleotide polymorphisms were detected in the exon2 and exon5 of LMP2 gene in 367 individuals. The LMP2 and MECL-1 genes putative protein included 219,274 amino acids, respectively. Alignment and phylogenetic of predicted porcine LMP2 and MECL-1 amino acid sequence with their homologies were analyzed. Tissues expression of LMP2 and MECL-1 mRNA were observed by real time quantitative PCR (Q-PCR) method, the results revealed MECL-1 expressed widely in all tissues, but LMP2 was not detected in muscle. The porcine MECL-1 gene was mapped to chromosome 6, closely linked to microsatellite SW1108 (LOD = 4.09, 84cR) by radiation hybrid panel.     

Concerted evolution within a trypsin gene cluster in Drosophila.

A cluster of four trypsin genes has previously been localized to cytological position 47D-F of the Drosophila melanogaster genome. One of these genes had been sequenced, and the presence of the other three genes was identified by cross-hybridization. Here, we present the DNA sequence of the entire genomic region encoding these four trypsin genes. In addition to the four previously inferred genes, we have identified a fifth trypsin-coding sequence located within this gene cluster. This new gene shows a high degree of sequence divergence (more than 30%) from the other four genes, although it retains all of the functional motifs that are characteristic of trypsin-coding sequences. In order to trace the molecular evolution of this gene cluster, we isolated and sequenced the homologous 7-kb region from the closely related species Drosophila erecta. A comparison of the DNA sequences between the two species provides strong evidence for the concerted evolution of some members of this gene family. Two genes within the cluster are evolving in concert, while a third gene appears to be evolving independently. The remaining two genes show an intermediate pattern of evolution. We propose a simple model, involving chromosome looping and gene conversion, to explain the relatively complex patterns of molecular evolution within this gene cluster.     

Molecular mapping of a recombination hotspot located in the second intron of the human TAP2 locus.

Recombination across the HLA class II region is not randomly distributed, as indicated by both strong linkage disequilibrium within the 100 kb encompassing the DRB1-DQA1-DQB1 loci and complete equilibrium between TAP1 and TAP2, the closest variant sites of which are < 15 kb. In an attempt to explain these observations, 39 novel polymorphic markers in a region encompassing the TAP, LMP, and DOB genes were used to delineate the site of crossover in 11 class II recombinant chromosomes. SSCP demonstrated that two recombination events occurred within an 850-bp interval in the second intron of TAP2, which separates the variant sites of TAP1 and TAP2. These data indicate the presence of a recombination hotspot, the first to be identified from the analysis of familial transmission in the human major histocompatibility complex. The region of crossover was cloned and sequenced from one of the recombinants, further defining the crossover site to a 138-bp segment nested within the 850-bp region. This represents the most precisely defined region of recombination in the human genome.     

Multispecies comparative analysis of a mammalian-specific genomic domain encoding secretory proteins

The mammalian-specific casein gene cluster comprises 3 or 4 evolutionarily related genes and 1 physically linked gene with a functional association. To gain a better understanding of the mechanisms regulating the entire casein cluster at the genomic level we initiated a multispecies comparative sequence analysis. Despite the high level of divergence at the coding level, these studies have identified uncharacterized family members within two species and the presence at orthologous positions of previously uncharacterized genes. Also the previous suggestion that the histatin/statherin gene family, located in this region, was primate specific was ruled out. All 11 genes identified in this region appear to encode secretory proteins. Conservation of a number of noncoding regions was observed; one coincides with an element previously suggested to be important for beta-casein gene expression in human and cow. The conserved regions might have biological importance for the regulation of genes in this genomic "neighborhood."     

Analysis of the genomic structure of the porcine CD1 gene cluster

CD1 is an MHC class I-like protein that presents lipid antigens to T cell receptors. We determined 470,187 bp of the genomic sequence encompassing the region encoding porcine CD1 genes. We identified 16 genes in this region and newly identified CD1A2, CD1B, CD1C, CD1D, and CD1E. Porcine CD1 genes were located in clusters between KIRREL and olfactory receptor (OR) genes, as observed in humans, although they were divided into two regions by a region encoding OR genes. Comparison of the genomic sequences of CD1 gene loci in pigs with other mammals showed that separation of the CD1 gene cluster by ORs was observed only in pigs. CD1A duplication in the porcine genome was estimated to have occurred after the divergence of the human and porcine. This analysis of the genomic sequence of the porcine CD1 family will contribute to our understanding of the evolution of mammalian CD1 genes.