The cellular landscape of the endochondral bone during the transition to extrauterine life.

The cellular complexity of the endochondral bone underlies its essential and pleiotropic roles during organismal life. While the adult bone has received significant attention, we still lack a deep understanding of the perinatal bone cellulome. Here, we have profiled the full composition of the murine endochondral bone at the single-cell level during the transition from fetal to newborn life and in comparison with the adult tissue, with particular emphasis on the mesenchymal compartment. The perinatal bone contains different fibroblastic clusters with blastema-like characteristics in organizing and supporting skeletogenesis, angiogenesis and hematopoiesis. Our data also suggest dynamic inter- and intra-compartment interactions, as well as a bone marrow milieu that seems prone to anti-inflammation, which we hypothesize is necessary to ensure the proper program of lymphopoiesis and the establishment of central and peripheral tolerance in early life. Our study provides an integrative roadmap for the future design of genetic and cellular functional assays to validate cellular interactions and lineage relationships within the perinatal bone.

The little skate genome and the evolutionary emergence of wing-like fins.

Skates are cartilaginous fish whose body plan features enlarged wing-like pectoral fins, enabling them to thrive in benthic environments1,2. However, the molecular underpinnings of this unique trait remain unclear. Here we investigate the origin of this phenotypic innovation by developing the little skate Leucoraja erinacea as a genomically enabled model. Analysis of a high-quality chromosome-scale genome sequence for the little skate shows that it preserves many ancestral jawed vertebrate features compared with other sequenced genomes, including numerous ancient microchromosomes. Combining genome comparisons with extensive regulatory datasets in developing fins-including gene expression, chromatin occupancy and three-dimensional conformation-we find skate-specific genomic rearrangements that alter the three-dimensional regulatory landscape of genes that are involved in the planar cell polarity pathway. Functional inhibition of planar cell polarity signalling resulted in a reduction in anterior fin size, confirming that this pathway is a major contributor to batoid fin morphology. We also identified a fin-specific enhancer that interacts with several hoxa genes, consistent with the redeployment of hox gene expression in anterior pectoral fins, and confirmed its potential to activate transcription in the anterior fin using zebrafish reporter assays. Our findings underscore the central role of genome reorganization and regulatory variation in the evolution of phenotypes, shedding light on the molecular origin of an enigmatic trait.

Expanded FLP toolbox for spatiotemporal protein degradation and transcriptomic profiling in Caenorhabditis elegans.

Control of gene expression in specific tissues and/or at certain stages of development allows the study and manipulation of gene function with high precision. Site-specific genome recombination by the flippase (FLP) and cyclization recombination (Cre) enzymes has proved particularly relevant. Joint efforts of many research groups have led to the creation of efficient FLP and Cre drivers to regulate gene expression in a variety of tissues in Caenorhabditis elegans. Here, we extend this toolkit by the addition of FLP lines that drive recombination specifically in distal tip cells, the somatic gonad, coelomocytes, and the epithelial P lineage. In some cases, recombination-mediated gene knockouts do not completely deplete protein levels due to persistence of long-lived proteins. To overcome this, we developed a spatiotemporally regulated degradation system for green fluorescent fusion proteins based on FLP-mediated recombination. Using 2 stable nuclear pore proteins, MEL-28/ELYS and NPP-2/NUP85 as examples, we report the benefit of combining tissue-specific gene knockout and protein degradation to achieve complete protein depletion. We also demonstrate that FLP-mediated recombination can be utilized to identify transcriptomes in a C. elegans tissue of interest. We have adapted RNA polymerase DamID for the FLP toolbox and by focusing on a well-characterized tissue, the hypodermis, we show that the vast majority of genes identified by RNA polymerase DamID are known to be expressed in this tissue. These tools allow combining FLP activity for simultaneous gene inactivation and transcriptomic profiling, thus enabling the inquiry of gene function in various complex biological processes.

Gene Regulatory Networks of Epidermal and Neural Fate Choice in a Chordate.

Neurons are a highly specialized cell type only found in metazoans. They can be scattered throughout the body or grouped together, forming ganglia or nerve cords. During embryogenesis, centralized nervous systems develop from the ectoderm, which also forms the epidermis. How pluripotent ectodermal cells are directed toward neural or epidermal fates, and to which extent this process is shared among different animal lineages, are still open questions. Here, by using micromere explants, we were able to define in silico the putative gene regulatory networks (GRNs) underlying the first steps of the epidermis and the central nervous system formation in the cephalochordate amphioxus. We propose that although the signal triggering neural induction in amphioxus (i.e., Nodal) is different from vertebrates, the main transcription factors implicated in this process are conserved. Moreover, our data reveal that transcription factors of the neural program seem to not only activate neural genes but also to potentially have direct inputs into the epidermal GRN, suggesting that the Nodal signal might also contribute to neural fate commitment by repressing the epidermal program. Our functional data on whole embryos support this result and highlight the complex interactions among the transcription factors activated by the signaling pathways that drive ectodermal cell fate choice in chordates.

Human prefoldin modulates co-transcriptional pre-mRNA splicing.

Prefoldin is a heterohexameric complex conserved from archaea to humans that plays a cochaperone role during the co-translational folding of actin and tubulin monomers. Additional functions of prefoldin have been described, including a positive contribution to transcription elongation and chromatin dynamics in yeast. Here we show that prefoldin perturbations provoked transcriptional alterations across the human genome. Severe pre-mRNA splicing defects were also detected, particularly after serum stimulation. We found impairment of co-transcriptional splicing during transcription elongation, which explains why the induction of long genes with a high number of introns was affected the most. We detected genome-wide prefoldin binding to transcribed genes and found that it correlated with the negative impact of prefoldin depletion on gene expression. Lack of prefoldin caused global decrease in Ser2 and Ser5 phosphorylation of the RNA polymerase II carboxy-terminal domain. It also reduced the recruitment of the CTD kinase CDK9 to transcribed genes, and the association of splicing factors PRP19 and U2AF65 to chromatin, which is known to depend on CTD phosphorylation. Altogether the reported results indicate that human prefoldin is able to act locally on the genome to modulate gene expression by influencing phosphorylation of elongating RNA polymerase II, and thereby regulating co-transcriptional splicing.

Novel miRNA-mRNA interactions conserved in essential cancer pathways.

Cancer is a complex disease in which unrestrained cell proliferation results in tumour development. Extensive research into the molecular mechanisms underlying tumorigenesis has led to the characterization of oncogenes and tumour suppressors that are key elements in cancer growth and progression, as well as that of other important elements like microRNAs. These genes and miRNAs appear to be constitutively deregulated in cancer. To identify signatures of miRNA-mRNA interactions potentially conserved in essential cancer pathways, we have conducted an integrative analysis of transcriptomic data, also taking into account methylation and copy number alterations. We analysed 18,605 raw transcriptome samples from The Cancer Genome Atlas covering 15 of the most common types of human tumours. From this global transcriptome study, we recovered known cancer-associated miRNA-targets and importantly, we identified new potential targets from miRNA families, also analysing the phenotypic outcomes of these genes/mRNAs in terms of survival. Further analyses could lead to novel approaches in cancer therapy.

DDRprot: a database of DNA damage response-related proteins.

The DNA Damage Response (DDR) signalling network is an essential system that protects the genome's integrity. The DDRprot database presented here is a resource that integrates manually curated information on the human DDR network and its sub-pathways. For each particular DDR protein, we present detailed information about its function. If involved in post-translational modifications (PTMs) with each other, we depict the position of the modified residue/s in the three-dimensional structures, when resolved structures are available for the proteins. All this information is linked to the original publication from where it was obtained. Phylogenetic information is also shown, including time of emergence and conservation across 47 selected species, family trees and sequence alignments of homologues. The DDRprot database can be queried by different criteria: pathways, species, evolutionary age or involvement in (PTM). Sequence searches using hidden Markov models can be also used.Database URL: http://ddr.cbbio.es.

Gene expression profiling in hearts of diabetic mice uncovers a potential role of estrogen-related receptor γ in diabetic cardiomyopathy.

Diabetic cardiomyopathy is characterized by an abnormal oxidative metabolism, but the underlying mechanisms remain to be defined. To uncover potential mechanisms involved in the pathophysiology of diabetic cardiomyopathy, we performed a gene expression profiling study in hearts of diabetic db/db mice. Diabetic hearts showed a gene expression pattern characterized by the up-regulation of genes involved in lipid oxidation, together with an abnormal expression of genes related to the cardiac contractile function. A screening for potential regulators of the genes differentially expressed in diabetic mice found that estrogen-related receptor γ (ERRγ) was increased in heart of db/db mice. Overexpression of ERRγ in cultured cardiomyocytes was sufficient to promote the expression of genes involved in lipid oxidation, increase palmitate oxidation and induce cardiomyocyte hypertrophy. Our findings strongly support a role for ERRγ in the metabolic alterations that underlie the development of diabetic cardiomyopathy.

Gene Expression Profiling Supports the Neural Crest Origin of Adult Rodent Carotid Body Stem Cells and Identifies CD10 as a Marker for Mesectoderm-Committed Progenitors.

Neural stem cells (NSCs) are promising tools for understanding nervous system plasticity and repair, but their use is hampered by the lack of markers suitable for their prospective isolation and characterization. The carotid body (CB) contains a population of peripheral NSCs, which support organ growth during acclimatization to hypoxia. We have set up CB neurosphere (NS) cultures enriched in differentiated neuronal (glomus) cells versus undifferentiated progenitors to investigate molecular hallmarks of cell classes within the CB stem cell (CBSC) niche. Microarray gene expression analysis in NS is compatible with CBSCs being neural crest derived-multipotent progenitor cells able to sustain CB growth upon exposure to hypoxia. Moreover, we have identified CD10 as a marker suitable for isolation of a population of CB mesectoderm-committed progenitor cells. CD10 + cells are resting in normoxia, and during hypoxia they are activated to proliferate and to eventually complete maturation into mesectodermal cells, thus participating in the angiogenesis necessary for CB growth. Our results shed light into the molecular and cellular mechanisms involved in CBSC fate choice, favoring a potential use of these cells for cell therapy. Stem Cells 2016;34:1637-1650.

Notch-regulated miR-223 targets the aryl hydrocarbon receptor pathway and increases cytokine production in macrophages from rheumatoid arthritis patients.

Evidence links aryl hydrocarbon receptor (AHR) activation to rheumatoid arthritis (RA) pathogenesis, although results are inconsistent. AHR agonists inhibit pro-inflammatory cytokine expression in macrophages, pivotal cells in RA aetiopathogenesis, which hints at specific circuits that regulate the AHR pathway in RA macrophages. We compared microRNA (miR) expression in CD14(+) cells from patients with active RA or with osteoarthritis (OA). Seven miR were downregulated and one (miR-223) upregulated in RA compared to OA cells. miR-223 upregulation correlated with reduced Notch3 and Notch effector expression in RA patients. Overexpression of the Notch-induced repressor HEY-1 and co-culture of healthy donor monocytes with Notch ligand-expressing cells showed direct Notch-mediated downregulation of miR-223. Bioinformatics predicted the AHR regulator ARNT (AHR nuclear translocator) as a miR-223 target. Pre-miR-223 overexpression silenced ARNT 3'UTR-driven reporter expression, reduced ARNT (but not AHR) protein levels and prevented AHR/ARNT-mediated inhibition of pro-inflammatory cytokine expression. miR-223 counteracted AHR/ARNT-induced Notch3 upregulation in monocytes. Levels of ARNT and of CYP1B1, an AHR/ARNT signalling effector, were reduced in RA compared to OA synovial tissue, which correlated with miR-223 levels. Our results associate Notch signalling to miR-223 downregulation in RA macrophages, and identify miR-223 as a negative regulator of the AHR/ARNT pathway through ARNT targeting.

Deconstruction of DNA methylation patterns during myogenesis reveals specific epigenetic events in the establishment of the skeletal muscle lineage.

The progressive restriction of differentiation potential from pluripotent embryonic stem cells (ESCs) to tissue-specific stem cells involves widespread epigenetic reprogramming, including modulation of DNA methylation patterns. Skeletal muscle stem cells are required for the growth, maintenance, and regeneration of skeletal muscle. To investigate the contribution of DNA methylation to the establishment of the myogenic program, we analyzed ESCs, skeletal muscle stem cells in proliferating (myoblasts) and differentiating conditions (myotubes), and mature myofibers. About 1.000 differentially methylated regions were identified during muscle-lineage determination and terminal differentiation, mainly located in gene bodies and intergenic regions. As a whole, myogenic stem cells showed a gain of DNA methylation, while muscle differentiation was accompanied by loss of DNA methylation in CpG-poor regions. Notably, the hypomethylated regions in myogenic stem cells were neighbored by enhancer-type chromatin, suggesting the involvement of DNA methylation in the regulation of cell-type specific enhancers. Interestingly, we demonstrated the hypomethylation of the muscle cell-identity Myf5 super-enhancer only in muscle cells. Furthermore, we observed that upstream stimulatory factor 1 binding to Myf5 super-enhancer occurs upon DNA demethylation in myogenic stem cells. Taken altogether, we characterized the unique DNA methylation signature of skeletal muscle stem cells and highlighted the importance of DNA methylation-mediated regulation of cell identity Myf5 super-enhancer during cellular differentiation.

Emergence and evolutionary analysis of the human DDR network: implications in comparative genomics and downstream analyses.

The DNA damage response (DDR) is a crucial signaling network that preserves the integrity of the genome. This network is an ensemble of distinct but often overlapping subnetworks, where different components fulfill distinct functions in precise spatial and temporal scenarios. To understand how these elements have been assembled together in humans, we performed comparative genomic analyses in 47 selected species to trace back their emergence using systematic phylogenetic analyses and estimated gene ages. The emergence of the contribution of posttranslational modifications to the complex regulation of DDR was also investigated. This is the first time a systematic analysis has focused on the evolution of DDR subnetworks as a whole. Our results indicate that a DDR core, mostly constructed around metabolic activities, appeared soon after the emergence of eukaryotes, and that additional regulatory capacities appeared later through complex evolutionary process. Potential key posttranslational modifications were also in place then, with interacting pairs preferentially appearing at the same evolutionary time, although modifications often led to the subsequent acquisition of new targets afterwards. We also found extensive gene loss in essential modules of the regulatory network in fungi, plants, and arthropods, important for their validation as model organisms for DDR studies.

Serine/threonine kinases and E2-ubiquitin conjugating enzymes in Planctomycetes: unexpected findings.

The regulation of signal transduction by phosphorylation and ubiquitination is essential to ensure the correct behavior of eukaryotic cells. We searched for protein families involved in such signaling in several eukaryotic species and in a limited set of prokaryotes, where two members of the Planctomycetes phylum were included as they exhibit eukaryote-like features (Gemmata obscuriglobus and Pirellula staleyi). We identified sequences homologous to eukaryotic serine/threonine kinases (STKs) and E2-ubiquitin conjugating enzymes in the two Planctomycetes species. To extend these analyses to the Planctomycetes/Verrucomicrobia/Chlamydia super-phylum, we performed comparative analyses using domains from kinases, phosphatases and GTPases that serve as signaling signatures, and we analyzed their distributions. We found substantial differences in kinome densities with regards to other prokaryote clades and among the groups in the Planctomycetes/Verrucomicrobia/Chlamydia super-phylum. In addition, we identified the presence of classic eukaryotic E2-conjugating ubiquitin proteins in prokaryotes, these having previously believed to exist only in eukaryotes. Our phylogenetic analyses of the STKs signature domains and E2-enzymes suggest the existence of horizontal gene transfer.

Functional analysis of the integration host factor site of the σ(54) Pu promoter of Pseudomonas putida by in vivo UV imprinting.

The integration host factor (IHF) of Pseudomonas putida connects cell growth to transcriptional activity of distinct promoters. The IHF site of the σ(54) promoter Pu of the TOL (m-xylene biodegradation) plasmid pWW0 of P. putida has been examined to define experimentally a relationship between occupation of the promoter by this factor, the biological activity of the protein and the tolerance of the target site to single-base changes through the bound DNA core sequence. The use of an in vivo high-intensity UV imprinting procedure to examine such an occupation of Pu by IHF allowed inspection of the interplay between the factor and cognate site variants under the physiologically relevant conditions of monocopy gene dosage. The resulting data were merged in a structural model for establishing key features of the IHF-DNA interaction. A functional consensus for first-order IHF binding was instrumental for a genome-wide survey of sequences with potential regulatory value. This search revealed that very few, if any, of the maximum 330 sites within intergenic regions were placed in locations controlling expression of central metabolic genes. It thus seems that the IHF regulon of P. putida has a degree of functional specialization that is not evenly distributed through all gene categories.

SIREs: searching for iron-responsive elements.

The iron regulatory protein/iron-responsive element regulatory system plays a crucial role in the post-transcriptional regulation of gene expression and its disruption results in human disease. IREs are cis-acting regulatory motifs present in mRNAs that encode proteins involved in iron metabolism. They function as binding sites for two related trans-acting factors, namely the IRP-1 and -2. Among cis-acting RNA regulatory elements, the IRE is one of the best characterized. It is defined by a combination of RNA sequence and structure. However, currently available programs to predict IREs do not show a satisfactory level of sensitivity and fail to detect some of the functional IREs. Here, we report an improved software for the prediction of IREs implemented as a user-friendly web server tool. The SIREs web server uses a simple data input interface and provides structure analysis, predicted RNA folds, folding energy data and an overall quality flag based on properties of well characterized IREs. Results are reported in a tabular format and as a schematic visual representation that highlights important features of the IRE. The SIREs (Search for iron-responsive elements) web server is freely available on the web at http://ccbg.imppc.org/sires/index.html.

Systemic approaches to biodegradation.

Biodegradation, the ability of microorganisms to remove complex chemicals from the environment, is a multifaceted process in which many biotic and abiotic factors are implicated. The recent accumulation of knowledge about the biochemistry and genetics of the biodegradation process, and its categorization and formalization in structured databases, has recently opened the door to systems biology approaches, where the interactions of the involved parts are the main subject of study, and the system is analysed as a whole. The global analysis of the biodegradation metabolic network is beginning to produce knowledge about its structure, behaviour and evolution, such as its free-scale structure or its intrinsic robustness. Moreover, these approaches are also developing into useful tools such as predictors for compounds' degradability or the assisted design of artificial pathways. However, it is the environmental application of high-throughput technologies from the genomics, metagenomics, proteomics and metabolomics that harbours the most promising opportunities to understand the biodegradation process, and at the same time poses tremendous challenges from the data management and data mining point of view.

Bionemo: molecular information on biodegradation metabolism.

Bionemo (http://bionemo.bioinfo.cnio.es) stores manually curated information about proteins and genes directly implicated in the Biodegradation metabolism. When possible, the database includes information on sequence, domains and structures for proteins; and sequence, regulatory elements and transcription units for genes. Thus, Bionemo is a unique resource that complements other biodegradation databases such as the University of Minessota Biocatalysis/Biodegradation Database, or Metarouter, which focus more on the biochemical aspects of biodegradation than in the nature of the biomolecules carrying out the reactions. Bionemo has been built by manually associating sequences database entries to biodegradation reactions, using the information extracted from published articles. Information on transcription units and their regulation was also extracted from the literature for biodegradation genes, and linked to the underlying biochemical network. In its current version, Bionemo contains sequence information for 324 reactions and transcription regulation information for more than 100 promoters and 100 transcription factors. The information in the Bionemo database is available via a web server and the full database is also downloadable as a PostgresSQL dump. To facilitate the programmatic use of the information contained in the database, an object-oriented Perl API is also provided.

The ancestral role of the phosphoenolpyruvate-carbohydrate phosphotransferase system (PTS) as exposed by comparative genomics.

The normal role of the phosphoenolpyruvate-carbohydrate phosphotransferase system (PTS) is phosphorylation and subsequent uptake of specific sugars. However, analysis of the distribution of PTS proteins in 206 genomes covering major bacterial groups indicates that the conventional function of PTS proteins as devices for carbohydrate phosphorylation and transport is an exception found in Enterobacteriacea, Vibrionales and Firmicutes, rather than a rule for all bacteria. Instead, available evidence suggests that a core set of C-responsive phosphotransferases have been evolutionarily drafted towards diversity of regulatory functions in response inter alia to the global economy of the C and N pools.

CARGO: a web portal to integrate customized biological information.

There is a huge quantity of information generated in Life Sciences, and it is dispersed in many databases and repositories. Despite the broad availability of the information, there is a great demand for methods that are able to look for, gather and display distributed data in a standardized and friendly way. CARGO (Cancer And Related Genes Online) is a configurable biological web portal designed as a tool to facilitate, integrate and visualize results from Internet resources, independently of their native format or access method. Through the use of small agents, called widgets, supported by a Rich Internet Application (RIA) paradigm based on AJAX, CARGO provides pieces of minimal, relevant and descriptive biological information. The tool is designed to be used by experimental biologists with no training in bioinformatics. In the current state, the system presents a list of human cancer genes. Available at http://cargo.bioinfo.cnio.es.

The phosphotransferase system formed by PtsP, PtsO, and PtsN proteins controls production of polyhydroxyalkanoates in Pseudomonas putida.

The genome of Pseudomonas putida KT2440 encodes five proteins of the phosphoenolpyruvate-carbohydrate phosphotransferase system. Two of these (FruA and FruB) form a dedicated system for fructose intake, while enzyme I(Ntr) (EI(Ntr); encoded by ptsP), NPr (ptsO), and EII(Ntr) (ptsN) act in concert to control the intracellular accumulation of polyhydroxyalkanoates, a typical product of carbon overflow.