scholarly journals Queryfuse is a sensitive algorithm for detection of gene-specific fusions

2020 ◽  
Author(s):  
Yuxiang Tan

ABSTRACTRecurrent chromosomal translocations, known as fusions, play important roles in carcinogenesis. They can serve as valuable diagnostic and therapeutic targets. RNA-seq is an ideal platform for detecting transcribed fusions, and computational methods have been developed to identify fusion transcripts from RNA-seq data. However, some transciptome realignment procedures for these methods are unnecessary, making this task computationally expensive and time consuming. Therefore, we have developed QueryFuse, a novel hypothesis-based algorithm that identifies gene-specific fusion from pre-aligned RNA-seq data. It is designed to help biologists quickly find and/or computationally validate fusions of interest, together with visualization and detailed properties of supporting reads. By aligning reads to Query genes at the pre-processing step with a more sensitive, memory intensive local aligner, QueryFuse can reduce alignment time and improve detection sensitivity.QueryFuse performed better or at comparable levels with two popular tools (deFuse and TopHatFusion) on both simulated and well-annotated cell-line datasets. Finally, using QueryFuse, we identified a novel fusion event with a potential therapeutic implication in clinical samples. Taken together, our results showed that QueryFuse is efficient and reliable for detecting gene-specific fusion events.

2021 ◽  
Vol 22 (3) ◽  
pp. 1146
Author(s):  
Reinhard Ullmann ◽  
Benjamin Valentin Becker ◽  
Simone Rothmiller ◽  
Annette Schmidt ◽  
Horst Thiermann ◽  
...  

Sulfur mustard (SM) is a chemical warfare agent that can damage DNA via alkylation and oxidative stress. Because of its genotoxicity, SM is cancerogenic and the progenitor of many chemotherapeutics. Previously, we developed an SM-resistant cell line via chronic exposure of the popular keratinocyte cell line HaCaT to increasing doses of SM over a period of 40 months. In this study, we compared the genomic landscape of the SM-resistant cell line HaCaT/SM to its sensitive parental line HaCaT in order to gain insights into genetic changes associated with continuous alkylation and oxidative stress. We established chromosome numbers by cytogenetics, analyzed DNA copy number changes by means of array Comparative Genomic Hybridization (array CGH), employed the genome-wide chromosome conformation capture technique Hi-C to detect chromosomal translocations, and derived mutational signatures by whole-genome sequencing. We observed that chronic SM exposure eliminated the initially prevailing hypotetraploid cell population in favor of a hyperdiploid one, which contrasts with previous observations that link polyploidization to increased tolerance and adaptability toward genotoxic stress. Furthermore, we observed an accumulation of chromosomal translocations, frequently flanked by DNA copy number changes, which indicates a high rate of DNA double-strand breaks and their misrepair. HaCaT/SM-specific single-nucleotide variants showed enrichment of C > A and T > A transversions and a lower rate of deaminated cytosines in the CpG dinucleotide context. Given the frequent use of HaCaT in toxicology, this study provides a valuable data source with respect to the original genotype of HaCaT and the mutational signatures associated with chronic alkylation and oxidative stress.


2021 ◽  
Vol 22 (3) ◽  
pp. 1407
Author(s):  
Hongxia Liu ◽  
Wang Zheng ◽  
Qianping Chen ◽  
Yuchuan Zhou ◽  
Yan Pan ◽  
...  

Nasopharyngeal carcinoma (NPC) is one of the most frequent head and neck malignant tumors and is majorly treated by radiotherapy. However, radiation resistance remains a serious obstacle to the successful treatment of NPC. The aim of this study was to discover the underlying mechanism of radioresistance and to elucidate novel genes that may play important roles in the regulation of NPC radiosensitivity. By using RNA-seq analysis of NPC cell line CNE2 and its radioresistant cell line CNE2R, lncRNA CASC19 was screened out as a candidate radioresistance marker. Both in vitro and in vivo data demonstrated that a high expression level of CASC19 was positively correlated with the radioresistance of NPC, and the radiosensitivity of NPC cells was considerably enhanced by knockdown of CASC19. The incidence of autophagy was enhanced in CNE2R in comparison with CNE2 and another NPC cell line HONE1, and silencing autophagy with LC3 siRNA (siLC3) sensitized NPC cells to irradiation. Furthermore, CASC19 siRNA (siCASC19) suppressed cellular autophagy by inhibiting the AMPK/mTOR pathway and promoted apoptosis through the PARP1 pathway. Our results revealed for the first time that lncRNA CASC19 contributed to the radioresistance of NPC by regulating autophagy. In significance, CASC19 might be a potential molecular biomarker and a new therapeutic target in NPC.


BMC Genomics ◽  
2021 ◽  
Vol 22 (1) ◽  
Author(s):  
Tracy M. Yamawaki ◽  
Daniel R. Lu ◽  
Daniel C. Ellwanger ◽  
Dev Bhatt ◽  
Paolo Manzanillo ◽  
...  

Abstract Background Elucidation of immune populations with single-cell RNA-seq has greatly benefited the field of immunology by deepening the characterization of immune heterogeneity and leading to the discovery of new subtypes. However, single-cell methods inherently suffer from limitations in the recovery of complete transcriptomes due to the prevalence of cellular and transcriptional dropout events. This issue is often compounded by limited sample availability and limited prior knowledge of heterogeneity, which can confound data interpretation. Results Here, we systematically benchmarked seven high-throughput single-cell RNA-seq methods. We prepared 21 libraries under identical conditions of a defined mixture of two human and two murine lymphocyte cell lines, simulating heterogeneity across immune-cell types and cell sizes. We evaluated methods by their cell recovery rate, library efficiency, sensitivity, and ability to recover expression signatures for each cell type. We observed higher mRNA detection sensitivity with the 10x Genomics 5′ v1 and 3′ v3 methods. We demonstrate that these methods have fewer dropout events, which facilitates the identification of differentially-expressed genes and improves the concordance of single-cell profiles to immune bulk RNA-seq signatures. Conclusion Overall, our characterization of immune cell mixtures provides useful metrics, which can guide selection of a high-throughput single-cell RNA-seq method for profiling more complex immune-cell heterogeneity usually found in vivo.


Oncotarget ◽  
2019 ◽  
Vol 10 (56) ◽  
pp. 5768-5779 ◽  
Author(s):  
Luana de O. Lopes ◽  
Jersey H. Maués ◽  
Hygor Ferreira-Fernandes ◽  
France K. Yoshioka ◽  
Severino C. Sousa Júnior ◽  
...  

2020 ◽  
Author(s):  
Viacheslav Mylka ◽  
Jeroen Aerts ◽  
Irina Matetovici ◽  
Suresh Poovathingal ◽  
Niels Vandamme ◽  
...  

ABSTRACTMultiplexing of samples in single-cell RNA-seq studies allows significant reduction of experimental costs, straightforward identification of doublets, increased cell throughput, and reduction of sample-specific batch effects. Recently published multiplexing techniques using oligo-conjugated antibodies or - lipids allow barcoding sample-specific cells, a process called ‘hashing’. Here, we compare the hashing performance of TotalSeq-A and -C antibodies, custom synthesized lipids and MULTI-seq lipid hashes in four cell lines, both for single-cell RNA-seq and single-nucleus RNA-seq. Hashing efficiency was evaluated using the intrinsic genetic variation of the cell lines. Benchmarking of different hashing strategies and computational pipelines indicates that correct demultiplexing can be achieved with both lipid- and antibody-hashed human cells and nuclei, with MULTISeqDemux as the preferred demultiplexing function and antibody-based hashing as the most efficient protocol on cells. Antibody hashing was further evaluated on clinical samples using PBMCs from healthy and SARS-CoV-2 infected patients, where we demonstrate a more affordable approach for large single-cell sequencing clinical studies, while simultaneously reducing batch effects.


2021 ◽  
Vol 12 ◽  
Author(s):  
Eric de Sousa ◽  
Joana R. Lérias ◽  
Antonio Beltran ◽  
Georgia Paraschoudi ◽  
Carolina Condeço ◽  
...  

Successful outcome of immune checkpoint blockade in patients with solid cancers is in part associated with a high tumor mutational burden (TMB) and the recognition of private neoantigens by T-cells. The quality and quantity of target recognition is determined by the repertoire of ‘neoepitope’-specific T-cell receptors (TCRs) in tumor-infiltrating lymphocytes (TIL), or peripheral T-cells. Interferon gamma (IFN-γ), produced by T-cells and other immune cells, is essential for controlling proliferation of transformed cells, induction of apoptosis and enhancing human leukocyte antigen (HLA) expression, thereby increasing immunogenicity of cancer cells. TCR αβ-dependent therapies should account for tumor heterogeneity and availability of the TCR repertoire capable of reacting to neoepitopes and functional HLA pathways. Immunogenic epitopes in the tumor-stroma may also be targeted to achieve tumor-containment by changing the immune-contexture in the tumor microenvironment (TME). Non protein-coding regions of the tumor-cell genome may also contain many aberrantly expressed, non-mutated tumor-associated antigens (TAAs) capable of eliciting productive anti-tumor immune responses. Whole-exome sequencing (WES) and/or RNA sequencing (RNA-Seq) of cancer tissue, combined with several layers of bioinformatic analysis is commonly used to predict possible neoepitopes present in clinical samples. At the ImmunoSurgery Unit of the Champalimaud Centre for the Unknown (CCU), a pipeline combining several tools is used for predicting private mutations from WES and RNA-Seq data followed by the construction of synthetic peptides tailored for immunological response assessment reflecting the patient’s tumor mutations, guided by MHC typing. Subsequent immunoassays allow the detection of differential IFN-γ production patterns associated with (intra-tumoral) spatiotemporal differences in TIL or peripheral T-cells versus TIL. These bioinformatics tools, in addition to histopathological assessment, immunological readouts from functional bioassays and deep T-cell ‘adaptome’ analyses, are expected to advance discovery and development of next-generation personalized precision medicine strategies to improve clinical outcomes in cancer in the context of i) anti-tumor vaccination strategies, ii) gauging mutation-reactive T-cell responses in biological therapies and iii) expansion of tumor-reactive T-cells for the cellular treatment of patients with cancer.


2021 ◽  
Author(s):  
Wu Biao ◽  
Yufeng Chen ◽  
Junlong Zhong ◽  
Shuping Zhong ◽  
Bin Wang ◽  
...  

Abstract Background: Rheumatoid arthritis (RA) is a common autoimmune disease that can occur at any age. If treatment is delayed, RA can seriously affect the patients’ quality of life. However, there is no diagnostic criteria for RA and the positive predictive value of the current biomarkers is moderate. Objective: to identify RA-associated susceptibility genes and explore their potential as a novel biomarker for diagnosis and evaluation of the prognosis of RA.Methods: Peripheral blood mononuclear cells (PBMCs) were collected from healthy human donors and RA patients. RNA-seq analyses were performed to identify the differentially expressed genes (DEGs) between RA and control samples. The PBMCs-mRNA in DEGs were further subjected to enrichment analysis. Furthermore, the hub genes and key modules associated with RA were screened by bioinformatics analyses. Then, the expression of hub genes in RA were assessed in mRNA expression profiles. Next, real time-quantitative PCR (RT-qPCR) analyses were performed to further confirm the expression of the hub genes from the PBMCs that collected from 47 patients with RA and 40 healthy controls. Finally, we evaluated the clinical characters for the candidate mRNAs.Results: RNA-seq analyses revealed the expression of 178 mRNAs from PBMCs were disregulated between the healthy controls and the RA patients. Bioinformatics analyses revealed 10 hub mRNAs. The top 3 significant functional modules screened from PPI network functionally were involved in DNA replication origin binding, chemokine activity, etc. After validating the 10 hub mRNAs in GSE93272 dataset and clinical samples, we identified 3 candidate mRNAs, including ASPM, DTL and RRM2. Among which, RRM2 showed great capacity in discriminating between remissive RA and active RA. Significant correlations were observed between DTL and IL-8, TNF-α, between RRM2 and CDAI, DAS-28, tender joints and swollen joints, respectively. The AUC values of ASPM, DTL and RRM2 were 0.654, 0.995 and 0.990, respectively.Conclusion: We successfully identified multiple candidate mRNAs associated with RA. RRM2 showed high diagnosis efficiency with the AUC of 0.990 (sensitivity=100%, specificity=97.5%). And RRM2 severed as an additional biomarker for evaluating disease activity. The findings provided a novel candidate biomarker for diagnosis and evaluation of the prognosis of RA.


2020 ◽  
Author(s):  
Romain Daveu ◽  
Caroline Hervet ◽  
Louane Sigrist ◽  
Davide Sassera ◽  
Aaron Jex ◽  
...  

AbstractWe studied a family of iflaviruses, a group of RNA viruses frequently found in arthropods, focusing on viruses associated with ticks. Our aim was to bring insight on the evolutionary dynamics of this group of viruses, which may interact with the biology of ticks. We explored systematically de novo RNA-Seq assemblies available for species of ticks which allowed to identify nine new genomes of iflaviruses. The phylogeny of virus sequences was not congruent with that of the tick hosts, suggesting recurrent host changes across tick genera along evolution. We identified five different variants with a complete or near-complete genome in Ixodes ricinus. These sequences were closely related, which allowed a fine-scale estimation of patterns of substitutions: we detected a strong excess of synonymous mutations suggesting evolution under strong positive selection. ISIV, a sequence found in the ISE6 cell line of Ixodes scapularis, was unexpectedly nearidentical with I. ricinus variants, suggesting a contamination of this cell line by I. ricinus material. Overall, our work constitutes a step in the understanding of the interactions between this family of viruses and ticks.


Blood ◽  
2019 ◽  
Vol 134 (Supplement_1) ◽  
pp. 3758-3758
Author(s):  
Jianping Li ◽  
Catalina Troche ◽  
Julia Hlavka Zhang ◽  
Jonathan Shrimp ◽  
Jacob S. Roth ◽  
...  

Despite improvements in chemotherapy that have increased the 5-year survival rates of pediatric ALL to close to 90%, 15-20% of patients may relapse with a very poor prognosis. Pediatric ALL patients, particularly those in relapse can harbor a specific point mutation (E1099K) in NSD2 (nuclear receptor binding SET domain protein 2) gene, also known as MMSET or WHSC1, which encodes a histone methyl transferase specific for H3K36me2. To understand the biology of mutant NSD2, we used CRISPR-Cas9 gene editing to disrupt the NSD2E1099K mutant allele in B-ALL cell lines (RCH-ACV and SEM) and T-ALL cell line (RPMI-8402) or insert the E1099K mutation into the NSD2WT T-ALL cell line (CEM) and B-ALL cell line (697). Cell lines in which the NSD2E1099K mutant allele is present display increased global levels of H3K36me2 and decreased H3K27me3. NSD2E1099Kcells demonstrate enhanced cell growth, colony formation and migration. NSD2E1099K mutant cell lines assayed by RNA-Seq exhibit an aberrant gene signature, mostly representing gene activation, with activation of signaling pathways, genes implicated in the epithelial mesenchymal transition and prominent expression of neural genes not generally found in hematopoietic tissues. Accordingly, NSD2E1099K cell lines showed prominent tropism to the central neural system in xenografts. To understand why this NSD2 mutations are identified prominently in children who relapse early from therapy for ALL, we performed high-throughput screening in our isogenic cell lines with the National Center for Advancing Translation Science (NCATS) Pharmaceutical Collection and other annotated chemical libraries and found that NSD2E1099K cells are resistant to glucocorticoids (GC) but not to other chemotherapeutic agents used to treat ALL such as vincristine, doxorubicin, cyclophosphamide, methotrexate, and 6-mercaptopurine. Accordingly, patient-derived-xenograft ALL cells with NSD2E1099K mutation were resistant to GC treatment. Reversion of NSD2E1099K mutation to NSD2WT restored GC sensitivity to both B- and T-ALL cell lines, which was accompanied by cell cycle arrest in G1 and induced-apoptosis. Furthermore, knock-in of the NSD2E1099K mutation conferred GC resistance to ALL cell lines by triggering cell cycle progression, proliferation and anti-apoptotic processes. Mice with NSD2E1099K xenografts were completely resistant to GC treatment while treatment of mice injected with isogenic NSD2WT cells led to significant tumor reduction and survival benefit. To illustrate these biological phenotypes and understand the molecular mechanism of GC resistance driven by NSD2E1099Kmutation, we investigated the GC-induced transcriptome, GC receptor (GR) binding sites and related epigenetic changes in isogenic ALL cell lines in response to GC treatment. RNA-Seq showed that GC transcriptional response was almost completely blocked in NSD2E1099K cells, especially in T-ALL cell lines, correlating with their lack of biological response. GC treatment activated apoptotic pathways and downregulated cell cycle and DNA repair pathways only in NSD2WT cells. The critical pro-apoptotic regulators BIM and BMF failed to be activated by GC in NSD2E1099K cells but were prominently activated when the NSD2 mutation was removed. Chromatin immunoprecipitation sequencing (ChIP-Seq) showed that, the NSD2E1099K mutation blocked the ability of GR and CTCF to bind most GC response elements (GREs) such as those within BIM and BMF. While GR binding in NSD2WT cells was accompanied by increased H3K27 acetylation and gene expression, this failed to occur in NSD2 mutant cells. Furthermore, we found that GR RNA and protein levels were repressed in ALL cells expressing NSD2E1099K and GC failed to induce GR expression in these cells. Paradoxically, while H3K27me3 levels were generally decreased in NSD2E1099K cells, we saw increased levels of H3K27me3 at the GRE within the GR gene body where GR itself and CTCF normally bind, suggesting a novel role for the polycomb repressive complex 2 and EZH2 inhibitors for this form of GC resistance. In conclusion, these studies demonstrate that NSD2E1099K mutation may play an important role in treatment failure of pediatric ALL relapse by interfering with the GR expression and its ability to bind and activate key target genes. Gene editing screens are being performed to understand how to overcome this resistance. Disclosures No relevant conflicts of interest to declare.


2021 ◽  
Author(s):  
Jasmijn A. Baaijens ◽  
Alessandro Zulli ◽  
Isabel M. Ott ◽  
Mary E. Petrone ◽  
Tara Alpert ◽  
...  

Effectively monitoring the spread of SARS-CoV-2 variants is essential to efforts to counter the ongoing pandemic. Wastewater monitoring of SARS-CoV-2 RNA has proven an effective and efficient technique to approximate COVID-19 case rates in the population. Predicting variant abundances from wastewater, however, is technically challenging. Here we show that by sequencing SARS-CoV-2 RNA in wastewater and applying computational techniques initially used for RNA-Seq quantification, we can estimate the abundance of variants in wastewater samples. We show by sequencing samples from wastewater and clinical isolates in Connecticut U.S.A. between January and April 2021 that the temporal dynamics of variant strains broadly correspond. We further show that this technique can be used with other wastewater sequencing techniques by expanding to samples taken across the United States in a similar timeframe. We find high variability in signal among individual samples, and limited ability to detect the presence of variants with clinical frequencies <10%; nevertheless, the overall trends match what we observed from sequencing clinical samples. Thus, while clinical sequencing remains a more sensitive technique for population surveillance, wastewater sequencing can be used to monitor trends in variant prevalence in situations where clinical sequencing is unavailable or impractical.


Sign in / Sign up

Export Citation Format

Share Document