RESUMEN
Carbapenem-resistant Acinetobacter baumannii (CRAB) has emerged as a major global pathogen with limited treatment options1. No new antibiotic chemical class with activity against A. baumannii has reached patients in over 50 years1. Here we report the identification and optimization of tethered macrocyclic peptide (MCP) antibiotics with potent antibacterial activity against CRAB. The mechanism of action of this molecule class involves blocking the transport of bacterial lipopolysaccharide from the inner membrane to its destination on the outer membrane, through inhibition of the LptB2FGC complex. A clinical candidate derived from the MCP class, zosurabalpin (RG6006), effectively treats highly drug-resistant contemporary isolates of CRAB both in vitro and in mouse models of infection, overcoming existing antibiotic resistance mechanisms. This chemical class represents a promising treatment paradigm for patients with invasive infections due to CRAB, for whom current treatment options are inadequate, and additionally identifies LptB2FGC as a tractable target for antimicrobial drug development.
Asunto(s)
Antibacterianos , Lipopolisacáridos , Proteínas de Transporte de Membrana , Animales , Humanos , Ratones , Acinetobacter baumannii/efectos de los fármacos , Acinetobacter baumannii/metabolismo , Antibacterianos/clasificación , Antibacterianos/farmacología , Farmacorresistencia Bacteriana Múltiple/efectos de los fármacos , Lipopolisacáridos/metabolismo , Pruebas de Sensibilidad Microbiana , Proteínas de Transporte de Membrana/metabolismo , Transporte Biológico/efectos de los fármacos , Modelos Animales de Enfermedad , Infecciones por Acinetobacter/tratamiento farmacológico , Infecciones por Acinetobacter/microbiología , Desarrollo de MedicamentosRESUMEN
MOTIVATION: Detection of random errors and systematic biases is a crucial step of a robust pipeline for processing high-throughput sequencing (HTS) data. Bioinformatics software tools capable of performing this task are available, either for general analysis of HTS data or targeted to a specific sequencing technology. However, most of the existing QC instruments only allow processing of one sample at a time. RESULTS: Qualimap 2 represents a next step in the QC analysis of HTS data. Along with comprehensive single-sample analysis of alignment data, it includes new modes that allow simultaneous processing and comparison of multiple samples. As with the first version, the new features are available via both graphical and command line interface. Additionally, it includes a large number of improvements proposed by the user community. AVAILABILITY AND IMPLEMENTATION: The implementation of the software along with documentation is freely available at http://www.qualimap.org. CONTACT: meyer@mpiib-berlin.mpg.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Asunto(s)
Genómica/métodos , Secuenciación de Nucleótidos de Alto Rendimiento/métodos , Control de Calidad , Alineación de Secuencia/normas , Análisis de Secuencia de ADN/métodos , Programas Informáticos , Algoritmos , HumanosRESUMEN
Next-generation sequencing (NGS) technologies are revolutionizing genome research, and in particular, their application to transcriptomics (RNA-seq) is increasingly being used for gene expression profiling as a replacement for microarrays. However, the properties of RNA-seq data have not been yet fully established, and additional research is needed for understanding how these data respond to differential expression analysis. In this work, we set out to gain insights into the characteristics of RNA-seq data analysis by studying an important parameter of this technology: the sequencing depth. We have analyzed how sequencing depth affects the detection of transcripts and their identification as differentially expressed, looking at aspects such as transcript biotype, length, expression level, and fold-change. We have evaluated different algorithms available for the analysis of RNA-seq and proposed a novel approach--NOISeq--that differs from existing methods in that it is data-adaptive and nonparametric. Our results reveal that most existing methodologies suffer from a strong dependency on sequencing depth for their differential expression calls and that this results in a considerable number of false positives that increases as the number of reads grows. In contrast, our proposed method models the noise distribution from the actual data, can therefore better adapt to the size of the data set, and is more effective in controlling the rate of false discoveries. This work discusses the true potential of RNA-seq for studying regulation at low expression ranges, the noise within RNA-seq data, and the issue of replication.
Asunto(s)
Algoritmos , Etiquetas de Secuencia Expresada , Perfilación de la Expresión Génica/métodos , Regulación de la Expresión Génica/fisiología , Modelos Genéticos , Análisis de Secuencia por Matrices de Oligonucleótidos/métodos , HumanosRESUMEN
A growing body of evidence suggests the non-protein coding human genome is of vital importance for human cell function. Besides small RNAs, the diverse class of long non-coding RNAs (lncRNAs) recently came into focus. However, their relevance for infection, a major evolutionary driving force, remains elusive. Using two commercially available microarray systems, namely NCode™ and Sureprint™ G3, we identified differential expression of 42 ncRNAs during influenza A virus (IAV) infection in human lung epithelial cells. This included several classes of lncRNAs, including large intergenic ncRNAs (lincRNAs). As analyzed by qRT-PCR, expression of one lincRNA, which we termed virus inducible lincRNA (VIN), is induced by several IAV strains (H1N1, H3N2, H7N7) as well as vesicular stomatitis virus. However, we did not observe an induction of VIN by influenza B virus, treatment with RNA mimics, or IFNß. Thus, VIN expression seems to be a specific response to certain viral infections. RNA fractionation and RNA-FISH experiments revealed that VIN is localized to the host cell nucleus. Most importantly, we show that abolition of VIN by RNA interference restricts IAV replication and viral protein synthesis, highlighting the relevance of this lincRNA for productive IAV infection. Our observations suggest that viral pathogens interfere with the non-coding portion of the human genome, thereby guaranteeing their successful propagation, and that the expression of VIN correlates with their virulence. Consequently, our study provides a novel approach for understanding virus pathogenesis in greater detail, which will enable future design of new antiviral strategies targeting the host's non-protein coding genome.
Asunto(s)
Virus de la Influenza A/fisiología , Riñón/virología , Pulmón/virología , ARN Largo no Codificante/genética , Virus de la Estomatitis Vesicular Indiana/fisiología , Animales , Línea Celular , Núcleo Celular/genética , Perros , Células Epiteliales/metabolismo , Células Epiteliales/virología , Regulación de la Expresión Génica , Humanos , Virus de la Influenza A/clasificación , Células de Riñón Canino Madin Darby , Análisis de Secuencia por Matrices de Oligonucleótidos , ARN Interferente Pequeño/farmacología , Proteínas Virales/metabolismo , Replicación ViralRESUMEN
MOTIVATION: The sequence alignment/map (SAM) and the binary alignment/map (BAM) formats have become the standard method of representation of nucleotide sequence alignments for next-generation sequencing data. SAM/BAM files usually contain information from tens to hundreds of millions of reads. Often, the sequencing technology, protocol and/or the selected mapping algorithm introduce some unwanted biases in these data. The systematic detection of such biases is a non-trivial task that is crucial to drive appropriate downstream analyses. RESULTS: We have developed Qualimap, a Java application that supports user-friendly quality control of mapping data, by considering sequence features and their genomic properties. Qualimap takes sequence alignment data and provides graphical and statistical analyses for the evaluation of data. Such quality-control data are vital for highlighting problems in the sequencing and/or mapping processes, which must be addressed prior to further analyses. AVAILABILITY: Qualimap is freely available from http://www.qualimap.org.
Asunto(s)
Secuenciación de Nucleótidos de Alto Rendimiento , Alineación de Secuencia/normas , Programas Informáticos , Algoritmos , Genómica , Humanos , Control de CalidadRESUMEN
AIM: We compared the expression of genes related to inflammatory and cytotoxic functions between MSI and MSS (HLA-class I-negative and HLA-class I-positive) colorectal cancers (CRCs), seeking evidence of differences in inflammatory mediators and cytotoxic T-cell responses. Twenty-two CRCs were divided into three study groups as a function of HLA class I expression and MSI phenotype: 8 MSI tumours, 6 MSS/HLA- tumours and 6 MSS/HLA+ tumours (controls). FINDINGS: A first comparison between eight MSI and six MSS/HLA-positive (control) cancers, based on microarray analysis on an Affymetrix(®) HG-U133-Plus-PM plate, identified 1974 differentially expressed genes (P < 0.05). We grouped genes in Gene Ontology functional categories: apoptotic programme (72 genes, P = 5.5·10(-3)), leucocyte activation (43 genes, P = 1.8·10(-5)), T-cell activation (24 genes, P = 6.3·10(-4)), inflammatory response (40 genes, 2.3·10(-2)) and cytokine production (10 genes, P = 1.9·10(-2)). Real-time PCR and immunohistochemical evaluation were used to validate the data, finding that increased mRNA levels of pro-inflammatory cytokines and cytotoxic mediators were associated with greater infiltration by CD8+T lymphocytes in the MSI group (P < 0.001). Finally, HLA-class I-negative tumours were not grouped together but rather in accordance with features of the gene expression profile of MSI or MSS tumours. As expected, genes associated with antigen processing machinery and MHC class I molecules (TAP2, B2m) were downregulated in MSS/HLA-class I-negative CRCs (n = 6) in comparison to controls. CONCLUSIONS: In conclusion, microarray and immunohistochemical data may be useful to comprehensively assess tumour-host interactions and differentiate MSI from MSS cancers. The two types of tumour, MSI/HLA-class I-negative and MSS/HLA-class I-negative, showed marked differences in the composition and intensity of infiltrating leucocytes, suggesting that their immune escape strategies involve distinct pathways.
Asunto(s)
Neoplasias Colorrectales/genética , Neoplasias Colorrectales/inmunología , Estudio de Asociación del Genoma Completo , Inestabilidad Genómica , Antígenos de Histocompatibilidad Clase I/genética , Escape del Tumor/inmunología , Adulto , Anciano , Anciano de 80 o más Años , Neoplasias Colorrectales/patología , Femenino , Perfilación de la Expresión Génica , Antígenos de Histocompatibilidad Clase I/inmunología , Humanos , Inmunohistoquímica , Masculino , Persona de Mediana Edad , Reacción en Cadena en Tiempo Real de la PolimerasaRESUMEN
MOTIVATION: The development of the omics technologies such as transcriptomics, proteomics and metabolomics has made possible the realization of systems biology studies where biological systems are interrogated at different levels of biochemical activity (gene expression, protein activity and/or metabolite concentration). An effective approach to the analysis of these complex datasets is the joined visualization of the disparate biomolecular data on the framework of known biological pathways. RESULTS: We have developed the Paintomics web server as an easy-to-use bioinformatics resource that facilitates the integrated visual analysis of experiments where transcriptomics and metabolomics data have been measured on different conditions for the same samples. Basically, Paintomics takes complete transcriptomics and metabolomics datasets, together with lists of significant gene or metabolite changes, and paints this information on KEGG pathway maps. AVAILABILITY: Paintomics is freely available at http://www.paintomics.org.
Asunto(s)
Perfilación de la Expresión Génica/métodos , Metabolómica/métodos , Programas Informáticos , Arabidopsis/genética , Arabidopsis/metabolismo , Biología Computacional , Internet , Biología de SistemasRESUMEN
Transcriptomic studies have revealed a large number of uncharacterized genes that are differentially expressed in biofilms, which may be important in regulating biofilm phenotypes such as resistance to antimicrobial agents. To identify biofilm genes of unknown function in P. aeruginosa, we made use of RNA-seq and selected 27 uncharacterized genes that were induced upon biofilm growth. Biofilms by respective mutants were subsequently analyzed for two biofilm characteristics, the biofilm architecture and drug susceptibility. The screen revealed 12 out of 27 genes to contribute to biofilm formation and 13 drug susceptibility, with 8 genes affecting both biofilm phenotypes. Amongst the genes affecting both biofilm phenotypes was PA2146, encoding a small hypothetical protein that exhibited some of the most substantial increases in transcript abundance during biofilm growth by P. aeruginosa PAO1 and clinical isolates. PA2146 is highly conserved in É£-proteobacteria. Inactivation of PA2146 affected both biofilm phenotypes in P. aeruginosa PAO1, with inactivation of homologs in Klebsiella pneumoniae and Escherichia coli having similar effects. Heterologous expression of PA2146 homologs complemented the P. aeruginosa ∆PA2146, suggesting that PA2146 homologs substitute for and play a similar role as PA2146 in P. aeruginosa.
Asunto(s)
Gammaproteobacteria , Regulación Bacteriana de la Expresión Génica , Proteínas Bacterianas/genética , Proteínas Bacterianas/metabolismo , Biopelículas , Tolerancia a Medicamentos , Gammaproteobacteria/metabolismo , Pseudomonas aeruginosaRESUMEN
Understanding the dynamics and mechanisms of acquired drug resistance across major classes of antibiotics and bacterial pathogens is of critical importance for the optimization of current anti-infective therapies and the development of novel ones. To systematically address this challenge, we developed a workflow combining experimental evolution in a morbidostat continuous culturing device with deep genomic sequencing of population samples collected in time series. This approach was applied to the experimental evolution of six populations of Escherichia coli BW25113 towards acquiring resistance to triclosan (TCS), an antibacterial agent in various consumer products. This study revealed the rapid emergence and expansion (up to 100% in each culture within 4 days) of missense mutations in the fabI gene, encoding enoyl-acyl carrier protein reductase, the known TCS molecular target. A follow-up analysis of isolated clones showed that distinct amino acid substitutions increased the drug IC90 in a 3-16-fold range, reflecting their proximity to the TCS-binding site. In contrast to other antibiotics, efflux-upregulating mutations occurred only rarely and with low abundance. Mutations in several other genes were detected at an earlier stage of evolution. Most notably, three distinct amino acid substitutions were mapped in the C-terminal periplasmic domain of CadC protein, an acid stress-responsive transcriptional regulator. While these mutations do not confer robust TCS resistance, they appear to play a certain, yet unknown, role in adaptation to relatively low drug pressure. Overall, the observed evolutionary trajectories suggest that the FabI enzyme is the sole target of TCS (at least up to the ~50 µm level), and amino acid substitutions in the TCS-binding site represent the main mechanism of robust TCS resistance in E. coli. This model study illustrates the potential utility of the established morbidostat-based approach for uncovering resistance mechanisms and target identification for novel drug candidates with yet unknown mechanisms of action.
Asunto(s)
Antibacterianos/farmacología , Bacterias/efectos de los fármacos , Bacterias/genética , Farmacorresistencia Bacteriana/genética , Genómica , Triclosán/farmacología , Sustitución de Aminoácidos , Escherichia coli/genética , Genes Bacterianos/genética , Mutación , Secuenciación Completa del GenomaRESUMEN
Resistance to the broad-spectrum antibiotic ciprofloxacin is detected at high rates for a wide range of bacterial pathogens. To investigate the dynamics of ciprofloxacin resistance development, we applied a comparative resistomics workflow for three clinically relevant species of Gram-negative bacteria: Escherichia coli, Acinetobacter baumannii, and Pseudomonas aeruginosa. We combined experimental evolution in a morbidostat with deep sequencing of evolving bacterial populations in time series to reveal both shared and unique aspects of evolutionary trajectories. Representative clone characterization by sequencing and MIC measurements enabled direct assessment of the impact of mutations on the extent of acquired drug resistance. In all three species, we observed a two-stage evolution: (i) early ciprofloxacin resistance reaching 4- to 16-fold the MIC for the wild type, commonly as a result of single mutations in DNA gyrase target genes (gyrA or gyrB), and (ii) additional genetic alterations affecting the transcriptional control of the drug efflux machinery or secondary target genes (DNA topoisomerase parC or parE). IMPORTANCE The challenge of spreading antibiotic resistance calls for systematic efforts to develop more "irresistible" drugs based on a deeper understanding of dynamics and mechanisms of antibiotic resistance acquisition. To address this challenge, we have established a comparative resistomics approach which combines experimental evolution in a continuous-culturing device, the morbidostat, with ultradeep sequencing of evolving microbial populations to identify evolutionary trajectories (mutations and genome rearrangements) leading to antibiotic resistance over a range of target pathogens. Here, we report the comparative resistomics study of three Gram-negative bacteria (Escherichia coli, Acinetobacter baumannii, and Pseudomonas aeruginosa), which revealed shared and species-specific aspects of the evolutionary landscape leading to robust resistance against the clinically important antibiotic ciprofloxacin. Despite some differences between morbidostat-deduced mutation profiles and those observed in clinical isolates of individual species, a cross-species comparative resistomics approach allowed us to recapitulate all types of clinically relevant ciprofloxacin resistance mechanisms. This observation supports the anticipated utility of this approach in guiding rational optimization of treatment regimens for current antibiotics and the development of novel antibiotics with minimized resistance propensities.
Asunto(s)
Antibacterianos/farmacología , Ciprofloxacina/farmacología , Farmacorresistencia Bacteriana/genética , Bacterias Gramnegativas/efectos de los fármacos , Bacterias Gramnegativas/genética , Sustitución de Aminoácidos , Bacterias Gramnegativas/clasificación , Pruebas de Sensibilidad Microbiana , Mutación/efectos de los fármacosRESUMEN
BACKGROUND: Transcription factors (TFs) control transcription by binding to specific regions of DNA called transcription factor binding sites (TFBSs). The identification of TFBSs is a crucial problem in computational biology and includes the subtask of predicting the location of known TFBS motifs in a given DNA sequence. It has previously been shown that, when scoring matches to known TFBS motifs, interdependencies between positions within a motif should be taken into account. However, this remains a challenging task owing to the fact that sequences similar to those of known TFBSs can occur by chance with a relatively high frequency. Here we present a new method for matching sequences to TFBS motifs based on intuitionistic fuzzy sets (IFS) theory, an approach that has been shown to be particularly appropriate for tackling problems that embody a high degree of uncertainty. RESULTS: We propose SCintuit, a new scoring method for measuring sequence-motif affinity based on IFS theory. Unlike existing methods that consider dependencies between positions, SCintuit is designed to prevent overestimation of less conserved positions of TFBSs. For a given pair of bases, SCintuit is computed not only as a function of their combined probability of occurrence, but also taking into account the individual importance of each single base at its corresponding position. We used SCintuit to identify known TFBSs in DNA sequences. Our method provides excellent results when dealing with both synthetic and real data, outperforming the sensitivity and the specificity of two existing methods in all the experiments we performed. CONCLUSIONS: The results show that SCintuit improves the prediction quality for TFs of the existing approaches without compromising sensitivity. In addition, we show how SCintuit can be successfully applied to real research problems. In this study the reliability of the IFS theory for motif discovery tasks is proven.
Asunto(s)
Biología Computacional/métodos , Elementos Reguladores de la Transcripción , Análisis de Secuencia de ADN/métodos , Factores de Transcripción/metabolismo , Algoritmos , Secuencia de Bases , Sitios de Unión , ADN/química , ADN/metabolismo , Factores de Transcripción/químicaRESUMEN
Chronic infection with Hepatitis B virus (HBV) is a major risk factor for the development of advanced liver disease including fibrosis, cirrhosis, and hepatocellular carcinoma (HCC). The relative contribution of virological factors to disease progression has not been fully defined and tools aiding the deconvolution of complex patient virus profiles is an unmet clinical need. Variable viral mutant signatures develop within individual patients due to the low-fidelity replication of the viral polymerase creating 'quasispecies' populations. Here we present the first comprehensive survey of the diversity of HBV quasispecies through ultra-deep sequencing of the complete HBV genome across two distinct European and Asian patient populations. Seroconversion to the HBV e antigen (HBeAg) represents a critical clinical waymark in infected individuals. Using a machine learning approach, a model was developed to determine the viral variants that accurately classify HBeAg status. Serial surveys of patient quasispecies populations and advanced analytics will facilitate clinical decision support for chronic HBV infection and direct therapeutic strategies through improved patient stratification.
Asunto(s)
ADN Viral , Variación Genética , Genoma Viral , Virus de la Hepatitis B/genética , Hepatitis B Crónica/virología , Aprendizaje Automático , Carcinoma Hepatocelular/virología , Progresión de la Enfermedad , Femenino , Humanos , Neoplasias Hepáticas/virología , Masculino , CuasiespeciesRESUMEN
Global transcriptomic analysis via RNA-seq is often hampered by the high abundance of ribosomal (r)RNA in bacterial cells. To remove rRNA and enrich coding sequences, subtractive hybridization procedures have become the approach of choice prior to RNA-seq, with their efficiency varying in a manner dependent on sample type and composition. Yet, despite an increasing number of RNA-seq studies, comparative evaluation of bacterial rRNA depletion methods has remained limited. Moreover, no such study has utilized RNA derived from bacterial biofilms, which have potentially higher rRNA:mRNA ratios and higher rRNA carryover during RNA-seq analysis. Presently, we evaluated the efficiency of three subtractive hybridization-based kits in depleting rRNA from samples derived from biofilm, as well as planktonic cells of the opportunistic human pathogen Pseudomonas aeruginosa. Our results indicated different rRNA removal efficiency for the three procedures, with the Ribo-Zero kit yielding the highest degree of rRNA depletion, which translated into enhanced enrichment of non-rRNA transcripts and increased depth of RNA-seq coverage. The results indicated that, in addition to improving RNA-seq sensitivity, efficient rRNA removal enhanced detection of low abundance transcripts via qPCR. Finally, we demonstrate that the Ribo-Zero kit also exhibited the highest efficiency when P. aeruginosa/Staphylococcus aureus co-culture RNA samples were tested.
Asunto(s)
Bacterias/genética , Perfilación de la Expresión Génica/métodos , Análisis de Secuencia de ARN/métodos , Fenómenos Fisiológicos Bacterianos , Biopelículas , Regulación Bacteriana de la Expresión Génica , Pseudomonas aeruginosa/genética , Pseudomonas aeruginosa/fisiología , ARN Bacteriano , Sensibilidad y Especificidad , Staphylococcus aureusRESUMEN
Heteroaryldihydropyrimidine (HAP) and sulfamoylbenzamide (SBA) are promising non-nucleos(t)ide HBV replication inhibitors. HAPs are known to promote core protein mis-assembly, but the molecular mechanism of abnormal assembly is still elusive. Likewise, the assembly status of core protein induced by SBA remains unknown. Here we show that SBA, unlike HAP, does not promote core protein mis-assembly. Interestingly, two reference compounds HAP_R01 and SBA_R01 bind to the same pocket at the dimer-dimer interface in the crystal structures of core protein Y132A hexamer. The striking difference lies in a unique hydrophobic subpocket that is occupied by the thiazole group of HAP_R01, but is unperturbed by SBA_R01. Photoaffinity labeling confirms the HAP_R01 binding pose at the dimer-dimer interface on capsid and suggests a new mechanism of HAP-induced mis-assembly. Based on the common features in crystal structures we predict that T33 mutations generate similar susceptibility changes to both compounds. In contrast, mutations at positions in close contact with HAP-specific groups (P25A, P25S, or V124F) only reduce susceptibility to HAP_R01, but not to SBA_R01. Thus, HAP and SBA are likely to have distinctive resistance profiles. Notably, P25S and V124F substitutions exist in low-abundance quasispecies in treatment-naïve patients, suggesting potential clinical relevance.
Asunto(s)
Benzamidas/farmacología , Virus de la Hepatitis B/fisiología , Pirimidinas/farmacología , Replicación Viral/efectos de los fármacos , Antivirales/farmacología , Benzamidas/química , Sitios de Unión , Cápside/metabolismo , Cristalografía por Rayos X , Replicación del ADN/efectos de los fármacos , Células Hep G2 , Virus de la Hepatitis B/efectos de los fármacos , Humanos , Ligandos , Mutación/genética , Etiquetas de Fotoafinidad , Pirimidinas/química , Proteínas Virales/químicaRESUMEN
Mannose binding lectin (MBL) plays important role in the innate immunity of human. Mutations in the MBL2 gene can significantly change the serum level of MBL, and consequently alter the susceptibility and progression of infectious disease. However, the association between the MBL2 profile and the HBV mutation and quasispecies complexity has not yet been reported. Our approach includes the study of the MBL2 gene genotype as well as ultra-deep sequencing of the HBV viruses obtained from the plasma of 50 treatment naïve patients with chronic HBV infection. We found that the liver function was better among patients within the high MBL2 group with respect to those within the medium/low MBL2 group. Likewise, the number of mutations in the HBV X gene as well as the viral quasispecies complexity were significantly higher in medium/low MBL2 production group. Nucleotide substitution rates were also higher within the medium/low MBL2 production group in all positions described to have an influence in liver cancer development, except for A1499G. In this work we show that the MBL2 profile may have an impact on the HBV X gene mutations as well as on viral quasispecies complexity.
Asunto(s)
Genes Virales , Haplotipos , Virus de la Hepatitis B/genética , Hepatitis B/inmunología , Lectina de Unión a Manosa/genética , Mutación , Adulto , Anciano , Progresión de la Enfermedad , Femenino , Predisposición Genética a la Enfermedad , Hepatitis B/genética , Humanos , Inmunidad Innata , Masculino , Lectina de Unión a Manosa/sangre , Persona de Mediana Edad , Reacción en Cadena en Tiempo Real de la Polimerasa/métodos , Adulto JovenRESUMEN
Analysis of fusion transcripts has become increasingly important due to their link with cancer development. Since high-throughput sequencing approaches survey fusion events exhaustively, several computational methods for the detection of gene fusions from RNA-seq data have been developed. This kind of analysis, however, is complicated by native trans-splicing events, the splicing-induced complexity of the transcriptome and biases and artefacts introduced in experiments and data analysis. There are a number of tools available for the detection of fusions from RNA-seq data; however, certain differences in specificity and sensitivity between commonly used approaches have been found. The ability to detect gene fusions of different types, including isoform fusions and fusions involving non-coding regions, has not been thoroughly studied yet. Here, we propose a novel computational toolkit called InFusion for fusion gene detection from RNA-seq data. InFusion introduces several unique features, such as discovery of fusions involving intergenic regions, and detection of anti-sense transcription in chimeric RNAs based on strand-specificity. Our approach demonstrates superior detection accuracy on simulated data and several public RNA-seq datasets. This improved performance was also evident when evaluating data from RNA deep-sequencing of two well-established prostate cancer cell lines. InFusion identified 26 novel fusion events that were validated in vitro, including alternatively spliced gene fusion isoforms and chimeric transcripts that include intergenic regions. The toolkit is freely available to download from http:/bitbucket.org/kokonech/infusion.
Asunto(s)
Biología Computacional/métodos , Fusión Génica/genética , Secuenciación de Nucleótidos de Alto Rendimiento , Proteínas de Fusión Oncogénica/genética , Algoritmos , Humanos , Neoplasias/genética , Proteínas de Fusión Oncogénica/aislamiento & purificación , Análisis de Secuencia de ARN/métodos , Programas Informáticos , Transcriptoma/genéticaRESUMEN
Current anti-influenza therapy depends on administering drugs soon after infection, which is often impractical. We assessed whether combinations of oseltamivir (a neuraminidase inhibitor) and T-705 (a nonspecific inhibitor of viral polymerases) could extend the window for treating lethal infection with highly pathogenic A(H5N1) influenza virus in mice. Combination therapy protected 100% of mice, even when delayed until 96 h postinoculation. Compared to animals receiving monotherapy, mice receiving combination therapy had reduced viral loads and restricted viral spread in lung tissues, limited lung damage, and decreased inflammatory cytokine production. Next-generation sequencing showed that virus populations in T-705-treated mice had greater genetic variability, with more frequent transversion events, than did populations in control and oseltamivir-treated mice, but no substitutions associated with resistance to oseltamivir or T-705 were detected. Thus, combination therapy extended the treatment window for A(H5N1) influenza infection in mice and should be considered for evaluation in a clinical setting.
Asunto(s)
Amidas/farmacología , Subtipo H5N1 del Virus de la Influenza A , Infecciones por Orthomyxoviridae/tratamiento farmacológico , Oseltamivir/farmacología , Pirazinas/farmacología , Animales , Perros , Quimioterapia Combinada , Femenino , Células de Riñón Canino Madin Darby , Ratones , Ratones Endogámicos BALB C , Infecciones por Orthomyxoviridae/metabolismo , Infecciones por Orthomyxoviridae/patologíaRESUMEN
Infection with the human pathogen Helicobacter pylori (H. pylori) is a major risk factor for gastric cancer. Since the bacterium exerts multiple genotoxic effects, we examined the circumstances of DNA damage accumulation and identified regions within the host genome with high susceptibility to H. pylori-induced damage. Infection impaired several DNA repair factors, the extent of which depends on a functional cagPAI. This leads to accumulation of a unique DNA damage pattern, preferentially in transcribed regions and proximal to telomeres, in both gastric cell lines and primary gastric epithelial cells. The observed pattern correlates with focal amplifications in adenocarcinomas of the stomach and partly overlaps with known cancer genes. We thus demonstrate an impact of a bacterial infection directed toward specific host genomic regions and describe underlying characteristics that make such regions more likely to acquire heritable changes during infection, which could contribute to cellular transformation.
Asunto(s)
Daño del ADN , Mucosa Gástrica/metabolismo , Infecciones por Helicobacter/genética , Mucosa Gástrica/microbiología , Amplificación de Genes , Genoma Humano , Infecciones por Helicobacter/patología , Humanos , Proteínas Tirosina Quinasas Receptoras/genética , Proteínas Tirosina Quinasas Receptoras/metabolismo , Telómero/genéticaRESUMEN
Eukaryotic gene control regions are known to be spread throughout non-coding DNA sequences which may appear distant from the gene promoter. Transcription factors are proteins that coordinately bind to these regions at transcription factor binding sites to regulate gene expression. Several tools allow to detect significant co-occurrences of closely located binding sites (cis-regulatory modules, CRMs). However, these tools present at least one of the following limitations: 1) scope limited to promoter or conserved regions of the genome; 2) do not allow to identify combinations involving more than two motifs; 3) require prior information about target motifs. In this work we present CisMiner, a novel methodology to detect putative CRMs by means of a fuzzy itemset mining approach able to operate at genome-wide scale. CisMiner allows to perform a blind search of CRMs without any prior information about target CRMs nor limitation in the number of motifs. CisMiner tackles the combinatorial complexity of genome-wide cis-regulatory module extraction using a natural representation of motif combinations as itemsets and applying the Top-Down Fuzzy Frequent- Pattern Tree algorithm to identify significant itemsets. Fuzzy technology allows CisMiner to better handle the imprecision and noise inherent to regulatory processes. Results obtained for a set of well-known binding sites in the S. cerevisiae genome show that our method yields highly reliable predictions. Furthermore, CisMiner was also applied to putative in-silico predicted transcription factor binding sites to identify significant combinations in S. cerevisiae and D. melanogaster, proving that our approach can be further applied genome-wide to more complex genomes. CisMiner is freely accesible at: http://genome2.ugr.es/cisminer. CisMiner can be queried for the results presented in this work and can also perform a customized cis-regulatory module prediction on a query set of transcription factor binding sites provided by the user.