Parent-specific copy number in paired tumor-normal studies using circular binary segmentation

Endocrine persistence in ER+ breast cancer is accompanied by metabolic vulnerability in oxidative phosphorylation

Despite adjuvant treatment with endocrine therapies, estrogen receptor-positive (ER+) breast cancers recur in a significant proportion of patients. Recurrences are attributable to clinically undetectable endocrine-tolerant persister cancer cells that retain tumor-forming potential. Therefore, strategies targeting such persister cells may prevent recurrent disease. Using CRISPR-Cas9 genome-wide knockout screening in ER+ breast cancer cells, we identified a survival mechanism involving metabolic reprogramming with reliance upon mitochondrial respiration in endocrine-tolerant persister cells. Quantitative proteomic profiling showed reduced levels of glycolytic proteins in persisters. Metabolic tracing of glucose revealed an energy-depleted state in persisters where oxidative phosphorylation was required to generate ATP. A phase II clinical trial was conducted to evaluate changes in mitochondrial markers in primary ER+/HER2-breast tumors induced by neoadjuvant endocrine therapy ( NCT04568616 ). In an analysis of tumor specimens from 32 patients, tumors exhibiting residual cell proliferation after aromatase inhibitor-induced estrogen deprivation with letrozole showed increased mitochondrial content. Genetic profiling and barcode lineage tracing showed that endocrine-tolerant persistence occurred stochastically without genetic predisposition. Mice bearing cell line- and patient-derived xenografts were used to measure the anti-tumor effects of mitochondrial complex I inhibition in the context of endocrine therapy. Pharmacological inhibition of complex I suppressed the tumor-forming potential of persisters and synergized with the anti-estrogen fulvestrant to induce regression of patient-derived xenografts. These findings indicate that mitochondrial metabolism is essential in endocrine-tolerant persister ER+ breast cancer cells and warrant the development of treatment strategies to leverage this vulnerability in the context of endocrine-sensitive disease.

Statement of Significance

Endocrine-tolerant persister cancer cells that survive endocrine therapy can cause recurrent disease. Persister cells exhibit increased energetic dependence upon mitochondria for survival and tumor re-growth potential.

High-Resolution Haplotyping of the PAH Gene Enables Early Gestation Noninvasive Prenatal Diagnosis of Phenylketonuria and Evolution Analysis of Recurrent Pathogenic Variations

BACKGROUND

The clinical performance of RHDO-based NIPD for PKU during early gestation remains under-evaluated. Furthermore, studies focused on SNP loci obtained by next-generation sequencing to analyze the genetic evolution of pathogenic variations in PKU is limited.

METHODS

Maternal peripheral blood, along with proband and paternal samples, was collected between 7 and 12 weeks of gestation. The PAH gene and surrounding high heterozygosity SNPs were targeted for enrichment and sequencing. Fetal genotypes were inferred using RHDO-based NIPD. High-resolution PAH haplotypes were used for the analysis of two common pathogenic variants in the Chinese population: c.728G>A and c.1238G>C.

RESULTS

Sixty one PKU families participated with an average fetal fraction of 6.08%. The median gestational age was 8+6 weeks. RHDO-based NIPD successfully identified fetal genotypes in 59 cases (96.72%, 59/62). Two cases failed because of insufficient informative SNPs. In addition, a recombination event was assessed in one fetus of 59 cases. Six, and three haplotypes were identified for c.728G>A(p.Arg243Gln) and c.1238G>C(p.Arg413Pro), respectively. Hap_3 and hap_8 were identified as the ancestral haplotypes for these pathogenic variants, with other haplotypes arising from mutations or recombination based on these ancestral haplotypes.

CONCLUSIONS

This study validates the feasibility of an RHDO-based assay for NIPD of PKU in early pregnancy and introduces its application in the demonstration of founder effects in recurrent pathogenic variations, offering new insights into the evolutionary analysis of PAH variations.

Deconstructing Intratumoral Heterogeneity through Multiomic and Multiscale Analysis of Serial Sections

Tumors may contain billions of cells, including distinct malignant clones and nonmalignant cell types. Clarifying the evolutionary histories, prevalence, and defining molecular features of these cells is essential for improving clinical outcomes, since intratumoral heterogeneity provides fuel for acquired resistance to targeted therapies. Here we present a statistically motivated strategy for deconstructing intratumoral heterogeneity through multiomic and multiscale analysis of serial tumor sections (MOMA). By combining deep sampling of IDH-mutant astrocytomas with integrative analysis of single-nucleotide variants, copy-number variants, and gene expression, we reconstruct and validate the phylogenies, spatial distributions, and transcriptional profiles of distinct malignant clones. By genotyping nuclei analyzed by single-nucleus RNA-seq for truncal mutations, we further show that commonly used algorithms for identifying cancer cells from single-cell transcriptomes may be inaccurate. We also demonstrate that correlating gene expression with tumor purity in bulk samples can reveal optimal markers of malignant cells and use this approach to identify a core set of genes that are consistently expressed by astrocytoma truncal clones, including AKR1C3, whose expression is associated with poor outcomes in several types of cancer. In summary, MOMA provides a robust and flexible strategy for precisely deconstructing intratumoral heterogeneity and clarifying the core molecular properties of distinct cellular populations in solid tumors.

First-trimester noninvasive prenatal diagnosis of seven facioscapulohumeral muscular dystrophy type 1 families using SNP-based amplicon sequencing: An earlier, rapid and safer way

The study is to explore the feasibility and value of SNP-based noninvasive prenatal diagnosis (NIPD) for facioscapulohumeral muscular dystrophy type 1 (FSHD1) in early pregnancy weeks. We prospectively collected seven FSHD1 families, with an average gestational age of 8+6. Among these seven couples, there were three affected FSHD1 mothers and four affected fathers. A multiplex-PCR panel comprising 402 amplicons was designed to selective enrich for highly heterozygous SNPs upstream of the DUX4 gene. Risk haplotype was constructed based on familial linkage analysis. Fetal genotypes were accurately inferred through relative haplotype dosage analysis using Bayes Factor. All tests were successfully completed in a single attempt, and no recombination events were detected. NIPD results were provided within a week, which is 4 weeks earlier than karyomapping and 7 weeks earlier than Bionano single-molecule optical mapping (BOM). Ultimately, five FSHD1 fetuses and two normal fetuses were successfully identified, with a 100% concordance rate with karyomapping and BOM. Therefore, SNP-based NIPD for FSHD1 was demonstrated to be feasible and accurate in early weeks of gestation, although the risk of recombination events cannot be completely eliminated. In the future, testing of more cases is still necessary to fully determine the clinical utility.

Pharmacogenomic analysis in adrenocortical carcinoma reveals genetic features associated with mitotane sensitivity and potential therapeutics

Background

Adrenocortical carcinoma (ACC) is an aggressive endocrine malignancy with limited therapeutic options. Treating advanced ACC with mitotane, the cornerstone therapy, remains challenging, thus underscoring the significance to predict mitotane response prior to treatment and seek other effective therapeutic strategies.

Objective

We aimed to determine the efficacy of mitotane via an in vitro assay using patient-derived ACC cells (PDCs), identify molecular biomarkers associated with mitotane response and preliminarily explore potential agents for ACC.

Methods

In vitro mitotane sensitivity testing was performed in 17 PDCs and high-throughput screening against 40 compounds was conducted in 8 PDCs. Genetic features were evaluated in 9 samples using exomic and transcriptomic sequencing.

Results

PDCs exhibited variable sensitivity to mitotane treatment. The median cell viability inhibition rate was 48.4% (IQR: 39.3-59.3%) and -1.2% (IQR: -26.4-22.1%) in responders (n=8) and non-responders (n=9), respectively. Median IC50 and AUC were remarkably lower in responders (IC50: 53.4 µM vs 74.7 µM, P<0.0001; AUC: 158.0 vs 213.5, P<0.0001). Genomic analysis revealed CTNNB1 somatic alterations were only found in responders (3/5) while ZNRF3 alterations only in non-responders (3/4). Transcriptomic profiling found pathways associated with lipid metabolism were upregulated in responder tumors whilst CYP27A1 and ABCA1 expression were positively correlated to in vitro mitotane sensitivity. Furthermore, pharmacologic analysis identified that compounds including disulfiram, niclosamide and bortezomib exhibited efficacy against PDCs.

Conclusion

ACC PDCs could be useful for testing drug response, drug repurposing and guiding personalized therapies. Our results suggested response to mitotane might be associated with the dependency on lipid metabolism. CYP27A1 and ABCA1 expression could be predictive markers for mitotane response, and disulfiram, niclosamide and bortezomib could be potential therapeutics, both warranting further investigation.

Deconstructing intratumoral heterogeneity through multiomic and multiscale analysis of serial sections

Tumors may contain billions of cells including distinct malignant clones and nonmalignant cell types. Clarifying the evolutionary histories, prevalence, and defining molecular features of these cells is essential for improving clinical outcomes, since intratumoral heterogeneity provides fuel for acquired resistance to targeted therapies. Here we present a statistically motivated strategy for deconstructing intratumoral heterogeneity through multiomic and multiscale analysis of serial tumor sections (MOMA). By combining deep sampling of IDH-mutant astrocytomas with integrative analysis of single-nucleotide variants, copy-number variants, and gene expression, we reconstruct and validate the phylogenies, spatial distributions, and transcriptional profiles of distinct malignant clones. By genotyping nuclei analyzed by single-nucleus RNA-seq for truncal mutations, we further show that commonly used algorithms for identifying cancer cells from single-cell transcriptomes may be inaccurate. We also demonstrate that correlating gene expression with tumor purity in bulk samples can reveal optimal markers of malignant cells and use this approach to identify a core set of genes that is consistently expressed by astrocytoma truncal clones, including AKR1C3, whose expression is associated with poor outcomes in several types of cancer. In summary, MOMA provides a robust and flexible strategy for precisely deconstructing intratumoral heterogeneity and clarifying the core molecular properties of distinct cellular populations in solid tumors.

Exploring factors impacting haplotype-based noninvasive prenatal diagnosis for single-gene recessive disorders

Haplotype-based noninvasive prenatal diagnosis (NIPD) is applicable for various recessive single-gene disorders in proband families. However, a comprehensive exploration of critical factors influencing the assay performance, such as fetal fraction, informative single nucleotide polymorphism (SNP) count, and recombination events, has yet to be performed. It is critical to identify key factors affecting NIPD performance, including its accuracy and success rate, and their impact on clinical diagnostics to guide clinical practice. We conducted a prospective study, recruiting 219 proband families with singleton pregnancies at risk for eight recessive single-gene disorders (Duchenne muscular dystrophy, spinal muscular atrophy, phenylketonuria, methylmalonic acidemia, hemophilia A, hemophilia B, non-syndromic hearing loss, and congenital adrenal hyperplasia) at 7-14 weeks of gestation. Haplotype-based NIPD was performed by evaluating the relative haplotype dosage (RHDO) in maternal circulation, and the results were validated via invasive prenatal diagnosis or newborn follow-ups. Among the 219 families, the median gestational age at first blood draw was 8+5 weeks. Initial testing succeeded for 190 families and failed for 29 due to low fetal fraction (16), insufficient informative SNPs (9), and homologous recombination near pathogenic variation (4). Among low fetal fraction families, successful testing was achieved for 11 cases after a redraw, while 5 remained inconclusive. Test failures linked to insufficient informative SNPs correlated with linkage disequilibrium near the genes, with F8 and MMUT exhibiting the highest associated failure rates (14.3% and 25%, respectively). Homologous recombination was relatively frequent around the DMD and SMN1 genes (8.8% and 4.8%, respectively) but led to detection failure in only 44.4% (4/9) of such cases. All NIPD results from the 201 successful families were consistent with invasive diagnostic findings or newborn follow-up. Fetal fraction, informative SNPs count, and homologous recombination are pivotal to NIPD performance. Redrawing blood effectively improves the success rate for low fetal fraction samples. However, informative SNPs count and homologous recombination rates vary significantly across genes, necessitating careful consideration in clinical practice. We have designed an in silico method based on linkage disequilibrium data to predict the number of informative SNPs. This can identify genomic regions where there might be an insufficient number of SNPs, thereby guiding panel design. With these factors properly accounted for, NIPD is highly accurate and reliable.

Genomic alterations driving precancerous to cancerous lesions in esophageal cancer development

Esophageal squamous cell carcinoma (ESCC) develops through a series of increasingly abnormal precancerous lesions. Previous studies have revealed the striking differences between normal esophageal epithelium and ESCC in copy number alterations (CNAs) and mutations in genes driving clonal expansion. However, due to limited data on early precancerous lesions, the timing of these transitions and which among them are prerequisites for malignant transformation remained unclear. Here, we analyze 1,275 micro-biopsies from normal esophagus, early and late precancerous lesions, and esophageal cancers to decipher the genomic alterations at each stage. We show that the frequency of TP53 biallelic inactivation increases dramatically in early precancerous lesion stage while CNAs and APOBEC mutagenesis substantially increase at late stages. TP53 biallelic loss is the prerequisite for the development of CNAs of genes in cell cycle, DNA repair, and apoptosis pathways, suggesting it might be one of the earliest steps initiating malignant transformation.

Haplotype-based noninvasive prenatal diagnosis of methylmalonic acidemia and the discovery of a recurrent pathogenic haplotype associated with c.609G>A

BACKGROUND

Early diagnosis and intervention are crucial for the prognosis of methylmalonic acidemia (MMA). However, research focused on early prenatal diagnosis of MMA is limited.

METHODS

A 161.89kb capture panel was designed for selectively enriching highly heterozygous SNPs. Fetal genotypes were inferred using relative haplotype dosage (RHDO) and Bayes factor, followed by invasive prenatal diagnosis (IPD) for validation. A core pathogenic haplotype associated with c.609G>A was identified based on the frequency differences between pathogenic and normal haplotypes.

RESULTS

We recruited 41 pregnancies at risk of MMA with a median gestational age of 8+2  weeks. The assay success rate of NIPD-MMA for maternal variants was 92.7% (38/41), and after incorporating the paternal result, the overall assay success rate reached 100% (41/41). All NIPD results were concordant with IPD. Notably, a core haplotype (hap_2), comprising 28 SNPs, demonstrates significant enrichment within pathogenic haplotypes bearing the c.609G>A variation. On average, c.609G>A carriers had 22.38 heterozygous loci within these 28 SNPs.

CONCLUSION

NIPD-MMA presents a viable choice for early, accurate, and safe prenatal diagnosis. Furthermore, the discovery of the recurrent core pathogenic haplotype provides a novel approach for haplotype phasing and has the potential for realizing proband-independent NIPD in the future.

Non-invasive prenatal testing (NIPT): Combination of copy number variant and gene analyses using an "in-house" target enrichment next generation sequencing-Solution for non-centralized NIPT laboratory?

OBJECTIVE

Recent studies have integrated copy number variant (CNV) and gene analysis using target enrichment. Here, we transferred this concept to our routine genetics laboratory, which is not linked to centralized non-invasive prenatal testing (NIPT) facilities.

METHOD

From a cohort of 100 pregnant women, 22 were selected for the analysis of maternal genomic DNA (gDNA) along with fetal cell-free DNA. Using targeted enrichment, 135 genes were analyzed, combined with aberrations of chromosomes 21, 18, 13, X, and Y. The data were subjected to specificity and sensitivity analyses, and correlated with the results from invasive testing methods.

RESULTS

The sensitivity/specificity was determined for the CNV analysis of chromosomes: 21 (80%/75%), 18 (-/82%), 13 (100%/67%), and Y (100%/100%). The gene detection was valid for maternal gDNA. However, for cell-free fetal DNA, it was not possible to determine the boundary between an artifact and a real sequence variant.

CONCLUSION

The target enrichment method combining CNV and gene detection seems feasible in a regular laboratory. However, this method can only be responsibly optimized with a sufficient number of controls and further validation on a strong bioinformatic background. The present results showed that NIPT should be performed in specialized centers, and that its introduction to isolated laboratories may not provide valid data.

PersonALL: a genetic scoring guide for personalized risk assessment in pediatric B-cell precursor acute lymphoblastic leukemia

BACKGROUND

Recurrent genetic lesions provide basis for risk assessment in pediatric acute lymphoblastic leukemia (ALL). However, current prognostic classifiers rely on a limited number of predefined sets of alterations.

METHODS

Disease-relevant copy number aberrations (CNAs) were screened genome-wide in 260 children with B-cell precursor ALL. Results were integrated with cytogenetic data to improve risk assessment.

RESULTS

CNAs were detected in 93.8% (n = 244) of the patients. First, cytogenetic profiles were combined with IKZF1 status (IKZF1normal, IKZF1del and IKZF1plus) and three prognostic subgroups were distinguished with significantly different 5-year event-free survival (EFS) rates, IKAROS-low (n = 215): 86.3%, IKAROS-medium (n = 27): 57.4% and IKAROS-high (n = 18): 37.5%. Second, contribution of genetic aberrations to the clinical outcome was assessed and an aberration-specific score was assigned to each prognostically relevant alteration. By aggregating the scores of aberrations emerging in individual patients, personalized cumulative values were calculated and used for defining four prognostic subgroups with distinct clinical outcomes. Two favorable subgroups included 60% of patients (n = 157) with a 5-year EFS of 96.3% (excellent risk, n = 105) and 87.2% (good risk, n = 52), respectively; while 40% of patients (n = 103) showed high (n = 74) or ultra-poor (n = 29) risk profile (5-year EFS: 67.4% and 39.0%, respectively).

CONCLUSIONS

PersonALL, our conceptually novel prognostic classifier considers all combinations of co-segregating genetic alterations, providing a highly personalized patient stratification.

Salivary gland cancer organoids are valid for preclinical genotype-oriented medical precision trials

Salivary gland cancers (SGCs) are heterogeneous tumors, and precision oncology represents a promising therapeutic approach; however, its impact on SGCs remains obscure. This study aimed to establish a translational model for testing molecular-targeted therapies by combining patient-derived organoids and genomic analyses of SGCs. We enrolled 29 patients, including 24 with SGCs and 5 with benign tumors. Resected tumors were subjected to organoid and monolayer cultures, as well as whole-exome sequencing. Organoid and monolayer cultures of SGCs were successfully established in 70.8% and 62.5% of cases, respectively. Organoids retained most histopathological and genetic profiles of their original tumors. In contrast, 40% of the monolayer-cultured cells did not harbor somatic mutations of their original tumors. The efficacy of molecular-targeted drugs tested on organoids depended on their oncogenic features. Organoids recapitulated the primary tumors and were useful for testing genotype-oriented molecular targeted therapy, which is valuable for precision medicine in patients with SGCs.

Oncogene Overlap Analysis of Circulating Cell-free Tumor DNA to Explore the Appropriate Criteria for Defining MET Copy Number-Driven Lung Cancer

INTRODUCTION

Defining clinically relevant MET amplification levels in non-small cell lung cancer (NSCLC) remains challenging. We hypothesize that oncogene overlap and MET amplicon size decline with increase in MET plasma copy number (pCN), thus enriching for MET-dependent states.

PATIENTS AND METHODS

We interrogated cell-free DNA NGS results of 16,782 patients with newly diagnosed advanced NSCLC to identify those with MET amplification as reported using Guardant360. Co-occurring genomic mutations and copy number alterations within each sample were evaluated. An exploratory method of adjusting for tumor fraction was also performed and amplicon size for MET was analyzed when available.

RESULTS

MET amplification was detected in 207 (1.2%) of samples. pCN ranged from 2.1 to 52.9. Of these, 43 (20.8%) had an overlapping oncogenic driver, including 23 (11.1%) METex14 skipping or other MET mutations. The degree of (non-MET) oncogene overlap decreased with increases in pCN. Patients with MET pCN ≥ 2.7 had lower rates of overlapping drivers compared to those with MET pCN < 2.7 (6.1% vs. 16.3%, P = .033). None of the 7 patients with pCN > 6.7 had an overlapping driver. After adjusting for tumor fraction, adjusted pCN (ApCN) was also lower for those with overlapping drivers than those without (median ApCN 4.9 vs. 7.3, P =.024). There was an inverse relationship between amplicon size and pCN.

CONCLUSIONS

We propose that a high MET pCN and/or ApCN, together with the absence of overlapping oncogenic drivers and small MET amplicon size, will enrich for patients most likely to derive benefit from MET targeted therapy.

Sensitization of Patient-Derived Colorectal Cancer Organoids to Photon and Proton Radiation by Targeting DNA Damage Response Mechanisms

Pathological complete response (pCR) has been correlated with overall survival in several cancer entities including colorectal cancer. Novel total neoadjuvant treatment (TNT) in rectal cancer has achieved pathological complete response in one-third of the patients. To define better treatment options for nonresponding patients, we used patient-derived organoids (PDOs) as avatars of the patient's tumor to apply both photon- and proton-based irradiation as well as single and combined chemo(radio)therapeutic treatments. While response to photon and proton therapy was similar, PDOs revealed heterogeneous responses to irradiation and different chemotherapeutic drugs. Radiotherapeutic response of the PDOs was significantly correlated with their ability to repair irradiation-induced DNA damage. The classical combination of 5-FU and irradiation could not sensitize radioresistant tumor cells. Ataxia-telangiectasia mutated (ATM) kinase was activated upon radiation, and by inhibition of this central sensor of DNA damage, radioresistant PDOs were resensitized. The study underlined the capability of PDOs to define nonresponders to irradiation and could delineate therapeutic approaches for radioresistant patients.

Molecular determinants of outcomes in relapsed or refractory mantle cell lymphoma treated with ibrutinib or temsirolimus in the MCL3001 (RAY) trial

Mantle cell lymphoma (MCL) is a rare, incurable lymphoma subtype characterized by heterogeneous outcomes. To better understand the clinical behavior and response to treatment, predictive biomarkers are needed. Using residual archived material from patients enrolled in the MCL3001 (RAY) study, we performed detailed analyses of gene expression and targeted genetic sequencing. This phase III clinical trial randomized patients with relapsed or refractory MCL to treatment with either ibrutinib or temsirolimus. We confirmed the prognostic capability of the gene expression proliferation assay MCL35 in this cohort treated with novel agents; it outperformed the simplified MCL International Prognostic Index in discriminating patients with different outcomes. Regardless of treatment arm, our data demonstrated that this assay captures the risk conferred by known biological factors, including increased MYC expression, blastoid morphology, aberrations of TP53, and truncated CCND1 3' untranslated region. We showed the negative impact of BIRC3 mutations/deletions on outcomes in this cohort and identified that deletion of chromosome 8p23.3 also negatively impacts survival. Restricted to patients with deletions/alterations in TP53, ibrutinib appeared to abrogate the deleterious impact on outcome. These data illustrate the potential to perform a molecular analysis of predictive biomarkers on routine patient samples that can meaningfully inform clinical practice.

Title: Noninvasive prenatal testing of hereditary colorectal cancer syndromes using cell-free DNA in maternal plasma

OBJECTIVE

This study aimed to establish a practical protocol for early noninvasive prenatal testing (NIPT) for fetuses at risk of Peutz-Jeghers syndrome (PJS) or familial adenomatous polyposis (FAP), two classical types of hereditary colorectal cancer syndromes, for risk evaluation and whole-life monitoring.

METHOD

Target enrichment was performed using hybridization probes coordinating the STK11 gene region and APC gene region, with 1,458 highly heterozygous SNPs included. Semitarget amplification random sequencing was used for large fragment deletion detection. For relative haplotype dosage (RHDO) analysis, haplotype construction was performed by SHAPEIT software, the CBS algorithm was used for recombination event calculation, and Bayes factor was used for the determination of whether the fetus was affected.

RESULTS

Haplotypes were successfully constructed in the nine recruited families with different pedigree characteristics, and the results for the RHDO analysis were consistent with the amniocentesis sampling detection results. The cell-free fetal DNA fraction can be detected as low as 2% in maternal plasma.

CONCLUSION

This is the first NIPT assay on hereditary colorectal cancer syndromes based upon RHDO analysis. This article is protected by copyright. All rights reserved.

Haplotype-Based Noninvasive Prenatal Diagnosis of 21 Families With Duchenne Muscular Dystrophy: Real-World Clinical Data in China

Noninvasive prenatal diagnosis (NIPD) of single-gene disorders has recently become the focus of clinical laboratories. However, reports on the clinical application of NIPD of Duchenne muscular dystrophy (DMD) are limited. This study aimed to evaluate the detection performance of haplotype-based NIPD of DMD in a real clinical environment. Twenty-one DMD families at 7-12 weeks of gestation were prospectively recruited. DNA libraries of cell-free DNA from the pregnant and genomic DNA from family members were captured using a custom assay for the enrichment of DMD gene exons and spanning single-nucleotide polymorphisms, followed by next-generation sequencing. Parental haplotype phasing was based on family linkage analysis, and fetal genotyping was inferred using the Bayes factor through target maternal plasma sequencing. Finally, the entire experimental process was promoted in the local clinical laboratory. We recruited 13 complete families, 6 families without paternal samples, and 2 families without probands in which daughter samples were collected. Two different maternal haplotypes were constructed based on family members in all 21 pedigrees at as early as 7 gestational weeks. Among the included families, the fetal genotypes of 20 families were identified at the first blood collection, and a second blood collection was performed for another family due to low fetal concentration. The NIPD result of each family was reported within 1 week. The fetal fraction in maternal cfDNA ranged from 1.87 to 11.68%. In addition, recombination events were assessed in two fetuses. All NIPD results were concordant with the findings of invasive prenatal diagnosis (chorionic villus sampling or amniocentesis). Exon capture and haplotype-based NIPD of DMD are regularly used for DMD genetic diagnosis, carrier screening, and noninvasive prenatal diagnosis in the clinic. Our method, haplotype-based early screening for DMD fetal genotyping via cfDNA sequencing, has high feasibility and accuracy, a short turnaround time, and is inexpensive in a real clinical environment.

Estimation of intrafamilial DNA contamination in family trio genome sequencing using deviation from Mendelian inheritance

With the increasing number of sequencing projects involving families, quality control tools optimized for family genome sequencing are needed. However, accurately quantifying contamination in a DNA mixture is particularly difficult when genetically related family members are the sources. We developed TrioMix, a maximum likelihood estimation (MLE) framework based on Mendel's law of inheritance, to quantify DNA mixture between family members in genome sequencing data of parent-offspring trios. TrioMix can accurately deconvolute any intrafamilial DNA contamination, including parent-offspring, sibling-sibling, parent-parent, and even multiple familial sources. In addition, TrioMix can be applied to detect genomic abnormalities that deviate from Mendelian inheritance patterns, such as uniparental disomy (UPD) and chimerism. A genome-wide depth and variant allele frequency plot generated by TrioMix facilitates tracing the origin of Mendelian inheritance deviations. We showed that TrioMix could accurately deconvolute genomes in both simulated and real data sets.

Determinants of anti-PD-1 response and resistance in clear cell renal cell carcinoma

ADAPTeR is a prospective, phase II study of nivolumab (anti-PD-1) in 15 treatment-naive patients (115 multiregion tumor samples) with metastatic clear cell renal cell carcinoma (ccRCC) aiming to understand the mechanism underpinning therapeutic response. Genomic analyses show no correlation between tumor molecular features and response, whereas ccRCC-specific human endogenous retrovirus expression indirectly correlates with clinical response. T cell receptor (TCR) analysis reveals a significantly higher number of expanded TCR clones pre-treatment in responders suggesting pre-existing immunity. Maintenance of highly similar clusters of TCRs post-treatment predict response, suggesting ongoing antigen engagement and survival of families of T cells likely recognizing the same antigens. In responders, nivolumab-bound CD8+ T cells are expanded and express GZMK/B. Our data suggest nivolumab drives both maintenance and replacement of previously expanded T cell clones, but only maintenance correlates with response. We hypothesize that maintenance and boosting of a pre-existing response is a key element of anti-PD-1 mode of action.

Whole chromosome loss and genomic instability in mouse embryos after CRISPR-Cas9 genome editing

Karyotype alterations have emerged as on-target complications from CRISPR-Cas9 genome editing. However, the events that lead to these karyotypic changes in embryos after Cas9-treatment remain unknown. Here, using imaging and single-cell genome sequencing of 8-cell stage embryos, we track both spontaneous and Cas9-induced karyotype aberrations through the first three divisions of embryonic development. We observe the generation of abnormal structures of the nucleus that arise as a consequence of errors in mitosis, including micronuclei and chromosome bridges, and determine their contribution to common karyotype aberrations including whole chromosome loss that has been recently reported after editing in embryos. Together, these data demonstrate that Cas9-mediated germline genome editing can lead to unwanted on-target side effects, including major chromosome structural alterations that can be propagated over several divisions of embryonic development.

The genomic profiling and MAMLD1 expression in human and canines with Cushing's disease

BACKGROUND

Cushing's disease (CD) is defined as hypercortisolemia caused by adrenocorticotropic hormone (ACTH)-secreting pituitary adenomas (corticotroph PA) that afflicts humans and dogs. In order to map common aberrant genomic features of CD between humans and dogs, we performed genomic sequencing and immunostaining on corticotroph PA.

METHODS

For inclusion, humans and dog were diagnosed with CD. Whole exome sequencing (WES) was conducted on 6 human corticotroph PA. Transcriptome RNA-Seq was performed on 6 human and 7 dog corticotroph PA. Immunohistochemistry (IHC) was complete on 31 human corticotroph PA. Corticotroph PA were compared with normal tissue and between species analysis were also performed.

RESULTS

Eight genes (MAMLD1, MNX1, RASEF, TBX19, BIRC5, TK1, GLDC, FAM131B) were significantly (P < 0.05) overexpressed across human and canine corticotroph PA. IHC revealed MAMLD1 to be positively (3+) expressed in the nucleus of ACTH-secreting tumor cells of human corticotroph PA (22/31, 70.9%), but absent in healthy human pituitary glands.

CONCLUSIONS

In this small exploratory cohort, we provide the first preliminary insights into profiling the genomic characterizations of human and dog corticotroph PA with respect to MAMLD1 overexpression, a finding of potential direct impact to CD microadenoma diagnosis. Our study also offers a rationale for potential use of the canine model in development of precision therapeutics.

Binding of the Treslin-MTBP Complex to Specific Regions of the Human Genome Promotes the Initiation of DNA Replication

The processes that control where higher eukaryotic cells initiate DNA replication throughout the genome are not understood clearly. In metazoans, the Treslin-MTBP complex mediates critical final steps in formation of the activated replicative helicase prior to initiation of replication. Here, we map the genome-wide distribution of the MTBP subunit of this complex in human cells. Our results indicate that MTBP binds to at least 30,000 sites in the genome. A majority of these sites reside in regions of open chromatin that contain transcriptional-regulatory elements (e.g., promoters, enhancers, and super-enhancers), which are known to be preferred areas for initiation of replication. Furthermore, many binding sites encompass two genomic features: a nucleosome-free DNA sequence (e.g., G-quadruplex DNA or AP-1 motif) and a nucleosome bearing histone marks characteristic of open chromatin, such as H3K4me2. Taken together, these findings indicate that Treslin-MTBP associates coordinately with multiple genomic signals to promote initiation of replication.

Kumagai and Dunphy show that Treslin-MTBP, activator of the replicative helicase, binds to at least 30,000 sites in the human genome. Many sites contain a nucleosome with active chromatin marks and nucleosome-free DNA (G-quadruplex or AP-1 site). Thus, Treslin-MTBP associates with multiple genomic elements to promote initiation of DNA replication.

Characterizing genetic intra-tumor heterogeneity across 2,658 human cancer genomes

Intra-tumor heterogeneity (ITH) is a mechanism of therapeutic resistance and therefore an important clinical challenge. However, the extent, origin, and drivers of ITH across cancer types are poorly understood. To address this, we extensively characterize ITH across whole-genome sequences of 2,658 cancer samples spanning 38 cancer types. Nearly all informative samples (95.1%) contain evidence of distinct subclonal expansions with frequent branching relationships between subclones. We observe positive selection of subclonal driver mutations across most cancer types and identify cancer type-specific subclonal patterns of driver gene mutations, fusions, structural variants, and copy number alterations as well as dynamic changes in mutational processes between subclonal expansions. Our results underline the importance of ITH and its drivers in tumor evolution and provide a pan-cancer resource of comprehensively annotated subclonal events from whole-genome sequencing data.

Genome-wide association study of germline copy number variations reveals an association with prostate cancer aggressiveness

Prostate cancer is a major health burden, being the second most commonly diagnosed malignancy in men worldwide. Overtreatment represents a major problem in prostate cancer therapy, leading to significant long-term quality-of-life effects for patients and a broad socio-ecological burden. Biomarkers that could facilitate risk stratification of prostate cancer aggressiveness at the time of diagnosis may help to guide clinical treatment decisions and reduce overtreatment. Previous research on genetic variations in prostate cancer has shown that germline copy number variations as well as somatic copy number alterations are commonly present in cancer patients, altering a greater portion of the cancer genome than any other type of genetic variation. To investigate the effect of germline copy number variations on cancer aggressiveness we have compared genome-wide screening data from genomic DNA isolated from the blood of 120 patients with aggressive prostate cancer, 231 patients with non-aggressive prostate cancer and 87 controls with benign prostatic hyperplasia from the Prostate Cancer Study of Austria biobank using the Affymetrix SNP 6.0 array. We could show that patients with an aggressive form of prostate cancer had a higher frequency of copy number variations [mean count of copy number segments (CNS) = 12.9, median count of CNS = 9] compared to patients with non-aggressive prostate cancer (mean count of CNS = 10.4, median count of CNS = 8) or control patients diagnosed with benign prostatic hyperplasia (mean count of CNS = 9.3, median count of CNS = 8). In general, we observed that copy number gain is a rarer event, compared to copy number loss within all three patient groups. Furthermore, we could show a significant effect of copy number losses located on chromosomes 8, 9 and 10 on prostate cancer aggressiveness (P = 0.040, P = 0.037 and P = 0.005, respectively). Applying a cross-validation analysis yielded an area under the curve of 0.63. Our study reports promising findings suggesting that copy number losses might play an important role in the establishment of novel biomarkers to predict prostate cancer aggressiveness at the time of diagnosis. Such markers could be used to facilitate risk stratification to reduce overtreatment of prostate cancer patients.

Mutational profiling of micro-dissected pre-malignant lesions from archived specimens

BACKGROUND

Systematic cancer screening has led to the increased detection of pre-malignant lesions (PMLs). The absence of reliable prognostic markers has led mostly to over treatment resulting in potentially unnecessary stress, or insufficient treatment and avoidable progression. Importantly, most mutational profiling studies have relied on PML synchronous to invasive cancer, or performed in patients without outcome information, hence limiting their utility for biomarker discovery. The limitations in comprehensive mutational profiling of PMLs are in large part due to the significant technical and methodological challenges: most PML specimens are small, fixed in formalin and paraffin embedded (FFPE) and lack matching normal DNA.

METHODS

Using test DNA from a highly degraded FFPE specimen, multiple targeted sequencing approaches were evaluated, varying DNA input amount (3-200 ng), library preparation strategy (BE: Blunt-End, SS: Single-Strand, AT: A-Tailing) and target size (whole exome vs. cancer gene panel). Variants in high-input DNA from FFPE and mirrored frozen specimens were used for PML-specific variant calling training and testing, respectively. The resulting approach was applied to profile and compare multiple regions micro-dissected (mean area 5 mm²) from 3 breast ductal carcinoma in situ (DCIS).

RESULTS

Using low-input FFPE DNA, BE and SS libraries resulted in 4.9 and 3.7 increase over AT libraries in the fraction of whole exome covered at 20x (BE:87%, SS:63%, AT:17%). Compared to high-confidence somatic mutations from frozen specimens, PML-specific variant filtering increased recall (BE:85%, SS:80%, AT:75%) and precision (BE:93%, SS:91%, AT:84%) to levels expected from sampling variation. Copy number alterations were consistent across all tested approaches and only impacted by the design of the capture probe-set. Applied to DNA extracted from 9 micro-dissected regions (8 PML, 1 normal epithelium), the approach achieved comparable performance, illustrated the data adequacy to identify candidate driver events (GATA3 mutations, ERBB2 or FGFR1 gains, TP53 loss) and measure intra-lesion genetic heterogeneity.

CONCLUSION

Alternate experimental and analytical strategies increased the accuracy of DNA sequencing from archived micro-dissected PML regions, supporting the deeper molecular characterization of early cancer lesions and achieving a critical milestone in the development of biology-informed prognostic markers and precision chemo-prevention strategies.

Evolution from adherent to suspension: systems biology of HEK293 cell line development

The need for new safe and efficacious therapies has led to an increased focus on biologics produced in mammalian cells. The human cell line HEK293 has bio-synthetic potential for human-like production attributes and is currently used for manufacturing of several therapeutic proteins and viral vectors. Despite the increased popularity of this strain we still have limited knowledge on the genetic composition of its derivatives. Here we present a genomic, transcriptomic and metabolic gene analysis of six of the most widely used HEK293 cell lines. Changes in gene copy and expression between industrial progeny cell lines and the original HEK293 were associated with cellular component organization, cell motility and cell adhesion. Changes in gene expression between adherent and suspension derivatives highlighted switching in cholesterol biosynthesis and expression of five key genes (RARG, ID1, ZIC1, LOX and DHRS3), a pattern validated in 63 human adherent or suspension cell lines of other origin.

Patient-derived cells from recurrent tumors that model the evolution of IDH-mutant glioma

Background

IDH-mutant lower-grade gliomas (LGGs) evolve under the selective pressure of therapy, but well-characterized patient-derived cells (PDCs) modeling evolutionary stages are lacking. IDH-mutant LGGs may develop therapeutic resistance associated with chemotherapy-driven hypermutation and malignant progression. The aim of this study was to establish and characterize PDCs, single-cell-derived PDCs (scPDCs), and xenografts (PDX) of IDH1-mutant recurrences representing distinct stages of tumor evolution.

Methods

We derived and validated cell cultures from IDH1-mutant recurrences of astrocytoma and oligodendroglioma. We used exome sequencing and phylogenetic reconstruction to examine the evolutionary stage represented by PDCs, scPDCs, and PDX relative to corresponding spatiotemporal tumor tissue and germline DNA. PDCs were also characterized for growth and tumor immortality phenotypes, and PDX were examined histologically.

Results

The integrated astrocytoma phylogeny revealed 2 independent founder clonal expansions of hypermutated (HM) cells in tumor tissue that are faithfully represented by independent PDCs. The oligodendroglioma phylogeny showed more than 4000 temozolomide-associated mutations shared among tumor samples, PDCs, scPDCs, and PDX, suggesting a shared monoclonal origin. The PDCs from both subtypes exhibited hallmarks of tumorigenesis, retention of subtype-defining genomic features, production of 2-hydroxyglutarate, and subtype-specific telomere maintenance mechanisms that confer tumor cell immortality. The oligodendroglioma PDCs formed infiltrative intracranial tumors with characteristic histology.

Conclusions

These PDCs, scPDCs, and PDX are unique and versatile community resources that model the heterogeneous clonal origins and functions of recurrent IDH1-mutant LGGs. The integrated phylogenies advance our knowledge of the complex evolution and immense mutational load of IDH1-mutant HM glioma.

Bioinformatics Tools and Resources for Cancer Immunotherapy Study

In recent years, cancer immunotherapy has emerged as a highly promising approach to treat patients with cancer, as the patient's own immune system is harnessed to attack cancer cells. However, the application of these approaches is still limited to a minority of patients with cancer and it is difficult to predict which patients will derive the greatest clinical benefit.One of the challenges faced by the biomedical community in the search of more effective biomarkers is the fact that translational research efforts involve collecting and accessing data at many different levels: from the type of material examined (e.g., cell line, animal models, clinical samples) to multiple data type (e.g., pharmacodynamic markers, genetic sequencing data) to the scale of a study (e.g., small preclinical study, moderate retrospective study on stored specimen sets, clinical trials with large cohorts).This chapter reviews several publicly available bioinformatics tools and data resources for high throughput molecular analyses applied to a range of data types, including those generated from microarray, whole-exome sequencing (WES), RNA-seq, DNA copy number, and DNA methylation assays, that are extensively used for integrative multidimensional data analysis and visualization.

A Unique Panel of Patient-Derived Cutaneous Squamous Cell Carcinoma Cell Lines Provides a Preclinical Pathway for Therapeutic Testing

BACKGROUND

Cutaneous squamous cell carcinoma (cSCC) incidence continues to rise with increasing morbidity and mortality, with limited treatment options for advanced disease. Future improvements in targeted therapy will rely on advances in genomic/transcriptomic understanding and the use of model systems for basic research. We describe here the panel of 16 primary and metastatic cSCC cell lines developed and characterised over the past three decades in our laboratory in order to provide such a resource for future preclinical research and drug screening.

METHODS

Primary keratinocytes were isolated from cSCC tumours and metastases, and cell lines were established. These were characterised using short tandem repeat (STR) profiling and genotyped by whole exome sequencing. Multiple in vitro assays were performed to document their morphology, growth characteristics, migration and invasion characteristics, and in vivo xenograft growth.

RESULTS

STR profiles of the cSCC lines allow the confirmation of their unique identity. Phylogenetic trees derived from exome sequence analysis of the matched primary and metastatic lines provide insight into the genetic basis of disease progression. The results of in vivo and in vitro analyses allow researchers to select suitable cell lines for specific experimentation.

CONCLUSIONS

There are few well-characterised cSCC lines available for widespread preclinical experimentation and drug screening. The described cSCC cell line panel provides a critical tool for in vitro and in vivo experimentation.

Directional allelic imbalance profiling and visualization from multi-sample data with RECUR

MOTIVATION

Genetic analysis of cancer regularly includes two or more samples from the same patient. Somatic copy number alterations leading to allelic imbalance (AI) play a critical role in cancer initiation and progression. Directional analysis and visualization of the alleles in imbalance in multi-sample settings allow for inference of recurrent mutations, providing insights into mutation rates, clonality and the genomic architecture and etiology of cancer.

RESULTS

The REpeat Chromosomal changes Uncovered by Reflection (RECUR) is an R application for the comparative analysis of AI profiles derived from SNP array and next-generation sequencing data. The algorithm accepts genotype calls and 'B allele' frequencies (BAFs) from at least two samples derived from the same individual. For a predefined set of genomic regions with AI, RECUR compares BAF values among samples. In the presence of AI, the expected value of a BAF can shift in two possible directions, reflecting an increased or decreased abundance of the maternal haplotype, relative to the paternal. The phenomenon of opposite haplotype shifts, or 'mirrored subclonal allelic imbalance', is a form of heterogeneity, and has been linked to clinico-pathological features of cancer. RECUR detects such genomic segments of opposite haplotypes in imbalance and plots BAF values for all samples, using a two-color scheme for intuitive visualization.

AVAILABILITY AND IMPLEMENTATION

RECUR is available as an R application. Source code and documentation are available at scheet.org.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Gene Co-expression Network and Copy Number Variation Analyses Identify Transcription Factors Associated With Multiple Myeloma Progression

Multiple myeloma (MM) has two clinical precursor stages of disease: monoclonal gammopathy of undetermined significance (MGUS) and smoldering multiple myeloma (SMM). However, the mechanism of progression is not well understood. Because gene co-expression network analysis is a well-known method for discovering new gene functions and regulatory relationships, we utilized this framework to conduct differential co-expression analysis to identify interesting transcription factors (TFs) in two publicly available datasets. We then used copy number variation (CNV) data from a third public dataset to validate these TFs. First, we identified co-expressed gene modules in two publicly available datasets each containing three conditions: normal, MGUS, and SMM. These modules were assessed for condition-specific gene expression, and then enrichment analysis was conducted on condition-specific modules to identify their biological function and upstream TFs. TFs were assessed for differential gene expression between normal and MM precursors, then validated with CNV analysis to identify candidate genes. Functional enrichment analysis reaffirmed known functional categories in MM pathology, the main one relating to immune function. Enrichment analysis revealed a handful of differentially expressed TFs between normal and either MGUS or SMM in gene expression and/or CNV. Overall, we identified four genes of interest (MAX, TCF4, ZNF148, and ZNF281) that aid in our understanding of MM initiation and progression.

Copy number alterations associated with clinical features in an underrepresented population with breast cancer

Background

As the most incident tumor among women worldwide, breast cancer is a heterogeneous disease. Tremendous efforts have been made to understand how tumor characteristics as histological type, molecular subtype, and tumor microenvironment collectively influence disease diagnosis to treatment, which impact outcomes. Differences between populations and environmental and cultural factors have impacts on the origin and evolution of the disease, as well as the therapeutic challenges that arise due to these factors. We, then, compared copy number variations (CNVs) in mucinous and nonmucinous luminal breast tumors from a Brazilian cohort to investigate major CNV imbalances in mucinous tumors versus non‐mucinous luminal tumors, taking into account their clinical and pathological features.

Methods

48 breast tumor samples and 48 matched control blood samples from Brazilian women were assessed for CNVs by chromosome microarray. Logistic regression and random forest models were used in order to assess CNVs in chromosomal regions from tumors.

Results

CNVs that were identified in chromosomes 1, 5, 8, 17, 19, and 21 classify tumors according to their histological type, ethnicity, disease stage, and familial history.

Conclusion

Copy number alterations described in this study provide a better understanding of the landscape of genomic aberrations in mucinous breast cancers that are associated with clinical features.

Defective homologous recombination DNA repair as therapeutic target in advanced chordoma

Chordomas are rare bone tumors with few therapeutic options. Here we show, using whole-exome and genome sequencing within a precision oncology program, that advanced chordomas (n = 11) may be characterized by genomic patterns indicative of defective homologous recombination (HR) DNA repair and alterations affecting HR-related genes, including, for example, deletions and pathogenic germline variants of BRCA2, NBN, and CHEK2. A mutational signature associated with HR deficiency was significantly enriched in 72.7% of samples and co-occurred with genomic instability. The poly(ADP-ribose) polymerase (PARP) inhibitor olaparib, which is preferentially toxic to HR-incompetent cells, led to prolonged clinical benefit in a patient with refractory chordoma, and whole-genome analysis at progression revealed a PARP1 p.T910A mutation predicted to disrupt the autoinhibitory PARP1 helical domain. These findings uncover a therapeutic opportunity in chordoma that warrants further exploration, and provide insight into the mechanisms underlying PARP inhibitor resistance.

Chordomas are rare bone tumors with limited therapeutic options. Here, the authors identify molecular alterations associated with defective homologous recombination DNA repair in advanced chordomas and report prolonged response in a patient treated with a PARP inhibitor, which later acquired resistance due to a newly gained PARP1 mutation.

Gliosarcoma Is Driven by Alterations in PI3K/Akt, RAS/MAPK Pathways and Characterized by Collagen Gene Expression Signature

Gliosarcoma is a very rare brain tumor reported to be a variant of glioblastoma (GBM), IDH-wildtype. While differences in molecular and histological features between gliosarcoma and GBM were reported, detailed information on the genetic background of this tumor is lacking. We intend to fill in this knowledge gap by the complex analysis of somatic mutations, indels, copy number variations, translocations and gene expression patterns in gliosarcomas. Using next generation sequencing, we determined somatic mutations, copy number variations (CNVs) and translocations in 10 gliosarcomas. Six tumors have been further subjected to RNA sequencing analysis and gene expression patterns have been compared to those of GBMs. We demonstrate that gliosarcoma bears somatic alterations in gene coding for PI3K/Akt (PTEN, PI3K) and RAS/MAPK (NF1, BRAF) signaling pathways that are crucial for tumor growth. Interestingly, the frequency of PTEN alterations in gliosarcomas was much higher than in GBMs. Aberrations of PTEN were the most frequent and occurred in 70% of samples. We identified genes differentially expressed in gliosarcoma compared to GBM (including collagen signature) and confirmed a difference in the protein level by immunohistochemistry. We found several novel translocations (including translocations in the RABGEF1 gene) creating potentially unfavorable combinations. Collected results on genetic alterations and transcriptomic profiles offer new insights into gliosarcoma pathobiology, highlight differences in gliosarcoma and GBM genetic backgrounds and point out to distinct molecular cues for targeted treatment.

Genomic analysis of the origins and evolution of multicentric diffuse lower-grade gliomas

Background

Rare multicentric lower-grade gliomas (LGGs) represent a unique opportunity to study the heterogeneity among distinct tumor foci in a single patient and to infer their origins and parallel patterns of evolution.

Methods

In this study, we integrate clinical features, histology, and immunohistochemistry for 4 patients with multicentric LGG, arising both synchronously and metachronously. For 3 patients we analyze the phylogeny of the lesions using exome sequencing, including one case with a total of 8 samples from the 2 lesions.

Results

One patient was diagnosed with multicentric isocitrate dehydrogenase 1 (IDH1) mutated diffuse astrocytomas harboring distinct IDH1 mutations, R132H and R132C; the latter mutation has been associated with Li-Fraumeni syndrome, which was subsequently confirmed in the patient's germline DNA and shown in additional cases with The Cancer Genome Atlas data. In another patient, phylogenetic analysis of synchronously arising grade II and grade III diffuse astrocytomas demonstrated a single shared mutation, IDH1 R132H, and revealed convergent evolution via non-overlapping mutations in ATRX and TP53. In 2 cases, there was divergent evolution of IDH1-mutated and 1p/19q-codeleted oligodendroglioma and IDH1-mutated and 1p/19q-intact diffuse astrocytoma, occurring synchronously in one case and metachronously in a second.

Conclusions

Each tumor in multicentric LGG cases may arise independently or may diverge very early in their development, presenting as genetically and histologically distinct tumors. Comprehensive sampling of these lesions can therefore significantly alter diagnosis and management. Additionally, somatic IDH1 R132C mutation in either multicentric or solitary LGG identifies unsuspected germline TP53 mutation, validating the limited number of published cases.

Quantitative analysis of somatically acquired and constitutive uniparental disomy in gastrointestinal cancers

Somatically acquired uniparental disomies (aUPDs) are frequent events in solid tumors and have been associated with cancer‐related genes. Studies assessing their functional consequences across several cancer types are therefore necessary. Here, we aimed at integrating aUPD profiles with the mutational status of cancer‐related genes in a tumor‐type specific manner. Using TCGA datasets for 1,032 gastrointestinal cancers, including colon (COAD), rectum (READ), stomach (STAD), esophageal adenocarcinoma (EAC) and esophageal squamous cell carcinoma (ESCC), we show a non‐random distribution of aUPD, suggesting the existence of a cancer‐specific landscape of aUPD events. Our analysis indicates that aUPD acts as a “second hit” in Knudson's model in order to achieve biallelic inactivation of tumor suppressor genes. In particular, APC, ARID1A and NOTCH1 were recurrently inactivated by the presence of homozygous mutation as a consequence of aUPD in COAD and READ, STAD and ESCC, respectively. Furthermore, while TP53 showed inactivation caused by aUPD at chromosome arm 17p across all tumor types, copy number losses at this genomic position were also frequent. By experimental and computationally inferring genome ploidy, we demonstrate that an increased number of aUPD events, both affecting the whole chromosome or segments of it, were present in highly aneuploid genomes compared to near‐diploid tumors. Finally, the presence of mosaic UPD was detected at a higher frequency in DNA extracted from peripheral blood lymphocytes of patients with colorectal cancer compared to healthy individuals. In summary, our study defines specific profiles of aUPD in gastrointestinal cancers and provides unequivocal evidence of their relevance in cancer.

What's new?

Somatically acquired uniparental disomies (aUPDs), in which two copies of a chromosome originate from the same parent, have been documented in various human cancers. Here, the authors examined the frequency of aUPDs in different gastrointestinal cancer types. Events involving aUPDs were found to occur at high incidence in gastrointestinal cancers and at increased frequency particularly in highly aneuploid genomes. The data also reveal a nonrandom distribution of aUPDs, with evidence of biallelic inactivation of tumor suppressor genes and activation of oncogenes in a tumor type‐specific manner. The findings suggest that aUPDs are functionally relevant in gastrointestinal malignancies.

Identifying disease-associated copy number variations by a doubly penalized regression model

Copy number variation (CNV) of DNA plays an important role in the development of many diseases. However, due to the irregularity and sparsity of the CNVs, studying the association between CNVs and a disease outcome or a trait can be challenging. Up to now, not many methods have been proposed in the literature for this problem. Most of the current researchers reply on an ad hoc two-stage procedure by first identifying CNVs in each individual genome and then performing an association test using these identified CNVs. This potentially leads to information loss and as a result a lower power to identify disease associated CNVs. In this article, we describe a new method that combines the two steps into a single coherent model to identify the common CNV across patients that are associated with certain diseases. We use a double penalty model to capture CNVs' association with both the intensities and the disease trait. We validate its performance in simulated datasets and a data example on platinum resistance and CNV in ovarian cancer genome.

Whole genome sequencing puts forward hypotheses on metastasis evolution and therapy in colorectal cancer

Incomplete understanding of the metastatic process hinders personalized therapy. Here we report the most comprehensive whole-genome study of colorectal metastases vs. matched primary tumors. 65% of somatic mutations originate from a common progenitor, with 15% being tumor- and 19% metastasis-specific, implicating a higher mutation rate in metastases. Tumor- and metastasis-specific mutations harbor elevated levels of BRCAness. We confirm multistage progression with new components ARHGEF7/ARHGEF33. Recurrently mutated non-coding elements include ncRNAs RP11-594N15.3, AC010091, SNHG14, 3’ UTRs of FOXP2, DACH2, TRPM3, XKR4, ANO5, CBL, CBLB, the latter four potentially dual protagonists in metastasis and efferocytosis-/PD-L1 mediated immunosuppression. Actionable metastasis-specific lesions include FAT1, FGF1, BRCA2, KDR, and AKT2-, AKT3-, and PDGFRA-3’ UTRs. Metastasis specific mutations are enriched in PI3K-Akt signaling, cell adhesion, ECM and hepatic stellate activation genes, suggesting genetic programs for site-specific colonization. Our results put forward hypotheses on tumor and metastasis evolution, and evidence for metastasis-specific events relevant for personalized therapy.

The evolution and genetic nature of metastatic lesions is not completely characterized. Here the authors perform a comprehensive whole-genome study of colorectal metastases in comparison to matched primary tumors and define a multistage progression model and metastasis-specific changes that, in part, are therapeutically actionable.

Metastatic adrenocortical carcinoma displays higher mutation rate and tumor heterogeneity than primary tumors

Adrenocortical cancer (ACC) is a rare cancer with poor prognosis and high mortality due to metastatic disease. All reported genetic alterations have been in primary ACC, and it is unknown if there is molecular heterogeneity in ACC. Here, we report the genetic changes associated with metastatic ACC compared to primary ACCs and tumor heterogeneity. We performed whole-exome sequencing of 33 metastatic tumors. The overall mutation rate (per megabase) in metastatic tumors was 2.8-fold higher than primary ACC tumor samples. We found tumor heterogeneity among different metastatic sites in ACC and discovered recurrent mutations in several novel genes. We observed 37–57% overlap in genes that are mutated among different metastatic sites within the same patient. We also identified new therapeutic targets in recurrent and metastatic ACC not previously described in primary ACCs.

Adrenocortical cancer (ACC) is a rarely diagnosed and aggressive cancer whose metastatic form has been scarcely studied. Here, the authors study primary and metastatic ACC to investigate genomic heterogeneity, discovering higher mutation rates in metastatic lesions and novel recurrent mutations.

Genome-Informed Targeted Therapy for Osteosarcoma

Osteosarcoma is a highly aggressive cancer for which treatment has remained essentially unchanged for more than 30 years. Osteosarcoma is characterized by widespread and recurrent somatic copy-number alterations (SCNA) and structural rearrangements. In contrast, few recurrent point mutations in protein-coding genes have been identified, suggesting that genes within SCNAs are key oncogenic drivers in this disease. SCNAs and structural rearrangements are highly heterogeneous across osteosarcoma cases, suggesting the need for a genome-informed approach to targeted therapy. To identify patient-specific candidate drivers, we used a simple heuristic based on degree and rank order of copy-number amplification (identified by whole-genome sequencing) and changes in gene expression as identified by RNA sequencing. Using patient-derived tumor xenografts, we demonstrate that targeting of patient-specific SCNAs leads to significant decrease in tumor burden, providing a road map for genome-informed treatment of osteosarcoma. SIGNIFICANCE: Osteosarcoma is treated with a chemotherapy regimen established 30 years ago. Although osteosarcoma is genomically complex, we hypothesized that tumor-specific dependencies could be identified within SCNAs. Using patient-derived tumor xenografts, we found a high degree of response for "genome-matched" therapies, demonstrating the utility of a targeted genome-informed approach.This article is highlighted in the In This Issue feature, p. 1.

The clinical implications of G1-G6 transcriptomic signature and 5-gene score in Korean patients with hepatocellular carcinoma

BACKGROUND

Efforts have been made to classify Hepatocellular Carcinoma (HCC) at surgically curable stages because molecular classification, which is prognostically informative, can accurately identify patients in need of additional early therapeutic interventions. Recently, HCC classification based French studies on the expression of 16 genes and 5 genes were proposed. In 16-gene classification, transcriptomic signatures (G1-G6) were used to classify HCC patients into clinical, genomic and pathway-specific subgroups. In 5-gene score classification, the good or poor prognosis of HCC patients was predicted. The patient's cohort in these studies was mainly from Caucasian and African populations. Here, we aimed to validate G1-G6 and 5-gene score signatures in 205 Korean HCC patients since genomic profiles of Korean patients are distinct from other regions.

METHODS

Integrated analyses using whole-exome sequencing, copy number variation and clinical data was performed against these two signatures to find statistical correlations. Kaplan-Meier, univariate and multivariate COX regression analysis were performed for Disease-Specific Survival (DSS) and Recurrence-Free Survival (RFS).

RESULTS

The G2 and G3 subgroups of transcriptomic signature were significantly associated with TP53 mutations while G5 and G6 subgroups were significantly associated with CTNNB1 mutations which is in concordance with original French studies. Similarly, the poor prognosis group of 5-gene score showed shorter DSS (p = 0.045) and early RFS (p = 0.023) as well as a significant association with microvascular invasion, tumor size (> 5 cm), elevated AFP levels, and RB1 mutations. However, the 5-gene score was not an independent prognostic factor for survival.

CONCLUSION

The G1-G6 and 5-gene signatures showed significant concordance between genetic profiles of Korean HCC patients and patients in original French studies. Thus, G1-G6 and 5-gene score signatures can be targeted as potential therapeutic biomarkers against HCC patients worldwide.

Spatial mutation patterns as markers of early colorectal tumor cell mobility

A growing body of evidence suggests that a subset of human cancers grows as single clonal expansions. In such a nearly neutral evolution scenario, it is possible to infer the early ancestral tree of a full-grown tumor. We hypothesized that early tree reconstruction can provide insights into the mobility phenotypes of tumor cells during their first few cell divisions. We explored this hypothesis by means of a computational multiscale model of tumor expansion incorporating the glandular structure of colorectal tumors. After calibrating the model to multiregional and single gland data from 19 human colorectal tumors using approximate Bayesian computation, we examined the role of early tumor cell mobility in shaping the private mutation patterns of the final tumor. The simulations showed that early cell mixing in the first tumor gland can result in side-variegated patterns where the same private mutations could be detected on opposite tumor sides. In contrast, absence of early mixing led to nonvariegated, sectional mutation patterns. These results suggest that the patterns of detectable private mutations in colorectal tumors may be a marker of early cell movement and hence the invasive and metastatic potential of the tumor at the start of the growth. In alignment with our hypothesis, we found evidence of early abnormal cell movement in 9 of 15 invasive colorectal carcinomas ("born to be bad"), but in none of 4 benign adenomas. If validated with a larger dataset, the private mutation patterns may be used for outcome prediction among screen-detected lesions with unknown invasive potential.

Genomic features of renal cell carcinoma with venous tumor thrombus

A venous tumor thrombus (VTT) is a potentially lethal complication of renal cell carcinoma (RCC) but virtually nothing is known about the underlying natural history. Based on our observation that venous thrombi contain significant numbers of viable tumor cells, we applied multiregion whole exome sequencing to a total of 37 primary tumor and VTT samples including normal tissue specimens from five consecutive patients. Our findings demonstrate mutational heterogeneity between primary tumor and VTT with 106 of 483 genes (22%) harboring functional SNVs and/or indels altered in either primary tumor or thrombus. Reconstruction of the clonal phylogeny showed clustering of tumor samples and VTT samples, respectively, in the majority of tumors. However, no new subclones were detected suggesting that pre-existing subclones of the primary tumor drive VTT formation. Importantly, we found several lines of evidence for “BRCAness” in a subset of tumors. These included mutations in genes that confer “BRCAness”, a mutational signature and an increase of small indels. Re-analysis of SNV calls from the TCGA KIRC-US cohort confirmed a high frequency of the “BRCAness” mutational signature AC3 in clear cell RCC. Our findings warrant further pre-clinical experiments and may lead to novel personalized therapies for RCC patients.

Novel variants identified with next-generation sequencing in Polish patients with cone-rod dystrophy

Purpose

The aim of this study was to identify the molecular genetic basis of cone-rod dystrophy in 18 unrelated families of Polish origin. Cone-rod dystrophy is one of the inherited retinal dystrophies, which constitute a highly heterogeneous group of disorders characterized by progressive dysfunction of photoreceptors and retinal pigment epithelium (RPE) cells.

Methods

The study group was composed of four groups of patients representing different Mendelian inheritance of the disease: autosomal dominant (AD), autosomal recessive (AR), X-linked recessive (XL), and autosomal recessive or X-linked recessive (AR/XL). The combined molecular strategy included Sanger sequencing of the RPGR-ORF15 gene (three families with XL and three families with the AR/XL mode of inheritance), mutation-specific microarray analysis of the ABCA4 gene (five families with the AR mode of inheritance and two families with the AR/XL mode of inheritance), targeted next-generation sequencing (NGS) of inherited retinal disease-associated (IRD) genes (seven families with the AD mode of inheritance and five families with the AR mode of inheritance), and whole exome sequencing, performed in select families who had been mutation-negative in the analysis with the targeted NGS panel (one family with the AD mode of inheritance, one family with the AR mode of inheritance, and two families with the AR/XL mode of inheritance).

Results

Based on this combined strategy, we managed to identify potentially causative variants in seven out of 18 families with CRD. Five of these variants are novel: c.3142_3143dupAA, p.(Glu1049Argfs*41) in the RPGR-ORF15 gene, two variants: c.1612delT, p.(Trp538Glyfs*15) and c.2389dupG, p.(Ile798Hisfs*20) in the PROM1 gene in one family, c.592A>C, p.(Ser198Arg) in the PRPH2 gene and the variant c.1691A>G, p.(Asp564Gly) in the ATF6 gene that we have already reported to be pathogenic. NGS on the IRD panel allowed the molecular basis of CRD to be identified in four out of 14 families with a total detection rate of 38%. WES allowed identification of the molecular genetic basis of CRD in one family.

Conclusions

This is the first report on the spectrum of disease genes and pathogenic variants causing CRD in the Polish population. The study presents five novel variants identified in four genes and therefore, broadens the spectrum of probable pathogenic variants associated with CRD.

TumorNext-Lynch-MMR: a comprehensive next generation sequencing assay for the detection of germline and somatic mutations in genes associated with mismatch repair deficiency and Lynch syndrome

The current algorithm for Lynch syndrome diagnosis is highly complex with multiple steps which can result in an extended time to diagnosis while depleting precious tumor specimens. Here we describe the analytical validation of a custom probe-based NGS tumor panel, TumorNext-Lynch-MMR, which generates a comprehensive genetic profile of both germline and somatic mutations that can accelerate and streamline the time to diagnosis and preserve specimen. TumorNext-Lynch-MMR can detect single nucleotide variants, small insertions and deletions in 39 genes that are frequently mutated in Lynch syndrome and colorectal cancer. Moreover, the panel provides microsatellite instability status and detects loss of heterozygosity in the five Lynch genes; MSH2, MSH6, MLH1, PMS2 and EPCAM. Clinical cases are described that highlight the assays ability to differentiate between somatic and germline mutations, precisely classify variants and resolve discordant cases.

The landscape of genomic alterations across childhood cancers

Pan-cancer analyses that examine commonalities and differences among various cancer types have emerged as a powerful way to obtain novel insights into cancer biology. Here we present a comprehensive analysis of genetic alterations in a pan-cancer cohort including 961 tumours from children, adolescents, and young adults, comprising 24 distinct molecular types of cancer. Using a standardized workflow, we identified marked differences in terms of mutation frequency and significantly mutated genes in comparison to previously analysed adult cancers. Genetic alterations in 149 putative cancer driver genes separate the tumours into two classes: small mutation and structural/copy-number variant (correlating with germline variants). Structural variants, hyperdiploidy, and chromothripsis are linked to TP53 mutation status and mutational signatures. Our data suggest that 7-8% of the children in this cohort carry an unambiguous predisposing germline variant and that nearly 50% of paediatric neoplasms harbour a potentially druggable event, which is highly relevant for the design of future clinical trials.

Universal informative CpG sites for inferring tumor purity from DNA methylation microarray data

Tumor purity is an intrinsic property of tumor samples and potentially has severe impact on many types of data analysis. We have previously developed a statistical method, InfiniumPurify, which could infer purity of a tumor sample given its tumor type (available in TCGA) or a set of informative CpG (iDMC) sites. However, in many clinical practices, researchers may focus on a specific type of tumor samples that is not included in TCGA, and samples which are too few to identify reliable iDMCs. This greatly restricts the application of InfiniumPurify in cancer research. In this paper, we proposed an updated version of InfiniumPurify (termed as uiInfiniumPurify) through identifying a universal set of iDMCs (uiDMCs) and redesigning the algorithm to determine hyper- and hypo-methylation status of each uiDMC. Through the application, we estimated tumor purities of 8830 tumor samples from TCGA. Result shows that our estimates are highly consistent with those by other available methods. Consequently, the updated uiInfiniumPurify, can be applied to a single sample (or a few samples) of interest whose tumor type is not included in TCGA. This characteristic will greatly broaden the application of uiInfiniumPurify in cancer research.

Integrative genomic and transcriptomic analysis of leiomyosarcoma

Leiomyosarcoma (LMS) is an aggressive mesenchymal malignancy with few therapeutic options. The mechanisms underlying LMS development, including clinically actionable genetic vulnerabilities, are largely unknown. Here we show, using whole-exome and transcriptome sequencing, that LMS tumors are characterized by substantial mutational heterogeneity, near-universal inactivation of TP53 and RB1, widespread DNA copy number alterations including chromothripsis, and frequent whole-genome duplication. Furthermore, we detect alternative telomere lengthening in 78% of cases and identify recurrent alterations in telomere maintenance genes such as ATRX, RBL2, and SP100, providing insight into the genetic basis of this mechanism. Finally, most tumors display hallmarks of "BRCAness", including alterations in homologous recombination DNA repair genes, multiple structural rearrangements, and enrichment of specific mutational signatures, and cultured LMS cells are sensitive towards olaparib and cisplatin. This comprehensive study of LMS genomics has uncovered key biological features that may inform future experimental research and enable the design of novel therapies.

The whole-genome landscape of medulloblastoma subtypes

Current therapies for medulloblastoma, a highly malignant childhood brain tumour, impose debilitating effects on the developing child, and highlight the need for molecularly targeted treatments with reduced toxicity. Previous studies have been unable to identify the full spectrum of driver genes and molecular processes that operate in medulloblastoma subgroups. Here we analyse the somatic landscape across 491 sequenced medulloblastoma samples and the molecular heterogeneity among 1,256 epigenetically analysed cases, and identify subgroup-specific driver alterations that include previously undiscovered actionable targets. Driver mutations were confidently assigned to most patients belonging to Group 3 and Group 4 medulloblastoma subgroups, greatly enhancing previous knowledge. New molecular subtypes were differentially enriched for specific driver events, including hotspot in-frame insertions that target KBTBD4 and ‘enhancer hijacking’ events that activate PRDM6. Thus, the application of integrative genomics to an extensive cohort of clinical samples derived from a single childhood cancer entity revealed a series of cancer genes and biologically relevant subtype diversity that represent attractive therapeutic targets for the treatment of patients with medulloblastoma.

Genomic analysis of 491 medulloblastoma samples, including methylation profiling of 1,256 cases, effectively assigns candidate drivers to most tumours across all molecular subgroups.

Medulloblastomas are highly malignant brain tumours that develop during childhood. Paul Northcott and colleagues analysed the whole-genome sequences of 491 medulloblastomas in order to characterize the genomic landscape across tumours and identify new drivers and mutational signatures. Their integrative genomic analyses, including methylation profiling of 1,256 medulloblastomas, identifies subgroup-specific driver mutations and suggests additional tumour subtypes. The authors assign driver mutations to a high proportion of the less well characterized Group 3 and Group 4, which together contribute to more than 60% of all medulloblastomas.