Alternative polyadenylation (APA) is a key post-transcriptional regulatory mechanism implicated in various diseases. Existing APA analysis tools are generally restricted to site detection and comparison, precluding differently expressed gene (DEG) analysis. Furthermore, standard RNA-seq-based DEG methods, though commonly used for gene expression profiling, demonstrate limited efficacy in identifying DEGs when directly applied to APA-seq datasets. To address this limitation, we developed APAdeg, a novel statistical method specifically tailored for DEG analysis of APA-seq data. APAdeg integrates both the total read count of a gene and the site-specific read counts within the gene into a generalized linear mixed model, thereby improving the efficiency of DEG detection. Benchmarking analyses on both simulated and empirical APA-seq data demonstrated that APAdeg consistently outperforms RNA-seq-based methods in DEG inference. Application of APAdeg to APA-seq data from distinct cancer types revealed that only a small proportion of DEGs exhibited significant changes in 3' untranslated region length, with an equally small proportion showing significant alterations in intronic APA usage. To facilitate widespread adoption, we implemented APAdeg as an R package. Collectively, APAdeg significantly enhances the accuracy of DEG analysis from APA-seq data, thereby advancing research into APA-mediated gene regulation in diseases and health.

Deciphering the free energy landscape of biomolecular structure space is crucial for understanding many complex molecular processes, such as protein-protein interaction, RNA folding, and protein folding. A major source of current dynamic structure data is Molecular Dynamics (MD) simulations. Several methods have been proposed to investigate the free energy landscape from MD data, but all of them rely on the assumption that kinetic similarity is associated with global geometric similarity, which may lead to unsatisfactory results. In this paper, we proposed a new method called Conditional Angle Partition Tree to reveal the hierarchical free energy landscape by correlating local geometric similarity with kinetic similarity. Its application on the benchmark alanine dipeptide MD data showed a much better performance than existing methods in exploring and understanding the free energy landscape. We also applied it to the MD data of Villin HP35. Our results are more reasonable on various aspects than those from other methods and very informative on the hierarchical structure of its energy landscape.

Exploring the relationship between factors of interest is a fundamental step for further analysis on various scientific problems such as understanding the genetic mechanism underlying specific disease, brain functional connectivity analysis. There are many methods proposed for association analysis and each has its own advantages, but none of them is suitable for all kinds of situations. This brings difficulties and confusions to practitioner on which one to use when facing a real problem. In this paper, we propose to combine power of different methods to detect associations in large data sets. It goes as combining the weaker to be stronger. Numerical results from simulation study and real data applications show that our new framework is powerful. Importantly, the framework can also be applied to other problems. Availability: The R script is available at https://jiangdata.github.io/resources/DM.zip.

Subcellular localization of microRNAs (miRNAs) is an important reflection of their biological functions. Considering the spatio-temporal specificity of miRNA subcellular localization, experimental detection techniques are expensive and time-consuming, which strongly motivates an efficient and economical computational method to predict miRNA subcellular localization. In this paper, we describe a computational framework, MiRLoc, to predict the subcellular localization of miRNAs. In contrast to existing methods, MiRLoc uses the functional similarity between miRNAs instead of sequence features and incorporates information about the subcellular localization of the corresponding target mRNAs. The results show that miRNA functional similarity data can be effectively used to predict miRNA subcellular localization, and that inclusion of subcellular localization information of target mRNAs greatly improves prediction performance.

N6-methyladenine (6mA) is an important DNA modification form associated with a wide range of biological processes. Identifying accurately 6mA sites on a genomic scale is crucial for under-standing of 6mA's biological functions. However, the existing experimental techniques for detecting 6mA sites are cost-ineffective, which implies the great need of developing new computational methods for this problem. In this paper, we developed, without requiring any prior knowledge of 6mA and manually crafted sequence features, a deep learning framework named Deep6mA to identify DNA 6mA sites, and its performance is superior to other DNA 6mA prediction tools. Specifically, the 5-fold cross-validation on a benchmark dataset of rice gives the sensitivity and specificity of Deep6mA as 92.96% and 95.06%, respectively, and the overall prediction accuracy is 94%. Importantly, we find that the sequences with 6mA sites share similar patterns across different species. The model trained with rice data predicts well the 6mA sites of other three species: Arabidopsis thaliana, Fragaria vesca and Rosa chinensis with a prediction accuracy over 90%. In addition, we find that (1) 6mA tends to occur at GAGG motifs, which means the sequence near the 6mA site may be conservative; (2) 6mA is enriched in the TATA box of the promoter, which may be the main source of its regulating downstream gene expression.