Dynamic regulatory module networks for inference of cell type-specific transcriptional networks.

Changes in transcriptional regulatory networks can significantly alter cell fate. To gain insight into transcriptional dynamics, several studies have profiled bulk multi-omic data sets with parallel transcriptomic and epigenomic measurements at different stages of a developmental process. However, integrating these data to infer cell type-specific regulatory networks is a major challenge. We present dynamic regulatory module networks (DRMNs), a novel approach to infer cell type-specific cis-regulatory networks and their dynamics. DRMN integrates expression, chromatin state, and accessibility to predict cis-regulators of context-specific expression, where context can be cell type, developmental stage, or time point, and uses multitask learning to capture network dynamics across linearly and hierarchically related contexts. We applied DRMNs to study regulatory network dynamics in three developmental processes, each showing different temporal relationships and measuring a different combination of regulatory genomic data sets: cellular reprogramming, liver dedifferentiation, and forward differentiation. DRMN identified known and novel regulators driving cell type-specific expression patterns, showing its broad applicability to examine dynamics of gene regulatory networks from linearly and hierarchically related multi-omic data sets.

Rapid changes in chromatin structure during dedifferentiation of primary hepatocytes in vitro.

Primary hepatocytes are widely used in the pharmaceutical industry to screen drug candidates for hepatotoxicity, but hepatocytes quickly dedifferentiate and lose their mature metabolic function in culture. Attempts have been made to better recapitulate the in vivo liver environment in culture, but the full spectrum of signals required to maintain hepatocyte function ex vivo remains elusive. To elucidate molecular changes that accompany, and may contribute to dedifferentiation of hepatocytes ex vivo, we performed lineage tracing and comprehensive profiling of alterations in their gene expression profiles and chromatin landscape during culture. First, using genetically tagged hepatocytes we demonstrate that expression of the fetal gene alpha-fetoprotein in cultured hepatocytes comes from cells that previously expressed the mature gene albumin, and not from a population of albumin-negative precursor cells, proving mature hepatocytes undergo true dedifferentiation in culture. Next we studied the dedifferentiation process in detail through bulk RNA-sequencing of hepatocytes cultured over an extended period. We identified three distinct phases of dedifferentiation: an early phase, where mature hepatocyte genes are rapidly downregulated in a matter of hours; a middle phase, where fetal genes are activated; and a late phase, where initially rare contaminating non-parenchymal cells proliferate, taking over the culture. Lastly, to better understand the signaling events that result in the rapid downregulation of mature genes in hepatocytes, we examined changes in chromatin accessibility in these cells during the first 24 h of culture using Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq). We find that drastic and rapid changes in chromatin accessibility occur immediately upon the start of culture. Using binding motif analysis of the areas of open chromatin sharing similar temporal profiles, we identify several candidate transcription factors potentially involved in the dedifferentiation of primary hepatocytes in culture.

Rational immunosilencing of a promiscuous T-cell epitope in the capsid of an adeno-associated virus.

Owing to the immunogenicity of adeno-associated viruses (AAVs), gene therapies using AAVs face considerable obstacles. Here, by leveraging ex vivo T-cell assays, the prediction of epitope binding to major histocompatibility complex class-II alleles, sequence-conservation analysis in AAV phylogeny and site-directed mutagenesis, we show that the replacement of amino acid residues in a promiscuous and most immunodominant T-cell epitope in the AAV9 capsid with AAV5 sequences abrogates the immune responses of peripheral blood mononuclear cells to the chimaeric vector while preserving its functions, potency, cellular specificity, transduction efficacy and biodistribution. This rational approach to the immunosilencing of capsid epitopes promiscuously binding to T cells may be applied to other AAV vectors and epitope regions.

Co-culture with mouse embryonic fibroblasts improves maintenance of metabolic function of human small hepatocyte progenitor cells.

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.

Reproducibility across single-cell RNA-seq protocols for spatial ordering analysis.

As newer single-cell protocols generate increasingly more cells at reduced sequencing depths, the value of a higher read depth may be overlooked. Using data from three different single-cell RNA-seq protocols that lend themselves to having either higher read depth (Smart-seq) or many cells (MARS-seq and 10X), we evaluate their ability to recapitulate biological signals in the context of spatial reconstruction. Overall, we find gene expression profiles after spatial reconstruction analysis are highly reproducible between datasets despite being generated by different protocols and using different computational algorithms. While UMI-based protocols such as 10X and MARS-seq allow for capturing more cells, Smart-seq's higher sensitivity and read-depth allow for analysis of lower expressed genes and isoforms. Additionally, we evaluate trade-offs for each protocol by performing subsampling analyses and find that optimizing the balance between sequencing depth and number of cells within a protocol is necessary for efficient use of resources. Our analysis emphasizes the importance of selecting a protocol based on the biological questions and features of interest.