Analysis of Complement C3 Gene Reveals Susceptibility to Severe Preeclampsia

Preeclampsia (PE) is a common vascular disease of pregnancy with genetic predisposition. Dysregulation of the complement system has been implicated, but molecular mechanisms are incompletely understood. In this study, we determined the potential linkage of severe PE to the most central complement gene, C3. Three cohorts of Finnish patients and controls were recruited for a genetic case-control study. Participants were genotyped using Sequenom genotyping and Sanger sequencing. Initially, we studied 259 Finnish patients with severe PE and 426 controls from the Southern Finland PE and the Finnish population-based PE cohorts. We used a custom-made single nucleotide polymorphism (SNP) genotyping assay consisting of 98 SNPs in 18 genes that encode components of the complement system. Following the primary screening, C3 was selected as the candidate gene and consequently Sanger sequenced. Fourteen SNPs from C3 were also genotyped by a Sequenom panel in 960 patients with severe PE and 705 controls, including already sequenced individuals. Three of the 43 SNPs observed within C3 were associated with severe PE: rs2287845 (p = 0.038, OR = 1.158), rs366510 (p = 0.039, OR = 1.158), and rs2287848 (p = 0.041, OR = 1.155). We also discovered 16 SNP haplotypes with extreme linkage disequilibrium in the middle of the gene with a protective (p = 0.044, OR = 0.628) or a predisposing (p = 0.011, OR = 2.110) effect to severe PE depending on the allele combination. Genetic variants associated with PE are located in key domains of C3 and could thereby influence the function of C3. This is, as far as we are aware, the first candidate gene in the complement system with an association to a clinically relevant PE subphenotype, severe PE. The result highlights a potential role for the complement system in the pathogenesis of PE and may help in defining prognostic and therapeutic subgroups of preeclamptic women.

Introduction

Preeclampsia (PE) is a serious vascular complication of pregnancy, which may lead to a life-threatening multi-organ dysfunction and a convulsive condition, eclampsia (1). PE affects 3–5% of pregnancies in all ethnic groups. The development and progression of the disease are unpredictable with delivery being the only effective cure.

Preeclampsia has been the subject of numerous genetic studies and several associating single nucleotide polymorphisms (SNPs) have been identified. Among the genes where associating SNPs have been described are genes linked to hypertension and vascular and metabolic disease (2–4), all diseases whose risk is increased in the later life of PE patients (5). Furthermore, genes encoding for proteins involved in the immunological processes have also been found to harbor SNPs that predispose patients to PE (6, 7).

Pregnancy is the ultimate immunological paradox, where the maternal immune system must accommodate to protect the mother and growing fetus from pathogens while allowing the semiallograft fetus to persist and thrive. PE is a vascular disease that involves poor placentation (8), especially in the severe and/or early onset (diagnosis or delivery <34 weeks of gestation) forms of the disease, where immunological mechanisms have been implicated (9, 10). Among the immunological effector mechanisms, inadequate control of the maternal complement system has been suggested to contribute to the etiology of PE (11, 12).

The complement system is a part of the immune system that is involved in generating inflammation and mediating the clearance of microbes and injured tissue materials. It can be activated by the classical, the lectin, or the alternative pathway, which proceed stepwise in a controlled cascade of interactions between surface-bound and soluble proteins in the serum. Complement C3 is the central component of all activation pathways. It is among the most ancient components of innate immunity that has evolved over 1,000 million years ago (13). Indeed, the functional domains of the human C3 are conserved in corals and Cnidarians (14, 15). The ancient evolutionary attribute of C3 and its abundance in the human serum indicate its important role as the key component of immunity against infection and in the discrimination between self and non-self (16).

C3 is a large protein formed by pair of disulfide-linked α- and β-chains and 13 individual domains. In shape, it is a “robot”-like molecule that has eight macroglobulin domains “the body,” a linker (LNK), the C3a anaphylatoxin, an arm-like region with the “C1r/s, UEGF, BMP1” (CUB), a thioester-containing domain (TED), an N-terminal domain (α'NT), and the “head” (C345C) linked to the body with an anchor (17). The domains are encoded by 41 exons of the C3 gene. When C3 is activated to C3b, an internal thioester bond is disrupted allowing covalent attachment of C3b to target surfaces. Subsequently, factor B binds to the MG2 and CUB domains of C3b (18). A C3 convertase, C3bBb, is formed when factor D activates C3b-bound factor B to breakdown product. Thereafter, C3bBb cleaves new C3 molecules to C3b to release anaphylatoxic C3a fragments to the circulation (19). The main inhibitors of C3 activation, factor H, decay accelerating factor (CD55), and CR1 bind partially to the factor B-binding site to prevent or disrupt C3bBb formation (17, 20).

Mao et al. showed that alternative pathway complement activation is the key mechanism for reproductive failure in complement inhibitor deficient (Crry^−/−) mice (21). Recently, it was shown that alternative complement pathway becomes activated also in human pregnancies, where severe PE develops (11). Successful treatment of a patient suffering from HELLP syndrome (hemolysis, elevated liver enzymes, low platelet count), a life-threatening complication of PE, by eculizumab, a targeted inhibitor of complement protein C5, demonstrated that the complement system could provide a promising target for drug development in severe PE (22). C5 is the initiator of the final stages of complement activation, i.e., the lytic terminal pathway.

We have looked for SNP association with severe PE among 18 genes coding for the complement system. The most promising associations were found in C3, where linkage both to individual SNPs and to a distinct haplotype, was observed. C3 was thus subsequently chosen for detailed capillary sequencing of its exons and promoter regions (PROMs) in women with severe PE and controls with non-PE pregnancies.

Materials and Methods

This case-control study was conducted at the Department of Medical and Clinical Genetics and Institute for Molecular Medicine Finland in the University of Helsinki. The number of women and methods of genotyping in each stage of the study are described in detail in the work flow chart (Figure 1). All subjects provided a written informed consent in accordance with the Declaration of Helsinki. Study protocols were approved by the local Ethical Committees, specifically, for the FINNPEC study, ethical approval has been obtained from the Coordinating Ethics Committee, Hospital District of Helsinki, and Uusimaa, for the Finnish population-based PE cohort was approved by the ethics committee of the Finnish Red Cross Blood Service and by the Ministry of Social Affairs and Health and for Southern Finland PE study was approved by the Ethics Committee of the Department of Obstetrics and Gynaecology at Helsinki University Central Hospital.

FIGURE 1

Figure 1. The workflow of C3 project with the analysis process and appropriate method in parallel. The single nucleotide polymorphisms (SNPs) that were assessed in each stage are identified by the respective color in the method panel in Figure 2. SPE, severe preeclampsia.

Selection of Study Subjects and Diagnostic Criteria

Study subjects were selected from three Finnish PE cohorts; The Finnish population-based PE cohort (23), The Southern Finland PE cohort (23), and the Finnish Genetics of PE Consortium (FINNPEC) cohort (24) (Table S1 in Supplementary Material). Only women with a singleton pregnancy were included in the study. Except for three women from the Southern Finland PE cohort, where data on subjective symptoms were not available, all women with PE met the criteria of severe PE (23). Severe PE was defined as having repeatedly maximum systolic blood pressure (BP) ≥ 160 mmHg and/or maximum diastolic BP ≥ 110 mmHg or proteinuria ≥5 g/day or significant subjective symptoms in a woman diagnosed with PE according to the American College of Obstetricians and Gynecologists (25).

The clinical characteristics of all women whose samples were used for sequencing and Sequenom genotyping are described in Table 1. Body mass index was defined as the pre-pregnancy weight in kilograms divided by height in meters squared (kg/m²). Pregnancy weight and height measures were obtained from the antenatal charts. Relative birth weight SD units (Z-score) were calculated according to Finnish standards (26).

TABLE 1

Table 1. Clinical characteristics of the Southern Finland cohort (N = 32, data missing for five individuals) and the Finnish Genetics of Preeclampsia (FINNPEC) cohort participants in C3 sequencing and genotyping.

Subjects in Complement Genotyping

The SNP genotyping was performed in 259 PE women and 426 non-PE controls. The PE women and controls were selected from the Finnish population-based PE cohort and the Southern Finland PE cohort with preference on severity of the disease.

Subjects in C3 Sequencing and Microsatellite Analysis

We performed C3 sequencing and the microsatellite analysis of the upstream regulatory motif in 32 severe PE women from the Southern Finland PE cohort (Table 1) and 95 non-PE controls from the FINNPEC cohort. These data were also used in the relative extended haplotype homozygosity (REHH) statistics.

Subjects in C3 Replication by Sequencing

Ninety-five women with severe PE from the FINNPEC cohort were selected for the second stage of sequencing, which involved re-sequencing the middle part of the gene indicated by blue in Figure 2. The replication sequencing data were used in the REHH statistics.

FIGURE 2

Figure 2. The haploblock structure within C3 determined by Sanger sequencing and replication studies. The intensifying gradient of red represents the increasing relative linkage disequilibrium between two variants. The associating Block 2 is identified by the black triangle in the middle of the gene and had a multiallelic D' score of 0.623. The single nucleotide polymorphisms (SNPs) that were assessed in each replication stage are identified by the respective color square that correspond to colors in the material panel in Figure 1. The independently associating SNPs are indicated by stars. All SNPs in this image are listed in Table S4 in Supplementary Material.

Subjects in Replication by Sequenom Genotyping

Fourteen SNPs were validated in 960 women with severe PE and 705 non-PE controls from the FINNPEC cohort (including all of the FINNPEC participants included in the initial exploratory Sequenom genotyping). These data were combined with sequenced genotypes of C3 in the initial phase and replication phase for the single SNP association analysis. The study subjects are described in Table 1.

Complement Genotyping

A custom SNP genotyping was performed using Sequenom’s MassARRAY MALDI-TOF Mass Spectrometry Compact platform and iPLEX Gold chemistry (Sequenom Inc., San Diego, CA, USA) with standard protocols as described elsewhere (27). Briefly, 18 genes coding for components of the complement system were chosen for genotyping and for each gene, SNPs with assumed relevance were selected based on published data on protein function, activity, or disease association (Table S2 in Supplementary Material). Furthermore, we included potentially functional nonsense and missense SNPs, and finally, we also included some intronic, promoter or 3' end SNPs as markers of association. A total of 98 SNPs were assayed in four multiplexes of 15–34 markers each. We focused on SNPs with minor allele frequencies >0.05 in European populations based on the HapMap data (28). Genotypes were analyzed using Sequenom’s MassARRAY Typer version 4.0 software.

Microsatellite Analysis

The size of the (CA)_n repeat polymorphism GF100472 in the promoter area of the C3 gene was determined by fragment analysis. The repeat area was amplified by PCR using FAM-labeled forward primer and non-labeled reverse primer (Table S3 in Supplementary Material) in PCR conditions of an initial denaturation at 95°C for 10 min, followed by 32 amplification cycles of 95°C for 30 s, 67°C for 30 s, and 72°C for 50 s and final extension at 72°C for 7 min. The sizes of the amplified fragments were determined using an automated capillary sequencer ABI3730xl DNA Analyzer (Applied Biosystems). GeneScan 500 LIZ Size Standard (Applied Biosystems) was used to size the fragment data. The number of CA repeats was determined with Gene Mapper v4.0 software (Applied Biosystems) and the repeat alleles were classified as short length (CA_8–10), medium length (CA_11–12), or long (CA₁₅) for the purpose of analysis.

C3 Sequencing

The exonic areas of the C3 gene including flanking intronic regions and potential splice sites and 650 kb PROM were sequenced using standard Sanger sequencing with primers detailed in S3. Amplitaq Gold (Applied Biosystems) enzyme was used in the reactions. Samples were purified from excess primers by digestion using 2.5 U of Shrimp Alkaline Phosphatase USB and 5 U of Exonuclease I (New England Biolabs) at 37°C for 60 min, followed by inactivation of 15 min at 80°C. Purified samples were prepared for sequencing using the BigDye 3.1 terminator (Applied Biosystems) as instructed by the manufacturer. The sequencing reaction was as follows: initiating step of 96°C for 1 min, 25 cycles of 96°C for 10 s, 53°C for 5 s, 60°C for 4 min. Sequence samples were purified with the Millipore Multiscreen plates (Millipore, USA) with Sephadex G-50 Superfine Sepharose (Amersham Biosciences, Sweden). Electrophoresis was performed with an ABI 3730 DNA Analyzer (Applied Biosystems) and base calling using the Sequence Analysis 5.2 software (Applied Biosystems). Initial analysis was carried out in Variant Reporter 1.0 software (Applied Biosystems) and the reported results were checked by Sequencher 4.1.4 software (Gene Codes, USA).

Replication by Sequenom Genotyping

Fourteen SNPs covering the length of the gene (indicated by red rectangles in Figure 2) were included in a single Sequenom iPLEX. The purpose of the Sequenom replication was to increase the sample size to gain reliability of analyses. The assay design and the genotyping were performed with Sequenom MassArray system at the FIMM Technology Centre, University of Helsinki. The Technology Centre performed routine quality control steps to ensure high quality of the genotyping.

Quality Control

Genotyping results from all methods were tested for deviations from the Hardy–Weinberg equilibrium (p < 0.05), and none were discovered in the controls. Data on three individuals were removed due to unresolved discrepancy between the sequenced and the genotyped results and data from 11 individuals were removed due to >10% failed genotyping by the Sequenom (Table S1 in Supplementary Material).

In Silico Analysis of Functionality

Ensemble Variant Effect Predictor and Human Splicing Finder 3.0 online softwares were used to assess the consequences of the five intronic SNPs of interest (29). RNAsnp by Center for non-coding RNA in Technology and Health (RTH) was the final online software that was applied to detect possible local RNA secondary structure changes, which might be introduced by exonic SNPs and could lead to changes in posttranscriptional processes of an otherwise functional gene (30, 31). Mode 2 of the program was used, as it is applicable for large mRNAs (>1,000 nucleotides). Following the suggested limits, a p-value <0.2 was considered indicative of an SNP induced change in secondary structure.

The Relative Extended Haplotype Homozygosity

To predict selection pressure toward the discovered haplotype and pinpoint its structure, we completed a REHH analysis using the associating rs2287845 as a focal SNP (32). REHH was conducted in R following the developer’s instructions to reveal the evolutionary selection pressures underlying the haplotype structure. All Sanger sequenced individuals (initial cohort and replication cohort N = 213) were used for REHH and missing genotypes were imputed using fastPHASE software (33, 34).

Statistical Analyses

The results of the complement SNP genotyping were analyzed in PLINK (35). The C3 sequencing results were analyzed for association in gPLINK and PLINK. Single SNP association to disease was evaluated by Fisher’s Exact test and results were confirmed in Haploview (36). Haplotype analysis was conducted with the Haploview software (36). Association analysis of individual SNPs and the discovered haplotype blocks was done using χ² test in the Haploview program.

Results

SNP Genotypes of Complement Genes

To analyze potential associations between PE and complement genes, we genotyped selected SNPs in 17 genes coding components of the complement system. No differences were observed in 64/72 SNPs between 259 PE women and 426 non-PE controls (data not shown). Out of the remaining eight SNPs, three associated to C3 in genotypic model analysis and in C3, rs2230204, and rs2230205 were associated after permutation. rs2230204 in the proximity of C3 exon 14 was most likely to have an independent effect (likelihood ratio test: χ² = 5.1, df = 1, p = 0.024).

C3 Sequences

In the C3 promoter, exons, and flanking introns, a total of 43 SNPs were identified in severe PE women, non-PE controls, or both (Figure 2). rs200046246 is located in a predicted transcription factor-binding site. It is a missense SNP that causes an amino acid change K779R. It was predicted by SIFT and POLYPHEN2 to be well tolerated and benign, apparently, because there is no change in the charge of the amino acid (both lysine/K and arginine/R are basic residues). Other variants were either synonymous or intronic. Six SNPs were independently associated to severe PE. rs190390034 (−39 from exon 2) had the strongest association with a predisposing effect with minor allele T [χ² = 7.72, p-value = 0.005; OR = 7.627 (95% CI 1.442, 40.350)] (Figure 2). While we had appropriate power to assess three of six SNPs with single SNP association: rs2287845 [minor allele frequency (MAF) = 0.426, p = 0.038, OR = 1.158, 95% CI = 1.009, 1.331], rs366510 (MAF = 0.426, p = 0.039, OR = 1.158, 95% CI = 1.008, 1.330), and rs2287848 (MAF = 0.426, p = 0.041, OR = 1.155, 95% CI = 1.006, 1.327), the associations of the remaining three SNPs are only suggestive (Table 2).

TABLE 2

Table 2. Six non-coding single nucleotide polymorphism (SNP) have independent allelic association to severe preeclampsia.

A set of alleles spanning from 5' intron proximal to exon 18 to the 5' intron proximal to exon 25 and consisting of 16 SNPs was found to be associated with severe PE. One haplotype showed an association in a protective (frequency in severe PE women and non-PE controls 0.275 and 0.366, respectively, p = 0.044), and another one in a predisposing manner (frequency in severe PE women and non-PE controls = 0.182 and 0.094, respectively, p = 0.011) (Table 3). Mapping the 16 SNPs contributing to the haplotype along the C3 gene showed that the tagging SNP rs2287845 is located 7bp in 3' direction of exon 22, which corresponds to the alpha-chain of the gene product at the edge of MG7 and directly before (5') the CUB domain (Figure 3). The proposed haplotype spans across the middle of the gene from the 3' intron of exon 18 situated at the 5' end of the ANA domain right before the α'NT domain and at 5' direction reaching past the CUB and into the TED domain. Three of the SNPs in the haplotype (rs406514, rs11569450, and rs432823) are located in the introns flanking exon 18, which codes for the α'NT domain. The regulatory microsatellite region was not associated with severe PE.

TABLE 3

Table 3. The region covered by haploblock 2 had two haplotypes of 16 SNPs with suggestive association to preeclampsia.

FIGURE 3

Figure 3. C3 is located in chromosome 19 in the reverse strand. The discovered single nucleotide polymorphism (SNP) associations, their positions, and associating haploblocks are depicted in (A). The associating SNPs and their locations in the domains of the functional C3 are shown in the schematic structure in (B). All associating SNPs are located in introns with the exception of rs183805948, which is located in the PROM. The domains are TED—thioester-containing domain (C3d); CUB—complement C1r/C1s, Uegf, Bmp1; MG1–7—macroglobulin domains 1–7; α'NT—N-terminal region of the cleaved α-chain in the linker domain; Ana—anaphylatoxic region (C3a); LNK—β-strand of the linker domain; PROM—promoter region.

In Silico Functional Analysis

Human Splicing Finder 3.0 found splicing motives to be influenced by four of the associating C3 SNPs: rs190390034, rs2287848, rs366510, and rs2241391. rs366510, rs2241391, and rs190390034 were predicted to create an exonic splicing enhancer (ESE) site within an intron, whereas rs2287848 was found to create a novel ESE site as well as to cause an alteration of an existing ESE site. In RNAsnp, rs190390034 had a p-value of 0.029 indicating a disruption in local RNA folding and rs2241391 had a p-value of 0.169 indicating a likely disruption in local RNA folding. The associating SNPs within the haploblock region (rs2287845, rs366510, and rs2287848) had RNAsnp p-values >0.2 indicating no significant structural change caused by these variants.

The Relative Extended Haplotype Homozygosity

The tagging SNP rs2287845 was associated with an extended haplotype to the 5' direction on the C3 gene suggesting that this structure with tight linkage disequilibrium in the middle of the gene results from a positive selection pressure (Figure 4).

FIGURE 4

Figure 4. The extended haplotype homozygosity (EHH) diagram of (A) women with severe preeclampsia (PE) and (B) non-PE controls using rs2287845 as a focal single nucleotide polymorphism (SNP) representative of the haploblock structure. The ancestral allele curve (in red) is close to 1 for a long stretch in non-PE controls particularly in the 5' upstream direction indicating that the haplotype surrounding the focal SNP (indicated by severed vertical line) is under selection in this cohort. Major allele (ancestral, in red) C, minor allele (derived, in blue) T.

Discussion

In this study, we have described genetic variants of the key complement component C3 in women with severe PE and control women with non-PE pregnancies. We identified a 16-nucleotide haplotype signature in the highly conserved middle region of the maternal C3 gene with three associating SNPs (rs2287845, rs366510, and rs2287848) that could influence susceptibility to severe PE.

The three C3 SNPs that associate with severe PE are located in the area of most intense linkage disequilibrium in the middle of the associating haploblock. The discovery of predisposing SNPs within the haploblock structure supports a possible functional role for the haploblock-encoded gene product. Furthermore, rs190390034 is an intronic variant 38 bp 5' downstream of exon 2 that is in tight linkage disequilibrium with the haploblock structure in the middle of the gene. rs190390034 has a perfect D' score of 100 for the tagging SNP rs2287845 (Figure 2). Among 43 observed SNPs, rs190390034 had the best allelic association with severe PE. While suggestive due to insufficient sample size, the independent association of rs190390034 with severe PE supports the role of the haploblock in severe PE. With the available genotypes for the remaining two of the six associating SNPs within C3 we were underpowered due to insufficient sample size (rs190390034) and rarity of the variants (rs183805948 and rs2241391). Therefore, the latter results should be interpreted with caution.

In support of the causality of changes in the haploblock 2 area of the C3 gene, two of the three associating SNPs within the haploblock described in our study have been implicated in prior studies in other phenotypes. rs2287845*C has been shown to be significantly associated with overall survival and prognosis in patients with early stage non-small cell lung cancer after surgical resection (37). rs366510 is a probable splicing variant that has been linked to asthma and related phenotypes in two independent studies (38, 39). None of these studies looked for a haplotype association within C3. Intronic variants causing C3 splicing mutations have been described in patients with C3 deficiency due to exon deletion from the C3 mRNA (40, 41). The current study is according to our knowledge the first one to show a disease association in C3 with a haplotype-based mechanism instead of the conventional single SNP association. Indeed, it has become apparent that while mutations leading to changes in amino acid sequences are readily detected, more subtle changes in gene regulatory elements are most likely accountable for much of the phenotypic variation we observe in complex diseases (42, 43). Such regulatory features remain cryptic in analyses at the translational and posttranslational level.

In other genes, intronic disease associations with causative splicing defects have been described (44, 45). In the field of reproductive immunology in a set-up similar to our study, an intronic haplotype was found to result in IL-10 secretion changes in women with idiopathic recurrent miscarriage (46). Furthermore, it was recently shown that mechanical nucleosome binding occurs even on top of genes (47). Therefore, it is possible that non-coding SNPs may have an important regulatory role, e.g., by influencing DNA folding.

Complement C3 plays a central role in a successful pregnancy. Inappropriate complement activation may play a role in the initial stages of PE pregnancies contributing to inadequate placentation or placental dysfunction. The anecdotal reports of success of eculizumab in the treatment of a full-blown disease indicate that the complement system is also involved in the later stages of the disease possibly by generating inflammation or tissue damage (22, 48). Problems may emerge if disturbances in the removal of ischemic or injured placental components by complement and phagocytes occur (49). Lack of functioning C3 in mice led to fewer pregnancies and to a higher fetal reabsorption rate, while fetal and placental weights were lower (50). On the other hand, C3-mediated over-activation of the complement system was shown to induce hypertension following placental ischemia in rats (51). Furthermore, complement activation at the feto-maternal interface of Crry^−/− mice that lack a key complement regulator was shown to cause fetal loss. The embryos were rescued when Crry^−/− mice were bred to C3^−/− mice (52). These observations underline the importance of balanced activation and regulation of the complement system for a healthy pregnancy.

Because C3 activation by the C3 convertases requires extensive conformational changes and translocation of the CUB/TED and α'NT domains (17), protein changes caused by variants in the middle of the C3 gene may hinder binding of factor B to C3b causing the C3 convertase to function inefficiently (Figure 3). If the haplotype described here has an effect on C3 function as suggested by its critical location, it is possible that C3 activation in the individuals with the protective haplotype is properly regulated. Thereby, the extravillous trophoblasts (EVTs) invading the maternal tissue during placentation would not encounter a vigorous complement attack. Thus, they could successfully remodel the uterine spiral arteries resulting in a non-PE pregnancy with a healthy blood flow and placental development (9). Concurrently, the predisposing haplotype may result in an increased level of complement activation as indicated by increased factor B (Bb) levels early in the pregnancy. Complement attack could compromise the EVT invasion and consequently the placental function resulting in an increased occurrence of severe PE (11, 53).

C3 promoter activity is dependent on the dinucleotide repeat polymorphism GF100472 such that the longer the CA repeat region, the lower the transcriptional activity of C3 (54). A shorter repeat has shown protective effect against mesial temporal lobe epilepsies and febrile seizures. However, in the present study, we did not find any indication of association of the CA-repeat to severe PE in a small patient cohort.

The functional polymorphism rs2230199 in C3 is known as the slow/fast mutation that influences C3 protein mobility in electrophoresis gels. C3F has been described as predisposing to PE (55), but we did not find any association of C3F with severe PE. Our result concurs with an early study of C3 allotypes that did not find association of C3F with PE (56).

The REHH analysis shows that the haploblock structure in the middle of C3 is tightest for the ancestral rs2287845 allele in non-PE controls suggesting that the structure results from a positive evolutionary selection pressure (Figure 4B). A similar pattern is observed for the ancestral allele of rs2287845 in PE women. However, in severe PE, the haploblock structure disintegrates noticeably sooner than in controls, indicating a loosened force of active selection (Figure 4A). It would follow that due to a possible regulatory feature caused by seemingly benign variants that have been introduced into the middle of C3, the risk for severe PE increases, while haplotypes with the ancestral genotype are protective from severe PE and, accordingly, under stronger positive natural selection.

The heterogeneity of PE is reflected in the comparison of our results to another recent study. Wu et al. found that rs698090 in MASP1 is associated to late-onset but not to early-onset PE, and nominally to severe PE in a Chinese population (57). In our initial genotyping, we did not find allelic or genotypic association of rs608090 to severe PE. It is possible that different complement pathways contribute to early-onset PE and late-onset PE and the mechanism of these varying associations merit more studies.

Targeting gene regulatory effects may provide new opportunities for PE risk assessment and diagnosis, maybe even future drug development (58). The reported results reveal significant differences between PE and healthy pregnant women but the roles of individual SNPs should be considered suggestive and treated with caution. With further studies to confirm our findings, assessing C3 genetic polymorphisms may be developed as a tool to find patients with the highest risk of severe PE.

Ethics Statement

All subjects provided a written informed consent in accordance with the Declaration of Helsinki. Study protocols were approved by the local Ethical Committees, specifically, for the FINNPEC study, ethical approval has been obtained from the Coordinating Ethics Committee, Hospital District of Helsinki and Uusimaa, for the Finnish population based pre-eclampsia cohort was approved by the ethics committee of the Finnish Red Cross Blood Service and by the Ministry of Social Affairs and Health and for Southern Finland pre-eclampsia study was approved by the Ethics Committee of the Department of Obstetrics and Gynaecology at Helsinki University Central Hospital.

Author Contributions

AIL designed the study with SM and HL. AIL and TK designed part of the primers for sequencing, AIL supervised laboratory work, sequenced all samples from patients and some controls, read the results, conducted association analyses of the C3 sequencing project, performed in silico analyses, and drafted the manuscript. TK conducted laboratory work for the regulatory regions branch of the study, analyzed these data, and participated in main data analysis. VH designed the complement genotyping Sequenom chip. PO analyzed the genotyping data. LK performed and interpreted the REHH analysis with AIL and PS. PV described patient cohorts. HL, SH, EK, JK, KK, and AP form the board of the FINNPEC cohort and are responsible for the clinical data and biological samples used in this study. LH provided samples and is with HL responsible for the clinical data of the Finnish population-based preeclampsia cohort. HL provided the samples and clinical data of The Southern Finland preeclampsia cohort. All the authors collaborated in drafting the manuscript and accepted the final version of the manuscript.

Conflict of Interest Statement

This research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Jaana Oikkonen Ph.D., who provided bioinformatics assistance in mapping the promoter and regulatory regions of C3. We thank Pekka Ellonen for designing part of the sequencing primers and Hanna Jarva for her expertise in complementology. Tia Aalto-Viljakainen, MD, described the replication dataset for which we are grateful. Riitta Lehtinen, Susanna Mehtälä, and Eija Kortelainen are thanked for their technical expertise. We are indebted to all the women, who participated in these studies.

Funding

This study was supported by Maud Kuistila memorial foundation (AIL), Oskar Öflund’s Foundation (AIL), Alfred Kordelin Foundation (AIL), Jane and Aatos Erkko Foundation (HL), The Academy of Finland (121196, 134957, and 278941, HL), Päivikki and Sakari Sohlberg Foundation (HL), Sigrid Jusélius Foundation (SM), Doctoral Programme in Biomedicine, Clinical Graduate Program (TK), Research Foundation of the University of Helsinki (TK) and Biomedicum Helsinki Foundation (TK). Finska Läkaresällskapet, Emil Aaltonen Foundation, Finnish Medical Foundation, University of Helsinki Funds, Special State Subsidy for Health Research (EVO funding), Novo Nordisk Foundation, Signe and Ane Gyllenberg Foundation, and Foundation for Pediatric Research contributed to the FINNPEC study. JK is a recipient of The Royal Society Wolfson Research Excellence award. Funders had no influence in methodology or data interpretation of our work.

Supplementary Material

The Supplementary Material for this article can be found online at http://journal.frontiersin.org/article/10.3389/fimmu.2017.00589/full#supplementary-material.

References

1. Roberts JM, Pearson G, Cutler J, Lindheimer M; NHLBI Working Group on Research on Hypertension During Pregnancy. Summary of the NHLBI working group on research on hypertension during pregnancy. Hypertension (2003) 41:437–45. doi: 10.1161/01.HYP.0000054981.03589.E9

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Andraweera PH, Gatford KL, Dekker GA, Leemaqz S, Jayasekara RW, Dissanayake VH, et al. The INSR rs2059806 single nucleotide polymorphism, a genetic risk factor for vascular and metabolic disease, associates with pre-eclampsia. Reprod Biomed Online (2017) 34(4):392–8. doi:10.1016/j.rbmo.2017.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Thomsen LC, McCarthy NS, Melton PE, Cadby G, Austgulen R, Nygard OK, et al. The antihypertensive MTHFR gene polymorphism rs17367504-G is a possible novel protective locus for preeclampsia. J Hypertens (2017) 35:132–9. doi:10.1097/HJH.0000000000001131

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Wan JP, Zhao H, Li T, Li CZ, Wang XT, Chen ZJ. The common variant rs11646213 is associated with preeclampsia in Han Chinese women. PLoS One (2013) 8:e71202. doi:10.1371/journal.pone.0071202

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Staff AC, Redman CW, Williams D, Leeson P, Moe K, Thilaganathan B, et al. Pregnancy and long-term maternal cardiovascular health: progress through harmonization of research cohorts and biobanks. Hypertension (2016) 67:251–60. doi:10.1161/HYPERTENSIONAHA.115.06357

CrossRef Full Text | Google Scholar

6. Best LG, Nadeau M, Davis K, Lamb F, Bercier S, Anderson CM. Genetic variants, immune function, and risk of pre-eclampsia among American Indians. Am J Reprod Immunol (2012) 67:152–9. doi:10.1111/j.1600-0897.2011.01076.x

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Quach K, Grover SA, Kenigsberg S, Librach CL. A combination of single nucleotide polymorphisms in the 3'untranslated region of HLA-G is associated with preeclampsia. Hum Immunol (2014) 75:1163–70. doi:10.1016/j.humimm.2014.10.009

CrossRef Full Text | Google Scholar

8. Pijnenborg R, Vercruysse L, Hanssens M. The uterine spiral arteries in human pregnancy: facts and controversies. Placenta (2006) 27:939–58. doi:10.1016/j.placenta.2005.12.006

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Lokki AI, Heikkinen-Eloranta J, Jarva H, Saisto T, Lokki ML, Laivuori H, et al. Complement activation and regulation in preeclamptic placenta. Front Immunol (2014) 5:312. doi:10.3389/fimmu.2014.00312

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Tang MX, Zhang YH, Hu L, Kwak-Kim J, Liao AH. CD14++ CD16+ HLA-DR+ monocytes in peripheral blood are quantitatively correlated with the severity of pre-eclampsia. Am J Reprod Immunol (2015) 74:116–22. doi:10.1111/aji.12389

CrossRef Full Text | Google Scholar

11. Hoffman MC, Rumer KK, Kramer A, Lynch AM, Winn VD. Maternal and fetal alternative complement pathway activation in early severe preeclampsia. Am J Reprod Immunol (2014) 71:55–60. doi:10.1111/aji.12162

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Boij R, Svensson J, Nilsson-Ekdahl K, Sandholm K, Lindahl TL, Palonek E, et al. Biomarkers of coagulation, inflammation, and angiogenesis are independently associated with preeclampsia. Am J Reprod Immunol (2012) 68:258–70. doi:10.1111/j.1600-0897.2012.01158.x

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Non-aka M, Kimura A. Genomic view of the evolution of the complement system. Immunogenetics (2006) 58:701–13. doi:10.1007/s00251-006-0142-1

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Non-aka M, Yoshizaki F. Evolution of the complement system. Mol Immunol (2004) 40:897–902. doi:10.1016/j.molimm.2003.10.009

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Non-aka M. Evolution of the complement system. Subcell Biochem (2014) 80:31–43. doi:10.1007/978-94-017-8881-6_3

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Meri S. Self-non-self discrimination by the complement system. FEBS Lett (2016) 590:2418–34. doi:10.1002/1873-3468.12284

CrossRef Full Text | Google Scholar

17. Janssen BJ, Christodoulidou A, McCarthy A, Lambris JD, Gros P. Structure of C3b reveals conformational changes that underlie complement activity. Nature (2006) 444:213–6. doi:10.1038/nature05172

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Forneris F, Ricklin D, Wu J, Tzekou A, Wallace RS, Lambris JD, et al. Structures of C3b in complex with factors B and D give insight into complement convertase formation. Science (2010) 330:1816–20. doi:10.1126/science.1195821

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Janssen BJ, Gros P. Structural insights into the central complement component C3. Mol Immunol (2007) 44:3–10. doi:10.1016/j.molimm.2006.06.017

CrossRef Full Text | Google Scholar

20. Taniguchi-Sidle A, Isenman DE. Interactions of human complement component C3 with factor B and with complement receptors type 1 (CR1, CD35) and type 3 (CR3, CD11b/CD18) involve an acidic sequence at the N-terminus of C3 alpha’-chain. J Immunol (1994) 153:5285–302.

PubMed Abstract | Google Scholar

21. Mao D, Wu X, Deppong C, Friend LD, Dolecki G, Nelson DM, et al. Negligible role of antibodies and C5 in pregnancy loss associated exclusively with C3-dependent mechanisms through complement alternative pathway. Immunity (2003) 19:813–22. doi:10.1016/S1074-7613(03)00321-2

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Burwick RM, Feinberg BB. Eculizumab for the treatment of preeclampsia/HELLP syndrome. Placenta (2013) 34:201–3. doi:10.1016/j.placenta.2012.11.014

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Lokki AI, Klemetti MM, Heino S, Hiltunen L, Heinonen S, Laivuori H. Association of the rs1424954 polymorphism of the ACVR2A gene with the risk of pre-eclampsia is not replicated in a Finnish study population. BMC Res Notes (2011) 4:545. doi:10.1186/1756-0500-4-545

CrossRef Full Text | Google Scholar

24. Jääskeläinen T, Heinonen S, Kajantie E, Kere J, Kivinen K, Pouta A, et al. Cohort profile: the Finnish Genetics of Pre-Eclampsia Consortium (FINNPEC). BMJ Open (2016) 6:e013148. doi:10.1136/bmjopen-2016-013148

PubMed Abstract | CrossRef Full Text | Google Scholar

25. ACOG Committee on Obstetric Practice. ACOG practice bulletin. Diagnosis and management of preeclampsia and eclampsia. Number 33, January 2002. American College of Obstetricians and Gynecologists. Int J Gynaecol Obstet (2002) 99:159–67. doi:10.1016/S0029-7844(01)01747-1

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Pihkala J, Hakala T, Voutilainen P, Raivio K. Characteristic of recent fetal growth curves in Finland. Duodecim (1989) 105:1540–6.

Google Scholar

27. Holmberg V, Onkamo P, Lahtela E, Lahermo P, Bedu-Addo G, Mockenhaupt FP, et al. Mutations of complement lectin pathway genes MBL2 and MASP2 associated with placental malaria. Malar J (2012) 11:61. doi:10.1186/1475-2875-11-61

PubMed Abstract | CrossRef Full Text | Google Scholar

28. International HapMap Consortium. The International HapMap Project. Nature (2003) 426:789–96. doi:10.1038/nature02168

CrossRef Full Text | Google Scholar

29. Desmet FO, Hamroun D, Lalande M, Collod-Beroud G, Claustres M, Beroud C. Human Splicing Finder: an online bioinformatics tool to predict splicing signals. Nucleic Acids Res (2009) 37:e67. doi:10.1093/nar/gkp215

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Sabarinathan R, Tafer H, Seemann SE, Hofacker IL, Stadler PF, Gorodkin J. The RNAsnp web server: predicting SNP effects on local RNA secondary structure. Nucleic Acids Res (2013) 41:W475–9. doi:10.1093/nar/gkt291

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Sabarinathan R, Tafer H, Seemann SE, Hofacker IL, Stadler PF, Gorodkin J. RNAsnp: efficient detection of local RNA secondary structure changes induced by SNPs. Hum Mutat (2013) 34:546–56. doi:10.1002/humu.22273

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Gautier M, Vitalis R. rehh: an R package to detect footprints of selection in genome-wide SNP data from haplotype structure. Bioinformatics (2012) 28:1176–7. doi:10.1093/bioinformatics/bts115

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Stephens M, Donnelly P. A comparison of Bayesian methods for haplotype reconstruction from population genotype data. Am J Hum Genet (2003) 73:1162–9. doi:10.1086/379378

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Scheet P, Stephens M. A fast and flexible statistical model for large-scale population genotype data: applications to inferring missing genotypes and haplotypic phase. Am J Hum Genet (2006) 78:629–44. doi:10.1086/502802

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet (2007) 81:559–75. doi:10.1086/519795

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics (2005) 21:263–5. doi:10.1093/bioinformatics/bth457

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Lee SY, Choi JE, Jeon HS, Choi YY, Lee WK, Lee EB, et al. A panel of genetic polymorphism for the prediction of prognosis in patients with early stage non-small cell lung cancer after surgical resection. PLoS One (2015) 10:e0140216. doi:10.1371/journal.pone.0140216

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Barnes KC, Grant AV, Baltadzhieva D, Zhang S, Berg T, Shao L, et al. Variants in the gene encoding C3 are associated with asthma and related phenotypes among African Caribbean families. Genes Immun (2006) 7:27–35. doi:10.1038/sj.gene.6364267

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Purwar R, Langer K, Werfel T. Polymorphisms within the C3 gene are associated with specific IgE levels to common allergens and super-antigens among atopic dermatitis patients. Exp Dermatol (2009) 18:30–4. doi:10.1111/j.1600-0625.2008.00759.x

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Silva KR, Fraga TR, Lucatelli JF, Grumach AS, Isaac L. Skipping of exon 27 in C3 gene compromises TED domain and results in complete human C3 deficiency. Immunobiology (2016) 221:641–9. doi:10.1016/j.imbio.2016.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Tsukamoto H, Horiuchi T, Kokuba H, Nagae S, Nishizaka H, Sawabe T, et al. Molecular analysis of a novel hereditary C3 deficiency with systemic lupus erythematosus. Biochem Biophys Res Commun (2005) 330:298–304. doi:10.1016/j.bbrc.2005.02.159

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Brest P, Lapaquette P, Souidi M, Lebrigand K, Cesaro A, Vouret-Craviari V, et al. A synonymous variant in IRGM alters a binding site for miR-196 and causes deregulation of IRGM-dependent xenophagy in Crohn’s disease. Nat Genet (2011) 43:242–5. doi:10.1038/ng.762

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Acedo A, Sanz DJ, Duran M, Infante M, Perez-Cabornero L, Miner C, et al. Comprehensive splicing functional analysis of DNA variants of the BRCA2 gene by hybrid minigenes. Breast Cancer Res (2012) 14:R87. doi:10.1186/bcr3202

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Nakamura K, Du L, Tunuguntla R, Fike F, Cavalieri S, Morio T, et al. Functional characterization and targeted correction of ATM mutations identified in Japanese patients with ataxia-telangiectasia. Hum Mutat (2012) 33:198–208. doi:10.1002/humu.21632

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Mancini C, Vaula G, Scalzitti L, Cavalieri S, Bertini E, Aiello C, et al. Megalencephalic leukoencephalopathy with subcortical cysts type 1 (MLC1) due to a homozygous deep intronic splicing mutation (c.895-226T > G) abrogated in vitro using an antisense morpholino oligonucleotide. Neurogenetics (2012) 13:205–14. doi:10.1007/s10048-012-0331-z

CrossRef Full Text | Google Scholar

46. Qaddourah RH, Magdoud K, Saldanha FL, Mahmood N, Mustafa FE, Mahjoub T, et al. IL-10 gene promoter and intron polymorphisms and changes in IL-10 secretion in women with idiopathic recurrent miscarriage. Hum Reprod (2014) 29:1025–34. doi:10.1093/humrep/deu043

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Eslami-Mossallam B, Schram RD, Tompitak M, van Noort J, Schiessel H. Multiplexing genetic and nucleosome positioning codes: a computational approach. PLoS One (2016) 11:e0156905. doi:10.1371/journal.pone.0156905

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Nadeau-Vallee M, Obari D, Palacios J, Brien ME, Duval C, Chemtob S, et al. Sterile inflammation and pregnancy complications: a review. Reproduction (2016) 6:R277–92. doi:10.1530/REP-16-0453

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Reichhardt MP, Jarva H, Lokki AI, Laivuori H; FINNPEC Study Group, Vuorela P, et al. The salivary scavenger and agglutinin (SALSA) in healthy and complicated pregnancy. PLoS One (2016) 11:e0147867. doi:10.1371/journal.pone.0147867

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Chow W, Lee Y, Wong P, Chung M, Lee K, Yeung WS. Complement 3 deficiency impairs early pregnancy in mice. Mol Reprod Dev (2009) 76:647–55. doi:10.1002/mrd.21013

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Lillegard KE, Johnson AC, Lojovich SJ, Bauer AJ, Marsh HC, Gilbert JS, et al. Complement activation is critical for placental ischemia-induced hypertension in the rat. Mol Immunol (2013) 56:91–7. doi:10.1016/j.molimm.2013.04.009

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Xu C, Mao D, Holers VM, Palanca B, Cheng AM, Molina H. A critical role for murine complement regulator crry in fetomaternal tolerance. Science (2000) 287:498–501. doi:10.1126/science.287.5452.498

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Velickovic I, Dalloul M, Wong KA, Bakare O, Schweis F, Garala M, et al. Complement factor B activation in patients with preeclampsia. J Reprod Immunol (2015) 109:94–100. doi:10.1016/j.jri.2014.12.002

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Jamali S, Salzmann A, Perroud N, Ponsole-Lenfant M, Cillario J, Roll P, et al. Functional variant in complement C3 gene promoter and genetic susceptibility to temporal lobe epilepsy and febrile seizures. PLoS One (2010) 5:e12740. doi:10.1371/journal.pone.0012740

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Rhim MS, Meddeb S, Kaabia O, Jalloul M, Sakouhi M, Jrzad BB, et al. C3F gene mutation is involved in the susceptibility to pre-eclampsia. Arch Gynecol Obstet (2015) 291:1023–7. doi:10.1007/s00404-014-3515-y

CrossRef Full Text | Google Scholar

56. Cooper DW, Deane EM, Marshall P, Gallery ED. C3 allotypes in pregnancy hypertension and eclampsia. Hum Hered (1988) 38:52–5. doi:10.1159/000153755

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Wu W, Yang H, Feng Y, Zhang P, Li S, Wang X, et al. Polymorphisms in complement genes and risk of preeclampsia in Taiyuan, China. Inflamm Res (2016) 65:837–45. doi:10.1007/s00011-016-0968-4

CrossRef Full Text | Google Scholar

58. Dai Y, Qiu Z, Diao Z, Shen L, Xue P, Sun H, et al. MicroRNA-155 inhibits proliferation and migration of human extravillous trophoblast derived HTR-8/SVneo cells via down-regulating cyclin D1. Placenta (2012) 33:824–9. doi:10.1016/j.placenta.2012.07.012

PubMed Abstract | CrossRef Full Text | Google Scholar