A new variant in the ZCCHC8 gene: diverse clinical phenotypes and expression in the lung.

Introduction: Pulmonary fibrosis is a severe disease which can be familial. A genetic cause can only be found in ∼40% of families. Searching for shared novel genetic variants may aid the discovery of new genetic causes of disease. Methods: Whole-exome sequencing was performed in 152 unrelated patients with a suspected genetic cause of pulmonary fibrosis from the St Antonius interstitial lung disease biobank. Variants of interest were selected by filtering for novel, potentially deleterious variants that were present in at least three unrelated pulmonary fibrosis patients. Results: The novel c.586G>A p.(E196K) variant in the ZCCHC8 gene was observed in three unrelated patients: two familial patients and one sporadic patient, who was later genealogically linked to one of the families. The variant was identified in nine additional relatives with pulmonary fibrosis and other telomere-related phenotypes, such as pulmonary arterial venous malformations, emphysema, myelodysplastic syndrome, acute myeloid leukaemia and dyskeratosis congenita. One family showed incomplete segregation, with absence of the variant in one pulmonary fibrosis patient who carried a PARN variant. The majority of ZCCHC8 variant carriers showed short telomeres in blood. ZCCHC8 protein was located in different lung cell types, including alveolar type 2 (AT2) pneumocytes, the culprit cells in pulmonary fibrosis. AT2 cells showed telomere shortening and increased DNA damage, which was comparable to patients with sporadic pulmonary fibrosis and those with pulmonary fibrosis carrying a telomere-related gene variant, respectively. Discussion: The ZCCHC8 c.586G>A variant confirms the involvement of ZCCHC8 in pulmonary fibrosis and short-telomere syndromes and underlines the importance of including the ZCCHC8 gene in diagnostic gene panels for these diseases.

Clonality, Antigen Recognition, and Suppression of CD8+ T Cells Differentially Affect Prognosis of Breast Cancer Subtypes.

PURPOSE: In breast cancer, response rates to immune therapies are generally low and differ significantly across molecular subtypes, urging a better understanding of immunogenicity and immune evasion. EXPERIMENTAL DESIGN: We interrogated large gene-expression data sets including 867 node-negative, treatment-naïve breast cancer patients (microarray data) and 347 breast cancer patients (whole-genome sequencing and transcriptome data) according to parameters of T cells as well as immune microenvironment in relation to patient survival. RESULTS: We developed a 109-immune gene signature that captures abundance of CD8 tumor-infiltrating lymphocytes (TIL) and is prognostic in basal-like, her2, and luminal B breast cancer, but not in luminal A or normal-like breast cancer. Basal-like and her2 are characterized by highest CD8 TIL abundance, highest T-cell clonality, highest frequencies of memory T cells, and highest antigenicity, yet only the former shows highest expression level of immune and metabolic checkpoints and highest frequency of myeloid suppressor cells. Also, luminal B shows a high antigenicity and T-cell clonality, yet a low abundance of CD8 TILs. In contrast, luminal A and normal-like both show a low antigenicity, and notably, a low and high abundance of CD8 TILs, respectively, which associates with T-cell influx parameters, such as expression of adhesion molecules. CONCLUSIONS: Collectively, our data argue that not only CD8 T-cell presence itself, but rather T-cell clonality, T-cell subset distribution, coinhibition, and antigen presentation reflect occurrence of a CD8 T-cell response in breast cancer subtypes, which have been aborted by distinct T-cell-suppressive mechanisms, providing a rationale for subtype-specific combination immune therapies.

Natural helix 9 mutants of PPARγ differently affect its transcriptional activity.

OBJECTIVE: The nuclear receptor PPARγ is the master regulator of adipocyte differentiation, distribution, and function. In addition, PPARγ induces terminal differentiation of several epithelial cell lineages, including colon epithelia. Loss-of-function mutations in PPARG result in familial partial lipodystrophy subtype 3 (FPDL3), a rare condition characterized by aberrant adipose tissue distribution and severe metabolic complications, including diabetes. Mutations in PPARG have also been reported in sporadic colorectal cancers, but the significance of these mutations is unclear. Studying these natural PPARG mutations provides valuable insights into structure-function relationships in the PPARγ protein. We functionally characterized a novel FPLD3-associated PPARγ L451P mutation in helix 9 of the ligand binding domain (LBD). Interestingly, substitution of the adjacent amino acid K450 was previously reported in a human colon carcinoma cell line. METHODS: We performed a detailed side-by-side functional comparison of these two PPARγ mutants. RESULTS: PPARγ L451P shows multiple intermolecular defects, including impaired cofactor binding and reduced RXRα heterodimerisation and subsequent DNA binding, but not in DBD-LBD interdomain communication. The K450Q mutant displays none of these functional defects. Other colon cancer-associated PPARγ mutants displayed diverse phenotypes, ranging from complete loss of activity to wildtype activity. CONCLUSIONS: Amino acid changes in helix 9 can differently affect LBD integrity and function. In addition, FPLD3-associated PPARγ mutations consistently cause intra- and/or intermolecular defects; colon cancer-associated PPARγ mutations on the other hand may play a role in colon cancer onset and progression, but this is not due to their effects on the most well-studied functional characteristics of PPARγ.

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.

Genomic profiling of CHEK2*1100delC-mutated breast carcinomas.

BACKGROUND: CHEK2*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. METHODS: We 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. RESULTS: High 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. CONCLUSIONS: In 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.

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

INTRODUCTION: 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. METHODS: For this purpose, genomic and gene expression profiles of basal-like BRCA1-mutated breast tumors (n = 27) were compared with basal-like familial BRCAX (non-BRCA1/2/CHEK2*1100delC) tumors (n = 14) in a familial cohort of 120 breast carcinomas. RESULTS: Genome wide copy number profiles of the BRCA1-mutated breast carcinomas in our data appeared heterogeneous. Gene expression analyses identified varying amounts of tumor infiltrating lymphocytes (TILs) as a major cause for this heterogeneity. Indeed, selecting tumors with relative low amounts of TILs, resulted in the identification of three known but also five previously unrecognized BRCA1-associated copy number aberrations. Moreover, these aberrations occurred with high frequencies in the BRCA1-mutated tumor samples. Using these regions it was possible to discriminate BRCA1-mutated from BRCAX breast carcinomas, and they were validated in two independent cohorts. To further substantiate our findings, we used flow cytometry to isolate cancer cells from formalin-fixed, paraffin-embedded, BRCA1-mutated triple negative breast carcinomas with estimated TIL percentages of 40% and higher. Genomic profiles of sorted and unsorted fractions were compared by shallow whole genome sequencing and confirm our findings. CONCLUSION: This study shows that genomic profiling of in particular basal-like, and thus BRCA1-mutated, breast carcinomas is severely affected by the presence of high numbers of TILs. Previous reports on genomic profiling of BRCA1-mutated breast carcinomas have largely neglected this. Therefore, our findings have direct consequences on the interpretation of published genomic data. Also, these findings could prove valuable in light of currently used genomic tools for assessing BRCA1-involvement in breast cancer patients and pathogenicity assessment of BRCA1 variants of unknown significance. The BRCA1-associated genomic aberrations identified in this study provide possible leads to a better understanding of BRCA1-associated oncogenesis.

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.

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.