Large Structural Variants Research Articles

How structural variants shape avian phenotypes: Lessons from model systems.

Despite receiving significant recent attention, the relevance of structural variation (SV) in driving phenotypic diversity remains understudied, although recent advances in long-read sequencing, bioinformatics and pangenomic approaches have enhanced SV detection. We review the role of SVs in shaping phenotypes in avian model systems, and identify some general patterns in SV type, length and their associated traits. We found that most of the avian SVs so far identified are short indels in chickens, which are frequently associated with changes in body weight and plumage colouration. Overall, we found that relatively short SVs are more frequently detected, likely due to a combination of their prevalence compared to large SVs, and a detection bias, stemming primarily from the widespread use of short-read sequencing and associated analytical methods. SVs most commonly involve non-coding regions, especially introns, and when patterns of inheritance were reported, SVs associated primarily with dominant discrete traits. We summarise several examples of phenotypic convergence across different species, mediated by different SVs in the same or different genes and different types of changes in the same gene that can lead to various phenotypes. Complex rearrangements and supergenes, which can simultaneously affect and link several genes, tend to have pleiotropic phenotypic effects. Additionally, SVs commonly co-occur with single-nucleotide polymorphisms, highlighting the need to consider all types of genetic changes to understand the basis of phenotypic traits. We end by summarising expectations for when long-read technologies become commonly implemented in non-model birds, likely leading to an increase in SV discovery and characterisation. The growing interest in this subject suggests an increase in our understanding of the phenotypic effects of SVs in upcoming years.

Abstract 4336: Homologous recombination repair deficiency, and not recurrence, determines mutational burden and clonal dynamics in high grade serous ovarian cancer

Abstract High grade serous ovarian cancer (HGSOC) is the most lethal gynecologic malignancy, killing more than 9,000 women each year in the United States. Nearly 80% of patients with HGSOC will experience recurrence within 5 years, but little is known about the mechanisms that drive this process. Up to 60% of the genome in an HGSOC tumor is impacted by structural variant mutations and pronounced intratumoral heterogeneity. Intratumor heterogeneity is believed to be a key feature of recurrence and therapy resistance with some evidence indicating diverse and complex mechanisms of the seeding of metastatic sites and recurrent tumors. We included 35 HGSOC patients for whom paired chemo-naive and chemoresistant tumors as well as germline DNA was available. Whole genome sequencing was used to identify somatic single nucleotide variants (SNV) insertion/deletion variants (indels), and structural variants (SVs) in each tumor. The number of somatic SVs and indels was higher in tumors with homologous recombination repair deficiency (HRD; P=0.000252 SVs, P=0.00385 Indel) but was not affected by recurrent status (P=0.826 SV, P=0.145 SNV/indel). Clonal composition and dynamics were measured using SNVs/indels as tumors progressed from chemo-naïve primary samples to recurrent chemo-resistant tumors. Surprisingly, few changes were observed in clonal abundance and complexity through progression. When this analysis was repeated with SVs homologous recombination repair proficient tumors tend to be polyclonal while HRD tumors tend to be monoclonal. Five denovo structural variant signatures were identified in our cohort and led to three distinct structural variant classes tumors defined by deletion, tumors defined by duplication, and tumors defined by copy number neutral changes (inversions, translocations and balanced complex events). Patients with tumors defined by deletions and copy number neutral changes carry pathogenic SVs correlating with reduced survival (P=0.00023), while DNA gains are less subject to pathogenic alterations. Each class displayed distinct regions of the chromosome that was frequently affected by large scale SV events (>5Mb), which was also observed in the HRP/HRD status. Although no regions were notably altered across recurrence, GO analysis revealed that recurrent tumors have a significantly reduced immune response, which was not seen in the primary tumors. Within a subset of tumors, we are plan to utilize Oxford Nanopore Technologies ultra long read sequencing which will allow for confirmation of SVs observed in the WGS as well as identification of more complex events that cannot be deciphered using WGS techniques. Citation Format: Michael A. Diaz, Nicole Gull, Pei-Chen Peng, Kate Lawrenson, Bobbie J. Rimel, Jenny Lester, Beth Karlan, Simon A. Gayther, Michelle R. Jones. Homologous recombination repair deficiency, and not recurrence, determines mutational burden and clonal dynamics in high grade serous ovarian cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 4336.

Common mitochondrial deletions in RNA-Seq: evaluation of bulk, single-cell, and spatial transcriptomic datasets

Common mitochondrial DNA (mtDNA) deletions are large structural variants in the mitochondrial genome that accumulate in metabolically active tissues with age and have been investigated in various diseases. We applied the Splice-Break2 pipeline (designed for high-throughput quantification of mtDNA deletions) to human RNA-Seq datasets and describe the methodological considerations for evaluating common deletions in bulk, single-cell, and spatial transcriptomics datasets. A robust evaluation of 1570 samples from 14 RNA-Seq studies showed: (i) the abundance of some common deletions detected in PCR-amplified mtDNA correlates with levels observed in RNA-Seq data; (ii) RNA-Seq library preparation method has a strong effect on deletion detection; (iii) deletions had a significant, positive correlation with age in brain and muscle; (iv) deletions were enriched in cortical grey matter, specifically in layers 3 and 5; and (v) brain regions with dopaminergic neurons (i.e., substantia nigra, ventral tegmental area, and caudate nucleus) had remarkable enrichment of common mtDNA deletions.

Characterisation of a LINE-1 Insertion in the RP1 Gene by Targeted Adaptive Nanopore Sequencing in a Family with Retinitis Pigmentosa

Retinitis pigmentosa (RP) is a group of inherited degenerative retinal disorders affecting more than 1.5 million people worldwide. For 30-50% of individuals with RP, the genetic cause remains unresolved by current clinical diagnostic gene panels. It is likely explained by variants in novel RP-associated genes or noncoding regulatory regions, or by complex genetic alterations such as large structural variants. Recent developments in long-read sequencing techniques have opened an opportunity for efficient analysis of complex genetic variants. We analysed a Finnish family with dominantly inherited RP affecting six individuals in three generations. Two affected individuals underwent a comprehensive clinical examination in combination with a clinical diagnostic gene panel, followed by whole exome sequencing in our laboratory. They exhibited typical signs of RP, yet initial sequence analysis found no causative variants. Reanalysis of the sequencing data detected a LINE-1 (L1) retrotransposon insertion of unknown size in exon 4 of the RP1 axonemal microtubule-associated (RP1) gene. The large chimeric L1 insertion that segregated with the disease was further characterised using targeted adaptive nanopore sequencing of RP1, allowing us to identify a 5.6 kb L1 transposable element insertion in RP1 as the cause of RP in this family with dominantly inherited RP.

Improved sensitivity for detection of pathogenic variants in familial NF2-related schwannomatosis

PurposeTo determine the impact of additional genetic screening techniques on the rate of detection of pathogenic variants leading to familial NF2-related schwannomatosis.MethodsWe conducted genetic screening of a cohort of 168...

Abstract PR005: Tracing the clonal dynamic of metastatic castration-resistant prostate cancer over immunotherapy using circulating tumor DNA

Abstract Human cancer tissues are known to release their own DNA into body fluids such as plasma. These DNA molecules derived from tumors, collectively known as circulating tumor DNA (ctDNA), contain genetic information about disease progression at the time of sampling. Importantly, ctDNA has recently been utilized for capturing the whole-body disease profile and when sampled at various time points over the treatment, can inform the dynamic of clonal heterogeneity and its contribution to therapy failures. Metastatic castration-resistant prostate cancer (mCRPC) is a class of cancerous diseases with a high degree of clonal heterogeneity and high mortality which is largely due to the frequent emergence of therapy resistance. Therefore, reconstructing the clonal history for mCRPC diseases using ctDNA has great potential in identification of novel and targetable molecular mechanisms that drive mCRPC progression. To assess this possibility, we hypothesize that tumor clones that are responsible for therapy resistance carry distinct genomic profiles, which are captured in plasma ctDNA and could be computationally resolved through genomic sequencing and clonal reconstruction. To address this question, we obtained plasma cell-free DNA (cfDNA), a mixture of ctDNA and other tissue-derived DNA, from 38 individuals treated with combination PD-L1 and PARP inhibition involved in a recent clinical trial (NCT02484404). Whole-genome sequencing was performed using cfDNA and fragmented buffy coat DNA as germline control. cfDNA samples with a tumor fraction of lower than 10% were excluded. Various computational strategies were used to model the clonal structures of diseases and estimate the temporal order in which small and large structural variants emerged. Through this approach, we observed a negative association between cfDNA tumor fraction and therapy response. Multiple clonal evolutionary patterns leading to therapy failures were observed. For example, both gain and loss of activating androgen receptor mutations were detected in different cases following the treatment, which represents two distinct mechanisms of evading immunotherapy targeting. As a result, changes in the relevant transcriptional activities were also detected through ctDNA fragment-based computational inference. Additionally, clonal persistence was detected in multiple cases that exhibited stable diseases upon the treatment, which involves loss-of-function mutations in p53, ATM, and CDK12. To estimate the potential contribution to plasma ctDNA by metastatic diseases, multi-region whole genome and exome sequencing data of the tumors that are available for some cases prior to the treatment were also analyzed. Clones resolved from analyzing ctDNA were computationally mapped to the primary disease or metastatic tissue with a shared set of mutations in a subset of cases and may explain the patterns of treatment response. Our novel findings from this study will fill in a critical gap of knowledge about the complex genetic mechanisms driving mCRPC immunotherapy resistance. Citation Format: Chennan Li, Anna Baj, Clara C. Y. Seo, Nicholas T. Terrigino, John R. Bright, S. Thomas Hennigan, Isaiah M. King, Scott Wilkinson, Shana Y. Trostel, William D. Figg, William L. Dahut, Jung Min Lee, David Y. Takeda, Fatima Karzai, Adam G. Sowalsky. Tracing the clonal dynamic of metastatic castration-resistant prostate cancer over immunotherapy using circulating tumor DNA [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Translating Cancer Evolution and Data Science: The Next Frontier; 2023 Dec 3-6; Boston, Massachusetts. Philadelphia (PA): AACR; Cancer Res 2024;84(3 Suppl_2):Abstract nr PR005.

Highly efficient CRISPR/Cas9-mediated exon skipping for recessive dystrophic epidermolysis bullosa.

Gene therapy based on the CRISPR/Cas9 system has emerged as a promising strategy for treating the monogenic fragile skin disorder recessive dystrophic epidermolysis bullosa (RDEB). With this approach problematic wounds could be grafted with gene edited, patient-specificskin equivalents. Precise gene editing using homology-directed repair (HDR) is the ultimate goal, however low efficiencies have hindered progress. Reframing strategies based on highly efficient non-homologous end joining (NHEJ) repair aimed at excising dispensable, mutation-harboring exons offer a promising alternative approach for restoring the COL7A1 open reading frame. To this end, we employed an exon skipping strategy using dual single guide RNA (sgRNA)/Cas9 ribonucleoproteins (RNPs) targeted at three novel COL7A1 exons (31, 68, and 109) containing pathogenic heterozygous mutations, and achieved exon deletion rates of up to 95%. Deletion of exon 31 in both primary human RDEB keratinocytes and fibroblasts resulted in the restoration of type VII collagen (C7), leading to increased cellular adhesion in vitro and accurate C7 deposition at the dermal-epidermal junction in a 3D skin model. Taken together, we extend the list of COL7A1 exons amenable to therapeutic deletion. As an incidental finding, we find that long-read Nanopore sequencing detected large on-target structural variants comprised of deletions up to >5 kb at a frequency of ~10%. Although this frequency may be acceptable given the high rates of intended editing outcomes, our data demonstrate that standard short-read sequencing may underestimate the full range of unexpected Cas9-mediated editing events.

Structural Elucidation of Sialylated O-Glycan Alditols Obtained from Mucins by Mass Spectrometry.

Acidic O-glycans having sialic acid and/or sulfate residue are abundantly expressed in intestinal mucins. However, structural elucidation of acidic O-glycans is a laborious and time-consuming task due to their large structural variations. Here, we describe a methodology of structural elucidation for sialylated O-glycan alditols from intestinal mucins using tandem mass spectroscopy. Methylesterification and mild periodate oxidation of sialylated O-glycan alditols assist mass analysis. This description includes the purification process of O-glycan alditols for structural analysis.

Light-driven dynamical tuning of the thermal conductivity in ferroelectrics.

Dynamical tuning of the thermal conductivity in crystals, κ, is critical for thermal management applications, as well as for energy harvesting and the development of novel phononic devices able to perform logic operations with phonons. Such a desired κ control can be achieved in functional materials that experience large structural and phonon variations as a result of field-induced phase transformations. However, this approach is only practical within reduced temperature intervals containing zero-bias phase transition points, since otherwise the necessary driving fields become excessively large and the materials' performances are detrimentally affected. Here, based on first-principles calculations, we propose an alternative strategy for dynamically tuning κ that is operative over broad temperature conditions and realizable in a wide class of materials. By shining light on the archetypal perovskite oxide KNbO3, we predict that ultrafast and reversible ferroelectric-to-paraelectric phase transitions are induced, yielding large and anisotropic κ variations (up to ≈30% at T = 300 K). These light-induced thermal transport shifts can take place at temperatures spanning several hundreds of kelvin and are essentially the result of anharmonic effects affecting the phonon lifetimes.

Genetic findings of individuals with inherited retinal diseases in Finland

Genetic research can benefit from population isolates like Finland, as they tend to have a high concentration of deleterious alleles on a small number of low‐frequency variants. To investigate gene findings from the patients with inherited retinal disorders (IRD) in Finland, retrospectively collected data from the patients with IRDs treated between January 1st of 2012, and December 31st of 2020 at the Department of Ophthalmology, Helsinki University Hospital. Of the 638 patients that underwent genetic testing, 79% received a likely genetic diagnosis (504/638). The genetic findings corresponded to autosomal recessive (AR) inheritance in 242 patients, X‐linked (XL) in 158, autosomal dominant (AD) in 94, and mitochondrial inheritance in 22 patients. The most common genetic diagnosis was X‐linked retinoschisis as 17% (84/504) harboured pathogenic variants in RS1. The most common retinitis pigmentosa genes with causative variants were CERKL (8.5%, 43/504), RPGR (5.8%, 24/504), and EYS (5%, 25/504). CERKL variant c.375C>G, p.(Cys125Trp) and EYS c.1155T>A, p.(Cys385Ter) are frequent founder mutations. The Stargardt disease is less common than in other Western countries, the pathogenic variants in the ABCA4 gene were found in 24 patients. Pathogenic variants were also found in CHM (26), and CLRN1 (16). The genetic landscape of IRDs in Finland differs from other Western countries due to its specific isolated population history. The high yield of genetic findings is possibly explained by the enriched founder mutations in the population. One‐fifth of the patients are lacking genetic diagnosis, which might be due to shortages in the testing methods, i.e., non‐coding variants and complex structural genomic variants are missing. To address this problem, we have been utilizing the whole genome and long‐read sequencing methods to uncover new pathogenic variants. In two families with IRDs, we have found large structural variants in the RP1 and PRPH2 genes using Nanopore long‐read sequencing. In the future, we try to utilize the advantages of the population isolate to find unknown genes causing IRDs.

Combined analysis of transposable elements and structural variation in maize genomes reveals genome contraction outpaces expansion.

Structural differences between genomes are a major source of genetic variation that contributes to phenotypic differences. Transposable elements, mobile genetic sequences capable of increasing their copy number and propagating themselves within genomes, can generate structural variation. However, their repetitive nature makes it difficult to characterize fine-scale differences in their presence at specific positions, limiting our understanding of their impact on genome variation. Domesticated maize is a particularly good system for exploring the impact of transposable element proliferation as over 70% of the genome is annotated as transposable elements. High-quality transposable element annotations were recently generated for de novo genome assemblies of 26 diverse inbred maize lines. We generated base-pair resolved pairwise alignments between the B73 maize reference genome and the remaining 25 inbred maize line assemblies. From this data, we classified transposable elements as either shared or polymorphic in a given pairwise comparison. Our analysis uncovered substantial structural variation between lines, representing both simple and complex connections between TEs and structural variants. Putative insertions in SNP depleted regions, which represent recently diverged identity by state blocks, suggest some TE families may still be active. However, our analysis reveals that within these recently diverged genomic regions, deletions of transposable elements likely account for more structural variation events and base pairs than insertions. These deletions are often large structural variants containing multiple transposable elements. Combined, our results highlight how transposable elements contribute to structural variation and demonstrate that deletion events are a major contributor to genomic differences.

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.

Benchmarking long-read genome sequence alignment tools for human genomics applications.

The utility of long-read genome sequencing platforms has been shown in many fields including whole genome assembly, metagenomics, and amplicon sequencing. Less clear is the applicability of long reads to reference-guided human genomics, which is the foundation of genomic medicine. Here, we benchmark available platform-agnostic alignment tools on datasets from nanopore and single-molecule real-time platforms to understand their suitability in producing a genome representation. For this study, we leveraged publicly-available data from sample NA12878 generated on Oxford Nanopore and sample NA24385 on Pacific Biosciences platforms. We employed state of the art sequence alignment tools including GraphMap2, long-read aligner (LRA), Minimap2, CoNvex Gap-cost alignMents for Long Reads (NGMLR), and Winnowmap2. Minimap2 and Winnowmap2 were computationally lightweight enough for use at scale, while GraphMap2 was not. NGMLR took a long time and required many resources, but produced alignments each time. LRA was fast, but only worked on Pacific Biosciences data. Each tool widely disagreed on which reads to leave unaligned, affecting the end genome coverage and the number of discoverable breakpoints. No alignment tool independently resolved all large structural variants (1,001-100,000 base pairs) present in the Database of Genome Variants (DGV) for sample NA12878 or the truthset for NA24385. These results suggest a combined approach is needed for LRS alignments for human genomics. Specifically, leveraging alignments from three tools will be more effective in generating a complete picture of genomic variability. It should be best practice to use an analysis pipeline that generates alignments with both Minimap2 and Winnowmap2 as they are lightweight and yield different views of the genome. Depending on the question at hand, the data available, and the time constraints, NGMLR and LRA are good options for a third tool. If computational resources and time are not a factor for a given case or experiment, NGMLR will provide another view, and another chance to resolve a case. LRA, while fast, did not work on the nanopore data for our cluster, but PacBio results were promising in that those computations completed faster than Minimap2. Due to its significant burden on computational resources and slow run time, Graphmap2 is not an ideal tool for exploration of a whole human genome generated on a long-read sequencing platform.

Platinum Talen-Mediated Knock-in of Single-Stranded DNA Templates Facilitates Manufacturing of Clinical-Scale Non-Virally Gene-Edited T Cells from Peripheral Blood

Background: During the past decade, technologies that manufacture gene-modified immune cells, such as chimeric antigen receptor T cells, have been rapidly evolved and introduced in human clinics. T cells reprogrammed to express antigen-specific T-cell receptors (TCR-T) are an up-and-coming tool to treat intractable infection or malignant neoplasm in a personalized manner. However, several caveats still exist in the production of TCR-T for clinical use, because competitive surface expression of endogenous and transduced TCRs might hamper desired expression of transduced TCR. Moreover, mispairing of TCR subunits could mediate unexpected harmful immune complications. Therefore, in this study, we aimed to develop a method to introduce any desired TCRs in human primary T cells independent of interference with endogenous TCR by way of targeted knockout of TRA and TRB gene loci using highly-efficient transcription activator-like effector nuclease (TALEN), named Platinum TALEN (Sakuma T, Sci Rep 2013). Methods: We designed and synthesized two pairs of 5'-Cap-1 Platinum TALEN mRNA targeting the gene sequences that encode constant regions of human TRA ( TRAC) and TRB ( TRBC). TCR knockout efficiencies by these TALEN pairs were assessed by flow cytometric analysis using anti-CD3/TCR antibodies. Next, we constructed single-stranded DNA (ssDNA) homology-directed repair templates for the knock-in of desired TCR transgene into the TALEN target sequence of the TRAC locus. As a model TCR, we selected a 1G4 clone that reacts with NY-ESO-1-derived peptide in an HLA-A*02-restricted manner. We delivered these TALEN mRNA pairs and 1G4-encoding ssDNA by electroporation into primary human T cells collected from healthy volunteer adult donors after 72 hours of anti-CD3/CD28 beads activation. Genome-edited T cells were subsequently cultured in Therapeak TM X-VIVO15 (Lonza) supplemented with 10% human AB serum in the presence of 2-mercaptoethanol and interleukin-2 followed by second anti-CD3/CD28 beads activation. Density of T cells were maintained approximately at 1.0-3.0 x10 6/mL and transferred to larger flasks with additional fresh media as appropriate according to cell counts throughout the culture.1G4-transduced T cells were detected and counted by HLA-A*02:01 NY-ESO-1 tetramers at 17-22 days after electroporation. Additionally, we performed optical genome mapping analysis (Bionano Genomics) of Platinum TALEN-treated cells to screen genome-editing-associated large chromosomal structural variations (SVs) including unintended translocation. Results: Platinum TALEN mRNA pairs targeting TRAC and TRBC yielded a high proportion (>70%) of CD3/TCR negative T cells without the cost of compromised cell viability. By co-electroporation of 1G4 ssDNA and these Platinum TALEN mRNA pairs, we successfully obtained an average of 10.5 (range 2.2-19.0) million (n=5)NY-ESO-1-tetramer-positive cells with high viability from 3 millionprimary T cells until 3 weeks after electroporation. By flow cytometric analysis, a proportion of 1G4-tetramer positive cells that express Vβ chains other than Vβ13 (Vβ of 1G4) was approximately 10%. The PD-1-positive fraction in the expanded 1G4-expressing T cells was less than 10%, while levels of TIM-3 expression were 40-60%. These tetramer-positive cells showed specific cytotoxicities against HLA-A*02-positive cells pulsed with cognate NY-ESO-1 peptides. In an optical genome mapping analysis of genomic DNA obtained from TALEN-treated T cells, we found unique SVs with relatively small sizes and low coverage but no large unique SVs with high coverage when compared with control genomic DNA. Additionally, no interchromosomal translocation was observed. Conclusions: We have successfully designed Platinum TALEN mRNA pairs targeting TRAC and TRBC which facilitate the production of genome-edited human primary T cells. We also developed an electroporation and expansion culture protocol that enables us to produce viable genome-edited 1G4-transduced T cells at a large cell dose equivalent to a clinical scale. Collectively, these results suggest that targeted orthotopic knock-in of antigen-specific TCR-encoding ssDNA into the TCR locus by the use of Platinum TALEN is a promising strategy that can be applied for the clinical manufacturing of therapeutic TCR-T cells.

A pangenome graph reference of 30 chicken genomes allows genotyping of large and complex structural variants

BackgroundThe red junglefowl, the wild outgroup of domestic chickens, has historically served as a reference for genomic studies of domestic chickens. These studies have provided insight into the etiology of traits of commercial importance. However, the use of a single reference genome does not capture diversity present among modern breeds, many of which have accumulated molecular changes due to drift and selection. While reference-based resequencing is well-suited to cataloging simple variants such as single-nucleotide changes and short insertions and deletions, it is mostly inadequate to discover more complex structural variation in the genome.MethodsWe present a pangenome for the domestic chicken consisting of thirty assemblies of chickens from different breeds and research lines.ResultsWe demonstrate how this pangenome can be used to catalog structural variants present in modern breeds and untangle complex nested variation. We show that alignment of short reads from 100 diverse wild and domestic chickens to this pangenome reduces reference bias by 38%, which affects downstream genotyping results. This approach also allows for the accurate genotyping of a large and complex pair of structural variants at the K feathering locus using short reads, which would not be possible using a linear reference.ConclusionsWe expect that this new paradigm of genomic reference will allow better pinpointing of exact mutations responsible for specific phenotypes, which will in turn be necessary for breeding chickens that meet new sustainability criteria and are resilient to quickly evolving pathogen threats.

Short-Read Sequencing Detects Large Structural Variants in Cancer Genomes.

Most large structural variants (SV) can be detected by short-read sequencing (SRS) of cancer genomes.

Accurate identification of structural variations from cancer samples.

Structural variations (SVs) are commonly found in cancer genomes. They can cause gene amplification, deletion and fusion, among other functional consequences. With an average read length of hundreds of kilobases, nano-channel-based optical DNA mapping is powerful in detecting large SVs. However, existing SV calling methods are not tailored for cancer samples, which have special properties such as mixed cell types and sub-clones. Here we propose the Cancer Optical Mapping for detecting Structural Variations (COMSV) method that is specifically designed for cancer samples. It shows high sensitivity and specificity in benchmark comparisons. Applying to cancer cell lines and patient samples, COMSV identifies hundreds of novel SVs per sample.

Unveiling local and global conformational changes and allosteric communications in SOD1 systems using molecular dynamics simulation and network analyses

BackgroundAmyotrophic lateral sclerosis (ALS) is a serious neurodegenerative disorder affecting nerve cells in the brain and spinal cord that is caused by mutations in the superoxide dismutase 1 (SOD1) enzyme. ALS-related mutations cause misfolding, dimerisation instability, and increased formation of aggregates. The underlying allosteric mechanisms, however, remain obscure as far as details of their fundamental atomistic structure are concerned. Hence, this gap in knowledge limits the development of novel SOD1 inhibitors and the understanding of how disease-associated mutations in distal sites affect enzyme activity.Methods: We combined microsecond-scale based unbiased molecular dynamics (MD) simulation with network analysis to elucidate the local and global conformational changes and allosteric communications in SOD1 Apo (unmetallated form), Holo, Apo_CallA (mutant and unmetallated form), and Holo_CallA (mutant form) systems. To identify hotspot residues involved in SOD1 signalling and allosteric communications, we performed network centrality, community network, and path analyses.Results: Structural analyses showed that unmetallated SOD1 systems and cysteine mutations displayed large structural variations in the catalytic sites, affecting structural stability. Inter- and intra H-bond analyses identified several important residues crucial for maintaining interfacial stability, structural stability, and enzyme catalysis. Dynamic motion analysis demonstrated more balanced atomic displacement and highly correlated motions in the Holo system. The rationale for structural disparity observed in the disulfide bond formation and R143 configuration in Apo and Holo systems were elucidated using distance and dihedral probability distribution analyses.Conclusion: Our study highlights the efficiency of combining extensive MD simulations with network analyses to unravel the features of protein allostery.

Targeted haplotyping in pharmacogenomics using Oxford Nanopore Technologies' adaptive sampling.

Pharmacogenomics (PGx) studies the impact of interindividual genomic variation on drug response, allowing the opportunity to tailor the dosing regimen for each patient. Current targeted PGx testing platforms are mainly based on microarray, polymerase chain reaction, or short-read sequencing. Despite demonstrating great value for the identification of single nucleotide variants (SNVs) and insertion/deletions (INDELs), these assays do not permit identification of large structural variants, nor do they allow unambiguous haplotype phasing for star-allele assignment. Here, we used Oxford Nanopore Technologies' adaptive sampling to enrich a panel of 1,036 genes with well-documented PGx relevance extracted from the Pharmacogenomics Knowledge Base (PharmGKB). By evaluating concordance with existing truth sets, we demonstrate accurate variant and star-allele calling for five Genome in a Bottle reference samples. We show that up to three samples can be multiplexed on one PromethION flow cell without a significant drop in variant calling performance, resulting in 99.35% and 99.84% recall and precision for the targeted variants, respectively. This work advances the use of nanopore sequencing in clinical PGx settings.

Most large structural variants in cancer genomes can be detected without long reads

Short-read sequencing is the workhorse of cancer genomics yet is thought to miss many structural variants (SVs), particularly large chromosomal alterations. To characterize missing SVs in short-read whole genomes, we analyzed ‘loose ends’—local violations of mass balance between adjacent DNA segments. In the landscape of loose ends across 1,330 high-purity cancer whole genomes, most large (>10-kb) clonal SVs were fully resolved by short reads in the 87% of the human genome where copy number could be reliably measured. Some loose ends represent neotelomeres, which we propose as a hallmark of the alternative lengthening of telomeres phenotype. These pan-cancer findings were confirmed by long-molecule profiles of 38 breast cancer and melanoma cases. Our results indicate that aberrant homologous recombination is unlikely to drive the majority of large cancer SVs. Furthermore, analysis of mass balance in short-read whole genome data provides a surprisingly complete picture of cancer chromosomal structure.