Long-read genome sequencing resolves a complex 13q structural variant associated with syndromic anophthalmia.

Microphthalmia, anophthalmia, and coloboma (MAC) are a heterogeneous spectrum of anomalous eye development and degeneration with genetic and environmental etiologies. Structural and copy number variants of chromosome 13 have been implicated in MAC; however, the specific loci involved in disease pathogenesis have not been well-defined. Herein we report a newborn with syndromic degenerative anophthalmia and a complex de novo rearrangement of chromosome 13q. Long-read genome sequencing improved the resolution and clinical interpretation of a duplication-triplication/inversion-duplication (DUP-TRP/INV-DUP) and terminal deletion. Sequence features at the breakpoint junctions suggested microhomology-mediated break-induced replication (MMBIR) of the maternal chromosome as the origin. Comparing this rearrangement to previously reported copy number alterations in 13q, we refine a putative dosage-sensitive critical region for MAC that might provide new insights into its molecular etiology.

An infant with congenital respiratory insufficiency and diaphragmatic paralysis: A novel BICD2 phenotype?

Monoallelic pathogenic variants in BICD2 are associated with autosomal dominant Spinal Muscular Atrophy Lower Extremity Predominant 2A and 2B (SMALED2A, SMALED2B). As part of the cellular vesicular transport, complex BICD2 facilitates the flow of constitutive secretory cargoes from the trans-Golgi network, and its dysfunction results in motor neuron loss. The reported phenotypes among patients with SMALED2A and SMALED2B range from a congenital onset disorder of respiratory insufficiency, arthrogryposis, and proximal or distal limb weakness to an adult-onset disorder of limb weakness and contractures. We report an infant with congenital respiratory insufficiency requiring mechanical ventilation, congenital diaphragmatic paralysis, decreased lung volume, and single finger camptodactyly. The infant displayed appropriate antigravity limb movements but had radiological, electrophysiological, and histopathological evidence of myopathy. Exome sequencing and long-read whole-genome sequencing detected a novel de novo BICD2 variant (NM_001003800.1:c.[1543G>A];[=]). This is predicted to encode p.(Glu515Lys); p.Glu515 is located in the coiled-coil 2 mutation hotspot. We hypothesize that this novel phenotype of diaphragmatic paralysis without clear appendicular muscle weakness and contractures of large joints is a presentation of BICD2-related disease.

Comprehensive genomic profiling of glioblastoma tumors, BTICs, and xenografts reveals stability and adaptation to growth environments.

Glioblastoma multiforme (GBM) is the most deadly brain tumor, and currently lacks effective treatment options. Brain tumor-initiating cells (BTICs) and orthotopic xenografts are widely used in investigating GBM biology and new therapies for this aggressive disease. However, the genomic characteristics and molecular resemblance of these models to GBM tumors remain undetermined. We used massively parallel sequencing technology to decode the genomes and transcriptomes of BTICs and xenografts and their matched tumors in order to delineate the potential impacts of the distinct growth environments. Using data generated from whole-genome sequencing of 201 samples and RNA sequencing of 118 samples, we show that BTICs and xenografts resemble their parental tumor at the genomic level but differ at the mRNA expression and epigenomic levels, likely due to the different growth environment for each sample type. These findings suggest that a comprehensive genomic understanding of in vitro and in vivo GBM model systems is crucial for interpreting data from drug screens, and can help control for biases introduced by cell-culture conditions and the microenvironment in mouse models. We also found that lack of MGMT expression in pretreated GBM is linked to hypermutation, which in turn contributes to increased genomic heterogeneity and requires new strategies for GBM treatment.

Slow Off-Rate Modified Aptamer (SOMAmer) Proteomic Analysis of Patient-Derived Malignant Glioma Identifies Distinct Cellular Proteomes.

Malignant gliomas derive from brain glial cells and represent >75% of primary brain tumors. This includes anaplastic astrocytoma (grade III; AS), the most common and fatal glioblastoma multiforme (grade IV; GBM), and oligodendroglioma (ODG). We have generated patient-derived AS, GBM, and ODG cell models to study disease mechanisms and test patient-centered therapeutic strategies. We have used an aptamer-based high-throughput SOMAscan® 1.3K assay to determine the proteomic profiles of 1307 different analytes. SOMAscan® proteomes of AS and GBM self-organized into closely adjacent proteomes which were clearly distinct from ODG proteomes. GBM self-organized into four proteomic clusters of which SOMAscan® cluster 4 proteome predicted a highly inter-connected proteomic network. Several up- and down-regulated proteins relevant to glioma were successfully validated in GBM cell isolates across different SOMAscan® clusters and in corresponding GBM tissues. Slow off-rate modified aptamer proteomics is an attractive analytical tool for rapid proteomic stratification of different malignant gliomas and identified cluster-specific SOMAscan® signatures and functionalities in patient GBM cells.

Improved structural variant interpretation for hereditary cancer susceptibility using long-read sequencing.

PURPOSE: Structural variants (SVs) may be an underestimated cause of hereditary cancer syndromes given the current limitations of short-read next-generation sequencing. Here we investigated the utility of long-read sequencing in resolving germline SVs in cancer susceptibility genes detected through short-read genome sequencing. METHODS: Known or suspected deleterious germline SVs were identified using Illumina genome sequencing across a cohort of 669 advanced cancer patients with paired tumor genome and transcriptome sequencing. Candidate SVs were subsequently assessed by Oxford Nanopore long-read sequencing. RESULTS: Nanopore sequencing confirmed eight simple pathogenic or likely pathogenic SVs, resolving three additional variants whose impact could not be fully elucidated through short-read sequencing. A recurrent sequencing artifact on chromosome 16p13 and one complex rearrangement on chromosome 5q35 were subsequently classified as likely benign, obviating the need for further clinical assessment. Variant configuration was further resolved in one case with a complex pathogenic rearrangement affecting TSC2. CONCLUSION: Our findings demonstrate that long-read sequencing can improve the validation, resolution, and classification of germline SVs. This has important implications for return of results, cascade carrier testing, cancer screening, and prophylactic interventions.

Hypermutation signature reveals a slippage and realignment model of translesion synthesis by Rev3 polymerase in cisplatin-treated yeast.

Gene-gene or gene-drug interactions are typically quantified using fitness as a readout because the data are continuous and easily measured in high throughput. However, to what extent fitness captures the range of other phenotypes that show synergistic effects is usually unknown. Using Saccharomyces cerevisiae and focusing on a matrix of DNA repair mutants and genotoxic drugs, we quantify 76 gene-drug interactions based on both mutation rate and fitness and find that these parameters are not connected. Independent of fitness defects, we identified six cases of synthetic hypermutation, where the combined effect of the drug and mutant on mutation rate was greater than predicted. One example occurred when yeast lacking RAD1 were exposed to cisplatin, and we characterized this interaction using whole-genome sequencing. Our sequencing results indicate mutagenesis by cisplatin in rad1Δ cells appeared to depend almost entirely on interstrand cross-links at GpCpN motifs. Interestingly, our data suggest that the following base on the template strand dictates the addition of the mutated base. This result differs from cisplatin mutation signatures in XPF-deficient Caenorhabditis elegans and supports a model in which translesion synthesis polymerases perform a slippage and realignment extension across from the damaged base. Accordingly, DNA polymerase ζ activity was essential for mutagenesis in cisplatin-treated rad1Δ cells. Together these data reveal the potential to gain new mechanistic insights from nonfitness measures of gene-drug interactions and extend the use of mutation accumulation and whole-genome sequencing analysis to define DNA repair mechanisms.

Clonal expansion and epigenetic reprogramming following deletion or amplification of mutant IDH1.

IDH1 mutation is the earliest genetic alteration in low-grade gliomas (LGGs), but its role in tumor recurrence is unclear. Mutant IDH1 drives overproduction of the oncometabolite d-2-hydroxyglutarate (2HG) and a CpG island (CGI) hypermethylation phenotype (G-CIMP). To investigate the role of mutant IDH1 at recurrence, we performed a longitudinal analysis of 50 IDH1 mutant LGGs. We discovered six cases with copy number alterations (CNAs) at the IDH1 locus at recurrence. Deletion or amplification of IDH1 was followed by clonal expansion and recurrence at a higher grade. Successful cultures derived from IDH1 mutant, but not IDH1 wild type, gliomas systematically deleted IDH1 in vitro and in vivo, further suggestive of selection against the heterozygous mutant state as tumors progress. Tumors and cultures with IDH1 CNA had decreased 2HG, maintenance of G-CIMP, and DNA methylation reprogramming outside CGI. Thus, while IDH1 mutation initiates gliomagenesis, in some patients mutant IDH1 and 2HG are not required for later clonal expansions.

FORGE Canada Consortium: outcomes of a 2-year national rare-disease gene-discovery project.

Inherited monogenic disease has an enormous impact on the well-being of children and their families. Over half of the children living with one of these conditions are without a molecular diagnosis because of the rarity of the disease, the marked clinical heterogeneity, and the reality that there are thousands of rare diseases for which causative mutations have yet to be identified. It is in this context that in 2010 a Canadian consortium was formed to rapidly identify mutations causing a wide spectrum of pediatric-onset rare diseases by using whole-exome sequencing. The FORGE (Finding of Rare Disease Genes) Canada Consortium brought together clinicians and scientists from 21 genetics centers and three science and technology innovation centers from across Canada. From nation-wide requests for proposals, 264 disorders were selected for study from the 371 submitted; disease-causing variants (including in 67 genes not previously associated with human disease; 41 of these have been genetically or functionally validated, and 26 are currently under study) were identified for 146 disorders over a 2-year period. Here, we present our experience with four strategies employed for gene discovery and discuss FORGE's impact in a number of realms, from clinical diagnostics to the broadening of the phenotypic spectrum of many diseases to the biological insight gained into both disease states and normal human development. Lastly, on the basis of this experience, we discuss the way forward for rare-disease genetic discovery both in Canada and internationally.

Compound heterozygous TRPV4 mutations in two siblings with a complex phenotype including severe intellectual disability and neuropathy.

TRPV4 encodes a polymodal calcium-permeable plasma membrane channel. Dominant pathogenic mutations in TRPV4 lead to a wide spectrum of abnormal phenotypes. This is the first report of biallelic TRPV4 mutations and we describe two compound heterozygous siblings presenting with a complex phenotype including severe neuromuscular involvement. In light of previously well described dominant inheritance for TRPV4-related neuromuscular disease, our study suggests a role for compound heterozygosity and loss-of-function as a potential novel disease mechanism for this group of disorders. Profound intellectual disability was also noted in both affected children, suggesting that TRPV4 may be necessary for normal brain development.

A novel mutation in EED associated with overgrowth.

In a patient suspected clinically to have Weaver syndrome, we ruled out mutations in EZH2 and NSD1, then identified a previously undescribed de novo mutation in EZH2's partner protein EED. Both proteins are members of the Polycomb Repressive Complex 2 that maintains gene silencing. On the basis of the similarities of the patient's phenotype to Weaver syndrome, which is caused by de novo mutations in EZH2, and on other lines of evidence including mouse Eed hypomorphs, we characterize this mutation as probably pathogenic for a Weaver-like overgrowth syndrome. This is the first report of overgrowth and related phenotypes associated with a constitutional mutation in human EED.

Diagnostic value of next-generation sequencing in an unusual sphenoid tumor.

Extraordinary advancements in sequencing technology have made what was once a decade-long multi-institutional endeavor into a methodology with the potential for practical use in a clinical setting. We therefore set out to examine the clinical value of next-generation sequencing by enrolling patients with incurable or ambiguous tumors into the Personalized OncoGenomics initiative at the British Columbia Cancer Agency whereby whole genome and transcriptome analyses of tumor/normal tissue pairs are completed with the ultimate goal of directing therapeutics. First, we established that the sequencing, analysis, and communication with oncologists could be completed in less than 5 weeks. Second, we found that cancer diagnostics is an area that can greatly benefit from the comprehensiveness of a whole genome analysis. Here, we present a scenario in which a metastasized sphenoid mass, which was initially thought of as an undifferentiated squamous cell carcinoma, was rediagnosed as an SMARCB1-negative rhabdoid tumor based on the newly acquired finding of homozygous SMARCB1 deletion. The new diagnosis led to a change in chemotherapy and a complete nodal response in the patient. This study also provides additional insight into the mutational landscape of an adult SMARCB1-negative tumor that has not been explored at a whole genome and transcriptome level.

Mitofusin 2 Variant Presenting With a Phenotype of Multiple System Atrophy of Cerebellar Subtype.

Objectives: To investigate the etiology of cerebellar ataxia in an adult male patient. Methods: We performed standard neurologic assessment and genome sequencing of a 62-year-old man with rapidly progressive balance and gait abnormalities. Results: The propositus exhibited cognitive dysfunction, mild appendicular bradykinesia, prominent appendicular ataxia, dysarthria, and hypomimia with minimal dysautonomic symptoms. Nerve conduction studies showed mild peripheral sensory neuropathy and normal motor nerve conduction velocities. Brain imaging showed progressive cerebellar atrophy and gliosis of the olivopontocerebellar fibers, characterized by T2 hyperintensity within the pons. Genetic testing revealed a likely pathogenic germline variant in MFN2 (NM_014874: c.[838C>T];[=], p.(R280C)) in the GTPase domain (G) interface; pathogenic variants of MFN2 typically cause hereditary sensory and motor neuropathy VI or Charcot-Marie-Tooth disease 2A. The presence of progressive ataxia, "hot cross bun" sign, and dysautonomia has been associated with multiple system atrophy, cerebellar type (MSA-C). Discussion: We describe progressive cerebellar ataxia in an individual with a deleterious variant in MFN2. Our findings suggest that pathogenic variants in MFN2 can result in a spectrum of phenotypes including cerebellar ataxia with cerebellar-pontine atrophy in the absence of significant neuropathy and in a manner closely resembling MSA-C.

Whole genome and transcriptome integrated analyses guide clinical care of pediatric poor prognosis cancers.

The role for routine whole genome and transcriptome analysis (WGTA) for poor prognosis pediatric cancers remains undetermined. Here, we characterize somatic mutations, structural rearrangements, copy number variants, gene expression, immuno-profiles and germline cancer predisposition variants in children and adolescents with relapsed, refractory or poor prognosis malignancies who underwent somatic WGTA and matched germline sequencing. Seventy-nine participants with a median age at enrollment of 8.8 y (range 6 months to 21.2 y) are included. Germline pathogenic/likely pathogenic variants are identified in 12% of participants, of which 60% were not known prior. Therapeutically actionable variants are identified by targeted gene report and whole genome in 32% and 62% of participants, respectively, and increase to 96% after integrating transcriptome analyses. Thirty-two molecularly informed therapies are pursued in 28 participants with 54% achieving a clinical benefit rate; objective response or stable disease ≥6 months. Integrated WGTA identifies therapeutically actionable variants in almost all tumors and are directly translatable to clinical care of children with poor prognosis cancers.

Ortholog identification in the presence of domain architecture rearrangement.

Ortholog identification is used in gene functional annotation, species phylogeny estimation, phylogenetic profile construction and many other analyses. Bioinformatics methods for ortholog identification are commonly based on pairwise protein sequence comparisons between whole genomes. Phylogenetic methods of ortholog identification have also been developed; these methods can be applied to protein data sets sharing a common domain architecture or which share a single functional domain but differ outside this region of homology. While promiscuous domains represent a challenge to all orthology prediction methods, overall structural similarity is highly correlated with proximity in a phylogenetic tree, conferring a degree of robustness to phylogenetic methods. In this article, we review the issues involved in orthology prediction when data sets include sequences with structurally heterogeneous domain architectures, with particular attention to automated methods designed for high-throughput application, and present a case study to illustrate the challenges in this area.

Defining the heterogeneity of unbalanced structural variation underlying breast cancer susceptibility by nanopore genome sequencing.

Germline structural variants (SVs) are challenging to resolve by conventional genetic testing assays. Long-read sequencing has improved the global characterization of SVs, but its sensitivity at cancer susceptibility loci has not been reported. Nanopore long-read genome sequencing was performed for nineteen individuals with pathogenic copy number alterations in BRCA1, BRCA2, CHEK2 and PALB2 identified by prior clinical testing. Fourteen variants, which spanned single exons to whole genes and included a tandem duplication, were accurately represented. Defining the precise breakpoints of SVs in BRCA1 and CHEK2 revealed unforeseen allelic heterogeneity and informed the mechanisms underlying the formation of recurrent deletions. Integrating read-based and statistical phasing further helped define extended haplotypes associated with founder alleles. Long-read sequencing is a sensitive method for characterizing private, recurrent and founder SVs underlying breast cancer susceptibility. Our findings demonstrate the potential for nanopore sequencing as a powerful genetic testing assay in the hereditary cancer setting.

macroH2A2 antagonizes epigenetic programs of stemness in glioblastoma.

Self-renewal is a crucial property of glioblastoma cells that is enabled by the choreographed functions of chromatin regulators and transcription factors. Identifying targetable epigenetic mechanisms of self-renewal could therefore represent an important step toward developing effective treatments for this universally lethal cancer. Here we uncover an epigenetic axis of self-renewal mediated by the histone variant macroH2A2. With omics and functional assays deploying patient-derived in vitro and in vivo models, we show that macroH2A2 shapes chromatin accessibility at enhancer elements to antagonize transcriptional programs of self-renewal. macroH2A2 also sensitizes cells to small molecule-mediated cell death via activation of a viral mimicry response. Consistent with these results, our analyses of clinical cohorts indicate that high transcriptional levels of this histone variant are associated with better prognosis of high-grade glioma patients. Our results reveal a targetable epigenetic mechanism of self-renewal controlled by macroH2A2 and suggest additional treatment approaches for glioblastoma patients.

Distribution and properties of the genes encoding the biosynthesis of the bacterial cofactor, pyrroloquinoline quinone.

Pyrroloquinoline quinone (PQQ) is a small, redox active molecule that serves as a cofactor for several bacterial dehydrogenases, introducing pathways for carbon utilization that confer a growth advantage. Early studies had implicated a ribosomally translated peptide as the substrate for PQQ production. This study presents a sequence- and structure-based analysis of the components of the pqq operon. We find the necessary components for PQQ production are present in 126 prokaryotes, most of which are Gram-negative and a number of which are pathogens. A total of five gene products, PqqA, PqqB, PqqC, PqqD, and PqqE, are identified as being obligatory for PQQ production. Three of the gene products in the pqq operon, PqqB, PqqC, and PqqE, are members of large protein superfamilies. By combining evolutionary conservation patterns with information from three-dimensional structures, we are able to differentiate the gene products involved in PQQ biosynthesis from those with divergent functions. The observed persistence of a conserved gene order within analyzed operons strongly suggests a role for protein-protein interactions in the course of cofactor biosynthesis. These studies propose previously unidentified roles for several of the gene products, as well as identifying possible new targets for antibiotic design and application.

Genomic analysis of the basal lineage fungus Rhizopus oryzae reveals a whole-genome duplication.

Rhizopus oryzae is the primary cause of mucormycosis, an emerging, life-threatening infection characterized by rapid angioinvasive growth with an overall mortality rate that exceeds 50%. As a representative of the paraphyletic basal group of the fungal kingdom called "zygomycetes," R. oryzae is also used as a model to study fungal evolution. Here we report the genome sequence of R. oryzae strain 99-880, isolated from a fatal case of mucormycosis. The highly repetitive 45.3 Mb genome assembly contains abundant transposable elements (TEs), comprising approximately 20% of the genome. We predicted 13,895 protein-coding genes not overlapping TEs, many of which are paralogous gene pairs. The order and genomic arrangement of the duplicated gene pairs and their common phylogenetic origin provide evidence for an ancestral whole-genome duplication (WGD) event. The WGD resulted in the duplication of nearly all subunits of the protein complexes associated with respiratory electron transport chains, the V-ATPase, and the ubiquitin-proteasome systems. The WGD, together with recent gene duplications, resulted in the expansion of multiple gene families related to cell growth and signal transduction, as well as secreted aspartic protease and subtilase protein families, which are known fungal virulence factors. The duplication of the ergosterol biosynthetic pathway, especially the major azole target, lanosterol 14alpha-demethylase (ERG11), could contribute to the variable responses of R. oryzae to different azole drugs, including voriconazole and posaconazole. Expanded families of cell-wall synthesis enzymes, essential for fungal cell integrity but absent in mammalian hosts, reveal potential targets for novel and R. oryzae-specific diagnostic and therapeutic treatments.