Linearly scaling direct method for accurately inverting sparse banded matrices

In many problems in Computational Physics and Chemistry, one finds a special kind of sparse matrices, termed "banded matrices". These matrices, which are defined as having non-zero entries only within a given distance from the main diagonal, need often to be inverted in order to solve the associated linear system of equations. In this work, we introduce a new O(n) algorithm for solving such a system, being n X n the size of the matrix. We produce the analytical recursive expressions that allow to directly obtain the solution, as well as the pseudocode for its computer implementation. Moreover, we review the different options for possibly parallelizing the method, we describe the extension to deal with matrices that are banded plus a small number of non-zero entries outside the band, and we use the same ideas to produce a method for obtaining the full inverse matrix. Finally, we show that the New Algorithm is competitive, both in accuracy and in numerical efficiency, when compared to a standard method based in Gaussian elimination. We do this using sets of large random banded matrices, as well as the ones that appear when one tries to solve the 1D Poisson equation by finite differences.Comment: 24 pages, 5 figures, submitted to J. Comp. Phy

Detection of cytomegalovirus in bronchoalveolar lavage fluid from immunocompromised patients with pneumonitis by viral culture and DNA quantification

Altres ajuts: acords transformatius de la UABTo compare the detection of human cytomegalovirus (HCMV) in bronchoalveolar lavage (BAL) fluid by viral culture and quantitative polymerase chain reaction (qPCR), and to establish a viral load threshold that can identify cases of HCMV replication indicative of pneumonitis. There is currently no universal viral load cut-off to differentiate between patients with and without pneumonitis, and the interpretation of qPCR results is challenging. 176 consecutive BAL samples from immunosuppressed hosts with signs and/or symptoms of respiratory infection were prospectively studied by viral culture and qPCR. Concordant results were obtained in 81.25% of the BAL samples. The rest were discordant, as only 34% of the qPCR-positive BAL samples were positive by culture. The median HCMV load was significantly higher in culture-positive than in culture-negative BAL samples (5038 vs 178 IU/mL). Using a cut-off value of 1258 IU/mL of HCMV in BAL, pneumonia was diagnosed with a sensitivity of 76%, a specificity of 100%, a VPP of 100% and VPN of 98%, and HCMV was isolated in 100% of the BAL cultures. We found that a qPCR-negative was a quick and reliable way of ruling out HCMV pneumonitis, but a positive result did not always indicate clinically significant replication in the lung. However, an HCMV load in BAL fluid of ≥ 1258 IU/mL was always associated with disease, whereas < 200 IU/mL rarely so

An exact expression to calculate the derivatives of position-dependent observables in molecular simulations with flexible constraints

In this work, we introduce an algorithm to compute the derivatives of physical observables along the constrained subspace when flexible constraints are imposed on the system (i.e., constraints in which the hard coordinates are fixed to configuration-dependent values). The presented scheme is exact, it does not contain any tunable parameter, and it only requires the calculation and inversion of a sub-block of the Hessian matrix of second derivatives of the function through which the constraints are defined. We also present a practical application to the case in which the sought observables are the Euclidean coordinates of complex molecular systems, and the function whose minimization defines the constraints is the potential energy. Finally, and in order to validate the method, which, as far as we are aware, is the first of its kind in the literature, we compare it to the natural and straightforward finite-differences approach in three molecules of biological relevance: methanol, N-methyl-acetamide and a tri-glycine peptideComment: 13 pages, 8 figures, published versio

Advancing Key Gaps in the Knowledge of Plasmodium vivax Cryptic Infections Using Humanized Mouse Models and Organs-on-Chips

Plasmodium vivax is the most widely distributed human malaria parasite representing 36.3% of disease burden in the South-East Asia region and the most predominant species in the region of the Americas. Recent estimates indicate that 3.3 billion of people are under risk of infection with circa 7 million clinical cases reported each year. This burden is certainly underestimated as the vast majority of chronic infections are asymptomatic. For centuries, it has been widely accepted that the only source of cryptic parasites is the liver dormant stages known as hypnozoites. However, recent evidence indicates that niches outside the liver, in particular in the spleen and the bone marrow, can represent a major source of cryptic chronic erythrocytic infections. The origin of such chronic infections is highly controversial as many key knowledge gaps remain unanswered. Yet, as parasites in these niches seem to be sheltered from immune response and antimalarial drugs, research on this area should be reinforced if elimination of malaria is to be achieved. Due to ethical and technical considerations, working with the liver, bone marrow and spleen from natural infections is very difficult. Recent advances in the development of humanized mouse models and organs-on-a-chip models, offer novel technological frontiers to study human diseases, vaccine validation and drug discovery. Here, we review current data of these frontier technologies in malaria, highlighting major challenges ahead to study P. vivax cryptic niches, which perpetuate transmission and burden

The expression level of BAALC -associated microRNA miR-3151 is an independent prognostic factor in younger patients with cytogenetic intermediate-risk acute myeloid leukemia

Acute myeloid leukemia (AML) is a heterogeneous disease whose prognosis is mainly related to the biological risk conferred by cytogenetics and molecular profiling. In elderly patients (⩾60 years) with normal karyotype AML miR-3151 have been identified as a prognostic factor. However, miR-3151 prognostic value has not been examined in younger AML patients. In the present work, we have studied miR-3151 alone and in combination with BAALC, its host gene, in a cohort of 181 younger intermediate-risk AML (IR-AML) patients. Patients with higher expression of miR-3151 had shorter overall survival (P =0.0025), shorter leukemia-free survival (P =0.026) and higher cumulative incidence of relapse (P =0.082). Moreover, in the multivariate analysis miR-3151 emerged as independent prognostic marker in both the overall series and within the unfavorable molecular prognostic category. Interestingly, the combined determination of both miR-3151 and BAALC improved this prognostic stratification, with patients with low levels of both parameters showing a better outcome compared with those patients harboring increased levels of one or both markers (P =0.003). In addition, we studied the microRNA expression profile associated with miR-3151 identifying a six-microRNA signature. In conclusion, the analysis of miR-3151 and BAALC expression may well contribute to an improved prognostic stratification of younger patients with IR-AML

Metastatic gallbladder adenocarcinoma with signet-ring cells: A case report

<p>Abstract</p> <p>Introduction</p> <p>Signet-ring cell carcinoma is a rare and aggressive variant of mucinous adenocarcinoma. Only a few cases of gallbladder adenocarcinoma with signet-ring cells have been reported and because of this there is a lack of knowledge about the behavior and biology of this pathology.</p> <p>Case presentation</p> <p>We present the case of a 63-year-old Arab man with gallbladder signet-ring cell adenocarcinoma. He had an elective cholecystectomy and refused chemotherapy. Two months later, a small hepatic metastatic nodule was found, and nine months later he presented with multiple metastases in the liver, lymphatic nodes, both pleuras, peritoneum and subcutaneous tissue.</p> <p>Conclusion</p> <p>The proliferation of signet-ring cells in a gallbladder adenocarcinoma worsens the prognosis of an already adverse neoplasm. New lines of treatment in chemotherapy, such as cisplatin, or new biological therapy, such as monoclonal antibody c-myc oncogene, should be encouraged to improve the survival and life quality of these oncologic patients.</p

Human Gene Coexpression Landscape: Confident Network Derived from Tissue Transcriptomic Profiles

This is an open-access article distributed under the terms of the Creative Commons Attribution License.[Background]: Analysis of gene expression data using genome-wide microarrays is a technique often used in genomic studies to find coexpression patterns and locate groups of co-transcribed genes. However, most studies done at global >omic> scale are not focused on human samples and when they correspond to human very often include heterogeneous datasets, mixing normal with disease-altered samples. Moreover, the technical noise present in genome-wide expression microarrays is another well reported problem that many times is not addressed with robust statistical methods, and the estimation of errors in the data is not provided. [Methodology/Principal Findings]: Human genome-wide expression data from a controlled set of normal-healthy tissues is used to build a confident human gene coexpression network avoiding both pathological and technical noise. To achieve this we describe a new method that combines several statistical and computational strategies: robust normalization and expression signal calculation; correlation coefficients obtained by parametric and non-parametric methods; random cross-validations; and estimation of the statistical accuracy and coverage of the data. All these methods provide a series of coexpression datasets where the level of error is measured and can be tuned. To define the errors, the rates of true positives are calculated by assignment to biological pathways. The results provide a confident human gene coexpression network that includes 3327 gene-nodes and 15841 coexpression-links and a comparative analysis shows good improvement over previously published datasets. Further functional analysis of a subset core network, validated by two independent methods, shows coherent biological modules that share common transcription factors. The network reveals a map of coexpression clusters organized in well defined functional constellations. Two major regions in this network correspond to genes involved in nuclear and mitochondrial metabolism and investigations on their functional assignment indicate that more than 60% are house-keeping and essential genes. The network displays new non-described gene associations and it allows the placement in a functional context of some unknown non-assigned genes based on their interactions with known gene families. [Conclusions/Significance]: The identification of stable and reliable human gene to gene coexpression networks is essential to unravel the interactions and functional correlations between human genes at an omic scale. This work contributes to this aim, and we are making available for the scientific community the validated human gene coexpression networks obtained, to allow further analyses on the network or on some specific gene associations. The data are available free online at http://bioinfow.dep.usal.es/coexpression/. © 2008 Prieto et al.Funding and grant support was provided by the Ministery of Health, Spanish Government (ISCiii-FIS, MSyC; Project reference PI061153) and by the Ministery of Education, Castilla-Leon Local Government (JCyL; Project reference CSI03A06).Peer Reviewe

mRNA accumulation in the Cajal bodies of the diplotene larch microsporocyte

In microsporocytes of the European larch, we demonstrated the presence of several mRNAs in spherical nuclear bodies. In the nuclei of microsporocytes, we observed up to 12 bodies, ranging from 0.5 to 6 μm in diameter, during the prophase of the first meiotic division. Our previous studies revealed the presence of polyadenylated RNA (poly(A) RNA) in these bodies, but did not confirm the presence of nascent transcripts or splicing factors of the SR family. The lack of these molecules precludes the bodies from being the sites of synthesis and early maturation of primary transcripts (Kołowerzo et al., Protoplasma 236:13–19, 2009). However, the bodies serve as sites for the accumulation of splicing machinery, including the Sm proteins and small nuclear RNAs. Characteristic ultrastructures and the molecular composition of the nuclear bodies, which contain poly(A) RNA, are indicative of Cajal bodies (CBs). Here, we demonstrated the presence of several housekeeping gene transcripts—α-tubulin, pectin methylesterase, peroxidase and catalase, ATPase, and inositol-3-phosphate synthase—in CBs. Additionally, we observed transcripts of the RNA polymerase II subunits RPB2 and RPB10 RNA pol II and the core spliceosome proteins mRNA SmD1, SmD2, and SmE. The co-localization of nascent transcripts and mRNAs indicates that mRNA accumulation/storage, particularly in CBs, occurs in the nucleus of microsporocytes

Meiotic chromosomes and nucleolar behavior in testicular cells of the grassland spittlebugs Deois flavopicta, Mahanarva fimbriolata and Notozulia entreriana (Hemiptera, Auchenorrhyncha)

Spittlebugs annually infest pastures and cause severe damage, representing a serious problem for the tropical American beef cattle industry. Spittlebugs are an important biotic constraint to forage production and there is a lack of cytogenetic data for this group of insects. For these reasons, we conducted this work, in which the spermatogenesis and nucleolar behavior of Deois flavopicta, Mahanarva fimbriolata and Notozulia entreriana were studied. The males possessed testes in the shape of a “bunch of grapes”; a variable number of testicular lobes per individual and polyploid nuclei composed of several heteropycnotic bodies. A heteropycnotic area was located in the periphery of the nucleus (prophase I); the chiasmata were terminal or interstitial; metaphases I were circular or linear and anaphase showed late migration of the sex chromosome. The chromosome complement had 2n = 19, except for N. entreriana (2n = 15); the spermatids were round with heteropycnotic material in the center and elongated with conspicuos chromatin. The analysis of testes after silver nitrate staining showed polyploid nuclei with three large and three smaller nucleolar bodies. Early prophase cells had an intensely stained nucleolar body located close to the chromatin and another less evident body located away from the chromatin. The nucleolar bodies disintegrated during diplotene. Silver staining occurred in two autosomes, in terminal and subterminal locations, the latter probably corresponding to the nucleolus organizer regions (NORs). The spermatids were round with a round nucleolar body and silver staining was observed in the medial and posterior region of the elongated part of the spermatid head