Specification testing for ordinary differential equation models with fixed design and applications to COVID-19 epidemic models.

Checking the models about the ongoing Coronavirus Disease 2019 (COVID-19) pandemic is an important issue. Some famous ordinary differential equation (ODE) models, such as the SIR and SEIR models have been used to describe and predict the epidemic trend. Still, in many cases, only part of the equations can be observed. A test is suggested to check possibly partially observed ODE models with a fixed design sampling scheme. The asymptotic properties of the test under the null, global and local alternative hypotheses are presented. Two new propositions about U-statistics with varying kernels based on independent but non-identical data are derived as essential tools. Some simulation studies are conducted to examine the performances of the test. Based on the available public data, it is found that the SEIR model, for modeling the data of COVID-19 infective cases in certain periods in Japan and Algeria, respectively, maybe not be appropriate by applying the proposed test.

Experimental study of the dispersion of cough-generated droplets from a person going up- or downstairs.

The dispersion of cough-generated droplets from a person going up- or downstairs was investigated through a laboratory experiment in a water tunnel. This experiment was carried out with a manikin mounted at inclination angles facing the incoming flow to mimic a person going up or down. Detailed velocity measurements and flow visualization were conducted in the water tunnel experiments. To investigate the influence of the initial position on the motion of particles, a virtual particle approach was adopted to simulate the dispersion of particles using the measured velocity field. Particle clustering, which is caused by the unsteadiness of the flow, was observed in both flow visualization and virtual particle simulation. For the case of going upstairs, particles are concentrated below the person's shoulder and move downward with a short travel distance. For the case of going downstairs, particles dispersing over the person's head advect over for a long distance. We also found that the motion of the particles is closely related to the initial position. According to the results in this study, suggestions for the prevention of respiratory infectious disease are made.

Dose-Dependent Variation of Synchronous Metabolites and Modules in a Yin/Yang Transformation Model of Appointed Ischemia Metabolic Networks.

AIM: Chinese medicine Danhong injection (DHI) is an effective pharmaceutical preparation for treating cerebral infarction. Our previous study shows that DHI can remarkably reduce the ischemic stroke-induced infarct volume in a dose-dependent manner, but the pharmacological mechanism of the DHI dose-dependent relationship is not clear. Therefore, the dose-dependent efficacy of DHI on cerebral ischemia and the underlying mechanisms were further investigated in this study. METHODS: A middle cerebral artery occlusion (MCAO) model was established and the rats were randomly divided into six groups: sham, vehicle, DHI dose-1, DHI dose-2, DHI dose-3, and DHI dose-4. Forty-one metabolites in serum were selected as candidate biomarkers of efficacy phenotypes by the Agilent 1290 rapid-resolution liquid chromatography system coupled with the Agilent 6550 Q-TOF MS system. Then, the metabolic networks in each group were constructed using the Weighted Correlation Network analysis (WGCNA). Moreover, the Yang and Yin transformation of six patterns (which are defined by up- and downregulation of metabolites) and synchronous modules divided from a synchronous network were used to dynamically analyze the mechanism of the drug's effectiveness. RESULTS: The neuroprotective effect of DHI has shown a dose-dependent manner, and the high-dose group (DH3 and DH4) effect is better. The entropy of the metabolic network and the Yin/Yang index both showed a consistent dose-response relationship. Seven dose-sensitive metabolites maintained constant inverse upregulation or downregulation in the four dose groups. Three synchronous modules for the DH1-DH4 full-course network were identified. Glycine, N-acetyl-L-glutamate, and tetrahydrofolate as a new emerging module appeared in DH2/DH3 and enriched in glutamine and glutamate metabolism-related pathways. CONCLUSION: This study takes the DHI metabolic network as an example to provide a new method for the discovery of multiple targets related to pharmacological effects. Our results show that the three conservative allosteric module nodes, taurine, L-tyrosine, and L-leucine, may be one of the basic mechanisms of DHI in the treatment of cerebral infarction, and the other three new emerging module nodes glyoxylate, L-glutamate, and L-valine may participate in the glutamine and glutamate metabolism pathway to improve the efficacy of DHI.

Effects of slope and speed of escalator on the dispersion of cough-generated droplets from a passenger.

During the pandemic of COVID-19, the public is encouraged to take stairs or escalators instead of elevators. However, the dispersion of respiratory droplets in these places, featured by slopes and human motion, is not well understood yet. It is consequently unclear whether the commonly recommended social-distancing guidelines are still appropriate in these scenarios. In this work, we analyze the dispersion of cough-generated droplets from a passenger riding an escalator with numerical simulations, focusing on the effects of the slope and speed of the escalator on the droplet dispersion. In the simulations, a one-way coupled Eulerian-Lagrangian approach is adopted, with the air-flow solved using the Reynolds-averaged Navier-Stokes method and the droplets modeled as passive Lagrangian particles. It is found that the slope alters the vertical concentration of the droplets in the passenger's wake significantly. The deflection of cough-generated jet and the wake flow behind the passenger drive the cough-generated droplets upwards when descending an escalator and downwards when ascending, resulting in both higher suspension height and larger spreading range of the viral droplets on a descending escalator than on an ascending one. These findings suggest that the present social-distancing guidelines may be inadequate on descending escalators and need further investigation.

The motion of respiratory droplets produced by coughing.

Coronavirus disease 2019 has become a global pandemic infectious respiratory disease with high mortality and infectiousness. This paper investigates respiratory droplet transmission, which is critical to understanding, modeling, and controlling epidemics. In the present work, we implemented flow visualization, particle image velocimetry, and particle shadow tracking velocimetry to measure the velocity of the airflow and droplets involved in coughing and then constructed a physical model considering the evaporation effect to predict the motion of droplets under different weather conditions. The experimental results indicate that the convection velocity of cough airflow presents the relationship t -0.7 with time; hence, the distance from the cougher increases by t 0.3 in the range of our measurement domain. Substituting these experimental results into the physical model reveals that small droplets (initial diameter D ≤ 100 µm) evaporate to droplet nuclei and that large droplets with D ≥ 500 µm and an initial velocity u 0 ≥ 5 m/s travel more than 2 m. Winter conditions of low temperature and high relative humidity can cause more droplets to settle to the ground, which may be a possible driver of a second pandemic wave in the autumn and winter seasons.

Defect-Rich Copper-doped Ruthenium Hollow Nanoparticles for Efficient Hydrogen Evolution Electrocatalysis in Alkaline Electrolyte.

It is of great importance to develop highly efficient and stable Pt-free catalysts for electrochemical hydrogen generation from water electrolysis. Here, monodisperse 7.5 nm copper-doped ruthenium hollow nanoparticles (NPs) with abundant defects and amorphous/crystalline hetero-phases were prepared and employed as efficient hydrogen evolution electrocatalysts in alkaline electrolyte. Specifically, these NPs only require a low overpotential of 25 mV to achieve a current density of 10 mA cm-2 in 1.0 M KOH and show acceptable stability after 2000 potential cycles, which represents one of the best Ru-based electrocatalysts for hydrogen evolution. Mechanism analysis indicates that Cu incorporation can modify the electronic structure of Ru shell, thereby optimizing the energy barrier for water adsorption and dissociation processes or H adsorption/desorption. Cu doping paired with the defect-rich and highly open hollow structure of the NPs greatly enhances hydrogen evolution activity.

Longitudinal study of BK Polyomavirus outcomes, risk factors, and kinetics in renal transplantation patients.

BACKGROUND: In the immunocompromised conditions following renal transplantation, BK virus can reactivate and cause BK virus associated nephropathy (BKVN). Increased BK viral loads and extended duration of infection have been linked to development of BKVN. The aim of this study was to observe the incidence of BKV infection and BKVN, and kinetics of infection and disease in renal transplantation recipients. METHODS: From 2014 to 2018, we conducted a longitudinal cohort observational study of 139 renal transplantation patients treated at a single clinic. Quantitative PCR assay was conducted to assess longitudinal BK viral loads. Analysis of patient clinical characteristics was performed to determine risk factors for BKV infection and associated disease. RESULTS: Of our cohort, 29 (20.9%) patients developed high BK viremia, and 7 (5.0%) developed biopsy-confirmed BKVN. Clinical parameters associated with diabetes (FBS, HbA1c) and hyperlipidemia (TG, TC, LDL-C) were found to be correlated with development of high BK viremia or BKVN. In 3 of 4 patients receiving intravenous immunoglobulin (IVIG) treatment, BK viral loads were reduced by at least 1 log within 2-3 months of administration. Significant differences were measured in BK viral loads and kidney function between BK viremic patients and BKVN patients by 3-9 months post-transplantation. CONCLUSIONS: We identified diabetes and hyperlipidemia as potential risk factors for development of high BK viremia and/or BKVN. IVIG was seen to be effective in reducing viral titers. The period 3-9 months post-transplantation was identified as important for development of BKVN from high BK viremia.

A Randomized Controlled Trial of Chinese Medicine on Nonmotor Symptoms in Parkinson's Disease.

Nonmotor symptoms (NMS) of Parkinson's disease (PD) have devastating impacts on both patients and their caregivers. Jiawei-Liujunzi Tang (JLT) has been used to treat some NMS of PD based on the Chinese medicine theory since Qing dynasty. Here we report a double-blind, randomized, placebo-controlled, add-on clinical trial aiming at evaluating the efficacy and safety of the JLT in treating NMS in PD patients. We randomly assigned 111 patients with idiopathic PD to receive either JLT or placebo for 32 weeks. Outcome measures were baseline to week 32 changes in Movement Disorder Society-Sponsored Revision of Unified PD Rating Scale (MDS-UPDRS) Parts I-IV and in NMS assessment scale for PD (NMSS). We observed improvements in the NMSS total score (p = 0.019), mood/cognition (p = 0.005), and reduction in hallucinations (p = 0.024). In addition, post hoc analysis showed a significant reduction in constipation (p < 0.001). However, there was no evidence of improvement in MDS-UPDRS Part I total score (p = 0.216) at week 32. Adverse events (AEs) were mild and comparable between the two groups. In conclusion, long-term administration of JLT is well tolerated and shows significant benefits in improving NMS including mood, cognition, and constipation.

Before and after - Nutritional transformation of dysmorphism in a case of Costello syndrome.

Costello syndrome is a type of RASopathy mapped to HRAS gene in chromosome 11, characterized by prenatal overgrowth, postnatal failure to thrive, classic facial gestalt and multisystem involvement including cardiomyopathy and intellectual disability. We present a 7 months old child with severe failure to thrive whose "subtle" facial dysmorphism at the time eluded clinical recognition of the syndrome. It was only with optimization of his nutritional status that dysmorphic features became more apparent, which affirmed the molecular diagnosis of Costello syndrome from exome sequencing. The case illustrated how drastic failure to thrive can be in Costello syndrome, and how nutritional status can transform dysmorphic features in a child. It also highlights the importance of serial dysmorphic evaluation in difficult cases.

Single-index varying coefficient model for functional responses.

Recently, massive functional data have been widely collected over space across a set of grid points in various imaging studies. It is interesting to correlate functional data with various clinical variables, such as age and gender, in order to address scientific questions of interest. The aim of this article is to develop a single-index varying coefficient (SIVC) model for establishing a varying association between functional responses (e.g., image) and a set of covariates. It enjoys several unique features of both varying-coefficient and single-index models. An estimation procedure is developed to estimate varying coefficient functions, the index function, and the covariance function of individual functions. The optimal integration of information across different grid points is systematically delineated and the asymptotic properties (e.g., consistency and convergence rate) of all estimators are examined. Simulation studies are conducted to assess the finite-sample performance of the proposed estimation procedure. Furthermore, our real data analysis of a white matter tract dataset obtained from the Alzheimer's Disease Neuroimaging Initiative (ADNI) study confirms the advantage and accuracy of SIVC model over the popular varying coefficient model.

Functional Nonlinear Mixed Effects Models for Longitudinal Image Data.

Motivated by studying large-scale longitudinal image data, we propose a novel functional nonlinear mixed effects modeling (FNMEM) framework to model the nonlinear spatial-temporal growth patterns of brain structure and function and their association with covariates of interest (e.g., time or diagnostic status). Our FNMEM explicitly quantifies a random nonlinear association map of individual trajectories. We develop an efficient estimation method to estimate the nonlinear growth function and the covariance operator of the spatial-temporal process. We propose a global test and a simultaneous confidence band for some specific growth patterns. We conduct Monte Carlo simulation to examine the finite-sample performance of the proposed procedures. We apply FNMEM to investigate the spatial-temporal dynamics of white-matter fiber skeletons in a national database for autism research. Our FNMEM may provide a valuable tool for charting the developmental trajectories of various neuropsychiatric and neurodegenerative disorders.

Covariance-enhanced discriminant analysis.

Linear discriminant analysis has been widely used to characterize or separate multiple classes via linear combinations of features. However, the high dimensionality of features from modern biological experiments defies traditional discriminant analysis techniques. Possible interfeature correlations present additional challenges and are often underused in modelling. In this paper, by incorporating possible interfeature correlations, we propose a covariance-enhanced discriminant analysis method that simultaneously and consistently selects informative features and identifies the corresponding discriminable classes. Under mild regularity conditions, we show that the method can achieve consistent parameter estimation and model selection, and can attain an asymptotically optimal misclassification rate. Extensive simulations have verified the utility of the method, which we apply to a renal transplantation trial.

Ultrahigh dimensional time course feature selection.

Statistical challenges arise from modern biomedical studies that produce time course genomic data with ultrahigh dimensions. In a renal cancer study that motivated this paper, the pharmacokinetic measures of a tumor suppressor (CCI-779) and expression levels of 12,625 genes were measured for each of 33 patients at 8 and 16 weeks after the start of treatments, with the goal of identifying predictive gene transcripts and the interactions with time in peripheral blood mononuclear cells for pharmacokinetics over the time course. The resulting data set defies analysis even with regularized regression. Although some remedies have been proposed for both linear and generalized linear models, there are virtually no solutions in the time course setting. As such, a novel GEE-based screening procedure is proposed, which only pertains to the specifications of the first two marginal moments and a working correlation structure. Different from existing methods that either fit separate marginal models or compute pairwise correlation measures, the new procedure merely involves making a single evaluation of estimating functions and thus is extremely computationally efficient. The new method is robust against the mis-specification of correlation structures and enjoys theoretical readiness, which is further verified via Monte Carlo simulations. The procedure is applied to analyze the aforementioned renal cancer study and identify gene transcripts and possible time-interactions that are relevant to CCI-779 metabolism in peripheral blood.

Chronic noise exposure causes persistence of tau hyperphosphorylation and formation of NFT tau in the rat hippocampus and prefrontal cortex.

The non-auditory effects of noise exposure on the central nervous system have been established both epidemiologically and experimentally. Chronic noise exposure (CNE) has been associated with tau hyperphosphorylation and Alzheimer's disease (AD)-like pathological changes. However, experimental evidence for these associations remains limited. The aim of the current study was to explore the effects of CNE [100 dB sound pressure level (SPL) white noise, 4 h/d×14 d] on tau phosphorylation in the rat hippocampus and the prefrontal cortex. Forty-eight male Wistar rats were randomly assigned to two groups: a noise-exposed group and a control group. The levels of radioimmunoprecipitation assay (RIPA)-soluble and RIPA-insoluble phosphorylated tau at Ser202, Ser396, Ser404, and Ser422 in the hippocampus and the prefrontal cortex were measured at different time points (days 0, 3, 7, and 14) after the end of the last noise exposure. Exposure to white noise for 14 consecutive days significantly increased the levels of tau phosphorylation at Ser202, Ser396, Ser404, and Ser422, the sites typically phosphorylated in AD brains, in the hippocampus and the prefrontal cortex. Tau hyperphosphorylation persisted for 7 to 14 d after the cessation of noise exposure. These alterations were also concomitant with the generation of pathological neurofibrillary tangle (NFT) tau 3, 7 and 14 d after the end of the stimulus. Furthermore, lasting increases in proteins involved in hyperphosphorylation, namely glycogen synthase kinase 3ß (GSK3ß) and protein phosphatase 2A (PP2A), were found to occur in close correspondence with increase in tau hyperphosphorylation. The results of this study show that CNE leads to long-lasting increases in non-NFT hyperphosphorylated tau and delayed formation of misfolded NFT tau in the hippocampus and the prefrontal cortex. Our results also provide evidence for the involvement of GSK3ß and PP2A in these processes.

A Two-Stage Estimation Method for Random Coefficient Differential Equation Models with Application to Longitudinal HIV Dynamic Data.

We propose a two-stage estimation method for random coefficient ordinary differential equation (ODE) models. A maximum pseudo-likelihood estimator (MPLE) is derived based on a mixed-effects modeling approach and its asymptotic properties for population parameters are established. The proposed method does not require repeatedly solving ODEs, and is computationally efficient although it does pay a price with the loss of some estimation efficiency. However, the method does offer an alternative approach when the exact likelihood approach fails due to model complexity and high-dimensional parameter space, and it can also serve as a method to obtain the starting estimates for more accurate estimation methods. In addition, the proposed method does not need to specify the initial values of state variables and preserves all the advantages of the mixed-effects modeling approach. The finite sample properties of the proposed estimator are studied via Monte Carlo simulations and the methodology is also illustrated with application to an AIDS clinical data set.

Asymmetric interaction and indeterminate fitness correlation between cooperative partners in the fig-fig wasp mutualism.

Empirical observations have shown that cooperative partners can compete for common resources, but what factors determine whether partners cooperate or compete remain unclear. Using the reciprocal fig-fig wasp mutualism, we show that nonlinear amplification of interference competition between fig wasps-which limits the fig wasps' ability to use a common resource (i.e. female flowers)-keeps the common resource unsaturated, making cooperation locally stable. When interference competition was manually prevented, the fitness correlation between figs and fig wasps went from positive to negative. This indicates that genetic relatedness or reciprocal exchange between cooperative players, which could create spatial heterogeneity or self-restraint, was not sufficient to maintain stable cooperation. Moreover, our analysis of field-collected data shows that the fitness correlation between cooperative partners varies stochastically, and that the mainly positive fitness correlation observed during the warm season shifts to a negative correlation during the cold season owing to an increase in the initial oviposition efficiency of each fig wasp. This implies that the discriminative sanction of less-cooperative wasps (i.e. by decreasing the egg deposition efficiency per fig wasp) but reward to cooperative wasps by fig, a control of the initial value, will facilitate a stable mutualism. Our finding that asymmetric interaction leading to an indeterminate fitness interaction between symbiont (i.e. cooperative actors) and host (i.e. recipient) has the potential to explain why conflict has been empirically observed in both well-documented intraspecific and interspecific cooperation systems.

Model-Free Feature Screening for Ultrahigh Dimensional Data.

With the recent explosion of scientific data of unprecedented size and complexity, feature ranking and screening are playing an increasingly important role in many scientific studies. In this article, we propose a novel feature screening procedure under a unified model framework, which covers a wide variety of commonly used parametric and semiparametric models. The new method does not require imposing a specific model structure on regression functions, and thus is particularly appealing to ultrahigh-dimensional regressions, where there are a huge number of candidate predictors but little information about the actual model forms. We demonstrate that, with the number of predictors growing at an exponential rate of the sample size, the proposed procedure possesses consistency in ranking, which is both useful in its own right and can lead to consistency in selection. The new procedure is computationally efficient and simple, and exhibits a competent empirical performance in our intensive simulations and real data analysis.

Efficacy of a Chinese herbal proprietary medicine (Hemp Seed Pill) for functional constipation.

OBJECTIVES: Functional constipation (FC) is a common clinical complaint. Despite a lack of consolidated evidence, Chinese herbal medicine (CHM) has become a popular alternative treatment for this condition. The aim of this study was to assess, with a rigidly designed study, the efficacy and safety of a CHM proprietary medicine, Hemp Seed Pill (HSP), in optimal dosage for treating FC. METHODS: This study comprised two parts: trial I, a dose determination study, and trial II, a placebo-controlled clinical study. In trial I, the optimal dosage of HSP was first determined from among three doses (2.5, 5.0, and 7.5 g b.i.d.). In trial II, a randomized double-blind study, the efficacy and safety of HSP for FC patients (Rome III criteria) in excessive syndrome as defined by traditional Chinese medicine (TCM) theory were compared with placebo. All participants in trials underwent a 2-week run-in, an 8-week treatment, and an 8-week follow-up. The primary end point was the responder rate for complete spontaneous bowel movement (CSBM) during treatment. Participants with a mean increase of CSBM ≧ 1/week compared with their baselines were defined as responders. Secondary outcome measures included responder rate during follow-up, individual and global symptom assessments, and reported adverse effects (AEs). RESULTS: The dose of 7.5 g b.i.d. showed better therapeutic effect than that of 2.5 and 5.0 g b.i.d. among 96 subjects (32 per arm) in trial I and was therefore selected for comparison with placebo in trial II. In trial II, 120 subjects were randomized into two arms (60 per arm). Responder rates for the HSP and placebo groups were 43.3 and 8.3% during treatment and 30.0 and 15.0% in the follow-up period, respectively (P<0.05). Those in the HSP group showed benefit in terms of increased CSBM, relief in the severity of constipation and straining of evacuation, and effective reduction in the use of rescue therapy when compared with placebo. No serious AE was reported. CONCLUSIONS: HSP (7.5 g b.i.d.) is safe and effective for alleviating FC for subjects in excessive syndrome. Optimal dose determination may be crucial for all CHM studies.

[Treatment of focal bone defect in postoperative nonunion with autologous red bone marrow injection].

OBJECTIVE: To observe the clinical effect of autologous red bone marrow injection in treating focal bone defect in postoperative nonunion. METHODS: Thirteen patients with focal bone defect in postoperative nonunion (7 cases in tibia, 2 cases in femur, 4 cases in humerus), including 8 males and 5 females with the mean age of 32.5-years-old (ranging from 15 to 60 years). The bone defects were treated with autologous red bone marrow injection (1 time per 2 weeks, 5 times in total) and the X-rays of AP and LP were observed. RESULTS: Thirteen patients were followed up from 6 to 12 months with an average of 7.5 months. According to results of X-ray pictures, 13 cases obtained bone defect recovered completely, and the average time of union was 4 months. CONCLUSION: Autologous red bone marrow injection has ascendancy such as less wound and clear clinical effect, which can accelerate bone healing and promotes functional recovery of limb. It is a good method to treat focal bone defect in postoperative nonunion.

[Risk factors for capillary leak syndrome in children with tetralogy after operation].

OBJECTIVE: To determine risk factors of capillary leak syndrome(CLS) in children with tetralogy after operation. METHODS: Clinical data were retrospectively collected and analyzed from 32 tetralogy cases with CLS and 50 cases without CLS(control group), who received operation under cardiopulmonary bypass (CBP) in our hospital from October 2002 to December 2006. Risk factors with statistical significance were screened with univariate logistic regression analysis, independent risk factors of CLS were determined with multivariate logistic regression analysis. Postoperative outcome was compared between CLS group and control group. RESULT: Logistic analysis showed that the risk factors for CLS were age(OR=6.783), duration of CBP(OR=4.756)and MGoon index (OR=3.826). There were statistical differences in injection of colloid, time of inotropic drugs and ventilation between CLS and control groups(P<0.01). Eight CLS cases underwent peritoneal dialysis and 2 CLS cases died. CONCLUSION: The risk factors of CLS in children with tetralogy after CBP are age, duration of CBP and MGoon index.