Robust high-dimensional precision matrix estimation

The dependency structure of multivariate data can be analyzed using the covariance matrix $\Sigma$. In many fields the precision matrix $\Sigma^{-1}$ is even more informative. As the sample covariance estimator is singular in high-dimensions, it cannot be used to obtain a precision matrix estimator. A popular high-dimensional estimator is the graphical lasso, but it lacks robustness. We consider the high-dimensional independent contamination model. Here, even a small percentage of contaminated cells in the data matrix may lead to a high percentage of contaminated rows. Downweighting entire observations, which is done by traditional robust procedures, would then results in a loss of information. In this paper, we formally prove that replacing the sample covariance matrix in the graphical lasso with an elementwise robust covariance matrix leads to an elementwise robust, sparse precision matrix estimator computable in high-dimensions. Examples of such elementwise robust covariance estimators are given. The final precision matrix estimator is positive definite, has a high breakdown point under elementwise contamination and can be computed fast

Comment on ``Passage Times for Unbiased Polymer Translocation through a Narrow Pore''

One of the most fundamental quantities associated with polymer translocation through a nanopore is the translocation time $\tau$ and its dependence on the chain length $N$. Our simulation results based on both the bond fluctuation Monte Carlo and Molecular Dynamics methods confirm the original prediction $\tau\sim N^{2\nu+1}$, which scales in the same manner as the Rouse relaxation time of the chain except for a larger prefactor, and invalidates other scaling claims.Comment: 1+pages, 1 Figure, Minor change

SUMOylation of AMPK alpha 1 by PIAS4 specifically regulates mTORC1 signalling

AMP-activated protein kinase (AMPK) inhibits several anabolic pathways such as fatty acid and protein synthesis, and identification of AMPK substrate specificity would be useful to understand its role in particular cellular processes and develop strategies to modulate AMPK activity in a substrate-specific manner. Here we show that SUMOylation of Z attenuates AMPK activation specifically towards mTORC1 signalling. SUMOylation is also important for rapid inactivation of AMPK, to allow prompt restoration of mTORC1 signalling. PIAS4 and its SUMO E3 ligase activity are specifically required for the AMPK alpha 1 SUMOylation and the inhibition of AMPK alpha 1 activity towards mTORC1 signalling. The activity of a SUMOylation-deficient AMPK alpha 1a mutant is higher than the wild type towards mTORC1 signalling when reconstituted in AMPKa-deficient cells. PIAS4 depletion reduced growth of breast cancer cells, specifically when combined with direct AMPK activator A769662, suggesting that inhibiting AMPK alpha 1 SUMOylation can be explored to modulate AMPK activation and thereby suppress cancer cell growth.Peer reviewe

Apatites in Gale Crater

ChemCam is an active remote sensing instrument suite that has operated successfully on MSL since landing Aug. 6th, 2012. It uses laser pulses to remove dust and to analyze rocks up to 7 m away. Laser-induced breakdown spectroscopy (LIBS) obtains emission spectra of materials ablated from the samples in electronically excited states. The intensities of the emission lines scale with the abundances of the related element. ChemCam is sensitive to most major rock-forming elements as well as to a set of minor and trace elements such as F, Cl, Li, P, Sr, Ba, and Rb. The measured chemical composition can then be used to infer the mineralogical composition of the ablated material. Here, we report a summary of inferred apatite detections along the MSL traverse at Gale Crater. We present the geologic settings of these findings and derive some interpretations about the formation conditions of apatite in time and space

Role of microstructure and surface defects on the dissolution kinetics of CeO2, a UO2 fuel analogue.

The release of radionuclides from spent fuel in a geological disposal facility is controlled by the surface mediated dissolution of UO2 in groundwater. In this study we investigate the influence of reactive surface sites on the dissolution of a synthesised CeO2 analogue for UO2 fuel. Dissolution was performed on: CeO2 annealed at high temperature, which eliminated intrinsic surface defects (point defects and dislocations); CeO2-x annealed in inert and reducing atmospheres to induce oxygen vacancy defects; and on crushed CeO2 particles of different size fractions. BET surface area measurements were used as an indicator of reactive surface site concentration. Cerium stoichiometry, determined using X-ray Photoelectron Spectroscopy (XPS) and supported by X-ray Diffraction (XRD) analysis, was used to determine oxygen vacancy concentration. Upon dissolution in nitric acid medium at 90°C, a quantifiable relationship was established between the concentration of high energy surface sites and CeO2 dissolution rate; the greater the proportion of intrinsic defects and oxygen vacancies, the higher the dissolution rate. Dissolution of oxygen vacancy-containing CeO2-x gave rise to rates that were an order of magnitude greater than for CeO2 with fewer oxygen vacancies. While enhanced solubility of Ce3+ influenced the dissolution, it was shown that replacement of vacancy sites by oxygen significantly affected the dissolution mechanism due to changes in the lattice volume and strain upon dissolution and concurrent grain boundary decohesion. These results highlight the significant influence of defect sites and grain boundaries on the dissolution kinetics of UO2 fuel analogues and reduce uncertainty in the long-term performance of spent fuel in geological disposal

Vibrational Spectra of a Mechanosensitive Channel

We report the simulated vibrational spectra of a mechanosensitive membrane channel in different gating states. Our results show that while linear absorption is insensitive to structural differences, linear dichroism and sum-frequency generation spectroscopies are sensitive to the orientation of the transmembrane helices, which is changing during the opening process. Linear dichroism cannot distinguish an intermediate structure from the closed structure, but sum-frequency generation can. In addition, we find that two-dimensional infrared spectroscopy can be used to distinguish all three investigated gating states of the mechanosensitive membrane channel.

Fluctuating lattice-Boltzmann model for complex fluids

We develop and test numerically a lattice-Boltzmann (LB) model for nonideal fluids that incorporates thermal fluctuations. The fluid model is a momentum-conserving thermostat, for which we demonstrate how the temperature can be made equal at all length scales present in the system by having noise both locally in the stress tensor and by shaking the whole system in accord with the local temperature. The validity of the model is extended to a broad range of sound velocities. Our model features a consistent coupling scheme between the fluid and solid molecular dynamics objects, allowing us to use the LB fluid as a heat bath for solutes evolving in time without external Langevin noise added to the solute. This property expands the applicability of LB models to dense, strongly correlated systems with thermal fluctuations and potentially nonideal equations of state. Tests on the fluid itself and on static and dynamic properties of a coarse-grained polymer chain under strong hydrodynamic interactions are used to benchmark the model. The model produces results for single-chain diffusion that are in quantitative agreement with theory.Peer reviewe

Genome-wide linkage analysis of 972 bipolar pedigrees using single-nucleotide polymorphisms.

Because of the high costs associated with ascertainment of families, most linkage studies of Bipolar I disorder (BPI) have used relatively small samples. Moreover, the genetic information content reported in most studies has been less than 0.6. Although microsatellite markers spaced every 10 cM typically extract most of the genetic information content for larger multiplex families, they can be less informative for smaller pedigrees especially for affected sib pair kindreds. For these reasons we collaborated to pool family resources and carried out higher density genotyping. Approximately 1100 pedigrees of European ancestry were initially selected for study and were genotyped by the Center for Inherited Disease Research using the Illumina Linkage Panel 12 set of 6090 single-nucleotide polymorphisms. Of the ~1100 families, 972 were informative for further analyses, and mean information content was 0.86 after pruning for linkage disequilibrium. The 972 kindreds include 2284 cases of BPI disorder, 498 individuals with bipolar II disorder (BPII) and 702 subjects with recurrent major depression. Three affection status models (ASMs) were considered: ASM1 (BPI and schizoaffective disorder, BP cases (SABP) only), ASM2 (ASM1 cases plus BPII) and ASM3 (ASM2 cases plus recurrent major depression). Both parametric and non-parametric linkage methods were carried out. The strongest findings occurred at 6q21 (non-parametric pairs LOD 3.4 for rs1046943 at 119 cM) and 9q21 (non-parametric pairs logarithm of odds (LOD) 3.4 for rs722642 at 78 cM) using only BPI and schizoaffective (SA), BP cases. Both results met genome-wide significant criteria, although neither was significant after correction for multiple analyses. We also inspected parametric scores for the larger multiplex families to identify possible rare susceptibility loci. In this analysis, we observed 59 parametric LODs of 2 or greater, many of which are likely to be close to maximum possible scores. Although some linkage findings may be false positives, the results could help prioritize the search for rare variants using whole exome or genome sequencing

MutSβ exceeds MutSα in dinucleotide loop repair

The target substrates of DNA mismatch recognising factors MutSalpha (MSH2+MSH6) and MutSbeta (MSH2+MSH3) have already been widely researched. However, the extent of their functional redundancy and clinical substance remains unclear. Mismatch repair (MMR)-deficient tumours are strongly associated with microsatellite instability (MSI) and the degree and type of MSI seem to be dependent on the MMR gene affected, and is linked to its substrate specificities. Deficiency in MSH2 and MSH6 is associated with both mononucleotide and dinucleotide repeat instability. Although no pathogenic MSH3 mutations have been reported, its deficiency is also suggested to cause low dinucleotide repeat instability