The mouse Gene Expression Database (GXD): 2017 update.

The Gene Expression Database (GXD; www.informatics.jax.org/expression.shtml) is an extensive and well-curated community resource of mouse developmental expression information. Through curation of the scientific literature and by collaborations with large-scale expression projects, GXD collects and integrates data from RNA in situ hybridization, immunohistochemistry, RT-PCR, northern blot and western blot experiments. Expression data from both wild-type and mutant mice are included. The expression data are combined with genetic and phenotypic data in Mouse Genome Informatics (MGI) and made readily accessible to many types of database searches. At present, GXD includes over 1.5 million expression results and more than 300 000 images, all annotated with detailed and standardized metadata. Since our last report in 2014, we have added a large amount of data, we have enhanced data and database infrastructure, and we have implemented many new search and display features. Interface enhancements include: a new Mouse Developmental Anatomy Browser; interactive tissue-by-developmental stage and tissue-by-gene matrix views; capabilities to filter and sort expression data summaries; a batch search utility; gene-based expression overviews; and links to expression data from other species.

A customized and versatile high-density genotyping array for the mouse.

We designed a high-density mouse genotyping array containing 623,124 single-nucleotide polymorphisms that captures the known genetic variation present in the laboratory mouse. The array also contains 916,269 invariant genomic probes targeted to functional elements and regions known to harbor segmental duplications. The array opens the door to the characterization of genetic diversity, copy-number variation, allele-specific gene expression and DNA methylation, and will extend the successes of human genome-wide association studies to the mouse.

Klf9 is a key feedforward regulator of the transcriptomic response to glucocorticoid receptor activity.

The zebrafish has recently emerged as a model system for investigating the developmental roles of glucocorticoid signaling and the mechanisms underlying glucocorticoid-induced developmental programming. To assess the role of the Glucocorticoid Receptor (GR) in such programming, we used CRISPR-Cas9 to produce a new frameshift mutation, GR369-, which eliminates all potential in-frame initiation codons upstream of the DNA binding domain. Using RNA-seq to ask how this mutation affects the larval transcriptome under both normal conditions and with chronic cortisol treatment, we find that GR mediates most of the effects of the treatment, and paradoxically, that the transcriptome of cortisol-treated larvae is more like that of larvae lacking a GR than that of larvae with a GR, suggesting that the cortisol-treated larvae develop GR resistance. The one transcriptional regulator that was both underexpressed in GR369- larvae and consistently overexpressed in cortisol-treated larvae was klf9. We therefore used CRISPR-Cas9-mediated mutation of klf9 and RNA-seq to assess Klf9-dependent gene expression in both normal and cortisol-treated larvae. Our results indicate that Klf9 contributes significantly to the transcriptomic response to chronic cortisol exposure, mediating the upregulation of proinflammatory genes that we reported previously.

Position-dependent motif characterization using non-negative matrix factorization.

MOTIVATION: Cis-acting regulatory elements are frequently constrained by both sequence content and positioning relative to a functional site, such as a splice or polyadenylation site. We describe an approach to regulatory motif analysis based on non-negative matrix factorization (NMF). Whereas existing pattern recognition algorithms commonly focus primarily on sequence content, our method simultaneously characterizes both positioning and sequence content of putative motifs. RESULTS: Tests on artificially generated sequences show that NMF can faithfully reproduce both positioning and content of test motifs. We show how the variation of the residual sum of squares can be used to give a robust estimate of the number of motifs or patterns in a sequence set. Our analysis distinguishes multiple motifs with significant overlap in sequence content and/or positioning. Finally, we demonstrate the use of the NMF approach through characterization of biologically interesting datasets. Specifically, an analysis of mRNA 3'-processing (cleavage and polyadenylation) sites from a broad range of higher eukaryotes reveals a conserved core pattern of three elements.

Systematic variation in mRNA 3'-processing signals during mouse spermatogenesis.

Gene expression and processing during mouse male germ cell maturation (spermatogenesis) is highly specialized. Previous reports have suggested that there is a high incidence of alternative 3'-processing in male germ cell mRNAs, including reduced usage of the canonical polyadenylation signal, AAUAAA. We used EST libraries generated from mouse testicular cells to identify 3'-processing sites used at various stages of spermatogenesis (spermatogonia, spermatocytes and round spermatids) and testicular somatic Sertoli cells. We assessed differences in 3'-processing characteristics in the testicular samples, compared to control sets of widely used 3'-processing sites. Using a new method for comparison of degenerate regulatory elements between sequence samples, we identified significant changes in the use of putative 3'-processing regulatory sequence elements in all spermatogenic cell types. In addition, we observed a trend towards truncated 3'-untranslated regions (3'-UTRs), with the most significant differences apparent in round spermatids. In contrast, Sertoli cells displayed a much smaller trend towards 3'-UTR truncation and no significant difference in 3'-processing regulatory sequences. Finally, we identified a number of genes encoding mRNAs that were specifically subject to alternative 3'-processing during meiosis and postmeiotic development. Our results highlight developmental differences in polyadenylation site choice and in the elements that likely control them during spermatogenesis.

CGDSNPdb: a database resource for error-checked and imputed mouse SNPs.

The Center for Genome Dynamics Single Nucleotide Polymorphism Database (CGDSNPdb) is an open-source value-added database with more than nine million mouse single nucleotide polymorphisms (SNPs), drawn from multiple sources, with genotypes assigned to multiple inbred strains of laboratory mice. All SNPs are checked for accuracy and annotated for properties specific to the SNP as well as those implied by changes to overlapping protein-coding genes. CGDSNPdb serves as the primary interface to two unique data sets, the 'imputed genotype resource' in which a Hidden Markov Model was used to assess local haplotypes and the most probable base assignment at several million genomic loci in tens of strains of mice, and the Affymetrix Mouse Diversity Genotyping Array, a high density microarray with over 600,000 SNPs and over 900,000 invariant genomic probes. CGDSNPdb is accessible online through either a web-based query tool or a MySQL public login. Database URL: http://cgd.jax.org/cgdsnpdb/

C. elegans sequences that control trans-splicing and operon pre-mRNA processing.

Many mRNAs in Caenorhabditis elegans are generated through a trans-splicing reaction that adds one of two classes of spliced leader RNA to an independently transcribed pre-mRNA. SL1 leaders are spliced mostly to pre-mRNAs from genes with outrons, intron-like sequences at the 5'-ends of the pre-mRNAs. In contrast, SL2 leaders are nearly exclusively trans-spliced to genes that occur downstream in polycistronic pre-mRNAs produced from operons. Operon pre-mRNA processing requires separation into individual transcripts, which is accomplished by 3'-processing of upstream genes and spliced leader trans-splicing to the downstream genes. We used a novel computational analysis, based on nonnegative matrix factorization, to identify and characterize significant differences in the cis-acting sequence elements that differentiate various types of functional site, including internal versus terminal 3'-processing sites, and SL1 versus SL2 trans-splicing sites. We describe several key elements, including the U-rich (Ur) element that couples 3'-processing with SL2 trans-splicing, and a novel outron (Ou) element that occurs upstream of SL1 trans-splicing sites. Finally, we present models of the distinct classes of trans-splicing reaction, including SL1 trans-splicing at the outron, SL2 trans-splicing in standard operons, competitive SL1-SL2 trans-splicing in operons with large intergenic separation, and SL1 trans-splicing in SL1-type operons, which have no intergenic separation.