|
|
||||||||




* Childrens Hospital of Johannes Gutenberg-University of Mainz, Mainz, Germany;
Institute of Immunology, University of Heidelberg, Heidelberg, Germany;
Department of Biology, Massachusetts Institute of Technology, Cambridge, MA 02139; and
Institute of Medical Microbiology and Hygiene and Microbiology, Johannes Gutenberg-University of Mainz, Mainz, Germany
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
and C5
. After activation by the C5 convertases, C5 is cleaved into C5a and C5b. C5a is the most potent complement-derived proinflammatory peptide (3, 4). It serves not only as a chemoattractant but also initiates multiple defense mechanisms in leukocytes. Through its multiple binding sites C5b initiates and directs the assembly of the membrane attack complex.
C5 is encoded by a single gene copy on chromosome 9q34.1 which contains 41 exons that span a genomic region of 79 kb. The gene encodes the open reading frames for C5
(exon 116), as well as for C5
(exons 1641). Both are transcribed as a single 6.0 kb mRNA and translated into a pre-C5 protein (in
-
orientation). These single chain precursors are processed into the mature form of C5 by the removal of the four internal residues to yield a two-chain protein (5).
Through activation and interaction with respective receptors on various immune cells, complement is closely linked to the adaptive immune response. A normally functioning complement system is also required for physiological tissue regeneration and repair. Complement defects increase the susceptibility to infection and are frequently associated with autoimmune disorders (6, 7). Patients with dysfunction of C5 display a propensity for severe recurrent infections, particularly with Neisseriae, the causative pathogen for meningitis and extragenital gonorrhea. Furthermore, C5 deficiency produces an anti-inflammatory phenotype that may also be involved in certain forms of asthma (8).
Hereditary C5 deficiency has been reported in several families of different ethnic backgrounds and from different geographic regions (7). To date, only two heterozygous nonsense mutations of the human C5 (hC5)3 gene have been reported to account for complete C5 deficiency. Whereas in two C5-deficient families of African-American origin a nonsense mutation located in exon 1 (Gln1 to STOP) was found, another mutation in exon 36 (Arg1458 to STOP) was identified in a third family. Because the affected individuals reported so far were heterozygous for the mutations and thus are believed to be able to express one normal C5 allele, it was speculated that a complete C5 deficiency requires additional molecular defects (9). The majority of recognized mutations which cause diseases occur in the coding sequences of genes with a direct impact on the encoded protein structures and their function. However, another probably large and yet unrecognized proportion of diseases are caused by genetic aberrations resulting from mutations that influence gene transcription, including mRNA processing, export, stability, as well as translational control (10). Based on the fact that alternatively spliced variants have a tissue-specific distribution, alternative splicing has been also considered as a reason for compound heterozygous C5 deficiency with only one open-reading frame being mutated. So far, none of the analyzed splice sites responsible for pre-mRNA processing harbored a mutation that could lead to incorrect splicing of the hC5 gene. However, recent findings indicate that exonic sequences not only encode for protein information, but also contain sequences essential for correct mRNA splicing. These specific short oligonucleotide sequences, termed exonic splicing enhancers (ESEs), increase pre-mRNA splicing and, therefore, play an important role in constitutive and alternative splicing. Such ESEs often activate specific splice sites to facilitate the recognition of an exon by the splicing apparatus and thus, result in exon skipping. Therefore, gene sequence analysis has to be extended to enable the determination of the impact of a sequence alteration on putative ESE motifs.
In this study, we report a new case with selective near total C5 deficiency as revealed by immunochemical and functional analyses. Sequencing of the hC5 coding region did not explain the lack of C5 protein because the mutation, a conservative amino acid exchange from an arginine to a lysine codon, was not expected to produce a functional defect. However, the identified sequence alteration that affects splicing signals (either by ESE alteration or possibly by changing the 5' splicing site itself), leading to exon skipping and a truncated transcript, is likely responsible for the observed phenotype.
| Materials and Methods |
|---|
|
|
|---|
The propositus, a 5-year-old boy, is the third child of consanguineous Turkish parents. In the first year of life, he underwent cardiac surgery because of coarctation of the aorta and a ventricular septal defect. His development proceeded normally until the age of 3 years when he was admitted to the hospital because of headache, fever, vomiting, and abdominal purpura. Meningitis was suspected but cultures of the cerebrospinal fluid were negative. Neisseria meningitides was isolated in blood cultures and after treatment with cefotaxime for 10 days he recovered completely. Eleven months later, the patient was admitted again because of headache, fever, vomiting, hypotension, and dehydration. Cultures of blood, cerebrospinal fluid, and urine were negative. Because of a pronounced leukocytosis and elevated C-reactive protein levels in serum (196 mg/L) the patient was treated with cefuroxime and could be discharged from the hospital after 11 days. Six weeks later, the patient presented with headache and nuchal rigidity but no other signs of meningitis. He was treated again with cefuroxime because of elevated serum levels of C-reactive protein (145 mg/L), IL-2R (1350 U/ml), and an erythrocyte sedimentation rate, even though cultures of blood and urine gave negative results. For comprehensive complement analysis, serum and plasma samples were obtained from venous blood collected from all members of the family. Samples were frozen at 70°C until assayed. After complement C5 deficiency was diagnosed, the patient was placed on long-term prophylaxis with oral penicillin. Since then, he had no further recurrence of severe infections. His parents and a brother had decreased serum levels of complement C5 (Table I), but no history of severe infections. All family members gave written informed consent to undergo molecular studies.
|
Hemolytic assays. Functional activity of the classical (CH50, Ref. 11) and the alternative pathway (AH50, Ref. 12) of complement was measured according to described procedures. Hemolytic activity of C5 was tested using C5-deficient human serum (Calbiochem Merck Biosciences). C6-C9 function was titrated with Ab-sensitized sheep erythrocytes (EA), preincubated with C6-deficient rabbit serum to collect EAC15. Plasma concentrations of C4 and C3 were determined by nephelometry (Behring). C3dg/C3d was measured by double-decker rocket immunoelectrophoresis using rabbit anti-C3c in the lower and rabbit anti-C3d Abs (DAKO) in the upper gel (13). C5 protein concentration was determined by ELISA: in brief, microtiter plates (Nunc) were coated with rabbit anti-human C5 IgG (DAKO). After blocking nonspecific binding sites with 1% (BSA in PBS), appropriate dilutions of the samples were added. Following incubation, C5 was detected by a monoclonal mouse anti-human C5 IgG (Quidel), followed by the appropriate peroxidase-labeled third Ab (Dianova). The reaction was visualized by the addition of o-phenylenediamine/H2O2 substrate. Purified C5 (Calbiochem) was used as standard. SC5b-9 ELISA was performed as described before (14). In some experiments, the patients and control serum samples were preactivated with zymosan to generate SC5b-9 as a measure of terminal activation pathway.
C5 reconstitution. C5 reconstitution assays were performed to investigate restoration of total hemolytic activity (CH50) and chemotactic activity by addition of incremental amounts of purified C5 (Sigma-Aldrich) to the patients serum.
SDS-PAGE and Western blot analysis. To further investigate the presence and protein structure of the patients C5, serum (native and immunoprecipitated with rabbit anti-C5 IgG; DAKO) was subjected to SDS-PAGE (15) under nondenaturing and denaturing conditions. Proteins were transferred onto a nylon membrane (16) and probed by chemoluminescence-enhanced Western blot and immunodetection of C5-specific bands with various monoclonal and polyclonal C5-specific Abs.
Chemotaxis assay. Chemotaxis of polymorphonuclear neutrophil (PMN) was evaluated by using a two-compartment, blind-well Boyden chamber (17). A total of 20 ml of heparinized venous blood from healthy donors was layered on a sodium metrizoate/dextran gradient (Polymorphoprep; Immuno) in HBSS. After centrifugation at 470 x g for 30 min, the PMN-rich cell layer was washed with HBSS (without Ca2+/Mg2+) and residual erythrocytes were removed by hypotonic lysis. Subsequently, the cell pellets (>98% purity) were resuspended in HBSS, containing 0.25% BSA (Serva), and adjusted to 1 x 106 cells/ml.
Viability of PMN was found to be >98% as judged from trypan blue staining. A total of 0.9 ml of the cell suspensions (1 x 106/ml) was added to the upper compartment of the chemotaxis chamber separated by a nitrocellulose membrane (3.0 µm pore size, 13 mm in diameter, 140 µm thickness; Sartorius) from the lower compartment, containing 0.2 ml of the chemoattractant or buffer as control. All samples and cell suspensions were in HBSS with 2 mM Ca2+/0.8 mM Mg2+ (HBSS2+) and 0.25% BSA, pH 7.35. To analyze restoration of chemotactic activity, one aliquot of the patients serum was resubstituted with C5 to receive a final concentration corresponding to that in normal human serum (75 µg/ml). After incubation for 75 min at 37°C, filters were removed, fixed in 100% propanol, stained with hematoxylin, and mounted on microscope slides.
Migration into the filter was quantified by the leading front method (18). Duplicate chambers were used in each experiment and four microscopic fields (x400) for each filter were semiautomatically evaluated by an image analyzer (VRZ-U2; Leitz). Buffer control (random cell movement) was subtracted from each sample, and data are presented as migration distance (µm/75 min).
Genetic analyses
Preparation of RNA and DNA amplification and sequencing the hC5 gene. Genomic DNA from all available family members was isolated from peripheral blood using the QIAamp DNA Blood Mini kit (Qiagen). Coding sequences of the hC5 gene were then amplified with primers flanking all 41 exons and the promoter region 5' of exon 1 (PCR conditions and primer sequences will be supplied upon request). The PCR products were purified by isopropanol precipitation and sequenced using an automated sequencing system (A377; Applied Biosystems). RNA was prepared from lymphocytes and isolated using the QIAamp RNA Blood Mini kit (Qiagen), avian myeloblastosis virus reverse transcriptase (Promega), and a mixture of oligo(dT) primers and random hexamers were used to synthesize the first strand of cDNA from lymphocyte-derived total RNA, according to the manufacturers protocols. Oligonucleotide primers for amplification and sequencing the part of the hC5 cDNA encompassing exons 912 were C5-9f (5'-CAGCAGTCAAAGAACTGTCA-3') and C5-12r (5'-TCAGATGTTCTCCCACTAGC-3'). With these primers, amplification of cDNA containing the wild-type allele is expected to result in a 540-bp fragment.
Sequence analysis.
Obtained sequences containing the exonic sequences, intron/exon boundaries, and putative promoter sites were compared with the published wild-type sequence (Accession no. NM_001735) using the SEQUENCHER (Gene Codes) computer program that allows detection of mutations/polymorphisms. Variations were then annotated with predicted ESEs using the RESCUE-ESE web server (
http://genes.mit.edu/burgelab/rescue-ese/
). The RESCUE method relies upon the assumption that sequences that enhance the recognition of a splice site will occur more frequently in the vicinity of suboptimal splice sites than consensus splice sites. Although the details of the RESCUE protocol have been described previously (19), the basic approach was to count hexamer frequencies in databases of exons and introns that have been sorted according to splice site strength. The splice site scores were used to separate the exon database into "weak" and "strong" sets. The weak exon dataset was composed of the exons that were adjacent to splice sites that score in the bottom quartile of all the splice sites tested while the strong exons had splice sites that scored in the top quartile. As each exon is flanked by a 3' splice site and a 5' splice site, there were a total of four new sets created from the original exon database. The t statistic was used to identify cases of hexamers that had significantly different counts in: 1) exons relative to introns and 2) exons with weak (low scoring) splice sites relative to exons with strong (high scoring) splice sites. A list of hexamers and t scores was received which enabled the identification of the subset of hexamers that were significantly enriched in sequences that were both exonic and adjacent to a weak splice site (19).
Wild-type and mutant sequences of C5 were scored for agreement to the splice site consensus using a maximum entropy model for human splice sites (20).
| Results |
|---|
|
|
|---|
Functional activity of CH50 and AH50 in the patients serum hemolytic activity was almost undetectable (Table I), suggesting a deficiency of one of the terminal complement components, C5-C9. Functional and immunochemical assays revealed a deficiency of the fifth component of the complement system (Table I), while C3, C4, and C6-C9 were normal. Normal plasma levels of the C3 split product C3dg/C3d excluded a secondary deficiency due to complement consumption. In vitro C5 reconstitution fully restored total hemolytic activity (CH50, Fig. 1).
|
C5 ELISA revealed only trace amounts of C5 in the patients serum (Table I), providing the basis for close to zero C5 and CH50 activities. Therefore, we were not able to further characterize C5 protein structure by SDS-PAGE despite application of sensitive chemoluminescence-enhanced Western blot analysis. Zymosan activation of the patients serum led to a generation of only 22 ng/ml SC5b-9 in comparison to
120 µg/ml obtained in normal human control serum, further supporting the presence of only trace amounts of functionally active C5 in the patients serum.
As expected, only negligible chemotactic activity could be generated in the patients serum upon yeast activation. The addition of purified C5 restored (C5a-mediated) chemotactic activity to a level observed in normal serum (Fig. 2).
|
Direct sequencing of the hC5 gene of the propositus which included all exons and flanking regions including the promoter region revealed only a substitution of the normal adenine 1115 with a guanine, according to the sequence of Haviland et al. (5) (NM_001735). This mutation predicts the replacement of a normal lysine (K) with an arginine (R) in codon 372 (K372R), both belonging to the group of alkaline amino acids. Whereas the propositus was homozygous for the mutation, both parents and an older brother were heterozygous carriers of K372R. A sister was homozygous for the wild type. The mutation could not be detected in fifty normal individuals from the same ethnic background.
C5 mutation affects a putative ESE
We used the RESCUE-ESE program (19) to evaluate whether the sequence alteration affects an ESE and to determine its possible impact on the quality of the ESE. Indeed, the nucleotide exchange was likely to affect but not to abolish an ESE (ATCAAG into CATCAG, underlining indicates the position of the exchanged nucleotide in respect of its position in the hexamer sequence of the putative ESE, Fig. 3 and Table II). The set of 238 hexamers that were identified as RESCUE-ESEs were isolated on the basis of their enrichment in exons that are adjacent to weak sites. The t statistic was used to identify hexamers that occurred significantly more frequently in exons (relative to introns) and, within this set of exons, the t statistic was used again to identify hexamers that occurred more frequently in exons that were defined by suboptimal splice sites than exons adjacent to consensus sites. Both these comparisons were arranged such that positive t scores reflect enrichments that are consistent with ESEs (i.e., ESEs are enriched in exons relative to introns; ESEs are enriched in weak exons relative to strong exons). To better understand this ESE alteration we examined the t scores of the set of six overlapping hexamers that were altered by this mutation. In other words, we compared the degree to which each hexamer in the wild-type and mutant allele of the sequence was enriched in the attributes that define ESEs (Fig. 3A). As shown in Fig. 3B, hexamers derived from wild-type sequence (blue squares) are distributed more toward the upper right (enriched in both the weak splice site and the exon set) than those of the mutant sequence (red squares). This difference is consistent with the idea that the wild-type sequence is more likely to possess ESE activity than the mutant sequence.
|
|
We then isolated total RNA from lymphocytes of all family members and amplified the region from exons 9 to 12 as cDNA. Indeed, in all individuals having at least one normal allele we detected the expected normal 540 bp band. However, the propositus did not have the expected band. Instead, he was homozygous for a 424 bp fragment, which was also present in all heterozygous family members harboring the K372R amino acid exchange (Fig. 4A). Direct sequencing from cDNA confirmed skipping of exon 10 (116 bp) in all individuals with the mutated allele. This skipping leads to a frame shift and a premature STOP codon of the remaining open reading frame, abolishing the template for the
- and the
-chain portions of the open reading frame (Fig. 4B). Thus, the sequence alteration could represent a noncanonical splicing mutation that alters the efficiency of an ESE, resulting in the loss of exon 10. Another explanation for this loss of splicing mutation is that the nucleotide substitution is close to the 5' exon/intron junction and may consequently disrupt the recognition of the 5' splice site. Applying a maximum entropy model of splice sites, the mutation reduces the agreement of the C5 5' splice site with the consensus (MAXENT score 7.24 is reduced to 2.55 with the identified mutation). However, this change does not appear to invalidate agggtagat as a 5' splice site, as this exact sequence is used as the 5' splice site for exon 4 of the human topoisomerase (DNA) II binding protein 1 (TOPBP1) (NM 007027).
|
| Discussion |
|---|
|
|
|---|
There is certain evidence that ESEs are comprised of diverse sequences and occur frequently within exons. Ominously, these latter studies predict that many human genetic diseases linked to mutations within exons might be caused by the inactivation of ESEs (21). If the ESE alteration has the proposed effect, it should result in abnormal splicing of the affected exon 10 in the mutated allele. We therefore sequenced cDNA directly which was transcribed from lymphocyte RNA and found skipping of exon 10 in all family members who carry the mutated allele. Because this mutation leads to a frame shift and a premature stop codon after 336 aa, the reduced level of C5 in the heterozygous carriers and the almost complete C5 deficiency in the patient can be explained. To rule out the possibility of intronic mutations, which may also affect splicing, we sequenced the branching site in intron 9. No sequence alteration was detected. This again argues for the ESE alteration being the mechanism for the C5 protein deficiency. As mutations that cause exon skipping are, in general, harmful, they are eliminated by natural selection. To prove that the A to G transition is such an ESE mutation we sequenced exon 10 of 50 healthy blood donors of the same ethnic background and did not detect the sequence alteration in any individual (data not shown). These findings strongly suggest that the C5 deficiency reported here is due to A1115G and that our observation does not reflect a common polymorphism in the Turkish population.
The position of the sequence alteration could also indicate that the mutation, rather than affecting the efficiency of an ESE, disrupts the recognition of the 5' splice site. Indeed, the mutation does reduce the agreement of the C5 exon 10 5' splice site according to the maximum entropy model of splice sites. However, the altered 5' splice site seems to be still efficient, since the identical sequence motive (agggtagat) is the functional 5' splice site for exon 4 of TOPBP1. This again argues for the nucleotide substitution being deleterious due to subsequent ESE alteration.
A growing number of disease-causing mutations have been shown to inactivate ESEs, thereby causing exon skipping. Well known ESE mutations affect the genes BRCA2, GH1, and CFTR. Many of the deleterious BRCA2 mutations have an effect on the splicing of the gene transcript. As in our case, they are typically located close to intron/exon boundaries. Fackenthal and coworkers (22) systematically examined several BRCA2 mutations for potential ESE disruption mutations and identified a Thr2722-to-Arg mutation. It segregated with the affected individuals of a family with breast cancer and disrupted three potential ESE sites. This mutation caused deleterious protein truncation due to exon skipping and suggested its analysis as a potentially useful method to elucidate the clinical significance of a subset of the hitherto unclassified variants of BRCA1 and BRCA2. Furthermore, in growth hormone type 1 (GH1) deficiency, three different mutations have been described (23, 24, 25) which are not located at the 5' splice site of an intron but alter splicing of GH1. All three mutations reside within sequences that resemble exonic and intronic splicing enhancers. Ryther and coworkers (26) showed that these mutations cause splicing failures due to a defect in exon definition and that splicing enhancer mutations which weaken the recognition of exon 3 produce variable amounts of the dominant-negative 17.5 kDa GH1 isoform. In a study presented by Aznarez et al. (27), the consequence of two cystic fibrosis (CF) disease-causing mutations, E656X and 2108delA, was analyzed with respect to the function of a putative ESE in exon 13 of the CF transmembrane conductance regulator (CFTR) gene. Using minigene constructs, it was demonstrated that both mutations caused aberrant splicing in a predicted manner supporting a role for the putative ESE sequence in pre-mRNA splicing. They concluded that the severity of CF could be modulated by changes in the fidelity of CFTR pre-mRNA splicing. Thus, in all three genes as well as in the current case, a particular single base exchange in an exon causes inappropriate skipping by inactivation of a critical ESE motif.
In this report we present the first case of a homozygous mutation in the hC5 gene, leading to an almost complete complement C5 deficiency. We demonstrate that an initially inconsiderable sequence alteration results in exon skipping by affecting an ESE, a mechanism that will shed a new light on the role of putative silent/conservative mutations in disease-associated genes. Our findings indicate that putative polymorphic sequence alterations which segregate with defects, should be reviewed in this respect. Whether such noncanonical ESE mutations may explain variable phenotypes in other diseases, as it was suggested by Aznarez and coworkers (27), needs to be investigated.
| Disclosures |
|---|
|
|
|---|
| Acknowledgments |
|---|
| Footnotes |
|---|
1 N.P., D.P., and M.K. contributed equally. ![]()
2 Address correspondence and reprint requests to Dr. Joachim Pohlenz, Childrens Hospital of Johannes Gutenberg-University of Mainz, Langenbeckstrasse 1, Building 109, D-55101 Mainz, Germany. E-mail address: pohlenz{at}kinder.klinik.uni-mainz.de ![]()
3 Abbreviations used in this paper: hC5, human C5; ESE, exonic splicing enhancer; CH50, hemolytic activity of the classical pathway; AH50, hemolytic activity of the alternative pathway; PMN, polymorphonuclear neutrophil; SNP, single nucleotide polymorphism; GH1, growth hormone type 1; CF, cystic fibrosis. ![]()
Received for publication September 23, 2004. Accepted for publication January 25, 2005.
| References |
|---|
|
|
|---|
-chain deficiency in two siblings. Immunopharmacology 38:215.[Medline]
Related articles in The JI:
This article has been cited by other articles:
![]() |
Z. Wang and C. B. Burge Splicing regulation: From a parts list of regulatory elements to an integrated splicing code RNA, May 1, 2008; 14(5): 802 - 813. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Zou, R. K. Gopalraj, J. Lok, H. Zhu, I-F. Ling, J. F. Simpson, H. M. Tucker, J. F. Kelly, S. G. Younkin, D. W. Dickson, et al. Sex-dependent association of a common low-density lipoprotein receptor polymorphism with RNA splicing efficiency in the brain and Alzheimer's disease Hum. Mol. Genet., April 1, 2008; 17(7): 929 - 935. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Zhu, H. M. Tucker, K. E. Grear, J. F. Simpson, A. K. Manning, L. A. Cupples, and S. Estus A common polymorphism decreases low-density lipoprotein receptor exon 12 splicing efficiency and associates with increased cholesterol Hum. Mol. Genet., July 15, 2007; 16(14): 1765 - 1772. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. P. Gorlov, M. Kimmel, and C. I. Amos Strength of the purifying selection against different categories of the point mutations in the coding regions of the human genome Hum. Mol. Genet., April 1, 2006; 15(7): 1143 - 1150. [Abstract] [Full Text] [PDF] |
||||
![]() |
D.-Q. Xu and W. Mattox Identification of a splicing enhancer in MLH1 using COMPARE, a new assay for determination of relative RNA splicing efficiencies Hum. Mol. Genet., January 15, 2006; 15(2): 329 - 336. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |