Skip to main content

Relationship between clinical features and droplet digital PCR copy number in non-HIV patients with pneumocystis pneumonia



Droplet digital PCR (ddPCR) is a novel assay to detect pneumocystis jjrovecii (Pj) which has been defined to be more sensitive than qPCR in recent studies. We aimed to explore whether clinical features of pneumocystis pneumonia (PCP) were associated with ddPCR copy numbers of Pj.


A total of 48 PCP patients were retrospectively included. Pj detection was implemented by ddPCR assay within 4 h. Bronchoalveolar fluid (BALF) samples were collected from 48 patients with molecular diagnosis as PCP via metagenomic next generation sequencing (mNGS) or quantitative PCR detection. Univariate and multivariate logistic regression were performed to screen out possible indicators for the severity of PCP. The patients were divided into two groups according to ddPCR copy numbers, and their clinical features were further analyzed.


Pj loading was a pro rata increase with serum (1,3)-beta-D glucan, D-dimmer, neutrophil percentage, procalcitonin and BALF polymorphonuclear leucocyte percentage, while negative correlation with albumin, PaO2/FiO2, BALF cell count, and BALF lymphocyte percentage. D-dimmer and ddPCR copy number of Pj were independent indicators for moderate/severe PCP patients with PaO2/FiO2 lower than 300. We made a ROC analysis of ddPCR copy number of Pj for PaO2/FiO2 index and grouped the patients according to the cut-off value (2.75). The high copy numbers group was characterized by higher level of inflammatory markers. Compared to low copy number group, there was lower level of the total cell count while higher level of polymorphonuclear leucocyte percentage in BALF in the high copy numbers group. Different from patients with high copy numbers, those with high copy numbers had a tendency to develop more severe complications and required advanced respiratory support.


The scenarios of patients infected with high ddPCR copy numbers of Pj showed more adverse clinical conditions. Pj loading could reflect the severity of PCP to some extent.

Strengths and Limitations

This is the first study to investigate the relationship between ddPCR copy number and the clinical characteristics of patients with PCP.

Droplet digital PCR was a more sensitve and accurate method than qPCR for detection of pathogens.

This was a retrospective study and we did not include more samples for it.

Peer Review reports


Pneumocystis jirovecii (Pj) is a ubiquitous opportunistic fungus that can cause life-threatening pneumocystis pneumonia (PCP) in immunocompromised hosts [1]. Owing to the application of highly active immunosuppressive therapy, PCP has become increasingly common among non-HIV-infected patients with hematologic disorders or autoimmune diseases (HIV- human immunodeficiency virus) [2, 3]. Patients may present with fever, cough, and dyspnea when infected with Pj [4, 5]. Severe cases can lead to respiratory failure and even death [6]. Therefore, rapid and accurate detection is critical in patients with PCP. Pathological diagnosis is the gold standard for PCP, and microscopy is far from optimal because of the high probability of false negatives. In recent years, molecular detection methods have gradually been developed to detect Pj. Polymerase chain reaction (PCR) assays, which mainly target the mitochondrial large subunit rRNA gene and the multicopy major surface glycoprotein gene family, make the detection of Pj more sensitive than microscopic detection [7, 8]. However, no commercial PCR kits have been approved for clinical use in China, and their sensitivity is unsatisfactory [9]. Metagenomic next-generation sequencing (mNGS) is a high-throughput molecular technology that can achieve unbiased pathogen detection in a single assay. Due to its high sensitivity, it has become increasingly popular in recent years [8] and can be used to identify different pathogens involved in infectious diseases [10, 11]. However, its use has been limited by the high cost.

Droplet digital PCR (ddPCR), based on water-oil emulsion droplet technology, is a novel technique that can accurately quantify trace specimens [12, 13]. The number of droplets was counted after amplification, according to which an estimation of template concentration was made based on Poisson’s Law of Small Numbers [14]. Compared with conventional molecular methods, ddPCR has the advantages of higher sensitivity and specificity, absolute quantification without a standard curve, better reproducibility, and better resistance to PCR inhibitors [13]. Many studies using ddPCR assays have shown high sensitivity and precision [15,16,17]. Currently, several studies have focused on ddPCR detection of Pj and achieved excellent results [14, 18]. However, data on ddPCR copy numbers and clinical features of PCP patients are limited. We hypothesized that the clinical manifestations in patients with PCP are related to ddPCR copy numbers. In this study, we used ddPCR to accurately quantify copy numbers of Pj in BALF samples of PCP patients. We analyzed the ddPCR copy numbers and clinical features of PCP in order to determine the relationship between them.

Materials and methods

Study participants

This was a retrospective study performed on 48 non-HIV patients with clinical PCP diagnosis from September 2016 to June 2022 at Peking University People’s Hospital. Bronchoscopy was conducted within the initial five days of anti-PCP treatment commencement. The BALF specimens were all collected from patients during their hospitalization period, as part of clinical indications for bronchoscopy examinations, and were the remaining samples after necessary laboratory analysis for submission. The study was approved by the Institutional Review Board of Peking University People’s Hospital (No. 2022PHB457). The included patients needed to meet the following criteria: [1] Have at least one of the following severe immunosuppression conditions: primary immunodeficiencies, solid organ or bone marrow transplantation, hematologic neoplasia, or other diseases that cause deterioration in cellular immunity, including but not limited to radiotherapy and chemotherapy for cancer, autoimmune diseases of the connective tissue, asthma, and acute exacerbations chronic obstructive pulmonary disease (AECOPD), which required treatment with immunosuppressive medications such as prednisone at a dose of 0.3 mg/kg/day or its equivalent for a duration exceeding two weeks [2]. Patients who were diagnosed as PCP; [3] Patients with both retained BALF specimens and molecular testing results. Patients whose imaging and clinical manifestations do not correspond to PCP, which can be explained by other pathogens or colonization are excluded from the study. Forty patients who didn’t perform molecular detection (24 patient) or whose BALF specimen weren’t retained (16 patients) were excluded (Fig. 1).


The clinical diagnosis of PCP requires meeting at least two of the following three criteria: respiratory symptoms including cough and/or dyspnoea, hypoxia, typical radiological picture such as ground glass opacity on a chest computer tomography scan, diffuse interstitial opacity on chest x-ray; receipt of a full course of treatment for PCP. A PCP case was considered if a patient meets the clinical diagnostic criteria and is either PCR or mNGS positive, or if they meet the clinical diagnostic criteria and show response to empirical treatment even with negative molecular testing. Colonization were described not fulfilling the clinical diagnostic criteria for PCP with positive PCR [8].

Data collection

We recorded the following clinical and biological data for each patient: molecular detection results (mNGS or qPCR), underlying disease (hematological malignancies, autoimmune disease, renal diseases, cancer, and other diseases), vital signs, comorbidities, radiological signs (obtained by X-ray analysis or computed tomography scan), laboratory tests (routine blood tests, liver function, renal function, electrolytes, PCT, CRP, coagulation function, PaO2/FiO2, serum (1,3)-beta-D glucan (BDG), cytologic analysis for BALF), treatments and outcomes (Hospital stay, immunosuppressive medication, complications, oxygen support assays, and clinical outcome). Patients who had incomplete medical records, lost BALF samples or no molecular detection were excluded. BALF specimens were collected for qPCR or mNGS detection. Arterial blood gas results were collected on the same day of tracheoscopy. BDG results were included before or after BALF specimen were collected within 7 days. The requirement for written informed consent was waived off due to the retrospective nature of the study.

Fig. 1
figure 1

Study flowchart

DNA extraction

BALF specimens from patients were collected according to standard operating procedures. For DNA extraction, we used 1.5mL microcentrifuge tubes to contain a 0.6mL sample attached to a horizontal platform on a vortex mixer. The samples were agitated vigorously at 2,800–3,200 rpm for 15 min. Then they were centrifuged and the supernatant were discarded. Cells were resuspended in 100uL nucleic acid extraction reagent and transferred into another microcentrifuge tube with 1 g 0.5 mm glass beads. The samples were agitated vigorously at 2,800–3,200 rpm for 8 min and then had metal bath for 5 min at 95 °C. Nucleic Acid Extraction Kit (Capital Bio Technology) was used according to the manufacturer’s recommendation. DNA concentrations were determined with NanoDrop ONEc spectrophotometer (Thermo Fisher Scientific, USA).

ddPCR detection of Pj

We used ddPCR via the Targeting One Digital PCR System (Targeting One, Beijing, China) for Pj detection. The primers and probes were designed according to previous research [14]. The sequences used were the mtLSU reference sequence of Pj and the human reference gene N-acetylglucosamine kinase (NAGK). The Forward primer of mtLSU reference sequence of Pj was GTATAGCACTGAATATCTCGAGGG. The Forward primer of mtLSU reference sequence of Pj was.

GAGCTTTAATTACTGTTCTGGGCT. The sequence of probe was FAM-TTCGACTATCTACCTTATCGC -MGB. The Forward primer of NAGK was AGATGCTGGGCAGACACATC. The Forward primer of NAGK was CCCACCTTCACTCCCACCT. The sequence of probe was HEX-AGCAGTGTTGCCCGAGATTG ACCC-BHQ1. The 30 µL PCR reaction mixture included 7.5 µL 4× SuperMix, 3 µL 10× primers and probes mix (the final concentrations of primers and probes were 600 and 300 nM, respectively), and up to 19.5 µL of the DNA template. We add 800ng genomic DNA sample as the template for PCR in the 30ul reaction mixture. The rest procedures were completed as described in the previous study. The reaction mixture of 30 µl was transferred into an 8-strip PCR tube and subjected to thermal cycling with the following conditions: initial denaturation at 95 °C for 10 min, followed by 45 amplification cycles involving denaturation at 94 °C for 30 s and annealing at 57 °C for 1 min. The reaction concluded with a final step at 12 °C. Following amplification, the amplified droplets within the 8-strip PCR tube were connected to a droplet detection chip and subsequently analyzed using the Chip Reader. When the copy number of the vector DNA found in that sample was greater than 5, the sample was considered positive and the measured copy number was reported.


After DNA extraction from BALF samples, PCR detection was performed. The specific detection process and synthesis of primers were based on previous literature [19]. The NGS detection process was also conducted following the procedures outlined in previous literature [20].

Statistical analysis

The severity of PCP was assessed based on the PaO2/FiO2 index when BALF specimen were collected and was categorized as mild (≥ 300 mmHg) or severe (< 300 mmHg) pneumonia. A log10-transformation (log10 [ddPCR copy number] and log10 [sequence number]) was applied to normalize the data. We made the univariate logistic regression analysis and found ddPCR copy number to be a risk factor for the severity of PCP according to PaO2/FiO2 index. Then we made a receiver operating characteristic (ROC) analysis and calculated the cutoff value. Finally, the patients were grouped into low copy numbers (n = 24) and high copy numbers (n = 24) groups according to the cutoff value. Finally, Continuous variables in a normal distribution were presented as the mean (± S.D.) while the median (25th–75th percentile) in a non-normal distribution. Continuous data with a non-normal distribution were analyzed by Wilcoxon rank-sum test while T test with a normal distribution. Categorical variables were analyzed by chi-squared test. A p-value < 0.05 was considered statistically significant. The correlations of different markers and ddPCR copy numbers were determined using the Spearman’s correlation for nonparametric data. As for sample size, we set 1-α to be 0.8, AUC to be 0.75, and Confidence Interval Width to be 0.2. Using he PASS 15 software, we determined that a sample size of 40 is necessary.


Baseline characteristics of PCP patients

A total of 48 patients were enrolled in this study. As shown in Table 1, the median age of the subjects was 53 years, and 52% were men. The primary diseases were autoimmune and hematological diseases, accounting for 90% of all cases. Approximately 40% of the patients had comorbidities such as diabetes, hypertension, hyperlipidemia, and cardiovascular disease. Most patients had cough, fever, and dyspnea. Patients with moderate-to-severe PCP had fewer cardiovascular diseases and higher heart rates and higher ddPCR copy numbers. There were no significant differences between the mild and moderate/severe PCP patients regarding sex, age, primary diseases, fever, highest temperature, cough, and dyspnea.

Table 1 Baseline characteristics of PCP patients

Assessment of the ddPCR-based detection of PCP

First, we performed the ddPCR assay. Twenty-three specimens were detected using mNGS, while 25 were detected using PCR. Among the 23 patients who underwent mNGS assays, two patients had negative ddPCR results, which were positive on mNGS, and one had a negative mNGS and positive ddPCR. Of the 25 patients who underwent PCR assays, two patients had negative results from the PCR assay but positive ddPCR results and one had a positive result from the PCR assay and a negative ddPCR. The detailed data are presented in Table 2. Eleven patients had data on mNGS sequence numbers, and we performed a correlation analysis with ddPCR copy numbers. As shown in Fig. 2, there was a significant positive correlation between the two.

Table 2 The discrepancies among different assays
Fig. 2
figure 2

Correlation analysis between the mNGS sequence number and the Pj load of ddPCR in BALF. The mNGS sequence number had a good correlation with the Pj load of ddPCR in BALF.

Correlations between the ddPCR copy number and clinical indicators

We performed a correlation analysis between copy numbers and different indicators (Fig. 3). The copy number was positively correlated with serum [1, 3]-β-D glucan, d-dimer, neutrophil percentage, procalcitonin (PCT), and BALF polymorphonuclear leukocyte percentage while negatively correlated with albumin, PaO2/FiO2, BALF cell count, and BALF lymphocyte percentage.

Fig. 3
figure 3

Correlation of different clinical indicators with the ddPCR copy number in 48 patients of PCP. We found that the ddPCR copy number had significant correlation with neutrophil percentage (a), PCT (b), ALB (c), d-dimmer (d), BDG (e), PaO2/FiO2 (f), BALF total cell count (g), BALF lymphocyte percentage (h), and BALF locbocyte percentage (i). MPN, polymorphonuclear leucocyte

Statistical analysis of risk factors for PCP severity and ROC analysis of ddPCR copy numbers

Independent indicators were introduced into univariate logistic regression analysis (Table 3). The results showed that high levels of neutrophil percentage, d-dimer, and lg (copy number) were associated with significantly high-risk ratios of low PaO2/FiO2, while high levels of serum albumin (ALB) were associated with significantly low-risk ratios. When the significant independent variables were integrated into a multivariate logistic regression analysis, only D-dimer [hazard ratio (95% confidence interval (CI)):1.004 (1.000-1.009)] and lg (copy number) [hazard ratio (95% CI):1.787 (1.036–3.082)] were found to be significant independent predictors of PCP severity (p = 0.034 and p = 0.037, respectively). To investigate whether the ddPCR copy numbers were related to the clinical features of PCP, we performed a ROC analysis between lg (copy number) and PaO2/FiO2 index (Fig. 4). The AUC was 0.794. Finally, we determined the cut-off value of lg (copy number) to be 2.85.

Fig. 4
figure 4

ROC curve analysis of ddPCR copy number for severity of PCP. AUC was 0.794 and cut-off value of lg (copy number) was 2.85

Table 3 Univariate and multivariate analysis of risk factors for PCP severity

Comparison of clinical manifestions in high and low ddPCR copy number groups

The 48 patients enrolled from September 2016 to June 2022 were divided into low and high copy number groups according to the cut-off value of the ROC analysis. Table 4 shows the baseline characteristics and laboratory results of the patients with PCP with low or high ddPCR copy numbers. More patients with other diseases, such as tumors or glomerular diseases, were in the high-copy number group. No statistical difference was observed between the two groups in terms of age, sex, comorbidities, temperature, heart rate, cough, or dyspnea. The high copy number group had higher levels of inflammatory markers, such as neutrophil percentage and PCT. Other indicators in blood showed that the high copy number group exhibited lower levels of platelet count and ALB but higher levels of BDG. The PaO2/FiO2 index was calculated and found to be significantly lower in the high copy number group. Compared to the low copy number group, there were lower total cell counts and higher levels of polymorphonuclear leukocytes (PMN) in BALF in the high copy number group. Next, we compared the treatment outcomes between the two groups. Compared to patients with high copy numbers, those with low copy numbers tend to develop more severe complications and require advanced respiratory support. Although the difference was not statistically significant, the median length of hospitalization was longer in the high copy number group. We did not find any difference in immunosuppressive therapy methods between the two groups.

Table 4 Clinical manifestions of PCP patients with low and high copies in hospital


To our knowledge, this is the first study to investigate the relationship between ddPCR copy number and the clinical characteristics of patients with PCP. In this study, we found that ddPCR copy number was closely related to the clinical features of PCP and could reflect the severity of PCP to an extent.

As a novel assay, ddPCR has revealed good diagnostic performance for PCP. A study compared ddPCR with quantitative PCR and found the same results for 95.12% of BALF samples, while the remaining samples tested positive using ddPCR but negative using qPCR. mNGS results corroborated these inconsistent results. The authors also found that the cycle threshold value of qPCR was highly correlated with the ddPCR copy number (R2 = 0.84) [14]. In another study, using a cut-off value of > 4.5 copies/uL, the ddPCR assay exhibited a sensitivity of 100% and specificity of 94.3% in non-HIV patients with PCP [18]. In our study, mNGS yielded highly consistent results for BALF specimens (R2 = 0.81). This indicated the diagnostic accuracy of the ddPCR system. Three samples tested negative by ddPCR, and the mNGS sequence numbers of the two were 1 and 3, respectively. This implied that the samples had a low fungal load. Considering the small sample size of BALF and the long storage time of these specimens, the expected effect of ddPCR could not be detected. Despite this, the novel assay exhibited a promising performance.

Since ddPCR enables precise quantification of nucleic acids, we performed Spearman correlation analysis between copy numbers and different clinical indicators. Neutrophil percentage, PCT, and d-dimer levels were positively correlated, while albumin was negatively correlated with ddPCR copy number. Neutrophils are important factors in pulmonary inflammation and can lead to lung damage when excessively activated [21]. PCT is a useful marker for discriminating between viral and bacterial pneumonia [22]. D-dimer can indicate excessive inflammation and vascular endothelial injury [23]. Hypoalbuminemia can reflect the severity of inflammation and is an independent risk factor for mortality [24]. Therefore, the high fungal burden in BALF was associated with severe inflammation and poor prognosis in patients with PCP to some extent. The oxygenation index was negatively correlated with the ddPCR copy number in our results. A low oxygenation index indicated a more severe clinical condition. This means that copy number was correlated with the severity of PCP. BDG is an important component of the cell wall of most fungi, including Pj, making it a significant diagnostic indicator of PCP [25]. Some studies have investigated whether serum BDG could be used to estimate the fungal burden, but the results were inconsistent in people with different diseases [26, 27]. A study used qPCR assays to detect Pj in BALF, performed correlation analysis for Ct value and serum BDG, and found a weak correlation between them [26]. From our findings, ddPCR copy number significantly correlated with BDG. This was because the ddPCR assay accurately detected the fungal load. The BALF total cell count and BALF lymphocyte percentage were negatively correlated with the ddPCR copy number, while the BALF PMN percentage was positively correlated with it. PCP patients are prone to a decrease in the total cell count and lymphocytes due to the use of immunosuppressive agents. Reduced CD8 + T cell levels can be used to predict the adverse development of PCP in immunosuppressed hosts [28]. The BALF PMN percentage increased due to the decreased lymphocyte percentage.

PCP is a life-threatening condition in immunosuppressed hosts. Many studies have investigated indicators for predicting the prognosis of PCP. One study found that platelet count and γ-globulin could be risk factors in the HIV group in the multivariate analysis, while albumin could be a risk factor in the non-HIV group in the univariate analysis [29]. In another study, mixed infection and elevation of NLR were considered independent risk factors for poor prognosis of non-HIV PCP patients diagnosed by mNGS [10]. We performed univariate and multivariate regression analyses for PaO2/FiO2 index and found ddPCR copy number and D-dimer to be risk factors. In contrast to previous studies with outcome measures of mortality, we defined an oxygenation index before tracheoscopy detection with less than 300 as the outcome. The two outcome indicators were consistent. For example, plasma interleukin (IL)-6 and IL-8 levels were significantly higher in PCP patients with lower oxygenation as well as in non-survivors [30]. Studies have evaluated the relationship between Ct-value and PCP severity but found no significant difference between the severe and mild PCP groups. Nonetheless, a trend of higher Ct values in the mild group was found [8]. The accurate quantification of ddPCR enabled us to further investigate their relationship. Individuals with a low PaO2/FiO2 index had significantly higher ddPCR copy numbers. Further analysis revealed that a high fungal load in the BALF can be an independent risk factor for disease severity. To deeply explore the relationship between ddPCR copy number and clinical features, we performed ROC analysis for the copy number and divided the patients into two groups according to the cut-off value. We normalized the data by applying a log10-transformation (log10-[ddPCR copy number]). Differences in laboratory results and treatments were observed between the groups. High copy numbers were associated with more severe clinical presentations, such as higher inflammatory markers, more advanced oxygen support, and more severe complications. Based on the above findings, we concluded that the ddPCR copy number could reflect the severity of PCP. The treatment outcomes also showed differences between the two groups. The high copy number group had a longer hospital stay and more severe complications than the low copy number group. Univariate analysis showed that mechanical ventilation could be a risk factor for patient mortality [31]. These results indicate that the clinical condition was worse in the high copy number group. Therefore, the ddPCR copy number may help to evaluate clinical severity in patients with PCP. Categorizing patients according to copy numbers offers valuable insights for enhancing the management of PCP cases. Those with elevated copy numbers necessitate swift and comprehensive anti-PCP treatment, coupled with vigilant monitoring of abnormal laboratory parameters. Furthermore, their susceptibility to advancing into respiratory failure underscores the critical importance of close surveillance. We can include more patients to further confirm this finding and investigate whether the ddPCR copy number of Pj in BALF can reflect disease mortality in future studies.

There were several limitations to our study, along with opportunities for further studies. Firstly, the sample size was not sufficiently large. Secondly, not all patients are eligible for bronchoscopy. This study has exclusively explored BALF specimens, but future investigations can encompass additional respiratory specimens such as sputum. This would allow for broader applicability of our findings. Thirdly, it’s worth noting that in this study, the majority of patients initiated empirical treatment before undergoing bronchoscopy. This factor could potentially influence our quantitative results.


In conclusion, this study addressed the clinical features and ddPCR copy numbers in patients with PCP. Our data suggest that ddPCR copy number has a good correlation with many laboratory indicators such as albumin, PCT, d-dimer, BDG, and BALF lymphocytes and can, to an extent, reflect PCP severity.

Data Availability

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding authors.



Bronchoalveolar fluid


Droplet digital PCR

Pj :

Pneumocystis jjrovecii


Metagenomic next generation sequencing


Pneumocystis pneumonia


Polymerase chain reaction






Confidence interval


(1,3)-beta-D glucan


Polymorphonuclear leukocytes


C-reactive protein




Aspartate transaminase


  1. Fortea JI, Cuadrado A, Puente Á, Álvarez Fernández P, Huelin P, Álvarez Tato C, et al. Is routine Prophylaxis Against Pneumocystis jirovecii needed in liver transplantation? A retrospective single-centre experience and current prophylaxis strategies in Spain. J Clin Med. 2020;9(11).

  2. Liu CJ, Lee TF, Ruan SY, Yu CJ, Chien JY, Hsueh PR. Clinical characteristics, treatment outcomes, and prognostic factors of Pneumocystis pneumonia in non-HIV-infected patients. Infect Drug Resist. 2019;12:1457–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Lee WS, Hsueh PR, Hsieh TC, Chen FL, Ou TY, Jean SS. Caspofungin salvage therapy in Pneumocystis jirovecii pneumonia. J Microbiol Immunol Infect. 2017;50(4):547–8.

    Article  CAS  PubMed  Google Scholar 

  4. Brakemeier S, Pfau A, Zukunft B, Budde K, Nickel P. Prophylaxis and treatment of Pneumocystis Jirovecii pneumonia after solid organ transplantation. Pharmacol Res. 2018;134:61–7.

    Article  CAS  PubMed  Google Scholar 

  5. Kim TO, Lee JK, Kwon YS, Kim YI, Lim SC, Kim MS, et al. Clinical characteristics and prognosis of patients with pneumocystis jirovecii pneumonia without a compromised illness. PLoS ONE. 2021;16(2):e0246296.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Enomoto T, Azuma A, Kohno A, Kaneko K, Saito H, Kametaka M, et al. Differences in the clinical characteristics of Pneumocystis jirovecii pneumonia in immunocompromized patients with and without HIV infection. Respirology. 2010;15(1):126–31.

    Article  PubMed  Google Scholar 

  7. Salsé M, Mercier V, Carles MJ, Lechiche C, Sasso M. Performance of the RealStar(®) pneumocystis jirovecii PCR kit for the diagnosis of Pneumocystis pneumonia. Mycoses. 2021;64(10):1230–7.

    Article  CAS  PubMed  Google Scholar 

  8. Hammarström H, Grankvist A, Broman I, Kondori N, Wennerås C, Gisslen M, et al. Serum-based diagnosis of Pneumocystis pneumonia by detection of Pneumocystis jirovecii DNA and 1,3-β-D-glucan in HIV-infected patients: a retrospective case control study. BMC Infect Dis. 2019;19(1):658.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Del Corpo O, Butler-Laporte G, Sheppard DC, Cheng MP, McDonald EG, Lee TC. Diagnostic accuracy of serum (1–3)-β-D-glucan for Pneumocystis jirovecii pneumonia: a systematic review and meta-analysis. Clin Microbiol Infect. 2020;26(9):1137–43.

    Article  CAS  PubMed  Google Scholar 

  10. Duan J, Gao J, Liu Q, Sun M, Liu Y, Tan Y, et al. Characteristics and prognostic factors of Non-HIV immunocompromised patients with Pneumocystis Pneumonia diagnosed by Metagenomics Next-Generation sequencing. Front Med (Lausanne). 2022;9:812698.

    Article  PubMed  Google Scholar 

  11. Goldberg B, Sichtig H, Geyer C, Ledeboer N, Weinstock GM. Making the Leap from Research Laboratory to Clinic: Challenges and Opportunities for Next-Generation sequencing in Infectious Disease Diagnostics. mBio. 2015;6(6):e01888–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Gaňová M, Zhang H, Zhu H, Korabečná M, Neužil P. Multiplexed digital polymerase chain reaction as a powerful diagnostic tool. Biosens Bioelectron. 2021;181:113155.

    Article  CAS  PubMed  Google Scholar 

  13. Kojabad AA, Farzanehpour M, Galeh HEG, Dorostkar R, Jafarpour A, Bolandian M, et al. Droplet digital PCR of viral DNA/RNA, current progress, challenges, and future perspectives. J Med Virol. 2021;93(7):4182–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Yi J, Wang N, Wu J, Tang Y, Zhang J, Zhu L, et al. Development of a Droplet Digital polymerase chain reaction for sensitive detection of Pneumocystis jirovecii in respiratory tract specimens. Front Med (Lausanne). 2021;8:761788.

    Article  PubMed  Google Scholar 

  15. Li H, Bai R, Zhao Z, Tao L, Ma M, Ji Z, et al. Application of droplet digital PCR to detect the pathogens of infectious diseases. Biosci Rep. 2018;38(6).

  16. Chen B, Xie Y, Zhang N, Li W, Liu C, Li D, et al. Evaluation of Droplet Digital PCR assay for the diagnosis of Candidemia in blood samples. Front Microbiol. 2021;12:700008.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Zheng Y, Jin J, Shao Z, Liu J, Zhang R, Sun R, et al. Development and clinical validation of a droplet digital PCR assay for detecting Acinetobacter baumannii and Klebsiella pneumoniae in patients with suspected bloodstream infections. Microbiologyopen. 2021;10(6):e1247.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Jitmuang A, Nititammaluk A, Boonsong T, Sarasombath PT, Sompradeekul S, Chayakulkeeree M. A novel droplet digital polymerase chain reaction for diagnosis of Pneumocystis pneumonia (PCP)-a clinical performance study and survey of sulfamethoxazole-trimethoprim resistant mutations. J Infect. 2021;83(6):701–8.

    Article  CAS  PubMed  Google Scholar 

  19. Wu Y, Wang F, Wang C, Tang X, Liu X, Li S, et al. Detection of Pneumocystis jirovecii and Toxoplasma gondii in patients with lung infections by a duplex qPCR assay. PLoS Negl Trop Dis. 2021;15(12):e0010025.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Liu L, Yuan M, Shi Y, Su X. Clinical performance of BAL Metagenomic Next-Generation sequence and serum (1,3)-β-D-Glucan for Differential diagnosis of Pneumocystis jirovecii Pneumonia and Pneumocystis jirovecii Colonisation. Front Cell Infect Microbiol. 2021;11:784236.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Bordon J, Aliberti S, Fernandez-Botran R, Uriarte SM, Rane MJ, Duvvuri P, et al. Understanding the roles of cytokines and neutrophil activity and neutrophil apoptosis in the protective versus deleterious inflammatory response in pneumonia. Int J Infect Dis. 2013;17(2):e76–83.

    Article  CAS  PubMed  Google Scholar 

  22. Giulia B, Luisa A, Concetta S, Bruna LS, Chiara B, Marcello C. Procalcitonin and community-acquired pneumonia (CAP) in children. Clin Chim Acta. 2015;451(Pt B):215–8.

    Article  CAS  PubMed  Google Scholar 

  23. Huang X, Li D, Liu F, Zhao D, Zhu Y, Tang H. Clinical significance of D-dimer levels in refractory Mycoplasma pneumoniae pneumonia. BMC Infect Dis. 2021;21(1):14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Zou XL, Feng DY, Wu WB, Yang HL, Zhang TT. Blood urea nitrogen to serum albumin ratio independently predicts 30-day mortality and severity in patients with Escherichia coli bacteraemia. Med Clin (Barc). 2021;157(5):219–25.

    Article  CAS  PubMed  Google Scholar 

  25. Watanabe T, Yasuoka A, Tanuma J, Yazaki H, Honda H, Tsukada K, et al. Serum (1–>3) beta-D-glucan as a noninvasive adjunct marker for the diagnosis of Pneumocystis pneumonia in patients with AIDS. Clin Infect Dis. 2009;49(7):1128–31.

    Article  CAS  PubMed  Google Scholar 

  26. Mercier T, Aissaoui N, Gits-Muselli M, Hamane S, Prattes J, Kessler HH, et al. Variable correlation between Bronchoalveolar Lavage Fluid Fungal load and Serum-(1,3)-β-d-Glucan in patients with Pneumocystosis-A Multicenter ECMM Excellence Center Study. J Fungi (Basel). 2020;6(4).

  27. Szvalb AD, Malek AE, Jiang Y, Bhatti MM, Wurster S, Kontoyiannis DP. Serum (1,3)-Beta-d-Glucan has suboptimal performance for the diagnosis of Pneumocystis jirovecii pneumonia in cancer patients and correlates poorly with respiratory burden as measured by quantitative PCR. J Infect. 2020;81(3):443–51.

    Article  CAS  PubMed  Google Scholar 

  28. Tang G, Tong S, Yuan X, Lin Q, Luo Y, Song H, et al. Using routine laboratory markers and immunological indicators for Predicting Pneumocystis jiroveci Pneumonia in Immunocompromised Patients. Front Immunol. 2021;12:652383.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Ewig S, Bauer T, Schneider C, Pickenhain A, Pizzulli L, Loos U, et al. Clinical characteristics and outcome of Pneumocystis carinii pneumonia in HIV-infected and otherwise immunosuppressed patients. Eur Respir J. 1995;8(9):1548–53.

    Article  CAS  PubMed  Google Scholar 

  30. Sun J, Su J, Xie Y, Yin MT, Huang Y, Xu L, et al. Plasma IL-6/IL-10 ratio and IL-8, LDH, and HBDH Level predict the severity and the risk of death in AIDS patients with Pneumocystis Pneumonia. J Immunol Res. 2016;2016:1583951.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Torres HA, Chemaly RF, Storey R, Aguilera EA, Nogueras GM, Safdar A, et al. Influence of type of cancer and hematopoietic stem cell transplantation on clinical presentation of Pneumocystis jiroveci pneumonia in cancer patients. Eur J Clin Microbiol Infect Dis. 2006;25(6):382–8.

    Article  CAS  PubMed  Google Scholar 

Download references


The authors would like to thank TargetingOne Corporation (Beijing, China) for the supply of the reagents to undertake this study.


This research was supported by Evaluation of droplet digital PCR for the detection of Pneumocystis jirovecii in alveolar lavage fluid and analysis of correlation between copy number and clinical characteristics and the National and Provincial Key Clinical Specialty Capacity Building Project 2020 (Department of the Respiratory Medicine).

Author information

Authors and Affiliations



Conceptualization, Gao ZC and Gong PH; methodology, Bian WJ and Xie Y; software, Yang ZW; validation, Shang Y; formal analysis, Bian WJ; investigation, Xie Y; resources, Zhao LL; data curation, Ma XQ; writing—original draft preparation, Bian WJ; writing—review and editing, Xie Y; visualization, He YK; supervision, Yu WY and Chen YW; project administration, Xi W, Yang DH, Wang F; funding acquisition, Gao ZC and Gong PH. All authors have read and agreed to the published version of the manuscript.All authors reviewed the manuscript.

Corresponding authors

Correspondence to Pihua Gong or Zhancheng Gao.

Ethics declarations

Competing interests

The authors declare no competing interests.

Ethics approval and consent to participate

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board of Peking University People’s Hospital (No. 2022PHB457). Informed consent was obtained from all subjects involved in the study.

Consent for publication

Not applicable.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Bian, W., Xie, Y., Shang, Y. et al. Relationship between clinical features and droplet digital PCR copy number in non-HIV patients with pneumocystis pneumonia. BMC Infect Dis 23, 833 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: