Benchmarking digital PCR partition classification methods with empirical and simulated duplex data.

Digital PCR (dPCR) is a highly accurate technique for the quantification of target nucleic acid(s). It has shown great potential in clinical applications, like tumor liquid biopsy and validation of biomarkers. Accurate classification of partitions based on end-point fluorescence intensities is crucial to avoid biased estimators of the concentration of the target molecules. We have evaluated many clustering methods, from general-purpose methods to specific methods for dPCR and flowcytometry, on both simulated and real-life data. Clustering method performance was evaluated by simulating various scenarios. Based on our extensive comparison of clustering methods, we describe the limits of these methods, and formulate guidelines for choosing an appropriate method. In addition, we have developed a novel method for simulating realistic dPCR data. The method is based on a mixture distribution of a Poisson point process and a skew-$t$ distribution, which enables the generation of irregularities of cluster shapes and randomness of partitions between clusters ('rain') as commonly observed in dPCR data. Users can fine-tune the model parameters and generate labeled datasets, using their own data as a template. Besides, the database of experimental dPCR data augmented with the labeled simulated data can serve as training and testing data for new clustering methods. The simulation method is available as an R Shiny app.

Digital PCR Partition Classification.

BACKGROUND: Partition classification is a critical step in the digital PCR data analysis pipeline. A range of partition classification methods have been developed, many motivated by specific experimental setups. An overview of these partition classification methods is lacking and their comparative properties are often unclear, likely impacting the proper application of these methods. CONTENT: This review provides a summary of all available digital PCR partition classification approaches and the challenges they aim to overcome, serving as a guide for the digital PCR practitioner wishing to apply them. We additionally discuss strengths and weaknesses of these methods, which can further guide practitioners in vigilant application of these existing methods. This review provides method developers with ideas for improving methods or designing new ones. The latter is further stimulated by our identification and discussion of application gaps in the literature, for which there are currently no or few methods available. SUMMARY: This review provides an overview of digital PCR partition classification methods, their properties, and potential applications. Ideas for further advances are presented and may bolster method development.

Triplex digital PCR assays for the quantification of intact proviral HIV-1 DNA.

The development of an HIV-1 cure is hampered by the existence of a persistent (latent) reservoir that contains a small proportion of replication-competent intact proviruses which refuels viral replication upon treatment discontinuation. Therefore, an accurate evaluation and quantification of these (intact) proviruses is essential to determine the efficacy of HIV-1 cure strategies which aim to eliminate this reservoir. Here, we present two triplex digital PCR assays which resulted from a combination of two existing methods, the IPDA (a 2-colour digital PCR based method) and Q4PCR assays (4 colour qPCR method), and tested the functionality on a three-colour digital PCR platform. In the present paper, we provide a step-by-step experimental protocol for these triplex digital PCR assays and validate their performance on a latently infected Jurkat cell-line model and HIV-1 patient samples. Our data demonstrates the potential and flexibility of increasing the number of subgenomic regions of HIV-1 within the IPDA to acquire sensitive detection of the HIV-1 reservoir while benefitting from the advantages of a dPCR setup.

Evaluating lncRNA Expression Patterns during HIV-1 Treatment Interruption.

Lately, the interest in long non-coding RNAs (lncRNAs) as potential drug targets and predictive markers in the context of HIV-1 has peaked, but their in vivo expression and regulation remains largely unexplored. Therefore, the present study examined lncRNA expression patterns during a clinical antiretroviral treatment interruption (ATI) trial. Peripheral blood mononuclear cells were isolated from ten patients at four timepoints: prior to ATI, 7-15 days after stop, at viral rebound and 3 months post antiretroviral therapy re-initiation. RNA was extracted and RT-qPCR on five known HIV-1-related lncRNAs (HEAL, MALAT1, NEAT1, GAS5 and NRON) was performed and correlated with HIV-1 and host marker expression. All lncRNAs correlated stronger with interferon stimulated genes (ISGs) than with HIV-1 reservoir and replication markers. However, one lncRNA, HEAL, showed significant upregulation at viral rebound during ATI compared to baseline and re-initiation of therapy (p = 0.0010 and p = 0.0094, respectively), following a similar viral-load-driven expression pattern to ISGs. In vitro knockdown of HEAL caused a significant reduction in HIV-1 infection levels, validating HEAL's importance for HIV-1 replication. We conclude that the HIV-1-promoting lncRNA HEAL is upregulated at viral rebound during ATI, most likely induced by viral cues.

On the whereabouts of SARS-CoV-2 in the human body: A systematic review.

Since SARS-CoV-2 appeared in the human population, the scientific community has scrambled to gather as much information as possible to find good strategies for the containment and treatment of this pandemic virus. Here, we performed a systematic review of the current (pre)published SARS-CoV-2 literature with a focus on the evidence concerning SARS-CoV-2 distribution in human tissues and viral shedding in body fluids. In addition, this evidence is aligned with published ACE2 entry-receptor (single cell) expression data across the human body to construct a viral distribution and ACE2 receptor body map. We highlight the broad organotropism of SARS-CoV-2, as many studies identified viral components (RNA, proteins) in multiple organs, including the pharynx, trachea, lungs, blood, heart, vessels, intestines, brain, male genitals and kidneys. This also implicates the presence of viral components in various body fluids such as mucus, saliva, urine, cerebrospinal fluid, semen and breast milk. The main SARS-CoV-2 entry receptor, ACE2, is expressed at different levels in multiple tissues throughout the human body, but its expression levels do not always correspond with SARS-CoV-2 detection, indicating that there is a complex interplay between virus and host. Together, these data shed new light on the current view of SARS-CoV-2 pathogenesis and lay the foundation for better diagnosis and treatment of COVID-19 patients.

Recurrent chromosomal imbalances provide selective advantage to human embryonic stem cells under enhanced replicative stress conditions.

Human embryonic stem cells (hESCs) and embryonal tumors share a number of common features, including a compromised G1/S checkpoint. Consequently, these rapidly dividing hESCs and cancer cells undergo elevated levels of replicative stress, inducing genomic instability that drives chromosomal imbalances. In this context, it is of interest that long-term in vitro cultured hESCs exhibit a remarkable high incidence of segmental DNA copy number gains, some of which are also highly recurrent in certain malignancies such as 17q gain (17q+). The selective advantage of DNA copy number changes in these cells has been attributed to several underlying processes including enhanced proliferation. We hypothesized that these recurrent chromosomal imbalances become rapidly embedded in the cultured hESCs through a replicative stress driven Darwinian selection process. To this end, we compared the effect of hydroxyurea-induced replicative stress vs normal growth conditions in an equally mixed cell population of isogenic euploid and 17q + hESCs. We could show that 17q + hESCs rapidly overtook normal hESCs. Our data suggest that recurrent chromosomal segmental gains provide a proliferative advantage to hESCs under increased replicative stress, a process that may also explain the highly recurrent nature of certain imbalances in cancer.

IRF7 and RNH1 are modifying factors of HIV-1 reservoirs: a genome-wide association analysis.

BACKGROUND: Combination antiretroviral treatment (cART) cannot eradicate HIV-1 from the body due to the establishment of persisting viral reservoirs which are not affected by therapy and reinitiate new rounds of HIV-1 replication after treatment interruption. These HIV-1 reservoirs mainly comprise long-lived resting memory CD4+ T cells and are established early after infection. There is a high variation in the size of these viral reservoirs among virally suppressed individuals. Identification of host factors that contribute to or can explain this observed variation could open avenues for new HIV-1 treatment strategies. METHODS: In this study, we conducted a genome-wide quantitative trait locus (QTL) analysis to probe functionally relevant genetic variants linked to levels of cell-associated (CA) HIV-1 DNA, CA HIV-1 RNA, and RNA:DNA ratio in CD4+ T cells isolated from blood from a cohort of 207 (Caucasian) people living with HIV-1 (PLHIV) on long-term suppressive antiretroviral treatment (median = 6.6 years). CA HIV-1 DNA and CA HIV-1 RNA levels were measured with corresponding droplet digital PCR (ddPCR) assays, and genotype information of 522,455 single-nucleotide variants was retrieved via the Infinium Global Screening array platform. RESULTS: The analysis resulted in one significant association with CA HIV-1 DNA (rs2613996, P < 5 × 10-8) and two suggestive associations with RNA:DNA ratio (rs7113204 and rs7817589, P < 5 × 10-7). Then, we prioritized PTDSS2, IRF7, RNH1, and DEAF1 as potential HIV-1 reservoir modifiers and validated that higher expressions of IRF7 and RNH1 were accompanied by rs7113204-G. Moreover, RNA:DNA ratio, indicating relative HIV-1 transcription activity, was lower in PLHIV carrying this variant. CONCLUSIONS: The presented data suggests that the amount of CA HIV-1 DNA and RNA:DNA ratio can be influenced through PTDSS2, RNH1, and IRF7 that were anchored by our genome-wide association analysis. Further, these observations reveal potential host genetic factors affecting the size and transcriptional activity of HIV-1 reservoirs and could indicate new targets for HIV-1 therapeutic strategies.

Digital Polymerase Chain Reaction for Assessment of Mutant Mitochondrial Carry-over after Nuclear Transfer for In Vitro Fertilization.

BACKGROUND: The quantification of mitochondrial DNA heteroplasmy for the diagnosis of mitochondrial disease or after mitochondrial donation, is performed mainly using next-generation sequencing strategies (NGS). Digital PCR (dPCR) has the potential to offer an accurate alternative for mutation load quantification. METHODS: We assessed the mutation load of 23 low-input human samples at the m.11778 locus, which is associated with Leber's hereditary optic neuropathy (LHON) using 2 droplet digital PCR platforms (Stilla Naica and Bio-Rad QX200) and the standard NGS strategy. Assay validation was performed by analyzing a titration series with mutation loads ranging from 50% to 0.01%. RESULTS: A good concordance in mutation rates was observed between both dPCR techniques and NGS. dPCR established a distinctly lower level of background noise compared to NGS. Minor alleles with mutation loads lower than 1% could still be detected, with standard deviations of the technical replicates varying between 0.07% and 0.44% mutation load. Although no significant systematic bias was observed when comparing dPCR and NGS, a minor proportional bias was detected. A slight overestimation of the minor allele was observed for the NGS data, most probably due to amplification and sequencing errors in the NGS workflow. CONCLUSION: dPCR has proven to be an accurate tool for the quantification of mitochondrial heteroplasmy, even for samples harboring a low mutation load (<1%). In addition, this alternative technique holds multiple benefits compared to NGS (e.g., less hands-on time, more straightforward data-analysis, and a lower up-front capital investment).

Evaluating predictive markers for viral rebound and safety assessment in blood and lumbar fluid during HIV-1 treatment interruption.

BACKGROUND: Validated biomarkers to evaluate HIV-1 cure strategies are currently lacking, therefore requiring analytical treatment interruption (ATI) in study participants. Little is known about the safety of ATI and its long-term impact on patient health. OBJECTIVES: ATI safety was assessed and potential biomarkers predicting viral rebound were evaluated. METHODS: PBMCs, plasma and CSF were collected from 11 HIV-1-positive individuals at four different timepoints during ATI (NCT02641756). Total and integrated HIV-1 DNA, cell-associated (CA) HIV-1 RNA transcripts and restriction factor (RF) expression were measured by PCR-based assays. Markers of neuroinflammation and neuronal injury [neurofilament light chain (NFL) and YKL-40 protein] were measured in CSF. Additionally, neopterin, tryptophan and kynurenine were measured, both in plasma and CSF, as markers of immune activation. RESULTS: Total HIV-1 DNA, integrated HIV-1 DNA and CA viral RNA transcripts did not differ pre- and post-ATI. Similarly, no significant NFL or YKL-40 increases in CSF were observed between baseline and viral rebound. Furthermore, markers of immune activation did not increase during ATI. Interestingly, the RFs SLFN11 and APOBEC3G increased after ATI before viral rebound. Similarly, Tat-Rev transcripts were increased preceding viral rebound after interruption. CONCLUSIONS: ATI did not increase viral reservoir size and it did not reveal signs of increased neuronal injury or inflammation, suggesting that these well-monitored ATIs are safe. Elevation of Tat-Rev transcription and induced expression of the RFs SLFN11 and APOBEC3G after ATI, prior to viral rebound, indicates that these factors could be used as potential biomarkers predicting viral rebound.

Interferon-Mediated Long Non-Coding RNA Response in Macrophages in the Context of HIV.

Interferons play a critical role in the innate immune response against a variety of pathogens, such as HIV-1. Recent studies have shown that long non-coding genes are part of a reciprocal feedforward/feedback relationship with interferon expression. They presumably contribute to the cell type specificity of the interferon response, such as the phenotypic and functional transition of macrophages throughout the immune response. However, no comprehensive understanding exists today about the IFN-lncRNA interplay in macrophages, also a sanctuary for latent HIV-1. Therefore, we completed a poly-A+ RNAseq analysis on monocyte-derived macrophages (MDMs) treated with members of all three types of IFNs (IFN-α, IFN-ε, IFN-γ or IFN-λ) and on macrophages infected with HIV-1, revealing an extensive non-coding IFN and/or HIV-1 response. Moreover, co-expression correlation with mRNAs was used to identify important (long) non-coding hub genes within IFN- or HIV-1-associated gene clusters. This study identified and prioritized IFN related hub lncRNAs for further functional validation.

Digital PCR as a tool to measure HIV persistence.

Although antiretroviral therapy is able to suppress HIV replication in infected patients, the virus persists and rebounds when treatment is stopped. In order to find a cure that can eradicate the latent reservoir, one must be able to quantify the persisting virus. Traditionally, HIV persistence studies have used real-time PCR (qPCR) to measure the viral reservoir represented by HIV DNA and RNA. Most recently, digital PCR is gaining popularity as a novel approach to nucleic acid quantification as it allows for absolute target quantification. Various commercial digital PCR platforms are nowadays available that implement the principle of digital PCR, of which Bio-Rad's QX200 ddPCR is currently the most used platform in HIV research. Quantification of HIV by digital PCR is proving to be a valuable improvement over qPCR as it is argued to have a higher robustness to mismatches between the primers-probe set and heterogeneous HIV, and forfeits the need for a standard curve, both of which are known to complicate reliable quantification. However, currently available digital PCR platforms occasionally struggle with unexplained false-positive partitions, and reliable segregation between positive and negative droplets remains disputed. Future developments and advancements of the digital PCR technology are promising to aid in the accurate quantification and characterization of the persistent HIV reservoir.

Short-Course Toll-Like Receptor 9 Agonist Treatment Impacts Innate Immunity and Plasma Viremia in Individuals With Human Immunodeficiency Virus Infection.

BACKGROUND.: Treatment with latency reversing agents (LRAs) enhances human immunodeficiency virus type 1 (HIV-1) transcription in vivo but leads to only modest reductions in the size of the reservoir, possibly due to insufficient immune-mediated elimination of infected cells. We hypothesized that a single drug molecule-a novel Toll-like receptor 9 (TLR9) agonist, MGN1703-could function as an enhancer of innate immunity and an LRA in vivo. METHODS.: We conducted a single-arm, open-label study in which 15 virologically suppressed HIV-1-infected individuals on antiretroviral therapy received 60 mg MGN1703 subcutaneously twice weekly for 4 weeks. We characterized plasmacytoid dendritic cell, natural killer (NK), and T-cell activation using flow cytometry on baseline and after 4 weeks of treatment. HIV-1 transcription was quantified by measuring plasma HIV-1 RNA during MGN1703 administration. RESULTS.: In accordance with the cell type-specific expression of TLR9, MGN1703 treatment led to pronounced activation of plasmacytoid dendritic cells and substantial increases in plasma interferon-α2 levels (P < .0001). Consistently, transcription of interferon-stimulated genes (eg, OAS1, ISG15, Mx1; each P < .0001) were upregulated in CD4+ T cells as demonstrated by RNA sequencing. Further, proportions of activated cytotoxic NK cells and CD8+ T cells increased significantly during MGN1703 dosing, suggesting an enhancement of cellular immune responses. In 6 of 15 participants, plasma HIV-1 RNA increased from <20 copies/mL to >1500 copies/mL (range, 21-1571 copies/mL) during treatment. CONCLUSIONS.: TLR9 agonist treatment in HIV infection has a dual potential by increasing HIV-1 transcription and enhancing cytotoxic NK cell activation, both of which are key outcomes in HIV-1 eradication therapy. CLINICAL TRIALS REGISTRATION.: NCT02443935.

Accurate quantification of episomal HIV-1 two-long terminal repeat circles by use of optimized DNA isolation and droplet digital PCR.

Episomal HIV-1 two-long terminal repeat (2-LTR) circles are considered markers for ongoing viral replication. Two sample processing procedures were compared to accurately quantify 2-LTR in patients by using droplet digital PCR (ddPCR). Here, we show that plasmid isolation with a spiked non-HIV plasmid for normalization enables more accurate 2-LTR quantification than genomic DNA isolation.

HIV-1 RNA and HIV-1 DNA persistence during suppressive ART with PI-based or nevirapine-based regimens.

OBJECTIVES: Whether ART regimens differ in their propensity to allow persistent HIV-1 detection remains unclear. To investigate this, we performed a cross-sectional study to characterize HIV-1 persistence in peripheral blood during suppressive therapy with NRTIs plus a PI or nevirapine. METHODS: Residual plasma HIV-1 RNA was quantified by real-time PCR. Cell-associated proviral total HIV-1 DNA, unspliced and multiply spliced HIV-1 RNA and 2-long terminal repeat (2-LTR) circles were quantified by digital PCR. RESULTS: Comparing PI with nevirapine recipients, residual plasma HIV-1 RNA detection rates were 47/80 (58.8%) versus 37/81 (45.7%), with median (IQR) levels of 4 (3-8) versus 4 (3-7) copies/mL (P = 0.207); detection was less likely with longer duration of suppressive ART (P = 0.020), independently of treatment. HIV-1 DNA was detected in all patients, with median levels of 2.3 (IQR 2.0-2.7) versus 2.5 (IQR 2.1-2.7) log10 copies/10(6) PBMCs, respectively; HIV-1 DNA levels were associated with pre-ART viral load (P = 0.004) and with residual HIV-1 RNA (P = 0.034), unspliced HIV-1 RNA (P = 0.001) and 2-LTR circles (P = 0.005), independently of treatment. CONCLUSIONS: No significant differences were revealed in levels of residual plasma HIV-1 RNA, total HIV-1 DNA or intracellular markers of ongoing virus replication (unspliced and multiply spliced HIV-1 RNA and 2-LTR circles) between treatment groups.

Robust regression methods for real-time polymerase chain reaction.

Current real-time polymerase chain reaction (PCR) data analysis methods implement linear least squares regression methods for primer efficiency estimation based on standard curve dilution series. This method is sensitive to outliers that distort the outcome and are often ignored or removed by the end user. Here, robust regression methods are shown to provide a reliable alternative because they are less affected by outliers and often result in more precise primer efficiency estimators than the linear least squares method.

ddpcRquant: threshold determination for single channel droplet digital PCR experiments.

Digital PCR is rapidly gaining interest in the field of molecular biology for absolute quantification of nucleic acids. However, the first generation of platforms still needs careful validation and requires a specific methodology for data analysis to distinguish negative from positive signals by defining a threshold value. The currently described methods to assess droplet digital PCR (ddPCR) are based on an underlying assumption that the fluorescent signal of droplets is normally distributed. We show that this normality assumption does not likely hold true for most ddPCR runs, resulting in an erroneous threshold. We suggest a methodology that does not make any assumptions about the distribution of the fluorescence readouts. A threshold is estimated by modelling the extreme values in the negative droplet population using extreme value theory. Furthermore, the method takes shifts in baseline fluorescence between samples into account. An R implementation of our method is available, allowing automated threshold determination for absolute ddPCR quantification using a single fluorescent reporter.

Replication competent virus as an important source of bias in HIV latency models utilizing single round viral constructs.

The central memory T cell (TCM) model forms a unique HIV-1 latency model based on primary cells that closely resemble in vivo TCM. The virus employed in this model is based on an engineered vector incapable of replication after initial infection. We show that despite this strategy, replication competent viral particles are released into the culture medium due to recombination between overlapping sequences of the env deleted HIV genome that is co-transfected with intact env. This finding emphasizes the need for careful data analysis and interpretation if similar constructs are employed and urges for additional caution during laboratory work.

Subpar reporting of pre-analytical variables in RNA-focused blood plasma studies.

Extracellular RNA (cell-free RNA; exRNA) from blood-derived liquid biopsies is an appealing, minimally invasive source of disease biomarkers. As pre-analytical variables strongly influence exRNA measurements, their reporting is essential for meaningful interpretation and replication of results. The aim of this review was to chart to what extent pre-analytical variables are documented, to pinpoint shortcomings and to improve future reporting. In total, 200 blood plasma exRNA studies published in 2018 or 2023 were reviewed for annotation of 22 variables associated with blood collection, plasma preparation, and RNA purification. Our results show that pre-analytical variables are poorly documented, with only three out of 22 variables described in over half of the publications. The percentage of variables reported ranged from 4.6% to 54.6% (mean 24.84%) in 2023 and from 4.6% to 57.1% (mean 28.60%) in 2018. Recommendations and guidelines (i.e., BRISQ, ASCO-CAP, BloodPAC, PPMPT, and CEN standards) have currently not resulted in improved reporting. In conclusion, our results highlight the lack of reporting pre-analytical variables in exRNA studies and advocate for a consistent use of available standards, endorsed by funders and journals.

Measuring DNA quality by digital PCR using probability calculations.

BACKGROUND: Accurate methods to assess DNA integrity are needed for many biomolecular methods. A multiplex digital PCR (dPCR) method designed for interspaced target sequences can be used to assess sequence integrity of large DNA strands. The ratio of single positive partitions versus double positive partitions is then used to calculate the sheared DNA strands. However, this simple calculation is only valid with low DNA concentration. We here describe a method based on probability calculations which enables DNA quality analysis in a large dynamic range of DNA concentrations. RESULTS: Known DNA integrity percentages were mimicked using artificial double stranded DNA in low, intermediate and high DNA concentration scenarios, respectively 600, 12500 and 30000 copies of DNA per reaction. At low concentrations both methods were similar. However, at the intermediate concentration (12500 copies per reaction) the ratio based method started producing a larger error than the proposed probability calculation method with a mean relative error of 20.7 and 16.7 for the Bruner and the proposed method respectively. At the high concentration (30000 copies per reaction) only the proposed method provided accurate measurements with a mean relative error of 60.9 and 9.3 for the ratio based and the proposed method respectively. Furthermore, while both methods have a bias, it is constant for the proposed method, while it decreases with the integrity of the DNA for the ratio based method. The probability calculation equation was extended to 4 dimensions and a proof of concept experiment was performed, the data suggested that the 4 dimensional equation is valid. SIGNIFICANCE AND NOVELTY: We here validate a method of estimating DNA integrity with dPCR using multiple probe combinations, allowing fast and flexible DNA integrity analysis. Additionally, we extend the method from 2 to 4 plex for more accurate DNA integrity measurements.