Highly accurate RNA and DNA sequencing: Application to longstanding questions in aging and cancer

Thesis (Ph.D.)--University of Washington, 2016-06

Presence of a novel small RNA-containing virus in a laboratory culture of the endoparasitic wasp Venturia canescens (Hymenoptera : Ichneumonidae)

Parasitoid wasps use a variety of mechanisms to alter their host's physiology to the benefit of the developing endoparasite inside the host larva. Association of certain wasps with viruses and virus-like particles (VLPs) that contribute to their success in parasitism is one of the fascinating evolutionary adaptations conferring active or passive protection for the endoparasite from the host immune system. Venturia canescens has been shown to produce VLPs that provide protection for the developing parasitoid egg inside the host, Ephestia kuehniella. Here, we report on the presence of a novel small RNA-containing virus from V. canescens, designated as VcSRV, occurring in the ovaries of the wasp. The virus particles are found together with VcVLPs in the lumen of the calyx region of the ovaries and are injected together with the egg and VcVLPs into E kuehniella larvae where they enter hemocytes. Alignment of the RNA-dependent RNA polymerase gene of VcSRV indicates that the virus most likely belongs to the recently described genus Iflavirus. (c) 2004 Elsevier Ltd. All rights reserved.

Comparative mechanistic studies of de novo RNA synthesis by flavivirus RNA-dependent RNA polymerases

Flavivirus protein NS5 harbors the RNA-dependent RNA polymerase (RdRp) activity. In contrast to the RdRps of hepaci- and pestiviruses, which belong to the same family of Flaviviridae, NS5 carries two activities, a methyltransferase (MTase) and a RdRp. RdRp domains of Dengue virus (DV) and West Nile virus (WNV) NS5 were purified in high yield relative to full-length NS5 and showed full RdRp activity. Steady-state enzymatic parameters were determined on homopolymeric template poly(rC). The presence of the MTase domain does not affect the RdRp activity. Flavivirus RdRp domains might bear more than one GTP binding site displaying positive cooperativity. The kinetics of RNA synthesis by four Flaviviridae RdRps were compared. In comparison to Hepatitis C RdRp, DV and WNV as well as Bovine Viral Diarrhea virus RdRps show less rate limitation by early steps of short-product fort-nation. This suggests that they display a higher conformational flexibility upon the transition from initiation to elongation. (c) 2006 Elsevier Inc. All rights reserved.

The dawn of the RNA World:toward functional complexity through ligation of random RNA oligomers

A main unsolved problem in the RNA World scenario for the origin of life is how a template-dependent RNA polymerase ribozyme emerged from short RNA oligomers obtained by random polymerization on mineral surfaces. A number of computational studies have shown that the structural repertoire yielded by that process is dominated by topologically simple structures, notably hairpin-like ones. A fraction of these could display RNA ligase activity and catalyze the assembly of larger, eventually functional RNA molecules retaining their previous modular structure: molecular complexity increases but template replication is absent. This allows us to build up a stepwise model of ligation- based, modular evolution that could pave the way to the emergence of a ribozyme with RNA replicase activity, step at which information-driven Darwinian evolution would be triggered. Copyright © 2009 RNA Society.

Regulation of ribosomal RNA expression across the lifespan is fine-tuned by maternal diet before implantation

Cells and organisms respond to nutrient deprivation by decreasing global rates of transcription, translation and DNA replication. To what extent such changes can be reversed is largely unknown. We examined the effect of maternal dietary restriction on RNA synthesis in the offspring. Low protein diet fed either throughout gestation or for the preimplantation period alone reduced cellular RNA content across fetal somatic tissues during challenge and increased it beyond controls in fetal and adult tissues after challenge release. Changes in transcription of ribosomal RNA, the major component of cellular RNA, were responsible for this phenotype as evidenced by matching alterations in RNA polymerase I density and DNA methylation at ribosomal DNA loci. Cellular levels of the ribosomal transcription factor Rrn3 mirrored the rRNA expression pattern. In cell culture experiments, Rrn3 overexpression reduced rDNA methylation and increased rRNA expression; the converse occurred after inhibition of Rrn3 activity. These observations define novel mechanism where poor nutrition before implantation irreversibly alters basal rates of rRNA transcription thereafter in a process mediated by rDNA methylation and Rrn3 factor.

A model for the modular evolution of rna addressing open questions on the origin of life

A main unsolved problem in the RNA world scenario for the origin of life is how a template-dependent RNA polymerase ribozyme emerged from short RNA oligomers generated by random polymerization of ribonucleotides (Joyce and Orgel 2006). Current estimates establish a minimum size about 165 nt long for such a ribozyme (Johnston et al. 2001), a length three to four times that of the longest RNA oligomers obtained by random polymerization on clay mineral surfaces (Huang and Ferris 2003, 2006). To overcome this gap, we have developed a stepwise model of ligation-based, modular evolution of RNA (Briones et al. 2009) whose main conceptual steps are summarized in Figure 1. This scenario has two main advantages with respect to previous hypotheses put forward for the origin of the RNA world: i) short RNA....

First complete genome sequence of a capsicum chlorosis tospovirus isolate from Australia with an unusually large S RNA intergenic region

The first complete genome sequence of capsicum chlorosis virus (CaCV) from Australia was determined using a combination of Illumina HiSeq RNA and Sanger sequencing technologies. Australian CaCV had a tripartite genome structure like other CaCV isolates. The large (L) RNA was 8913 nucleotides (nt) in length and contained a single open reading frame (ORF) of 8634 nt encoding a predicted RNA-dependent RNA polymerase (RdRp) in the viral-complementary (vc) sense. The medium (M) and small (S) RNA segments were 4846 and 3944 nt in length, respectively, each containing two non-overlapping ORFs in ambisense orientation, separated by intergenic regions (IGR). The M segment contained ORFs encoding the predicted non-structural movement protein (NSm; 927 nt) and precursor of glycoproteins (GP; 3366 nt) in the viral sense (v) and vc strand, respectively, separated by a 449-nt IGR. The S segment coded for the predicted nucleocapsid (N) protein (828 nt) and non-structural suppressor of silencing protein (NSs; 1320 nt) in the vc and v strand, respectively. The S RNA contained an IGR of 1663 nt, being the largest IGR of all CaCV isolates sequenced so far. Comparison of the Australian CaCV genome with complete CaCV genome sequences from other geographic regions showed highest sequence identity with a Taiwanese isolate. Genome sequence comparisons and phylogeny of all available CaCV isolates provided evidence for at least two highly diverged groups of CaCV isolates that may warrant re-classification of AIT-Thailand and CP-China isolates as unique tospoviruses, separate from CaCV.

STRUCTURAL ANALYSIS OF RIBOSOME BINDING ELEMENTS OF THE 3´ UTR IN TWO POSITIVE-STRAND PLANT RNA VIRUSES

Turnip crinkle virus (TCV) and Pea enation mosaic virus (PEMV) are two positive (+)-strand RNA viruses that are used to investigate the regulation of translation and replication due to their small size and simple genomes. Both viruses contain cap-independent translation elements (CITEs) within their 3´ untranslated regions (UTRs) that fold into tRNA-shaped structures (TSS) according to nuclear magnetic resonance and small angle x-ray scattering analysis (TCV) and computational prediction (PEMV). Specifically, the TCV TSS can directly associate with ribosomes and participates in RNA-dependent RNA polymerase (RdRp) binding. The PEMV kissing-loop TSS (kl-TSS) can simultaneously bind to ribosomes and associate with the 5´ UTR of the viral genome. Mutational analysis and chemical structure probing methods provide great insight into the function and secondary structure of the two 3´ CITEs. However, lack of 3-D structural information has limited our understanding of their functional dynamics. Here, I report the folding dynamics for the TCV TSS using optical tweezers (OT), a single molecule technique. My study of the unfolding/folding pathways for the TCV TSS has provided an unexpected unfolding pathway, confirmed the presence of Ψ3 and hairpin elements, and suggested an interconnection between the hairpins and pseudoknots. In addition, this study has demonstrated the importance of the adjacent upstream adenylate-rich sequence for the formation of H4a/Ψ3 along with the contribution of magnesium to the stability of the TCV TSS. In my second project, I report on the structural analysis of the PEMV kl-TSS using NMR and SAXS. This study has re-confirmed the base-pair pattern for the PEMV kl-TSS and the proposed interaction of the PEMV kl-TSS with its interacting partner, hairpin 5H2. The molecular envelope of the kl-TSS built from SAXS analysis suggests the kl-TSS has two functional conformations, one of which has a different shape from the previously predicted tRNA-shaped form. Along with applying biophysical methods to study the structural folding dynamics of RNAs, I have also developed a technique that improves the production of large quantities of recombinant RNAs in vivo for NMR study. In this project, I report using the wild-type and mutant E.coli strains to produce cost-effective, site-specific labeled, recombinant RNAs. This technique was validated with four representative RNAs of different sizes and complexity to produce milligram amounts of RNAs. The benefit of using site-specific labeled RNAs made from E.coli was demonstrated with several NMR techniques.

Shared active site architecture between archaeal PolD and multi-subunit RNA polymerases revealed by X-ray crystallography

Archaeal replicative DNA polymerase D (PolD) constitute an atypical class of DNA polymerases made of a proofreading exonuclease subunit (DP1) and a larger polymerase catalytic subunit (DP2), both with unknown structures. We have determined the crystal structures of Pyrococcus abyssi DP1 and DP2 at 2.5 and 2.2 angstrom resolution, respectively, revealing a catalytic core strikingly different from all other known DNA polymerases (DNAPs). Rather, the PolD DP2 catalytic core has the same 'double-psi beta-barrel' architecture seen in the RNA polymerase (RNAP) superfamily, which includes multi-subunit transcriptases of all domains of life, homodimeric RNA-silencing pathway RNAPs and atypical viral RNAPs. This finding bridges together, in non-viral world, DNA transcription and DNA replication within the same protein superfamily. This study documents further the complex evolutionary history of the DNA replication apparatus in different domains of life and proposes a classification of all extant DNAPs.

The role of non-coding RNA in the development of pulmonary arterial hypertension

Pulmonary arterial hypertension (PAH) is a progressive disease of the small pulmonary arteries, characterised by pulmonary vascular remodelling due to excessive proliferation and resistance to apoptosis of pulmonary artery endothelial cells (PAECs) and pulmonary artery smooth muscle cells (PASMCs). The increased pulmonary vascular resistance and elevated pulmonary artery pressures result in right heart failure and premature death. Germline mutations of the bone morphogenetic protein receptor-2 (bmpr2) gene, a receptor of the transforming growth factor beta (TGF-β) superfamily, account for approximately 75%-80% of the cases of heritable form of PAH (HPAH) and 20% of sporadic cases or idiopathic PAH (IPAH). IPAH patients without known bmpr2 mutations show reduced expression of BMPR2. However only ~ 20% of bmpr2-mutation carriers will develop the disease, due to an incomplete penetrance, thus the need for a ‘second hit’ including other genetic and/or environmental factors is accepted. Diagnosis of PAH occurs most frequently when patients have reached an advanced stage of disease. Although modern PAH therapies can markedly improve a patient’s symptoms and slow the rate of clinical deterioration, the mortality rate from PAH remains unacceptably high. Therefore, the development of novel therapeutic approaches is required for the treatment of this multifaceted disease. Noncoding RNAs (ncRNAs) include microRNAs (miRNAs) and long noncoding RNAs (lncRNAs). MiRNAs are ~ 22 nucleotide long and act as negative regulators of gene ex-pression via degradation or translational inhibition of their target mRNAs. Previous studies showed extensive evidence for the role of miRNAs in the development of PAH. LncRNAs are transcribed RNA molecules greater than 200 nucleotides in length. Similar to classical mRNA, lncRNAs are translated by RNA polymerase II and are generally alternatively spliced and polyadenylated. LncRNAs are highly versatile and function to regulate gene expression by diverse mechanisms. Unlike miRNAs, which exhibit well-defined actions in negatively regulating gene expression via the 3’-UTR of mRNAs, lncRNAs play more diverse and unpredictable regulatory roles. Although a number of lncRNAs have been intensively investigated in the cancer field, studies of the role of lncRNAs in vascular diseases such as PAH are still at a very early stage. The aim of this study was to investigate the involvement of specific ncRNAs in the development of PAH using experimental animal models and cell culture. The first ncRNA we focused on was miR-143, which is up-regulated in the lung and right ventricle tissues of various animal models of PH, as well as in the lungs and PASMCs of PAH patients. We show that genetic ablation of miR-143 is protective against the development of chronic hypoxia induced PH in mice, assessed via measurement of right ventricular systolic pressure (RVSP), right ventricular hypertrophy (RVH) and pulmonary vascular remodelling. We further report that knockdown of miR-143-3p in WT mice via anti-miR-143-3p administration prior to exposure of mice to chronic hypoxia significantly decreases certain indices of PH (RVSP) although no significant changes in RVH and pulmo-nary vascular remodelling were observed. However, a reversal study using antimiR-143-3p treatment to modulate miR-143-3p demonstrated a protective effect on RVSP, RVH, and muscularisation of pulmonary arteries in the mouse chronic hypoxia induced PH model. In vitro experiments showed that miR-143-3p overexpression promotes PASMC migration and inhibits PASMC apoptosis, while knockdown miR-143-3p elicits the opposite effect, with no effects observed on cellular proliferation. Interestingly, miR-143-3p-enriched exosomes derived from PASMCs mediated cell-to-cell communication between PASMCs and PAECs, contributing to the pro-migratory and pro-angiogenic phenotype of PAECs that underlies the pathogenesis of PAH. Previous work has shown that miR-145-5p expression is upregulated in the chronic hypoxia induced mouse model of PH, as well as in PAH patients. Genetic ablation and pharmacological inhibition (subcutaneous injection) of miR-145-5p exert a protective against the de-velopment of PAH. In order to explore the potential for alternative, more lung targeted delivery strategies, miR-145-5p expression was inhibited in WT mice using intranasal-delivered antimiR-145-5p both prior to and post exposure to chronic hypoxia. The decreased expression of miR-145-5p in lung showed no beneficial effect on the development of PH compared with control antimiRNA treated mice exposed to chronic hypoxia. Thus, miR-143-3p modulated both cellular and exosome-mediated responses in pulmonary vascular cells, while the inhibition of miR-143-3p prevented the development of experimental pulmonary hypertension. We focused on two lncRNAs in this project: Myocardin-induced Smooth Muscle Long noncoding RNA, Inducer of Differentiation (MYOSLID) and non-annotated Myolnc16, which were identified from RNA sequencing studies in human coronary artery smooth muscle cells (HCASMCs) that overexpress myocardin. MYOSLID was significantly in-creased in PASMCs from patients with IPAH compared to healthy controls and increased in circulating endothelial progenitor cells (EPCs) from bmpr2 mutant PAH patients. Exposure of PASMCs to hypoxia in vitro led to a significant upregulation in MYOSLID expres-sion. MYOSLID expression was also induced by treatment of PASMC with BMP4, TGF-β and PDGF, which are known to be triggers of PAH in vitro. Small interfering RNA (siR-NA)-mediated knockdown MYOSLID inhibited migration and induced cell apoptosis without affecting cell proliferation and upregulated several genes in the BMP pathway in-cluding bmpr1α, bmpr2, id1, and id3. Modulation of MYOSLID also affected expression of BMPR2 at the protein level. In addition, MYOSLID knockdown affected the BMP-Smad and BMP-non-Smad signalling pathways in PASMCs assessed by phosphorylation of Smad1/5/9 and ERK1/2, respectively. In PAECs, MYOSLID expression was also induced by hypoxia exposure, VEGF and FGF2 treatment. In addition, MYOSLID knockdown sig-nificantly decreased the proliferation of PAECs. Thus, MYOSLID may be a novel modulator in pulmonary vascular cell functions, likely through the BMP-Smad and –non-Smad pathways. Treatment of PASMCs with inflammatory cytokines (IL-1 and TNF-α) significantly in-duced the expression of Myolnc16 at a very early time point. Knockdown of Myolnc16 in vitro decreased the expression of il-6, and upregulated the expression of il-1 and il-8 in PASMCs. Moreover, the expression levels of chemokines (cxcl1, cxcl6 and cxcl8) were sig-nificantly decreased with Myolnc16 knockdown. In addition, Myolnc16 knockdown decreased the MAP kinase signalling pathway assessed by phosphorylation of ERK1/2 and p38 MAPK and inhibited cell migration and proliferation in PASMCs. Thus, Myolnc16 may a novel modulator of PASMCs functions through anti-inflammatory signalling pathways. In summary, in this thesis we have demonstrated how miR-143-3p plays a protective role in the development of PH both in vivo animal models and patients, as well as in vitro cell cul-ture. Moreover, we have showed the role of two novel lncRNAs in pulmonary vascular cells. These ncRNAs represent potential novel therapeutic targets for the treatment of PAH with further work addressing to investigate the target genes, and the pathways modulated by these ncRNAs during the development of PAH.

Light controlled synthesis of nucleic acids

Dissertação apresentada para obtenção do grau de Doutor em Bioquímica - especialidade Biotecnologia, pela Universidade Nova de Lisboa,Faculdade de Ciências e Tecnologia

Towards modification of Medicago truncatula epigenome: genome editing with engineered nucleases

Periodic drought is the primary limitation of plant growth and crop yield. The rise of water demand caused by the increase in world population and climate change, leads to one of the biggest challenges of modern agriculture: to increase food and feed production. De novo DNA methylation is a process regulated by small interfering RNA (siRNAs), which play a role in plant response and adaptation to abiotic stress. In the particular case of water deficit, growing evidences suggest a link between the siRNA pathways and drought response in the model legume Medicago truncatula. As a first step to understand the role of DNA methylation under water stress, we have set up several bioinformatics and molecular methodologies allowing the design of Clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 systems and the assembly of TALENs (transcription activator-like effector nucleases), to target both dicer-like 3 (MtDCL3) and RNA-Dependent RNA polymerase (MtRDR2), enzymes of the RNA-directed DNA methylation pathway. TALENs efficiency was evaluated prior to plant transformation by a yeast-based assay using two different strategies to test TALENs activity: Polyacrylamide gel electrophoresis (PAGE) and Single strand conformation polymorphisms (SSCP). In this assay, yeast cells triple transformation emerged as good and rapid alternative to laborious yeast mating strategies. PAGE analysis might be a valuable tool to test TALENs efficacy in vivo if we could increase TALENs activity. SSCP-based approach proved to be ineffective due to the generation of several false positives. TALENs and CRISPR/Cas9 system constructed and designed in this work will in the future certainly enable the successful disruption of DCL3 and RDR2 genes and shed the light on the relationship between plant stress resistance and epigenetic regulation mediated by siRNAs in M.truncatula.

Reverse-transcriptase-associated RNaseH activity mediates template switching during reverse transcription in vitro.

During the first steps of reverse transcription of the retroviral genome, sequences present at the extremities of the RNA are used to reconstitute a host cell PolII promoter. The assembly of the promoter occurs by template switching, which takes advantage of a direct repeat at the ends of the RNA molecule. These steps are catalysed by the viral reverse transcriptase, which carries an intrinsic RNaseH activity that is probably also involved therein. To study the role of the RNaseH activity in this first template-switching event, an in vitro system has been developed based on primer extensions of synthetic RNAs. When an RNA was reverse transcribed with wild-type reverse transcriptase in the presence of a second RNA the 3' part of which was repeated at the 5' end of the first one, extension products could be observed corresponding to a chimeric cDNA comprising both RNA species. This template switching could not be detected when a mutant reverse transcriptase lacking the RNaseH activity was used. The results show that the RNaseH activity is needed to remove the 5' RNA sequences from the cDNA:RNA hybrid thereby enabling its translocation to another RNA containing an appropriate complementary target sequence.

Telomerase activity and human telomerase reverse transcriptase mRNA expression in soft tissue tumors: correlation with grade, histology, and proliferative activity.

Telomerase activity (TA) is detected in most human cancers but, with few exceptions, not in normal somatic cells. Little is known about TA in soft tissue tumors. We have examined a series of benign and malignant soft tissue tumors for TA using the telomerase repeat amplification protocol assay. Analysis of the expression of the human telomerase reverse transcriptase was also carried out using RT-PCR. TA was undetectable in benign lesions (15 of 15) and low-grade sarcomas (6 of 6) and was detectable in 50% (19 of 38) of intermediate-/high-grade sarcomas. Although the presence of TA in soft tissue tumors is synonymous with malignancy, it is neither a reliable method in making the distinction between reactive/benign and malignant (especially low-grade) lesions nor a reliable marker of tumor aggressiveness. Leiomyosarcomas and storiform/pleomorphic malignant fibrous histiocytomas rarely showed TA, irrespective of their grade. A strong correlation between human telomerase reverse transcriptase mRNA expression and TA was observed, supporting the close relationship between both parameters. No significant relationship was observed between proliferative activity (as assessed by MIB-1 immunolabeling) and TA. We verified that the absence of telomerase expression was not due to the presence of telomerase inhibitors and therefore alternative mechanism(s) for cell immortalization, yet to be determined, seem to be involved in the development and/or maintenance of some soft tissue sarcomas.

Negative control of quorum sensing by RpoN (sigma54) in Pseudomonas aeruginosa PAO1.

In Pseudomonas aeruginosa PAO1, the expression of several virulence factors such as elastase, rhamnolipids, and hydrogen cyanide depends on quorum-sensing regulation, which involves the lasRI and rhlRI systems controlled by N-(3-oxododecanoyl)-L-homoserine lactone and N-butyryl-L-homoserine lactone, respectively, as signal molecules. In rpoN mutants lacking the transcription factor sigma(54), the expression of the lasR and lasI genes was elevated at low cell densities, whereas expression of the rhlR and rhlI genes was markedly enhanced throughout growth by comparison with the wild type and the complemented mutant strains. As a consequence, the rpoN mutants had elevated levels of both signal molecules and overexpressed the biosynthetic genes for elastase, rhamnolipids, and hydrogen cyanide. The quorum-sensing regulatory protein QscR was not involved in the negative control exerted by RpoN. By contrast, in an rpoN mutant, the expression of the gacA global regulatory gene was significantly increased during the entire growth cycle, whereas another global regulatory gene, vfr, was downregulated at high cell densities. In conclusion, it appears that GacA levels play an important role, probably indirectly, in the RpoN-dependent modulation of the quorum-sensing machinery of P. aeruginosa.