Variability in Cerebrospinal Fluid MicroRNAs Through Life.

The development of the human brain starts in the first weeks of embryo differentiation. However, there are many relevant neurodevelopmental processes that take place after birth and during lifespan. Such a fine and changing scenario requires the coordinated expression of thousands of genes to achieve the proper specialization and inter-connectivity. In this context, microRNAs (miRNAs), which can modulate mRNA stability and translation, are gaining recognition for their involvement in both brain development and neurodevelopmental disorders. Therefore, cerebrospinal fluid (CSF) miRNAs should be perfectly differentiated in relevant age periods. In this study, we aimed to highlight the biological variability of miRNA expression in the CSF throughout life, which is also crucial for biomarker discovery in CNS pathologies, especially in children, where they are desperately needed. We analyzed the CSF microRNAome of 14 healthy children (aged 0-7.4 years) by smallRNA-Seq and compared it with previously published data in adults (N = 7) and elders (N = 11). miR-423-5p and miR-22-3p were overexpressed in the < 1 and > 3 years groups, respectively. Additionally, we detected 18 miRNAs that reached their highest peak of expression at different time-points during the lifespan and sets of miRNAs that were exclusively expressed in a specific age group. On the contrary, miR-191-5p showed stable expression in CSF from the first year of life. Our results remark the complex differential miRNA expression profile that can be observed through life, which underlines the need for including appropriate age-matched controls when the expression of CSF miRNAs is analyzed in different pathological contexts. Graphical abstract.

Phosphoinositide 3-Kinase-Regulated Pericyte Maturation Governs Vascular Remodeling.

BACKGROUND: Pericytes regulate vessel stabilization and function, and their loss is associated with diseases such as diabetic retinopathy or cancer. Despite their physiological importance, pericyte function and molecular regulation during angiogenesis remain poorly understood. METHODS: To decipher the transcriptomic programs of pericytes during angiogenesis, we crossed Pdgfrb(BAC)-CreERT2 mice into RiboTagflox/flox mice. Pericyte morphological changes were assessed in mural cell-specific R26-mTmG reporter mice, in which low doses of tamoxifen allowed labeling of single-cell pericytes at high resolution. To study the role of phosphoinositide 3-kinase (PI3K) signaling in pericyte biology during angiogenesis, we used genetic mouse models that allow selective inactivation of PI3Kα and PI3Kß isoforms and their negative regulator phosphate and tensin homolog deleted on chromosome 10 (PTEN) in mural cells. RESULTS: At the onset of angiogenesis, pericytes exhibit molecular traits of cell proliferation and activated PI3K signaling, whereas during vascular remodeling, pericytes upregulate genes involved in mature pericyte cell function, together with a remarkable decrease in PI3K signaling. Immature pericytes showed stellate shape and high proliferation, and mature pericytes were quiescent and elongated. Unexpectedly, we demonstrate that PI3Kß, but not PI3Kα, regulates pericyte proliferation and maturation during vessel formation. Genetic PI3Kß inactivation in pericytes triggered early pericyte maturation. Conversely, unleashing PI3K signaling by means of PTEN deletion delayed pericyte maturation. Pericyte maturation was necessary to undergo vessel remodeling during angiogenesis. CONCLUSIONS: Our results identify new molecular and morphological traits associated with pericyte maturation and uncover PI3Kß activity as a checkpoint to ensure appropriate vessel formation. In turn, our results may open new therapeutic opportunities to regulate angiogenesis in pathological processes through the manipulation of pericyte PI3Kß activity.

Intestinal epithelial deletion of the glucocorticoid receptor NR3C1 alters expression of inflammatory mediators and barrier function.

Glucocorticoids (GCs) are important hormones involved in the regulation of multiple physiologic functions. GCs are also widely used in anti-inflammatory/immunosuppressant drugs. GCs are synthesized by the adrenal cortex as part of the hypothalamus-pituitary-adrenal axis and also by intestinal epithelial cells, among other peripheral sites. GCs are one of the main therapy choices for the exacerbations of inflammatory bowel disease, but they are not useful to prolong remission, and development of tolerance with secondary treatment failure is frequent. Thus, GC actions at the intestinal epithelial level are of great importance, both physiologically and pharmacologically. We generated a tamoxifen-inducible nuclear receptor subfamily 3 group C member 1 (NR3C1)ΔIEC mouse model to study the effects of GCs on epithelial cells in vivo. Nr3c1 deletion in epithelial cells of the small intestine and colon was associated with limited colonic inflammation at 1 wk postdeletion, involving augmented epithelial proliferation and mucus production, plus local and systemic immune/inflammatory changes. This phenotype regressed substantially, but not completely, after 2 wk. The mechanism may involve augmented inflammatory signaling by epithelial cells or defective barrier function. We conclude that the epithelial GC receptor plays a significant role in colonic homeostasis in basal conditions, but its deficiency can be compensated in the short term. Future studies are required to assess the impact of Nr3c1 deletion in other conditions such as experimental colitis.-Aranda, C. J., Arredondo-Amador, M., Ocón, B., Lavín, J. L., Aransay, A. M., Martínez-Augustin, O., Sánchez de Medina, F. Intestinal epithelial deletion of the glucocorticoid receptor NR3C1 alters expression of inflammatory mediators and barrier function.

A Comprehensive Study of Vesicular and Non-Vesicular miRNAs from a Volume of Cerebrospinal Fluid Compatible with Clinical Practice.

Cerebrospinal fluid (CSF) microRNAs (miRNAs) have emerged as potential biomarkers for minimally invasive diagnosis of central nervous system malignancies. However, despite significant advances in recent years, this field still suffers from poor data reproducibility. This is especially true in cases of infants, considered a new subject group. Implementing efficient methods to study miRNAs from clinically realistic CSF volumes is necessary for the identification of new biomarkers. Methods: We compared six protocols for characterizing miRNAs, using 200-µL CSF from infants (aged 0-7). Four of the methods employed extracellular vesicle (EV) enrichment step and the other two obtained the miRNAs directly from cleared CSF. The efficiency of each method was assessed using real-time PCR and small RNA sequencing. We also determined the distribution of miRNAs among different CSF shuttles, using size-exclusion chromatography. Results: We identified 281 CSF miRNAs from infants. We demonstrated that the miRNAs could be efficiently detected using only 200 µL of biofluid in case of at least two of the six methods. In the exosomal fraction, we found 12 miRNAs that might be involved in neurodevelopment. Conclusion: The Norgen and Invitrogen protocols appear suitable for the analysis of a large number of miRNAs using small CSF samples.

Cluster Locator, online analysis and visualization of gene clustering.

Summary: Genes sharing functions, expression patterns or quantitative traits are not randomly distributed along eukaryotic genomes. In order to study the distribution of genes that share a given feature, we present Cluster Locator, an online analysis and visualization tool. Cluster Locator determines the number, size and position of all the clusters formed by the protein-coding genes on a list according to a given maximum gap, the percentage of gene clustering of the list and its statistical significance. The output includes a visual representation of the distribution of genes and gene clusters along the reference genome. Availability and implementation: Cluster Locator is freely available at http://clusterlocator.bnd.edu.uy/. Supplementary information: Supplementary data are available at Bioinformatics online.

Erratum: The metabolic co-regulator PGC1α suppresses prostate cancer metastasis.

This corrects the article DOI: 10.1038/ncb3357.

A synbiotic composed of Lactobacillus fermentum CECT5716 and FOS prevents the development of fatty acid liver and glycemic alterations in rats fed a high fructose diet associated with changes in the microbiota.

We investigated the effect of a high fructose diet (HFD) on Sprague Dawley rats and the impact of a synbiotic composed of Lactobacillus fermentum CECT5716 and fructooligosaccharides. Feeding the HFD for 5 weeks resulted in liver steatosis and insulin resistance but not obesity. These changes were associated with increased production of short-chain fatty acids and increased Bacteroidetes in feces, with an augmented Bacteroidetes/Firmicutes ratio, among other changes in the microbiota. In addition, barrier function was weakened, with increased LPS plasma levels. These data are consistent with increased fructose availability in the distal gut due to saturation of absorptive mechanisms, leading to dysbiosis, endotoxemia, hepatic steatosis, and insulin resistance. Treatment with the synbiotic prevented some of the pathological effects, so that treated rats did not develop steatosis or systemic inflammation, while dysbiosis and barrier function were greatly ameliorated. In addition, the synbiotic had hypolipidemic effects. The synbiotic composed by L. fermentum CECT5716 and fructooligosaccharides has beneficial effects in a model of metabolic syndrome induced by a HFD, suggesting it might be clinically useful in this type of condition, particularly considering that high fructose intake has been related to metabolic syndrome in humans.

The metabolic co-regulator PGC1α suppresses prostate cancer metastasis.

Cellular transformation and cancer progression is accompanied by changes in the metabolic landscape. Master co-regulators of metabolism orchestrate the modulation of multiple metabolic pathways through transcriptional programs, and hence constitute a probabilistically parsimonious mechanism for general metabolic rewiring. Here we show that the transcriptional co-activator peroxisome proliferator-activated receptor gamma co-activator 1α (PGC1α) suppresses prostate cancer progression and metastasis. A metabolic co-regulator data mining analysis unveiled that PGC1α is downregulated in prostate cancer and associated with disease progression. Using genetically engineered mouse models and xenografts, we demonstrated that PGC1α opposes prostate cancer progression and metastasis. Mechanistically, the use of integrative metabolomics and transcriptomics revealed that PGC1α activates an oestrogen-related receptor alpha (ERRα)-dependent transcriptional program to elicit a catabolic state and metastasis suppression. Importantly, a signature based on the PGC1α-ERRα pathway exhibited prognostic potential in prostate cancer, thus uncovering the relevance of monitoring and manipulating this pathway for prostate cancer stratification and treatment.

Global gene expression shift during the transition from early neural development to late neuronal differentiation in Drosophila melanogaster.

Regulation of transcription is one of the mechanisms involved in animal development, directing changes in patterning and cell fate specification. Large temporal data series, based on microarrays across the life cycle of the fly Drosophila melanogaster, revealed the existence of groups of genes which expression increases or decreases temporally correlated during the life cycle. These groups of genes are enriched in different biological functions. Here, instead of searching for temporal coincidence in gene expression using the entire genome expression data, we searched for temporal coincidence in gene expression only within predefined catalogues of functionally related genes and investigated whether a catalogue's expression profile can be used to generate larger catalogues, enriched in genes necessary for the same function. We analyzed the expression profiles from genes already associated with early neurodevelopment and late neurodifferentiation, at embryonic stages 16 and 17 of Drosophila life cycle. We hypothesized that during this interval we would find global downregulation of genes important for early neuronal development together with global upregulation of genes necessary for the final differentiation of neurons. Our results were consistent with this hypothesis. We then investigated if the expression profile of gene catalogues representing particular processes of neural development matched the temporal sequence along which these processes occur. The profiles of genes involved in patterning, neurogenesis, axogenesis or synaptic transmission matched the prediction, with largest transcript values at the time when the corresponding biological process takes place in the embryo. Furthermore, we obtained catalogues enriched in genes involved in temporally matching functions by performing a genome-wide systematic search for genes with their highest expression levels at the corresponding embryonic intervals. These findings imply the use of gene expression data in combination with known biological information to predict the involvement of functionally uncharacterized genes in particular biological events.

S-adenosylmethionine levels regulate the schwann cell DNA methylome.

Axonal myelination is essential for rapid saltatory impulse conduction in the nervous system, and malformation or destruction of myelin sheaths leads to motor and sensory disabilities. DNA methylation is an essential epigenetic modification during mammalian development, yet its role in myelination remains obscure. Here, using high-resolution methylome maps, we show that DNA methylation could play a key gene regulatory role in peripheral nerve myelination and that S-adenosylmethionine (SAMe), the principal methyl donor in cytosine methylation, regulates the methylome dynamics during this process. Our studies also point to a possible role of SAMe in establishing the aberrant DNA methylation patterns in a mouse model of diabetic neuropathy, implicating SAMe in the pathogenesis of this disease. These critical observations establish a link between SAMe and DNA methylation status in a defined biological system, providing a mechanism that could direct methylation changes during cellular differentiation and in diverse pathological situations.

Bivariate segmentation of SNP-array data for allele-specific copy number analysis in tumour samples.

BACKGROUND: SNP arrays output two signals that reflect the total genomic copy number (LRR) and the allelic ratio (BAF), which in combination allow the characterisation of allele-specific copy numbers (ASCNs). While methods based on hidden Markov models (HMMs) have been extended from array comparative genomic hybridisation (aCGH) to jointly handle the two signals, only one method based on change-point detection, ASCAT, performs bivariate segmentation. RESULTS: In the present work, we introduce a generic framework for bivariate segmentation of SNP array data for ASCN analysis. For the matter, we discuss the characteristics of the typically applied BAF transformation and how they affect segmentation, introduce concepts of multivariate time series analysis that are of concern in this field and discuss the appropriate formulation of the problem. The framework is implemented in a method named CnaStruct, the bivariate form of the structural change model (SCM), which has been successfully applied to transcriptome mapping and aCGH. CONCLUSIONS: On a comprehensive synthetic dataset, we show that CnaStruct outperforms the segmentation of existing ASCN analysis methods. Furthermore, CnaStruct can be integrated into the workflows of several ASCN analysis tools in order to improve their performance, specially on tumour samples highly contaminated by normal cells.

Whole transcriptome analysis of a reversible neurodegenerative process in Drosophila reveals potential neuroprotective genes.

BACKGROUND: Neurodegenerative diseases are progressive and irreversible and they can be initiated by mutations in specific genes. Spalt-like genes (Sall) encode transcription factors expressed in the central nervous system. In humans, SALL mutations are associated with hereditary syndromes characterized by mental retardation, sensorineural deafness and motoneuron problems, among others. Drosophila sall mutants exhibit severe neurodegeneration of the central nervous system at embryonic stage 16, which surprisingly reverts later in development at embryonic stage 17, suggesting a potential to recover from neurodegeneration. We hypothesize that this recovery is mediated by a reorganization of the transcriptome counteracting SALL lost. To identify genes associated to neurodegeneration and neuroprotection, we used mRNA-Seq to compare the transcriptome of Drosophila sall mutant and wild type embryos from neurodegeneration and reversal stages. RESULTS: Neurodegeneration stage is associated with transcriptional changes in 220 genes, of which only 5% were already described as relevant for neurodegeneration. Genes related to the groups of Redox, Lifespan/Aging and Mitochondrial diseases are significantly represented at this stage. By contrast, neurodegeneration reversal stage is associated with significant changes in 480 genes, including 424 not previously associated with neuroprotection. Immune response and Salt stress are the most represented groups at this stage. CONCLUSIONS: We identify new genes associated to neurodegeneration and neuroprotection by using an mRNA-Seq approach. The strong homology between Drosophila and human genes raises the possibility to unveil novel genes involved in neurodegeneration and neuroprotection also in humans.

Comparison of methods to detect copy number alterations in cancer using simulated and real genotyping data.

BACKGROUND: The detection of genomic copy number alterations (CNA) in cancer based on SNP arrays requires methods that take into account tumour specific factors such as normal cell contamination and tumour heterogeneity. A number of tools have been recently developed but their performance needs yet to be thoroughly assessed. To this aim, a comprehensive model that integrates the factors of normal cell contamination and intra-tumour heterogeneity and that can be translated to synthetic data on which to perform benchmarks is indispensable. RESULTS: We propose such model and implement it in an R package called CnaGen to synthetically generate a wide range of alterations under different normal cell contamination levels. Six recently published methods for CNA and loss of heterozygosity (LOH) detection on tumour samples were assessed on this synthetic data and on a dilution series of a breast cancer cell-line: ASCAT, GAP, GenoCNA, GPHMM, MixHMM and OncoSNP. We report the recall rates in terms of normal cell contamination levels and alteration characteristics: length, copy number and LOH state, as well as the false discovery rate distribution for each copy number under different normal cell contamination levels.Assessed methods are in general better at detecting alterations with low copy number and under a little normal cell contamination levels. All methods except GPHMM, which failed to recognize the alteration pattern in the cell-line samples, provided similar results for the synthetic and cell-line sample sets. MixHMM and GenoCNA are the poorliest performing methods, while GAP generally performed better. This supports the viability of approaches other than the common hidden Markov model (HMM)-based. CONCLUSIONS: We devised and implemented a comprehensive model to generate data that simulate tumoural samples genotyped using SNP arrays. The validity of the model is supported by the similarity of the results obtained with synthetic and real data. Based on these results and on the software implementation of the methods, we recommend GAP for advanced users and GPHMM for a fully driven analysis.

The RNA-binding protein human antigen R controls global changes in gene expression during Schwann cell development.

An important prerequisite to myelination in peripheral nerves is the establishment of one-to-one relationships between axons and Schwann cells. This patterning event depends on immature Schwann cell proliferation, apoptosis, and morphogenesis, which are governed by coordinated changes in gene expression. Here, we found that the RNA-binding protein human antigen R (HuR) was highly expressed in immature Schwann cells, where genome-wide identification of its target mRNAs in vivo in mouse sciatic nerves using ribonomics showed an enrichment of functionally related genes regulating these processes. HuR coordinately regulated expression of several genes to promote proliferation, apoptosis, and morphogenesis in rat Schwann cells, in response to NRG1, TGFß, and laminins, three major signals implicated in this patterning event. Strikingly, HuR also binds to several mRNAs encoding myelination-related proteins but, contrary to its typical function, negatively regulated their expression, likely to prevent ectopic myelination during development. These functions of HuR correlated with its abundance and subcellular localization, which were regulated by different signals in Schwann cells.

Distinct roles for Wnt-4 and Wnt-11 during retinoic acid-induced neuronal differentiation.

Retinoic acid and Wnt/ß-catenin signals play important roles during neuronal differentiation but less is known about noncanonical Wnt signals in this context. We examined retinoic acid and Wnt signaling in two human embryonal carcinoma cell lines, NTERA-2 (clone D1), which undergoes neuronal differentiation in response to retinoic acid, and 2102Ep, which does not. Retinoic acid treatment inhibited ß-catenin/Tcf activity in NTERA-2 cells but not in 2102Ep cells. Inhibition occurred downstream of ß-catenin but did not involve competition between retinoic acid receptors and ß-catenin for binding to p300 or Tcf-4. Ectopic expression of FZD1 partially restored inhibition in 2102Ep cells, suggesting the involvement of Wnt ligands. Retinoic acid treatment of NTERA-2 cells induced the expression of Wnt-4 and Wnt-11, both of which were able to inhibit ß-catenin/Tcf activity. Wnt-4 and Wnt-11 were found at cell borders in islands of cells that expressed OCT4 and GFAP and were predominantly negative for Nestin, PAX6, and GATA6. Gene silencing of Wnt-4, but not Wnt-11, reduced retinoic acid downregulation of OCT4 and Nanog and upregulation of PAX6, ASCL1, HOXC5, and NEUROD1, suggesting that Wnt-4 promotes early neuronal differentiation. Gene expression analysis of NTERA-2 cells stably overexpressing Wnt-11 suggested that Wnt-11 potentiates retinoic acid induction of early neurogenesis. Consistent with this, overexpression of Wnt-11 maintained a population of proliferating progenitor cells in cultures treated with retinoic acid for several weeks. These observations highlight the distinct roles of two noncanonical Wnts during the early stages of retinoic acid-induced neuronal differentiation.

High-density SNP genotyping detects homogeneity of Spanish and French Basques, and confirms their genomic distinctiveness from other European populations.

A recent study reported that Basques do not constitute a genetically distinct population, and that Basques from Spanish and French provinces do not show significant genetic similarity. These conclusions disagree with numerous previous studies, and are not consistent with the historical and linguistic evidence that supports the distinctiveness of Basques. In order to further investigate this controversy, we have genotyped 83 Spanish Basque individuals and used these data to infer population structure based on more than 60,000 single nucleotide polymorphisms of several European populations. Here, we present the first high-throughput analysis including Basques from Spanish and French provinces, and show that all Basques constitute a homogeneous group that can be clearly differentiated from other European populations.

Genetic diversity of Toscana virus.

Distribution of Toscana virus (TOSV) is evolving with climate change, and pathogenicity may be higher in nonexposed populations outside areas of current prevalence (Mediterranean Basin). To characterize genetic diversity of TOSV, we determined the coding sequences of isolates from Spain and France. TOSV is more diverse than other well-studied phleboviruses (e.g.,Rift Valley fever virus).

Combined functional and positional gene information for the identification of susceptibility variants in celiac disease.

BACKGROUND & AIMS: Celiac disease is a complex, immune-mediated disorder of the intestinal mucosa with a strong genetic component. HLA-DQ2 is the major determinant of risk, but other minor genes, still to be identified, also are involved. METHODS: We designed a strategy that combines gene expression profiling of intestinal biopsy specimens, linkage region information, and different bioinformatics tools for the selection of potentially regulatory single-nucleotide polymorphisms (SNPs) involved in the disease. We selected 361 SNPs from 71 genes that fulfilled stringent functional (changes in expression level) and positional criteria (located in regions that have been linked to the disease, other than HLA). These polymorphisms were genotyped in 262 celiac patients and 214 controls. RESULTS: We detected strong evidence of association with several SNPs (the most significant were rs6747096, P = 2.38 x 10(-5); rs7040561, P = 6.55 x 10(-5); and rs458046, P = 1.35 x 10(-4)) that pinpoint novel candidate determinants of predisposition to the disease in previously identified linkage regions (eg, SERPINE2 in 2q33, and PBX3 or PPP6C in 9q34). CONCLUSIONS: Our study shows that the combination of function and position is a valid strategy for the genetic dissection of complex traits.