Joint tensor modeling of single cell 3D genome and epigenetic data with Muscle.

Emerging single cell technologies that simultaneously capture long-range interactions of genomic loci together with their DNA methylation levels are advancing our understanding of three-dimensional genome structure and its interplay with the epigenome at the single cell level. While methods to analyze data from single cell high throughput chromatin conformation capture (scHi-C) experiments are maturing, methods that can jointly analyze multiple single cell modalities with scHi-C data are lacking. Here, we introduce Muscle, a semi-nonnegative joint decomposition of Multiple single cell tensors, to jointly analyze 3D conformation and DNA methylation data at the single cell level.

Whole genome methylation sequencing in blood from persons with mild cognitive impairment and dementia due to Alzheimer's disease identifies cognitive status.

Introduction: Whole genome methylation sequencing (WGMS) in blood identifies differential DNA methylation in persons with late-onset dementia due to Alzheimer's disease (AD) but has not been tested in persons with mild cognitive impairment (MCI).

Methods: We used WGMS to compare DNA methylation levels at 25,244,219 CpG loci in 382 blood samples from 99 persons with MCI, 109 with AD, and 174 who are cognitively unimpaired (CU).

Results: WGMS identified 9756 differentially methylated positions (DMPs) in persons with MCI, including 1743 differentially methylated genes encoding proteins in biological pathways related to synapse organization, dendrite development, and ion transport.

Harnessing the regenerative potential of interleukin11 to enhance heart repair.

Balancing between regenerative processes and fibrosis is crucial for heart repair, yet strategies regulating this balance remain a barrier to developing therapies. The role of Interleukin 11 (IL11) in heart regeneration remains controversial, as both regenerative and fibrotic functions have been reported. We uncovered that il11a, an Il11 homolog in zebrafish, can trigger robust regenerative programs in zebrafish hearts, including cardiomyocytes proliferation and coronary expansion, even in the absence of injury.

Integrative annotation scores of variants for impact on RNA binding protein activities.

Motivation: The ENCODE project generated a large collection of eCLIP-seq RNA binding protein (RBP) profiling data with accompanying RNA-seq transcriptomes of shRNA knockdown of RBPs. These data could have utility in understanding the functional impact of genetic variants, however their potential has not been fully exploited. We implement INCA (Integrative annotation scores of variants for impact on RBP activities) as a multi-step genetic variant scoring approach that leverages the ENCODE RBP data together with ClinVar and integrates multiple computational approaches to aggregate evidence.

Harnessing the regenerative potential of to enhance heart repair.

Balancing between regenerative processes and fibrosis is crucial for heart repair, yet strategies regulating this balance remain a barrier to developing therapies. While Interleukin11 (IL11) is known as a fibrotic factor, its contribution to heart regeneration is poorly understood. We uncovered that homolog in zebrafish, can trigger robust regenerative programs in zebrafish hearts, including cardiomyocytes proliferation and coronary expansion, even in the absence of injury.

PUF partner interactions at a conserved interface shape the RNA-binding landscape and cell fate in Caenorhabditis elegans.

Protein-RNA regulatory networks underpin much of biology. C. elegans FBF-2, a PUF-RNA-binding protein, binds over 1,000 RNAs to govern stem cells and differentiation.

FlyVISTA, an Integrated Machine Learning Platform for Deep Phenotyping of Sleep in .

Animal behavior depends on internal state. While subtle movements can signify significant changes in internal state, computational methods for analyzing these "microbehaviors" are lacking. Here, we present FlyVISTA, a machine-learning platform to characterize microbehaviors in freely-moving flies, which we use to perform deep phenotyping of sleep.

Whole genome methylation sequencing in blood identifies extensive differential DNA methylation in late-onset dementia due to Alzheimer's disease.

Introduction: DNA microarray-based studies report differentially methylated positions (DMPs) in blood between late-onset dementia due to Alzheimer's disease (AD) and cognitively unimpaired individuals, but interrogate < 4% of the genome.

Methods: We used whole genome methylation sequencing (WGMS) to quantify DNA methylation levels at 25,409,826 CpG loci in 281 blood samples from 108 AD and 173 cognitively unimpaired individuals.

Results: WGMS identified 28,038 DMPs throughout the human methylome, including 2707 differentially methylated genes (e.

MuDCoD: multi-subject community detection in personalized dynamic gene networks from single-cell RNA sequencing.

Motivation: With the wide availability of single-cell RNA-seq (scRNA-seq) technology, population-scale scRNA-seq datasets across multiple individuals and time points are emerging. While the initial investigations of these datasets tend to focus on standard analysis of clustering and differential expression, leveraging the power of scRNA-seq data at the personalized dynamic gene co-expression network level has the potential to unlock subject and/or time-specific network-level variation, which is critical for understanding phenotypic differences. Community detection from co-expression networks of multiple time points or conditions has been well-studied; however, none of the existing settings included networks from multiple subjects and multiple time points simultaneously.

Debiased personalized gene coexpression networks for population-scale scRNA-seq data.

Population-scale single-cell RNA-seq (scRNA-seq) data sets create unique opportunities for quantifying expression variation across individuals at the gene coexpression network level. Estimation of coexpression networks is well established for bulk RNA-seq; however, single-cell measurements pose novel challenges owing to technical limitations and noise levels of this technology. Gene-gene correlation estimates from scRNA-seq tend to be severely biased toward zero for genes with low and sparse expression.

Secondary bile acids function through the vitamin D receptor in myeloid progenitors to promote myelopoiesis.

Metabolic products of the microbiota can alter hematopoiesis. However, the contribution and site of action of bile acids is poorly understood. Here, we demonstrate that the secondary bile acids, deoxycholic acid (DCA) and lithocholic acid (LCA), increase bone marrow myelopoiesis.

Dozer: Debiased personalized gene co-expression networks for population-scale scRNA-seq data.

Population-scale single cell RNA-seq (scRNA-seq) datasets create unique opportunities for quantifying expression variation across individuals at the gene co-expression network level. Estimation of co-expression networks is well-established for bulk RNA-seq; however, single-cell measurements pose novel challenges due to technical limitations and noise levels of this technology. Gene-gene correlation estimates from scRNA-seq tend to be severely biased towards zero for genes with low and sparse expression.

AdaLiftOver: high-resolution identification of orthologous regulatory elements with Adaptive liftOver.

Motivation: Elucidating functionally similar orthologous regulatory regions for human and model organism genomes is critical for exploiting model organism research and advancing our understanding of results from genome-wide association studies (GWAS). Sequence conservation is the de facto approach for finding orthologous non-coding regions between human and model organism genomes. However, existing methods for mapping non-coding genomic regions across species are challenged by the multi-mapping, low precision, and low mapping rate issues.

Pathogenic human variant that dislocates GATA2 zinc fingers disrupts hematopoietic gene expression and signaling networks.

Although certain human genetic variants are conspicuously loss of function, decoding the impact of many variants is challenging. Previously, we described a patient with leukemia predisposition syndrome (GATA2 deficiency) with a germline GATA2 variant that inserts 9 amino acids between the 2 zinc fingers (9aa-Ins). Here, we conducted mechanistic analyses using genomic technologies and a genetic rescue system with Gata2 enhancer-mutant hematopoietic progenitor cells to compare how GATA2 and 9aa-Ins function genome-wide.

Joint tensor modeling of single cell 3D genome and epigenetic data with Muscle.

Emerging single cell technologies that simultaneously capture long-range interactions of genomic loci together with their DNA methylation levels are advancing our understanding of three-dimensional genome structure and its interplay with the epigenome at the single cell level. While methods to analyze data from single cell high throughput chromatin conformation capture (scHi-C) experiments are maturing, methods that can jointly analyze multiple single cell modalities with scHi-C data are lacking. Here, we introduce Muscle, a semi-nonnegative joint decomposition of ltiple ingle el tnsors, to jointly analyze 3D conformation and DNA methylation data at the single cell level.

Normalization and de-noising of single-cell Hi-C data with BandNorm and scVI-3D.

Single-cell high-throughput chromatin conformation capture methodologies (scHi-C) enable profiling of long-range genomic interactions. However, data from these technologies are prone to technical noise and biases that hinder downstream analysis. We develop a normalization approach, BandNorm, and a deep generative modeling framework, scVI-3D, to account for scHi-C specific biases.

Regulatory Architecture of the RCA Gene Cluster Captures an Intragenic TAD Boundary, CTCF-Mediated Chromatin Looping and a Long-Range Intergenic Enhancer.

The Regulators of Complement Activation (RCA) gene cluster comprises several tandemly arranged genes with shared functions within the immune system. RCA members, such as complement receptor 2 (), are well-established susceptibility genes in complex autoimmune diseases. Altered expression of RCA genes has been demonstrated at both the functional and genetic level, but the mechanisms underlying their regulation are not fully characterised.

scGAD: single-cell gene associating domain scores for exploratory analysis of scHi-C data.

Summary: Quantitative tools are needed to leverage the unprecedented resolution of single-cell high-throughput chromatin conformation (scHi-C) data and integrate it with other single-cell data modalities. We present single-cell gene associating domain (scGAD) scores as a dimension reduction and exploratory analysis tool for scHi-C data. scGAD enables summarization at the gene unit while accounting for inherent gene-level genomic biases.

Gene by environment interaction mouse model reveals a functional role for 5-hydroxymethylcytosine in neurodevelopmental disorders.

Mouse knockouts of show altered neurodevelopmental behavior, deficits in striatal GABAergic signaling, and a genome-wide disruption of an environmentally sensitive DNA methylation modification (5-hydroxymethylcytosine [5hmC]) in the orthologs of a significant number of genes implicated in human neurodevelopmental disorders. We tested adult heterozygous mice ( ; lacking behavioral or neuropathological abnormalities) subjected to a prenatal stress and found that prenatally stressed female mice show repetitive behaviors and altered sociability, similar to the homozygote phenotype. Genomic profiling revealed disruptions in hippocampal and striatal 5hmC levels that are correlated to altered transcript levels of genes linked to these phenotypes (e.

MLG: multilayer graph clustering for multi-condition scRNA-seq data.

Single-cell transcriptome sequencing (scRNA-seq) enabled investigations of cellular heterogeneity at exceedingly higher resolutions. Identification of novel cell types or transient developmental stages across multiple experimental conditions is one of its key applications. Linear and non-linear dimensionality reduction for data integration became a foundational tool in inference from scRNA-seq data.

Identification of direct transcriptional targets of NFATC2 that promote β cell proliferation.

The transcription factor NFATC2 induces β cell proliferation in mouse and human islets. However, the genomic targets that mediate these effects have not been identified. We expressed active forms of Nfatc2 and Nfatc1 in human islets.

PKM2-TMEM33 axis regulates lipid homeostasis in cancer cells by controlling SCAP stability.

The pyruvate kinase M2 isoform (PKM2) is preferentially expressed in cancer cells to regulate anabolic metabolism. Although PKM2 was recently reported to regulate lipid homeostasis, the molecular mechanism remains unclear. Herein, we discovered an ER transmembrane protein 33 (TMEM33) as a downstream effector of PKM2 that regulates activation of SREBPs and lipid metabolism.