The landscape of persistent human DNA viruses in femoral bone

Presence of herpesviruses, parvoviruses, and polyomaviruses in sinonasal lymphoma

PURPOSE

Sinonasal lymphoma (SL) is a rare lymphatic neoplasm of the nasal cavities, paranasal sinuses and nasopharynx. Whereas some risk factors for SL subtypes have been identified, their aetiology is unknown. Along with other predisposing factors, the viral association of lymphomas, such as Epstein-Barr virus (EBV) and Burkitt and Hodgkin lymphomas, is well-established. Modern molecular biology techniques have enabled the discovery of novel human viruses, exemplified by the protoparvovirus cutavirus (CuV), associated with cutaneous T-cell lymphoma. These findings, and the anatomical location of the sinonasal tract with its rich microbiome and infectious agents, justify in-depth studies among SL.

METHODS

We analysed the presence of 20 viruses of Orthoherpesviridae, Parvoviridae, and Polyomaviridae by qPCR in 24 SL tumours. We performed RNAscope in situ hybridisation (RISH) to localize the viruses. Parvovirus-specific IgG was analysed by enzyme immunoassay and targeted next-generation sequencing (NGS) was applied to detect CuV in plasma.

RESULTS

We detected viral DNA in 15/24 (63%) tumours; nine of EBV, six of human herpesvirus (HHV) -7, four each of HHV-6B and parvovirus B19, two of cytomegalovirus, and one each of CuV and Merkel-cell polyomavirus. We found tumours with up to four viruses per tumour, and localized CuV and EBV DNAs by RISH. Two of the ten plasma samples exhibited CuV IgG, and one plasma sample demonstrated CuV viremia by NGS.

CONCLUSION

Viruses were frequent findings in SL. The EBV detection rate was high in diffuse large B-cell lymphoma, and co-detections with other viruses were prevalent.

Hybrid-Capture Target Enrichment in Human Pathogens: Identification, Evolution, Biosurveillance, and Genomic Epidemiology

High-throughput sequencing (HTS) has revolutionised the field of pathogen genomics, enabling the direct recovery of pathogen genomes from clinical and environmental samples. However, pathogen nucleic acids are often overwhelmed by those of the host, requiring deep metagenomic sequencing to recover sufficient sequences for downstream analyses (e.g., identification and genome characterisation). To circumvent this, hybrid-capture target enrichment (HC) is able to enrich pathogen nucleic acids across multiple scales of divergences and taxa, depending on the panel used. In this review, we outline the applications of HC in human pathogens-bacteria, fungi, parasites and viruses-including identification, genomic epidemiology, antimicrobial resistance genotyping, and evolution. Importantly, we explored the applicability of HC to clinical metagenomics, which ultimately requires more work before it is a reliable and accurate tool for clinical diagnosis. Relatedly, the utility of HC was exemplified by COVID-19, which was used as a case study to illustrate the maturity of HC for recovering pathogen sequences. As we unravel the origins of COVID-19, zoonoses remain more relevant than ever. Therefore, the role of HC in biosurveillance studies is also highlighted in this review, which is critical in preparing us for the next pandemic. We also found that while HC is a popular tool to study viruses, it remains underutilised in parasites and fungi and, to a lesser extent, bacteria. Finally, weevaluated the future of HC with respect to bait design in the eukaryotic groups and the prospect of combining HC with long-read HTS.

Exploring the relationship between anellovirus load and clinical variables in hospitalized COVID-19 patients: Implications for immune activation and inflammation

Objectives

Anelloviruses have been linked with host-immunocompetence and inflammation. Here, we studied the anellovirus load in hospitalized COVID-19 patients.

Methods

We collected samples of patients recruited in the DAWN-Plasma trial that received convalescent plasma (CP) therapy (four plasma units) combined with standard of care (SOC) or SOC of alone. Plasma samples were collected on day 0 and 6 of hospitalization and we quantified anellovirus load. With multivariate models, clinical variables were associated with changes in anellovirus load.

Results

Samples were collected on day 0 and 6 of 150 patients (103 CP + SOC and 47 SOC). Anellovirus load was higher on day 0 compared to day 6 and we found a significant drop in SOC patients. Patients receiving immunosuppressive drug had a lower anellovirus load (coefficient: 1.021, 95% confidence interval [CI] 0.270-1.772, P = 0.008), while patients admitted to the emergency room displayed a higher abundance on day 0 (1.308, 95% CI 0.443-2.173, P = 0.003). Unspecific markers of inflammation and organ damage, D-dimer (0.001, 95% CI <0.001-0.001, P = 0.001) and lactate dehydrogenase (0.002, 95% CI 0.001-0.004, P = 0.044), were positively associated with anellovirus load. Finally, anellovirus load on day 0 (-39.9, 95% CI -75.72 to -4.27, P = 0.029) was negatively associated with SARS-CoV-2 antibody response on day.

Conclusion

The results showed associations between clinical variables and anellovirus load in COVID-19 patients. Many variables share properties related to host immunocompetence or inflammation. Therefore, we expect that anellovirus abundance displays the net state of immune activation.

Infectious disease in the Pleistocene: Old friends or old foes?

The impact of endemic and epidemic disease on humans has traditionally been seen as a comparatively recent historical phenomenon associated with the Neolithisation of human groups, an increase in population size led by sedentarism, and increasing contact with domesticated animals as well as species occupying opportunistic symbiotic and ectosymbiotic relationships with humans. The orthodox approach is that Neolithisation created the conditions for increasing population size able to support a reservoir of infectious disease sufficient to act as selective pressure. This orthodoxy is the result of an overly simplistic reliance on skeletal data assuming that no skeletal lesions equated to a healthy individual, underpinned by the assumption that hunter-gatherer groups were inherently healthy while agricultural groups acted as infectious disease reservoirs. The work of van Blerkom, Am. J. Phys. Anthropol., vol. suppl 37 (2003), Wolfe et al., Nature, vol. 447 (2007) and Houldcroft and Underdown, Am. J. Phys. Anthropol., vol. 160, (2016) has changed this landscape by arguing that humans and pathogens have long been fellow travelers. The package of infectious diseases experienced by our ancient ancestors may not be as dissimilar to modern infectious diseases as was once believed. The importance of DNA, from ancient and modern sources, to the study of the antiquity of infectious disease, and its role as a selective pressure cannot be overstated. Here we consider evidence of ancient epidemic and endemic infectious diseases with inferences from modern and ancient human and hominin DNA, and from circulating and extinct pathogen genomes. We argue that the pandemics of the past are a vital tool to unlock the weapons needed to fight pandemics of the future.

The epidemiological and genetic characteristics of human parvovirus B19 in patients with febrile rash illnesses in China

To understand the epidemiological and genetic characteristics of B19V, a multiple-province surveillance of patients with febrile rash illnesses (FRIs) were conducted in China during 2009 ~ 2021. The clinical specimens of 3,820 FRI patients were collected and tested for B19V DNA. A total of 99 (2.59%) patients were positive for B19V, and 49 (49.49%) were children under 5 years old. B19V infections occurred throughout the year without obvious seasonal pattern. Ten NS1-VP1u sequences and seven genome sequences were obtained in this study, identified as subgenotype 1a. Combined with the globally representative genome sequences, no temporal and geographic clustering trends of B19V were observed, and there was no significant correlation between B19V sequences and clinical manifestations. The evolutionary rate of the B19V genome was 2.30 × 10-4 substitutions/site/year. The number of negative selection sites was higher than that of positive selection sites. It was the first to comprehensively describe the prevalence patterns and evolutionary characteristics of B19V in FRI patients in China. B19V played the role in FRI patients. Children under 5 years old were the main population of B19V infection. Subgenotype 1a was prevalent in FRI patients in China. B19V showed a high mutation rate, while negative selection acted on the genome.

Characterization of the Human Blood Virome in Iranian Multiple Transfused Patients

Blood transfusion safety is an essential element of public health. Current blood screening strategies rely on targeted techniques that could miss unknown or unexpected pathogens. Recent studies have demonstrated the presence of a viral community (virobiota/virome) in the blood of healthy individuals. Here, we characterized the blood virome in patients frequently exposed to blood transfusion by using Illumina metagenomic sequencing. The virome of these patients was compared to viruses present in healthy blood donors. A total number of 155 beta-thalassemia, 149 hemodialysis, and 100 healthy blood donors were pooled with five samples per pool. Members of the Anelloviridae and Flaviviridae family were most frequently observed. Interestingly, samples of healthy blood donors harbored traces of potentially pathogenic viruses, including adeno-, rota-, and Merkel cell polyomavirus. Viruses of the Anelloviridae family were most abundant in the blood of hemodialysis patients and displayed a higher anellovirus richness. Pegiviruses (Flaviviridae) were only observed in patient populations. An overall trend of higher eukaryotic read abundance in both patient groups was observed. This might be associated with increased exposure through blood transfusion. Overall, the findings in this study demonstrated the presence of various viruses in the blood of Iranian multiple-transfused patients and healthy blood donors.

Herpesviruses, polyomaviruses, parvoviruses, papillomaviruses, and anelloviruses in vestibular schwannoma

Etiology of vestibular schwannoma (VS) is unknown. Viruses can infect and reside in neural tissues for decades, and new viruses with unknown tumorigenic potential have been discovered. The presence of herpesvirus, polyomavirus, parvovirus, and anellovirus DNA was analyzed by quantitative PCR in 46 formalin-fixed paraffin-embedded VS samples. Five samples were analyzed by targeted next-generation sequencing. Viral DNA was detected altogether in 24/46 (52%) tumor samples, mostly representing anelloviruses (46%). Our findings show frequent persistence of anelloviruses, considered normal virome, in VS. None of the other viruses showed an extensive presence, thereby suggesting insignificant role in VS.

The complexity landscape of viral genomes

BACKGROUND

Viruses are among the shortest yet highly abundant species that harbor minimal instructions to infect cells, adapt, multiply, and exist. However, with the current substantial availability of viral genome sequences, the scientific repertory lacks a complexity landscape that automatically enlights viral genomes' organization, relation, and fundamental characteristics.

RESULTS

This work provides a comprehensive landscape of the viral genome's complexity (or quantity of information), identifying the most redundant and complex groups regarding their genome sequence while providing their distribution and characteristics at a large and local scale. Moreover, we identify and quantify inverted repeats abundance in viral genomes. For this purpose, we measure the sequence complexity of each available viral genome using data compression, demonstrating that adequate data compressors can efficiently quantify the complexity of viral genome sequences, including subsequences better represented by algorithmic sources (e.g., inverted repeats). Using a state-of-the-art genomic compressor on an extensive viral genomes database, we show that double-stranded DNA viruses are, on average, the most redundant viruses while single-stranded DNA viruses are the least. Contrarily, double-stranded RNA viruses show a lower redundancy relative to single-stranded RNA. Furthermore, we extend the ability of data compressors to quantify local complexity (or information content) in viral genomes using complexity profiles, unprecedently providing a direct complexity analysis of human herpesviruses. We also conceive a features-based classification methodology that can accurately distinguish viral genomes at different taxonomic levels without direct comparisons between sequences. This methodology combines data compression with simple measures such as GC-content percentage and sequence length, followed by machine learning classifiers.

CONCLUSIONS

This article presents methodologies and findings that are highly relevant for understanding the patterns of similarity and singularity between viral groups, opening new frontiers for studying viral genomes' organization while depicting the complexity trends and classification components of these genomes at different taxonomic levels. The whole study is supported by an extensive website (https://asilab.github.io/canvas/) for comprehending the viral genome characterization using dynamic and interactive approaches.

Torque Teno Virus Primary Infection Kinetics in Early Childhood

Human torque teno viruses (TTVs) are a diverse group of small nonenveloped viruses with circular, single-stranded DNA genomes. These elusive anelloviruses are harbored in the blood stream of most humans and have thus been considered part of the normal flora. Whether the primary infection as a rule take(s) place before or after birth has been debated. The aim of our study was to determine the time of TTV primary infection and the viral load and strain variations during infancy and follow-up for up to 7 years. TTV DNAs were quantified in serial serum samples from 102 children by a pan-TTV quantitative PCR, and the amplicons from representative time points were cloned and sequenced to disclose the TTV strain diversity. We detected an unequivocal rise in TTV-DNA prevalence, from 39% at 4 months of age to 93% at 2 years; all children but one, 99%, became TTV-DNA positive before age 4 years. The TTV-DNA quantities ranged from 5 × 10¹ to 4 × 10⁷ copies/mL, both within and between the children. In conclusion, TTV primary infections occur mainly after birth, and increase during the first two years with high intra- and interindividual variation in both DNA quantities and virus strains.

Detection of Low-Copy Human Virus DNA upon Prolonged Formalin Fixation

Formalin fixation, albeit an outstanding method for morphological and molecular preservation, induces DNA damage and cross-linking, which can hinder nucleic acid screening. This is of particular concern in the detection of low-abundance targets, such as persistent DNA viruses. In the present study, we evaluated the analytical sensitivity of viral detection in lung, liver, and kidney specimens from four deceased individuals. The samples were either frozen or incubated in formalin (±paraffin embedding) for up to 10 days. We tested two DNA extraction protocols for the control of efficient yields and viral detections. We used short-amplicon qPCRs (63–159 nucleotides) to detect 11 DNA viruses, as well as hybridization capture of these plus 27 additional ones, followed by deep sequencing. We observed marginally higher ratios of amplifiable DNA and scantly higher viral genoprevalences in the samples extracted with the FFPE dedicated protocol. Based on the findings in the frozen samples, most viruses were detected regardless of the extended fixation times. False-negative calls, particularly by qPCR, correlated with low levels of viral DNA (<250 copies/million cells) and longer PCR amplicons (>150 base pairs). Our data suggest that low-copy viral DNAs can be satisfactorily investigated from FFPE specimens, and encourages further examination of historical materials.

Survey of Viral Reactivations in Elite Athletes: A Case-Control Study

Exercise-induced immune perturbations have been proposed to increase susceptibility to viral infections. We investigated the replication of persisting viruses as indicators of immune function in elite cross-country skiers after ten months of sustained high-performance exercise. The viruses evaluated, nine human herpesviruses (HHVs) and torque teno virus (TTV), are typically restrained in health but replicate actively in immunosuppressed individuals. We collected sera from 27 Finnish elite cross-country skiers at the end of the competition's season and 27 matched controls who perform moderate exercise. We quantified all the HHVs and-TTV via highly sensitive qPCRs. To verify equal past exposures between the groups, we assessed the IgG antibody prevalences toward HHV-4 (Epstein-Barr virus, EBV) and HHV-5 (human cytomegalovirus, HCMV). We found equal TTV DNA prevalences in athletes (63%) and controls (63%) and loads with respective geometric means of 1.7 × 103 and 1.2 × 103 copies/mL of serum. Overall, the copy numbers were low and consistent with those of healthy individuals. Neither of the groups presented with herpesvirus viremia despite similar past exposures to HHVs (seroprevalences of EBV 70% vs. 78% and HCMV 52% vs. 44% in athletes and controls, respectively). We found no evidence of increased replication of persistent viruses in elite athletes, arguing against impaired viral immunity due to high-performance exercise.

The Human Bone Marrow Is Host to the DNAs of Several Viruses

The long-term impact of viruses residing in the human bone marrow (BM) remains unexplored. However, chronic inflammatory processes driven by single or multiple viruses could significantly alter hematopoiesis and immune function. We performed a systematic analysis of the DNAs of 38 viruses in the BM. We detected, by quantitative PCRs and next-generation sequencing, viral DNA in 88.9% of the samples, up to five viruses in one individual. Included were, among others, several herpesviruses, hepatitis B virus, Merkel cell polyomavirus and, unprecedentedly, human papillomavirus 31. Given the reactivation and/or oncogenic potential of these viruses, their repercussion on hematopoietic and malignant disorders calls for careful examination. Furthermore, the implications of persistent infections on the engraftment, regenerative capacity, and outcomes of bone marrow transplantation deserve in-depth evaluation.

A semi-automatic methodology for analysing distributed and private biobanks

Privacy issues limit the analysis and cross-exploration of most distributed and private biobanks, often raised by the multiple dimensionality and sensitivity of the data associated with access restrictions and policies. These characteristics prevent collaboration between entities, constituting a barrier to emergent personalized and public health challenges, namely the discovery of new druggable targets, identification of disease-causing genetic variants, or the study of rare diseases. In this paper, we propose a semi-automatic methodology for the analysis of distributed and private biobanks. The strategies involved in the proposed methodology efficiently enable the creation and execution of unified genomic studies using distributed repositories, without compromising the information present in the datasets. We apply the methodology to a case study in the current Covid-19, ensuring the combination of the diagnostics from multiple entities while maintaining privacy through a completely identical procedure. Moreover, we show that the methodology follows a simple, intuitive, and practical scheme.

Efficient DNA sequence compression with neural networks

BACKGROUND

The increasing production of genomic data has led to an intensified need for models that can cope efficiently with the lossless compression of DNA sequences. Important applications include long-term storage and compression-based data analysis. In the literature, only a few recent articles propose the use of neural networks for DNA sequence compression. However, they fall short when compared with specific DNA compression tools, such as GeCo2. This limitation is due to the absence of models specifically designed for DNA sequences. In this work, we combine the power of neural networks with specific DNA models. For this purpose, we created GeCo3, a new genomic sequence compressor that uses neural networks for mixing multiple context and substitution-tolerant context models.

FINDINGS

We benchmark GeCo3 as a reference-free DNA compressor in 5 datasets, including a balanced and comprehensive dataset of DNA sequences, the Y-chromosome and human mitogenome, 2 compilations of archaeal and virus genomes, 4 whole genomes, and 2 collections of FASTQ data of a human virome and ancient DNA. GeCo3 achieves a solid improvement in compression over the previous version (GeCo2) of $2.4\%$, $7.1\%$, $6.1\%$, $5.8\%$, and $6.0\%$, respectively. To test its performance as a reference-based DNA compressor, we benchmark GeCo3 in 4 datasets constituted by the pairwise compression of the chromosomes of the genomes of several primates. GeCo3 improves the compression in $12.4\%$, $11.7\%$, $10.8\%$, and $10.1\%$ over the state of the art. The cost of this compression improvement is some additional computational time (1.7-3 times slower than GeCo2). The RAM use is constant, and the tool scales efficiently, independently of the sequence size. Overall, these values outperform the state of the art.

CONCLUSIONS

GeCo3 is a genomic sequence compressor with a neural network mixing approach that provides additional gains over top specific genomic compressors. The proposed mixing method is portable, requiring only the probabilities of the models as inputs, providing easy adaptation to other data compressors or compression-based data analysis tools. GeCo3 is released under GPLv3 and is available for free download at https://github.com/cobilab/geco3.