Genomic structure and evolution of the human pepsinogen A multigene family

Summary A human cosmid library was screened with a pepsinogen A (PGA) cDNA probe, yielding 18 clones with (parts of) one, two or three PGA genes. By aligning these cosmids a restriction map of a PGA gene quadruplet was obtained in which the four genes are arranged in a highly ordered fashion in a head-to-tail orientation. Using the length in kilobases of the large polymorphic EcoRI fragment of the PGA genes, this quadruplet can be described as 15.0-12.0-12.0-16.6. An AvaII polymorphism allowed us to identify the two PGA haplotypes of the individual whose DNA had been cloned in the cosmid library to be a gene triplet and a gene quadruplet. By comparing the restriction maps of the central 12.0 genes in these multiplets to those of the flanking 15.0 and 16.6 genes, we postulate that these central genes arose from unequal but homologous crossing over between two 15.0–16.6 gene pairs. This hypothesis provides for the creation of a variety of haplotypes by additional cross overs and mutations. Southern blots of family and population material supports the existance of at least five common PGA haplotypes, including a single-gene haplotype, giving rise to a large number of different EcoRI patterns. The single PGA gene is probably the reciprocal crossing over product. Comparison between the DNA and protein polymorphisms suggests further micro-heterogeneity in the different PGA haplotypes.     

Assignment of the pepsinogen gene complex (PGA) to human chromosome region 11q13 by in situ hybridization

The genes coding for human pepsinogen (PGA3, PGA4, and PGA5) were assigned to chromosome region 11q13 by in situ hybridization. Previously we localized the PGA gene complex to a centromeric region of chromosome 11 (p11----q13) by Southern blot analysis of mouse-human somatic cell hybrids. Our in situ hybridization results confirm this assignment and further localize the genes to a smaller region on the long arm.     

Family and population studies on the human pepsinogen A multigene family

Summary Human pepsinogen A (PGA) displays highly polymorphic isozymogen patterns after polyacrylamide gel electrophoresis and activity staining. The patterns differ with respect to the presence and the relative intensity of the individual fractions. Family studies strongly suggest that these isozymogen patterns are encoded by allelic haplotypes, encompassing different numbers and types of PGA genes. In this paper, we confirm the essential features of this multigene model. We establish the relationship between the haplotypes and the corresponding isozymogen patterns by determination of the PGA polymorphism at both the DNA and the protein level in 117 Dutch individuals, 60 of whom were unrelated. The combination of HindIII and EcoRI restriction fragment length polymorphisms (RFLPs) has enabled us to define different haplotypes, which are shown to segregate within families. Most genes are characterized by their specific EcoRI fragments. The HindIII RFLP is in strong linkage disequilibrium with PGA genes showing strong expression of the relevant isozymogen. Although a general picture of the relationship between genotypes and phenotypes is emerging, there are exceptions, suggesting that rare haplotypes evolve by unique crossover events.     

Molecular cloning of a pair of human pepsinogen A genes which differ by a Glu→Lys mutation in the activation peptide

Summary Three human cosmid clones containing pepsinogen A (PGA) encoding sequences were isolated from a genomic bank derived from a single individual. One cosmid contains two PGA genes in tandem in a head-to-tail orientation, while the other two cosmids each contain a single PGA gene. The three cosmids were characterized by restriction mapping and sequence analysis (exons 1 and 2 and flanking regions). As judged from these data, three of the four PGA genes isolated appear to be nearly identical, but one of the tandem genes is clearly different from the other genes. The first exon of all four genes codes for the same amino acid sequence. However, in the second exon of one of the tandem genes we found a nucleotide substitution giving rise to a GluLys substitution of the 43rd amino acid residue of the activation peptide, leading to a charge difference of the corresponding isozymogens. The presence of two distinct PGA genes in the isolated gene pair conclusively proves the multigene structure of the PGA system. These genes might be responsible for at least part of the electrophoretic polymorphism at the protein level.     

Nucleotide sequence comparison of five human pepsinogen A (PGA) genes: evolution of the PGA multigene family

To unravel the genetic basis for the pepsinogen A (PGA) protein polymorphism, we have isolated and characterized a number of PGA genes, distinguishable by polymorphic EcoRI fragments of 12.0, 15.0, and 16.6 kb. Using a HindIII or AvaII polymorphism, we can discriminate between different 15.0 (15.0 and 15.0*) and 12.0 (12.0s and 12.0l) genes, respectively. The coding sequences of a 15.0 and a 16.6 gene were determined, together with considerable stretches of the 5'- and 3'-flanking regions and introns. The genes were demonstrated to encode Pg5 and Pg4, respectively. Because substitutions in codons 43 and 207 appeared to be critical in the determination of the encoded proteins, we sequenced only these regions in the two 12.0 genes and the 15.0* gene. On the basis of these partial sequences, we assume that these genes encode Pg3. In the evolutionary model of the PGA gene cluster presented here, the 12.0 genes arose by an unequal, but homologous crossover. The results of sequence analysis of the second intron of the 12.0s, 12.0l, 15.0, and 16.6 genes suggest that the two 12.0 genes have arisen from two different crossover events.     

Human pepsinogen A (PGA): an informative gene complex located at 11q13

Summary Human pepsinogen (PGA) exhibits extensive polymorphism that can be detected both at the protein and the DNA level. We describe here two restriction fragment length polymorphisms, EcoRI and BglII, which provide for the detection of three of the most common PGA haplotypes (A, B, and C) in the United States population. The relationship of these polymorphisms to each PGA haplotype was determined by analysis of DNA from individuals exhibiting the corresponding protein phenotypes and by analysis of a series of human × mouse somatic cell hybrids containing the individual chromosome 11 homologous from heterozygous individuals exhibiting the AB and AC protein phenotypes. The use of the BglII polymorphism in combination with previously described EcoRI polymorphism provides a very informative marker of 11q13.     

Characterization of the Bacillus subtilis ywsC gene, involved in gamma-polyglutamic acid production.

The genes required for gamma-polyglutamic acid (PGA) production were cloned from Bacillus subtilis IFO16449, a strain isolated from fermented soybeans. There were four open reading frames in the cloned 4.2-kb DNA fragment, and they were almost identical to those in the ywsC and ywtABC genes of B. subtlis 168. Northern blot analysis showed that the four genes constitute an operon. Three genes, ywsC, ywtA, and ywtB, were disrupted to determine which gene plays a central role in PGA biosynthesis. No PGA was produced in Delta ywsC and Delta ywtA strains, indicating that both of these genes are essential for PGA production. To clarify the function of the YwsC protein, histidine-tagged YwsC (YwsC-His) was produced in the Delta ywsC strain and purified from the lysozyme-treated lysate of the transformant by Ni-nitrilotriacetic acid affinity chromatography. Western blot analysis revealed that the YwsC-His protein consists of two subunits, the 44-kDa and 33-kDa proteins, which are encoded by in-phase overlapping in the ywsC gene. (14)C-labeled PGA was synthesized by the purified proteins from L-[(14)C]-glutamate in the presence of ATP and MnCl(2), through an acylphosphate intermediate, indicating that the ywsC gene encodes PGA synthetase (EC 6.3.2), a crucial enzyme in PGA biosynthesis.     

Localization of the photoreceptor gene ROM1 to human chromosome 11 and mouse chromosome 19: sublocalization to human 11q13 between PGA and PYGM.

Rom-1 is a retinal integral membrane protein that, together with the product of the human retinal degeneration slow gene (RDS), defines a photoreceptor-specific protein family. The gene for rom-1 (HGM symbol: ROM1) has been assigned to human chromosome 11 and mouse chromosome 19 by Southern blot analysis of somatic cell hybrid DNAs. ROM1 was regionally sublocalized to human 11p13-11q13 by using three mouse-human somatic cell hybrids; in situ hybridization refined the sublocalization to human 11q13. Analysis of somatic cell hybrids suggested that the most likely localization of ROM1 is in the approximately 2-cM interval between human PGA (human pepsinogen A) and PYGM (muscle glycogen phosphorylase). ROM1 appears to be a new member of a conserved syntenic group whose members include such genes as CD5, CD20, and OSBP (oxysterol-binding protein), on human chromosome 11 and mouse chromosome 19. Localization of the ROM1 gene will permit the examination of its linkage to hereditary retinopathies in man and mouse.     

Human pepsinogen C (progastricsin). Isolation of cDNA clones, localization to chromosome 6, and sequence homology with pepsinogen A

The entire pepsinogen C (PGC) coding sequence was determined by analysis of a series of five overlapping cDNA clones identified in a library constructed from human gastric mucosa poly(A+) RNA. A partial cDNA clone was initially identified using a 256-fold degenerate oligonucleotide probe for amino acid residues 4-12 of pepsin C, and subsequently 4 additional clones were identified upon rescreening with a probe complementary to the 5' region of the original cDNA clone. Northern analysis of gastric mucosa poly(A+) RNA with a PGC cDNA probe revealed an mRNA 1.5-kilobase species that was indistinguishable from that detected with a human pepsinogen A (PGA) cDNA probe. In contrast, the PGC and PGA cDNA probes detected distinct genomic restriction fragments indicating there was no detectable cross-hybridization under high stringency conditions. The PGC gene was localized to human chromosome 6 by analysis of a panel of human x mouse somatic cell hybrids. The regions containing the active site aspartyl groups of PGC are conserved in relationship to several other aspartic proteinases. We propose that the absence of detectable immunologic cross-reactivity between the two groups of human pepsinogens, A and C, results from divergent evolution of sequences located on the surface of the zymogens in contrast to the strongly conserved active site regions located within the binding cleft of the enzymes that are inaccessible for antigenic recognition.     

Haplotype diversity across 100 candidate genes for inflammation, lipid metabolism, and blood pressure regulation in two populations

Recent studies have suggested that a significant fraction of the human genome is contained in blocks of strong linkage disequilibrium, ranging from ~5 to >100 kb in length, and that within these blocks a few common haplotypes may account for >90% of the observed haplotypes. Furthermore, previous studies have suggested that common haplotypes in candidate genes are generally shared across populations and represent the majority of chromosomes in each population. The conclusions drawn from these preliminary studies, however, are based on an incomplete knowledge of the variation in the regions examined. To bridge this gap in knowledge, we have completely resequenced 100 candidate genes in a population of African descent and one of European descent. Although these genes have been well studied because of their medical importance, we demonstrate that a large amount of sequence variation has not yet been described. We also report that the average number of inferred haplotypes per gene, when complete data is used, is higher than in previous reports and that the number and proportion of all haplotypes represented by common haplotypes per gene is variable. Furthermore, we demonstrate that haplotypes shared between the two populations constitute only a fraction of the total number of haplotypes observed and that these shared haplotypes represent fewer of the African-descent chromosomes than was expected from previous studies. Finally, we show that restricting variation discovery to coding regions does not adequately describe all common haplotypes or the true haplotype block structure observed when all common variation is used to infer haplotypes. These data, derived from complete knowledge of genetic variation in these genes, suggest that the haplotype architecture of candidate genes across the human genome is more complex than previously suggested, with important implications for candidate gene and genomewide association studies.     

Human pepsinogen C (progastricsin) polymorphism: evidence for a single locus located at 6p21.1-pter

A series of six clones containing the entire human pepsinogen C gene (PGC) was identified in a cosmid vector library by using cDNA and oligonucleotide probes. The 10.7-kb PGC gene includes nine exons and exhibits a high degree of sequence identity (60%) with the functionally related pepsinogen A genes. The predicted amino acid sequence was identical with the partial amino-terminal and carboxyl-terminal sequences of purified pepsinogen C. An informative restriction fragment length polymorphism was detected with several restriction enzymes and involved an insertion or deletion of 100 bp of intron sequence located between exons 7 and 8. Evidence that there is only a single PGC gene in humans is presented. The PGC gene and the prolactin gene were regionally localized to 6p21.1-pter by analysis of mouse X human somatic cell hybrids.     

Evolutionary signatures of common human cis-regulatory haplotypes

Variation in gene expression may give rise to a significant fraction of inter-individual phenotypic variation. Studies searching for the underlying genetic controls for such variation have been conducted in model organisms and humans in recent years. In our previous effort of assessing conserved underlying haplotype patterns across ethnic populations, we constructed common haplotypes using SNPs having conserved linkage disequilibrium (LD) across ethnic populations. These common haplotypes cluster into a simple evolutionary structure based on their frequencies, defining only up to three conserved clusters termed 'haplotype frameworks'. One intriguing preliminary finding was that a significant portion of reported variants strongly associated with cis-regulation tags these globally conserved haplotype frameworks. Here we expand the investigation by collecting genes showing stringently determined cis-association between genotypes and expression phenotypes from major studies. We conducted phylogenetic analysis of current major haplotypes along with the corresponding haplotypes derived from chimpanzee reference sequences. Our analysis reveals that, for the vast majority of such cis-regulatory genes, the tagging SNPs showing the strongest association also tag the haplotype lineages directly separated from ancestry, inferred from either chimpanzee reference sequences or the allele frequency-derived haplotype frameworks, suggesting that the differentially expressed phenotypes were evolved relatively early in human history. Such evolutionary signatures provide keys for a more effective identification of globally-conserved candidate regulatory haplotypes across human genes in future epidemiologic and pharmacogenetic studies.     

Identification of pepsinogen gene in the genome of stomachless fish, Takifugu rubripes

Pepsinogen is a precursor of pepsin, a gastric specific protease belonging to the aspartic proteinase family. In teleosts, several species, such as zebrafish and puffer, have independently lost gastric glands. So whether puffer have pepsinogen gene or not is an interesting issue. A search of GSS database for pufferfish, Takifugu rubripes, revealed five different aspartic proteinase genes in its genome. One of them (pPep) has typical pepsinogen structure and belongs to the fish pepsinogen cluster by phylogenic analysis. The pPep antisense probe hybridized to the gastric glands of flounder stomach. Therefore, we concluded that pPep is a pufferfish pepsinogen. The pufferfish pepsinogen mRNA was not expressed in the digestive organs but specifically in the skin. We speculated that while Tetraodontiformes evolutionarily lost gastric glands, pufferfish pepsinogen acquired an alternative function to food digestion.     

Gene conversion rapidly generates major histocompatibility complex diversity in recently founded bird populations

Population bottlenecks can restrict variation at functional genes, reducing the ability of populations to adapt to new and changing environments. Understanding how populations generate adaptive genetic variation following bottlenecks is therefore central to evolutionary biology. Genes of the major histocompatibility complex (MHC) are ideal models for studying adaptive genetic variation due to their central role in pathogen recognition. While de novo MHC sequence variation is generated by point mutation, gene conversion can generate new haplotypes by transferring sections of DNA within and across duplicated MHC loci. However, the extent to which gene conversion generates new MHC haplotypes in wild populations is poorly understood. We developed a 454 sequencing protocol to screen MHC class I exon 3 variation across all 13 island populations of Berthelot's pipit (Anthus berthelotii). We reveal that just 11-15 MHC haplotypes were retained when the Berthelot's pipit dispersed across its island range in the North Atlantic ca. 75,000 years ago. Since then, at least 26 new haplotypes have been generated in situ across populations. We show that most of these haplotypes were generated by gene conversion across divergent lineages, and that the rate of gene conversion exceeded that of point mutation by an order of magnitude. Gene conversion resulted in significantly more changes at nucleotide sites directly involved with pathogen recognition, indicating selection for functional variants. We suggest that the creation of new variants by gene conversion is the predominant mechanism generating MHC variation in genetically depauperate populations, thus allowing them to respond to pathogenic challenges.     

Haplotypes at ATM identify coding-sequence variation and indicate a region of extensive linkage disequilibrium

Genetic variation in the human population may lead to functional variants of genes that contribute to risk for common chronic diseases such as cancer. In an effort to detect such possible predisposing variants, we constructed haplotypes for a candidate gene and tested their efficacy in association studies. We developed haplotypes consisting of 14 biallelic neutral-sequence variants that span 142 kb of the ATM locus. ATM is the gene responsible for the autosomal recessive disease ataxia-telangiectasia (AT). These ATM noncoding single-nucleotide polymorphisms (SNPs) were genotyped in nine CEPH families (89 individuals) and in 260 DNA samples from four different ethnic origins. Analysis of these data with an expectation-maximization algorithm revealed 22 haplotypes at this locus, with three major haplotypes having frequencies > or = .10. Tests for recombination and linkage disequilibrium (LD) show reduced recombination and extensive LD at the ATM locus, in all four ethnic groups studied. The most striking example was found in the study population of European ancestry, in which no evidence for recombination could be discerned. The potential of ATM haplotypes for detection of genetic variants through association studies was tested by analysis of 84 individuals carrying one of three ATM coding SNPs. Each coding SNP was detected by association with an ATM haplotype. We demonstrate that association studies with haplotypes for candidate genes have significant potential for the detection of genetic backgrounds that contribute to disease.     

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.     

Different patterns of evolution in the centromeric and telomeric regions of group A and B haplotypes of the human killer cell Ig-like receptor locus

The fast evolving human KIR gene family encodes variable lymphocyte receptors specific for polymorphic HLA class I determinants. Nucleotide sequences for 24 representative human KIR haplotypes were determined. With three previously defined haplotypes, this gave a set of 12 group A and 15 group B haplotypes for assessment of KIR variation. The seven gene-content haplotypes are all combinations of four centromeric and two telomeric motifs. 2DL5, 2DS5 and 2DS3 can be present in centromeric and telomeric locations. With one exception, haplotypes having identical gene content differed in their combinations of KIR alleles. Sequence diversity varied between haplotype groups and between centromeric and telomeric halves of the KIR locus. The most variable A haplotype genes are in the telomeric half, whereas the most variable genes characterizing B haplotypes are in the centromeric half. Of the highly polymorphic genes, only the 3DL3 framework gene exhibits a similar diversity when carried by A and B haplotypes. Phylogenetic analysis and divergence time estimates, point to the centromeric gene-content motifs that distinguish A and B haplotypes having emerged ~6 million years ago, contemporaneously with the separation of human and chimpanzee ancestors. In contrast, the telomeric motifs that distinguish A and B haplotypes emerged more recently, ~1.7 million years ago, before the emergence of Homo sapiens. Thus the centromeric and telomeric motifs that typify A and B haplotypes have likely been present throughout human evolution. The results suggest the common ancestor of A and B haplotypes combined a B-like centromeric region with an A-like telomeric region.