Search

Using genome-wide transcriptome analysis by RNA sequencing of first trimester plasma RNA, we tested whether the identification of pregnancies at risk of developing pre-eclampsia with or without preterm birth or growth restriction is possible between weeks 9-14, prior to the appearance of clinical symptoms. We implemented a metaheuristic approach in the self-learning SVM algorithm for differential gene expression analysis of normal pregnancies (n = 108), affected pregnancies (n = 34) and non-pregnant controls (n = 19). Presymptomatic candidate markers for affected pregnancies were validated by RT-qPCR in first trimester samples (n = 34) from an independent cohort. PRKG1 was significantly downregulated in a subset of pregnancies with birth weights below the 10thpercentile as shared symptom. The NRIP1/ZEB2 ratio was found to be upregulated in pregnancies with pre-eclampsia or trisomy 21. Complementary quantitative analysis of both the linear and circular forms of NRIP1 permitted discrimination between pre-eclampsia and trisomy 21. Pre-eclamptic pregnancies showed an increase in linear NRIP1 compared to circular NRIP1, while trisomy 21 pregnancies did not.

Copy number variations (CNVs) are major causative contributors both in the genesis of genetic diseases and human neoplasias. While "High-Throughput" sequencing technologies are increasingly becoming the primary choice for genomic screening analysis, their ability to efficiently detect CNVs is still heterogeneous and remains to be developed. The aim of this white paper is to provide a guiding framework for the future contributions of ELIXIR's recently established h uman CNV Community, with implications beyond human disease diagnostics and population genomics. This white paper is the direct result of a strategy meeting that took place in September 2018 in Hinxton (UK) and involved representatives of 11 ELIXIR Nodes. The meeting led to the definition of priority objectives and tasks, to address a wide range of CNV-related challenges ranging from detection and interpretation to sharing and training. Here, we provide suggestions on how to align these tasks within the ELIXIR Platforms strategy, and on how to frame the activities of this new ELIXIR Community in the international context., Competing Interests: No competing interests were disclosed., (Copyright: © 2020 Salgado D et al.)

Objectives: The majority of patients with locally advanced larynx or hypopharynx squamous cell carcinoma are treated with organ-preserving chemoradiotherapy (CRT). Clinical outcome following CRT varies greatly. We hypothesized that tumor microRNA (miRNA) expression is predictive for outcome following CRT., Methods: Next-generation sequencing (NGS) miRNA profiling was performed on 37 formalin-fixed paraffin-embedded (FFPE) tumor samples. Patients with a recurrence-free survival (RFS) of less than 2 years and patients with late/no recurrence within 2 years were compared by differential expression analysis. Tumor-specific miRNAs were selected based on normal mucosa miRNA expression data from The Cancer Genome Atlas database. A model was constructed to predict outcome using group-regularized penalized logistic ridge regression. Candidate miRNAs were validated by RT-qPCR in the initial sample set as well as in 46 additional samples., Results: Thirteen miRNAs were differentially expressed (p < 0.05, FDR < 0.1) according to outcome group. Initial class prediction in the NGS cohort (n = 37) resulted in a model combining five miRNAs and disease stage, able to predict CRT outcome with an area under the curve (AUC) of 0.82. In the RT-qPCR cohort (n = 83), 25 patients (30%) experienced early recurrence (median RFS 8 months; median follow-up 42 months). Class prediction resulted in a model combining let-7i-5p, miR-192-5p and disease stage, able to discriminate patients with good versus poor clinical outcome (AUC:0.80)., Conclusion: The combined miRNA expression and disease stage prediction model for CRT outcome is superior to using either factor alone. This study indicates NGS miRNA profiling using FFPE specimens is feasible, resulting in clinically relevant biomarkers., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

Background: Patients with longstanding inflammatory bowel disease (IBD; ie, ulcerative colitis and Crohn's disease) have an increased risk of colorectal cancer (CRC). Due to ongoing inflammation, IBD-associated dysplastic lesions can develop. These lesions have an increased risk to progress to cancer compared with sporadic adenomas, which are also found in these patients. Differentiating between these 2 types of dysplasia remains challenging, both clinically and histologically, while treatment strategies may differ. Therefore, the aim of this study was to investigate molecular alterations associated with colorectal dysplasia to cancer progression in IBD and evaluate to what extent these alterations differ from sporadic adenomas., Methods: DNA copy number aberrations and mutation analyses of 48 genes were performed by next-generation sequencing in 43 IBD-associated dysplastic lesions, 30 of which were dysplastic and 13 of which were cancers. Results were compared with existing DNA copy number and mutation data from 118 sporadic adenomas and 24 sporadic cancers., Results: Inflammatory bowel disease-associated dysplastic lesions harbor patterns of DNA copy number aberrations comparable to carcinomas, which are rare in sporadic adenomas. TP53 mutation was the most frequent mutation observed in IBD-associated dysplastic lesions and in cancers. FBXW7 was mutated significantly more often in IBD-associated dysplastic lesions than in sporadic adenomas., Conclusions: Inflammatory bowel disease-associated dysplastic lesions show more DNA copy number aberrations than sporadic adenomas. TP53 and FBXW7 mutations appear to be involved in the development of IBD-associated dysplastic lesions and cancer. These findings indicate that IBD-associated dysplastic lesions are more genomically unstable, possibly reflecting a faster progression toward cancer., (© 2019 Crohn’s & Colitis Foundation. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Multiple tumors in patients are frequently diagnosed, either synchronous or metachronous. The distinction between a second primary and a metastasis is important for treatment. Chromosomal DNA copy number aberrations (CNA) patterns are highly unique to specific tumors. The aim of this study was to assess genome-wide CNA-patterns as method to identify clonally related tumors in a prospective cohort of patients with synchronous or metachronous tumors, with at least one intrapulmonary tumor. In total, 139 tumor pairs from 90 patients were examined: 35 synchronous and 104 metachronous pairs. Results of CNA were compared to histological type, clinicopathological methods (Martini-Melamed-classification (MM) and ACCP-2013-criteria), and, if available, EGFR- and KRAS-mutation analysis. CNA-results were clonal in 74 pairs (53%), non-clonal in 33 pairs (24%), and inconclusive in 32 pairs (23%). Histological similarity was found in 130 pairs (94%). Concordance between histology and conclusive CNA-results was 69% (74 of 107 pairs: 72 clonal and two non-clonal). In 31 of 103 pairs with similar histology, genetics revealed non-clonality. In two out of four pairs with non-matching histology, genetics revealed clonality. The subgroups of synchronous and metachronous pairs showed similar outcome for the comparison of histological versus CNA-results. MM-classification and ACCP-2013-criteria, applicable on 34 pairs, and CNA-results were concordant in 50% and 62% respectively. Concordance between mutation matching and conclusive CNA-results was 89% (8 of 9 pairs: six clonal and two non-clonal). Interestingly, in one patient both tumors had the same KRAS mutation, but the CNA result was non-clonal. In conclusion, although some concordance between histological comparison and CNA profiling is present, arguments exist to prefer extensive molecular testing to determine whether a second tumor is a metastasis or a second primary., Competing Interests: The authors have declared that no competing interests exist.

Aims: Outcomes of colorectal cancer (CRC) treatment and survival have steadily improved during the past decades, accompanied by an increased risk of developing second primary tumours and metastatic tumours at unusual sites. Metastatic CRC can show mucosal colonisation, thereby mimicking a second primary tumour. This potential confusion could lead to incorrect diagnosis and consequently inadequate treatment of the patient. The aim of this study was to differentiate between metastatic CRC and a second primary (gallbladder cancer, GBC) using a combination of standard histopathology and molecular techniques., Methods and Results: Ten consecutive patients with both CRC and GBC were identified in our region using the Dutch National Pathology Archive (PALGA). Two patients served as negative controls. Histology of GBC was reviewed by nine pathologists. A combination of immunohistochemistry, microsatellite analysis, genomewide DNA copy number analysis and targeted somatic mutation analysis was used to aid in differential diagnosis. In two patients, CRC and GBC were clonally related, as confirmed by somatic mutation analysis. For one case, this was confirmed by genomewide DNA copy number analysis. However, in both cases, pathologists initially considered the GBC as a second primary tumour., Conclusions: Metastatic CRC displaying mucosal colonisation is often misinterpreted as a second primary tumour. A combination of traditional histopathology and molecular techniques improves this interpretation, and lowers the risk of inadequate treatment., (© 2019 The Authors. Histopathology Published by John Wiley & Sons Ltd.)

Summary: Chromosomal copy number aberrations can be efficiently detected and quantified using low-coverage whole-genome sequencing, but analysis is hampered by the lack of knowledge on absolute DNA copy numbers and tumor purity. Here, we describe an analytical tool for Absolute Copy number Estimation, ACE, which scales relative copy number signals from chromosomal segments to optimally fit absolute copy numbers, without the need for additional genetic information, such as SNP data. In doing so, ACE derives an estimate of tumor purity as well. ACE facilitates analysis of large numbers of samples, while maintaining the flexibility to customize models and generate output of single samples., Availability and Implementation: ACE is freely available via www.bioconductor.org and at www.github.com/tgac-vumc/ACE., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2018. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

An important factor in quality control of non-invasive prenatal screening (NIPS) or testing (NIPT) is a sufficient percentage of fetal DNA to avoid false-negative results. Here we evaluate 14,379 shallow whole-genome sequenced diagnostic NIPS samples, as well as negative controls, for both technical and biological factors that can influence fetal fraction and its assessment. Technically, bioinformatics analyses can have a profound impact on fetal fraction determination. We found best performance for fetal fraction determination with the Y chromosome based tool DEFRAG for male fetuses and the count based tool SeqFF for female fetuses. Biologically, gestational age of up to 21 weeks and maternal age had no influence on fetal fraction, while an increase in weight and BMI had a negative influence on fetal fraction. While a trend was observed, no statistically significant difference in fetal fraction was found between trisomy and normal samples. Overall, these results confirm the influence of biological factors and give insight into technical factors that can affect fetal fractions in NIPS.

Offering self-sampling for HPV testing improves the effectiveness of current cervical screening programs by increasing population coverage. Molecular markers directly applicable on self-samples are needed to stratify HPV-positive women at risk of cervical cancer (so-called triage) and to avoid over-referral and overtreatment. Deregulated microRNAs (miRNAs) have been implicated in the development of cervical cancer, and represent potential triage markers. However, it is unknown whether deregulated miRNA expression is reflected in self-samples. Our study is the first to establish genome-wide miRNA profiles in HPV-positive self-samples to identify miRNAs that can predict the presence of CIN3 and cervical cancer in self-samples. Small RNA sequencing (sRNA-Seq) was conducted to determine genome-wide miRNA expression profiles in 74 HPV-positive self-samples of women with and without cervical precancer (CIN3). The optimal miRNA marker panel for CIN3 detection was determined by GRridge, a penalized method on logistic regression. Six miRNAs were validated by qPCR in 191 independent HPV-positive self-samples. Classification of sRNA-Seq data yielded a 9-miRNA marker panel with a combined area under the curve (AUC) of 0.89 for CIN3 detection. Validation by qPCR resulted in a combined AUC of 0.78 for CIN3+ detection. Our study shows that deregulated miRNA expression associated with CIN3 and cervical cancer development can be detected by sRNA-Seq in HPV-positive self-samples. Validation by qPCR indicates that miRNA expression analysis offers a promising novel molecular triage strategy for CIN3 and cervical cancer detection applicable to self-samples., (© 2018 The Authors. International Journal of Cancer published by John Wiley & Sons Ltd on behalf of UICC.)

Vulvar squamous cell carcinoma (VSCC) and precancerous vulvar intraepithelial neoplasia (VIN) can develop through human papillomavirus (HPV)-dependent and -independent pathways, indicating a heterogeneous disease. Only a minority of VIN progress, but current clinicopathological classifications are insufficient to predict the cancer risk. Here we analyzed copy number alterations (CNA) to assess the molecular heterogeneity of vulvar lesions in relation to HPV and cancer risk. HPV-status and CNA by means of whole-genome next-generation shallow-sequencing were assessed in VSCC and VIN. The latter included VIN of women with associated VSCC (VIN VSCC ) and women who did not develop VSCC during follow-up (VIN no VSCC ). HPV-testing resulted in 41 HPV-positive (16 VIN VSCC , 14 VIN no VSCC , and 11 VSCC) and 24 HPV-negative (11 VIN VSCC and 13 VSCC) lesions. HPV-positive and -negative VSCC showed a partially overlapping pattern of recurrent CNA, including frequent gains of 3q and 8q. In contrast, amplification of 11q13/cyclinD1 was exclusively found in HPV-negative lesions. HPV-negative VIN VSCC had less CNA than HPV-negative VSCC (P = .009), but shared chromosome 8 alterations. HPV-positive VIN no VSCC had less CNA than VIN VSCC (P = .022). Interestingly, 1pq gain was detected in 81% of HPV-positive VIN VSCC and only in 21% of VIN no VSCC (P = .001). In conclusion, HPV-dependent and -independent vulvar carcinogenesis is characterized by distinct CNA patterns at the VIN stage, while more comparable patterns are present at the cancer stage. Cancer risk in VIN seems to be reflected by the extent of CNA, in particular chromosome 1 gain in HPV-positive cases., (© 2018 The Authors. Cancer Medicine published by John Wiley & Sons Ltd.)

The genomic complexity of profound copy number aberrations has prevented effective molecular stratification of ovarian cancers. Here, to decode this complexity, we derived copy number signatures from shallow whole-genome sequencing of 117 high-grade serous ovarian cancer (HGSOC) cases, which were validated on 527 independent cases. We show that HGSOC comprises a continuum of genomes shaped by multiple mutational processes that result in known patterns of genomic aberration. Copy number signature exposures at diagnosis predict both overall survival and the probability of platinum-resistant relapse. Measurement of signature exposures provides a rational framework to choose combination treatments that target multiple mutational processes.

Oncolytic adenoviruses represent a novel class of anticancer agents. Their efficacy in killing cancer cells is variable, suggesting that there is room for improvement. Host miRNAs have been shown to play important roles in susceptibility of cells to replication of different viruses. This study investigated if adenovirus replication in human prostate cancer cells is influenced by host cell miRNA expression. To this end, human miRNA expression in response to adenovirus infection was analyzed, and functional screens for lytic adenovirus replication were performed using synthetic miRNA mimic and inhibitor libraries. Adenovirus infection generally reduced miRNA expression. On top of this nonspecific interference with miRNA biogenesis, a set of miRNAs, including in particular miR-222, was found specifically reduced. Another set of miRNAs was found to promote adenovirus-induced death of prostate cancer cells. In most cases, this did not stimulate adenovirus propagation. The exception was miR-26b. Overexpression of miR-26b inhibited adenovirus-induced NF-κB activation, augmented adenovirus-mediated cell death, increased adenovirus progeny release, and promoted adenovirus propagation and spread in several human prostate cancer cell lines. This suggests that miR-26b is particularly useful to be combined with oncolytic adenovirus for more effective treatment of prostate cancer.

Tumor profiling of DNA alterations, i.e. gene point mutations, somatic copy number aberrations (CNAs) and structural variants (SVs), improves insight into the molecular pathology of cancer and clinical outcome. Here, associations between genomic aberrations and disease recurrence in stage II and III colon cancers were investigated. A series of 114 stage II and III microsatellite stable colon cancer samples were analyzed by high-resolution array-comparative genomic hybridization (array-CGH) to detect CNAs and CNA-associated chromosomal breakpoints (SVs). For 60 of these samples mutation status of APC, TP53, KRAS, PIK3CA, FBXW7, SMAD4, BRAF and NRAS was determined using targeted massive parallel sequencing. Loss of chromosome 18q12.1-18q12.2 occurred more frequently in tumors that relapsed than in relapse-free tumors (p < 0.001; FDR = 0.13). In total, 267 genes were recurrently affected by SVs (FDR < 0.1). CNAs and SVs were not associated with disease-free survival (DFS). Mutations in APC and TP53 were associated with increased CNAs. APC mutations were associated with poor prognosis in (5-fluorouracil treated) stage III colon cancers (p = 0.005; HR = 4.1), an effect that was further enhanced by mutations in MAPK pathway (KRAS, NRAS, BRAF) genes. We conclude that among multiple genomic alterations in CRC, strongest associations with clinical outcome were observed for common mutations in APC.

Chromosomal instability is a hallmark of cancer, but mitotic regulators are rarely mutated in tumors. Mutations in the condensin complexes, which restructure chromosomes to facilitate segregation during mitosis, are significantly enriched in cancer genomes, but experimental evidence implicating condensin dysfunction in tumorigenesis is lacking. We report that mice inheriting missense mutations in a condensin II subunit (Caph2 nes ) develop T-cell lymphoma. Before tumors develop, we found that the same Caph2 mutation impairs ploidy maintenance to a different extent in different hematopoietic cell types, with ploidy most severely perturbed at the CD4 + CD8 + T-cell stage from which tumors initiate. Premalignant CD4 + CD8 + T cells show persistent catenations during chromosome segregation, triggering DNA damage in diploid daughter cells and elevated ploidy. Genome sequencing revealed that Caph2 single-mutant tumors are near diploid but carry deletions spanning tumor suppressor genes, whereas P53 inactivation allowed Caph2 mutant cells with whole-chromosome gains and structural rearrangements to form highly aggressive disease. Together, our data challenge the view that mitotic chromosome formation is an invariant process during development and provide evidence that defective mitotic chromosome structure can promote tumorigenesis., (© 2016 Woodward et al.; Published by Cold Spring Harbor Laboratory Press.)

Aim: Gangliogliomas (GGs) and dysembryoplastic neuroepithelial tumours (DNTs) represent the most common histological entities within the spectrum of glioneuronal tumours (GNTs). The wide variability of morphological features complicates histological classification, including discrimination from prognostically distinct diffuse low-grade astrocytomas (AIIs). This study was performed to increase our understanding of these tumours., Methods: We studied chromosomal copy number aberrations (CNAs) by genome-wide sequencing in a large cohort of GNTs and linked these to comprehensive histological analysis and clinical characteristics. One hundred fourteen GNTs were studied: 50 GGs and 64 DNTs. Also, a data set of CNAs from 38 diffuse AIIs was included., Results: The most frequent CNAs in both GGs and DNTs were gains at chromosomes 5 and 7, often concurrent, and gain at chromosome 6. None of the CNAs was linked to histological subtype, immunohistochemical features or to clinical characteristics. Comparison of AIIs and diffuse GNTs revealed that gain at whole chromosome 5 is only observed in GNTs. CNA patterns indicative of chromothripsis were detected in three GNTs., Conclusion: We conclude that GNTs with diverse morphologies share molecular features, and our findings support the need to improve classification and differential diagnosis of tumour entities within the spectrum of GNTs, as well as their distinction from other gliomas., (© 2015 British Neuropathological Society.)

Somatic mutations in cervical intraepithelial neoplasia (CIN) are largely unknown. Here, we profiled 35 cervical carcinomas and 23 CIN grade 2/3 (CIN2/3) for mutations in 48 cancer-related genes using a Next Generation Sequencing-based cancer panel. PIK3CA exon 9 was the most frequently mutated locus in cervical carcinoma and the only mutated locus detected in CIN2/3. These PIK3CA exon 9 mutation findings were verified in a large, independent series (n = 647) covering all stages of cervical carcinogenesis using high resolution melting-guided Sanger sequencing. PIK3CA exon 9 mutation frequency was 37.1% (13/35; 95%CI 21.2-54.0%) in cervical carcinoma, and 2.4% (5/209; 95%CI 0.5-4.7%) in CIN3. No PIK3CA exon 9 mutations were detected in CIN2 (0/144), CIN1 (0/154) and normal cervix (0/105). In a third series of 46 CIN2/3 lesions from women with a known 5-year history of preceding high-risk human papillomavirus (hrHPV) infection, detection of PIK3CA exon 9 mutation was confined to 2 (5.4%; 95%CI 0.0-13.2%) CIN3 lesions with preceding hrHPV infection ≥5 years, and was absent in those with a short duration (<5 years) of preceding hrHPV infection. In conclusion, somatic mutation in PIK3CA represents a late event during cervical carcinogenesis, detected in a substantial subset of cervical carcinoma, but only in a minority of CIN3.

Background: Cancer is caused by somatic DNA alterations such as gene point mutations, DNA copy number aberrations (CNA) and structural variants (SVs). Genome-wide analyses of SVs in large sample series with well-documented clinical information are still scarce. Consequently, the impact of SVs on carcinogenesis and patient outcome remains poorly understood. This study aimed to perform a systematic analysis of genes that are affected by CNA-associated chromosomal breaks in colorectal cancer (CRC) and to determine the clinical relevance of recurrent breakpoint genes., Methods: Primary CRC samples of patients with metastatic disease from CAIRO and CAIRO2 clinical trials were previously characterized by array-comparative genomic hybridization. These data were now used to determine the prevalence of CNA-associated chromosomal breaks within genes across 352 CRC samples. In addition, mutation status of the commonly affected APC, TP53, KRAS, PIK3CA, FBXW7, SMAD4, BRAF and NRAS genes was determined for 204 CRC samples by targeted massive parallel sequencing. Clinical relevance was assessed upon stratification of patients based on gene mutations and gene breakpoints that were observed in >3% of CRC cases., Results: In total, 748 genes were identified that were recurrently affected by chromosomal breaks (FDR <0.1). MACROD2 was affected in 41% of CRC samples and another 169 genes showed breakpoints in >3% of cases, indicating that prevalence of gene breakpoints is comparable to the prevalence of well-known gene point mutations. Patient stratification based on gene breakpoints and point mutations revealed one CRC subtype with very poor prognosis., Conclusions: We conclude that CNA-associated chromosomal breaks within genes represent a highly prevalent and clinically relevant subset of SVs in CRC.

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., (Copyright © 2015 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserved.)

Detection of DNA copy number aberrations by shallow whole-genome sequencing (WGS) faces many challenges, including lack of completion and errors in the human reference genome, repetitive sequences, polymorphisms, variable sample quality, and biases in the sequencing procedures. Formalin-fixed paraffin-embedded (FFPE) archival material, the analysis of which is important for studies of cancer, presents particular analytical difficulties due to degradation of the DNA and frequent lack of matched reference samples. We present a robust, cost-effective WGS method for DNA copy number analysis that addresses these challenges more successfully than currently available procedures. In practice, very useful profiles can be obtained with ∼0.1× genome coverage. We improve on previous methods by first implementing a combined correction for sequence mappability and GC content, and second, by applying this procedure to sequence data from the 1000 Genomes Project in order to develop a blacklist of problematic genome regions. A small subset of these blacklisted regions was previously identified by ENCODE, but the vast majority are novel unappreciated problematic regions. Our procedures are implemented in a pipeline called QDNAseq. We have analyzed over 1000 samples, most of which were obtained from the fixed tissue archives of more than 25 institutions. We demonstrate that for most samples our sequencing and analysis procedures yield genome profiles with noise levels near the statistical limit imposed by read counting. The described procedures also provide better correction of artifacts introduced by low DNA quality than prior approaches and better copy number data than high-resolution microarrays at a substantially lower cost., (© 2014 Scheinin et al.; Published by Cold Spring Harbor Laboratory Press.)

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.

Purpose: Next generation DNA sequencing (NGS) holds promise for diagnostic applications, yet implementation in routine molecular pathology practice requires performance evaluation on DNA derived from routine formalin-fixed paraffin-embedded (FFPE) tissue specimens. The current study presents a comprehensive analysis of TruSeq Amplicon Cancer Panel-based NGS using a MiSeq Personal sequencer (TSACP-MiSeq-NGS) for somatic mutation profiling., Methods: TSACP-MiSeq-NGS (testing 212 hotspot mutation amplicons of 48 genes) and a data analysis pipeline were evaluated in a retrospective learning/test set approach (n = 58/n = 45 FFPE-tumor DNA samples) against 'gold standard' high-resolution-melting (HRM)-sequencing for the genes KRAS, EGFR, BRAF and PIK3CA. Next, the performance of the validated test algorithm was assessed in an independent, prospective cohort of FFPE-tumor DNA samples (n = 75)., Results: In the learning set, a number of minimum parameter settings was defined to decide whether a FFPE-DNA sample is qualified for TSACP-MiSeq-NGS and for calling mutations. The resulting test algorithm revealed 82% (37/45) compliance to the quality criteria and 95% (35/37) concordant assay findings for KRAS, EGFR, BRAF and PIK3CA with HRM-sequencing (kappa = 0.92; 95% CI = 0.81-1.03) in the test set. Subsequent application of the validated test algorithm to the prospective cohort yielded a success rate of 84% (63/75), and a high concordance with HRM-sequencing (95% (60/63); kappa = 0.92; 95% CI = 0.84-1.01). TSACP-MiSeq-NGS detected 77 mutations in 29 additional genes., Conclusion: TSACP-MiSeq-NGS is suitable for diagnostic gene mutation profiling in oncopathology.

Functional biomolecules, including small noncoding RNAs (ncRNAs), are released and transmitted between mammalian cells via extracellular vesicles (EVs), including endosome-derived exosomes. The small RNA composition in cells differs from exosomes, but underlying mechanisms have not been established. We generated small RNA profiles by RNA sequencing (RNA-seq) from a panel of human B cells and their secreted exosomes. A comprehensive bioinformatics and statistical analysis revealed nonrandomly distributed subsets of microRNA (miRNA) species between B cells and exosomes. Unexpectedly, 3' end adenylated miRNAs are relatively enriched in cells, whereas 3' end uridylated isoforms appear overrepresented in exosomes, as validated in naturally occurring EVs isolated from human urine samples. Collectively, our findings suggest that posttranscriptional modifications, notably 3' end adenylation and uridylation, exert opposing effects that may contribute, at least in part, to direct ncRNA sorting into EVs., (Copyright © 2014 The Authors. Published by Elsevier Inc. All rights reserved.)

Background: The disease course of patients with diffuse low-grade glioma is notoriously unpredictable. Temporal and spatially distinct samples may provide insight into the evolution of clinically relevant copy number aberrations (CNAs). The purpose of this study is to identify CNAs that are indicative of aggressive tumor behavior and can thereby complement the prognostically favorable 1p/19q co-deletion., Results: Genome-wide, 50 base pair single-end sequencing was performed to detect CNAs in a clinically well-characterized cohort of 98 formalin-fixed paraffin-embedded low-grade gliomas. CNAs are correlated with overall survival as an endpoint. Seventy-five additional samples from spatially distinct regions and paired recurrent tumors of the discovery cohort were analyzed to interrogate the intratumoral heterogeneity and spatial evolution. Loss of 10q25.2-qter is a frequent subclonal event and significantly correlates with an unfavorable prognosis. A significant correlation is furthermore observed in a validation set of 126 and confirmation set of 184 patients. Loss of 10q25.2-qter arises in a longitudinal manner in paired recurrent tumor specimens, whereas the prognostically favorable 1p/19q co-deletion is the only CNA that is stable across spatial regions and recurrent tumors., Conclusions: CNAs in low-grade gliomas display extensive intratumoral heterogeneity. Distal loss of 10q is a late onset event and a marker for reduced overall survival in low-grade glioma patients. Intratumoral heterogeneity and higher frequencies of distal 10q loss in recurrences suggest this event is involved in outgrowth to the recurrent tumor.

The somatic mutation burden in healthy white blood cells (WBCs) is not well known. Based on deep whole-genome sequencing, we estimate that approximately 450 somatic mutations accumulated in the nonrepetitive genome within the healthy blood compartment of a 115-yr-old woman. The detected mutations appear to have been harmless passenger mutations: They were enriched in noncoding, AT-rich regions that are not evolutionarily conserved, and they were depleted for genomic elements where mutations might have favorable or adverse effects on cellular fitness, such as regions with actively transcribed genes. The distribution of variant allele frequencies of these mutations suggests that the majority of the peripheral white blood cells were offspring of two related hematopoietic stem cell (HSC) clones. Moreover, telomere lengths of the WBCs were significantly shorter than telomere lengths from other tissues. Together, this suggests that the finite lifespan of HSCs, rather than somatic mutation effects, may lead to hematopoietic clonal evolution at extreme ages.

Background: Complex designs are common in (observational) clinical studies. Sequencing data for such studies are produced more and more often, implying challenges for the analysis, such as excess of zeros, presence of random effects and multi-parameter inference. Moreover, when sample sizes are small, inference is likely to be too liberal when, in a Bayesian setting, applying a non-appropriate prior or to lack power when not carefully borrowing information across features., Results: We show on microRNA sequencing data from a clinical cancer study how our software ShrinkBayes tackles the aforementioned challenges. In addition, we illustrate its comparatively good performance on multi-parameter inference for groups using a data-based simulation. Finally, in the small sample size setting, we demonstrate its high power and improved FDR estimation by use of Gaussian mixture priors that include a point mass., Conclusion: ShrinkBayes is a versatile software package for the analysis of count-based sequencing data, which is particularly useful for studies with small sample sizes or complex designs.

Background: Patients with human papillomavirus (HPV)-positive oropharyngeal squamous cell carcinomas (OPSCCs) have a better prognosis than patients with HPV-negative OPSCCs. Important factors contributing to this better prognosis are relatively low numbers of local/regional recurrences (LRRs) and second primary tumors (SPTs) in patients with HPV-positive OPSCC. These low numbers may be explained in addition by the absence of a 'field cancerization' effect, which is a cause of LRRs and SPTs in patients with HPV-negative OPSCC. We aimed to detect a possible 'field effect' in patients with HPV-positive OPSCC. As HPV is involved in the early stage of carcinogenesis in OPSCCs, its presence is considered a reliable marker for the detection of such a field effect. Therefore, the presence of transcriptionally active HPV was analyzed in the mucosa surrounding HPV-positive OPSCCs., Methods: We included 20 patients who were surgically treated for an HPV-positive OPSCC in the period 2000-2006. Of each patient, the formalin-fixed paraffin-embedded tumor sample and all available resection margins were collected. In total, 97 resection margins were investigated with an average of five resection margins per tumor. All samples were analyzed for the presence of tumor and the presence of transcriptionally active HPV by HPV16-E6-mRNA detection., Results: All tumors were HPV16-E6-mRNA positive. HPV16-E6-mRNA could be detected in the resection margins that contained tumor (n = 6). All tumor-negative resection margins (n = 91) scored negative for HPV16-E6-mRNA., Conclusions: In conclusion, transcriptional active HPV could not be detected in the mucosa surrounding an HPV-positive OPSCC, which suggests the absence of field effect. This observation may explain the lower number of LRRs and SPTs in HPV-positive patients., (© 2013 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

Background: Mosaic IDH1 mutations are described as the cause of metaphyseal chondromatosis with increased urinary excretion of D-2-hydroxyglutarate (MC-HGA), and mutations in IDH2 as the cause of D-2-hydroxyglutaric aciduria (D-2HGA) type II. Mosaicism for IDH2 mutations has not previously been reported as a cause of D-2HGA. Here we describe three cases: one MC-HGA case with IDH1 mosaic mutations, and two D-2HGA type II cases. In one D-2HGA case we identified mosaicism for an IDH2 mutation as the genetic cause of this disorder; the other D-2HGA case was caused by a heterozygous IDH2 mutation, while the unaffected mother was a mosaic carrier., Methods: We performed amplicon deep sequencing using the 454 GS Junior platform, next to Sanger sequencing, to identify and confirm mosaicism of IDH1 or IDH2 mutations in MC-HGA or D-2HGA, respectively., Results and Conclusions: We identified different mutant allele percentages in DNA samples derived from different tissues (blood vs fibroblasts). Furthermore, we found that mutant allele percentages of IDH1 decreased after more passages had occurred in fibroblast cell cultures. We describe a method for the detection and validation of mosaic mutations in IDH1 and IDH2, making quantification with laborious cloning techniques obsolete.

Objectives: Epidermal growth factor receptor (EGFR) mutations have been identified in lung adenocarcinomas and are associated with high response chance to EGFR tyrosine kinase inhibitors. EGFR mutations can be detected in tumour tissue, cytology specimens and blood from lung cancer patients. Thus far, EGFR mutation analysis has not been systematically demonstrated for sputum samples. The aim of the present study was to determine whether EGFR mutation analysis is attainable on sputum samples, employing different assays in a multicenter study., Materials and Methods: Sputum DNA from 10 lung cancer patients with confirmed EGFR mutation in their tumour tissue, 10 lung cancer patients without evidence of an EGFR mutation, and 10 patients with chronic obstructive pulmonary disease (COPD) was used for mutation analysis by Cycleave PCR, COLD-PCR, PangaeaBiotech SL Technology (PST), and High Resolution Melting, respectively. Targeted resequencing (TruSeq Amplicon Cancer Panel) and droplet digital PCR were additionally performed on the 10 samples with EGFR mutation., Results: Dependent on the assay, EGFR mutations could be detected in 30-50% of the sputum samples of patients with EGFR mutations. The different techniques revealed consistent results, with slightly higher sensitivity for PST. Neither the lung cancer patients without EGFR mutation nor the COPD controls tested positive for EGFR mutations in their sputum samples, indicating high clinical specificity of all assays., Conclusion: EGFR mutations can be detected in sputum samples from patients with EGFR-mutated non-small cell lung cancer, which may replace biopsy procedure for some patients., (Copyright © 2013 Elsevier Ireland Ltd. All rights reserved.)

Purpose: A gene expression signature, predictive for local recurrence after breast-conserving treatment, has previously been identified from a series of 165 young patients with breast cancer. We evaluated this signature on both another platform and an independent series, compared its performance with other published gene-sets, and investigated the gene expression profile of a larger data set., Experimental Design: Gene expression tumor profiles were obtained on 148 of the initial 165 Dutch patients and on an independent validation series of 195 French patients. Both unsupervised and supervised classifications were used to study the gene expression profile of the 343 breast cancers and to identify subgroups that differ for their risk of local recurrence., Results: The previous local recurrence signature was validated across platforms. However, when applied to the French patients, the signature did not reproduce its reported performance and did not better classify the patients than other published gene sets. Hierarchical clustering of all 343 breast cancers did not show any grouping reflecting local recurrence status. Genes related to proliferation were found differentially expressed between patients with or without local recurrence only in triple-negative tumors. Supervised classification revealed no significant gene set predictive for local recurrence or able to outperform classification based on clinical variables., Conclusions: Although the previously identified local recurrence signature was robust on another platform, we were neither able to validate it on an independent data set, nor able to define a strong gene expression classifier for local recurrence using a larger data set. We conclude that there are no significant differences in gene expression pattern in tumors from patients with and without local recurrence after breast-conserving treatment.

Fanconi anemia (FA) is a rare genetic instability syndrome characterized by developmental defects, bone marrow failure, and a high cancer risk. Fifteen genetic subtypes have been distinguished. The majority of patients (≈85%) belong to the subtypes A (≈60%), C (≈15%) or G (≈10%), while a minority (≈15%) is distributed over the remaining 12 subtypes. All subtypes seem to fit within the "classical" FA phenotype, except for D1 and N patients, who have more severe clinical symptoms. Since FA patients need special clinical management, the diagnosis should be firmly established, to exclude conditions with overlapping phenotypes. A valid FA diagnosis requires the detection of pathogenic mutations in a FA gene and/or a positive result from a chromosomal breakage test. Identification of the pathogenic mutations is also important for adequate genetic counselling and to facilitate prenatal or preimplantation genetic diagnosis. Here we describe and validate a comprehensive protocol for the molecular diagnosis of FA, based on massively parallel sequencing. We used this approach to identify BRCA2, FANCD2, FANCI and FANCL mutations in novel unclassified FA patients.

Dynamic modification of histone proteins plays a key role in regulating gene expression. However, histones themselves can also be dynamic, which potentially affects the stability of histone modifications. To determine the molecular mechanisms of histone turnover, we developed a parallel screening method for epigenetic regulators by analyzing chromatin states on DNA barcodes. Histone turnover was quantified by employing a genetic pulse-chase technique called RITE, which was combined with chromatin immunoprecipitation and high-throughput sequencing. In this screen, the NuB4/HAT-B complex, containing the conserved type B histone acetyltransferase Hat1, was found to promote histone turnover. Unexpectedly, the three members of this complex could be functionally separated from each other as well as from the known interacting factor and histone chaperone Asf1. Thus, systematic and direct interrogation of chromatin structure on DNA barcodes can lead to the discovery of genes and pathways involved in chromatin modification and dynamics., Competing Interests: The authors have declared that no competing interests exist.

Genome-wide RNA expression profiling has yielded tumor subtypes with strong predictive or prognostic value for a wide variety of cancers. Recently, for breast cancer two RNA expression classifiers have been adopted by the World Health Organization (WHO) and approved by the Food and Drug Administration (FDA). Also on the basis of DNA copy number profiles, tumor subtypes with different prognosis have been described, but have not yet led to clinical implementation. The genomic revolution caused by next generation sequencing of DNA samples presents additional mutation, balanced translocations, single-nucleotide polymorphisms (SNP), and copy neutral loss of heterozygosity data simultaneously. We foresee a further boost of the potential of DNA profiling in the clinic when these multidimensional DNA factors will be implemented. Here we evaluate the current stratification power with DNA copy numbers. In a training and validation approach using data of 400 published breast cancer samples, we show that a DNA copy number classifier accurately classifies RNA expression subtypes. We consider this an important step forward for clinical implementation of genomic subtyping using DNA and discuss the extra dimensions upcoming techniques will bring to the DNA palette., (©2011 AACR.)

MicroRNAs (miRNA) regulate many genes critical for tumorigenesis. We profiled miRNAs from 11 normal breast tissues, 17 noninvasive, 151 invasive breast carcinomas, and 6 cell lines by in-house-developed barcoded Solexa sequencing. miRNAs were organized in genomic clusters representing promoter-controlled miRNA expression and sequence families representing seed sequence-dependent miRNA target regulation. Unsupervised clustering of samples by miRNA sequence families best reflected the clustering based on mRNA expression available for this sample set. Clustering and comparative analysis of miRNA read frequencies showed that normal breast samples were separated from most noninvasive ductal carcinoma in situ and invasive carcinomas by increased miR-21 (the most abundant miRNA in carcinomas) and multiple decreased miRNA families (including miR-98/let-7), with most miRNA changes apparent already in the noninvasive carcinomas. In addition, patients that went on to develop metastasis showed increased expression of mir-423, and triple-negative breast carcinomas were most distinct from other tumor subtypes due to upregulation of the mir~17-92 cluster. However, absolute miRNA levels between normal breast and carcinomas did not reveal any significant differences. We also discovered two polymorphic nucleotide variations among the more abundant miRNAs miR-181a (T19G) and miR-185 (T16G), but we did not identify nucleotide variations expected for classical tumor suppressor function associated with miRNAs. The differentiation of tumor subtypes and prediction of metastasis based on miRNA levels is statistically possible but is not driven by deregulation of abundant miRNAs, implicating far fewer miRNAs in tumorigenic processes than previously suggested., (©2011 AACR.)

Gene expression is increasingly used for diagnostic, prognostic, and therapeutic purposes in clinical practice. We tested the hypothesis that volatile anesthetics (VA) affect gene expression of tumor cells. Cells from the neuronal cell line SH-SY5Y and from the breast cell line MCF-7 were exposed ex vivo to enflurane, isoflurane, desflurane, halothane, sevoflurane, or nitrous oxide. Microarray gene expression profiles were studied. We observed significant differences in gene expression levels of cell cultures and response in time when exposed to different VA. Some genes used for predictive genetic fingerprints for breast cancer were affected by VA. Our findings suggest that VA modulate gene expression in breast and brain tumor cell cultures in a unique and time-dependent manner.

The mechanism by which the immune system produces effector and memory T cells is largely unclear. To allow a large-scale assessment of the development of single naive T cells into different subsets, we have developed a technology that introduces unique genetic tags (barcodes) into naive T cells. By comparing the barcodes present in antigen-specific effector and memory T cell populations in systemic and local infection models, at different anatomical sites, and for TCR-pMHC interactions of different avidities, we demonstrate that under all conditions tested, individual naive T cells yield both effector and memory CD8+ T cell progeny. This indicates that effector and memory fate decisions are not determined by the nature of the priming antigen-presenting cell or the time of T cell priming. Instead, for both low and high avidity T cells, individual naive T cells have multiple fates and can differentiate into effector and memory T cell subsets.

The cyclin dependent kinase (CDK) inhibitors p15, p16, p21, and p27 are frequently deleted, silenced, or downregulated in many malignancies. Inactivation of CDK inhibitors predisposes mice to tumor development, showing that these genes function as tumor suppressors. Here, we describe high-throughput murine leukemia virus insertional mutagenesis screens in mice that are deficient for one or two CDK inhibitors. We retrieved 9,117 retroviral insertions from 476 lymphomas to define hundreds of loci that are mutated more frequently than expected by chance. Many of these loci are skewed toward a specific genetic context of predisposing germline and somatic mutations. We also found associations between these loci with gender, age of tumor onset, and lymphocyte lineage (B or T cell). Comparison of retroviral insertion sites with single nucleotide polymorphisms associated with chronic lymphocytic leukemia revealed a significant overlap between the datasets. Together, our findings highlight the importance of genetic context within large-scale mutation detection studies, and they show a novel use for insertional mutagenesis data in prioritizing disease-associated genes that emerge from genome-wide association studies.

The magnitude of antigen-specific CD8+ T cell responses is not fixed but correlates with the severity of infection. Although by definition T cell response size is the product of both the capacity to recruit naïve T cells (clonal selection) and their subsequent proliferation (clonal expansion), it remains undefined how these two factors regulate antigen-specific T cell responses. We determined the relative contribution of recruitment and expansion by labeling naïve T cells with unique genetic tags and transferring them into mice. Under disparate infection conditions with different pathogens and doses, recruitment of antigen-specific T cells was near constant and close to complete. Thus, naïve T cell recruitment is highly efficient, and the magnitude of antigen-specific CD8+ T cell responses is primarily controlled by clonal expansion.

T cells, as well as other cell types, are composed of phenotypically and functionally distinct subsets. However, for many of these populations it is unclear whether they develop from common or separate progenitors. To address such issues, we developed a novel approach, termed cellular barcoding, that allows the dissection of lineage relationships. We demonstrate that the labeling of cells with unique identifiers coupled to a microarray-based detection system can be used to analyze family relationships between the progeny of such cells. To exemplify the potential of this technique, we studied migration patterns of families of antigen-specific CD8(+) T cells in vivo. We demonstrate that progeny of individual T cells rapidly seed independent lymph nodes and that antigen-specific CD8(+) T cells present at different effector sites are largely derived from a common pool of precursors. These data show how locally primed T cells disperse and provide a technology for kinship analysis with wider utility.

p53 and p19(ARF) are tumor suppressors frequently mutated in human tumors. In a high-throughput screen in mice for mutations collaborating with either p53 or p19(ARF) deficiency, we identified 10,806 retroviral insertion sites, implicating over 300 loci in tumorigenesis. This dataset reveals 20 genes that are specifically mutated in either p19(ARF)-deficient, p53-deficient or wild-type mice (including Flt3, mmu-mir-106a-363, Smg6, and Ccnd3), as well as networks of significant collaborative and mutually exclusive interactions between cancer genes. Furthermore, we found candidate tumor suppressor genes, as well as distinct clusters of insertions within genes like Flt3 and Notch1 that induce mutants with different spectra of genetic interactions. Cross species comparative analysis with aCGH data of human cancer cell lines revealed known and candidate oncogenes (Mmp13, Slamf6, and Rreb1) and tumor suppressors (Wwox and Arfrp2). This dataset should prove to be a rich resource for the study of genetic interactions that underlie tumorigenesis.

Eberwine(-like) amplification of mRNA adds distinct 6-10 bp nucleotide stretches to the 5' end of amplified RNA transcripts. Analysis of over six thousand microarrays reveals that probes containing motifs complementary to these stretches are associated with aberrantly high signals up to a hundred fold the signal observed in unaffected probes. This is not observed when total RNA is used as target source. Different T7 primer sequences are used in different laboratories and platforms and consequently different T7 primer bias is observed in different datasets. This will hamper efforts to compare data sets across platforms.