Search

Summary: Properly and effectively managing reference datasets is an important task for many bioinformatics analyses. Refgenie is a reference asset management system that allows users to easily organize, retrieve and share such datasets. Here, we describe the integration of refgenie into the Galaxy platform. Server administrators are able to configure Galaxy to make use of reference datasets made available on a refgenie instance. In addition, a Galaxy Data Manager tool has been developed to provide a graphical interface to refgenie's remote reference retrieval functionality. A large collection of reference datasets has also been made available using the CVMFS (CernVM File System) repository from GalaxyProject.org, with mirrors across the USA, Canada, Europe and Australia, enabling easy use outside of Galaxy., Availability and Implementation: The ability of Galaxy to use refgenie assets was added to the core Galaxy framework in version 22.01, which is available from https://github.com/galaxyproject/galaxy under the Academic Free License version 3.0. The refgenie Data Manager tool can be installed via the Galaxy ToolShed, with source code managed at https://github.com/BlankenbergLab/galaxy-tools-blankenberg/tree/main/data_managers/data_manager_refgenie_pull and released using an MIT license. Access to existing data is also available through CVMFS, with instructions at https://galaxyproject.org/admin/reference-data-repo/. No new data were generated or analyzed in support of this research., (© The Author(s) 2022. Published by Oxford University Press.)

Background: The Public Health Alliance for Genomic Epidemiology (PHA4GE) (https://pha4ge.org) is a global coalition that is actively working to establish consensus standards, document and share best practices, improve the availability of critical bioinformatics tools and resources, and advocate for greater openness, interoperability, accessibility, and reproducibility in public health microbial bioinformatics. In the face of the current pandemic, PHA4GE has identified a need for a fit-for-purpose, open-source SARS-CoV-2 contextual data standard., Results: As such, we have developed a SARS-CoV-2 contextual data specification package based on harmonizable, publicly available community standards. The specification can be implemented via a collection template, as well as an array of protocols and tools to support both the harmonization and submission of sequence data and contextual information to public biorepositories., Conclusions: Well-structured, rich contextual data add value, promote reuse, and enable aggregation and integration of disparate datasets. Adoption of the proposed standard and practices will better enable interoperability between datasets and systems, improve the consistency and utility of generated data, and ultimately facilitate novel insights and discoveries in SARS-CoV-2 and COVID-19. The package is now supported by the NCBI's BioSample database., (© The Author(s) 2022. Published by Oxford University Press GigaScience.)

Influenza viruses continually evolve to evade population immunity, and the different lineages are assigned into clades based on shared mutations. We have developed a publicly available computational workflow, the Influenza Classification Suite, for rapid clade mapping of sequenced influenza viruses. This suite provides a user-friendly workflow implemented in Galaxy to automate clade calling and antigenic site extraction. Workflow input includes clade definition and amino acid index array files, which can be customized to identify any clades of interest. The Influenza Classification Suite provides rapid, high-resolution understanding of circulating influenza strain evolution to inform influenza vaccine effectiveness and the need for potential vaccine reformulation., (© 2020 The Authors. Influenza and Other Respiratory Viruses Published by John Wiley & Sons Ltd.)

MLST (multi-locus sequence typing) is a classic technique for genotyping bacteria, widely applied for pathogen outbreak surveillance. Traditionally, MLST is based on identifying sequence types from a small number of housekeeping genes. With the increasing availability of whole-genome sequencing data, MLST methods have evolved towards larger typing schemes, based on a few hundred genes [core genome MLST (cgMLST)] to a few thousand genes [whole genome MLST (wgMLST)]. Such large-scale MLST schemes have been shown to provide a finer resolution and are increasingly used in various contexts such as hospital outbreaks or foodborne pathogen outbreaks. This methodological shift raises new computational challenges, especially given the large size of the schemes involved. Very few available MLST callers are currently capable of dealing with large MLST schemes. We introduce MentaLiST, a new MLST caller, based on a k-mer voting algorithm and written in the Julia language, specifically designed and implemented to handle large typing schemes. We test it on real and simulated data to show that MentaLiST is faster than any other available MLST caller while providing the same or better accuracy, and is capable of dealing with MLST schemes with up to thousands of genes while requiring limited computational resources. MentaLiST source code and easy installation instructions using a Conda package are available at https://github.com/WGS-TB/MentaLiST.

Ultrasensitive methods for rare allele detection are critical to leverage the full potential offered by liquid biopsies. Here, we describe a novel molecular barcoding method for the precise detection and quantification of circulating tumor DNA (ctDNA). The major benefits of our design include straightforward and cost-effective production of barcoded adapters to tag individual DNA molecules before PCR and sequencing, and better control over cross-contamination between experiments. We validated our approach in a cohort of 24 patients with a broad spectrum of cancer diagnoses by targeting and quantifying single-nucleotide variants (SNVs), indels and genomic rearrangements in plasma samples. By using personalized panels targeting a priori known mutations, we demonstrate comprehensive error-suppression capabilities for SNVs and detection thresholds for ctDNA below 0.1%. We also show that our semi-degenerate barcoded adapters hold promise for noninvasive genotyping in the absence of tumor biopsies and monitoring of minimal residual disease in longitudinal plasma samples. The benefits demonstrated here include broad applicability, flexibility, affordability and reproducibility in the research and clinical settings.

Background: Inherited arrhythmia syndromes are responsible for a significant portion of autopsy-negative sudden unexpected death (SUD) cases, but molecular autopsy used to identify potentially causal variants is not routinely included in SUD investigations. We collaborated with a medical examiner's office to assist in finding a diagnosis for their autopsy-negative child SUD cases., Methods and Results: 191 child SUD cases (<5 years of age) were selected for analyses. Our next generation sequencing panel incorporated 38 inherited arrhythmia syndrome candidate genes and another 33 genes not previously investigated for variants that may underlie SUDY pathophysiology. Overall, we identified 11 potentially causal disease-associated variants in 12 cases, for an overall yield of 6.3%. We also identified 31 variants of uncertain significance in 36 cases and 16 novel variants predicted to be pathogenic in silico in 15 cases. The disease-associated variants were reported to the medical examiner to notify surviving relatives and recommend clinical assessment., Conclusions: We have identified variants that may assist in the diagnosis of at least 6.3% of autopsy-negative child SUD cases and reduce risk of future SUD in surviving relatives. We recommend a cautious approach to variant interpretation. We also suggest inclusion of cardiomyopathy genes as well as other candidate SUD genes in molecular autopsy analyses., (© 2017 American Heart Association, Inc.)

Background: A plethora of options to detect mutations in tumor-derived DNA currently exist but each suffers limitations in analytical sensitivity, cost, or scalability. Droplet digital PCR (ddPCR) is an appealing technology for detecting the presence of specific mutations based on a priori knowledge and can be applied to tumor biopsies, including formalin-fixed paraffin embedded (FFPE) tissues. More recently, ddPCR has gained popularity in its utility in quantifying circulating tumor DNA., Methods: We have developed a suite of novel ddPCR assays for detecting recurrent mutations that are prevalent in common B-cell non-Hodgkin lymphomas (NHLs), including diffuse large B-cell lymphoma, follicular lymphoma, and lymphoplasmacytic lymphoma. These assays allowed the differentiation and counting of mutant and wild-type molecules using one single hydrolysis probe. We also implemented multiplexing that allowed the simultaneous detection of distinct mutations and an "inverted" ddPCR assay design, based on employing probes matching wild-type alleles, capable of detecting the presence of multiple single nucleotide polymorphisms., Results: The assays successfully detected and quantified somatic mutations commonly affecting enhancer of zeste 2 polycomb repressive complex 2 subunit (EZH2) (Y641) and signal transducer and activator of transcription 6 (STAT6) (D419) hotspots in fresh tumor, FFPE, and liquid biopsies. The "inverted" ddPCR approach effectively reported any single nucleotide variant affecting either of these 2 hotspots as well. Finally, we could effectively multiplex hydrolysis probes targeting 2 additional lymphoma-related hotspots: myeloid differentiation primary response 88 (MYD88; L265P) and cyclin D3 (CCND3; I290R)., Conclusions: Our suite of ddPCR assays provides sufficient analytical sensitivity and specificity for either the invasive or noninvasive detection of multiple recurrent somatic mutations in B-cell NHLs., (© 2016 American Association for Clinical Chemistry.)

The majority of diffuse large B-cell lymphoma (DLBCL) tumors contain mutations in histone-modifying enzymes (HMEs), indicating a potential therapeutic benefit of histone deacetylase inhibitors (HDIs), and preclinical data suggest that HDIs augment the effect of rituximab. In this randomized phase 2 study, we evaluated the response rate and toxicity of panobinostat, a pan-HDI administered 30 mg orally 3 times weekly, with or without rituximab, in 40 patients with relapsed or refractory de novo (n = 27) or transformed (n = 13) DLBCL. Candidate genes and whole exomes were sequenced in relapse tumor biopsies to search for molecular correlates, and these data were used to quantify circulating tumor DNA (ctDNA) in serial plasma samples. Eleven of 40 patients (28%) responded to panobinostat (95% confidence interval [CI] 14.6-43.9) and rituximab did not increase responses. The median duration of response was 14.5 months (95% CI 9.4 to "not reached"). At time of data censoring, 6 of 11 patients had not progressed. Of the genes tested for mutations, only those in MEF2B were significantly associated with response. We detected ctDNA in at least 1 plasma sample from 96% of tested patients. A significant increase in ctDNA at day 15 relative to baseline was strongly associated with lack of response (sensitivity 71.4%, specificity 100%). We conclude that panobinostat induces very durable responses in some patients with relapsed DLBCL, and early responses can be predicted by mutations in MEF2B or a significant change in ctDNA level at 15 days after treatment initiation. This clinical trial was registered at www.ClinicalTrials.gov (#NCT01238692)., (© 2016 by The American Society of Hematology.)

Purpose: Relapsed or refractory diffuse large B-cell lymphoma (rrDLBCL) is fatal in 90% of patients, and yet little is known about its biology., Experimental Design: Using exome sequencing, we characterized the mutation profiles of 38 rrDLBCL biopsies obtained at the time of progression after immunochemotherapy. To identify genes that may be associated with relapse, we compared the mutation frequency in samples obtained at relapse to an unrelated cohort of 138 diagnostic DLBCLs and separately amplified specific mutations in their matched diagnostic samples to identify clonal expansions., Results: On the basis of a higher frequency at relapse and evidence for clonal selection, TP53, FOXO1, MLL3 (KMT2C), CCND3, NFKBIZ, and STAT6 emerged as top candidate genes implicated in therapeutic resistance. We observed individual examples of clonal expansions affecting genes whose mutations had not been previously associated with DLBCL including two regulators of NF-κB: NFKBIE and NFKBIZ We detected mutations that may be affect sensitivity to novel therapeutics, such as MYD88 and CD79B mutations, in 31% and 23% of patients with activated B-cell-type of rrDLBCL, respectively. We also identified recurrent STAT6 mutations affecting D419 in 36% of patients with the germinal center B (GCB) cell rrDLBCL. These were associated with activated JAK/STAT signaling, increased phospho-STAT6 protein expression and increased expression of STAT6 target genes., Conclusions: This work improves our understanding of therapeutic resistance in rrDLBCL and has identified novel therapeutic opportunities especially for the high-risk patients with GCB-type rrDLBCL. Clin Cancer Res; 22(9); 2290-300. ©2015 AACR., (©2015 American Association for Cancer Research.)

Non-Hodgkin lymphomas (NHLs) are the most common cancer to affect pet dogs. In contrast to the many genes whose mutation contributes to lymphomagenesis in humans, relatively little is known about the acquired genetic alterations that lead to canine B-cell lymphomas (cBCLs). We performed a survey of 84 canine NHL tumors to identify genes affected by somatic point mutations. We found mutations affecting TRAF3, which encodes a negative regulator of nuclear factor (NF)-κB, to be a common feature of cBCLs, with mutations observed in 44% of tumors including a combination of somatic and rare germ-line variants. Overall, 30% of the tumors contained ≥1 somatic TRAF3 mutation. The majority of mutations are predicted to cause loss of TRAF3 protein including those impacting reading frame and splicing. To determine whether TRAF3 loss might be relevant to human NHL, we also analyzed 148 human diffuse large B-cell lymphoma (DLBCL) tumors and identified loss of TRAF3 as a common event, affecting ∼9% of DLBCLs, and reduced expression of TRAF3 among deleted cases. This study implicates mutations affecting NF-κB activity as a novel genetic commonality between human and canine NHLs and supports the potential utility of cBCLs with mutated TRAF3 as a model of the more aggressive activated B-cell subgroup of DLBCL., (© 2015 by The American Society of Hematology.)

Background: Human papilloma virus (HPV) prevalence studies performed in different regions and population groups across Canada would inform public health decisions regarding implementation of anti-HPV vaccines., Methods: A total of 8,700 liquid-based specimens from 8,660 women aged 13-86 from throughout British Columbia were collected. DNA was isolated from 4,980 of these samples and assessed for HPV prevalence and type distribution. HPV was detected by PCR analysis using tagged GP5+/6+ consensus primers to amplify the L1 region of HPV; typing was done by bi-directional sequencing of PCR products., Results: Overall HPV prevalence was 16.8% (age adjusted 15.5%). Prevalence of high-risk HPV was 13.9, and 10.7% of samples contained HPV16. HPV prevalence was highest in the youngest group of women (<20 years). One-third of HPV positive samples contained more than one HPV type. Percentages of low-grade (LGIL) and high-grade intraepithelial lesions (HGIL) containing high-risk HPV are 52.3 and 79.4%, respectively., Conclusions: Overall HPV prevalence in this study is within the range of estimates from other studies. The prevalence of HPV16 is higher than what is found in other Canadian and international studies. HPV16 and HPV18 compose a majority of the high-risk virus in this study. Use of current HPV vaccines could considerably reduce HPV-related conditions including cervical cancer and procedures such as colposcopy.

Telomeres protect against DNA degradation at the ends of linear chromosomes. The number of telomere repeats is reduced over time in human aging. Using flow FISH we have assessed telomere length in 134 exceptionally healthy seniors aged 85 or older who have never been diagnosed with cancer, cardiovascular disease, major pulmonary disease, diabetes or Alzheimer disease (the 'Super-seniors') and 47 randomly-ascertained mid-life individuals aged 40-50 years. We compared their telomere lengths to a reference interval based on 400 individuals aged 1-100 years and show that Super-seniors do not have exceptionally long telomeres for their age. Consistent with the known trend of telomere shortening over time; however, they have shorter telomeres than the younger control group. Furthermore, we show that variability in telomere length was lower in the Super-seniors than in the mid-life controls or the reference data. Reduced telomere length variation was observed for lymphocytes, CD45RA-positive T-cells and memory T-cells. These results suggest that individuals, some types of their somatic cells, or both, may be selected for an optimal rather than extreme telomere length. Selection of individuals and/or cells that have an optimal telomere repeat length could contribute to disease resistance and promote healthy aging.