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.     

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.     

Theoretical population genetics of repeated genes forming a multigene family

The evolution of repeated genes forming a multigene family in a finite population is studied with special reference to the probability of gene identity, i.e., the identity probability of two gene units chosen from the gene family. This quantity is called clonality and is defined as the sum of squares of the frequencies of gene lineages in the family. The multigene family is undergoing continuous unequal somatic crossing over, ordinary interchromosomal crossing over, mutation and random frequency drift. Two measures of clonality are used: clonality within one chromosome and that between two different chromosomes. The equilibrium properties of the means, the variances and the covariance of the two measures of clonality are investigated by using the diffusion equation method under the assumption of constant number of gene units in the multigene family. Some models of natural selection based on clonality are considered. The possible significance of the variance and covariance of clonality among the chromosomes on the adaptive differentiation of gene families such as those producing antibodies is discussed.     

Isolation of a multigene family containing human alpha-tubulin sequences

The boundaries of the origin of polyoma DNA replication have been analyzed using a set of deletion mutants. The majority of these had small deletions, 5 to 30 base-pairs in size, which together removed most of the non-translated sequences of the genome. The phenotype of the mutants was characterized by analysis of infectivity, transforming ability and DNA synthesis. All mutants with reduced or abolished infectivity had corresponding defects of viral DNA synthesis. The effect of the deletion was cis-acting, since the replication of the mutants was not stimulated by the presence of wild-type DNA. Deletions causing a reduction of DNA synthesis were found at two sites. The first at the 32 base-pair inverted repeat sequence and the neighbouring A · T tract previously implicated in the initiation of DNA synthesis, and the second close to the late genes. The two sites were separated by at least 60 base-pairs of non-essential DNA. Only one mutant with a deletion at the second site was unable to express early gene functions.The mutants were constructed by linearization, shortening and recircularization of polyoma DNA inserted into the plasmid pBR322. The mutagenesis was directed at restriction endonuclease BglI or PvuII cleavage sites. The BglI-directed mutagenesis was focussed to polyoma DNA by using as a vector a derivative of pBR322 resistant to cleavage by BglI.     

Evolution of the regulators of G-protein signaling multigene family in mouse and human

The regulators of G-protein signaling (RGS) proteins are important regulatory and structural components of G-protein coupled receptor complexes. RGS proteins are GTPase activating proteins (GAPs) of Gi-and Gq-class Galpha proteins, and thereby accelerate signaling kinetics and termination. Here, we mapped the chromosomal positions of all 21 Rgs genes in mouse, and determined human RGS gene structures using genomic sequence from partially assembled bacterial artificial chromosomes (BACs) and Celera fragments. In mice and humans, 18 of 21 RGS genes are either tandemly duplicated or tightly linked to genes encoding other components of G-protein signaling pathways, including Galpha, Ggamma, receptors (GPCR), and receptor kinases (GPRK). A phylogenetic tree revealed seven RGS gene subfamilies in the yeast and metazoan genomes that have been sequenced. We propose that similar systematic analyses of all multigene families from human and other mammalian genomes will help complete the assembly and annotation of the human genome sequence.     

Interpretation of polymorphic DNA patterns in the human alpha-amylase multigene family.

Previous molecular studies have clearly shown that the human amylase locus has a very complicated structure. Multiple salivary and pancreatic amylase genes are present on haplotypes with variable numbers of genes. To study the population heterogeneity, human genomic DNA from family members and random individuals was digested with a number of different restriction enzymes and hybridized with probes representing various parts of the human pancreatic amylase cDNA. The complex patterns obtained were, in most cases, compatible with predictions from the restriction enzyme maps of cloned human amylase genes. With some enzymes deviations from the predicted intensities of the bands associated with the pancreatic amylase gene AMY2A were observed. These findings can be explained by unequal homologous crossovers between AMY2A and AMY1A, resulting in haplotypes with one gene less or one gene more than the haplotypes described thus far. Moreover, a very complicated TaqI polymorphism was found that can be explained by homologous crossovers between different salivary amylase genes. Because some salivary amylase genes have an inverted orientation with respect to the others, these data provide evidence for the occurrence of intrachromosomal, homologous crossovers, as proposed by us previously (P. C. Groot et al., 1990, Genomics 8: 97-105).     

A multigene family encoding R-SNAREs in the ciliate Paramecium tetraurelia

SNARE proteins (soluble N-ethylmaleimide-sensitive factor attachment protein receptors) mediate membrane interactions and are conventionally divided into Q-SNAREs and R-SNAREs according to the possession of a glutamine or arginine residue at the core of their SNARE domain. Here, we describe a set of R-SNAREs from the ciliate Paramecium tetraurelia consisting of seven families encoded by 12 genes that are expressed simultaneously. The complexity of the endomembrane system in Paramecium can explain this high number of genes. All P. tetraurelia synaptobrevins (PtSybs) possess a SNARE domain and show homology to the Longin family of R-SNAREs such as Ykt6, Sec22 and tetanus toxin-insensitive VAMP (TI-VAMP). We localized four exemplary PtSyb subfamilies with GFP constructs and antibodies on the light and electron microscopic level. PtSyb1-1, PtSyb1-2 and PtSyb3-1 were found in the endoplasmic reticulum, whereas PtSyb2 is localized exclusively in the contractile vacuole complex. PtSyb6 was found cytosolic but also resides in regularly arranged structures at the cell cortex (parasomal sacs), the cytoproct and oral apparatus, probably representing endocytotic compartments. With gene silencing, we showed that the R-SNARE of the contractile vacuole complex, PtSyb2, functions to maintain structural integrity as well as functionality of the osmoregulatory system but also affects cell division.     

A human RNA polymerase II subunit is encoded by a recently generated multigene family

Background  

The sequences encoding the yeast RNA polymerase II (RPB) subunits are single copy genes.     

A human RNA polymerase II subunit is encoded by a recently generated multigene family

Background  

Some origins in eukaryotic chromosomes fire more frequently than others. In the fission yeast, Schizosaccharomyces pombe, the relative firing frequencies of the three origins clustered 4-8 kbp upstream of the ura4 gene are controlled by a replication enhancer - an element that stimulates nearby origins in a relatively position-and orientation-independent fashion. The important sequence motifs within this enhancer were not previously localized.     

A multigene framework for polychaete phylogenetic studies

Producing a robust phylogenetic reconstruction for Polychaeta using either morphological or molecular data sets has proven very difficult. There remain many conflicts between morphological analyses and hypotheses based on DNA data, the latter principally derived from 18S rRNA sequences. For the present study a data set covering a broad range of polychaete diversity was assembled, including 38 new sequences from 21 species. Besides available 18S rRNA data, five additional gene segments were examined: the D1 and D9-10 expansion regions of 28S rRNA, histone H3, snU2 RNA and cytochrome c oxidase subunit I. Maximum parsimony, maximum likelihood and Bayesian analyses were conducted.Annelida and Mollusca were reciprocally monophyletic in maximum likelihood analyses, but Polychaeta included a cephalopod in maximum parsimony analyses, and a patellogastropod in Bayesian analyses. When rooted on the Mollusca, optimal topologies from maximum likelihood analyses showed a recognisable basal group of taxa, including Oweniidae, Chaetopteridae and Amphinomidae. The six studied phyllodocidan families plus Orbiniidae (as the sister group of the scale-worms) formed the next most basal group. All analyses support the inclusion of Echiura, Clitellata and Siboglinidae within polychaetes. Bayesian analyses show Echiura as the sister group of Capitellidae, in agreement with previous 18S rRNA results, In contrast, Echiura formed the sister group to Trichobranchidae in maximum likelihood and maximum parsimony analyses.Supra-familial groupings consistently recovered within Polychaeta in the analyses are: (i) Terebellida without Ampharetidae; (ii) Scolecida (excepting Orbiniidae); (iii) Eunicidae, Lumbrineridae and Clitellata; and (iv) “Cirratuliformia” (including Sternaspidae) plus Sabellidae, Serpulidae and Spionidae.     

Structure and evolution of the apolipoprotein multigene family

We present the complementary DNA and deduced amino acid sequence of rat apolipoprotein A-II (apoA-II), and the results of a detailed statistical analysis of the nucleotide and amino acid sequences of all the apolipoprotein gene sequences published to date: namely, those of human and rat apoA-I, apoA-II and apoE, rat apoA-IV, and human apoC-I, C-II and C-III. Our results indicate that the apolipoprotein genes have very similar genomic structures, each having a total of three introns at the same locations. Using the exon/intron junctions as reference points, we have obtained an alignment of the coding regions of all the genes studied. It appears that the mature peptide regions of these genes are almost completely made up of tandem repeats of 11 codons. The part of mature peptide region encoded by exon 3 contains a common block of 33 codons, whereas the part encoded by exon 4 contains a much more variable number of internal repeats of 11 codons. These genes have apparently evolved from a primordial gene through multiple partial (internal) and complete gene duplications. On the basis of the degree of homology of the various sequences, and the pattern of the internal repeats in these genes, we propose an evolutionary tree for the apolipoprotein genes and give rough estimates of the divergence times between these genes. Our results show that apoA-II has evolved extremely rapidly and that apoA-I and apoE also have evolved at high rates but some regions are better conserved than the others. The rate of evolution of individual regions seems to be related to the stringency of their functional requirements.     

Adenylate deaminase. A multigene family in humans and rats

Multiple AMP deaminase (AMP-D) isoforms have been found in vertebrates, and tissue-specific inherited deficiencies of AMP-D have been described in two unrelated clinical syndromes suggesting there may be more than one AMP-D gene in higher eukaryotes. Using a newly isolated cDNA cloned from an adult rat brain library and a previously reported cDNA cloned from adult rat skeletal muscle, two linked AMP-D genes have been identified in rat and man. ampd1 is expressed at high levels in skeletal muscle of the adult rat. ampd2 is the predominant gene expressed in non-muscle tissues and smooth muscle of the adult rat, and it is also the predominant gene expressed in embryonic muscle and undifferentiated myoblasts. Both genes are expressed in cardiac muscle of the adult rat. The peptides encoded by these two genes have distinct immunological properties. The conservation of nucleotide sequence and exon/intron boundaries in these two genes suggests they arose by duplication of a common primordial gene around 150 million years ago.     

Predicted expansion of the claudin multigene family

Claudins (Cldn) are essential membrane proteins of tight junctions (TJs), which form the paracellular permselective barrier. They are produced by a multi-gene family of 24 reported members in mouse and human. Based on a comprehensive search combined with phylogenetic analyses, we identified three novel claudins (claudin-25, -26, and -27). Quantitative RT-PCR revealed that the three novel claudins were expressed in a tissue- and/or developmental stage-dependent manner. Claudins-25 and -26, but not claudin-27, were immunofluorescently localized to TJs when exogenously expressed in cultured MDCK and Eph epithelial cell lines. These findings expand the claudin family to include at least 27 members.     

Structures of two HaeIII-type genes in the human salivary proline-rich protein multigene family

Two members of the human salivary proline-rich protein (PRP) multigene family have been isolated and completely sequenced. These PRP genes, PRH1 and PRH2, are of the HaeIII-type subfamily and code for acidic PRP proteins. Both genes are approximately 3.5 kilobase pairs (kb) in length and contain four exons. Exon 3 encodes the proline-rich part of the protein and includes five 63-base pair (bp) repeats. CAT and ATA boxes and several possible enhancer sequences occur in a 1-kb region 5' to exon 1. Two sets of repeats occur in the sequenced region in addition to the 63-bp repeats: one pair of about 140 bp flanks 500 bp of DNA in the first intervening sequence, and the other pair of 72 bp is tandemly repeated 1.4 kb 5' to the PRH1 gene. The 4-kb region of sequenced DNA from PRH1 differs by an average of 8.7% from the same region in PRH2, but the nucleotide sequences of the exon 3 of the two genes differ by only 0.2%. This result suggests the occurrence of a recent gene conversion event. The regions containing the 5-fold repeated sequences of 63 bp are identical in the two genes, PRH1 and PRH2. A comparison of the human HaeIII and BstNI subfamily repeats and a comparison of the human, mouse, and rat repeats suggest that the individual repeats have evolved in a concerted fashion within each gene and within the PRP gene family as a whole.     

A multigene family encodes the human cytomegalovirus glycoprotein complex gcII (gp47-52 complex).

The HXLF (HindIII-X left reading frame) gene family is a group of five genes that share one or two regions of homology and are arranged in tandem within the short unique component of the human cytomegalovirus genome (K. Weston and B.G. Barrell, J. Mol. Biol. 192:177-208, 1986). These genes were cloned into an SP6 expression vector in both the sense and antisense orientations. An abundant 1.62-kilobase (kb) bicistronic mRNA, predicted to originate from HXLF1 and HXLF2, was detected in the cytoplasm of infected human fibroblast cells by Northern (RNA) blot analysis. Less abundant RNAs of 1.0 and 0.8 kb, predicted to originate from the HXLF5 and HXLF2 genes, respectively, were also detected. Monocistronic, bicistronic, and polycistronic RNAs synthesized in vitro by using SP6 polymerase were translated in rabbit reticulocyte lysates with or without canine pancreatic microsomal membranes. The HXLF1 or the HXLF1 and HXLF2 translation products were detected when the above mRNAs were used. The HXLF3, HXLF4, and HXLF5 gene products were not detected by in vitro translation of the SP6-derived polycistronic mRNA. Nonglycosylated or glycosylated HXLF1 and HXLF2 gene products were immunoprecipitated by monoclonal antibody 9E10, which is specific for a virion envelope glycoprotein complex designated gcII (gp47-52 complex). In addition, the monoclonal antibody 9E10 immunoprecipitated a diffuse glycoprotein band, designated gp47-52, from HCMV-infected cell lysates. The amino acid composition of gp47-52 purified from viron envelopes has the highest similarity to the predicted amino acid composition of the HXLF1 plus HXLF2 open reading frames, but it is more similar to HXLF2 than to HXLF1. The Northern blot results imply that gp47-52 is synthesized predominantly from the abundant 1.62-kb bicistronic mRNA encoded by the HXLF1 and HXLF2 genes. However, the glycoprotein could also be synthesized by the monocistronic 0.8-kb mRNA encoded by the HXLF2 gene as well as by the mRNAs predicted from the other HXLF genes.