Vascular endothelial cells in vivo are exquisitely regulated by their local environment, which is disrupted or absent when using methods such as FACS sorting of cells isolated from animals or in vitro cell culture. Here, we profile the gene expression patterns of undisturbed endothelial cells in living animals using a novel “AngioTag” zebrafish transgenic line that permits isolation of actively translating mRNAs from endothelial cells in their native environment. This transgenic line uses the endothelial cell-specific kdrl promoter to drive expression of an epitope tagged Rpl10a 60 S ribosomal subunit protein, allowing for Translating Ribosome Affinity Purification (TRAP) of actively translating endothelial cell mRNAs. By performing TRAP-RNAseq on AngioTag animals, we demonstrate strong enrichment of endothelial-specific genes and have uncovered both novel endothelial genes and unique endothelial gene expression signatures for different adult organs. Finally, we generated a versatile “UAS: RiboTag” transgenic line to allow a wider array of different zebrafish cell and tissue types to be examined using TRAP-RNAseq methods. These new tools offer an unparalleled resource to study the molecular identity of cells in their normal in vivo context.

Introduction The majority of studies on oxidative phosphorylation in immune cells have been performed in mouse models, necessitating human translation. To understand the impact of oxidative phosphorylation (OXPHOS) deficiency on human immunity, we studied children with primary mitochondrial disease (MtD). Methods scRNAseq analysis of peripheral blood mononuclear cells was performed on matched children with MtD (N = 4) and controls (N = 4). To define B cell function we performed phage display immunoprecipitation sequencing on a cohort of children with MtD (N = 19) and controls (N = 16). Results Via scRNAseq, we found marked reductions in select populations involved in the humoral immune response, especially antigen presenting cells, B cell and plasma populations, with sparing of T cell populations. MTRNR2L8 , a marker of bioenergetic stress, was significantly elevated in populations that were most depleted. mir4485 , a miRNA contained in the intron of MTRNR2L8 , was co-expressed. Knockdown studies of mir4485 demonstrated its role in promoting survival by modulating apoptosis. To determine the functional consequences of our findings on humoral immunity, we studied the antiviral antibody repertoire in children with MtD and controls using phage display and immunoprecipitation sequencing. Despite similar viral exposomes, MtD displayed antiviral antibodies with less robust fold changes and limited polyclonality. Discussion Overall, we show that children with MtD display perturbations in the B cell repertoire which may impact humoral immunity and the ability to clear viral infections.

Modulation of metabolic flux through pyruvate dehydrogenase complex (PDC) plays an important role in T cell activation and differentiation. PDC sits at the transition between glycolysis and the tricarboxylic acid cycle and is a major producer of acetyl-CoA, marking it as a potential metabolic and epigenetic node To understand the role of pyruvate dehydrogenase complex in T cell differentiation, we generated mice deficient in T cell pyruvate dehydrogenase E1A (Pdha) subunit using a CD4-cre recombinase-based strategy. Herein, we show that genetic ablation of PDC activity in T cells (TPdh-/-) leads to marked perturbations in glycolysis, the tricarboxylic acid cycle, and OXPHOS. TPdh-/- T cells became dependent upon substrate level phosphorylation via glycolysis, secondary to depressed OXPHOS. Due to the block of PDC activity, histone acetylation was also reduced, including H3K27, a critical site for CD8+ TM differentiation. Transcriptional and functional profiling revealed abnormal CD8+ TM differentiation in vitro. Collectively, our data indicate that PDC integrates the metabolome and epigenome in CD8+ memory T cell differentiation. Targeting this metabolic and epigenetic node can have widespread ramifications on cellular function.

Acute myeloid leukemia (AML) is characterized by recurring chromosomal abnormalities that encode oncogenic fusion proteins. Approximately 10% of AML cases are associated with a chromosome 16 inversion [inv(16)(p13q22)] or translocation t(16;16)(p13q22). These chromosomal rearrangements results in the formation of the fusion oncogene CBFB-MYH11 and ultimately the fusion protein CBFβ-SMMHC. Although CBFβ-SMMHC was initially considered a dominant negative repressor of RUNX1, we have recently shown that it works together with RUNX1 to activate gene expression through direct target gene binding. RUNX1 is known to regulate gene expression at both transcriptional and epigenetic levels. Aberrant DNA methylation patterns have been reported across different cancer types including AML. Recent studies have shown that DNA methylation negatively influences RUNX1 DNA binding. We hypothesize that changes to DNA methylation and consequently dysregulated RUNX1 DNA binding contribute to the pathogenesis of inv(16) AML. In this study, we examined the epigenetic and transcriptomic landscapes of CBFB-MYH11 knockin mice and assessed the contribution of the resulting fusion protein CBFβ-SMMHC to RUNX1 DNA target binding. To assess the effects of DNA methylation on RUNX1 DNA binding in vitro, we initially performed fluorescence polarization and electrophoretic mobility shift assays using unmethylated and methylated RUNX1 DNA target sequences. RUNX1, RUNX1/CBFβ and RUNX1/CBFβ-SMMHC showed significantly lower binding affinity for methylated DNA compared to unmethylated DNA probes. To determine genome-wide DNA methylation changes, we performed enzymatic methyl sequencing on sorted Lin -Sca1 -c-Kit + bone marrow cells from control and pre-leukemic inv(16) mice. Inv(16) mice showed a significant genome-wide DNA hypermethylation compared to controls (13,478 hyper- vs 705 hypomethylated regions; q-value < 0.05), with 14% of the differentially methylated regions (DMRs) lying in promoters. Genes associated with these DMRs were enriched for gene ontology (GO) processes related to hematopoiesis, such as myeloid cell homeostasis and erythrocyte differentiation, and DNA replication, such as mitotic DNA replication initiation, among others. In addition, several RUNX1 target genes showed promoter hypermethylation, such as Klf1 or the thrombopoietin receptor Mpl. To understand the impact of aberrant DNA methylation on gene expression, we performed RNA-seq using the same cell population. Gene expression changes were observed in inv(16) mice compared to controls, with 741 genes being upregulated and 278 downregulated from a total of 1019 differentially expressed genes (DEG; absolute fold-change ≥ 2, q-value < 0.05). Enriched GO processes identified from the DEGs were associated with immune responses, such as response to interferon-alpha (and beta) and leukocyte activation involved in inflammatory response. We then intersected the DMR (focusing on promoters) and DEG lists to identify genes whose expression might be directly influenced by DNA methylation changes in inv(16) leukemia. Hypermethylation of promoters was associated with 102 genes, with 67% of them showing downregulation of gene expression in inv(16) mice, whereas hypomethylation was associated with 18 genes with 94% showing upregulation of gene expression. These results show that in inv(16) mice, promoter hypermethylation is associated with downregulation of gene expression, and hypomethylation with increased expression. We also confirmed that several RUNX1 target genes showing promoter hypermethylation also had reduced gene expression in inv(16) mice. The initial analysis suggests that CBFβ-SMMHC expression significantly affects the methylation landscape in mice, with important RUNX1 target genes exhibiting promoter hypermethylation and subsequent downregulation of gene expression. We are currently performing chromatin immunocleavage sequencing using the same cell population to determine RUNX1 binding in vivo. Overall, this study explores a novel regulatory mechanism of RUNX1 function through DNA methylation, and ultimately a new role for RUNX1 in the pathogenesis of inv(16) AML.

Despite advances in identifying rare and common genetic variants conferring risk for ADHD, the lack of a transcriptomic understanding of cortico-striatal brain circuitry has stymied a molecular mechanistic understanding of this disorder. To address this gap, we mapped the transcriptome of the caudate nucleus and anterior cingulate cortex in post-mortem tissue from 60 individuals with and without ADHD. Significant differential expression of genes was found in the anterior cingulate cortex and, to a lesser extent, the caudate. Significant downregulation emerged of neurotransmitter gene pathways, particularly glutamatergic, in keeping with models that implicate these neurotransmitters in ADHD. Consistent with the genetic overlap between mental disorders, correlations were found between the cortico-striatal transcriptomic changes seen in ADHD and those seen in other neurodevelopmental and mood disorders. This transcriptomic evidence points to cortico-striatal neurotransmitter anomalies in the pathogenesis of ADHD, consistent with current models of the disorder.

Using adult zebrafish inner ears as a model for sensorineural regeneration, we ablated the mechanosensory receptors and characterized the single-cell epigenome and transcriptome at consecutive time points during hair cell regeneration. We utilized deep learning on the regeneration-induced open chromatin sequences and identified cell-specific transcription factor (TF) motif patterns. Enhancer activity correlated with gene expression and identified potential gene regulatory networks. A pattern of overlapping Sox- and Six-family TF gene expression and binding motifs was detected, suggesting a combinatorial program of TFs driving regeneration and cell identity. Pseudotime analysis of single-cell transcriptomic data suggested that support cells within the sensory epithelium changed cell identity to a "progenitor" cell population that could differentiate into hair cells. We identified a 2.6 kb DNA enhancer upstream of the sox2 promoter that, when deleted, showed a dominant phenotype that resulted in a hair-cell-regeneration-specific deficit in both the lateral line and adult inner ear.

ABSTRACTVascular endothelial cellsin vivoare exquisitely regulated by their local environment, which is disrupted or absent when using methods such as FACS sorting of cells isolated from animals orin vitrocell culture. Here, we profile the gene expression patterns of undisturbed endothelial cells in living animals using a novel “AngioTag” zebrafish transgenic line that permits isolation of actively translating mRNAs from endothelial cells in their native environment. This transgenic line uses the endothelial cell-specifickdrlpromoter to drive expression of an epitope tagged Rpl10a 60S ribosomal subunit protein, allowing for Translating Ribosome Affinity Purification (TRAP) of actively translating endothelial cell mRNAs. By performing TRAP-RNAseq on AngioTag animals, we demonstrate strong enrichment of endothelial specific genes and uncover novel endothelial genes and unique endothelial gene expression signatures for different adult organs. Finally, we generated a versatile “UAS:RiboTag” transgenic line to allow a wider array of different zebrafish cell and tissue types to be examined using TRAP-RNAseq methods. These new tools offer an unparalleled resource to study the molecular identity of cells in their normalin vivocontext.

Single-cell RNA sequencing (scRNA-seq) of human primary tissues is a rapidly emerging tool for investigating human health and disease at the molecular level. However, optimal processing of solid tissues presents a number of technical and logistical challenges, especially for tissues that are only available at autopsy, which includes pancreatic islets, a tissue that is highly relevant to diabetes. To assess the possible effects of different sample preparation protocols on fresh islet samples, we performed a detailed comparison of scRNA-seq data generated with islets isolated from a human donor but processed according to four treatment strategies, including fixation and cryopreservation. We found significant and reproducible differences in the proportion of cell types identified, and more minor effects on cell-specific patterns of gene expression. Fresh islets from a second donor confirmed gene expression signatures of alpha and beta subclusters. These findings may well apply to other tissues, emphasizing the need for careful consideration when choosing processing methods, comparing results between different studies, and/or interpreting data in the context of multiple cell types from preserved tissue.

The MAPK/ERK signaling pathway is an essential regulator of numerous cell processes that are crucial for normal development as well as cancer progression. While much is known regarding MAPK/ERK signal conveyance from the cell membrane to the nucleus, the transcriptional and epigenetic mechanisms that govern gene expression downstream of MAPK signaling are not fully elucidated. This study employed an integrated epigenome analysis approach to interrogate the effects of MAPK/ERK pathway inhibition on the global transcriptome, the active chromatin landscape, and protein–DNA interactions in 501mel melanoma cells. Treatment of these cells with the small-molecule MEK inhibitor AZD6244 induces hyperpigmentation, widespread gene expression changes including alteration of genes linked to pigmentation, and extensive epigenomic reprogramming of transcriptionally distinct regulatory regions associated with the active chromatin mark H3K27ac. Regulatory regions with differentially acetylated H3K27ac regions following AZD6244 treatment are enriched in transcription factor binding motifs of ETV/ETS and ATF family members as well as the lineage-determining factors MITF and SOX10. H3K27ac-dense enhancer clusters known as super-enhancers show similar transcription factor motif enrichment, and furthermore, these super-enhancers are associated with genes encoding MITF, SOX10, and ETV/ETS proteins. Along with genome-wide resetting of the active enhancer landscape, MEK inhibition also results in widespread SOX10 recruitment throughout the genome, including increased SOX10 binding density at H3K27ac-marked enhancers. Importantly, these MEK inhibitor-responsive enhancers marked by H3K27ac and occupied by SOX10 are located near melanocyte lineage-specific and pigmentation genes and overlap numerous human SNPs associated with pigmentation and melanoma phenotypes, highlighting the variants located within these regions for prioritization in future studies. These results reveal the epigenetic reprogramming underlying the re-activation of melanocyte pigmentation and developmental transcriptional programs in 501mel cells in response to MEK inhibition and suggest extensive involvement of a MEK-SOX10 axis in the regulation of these processes. The dynamic chromatin changes identified here provide a rich genomic resource for further analyses of the molecular mechanisms governing the MAPK pathway in pigmentation- and melanocyte-associated diseases.

Vascular endothelial cells are exquisitely regulated by their local environment, which is absent or disrupted when using methods such as FACS or cell culture to study native signaling pathways. Here, we profile the gene expression patterns of undisturbed endothelial cells in living animals using a novel ‘AngioTag’ zebrafish transgenic line that permits isolation of actively translating mRNAs from endothelial cells in their native environment. This transgenic line uses the endothelial cell-specific promoter to drive expression of an epitope tagged Rpl10a 60S ribosomal subunit protein, allowing for Translating Ribosome Affinity Purification (TRAP) of actively translating endothelial cell mRNAs. By performing TRAP-RNAseq on AngioTag animals, we demonstrate strong enrichment of endothelial specific genes and uncover novel endothelially expressed genes. Additionally, we generated a ‘UAS:RiboTag’ transgenic line to allow for the study of a wider array of zebrafish cell and tissue types using TRAP-RNAseq methods. This new tool offers an unparalleled resource to study cause and effect relationships in the context of gene loss or gain of function .

It is well known that development of prostate cancer (PC) can be attributed to somatic mutations of the genome, acquired within proto-oncogenes or tumor-suppressor genes. What is less well understood is how germline variation contributes to disease aggressiveness in PC patients. To map germline modifiers of aggressive neuroendocrine PC, we generated a genetically diverse F2 intercross population using the transgenic TRAMP mouse model and the wild-derived WSB/EiJ (WSB) strain. The relevance of germline modifiers of aggressive PC identified in these mice was extensively correlated in human PC datasets and functionally validated in cell lines.

Melanocyte stem cells (McSCs) and mouse models of hair graying serve as useful systems to uncover mechanisms involved in stem cell self-renewal and the maintenance of regenerating tissues. Interested in assessing genetic variants that influence McSC maintenance, we found previously that heterozygosity for the melanogenesis associated transcription factor, Mitf, exacerbates McSC differentiation and hair graying in mice that are predisposed for this phenotype. Based on transcriptome and molecular analyses of Mitfmi-vga9/+ mice, we report a novel role for MITF in the regulation of systemic innate immune gene expression. We also demonstrate that the viral mimic poly(I:C) is sufficient to expose genetic susceptibility to hair graying. These observations point to a critical suppressor of innate immunity, the consequences of innate immune dysregulation on pigmentation, both of which may have implications in the autoimmune, depigmenting disease, vitiligo.

In 2018, approximately 165,000 new prostate cancer (PC) cases will be diagnosed, and over 29,000 men will succumb to PC in the U.S. alone. The means of assessing outcome in the clinic are inaccurate, and there is a pressing need to more precisely identify men at risk of aggressive PC. We previously identified HIST1H1A as a susceptibility gene for aggressive PC. HIST1H1A encodes H1.1, a member of the linker histone family that is involved in chromatin organization and compaction. To understand the molecular basis of aggressive PC, we have characterized how germline variation modulates susceptibility to neuroendocrine differentiation, which is a form of aggressive PC. Immunohistochemistry studies revealed that HIST1H1A is over-expressed in normal human prostate tissue compared to prostate adenocarcinoma. Functional characterization of HIST1H1A in prostate LNCaP cells indicated that HIST1HA over-expression increased cell growth, as well as the expression of neuroendocrine and epithelial-to-mesenchymal markers in vitro. Assay for Transposase-Accessible Chromatin (ATAC-seq), which is used to assess chromatin compaction and thus the transcriptional availability of individual genomic regions, demonstrated that H1.1 plays a prominent role in modulating Wnt signaling pathway genes, which are implicated in prostate tumorigenesis. These results demonstrate that HIST1H1A is a modulator of aggressive PC susceptibility.

It is unclear how standing genetic variation affects the prognosis of prostate cancer patients. To provide one controlled answer to this problem, we crossed a dominant, penetrant mouse model of prostate cancer to Diversity Outbred mice, a collection of animals that carries over 40 million SNPs. Integration of disease phenotype and SNP variation data in 493 F1 males identified a metastasis modifier locus on Chromosome 8 (LOD = 8.42); further analysis identified the genes Rwdd4, Cenpu, and Casp3 as functional effectors of this locus. Accordingly, analysis of over 5,300 prostate cancer patient samples revealed correlations between the presence of genetic variants at these loci, their expression levels, cancer aggressiveness, and patient survival. We also observed that ectopic overexpression of RWDD4 and CENPU increased the aggressiveness of two human prostate cancer cell lines. In aggregate, our approach demonstrates how well-characterized genetic variation in mice can be harnessed in conjunction with systems genetics approaches to identify and characterize germline modifiers of human disease processes.

Mutations in the gene encoding transcription factor TFAP2A result in pigmentation anomalies in model organisms and premature hair graying in humans. However, the pleiotropic functions of TFAP2A and its redundantly-acting paralogs have made the precise contribution of TFAP2-type activity to melanocyte differentiation unclear. Defining this contribution may help to explain why TFAP2A expression is reduced in advanced-stage melanoma compared to benign nevi. To identify genes with TFAP2A-dependent expression in melanocytes, we profile zebrafish tissue and mouse melanocytes deficient in Tfap2a, and find that expression of a small subset of genes underlying pigmentation phenotypes is TFAP2A-dependent, including Dct, Mc1r, Mlph, and Pmel. We then conduct TFAP2A ChIP-seq in mouse and human melanocytes and find that a much larger subset of pigmentation genes is associated with active regulatory elements bound by TFAP2A. These elements are also frequently bound by MITF, which is considered the "master regulator" of melanocyte development. For example, the promoter of TRPM1 is bound by both TFAP2A and MITF, and we show that the activity of a minimal TRPM1 promoter is lost upon deletion of the TFAP2A binding sites. However, the expression of Trpm1 is not TFAP2A-dependent, implying that additional TFAP2 paralogs function redundantly to drive melanocyte differentiation, which is consistent with previous results from zebrafish. Paralogs Tfap2a and Tfap2b are both expressed in mouse melanocytes, and we show that mouse embryos with Wnt1-Cre-mediated deletion of Tfap2a and Tfap2b in the neural crest almost completely lack melanocytes but retain neural crest-derived sensory ganglia. These results suggest that TFAP2 paralogs, like MITF, are also necessary for induction of the melanocyte lineage. Finally, we observe a genetic interaction between tfap2a and mitfa in zebrafish, but find that artificially elevating expression of tfap2a does not increase levels of melanin in mitfa hypomorphic or loss-of-function mutants. Collectively, these results show that TFAP2 paralogs, operating alongside lineage-specific transcription factors such as MITF, directly regulate effectors of terminal differentiation in melanocytes. In addition, they suggest that TFAP2A activity, like MITF activity, has the potential to modulate the phenotype of melanoma cells.

SOX10 is required for melanocyte development and maintenance, and has been linked to melanoma initiation and progression. However, the molecular mechanisms by which SOX10 guides the appropriate gene expression programs necessary to promote the melanocyte lineage are not fully understood. Here we employ genetic and epigenomic analysis approaches to uncover novel genomic targets and previously unappreciated molecular roles of SOX10 in melanocytes. Through global analysis of SOX10-binding sites and epigenetic characteristics of chromatin states, we uncover an extensive catalog of SOX10 targets genome-wide. Our findings reveal that SOX10 predominantly engages ‘open’ chromatin regions and binds to distal regulatory elements, including novel and previously known melanocyte enhancers. Integrated chromatin occupancy and transcriptome analysis suggest a role for SOX10 in both transcriptional activation and repression to regulate functionally distinct classes of genes. We demonstrate that distinct epigenetic signatures and cis-regulatory sequence motifs predicted to bind putative co-regulatory transcription factors define SOX10-activated and SOX10-repressed target genes. Collectively, these findings uncover a central role of SOX10 as a global regulator of gene expression in the melanocyte lineage by targeting diverse regulatory pathways.

Metastasis is the primary cause of death in breast cancer. Earlier studies using a mammary tumorigenesis mouse model identified Necdin (Ndn) as a germline modifier of metastasis. Differential expression of Ndn induces a gene-expression signature that predicts prognosis in human breast cancer. Additionally, a non-synonymous germline single nucleotide polymorphism (T50C; V17A) in Ndn distinguishes mouse strains with differing metastatic capacities. To better understand how hereditary factors influence metastasis in breast cancer, we characterized NDN-mediated transcription. Haplotype analysis in a well-characterized breast cancer cohort revealed that NDN germline variation is associated with both NDN expression levels and patient outcome. To examine the role of NDN in mammary tumor metastasis and transcriptional regulation, mouse mammary tumor cell lines stably over-expressing either the wildtype 50T or variant 50C Ndn allele were generated. Cells over-expressing Ndn 50T, but not Ndn 50C, exhibited significant decrease in cell invasiveness and pulmonary metastases compared to control cells. Transcriptome analyses identified a 71-gene expression signature that distinguishes cells over-expressing the two Ndn allelic variants. Furthermore, ChIP assays revealed c-Myc, a target gene of NDN, to be differentially regulated by the allelic variants. These data demonstrate that NDN and the T50C allele regulate gene expression and metastasis efficiency.

Retroviruses integrate into the host genome in patterns specific to each virus. Understanding the causes of these patterns can provide insight into viral integration mechanisms, pathology and genome evolution, and is critical to the development of safe gene therapy vectors. We generated murine leukemia virus integrations in human HepG2 and K562 cells and subjected them to second-generation sequencing, using a DNA barcoding technique that allowed us to quantify independent integration events. We characterized >3,700,000 unique integration events in two ENCODE-characterized cell lines. We find that integrations were most highly enriched in a subset of strong enhancers and active promoters. In both cell types, approximately half the integrations were found in <2% of the genome, demonstrating genomic influences even narrower than previously believed. The integration pattern of murine leukemia virus appears to be largely driven by regions that have high enrichment for multiple marks of active chromatin; the combination of histone marks present was sufficient to explain why some strong enhancers were more prone to integration than others. The approach we used is applicable to analyzing the integration pattern of any exogenous element and could be a valuable preclinical screen to evaluate the safety of gene therapy vectors.