{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,7,1]],"date-time":"2026-07-01T14:44:35Z","timestamp":1782917075574,"version":"3.54.5"},"reference-count":37,"publisher":"MDPI AG","issue":"13","license":[{"start":{"date-parts":[[2026,7,1]],"date-time":"2026-07-01T00:00:00Z","timestamp":1782864000000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/501100012166","name":"National Key R&D Program of China","doi-asserted-by":"publisher","award":["2024YFC2309905"],"award-info":[{"award-number":["2024YFC2309905"]}],"id":[{"id":"10.13039\/501100012166","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["Diagnostics"],"abstract":"<jats:p>Background\/Objectives: RNA virus metagenomic sequencing is a core technology for emerging infectious disease prevention and control, as well as for rapid pathogen identification. However, two major bottlenecks hinder its clinical application: the low fraction of informative sequencing reads caused by host rRNA contamination, and insufficient viral genome coverage. This study aimed to optimize the experimental parameters of RNA virus metagenomic sequencing, address the above bottlenecks, and establish a standardized workflow. Methods: Forty-five clinically positive samples (20 influenza virus-positive; 25 SARS-CoV-2-positive) were investigated in three parallel comparative experiments: rRNA depletion versus no depletion; probe-mediated RNase H digestion versus rRNA blocking; and two fragmentation timing strategies (fragmentation before versus after reverse transcription). Sequencing was performed on the GeneMind platform, and key performance metrics were systematically analyzed. Results: Following rRNA depletion, the host sequence proportion in the influenza virus and SARS-CoV-2 samples decreased from 39.5 to 90.5% to 3.6 to 32.2%, while the 10\u00d7 genomic coverage increased from 0 to 99.4% to 98.1 to 100.0%. The proportion of host sequences captured by probe capture depletion (0.3\u201316.2%) was significantly (p &lt; 0.05) lower than that captured by rRNA blocking module (14.3\u201392.3%). No significant differences were observed in the 10\u00d7 genomic coverage (96.5\u2013100.0%) or the fraction of effective viral reads between the two fragmentation strategies (p &gt; 0.05). rRNA depletion is key to improving library quality, with post-capture probe digestion being optimal. Conclusions: The suggested optimization process will enhance sequencing efficiency and support the standardization of clinical RNA virus identification.<\/jats:p>","DOI":"10.3390\/diagnostics16132065","type":"journal-article","created":{"date-parts":[[2026,7,1]],"date-time":"2026-07-01T12:41:38Z","timestamp":1782909698000},"page":"2065","update-policy":"https:\/\/doi.org\/10.3390\/mdpi_crossmark_policy","source":"Crossref","is-referenced-by-count":0,"title":["Optimization of Metagenomic Library Construction for Influenza A Virus and SARS-CoV-2: Systematic Comparison of rRNA Depletion Strategies and Fragmentation Orders"],"prefix":"10.3390","volume":"16","author":[{"ORCID":"https:\/\/orcid.org\/0000-0001-7549-8286","authenticated-orcid":false,"given":"Yi","family":"Sun","sequence":"first","affiliation":[{"name":"Zhejiang Key Laboratory of Public Health Detection and Pathogenesis Research, Department of Microbiology, Zhejiang Provincial Center for Disease Control and Prevention, Hangzhou 310051, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Feng","family":"Wang","sequence":"additional","affiliation":[{"name":"Zhejiang Key Laboratory of Public Health Detection and Pathogenesis Research, Department of Microbiology, Zhejiang Provincial Center for Disease Control and Prevention, Hangzhou 310051, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Lingfeng","family":"Mao","sequence":"additional","affiliation":[{"name":"Hangzhou Baiyi Technology Co., Ltd., Hangzhou 310051, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Wenjun","family":"Lu","sequence":"additional","affiliation":[{"name":"Hangzhou Baiyi Technology Co., Ltd., Hangzhou 310051, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Hao","family":"Wu","sequence":"additional","affiliation":[{"name":"Hangzhou Baiyi Technology Co., Ltd., Hangzhou 310051, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Haiyan","family":"Mao","sequence":"additional","affiliation":[{"name":"Zhejiang Key Laboratory of Public Health Detection and Pathogenesis Research, Department of Microbiology, Zhejiang Provincial Center for Disease Control and Prevention, Hangzhou 310051, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Yanjun","family":"Zhang","sequence":"additional","affiliation":[{"name":"Zhejiang Key Laboratory of Public Health Detection and Pathogenesis Research, Department of Microbiology, Zhejiang Provincial Center for Disease Control and Prevention, Hangzhou 310051, China"}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"1968","published-online":{"date-parts":[[2026,7,1]]},"reference":[{"key":"ref_1","doi-asserted-by":"crossref","unstructured":"Luong, Q.X.T., Hoang, P.T., Ho, P.T., Ayun, R.Q., Lee, T.K., and Lee, S. (2025). Potential Broad-Spectrum Antiviral Agents: A Key Arsenal Against Newly Emerging and Reemerging Respiratory RNA Viruses. Int. J. Mol. Sci., 26.","DOI":"10.3390\/ijms26041481"},{"key":"ref_2","doi-asserted-by":"crossref","unstructured":"Paremskaia, A.I., Volchkov, P.Y., and Deviatkin, A.A. (2024). IAVCP (Influenza A Virus Consensus and Phylogeny): Automatic Identification of the Genomic Sequence of the Influenza A Virus from High-Throughput Sequencing Data. Viruses, 16.","DOI":"10.3390\/v16060873"},{"key":"ref_3","doi-asserted-by":"crossref","unstructured":"Gilbert-Girard, S., Piret, J., Carbonneau, J., H\u00e9naut, M., Goyette, N., and Boivin, G. (2024). Viral interference between severe acute respiratory syndrome coronavirus 2 and influenza A viruses. PLoS Pathog., 20.","DOI":"10.1101\/2024.02.02.578538"},{"key":"ref_4","doi-asserted-by":"crossref","first-page":"565","DOI":"10.1016\/S0140-6736(20)30251-8","article-title":"Genomic characterisation and epidemiology of 2019 novel coronavirus: Implications for virus origins and receptor binding","volume":"395","author":"Lu","year":"2020","journal-title":"Lancet"},{"key":"ref_5","doi-asserted-by":"crossref","first-page":"265","DOI":"10.1038\/s41586-020-2008-3","article-title":"A new coronavirus associated with human respiratory disease in China","volume":"579","author":"Wu","year":"2020","journal-title":"Nature"},{"key":"ref_6","doi-asserted-by":"crossref","first-page":"270","DOI":"10.1038\/s41586-020-2012-7","article-title":"A pneumonia outbreak associated with a new coronavirus of probable bat origin","volume":"579","author":"Zhou","year":"2020","journal-title":"Nature"},{"key":"ref_7","doi-asserted-by":"crossref","unstructured":"Heom, K.A., Wangsanuwat, C., Butkovich, L.V., Tam, S.C., Rowe, A.R., O\u2019Malley, M.A., and Dey, S.S. (2023). Targeted rRNA depletion enables efficient mRNA sequencing in diverse bacterial species and complex co-cultures. mSystems, 8.","DOI":"10.1128\/msystems.00281-23"},{"key":"ref_8","doi-asserted-by":"crossref","first-page":"16","DOI":"10.1016\/j.slast.2022.09.004","article-title":"Improving an rRNA depletion protocol with statistical design of experiments","volume":"28","author":"David","year":"2023","journal-title":"SLAS Technol."},{"key":"ref_9","doi-asserted-by":"crossref","first-page":"e5","DOI":"10.1093\/nar\/gkaa1072","article-title":"Optimized design of antisense oligomers for targeted rRNA depletion","volume":"49","author":"Phelps","year":"2021","journal-title":"Nucleic Acids Res."},{"key":"ref_10","doi-asserted-by":"crossref","first-page":"807","DOI":"10.1038\/nmeth.1507","article-title":"Validation of two ribosomal RNA removal methods for microbial metatranscriptomics","volume":"7","author":"He","year":"2010","journal-title":"Nat. Methods"},{"key":"ref_11","doi-asserted-by":"crossref","first-page":"i884","DOI":"10.1093\/bioinformatics\/bty560","article-title":"fastp: An ultra-fast all-in-one FASTQ preprocessor","volume":"34","author":"Chen","year":"2018","journal-title":"Bioinformatics"},{"key":"ref_12","unstructured":"Li, H. (2013). Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv."},{"key":"ref_13","doi-asserted-by":"crossref","first-page":"giab008","DOI":"10.1093\/gigascience\/giab008","article-title":"Twelve years of SAMtools and BCFtools","volume":"10","author":"Danecek","year":"2021","journal-title":"GigaScience"},{"key":"ref_14","doi-asserted-by":"crossref","first-page":"257","DOI":"10.1186\/s13059-019-1891-0","article-title":"Improved metagenomic analysis with Kraken 2","volume":"20","author":"Wood","year":"2019","journal-title":"Genome Biol."},{"key":"ref_15","doi-asserted-by":"crossref","first-page":"e104","DOI":"10.7717\/peerj-cs.104","article-title":"Bracken: Estimating species abundance in metagenomics data","volume":"3","author":"Lu","year":"2017","journal-title":"PeerJ Comput. Sci."},{"key":"ref_16","doi-asserted-by":"crossref","unstructured":"Tierney, L. (2012). The R Statistical Computing Environment. Proceedings of the Statistical Challenges in Modern Astronomy V, Springer.","DOI":"10.1007\/978-1-4614-3520-4_41"},{"key":"ref_17","doi-asserted-by":"crossref","first-page":"135","DOI":"10.1177\/096228029900800204","article-title":"Measuring agreement in method comparison studies","volume":"8","author":"Bland","year":"1999","journal-title":"Stat. Methods Med. Res."},{"key":"ref_18","doi-asserted-by":"crossref","first-page":"239","DOI":"10.1016\/j.jare.2019.12.005","article-title":"A reverse transcriptase-mediated ribosomal RNA depletion (RTR2D) strategy for the cost-effective construction of RNA sequencing libraries","volume":"24","author":"Zeng","year":"2020","journal-title":"J. Adv. Res."},{"key":"ref_19","doi-asserted-by":"crossref","first-page":"631","DOI":"10.1038\/s41576-019-0150-2","article-title":"RNA sequencing: The teenage years","volume":"20","author":"Stark","year":"2019","journal-title":"Nat. Rev. Genet."},{"key":"ref_20","doi-asserted-by":"crossref","first-page":"111","DOI":"10.1093\/ajcp\/aqac135","article-title":"Application of CRISPR-Based Human and Bacterial Ribosomal RNA Depletion for SARS-CoV-2 Shotgun Metagenomic Sequencing","volume":"159","author":"Clemons","year":"2023","journal-title":"Am. J. Clin. Pathol."},{"key":"ref_21","doi-asserted-by":"crossref","first-page":"12436","DOI":"10.1093\/nar\/gkaa1053","article-title":"Genome-wide mapping of SARS-CoV-2 RNA structures identifies therapeutically-relevant elements","volume":"48","author":"Manfredonia","year":"2020","journal-title":"Nucleic Acids Res."},{"key":"ref_22","doi-asserted-by":"crossref","first-page":"1128","DOI":"10.1038\/s41467-022-28603-2","article-title":"Secondary structural ensembles of the SARS-CoV-2 RNA genome in infected cells","volume":"13","author":"Lan","year":"2022","journal-title":"Nat. Commun."},{"key":"ref_23","unstructured":"Shaw, M., and Palese, P. (2007). Orthomyxoviridae: The viruses and their replication. Fields Virol., 1647\u20131689."},{"key":"ref_24","doi-asserted-by":"crossref","first-page":"425","DOI":"10.1038\/s41579-021-00542-7","article-title":"Influenza virus and SARS-CoV-2: Pathogenesis and host responses in the respiratory tract","volume":"19","author":"Flerlage","year":"2021","journal-title":"Nat. Rev. Microbiol."},{"key":"ref_25","doi-asserted-by":"crossref","first-page":"141","DOI":"10.1038\/s41579-020-00459-7","article-title":"Characteristics of SARS-CoV-2 and COVID-19","volume":"19","author":"Hu","year":"2021","journal-title":"Nat. Rev. Microbiol."},{"key":"ref_26","doi-asserted-by":"crossref","first-page":"105499","DOI":"10.1016\/j.jcv.2023.105499","article-title":"Comparison of the Ct-values for genomic and subgenomic SARS-CoV-2 RNA reveals limited predictive value for the presence of replication competent virus","volume":"165","author":"Roesmann","year":"2023","journal-title":"J. Clin. Virol."},{"key":"ref_27","doi-asserted-by":"crossref","first-page":"e02102-20","DOI":"10.1128\/JVI.02102-20","article-title":"Segment-Specific Kinetics of mRNA, cRNA, and vRNA Accumulation during Influenza Virus Infection","volume":"95","author":"Phan","year":"2021","journal-title":"J. Virol."},{"key":"ref_28","doi-asserted-by":"crossref","unstructured":"Daroch, N., Kannan, S.K., Srikantaiah, V., Mishra, R., and Ishtiaq, F. (2025). Respiratory virus infection dynamics and genomic surveillance to detect seasonal influenza subtypes in wastewater: A longitudinal study in Bengaluru, India. PLoS Glob. Public Health, 5.","DOI":"10.1371\/journal.pgph.0004640"},{"key":"ref_29","doi-asserted-by":"crossref","unstructured":"Herbert, Z.T., Kershner, J.P., Butty, V.L., Thimmapuram, J., Choudhari, S., Alekseyev, Y.O., Fan, J., Podnar, J.W., Wilcox, E., and Gipson, J. (2018). Cross-site comparison of ribosomal depletion kits for Illumina RNAseq library construction. BMC Genom., 19.","DOI":"10.1186\/s12864-018-4585-1"},{"key":"ref_30","doi-asserted-by":"crossref","unstructured":"Zhao, W., He, X., Hoadley, K.A., Parker, J.S., Hayes, D.N., and Perou, C.M. (2014). Comparison of RNA-Seq by poly (A) capture, ribosomal RNA depletion, and DNA microarray for expression profiling. BMC Genom., 15.","DOI":"10.1186\/1471-2164-15-419"},{"key":"ref_31","doi-asserted-by":"crossref","first-page":"12281","DOI":"10.1038\/s41598-019-48692-2","article-title":"Efficient and specific oligo-based depletion of rRNA","volume":"9","author":"Kraus","year":"2019","journal-title":"Sci. Rep."},{"key":"ref_32","doi-asserted-by":"crossref","first-page":"7","DOI":"10.1186\/s12575-023-00193-3","article-title":"Blocking Abundant RNA Transcripts by High-Affinity Oligonucleotides during Transcriptome Library Preparation","volume":"25","author":"Everaert","year":"2023","journal-title":"Biol. Proced. Online"},{"key":"ref_33","doi-asserted-by":"crossref","unstructured":"Wangsanuwat, C., Heom, K.A., Liu, E., O\u2019Malley, M.A., and Dey, S.S. (2020). Efficient and cost-effective bacterial mRNA sequencing from low input samples through ribosomal RNA depletion. BMC Genom., 21.","DOI":"10.1186\/s12864-020-07134-4"},{"key":"ref_34","doi-asserted-by":"crossref","first-page":"101","DOI":"10.1186\/s40249-025-01371-z","article-title":"SMART-RNA-Metavirome: A practical RNA metavirome platform compatible with high-throughput sequencing of both short and long reads","volume":"14","author":"Liu","year":"2025","journal-title":"Infect. Dis. Poverty"},{"key":"ref_35","doi-asserted-by":"crossref","first-page":"105695","DOI":"10.1016\/j.jcv.2024.105695","article-title":"Multicenter benchmarking of short and long read wet lab protocols for clinical viral metagenomics","volume":"173","author":"Huber","year":"2024","journal-title":"J. Clin. Virol."},{"key":"ref_36","doi-asserted-by":"crossref","first-page":"25","DOI":"10.1016\/j.ymeth.2013.03.009","article-title":"Zinc-mediated RNA fragmentation allows robust transcript reassembly upon whole transcriptome RNA-Seq","volume":"63","author":"Wery","year":"2013","journal-title":"Methods"},{"key":"ref_37","doi-asserted-by":"crossref","unstructured":"Puchta, M., Boczkowska, M., and Groszyk, J. (2020). Low RIN Value for RNA-Seq Library Construction from Long-Term Stored Seeds: A Case Study of Barley Seeds. Genes, 11.","DOI":"10.3390\/genes11101190"}],"container-title":["Diagnostics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/www.mdpi.com\/2075-4418\/16\/13\/2065\/pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,7,1]],"date-time":"2026-07-01T13:52:08Z","timestamp":1782913928000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.mdpi.com\/2075-4418\/16\/13\/2065"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2026,7,1]]},"references-count":37,"journal-issue":{"issue":"13","published-online":{"date-parts":[[2026,7]]}},"alternative-id":["diagnostics16132065"],"URL":"https:\/\/doi.org\/10.3390\/diagnostics16132065","relation":{},"ISSN":["2075-4418"],"issn-type":[{"value":"2075-4418","type":"electronic"}],"subject":[],"published":{"date-parts":[[2026,7,1]]}}}