lmerSeq: an R package for analyzing transformed RNA-Seq data with linear mixed effects models.

BMC Bioinformatics

Center for Genes, Environment and Health, National Jewish Health, 1400 Jackson St, Denver, CO, 80206, USA.

Published: November 2022


Category Ranking

98%

Total Visits

921

Avg Visit Duration

2 minutes

Citations

20

Article Abstract

Background: Studies that utilize RNA Sequencing (RNA-Seq) in conjunction with designs that introduce dependence between observations (e.g. longitudinal sampling) require specialized analysis tools to accommodate this additional complexity. This R package contains a set of utilities to fit linear mixed effects models to transformed RNA-Seq counts that properly account for this dependence when performing statistical analyses.

Results: In a simulation study comparing lmerSeq and two existing methodologies that also work with transformed RNA-Seq counts, we found that lmerSeq was comprehensively better in terms of nominal error rate control and statistical power.

Conclusions: Existing R packages for analyzing transformed RNA-Seq data with linear mixed models are limited in the variance structures they allow and/or the transformation methods they support. The lmerSeq package offers more flexibility in both of these areas and gave substantially better results in our simulations.

Download full-text PDF

Source
http://www.ncbi.nlm.nih.gov/pmc/articles/PMC9670578PMC
http://dx.doi.org/10.1186/s12859-022-05019-9DOI Listing

Publication Analysis

Top Keywords

transformed rna-seq
16
linear mixed
12
lmerseq package
8
analyzing transformed
8
rna-seq data
8
data linear
8
mixed effects
8
effects models
8
rna-seq counts
8
rna-seq
5

Similar Publications

Background: Primary Membranous Nephropathy (PMN) is characterized by dysregulated immune responses, with B cells playing critical roles in disease pathogenesis. However, the immunopathogenic mechanisms underlying B cell involvement in PMN remain elusive.

Methods: We employed single-cell RNA sequencing on peripheral blood mononuclear cell samples (PBMC) obtained from 6 patients with PMN and 3 healthy controls (NC) to explore the transformation of B cells and their interaction with immune cells.

View Article and Find Full Text PDF

Transcription factors (TFs) regulate gene expression by interacting with DNA in a sequence-specific manner. High-throughput in vitro technologies, such as protein-binding microarrays and HT-SELEX (high-throughput systematic evolution of ligands by exponential enrichment), have revealed the DNA-binding specificities of hundreds of TFs. However, they have limited ability to reliably identify lower-affinity DNA binding sites, which are increasingly recognized as important for precise spatiotemporal control of gene expression.

View Article and Find Full Text PDF

Multiscale identification of DNASE1L3 as a key target in MASLD progression to hepatocellular carcinoma.

Pathol Res Pract

August 2025

Institute of Interdisciplinary Integrative Medicine Research, Shanghai University of Traditional Chinese Medicine, Shanghai 201203, China. Electronic address:

Background: The rising incidence of hepatocellular carcinoma (HCC) linked to metabolic dysfunction-associated steatotic liver disease (MASLD) underscores the need to identify key drivers of malignant transformation. This study aimed to discover critical genes governing MASLD-to-HCC progression for early intervention.

Methods: Bulk and single-cell transcriptomic data were analyzed using pseudotime algorithms (Mfuzz, Monocle3) to pinpoint progression-related genes.

View Article and Find Full Text PDF

Bulk RNA-seq deconvolution typically uses single-cell RNA-sequencing (scRNA-seq) references, but some cell types are only detectable through single-nucleus RNA sequencing (snRNA-seq). Because snRNA-seq captures nuclear, but not cytoplasmic, transcripts, direct use as a reference could reduce deconvolution accuracy. Here, we systematically benchmark strategies to integrate both modalities, focusing on transformations and gene-filtering approaches that harmonize snRNA-seq with scRNA-seq references.

View Article and Find Full Text PDF

Unlabelled: Cysteine proteases (CPs), a pivotal class of proteolytic enzymes ubiquitously distributed across plant genomes, play critical roles in plant development, senescence, and immune responses. However, systematic investigations of CPs in maize ( L.) remain limited.

View Article and Find Full Text PDF