Block CUR: Decomposing Matrices using Groups of Columns

A common problem in large-scale data analysis is to approximate a matrix using a combination of specifically sampled rows and columns, known as CUR decomposition. Unfortunately, in many real-world environments, the ability to sample specific individual rows or columns of the matrix is limited by either system constraints or cost. In this paper, we consider matrix approximation by sampling predefined \emph{blocks} of columns (or rows) from the matrix. We present an algorithm for sampling useful column blocks and provide novel guarantees for the quality of the approximation. This algorithm has application in problems as diverse as biometric data analysis to distributed computing. We demonstrate the effectiveness of the proposed algorithms for computing the Block CUR decomposition of large matrices in a distributed setting with multiple nodes in a compute cluster, where such blocks correspond to columns (or rows) of the matrix stored on the same node, which can be retrieved with much less overhead than retrieving individual columns stored across different nodes. In the biometric setting, the rows correspond to different users and columns correspond to users' biometric reaction to external stimuli, {\em e.g.,}~watching video content, at a particular time instant. There is significant cost in acquiring each user's reaction to lengthy content so we sample a few important scenes to approximate the biometric response. An individual time sample in this use case cannot be queried in isolation due to the lack of context that caused that biometric reaction. Instead, collections of time segments ({\em i.e.,} blocks) must be presented to the user. The practical application of these algorithms is shown via experimental results using real-world user biometric data from a content testing environment.Comment: shorter version to appear in ECML-PKDD 201

A Novel Unsupervised Method to Identify Genes Important in the Anti-viral Response: Application to Interferon/Ribavirin in Hepatitis C Patients

Background: Treating hepatitis C with interferon/ribavirin results in a varied response in terms of decrease in viral titer and ultimate outcome. Marked responders have a sharp decline in viral titer within a few days of treatment initiation, whereas in other patients there is no effect on the virus (poor responders). Previous studies have shown that combination therapy modifies expression of hundreds of genes in vitro and in vivo. However, identifying which, if any, of these genes have a role in viral clearance remains challenging. Aims: The goal of this paper is to link viral levels with gene expression and thereby identify genes that may be responsible for early decrease in viral titer. Methods: Microarrays were performed on RNA isolated from PBMC of patients undergoing interferon/ribavirin therapy. Samples were collected at pre-treatment (day 0), and 1, 2, 7, 14 and 28 days after initiating treatment. A novel method was applied to identify genes that are linked to a decrease in viral titer during interferon/ribavirin treatment. The method uses the relationship between inter-patient gene expression based proximities and inter-patient viral titer based proximities to define the association between microarray gene expression measurements of each gene and viral-titer measurements. Results: We detected 36 unique genes whose expressions provide a clustering of patients that resembles viral titer based clustering of patients. These genes include IRF7, MX1, OASL and OAS2, viperin and many ISG's of unknown function. Conclusion: The genes identified by this method appear to play a major role in the reduction of hepatitis C virus during the early phase of treatment. The method has broad utility and can be used to analyze response to any group of factors influencing biological outcome such as antiviral drugs or anti-cancer agents where microarray data are available. © 2007 Brodsky et al

Radial Growth of Qilian Juniper on the Northeast Tibetan Plateau and Potential Climate Associations

There is controversy regarding the limiting climatic factor for tree radial growth at the alpine treeline on the northeastern Tibetan Plateau. In this study, we collected 594 increment cores from 331 trees, grouped within four altitude belts spanning the range 3550 to 4020 m.a.s.l. on a single hillside. We have developed four equivalent ring-width chronologies and shown that there are no significant differences in their growth-climate responses during 1956 to 2011 or in their longer-term growth patterns during the period AD 1110–2011. The main climate influence on radial growth is shown to be precipitation variability. Missing ring analysis shows that tree radial growth at the uppermost treeline location is more sensitive to climate variation than that at other elevations, and poor tree radial growth is particularly linked to the occurrence of serious drought events. Hence water limitation, rather than temperature stress, plays the pivotal role in controlling the radial growth of Sabina przewalskii Kom. at the treeline in this region. This finding contradicts any generalisation that tree-ring chronologies from high-elevation treeline environments are mostly indicators of temperature changes

Engineering a BCR-ABL-activated caspase for the selective elimination of leukemic cells.

Increased understanding of the precise molecular mechanisms involved in cell survival and cell death signaling pathways offers the promise of harnessing these molecules to eliminate cancer cells without damaging normal cells. Tyrosine kinase oncoproteins promote the genesis of leukemias through both increased cell proliferation and inhibition of apoptotic cell death. Although tyrosine kinase inhibitors, such as the BCR-ABL inhibitor imatinib, have demonstrated remarkable efficacy in the clinic, drug-resistant leukemias emerge in some patients because of either the acquisition of point mutations or amplification of the tyrosine kinase, resulting in a poor long-term prognosis. Here, we exploit the molecular mechanisms of caspase activation and tyrosine kinase/adaptor protein signaling to forge a unique approach for selectively killing leukemic cells through the forcible induction of apoptosis. We have engineered caspase variants that can directly be activated in response to BCR-ABL. Because we harness, rather than inhibit, the activity of leukemogenic kinases to kill transformed cells, this approach selectively eliminates leukemic cells regardless of drug-resistant mutations

Liver-Targeting of Interferon-Alpha with Tissue-Specific Domain Antibodies

PMCID: PMC3581439This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Broadband gradient impedance matching using an acoustic metamaterial for ultrasonic transducers

2016-2017 > Academic research: refereed > Publication in refereed journal201804_a bcmaVersion of RecordPublishe

Assisted evolution enables HIV-1 to overcome a high trim5α-imposed genetic barrier to rhesus macaque tropism

Diversification of antiretroviral factors during host evolution has erected formidable barriers to cross-species retrovirus transmission. This phenomenon likely protects humans from infection by many modern retroviruses, but it has also impaired the development of primate models of HIV-1 infection. Indeed, rhesus macaques are resistant to HIV-1, in part due to restriction imposed by the TRIM5α protein (rhTRIM5α). Initially, we attempted to derive rhTRIM5α-resistant HIV-1 strains using two strategies. First, HIV-1 was passaged in engineered human cells expressing rhTRIM5α. Second, a library of randomly mutagenized capsid protein (CA) sequences was screened for mutations that reduced rhTRIM5α sensitivity. Both approaches identified several individual mutations in CA that reduced rhTRIM5α sensitivity. However, neither approach yielded mutants that were fully resistant, perhaps because the locations of the mutations suggested that TRIM5α recognizes multiple determinants on the capsid surface. Moreover, even though additive effects of various CA mutations on HIV-1 resistance to rhTRIM5α were observed, combinations that gave full resistance were highly detrimental to fitness. Therefore, we employed an 'assisted evolution' approach in which individual CA mutations that reduced rhTRIM5α sensitivity without fitness penalties were randomly assorted in a library of viral clones containing synthetic CA sequences. Subsequent passage of the viral library in rhTRIM5α-expressing cells resulted in the selection of individual viral species that were fully fit and resistant to rhTRIM5α. These viruses encoded combinations of five mutations in CA that conferred complete or near complete resistance to the disruptive effects of rhTRIM5α on incoming viral cores, by abolishing recognition of the viral capsid. Importantly, HIV-1 variants encoding these CA substitutions and SIVmac239 Vif replicated efficiently in primary rhesus macaque lymphocytes. These findings demonstrate that rhTRIM5α is difficult to but not impossible to evade, and doing so should facilitate the development of primate models of HIV-1 infection

TEMPRANILLO is a regulator of juvenility in plants

Many plants are incapable of flowering in inductive daylengths during the early juvenile vegetative phase (JVP). Arabidopsis mutants with reduced expression of TEMPRANILLO (TEM), a repressor of FLOWERING LOCUS T (FT) had a shorter JVP than wild-type plants. Reciprocal changes in mRNA expression of TEM and FT were observed in both Arabidopsis and antirrhinum, which correlated with the length of the JVP. FT expression was induced just prior to the end of the JVP and levels of TEM1 mRNA declined rapidly at the time when FT mRNA levels were shown to increase. TEM orthologs were isolated from antirrhinum (AmTEM) and olive (OeTEM) and were expressed most highly during their juvenile phase. AmTEM functionally complemented AtTEM1 in the tem1 mutant and over-expression of AmTEM prolonged the JVP through repression of FT and CONSTANS (CO). We propose that TEM may have a general role in regulating JVP in herbaceous and woody species

Spin-Imbalance in a One-Dimensional Fermi Gas

Superconductivity and magnetism generally do not coexist. Changing the relative number of up and down spin electrons disrupts the basic mechanism of superconductivity, where atoms of opposite momentum and spin form Cooper pairs. Nearly forty years ago Fulde and Ferrell and Larkin and Ovchinnikov proposed an exotic pairing mechanism (FFLO) where magnetism is accommodated by formation of pairs with finite momentum. Despite intense theoretical and experimental efforts, however, polarized superconductivity remains largely elusive. Here we report experimental measurements of density profiles of a two spin mixture of ultracold 6Li atoms trapped in an array of one dimensional (1D) tubes, a system analogous to electrons in 1D wires. At finite spin imbalance, the system phase separates with an inverted phase profile in comparison to the three-dimensional case. In 1D we find a partially polarized core surrounded by wings composed of either a completely paired BCS superfluid or a fully polarized Fermi gas, depending on the degree of polarization. Our observations are in quantitative agreement with theoretical calculations in which the partially polarized phase is found to be a 1D analogue of the FFLO state. This study demonstrates how ultracold atomic gases in 1D may be used to create non-trivial new phases of matter, and also paves the way for direct observation and further study of the FFLO phase.Comment: 30 pages, 7 figure

Selective Constraints on Amino Acids Estimated by a Mechanistic Codon Substitution Model with Multiple Nucleotide Changes

Empirical substitution matrices represent the average tendencies of substitutions over various protein families by sacrificing gene-level resolution. We develop a codon-based model, in which mutational tendencies of codon, a genetic code, and the strength of selective constraints against amino acid replacements can be tailored to a given gene. First, selective constraints averaged over proteins are estimated by maximizing the likelihood of each 1-PAM matrix of empirical amino acid (JTT, WAG, and LG) and codon (KHG) substitution matrices. Then, selective constraints specific to given proteins are approximated as a linear function of those estimated from the empirical substitution matrices. Akaike information criterion (AIC) values indicate that a model allowing multiple nucleotide changes fits the empirical substitution matrices significantly better. Also, the ML estimates of transition-transversion bias obtained from these empirical matrices are not so large as previously estimated. The selective constraints are characteristic of proteins rather than species. However, their relative strengths among amino acid pairs can be approximated not to depend very much on protein families but amino acid pairs, because the present model, in which selective constraints are approximated to be a linear function of those estimated from the JTT/WAG/LG/KHG matrices, can provide a good fit to other empirical substitution matrices including cpREV for chloroplast proteins and mtREV for vertebrate mitochondrial proteins. The present codon-based model with the ML estimates of selective constraints and with adjustable mutation rates of nucleotide would be useful as a simple substitution model in ML and Bayesian inferences of molecular phylogenetic trees, and enables us to obtain biologically meaningful information at both nucleotide and amino acid levels from codon and protein sequences.Comment: Table 9 in this article includes corrections for errata in the Table 9 published in 10.1371/journal.pone.0017244. Supporting information is attached at the end of the article, and a computer-readable dataset of the ML estimates of selective constraints is available from 10.1371/journal.pone.001724