Single-cell multiomics reveals the interplay of clonal evolution and cellular plasticity in hepatoblastoma.

Hepatoblastomas (HB) display heterogeneous cellular phenotypes that influence the clinical outcome, but the underlying mechanisms are poorly understood. Here, we use a single-cell multiomic strategy to unravel the molecular determinants of this plasticity. We identify a continuum of HB cell states between hepatocytic (scH), liver progenitor (scLP) and mesenchymal (scM) differentiation poles, with an intermediate scH/LP population bordering scLP and scH areas in spatial transcriptomics. Chromatin accessibility landscapes reveal the gene regulatory networks of each differentiation pole, and the sequence of transcription factor activations underlying cell state transitions. Single-cell mapping of somatic alterations reveals the clonal architecture of each tumor, showing that each genetic subclone displays its own range of cellular plasticity across differentiation states. The most scLP subclones, overexpressing stem cell and DNA repair genes, proliferate faster after neo-adjuvant chemotherapy. These results highlight how the interplay of clonal evolution and epigenetic plasticity shapes the potential of HB subclones to respond to chemotherapy.

Molecular characterization of corona radiata cells from patients with diminished ovarian reserve using microarray and microfluidic-based gene expression profiling.

BACKGROUND: Diminished ovarian reserve (DOR) is one of the causes of infertility in young women. In this prospective study, gene expression profiling (GEP) of corona radiata cells (CRC) was performed to identify genes deregulated in DOR patients. METHODS: Microarray-based GEP of CRC isolated from eight women undergoing IVF was performed to identify genes differentially expressed between patients with normal ovarian reserve and DOR patients. Microfluidic-based quantitative RT-PCR assays were used to validate selected transcripts on 40 independent patients. A principal component analysis was used to identify more homogeneous subgroups of DOR patients. In silico analyses focusing on cis-regulation were performed to refine the interactions between patient's biological characteristics and their GEP. RESULTS: Forty-eight transcripts were differentially expressed, including CXXC finger protein 5 (CXXC5), forkhead box C1 (FOXC1) (down-regulated in DOR) as well as connective tissue growth factor (CTGF), follistatin-like 3 (FSTL3), prostaglandin-endoperoxide synthase 2 (PTGS2) and suppressor of cytokine signaling 2 (SOCS2) (up-regulated in DOR). According to these transcripts, two DOR patients' subgroups (DOR Gr1 and Gr2) were identified. In DOR Gr2 patients, C-terminal domain 2 (CITED2), CTGF, growth arrest-specific 1 (GAS1), insulin receptor substrate 2 (IRS2), PTGS2, SOCS2 and Versican (VCAN) were expressed at significantly higher levels and CXXC5, FOXC1, guanylate-binding protein 2 (GBP2) and zinc finger MIZ-domain containing 1 (ZMIZ1) at significantly lower levels. Higher baseline estradiol (E(2)) levels were observed in DOR Gr2 patients (P < 0.006). The in silico analyses suggested that all 11 genes differentially expressed between DOR Gr1 and DOR Gr2 subgroups could be transcriptional targets of estrogen. CONCLUSIONS: Despite small sample size limitations, 12 genes deregulated in the CRC of DOR patients were identified, which could be involved in DOR pathogenesis. A DOR patient's subgroup with high baseline E(2) levels and deregulated estrogen-responsive genes was also identified.

An Accessible and Unique Insight into Metastasis Mutational Content Through Whole-exome Sequencing of Circulating Tumor Cells in Metastatic Prostate Cancer.

BACKGROUND: Genomic analysis of circulating tumor cells (CTCs) could provide a unique and accessible representation of tumor diversity but remains hindered by technical challenges associated with CTC rarity and heterogeneity. OBJECTIVE: To evaluate CTCs as surrogate samples for genomic analyses in metastatic castration-resistant prostate cancer (mCRPC). DESIGN, SETTING, AND PARTICIPANTS: Three isolation strategies (filter laser-capture microdissection, self-seeding microwell chips, and fluorescence-activated cell sorting) were developed to capture CTCs with various epithelial and mesenchymal phenotypes and isolate them at the single-cell level. Whole-genome amplification (WGA) and WGA quality control were performed on 179 CTC samples, matched metastasis biopsies, and negative controls from 11 patients. All patients but one were pretreated with enzalutamide or abiraterone. Whole-exome sequencing (WES) of 34 CTC samples, metastasis biopsies, and negative controls were performed for seven patients. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: WES of CTCs was rigorously qualified in terms of percentage coverage at 10× depth, allelic dropout, and uncovered regions. Shared somatic mutations between CTCs and matched metastasis biopsies were identified. A customized approach based on determination of mutation rates for CTC samples was developed for identification of CTC-exclusive mutations. RESULTS AND LIMITATIONS: Shared mutations were mostly detected in epithelial CTCs and were recurrent. For two patients for whom a deeper analysis was performed, a few CTCs were sufficient to represent half to one-third of the mutations in the matched metastasis biopsy. CTC-exclusive mutations were identified in both epithelial and nonepithelial CTCs and affected cytoskeleton, invasion, DNA repair, and cancer-driver genes. Some 41% of CTC-exclusive mutations had a predicted deleterious impact on protein function. Phylogenic relationships between CTCs with distinct phenotypes were evidenced. CONCLUSIONS: CTCs can provide unique insight into metastasis mutational diversity and reveal undiagnosed genomic aberrations in matched metastasis biopsies. PATIENT SUMMARY: Our results demonstrate the clinical potential of circulating tumor cells to provide insight into metastatic events that could be critical to target using precision medicine.

Robust physical methods that enrich genomic regions identical by descent for linkage studies: confirmation of a locus for osteogenesis imperfecta.

BACKGROUND: The monogenic disease osteogenesis imperfecta (OI) is due to single mutations in either of the collagen genes ColA1 or ColA2, but within the same family a given mutation is accompanied by a wide range of disease severity. Although this phenotypic variability implies the existence of modifier gene variants, genome wide scanning of DNA from OI patients has not been reported. Promising genome wide marker-independent physical methods for identifying disease-related loci have lacked robustness for widespread applicability. Therefore we sought to improve these methods and demonstrate their performance to identify known and novel loci relevant to OI. RESULTS: We have improved methods for enriching regions of identity-by-descent (IBD) shared between related, afflicted individuals. The extent of enrichment exceeds 10- to 50-fold for some loci. The efficiency of the new process is shown by confirmation of the identification of the Col1A2 locus in osteogenesis imperfecta patients from Amish families. Moreover the analysis revealed additional candidate linkage loci that may harbour modifier genes for OI; a locus on chromosome 1q includes COX-2, a gene implicated in osteogenesis. CONCLUSION: Technology for physical enrichment of IBD loci is now robust and applicable for finding genes for monogenic diseases and genes for complex diseases. The data support the further investigation of genetic loci other than collagen gene loci to identify genes affecting the clinical expression of osteogenesis imperfecta. The discrimination of IBD mapping will be enhanced when the IBD enrichment procedure is coupled with deep resequencing.

Technical Validation of a Reverse-Transcription Quantitative Polymerase Chain Reaction In Vitro Diagnostic Test for the Determination of MiR-31-3p Expression Levels in Formalin-Fixed Paraffin-Embedded Metastatic Colorectal Cancer Tumor Specimens.

MiR-31-3p expression has been shown to be a predictive biomarker for response to anti-epithelial growth factor receptor therapy in patients with RAS wild-type metastatic colorectal cancer (mCRC). To aid in the quantification of miR-31-3p expression in formalin-fixed paraffin-embedded (FFPE) primary tumor samples from patients with mCRC, a reverse-transcription quantitative polymerase chain reaction (RT-qPCR) assay was developed and validated. Assay development included the identification of a microRNA reference standard and the determination of an appropriate relative quantification cutoff for differentiating low versus high miR-31-3p expression. Sample specimens for the validation studies included both FFPE slides and shavings. Polymerase chain reaction (PCR) efficiency and linearity, analytical sensitivity and specificity, assay robustness, reproducibility, and accuracy were demonstrated across a number of test conditions and differing quantitative PCR platforms. The data from this study provide evidence as to the feasibility of quantifying the expression of miR-31-3p from FFPE tumor tissue using a standardized RT-qPCR assay.

Whole exome sequencing for determination of tumor mutation load in liquid biopsy from advanced cancer patients.

Tumor mutation load (TML) has been proposed as a biomarker of patient response to immunotherapy in several studies. TML is usually determined by tumor biopsy DNA (tDNA) whole exome sequencing (WES), therefore TML evaluation is limited by informative biopsy availability. Circulating cell free DNA (cfDNA) provided by liquid biopsy is a surrogate specimen to biopsy for molecular profiling. Nevertheless performing WES on DNA from plasma is technically challenging and the ability to determine tumor mutation load from liquid biopsies remains to be demonstrated. In the current study, WES was performed on cfDNA from 32 metastatic patients of various cancer types included into MOSCATO 01 (NCT01566019) and/or MATCHR (NCT02517892) molecular triage trials. Results from targeted gene sequencing (TGS) and WES performed on cfDNA were compared to results from tumor tissue biopsy. In cfDNA samples, WES mutation detection sensitivity was 92% compared to targeted sequencing (TGS). When comparing cfDNA-WES to tDNA-WES, mutation detection sensitivity was 53%, consistent with previously published prospective study comparing cfDNA-TGS to tDNA-TGS. For samples in which presence of tumor DNA was confirmed in cfDNA, tumor mutation load from liquid biopsy was correlated with tumor biopsy. Taken together, this study demonstrated that liquid biopsy may be applied to determine tumor mutation load. Qualification of liquid biopsy for interpretation is a crucial point to use cfDNA for mutational load estimation.

Massively parallel DNA sequencing from routinely processed cytological smears.

BACKGROUND: Data generated by next-generation sequencing technologies have a pivotal role in precision medicine. These high-throughput techniques are preferentially performed on fresh tissue, but there is an increasing need for protocols adapted to materials derived from formalin-fixed, paraffin-embedded tissue and cytology specimens. METHODS: The aim of this work was to show that cytological material collected from archival smears processed for routine diagnoses could be used for massively parallel sequencing and array-based genomic analysis for further studies. RESULTS: As a proof of concept, data obtained from May-Grünwald Giemsa- and Diff-Quik-stained archival smears were shown to be in keeping with those obtained from matched frozen controls. CONCLUSIONS: The quality of DNA extracted from routinely processed smears is compatible with the multitargeted sequencing of a large series of genes of interest with methods such as array-based genomic analysis and whole-exome sequencing.

An innovative strategy for the molecular diagnosis of Usher syndrome identifies causal biallelic mutations in 93% of European patients.

Usher syndrome (USH), the most prevalent cause of hereditary deafness-blindness, is an autosomal recessive and genetically heterogeneous disorder. Three clinical subtypes (USH1-3) are distinguishable based on the severity of the sensorineural hearing impairment, the presence or absence of vestibular dysfunction, and the age of onset of the retinitis pigmentosa. A total of 10 causal genes, 6 for USH1, 3 for USH2, and 1 for USH3, and an USH2 modifier gene, have been identified. A robust molecular diagnosis is required not only to improve genetic counseling, but also to advance gene therapy in USH patients. Here, we present an improved diagnostic strategy that is both cost- and time-effective. It relies on the sequential use of three different techniques to analyze selected genomic regions: targeted exome sequencing, comparative genome hybridization, and quantitative exon amplification. We screened a large cohort of 427 patients (139 USH1, 282 USH2, and six of undefined clinical subtype) from various European medical centers for mutations in all USH genes and the modifier gene. We identified a total of 421 different sequence variants predicted to be pathogenic, about half of which had not been previously reported. Remarkably, we detected large genomic rearrangements, most of which were novel and unique, in 9% of the patients. Thus, our strategy led to the identification of biallelic and monoallelic mutations in 92.7% and 5.8% of the USH patients, respectively. With an overall 98.5% mutation characterization rate, the diagnosis efficiency was substantially improved compared with previously reported methods.

Functional variants of POC5 identified in patients with idiopathic scoliosis.

Idiopathic scoliosis (IS) is a spine deformity that affects approximately 3% of the population. The underlying causes of IS are not well understood, although there is clear evidence that there is a genetic component to the disease. Genetic mapping studies suggest high genetic heterogeneity, but no IS disease-causing gene has yet been identified. Here, genetic linkage analyses combined with exome sequencing identified a rare missense variant (p.A446T) in the centriolar protein gene POC5 that cosegregated with the disease in a large family with multiple members affected with IS. Subsequently, the p.A446T variant was found in an additional set of families with IS and in an additional 3 cases of IS. Moreover, POC5 variant p.A455P was present and linked to IS in one family and another rare POC5 variant (p.A429V) was identified in an additional 5 cases of IS. In a zebrafish model, expression of any of the 3 human IS-associated POC5 variant mRNAs resulted in spine deformity, without affecting other skeletal structures. Together, these findings indicate that mutations in the POC5 gene contribute to the occurrence of IS.

SDH mutations establish a hypermethylator phenotype in paraganglioma.

Paragangliomas are neuroendocrine tumors frequently associated with mutations in RET, NF1, VHL, and succinate dehydrogenase (SDHx) genes. Methylome analysis of a large paraganglioma cohort identified three stable clusters, associated with distinct clinical features and mutational status. SDHx-related tumors displayed a hypermethylator phenotype, associated with downregulation of key genes involved in neuroendocrine differentiation. Succinate accumulation in SDH-deficient mouse chromaffin cells led to DNA hypermethylation by inhibition of 2-OG-dependent histone and DNA demethylases and established a migratory phenotype reversed by decitabine treatment. Epigenetic silencing was particularly severe in SDHB-mutated tumors, potentially explaining their malignancy. Finally, inactivating FH mutations were identified in the only hypermethylated tumor without SDHx mutations. These findings emphasize the interplay between the Krebs cycle, epigenomic changes, and cancer.

Association of TALS developmental disorder with defect in minor splicing component U4atac snRNA.

The spliceosome, a ribonucleoprotein complex that includes proteins and small nuclear RNAs (snRNAs), catalyzes RNA splicing through intron excision and exon ligation to produce mature messenger RNAs, which, in turn serve as templates for protein translation. We identified four point mutations in the U4atac snRNA component of the minor spliceosome in patients with brain and bone malformations and unexplained postnatal death [microcephalic osteodysplastic primordial dwarfism type 1 (MOPD 1) or Taybi-Linder syndrome (TALS); Mendelian Inheritance in Man ID no. 210710]. Expression of a subgroup of genes, possibly linked to the disease phenotype, and minor intron splicing were affected in cell lines derived from TALS patients. Our findings demonstrate a crucial role of the minor spliceosome component U4atac snRNA in early human development and postnatal survival.