DOT1L regulates chromatin reorganization and gene expression during sperm differentiation.

Spermatozoa have a unique genome organization. Their chromatin is almost completely devoid of histones and is formed instead of protamines, which confer a high level of compaction and preserve paternal genome integrity until fertilization. Histone-to-protamine transition takes place in spermatids and is indispensable for the production of functional sperm. Here, we show that the H3K79-methyltransferase DOT1L controls spermatid chromatin remodeling and subsequent reorganization and compaction of the spermatozoon genome. Using a mouse model in which Dot1l is knocked-out (KO) in postnatal male germ cells, we found that Dot1l-KO sperm chromatin is less compact and has an abnormal content, characterized by the presence of transition proteins, immature protamine 2 forms and a higher level of histones. Proteomic and transcriptomic analyses performed on spermatids reveal that Dot1l-KO modifies the chromatin prior to histone removal and leads to the deregulation of genes involved in flagellum formation and apoptosis during spermatid differentiation. As a consequence of these chromatin and gene expression defects, Dot1l-KO spermatozoa have less compact heads and are less motile, which results in impaired fertility.

Molecular-based targeted therapies in patients with hepatocellular carcinoma and hepato-cholangiocarcinoma refractory to atezolizumab/bevacizumab.

BACKGROUND & AIMS: The "French Medicine Genomic program 2025" has been designed to give patients with cancers that are refractory to systemic treatments access to off-label therapies adapted to their specific genomic profile. Herein, we reported the results of this program in patients with advanced hepatocellular carcinoma (HCC) and hepato-cholangiocarcinoma (H-CCK). METHODS: In one center, all patients with HCC or H-CCK who progressed under atezolizumab/bevacizumab with available tumor frozen samples benefited from whole-genome/-exome and RNA sequencing. Targeted therapies were matched to genomic alterations following the recommendations of a molecular tumor board and radiological response and overall survival were assessed. RESULTS: Among 135 patients with HCC and H-CCK treated by atezolizumab/bevacizumab, 20 patients benefited from genomic analysis after progression (16 HCC; 4 H-CCK). Nineteen patients had analyzable data, 70% were male, median age was 57 years, 65% had metastatic disease and 45% had vascular invasion. Among these 19 patients, 14 patients (76%) harbored at least one actionable genomic alteration and 9/14 received an adapted targeted therapy (45%). One patient with H-CCK showing CDK4 amplification was treated with palbociclib and achieved a partial radiological response for 16 months. Another patient with H-CCK, high HER2 overexpression and a high homologous recombination score was treated with trastuzumab/olaparib and had stable disease. One patient with an HCC and bi-allelic inactivation of TSC2 achieved a complete radiological response under everolimus. The remaining six treated patients (all HCC) had progressive disease, including three patients treated with trametinib, two with everolimus and one with olaparib. CONCLUSION: Molecular-based guided therapy is feasible in patients with HCC/H-CCK progressing under atezolizumab/bevacizumab and may be useful in a small subset of patients. IMPACT AND IMPLICATIONS: The use of whole-genome/-exome and RNA sequencing in clinical practice has not been reported in patients with hepatocellular carcinoma and hepato-cholangiocarcinoma. Herein, we performed a pilot study which suggested that whole-genome/-exome and RNA sequencing is feasible on tumor biopsies from patients refractory to atezolizumab/bevacizumab, with a small subset of patients exhibiting at least one actionable genomic alteration and receiving an adapted targeted therapy. This proof-of-concept study suggests that this clinical strategy could benefit a small subset of patients. Finally, validation of this approach will be required in a larger cohort of patients.

Contribution of whole genome sequencing in the molecular diagnosis of mosaic partial deletion of the NF1 gene in neurofibromatosis type 1.

Neurofibromatosis type 1 (NF1) is an autosomal dominant disease with complete penetrance but highly variable expressivity. In most patients, Next Generation Sequencing (NGS) technologies allow the identification of a loss-of-function pathogenic variant in the NF1 gene, a negative regulator of the RAS-MAPK pathway. We describe the 5-year diagnosis wandering of a patient with a clear NF1 clinical diagnosis, but no molecular diagnosis using standard molecular technologies. The patient presented with a typical NF1 phenotype but NF1 targeted NGS, NF1 transcript analysis, MLPA, and array comparative genomic hybridization failed to reveal a genetic aberration. After 5 years of unsuccessful investigations, trio WGS finally identified a de novo mosaic (VAF ~ 14%) 24.6 kb germline deletion encompassing the promoter and first exon of NF1. This case report illustrates the relevance of WGS to detect structural variants including copy number variants that would be missed by alternative approaches. The identification of the causal pathogenic variant allowed a tailored genetic counseling with a targeted non-invasive prenatal diagnosis by detecting the deletion in plasmatic cell-free DNA from the proband's pregnant partner. This report clearly highlights the need to make WGS a clinically accessible test, offering a tremendous opportunity to identify a molecular diagnosis for otherwise unsolved cases.

Biallelic variants in ZNF142 lead to a syndromic neurodevelopmental disorder.

Biallelic variants of the gene encoding for the zinc-finger protein 142 (ZNF142) have recently been associated with intellectual disability (ID), speech impairment, seizures, and movement disorders in nine individuals from five families. In this study, we obtained phenotype and genotype information of 26 further individuals from 16 families. Among the 27 different ZNF142 variants identified in the total of 35 individuals only four were missense. Missense variants may give a milder phenotype by changing the local structure of ZF motifs as suggested by protein modeling; but this correlation should be validated in larger cohorts and pathogenicity of the missense variants should be investigated with functional studies. Clinical features of the 35 individuals suggest that biallelic ZNF142 variants lead to a syndromic neurodevelopmental disorder with mild to moderate ID, varying degrees of delay in language and gross motor development, early onset seizures, hypotonia, behavioral features, movement disorders, and facial dysmorphism. The differences in symptom frequencies observed in the unpublished individuals compared to those of published, and recognition of previously underemphasized facial features are likely to be due to the small sizes of the previous cohorts, which underlines the importance of larger cohorts for the phenotype descriptions of rare genetic disorders.

Circulating tumor DNA and circulating tumor cells in metastatic triple negative breast cancer patients.

Circulating tumor DNA (ctDNA) is a new circulating tumor biomarker which might be used as a prognostic biomarker in a way similar to circulating tumor cells (CTCs). Here, we used the high prevalence of TP53 mutations in triple negative breast cancer (TNBC) to compare ctDNA and CTC detection rates and prognostic value in metastatic TNBC patients. Forty patients were enrolled before starting a new line of treatment. TP53 mutations were characterized in archived tumor tissues and in plasma DNA using two next generation sequencing (NGS) platforms in parallel. Archived tumor tissue was sequenced successfully for 31/40 patients. TP53 mutations were found in 26/31 (84%) of tumor samples. The same mutation was detected in the matched plasma of 21/26 (81%) patients with an additional mutation found only in the plasma for one patient. Mutated allele fractions ranged from 2 to 70% (median 5%). The observed correlation between the two NGS approaches (R(2) = 0.903) suggested that ctDNA levels data were quantitative. Among the 27 patients with TP53 mutations, CTC count was ≥1 in 19 patients (70%) and ≥5 in 14 patients (52%). ctDNA levels had no prognostic impact on time to progression (TTP) or overall survival (OS), whereas CTC numbers were correlated with OS (p = 0.04) and marginally with TTP (p = 0.06). Performance status and elevated LDH also had significant prognostic impact. Here, absence of prognostic impact of baseline ctDNA level suggests that mechanisms of ctDNA release in metastatic TNBC may involve, beyond tumor burden, biological features that do not dramatically affect patient outcome.

A metabolic signature for NADSYN1-dependent congenital NAD deficiency disorder.

Nicotinamide adenine dinucleotide (NAD) is essential for embryonic development. To date, biallelic loss-of-function variants in 3 genes encoding nonredundant enzymes of the NAD de novo synthesis pathway - KYNU, HAAO, and NADSYN1 - have been identified in humans with congenital malformations defined as congenital NAD deficiency disorder (CNDD). Here, we identified 13 further individuals with biallelic NADSYN1 variants predicted to be damaging, and phenotypes ranging from multiple severe malformations to the complete absence of malformation. Enzymatic assessment of variant deleteriousness in vitro revealed protein domain-specific perturbation, complemented by protein structure modeling in silico. We reproduced NADSYN1-dependent CNDD in mice and assessed various maternal NAD precursor supplementation strategies to prevent adverse pregnancy outcomes. While for Nadsyn1+/- mothers, any B3 vitamer was suitable to raise NAD, preventing embryo loss and malformation, Nadsyn1-/- mothers required supplementation with amidated NAD precursors (nicotinamide or nicotinamide mononucleotide) bypassing their metabolic block. The circulatory NAD metabolome in mice and humans before and after NAD precursor supplementation revealed a consistent metabolic signature with utility for patient identification. Our data collectively improve clinical diagnostics of NADSYN1-dependent CNDD, provide guidance for the therapeutic prevention of CNDD, and suggest an ongoing need to maintain NAD levels via amidated NAD precursor supplementation after birth.

Nebula--a web-server for advanced ChIP-seq data analysis.

MOTIVATION: ChIP-seq consists of chromatin immunoprecipitation and deep sequencing of the extracted DNA fragments. It is the technique of choice for accurate characterization of the binding sites of transcription factors and other DNA-associated proteins. We present a web service, Nebula, which allows inexperienced users to perform a complete bioinformatics analysis of ChIP-seq data. RESULTS: Nebula was designed for both bioinformaticians and biologists. It is based on the Galaxy open source framework. Galaxy already includes a large number of functionalities for mapping reads and peak calling. We added the following to Galaxy: (i) peak calling with FindPeaks and a module for immunoprecipitation quality control, (ii) de novo motif discovery with ChIPMunk, (iii) calculation of the density and the cumulative distribution of peak locations relative to gene transcription start sites, (iv) annotation of peaks with genomic features and (v) annotation of genes with peak information. Nebula generates the graphs and the enrichment statistics at each step of the process. During Steps 3-5, Nebula optionally repeats the analysis on a control dataset and compares these results with those from the main dataset. Nebula can also incorporate gene expression (or gene modulation) data during these steps. In summary, Nebula is an innovative web service that provides an advanced ChIP-seq analysis pipeline providing ready-to-publish results. AVAILABILITY: Nebula is available at http://nebula.curie.fr/ SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A Melanoma-Tailored Next-Generation Sequencing Panel Coupled with a Comprehensive Analysis to Improve Routine Melanoma Genotyping.

BACKGROUND: Tumor molecular deciphering is crucial in clinical management. Pan-cancer next-generation sequencing panels have moved towards exhaustive molecular characterization. However, because of treatment resistance and the growing emergence of pharmacological targets, tumor-specific customized panels are needed to guide therapeutic strategies. OBJECTIVE: The objective of this study was to present such a customized next-generation sequencing panel in melanoma. METHODS: Melanoma patients with somatic molecular profiling performed as part of routine care were included. High-throughput sequencing was performed with a melanoma tailored next-generation sequencing panel of 64 genes involved in molecular classification, prognosis, theranostic, and therapeutic resistance. Single nucleotide variants and copy number variations were screened, and a comprehensive molecular analysis identified clinically relevant alterations. RESULTS: Four hundred and twenty-one melanoma cases were analyzed (before any treatment initiation for 94.8% of patients). After bioinformatic prioritization, we uncovered 561 single nucleotide variants, 164 copy number variations, and four splice-site mutations. At least one alteration was detected in 368 (87.4%) lesions, with BRAF, NRAS, CDKN2A, CCND1, and MET as the most frequently altered genes. Among patients with BRAFV600 mutated melanoma, 44.5% (77 of 173) harbored at least one concurrent alteration driving potential resistance to mitogen-activated protein kinase inhibitors. In patients with RAS hotspot mutated lesions and in patients with neither BRAFV600 nor RAS hotspot mutations, alterations constituting potential pharmacological targets were found in 56.9% (66 of 116) and 47.7% (63 of 132) of cases, respectively. CONCLUSIONS: Our tailored next-generation sequencing assay coupled with a comprehensive analysis may improve therapeutic management in a significant number of patients with melanoma. Updating such a panel and implementing multi-omic approaches will further enhance patients' clinical management.

Heterogeneity of neuroblastoma cell identity defined by transcriptional circuitries.

Neuroblastoma is a tumor of the peripheral sympathetic nervous system, derived from multipotent neural crest cells (NCCs). To define core regulatory circuitries (CRCs) controlling the gene expression program of neuroblastoma, we established and analyzed the neuroblastoma super-enhancer landscape. We discovered three types of identity in neuroblastoma cell lines: a sympathetic noradrenergic identity, defined by a CRC module including the PHOX2B, HAND2 and GATA3 transcription factors (TFs); an NCC-like identity, driven by a CRC module containing AP-1 TFs; and a mixed type, further deconvoluted at the single-cell level. Treatment of the mixed type with chemotherapeutic agents resulted in enrichment of NCC-like cells. The noradrenergic module was validated by ChIP-seq. Functional studies demonstrated dependency of neuroblastoma with noradrenergic identity on PHOX2B, evocative of lineage addiction. Most neuroblastoma primary tumors express TFs from the noradrenergic and NCC-like modules. Our data demonstrate a previously unknown aspect of tumor heterogeneity relevant for neuroblastoma treatment strategies.

Bioinformatics for precision medicine in oncology: principles and application to the SHIVA clinical trial.

Precision medicine (PM) requires the delivery of individually adapted medical care based on the genetic characteristics of each patient and his/her tumor. The last decade witnessed the development of high-throughput technologies such as microarrays and next-generation sequencing which paved the way to PM in the field of oncology. While the cost of these technologies decreases, we are facing an exponential increase in the amount of data produced. Our ability to use this information in daily practice relies strongly on the availability of an efficient bioinformatics system that assists in the translation of knowledge from the bench towards molecular targeting and diagnosis. Clinical trials and routine diagnoses constitute different approaches, both requiring a strong bioinformatics environment capable of (i) warranting the integration and the traceability of data, (ii) ensuring the correct processing and analyses of genomic data, and (iii) applying well-defined and reproducible procedures for workflow management and decision-making. To address the issues, a seamless information system was developed at Institut Curie which facilitates the data integration and tracks in real-time the processing of individual samples. Moreover, computational pipelines were developed to identify reliably genomic alterations and mutations from the molecular profiles of each patient. After a rigorous quality control, a meaningful report is delivered to the clinicians and biologists for the therapeutic decision. The complete bioinformatics environment and the key points of its implementation are presented in the context of the SHIVA clinical trial, a multicentric randomized phase II trial comparing targeted therapy based on tumor molecular profiling versus conventional therapy in patients with refractory cancer. The numerous challenges faced in practice during the setting up and the conduct of this trial are discussed as an illustration of PM application.

Breakpoint features of genomic rearrangements in neuroblastoma with unbalanced translocations and chromothripsis.

Neuroblastoma is a pediatric cancer of the peripheral nervous system in which structural chromosome aberrations are emblematic of aggressive tumors. In this study, we performed an in-depth analysis of somatic rearrangements in two neuroblastoma cell lines and two primary tumors using paired-end sequencing of mate-pair libraries and RNA-seq. The cell lines presented with typical genetic alterations of neuroblastoma and the two tumors belong to the group of neuroblastoma exhibiting a profile of chromothripsis. Inter and intra-chromosomal rearrangements were identified in the four samples, allowing in particular characterization of unbalanced translocations at high resolution. Using complementary experiments, we further characterized 51 rearrangements at the base pair resolution that revealed 59 DNA junctions. In a subset of cases, complex rearrangements were observed with templated insertion of fragments of nearby sequences. Although we did not identify known particular motifs in the local environment of the breakpoints, we documented frequent microhomologies at the junctions in both chromothripsis and non-chromothripsis associated breakpoints. RNA-seq experiments confirmed expression of several predicted chimeric genes and genes with disrupted exon structure including ALK, NBAS, FHIT, PTPRD and ODZ4. Our study therefore indicates that both non-homologous end joining-mediated repair and replicative processes may account for genomic rearrangements in neuroblastoma. RNA-seq analysis allows the identification of the subset of abnormal transcripts expressed from genomic rearrangements that may be involved in neuroblastoma oncogenesis.