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.

Relationship between CD4 T cell turnover, cellular differentiation and HIV persistence during ART.

The precise role of CD4 T cell turnover in maintaining HIV persistence during antiretroviral therapy (ART) has not yet been well characterized. In resting CD4 T cell subpopulations from 24 HIV-infected ART-suppressed and 6 HIV-uninfected individuals, we directly measured cellular turnover by heavy water labeling, HIV reservoir size by integrated HIV-DNA (intDNA) and cell-associated HIV-RNA (caRNA), and HIV reservoir clonality by proviral integration site sequencing. Compared to HIV-negatives, ART-suppressed individuals had similar fractional replacement rates in all subpopulations, but lower absolute proliferation rates of all subpopulations other than effector memory (TEM) cells, and lower plasma IL-7 levels (p = 0.0004). Median CD4 T cell half-lives decreased with cell differentiation from naïve to TEM cells (3 years to 3 months, p<0.001). TEM had the fastest replacement rates, were most highly enriched for intDNA and caRNA, and contained the most clonal proviral expansion. Clonal proviruses detected in less mature subpopulations were more expanded in TEM, suggesting that they were maintained through cell differentiation. Earlier ART initiation was associated with lower levels of intDNA, caRNA and fractional replacement rates. In conclusion, circulating integrated HIV proviruses appear to be maintained both by slow turnover of immature CD4 subpopulations, and by clonal expansion as well as cell differentiation into effector cells with faster replacement rates.

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.

Auditory representation of learned sound sequences in motor regions of the macaque brain.

Human speech production requires the ability to couple motor actions with their auditory consequences. Nonhuman primates might not have speech because they lack this ability. To address this question, we trained macaques to perform an auditory-motor task producing sound sequences via hand presses on a newly designed device ("monkey piano"). Catch trials were interspersed to ascertain the monkeys were listening to the sounds they produced. Functional MRI was then used to map brain activity while the animals listened attentively to the sound sequences they had learned to produce and to two control sequences, which were either completely unfamiliar or familiar through passive exposure only. All sounds activated auditory midbrain and cortex, but listening to the sequences that were learned by self-production additionally activated the putamen and the hand and arm regions of motor cortex. These results indicate that, in principle, monkeys are capable of forming internal models linking sound perception and production in motor regions of the brain, so this ability is not special to speech in humans. However, the coupling of sounds and actions in nonhuman primates (and the availability of an internal model supporting it) seems not to extend to the upper vocal tract, that is, the supralaryngeal articulators, which are key for the production of speech sounds in humans. The origin of speech may have required the evolution of a "command apparatus" similar to the control of the hand, which was crucial for the evolution of tool use.

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.

ADAMTS9 and ADAMTS20 are differentially affected by loss of B3GLCT in mouse model of Peters plus syndrome.

Peters plus syndrome (MIM #261540 PTRPLS), characterized by defects in eye development, prominent forehead, hypertelorism, short stature and brachydactyly, is caused by mutations in the ß3-glucosyltransferase (B3GLCT) gene. Protein O-fucosyltransferase 2 (POFUT2) and B3GLCT work sequentially to add an O-linked glucose ß1-3fucose disaccharide to properly folded thrombospondin type 1 repeats (TSRs). Forty-nine proteins are predicted to be modified by POFUT2, and nearly half are members of the ADAMTS superfamily. Previous studies suggested that O-linked fucose is essential for folding and secretion of POFUT2-modified proteins and that B3GLCT-mediated extension to the disaccharide is essential for only a subset of targets. To test this hypothesis and gain insight into the origin of PTRPLS developmental defects, we developed and characterized two mouse B3glct knockout alleles. Using these models, we tested the role of B3GLCT in enabling function of ADAMTS9 and ADAMTS20, two highly conserved targets whose functions are well characterized in mouse development. The mouse B3glct mutants developed craniofacial and skeletal abnormalities comparable to PTRPLS. In addition, we observed highly penetrant hydrocephalus, white spotting and soft tissue syndactyly. We provide strong genetic and biochemical evidence that hydrocephalus and white spotting in B3glct mutants resulted from loss of ADAMTS20, eye abnormalities from partial reduction of ADAMTS9 and cleft palate from loss of ADAMTS20 and partially reduced ADAMTS9 function. Combined, these results provide compelling evidence that ADAMTS9 and ADAMTS20 were differentially sensitive to B3GLCT inactivation and suggest that the developmental defects in PTRPLS result from disruption of a subset of highly sensitive POFUT2/B3GLCT targets such as ADAMTS20.

Optical mapping reveals a higher level of genomic architecture of chained fusions in cancer.

Genomic rearrangements are common in cancer, with demonstrated links to disease progression and treatment response. These rearrangements can be complex, resulting in fusions of multiple chromosomal fragments and generation of derivative chromosomes. Although methods exist for detecting individual fusions, they are generally unable to reconstruct complex chained events. To overcome these limitations, we adopted a new optical mapping approach, allowing megabase-length genome maps to be reconstructed and rearranged genomes to be visualized without loss of integrity. Whole-genome mapping (Bionano Genomics) of a well-studied highly rearranged liposarcoma cell line resulted in 3338 assembled consensus genome maps, including 72 fusion maps. These fusion maps represent 112.3 Mb of highly rearranged genomic regions, illuminating the complex architecture of chained fusions, including content, order, orientation, and size. Spanning the junction of 147 chromosomal translocations, we found a total of 28 Mb of interspersed sequences that could not be aligned to the reference genome. Traversing these interspersed sequences using short-read sequencing breakpoint calls, we were able to identify and place 399 sequencing fragments within the optical mapping gaps, thus illustrating the complementary nature of optical mapping and short-read sequencing. We demonstrate that optical mapping provides a powerful new approach for capturing a higher level of complex genomic architecture, creating a scaffold for renewed interpretation of sequencing data of particular relevance to human cancer.

GRIDSS: sensitive and specific genomic rearrangement detection using positional de Bruijn graph assembly.

The identification of genomic rearrangements with high sensitivity and specificity using massively parallel sequencing remains a major challenge, particularly in precision medicine and cancer research. Here, we describe a new method for detecting rearrangements, GRIDSS (Genome Rearrangement IDentification Software Suite). GRIDSS is a multithreaded structural variant (SV) caller that performs efficient genome-wide break-end assembly prior to variant calling using a novel positional de Bruijn graph-based assembler. By combining assembly, split read, and read pair evidence using a probabilistic scoring, GRIDSS achieves high sensitivity and specificity on simulated, cell line, and patient tumor data, recently winning SV subchallenge #5 of the ICGC-TCGA DREAM8.5 Somatic Mutation Calling Challenge. On human cell line data, GRIDSS halves the false discovery rate compared to other recent methods while matching or exceeding their sensitivity. GRIDSS identifies nontemplate sequence insertions, microhomologies, and large imperfect homologies, estimates a quality score for each breakpoint, stratifies calls into high or low confidence, and supports multisample analysis.

Neural entrainment is associated with subjective groove and complexity for performed but not mechanical musical rhythms.

Both movement and neural activity in humans can be entrained by the regularities of an external stimulus, such as the beat of musical rhythms. Neural entrainment to auditory rhythms supports temporal perception, and is enhanced by selective attention and by hierarchical temporal structure imposed on rhythms. However, it is not known how neural entrainment to rhythms is related to the subjective experience of groove (the desire to move along with music or rhythm), the perception of a regular beat, the perception of complexity, and the experience of pleasure. In two experiments, we used musical rhythms (from Steve Reich's Clapping Music) to investigate whether rhythms that are performed by humans (with naturally variable timing) and rhythms that are mechanical (with precise timing), elicit differences in (1) neural entrainment, as measured by inter-trial phase coherence, and (2) subjective ratings of the complexity, preference, groove, and beat strength of rhythms. We also combined results from the two experiments to investigate relationships between neural entrainment and subjective perception of musical rhythms. We found that mechanical rhythms elicited a greater degree of neural entrainment than performed rhythms, likely due to the greater temporal precision in the stimulus, and the two types only elicited different ratings for some individual rhythms. Neural entrainment to performed rhythms, but not to mechanical ones, correlated with subjective desire to move and subjective complexity. These data, therefore, suggest multiple interacting influences on neural entrainment to rhythms, from low-level stimulus properties to high-level cognition and perception.

Digital PCR of Genomic Rearrangements for Monitoring Circulating Tumour DNA.

Identifying circulating tumour DNA (ctDNA) for monitoring of cancer therapy is dependent on the development of readily designed, sensitive cancer-specific DNA markers. Genomic rearrangements that are present in the vast majority of cancers provide such markers.Tumour DNA isolated from two fresh-frozen lung tumours underwent whole genome sequencing. Genomic rearrangements were detected using a new computational algorithm, GRIDSS. Four genomic rearrangements from each tumour were chosen for further study using rearrangement-specific primers. Six of the eight rearrangements tested were identified as tumour-specific, the remaining two were present in the germline. ctDNA was quantified using digital PCR for the tumour genomic rearrangements in patient blood. Interestingly, one of the patients had no detectable ctDNA either prior to or post surgery although the rearrangements were readily detectable in the tumour DNA.This study demonstrates the feasibility of using digital PCR based on genomic rearrangements for the monitoring of minimal residual disease. In addition, whole genome sequencing provided further information enabling therapeutic choices including the identification of a cryptic EGFR exon 19 deletion in one patient and the identification of a high somatic mutation load in the other patient. This approach can be used as a model for all cancers with rearranged genomes.

Overcoming Clinical Resistance to EZH2 Inhibition Using Rational Epigenetic Combination Therapy.

Epigenetic dependencies have become evident in many cancers. On the basis of 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. On the basis of 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. SIGNIFICANCE: Genomic studies of patient epithelioid sarcomas and rhabdoid tumors identify mutations converging on a common pathway for response to EZH2 inhibition. Resistance mutations decouple drug-induced differentiation from cell-cycle control. We identify an epigenetic combination strategy to overcome resistance and improve durability of response, supporting its investigation in clinical trials. See related commentary by Paolini and Souroullas, p. 903. This article is featured in Selected Articles from This Issue, p. 897.

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.

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.

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.

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.