miRNA profiling in leaf and cork tissues of Quercus suber reveals novel miRNAs and tissue-specific expression patterns

The differentiation of cork (phellem) cells from the phellogen (cork cambium) is a secondary growth process observed in the cork oak tree conferring a unique ability to produce a thick layer of cork. At present, the molecular regulators of phellem differentiation are unknown. The previously documented involvement of microRNAs (miRNAs) in the regulation of developmental processes, including secondary growth, motivated the search for these regulators in cork oak tissues. We performed deep sequencing of the small RNA fraction obtained from cork oak leaves and differentiating phellem. RNA sequences with lengths of 19–25 nt derived from the two libraries were analysed, leading to the identification of 41 families of conserved miRNAs, of which the most abundant were miR167, miR165/166, miR396 and miR159. Thirty novel miRNA candidates were also unveiled, 11 of which were unique to leaves and 13 to phellem. Northern blot detection of a set of conserved and novel miRNAs confirmed their differential expression profile. Prediction and analysis of putative miRNA target genes provided clues regarding processes taking place in leaf and phellem tissues, but further experimental work will be needed for functional characterization. In conclusion, we here provide a first characterization of the miRNA population in a Fagacea species, and the comparative analysis of miRNA expression in leaf and phellem libraries represents an important step to uncovering specific regulatory networks controlling phellem differentiation.     

Genome-wide methylation profiling identifies novel methylated genes in neuroblastoma tumors

Neuroblastoma is a very heterogeneous tumor of childhood. The clinical spectra range from very aggressive metastatic disease to spontaneous regression, even without therapy. Aberrant DNA methylation pattern is a common feature of most cancers. For neuroblastoma, it has been demonstrated both for single genes as well as genome-wide, where a so-called methylator phenotype has been described. Here, we present a study using Illumina 450K methylation arrays on 60 neuroblastoma tumors. We show that aggressive tumors, characterized by International Neuroblastoma Risk Group (INRG) as stage M, are hypermethylated compared to low-grade tumors. On the contrary, INRG stage L tumors display more non-CpG methylation. The genes with the highest number of hypermethylated CpG sites in INRG M tumors are TERT, PCDHGA4, DLX5, and DLX6-AS1. Gene ontology analysis showed a representation of neuronal tumor relevant gene functions among the differentially methylated genes. For validation, we used a set of independent tumors previously analyzed with the Illumina 27K methylation arrays, which confirmed the differentially methylated sites. Top candidate genes with aberrant methylation were analyzed for altered gene expression through the R2 platform (http://r2.amc.nl), and for correlations between methylation and gene expression in a public dataset. Altered expression in nonsurvivors was found for the genes B3GALT4 and KIAA1949, CLIC5, DLX6-AS, TERT, and PIRT, and strongest correlations were found for TRIM36, KIAA0513, and PIRT. Our data indicate that methylation profiling can be used for patient stratification and informs on epigenetically deregulated genes with the potential of increasing our knowledge about the underlying mechanisms of tumor development.     

Global profiling of miRNAs and the hairpin precursors: insights into miRNA processing and novel miRNA discovery

MicroRNAs (miRNAs) constitute an important class of small regulatory RNAs that are derived from distinct hairpin precursors (pre-miRNAs). In contrast to mature miRNAs, which have been characterized in numerous genome-wide studies of different organisms, research on global profiling of pre-miRNAs is limited. Here, using massive parallel sequencing, we have performed global characterization of both mouse mature and precursor miRNAs. In total, 87 369 704 and 252 003 sequencing reads derived from 887 mature and 281 precursor miRNAs were obtained, respectively. Our analysis revealed new aspects of miRNA/pre-miRNA processing and modification, including eight Ago2-cleaved pre-miRNAs, eight new instances of miRNA editing and exclusively 5′ tailed mirtrons. Furthermore, based on the sequences of both mature and precursor miRNAs, we developed a miRNA discovery pipeline, miRGrep, which does not rely on the availability of genome reference sequences. In addition to 239 known mouse pre-miRNAs, miRGrep predicted 41 novel ones with high confidence. Similar as known ones, the mature miRNAs derived from most of these novel loci showed both reduced abundance following Dicer knockdown and the binding with Argonaute2. Evaluation on data sets obtained from Caenorhabditis elegans and Caenorhabditis sp.11 demonstrated that miRGrep could be widely used for miRNA discovery in metazoans, especially in those without genome reference sequences.     

Integrated miRNA and mRNA expression profiling of mouse mammary tumor models identifies miRNA signatures associated with mammary tumor lineage

Background  

MicroRNAs (miRNAs) are small, non-coding, endogenous RNAs involved in regulating gene expression and protein translation. miRNA expression profiling of human breast cancers has identified miRNAs related to the clinical diversity of the disease and potentially provides novel diagnostic and prognostic tools for breast cancer therapy. In order to further understand the associations between oncogenic drivers and miRNA expression in sub-types of breast cancer, we performed miRNA expression profiling on mammary tumors from eight well-characterized genetically engineered mouse (GEM) models of human breast cancer, including MMTV-H-Ras, -Her2/neu, -c-Myc, -PymT, -Wnt1 and C3(1)/SV40 T/t-antigen transgenic mice, BRCA1 ^fl/fl ;p53 ^+/-;MMTV-cre knock-out mice and the p53 ^fl/fl ;MMTV-cre transplant model.     

Metabolic and miRNA profiling of TMV infected plants reveals biphasic temporal changes

Plant viral infections induce changes including gene expression and metabolic components. Identification of metabolites and microRNAs (miRNAs) differing in abundance along infection may provide a broad view of the pathways involved in signaling and defense that orchestrate and execute the response in plant-pathogen interactions. We used a systemic approach by applying both liquid and gas chromatography coupled to mass spectrometry to determine the relative level of metabolites across the viral infection, together with a miRs profiling using a micro-array based procedure. Systemic changes in metabolites were characterized by a biphasic response after infection. The first phase, detected at one dpi, evidenced the action of a systemic signal since no virus was detected systemically. Several of the metabolites increased at this stage were hormone-related. miRs profiling after infection also revealed a biphasic alteration, showing miRs alteration at 5 dpi where no virus was detected systemically and a late phase correlating with virus accumulation. Correlation analyses revealed a massive increase in the density of correlation networks after infection indicating a complex reprogramming of the regulatory pathways, either in response to the plant defense mechanism or to the virus infection itself. Our data propose the involvement of a systemic signaling on early miRs alteration.     

Massive analysis of cDNA Ends (MACE) and miRNA expression profiling identifies proatherogenic pathways in chronic kidney disease

Epigenetic dysregulation contributes to the high cardiovascular disease burden in chronic kidney disease (CKD) patients. Although microRNAs (miRNAs) are central epigenetic regulators, which substantially affect the development and progression of cardiovascular disease (CVD), no data on miRNA dysregulation in CKD-associated CVD are available until now. We now performed high-throughput miRNA sequencing of peripheral blood mononuclear cells from ten clinically stable hemodialysis (HD) patients and ten healthy controls, which allowed us to identify 182 differentially expressed miRNAs (e.g., miR-21, miR-26b, miR-146b, miR-155). To test biological relevance, we aimed to connect miRNA dysregulation to differential gene expression. Genome-wide gene expression profiling by MACE (Massive Analysis of cDNA Ends) identified 80 genes to be differentially expressed between HD patients and controls, which could be linked to cardiovascular disease (e.g., KLF6, DUSP6, KLF4), to infection / immune disease (e.g., ZFP36, SOCS3, JUND), and to distinct proatherogenic pathways such as the Toll-like receptor signaling pathway (e.g., IL1B, MYD88, TICAM2), the MAPK signaling pathway (e.g., DUSP1, FOS, HSPA1A), and the chemokine signaling pathway (e.g., RHOA, PAK1, CXCL5). Formal interaction network analysis proved biological relevance of miRNA dysregulation, as 68 differentially expressed miRNAs could be connected to 47 reciprocally expressed target genes. Our study is the first comprehensive miRNA analysis in CKD that links dysregulated miRNA expression with differential expression of genes connected to inflammation and CVD. After recent animal data suggested that targeting miRNAs is beneficial in experimental CVD, our data may now spur further research in the field of CKD-associated human CVD.     

Global gene expression profiling reveals similarities and differences among mouse pluripotent stem cells of different origins and strains

Pluripotent stem cell lines with similar phenotypes can be derived from both blastocysts (embryonic stem cells, ESC) and primordial germ cells (embryonic germ cells, EGC). Here, we present a compendium DNA microarray analysis of multiple mouse ESCs and EGCs from different genetic backgrounds (strains 129 and C57BL/6) cultured under standard conditions and in differentiation-promoting conditions by the withdrawal of Leukemia Inhibitory Factor (LIF) or treatment with retinoic acid (RA). All pluripotent cell lines showed similar gene expression patterns, which separated them clearly from other tissue stem cells with lower developmental potency. Differences between pluripotent lines derived from different sources (ESC vs. EGC) were smaller than differences between lines derived from different mouse strains (129 vs. C57BL/6). Even in the differentiation-promoting conditions, these pluripotent cells showed the same general trends of gene expression changes regardless of their origin and genetic background. These data indicate that ESCs and EGCs are indistinguishable based on global gene expression patterns alone. On the other hand, a detailed comparison between a group of ESC lines and a group of EGC lines identified 20 signature genes whose average expression levels were consistently higher in ESC lines, and 84 signature genes whose average expression levels were consistently higher in EGC lines, irrespective of mouse strains. Similar analysis identified 250 signature genes whose average expression levels were consistently higher in a group of 129 cell lines, and 337 signature genes whose average expression levels were consistently higher in a group of C57BL/6 cell lines. Although none of the genes was exclusively expressed in either ESCs versus EGCs or 129 versus C57BL/6, in combination these signature genes provide a reliable separation and identification of each cell type. Differentiation-promoting conditions also revealed some minor differences between the cell lines. For example, in the presence of RA, EGCs showed a lower expression of muscle- and cardiac-related genes and a higher expression of gonad-related genes than ESCs. Taken together, the results provide a rich source of information about the similarities and differences between ESCs and EGCs as well as 129 lines and C57BL/6 lines. Such information will be crucial to our understanding of pluripotent stem cells. The results also underscore the importance of studying multiple cell lines from different strains when making comparisons based on gene expression analysis.