Analysis of gene expression has undergone a technological revolution. What was impossible 6 years ago is now routine. High-throughput DNA sequencing machines capable of generating hundreds of millions of reads allow, indeed force, a major revision toward the study of the genome’s functional output—the transcriptome. This thesis examines the history of DNA sequencing, measurement of gene expression by sequencing, isoform complexity driven by alternative splicing and mammalian piRNA precursor biogenesis. Examination of these topics is framed around development of a novel RNA-templated DNA-DNA ligation assay (SeqZip) that allows for efficient analysis of abundant, complex, and functional long RNAs. The discussion focuses on the future of transcriptome analysis, development and applications of SeqZip, and challenges presented to biomedical researchers by extremely large and rich datasets.

Transposon silencing is required to maintain genome stability. The non-coding piRNAs effectively suppress of transposon activity during germline development. In the Drosophila female germline, long precursors of piRNAs are transcribed from discrete heterochromatic clusters and then processed into primary piRNAs in the perinuclear nuage. However, the detailed mechanism of piRNA biogenesis, specifically how the nuclear and cytoplasmic processes are connected, is not well understood. The nuclear DEAD box protein UAP56 has been previously implicated in protein-coding gene transcript splicing and export. I have identified a novel function of UAP56 in piRNA biogenesis. In Drosophila egg chambers, UAP56 co-localizes with the cluster-associated HP1 variant Rhino. Nuage is a germline-specific perinuclear structure rich in piRNA biogenesis proteins, including Vasa, a DEAD box with an established role in piRNA production. Vasa-containing nuage granules localize directly across the nuclear envelope from cluster foci containing UAP56 and Rhino, and cluster transcripts immunoprecipitate with both Vasa and UAP56. Significantly, a charge-substitution mutation that alters a conserved surface residue in UAP56 disrupts co-localization with Rhino, germline piRNA production, transposon silencing, and perinuclear localization of Vasa. I therefore propose that UAP56 and Vasa function in a piRNA-processing compartment that spans the nuclear envelope.

Small RNAs are single-stranded, 18–36 nucleotide RNAs that can be categorized as miRNA, siRNA, and piRNA. miRNA are expressed ubiquitously in tissues and at particular developmental stages. They fine-tune gene expression by regulating the stability and translation of mRNAs. piRNAs are mainly expressed in the animal gonads and their major function is repressing transposable elements to ensure the faithful transfer of genetic information from generation to generation. My thesis research focused on the biogenesis of miRNAs and piRNAs using both experimental and computational strategies. The biogenesis of miRNAs involves sequential processing of their precursors by the RNase III enzymes Drosha and Dicer to generate miRNA/miRNA* duplexes, which are subsequently loaded into Argonaute proteins to form the RNA-induced silencing complex (RISC). We discovered that, after assembled into Ago1, more than a quarter of Drosophila miRNAs undergo 3′ end trimming by the 3′-to-5′ exoribonuclease Nibbler. Such trimming occurs after removal of the miRNA* strand from pre-RISC and may be the final step in RISC assembly, ultimately enhancing target messenger RNA repression. Moreover, by developing a specialized Burrow-Wheeler Transform based short reads aligner, we discovered that in the absence of Nibbler a subgroup of miRNAs undergoes increased tailing—non-templated nucleotide addition to their 3′ ends, which are usually associated with miRNA degradation. Therefore, the 3′ trimming by Nibbler might increase miRNA stability by protecting them from degradation. In Drosophila germ line, piRNAs associate with three PIWI-clade Argonaute proteins, Piwi, Aub, and Ago3. piRNAs bound by Aub and Ago3 are generated by reciprocal cleavages of sense and antisense transposon transcripts (a.k.a., the “Ping-Pong” cycle), which amplifies piRNA abundance and degrades transposon transcripts in the cytoplasm. On the other hand, Piwi and its associated piRNA repress the transcription of transposons in the nucleus. We discovered that Aub- and Ago3-mediated transposon RNA cleavage not only generates piRNAs bound to each other, but also produces substrates for the endonuclease Zucchini, which processively cleaves those substrates in a periodicity of ~26 nt and generates piRNAs that predominantly load into Piwi. Without Aub or Ago3, the abundance of Piwi-bound piRNAs drops and transcriptional silencing is compromised. Our discovery revises the current model of piRNA biogenesis.

PIWI-interacting RNAs (piRNA) are a group of 23–35 nucleotide (nt) short RNAs that protect animal gonads from transposon activities. In Drosophila germ line, piRNAs can be categorized into two different categories— primary and secondary piRNAs— based on their origins. Primary piRNAs, generated from transcripts of specific genomic regions called piRNA clusters, which are enriched in transposon fragments that are unlikely to retain transposition activity. The transcription and maturation of primary piRNAs from those cluster transcripts are poorly understood. After being produced, a group of primary piRNAs associates Piwi proteins and directs them to repress transposons at the transcriptional level in the nucleus. Other than their direct role in repressing transposons, primary piRNAs can also initiate the production of secondary piRNA. piRNAs with such function are loaded in a second PIWI protein named Aubergine (Aub). Similar to Piwi, Aub is guided by piRNAs to identify its targets through base-pairing. Differently, Aub functions in the cytoplasm by cleaving transposon mRNAs. The 5' cleavage products are not degraded but loaded into the third PIWI protein Argonaute3 (Ago3). It is believed that an unidentified nuclease trims the 3' ends of those cleavage products to 23–29 nt, becoming mature piRNAs remained in Ago3. Such piRNAs whose 5' ends are generated by another PIWI protein are named secondary piRNAs. Intriguingly, secondary piRNAs loaded into Ago3 also cleave transposon mRNA or piRNA cluster transcripts and produce more secondary piRNAs loaded into Aub. This reciprocal feed-forward loop, named the “Ping-Pong cycle”, amplified piRNA abundance. By dissecting and analyzing data from large-scale deep sequencing of piRNAs and transposon transcripts, my dissertation research elucidates the biogenesis of germline piRNAs in Drosophila. How primary piRNAs are processed into mature piRNAs remains enigmatic. I discover that primary piRNA signal on the genome display a fixed periodicity of ~26 nt. Such phasing depends on Zucchini, Armitage and some other primary piRNA pathway components. Further analysis suggests that secondary piRNAs bound to Ago3 can initiate phased primary piRNA production from cleaved transposon RNAs. The first ~26 nt becomes a secondary piRNA that bind Aub while the subsequent piRNAs bind Piwi, allowing piRNAs to spread beyond the site of RNA cleavage. This discovery adds sequence diversity to the piRNA pool, allowing adaptation to changes in transposon sequence. We further find that most Piwi-associated piRNAs are generated from the cleavage products of Ago3, instead of being processed from piRNA cluster transcripts as the previous model suggests. The cardinal function of Ago3 is to produce antisense piRNAs that direct transcriptional silencing by Piwi, rather to make piRNAs that guide post-transcriptional silencing by Aub. Although Ago3 slicing is required to efficiently trigger phased piRNA production, an alternative, slicing-independent pathway suffices to generate Piwi-bound piRNAs that repress transcription of a subset of transposon families. The alternative pathway may help flies silence newly acquired transposons for which they lack extensively complementary piRNAs. The Ping-Pong model depicts that first ten nucleotides of Aub-bound piRNAs are complementary to the first ten nt of Ago3-bound piRNAs. Supporting this view, piRNAs bound to Aub typically begin with Uridine (1U), while piRNAs bound to Ago3 often have adenine at position 10 (10A). Furthermore, the majority of Ping-Pong piRNAs form this 1U:10A pair. The Ping-Pong model proposes that the 10A is a consequence of 1U. By statistically quantifying those target piRNAs not paired to g1U, we discover that 10A is not directly caused by 1U. Instead, fly Aub as well as its homologs, Siwi in silkmoth and MILI in mice, have an intrinsic preference for adenine at the t1 position of their target RNAs. On the other hand, this t1A (and g10A after loading) piRNA directly give rise to 1U piRNA in the next Ping-Pong cycle, maximizing the affinity between piRNAs and PIWI proteins.

Os peptídeos antimicrobianos desempenham papéis protetores críticos em uma gama de doenças humanas, incluindo o câncer. Vários estudos demonstraram funções - tais como proliferação, angiogênese, apoptose e imunomodulação - desses peptídeos em vias cancerígenas cruciais. Investigamos o papel do Peptídeo antimicrobiano LL-37 sobre stemness em câncer de mama (SKBR3) e células de melanoma (A375). Análise por PCR array da expressão diferencial de genes em SKBR3 e A375 com knockdown por siRNA para o mRNA de LL-37 revelou uma regulação negativa de genes relacionados com stemness, incluindo transcriptase reversa da telomerase, forkhead box D3 e para o fator indiferenciado de transcrição de células embrionárias 1, notavelmente em células de câncer de mama.Além disso, as células SKBR3 com knockdown para a expressão de LL-37 mostraram uma diminuição da produção de oncosferas em comparação com controles negativos, enquanto as células A375 exibiram uma produção aumentada. Tomados em conjunto, nossos achados indicam um papel para LL- 37 em stemness, dependendo do tipo de celular analisado / Antimicrobial peptides play critical protective roles in a range of human diseases, including cancer. Multiple studies have demonstrated functions -- such as proliferation, angiogenesis, apoptosis and immunomodulation -- of these peptides in crucial cancer pathways. We investigated the role of the antimicrobial peptide LL-37 on stemness in breast cancer (SKBR3) and melanoma cells (A375). PCR array analysis of differential gene expression in SKBR3 and A375 cancer cell lines downregulated for LL-37 expression by siRNA revealed downregulation of genes related to stemness, including telomerase reverse transcriptase, forkhead box D3 and undifferentiated embryonic cell transcription factor 1, remarkably in breast cancer cells. Furthermore, SKBR3 cells knocked down for LL-37 expression showed a decreased production of oncospheres in comparison with negative controls, while A375 cells exhibited increased production. Taken collectively, our findings indicate a role for LL-37 in cancer cell stemness depending on the cell type

Os peptídeos antimicrobianos desempenham papéis protetores críticos em uma gama de doenças humanas, incluindo o câncer. Vários estudos demonstraram funções - tais como proliferação, angiogênese, apoptose e imunomodulação - desses peptídeos em vias cancerígenas cruciais. Investigamos o papel do Peptídeo antimicrobiano LL-37 sobre stemness em câncer de mama (SKBR3) e células de melanoma (A375). Análise por PCR array da expressão diferencial de genes em SKBR3 e A375 com knockdown por siRNA para o mRNA de LL-37 revelou uma regulação negativa de genes relacionados com stemness, incluindo transcriptase reversa da telomerase, forkhead box D3 e para o fator indiferenciado de transcrição de células embrionárias 1, notavelmente em células de câncer de mama.Além disso, as células SKBR3 com knockdown para a expressão de LL-37 mostraram uma diminuição da produção de oncosferas em comparação com controles negativos, enquanto as células A375 exibiram uma produção aumentada. Tomados em conjunto, nossos achados indicam um papel para LL- 37 em stemness, dependendo do tipo de celular analisado / Antimicrobial peptides play critical protective roles in a range of human diseases, including cancer. Multiple studies have demonstrated functions -- such as proliferation, angiogenesis, apoptosis and immunomodulation -- of these peptides in crucial cancer pathways. We investigated the role of the antimicrobial peptide LL-37 on stemness in breast cancer (SKBR3) and melanoma cells (A375). PCR array analysis of differential gene expression in SKBR3 and A375 cancer cell lines downregulated for LL-37 expression by siRNA revealed downregulation of genes related to stemness, including telomerase reverse transcriptase, forkhead box D3 and undifferentiated embryonic cell transcription factor 1, remarkably in breast cancer cells. Furthermore, SKBR3 cells knocked down for LL-37 expression showed a decreased production of oncospheres in comparison with negative controls, while A375 cells exhibited increased production. Taken collectively, our findings indicate a role for LL-37 in cancer cell stemness depending on the cell type

Skeletal muscle is a complicated and heterogeneous striated muscle tissue that serves critical mechanical and metabolic functions in the organism. The process of generating skeletal muscle, myogenesis, is elaborately coordinated by members of the protein kinase family, which transmit diverse signals initiated by extracellular stimuli to myogenic transcriptional hierarchy in muscle cells. Mitogen-activated protein kinases (MAPKs) including p38 MAPK, c-Jun N terminal kinase (JNK) and extracellular signal-regulated protein kinase (ERK) are components of serine/threonine protein kinase cascades that play important roles in skeletal muscle differentiation. The exploration of MAPK upstream kinases identified mitogen activated protein kinase kinase kinase kinase 4 (MAP4K4), a serine/threonine protein kinase that modulates p38 MAPK, JNK and ERK activities in multiple cell lines. Our lab further discovered that Map4k4 regulates peroxisome proliferator-activated receptor γ (PPARγ) translation in cultured adipocytes through inactivating mammalian target of rapamycin (mTOR), which controls skeletal muscle differentiation and hypotrophy in kinase-dependent and -independent manners. These findings suggest potential involvement of Map4k4 in skeletal myogenesis. Therefore, for the first part of my thesis, I characterize the role of Map4k4 in skeletal muscle differentiation in cultured muscle cells. Here I show that Map4k4 functions as a myogenic suppressor mainly at the early stage of skeletal myogenesis with a moderate effect on myoblast fusion during late-stage muscle differentiation. In agreement, Map4k4 expression and protein kinase activity are declined with myogenic differentiation. The inhibitory effect of Map4k4 on skeletal myogenesis requires its kinase activity. Surprisingly, none of the identified Map4k4 downstream effectors including p38 MAPK, JNK and ERK is involved in the Map4k4-mediated myogenic differentiation. Instead, expression of myogenic regulatory factor Myf5, a positive mediator of skeletal muscle differentiation is transiently regulated by Map4k4 to partially control skeletal myogenesis. Mechanisms by which Map4k4 modulates Myf5 amount have yet to be determined. In the second part of my thesis, I assess the relationship between Map4k4 and IGF-mediated signaling pathways. Although siRNA-mediated silencing of Map4k4 results in markedly enhanced myotube formation that is identical to the IGF-induced muscle hypertrophic phenotype, and Map4k4 regulates IGF/Akt signaling downstream effector mTOR in cultured adipocytes, Map4k4 appears not to be involved in the IGF-mediated ERK1/2 signaling axis and the IGF-mediated Akt signaling axis in C2C12 myoblasts. Furthermore, Map4k4 does not affect endogenous Akt signaling or mTOR activity during C2C12 myogenic differentiation. The results presented here not only identify Map4k4 as a novel suppressor of skeletal muscle differentiation, but also add to our knowledge of Map4k4 action on multiple signaling pathways in muscle cells during skeletal myogenesis. The effects that Map4k4 exerts on myoblast differentiation, fusion and Myf5 expression implicate Map4k4 as a potential drug target for muscle mass growth, skeletal muscle regeneration and muscular dystrophy.