Routine use of clinical exome-based next-generation sequencing for evaluation of patients with thrombotic microangiopathies.

Next-generation sequencing is increasingly used for clinical evaluation of patients presenting with thrombotic microangiopathies because it allows for simultaneous interrogation of multiple complement and coagulation pathway genes known to be associated with disease. However, the diagnostic yield is undefined in routine clinical practice. Historic studies relied on case-control cohorts, did not apply current guidelines for variant pathogenicity assessment, and used targeted gene enrichment combined with next-generation sequencing. A clinically enhanced exome, targeting ~54 Mb, was sequenced for 73 patients. Variant analysis and interpretation were performed on genes with biological relevance in thrombotic microangiopathy (C3,CD46, CFB, CFH, CFI, DGKE, and THBD). CFHR3-CFHR1 deletion status was also assessed using multiplex ligation-dependent probe amplification. Variants were classified using American College of Medical Genetics and Genomics guidelines. We identified 5 unique novel and 14 unique rare variants in 25% (18/73) of patients, including a total of 5 pathogenic, 4 likely pathogenic, and 15 variants of uncertain clinical significance. Nine patients had homozygous deletions in CFHR3-CFHR1. The diagnostic yield, defined as the presence of a pathogenic variant, likely pathogenic variant or homozygous deletion of CFHR3-CFHR1, was 25% for all patients tested. Variants of uncertain clinical significance were identified in 21% (15/73) of patients.These results illustrate the expected diagnositic yield in the setting of thrombotic microangiopathies through the application of standardized variant interpretation, and highlight the utility of such an approach. Sequencing a clinically enhanced exome to enable targeted, disease-specific variant analysis is a viable approach. The moderate rate of variants of uncertain clinical significance highlights the paucity of data surrounding the variants in our cohort and illustrates the need for expanded variant curation resources to aid in thrombotic microangiopathy-related disease variant classification.

Thrombotic microangiopathies, defined by microangiopathic hemolytic anemia, microvascular thrombosis, and organ damage, are a diverse group of disorders with varied etiologies (1, 2). Atypical hemolytic uremic syndrome is an event-triggered thrombotic microangiopathy due to over-activation of the alternative complement pathway predominantly involving the kidneys (1, 3). Hereditary and acquired genetic variants in the alternative complement and coagulation pathways have been identified in association with disease (1, 4–9). Due to the ability to evaluate multiple genes simultaneously for known and novel mutations, massively parallel, or ‘next-generation’ sequencing, is poised to play a critical role in identifying patients with thrombotic microangiopathy-related genetic variants that have implications for diagnosis, therapy, and prognosis (3, 7, 10–14). Clinical practice guidelines recommend genetic testing for all patients presenting with atypical hemolytic uremic syndrome to guide prognosis and treatment (15–17).

Historically, when a genetic etiology was considered, single-gene analysis was used to identify variants associated with clinical phenotypes. As the number of causative genes grows, this iterative process becomes costly and time-consuming. The advent of rapid, economical, and efficient next-generation sequencing technologies has revolutionized medical practice. Since completion of the initial draft of the human genome in 2000, the cost and time required to sequence the human genome have decreased dramatically (18). Additional improvements in accuracy and automation allowed next-generation sequencing to become a widespread clinical diagnostic tool for pathogenic variant detection. Several reports described the utility of targeted gene enrichment combined with next-generation sequencing in renal disease evaluation (14, 19–23). Using this approach, investigators reported a positive genetic diagnosis in 43%–61% of patients sequenced for atypical hemolytic uremic syndrome (3, 4, 7, 24). However, none of these studies used application of current variant interpretation standards as defined by the American College of Medical Genetics and Genomics/American Molecular Pathology Association (ACMG/AMP) (25).

The development of the clinical exome holds significant promise for genetic testing in patients with constitutional disease. The clinical research exome has become a dominant exome-capture reagent for clinical next-generation sequencing in the assessment of constitutional disease because it provides for enrichment of greater than 4,600 genes known to be associated with complex and Mendelian disease as curated from databases including Online Mendelian Inheritance in Man, ClinVar, and the Human Gene Mutation Database. This represents a cost-effective technique for generation of hundreds of disease-specific gene panels that is readily adaptable as clinical knowledge advances. With approved human research protocols, variants across the exome may be evaluated for discovery and other clinical research investigations. To our knowledge, no prior studies have utilized clinical exome sequencing followed by targeted, disease specific bioinformatic analysis and reporting in the evaluation of patients with thrombotic microangiopathies.

This study is the first to report the diagnostic yield of next-generation sequencing testing using a clinical exome-based strategy and rigorous application of current clinical variant interpretation standards for patients with thrombotic microangiopathies in routine clinical practice. Our results demonstrate that this approach identifies a genetic cause, defined as identification of pathogenic variants or likely pathogenic variants in 12% of all patients referred for thrombotic microangiopathy genetic testing. The diagnostic yield increases to 25% by including detection of homozygous deletion of CFHR3-CFHR1. This is significant as most cases in this cohort referred for clinical sequencing represent singletons with little clinical data provided.

Materials and Methods

Laboratory Workflow

This study was granted exempt status from the Washington University School of Medicine institutional review board. All tests were requested by licensed physicians (Figure 1). The data presented represent consecutive samples submitted to Genomics and Pathology Services at Washington University in Saint Louis (GPS@WUSTL) between February 1, 2015 and February 29, 2016. The laboratory is College of American Pathologists (CAP)-accredited and Clinical Laboratory Improvement Amendments (CLIA)-certified. Test methodology was two-fold, encompassing both next-generation sequencing for genetic variant detection, and multiplex ligation-dependent probe amplification for detection of deletion of CFHR1-CFHR3. The submitted specimen for testing represented 2-5mls of peripheral blood obtained in a lavender-top EDTA tube with genomic DNA extracted manually using the QIAamp DNA blood mini kit (Qiagen; Valencia, CA). Deletion analysis was performed on extracted DNA as a send-out test to Cincinnati Children’s Hospital. Briefly, oligonucleotide probes hybridize adjacent regions of the CFHR1 and CFHR3 genes. Only probes that are adjacently hybridized are able to be ligated, amplified and quantified using capillary electrophoresis. In the presence of CFHR1-CFHR3 deletions, probe amplification will be diminished or absent (26, 27). Library preparation for the next-generation sequencing assay began by fragmenting DNA to approximately 200 base pairs (bp) by ultrasonication, followed by end repair, A-tailing, and ligation to sequencing adapters. Target capture was performed using the SureSelectXT Clinical Research Exome (Agilent Technologies; Santa Clara, CA) encompassing 54Mb of target space. Libraries were sequenced using an Illumina (San Diego, CA) HiSeq2500 with paired 2x101-bp reads. Analytic sensitivity, specificity, and reproducibility were established per the CAP next-generation sequencing-testing checklist (28).

Figure 1

Flow chart depicting the thrombotic microangiopathy genetic testing results by numbers of patients with variants, variant interpretation and variant distribution among complement and coagulation pathway genes.

Variant Annotation and Reporting

Variant classifications were based on standards and guidelines presented as a joint consensus recommendation published by the American College of Medical Genetics and Genomics and the Association of Molecular Pathology (25). Variant calls were reported using Human Genome Variation Society nomenclature (http://www.hgvs.org/mutnomen), and variant attributes were examined using sequence variation databases including the Exome Aggregation Consortium (ExAC) (exac.broadinstitute.org, encompassing data from the 1000 Genomes project, and NHLBI-GO Exome Sequencing Project, among others, along with annotated data from the Single Nucleotide Polymorphism Database (v135) (http://www.ncbi.nlm.nih.gov/projects/SNP/)), the FH atypical hemolytic uremic syndrome mutation database (http://www.fh-hus.org/), ClinVar (http://www.ncbi.nlm.nih.gov/clinvar/), the HGMD public resource (http://www.hgmd.cf.ac.uk/ac/index.php), and an internally curated clinical-grade database of variants and interpretations housed in the Clinical Genomicist Workspace (PierianDx, St. Louis, MO) (29, 30). Variants were classified on the basis of the sum of evidence surrounding the genomic alteration: level 1-pathogenic, level 2-likely pathogenic, level 3-variant of uncertain significance, level 4-likely benign, level 5-benign (25). All results were reviewed by a pathologist with subspecialty boards in Molecular Genetics from the American Board of Pathology, or a clinical laboratory geneticist certified in Clinical Molecular Genetics and Clinical Cytogenetics from the American Board of Medical Genetics and Genomics prior to release to the patient’s medical record.

Statistics

Statistical analyses were performed using GraphPad software. Categorical variables were examined using the Fisher’s exact test, two-tailed. Continuous variables were evaluated using a Mann-Whitney test.

Results

Assay design and validation

The thrombotic microangiopathy genetic test is performed using the Agilent SureSelectXT Clinical Research exome with the reported gene set encompassing CFH, CFI, CFB, C3, CD46, DGKE, and THBD. In total, the capture reagent targets approximately 54Mb of exonic and selected intronic regions for sequencing. Bioinformatic filtering is used to restrict the analyzed and reported gene set. Validation was performed using 28 genomic DNA samples including 3 HapMap cell lines provided by the National Institute of General Medical Sciences-Coriell cell repositories, 5 blinded peripheral blood samples from individuals without a history of renal disease, 17 blinded patient samples harboring at least one known pathogenic genetic variant, a single unblinded patient sample, and 2 production patient samples. All analyses were based on human reference sequence UCSC build hg19; NCBI build 37.2. Analytic sensitivity, specificity, and positive predictive value to detect single base-pair variation in coding regions throughout the entire Agilent SureSelectXT Clinical Research exome were approximately 96.6%, 100.0%, and 99.3%, respectively, as determined by comparison of genetic sequence called by this assay using HapMap DNA sample NA12878 to the high-confidence genotypes reported at those positions by Complete Genomics (www.completegenomics.com). Analysis of the subset of thrombotic microangiopathy disease-associated genes showed sensitivity, specificity and positive predictive value all approaching 100%. Single nucleotide variants and small insertion/deletion events were called using SAMtools mpileup (31, 32).

Patient Characteristics

Seventy-three patients underwent exome-based thrombotic microangiopathy genetic testing. The demographic characteristics are summarized in Table 1. The average age at the time of testing was 36.7 +/− 16.5 years. Forty-six females and twenty-seven males were tested. Patients were classified based on the referring clinician’s clinical diagnosis which included 61 atypical hemolytic uremic syndrome, 4 thrombotic thrombocytopenic purpura, and 8 thrombotic microangiopathy.

Table 1

Patient demographics

Diagnosis	Number of patients	Age (years)	Gender	Race
aHUS/HUS	61	35.5 +/− 16.2	22M:39F	37 Caucasian 5 Hispanic 6 African American 4 Asian 1 Mediterranean 3 Mixed 5 Unknown
TTP	4	47.8 +/− 3.3	2M:2F	3 Caucasian 1 African American
TMA	8	40.9 +/− 21	3M:5F	5 Caucasian 1 Hispanic 1 African American 1 Asian
All patients	73	36.7 +/− 16.5	27M:46F	45 Caucasian 6 Hispanic 8 African American 5 Asian 1 Mediterranean 3 Mixed 5 Unknown

aHUS – atypical hemolytic uremic syndrome; HUS – hemolytic uremic syndrome; TTP – thrombotic thrombocytopenic purpura; TMA – thrombotic microangiopathy

Variant Interpretation

Prior to variant calling and interpretation, all samples underwent detailed quality control analysis (28). A quality control metrics report is generated detailing the total number of reads, percent mapped to the genome, on-target, on-target unique, mean mapping quality, and depth of coverage. Genomic positions failing to meet a coverage depth of at least 10x, or those failing to meet a variant quality call of 30 were filtered out of analysis. All variants were initially filtered for a population minor allele frequency <5%; this cut-off was chosen since this represents stand-alone evidence for a benign classification according to ACMG/AMP guidelines (25). All variants with a minor allele frequency <5% were analyzed and scored as pathogenic, likely pathogenic, uncertain clinical significance, likely benign, or benign based on the scheme outlined by the ACMG/AMP (25). To facilitate classification, a variant assessment worksheet with the ACMG/AMP criteria is utilized for each variant and archived as part of the permanent testing record.

Genetic Variants

The genetic variants are summarized in Table 2. Of the seventy-three patients sequenced, 27% (20/73) harbored 24 variants classified as variants of uncertain clinical significance, pathogenic, or likely pathogenic and 12% (9/73) had homozygous deletions in CFHR3-CFHR1. Among the 24 single nucleotide variants are five pathogenic variants, four likely pathogenic variants, and fifteen variants of uncertain clinical significance (Figure 1). Six variants (5 of which are unique) are novel, having not been reported in the literature or deposited in existing population or disease-associated variant databases (ExAC, NHLBI, FH-aHUS, dbSNP). Details of the pathogenic/likely pathogenic variants are summarized in Table 3 and variants of uncertain clinical significance in Table 4. Six patients carried multiple genomic alterations: greater than one single nucleotide variant and/or harboring a homozygous CFHR3-CFHR1 deletion (Figure 2). The patients identified with pathogenic or likely pathogenic variants tended to be younger at the time of testing, 29.2 +/− 14.5 years compared with 37.8 +/−16.6 years. However, this difference was not statistically significant.

Table 2

Genetic variants meeting criteria for pathogenic, likely pathogenic or VUS identified in patients referred for TMA genetic testing.

Patient	Gene	Interpretation	Variant (coding)	Variant (genomic)	Type	Clinical Diagnosis
65	C3	LP	p.R1042L	chr19:g.6694471C>A	Missense	aHUS
7	C3	LP	p.R161W	chr19:g.6718128G>A	Missense	aHUS
3	C3	VUS	p.T647M^	chr19:g.6707846G>A	Missense	aHUS
40	C3	VUS	p.Q185H	chr19:g.6714407C>G	Missense	aHUS
80	C3	VUS	p.T1383N	chr19:g.6684423G>T	Missense	aHUS
78	C3	VUS	p.E1636G	chr19:g.6677978T>C	Missense	aHUS
12	CD46	P	c.287-2A>G	chr1:g.207930883A>G	Splice site	HUS
35	CD46	P	c.287-2A>G	chr1:g.207930883A>G	Splice site	Post-partum TMA
25	CD46	P	c.1127+2T>G^	chr1:g.207959029T>G	Splice site	TTP
5	CD46	LP	p.C35Y	chr1:g.207930365G>A	Missense	HUS
38	CD46	VUS	p.P324L	chr1:g.207943690C>T	Missense	aHUS s/p renal transplant
21	CFH	P	c.619+1G>A^	chr1:g.196646798G>A	Splice site	aHUS
84	CFH	VUS	p.I372V^	chr1:g.196658699A>G	Missense	TMA, loss of 2 kidney transplants
31	CFH	VUS	p.Q950H	chr1:g.196709816G>T	Missense	aHUS
12	CFH	VUS	p.G918E^	chr1:g.196706761G>A	Missense	HUS
35	CFH	VUS	p.G918E^	chr1:g.196706761G>A	Missense	Post-partum TMA
69	CFI	LP	p.I370N^	chr4:g.110670413A>T	Missense	aHUS
30	CFI	VUS	p.I306V	chr4:g.110673648T>C	Missense	HUS
45	CFI	VUS	p.D403N	chr4:g.110667600C>T	Missense	aHUS
79	DGKE	P	p.W322*	chr17:g.54926134G>A	Nonsense	aHUS
27	THBD	VUS	p.P495S	chr20:g.23028659G>A	Missense	HUS
79	THBD	VUS	p.A43T	chr20:g.23030015C>T	Missense	aHUS
42	THBD	VUS	p.A43T	chr20:g.23030015C>T	Missense	aHUS
65	THBD	VUS	p.A43T	chr20:g.23030015C>T	Missense	aHUS

TMA – Thrombotic microangiopathy, HUS – hemolytic uremic syndrome, aHUS – atypical hemolytic uremic syndrome, P – pathogenic, LP – likely pathogenic, VUS – variant of uncertain clinical significance

Novel variant

Table 4

Variants of Uncertain Clinical Significance

Disease	Gene	Variant (coding)	SIFT	PolyPhen	PROVEAN	MAF*	GERP	dbSNP
aHUS	CD46	p.P324L	0.094	0.998	−2.10	0.0006365	−0.557	rs41317833
TMA	CFH	p.I372V^	0.025	0.981	−0.72	Not present	−7.48	No
aHUS	CFH	p.Q950H	0.006	0.986	−2.98	0.003583	−4.92	rs149474608
HUS & TMA	CFH	p.G918E^#	0.042	0.992	−5.90	Not present	6.04	No
HUS	CFI	p.I306V	0.458	0.012	0.08	0.0004798	−2.55	rs113273712
aHUS	CFI	p.D403N	1	0.003	1.45	2.48e-05	0.0847	rs139881195
HUS	C3	p.T647M^	0.018	0.321	−1.13	Not present	2.63	No
aHUS	C3	p.Q185H	0.065	0.033	−1.36	Not present	−1.4	No
aHUS	C3	p.T1383N	0.036	0.63	−1.62	9.06e-05	1.59	rs139100972
aHUS	C3	p.E1636G	0.002	0.967	−5.55	Not present	3.95	No
HUS	THBD	p.P495S	0.146	0.273	−2.74	0.000575	−0.884	rs1800578
aHUS x 3 patients	THBD	p.A43T	0.533	0.002	−0.39	0.00343	1.61	rs1800576

TMA – Thrombotic microangiopathy, HUS – hemolytic uremic syndrome, aHUS – atypical hemolytic uremic syndrome

Identified in association with a pathogenic variant in CD46

Refers to the overall minor allele frequency (MAF) in the ExAC database

Novel variant

Figure 2

Summary of thrombotic microangiopathy -genetic testing results. Patients are grouped according to clinical diagnosis. Genes are grouped according to function. Age represents the age at the time of testing. AA – African American, CA – Caucasian, NA – Native American, H – Hispanic, A – Asian, Unk – Unknown, aHUS – atypical hemolytic uremic syndrome, HUS – hemolytic uremic syndrome, TMA – thrombotic microangiopathy, TTP – thrombotic thrombocytopenic purpura, VUS – variant of uncertain clinical significance.

Uniquely occuring pathogenic and likely pathogenic variants (n=8) classified by the ACMG/AMP criteria were compared with unique variants of uncertain clinical significance (n=12) classified using the same criteria. Pathogenic and likely pathogenic variants tended to have lower minor allele frequencies across all populations based on the exome aggregation consortium (ExAC) database, 1.11e-05 ± 2.18e-05 compared with variants of uncertain clinical significance, 7.35e-04 ± 1.32e-03 (p=0.11) (30). The genomic evolutionary rate profiling (GERP++) scores quantifying evolutionary constraint differed significantly between the variant classifications, 4.17 ± 1.86 for pathogenic/likely pathogenic versus -0.16 ± 3.73 for variants of uncertain clinical significance (p=0.004) (33). In silico prediction tools were compared for the missense variants classified as pathogenic/likely pathogenic (n=4) and variants of uncertain clinical significance (n=12). SIFT (release 63), PolyPhen-2 (version 2.2.2, r394), and the protein variation effect analyzer (PROVEAN, version 1.1.5) algorithms showed statistically significant differences between the pathogenic/likely pathogenic variantscompared with variants of uncertain clinical significance, −1.80 ± 2.10 (p=0.001, p=0.03, and p=0.004, respectively) (34).

CFH and CD46 were most commonly identified as harboring single nucleotide variants, each occurring in 7% (5/73) of patients (Figures 1 and 2). The five patients with CFH variants had four distinct single nucleotide variants, three classified as variants of uncertain clinical significance and one as pathogenic (Tables 3 and 4). Two were previously reported (24, 35–39). One patient with the p.Q950H variant of uncertain clinical significance also showed homozygous deletion of CFHR3-CFHR1. Two related patients carrying the p.G918E CFH variant also harbored the c.287-2A>G splice-site pathogenic variant in CD46. CD46 contained four distinct variants (Figures 1 and 2). Of the four variants, two were classified as pathogenic, one likely pathogenic, and one of uncertain clinical significance (Tables 3 and 4).

Table 3

Pathogenic and Likely Pathogenic Variants

Diagnosis	Gene	Variant (coding)	SIFT	PolyPhen	PROVEAN	MAF*	GERP	dbSNP	Human Splicing Finder
HUS & TMA	CD46	c.287-2A>G	n/a	n/a	n/a	3.121e-05	4.07	No	Positive#
TMA	CD46	c.1127+2T>G^	n/a	n/a	n/a	0	2.96	No	Positive#
HUS	CD46	p.C35Y	0	1	−10.16	0	3.72	No	n/a
aHUS	CFH	c.619+1G>A^	n/a	n/a	n/a	0	5.21	No	Positive#
aHUS	CFI	p.I370N^	0	0.998	−6.09	0	5.71	No	n/a
HUS	C3	p.R161W	0.001	0.952	−5.06	0	0.356	No	n/a
aHUS	C3	p.R1042L	0	0.996	−6.16	0	5.76	No	n/a
aHUS	DGKE	p.W332*	n/a	n/a	n/a	5.766e-05	5.59	rs138924661	n/a

TMA – Thrombotic microangiopathy, HUS – hemolytic uremic syndrome, aHUS – atypical hemolytic uremic syndrome

Refers to the overall minor allele frequency (MAF) in the ExAC database

Alteration of wild-type acceptor site, most probably affecting splicing

Novel variant

Three distinct variants were found in CFI, one likely pathogenic and two of uncertain clinical significance (Figure 2, Table 3). The two variants of uncertain clinical significance have been previously reported (40, 41). The patient with the previously reported p.I306V variant also harbored a homozygous deletion in CFHR3-CFHR1.

Six distinct variants were identified in C3, two likely pathogenic and four variants of uncertain clinical significance (Figure 2, Table 4). Five were previously reported in the literature and/or the FH atypical hemolytic uremic syndrome database (4, 24, 42). One patient with a likely pathogenic variant, p.R1042L, also carried the p.A43T variant of uncertain clinical significance in THBD.

THBD harbored variants in four patients (Figure 2). Three patients carried the p.A43T variant and one carried the p.P495S variant, both of uncertain clinical significance. Both variants have been previously reported and have low minor allele frequency (<0.4%) (43). As previously mentioned, one patient with a p.A43T variant also harbored a likely pathogenic variant in C3. Another patient with the p.A43T variant carried a pathogenic variant in DGKE.

The only variant identified in DGKE, p.W322* has been previously described in the literature as pathogenic (Figure 2, Tables 2 and 3) (44).

Diagnostic Yield

Current ACMG/AMP guidelines indicate that variants classified as of uncertain clinical significance are not to be used in clinical decision making (25). Therefore, these are not considered a positive genetic diagnosis in this study for purposes of measuring clinical utility or diagnostic yield. In patients with a clinical diagnosis of thrombotic microangiopathy, atypical hemolytic uremic syndrome, hemolytic uremic syndrome, or thrombotic thrombocytopenic purpura (n=73), 9 had pathogenic or likely pathogenic variants and 9 had homozygous deletions of CFHR3-CFHR1, representing a diagnostic yield of 25% (18/73). Isolated variants of uncertain clinical significance were identified in 10 (14%) patients..

Discussion

We report our one-year experience using exome-based next-generation sequencing in routine diagnostic evaluation of patients referred for thrombotic microangiopathy genetic testing. Our results showed pathogenic or likely pathogenic variants in 12% (9/73) and variants of uncertain clinical significance in 21% (15/73) of all patients tested. This is comparable to the rate of pathogenic and likely pathogenic variants identified in a previous large series of patients sequenced for thrombotic microangiopathy (3). In that particular study, 193 patients were sequenced, and 23 pathogenic or likely pathogenic variants were identified for a diagnostic yield of 12% (3). In our current study, when homozygous deletions of CFHR3-CFHR1 are included, the overall diagnostic yield increases to 21%. Similarly, in the large series by Bu et al., including the rate of homozygous deletions of CFHR3-CFHR1 increases their diagnostic yield to 21% (3). The current study shows that exome-based, targeted next-generation sequencing sequencing in a CAP/CLIA certified laboratory is an efficient and reliable approach for identification of diagnostic variants in patients with thrombotic microangiopathy. The positive diagnostic rate is comparable to a similar series that used non-exome based next-generation sequencing sequencing techniques (3).

Genetic sequencing has the potential to uncover single nucleotide variants of clinical significance that have not been previously described in the literature or reported in existing databases. A significant finding of our study is the identification of 5 such novel variants. Three of these are pathogenic or likely pathogenic based on standardized interpretation guidelines (25). In addition to revealing new pathogenic variants, comprehensive gene sequencing reveals novel variants, and currently undefinable variants, classified as variants of uncertain clinical significance.

In all patients referred for thrombotic microangiopathy genetic testing, variants of uncertain clinical significance represent the largest category of single nucleotide variants identified at 21% (15/73). However, this is significantly lower than the rate of variants of uncertain clinical significance previously reported by Bu et al. of 35% (68/193) (p=0.026) (3). A portion of this discrepancy may reflect the absence of CFHR5 and PLG in our reportable gene-set. However, this discordance may be further attributed to differences in criteria applied to classify variants between the two studies. In the current study, comparison of variant classifications showed significantly higher GERP++ scores and lower minor allele frequency in pathogenic/likely pathogenic variants compared with variants of uncertain clinical significance. Likewise, in silico prediction tools showed significant differences between these two groups.

In the current series, patients referred for thrombotic microangiopathy genetic sequencing most commonly harbor single nucleotide variants in CFH, but variants identified in CD46 are most likely to be classified as pathogenic due to the higher incidence of splice site variants identified in the gene. Homozygous CFHR3-CFHR1 deletions are identified in 12% of patients, indicating an increased risk for the development of anti-complement Factor H antibodies (45).

Resolving variants of uncertain clinical significance is challenging for clinicians due to less than definitive information on the variant and its application to patient management. In such settings, discussions with a clinical team composed of molecular pathologists, geneticists, genetic counselors and nephrologists may be useful. Current recommendations indicate that variants of uncertain clinical significance should not be used in clinical decision making (25). Therefore, to further define associations between gene variants and pathogenesis, genetic findings need to be correlated with functional studies, clinical outcome, pathology findings, familial cosegregation, and treatment response (25). Renal genetic analysis is challenging due to incomplete penetrance, variable expressivity, and complex modes of inheritance. It is anticipated that as more laboratories describe genetic variants, further studies will be performed to illuminate their significance.

In addition to follow-up analyses, our data suggest a need for expanded variant curation resources in the setting of renal disease including a centralized database to deposit annotated disease-associated variants identified as part of routine clinical care. Presently, the ClinVar database serves as a public archive documenting the report of relationships between human phenotype and genetic variation. However the paucity of renal disease associated variation makes this resource limited in function. Expansion of the ClinVar database and other renal genetic knowledgebases would be invaluable to assist laboratories in assigning significance to individual variants with little published information. Such a database needs accurate clinical information in addition to variant information and could serve as a source for other investigators to utilize in attempts to further characterize unique genetic findings.

New disease gene associations continue to emerge. Chong et al. report 915 new disease genes over the last four years (46). It is therefore important to design a test readily adaptable to new discoveries without incurring significant cost and time. Utilizing a unique panel-based next-generation sequencing test requires re-validation for each new gene target added. In contrast, building a test on an exome backbone allows for simpler and more rapid validation of new genes, making this approach nimble and cost-effective. Using this approach, clinically relevant genes are reported in the medical record, while the remaining exome sequence is captured and bioinformatically masked. Thus, these data can be leveraged as new biologic gene targets emerge. Under appropriate consent/human research protections, these data also provide a rich resource of research material to investigate additional genes associated with thrombotic microangiopathy and potential genetic modifiers contributing to states of inappropriate complement activation.

The study is limited by the small sample size, an inherent problem when studying rare diseases such as atypical hemolytic uremic syndrome. Since many of the cases are referred from outside institutions, complete clinical data are not available for majority of patients. Reliance on the referring clinician’s diagnostic impression may artificially lower the diagnostic yield. Nonetheless, this reflects the results clinicians can expect in actual clinical practice. This study describes the use of clinical exome-based next-generation sequencing testing with targeted bioinformatics analysis and reporting for routine care of patients referred for thrombotic microangiopathy genetic testing. Our results highlight the inherent difficulty in applying the current ACMG/AMP guidelines in interpreting single nucleotide variants due to the paucity of available data in the renal disease setting. The results show a lower diagnostic yield than previously reported and a high number of variants of uncertain clinical significance. The decrease in diagnostic yield appears to be due to the absence of variants of uncertain clinical significance for inclusion as a positive diagnostic finding. The current study demonstrates that exome-based next-generation sequencing testing in patients with thrombotic microangiopathies can be applied in a clinical setting and provide meaningful results to assist in clinical care of patients.