Submaximal cardiopulmonary thresholds on a robotics-assisted tilt table, a cycle and a treadmill: a comparative analysis

BACKGROUND: The robotics-assisted tilt table (RATT), including actuators for tilting and cyclical leg movement, is used for rehabilitation of severely disabled neurological patients. Following further engineering development of the system, i.e. the addition of force sensors and visual bio-feedback, patients can actively participate in exercise testing and training on the device. Peak cardiopulmonary performance parameters were previously investigated, but it also important to compare submaximal parameters with standard devices. The aim of this study was to evaluate the feasibility of the RATT for estimation of submaximal exercise thresholds by comparison with a cycle ergometer and a treadmill. METHODS: 17 healthy subjects randomly performed six maximal individualized incremental exercise tests, with two tests on each of the three exercise modalities. The ventilatory anaerobic threshold (VAT) and respiratory compensation point (RCP) were determined from breath-by-breath data. RESULTS: VAT and RCP on the RATT were lower than the cycle ergometer and the treadmill: oxygen uptake (V'O2) at VAT was [mean (SD)] 1.2 (0.3), 1.5 (0.4) and 1.6 (0.5) L/min, respectively (p < 0.001); V'O2 at RCP was 1.7 (0.4), 2.3 (0.8) and 2.6 (0.9) L/min, respectively (p = 0.001). High correlations for VAT and RCP were found between the RATT vs the cycle ergometer and RATT vs the treadmill (R on the range 0.69-0.80). VAT and RCP demonstrated excellent test-retest reliability for all three devices (ICC from 0.81 to 0.98). Mean differences between the test and retest values on each device were close to zero. The ventilatory equivalent for O2 at VAT for the RATT and cycle ergometer were similar and both were higher than the treadmill. The ventilatory equivalent for CO2 at RCP was similar for all devices. Ventilatory equivalent parameters demonstrated fair-to-excellent reliability and repeatability. CONCLUSIONS: It is feasible to use the RATT for estimation of submaximal exercise thresholds: VAT and RCP on the RATT were lower than the cycle ergometer and the treadmill, but there were high correlations between the RATT vs the cycle ergometer and vs the treadmill. Repeatability and test-retest reliability of all submaximal threshold parameters from the RATT were comparable to those of standard devices.

Cell cycle reactivation of cochlear progenitor cells in neonatal FUCCI mice by a GSK3 small molecule inhibitor.

Due to the lack of regenerative capacity of the mammalian auditory epithelium, sensory hair cell loss results in permanent hearing deficit. Nevertheless, a population of tissue resident stem/progenitor cells has been recently described. Identification of methods to trigger their activity could lead to exploitation of their potential therapeutically. Here we validate the use of transgenic mice reporting cell cycle progression (FUCCI), and stemness (Lgr5-GFP), as a valuable tool to identify regulators of cell cycle re-entry of supporting cells within the auditory epithelium. The small molecule compound CHIR99021 was used to inhibit GSK3 activity. This led to a significant increase in the fraction of proliferating sphere-forming cells, labeled by the FUCCI markers and in the percentage of Lgr5-GFP + cells, as well as a selective increase in the fraction of S-G2-M cells in the Lgr5 + population. Using whole mount cultures of the organ of Corti we detected a statistically significant increment in the fraction of proliferating Sox2 supporting cells after CHIR99021 treatment, but only rarely appearance of novel MyoVIIa+/Edu + hair cells. In conclusion, these tools provide a robust mean to identify novel regulators of auditory organ regeneration and to clarify the contribution of stem cell activity.

Sleep-Wake Cycle Dysfunction in the TgCRND8 Mouse Model of Alzheimer's Disease: From Early to Advanced Pathological Stages.

In addition to cognitive decline, individuals affected by Alzheimer's disease (AD) can experience important neuropsychiatric symptoms including sleep disturbances. We characterized the sleep-wake cycle in the TgCRND8 mouse model of AD, which overexpresses a mutant human form of amyloid precursor protein resulting in high levels of β-amyloid and plaque formation by 3 months of age. Polysomnographic recordings in freely-moving mice were conducted to study sleep-wake cycle architecture at 3, 7 and 11 months of age and corresponding levels of β-amyloid in brain regions regulating sleep-wake states were measured. At all ages, TgCRND8 mice showed increased wakefulness and reduced non-rapid eye movement (NREM) sleep during the resting and active phases. Increased wakefulness in TgCRND8 mice was accompanied by a shift in the waking power spectrum towards fast frequency oscillations in the beta (14-20 Hz) and low gamma range (20-50 Hz). Given the phenotype of hyperarousal observed in TgCRND8 mice, the role of noradrenergic transmission in the promotion of arousal, and previous work reporting an early disruption of the noradrenergic system in TgCRND8, we tested the effects of the alpha-1-adrenoreceptor antagonist, prazosin, on sleep-wake patterns in TgCRND8 and non-transgenic (NTg) mice. We found that a lower dose (2 mg/kg) of prazosin increased NREM sleep in NTg but not in TgCRND8 mice, whereas a higher dose (5 mg/kg) increased NREM sleep in both genotypes, suggesting altered sensitivity to noradrenergic blockade in TgCRND8 mice. Collectively our results demonstrate that amyloidosis in TgCRND8 mice is associated with sleep-wake cycle dysfunction, characterized by hyperarousal, validating this model as a tool towards understanding the relationship between β-amyloid overproduction and disrupted sleep-wake patterns in AD.

Regulation of histone mRNA in the unperturbed cell cycle: evidence suggesting control at two posttranscriptional steps.

The levels of histone mRNA increase 35-fold as selectively detached mitotic CHO cells progress from mitosis through G1 and into S phase. Using an exogenous gene with a histone 3' end which is not sensitive to transcriptional or half-life regulation, we show that 3' processing is regulated as cells progress from G1 to S phase. The half-life of histone mRNA is similar in G1- and S-phase cells, as measured after inhibition of transcription by actinomycin D (dactinomycin) or indirectly after stabilization by the protein synthesis inhibitor cycloheximide. Taken together, these results suggest that the change in histone mRNA levels between G1- and S-phase cells must be due to an increase in the rate of biosynthesis, a combination of changes in transcription rate and processing efficiency. In G2 phase, there is a rapid 35-fold decrease in the histone mRNA concentration which our results suggest is due primarily to an altered stability of histone mRNA. These results are consistent with a model for cell cycle regulation of histone mRNA levels in which the effects on both RNA 3' processing and transcription, rather than alterations in mRNA stability, are the major mechanisms by which low histone mRNA levels are maintained during G1.

Three-dimensional culture and characterization of mononuclear cells from human bone marrow

BACKGROUND AIMS The diverse phenotypic changes and clinical and economic disadvantages associated with the monolayer expansion of bone marrow-derived mesenchymal stromal cells (MSCs) have focused attention on the development of one-step intraoperative cells therapies and homing strategies. The mononuclear cell fraction of bone marrow, inclusive of discrete stem cell populations, is not well characterized, and we currently lack suitable cell culture systems in which to culture and investigate the behavior of these cells. METHODS Human bone marrow-derived mononuclear cells were cultured within fibrin for 2 weeks with or without fibroblast growth factor-2 supplementation. DNA content and cell viability of enzymatically retrieved cells were determined at days 7 and 14. Cell surface marker profiling and cell cycle analysis were performed by means of multi-color flow cytometry and a 5-ethynyl-2'-deoxyuridine incorporation assay, respectively. RESULTS Total mononuclear cell fractions, isolated from whole human bone marrow, was successfully cultured in fibrin gels for up to 14 days under static conditions. Discrete niche cell populations including MSCs, pericytes and hematopoietic stem cells were maintained in relative quiescence for 7 days in proportions similar to that in freshly isolated cells. Colony-forming unit efficiency of enzymatically retrieved MSCs was significantly higher at day 14 compared to day 0; and in accordance with previously published works, it was fibroblast growth factor-2-dependant. CONCLUSIONS Fibrin gels provide a simple, novel system in which to culture and study the complete fraction of bone marrow-derived mononuclear cells and may support the development of improved bone marrow cell-based therapies.

Complete genome of a Puumala virus strain from Central Europe

Puumala virus (PUUV) is one of the predominant hantavirus species in Europe causing mild to moderate cases of haemorrhagic fever with renal syndrome. Parts of Lower Saxony in north-western Germany are endemic for PUUV infections. In this study, the complete PUUV genome sequence of a bank vole-derived tissue sample from the 2007 outbreak was determined by a combined primer-walking and RNA ligation strategy. The S, M and L genome segments were 1,828, 3,680 and 6,550 nucleotides in length, respectively. Sliding-window analyses of the nucleotide sequences of all available complete PUUV genomes indicated a non-homogenous distribution of variability with hypervariable regions located at the 3′-ends of the S and M segments. The overall similarity of the coding genome regions to the other PUUV strains ranged between 80.1 and 84.7 % at the level of the nucleotide sequence and between 89.5 and 98.1 % for the deduced amino acid sequences. In comparison to the phylogenetic trees of the complete coding sequences, trees based on partial segments revealed a general drop in phylogenetic support and a lower resolution. The Astrup strain S and M segment sequences showed the highest similarity to sequences of strains from geographically close sites in the Osnabrück Hills region. In conclusion, a primer-walking-mediated strategy resulted in the determination of the first complete nucleotide sequence of a PUUV strain from Central Europe. Different levels of variability along the genome provide the opportunity to choose regions for analyses according to the particular research question, e.g., large-scale phylogenetics or within-host evolution.

Variability of mesospheric water vapor above Bern in relation to the 27-day solar rotation cycle

Many studies investigated solar–terrestrial responses (thermal state, O₃ , OH, H₂O) with emphasis on the tropical upper atmosphere. In this paper the Focus is switched to water vapor in the mesosphere at a mid-latitudinal location. Eight years of water vapor profile measurements above Bern (46.88°N/7.46°E) are investigated to study oscillations with the Focus on periods between 10 and 50 days. Different spectral analyses revealed prominent features in the 27-day oscillation band, which are enhanced in the upper mesosphere (above 0.1 hPa, ∼64 km) during the rising sun spot activity of solar cycle 24. Local as well as zonal mean Aura MLS observations Support these results by showing a similar behavior. The relationship between mesospheric water and the solar Lyman-α flux is studied by comparing thesi-milarity of their temporal oscillations. The H₂O oscillation is negatively correlated to solar Lyman-α oscillation with a correlation coefficient of up to −0.3 to −0.4, and the Phase lag is 6–10 days at 0.04 hPa. The confidence level of the correlation is ≥99%. This finding supports the assumption that the 27-day oscillation in Lyman-α causes a periodical photo dissociation loss in mesospheric water. Wavelet power spectra, cross-wavelet transform and wavelet coherence analysis (WTC)complete our study. More periods of high common wavelet power of H₂O and solar Lyman-α are present when amplitudes of the Lyman-α flux increase. Since this is not a measure of physical correlation a more detailed view on WTC is necessary, where significant (two sigma level)correlations occur intermittently in the 27 and 13-day band with variable Phase lock behavior. Large Lyman-α oscillations appeared after the solar super storm in July 2012 and the H₂O oscillations show a well pronounced anticorrelation. The competition between advective transport and photo dissociation loss of mesospheric water vapor may explain the sometimes variable Phase relationship of mesospheric H₂O and solar Lyman-α oscillations. Generally, the WTC analysis indicates that solar variability causes observable photochemical and dynamical processes in the mid-latitude mesosphere.

Sulfate assimilation in higher plants. Characterization of a stable intermediate in the adenosine 5′-phosphosulfate reductase reaction

The enzyme catalysing the reduction of adenosine 5′-phosphosulfate (AdoPS) to sulfite in higher plants, AdoPS reductase, is considered to be the key enzyme of assimilatory sulfate reduction. In order to address its reaction mechanism, the APR2 isoform of this enzyme from Arabidopsis thaliana was overexpressed in Escherichia coli and purified to homogeneity. Incubation of the enzyme with [35S]AdoPS at 4 °C resulted in radioactive labelling of the protein. Analysis of APR2 tryptic peptides revealed 35SO2–3 bound to Cys248, the only Cys conserved between AdoPS and prokaryotic phosphoadenosine 5′-phosphosulfate reductases. Consistent with this result, radioactivity could be released from the protein by incubation with thiols, inorganic sulfide and sulfite. The intermediate remained stable, however, after incubation with sulfate, oxidized glutathione or AdoPS. Because truncated APR2, missing the thioredoxin-like C-terminal part, could be labelled even at 37 °C, and because this intermediate was more stable than the complete protein, we conclude that the thioredoxin-like domain was required to release the bound SO2–3 from the intermediate. Taken together, these results demonstrate for the first time the binding of 35SO2–3 from [35S]AdoPS to AdoPS reductase and its subsequent release, and thus contribute to our understanding of the molecular mechanism of AdoPS reduction in plants.

Chromium uptake and adsorption in marine phytoplankton - Implications for the marine chromium cycle

Using the radioisotope 51Cr, we investigated the controls of cellular Cr accumulation in an array of marine phytoplankton grown in environmentally relevant Cr concentrations (1–10 nM). Given the affinity of Cr(III) for amorphous Fe-hydroxide mineral surfaces, and the formation of these mineral phases on the outside of phytoplankton cells, extracellular Cr was monitored in a model diatom species (Thalassiosira weissflogii) as extracellular Fe concentrations varied. Extracellular Cr in T. weissflogii increased with increasing extracellular Fe, demonstrating that Cr may be removed from seawater via extracellular adsorption to phytoplankton. Short-term Cr(VI) and Cr(III) uptake experiments performed with T. weissflogii demonstrated that Cr(III) was the primary oxidation state adsorbing to cells and being internalized by them. Cellular Cr:C ratios (<0.5 μmol Cr mol C−1) of the eight phytoplankton species surveyed were significantly lower than previously reported Cr:C ratios in marine particles with a high biogenic component (10–300 μmol Cr mol C−1). This indicates that Cr(III) likely accumulates in marine particles due to uptake and/or adsorption. Mass balance calculations demonstrate that surface water Cr deficits can be explained via loss of Cr(III) to exported particles, thereby providing a mechanism to account for the nutrient depth profile for Cr in modern seawater. Given the large fractionation of stable Cr isotopes during Cr(VI) reduction, Cr(III) associated with exported organic carbon is likely enriched in lighter isotopes. Most sedimentary Cr isotope studies have thus far neglected internal fractionating processes in the marine Cr cycle, but our data indicate that loss of Cr to exported particles may be traced in the sedimentary d53Cr record.

TAC102 Is a Novel Component of the Mitochondrial Genome Segregation Machinery in Trypanosomes

Trypanosomes show an intriguing organization of their mitochondrial DNA into a catenated network, the kinetoplast DNA (kDNA). While more than 30 proteins involved in kDNA replication have been described, only few components of kDNA segregation machinery are currently known. Electron microscopy studies identified a high-order structure, the tripartite attachment complex (TAC), linking the basal body of the flagellum via the mitochondrial membranes to the kDNA. Here we describe TAC102, a novel core component of the TAC, which is essential for proper kDNA segregation during cell division. Loss of TAC102 leads to mitochondrial genome missegregation but has no impact on proper organelle biogenesis and segregation. The protein is present throughout the cell cycle and is assembled into the newly developing TAC only after the pro-basal body has matured indicating a hierarchy in the assembly process. Furthermore, we provide evidence that the TAC is replicated de novo rather than using a semi-conservative mechanism. Lastly, we demonstrate that TAC102 lacks an N-terminal mitochondrial targeting sequence and requires sequences in the C-terminal part of the protein for its proper localization.

Seasonal Distributions of Satellite-Measured Phytoplankton Pigment Concentration Along the Chilean Coast

Five years (1979-1983) of Coastal Zone Color Scanner satellite ocean color data are used to examine seasonal patterns of phytoplankton pigment concentration along the Chilean coast from 20 degrees S to 45 degrees S. Four kilometer resolution, 2-4 day composites document the presence of filaments of elevated pigment concentration extending offshore throughout the study area, with maximum offshore extension at higher latitudes. In three years, 1979, 1981, and 1983, sufficient data exist in monthly composites to allow recreation of portions of the seasonal cycle. Data in 1979 are the most complete. Near-shore concentrations and cross-shelf extension of pigment concentrations in 1979 are maximum in austral winter throughout the study area and minimum in summer. Available data from 1981 and 1983 are consistent with this temporal pattern but with concentrations approximately double those of 1979. Seasonal, spatial patterns within 10 km of shore and 50 km offshore indicate a latitudinal discontinuity both in absolute concentration and in the magnitude of the seasonal cycle at approximately 33 degrees S in both 1979 and in the climatological time series. The discontinuity is strongest ill fall-winter and weakest in summer. South of this latitude, concentrations are relatively high (2-3 mg m(-3) in 1979), a strong seasonal cycle is present, and patterns 50 km offshore are correlated with those within 10 km of shore. North of 33 degrees S, concentrations are < 1.5 mg m(-3) (in 1979), and the seasonal cycle within 10 km of shore is present but much weaker and less obviously correlated with that 50 km offshore. The seasonal cycle of pigment concentrations is 180 degrees out of phase with monthly averaged upwelling favorable winds. Noncoincident Pathfinder sea surface temperature data show that over most latitudes, coastal low surface temperatures lag wind forcing by 1-2 months, but these too are out of phase with the pigment seasonal cycle. These data point to control of pigment patterns along the Chilean coast by the interaction of upwelling with circulation patterns unconnected to local wind forcing.

The Impact of Whaling on the Ocean Carbon Cycle: Why Bigger Was Better

Background: Humans have reduced the abundance of many large marine vertebrates, including whales, large fish, and sharks, to only a small percentage of their pre-exploitation levels. Industrial fishing and whaling also tended to preferentially harvest the largest species and largest individuals within a population. We consider the consequences of removing these animals on the ocean's ability to store carbon. Methodology/Principal Findings: Because body size is critical to our arguments, our analysis focuses on populations of baleen whales. Using reconstructions of pre-whaling and modern abundances, we consider the impact of whaling on the amount of carbon stored in living whales and on the amount of carbon exported to the deep sea by sinking whale carcasses. Populations of large baleen whales now store 9.1 x 10(6) tons less carbon than before whaling. Some of the lost storage has been offset by increases in smaller competitors; however, due to the relative metabolic efficiency of larger organisms, a shift toward smaller animals could decrease the total community biomass by 30% or more. Because of their large size and few predators, whales and other large marine vertebrates can efficiently export carbon from the surface waters to the deep sea. We estimate that rebuilding whale populations would remove 1.6 x 10(5) tons of carbon each year through sinking whale carcasses. Conclusions/Significance: Even though fish and whales are only a small portion of the ocean's overall biomass, fishing and whaling have altered the ocean's ability to store and sequester carbon. Although these changes are small relative to the total ocean carbon sink, rebuilding populations of fish and whales would be comparable to other carbon management schemes, including ocean iron fertilization.

Molecular consequences of the loss of 3'→5' exonuclease activity of DNA polymerases delta and epsilon

The DNA replication polymerases δ and ϵ have an inherent proofreading mechanism in the form of a 3'→5' exonuclease. Upon recognition of errant deoxynucleotide incorporation into DNA, the nascent primer terminus is partitioned to the exonuclease active site where the incorrectly paired nucleotide is excised before resumption of polymerization. The goal of this project was to identify the cellular and molecular consequences of an exonuclease deficiency. The proofreading capability of model system MEFs with EXOII mutations was abolished without altering polymerase function.^ It was hypothesized that 3'→5' exonucleases of polymerases δ and ϵ are critical for prevention of replication stress and important for sensitization to nucleoside analogs. To test this hypothesis, two aims were formulated: Determine the effect of the exonuclease active site mutation on replication related molecular signaling and identify the molecular consequences of an exonuclease deficiency when replication is challenged with nucleoside analogs.^ Via cell cycle studies it was determined that larger populations of exonuclease deficient cells are in the S-phase. There was an increase in levels of replication proteins, cell population growth and DNA synthesis capacity without alteration in cell cycle progression. These findings led to studies of proteins involved in checkpoint activation and DNA damage sensing. Finally, collective modifications at the level of DNA replication likely affect the strand integrity of DNA at the chromosomal level.^ Gemcitabine, a DNA directed nucleoside analog is a substrate of polymerases δ and ϵ and exploits replication to become incorporated into DNA. Though accumulation of gemcitabine triphosphate was similar in all cell types, incorporation into DNA and rates of DNA synthesis were increased in exonuclease defective cells and were not consistent with clonogenic survival. This led to molecular signaling investigations which demonstrated an increase in S-phase cells and activation of a DNA damage response upon gemcitabine treatment.^ Collectively, these data indicate that the loss of exonuclease results in a replication stress response that is likely required to employ other repair mechanisms to remove unexcised mismatches introduced into DNA during replication. When challenged with nucleoside analogs, this ongoing stress response coupled with repair serves as a resistance mechanism to cell death.^

SYSTEMIC TOXICITY IN RODENTS FOLLOWING ACUTE AND SUBCHRONIC INHALATION OF FORMALDEHYDE (FLOW CYTOMETRY, CELL CYCLE KINETICS, ALKALINE ELUTION)

Systemic toxicity was evaluated in Sprague-Dawley (SD) rats and A-strain mice exposed to HCHO inhalation at 0, 0.5, 3, or 15 ppm for six hours/day, five days/week for up to 24 weeks. Toxicity was measured by flow cytometry to detect changes in cell cycle RNA and DNA content and by alkaline elution to detect DNA protein cross-link (DPC) formation.^ A G(,2)M block was detected in SD rat marrow following one week of exposure to 0.5, 3, or 15 ppm HCHO, but this block did not persist. No effect was noticed in mouse marrow. Only a minimal increase in RNA content was detected in rat or mouse marrow while exfoliated lung cells showed a significant increase in RNA activity after one week of exposure.^ Acute exposure in SD rats for four hours/day for one or three days at 150 ppm showed an increase in RNA activity in exfoliated lung cells but not in the marrow after one day. On the third day, dead cells were detected in exfoliated lung cells.^ In alkaline elution studies, no DPC were detected in marrow of SD rats after 24 weeks exposure up to 15 ppm. During acute exposures, a dose response relationship was detected in SD rat exfoliated lung cells which yielded cross-linking factors of 0.954, 1.237, and 1.417 following a four hour exposure to 15, 50, or 150 ppm, respectively. No DPC were detected in the marrow at 150 ppm. In vitro exposures to HCHO of CHO and SHE cells and rat marrow cells revealed the production of DPC and DNA-DNA cross-links.^ Cytoxan treatment of SD rats was used to provide positive controls for flow cytometry and alkaline elution. A drastic reduction in RNA content and cycling cells occurred one day following treatment. After four days, RNA content was greatly increased; and on day eleven the marrow had regenerated. DPCs were detected in both the marrow and the exfoliated lung cells.^ The lack of significant responses in SD rats and A-strain mice below 15 ppm HCHO is explainable by host defense mechanisms. Apparently, the mucociliary apparatus and enzymatic detoxification are sufficient to reduce systemic toxicity to low level concentrations of formaldehyde. ^

Motivational interviewing and feedback for college drinkers: Handling missing values with complete case analysis and multiple imputation

Objective: In this secondary data analysis, three statistical methodologies were implemented to handle cases with missing data in a motivational interviewing and feedback study. The aim was to evaluate the impact that these methodologies have on the data analysis. ^ Methods: We first evaluated whether the assumption of missing completely at random held for this study. We then proceeded to conduct a secondary data analysis using a mixed linear model to handle missing data with three methodologies (a) complete case analysis, (b) multiple imputation with explicit model containing outcome variables, time, and the interaction of time and treatment, and (c) multiple imputation with explicit model containing outcome variables, time, the interaction of time and treatment, and additional covariates (e.g., age, gender, smoke, years in school, marital status, housing, race/ethnicity, and if participants play on athletic team). Several comparisons were conducted including the following ones: 1) the motivation interviewing with feedback group (MIF) vs. the assessment only group (AO), the motivation interviewing group (MIO) vs. AO, and the intervention of the feedback only group (FBO) vs. AO, 2) MIF vs. FBO, and 3) MIF vs. MIO.^ Results: We first evaluated the patterns of missingness in this study, which indicated that about 13% of participants showed monotone missing patterns, and about 3.5% showed non-monotone missing patterns. Then we evaluated the assumption of missing completely at random by Little's missing completely at random (MCAR) test, in which the Chi-Square test statistic was 167.8 with 125 degrees of freedom, and its associated p-value was p=0.006, which indicated that the data could not be assumed to be missing completely at random. After that, we compared if the three different strategies reached the same results. For the comparison between MIF and AO as well as the comparison between MIF and FBO, only the multiple imputation with additional covariates by uncongenial and congenial models reached different results. For the comparison between MIF and MIO, all the methodologies for handling missing values obtained different results. ^ Discussions: The study indicated that, first, missingness was crucial in this study. Second, to understand the assumptions of the model was important since we could not identify if the data were missing at random or missing not at random. Therefore, future researches should focus on exploring more sensitivity analyses under missing not at random assumption.^