VEGF-C/VEGFR-3 signaling plays a central role in lymphatic development, regulating the budding of lymphatic progenitor cells from embryonic veins and maintaining the expression of PROX1 during later developmental stages. However, how VEGFR-3 activation translates into target gene expression is still not completely understood. We used cap analysis of gene expression (CAGE) RNA sequencing to characterize the transcriptional changes invoked by VEGF-C in LECs and to identify the transcription factors (TFs) involved. We found that MAFB, a TF involved in differentiation of various cell types, is rapidly induced and activated by VEGF-C. MAFB induced expression of PROX1 as well as other TFs and markers of differentiated LECs, indicating a role in the maintenance of the mature LEC phenotype. Correspondingly, E14.5 Mafb(-/-) embryos showed impaired lymphatic patterning in the skin. This suggests that MAFB is an important TF involved in lymphangiogenesis.

To understand biological processes and diseases, it is crucial to unravel the concerted interplay of transcription factors (TFs), microRNAs (miRNAs) and their targets within regulatory networks and fundamental sub-networks. An integrative computational resource generating a comprehensive view of these regulatory molecular interactions at a genome-wide scale would be of great interest to biologists, but is not available to date.

The Pax6 transcription factor is expressed in cerebellar granule cells and when mutated, as in the Sey/Sey mouse, produces granule cells with disturbed survival and migration and with defects in neurite extension. The impact of Pax6 on other genes in the context of cerebellar development has not been identified. In this study, we performed transcriptome comparisons between wildtype and Pax6-null whole cerebellar tissue at embryonic day (E) 13.5, 15.5 and 18.5 using Affymetrix arrays (U74Av2). Statistical analyses identified 136 differentially regulated transcripts (FDR 0.05, 1.2-fold change cutoff) over time in Pax6-null cerebellar tissue. In parallel we examined the Math1-null granuloprival cerebellum and identified 228 down-regulated transcripts (FDR 0.05, 1.2-fold change cutoff). The intersection of these two microarray datasets produced a total of 21 differentially regulated transcripts. For a subset of the identified transcripts, we used qRT-PCR to validate the microarray data and demonstrated the expression in the rhombic lip lineage and differential expression in Pax6-null cerebellum with in situ hybridisation analysis. The candidate genes identified in this way represent direct or indirect Pax6-downstream genes involved in cerebellar development.

Phosphatidylinositol glycan biosynthesis class A protein (PIGA) is one of the enzymes involved in the biosynthesis of glycosylphosphatidylinositol (GPI) anchor proteins, which function as enzymes, adhesion molecules, complement regulators and co-receptors in signal transduction pathways. Until recently, only somatic PIGA mutations had been reported in patients with paroxysmal nocturnal hemoglobinuria (PNH), while germline mutations had not been observed, and were suspected to result in lethality. However, in just two years, whole exome sequencing (WES) analyses have identified germline PIGA mutations in male patients with XLIDD (X-linked intellectual developmental disorder) with a wide spectrum of clinical presentations.

Although it is generally accepted that cellular differentiation requires changes to transcriptional networks, dynamic regulation of promoters and enhancers at specific sets of genes has not been previously studied en masse. Exploiting the fact that active promoters and enhancers are transcribed, we simultaneously measured their activity in 19 human and 14 mouse time courses covering a wide range of cell types and biological stimuli. Enhancer RNAs, then messenger RNAs encoding transcription factors, dominated the earliest responses. Binding sites for key lineage transcription factors were simultaneously overrepresented in enhancers and promoters active in each cellular system. Our data support a highly generalizable model in which enhancer transcription is the earliest event in successive waves of transcriptional change during cellular differentiation or activation.

We describe a woman who presented with cataracts, optic atrophy, lipodystrophy/lipoatrophy, and peripheral neuropathy. Exome sequencing identified a c.235C > G p.(Leu79Val) variant in the optic atrophy 3 (OPA3) gene that was confirmed to be de novo. This report expands the severity of the phenotypic spectrum of autosomal dominant OPA3 mutations.

The mammalian inner ear (IE) subserves auditory and vestibular sensations via highly specialized cells and proteins. Sensory receptor hair cells (HCs) are necessary for transducing mechanical inputs and stimulating sensory neurons by using a host of known and as yet unknown protein machinery. To understand the protein composition of these unique postmitotic cells, in which irreversible protein degradation or damage can lead to impaired hearing and balance, we analyzed IE samples by tandem mass spectrometry to generate an unbiased, shotgun-proteomics view of protein identities and abundances. By using Pou4f3/eGFP-transgenic mice in which HCs express GFP driven by Pou4f3, we FACS purified a population of HCs to analyze and compare the HC proteome with other IE subproteomes from sensory epithelia and whole IE. We show that the mammalian HC proteome comprises hundreds of uniquely or highly expressed proteins. Our global proteomic analysis of purified HCs extends the existing HC transcriptome, revealing previously undetected gene products and isoform-specific protein expression. Comparison of our proteomic data with mouse and human databases of genetic auditory/vestibular impairments confirms the critical role of the HC proteome for normal IE function, providing a cell-specific pool of candidates for novel, important HC genes. Several proteins identified exclusively in HCs by proteomics and verified by immunohistochemistry map to human genetic deafness loci, potentially representing new deafness genes.

High-grade serous ovarian cancer (HGSC) exhibits extensive malignant clonal diversity with widespread but non-random patterns of disease dissemination. We investigated whether local immune microenvironment factors shape tumor progression properties at the interface of tumor-infiltrating lymphocytes (TILs) and cancer cells. Through multi-region study of 212 samples from 38 patients with whole-genome sequencing, immunohistochemistry, histologic image analysis, gene expression profiling, and T and B cell receptor sequencing, we identified three immunologic subtypes across samples and extensive within-patient diversity. Epithelial CD8+ TILs negatively associated with malignant diversity, reflecting immunological pruning of tumor clones inferred by neoantigen depletion, HLA I loss of heterozygosity, and spatial tracking between T cell and tumor clones. In addition, combinatorial prognostic effects of mutational processes and immune properties were observed, illuminating how specific genomic aberration types associate with immune response and impact survival. We conclude that within-patient spatial immune microenvironment variation shapes intraperitoneal malignant spread, provoking new evolutionary perspectives on HGSC clonal dispersion.

Evolutionarily conserved RFX transcription factors (TFs) regulate their target genes through a DNA sequence motif called the X-box. Thereby they regulate cellular specialization and terminal differentiation. Here, we provide a comprehensive analysis of all the eight human RFX genes (RFX1-8), their spatial and temporal expression profiles, potential upstream regulators and target genes.

MEDLINE®/PubMed® indexes over 20 million biomedical articles, providing curated annotation of its contents using a controlled vocabulary known as Medical Subject Headings (MeSH). The MeSH vocabulary, developed over 50+ years, provides a broad coverage of topics across biomedical research. Distilling the essential biomedical themes for a topic of interest from the relevant literature is important to both understand the importance of related concepts and discover new relationships.

The global effort to annotate the non-coding portion of the human genome relies heavily on chromatin immunoprecipitation data generated with high-throughput DNA sequencing (ChIP-seq). ChIP-seq is generally successful in detailing the segments of the genome bound by the immunoprecipitated transcription factor (TF), however almost all datasets contain genomic regions devoid of the canonical motif for the TF. It remains to be determined if these regions are related to the immunoprecipitated TF or whether, despite the use of controls, there is a portion of peaks that can be attributed to other causes.

JASPAR (http://jaspar.genereg.net) is the largest open-access database of matrix-based nucleotide profiles describing the binding preference of transcription factors from multiple species. The fifth major release greatly expands the heart of JASPAR-the JASPAR CORE subcollection, which contains curated, non-redundant profiles-with 135 new curated profiles (74 in vertebrates, 8 in Drosophila melanogaster, 10 in Caenorhabditis elegans and 43 in Arabidopsis thaliana; a 30% increase in total) and 43 older updated profiles (36 in vertebrates, 3 in D. melanogaster and 4 in A. thaliana; a 9% update in total). The new and updated profiles are mainly derived from published chromatin immunoprecipitation-seq experimental datasets. In addition, the web interface has been enhanced with advanced capabilities in browsing, searching and subsetting. Finally, the new JASPAR release is accompanied by a new BioPython package, a new R tool package and a new R/Bioconductor data package to facilitate access for both manual and automated methods.

Sialic acid storage diseases are neurodegenerative disorders characterized by accumulation of sialic acid in the lysosome. These disorders are caused by mutations in SLC17A5, the gene encoding sialin, a sialic acid transporter located in the lysosomal membrane. The most common form of sialic acid storage disease is the slowly progressive Salla disease, presenting with hypotonia, ataxia, epilepsy, nystagmus and findings of cerebral and cerebellar atrophy. Hypomyelination and corpus callosum hypoplasia are typical as well. We report a 16 year-old boy with an atypically mild clinical phenotype of sialic acid storage disease characterized by psychomotor retardation and a mixture of spasticity and rigidity but no ataxia, and only weak features of hypomyelination and thinning of corpus callosum on MRI of the brain.

In this study, we developed and evaluated a novel text-mining approach, using non-negative tensor factorization (NTF), to simultaneously extract and functionally annotate transcriptional modules consisting of sets of genes, transcription factors (TFs), and terms from MEDLINE abstracts. A sparse 3-mode term × gene × TF tensor was constructed that contained weighted frequencies of 106,895 terms in 26,781 abstracts shared among 7,695 genes and 994 TFs. The tensor was decomposed into sub-tensors using non-negative tensor factorization (NTF) across 16 different approximation ranks. Dominant entries of each of 2,861 sub-tensors were extracted to form term-gene-TF annotated transcriptional modules (ATMs). More than 94% of the ATMs were found to be enriched in at least one KEGG pathway or GO category, suggesting that the ATMs are functionally relevant. One advantage of this method is that it can discover potentially new gene-TF associations from the literature. Using a set of microarray and ChIP-Seq datasets as gold standard, we show that the precision of our method for predicting gene-TF associations is significantly higher than chance. In addition, we demonstrate that the terms in each ATM can be used to suggest new GO classifications to genes and TFs. Taken together, our results indicate that NTF is useful for simultaneous extraction and functional annotation of transcriptional regulatory networks from unstructured text, as well as for literature based discovery. A web tool called Transcriptional Regulatory Modules Extracted from Literature (TREMEL), available at http://binf1.memphis.edu/tremel, was built to enable browsing and searching of ATMs.

MicroRNAs (miRNAs) can group together along the human genome to form stable secondary structures made of several hairpins hosting miRNAs in their stems. The few known examples of such structures are all involved in cancer development. A large scale computational analysis of human chromosomes crossing sequence analysis and deep sequencing data revealed the presence of >400 structural clusters of miRNAs in the human genome. An a posteriori analysis validates predictions as bona fide miRNAs. A functional analysis of structural clusters position along the chromosomes co-localizes them with genes involved in several key cellular processes like immune systems, sensory systems, signal transduction and development. Immune systems diseases, infectious diseases and neurodegenerative diseases are characterized by genes that are especially well organized around structural clusters of miRNAs. Target genes functional analysis strongly supports a regulatory role of most predicted miRNAs and, notably, a strong involvement of predicted miRNAs in the regulation of cancer pathways. This analysis provides new fundamental insights on the genomic organization of miRNAs in human chromosomes.

In the central nervous system (CNS), myelin is produced from spirally-wrapped oligodendrocyte plasma membrane and, as exemplified by the debilitating effects of inherited or acquired myelin abnormalities in diseases such as multiple sclerosis, it plays a critical role in nervous system function. Myelin sheath production coincides with rapid up-regulation of numerous genes. The complexity of their subsequent expression patterns, along with recently recognized heterogeneity within the oligodendrocyte lineage, suggest that the regulatory networks controlling such genes drive multiple context-specific transcriptional programs. Conferring this nuanced level of control likely involves a large repertoire of interacting transcription factors (TFs). Here, we combined novel strategies of computational sequence analyses with in vivo functional analysis to establish a TF network model of coordinate myelin-associated gene transcription. Notably, the network model captures regulatory DNA elements and TFs known to regulate oligodendrocyte myelin gene transcription and/or oligodendrocyte development, thereby validating our approach. Further, it links to numerous TFs with previously unsuspected roles in CNS myelination and suggests collaborative relationships amongst both known and novel TFs, thus providing deeper insight into the myelin gene transcriptional network.

Fragile X syndrome (FXS) is the most common inherited form of intellectual disability in males and the most common genetic cause of autism. Although executive dysfunction is consistently found in humans with FXS, evidence of executive dysfunction in Fmr1 KO mice, a mouse model of FXS, has been inconsistent. One possible explanation for this is that executive dysfunction in Fmr1 KO mice, similar to humans with FXS, is only evident when cognitive demands are high. Using touchscreen operant conditioning chambers, male Fmr1 KO mice and their male wildtype littermates were tested on the acquisition of a pairwise visual discrimination followed by four serial reversals of the response rule. We assessed reversal learning performance under two different conditions. In the first, the correct stimulus was salient and the incorrect stimulus was non-salient. In the second and more challenging condition, the incorrect stimulus was salient and the correct stimulus was non-salient; this increased cognitive load by introducing conflict between sensory-driven (i.e., bottom-up) and task-dependent (i.e., top-down) signals. Fmr1 KOs displayed two distinct impairments relative to wildtype littermates. First, Fmr1 KOs committed significantly more learning-type errors during the second reversal stage, but only under high cognitive load. Second, during the first reversal stage, Fmr1 KOs committed significantly more attempts to collect a reward during the timeout following an incorrect response. These findings indicate that Fmr1 KO mice display executive dysfunction that, in some cases, is only evident under high cognitive load.

We developed Ulysses as a user-oriented system that uses a process called Interolog Analysis for the parallel analysis and display of protein interactions detected in various species. Ulysses was designed to perform such Interolog Analysis by the projection of model organism interaction data onto homologous human proteins, and thus serves as an accelerator for the analysis of uncharacterized human proteins. The relevance of projections was assessed and validated against published reference collections. All source code is freely available, and the Ulysses system can be accessed via a web interface http://www.cisreg.ca/ulysses.

Evolutionarily conserved sequences within or adjoining orthologous genes often serve as critical cis-regulatory regions. Recent studies have identified long, non-coding genomic regions that are perfectly conserved between human and mouse, termed ultra-conserved regions (UCRs). Here, we focus on UCRs that cluster around genes involved in early vertebrate development; genes conserved over 450 million years of vertebrate evolution.

Unravelling regulatory programs governed by transcription factors (TFs) is fundamental to understanding biological systems. TFCat is a catalog of mouse and human TFs based on a reliable core collection of annotations obtained by expert review of the scientific literature. The collection, including proven and homology-based candidate TFs, is annotated within a function-based taxonomy and DNA-binding proteins are organized within a classification system. All data and user-feedback mechanisms are available at the TFCat portal (http://www.tfcat.ca).