Search

Abstract Inversions are balanced structural variants that often remain undetected in genetic diagnostics. We present a female proband with a de novo Chromosome 15 paracentric inversion, disrupting MEIS2 and NUSAP1. The inversion was detected by short-read genome sequencing and confirmed with adaptive long-read sequencing. The breakpoint junction analysis revealed a 96 bp (bp) deletion and an 18 bp insertion in the two junctions, suggesting that the rearrangement arose through a replicative error. Transcriptome sequencing of cultured fibroblasts revealed normal MEIS2 levels and 0.61-fold decreased expression of NUSAP1. Furthermore, three fusion transcripts were detected and confirmed by Sanger sequencing. Heterozygous loss of MEIS2 (MIM# 600987) is associated with a cleft palate, heart malformations, and intellectual impairment, which overlap with the clinical symptoms observed in the proband. The observed fusion transcripts are likely non-functional, and MEIS2 haploinsufficiency is the likely disease causative mechanism. Altogether, this study’s findings illustrate the importance of including inversions in rare disease diagnostic testing and highlight the value of long read sequencing for the validation and characterization of such variants.

Abstract Long-read genome sequencing (lrGS) is a promising method in genetic diagnostics. Here we investigate the potential of lrGS to detect a disease-associated chromosomal translocation between 17p13 and the 19 centromere. We constructed two sets of phased and non-phased de novo assemblies; (i) based on lrGS only and (ii) hybrid assemblies combining lrGS with optical mapping using lrGS reads with a median coverage of 34X. Variant calling detected both structural variants (SVs) and small variants and the accuracy of the small variant calling was compared with those called with short-read genome sequencing (srGS). The de novo and hybrid assemblies had high quality and contiguity with N50 of 62.85 Mb, enabling a near telomere to telomere assembly with less than a 100 contigs per haplotype. Notably, we successfully identified the centromeric breakpoint of the translocation. A concordance of 92% was observed when comparing small variant calling between srGS and lrGS. In summary, our findings underscore the remarkable potential of lrGS as a comprehensive and accurate solution for the analysis of SVs and small variants. Thus, lrGS could replace a large battery of genetic tests that were used for the diagnosis of a single symptomatic translocation carrier, highlighting the potential of lrGS in the realm of digital karyotyping.

Loss-of-function variants in the FLG gene have been identified as the strongest cause of susceptibility to atopic dermatitis (AD) in Europeans and Asians. However, very little is known about the genetic etiology behind AD in African populations, where the prevalence of AD is notably high. We sought to investigate the genetic origins of AD by performing whole-genome sequencing in an Ethiopian family with 12 individuals and several affected in different generations. We identified 2 variants within FLG2 (p.D13Y) and NOD2 (p.A918S) genes cosegregating with AD in the affected individuals. Further genotyping analyses in both Ethiopian and Swedish AD cases and controls revealed a significant association with the FLG2 variant (p.D13Y, P < .0013) only in the Ethiopian cohort. However, the NOD2 variant (p.A918S) did not show any association in our Ethiopian cohort. Instead, 2 previously recognized NOD2 variants (p.A849V, P < .0085 and p.G908R, P < .0036) were significantly associated with AD in our Ethiopian cohort. Our study indicates that the FLG2 and NOD2 genes might be important in the etiology of AD in Ethiopians. Additional genetic and functional studies are needed to confirm the role of these genes and the associated variants into the development of AD.

Abstract Psoriatic arthritis mutilans (PAM) is the rarest and most severe form of psoriatic arthritis, characterized by erosions of the small joints and osteolysis leading to joint disruption. Despite its severity, the underlying mechanisms are unknown, and no susceptibility genes have hitherto been identified. We aimed to investigate the genetic basis of PAM by performing massive parallel sequencing in sixty-one patients from the PAM Nordic cohort. We found rare variants in the NADPH oxidase 4 (NOX4) in four patients. In silico predictions show that the identified variants are potentially damaging. NOXs are the only enzymes producing reactive oxygen species (ROS). NOX4 is specifically involved in the differentiation of osteoclasts, the cells implicated in bone resorption. Functional follow-up studies using cell culture, zebrafish models, and measurement of ROS in patients uncovered that these NOX4 variants increase ROS levels both in vitro and in vivo. We propose NOX4 as the first candidate susceptibility gene for PAM. Our study links high levels of ROS caused by NOX4 variants to the development of PAM, offering a potential therapeutic target.

Abstract 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.

Abstract Long-read sequencing has dramatically increased our understanding of human genome variation. Here, we demonstrate that long-read technology can give new insights into the genomic architecture of individual cells. Clonally expanded CD8+ T-cells from a human donor were subjected to droplet-based multiple displacement amplification (dMDA) to generate long molecules with reduced bias. PacBio sequencing generated up to 40% genome coverage per single-cell, enabling detection of single nucleotide variants (SNVs), structural variants (SVs), and tandem repeats, also in regions inaccessible by short reads. 28 somatic SNVs were detected, including one case of mitochondrial heteroplasmy. 5473 high-confidence SVs/cell were discovered, a sixteen-fold increase compared to Illumina-based results from clonally related cells. Single-cell de novo assembly generated a genome size of up to 598 Mb and 1762 (12.8%) complete gene models. In summary, our work shows the promise of long-read sequencing toward characterization of the full spectrum of genetic variation in single cells.

Purpose: Although chromosome 21 is the smallest human chromosome, it is highly relevant in the pathogenicity of both cancer and congenital diseases, including Alzheimer disease and trisomy 21 (Down syndrome). In addition, cases with rare structural variants (SVs) of chromosome 21 have been reported. These events vary in size and include large chromosomal events, such as ring chromosomes and small partial aneuploidies. The p-arm of the acrocentric chromosome 21 was devoid of reference genomic sequence in GRCh37 and GRCh38, which hampered our ability to solve genomic rearrangements and find the mechanism of formation of disease-causing SVs. We hypothesize that conserved satellite structures and segmental duplications located on the p-arm play an important role in the formation of complex SVs involving chromosome 21. Methods: Three cases with complex chromosome 21 rearrangements were studied with a combination of short-read and long-read genome sequencing, as well as optical genome mapping. The data were aligned to the T2T-CHM13 assembly. Results: We were able to resolve all 3 complex chromosome 21 rearrangements in which 15, 8, and 26 breakpoints were identified, respectively. By comparing the identified SV breakpoints, we were able to pinpoint a region between 21p13 and 21p12 that appears to be frequently involved in chromosome 21 rearrangements. Importantly, we observed acrocentric satellite DNA at several breakpoint junctions suggesting an important role for those elements in the formation of complex SVs. Conclusion: Taken together, our results provide further insights into the architecture and underlying mechanisms of complex rearrangements on acrocentric chromosomes.

IntroductionThe 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.MethodsFor 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.ResultsBoth 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.DiscussionThe 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.

FOXC1 is a ubiquitously expressed forkhead transcription factor that plays a critical role during early development. Germline pathogenic variants in FOXC1 are associated with anterior segment dysgenesis and Axenfeld-Rieger syndrome (ARS, #602482), an autosomal dominant condition with ophthalmologic anterior segment abnormalities, high risk for glaucoma and extraocular findings including distinctive facial features, as well as dental, skeletal, audiologic, and cardiac anomalies. De Hauwere syndrome is an ultrarare condition previously associated with 6p microdeletions and characterized by anterior segment dysgenesis, joint instability, short stature, hydrocephalus, and skeletal abnormalities. Here, we report clinical findings of two unrelated adult females with FOXC1 haploinsufficiency who have ARS and skeletal abnormalities. Final molecular diagnoses of both patients were achieved using genome sequencing. Patient 1 had a complex rearrangement involving a 4.9 kB deletion including FOXC1 coding region (Hg19; chr6:1,609,721-1,614,709), as well as a 7 MB inversion (Hg19; chr6:1,614,710-8,676,899) and a second deletion of 7.1 kb (Hg19; chr6:8,676,900-8,684,071). Patient 2 had a heterozygous single nucleotide deletion, resulting in a frameshift and a premature stop codon in FOXC1 (NM_001453.3): c.467del, p.(Pro156Argfs*25). Both individuals had moderate short stature, skeletal abnormalities, anterior segment dysgenesis, glaucoma, joint laxity, pes planovalgus, dental anomalies, hydrocephalus, distinctive facial features, and normal intelligence. Skeletal surveys revealed dolichospondyly, epiphyseal hypoplasia of femoral and humeral heads, dolichocephaly with frontal bossin gand gracile long bones. We conclude that haploinsufficiency of FOXC1 causes ARS and a broad spectrum of symptoms with variable expressivity that at its most severe end also includes a phenotype overlapping with De Hauwere syndrome.

The vast majority of the human genome is non-coding. There is a diversity of non-coding features, some of which have functional importance. Although the non-coding regions constitute the majority of the genome, they remain understudied, and for a long time, these regions have been referred to as junk DNA. Pseudogenes are one of these features. A pseudogene is a non-functional copy of a protein-coding gene. Pseudogenes may arise through a variety of genetic mechanisms. Processed pseudogenes are formed through reverse transcription of mRNA by LINE elements, after which the cDNA is integrated into the genome. Processed pseudogenes are known to be variable across populations; however, the variability and distribution remains unknown. Herein, we apply a custom-designed processed pseudogene pipeline on the whole genome sequencing data of 3,500 individuals; 2,500 individuals from the thousand genomes dataset, as well as 1,000 Swedish individuals. Through these analyses, we discover over 3,000 pseudogenes missing from the GRCh38 reference. Utilising our pipeline, we position 74% of the detected processed pseudogenes—allowing for analyses of formation. Notably, we find that common structural variant callers, such as Delly, classify the processed pseudogenes as deletion events, which are later predicted to be truncating variants. By compiling lists of non-reference processed pseudogenes and their frequencies, we find a great variability of pseudogenes; indicating that non-reference processed pseudogenes may be useful for DNA testing and as population-specific markers. In summary, our findings highlight a great diversity of processed pseudogenes, that processed pseudogenes are actively formed in the human genome; and that our pipeline may be used to reduce false positive structural variation caused by the misalignment and subsequent misclassification of non-reference processed pseudogenes.

Abstract Urinary bladder wound healing relies on multiple biological events that are finely tuned in a spatial–temporal manner. MicroRNAs are small non-coding RNA molecules with regulatory functions. We hypothesized that microRNAs are important molecules in the coordination of normal urinary bladder wound healing. We aimed at identifying microRNAs expressed during bladder wound healing using Affymetrix global array for microRNA profiling of the rodent urinary bladder during healing of a surgically created wound. Results were validated in the rat bladders by real-time PCR (RT-PCR) using three of the differentially expressed (DE) microRNAs. The model was thereafter validated in human cells, by measuring the expression of eight of the DE microRNAs upon in vitro wound-healing assays in primary urothelial cells. Our results indicated that 508 (40%) of all rodent microRNAs were expressed in the urinary bladder during wound healing. Thirteen of these microRNAs (1%) were DE (false discovery rate (FDR) 0.25) in wounded compared to non-wounded bladders. Bioinformatic analyses helped us to identify target molecules for the DE microRNAs, and biological pathways involved in tissue repair. All data are made available in an open-access database for other researchers to explore.

Abstract We identified an autosomal dominant progranulin mutation carrier without symptoms of dementia in her lifetime (Reduced Penetrance Mutation Carrier, RedPenMC). This resistance to develop expected pathology presents a unique opportunity to interrogate neurodegenerative mechanisms. We performed multimodal single-nuclei analyses of post-mortem frontal cortex from RedPenMC, including transcriptomics and global levels of chromatin marks. RedPenMC had an increased ratio of GRN-expressing microglia, higher levels of activating histone mark H3k4me3 in microglia and lower levels of the repressive chromatin marks H3k9me1 and H3k9me3 in the frontal cortex than her affected mutation carrier son and evidence of higher protein levels of progranulin in both plasma and brain homogenates. Although the study is limited to one case, the results support that restoring brain progranulin levels may be sufficient to escape neurodegeneration and FTD. In addition to previously identified modifier genes, it is possible that epigenetic marks may contribute to the increased progranulin expression in cases of reduced penetrance. These findings may stimulate similar follow-up studies and new therapeutic approaches.

ABSTRACT 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.

Abstract 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.

Abstract Background Exome and genome sequencing is becoming the method of choice for rare disease diagnostics. One of the key challenges remaining is distinguishing the disease causing variants from the benign background variation. After analysis and annotation of the sequencing data there are typically thousands of candidate variants requiring further investigation. One of the most effective and least biased ways to reduce this number is to assess the rarity of a variant in any population. Currently, there are a number of reliable sources of information for major population frequencies when considering single nucleotide variants (SNVs) and small insertion and deletions (INDELs), with gnomAD as the most prominent public resource available. However, local variation or frequencies in sub-populations may be underrepresented in these public resources. In contrast, for structural variation (SV), the background frequency in the general population is more or less unknown mostly due to challenges in calling SVs in a consistent way. Keeping track of local variation is one way to overcome these problems and significantly reduce the number of potential disease causing variants retained for manual inspection, both for SNVs and SVs. Results Here, we present loqusdb, a tool to solve the challenge of keeping track of any type of variant observations from genome sequencing data. Loqusdb was designed to handle a large flow of samples and unlike other solutions, samples can be added continuously to the database without rebuilding it, facilitating improvements and additions. We assessed the added value of a local observations database using 98 samples annotated with information from a background of 888 unrelated individuals. Conclusions We show both how powerful SV analysis can be when filtering for population frequencies and how the number of apparently rare SNVs/INDELs can be reduced by adding local population information even after annotating the data with other large frequency databases, such as gnomAD. In conclusion, we show that a local frequency database is an attractive, and a necessary addition to the publicly available databases that facilitate the analysis of exome and genome data in a clinical setting.

Abstract 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.

Abstract 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.

Induced pluripotent stem cells (iPSCs) from patients are an attractive disease model to study tissues with poor accessibility such as the brain. Using this approach, we and others have shown that trisomy 21 results in genome-wide transcriptional dysregulations. The effects of loss of genes on chromosome 21 is much less characterized. Here, we use patient-derived neural cells from an individual with neurodevelopmental delay and a ring chromosome 21 with two deletions spanning 3.8 Mb at the terminal end of 21q22.3, containing 60 protein-coding genes. To investigate the molecular perturbations of the partial monosomy on neural cells, we established patient-derived iPSCs from fibroblasts retaining the ring chromosome 21, and we then induced iPSCs into neuroepithelial stem cells. RNA-Seq analysis of NESCs with the ring chromosome revealed downregulation of 18 genes within the deleted region together with global transcriptomic dysregulations when compared to euploid NESCs. Since the deletions on chromosome 21 represent a genetic “contrary” to trisomy of the corresponding region, we further compared the dysregulated transcriptomic profile in with that of two NESC lines with trisomy 21. The analysis revealed opposed expression changes for 23 genes on chromosome 21 as well as 149 non-chromosome 21 genes. Taken together, our results bring insights into the effects on the global and chromosome 21 specific gene expression from a partial monosomy of chromosome 21qter during early neuronal differentiation.

Chromoanagenesis is a descriptive term that encompasses classes of catastrophic mutagenic processes that generate localized and complex chromosome rearrangements in both somatic and germline genomes. Herein, we describe a 5-year-old female presenting with a constellation of clinical features consistent with a clinical diagnosis of Coffin–Siris syndrome 1 (CSS1). Initial G-banded karyotyping detected a 90-Mb pericentric and a 47-Mb paracentric inversion on a single chromosome. Subsequent analysis of short-read whole-genome sequencing data and genomic optical mapping revealed additional inversions, all clustered on chromosome 6, one of them disrupting ARID1B for which haploinsufficiency leads to the CSS1 disease trait (MIM:135900). The aggregate structural variant data show that the resolved, the resolved derivative chromosome architecture presents four de novo inversions, one pericentric and three paracentric, involving six breakpoint junctions in what appears to be a shuffling of genomic material on this chromosome. Each junction was resolved to nucleotide-level resolution with mutational signatures suggestive of non-homologous end joining. The disruption of the gene ARID1B is shown to occur between the fourth and fifth exon of the canonical transcript with subsequent qPCR studies confirming a decrease in ARID1B expression in the patient versus healthy controls. Deciphering the underlying genomic architecture of chromosomal rearrangements and complex structural variants may require multiple technologies and can be critical to elucidating the molecular etiology of a patient’s clinical phenotype or resolving unsolved Mendelian disease cases.

Abstract Background Since different types of genetic variants, from single nucleotide variants (SNVs) to large chromosomal rearrangements, underlie intellectual disability, we evaluated the use of whole-genome sequencing (WGS) rather than chromosomal microarray analysis (CMA) as a first-line genetic diagnostic test. Methods We analyzed three cohorts with short-read WGS: (i) a retrospective cohort with validated copy number variants (CNVs) (cohort 1, n = 68), (ii) individuals referred for monogenic multi-gene panels (cohort 2, n = 156), and (iii) 100 prospective, consecutive cases referred to our center for CMA (cohort 3). Bioinformatic tools developed include FindSV, SVDB, Rhocall, Rhoviz, and vcf2cytosure. Results First, we validated our structural variant (SV)-calling pipeline on cohort 1, consisting of three trisomies and 79 deletions and duplications with a median size of 850 kb (min 500 bp, max 155 Mb). All variants were detected. Second, we utilized the same pipeline in cohort 2 and analyzed with monogenic WGS panels, increasing the diagnostic yield to 8%. Next, cohort 3 was analyzed by both CMA and WGS. The WGS data was processed for large (> 10 kb) SVs genome-wide and for exonic SVs and SNVs in a panel of 887 genes linked to intellectual disability as well as genes matched to patient-specific Human Phenotype Ontology (HPO) phenotypes. This yielded a total of 25 pathogenic variants (SNVs or SVs), of which 12 were detected by CMA as well. We also applied short tandem repeat (STR) expansion detection and discovered one pathologic expansion in ATXN7. Finally, a case of Prader-Willi syndrome with uniparental disomy (UPD) was validated in the WGS data. Important positional information was obtained in all cohorts. Remarkably, 7% of the analyzed cases harbored complex structural variants, as exemplified by a ring chromosome and two duplications found to be an insertional translocation and part of a cryptic unbalanced translocation, respectively. Conclusion The overall diagnostic rate of 27% was more than doubled compared to clinical microarray (12%). Using WGS, we detected a wide range of SVs with high accuracy. Since the WGS data also allowed for analysis of SNVs, UPD, and STRs, it represents a powerful comprehensive genetic test in a clinical diagnostic laboratory setting.

BackgroundAnalysis 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.MethodsDNA 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.ResultsCancer-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.ConclusionsWe 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.

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/.

Optic nerve hypoplasia (ONH) is a congenital malformation with a reduced number of retinal ganglion cell axons in a thin optic nerve. It is a common cause of visual impairment in children and ONH is associated with neurodevelopmental disorders, pituitary hormone deficiencies, and brain malformations. In most cases, the aetiology is unknown, but both environmental factors and genetic causes have been described. This study aimed to identify genetic variants underlying ONH in a well-characterised cohort of individuals with ONH. We performed array comparative genomic hybridization and whole genome sequencing in 29 individuals with ONH. Rare variants were verified by Sanger sequencing and inheritance was assessed in parental samples. We identified 11 rare single nucleotide variants (SNVs) in ten individuals, including a homozygous variant in KIF7 (previously associated with Joubert syndrome), a heterozygous de novo variant in COL4A1 (previously described in an individual with porencephaly), and a homozygous variant in COL4A2. In addition, one individual harboured a heterozygous variant in OPA1 and a heterozygous variant in COL4A1, both were inherited and assessed as variants of unknown clinical significance. Finally, a heterozygous deletion of 341 kb involving exons 7-18 of SOX5 (associated with Lamb-Schaffer syndrome) was identified in one individual. The overall diagnostic yield of pathogenic or likely pathogenic variants in individuals with ONH using whole genome sequencing was 4/29 (14%). Our results show that there is a genetic heterogeneity in ONH and indicate that genetic causes of ONH are not rare. We conclude that genetic testing is valuable in a substantial proportion of the individuals with ONH, especially in cases with non-isolated ONH.

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/.

Complex chromosomal rearrangements (CCRs) are rearrangements involving more than two chromosomes or more than two breakpoints. Whole genome sequencing (WGS) allows for outstanding high resolution characterization on the nucleotide level in unique sequences of such rearrangements, but problems remain for mapping breakpoints in repetitive regions of the genome, which are known to be prone to rearrangements. Hence, multiple complementary WGS experiments are sometimes needed to solve the structures of CCRs. We have studied three individuals with CCRs: Case 1 and Case 2 presented with de novo karyotypically balanced, complex interchromosomal rearrangements (46,XX,t(2;8;15)(q35;q24.1;q22) and 46,XY,t(1;10;5)(q32;p12;q31)), and Case 3 presented with a de novo, extremely complex intrachromosomal rearrangement on chromosome 1. Molecular cytogenetic investigation revealed cryptic deletions in the breakpoints of chromosome 2 and 8 in Case 1, and on chromosome 10 in Case 2, explaining their clinical symptoms. In Case 3, 26 breakpoints were identified using WGS, disrupting five known disease genes. All rearrangements were subsequently analyzed using optical maps, linked-read WGS, and short-read WGS. In conclusion, we present a case series of three unique de novo CCRs where we by combining the results from the different technologies fully solved the structure of each rearrangement. The power in combining short-read WGS with long-molecule sequencing or optical mapping in these unique de novo CCRs in a clinical setting is demonstrated.

Clustered copy number variants (CNVs) as detected by chromosomal microarray analysis (CMA) are often reported as germline chromothripsis. However, such cases might need further investigations by massive parallel whole genome sequencing (WGS) in order to accurately define the underlying complex rearrangement, predict the occurrence mechanisms and identify additional complexities. Here, we utilized WGS to delineate the rearrangement structure of 21 clustered CNV carriers first investigated by CMA and identified a total of 83 breakpoint junctions (BPJs). The rearrangements were further sub-classified depending on the patterns observed: I) Cases with only deletions (n = 8) often had additional structural rearrangements, such as insertions and inversions typical to chromothripsis; II) cases with only duplications (n = 7) or III) combinations of deletions and duplications (n = 6) demonstrated mostly interspersed duplications and BPJs enriched with microhomology. In two cases the rearrangement mutational signatures indicated both a breakage-fusion-bridge cycle process and haltered formation of a ring chromosome. Finally, we observed two cases with Alu- and LINE-mediated rearrangements as well as two unrelated individuals with seemingly identical clustered CNVs on 2p25.3, possibly a rare European founder rearrangement. In conclusion, through detailed characterization of the derivative chromosomes we show that multiple mechanisms are likely involved in the formation of clustered CNVs and add further evidence for chromoanagenesis mechanisms in both "simple" and highly complex chromosomal rearrangements. Finally, WGS characterization adds positional information, important for a correct clinical interpretation and deciphering mechanisms involved in the formation of these rearrangements.

Copy number variations (CNVs) within the human genome have been linked to a diversity of inherited diseases and phenotypic traits. The currently used methodology to measure copy numbers has limited resolution and/or precision, especially for regions with more than 4 copies. Whole genome sequencing (WGS) offers an alternative data source to allow for the detection and characterization of the copy number across different genomic regions in a single experiment. A plethora of tools have been developed to utilize WGS data for CNV detection. None of these tools are designed specifically to accurately estimate copy numbers of complex regions in a small cohort or clinical setting. Herein, we present AMYCNE (automatic modeling functionality for copy number estimation), a CNV analysis tool using WGS data. AMYCNE is multifunctional and performs copy number estimation of complex regions, annotation of VCF files, and CNV detection on individual samples. The performance of AMYCNE was evaluated using AMY1A ddPCR measurements from 86 unrelated individuals. In addition, we validated the accuracy of AMYCNE copy number predictions on two additional genes (FCGR3A and FCGR3B) using datasets available through the 1000 genomes consortium. Finally, we simulated levels of mosaic loss and gain of chromosome X and used this dataset for benchmarking AMYCNE. The results show a high concordance between AMYCNE and ddPCR, validating the use of AMYCNE to measure tandem AMY1 repeats with high accuracy. This opens up new possibilities for the use of WGS for accurate copy number determination of other complex regions in the genome in small cohorts or single individuals.

Reliable detection of large structural variation ( > 1000 bp) is important in both rare and common genetic disorders. Whole genome sequencing (WGS) is a technology that may be used to identify a large proportion of the genomic structural variants (SVs) in an individual in a single experiment. Even though SV callers have been extensively used in research to detect mutations, the potential usage of SV callers within routine clinical diagnostics is still limited. One well known, but not well-addressed problem is the large number of benign variants and reference errors present in the human genome that further complicates analysis. Even though there is a wide range of SV-callers available, the number of callers that allow detection of the entire spectra of SV at a low computational cost is still relatively limited.

Reliable detection of large structural variation ( > 1000 bp) is important in both rare and common genetic disorders. Whole genome sequencing (WGS) is a technology that may be used to identify a large proportion of the genomic structural variants (SVs) in an individual in a single experiment. Even though SV callers have been extensively used in research to detect mutations, the potential usage of SV callers within routine clinical diagnostics is hindered by high computational costs, usage of non-standard output format, and limited support for the various sequencing platforms and libraries. Another well known, but not well-addressed problem is the large number of benign variants and reference errors present in the human genome that further complicates analysis. Here we present TIDDIT, a time efficient variant caller, that uses discordant read pairs as well as the depth of coverage and split reads to detect and classify a large spectrum of SVs. As part of the software suite, TIDDIT also includes a database functionality that enables filtering for rare variants and reduces the number of false positive calls and background noise. Benchmarked against five state-of-the-art SV callers, TIDDIT performs at an equal/superior level while using only 2 CPU hours per sample. Thanks to its speed, sensitivity, flexibility and ability to easily detect variants on a wide range of WGS library types, TIDDIT solves many of the problems that are currently hindering the utilization of WGS for SV calling in clinical settings.