HTSinfer: inferring metadata from bulk Illumina RNA-Seq libraries
Language English Country England, Great Britain Media print
Document type Journal Article
PubMed
39969909
PubMed Central
PMC11889452
DOI
10.1093/bioinformatics/btaf076
PII: 8024254
Knihovny.cz E-resources
- MeSH
- Gene Library MeSH
- Metadata * MeSH
- Sequence Analysis, RNA * methods MeSH
- RNA-Seq * methods MeSH
- Software * MeSH
- Computational Biology methods MeSH
- High-Throughput Nucleotide Sequencing * methods MeSH
- Publication type
- Journal Article MeSH
SUMMARY: The Sequencing Read Archive is one of the largest and fastest-growing repositories of sequencing data, containing tens of petabytes of sequenced reads. Its data is used by a wide scientific community, often beyond the primary study that generated them. Such analyses rely on accurate metadata concerning the type of experiment and library, as well as the organism from which the sequenced reads were derived. These metadata are typically entered manually by contributors in an error-prone process, and are frequently incomplete. In addition, easy-to-use computational tools that verify the consistency and completeness of metadata describing the libraries to facilitate data reuse, are largely unavailable. Here, we introduce HTSinfer, a Python-based tool to infer metadata directly and solely from bulk RNA-sequencing data generated on Illumina platforms. HTSinfer leverages genome sequence information and diagnostic genes to rapidly and accurately infer the library source and library type, as well as the relative read orientation, 3' adapter sequence and read length statistics. HTSinfer is written in a modular manner, published under a permissible free and open-source license and encourages contributions by the community, enabling easy addition of new functionalities, e.g. for the inference of additional metrics, or the support of different experiment types or sequencing platforms. AVAILABILITY AND IMPLEMENTATION: HTSinfer is released under the Apache License 2.0. Latest code is available via GitHub at https://github.com/zavolanlab/htsinfer, while releases are published on Bioconda. A snapshot of the HTSinfer version described in this article was deposited at Zenodo at 10.5281/zenodo.13985958.
Biozentrum University of Basel Basel 4056 Switzerland
Institute of Computer Science Masaryk University Brno 60200 Czech Republic
See more in PubMed
Bai J, Bandla C, Guo J et al. BioContainers registry: Searching bioinformatics and proteomics tools, packages, and containers. J Proteome Res 2021;20:2056–61. PubMed PMC
Balajti M, Kandari R, Jurič B et al. HTSinfer: Supplementary materials. Zenodo. 2024a. 10.5281/zenodo.14177729
Balajti M, Kanitz A, Kandhari R et al. zavolanlab/htsinfer: v1.0.0-rc.1. Zenodo. 2024b. 10.5281/zenodo.13985958
Ban N, Beckmann R, Cate JHD et al. A new system for naming ribosomal proteins. Curr Opin Struct Biol 2014;24:165–9. PubMed PMC
Bray NL, Pimentel H, Melsted P et al. Near-optimal probabilistic RNA-seq quantification. Nat Biotechnol 2016;34:525–7. PubMed
Dobin A, Davis CA, Schlesinger F et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics 2013;29:15–21. PubMed PMC
European Organization For Nuclear Research. OpenAIRE. Zenodo. 2013. https://doi.org/10.25495/7GXK-RD71
Grüning B, Dale R, Sjödin A et al. Bioconda: sustainable and comprehensive software distribution for the life sciences. Nat Methods 2018;15:475–6. PubMed PMC
Harrison PW, Amode MR, Austine-Orimoloye O et al. Ensembl 2024. Nucleic Acids Res 2024;52:D891–9. PubMed PMC
Jiménez RC, Kuzak M, Alhamdoosh M et al. Four simple recommendations to encourage best practices in research software. F1000Res 2017;6:876. PubMed PMC
Katsantoni M, Gypas F, Herrmann CJ et al. ZARP: A user-friendly and versatile RNA-seq analysis workflow. F1000Res 2024;13:533.
Katz KS, Shutov O, Lapoint R et al. STAT: A fast, scalable, minhash-based k-mer tool to assess Sequence Read Archive next-generation sequence submissions. Genome Biol 2021;22:270. PubMed PMC
Katz Y, Wang ET, Airoldi EM et al. Analysis and design of RNA sequencing experiments for identifying isoform regulation. Nat Methods 2010;7:1009–15. PubMed PMC
Li H, Handsaker B, Wysoker A et al. The sequence alignment/map format and samtools. Bioinformatics 2009;25:2078–9. PubMed PMC
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet j 2011;17:10.
Patro R, Duggal G, Love MI et al. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods 2017;14:417–9. PubMed PMC
Petibon C, Malik Ghulam M, Catala M et al. Regulation of ribosomal protein genes: An ordered anarchy. Wiley Interdiscip Rev RNA 2021;12:e1632. PubMed PMC
Sayers EW, Bolton EE, Brister JR et al. Database resources of the national center for biotechnology information. Nucleic Acids Res 2022;50:D20–6. PubMed PMC
Scarpin MR, Busche M, Martinez RE et al. An updated nomenclature for plant ribosomal protein genes. Plant Cell 2023;35:640–3. PubMed PMC
Thorrez L, Van Deun K, Tranchevent L-C et al. Using ribosomal protein genes as reference: A tale of caution. PLoS One 2008;3:e1854. PubMed PMC
Wang L, Nie J, Sicotte H et al. Measure transcript integrity using RNA-seq data. BMC Bioinformatics 2016;17:58. PubMed PMC
Wik E, Olin H, Haughey C et al. GUESSmyLT: Software to guess the RNA-seq library type of paired and single end read files. JOSS 2019;4:1344.