Trichoderma harzianum expressed sequence tags for identification of genes with putative roles in mycoparasitism against Fusarium solani

The plant pathogen Fusarium solani causes a disease root rot of common bean (Phaseolus vulgaris) resulting in great losses of yield in irrigated areas of the Southeast and Midwest regions of Brazil. Species of the genus Trichoderma have been used in the biological control of this pathogen as an alternative to chemical control. To gain new insights into the biocontrol mechanism used by Trichoderma harzianum against the phytopathogenic fungus, Fusarium solani, we performed a transcriptome analysis using expressed sequence tags (ESTs) and quantitative real-time PCR (RT-qPCR) approaches. A cDNA library from T. harzianum mycelium (isolate ALL42) grown on cell walls of F. solani (CWFS) was constructed and analyzed. A total of 2927 high quality sequences were selected from 3845 and 37.7% were identified as unique genes. The Gene Ontology analysis revealed that the majority of the annotated genes are involved in metabolic processes (80.9%), followed by cellular process (73.7%). We tested twenty genes that encode proteins with potential role in biological control. RT-qPCR analysis showed that none of these genes were expressed when T. harzianum was challenged with itself. These genes showed different patterns of expression during in vitro interaction between T. harzianum and F. solani. (C) 2012 Elsevier Inc. All rights reserved.

Dissecting Phaseolus vulgaris Innate Immune System against Colletotrichum lindemuthianum Infection

Background: The genus Colletotrichum is one of the most economically important plant pathogens, causing anthracnose on a wide range of crops including common beans (Phaseolus vulgaris L.). Crop yield can be dramatically decreased depending on the plant cultivar used and the environmental conditions. This study aimed to identify potential genetic components of the bean immune system to provide environmentally friendly control measures against this fungus. Methodology and Principal Findings: As the common bean is not amenable to reverse genetics to explore functionality and its genome is not fully curated, we used putative Arabidopsis orthologs of bean expressed sequence tag (EST) to perform bioinformatic analysis and experimental validation of gene expression to identify common bean genes regulated during the incompatible interaction with C. lindemuthianum. Similar to model pathosystems, Gene Ontology (GO) analysis indicated that hormone biosynthesis and signaling in common beans seem to be modulated by fungus infection. For instance, cytokinin and ethylene responses were up-regulated and jasmonic acid, gibberellin, and abscisic acid responses were down-regulated, indicating that these hormones may play a central role in this pathosystem. Importantly, we have identified putative bean gene orthologs of Arabidopsis genes involved in the plant immune system. Based on experimental validation of gene expression, we propose that hypersensitive reaction as part of effector-triggered immunity may operate, at least in part, by down-regulating genes, such as FLS2-like and MKK5-like, putative orthologs of the Arabidopsis genes involved in pathogen perception and downstream signaling. Conclusions/Significance: We have identified specific bean genes and uncovered metabolic processes and pathways that may be involved in the immune response against pathogens. Our transcriptome database is a rich resource for mining novel defense-related genes, which enabled us to develop a model of the molecular components of the bean innate immune system regulated upon pathogen attack.

The Aurora A and B kinases are up-regulated in bone marrow-derived chronic lymphocytic leukemia cells and represent potential therapeutic targets

Background The malignant B cells in chronic lymphocytic leukemia receive signals from the bone marrow and lymph node microenvironments which regulate their survival and proliferation. Characterization of these signals and the pathways that propagate them to the interior of the cell is important for the identification of novel potential targets for therapeutic intervention. Design and Methods We compared the gene expression profiles of chronic lymphocytic leukemia B cells purified from bone marrow and peripheral blood to identify genes that are induced by the bone marrow microenvironment. Two of the differentially expressed genes were further studied in cell culture experiments and in an animal model to determine whether they could represent appropriate therapeutic targets in chronic lymphocytic leukemia. Results Functional classification analysis revealed that the majority of differentially expressed genes belong to gene ontology categories related to cell cycle and mitosis. Significantly up-regulated genes in bone marrow-derived tumor cells included important cell cycle regulators, such as Aurora A and B, survivin and CDK6. Down-regulation of Aurora A and B by RNA interference inhibited proliferation of chronic lymphocytic leukemia-derived cell lines and induced low levels of apoptosis. A similar effect was observed with the Aurora kinase inhibitor VX-680 in primary chronic lymphocytic leukemia cells that were induced to proliferate by CpG-oligonucleotides and interleukin-2. Moreover, VX-680 significantly blocked leukemia growth in a mouse model of chronic lymphocytic leukemia. Conclusions Aurora A and B are up-regulated in proliferating chronic lymphocytic leukemia cells and represent potential therapeutic targets in this disease.

Changes in the expression of proteins associated with aerobic glycolysis and cell migration are involved in tumorigenic ability of two glioma cell lines

Background: The most frequent and malignant brain cancer is glioblastoma multiforme (GBM). In gliomas, tumor progression and poor prognosis are associated with the tumorigenic ability of the cells. U87MG cells (wild-type p53) are known to be tumorigenic in nude mice, but T98G cells (mutant p53) are not tumorigenic. We investigated the proteomic profiling of these two cell lines in order to gain new insights into the mechanisms that may be involved in tumorigenesis. Results: We found 24 differentially expressed proteins between T98G and U87MG cells. Gene Ontology supports the notion that over-representation of differentially expressed proteins is involved in glycolysis, cell migration and stress oxidative response. Among those associated with the glycolysis pathway, TPIS and LDHB are up-regulated in U87MG cells. Measurement of glucose consumption and lactate production suggests that glycolysis is more effective in U87MG cells. On the other hand, G6PD expression was 3-fold higher in T98G cells and this may indicate a shift to the pentose-phosphate pathway. Moreover, GRP78 expression was also three-fold higher in T98G than in U87MG cells. Under thapsigargin treatment both cell lines showed increased GRP78 expression and the effect of this agent was inversely correlated to cell migration. Quantitative RT-PCR and immunohistochemistry of GRP78 in patient samples indicated a higher level of expression of GRP78 in grade IV tumors compared to grade I and non-neoplastic tissues, respectively. Conclusions: Taken together, these results suggest an important role of proteins involved in key functions such as glycolysis and cell migration that may explain the difference in tumorigenic ability between these two glioma cell lines and that may be extrapolated to the differential aggressiveness of glioma tumors.

Proteomic analysis of Porphyromonas gingivalis exposed to nicotine and cotinine

Cogo K, de Andrade A, Labate CA, Bergamaschi CC, Berto LA, Franco GCN, Goncalves RB, Groppo FC. Proteomic analysis of Porphyromonas gingivalis exposed to nicotine and cotinine. J Periodont Res 2012; 47: 766775. (c) 2012 John Wiley & Sons A/S Background and Objective: Smokers are more predisposed than nonsmokers to infection with Porphyromonas gingivalis, one of the most important pathogens involved in the onset and development of periodontitis. It has also been observed that tobacco, and tobacco derivatives such as nicotine and cotinine, can induce modifications to P. gingivalis virulence. However, the effect of the major compounds derived from cigarettes on expression of protein by P.gingivalis is poorly understood. Therefore, this study aimed to evaluate and compare the effects of nicotine and cotinine on the P.gingivalis proteomic profile. Material and Methods: Total proteins of P gingivalis exposed to nicotine and cotinine were extracted and separated by two-dimensional electrophoresis. Proteins differentially expressed were successfully identified through liquid chromatography-mass spectrometry and primary sequence databases using MASCOT search engine, and gene ontology was carried out using DAVID tools. Results: Of the approximately 410 protein spots that were reproducibly detected on each gel, 23 were differentially expressed in at least one of the treatments. A particular increase was seen in proteins involved in metabolism, virulence and acquisition of peptides, protein synthesis and folding, transcription and oxidative stress. Few proteins showed significant decreases in expression; those that did are involved in cell envelope biosynthesis and proteolysis and also in metabolism. Conclusion: Our results characterized the changes in the proteome of P.gingivalis following exposure to nicotine and cotinine, suggesting that these substances may modulate, with minor changes, protein expression. The present study is, in part, a step toward understanding the potential smokepathogen interaction that may occur in smokers with periodontitis.

Pilot survey of expressed sequence tags (ESTs) from the asexual blood stages of Plasmodium vivax in human patients

Abstract Background Plasmodium vivax is the most widely distributed human malaria, responsible for 70–80 million clinical cases each year and large socio-economical burdens for countries such as Brazil where it is the most prevalent species. Unfortunately, due to the impossibility of growing this parasite in continuous in vitro culture, research on P. vivax remains largely neglected. Methods A pilot survey of expressed sequence tags (ESTs) from the asexual blood stages of P. vivax was performed. To do so, 1,184 clones from a cDNA library constructed with parasites obtained from 10 different human patients in the Brazilian Amazon were sequenced. Sequences were automatedly processed to remove contaminants and low quality reads. A total of 806 sequences with an average length of 586 bp met such criteria and their clustering revealed 666 distinct events. The consensus sequence of each cluster and the unique sequences of the singlets were used in similarity searches against different databases that included P. vivax, Plasmodium falciparum, Plasmodium yoelii, Plasmodium knowlesi, Apicomplexa and the GenBank non-redundant database. An E-value of <10-30 was used to define a significant database match. ESTs were manually assigned a gene ontology (GO) terminology Results A total of 769 ESTs could be assigned a putative identity based upon sequence similarity to known proteins in GenBank. Moreover, 292 ESTs were annotated and a GO terminology was assigned to 164 of them. Conclusion These are the first ESTs reported for P. vivax and, as such, they represent a valuable resource to assist in the annotation of the P. vivax genome currently being sequenced. Moreover, since the GC-content of the P. vivax genome is strikingly different from that of P. falciparum, these ESTs will help in the validation of gene predictions for P. vivax and to create a gene index of this malaria parasite.

Genome mapping and expression analyses of human intronic noncoding RNAs reveal tissue-specific patterns and enrichment in genes related to regulation of transcription

Abstract Background RNAs transcribed from intronic regions of genes are involved in a number of processes related to post-transcriptional control of gene expression. However, the complement of human genes in which introns are transcribed, and the number of intronic transcriptional units and their tissue expression patterns are not known. Results A survey of mRNA and EST public databases revealed more than 55,000 totally intronic noncoding (TIN) RNAs transcribed from the introns of 74% of all unique RefSeq genes. Guided by this information, we designed an oligoarray platform containing sense and antisense probes for each of 7,135 randomly selected TIN transcripts plus the corresponding protein-coding genes. We identified exonic and intronic tissue-specific expression signatures for human liver, prostate and kidney. The most highly expressed antisense TIN RNAs were transcribed from introns of protein-coding genes significantly enriched (p = 0.002 to 0.022) in the 'Regulation of transcription' Gene Ontology category. RNA polymerase II inhibition resulted in increased expression of a fraction of intronic RNAs in cell cultures, suggesting that other RNA polymerases may be involved in their biosynthesis. Members of a subset of intronic and protein-coding signatures transcribed from the same genomic loci have correlated expression patterns, suggesting that intronic RNAs regulate the abundance or the pattern of exon usage in protein-coding messages. Conclusion We have identified diverse intronic RNA expression patterns, pointing to distinct regulatory roles. This gene-oriented approach, using a combined intron-exon oligoarray, should permit further comparative analysis of intronic transcription under various physiological and pathological conditions, thus advancing current knowledge about the biological functions of these noncoding RNAs.

Genes involved in thoracic exoskeleton formation during the pupal-to-adult molt in a social insect model, Apis mellifera

Background: The insect exoskeleton provides shape, waterproofing, and locomotion via attached somatic muscles. The exoskeleton is renewed during molting, a process regulated by ecdysteroid hormones. The holometabolous pupa transforms into an adult during the imaginal molt, when the epidermis synthe3sizes the definitive exoskeleton that then differentiates progressively. An important issue in insect development concerns how the exoskeletal regions are constructed to provide their morphological, physiological and mechanical functions. We used whole-genome oligonucleotide microarrays to screen for genes involved in exoskeletal formation in the honeybee thoracic dorsum. Our analysis included three sampling times during the pupal-to-adult molt, i.e., before, during and after the ecdysteroid-induced apolysis that triggers synthesis of the adult exoskeleton. Results: Gene ontology annotation based on orthologous relationships with Drosophila melanogaster genes placed the honeybee differentially expressed genes (DEGs) into distinct categories of Biological Process and Molecular Function, depending on developmental time, revealing the functional elements required for adult exoskeleton formation. Of the 1,253 unique DEGs, 547 were upregulated in the thoracic dorsum after apolysis, suggesting induction by the ecdysteroid pulse. The upregulated gene set included 20 of the 47 cuticular protein (CP) genes that were previously identified in the honeybee genome, and three novel putative CP genes that do not belong to a known CP family. In situ hybridization showed that two of the novel genes were abundantly expressed in the epidermis during adult exoskeleton formation, strongly implicating them as genuine CP genes. Conserved sequence motifs identified the CP genes as members of the CPR, Tweedle, Apidermin, CPF, CPLCP1 and Analogous-to-Peritrophins families. Furthermore, 28 of the 36 muscle-related DEGs were upregulated during the de novo formation of striated fibers attached to the exoskeleton. A search for cis-regulatory motifs in the 5′-untranslated region of the DEGs revealed potential binding sites for known transcription factors. Construction of a regulatory network showed that various upregulated CP- and muscle-related genes (15 and 21 genes, respectively) share common elements, suggesting co-regulation during thoracic exoskeleton formation. Conclusions: These findings help reveal molecular aspects of rigid thoracic exoskeleton formation during the ecdysteroid-coordinated pupal-to-adult molt in the honeybee.

Transcriptional profiling of Saccharomyces cerevisiae exposed to propolis

Abstract Background Propolis is a natural product of plant resins collected by honeybees (Apis mellifera) from various plant sources. Our previous studies indicated that propolis sensitivity is dependent on the mitochondrial function and that vacuolar acidification and autophagy are important for yeast cell death caused by propolis. Here, we extended our understanding of propolis-mediated cell death in the yeast Saccharomyces cerevisiae by applying systems biology tools to analyze the transcriptional profiling of cells exposed to propolis. Methods We have used transcriptional profiling of S. cerevisiae exposed to propolis. We validated our findings by using real-time PCR of selected genes. Systems biology tools (physical protein-protein interaction [PPPI] network) were applied to analyse the propolis-induced transcriptional bevavior, aiming to identify which pathways are modulated by propolis in S. cerevisiae and potentially influencing cell death. Results We were able to observe 1,339 genes modulated in at least one time point when compared to the reference time (propolis untreated samples) (t-test, p-value 0.01). Enrichment analysis performed by Gene Ontology (GO) Term finder tool showed enrichment for several biological categories among the genes up-regulated in the microarray hybridization such as transport and transmembrane transport and response to stress. Real-time RT-PCR analysis of selected genes showed by our microarray hybridization approach was capable of providing information about S. cerevisiae gene expression modulation with a considerably high level of confidence. Finally, a physical protein-protein (PPPI) network design and global topological analysis stressed the importance of these pathways in response of S. cerevisiae to propolis and were correlated with the transcriptional data obtained thorough the microarray analysis. Conclusions In summary, our data indicate that propolis is largely affecting several pathways in the eukaryotic cell. However, the most prominent pathways are related to oxidative stress, mitochondrial electron transport chain, vacuolar acidification, regulation of macroautophagy associated with protein target to vacuole, cellular response to starvation, and negative regulation of transcription from RNA polymerase II promoter. Our work emphasizes again the importance of S. cerevisiae as a model system to understand at molecular level the mechanism whereby propolis causes cell death in this organism at the concentration herein tested. Our study is the first one that investigates systematically by using functional genomics how propolis influences and modulates the mRNA abundance of an organism and may stimulate further work on the propolis-mediated cell death mechanisms in fungi.

New insights in osteogenic differentiation revealed by mass spectrometric assessment of phosphorylated substrates in murine skin mesenchymal cells

Abstract Background Bone fractures and loss represent significant costs for the public health system and often affect the patients quality of life, therefore, understanding the molecular basis for bone regeneration is essential. Cytokines, such as IL-6, IL-10 and TNFα, secreted by inflammatory cells at the lesion site, at the very beginning of the repair process, act as chemotactic factors for mesenchymal stem cells, which proliferate and differentiate into osteoblasts through the autocrine and paracrine action of bone morphogenetic proteins (BMPs), mainly BMP-2. Although it is known that BMP-2 binds to ActRI/BMPR and activates the SMAD 1/5/8 downstream effectors, little is known about the intracellular mechanisms participating in osteoblastic differentiation. We assessed differences in the phosphorylation status of different cellular proteins upon BMP-2 osteogenic induction of isolated murine skin mesenchymal stem cells using Triplex Stable Isotope Dimethyl Labeling coupled with LC/MS. Results From 150 μg of starting material, 2,264 proteins were identified and quantified at five different time points, 235 of which are differentially phosphorylated. Kinase motif analysis showed that several substrates display phosphorylation sites for Casein Kinase, p38, CDK and JNK. Gene ontology analysis showed an increase in biological processes related with signaling and differentiation at early time points after BMP2 induction. Moreover, proteins involved in cytoskeleton rearrangement, Wnt and Ras pathways were found to be differentially phosphorylated during all timepoints studied. Conclusions Taken together, these data, allow new insights on the intracellular substrates which are phosphorylated early on during differentiation to BMP2-driven osteoblastic differentiation of skin-derived mesenchymal stem cells.

Expression analysis and in silico characterization of intronic long noncoding RNAs in renal cell carcinoma: emerging functional associations

Abstract Background Intronic and intergenic long noncoding RNAs (lncRNAs) are emerging gene expression regulators. The molecular pathogenesis of renal cell carcinoma (RCC) is still poorly understood, and in particular, limited studies are available for intronic lncRNAs expressed in RCC Methods Microarray experiments were performed with custom-designed arrays enriched with probes for lncRNAs mapping to intronic genomic regions. Samples from 18 primary RCC tumors and 11 nontumor adjacent matched tissues were analyzed. Meta-analyses were performed with microarray expression data from three additional human tissues (normal liver, prostate tumor and kidney nontumor samples), and with large-scale public data for epigenetic regulatory marks and for evolutionarily conserved sequences. Results A signature of 29 intronic lncRNAs differentially expressed between RCC and nontumor samples was obtained (false discovery rate (FDR) <5%). A signature of 26 intronic lncRNAs significantly correlated with the RCC five-year patient survival outcome was identified (FDR <5%, p-value ≤0.01). We identified 4303 intronic antisense lncRNAs expressed in RCC, of which 22% were significantly (p <0.05) cis correlated with the expression of the mRNA in the same locus across RCC and three other human tissues. Gene Ontology (GO) analysis of those loci pointed to 'regulation of biological processes’ as the main enriched category. A module map analysis of the protein-coding genes significantly (p <0.05) trans correlated with the 20% most abundant lncRNAs, identified 51 enriched GO terms (p <0.05). We determined that 60% of the expressed lncRNAs are evolutionarily conserved. At the genomic loci containing the intronic RCC-expressed lncRNAs, a strong association (p <0.001) was found between their transcription start sites and genomic marks such as CpG islands, RNA Pol II binding and histones methylation and acetylation. Conclusion Intronic antisense lncRNAs are widely expressed in RCC tumors. Some of them are significantly altered in RCC in comparison with nontumor samples. The majority of these lncRNAs is evolutionarily conserved and possibly modulated by epigenetic modifications. Our data suggest that these RCC lncRNAs may contribute to the complex network of regulatory RNAs playing a role in renal cell malignant transformation.

Profilo genomico e di espressione dell'Osteosarcoma. Identificazione dei geni con alto grado di correlazione tra numero di copie ed espressione

To identify the regions of recurrent copy number abnormality in osteosarcoma and their effect on gene expression, we performed an integrated genome-wide high-resolution array CGH (aCGH) and gene expression profiling analysis on 40 human OS tissues and 12 OS cell lines. This analysis identified several recurrent chromosome regions that contain genes that show a gene dosage effect on gene expression. A further search, performed on those genes that were over-expressed and localized in the frequently amplified chromosomal regions, greatly reduced the number of candidate genes while their characterization using gene ontology (GO) analysis suggests the importance of the deregulation of the G1-to-S phase in the development of the disease. We also identified frequent deletions on 3q in the vicinity of LSAMP and performed a fine mapping analysis of the breakpoints. We precisely mapped the breakpoints in several instances and demonstrated that the majority do not involve the LSAMP gene itself, and that they appear to form by a process of non-homologous end joining. In addition, aCGH analysis revealed frequent gains of IGF1R that were highly correlated with its protein level. Blockade of IGF1R in OS cell lines with high copy number gain led to growth inhibition suggesting that IGF1R may be a viable drug target in OS, particularly in patients with copy number driven overexpression of this receptor.

A clustering method for robust and reliable large scale functional and structural protein sequence annotation

Bioinformatics, in the last few decades, has played a fundamental role to give sense to the huge amount of data produced. Obtained the complete sequence of a genome, the major problem of knowing as much as possible of its coding regions, is crucial. Protein sequence annotation is challenging and, due to the size of the problem, only computational approaches can provide a feasible solution. As it has been recently pointed out by the Critical Assessment of Function Annotations (CAFA), most accurate methods are those based on the transfer-by-homology approach and the most incisive contribution is given by cross-genome comparisons. In the present thesis it is described a non-hierarchical sequence clustering method for protein automatic large-scale annotation, called “The Bologna Annotation Resource Plus” (BAR+). The method is based on an all-against-all alignment of more than 13 millions protein sequences characterized by a very stringent metric. BAR+ can safely transfer functional features (Gene Ontology and Pfam terms) inside clusters by means of a statistical validation, even in the case of multi-domain proteins. Within BAR+ clusters it is also possible to transfer the three dimensional structure (when a template is available). This is possible by the way of cluster-specific HMM profiles that can be used to calculate reliable template-to-target alignments even in the case of distantly related proteins (sequence identity < 30%). Other BAR+ based applications have been developed during my doctorate including the prediction of Magnesium binding sites in human proteins, the ABC transporters superfamily classification and the functional prediction (GO terms) of the CAFA targets. Remarkably, in the CAFA assessment, BAR+ placed among the ten most accurate methods. At present, as a web server for the functional and structural protein sequence annotation, BAR+ is freely available at http://bar.biocomp.unibo.it/bar2.0.

Transkriptom des Meristemgewebes des Kolumnarapfels

Kolumnare Apfelbäume (Malus x domestica) stellen aufgrund ihres auffälligen Phänotyps eine ökonomisch interessante Wuchsform dar. Diese extreme Form des Kurztriebwuchses zeichnet sich durch einen insgesamt sehr schlanken, säulenförmigen Habitus aus, welcher eine dichte Pflanzung und damit einhergehend Ertragssteigerungen im Vergleich zu normalwüchsigen Bäumen ermöglicht. Verursacht wird der Phänotyp durch die Anwesenheit eines einzelnen, dominanten Allels des Columnar (Co)-Gens. Bis auf die approximative Lokalisation des Gens auf Chromosom 10 ist über mögliche Identität und Funktion bislang nichts bekannt.rnIn der vorliegenden Arbeit wurde ein erster Versuch unternommen, mit Hilfe von Next Generation Sequencing (NGS) Technologien und RNA-Seq Einblicke in das Transkriptom des Sprossapikalmeristems (SAM) kolumnarer Apfelbäume zu gewinnen. So konnte gezeigt werden, dass unabhängig vom Zeitpunkt der Entnahme des Materials mehrere hundert Gene differentiell reguliert werden. Diese lassen sich funktional in mehrere überrepräsentierte Kategorien gruppieren, von denen sich einige wiederum mit dem kolumnaren Phänotyp assoziieren lassen. Durch den Einsatz weiterer Expressionsstudien (Microarrays, qRT-PCR) konnten frühere Ergebnisse bezüglich des Hormonhaushalts auf Genebene bestätigt und neue Erkenntnisse gewonnen werden, die eine mögliche Erklärung für den Phänotyp darstellen. Weiterhin ergab der Vergleich aller durchgeführten Expressionsstudien eine Anreicherung signifikant differentiell regulierter Gene auf Chromosom 10, was auf einen „selective sweep“ hindeutet. Eine potentielle epigenetische Regulation dieser Gene durch das Genprodukt von Co könnte daher möglich sein. Mehr als die Hälfte dieser Gene lassen sich darüber hinaus aufgrund ihrer Funktion direkt mit dem kolumnaren Phänotyp assoziieren.rnDiese Ergebnisse zeigen, dass die Anwesenheit des Co-Allels massive Veränderungen in der Genregulation des SAMs mit sich bringt, wobei einige dieser differentiell regulierten Gene mit großer Wahrscheinlichkeit an der Etablierung des kolumnaren Phänotyps beteiligt sind. Auch wenn die Funktion des Co-Genproduktes nicht abschließend geklärt werden konnte, sind doch anhand der Resultate schlüssige Hypothesen diesbezüglich möglich.rn

Genexpressionsanalyse boviner mesenchymaler Stammzellen und deren in-vitro-differenzierten Folgelinien

Mesenchymale Stamzellen (MSC) sind Vertreter der adulten Stammzellen. Sie bergen durch ihre große Plastizität ein immenses Potential für die klinische Nutzung in Form von Stammzelltherapien. Zellen dieses Typs kommen vornehmlich im Knochenmark der großen Röhrenknochen vor und können zu Knochen, Knorpel und Fettzellen differenzieren. MSC leisten einen wichtigen Beitrag im Rahmen regenerativer Prozesse, beispielsweise zur Heilung von Frakturen. Breite Studien demonstrieren bereits jetzt auch bei komplexeren Erkrankungen (z.B. Osteoporose) therapeutisch vielversprechende Einsatzmöglichkeiten. Oft kommen hierbei aus MSC gezielt differenzierte Folgelinien aus Zellkulturen zum Einsatz. Dies bedingt eine kontrollierte Steuerung der Differenzierungsprozesse in vitro. Der Differenzierung einer Stammzelle liegt eine komplexe Veränderung ihrer Genexpression zugrunde. Genexpressionsmuster zur Erhaltung und Proliferation der Stammzellen müssen durch solche, die der linienspezifischen Differenzierung dienen, ersetzt werden. Die mit der Differenzierung einhergehende, transkriptomische Neuausrichtung ist für das Verständnis der Prozesse grundlegend und wurde bislang nur unzureichend untersucht. Ziel der vorliegenden Arbeit ist eine transkriptomweite und vergleichende Genexpressionsanalyse Mesenchymaler Stammzellen und deren in vitro differenzierten Folgelinien mittels Plasmid - DNA Microarrays und Sequenziertechniken der nächsten Generation (RNA-Seq, Illumina Plattform). In dieser Arbeit diente das Hausrind (Bos taurus) als Modellorganismus, da es genetisch betrachtet eine hohe Ähnlichkeit zum Menschen aufweist und Knochenmark als Quelle von MSC gut verfügbar ist. Primärkulturen Mesenchymaler Stammzellen konnten aus dem Knochenmark von Rindern erfolgreich isoliert werden. Es wurden in vitro Zellkultur - Versuche durchgeführt, um die Zellen zu Osteoblasten, Chondrozyten und Adipozyten zu differenzieren. Zur Genexpressionsanalyse wurde RNA aus jungen MSC und einer MSC Langzeitkultur („alte MSC“), sowie aus den differenzierten Zelllinien isoliert und für nachfolgende Experimente wo nötig amplifiziert. Der Erfolg der Differenzierungen konnte anhand der Genexpression von spezifischen Markergenen und mittels histologischer Färbungen belegt werden. Hierbei zeigte sich die Differenzierung zu Osteoblasten und Adipozyten erfolgreich, während die Differenzierung zu Chondrozyten trotz diverser Modifikationen am Protokoll nicht erfolgreich durchgeführt werden konnte. Eine vergleichende Hybridisierung zur Bestimmung differentieller Genexpression (MSC vs. Differenzierung) mittels selbst hergestellter Plasmid - DNA Microarrays ergab für die Osteogenese mit Genen wie destrin und enpp1, für die undifferenzierten MSC mit dem Gen sema3c neue Kandidatengene, deren biologische Funktion aufzuklären in zukünftigen Experimenten vielversprechende Ergebnisse liefern sollte. Die Analyse der transkriptomweiten Genexpression mittels NGS lieferte einen noch umfangreicheren Einblick ins Differenzierungsgeschehen. Es zeigte sich eine hohe Ähnlichkeit im Expressionsprofil von jungen MSC und Adipozyten, sowie zwischen den Profilen der alten MSC (eine Langzeitkultur) und Osteoblasten. Die alten MSC wiesen deutliche Anzeichen für eine spontane Differenzierung in die osteogene Richtung auf. Durch Analyse der 100 am stärksten exprimierten Gene jeder Zelllinie ließen sich für junge MSC und Adipozyten besonders Gene der extrazellulären Matrix (z.B col1a1,6 ; fn1 uvm.) auffinden. Sowohl Osteoblasten, als auch die alten MSC exprimieren hingegen verstärkt Gene mit Bezug zur oxidativen Phosphorylierung, sowie ribosomale Proteine. Eine Betrachtung der differentiellen Genexpression (junge MSC vs. Differenzierung) mit anschließender Pathway Analyse und Genontologie Anreicherungsstatistik unterstützt diese Ergebnisse vor allem bei Osteoblasten, wo nun jedoch zusätzlich auch Gene zur Regulation der Knochenentwicklung und Mineralisierung in den Vordergrund treten. Für Adipozyten konnte mit Genen des „Jak-STAT signaling pathway“, der Fokalen Adhäsion, sowie Genen des „Cytokine-cytokine receptor interaction pathway“ sehr spannende Einsichten in die Biologie dieses Zelltyps erlangt werden, die sicher weiterer Untersuchungen bedürfen. In undifferenzierten MSC konnte durch differentielle Genexpressionsanalyse die Rolle des nicht kanonischen Teils des WNT Signalweges als für die Aufrechterhaltung des Stammzellstatus potentiell äußerst einflussreich ermittelt werden. Die hier diskutierten Ergebnisse zeigen beispielhaft, dass besonders mittels Genexpressionsanalyse im Hochdurchsatzverfahren wertvolle Einblicke in die komplexe Biologie der Stammzelldifferenzierung möglich sind. Als Grundlage für nachfolgende Arbeiten konnten interessante Gene ermittelt und Hypothesen zu deren Einfluss auf Stammzelleigenschaften und Differenzierungsprozesse aufgestellt werden. Um einen besseren Einblick in den Differenzierungsverlauf zu ermöglichen, könnten künftig NGS Analysen zu unterschiedlichen Differenzierungszeitpunkten durchgeführt werden. Zudem wären weitere Anstrengungen zur erfolgreichen Etablierung der chondrogenen Differenzierung zur vollständigen Analyse der Genexpression des trilinearen Differenzierungspotentials von MSC wünschenswert.