Indexing k-mers in linear space for quality value compression.

Many bioinformatics tools heavily rely on [Formula: see text]-mer dictionaries to describe the composition of sequences and allow for faster reference-free algorithms or look-ups. Unfortunately, naive [Formula: see text]-mer dictionaries are very memory-inefficient, requiring very large amount of storage space to save each [Formula: see text]-mer. This problem is generally worsened by the necessity of an index for fast queries. In this work, we discuss how to build an indexed linear reference containing a set of input [Formula: see text]-mers and its application to the compression of quality scores in FASTQ files. Most of the entropies of sequencing data lie in the quality scores, and thus they are difficult to compress. Here, we present an application to improve the compressibility of quality values while preserving the information for SNP calling. We show how a dictionary of significant [Formula: see text]-mers, obtained from SNP databases and multiple genomes, can be indexed in linear space and used to improve the compression of quality value. Availability: The software is freely available at https://github.com/yhhshb/yalff

Effects of local hypothermia-rewarming on physiology, metabolism and inflammation of acutely injured human spinal cord.

In five patients with acute, severe thoracic traumatic spinal cord injuries (TSCIs), American spinal injuries association Impairment Scale (AIS) grades A-C, we induced cord hypothermia (33 °C) then rewarming (37 °C). A pressure probe and a microdialysis catheter were placed intradurally at the injury site to monitor intraspinal pressure (ISP), spinal cord perfusion pressure (SCPP), tissue metabolism and inflammation. Cord hypothermia-rewarming, applied to awake patients, did not cause discomfort or neurological deterioration. Cooling did not affect cord physiology (ISP, SCPP), but markedly altered cord metabolism (increased glucose, lactate, lactate/pyruvate ratio (LPR), glutamate; decreased glycerol) and markedly reduced cord inflammation (reduced IL1β, IL8, MCP, MIP1α, MIP1β). Compared with pre-cooling baseline, rewarming was associated with significantly worse cord physiology (increased ICP, decreased SCPP), cord metabolism (increased lactate, LPR; decreased glucose, glycerol) and cord inflammation (increased IL1β, IL8, IL4, IL10, MCP, MIP1α). The study was terminated because three patients developed delayed wound infections. At 18-months, two patients improved and three stayed the same. We conclude that, after TSCI, hypothermia is potentially beneficial by reducing cord inflammation, though after rewarming these benefits are lost due to increases in cord swelling, ischemia and inflammation. We thus urge caution when using hypothermia-rewarming therapeutically in TSCI

Brachytherapy of stage II mobile tongue carcinoma. Prediction of local control and QOL

BACKGROUND: There is no consensus as to the prognostic model for brachytherapy of tongue carcinoma. This study was designed to evaluate the prognostic factors for local control based on a large population under a unified treatment policy. RESULTS: Between 1970 and 1998, 433 patients with stage II tongue squamous cell carcinoma were treated by low-dose-rate brachytherapy. This series included 277 patients treated with a linear source with a minimum follow-up of 3 years. A spacer was introduced in 1987. The primary local control rates were 85.6%. CONCLUSION: In the multivariate analysis, an invasive growth pattern was a significant factor for local recurrence. The disease-related survival was influenced by old age and an invasive growth pattern. A spacer lowered mandibular bone complications. The growth pattern was the most important factor for recurrence. Brachytherapy was associated with a high cure rate and the use of spacers brought about good quality of life (QOL)

Functional Relations in Stokes Multipliers and Solvable Models related to U_q(A^{(1)}_n)

Recently, Dorey and Tateo have investigated functional relations among Stokes multipliers for a Schr{\"o}dinger equation (second order differential equation) with a polynomial potential term in view of solvable models. Here we extend their studies to a restricted case of n+1-th order linear differential equations.Comment: 20 pages, some explanations improved, To appear in J. Phys.

Internal states of model isotropic granular packings. III. Elastic properties

In this third and final paper of a series, elastic properties of numerically simulated isotropic packings of spherical beads assembled by different procedures and subjected to a varying confining pressure P are investigated. In addition P, which determines the stiffness of contacts by Hertz's law, elastic moduli are chiefly sensitive to the coordination number, the possible values of which are not necessarily correlated with the density. Comparisons of numerical and experimental results for glass beads in the 10kPa-10MPa range reveal similar differences between dry samples compacted by vibrations and lubricated packings. The greater stiffness of the latter, in spite of their lower density, can hence be attributed to a larger coordination number. Voigt and Reuss bounds bracket bulk modulus B accurately, but simple estimation schemes fail for shear modulus G, especially in poorly coordinated configurations under low P. Tenuous, fragile networks respond differently to changes in load direction, as compared to load intensity. The shear modulus, in poorly coordinated packings, tends to vary proportionally to the degree of force indeterminacy per unit volume. The elastic range extends to small strain intervals, in agreement with experimental observations. The origins of nonelastic response are discussed. We conclude that elastic moduli provide access to mechanically important information about coordination numbers, which escape direct measurement techniques, and indicate further perspectives.Comment: Published in Physical Review E 25 page

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

新型コロナウイルス感染者の隔離短縮は可能か？. 京都大学プレスリリース. 2022-08-24.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

Linear-time protein 3-D structure searching with insertions and deletions

<p>Abstract</p> <p>Background</p> <p>Two biomolecular 3-D structures are said to be similar if the RMSD (root mean square deviation) between the two molecules' sequences of 3-D coordinates is less than or equal to some given constant bound. Tools for searching for similar structures in biomolecular 3-D structure databases are becoming increasingly important in the structural biology of the post-genomic era.</p> <p>Results</p> <p>We consider an important, fundamental problem of reporting all substructures in a 3-D structure database of chain molecules (such as proteins) which are similar to a given query 3-D structure, with consideration of indels (<it>i.e.</it>, insertions and deletions). This problem has been believed to be very difficult but its exact computational complexity has not been known. In this paper, we first prove that the problem in unbounded dimensions is NP-hard. We then propose a new algorithm that dramatically improves the average-case time complexity of the problem in 3-D in case the number of indels <it>k </it>is bounded by a constant. Our algorithm solves the above problem for a query of size <it>m </it>and a database of size <it>N </it>in average-case <it>O</it>(<it>N</it>) time, whereas the time complexity of the previously best algorithm was <it>O</it>(<it>Nm</it>^<it>k</it>+1).</p> <p>Conclusions</p> <p>Our results show that although the problem of searching for similar structures in a database based on the RMSD measure with indels is NP-hard in the case of unbounded dimensions, it can be solved in 3-D by a simple average-case linear time algorithm when the number of indels is bounded by a constant.</p

A buccal mucosa carcinoma treated with high dose rate brachytherapy

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/135251/1/acm20008.pd