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.

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.

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.

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.

Targeted plasma proteomics identifies MICA and IL1R1 proteins associated with HIV-1 reservoir size.

HIV-1 reservoir shows high variability in size and activity among virally suppressed individuals. Differences in the size of the viral reservoir may relate to differences in plasma protein concentrations. We tested whether plasma protein expression levels are associated with levels of cell-associated (CA) HIV-1 DNA and RNA in 211 virally suppressed people living with HIV (PLHIV). Plasma concentrations of FOLR1, IL1R1, MICA, and FETUB showed a positive association with CA HIV-1 RNA and DNA. Moreover, SNPs in close proximity to IL1R1 and MICA genes were found to influence the levels of CA HIV-1 RNA and DNA. We found a difference in mRNA expression of the MICA gene in homozygotes carrying the rs9348866-A allele compared to the ones carrying the G allele (p < 0.005). Overall, our findings pinpoint plasma proteins that could serve as potential targets for therapeutic interventions to lower or even eradicate cells containing CA HIV-1 RNA and DNA in PLHIV.

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.

People with HIV have higher percentages of circulating CCR5+ CD8+ T cells and lower percentages of CCR5+ regulatory T cells.

CCR5 is the main HIV co-receptor. We aimed to (1) compare CCR5 expression on immune cells between people living with HIV (PLHIV) using combination antiretroviral therapy (cART) and HIV-uninfected controls, (2) relate CCR5 expression to viral reservoir size and (3) assess determinants of CCR5 expression. This cross-sectional study included 209 PLHIV and 323 controls. Percentages of CCR5+ cells (%) and CCR5 mean fluorescence intensity assessed by flow cytometry in monocytes and lymphocyte subsets were correlated to host factors, HIV-1 cell-associated (CA)-RNA and CA-DNA, plasma inflammation markers and metabolites. Metabolic pathways were identified. PLHIV displayed higher percentages of CCR5+ monocytes and several CD8+ T cell subsets, but lower percentages of CCR5+ naive CD4+ T cells and regulatory T cells (Tregs). HIV-1 CA-DNA and CA-RNA correlated positively with percentages of CCR5+ lymphocytes. Metabolome analysis revealed three pathways involved in energy metabolism associated with percentage of CCR5+ CD8+ T cells in PLHIV. Our results indicate that CCR5 is differently expressed on various circulating immune cells in PLHIV. Hence, cell-trafficking of CD8+ T cells and Tregs may be altered in PLHIV. Associations between energy pathways and percentage of CCR5+ CD8+ T cells in PLHIV suggest higher energy demand of these cells in PLHIV.

A 3'-end capture sequencing method for high-throughput targeted gene expression profiling.

Molecular phenotyping through shallow 3'-end RNA-sequencing workflows is increasingly applied in the context of large-scale chemical or genetic perturbation screens to study disease biology or support drug discovery. While these workflows enable accurate quantification of the most abundant genes, they are less effective for applications that require expression profiling of low abundant transcripts, like long noncoding RNAs (lncRNAs), or selected gene panels. To tackle these issues, we describe a workflow combining 3'-end library preparation with 3'-end hybrid capture probes and shallow RNA-sequencing for cost-effective, targeted quantification of subsets of (low abundant) genes across hundreds to thousands of samples. To assess the performance of the method, we designed a capture probe set for more than 100 mRNA and lncRNA target genes and applied the workflow to a cohort of 360 samples. When compared to standard 3'-end RNA-sequencing, 3'-end capture sequencing resulted in a more than 200-fold enrichment of target gene abundance while conserving relative intergene and intersample abundances. 3'-end RNA capture sequencing enables accurate targeted gene expression profiling at extremely shallow sequencing depth.

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.

Organ-specific genome diversity of replication-competent SARS-CoV-2.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is not always confined to the respiratory system, as it impacts people on a broad clinical spectrum from asymptomatic to severe systemic manifestations resulting in death. Further, accumulation of intra-host single nucleotide variants during prolonged SARS-CoV-2 infection may lead to emergence of variants of concern (VOCs). Still, information on virus infectivity and intra-host evolution across organs is sparse. We report a detailed virological analysis of thirteen postmortem coronavirus disease 2019 (COVID-19) cases that provides proof of viremia and presence of replication-competent SARS-CoV-2 in extrapulmonary organs of immunocompromised patients, including heart, kidney, liver, and spleen (NCT04366882). In parallel, we identify organ-specific SARS-CoV-2 genome diversity and mutations of concern N501Y, T1027I, and Y453F, while the patient had died long before reported emergence of VOCs. These mutations appear in multiple organs and replicate in Vero E6 cells, highlighting their infectivity. Finally, we show two stages of fatal disease evolution based on disease duration and viral loads in lungs and plasma. Our results provide insights about the pathogenesis and intra-host evolution of SARS-CoV-2 and show that COVID-19 treatment and hygiene measures need to be tailored to specific needs of immunocompromised patients, even when respiratory symptoms cease.

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.

The RNA Atlas expands the catalog of human non-coding RNAs.

Existing compendia of non-coding RNA (ncRNA) are incomplete, in part because they are derived almost exclusively from small and polyadenylated RNAs. Here we present a more comprehensive atlas of the human transcriptome, which includes small and polyA RNA as well as total RNA from 300 human tissues and cell lines. We report thousands of previously uncharacterized RNAs, increasing the number of documented ncRNAs by approximately 8%. To infer functional regulation by known and newly characterized ncRNAs, we exploited pre-mRNA abundance estimates from total RNA sequencing, revealing 316 microRNAs and 3,310 long non-coding RNAs with multiple lines of evidence for roles in regulating protein-coding genes and pathways. Our study both refines and expands the current catalog of human ncRNAs and their regulatory interactions. All data, analyses and results are available for download and interrogation in the R2 web portal, serving as a basis for future exploration of RNA biology and function.

The Architecture of Circulating Immune Cells Is Dysregulated in People Living With HIV on Long Term Antiretroviral Treatment and Relates With Markers of the HIV-1 Reservoir, Cytomegalovirus, and Microbial Translocation.

Long-term changes in the immune system of successfully treated people living with HIV (PLHIV) remain incompletely understood. In this study, we assessed 108 white blood cell (WBC) populations in a cohort of 211 PLHIV on stable antiretroviral therapy and in 56 HIV-uninfected controls using flow cytometry. We show that marked differences exist in T cell maturation and differentiation between PLHIV and HIV-uninfected controls: PLHIV had reduced percentages of CD4+ T cells and naïve T cells and increased percentages of CD8+ T cells, effector T cells, and T helper 17 (Th17) cells, together with increased Th17/regulatory T cell (Treg) ratios. PLHIV also exhibited altered B cell maturation with reduced percentages of memory B cells and increased numbers of plasmablasts. Determinants of the T and B cell composition in PLHIV included host factors (age, sex, and smoking), markers of the HIV reservoir, and CMV serostatus. Moreover, higher circulating Th17 percentages were associated with higher plasma concentrations of interleukin (IL) 6, soluble CD14, the gut homing chemokine CCL20, and intestinal fatty acid binding protein (IFABP). The changes in circulating lymphocytes translated into functional changes with reduced interferon (IFN)- γ responses of peripheral blood mononuclear cells to stimulation with Candida albicans and Mycobacterium tuberculosis. In conclusion, this comprehensive analysis confirms the importance of persistent abnormalities in the number and function of circulating immune cells in PLHIV on stable treatment.

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

Chronic HIV infection induces transcriptional and functional reprogramming of innate immune cells.

Chronic inflammation and immune dysfunction play a key role in the development of non-AIDS-related comorbidities. The aim of our study was to characterize the functional phenotype of immune cells in people living with HIV (PLHIV). We enrolled a cross-sectional cohort study of PLHIV on stable antiretroviral therapy and healthy controls. We assessed ex vivo cytokine production capacity and transcriptomics of monocytes and T cells upon bacterial, fungal, and viral stimulation. PLHIV exhibited an exacerbated proinflammatory profile in monocyte-derived cytokines, but not in lymphocyte-derived cytokines. Particularly, the production of the IL-1ß to imiquimod, E. coli LPS, and Mycobacterium tuberculosis was increased, and this production correlated with plasma concentrations of high-sensitivity C-reactive protein and soluble CD14. This increase in monocyte responsiveness remained stable over time in subsequent blood sampling after more than 1 year. Transcriptome analyses confirmed priming of the monocyte IL-1ß pathway, consistent with a monocyte-trained immunity phenotype. Increased plasma concentrations of ß-glucan, a well-known inducer of trained immunity, were associated with increased innate cytokine responses. Monocytes of PLHIV exhibited a sustained proinflammatory immune phenotype with priming of the IL-1ß pathway. Training of the innate immune system in PLHIV likely plays a role in long-term HIV complications and provides a promising therapeutic target for inflammation-related comorbidities.

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.

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.

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.