Assessment of fight outcome is needed to activate socially driven transcriptional changes in the zebrafish brain

Group living animals must be able to express different behavior profiles depending on their social status. Therefore, the same genotype may translate into different behavioral phenotypes through socially driven differential gene expression. However, how social information is translated into a neurogenomic response and what are the specific cues in a social interaction that signal a change in social status are questions that have remained unanswered. Here, we show for the first time, to our knowledge, that the switch between status-specific neurogenomic states relies on the assessment of fight outcome rather than just on self- or opponent-only assessment of fighting ability. For this purpose, we manipulated the perception of fight outcome in male zebrafish and measured its impact on the brain transcriptome using a zebrafish whole genome gene chip. Males fought either a real opponent, and a winner and a loser were identified, or their own image on a mirror, in which case, despite expressing aggressive behavior, males did not experience either a victory or a defeat. Massive changes in the brain transcriptome were observed in real opponent fighters, with losers displaying both a higher number of differentially expressed genes and of coexpressed gene modules than winners. In contrast, mirror fighters expressed a neurogenomic state similar to that of noninteracting fish. The genes that responded to fight outcome included immediate early genes and genes involved in neuroplasticity and epigenetic modifications. These results indicate that, even in cognitively simple organisms such as zebrafish, neurogenomic responses underlying changes in social status rely on mutual assessment of fighting ability.

Using the R package Shiny to create web applications that facilitate quantitative gene expression analysis

Thesis (Master's)--University of Washington, 2016-06

Statistical Hurdle Models for Single Cell Gene Expression: Differential Expression and Graphical Modeling

Thesis (Ph.D.)--University of Washington, 2016-06

Phosphoregulators: Protein kinases and protein phosphatases of mouse

With the completion of the human and mouse genome sequences, the task now turns to identifying their encoded transcripts and assigning gene function. In this study, we have undertaken a computational approach to identify and classify all of the protein kinases and phosphatases present in the mouse gene complement. A nonredundant set of these sequences was produced by mining Ensembl gene predictions and publicly available cDNA sequences with a panel of InterPro domains. This approach identified 561 candidate protein kinases and 162 candidate protein phosphatases. This cohort was then analyzed using TribeMCL protein sequence similarity clustering followed by CLUSTALV alignment and hierarchical tree generation. This approach allowed us to (1) distinguish between true members of the protein kinase and phosphatase families and enzymes of related biochemistry, (2) determine the structure of the families, and (3) suggest functions for previously uncharacterized members. The classifications obtained by this approach were in good agreement with previous schemes and allowed us to demonstrate domain associations with a number of clusters. Finally, we comment on the complementary nature of cDNA and genome-based gene detection and the impact of the FANTOM2 transcriptome project.

Continued discovery of transcriptional units expressed in cells of the mouse mononuclear phagocyte lineage

The current RIKEN transcript set represents a significant proportion of the mouse transcriptome but transcripts expressed in the innate and acquired immune systems are poorly represented. In the present study we have assessed the complexity of the transcriptome expressed in mouse macrophages before and after treatment with lipopolysaccharide, a global regulator of macrophage gene expression, using existing RIKEN 19K arrays. By comparison to array profiles of other cells and tissues, we identify a large set of macrophage-enriched genes, many of which have obvious functions in endocytosis and phagocytosis. In addition, a significant number of LPS-inducible genes were identified. The data suggest that macrophages are a complex source of mRNA for transcriptome studies. To assess complexity and identify additional macrophage expressed genes, cDNA libraries were created from purified populations of macrophage and dendritic cells, a functionally related cell type. Sequence analysis revealed a high incidence of novel mRNAs within these cDNA libraries. These studies provide insights into the depths of transcriptional complexity still untapped amongst products of inducible genes, and identify macrophage and dendritic cell populations as a starting point for sampling the inducible mammalian transcriptome.

The functional genomics of noncoding RNA

Large numbers of noncoding RNA transcripts (ncRNAS) are being revealed by complementary DNA cloning and genome tiling array studies in animals. The big and as yet largely unanswered question is whether these transcripts are relevant. A paper by Willingham et al. shows the way forward by developing a strategy for large-scale functional screening of ncRNAs, involving small interfering RNA knockdowns in cell-based screens, which identified a previously unidentified ncRNA repressor of the transcription factor NFAT. It appears likely that ncRNAs constitute a critical hidden layer of gene regulation in complex organisms, the understanding of which requires new approaches in functional genomics.

In vitro and in silico analysis of signal peptides from the human blood fluke, Schistosoma mansoni

Proteins secreted by and anchored on the surfaces of parasites are in intimate contact with host tissues. The transcriptome of infective cercariae of the blood fluke, Schistosoma mansoni, was screened using signal sequence trap to isolate cDNAs encoding predicted proteins with an N-terminal signal peptide. Twenty cDNA fragments were identified, most of which contained predicted signal peptides or transmembrane regions, including a novel putative seven-transmembrane receptor and a membrane-associated mitogen-activated protein kinase. The developmental expression pattern within different life-cycle stages ranged from ubiquitous to a transcript that was highly upregulated in the cercaria. A bioinformatics-based comparison of 100 signal peptides from each of schistosomes, humans, a parasitic nematode and Escherichia coli showed that differences in the sequence composition of signal peptides, notably the residues flanking the predicted cleavage site, might account for the negative bias exhibited in the processing of schistosome signal peptides in mammalian cells. (c) 2005 Federation of European Microbiological Societies. Published by Elsevier B.V. All rights reserved.

Transcription factor Tfec contributes to the IL-4-inducible expression of a small group of genes in mouse macrophages including the granulocyte colony-stimulating factor receptor

Expression of the mouse transcription factor EC (Tfec) is restricted to the myeloid compartment, suggesting a function for Tfec in the development or function of these cells. However, mice lacking Tfec develop normally, indicating a redundant role for Tfec in myeloid cell development. We now report that Tfec is specifically induced in bone marrow-derived macrophages upon stimulation with the Th2 cytokines, IL-4 and IL-13, or LPS. LPS induced a rapid and transient up-regulation of Tfec mRNA expression and promoter activity, which was dependent on a functional NF-kappa B site. IL-4, however, induced a rapid, but long-lasting, increase in Tfec mRNA, which, in contrast to LPS stimulation, also resulted in detectable levels of Tfec protein. IL-4-induced transcription of Tfec was absent in macrophages lacking Stat6, and its promoter depended on two functional Stat6-binding sites. A global comparison of IL-4-induced genes in both wild-type and Tfec mutant macrophages revealed a surprisingly mild phenotype with only a few genes affected by Tfec deficiency. These included the G-CSFR (Csf3r) gene that was strongly up-regulated by IL-4 in wild-type macrophages and, to a lesser extent, in Tfec mutant macrophages. Our study also provides a general definition of the transcriptome in alternatively activated mouse macrophages and identifies a large number of novel genes characterizing this cell type.

Rapid evolution of noncoding RNAs: lack of conservation does not mean lack of function

The mammalian transcriptome contains many nonprotein-coding RNAs (ncRNAs), but most of these are of unclear significance and lack strong sequence conservation, prompting suggestions that they might be non-functional. However, certain long functional ncRNAs such as Air and Xist are also poorly conserved. In this article, we systematically analyzed the conservation of several groups of functional ncRNAs, including miRNAs, snoRNAs and longer ncRNAs whose function has been either documented or confidently predicted. As expected, miRNAs and snoRNAs were highly conserved. By contrast, the longer functional non-micro, non-sno ncRNAs were much less conserved with many displaying rapid sequence evolution. Our findings suggest that longer ncRNAs are under the influence of different evolutionary constraints and that the lack of conservation displayed by the thousands of candidate ncRNAs does not necessarily signify an absence of function.

Transcriptional profile reveals altered hepatic lipid and cholesterol metabolism in hyposulfatemic NaS1 null mice

Sulfate plays an essential role in human growth and development, and its circulating levels are maintained by the renal Na+-SO42- cotransporter, NaS1. We previously generated a NaS1 knockout ( Nas1(-/-)) mouse, an animal model for hyposulfatemia, that exhibits reduced growth and liver abnormalities including hepatomegaly. In this study, we investigated the hepatic gene expression profile of Nas1(-/-) mice using oligonucleotide microarrays. The mRNA expression levels of 92 genes with known functional roles in metabolism, cell signaling, cell defense, immune response, cell structure, transcription, or protein synthesis were increased ( n = 51) or decreased ( n = 41) in Nas1(-/-) mice when compared with Nas1(-/-) mice. The most upregulated transcript levels in Nas1(-/-) mice were found for the sulfotransferase genes, Sult3a1 ( approximate to 500% increase) and Sult2a2 ( 100% increase), whereas the metallothionein-1 gene, Mt1, was among the most downregulated genes ( 70% decrease). Several genes involved in lipid and cholesterol metabolism, including Scd1, Acly, Gpam, Elov16, Acsl5, Mvd, Insig1, and Apoa4, were found to be upregulated ( >= 30% increase) in Nas1(+/+) mice. In addition, Nas1(+/+) mice exhibited increased levels of hepatic lipid ( approximate to 16% increase), serum cholesterol ( approximate to 20% increase), and low-density lipoprotein ( approximate to 100% increase) and reduced hepatic glycogen ( approximate to 50% decrease) levels. In conclusion, these data suggest an altered lipid and cholesterol metabolism in the hyposulfatemic Nas1(-/-) mouse and provide new insights into the metabolic state of the liver in Nas1(-/-) mice.

Differential use of signal peptides and membrane domains is a common occurrence in the protein output of transcriptional units

Membrane organization describes the orientation of a protein with respect to the membrane and can be determined by the presence, or absence, and organization within the protein sequence of two features: endoplasmic reticulum signal peptides and alpha-helical transmembrane domains. These features allow protein sequences to be classified into one of five membrane organization categories: soluble intracellular proteins, soluble secreted proteins, type I membrane proteins, type II membrane proteins, and multi- spanning membrane proteins. Generation of protein isoforms with variable membrane organizations can change a protein's subcellular localization or association with the membrane. Application of MemO, a membrane organization annotation pipeline, to the FANTOM3 Isoform Protein Sequence mouse protein set revealed that within the 8,032 transcriptional units ( TUs) with multiple protein isoforms, 573 had variation in their use of signal peptides, 1,527 had variation in their use of transmembrane domains, and 615 generated protein isoforms from distinct membrane organization classes. The mechanisms underlying these transcript variations were analyzed. While TUs were identified encoding all pairwise combinations of membrane organization categories, the most common was conversion of membrane proteins to soluble proteins. Observed within our highconfidence set were 156 TUs predicted to generate both extracellular soluble and membrane proteins, and 217 TUs generating both intracellular soluble and membrane proteins. The differential use of endoplasmic reticulum signal peptides and transmembrane domains is a common occurrence within the variable protein output of TUs. The generation of protein isoforms that are targeted to multiple subcellular locations represents a major functional consequence of transcript variation within the mouse transcriptome.

Experimental validation of the regulated expression of large numbers of non-coding RNAs from the mouse genome

Recent large-scale analyses of mainly full-length cDNA libraries generated from a variety of mouse tissues indicated that almost half of all representative cloned sequences did flat contain ail apparent protein-coding sequence, and were putatively derived from non-protein-coding RNA (ncRNA) genes. However, many of these clones were singletons and the majority were unspliced, raising the possibility that they may be derived from genomic DNA or unprocessed pre-rnRNA contamination during library construction, or alternatively represent nonspecific transcriptional noise. Here we Show, using reverse transcriptase-dependent PCR, microarray, and Northern blot analyses, that many of these clones were derived from genuine transcripts Of unknown function whose expression appears to be regulated. The ncRNA transcripts have larger exons and fewer introns than protein-coding transcripts. Analysis of the genomic landscape around these sequences indicates that some cDNA clones were produced not from terminal poly(A) tracts but internal priming sites within longer transcripts, only a minority of which is encompassed by known genes. A significant proportion of these transcripts exhibit tissue-specific expression patterns, as well as dynamic changes in their expression in macrophages following lipopolysaccharide Stimulation. Taken together, the data provide strong support for the conclusion that ncRNAs are an important, regulated component of the mammalian transcriptome.

Non-coding RNAs in the nervous system

Increasing evidence suggests that the development and function of the nervous system is heavily dependent on RNA editing and the intricate spatiotemporal expression of a wide repertoire of non-coding RNAs, including micro RNAs, small nucleolar RNAs and longer non-coding RNAs. Non-coding RNAs may provide the key to understanding the multi-tiered links between neural development, nervous system function, and neurological diseases.

Transcript annotation in FANTOM3: Mouse gene catalog based on physical cDNAs

T he international FANTOM consortium aims to produce a comprehensive picture of the mammalian transcriptome, based upon an extensive cDNA collection and functional annotation of full-length enriched cDNAs. The previous dataset, FANTOM(2), comprised 60,770 full- length enriched cDNAs. Functional annotation revealed that this cDNA dataset contained only about half of the estimated number of mouse protein- coding genes, indicating that a number of cDNAs still remained to be collected and identified. To pursue the complete gene catalog that covers all predicted mouse genes, cloning and sequencing of full- length enriched cDNAs has been continued since FANTOM2. In FANTOM3, 42,031 newly isolated cDNAs were subjected to functional annotation, and the annotation of 4,347 FANTOM2 cDNAs was updated. To accomplish accurate functional annotation, we improved our automated annotation pipeline by introducing new coding sequence prediction programs and developed a Web- based annotation interface for simplifying the annotation procedures to reduce manual annotation errors. Automated coding sequence and function prediction was followed with manual curation and review by expert curators. A total of 102,801 full- length enriched mouse cDNAs were annotated. Out of 102,801 transcripts, 56,722 were functionally annotated as protein coding ( including partial or truncated transcripts), providing to our knowledge the greatest current coverage of the mouse proteome by full- length cDNAs. The total number of distinct non- protein- coding transcripts increased to 34,030. The FANTOM3 annotation system, consisting of automated computational prediction, manual curation, and. nal expert curation, facilitated the comprehensive characterization of the mouse transcriptome, and could be applied to the transcriptomes of other species.

The Pineapple EST sequencing and microarray project

We have initiated an EST sequencing project to survey a range of expressed sequences from green fruit, yellow fruit, roots, and root-knot nematode infected root/gall tissues. In total, 5681 edited EST sequences were retrieved. Clone redundancy was high in the fruit libraries, with the combined fruit 1548 clone sequences clustering into just 634 contigs comprising 191 consensus sequences and 443 singletons. Half of all fruit EST clone sequences clustered within approximately 14 and 9% of contigs from green unripe and yellow ripe libraries respectively, indicating that a small subset of genes dominates the majority of the transcriptome. The root and root/gall libraries had lower levels of redundancy than the fruit libraries. Half of the root/gall ESTs clustered within approximately 40% of all contigs, indicating the roots possess a more complex transcriptome. Contig assembly and cluster analysis revealed major differences in the abundant gene sequences expressed between the unripe green and the ripe yellow fruit tissues, or gene sequences expressed between the weeks 1-4 and weeks 5-10 nematode infected gall vascular cylinder libraries.