Jeroen van Rooij

Quantifying prion disease penetrance using large population control cohorts

Large genomic reference data sets reveal a spectrum of pathogenicity in the prion protein gene and provide genetic validation for a therapeutic strategy in prion disease. Share trumps rare No longer just buzz words, “patient empowerment” and “data sharing” are enabling breakthrough research on rare genetic diseases. Although more than 100,000 genetic variants are believed to drive disease in humans, little is known about penetrance—the probability that a mutation will actually cause disease in the carrier. This conundrum persists because small sample sizes breed imperfect alliance estimates between mutations and disease risk. Now, a patient-turned-scientist joined with a large bioinformatics team to analyze vast amounts of shared data—from the Exome Aggregation Consortium and the 23andMe database—to provide insights into genetic-variant penetrance and possible treatment approaches for a rare, fatal genetic prion disease. More than 100,000 genetic variants are reported to cause Mendelian disease in humans, but the penetrance—the probability that a carrier of the purported disease-causing genotype will indeed develop the disease—is generally unknown. We assess the impact of variants in the prion protein gene (PRNP) on the risk of prion disease by analyzing 16,025 prion disease cases, 60,706 population control exomes, and 531,575 individuals genotyped by 23andMe Inc. We show that missense variants in PRNP previously reported to be pathogenic are at least 30 times more common in the population than expected on the basis of genetic prion disease prevalence. Although some of this excess can be attributed to benign variants falsely assigned as pathogenic, other variants have genuine effects on disease susceptibility but confer lifetime risks ranging from <0.1 to ~100%. We also show that truncating variants in PRNP have position-dependent effects, with true loss-of-function alleles found in healthy older individuals, a finding that supports the safety of therapeutic suppression of prion protein expression.

Disease variants alter transcription factor levels and methylation of their binding sites

Most disease-associated genetic variants are noncoding, making it challenging to design experiments to understand their functional consequences. Identification of expression quantitative trait loci (eQTLs) has been a powerful approach to infer the downstream effects of disease-associated variants, but most of these variants remain unexplained. The analysis of DNA methylation, a key component of the epigenome, offers highly complementary data on the regulatory potential of genomic regions. Here we show that disease-associated variants have widespread effects on DNA methylation in trans that likely reflect differential occupancy of trans binding sites by cis-regulated transcription factors. Using multiple omics data sets from 3,841 Dutch individuals, we identified 1,907 established trait-associated SNPs that affect the methylation levels of 10,141 different CpG sites in trans (false discovery rate (FDR) < 0.05). These included SNPs that affect both the expression of a nearby transcription factor (such as NFKB1, CTCF and NKX2-3) and methylation of its respective binding site across the genome. Trans methylation QTLs effectively expose the downstream effects of disease-associated variants.

Identification of context-dependent expression quantitative trait loci in whole blood

Genetic risk factors often localize to noncoding regions of the genome with unknown effects on disease etiology. Expression quantitative trait loci (eQTLs) help to explain the regulatory mechanisms underlying these genetic associations. Knowledge of the context that determines the nature and strength of eQTLs may help identify cell types relevant to pathophysiology and the regulatory networks underlying disease. Here we generated peripheral blood RNA–seq data from 2,116 unrelated individuals and systematically identified context-dependent eQTLs using a hypothesis-free strategy that does not require previous knowledge of the identity of the modifiers. Of the 23,060 significant cis-regulated genes (false discovery rate (FDR) ≤ 0.05), 2,743 (12%) showed context-dependent eQTL effects. The majority of these effects were influenced by cell type composition. A set of 145 cis-eQTLs depended on type I interferon signaling. Others were modulated by specific transcription factors binding to the eQTL SNPs.

Effect of oral acyclovir after penetrating keratoplasty for herpetic keratitis: a placebo-controlled multicenter trial.

OBJECTIVE To determine the prophylactic effect of oral acyclovir on the recurrence rate of herpetic eye disease after penetrating keratoplasty. DESIGN A randomized, double-masked, placebo-controlled multicenter trial. PARTICIPANTS Sixty-eight consecutive patients (68 eyes) with corneal opacities due to herpetic eye disease who underwent penetrating keratoplasty. INTERVENTION Oral acyclovir 400 mg twice daily or placebo tablets for 6 months. MAIN OUTCOME MEASURES The recurrence rate of herpetic eye disease-related events and rejection episodes, proven by viral cell culture or polymerase chain reaction. RESULTS During the 2-year follow-up period, there were 3 culture-proven herpetic eye disease recurrences in the acyclovir group and 9 in the placebo group. Lifetime survival analysis of the probability of remaining free from recurrence revealed a significantly reduced risk of recurrent herpetic disease in the acyclovir-treated group. CONCLUSION This study suggests that oral acyclovir effectively prevents herpes-related recurrences after penetrating keratoplasty in herpetic eye disease.

Age-related risk factors, culture outcomes, and prognosis in patients admitted with infectious keratitis to two Dutch tertiary referral centers.

Purpose: To assess age-related risk factors (RFs), microbiologic profile, and prognosis of infectious keratitis and create guidelines for prevention and treatment. Methods: Retrospective review of patients with infectious keratitis admitted to 2 Dutch tertiary referral centers from January 2002 to December 2004. Results: Forty-nine patients were admitted to the Academic Medical Center (Amsterdam) and 107 to the Rotterdam Eye Hospital. Mean age was 56.6 ± 24.4 (SD) years; 49.4% were ≥60 years of age. The most common RFs among the elderly were systemic illness (36.4%), ocular surgery (33.8%), topical steroids (26%), blepharitis (20.8%), and herpetic eye disease (28.6%). This was significantly different from the most common RFs among younger patients (contact lens wear, 62.7%; χ2, P = 0.000). Gram-negative infections predominated (52.3%) and were more prevalent among younger patients (χ2, P = 0.000). Gram-positive infections prevailed among the elderly. Untreated patients had higher culture positive rates (68.7%) than patients treated with antibiotics before culturing (41.3%; χ2, P = 0.001). Elderly patients had a higher risk of perforations than younger patients (27.6% vs. 9.9%), a worse prognosis (mean VA, 6/30 vs. 6/10), and more often needed surgery (57.1% vs. 23.4%; P < 0.005 in all cases). Conclusions: Infectious keratitis is a more severe disease in elderly than in younger patients with more complications and a worse prognosis. Elderly patients have multiple and more diverse risk factors, making prevention difficult. Prevention should aim at minimizing topical steroid use and controlling blepharitis, ocular surface disease, and herpetic eye disease. Initial antibiotic treatment should include sufficient coverage of Gram-positive pathogens.

Novel Genetic Variants for Cartilage Thickness and Hip Osteoarthritis.

Osteoarthritis is one of the most frequent and disabling diseases of the elderly. Only few genetic variants have been identified for osteoarthritis, which is partly due to large phenotype heterogeneity. To reduce heterogeneity, we here examined cartilage thickness, one of the structural components of joint health. We conducted a genome-wide association study of minimal joint space width (mJSW), a proxy for cartilage thickness, in a discovery set of 13,013 participants from five different cohorts and replication in 8,227 individuals from seven independent cohorts. We identified five genome-wide significant (GWS, P≤5·0×10−8) SNPs annotated to four distinct loci. In addition, we found two additional loci that were significantly replicated, but results of combined meta-analysis fell just below the genome wide significance threshold. The four novel associated genetic loci were located in/near TGFA (rs2862851), PIK3R1 (rs10471753), SLBP/FGFR3 (rs2236995), and TREH/DDX6 (rs496547), while the other two (DOT1L and SUPT3H/RUNX2) were previously identified. A systematic prioritization for underlying causal genes was performed using diverse lines of evidence. Exome sequencing data (n = 2,050 individuals) indicated that there were no rare exonic variants that could explain the identified associations. In addition, TGFA, FGFR3 and PIK3R1 were differentially expressed in OA cartilage lesions versus non-lesioned cartilage in the same individuals. In conclusion, we identified four novel loci (TGFA, PIK3R1, FGFR3 and TREH) and confirmed two loci known to be associated with cartilage thickness.The identified associations were not caused by rare exonic variants. This is the first report linking TGFA to human OA, which may serve as a new target for future therapies.

Lack of evidence for a role of genetic variation in TMEM230 in the risk for Parkinson's disease in the Caucasian population

Mutations in TMEM230 have recently been associated to Parkinsons disease (PD). To further understand the role of this gene in the Caucasian population, we interrogated our large repository of next generation sequencing data from unrelated PD cases and controls, as well as multiplex families with autosomal dominant PD. We identified 2 heterozygous missense variants in 2 unrelated PD cases and not in our control database (p.Y106H and p.I162V), and a heterozygous missense variant in 2 PD cases from the same family (p.A163T). However, data presented herein is not sufficient to support the role of any of these variants in PD pathology. A series of unified sequence kernel association tests also failed to show a cumulative effect of rare variation in this gene on the risk of PD in the general Caucasian population. Further evaluation of genetic data from different populations is needed to understand the genetic role of TMEM230 in PD etiology.

Uncompromised 10-year survival of oldest old carrying somatic mutations in DNMT3A and TET2.

To the editor: Recent large-scale sequencing studies report recurrent somatic mutations in the blood of elderly individuals in genes previously linked to clonal expansion of hematopoietic stem cells.[1][1][⇓][2][⇓][3]-[4][4] Particularly for DNMT3A and TET2 , a steep age-associated increase in

Characterization of pathogenic SORL1 genetic variants for association with Alzheimer’s disease: a clinical interpretation strategy

Accumulating evidence suggests that genetic variants in the SORL1 gene are associated with Alzheimer disease (AD), but a strategy to identify which variants are pathogenic is lacking. In a discovery sample of 115 SORL1 variants detected in 1908 Dutch AD cases and controls, we identified the variant characteristics associated with SORL1 variant pathogenicity. Findings were replicated in an independent sample of 103 SORL1 variants detected in 3193 AD cases and controls. In a combined sample of the discovery and replication samples, comprising 181 unique SORL1 variants, we developed a strategy to classify SORL1 variants into five subtypes ranging from pathogenic to benign. We tested this pathogenicity screen in SORL1 variants reported in two independent published studies. SORL1 variant pathogenicity is defined by the Combined Annotation Dependent Depletion (CADD) score and the minor allele frequency (MAF) reported by the Exome Aggregation Consortium (ExAC) database. Variants predicted strongly damaging (CADD score >30), which are extremely rare (ExAC-MAF <1 × 10−5) increased AD risk by 12-fold (95% CI 4.2–34.3; P=5 × 10−9). Protein-truncating SORL1 mutations were all unknown to ExAC and occurred exclusively in AD cases. More common SORL1 variants (ExAC-MAF≥1 × 10−5) were not associated with increased AD risk, even when predicted strongly damaging. Findings were independent of gender and the APOE-ɛ4 allele. High-risk SORL1 variants were observed in a substantial proportion of the AD cases analyzed (2%). Based on their effect size, we propose to consider high-risk SORL1 variants next to variants in APOE, PSEN1, PSEN2 and APP for personalized risk assessments in clinical practice.

Hypothesis-free identification of modulators of genetic risk factors

Genetic risk factors often localize in non-coding regions of the genome with unknown effects on disease etiology. Expression quantitative trait loci (eQTLs) help to explain the regulatory mechanisms underlying the association of genetic risk factors with disease. More mechanistic insights can be derived from knowledge of the context, such as cell type or the activity of signaling pathways, influencing the nature and strength of eQTLs. Here, we generated peripheral blood RNA-seq data from 2,116 unrelated Dutch individuals and systematically identified these context-dependent eQTLs using a hypothesis-free strategy that does not require prior knowledge on the identity of the modifiers. Out of the 23,060 significant cis-regulated genes (false discovery rate < 0.05), 2,743 genes (12%) show context-dependent eQTL effects. The majority of those were influenced by cell type composition, revealing eQTLs that are particularly strong in cell types such as CD4+ T-cells, erythrocytes, and even lowly abundant eosinophils. A set of 145 cis-eQTLs were influenced by the activity of the type I interferon signaling pathway and we identified several cis-eQTLs that are modulated by specific transcription factors that bind to the eQTL SNPs. This demonstrates that large-scale eQTL studies in unchallenged individuals can complement perturbation experiments to gain better insight in regulatory networks and their stimuli.