Resistance to paclitaxel is associated with a variant of the gene BCL2 in multiple tumor types.

Paclitaxel, the most commonly used form of chemotherapy, is utilized in curative protocols in different types of cancer. The response to treatment differs among patients. Biological interpretation of a mechanism to explain this personalized response is still unavailable. Since paclitaxel is known to target BCL2 and TUBB1, we used pan-cancer genomic data from hundreds of patients to show that a single-nucleotide variant in the BCL2 sequence can predict a patient's response to paclitaxel. Here, we show a connection between this BCL2 genomic variant, its transcript structure, and protein abundance. We demonstrate these findings in silico, in vitro, in formalin-fixed paraffin-embedded (FFPE) tissue, and in patient lymphocytes. We show that tumors with the specific variant are more resistant to paclitaxel. We also show that tumor and normal cells with the variant express higher levels of BCL2 protein, a phenomenon that we validated in an independent cohort of patients. Our results indicate BCL2 sequence variations as determinants of chemotherapy resistance. The knowledge of individual BCL2 genomic sequences prior to the choice of chemotherapy may improve patient survival. The current work also demonstrates the benefit of community-wide, integrative omics data sources combined with in-lab experimentation and validation sets.

Preferential translation of Hsp83 in Leishmania requires a thermosensitive polypyrimidine-rich element in the 3' UTR and involves scanning of the 5' UTR.

Heat shock proteins (HSPs) provide a useful system for studying developmental patterns in the digenetic Leishmania parasites, since their expression is induced in the mammalian life form. Translation regulation plays a key role in control of protein coding genes in trypanosomatids, and is directed exclusively by elements in the 3' untranslated region (UTR). Using sequential deletions of the Leishmania Hsp83 3' UTR (888 nucleotides [nt]), we mapped a region of 150 nt that was required, but not sufficient for preferential translation of a reporter gene at mammalian-like temperatures, suggesting that changes in RNA structure could be involved. An advanced bioinformatics package for prediction of RNA folding (UNAfold) marked the regulatory region on a highly probable structural arm that includes a polypyrimidine tract (PPT). Mutagenesis of this PPT abrogated completely preferential translation of the fused reporter gene. Furthermore, temperature elevation caused the regulatory region to melt more extensively than the same region that lacked the PPT. We propose that at elevated temperatures the regulatory element in the 3' UTR is more accessible to mediators that promote its interaction with the basal translation components at the 5' end during mRNA circularization. Translation initiation of Hsp83 at all temperatures appears to proceed via scanning of the 5' UTR, since a hairpin structure abolishes expression of a fused reporter gene.

In silico whole-genome screening for cancer-related single-nucleotide polymorphisms located in human mRNA untranslated regions.

BACKGROUND: A promising application of the huge amounts of genetic data currently available lies in developing a better understanding of complex diseases, such as cancer. Analysis of publicly available databases can help identify potential candidates for genes or mutations specifically related to the cancer phenotype. In spite of their huge potential to affect gene function, no systematic attention has been paid so far to the changes that occur in untranslated regions of mRNA. RESULTS: In this study, we used Expressed Sequence Tag (EST) databases as a source for cancer-related sequence polymorphism discovery at the whole-genome level. Using a novel computational procedure, we focused on the identification of untranslated region (UTR)-localized non-coding Single Nucleotide Polymorphisms (UTR-SNPs) significantly associated with the tumoral state. To explore possible relationships between genetic mutation and phenotypic variation, bioinformatic tools were used to predict the potential impact of cancer-associated UTR-SNPs on mRNA secondary structure and UTR regulatory elements. We provide a comprehensive and unbiased description of cancer-associated UTR-SNPs that may be useful to define genotypic markers or to propose polymorphisms that can act to alter gene expression levels. Our results suggest that a fraction of cancer-associated UTR-SNPs may have functional consequences on mRNA stability and/or expression. CONCLUSION: We have undertaken a comprehensive effort to identify cancer-associated polymorphisms in untranslated regions of mRNA and to characterize putative functional UTR-SNPs. Alteration of translational control can change the expression of genes in tumor cells, causing an increase or decrease in the concentration of specific proteins. Through the description of testable candidates and the experimental validation of a number of UTR-SNPs discovered on the secreted protein acidic and rich in cysteine (SPARC) gene, this report illustrates the utility of a cross-talk between in silico transcriptomics and cancer genetics.

Primordia vita. Deconvolution from modern sequences.

Evolution of the triplet code is reconstructed on the basis of consensus temporal order of appearance of amino acids. Several important predictions are confirmed by computational sequence analyses. The earliest amino acids, alanine and glycine, have been encoded by GCC and GGC codons, as today. They were succeeded, respectively, by A- and G-series of amino acids, encoded by pyrimidine-central and purine-central codons. The length of the earliest proteins is estimated to be 6-7 residues. The earliest mRNAs were short G+C-rich molecules. These short sequences could have formed hairpins. This is confirmed by analysis of modern prokaryotic mRNA sequences. Predominant size of detected ancient hairpins also corresponds to 6-7 amino acids, as above. Vestiges of last common ancestor can be found in extant proteins in form of entirely conserved short sequences of size six to nine residues present in all or almost all sequenced prokaryotic proteomes (omnipresent motifs). The functions of the topmost conserved octamers are not involved in the basic elementary syntheses. This suggests an initial abiotic supply of amino acids, bases and sugars.

Deleterious mutation prediction in the secondary structure of RNAs.

Methods for computationally predicting deleterious mutations have recently been investigated for proteins, mainly by probabilistic estimations in the context of genomic research for identifying single nucleotide polymorphisms that can potentially affect protein function. It has been demonstrated that in cases where a few homologs are available, ab initio predicted structures modeled by the Rosetta method can become useful for including structural information to improve the deleterious mutation prediction methods for proteins. In the field of RNAs where very few homologs are available at present, this analogy can serve as a precursor to investigate a deleterious mutation prediction approach that is based on RNA secondary structure. When attempting to develop models for the prediction of deleterious mutations in RNAs, useful structural information is available from folding algorithms that predict the secondary structure of RNAs, based on energy minimization. Detecting mutations with desired structural effects among all possible point mutations may then be valuable for the prediction of deleterious mutations that can be tested experimentally. Here, a method is introduced for the prediction of deleterious mutations in the secondary structure of RNAs. The mutation prediction method, based on subdivision of the initial structure into smaller substructures and construction of eigenvalue tables, is independent of the folding algorithms but relies on their success to predict the folding of small RNA structures. Application of this method to predict mutations that may cause structural rearrangements, thereby disrupting stable motifs, is given for prokaryotic transcription termination in the thiamin pyrophosphate and S-adenosyl-methionine induced riboswitches. Ribo switches are mRNA structures that have recently been found to regulate transcription termination or translation initiation in bacteria by conformation rearrangement in response to direct metabolite binding. Predicting deleterious mutations on riboswitches may succeed to systematically intervene in bacterial genetic control.

Sequence-dependent solution structure and motions of 13 TATA/TBP (TATA-box binding protein) complexes.

The TATA element is a well-known example of a DNA promoter sequence recognized by the TATA box binding protein (TBP) through its intrinsic motion and deformability. Although TBP recognizes the TATA element octamer unusually (through the minor groove, which lacks the distinctive features of the major groove), single base-pair replacements alter transcriptional activity. Recent crystallographic experiments have suggested that TATA/TBP complexes differing by a single base pair retain substantial structural similarity despite their functional differences in activating transcription. To investigate the subtle role of sequence-dependent motion within the TATA element and certain aspects of its effect on assembly of the transcriptional complex, we examine 5-ns dynamics trajectories of 13 variant TATA/TBP complexes differing from each other by a single base pair. They include the wild-type (WT) adenovirus 2 major late promoter (AdMLP) TATA element, TATAAAAG (the octamer specifies positions -31 to -24 with respect to the transcription initiation site), and the variants A31 (i.e., AATAAAAG), T30, A29, C29, G28, T28, T27, G26, T26, C25, T25, and T24. Our simulated TATA/TBP complexes develop sequence-dependent structure and motion trends that may lead to favorable orientations for high-activity variants (with respect to binding TFIIA, TFIIB, and other transcription factors), while conversely, accelerate dissociation of low-activity TATA/TBP complexes. The motions that promote favorable geometries for preinitiation complexes include small rotations between TBP's N- and C-terminal domains, sense strand DNA backbone "slithering," and rotations in TBP's H2 and H2' helices. Low-activity variants tend to translate the H1 and H1' helices and withdraw the intercalating phenylalanines. These cumulative DNA and protein motions lead to a spatial spread of complex orientations up to 4 A; this is associated with an overall bend of the variant TATA/TBP complexes that spans 93 degrees to 110 degrees (107 degrees for the crystal reference). Taken together, our analyses imply larger differences when these local structural and bending changes are extended to longer DNA (upstream and downstream) and suggest that specific local TATA/TBP motions (e.g., shifts in TBP helices and TATA bases and backbone) play a role in modulating the formation and maintenance of the transcription initiation complex.