Search

Reproducibility of computational research is often challenging despite established guidelines and best practices. Translating these guidelines into practical applications remains difficult. Here, we present ENCORE, an approach to enhance transparency and reproducibility by guiding researchers in how to structure and document a computational project. ENCORE builds on previous efforts in computational reproducibility and integrates all project components into a standardized file system structure. It utilizes pre-defined files as documentation templates, leverages GitHub for software versioning, and includes an HTML-based navigator. ENCORE is designed to be agnostic to the type of computational project, data, programming language, and ICT infrastructure, and does not rely on specific software tools. We also share our group's experience using ENCORE, highlighting that the most significant challenge to the routine adoption of approaches like ours is the lack of incentives to motivate researchers to dedicate sufficient time and effort to ensure reproducibility., (© 2024. The Author(s).)

The reconstruction of clonal families (CFs) in B-cell receptor (BCR) repertoire analysis is a crucial step to understand the adaptive immune system and how it responds to antigens. The BCR repertoire of an individual is formed throughout life and is diverse due to several factors such as gene recombination and somatic hypermutation. The use of Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) using next generation sequencing enabled the generation of full BCR repertoires that also include rare CFs. The reconstruction of CFs from AIRR-seq data is challenging and several approaches have been developed to solve this problem. Currently, most methods use the heavy chain (HC) only, as it is more variable than the light chain (LC). CF reconstruction options include the definition of appropriate sequence similarity measures, the use of shared mutations among sequences, and the possibility of reconstruction without preliminary clustering based on V- and J-gene annotation. In this study, we aimed to systematically evaluate different approaches for CF reconstruction and to determine their impact on various outcome measures such as the number of CFs derived, the size of the CFs, and the accuracy of the reconstruction. The methods were compared to each other and to a method that groups sequences based on identical junction sequences and another method that only determines subclones. We found that after accounting for data set variability, in particular sequencing depth and mutation load, the reconstruction approach has an impact on part of the outcome measures, including the number of CFs. Simulations indicate that unique junctions and subclones should not be used as substitutes for CF and that more complex methods do not outperform simpler methods. Also, we conclude that different approaches differ in their ability to correctly reconstruct CFs when not considering the LC and to identify shared CFs. The results showed the effect of different approaches on the reconstruction of CFs and highlighted the importance of choosing an appropriate method., (© 2024. The Author(s).)

Objectives: To characterize the T cell receptor (TCRβ) repertoire in peripheral blood and muscle tissues of treatment naïve patients with newly diagnosed idiopathic inflammatory myopathies (IIMs)., Methods: High throughput RNA sequencing of the TCRβ chain was performed in peripheral blood and muscle tissue in twenty newly-diagnosed treatment-naïve IIM patients (9 DM, 5 NM/OM, 5 IMNM and 1 ASyS) and healthy controls. Results thereof were correlated with markers of disease activity., Results: Muscle tissue of IIM patients shows more expansion of TCRβ clones and decreased diversity when compared to peripheral blood of IIM as well as healthy controls (both p=0.0001). Several expanded TCRβ clones in muscle are tissue restricted and cannot be retrieved in peripheral blood. These clones have significantly longer CDR3 regions when compared to clones (also) found in circulation (p=0.0002), while their CDR3 region is more hydrophobic (p<0.01). Network analysis shows that clonal TCRβ signatures are shared between patients. Increased clonal expansion in muscle tissue is significantly correlated with increased CK levels (p=0.03), while it tends to correlate with decreased muscle strength (p=0.08)., Conclusion: Network analysis of clones in muscle of IIM patients shows shared clusters of sequences across patients. Muscle-restricted CDR3 TCRβ clones show specific structural features in their T cell receptor. Our results indicate that clonal TCRβ expansion in muscle tissue might be associated with disease activity. Collectively, these findings support a role for specific clonal T cell responses in muscle tissue in the pathogenesis of the IIM subtypes studied., Competing Interests: The 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., (Copyright © 2024 Anang, Walter, Lim, Niewold, van der Weele, Aronica, Eftimov, Raaphorst, van Schaik, van Kampen, van der Kooi and de Vries.)

Objective: To unravel B-cell receptor (BcR) characteristics in muscle tissues and peripheral blood and gain more insight into BcR repertoire changes in peripheral blood in idiopathic inflammatory myopathies (IIMs), and study how this correlates to the clinical response to IVIG., Methods: Nineteen treatment-naive patients with newly diagnosed IIM were prospectively treated with IVIG monotherapy. RNA-based BcR repertoire sequencing was performed in muscle biopsies collected before, and in peripheral blood (PB) collected before and nine weeks after IVIG treatment. Results were correlated to patients' clinical improvement based on the total improvement score (TIS)., Results: Prior to IVIG treatment, BcR clones found in muscle tissue could be retrieved in peripheral blood. Nine weeks after IVIG treatment, new patient-specific dominant BcR clones appeared in peripheral blood while pre-treatment dominant BcR clones disappeared. The cumulative frequency of all dominant BcR clones before treatment was significantly higher in individuals who responded to IVIG compared with those who did not respond to IVIG, and correlated with a higher CK. During follow-up, a decrease in the cumulative frequency of all dominant clones correlated with a higher TIS., Conclusion: In treatment-naive patients with newly diagnosed IIM, muscle tissue and peripheral blood share expanded BcR clones. In our study a higher cumulative frequency of dominant BcR clones in blood before treatment was associated with a higher CK and better treatment response, suggesting that response to IVIG may depend on the composition of the pre-treatment BcR repertoire., (© The Author(s) 2022. Published by Oxford University Press on behalf of the British Society for Rheumatology.)

A significant proportion of multiple sclerosis (MS) patients treated with interferon beta-1a (Rebif™) develop anti-drug antibodies (ADA) with a negative impact on treatment efficacy. We hypothesized that high-throughput B-cell receptor (BCR) repertoire analysis could be used to predict and monitor ADA development. To study this we analyzed 228 peripheral blood samples from 68 longitudinally followed patients starting on interferon beta-1a. Our results show that whole blood BCR analysis does not reflect, and does not predict ADA development in MS patients treated with interferon beta-1a. We propose that BCR analysis of phenotypically selected cell subsets or tissues might be more informative., (Copyright © 2022. Published by Elsevier B.V.)

Delineating the origins and properties of antibodies elicited by SARS-CoV-2 infection and vaccination is critical for understanding their benefits and potential shortcomings. Therefore, we investigate the SARS-CoV-2 spike (S)-reactive B cell repertoire in unexposed individuals by flow cytometry and single-cell sequencing. We show that ∼82% of SARS-CoV-2 S-reactive B cells harbor a naive phenotype, which represents an unusually high fraction of total human naive B cells (∼0.1%). Approximately 10% of these naive S-reactive B cells share an IGHV1-69/IGKV3-11 B cell receptor pairing, an enrichment of 18-fold compared to the complete naive repertoire. Following SARS-CoV-2 infection, we report an average 37-fold enrichment of IGHV1-69/IGKV3-11 B cell receptor pairing in the S-reactive memory B cells compared to the unselected memory repertoire. This class of B cells targets a previously undefined non-neutralizing epitope on the S2 subunit that becomes exposed on S proteins used in approved vaccines when they transition away from the native pre-fusion state because of instability. These findings can help guide the improvement of SARS-CoV-2 vaccines., (© 2022. The Author(s).)

A major goal of current HIV-1 vaccine design efforts is to induce broadly neutralizing antibodies (bNAbs). The VH1-2-derived bNAb IOMA directed to the CD4-binding site of the HIV-1 envelope glycoprotein is of interest because, unlike the better-known VH1-2-derived VRC01-class bNAbs, it does not require a rare short light chain complementarity-determining region 3 (CDRL3). Here, we describe three IOMA-class NAbs, ACS101-103, with up to 37% breadth, that share many characteristics with IOMA, including an average-length CDRL3. Cryo-electron microscopy revealed that ACS101 shares interactions with those observed with other VH1-2 and VH1-46-class bNAbs, but exhibits a unique binding mode to residues in loop D. Analysis of longitudinal sequences from the patient suggests that a transmitter/founder-virus lacking the N276 glycan might have initiated the development of these NAbs. Together these data strengthen the rationale for germline-targeting vaccination strategies to induce IOMA-class bNAbs and provide a wealth of sequence and structural information to support such strategies., (© 2022. The Author(s).)

Background: In patients with rheumatoid arthritis (RA) different joints were shown to share the same dominant T-cell clones, suggesting shared characteristics of the inflammatory process and indicating that strategies to selectively target the antigen receptor might be feasible. Since T- and B-lymphocytes closely interact in adaptive responses, we analysed to what extent different joints also share dominant B-cell clones., Methods: In 11 RA patients, quantitative B-cell receptor (BCR) repertoire analysis was performed in simultaneously obtained samples from inflamed synovial tissue (ST) from distinct locations within one joint, from multiple joints, from synovial fluid (SF) and peripheral blood (PB)., Results: ST biopsies from different locations in the same joint showed clear overlap in the top-25 dominant BCR clones (16.7%, SD 12.5), in the same range as the overlap between ST and SF in the same joint (8.0%, SD 8.8) and the overlap between ST-ST between different joints (9.1%, SD 8.2), but clearly higher than the overlap between ST and PB (1.7%, SD 2.4; p<0.05) and SF and PB (2.7%, SD 4.1; p<0.05). Interestingly, these figures were substantially lower than the overlap observed in previous T-cell clonality studies., Conclusions: We conclude that in RA BCR clonal responses may be more localized than TCR clonal responses, pointing to antigen-selective influx, proliferation and/or maturation of B-cells. B lineage cells in the SF may adequately represent the dominant BCR clones of the ST, which is in contrast to T-cells. Collectively, the presence of shared B- and especially T-cells in different joints from the same patient suggests that approaches might be feasible that aim to develop antigen-receptor specific targeting of lymphocyte clones in RA as an alternative to more generalized immunosuppressive strategies., Competing Interests: The 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., (Copyright © 2022 Musters, Balzaretti, van Schaik, Jongejan, van der Weele, Tas, van Kampen and de Vries.)

Objective: To comparatively analyse the aberrant affinity maturation of the antinuclear and rheumatoid factor (RF) B cell repertoires in blood and tissues of patients with Sjögren's syndrome (SjS) using an integrated omics workflow., Methods: Peptide sequencing of anti-Ro60, anti-Ro52, anti-La and RF was combined with B cell repertoire analysis at the DNA, RNA and single cell level in blood B cell subsets, affected salivary gland and extranodal marginal zone lymphomas of mucosa-associated lymphoid tissue (MALT) of patients with SjS., Results: Affected tissues contained anti-Ro60, anti-Ro52, anti-La and RF clones as a small part of a polyclonal infiltrate. Anti-Ro60, anti-La and anti-Ro52 clones outnumbered RF clones. MALT lymphoma tissues contained monoclonal RF expansions. Autoreactive clones were not selected from a restricted repertoire in a circulating B cell subset. The antinuclear antibody (ANA) repertoires displayed similar antigen-dependent and immunoglobulin (Ig) G1-directed affinity maturation. RF clones displayed antigen-dependent, IgM-directed and more B cell receptor integrity-dependent affinity maturation. This coincided with extensive intra-clonal diversification in RF-derived lymphomas. Regeneration of clinical disease manifestations after rituximab coincided with large RF clones, which not necessarily belonged to the lymphoma clone, that displayed continuous affinity maturation and intra-clonal diversification., Conclusion: The ANA and RF repertoires in patients with SjS display tissue-restricted, antigen-dependent and divergent affinity maturation. Affinity maturation of RF clones deviates further during RF clone derived lymphomagenesis and during regeneration of the autoreactive repertoire after temporary disruption by rituximab. These data give insight into the molecular mechanisms of autoreactive inflammation in SjS, assist MALT lymphoma diagnosis and allow tracking its response to rituximab., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2022. Re-use permitted under CC BY. Published by BMJ.)

High-throughput T-cell receptor repertoire sequencing constitutes a powerful tool to study T cell responses at the clonal level. However, it does not give information on the functional phenotype of the responding clones and lacks a statistical framework for quantitative evaluation. To overcome this, we combined datasets from different experiments, all starting from the same blood samples. We used a novel, sensitive, UMI-based protocol to perform repertoire analysis on experimental replicates. Applying established bioinformatic routines for transcriptomic expression analysis we explored the dynamics of antigen-induced clonal expansion after in vitro stimulation, identified antigen-responsive clones, and confirmed their activation status using the expression of activation markers upon antigen re-challenge. We demonstrate that the addition of IL-4 after antigen stimulation drives the expansion of T cell clones encoding unique receptor sequences. We show that our approach represents a scalable, high-throughput immunological tool, which can be used to identify and characterize antigen-responsive T cells at clonal level., Competing Interests: The 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., (Copyright © 2021 Pollastro, de Bourayne, Balzaretti, Jongejan, van Schaik, Niewold, van Kampen, Maillère and de Vries.)

The genetic background of childhood body mass index (BMI), and the extent to which the well-known associations of childhood BMI with adult diseases are explained by shared genetic factors, are largely unknown. We performed a genome-wide association study meta-analysis of BMI in 61,111 children aged between 2 and 10 years. Twenty-five independent loci reached genome-wide significance in the combined discovery and replication analyses. Two of these, located near NEDD4L and SLC45A3, have not previously been reported in relation to either childhood or adult BMI. Positive genetic correlations of childhood BMI with birth weight and adult BMI, waist-to-hip ratio, diastolic blood pressure and type 2 diabetes were detected (Rg ranging from 0.11 to 0.76, P-values <0.002). A negative genetic correlation of childhood BMI with age at menarche was observed. Our results suggest that the biological processes underlying childhood BMI largely, but not completely, overlap with those underlying adult BMI. The well-known observational associations of BMI in childhood with cardio-metabolic diseases in adulthood may reflect partial genetic overlap, but in light of previous evidence, it is also likely that they are explained through phenotypic continuity of BMI from childhood into adulthood., Competing Interests: I have read the journal's policy and the authors of this manuscript have the following competing interests: MMcC: The views expressed in this article are those of the author(s) and not necessarily those of the NHS, the NIHR, or the Department of Health. He serves on advisory panels for Pfizer, NovoNordisk, Zoe Global; has received honoraria from Merck, Pfizer, NovoNordisk and Eli Lilly; has stock options in Zoe Global; has received research funding from Abbvie, Astra Zeneca, Boehringer Ingelheim, Eli Lilly, Janssen, Merck, NovoNordisk, Pfizer, Roche, Sanofi Aventis, Servier & Takeda. MS receives funding from Pfizer Inc. for a project not related to this research. IB and spouse own stock in GlaxoSmithKline and Incyte Corp. ZE and CDB currently serve on the editorial board of PLOS Genetics. AC reports personal fees from Novartis, personal fees from Thermo Fisher Scientific, personal fees from Philips, personal fees from Sanofi, personal fees from Stallergenes Greer, outside the submitted work. KC in involved in consultancy for Danone Research, LNC-therapeutic and Confo-therapeutic but no personal funding is received and this activity not linked to the present research.

There is limited evidence on association between adherence to the Dietary Approaches to Stop Hypertension (DASH diet) and a lower blood pressure (BP) in children. In a population-based cohort study, among 1068 Dutch children aged 5 to 7, we evaluated the association between a DASH-type diet, 29 known genetic variants incorporated in a genetic risk score, and their interaction on BP. We calculated DASH score based on the food intake data measured through a validated 71-item food frequency questionnaire. In our sample, DASH score ranged from 9 (low adherence to the DASH diet) to 33 (median=21), and genetic score ranged from 18 (low genetic risk on high BP) to 41 (median=29). After adjustment for covariates, each 10 unit increase in DASH score was associated with a lower systolic BP of 0.7 mm Hg ( P =0.033). DASH score was negatively associated with hypertension (odds ratio=0.96 [0.92-0.99], P =0.044). Similarly, each SD increment in genetic score was associated with 0.5 mm Hg higher diastolic BP ( P =0.002). We found a positive interaction between low DASH score and high genetic score on diastolic BP adjusted for BP risk factors (β=1.52, P interaction =0.019 in additive scale and β=0.03, P interaction =0.021 in multiplicative scale). Our findings show that adherence to the DASH-type diet, as well as a low (adult-derived) genetic risk profile for BP, is associated with lower BP in children and that the genetic basis of BP phenotypes at least partly overlaps between adults and children. In addition, we found evidence of a gene-diet interaction on BP in children.

Objective: To gain more insight into the dynamics of lymphocyte depletion and develop new predictors of clinical response to rituximab in rheumatoid arthritis (RA)., Methods: RNA-based next-generation sequencing was used to analyse the B cell receptor (BCR) repertoire in peripheral blood and synovial tissue samples collected from 24 seropositive patients with RA treated with rituximab. Clonal expansion, mutation load and clonal overlap were assessed in samples collected before, at week 4 and at week 16 or 24 after treatment and correlated to the patients' clinical response., Results: After 4 weeks of rituximab-induced B cell depletion, the peripheral blood BCR repertoire of treated patients consisted of fewer, more dominant and more mutated BCR clones. No significant changes in the synovial tissue BCR repertoire were detected until week 16 post-treatment, when a reduced clonal overlap with baseline and an increased mutation load were observed. In patients who were non-responders at month 3 (n=5) using the European League Against Rheumatism response criteria, peripheral blood samples taken at week 4 after rituximab treatment showed more dominant clones compared with moderate responders (n=9) (median (IQR): 36 (27-52) vs 18 (16-26); p<0.01) and more clonal overlap with the baseline (median (IQR): 5% (2%-20%) vs 0% (0%-0%); p≤0.01)., Conclusion: Significant changes in BCR clonality are observed in peripheral blood of patients 4 weeks after rituximab treatment, while changes in synovial tissue were observed at later time points. Incomplete depletion of the dominant baseline peripheral blood BCR repertoire in the first month of treatment might predict clinical non-response at 3 months., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2019. Re-use permitted under CC BY. Published by BMJ.)

Although hundreds of genome-wide association studies-implicated loci have been reported for adult obesity-related traits, less is known about the genetics specific for early-onset obesity and with only a few studies conducted in non-European populations to date. Searching for additional genetic variants associated with childhood obesity, we performed a trans-ancestral meta-analysis of 30 studies consisting of up to 13 005 cases (≥95th percentile of body mass index (BMI) achieved 2-18 years old) and 15 599 controls (consistently <50th percentile of BMI) of European, African, North/South American and East Asian ancestry. Suggestive loci were taken forward for replication in a sample of 1888 cases and 4689 controls from seven cohorts of European and North/South American ancestry. In addition to observing 18 previously implicated BMI or obesity loci, for both early and late onset, we uncovered one completely novel locus in this trans-ancestral analysis (nearest gene, METTL15). The variant was nominally associated with only the European subgroup analysis but had a consistent direction of effect in other ethnicities. We then utilized trans-ancestral Bayesian analysis to narrow down the location of the probable causal variant at each genome-wide significant signal. Of all the fine-mapped loci, we were able to narrow down the causative variant at four known loci to fewer than 10 single nucleotide polymorphisms (SNPs) (FAIM2, GNPDA2, MC4R and SEC16B loci). In conclusion, an ethnically diverse setting has enabled us to both identify an additional pediatric obesity locus and further fine-map existing loci., (© The Author(s) 2019. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Many autoimmune diseases are hallmarked by autoreactive B and plasma cell responses that are directly or indirectly involved in disease pathogenesis. These B-cell responses show large variability between diseases, both in terms of the secreted autoantibody repertoire and the dynamics and characteristics of the underlying B-cell responses. Hence, different mechanisms have been proposed to explain the emergence of autoreactive B cells in an otherwise self-tolerant immune system. Notably, most mechanistic insights have been obtained from murine studies using models harboring genetic modifications of B and T cells. Given recent technological advances that have rendered autoreactive human B cells accessible for analysis, we here discuss the phenomenon of extensive N -glycosylation of the B-cell receptor (BCR) variable domain of a prototypic human autoreactive B-cell response and its potential role in the generation of autoimmunity. Anti-citrullinated protein antibodies (ACPA) hallmark the most disease-specific autoimmune response in Rheumatoid Arthritis (RA). Remarkably, ACPA-IgG are heavily N -glycosylated in the variable domain due to somatic mutations that generate abundant N -glycosylation consensus sequences. These sites, obtained from full-length BCR sequences of ACPA-expressing B cells from 12 ACPA-positive RA patients, were here analyzed in detail. Sites that required a single nucleotide mutation to be generated were defined as single somatic hypermutation (s-SHM) sites, whereas sites requiring multiple mutations were defined as m-SHM sites. IgG sequences of 12 healthy donors were used as control. Computational modeling of the germinal center reaction (CLONE algorithm) was used with the germline counterparts of ACPA-IgG heavy chain (HC) sequences to simulate the germinal center response. Our analyses revealed an abundance of N -glycosylation sites in ACPA-IgG HC that frequently required multiple mutations and predominated in specific positions. Based on these data, and taking into account recent insights into the dynamics of the ACPA-response during disease development, we here discuss the hypothesis that N -glycosylation sites in ACPA-IgG variable domains could lead to alternative, possibly antibody affinity-independent selection forces. Presumably, this occurs during germinal center responses allowing these B cells to escape from putative tolerance checkpoints, thereby driving autoreactive B cell development in the pathogenesis of RA., (Copyright © 2019 Vergroesen, Slot, van Schaik, Koning, Rispens, van Kampen, Toes and Scherer.)

Birth weight variation is influenced by fetal and maternal genetic and non-genetic factors, and has been reproducibly associated with future cardio-metabolic health outcomes. In expanded genome-wide association analyses of own birth weight (n = 321,223) and offspring birth weight (n = 230,069 mothers), we identified 190 independent association signals (129 of which are novel). We used structural equation modeling to decompose the contributions of direct fetal and indirect maternal genetic effects, then applied Mendelian randomization to illuminate causal pathways. For example, both indirect maternal and direct fetal genetic effects drive the observational relationship between lower birth weight and higher later blood pressure: maternal blood pressure-raising alleles reduce offspring birth weight, but only direct fetal effects of these alleles, once inherited, increase later offspring blood pressure. Using maternal birth weight-lowering genotypes to proxy for an adverse intrauterine environment provided no evidence that it causally raises offspring blood pressure, indicating that the inverse birth weight-blood pressure association is attributable to genetic effects, and not to intrauterine programming.

Genetic and immunological evidence clearly points to a role for T cells in the pathogenesis of rheumatoid arthritis (RA). Selective targeting of such disease-associated T cell clones might be highly effective while having few side effects. However, such selective targeting may only be feasible if the same T cell clones dominate the immune response at different sites of inflammation. We leveraged high-throughput technology to quantitatively assess whether different T cell clones dominate the inflammatory infiltrate at various sites of inflammation in this prototypic autoimmune disease. In 13 RA patients, we performed quantitative next-generation sequencing-based human TCRβ repertoire analysis in simultaneously obtained samples from inflamed synovial tissue (ST) from distinct locations within one joint, from multiple joints, and from synovial fluid (SF) and peripheral blood (PB). Identical TCRβ clones dominate inflammatory responses in ST samples taken from different locations within a single joint and when sampled in different joints. Although overall ST-SF overlap was comparable to higher ST-ST values, the overlap in dominant TCRβ clones in ST-SF comparisons was much lower than ST-ST and comparable to the low ST-PB overlap. In individual RA patients, a limited number of TCRβ clones dominate the immune response in the inflamed ST regardless of the location within a joint and which joint undergoes biopsy; in contrast, there is limited overlap of ST with SF or PB TCR repertoires. This limited breadth of the T cell response in ST of the individual RA patient indicates that development of immunotherapies that selectively modulate dominant T cell responses might be feasible., (Copyright © 2018 by The American Association of Immunologists, Inc.)

Virus discovery combining sequence unbiased amplification with next generation sequencing is now state-of-the-art. We have previously determined that the performance of the unbiased amplification technique which is operational at our institute, VIDISCA-454, is efficient when respiratory samples are used as input. The performance of the assay is, however, not known for other clinical materials like blood or stool samples. Here, we investigated the sensitivity of VIDISCA-454 with feces-suspensions and serum samples that are positive and that have been quantified for norovirus and human immunodeficiency virus type 1, respectively. The performance of VIDISCA-454 in serum samples was equal to its performance in respiratory material, with an estimated lower threshold of 1,000 viral genome copies. The estimated threshold in feces-suspension is around 200,000 viral genome copies. The decreased sensitivity in feces suspension is mainly due to sequences that share no recognizable identity with known sequences. Most likely these sequences originate from bacteria and phages which are not completely sequenced.