Search

Supplementary files to support the manuscript, including antibody and cell line mutational information and confirmation/setup information for different experiments.

Melanoma is one of the most serious forms of skin cancer, and its increasing incidence coupled with nonlasting therapeutic options for metastatic disease highlights the need for additional novel approaches for its management. In this study, we determined the potential interactions between polo-like kinase 1 (PLK1, a serine/threonine kinase involved in mitotic regulation) and NOTCH1 (a type I transmembrane protein deciding cell fate during development) in melanoma. Employing an in-house human melanoma tissue microarray (TMA) containing multiple cases of melanomas and benign nevi, coupled with high-throughput, multispectral quantitative fluorescence imaging analysis, we found a positive correlation between PLK1 and NOTCH1 in melanoma. Furthermore, The Cancer Genome Atlas database analysis of patients with melanoma showed an association of higher mRNA levels of PLK1 and NOTCH1 with poor overall, as well as disease-free, survival. Next, utilizing small-molecule inhibitors of PLK1 and NOTCH (BI 6727 and MK-0752, respectively), we found a synergistic antiproliferative response of combined treatment in multiple human melanoma cells. To determine the molecular targets of the overall and synergistic responses of combined PLK1 and NOTCH inhibition, we conducted RNA-sequencing analysis employing a unique regression model with interaction terms. We identified the modulations of several key genes relevant to melanoma progression/metastasis, including MAPK, PI3K, and RAS, as well as some new genes such as Apobec3G, BTK, and FCER1G, which have not been well studied in melanoma. In conclusion, our study demonstrated a synergistic antiproliferative response of concomitant targeting of PLK1 and NOTCH in melanoma, unraveling a potential novel therapeutic approach for detailed preclinical/clinical evaluation.

The RNA-regulatory exosome complex (EC) posttranscriptionally and cotranscriptionally processes and degrades RNAs in a context-dependent manner. Although the EC functions in diverse cell types, its contributions to stem and progenitor cell development are not well understood. Previously, we demonstrated that the transcriptional regulator of erythrocyte development, GATA1, represses EC subunit genes, and the EC maintains erythroid progenitors in vitro. To determine if this mechanism operates in vivo, we used the hematopoietic-specific Vav1-Cre and “conditional by inversion” mouse system to ablate Exosc3, encoding an EC structural subunit. Although Exosc3C/C Cre+ embryos developed normally until embryonic day 14.5, Exosc3 ablation was embryonic lethal and severely reduced erythromyeloid progenitor activity. RNA sequencing analysis of Exosc3-ablated burst-forming unit-erythroid revealed elevated transcripts encoding multiple proapoptotic factors, and the mutant erythroid progenitors exhibited increased apoptosis. We propose that the EC controls an ensemble of apoptosis-regulatory RNAs, thereby promoting erythroid progenitor survival and developmental erythropoiesis in vivo.

Cellular differentiation requires vast remodeling of transcriptomes, and therefore machinery mediating remodeling controls differentiation. Relative to transcriptional mechanisms governing differentiation, post-transcriptional processes are less well understood. As an important post-transcriptional determinant of transcriptomes, the RNA exosome complex (EC) mediates processing and/or degradation of select RNAs. During erythropoiesis, the erythroid transcription factor GATA1 represses EC subunit genes. Depleting EC structural subunits prior to GATA1-mediated repression is deleterious to erythroid progenitor cells. To assess the importance of the EC catalytic subunits Dis3 and Exosc10 in this dynamic process, we asked if these subunits function non-redundantly to control erythropoiesis. Dis3 or Exosc10 depletion in primary murine hematopoietic progenitor cells reduced erythroid progenitors and their progeny, while sparing myeloid cells. Dis3 loss severely compromised erythroid progenitor and erythroblast survival, rendered erythroblasts hypersensitive to apoptosis-inducing stimuli and induced γ-H2AX, indicative of DNA double-stranded breaks. Dis3 loss-of-function phenotypes were more severe than those caused by Exosc10 depletion. We innovated a genetic rescue system to compare human Dis3 with multiple myeloma-associated Dis3 mutants S447R and R750K, and only wild type Dis3 was competent to rescue progenitors. Thus, Dis3 establishes a disease mutation-sensitive, cell type-specific survival mechanism to enable a differentiation program.

Memory T cells underpin vaccine-induced immunity but are not yet fully understood. To distinguish features of memory cells that confer protective immunity, we used single cell transcriptome analysis to compare antigen-specific CD4+T cells recalled to lungs of mice that received a protective or nonprotective subunit vaccine followed by challenge with a fungal pathogen. We unexpectedly found populations specific to protection that expressed a strong type I interferon response signature, whose distinctive transcriptional signature appeared unconventionally dependent on IFN-γ receptor. We also detected a unique population enriched in protection that highly expressed the gene for the natural killer cell marker NKG7. Lastly, we detected differences in TCR gene use and in Th1- and Th17-skewed responses after protective and nonprotective vaccine, respectively, reflecting heterogeneousIfng- andIl17a-expressing populations. Our findings highlight key features of transcriptionally diverse and distinctive antigen-specific T cells associated with protective vaccine-induced immunity.

Island populations repeatedly evolve extreme body sizes, but the genomic basis of this pattern remains largely unknown. To understand how organisms on islands evolve gigantism, we compared genome-wide patterns of gene expression in Gough Island mice, the largest wild house mice in the world, and mainland mice from the WSB/EiJ wild-derived inbred strain. We used RNA-seq to quantify differential gene expression in three key metabolic organs: gonadal adipose depot, hypothalamus, and liver. Between 4,000 and 8,800 genes were significantly differentially expressed across the evaluated organs, representing between 20% and 50% of detected transcripts, with 20% or more of differentially expressed transcripts in each organ exhibiting expression fold changes of at least 2×. A minimum of 73 candidate genes for extreme size evolution, including Irs1 and Lrp1, were identified by considering differential expression jointly with other data sets: 1) genomic positions of published quantitative trait loci for body weight and growth rate, 2) whole-genome sequencing of 16 wild-caught Gough Island mice that revealed fixed single-nucleotide differences between the strains, and 3) publicly available tissue-specific regulatory elements. Additionally, patterns of differential expression across three time points in the liver revealed that Arid5b potentially regulates hundreds of genes. Functional enrichment analyses pointed to cell cycling, mitochondrial function, signaling pathways, inflammatory response, and nutrient metabolism as potential causes of weight accumulation in Gough Island mice. Collectively, our results indicate that extensive gene regulatory evolution in metabolic organs accompanied the rapid evolution of gigantism during the short time house mice have inhabited Gough Island.

The second messenger nucleotide ppGpp dramatically alters gene expression in bacteria to adjust cellular metabolism to nutrient availability. ppGpp binds to two sites on RNA polymerase (RNAP) in Escherichia coli , but it has also been reported to bind to many other proteins. To determine the role of the RNAP binding sites in the genome-wide effects of ppGpp on transcription, we used RNA-seq to analyze transcripts produced in response to elevated ppGpp levels in strains with/without the ppGpp binding sites on RNAP. We examined RNAs rapidly after ppGpp production without an accompanying nutrient starvation. This procedure enriched for direct effects of ppGpp on RNAP rather than for indirect effects on transcription resulting from starvation-induced changes in metabolism or on secondary events from the initial effects on RNAP. The transcriptional responses of all 757 genes identified after 5 minutes of ppGpp induction depended on ppGpp binding to RNAP. Most (>75%) were not reported in earlier studies. The regulated transcripts encode products involved not only in translation but also in many other cellular processes. In vitro transcription analysis of more than 100 promoters from the in vivo dataset identified a large collection of directly regulated promoters, unambiguously demonstrated that most effects of ppGpp on transcription in vivo were direct, and allowed comparison of DNA sequences from inhibited, activated, and unaffected promoter classes. Our analysis greatly expands our understanding of the breadth of the stringent response and suggests promoter sequence features that contribute to the specific effects of ppGpp.

Summary Cell type annotation is a fundamental task in the analysis of single-cell RNA-sequencing data. In this work, we present CellO, a machine learning-based tool for annotating human RNA-seq data with the Cell Ontology. CellO enables accurate and standardized cell type classification of cell clusters by considering the rich hierarchical structure of known cell types. Furthermore, CellO comes pre-trained on a comprehensive data set of human, healthy, untreated primary samples in the Sequence Read Archive. CellO's comprehensive training set enables it to run out of the box on diverse cell types and achieves competitive or even superior performance when compared to existing state-of-the-art methods. Lastly, CellO's linear models are easily interpreted, thereby enabling exploration of cell-type-specific expression signatures across the ontology. To this end, we also present the CellO Viewer: a web application for exploring CellO's models across the ontology., Graphical abstract, Highlights • CellO hierarchically annotates single-cell RNA-seq data using the Cell Ontology • CellO is pre-trained on a comprehensive data set comprising diverse cell types • CellO achieves superior or comparable performance with existing methods • The CellO Viewer is a novel web application for exploring CellO's trained models, Genomics; Classification of Bioinformatical Subject; Genomic Analysis

Tissue damage induces rapid recruitment of leukocytes and changes in the transcriptional landscape that influence wound healing. However, the cell-type specific transcriptional changes that influence leukocyte function and tissue repair have not been well characterized. Here, we employed translating ribosome affinity purification (TRAP) and RNA sequencing, TRAP-seq, in larval zebrafish to identify genes differentially expressed in neutrophils, macrophages, and epithelial cells in response to wounding. We identified the complement pathway and c3a.1, homologous to the C3 component of human complement, as significantly increased in neutrophils in response to wounds. c3a.1−/− zebrafish larvae have impaired neutrophil directed migration to tail wounds with an initial lag in recruitment early after wounding. Moreover, c3a.1−/− zebrafish larvae have impaired recruitment to localized bacterial infections and reduced survival that is, at least in part, neutrophil mediated. Together, our findings support the power of TRAP-seq to identify cell type specific changes in gene expression that influence neutrophil behavior in response to tissue damage.

Melanoma is one of the most serious forms of skin cancer, and its increasing incidence coupled with nonlasting therapeutic options for metastatic disease highlights the need for additional novel approaches for its management. In this study, we determined the potential interactions between polo-like kinase 1 (PLK1, a serine/threonine kinase involved in mitotic regulation) and NOTCH1 (a type I transmembrane protein deciding cell fate during development) in melanoma. Employing an in-house human melanoma tissue microarray (TMA) containing multiple cases of melanomas and benign nevi, coupled with high-throughput, multispectral quantitative fluorescence imaging analysis, we found a positive correlation between PLK1 and NOTCH1 in melanoma. Furthermore, The Cancer Genome Atlas database analysis of patients with melanoma showed an association of higher mRNA levels of PLK1 and NOTCH1 with poor overall, as well as disease-free, survival. Next, utilizing small-molecule inhibitors of PLK1 and NOTCH (BI 6727 and MK-0752, respectively), we found a synergistic antiproliferative response of combined treatment in multiple human melanoma cells. To determine the molecular targets of the overall and synergistic responses of combined PLK1 and NOTCH inhibition, we conducted RNA-sequencing analysis employing a unique regression model with interaction terms. We identified the modulations of several key genes relevant to melanoma progression/metastasis, including MAPK, PI3K, and RAS, as well as some new genes such as Apobec3G, BTK, and FCER1G, which have not been well studied in melanoma. In conclusion, our study demonstrated a synergistic antiproliferative response of concomitant targeting of PLK1 and NOTCH in melanoma, unraveling a potential novel therapeutic approach for detailed preclinical/clinical evaluation.

Following acute injury, leukocytes rapidly infiltrate into tissues. For efficient recruitment, leukocytes must sense and respond to signals from both from the damaged tissue and from one another. However, the cell type specific transcriptional changes that influence leukocyte recruitment and wound healing have not been well characterized. In this study, we performed a large-scale translating ribosome affinity purification (TRAP) and RNA sequencing screen in larval zebrafish to identify genes differentially expressed by neutrophils, macrophages, and epithelial cells in the context of wounding. We identified the complement pathway andc3a.1, homologous to the C3A component of human complement, as significantly increased in neutrophils in response to a wound. We report thatc3a.1−/−zebrafish larvae have impaired neutrophil responses to both tail wounds and localized bacterial infections, as well as increased susceptibility to infection due to a neutrophil-intrinsic function of C3A. We further show that C3A enhances migration of human primary neutrophils to IL-8 and thatc3a.1−/−larvae have impaired neutrophil migrationin vivo, and a decrease in neutrophil directed migration speed early after wounding. Together, our findings suggest a role for C3A in mediating efficient neutrophil migration to damaged tissues and support the power of TRAP to identify cell-specific changes in gene expression associated with wound-associated inflammation.

Motivation The NCBI’s Sequence Read Archive (SRA) promises great biological insight if one could analyze the data in the aggregate; however, the data remain largely underutilized, in part, due to the poor structure of the metadata associated with each sample. The rules governing submissions to the SRA do not dictate a standardized set of terms that should be used to describe the biological samples from which the sequencing data are derived. As a result, the metadata include many synonyms, spelling variants and references to outside sources of information. Furthermore, manual annotation of the data remains intractable due to the large number of samples in the archive. For these reasons, it has been difficult to perform large-scale analyses that study the relationships between biomolecular processes and phenotype across diverse diseases, tissues and cell types present in the SRA. Results We present MetaSRA, a database of normalized SRA human sample-specific metadata following a schema inspired by the metadata organization of the ENCODE project. This schema involves mapping samples to terms in biomedical ontologies, labeling each sample with a sample-type category, and extracting real-valued properties. We automated these tasks via a novel computational pipeline. Availability and implementation The MetaSRA is available at metasra.biostat.wisc.edu via both a searchable web interface and bulk downloads. Software implementing our computational pipeline is available at http://github.com/deweylab/metasra-pipeline Supplementary information Supplementary data are available at Bioinformatics online.

Whole-genome alignment (WGA) is the prediction of evolutionary relationships at the nucleotide level between two or more genomes. It combines aspects of both colinear sequence alignment and gene orthology prediction, and is typically more challenging to address than either of these tasks due to the size and complexity of whole genomes. Despite the difficulty of this problem, numerous methods have been developed for its solution because WGAs are valuable for genome-wide analyses, such as phylogenetic inference, genome annotation, and function prediction. In this chapter, we discuss the meaning and significance of WGA and present an overview of the methods that address it. We also examine the problem of evaluating whole-genome aligners and offer a set of methodological challenges that need to be tackled in order to make the most effective use of our rapidly growing databases of whole genomes.

Publicly available RNA-seq data is routinely used for retrospective analysis to elucidate new biology. Novel transcript discovery enabled by joint examination of large collections of RNA-seq datasets has emerged as one such analysis. Current methods for transcript discovery rely on a ‘2-Step’ approach where the first step encompasses building transcripts from individual datasets, followed by the second step that merges predicted transcripts across datasets. To increase the power of transcript discovery from large collections of RNA-seq datasets, we developed a novel ‘1-Step’ approach named Pooling RNA-seq and Assembling Models (PRAM) that builds transcript models from pooled RNA-seq datasets. We demonstrate in a computational benchmark that ‘1-Step’ outperforms ‘2-Step’ approaches in predicting overall transcript structures and individual splice junctions, while performing competitively in detecting exonic nucleotides. Applying PRAM to 30 human ENCODE RNA-seq datasets identified unannotated transcripts with epigenetic and RAMPAGE signatures similar to those of recently annotated transcripts. In a case study, we discovered and experimentally validated new transcripts through the application of PRAM to mouse hematopoietic RNA-seq datasets. Notably, we uncovered new transcripts that share a differential expression pattern with a neighboring genePik3cgimplicated in human hematopoietic phenotypes, and we provided evidence for the conservation of this relationship in human. PRAM is implemented as an R/Bioconductor package and is available athttps://bioconductor.org/packages/pram.

The second messenger nucleotide ppGpp dramatically alters gene expression in bacteria to adjust cellular metabolism to nutrient availability. ppGpp binds to two sites on RNA polymerase (RNAP) in Escherichia coli, but it has also been reported to bind to many other proteins. To determine the role of the RNAP binding sites in the genome-wide effects of ppGpp on transcription, we used RNA-seq to analyze transcripts produced in response to elevated ppGpp levels in strains with/without the ppGpp binding sites on RNAP. We examined RNAs rapidly after ppGpp production without an accompanying nutrient starvation. This procedure enriched for direct effects of ppGpp on RNAP rather than for indirect effects on transcription resulting from starvation-induced changes in metabolism or on secondary events from the initial effects on RNAP. The transcriptional responses of all 757 genes identified after 5 minutes of ppGpp induction depended on ppGpp binding to RNAP. Most (>75%) were not reported in earlier studies. The regulated transcripts encode products involved not only in translation but also in many other cellular processes. In vitro transcription analysis of more than 100 promoters from the in vivo dataset identified a large collection of directly regulated promoters, unambiguously demonstrated that most effects of ppGpp on transcription in vivo were direct, and allowed comparison of DNA sequences from inhibited, activated, and unaffected promoter classes. Our analysis greatly expands our understanding of the breadth of the stringent response and suggests promoter sequence features that contribute to the specific effects of ppGpp.

Whole-genome alignment (WGA) is the prediction of evolutionary relationships at the nucleotide level between two or more genomes. It combines aspects of both colinear sequence alignment and gene orthology prediction, and is typically more challenging to address than either of these tasks due to the size and complexity of whole genomes. Despite the difficulty of this problem, numerous methods have been developed for its solution because WGAs are valuable for genome-wide analyses, such as phylogenetic inference, genome annotation, and function prediction. In this chapter, we discuss the meaning and significance of WGA and present an overview of the methods that address it. We also examine the problem of evaluating whole-genome aligners and offer a set of methodological challenges that need to be tackled in order to make the most effective use of our rapidly growing databases of whole genomes.

RNA-seq is currently the technology of choice for global measurement of transcript abundances in cells. Despite its successes, isoform-level quantification remains difficult because short RNA-seq reads are often compatible with multiple alternatively spliced isoforms. Existing methods rely heavily on uniquely mapping reads, which are not available for numerous isoforms that lack regions of unique sequence. To improve quantification accuracy in such difficult cases, we developed a novel computational method, prior-enhanced RSEM (pRSEM), which uses a complementary data type in addition to RNA-seq data. We found that ChIP-seq data of RNA polymerase II and histone modifications were particularly informative in this approach. In qRT-PCR validations, pRSEM was shown to be superior than competing methods in estimating relative isoform abundances within or across conditions. Data-driven simulations suggested that pRSEM has a greatly decreased false-positive rate at the expense of a small increase in false-negative rate. In aggregate, our study demonstrates that pRSEM transforms existing capacity to precisely estimate transcript abundances, especially at the isoform level.

Melanoma is one of the most brutal forms of skin cancer, and its increasing incidence coupled with non-lasting therapeutic options for metastatic tumor highlight the need for additional strategies for the management of this neoplasm. Using tissue microarray analysis, we previously found that the expression of polo-like kinase 1 (Plk1, a serine/threonine kinase involved in mitotic regulation) and Notch1 (a type I transmembrane protein deciding cell fate during development) were positively correlated in melanoma (Cancer Res 2018; 78 [13 Suppl]: Abstract nr 2530), and their combined inhibition resulted in a synergistic anti-proliferative response in human melanoma cells (Cancer Res 2019; 79 [13 Suppl]: Abstract nr 302). In this study, to determine the possible mechanisms behind this observed synergism, we used RNA-seq technology to obtain the differential gene expression following treatment of SK-MEL-2 human metastatic melanoma cells with Plk1 inhibitor volasertib (BI6727, 20 nM) and Notch1 inhibitor MK-0752 (100 μM) for 48 h. After data pre-processing by RSEM algorithm, the DESeq2 package was implemented to identify differentially expressed genes (DEGs, |log2-fold change| >= 1, false positive rate ⇐ 0.05) when comparing the individual and combined treatments to vehicle (DMSO), as well as the interaction between volasertib:MK-0752. As a result, we identified 909 DEGs from volasertib treatment, 675 DEGs from MK-0752 treatment, 2142 genes from the combined treatment of volasertib and MK-0752, as well as 304 DEGs from the interaction of volasertib and MK-0752. In addition, employing GOstats and KEGGprofile packages in R programming, we conducted Gene Ontology (GO) and KEGG pathway analysis of the various DEGs. In GO analysis (counts >= 2, p ⇐ 10−5), we identified 202 downregulated GO terms affected by the combined inhibition of Plk1 and Notch1, including metabolism, cell proliferation, and migration. In KEGG pathway analysis, the combined inhibition of Plk1 and Notch1 was found to be associated with downregulation of several pathways shared with single drug treatments, such as PI3K-Akt, extracellular matrix receptor interaction, and protein digestion and absorption, as well as some novel pathways that were only affected by combined treatment, such as MAPK, Ras, and Rap1 pathways. Interestingly, our analysis predicted that the combined inhibition of Plk1 and Notch may make the melanoma cells more sensitive to immune responses. Overall, our data demonstrated that not only does targeting both Plk1 and Notch1 signaling pathways alters multiple melanoma progression pathways, but it may also potentially result in an increased sensitivity to other therapeutic targets, such as immune checkpoint blockade. However, these mechanistic findings need to be validated further in other relevant in vitro and in vivo models. Citation Format: Shengqin Su, Gagan Chhabra, Mary A. Ndiaye, Chandra K. Singh, Colin N. Dewey, Nihal Ahmad. RNA-seq analysis of differential gene expression in melanoma cells after combined inhibition of Plk1 and Notch [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 222.

Thousands of cis-elements in genomes are predicted to have vital functions. Although conservation, activity in surrogate assays, polymorphisms, and disease mutations provide functional clues, deletion from endogenous loci constitutes the gold-standard test. A GATA-2-binding, Gata2 intronic cis-element (+9.5) required for hematopoietic stem cell genesis in mice is mutated in a human immunodeficiency syndrome. Because +9.5 is the only cis-element known to mediate stem cell genesis, we devised a strategy to identify functionally comparable enhancers ("+9.5-like") genome-wide. Gene editing revealed +9.5-like activity to mediate GATA-2 occupancy, chromatin opening, and transcriptional activation. A +9.5-like element resided in Samd14, which encodes a protein of unknown function. Samd14 increased hematopoietic progenitor levels/activity and promoted signaling by a pathway vital for hematopoietic stem/progenitor cell regulation (stem cell factor/c-Kit), and c-Kit rescued Samd14 loss-of-function phenotypes. Thus, the hematopoietic stem/progenitor cell cistrome revealed a mediator of a signaling pathway that has broad importance for stem/progenitor cell biology.

Motivation: With improvements in next-generation sequencing technologies and reductions in price, ordered RNA-seq experiments are becoming common. Of primary interest in these experiments is identifying genes that are changing over time or space, for example, and then characterizing the specific expression changes. A number of robust statistical methods are available to identify genes showing differential expression among multiple conditions, but most assume conditions are exchangeable and thereby sacrifice power and precision when applied to ordered data. Results: We propose an empirical Bayes mixture modeling approach called EBSeq-HMM. In EBSeq-HMM, an auto-regressive hidden Markov model is implemented to accommodate dependence in gene expression across ordered conditions. As demonstrated in simulation and case studies, the output proves useful in identifying differentially expressed genes and in specifying gene-specific expression paths. EBSeq-HMM may also be used for inference regarding isoform expression. Availability and implementation: An R package containing examples and sample datasets is available at Bioconductor. Contact: kendzior@biostat.wisc.edu Supplementary information: Supplementary data are available at Bioinformatics online.

An enhancer with amalgamated E-box and GATA motifs (+9.5) controls expression of the regulator of hematopoiesis GATA-2. While similar GATA-2-occupied elements are common in the genome, occupancy does not predict function, and GATA-2-dependent genetic networks are incompletely defined. A "+9.5-like" element resides in an intron of Samd14 (Samd14-Enh) encoding a sterile alpha motif (SAM) domain protein. Deletion of Samd14-Enh in mice strongly decreased Samd14 expression in bone marrow and spleen. Although steady-state hematopoiesis was normal, Samd14-Enh-/- mice died in response to severe anemia. Samd14-Enh stimulated stem cell factor/c-Kit signaling, which promotes erythrocyte regeneration. Anemia activated Samd14-Enh by inducing enhancer components and enhancer chromatin accessibility. Thus, a GATA-2/anemia-regulated enhancer controls expression of an SAM domain protein that confers survival in anemia. We propose that Samd14-Enh and an ensemble of anemia-responsive enhancers are essential for erythrocyte regeneration in stress erythropoiesis, a vital process in pathologies, including β-thalassemia, myelodysplastic syndrome, and viral infection.

The vertebrate retina develops in close proximity to the forebrain and neural crest-derived cartilages of the face and jaw. Coloboma, a congenital eye malformation, is associated with aberrant forebrain development (holoprosencephaly) and with craniofacial defects (frontonasal dysplasia) in humans, suggesting a critical role for cross-lineage interactions during retinal morphogenesis. ZIC2, a zinc-finger transcription factor, is linked to human holoprosencephaly. We have previously used morpholino assays to show zebrafish zic2 functions in the developing forebrain, retina and craniofacial cartilage. We now report that zebrafish with genetic lesions in zebrafish zic2 orthologs, zic2a and zic2b, develop with retinal coloboma and craniofacial anomalies. We demonstrate a requirement for zic2 in restricting pax2a expression and show evidence that zic2 function limits Hh signaling. RNA-seq transcriptome analysis identified an early requirement for zic2 in periocular neural crest as an activator of alx1, a transcription factor with essential roles in craniofacial and ocular morphogenesis in human and zebrafish. Collectively, these data establish zic2 mutant zebrafish as a powerful new genetic model for in-depth dissection of cell interactions and genetic controls during craniofacial complex development.

MotivationThe NCBI’s Sequence Read Archive (SRA) promises great biological insight if one could analyze the data in the aggregate; however, the data remain largely underutilized, in part, due to the poor structure of the metadata associated with each sample. The rules governing submissions to the SRA do not dictate a standardized set of terms that should be used to describe the biological samples from which the sequencing data are derived. As a result, the metadata include many synonyms, spelling variants, and references to outside sources of information. Furthermore, manual annotation of the data remains intractable due to the large number of samples in the archive. For these reasons, it has been difficult to perform large-scale analyses that study the relationships between biomolecular processes and phenotype across diverse diseases, tissues, and cell types present in the SRA.ResultsWe present MetaSRA, a database of normalized SRA sample-specific metadata following a schema inspired by the metadata organization of the ENCODE project. This schema involves mapping samples to terms in biomedical ontologies, labeling each sample with a sample-type category, and extracting real-valued properties. We automated these tasks via a novel computational pipeline.AvailabilityThe MetaSRA database is available at http://deweylab.biostat.wisc.edu/metasra. Software implementing our computational pipeline is available at https://github.com/deweylab/metasra-pipeline.Contactcdewey@biostat.wisc.edu

Cis-element requirement for the emergence of HSCs in the AGM and for hemogenic endothelium to generate HSC-containing c-Kit+ cell clusters., The generation of hematopoietic stem cells (HSCs) from hemogenic endothelium within the aorta, gonad, mesonephros (AGM) region of the mammalian embryo is crucial for development of the adult hematopoietic system. We described a deletion of a Gata2 cis-element (+9.5) that depletes fetal liver HSCs, is lethal at E13–14 of embryogenesis, and is mutated in an immunodeficiency that progresses to myelodysplasia/leukemia. Here, we demonstrate that the +9.5 element enhances Gata2 expression and is required to generate long-term repopulating HSCs in the AGM. Deletion of the +9.5 element abrogated the capacity of hemogenic endothelium to generate HSC-containing clusters in the aorta. Genomic analyses indicated that the +9.5 element regulated a rich ensemble of genes that control hemogenic endothelium and HSCs, as well as genes not implicated in hematopoiesis. These results reveal a mechanism that controls stem cell emergence from hemogenic endothelium to establish the adult hematopoietic system.

The Xenopus Cripto-1 protein is confined to the cells of the animal hemisphere during early embryogenesis where it regulates the formation of anterior structures. Cripto-1 protein accumulates only in animal cells because cripto-1 mRNA in cells of the vegetal hemisphere is translationally repressed. Here, we show that the RNA binding protein, Bicaudal-C (Bic-C), functioned directly in this vegetal cell-specific repression. While Bic-C protein is normally confined to vegetal cells, ectopic expression of Bic-C in animal cells repressed a cripto-1 mRNA reporter and associated with endogenous cripto-1 mRNA. Repression by Bic-C required its N-terminal domain, comprised of multiple KH motifs, for specific binding to relevant control elements within the cripto-1 mRNA and a functionally separable C-terminal translation repression domain. Bic-C-mediated repression required the 5′ CAP and translation initiation factors, but not a poly(A) tail or the conserved SAM domain within Bic-C. Bic-C-directed immunoprecipitation followed by deep sequencing of associated mRNAs identified multiple Bic-C-regulated mRNA targets, including cripto-1 mRNA, providing new insights and tools for understanding the role of Bic-C in vertebrate development.

Motivation: Alternative splicing and other processes that allow for different transcripts to be derived from the same gene are significant forces in the eukaryotic cell. RNA-Seq is a promising technology for analyzing alternative transcripts, as it does not require prior knowledge of transcript structures or genome sequences. However, analysis of RNA-Seq data in the presence of genes with large numbers of alternative transcripts is currently challenging due to efficiency, identifiability and representation issues. Results: We present RNA-Seq models and associated inference algorithms based on the concept of probabilistic splice graphs, which alleviate these issues. We prove that our models are often identifiable and demonstrate that our inference methods for quantification and differential processing detection are efficient and accurate. Availability: Software implementing our methods is available at http://deweylab.biostat.wisc.edu/psginfer. Contact: ude.csiw.tatsoib@yewedc Supplementary information: Supplementary data are available at Bioinformatics online.

Titin, a sarcomeric protein expressed primarily in striated muscles, is responsible for maintaining the structure and biomechanical properties of muscle cells. Cardiac titin undergoes developmental size reduction from 3.7 megadaltons in neonates to primarily 2.97 megadaltons in the adult. This size reduction results from gradually increased exon skipping between exons 50 and 219 of titin mRNA. Our previous study reported that Rbm20 is the splicing factor responsible for this process. In this work, we investigated its molecular mechanism. We demonstrate that Rbm20 mediates exon skipping by binding to titin pre-mRNA to repress the splicing of some regions; the exons/introns in these Rbm20-repressed regions are ultimately skipped. Rbm20 was also found to mediate intron retention and exon shuffling. The two Rbm20 speckles found in nuclei from muscle tissues were identified as aggregates of Rbm20 protein on the partially processed titin pre-mRNAs. Cooperative repression and alternative 3' splice site selection were found to be used by Rbm20 to skip different subsets of titin exons, and the splicing pathway selected depended on the ratio of Rbm20 to other splicing factors that vary with tissue type and developmental age.

This report of independent genome sequences of two natural populations of Drosophila melanogaster (37 from North America and 6 from Africa) provides unique insight into forces shaping genomic polymorphism and divergence. Evidence of interactions between natural selection and genetic linkage is abundant not only in centromere- and telomere-proximal regions, but also throughout the euchromatic arms. Linkage disequilibrium, which decays within 1 kbp, exhibits a strong bias toward coupling of the more frequent alleles and provides a high-resolution map of recombination rate. The juxtaposition of population genetics statistics in small genomic windows with gene structures and chromatin states yields a rich, high-resolution annotation, including the following: (1) 5′- and 3′-UTRs are enriched for regions of reduced polymorphism relative to lineage-specific divergence; (2) exons overlap with windows of excess relative polymorphism; (3) epigenetic marks associated with active transcription initiation sites overlap with regions of reduced relative polymorphism and relatively reduced estimates of the rate of recombination; (4) the rate of adaptive nonsynonymous fixation increases with the rate of crossing over per base pair; and (5) both duplications and deletions are enriched near origins of replication and their density correlates negatively with the rate of crossing over. Available demographic models of X and autosome descent cannot account for the increased divergence on the X and loss of diversity associated with the out-of-Africa migration. Comparison of the variation among these genomes to variation among genomes from D. simulans suggests that many targets of directional selection are shared between these species.

Motivation: RNA-Seq is a promising new technology for accurately measuring gene expression levels. Expression estimation with RNA-Seq requires the mapping of relatively short sequencing reads to a reference genome or transcript set. Because reads are generally shorter than transcripts from which they are derived, a single read may map to multiple genes and isoforms, complicating expression analyses. Previous computational methods either discard reads that map to multiple locations or allocate them to genes heuristically. Results: We present a generative statistical model and associated inference methods that handle read mapping uncertainty in a principled manner. Through simulations parameterized by real RNA-Seq data, we show that our method is more accurate than previous methods. Our improved accuracy is the result of handling read mapping uncertainty with a statistical model and the estimation of gene expression levels as the sum of isoform expression levels. Unlike previous methods, our method is capable of modeling non-uniform read distributions. Simulations with our method indicate that a read length of 20–25 bases is optimal for gene-level expression estimation from mouse and maize RNA-Seq data when sequencing throughput is fixed. Availability: An initial C++ implementation of our method that was used for the results presented in this article is available at http://deweylab.biostat.wisc.edu/rsem. Contact: cdewey@biostat.wisc.edu Supplementary information: Supplementary data are available at Bioinformatics on

Non-coding DNA elements differentially control stem and progenitor cell transitions required for development., Cis-element encyclopedias provide information on phenotypic diversity and disease mechanisms. Although cis-element polymorphisms and mutations are instructive, deciphering function remains challenging. Mutation of an intronic GATA motif (+9.5) in GATA2, encoding a master regulator of hematopoiesis, underlies an immunodeficiency associated with myelodysplastic syndrome (MDS) and acute myeloid leukemia (AML). Whereas an inversion relocalizes another GATA2 cis-element (−77) to the proto-oncogene EVI1, inducing EVI1 expression and AML, whether this reflects ectopic or physiological activity is unknown. We describe a mouse strain that decouples −77 function from proto-oncogene deregulation. The −77−/− mice exhibited a novel phenotypic constellation including late embryonic lethality and anemia. The −77 established a vital sector of the myeloid progenitor transcriptome, conferring multipotentiality. Unlike the +9.5−/− embryos, hematopoietic stem cell genesis was unaffected in −77−/− embryos. These results illustrate a paradigm in which cis-elements in a locus differentially control stem and progenitor cell transitions, and therefore the individual cis-element alterations cause unique and overlapping disease phenotypes.