Feasibility to use whole-genome sequencing as a sole diagnostic method to detect genomic aberrations in pediatric B-cell acute lymphoblastic leukemia.

Introduction: The suitability of whole-genome sequencing (WGS) as the sole method to detect clinically relevant genomic aberrations in B-cell acute lymphoblastic leukemia (ALL) was investigated with the aim of replacing current diagnostic methods. Methods: For this purpose, we assessed the analytical performance of 150 bp paired-end WGS (90x leukemia/30x germline). A set of 88 retrospective B-cell ALL samples were selected to represent established ALL subgroups as well as ALL lacking stratifying markers by standard-of-care (SoC), so-called B-other ALL. Results: Both the analysis of paired leukemia/germline (L/N)(n=64) as well as leukemia-only (L-only)(n=88) detected all types of aberrations mandatory in the current ALLTogether trial protocol, i.e., aneuploidies, structural variants, and focal copy-number aberrations. Moreover, comparison to SoC revealed 100% concordance and that all patients had been assigned to the correct genetic subgroup using both approaches. Notably, WGS could allocate 35 out of 39 B-other ALL samples to one of the emerging genetic subgroups considered in the most recent classifications of ALL. We further investigated the impact of high (90x; n=58) vs low (30x; n=30) coverage on the diagnostic yield and observed an equally perfect concordance with SoC; low coverage detected all relevant lesions. Discussion: The filtration of the WGS findings with a short list of genes recurrently rearranged in ALL was instrumental to extract the clinically relevant information efficiently. Nonetheless, the detection of DUX4 rearrangements required an additional customized analysis, due to multiple copies of this gene embedded in the highly repetitive D4Z4 region. We conclude that the diagnostic performance of WGS as the standalone method was remarkable and allowed detection of all clinically relevant genomic events in the diagnostic setting of B-cell ALL.

Insights into cellular behavior and micromolecular communication in urothelial micrografts.

Autologous micrografting is a technique currently applied within skin wound healing, however, the potential use for surgical correction of other organs with epithelial lining, including the urinary bladder, remains largely unexplored. Currently, little is known about the micrograft expansion potential and the micromolecular events that occur in micrografted urothelial cells. In this study, we aimed to evaluate the proliferative potential of different porcine urothelial micrograft sizes in vitro, and, furthermore, to explore how urothelial micrografts communicate and which microcellular events are triggered. We demonstrated that increased tissue fragmentation subsequently potentiated the yield of proliferative cells and the cellular expansion potential, which confirms, that the micrografting principles of skin epithelium also apply to uroepithelium. Furthermore, we targeted the expression of the extracellular signal-regulated kinase (ERK) pathway and demonstrated that ERK activation occurred predominately at the micrograft borders and that ERK inhibition led to decreased urothelial migration and proliferation. Finally, we successfully isolated extracellular vesicles from the micrograft culture medium and evaluated their contents and relevance within various enriched biological processes. Our findings substantiate the potential of applying urothelial micrografting in future tissue-engineering models for reconstructive urological surgery, and, furthermore, highlights certain mechanisms as potential targets for future wound healing treatments.

Mosaic Deletions of Known Genes Explain Skeletal Dysplasias With High and Low Bone Mass.

Mosaicism, a state in which an individual has two or more genetically distinct populations of cells in the body, can be difficult to detect because of either mild or atypical clinical presentation and limitations in the commonly used detection methods. Knowledge of the role of mosaicism is limited in many skeletal disorders, including osteopathia striata with cranial sclerosis (OSCS) and cleidocranial dysplasia (CCD). We used whole-genome sequencing (WGS) with coverage >40× to identify the genetic causes of disease in two clinically diagnosed patients. In a female patient with OSCS, we identified a mosaic 7-nucleotide frameshift deletion in exon 2 of AMER1, NM_152424.4:c.855_861del:p.(His285Glnfs*7), affecting 8.3% of the WGS reads. In a male patient with CCD, approximately 34% of the WGS reads harbored a 3710-basepair mosaic deletion, NC_000006.11:g.45514471_45518181del, starting in intron 8 of RUNX2 and terminating in the 3' untranslated region. Droplet digital polymerase chain reaction was used to validate these deletions and quantify the absolute level of mosaicism in each patient. Although constitutional variants in AMER1 and RUNX2 are a known cause of OSCS and CCD, respectively, the mosaic changes here reported have not been described previously. Our study indicates that mosaicism should be considered in unsolved cases of skeletal dysplasia and should be investigated with comprehensive and sensitive detection methods. © 2022 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Multi-omics analysis reveals multiple mechanisms causing Prader-Willi like syndrome in a family with a X;15 translocation.

Prader-Willi syndrome (PWS; MIM# 176270) is a neurodevelopmental disorder caused by the loss of expression of paternally imprinted genes within the PWS region located on 15q11.2. It is usually caused by either maternal uniparental disomy of chromosome 15 (UPD15) or 15q11.2 recurrent deletion(s). Here, we report a healthy carrier of a balanced X;15 translocation and her two daughters, both with the karyotype 45,X,der(X)t(X;15)(p22;q11.2),-15. Both daughters display symptoms consistent with haploinsufficiency of the SHOX gene and PWS. We explored the architecture of the derivative chromosomes and investigated effects on gene expression in patient-derived neural cells. First, a multiplex ligation-dependent probe amplification methylation assay was used to determine the methylation status of the PWS-region revealing maternal UPD15 in daughter 2, explaining her clinical symptoms. Next, short read whole genome sequencing and 10X genomics linked read sequencing was used to pinpoint the exact breakpoints of the translocation. Finally, we performed transcriptome sequencing on neuroepithelial stem cells from the mother and from daughter 1 and observed biallelic expression of genes in the PWS region (including SNRPN) in daughter 1. In summary, our multi-omics analysis highlights two different PWS mechanisms in one family and provide an example of how structural variation can affect imprinting through long-range interactions.

Linked-read whole-genome sequencing resolves common and private structural variants in multiple myeloma.

Multiple myeloma (MM) is an incurable and aggressive plasma cell malignancy characterized by a complex karyotype with multiple structural variants (SVs) and copy-number variations (CNVs). Linked-read whole-genome sequencing (lrWGS) allows for refined detection and reconstruction of SVs by providing long-range genetic information from standard short-read sequencing. This makes lrWGS an attractive solution for capturing the full genomic complexity of MM. Here we show that high-quality lrWGS data can be generated from low numbers of cells subjected to fluorescence-activated cell sorting (FACS) without DNA purification. Using this protocol, we analyzed MM cells after FACS from 37 patients with MM using lrWGS. We found high concordance between lrWGS and fluorescence in situ hybridization (FISH) for the detection of recurrent translocations and CNVs. Outside of the regions investigated by FISH, we identified >150 additional SVs and CNVs across the cohort. Analysis of the lrWGS data allowed for resolution of the structure of diverse SVs affecting the MYC and t(11;14) loci, causing the duplication of genes and gene regulatory elements. In addition, we identified private SVs causing the dysregulation of genes recurrently involved in translocations with the IGH locus and show that these can alter the molecular classification of MM. Overall, we conclude that lrWGS allows for the detection of aberrations critical for MM prognostics and provides a feasible route for providing comprehensive genetics. Implementing lrWGS could provide more accurate clinical prognostics, facilitate genomic medicine initiatives, and greatly improve the stratification of patients included in clinical trials.

A somatic UBA2 variant preceded ETV6-RUNX1 in the concordant BCP-ALL of monozygotic twins.

Genetic analysis of leukemic clones in monozygotic twins with concordant acute lymphoblastic leukemia (ALL) has proved a unique opportunity to gain insight into the molecular phylogenetics of leukemogenesis. Using whole-genome sequencing, we characterized constitutional and somatic single nucleotide variants/insertion-deletions (indels) and structural variants in a monozygotic twin pair with concordant ETV6-RUNX1+ B-cell precursor ALL (BCP-ALL). In addition, digital PCR (dPCR) was applied to evaluate the presence of and quantify selected somatic variants at birth, diagnosis, and remission. A shared somatic complex rearrangement involving chromosomes 11, 12, and 21 with identical fusion sequences in leukemias of both twins offered direct proof of a common clonal origin. The ETV6-RUNX1 fusion detected at diagnosis was found to originate from this complex rearrangement. A shared somatic frameshift deletion in UBA2 was also identified in diagnostic samples. In addition, each leukemia independently acquired analogous deletions of 3 genes recurrently targeted in BCP-ALLs (ETV6, ATF7IP, and RAG1/RAG2), providing evidence of a convergent clonal evolution only explained by a strong concurrent selective pressure. Quantification of the UBA2 deletion by dPCR surprisingly indicated it persisted in remission. This, for the first time to our knowledge, provided evidence of a UBA2 variant preceding the well-established initiating event ETV6-RUNX1. Further, we suggest the UBA2 deletion exerted a leukemia predisposing effect and that its essential role in Small Ubiquitin-like Modifier (SUMO) attachment (SUMOylation), regulating nearly all physiological and pathological cellular processes such as DNA-repair by nonhomologous end joining, may hold a mechanistic explanation for the predisposition.

Genomic profile - a possible diagnostic and prognostic marker in upper tract urothelial carcinoma.

OBJECTIVES: To investigate gene alterations as diagnostic and prognostic markers in upper tract urothelial carcinoma (UTUC). PATIENTS AND METHODS: Patients with UTUC who underwent nephroureterectomy between 2005 and 2012 were followed until November 2020. DNA was extracted from paraffin-embedded tumour tissue. Next-generation sequencing using a 388-gene panel was performed. First a blinded analysis using principal component analysis and hierarchical clustering was used to search for patterns of mutations. Then a comparative analysis using analysis of variance (ANOVA) was used to search for mutations enriched in groups of various grades, stages, and survival. In addition, careful manual annotation was used to identify pathogenic mutations over-represented in tumours of high grade/stage and/or poor survival. RESULTS: A total of 39 patients were included. All tumour stages and grades were represented in the cohort. The median follow-up was 10.6 years. In all, 11 patients died from UTUC during the follow-up. Tumour mutational burden showed a statistically significant correlation with stage, grade, and stage + grade. Grade 1, Grade 2, and Grade 3 tumours had different mutational patterns. Patients who died from UTUC had pathogenic mutations in specific genes e.g. tumour protein p53 (TP53) and HRas proto-oncogene, GTPase (HRAS). Patients with Ta Grade 1 tumours with a known pathogenic fibroblast growth factor receptor 3 (FGFR3) mutation did not die from UTUC. CONCLUSION: The genetic analysis was highly concordant with histopathological features and added prognostic information in some cases. Thus, results from genomic profiling may contribute to the choice of treatment and follow-up regimens in the future.

High diagnostic yield in skeletal ciliopathies using massively parallel genome sequencing, structural variant screening and RNA analyses.

Skeletal ciliopathies are a heterogenous group of disorders with overlapping clinical and radiographic features including bone dysplasia and internal abnormalities. To date, pathogenic variants in at least 30 genes, coding for different structural cilia proteins, are reported to cause skeletal ciliopathies. Here, we summarize genetic and phenotypic features of 34 affected individuals from 29 families with skeletal ciliopathies. Molecular diagnostic testing was performed using massively parallel sequencing (MPS) in combination with copy number variant (CNV) analyses and in silico filtering for variants in known skeletal ciliopathy genes. We identified biallelic disease-causing variants in seven genes: DYNC2H1, KIAA0753, WDR19, C2CD3, TTC21B, EVC, and EVC2. Four variants located in non-canonical splice sites of DYNC2H1, EVC, and KIAA0753 led to aberrant splicing that was shown by sequencing of cDNA. Furthermore, CNV analyses showed an intragenic deletion of DYNC2H1 in one individual and a 6.7 Mb de novo deletion on chromosome 1q24q25 in another. In five unsolved cases, MPS was performed in family setting. In one proband we identified a de novo variant in PRKACA and in another we found a homozygous intragenic deletion of IFT74, removing the first coding exon and leading to expression of a shorter message predicted to result in loss of 40 amino acids at the N-terminus. These findings establish IFT74 as a new skeletal ciliopathy gene. In conclusion, combined single nucleotide variant, CNV and cDNA analyses lead to a high yield of genetic diagnoses (90%) in a cohort of patients with skeletal ciliopathies.

Integration of whole genome sequencing into a healthcare setting: high diagnostic rates across multiple clinical entities in 3219 rare disease patients.

BACKGROUND: We report the findings from 4437 individuals (3219 patients and 1218 relatives) who have been analyzed by whole genome sequencing (WGS) at the Genomic Medicine Center Karolinska-Rare Diseases (GMCK-RD) since mid-2015. GMCK-RD represents a long-term collaborative initiative between Karolinska University Hospital and Science for Life Laboratory to establish advanced, genomics-based diagnostics in the Stockholm healthcare setting. METHODS: Our analysis covers detection and interpretation of SNVs, INDELs, uniparental disomy, CNVs, balanced structural variants, and short tandem repeat expansions. Visualization of results for clinical interpretation is carried out in Scout-a custom-developed decision support system. Results from both singleton (84%) and trio/family (16%) analyses are reported. Variant interpretation is done by 15 expert teams at the hospital involving staff from three clinics. For patients with complex phenotypes, data is shared between the teams. RESULTS: Overall, 40% of the patients received a molecular diagnosis ranging from 19 to 54% for specific disease groups. There was heterogeneity regarding causative genes (n = 754) with some of the most common ones being COL2A1 (n = 12; skeletal dysplasia), SCN1A (n = 8; epilepsy), and TNFRSF13B (n = 4; inborn errors of immunity). Some causative variants were recurrent, including previously known founder mutations, some novel mutations, and recurrent de novo mutations. Overall, GMCK-RD has resulted in a large number of patients receiving specific molecular diagnoses. Furthermore, negative cases have been included in research studies that have resulted in the discovery of 17 published, novel disease-causing genes. To facilitate the discovery of new disease genes, GMCK-RD has joined international data sharing initiatives, including ClinVar, UDNI, Beacon, and MatchMaker Exchange. CONCLUSIONS: Clinical WGS at GMCK-RD has provided molecular diagnoses to over 1200 individuals with a broad range of rare diseases. Consolidation and spread of this clinical-academic partnership will enable large-scale national collaboration.

Cell-free tumour DNA analysis detects copy number alterations in gastro-oesophageal cancer patients.

BACKGROUND: Analysis of cell-free tumour DNA, a liquid biopsy, is a promising biomarker for cancer. We have performed a proof-of principle study to test the applicability in the clinical setting, analysing copy number alterations (CNAs) in plasma and tumour tissue from 44 patients with gastro-oesophageal cancer. METHODS: DNA was isolated from blood plasma and a tissue sample from each patient. Array-CGH was applied to the tissue DNA. The cell-free plasma DNA was sequenced by low-coverage whole-genome sequencing using a clinical pipeline for non-invasive prenatal testing. WISECONDOR and ichorCNA, two bioinformatic tools, were used to process the output data and were compared to each other. RESULTS: Cancer-associated CNAs could be seen in 59% (26/44) of the tissue biopsies. In the plasma samples, a targeted approach analysing 61 regions of special interest in gastro-oesophageal cancer detected cancer-associated CNAs with a z-score >5 in 11 patients. Broadening the analysis to a whole-genome view, 17/44 patients (39%) had cancer-associated CNAs using WISECONDOR and 13 (30%) using ichorCNA. Of the 26 patients with tissue-verified cancer-associated CNAs, 14 (54%) had corresponding CNAs in plasma. Potentially clinically actionable amplifications overlapping the genes VEGFA, EGFR and FGFR2 were detected in the plasma from three patients. CONCLUSIONS: We conclude that low-coverage whole-genome sequencing without prior knowledge of the tumour alterations could become a useful tool for cell-free tumour DNA analysis of total CNAs in plasma from patients with gastro-oesophageal cancer.

Novel pathogenic genomic variants leading to autosomal dominant and recessive Robinow syndrome.

Robinow syndrome (RS) is a genetically heterogeneous disorder characterized by skeletal dysplasia and a distinctive facial appearance. Previous studies have revealed locus heterogeneity with rare variants in DVL1, DVL3, FZD2, NXN, ROR2, and WNT5A underlying the etiology of RS. The aforementioned "Robinow-associated genes" and their gene products all play a role in the WNT/planar cell polarity signaling pathway. We performed gene-targeted Sanger sequencing, exome sequencing, genome sequencing, and array comparative genomic hybridization on four subjects with a clinical diagnosis of RS who had not had prior DNA testing. Individuals in our cohort were found to carry pathogenic or likely pathogenic variants in three RS related genes: DVL1, ROR2, and NXN. One subject was found to have a nonsense variant (c.817C > T [p.Gln273*]) in NXN in trans with an ~1 Mb telomeric deletion on chromosome 17p containing NXN, which supports our contention that biallelic NXN variant alleles are responsible for a novel autosomal recessive RS locus. These findings provide increased understanding of the role of WNT signaling in skeletal development and maintenance. These data further support the hypothesis that dysregulation of the noncanonical WNT pathway in humans gives rise to RS.

Sarek: A portable workflow for whole-genome sequencing analysis of germline and somatic variants.

Whole-genome sequencing (WGS) is a fundamental technology for research to advance precision medicine, but the limited availability of portable and user-friendly workflows for WGS analyses poses a major challenge for many research groups and hampers scientific progress. Here we present Sarek, an open-source workflow to detect germline variants and somatic mutations based on sequencing data from WGS, whole-exome sequencing (WES), or gene panels. Sarek features (i) easy installation, (ii) robust portability across different computer environments, (iii) comprehensive documentation, (iv) transparent and easy-to-read code, and (v) extensive quality metrics reporting. Sarek is implemented in the Nextflow workflow language and supports both Docker and Singularity containers as well as Conda environments, making it ideal for easy deployment on any POSIX-compatible computers and cloud compute environments. Sarek follows the GATK best-practice recommendations for read alignment and pre-processing, and includes a wide range of software for the identification and annotation of germline and somatic single-nucleotide variants, insertion and deletion variants, structural variants, tumour sample purity, and variations in ploidy and copy number. Sarek offers easy, efficient, and reproducible WGS analyses, and can readily be used both as a production workflow at sequencing facilities and as a powerful stand-alone tool for individual research groups. The Sarek source code, documentation and installation instructions are freely available at https://github.com/nf-core/sarek and at https://nf-co.re/sarek/.

Rare variants in dynein heavy chain genes in two individuals with situs inversus and developmental dyslexia: a case report.

BACKGROUND: Developmental dyslexia (DD) is a neurodevelopmental learning disorder with high heritability. A number of candidate susceptibility genes have been identified, some of which are linked to the function of the cilium, an organelle regulating left-right asymmetry development in the embryo. Furthermore, it has been suggested that disrupted left-right asymmetry of the brain may play a role in neurodevelopmental disorders such as DD. However, it is unknown whether there is a common genetic cause to DD and laterality defects or ciliopathies. CASE PRESENTATION: Here, we studied two individuals with co-occurring situs inversus (SI) and DD using whole genome sequencing to identify genetic variants of importance for DD and SI. Individual 1 had primary ciliary dyskinesia (PCD), a rare, autosomal recessive disorder with oto-sino-pulmonary phenotype and SI. We identified two rare nonsynonymous variants in the dynein axonemal heavy chain 5 gene (DNAH5): a previously reported variant c.7502G > C; p.(R2501P), and a novel variant c.12043 T > G; p.(Y4015D). Both variants are predicted to be damaging. Ultrastructural analysis of the cilia revealed a lack of outer dynein arms and normal inner dynein arms. MRI of the brain revealed no significant abnormalities. Individual 2 had non-syndromic SI and DD. In individual 2, one rare variant (c.9110A > G;p.(H3037R)) in the dynein axonemal heavy chain 11 gene (DNAH11), coding for another component of the outer dynein arm, was identified. CONCLUSIONS: We identified the likely genetic cause of SI and PCD in one individual, and a possibly significant heterozygosity in the other, both involving dynein genes. Given the present evidence, it is unclear if the identified variants also predispose to DD and further studies into the association between laterality, ciliopathies and DD are needed.

pyCancerSig: subclassifying human cancer with comprehensive single nucleotide, structural and microsatellite mutational signature deconstruction from whole genome sequencing.

BACKGROUND: DNA damage accumulates over the course of cancer development. The often-substantial amount of somatic mutations in cancer poses a challenge to traditional methods to characterize tumors based on driver mutations. However, advances in machine learning technology can take advantage of this substantial amount of data. RESULTS: We developed a command line interface python package, pyCancerSig, to perform sample profiling by integrating single nucleotide variation (SNV), structural variation (SV) and microsatellite instability (MSI) profiles into a unified profile. It also provides a command to decipher underlying cancer processes, employing an unsupervised learning technique, Non-negative Matrix Factorization, and a command to visualize the results. The package accepts common standard file formats (vcf, bam). The program was evaluated using a cohort of breast- and colorectal cancer from The Cancer Genome Atlas project (TCGA). The result showed that by integrating multiple mutations modes, the tool can correctly identify cases with known clear mutational signatures and can strengthen signatures in cases with unclear signal from an SNV-only profile. The software package is available at https://github.com/jessada/pyCancerSig. CONCLUSIONS: pyCancerSig has demonstrated its capability in identifying known and unknown cancer processes, and at the same time, illuminates the association within and between the mutation modes.

Alu-Alu mediated intragenic duplications in IFT81 and MATN3 are associated with skeletal dysplasias.

Skeletal dysplasias are a diverse group of rare Mendelian disorders with clinical and genetic heterogeneity. Here, we used targeted copy number variant (CNV) screening and identified intragenic exonic duplications, formed through Alu-Alu fusion events, in two individuals with skeletal dysplasia and negative exome sequencing results. First, we detected a homozygous tandem duplication of exon 9 and 10 in IFT81 in a boy with Jeune syndrome, or short-rib thoracic dysplasia (SRTD) (MIM# 208500). Western blot analysis did not detect any wild-type IFT81 protein in fibroblasts from the patient with the IFT81 duplication, but only a shorter isoform of IFT81 that was also present in the normal control samples. Complementary zebrafish studies suggested that loss of full-length IFT81 protein but expression of a shorter form of IFT81 protein affects the phenotype while being compatible with life. Second, a de novo tandem duplication of exons 2 to 5 in MATN3 was identified in a girl with multiple epiphyseal dysplasia (MED) type 5 (MIM# 607078). Our data highlights the importance of detection and careful characterization of intragenic duplication CNVs, presenting them as a novel and very rare genetic mechanism in IFT81-related Jeune syndrome and MATN3-related MED.

Targeted copy number screening highlights an intragenic deletion of WDR63 as the likely cause of human occipital encephalocele and abnormal CNS development in zebrafish.

Congenital malformations affecting the neural tube can present as isolated malformations or occur in association with other developmental abnormalities and syndromes. Using high-resolution copy number screening in 66 fetuses with neural tube defects, we identified six fetuses with likely pathogenic mutations, three aneuploidies (one trisomy 13 and two trisomy 18) and three deletions previously reported in NTDs (one 22q11.2 deletion and two 1p36 deletions) corresponding to 9% of the cohort. In addition, we identified five rare deletions and two duplications of uncertain significance including a rare intragenic heterozygous in-frame WDR63 deletion in a fetus with occipital encephalocele. Whole genome sequencing verified the deletion and excluded known pathogenic variants. The deletion spans exons 14-17 resulting in the expression of a protein missing the third and fourth WD-repeat domains. These findings were supported by CRISPR/Cas9-mediated somatic deletions in zebrafish. Injection of two different sgRNA-pairs targeting relevant intronic regions resulted in a deletion mimicking the human deletion and a concomitant increase of abnormal embryos with body and brain malformations (41%, n = 161 and 62%, n = 224, respectively), including a sac-like brain protrusion (7% and 9%, P < 0.01). Similar results were seen with overexpression of RNA encoding the deleted variant in zebrafish (total abnormal; 46%, n = 255, P < 0.001) compared with the overexpression of an equivalent amount of wild-type RNA (total abnormal; 3%, n = 177). We predict the in-frame WDR63 deletion to result in a dominant negative or gain-of-function form of WDR63. These are the first findings supporting a role for WDR63 in encephalocele formation.