Evolutionary relationships among proteins encoded by the regulator of complement activation gene cluster

Evolutionary relationships among members of the regulator of complement activation (RCA) gene cluster were analyzed using neighbor-joining and parsimony methods of phylogenetic tree inference. We investigated the structural and functional similarities among short consensus repeats (SCRs) of the following human proteins: the alpha chain of the C4b-binding protein (C4bpalpha), factor H (FH), factor H-related proteins (FHR-1 through FHR-4), complement receptors type 1 (CR1) and type 2 (CR2), the CR1-like protein (CR1L), membrane cofactor protein (MCP), decay accelerating factor (DAF), and the sand bass proteins, the cofactor protein (SBP1) and its homolog, the cofactor-related protein (SBCRP-1). Also included are the beta chain of the human C4b-binding protein (C4bpbeta) and the b subunit of human blood-clotting factor XIII (FXIIIb). Our results indicate that the human plasma complement regulators, FH and C4bpalpha, fall into two distinct groups on the basis of their sequence divergence. Homology among RCA proteins is in agreement with their chromosomal location, with the exception of C4bpbeta. The evolutionary relationships among individual short consensus repeats are confirmed by the exon/intron structure of the RCA members. Structural similarities among repeats of the RCA proteins correlate with their functional activities and demonstrate the importance of the N-terminal SCRs.     

Structure of the gene for human complement protein decay accelerating factor.

Decay accelerating factor (DAF) is a glycophospholipid-anchored membrane protein that is part of the regulators of complement activation (RCA) gene family located on human chromosome 1, band q32. These proteins, beginning at their amino terminus, consist largely of multiple copies of an approximately 60 amino acid short consensus repeat (SCR). A DAF cDNA clone was used to identify overlapping bacteriophage genomic clones. The human DAF gene spans approximately 40 kb and consists of 11 exons. The length of these exons and introns varies considerably, with the exons ranging from 21 to 956 bp and the introns ranging from approximately 0.5 to 19.8 kb. SCR I, II, and IV are all encoded by single exons; however, SCR III is encoded by two separate exons, with the splice junction occurring after the second nucleotide of the codon for the glycine residue at position 34 of the consensus sequence. This feature has also been found in CR1, CR2, membrane cofactor protein, and murine factor H. Following the SCR in DAF is a 76 amino acid serine/threonine-rich domain encoded on three separate exons. Exon 10 encodes the Alu family sequence that has been found as an insert in a minor class of DAF cDNA, thus indicating that this mRNA arises by standard alternative splicing. The last DAF exon, which comes after the largest intron of 19.8 kb, encodes the hydrophobic carboxy terminus and the 3'UT region. The nature of the signal that directs posttranslational attachment of a glycophospholipid anchor to DAF is not known, but that signal is apparently spread over three exons and greater than 20 kb. An analysis of the DAF gene provides additional evidence for the common evolutionary heritage of the RCA gene family. The exon/intron structure of this gene will facilitate experiments aimed at understanding the functions of the various domains of DAF.     

Organization, sequence and nuclease hypersensitivity of repetitive elements flanking the chicken apoVLDLII gene: extended sequence similarity to elements flanking the chicken vitellogenin gene.

Analysis of nuclease hypersensitivity of regions flanking the estrogen-dependent, chicken apoVLDLII gene has revealed an hepatic, DNaseI hypersensitive site whose sensitivity is influenced by both the developmental stage and sex of the bird. The site is located 3.0kb upstream from the gene, in a block of middle repetitive elements. Contact hybridization studies indicate that the block consists of contiguous copies of two elements with reiteration frequencies of 500-1000 and 10,000-30,000 copies per haploid genome. Sequencing of 1.8kb spanning the repeats has revealed that the higher frequency element is a member of the CR1 family. The adjacent lower frequency repeat can also be found next to another member of the CR1 family located in the 3' flanking region of the vitellogenin gene. The hypersensitive site has been mapped to one of the two most highly conserved regions of the CR1 element. This region displays homology with a silencer sequence recently identified in a CR1 element flanking the chicken lysozyme gene.     

A new repetitive element of the CR1 family downstream of the chicken vitellogenin gene.

We have analyzed a repetitive DNA sequence found in the 3'-flanking region of the chicken vitellogenin gene. By its sequence, the repetitive DNA has been identified as a hitherto unreported member of the chicken CR1 family of repetitive elements. The CR1 sequence displays the structural characteristics of a long terminal repeat located at the 3' end of an avian retrovirus. The CR1 element lies 2.2 kb downstream of the vitellogenin gene and 'points' away from the gene rather than toward it. In this respect, this element differs from other CR1 repeats. The CR1 element is embedded in a region showing changes in chromatin structure implying a potential role for this sequence in determining the structural state of the local chromatin.     

Characterization of CR1- and membrane cofactor protein-like proteins of two primates

We have identified and characterized C3b binding proteins of two primates, orangutan (Pongo pygmaeus) and gorilla (Gorilla gorilla). Detergent solubilized 125I surface-labeled E and PBMC were subjected to affinity chromatography with homologous or human iC3/C3b. These ligands bound a 225,000 single chain protein from orangutan E and PBMC and a 220,000 protein from gorilla E. Proteins of the same Mr were immunoprecipitated by a rabbit polyclonal and two murine mAb to the human CR1 (CD35). The C3b binding protein of gorilla E aligned with that of the common human CR1 polymorphic size variant. Human or orangutan iC3 was also a ligand for a surface-labeled protein doublet of 59,000 and 65,000 from orangutan E. The doublet pattern and mol wts are similar to membrane cofactor protein (or CD46). Further, this doublet was immunoprecipitated by a mAb to human MCP. The MCP-like protein doublet was not isolated from gorilla or human E. Decay accelerating factor (DAF) of orangutan E was also identified and was structurally and antigenically distinct from the MCP-like protein. Orangutan or gorilla E preparations were a cofactor for the cleavage of human iC3 by human factor I and produced the same cleavage fragments as human CR1. Cofactor activity of orangutan E was partially inhibited by preclearance of CR1 and more completely inhibited by preclearance of MCP. Cofactor activity of gorilla E was inhibited by coincubation with a monoclonal antibody to human CR1. These data indicate that the orangutan and gorilla high m.w. proteins are equivalent to human CR1. The orangutan E membrane protein doublet with m.w. of 59,000 and 65,000 possesses biochemical, antigenic, and functional properties of human membrane cofactor protein.     

The human immunoglobulin κ locus on yeast artificial chromosomes (YACs)

The human immunoglobulin κ locus is a duplicated structure. Contigs of 600 kb with 40 Vκ genes and 440 kb with 36 Vκ genes had been established for the Cκ proximal (p) and distal (d) copies, respectively. In addition the human genome contains more than 24 dispersed Vκ genes, called orphons. In the present study, 22 κ-locus derived YACs were analyzed in detail, while 30 orphon-derived YACs were characterized only with respect to some parameters. The κ-locus derived YACs allowed three gaps to be closed which previously could not be bridged by cosmid and phage λ cloning. At the 5′ side, the p contig was extended in the YACs by 50 kb and the d contig by 16 kb. At the 3′side, the d contig was extended by 11.5 kb. Beyond the 3′ end of the d contig a new Vκ gene was found, which is located, according to pulsed field gel electrophoresis (PFGE) experiments, at a distance of at least 140 kb from the last Vκ gene of the contig. This Vκ gene, which was termed Z0, occurred on three YACs, albeit at distances smaller than 140 kb; this was probably due to deletions in the YACs caused by abundant repetitive sequences at the borders of the locus. According to its sequence and to the restriction map of its surroundings, Z0 is an orphon gene of the so-called Z family, of which several members are known to be dispersed throughout the genome. The possibility that Z0 has been the parent of the other Z orphons is discussed.     

Characterization of the human laminin beta2 chain locus (LAMB2): linkage to a gene containing a nonprocessed, transcribed LAMB2-like pseudogene (LAMB2L) and to the gene encoding glutaminyl tRNA synthetase (QARS).

The laminin beta2 chain is an important constituent of certain kidney and muscle basement membranes. We have generated a detailed physical map of a 110-kb genomic DNA segment surrounding the human laminin beta2 chain gene (LAMB2) on chromosome 3p21.3-->p21.2, a region paralogous with the chromosome 7q22-->q31 region that contains the laminin beta1 chain gene locus (LAMB1). Several CpG islands and a novel polymorphic microsatellite marker (D3S4594) were identified. The 3' end of LAMB2 lies 16 kb from the 5' end of the glutaminyl tRNA synthetase gene (QARS). About 20 kb upstream of LAMB2 we found a gene encoding a transcribed, non-processed LAMB2-like pseudogene (LAMB2L). The sequence of 1.75 kb of genomic DNA at the 3' end of LAMB2L was similar to exons 8-12 of the laminin beta2 chain gene. The LAMB2L-LAMB2-QARS cluster lies telomeric to the gene encoding the laminin-binding protein dystroglycan (DAG1).     

DNase I sensitive domain of the gene coding for the glycolytic enzyme glyceraldehyde-3-phosphate dehydrogenase

We have cloned a 36-kilobase segment of chicken DNA containing the gene coding for glyceraldehyde-3-phosphate dehydrogenase [GAPDH (EC 1.2.1.12)], a glycolytic enzyme which is expressed constitutively in all cell types. Using defined segments of this cloned DNA as probes, we have determined the DNase I sensitive domain of the GAPDH natural gene in the hen oviduct. When nuclei isolated from hen oviduct were treated with DNase I under conditions known to preferentially degrade actively transcribed genes (i.e., 15-20% of the DNA rendered perchloric acid soluble), a region of approximately 12 kilobase(s) (kb) containing the GAPDH coding sequences and flanking DNA was found to be highly susceptible to digestion by DNase I. Approximately 4 kb downstream from the end of the coding sequences, there was an abrupt transition from the DNase I sensitive or "open" configuration to the resistant or "closed" configuration. The chromatin then remained in a closed conformation for at least 10 kb further downstream. On the 5' side of the gene, the transition from a sensitive to a resistant configuration was located about 4 kb upstream from the gene. In addition, we have localized two repeated sequences in the area of DNA that was cloned. One of these is of the CR1 family of middle repetitive elements. It is located about 18 kb 3' to the gene and as such lies well outside of the DNase I sensitive region which encompasses GAPDH. The other repetitive element is of an uncharacterized family. It is located upstream from the gene and appears to be within a region of transition from the DNase I sensitive to resistant states.     

Clustered cadherin genes: a sequence-ready contig for the desmosomal cadherin locus on human chromosome 18

We describe the assembly of a cosmid and PAC contig of approximately 700 kb on human chromosome 18q12 spanning the DSC and DSG genes coding for the desmocollins and desmogleins. These are members of the cadherin superfamily of calcium-dependent cell adhesion proteins present in the desmosome type of cell junction found especially in epithelial cells. They provide the strong cell-cell adhesion generated by this type of cell junction for which expression of both a desmocollin and a desmoglein is required. In the autoimmune skin diseases pemphigus foliaceous and pemphigus vulgaris (PV), where the autoantigens are, respectively, encoded by the DSG1 and DSG3 genes, severe areas of acantholysis (cell separation), potentially life-threatening in the case of PV, are evident. Dominant mutations in the DSG1 gene causing striate palmoplantar keratoderma result in hyperkeratosis of the skin on the parts of the body where pressure and abrasion are greatest, viz., on the palms and soles. These genes are also candidate tumor suppressor genes in squamous cell carcinomas and other epithelial cancers. We have screened two chromosome 18-specific cosmid libraries by hybridization with previously isolated YAC clones and DSC and DSG cDNAs, and a whole genome PAC library, both by hybridization with the YACs and by screening by PCR using cDNA sequences and YAC end sequence. The contigs were extended by further PCR screens using STSs generated by vectorette walking from the ends of the cosmids and PACs, together with sequence from PAC ends. Despite screening of two libraries, the cosmid contig still had four gaps. The PAC contig filled these gaps and in fact covered the whole locus. The positions of 45 STSs covering the whole of this region are presented. The desmocollin and desmoglein genes, which are about 30-35 kb in size, are quite well separated at approximately 20-30 kb apart and are arranged in two clusters, one DSC cluster and one DSG cluster, which are transcribed outward from the interlocus region. The order of the genes is correlated with the spatial order of gene expression in the developing mouse embryo, and this, and previous transgenic experiments, suggests that long-range genetic elements that coordinate expression of these genes may be present. The complete bacterial clone contig described in this paper is thus a resource not only for future sequencing but also for investigations into the control of expression of these clustered genes.     

Junctions between genes in the haptoglobin gene cluster of primates.

To investigate the nature of the recombination that generated the haptoglobin three-gene cluster in Old World primates, we sequenced the region between the second gene (HPR) and the third gene (HPP) in chimpanzees (15 kb), as well as the region 3' to the cluster in humans (14 kb). Comparison to the previously sequenced human haptoglobin (HP) and HPR genes showed that the junction point between HP and HPR in humans (junction 1) was not identical to the junction point between the HPR and HPP genes of the chimpanzee (junction 2). An Alu sequence was found at each junction, but both Alu sequences lacked short direct repeats of the flanking genomic DNA. The lack of direct repeats implies that both junction Alu sequences are the products of recombination between different Alu elements. In addition, other insertion and deletion events are clustered in the regions near the junction Alu sequences. The observation that Alu sequences define the junctions between genes in the haptoglobin gene cluster emphasizes the importance of Alu sequences in the evolution of multigene families.     

Structure and characterization of a murine chromosomal fragment containing the interferon beta gene

In this paper we report on the cloning and characterization of the murine interferon (IFN) beta gene. We have isolated and sequenced a 2.8 kb genomic fragment containing the murine IFN beta gene flanked by 1.2 kb 5' and 1 kb 3' untranslated regions (1 kb = 10(3) base-pairs). The mRNA cap site has been defined. An extensive analysis of the flanking sequence is provided and points out striking features such as: the presence of A + T-rich motifs characteristic of transiently expressed mRNAs, and homologies to repetitive R-type element flanks and to hormone-responsive elements. Comparison of the MuIFN beta 5' flanking region with those from other species reveals similarities in the sequences required for the regulated expression of such inducible genes. Computer analysis of the 130 base-pairs preceding the cap site has revealed TGAAAG motifs and shows that the presence of such elements and their permutants have biological significance, according to statistical calculations. Thus, the comparison between the mouse promoter reported here and the promoters from other species highlights the region containing the hexanucleotide blocks, which is strongly conserved.     

Molecular cloning and characterization of the human beta-like globin gene cluster

The genes encoding human embryonic (epsilon), fetal (G gamma, A gamma) and adult (delta, beta) beta-like globin polypeptides were isolated as a set of overlapping cloned DNA fragments from bacteriophage lambda libraries of high molecular weight (15-20 kb) chromosomal DNA. The 65 kb of DNA represented in these overlapping clones contains the genes for all five beta-like polypeptides, including the embryonic epsilon-globin gene, for which the chromosomal location was previously unknown. All five genes are transcribed from the same DNA strand and are arranged in the order 5'-epsilon-(13.3 kb)-G gamma-(3.5 kb)-A gamma-(13.9 kb)-delta-(5.4 kb)-beta-3'. Thus the genes are positioned on the chromosome in the order of their expression during development. In addition to the five known beta-like globin genes, we have detected two other beta-like globin sequences which do not correspond to known polypeptides. One of these sequences has been mapped to the A gamma-delta intergenic region while the other is located 6-9 kb 5' to the epsilon gene. Cross hybridization experiments between the intergenic sequences of the gene cluster have revealed a nonglobin repeat sequence (*) which is interspersed with the globin genes in the following manner: 5'-**epsilon-*G gamma-A gamma*-**delta-beta*-3'. Fine structure mapping of the region located 5' to the delta-globin gene revealed two repeats with a maximum size of 400 bp, which are separated by approximately 700 bp of DNA not repeated within the cluster. Preliminary experiments indicate that this repeat family is also repeated many times in the human genome.     

An intronic signal for alternative splicing in the human genome

An important level at which the expression of programmed cell death (PCD) genes is regulated is alternative splicing. Our previous work identified an intronic splicing regulatory element in caspase-2 (casp-2) gene. This 100-nucleotide intronic element, In100, consists of an upstream region containing a decoy 3' splice site and a downstream region containing binding sites for splicing repressor PTB. Based on the signal of In100 element in casp-2, we have detected the In100-like sequences as a family of sequence elements associated with alternative splicing in the human genome by using computational and experimental approaches. A survey of human genome reveals the presence of more than four thousand In100-like elements in 2757 genes. These In100-like elements tend to locate more frequent in intronic regions than exonic regions. EST analyses indicate that the presence of In100-like elements correlates with the skipping of their immediate upstream exons, with 526 genes showing exon skipping in such a manner. In addition, In100-like elements are found in several human caspase genes near exons encoding the caspase active domain. RT-PCR experiments show that these caspase genes indeed undergo alternative splicing in a pattern predicted to affect their functional activity. Together, these results suggest that the In100-like elements represent a family of intronic signals for alternative splicing in the human genome.