Interfaces of Malignant and Immunologic Clonal Dynamics in Ovarian Cancer.

High-grade serous ovarian cancer (HGSC) exhibits extensive malignant clonal diversity with widespread but non-random patterns of disease dissemination. We investigated whether local immune microenvironment factors shape tumor progression properties at the interface of tumor-infiltrating lymphocytes (TILs) and cancer cells. Through multi-region study of 212 samples from 38 patients with whole-genome sequencing, immunohistochemistry, histologic image analysis, gene expression profiling, and T and B cell receptor sequencing, we identified three immunologic subtypes across samples and extensive within-patient diversity. Epithelial CD8+ TILs negatively associated with malignant diversity, reflecting immunological pruning of tumor clones inferred by neoantigen depletion, HLA I loss of heterozygosity, and spatial tracking between T cell and tumor clones. In addition, combinatorial prognostic effects of mutational processes and immune properties were observed, illuminating how specific genomic aberration types associate with immune response and impact survival. We conclude that within-patient spatial immune microenvironment variation shapes intraperitoneal malignant spread, provoking new evolutionary perspectives on HGSC clonal dispersion.

Obtaining genetics insights from deep learning via explainable artificial intelligence.

Artificial intelligence (AI) models based on deep learning now represent the state of the art for making functional predictions in genomics research. However, the underlying basis on which predictive models make such predictions is often unknown. For genomics researchers, this missing explanatory information would frequently be of greater value than the predictions themselves, as it can enable new insights into genetic processes. We review progress in the emerging area of explainable AI (xAI), a field with the potential to empower life science researchers to gain mechanistic insights into complex deep learning models. We discuss and categorize approaches for model interpretation, including an intuitive understanding of how each approach works and their underlying assumptions and limitations in the context of typical high-throughput biological datasets.

JASPAR 2024: 20th anniversary of the open-access database of transcription factor binding profiles.

JASPAR (https://jaspar.elixir.no/) is a widely-used open-access database presenting manually curated high-quality and non-redundant DNA-binding profiles for transcription factors (TFs) across taxa. In this 10th release and 20th-anniversary update, the CORE collection has expanded with 329 new profiles. We updated three existing profiles and provided orthogonal support for 72 profiles from the previous release's UNVALIDATED collection. Altogether, the JASPAR 2024 update provides a 20% increase in CORE profiles from the previous release. A trimming algorithm enhanced profiles by removing low information content flanking base pairs, which were likely uninformative (within the capacity of the PFM models) for TFBS predictions and modelling TF-DNA interactions. This release includes enhanced metadata, featuring a refined classification for plant TFs' structural DNA-binding domains. The new JASPAR collections prompt updates to the genomic tracks of predicted TF binding sites (TFBSs) in 8 organisms, with human and mouse tracks available as native tracks in the UCSC Genome browser. All data are available through the JASPAR web interface and programmatically through its API and the updated Bioconductor and pyJASPAR packages. Finally, a new TFBS extraction tool enables users to retrieve predicted JASPAR TFBSs intersecting their genomic regions of interest.

OnTarget: in silico design of MiniPromoters for targeted delivery of expression.

MiniPromoters, or compact promoters, are short DNA sequences that can drive expression in specific cells and tissues. While broadly useful, they are of high relevance to gene therapy due to their role in enabling precise control of where a therapeutic gene will be expressed. Here, we present OnTarget (http://ontarget.cmmt.ubc.ca), a webserver that streamlines the MiniPromoter design process. Users only need to specify a gene of interest or custom genomic coordinates on which to focus the identification of promoters and enhancers, and can also provide relevant cell-type-specific genomic evidence (e.g. accessible chromatin regions, histone modifications, etc.). OnTarget combines the provided data with internal data to identify candidate promoters and enhancers and design MiniPromoters. To illustrate the utility of OnTarget, we designed and characterized two MiniPromoters targeting different cell populations relevant to Parkinson Disease.

The genome atlas: navigating a new era of reference genomes.

The reference genome serves two distinct purposes within the field of genomics. First, it provides a persistent structure against which findings can be reported, allowing for universal knowledge exchange between users. Second, it reduces the computational costs and time required to process genomic data by creating a scaffold that can be relied upon by analysis software. Here, we posit that current efforts to extend the linear reference to a graph-based structure while trying to fulfil both of these purposes concurrently will face a trade-off between comprehensiveness and computational efficiency. In this article, we explore how the reference genome is used and suggest an alternative structure, The Genome Atlas (TGA), to fulfil the bipartite role of the reference genome.

CIAO1 and MMS19 deficiency: A lethal neurodegenerative phenotype caused by cytosolic Fe-S cluster protein assembly disorders.

PURPOSE: The functionality of many cellular proteins depends on cofactors; yet, they have only been implicated in a minority of Mendelian diseases. Here, we describe the first 2 inherited disorders of the cytosolic iron-sulfur protein assembly system. METHODS: Genetic testing via genome sequencing was applied to identify the underlying disease cause in 3 patients with microcephaly, congenital brain malformations, progressive developmental and neurologic impairments, recurrent infections, and a fatal outcome. Studies in patient-derived skin fibroblasts and zebrafish models were performed to investigate the biochemical and cellular consequences. RESULTS: Metabolic analysis showed elevated uracil and thymine levels in body fluids but no pathogenic variants in DPYD, encoding dihydropyrimidine dehydrogenase. Genome sequencing identified compound heterozygosity in 2 patients for missense variants in CIAO1, encoding cytosolic iron-sulfur assembly component 1, and homozygosity for an in-frame 3-nucleotide deletion in MMS19, encoding the MMS19 homolog, cytosolic iron-sulfur assembly component, in the third patient. Profound alterations in the proteome, metabolome, and lipidome were observed in patient-derived fibroblasts. We confirmed the detrimental effect of deficiencies in CIAO1 and MMS19 in zebrafish models. CONCLUSION: A general failure of cytosolic and nuclear iron-sulfur protein maturation caused pleiotropic effects. The critical function of the cytosolic iron-sulfur protein assembly machinery for antiviral host defense may well explain the recurrent severe infections occurring in our patients.

Human complete NFAT1 deficiency causes a triad of joint contractures, osteochondromas, and B-cell malignancy.

The discovery of humans with monogenic disorders has a rich history of generating new insights into biology. Here we report the first human identified with complete deficiency of nuclear factor of activated T cells 1 (NFAT1). NFAT1, encoded by NFATC2, mediates calcium-calcineurin signals that drive cell activation, proliferation, and survival. The patient is homozygous for a damaging germline NFATC2 variant (c.2023_2026delTACC; p.Tyr675Thrfs∗18) and presented with joint contractures, osteochondromas, and recurrent B-cell lymphoma. Absence of NFAT1 protein in chondrocytes caused enrichment in prosurvival and inflammatory genes. Systematic single-cell-omic analyses in PBMCs revealed an environment that promotes lymphomagenesis with accumulation of naïve B cells (enriched for oncogenic signatures MYC and JAK1), exhausted CD4+ T cells, impaired T follicular helper cells, and aberrant CD8+ T cells. This work highlights the pleiotropic role of human NFAT1, will empower the diagnosis of additional patients with NFAT1 deficiency, and further defines the detrimental effects associated with long-term use of calcineurin inhibitors.

JASPAR 2022: the 9th release of the open-access database of transcription factor binding profiles.

JASPAR (http://jaspar.genereg.net/) is an open-access database containing manually curated, non-redundant transcription factor (TF) binding profiles for TFs across six taxonomic groups. In this 9th release, we expanded the CORE collection with 341 new profiles (148 for plants, 101 for vertebrates, 85 for urochordates, and 7 for insects), which corresponds to a 19% expansion over the previous release. We added 298 new profiles to the Unvalidated collection when no orthogonal evidence was found in the literature. All the profiles were clustered to provide familial binding profiles for each taxonomic group. Moreover, we revised the structural classification of DNA binding domains to consider plant-specific TFs. This release introduces word clouds to represent the scientific knowledge associated with each TF. We updated the genome tracks of TFBSs predicted with JASPAR profiles in eight organisms; the human and mouse TFBS predictions can be visualized as native tracks in the UCSC Genome Browser. Finally, we provide a new tool to perform JASPAR TFBS enrichment analysis in user-provided genomic regions. All the data is accessible through the JASPAR website, its associated RESTful API, the R/Bioconductor data package, and a new Python package, pyJASPAR, that facilitates serverless access to the data.

Deregulated Regulators: Disease-Causing cis Variants in Transcription Factor Genes.

Whole-genome sequencing is accelerating identification of noncoding variants that disrupt gene expression, although reports of such regulatory variants implicated in disease remain rare. A notable subset of described variants affect transcription factor (TF) genes and other master regulators in cis through dosage effects. From the literature, we compiled 46 regulatory variants linked to 40 TF genes implicated in rare diseases. We discuss the genomic geography of these variants and the evidence presented for their potential pathogenicity. To help advance research on candidate disease variants into the literature, we introduce an evidence framework specific to regulatory variants, which are under-represented in current variant classification guidelines. The clinical research interpretation of patient genomes may be advanced by considering regulatory variants, particularly those that deregulate TF genes.

RevUP: an online scoring system for regulatory variants implicated in rare diseases.

SUMMARY: To address the difficulty in assessing the implication of regulatory variants in diseases, a scoring scheme previously published allows the calculation of the Regulatory Variant Evidence score (RVE-score). The score represents the accumulated evidence for a causative role of a regulatory variant in a disease. Regulatory Evidence for Variants Underlying Phenotypes was built to calculate the RVE-score of regulatory variants, based on the 24 criteria, with a hybrid approach combining information retrieved from public databases and user input. AVAILABILITY AND IMPLEMENTATION: RevUP is freely available at http://www.revup-classifier.ca. The source code is available at https://github.com/wassermanlab/revup. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Targeting AXL kinase sensitizes leukemic stem and progenitor cells to venetoclax treatment in acute myeloid leukemia.

The abundance of genetic abnormalities and phenotypic heterogeneities in acute myeloid leukemia (AML) poses significant challenges to the development of improved treatments. Here, we demonstrated that a key growth arrest-specific gene 6/AXL axis is highly activated in cells from patients with AML, particularly in stem/progenitor cells. We developed a potent selective AXL inhibitor that has favorable pharmaceutical properties and efficacy against preclinical patient-derived xenotransplantation (PDX) models of AML. Importantly, inhibition of AXL sensitized AML stem/progenitor cells to venetoclax treatment, with strong synergistic effects in vitro and in PDX models. Mechanistically, single-cell RNA-sequencing and functional validation studies uncovered that AXL inhibition, alone or in combination with venetoclax, potentially targets intrinsic metabolic vulnerabilities of AML stem/progenitor cells and shows a distinct transcriptomic profile and inhibits mitochondrial oxidative phosphorylation. Inhibition of AXL or BCL-2 also differentially targets key signaling proteins to synergize in leukemic cell killing. These findings have a direct translational impact on the treatment of AML and other cancers with high AXL activity.

Linkage analysis identifies an isolated strabismus locus at 14q12 overlapping with FOXG1 syndrome region.

Strabismus is a common condition, affecting 1%-4% of individuals. Isolated strabismus has been studied in families with Mendelian inheritance patterns. Despite the identification of multiple loci via linkage analyses, no specific genes have been identified from these studies. The current study is based on a seven-generation family with isolated strabismus inherited in an autosomal dominant manner. A total of 13 individuals from a common ancestor have been included for linkage analysis. Among these, nine are affected and four are unaffected. A single linkage signal has been identified at an 8.5 Mb region of chromosome 14q12 with a multipoint LOD (logarithm of the odds) score of 4.69. Disruption of this locus is known to cause FOXG1 syndrome (or congenital Rett syndrome; OMIM #613454 and *164874), in which 84% of affected individuals present with strabismus. With the incorporation of next-generation sequencing and in-depth bioinformatic analyses, a 4 bp non-coding deletion was prioritised as the top candidate for the observed strabismus phenotype. The deletion is predicted to disrupt regulation of FOXG1, which encodes a transcription factor of the Forkhead family. Suggestive of an autoregulation effect, the disrupted sequence matches the consensus FOXG1 and Forkhead family transcription factor binding site and has been observed in previous ChIP-seq studies to be bound by Foxg1 in early mouse brain development. Future study of this specific deletion may shed light on the regulation of FOXG1 expression and may enhance our understanding of the mechanisms contributing to strabismus and FOXG1 syndrome.

Bi-allelic GOT2 Mutations Cause a Treatable Malate-Aspartate Shuttle-Related Encephalopathy.

Early-infantile encephalopathies with epilepsy are devastating conditions mandating an accurate diagnosis to guide proper management. Whole-exome sequencing was used to investigate the disease etiology in four children from independent families with intellectual disability and epilepsy, revealing bi-allelic GOT2 mutations. In-depth metabolic studies in individual 1 showed low plasma serine, hypercitrullinemia, hyperlactatemia, and hyperammonemia. The epilepsy was serine and pyridoxine responsive. Functional consequences of observed mutations were tested by measuring enzyme activity and by cell and animal models. Zebrafish and mouse models were used to validate brain developmental and functional defects and to test therapeutic strategies. GOT2 encodes the mitochondrial glutamate oxaloacetate transaminase. GOT2 enzyme activity was deficient in fibroblasts with bi-allelic mutations. GOT2, a member of the malate-aspartate shuttle, plays an essential role in the intracellular NAD(H) redox balance. De novo serine biosynthesis was impaired in fibroblasts with GOT2 mutations and GOT2-knockout HEK293 cells. Correcting the highly oxidized cytosolic NAD-redox state by pyruvate supplementation restored serine biosynthesis in GOT2-deficient cells. Knockdown of got2a in zebrafish resulted in a brain developmental defect associated with seizure-like electroencephalography spikes, which could be rescued by supplying pyridoxine in embryo water. Both pyridoxine and serine synergistically rescued embryonic developmental defects in zebrafish got2a morphants. The two treated individuals reacted favorably to their treatment. Our data provide a mechanistic basis for the biochemical abnormalities in GOT2 deficiency that may also hold for other MAS defects.

Glutaminase Deficiency Caused by Short Tandem Repeat Expansion in GLS.

We report an inborn error of metabolism caused by an expansion of a GCA-repeat tract in the 5' untranslated region of the gene encoding glutaminase (GLS) that was identified through detailed clinical and biochemical phenotyping, combined with whole-genome sequencing. The expansion was observed in three unrelated patients who presented with an early-onset delay in overall development, progressive ataxia, and elevated levels of glutamine. In addition to ataxia, one patient also showed cerebellar atrophy. The expansion was associated with a relative deficiency of GLS messenger RNA transcribed from the expanded allele, which probably resulted from repeat-mediated chromatin changes upstream of the GLS repeat. Our discovery underscores the importance of careful examination of regions of the genome that are typically excluded from or poorly captured by exome sequencing.

Demonstrating the utility of flexible sequence queries against indexed short reads with FlexTyper.

Across the life sciences, processing next generation sequencing data commonly relies upon a computationally expensive process where reads are mapped onto a reference sequence. Prior to such processing, however, there is a vast amount of information that can be ascertained from the reads, potentially obviating the need for processing, or allowing optimized mapping approaches to be deployed. Here, we present a method termed FlexTyper which facilitates a "reverse mapping" approach in which high throughput sequence queries, in the form of k-mer searches, are run against indexed short-read datasets in order to extract useful information. This reverse mapping approach enables the rapid counting of target sequences of interest. We demonstrate FlexTyper's utility for recovering depth of coverage, and accurate genotyping of SNP sites across the human genome. We show that genotyping unmapped reads can correctly inform a sample's population, sex, and relatedness in a family setting. Detection of pathogen sequences within RNA-seq data was sensitive and accurate, performing comparably to existing methods, but with increased flexibility. We present two examples of ways in which this flexibility allows the analysis of genome features not well-represented in a linear reference. First, we analyze contigs from African genome sequencing studies, showing how they distribute across families from three distinct populations. Second, we show how gene-marking k-mers for the killer immune receptor locus allow allele detection in a region that is challenging for standard read mapping pipelines. The future adoption of the reverse mapping approach represented by FlexTyper will be enabled by more efficient methods for FM-index generation and biology-informed collections of reference queries. In the long-term, selection of population-specific references or weighting of edges in pan-population reference genome graphs will be possible using the FlexTyper approach. FlexTyper is available at https://github.com/wassermanlab/OpenFlexTyper.

JASPAR 2020: update of the open-access database of transcription factor binding profiles.

JASPAR (http://jaspar.genereg.net) is an open-access database of curated, non-redundant transcription factor (TF)-binding profiles stored as position frequency matrices (PFMs) for TFs across multiple species in six taxonomic groups. In this 8th release of JASPAR, the CORE collection has been expanded with 245 new PFMs (169 for vertebrates, 42 for plants, 17 for nematodes, 10 for insects, and 7 for fungi), and 156 PFMs were updated (125 for vertebrates, 28 for plants and 3 for insects). These new profiles represent an 18% expansion compared to the previous release. JASPAR 2020 comes with a novel collection of unvalidated TF-binding profiles for which our curators did not find orthogonal supporting evidence in the literature. This collection has a dedicated web form to engage the community in the curation of unvalidated TF-binding profiles. Moreover, we created a Q&A forum to ease the communication between the user community and JASPAR curators. Finally, we updated the genomic tracks, inference tool, and TF-binding profile similarity clusters. All the data is available through the JASPAR website, its associated RESTful API, and through the JASPAR2020 R/Bioconductor package.

GeneBreaker: Variant simulation to improve the diagnosis of Mendelian rare genetic diseases.

Mendelian rare genetic diseases affect 5%-10% of the population, and with over 5300 genes responsible for ∼7000 different diseases, they are challenging to diagnose. The use of whole-genome sequencing (WGS) has bolstered the diagnosis rate significantly. The effective use of WGS relies on the ability to identify the disrupted gene responsible for disease phenotypes. This process involves genomic variant calling and prioritization, and is the beneficiary of improvements to sequencing technology, variant calling approaches, and increased capacity to prioritize genomic variants with potential pathogenicity. As analysis pipelines continue to improve, careful testing of their efficacy is paramount. However, real-life cases typically emerge anecdotally, and utilization of clinically sensitive and identifiable data for testing pipeline improvements is regulated and limiting. We identified the need for a gene-based variant simulation framework that can create mock rare disease scenarios, utilizing known pathogenic variants or through the creation of novel gene-disrupting variants. To fill this need, we present GeneBreaker, a tool that creates synthetic rare disease cases with utility for benchmarking variant calling approaches, testing the efficacy of variant prioritization, and as an educational mechanism for training diagnostic practitioners in the expanding field of genomic medicine. GeneBreaker is freely available at http://GeneBreaker.cmmt.ubc.ca.

Human MiniPromoters for ocular-rAAV expression in ON bipolar, cone, corneal, endothelial, Müller glial, and PAX6 cells.

Small and cell-type restricted promoters are important tools for basic and preclinical research, and clinical delivery of gene therapies. In clinical gene therapy, ophthalmic trials have been leading the field, with over 50% of ocular clinical trials using promoters that restrict expression based on cell type. Here, 19 human DNA MiniPromoters were bioinformatically designed for rAAV, tested by neonatal intravenous delivery in mouse, and successful MiniPromoters went on to be tested by intravitreal, subretinal, intrastromal, and/or intravenous delivery in adult mouse. We present promoter development as an overview for each cell type, but only show results in detail for the recommended MiniPromoters: Ple265 and Ple341 (PCP2) ON bipolar, Ple349 (PDE6H) cone, Ple253 (PITX3) corneal stroma, Ple32 (CLDN5) endothelial cells of the blood-retina barrier, Ple316 (NR2E1) Müller glia, and Ple331 (PAX6) PAX6 positive. Overall, we present a resource of new, redesigned, and improved MiniPromoters for ocular gene therapy that range in size from 784 to 2484 bp, and from weaker, equal, or stronger in strength relative to the ubiquitous control promoter smCBA. All MiniPromoters will be useful for therapies involving small regulatory RNA and DNA, and proteins ranging from 517 to 1084 amino acids, representing 62.9-90.2% of human proteins.

Secondary biogenic amine deficiencies: genetic etiology, therapeutic interventions, and clinical effects.

Monoamine neurotransmitter disorders present predominantly with neurologic features, including dystonic or dyskinetic cerebral palsy and movement disorders. Genetic conditions that lead to secondary defects in the synthesis, catabolism, transport, and metabolism of biogenic amines can lead to neurotransmitter abnormalities, which can present with similar features. Eleven patients with secondary neurotransmitter abnormalities were enrolled between 2011 and 2015. All patients underwent research-based whole exome and/or whole genome sequencing (WES/WGS). A trial of treatment with levodopa/carbidopa and 5-hydroxytryptophan was initiated. In six families with abnormal neurotransmitter profiles and neurological phenotypes, variants in known disease-causing genes (KCNJ6, SCN2A, CSTB in 2 siblings, NRNX1, KIF1A and PAK3) were identified, while one patient had a variant of uncertain significance in a candidate gene (DLG4) that may explain her phenotype. In 3 patients, no compelling candidate genes were identified. A trial of neurotransmitter replacement therapy led to improvement in motor and behavioral symptoms in all but two patients. The patient with KCNJ6 variant did not respond to L-dopa therapy, but rather experienced increased dyskinetic movements even at low dose of medication. The patient's symptoms harboring the NRNX1 deletion remained unaltered. This study demonstrates the utility of genome-wide sequencing in further understanding the etiology and pathophysiology of neurometabolic conditions, and the potential of secondary neurotransmitter deficiencies to serve as novel therapeutic targets. As there was a largely favorable response to therapy in our case series, a careful trial of neurotransmitter replacement therapy should be considered in patients with cerebrospinal fluid (CSF) monoamines below reference range.

Gene expression models based on transcription factor binding events confer insight into functional cis-regulatory variants.

MOTIVATION: Deciphering the functional roles of cis-regulatory variants is a critical challenge in genome analysis and interpretation. It has been hypothesized that altered transcription factor (TF) binding events are a central mechanism by which cis-regulatory variants impact gene expression levels. However, we lack a computational framework to understand and quantify such mechanistic contributions. RESULTS: We present TF2Exp, a gene-based framework to predict the impact of altered TF-binding events on gene expression levels. Using data from lymphoblastoid cell lines, TF2Exp models were applied successfully to predict the expression levels of 3196 genes. Alterations within DNase I hypersensitive, CTCF-bound and tissue-specific TF-bound regions were the greatest contributing features to the models. TF2Exp models performed as well as models based on common variants, both in cross-validation and external validation. Combining TF alteration and common variant features can further improve model performance. Unlike variant-based models, TF2Exp models have the unique advantage to evaluate the functional impact of variants in linkage disequilibrium and uncommon variants. We find that adding TF-binding events altered only by uncommon variants could increase the number of predictable genes (R2 > 0.05). Taken together, TF2Exp represents a key step towards interpreting the functional roles of cis-regulatory variants in the human genome. AVAILABILITY AND IMPLEMENTATION: The code and model training results are publicly available at https://github.com/wqshi/TF2Exp. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.