TIDDIT, an efficient and comprehensive structural variant caller for massive parallel sequencing data

Understanding the epidemiology and pathogenesis of Mycobacterium tuberculosis with non-redundant pangenome of epidemic strains in China

Tuberculosis is a major public health threat resulting in more than one million lives lost every year. Many challenges exist to defeat this deadly infectious disease which address the importance of a thorough understanding of the biology of the causative agent Mycobacterium tuberculosis (MTB). We generated a non-redundant pangenome of 420 epidemic MTB strains from China including 344 Lineage 2 strains, 69 Lineage 4 strains, six Lineage 3 strains, and one Lineage 1 strain. We estimate that MTB strains have a pangenome of 4,278 genes encoding 4,183 proteins, of which 3,438 are core genes. However, due to 99,694 interruptions in 2,447 coding genes, we can only confidently confirm 1,651 of these genes are translated in all samples. Of these interruptions, 67,315 (67.52%) could be classified by various genetic variations detected by currently available tools, and more than half of them are due to structural variations, mostly small indels. Assuming a proportion of these interruptions are artifacts, the number of active core genes would still be much lower than 3,438. We further described differential evolutionary patterns of genes under the influences of selective pressure, population structure and purifying selection. While selective pressure is ubiquitous among these coding genes, evolutionary adaptations are concentrated in 1,310 genes. Genes involved in cell wall biogenesis are under the strongest selective pressure, while the biological process of disruption of host organelles indicates the direction of the most intensive positive selection. This study provides a comprehensive view on the genetic diversity and evolutionary patterns of coding genes in MTB which may deepen our understanding of its epidemiology and pathogenicity.

Single-Cell RNA Sequencing of Ewing Sarcoma Tumors Demonstrates Transcriptional Heterogeneity and Clonal Evolution

PURPOSE

Ewing sarcoma is the second most common bone cancer in children, accounting for 2% of pediatric cancer diagnoses. Patients who present with metastatic disease at the time of diagnosis have a dismal prognosis compared with the >70% 5-year survival of those with localized disease. Novel therapeutic approaches that can impact metastatic disease are desperately needed, as well as a deeper understanding of the heterogeneity of Ewing sarcoma tumors.

EXPERIMENTAL DESIGN

In this study, we utilized single-cell RNA sequencing to characterize the transcriptional landscape of primary Ewing sarcoma tumors and the surrounding tumor microenvironment in a cohort of seven untreated patients with Ewing sarcoma, as well as in circulating tumor cells (CTC). A potential CTC therapeutic target was evaluated through immunofluorescence of fixed CTCs from a separate cohort.

RESULTS

Primary tumor samples demonstrate a heterogeneous transcriptional landscape with several conserved gene expression programs, including those composed of genes related to proliferation and Ewing sarcoma gene targets, which were found to correlate with overall survival. Copy-number analysis identified subclonal evolution within patients prior to treatment. Analyses of the immune microenvironment reveal an immunosuppressive microenvironment with complex intercellular communication among the tumor and immune cells. Single-cell RNA sequencing and immunofluorescence of CTCs at the time of diagnosis identified TSPAN8 as a potential therapeutic target.

CONCLUSIONS

Ewing sarcoma tumors demonstrate significant transcriptional heterogeneity as well as a complex immunosuppressive microenvironment. This work evaluates several proposed targets that warrant further exploration as novel therapeutic strategies.

Transcription-replication conflicts drive R-loop-dependent nucleosome eviction and require DOT1L activity for transcription recovery

Progressing transcription and replication machineries profoundly impact their underlying chromatin template. Consequently, transcription-replication conflict (TRC) sites are vulnerable to chromatin and epigenome alterations, provoking genome instability. Here, we engineered an inducible TRC reporter system using a genome-integrated R-loop-prone sequence and characterized the dynamic changes of the local chromatin structure inflicted by TRCs, leading to reduced nucleosome occupancy and replication fork blockage. Strikingly, inducing a small number of TRCs on the genome results in a measurable global replication stress response. Furthermore, we find a TRC-dependent increase in H3K79 methylation specifically at the R-loop forming TRC site. Accordingly, inhibition of the H3K79 methyltransferase DOT1L leads to reduced transcriptional output and an exacerbated DNA damage response, suggesting that deposition of this mark is required for effective transcription recovery and resolution of TRCs. Our work shows the molecular dynamics and reveals a specific epigenetic modifier bookmarking TRC sites, relevant to cancer and other diseases.

A de novo, mosaic and complex chromosome 21 rearrangement causes APP triplication and familial autosomal dominant early onset Alzheimer disease

Copy number variation (CNV) of the amyloid-β precursor protein gene (APP) is a known cause of autosomal dominant Alzheimer disease (ADAD), but de novo genetic variants causing ADAD are rare. We report a mother and daughter with neuropathologically confirmed definite Alzheimer disease (AD) and extensive cerebral amyloid angiopathy (CAA). Copy number analysis identified an increased number of APP copies and genome sequencing (GS) revealed the underlying complex genomic rearrangement (CGR) including a triplication of APP with two unique breakpoint junctions (BPJs). The mosaic state in the mother had likely occurred de novo. Digital droplet PCR (ddPCR) on 42 different tissues, including 17 different brain regions, showed the derivative chromosome at varying mosaic levels (20-96%) in the mother who had symptom onset at age 58 years. In contrast, the derivative chromosome was present in all analyzed cells in the daughter whose symptom onset was at 34 years. This study reveals the architecture of a de novo CGR causing APP triplication and ADAD with a striking difference in age at onset between the fully heterozygous daughter compared to the mosaic mother. The GS analysis identified the complexity of the CGR illustrating its usefulness in identifying structural variants (SVs) in neurodegenerative disorders.

Revealing genes related teat number traits via genetic variation in Yorkshire pigs based on whole-genome sequencing

BACKGROUND

Teat number is one of the most important indicators to evaluate the lactation performance of sows, and increasing the teat number has become an important method to improve the economic efficiency of farms. Therefore, it is particularly important to deeply analyze the genetic mechanism of teat number traits in pigs. In this study, we detected Single Nucleotide Ploymorphism (SNP), Insertion-Deletion (InDel) and Structural variant (SV) by high-coverage whole-genome resequencing data, and selected teat number at birth and functional teat number as two types of teat number traits for genome-wide association study (GWAS) to reveal candidate genes associated with pig teat number traits.

RESULTS

In this study, we used whole genome resequencing data from 560 Yorkshire sows to detect SNPs, InDels and SVs, and performed GWAS for the traits of born teat number and functional teat number, and detected a total of 85 significant variants and screened 214 candidate genes, including HEG1, XYLT1, SULF1, MUC13, VRTN, RAP1A and NPVF. Among them, HEG1 and XYLT1 were the new candidate genes in this study. The co-screening and population validation of multiple traits suggested that HEG1 may have a critical effect on the born teat number.

CONCLUSION

Our study shows that more candidate genes associated with pig teat number traits can be identified by GWAS with different variant types. Through large population validation, we found that HEG1 may be a new key candidate gene affecting pig teat number traits. In conclusion, the results of this study provide new information for exploring the genetic mechanisms affecting pig teat number traits and genetic improvement of pigs.

A national long-read sequencing study on chromosomal rearrangements uncovers hidden complexities

Clinical genetic laboratories often require a comprehensive analysis of chromosomal rearrangements/structural variants (SVs), from large events like translocations and inversions to supernumerary ring/marker chromosomes and small deletions or duplications. Understanding the complexity of these events and their clinical consequences requires pinpointing breakpoint junctions and resolving the derivative chromosome structure. This task often surpasses the capabilities of short-read sequencing technologies. In contrast, long-read sequencing techniques present a compelling alternative for clinical diagnostics. Here, Genomic Medicine Sweden-Rare Diseases has explored the utility of HiFi Revio long-read genome sequencing (lrGS) for digital karyotyping of SVs nationwide. The 16 samples from 13 families were collected from all Swedish healthcare regions. Prior investigations had identified 16 SVs, ranging from simple to complex rearrangements, including inversions, translocations, and copy number variants. We have established a national pipeline and a shared variant database for variant calling and filtering. Using lrGS, 14 of the 16 known SVs are detected. Of these, 13 are mapped at nucleotide resolution, and one complex rearrangement is only visible by read depth. Two Chromosome 21 rearrangements, one mosaic, remain undetected. Average read lengths are 8.3-18.8 kb with coverage exceeding 20× for all samples. De novo assembly results in a limited number of phased contigs per individual (N50 6-86 Mb), enabling direct characterization of the chromosomal rearrangements. In a national pilot study, we demonstrate the utility of HiFi Revio lrGS for analyzing chromosomal rearrangements. Based on our results, we propose a 5-year plan to expand lrGS use for rare disease diagnostics in Sweden.

Resolving complex duplication variants in autism spectrum disorder using long-read genome sequencing

Rare or de novo structural variation, primarily in the form of copy number variants, is detected in 5%-10% of autism spectrum disorder (ASD) families. While complex structural variants involving duplications can generally be detected using microarray or short-read genome sequencing (GS), these methods frequently fail to characterize breakpoints at nucleotide resolution, requiring additional molecular methods for validation and fine-mapping. Here, we use Oxford Nanopore Technologies PromethION long-read GS to characterize complex genomic rearrangements (CGRs) involving large duplications that segregate with ASD in five families. In total, we investigated 13 CGR carriers and were able to resolve all breakpoint junctions at nucleotide resolution. While all breakpoints were identified, the precise genomic architecture of one rearrangement remained unresolved with three different potential structures. The findings in two families include potential fusion genes formed through duplication rearrangements, involving IL1RAPL1-DMD and SUPT16H-CHD8 In two of the families originating from the same geographical region, an identical rearrangement involving ANK2 was identified, which likely represents a founder variant. In addition, we analyze methylation status directly from the long-read data, allowing us to assess the activity of rearranged genes and regulatory regions. Investigation of methylation across the CGRs reveals aberrant methylation status in carriers across a rearrangement affecting the CREBBP locus. In aggregate, our results demonstrate the utility of nanopore sequencing to pinpoint CGRs associated with ASD in five unrelated families, and highlight the importance of a gene-centric description of disease-associated complex chromosomal rearrangements.

Leveraging the T2T assembly to resolve rare and pathogenic inversions in reference genome gaps

Chromosomal inversions (INVs) are particularly challenging to detect due to their copy-number neutral state and association with repetitive regions. Inversions represent about 1/20 of all balanced structural chromosome aberrations and can lead to disease by gene disruption or altering regulatory regions of dosage-sensitive genes in cis Short-read genome sequencing (srGS) can only resolve ∼70% of cytogenetically visible inversions referred to clinical diagnostic laboratories, likely due to breakpoints in repetitive regions. Here, we study 12 inversions by long-read genome sequencing (lrGS) (n = 9) or srGS (n = 3) and resolve nine of them. In four cases, the inversion breakpoint region was missing from at least one of the human reference genomes (GRCh37, GRCh38, T2T-CHM13) and a reference agnostic analysis was needed. One of these cases, an INV9 mappable only in de novo assembled lrGS data using T2T-CHM13 disrupts EHMT1 consistent with a Mendelian diagnosis (Kleefstra syndrome 1; MIM#610253). Next, by pairwise comparison between T2T-CHM13, GRCh37, and GRCh38, as well as the chimpanzee and bonobo, we show that hundreds of megabases of sequence are missing from at least one human reference, highlighting that primate genomes contribute to genomic diversity. Aligning population genomic data to these regions indicated that these regions are variable between individuals. Our analysis emphasizes that T2T-CHM13 is necessary to maximize the value of lrGS for optimal inversion detection in clinical diagnostics. These results highlight the importance of leveraging diverse and comprehensive reference genomes to resolve unsolved molecular cases in rare diseases.

CDK4 is co-amplified with either TP53 promoter gene fusions or MDM2 through distinct mechanisms in osteosarcoma

Amplification of the MDM2 and CDK4 genes on chromosome 12 is commonly associated with low-grade osteosarcomas. In this study, we conducted high-resolution genomic and transcriptomic analyses on 33 samples from 25 osteosarcomas, encompassing both high- and low-grade cases with MDM2 and/or CDK4 amplification. We discerned four major subgroups, ranging from nearly intact genomes to heavily rearranged ones, each harbouring CDK4 and MDM2 amplification or CDK4 amplification with TP53 structural alterations. While amplicons involving MDM2 exhibited signs of an initial chromothripsis event, no evidence of chromothripsis was found in TP53-rearranged cases. Instead, the initial disruption of the TP53 locus led to co-amplification of the CDK4 locus. Additionally, we observed recurring promoter swapping events involving the regulatory regions of the FRS2, PLEKHA5, and TP53 genes. These events resulted in ectopic expression of partner genes, with the ELF1 gene being upregulated by the FRS2 and TP53 promoter regions in two distinct cases.

IDeRare: a lightweight and extensible open-source phenotype and exome analysis pipeline for germline rare disease diagnosis

Objectives

Diagnosing rare diseases is an arduous and challenging process in clinical settings, resulting in the late discovery of novel variants and referral loops. To help clinicians, we built IDeRare pipelines to accelerate phenotype-genotype analysis for patients with suspected rare diseases.

Materials and Methods

IDeRare pipeline is separated into phenotype and genotype parts. The phenotype utilizes our handmade Python library, while the genotype part utilizes command line (bash) and Python script to combine bioinformatics executable and Docker image.

Results

We described various implementations of IDeRare phenotype and genotype parts with real-world clinical and exome data using IDeRare, accelerating the terminology conversion process and giving insight on the diagnostic pathway based on disease linkage analysis until exome analysis and HTML-based reporting for clinicians.

Conclusion

IDeRare is freely available under the BSD-3 license, obtainable via GitHub. The portability of IDeRare pipeline could be easily implemented for semi-technical users and extensible for advanced users.

The impact of next-generation sequencing for diagnosis and disease understanding of myeloid malignancies

INTRODUCTION

Defining the chromosomal and molecular changes associated with myeloid neoplasms (MNs) optimizes clinical care through improved diagnosis, prognosis, treatment planning, and patient monitoring. This review will concisely describe the techniques used to profile MNs clinically today, with descriptions of challenges and emerging approaches that may soon become standard-of-care.

AREAS COVERED

In this review, the authors discuss molecular assessment of MNs using non-sequencing techniques, including conventional cytogenetic analysis, fluorescence in situ hybridization, chromosomal genomic microarray testing; as well as DNA- or RNA-based next-generation sequencing (NGS) assays; and sequential monitoring via digital PCR or measurable residual disease assays. The authors explain why distinguishing somatic from germline alleles is critical for optimal management. Finally, they introduce emerging technologies, such as long-read, whole exome/genome, and single-cell sequencing, which are reserved for research purposes currently but will become clinical tests soon.

EXPERT OPINION

The authors describe challenges to the adoption of comprehensive genomic tests for those in resource-constrained environments and for inclusion into clinical trials. In the future, all aspects of patient care will likely be influenced by the adaptation of artificial intelligence and mathematical modeling, fueled by rapid advances in telecommunications.

A combination of long- and short-read genomics reveals frequent p-arm breakpoints within chromosome 21 complex genomic rearrangements

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.

Scalable and efficient DNA sequencing analysis on different compute infrastructures aiding variant discovery

DNA variation analysis has become indispensable in many aspects of modern biomedicine, most prominently in the comparison of normal and tumor samples. Thousands of samples are collected in local sequencing efforts and public databases requiring highly scalable, portable, and automated workflows for streamlined processing. Here, we present nf-core/sarek 3, a well-established, comprehensive variant calling and annotation pipeline for germline and somatic samples. It is suitable for any genome with a known reference. We present a full rewrite of the original pipeline showing a significant reduction of storage requirements by using the CRAM format and runtime by increasing intra-sample parallelization. Both are leading to a 70% cost reduction in commercial clouds enabling users to do large-scale and cross-platform data analysis while keeping costs and CO2 emissions low. The code is available at https://nf-co.re/sarek.

A Rare Case of De Novo Beta-Thalassemia Diagnosed by Whole-Genome Sequencing in an Ethnically Danish Newborn

In 2020, a 2-month-old ethnically Danish girl was diagnosed with β-thalassemia after presenting with persistent jaundice. The peripheral blood smear showed significant aniso- and poikilocytosis, increased number of reticulocytes and erythroblastosis. Trio analysis of the index patient and both parents was performed by whole-genome sequencing. Here, both parents were found normal, however the analysis revealed an apparently de novo HBB:c.444A > C variant in the child. The child has recently been discharged three months after a successful bone marrow transplantation with a matched sibling-donor.

Long-read sequencing and optical mapping generates near T2T assemblies that resolves a centromeric translocation

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.

Rare coding variants in NOX4 link high ROS levels to psoriatic arthritis mutilans

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.

Whole genome case-control study of central nervous system toxicity due to antimicrobial drugs

A genetic predisposition to central nervous system (CNS) toxicity induced by antimicrobial drugs (antibiotics, antivirals, antifungals, and antiparasitic drugs) has been suspected. Whole genome sequencing of 66 cases and 833 controls was performed to investigate whether antimicrobial drug-induced CNS toxicity was associated with genetic variation. The primary objective was to test whether antimicrobial-induced CNS toxicity was associated with seventeen efflux transporters at the blood-brain barrier. In this study, variants or structural elements in efflux transporters were not significantly associated with CNS toxicity. Secondary objectives were to test whether antimicrobial-induced CNS toxicity was associated with genes over the whole genome, with HLA, or with structural genetic variation. Uncommon variants in and close to three genes were significantly associated with CNS toxicity according to a sequence kernel association test combined with an optimal unified test (SKAT-O). These genes were LCP1 (q = 0.013), RETSAT (q = 0.013) and SFMBT2 (q = 0.035). Two variants were driving the LCP1 association: rs6561297 (p = 1.15x10-6, OR: 4.60 [95% CI: 2.51-8.46]) and the regulatory variant rs10492451 (p = 1.15x10-6, OR: 4.60 [95% CI: 2.51-8.46]). No common genetic variant, HLA-type or structural variation was associated with CNS toxicity. In conclusion, CNS toxicity due to antimicrobial drugs was associated with uncommon variants in LCP1, RETSAT and SFMBT2.

Disruption of the TP53 locus in osteosarcoma leads to TP53 promoter gene fusions and restoration of parts of the TP53 signalling pathway

TP53 is the most frequently mutated gene in human cancer. This gene shows not only loss-of-function mutations but also recurrent missense mutations with gain-of-function activity. We have studied the primary bone malignancy osteosarcoma, which harbours one of the most rearranged genomes of all cancers. This is odd since it primarily affects children and adolescents who have not lived the long life thought necessary to accumulate massive numbers of mutations. In osteosarcoma, TP53 is often disrupted by structural variants. Here, we show through combined whole-genome and transcriptome analyses of 148 osteosarcomas that TP53 structural variants commonly result in loss of coding parts of the gene while simultaneously preserving and relocating the promoter region. The transferred TP53 promoter region is fused to genes previously implicated in cancer development. Paradoxically, these erroneously upregulated genes are significantly associated with the TP53 signalling pathway itself. This suggests that while the classical tumour suppressor activities of TP53 are lost, certain parts of the TP53 signalling pathway that are necessary for cancer cell survival and proliferation are retained. In line with this, our data suggest that transposition of the TP53 promoter is an early event that allows for a new normal state of genome-wide rearrangements in osteosarcoma. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Osteosarcomas With Few Chromosomal Alterations or Adult Onset Are Genetically Heterogeneous

Osteosarcoma is the most common primary bone malignancy, often detected in children and adolescents and commonly associated with TP53 alterations along with a high number of chromosomal rearrangements. However, osteosarcoma can affect patients of any age, and some tumors display less genetic complexity. Besides TP53 variants, data on key driving mutations are lacking for many osteosarcomas, particularly those affecting adults. To detect osteosarcoma-specific alterations, we screened transcriptomic and genomic sequencing and copy number data from 150 bone tumors originally diagnosed as osteosarcomas. To increase the precision in gene fusion detection, we developed a bioinformatic tool denoted as NAFuse, which extracts gene fusions that are verified at both the genomic and transcriptomic levels. Apart from the already reported genetic subgroups of osteosarcoma with TP53 structural variants, or MDM2 and/or CDK4 amplification, we did not identify any recurrent genetic driver that signifies the remaining cases. Among the plethora of mutations identified, we found genetic alterations characteristic of, or similar to, those of other bone and soft tissue tumors in 8 cases. These mutations were found in tumors with relatively few other genetic alterations or in adults. Due to the lack of clinical context and available tissue, we can question the diagnosis only on a genetic basis. However, our findings support the notion that osteosarcomas with few chromosomal alterations or adult onset seem genetically distinct from conventional osteosarcomas of children and adolescents.

Rare Deletions or Large Duplications Contribute to Genetic Variation in Patients with Severe Tinnitus and Meniere Disease

Meniere disease (MD) is a debilitating disorder of the inner ear defined by sensorineural hearing loss (SNHL) associated with episodes of vertigo and tinnitus. Severe tinnitus, which occurs in around 1% of patients, is a multiallelic disorder associated with a burden of rare missense single nucleotide variants in synaptic genes. Rare structural variants (SVs) may also contribute to MD and severe tinnitus. In this study, we analyzed exome sequencing data from 310 MD Spanish patients and selected 75 patients with severe tinnitus based on a Tinnitus Handicap Inventory (THI) score > 68. Three rare deletions were identified in two unrelated individuals overlapping the ERBB3 gene in the positions: NC_000012.12:g.56100028_56100172del, NC_000012.12:g.56100243_56101058del, and NC_000012.12:g.56101359_56101526del. Moreover, an ultra-rare large duplication was found covering the AP4M1, COPS6, MCM7, TAF6, MIR106B, MIR25, and MIR93 genes in another two patients in the NC_000007.14:g.100089053_100112257dup region. All the coding genes exhibited expression in brain and inner ear tissues. These results confirm the contribution of large SVs to severe tinnitus in MD and pinpoint new candidate genes to get a better molecular understanding of the disease.

Structural insights into coordinating 5S RNP rotation with ITS2 pre-RNA processing during ribosome formation

The rixosome defined in Schizosaccharomyces pombe and humans performs diverse roles in pre-ribosomal RNA processing and gene silencing. Here, we isolate and describe the conserved rixosome from Chaetomium thermophilum, which consists of two sub-modules, the sphere-like Rix1-Ipi3-Ipi1 and the butterfly-like Las1-Grc3 complex, connected by a flexible linker. The Rix1 complex of the rixosome utilizes Sda1 as landing platform on nucleoplasmic pre-60S particles to wedge between the 5S rRNA tip and L1-stalk, thereby facilitating the 180° rotation of the immature 5S RNP towards its mature conformation. Upon rixosome positioning, the other sub-module with Las1 endonuclease and Grc3 polynucleotide-kinase can reach a strategic position at the pre-60S foot to cleave and 5' phosphorylate the nearby ITS2 pre-rRNA. Finally, inward movement of the L1 stalk permits the flexible Nop53 N-terminus with its AIM motif to become positioned at the base of the L1-stalk to facilitate Mtr4 helicase-exosome participation for completing ITS2 removal. Thus, the rixosome structure elucidates the coordination of two central ribosome biogenesis events, but its role in gene silencing may adapt similar strategies.

Quartet DNA reference materials and datasets for comprehensively evaluating germline variant calling performance

BACKGROUND

Genomic DNA reference materials are widely recognized as essential for ensuring data quality in omics research. However, relying solely on reference datasets to evaluate the accuracy of variant calling results is incomplete, as they are limited to benchmark regions. Therefore, it is important to develop DNA reference materials that enable the assessment of variant detection performance across the entire genome.

RESULTS

We established a DNA reference material suite from four immortalized cell lines derived from a family of parents and monozygotic twins. Comprehensive reference datasets of 4.2 million small variants and 15,000 structural variants were integrated and certified for evaluating the reliability of germline variant calls inside the benchmark regions. Importantly, the genetic built-in-truth of the Quartet family design enables estimation of the precision of variant calls outside the benchmark regions. Using the Quartet reference materials along with study samples, batch effects are objectively monitored and alleviated by training a machine learning model with the Quartet reference datasets to remove potential artifact calls. Moreover, the matched RNA and protein reference materials and datasets from the Quartet project enables cross-omics validation of variant calls from multiomics data.

CONCLUSIONS

The Quartet DNA reference materials and reference datasets provide a unique resource for objectively assessing the quality of germline variant calls throughout the whole-genome regions and improving the reliability of large-scale genomic profiling.

GeneToCN: an alignment-free method for gene copy number estimation directly from next-generation sequencing reads

Genomes exhibit large regions with segmental copy number variation, many of which include entire genes and are multiallelic. We have developed a computational method GeneToCN that counts the frequencies of gene-specific k-mers in FASTQ files and uses this information to infer copy number of the gene. We validated the copy number predictions for amylase genes (AMY1, AMY2A, AMY2B) using experimental data from digital droplet PCR (ddPCR) on 39 individuals and observed a strong correlation (R = 0.99) between GeneToCN predictions and experimentally determined copy numbers. An additional validation on FCGR3 genes showed a higher concordance for FCGR3A compared to two other methods, but reduced accuracy for FCGR3B. We further tested the method on three different genomic regions (SMN, NPY4R, and LPA Kringle IV-2 domain). Predicted copy number distributions of these genes in a set of 500 individuals from the Estonian Biobank were in good agreement with the previously published studies. In addition, we investigated the possibility to use GeneToCN on sequencing data generated by different technologies by comparing copy number predictions from Illumina, PacBio, and Oxford Nanopore data of the same sample. Despite the differences in variability of k-mer frequencies, all three sequencing technologies give similar predictions with GeneToCN.

Precision medicine in rare diseases: What is next?

Molecular diagnostics is a cornerstone of modern precision medicine, broadly understood as tailoring an individual's treatment, follow-up, and care based on molecular data. In rare diseases (RDs), molecular diagnoses reveal valuable information about the cause of symptoms, disease progression, familial risk, and in certain cases, unlock access to targeted therapies. Due to decreasing DNA sequencing costs, genome sequencing (GS) is emerging as the primary method for precision diagnostics in RDs. Several ongoing European initiatives for precision medicine have chosen GS as their method of choice. Recent research supports the role for GS as first-line genetic investigation in individuals with suspected RD, due to its improved diagnostic yield compared to other methods. Moreover, GS can detect a broad range of genetic aberrations including those in noncoding regions, producing comprehensive data that can be periodically reanalyzed for years to come when further evidence emerges. Indeed, targeted drug development and repurposing of medicines can be accelerated as more individuals with RDs receive a molecular diagnosis. Multidisciplinary teams in which clinical specialists collaborate with geneticists, genomics education of professionals and the public, and dialogue with patient advocacy groups are essential elements for the integration of precision medicine into clinical practice worldwide. It is also paramount that large research projects share genetic data and leverage novel technologies to fully diagnose individuals with RDs. In conclusion, GS increases diagnostic yields and is a crucial step toward precision medicine for RDs. Its clinical implementation will enable better patient management, unlock targeted therapies, and guide the development of innovative treatments.

Long-read whole-genome analysis of human single cells

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.

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.

Transposable element insertions in 1000 Swedish individuals

The majority of rare diseases are genetic, and regardless of advanced high-throughput genomics-based investigations, 60% of patients remain undiagnosed. A major factor limiting our ability to identify disease-causing alterations is a poor understanding of the morbid and normal human genome. A major genomic contributor of which function and distribution remain largely unstudied are the transposable elements (TE), which constitute 50% of our genome. Here we aim to resolve this knowledge gap and increase the diagnostic yield of rare disease patients investigated with clinical genome sequencing. To this end we characterized TE insertions in 1000 Swedish individuals from the SweGen dataset and 2504 individuals from the 1000 Genomes Project (1KGP), creating seven population-specific TE insertion databases. Of note, 66% of TE insertions in SweGen were present at >1% in the 1KGP databases, proving that most insertions are common across populations. Focusing on the rare TE insertions, we show that even though ~0.7% of those insertions affect protein coding genes, they rarely affect known disease casing genes (<0.1%). Finally, we applied a TE insertion identification workflow on two clinical cases where disease causing TE insertions were suspected and could verify the presence of pathogenic TE insertions in both. Altogether we demonstrate the importance of TE insertion detection and highlight possible clinical implications in rare disease diagnostics.

Whole-genome informed circulating tumor DNA analysis by multiplex digital PCR for disease monitoring in B-cell lymphomas: a proof-of-concept study

Introduction

Analyzing liquid biopsies for tumor-specific aberrations can facilitate detection of measurable residual disease (MRD) during treatment and at follow-up. In this study, we assessed the clinical potential of using whole-genome sequencing (WGS) of lymphomas at diagnosis to identify patient-specific structural (SVs) and single nucleotide variants (SNVs) to enable longitudinal, multi-targeted droplet digital PCR analysis (ddPCR) of cell-free DNA (cfDNA).

Methods

In 9 patients with B-cell lymphoma (diffuse large B-cell lymphoma and follicular lymphoma), comprehensive genomic profiling at diagnosis was performed by 30X WGS of paired tumor and normal specimens. Patient-specific multiplex ddPCR (m-ddPCR) assays were designed for simultaneous detection of multiple SNVs, indels and/or SVs, with a detection sensitivity of 0.0025% for SV assays and 0.02% for SNVs/indel assays. M-ddPCR was applied to analyze cfDNA isolated from serially collected plasma at clinically critical timepoints during primary and/or relapse treatment and at follow-up.

Results

A total of 164 SNVs/indels were identified by WGS including 30 variants known to be functionally relevant in lymphoma pathogenesis. The most frequently mutated genes included KMT2D, PIM1, SOCS1 and BCL2. WGS analysis further identified recurrent SVs including t(14;18)(q32;q21) (IGH::BCL2), and t(6;14)(p25;q32) (IGH::IRF4). Plasma analysis at diagnosis showed positive circulating tumor DNA (ctDNA) levels in 88% of patients and the ctDNA burden correlated with baseline clinical parameters (LDH and sedimentation rate, p-value <0.01). While clearance of ctDNA levels after primary treatment cycle 1 was observed in 3/6 patients, all patients analyzed at final evaluation of primary treatment showed negative ctDNA, hence correlating with PET-CT imaging. One patient with positive ctDNA at interim also displayed detectable ctDNA (average variant allele frequency (VAF) 6.9%) in the follow-up plasma sample collected 2 years after final evaluation of primary treatment and 25 weeks before clinical manifestation of relapse.

Conclusion

In summary, we demonstrate that multi-targeted cfDNA analysis, using a combination of SNVs/indels and SVs candidates identified by WGS analysis, provides a sensitive tool for MRD monitoring and can detect lymphoma relapse earlier than clinical manifestation.

Discovery of non-reference processed pseudogenes in the Swedish population

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.

Multi-organ landscape of therapy-resistant melanoma

Metastasis and failure of present-day therapies represent the most common causes of mortality in patients with cutaneous melanoma. To identify the underlying genetic and transcriptomic landscapes, in this study we analyzed multi-organ metastases and tumor-adjacent tissues from 11 rapid autopsies after treatment with MAPK inhibitor (MAPKi) and/or immune checkpoint blockade (ICB) and death due to acquired resistance. Either treatment elicits shared genetic alterations that suggest immune-evasive, cross-therapy resistance mechanisms. Large, non-clustered deletions, inversions and inter-chromosomal translocations dominate rearrangements. Analyzing data from separate melanoma cohorts including 345 therapy-naive patients and 35 patients with patient-matched pre-treatment and post-acquired resistance tumor samples, we performed cross-cohort analyses to identify MAPKi and ICB as respective contributors to gene amplifications and deletions enriched in autopsy versus therapy-naive tumors. In the autopsy cohort, private/late mutations and structural variants display shifted mutational and rearrangement signatures, with MAPKi specifically selecting for signatures of defective homologous-recombination, mismatch and base-excision repair. Transcriptomic signatures and crosstalks with tumor-adjacent macroenvironments nominated organ-specific adaptive pathways. An immune-desert, CD8+-macrophage-biased archetype, T-cell exhaustion and type-2 immunity characterized the immune contexture. This multi-organ analysis of therapy-resistant melanoma presents preliminary insights with potential to improve therapeutic strategies.

The cost-effectiveness of whole genome sequencing in neurodevelopmental disorders

Whole genome sequencing (WGS) has the potential to be a comprehensive genetic test, especially relevant for individuals with neurodevelopmental disorders, syndromes and congenital malformations. However, the cost consequences of using whole genome sequencing as a first-line genetic test for these individuals are not well understood. The study objective was to compare the healthcare costs and diagnostic yield when WGS is performed as the first-line test instead of chromosomal microarray analysis (CMA). Two cohorts were analyzed retrospectively using register data, cohort CMA (418 patients referred for CMA at the department of Clinical Genetics, Karolinska University Hospital, during 2015) and cohort WGS (89 patients included in a WGS-first prospective study in 2017). The analysis compared healthcare consumption over a 2-year period after referral for genetic testing, the diagnostic yield over a 2- and 3-year period after referral was also compiled. The mean healthcare cost per patient in cohort WGS was $2,339 lower compared to cohort CMA ($ - 2339, 95% CI - 12,238-7561; P = 0.64) including higher costs for genetic investigations ($1065, 95% CI 834-1295; P < 0.001) and lower costs for outpatient care ($ - 2330, 95% CI - 3992 to (- 669); P = 0.006). The diagnostic yield was 23% higher for cohort WGS (cohort CMA 20.1%, cohort WGS 24.7%) (0.046, 95% CI - 0.053-0.145; P = 0.36). WGS as a first-line diagnostic test for individuals with neurodevelopmental disorders is associated with statistically non-significant lower costs and higher diagnostic yield compared with CMA. This indicates that prioritizing WGS over CMA in health care decision making will yield positive expected outcomes as well as showing a need for further research.

Blocking Genomic Instability Prevents Acquired Resistance to MAPK Inhibitor Therapy in Melanoma

Blocking cancer genomic instability may prevent tumor diversification and escape from therapies. We show that, after MAPK inhibitor (MAPKi) therapy in patients and mice bearing patient-derived xenografts (PDX), acquired resistant genomes of metastatic cutaneous melanoma specifically amplify resistance-driver, nonhomologous end-joining (NHEJ), and homologous recombination repair (HRR) genes via complex genomic rearrangements (CGR) and extrachromosomal DNAs (ecDNA). Almost all sensitive and acquired-resistant genomes harbor pervasive chromothriptic regions with disproportionately high mutational burdens and significant overlaps with ecDNA and CGR spans. Recurrently, somatic mutations within ecDNA and CGR amplicons enrich for HRR signatures, particularly within acquired resistant tumors. Regardless of sensitivity or resistance, breakpoint-junctional sequence analysis suggests NHEJ as critical to double-stranded DNA break repair underlying CGR and ecDNA formation. In human melanoma cell lines and PDXs, NHEJ targeting by a DNA-PKCS inhibitor prevents/delays acquired MAPKi resistance by reducing the size of ecDNAs and CGRs early on combination treatment. Thus, targeting the causes of genomic instability prevents acquired resistance.

SIGNIFICANCE

Acquired resistance often results in heterogeneous, redundant survival mechanisms, which challenge strategies aimed at reversing resistance. Acquired-resistant melanomas recurrently evolve resistance-driving and resistance-specific amplicons via ecDNAs and CGRs, thereby nominating chromothripsis-ecDNA-CGR biogenesis as a resistance-preventive target. Specifically, targeting DNA-PKCS/NHEJ prevents resistance by suppressing ecDNA/CGR rearrangements in MAPKi-treated melanomas. This article is highlighted in the In This Issue feature, p. 799.

In Silico Identification and Functional Characterization of Genetic Variations across DLBCL Cell Lines

Diffuse large B-cell lymphoma (DLBCL) is the most common form of non-Hodgkin lymphoma and frequently develops through the accumulation of several genetic variations. With the advancement in high-throughput techniques, in addition to mutations and copy number variations, structural variations have gained importance for their role in genome instability leading to tumorigenesis. In this study, in order to understand the genetics of DLBCL pathogenesis, we carried out a whole-genome mutation profile analysis of eleven human cell lines from germinal-center B-cell-like (GCB-7) and activated B-cell-like (ABC-4) subtypes of DLBCL. Analysis of genetic variations including small sequence variants and large structural variations across the cell lines revealed distinct variation profiles indicating the heterogeneous nature of DLBCL and the need for novel patient stratification methods to design potential intervention strategies. Validation and prognostic significance of the variants was assessed using annotations provided for DLBCL samples in cBioPortal for Cancer Genomics. Combining genetic variations revealed new subgroups between the subtypes and associated enriched pathways, viz., PI3K-AKT signaling, cell cycle, TGF-beta signaling, and WNT signaling. Mutation landscape analysis also revealed drug-variant associations and possible effectiveness of known and novel DLBCL treatments. From the whole-genome-based mutation analysis, our findings suggest putative molecular genetics of DLBCL lymphomagenesis and potential genomics-driven precision treatments.

Independent Evolution of Sex Chromosomes and Male Pregnancy-Related Genes in Two Seahorse Species

Unlike birds and mammals, many teleosts have homomorphic sex chromosomes, and changes in the chromosome carrying the sex-determining locus, termed "turnovers", are common. Recent turnovers allow studies of several interesting questions. One question is whether the new sex-determining regions evolve to become completely non-recombining, and if so, how and why. Another is whether (as predicted) evolutionary changes that benefit one sex accumulate in the newly sex-linked region. To study these questions, we analyzed the genome sequences of two seahorse species of the Syngnathidae, a fish group in which many species evolved a unique structure, the male brood pouch. We find that both seahorse species have XY sex chromosome systems, but their sex chromosome pairs are not homologs, implying that at least one turnover event has occurred. The Y-linked regions occupy 63.9% and 95.1% of the entire sex chromosome of the two species and do not exhibit extensive sequence divergence with their X-linked homologs. We find evidence for occasional recombination between the extant sex chromosomes that may account for their homomorphism. We argue that these Y-linked regions did not evolve by recombination suppression after the turnover, but by the ancestral nature of the low crossover rates in these chromosome regions. With such an ancestral crossover landscape, a turnover can instantly create an extensive Y-linked region. Finally, we test for adaptive evolution of male pouch-related genes after they became Y-linked in the seahorse.

CONGA: Copy number variation genotyping in ancient genomes and low-coverage sequencing data

To date, ancient genome analyses have been largely confined to the study of single nucleotide polymorphisms (SNPs). Copy number variants (CNVs) are a major contributor of disease and of evolutionary adaptation, but identifying CNVs in ancient shotgun-sequenced genomes is hampered by typical low genome coverage (<1×) and short fragments (<80 bps), precluding standard CNV detection software to be effectively applied to ancient genomes. Here we present CONGA, tailored for genotyping CNVs at low coverage. Simulations and down-sampling experiments suggest that CONGA can genotype deletions >1 kbps with F-scores >0.75 at ≥1×, and distinguish between heterozygous and homozygous states. We used CONGA to genotype 10,002 outgroup-ascertained deletions across a heterogenous set of 71 ancient human genomes spanning the last 50,000 years, produced using variable experimental protocols. A fraction of these (21/71) display divergent deletion profiles unrelated to their population origin, but attributable to technical factors such as coverage and read length. The majority of the sample (50/71), despite originating from nine different laboratories and having coverages ranging from 0.44×-26× (median 4×) and average read lengths 52-121 bps (median 69), exhibit coherent deletion frequencies. Across these 50 genomes, inter-individual genetic diversity measured using SNPs and CONGA-genotyped deletions are highly correlated. CONGA-genotyped deletions also display purifying selection signatures, as expected. CONGA thus paves the way for systematic CNV analyses in ancient genomes, despite the technical challenges posed by low and variable genome coverage.

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.

Approaches to long-read sequencing in a clinical setting to improve diagnostic rate

Over the past decade, advances in genetic testing, particularly the advent of next-generation sequencing, have led to a paradigm shift in the diagnosis of molecular diseases and disorders. Despite our present collective ability to interrogate more than 90% of the human genome, portions of the genome have eluded us, resulting in stagnation of diagnostic yield with existing methodologies. Here we show how application of a new technology, long-read sequencing, has the potential to improve molecular diagnostic rates. Whole genome sequencing by long reads was able to cover 98% of next-generation sequencing dead zones, which are areas of the genome that are not interpretable by conventional industry-standard short-read sequencing. Through the ability of long-read sequencing to unambiguously call variants in these regions, we discovered an immunodeficiency due to a variant in IKBKG in a subject who had previously received a negative genome sequencing result. Additionally, we demonstrate the ability of long-read sequencing to detect small variants on par with short-read sequencing, its superior performance in identifying structural variants, and thirdly, its capacity to determine genomic methylation defects in native DNA. Though the latter technical abilities have been demonstrated, we demonstrate the clinical application of this technology to successfully identify multiple types of variants using a single test.

Genome sequencing is a sensitive first-line test to diagnose individuals with intellectual disability

PURPOSE

Individuals with intellectual disability (ID) and/or neurodevelopment disorders (NDDs) are currently investigated with several different approaches in clinical genetic diagnostics.

METHODS

We compared the results from 3 diagnostic pipelines in patients with ID/NDD: genome sequencing (GS) first (N = 100), GS as a secondary test (N = 129), or chromosomal microarray (CMA) with or without FMR1 analysis (N = 421).

RESULTS

The diagnostic yield was 35% (GS-first), 26% (GS as a secondary test), and 11% (CMA/FMR1). Notably, the age of diagnosis was delayed by 1 year when GS was performed as a secondary test and the cost per diagnosed individual was 36% lower with GS first than with CMA/FMR1. Furthermore, 91% of those with a negative result after CMA/FMR1 analysis (338 individuals) have not yet been referred for additional genetic testing and remain undiagnosed.

CONCLUSION

Our findings strongly suggest that genome analysis outperforms other testing strategies and should replace traditional CMA and FMR1 analysis as a first-line genetic test in individuals with ID/NDD. GS is a sensitive, time- and cost-effective method that results in a confirmed molecular diagnosis in 35% of all referred patients.

Bioinformatics: From NGS Data to Biological Complexity in Variant Detection and Oncological Clinical Practice

The use of next-generation sequencing (NGS) techniques for variant detection has become increasingly important in clinical research and in clinical practice in oncology. Many cancer patients are currently being treated in clinical practice or in clinical trials with drugs directed against specific genomic alterations. In this scenario, the development of reliable and reproducible bioinformatics tools is essential to derive information on the molecular characteristics of each patient’s tumor from the NGS data. The development of bioinformatics pipelines based on the use of machine learning and statistical methods is even more relevant for the determination of complex biomarkers. In this review, we describe some important technologies, computational algorithms and models that can be applied to NGS data from Whole Genome to Targeted Sequencing, to address the problem of finding complex cancer-associated biomarkers. In addition, we explore the future perspectives and challenges faced by bioinformatics for precision medicine both at a molecular and clinical level, with a focus on an emerging complex biomarker such as homologous recombination deficiency (HRD).

Patient-Specific Assays Based on Whole-Genome Sequencing Data to Measure Residual Disease in Children With Acute Lymphoblastic Leukemia: A Proof of Concept Study

Risk-adapted treatment in acute lymphoblastic leukemia (ALL) relies on genetic information and measurable residual disease (MRD) monitoring. In this proof of concept study, DNA from diagnostic bone marrow (BM) of six children with ALL, without stratifying genetics or central nervous system (CNS) involvement, underwent whole-genome sequencing (WGS) to identify structural variants (SVs) in the leukemic blasts. Unique sequences generated by SVs were targeted with patient-specific droplet digital PCR (ddPCR) assays. Genomic DNA (gDNA) from BM and cell-free DNA (cfDNA) from plasma and cerebrospinal fluid (CSF) were analyzed longitudinally. WGS with 30× coverage enabled target identification in all cases. Limit of quantifiability (LoQ) and limit of detection (LoD) for the ddPCR assays (n = 15) were up to 10^-5 and 10^-6, respectively. All targets were readily detectable in a multiplexed ddPCR with minimal DNA input (1 ng of gDNA) at a 10^-1 dilution, and targets for half of the patients were also detectable at a 10^-2 dilution. The level of MRD in BM at end of induction and end of consolidation block 1 was in a comparable range between ddPCR and clinical routine methods for samples with detectable residual disease, although our approach consistently detected higher MRD values for patients with B-cell precursor ALL. Additionally, several samples with undetectable MRD by flow cytometry were MRD-positive by ddPCR. In plasma, the level of leukemic targets decreased in cfDNA over time following the MRD level detected in BM. cfDNA was successfully extracted from all diagnostic CSF samples (n = 6), and leukemic targets were detected in half of these. The results suggest that our approach to design molecular assays, together with ddPCR quantification, is a technically feasible option for accurate MRD quantification and that cfDNA may contribute valuable information regarding MRD and low-grade CNS involvement.

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.

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.

svBreak: A New Approach for the Detection of Structural Variant Breakpoints Based on Convolutional Neural Network

Structural variation (SV) is an important type of genome variation and confers susceptibility to human cancer diseases. Systematic analysis of SVs has become a crucial step for the exploration of mechanisms and precision diagnosis of cancers. The central point is how to accurately detect SV breakpoints by using next-generation sequencing (NGS) data. Due to the cooccurrence of multiple types of SVs in the human genome and the intrinsic complexity of SVs, the discrimination of SV breakpoint types is a challenging task. In this paper, we propose a convolutional neural network- (CNN-) based approach, called svBreak, for the detection and discrimination of common types of SV breakpoints. The principle of svBreak is that it extracts a set of SV-related features for each genome site from the sequencing reads aligned to the reference genome and establishes a data matrix where each row represents one site and each column represents one feature and then adopts a CNN model to analyze such data matrix for the prediction of SV breakpoints. The performance of the proposed approach is tested via simulation studies and application to a real sequencing sample. The experimental results demonstrate the merits of the proposed approach when compared with existing methods. Thus, svBreak can be expected to be a supplementary approach in the field of SV analysis in human tumor genomes.

Case Report: Whole genome sequencing identifies CCDC88C as a novel JAK2 fusion partner in pediatric T-cell acute lymphoblastic leukemia

In the present report, we applied whole genome sequencing (WGS) to genetically characterize a case of pediatric T-cell acute lymphoblastic leukemia (ALL) refractory to standard therapy. WGS identified a novel JAK2 fusion, with CCDC88C as a partner. CCDC88C encodes a protein part of the Wnt signaling pathway and has previously been described in hematological malignancies as fusion partner to FLT3 and PDGFRB. The novel CCDC88C::JAK2 fusion gene results in a fusion transcript, predicted to produce a hybrid protein, which retains the kinase domain of JAK2 and is expected to respond to JAK2 inhibitors. This report illustrates the potential of WGS in the diagnostic setting of ALL.

Massive parallel sequencing in individuals with multiple primary tumours reveals the benefit of re-analysis

Multiple primary cancers, defined as three or more primary tumours, are rare, and there are few genetic studies concerning them. There is a need for increased knowledge on the heritability of multiple primary cancers and genotype-phenotype correlations. We have performed whole-genome/exome sequencing (WGS/WES) in ten individuals with three or more primary tumours, with no previous findings on standard clinical genetic investigations. In one individual with a clinical diagnosis of MEN1, a likely pathogenic cryptic splice site variant was detected in the MEN1 gene. The variant (c.654C > A) is synonymous but we showed in a cDNA analysis that it affects splicing and leads to a frameshift, with the theoretical new amino acid sequence p.(Gly219Glufs*13). In one individual with metachronous colorectal cancers, ovarian cancer, endometrial cancer and chronic lymphocytic leukaemia, we found a likely pathogenic variant in the MLH1 gene (c.27G > A), and two risk factor variants in the genes CHEK2 and HOXB13. The MLH1 variant is synonymous but has previously been shown to be associated to constitutional low-grade hypermethylation of the MLH1 promoter, and segregates with disease in families with colorectal and endometrial cancer. No pathogenic single nucleotide or structural variants were detected in the remaining eight individuals in the study. The pathogenic variants found by WGS/WES were in genes already sequenced by Sanger sequencing and WES in the clinic, without any findings. We conclude that, in individuals with an unequivocal clinical diagnosis of a specific hereditary cancer syndrome, where standard clinical testing failed to detect a causative variant, re-analysis may lead to a diagnosis.

Deciphering complex genome rearrangements in C. elegans using short-read whole genome sequencing

Genomic rearrangements cause congenital disorders, cancer, and complex diseases in human. Yet, they are still understudied in rare diseases because their detection is challenging, despite the advent of whole genome sequencing (WGS) technologies. Short-read (srWGS) and long-read WGS approaches are regularly compared, and the latter is commonly recommended in studies focusing on genomic rearrangements. However, srWGS is currently the most economical, accurate, and widely supported technology. In Caenorhabditis elegans (C. elegans), such variants, induced by various mutagenesis processes, have been used for decades to balance large genomic regions by preventing chromosomal crossover events and allowing the maintenance of lethal mutations. Interestingly, those chromosomal rearrangements have rarely been characterized on a molecular level. To evaluate the ability of srWGS to detect various types of complex genomic rearrangements, we sequenced three balancer strains using short-read Illumina technology. As we experimentally validated the breakpoints uncovered by srWGS, we showed that, by combining several types of analyses, srWGS enables the detection of a reciprocal translocation (eT1), a free duplication (sDp3), a large deletion (sC4), and chromoanagenesis events. Thus, applying srWGS to decipher real complex genomic rearrangements in model organisms may help designing efficient bioinformatics pipelines with systematic detection of complex rearrangements in human genomes.

Detection of somatic structural variants from short-read next-generation sequencing data

Somatic structural variants (SVs), which are variants that typically impact >50 nucleotides, play a significant role in cancer development and evolution but are notoriously more difficult to detect than small variants from short-read next-generation sequencing (NGS) data. This is due to a combination of challenges attributed to the purity of tumour samples, tumour heterogeneity, limitations of short-read information from NGS and sequence alignment ambiguities. In spite of active development of SV detection tools (callers) over the past few years, each method has inherent advantages and limitations. In this review, we highlight some of the important factors affecting somatic SV detection and compared the performance of seven commonly used SV callers. In particular, we focus on the extent of change in sensitivity and precision for detecting different SV types and size ranges from samples with differing variant allele frequencies and sequencing depths of coverage. We highlight the reasons for why some SV callers perform well in some settings but not others, allowing our evaluation findings to be extended beyond the seven SV callers examined in this paper. As the importance of large SVs become increasingly recognized in cancer genomics, this paper provides a timely review on some of the most impactful factors influencing somatic SV detection that should be considered when choosing SV callers.

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.

Ethnic-specific association of amylase gene copy number with adiposity traits in a large Middle Eastern biobank

Studies assessing the impact of amylase genes copy number (CN) on adiposity report conflicting findings in different global populations, likely reflecting the impact of ancestral and ethnic-specific environment and lifestyle on selection at the amylase loci. Here, we leverage population size and detailed adiposity measures from a large population biobank to resolve confounding effects and determine the relationship between salivary (AMY1) and pancreatic (AMY2A) amylase genes CN and adiposity in 2935 Qatari individuals who underwent whole-genome sequencing (WGS) as part of the Qatar Genome Programme. We observe a negative association between AMY1 CNs and trunk fat percentage in the Qatari population (P = 7.50 × 10^-3) and show that Qataris of Arab descent have significantly lower CN at AMY1 (P = 1.32 × 10^-10) as well as less favorable adiposity and metabolic profiles (P < 1.34 × 10^-8) than Qataris with Persian ancestry. Indeed, lower AMY1 CN was associated with increased total and trunk fat percentages in Arabs (P < 4.60 × 10^-3) but not in Persians. Notably, overweight and obese Persians reported a significant trend towards dietary restraint following weight gain compared to Arabs (P = 4.29 × 10^-5), with AMY1 CN showing negative association with dietary self-restraint (P = 3.22 × 10^-3). This study reports an association between amylase gene CN and adiposity traits in a large Middle Eastern population. Importantly, we leverage rich biobank data to demonstrate that the strength of this association varies with ethnicity, and may be influenced by population-specific behaviors that also contribute to adiposity traits.

A Novel Mutation of the NARROW LEAF 1 Gene Adversely Affects Plant Architecture in Rice (Oryza sativa L.)

Plant architecture is critical for enhancing the adaptability and productivity of crop plants. Mutants with an altered plant architecture allow researchers to elucidate the genetic network and the underlying mechanisms. In this study, we characterized a novel nal1 rice mutant with short height, small panicle, and narrow and thick deep green leaves that was identified from a cross between a rice cultivar and a weedy rice accession. Bulked segregant analysis coupled with genome re-sequencing and cosegregation analysis revealed that the overall mutant phenotype was caused by a 1395-bp deletion spanning over the last two exons including the transcriptional end site of the nal1 gene. This deletion resulted in chimeric transcripts involving nal1 and the adjacent gene, which were validated by a reference-guided assembly of transcripts followed by PCR amplification. A comparative transcriptome analysis of the mutant and the wild-type rice revealed 263 differentially expressed genes involved in cell division, cell expansion, photosynthesis, reproduction, and gibberellin (GA) and brassinosteroids (BR) signaling pathways, suggesting the important regulatory role of nal1. Our study indicated that nal1 controls plant architecture through the regulation of genes involved in the photosynthetic apparatus, cell cycle, and GA and BR signaling pathways.