The Impact of Blood Sample Processing on Ribonucleic Acid (RNA) Sequencing.

In gene quantification and expression analysis, issues with sample selection and processing can be serious, as they can easily introduce irrelevant variables and lead to ambiguous results. This study aims to investigate the extent and mechanism of the impact of sample selection and processing on ribonucleic acid (RNA) sequencing. RNA from PBMCs and blood samples was investigated in this study. The integrity of this RNA was measured under different storage times. All the samples underwent high-throughput sequencing for comprehensive evaluation. The differentially expressed genes and their potential functions were analyzed after the samples were placed at room temperature for 0h, 4h and 8h, and different feature changes in these samples were also revealed. The sequencing results showed that the differences in gene expression were higher with an increased storage time, while the total number of genes detected did not change significantly. There were five genes showing gradient patterns over different storage times, all of which were protein-coding genes that had not been mentioned in previous studies. The effect of different storage times on seemingly the same samples was analyzed in this present study. This research, therefore, provides a theoretical basis for the long-term consideration of whether sample processing should be adequately addressed.

Development of a Point-of-Care Microfluidic RNA Extraction Slide for Gene Expression Diagnosis after Irradiation.

In times of war, radiological/nuclear emergency scenarios have become a reemphasized threat. However, there are challenges in transferring whole-blood samples to laboratories for specialized diagnostics using RNA. This project aims to miniaturize the process of unwieldy conventional RNA extraction with its stationed technical equipment using a microfluidic-based slide (MBS) for point-of-care diagnostics. The MBS is thought to be a preliminary step toward the development of a so-called lab-on-a-chip microfluidic device. A MBS would enable early and fast field care combined with gene expression (GE) analysis for the prediction of hematologic acute radiation syndrome (HARS) severity or identification of RNA microbes. Whole blood samples from ten healthy donors were irradiated with 0, 0.5 and 4 Gy, simulating different ARS severity degrees. RNA quality and quantity of a preliminary MBS was compared with a conventional column-based (CB) RNA extraction method. GE of four HARS severity-predicting radiation-induced genes (FDXR, DDB2, POU2AF1 and WNT3) was examined employing qRT-PCR. Compared to the CB method, twice as much total RNA from whole blood could be extracted using the MBS (6.6 ± 3.2 µg vs. 12.0 ± 5.8 µg) in half of the extraction time, and all MBS RNA extracts appeared DNA-free in contrast to the CB method (30% were contaminated with DNA). Using MBS, RNA quality [RNA integrity number equivalent (RINe)] values decreased about threefold (3.3 ± 0.8 vs. 9.0 ± 0.4), indicating severe RNA degradation, while expected high-quality RINe ≥ 8 were found using column-based method. However, normalized cycle threshold (Ct) values, as well as radiation-induced GE fold-changes appeared comparable for all genes utilizing both methods, indicating that no RNA degradation took place. In summary, the preliminary MBS showed promising features such as: 1. halving the RNA extraction time without the burden of heavy technical equipment (e.g., a centrifuge); 2. absence of DNA contamination in contrast to CB RNA extraction; 3. reduction in blood required, because of twice the biological output of RNA; and 4. equal GE performance compared to CB, thus, increasing its appeal for later semi-automatic parallel field applications.

HIV-DNA decrease during treatment in primary HIV-1 infection with three different drug regimens: Italian Network of Acute HIV Infection (INACTION) clinical trial.

As the introduction of antiretroviral therapy (ART) during primary HIV-1 infection (PHI) could restrict the establishment of HIV reservoirs, we aimed to assess the effect of three different ART regimens on HIV-DNA load in people living with HIV (PLWH), who started ART in PHI. Randomized, open-label, multicentric study, including subjects in PHI (defined as an incomplete HIV-1 Western blot and detectable plasma HIV-RNA) in the Italian Network of Acute HIV Infection cohort. Participants were randomly assigned (10:10:8) to a fixed-dose combination of tenofovir alafenamide fumarate (TAF) 10 mg plus emtricitabine (FTC) 200 mg, darunavir 800 mg, and cobicistat 150 mg once daily (group A), or TAF 25 mg plus FTC 200 mg, dolutegravir 50 mg once daily (group B), or an intensified four-drug regimen (TAF 10 mg plus FTC 200 mg, dolutegravir 50 mg, darunavir 800 mg, and cobicistat 150 mg once daily) (group C). The primary endpoint was the decrease of HIV-DNA copies/106 peripheral blood mononuclear cells (PBMCs) at weeks (W) 12 and 48. Secondary endpoints were increased in CD4+ cells and in CD4+/CD8+ ratio and percentage of PLWH reaching undetectable HIV-RNA. HIV-DNA was quantified by Droplet Digital PCR (Biorad QX100) and normalized to RPP30 reference gene. This study was registered in ClinicalTrials.gov (number NCT04225325). Among 78 participants enrolled, 30 were randomized to group 1, 28 to group 2, and 20 to group 3. At baseline, median CD4+ count was 658/µL (476-790), HIV-RNA 5.37 (4.38, 6.12) log10 copies/mL, without statistical difference in their change among groups at weeks 12 and 48 (p = 0.432 and 0.234, respectively). The trial was prematurely discontinued for slow accrual and for COVID-19 pandemic-associated restrictions. In the per-protocol analysis, PLWH (n = 72) with undetectable viral load was 54.3% at W12 and 86.4% at W48. Interestingly, the CD4/CD8 ratio progressively increased over time, up to normalization in almost half of the cohort by week 48, despite a deflection in group 3; no difference was observed by the Fiebig stage (I-III vs. IV-VI). HIV-DNA decreased from 4.46 (4.08, 4.81) log10 copies/106 PBMCs to 4.22 (3.79, 4.49) at week 12, and 3.87 (3.46, 4.34) at week 48, without difference among groups. At multivariable analysis, HIV-DNA delta at W48 was associated only with the increase of CD4+ count by 100 cells/mm3 but not with the Fiebig stage, the CD4+/CD8+ ratio, and treatment arm, despite a higher decrease in group 3. Six adverse events were recorded during our study, which did not cause any withdrawal from the study. We observed a decrease in HIV-DNA from baseline to W48 in PLWH treated during PHI, associated with an increase in CD4+ count, unrelated to the treatment arm.

Early prediction of preeclampsia in pregnancy with cell-free RNA.

Liquid biopsies that measure circulating cell-free RNA (cfRNA) offer an opportunity to study the development of pregnancy-related complications in a non-invasive manner and to bridge gaps in clinical care1-4. Here we used 404 blood samples from 199 pregnant mothers to identify and validate cfRNA transcriptomic changes that are associated with preeclampsia, a multi-organ syndrome that is the second largest cause of maternal death globally5. We find that changes in cfRNA gene expression between normotensive and preeclamptic mothers are marked and stable early in gestation, well before the onset of symptoms. These changes are enriched for genes specific to neuromuscular, endothelial and immune cell types and tissues that reflect key aspects of preeclampsia physiology6-9, suggest new hypotheses for disease progression and correlate with maternal organ health. This enabled the identification and independent validation of a panel of 18 genes that when measured between 5 and 16 weeks of gestation can form the basis of a liquid biopsy test that would identify mothers at risk of preeclampsia long before clinical symptoms manifest themselves. Tests based on these observations could help predict and manage who is at risk for preeclampsia-an important objective for obstetric care10,11.

Plasma cell-free RNA characteristics in COVID-19 patients.

The pathogenesis of COVID-19 is still elusive, which impedes disease progression prediction, differential diagnosis, and targeted therapy. Plasma cell-free RNAs (cfRNAs) carry unique information from human tissue and thus could point to resourceful solutions for pathogenesis and host-pathogen interactions. Here, we performed a comparative analysis of cfRNA profiles between COVID-19 patients and healthy donors using serial plasma. Analyses of the cfRNA landscape, potential gene regulatory mechanisms, dynamic changes in tRNA pools upon infection, and microbial communities were performed. A total of 380 cfRNA molecules were up-regulated in all COVID-19 patients, of which seven could serve as potential biomarkers (AUC > 0.85) with great sensitivity and specificity. Antiviral (NFKB1A, IFITM3, and IFI27) and neutrophil activation (S100A8, CD68, and CD63)-related genes exhibited decreased expression levels during treatment in COVID-19 patients, which is in accordance with the dynamically enhanced inflammatory response in COVID-19 patients. Noncoding RNAs, including some microRNAs (let 7 family) and long noncoding RNAs (GJA9-MYCBP) targeting interleukin (IL6/IL6R), were differentially expressed between COVID-19 patients and healthy donors, which accounts for the potential core mechanism of cytokine storm syndromes; the tRNA pools change significantly between the COVID-19 and healthy group, leading to the accumulation of SARS-CoV-2 biased codons, which facilitate SARS-CoV-2 replication. Finally, several pneumonia-related microorganisms were detected in the plasma of COVID-19 patients, raising the possibility of simultaneously monitoring immune response regulation and microbial communities using cfRNA analysis. This study fills the knowledge gap in the plasma cfRNA landscape of COVID-19 patients and offers insight into the potential mechanisms of cfRNAs to explain COVID-19 pathogenesis.

RNA profiles reveal signatures of future health and disease in pregnancy.

Maternal morbidity and mortality continue to rise, and pre-eclampsia is a major driver of this burden1. Yet the ability to assess underlying pathophysiology before clinical presentation to enable identification of pregnancies at risk remains elusive. Here we demonstrate the ability of plasma cell-free RNA (cfRNA) to reveal patterns of normal pregnancy progression and determine the risk of developing pre-eclampsia months before clinical presentation. Our results centre on comprehensive transcriptome data from eight independent prospectively collected cohorts comprising 1,840 racially diverse pregnancies and retrospective analysis of 2,539 banked plasma samples. The pre-eclampsia data include 524 samples (72 cases and 452 non-cases) from two diverse independent cohorts collected 14.5 weeks (s.d., 4.5 weeks) before delivery. We show that cfRNA signatures from a single blood draw can track pregnancy progression at the placental, maternal and fetal levels and can robustly predict pre-eclampsia, with a sensitivity of 75% and a positive predictive value of 32.3% (s.d., 3%), which is superior to the state-of-the-art method2. cfRNA signatures of normal pregnancy progression and pre-eclampsia are independent of clinical factors, such as maternal age, body mass index and race, which cumulatively account for less than 1% of model variance. Further, the cfRNA signature for pre-eclampsia contains gene features linked to biological processes implicated in the underlying pathophysiology of pre-eclampsia.

Identification of Characteristic Genes in Whole Blood of Intervertebral Disc Degeneration Patients by Weighted Gene Coexpression Network Analysis (WGCNA).

Intervertebral disc degeneration (IDD) is a major cause of lower back pain. However, to date, the molecular mechanism of the IDD remains unclear. Gene expression profiles and clinical traits were downloaded from the Gene Expression Omnibus (GEO) database. Firstly, weighted gene coexpression network analysis (WGCNA) was used to screen IDD-related genes. Moreover, least absolute shrinkage and selection operator (LASSO) logistic regression and support vector machine (SVM) algorithms were used to identify characteristic genes. Furthermore, we further investigated the immune landscape by the Cell-type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT) algorithm and the correlations between key characteristic genes and infiltrating immune cells. Finally, a competing endogenous RNA (ceRNA) network was established to show the regulatory mechanisms of characteristic genes. A total of 2458 genes were identified by WGCNA, and 48 of them were disordered. After overlapping the genes obtained by LASSO and SVM-RFE algorithms, genes including LINC01347, ASAP1-IT1, lnc-SEPT7L-1, B3GNT8, CHRNB3, CLEC4F, LOC102724000, SERINC2, and LOC102723649 were identified as characteristic genes of IDD. Moreover, differential analysis further identified ASAP1-IT1 and SERINC2 as key characteristic genes. Furthermore, we found that the expression of both ASAP1-IT1 and SERINC2 was related to the proportions of T cells gamma delta and Neutrophils. Finally, a ceRNA network was established to show the regulatory mechanisms of ASAP1-IT1 and SERINC2. In conclusion, the present study identified ASAP1-IT1 and SERINC2 as the key characteristic genes of IDD through integrative bioinformatic analyses, which may contribute to the diagnosis and treatment of IDD.

Liquid biopsies in epilepsy: biomarkers for etiology, diagnosis, prognosis, and therapeutics.

Epilepsy is one of the most common diseases of the central nervous system, impacting nearly 50 million people around the world. Heterogeneous in nature, epilepsy presents in children and adults alike. Currently, surgery is one treatment approach that can completely cure epilepsy. However, not all individuals are eligible for surgical procedures or have successful outcomes. In addition to surgical approaches, antiepileptic drugs (AEDs) have also allowed individuals with epilepsy to achieve freedom from seizures. Others have found treatment through nonpharmacologic approaches such as vagus nerve stimulation, or responsive neurostimulation. Difficulty in accessing samples of human brain tissue along with advances in sequencing technology have driven researchers to investigate sampling liquid biopsies in blood, serum, plasma, and cerebrospinal fluid within the context of epilepsy. Liquid biopsies provide minimal or non-invasive sample collection approaches and can be assayed relatively easily across multiple time points, unlike tissue-based sampling. Various efforts have investigated circulating nucleic acids from these samples including microRNAs, cell-free DNA, transfer RNAs, and long non-coding RNAs. Here, we review nucleic acid-based liquid biopsies in epilepsy to improve understanding of etiology, diagnosis, prediction, and therapeutic monitoring.

Splicing machinery is impaired in rheumatoid arthritis, associated with disease activity and modulated by anti-TNF therapy.

OBJECTIVES: To characterise splicing machinery (SM) alterations in leucocytes of patients with rheumatoid arthritis (RA), and to assess its influence on their clinical profile and therapeutic response. METHODS: Leucocyte subtypes from 129 patients with RA and 29 healthy donors (HD) were purified, and 45 selected SM elements (SME) were evaluated by quantitative PCR-array based on microfluidic technology (Fluidigm). Modulation by anti-tumour necrosis factor (TNF) therapy and underlying regulatory mechanisms were assessed. RESULTS: An altered expression of several SME was found in RA leucocytes. Eight elements (SNRNP70, SNRNP200, U2AF2, RNU4ATAC, RBM3, RBM17, KHDRBS1 and SRSF10) were equally altered in all leucocytes subtypes. Logistic regressions revealed that this signature might: discriminate RA and HD, and anti-citrullinated protein antibodies (ACPAs) positivity; classify high-disease activity (disease activity score-28 (DAS28) >5.1); recognise radiological involvement; and identify patients showing atheroma plaques. Furthermore, this signature was altered in RA synovial fluid and ankle joints of K/BxN-arthritic mice. An available RNA-seq data set enabled to validate data and identified distinctive splicing events and splicing variants among patients with RA expressing high and low SME levels. 3 and 6 months anti-TNF therapy reversed their expression in parallel to the reduction of the inflammatory profile. In vitro, ACPAs modulated SME, at least partially, by Fc Receptor (FcR)-dependent mechanisms. Key inflammatory cytokines further altered SME. Lastly, induced SNRNP70-overexpression and KHDRBS1-overexpression reversed inflammation in lymphocytes, NETosis in neutrophils and adhesion in RA monocytes and influenced activity of RA synovial fibroblasts. CONCLUSIONS: Overall, we have characterised for the first time a signature comprising eight dysregulated SME in RA leucocytes from both peripheral blood and synovial fluid, linked to disease pathophysiology, modulated by ACPAs and reversed by anti-TNF therapy.

Elevated N6-Methyladenosine RNA Levels in Peripheral Blood Immune Cells: A Novel Predictive Biomarker and Therapeutic Target for Colorectal Cancer.

Effective biomarkers for the diagnosis of colorectal cancer (CRC) are essential for improving prognosis. Imbalance in regulation of N6-methyladenosine (m6A) RNA has been associated with a variety of cancers. However, whether the m6A RNA levels of peripheral blood can serve as a diagnostic biomarker for CRC is still unclear. In this research, we found that the m6A RNA levels of peripheral blood immune cells were apparently elevated in the CRC group compared with those in the normal controls (NCs) group. Furthermore, the m6A levels arose as CRC progressed and metastasized, while these levels decreased after treatment. The area under the curve (AUC) of the m6A levels was 0.946, which was significantly higher than the AUCs for carcinoembryonic antigen (CEA; 0.817), carbohydrate antigen 125 (CA125; 0.732), and carbohydrate antigen 19-9 (CA19-9; 0.771). Moreover, the combination of CEA, CA125, and CA19-9 with m6A levels improved the AUC to 0.977. Bioinformatics and qRT-PCR analysis further confirmed that the expression of m6A modifying regulator IGF2BP2 was markedly elevated in peripheral blood of CRC patients. Gene set variation analysis (GSVA) implied that monocyte was the most abundant m6A-modified immune cell type in CRC patients' peripheral blood. Additionally, m6A modifications were negatively related to the immune response of monocytes. In conclusion, our results revealed that m6A RNA of peripheral blood immune cells was a prospective non-invasive diagnostic biomarker for CRC patients and might provide a valuable therapeutic target.

Exploration of differentially-expressed exosomal mRNAs, lncRNAs and circRNAs from serum samples of gallbladder cancer and xantho-granulomatous cholecystitis patients.

Gallbladder cancer (GBC) is the most common biliary tract malignancy worldwide. Although a growing number of studies have explored the mechanism of GBC, thus far, few molecules have been discovered that can be utilized as specific biomarkers for the early diagnosis and therapeutic treatment of GBC. Recent studies have shown that exosomes not only participate in the progression of tumors, but also carry specific information that can define multiple cancer types. The present study investigated the expression profiles of coding (or messenger) ribonucleic acids (mRNAs) and non-coding RNAs (ncRNAs, including long non-coding RNAs [lncRNAs] and circular RNAs [circRNAs]) in plasma-derived exosomes from GBC patients. Using high-throughput RNA sequencing and subsequent bioinformatic analysis, a number of differentially expressed (DE) mRNAs, lncRNAs, and circRNAs were identified in GBC exosomes, compared to their expressions in xantho-granulomatous cholecystitis (XGC) exosomes. Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genome (KEGG) analyses were then conducted to investigate the potential functions of these DE RNAs. Furthermore, the interaction networks and competing endogenous RNA networks of these DE RNAs and their target genes were investigated, revealing a complex regulatory network among mRNAs and ncRNAs. In summary, this study demonstrates the diagnostic value of plasma-derived exosomes in GBC and provides a new perspective on the mechanism of GBC.

Assembly of long L-RNA by native RNA ligation.

Due to their intrinsic nuclease resistance, L-oligonucleotides are being increasingly utilized in the development of molecular tools and sensors. Yet, it remains challenging to synthesize long L-oligonucleotides, potential limiting future applications. Herein, we report straightforward and versitile approach to assemble long L-RNAs from two or more shorter fragments using T4 RNA ligase 1. We show that this approach is compatible with the assembly of several classes of functional L-RNA, which we highlight by generating a 124 nt L-RNA biosensor that functions in serum.

Eight-year longitudinal study of whole blood gene expression profiles in individuals undergoing long-term medical follow-up.

Blood circulates throughout the body via the peripheral tissues, contributes to host homeostasis and maintains normal physiological functions, in addition to responding to lesions. Previously, we revealed that gene expression analysis of peripheral blood cells is a useful approach for assessing diseases such as diabetes mellitus and cancer because the altered gene expression profiles of peripheral blood cells can reflect the presence and state of diseases. However, no chronological assessment of whole gene expression profiles has been conducted. In the present study, we collected whole blood RNA from 61 individuals (average age at registration, 50 years) every 4 years for 8 years and analyzed gene expression profiles using a complementary DNA microarray to examine whether these profiles were stable or changed over time. We found that the genes with very stable expression were related mostly to immune system pathways, including antigen cell presentation and interferon-related signaling. Genes whose expression was altered over the 8-year study period were principally involved in cellular machinery pathways, including development, signal transduction, cell cycle, apoptosis, and survival. Thus, this chronological examination study showed that the gene expression profiles of whole blood can reveal unmanifested physiological changes.

Optimized protocol for the extraction of RNA and DNA from frozen whole blood sample stored in a single EDTA tube.

Cryopreservation of whole blood is useful for DNA collection, and clinical and basic research. Blood samples in ethylenediaminetetraacetic acid disodium salt (EDTA) tubes stored at - 80 °C are suitable for DNA extraction, but not for high-quality RNA extraction. Herein, a new methodology for high-quality RNA extraction from human blood samples is described. Quickly thawing frozen whole blood on aluminum blocks at room temperature could minimize RNA degradation, and improve RNA yield and quality compared with thawing the samples in a 37 °C water bath. Furthermore, the use of the NucleoSpin RNA kit increased RNA yield by fivefold compared with the PAXgene Blood RNA Kit. Thawing blood samples on aluminum blocks significantly increased the DNA yield by ~ 20% compared with thawing in a 37 °C water bath or on ice. Moreover, by thawing on aluminum blocks and using the NucleoSpin RNA and QIAamp DNA Blood kits, the extraction of RNA and DNA of sufficient quality and quantity was achieved from frozen EDTA whole blood samples that were stored for up to 8.5 years. Thus, extracting RNA from frozen whole blood in EDTA tubes after long-term storage is feasible. These findings may help advance gene expression analysis, as well as biomarker research for various diseases.

A blood RNA transcriptome signature for COVID-19.

BACKGROUND: COVID-19 is a respiratory viral infection with unique features including a more chronic course and systemic disease manifestations including multiple organ involvement; and there are differences in disease severity between ethnic groups. The immunological basis for disease has not been fully characterised. Analysis of whole-blood RNA expression may provide valuable information on disease pathogenesis. METHODS: We studied 45 patients with confirmed COVID-19 infection within 10 days from onset of illness and a control group of 19 asymptomatic healthy volunteers with no known exposure to COVID-19 in the previous 14 days. Relevant demographic and clinical information was collected and a blood sample was drawn from all participants for whole-blood RNA sequencing. We evaluated differentially-expressed genes in COVID-19 patients (log2 fold change ≥ 1 versus healthy controls; false-discovery rate < 0.05) and associated protein pathways and compared these to published whole-blood signatures for respiratory syncytial virus (RSV) and influenza. We developed a disease score reflecting the overall magnitude of expression of internally-validated genes and assessed the relationship between the disease score and clinical disease parameters. RESULTS: We found 135 differentially-expressed genes in the patients with COVID-19 (median age 35 years; 82% male; 36% Chinese, 53% South Asian ethnicity). Of the 117 induced genes, 14 were found in datasets from RSV and 40 from influenza; 95 genes were unique to COVID-19. Protein pathways were mostly generic responses to viral infections, including apoptosis by P53-associated pathway, but also included some unique pathways such as viral carcinogenesis. There were no major qualitative differences in pathways between ethnic groups. The composite gene-expression score was correlated with the time from onset of symptoms and nasal swab qPCR CT values (both p < 0.01) but was not related to participant age, gender, ethnicity or the presence or absence of chest X-ray abnormalities (all p > 0.05). CONCLUSIONS: The whole-blood transcriptome of COVID-19 has overall similarity with other respiratory infections but there are some unique pathways that merit further exploration to determine clinical relevance. The approach to a disease score may be of value, but needs further validation in a population with a greater range of disease severity.

Fit-for-purpose quantitative liquid biopsy based droplet digital PCR assay development for detection of programmed cell death ligand-1 (PD-L1) RNA expression in PAXgene blood samples.

Development of a clinically applicable liquid biopsy-based test for PD-L1 mRNA expression would be beneficial in providing complementary evidence to current immunohistochemistry assays. Hence, we report the development of a fit-for-purpose assay for detection of blood PD-L1 mRNA expression using droplet digital polymerase chain reaction (ddPCR). TaqMan® assays were selected based on coverage of the PD-L1 gene and were tested for linearity and efficiency using real-time quantitative PCR. Four reference genes were analyzed in positive control cell line (A549 treated with interferon gamma, [IFN γ]) genomic DNA. The PD-L1 primer/probe sets were also evaluated in ddPCR for limit of blank, limit of detection, and precision. Finally, thirty-five healthy volunteer samples were evaluated to establish a baseline level of PD-L1 expression. In ddPCR, the limit of blank was determined to be 0 copies and the limit of detection was determined to be less than or equal to 19 copies of PD-L1. The average intra-run coefficient of variation in the ddPCR assay was 7.44% and average inter-run CV was 7.70%. Treatment of A549 cells with IFN γ resulted in a 6.7-fold increase in PD-L1 expression (21,580 copies in untreated cDNA versus 145,000 copies in treated cDNA). Analysis of healthy human samples yielded a median value of 1659 PD-L1 copies/µL with a range of 768-7510 copies/µL. The assay was transferred to an external service provider and results from our in-house experiments and those conducted externally shows a correlation of 0.994. In conclusion, a fit-for-purpose liquid biopsy-based, purely quantitative ddPCR assay for the detection of PD-L1 mRNA expression was developed and validated using PAXgene RNA blood samples. Linearity, reproducibility, limit of blank and limit of detection were measured and deemed suitable for clinical application. This ultra-sensitive liquid biopsy ddPCR assay has promising clinical potential in screening, longitudinal monitoring and disease progression detection.