Leucine-rich repeat containing 10 (LRRC10) is a cardiomyocyte-specific protein, but its role in cardiac biology is little understood. Recently Lrrc10 was identified as required for endogenous cardiac regeneration in zebrafish; however, whether LRRC10 plays a role in mammalian heart regeneration remains unclear. In this study, we demonstrate that Lrrc10-/- knockout mice exhibit a loss of the neonatal mouse regenerative response, marked by reduced cardiomyocyte cytokinesis and increased cardiomyocyte binucleation. Interestingly, LRRC10 deletion disrupts the regenerative transcriptional landscape of the regenerating neonatal mouse heart. Remarkably, cardiac overexpression of LRRC10 restores cardiomyocyte cytokinesis, increases cardiomyocyte mononucleation, and the cardiac regenerative capacity of Lrrc10-/- mice. Our results are consistent with a model in which LRRC10 is required for cardiomyocyte cytokinesis as well as regulation of the transcriptional landscape during mammalian heart regeneration.

No abstract available

Cell type-specific gene expression patterns are outputs of transcriptional gene regulatory networks (GRNs) that connect transcription factors and signaling proteins to target genes. Single-cell technologies such as single cell RNA-sequencing (scRNA-seq) and single cell Assay for Transposase-Accessible Chromatin using sequencing (scATAC-seq), can examine cell-type specific gene regulation at unprecedented detail. However, current approaches to infer cell type-specific GRNs are limited in their ability to integrate scRNA-seq and scATAC-seq measurements and to model network dynamics on a cell lineage. To address this challenge, we have developed single-cell Multi-Task Network Inference (scMTNI), a multi-task learning framework to infer the GRN for each cell type on a lineage from scRNA-seq and scATAC-seq data. Using simulated and real datasets, we show that scMTNI is a broadly applicable framework for linear and branching lineages that accurately infers GRN dynamics and identifies key regulators of fate transitions for diverse processes such as cellular reprogramming and differentiation.

Histone modifying enzymes play a central role in maintaining cell identity by establishing a conducive chromatin environment for lineage specific transcription factor activity. Pluripotent embryonic stem cell (ESC) identity is characterized by a lower abundance of gene repression associated histone modifications that enables rapid response to differentiation cues. The KDM3 family of histone demethylases removes the repressive histone H3 lysine 9 dimethylation (H3K9me2). Here we uncover a surprising role for the KDM3 proteins in the maintenance of the pluripotent state through post-transcriptional regulation. We find that KDM3A and KDM3B interact with RNA processing factors such as EFTUD2 and PRMT5. Acute selective degradation of the endogenous KDM3A and KDM3B proteins resulted in altered splicing independent of H3K9me2 status or catalytic activity. These splicing changes partially resemble the splicing pattern of the more blastocyst-like ground state of pluripotency and occurred in important chromatin and transcription factors such as Dnmt3b, Tbx3 and Tcf12. Our findings reveal non-canonical roles of histone demethylating enzymes in splicing to regulate cell identity.

Mechanisms of haematopoietic and cardiac patterning remain poorly understood. Here we show that the BMP and Wnt signalling pathways are integrated in an endoglin (Eng)-dependent manner in cardiac and haematopoietic lineage specification. Eng is expressed in early mesoderm and marks both haematopoietic and cardiac progenitors. In the absence of Eng, yolk sacs inappropriately express the cardiac marker, Nkx2.5. Conversely, high levels of Eng in vitro and in vivo increase haematopoiesis and inhibit cardiogenesis. Levels of Eng determine the activation of both BMP and Wnt pathways, which are integrated downstream of Eng by phosphorylation of Smad1 by Gsk3. By interrogating Eng-dependent Wnt-mediated transcriptional changes, we identify Jdp2 as a key Eng-dependent Wnt target, sufficient to establish haematopoietic fate in early mesoderm when BMP and Wnt crosstalk is disturbed. These studies provide mechanistic insight into the integration of BMP and Wnt signalling in the establishment of haematopoietic and cardiac progenitors during embryogenesis.

Regulation of gene expression is an important mechanism through which genetic variation can affect complex traits. A substantial portion of gene expression variation can be explained by both local (cis) and distal (trans) genetic variation. Much progress has been made in uncovering cis-acting expression quantitative trait loci (cis-eQTL), but trans-eQTL have been more difficult to identify and replicate. Here we take advantage of our ability to predict the cis component of gene expression coupled with gene mapping methods such as PrediXcan to identify high confidence candidate trans-acting genes and their targets. That is, we correlate the cis component of gene expression with observed expression of genes in different chromosomes. Leveraging the shared cis-acting regulation across tissues, we combine the evidence of association across all available Genotype-Tissue Expression Project tissues and find 2,356 trans-acting/target gene pairs with high mappability scores. Reassuringly, trans-acting genes are enriched in transcription and nucleic acid binding pathways and target genes are enriched in known transcription factor binding sites. Interestingly, trans-acting genes are more significantly associated with selected complex traits and diseases than target or background genes, consistent with percolating trans effects. Our scripts and summary statistics are publicly available for future studies of trans-acting gene regulation.

As yeast cultures enter stationary phase in rich, glucose-based medium, differentiation of two major subpopulations of cells, termed quiescent and nonquiescent, is observed. Differences in mRNA abundance between exponentially growing and stationary-phase cultures and quiescent and nonquiescent cells are known, but little was known about protein abundance in these cells. To measure protein abundance in exponential and stationary-phase cultures, the yeast GFP-fusion library (4159 strains) was examined during exponential and stationary phases, using high-throughput flow cytometry (HyperCyt). Approximately 5% of proteins in the library showed twofold or greater changes in median fluorescence intensity (abundance) between the two conditions. We examined 38 strains exhibiting two distinct fluorescence-intensity peaks in stationary phase and determined that the two fluorescence peaks distinguished quiescent and nonquiescent cells, the two major subpopulations of cells in stationary-phase cultures. GFP-fusion proteins in this group were more abundant in quiescent cells, and half were involved in mitochondrial function, consistent with the sixfold increase in respiration observed in quiescent cells and the relative absence of Cit1p:GFP in nonquiescent cells. Finally, examination of quiescent cell-specific GFP-fusion proteins revealed symmetry in protein accumulation in dividing quiescent and nonquiescent cells after glucose exhaustion, leading to a new model for the differentiation of these cells.

We mapped Polycomb-associated H3K27 trimethylation (H3K27me3) and Trithorax-associated H3K4 trimethylation (H3K4me3) across the whole genome in human embryonic stem (ES) cells. The vast majority of H3K27me3 colocalized on genes modified with H3K4me3. These commodified genes displayed low expression levels and were enriched in developmental function. Another significant set of genes lacked both modifications and was also expressed at low levels in ES cells but was enriched for gene function in physiological responses rather than development. Commodified genes could change expression levels rapidly during differentiation, but so could a substantial number of genes in other modification categories. SOX2, POU5F1, and NANOG, pluripotency-associated genes, shifted from modification by H3K4me3 alone to colocalization of both modifications as they were repressed during differentiation. Our results demonstrate that H3K27me3 modifications change during early differentiation, both relieving existing repressive domains and imparting new ones, and that colocalization with H3K4me3 is not restricted to pluripotent cells.

The majority of bacteria and archaea rely on CRISPR-Cas systems for RNA-guided, adaptive immunity against mobile genetic elements. The Cas9 family of type II CRISPR-associated DNA endonucleases generates programmable double strand breaks in the CRISPR-complementary DNA targets flanked by the PAM motif. Nowadays, CRISPR-Cas9 provides a set of powerful tools for precise genome manipulation in eukaryotes and prokaryotes. Recently, a few Cas9 orthologs have been reported to possess intrinsic CRISPR-guided, sequence-specific ribonuclease activities. These discoveries fundamentally expanded the targeting capability of CRISPR-Cas9 systems, and promise to provide new CRISPR tools to manipulate specific cellular RNA transcripts. Here we present a detailed method for the biochemical characterization of Cas9's RNA-targeting potential.

Genome-wide association studies (GWAS) have hitherto identified several germline variants associated with cancer susceptibility, but the molecular functions of these risk modulators remain largely uncharacterized. Recent studies have begun to uncover the regulatory potential of noncoding GWAS SNPs using epigenetic information in corresponding cancer cell types and matched normal tissues. However, this approach does not explore the potential effect of risk germline variants on other important cell types that constitute the microenvironment of tumor or its precursor. This paper presents evidence that the breast-cancer-associated variant rs3903072 may regulate the expression of CTSW in tumor-infiltrating lymphocytes. CTSW is a candidate tumor-suppressor gene, with expression highly specific to immune cells and also positively correlated with breast cancer patient survival. Integrative analyses suggest a putative causative variant in a GWAS-linked enhancer in lymphocytes that loops to the 3' end of CTSW through three-dimensional chromatin interaction. Our work thus poses the possibility that a cancer-associated genetic variant could regulate a gene not only in the cell of cancer origin but also in immune cells in the microenvironment, thereby modulating the immune surveillance by T lymphocytes and natural killer cells and affecting the clearing of early cancer initiating cells.

We performed RNA-seq and high-resolution mass spectrometry on 128 blood samples from COVID-19-positive and COVID-19-negative patients with diverse disease severities and outcomes. Quantified transcripts, proteins, metabolites, and lipids were associated with clinical outcomes in a curated relational database, uniquely enabling systems analysis and cross-ome correlations to molecules and patient prognoses. We mapped 219 molecular features with high significance to COVID-19 status and severity, many of which were involved in complement activation, dysregulated lipid transport, and neutrophil activation. We identified sets of covarying molecules, e.g., protein gelsolin and metabolite citrate or plasmalogens and apolipoproteins, offering pathophysiological insights and therapeutic suggestions. The observed dysregulation of platelet function, blood coagulation, acute phase response, and endotheliopathy further illuminated the unique COVID-19 phenotype. We present a web-based tool (covid-omics.app) enabling interactive exploration of our compendium and illustrate its utility through a machine learning approach for prediction of COVID-19 severity.

Recent advances in consortium-scale genome-wide association studies (GWAS) have highlighted the involvement of common genetic variants in autism spectrum disorder (ASD), but our understanding of their etiologic roles, especially the interplay with rare variants, is incomplete. In this work, we introduce an analytical framework to quantify the transmission disequilibrium of genetically regulated gene expression from parents to offspring. We applied this framework to conduct a transcriptome-wide association study (TWAS) on 7,805 ASD proband-parent trios, and replicated our findings using 35,740 independent samples. We identified 31 associations at the transcriptome-wide significance level. In particular, we identified POU3F2 (p = 2.1E-7), a transcription factor mainly expressed in developmental brain. Gene targets regulated by POU3F2 showed a 2.7-fold enrichment for known ASD genes (p = 2.0E-5) and a 2.7-fold enrichment for loss-of-function de novo mutations in ASD probands (p = 7.1E-5). These results provide a novel connection between rare and common variants, whereby ASD genes affected by very rare mutations are regulated by an unlinked transcription factor affected by common genetic variations.

High-throughput chromosome conformation capture assays, such as Hi-C, have shown that the genome is organized into organizational units such as topologically associating domains (TADs), which can impact gene regulatory processes. The sparsity of Hi-C matrices poses a challenge for reliable detection of these units. We present GRiNCH, a constrained matrix-factorization-based approach for simultaneous smoothing and discovery of TADs from sparse contact count matrices. GRiNCH shows superior performance against seven TAD-calling methods and three smoothing methods. GRiNCH is applicable to multiple platforms including SPRITE and HiChIP and can predict novel boundary factors with potential roles in genome organization.

Derivation and culture of small hepatocyte progenitor cells (SHPCs) capable of proliferating in vitro has been described in rodents and recently in humans. These cells are capable of engrafting in injured livers, however, they display de-differentiated morphology and reduced xenobiotic metabolism activity in culture over passages. Here we report that SHPCs derived from adult primary human hepatocytes (PHHs) and cultured on mouse embryonic fibroblasts (MEFs) not only display differentiated morphology and exhibit gene expression profiles similar to adult PHHs, but importantly, they retain their phenotype over several passages. Further, unlike previous reports, where extensive manipulations of culture conditions are required to convert SHPCs to metabolically functional hepatocytes, SHPCs in our co-culture system maintain expression of xenobiotic metabolism-associated genes. We show that SHPCs in co-culture are able to perform xenobiotic metabolism at rates equal to their parent PHHs as evidenced by the metabolism of acetaminophen to all of its major metabolites. In summary, we present an improved co-culture system that allows generation of SHPCs from adult PHHs that maintain their differentiated phenotype over multiple passages. Our findings would be useful for expansion of limited PHHs for use in studies of drug metabolism and toxicity testing.

During an immune response, macrophages systematically rewire their metabolism in specific ways to support their diversve functions. However, current knowledge of macrophage metabolism is largely concentrated on central carbon metabolism. Using multi-omics analysis, we identified nucleotide metabolism as one of the most significantly rewired pathways upon classical activation. Further isotopic tracing studies revealed several major changes underlying the substantial metabolomic alterations: 1) de novo synthesis of both purines and pyrimidines is shut down at several specific steps; 2) nucleotide degradation activity to nitrogenous bases is increased but complete oxidation of bases is reduced, causing a great accumulation of nucleosides and bases; and 3) cells gradually switch to primarily relying on salvaging the nucleosides and bases for maintaining most nucleotide pools. Mechanistically, the inhibition of purine nucleotide de novo synthesis is mainly caused by nitric oxide (NO)-driven inhibition of the IMP synthesis enzyme ATIC, with NO-independent transcriptional downregulation of purine synthesis genes augmenting the effect. The inhibition of pyrimidine nucleotide de novo synthesis is driven by NO-driven inhibition of CTP synthetase (CTPS) and transcriptional downregulation of thymidylate synthase (TYMS). For the rewiring of degradation, purine nucleoside phosphorylase (PNP) and uridine phosphorylase (UPP) are transcriptionally upregulated, increasing nucleoside degradation activity. However, complete degradation of purine bases by xanthine oxidoreductase (XOR) is inhibited by NO, diverting flux into nucleotide salvage. Inhibiting the activation-induced switch from nucleotide de novo synthesis to salvage by knocking out the purine salvage enzyme hypoxanthine-guanine phosporibosyl transferase (Hprt) significantly alters the expression of genes important for activated macrophage functions, suppresses macrophage migration, and increases pyroptosis. Furthermore, knocking out Hprt or Xor increases proliferation of the intracellular parasite Toxoplasma gondii in macrophages. Together, these studies comprehensively reveal the characteristics, the key regulatory mechanisms, and the functional importance of the dynamic rewiring of nucleotide metabolism in classically activated macrophages.

Sus scrofa domesticus (pig) has served as a superb large mammalian model for biomedical studies because of its comparable physiology and organ size to humans. The derivation of transgene-free porcine induced pluripotent stem cells (PiPSCs) will, therefore, benefit the development of porcine-specific models for regenerative biology and its medical applications. In the past, this effort has been hampered by a lack of understanding of the signaling milieu that stabilizes the porcine pluripotent state in vitro. Here, we report that transgene-free PiPSCs can be efficiently derived from porcine fibroblasts by episomal vectors along with microRNA-302/367 using optimized protocols tailored for this species. PiPSCs can be differentiated into derivatives representing the primary germ layers in vitro and can form teratomas in immunocompromised mice. Furthermore, the transgene-free PiPSCs preserve intrinsic species-specific developmental timing in culture, known as developmental allochrony. This is demonstrated by establishing a porcine in vitro segmentation clock model that, for the first time, displays a specific periodicity at ∼3.7 h, a timescale recapitulating in vivo porcine somitogenesis. We conclude that the transgene-free PiPSCs can serve as a powerful tool for modeling development and disease and developing transplantation strategies. We also anticipate that they will provide insights into conserved and unique features on the regulations of mammalian pluripotency and developmental timing mechanisms.

Cellular gene expression changes throughout a dynamic biological process, such as differentiation. Pseudotimes estimate cells' progress along a dynamic process based on their individual gene expression states. Ordering the expression data by pseudotime provides information about the underlying regulator-gene interactions. Because the pseudotime distribution is not uniform, many standard mathematical methods are inapplicable for analyzing the ordered gene expression states. Here we present single-cell inference of networks using Granger ensembles (SINGE), an algorithm for gene regulatory network inference from ordered single-cell gene expression data. SINGE uses kernel-based Granger causality regression to smooth irregular pseudotimes and missing expression values. It aggregates predictions from an ensemble of regression analyses to compile a ranked list of candidate interactions between transcriptional regulators and target genes. In two mouse embryonic stem cell differentiation datasets, SINGE outperforms other contemporary algorithms. However, a more detailed examination reveals caveats about poor performance for individual regulators and uninformative pseudotimes.

Our inability to derive the neuronal diversity that comprises the posterior central nervous system (pCNS) using human pluripotent stem cells (hPSCs) poses an impediment to understanding human neurodevelopment and disease in the hindbrain and spinal cord. Here, we establish a modular, monolayer differentiation paradigm that recapitulates both rostrocaudal (R/C) and dorsoventral (D/V) patterning, enabling derivation of diverse pCNS neurons with discrete regional specificity. First, neuromesodermal progenitors (NMPs) with discrete HOX profiles are converted to pCNS progenitors (pCNSPs). Then, by tuning D/V signaling, pCNSPs are directed to locomotor or somatosensory neurons. Expansive single-cell RNA-sequencing (scRNA-seq) analysis coupled with a novel computational pipeline allowed us to detect hundreds of transcriptional markers within region-specific phenotypes, enabling discovery of gene expression patterns across R/C and D/V developmental axes. These findings highlight the potential of these resources to advance a mechanistic understanding of pCNS development, enhance in vitro models, and inform therapeutic strategies.

Nitrogen is one of the most inaccessible plant nutrients, but certain species have overcome this limitation by establishing symbiotic interactions with nitrogen-fixing bacteria in the root nodule. This root-nodule symbiosis (RNS) is restricted to species within a single clade of angiosperms, suggesting a critical, but undetermined, evolutionary event at the base of this clade. To identify putative regulatory sequences implicated in the evolution of RNS, we evaluated the genomes of 25 species capable of nodulation and identified 3091 conserved noncoding sequences (CNS) in the nitrogen-fixing clade (NFC). We show that the chromatin accessibility of 452 CNS correlates significantly with the regulation of genes responding to lipochitooligosaccharides in Medicago truncatula. These included 38 CNS in proximity to 19 known genes involved in RNS. Five such regions are upstream of MtCRE1, Cytokinin Response Element 1, required to activate a suite of downstream transcription factors necessary for nodulation in M. truncatula. Genetic complementation of an Mtcre1 mutant showed a significant decrease of nodulation in the absence of the five CNS, when they are driving the expression of a functional copy of MtCRE1. CNS identified in the NFC may harbor elements required for the regulation of genes controlling RNS in M. truncatula.

Changes in the three-dimensional (3D) structure of the genome are an emerging hallmark of cancer. Cancer-associated copy number variants and single nucleotide polymorphisms promote rewiring of chromatin loops, disruption of topologically associating domains (TADs), active/inactive chromatin state switching, leading to oncogene expression and silencing of tumor suppressors. However, little is known about 3D changes during cancer progression to a chemotherapy-resistant state. We integrated chromatin conformation capture (Hi-C), RNA-seq, and whole-genome sequencing obtained from triple-negative breast cancer patient-derived xenograft primary tumors (UCD52) and carboplatin-resistant samples and found increased short-range (< 2 Mb) interactions, chromatin looping, formation of TAD, chromatin state switching into a more active state, and amplification of ATP-binding cassette transporters. Transcriptome changes suggested the role of long-noncoding RNAs in carboplatin resistance. Rewiring of the 3D genome was associated with TP53, TP63, BATF, FOS-JUN family of transcription factors and led to activation of aggressiveness-, metastasis- and other cancer-related pathways. Integrative analysis highlighted increased ribosome biogenesis and oxidative phosphorylation, suggesting the role of mitochondrial energy metabolism. Our results suggest that 3D genome remodeling may be a key mechanism underlying carboplatin resistance.