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.

Digital PCR: a tool in clostridial mutant selection and detection.

The ClosTron mutagenesis system has enabled researchers to efficiently edit the clostridial genome. Since site-specific insertion of the mobile ClosTron insert may cause errors, validation is key. In this paper we describe the use of digital PCR (dPCR) as an alternative tool in selecting clostridial mutant strains. Clostridium perfringens chitinase mutant strains were constructed in which the mobile ClosTron intron was inserted into one of the chitinase genes. On-target insertion of the mobile intron was validated through conventional PCR. In order to confirm the absence of off-target insertions, dPCR was used to determine the amount of the ClosTron intron as well as the amount of a reference gene, located in close proximity to the interrupted gene. Subsequently, mutant strains containing an equivalent amount of both genes were selected as these do not contain additional off-target mobile ClosTron inserts. The outcome of this selection procedure was confirmed through a validated PCR-based approach. In addition to its application in mutant selection, dPCR can be used in other aspects of clostridial research, such as the distinction and easy quantification of different types of strains (wildtype vs. mutant) in complex matrices, such as faecal samples, a process in which other techniques are hampered by bacterial overgrowth (plating) or inhibition by matrix contaminants (qPCR). This research demonstrates that dPCR is indeed a high-throughput method in the selection of clostridial insertion mutants as well as a robust and accurate tool in distinguishing between wildtype and mutant C. perfringens strains, even in a complex matrix such as faeces. KEY POINTS: • Digital PCR as an alternative in ClosTron mutant selection • Digital PCR is an accurate tool in bacterial quantification in a complex matrix • Digital PCR is an alternative tool with great potential to microbiological research.

Techniques used to assess intussusceptive angiogenesis: A systematic review.

Intussusceptive angiogenesis (IA) is an important physiological form of angiogenesis in which an existing vessel splits in two by the formation of an intraluminal tissue pillar. The presence of these intraluminal pillars form the hallmark of ongoing IA in growing vascular beds. However, their visualization is technically challenging. The goal of this systematic review was to investigate which techniques are being used to identify intraluminal pillars and to formulate important points to keep in mind when studying IA. A systematic literature search resulted in 154 evaluated articles of which the majority (65%) provided sufficient data to unambiguously demonstrate the presence of intraluminal pillars. Scanning electron microscopy imaging of vascular corrosion casts and serial sectioning of ultrathin sections are the most used techniques. New methods such as serial block face scanning electron microscopy and micro computed tomography (µCT) are gaining importance. Moreover, our results indicate that IA was studied in a variety of animals and tissues. IA is a biologically very relevant form of angiogenesis. Techniques to visualize intraluminal pillars need to have a minimal resolution of 1 µm and should provide information on the 3D-nature of the pillars. Optimally, several techniques are combined to demonstrate ongoing IA.

Paternally expressed retrotransposon Gag-like 1 gene, RTL1, is one of the crucial elements for placental angiogenesis in horses†.

RTL1 (retrotransposon Gag-like 1) is an essential gene in the development of the human and murine placenta. Several fetal and placental abnormalities such as intrauterine growth restriction (IUGR) and hydrops conditions have been associated with altered expression of this gene. However, the function of RTL1 has not been identified. RTL1 is located on a highly conserved region in eutherian mammals. Therefore, the genetic and molecular analysis in horses could hold important implications for other species, including humans. Here, we demonstrated that RTL1 is paternally expressed and is localized within the endothelial cells of the equine (Equus caballus) chorioallantois. We developed an equine placental microvasculature primary cell culture and demonstrated that RTL1 knockdown leads to loss of the sprouting ability of these endothelial cells. We further demonstrated an association between abnormal expression of RTL1 and development of hydrallantois. Our data suggest that RTL1 may be essential for placental angiogenesis, and its abnormal expression can lead to placental insufficiency. This placental insufficiency could be the reason for IUGR and hydrops conditions reported in other species, including humans.

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).

Development of a 3':5' digital PCR assay to determine horse mRNA integrity.

Accurate tools to measure RNA integrity are essential to obtain reliable gene expression data. The reverse transcription quantitative PCR (RT-qPCR) based 3':5' assay permits a direct determination of messenger RNA (mRNA) integrity. However, the use of standard curves and the possible effect of PCR inhibitors make this method cumbersome and prone to variation, especially in small samples. Here we developed a triplex digital PCR (dPCR) 3':5' assay for assessing RNA integrity in equine samples as rapid and simple alternative to RT-qPCR. This dPCR assay not only provides a straight forward analysis of the mRNA integrity, but also of its quantity.

Evaluation of HIV-1 reservoir levels as possible markers for virological failure during boosted darunavir monotherapy.

BACKGROUND: The gold standard for HIV-1 treatment is to administer triple antiretroviral therapy, but a shift to simplified regimens is being explored. Boosted darunavir monotherapy can be considered for patients who are for specific reasons not good candidates for dual or triple therapy. Still, a number of patients fail virologically or need to switch treatment. OBJECTIVES: To identify predictive markers for those patients that are more likely to sustain virological control under monotherapy, virological and immunological markers were explored in HIV-1-positive patients that experienced virological failure on ritonavir-boosted darunavir monotherapy in the PROTEA trial. METHODS: As a retrospective nested study of the PROTEA study (NCT01448707), we analysed 77 HIV-1-infected patients who were on darunavir/ritonavir 800/100 mg monotherapy up to 96 weeks. Patients were appointed to three distinct cohorts based on viral loads (VLs): (i) undetectable VL after 96 weeks; (ii) very-low-level viraemia (5-39 copies/mL); and (iii) failing treatment. Total HIV-1 DNA, integrated HIV-1 DNA and 2-long terminal repeat circular HIV-1 DNA (2LTR circles) were measured in PBMCs at baseline, week 48 and week 96. RESULTS: Total HIV-1 DNA and integrated HIV-1 DNA at baseline differed significantly between patients who experienced virological failure on monotherapy (P < 0.01 and P < 0.001). Although a higher level of HIV-1 DNA was measured in failures, this marker by itself does not provide enough predictive value to prospectively predict virological failure in patients on monotherapy. CONCLUSIONS: HIV-1 reservoir markers correlate with therapy failure in ritonavir-boosted darunavir monotherapy. However, their role as a predictive marker combined with other markers in a routine clinical setting should be further explored.

Integrated and Total HIV-1 DNA Predict Ex Vivo Viral Outgrowth.

The persistence of a reservoir of latently infected CD4 T cells remains one of the major obstacles to cure HIV. Numerous strategies are being explored to eliminate this reservoir. To translate these efforts into clinical trials, there is a strong need for validated biomarkers that can monitor the reservoir over time in vivo. A comprehensive study was designed to evaluate and compare potential HIV-1 reservoir biomarkers. A cohort of 25 patients, treated with suppressive antiretroviral therapy was sampled at three time points, with median of 2.5 years (IQR: 2.4-2.6) between time point 1 and 2; and median of 31 days (IQR: 28-36) between time point 2 and 3. Patients were median of 6 years (IQR: 3-12) on ART, and plasma viral load (<50 copies/ml) was suppressed for median of 4 years (IQR: 2-8). Total HIV-1 DNA, unspliced (us) and multiply spliced HIV-1 RNA, and 2LTR circles were quantified by digital PCR in peripheral blood, at 3 time points. At the second time point, a viral outgrowth assay (VOA) was performed, and integrated HIV-1 DNA and relative mRNA expression levels of HIV-1 restriction factors were quantified. No significant change was found for long- and short-term dynamics of all HIV-1 markers tested in peripheral blood. Integrated HIV-1 DNA was associated with total HIV-1 DNA (p<0.001, R² = 0.85), us HIV-1 RNA (p = 0.029, R² = 0.40), and VOA (p = 0.041, R2 = 0.44). Replication-competent virus was detected in 80% of patients by the VOA and it correlated with total HIV-1 DNA (p = 0.039, R² = 0.54). The mean quantification difference between Alu-PCR and VOA was 2.88 log10, and 2.23 log10 between total HIV-1 DNA and VOA. The levels of usHIV-1 RNA were inversely correlated with mRNA levels of several HIV-1 restriction factors (TRIM5α, SAMHD1, MX2, SLFN11, pSIP1). Our study reveals important correlations between the viral outgrowth and total and integrated HIV-1 DNA measures, suggesting that the total pool of HIV-1 DNA may predict the size of the replication-competent virus in ART suppressed patients.

Intussusceptive Pillar Formation in Developing Porcine Glomeruli.

BACKGROUND/AIMS: Intussusceptive angiogenesis (IA) is a dynamic process which contributes to vascular expansion and remodeling. Intraluminal pillars have long been the distinctive structural indicator of IA. However, the mechanism of their formation has not been fully elucidated. METHODS: Using light and electron microscopy, we studied intussusceptive vascular growth in the developing porcine metanephric kidney. RESULTS: We observed intraluminal pillars formed by endothelial cells in the vasculature of developing glomeruli. Their diameter was < 2.5 µm, consistent with the diameter of nascent pillars. TEM revealed that the majority of these pillars consisted only of endothelium. However, a central core of extracellular matrix (ECM) covered by endothelium, reminiscent of a more mature intussusceptive pillar, was also found in the lumen of a glomerular capillary. Perivascular cells or pericytes were not involved in the pillar structure during these stages of formation. CONCLUSION: This study shows ECM presence in a mature intussusceptive pillar without any perivascular cell involvement in the structure. This leads to the hypothesis that ECM deposition precedes the participation of these cells in the formation of intraluminal pillars during IA in porcine metanephric glomerular capillaries.

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.

HIV integration sites in latently infected cell lines: evidence of ongoing replication.

BACKGROUND: Assessing the location and frequency of HIV integration sites in latently infected cells can potentially inform our understanding of how HIV persists during combination antiretroviral therapy. We developed a novel high throughput sequencing method to evaluate HIV integration sites in latently infected cell lines to determine whether there was virus replication or clonal expansion in these cell lines observed as multiple integration events at the same position. RESULTS: We modified a previously reported method using random DNA shearing and PCR to allow for high throughput robotic processing to identify the site and frequency of HIV integration in latently infected cell lines. Latently infected cell lines infected with intact virus demonstrated multiple distinct HIV integration sites (28 different sites in U1, 110 in ACH-2 and 117 in J1.1 per 150,000 cells). In contrast, cell lines infected with replication-incompetent viruses (J-Lat cells) demonstrated single integration sites. Following in vitro passaging of the ACH-2 cell line, we observed a significant increase in the frequency of unique HIV integration sites and there were multiple mutations and large deletions in the proviral DNA. When the ACH-2 cell line was cultured with the integrase inhibitor raltegravir, there was a significant decrease in the number of unique HIV integration sites and a transient increase in the frequency of 2-LTR circles consistent with virus replication in these cells. CONCLUSION: Cell lines latently infected with intact HIV demonstrated multiple unique HIV integration sites indicating that these cell lines are not clonal and in the ACH-2 cell line there was evidence of low level virus replication. These findings have implications for the use of latently infected cell lines as models of HIV latency and for the use of these cells as standards.

The need for transparency and good practices in the qPCR literature.

Two surveys of over 1,700 publications whose authors use quantitative real-time PCR (qPCR) reveal a lack of transparent and comprehensive reporting of essential technical information. Reporting standards are significantly improved in publications that cite the Minimum Information for Publication of Quantitative Real-Time PCR Experiments (MIQE) guidelines, although such publications are still vastly outnumbered by those that do not.

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.

Evaluation of normalization strategies used in real-time quantitative PCR experiments in HepaRG cell line studies.

BACKGROUND: The HepaRG cell line is widely used as an alternative for primary human hepatocytes for numerous applications, including drug screening, and is progressively gaining importance as a human-relevant cell source. Consequently, increasing numbers of experiments are being performed with this cell line, including real-time quantitative PCR (RT-qPCR) experiments for gene expression studies. CONTENT: When RT-qPCR experiments are performed, results are reliable only when attention is paid to several critical aspects, including a proper normalization strategy. Therefore, in 2011 we determined the most optimal reference genes for gene expression studies in the HepaRG cell system, according to the MIQE (Minimum Information for Publication of Quantitative Real-Time PCR Experiments) guidelines. This study additionally provided clear evidence that the use of a single reference gene [glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal protein S18 (RPS18), or actin, beta (ACTB)] was insufficient for normalization in HepaRG cells. Our screening of relevant studies published after our study suggested that the findings of our study were completely ignored. SUMMARY: In none of the 24 reviewed studies was a proper normalization method used. Only 1 reference gene was included for normalization in 21 out of the 24 reported studies we screened, with RPS18 and GAPDH used most frequently, followed by hypoxanthine phosphoribosyltransferase 1 (HPRT1), glutathione synthetase (GSS) (hGus), ß-2 microglobin (B2M), and acidic ribosomal phosphoprotein P0 (36B4). For 2 studies the use of multiple reference genes (2 and 3) was reported, but these had not been prevalidated for expression stability in HepaRG cells. In 1 study, there was no evidence that any reference gene had been used. Current RT-qPCR gene expression studies in HepaRG cells are being performed without adequate consideration or evaluation of reference genes. Such studies can yield erroneous and biologically irrelevant results.