Epigenetic targeting of PGBD5-dependent DNA damage in SMARCB1-deficient sarcomas.

Despite the potential of targeted epigenetic therapies, most cancers do not respond to current epigenetic drugs. The Polycomb repressive complex EZH2 inhibitor tazemetostat was recently approved for the treatment of SMARCB1-deficient epithelioid sarcomas, based on the functional antagonism between PRC2 and loss of SMARCB1. Through the analysis of tazemetostat-treated patient tumors, we recently defined key principles of their response and resistance to EZH2 epigenetic therapy. Here, using transcriptomic inference from SMARCB1-deficient tumor cells, we nominate the DNA damage repair kinase ATR as a target for rational combination EZH2 epigenetic therapy. We show that EZH2 inhibition promotes DNA damage in epithelioid and rhabdoid tumor cells, at least in part via its induction of the transposase-derived PGBD5. We leverage this collateral synthetic lethal dependency to target PGBD5-dependent DNA damage by inhibition of ATR but not CHK1 using elimusertib. Consequently, combined EZH2 and ATR inhibition improves therapeutic responses in diverse patient-derived epithelioid and rhabdoid tumors in vivo. This advances a combination epigenetic therapy based on EZH2-PGBD5 synthetic lethal dependency suitable for immediate translation to clinical trials for patients.

Childhood cancer mutagenesis caused by transposase-derived PGBD5.

Genomic rearrangements are a hallmark of most childhood tumors, including medulloblastoma, one of the most common brain tumors in children, but their causes remain largely unknown. Here, we show that PiggyBac transposable element derived 5 (Pgbd5) promotes tumor development in multiple developmentally accurate mouse models of Sonic Hedgehog (SHH) medulloblastoma. Most Pgbd5-deficient mice do not develop tumors, while maintaining normal cerebellar development. Ectopic activation of SHH signaling is sufficient to enforce cerebellar granule cell progenitor-like cell states, which exhibit Pgbd5-dependent expression of distinct DNA repair and neurodevelopmental factors. Mouse medulloblastomas expressing Pgbd5 have increased numbers of somatic structural DNA rearrangements, some of which carry PGBD5-specific sequences at their breakpoints. Similar sequence breakpoints recurrently affect somatic DNA rearrangements of known tumor suppressors and oncogenes in medulloblastomas in 329 children. This identifies PGBD5 as a medulloblastoma mutator and provides a genetic mechanism for the generation of oncogenic DNA rearrangements in childhood cancer.

Commonality and variation in mental representations of music revealed by a cross-cultural comparison of rhythm priors in 15 countries.

Music is present in every known society but varies from place to place. What, if anything, is universal to music cognition? We measured a signature of mental representations of rhythm in 39 participant groups in 15 countries, spanning urban societies and Indigenous populations. Listeners reproduced random 'seed' rhythms; their reproductions were fed back as the stimulus (as in the game of 'telephone'), such that their biases (the prior) could be estimated from the distribution of reproductions. Every tested group showed a sparse prior with peaks at integer-ratio rhythms. However, the importance of different integer ratios varied across groups, often reflecting local musical practices. Our results suggest a common feature of music cognition: discrete rhythm 'categories' at small-integer ratios. These discrete representations plausibly stabilize musical systems in the face of cultural transmission but interact with culture-specific traditions to yield the diversity that is evident when mental representations are probed across many cultures.

Overcoming clinical resistance to EZH2 inhibition using rational epigenetic combination therapy.

Epigenetic dependencies have become evident in many cancers. Based on antagonism between BAF/SWI/SNF and PRC2 in SMARCB1-deficient sarcomas, we recently completed the clinical trial of the EZH2 inhibitor tazemetostat. However, the principles of tumor response to epigenetic therapy in general, and tazemetostat in particular, remain unknown. Using functional genomics and diverse experimental models, we define molecular mechanisms of tazemetostat resistance in SMARCB1-deficient tumors. We found distinct acquired mutations that converge on the RB1/E2F axis and decouple EZH2-dependent differentiation and cell cycle control. This allows tumor cells to escape tazemetostat-induced G1 arrest, suggests a general mechanism for effective therapy, and provides prospective biomarkers for therapy stratification, including PRICKLE1. Based on this, we develop a combination strategy to circumvent tazemetostat resistance using bypass targeting of AURKB. This offers a paradigm for rational epigenetic combination therapy suitable for translation to clinical trials for epithelioid sarcomas, rhabdoid tumors, and other epigenetically dysregulated cancers.

Comprehensive and accurate genome analysis at scale using DRAGEN accelerated algorithms.

Research and medical genomics require comprehensive and scalable solutions to drive the discovery of novel disease targets, evolutionary drivers, and genetic markers with clinical significance. This necessitates a framework to identify all types of variants independent of their size (e.g., SNV/SV) or location (e.g., repeats). Here we present DRAGEN that utilizes novel methods based on multigenomes, hardware acceleration, and machine learning based variant detection to provide novel insights into individual genomes with ~30min computation time (from raw reads to variant detection). DRAGEN outperforms all other state-of-the-art methods in speed and accuracy across all variant types (SNV, indel, STR, SV, CNV) and further incorporates specialized methods to obtain key insights in medically relevant genes (e.g., HLA, SMN, GBA). We showcase DRAGEN across 3,202 genomes and demonstrate its scalability, accuracy, and innovations to further advance the integration of comprehensive genomics for research and medical applications.

A transposase-derived gene required for human brain development.

DNA transposable elements and transposase-derived genes are present in most living organisms, including vertebrates, but their function is largely unknown. PiggyBac Transposable Element Derived 5 (PGBD5) is an evolutionarily conserved vertebrate DNA transposase-derived gene with retained nuclease activity in cells. Vertebrate brain development is known to be associated with prominent neuronal cell death and DNA breaks, but their causes and functions are not well understood. Here, we show that PGBD5 contributes to normal brain development in mice and humans, where its deficiency causes disorder of intellectual disability, movement and seizures. In mice, Pgbd5 is required for the developmental induction of post-mitotic DNA breaks and recurrent somatic genome rearrangements in neurons. Together, these studies nominate PGBD5 as the long-hypothesized neuronal DNA nuclease required for brain function in mammals.

Consequences of Contracting COVID-19 or Taking the COVID-19 Vaccine for Individuals with a History of Lyme Disease.

Individuals with Lyme disease can be very symptomatic. This survey compares the burden of illness for individuals with a history of Lyme disease (HLD) with individuals with a HLD who have either contracted COVID-19 or who have taken the COVID-19 vaccine. The findings describe the relative symptom burden among these three groups using a cross-sectional descriptive survey investigating the burden of Lyme disease in a pandemic. The survey includes the General Symptom Questionnaire-30 (GSQ-30), a brief self-report scale designed to assess the symptom burden in Lyme disease (LD). The results of this survey show that the overall burden of illness among individuals with HLD is not significantly different after contracting COVID-19 or after COVID-19 vaccination. A new survey will be needed to better understand why one in five individuals with a HLD reported long COVID after contracting COVID-19. These results should help clinicians and their patients to discuss the consequences of contracting a COVID-19 infection or being vaccinated against COVID-19.

Overcoming clinical resistance to EZH2 inhibition using rational epigenetic combination therapy.

Essential epigenetic dependencies have become evident in many cancers. Based on the functional antagonism between BAF/SWI/SNF and PRC2 in SMARCB1-deficient sarcomas, we and colleagues recently completed the clinical trial of the EZH2 inhibitor tazemetostat. However, the principles of tumor response to epigenetic therapy in general, and tazemetostat in particular, remain unknown. Using functional genomics of patient tumors and diverse experimental models, we sought to define molecular mechanisms of tazemetostat resistance in SMARCB1-deficient sarcomas and rhabdoid tumors. We found distinct classes of acquired mutations that converge on the RB1/E2F axis and decouple EZH2-dependent differentiation and cell cycle control. This allows tumor cells to escape tazemetostat-induced G1 arrest despite EZH2 inhibition, and suggests a general mechanism for effective EZH2 therapy. This also enables us to develop combination strategies to circumvent tazemetostat resistance using cell cycle bypass targeting via AURKB, and synthetic lethal targeting of PGBD5-dependent DNA damage repair via ATR. This reveals prospective biomarkers for therapy stratification, including PRICKLE1 associated with tazemetostat resistance. In all, this work offers a paradigm for rational epigenetic combination therapy suitable for immediate translation to clinical trials for epithelioid sarcomas, rhabdoid tumors, and other epigenetically dysregulated cancers.

The complexity-aesthetics relationship for musical rhythm is more fixed than flexible: Evidence from children and expert dancers.

The urge to move to music (groove) depends in part on rhythmic syncopation in the music. For adults, the syncopation-groove relationship has an inverted-U shape: listeners want to move most to rhythms that have some, but not too much, syncopation. However, we do not know whether the syncopation-groove relationship is relatively sensitive to, or resistant to, a listener's experience. In two sets of experiments, we tested whether the syncopation-groove relationship is affected by dance experience or changes through development in childhood. Dancers and nondancers rated groove for 50 rhythmic patterns varying in syncopation. Dancers' and nondancers' ratings did not differ (and Bayesian tests provided substantial evidence that they were equivalent) in terms of mean groove and the optimal level of syncopation. Similarly, ballet and hip-hop dancers' syncopation-groove relationships did not differ. However, dancers had more robust syncopation-groove relationships (higher goodness-of-fit) than nondancers. Children (3-6 years old) completed two tasks to assess their syncopation-groove relationships: In a 2-alternative-forced choice task, children compared rhythms from 2 of 3 possible levels of syncopation (low, medium, and high) and chose which rhythm in a pair was better for dancing. In a dance task, children danced to the same rhythms. Results from both tasks indicated that for children, as for adults, medium syncopation rhythms elicit more groove than low syncopation rhythms. A follow-up experiment replicated the 2-alternative-forced choice task results. Taken together, the results suggest the optimal level of syncopation for groove is resistant to experience, although experience may affect the robustness of the inverted-U relationship. RESEARCH HIGHLIGHTS: In Experiment 1, dancers and nondancers rated groove (the urge to move) for musical rhythms, demonstrating the same inverted-U relationships between syncopation and groove. In Experiment 2, children and adults both chose rhythms with moderate syncopation more than low syncopation as more groove-inducing or better for dancing. Children also danced more for moderate than low syncopation, showing a close perception-behavior relationship across tasks. Similarities in the syncopation-groove relationship regardless of dance training and age suggest that this perceptual and behavioral groove response to rhythmic complexity may be quite resistant to experience.

Collective dynamics support group drumming, reduce variability, and stabilize tempo drift.

Humans are social animals who engage in a variety of collective activities requiring coordinated action. Among these, music is a defining and ancient aspect of human sociality. Human social interaction has largely been addressed in dyadic paradigms, and it is yet to be determined whether the ensuing conclusions generalize to larger groups. Studied more extensively in non-human animal behavior, the presence of multiple agents engaged in the same task space creates different constraints and possibilities than in simpler dyadic interactions. We addressed whether collective dynamics play a role in human circle drumming. The task was to synchronize in a group with an initial reference pattern and then maintain synchronization after it was muted. We varied the number of drummers from solo to dyad, quartet, and octet. The observed lower variability, lack of speeding up, smoother individual dynamics, and leader-less inter-personal coordination indicated that stability increased as group size increased, a sort of temporal wisdom of crowds. We propose a hybrid continuous-discrete Kuramoto model for emergent group synchronization with a pulse-based coupling that exhibits a mean field positive feedback loop. This research suggests that collective phenomena are among the factors that play a role in social cognition.

Undetectable very-low frequency sound increases dancing at a live concert.

Does low frequency sound (bass) make people dance more? Music that makes people want to move tends to have more low frequency sound, and bass instruments typically provide the musical pulse that people dance to1. Low pitches confer advantages in perception and movement timing, and elicit stronger neural responses for timing compared to high pitches2, suggesting superior sensorimotor communication. Low frequency sound is processed via vibrotactile3 and vestibular4 (in addition to auditory) pathways, and stimulation of these non-auditory modalities in the context of music can increase ratings of groove (the pleasurable urge to move to music)3, and modulate musical rhythm perception4. Anecdotal accounts describe intense physical and psychological effects of low frequencies, especially in electronic dance music5, possibly reflecting effects on physiological arousal. We do not, however, know if these associations extend to direct causal effects of low frequencies in complex, real-world, social contexts like dancing at concerts, or if low frequencies that are not consciously detectable can affect behaviour. We tested whether non-auditory low-frequency stimulation would increase audience dancing by turning very-low frequency (VLF) speakers on and off during a live electronic music concert and measuring audience members' movements using motion-capture. Movement increased when VLFs were present, and because the VLFs were below or near auditory thresholds (and a subsequent experiment suggested they were undetectable), we believe this represents an unconscious effect on behaviour, possibly via vestibular and/or tactile processing.

StructuralVariantAnnotation: a R/Bioconductor foundation for a caller-agnostic structural variant software ecosystem.

SUMMARY: StructuralVariantAnnotation is an R/Bioconductor package that provides a framework for decoupling downstream analysis of structural variant breakpoints from upstream variant calling methods. It standardizes the representational format from BEDPE, or any of the three different notations supported by VCF into a breakpoint GRanges data structure suitable for use by the wider Bioconductor ecosystem. It handles both transitive breakpoints and duplication/insertion notational differences of identical variants-both common scenarios when comparing short/long read-based call sets that confound downstream analysis. StructuralVariantAnnotation provides the caller-agnostic foundation needed for a R/Bioconductor ecosystem of structural variant annotation, classification and interpretation tools able to handle both simple and complex genomic rearrangements. AVAILABILITY AND IMPLEMENTATION: StructuralVariantAnnotation is implemented in R and available for download as the Bioconductor StructuralVariantAnnotation package. Details can be found at https://www.bioconductor.org/packages/release/bioc/html/StructuralVariantAnnotation.html. It has been released under a GPL license.

O-fucosylation of thrombospondin type 1 repeats is essential for ECM remodeling and signaling during bone development.

Many extracellular matrix (ECM) associated proteins that influence ECM properties have Thrombospondin type 1 repeats (TSRs) which are modified with O-linked fucose. The O-fucose is added in the endoplasmic reticulum to folded TSRs by the enzyme Protein O-fucosyltransferase-2 (POFUT2) and is proposed to promote efficient trafficking of substrates. The importance of this modification for function of TSR-proteins is underscored by the early embryonic lethality of mouse embryos lacking Pofut2. To overcome early lethality and investigate the impact of the Pofut2 knockout on the secretion of POFUT2 substrates and on extracellular matrix properties in vivo, we deleted Pofut2 in the developing limb mesenchyme using Prrx1-Cre recombinase. Loss of Pofut2 in the limb mesenchyme caused significant shortening of the limbs, long bones and tendons and stiff joint resembling the musculoskeletal dysplasias in human and in mice with mutations in ADAMTS or ADAMTSL proteins. Limb shortening was evident at embryonic day 14.5 where loss of O-fucosylation led to an accumulation of fibrillin 2 (FBN2), decreased BMP and IHH signaling, and increased TGF-ß signaling. Consistent with these changes we saw a decrease in the size of the hypertrophic zone with lower levels of Collagen-X. Unexpectedly, we observed minimal effects of the Pofut2 knockout on secretion of two POFUT2 substrates, CCN2 or ADAMTS17, in the developing bone. In contrast, CCN2 and two other POFUT2 substrates important for bone development, ADAMTS6 and 10, showed a decrease in secretion from POFUT2-null HEK293T cells in vitro. These combined results suggest that the impact of the Pofut2 mutation is cell-type specific. In addition, these observations raise the possibility that the O-fucose modification on TSRs extends beyond promoting efficient trafficking of POFUT2 substrates and has the potential to influence their function in the extracellular environment.

Unscrambling cancer genomes via integrated analysis of structural variation and copy number.

Complex somatic genomic rearrangements and copy number alterations are hallmarks of nearly all cancers. We have developed an algorithm, LINX, to aid interpretation of structural variant and copy number data derived from short-read, whole-genome sequencing. LINX classifies raw structural variant calls into distinct events and predicts their effect on the local structure of the derivative chromosome and the functional impact on affected genes. Visualizations facilitate further investigation of complex rearrangements. LINX allows insights into a diverse range of structural variation events and can reliably detect pathogenic rearrangements, including gene fusions, immunoglobulin enhancer rearrangements, intragenic deletions, and duplications. Uniquely, LINX also predicts chained fusions that we demonstrate account for 13% of clinically relevant oncogenic fusions. LINX also reports a class of inactivation events that we term homozygous disruptions that may be a driver mutation in up to 9% of tumors and may frequently affect PTEN, TP53, and RB1.

A multi-platform reference for somatic structural variation detection.

Accurate detection of somatic structural variation (SV) in cancer genomes remains a challenging problem. This is in part due to the lack of high-quality, gold-standard datasets that enable the benchmarking of experimental approaches and bioinformatic analysis pipelines. Here, we performed somatic SV analysis of the paired melanoma and normal lymphoblastoid COLO829 cell lines using four different sequencing technologies. Based on the evidence from multiple technologies combined with extensive experimental validation, we compiled a comprehensive set of carefully curated and validated somatic SVs, comprising all SV types. We demonstrate the utility of this resource by determining the SV detection performance as a function of tumor purity and sequence depth, highlighting the importance of assessing these parameters in cancer genomics projects. The truth somatic SV dataset as well as the underlying raw multi-platform sequencing data are freely available and are an important resource for community somatic benchmarking efforts.

svaRetro and svaNUMT: modular packages for annotating retrotransposed transcripts and nuclear integration of mitochondrial DNA in genome sequencing data.

Nuclear integration of mitochondrial genomes and retrocopied transcript insertion are biologically important but often-overlooked aspects of structural variant (SV) annotation. While tools for their detection exist, these typically rely on reanalysis of primary data using specialised detectors rather than leveraging calls from general purpose structural variant callers. Such reanalysis potentially leads to additional computational expense and does not take advantage of advances in general purpose structural variant calling. Here, we present svaRetro and svaNUMT; R packages that provide functions for annotating novel genomic events, such as nonreference retrocopied transcripts and nuclear integration of mitochondrial DNA. The packages were developed to work within the Bioconductor framework. We evaluate the performance of these packages to detect events using simulations and public benchmarking datasets, and annotate processed transcripts in a public structural variant database. svaRetro and svaNUMT provide modular, SV-caller agnostic tools for downstream annotation of structural variant calls.

An international virtual hackathon to build tools for the analysis of structural variants within species ranging from coronaviruses to vertebrates.

In October 2020, 62 scientists from nine nations worked together remotely in the Second Baylor College of Medicine & DNAnexus hackathon, focusing on different related topics on Structural Variation, Pan-genomes, and SARS-CoV-2 related research.   The overarching focus was to assess the current status of the field and identify the remaining challenges. Furthermore, how to combine the strengths of the different interests to drive research and method development forward. Over the four days, eight groups each designed and developed new open-source methods to improve the identification and analysis of variations among species, including humans and SARS-CoV-2. These included improvements in SV calling, genotyping, annotations and filtering. Together with advancements in benchmarking existing methods. Furthermore, groups focused on the diversity of SARS-CoV-2. Daily discussion summary and methods are available publicly at  https://github.com/collaborativebioinformatics provides valuable insights for both participants and the research community.

GRIDSS2: comprehensive characterisation of somatic structural variation using single breakend variants and structural variant phasing.

GRIDSS2 is the first structural variant caller to explicitly report single breakends-breakpoints in which only one side can be unambiguously determined. By treating single breakends as a fundamental genomic rearrangement signal on par with breakpoints, GRIDSS2 can explain 47% of somatic centromere copy number changes using single breakends to non-centromere sequence. On a cohort of 3782 deeply sequenced metastatic cancers, GRIDSS2 achieves an unprecedented 3.1% false negative rate and 3.3% false discovery rate and identifies a novel 32-100 bp duplication signature. GRIDSS2 simplifies complex rearrangement interpretation through phasing of structural variants with 16% of somatic calls phasable using paired-end sequencing.

VIRUSBreakend: Viral Integration Recognition Using Single Breakends.

MOTIVATION: Integration of viruses into infected host cell DNA can cause DNA damage and disrupt genes. Recent cost reductions and growth of whole genome sequencing has produced a wealth of data in which viral presence and integration detection is possible. While key research and clinically relevant insights can be uncovered, existing software has not achieved widespread adoption, limited in part due to high computational costs, the inability to detect a wide range of viruses, as well as precision and sensitivity. RESULTS: Here, we describe VIRUSBreakend, a high-speed tool that identifies viral DNA presence and genomic integration. It utilizes single breakends, breakpoints in which only one side can be unambiguously placed, in a novel virus-centric variant calling and assembly approach to identify viral integrations with high sensitivity and a near-zero false discovery rate. VIRUSBreakend detects viral integrations anywhere in the host genome including regions such as centromeres and telomeres unable to be called by existing tools. Applying VIRUSBreakend to a large metastatic cancer cohort, we demonstrate that it can reliably detect clinically relevant viral presence and integration including HPV, HBV, MCPyV, EBV and HHV-8. AVAILABILITY AND IMPLEMENTATION: VIRUSBreakend is part of the Genomic Rearrangement IDentification Software Suite (GRIDSS). It is available under a GPLv3 license from https://github.com/PapenfussLab/VIRUSBreakend. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Evolution of late-stage metastatic melanoma is dominated by aneuploidy and whole genome doubling.

Although melanoma is initiated by acquisition of point mutations and limited focal copy number alterations in melanocytes-of-origin, the nature of genetic changes that characterise lethal metastatic disease is poorly understood. Here, we analyze the evolution of human melanoma progressing from early to late disease in 13 patients by sampling their tumours at multiple sites and times. Whole exome and genome sequencing data from 88 tumour samples reveals only limited gain of point mutations generally, with net mutational loss in some metastases. In contrast, melanoma evolution is dominated by whole genome doubling and large-scale aneuploidy, in which widespread loss of heterozygosity sculpts the burden of point mutations, neoantigens and structural variants even in treatment-naïve and primary cutaneous melanomas in some patients. These results imply that dysregulation of genomic integrity is a key driver of selective clonal advantage during melanoma progression.