Amadeus T. Heinz1,2
Silke Grumaz3
Christoph Slavetinsky4Michaela Döring1Manon Queudeville5
Rupert Handgretinger6
Martin Ebinger1*- 1Department of Pediatric Hematology and Oncology, University Children´s Hospital Tuebingen, Tuebingen, Germany
- 2Stuttgart Cancer Center, Zentrum für Kinder-, Jugend- und Frauenmedizin (Olgahospital), Pädiatrie 5 (Pädiatrische Onkologie, Hämatologie, Immunologie), Klinikum der Landeshauptstadt Stuttgart, Stuttgart, Germany
- 3Noscendo GmbH, Duisburg, Germany
- 4Department of Pediatric Surgery and Urology, University Children´s Hospital Tuebingen, Tuebingen, Germany
- 5Department for Stem Cell Transplantation and Immunology, Klinikum Hamburg-Eppendorf (UKE), Hamburg, Germany
- 6Abu Dhabi Stem Cells Center, Abu Dhabi, United Arab Emirates
The etiology of pediatric acute lymphatic leukemia (ALL) is still unclear. Whole-metagenome shotgun sequencing of bone marrow samples in patients with treatment-naïve ALL (n=6) was performed for untargeted investigation of bacterial and viral DNA. The control group consisted of healthy children (n=4) and children with non-oncologic diseases (n=2) undergoing bone marrow sampling. Peripheral blood of all participants was investigated at the same time. After bioinformatical elimination of potential contaminants by comparison with the employed controls, no significant amounts of microbial or viral DNA were identified.
1 Introduction
Acute lymphoblastic leukemia (ALL) is the most common childhood cancer (Pui, 1995). Nevertheless, its etiology is still unknown. Derived from epidemiological studies, infection-related hypotheses have frequently been discussed in the past (Greaves and Alexander, 1993; Kinlen et al., 1995). However, no specific causative agent has been identified, so far (O’Connor and Boneva, 2007). Using next-generation sequencing (NGS) techniques for untargeted amplification of the whole viral and bacterial metagenome in pre-treatment bone marrow samples, a higher viral load has been identified in pediatric patients with ALL compared to acute myeloid leukemia (AML) (Francis et al., 2014). To further investigate the infection-related hypothesis in pediatric ALL using novel sequencing techniques, we performed a pilot trial to search for microbial and viral DNA using whole metagenome sequencing (WMS) of peripheral blood and bone marrow samples of six patients with treatment-naïve ALL and six children without oncologic diseases acting as negative controls.
2 Methods
2.1 Collection of samples
All samples were collected during planned bone marrow punctures of the iliac crest after disinfection of the skin and under sterile conditions using an aspiration cannula. Bone marrow punctures were scheduled in six patients for initial diagnosis of ALL, for bone marrow donation in four healthy participants and for other reasons in two participants (diagnostic purpose in one patient with severe aplastic anemia (SAA) and gathering of an autologous backup ahead of bone marrow transplantation in one patient suffering from sickle cell anemia). All work was conducted within the formal approval of the institutional ethics comittee (University Tuebingen, approval number 561/2020BO1, gathered 11/2020) and in accordance with the Declaration of Helsinki. Blood and bone marrow samples were collected after obtaining written informed consent of participant’s guardians and informing study participants according to the corresponding age group.
The sample collection for this study was performed in the midst or the end of the bone marrow aspiration, not at the beginning, to minimize the risk of contamination through commensal bacteria of the skin. A total of five milliliters of bone marrow was drawn, plus one additional milliliter of bone marrow for microbial cultivation as contamination control. Five milliliters of peripheral blood were collected through a peripheral venous blood draw or a central venous catheter (Hickman line). In patients below three years of age, two milliliters of blood were taken. Collection of bone marrow and peripheral blood was performed simultaneously, whenever possible.
The collected sample material was placed in DNA stabilizer tubes (Streck Cell-free DNA BCT, CE-IVD, Streck, La Vista, NE, USA) and sent to the laboratory at room temperature. After centrifugation for separation of plasma and peripheral blood cells, the plasma was stored at -80°C until further processing. Bone marrow samples were treated in the same manner.
2.2 Metagenomic next-generation sequencing
Blood samples and bone marrow samples were processed in a blinded manner and therefore processed identically. Cells were separated from the supernatant by centrifugation at 1,600 x g for 10 min at 4°C and the supernatant was transferred to a fresh reaction tube. Then, a second centrifugation step at 16,000 x g for 10 min at 4°C was performed, supernatants were again transferred, and aliquots were further stored at -80°C. Nucleic acid isolation, quality controls and library preparation were carried out as previously described (Grumaz et al., 2016). NGS was performed on a Illumina © NextSeq 550 instrument. Adequate positive and negative controls accompanied all laboratory and sequencing procedures. Raw sequencing data were subjected to various QC controls comprising PHRED-Score filtering, adapter trimming, complexity filtering as well as k-mer based classification and contamination screening. To pass the quality filter, read quality needed to surpass a PHRED score of 20 and achieve a minimal length of 50 bp after quality control. All data generated were analyzed using Noscendo’s DISQVER platform, integrating the CE-IVD for pathogen detection assay from blood. The DISQVER platform comprises a curated microbial genome reference database of over 16,000 microbial species covering >1,500 pathogens and can detect bacteria, DNA-viruses, fungi, and parasites while differentiating contaminating commensals from infective agents via run controls and dedicated control collectives. For each pathogen type (bacterium, virus, fungus, parasite) a dedicated minimum threshold for reporting was set, that is further refined by incorporating data for each species individually regarding their observation in control cohorts. The limit of detection (LOD) was determined by bioinformatic and wetlab analytical validation procedures, which are further described in Supplementary Material 1. For the PCA analysis, raw sequencing data after removal of human sequences and after classification of the non-human fraction was used for each sample. Read counts were log10 transformed prior to PCA analysis. The positive control was used after bioinformatic exclusion of the known, spiked-in microbial sequences, in order to not bias the results by taking these into account.
2.3 Statistical methods
Statistical analyses were conducted using SPSS 27 (Armonk, New York) and RStudio 2022.12.0 (Boston, MA). Descriptive statistics are reported as median [range], if not otherwise specified. For comparison of means between cfDNA levels in control and ALL cohorts, the Wilcoxon rank sum test was used as values were not-normally distributed. For comparison of read mapping proportions between the different groups, the Kruskal-Wallis test was used.
3 Results
3.1 Participant characteristics
As shown in Table 1, participants were well balanced concerning gender and age. Common type ALL was the most frequent diagnosis in the ALL group, whereas controls were mainly healthy children and young adolescents willing to donate bone marrow for a sibling.
Table 1 Participant characteristics.
In the ALL group, four out of six patients received antibiotic treatment ahead of BMA. In three patients, intravenous Piperacillin/Tazobactam or Cefuroxime plus Clindamycin was administered just one day ahead of BMA. The remaining case with fever and mild pancytopenia was treated intravenously for six days with Cefuroxime and then discharged. Fourteen days later, ALL was diagnosed via BMA.
In the control group, only the two patients suffering from SAA or sickle cell anemia received antibiotics ahead of BMA: Piperacillin/Tazobactam in the SAA patient two days ahead of BMA, and long-term Penicillin prophylaxis in the sickle cell patient.
3.2 Results of WMS
In total, an mean of 4.8 ng/µl of cell free DNA (cfDNA) could be isolated in the bone marrow of the ALL patients, compared to 2 ng/µl in the control group (Table 1, p=0.18). In peripheral blood, mean cfDNA concentration was 2.1 ng/µl in ALL and 0.3 ng/µl in the control group (p=0.08).
The cfDNA samples were processed in 6 different batches with different reagent lots due to logistical requirements. The batches separate into different clusters (Figure 1A) in contrast to the plot showing the different groups (patients, controls, synthetic sequencing run controls, Figure 1B) where no separation based on any of these groups was possible. In the read mapping statistics regarding the proportion of human reads versus non-human reads (Table 2), no significant differences between the groups were discernible (Figure 2A), with the exception of the synthetic control that shows a statistically significant increase in the non-human proportion, due to the synthetic microbial sequences (Figure 2B). In blood and bone marrow samples of patients or controls, no significant amounts of viral DNA were identified. DNA fragments for Torque Teno virus were found in two patients in the raw data, but below LOD and reporting thresholds.
Figure 1 PCA plot visualizing the clustering of datasets based on the six investigated batches (A) compared to the distribution in the PCA plot based on the different groups (patients, controls and synthetic run controls) (B).
Table 2 Read mapping statistics of whole metagenome shotgun sequencing for all patient samples (n=23) in comparison to synthetic control samples.
Figure 2 Comparison of mapped read proportions between the cohort of patients with newly diagnosed acute lymphatic leukemia (patient), the cohort of healthy donors and patients with non-oncologic diseases (control), and synthetic quality controls for the wetlab and sequencing process (seq-control) for human read fraction (A) and non-human read fraction (B). **: p≤0.01, ****: p≤0.0001, ns: non significant.
4 Discussion
In this pilot trial, WMS of bone marrow samples in pediatric patients with treatment-naïve ALL revealed no specific infectious agents and no differences in the viral or bacterial metagenome when compared to peripheral blood of the same patients (internal validation) or in comparison to peripheral blood and BM of healthy children or children with non-oncologic diseases (external validation). As a limitation, only DNA viruses could have been identified, which could explain the opposing results of a higher viral infection load in ALL previously reported (Francis et al., 2014) where RNA has been amplified.
The absence of a discernible bacterial metagenome in the peripheral blood and bone marrow samples after implementation of rigid controls, could also be linked to the recent discussion about the hypothesized “human blood microbiome”, which has been postulated by several small trials using NGS techniques on peripheral blood samples (McLaughlin et al., 2002; Moriyama et al., 2008; Païssé et al., 2016), reviewed by Castillo et al. (2019) A frequent issue that arises from the combination of high-throughput sequencing and low-biomass samples (Castillo et al., 2019) is a contamination of the used DNA extraction kits with commensal bacteria (“kitome”) as well as cross-contamination between different samples in the laboratory (“splashome”) that might emulate a “microbiome”. This has been demonstrated by Olomu et al. (2020), who opposed the findings of a hypothetical placental microbiome (Bassols et al., 2016; Zheng et al., 2017) by adjusting for these confounders. Just recently, the idea of a “human blood microbiome” has been rendered unlikely by the findings of a population-based study including nearly 10,000 healthy subjects (Tan et al., 2023). However, the small cohort size of our pilot study as well as the fact that several patients received antibiotic treatment ahead of inclusion into this trial have to be acknowledged as important limitations in interpreting our data.
As a conclusion, further unbiased investigation of a potential, infectious-related etiology of pediatric ALL using untargeted DNA amplification methods should always implement rigid controls against potential confounders, and might consider the detection of viral DNA and RNA in larger sample sizes.
Data availability statement
The data presented in the study are deposited in the NCBI BioProject repository, accession number PRJNA1123915.
Ethics statement
The studies involving humans were approved by Ethical committee of Medical Faculty, University Tuebingen, approval number 561/2020BO1, gathered 11/2020. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin.
Author contributions
ATH: Writing – review & editing, Writing – original draft, Validation, Investigation, Formal analysis, Data curation, Conceptualization. SG: Writing – review & editing, Writing – original draft, Validation, Software, Methodology, Investigation, Funding acquisition, Formal analysis, Data curation, Conceptualization. CS: Writing – review & editing, Writing – original draft, Validation, Methodology, Conceptualization. MD: Writing – review & editing, Data curation. MQ: Writing – review & editing, Data curation. RH: Writing – review & editing, Supervision, Project administration, Conceptualization. ME: Writing – review & editing, Supervision, Project administration, Data curation, Conceptualization.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. Sequencing of blood and bone marrow samples was performed and funded by the Noscendo GmbH, Königstr. 34, 47198 Duisburg, Germany. We kindly acknowledge the support by the “Stiftung für krebskranke Kinder Tübingen e.V.” and the Open Access Publishing Fund of University of Tübingen.
Conflict of interest
Author SG was employed by the company Noscendo GmbH.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2024.1355787/full#supplementary-material
References
Bassols, J., Serino, M., Carreras-Badosa, G., Burcelin, R., Blasco-Baque, V., Lopez-Bermejo, A., et al. (2016). Gestational diabetes is associated with changes in placental microbiota and microbiome. Pediatr. Res. 80, 777–784. doi: 10.1038/pr.2016.155
Castillo, D. J., Rifkin, R. F., Cowan, D. A., Potgieter, M. (2019). The healthy human blood microbiome: fact or fiction? Front. Cell Infect. Microbiol. 9, 148–148. doi: 10.3389/fcimb.2019.00148
Francis, S. S., Wendt, J., Wallace, A. D., Riley, L. W., Liu, F., Metayer, C., et al. (2014). Bone marrow microbiome: metagenomic comparsion of childhood acute lymphoblastic and acute myeloid leukamias. Blood 124, 3771–3771. doi: 10.1182/blood.V124.21.3771.3771
Greaves, M., Alexander, F. (1993). An infectious etiology for common acute lymphoblastic leukemia in childhood? Leukemia 7, 349–360.
Grumaz, S., Stevens, P., Grumaz, C., Decker, S. O., Weigand, M. A., Hofer, S., et al. (2016). Next-generation sequencing diagnostics of bacteremia in septic patients. Genome Med. 8 (1), 73. doi: 10.1186/s13073-016-0326-8
Kinlen, L., Dickson, M., Stiller, C. (1995). Childhood leukaemia and non-Hodgkin's lymphoma near large rural construction sites, with a comparison with Sellafield nuclear site. BMJ 310, 763–768. doi: 10.1136/bmj.310.6982.763
McLaughlin, R. W., Vali, H., Lau, P. C. K., Palfree, R. G. E., Ciccio, A. D., Sirois, M., et al. (2002). Are there naturally occurring pleomorphic bacteria in the blood of healthy humans? J. Clin. Microbiol. 40, 4771–4775. doi: 10.1128/JCM.40.12.4771-4775.2002
Moriyama, K., Ando, C., Tashiro, K., Kuhara, S., Okamura, S., Nakano, S., et al. (2008). Polymerase chain reaction detection of bacterial 16S rRNA gene in human blood. Microbiol. Immunol. 52, 375–382. doi: 10.1111/j.1348-0421.2008.00048.x
O’Connor, S. M., Boneva, R. S. (2007). Infectious etiologies of childhood leukemia: plausibility and challenges to proof. Environ. Health Perspectives. 115, 146–150. doi: 10.1289/ehp.9024
Olomu, I. N., Pena-Cortes, L. C., Long, R. A., Vyas, A., Krichevskiy, O., Luellwitz, R., et al. (2020). Elimination of “kitome” and “splashome” contamination results in lack of detection of a unique placental microbiome. BMC Microbiol. 20, 1–19. doi: 10.1186/s12866-020-01839-y
Païssé, S., Valle, C., Servant, F., Courtney, M., Burcelin, R., Amar, J., et al. (2016). Comprehensive description of blood microbiome from healthy donors assessed by 16S targeted metagenomic sequencing. Transfusion 56, 1138–1147. doi: 10.1111/trf.13477
Pui, C.-H. (1995). Childhood leukemias. New Engl. J. Med. 332 (24), 1618–1630. doi: 10.1056/NEJM199506153322407
Tan, C. C. S., Ko, K. K. K., Chen, H., Liu, J., Loh, M., Chia, M., et al. (2023). No evidence for a common blood microbiome based on a population study of 9,770 healthy humans. Nat. Microbiol. 8 (5), 973–985. doi: 10.1038/s41564-023-01350-w
Keywords: high-throughput sequencing, ALL, shotgun sequencing, infectious etiology, microbiome
Citation: Heinz AT, Grumaz S, Slavetinsky C, Döring M, Queudeville M, Handgretinger R and Ebinger M (2024) No evidence on infectious DNA-based agents in pediatric acute lymphoblastic leukemia using whole metagenome shotgun sequencing. Front. Cell. Infect. Microbiol. 14:1355787. doi: 10.3389/fcimb.2024.1355787
Received: 14 December 2023; Accepted: 24 May 2024;
Published: 21 June 2024.
Edited by:
Spiros Vlahopoulos, University of Athens, GreeceReviewed by:
Lingyun Zou, Chongqing Emergency Medical Center, ChinaJie Liu, Qingdao University, China
Copyright © 2024 Heinz, Grumaz, Slavetinsky, Döring, Queudeville, Handgretinger and Ebinger. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Martin Ebinger, martin.ebinger@med.uni-tuebingen.de