Distinct fragmentation patterns of circulating viral cell-free DNA in 83,552 non-invasive prenatal testing samples

Abstract

Aim: The fragmentation characteristics of cell-free DNA (cfDNA) are informative biomarkers in liquid biopsies, including non-invasive prenatal testing (NIPT), as they provide insights into the origins of the cfDNA. Viral infections by DNA viruses can contribute to the available cfDNA in these samples. Here, we characterize the fragment size distribution of viral cfDNA fragments obtained from available anonymous NIPT samples.

Methods: A viral database of 224 DNA viruses was generated from the NCBI RefSeq viral database. Paired-end cfDNA sequencing reads from 83.522 NIPT samples that did not map to any of the human chromosomes, or mitochondrial DNA of the human reference genome build GRCh38 (excluding alternative and unplaced contigs) were remapped to the generated viral database. Reads mapping to the 14 most abundant DNA viruses were selected, and fragment size distributions were analyzed in detail.

Results: Distinct fragmentation patterns were identified for several DNA viruses, most likely due to differences in viral tropism, chromatinization (binding of nucleosomes), and the topology of the viral DNA. In high viral load parvo B19 positive samples, the fragment size distribution differed between samples, potentially reflecting the state of the infection.

Conclusion: These findings outline the potential for liquid biopsies to elucidate the dynamics behind the viral infection, which may potentially have various clinical applications. Our data provide preliminary insights on the use of fragmentomics of viral cfDNA to distinguish between reactivation, reinfection, and primary infection and monitoring the state of viral infections.

Keywords

Virome / fragmentomics / DNA degradation / epigenetics / viral cell-free DNA / NIPT

References

[1] Lo YMD,Jiang P.Epigenetics, fragmentomics, and topology of cell-free DNA in liquid biopsies.Science2021;372:eaaw3616

[2] Liu S,Chen F.Genomic analyses from non-invasive prenatal testing reveal genetic associations, patterns of viral infections, and Chinese population history.Cell2018;175:347-359.e14

[3] Lo YM,Sun H.Maternal plasma DNA sequencing reveals the genome-wide genetic and mutational profile of the fetus.Sci Transl Med2010;2:61ra91

[4] Kohabir K,Sistermans EA.Fragmentomic cfDNA patterns in noninvasive prenatal testing and beyond.J Biomed Transl Res2021;7:38-47

[5] Yu SC,Zheng YW.Size-based molecular diagnostics using plasma DNA for noninvasive prenatal testing.Proc Natl Acad Sci U S A2014;111:8583-8 PMCID:PMC4060699

[6] der Pol Y, Mouliere F. Toward the early detection of cancer by decoding the epigenetic and environmental fingerprints of cell-free DNA.Cancer Cell2019;36:350-68

[7] Tsai K.Epigenetic and epitranscriptomic regulation of viral replication.Nat Rev Microbiol2020;18:559-70 PMCID:PMC7291935

[8] Zalckvar E,Tillo D.Nucleosome maps of the human cytomegalovirus genome reveal a temporal switch in chromatin organization linked to a major IE protein.Proc Natl Acad Sci U S A2013;110:13126-31 PMCID:PMC3740854

[9] Oh J,Chen EZ.Genome wide nucleosome mapping for HSV-1 shows nucleosomes are deposited at preferred positions during lytic infection.PLoS One2015;10:e0117471 PMCID:PMC4339549

[10] Backer MW, Russcher A, Kroes AC, Koppelman MH, Lanfermeijer M, Zaaijer HL. Detection of parvovirus B19 DNA in blood: viruses or DNA remnants?.J Clin Virol2016;84:19-23

[11] Ryan JL,Swinnen LJ.Epstein-Barr virus (EBV) DNA in plasma is not encapsidated in patients with EBV-related malignancies.Diagn Mol Pathol2004;13:61-8

[12] Lam WKJ,Chan KCA.Sequencing-based counting and size profiling of plasma Epstein-Barr virus DNA enhance population screening of nasopharyngeal carcinoma.Proc Natl Acad Sci U S A2018;115:E5115-24 PMCID:PMC5984543

[13] Lam WKJ,Chan KCA.Methylation analysis of plasma DNA informs etiologies of Epstein-Barr virus-associated diseases.Nat Commun2019;10:3256 PMCID:PMC6646310

[14] Peddu V,Casto AM.High-resolution profiling of human cytomegalovirus cell-free DNA in human plasma highlights its exceptionally fragmented nature.Sci Rep2020;10:3734 PMCID:PMC7048871

[15] Bailey TL,Buske FA.MEME SUITE: tools for motif discovery and searching.Nucleic Acids Res2009;37:W202-8 PMCID:PMC2703892

[16] Huerta-Cepas J,Bork P.ETE 3: reconstruction, analysis, and visualization of phylogenomic data.Mol Biol Evol2016;33:1635-8 PMCID:PMC4868116

[17] Ondov BD,Melsted P.Mash: fast genome and metagenome distance estimation using MinHash.Genome Biol2016;17:132 PMCID:PMC4915045

[18] Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. Available from: https://arxiv.org/abs/1303.3997. [Last accessed on 28 Sep 2021]

[19] Li H,Wysoker A.1000 Genome Project Data Processing SubgroupThe sequence alignment/map format and SAMtools.Bioinformatics2009;25:2078-9 PMCID:PMC2723002

[20] Li H.Fast and accurate long-read alignment with Burrows-Wheeler transform.Bioinformatics2010;26:589-95 PMCID:PMC2828108

[21] Faust GG.SAMBLASTER: fast duplicate marking and structural variant read extraction.Bioinformatics2014;30:2503-5 PMCID:PMC4147885

[22] Kim SK,Geis J.Determination of fetal DNA fraction from the plasma of pregnant women using sequence read counts.Prenat Diagn2015;35:810-5

[23] Fritz M, Leinonen R, Cochrane G, Birney E. Efficient storage of high throughput DNA sequencing data using reference-based compression.Genome Res2011;21:734-40 PMCID:PMC3083090

[24] Pantry SN.Latency, Integration, and reactivation of human herpesvirus-6.Viruses2017;9:194 PMCID:PMC5537686

[25] Cook L,Boonyaratanakornkit J,Sam SS.Does size matter?.J Clin Microbiol2018;56:e01061-18 PMCID:PMC6258844