Modeling and predicting individual variation in COVID-19 vaccine-elicited antibody response in the general population.

As we learned during the COVID-19 pandemic, vaccines are one of the most important tools in infectious disease control. To date, an unprecedentedly large volume of high-quality data on COVID-19 vaccinations have been accumulated. For preparedness in future pandemics beyond COVID-19, these valuable datasets should be analyzed to best shape an effective vaccination strategy. We are collecting longitudinal data from a community-based cohort in Fukushima, Japan, that consists of 2,407 individuals who underwent serum sampling two or three times after a two-dose vaccination with either BNT162b2 or mRNA-1273. Using the individually reconstructed time courses of the vaccine-elicited antibody response based on mathematical modeling, we first identified basic demographic and health information that contributed to the main features of the antibody dynamics, i.e., the peak, the duration, and the area under the curve. We showed that these three features of antibody dynamics were partially explained by underlying medical conditions, adverse reactions to vaccinations, and medications, consistent with the findings of previous studies. We then applied to these factors a recently proposed computational method to optimally fit an "antibody score", which resulted in an integer-based score that can be used as a basis for identifying individuals with higher or lower antibody titers from basic demographic and health information. The score can be easily calculated by individuals themselves or by medical practitioners. Although the sensitivity of this score is currently not very high, in the future, as more data become available, it has the potential to identify vulnerable populations and encourage them to get booster vaccinations. Our mathematical model can be extended to any kind of vaccination and therefore can form a basis for policy decisions regarding the distribution of booster vaccines to strengthen immunity in future pandemics.

Multiscale modeling of HBV infection integrating intra- and intercellular viral propagation to analyze extracellular viral markers.

Chronic infection with hepatitis B virus (HBV) is caused by the persistence of closed circular DNA (cccDNA) in the nucleus of infected hepatocytes. Despite available therapeutic anti-HBV agents, eliminating the cccDNA remains challenging. Thus, quantifying and understanding the dynamics of cccDNA are essential for developing effective treatment strategies and new drugs. However, such study requires repeated liver biopsy to measure the intrahepatic cccDNA, which is basically not accepted because liver biopsy is potentially morbid and not common during hepatitis B treatment. We here aimed to develop a noninvasive method for quantifying cccDNA in the liver using surrogate markers in peripheral blood. We constructed a multiscale mathematical model that explicitly incorporates both intracellular and intercellular HBV infection processes. The model, based on age-structured partial differential equations, integrates experimental data from in vitro and in vivo investigations. By applying this model, we roughly predicted the amount and dynamics of intrahepatic cccDNA within a certain range using specific viral markers in serum samples, including HBV DNA, HBsAg, HBeAg, and HBcrAg. Our study represents a significant step towards advancing the understanding of chronic HBV infection. The noninvasive quantification of cccDNA using our proposed method holds promise for improving clinical analyses and treatment strategies. By comprehensively describing the interactions of all components involved in HBV infection, our multiscale mathematical model provides a valuable framework for further research and the development of targeted interventions.

Isolation may select for earlier and higher peak viral load but shorter duration in SARS-CoV-2 evolution.

During the COVID-19 pandemic, human behavior change as a result of nonpharmaceutical interventions such as isolation may have induced directional selection for viral evolution. By combining previously published empirical clinical data analysis and multi-level mathematical modeling, we find that the SARS-CoV-2 variants selected for as the virus evolved from the pre-Alpha to the Delta variant had earlier and higher peak in viral load dynamics but a shorter duration of infection. Selection for increased transmissibility shapes the viral load dynamics, and the isolation measure is likely to be a driver of these evolutionary transitions. In addition, we show that a decreased incubation period and an increased proportion of asymptomatic infection are also positively selected for as SARS-CoV-2 mutated to adapt to human behavior (i.e., Omicron variants). The quantitative information and predictions we present here can guide future responses in the potential arms race between pandemic interventions and viral evolution.

Analysis of the risk and pre-emptive control of viral outbreaks accounting for within-host dynamics: SARS-CoV-2 as a case study.

In the era of living with COVID-19, the risk of localised SARS-CoV-2 outbreaks remains. Here, we develop a multiscale modelling framework for estimating the local outbreak risk for a viral disease (the probability that a major outbreak results from a single case introduced into the population), accounting for within-host viral dynamics. Compared to population-level models previously used to estimate outbreak risks, our approach enables more detailed analysis of how the risk can be mitigated through pre-emptive interventions such as antigen testing. Considering SARS-CoV-2 as a case study, we quantify the within-host dynamics using data from individuals with omicron variant infections. We demonstrate that regular antigen testing reduces, but may not eliminate, the outbreak risk, depending on characteristics of local transmission. In our baseline analysis, daily antigen testing reduces the outbreak risk by 45% compared to a scenario without antigen testing. Additionally, we show that accounting for heterogeneity in within-host dynamics between individuals affects outbreak risk estimates and assessments of the impact of antigen testing. Our results therefore highlight important factors to consider when using multiscale models to design pre-emptive interventions against SARS-CoV-2 and other viruses.

Multiscale modeling of HBV infection integrating intra- and intercellular viral propagation for analyzing extracellular viral markers.

Chronic infection of hepatitis B virus (HBV) is caused by the persistence of closed circular DNA (cccDNA) in the nucleus of infected hepatocytes. Despite available therapeutic anti-HBV agents, eliminating the cccDNA remains challenging. The quantifying and understanding dynamics of cccDNA are essential for developing effective treatment strategies and new drugs. However, it requires a liver biopsy to measure the intrahepatic cccDNA, which is basically not accepted because of the ethical aspect. We here aimed to develop a non-invasive method for quantifying cccDNA in the liver using surrogate markers present in peripheral blood. We constructed a multiscale mathematical model that explicitly incorporates both intracellular and intercellular HBV infection processes. The model, based on age-structured partial differential equations (PDEs), integrates experimental data from in vitro and in vivo investigations. By applying this model, we successfully predicted the amount and dynamics of intrahepatic cccDNA using specific viral markers in serum samples, including HBV DNA, HBsAg, HBeAg, and HBcrAg. Our study represents a significant step towards advancing the understanding of chronic HBV infection. The non-invasive quantification of cccDNA using our proposed methodology holds promise for improving clinical analyses and treatment strategies. By comprehensively describing the interactions of all components involved in HBV infection, our multiscale mathematical model provides a valuable framework for further research and the development of targeted interventions.

OptEF-BCI: An Optimization-Based Hybrid EEG and fNIRS-Brain Computer Interface.

Multimodal data fusion (electroencephalography (EEG) and functional near-infrared spectroscopy (fNIRS)) has been developed as an important neuroimaging research field in order to circumvent the inherent limitations of individual modalities by combining complementary information from other modalities. This study employed an optimization-based feature selection algorithm to systematically investigate the complementary nature of multimodal fused features. After preprocessing the acquired data of both modalities (i.e., EEG and fNIRS), the temporal statistical features were computed separately with a 10 s interval for each modality. The computed features were fused to create a training vector. A wrapper-based binary enhanced whale optimization algorithm (E-WOA) was used to select the optimal/efficient fused feature subset using the support-vector-machine-based cost function. An online dataset of 29 healthy individuals was used to evaluate the performance of the proposed methodology. The findings suggest that the proposed approach enhances the classification performance by evaluating the degree of complementarity between characteristics and selecting the most efficient fused subset. The binary E-WOA feature selection approach showed a high classification rate (94.22 ± 5.39%). The classification performance exhibited a 3.85% increase compared with the conventional whale optimization algorithm. The proposed hybrid classification framework outperformed both the individual modalities and traditional feature selection classification (p < 0.01). These findings indicate the potential efficacy of the proposed framework for several neuroclinical applications.

Increased flexibility of the SARS-CoV-2 RNA-binding site causes resistance to remdesivir.

Mutations continue to accumulate within the SARS-CoV-2 genome, and the ongoing epidemic has shown no signs of ending. It is critical to predict problematic mutations that may arise in clinical environments and assess their properties in advance to quickly implement countermeasures against future variant infections. In this study, we identified mutations resistant to remdesivir, which is widely administered to SARS-CoV-2-infected patients, and discuss the cause of resistance. First, we simultaneously constructed eight recombinant viruses carrying the mutations detected in in vitro serial passages of SARS-CoV-2 in the presence of remdesivir. We confirmed that all the mutant viruses didn't gain the virus production efficiency without remdesivir treatment. Time course analyses of cellular virus infections showed significantly higher infectious titers and infection rates in mutant viruses than wild type virus under treatment with remdesivir. Next, we developed a mathematical model in consideration of the changing dynamic of cells infected with mutant viruses with distinct propagation properties and defined that mutations detected in in vitro passages canceled the antiviral activities of remdesivir without raising virus production capacity. Finally, molecular dynamics simulations of the NSP12 protein of SARS-CoV-2 revealed that the molecular vibration around the RNA-binding site was increased by the introduction of mutations on NSP12. Taken together, we identified multiple mutations that affected the flexibility of the RNA binding site and decreased the antiviral activity of remdesivir. Our new insights will contribute to developing further antiviral measures against SARS-CoV-2 infection.

Potential Anti-Mpox Virus Activity of Atovaquone, Mefloquine, and Molnupiravir, and Their Potential Use as Treatments.

BACKGROUND: Mpox virus (MPXV) is a zoonotic orthopoxvirus and caused an outbreak in 2022. Although tecovirimat and brincidofovir are approved as anti-smallpox drugs, their effects in mpox patients have not been well documented. In this study, by a drug repurposing approach, we identified potential drug candidates for treating mpox and predicted their clinical impacts by mathematical modeling. METHODS: We screened 132 approved drugs using an MPXV infection cell system. We quantified antiviral activities of potential drug candidates by measuring intracellular viral DNA and analyzed the modes of action by time-of-addition assay and electron microscopic analysis. We further predicted the efficacy of drugs under clinical concentrations by mathematical simulation and examined combination treatment. RESULTS: Atovaquone, mefloquine, and molnupiravir exhibited anti-MPXV activity, with 50% inhibitory concentrations of 0.51-5.2 µM, which was more potent than cidofovir. Whereas mefloquine was suggested to inhibit viral entry, atovaquone and molnupiravir targeted postentry processes. Atovaquone was suggested to exert its activity through inhibiting dihydroorotate dehydrogenase. Combining atovaquone with tecovirimat enhanced the anti-MPXV effect of tecovirimat. Quantitative mathematical simulations predicted that atovaquone can promote viral clearance in patients by 7 days at clinically relevant drug concentrations. CONCLUSIONS: These data suggest that atovaquone would be a potential candidate for treating mpox.

Smoking enhances the expression of angiotensin-converting enzyme 2 involved in the efficiency of severe acute respiratory syndrome coronavirus 2 infection.

Smoking is one of the risk factors most closely related to the severity of coronavirus disease 2019 (COVID-19). However, the relationship between smoking history and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infectivity is unknown. In this study, we evaluated the ACE2 expression level in the lungs of current smokers, ex-smokers, and nonsmokers. The ACE2 expression level of ex-smokers who smoked cigarettes until recently (cessation period shorter than 6 months) was higher than that of nonsmokers and ex-smokers with a long history of nonsmoking (cessation period longer than 6 months). We also showed that the efficiency of SARS-CoV-2 infection was enhanced in a manner dependent on the angiotensin-converting enzyme 2 (ACE2) expression level. Using RNA-seq analysis on the lungs of smokers, we identified that the expression of inflammatory signaling genes was correlated with ACE2 expression. Notably, with increasing duration of smoking cessation among ex-smokers, not only ACE2 expression level but also the expression levels of inflammatory signaling genes decreased. These results indicated that smoking enhances the expression levels of ACE2 and inflammatory signaling genes. Our data suggest that the efficiency of SARS-CoV-2 infection is enhanced by smoking-mediated upregulation of ACE2 expression level.

Artificial neural network models: implementation of functional near-infrared spectroscopy-based spontaneous lie detection in an interactive scenario.

Deception is an inevitable occurrence in daily life. Various methods have been used to understand the mechanisms underlying brain deception. Moreover, numerous efforts have been undertaken to detect deception and truth-telling. Functional near-infrared spectroscopy (fNIRS) has great potential for neurological applications compared with other state-of-the-art methods. Therefore, an fNIRS-based spontaneous lie detection model was used in the present study. We interviewed 10 healthy subjects to identify deception using the fNIRS system. A card game frequently referred to as a bluff or cheat was introduced. This game was selected because its rules are ideal for testing our hypotheses. The optical probe of the fNIRS was placed on the subject's forehead, and we acquired optical density signals, which were then converted into oxy-hemoglobin and deoxy-hemoglobin signals using the Modified Beer-Lambert law. The oxy-hemoglobin signal was preprocessed to eliminate noise. In this study, we proposed three artificial neural networks inspired by deep learning models, including AlexNet, ResNet, and GoogleNet, to classify deception and truth-telling. The proposed models achieved accuracies of 88.5%, 88.0%, and 90.0%, respectively. These proposed models were compared with other classification models, including k-nearest neighbor, linear support vector machines (SVM), quadratic SVM, cubic SVM, simple decision trees, and complex decision trees. These comparisons showed that the proposed models performed better than the other state-of-the-art methods.

Fractional transit compartment model for describing drug delayed response to tumors using Mittag-Leffler distribution on age-structured PKPD model.

The response of a cell population is often delayed relative to drug injection, and individual cells in a population of cells have a specific age distribution. The application of transit compartment models (TCMs) is a common approach for describing this delay. In this paper, we propose a TCM in which damaged cells caused by a drug are given by a single fractional derivative equation. This model describes the delay as a single equation composed of fractional and ordinary derivatives, instead of a system of ODEs expressed in multiple compartments, applicable to the use of the PK concentration in the model. This model tunes the number of compartments in the existing model and expresses the delay in detail by estimating an appropriate fractional order. We perform model robustness, sensitivity analysis, and change of parameters based on the amount of data. Additionally, we resolve the difficulty in parameter estimation and model simulation using a semigroup property, consisting of a system with a mixture of fractional and ordinary derivatives. This model provides an alternative way to express the delays by estimating an appropriate fractional order without determining the pre-specified number of compartments.

Designing isolation guidelines for COVID-19 patients with rapid antigen tests.

Appropriate isolation guidelines for COVID-19 patients are warranted. Currently, isolating for fixed time is adopted in most countries. However, given the variability in viral dynamics between patients, some patients may no longer be infectious by the end of isolation, whereas others may still be infectious. Utilizing viral test results to determine isolation length would minimize both the risk of prematurely ending isolation of infectious patients and the unnecessary individual burden of redundant isolation of noninfectious patients. In this study, we develop a data-driven computational framework to compute the population-level risk and the burden of different isolation guidelines with rapid antigen tests (i.e., lateral flow tests). Here, we show that when the detection limit is higher than the infectiousness threshold values, additional consecutive negative results are needed to ascertain infectiousness status. Further, rapid antigen tests should be designed to have lower detection limits than infectiousness threshold values to minimize the length of prolonged isolation.

Estimation of timing of infection from longitudinal SARS-CoV-2 viral load data: mathematical modelling study.

BACKGROUND: Multiple waves of the COVID-19 epidemic have hit most countries by the end of 2021. Most of those waves are caused by emergence and importation of new variants. To prevent importation of new variants, combination of border control and contact tracing is essential. However, the timing of infection inferred by interview is influenced by recall bias and hinders the contact tracing process. METHODS: We propose a novel approach to infer the timing of infection, by employing a within-host model to capture viral load dynamics after the onset of symptoms. We applied this approach to ascertain secondary transmission which can trigger outbreaks. As a demonstration, the 12 initial reported cases in Singapore, which were considered as imported because of their recent travel history to Wuhan, were analyzed to assess whether they are truly imported. RESULTS: Our approach suggested that 6 cases were infected prior to the arrival in Singapore, whereas other 6 cases might have been secondary local infection. Three among the 6 potential secondary transmission cases revealed that they had contact history to previously confirmed cases. CONCLUSIONS: Contact trace combined with our approach using viral load data could be the key to mitigate the risk of importation of new variants by identifying cases as early as possible and inferring the timing of infection with high accuracy.

Different efficacies of neutralizing antibodies and antiviral drugs on SARS-CoV-2 Omicron subvariants, BA.1 and BA.2.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron subvariant BA.2 has spread in many countries, replacing the earlier Omicron subvariant BA.1 and other variants. Here, using a cell culture infection assay, we quantified the intrinsic sensitivity of BA.2 and BA.1 compared with other variants of concern, Alpha, Gamma, and Delta, to five approved-neutralizing antibodies and antiviral drugs. Our assay revealed the diverse sensitivities of these variants to antibodies, including the loss of response of both BA.1 and BA.2 to casirivimab and of BA.1 to imdevimab. In contrast, EIDD-1931 and nirmatrelvir showed a more conserved activities to these variants. The viral response profile combined with mathematical analysis estimated differences in antiviral effects among variants in the clinical concentrations. These analyses provide essential evidence that gives insight into variant emergence's impact on choosing optimal drug treatment.

First-principles study on structural, electronic, magnetic and thermodynamic properties of lithium ferrite LiFe5O8.

Lithium ferrite, LiFe5O8 (LFO), has attracted great attention for various applications, and there has been extensive experimental studies on its material properties and applications. However, no systematic theoretical study has yet been reported, so understanding of its material properties at the atomic scale is still required. In this work, we present a comprehensive investigation into the structural, electronic, magnetic and thermodynamic properties of LFO using first-principles calculations. We demonstrate that the ordered α-phase with ferrimagnetic spin configuration is energetically favourable among various crystalline phases with different magnetic configurations. By applying the DFT + U approach with U = 4 eV, we reproduce the lattice constant, band gap energy, and total magnetization in good agreement with experiments, emphasizing the importance of considering strong correlation and spin-polarization effects originating from the 3d states of Fe atoms. We calculated the phonon dispersions of LFO with ferrimagnetic and non-magnetic states, and subsequently evaluated the Gibbs free energy differences between the two states, plotting the P-T diagram for thermodynamic stability of the ferrimagnetic against non-magnetic state. From the P-T diagram, the Curie temperature is found to be ∼925 K at the normal condition and gradually increase with increasing pressure. Our calculations explain the experimental observations for material properties of LFO, providing a comprehensive understanding of the underlying mechanism and useful guidance for enhancing performance of LFO-based devices.

Designing isolation guidelines for COVID-19 patients utilizing rapid antigen tests: a simulation study using viral dynamics models.

Appropriate isolation guidelines for COVID-19 patients are warranted. Currently, isolating for fixed time is adapted in most countries. However, given the variability in viral dynamics between patients, some patients may no longer be infectious by the end of isolation (thus they are redundantly isolated), whereas others may still be infectious. Utilizing viral test results to determine ending isolation would minimize both the risk of ending isolation of infectious patients and the burden due to redundant isolation of noninfectious patients. In our previous study, we proposed a computational framework using SARS-CoV-2 viral dynamics models to compute the risk and the burden of different isolation guidelines with PCR tests. In this study, we extend the computational framework to design isolation guidelines for COVID-19 patients utilizing rapid antigen tests. Time interval of tests and number of consecutive negative tests to minimize the risk and the burden of isolation were explored. Furthermore, the approach was extended for asymptomatic cases. We found the guideline should be designed considering various factors: the infectiousness threshold values, the detection limit of antigen tests, symptom presence, and an acceptable level of releasing infectious patients. Especially, when detection limit is higher than the infectiousness threshold values, more consecutive negative results are needed to ascertain loss of infectiousness. To control the risk of releasing of infectious individuals under certain levels, rapid antigen tests should be designed to have lower detection limits than infectiousness threshold values to minimize the length of prolonged isolation, and the length of prolonged isolation increases when the detection limit is higher than the infectiousness threshold values, even though the guidelines are optimized for given conditions.

Incomplete antiviral treatment may induce longer durations of viral shedding during SARS-CoV-2 infection.

The duration of viral shedding is determined by a balance between de novo infection and removal of infected cells. That is, if infection is completely blocked with antiviral drugs (100% inhibition), the duration of viral shedding is minimal and is determined by the length of virus production. However, some mathematical models predict that if infected individuals are treated with antiviral drugs with efficacy below 100%, viral shedding may last longer than without treatment because further de novo infections are driven by entry of the virus into partially protected, uninfected cells at a slower rate. Using a simple mathematical model, we quantified SARS-CoV-2 infection dynamics in non-human primates and characterized the kinetics of viral shedding. We counterintuitively found that treatments initiated early, such as 0.5 d after virus inoculation, with intermediate to relatively high efficacy (30-70% inhibition of virus replication) yield a prolonged duration of viral shedding (by about 6.0 d) compared with no treatment.

Detection of significant antiviral drug effects on COVID-19 with reasonable sample sizes in randomized controlled trials: A modeling study.

BACKGROUND: Development of an effective antiviral drug for Coronavirus Disease 2019 (COVID-19) is a global health priority. Although several candidate drugs have been identified through in vitro and in vivo models, consistent and compelling evidence from clinical studies is limited. The lack of evidence from clinical trials may stem in part from the imperfect design of the trials. We investigated how clinical trials for antivirals need to be designed, especially focusing on the sample size in randomized controlled trials. METHODS AND FINDINGS: A modeling study was conducted to help understand the reasons behind inconsistent clinical trial findings and to design better clinical trials. We first analyzed longitudinal viral load data for Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) without antiviral treatment by use of a within-host virus dynamics model. The fitted viral load was categorized into 3 different groups by a clustering approach. Comparison of the estimated parameters showed that the 3 distinct groups were characterized by different virus decay rates (p-value < 0.001). The mean decay rates were 1.17 d-1 (95% CI: 1.06 to 1.27 d-1), 0.777 d-1 (0.716 to 0.838 d-1), and 0.450 d-1 (0.378 to 0.522 d-1) for the 3 groups, respectively. Such heterogeneity in virus dynamics could be a confounding variable if it is associated with treatment allocation in compassionate use programs (i.e., observational studies). Subsequently, we mimicked randomized controlled trials of antivirals by simulation. An antiviral effect causing a 95% to 99% reduction in viral replication was added to the model. To be realistic, we assumed that randomization and treatment are initiated with some time lag after symptom onset. Using the duration of virus shedding as an outcome, the sample size to detect a statistically significant mean difference between the treatment and placebo groups (1:1 allocation) was 13,603 and 11,670 (when the antiviral effect was 95% and 99%, respectively) per group if all patients are enrolled regardless of timing of randomization. The sample size was reduced to 584 and 458 (when the antiviral effect was 95% and 99%, respectively) if only patients who are treated within 1 day of symptom onset are enrolled. We confirmed the sample size was similarly reduced when using cumulative viral load in log scale as an outcome. We used a conventional virus dynamics model, which may not fully reflect the detailed mechanisms of viral dynamics of SARS-CoV-2. The model needs to be calibrated in terms of both parameter settings and model structure, which would yield more reliable sample size calculation. CONCLUSIONS: In this study, we found that estimated association in observational studies can be biased due to large heterogeneity in viral dynamics among infected individuals, and statistically significant effect in randomized controlled trials may be difficult to be detected due to small sample size. The sample size can be dramatically reduced by recruiting patients immediately after developing symptoms. We believe this is the first study investigated the study design of clinical trials for antiviral treatment using the viral dynamics model.

Revisiting the guidelines for ending isolation for COVID-19 patients.

Since the start of the COVID-19 pandemic, two mainstream guidelines for defining when to end the isolation of SARS-CoV-2-infected individuals have been in use: the one-size-fits-all approach (i.e. patients are isolated for a fixed number of days) and the personalized approach (i.e. based on repeated testing of isolated patients). We use a mathematical framework to model within-host viral dynamics and test different criteria for ending isolation. By considering a fixed time of 10 days since symptom onset as the criterion for ending isolation, we estimated that the risk of releasing an individual who is still infectious is low (0-6.6%). However, this policy entails lengthy unnecessary isolations (4.8-8.3 days). In contrast, by using a personalized strategy, similar low risks can be reached with shorter prolonged isolations. The obtained findings provide a scientific rationale for policies on ending the isolation of SARS-CoV-2-infected individuals.

Time variation in the probability of failing to detect a case of polymerase chain reaction testing for SARS-CoV-2 as estimated from a viral dynamics model.

Viral tests including polymerase chain reaction (PCR) tests are recommended to diagnose COVID-19 infection during the acute phase of infection. A test should have high sensitivity; however, the sensitivity of the PCR test is highly influenced by viral load, which changes over time. Because it is difficult to collect data before the onset of symptoms, the current literature on the sensitivity of the PCR test before symptom onset is limited. In this study, we used a viral dynamics model to track the probability of failing to detect a case of PCR testing over time, including the presymptomatic period. The model was parametrized by using longitudinal viral load data collected from 30 hospitalized patients. The probability of failing to detect a case decreased toward symptom onset, and the lowest probability was observed 2 days after symptom onset and increased afterwards. The probability on the day of symptom onset was 1.0% (95% CI: 0.5 to 1.9) and that 2 days before symptom onset was 60.2% (95% CI: 57.1 to 63.2). Our study suggests that the diagnosis of COVID-19 by PCR testing should be done carefully, especially when the test is performed before or way after symptom onset. Further study is needed of patient groups with potentially different viral dynamics, such as asymptomatic cases.