Efficient computation of subset selection probabilities with application to Cox regression

An efficient method is presented to compute the probabilityof selection of a specified subset from the set of all subsetsof a fixed size where the subsets are taken from a populationwhose units have varying individual probabilities of selection.The problem is motivated by the computation of the exact marginallikelihood for the Cox proportional hazards model.     

The neuronal dynamics underlying cognitive flexibility in set shifting tasks

The ability to switch attention from one aspect of an object to another or in other words to switch the “attentional set” as investigated in tasks like the “Wisconsin Card Sorting Test” is commonly referred to as cognitive flexibility. In this work we present a biophysically detailed neurodynamical model which illustrates the neuronal base of the processes related to this cognitive flexibility. For this purpose we conducted behavioral experiments which allow the combined evaluation of different aspects of set shifting tasks: uninstructed set shifts as investigated in Wisconsin-like tasks, effects of stimulus congruency as investigated in Stroop-like tasks and the contribution of working memory as investigated in “Delayed-Match-to-Sample” tasks. The work describes how general experimental findings are usable to design the architecture of a biophysical detailed though minimalistic model with a high orientation on neurobiological findings and how, in turn, the simulations support experimental investigations. The resulting model is able to account for experimental and individual response times and error rates and enables the switch of attention as a system inherent model feature: The switching process suggested by the model is based on the memorization of the visual stimuli and does not require any synaptic learning. The operation of the model thus demonstrates with at least a high probability the neuronal dynamics underlying a key component of human behavior: the ability to adapt behavior according to context requirements—cognitive flexibility. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Action Editor: Peter Dayan     

Computational analysis reveals abundance of potential glycoproteins in Archaea, Bacteria and Eukarya

Glycosylation is the most common type of post-translational modification (PTM) and is known to affect protein stability, folding and activity. Inactivity of enzymes mediating glycosylation can result in serious disorders including colon cancer and brain disorders. Out of five main types of glycosylation, N-linked glycosylation is most abundant and characterized by the addition of a sugar group to an Asparagine residue at the N-X-S/T motif. Enzyme mediating such transfer is known as oligosaccharyl transferase (OST). It has been hypothesized before that a significant number of proteins serve as glycoproteins. In this study, we used programming implementations of Python to statistically quantify the representation of glycoproteins by scanning all the available proteome sequence data at ExPASy server for the presence of glycoproteins and also the enzyme which plays critical role in glycosylation i.e. OST. Our results suggest that more than 50% of the proteins carry N-X-S/T motif i.e. they could be potential glycoproteins. Furthermore, approximately 28-36% (1/3) of proteins possesses signature motifs which are characteristic features of enzyme OST. Quantifying this bias individually reveals that both the number of proteins tagged with N-X-S/T motif and the average number of motifs per protein is significantly higher in case of eukaryotes when compared to prokaryotes. In the light of these results we conclude that there is a significant bias in the representation of glycoproteins in the proteomes of all species and is manifested substantially in eukaryotes and claim for glycosylation to be the most common and ubiquitous PTM in cells, especially in eukaryotes.     

Use of in silico tools for classification of novel missense mutations identified in dystrophin gene in developing countries

DMD gene which is composed of 79 exons is the largest known gene located on X chromosome (Xp21). Point mutations in the dystrophin gene are responsible for 30–35% of cases with DMD/BMD. Mutation analysis of all the exons of the DMD gene is costly in developing countries, therefore, a few of the exons are selected to be analyzed routinely in clinical laboratories. In this study, direct sequencing was used for detection of point mutations in 10 exons of dystrophin gene in patients affected with DMD without detectable large rearrangements. Freely available programs were used to predict the damaging effects of the mutations. Point mutations were successfully detected in three patients. Three novel mutations, two missense mutations located on nonconservative domains and a single nucleotide deletion, were detected. Missense mutations were predicted to change splicing efficiency. Detection of point mutations by DNA analysis followed by prediction of the pathogenecity by using bioinformatic tool might be an asset to provide proper diagnosis or genetic counseling to patients and their family.     

Current and future applications of flow cytometry in aquatic microbiology

Flow cytometry has become a valuable tool in aquatic and environmental microbiology that combines direct and rapid assays to determine numbers, cell size distribution and additional biochemical and physiological characteristics of individual cells, revealing the heterogeneity present in a population or community. Flow cytometry exhibits three unique technical properties of high potential to study the microbiology of aquatic systems: (i) its tremendous velocity to obtain and process data; (ii) the sorting capacity of some cytometers, which allows the transfer of specific populations or even single cells to a determined location, thus allowing further physical, chemical, biological or molecular analysis; and (iii) high-speed multiparametric data acquisition and multivariate data analysis. Flow cytometry is now commonly used in aquatic microbiology, although the application of cell sorting to microbial ecology and quantification of heterotrophic nanoflagellates and viruses is still under development. The recent development of laser scanning cytometry also provides a new way to further analyse sorted cells or cells recovered on filter membranes or slides. The main infrastructure limitations of flow cytometry are: cost, need for skilled and well-trained operators, and adequate refrigeration systems for high-powered lasers and cell sorters. The selection and obtaining of the optimal fluorochromes, control microorganisms and validations for a specific application may sometimes be difficult to accomplish.     

Oligo-(R)-3-hydroxybutyrate modification of sorting signal enables pore formation by Escherichia coli OmpA

The outer membrane protein A (OmpA) of Escherichia coli is a well-known model for protein targeting and protein folding. Wild-type OmpA, isolated either from cytoplasmic inclusion bodies or from outer membranes, forms narrow pores of ∼ 80 pS in planar lipid bilayers at room temperature. The pores are well structured with narrow conductance range when OmpA is isolated using lithium dodecyl sulfate (LDS) or RapiGest surfactant but display irregular conductance when OmpA is isolated with urea or guanidine hydrochloride. Previous studies have shown that serine residues S163 and S167 of the sorting signal of OmpA (residues 163-169), i.e., the essential sequence for outer membrane incorporation, are covalently modified by oligomers of (R)-3-hydroxybutyrate (cOHB). Here we find that single-mutants S163 and S167 of OmpA, which still contain cOHB on one serine of the sorting signal, form narrow pores in planar lipid bilayers at room temperature with lower and more irregular conductance than wild-type OmpA, whereas double mutants S163:S167 and S163:V166 of OmpA, with no cOHB on the sorting signal, are unable to form stable pores in planar lipid bilayers. Our results indicate that modification of serines in the sorting signal of OmpA by cOHB in the cytoplasm enables OmpA to incorporate into lipid bilayers at room temperature as a narrow pore. They further suggest that cOHB modification may be an important factor in protein targeting and protein folding.     

BLOC-2, AP-3, and AP-1 proteins function in concert with Rab38 and Rab32 proteins to mediate protein trafficking to lysosome-related organelles

Lysosome-related organelles (LROs) are synthesized in specialized cell types where they largely coexist with conventional lysosomes. Most of the known cellular transport machinery involved in biogenesis are ubiquitously expressed and shared between lysosomes and LROs. Examples of common components are the adaptor protein complex-3 (AP-3) and biogenesis of lysosome-related organelle complex (BLOC)-2. These protein complexes control sorting and transport of newly synthesized integral membrane proteins from early endosomes to both lysosomes and LROs such as the melanosome. However, it is unknown what factors cooperate with the ubiquitous transport machinery to mediate transport to LROs in specialized cells. Focusing on the melanosome, we show that the ubiquitous machinery interacts with cell type-specific Rab proteins, Rab38 and Rab32, to facilitate transport to the maturing organelle. BLOC-2, AP-3, and AP-1 coimmunoprecipitated with Rab38 and Rab32 from MNT-1 melanocytic cell extracts. BLOC-2, AP-3, AP-1, and clathrin partially colocalized with Rab38 and Rab32 by confocal immunofluorescence microscopy in MNT-1 cells. Rab38- and Rab32-deficient MNT-1 cells displayed abnormal trafficking and steady state levels of known cargoes of the BLOC-2, AP-3, and AP-1 pathways, the melanin-synthesizing enzymes tyrosinase and tyrosinase-related protein-1. These observations support the idea that Rab38 and Rab32 are the specific factors that direct the ubiquitous machinery to mediate transport from early endosomes to maturing LROs. Additionally, analysis of tyrosinase-related protein-2 and total melanin production indicates that Rab32 has unique functions that cannot be carried out by Rab38 in melanosome biogenesis.     

单细胞转录组测序数据分析算法在干细胞研究中的应用

细胞的转录组决定其生理状态,每个细胞的转录组都是唯一的。借助单细胞转录组测序可分析单个干细胞的转录组特征,通过进一步的运算方法可以根据转录组特征对细胞进行细胞状态测定以及系谱分化特征的重建,在干细胞及组织发育研究中发挥了强大的作用,推动了其快速发展,加速了对干细胞分化及组织发育的相关过程及调控路径的认识。尤其是在干细胞领域的应用,得益于新算法的发展,单细胞转录组测序分析可用来阐述干细胞的起源、异质性,尤其是对干细胞的分化过程进行连续观察。本文主要对应用于干细胞分化相关研究的单细胞转录组测序数据新的算法及其应用进行了综述。     

超声图像处理中Snake模型研究

Snake模型是一种基于高层信息的有效目标轮廓提取算法,其优点是作用过程及最后结果的目标轮廓是一条完整的曲线,从而引起广泛的关注。鉴于医学超声图像的信噪比较低,用经典的边缘提取算法无法得到较好的结果,因此人们将Snake模型进行了各种各样的改进,并且越来越多地将它运用到医学超声图像处理中来。本文对乳腺超声图像进行阈值分割、形态滤波等一系列预处理后,将改进的Snake模型对乳腺超声图像进行肿瘤的边缘提取,得到了比较好的结果。