Abstract Ascoviruses (AVs) are members of the family Ascoviridae that are transmitted by female endoparasitic wasps and cause lethal infection in lepidopteran insects. AVs possess large double stranded DNA genomes ranging from 116-186 kbp. Recently, genomes of four AV species have been completely sequenced and have revealed important genes potentially needed for virus DNA replication and infection. Phylogenetic analyses of several of these genes indicate that AVs are closely related to iridoviruses and likely evolved from them. Two unique features, mode of transmission and cytopathology which involves cleavage of cells into virus-containing vesicles, make AVs different from other insect pathogenic viruses. During this decade, tremendous advancements in the study of RNA silencing mechanisms have openned a new dimension in virology. It is now evident that viruses reshape the cellular environment by reprogramming host RNA silencing machinery. The process of RNA silencing involves small non-coding RNAs, which with the help of nuclease-containing regulatory proteins bind to complementary messenger RNA (mRNA) targets, resulting in inhibition of gene expression. This sophisticated style of gene regulation has attained a fundamental status in living organisms, since RNA silencing has been revealed to be ubiquitous from viruses to prokaryotes to eukaryotes. Two main categories of small RNAs, short interfering RNA (siRNA) and microRNA (miRNA), have been defined as major players in RNA silencing. Interestingly, viral genomes like that of their hosts, encode miRNAs that can be used during virus invasion to manipulate host genes as well as miRNA biogenesis. Here, we report on the identification of the first insect virus miRNA (HvAVmiR- 1) derived from the major capsid protein (MCP) gene of Heliothis virescens ascovirus 7 (HvAV3e). HvAV-miR-1 expression was found to be strictly regulated and specifically detected from 96 h post-infection. HvAV-miR-1 expression coincides with a marked reduction of the expression of HvAV3e DNA polymerase I, which is a predicted target. Ectopic expression of the full-length and truncated versions of MCP retaining the miRNA sequence significantly reduced DNA polymerase I transcript levels and inhibited viral replication. Our results indicate that HvAV-miR-1 directs degradation of DNA polymerase I transcripts and regulates replication of HvAV3e. Further, we investigated changes in the expression levels of host miRNAs upon HvAV3e infection in an insect cell line derived from Helicoverpa zea fat body and investigated the role of a host miRNA, Hz-miR24, in the hostvirus system. It was found that Hz-miR24 is differentially expressed following virus infection, with an increase in its expression levels late in infection. Functional analyses demonstrated that Hz-miR24 targets viral DNA-dependent RNA polymerase and its β subunit mRNAs. This was confirmed using ectopic expression of Hz-miR24 and a green fluorescent protein-based reporter system. Expression of the target gene was substantially enhanced in cells transfected with a synthesized inhibitor of Hz-miR24. These findings suggest that ascoviruses encode their own miRNA(s) and concurrently manipulate host miRNAs that in turn regulate the expression of their genes at specific time points after infection. In connection to RNA silencing, we characterized a ribonuclease III (RNase III) protein encoded by HvAV3e. We found that RNase III protein was functional in vivo as well as in vitro and catalyzed long and short double stranded RNAs. Expression analyses during virus infection revealed autoregulation of this protein by degradation of its RNA transcripts. Moreover, RNase III protein was found to be involved in suppression of RNA silencing and essential for virus DNA replication and infection. Finally, we studied another ascoviral 8 protein, a putative inhibitor of apoptosis (IAP), which was found to be essential for virus DNA replication and pathology. Further, despite inhibition of apoptosis by HvAV3e, the IAP-like protein was found dispensable for the inhibition of replication. In conclusion, for successful invasion and attenuation of host antiviral responses, ascoviruses seem to utilize viral encoded proteins as well as miRNAs. Since the genomes of these viruses have only recently been sequenced, the role of many of the encoded genes essential for pathogenesis and manipulation of antiviral defence mechanisms remains to be eluciated.

Abstract Ascoviruses (AVs) are members of the family Ascoviridae that are transmitted by female endoparasitic wasps and cause lethal infection in lepidopteran insects. AVs possess large double stranded DNA genomes ranging from 116-186 kbp. Recently, genomes of four AV species have been completely sequenced and have revealed important genes potentially needed for virus DNA replication and infection. Phylogenetic analyses of several of these genes indicate that AVs are closely related to iridoviruses and likely evolved from them. Two unique features, mode of transmission and cytopathology which involves cleavage of cells into virus-containing vesicles, make AVs different from other insect pathogenic viruses. During this decade, tremendous advancements in the study of RNA silencing mechanisms have openned a new dimension in virology. It is now evident that viruses reshape the cellular environment by reprogramming host RNA silencing machinery. The process of RNA silencing involves small non-coding RNAs, which with the help of nuclease-containing regulatory proteins bind to complementary messenger RNA (mRNA) targets, resulting in inhibition of gene expression. This sophisticated style of gene regulation has attained a fundamental status in living organisms, since RNA silencing has been revealed to be ubiquitous from viruses to prokaryotes to eukaryotes. Two main categories of small RNAs, short interfering RNA (siRNA) and microRNA (miRNA), have been defined as major players in RNA silencing. Interestingly, viral genomes like that of their hosts, encode miRNAs that can be used during virus invasion to manipulate host genes as well as miRNA biogenesis. Here, we report on the identification of the first insect virus miRNA (HvAVmiR- 1) derived from the major capsid protein (MCP) gene of Heliothis virescens ascovirus 7 (HvAV3e). HvAV-miR-1 expression was found to be strictly regulated and specifically detected from 96 h post-infection. HvAV-miR-1 expression coincides with a marked reduction of the expression of HvAV3e DNA polymerase I, which is a predicted target. Ectopic expression of the full-length and truncated versions of MCP retaining the miRNA sequence significantly reduced DNA polymerase I transcript levels and inhibited viral replication. Our results indicate that HvAV-miR-1 directs degradation of DNA polymerase I transcripts and regulates replication of HvAV3e. Further, we investigated changes in the expression levels of host miRNAs upon HvAV3e infection in an insect cell line derived from Helicoverpa zea fat body and investigated the role of a host miRNA, Hz-miR24, in the hostvirus system. It was found that Hz-miR24 is differentially expressed following virus infection, with an increase in its expression levels late in infection. Functional analyses demonstrated that Hz-miR24 targets viral DNA-dependent RNA polymerase and its β subunit mRNAs. This was confirmed using ectopic expression of Hz-miR24 and a green fluorescent protein-based reporter system. Expression of the target gene was substantially enhanced in cells transfected with a synthesized inhibitor of Hz-miR24. These findings suggest that ascoviruses encode their own miRNA(s) and concurrently manipulate host miRNAs that in turn regulate the expression of their genes at specific time points after infection. In connection to RNA silencing, we characterized a ribonuclease III (RNase III) protein encoded by HvAV3e. We found that RNase III protein was functional in vivo as well as in vitro and catalyzed long and short double stranded RNAs. Expression analyses during virus infection revealed autoregulation of this protein by degradation of its RNA transcripts. Moreover, RNase III protein was found to be involved in suppression of RNA silencing and essential for virus DNA replication and infection. Finally, we studied another ascoviral 8 protein, a putative inhibitor of apoptosis (IAP), which was found to be essential for virus DNA replication and pathology. Further, despite inhibition of apoptosis by HvAV3e, the IAP-like protein was found dispensable for the inhibition of replication. In conclusion, for successful invasion and attenuation of host antiviral responses, ascoviruses seem to utilize viral encoded proteins as well as miRNAs. Since the genomes of these viruses have only recently been sequenced, the role of many of the encoded genes essential for pathogenesis and manipulation of antiviral defence mechanisms remains to be eluciated.

Double-stranded RNA (dsRNA), a double helix formed by two antiparallel complementary RNA strands, is a unique structure with a variety of biological effects. dsRNA can be introduced into the cell from exogenous sources or it can be produced endogenously. There are four basic mechanisms producing dsRNA: inverted repeat transcription, convergent transcription, pairing of sense and antisense RNAs produced in trans, and RNA dependent RNA polymerase-mediated synthesis dsRNA. Different mechanisms of production determine additional structural features of dsRNA, such as dsRNA termini, mismatches etc. These features may affect cellular response to dsRNA. Recognition of dsRNA can trigger several responses that act in sequence-specific or sequence-independent manners. The main sequence- specific response triggered by dsRNA is RNA interference (RNAi) is. Our laboratory has been studying mechanism of induction of RNAi in mammalian cells using one specific type of long dsRNA expression system. The dsRNA used in these experiments formed hairpin structure with long 5' and 3' single-strand RNA overhangs. We hypothesized that other dsRNA substrates might be more efficient than the one used in mammalian RNAi experiments since 2002. Accordingly, the main aim of my thesis was to compare efficiency of different dsRNA...

Oncolytic viruses are novel candidates for cancer therapy and their efficacy relies on their capacity to overcome the host’s anti-viral barriers. In mammalian cells, the anti-viral response involves a protein-signaling cascade known as the interferon pathway, which alerts the immune system and limits the propagation of infection. Given that most cancer cells have defects in this pathway, they are susceptible to viral infection and responsive to oncolytic virotherapy. For reasons that remain unknown, many cancers are still refractory to oncolytic viruses, which suggests the existence of additional antiviral mechanisms. In this study, we investigate the potential involvement of an alternative antiviral pathway in cancer cells. Given that insects and plants rely on the RNA silencing pathway for their anti-viral protection, we investigated the presence of a similar mechanism in cancer cells. We found viral genome-derived small RNAs in various cancer cell lines upon infection, which is indicative of an RNA-mediated antiviral response. Also, various viruses encode suppressors of the RNA interference pathway. To determine if an oncolytic virus could benefit from such a factor, we engineered an oncolytic virus variant to encode the Nodamura virus B2 protein, a known inhibitor of RNA silencing-mediated immune responses. Using this virus, we observed enhanced cytotoxicity in 33 out of the 38 human cancer cell lines tested. Furthermore, our results show inhibition of viral genome cleavage and altered microRNA processing by our B2-expressing oncolytic virus. Taken together, our data suggests the blockade of RNA silencing antiviral pathways and/or antiviral microRNA processing improves the efficacy of our B2-encoding virus in a cell-line specific manner. Overall, our results establish the improved potential of our novel virus therapy and demonstrate for the first time the involvement of RNA pathways in the antiviral defense of cancer cells.

No description available.

No description available.

Ο καρκίνος του πνεύμονα είναι η πιο συχνή αιτία θανάτου από καρκίνο παγκοσμίως. Είναι γνωστό ότι ο καρκίνος του πνεύμονα είναι διαδικασία πολλαπλών σταδίων, στην οποία ενέχονται γενετικοί και επιγενετικοί μηχανισμοί. Ενεργοποίηση ογκογονιδίων συμβαίνει σε όλους τους βρογχοπνευμονικούς καρκίνους με αποτέλεσμα την αύξηση των μιτογόνων σημάτων. Στον καρκίνο του πνεύμονα τα πιο συχνά ενεργοποιημένα ογκογονίδια είναι τα EGFR, ERbB2, MYC, KRAS, MET, CCND1, CDK4, EML4-ALK fusion, και BCL2. Επίσης, η απώλεια ογκοκατασταλτικών γονιδίων είναι ιδιαίτερα σημαντική στην πνευμονική καρκινογένεση και είναι συνήθως αποτέλεσμα απενεργοποίησης και των δυο αλληλόμορφων. Συχνά απενεργοποιημένα ογκοκατασταλτικά γονίδια στον καρκίνο του πνεύμονα είναι TP53, RB1, STK11, CDKN2A, FHIT και PTEN. Οι επιγενετικοί μηχανισμοί περιλαμβάνουν την μεθυλίωση του DNA, την τροποποίηση των ιστονών και τη ρύθμιση της γονιδιακής έκφρασης μέσω των microRNAs. Τα microRNAs είναι μικρά, μη κωδικοποιούντα μόρια RNA που εμπλέκονται στην αρνητική μετα-μεταγραφική ρύθμιση της έκφρασης των γονιδίων. Μελέτες έχουν αποδείξει το ρόλο των miRNAs στην φυσιολογική πνευμονική ανάπτυξη και ομοιόσταση αλλά και τον ενεργό ρολό τους στην παθογένεια πνευμονικών νοσημάτων όπως είναι ο καρκίνος του πνεύμονα. Η δημιουργία ωρίμων, λειτουργικών microRNAs απαιτεί τη συντονισμένη δράση μιας ομάδας πρωτεϊνών που στο σύνολο τους απαρτίζουν το μηχανισμό ρύθμισης των microRNA (microRNA machinery). Ο μηχανισμός ελέγχου των microRNA ρυθμίζει μέσω των παραγομένων microRNAs την έκφραση πολλών ογκοκατασταλτικών γονιδίων και ογκογονιδίων. Κύρια συστατικά του μηχανισμού ρύθμισης των microRNA είναι οι ριβονουκλεάσες Drosha, Dicer και AGO2. Σκοπός της παρούσας διατριβής ήταν η μελέτη της κυτταρικής εντόπισης και έκφρασης των συστατικών του μηχανισμού ρύθμισης των microRNA, Drosha, Dicer και AGO2, στον μη μικροκυτταρικό καρκίνο του πνεύμονα. Συγκεκριμένα, ελέγχθηκε η κυτταρική εντόπιση των Drosha, Dicer και AGO2 στις κυτταρικές σειρές καρκίνου του πνεύμονα A549, H23, H358, H661, HCC827 με τη μέθοδο του ανοσοφθορισμού. Στις ίδιες κυτταρικές σειρές, μελετήθηκαν τα κυτταρικά επίπεδα των πρωτεϊνών Drosha, Dicer και AGO2 με την μέθοδο της SDS-PAGE και του ανοσοαποτυπώματος. H έκφραση των πρωτεϊνών αυτών μελετήθηκε σε ιστολογικές τομές παραφίνης μη μικροκυτταρικού καρκίνου του πνεύμονα NSCLC με την μέθοδο της ανοσοϊστοχημείας. Επιπλέον συσχετίσαμε τα επίπεδα της ανοσοϊστοχημικής χρώσης αυτών των ριβονουκλεασών με κλινικοπαθολογοανατομικές παραμέτρους. Η παρούσα εργασία είναι η πρώτη που μελετά την κυτταρική εντόπιση της Drosha in vitro και σε ιστούς από ανθρώπινο καρκίνο του πνεύμονα. Τα επίπεδα ανοσοέκφρασης της Drosha ήταν στατιστικά χαμηλότερα στα νεοπλασματικά κύτταρα NSCLC, σε σχέση με τα φυσιολογικά. Επίσης, τα κυτταρικά επίπεδα της Drosha ήταν στατιστικά χαμηλότερα στα NSCLC σταδίου Ι σε σχέση με το φυσιολογικό ιστό. Όμως, στατιστικά σημαντική διαφορά δεν προέκυψε από την σύγκριση καρκινικών ιστών μεταξύ τους κατά ιστολογικό τύπο, στάδιο νόσου και βαθμό κακοήθειας. Τα ευρήματα αυτά υποδηλώνουν συμμετοχή της ριβονουκλεάσης Drosha στην πνευμονική κακοήθη εξαλλαγή και στην παθογένεια του NSCLC αλλά όχι στην εξέλιξη της νόσου. Η παρούσα εργασία είναι η πρώτη που μελετά την κυτταρική εντόπιση της Dicer in vitro και σε ιστούς από ανθρώπινο καρκίνο του πνεύμονα. Τα πειράματα ανοσοϊστοχημείας, ανέδειξαν ότι τα επίπεδα ανοσοέκφρασης της Dicer ήταν στατιστικά χαμηλότερα στα NSCLC σταδίου Ι σε σχέση με το φυσιολογικό ιστό (p=0,040). Μάλιστα, παρατηρήθηκε στατιστικά σημαντική διαφορά στην ανοσοέκφραση της Dicer στην σύγκριση των τριών σταδίων μεταξύ τους (p=0,049) και αυτό το εύρημα παρουσιάζεται για πρώτη φορά στη βιβλιογραφία από την παρούσα μελέτη. Όμως, τα κυτταρικά επίπεδα αυτής της πρωτεΐνης δεν σχετίζονται με τον ιστολογικό τύπο αλλά και το βαθμό της κακοήθειας. Τα ευρήματα μας αυτά εισηγούνται τη συμμετοχή της ριβονουκλεάσης Dicer στην πνευμονική καρκινογένεση και στην εξέλιξη της νόσου. Τέλος, τα κυτταρικά επίπεδα της ενδονουκλεάσης AGO2 είναι στατιστικά χαμηλότερα στα πνευμονικά νεοπλασματικά κύτταρα σε σχέση με τα φυσιολογικά. Η πρωτεϊνική έκφραση των κυτταρικών σειρών NSCLC παρουσίασε σχεδόν ομοιόμορφη κατανομή. Μάλιστα, και για την πρωτεΐνη AGO2 τα επίπεδα ανοσοέκφρασης είναι στατιστικά χαμηλότερα στα NSCLC σταδίου Ι σε σχέση με το φυσιολογικό ιστό (p=0,000). Όμως, παρατηρήθηκε ότι τα κυτταρικά επίπεδα αυτής της πρωτεΐνης δεν σχετίζονται με τον ιστολογικό τύπο, το στάδιο της νόσου αλλά και το βαθμό της κακοήθειας. Το γεγονός αυτό ενισχύει την άποψη ότι η AGO2 συμμετέχει στην παθοβιολογία του NSCLC αλλά πιθανά όχι στην εξέλιξη της νόσου. Εάν αποδειχθεί σημαντική η συμμετοχή του μηχανισμού ρύθμισης των microRNA στην παθογένεια της πνευμονικής κακοήθειας, θα υπάρξει η δυνατότητα να χρησιμοποιηθούν για την δημιουργία υποομάδων («μοριακά πορτραίτα») του καρκίνου του πνεύμονα, οι οποίες να έχουν προγνωστική αλλά και θεραπευτική αξία (στοχευμένες θεραπείες). / Lung cancer is the leading cause of cancer related death worldwide. Decades of research have contributed to our understanding that lung cancer is a multistep process involving genetic and epigenetic alterations. Oncogene activation occurs in all lung cancers, resulting in persistent upregulation of mitogenic signals. In lung cancer commonly activated oncogenes are EGFR, ERbB2, MYC, KRAS, MET, CCND1, CDK4, EML4-ALK fusion, and BCL2. Loss of tumor suppressor gene (TSG) function is also important in lung carcinogenesis and usually results from silencing of both alleles. Commonly unactivated TSGs in lung cancer are TP53, RB1, STK11, CDKN2A, FHIT and PTEN. Epigenetic alterations include DNA methylation, histone modification and microRNA regulation of gene expression. MicroRNAs are small non-protein encoding RNAs, responsible for the negative post transcriptional regulation of gene expression. Studies have shown the role of microRNAs in normal pulmonary development and homeostasis but also in the pathogenesis of multiple lung diseases including lung cancer. The biogenesis of mature and functional microRNAs requires the orchestrated action of a group of proteins, collectively refered to as miRNA machinery. The miRNA machinery regulates the expression of many TSGs and oncogenes in a miRNA guided fashion. Drosha, Dicer and AGO2 are main components of the miRNA machinery. Our study adressed the cellular localization and protein levels of Drosha, Dicer and AGO2, components of the miRNA machinery, in NSCLC cell lines, and in NSCLC FFPE tissue sections. We employed immunofluorescence and Western blot analysis in five NSCLC cell lines and immunohistochemistry on FFPE NSCLC tissue sections. Staining intensity of the FFPE tissues was correlated with clinicopathological parameters. Altered Drosha cellular distribution was evident in neoplasia. The staining intensity of Drosha (p=0,03) was significantly lower in neoplastic tissues compared to normal tissues. When we compared neoplastic tissue stage I with normal tissues, Drosha’s staining intensity (p=0,002) was significantly lower. Drosha, protein levels were not significantly associated with age, tumor histology, grade or stage. Altered Dicer nuclear distribution was evident in lung neoplasia. The staining intensity of Dicer was significantly lower in neoplastic tissues stage I compared to normal tissues (p=0,04). Dicer’s protein levels in FFPE tissues were significantly associated with tumor stage (p=0,049). AGO2 excibited physiological cytoplasmic distribution in lung neoplasia. The staining intensity of AGO2 was significantly lower in neoplastic tissues compared to normal tissues (p=0,000). When we compared neoplastic tissue stage I with normal tissues, AGO2 staining intensity (p=0,000) was significantly lower. AGO2 protein levels were not significantly associated with age, tumor histology, grade or stage. Our findings provide evidence that the miRNA machinery components Drosha, Dicer and AGO2 are involved in lung carcinogenesis but only Dicer is implicated in cancer progression. The expression levels of the miRNA processing components might contribute to improved cancerous molecular portraits for achieving personalized medicine, the selection of patient-tailored treatment regimens.

Les microARN (miARN) sont des petits ARN non-codants pour des protéines qui permettent d’inhiber la traduction d’ARN messagers. Pour obtenir un miARN, un gène de miARN passe à travers une voie de maturation dans laquelle il sera coupé à deux reprises par les enzymes Drosha et Dicer. Pour interagir avec les enzymes, les gènes de miARN possèdent une structure générale en tige-boucle. Cependant, les détails de cette structure sont encore peu connus. L’objectif principal de ce projet était d’établir s’il y a une relation entre l’efficacité de maturation et la dynamique de la structure. Pour cela, l’efficacité de maturation de plusieurs variants de miARN a été évaluée par Northern Blot. La dynamique de la structure a été mesurée par un programme informatique à partir de l’information de la séquence. Une corrélation de 0,74 avec une valeur p de 0,02206 a été obtenue entre les dynamiques et les ratios d’efficacités de maturation de miARN. Cette corrélation est supérieure à celle obtenue basée sur l’énergie libre des structures prédites les plus stables qui n’atteignent pas 0,6. Les mutants de miR128-1 et miR188 ont été découverts comme diminuant la maturation. De plus, les mutants de miR125a, miR188 et miR330 affectent le site de clivage de Drosha. Une meilleure connaissance de la dynamique de l’ARN impliquée dans la maturation permettrait de définir l’impact des mutations dans les séquences de miARN ou encore de prédire les séquences pouvant générer des miARN. / MicroRNAs (miRNAs) are small non-coding RNAs, which can inhibit target messenger RNAs translation. In order to obtain a miRNA, two enzymes, Drosha and Dicer cut the gene of miRNA. The RNA interacts with the proteins by its general hairpin structure. However, the details of the structure are still missing. The objective of this project is to establish if there is a relation between the efficiency of maturation and the RNA’s structural dynamics. In order to do this, the maturation efficiency of miRNA variants is measured by Northern Blot. The structural dynamics is measured by a program assessing the information of the sequence. The correlation between the dynamics and the maturation efficiency of the miRNA is 0.74 with a p-value of 0.02206. This correlation is superior to those based on free energy, which does not reach 0.6. The tested mutants of miR128-1 and miR188 have inhibited maturation; also, those of miR125a, miR188 and miR330 have modified the cleavage site of Drosha. A better knowledge of the dynamic structure involved in maturation would help define the impact of miRNA mutation or to predict sequences that are able to generate miRNAs.

Host organisms evolved antiviral responses, which can recognize the viral infection and deal with it. One of the frequent signs of viral infection in a cell is appearance of double-stranded RNA (dsRNA). One of the pathways responding to dsRNA is RNA interference (RNAi), which functions as the key antiviral defence system in invertebrates and plants. Mammals, however, utilize for antiviral defence a different dsRNA-sensing pathway called the interferon response. RNAi functions only in mammalian oocytes and early embryonal stages although its enzymatic machinery is present in all somatic cells, where it is employed in the microRNA pathway. A previous study indicated that the functionality of RNAi in mouse oocytes functions due to an oocyte-specific isoform of protein Dicer (DicerO ), which is truncated at the N-terminus. In my thesis, I aimed to assess whether DicerO processes RNAi substrates more efficiently in vitro than the full-length Dicer (DicerS ), which is found in somatic cells. Therefore, I developed Dicer purification protocol for obtaining both recombinant mouse Dicer isoforms of high purity. I examined their activity in a non-radioactive cleavage assay using RNA substrates with structural features characteristic of RNAi substrates. My results suggest that recombinant DicerO and DicerS do not...

microRNAs (miRNAs) are small (20-23 nt), non-coding single stranded RNA molecules that play an important role in post-transcriptional regulation of protein-coding genes. miRNAs have been found in all animal lineages, and have been implicated as critical regulators during development in multiple species. The echinoderms, Strongylocentrotus purpuratus (sea urchin) and Patiria miniata (sea star) are excellent model organisms for studying development due to their well-characterized transcriptional gene networks, ease of working with their embryos in the laboratory and phylogenetic position as invertebrate deuterostomes. Literature on miRNAs in echinoderm embryogenesis is limited. It has been shown that RNAi genes are developmentally expressed and regulated in sea urchin embryos, but no study in the sea urchin has examined the expression of miRNAs. The goal of my work has been to study miRNA regulation in echinoderm developmental gene networks. I have identified developmentally regulated miRNAs in sea urchin and sea star embryos, using a combination of computational and wet lab experimental techniques. I developed a probabilistic model (named HHMMiR) based on hierarchical hidden Markov models (HHMMs) to classify genomic hairpins into miRNA precursors and random stem-loop structures. I then extended this model to make an efficient decoder by introduction of explicit state duration densities. We used the Illumina Genome Analyzer to sequence small RNA libraries in mixed stage population of embryos from one to three days after fertilization of S. purpuratus and P. miniata. We developed a computational pipeline for analysis of these miRNAseq data to reveal the miRNA populations in both species, and study their differential expression. We also used northern blots and whole mount in situ hybridization experimental techniques to study the temporal and spatial expression patterns of some of these miRNAs in sea urchin embryos. By knocking down the major components of the miRNA biogenesis pathway, we studied the global effects of miRNAs on embryo morphology and differentiation genes. The biogenesis genes selected for this purpose are the RNAse III enzyme, Dicer and Argonaute. Dicer is necessary for the processing of mature miRNAs from hairpin structures while Ago is a necessary part of the RISC (RNA interference silencing complex) assembly, which is required for the miRNA to hybridize to its target mRNA site. Knocking down these genes hinders normal development of the sea urchin embryo and leads to loss of the larval skeleton, a novel phenotype not seen in sea stars, as well as abnormal gastrulation. Comparison of differentiation gene marker expression between control and Ago knocked down sea urchin embryos shows interesting patterns of expansion and suppression of adjoining some embryonic territories, while ingression of larval skeletogenesis progenitors does not occur.