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: 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: 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: Residual resistivity ratio (RRR) and thermal conductivity of metal matrix in metal/superconductor composite wires are important parameters for designing superconducting magnets. However, the resistivity of silver in reacted Ag/Bi-2212 wires has yet to be determined over temperature range from 4.2 K to 80 K because Bi-2212 filaments have a critical transition temperature T c of ∼ 80 K, and because it is unknown whether the RRR of Ag/Bi-2212 degrades with Cu diffusing from Bi-2212 filaments into silver sheathes at elevated temperatures and to what degree it varies with heat treatment. We measured the resistivity of stand-alone Ag and AgMg (Ag-0.2 wt.% Mg) wires as well as the resistivity of Ag and Ag-0.2 wt.% Mg in Ag/Bi- 2212 round wires reacted in 1 bar oxygen at 890 °C for 1, 8, 24 and 48 hours and quickly cooled to room temperature. The heat treatment was designed to reduce the critical current I c of Bi- 2212 wires to nearly zero while allowing Cu lo...

Description: N 6 -methyladenosine (m 6 A) is a prevalent RNA methylation modification involved in the regulation of degradation, subcellular localization, splicing and local conformation changes of RNA transcripts. High-throughput experiments have demonstrated that only a small fraction of the m 6 A consensus motifs in mammalian transcriptomes are modified. Therefore, accurate identification of RNA m 6 A sites becomes emergently important. For the above purpose, here a computational predictor of mammalian m 6 A site named SRAMP is established. To depict the sequence context around m 6 A sites, SRAMP combines three random forest classifiers that exploit the positional nucleotide sequence pattern, the K-nearest neighbor information and the position-independent nucleotide pair spectrum features, respectively. SRAMP uses either genomic sequences or cDNA sequences as its input. With either kind of input sequence, SRAMP achieves competitive performance in both cross-validation tests and rigorous independent benchmarking tests. Analyses of the informative features and overrepresented rules extracted from the random forest classifiers demonstrate that nucleotide usage preferences at the distal positions, in addition to those at the proximal positions, contribute to the classification. As a public prediction server, SRAMP is freely available at http://www.cuilab.cn/sramp/ .

Description: It has been shown that a realistic level of magnetization of dense molecular cloud cores can suppress the formation of a rotationally supported disc (RSD) through catastrophic magnetic braking in the axisymmetric ideal MHD limit. In this study, we present conditions for the formation of RSDs through non-ideal MHD effects computed self-consistently from an equilibrium chemical network. We find that removing from the standard MRN distribution the large population of very small grains (VSGs) of ~ 10 Å to few 100 Å that dominate the coupling of the bulk neutral matter to the magnetic field increases the ambipolar diffusivity by ~ 1–2 orders of magnitude at densities below 10 10 /cm –3 . The enhanced ambipolar diffusion (AD) in the envelope reduces the amount of magnetic flux dragged by the collapse into the circumstellar disc-forming region. Therefore, magnetic braking is weakened and more angular momentum can be retained. With continuous high angular momentum inflow, RSDs of tens of au are able to form, survive, and even grow in size, depending on other parameters including cosmic ray ionization rate, magnetic field strength, and rotation speed. Some discs become self-gravitating and evolve into rings in our 2D (axisymmetric) simulations, which have the potential to fragment into (close) multiple systems in 3D. We conclude that disc formation in magnetized cores is highly sensitive to chemistry, especially to grain sizes. A moderate grain coagulation/growth to remove the large population of VSGs, either in the prestellar phase or during free-fall collapse, can greatly promote AD and help formation of tens of au RSDs.