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

Single variant, yet "double trouble": TSC and KBG syndrome because of a large de novo inversion

Despite the advances in high-throughput sequencing, many rare disease patients remain undiagnosed. In particular, the patients with well-defined clinical phenotypes and established clinical diagnosis, yet missing or partial genetic diagnosis, may hold a clue to more complex genetic mechanisms of a disease that could be missed by available clinical tests. Here, we report a patient with a clinical diagnosis of Tuberous sclerosis, combined with unusual secondary features, but negative clinical tests including TSC1 and TSC2 Short-read whole-genome sequencing combined with advanced bioinformatics analyses were successful in uncovering a de novo pericentric 87-Mb inversion with breakpoints in TSC2 and ANKRD11, which explains the TSC clinical diagnosis, and confirms a second underlying monogenic disorder, KBG syndrome. Our findings illustrate how complex variants, such as large inversions, may be missed by clinical tests and further highlight the importance of well-defined clinical diagnoses in uncovering complex molecular mechanisms of a disease, such as complex variants and "double trouble" effects.

Real-world diagnostic outcomes and cost-effectiveness of genome-wide sequencing for developmental and seizure disorders: Evidence from Canada

PURPOSE

To determine real-world diagnostic rates, cost trajectories, and cost-effectiveness of exome sequencing (ES) and genome sequencing (GS) for children with developmental and/or seizure disorders in British Columbia, Canada.

METHODS

Based on medical records review, we estimated real-world costs and outcomes for 491 patients who underwent standard of care (SOC) diagnostic testing at British Columbia Children's Hospital. Results informed a state-transition Markov model examining cost-effectiveness of 3 competing diagnostic strategies: (1) SOC with last-tier access to ES, (2) streamlined ES access, and (3) first-tier GS.

RESULTS

Through SOC, 49.4% (95% CI: 40.6, 58.2) of patients were diagnosed at an average cost of C$11,683 per patient (95% CI: 9200, 14,166). Compared with SOC, earlier ES or GS access yielded similar or improved diagnostic rates and shorter times to genetic diagnosis, with 94% of simulations demonstrating cost savings for streamlined ES and 60% for first-tier GS. Net benefit from the perspective of the health care system was C$2956 (95% CI: -608, 6519) for streamlined ES compared with SOC.

CONCLUSION

Using real-world data, we found earlier access to ES may yield more rapid genetic diagnosis of childhood developmental and seizure disorders and cost savings compared with current practice in a Canadian health care system.

Molecular basis of essential and morphological variations across 12 balancer strains in C. elegans

Balancers are primarily chromosomal rearrangements that allow lethal or sterile mutations to be stably maintained as heterozygotes. Strains with balanced lethal/sterile mutations are available at the Caenorhabditis Genetics Center. Such strains harbor morphological markers, with associated molecular changes, that are in trans to the balancer. In many cases, only the genetic position (in cM) has been described for balanced mutations or morphological markers. We used short-read whole genome sequencing to uncover the genomic position of those variants (balanced mutations and linked markers) and predicted their effect. We investigated 12 different strains, and characterized at a molecular level 12 variants.

A novel FAME1 repeat configuration in a European family identified using a combined genomics approach

Familial adult myoclonic epilepsy (FAME) is an adult-onset neurological disease characterized by cortical tremor, myoclonus, and seizures due to a pentanucleotide repeat expansion: a combination of pathogenic TTTCA expansion associated with a TTTTA repeat in introns of six different genes. Repeat-primed PCR (RP-PCR) is an inexpensive test for expansions at known loci. The analysis of the SAMD12 locus revealed that the repeats have different size, configuration, and composition. The TTTCA repeats can be very long (>1000 repeats) but also very short (14 being the shortest identified). Here, we report siblings of European descent with the clinical diagnosis of FAME yet a negative RP-PCR test. Using short-read genome sequencing, we identified the pentanucleotide expansion in intron 4 of SAMD12, which was confirmed by CRIPSR-Cas9-mediated enrichment and long-read sequencing to be of (TTTTA)_~879 (TTTCA)₃ (TTTTA)₇ (TTTCA)₇ configuration. Our finding is the first to associate the SAMD12 locus in European patients with FAME and currently represents the shortest identified TTTCA expansion. Our results suggest that the SAMD12 locus should be tested in patients with suspected FAME independent of ethnicity. Furthermore, RP-PCR may miss the underlying mutation, and genome sequencing may be needed to confirm the pathogenic repeat.

The Power of Clinical Diagnosis for Deciphering Complex Genetic Mechanisms in Rare Diseases

Complex genetic disease mechanisms, such as structural or non-coding variants, currently pose a substantial difficulty in frontline diagnostic tests. They thus may account for most unsolved rare disease patients regardless of the clinical phenotype. However, the clinical diagnosis can narrow the genetic focus to just a couple of genes for patients with well-established syndromes defined by prominent physical and/or unique biochemical phenotypes, allowing deeper analyses to consider complex genetic origin. Then, clinical-diagnosis-driven genome sequencing strategies may expedite the development of testing and analytical methods to account for complex disease mechanisms as well as to advance functional assays for the confirmation of complex variants, clinical management, and the development of new therapies.

Genome sequencing of C. elegans balancer strains reveals previously unappreciated complex genomic rearrangements

Genetic balancers in Caenorhabditis elegans are complex variants that allow lethal or sterile mutations to be stably maintained in a heterozygous state by suppressing crossover events. Balancers constitute an invaluable tool in the C. elegans scientific community and have been widely used for decades. The first/traditional balancers were created by applying X-rays, UV, or gamma radiation on C. elegans strains, generating random genomic rearrangements. Their structures have been mostly explored with low-resolution genetic techniques (e.g., fluorescence in situ hybridization or PCR), before genomic mapping and molecular characterization through sequencing became feasible. As a result, the precise nature of most chromosomal rearrangements remains unknown, whereas, more recently, balancers have been engineered using the CRISPR-Cas9 technique for which the structure of the chromosomal rearrangement has been predesigned. Using short-read whole-genome sequencing (srWGS) and tailored bioinformatic analyses, we previously interpreted the structure of four chromosomal balancers randomly created by mutagenesis processes. Here, we have extended our analyses to five CRISPR-Cas9 balancers and 17 additional traditional balancing rearrangements. We detected and experimentally validated their breakpoints and have interpreted the balancer structures. Many of the balancers were found to be more intricate than previously described, being composed of complex genomic rearrangements (CGRs) such as chromoanagenesis-like events. Furthermore, srWGS revealed additional structural variants and CGRs not known to be part of the balancer genomes. Altogether, our study provides a comprehensive resource of complex genomic variations in C. elegans and highlights the power of srWGS to study the complexity of genomes by applying tailored analyses.

Different structural variant prediction tools yield considerably different results in Caenorhabditis elegans

The accurate characterization of structural variation is crucial for our understanding of how large chromosomal alterations affect phenotypic differences and contribute to genome evolution. Whole-genome sequencing is a popular approach for identifying structural variants, but the accuracy of popular tools remains unclear due to the limitations of existing benchmarks. Moreover, the performance of these tools for predicting variants in non-human genomes is less certain, as most tools were developed and benchmarked using data from the human genome. To evaluate the use of long-read data for the validation of short-read structural variant calls, the agreement between predictions from a short-read ensemble learning method and long-read tools were compared using real and simulated data from Caenorhabditis elegans. The results obtained from simulated data indicate that the best performing tool is contingent on the type and size of the variant, as well as the sequencing depth of coverage. These results also highlight the need for reference datasets generated from real data that can be used as 'ground truth' in benchmarks.

Epidemiology of Δ8THC-Related Carcinogenesis in USA: A Panel Regression and Causal Inferential Study

The use of Δ8THC is increasing at present across the USA in association with widespread cannabis legalization and the common notion that it is “legal weed”. As genotoxic actions have been described for many cannabinoids, we studied the cancer epidemiology of Δ8THC. Data on 34 cancer types was from the Centers for Disease Control Atlanta Georgia, substance abuse data from the Substance Abuse and Mental Health Services Administration, ethnicity and income data from the U.S. Census Bureau, and cannabinoid concentration data from the Drug Enforcement Agency, were combined and processed in R. Eight cancers (corpus uteri, liver, gastric cardia, breast and post-menopausal breast, anorectum, pancreas, and thyroid) were related to Δ8THC exposure on bivariate testing, and 18 (additionally, stomach, Hodgkins, and Non-Hodgkins lymphomas, ovary, cervix uteri, gall bladder, oropharynx, bladder, lung, esophagus, colorectal cancer, and all cancers (excluding non-melanoma skin cancer)) demonstrated positive average marginal effects on fully adjusted inverse probability weighted interactive panel regression. Many minimum E-Values (mEVs) were infinite. p-values rose from 8.04 × 10⁻⁷⁸. Marginal effect calculations revealed that 18 Δ8THC-related cancers are predicted to lead to a further 8.58 cases/100,000 compared to 7.93 for alcoholism and −8.48 for tobacco. Results indicate that between 8 and 20/34 cancer types were associated with Δ8THC exposure, with very high effect sizes (mEVs) and marginal effects after adjustment exceeding tobacco and alcohol, fulfilling the epidemiological criteria of causality and suggesting a cannabinoid class effect. The inclusion of pediatric leukemias and testicular cancer herein demonstrates heritable malignant teratogenesis.

Case Report: Biallelic Loss of Function ATM due to Pathogenic Synonymous and Novel Deep Intronic Variant c.1803-270T > G Identified by Genome Sequencing in a Child With Ataxia–Telangiectasia

Ataxia–telangiectasia (AT) is a complex neurodegenerative disease with an increased risk for bone marrow failure and malignancy. AT is caused by biallelic loss of function variants in ATM, which encodes a phosphatidylinositol 3-kinase that responds to DNA damage. Herein, we report a child with progressive ataxia, chorea, and genome instability, highly suggestive of AT. The clinical ataxia gene panel identified a maternal heterozygous synonymous variant (NM_000051.3: c.2250G > A), previously described to result in exon 14 skipping. Subsequently, trio genome sequencing led to the identification of a novel deep intronic variant [NG_009830.1(NM_000051.3): c.1803-270T > G] inherited from the father. Transcript analyses revealed that c.1803-270T > G results in aberrant inclusion of 56 base pairs of intron 11. In silico tests predicted a premature stop codon as a consequence, suggesting non-functional ATM; and DNA repair analyses confirmed functional loss of ATM. Our findings highlight the power of genome sequencing, considering deep intronic variants in undiagnosed rare disease patients.

Whole genome sequencing facilitates intragenic variant interpretation following modifier screening in C. elegans

Background

Intragenic modifiers (in-phase, second-site variants) are known to have dramatic effects on clinical outcomes, affecting disease attributes such as severity or age of onset. However, despite their clinical importance, the focus of many genetic screens in model systems is on the discovery of extragenic variants, with many labs still relying upon more traditional methods to identify modifiers. However, traditional methods such as PCR and Sanger sequencing can be time-intensive and do not permit a thorough understanding of the intragenic modifier effects in the context of non-isogenic genomic backgrounds.

Results

Here, we apply high throughput approaches to identify and understand intragenic modifiers using Caenorhabditis elegans. Specifically, we applied whole genome sequencing (WGS) to a mutagen-induced forward genetic screen to identify intragenic suppressors of a temperature-sensitive zyg-1(it25) allele in C. elegans. ZYG-1 is a polo kinase that is important for centriole function and cell divisions, and mutations that truncate its human orthologue, PLK4, have been associated with microcephaly. Combining WGS and CRISPR/Cas9, we rapidly identify intragenic modifiers, show that these variants are distributed non-randomly throughout zyg-1 and that genomic context plays an important role on phenotypic outcomes.

Conclusions

Ultimately, our work shows that WGS facilitates high-throughput identification of intragenic modifiers in clinically relevant genes by reducing hands-on research time and overall costs and by allowing thorough understanding of the intragenic phenotypic effects in the context of different genetic backgrounds.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08142-8.