A Transcriptionally Distinct CXCL13+CD103+CD8+ T-cell Population Is Associated with B-cell Recruitment and Neoantigen Load in Human Cancer.

The chemokine CXCL13 mediates recruitment of B cells to tumors and is essential for the formation of tertiary lymphoid structures (TLSs). TLSs are thought to support antitumor immunity and are associated with improved prognosis. However, it remains unknown whether TLSs are formed in response to the general inflammatory character of the tumor microenvironment, or rather, are induced by (neo)antigen-specific adaptive immunity. We here report on the finding that the TGFß-dependent CD103+CD8+ tumor-infiltrating T-cell (TIL) subpopulation expressed and produced CXCL13. Accordingly, CD8+ T cells from peripheral blood activated in the presence of TGFß upregulated CD103 and secreted CXCL13. Conversely, inhibition of TGFß receptor signaling abrogated CXCL13 production. CXCL13+CD103+CD8+ TILs correlated with B-cell recruitment, TLSs, and neoantigen burden in six cohorts of human tumors. Altogether, our findings indicated that TGFß plays a noncanonical role in coordinating immune responses against human tumors and suggest a potential role for CXCL13+CD103+CD8+ TILs in mediating B-cell recruitment and TLS formation in human tumors.

The emerging landscape of dynamic DNA methylation in early childhood.

BACKGROUND: DNA methylation has been found to associate with disease, aging and environmental exposure, but it is unknown how genome, environment and disease influence DNA methylation dynamics in childhood. RESULTS: By analysing 538 paired DNA blood samples from children at birth and at 4-5 years old and 726 paired samples from children at 4 and 8 years old from four European birth cohorts using the Illumina Infinium Human Methylation 450 k chip, we have identified 14,150 consistent age-differential methylation sites (a-DMSs) at epigenome-wide significance of p < 1.14 × 10-7. Genes with an increase in age-differential methylation were enriched in pathways related to 'development', and were more often located in bivalent transcription start site (TSS) regions, which can silence or activate expression of developmental genes. Genes with a decrease in age-differential methylation were involved in cell signalling, and enriched on H3K27ac, which can predict developmental state. Maternal smoking tended to decrease methylation levels at the identified da-DMSs. We also found 101 a-DMSs (0.71%) that were regulated by genetic variants using cis-differential Methylation Quantitative Trait Locus (cis-dMeQTL) mapping. Moreover, a-DMS-associated genes during early development were significantly more likely to be linked with disease. CONCLUSION: Our study provides new insights into the dynamic epigenetic landscape of the first 8 years of life.

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.

Towards sustained silencing of HER2/neu in cancer by epigenetic editing.

UNLABELLED: 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.

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.

Wnt signaling and Dupuytren's disease.

BACKGROUND: Dupuytren's 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 Dupuytren's 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 Dupuytren's 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 Dupuytren's disease and 8 SNPs at three loci. Tests of replication and joint analysis of all data from 2325 patients with Dupuytren's 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 Dupuytren's 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 Dupuytren's disease.

A chronic obstructive pulmonary disease related signature in squamous cell lung cancer.

The epidemiological relationship between squamous cell lung cancer (SCC) and chronic obstructive pulmonary disease (COPD), both smoking-related diseases, suggests that they have also a genetic basis. We compared 35 SCC patients with and without COPD with whole-genome gene expression profiles of laser microdissected tissue. Validation of differential expression was performed for 25 genes using quantitative (q)RT-PCR. Subsequently, we performed array-based CGH on the same tumor samples. We found that 374 probes were differentially expressed in SCC from patients with and without COPD. Forty-four probes were derived from genes with mitochondrial functions and 34 probes were located on 5q. All these probes showed a lower expression level in SCC from non-COPD patients. For a random selection of 25 mitochondrial and 5q genes, we confirmed the differential expression by qRT-PCR. Loss of 3p, 5q and 9p was observed, via array-CGH, to be more frequent in SCC from non-COPD patients as compared to SCC from COPD patients. Combination of chromosomal aberrations and the location of the differentially expressed genes showed significant association for loss of the 5q31.2-31.3 region and reduced expression of the 5q genes. In conclusion, a more frequent loss of 5q and a low expression of genes located on 5q in SCC tumors of non-COPD patients compared to tumors from COPD patients was identified suggesting that different oncogenetic pathways are operational in patients with and without COPD.

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.

Global correlation of genome and transcriptome changes in classical Hodgkin lymphoma.

To identify genes involved in the pathogenesis of classical Hodgkin lymphoma (cHL), we performed serial analysis of gene expression (SAGE) and array-based comparative genomic hybridization (aCGH). Comparison of SAGE libraries of cHL cell lines L428 and L1236 with that of germinal centre B cells revealed consistent overexpression of only 14 genes. In contrast, 141 genes were downregulated in both cHL cell lines, including many B cell and HLA genes. aCGH revealed gain of 2p, 7p, 9p, 11q and Xq and loss of 4q and 11q. Eighteen percent of the differentially expressed genes mapped to regions with loss or gain and a good correlation was observed between underexpression and loss or overexpression and gain of DNA. Remarkably, gain of 2p and 9p did not correlate with increased expression of the proposed target genes c-REL and JAK2. Downregulation of many genes within the HLA region also did not correlate with loss of DNA. FSCN1 and IRAK1 mapping at genomic loci (7p and Xq) that frequently showed gain were overexpressed in cHL cell lines and might be involved in the pathogenesis of cHL.

Functional analysis of lung tumor suppressor activity at 3p21.3.

The early and frequent occurrence of deletions at 3p21.3 in lung cancer has led to the consideration of this chromosomal region as a lung cancer (LUCA) critical region with tumor suppressor activity. We covered this 19 genes-containing region with overlapping P1 artificial chromosomes (PACs), in which genes are likely accompanied by their own promoters or other regulatory sequences. With these PACs we transfected cells from a small cell lung cancer (SCLC) cell line which readily caused tumors in nude mice. Per PAC we selected two cell clones with a low number of PAC copies integrated at a single genomic site. The selected clones were s.c. injected into nude mice to investigate whether the integrated genes suppressed the tumor-inducing capacity of the original SCLC cell line. We could demonstrate PAC-specific gene expression in the transfected cells. All of the PAC integration sites were different. It appeared that introduction of a PAC or even an empty PAC vector causes some chromosomal instability, which in principle may either promote or inhibit cell growth. However, both cell clones with integration of the same PAC from the centromeric part of the LUCA region in different genomic sites were the sole pair of clones that caused smaller tumors than did the original SCLC cell line. This suggests that rather than the induced chromosomal instability, the DNA sequence of that PAC, which in addition to two protein-encoding genes contains at least one potential miRNA gene, is responsible for the tumor suppressor activity.

Array comparative genomic hybridization reveals a very high frequency of deletions of the long arm of chromosome 6 in testicular lymphoma.

Despite the fact that numerous studies have been performed on diffuse large B-cell lymphoma (DLBCL), only few have concerned extranodal lymphomas occurring in the testis. We performed a cytogenetic and molecular study of 17 testicular non-Hodgkin lymphomas, of which 14 were proven primary DLBCL of the testis. Cytogenetic analysis revealed in 8 out of 11 evaluable cases a structural abnormality of the long arm of chromosome 6, with deletion or addition of material of unknown origin, and with breakpoints spanning the region 6q12-6q23. The cytogenetic findings were confirmed by fluorescent in situ hybridization (FISH) with a chromosome 6 painting probe. Using array based-comparative genomic hybridization on 16 evaluable cases, including 5 cases not tested by cytogenetics or FISH, 14 (88%) showed chromosome 6q deletions. We identified two regions of minimal deletion (RMD), at 104-113 Mb (6q16.3-q21) and 137.5-138.8 Mb (6q23.3), respectively. In one case, we observed a 2.7 Mb homozygous deletion ranging from 135.3 to 138.0 Mb that partly overlapped with the RMD at 6q23.3. Our study indicates that 6q deletions play a major pathogenetic role in DLBCL of the testis and that many of these deletions are part of unbalanced translocations.

A DGGE system for comprehensive mutation screening of BRCA1 and BRCA2: application in a Dutch cancer clinic setting.

Rapid and reliable identification of deleterious changes in the breast cancer genes BRCA1 and BRCA2 has become one of the major issues in most DNA services laboratories. To rapidly detect all possible changes within the coding and splice site determining sequences of the breast cancer genes, we established a semiautomated denaturing gradient gel electrophoresis (DGGE) mutation scanning system. All exons of both genes are covered by the DGGE scan, comprising 120 amplicons. We use a semiautomated approach, amplifying all individual amplicons with the same PCR program, after which the amplicons are pooled. DGGE is performed using three slightly different gel conditions. Validation was performed using DNA samples with known sequence variants in 107 of the 120 amplicons; all variants were detected. This DGGE mutation scanning, in combination with a PCR test for two Dutch founder deletions in BRCA1 was then applied in 431 families in which 52 deleterious changes and 70 unclassified variants were found. Fifteen unclassified variants were not reported before. The system was easily adopted by five other laboratories, where in another 3,593 families both exons 11 were analyzed by the protein truncation test (PTT) and the remaining exons by DGGE. In total, a deleterious change (nonsense, frameshift, splice-site mutation, or large deletion) was found in 661 families (16.4%), 462 in BRCA1 (11.5%), 197 in BRCA2 (4.9%), and in two index cases a deleterious change in both BRCA1 and BRCA2 was identified. Eleven deleterious changes in BRCA1 and 36 in BRCA2 had not been reported before. In conclusion, this DGGE mutation screening method for BRCA1 and BRCA2 is proven to be highly sensitive and is easy to adopt, which makes screening of large numbers of patients feasible. The results of screening of BRCA1 and BRCA2 in more than 4,000 families present a valuable overview of mutations in the Dutch population.

BCL6 alternative breakpoint region break and homozygous deletion of 17q24 in the nodular lymphocyte predominance type of Hodgkin's lymphoma-derived cell line DEV.

DEV is the only cell line derived from nodular lymphocyte predominance type of Hodgkin's lymphoma (NLPHL); however, a comprehensive report about the genetic and immunophenotypic profile of this unique cell line is lacking. We analyzed DEV with respect to immunophenotype and genetic aberrations. The immunostaining revealed positivity for CD45, CD20, CD22, CD79a, IgA2, CD80, CD86, CD74, and BCL6. Cytogenetically, DEV has complex chromosome 3 translocations involving chromosomes 7, 14, and 22. A detailed analysis of the 3q27 breakpoint of the der(3)t(3;14)(p14;q32)t(3;22)(q27;q11.2) revealed a break in the BCL6 alternative breakpoint region. Using array comparative genomic hybridization, a 3-megabase homozygous deletion at 17q24.1-24.2 was identified. Fluorescence in situ hybridization indicated the presence of 2 chromosome 17 homologues, each of which carried a small interstitial deletion. Eight microsatellite markers flanking the homozygously deleted region all showed a homozygous pattern suggesting loss of one of the parental alleles. D17S1809 and D17S1816 could not be amplified using DEV DNA, in keeping with a location within the homozygously deleted segment. In conclusion, DEV has an immunophenotype that is consistent with the neoplastic cells of NLPHL cases, the lymphocytic and histiocytic cells. We demonstrated involvement of the BCL6 gene based on the presence of a breakpoint in the alternative breakpoint region and nuclear staining for BCL6 protein and identified a homozygously deleted region at 17q24.

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.

A substantial proportion of microsatellite-unstable colon tumors carry TP53 mutations while not showing chromosomal instability.

Chromosomal instability in colon tumors implies the presence of numerical and structural chromosome aberrations and is further characterized by the absence of microsatellite instability and the occurrence of KRAS and/or TP53 mutations. In a previous screening of 194 colon tumors for both microsatellite instability and TP53 mutation, we found 25 microsatellite-unstable tumors, in 9 (36%) of which, presumed to be chromosomally stable, there were TP53 mutations. This prompted us to investigate whether a TP53 mutation in these microsatellite-unstable tumors would be an indicator of chromosomal instability, that is, whether this would be a category of tumors showing both microsatellite and chromosomal instability. For chromosomal instability assessment, we performed array-comparative genomic hybridization analysis of tumor and control DNA extracted from formalin-fixed, paraffin-embedded stage III colon tumor specimens. The array consisted of 435 subtelomere-specific BACs. We compared all but one (whose DNA was of bad quality) of the microsatellite-unstable TP53 mutation-containing tumors (8) with a similarly sized group of microsatellite-unstable tumors without TP53 mutation (11). Microsatellite-unstable tumors with a TP53 mutation showed on average 0.9 aberrations (range 0-3) when assessed with this array system. Those without a TP53 mutation showed on average 0.7 aberrations (range 0-2). Thus, microsatellite-unstable tumors showed few chromosomal abnormalities regardless of TP53 mutation status. Because, in our study, the microsatellite-stable tumors had on average 3.4chromosomal abnormalities (range 0-7), a clear difference exists between microsatellite-unstable and -stable tumors. Because a substantial proportion of microsatellite-unstable colon tumors carry a TP53 mutation while showing relativelyfewchromosomal aberrations, a TP53 mutation in these tumors cannot be considered to be an indicator of chromosomal instability.