Up-regulation of cholesterol synthesis pathways and limited neurodegeneration in a knock-in Sod1 mutant mouse model of ALS.

Amyotrophic lateral sclerosis (ALS) is a severe neurodegenerative disorder affecting brain and spinal cord motor neurons. Mutations in the copper/zinc superoxide dismutase gene ( SOD1 ) are associated with ∼20% of inherited and 1-2% of sporadic ALS cases. Much has been learned from mice expressing transgenic copies of mutant SOD1, which typically involve high-level transgene expression, thereby differing from ALS patients expressing one mutant gene copy. To generate a model that more closely represents patient gene expression, we created a knock-in point mutation (G85R, a human ALS-causing mutation) in the endogenous mouse Sod1 gene, leading to mutant SOD1 G85R protein expression. Heterozygous Sod1 G85R mutant mice resemble wild type, whereas homozygous mutants have reduced body weight and lifespan, a mild neurodegenerative phenotype, and express very low mutant SOD1 protein levels with no detectable SOD1 activity. Homozygous mutants exhibit partial neuromuscular junction denervation at 3-4 months of age. Spinal cord motor neuron transcriptome analyses of homozygous Sod1 G85R mice revealed up-regulation of cholesterol synthesis pathway genes compared to wild type. Transcriptome and phenotypic features of these mice are similar to Sod1 knock-out mice, suggesting the Sod1 G85R phenotype is largely driven by loss of SOD1 function. By contrast, cholesterol synthesis genes are down-regulated in severely affected human TgSOD1 G93A transgenic mice at 4 months. Our analyses implicate dysregulation of cholesterol or related lipid pathway genes in ALS pathogenesis. The Sod1 G85R knock-in mouse is a useful ALS model to examine the importance of SOD1 activity in control of cholesterol homeostasis and motor neuron survival. SIGNIFICANCE STATEMENT: Amyotrophic lateral sclerosis is a devastating disease involving the progressive loss of motor neurons and motor function for which there is currently no cure. Understanding biological mechanisms leading to motor neuron death is critical for developing new treatments. Using a new knock-in mutant mouse model carrying a Sod1 mutation that causes ALS in patients, and in the mouse, causes a limited neurodegenerative phenotype similar to Sod1 loss-of-function, we show that cholesterol synthesis pathway genes are up-regulated in mutant motor neurons, whereas the same genes are down-regulated in transgenic SOD1 mice with a severe phenotype. Our data implicate dysregulation of cholesterol or other related lipid genes in ALS pathogenesis and provide new insights that could contribute to strategies for disease intervention.

Functional Annotation Routines Used by ABRF Bioinformatics Core Facilities - Observations, Comparisons, and Considerations.

The functional annotation of gene lists is a common analysis routine required for most genomics experiments, and bioinformatics core facilities must support these analyses. In contrast to methods such as the quantitation of RNA-Seq reads or differential expression analysis, our research group noted a lack of consensus in our preferred approaches to functional annotation. To investigate this observation, we selected 4 experiments that represent a range of experimental designs encountered by our cores and analyzed those data with 6 tools used by members of the Association of Biomolecular Resource Facilities (ABRF) Genomic Bioinformatics Research Group (GBIRG). To facilitate comparisons between tools, we focused on a single biological result for each experiment. These results were represented by a gene set, and we analyzed these gene sets with each tool considered in our study to map the result to the annotation categories presented by each tool. In most cases, each tool produces data that would facilitate identification of the selected biological result for each experiment. For the exceptions, Fisher's exact test parameters could be adjusted to detect the result. Because Fisher's exact test is used by many functional annotation tools, we investigated input parameters and demonstrate that, while background set size is unlikely to have a significant impact on the results, the numbers of differentially expressed genes in an annotation category and the total number of differentially expressed genes under consideration are both critical parameters that may need to be modified during analyses. In addition, we note that differences in the annotation categories tested by each tool, as well as the composition of those categories, can have a significant impact on results.

CRISPR-induced exon skipping of ß-catenin reveals tumorigenic mutants driving distinct subtypes of liver cancer.

CRISPR/Cas9-driven cancer modeling studies are based on the disruption of tumor suppressor genes by small insertions or deletions (indels) that lead to frame-shift mutations. In addition, CRISPR/Cas9 is widely used to define the significance of cancer oncogenes and genetic dependencies in loss-of-function studies. However, how CRISPR/Cas9 influences gain-of-function oncogenic mutations is elusive. Here, we demonstrate that single guide RNA targeting exon 3 of Ctnnb1 (encoding ß-catenin) results in exon skipping and generates gain-of-function isoforms in vivo. CRISPR/Cas9-mediated exon skipping of Ctnnb1 induces liver tumor formation in synergy with YAPS127A in mice. We define two distinct exon skipping-induced tumor subtypes with different histological and transcriptional features. Notably, ectopic expression of two exon-skipped ß-catenin transcript isoforms together with YAPS127A phenocopies the two distinct subtypes of liver cancer. Moreover, we identify similar CTNNB1 exon-skipping events in patients with hepatocellular carcinoma. Collectively, our findings advance our understanding of ß-catenin-related tumorigenesis and reveal that CRISPR/Cas9 can be repurposed, in vivo, to study gain-of-function mutations of oncogenes in cancer. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

AAV-delivered suppressor tRNA overcomes a nonsense mutation in mice.

Gene therapy is a potentially curative medicine for many currently untreatable diseases, and recombinant adeno-associated virus (rAAV) is the most successful gene delivery vehicle for in vivo applications1-3. However, rAAV-based gene therapy suffers from several limitations, such as constrained DNA cargo size and toxicities caused by non-physiological expression of a transgene4-6. Here we show that rAAV delivery of a suppressor tRNA (rAAV.sup-tRNA) safely and efficiently rescued a genetic disease in a mouse model carrying a nonsense mutation, and effects lasted for more than 6 months after a single treatment. Mechanistically, this was achieved through a synergistic effect of premature stop codon readthrough and inhibition of nonsense-mediated mRNA decay. rAAV.sup-tRNA had a limited effect on global readthrough at normal stop codons and did not perturb endogenous tRNA homeostasis, as determined by ribosome profiling and tRNA sequencing, respectively. By optimizing the AAV capsid and the route of administration, therapeutic efficacy in various target tissues was achieved, including liver, heart, skeletal muscle and brain. This study demonstrates the feasibility of developing a toolbox of AAV-delivered nonsense suppressor tRNAs operating on premature termination codons (AAV-NoSTOP) to rescue pathogenic nonsense mutations and restore gene function under endogenous regulation. As nonsense mutations account for 11% of pathogenic mutations, AAV-NoSTOP can benefit a large number of patients. AAV-NoSTOP obviates the need to deliver a full-length protein-coding gene that may exceed the rAAV packaging limit, elicit adverse immune responses or cause transgene-related toxicities. It therefore represents a valuable addition to gene therapeutics.

Type I IFN-Driven Immune Cell Dysregulation in Rat Autoimmune Diabetes.

Type 1 diabetes is a chronic autoimmune disease, characterized by the immune-mediated destruction of insulin-producing ß cells of pancreatic islets. Essential components of the innate immune antiviral response, including type I IFN and IFN receptor (IFNAR)-mediated signaling pathways, likely contribute to human type 1 diabetes susceptibility. We previously showed that LEW.1WR1 Ifnar1 -/- rats have a significant reduction in diabetes frequency following Kilham rat virus (KRV) infection. To delineate the impact of IFNAR loss on immune cell populations in KRV-induced diabetes, we performed flow cytometric analysis in spleens from LEW.1WR1 wild-type (WT) and Ifnar1 -/- rats after viral infection but before the onset of insulitis and diabetes. We found a relative decrease in CD8+ T cells and NK cells in KRV-infected LEW.1WR1 Ifnar1 -/- rats compared with KRV-infected WT rats; splenic regulatory T cells were diminished in WT but not Ifnar1 -/- rats. In contrast, splenic neutrophils were increased in KRV-infected Ifnar1 -/- rats compared with KRV-infected WT rats. Transcriptional analysis of splenic cells from KRV-infected rats confirmed a reduction in IFN-stimulated genes in Ifnar1 -/- compared with WT rats and revealed an increase in transcripts related to neutrophil chemotaxis and MHC class II. Single-cell RNA sequencing confirmed that MHC class II transcripts are increased in monocytes and macrophages and that numerous types of splenic cells harbor KRV. Collectively, these findings identify dynamic shifts in innate and adaptive immune cells following IFNAR disruption in a rat model of autoimmune diabetes, providing insights toward the role of type I IFNs in autoimmunity.

Dietary suppression of MHC class II expression in intestinal epithelial cells enhances intestinal tumorigenesis.

Little is known about how interactions of diet, intestinal stem cells (ISCs), and immune cells affect early-stage intestinal tumorigenesis. We show that a high-fat diet (HFD) reduces the expression of the major histocompatibility complex class II (MHC class II) genes in intestinal epithelial cells, including ISCs. This decline in epithelial MHC class II expression in a HFD correlates with reduced intestinal microbiome diversity. Microbial community transfer experiments suggest that epithelial MHC class II expression is regulated by intestinal flora. Mechanistically, pattern recognition receptor (PRR) and interferon-gamma (IFNγ) signaling regulates epithelial MHC class II expression. MHC class II-negative (MHC-II-) ISCs exhibit greater tumor-initiating capacity than their MHC class II-positive (MHC-II+) counterparts upon loss of the tumor suppressor Apc coupled with a HFD, suggesting a role for epithelial MHC class II-mediated immune surveillance in suppressing tumorigenesis. ISC-specific genetic ablation of MHC class II increases tumor burden cell autonomously. Thus, HFD perturbs a microbiome-stem cell-immune cell interaction that contributes to tumor initiation in the intestine.

Reduced Neurog3 Gene Dosage Shifts Enteroendocrine Progenitor Towards Goblet Cell Lineage in the Mouse Intestine.

BACKGROUND & AIMS: Transient expression of Neurog3 commits intestinal secretory progenitors to become enteroendocrine-biased progenitors and hence drive enteroendocrine differentiation. Loss of Neurog3 in mouse resulted in the depletion of intestinal enteroendocrine cells (EECs) and an increase in goblet cells. Earlier studies in developing mouse pancreas identified a role of Neurog3 gene dosage in endocrine and exocrine cell fate allocation. We aimed to determine whether Neurog3 gene dosage controls fate choice of enteroendocrine progenitors. METHODS: We acquired mutant Neurog3 reporter mice carrying 2, 1, or null Neurog3 alleles to study Neurog3 gene dosage effect by lineage tracing. Cell types arising from Neurog3+ progenitors were determined by immunohistochemistry using antibodies against intestinal lineage-specific markers. RNA sequencing of sorted Neurog3+/+, Neurog3+/-, or bulk intestinal cells were performed and differentially expressed genes were analyzed. RESULTS: We identified 2731 genes enriched in sorted Neurog3+/+-derived cells in the Neurog3+/+EYFP mouse intestine when compared with bulk duodenum epithelial cells. In the intestine of Neurog3+/-EGFP heterozygous mouse, we observed a 63% decrease in EEC numbers. Many Neurog3-derived cells stained for goblet marker Mucin 2. RNA sequencing of sorted Neurog3+/- cells uncovered enriched expression of genes characteristic for both goblet and enteroendocrine cells, indicating the mixed lineages arose from Neurog3+ progenitors. Consistent with this hypothesis, deletion of both Neurog3 alleles resulted in the total absence of EECs. All Neurog3+-derived cells stained for Mucin 2. CONCLUSIONS: We identified that the fate of Neurog3+ enteroendocrine progenitors is dependent on Neurog3 gene dosage. High Neurog3 gene dosage enforces the commitment of secretory progenitors to an EE lineage, while constraining their goblet cell lineage potential. Transcriptome profiling data was deposited to Gene Ontology omnibus, accession number: GSE149203.

A Translational Model for Venous Thromboembolism: MicroRNA Expression in Hibernating Black Bears.

BACKGROUND: Hibernating American black bears have significantly different clotting parameters than their summer active counterparts, affording them protection against venous thromboembolism during prolonged periods of immobility. We sought to evaluate if significant differences exist between the expression of microRNAs in the plasma of hibernating black bears compared with their summer active counterparts, potentially contributing to differences in hemostasis during hibernation. MATERIALS AND METHODS: MicroRNA sequencing was assessed in plasma from 21 American black bears in summer active (n = 11) and hibernating states (n = 10), and microRNA signatures during hibernating and active state were established using both bear and human genome. MicroRNA targets were predicted using messenger RNA (mRNA) transcripts from black bear kidney cells. In vitro studies were performed to confirm the relationship between identified microRNAs and mRNA expression, using artificial microRNA and human liver cells. RESULTS: Using the bear genome, we identified 15 microRNAs differentially expressed in the plasma of hibernating black bears. Of these microRNAs, three were significantly downregulated (miR-141-3p, miR-200a-3p, and miR-200c-3p), were predicted to target SERPINC1, the gene for antithrombin, and demonstrated regulatory control of the gene mRNA expression in cell studies. CONCLUSIONS: Our findings suggest that the hibernating black bears' ability to maintain hemostasis and achieve protection from venous thromboembolism during prolonged periods of immobility may be due to changes in microRNA signatures and possible upregulation of antithrombin expression.

An improved zebrafish transcriptome annotation for sensitive and comprehensive detection of cell type-specific genes.

The zebrafish is ideal for studying embryogenesis and is increasingly applied to model human disease. In these contexts, RNA-sequencing (RNA-seq) provides mechanistic insights by identifying transcriptome changes between experimental conditions. Application of RNA-seq relies on accurate transcript annotation for a genome of interest. Here, we find discrepancies in analysis from RNA-seq datasets quantified using Ensembl and RefSeq zebrafish annotations. These issues were due, in part, to variably annotated 3' untranslated regions and thousands of gene models missing from each annotation. Since these discrepancies could compromise downstream analyses and biological reproducibility, we built a more comprehensive zebrafish transcriptome annotation that addresses these deficiencies. Our annotation improves detection of cell type-specific genes in both bulk and single cell RNA-seq datasets, where it also improves resolution of cell clustering. Thus, we demonstrate that our new transcriptome annotation can outperform existing annotations, providing an important resource for zebrafish researchers.

An atlas of cell types in the mouse epididymis and vas deferens.

Following testicular spermatogenesis, mammalian sperm continue to mature in a long epithelial tube known as the epididymis, which plays key roles in remodeling sperm protein, lipid, and RNA composition. To understand the roles for the epididymis in reproductive biology, we generated a single-cell atlas of the murine epididymis and vas deferens. We recovered key epithelial cell types including principal cells, clear cells, and basal cells, along with associated support cells that include fibroblasts, smooth muscle, macrophages and other immune cells. Moreover, our data illuminate extensive regional specialization of principal cell populations across the length of the epididymis. In addition to region-specific specialization of principal cells, we find evidence for functionally specialized subpopulations of stromal cells, and, most notably, two distinct populations of clear cells. Our dataset extends on existing knowledge of epididymal biology, and provides a wealth of information on potential regulatory and signaling factors that bear future investigation.

Bioinformatics Core Survey Highlights the Challenges Facing Data Analysis Facilities.

Over the last decade, the cost of -omics data creation has decreased 10-fold, whereas the need for analytical support for those data has increased exponentially. Consequently, bioinformaticians face a second wave of challenges: novel applications of existing approaches (e.g., single-cell RNA sequencing), integration of -omics data sets of differing size and scale (e.g., spatial transcriptomics), as well as novel computational and statistical methods, all of which require more sophisticated pipelines and data management. Nonetheless, bioinformatics cores are often asked to operate under primarily a cost-recovery model, with limited institutional support. Seeing the need to assess bioinformatics core operations, the Association of Biomolecular Resource Facilities Genomics Bioinformatics Research Group conducted a survey to answer questions about staffing, services, financial models, and challenges to better understand the challenges bioinformatics core facilities are currently faced with and will need to address going forward. Of the respondent groups, we chose to focus on the survey data from smaller cores, which made up the majority. Although all cores indicated similar challenges in terms of changing technologies and analysis needs, small cores tended to have the added challenge of funding their operations largely through cost-recovery models with heavy administrative burdens.

DolphinNext: a distributed data processing platform for high throughput genomics.

BACKGROUND: The emergence of high throughput technologies that produce vast amounts of genomic data, such as next-generation sequencing (NGS) is transforming biological research. The dramatic increase in the volume of data, the variety and continuous change of data processing tools, algorithms and databases make analysis the main bottleneck for scientific discovery. The processing of high throughput datasets typically involves many different computational programs, each of which performs a specific step in a pipeline. Given the wide range of applications and organizational infrastructures, there is a great need for highly parallel, flexible, portable, and reproducible data processing frameworks. Several platforms currently exist for the design and execution of complex pipelines. Unfortunately, current platforms lack the necessary combination of parallelism, portability, flexibility and/or reproducibility that are required by the current research environment. To address these shortcomings, workflow frameworks that provide a platform to develop and share portable pipelines have recently arisen. We complement these new platforms by providing a graphical user interface to create, maintain, and execute complex pipelines. Such a platform will simplify robust and reproducible workflow creation for non-technical users as well as provide a robust platform to maintain pipelines for large organizations. RESULTS: To simplify development, maintenance, and execution of complex pipelines we created DolphinNext. DolphinNext facilitates building and deployment of complex pipelines using a modular approach implemented in a graphical interface that relies on the powerful Nextflow workflow framework by providing 1. A drag and drop user interface that visualizes pipelines and allows users to create pipelines without familiarity in underlying programming languages. 2. Modules to execute and monitor pipelines in distributed computing environments such as high-performance clusters and/or cloud 3. Reproducible pipelines with version tracking and stand-alone versions that can be run independently. 4. Modular process design with process revisioning support to increase reusability and pipeline development efficiency. 5. Pipeline sharing with GitHub and automated testing 6. Extensive reports with R-markdown and shiny support for interactive data visualization and analysis. CONCLUSION: DolphinNext is a flexible, intuitive, web-based data processing and analysis platform that enables creating, deploying, sharing, and executing complex Nextflow pipelines with extensive revisioning and interactive reporting to enhance reproducible results.

HIV-1-induced cytokines deplete homeostatic innate lymphoid cells and expand TCF7-dependent memory NK cells.

Human immunodeficiency virus 1 (HIV-1) infection is associated with heightened inflammation and excess risk of cardiovascular disease, cancer and other complications. These pathologies persist despite antiretroviral therapy. In two independent cohorts, we found that innate lymphoid cells (ILCs) were depleted in the blood and gut of people with HIV-1, even with effective antiretroviral therapy. ILC depletion was associated with neutrophil infiltration of the gut lamina propria, type 1 interferon activation, increased microbial translocation and natural killer (NK) cell skewing towards an inflammatory state, with chromatin structure and phenotype typical of WNT transcription factor TCF7-dependent memory T cells. Cytokines that are elevated during acute HIV-1 infection reproduced the ILC and NK cell abnormalities ex vivo. These results show that inflammatory cytokines associated with HIV-1 infection irreversibly disrupt ILCs. This results in loss of gut epithelial integrity, microbial translocation and memory NK cells with heightened inflammatory potential, and explains the chronic inflammation in people with HIV-1.

Genomic Characterization of Endothelial Enhancers Reveals a Multifunctional Role for NR2F2 in Regulation of Arteriovenous Gene Expression.

RATIONALE: Significant progress has revealed transcriptional inputs that underlie regulation of artery and vein endothelial cell fates. However, little is known concerning genome-wide regulation of this process. Therefore, such studies are warranted to address this gap. OBJECTIVE: To identify and characterize artery- and vein-specific endothelial enhancers in the human genome, thereby gaining insights into mechanisms by which blood vessel identity is regulated. METHODS AND RESULTS: Using chromatin immunoprecipitation and deep sequencing for markers of active chromatin in human arterial and venous endothelial cells, we identified several thousand artery- and vein-specific regulatory elements. Computational analysis revealed that NR2F2 (nuclear receptor subfamily 2, group F, member 2) sites were overrepresented in vein-specific enhancers, suggesting a direct role in promoting vein identity. Subsequent integration of chromatin immunoprecipitation and deep sequencing data sets with RNA sequencing revealed that NR2F2 regulated 3 distinct aspects related to arteriovenous identity. First, consistent with previous genetic observations, NR2F2 directly activated enhancer elements flanking cell cycle genes to drive their expression. Second, NR2F2 was essential to directly activate vein-specific enhancers and their associated genes. Our genomic approach further revealed that NR2F2 acts with ERG (ETS-related gene) at many of these sites to drive vein-specific gene expression. Finally, NR2F2 directly repressed only a small number of artery enhancers in venous cells to prevent their activation, including a distal element upstream of the artery-specific transcription factor, HEY2 (hes related family bHLH transcription factor with YRPW motif 2). In arterial endothelial cells, this enhancer was normally bound by ERG, which was also required for arterial HEY2 expression. By contrast, in venous endothelial cells, NR2F2 was bound to this site, together with ERG, and prevented its activation. CONCLUSIONS: By leveraging a genome-wide approach, we revealed mechanistic insights into how NR2F2 functions in multiple roles to maintain venous identity. Importantly, characterization of its role at a crucial artery enhancer upstream of HEY2 established a novel mechanism by which artery-specific expression can be achieved.

Evolutionarily conserved pachytene piRNA loci are highly divergent among modern humans.

In the fetal mouse testis, PIWI-interacting RNAs (piRNAs) guide PIWI proteins to silence transposons but, after birth, most post-pubertal pachytene piRNAs map to the genome uniquely and are thought to regulate genes required for male fertility. In the human male, the developmental classes, precise genomic origins and transcriptional regulation of postnatal piRNAs remain undefined. Here, we demarcate the genes and transcripts that produce postnatal piRNAs in human juvenile and adult testes. As in the mouse, human A-MYB drives transcription of both pachytene piRNA precursor transcripts and messenger RNAs encoding piRNA biogenesis factors. Although human piRNA genes are syntenic to those in other placental mammals, their sequences are poorly conserved. In fact, pachytene piRNA loci are rapidly diverging even among modern humans. Our findings suggest that, during mammalian evolution, pachytene piRNA genes are under few selective constraints. We speculate that pachytene piRNA diversity may provide a hitherto unrecognized driver of reproductive isolation.

Recovery of viable endocrine-specific cells and transcriptomes from human pancreatic islet-engrafted mice.

Human pancreatic islets engrafted into immunodeficient mice serve as an important model for in vivo human diabetes studies. Following engraftment, islet function can be monitored in vivo by measuring circulating glucose and human insulin; however, it will be important to recover viable cells for more complex graft analyses. Moreover, RNA analyses of dissected grafts have not distinguished which hormone-specific cell types contribute to gene expression. We developed a method for recovering live cells suitable for fluorescence-activated cell sorting from human islets engrafted in mice. Although yields of recovered islet cells were relatively low, the ratios of bulk-sorted ß, α, and δ cells and their respective hormone-specific RNA-Seq transcriptomes are comparable pretransplant and posttransplant, suggesting that the cellular characteristics of islet grafts posttransplant closely mirror the original donor islets. Single-cell RNA-Seq transcriptome analysis confirms the presence of appropriate ß, α, and δ cell subsets. In addition, ex vivo perifusion of recovered human islet grafts demonstrated glucose-stimulated insulin secretion. Viable cells suitable for patch-clamp analysis were recovered from transplanted human embryonic stem cell-derived ß cells. Together, our functional and hormone-specific transcriptome analyses document the broad applicability of this system for longitudinal examination of human islet cells undergoing developmental/metabolic/pharmacogenetic manipulation in vivo and may facilitate the discovery of treatments for diabetes.

Diverse repertoire of human adipocyte subtypes develops from transcriptionally distinct mesenchymal progenitor cells.

Single-cell sequencing technologies have revealed an unexpectedly broad repertoire of cells required to mediate complex functions in multicellular organisms. Despite the multiple roles of adipose tissue in maintaining systemic metabolic homeostasis, adipocytes are thought to be largely homogenous with only 2 major subtypes recognized in humans so far. Here we report the existence and characteristics of 4 distinct human adipocyte subtypes, and of their respective mesenchymal progenitors. The phenotypes of these distinct adipocyte subtypes are differentially associated with key adipose tissue functions, including thermogenesis, lipid storage, and adipokine secretion. The transcriptomic signature of "brite/beige" thermogenic adipocytes reveals mechanisms for iron accumulation and protection from oxidative stress, necessary for mitochondrial biogenesis and respiration upon activation. Importantly, this signature is enriched in human supraclavicular adipose tissue, confirming that these cells comprise thermogenic depots in vivo, and explain previous findings of a rate-limiting role of iron in adipose tissue browning. The mesenchymal progenitors that give rise to beige/brite adipocytes express a unique set of cytokines and transcriptional regulators involved in immune cell modulation of adipose tissue browning. Unexpectedly, we also find adipocyte subtypes specialized for high-level expression of the adipokines adiponectin or leptin, associated with distinct transcription factors previously implicated in adipocyte differentiation. The finding of a broad adipocyte repertoire derived from a distinct set of mesenchymal progenitors, and of the transcriptional regulators that can control their development, provides a framework for understanding human adipose tissue function and role in metabolic disease.

CRISPR-SONIC: targeted somatic oncogene knock-in enables rapid in vivo cancer modeling.

CRISPR/Cas9 has revolutionized cancer mouse models. Although loss-of-function genetics by CRISPR/Cas9 is well-established, generating gain-of-function alleles in somatic cancer models is still challenging because of the low efficiency of gene knock-in. Here we developed CRISPR-based Somatic Oncogene kNock-In for Cancer Modeling (CRISPR-SONIC), a method for rapid in vivo cancer modeling using homology-independent repair to integrate oncogenes at a targeted genomic locus. Using a dual guide RNA strategy, we integrated a plasmid donor in the 3'-UTR of mouse ß-actin, allowing co-expression of reporter genes or oncogenes from the ß-actin promoter. We showed that knock-in of oncogenic Ras and loss of p53 efficiently induced intrahepatic cholangiocarcinoma in mice. Further, our strategy can generate bioluminescent liver cancer to facilitate tumor imaging. This method simplifies in vivo gain-of-function genetics by facilitating targeted integration of oncogenes.

HLA Class II Antigen Processing and Presentation Pathway Components Demonstrated by Transcriptome and Protein Analyses of Islet ß-Cells From Donors With Type 1 Diabetes.

Type 1 diabetes studies consistently generate data showing islet ß-cell dysfunction and T cell-mediated anti-ß-cell-specific autoimmunity. To explore the pathogenesis, we interrogated the ß-cell transcriptomes from donors with and without type 1 diabetes using both bulk-sorted and single ß-cells. Consistent with immunohistological studies, ß-cells from donors with type 1 diabetes displayed increased Class I transcripts and associated mRNA species. These ß-cells also expressed mRNA for Class II and Class II antigen presentation pathway components, but lacked the macrophage marker CD68. Immunohistological study of three independent cohorts of donors with recent-onset type 1 diabetes showed Class II protein and its transcriptional regulator Class II MHC trans-activator protein expressed by a subset of insulin+CD68- ß-cells, specifically found in islets with lymphocytic infiltrates. ß-Cell surface expression of HLA Class II was detected on a portion of CD45-insulin+ ß-cells from donors with type 1 diabetes by immunofluorescence and flow cytometry. Our data demonstrate that pancreatic ß-cells from donors with type 1 diabetes express Class II molecules on selected cells with other key genes in those pathways and inflammation-associated genes. ß-Cell expression of Class II molecules suggests that ß-cells may interact directly with islet-infiltrating CD4+ T cells and may play an immunopathogenic role.

DEBrowser: interactive differential expression analysis and visualization tool for count data.

BACKGROUND: Sequencing data has become a standard measure of diverse cellular activities. For example, gene expression is accurately measured by RNA sequencing (RNA-Seq) libraries, protein-DNA interactions are captured by chromatin immunoprecipitation sequencing (ChIP-Seq), protein-RNA interactions by crosslinking immunoprecipitation sequencing (CLIP-Seq) or RNA immunoprecipitation (RIP-Seq) sequencing, DNA accessibility by assay for transposase-accessible chromatin (ATAC-Seq), DNase or MNase sequencing libraries. The processing of these sequencing techniques involves library-specific approaches. However, in all cases, once the sequencing libraries are processed, the result is a count table specifying the estimated number of reads originating from each genomic locus. Differential analysis to determine which loci have different cellular activity under different conditions starts with the count table and iterates through a cycle of data assessment, preparation and analysis. Such complex analysis often relies on multiple programs and is therefore a challenge for those without programming skills. RESULTS: We developed DEBrowser as an R bioconductor project to interactively visualize every step of the differential analysis, without programming. The application provides a rich and interactive web based graphical user interface built on R's shiny infrastructure. DEBrowser allows users to visualize data with various types of graphs that can be explored further by selecting and re-plotting any desired subset of data. Using the visualization approaches provided, users can determine and correct technical variations such as batch effects and sequencing depth that affect differential analysis. We show DEBrowser's ease of use by reproducing the analysis of two previously published data sets. CONCLUSIONS: DEBrowser is a flexible, intuitive, web-based analysis platform that enables an iterative and interactive analysis of count data without any requirement of programming knowledge.