47 research outputs found
Computational analysis of next generation sequencing data : from transcription start sites in bacteria to human non-coding RNAS
No full textThe advent of next generation sequencing (NGS) technologies has revolutionized the field of molecular biology by providing a wealth of sequence data. “Transcriptomics”, which aims to identify and annotate the complete set of RNA molecules transcribed from a genome, is one of the main applications of these high-throughput methods. Special attention has been paid in determining the exact position of the 5’ ends of RNA transcripts, the transcription start sites (TSSs), and subsequently in identifying the regulatory motifs that are ultimately responsible for governing gene expression. Recently, a novel experimental approach termed dRNA-seq has emerged which enables TSS identification in prokaryotic genomes at a genome-wide scale. While the experimental procedure has reached a point of maturity, the computational downstream analysis of dRNA-seq data is still in its infancy. Analysis of dRNA-seq data was previously done manually, a tedious task that is prone to errors and biases. In order to automate this process we developed a computational tool for accurate and systematic analysis of dRNA-seq data to identify the TSSs genome-wide. In particular, we used a Bayesian framework for TSS calling and a Hidden Markov Model to infer the canonical motifs in the promoter regions of TSSs in order to further capture TSSs that show low evidence of expression. In a second contribution, we exploited the power of next generation sequencing to identify and characterize the expression and processing mechanisms of snoRNAs. SnoRNAs are a particular class of non-protein coding RNAs whose main function is post-transcriptional modification of other non-protein coding RNAs. SnoRNAs carry out their function as part of ribonucleoprotein complexes (RNPs). In order to gain insights into these protein-RNA interactions, we used a technique called PAR-CLIP (Photoactivatable-Ribonucleoside-Enhanced Crosslinking and Immunoprecipitation) that allows the identification of protein-RNA contacts at nucleotide resolution. Using PAR-CLIP data, we were able to demonstrate that snoRNAs undergo precise processing and that many loci in the human genome generate snoRNA-like transcripts whose evolutionary conservation and expression are considerably lower than currently catalogued snoRNAs. Finally, we set out to use small RNA-seq data from the ENCODE project to construct a comprehensive catalog of genomic loci that give rise to snoRNAs. In addition we expanded the current catalog of human snoRNAs and studied the plasticity of snoRNA expression across different cell types. Our analysis confirmed prior observations that several snoRNAs show cell type specific expression, mainly in neurons. A more striking observation was that snoRNA expression appears to be strongly dysregulated in cancers which could lead to the identification of novel biomarkers
How Prevalent is Functional Alternative Splicing in the
No full textthis article can be found at doi: 10.1016/j.tig.2003.12.004 Corresponding author: Gil Ast ([email protected]
Assessing the number of ancestral alternatively spliced exons in the human genome
No full textAbstract Background It is estimated that between 35% and 74% of all human genes undergo alternative splicing. However, as a gene that undergoes alternative splicing can have between one and dozens of alternative exons, the number of alternatively spliced genes by itself is not informative enough. An additional parameter, which was not addressed so far, is therefore the number of human exons that undergo alternative splicing. We have previously described an accurate machine-learning method allowing the detection of conserved alternatively spliced exons without using ESTs, which relies on specific features of the exon and its genomic vicinity that distinguish alternatively spliced exons from constitutive ones. Results In this study we use the above-described approach to calculate that 7.2% (± 1.1%) of all human exons that are conserved in mouse are alternatively spliced in both species. Conclusion This number is the first estimation for the extent of ancestral alternatively spliced exons in the human genome.</p
Understanding the functionality of transcript diversity
No full textRecent years have seen a huge increase in the amount of genomic DNA
being sequenced from a wide variety of organisms, giving us an unprecedented
insight into the molecular diversity seen in nature. As a
result a host of methods have been developed, both experimental and
computational, to understand the functional significance of such diversity
and how it relates to organismal and environmental complexity.
In this thesis I use comparative approaches to explore two areas of
molecular biology where there is evidence for large amounts of transcript
diversity. Firstly, I explore the unprecedented view of microbial
sequence diversity offered by metagenomic sequencing projects, using
sequence similarity and adapted genomic context methods to quantify
the amount of functional novelty in these samples. Secondly, I look
at the transcript diversity generated by alternative splicing. I develop
methods to detect and visualise alternative splicing events and apply
these to the detection of conserved alternative splicing events
Alu-containing exons are alternatively spliced
No full textAlu repetitive elements are found in ∼1.4 million copies in the human genome, comprising more than one-tenth of it. Numerous studies describe exonizations of Alu elements, that is, splicing-mediated insertions of parts of Alu sequences into mature mRNAs. To study the connection between the exonization of Alu elements and alternative splicing, we used a database of ESTs and cDNAs aligned to the human genome. We compiled two exon sets, one of 1176 alternatively spliced internal exons, and another of 4151 constitutively spliced internal exons. Sixty one alternatively spliced internal exons (5.2%) had a significant BLAST hit to an Alu sequence, but none of the constitutively spliced internal exons had such a hit. The vast majority (84%) of the Alu-containing exons that appeared within the coding region of mRNAs caused a frame-shift or a premature termination codon. Alu-containing exons were included in transcripts at lower frequencies than alternatively spliced exons that do not contain an Alu sequence. These results indicate that internal exons that contain an Alu sequence are predominantly, if not exclusively, alternatively spliced. Presumably, evolutionary events that cause a constitutive insertion of an Alu sequence into an mRNA are deleterious and selected against. Alu elements are short interspersed elements (SINEs), typically 300 nucleotides long, which account for>10 % of the huma
MicroRNA-gene association as a prognostic biomarker in cancer exposes disease mechanisms.
No full textThe transcriptional networks that regulate gene expression and modifications to this network are at the core of the cancer phenotype. MicroRNAs, a well-studied species of small non-coding RNA molecules, have been shown to have a central role in regulating gene expression as part of this transcriptional network. Further, microRNA deregulation is associated with cancer development and with tumor progression. Glioblastoma Multiform (GBM) is the most common, aggressive and malignant primary tumor of the brain and is associated with one of the worst 5-year survival rates among all human cancers. To study the transcriptional network and its modifications in GBM, we utilized gene expression, microRNA sequencing, whole genome sequencing and clinical data from hundreds of patients from different datasets. Using these data and a novel microRNA-gene association approach we introduce, we have identified unique microRNAs and their associated genes. This unique behavior is composed of the ability of the quantifiable association of the microRNA and the gene expression levels, which we show stratify patients into clinical subgroups of high statistical significance. Importantly, this stratification goes unobserved by other methods and is not affiliated by other subsets or phenotypes within the data. To investigate the robustness of the introduced approach, we demonstrate, in unrelated datasets, robustness of findings. Among the set of identified microRNA-gene associations, we closely study the example of MAF and hsa-miR-330-3p, and show how their co-behavior stratifies patients into prognosis clinical groups and how whole genome sequences tells us more about a specific genomic variation as a possible basis for patient variances. We argue that these identified associations may indicate previously unexplored specific disease control mechanisms and may be used as basis for further study and for possible therapeutic intervention
Lipid bilayer composition influences small multidrug transporters
No full textBackground:
Membrane proteins are influenced by their surrounding lipids. We investigate the effect of bilayer composition on the membrane transport activity of two members of the small multidrug resistance family; the Escherichia coli transporter, EmrE and the Mycobacterium tuberculosis, TBsmr. In particular we address the influence of phosphatidylethanolamine and anionic lipids on the activity of these multidrug transporters. Phosphatidylethanolamine lipids are native to the membranes of both transporters and also alter the lateral pressure profile of a lipid bilayer. Lipid bilayer lateral pressures affect membrane protein insertion, folding and activity and have been shown to influence reconstitution, topology and activity of membrane transport proteins.
Results:
Both EmrE and TBsmr are found to exhibit a similar dependence on lipid composition, with phosphatidylethanolamine increasing methyl viologen transport. Anionic lipids also increase transport for both EmrE and TBsmr, with the proteins showing a preference for their most prevalent native anionic lipid headgroup; phosphatidylglycerol for EmrE and phosphatidylinositol for TBsmr.
Conclusion:
These findings show that the physical state of the membrane modifies drug transport and that substrate translocation is dependent on in vitro lipid composition. Multidrug transport activity seems to respond to alterations in the lateral forces exerted upon the transport proteins by the bilayer
AluGene: a database of Alu elements incorporated within proteincoding genes
No full textAlu elements are short interspersed elements (SINEs) ~300 nucleotides in length. More than 1 million Alus are found in the human genome. Despite their being genetically functionless, recent ®ndings suggest that Alu elements may have a broad evolutionary impact by affecting gene struc-tures, protein sequences, splicing motifs and expression patterns. Because of these effects, com-piling a genomic database of Alu sequences that reside within protein-coding genes seemed a useful enterprise. Presently, such data are limited since the structural and positional information on genes and Alu sequences are scattered throughout incom-patible and unconnected databases. AluGen
Time to Diagnosis and Persistence: The Two Major Determinants of Effective Tuberculosis Control
No full textThe greatest challenge confronting effective tuberculosis (TB) eradication is the time to diagnosis, and duration of treatment of chronically infected individuals which represent a pool of infection. In an attempt to help limit the spread of TB in New Zealand, a fast SNP based diagnostic test was developed, to quickly identify the highly transmissible and virulent endemic Rangipo strain. The role of VapBC toxin-antitoxin systems in M. tuberculosis has been the subject of great interest recently, due to their expanded number in the genome and links with virulence and the regulation of cell growth in response to environmental stress. Their ability to regulate growth under adverse conditions for presumed survival advantages possibly leading to dormancy or persistence, make them ideal candidates for the development of new M. tuberculosis treatments. To establish differential expression of vapC, and therefore identify possible pathways and functions of the VapBC proteins, RT-qPCR was used to assess the expression levels of vapB and vapC in M. smegmatis under conditions of stress. No consistant changes in vapC mRNA levels were observed, resulting in the hypothesis that it is not the transcriptional differences which are important in the regulation of VapC, but post-transcriptional factors. In order to investigate the function(s) of M. tuberculosis VapBCs, these VapBC proteins were expressed and purified in M. smegmatis, and the VapC toxin tested for RNase activity. The purification, expression, RNase testing and bioinformatic analysis of M. tuberculosis VapCs suggested that VapCRv2530c, VapCRv0065 and VapCRv0617 may all target the same recognition sequence, UA*GG. Bioinformatic analysis revealed an abundance of this target sequence in horizontal gene transfer and TA genes, raising the possibility that VapC toxins could be functioning as selfish elements, or initiating transcriptional regulation cascades when a rapid change in the proteomic response and metabolic state of the cell is required. It is intriguing that the three M. tuberculosis VapC proteins tested thus far appear to target the same recognition sequence, possibly suggesting that all 47 VapCs are RNases and are targeting the same sequence. Alternatively; VapCs may belong to sub-groups targeting different sequences, allowing M. tuberculosis to exude both gross and fine metabolic control; or, they may share the same target, but are regulated by different activators triggered in response to different environmental stimuli
