Abstract
Pathologic thrombosis is a major cause of mortality. Hemolytic-uremic syndrome (HUS) features episodes of small vessel thrombosis resulting in microangiopathic hemolytic anemia, thrombocytopenia and renal failure1. Atypical HUS (aHUS) can result from genetic or autoimmune factors2 that lead to pathologic complement cascade activation3. By exome sequencing we identify recessive mutations in DGKE (diacylglycerol kinase epsilon) that co-segregate with aHUS in 9 unrelated kindreds, defining a distinctive Mendelian disease. Affected patients present with aHUS before age 1, have persistent hypertension, hematuria and proteinuria (sometimes nephrotic range), and develop chronic kidney disease with age. DGKE is found in endothelium, platelets, and podocytes. Arachidonic acid-containing diacylglycerols (DAG) activate protein kinase C, which promotes thrombosis. DGKE normally inactivates DAG signaling. We infer that loss of DGKE function results in a pro-thrombotic state. These findings identify a new mechanism of pathologic thrombosis and kidney failure and have immediate implications for treatment of aHUS patients.
Mendelian forms of atypical HUS (aHUS) have implicated mutations in genes of the complement cascade, including complement factors B (CFB), H (CFH), and I (CFI), complement component 3 (C3), membrane cofactor protein (MCP) and thrombomodulin (THBD)2. All of these mutations result in unrestricted complement activation3 Most are transmitted as autosomal dominant traits with markedly reduced penetrance; only recessive mutations in CFH and MCP show apparently high penetrance2. Nonetheless, nearly half of aHUS patients without secondary causes have no discernable genetic or autoimmune abnormality4.
We studied two unrelated families (kindreds 1 and 2), each with two siblings diagnosed with aHUS in infancy and unaffected unrelated parents. There were no pathogenic mutations in known aHUS genes nor anti-CFH antibodies (Supplementary Table 1). All four presented between 4 and 8 months of age with microangiopathic hemolytic anemia, thrombocytopenia and acute renal failure (Table 1 and Supplementary Table 2). Three had renal biopsies before age 3, all with pathology demonstrating chronic thrombotic microangiopathy (Table 1 and Fig. 1a-d). We performed exome sequencing of these 4 affected subjects (Supplementary Table 3). High quality variations from the reference sequence were called, their impact on encoded proteins determined and allele frequencies estimated.
Table 1.
Subject ID |
Age at Dx (yr), Gender |
Hgb nadir nl>10.3 g/dL) |
Evidence of MAb |
PLT (nl<220 x109/L) |
Cr (nl<0.2 mg/dL) |
Dialysis at onset |
Dx of first renal Bx (age, yr) |
DGKE mutation 1 |
DGKE mutation 2 |
---|---|---|---|---|---|---|---|---|---|
1-3 | 0.7, M | 7.2 | Yes | 50 | 2.84 | Yes | cTMA (2) | p.Trp322* | p.Trp322* |
1-4 | 0.3, F | 5.7 | Yes | 36 | 1.89 | Yes | nd | p.Trp322* | p.Trp322* |
2-5c | 0.6, F | 7.3 | Yes | 35 | 5.32 | No | cTMA (1) | p.Arg63Pro | p.Val163Serfs*3 |
2-7c | 0.3, F | 6.8 | Yes | 168 | 1.35 | No | cTMA (1.5) | p.Arg63Pro | p.Val163Serfs*3 |
3-3c | 0.5, F | 8.4 | Yes | 132 | 0.62 | No | cTMA (0.7) | p.Trp322* | p.Ser11* |
4-1c | 0.3, F | 3.7 | Yes | 390 | 8.51 | Yes | cTMA (1) | p.Trp322* | p.Trp322* |
5-3c | 0.6, M | 7.1 | Yes | 88 | 5.09 | Yes | cTMA (1) | p.Trp322* | p.Trp158Leufs*8 |
6-3 | 0.5, F | 9.0 | Yes | 99 | 2.52 | Yes | cTMA (1) | p.Gln334* | p.Gln334* |
7-3 | 0.9, M | 4.9 | Yes | 125 | 6.79 | Yes | cTMA (2) | p.IVS5-1 | p.IVS5-1 |
8-3 | 0.7, M | 5.0 | Yes | 214 | 7.19 | Yes | cTMA (1e) | p.Trp322* | p.Trp322* |
9-3 | 0.3, M | 6.4 | Yes | 33 | 4.21 | Yes | cTMA (7) | p.Arg273Pro | p.Arg273Pro |
9-4 | 0.8, M | 8.2 | nd | 57 | 0.6 | No | nd | p.Arg273Pro | p.Arg273Pro |
9-6 | 0.3, F | 7.3 | Yes | 32 | 2.21 | Yes | nd | p.Arg273Pro | p.Arg273Pro |
Bx, biopsy; Cr, creatinine; Dx, diagnosis; cTMA, chronic thrombotic microangiopathy; Hgb, hemoglobin; MA, microangiopathy; nd, not done; nl, normal; PLT, platelet count.
Reference values for infants <1 year of age32.
Microangiopathy diagnosed by high lactate dehydrogenase, low haptoglobin and/or schistocytes on blood smear (details in Supplementary Table 2).
Kindreds previously reported33.
Patient had many normal platelet counts with acute aHUS.
Global renal cortical ischemia and malignant nephroangiosclerosis noted in both kidneys.
We posited autosomal recessive transmission in these families and sought genes with rare homozygous or compound heterozygous variants (minor allele frequency < 1%, and homozygous/compound heterozygous genotypes not previously seen in databases) that were shared by both affected subjects (Supplementary Table 4). In kindred 1, there was a single novel homozygous variant shared by both affected subjects, and there was one novel shared compound heterozygous genotype in kindred 2. These novel genotypes occurred in the same gene, diacylglycerol kinase epsilon (DGKE). The former was a homozygous premature termination codon (p.Trp322*) while the latter genotype was a compound heterozygote for a frameshift (p.Val163Serfs*3) and a missense mutation (p.Arg63Pro); 3 unaffected siblings had zero or one of these variants (Fig. 2, Table 1 and Supplementary Fig. 1a).
To extend these findings, we sequenced DGKE in 47 additional unrelated probands with pediatric-onset aHUS and 36 adult-onset aHUS probands in whom mutation in known aHUS-associated genes or anti-CFH antibodies were not found (Supplementary Table 1). The results identified 6 additional index cases, harboring rare homozygous or compound heterozygous DGKE variants, all in pediatric-onset cases (Fig. 2, Table 1, and Supplementary Fig. 1a). Parental samples, available for all but one kindred, were heterozygous for one of the mutations with the exception of kindred 5, in which one mutation was apparently de novo. Additionally, kindred 9, independently ascertained in Germany with three affected subjects, showed complete linkage to the DGKE locus (LOD score 2.53; Supplementary Fig. 1b) and sequencing of all exons in the interval identified a homozygous DGKE p.Arg273Pro mutation (Fig. 2 and Supplementary Fig. 1a). These 9 patients all met clinical criteria for aHUS at presentation (Table 1 and Supplementary Table 2). Six had renal biopsies before age 2, all read as chronic thrombotic microangiopathy (Table 1 and Fig. 1e-g).
Collectively, the rare DGKE variants found in the 9 kindreds included 3 different premature termination codons, 2 frameshift mutations, 1 splice donor site mutation and two missense mutations that occur at conserved positions (Fig. 2 and Supplementary Fig. 1c). Only one of these variants, p.Trp322*, was previously seen among 8,475 subjects from NHLBI or Yale exome databases; this variant was heterozygous in two people of European ancestry. p.Trp322* was found in five apparently unrelated aHUS subjects of European ancestry, and was homozygous in three. These three subjects shared an identical and extremely rare haplotype spanning no more than 400 kb at the DGKE locus (Supplementary Fig. 2 and Supplementary Table 5). This indicates a common ancestry for the mutation in each family, with the last common ancestor estimated to have occurred 53 generations ago (95% confidence interval, 33-73; Supplementary Fig. 3). The remote shared ancestry of the mutation is consistent with these 3 families not being closely related.
Twenty-two percent of siblings of index cases in these families (4/18) had aHUS, consistent with recessive transmission with high penetrance. Moreover, rare DGKE variants precisely cosegregate with aHUS in these families, yielding a LOD score of 8.9 (likelihood ratio of 7.9 × 108 in favor of linkage) for complete linkage of these variants to aHUS under a recessive model with rare phenocopies and high penetrance. Similarly, we found no homozygous/compound heterozygous damaging variants (premature termination, frameshift or splice site variants) among 8,475 subjects in NHLBI and Yale exome databases (Supplementary Fig. 4). The association of rare recessive genotypes with aHUS is extremely strong (P = 2 × 10−16, Fisher’s exact test). Together, these genetic findings unequivocally establish recessive loss of function mutations in DGKE as the cause of aHUS in these families.
Patients with DGKE mutations all presented with aHUS in the first year of life (mean 0.5 years, range 0.3 – 0.9; Table 1). DGKE mutations were found in 9 of 22 pediatric aHUS patients with disease onset < 1 year of age, and 0 of 28 diagnosed after age 1 (P = 2 × 10−4). DGKE was a frequent cause of aHUS in the first year of life (13 of 49 aHUS cases, 27%) and accounted for 50% of familial disease in this age group (3 of 6 kindreds). This uniformly early age of onset defines a distinct subgroup of aHUS (Fig. 3a). The clinical course in DGKE-mutant subjects featured relapsing episodes of HUS before age 5 (Fig. 3b).
Abnormal complement activation is a feature of all previously described forms of aHUS3. Because DGKE encodes an intracellular enzyme, it is not obvious that complement activation plays a role in DGKE-mutant patients. Detailed assessment of the complement system revealed no compelling abnormality in any subject (Supplementary Table 1). Moreover, two patients with DGKE mutations had HUS relapses while on anti-complement therapy (eculizumab and fresh frozen plasma infusions, respectively) which are believed efficacious in patients with complement defects5 (Supplementary Fig. 5 and Supplementary Table 2).
Among patients whose renal function recovered after the onset of aHUS, hypertension, microhematuria and proteinuria persisted in all but one (subject 5-3). In contrast, most patients with other aHUS subtypes have no residual renal abnormalities between episodes. Progression to chronic kidney disease (CKD) stages 4 and 5 was common by the second decade of life, long after the last acute episode of HUS (Fig. 3c). Interesingly, three patients developed nephrotic syndrome 3-5 years after disease onset, a very rare event in other forms of HUS. Repeat renal biopsies in these patients were read as chronic thrombotic microangiopathy (Supplementary Table 2). This is not likely a genotypic effect, because two affected siblings have not developed nephrotic syndrome. Treating physicians did not report extra-renal phenotypes.
Three subjects received cadaveric renal transplantation at 2, 19 and 21 years using standard pre- and post-transplantation protocols (Supplementary Table 2). Two allografts have survived 2 and 4 years, while the other failed after 6 years due to chronic rejection. Importantly, there have been no HUS recurrences post-transplant. This contrasts with aHUS with defects in the soluble complement cascade in which recurrent HUS is very common and graft failure almost invariably occurs without anti-complement therapy2.
DGKE was first cloned from a human endothelial cell line6; studies have noted high expression in testes and little expression elsewhere6. We examined DGKE expression by western blotting in human endothelial cells and platelets, the major cell types involved in thrombosis. We found DGKE expression in both (Fig. 4a-b). Moreover, staining of normal human kidney revealed DGKE in endothelium of glomerular capillaries and podocytes (Fig. 4c and Supplementary Fig. 7). Similar results were seen in rat, with colocalization with WT1, a podocyte marker (Supplementary Fig. 8). The specificity of this staining was established by parallel study of a renal biopsy of subject 2-7, a compound heterozygote for a frameshift mutation and a missense mutation in a beta-pleated sheet in the first DGKE C1 domain. Virtually no DGKE was seen in glomeruli of this patient (Fig. 4d and Supplementary Fig. 7a), while other antisera, such as anti-CD34, show normal staining (Supplementary Fig. 7c).
These findings establish recessive loss of function mutations in DGKE as a frequent cause of aHUS in the first year of life. Consanguinity, recurrence among siblings, persistence of hypertension, microhematuria and proteinuria (especially if in the nephrotic range) along with the absence of complement abnormalities, suggest the diagnosis of DKGE nephropathy. The molecular diagnosis of DGKE mutation should be straighforward.
DGKE is the first gene implicated in aHUS that is not an integral component of the complement cascade, raising the question of the pathophysiologic mechanism. DGKE preferentially phosphorylates arachidonic acid-containing diacylglycerol (AADAG) to the corresponding phosphatidic acid (PA)7. AADAG is a major signaling molecule of the diacylglycerol family, produced by hydrolysis of phosphatidylinositol 4,5-bisphosphate (PIP2) by phospholipase C (PLC) in response to cell surface receptor signaling8 (Supplementary Fig. 9). AADAGs activate protein kinase C (PKC)9. In endothelial cells, PKC increases production of various pro-thrombotic (von Willebrand factor10, plasminogen activator inhibitor-111, platelet-activating factor12, and tissue factor13) and antithrombotic factors such as tissue-type plasminogen activator14 (Supplementary Fig. 10a). The factors determining the balance between these pro- and anti-thrombotic factors are poorly understood. AADAG-dependent PKC signaling also drives thrombin-induced platelet activation15 (Supplementary Fig. 10b). Phosphorylation of AADAG to phosphatidic acid by DGKE terminates AADAG signaling. It is therefore plausible that loss of DGKE results in sustained AADAG signaling16, resulting in a prothrombotic state. This mechanism is supported by experiments with R5902217, a small molecule inhibitor of several DGK’s including DGKE. This produces platelet activation18 and inhibition of endothelial prostacyclin production19.
A similar mechanism may pertain to podocytes, where DAGs modifies slit diaphragm function20, including inducing endocytosis of nephrin21, an effect that could contribute to proteinuria and kidney failure. Additionally, VEGF signaling is essential for podocyte and renal endothelial cell survival, and loss of signaling in renal endothelium produces thrombotic microangiopathy22,23. Importantly, PKC-dependent downregulation of VEGFR2 has been reported in podocyte24 and endothelial cells25. This interaction suggests a potential mechanism for the pronounced renal effects of this form of microangiopathy (Supplementary Fig. 10a). Further work will be required to determine the detailed biochemical mechanism(s) linking DGKE deficiency to aHUS.
The specific triggers that account for the episodic nature of acute attacks of aHUS are poorly understood. Dgke-null mice were not reported to have a thrombotic phenotype, perhaps because such an inciting factor was missing26 or intrinsic species differences.
The universal findings of hypertension, microhematuria and proteinuria, and the unique finding of nephrotic syndrome among aHUS patients with DGKE mutations suggests they might play a role in other kidney diseases with similar glomerular phenotypes, such as systemic lupus erythematosis-associated glomerulonephritis27, severe preeclampsia/HELLP syndrome28, or membranoproliferative glomerulonephritis (MPGN)29. It is presently unclear whether the hypertension, microhematuria and proteinuria found in patients with DGKE mutations require prior acute episodes of HUS or whether it could occur in their absence.
In this regard, it is of great interest that during review of this manuscript a paper reporting 3 families with an early-onset MPGN-like syndrome featuring proteinuria and renal failure were described with recessive DGKE mutations. None of the affected patients were noted to have acute episodes suggesting aHUS, however histologic features of glomerular microangiopathy were noted. These authors also presented experimental evidence that loss of DGKE in podocytes might contribute to proteinuria and renal damage via aberrant activation of TRPC630.
Lastly, these findings have immediate implications for the treatment of aHUS. Current guidelines recommend that patients with aHUS be treated with eculizumab, an anti-C5 antibody that inhibits the complement cascade, or plasma therapy5. The absence of evidence linking DGKE deficiency to the complement cascade and relapses of acute HUS in these patients while receiving these therapies suggests they may not benefit DGKE patients. Moreover, unlike patients with soluble complement defects, it appears that renal transplantation can be efficacious and safe in patients with aHUS due to DGKE mutations, underscoring the importance of diagnosing DGKE nephropathy.
URLs. BLAST, http://blast.ncbi.nlm.nih.gov/Blast.cgi; ClustalW2, http://www.ebi.ac.uk/Tools/msa/clustalw2/; NCBI protein, http://www.ncbi.nlm.nih.gov/protein; NHLBI GO Exome Sequencing Project (ESP), https://esp.gs.washington.edu/drupal/; UCSC Genome browser, http://genome.ucsc.edu/; dbSNP, http://www.ncbi.nlm.nih.gov/snp.
METHODS
Methods and any associated references are available in the online version of the paper.
Accession codes
The DGKE variants described were deposited in dbSNP under batch accession number 1058996. mRNA and protein sequences are available at NCBI under the following accession numbers: human DGKE, NM_003647.2 and NP_003638.1, cow DGKE, NP_001179859.1, mouse Dgke, NP_062378.1, Xenopus dgke, NP_001087580.1, zebrafish dgke NP_001165699.1, fruitfly dgke, NP_725228.1, and worm dgk-2, NP_001024679.1. The reference protein sequences for pig (NP_001161117; 428 amino acids) and rat (NP_001034430; 407 amino acids) were ~70% shorter than expected. Full-length pig DGKE protein sequence (564 amino acids) was derived by translating the cDNA AK400222, while full-length rat DGKE protein (567 amino acids) was reassembled from the mRNA NM_001039341 by removing a short retained intron.
Detailed methods
aHUS kindreds
The French atypical HUS (aHUS) cohort34 includes 139 patients from 131 unrelated kindreds with pediatric-onset aHUS and 36 unrelated subjects with adult-onset aHUS. These patients were recruited from 24 participating pediatric nephrology centers from France, and 1 from Belgium; 36 adult nephrology centers also enroll patients. This cohort includes 55 patients from 49 kindreds without mutation in any of the known aHUS genes or anti-CFH antibodies. One additional kindred including 3 affected subjects without mutations in known aHUS genes was ascertained in Dubai and Germany. IRB protocols were approved at all sites involved in the study and all subjects studied provided informed consent.
For all patients, relevant data were abstracted from detailed medical records. aHUS cases were defined using the classic diagnostic triad: microangiopathic hemolytic anemia (hemoglobin < 10 g/dL, with lactate dehydrogenase level > 250 IU/L, haptoglobin < 0.4 g/L or presence of schizocytes on blood smear), thrombocytopenia (platelets count lower than the lower limit of normal for age) and renal failure (serum creatinine level higher than the upper limit of normal for age). If one of the elements of the diagnostic triad was missing, clear evidence of thrombotic microangiopathic (TMA) lesions observed on a kidney biopsy was necessary to substantiate an aHUS diagnosis. Remission was defined by normalization of platelet and LDH levels, and relapse was defined by recurrence of microangiopathic hemolytic anemia, thrombocytopenia and/or a 25% increase in serum creatinine after at least two weeks remission. We specifically excluded patients with HUS caused by Shiga-toxin producing Escherichia coli and patients with secondary causes of HUS, such as drug exposure, autoimmune diseases, infections (Streptococcus pneumoniae, or human immunodeficiency virus), bone marrow or solid organ transplantation, and cobalamin deficiency.
Assessment of aHUS susceptibility factors (proteins and genes)
Assessment of complement system function was done at the complement laboratory at Hôpital Européen Georges-Pompidou (Paris, France), the National reference center for the evaluation of complement disorders. Blood samples from aHUS patients were collected for investigations of complement function and genetic analyses. Similar studies were performed in parallel on samples from 200 unrelated healthy French subjects to generate reference levels and validate any rare variants identified. Extensive investigations of the complement system were performed as described35. Briefly, the plasma concentrations of factors H and I (CFH, CFI) were measured by ELISA while levels of factors 3, 4 and B (C3, C4 and CFB) were measured by nephelometry. MCP surface expression was analyzed on granulocytes using anti-MCP phycoerythrin-conjugated antibodies (Serotec, UK). All patients were screened for anti-CFH antibodies and ADAMTS13 deficiency. All coding sequences of the CFH, CFI, MCP, C3, CFB and thrombomodulin (THBD) genes were sequenced as described35. Screening for unequal crossing over between homologous genes at the CFH, CFHR1, and CFHR3 loci was done with multiplex ligation-dependent probe amplification (MLPA) from MRC Holland.
DNA sequencing and analysis
Genomic DNA from 2 aHUS index cases and 2 affected siblings was prepared and subjected to exome capture using NimbleGen 2.1M human exome capture arrays (Life technologies) followed by next generation sequencing on the Illumina sequencing platform as previously described36. Illumina’s processing software ELAND (CASAVA 1.8.2) was used to map reads to the human reference genome (build 19), and SAMtools37 was used to call single nucleotide variants and insertion/deletion at targeted bases. Variants with minor allele frequencies < 1% in the Yale (1,972 European subjects), NHLBI GO (4,300 European and 2,202 African American subjects; last accessed November, 2012), dbSNP (version 135) or 1000 Genomes (1,094 subjects of various ethnicities; May, 2011 data release) databases were selected and annotated for impact on the encoded protein and for conservation of the reference base and amino acid among orthologs across phylogeny. Variants of interest were verified by direct Sanger sequencing. Nomenclature of sequence DGKE variants is based on NCBI Reference Sequence NM_003647.
Genomic DNA of two affected members of kindred 9 underwent targeted enrichment of all exons from the interval 49.4-60.3 Mb on chromosome 17 based on evidence of linkage in this family (see below) followed by next generation sequencing. The targeted capture reagent was prepared by Roche NimbleGen (Madison, WI). Analysis prioritized homozygous protein-altering variants in the linked interval that were not present in the Ensembl SNP database, release 54.
To search for DGKE mutations in additional patients with aHUS the complete coding region of DGKE was amplified by PCR using specific oligonucleotide primer pairs for each of the 11 exons and subjected to direct Sanger sequencing. Samples from 48 index patients with aHUS but without mutation in the known aHUS genes and 38 index patients with heterozygous variants in known aHUS genes were sequenced.
DGKE variants were identified in exome data from subjects not known to have aHUS from Yale and NHLBI exome databases; potential compound heterozygous variants could only be identified in the Yale database since variants in NHLBI are not linked to genotypes.
Origin of DGKE p.Trp322*
Single nucleotide polymorphisms (SNP) flanking DGKE were selected using Haploview 4.238. All had minor allele frequency > 10% and no violation of Hardy-Weinberg equilibrium (p > 0.05) among European subjects (CEU). SNPs were genotyped in 3 apparently unrelated subjects homozygous for p.Trp322* by PCR and direct Sanger sequencing. The boundaries of the homozygous segments were determined for each family by genotyping additional proximal and distal tag SNPs until heterozygous positions were recorded. This defined a shared homozygous haplotype comprising 16 SNPs no more than 400 kilobases in length.
The population frequency of the shared haplotype was estimated from the frequency of the haplotype in each block of linkage disequilibrium (LD); a frequency of 0.1% was assigned to haplotypes not previously detected in HapMap CEU dataset.
Given the rarity of the haplotype bearing DGKE p.Trp322* in CEU subjects, we examined the ancestry of the siblings from kindred 1 who were homozygous for this mutation by principal component analysis (PCA). All tag SNP genotypes were extracted from exome sequence data using a Perl script that produces PLINK-compatible output files. These SNPs were combined with HapMap data (Phase II release, 2010-08-18). Tag SNPs extracted using PLINK’s LD-based SNP pruning algorithm39 were used as inputs to perform PCA with EIGENSTRAT software (version 3.0)40.
DMLE+2.3 software was used to estimate the age of the most common ancestor carrying DGKE p.Trp322*41. DMLE+2.3 uses a Bayesian approach to infer mutation age of a given variant based on observed LD data from polymorphic markers located within the shared segment. Data for six polymorphic loci spanning the shared segment were entered as genotypes. To minimize bias from consanguinity, we treated the data from these 3 kindreds as 3 independent chromosomes; however, results were similar when data from 6 chromosomes were used. Results were also similar whether analyses were performed under a recessive or dominant genetic model. The three unrelated families harboring DGKEW322* are all from western European ancestry. We estimated western Europe’s population growth rate (PGR) using census data from France with the equation PGR = LN(T1/T0)/g, where T1 and T0 are the French populations in the years 2009 (62.5M) and 1806 (21M), respectively, and g is the number of generations during this period (g = 8.1 assuming ~25 years per generation). With these parameters, PGR is estimated at 0.09. When using T0 from the 1954 census data (42.8M; g = 2.2), PGR is 0.17. To calculate the proportion of sampled chromosomes (PSC), we used the following values: number of chromosomes for DGKE nephropathy patients harboring p.Trp322* (n = 8), estimated minor allele frequency in CEU subjects (MAF ≤ 0.01%, based on empiric data from control exomes), and Western Europe’s population (P = 200 millions; PSC = n/[MAF×P×2 chromosomes]). A sensitivity analysis was performed with various combinations of PSC and PGR values centered around the estimates described above. The software ESTIAGE was used to ascertain whether results obtained with DMLE+2.3 were plausible42. ESTIAGE uses a likelihood-based method to estimate mutation age from information related to polymorphic markers located within or at the boundaries of the shared homozygous segment. Twenty-four SNPs spanning the DGKE locus were used for input. Mutation rate was set to 2×10−8. Examples of input files for both softwares are available upon request.
Analysis of linkage and association
Genome-wide analysis of linkage analysis in kindred 9 was performed by comparing inheritance of SNP genotypes (Affymetrix GeneChip® Human Mapping 250K SNP NspI Array) to the inheritance of aHUS in all pedigree members specifying aHUS as a rare recessive trait with high penetrance and rare phenocopies, as described. Multipoint LOD scores were calculated using ALLEGRO43. Haplotypes were reconstructed with ALLEGRO and presented graphically with HaploPainter44.
The significance of rare DGKE variants in all aHUS kindreds was assessed by comparing their segregation to the inheritance of aHUS in the kindreds. Parametric LOD scores were calculated specifying aHUS in kindreds with DGKE mutations as an autosomal recessive trait with complete penetrance and zero phenocopies. Fisher’s exact test was used to compare the prevalence of recessive DGKE variants among non-consanguineous index cases of pediatric-onset aHUS to the corresponding prevalence in 8,475 control exomes from NHLBI and Yale.
DGKE protein expression in platelets and endothelial cells
Venous blood was drawn from 3 healthy adult human volunteers and from 3 wild type C57/Bl6 adult mice at 8- and 18-weeks of age. Platelets were prepared from blood by differential centrifugation as described45. The cytoplasmic and membrane fractions of total protein extracts were harvested by subcellular protein fractionation (Thermo Scientific). Alternatively, whole platelet lysates were prepared by lysis with NP-40 in the presence of protease and phosphatase inhibitors. Whole cell protein extracts from human embryonic venous endothelial cells (HUVECs) were prepared in an analogous fashion.
For western blotting, 50 μg of cytoplasmic and membrane platelet extracts, and total platelet and endothelial protein extracts were subjected to SDS-PAGE and analysed by Western blotting. After blocking, the membranes were probed with the following primary antibodies: for mouse extracts, rabbit polyclonal anti-DGKE (1:1000; Abcam); for human extracts, mouse monoclonal anti-DGKE (1:1000; Sigma) or mouse monoclonal anti-DGKE (1:1000; R&D). Anti-β-tubulin (1:500; Santa Cruz) and anti-Na,K-ATPase (1:500; Cell Signaling) primary antibodies were also used as loading controls and to assess the purity of the cytoplasmic and membrane fractions, respectively. Secondary antibodies linked to horseradish peroxidase (1:5000; Thermo Scientific) and incubated with the membrane for 1 hour at room temperature.
Microscopy
We stained kidney sections from two adult subjects with peri-renal tumors and from aHUS patient 2-7. A standard DAB protocol for renal tissue was used for immunohistochemistry using anti-DGKE and anti-CD34 antibodies. Briefly, paraffin was removed with xylene/alcohol and heat-induced antigen retrieval was done with 10mM citrate, pH 6.0. Endogenous peroxidase activity was quenched with hydrogen peroxide (H2O2). Slides were washed twice with TBS-Tween in between each of the following steps. After blocking with 5% human serum in TBS, the slides were incubated with primary antibodies for 1 h, biotin-labeled secondary antibody for 1 h, and then streptavidin-horseradish peroxidase (Dako, P0397) for 30 minutes. DAB reagent (Dako) was then applied for 5 or 30 minutes and slides were then washed in water. Hematoxylin counterstain was applied to slides before mounting. The primary antibodies included rabbit polyclonal anti-DGKE (1:50; Novus Biologicals) and mouse monoclonal anti-CD34 (1:50; DAKO). The immunoperoxidase protocol used for Supplementary Fig. 7d-e was similar except for the following modifications: the primary antibody was mouse monoclonal anti-DGKE (1:30; R&D) was applied overnight, at 4°C; Tris-EDTA pH9.0 was used for antigen retrieval; PBS-based blocking solution contained 4% bovine serum albumin, 10% normal goat serum, and 0.1% Triton X-100.
Kidney specimens were obtained from adult wild type rats following infusion of PBS and 4% paraformaldehyde and equilibration overnight in 30% sucrose solution prior to freezing. Cryostat sections were prepared from frozen tissue (thickness, 10 microns) and fixed using the Nakane protocol. Slides were washed three times with TBS in between each step described below. After permeabilization with 0.1% Triton X-100, slides were blocked with 10% goat serum, 1% bovine serum albumin and 0.1% Triton X-100 diluted in TBS at room temperature for 1 h. The slides were incubated with primary antibodies at 4°C overnight, and then with secondary antibodies for 1 h at room temperature. Slides were mounted with DAPI nuclear counterstain (Vector Laboratories). Primary antibodies used were directed against DGKE labeled secondary antibodies (Invitrogen) were used at 1:200.
Supplementary Material
ACKNOWLEDGEMENTS
We thank the aHUS subjects, their families, and the health care professionals whose participation made this study possible; Junhui Zhang, Carol Nelson-Williams, SueAnn Mentone, Delphine Beury and others members of the complement laboratory at Hôpital Européen Georges-Pompidou for technical support; the staff of the Yale Center for Genome Analysis for exome production; Shuta Ishibe, Shigeru Shibata, Ute Scholl, Marie-Agnes Dragon-Durey, Lubka Roumenina, Michal Malina, Quentin Vincent and Laurent Abel for helpful discussions; Diane Damotte for the anti-CD34 staining. This work was supported by NIH grants U54 HG006504 01 (Yale Center for Mendelian Genomics), P30 DK079310 05 (Yale O’Brien Center for Kidney Research) and UL1TR00142 07 (Yale Center for Translational Science Award), by grants from the Délégation Régionale à la Recherche Clinique, Assistance Publique – Hôpitaux de Paris to V.F.B., such as Programme Hospitalier de Recherche Clinique, (AOM08198), and the Association pour l’Information et la Recherche dans les maladies Rénales génétiques (AIRG France), and a grant from EURenOmics (2012-305608) to F.S. and V.F.B.. M.L. is the recipient of a Kidney Research Scientist Core Education and National Training (KRESCENT) Program Post-Doctoral Fellowship Award from the Kidney Foundation of Canada and is a member of the Investigative Medicine Ph.D. program at Yale University School of Medicine.
Footnotes
AUTHOR’S CONTRIBUTIONS: M.L., V.F.B., and R.P.L. designed experiments and analyzed data. M.C. and R.P.L. developed the exome analysis protocol. S.M.M., J.D.O., J.A. and H.T. directed the exome capture, DNA sequencing infrastructure and information technology. M.C., M.L., R.P.L., G.N. and P.N. performed bioinformatic and statistical analyses. M.L., W.J. and R.P.L. analyzed age of shared mutation. Sanger sequencing was done by M.L., W.J. and V.F.B.. W.H.T., J.H. and F. F. performed Western blotting experiments. M.L. performed the immunofluorescence studies. M.LQ. and F.F. performed the immunohistochemistry studies on human kidneys. F.Z., S.T., F.N., F.M., D.M., G.D., V.B., B.L., L.C., M.A.M., E.S. and C.L. ascertained and evaluated aHUS patients. C.L., V.F.B and F.S. recruited aHUS patients. N.R.L., G.W M., and M.C.G. provided renal pathology expertise. M.L., V.F.B. and R.P.L. wrote the manuscript.
COMPETING FINANCIAL INTERESTS F.S., F.F., C.L. and V.F.B. have received fees from Alexion Pharmaceuticals for invited lectures and are members of an expert board supported by Alexion Pharmaceuticals. C.L. is an unpaid coordinator for France for the clinical trial “Eculizumab in atypical HUS”. P.N. is a founder, CEO, and shareholder of ATLAS Biolabs GmbH, a service provider for genomic analyses. The authors declare no other competing financial interests.
References
- 1.Neild GH. Haemolytic-uraemic syndrome in practice. Lancet. 1994;343:398–401. doi: 10.1016/s0140-6736(94)91228-9. [DOI] [PubMed] [Google Scholar]
- 2.Loirat C, Frémeaux-Bacchi V. Atypical hemolytic uremic syndrome. Orphanet J. Rare Dis. 2011;6:60. doi: 10.1186/1750-1172-6-60. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 3.Noris M, Mescia F, Remuzzi G. STEC-HUS, atypical HUS and TTP are all diseases of complement activation. Nat. Rev. Nephrol. 2012;8:622–633. doi: 10.1038/nrneph.2012.195. [DOI] [PubMed] [Google Scholar]
- 4.Noris M, Remuzzi G. Atypical hemolytic–uremic syndrome. N. Engl. J. Med. 2009;361:1676–1687. doi: 10.1056/NEJMra0902814. [DOI] [PubMed] [Google Scholar]
- 5.Zuber J, Fakhouri F, Roumenina LT, Loirat C, Frémeaux-Bacchi V. Use of eculizumab for atypical haemolytic uraemic syndrome and C3 glomerulopathies. Nat. Rev. Nephrol. 2012;8:643–657. doi: 10.1038/nrneph.2012.214. [DOI] [PubMed] [Google Scholar]
- 6.Tang W, Bunting M, Zimmerman GA, McIntyre TM, Prescott SM. Molecular cloning of a novel human diacylglycerol kinase highly selective for arachidonate-containing substrates. J. Biol. Chem. 1996;271:10237–10241. [PubMed] [Google Scholar]
- 7.Shulga YV, Topham MK, Epand RM. Regulation and functions of diacylglycerol kinases. Chem. Rev. 2011;111:6186–6208. doi: 10.1021/cr1004106. [DOI] [PubMed] [Google Scholar]
- 8.Rhee SG. Regulation of phosphoinositide-specific phospholipase C. Annu. Rev. Biochem. 2001;70:281–312. doi: 10.1146/annurev.biochem.70.1.281. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 9.Pettitt TR, et al. Diacylglycerol and phosphatidate generated by phospholipases C and D, respectively, have distinct fatty acid compositions and functions. Phospholipase D-derived diacylglycerol does not activate protein kinase C in porcine aortic endothelial cells. J. Biol. Chem. 1997;272:17354–17359. doi: 10.1074/jbc.272.28.17354. [DOI] [PubMed] [Google Scholar]
- 10.Carew MA, Paleolog EM, Pearson JD. The roles of protein kinase C and intracellular Ca2+ in the secretion of von Willebrand factor from human vascular endothelial cells. Biochem. J. 1992;286:631–635. doi: 10.1042/bj2860631. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 11.Ren S, Shatadal S, Shen GX. Protein kinase C-beta mediates lipoprotein-induced generation of PAI-1 from vascular endothelial cells. Am. J. Physiol. Endocr. Metab. 2000;278:E656–662. doi: 10.1152/ajpendo.2000.278.4.E656. [DOI] [PubMed] [Google Scholar]
- 12.Whatley RE, et al. The regulation of platelet-activating factor production in endothelial cells. The role of calcium and protein kinase C. J. Biol. Chem. 1989;264:6325–6333. [PubMed] [Google Scholar]
- 13.Herbert JM, Savi P, Laplace MC, Dumas A, Dol F. Chelerythrine, a selective protein kinase C inhibitor, counteracts pyrogen-induced expression of tissue factor without effect on thrombomodulin down-regulation in endothelial cells. Thromb. Res. 1993;71:487–493. doi: 10.1016/0049-3848(93)90122-5. [DOI] [PubMed] [Google Scholar]
- 14.Levin EG, Marotti KR, Santell L. Protein kinase C and the stimulation of tissue plasminogen activator release from human endothelial cells. Dependence on the elevation of messenger RNA. J. Biol. Chem. 1989;264:16030–16036. [PubMed] [Google Scholar]
- 15.Offermans S. Activation of platelet function through G protein–coupled receptors. Circ. Res. 2006;15:1293–1304. doi: 10.1161/01.RES.0000251742.71301.16. [DOI] [PubMed] [Google Scholar]
- 16.Pettitt TR, Wakelam MJ. Diacylglycerol kinase epsilon, but not zeta, selectively removes polyunsaturated diacylglycerol, inducing altered protein kinase C distribution in vivo. J. Biol. Chem. 1999;274:36181–36186. doi: 10.1074/jbc.274.51.36181. [DOI] [PubMed] [Google Scholar]
- 17.Yada Y, Ozeki T, Kanoh H, Nozawa Y. Purification and characterization of cytosolic diacylglycerol kinases of human platelets. J. Biol. Chem. 1990;265:19237–19243. [PubMed] [Google Scholar]
- 18.Nunn DL, Watson SP. A diacylglycerol kinase inhibitor, R59022, potentiates secretion by and aggregation of thrombin-stimulated human platelets. Biochem. J. 1987;243:809–813. doi: 10.1042/bj2430809. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 19.De Nucci G, Gryglewski RJ, Warner TD, Vane JR. Receptor-mediated release of endothelium-derived relaxing factor and prostacyclin from bovine aortic endothelial cells is coupled. Proc. Natl. Acad. Sci. U S A. 1988;85:2334–2338. doi: 10.1073/pnas.85.7.2334. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 20.Hofmann T, et al. Direct activation of human TRPC6 and TRPC3 channels by diacylglycerol. Nature. 1999;397:259–263. doi: 10.1038/16711. [DOI] [PubMed] [Google Scholar]
- 21.Quack I, et al. PKC mediates beta-arrestin2-dependent nephrin endocytosis in hyperglycemia. J. Biol. Chem. 2011;286:12959–12970. doi: 10.1074/jbc.M110.204024. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 22.Eremina V, et al. VEGF inhibition and renal thrombotic microangiopathy. N. Engl. J. Med. 2008;358:1129–1136. doi: 10.1056/NEJMoa0707330. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 23.Sison K, et al. Glomerular structure and function require paracrine, not autocrine, VEGF-VEGFR-2 signaling. J. Am. Soc. Nephrol. 2010;21:1691–1701. doi: 10.1681/ASN.2010030295. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 24.Hoshi S, Nomoto K, Kuromitsu J, Tomari S, Nagata M. High glucose induced VEGF expression via PKC and ERK in glomerular podocytes. Biochem. Biophys. Res. Commun. 2002;290:177–184. doi: 10.1006/bbrc.2001.6138. [DOI] [PubMed] [Google Scholar]
- 25.Rask-Madsen C, King GL. Differential regulation of VEGF signaling by PKC-alpha and PKC-epsilon in endothelial cells. Arterioscl. Throm. Vasc. 2008;28:919–924. doi: 10.1161/ATVBAHA.108.162842. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 26.Rodriguez de Turco EB, et al. Diacylglycerol kinase epsilon regulates seizure susceptibility and long-term potentiation through arachidonoyl-inositol lipid signaling. Proc. Natl. Acad. Sci. U S A. 2001;98:4740–4745. doi: 10.1073/pnas.081536298. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 27.Lansigan F, Isufi I, Tagoe CE. Microangiopathic haemolytic anaemia resembling thrombotic thrombocytopenic purpura in systemic lupus erythematosus: the role of ADAMTS13. Rheumatology. 2011;50:824–829. doi: 10.1093/rheumatology/keq395. [DOI] [PubMed] [Google Scholar]
- 28.Ganesan C, Maynard SE. Acute kidney injury in pregnancy: the thrombotic microangiopathies. J. Nephrol. 2011;24:554–563. doi: 10.5301/JN.2011.6250. [DOI] [PubMed] [Google Scholar]
- 29.Skerka C, et al. Autoimmune forms of thrombotic micorangiopathy and membranoproliferative glomerulonephritis: Indications for a disease spectrum and common pathogenic principles. Mol. Immunol. 2009;46:2801–2807. doi: 10.1016/j.molimm.2009.05.018. [DOI] [PubMed] [Google Scholar]
- 30.Ozaltin F, et al. DGKE variants cause a glomerular microangiopathy that mimics membranoproliferative GN. J. Am. Soc. Nephrol. 2012 doi: 10.1681/ASN.2012090903. doi: 10.1681/ASN.2012090903. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 31.Nørholm MHH, Shulga YV, Aoki S, Epand RM, von Heijne G. Flanking residues help determine whether a hydrophobic segment adopts a monotopic or bitopic topology in the endoplasmic reticulum membrane. J. Biol. Chem. 2011;286:25284–25290. doi: 10.1074/jbc.M111.244616. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 32.Soldin SJ, Brugnara C, Wong EC, editors. Pediatric Reference Intervals. AACC Press; Washington, DC: 2005. [Google Scholar]
- 33.Sellier-Leclerc A-L, et al. Differential impact of complement mutations on clinical characteristics in atypical hemolytic uremic syndrome. J. Am. Soc. Nephrol. 2007;18:2392–2400. doi: 10.1681/ASN.2006080811. [DOI] [PubMed] [Google Scholar]
- 34.Frémeaux-Bacchi V, et al. Genetics and outcome of atypical hemolytic uremic syndrome: a nationwide French series comparing children and adults. Clin. J. Am. Soc. Nephrol. 2013 doi: 10.2215/CJN.04760512. doi: 10.2215/CJN.04760512. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 35.Roumenina LT, et al. Alternative complement pathway assessment in patients with atypical HUS. J. Immunol. Methods. 2011;365:8–26. doi: 10.1016/j.jim.2010.12.020. [DOI] [PubMed] [Google Scholar]
- 36.Boyden LM, et al. Mutations in kelch-like 3 and cullin 3 cause hypertension and electrolyte abnormalities. Nature. 2012;482:98–102. doi: 10.1038/nature10814. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 37.Li H, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–2079. doi: 10.1093/bioinformatics/btp352. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 38.Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21:263–265. doi: 10.1093/bioinformatics/bth457. [DOI] [PubMed] [Google Scholar]
- 39.Purcell S, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 2007;81:559–575. doi: 10.1086/519795. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 40.Price AL, et al. Principal components analysis corrects for stratification in genome-wide association studies. Nat. Genet. 2006;38:904–909. doi: 10.1038/ng1847. [DOI] [PubMed] [Google Scholar]
- 41.Reeve JP, Rannala B. DMLE+: Bayesian linkage disequilibrium gene mapping. Bioinformatics. 2002;18:894–895. doi: 10.1093/bioinformatics/18.6.894. [DOI] [PubMed] [Google Scholar]
- 42.Genin E, Tullio-Pelet A, Begeot F, Lyonnet S, Abel L. Estimating the age of rare disease mutations: the example of Triple-A syndrome. J. Med. Genet. 2004;41:445–449. doi: 10.1136/jmg.2003.017962. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 43.Gudbjartsson DF, Jonasson K, Frigge ML, Kong A. Allegro, a new computer program for multipoint linkage analysis. Nat. Genet. 2000;25:12–13. doi: 10.1038/75514. [DOI] [PubMed] [Google Scholar]
- 44.Thiele H, Nürnberg P. HaploPainter: a tool for drawing pedigrees with complex haplotypes. Bioinformatics. 2005;21:1730–1732. doi: 10.1093/bioinformatics/bth488. [DOI] [PubMed] [Google Scholar]
- 45.Tang WH, et al. Glucose and collagen regulate human platelet activity through aldose reductase induction of thromboxane. J. Clin. Invest. 2011;121:4462–4476. doi: 10.1172/JCI59291. [DOI] [PMC free article] [PubMed] [Google Scholar]
Associated Data
This section collects any data citations, data availability statements, or supplementary materials included in this article.