ALBERT

Description: Motivation: RNA-seq technology has been widely adopted as an attractive alternative to microarray-based methods to study global gene expression. However, robust statistical tools to analyze these complex datasets are still lacking. By grouping genes with similar expression profiles across treatments, cluster analysis provides insight into gene functions and networks, and hence is an important technique for RNA-seq data analysis. Results: In this manuscript, we derive clustering algorithms based on appropriate probability models for RNA-seq data. An expectation-maximization algorithm and another two stochastic versions of expectation-maximization algorithms are described. In addition, a strategy for initialization based on likelihood is proposed to improve the clustering algorithms. Moreover, we present a model-based hybrid-hierarchical clustering method to generate a tree structure that allows visualization of relationships among clusters as well as flexibility of choosing the number of clusters. Results from both simulation studies and analysis of a maize RNA-seq dataset show that our proposed methods provide better clustering results than alternative methods such as the K-means algorithm and hierarchical clustering methods that are not based on probability models. Availability and implementation: An R package, MBCluster.Seq, has been developed to implement our proposed algorithms. This R package provides fast computation and is publicly available at http://www.r-project.org . Contact: sy@swufe.edu.cn ; pliu@iastate.edu Supplementary Information: Supplementary data are available at Bioinformatics online.

Description: Author(s): Q. Fu, K. Zhu, Y. R. Lu, M. J. Easton, S. L. Gao, Z. Wang, F. J. Jia, H. P. Li, P. P. Gan, and Y. He A deuteron radio-frequency quadrupole (RFQ) is being built by the RFQ group at Peking University. It is a very compact high-current RFQ, operating at 162.5 MHz in continuous-wave mode. By optimizing the beam dynamics design, our simulations reached 98% transmission efficiency for acceleration of the... [Phys. Rev. Accel. Beams 20, 120101] Published Fri Dec 22, 2017

Description: Using data taken as part of the Bluedisk project, we study the connection between neutral hydrogen (H i ) in the environment of spiral galaxies and that in the galaxies themselves. We measure the total H i mass present in the environment in a statistical way by studying the distribution of noise peaks in the H i data cubes obtained for 40 galaxies observed with Westerbork Synthesis Radio Telescope. We find that galaxies whose H i mass fraction is high relative to standard scaling relations have an excess H i mass in the surrounding environment as well. Gas in the environment consists of gas clumps which are individually below the detection limit of our H i data. These clumps may be hosted by small satellite galaxies and/or be the high-density peaks of a more diffuse gas distribution in the intergalactic medium. We interpret this result as an indication for a picture in which the ${\rm H}\,{\small {I}}$ -rich central galaxies accrete gas from an extended gas reservoir present in their environment.

Description: We report new imaging polarimetry observations of the Galactic compact H ii region K3-50 using CanariCam at the Gran Telescopio Canarias. We use a standard polarimetric analysis technique, first outlined by Aitken, to decompose the observed polarization images centred at 8.7, 10.3, and 12.5 μm into the emissive and absorptive components from silicate grains that are aligned with the local magnetic field. These components reveal the spatially resolved magnetic field structures across the mid-infrared emission area of K3-50. We examine these structures and show that they are consistent with previously observed features and physical models of K3-50, such as the molecular torus and the ionized outflow. We propose a 3D geometry for all the structures seen at different wavelengths. We also compute relevant physical quantities in order to estimate the associated magnetic field strengths that would be implied under various physical assumptions. We compare these results with magnetohydrodynamic simulations of protostar formation that predict the magnetic field strength and configuration. We find that the magnetic field may be dynamically important in the innermost 0.2 pc of the molecular torus, but that the torus is more likely to be rotationally supported against gravity outside this radius. Similarly, magnetic fields are unlikely to dominate the global physics of the ionized outflow, but they may be important in helping confine the flow near the cavity wall in some locations. Ours is the first application of the Aitken technique to spatially resolved magnetic field structures in multiple layers along the line of sight, effectively a method of ‘polarization tomography’.

Description: Author(s): Q. Zhang, P. Li, Y. Wen, C. Zhao, J. W. Zhang, A. Manchon, W. B. Mi, Y. Peng, and X. X. Zhang To understand the interfacial scattering effect on the anomalous Hall effect (AHE), we prepared multilayers of ( F e ( 36 / n ) nm / A u ( 12 / n ) nm ) n using an e-beam evaporator. This structure design allowed us to investigate the effect of interfacial scattering on the AHE, while keeping the samples' thickness an… [Phys. Rev. B 94, 024428] Published Fri Jul 22, 2016

Description: Multiple sclerosis (MS) is the most common autoimmune disease of the central nervous system (CNS). It is characterized by the infiltration of autoreactive immune cells into the CNS, which target the myelin sheath, leading to the loss of neuronal function. Although it is accepted that MS is a multifactorial disorder with both genetic and environmental factors influencing its development and course, the molecular pathogenesis of MS has not yet been fully elucidated. Here, we studied the longitudinal gene expression profiles of whole-blood RNA from a cohort of 195 MS patients and 66 healthy controls. We analyzed these transcriptomes at both the individual transcript and the biological pathway level. We found 62 transcripts to be significantly up-regulated in MS patients; the expression of 11 of these genes was counter-regulated by interferon treatment, suggesting partial restoration of a ‘healthy’ gene expression profile. Global pathway analyses linked the proteasome and Wnt signaling to MS disease processes. Since genotypes from a subset of individuals were available, we were able to identify expression quantitative trait loci (eQTL), a number of which involved two genes of the MS gene signature. However, all these eQTL were also present in healthy controls. This study highlights the challenge posed by analyzing transcripts from whole blood and how these can be mitigated by using large, well-characterized cohorts of patients with longitudinal follow-up and multi-modality measurements.

Description: Thousands of novel transcripts have been identified using deep transcriptome sequencing. This discovery of large and ‘hidden’ transcriptome rejuvenates the demand for methods that can rapidly distinguish between coding and noncoding RNA. Here, we present a novel alignment-free method, Coding Potential Assessment Tool (CPAT), which rapidly recognizes coding and noncoding transcripts from a large pool of candidates. To this end, CPAT uses a logistic regression model built with four sequence features: open reading frame size, open reading frame coverage, Fickett TESTCODE statistic and hexamer usage bias. CPAT software outperformed (sensitivity: 0.96, specificity: 0.97) other state-of-the-art alignment-based software such as Coding-Potential Calculator (sensitivity: 0.99, specificity: 0.74) and Phylo Codon Substitution Frequencies (sensitivity: 0.90, specificity: 0.63). In addition to high accuracy, CPAT is approximately four orders of magnitude faster than Coding-Potential Calculator and Phylo Codon Substitution Frequencies, enabling its users to process thousands of transcripts within seconds. The software accepts input sequences in either FASTA- or BED-formatted data files. We also developed a web interface for CPAT that allows users to submit sequences and receive the prediction results almost instantly.

Description: Dilated cardiomyopathy (DCM) is a leading cause of heart failure. In families with autosomal-dominant DCM, heterozygous missense mutations were identified in RNA-binding motif protein 20 ( RBM20 ), a spliceosome protein induced during early cardiogenesis. Dermal fibroblasts from two unrelated patients harboring an RBM20 R636S missense mutation were reprogrammed to human induced pluripotent stem cells (hiPSCs) and differentiated to beating cardiomyocytes (CMs). Stage-specific transcriptome profiling identified differentially expressed genes ranging from angiogenesis regulator to embryonic heart transcription factor as initial molecular aberrations. Furthermore, gene expression analysis for RBM20- dependent splice variants affected sarcomeric ( TTN and LDB3 ) and calcium (Ca 2+ ) handling ( CAMK2D and CACNA1C ) genes. Indeed, RBM20 hiPSC-CMs exhibited increased sarcomeric length ( RBM20 : 1.747 ± 0.238 µm versus control: 1.404 ± 0.194 µm; P 〈 0.0001) and decreased sarcomeric width ( RBM20 : 0.791 ± 0.609 µm versus control: 0.943 ± 0.166 µm; P 〈 0.0001). Additionally, CMs showed defective Ca 2+ handling machinery with prolonged Ca 2+ levels in the cytoplasm as measured by greater area under the curve ( RBM20 : 814.718 ± 94.343 AU versus control: 206.941 ± 22.417 AU; P 〈 0.05) and higher Ca 2+ spike amplitude ( RBM20 : 35.281 ± 4.060 AU versus control:18.484 ± 1.518 AU; P 〈 0.05). β-adrenergic stress induced with 10 µ m norepinephrine demonstrated increased susceptibility to sarcomeric disorganization ( RBM20 : 86 ± 10.5% versus control: 40 ± 7%; P 〈 0.001). This study features the first hiPSC model of RBM20 familial DCM. By monitoring human cardiac disease according to stage-specific cardiogenesis, this study demonstrates RBM20 familial DCM is a developmental disorder initiated by molecular defects that pattern maladaptive cellular mechanisms of pathological cardiac remodeling. Indeed, hiPSC-CMs recapitulate RBM20 familial DCM phenotype in a dish and establish a tool to dissect disease-relevant defects in RBM20 splicing as a global regulator of heart function.

Description: This paper is concerned with travelling wave solutions of periodic monostable differential equations on a two-dimensional lattice. The existence, stability and uniqueness of supercritical travelling wave solutions are obtained by comparison principle.

Description: Author(s): J. M. Yao (尧江明), E. F. Zhou (周恩付), and Z. P. Li (李志攀) We report the first beyond-mean-field study of low-lying parity-doublet states in Ra 224 by extending the multireference relativistic energy density functional method to include dynamical correlations related to symmetry restoration and quadrupole-octupole shape fluctuation with a generator coordinat… [Phys. Rev. C 92, 041304(R)] Published Tue Oct 20, 2015