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.

Functionality of the non-coding transcripts encoded by the human genome is the coveted goal of the modern genomics research. While commonly relied on the classical methods of forward genetics, integration of different genomics datasets in a global Systems Biology fashion presents a more productive avenue of achieving this very complex aim. Here we report application of a Systems Biology-based approach to dissect functionality of a newly identified vast class of very long intergenic non-coding (vlinc) RNAs. Using highly quantitative FANTOM5 CAGE dataset, we show that these RNAs could be grouped into 1542 novel human genes based on analysis of insulators that we show here indeed function as genomic barrier elements. We show that vlinc RNAs genes likely function in cisto activate nearby genes. This effect while most pronounced in closely spaced vlinc RNA-gene pairs can be detected over relatively large genomic distances. Furthermore, we identified 101 vlinc RNA genes likely involved in early embryogenesis based on patterns of their expression and regulation. We also found another 109 such genes potentially involved in cellular functions also happening at early stages of development such as proliferation, migration and apoptosis. Overall, we show that Systems Biology-based methods have great promise for functional annotation of non-coding RNAs.

Mesothelioma is a highly malignant tumor that is primarily caused by occupational or environmental exposure to asbestos fibers. Despite worldwide restrictions on asbestos usage, further cases are expected as diagnosis is typically 20-40 years after exposure. Once diagnosed there is a very poor prognosis with a median survival rate of 9 months. Considering this the development of early pre clinical diagnostic markers may help improve clinical outcomes.

The analysis of CAGE (Cap Analysis of Gene Expression) time-course has been proposed by the FANTOM5 Consortium to extend the understanding of the sequence of events facilitating cell state transition at the level of promoter regulation. To identify the most prominent transcriptional regulations induced by growth factors in human breast cancer, we apply here the Complexity Invariant Dynamic Time Warping motif EnRichment (CIDER) analysis approach to the CAGE time-course datasets of MCF-7 cells stimulated by epidermal growth factor (EGF) or heregulin (HRG). We identify a multi-level cascade of regulations rooted by the Serum Response Factor (SRF) transcription factor, connecting the MAPK-mediated transduction of the HRG stimulus to the negative regulation of the MAPK pathway by the members of the DUSP family phosphatases. The finding confirms the known primary role of FOS and FOSL1, members of AP-1 family, in shaping gene expression in response to HRG induction. Moreover, we identify a new potential regulation of DUSP5 and RARA (known to antagonize the transcriptional regulation induced by the estrogen receptors) by the activity of the AP-1 complex, specific to HRG response. The results indicate that a divergence in AP-1 regulation determines cellular changes of breast cancer cells stimulated by ErbB receptors.

The generation of myeloid cells from their progenitors is regulated at the level of transcription by combinatorial control of key transcription factors influencing cell-fate choice. To unravel the global dynamics of this process at the transcript level, we generated transcription profiles for 91 human cell types of myeloid origin by use of CAGE profiling. The CAGE sequencing of these samples has allowed us to investigate diverse aspects of transcription control during myelopoiesis, such as identification of novel transcription factors, miRNAs, and noncoding RNAs specific to the myeloid lineage. We further reconstructed a transcription regulatory network by clustering coexpressed transcripts and associating them with enriched cis-regulatory motifs. With the use of the bidirectional expression as a proxy for enhancers, we predicted over 2000 novel enhancers, including an enhancer 38 kb downstream of IRF8 and an intronic enhancer in the KIT gene locus. Finally, we highlighted relevance of these data to dissect transcription dynamics during progressive maturation of granulocyte precursors. A multifaceted analysis of the myeloid transcriptome is made available (www.myeloidome.roslin.ed.ac.uk). This high-quality dataset provides a powerful resource to study transcriptional regulation during myelopoiesis and to infer the likely functions of unannotated genes in human innate immunity.

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.

MicroRNAs (miRNAs) modulate the post-transcriptional regulation of target genes and are related to biology of complex human traits, but genetic landscape of miRNAs remains largely unknown. Given the strikingly tissue-specific miRNA expression profiles, we here expand a previous method to quantitatively evaluate enrichment of genome-wide association study (GWAS) signals on miRNA-target gene networks (MIGWAS) to further estimate tissue-specific enrichment. Our approach integrates tissue-specific expression profiles of miRNAs (∼1800 miRNAs in 179 cells) with GWAS to test whether polygenic signals enrich in miRNA-target gene networks and whether they fall within specific tissues. We applied MIGWAS to 49 GWASs (nTotal = 3 520 246), and successfully identified biologically relevant tissues. Further, MIGWAS could point miRNAs as candidate biomarkers of the trait. As an illustrative example, we performed differentially expressed miRNA analysis between rheumatoid arthritis (RA) patients and healthy controls (n = 63). We identified novel biomarker miRNAs (e.g. hsa-miR-762) by integrating differentially expressed miRNAs with MIGWAS results for RA, as well as novel associated loci with significant genetic risk (rs56656810 at MIR762 at 16q11; n = 91 482, P = 3.6 × 10-8). Our result highlighted that miRNA-target gene network contributes to human disease genetics in a cell type-specific manner, which could yield an efficient screening of miRNAs as promising biomarkers.

Activated mouse macrophages metabolize arginine via NO synthase (NOS2) to produce NO as an antimicrobial effector. Published gene expression datasets provide little support for the activation of this pathway in human macrophages. Generation of NO requires the coordinated regulation of multiple genes. We have generated RNA-sequencing data from bone marrow-derived macrophages from representative rodent (rat), monogastric (pig and horse), and ruminant (sheep, goat, cattle, and water buffalo) species, and analyzed the expression of genes involved in arginine metabolism in response to stimulation with LPS. In rats, as in mice, LPS strongly induced Nos2, the arginine transporter Slc7a2, arginase 1 (Arg1), GTP cyclohydrolase (Gch1), and argininosuccinate synthase (Ass1). None of these responses was conserved across species. Only cattle and water buffalo showed substantial NOS2 induction. The species studied also differed in expression and regulation of arginase (ARG2, rather than ARG1), and amino acid transporters. Variation between species was associated with rapid promoter evolution. Differential induction of NOS2 and ARG2 between the ruminant species was associated with insertions of the Bov-A2 retrotransposon in the promoter region. Bov-A2 was shown to possess LPS-inducible enhancer activity in transfected RAW264.7 macrophages. Consistent with a function in innate immunity, NO production and arginine metabolism vary greatly between species and differences may contribute to pathogen host restriction.

Deciphering the most common modes by which chromatin regulates transcription, and how this is related to cellular status and processes is an important task for improving our understanding of human cellular biology. The FANTOM5 and ENCODE projects represent two independent large scale efforts to map regulatory and transcriptional features to the human genome. Here we investigate chromatin features around a comprehensive set of transcription start sites in four cell lines by integrating data from these two projects.

The availability of fast alignment-free algorithms has greatly reduced the computational burden of RNA-seq processing, especially for relatively poorly assembled genomes. Using these approaches, previous RNA-seq datasets could potentially be processed and integrated with newly sequenced libraries. Confounding factors in such integration include sequencing depth and methods of RNA extraction and selection. Different selection methods (typically, either polyA-selection or rRNA-depletion) omit different RNAs, resulting in different fractions of the transcriptome being sequenced. In particular, rRNA-depleted libraries sample a broader fraction of the transcriptome than polyA-selected libraries. This study aimed to develop a systematic means of accounting for library type that allows data from these two methods to be compared.

Eimeria species cause the intestinal disease coccidiosis, most notably in poultry. While the direct impact of coccidiosis on animal health and welfare is clear, its influence on the enteric microbiota and by-stander effects on chicken health and production remains largely unknown, with the possible exception of Clostridium perfringens (necrotic enteritis). This study evaluated the composition and structure of the caecal microbiome in the presence or absence of a defined Eimeria tenella challenge infection in Cobb500 broiler chickens using 16S rRNA amplicon sequencing. The severity of clinical coccidiosis in individual chickens was quantified by caecal lesion scoring and microbial changes associated with different lesion scores identified. Following E. tenella infection the diversity of taxa within the caecal microbiome remained largely stable. However, infection induced significant changes in the abundance of some microbial taxa. The greatest changes were detected in birds displaying severe caecal pathology; taxa belonging to the order Enterobacteriaceae were increased, while taxa from Bacillales and Lactobacillales were decreased with the changes correlated with lesion severity. Significantly different profiles were also detected in infected birds which remained asymptomatic (lesion score 0), with taxa belonging to the genera Bacteroides decreased and Lactobacillus increased. Many differential taxa from the order Clostridiales were identified, with some increasing and others decreasing in abundance in Eimeria-infected animals. The results support the view that caecal microbiome dysbiosis associated with Eimeria infection contributes to disease pathology, and could be a target for intervention to mitigate the impact of coccidiosis on poultry productivity and welfare. This work highlights that E. tenella infection has a significant impact on the abundance of some caecal bacteria with notable differences detected between lesion score categories emphasising the importance of accounting for differences in caecal lesions when investigating the relationship between E. tenella and the poultry intestinal microbiome.

The FANTOM5 expression atlas is a quantitative measurement of the activity of nearly 200,000 promoter regions across nearly 2,000 different human primary cells, tissue types and cell lines. Generation of this atlas was made possible by the use of CAGE, an experimental approach to localise transcription start sites at single-nucleotide resolution by sequencing the 5' ends of capped RNAs after their conversion to cDNAs. While 50% of CAGE-defined promoter regions could be confidently associated to adjacent transcriptional units, nearly 100,000 promoter regions remained gene-orphan. To address this, we used the CAGEscan method, in which random-primed 5'-cDNAs are paired-end sequenced. Pairs starting in the same region are assembled in transcript models called CAGEscan clusters. Here, we present the production and quality control of CAGEscan libraries from 56 FANTOM5 RNA sources, which enhances the FANTOM5 expression atlas by providing experimental evidence associating core promoter regions with their cognate transcripts.

Colony-stimulating factor (CSF)-1 controls the survival, proliferation, and differentiation of macrophages, which are recognized as scavengers and agents of the innate and the acquired immune systems. Because of their plasticity, macrophages are endowed with many other essential roles during development and tissue homeostasis. We present evidence that CSF-1 plays an important trophic role in postnatal organ growth and kidney repair. Notably, the injection of CSF-1 postnatally enhanced kidney weight and volume and was associated with increased numbers of tissue macrophages. Moreover, CSF-1 promotes postnatal renal repair in mice after ischemia-reperfusion injury by recruiting and influencing macrophages toward a reparative state. CSF-1 treatment rapidly accelerated renal repair with tubular epithelial cell replacement, attenuation of interstitial fibrosis, and functional recovery. Analysis of macrophages from CSF-1-treated kidneys showed increased expression of insulin-like growth factor-1 and anti-inflammatory genes that are known CSF-1 targets. Taken together, these data suggest that CSF-1 is important in kidney growth and the promotion of endogenous repair and resolution of inflammatory injury.

The ten mouse and six human members of the Schlafen (Slfn) gene family all contain an AAA domain. Little is known of their function, but previous studies suggest roles in immune cell development. In this report, we assessed Slfn regulation and function in macrophages, which are key cellular regulators of innate immunity.

The number of mammalian transcripts identified by full-length cDNA projects and genome sequencing projects is increasing remarkably. Clustering them into a strictly nonredundant and comprehensive set provides a platform for functional analysis of the transcriptome and proteome, but the quality of the clustering and predictive usefulness have previously required manual curation to identify truncated transcripts and inappropriate clustering of closely related sequences. A Representative Transcript and Protein Sets (RTPS) pipeline was previously designed to identify the nonredundant and comprehensive set of mouse transcripts based on clustering of a large mouse full-length cDNA set (FANTOM2). Here we propose an alternative method that is more robust, requires less manual curation, and is applicable to other organisms in addition to mouse. RTPSs of human, mouse, and rat have been produced by this method and used for validation. Their comprehensiveness and quality are discussed by comparison with other clustering approaches. The RTPSs are available at .

The beta-arrestins (ARRB1 and ARRB2) regulate G-protein coupled receptor (GPCR) dependent- and independent-signaling pathways and are ubiquitously expressed. Here we show that ARRB2 mRNA and protein expression is enriched in macrophages, and that it regulates complement C1q expression and cell survival. Basal and Toll-like receptor (TLR) inducible expression of mRNAs encoding the complement subcomponents C1qa, C1qb and C1qc was greatly reduced in bone marrow-derived macrophages (BMM) from ARRB2-deficient, but not ARRB1-deficient mice, while factor-independent survival of ARRB2(-/-) BMM was enhanced compared to wildtype BMM. TatARRB2(23), a cell-permeable peptide that contains the MAPK JNK-binding motif from within the ARRB2 C-domain, impaired ARRB2 interaction with JNK3, down-regulated C1q expression and permitted factor-independent survival in BMM, thus suggesting that this peptide antagonises ARRB2 function in macrophages. In addition, TatARRB2(23) transiently activated the phosphorylation of JNK and ERK, but not p38 in BMM. These data imply that ARRB2 acts to limit JNK/ERK activation and survival in macrophages, but is required for basal and TLR-inducible complement C1q expression. Given that loss of C1q function is strongly associated with the development of systemic lupus erythematosus, ARRB2 may act to limit the development of autoimmune disease.

Gene regulatory networks in living cells are controlled by the interaction of multiple cell type-specific transcription regulators with DNA binding sites in target genes. Interferon regulatory factor 8 (IRF8), also known as interferon consensus sequence binding protein (ICSBP), is a transcription factor expressed predominantly in myeloid and lymphoid cell lineages. To find the functional direct target genes of IRF8, the gene expression profiles of siRNA knockdown samples and genome-wide binding locations by ChIP-chip were analyzed in THP-1 myelomonocytic leukemia cells. Consequently, 84 genes were identified as functional direct targets. The ETS family transcription factor PU.1, also known as SPI1, binds to IRF8 and regulates basal transcription in macrophages. Using the same approach, we identified 53 direct target genes of PU.1; these overlapped with 19 IRF8 targets. These 19 genes included key molecules of IFN signaling such as OAS1 and IRF9, but excluded other IFN-related genes amongst the IRF8 functional direct target genes. We suggest that IRF8 and PU.1 can have both combined, and independent actions on different promoters in myeloid cells.

Genetic variants underlying complex traits, including disease susceptibility, are enriched within the transcriptional regulatory elements, promoters and enhancers. There is emerging evidence that regulatory elements associated with particular traits or diseases share similar patterns of transcriptional activity. Accordingly, shared transcriptional activity (coexpression) may help prioritise loci associated with a given trait, and help to identify underlying biological processes. Using cap analysis of gene expression (CAGE) profiles of promoter- and enhancer-derived RNAs across 1824 human samples, we have analysed coexpression of RNAs originating from trait-associated regulatory regions using a novel quantitative method (network density analysis; NDA). For most traits studied, phenotype-associated variants in regulatory regions were linked to tightly-coexpressed networks that are likely to share important functional characteristics. Coexpression provides a new signal, independent of phenotype association, to enable fine mapping of causative variants. The NDA coexpression approach identifies new genetic variants associated with specific traits, including an association between the regulation of the OCT1 cation transporter and genetic variants underlying circulating cholesterol levels. NDA strongly implicates particular cell types and tissues in disease pathogenesis. For example, distinct groupings of disease-associated regulatory regions implicate two distinct biological processes in the pathogenesis of ulcerative colitis; a further two separate processes are implicated in Crohn's disease. Thus, our functional analysis of genetic predisposition to disease defines new distinct disease endotypes. We predict that patients with a preponderance of susceptibility variants in each group are likely to respond differently to pharmacological therapy. Together, these findings enable a deeper biological understanding of the causal basis of complex traits.

The majority of chickens in sub-Saharan Africa are indigenous ecotypes, well adapted to the local environment and raised in scavenging production systems. Although they are generally resilient to disease challenge, routine vaccination and biosecurity measures are rarely applied and infectious diseases remain a major cause of mortality and reduced productivity. Management and genetic improvement programmes are hampered by lack of routine data recording. Selective breeding based on genomic technologies may provide the means to enhance sustainability. In this study, we investigated the genetic architecture of antibody response to four major infectious diseases [infectious bursal disease (IBDV), Marek's disease (MDV), fowl typhoid (SG), fowl cholera (PM)] and resistance to Eimeria and cestode parasitism, along with two production traits [body weight and body condition score (BCS)] in two distinct indigenous Ethiopian chicken ecotypes. We conducted variance component analyses, genome-wide association studies, and pathway and selective sweep analyses.

Lysyl hydroxylation and arginyl demethylation are post-translational events that are important for many cellular processes. The jumonji domain containing protein 6 (JMJD6) has been reported to catalyze both lysyl hydroxylation and arginyl demethylation on diverse protein substrates. It also interacts directly with RNA. This review summarizes knowledge of JMJD6 functions that have emerged in the last 15 years and considers how a single Jumonji C (JmjC) domain-containing enzyme can target so many different substrates. New links and synergies between the three main proposed functions of Jmjd6 in histone demethylation, promoter proximal pause release of polymerase II and RNA splicing are discussed. The physiological context of the described molecular functions is considered and recently described novel roles for JMJD6 in cancer and immune biology are reviewed. The increased knowledge of JMJD6 functions has wider implications for our general understanding of the JmjC protein family of which JMJD6 is a member.