Computational analysis of expressed sequence tags for understanding gene regulation.

High-throughput sequencing has provided a myriad of genetic data for thousands of organisms. Computational analysis of one data type, expressed sequence tags (ESTs) yields insight into gene expression, alternative splicing, tissue specificity gene functionality and the detection and differentiation of pseudogenes. Two computational methods have been developed to analyze alternative splicing events and to detect and characterize pseudogenes using ESTs. A case study of rat phosphodiesterase 4 (PDE4) genes yielded more than twenty-five previously unreported isoforms. These were experimentally verified through wet lab collaboration and found to be tissue specific. In addition, thirteen cytochrome-like gene and pseudogene sequences from the human genome were analyzed for pseudogene properties. Of the thirteen sequences, one was identified as the actual cytochrome gene, two were found to be non-cytochrome-related sequences, and eight were determined to be pseudogenes. The remaining two sequences were identified to be duplicates. As a precursor to applying the two new methods, the efficiency of three BLAST algorithms (NCBI BLAST, WU BLAST and mpiBLAST) were examined for comparing large numbers of short sequences (ESTs) to fewer large sequences (genomic regions). In general, WU BLAST was found to be the most efficient sequence comparison tool. These approaches illustrate the power of ESTs in understanding gene expression. Efficient computational analysis of ESTs (such as the two tools described) will be vital to understanding the complexity of gene expression as more high-throughput EST data is made available via advances in molecular sequencing technologies, such as the current next-generation approaches

HIV-1 Vpr N-terminal tagging affects alternative splicing of the viral genome

To facilitate studies on Vpr function in replicating HIV-1, we aimed to tag the protein in an infectious virus. First we showed that N-, but not C-terminal HA/FLAG tagging of Vpr protein preserves Vpr cytopathicity. Cloning the tags into proviral DNA however ablated viral production and replication. By construction of additional viral variants we could show this defect was not protein-but RNA-dependent and sequence specific, and characterized by oversplicing of the genomic RNA. Simulation of genomic RNA folding suggested that introduction of the tag sequence induced an alternative folding structure in a region enriched in splice sites and splicing regulatory sequences. In silico predictions identified the HA/His(6)-Vpr tagging in HIV-1 to affect mRNA folding less than HA/FLAG-Vpr tagging. In vitro infectivity and mRNA splice pattern improved but did not reach wild-type values. Thus, sequence-specific insertions may interfere with mRNA splicing, possibly due to altered RNA folding. Our results point to the complexity of viral RNA genome sequence interactions. This should be taken into consideration when designing viral manipulation strategies, for both research as for biological interventions

Analysis of the EIAV Rev-Responsive Element (RRE) Reveals a Conserved RNA Motif Required for High Affinity Rev Binding in Both HIV-1 and EIAV

A cis-acting RNA regulatory element, the Rev-responsive element (RRE), has essential roles in replication of lentiviruses, including human immunodeficiency virus (HIV-1) and equine infection anemia virus (EIAV). The RRE binds the viral trans-acting regulatory protein, Rev, to mediate nucleocytoplasmic transport of incompletely spliced mRNAs encoding viral structural genes and genomic RNA. Because of its potential as a clinical target, RRE-Rev interactions have been well studied in HIV-1; however, detailed molecular structures of Rev-RRE complexes in other lentiviruses are still lacking. In this study, we investigate the secondary structure of the EIAV RRE and interrogate regulatory protein-RNA interactions in EIAV Rev-RRE complexes. Computational prediction and detailed chemical probing and footprinting experiments were used to determine the RNA secondary structure of EIAV RRE-1, a 555 nt region that provides RRE function in vivo. Chemical probing experiments confirmed the presence of several predicted loop and stem-loop structures, which are conserved among 140 EIAV sequence variants. Footprinting experiments revealed that Rev binding induces significant structural rearrangement in two conserved domains characterized by stable stem-loop structures. Rev binding region-1 (RBR-1) corresponds to a genetically-defined Rev binding region that overlaps exon 1 of the EIAV rev gene and contains an exonic splicing enhancer (ESE). RBR-2, characterized for the first time in this study, is required for high affinity binding of EIAV Rev to the RRE. RBR-2 contains an RNA structural motif that is also found within the high affinity Rev binding site in HIV-1 (stem-loop IIB), and within or near mapped RRE regions of four additional lentiviruses. The powerful integration of computational and experimental approaches in this study has generated a validated RNA secondary structure for the EIAV RRE and provided provocative evidence that high affinity Rev binding sites of HIV-1 and EIAV share a conserved RNA structural motif. The presence of this motif in phylogenetically divergent lentiviruses suggests that it may play a role in highly conserved interactions that could be targeted in novel anti-lentiviral therapies

Structure and Function of Lentiviral Genomic and Messenger RNA

The positive sense lentiviral RNA genome is packaged within the virus as a dimer of two single strands. The RNA of primate lentiviruses human immunodeficiency virus (HIV-1) and simian immunodeficiency virus (SIVmac239) are distantly related and the secondary structures of these viral RNAs share many known biological functions. Using selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE), I present an analysis of the secondary structure of ex virio genomic SIVmac239 RNA in relation to that of HIV-1 as well as an investigation into the secondary structure of the various in vitro mRNA species of HIV-1 resolved using the SHAPE technique. First, I describe a SHAPE-derived model of SIVmac239 genomic RNA structure. When compared to that of HIV-1, I find very few conserved structural regions outside the previously studied functional structures. I observe that this is due to the flexible nature of the adenosine-rich lentiviral genome. The structures that are conserved are located in regions with high guanosine concentration, forming more stable pairing interactions. These results suggest that lentiviral genomic RNA structure is flexible and metastable unless held by stronger pairing interactions that seem to persist through the course of viral evolution. The lentiviral genomic RNA structures that I have studied do share a few common base pairs, including a small stem-loop at the site of the first splice acceptor (SA1). In the second part, I describe the effect of mutating this structure on viral replication and on the splicing profile of the viral mRNA. To further investigate viral splicing regulation, I determined the SHAPE-derived structures of the most abundant mRNA variants for all of the protein products of HIV-1. Results reveal local interactions that form at regulatory regions in the viral transcripts. Because RNA is an important feature throughout the lentiviral replication cycle, a greater understanding into the role of RNA in various aspects of viral replication will increase comprehension toward the complex biology of infection. This analysis provides insight into evolutionary conservation of RNA structures that play functional roles and may be possible targets for novel factors as part of a broad spectrum of viral inhibitory agents.Doctor of Philosoph

Comparison of SIV and HIV-1 Genomic RNA Structures Reveals Impact of Sequence Evolution on Conserved and Non-Conserved Structural Motifs

RNA secondary structure plays a central role in the replication and metabolism of all RNA viruses, including retroviruses like HIV-1. However, structures with known function represent only a fraction of the secondary structure reported for HIV-1NL4-3. One tool to assess the importance of RNA structures is to examine their conservation over evolutionary time. To this end, we used SHAPE to model the secondary structure of a second primate lentiviral genome, SIVmac239, which shares only 50% sequence identity at the nucleotide level with HIV-1NL4-3. Only about half of the paired nucleotides are paired in both genomic RNAs and, across the genome, just 71 base pairs form with the same pairing partner in both genomes. On average the RNA secondary structure is thus evolving at a much faster rate than the sequence. Structure at the Gag-Pro-Pol frameshift site is maintained but in a significantly altered form, while the impact of selection for maintaining a protein binding interaction can be seen in the conservation of pairing partners in the small RRE stems where Rev binds. Structures that are conserved between SIVmac239 and HIV-1NL4-3 also occur at the 5′ polyadenylation sequence, in the plus strand primer sites, PPT and cPPT, and in the stem-loop structure that includes the first splice acceptor site. The two genomes are adenosine-rich and cytidine-poor. The structured regions are enriched in guanosines, while unpaired regions are enriched in adenosines, and functionaly important structures have stronger base pairing than nonconserved structures. We conclude that much of the secondary structure is the result of fortuitous pairing in a metastable state that reforms during sequence evolution. However, secondary structure elements with important function are stabilized by higher guanosine content that allows regions of structure to persist as sequence evolution proceeds, and, within the confines of selective pressure, allows structures to evolve

Characterization of the longitudinal HIV-1 quasispecies evolution in HIV-1 infected individuals co-infected with Mycobacterium tuberculosis

One of the earliest and most striking observations made about HIV is the extensive genetic variation that the virus has within individual hosts, particularly in the hypervariable regions of the env gene which is divided into 5 variable regions (V1-V5) and 5 more constant (C1-C5) regions. HIV evolves at any time over the course of an individual’s infection and infected individuals harbours a population of genetically related but non-identical viruses that are under constant change and ready to adapt to changes in their environment. These genetically heterogeneous populations of closely related genomes are called quasispecies [65]. Tuberculosis or tubercle forming disease is an acute and/or chronic bacterial infection that primarily attacks the lungs, but which may also affect the kidneys, bones, lymph nodes, and brain. The disease is caused by Mycobacterium tuberculosis (MTB), a slow growing rod-shaped, acid fast bacterium. It is transmitted from person to person through inhalation of bacteria-carrying air droplets. Worldwide, one person out of three is infected with Mycobacterium tuberculosis – two billion people in total. TB currently holds the seventh place in the global ranking of causes of death [73]. In 2008, there were an estimated 9.4 (range, 8.9–9.9 million) million incident cases (equivalent to 139 cases per 100 000 population) of TB globally [75]. A complex biological interplay occurs between M. tuberculosis and HIV in coinfected host that results in the worsening of both pathologies. HIV promotes progression of M. tuberculosis either by endogenous reactivation or exogenous reinfection [77, 78] and, the course of HIV-1 infection is accelerated subsequent to the development of TB [80]. Active TB is associated with an increase in intra-patient HIV-1 diversity both systemically and at the infected lung sites [64,122]. The sustainability or reversal of the HIV-1 quasispecies heterogeneity after TB treatment is not known. Tetanus toxoid vaccinated HIV-1 infected patients developed a transient increase in HIV-1 heterogeneity which was reversed after few weeks [121]. Emergence of a heterogeneous HIV-1 population within a patient may be one of the mechanisms to escape strong immune or drug pressure [65,128]. The existence of better fitting and/or immune escape HIV-variants can lead to an increase in HIV-1 replication [129,130]. It might be that TB favourably selected HIV-1 variants which are sources for consistent HIV-1 replication. Understanding the mechanisms underlying the impacts of TB on HIV-1 is essential for the development of effective measures to reduce TB related morbidity and mortality in HIV-1 infected individuals. In the present study we studied whether the increase in HIV-1 quasispecies diversity during active TB is reversed or preserved throughout the course of antituberculous chemotherapy. For this purpose Two time point HIV-1 quasispecies were evaluated by comparing HIV-1 infected patients with active tuberculosis (HIV-1/TB) and HIV-1 infected patients without tuberculosis (HIV-1/non TB). Plasma samples were obtained from the Frankfurt HIV cohort and HIV-1 RNA was isolated. C2V5 env was amplified by PCR and molecular cloning was performed. Eight to twenty five clones were sequenced from each patient. Various phylogenetic analyses were performed including tree inferences, intra-patient viral diversity and divergence, selective pressure, co-receptor usage prediction and two time point identity of quasispecies comparison using Mantel’s test. We found out from this study that: 1) Active TB sustains HIV-1 quasispecies diversity for longer period 2. Active TB increases the rate of HIV-1 divergence 3) TB might slow down evolution of X4 variants And we concluded that active TB has an impact on HIV-1 viral diversity and divergence over time. The influence of active TB on longitudinal evolution of HIV- 1 may be predominant for R5 viruses. The use of CCR5-coreceptor inhibitors for HIV-1/TB patients as therapeutic approach needs further investigation.Eine der ersten und überraschenden Beobachtungen, welche bei der Analyse des HI-Virus gemacht wurden ist seine ausgeprägte Genetische Variabilität besonders die hypervariable Region des env Genes betreffen. Dieses wird in 5 variable Regionen (V1-V5) sowie 5 stärker konservierte Regionen (C1-C5) unterteilt. HIV wandelt sich zu jedem Zeitpunkt im Verlauf der Infektion und jedes infizierte Individuum ist Träger einer Population von genetisch verwandten jedoch nicht identischen Viren, welche sich kontinuierlich verändern und an die Erfordernisse innerhalb der Umgebung anpassen. Diese genetisch heterogenen, jedoch eng verwandten Populationen werden Quasispecies genannt. Tuberkulose ist eine mykobakterielle Infektion, welche sowohl akute als auch chronische Verläufe zeigt. Neben den Lungen als primärem Manifestationsort können auch die Nieren, Knochen und andere Organe befallen sein. Eine von drei Personen weltweit ist mit Mycobacterium tuberculosis infiziert, insgesamt 2 Milliarden Menschen. In HIV/TB Co-Inifzierten Menschen entsteht ein komplexes Zusammenspiel zwischen HIV und M. tuberculosis, welches zu einer Verschlechterung beider Krankheitsbilder führt. HIV führt durch endogene Rekativierung oder exogene Re-Infektion zu einer Progression der Tuberkulose, welche im weiteren Verlauf die Krankheitsprogression von HIV beschleunigt. Sowohl Morbidität als auch Mortalität sind in HIV-1/TB Co-Infizierten Menschen erhöht. Aktive Lungentuberkulose und Miliartuberkulose gehen mit dem Anstieg der Diversifität der HIV Viren innerhalb eines Wirtes einher. Wie lange diese erhöhte Heterogenität der HIV Quasispecies nach der erfolgreichen Behandlung einer Tuberkulose bestehen bleibt ist bisher noch unklar. Das Verständnis des dem Zusammenspiel von HIV und TB zugrundeliegenden Mechanismus ist essentiell für die Entwicklung von effektiven Massnahmen zur Senkung der Morbidität und Mortalität in HIV/TB Co-infizierten Menschen. Die gegenwärtige Forschungsarbeit folgte daher der Frage, ob wärend einer aktiven TB Infektion eine Zunahme der Diversität der HIV-1 Quasispecies zu beobachten ist und ob diese Diversität während einer TB Therapie erhalten bleibt oder sich zurück bildet. Hierfür wurden die HIV-1 Quasispecies zu zwei Zeitpunkten untersucht, wobei Proben von HIV-1 infizierten Patienten mit aktiver Tuberkulose (HIV-1/TB) und HIV infizierte Patienten ohne Tuberkulose (HIV-1/non TB) verglichen wurden. Aus Plasmaproben der Frankfurter HIV Cohorte wurde HIV-1 RNA isoliert. C2V5 env wurde durch PCR amplifiziert und molekular cloniert. Acht bis fünfundzwanzig Clone wurden für jeden Patienten sequenziert. Mehrere phylogenetische Analysen wurden durchgeführt, welche tree inferences, Intra-Patienten- und virale Diversität und Divergenz, Selektionsdruckanalysen, Vorhersage der Co-Rezeptornutzung sowie Zweipunktanalysen der Identität von Quasispecies mit Hilfe des Mantel’s Test miteinschlossen. Die Analysen ergaben die folgenden Ergebnisse: 1) Eine aktive TB erhält die Diversität von HIV-1 Quasispecies über einen längeren Zeitraum. 2. Eine aktive TB verstärkt die HIV -1 Divergenz 3) TB könnte zu einer langsameren Evolution von X4 Varianten führen. Schlussfolgerung: eine aktive TB beeinflusst die Entwicklung der viralen Diversität und Divergenz von HIV-1 im Verlauf der Krankheit. Der Einfluss der aktiven TB auf die longitudinale Evolution von HIV-1 könnte insbesondere R5 Viren betreffen. Der Einsatz von CCR5-Corezeptor Inhibitoren in HIV-1/TB coinifizerten Patienten sollte daher in Langzeitstudien untersucht werden

Mining Functional Elements in Messenger RNAs: Overview, Challenges, and Perspectives

Eukaryotic messenger RNA (mRNA) contains not only protein-coding regions but also a plethora of functional cis-elements that influence or coordinate a number of regulatory aspects of gene expression, such as mRNA stability, splicing forms, and translation rates. Understanding the rules that apply to each of these element types (e.g., whether the element is defined by primary or higher-order structure) allows for the discovery of novel mechanisms of gene expression as well as the design of transcripts with controlled expression. Bioinformatics plays a major role in creating databases and finding non-evident patterns governing each type of eukaryotic functional element. Much of what we currently know about mRNA regulatory elements in eukaryotes is derived from microorganism and animal systems, with the particularities of plant systems lagging behind. In this review, we provide a general introduction to the most well-known eukaryotic mRNA regulatory motifs (splicing regulatory elements, internal ribosome entry sites, iron-responsive elements, AU-rich elements, zipcodes, and polyadenylation signals) and describe available bioinformatics resources (databases and analysis tools) to analyze eukaryotic transcripts in search of functional elements, focusing on recent trends in bioinformatics methods and tool development. We also discuss future directions in the development of better computational tools based upon current knowledge of these functional elements. Improved computational tools would advance our understanding of the processes underlying gene regulations. We encourage plant bioinformaticians to turn their attention to this subject to help identify novel mechanisms of gene expression regulation using RNA motifs that have potentially evolved or diverged in plant species