Maarten P. G. Massink

Loss-of-Function Mutations in the WNT Co-receptor LRP6 Cause Autosomal-Dominant Oligodontia

Tooth agenesis is one of the most common developmental anomalies in man. Oligodontia, a severe form of tooth agenesis, occurs both as an isolated anomaly and as a syndromal feature. We performed exome sequencing on 20 unrelated individuals with apparent non-syndromic oligodontia and failed to detect mutations in genes previously associated with oligodontia. In three of the probands, we detected heterozygous variants in LRP6, and sequencing of additional oligodontia-affected individuals yielded one additional mutation in LRP6. Three mutations (c.1144_1145dupAG [p.Ala383Glyfs(∗)8], c.1779dupT [p.Glu594(∗)], and c.2224_2225dupTT [p.Leu742Phefs(∗)7]) are predicted to truncate the protein, whereas the fourth (c.56C>T [p.Ala19Val]) is a missense variant of a conserved residue located at the cleavage site of the proteins signal peptide. All four affected individuals harboring a LRP6 mutation had a family history of tooth agenesis. LRP6 encodes a transmembrane cell-surface protein that functions as a co-receptor with members from the Frizzled protein family in the canonical Wnt/β-catenin signaling cascade. In this same pathway, WNT10A was recently identified as a major contributor in the etiology of non-syndromic oligodontia. We show that the LRP6 missense variant (c.56C>T) results in altered glycosylation and improper subcellular localization of the protein, resulting in abrogated activation of the Wnt pathway. Our results identify LRP6 variants as contributing to the etiology of non-syndromic autosomal-dominant oligodontia and suggest that this gene is a candidate for screening in DNA diagnostics.

Proteomics of Mouse BRCA1-deficient Mammary Tumors Identifies DNA Repair Proteins with Potential Diagnostic and Prognostic Value in Human Breast Cancer

Breast cancer 1, early onset (BRCA1) hereditary breast cancer, a type of cancer with defects in the homology-directed DNA repair pathway, would benefit from the identification of proteins for diagnosis, which might also be of potential use as screening, prognostic, or predictive markers. Sporadic breast cancers with defects in the BRCA1 pathway might also be diagnosed. We employed proteomics based on one-dimensional gel electrophoresis in combination with nano-LC-MS/MS and spectral counting to compare the protein profiles of mammary tumor tissues of genetic mouse models either deficient or proficient in BRCA1. We identified a total of 3,545 proteins, of which 801 were significantly differentially regulated between the BRCA1-deficient and -proficient breast tumors. Pathway and protein complex analysis identified DNA repair and related functions as the major processes associated with the up-regulated proteins in the BRCA1-deficient tumors. In addition, by selecting highly connected nodes, we identified a BRCA1 deficiency signature of 45 proteins that enriches for homology-directed DNA repair deficiency in human gene expression breast cancer data sets. This signature also exhibits prognostic power across multiple data sets, with optimal performance in a data set enriched in tumors deficient in homology-directed DNA repair. In conclusion, by comparing mouse proteomes from BRCA1-proficient and -deficient mammary tumors, we were able to identify several markers associated with BRCA1 deficiency and a prognostic signature for human breast cancer deficient in homology-directed DNA repair.

High resolution SNP array profiling identifies variability in retinoblastoma genome stability

Both hereditary and nonhereditary retinoblastoma (Rb) are commonly initiated by loss of both copies of the retinoblastoma tumor suppressor gene (RB1), while additional genomic changes are required for tumor initiation and progression. Our aim was to determine whether there is genomic heterogeneity between different clinical Rb subtypes. Therefore, 21 Rb tumors from 11 hereditary patients and 10 nonhereditary Rb patients were analyzed using high‐resolution single nucleotide polymorphism (SNP) arrays and gene losses and gains were validated with Multiplex Ligation‐dependent Probe Amplification. In these tumors only a few focal aberrations were detected. The most frequent was a focal gain on chromosome 2p24.3, the minimal region of gain encompassing the oncogene MYCN. The genes BAZ1A, OTX2, FUT8, and AKT1 were detected in four focal regions on chromosome 14 in one nonhereditary Rb. There was a large difference in number of copy number aberrations between tumors. A subset of nonhereditary Rbs turned out to be the most genomic unstable, while especially very young patients with hereditary Rb display stable genomes. Established Rb copy number aberrations, including gain of chromosome arm 1q and loss of chromosome arm 16q, turned out to be preferentially associated with the nonhereditary Rbs with later age of diagnosis. In contrast, copy number neutral loss of heterozygosity was detected mainly on chromosome 13, where RB1 resides, irrespective of hereditary status or age. Focal amplifications and deletions and copy number neutral loss of heterozygosity besides chromosome 13 appear to be rare events in retinoblastoma.

Proper genomic profiling of (BRCA1‐mutated) basal‐like breast carcinomas requires prior removal of tumor infiltrating lymphocytes

BRCA1‐mutated breast carcinomas may have distinct biological features, suggesting the involvement of specific oncogenic pathways in tumor development. The identification of genomic aberrations characteristic for BRCA1‐mutated breast carcinomas could lead to a better understanding of BRCA1‐associated oncogenic events and could prove valuable in clinical testing for BRCA1‐involvement in patients.

Genomic profiling of CHEK2*1100delC-mutated breast carcinomas

BackgroundCHEK2*1100delC is a moderate-risk breast cancer susceptibility allele with a high prevalence in the Netherlands. We performed copy number and gene expression profiling to investigate whether CHEK2*1100delC breast cancers harbor characteristic genomic aberrations, as seen for BRCA1 mutated breast cancers.MethodsWe performed high-resolution SNP array and gene expression profiling of 120 familial breast carcinomas selected from a larger cohort of 155 familial breast tumors, including BRCA1, BRCA2, and CHEK2 mutant tumors. Gene expression analyses based on a mRNA immune signature was used to identify samples with relative low amounts of tumor infiltrating lymphocytes (TILs), which were previously found to disturb tumor copy number and LOH (loss of heterozygosity) profiling. We specifically compared the genomic and gene expression profiles of CHEK2*1100delC breast cancers (n = 14) with BRCAX (familial non-BRCA1/BRCA2/CHEK2*1100delC mutated) breast cancers (n = 34) of the luminal intrinsic subtypes for which both SNP-array and gene expression data is available.ResultsHigh amounts of TILs were found in a relatively small number of luminal breast cancers as compared to breast cancers of the basal-like subtype. As expected, these samples mostly have very few copy number aberrations and no detectable regions of LOH. By unsupervised hierarchical clustering of copy number data we observed a great degree of heterogeneity amongst the CHEK2*1100delC breast cancers, comparable to the BRCAX breast cancers. Furthermore, copy number aberrations were mostly seen at low frequencies in both the CHEK2*1100delC and BRCAX group of breast cancers. However, supervised class comparison identified copy number loss of chromosomal arm 1p to be associated with CHEK2*1100delC status.ConclusionsIn conclusion, in contrast to basal-like BRCA1 mutated breast cancers, no apparent specific somatic copy number aberration (CNA) profile for CHEK2*1100delC breast cancers was found. With the possible exception of copy number loss of chromosomal arm 1p in a subset of tumors, which might be involved in CHEK2 tumorigenesis. This difference in CNAs profiles might be explained by the need for BRCA1-deficient tumor cells to acquire survival factors, by for example specific copy number aberrations, to expand. Such factors may not be needed for breast tumors with a defect in a non-essential gene such as CHEK2.

Somatic genomic alterations in retinoblastoma beyond RB1 are rare and limited to copy number changes

Retinoblastoma is a rare childhood cancer initiated by RB1 mutation or MYCN amplification, while additional alterations may be required for tumor development. However, the view on single nucleotide variants is very limited. To better understand oncogenesis, we determined the genomic landscape of retinoblastoma. We performed exome sequencing of 71 retinoblastomas and matched blood DNA. Next, we determined the presence of single nucleotide variants, copy number alterations and viruses. Aside from RB1, recurrent gene mutations were very rare. Only a limited fraction of tumors showed BCOR (7/71, 10%) or CREBBP alterations (3/71, 4%). No evidence was found for the presence of viruses. Instead, specific somatic copy number alterations were more common, particularly in patients diagnosed at later age. Recurrent alterations of chromosomal arms often involved less than one copy, also in highly pure tumor samples, suggesting within-tumor heterogeneity. Our results show that retinoblastoma is among the least mutated cancers and signify the extreme sensitivity of the childhood retina for RB1 loss. We hypothesize that retinoblastomas arising later in retinal development benefit more from subclonal secondary alterations and therefore, these alterations are more selected for in these tumors. Targeted therapy based on these subclonal events might be insufficient for complete tumor control.

A Meta-Analysis of Retinoblastoma Copy Numbers Refines the List of Possible Driver Genes Involved in Tumor Progression.

Background While RB1 loss initiates retinoblastoma development, additional somatic copy number alterations (SCNAs) can drive tumor progression. Although SCNAs have been identified with good concordance between studies at a cytoband resolution, accurate identification of single genes for all recurrent SCNAs is still challenging. This study presents a comprehensive meta-analysis of genome-wide SCNAs integrated with gene expression profiling data, narrowing down the list of plausible retinoblastoma driver genes. Methods We performed SCNA profiling of 45 primary retinoblastoma samples and eight retinoblastoma cell lines by high-resolution microarrays. We combined our data with genomic, clinical and histopathological data of ten published genome-wide SCNA studies, which strongly enhanced the power of our analyses (N = 310). Results Comprehensive recurrence analysis of SCNAs in all studies integrated with gene expression data allowed us to reduce candidate gene lists for 1q, 2p, 6p, 7q and 13q to a limited gene set. Besides the well-established driver genes RB1 (13q-loss) and MYCN (2p-gain) we identified CRB1 and NEK7 (1q-gain), SOX4 (6p-gain) and NUP205 (7q-gain) as novel retinoblastoma driver candidates. Depending on the sample subset and algorithms used, alternative candidates were identified including MIR181 (1q-gain) and DEK (6p gain). Remarkably, our study showed that copy number gains rarely exceeded change of one copy, even in pure tumor samples with 100% homozygosity at the RB1 locus (N = 34), which is indicative for intra-tumor heterogeneity. In addition, profound between-tumor variability was observed that was associated with age at diagnosis and differentiation grades. Interpretation Since focal alterations at commonly altered chromosome regions were rare except for 2p24.3 (MYCN), further functional validation of the oncogenic potential of the described candidate genes is now required. For further investigations, our study provides a refined and revised set of candidate retinoblastoma driver genes.

Coregulation of FANCA and BRCA1 in human cells

Fanconi anemia (FA) is a genetically heterogeneous syndrome associated with increased cancer predisposition. The underlying genes govern the FA pathway which functions to protect the genome during the S-phase of the cell cycle. While upregulation of FA genes has been linked to chemotherapy resistance, little is known about their regulation in response to proliferative stimuli. The purpose of this study was to examine how FA genes are regulated, especially in relation to the cell cycle, in order to reveal their possible participation in biochemical networks. Expression of 14 FA genes was monitored in two human cell-cycle models and in two RB1/E2F pathway-associated primary cancers, retinoblastoma and basal breast cancer. In silico studies were performed to further evaluate coregulation and identify connected networks and diseases. Only FANCA was consistently induced over 2-fold; FANCF failed to exhibit any regulatory fluctuations. Two tools exploiting public data sets indicated coregulation of FANCA with BRCA1. Upregulation of FANCA and BRCA1 correlated with upregulation of E2F3. Genes coregulated with both FANCA and BRCA1 were enriched for MeSH-Term id(s) genomic instability, microcephaly, and Bloom syndrome, and enriched for the cellular component centrosome. The regulation of FA genes appears highly divergent. In RB1-linked tumors, upregulation of FA network genes was associated with reduced expression of FANCF. FANCA and BRCA1 may jointly act in a subnetwork - supporting vital function(s) at the subcellular level (centrosome) as well as at the level of embryonic development (mechanisms controlling head circumference).

Variants in members of the cadherin-catenin complex, CDH1 and CTNND1, cause blepharocheilodontic syndrome

Blepharocheilodontic syndrome (BCDS) consists of lagophthalmia, ectropion of the lower eyelids, distichiasis, euryblepharon, cleft lip/palate and dental anomalies and has autosomal dominant inheritance with variable expression. We identified heterozygous variants in two genes of the cadherin–catenin complex, CDH1, encoding E-cadherin, and CTNND1, encoding p120 catenin delta1 in 15 of 17 BCDS index patients, as was recently described in a different publication. CDH1 plays an essential role in epithelial cell adherence; CTNND1 binds to CDH1 and controls the stability of the complex. Functional experiments in zebrafish and human cells showed that the CDH1 variants impair the cell adhesion function of the cadherin–catenin complex in a dominant-negative manner. Variants in CDH1 have been linked to familial hereditary diffuse gastric cancer and invasive lobular breast cancer; however, no cases of gastric or breast cancer have been reported in our BCDS cases. Functional experiments reported here indicated the BCDS variants comprise a distinct class of CDH1 variants. Altogether, we identified the genetic cause of BCDS enabling DNA diagnostics and counseling, in addition we describe a novel class of dominant negative CDH1 variants.

Intestinal Failure and Aberrant Lipid Metabolism in Patients With DGAT1 Deficiency

Background & Aims Congenital diarrheal disorders are rare inherited intestinal disorders characterized by intractable, sometimes life-threatening, diarrhea and nutrient malabsorption; some have been associated with mutations in diacylglycerol-acyltransferase 1 (DGAT1), which catalyzes formation of triacylglycerol from diacylglycerol and acyl-CoA. We investigated the mechanisms by which DGAT1 deficiency contributes to intestinal failure using patient-derived organoids. Methods We collected blood samples from 10 patients, from 6 unrelated pedigrees, who presented with early-onset severe diarrhea and/or vomiting, hypoalbuminemia, and/or (fatal) protein-losing enteropathy with intestinal failure; we performed next-generation sequencing analysis of DNA from 8 patients. Organoids were generated from duodenal biopsies from 3 patients and 3 healthy individuals (controls). Caco-2 cells and patient-derived dermal fibroblasts were transfected or transduced with vectors that express full-length or mutant forms of DGAT1 or full-length DGAT2. We performed CRISPR/Cas9-guided disruption of DGAT1 in control intestinal organoids. Cells and organoids were analyzed by immunoblot, immunofluorescence, flow cytometry, chromatography, quantitative real-time polymerase chain reaction, and for the activity of caspases 3 and 7. Results In the 10 patients, we identified 5 bi-allelic loss-of-function mutations in DGAT1. In patient-derived fibroblasts and organoids, the mutations reduced expression of DGAT1 protein and altered triacylglycerol metabolism, resulting in decreased lipid droplet formation after oleic acid addition. Expression of full-length DGAT2 in patient-derived fibroblasts restored formation of lipid droplets. Organoids derived from patients with DGAT1 mutations were more susceptible to lipid-induced cell death than control organoids. Conclusions We identified a large cohort of patients with congenital diarrheal disorders with mutations in DGAT1 that reduced expression of its product; dermal fibroblasts and intestinal organoids derived from these patients had altered lipid metabolism and were susceptible to lipid-induced cell death. Expression of full-length wildtype DGAT1 or DGAT2 restored normal lipid metabolism in these cells. These findings indicate the importance of DGAT1 in fat metabolism and lipotoxicity in the intestinal epithelium. A fat-free diet might serve as the first line of therapy for patients with reduced DGAT1 expression. It is important to identify genetic variants associated with congenital diarrheal disorders for proper diagnosis and selection of treatment strategies.