Principles of regulatory information conservation between mouse and human

LINE retrotransposons characterize mammalian tissue-specific and evolutionarily dynamic regulatory regions

BACKGROUND

To investigate the mechanisms driving regulatory evolution across tissues, we experimentally mapped promoters, enhancers, and gene expression in the liver, brain, muscle, and testis from ten diverse mammals.

RESULTS

The regulatory landscape around genes included both tissue-shared and tissue-specific regulatory regions, where tissue-specific promoters and enhancers evolved most rapidly. Genomic regions switching between promoters and enhancers were more common across species, and less common across tissues within a single species. Long Interspersed Nuclear Elements (LINEs) played recurrent evolutionary roles: LINE L1s were associated with tissue-specific regulatory regions, whereas more ancient LINE L2s were associated with tissue-shared regulatory regions and with those switching between promoter and enhancer signatures across species.

CONCLUSIONS

Our analyses of the tissue-specificity and evolutionary stability among promoters and enhancers reveal how specific LINE families have helped shape the dynamic mammalian regulome.

Predicting unrecognized enhancer-mediated genome topology by an ensemble machine learning model

Transcriptional enhancers commonly work over long genomic distances to precisely regulate spatiotemporal gene expression patterns. Dissecting the promoters physically contacted by these distal regulatory elements is essential for understanding developmental processes as well as the role of disease-associated risk variants. Modern proximity-ligation assays, like HiChIP and ChIA-PET, facilitate the accurate identification of long-range contacts between enhancers and promoters. However, these assays are technically challenging, expensive, and time-consuming, making it difficult to investigate enhancer topologies, especially in uncharacterized cell types. To overcome these shortcomings, we therefore designed LoopPredictor, an ensemble machine learning model, to predict genome topology for cell types which lack long-range contact maps. To enrich for functional enhancer-promoter loops over common structural genomic contacts, we trained LoopPredictor with both H3K27ac and YY1 HiChIP data. Moreover, the integration of several related multi-omics features facilitated identifying and annotating the predicted loops. LoopPredictor is able to efficiently identify cell type–specific enhancer-mediated loops, and promoter–promoter interactions, with a modest feature input requirement. Comparable to experimentally generated H3K27ac HiChIP data, we found that LoopPredictor was able to identify functional enhancer loops. Furthermore, to explore the cross-species prediction capability of LoopPredictor, we fed mouse multi-omics features into a model trained on human data and found that the predicted enhancer loops outputs were highly conserved. LoopPredictor enables the dissection of cell type–specific long-range gene regulation and can accelerate the identification of distal disease-associated risk variants.

The 4-NQO mouse model: An update on a well-established in vivo model of oral carcinogenesis

The early detection and management of oral premalignant lesions (OPMDs) improve their outcomes. Animal models that mimic histological and biological processes of human oral carcinogenesis may help to improve the identification of OPMD at-risk of progression into oral squamous cell carcinoma and to develop preventive strategies for the entire field of cancerization. No animal model is perfectly applicable for investigating human oral carcinogenesis. However, the 4-nitroquinoline 1-oxide (4-NQO) mouse model is well established and mimics several morphological, histological, genomic and molecular features of human oral carcinogenesis. Some of the reasons for the success of this model include its reproducible experimental conditions with limited variation, the possibility of realizing longitudinal studies with invasive intervention or gene manipulation, and sample availability for all stages of oral carcinogenesis, especially premalignant lesions. Moreover, the role of histological and molecular alterations in the field of cancerization (i.e., macroscopically healthy mucosa exposed to a carcinogen) during oral carcinogenesis can be easily explored using this model. In this review, we discuss the advantages and drawbacks of this model for studying human oral carcinogenesis. In summary, the 4-NQO-induced murine oral cancer model is relevant for investigating human oral carcinogenesis, including the immune microenvironment, and for evaluating therapeutic and chemoprevention agents.

Functional interrogation of HOXA9 regulome in MLLr leukemia via reporter-based CRISPR/Cas9 screen

Aberrant HOXA9 expression is a hallmark of most aggressive acute leukemias, notably those with KMT2A (MLL) gene rearrangements. HOXA9 overexpression not only predicts poor diagnosis and outcome but also plays a critical role in leukemia transformation and maintenance. However, our current understanding of HOXA9 regulation in leukemia is limited, hindering development of therapeutic strategies. Here, we generated the HOXA9-mCherry knock-in reporter cell lines to dissect HOXA9 regulation. By utilizing the reporter and CRISPR/Cas9 screens, we identified transcription factors controlling HOXA9 expression, including a novel regulator, USF2, whose depletion significantly down-regulated HOXA9 expression and impaired MLLr leukemia cell proliferation. Ectopic expression of Hoxa9 rescued impaired leukemia cell proliferation upon USF2 loss. Cut and Run analysis revealed the direct occupancy of USF2 at HOXA9 promoter in MLLr leukemia cells. Collectively, the HOXA9 reporter facilitated the functional interrogation of the HOXA9 regulome and has advanced our understanding of the molecular regulation network in HOXA9-driven leukemia.

Tumor suppressor p53: from engaging DNA to target gene regulation

The p53 transcription factor confers its potent tumor suppressor functions primarily through the regulation of a large network of target genes. The recent explosion of next generation sequencing protocols has enabled the study of the p53 gene regulatory network (GRN) and underlying mechanisms at an unprecedented depth and scale, helping us to understand precisely how p53 controls gene regulation. Here, we discuss our current understanding of where and how p53 binds to DNA and chromatin, its pioneer-like role, and how this affects gene regulation. We provide an overview of the p53 GRN and the direct and indirect mechanisms through which p53 affects gene regulation. In particular, we focus on delineating the ubiquitous and cell type-specific network of regulatory elements that p53 engages; reviewing our understanding of how, where, and when p53 binds to DNA and the mechanisms through which these events regulate transcription. Finally, we discuss the evolution of the p53 GRN and how recent work has revealed remarkable differences between vertebrates, which are of particular importance to cancer researchers using mouse models.

Mice Are Not Humans: The Case of p53

Mouse models represent powerful tools in cancer research that have influenced much of our understanding of tumor development and our current anticancer strategies. Yet, central signaling pathways differ considerably between mouse and human, including gene regulatory networks (GRNs) of the most important tumor suppressor p53 and its oncogenic sibling p63.

Evaluating the probability of CRISPR-based gene drive contaminating another species

The probability D that a given clustered regularly interspaced short palindromic repeats (CRISPR)-based gene drive element contaminates another, nontarget species can be estimated by the following Drive Risk Assessment Quantitative Estimate (DRAQUE) Equation: D = h y b + t r a n s f × e x p r e s s × c u t × f l a n k × i m m u n e × n o n e x t i n c t with hyb = probability of hybridization between the target species and a nontarget species; transf = probability of horizontal transfer of a piece of DNA containing the gene drive cassette from the target species to a nontarget species (with no hybridization); express = probability that the Cas9 and guide RNA genes are expressed; cut = probability that the CRISPR-guide RNA recognizes and cuts at a DNA site in the new host; flank = probability that the gene drive cassette inserts at the cut site; immune = probability that the immune system does not reject Cas9-expressing cells; nonextinct = probability of invasion of the drive within the population. We discuss and estimate each of the seven parameters of the equation, with particular emphasis on possible transfers within insects, and between rodents and humans. We conclude from current data that the probability of a gene drive cassette to contaminate another species is not insignificant. We propose strategies to reduce this risk and call for more work on estimating all the parameters of the formula.

Supervised enhancer prediction with epigenetic pattern recognition and targeted validation

Enhancers are important noncoding elements, but they have been traditionally hard to characterize experimentally. The development of massively parallel assays allows the characterization of large numbers of enhancers for the first time. Here, we developed a framework using Drosophila STARR-seq to create shape-matching filters based on meta-profiles of epigenetic features. We integrated these features with supervised machine-learning algorithms to predict enhancers. We further demonstrated our model could be transferred to predict enhancers in mammals. We comprehensively validated the predictions using a combination of in vivo and in vitro approaches, involving transgenic assays in mouse and transduction-based reporter assays in human cell lines (153 enhancers in total). The results confirmed our model can accurately predict enhancers in different species without re-parameterization. Finally, we examined the transcription-factor binding patterns at predicted enhancers versus promoters. We demonstrated that these patterns enable the construction of a secondary model effectively discriminating between enhancers and promoters.

Comparative transcriptomics of primary cells in vertebrates

Gene expression profiles in homologous tissues have been observed to be different between species, which may be due to differences between species in the gene expression program in each cell type, but may also reflect differences in cell type composition of each tissue in different species. Here, we compare expression profiles in matching primary cells in human, mouse, rat, dog, and chicken using Cap Analysis Gene Expression (CAGE) and short RNA (sRNA) sequencing data from FANTOM5. While we find that expression profiles of orthologous genes in different species are highly correlated across cell types, in each cell type many genes were differentially expressed between species. Expression of genes with products involved in transcription, RNA processing, and transcriptional regulation was more likely to be conserved, while expression of genes encoding proteins involved in intercellular communication was more likely to have diverged during evolution. Conservation of expression correlated positively with the evolutionary age of genes, suggesting that divergence in expression levels of genes critical for cell function was restricted during evolution. Motif activity analysis showed that both promoters and enhancers are activated by the same transcription factors in different species. An analysis of expression levels of mature miRNAs and of primary miRNAs identified by CAGE revealed that evolutionary old miRNAs are more likely to have conserved expression patterns than young miRNAs. We conclude that key aspects of the regulatory network are conserved, while differential expression of genes involved in cell-to-cell communication may contribute greatly to phenotypic differences between species.

Abeliophyllum distichum Nakai alleviates postmenopausal osteoporosis in ovariectomized rats and prevents RANKL-induced osteoclastogenesis in vitro

ETHNOPHARMACOLOGICAL RELEVANCE

Abeliophyllum distichum Nakai (AD), called Miseon, is one of Korea's monotypic endemic species. As a folk remedy, the AD has been used to treat inflammatory disease, stomachaches, diarrhea, and gynecologic disease in Korea. Some researchers have reported that the AD has anti-cancer, anti-inflammatory, and anti-oxidant effects. But the protective effect of AD leaf for osteoporosis has not been reported yet.

AIM OF THE STUDY

This study aimed to analyze the effects and mechanism of AD-ethyl acetate fraction (EA) extract on the osteoporosis, one of the gynecologic disease.

MATERIALS AND METHODS

The RAW 264.7 cells were used as a model for RANKL-induced osteoclastogenesis. We measured the TRAcP activity, expressions of NFATc1, c-fos, and MAPK to investigate the effect of AD-EA. OVX-induced osteoporosis rat model was used as menopausal osteoporosis. After both ovaries were removed through a surgical procedure, and AD-EA or 17b-estradiol was orally administered for 8 weeks. BMD of femurs was measured as well as the bone morphometric parameter, such as BV/TV, trabecular thickness, number and surface using a micro CT.

RESULTS

AD-EA significantly inhibited TRAcP activity, actin ring formation, pit formation and the expressions of osteoclast-related genes in a dose-dependent manner through the inhibition of the MAPK and c-fos/NFATc1 pathway. In addition, low dose administration of AD-EA improved the deterioration of trabecular bone microarchitecture caused by OVX through the inhibition of bone resorption by TRAcP and CTK.

CONCLUSIONS

These results suggest that AD-EA may contribute to the therapy of osteoporosis caused by menopause in women.

Discovery of regulatory noncoding variants in individual cancer genomes by using cis-X

We developed cis-X, a computational method for discovering regulatory noncoding variants in cancer by integrating whole-genome and transcriptome sequencing data from a single cancer sample. cis-X first finds aberrantly cis-activated genes that exhibit allele-specific expression accompanied by an elevated outlier expression. It then searches for causal noncoding variants that may introduce aberrant transcription factor binding motifs or enhancer hijacking by structural variations. Analysis of 13 T-lineage acute lymphoblastic leukemias identified a recurrent intronic variant predicted to cis-activate the TAL1 oncogene, a finding validated in vivo by chromatin immunoprecipitation sequencing of a patient-derived xenograft. Candidate oncogenes include the prolactin receptor PRLR activated by a focal deletion that removes a CTCF-insulated neighborhood boundary. cis-X may be applied to pediatric and adult solid tumors that are aneuploid and heterogeneous. In contrast to existing approaches, which require large sample cohorts, cis-X enables the discovery of regulatory noncoding variants in individual cancer genomes.

The contribution of non-coding regulatory elements to cardiovascular disease

Cardiovascular disease collectively accounts for a quarter of deaths worldwide. Genome-wide association studies across a range of cardiovascular traits and pathologies have highlighted the prevalence of common non-coding genetic variants within candidate loci. Here, we review genetic, epigenomic and molecular approaches to investigate the contribution of non-coding regulatory elements in cardiovascular biology. We then discuss recent insights on the emerging role of non-coding variation in predisposition to cardiovascular disease, with a focus on novel mechanistic examples from functional genomics studies. Lastly, we consider the clinical significance of these findings at present, and some of the current challenges facing the field.

Perspectives on ENCODE

The Encylopedia of DNA Elements (ENCODE) Project launched in 2003 with the long-term goal of developing a comprehensive map of functional elements in the human genome. These included genes, biochemical regions associated with gene regulation (for example, transcription factor binding sites, open chromatin, and histone marks) and transcript isoforms. The marks serve as sites for candidate cis-regulatory elements (cCREs) that may serve functional roles in regulating gene expression1. The project has been extended to model organisms, particularly the mouse. In the third phase of ENCODE, nearly a million and more than 300,000 cCRE annotations have been generated for human and mouse, respectively, and these have provided a valuable resource for the scientific community.

CHANGE-seq reveals genetic and epigenetic effects on CRISPR-Cas9 genome-wide activity

Current methods can illuminate the genome-wide activity of CRISPR-Cas9 nucleases, but are not easily scalable to the throughput needed to fully understand the principles that govern Cas9 specificity. Here we describe ‘circularization for high-throughput analysis of nuclease genome-wide effects by sequencing’ (CHANGE-seq), a scalable, automatable tagmentation-based method for measuring the genome-wide activity of Cas9 in vitro. We applied CHANGE-seq to 110 sgRNA targets across 13 therapeutically relevant loci in human primary T-cells and identified 201,934 off-target sites, enabling the training of a machine learning model to predict off-target activity. Comparing matched genome-wide off-target, chromatin modification and accessibility, and transcriptional data, we found that cellular off-target activity was two to four times more likely to occur near active promoters, enhancers, and transcribed regions. Finally, CHANGE-seq analysis of 6 targets across 8 individual genomes revealed that human single-nucleotide variation had significant effects on activity at ~15.2% of off-target sites analyzed. CHANGE-seq is a simplified, sensitive, and scalable approach to understanding the specificity of genome editors.

RECAP reveals the true statistical significance of ChIP-seq peak calls

Motivation

Chromatin Immunopreciptation (ChIP)-seq is used extensively to identify sites of transcription factor binding or regions of epigenetic modifications to the genome. A key step in ChIP-seq analysis is peak calling, where genomic regions enriched for ChIP versus control reads are identified. Many programs have been designed to solve this task, but nearly all fall into the statistical trap of using the data twice—once to determine candidate enriched regions, and again to assess enrichment by classical statistical hypothesis testing. This double use of the data invalidates the statistical significance assigned to enriched regions, thus the true significance or reliability of peak calls remains unknown.

Results

Using simulated and real ChIP-seq data, we show that three well-known peak callers, MACS, SICER and diffReps, output biased P-values and false discovery rate estimates that can be many orders of magnitude too optimistic. We propose a wrapper algorithm, RECAP, that uses resampling of ChIP-seq and control data to estimate a monotone transform correcting for biases built into peak calling algorithms. When applied to null hypothesis data, where there is no enrichment between ChIP-seq and control, P-values recalibrated by RECAP are approximately uniformly distributed. On data where there is genuine enrichment, RECAP P-values give a better estimate of the true statistical significance of candidate peaks and better false discovery rate estimates, which correlate better with empirical reproducibility. RECAP is a powerful new tool for assessing the true statistical significance of ChIP-seq peak calls.

Availability and implementation

The RECAP software is available through www.perkinslab.ca or on github at https://github.com/theodorejperkins/RECAP.

Supplementary information

Supplementary data are available at Bioinformatics online.

TETs compete with DNMT3 activity in pluripotent cells at thousands of methylated somatic enhancers

Mammalian cells stably maintain high levels of DNA methylation despite expressing both positive (DNMT3A/B) and negative (TET1-3) regulators. Here, we analyzed the independent and combined effects of these regulators on the DNA methylation landscape using a panel of knockout human embryonic stem cell (ESC) lines. The greatest impact on global methylation levels was observed in DNMT3-deficient cells, including reproducible focal demethylation at thousands of normally methylated loci. Demethylation depends on TET expression and occurs only when both DNMT3s are absent. Dynamic loci are enriched for hydroxymethylcytosine and overlap with subsets of putative somatic enhancers that are methylated in ESCs and can be activated upon differentiation. We observe similar dynamics in mouse ESCs that were less frequent in epiblast stem cells (EpiSCs) and scarce in somatic tissues, suggesting a conserved pluripotency-linked mechanism. Taken together, our data reveal tightly regulated competition between DNMT3s and TETs at thousands of somatic regulatory sequences within pluripotent cells.

GCN5 acetylation is required for craniofacial chondrocyte maturation

Development of the craniofacial structures requires the precise differentiation of cranial neural crest cells into osteoblasts or chondrocytes. Here, we explore the epigenetic and non-epigenetic mechanisms that are required for the development of craniofacial chondrocytes. We previously demonstrated that the acetyltransferase activity of the highly conserved acetyltransferase GCN5, or KAT2A, is required for murine craniofacial development. We show that Gcn5 is required cell autonomously in the cranial neural crest. Moreover, GCN5 is required for chondrocyte development following the arrival of the cranial neural crest within the pharyngeal arches. Using a combination of in vivo and in vitro inhibition of GCN5 acetyltransferase activity, we demonstrate that GCN5 is a potent activator of chondrocyte maturation, acting to control chondrocyte maturation and size increase during pre-hypertrophic maturation to hypertrophic chondrocytes. Rather than acting as an epigenetic regulator of histone H3K9 acetylation, our findings suggest GCN5 primarily acts as a non-histone acetyltransferase to regulate chondrocyte development. Here, we investigate the contribution of GCN5 acetylation to the activity of the mTORC1 pathway. Our findings indicate that GCN5 acetylation is required for activation of this pathway, either via direct activation of mTORC1 or through indirect mechanisms. We also investigate one possibility of how mTORC1 activity is regulated through RAPTOR acetylation, which is hypothesized to enhance mTORC1 downstream phosphorylation. This study contributes to our understanding of the specificity of acetyltransferases, and the cell type specific roles in which these enzymes function.

An integrative view of the regulatory and transcriptional landscapes in mouse hematopoiesis

Thousands of epigenomic data sets have been generated in the past decade, but it is difficult for researchers to effectively use all the data relevant to their projects. Systematic integrative analysis can help meet this need, and the VISION project was established for validated systematic integration of epigenomic data in hematopoiesis. Here, we systematically integrated extensive data recording epigenetic features and transcriptomes from many sources, including individual laboratories and consortia, to produce a comprehensive view of the regulatory landscape of differentiating hematopoietic cell types in mouse. By using IDEAS as our integrative and discriminative epigenome annotation system, we identified and assigned epigenetic states simultaneously along chromosomes and across cell types, precisely and comprehensively. Combining nuclease accessibility and epigenetic states produced a set of more than 200,000 candidate cis-regulatory elements (cCREs) that efficiently capture enhancers and promoters. The transitions in epigenetic states of these cCREs across cell types provided insights into mechanisms of regulation, including decreases in numbers of active cCREs during differentiation of most lineages, transitions from poised to active or inactive states, and shifts in nuclease accessibility of CTCF-bound elements. Regression modeling of epigenetic states at cCREs and gene expression produced a versatile resource to improve selection of cCREs potentially regulating target genes. These resources are available from our VISION website to aid research in genomics and hematopoiesis.

Toward Development of the Male Pill: A Decade of Potential Non-hormonal Contraceptive Targets

With the continued steep rise of the global human population, and the paucity of safe and practical contraceptive options available to men, the need for development of effective and reversible non-hormonal methods of male fertility control is widely recognized. Currently there are several contraceptive options available to men, however, none of the non-hormonal alternatives have been clinically approved. To advance progress in the development of a safe and reversible contraceptive for men, further identification of novel reproductive tract-specific druggable protein targets is required. Here we provide an overview of genes/proteins identified in the last decade as specific or highly expressed in the male reproductive tract, with deletion phenotypes leading to complete male infertility in mice. These phenotypes include arrest of spermatogenesis and/or spermiogenesis, abnormal spermiation, abnormal spermatid morphology, abnormal sperm motility, azoospermia, globozoospermia, asthenozoospermia, and/or teratozoospermia, which are all desirable outcomes for a novel male contraceptive. We also consider other associated deletion phenotypes that could impact the desirability of a potential contraceptive. We further discuss novel contraceptive targets underscoring promising leads with the objective of presenting data for potential druggability and whether collateral effects may exist from paralogs with close sequence similarity.

The Role of Chromatin Accessibility in cis-Regulatory Evolution

Transcription factor (TF) binding is determined by sequence as well as chromatin accessibility. Although the role of accessibility in shaping TF-binding landscapes is well recorded, its role in evolutionary divergence of TF binding, which in turn can alter cis-regulatory activities, is not well understood. In this work, we studied the evolution of genome-wide binding landscapes of five major TFs in the core network of mesoderm specification, between Drosophila melanogaster and Drosophila virilis, and examined its relationship to accessibility and sequence-level changes. We generated chromatin accessibility data from three important stages of embryogenesis in both Drosophila melanogaster and Drosophila virilis and recorded conservation and divergence patterns. We then used multivariable models to correlate accessibility and sequence changes to TF-binding divergence. We found that accessibility changes can in some cases, for example, for the master regulator Twist and for earlier developmental stages, more accurately predict binding change than is possible using TF-binding motif changes between orthologous enhancers. Accessibility changes also explain a significant portion of the codivergence of TF pairs. We noted that accessibility and motif changes offer complementary views of the evolution of TF binding and developed a combined model that captures the evolutionary data much more accurately than either view alone. Finally, we trained machine learning models to predict enhancer activity from TF binding and used these functional models to argue that motif and accessibility-based predictors of TF-binding change can substitute for experimentally measured binding change, for the purpose of predicting evolutionary changes in enhancer activity.

Multi-species annotation of transcriptome and chromatin structure in domesticated animals

BACKGROUND

Comparative genomics studies are central in identifying the coding and non-coding elements associated with complex traits, and the functional annotation of genomes is a critical step to decipher the genotype-to-phenotype relationships in livestock animals. As part of the Functional Annotation of Animal Genomes (FAANG) action, the FR-AgENCODE project aimed to create reference functional maps of domesticated animals by profiling the landscape of transcription (RNA-seq), chromatin accessibility (ATAC-seq) and conformation (Hi-C) in species representing ruminants (cattle, goat), monogastrics (pig) and birds (chicken), using three target samples related to metabolism (liver) and immunity (CD4+ and CD8+ T cells).

RESULTS

RNA-seq assays considerably extended the available catalog of annotated transcripts and identified differentially expressed genes with unknown function, including new syntenic lncRNAs. ATAC-seq highlighted an enrichment for transcription factor binding sites in differentially accessible regions of the chromatin. Comparative analyses revealed a core set of conserved regulatory regions across species. Topologically associating domains (TADs) and epigenetic A/B compartments annotated from Hi-C data were consistent with RNA-seq and ATAC-seq data. Multi-species comparisons showed that conserved TAD boundaries had stronger insulation properties than species-specific ones and that the genomic distribution of orthologous genes in A/B compartments was significantly conserved across species.

CONCLUSIONS

We report the first multi-species and multi-assay genome annotation results obtained by a FAANG project. Beyond the generation of reference annotations and the confirmation of previous findings on model animals, the integrative analysis of data from multiple assays and species sheds a new light on the multi-scale selective pressure shaping genome organization from birds to mammals. Overall, these results emphasize the value of FAANG for research on domesticated animals and reinforces the importance of future meta-analyses of the reference datasets being generated by this community on different species.

Discordant gene responses to radiation in humans and mice and the role of hematopoietically humanized mice in the search for radiation biomarkers

The mouse (Mus musculus) is an extensively used model of human disease and responses to stresses such as ionizing radiation. As part of our work developing gene expression biomarkers of radiation exposure, dose, and injury, we have found many genes are either up-regulated (e.g. CDKN1A, MDM2, BBC3, and CCNG1) or down-regulated (e.g. TCF4 and MYC) in both species after irradiation at ~4 and 8 Gy. However, we have also found genes that are consistently up-regulated in humans and down-regulated in mice (e.g. DDB2, PCNA, GADD45A, SESN1, RRM2B, KCNN4, IFI30, and PTPRO). Here we test a hematopoietically humanized mouse as a potential in vivo model for biodosimetry studies, measuring the response of these 14 genes one day after irradiation at 2 and 4 Gy, and comparing it with that of human blood irradiated ex vivo, and blood from whole body irradiated mice. We found that human blood cells in the hematopoietically humanized mouse in vivo environment recapitulated the gene expression pattern expected from human cells, not the pattern seen from in vivo irradiated normal mice. The results of this study support the use of hematopoietically humanized mice as an in vivo model for screening of radiation response genes relevant to humans.

Which came first, the transcriptional regulator or its target genes? An evolutionary perspective into the construction of eukaryotic regulons

Eukaryotic regulons are regulatory units formed by a set of genes under the control of the same transcription factor (TF). Despite the functional plasticity, TFs are highly conserved and recognize the same DNA sequences in different organisms. One of the main factors that confer regulatory specificity is the distribution of the binding sites of the TFs along the genome, allowing the configuration of different transcriptional regulatory networks (TRNs) from the same regulator. A similar scenario occurs between tissues of the same organism, where a TRN can be rewired by epigenetic factors, modulating the accessibility of the TF to its binding sites. In this article we discuss concepts that can help to formulate testable hypotheses about the construction of regulons, exploring the presence and absence of the elements that form a TRN throughout the evolution of an ancestral lineage. This article is part of a Special Issue entitled: Transcriptional Profiles and Regulatory Gene Networks edited by Dr. Federico Manuel Giorgi and Dr. Shaun Mahony.

Accelerated evolution at chaperone promoters among Antarctic notothenioid fishes

BACKGROUND

Antarctic fishes of the Notothenioidei suborder constitutively upregulate multiple inducible chaperones, a highly derived adaptation that preserves proteostasis in extreme cold, and represent a system for studying the evolution of gene frontloading. We screened for Hsf1-binding sites, as Hsf1 is a master transcription factor of the heat shock response, and highly-conserved non-coding elements within proximal promoters of chaperone genes across 10 Antarctic notothens, 2 subpolar notothens, and 17 perciform fishes. We employed phylogenetic models of molecular evolution to determine whether (i) changes in motifs associated with Hsf1-binding and/or (ii) relaxed purifying selection or exaptation at ancestral cis-regulatory elements coincided with the evolution of chaperone frontloading in Antarctic notothens.

RESULTS

Antarctic notothens exhibited significantly fewer Hsf1-binding sites per bp at chaperone promoters than subpolar notothens and Serranoidei, the most closely-related suborder to Notothenioidei included in this study. 90% of chaperone promoters exhibited accelerated substitution rates among Antarctic notothens relative to other perciformes. The proportion of bases undergoing accelerated evolution (i) was significantly greater in Antarctic notothens than in subpolar notothens and Perciformes in 70% of chaperone genes and (ii) increased among bases that were more conserved among perciformes. Lastly, we detected evidence of relaxed purifying selection and exaptation acting on ancestrally conserved cis-regulatory elements in the Antarctic notothen lineage and its major branches.

CONCLUSION

A large degree of turnover has occurred in Notothenioidei at chaperone promoter regions that are conserved among perciform fishes following adaptation to the cooling of the Southern Ocean. Additionally, derived reductions in Hsf1-binding site frequency suggest cis-regulatory modifications to the classical heat shock response. Of note, turnover events within chaperone promoters were less frequent in the ancestral node of Antarctic notothens relative to younger Antarctic lineages. This suggests that cis-regulatory divergence at chaperone promoters may be greater between Antarctic notothen lineages than between subpolar and Antarctic clades. These findings demonstrate that strong selective forces have acted upon cis-regulatory elements of chaperone genes among Antarctic notothens.

Building gene regulatory networks from scATAC-seq and scRNA-seq using Linked Self Organizing Maps

Rapid advances in single-cell assays have outpaced methods for analysis of those data types. Different single-cell assays show extensive variation in sensitivity and signal to noise levels. In particular, scATAC-seq generates extremely sparse and noisy datasets. Existing methods developed to analyze this data require cells amenable to pseudo-time analysis or require datasets with drastically different cell-types. We describe a novel approach using self-organizing maps (SOM) to link scATAC-seq regions with scRNA-seq genes that overcomes these challenges and can generate draft regulatory networks. Our SOMatic package generates chromatin and gene expression SOMs separately and combines them using a linking function. We applied SOMatic on a mouse pre-B cell differentiation time-course using controlled Ikaros over-expression to recover gene ontology enrichments, identify motifs in genomic regions showing similar single-cell profiles, and generate a gene regulatory network that both recovers known interactions and predicts new Ikaros targets during the differentiation process. The ability of linked SOMs to detect emergent properties from multiple types of highly-dimensional genomic data with very different signal properties opens new avenues for integrative analysis of heterogeneous data.

Evolutionary "Experiments" in Symbiosis: The Study of Model Animals Provides Insights into the Mechanisms Underlying the Diversity of Host-Microbe Interactions

Current work in experimental biology revolves around a handful of animal species. Studying only a few organisms limits science to the answers that those organisms can provide. Nature has given us an overwhelming diversity of animals to study, and recent technological advances have greatly accelerated the ability to generate genetic and genomic tools to develop model organisms for research on host-microbe interactions. With the help of such models the authors therefore hope to construct a more complete picture of the mechanisms that underlie crucial interactions in a given metaorganism (entity consisting of a eukaryotic host with all its associated microbial partners). As reviewed here, new knowledge of the diversity of host-microbe interactions found across the animal kingdom will provide new insights into how animals develop, evolve, and succumb to the disease.

Epigenomes in Cardiovascular Disease

Supplemental Digital Content is available in the text.

If unifying principles could be revealed for how the same genome encodes different eukaryotic cells and for how genetic variability and environmental input are integrated to impact cardiovascular health, grand challenges in basic cell biology and translational medicine may succumb to experimental dissection. A rich body of work in model systems has implicated chromatin-modifying enzymes, DNA methylation, noncoding RNAs, and other transcriptome-shaping factors in adult health and in the development, progression, and mitigation of cardiovascular disease. Meanwhile, deployment of epigenomic tools, powered by next-generation sequencing technologies in cardiovascular models and human populations, has enabled description of epigenomic landscapes underpinning cellular function in the cardiovascular system. This essay aims to unpack the conceptual framework in which epigenomes are studied and to stimulate discussion on how principles of chromatin function may inform investigations of cardiovascular disease and the development of new therapies.

The Applicability of Mouse Models to the Study of Human Disease

The laboratory mouse Mus musculus has long been used as a model organism to test hypotheses and treatments related to understanding the mechanisms of disease in humans; however, for these experiments to be relevant, it is important to know the complex ways in which mice are similar to humans and, crucially, the ways in which they differ. In this chapter, an in-depth analysis of these similarities and differences is provided to allow researchers to use mouse models of human disease and primary cells derived from these animal models under the most appropriate and meaningful conditions.

Although there are considerable differences between mice and humans, particularly regarding genetics, physiology, and immunology, a more thorough understanding of these differences and their effects on the function of the whole organism will provide deeper insights into relevant disease mechanisms and potential drug targets for further clinical investigation. Using specific examples of mouse models of human lung disease, i.e., asthma, chronic obstructive pulmonary disease, and pulmonary fibrosis, this chapter explores the most salient features of mouse models of human disease and provides a full assessment of the advantages and limitations of these models, focusing on the relevance of disease induction and their ability to replicate critical features of human disease pathophysiology and response to treatment. The chapter concludes with a discussion on the future of using mice in medical research with regard to ethical and technological considerations.

Dynamic motif occupancy (DynaMO) analysis identifies transcription factors and their binding sites driving dynamic biological processes

Biological processes are usually associated with genome-wide remodeling of transcription driven by transcription factors (TFs). Identifying key TFs and their spatiotemporal binding patterns are indispensable to understanding how dynamic processes are programmed. However, most methods are designed to predict TF binding sites only. We present a computational method, dynamic motif occupancy analysis (DynaMO), to infer important TFs and their spatiotemporal binding activities in dynamic biological processes using chromatin profiling data from multiple biological conditions such as time-course histone modification ChIP-seq data. In the first step, DynaMO predicts TF binding sites with a random forests approach. Next and uniquely, DynaMO infers dynamic TF binding activities at predicted binding sites using their local chromatin profiles from multiple biological conditions. Another landmark of DynaMO is to identify key TFs in a dynamic process using a clustering and enrichment analysis of dynamic TF binding patterns. Application of DynaMO to the yeast ultradian cycle, mouse circadian clock and human neural differentiation exhibits its accuracy and versatility. We anticipate DynaMO will be generally useful for elucidating transcriptional programs in dynamic processes.

Fine-Resolution Mapping of TF Binding and Chromatin Interactions

Transcription factor (TF) binding to DNA is crucial for transcriptional regulation. There are multiple methods for mapping such binding. These methods balance between input requirements, spatial resolution, and compatibility with high-throughput automation. Here, we describe SLIM-ChIP (short-fragment-enriched, low-input, indexed MNase ChIP), which combines enzymatic fragmentation of chromatin and on-bead indexing to address these desiderata. SLIM-ChIP reproduces a high-resolution binding map of yeast Reb1 comparable with existing methods, yet with less input material and full compatibility with high-throughput procedures. We demonstrate the robustness and flexibility of SLIM-ChIP by probing additional factors in yeast and mouse. Finally, we show that SLIM-ChIP provides information on the chromatin landscape surrounding the bound transcription factor. We identify a class of Reb1 sites where the proximal −1 nucleosome tightly interacts with Reb1 and maintains unidirectional transcription. SLIM-ChIP is an attractive solution for mapping DNA binding proteins and charting the surrounding chromatin occupancy landscape at a single-cell level.

•

SLIM-ChIP is a low-input, robust, high-resolution, automatable TF mapping protocol

•

SLIM-ChIP is applicable to a range of TFs from yeast to mammals

•

DNA fragments from SLIM-ChIP provide targeted footprinting at the bound site

•

RSC-mediated Reb1-nucleosome interactions affect promoter directionality

DNA methylation profiling allows for characterization of atrial and ventricular cardiac tissues and hiPSC-CMs

Background

Cardiac disease modelling using human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CM) requires thorough insight into cardiac cell type differentiation processes. However, current methods to discriminate different cardiac cell types are mostly time-consuming, are costly and often provide imprecise phenotypic evaluation. DNA methylation plays a critical role during early heart development and cardiac cellular specification. We therefore investigated the DNA methylation pattern in different cardiac tissues to identify CpG loci for further cardiac cell type characterization.

Results

An array-based genome-wide DNA methylation analysis using Illumina Infinium HumanMethylation450 BeadChips led to the identification of 168 differentially methylated CpG loci in atrial and ventricular human heart tissue samples (n = 49) from different patients with congenital heart defects (CHD). Systematic evaluation of atrial-ventricular DNA methylation pattern in cardiac tissues in an independent sample cohort of non-failing donor hearts and cardiac patients using bisulfite pyrosequencing helped us to define a subset of 16 differentially methylated CpG loci enabling precise characterization of human atrial and ventricular cardiac tissue samples. This defined set of reproducible cardiac tissue-specific DNA methylation sites allowed us to consistently detect the cellular identity of hiPSC-CM subtypes.

Conclusion

Testing DNA methylation of only a small set of defined CpG sites thus makes it possible to distinguish atrial and ventricular cardiac tissues and cardiac atrial and ventricular subtypes of hiPSC-CMs. This method represents a rapid and reliable system for phenotypic characterization of in vitro-generated cardiomyocytes and opens new opportunities for cardiovascular research and patient-specific therapy.

Electronic supplementary material

The online version of this article (10.1186/s13148-019-0679-0) contains supplementary material, which is available to authorized users.

Distribution of Initiation Times Reveals Mechanisms of Transcriptional Regulation in Single Cells

Transcription is the dominant point of control of gene expression. Biochemical studies have revealed key molecular components of transcription and their interactions, but the dynamics of transcription initiation in cells is still poorly understood. This state of affairs is being remedied with experiments that observe transcriptional dynamics in single cells using fluorescent reporters. Quantitative information about transcription initiation dynamics can also be extracted from experiments that use electron micrographs of RNA polymerases caught in the act of transcribing a gene (Miller spreads). Inspired by these data, we analyze a general stochastic model of transcription initiation and elongation and compute the distribution of transcription initiation times. We show that different mechanisms of initiation leave distinct signatures in the distribution of initiation times that can be compared to experiments. We analyze published data from micrographs of RNA polymerases transcribing ribosomal RNA genes in Escherichia coli and compare the observed distributions of interpolymerase distances with the predictions from previously hypothesized mechanisms for the regulation of these genes. Our analysis demonstrates the potential of measuring the distribution of time intervals between initiation events as a probe for dissecting mechanisms of transcription initiation in live cells.

Shaping the human brain: evolutionary cis-regulatory plasticity drives changes in synaptic activity-controlled adaptive gene expression

Neuronal activity-induced gene expression programs involved in synaptic structure- and plasticity-related functions are similar in mice and humans, yet bear distinct features. These include gains or losses of activity-responsiveness of certain genes and differences in gene induction profiles. Here, we discuss a possible origin of dissimilarities in activity-regulated transcription between species. We highlight that while synapse-to-nucleus signalling pathways are evolutionarily conserved, cis-regulatory plasticity has been driving species-specific remodelling of the activity-controlled enhancer landscape, thereby affecting gene regulation. In particular, evolutionary rearrangements of transcription factor binding site placements together with potential species-dependent developmental stage- and/or cell type-specific epigenetic and other trans-acting mechanisms are most likely at least in part accountable for between-species diversity in activity-regulated transcription. It is conceivable that cis-regulatory plasticity may have equipped the synaptic activity-driven adaptive gene program in human neurons with unique, species-specific qualities.

Identification and Conservation Analysis of Cis-Regulatory Elements in Pig Liver

The liver plays a key role in metabolism and affects pig production. However, the functional annotation of noncoding regions of the pig liver remains poorly understood. We revealed the landscape of cis-regulatory elements and their functional characterization in pig liver. We identified 102,373 cis-regulatory elements in the pig liver, including enhancers, promoters, super-enhancers, and broad H3K4me3 domains, and highlighted 26 core transcription regulatory factors in the pig liver as well. We found similarity of cis-regulatory elements among those of pigs, humans, and cattle. Despite the low proportion of functionally conserved enhancers (~30%) between pig and human liver tissue, ~78% of the pig liver enhancer orthologues sequence could play an enhancer role in other human tissues. Additionally, we observed that the ratio of consistent super-enhancer-associated genes was significantly higher than the ratio of functionally conserved super-enhancers. Approximately 54% of the core regulation factors driven by super-enhancers were consistent across the liver from these three species. Our pig liver annotation and functional characterization studies provide a system and resource for noncoding annotation for future gene regulatory studies in pigs. Furthermore, our study also showed the high level functional conservation of cis-regulatory elements in mammals; it also improved our understanding of regulation function of mammal cis-regulatory elements.

Conservation and divergence of the p53 gene regulatory network between mice and humans

Understanding the p53 tumor suppressor pathway remains crucial for the design of anticancer strategies. Studies in human tumors and mouse models help to unravel the molecular mechanisms that underlie the p53 signaling pathway. Yet, the p53 gene regulatory network (GRN) is not the same in mice and humans. The comparison of the regulatory networks of p53 in mice and humans reveals that gene up- and down-regulation by p53 are distinctly affected during evolution. Importantly, gene up-regulation by p53 underwent more rapid evolution and gene down-regulation has been evolutionarily constrained. This difference stems from the two major mechanisms employed by p53 to regulate gene expression: up-regulation through direct p53 target gene binding and indirect down-regulation through the p53-p21-DREAM pathway. More than 1000 genes have been identified to differ in their p53-dependent expression between mice and humans. Analysis of p53 gene expression profiles and p53 binding data reveal that turnover of p53 binding sites is the major mechanism underlying extensive variation in p53-dependent gene up-regulation. Only a core set of high-confidence genes appears to be directly regulated by p53 in both species. In contrast to up-regulation, p53-induced down-regulation is well conserved between mice and humans and controls cell cycle genes. Here a curated data set is provided that extends the previously established web-atlas at www.targetgenereg.org to assess the p53 response of any human gene of interest and its mouse ortholog. Taken together, the analysis reveals a limited translation potential from mouse models to humans for the p53 GRN.

Comparison of reprogramming factor targets reveals both species-specific and conserved mechanisms in early iPSC reprogramming

BACKGROUND

Both human and mouse fibroblasts can be reprogrammed to pluripotency with Oct4, Sox2, Klf4, and c-Myc (OSKM) transcription factors. While both systems generate pluripotency, human reprogramming takes considerably longer than mouse.

RESULTS

To assess additional similarities and differences, we sought to compare the binding of the reprogramming factors between the two systems. In human fibroblasts, the OSK factors initially target many more closed chromatin sites compared to mouse. Despite this difference, the intra- and intergenic distribution of target sites, target genes, primary binding motifs, and combinatorial binding patterns between the reprogramming factors are largely shared. However, while many OSKM binding events in early mouse cell reprogramming occur in syntenic regions, only a limited number is conserved in human.

CONCLUSIONS

Our findings suggest similar general effects of OSKM binding across these two species, even though the detailed regulatory networks have diverged significantly.

Genome-wide profiling of histone H3K27 acetylation featured fatty acid signalling in pancreatic beta cells in diet-induced obesity in mice

AIMS/HYPOTHESIS

Epigenetic regulation of gene expression has been implicated in the pathogenesis of obesity and type 2 diabetes. However, detailed information, such as key transcription factors in pancreatic beta cells that mediate environmental effects, is not yet available.

METHODS

To analyse genome-wide cis-regulatory profiles and transcriptome of pancreatic islets derived from a diet-induced obesity (DIO) mouse model, we conducted chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-Seq) of histone H3 lysine 27 acetylation (histone H3K27ac) and high-throughput RNA sequencing. Transcription factor-binding motifs enriched in differential H3K27ac regions were examined by de novo motif analysis. For the predicted transcription factors, loss of function experiments were performed by transfecting specific siRNA in INS-1, a rat beta cell line, with and without palmitate treatment. Epigenomic and transcriptional changes of possible target genes were evaluated by ChIP and quantitative RT-PCR.

RESULTS

After long-term feeding with a high-fat diet, C57BL/6J mice were obese and mildly glucose intolerant. Among 39,350 islet cis-regulatory regions, 13,369 and 4610 elements showed increase and decrease in ChIP-Seq signals, respectively, significantly associated with global change in gene expression. Remarkably, increased H3K27ac showed a distinctive genomic localisation, mainly in the proximal-promoter regions, revealing enriched elements for nuclear respiratory factor 1 (NRF1), GA repeat binding protein α (GABPA) and myocyte enhancer factor 2A (MEF2A) by de novo motif analysis, whereas decreased H3K27ac was enriched for v-maf musculoaponeurotic fibrosarcoma oncogene family protein K (MAFK), a known negative regulator of beta cells. By siRNA-mediated knockdown of NRF1, GABPA or MEF2A we found that INS-1 cells exhibited downregulation of fatty acid β-oxidation genes in parallel with decrease in the associated H3K27ac. Furthermore, in line with the epigenome in DIO mice, palmitate treatment caused increase in H3K27ac and induction of β-oxidation genes; these responses were blunted when NRF1, GABPA or MEF2A were suppressed.

CONCLUSIONS/INTERPRETATION

These results suggest novel roles for DNA-binding proteins and fatty acid signalling in obesity-induced epigenomic regulation of beta cell function.

DATA AVAILABILITY

The next-generation sequencing data in the present study were deposited at ArrayExpress. RNA-Seq: Dataset name: ERR2538129 (Control), ERR2538130 (Diet-induced obesity) Repository name and number: E-MTAB-6718 - RNA-Seq of pancreatic islets derived from mice fed a long-term high-fat diet against chow-fed controls. ChIP-Seq: Dataset name: ERR2538131 (Control), ERR2538132 (Diet-induced obesity) Repository name and number: E-MTAB-6719 - H3K27ac ChIP-Seq of pancreatic islets derived from mice fed a long-term high-fat diet (HFD) against chow-fed controls.

Regulatory Architecture of the Neuronal Cacng2/Tarpγ2 Gene Promoter: Multiple Repressive Domains, a Polymorphic Regulatory Short Tandem Repeat, and Bidirectional Organization with Co-regulated lncRNAs

CACNG2 (TARPγ2, Stargazin) is a multi-functional regulator of excitatory neurotransmission and has been implicated in the pathological processes of several brain diseases. Cacng2 function is dependent upon expression level, but currently, little is known about the molecular mechanisms that control expression of this gene. To address this deficit and investigate disease-related gene variants, we have cloned and characterized the rat Cacng2 promoter and have defined three major features: (i) multiple repressive domains that include an array of RE-1 silencing transcription factor (REST) elements, and a calcium regulatory element-binding factor (CaRF) element, (ii) a (poly-GA) short tandem repeat (STR), and (iii) bidirectional organization with expressed lncRNAs. Functional activity of the promoter was demonstrated in transfected neuronal cell lines (HT22 and PC12), but although selective removal of REST and CaRF domains was shown to enhance promoter-driven transcription, the enhanced Cacng2 promoter constructs were still about fivefold weaker than a comparable rat Synapsin-1 promoter sequence. Direct evidence of REST activity at the Cacng2 promoter was obtained through co-transfection with an established dominant-negative REST (DNR) construct. Investigation of the GA-repeat STR revealed polymorphism across both animal strains and species, and size variation was also observed in absence epilepsy disease model cohorts (Genetic Absence Epilepsy Rats, Strasbourg [GAERS] and non-epileptic control [NEC] rats). These data provide evidence of a genotype (STR)-phenotype correlation that may be unique with respect to proximal gene regulatory sequence in the demonstrated absence of other promoter, or 3' UTR variants in GAERS rats. However, although transcriptional regulatory activity of the STR was demonstrated in further transfection studies, we did not find a GAERS vs. NEC difference, indicating that this specific STR length variation may only be relevant in the context of other (Cacna1h and Kcnk9) gene variants in this disease model. Additional studies revealed further (bidirectional) complexity at the Cacng2 promoter, and we identified novel, co-regulated, antisense rat lncRNAs that are paired with Cacng2 mRNA. These studies have provided novel insights into the organization of a synaptic protein gene promoter, describing multiple repressive and modulatory domains that can mediate diverse regulatory inputs.

Development and Evolutionary Constraints in Animals

We review the evolutionary importance of developmental mechanisms in constraining evolutionary changes in animals—in other words, developmental constraints. We focus on hard constraints that can act on macroevolutionary timescales. In particular, we discuss the causes and evolutionary consequences of the ancient metazoan constraint that differentiated cells cannot divide and constraints against changes of phylotypic stages in vertebrates and other higher taxa. We conclude that in all cases these constraints are caused by complex and highly controlled global interactivity of development, the disturbance of which has grave consequences. Mutations that affect such global interactivity almost unavoidably have many deleterious pleiotropic effects, which will be strongly selected against and will lead to long-term evolutionary stasis. The discussed developmental constraints have pervasive consequences for evolution and critically restrict regeneration capacity and body plan evolution.

Correlated Evolution of Two Copulatory Organs via a Single cis-Regulatory Nucleotide Change

Diverse traits often covary between species [1-3]. The possibility that a single mutation could contribute to the evolution of several characters between species [3] is rarely investigated as relatively few cases are dissected at the nucleotide level. Drosophila santomea has evolved additional sex comb sensory teeth on its legs and has lost two sensory bristles on its genitalia. We present evidence that a single nucleotide substitution in an enhancer of the scute gene contributes to both changes. The mutation alters a binding site for the Hox protein Abdominal-B in the developing genitalia, leading to bristle loss, and for another factor in the developing leg, leading to bristle gain. Our study suggests that morphological evolution between species can occur through a single nucleotide change affecting several sexually dimorphic traits. VIDEO ABSTRACT.

Contemporary Circulating Enterovirus D68 Strains Have Acquired the Capacity for Viral Entry and Replication in Human Neuronal Cells

Since the EV-D68 outbreak during the summer of 2014, evidence of a causal link to a type of limb paralysis (AFM) has been mounting. In this article, we describe a neuronal cell culture model (SH-SY5Y cells) in which a subset of contemporary 2014 outbreak strains of EV-D68 show infectivity in neuronal cells, or neurotropism. We confirmed the difference in neurotropism in vitro using primary human neuron cell cultures and in vivo with a mouse paralysis model. Using the SH-SY5Y cell model, we determined that a barrier to viral entry is at least partly responsible for neurotropism. SH-SY5Y cells may be useful in determining if specific EV-D68 genetic determinants are associated with neuropathogenesis, and replication in this cell line could be used as rapid screening tool for identification of neurotropic EV-D68 strains. This may assist with better understanding of pathogenesis and epidemiology and with the development of potential therapies.

Enterovirus D68 (EV-D68) has historically been associated with respiratory illnesses. However, in the summers of 2014 and 2016, EV-D68 outbreaks coincided with a spike in polio-like acute flaccid myelitis/paralysis (AFM/AFP) cases. This raised concerns that EV-D68 could be the causative agent of AFM during these recent outbreaks. To assess the potential neurotropism of EV-D68, we utilized the neuroblastoma-derived neuronal cell line SH-SY5Y as a cell culture model to determine if differential infection is observed for different EV-D68 strains. In contrast to HeLa and A549 cells, which support viral infection of all EV-D68 strains tested, SH-SY5Y cells only supported infection by a subset of contemporary EV-D68 strains, including isolates from the 2014 outbreak. Viral replication and infectivity in SH-SY5Y were assessed using multiple assays: virus production, cytopathic effects, cellular ATP release, and VP1 capsid protein production. Similar differential neurotropism was also observed in differentiated SH-SY5Y cells, primary human neuron cultures, and a mouse paralysis model. Using the SH-SY5Y cell culture model, we determined that barriers to viral binding and entry were at least partly responsible for the differential infectivity phenotype. Transfection of genomic RNA into SH-SY5Y generated virions for all EV-D68 isolates, but only a single round of replication was observed from strains that could not directly infect SH-SY5Y. In addition to supporting virus replication and other functional studies, this cell culture model may help identify the signatures of virulence to confirm epidemiological associations between EV-D68 strains and AFM and allow for the rapid identification and characterization of emerging neurotropic strains.

Nature and nurture of tissue-specific macrophage phenotypes

Macrophages are key players in immunity and tissue homeostasis but can also contribute to a diverse range of human diseases, including cardiovascular diseases. Enhancers, cis-acting DNA elements regulating gene activity, have been shown to be crucial for control of macrophage development and function. The selection and activities of macrophage-specific enhancers are regulated by the combined actions of lineage determining transcription factors (LDTFs) and signal dependent transcription factors (SDTFs) that are specified by developmental origin and tissue-specific signals. As a consequence, each tissue resident macrophage population adopts a distinct phenotype. In this review, we discuss recent work on how environmental factors affect the activation status of enhancers and can lead to long-lasting epigenetic changes resulting in innate immune memory. Furthermore, we discuss how non-coding genetic variation affects gene expression by altering transcription factor binding through local and domain-wide mechanisms. These findings have implications for interpretation of non-coding risk alleles that are associated with human disease and efforts to target macrophages for therapeutic purposes.

Prediction of gene regulatory enhancers across species reveals evolutionarily conserved sequence properties

Genomic regions with gene regulatory enhancer activity turnover rapidly across mammals. In contrast, gene expression patterns and transcription factor binding preferences are largely conserved between mammalian species. Based on this conservation, we hypothesized that enhancers active in different mammals would exhibit conserved sequence patterns in spite of their different genomic locations. To investigate this hypothesis, we evaluated the extent to which sequence patterns that are predictive of enhancers in one species are predictive of enhancers in other mammalian species by training and testing two types of machine learning models. We trained support vector machine (SVM) and convolutional neural network (CNN) classifiers to distinguish enhancers defined by histone marks from the genomic background based on DNA sequence patterns in human, macaque, mouse, dog, cow, and opossum. The classifiers accurately identified many adult liver, developing limb, and developing brain enhancers, and the CNNs outperformed the SVMs. Furthermore, classifiers trained in one species and tested in another performed nearly as well as classifiers trained and tested on the same species. We observed similar cross-species conservation when applying the models to human and mouse enhancers validated in transgenic assays. This indicates that many short sequence patterns predictive of enhancers are largely conserved. The sequence patterns most predictive of enhancers in each species matched the binding motifs for a common set of TFs enriched for expression in relevant tissues, supporting the biological relevance of the learned features. Thus, despite the rapid change of active enhancer locations between mammals, cross-species enhancer prediction is often possible. Our results suggest that short sequence patterns encoding enhancer activity have been maintained across more than 180 million years of mammalian evolution.

Human intraepithelial lymphocytes

The location of intraepithelial lymphocytes (IEL) between epithelial cells, their effector memory, cytolytic and inflammatory phenotype positions them to kill infected epithelial cells and protect the intestine against pathogens. Human TCRαβ+CD8αβ+ IEL have the dual capacity to recognize modified self via natural killer (NK) receptors (autoreactivity) as well as foreign antigen via the T cell receptor (TCR), which is accomplished in mouse by two cell subsets, the naturally occurring TCRαβ+CD8αα+ and adaptively induced TCRαβ+CD8αβ+ IEL subsets, respectively. The private/oligoclonal nature of the TCR repertoire of both human and mouse IEL suggests local environmental factors dictate the specificity of IEL responses. The line between sensing of foreign antigens and autoreactivity is blurred for IEL in celiac disease, where recognition of stress ligands by induced activating NK receptors in conjunction with inflammatory signals such as IL-15 can result in low-affinity TCR/non-cognate antigen and NK receptor/stress ligand interactions triggering destruction of intestinal epithelial cells.

From lipid locus to drug target through human genomics

In the last decade, over 175 genetic loci have robustly been associated to levels of major circulating blood lipids. Most loci are specific to one or two lipids, whereas some (SUGP1, ZPR1, TRIB1, HERPUD1, and FADS1) are associated to all. While exposing the polygenic architecture of circulating lipids and the underpinnings of dyslipidaemia, these genome-wide association studies (GWAS) have provided further evidence of the critical role that lipids play in coronary heart disease (CHD) risk, as indicated by the 2.7-fold enrichment for macrophage gene expression in atherosclerotic plaques and the association of 25 loci (such as PCSK9, APOB, ABCG5-G8, KCNK5, LPL, HMGCR, NPC1L1, CETP, TRIB1, ABO, PMAIP1-MC4R, and LDLR) with CHD. These GWAS also confirmed known and commonly used therapeutic targets, including HMGCR (statins), PCSK9 (antibodies), and NPC1L1 (ezetimibe). As we head into the post-GWAS era, we offer suggestions for how to move forward beyond genetic risk loci, towards refining the biology behind the associations and identifying causal genes and therapeutic targets. Deep phenotyping through lipidomics and metabolomics will refine and increase the resolution to find causal and druggable targets, and studies aimed at demonstrating gene transcriptional and regulatory effects of lipid associated loci will further aid in identifying these targets. Thus, we argue the need for deeply phenotyped, large genetic association studies to reduce costs and failures and increase the efficiency of the drug discovery pipeline. We conjecture that in the next decade a paradigm shift will tip the balance towards a data-driven approach to therapeutic target development and the application of precision medicine where human genomics takes centre stage.