Detection of Genomic Structural Variations Associated with Drug Sensitivity and Resistance in Acute Leukemia.

Acute leukemia is a particularly problematic collection of hematological cancers, and, while somewhat rare, the survival rate of patients is typically abysmal without bone marrow transplantation. Furthermore, traditional chemotherapies used as standard-of-care for patients cause significant side effects. Understanding the evolution of leukemia to identify novel targets and, therefore, drug treatment regimens is a significant medical need. Genomic rearrangements and other structural variations (SVs) have long been known to be causative and pathogenic in multiple types of cancer, including leukemia. These SVs may be involved in cancer initiation, progression, clonal evolution, and drug resistance, and a better understanding of SVs from individual patients may help guide therapeutic options. Here, we show the utilization of optical genome mapping (OGM) to detect known and novel SVs in the samples of patients with leukemia. Importantly, this technology provides an unprecedented level of granularity and quantitation unavailable to other current techniques and allows for the unbiased detection of novel SVs, which may be relevant to disease pathogenesis and/or drug resistance. Coupled with the chemosensitivities of these samples to FDA-approved oncology drugs, we show how an impartial integrative analysis of these diverse datasets can be used to associate the detected genomic rearrangements with multiple drug sensitivity profiles. Indeed, an insertion in the gene MUSK is shown to be associated with increased sensitivity to the clinically relevant agent Idarubicin, while partial tandem duplication events in the KMT2A gene are related to the efficacy of another frontline treatment, Cytarabine.

Detection of mosaic and population-level structural variants with Sniffles2.

Calling structural variations (SVs) is technically challenging, but using long reads remains the most accurate way to identify complex genomic alterations. Here we present Sniffles2, which improves over current methods by implementing a repeat aware clustering coupled with a fast consensus sequence and coverage-adaptive filtering. Sniffles2 is 11.8 times faster and 29% more accurate than state-of-the-art SV callers across different coverages (5-50×), sequencing technologies (ONT and HiFi) and SV types. Furthermore, Sniffles2 solves the problem of family-level to population-level SV calling to produce fully genotyped VCF files. Across 11 probands, we accurately identified causative SVs around MECP2, including highly complex alleles with three overlapping SVs. Sniffles2 also enables the detection of mosaic SVs in bulk long-read data. As a result, we identified multiple mosaic SVs in brain tissue from a patient with multiple system atrophy. The identified SV showed a remarkable diversity within the cingulate cortex, impacting both genes involved in neuron function and repetitive elements.

Structural variant analysis of a cancer reference cell line sample using multiple sequencing technologies.

BACKGROUND: The cancer genome is commonly altered with thousands of structural rearrangements including insertions, deletions, translocation, inversions, duplications, and copy number variations. Thus, structural variant (SV) characterization plays a paramount role in cancer target identification, oncology diagnostics, and personalized medicine. As part of the SEQC2 Consortium effort, the present study established and evaluated a consensus SV call set using a breast cancer reference cell line and matched normal control derived from the same donor, which were used in our companion benchmarking studies as reference samples. RESULTS: We systematically investigated somatic SVs in the reference cancer cell line by comparing to a matched normal cell line using multiple NGS platforms including Illumina short-read, 10X Genomics linked reads, PacBio long reads, Oxford Nanopore long reads, and high-throughput chromosome conformation capture (Hi-C). We established a consensus SV call set of a total of 1788 SVs including 717 deletions, 230 duplications, 551 insertions, 133 inversions, 146 translocations, and 11 breakends for the reference cancer cell line. To independently evaluate and cross-validate the accuracy of our consensus SV call set, we used orthogonal methods including PCR-based validation, Affymetrix arrays, Bionano optical mapping, and identification of fusion genes detected from RNA-seq. We evaluated the strengths and weaknesses of each NGS technology for SV determination, and our findings provide an actionable guide to improve cancer genome SV detection sensitivity and accuracy. CONCLUSIONS: A high-confidence consensus SV call set was established for the reference cancer cell line. A large subset of the variants identified was validated by multiple orthogonal methods.

ATTCT and ATTCC repeat expansions in the ATXN10 gene affect disease penetrance of spinocerebellar ataxia type 10.

Spinocerebellar ataxia type 10 (SCA10) is an autosomal-dominant disorder caused by an expanded pentanucleotide repeat in the ATXN10 gene. This repeat expansion, when fully penetrant, has a size of 850-4,500 repeats. It has been shown that the repeat composition can be a modifier of disease, e.g., seizures. Here, we describe a Mexican kindred in which we identified both pure (ATTCT)n and mixed (ATTCT)n-(ATTCC)n expansions in the same family. We used amplification-free targeted sequencing and optical genome mapping to decipher the composition of these repeat expansions. We found a considerable degree of mosaicism of the repeat expansion. This mosaicism was confirmed in skin fibroblasts from individuals with ATXN10 expansions with RNAScope in situ hybridization. All affected family members with the mixed ATXN10 repeat expansion showed typical clinical signs of spinocerebellar ataxia and epilepsy. In contrast, individuals with the pure ATXN10 expansion present with Parkinson's disease or are unaffected, even in individuals more than 20 years older than the average age at onset for SCA10. Our findings suggest that the pure (ATTCT)n expansion is non-pathogenic, while repeat interruptions, e.g., (ATTCC)n, are necessary to cause SCA10. This mechanism has been recently described for several other repeat expansions including SCA31 (BEAN1), SCA37 (DAB1), and three loci for benign adult familial myoclonic epilepsy BAFME (SAMD12, TNRC6A, RAPGEF2). Therefore, long-read sequencing and optical genome mapping of the entire genomic structure of repeat expansions are critical for clinical practice and genetic counseling, as variations in the repeat can affect disease penetrance, symptoms, and disease trajectory.

Screening and characterization of polyhydroxyalkanoate granules, and phylogenetic analysis of polyhydroxyalkanoate synthase gene PhaC in cyanobacteria.

Using Nile Red and BODIPY 493/503 dye-staining and fluorescence microscopy, twenty cyanobacterial strains, including ten commercially available strains and ten environmental isolates from estuaries, freshwater ponds, and lagoons, were screened for the accumulation of ecologically important and potentially biotechnologically significant carbon storage granules such as polyhydroxyalkanoates (PHA). Dye-staining granules were observed in six strains. Three Synechocystis, spp. strains WHSYN, LSNM, and CGF-1, and a Phormidium-like sp. CGFILA were isolated from environmental sources and found to produce granules of polyhydroxyalkanoate (PHA) according to PHA synthase gene (phaC) PCR screening and 1 H NMR analyses. The environmental isolate, Nodularia sp. Las Olas and commercially available Phormidium cf. iriguum CCALA 759 displayed granules but screened negative for PHA according to phaC PCR and 1 H NMR analyses. Partial polyhydroxyalkanoate synthase subunit C (phaC) and 16S rRNA gene sequences obtained from the PHA-accumulating strains and analyzed alongside publicly available phaC, phaE, 16S rRNA, and 23S rRNA data help in understanding the distribution and evolutionary history of PHA biosynthesis within the phylum Cyanobacteria. The data show that the presence of phaC is highly conserved within the genus Synechocystis, and present in at least one isolate of Phormidium. Maximum likelihood analyses and cophylogenetic modeling of PHA synthase gene sequences provide evidence of a recent horizontal gene transfer event between distant genera of cyanobacteria related to Pleurocapsa sp. PCC 7327 and Phormidium-like sp. CGFILA. These findings will help guide additional screening for PHA producers, and may explain why some Phormidium species produce PHAs, while others do not.

Unbiased optical mapping of telomere-integrated endogenous human herpesvirus 6.

Next-generation sequencing technologies allowed sequencing of thousands of genomes. However, there are genomic regions that remain difficult to characterize, including telomeres, centromeres, and other low-complexity regions, as well as transposable elements and endogenous viruses. Human herpesvirus 6A and 6B (HHV-6A and HHV-6B) are closely related viruses that infect most humans and can integrate their genomes into the telomeres of infected cells. Integration also occurs in germ cells, meaning that the virus can be inherited and result in individuals harboring the virus in every cell of their body. The integrated virus can reactivate and cause disease in humans. While it is well established that the virus resides in the telomere region, the integration locus is poorly defined due to the low sequence complexity (TTAGGG)n of telomeres that cannot be easily resolved through sequencing. We therefore employed genome imaging of the integrated HHV-6A and HHV-6B genomes using whole-genome optical site mapping technology. Using this technology, we identified which chromosome arm harbors the virus genome and obtained a high-resolution map of the integration loci of multiple patients. Surprisingly, this revealed long telomere sequences at the virus-subtelomere junction that were previously missed using PCR-based approaches. Contrary to what was previously thought, our technique revealed that the telomere lengths of chromosomes harboring the integrated virus genome were comparable to the other chromosomes. Taken together, our data shed light on the genetic structure of the HHV-6A and HHV-6B integration locus, demonstrating the utility of optical mapping for the analysis of genomic regions that are difficult to sequence.