Skip to main content

ORIGINAL RESEARCH article

Front. Immunol., 15 November 2022
Sec. Microbial Immunology

Circulating mitochondrial cell-free DNA dynamics in patients with mycobacterial pulmonary infections: Potential for a novel biomarker of disease

Sheng-Wei Pan,,Sheng-Wei Pan1,2,3Rehan R. SyedRehan R. Syed4Donald G. CatanzaroDonald G. Catanzaro5Mei-Lin Ho,Mei-Lin Ho6,7Chin-Chung Shu,&#x;Chin-Chung Shu8,9†Tsung-Yeh Tsai,Tsung-Yeh Tsai1,2Yen-Han Tseng,Yen-Han Tseng1,2Jia-Yih Feng,*&#x;Jia-Yih Feng1,2*†Yuh-Min Chen,Yuh-Min Chen1,2Wei-Juin SuWei-Juin Su10Antonino CatanzaroAntonino Catanzaro3Timothy C. RodwellTimothy C. Rodwell3
  • 1Department of Chest Medicine, Taipei Veterans General Hospital, Taipei, Taiwan
  • 2School of Medicine, National Yang Ming Chiao Tung University, Taipei, Taiwan
  • 3Division of Pulmonary, Critical Care and Sleep Medicine, Department of Medicine, University of California San Diego, La Jolla, CA, United States
  • 4Division of Infectious Diseases and Global Public Health, University of California San Diego, La Jolla, CA, United States
  • 5Department of Biological Sciences, University of Arkansas, Fayetteville, AR, United States
  • 6Department of Chemistry, Soochow University, Taipei, Taiwan
  • 7Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA, United States
  • 8Department of Internal Medicine, National Taiwan University Hospital, Taipei, Taiwan
  • 9College of Medicine, National Taiwan University, Taipei, Taiwan
  • 10Division of Chest Medicine, China Medical University Hospital, Taipei Branch, Taipei, Taiwan

Objectives: Human mitochondrial cell-free DNA (Mt-cfDNA) may serve as a useful biomarker for infectious processes. We investigated Mt-cfDNA dynamics in patients with pulmonary mycobacterial infections to determine if this novel biomarker could be used to differentiate disease states and severity.

Methods: Patients with pulmonary tuberculosis (PTB), latent tuberculosis infection (LTBI), and nontuberculous mycobacterial-lung disease (NTM-LD) were enrolled at a tertiary care hospital in Taiwan between June 2018 and August 2021. Human Mt-cfDNA and nuclear-cfDNA (Nu-cfDNA) copy numbers were estimated by quantitative polymerase chain reaction. Variables associated with PTB and 2-month sputum culture-positivity, indicating poor treatment response, were assessed using logistic regression.

Results: Among 97 patients with PTB, 64 with LTBI, and 51 with NTM-LD, Mt-cfDNA levels were higher in patients with PTB than in LTBI (p=0.001) or NTM-LD (p=0.006). In the Mycobacterium tuberculosis-infected population, Mt-cfDNA levels were highest in smear-positive PTB patients, followed by smear-negative PTB (p<0.001), and were lowest in LTBI persons (p=0.009). A Mt-cfDNA, but not Nu-cfDNA, level higher than the median helped differentiate culture-positive PTB from culture-negative PTB and LTBI (adjusted OR 2.430 [95% CI 1.139–5.186], p=0.022) and differentiate PTB from NTM-LD (adjusted OR 4.007 [1.382–12.031], p=0.011). Mt-cfDNA levels decreased after 2 months of treatment in PTB patients (p=0.010). A cutoff Mt-cfDNA level greater than 62.62 x 106 copies/μL-plasma was associated with a 10-fold risk of 2-month culture-positivity (adjusted OR 9.691 [1.046–89.813], p=0.046).

Conclusion: Elevated Mt-cfDNA levels were associated with PTB disease and failed sputum conversion at 2 months in PTB patients, and decreased after treatment.

Introduction

An estimated 10 million people infected with Mycobacterium tuberculosis (MTB) develop active pulmonary tuberculosis disease (PTB) yearly (1), and up to one third of the population has latent tuberculosis infection (LTBI) at any time, with a 5–10% lifelong risk of developing PTB (2). LTBI, subclinical TB and active TB disease comprise the spectrum of MTB infection. The ability to identify individuals with PTB, differentiate LTBI from PTB and other mycobacterial infections, and to verify treatment success, is key both for individual patient outcomes and for preventing MTB transmission to others. Sputum culture positivity for MTB is one of the primary reference standards for TB diagnosis, just as the conversion from sputum positive to negative after 2 months of intensive-phase treatment is considered a reference indicator of therapeutic response. However, the slow growth of MTB cultures contributes to significant diagnostic delays, making evaluating sputum culture conversion both time-consuming and laborious. Easy-to-measure host inflammatory biomarkers that can differentiate PTB from other diagnoses or that can be used to support assessments of TB treatment response, would be transformative (3). Levels of c-reactive protein (CRP) and blood monocyte-to-lymphocyte ratio (MLR) may help to identify PTB in high-risk populations, and both decrease during treatment, but their changes are poorly predictive of sputum culture conversion at 2 months (38). Interferon-γ release assays (IGRA) measure interferon-γ response to MTB-specific antigens and are useful for identifying individuals with LTBI, however, they cannot be used to differentiate LTBI from PTB reliably, and are not useful for monitoring response to treatment (9, 10).

Cell-free DNA (cfDNA) is circulating fragmented DNA in blood that is actively secreted by living cells or released from dying cells and destroyed microorganisms and other pathogens (11, 12). Interest in cfDNA as a biomarker of altered physiologic and pathogenic states has prompted research into its diagnostic value in pregnancy and organ transplant monitoring and cancer diagnosis (1315). Elevated levels of host-derived cfDNA can be also be predictive of a poor prognosis in infection (1618). Both nuclear cfDNA (Nu-cfDNA) and mitochondrial cfDNA (Mt-cfDNA) are host-derived cfDNA, with the former originated from the nucleus and the later from mitochondria distributed throughout the cytoplasm in a cell. They can be released by activated inflammatory cells and only Mt-cfDNA can be released by platelets during infection (19, 20). Studies have shown that both plasma Nu-cfDNA and Mt-cfDNA levels correlate with disease severity and mortality in patients with sepsis (18, 21, 22). Prior research has also demonstrated that among HIV-positive individuals, plasma Mt-cfDNA levels were higher in those with TB and severe inflammation than in those without disease (23). Wiens et al. reported that macrophages infected with more virulent MTB released more mitochondrial DNA into the cytosol, suggesting mitochondrial DNA could be a biomarker indicating innate immune responses in TB (24). Studies have also demonstrated that MTB can cause mitochondrial damage in macrophages, leading to the release of mitochondrial DNA (25). Yet, it remains unclear if the release of Mt-cfDNA from inflammatory cells into the extracellular space, presumably resulting in high circulating Mt-cfDNA level, is related to burden of MTB nor if observed levels are different for LTBI or PTB and other pulmonary mycobacterial diseases. It is also not yet known if elevated Mt-cfDNA levels at PTB diagnosis reverts to baseline after PTB disease is effectively treated.

We hypothesized that plasma Mt-cfDNA is associated with MTB bacterial burden and inflammatory response in TB. In order to explore these ideas we a) examined the association between plasma Mt-cfDNA levels in patients with and without culture-positive PTB; b) evaluated changes in plasma Mt-cfDNA level during PTB treatment and c) evaluated the potential predictive value of Mt-cfDNA levels on subsequent failure of sputum conversion after two months of treatment. Finally, we measured levels of Mt-cfDNA in supernatants from a cell line stimulated with MTB early secretory antigenic target (ESAT)-6 antigen to develop potential supporting evidence for the causative role of MTB on Mt-cfDNA release from macrophages.

Methods

Study design and enrollment

This prospective observational study was conducted at the Taipei Veterans General Hospital in Taipei City, Taiwan, from June 2018 to August 2021. The cohort was composed of patients with LTBI, PTB, and nontuberculous mycobacterial-lung disease (NTM-LD). Adult outpatients with PTB were enrolled if 1) their respiratory samples were culture-positive for MTB or 2) they were culture-negative but a lung biopsy showed compatible pathologic findings and tested positive for the MTB-specific IS6110 gene (6). Asymptomatic persons with a history of close contact in the past six months with another PTB patient (prior to enrollment) were included if they had normal chest radiography and IGRA-positivity, indicating LTBI (26). As comparators, patients with NTM-LD were enrolled if they had positive NTM cultures from respiratory samples and met the NTM-LD diagnostic criteria for active disease (27, 28). The exclusion criteria were 1) HIV infection in PTB patients tested at enrollment and in NTM-LD patients assessed by reviewing claim data, 2) concomitant active cancer, and 3) mixed MTB and NTM-LD infection. Written informed consents were obtained from participants (Taipei Veterans General Hospital Institutional Review Board Nos. 2017-12-001C, 2018-10-017A, 2019-07-003C, 2020-07-009AC, 2021-01-010BC).

Blood sampling for cfDNA test

Peripheral blood samples were obtained from all participants at baseline before treatment. In patients with PTB, blood tests were performed again after two months of standard intensive-phase anti-TB therapy. For plasma preparation, 8–10 mL of whole blood was collected with K2-EDTA Vacutainer tubes (BD, Franklin Lakes, NJ, USA) and centrifuged for 10 minutes at 1500 rpm within 2 hours of collection (29). Using the QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany), cfDNA was extracted from a 400 µL plasma sample with an elution volume of 50 µL and stored in a microcentrifuge tube at −80°C until real-time quantitative polymerase chain reaction (qPCR) was performed in batches (30, 31).

Real-time qPCR with absolute quantification

To estimate the Mt-cfDNA and Nu-cfDNA copy numbers in plasma samples, the levels of a noncoding human mitochondrial genome region (hMito) and a nuclear DNA, human beta2 microglobulin (hB2M), were measured by qPCR (20). To calculate the absolute copy number of the target genes from the cycle threshold (Ct) value, a reference 5-point standard curve was constructed in duplicate on each plate with known copies of hMito and hB2M amplicons as previously reported (20). The levels of hMito and hB2M in plasma were expressed as copy numbers per μL of the original plasma sample (copies/μL plasma) after adjustment for the dilution factor of 1.25 (see Supplement Appendix 1 for details).

Clinical data collection and the outcome

Patient clinical data were collected at baseline for patients with PTB and NTM-LD and at 2 months among those who underwent anti-TB treatment. Data collected included body mass index (BMI), smoking status, medical comorbidities, radiographic score and findings (32), and sputum acid-fast bacilli cultures results and smear grades. Laboratory data including CRP and complete blood count with lymphocytes, monocytes, and platelet counts were also recorded. Persons with LTBI and patients with NTM-LD were confirmed to have no microbiologic evidence of PTB during a follow-up period of at least 3 months. Successful sputum culture conversion from positive to negative after two months of treatment was the outcome measure in patients with culture-positive PTB.

Measuring cfDNA in macrophages stimulated by MTB antigens

To try and establish a more direct link between MTB infection and the release of Mt-cfDNA and Nu-cfDNA from macrophages, we also examined the supernatant of monocytes stimulated with MTB antigens. Using a human leukemic cell line (THP-1) (TIB-202, ATCC, Manassas, VA, USA), THP-1 cells were cultured and treated with phorbol 12-myristate 13-acetate (PMA, 20 ng/ml) (P8139, Sigma-Aldrich, St. Louis, MO, USA) for 24 hours, followed by PMA removal and rest for 72 hours (33, 34). Differentiated cells were then stimulated by MTB-specific recombinant ESAT-6 protein (Recombinant Protein Reference Standard, NR-49424) (BEI Resources, Manassas, VA, USA) (34) to mimic a TB infection. After ESAT-6 protein was used to stimulate cells at concentrations of 0, 0.1, 1.0, and 2.0 µg/mL for 24 hours, culture supernatants were harvested for qPCR targeting Mt-cfDNA and Nu-cfDNA with absolute quantification.

Statistical analysis

The Mann–Whitney U-test and Student’s t-test were used to compare continuous variables between PTB, LTB and NTM-LD groups. A paired t-test was performed to compare the cfDNA levels before and after treatment for PTB. The distribution of a categorical variables between groups were compared by the χ2 test and Fisher’s exact test. Linear correlations between cfDNA levels and interested continuous variables were evaluated using Pearson correlation analysis. Logistic regression analysis was used to identify factors associated with culture-positive PTB and sputum conversion. Odds ratios (ORs) and 95% confidence intervals (CIs) were calculated using logistic regression. All variables with a p-value less than or equal to 0.05 in a univariate analysis were entered in a multivariable analysis of independent factors. Areas under the receiver operating characteristic (ROC) curves (AUCs) for variables of interest were calculated and if the lower boundary of a 95% CI for the AUC was above 0.5 the optimal cut-off value was determined (35, 36). Analyses were completed using GraphPad Prism version 6.0 (GraphPad Software Inc, CA, USA) and SPSS version 18.0 (SPSS Inc, IL, USA).

Results

Characteristics of participants

Sixty-four TB contacts with LTBI, 97 patients with PTB, and 51 with NTM-LD were enrolled over the 3-year study period (Figure 1). Among PTB patients, 73 were culture-positive while 24 were culture-negative. Of the culture-positive cases, 34 (46.6%) were smear-positive and 12 (16.4%) had isoniazid-resistant MTB isolates. No isolates were rifampicin resistant. The platelet count, monocyte count, MLR, and CRP levels in the PTB group were all higher than those in the LTBI group (all p<0.05) but similar to those in the NTM-LD group (Table 1). Although sputum smear grades were similar between the PTB and NTM-LD groups, the rates of hemoptysis, bilateral lung disease, and radiographic scores were higher in the NTM-LD group.

FIGURE 1
www.frontiersin.org

Figure 1 Flow chart of participant enrollment. C+, culture-positive; C-, culture-negative; S+, smear-positive; S-, smear-negative; LTBI, latent tuberculosis infection; NTM-LD, nontuberculous mycobacterial lung disease; PTB, pulmonary tuberculosis.

TABLE 1
www.frontiersin.org

Table 1 Demographic and clinical characteristics of patients with latent tuberculosis infection (LTBI), pulmonary tuberculosis (PTB), and nontuberculous mycobacterial lung disease (NTM-LD) (N=212).

Mt-cfDNA and Nu-cfDNA levels in groups

Median levels of Mt-cfDNA and Nu-cfDNA were 45.44 x106 copies/μL-plasma (interquartile range [IQR] 4.52 x106–157.73 x106) and 22.25 x104 copies/μL-plasma (IQR 2.17 x104–90.45 x104), respectively. Mt-cfDNA levels were higher in the PTB group (87.22 x106 copies/μL-plasma [IQR 12.20 x106–252.24 x106]) than in the LTBI group (20.63 x106 copies/μL-plasma [IQR 1.15 x106–84.68 x106]) (p=0.001, Figure 2A). Similarly, Nu-cfDNA levels were higher in the PTB group (37.01 x104 copies/μL-plasma [IQR 5.51 x104–126.44 x104]) than in the LTBI group (9.18 x104 copies/μL-plasma [IQR 0.57 x104–37.53 x104]) (p<0.001, Figure 2B). The median Mt-cfDNA level in the PTB group was also higher than that in the NTM-LD group (19.32 x106 copies/μL-plasma [IQR 0.89 x106–115.07 x106]) (p=0.006, Figure 2A) but the Nu-cfDNA was not (30.27 x104 copies/μL-plasma [IQR 2.22 x104–113.94 x104] in NTM-LD) (p=0.355, Figure 2B). Among patients infected with MTB, including PTB and LTBI but excluding NTM-LD, the level of Mt-cfDNA was higher in smear-positive PTB patients than in smear-negative PTB patients (p<0.001), but higher in smear-negative PTB patients than in persons with LTBI (p=0.009) (Figure 2C). This trend of association was not evident with Nu-cfDNA levels (Figure 2D). For patients with NTM-LD, Mt-cfDNA levels were not different between smear-positive and smear-negative cases (p=0.523). Overall, and even only among the PTB patients, Mt-cfDNA level was correlated with Nu-cfDNA, platelet count, monocyte count, and CRP (all p<0.01) and was not associated with age, gender, BMI, lymphocyte count, or MLR (see Supplement Table S1).

FIGURE 2
www.frontiersin.org

Figure 2 Human mitochondrial cell-free DNA (Mt-cfDNA) (A) and nuclear cell-free DNA (Nu-cfDNA) (B) levels in persons with latent tuberculosis infection (LTBI), pulmonary tuberculosis (PTB) and nontuberculous mycobacterial lung disease (NTM-LD). The levels of Mt-cfDNA (C) and Nu-cfDNA (D) in LTBI, smear-negative PTB (n=63, including 24 pathology-confirmed and 39 culture-positive cases) and smear-positive PTB subgroups. Data are presented as medians with interquartile ranges, and with each dot representing a value.

Factors associated with culture-positive PTB

In the MTB-infected cohort consisting of persons with LTBI and PTB but not NTM-LD (Figure 1), factors associated with culture-positive PTB in the multivariable analysis included BMI, platelet count, MLR, and Mt-cfDNA level higher than the median value of 45.44 x106 copies/μL-plasma (adjusted OR 2.430, 95% CI 1.139–5.186, p=0.022) (Table 2). In the mycobacterial culture-positive cohort consisting of culture-positive PTB and NTM-LD patients, the independent factors associated with culture-positive PTB were sex, platelet count, radiographic score, and Mt-cfDNA level higher than 45.44 x106 copies/μL-plasma (adjusted OR 4.007, 95% CI 1.382–12.031, p=0.011) (Table 2).

TABLE 2
www.frontiersin.org

Table 2 Logistic regression model of factors associated with the diagnosis of culture-positive pulmonary tuberculosis (PTB).

Sputum conversion and Mt-cfDNA changes in the PTB group

Twelve (16.4%) of the 73 initially culture-positive PTB patients were still culture-positive for MTB after 2 months of treatments. Two-month culture-positivity was associated with sputum smear grade and radiographic score at enrollment (Figure 3A). In ROC curve analyses, the AUC for smear grade and radiographic scores to discriminate persistent culture-positivity after 2 months of anti-TB therapy were 0.716 (95% CI 0.569–0.862) and 0.736 (95% CI 0.586–0.885), respectively, while that for Mt-cfDNA level was 0.658 (95% CI 0.507–0.810). The AUCs for Nu-cfDNA level, platelet count, monocyte count, MLR, and CPR were 0.634 (95% CI 0.480–0.788), 0.536 (95% CI 0.325–0.746), 0.663 (95% CI 0.488–0.839), 0.582 (95% CI 0.416–0.748), and 0.643 (95% CI 0.490–0.795), respectively, all with 95% CI lower bounds below 0.5. Using the optimal cutoff value of 62.62 x 106 copies/μL-plasma to dichotomize Mt-cfDNA level, the sensitivity, specificity, positive and negative predictive values of high Mt-cfDNA to predict two-month culture-positivity were 91.67%, 49.18%, 26.19% and 96.77%, respectively. An Mt-cfDNA level above this level was associated with a 10-fold risk of persistent culture-positivity (crude OR 10.645, 95% CI 1.293–87.607, p=0.028). The results were similar after adjusting for smear grade and radiographic score (adjusted OR 9.691, 95% CI 1.046–89.813, p=0.046) (Table 3).

FIGURE 3
www.frontiersin.org

Figure 3 The rates and numbers (n/N) of 2-month culture-positivity in PTB patients are shown stratified by initial sputum smear grade, initial radiographic score, and dynamic change of Mt-cfDNA levels (A).The levels of human mitochondrial cell-free DNA (Mt-cfDNA) (B) and nuclear cell-free DNA (Nu-cfDNA) (C) at baseline and at 2 months of therapy in patients with pulmonary tuberculosis (PTB). Here, Mt-cfDNA levels at 2 months in the 59 patients with PTB were normally distributed and a paired t test was used. Data are presented as medians with interquartile ranges, and with each dot representing a value. Low↓, low Mt-cfDNA level (<62.62 x 106 copies/μL-plasma) at baseline with a further decrease at 2 months; Low↑, low Mt-cfDNA level at baseline with a subsequent increase at 2 months; High↓, high Mt-cfDNA level (> 62.62 x 106 copies/μL-plasma) at baseline with a subsequent decrease at 2 months; High↑, high Mt-cfDNA level at baseline with a further increase at 2 months.

TABLE 3
www.frontiersin.org

Table 3 Logistic regression model of factors associated with persistent culture-positive at 2 months of treatment in patients with PTB (N=73).

Among the PTB patients with follow-up blood tests (n=59, including 48 with initial culture-positivity), the median Mt-cfDNA level was lower after 2 months of treatment when compared to baseline level (75.31 x106 copies/μL-plasma [IQR 25.53 x106–197.93 x106] versus 97.78 x106 copies/μL-plasma [IQR 15.72 x106–314.22 x106], p=0.010 using paired t-test) but the Nu-cfDNA level was not statistically different (Figure 3B, C). We also examined Mt-cfDNA changes in 48 PTB patients with initial culture-positivity and used Mt-cfDNA dynamics to stratify patients (Figure 3A). In those patients with initially low levels of Mt-cfDNA (less than 62.62 x 106 copies/μL-plasma), all were culture negative at two months regardless of subsequent changes in Mt-cfDNA levels, while in those with initially high levels that subsequently decreased (starting above 62.62 x 106 copies/μL-plasma), 16% (4/25) were culture positive at two months, and in those with initially high and subsequently increasing levels of Mt-cfDNA, 40% (2/5) were still culture positive at two months (p=0.023 for trend) (Table 4).

TABLE 4
www.frontiersin.org

Table 4 Characteristics of pulmonary tuberculosis (PTB) patients stratified by Mt-cfDNA dynamic change at 2 months of treatment (N=59).

Mt-cfDNA in supernatant from stimulated monocytes

We quantified levels of Mt-cfDNA and Nu-cfDNA in supernatants from THP-1 cells (primary monocyte differentiated macrophages) stimulated with MTB ESAT-6 antigens. We observed a significantly higher mean Mt-cfDNA level (3.00 x106 ± 0.87 x106 copies/μL-plasma) in the supernatant after stimulation with high-dose (2.0 µg/mL) ESAT-6 protein than in the supernatants with no MTB antigen stimulation (1.18 x106 ± 0.50 x106 copies/μL-plasma) (p=0.011). But there was no significant difference in Nu-cfDNA level between the MTB antigen stimulated and non-stimulated samples (Figure 4A, B).

FIGURE 4
www.frontiersin.org

Figure 4 Measurement of mitochondrial cell-free DNA (Mt-cfDNA) (A) and nuclear cell-free DNA (Nu-cfDNA) (B) in supernatants from phorbol 12-myristate 13-acetate (PMA)-treated human leukemic cell culture after stimulation by recombinant 6 kDa early secretory antigenic target (ESAT-6) protein (at concentrations of 0, 0.1, 1.0 and 2.0 µg/mL). Data are presented as the mean ± standard deviation, and n=4 at every concentration level.

Discussion

In our study of patients with MTB and NTM infections we found that Mt-cfDNA levels were elevated in patients with PTB as compared to those with LTBI and NTM-LD, despite the greater clinical severity of disease among the NTM-LD patients. Notably, plasma Mt-cfDNA was positively associated with inflammatory markers and smear-positivity (a proxy for bacterial burden) in PTB patients, indicating a potential association of this biomarker with disease severity. For patients with culture positive PTB, high and increasing levels of Mt-cfDNA also appeared to have some potential for predicting persistent culture-positivity after 2 months of treatment, but more importantly, low Mt-cfDNA levels at start of treatment was highly predictive of culture negativity after two months of treatment. Furthermore, while evaluating the plausibility of our proposed linkage between plasma Mt-cfDNA and MTB infection, we observed stimulation of THP-1 cells with MTB antigen ESAT-6 is associated with Mt-cfDNA release into the extracellular space of inflammatory cells, supporting the hypothesis for the origin of the increased plasma Mt-cfDNA levels in PTB.

Identifying individuals with infectious PTB and monitoring treatment response remain critical aspects of TB control worldwide. Of the currently available blood biomarkers for this purpose, CRP performs better than symptom-based methods as a screening tool for PTB in people living with HIV (4), but neither its baseline level nor change over time is predictive of MTB sputum conversion, which despite its limitations, is still considered the reference standard for treatment success (5). Elevated monocyte and platelet counts are correlated with inflammation in PTB and trend down after treatment, but they cannot reliably predict treatment outcomes (3739). MLR may have utility in identifying TB disease in HIV-infected children and for differentiating active TB from LTBI in the general population (68), but persistently high MLRs have not been found to associate 2-month culture-positivity (7, 8). The present study, while limited in size and scope, indicates that Mt-cfDNA levels are elevated in individuals with PTB compared to those with LTBI which reinforces the potential of Mt-cfDNA as a novel biomarker for differentiating PTB from LTBI, possibly in conjunction with existing assays such as the IGRAs. In PTB patients, Mt-cfDNA level was positively correlated with CRP, monocyte count, and platelet count, and was higher in smear-positive patients than in smear-negative patients, suggesting that Mt-cfDNA is correlated with both inflammatory processes and bacterial burden. Our findings are supported by prior research showing elevated Mt-cfDNA to be associated with TB-related immune reconstitution inflammatory syndrome and interleukin-18, a proinflammatory cytokine, in HIV-infected patients (23). From a clinical standpoint, when excluding MTB smear-positive patients who can more easily be diagnosed with PTB using existing diagnostics, Mt-cfDNA level was also statistically different between smear-negative PTB patients and the contacts with LTBI, suggesting a potential role for Mt-cfDNA for differentiating low-bacterial-burden PTB from IGRA positive LTBI patients. However, the presence of outliers with low Mt-cfDNA levels in the smear-negative PTB subgroup suggests that Mt-cfDNA alone might not have sufficient sensitivity in the clinical setting to differentiate all smear-negative PTB from LTBI. However, it may nonetheless play complementary role in the medical decision making model when evaluating these patients. Based on the totality of our results, we believe that human origin Mt-cfDNA has the potential to be used as an indicator of the clinical stage and severity of MTB infection, reflecting the clinical spectrum of MTB infection from LTBI to smear-negative and smear-positive PTB, and its use as a standalone or complementary biomarker in the evaluation of patients with suspected PTB and LTBI warrants further investigation.

Importantly, our study demonstrates that Mt-cfDNA may also have some utility in treatment monitoring for patients with PTB. A high baseline Mt-cfDNA level was found to be moderately predictive of 2-month culture-positivity in PTB patients and it was observed that Mt-cfDNA levels decline after 2 months of treatment. Our observed rate of persistent culture-positivity of 16.4% at 2 months is in line with previous reports (35, 40). Notably, using changes in Mt-cfDNA to stratify patients, we noted that the rate of persistent 2-month culture-positivity was highest in the PTB subgroup with initially high Mt-cfDNA and subsequently increasing levels at 2 months while the rates were lower in other subgroups, with a statistically significant trend. We hypothesize that Mt-cfDNA could be an integrated surrogate marker for bacterial load and inflammation in PTB, which are important factors associated with decreased rates of on-treatment culture-conversion (3, 41, 42). Similarly, decreased Mt-cfDNA levels reflect decreased bacterial load and reduced inflammation following treatment, which could have clinical relevance in indicating successful treatment. To this end, while the exact cause of increased circulating Mt-cfDNA levels is likely complex and multidimensional, our proof-of-concept study suggests the potential value of Mt-cfDNA monitoring the treatment response of PTB.

Human cfDNA originates from cell apoptosis, neutrophil extracellular traps, and the death of inflammatory cells in response to microbes, as well as from active cellular secretion through extracellular vehicles (12). In comparison to Nu-cfDNA, Mt-cfDNA is a thousand-fold higher in copy numbers because every cell, except red blood cells, contains abundant mitochondria with DNA. We found Mt-cfDNA levels better correlate with culture-positive PTB and predict 2-month culture-positivity than Nu-cfDNA. These findings may be mechanistically supported by the relationship we observed between ESAT-6 stimulation and Mt-cfDNA levels in supernatants from activated THP-1 cells. However, human Mt-cfDNA may be elevated by various physical or psychiatric stresses and pathologic insults (17, 43, 44). To investigate the specificity of Mt-cfDNA as a biomarker of PTB, we also included patients with NTM-LD in the study and found that high Mt-cfDNA facilitated the differentiation of culture-positive PTB from NTM-LD. Interestingly, although the NTM-LD group had higher rates of hemoptysis and higher radiographic scores indicating more severe disease among the NTM patients, elevated Mt-cfDNA levels still differentiated PTB from NTM-LD suggesting the possibility of a unique association between PTB disease and Mt-cfDNA level that goes beyond lung pathology. Using MTB strains with different virulence to infect bone marrow-derived macrophages, Wiens et al. reported the strain-dependent infection-related release of mitochondrial DNA into the cytosol (24). Because the levels of several inflammatory markers and bacterial load were similar in the PTB and NTM-LD groups, the differentially higher level of Mt-cfDNA in PTB patients may indicate a specific role for MTB in causing the release of Mt-cfDNA that deserves further investigation.

The present study has some limitations. First, the timing of laboratory testing during the intensive phase of treatment in PTB patients was limited to only two time points and the number of cases eligible for persistence culture-positivity assessment was relatively small. Large-scale studies are warranted to validate the association between dynamic changes in Mt-cfDNA during follow-up and bacterial burden during and after tuberculosis treatment. Second, the enrolled participants were clinically stable HIV-uninfected outpatients at a single institute, so the generalizability of our findings to other populations of different conditions on diverse treatment regimens remains unknown. Further studies on the performance of Mt-cfDNA in differentiating mycobacterial diagnosis and predicting outcomes among PTB and NTM-LD patients with HIV co-infection are needed. In addition, confounders such as diabetes severity and antihyperglycemic medications might affect the treatment outcome in PTB patients but were not available for analysis (45). Our study is also limited by the lack of a healthy control group. In one recently published study, the authors disclosed that patients with multidrug-resistant TB had increased mitochondrial DNA copy numbers in blood leukocytes when compared to healthy individuals, and proposed this finding was a result of oxidative stress in TB patients (46), which may be one of the plausible mechanisms explaining the elevated extracellular plasma Mt-cfDNA levels in our patients with PTB. Further investigations incorporating healthy controls may help confirm the potential role of Mt-cfDNA in differentiating the spectrum ranging from an uninfected population, to LTBI and PTB. Finally, the plasma samples were not double centrifuged and, thus, the level of cfDNA may have been overestimated (31). However, our results should be robust because, using the same protocol, all intergroup and post-treatment comparisons consistently demonstrated that high levels of Mt-cfDNA were associated with PTB-related inflammation and bacterial burden. Nevertheless, the potential impact of differences in cfDNA measurement protocols among different studies should be acknowledged.

In conclusion, we found that Mt-cfDNA level was correlated with inflammatory markers and smear-positivity in PTB. Importantly, Mt-cfDNA, but not Nu-cfDNA, was associated with culture-positive PTB, was associated with sputum culture-positivity at 2 months of therapy, and decreased after treatment. Further studies are warranted to evaluate the potential complementary role of Mt-DNA in PTB identification and treatment response monitoring.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation. Requests to access these datasets should be directed to Sheng-Wei Pan, sanweipan@gmail.com.

Ethics statement

The studies involving human participants were reviewed and approved by Taipei Veterans General Hospital Institutional Review Board Nos. 2017-12-001C, 2018-10-017A, 2019-07-003C, 2020-07-009AC, 2021-01-010BC. The patients/participants provided their written informed consent to participate in this study.

Author contributions

S-WP, T-YT, J-YF, and W-JS collected data. All authors performed data analysis. S-WP, RS, DC, M-LH, and C-CS wrote the manuscript. J-YF and TR provided critical revisions for the paper. S-WP and J-YF were the guarantors of the paper, taking responsibility for the integrity of the work as a whole, from inception to published article. All authors contributed to the article and approved the submitted version.

Funding

S-WP was funded by grants from Taipei Veterans General Hospital [V110C-040] and the Ministry of Science and Technology [107-2314-B-075-057, 108-2314-B-075-001, 109-2314-B-075-094, 110-2314-B-075-077, 111-2314-B-075-059] in Taiwan. TR, DC and AC were funded by NIH NIAID grant # R01AI137681.

Acknowledgments

The authors thank the Medical Science and Technology Building of the Taipei Veterans General Hospital and the Stein clinical research building of the University of California, San Diego, for providing research facilities. The authors also would like to thank Dr. Peter G Chiles for providing help in the experiment consultation.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.1040947/full#supplementary-material

References

1. World Health Organization. Global tuberculosis report (2021). Available at: https://www.who.int/teams/global-tuberculosis-programme/tb-reports.

Google Scholar

2. Comstock GW, Livesay VT, W. SF. The prognosis of a positive tuberculin reaction in childhood and adolescence. Am J Epidemiol (1974) 99:131–8. doi: 10.1093/oxfordjournals.aje.a121593

CrossRef Full Text | Google Scholar

3. Zimmer AJ, Lainati F, Aguilera Vasquez N, Chedid C, McGrath S, Benedetti A, et al. Biomarkers that correlate with active pulmonary tuberculosis treatment response: A systematic review and meta-analysis. J Clin Microbiol (2022) 60:e0185921. doi: 10.1128/jcm.01859-21

CrossRef Full Text | Google Scholar

4. Yoon C, Semitala FC, Atuhuzuma E, Katende J, Mwebe S, Asege L, et al. Point-of-care c-reactive protein-based tuberculosis screening for people living with HIV: a diagnostic accuracy study. Lancet Infect Dis (2017) 17:1285–92. doi: 10.1016/S1473-3099(17)30488-7

CrossRef Full Text | Google Scholar

5. Visser ME, Stead MC, Walzl G, Warren R, Schomaker M, Grewal HM, et al. Baseline predictors of sputum culture conversion in pulmonary tuberculosis: importance of cavities, smoking, time to detection and W-Beijing genotype. PloS One (2012) 7:e29588. doi: 10.1371/journal.pone.0029588

CrossRef Full Text | Google Scholar

6. Pan SW, Su WJ, Chan YJ, Chuang FY, Feng JY, C. YM. Mycobacterium tuberculosis-derived circulating cell-free DNA in patients with pulmonary tuberculosis and persons with latent tuberculosis infection. PloS One (2021) 16:0253879. doi: 10.1371/journal.pone.0253879

CrossRef Full Text | Google Scholar

7. Buttle TS, Hummerstone CY, Billahalli T, Ward RJB, Barnes KE, Marshall NJ, et al. The monocyte-to-lymphocyte ratio: Sex-specific differences in the tuberculosis disease spectrum, diagnostic indices and defining normal ranges. PloS One (2021) 18:e0247745. doi: 10.1371/journal.pone.0247745

CrossRef Full Text | Google Scholar

8. Choudhary RK, Wall KM, Njuguna I, Pavlinac PB, LaCourse SM, Otieno V, et al. Monocyte-to-Lymphocyte ratio is associated with tuberculosis disease and declines with anti-TB treatment in HIV-infected children. J Acquir Immune Defic Syndr (2019) 80:174–81. doi: 10.1097/QAI.0000000000001893

CrossRef Full Text | Google Scholar

9. Carranza C, Pedraza-Sanchez S, de Oyarzabal-Mendez E. Diagnosis for latent tuberculosis infection: New alternatives. Front Immunol (2020) 11:2006. doi: 10.3389/fimmu.2020.02006

CrossRef Full Text | Google Scholar

10. Johnson JL, Geldenhuys H TA, Thiel BA, Suliman S, Pienaar B, Chheng P, et al. Effect of isoniazid therapy for latent TB infection on QuantiFERON-TB gold in-tube responses in adults with positive tuberculin skin test results in a high TB incidence area: a controlled study. Chest (2014) 145:612–7. doi: 10.1378/chest.13-1232

CrossRef Full Text | Google Scholar

11. Fernández-Carballo BL, Broger T, Wyss R, Banaei N. Toward the development of a circulating free DNA-based In vitro diagnostic test for infectious diseases: a review of evidence for tuberculosis. J Clin Microbiol (2019) 57:e01234–18. doi: 10.1128/JCM.01234-18

CrossRef Full Text | Google Scholar

12. Duvvuri B. Cell-free DNA as a biomarker in autoimmune rheumatic diseases. Front Immunol (2019) 10:502. doi: 10.3389/fimmu.2019.00502

CrossRef Full Text | Google Scholar

13. Ordoñez E, Rueda L, Cañadas MP, Fuster C. Evaluation of sample stability and automated DNA extraction for fetal sex determination using cell-free fetal DNA in maternal plasma. BioMed Res Int (2013) 2013:195363. doi: 10.1155/2013/195363

CrossRef Full Text | Google Scholar

14. De Vlaminck I, Valantine HA, Snyder TM, Strehl C, Cohen G, Luikart H, et al. Circulating cell-free DNA enables noninvasive diagnosis of heart transplant rejection. Sci Transl Med (2014) 6:241ra77. doi: 10.1126/scitranslmed.3007803

CrossRef Full Text | Google Scholar

15. Rolfo C, Mack PC, Scagliotti GV, Baas P, Barlesi F, Bivona TG, et al. Liquid biopsy for advanced non-small cell lung cancer (NSCLC): A statement paper from the IASLC. J Thorac Oncol (2018) 13:1248–68. doi: 10.1016/j.jtho.2018.05.030

CrossRef Full Text | Google Scholar

16. Zhang Y, Song R, Shen Y, Zhao Y, Zhao Z, Fan T, et al. High levels of circulating cell-free DNA are associated with a poor prognosis in patients with severe fever with thrombocytopenia syndrome. Clin Infect Dis (2020) 70:1941–9. doi: 10.1093/cid/ciz553

CrossRef Full Text | Google Scholar

17. Scozzi D, Cano M, Ma L, Zhou D, Zhu JH, O’Halloran JA, et al. Circulating mitochondrial DNA is an early indicator of severe illness and mortality from COVID-19. JCI Insight (2021) 6:e143299. doi: 10.1172/jci.insight.143299

CrossRef Full Text | Google Scholar

18. Kung CT, Hsiao SY, Tsai TC, Su CM, Chang WN, Huang CR, et al. Plasma nuclear and mitochondrial DNA levels as predictors of outcome in severe sepsis patients in the emergency room. J Transl Med (2012) 10:130. doi: 10.1186/1479-5876-10-130

CrossRef Full Text | Google Scholar

19. Boudreau LH, Duchez AC, Cloutier N, Soulet D, Martin N, Bollinger J, et al. Platelets release mitochondria serving as substrate for bactericidal group IIA-secreted phospholipase A2 to promote inflammation. Blood (2014) 124:2173–83. doi: 10.1182/blood-2014-05-573543

CrossRef Full Text | Google Scholar

20. Rosa HS, Ajaz S, Gnudi L, M. AN. A case for measuring both cellular and cell-free mitochondrial DNA as a disease biomarker in human blood. FASEB J (2020) 34:12278–88. doi: 10.1096/fj.202000959RR

CrossRef Full Text | Google Scholar

21. Nakahira K, Kyung SY, Rogers AJ, Gazourian L, Youn S, Massaro AF, et al. Circulating mitochondrial DNA in patients in the ICU as a marker of mortality: derivation and validation. PloS Med (2013) 10:e1001577. doi: 10.1371/journal.pmed.1001577

CrossRef Full Text | Google Scholar

22. Yan HP, Li M, Lu XL, Zhu YM, Ou-Yang WX, Xiao ZH, et al. Use of plasma mitochondrial DNA levels for determining disease severity and prognosis in pediatric sepsis: a case control study. BMC Pediatr (2018) 18:267. doi: 10.1186/s12887-018-1239-z

CrossRef Full Text | Google Scholar

23. Tan HY, Yong YK, Shankar EM, Paukovics G, Ellegård R, Larsson M, et al. Aberrant inflammasome activation characterizes tuberculosis-associated immune reconstitution inflammatory syndrome. J Immunol (2016) 196:4052–63. doi: 10.4049/jimmunol.1502203

CrossRef Full Text | Google Scholar

24. Wiens KE, JD E. The mechanism for type I interferon induction by mycobacterium tuberculosis is bacterial strain-dependent. PloS Pathog (2016) 12:e1005809. doi: 10.1371/journal.ppat.1005809

CrossRef Full Text | Google Scholar

25. Olson GS, Murray TA, Jahn AN, Mai D, Diercks AH, Gold ES. Type I interferon decreases macrophage energy metabolism during mycobacterial infection. Cell Rep (2021) 35:109195. doi: 10.1016/j.celrep.2021.109195

CrossRef Full Text | Google Scholar

26. Pan SW, Kou YR, Hu TM, Wu YC, Lee YC, Feng JY, et al. Assessment of latent tuberculosis infection in psychiatric inpatients: A survey after tuberculosis outbreaks. J Microbiol Immunol Infect (2016) 49:575–83. doi: 10.1016/j.jmii.2015.10.005

CrossRef Full Text | Google Scholar

27. Pan SW, Shu CC, Lee CC, Feng JY, Chan YJ, Chen YM, et al. Role of soluble T-cell immunoglobulin mucin domain-3 in differentiating nontuberculous mycobacterial lung disease from pulmonary colonization. Arch Bronconeumol (2022) 58:547–553. doi: 10.1016/j.arbres.2021.01.035

CrossRef Full Text | Google Scholar

28. Griffith DE AT, Brown-Elliott BA, Catanzaro A, Daley C, Gordin F, Holland SM, et al. An official ATS/IDSA statement: diagnosis, treatment, and prevention of nontuberculous mycobacterial diseases. Am J Respir Crit Care Med (2007) 175:367–416. doi: 10.1164/rccm.200604-571ST

CrossRef Full Text | Google Scholar

29. Ushio R, Yamamoto M, Nakashima K, Watanabe H, Nagai K, Shibata Y, et al. Digital PCR assay detection of circulating mycobacterium tuberculosis DNA in pulmonary tuberculosis patient plasma. Tuberculosis (Edinb) (2016) 99:47–53. doi: 10.1016/j.tube.2016.04.004

CrossRef Full Text | Google Scholar

30. Mojtabanezhad Shariatpanahi A, Rokni P, Shahabi E, Varshoee Tabrizi F, K. MA. Simple and cost-effective laboratory methods to evaluate and validate cell-free DNA isolation. BMC Res Notes (2018) 11:757. doi: 10.1186/s13104-018-3866-8

CrossRef Full Text | Google Scholar

31. Vrablicova Z, Tomova K, Tothova L, Babickova J, Gromova B, Konecna B, et al. Nuclear and mitochondrial circulating cell-free DNA is increased in patients with inflammatory bowel disease in clinical remission. Front Med (Lausanne) (2020) 7:593316. doi: 10.3389/fmed.2020.593316

CrossRef Full Text | Google Scholar

32. Snider GL, Doctor L, Demas TA, S. AR. Obstructive airway disease in patients with treated pulmonary tuberculosis. Am Rev Respir Dis (1971) 103:625–40. doi: 10.1164/arrd.1971.103.5.625

CrossRef Full Text | Google Scholar

33. Starr T, Bauler TJ, Malik-Kale P. The phorbol 12-myristate-13-acetate differentiation protocol is critical to the interaction of THP-1 macrophages with salmonella typhimurium. PloS One (2018) 13:e0193601. doi: 10.1371/journal.pone.0193601

CrossRef Full Text | Google Scholar

34. Pan SW, Shu CC, Huang JR, Lee CC, Tseng YH, Hung JJ, et al. PD-L1 expression in monocytes correlates with bacterial burden and treatment outcomes in active pulmonary tuberculosis. Int J Mol Sci (2022) 23:1619. doi: 10.3390/ijms23031619

CrossRef Full Text | Google Scholar

35. Lee MR, Tsai CJ, Wang WJ, Chuang TY, Yang CM, Chang LY, et al. Plasma biomarkers can predict treatment response in tuberculosis patients: A prospective observational study. Med (Baltimore) (2015) 94:e1628. doi: 10.1097/MD.0000000000001628

CrossRef Full Text | Google Scholar

36. Ravibabu K, Tapu B, Ravi PJ, Vinay KA. Assessment of diagnostic accuracy and optimal cut points of blood lead levels on serum proteins among workers exposed to Pb at a lead battery plant. Int J Med Biochem (2019) 2:81–7. doi: 10.14744/ijmb.2019.41736

CrossRef Full Text | Google Scholar

37. Unsal E, Aksaray S, Köksal D. Potential role of interleukin 6 in reactive thrombocytosis and acute phase response in pulmonary tuberculosis. Postgrad Med J (2005) 81:604–7. doi: 10.1136/pgmj.2004.030544

CrossRef Full Text | Google Scholar

38. Turken O, Kunter E, Sezer M, Solmazgul E, Cerrahoglu K, Bozkanat E, et al. Hemostatic changes in active pulmonary tuberculosis. Int J Tuberc Lung Dis (2002) 6:927–32.

Google Scholar

39. Ştefanescu S, Cocoş R, Turcu-Stiolica A, Mahler B, Meca AD, Giura AMC, et al. Evaluation of prognostic significance of hematological profiles after the intensive phase treatment in pulmonary tuberculosis patients from Romania. PloS One (2021) 16:e0249301.

Google Scholar

40. Lee HY, Chae KO, Lee CH, Choi SM, Lee J, Park YS, et al. Culture conversion rate at 2 months of treatment according to diagnostic methods among patients with culture-positive pulmonary tuberculosis. PloS One (2014) 9:e103768. doi: 10.1371/journal.pone.0103768

CrossRef Full Text | Google Scholar

41. Sabiiti W, Azam K, Farmer ECW, Kuchaka D, Mtafya B, Bowness R, et al. Tuberculosis bacillary load, an early marker of disease severity: the utility of tuberculosis molecular bacterial load assay. Thorax (2020) 75:606–8. doi: 10.1136/thoraxjnl-2019-214238

CrossRef Full Text | Google Scholar

42. Su WJ, Feng JY, Chiu YC, Huang SF, L. YC. Role of 2-month sputum smears in predicting culture conversion in pulmonary tuberculosis. Eur Respir J (2011) 37:376–83. doi: 10.1183/09031936.00007410

CrossRef Full Text | Google Scholar

43. Simmons JD, Lee YL, Mulekar S, Kuck JL, Brevard SB, Gonzalez RP, et al. Elevated levels of plasma mitochondrial DNA DAMPs are linked to clinical outcome in severely injured human subjects. Ann Surg (2013) 258:591–6. doi: 10.1097/SLA.0b013e3182a4ea46

CrossRef Full Text | Google Scholar

44. Trumpff C, Marsland AL, Basualto-Alarcón C, Martin JL, Carroll JE, Sturm G, et al. Acute psychological stress increases serum circulating cell-free mitochondrial DNA. Psychoneuroendocrinology (2019) 106:268–76. doi: 10.1016/j.psyneuen.2019.03.026

CrossRef Full Text | Google Scholar

45. Degner NR, Wang JY, Golub JE, K. PC. Metformin use reverses the increased mortality associated with diabetes mellitus during tuberculosis treatment. Clin Infect Dis (2018) 66:198–205. doi: 10.1093/cid/cix819

CrossRef Full Text | Google Scholar

46. Freimane L, Barkane L, Igumnova V, Kivrane A, Zole E, Ranka R. Telomere length and mitochondrial DNA copy number in multidrug-resistant tuberculosis. Tuberculosis (Edinburgh Scotland) (2021) 131:102144. doi: 10.1016/j.tube.2021.102144

CrossRef Full Text | Google Scholar

Keywords: cell-free DNA, mitochondria, pulmonary tuberculosis, treatment response monitoring, biomarker

Citation: Pan S-W, Syed RR, Catanzaro DG, Ho M-L, Shu C-C, Tsai T-Y, Tseng Y-H, Feng J-Y, Chen Y-M, Su W-J, Catanzaro A and Rodwell TC (2022) Circulating mitochondrial cell-free DNA dynamics in patients with mycobacterial pulmonary infections: Potential for a novel biomarker of disease. Front. Immunol. 13:1040947. doi: 10.3389/fimmu.2022.1040947

Received: 09 September 2022; Accepted: 31 October 2022;
Published: 15 November 2022.

Edited by:

Esaki M. Shankar, Central University of Tamil Nadu, India

Reviewed by:

Vignesh Ramachandran, University of Kuala Lumpur, Malaysia
Vladimir M Pisarev, Federal Research and Clinical Center of Intensive Care Medicine and Rehabilitation, Russia

Copyright © 2022 Pan, Syed, Catanzaro, Ho, Shu, Tsai, Tseng, Feng, Chen, Su, Catanzaro and Rodwell. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jia-Yih Feng, peterofeng@gmail.com

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.