Search

Copyright © 2020 American Chemical Society. Host protein folding stress responses can play important roles in RNA virus replication and evolution. Prior work suggested a complicated interplay between the cytosolic proteostasis stress response, controlled by the transcriptional master regulator heat shock factor 1 (HSF1), and human immunodeficiency virus-1 (HIV-1). We sought to uncouple HSF1 transcription factor activity from cytotoxic proteostasis stress and thereby better elucidate the proposed role(s) of HSF1 in the HIV-1 lifecycle. To achieve this objective, we used chemical genetic, stress-independent control of HSF1 activity to establish whether and how HSF1 influences HIV-1 replication. Stress-independent HSF1 induction decreased both the total quantity and infectivity of HIV-1 virions. Moreover, HIV-1 was unable to escape HSF1-mediated restriction over the course of several serial passages. These results clarify the interplay between the host's heat shock response and HIV-1 infection and motivate continued investigation of chaperones as potential antiviral therapeutic targets.

Targeted mutagenesis mediated by nucleotide base deaminase–T7 RNA polymerase fusions has recently emerged as a novel and broadly useful strategy to power genetic diversification in the context of in vivo directed evolution campaigns. Here, we expand the utility of this approach by introducing a highly active adenosine deaminase–T7 RNA polymerase fusion protein (eMutaT7A→G), resulting in higher mutation frequencies to enable more rapid directed evolution. We also assess the benefits and potential downsides of using this more active mutator. We go on to show in Escherichia coli that adenosine deaminase-bearing mutators (MutaT7A→G or eMutaT7A→G) can be employed in tandem with a cytidine deaminase-bearing mutator (MutaT7C→T) to introduce all possible transition mutations simultaneously. We illustrate the efficacy of this in vivo mutagenesis approach by exploring mutational routes to antibacterial drug resistance. This work sets the stage for general application of optimized MutaT7 tools able to induce all types of transition mutations during in vivo directed evolution campaigns across diverse organisms.

SUMMARYMalaria-causing P. vivax parasites can linger in the human liver for weeks to years, and then reactivate to cause recurrent blood-stage infection. While an important target for malaria eradication, little is known about the molecular features of the replicative and non-replicative states of intracellular P. vivax parasites, or the human host-cell responses to them. Here, we leverage a bioengineered human microliver platform to culture Thai clinical isolates of P. vivax in primary human hepatocytes and conduct transcriptional profiling of infected cultures. By coupling enrichment strategies with bulk and single-cell analyses, we captured both parasite and host transcripts in individual hepatocytes throughout the infection course. We defined host- and state-dependent transcriptional signatures and identified previously unappreciated populations of replicative and non-replicative parasites, sharing features with sexual transmissive forms. We found that infection suppresses transcription of key hepatocyte function genes, and that P. vivax elicits an innate immune response that can be manipulated to control infection. Our work provides an extendible framework and resource for understanding host-parasite interactions and reveals new insights into the biology of malaria dormancy and transmission.

Hematopoiesis is a process of constitutive regeneration whereby hematopoietic stem and progenitor cells (HSPCs) replenish mature blood cells. During maturation and aging, HSPCs shift their output to support the demands of prenatal development and postnatal maturation both at homeostasis and in response to stress. How HSPC ontogeny changes throughout life is unknown; studies to date have largely focused on specific individual ages, particularly at single cell resolution. Here, we performed single cell RNA-seq of human HSPCs from early prenatal development into mature adulthood. We observed shifts in HSPC transcriptional states and differentiation trajectories over time. We identified age-specific gene expression patterns throughout human maturation and developed methods for identifying, prospectively purifying, and functionally validating age-specific HSC states. Together, our findings define the temporal maturation of human HSPCs and uncover principles applicable to age-biased blood diseases.SummarySingle cell RNA sequencing reveals that the mechanisms of human hematopoietic stem and progenitor cell (HSPC) fate commitment change over a lifetime from gestation to mature adulthood.

International audience; Malaria-causing Plasmodium vivax parasites can linger in the human liver for weeks to years and reactivate to cause recurrent blood-stage infection. Although they are an important target for malaria eradication, little is known about the molecular features of replicative and non-replicative intracellular liver-stage parasites and their host cell dependence. Here, we leverage a bioengineered human microliver platform to culture patient-derived P. vivax parasites for transcriptional profiling. Coupling enrichment strategies with bulk and single-cell analyses, we capture both parasite and host transcripts in individual hepatocytes throughout the course of infection. We define host- and state-dependent transcriptional signatures and identify unappreciated populations of replicative and non-replicative parasites that share features with sexual transmissive forms. We find that infection suppresses the transcription of key hepatocyte function genes and elicits an anti-parasite innate immune response. Our work provides a foundation for understanding host-parasite interactions and reveals insights into the biology of P. vivax dormancy and transmission.

The mammalian brain is complex, with multiple cell types performing a variety of diverse functions, but exactly how each cell type is affected in aging remains largely unknown. Here we performed a single-cell transcriptomic analysis of young and old mouse brains. We provide comprehensive datasets of aging-related genes, pathways and ligand-receptor interactions in nearly all brain cell types. Our analysis identified gene signatures that vary in a coordinated manner across cell types and gene sets that are regulated in a cell-type specific manner, even at times in opposite directions. These data reveal that aging, rather than inducing a universal program, drives a distinct transcriptional course in each cell population, and they highlight key molecular processes, including ribosome biogenesis, underlying brain aging. Overall, these large-scale datasets (accessible online at https://portals.broadinstitute.org/single_cell/study/aging-mouse-brain ) provide a resource for the neuroscience community that will facilitate additional discoveries directed towards understanding and modifying the aging process.

The sequence space accessible to evolving proteins can be enhanced by cellular chaperones that assist biophysically defective clients in navigating complex folding landscapes. It is also possible, however, for proteostasis mechanisms that promote strict quality control to greatly constrain accessible protein sequence space. Unfortunately, most efforts to understand how proteostasis mechanisms influence evolution rely on artificial inhibition or genetic knockdown of specific chaperones. The few experiments that perturb quality control pathways also generally modulate the levels of only individual quality control factors. Here, we use chemical genetic strategies to tune proteostasis networks via natural stress response pathways that regulate levels of entire suites of chaperones and quality control mechanisms. Specifically, we upregulate the unfolded protein response (UPR) to test the hypothesis that the host endoplasmic reticulum (ER) proteostasis network shapes the sequence space accessible to human immunodeficiency virus-1 (HIV) envelope (Env) protein. Elucidating factors that enhance or constrain Env sequence space is critical because Env evolves extremely rapidly, yielding HIV strains with antibody and drug escape mutations. We find that UPR-mediated upregulation of ER proteostasis factors, particularly those controlled by the IRE1-XBP1s UPR arm, globally reduces Env mutational tolerance. Conserved, functionally important Env regions exhibit the largest decreases in mutational tolerance upon XBP1s activation. This phenomenon likely reflects strict quality control endowed by XBP1s-mediated remodeling of the ER proteostasis environment. Intriguingly and in contrast, specific regions of Env, including regions targeted by broadly neutralizing antibodies, display enhanced mutational tolerance when XBP1s is activated, hinting at a role for host proteostasis network hijacking in potentiating antibody escape. These observations reveal a key function for proteostasis networks in decreasing instead of expanding the sequence space accessible to client proteins, while also demonstrating that the host ER proteostasis network profoundly shapes the mutational tolerance of Env in ways that could have important consequences for HIV adaptation.

Osteogenesis imperfecta (OI) is typically caused by autosomal dominant mutations in genes encoding collagen type-I, most commonly resulting in Gly→Ser triple-helical domain substitutions that disrupt collagen folding and/or stability. Here, we test the hypothesis that upregulating the endoplasmic reticulum (ER) proteo-stasis network via the unfolded protein response (UPR) can improve the folding and secretion of the clinically severe, prototypical OI-causingCOL1A1p.G425S collagen-α1(I) variant. We first show that small molecules that activate the entire UPR by causing global ER protein misfolding stress severely ablate collagen-I secretion from both G425S Colα1(I)- and wild-type (WT) Colα1(I)-expressing primary fibroblasts. In contrast, stress-independent, specific induction of just the UPR’s XBP1s transcriptional response can enhance collagen-I secretion from G425S Colα1(I) patient primary fibroblasts up to ~300% of basal levels. Notably, the effect is selective – collagen-I secretion from WT Colα1(I)-expressing healthy donor primary fibroblasts is unaltered by XBP1s. XBP1s pathway activation appears to post-translationally enhance the folding/assembly and secretion of G425S Colα1(I), as only modest impacts on collagen-I transcription or synthesis are observed. Consistent with this notion, we find that the stable, triple-helical collagen-I secreted by XBP1s-activated G425S α1(I) patient fibroblasts includes a higher proportion of the mutant α1(I) polypeptide than the collagen-I secreted under basal ER proteostasis conditions. We note that consistent reproducibility of these results is dependent on as yet unascertained experimental variables. Still, these promising observations suggest the potential for ER proteo-stasis network modulation to improve mutant collagen proteostasis in the collagenopathies, motivating further investigation of the effect’s generality, underlying mechanism, and potential therapeutic benefits.

The sequence space accessible to evolving proteins can be enhanced by cellular chaperones that assist biophysically defective clients in navigating complex folding landscapes. It is also possible, at least in theory, for proteostasis mechanisms that promote strict quality control to greatly constrain accessible protein sequence space. Unfortunately, most efforts to understand how proteostasis mechanisms influence evolution rely on artificial inhibition or genetic knockdown of specific chaperones. The few experiments that perturb quality control pathways also generally modulate the levels of only individual quality control factors. Here, we use chemical genetic strategies to tune proteostasis networks via natural stress response pathways that regulate the levels of entire suites of chaperones and quality control mechanisms. Specifically, we upregulate the unfolded protein response (UPR) to test the hypothesis that the host endoplasmic reticulum (ER) proteostasis network shapes the sequence space accessible to human immunodeficiency virus-1 (HIV-1) envelope (Env) protein. Elucidating factors that enhance or constrain Env sequence space is critical because Env evolves extremely rapidly, yielding HIV strains with antibody- and drug-escape mutations. We find that UPR-mediated upregulation of ER proteostasis factors, particularly those controlled by the IRE1-XBP1s UPR arm, globally reduces Env mutational tolerance. Conserved, functionally important Env regions exhibit the largest decreases in mutational tolerance upon XBP1s induction. Our data indicate that this phenomenon likely reflects strict quality control endowed by XBP1s-mediated remodeling of the ER proteostasis environment. Intriguingly, and in contrast, specific regions of Env, including regions targeted by broadly neutralizing antibodies, display enhanced mutational tolerance when XBP1s is induced, hinting at a role for host proteostasis network hijacking in potentiating antibody escape. These observations reveal a key function for proteostasis networks in decreasing instead of expanding the sequence space accessible to client proteins, while also demonstrating that the host ER proteostasis network profoundly shapes the mutational tolerance of Env in ways that could have important consequences for HIV adaptation.

Summary To gain insights into the molecular mechanisms and pathways involved in the activation of γ-herpesvirus (MHV68)-specific T cell receptor transnuclear (TN) CD8+ T cells, we performed a comprehensive transcriptomic analysis. Upon viral infection, we observed differential expression of several thousand transcripts encompassing various networks and pathways in activated TN cells compared with their naive counterparts. Activated cells highly upregulated galectin-3. We therefore explored the role of galectin-3 in influencing anti-MHV68 immunity. Galectin-3 was recruited at the immunological synapse during activation of CD8+ T cells and helped constrain their activation. The localization of galectin-3 to immune synapse was evident during the activation of both naive and memory CD8+ T cells. Galectin-3 knockout mice mounted a stronger MHV68-specific CD8+ T cell response to the majority of viral epitopes and led to better viral control. Targeting intracellular galectin-3 in CD8+ T cells may therefore serve to enhance response to efficiently control infections., Graphical Abstract, Highlights • A comprehensive transcriptomic analysis on γ-HV-expanded TN cells was performed • CD8 T cells upregulated Gal-3 that migrated intracellularly toward immune synapse • Stimulated naive and memory cells displayed galectin-3 recruitment to immune synapse • γ−HV-infected Gal-3 KO mice mounted better CD8+ T cells and controlled virus better, Immune Response; Immunology; Transcriptomics; Virology

We describe the preparation, evaluation, and application of an S100A12 protein-conjugated solid support, hereafter the "A12-resin", that can remove 99% of Zn(II) from complex biological solutions without significantly perturbing the concentrations of other metal ions. The A12-resin can be applied to selectively deplete Zn(II) from diverse tissue culture media and from other biological fluids, including human serum. To further demonstrate the utility of this approach, we investigated metabolic, transcriptomic, and metallomic responses of HEK293 cells cultured in medium depleted of Zn(II) using S100A12. The resulting data provide insight into how cells respond to acute Zn(II) deficiency. We expect that the A12-resin will facilitate interrogation of disrupted Zn(II) homeostasis in biological settings, uncovering novel roles for Zn(II) in biology.

Advances in next-generation sequencing technologies have allowed RNA sequencing to become an increasingly time efficient, cost-effective, and accessible tool for genomic research. We present here an automated and miniaturized workflow for RNA library preparation that minimizes reagent usage and processing time required per sample to generate Illumina compatible libraries for sequencing. The reduced-volume libraries show similar behavior to full-scale libraries with comparable numbers of genes detected and reproducible clustering of samples.

The hematopoietic system continuously generates and replenishes the supply of circulating blood cells from embryonic life throughout the entirety of human lifespan. Studies in mouse development have shown that the repertoire of mature blood cell types produced changes dramatically during development and aging, with hematopoietic stem and progenitor cells (HSPCs) adapting their output to meet age-specific physiologic needs. In humans, it is presumed that age-dependent changes in the production of mature blood lineages underlie the tendency of blood disorders to skew toward certain ages of onset. The observations that mature cell output changes throughout life but mechanisms of terminal hematopoietic differentiation within each lineage remain consistent suggest that age-specific hematopoietic states are programmed at the level of HSPCs. Although the developmental changes occurring in mouse hematopoiesis are well documented, the specific changes in human HSPC ontogeny occurring during prenatal development and postnatal aging from newborn, through childhood, and into adulthood are completely unknown. We hypothesized that temporal changes in human hematopoiesis are mediated by age-specific, occasionally transient, HSPC states and that mechanisms of HSPC lineage commitment change over time in order to meet the changing physiologic demands of the developing and aging human. To test this hypothesis, we comprehensively profiled human HSPC cell states from human fetal hematopoiesis through adulthood using single cell RNA sequencing (scRNAseq). We obtained CD34+ HSPCs from 14 different human donors covering a range of ages from first and second trimester fetal liver, umbilical cord blood, and pediatric and adult bone marrow. We obtained high quality sequencing data on a total of 38,873 individual HSPCs after filtering out apoptotic cells, hepatocytes, stromal cells, and endothelial cells. We then identified differentially expressed genes and performed dimensionality reduction and uniform manifold approximation and projection followed by Louvain clustering to identify 32 distinct cell types encompassing primitive stem and multipotent progenitor cells as well as committed progenitors in the myeloid, erythroid and lymphoid lineages. We found age-specific alterations in hematopoietic differentiation trajectories, particularly in the myeloid lineages. Additionally, we discovered an HSC, emerging in mid-gestation and diminishing at birth, with a characteristic immunophenotype and megakaryocyte (Meg) differentiation bias. Differential gene expression analysis of the Meg-biased fetal HSC identified increased expression of the MYB transcription factor relative to other HSCs, potentially illuminating a mechanistic role for MYB in driving megakaryocyte-erythroid progenitor (MEP) differentiation of fetal HSCs, that is distinct from the role of MYB in conferring erythroid differentiation bias at the MEP stage itself. Finally, we used this atlas of human developmental hematopoiesis to map lineage commitment and progenitor states in leukemia, highlighting the translational applicability of this resource. In summary, we have compiled the first comprehensive atlas of HSPCs across human development and aging. This resource allowed us to identify age-specific differentiation trajectories in human hematopoiesis and enabled identification of a Meg-biased fetal-specific HSC. Our research reveals novel mechanisms of maturation and aging of the human hematopoietic system, uncovers transient HSPC states and differentiation trajectories, and establishes a framework for interrogating the differentiation and maturation states of human leukemias that can likely be applied to other blood diseases. As a resource, we expect that this atlas will broadly impact the study of human hematopoietic development and aging, developmental immunology, and the pathophysiology of age-biased blood diseases. Disclosures No relevant conflicts of interest to declare.

The nematode Caenorhabditis elegans has emerged as a genetically tractable animal host in which to study evolutionarily conserved mechanisms of innate immune signaling. We previously showed that the PMK-1 p38 mitogen-activated protein kinase (MAPK) pathway regulates innate immunity of C. elegans through phosphorylation of the CREB/ATF bZIP transcription factor, ATF-7. Here, we have undertaken a genomic analysis of the transcriptional response of C. elegans to infection by Pseudomonas aeruginosa, combining genome-wide expression analysis by RNA-seq with ATF-7 chromatin immunoprecipitation followed by sequencing (ChIP-Seq). We observe that PMK-1-ATF-7 activity regulates a majority of all genes induced by pathogen infection, and observe ATF-7 occupancy in regulatory regions of pathogen-induced genes in a PMK-1-dependent manner. Moreover, functional analysis of a subset of these ATF-7-regulated pathogen-induced target genes supports a direct role for this transcriptional response in host defense. The genome-wide regulation through PMK-1– ATF-7 signaling reveals global control over the innate immune response to infection through a single transcriptional regulator in a simple animal host.Author SummaryInnate immunity is the first line of defense against invading microbes across metazoans. Caenorhabditis elegans lacks adaptive immunity and is therefore particularly dependent on mounting an innate immune response against pathogens. A major component of this response is the conserved PMK-1/p38 MAPK signaling cascade, the activation of which results in phosphorylation of the bZIP transcription factor ATF-7. Signaling via PMK-1 and ATF-7 causes broad transcriptional changes including the induction of many genes that are predicted to have antimicrobial activity including C-type lectins and lysozymes. In this study, we show that ATF-7 directly regulates the majority of innate immune response genes upon pathogen infection of C. elegans, and demonstrate that many ATF-7 targets function to promote pathogen resistance.

The mammalian brain is complex, with multiple cell types performing a variety of diverse functions, but exactly how each cell type is affected in aging remains largely unknown. Here we performed a single-cell transcriptomic analysis of young and old mouse brains. We provide comprehensive datasets of aging-related genes, pathways and ligand-receptor interactions in nearly all brain cell types. Our analysis identified gene signatures that vary in a coordinated manner across cell types and gene sets that are regulated in a cell-type specific manner, even at times in opposite directions. These data reveal that aging, rather than inducing a universal program, drives a distinct transcriptional course in each cell population, and they highlight key molecular processes, including ribosome biogenesis, underlying brain aging. Overall, these large-scale datasets (accessible online at https://portals.broadinstitute.org/single_cell/study/aging-mouse-brain ) provide a resource for the neuroscience community that will facilitate additional discoveries directed towards understanding and modifying the aging process.

Iron and heme play central roles in red blood cell production. However, the mechanisms by which iron and heme levels coordinate erythropoiesis remain incompletely understood. HRI is a heme-regulated kinase that controls translation by phosphorylating eIF2α. Here, we investigate the global impact of iron, heme and HRI on protein translation in vivo in murine primary erythroblasts using ribosome profiling. By defining the underlying changes in translation during iron and HRI deficiencies, we validate known regulators of this process, including Atf4, and identify novel pathways such as co-regulation of ribosomal protein mRNA translation. Surprisingly, we found that heme and HRI pathways, but not iron-regulated pathways, mediate the major protein translational and transcriptional responses to iron deficiency in erythroblasts in vivo and thereby identify previously unappreciated regulators of erythropoiesis. Our genome-wide study uncovers the major impact of the HRI-mediated integrated stress response for the adaptation to iron deficiency anemia.

The transcription factor Max is a basic-helix-loop-helix leucine zipper (bHLHLZ) protein that forms homodimers or interacts with other bHLHLZ proteins, including Myc and Mxd proteins. Among this dynamic network of interactions, the Myc/Max heterodimer has crucial roles in regulating normal cellular processes, but its transcriptional activity is deregulated in a majority of human cancers. Despite this significance, the arsenal of high-quality chemical probes to interrogate these proteins remains limited. We used small molecule microarrays to identify compounds that bind Max in a mechanistically unbiased manner. We discovered the asymmetric polycyclic lactam, KI-MS2-008, which stabilizes the Max homodimer while reducing Myc protein and Myc-regulated transcript levels. KI-MS2-008 also decreases viable cancer cell growth in a Myc-dependent manner and suppresses tumor growth in vivo. This approach demonstrates the feasibility of modulating Max with small molecules and supports altering Max dimerization as an alternative approach to targeting Myc., National Cancer Institute (Grant R01-CA160860), National Cancer Institute (Grant P30-CA14051), National Cancer Institute (Grant U01-CA176152), National Cancer Institute (Grant CA170378PQ2), National Institutes of Health (Grant CA170378PQ2)

Background Ribosomal RNA (rRNA) comprises at least 90% of total RNA extracted from mammalian tissue or cell line samples. Informative transcriptional profiling using massively parallel sequencing technologies requires either enrichment of mature poly-adenylated transcripts or targeted depletion of the rRNA fraction. The latter method is of particular interest because it is compatible with degraded samples such as those extracted from FFPE and also captures transcripts that are not poly-adenylated such as some non-coding RNAs. Here we provide a cross-site study that evaluates the performance of ribosomal RNA removal kits from Illumina, Takara/Clontech, Kapa Biosystems, Lexogen, New England Biolabs and Qiagen on intact and degraded RNA samples. Results We find that all of the kits are capable of performing significant ribosomal depletion, though there are differences in their ease of use. All kits were able to remove ribosomal RNA to below 20% with intact RNA and identify ~ 14,000 protein coding genes from the Universal Human Reference RNA sample at >1FPKM. Analysis of differentially detected genes between kits suggests that transcript length may be a key factor in library production efficiency. Conclusions These results provide a roadmap for labs on the strengths of each of these methods and how best to utilize them. Keywords: RNAseqr; RNA depletion; Illumina; NGS; ABRF; Transcriptomics, National Cancer Institute (U.S.) (Grant P30-CA14051), National Institute of Environmental Health Sciences (Grant P30-ES002109)

Rationale: Neonatal mice have the capacity to regenerate their hearts in response to injury, but this potential is lost after the first week of life. The transcriptional changes that underpin mammalian cardiac regeneration have not been fully characterized at the molecular level. Objective: The objectives of our study were to determine whether myocytes revert the transcriptional phenotype to a less differentiated state during regeneration and to systematically interrogate the transcriptional data to identify and validate potential regulators of this process. Methods and Results: We derived a core transcriptional signature of injury-induced cardiac myocyte (CM) regeneration in mouse by comparing global transcriptional programs in a dynamic model of in vitro and in vivo CM differentiation, in vitro CM explant model, as well as a neonatal heart resection model. The regenerating mouse heart revealed a transcriptional reversion of CM differentiation processes, including reactivation of latent developmental programs similar to those observed during destabilization of a mature CM phenotype in the explant model. We identified potential upstream regulators of the core network, including interleukin 13, which induced CM cell cycle entry and STAT6/STAT3 signaling in vitro. We demonstrate that STAT3/periostin and STAT6 signaling are critical mediators of interleukin 13 signaling in CMs. These downstream signaling molecules are also modulated in the regenerating mouse heart. Conclusions: Our work reveals new insights into the transcriptional regulation of mammalian cardiac regeneration and provides the founding circuitry for identifying potential regulators for stimulating heart regeneration. Keywords: cardiac myocyte; gene expression; growth factors/cytokines; myogenesis regeneration, National Institutes of Health (U.S.) (Grant F32HL104913), National Institutes of Health (U.S.) (Grant K99HL122514), National Institutes of Health (U.S.) (Grant U01HL098179), Natioanal Science Foundation (U.S.) (Grant CBET-0939511)

Predicting and constraining RNA virus evolution require understanding the molecular factors that define the mutational landscape accessible to these pathogens. RNA viruses typically have high mutation rates, resulting in frequent production of protein variants with compromised biophysical properties. Their evolution is necessarily constrained by the consequent challenge to protein folding and function. We hypothesized that host proteostasis mechanisms may be significant determinants of the fitness of viral protein variants, serving as a critical force shaping viral evolution. Here, we test that hypothesis by propagating influenza in host cells displaying chemically-controlled, divergent proteostasis environments. We find that both the nature of selection on the influenza genome and the accessibility of specific mutational trajectories are significantly impacted by host proteostasis. These findings provide new insights into features of host-pathogen interactions that shape viral evolution, and into the potential design of host proteostasis-targeted antiviral therapeutics that are refractory to resistance., National Institutes of Health (U.S.) (Award 1DP2GM119162), National Institutes of Health (U.S.) (Grant P30-ES002109)

Motivation: Gene expression profiling using RNA-seq is a powerful technique for screening RNA species’ landscapes and their dynamics in an unbiased way. While several advanced methods exist for differential expression analysis of RNA-seq data, proper tools to anal.yze RNA-seq time-course have not been proposed. Results: In this study, we use RNA-seq to measure gene expression during the early human T helper 17 (Th17) cell differentiation and T-cell activation (Th0). To quantify Th17-specific gene expression dynamics, we present a novel statistical methodology, DyNB, for analyzing time-course RNA-seq data. We use non-parametric Gaussian processes to model temporal correlation in gene expression and combine that with negative binomial likelihood for the count data. To account for experiment-specific biases in gene expression dynamics, such as differences in cell differentiation efficiencies, we propose a method to rescale the dynamics between replicated measurements. We develop an MCMC sampling method to make inference of differential expression dynamics between conditions. DyNB identifies several known and novel genes involved in Th17 differentiation. Analysis of differentiation efficiencies revealed consistent patterns in gene expression dynamics between different cultures. We use qRT-PCR to validate differential expression and differentiation efficiencies for selected genes. Comparison of the results with those obtained via traditional timepoint-wise analysis shows that time-course analysis together with time rescaling between cultures identifies differentially expressed genes which would not otherwise be detected. Availability: An implementation of the proposed computational methods will be available at http://research.ics.aalto.fi/csb/software/, Academy of Finland (Centre of Excellence in Moleculary Systems Immunology and Physiology Research (2012-2017) Grant 135320), Seventh Framework Programme (European Commission) (Grant EC-FP7-SYBILLA-201106), EU ERASysBio ERA-NET, Sigrid Juslius Foundation, FICS Graduate School

The unique relapsing nature of Plasmodium vivax infection is a major barrier to malaria eradication. Upon infection, dormant liver-stage forms, hypnozoites, linger for weeks to months and then relapse to cause recurrent blood-stage infection. Very little is known about hypnozoite biology; definitive biomarkers are lacking and in vitro platforms that support phenotypic studies are needed. Here, we recapitulate the entire liver stage of P. vivax in vitro, using a multiwell format that incorporates micropatterned primary human hepatocyte co-cultures (MPCCs). MPCCs feature key aspects of P. vivax biology, including establishment of persistent small forms and growing schizonts, merosome release, and subsequent infection of reticulocytes. We find that the small forms exhibit previously described hallmarks of hypnozoites, and we pilot MPCCs as a tool for testing candidate anti-hypnozoite drugs. Finally, we employ a hybrid capture strategy and RNA sequencing to describe the hypnozoite transcriptome and gain insight into its biology. Plasmodium vivax hypnozoites are difficult to study due to the lack of human liver platforms. Gural et al. recapitulated the entire liver stage of P. vivax in vitro, including formation and reactivation of hypnozoites and release of merosomes. Hybrid capture followed by RNA-seq revealed a first look into the hypnozoite transcriptome. Keywords: malaria; plasmodium vivax; liver stage; hypnozoite; RNA-seq; transcriptome; RNA capture; micropatterning; primary human hepatocytes; radical cure, Bill & Melinda Gates Foundation (Grant OPP1023607), National Cancer Institute (U.S.) (Grant P30-CA14051)

OBJECTIVE—The unfolding of type 1 diabetes involves a number of steps: defective immunological tolerance, priming of anti-islet autoimmunity, and destruction of insulin-producing β-cells. A number of genetic loci contribute to susceptibility to type 1 diabetes, but it is unclear which stages of the disease are influenced by the different loci. Here, we analyzed the frequency of type 1 diabetes–risk alleles among individuals from the Diabetes Prevention Trial–Type 1 (DPT-1) clinical trial, which tested a preventive effect of insulin in at-risk relatives of diabetic individuals, all of which presented with autoimmune manifestations but only one-third of which eventually progressed to diabetes. RESEARCH DESIGN AND METHODS—In this study, 708 individuals randomized into DPT-1 were genotyped for 37 single nucleotide polymorphisms in diabetes susceptibility loci. RESULTS—Susceptibility alleles at loci expected to influence immunoregulation (PTPN22, CTLA4, and IL2RA) did not differ between progressors and nonprogressors but were elevated in both groups relative to general population frequencies, as was the INS promoter variant. In contrast, HLA DQB1*0302 and DQB1*0301 differed significantly in progressors versus nonprogressors (DQB*0302, 42.6 vs. 34.7%, P = 0.0047; DQB*0301, 8.6 vs. 14.3%, P = 0.0026). Multivariate analysis of the factors contributing to progression demonstrated that initial titers of anti-insulin autoantibodies (IAAs) could account for some (P = 0.0016) but not all of this effect on progression (P = 0.00038 for the independent effect of the number of DQB*0302 alleles). The INS-23 genotype was most strongly associated with anti-IAAs (median IAA levels in TT individuals, 60 nU/ml; AT, 121; and AA, 192; P = 0.000037) and only suggestively to the outcome of oral insulin administration. CONCLUSIONS—With the exception of HLA, most susceptibility loci tested condition the risk of autoimmunity rather than the risk of failed immunoregulation that results in islet destruction. Future clinical trials might consider genotyping INS-23 in addition to HLA alleles as disease/treatment response modifier.

Evidence from animal models suggests that natural killer (NK) cells can be important players in the development of type 1 diabetes, although data in humans are still sparse. We studied the frequency and activation state of blood NK cells at different stages of human type 1 diabetes, and whether genetic or phenotypic NK cell peculiarities could be associated with an early onset of diabetes. The onset period is marked by a slight reduction in blood NK cells, but these are unusually activated in some patients (γ-interferon expression). This activation status does not correlate, however, with a particularly young age at onset. In contrast, NK cells in patients with long-standing type 1 diabetes had a markedly lower expression of p30/p46 NK-activating receptor molecules compared with those of control subjects. A slightly decreased expression of NKG2D in all type 1 diabetic patients relative to control subjects was observed, independent of the duration of disease, parallel to prior observations in the NOD mouse. Finally, type 1 diabetic patients had an increased frequency of KIR gene haplotypes that include the activating KIR2DS3 gene, with a genetic interaction between the KIR and HLA complexes. The reduced activation of NK cells in individuals with long-standing type 1 diabetes would seem to be a consequence rather than a cause, but other peculiarities may relate to type 1 diabetes pathogenesis.

The capacity to repair different types of DNA damage varies among individuals, making them more or less susceptible to the detrimental health consequences of damage exposures. Current methods for measuring DNA repair capacity (DRC) are relatively labor intensive, often indirect, and usually limited to a single repair pathway. Here, we describe a fluorescence-based multiplex flow-cytometric host cell reactivation assay (FM-HCR) that measures the ability of human cells to repair plasmid reporters, each bearing a different type of DNA damage or different doses of the same type of DNA damage. FM-HCR simultaneously measures repair capacity in any four of the following pathways: nucleotide excision repair, mismatch repair, base excision repair, nonhomologous end joining, homologous recombination, and methylguanine methyltransferase. We show that FM-HCR can measure interindividual DRC differences in a panel of 24 cell lines derived from genetically diverse, apparently healthy individuals, and we show that FM-HCR may be used to identify inhibitors or enhancers of DRC. We further develop a next-generation sequencing-based HCR assay (HCR-Seq) that detects rare transcriptional mutagenesis events due to lesion bypass by RNA polymerase, providing an added dimension to DRC measurements. FM-HCR and HCR-Seq provide powerful tools for exploring relationships among global DRC, disease susceptibility, and optimal treatment.

Objective It has been suggested that polymorphisms in IL1 are correlated with severe and/or erosive rheumatoid arthritis (RA), but the implicated alleles have differed among studies. The aim of this study was to perform a broad and well-powered search for association between allelic polymorphism in IL1A and IL1B and the susceptibility to or severity of RA. Methods Key coding and regulatory regions in IL1A and IL1B were sequenced in 24 patients with RA, revealing 4 novel single-nucleotide polymorphisms (SNPs) in IL1B. These and a comprehensive set of 24 SNPs tagging most of the underlying genetic diversity were genotyped in 3 independent RA case–control sample sets and 1 longitudinal RA cohort, totaling 3,561 patients and 3,062 control subjects. Results No fully significant associations were observed. Analysis of the discovery case–control sample sets indicated a potential association of IL1B promoter region SNPs with susceptibility to RA (for RA3/A, odds ratio [OR] 1.27, P = 0.0021) or with the incidence of radiographic erosions (for RA4/C, OR 1.56, P = 0.036), but these findings were not replicated in independent case–control samples. No association with rheumatoid factor, anti–cyclic citrullinated peptide, or the Disease Activity Score in 28 joints was found. None of the associations previously observed in other studies were replicated here. Conclusion In spite of a broad and highly powered study, we observed no robust, reproducible association between IL1A/B variants and the susceptibility to or severity of RA in white individuals of European descent. Our results provide evidence that, in the majority of cases, polymorphism in IL1A and IL1B is not a major contributor to genetic susceptibility to RA.

To investigate the effects of Lenercept , a recombinant soluble TNF receptor p55 fused to an immunoglobulin heavy chain IgG1, on the balance of pro- and anti-inflammatory mediators in sepsis.Post hoc analysis of a subgroup of patients enrolled in a multicenter phase III, prospective, double-blind, placebo-controlled, randomized study of Lenercept in severe sepsis.Surgical and medical intensive care units, and postoperative recovery room of a tertiary care teaching hospital.A total of 57 patients were enrolled in the multicenter study in our center.Septic patients were randomly assigned to receive either Lenercept 0.125 mg/kg or placebo. The patients were followed for up to 28 days after randomization.Circulating levels of TNF-alpha, IL-6, TNFsR75 and IL-1Ra were measured before and after treatment. The two groups were comparable with regard to age, gender and diagnosis distribution. The total level of TNF-alpha increased significantly in treated patients, compared to patients receiving placebo. The levels of the other inflammatory mediators did not differ between the two groupsLenercept -treated patients experienced a protracted TNF-alpha half-life, leading to higher total TNF-alpha levels throughout the study. However, the treatment had no effects on anti-inflammatory mediators. Therefore, peripheral inflammatory processes might not have been significantly modified by the treatment. This might account for the lack of efficacy this treatment in septic patients

Polymerase colony (polony) technology amplifies multiple individual DNA molecules within a thin acrylamide gel attached to a microscope slide. Each DNA molecule included in the reaction produces an immobilized colony of double-stranded DNA. We genotype these polonies by performing single base extensions with dye-labeled nucleotides, and we demonstrate the accurate quantitation of two allelic variants. We also show that polony technology can determine the phase, or haplotype, of two single- nucleotide polymorphisms (SNPs) by coamplifying distally located targets on a single chromosomal fragment. We correctly determine the genotype and phase of three different pairs of SNPs. In one case, the distance between the two SNPs is 45 kb, the largest distance achieved to date without separating the chromosomes by cloning or somatic cell fusion. The results indicate that polony genotyping and haplotyping may play an important role in understanding the structure of genetic variation.

Polycomb Repressive Complex 2 (PRC2) catalyzes histone H3 lysine 27 tri-methylation (H3K27me3), an epigenetic modification associated with gene repression. H3K27me3 is enriched at the promoters of a large cohort of developmental genes in embryonic stem cells (ESCs). Loss of H3K27me3 leads to a failure of ESCs to properly differentiate, making it difficult to determine the precise roles of PRC2 during lineage commitment. Moreover, while studies suggest that PRC2 prevents DNA methylation, how these two epigenetic regulators coordinate to regulate lineage programs is poorly understood. Using several PRC2 mutant ESC lines that maintain varying levels of H3K27me3, we found that partial maintenance of H3K27me3 allowed for proper temporal activation of lineage genes during directed differentiation of ESCs to spinal motor neurons (SMNs). In contrast, genes that function to specify other lineages failed to be repressed in these cells, suggesting that PRC2 is also necessary for lineage fidelity. We also found that loss of H3K27me3 leads to a modest gain in DNA methylation at PRC2 target regions in both ESCs and in SMNs. Our study demonstrates a critical role for PRC2 in safeguarding lineage decisions and in protecting genes against inappropriate DNA methylation., National Cancer Institute (U.S.) (Cancer Center Support (Core) Grant P30-CA14051), National Institutes of Health (U.S.) (Training Grant T 32 GM007287), Smith Family Foundation (Contract LTR DATED 11/6/09)

The histone H2A variant H2A.Z is essential for embryonic development and for proper control of developmental gene expression programs in embryonic stem cells (ESCs). Divergent regions of amino acid sequence of H2A.Z likely determine its functional specialization compared to core histone H2A. For example, H2A.Z contains three divergent residues in the essential C-terminal acidic patch that reside on the surface of the histone octamer as an uninterrupted acidic patch domain; however, we know little about how these residues contribute to chromatin structure and function. Here, we show that the divergent amino acids Gly92, Asp97, and Ser98 in the H2A.Z C-terminal acidic patch (H2A.Z[superscript AP3]) are critical for lineage commitment during ESC differentiation. H2A.Z is enriched at most H3K4me3 promoters in ESCs including poised, bivalent promoters that harbor both activating and repressive marks, H3K4me3 and H3K27me3 respectively. We found that while H2A.Z[superscript AP3] interacted with its deposition complex and displayed a highly similar distribution pattern compared to wild-type H2A.Z, its enrichment levels were reduced at target promoters. Further analysis revealed that H2A.Z[superscript AP3] was less tightly associated with chromatin, suggesting that the mutant is more dynamic. Notably, bivalent genes in H2A.Z[superscript AP3] ESCs displayed significant changes in expression compared to active genes. Moreover, bivalent genes in H2A.Z[superscript AP3] ESCs gained H3.3, a variant associated with higher nucleosome turnover, compared to wild-type H2A.Z. We next performed single cell imaging to measure H2A.Z dynamics. We found that H2A.Z[superscript AP3] displayed higher mobility in chromatin compared to wild-type H2A.Z by fluorescent recovery after photobleaching (FRAP). Moreover, ESCs treated with the transcriptional inhibitor flavopiridol resulted in a decrease in the H2A.Z[superscript AP3] mobile fraction and an increase in its occupancy at target genes indicating that the mutant can be properly incorporated into chromatin. Collectively, our work suggests that the divergent residues in the H2A.Z acidic patch comprise a unique domain that couples control of chromatin dynamics to the regulation of developmental gene expression patterns during lineage commitment., Massachusetts Life Sciences Center (David H. Koch Institute for Integrative Cancer Research at MIT Core Grant P30-CA14051), National Science Foundation (U.S.). Emergent Behaviors of Integrated Cellular Systems (Grant CBET-0939511), MIT Faculty Start-up Fund, Massachusetts Institute of Technology. Computational and Systems Biology Initiative (Merck & Co. Postdoctoral Fellowship)