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

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

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.

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.

Aplicações biotecnológicas do supressor de silenciamento de cucumber mosaic virus

Tese de Doutoramento, Ciências Agrárias, Faculdade de Ciências e Tecnologia, Universidade do Algarve, 2015

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.

A coronavirus detected in the vampire bat Desmodus rotundus

This article reports on the identification of a group 2 coronavirus (BatCoV DR/2007) in a Desmodus rotundus vampire bat in Brazil. Phylogenetic analysis of ORF1b revealed that BatCoV DR/2007 originates from a unique lineage in the archetypical group 2 coronaviruses, as described for bat species elsewhere with putative importance in Public Health.

The novel genome organization of the insect picorna-like virus Drosophila C virus suggests this virus belongs to a previously undescribed virus family

The complete nucleotide sequence of the genomic RNA from the insect picorna-like virus Drosophila C virus (DCV) was determined. The DCV sequence predicts a genome organization different to that of other RNA virus families whose sequences are known. The single-stranded positive-sense genomic RNA is 9264 nucleotides in length and contains two large open reading frames (ORFs) which are separated by 191 nucleotides. The 5' ORF contains regions of similarities with the RNA-dependent RNA polymerase, helicase and protease domains of viruses from the picornavirus, comovirus and sequivirus families. The 3' ORF encodes the capsid proteins as confirmed by N-terminal sequence analysis of these proteins. The capsid protein coding region is unusual in two ways: firstly the cistron appears to lack an initiating methionine and secondly no subgenomic RNA is produced, suggesting that the proteins may be translated through internal initiation of translation from the genomic length RNA. The finding of this novel genome organization for DCV shows that this virus is not a member of the Picornaviridae as previously thought, but belongs to a distinct and hitherto unrecognized virus family.

Molecular characterization of the toxic cyanobacterium Cylindrospermopsis raciborskii and design of a species-specific PCR

Cylindrospermopsis raciborskii is a toxic-bloom-forming cyanobacterium that is commonly found in tropical to subtropical climatic regions worldwide, but it is also recognized as a common component of cyanobacterial communities in temperate climates. Genetic profiles of C. raciborskii were examined in 19 cultured isolates originating from geographically diverse regions of Australia and represented by two distinct morphotypes. A 609-bp region of rpoC1, a DNA-dependent RNA polymerase gene, was amplified by PCR from these isolates with cyanobacterium-specific primers. Sequence analysis revealed that all isolates belonged to the same species, including morphotypes with straight or coiled trichomes. Additional rpoC1 gene sequences obtained for a range of cyanobacteria highlighted clustering of C. raciborskii with other heterocyst-producing cyanobacteria (orders Nostocales and Stigonematales). In contrast, randomly amplified polymorphic DNA and short tandemly repeated repetitive sequence profiles revealed a greater level of genetic heterogeneity among C. raciborskii isolates than did rpoC1 gene analysis, and unique band profiles were also found among each of the cyanobacterial genera examined. A PCR test targeting a region of the rpoC1 gene unique to C. raciborskii was developed for the specific identification of C. raciborskii from both purified genomic DNA and environmental samples. The PCR was evaluated with a number of cyanobacterial isolates, but a PCR-positive result was only achieved with C, raciborskii. This method provides an accurate alternative to traditional morphological identification of C. raciborskii.

Indirect immune detection of ecdysone receptor (EcR) during the formation of DNA puffs in Bradysia hygida (Diptera, Sciaridae)

Gene amplification occurs in Bradysia hygida salivary glands, at the end of the fourth larval instar. The hormone 20-hydroxyecdysone (20E) triggers this process, which results in DNA puff formation. Amplified genes are activated in two distinct groups. The activity of the first group is dependent on high levels of 20E, while the second group needs low hormone levels. Consequently, the salivary glands of B. hygida constitute an interesting biological model to study how 20E, and its receptors, affect gene amplification and activity. We produced polyclonal antibodies against B. hygida EcR (BhEcR). In western blots a polypeptide of about 66 kDa was detected in salivary gland extracts. The antibodies were also used for indirect immune-localization of BhEcR in polytene chromosomes. RNA-polymerase II was also immune-detected. We did not detect the receptor in chromosome C where the first and second groups of DNA puffs form during DNA puff anlage formation, but it was present during puff expansion. During the active phase of both groups of DNA puffs, RNA polymerase II co-localized with BhEcR. After puff regression, these antigens were not detected. Apparently, EcR plays a direct role in the transcription of amplified genes, but its role in gene amplification remains enigmatic.

Suppression of Myc oncogenic activity by ribosomal protein haploinsufficiency

The Myc oncogene regulates the expression of several components of the protein synthetic machinery, including ribosomal proteins, initiation factors of translation, RNA polymerase III and ribosomal DNA(1,2). Whether and how increasing the cellular protein synthesis capacity affects the multistep process leading to cancer remains to be addressed. Here we use ribosomal protein heterozygote mice as a genetic tool to restore increased protein synthesis in E mu-Myc/+ transgenic mice to normal levels, and show that the oncogenic potential of Myc in this context is suppressed. Our findings demonstrate that the ability of Myc to increase protein synthesis directly augments cell size and is sufficient to accelerate cell cycle progression independently of known cell cycle targets transcriptionally regulated by Myc. In addition, when protein synthesis is restored to normal levels, Myc- overexpressing precancerous cells are more efficiently eliminated by programmed cell death. Our findings reveal a new mechanism that links increases in general protein synthesis rates downstream of an oncogenic signal to a specific molecular impairment in the modality of translation initiation used to regulate the expression of selective messenger RNAs. We show that an aberrant increase in cap- dependent translation downstream of Myc hyperactivation specifically impairs the translational switch to internal ribosomal entry site ( IRES)- dependent translation that is required for accurate mitotic progression. Failure of this translational switch results in reduced mitotic- specific expression of the endogenous IRES- dependent form of Cdk11 ( also known as Cdc21 and PITSLRE)(3-5), which leads to cytokinesis defects and is associated with increased centrosome numbers and genome instability in E mu-Myc/+ mice. When accurate translational control is re- established in E mu-Myc/+ mice, genome instability is suppressed. Our findings demonstrate how perturbations in translational control provide a highly specific outcome for gene expression, genome stability and cancer initiation that have important implications for understanding the molecular mechanism of cancer formation at the post- genomic level.

High rates of molecular evolution in hantaviruses

Hantaviruses are rodent-borne Bunyaviruses that infect the Arvicolinae, Murinae, and Sigmodontinae subfamilies of Muridae. The rate of molecular evolution in the hantaviruses has been previously estimated at approximately 10(-7) nucleotide substitutions per site, per year (substitutions/site/year), based on the assumption of codivergence and hence shared divergence times with their rodent hosts. If substantiated, this would make the hantaviruses among the slowest evolving of all RNA viruses. However, as hantaviruses replicate with an RNA-dependent RNA polymerase, with error rates in the region of one mutation per genome replication, this low rate of nucleotide substitution is anomalous. Here, we use a Bayesian coalescent approach to estimate the rate of nucleotide substitution from serially sampled gene sequence data for hantaviruses known to infect each of the 3 rodent subfamilies: Araraquara virus ( Sigmodontinae), Dobrava virus ( Murinae), Puumala virus ( Arvicolinae), and Tula virus ( Arvicolinae). Our results reveal that hantaviruses exhibit shortterm substitution rates of 10(-2) to 10(-4) substitutions/site/year and so are within the range exhibited by other RNA viruses. The disparity between this substitution rate and that estimated assuming rodent-hantavirus codivergence suggests that the codivergence hypothesis may need to be reevaluated.

Complementation analysis of the flavivirus Kunjin NS3 and NS5 proteins defines the minimal regions essential for formation of a replication complex and shows a requirement of NS3 in cis for virus assembly

We have previously reported successful trans-complementation of defective Kunjin virus genomic RNAs with a range of large lethal deletions in the nonstructural genes NSI, NS3, and NS5 (A. A. Khromykh et al., J. Virol. 74:3253-3263, 2000). In this study we have mapped further the minimal region in the NS5 gene essential for efficient trans-complementation of genome-length RNAs in repBHK cells to the first 316 of the 905 codons. To allow amplification and easy detection of complemented defective RNAs with deletions apparently affecting virus assembly, we have developed a dual replicon complementation system. In this system defective replicon RNAs with a deletion(s) in the nonstructural genes also encoded the puromycin resistance gene (PAC gene) and the reporter gene for beta-galactosidase (beta-Gal). Complementation of these defective replicon RNAs in repBHK cells resulted in expression of PAC and beta-Gal which allowed establishment of cell lines stably producing replicating defective RNAs by selection with puromycin and comparison of replication efficiencies of complemented defective RNAs by beta-Gal assay. Using this system we demonstrated that deletions in the C-terminal 434 codons of NS3 (codons 178 to 611) were complemented for RNA replication, while any deletions in the first 178 codons were not. None of the genome-length RNAs containing deletions in NS3 shown to be complementable for RNA replication produced secreted defective viruses during complementation in repBHK cells. In contrast, structural proteins produced from these complemented defective RNAs were able to package helper replicon RNA. The results define minimal regions in the NS3 and NS5 genes essential for the formation of complementable replication complex and show a requirement of NS3 in cis for virus assembly.

Human immunodeficiency virus type 1 reverse transcription is stimulated by Tat from other lentiviruses

The tat gene is required by HIV-1 for efficient reverse transcription and this function of Tat can be distinguished from its role in transcription by RNA polymerase II using tat point mutations that abrogate each function independently The mechanism of Tat's role in reverse transcription, however, is not known, nor is it known whether this role is conserved among trans-activating factors in other retroviruses. Here we examine the abilities of heterologous viral trans-activating proteins from jembrana disease virus (jTat), HIV-2 (Tat2), and equine infectious anemia virus (eTat) to substitute for HIV-1 Tat (Tat1) and restore reverse transcription in HIV-1 carrying an inactivated tat gene. Natural endogenous reverse transcription assays showed that trans-activators from some retroviruses (Tat2 and jTat, but not eTat) could substitute for Tat1 in complementation of HIV-1 reverse transcription. Finally, we show that Y47 is critical for Tat1 to function in reverse transcription, but not HIV-1 gene expression. We mutated the homologous position in jTat to H62Y and found it did not improve its ability to stimulate reverse transcription, but an H62A mutation did inhibit jTat complementation. These data highlight the finding that the role of Tat in reverse transcription is not related to trans-activation and demonstrate that other tat genes conserve this function. (C) 2002 Elsevier Science (USA).

The first strand transfer reaction of HIV-1 reverse transcription is more efficient in infected cells than in cell-free natural endogenous reverse transcription reactions

Background: In the presence of dNTPs, intact HIV-1 virions are capable of reverse transcribing at least part of their genome, a process known as natural endogenous reverse transcription (NERT). PCR analysis of virion DNA produced by NERT revealed that the first strand transfer reaction (1stST) was inefficient in intact virions, with minus strand (-) strong stop DNA (ssDNA) copy numbers up to 200 times higher than post-1stST products measured using primers in U3 and U5. This was in marked contrast to the efficiency of 1stST observed in single-round cell infection assays, in which (-) ssDNA and U3-U5 copy numbers were indistinguishable. Objectives: To investigate the reasons for the discrepancy in first strand transfer efficiency between intact cell-free virus and the infection process. Study design: Alterations of both NERT reactions and the conditions of cell infection were used to test whether uncoating and/or entry play a role in the discrepancy in first strand transfer efficiency. Results and Conclusions: The difference in 1stST efficiency could not be attributed simply to viral uncoating, since addition of very low concentrations of detergent to NERT reactions removed the viral envelope without disrupting the reverse transcription complex, and these conditions resulted in no improvement in 1stST efficiency. Virus pseudotyped with surface glycoproteins from either vesicular stomatitis virus or amphotrophic murine leukaemia virus also showed low levels of 1stST in low detergent NERT assays and equivalent levels of (-) ssDNA and 1stST in single-round infections of cells, demonstrating that the gp120-mediated infection process did not select for virions capable of carrying out 1stST. These data indicate that a post-entry event or factor may be involved in efficient HIV-1 reverse transcription in vivo. (C) 2002 Elsevier Science B.V. All rights reserved.

Novos circuitos regulatórios de controlo espacio-temporal do desenvolvimento em bacillus subtilis

RESUMO: A esporulação em Bacillus subtilis é controlada por uma cascata de factores sigma da polimerase do RNA. F e E controlam os estágios precoces do desenvolvimento no pré-esporo e na célula mãe, respectivamente. Numa fase intermédia da diferenciação, quando a célula mãe acaba por envolver o pré-esporo, F é substituído por G e E é substituído por K. Vários mecanismos asseguram que a actividade dos diferentes factores sigma seja confinada a uma janela temporal precisa na célula adequada. Neste estudo, investigámos a função de um factor anti-G, designado por CsfB. Mostramos que para além da sua função de inibição da actividade do factor G em células pré-divisionais, CsfB é também necessário na célula mãe num estágio tardio do desenvolvimento. Mostramos que a expressão de csfB é activada na célula mãe a partir de um promotor dependente de K. Contudo, demonstramos que CsfB interage directamente com E e não com K, e que CsfB é suficiente para inibir a actividade transcricional dependente de E em células vegetativas de B. subtilis. Propomos que CsfB contribui para reduzir o período dependente de E, na linha de expressão genética da célula mãe, desse modo reduzindo a sobreposição entre os regulões E e K e aumentado a fidelidade do processo de desenvolvimento. Uma segunda proteína, YabK, partilha semelhança estrutural com CsfB. YabK é produzida no pré-esporo sob o comando de F, e é necessária para a esporulação. YabK contribui para a transição F/G no programa genético do pré-esporo, porque uma mutação que torna F sensível a CsfB ultrapassa parcialmente a função de YabK na esporulação. No entanto, YabK e CsfB funcionam por mecanismos diferentes, uma vez que YabK não liga directamente a F.---------ABSTRACT: Gene expression during spore development in Bacillus subtilis is governed by a cascade of RNA polymerase sigma factors. F and E control the early stages of development in the forespore and in the mother cell, respectively. At an intermediate stage in the differentiation process, when the larger mother cell finishes engulfment of the smaller forespore, F is replaced by G and E is replaced by K. Several mechanisms ensure the proper timing of activation of the cell type-specific sigma factors. Here, we have investigated the funtion of an anti-sigma G factor, called CsfB. We show here that in addition to its role in inhibiting G in pre-divisional cells, CsfB is also required in the mother cell at a late stage in development. We show that the expression of csfB is activated in the mother cell from a K-specific promoter. However, we demonstrate that CsfB binds directly to E but not to K in a yeast two-hybrid assay, and that CsfB is sufficient to inhibit E-dependent transcriptional activity in vegetative cells of B. subtilis. We posit that CsfB contributes to shutting off the early, E-controlled period in the mother cell line of gene expression, thus reducing the overlap between deployment of the E and K regulons and increasing the fidelity of the developmental process. A second protein, YabK, shares structural similarity with CsfB. YabK is produced in the forespore under F control, and is required for efficient sporulation. YabK contributes to the transition from the F- to the G-dependent period of gene expression, because a mutation that renders F sensitive to CsfB partially bypasses the need for YabK. Yet, YabK and CsfB must function in the control of sigma factor activity by different mechanisms because YabK does not bind directly to F.