Pieter van der Vlies

Lessons learned from whole exome sequencing in multiplex families affected by a complex genetic disorder, intracranial aneurysm

Genetic risk factors for intracranial aneurysm (IA) are not yet fully understood. Genomewide association studies have been successful at identifying common variants; however, the role of rare variation in IA susceptibility has not been fully explored. In this study, we report the use of whole exome sequencing (WES) in seven densely-affected families (45 individuals) recruited as part of the Familial Intracranial Aneurysm study. WES variants were prioritized by functional prediction, frequency, predicted pathogenicity, and segregation within families. Using these criteria, 68 variants in 68 genes were prioritized across the seven families. Of the genes that were expressed in IA tissue, one gene (TMEM132B) was differentially expressed in aneurysmal samples (n=44) as compared to control samples (n=16) (false discovery rate adjusted p-value=0.023). We demonstrate that sequencing of densely affected families permits exploration of the role of rare variants in a relatively common disease such as IA, although there are important study design considerations for applying sequencing to complex disorders. In this study, we explore methods of WES variant prioritization, including the incorporation of unaffected individuals, multipoint linkage analysis, biological pathway information, and transcriptome profiling. Further studies are needed to validate and characterize the set of variants and genes identified in this study.

Wnt Signaling and Dupuytren's Disease

BACKGROUND Dupuytrens disease is a benign fibromatosis of the hands and fingers that leads to flexion contractures. We hypothesized that multiple genetic and environmental factors influence susceptibility to this disease and sought to identify susceptibility genes to better understand its pathogenesis. METHODS We conducted a genomewide association study of 960 Dutch persons with Dupuytrens disease and 3117 controls (the discovery set) to test for association between the disease and genetic markers. We tested the 35 single-nucleotide polymorphisms (SNPs) most strongly associated with Dupuytrens disease (P<1×10(-4)) in the discovery set in three additional, independent case series comprising a total of 1365 affected persons and 8445 controls from Germany, the United Kingdom, and The Netherlands. RESULTS Initially, we observed a significant genomewide association between Dupuytrens disease and 8 SNPs at three loci. Tests of replication and joint analysis of all data from 2325 patients with Dupuytrens disease and 11,562 controls yielded an association with 11 SNPs from nine different loci (P<5.0×10(-8)). Six of these loci contain genes known to be involved in the Wnt-signaling pathway: WNT4 (rs7524102) (P=2.8×10(-9); odds ratio, 1.28), SFRP4 (rs16879765) (P=5.6×10(-39); odds ratio, 1.98), WNT2 (rs4730775) (P=3.0×10(-8); odds ratio, 0.83), RSPO2 (rs611744) (P=7.9×10(-15); odds ratio, 0.75), SULF1 (rs2912522) (P=2.0×10(-13); odds ratio, 0.72), and WNT7B (rs6519955) (P=3.2×10(-33); odds ratio, 1.54). CONCLUSIONS This study implicates nine different loci involved in genetic susceptibility to Dupuytrens disease. The fact that six of these nine loci harbor genes encoding proteins in the Wnt-signaling pathway suggests that aberrations in this pathway are key to the process of fibromatosis in Dupuytrens disease.

Targeted sequencing by proximity ligation for comprehensive variant detection and local haplotyping.

Despite developments in targeted gene sequencing and whole-genome analysis techniques, the robust detection of all genetic variation, including structural variants, in and around genes of interest and in an allele-specific manner remains a challenge. Here we present targeted locus amplification (TLA), a strategy to selectively amplify and sequence entire genes on the basis of the crosslinking of physically proximal sequences. We show that, unlike other targeted re-sequencing methods, TLA works without detailed prior locus information, as one or a few primer pairs are sufficient for sequencing tens to hundreds of kilobases of surrounding DNA. This enables robust detection of single nucleotide variants, structural variants and gene fusions in clinically relevant genes, including BRCA1 and BRCA2, and enables haplotyping. We show that TLA can also be used to uncover insertion sites and sequences of integrated transgenes and viruses. TLA therefore promises to be a useful method in genetic research and diagnostics when comprehensive or allele-specific genetic information is needed.

Targeted exome sequencing in clear cell renal cell carcinoma tumors suggests aberrant chromatin regulation as a crucial step in ccRCC development

Clear cell renal cell carcinomas are characterized by 3p loss, and by inactivation of Von Hippel Lindau (VHL), a tumorsuppressor gene located at 3p25. Recently, SETD2, located at 3p21, was identified as a new candidate ccRCC tumor‐suppressor gene. The combined mutational frequency in ccRCC tumors of VHL and SETD2 suggests that there are still undiscovered tumor‐suppressor genes on 3p. We screened all genes on 3p for mutations in 10 primary ccRCC tumors using exome‐sequencing. We identified inactivating mutations in VHL, PBRM1, and BAP1. Sequencing of PBRM1 in ccRCC‐derived cell lines confirmed its frequent inactivation in ccRCC. PBRM1 encodes for BAF180, the chromatin targeting subunit of the SWI/SNF complex. BAP1 encodes for BRCA1 associated protein‐1, involved in histone deubiquitination. Taken together, the accumulating data suggest an important role for aberrant chromatin regulation in ccRCC development. Hum Mutat 33:1059–1062, 2012.

Towards Sustained Silencing of HER2/neu in Cancer By Epigenetic Editing

The human epidermal growth factor receptor-2 (HER2/neu/ERBB2) is overexpressed in several cancer types. Although therapies targeting the HER2/neu protein result in inhibition of cell proliferation, the anticancer effect might be further optimized by limiting HER2/neu expression at the DNA level. Towards this aim, epigenetic editing was performed to suppress HER2/neu expression by inducing epigenetic silencing marks on the HER2/neu promoter.HER2/neu expression and HER2/neu promoter epigenetic modification status were determined in a panel of ovarian and breast cancer cell lines. HER2/neu-overexpressing cancer cells were transduced to express a zinc finger protein (ZFP), targeting the HER2/neugene, fused to histone methyltransferases (G9a, SUV39-H1)/super KRAB domain (SKD). Epigenetic assessment of the HER2/neu promoter showed that HER2/neu-ZFP fused to G9a efficiently induced the intended silencing histone methylation mark (H3K9me2). Importantly, H3K9me2 induction was associated with a dramatic downregulation of HER2/neu expression in HER2/neu- overexpressing cells. Downregulation by SKD, traditionally considered transient in nature, was associated with removal of the histone acetylation mark (H3ac). The downregulation of HER2/neu by induced H3K9 methylation and/or reduced H3 acetylation was sufficient to effectively inhibit cellular metabolic activity and clonogenicity. Furthermore, genome-wide analysis indicated preferential binding of the ZFP to its target sequence. These results not only show that H3K9 methylation can be induced but also that this epigenetic mark was instructive in promoting downregulation of HER2/neu expression. Implications: Epigenetic editing provides a novel (synergistic) approach to modulate expression of oncogenes. Mol Cancer Res; 11(9); 1029–39. ©2013 AACR.

FISH and array-CGH analysis of a complex chromosome 3 aberration suggests that loss of CNTN4 and CRBN contributes to mental retardation in 3pter deletions.

Imbalances of 3p telomeric sequences cause 3p− and trisomy 3p syndrome, respectively, showing distinct, but also shared clinical features. No causative genes have been identified in trisomy 3p patients, but for the 3p− syndrome, there is growing evidence that monosomy for one or more of four genes at 3pter, CHL1, CNTN4, CRBN, and MEGAP/srGAP3, may play a causative role. We describe here an analysis of a complex chromosome 3p aberration in a severely mentally retarded patient that revealed two adjacent segments with different copy number gains and a distal deletion. The deletion in this patient included the loci for CHL1, CNTN4, and CRBN, and narrowed the critical segment associated with the 3p− syndrome to 1.5 Mb, including the loci for CNTN4 and CRBN. We speculate that the deletion contributes more to this patients phenotype than the gains that were observed. We suggest that 3p− syndrome associated features are primarily caused by loss of CNTN4 and CRBN, with loss of CHL1 probably having an additional detrimental effect on the cognitive functioning of the present patient.

Genomic aberrations in squamous cell lung carcinoma related to lymph node or distant metastasis

About 50% of patients presenting with resectable lung cancer develop distant metastases within 5 years. Genomic markers predicting metastatic behaviour of squamous cell lung carcinoma (SCC) are currently underexposed. We analyzed a cohort of patients with primary SCC using array-based comparative genomic hybridization (aCGH) to identify which genomic aberrations are related to metastatic behaviour. The cohort consisted of 34 patients with a follow-up of at least 5 years, 8 with metastases in regional lymph nodes only and 26 patients without any metastases at the time of surgery. Eleven of the latter 26 developed metastases in distant organs within 3 years after surgery. Copy number changes observed in at least 40% of all SCC included gains at chromosomal arms 3q, 5p, 8q, 19q, 20p, 22q and losses at 3p, 4p, 4q, 5q, 8p and 9p. High copy number amplifications were observed at 2p15-p16, 3q24-q29, 8p11-p12, 8q23-q24, and 12p12, containing candidate oncogenes such as BCL11A, REL, ECT2, PIK3CA, ADAM9, MYC and KRAS. Amplification of 2p15-p16 is a novel finding in SCC. Another novel finding is the homozygous deletion observed at 4q33-34.1 in 15% of the SCC cases. Gains at 7q36, 8p12, 10q22, 12p12, loss at 4p14 and the homozygous deletions at 4q occurred significantly more frequent in SCC from patients with lymph node metastases only. SCC from patients with distant metastases showed a significantly higher gain frequency at 8q22-q24 and loss at 8p23 and 13q21, and a significantly lower gain frequency at 2p12 and 2p16 and loss at 11q25 compared with SCC from patients without metastases. Of these, gains at 7q, 8p and 10q were restricted to SCC with lymph node metastasis and gain at 8q was restricted to patients with distant metastasis. Two genomic aberrations, i.e. loss of 4p and gain of 19q12 were observed more frequently in SCC with only lymph node metastases as compared to SCC with distant metastases. In conclusion, we identified genomic aberrations in primary SCC that were related to lymph node or distant metastases.

High expression of calcium‐binding proteins, S100A10, S100A11 and CALM2 in anaplastic large cell lymphoma

Anaplastic large cell lymphomas (ALCL) are characterised by the presence of CD30‐positive large cells, which usually are of T‐cell type. Based on the presence or absence of translocations involving the anaplastic lymphoma kinase (ALK) locus, ALCL cases can be divided into two groups. To gain more insight in the biology of ALCL, we applied serial analysis of gene expression (SAGE) on the Karpas299 cell line and identified 25 up‐ and 19 downregulated genes. Comparison of the differentially expressed genes with DNA copy number changes in Karpas299 revealed that two overexpressed genes, S100A10 and S100A11, were located in an amplicon suggesting that the increased mRNA levels were caused by DNA amplification. Quantitative reverse transcription polymerase chain reaction on 5 ALCL cell lines and 12 ALCL tissues confirmed the SAGE data for 13 out of 14 up‐ and one out of four downregulated genes. Immunohistochemical staining confirmed the presence of S100A10, a calcium‐binding protein, in three out of five ALK+ and all 7 ALK− ALCL cases. S100A11 staining was confirmed in all ALK+ and six of seven ALK− ALCL cases. Three of the upregulated genes represented calcium‐binding proteins, which suggest that altered intracellular signaling might be associated with the oncogenesis of ALCL.

Hypoxia and Complement-and-Coagulation Pathways in the Deceased Organ Donor as the Major Target for Intervention to Improve Renal Allograft Outcome.

Background In the last few decades, strategies to improve allograft survival after kidney transplantation have been directed to recipient-dependent mechanisms of renal injury. In contrast, no such efforts have been made to optimize organ quality in the donor. Optimizing deceased donor kidney quality opens new possibilities to improve renal allograft outcome. Methods A total of 554 kidney biopsies were taken from donation after brain death (DBD) and donation after cardiac death (DCD) kidneys before donation, after cold ischemia and after reperfusion. Healthy living donor kidney biopsies served as controls. Transcriptomics was performed by whole genome microarray analyses followed by functional pathway analyses. Results Before organ retrieval and before cessation of blood circulation, metabolic pathways related to hypoxia and complement-and-coagulation cascades were the major pathways enhanced in DBD donors. Similar pathways were also enriched in DCD donors after the first warm ischemia time. Shortly after reperfusion of DCD grafts, pathways related to prolonged and worsening deprivation of oxygen were associated with delayed graft function in the recipient. Conclusion In conclusion, this large deceased donor study shows enrichment of hypoxia and complement-and-coagulation pathways already in DBD donors before cessation of blood flow, before organ retrieval. Therefore, future intervention therapies should target hypoxia and complement-and-coagulation cascades in the donor to improve renal allograft outcome in the recipient.

Systematic identification of tRNAome and its dynamics in Lactococcus lactis

Transfer RNAs (tRNAs) through their abundance and modification pattern significantly influence protein translation. Here, we present a systematic analysis of the tRNAome of Lactococcus lactis. Using the next‐generation sequencing approach, we identified 40 tRNAs which carry 16 different post‐transcriptional modifications as revealed by mass spectrometry analysis. While small modifications are located in the tRNA body, hypermodified nucleotides are mainly present in the anticodon loop, which through wobbling expand the decoding potential of the tRNAs. Using tRNA‐based microarrays, we also determined the dynamics in tRNA abundance upon changes in the growth rate and heterologous protein overexpression stress. With a fourfold increase in the growth rate, the relative abundance of tRNAs cognate to low abundance codons decrease, while the tRNAs cognate to major codons remain mostly unchanged. Significant changes in the tRNA abundances are observed upon protein overexpression stress, which does not correlate with the codon usage of the overexpressed gene but rather reflects the altered expression of housekeeping genes.