Ribosomal frameshifting and transcriptional slippage: from genetic steganography and cryptography to adventitious use

Genetic decoding is not ‘frozen’ as was earlier thought, but dynamic. One facet of this is frameshifting that often results in synthesis of a C-terminal region encoded by a new frame. Ribosomal frameshifting is utilized for the synthesis of additional products, for regulatory purposes and for translational ‘correction’ of problem or ‘savior’ indels. Utilization for synthesis of additional products occurs prominently in the decoding of mobile chromosomal element and viral genomes. One class of regulatory frameshifting of stable chromosomal genes governs cellular polyamine levels from yeasts to humans. In many cases of productively utilized frameshifting, the proportion of ribosomes that frameshift at a shift-prone site is enhanced by specific nascent peptide or mRNA context features. Such mRNA signals, which can be 5′ or 3′ of the shift site or both, can act by pairing with ribosomal RNA or as stem loops or pseudoknots even with one component being 4 kb 3′ from the shift site. Transcriptional realignment at slippage-prone sequences also generates productively utilized products encoded trans-frame with respect to the genomic sequence. This too can be enhanced by nucleic acid structure. Together with dynamic codon redefinition, frameshifting is one of the forms of recoding that enriches gene expression.

Étude des mécanismes non-conventionnels de traduction chez le virus de l'immunodéficience humaine de type 1 et le virus de l'hépatite C

Thèse numérisée par la Direction des bibliothèques de l'Université de Montréal.

Cryptic splicing events in the iron transporter ABCB7 and other key target genes in SF3B1-mutant myelodysplastic syndromes

The splicing factor SF3B1 is the most frequently mutated gene in myelodysplastic syndromes (MDS), and is strongly associated with the presence of ring sideroblasts (RS). We have performed a systematic analysis of cryptic splicing abnormalities from RNA sequencing data on hematopoietic stem cells (HSCs) of SF3B1-mutant MDS cases with RS. Aberrant splicing events in many downstream target genes were identified and cryptic 3' splice site usage was a frequent event in SF3B1-mutant MDS. The iron transporter ABCB7 is a well-recognized candidate gene showing marked downregulation in MDS with RS. Our analysis unveiled aberrant ABCB7 splicing, due to usage of an alternative 3' splice site in MDS patient samples, giving rise to a premature termination codon in the ABCB7 mRNA. Treatment of cultured SF3B1-mutant MDS erythroblasts and a CRISPR/Cas9-generated SF3B1-mutant cell line with the nonsense-mediated decay (NMD) inhibitor cycloheximide showed that the aberrantly spliced ABCB7 transcript is targeted by NMD. We describe cryptic splicing events in the HSCs of SF3B1-mutant MDS, and our data support a model in which NMD-induced downregulation of the iron exporter ABCB7 mRNA transcript resulting from aberrant splicing caused by mutant SF3B1 underlies the increased mitochondrial iron accumulation found in MDS patients with RS.Leukemia advance online publication, 17 June 2016; doi:10.1038/leu.2016.149.

External quality assessment of BRAF molecular analysis in melanoma

The availability of BRAF inhibitors has given metastatic melanoma patients an effective new treatment choice and molecular testing to determine the presence or absence of a BRAF codon 600 mutation is pivotal in the clinical management of these patients. This molecular test must be performed accurately and appropriately to ensure that the patient receives the most suitable treatment in a timely manner. Laboratories have introduced such testing; however, some experience low sample throughput making it critical that an external quality assurance programme is available to help promote a high standard of testing, reporting and provide an educational aspect for BRAF molecular testing. Laboratories took part in three rounds of external quality assessment (EQA) during a 12-month period giving participants a measure of the accuracy of genotyping, clinical interpretation of the result and experience in testing a range of different samples. Formalin fixed paraffin embedded tissue sections from malignant melanoma patients were distributed to participants for BRAF molecular testing. The standard of testing was generally high but distribution of a mutation other than the most common, p.(Val600Glu), highlighted concerns with detection or reporting of the presence of rarer mutations. The main issues raised in the interpretation of the results were the importance of clear unambiguous interpretation of the result tailored to the patient and the understanding that the treatment is different from that given to other stratified medicine programmes. The variability in reporting and wide range of methodologies used indicate a continuing need for EQA in this field.

Saccharomycotin transcriptomics by next-generation sequencing

The non-standard decoding of the CUG codon in Candida cylindracea raises a number of questions about the evolutionary process of this organism and other species Candida clade for which the codon is ambiguous. In order to find some answers we studied the transcriptome of C. cylindracea, comparing its behavior with that of Saccharomyces cerevisiae (standard decoder) and Candida albicans (ambiguous decoder). The transcriptome characterization was performed using RNA-seq. This approach has several advantages over microarrays and its application is booming. TopHat and Cufflinks were the software used to build the protocol that allowed for gene quantification. About 95% of the reads were mapped on the genome. 3693 genes were analyzed, of which 1338 had a non-standard start codon (TTG/CTG) and the percentage of expressed genes was 99.4%. Most genes have intermediate levels of expression, some have little or no expression and a minority is highly expressed. The distribution profile of the CUG between the three species is different, but it can be significantly associated to gene expression levels: genes with fewer CUGs are the most highly expressed. However, CUG content is not related to the conservation level: more and less conserved genes have, on average, an equal number of CUGs. The most conserved genes are the most expressed. The lipase genes corroborate the results obtained for most genes of C. cylindracea since they are very rich in CUGs and nothing conserved. The reduced amount of CUG codons that was observed in highly expressed genes may be due, possibly, to an insufficient number of tRNA genes to cope with more CUGs without compromising translational efficiency. From the enrichment analysis, it was confirmed that the most conserved genes are associated with basic functions such as translation, pathogenesis and metabolism. From this set, genes with more or less CUGs seem to have different functions. The key issues on the evolutionary phenomenon remain unclear. However, the results are consistent with previous observations and shows a variety of conclusions that in future analyzes should be taken into consideration, since it was the first time that such a study was conducted.

Alternate Conformations Regulate Ribosomal Recoding in a Positive-sense RNA Virus

Positive-sense RNA viruses are important animal, plant, insect and bacteria pathogens and constitute the largest group of RNA viruses. Due to the relatively small size of their genomes, these viruses have evolved a variety of non-canonical translation mechanisms to optimize coding capacity expanding their proteome diversity. One such strategy is codon redefinition or recoding. First described in viruses, recoding is a programmed translation event in which codon alterations are context dependent. Recoding takes place in a subset of messenger RNA (mRNAs) with some products reflecting new, and some reflecting standard, meanings. The ratio between the two is both critical and highly regulated. While a variety of recoding mechanisms have been documented, (ribosome shunting, stop-carry on, termination-reinitiation, and translational bypassing), the two most extensively employed by RNA viruses are Programmed Ribosomal Frameshifting (PRF) and Programmed Ribosomal Readthrough (PRT). While both PRT and PRF subvert normal decoding for expression of C-terminal extension products, the former involves an alteration of reading frame, and the latter requires decoding of a non-sense codon. Both processes occur at a low but defined frequency, and both require Recoding Stimulatory Elements (RSE) for regulation and optimum functionality. These stimulatory signals can be embedded in the RNA in the form of sequence or secondary structure, or trans-acting factors outside the mRNA such as proteins or micro RNAs (miRNA). Despite 40+ years of study, the precise mechanisms by which viral RSE mediate ribosome recoding for the synthesis of their proteins, or how the ratio of these products is maintained, is poorly defined. This study reveals that in addition to a long distance RNA:RNA interaction, three alternate conformations and a phylogenetically conserved pseudoknot regulate PRT in the carmovirus Turnip crinkle virus (TCV).

Development of a novel platform for high-throughput gene design and artificial gene synthesis to produce large libraries of recombinant venom peptides for drug discovery

Tese de Doutoramento em Ciências Veterinárias na Especialidade de Ciências Biológicas e Biomédicas

Targeted gene modification in Fragaria vesca mediated by CRISPR/Cas9 system

Genome editing is becoming an important biotechnological tool for gene function analysis and crop improvement, being the CRISPR-Cas9 (Clustered Regularly Interspaced Short Palindromic Repeat-CRISPR associated protein 9) system the most widely used. The natural CRISPR/Cas9 system has been reduced to two components: a single-guide RNA (sgRNA) for target recognition via RNA-DNA base pairing, which is commonly expressed using a promoter for small-RNAs (U6 promoter), and the Cas9 endonuclease for DNA cleavage (1). To validate the CRISPR/Cas9 system in strawberry plants, we designed two sgRNAs directed against the floral homeotic gene APETALA3 (sgRNA-AP3#1 and sgRNA-AP3#2). This gene was selected because ap3 mutations induce clear developmental phenotypes in which petals and stamens are missing or partially converted to sepals and carpels respectively (2). In this work, we used two different U6 promoters to drive the sgRNA-AP3s expression: AtU6-26 from Arabidopsis (4), and a U6 promoter from Fragaria vesca (FvU6) (this work). We also tested two different coding sequences of Cas9: a human- (hSpCas9) (3) and a plant-codon optimized (pSpCas9) (this work). Transient expression experiments using both CRISPR/Cas9 systems (AtU6-26:sgRNA-AP3#1_35S:hSpCas9_AtU6-26:sgRNA-AP3#2 and FvU6:sgRNA-AP3#1_35S:pSpCas9_FvU6:sgRNA-AP3#2) were performed infiltrating Agrobacterium tumefaciens into F. vesca fruits. PCR amplification and sequencing analyses across the target sites showed a deletion of 188-189 bp corresponding to the region comprised between the two cutting sites of Cas9, confirming that the CRISPR/Cas9 system is functional in F. vesca. Remarkably, the two systems showed different mutagenic efficiency that could be related to differences in expression of the U6 promoters as well as differences in the Cas9 transcripts stability and translation. Stable transformants for both F. vesca (2n) and Fragaria X anannassa (8n) are currently being established to test whether is possible to obtain heritable homozygous mutants derived from CRISPR/Cas9 strategies in strawberry. Thus, our work offers a promising tool for genome editing and gene functional analysis in strawberry. This tool might represent a more efficient alternative to the sometimes inefficient RNAi silencing methods commonly used in this species.

Hormone-Mediated Gene-Specific Translation Regulation

The central role of translation regulation in the control of critical cellular processes has long been recognized. Yet the systematic exploration of quantitative changes in translation at a genome-wide scale in response to specific stimuli has only recently become technically feasible. Using a genetic approach, we have identified new Arabidopsis weak-ethylene insensitive mutants that also display defects in translation, which suggested the existence of a previously unknown molecular module involved in ethylene-mediated translation regulation of components of this signaling pathway. To explore this link in detail, we implemented for Arabidopsis the ribosome-footprinting technology, which enables the study of translation at a whole-genome level at single codon resolution[1]. Using ribosome-footprinting we examined the effects of short exposure to ethylene on the Arabidopsis translatome looking for ethylene-triggered changes in translation rates that could not be explained by changes in transcript levels. The results of this research, in combination with the characterization of a subset of the aforementioned weak-ethylene insensitive mutants that are defective in the UPF genes (core-components of the nonsense-mediated mRNA decay machinery), uncovered a translation-based branch of the ethylene signaling pathway[2]. In the presence of ethylene, translation of a negative regulator of ethylene signaling EBF2 is repressed, despite induced transcription of this gene. These translational effects of ethylene require the long 3´UTR of EBF2 (3´EBF2), which is recognized by the C-terminal end of the key ethylene-signaling protein EIN2 (EIN2C) in the cytoplasm once EIN2C is released from the ER-membrane by proteolytic cleavage. EIN2C binds the 3´EBF2, recruits the UPF proteins and moves to P-bodies, where the translation of EBF2 in inhibited despite its mRNA accumulation. Once the ethylene signal is withdrawn, the translation of the stored EBF2 mRNAs is resumed, thus rapidly dampening the ethylene response. These findings represent a mechanistic paradigm of gene-specific regulation of translation in response to a key growth regulator. Translation regulatory elements can be located in both 3′ and 5′ UTRs. We are now focusing on the ead1 and ead2 mutants, another set of ethylene-signaling mutants defective in translational regulation. Ribosome-footprinting on the ead1 mutant revealed an accumulation of translating ribosomes in the 5´UTRs of uORF-containing genes and reduction in the levels of ribosomes in the main ORF. The mutant is also impaired in the translation of GFP when this reporter is fused to WT 5´UTR of potential EAD1 targets but not when GFP is fused to the uORF-less versions of the same 5´UTRs. Our hypothesis is that EAD1/2 work as a complex that is required for the efficient translation of mRNAs that have common structural (complex 5´UTR with uORFs) and functional (regulation of key cellular processes) features. We are working towards the identification of the conditions where the EAD1 regulation of translation is required. [1] Ingolia, N. et al. (2009) Genome-Wide Analysis in Vivo of Translation with Nucleotide Resolution Using Ribosome Profiling. Science, 324; 218-222 [2] Merchante, C. et al. (2015) Gene-Specific Translation Regulation Mediated by the Hormone-Signaling Molecule EIN2. Cell, 163(3): 684-697

Charcot Marie Tooth disease (CMT4A) due to GDAP1 mutation: report of a colombian family

Background: Mutations of GDAP1 gene cause autosomal dominant and autosomal recessive Charcot-Marie-Tooth disease and more than 40 different mutations have been reported. The recessive Q163X mutation has been described in patients of Spanish ancestry, and a founder mutation in South American patients, originating in Spain has been demonstrated. Objective: We describe physical and histological features, and the molecular impact of mutation Q163X in a Colombian family. Methods: We report two female patients, daughters of consanguineous parents, with onset of symptoms within the first two years of life, developing severe functional impairment, without evidence of dysmorphic features, hoarseness or diaphragmatic paralysis. Electrophysiology tests showed a sensory and motor neuropathy with axonal pattern. Sequencing of GDAP1 gene was requested and the study identified a homozygous point mutation (c.487 C>T) in exon 4, resulting in a premature stop codon (p.Q163X). This result confirms the diagnosis of Charcot-Marie-Tooth disease, type 4A. Results: The patients were referred to Physical Medicine and Rehabilitation service, in order to be evaluated for ambulation assistance. They have been followed by Pulmonology service, for pulmonary function assessment and diaphragmatic paralysis evaluation. Genetic counseling was offered. The study of the genealogy of the patient, phenotypic features, and electrophysiological findings must be included as valuable tools in the clinical approach of the patient with Charcot-Marie-Tooth disease, in order to define a causative mutation. In patients of South American origin, the presence of GDAP1 gene mutations should be considered, especially the Q163X mutation, as the cause of CMT4A disease.

Charcot Marie Tooth disease (CMT4A) due to GDAP1 mutation: report of a colombian family

Background: Mutations of GDAP1 gene cause autosomal dominant and autosomal recessive Charcot-Marie-Tooth disease and more than 40 different mutations have been reported. The recessive Q163X mutation has been described in patients of Spanish ancestry, and a founder mutation in South American patients, originating in Spain has been demonstrated. Objective: we describe physical and histological features, and the molecular impact of mutation Q163X in a Colombian family. Methods: We report two female patients, daughters of consanguineous parents, with onset of symptoms within the first two years of life, developing severe functional impairment, without evidence of dysmorphic features, hoarseness or diaphragmatic paralysis. Electrophysiology tests showed a sensory and motor neuropathy with axonal pattern. Sequencing of GDAP1 gene was requested and the study identified a homozygous point mutation (c.487 C>T) in exon 4, resulting in a premature stop codon (p.Q163X). This result confirms the diagnosis of Charcot-Marie-Tooth disease, type 4A. Results: The patients were referred to Physical Medicine and Rehabilitation service, in order to be evaluated for ambulation assistance. They have been followed by Pulmonology service, for pulmonary function assessment and diaphragmatic paralysis evaluation. Genetic counseling was offered. The study of the genealogy of the patient, phenotypic features, and electrophysiological findings must be included as valuable tools in the clinical approach of the patient with Charcot-Marie-Tooth disease, in order to define a causative mutation. In patients of South American origin, the presence of GDAP1 gene mutations should be considered, especially the Q163X mutation, as the cause of CMT4A disease.

Metabolic consequences of knocking out UGT85B1, the gene encoding the glucosyltransferase required for synthesis of dhurrin in Sorghum bicolor (L. Moench)

Many important food crops produce cyanogenic glucosides as natural defense compounds to protect against herbivory or pathogen attack. It has also been suggested that these nitrogen-based secondary metabolites act as storage reserves of nitrogen. In sorghum, three key genes, CYP79A1, CYP71E1 and UGT85B1, encode two Cytochrome P450s and a glycosyltransferase, respectively, the enzymes essential for synthesis of the cyanogenic glucoside dhurrin. Here, we report the use of targeted induced local lesions in genomes (TILLING) to identify a line with a mutation resulting in a premature stop codon in the N-terminal region of UGT85B1. Plants homozygous for this mutation do not produce dhurrin and are designated tcd2 (totally cyanide deficient 2) mutants. They have reduced vigor, being dwarfed, with poor root development and low fertility. Analysis using liquid chromatography-mass spectrometry (LC-MS) shows that tcd2 mutants accumulate numerous dhurrin pathway-derived metabolites, some of which are similar to those observed in transgenic Arabidopsis expressing the CYP79A1 and CYP71E1 genes. Our results demonstrate that UGT85B1 is essential for formation of dhurrin in sorghum with no co-expressed endogenous UDP-glucosyltransferases able to replace it. The tcd2 mutant suffers from self-intoxication because sorghum does not have a feedback mechanism to inhibit the initial steps of dhurrin biosynthesis when the glucosyltransferase activity required to complete the synthesis of dhurrin is lacking. The LC-MS analyses also revealed the presence of metabolites in the tcd2 mutant which have been suggested to be derived from dhurrin via endogenous pathways for nitrogen recovery, thus indicating which enzymes may be involved in such pathways.

Novel ciliate genetic code variants including the reassignment of all three stop codons to sense codons in Condylostoma magnum

mRNA translation in many ciliates utilizes variant genetic codes where stop codons are reassigned to specify amino acids. To characterize the repertoire of ciliate genetic codes, we analyzed ciliate transcriptomes from marine environments. Using codon substitution frequencies in ciliate protein-coding genes and their orthologs, we inferred the genetic codes of 24 ciliate species. Nine did not match genetic code tables currently assigned by NCBI. Surprisingly, we identified a novel genetic code where all three standard stop codons (TAA, TAG, and TGA) specify amino acids in Condylostoma magnum. We provide evidence suggesting that the functions of these codons in C. magnum depend on their location within mRNA. They are decoded as amino acids at internal positions, but specify translation termination when in close proximity to an mRNA 3' end. The frequency of stop codons in protein coding sequences of closely related Climacostomum virens suggests that it may represent a transitory state.mRNA translation in many ciliates utilizes variant genetic codes where stop codons are reassigned to specify amino acids. To characterize the repertoire of ciliate genetic codes, we analyzed ciliate transcriptomes from marine environments. Using codon substitution frequencies in ciliate protein-coding genes and their orthologs, we inferred the genetic codes of 24 ciliate species. Nine did not match genetic code tables currently assigned by NCBI. Surprisingly, we identified a novel genetic code where all three standard stop codons (TAA, TAG, and TGA) specify amino acids in Condylostoma magnum. We provide evidence suggesting that the functions of these codons in C. magnum depend on their location within mRNA. They are decoded as amino acids at internal positions, but specify translation termination when in close proximity to an mRNA 3' end. The frequency of stop codons in protein coding sequences of closely related Climacostomum virens suggests that it may represent a transitory state.

Comparative survey of the relative impact of mRNA features on local ribosome profiling read density

Ribosome profiling (Ribo-seq), a promising technology for exploring ribosome decoding rates, is characterized by the presence of infrequent high peaks in ribosome footprint density and by long alignment gaps. Here, to reduce the impact of data heterogeneity we introduce a simple normalization method, Ribo-seq Unit Step Transformation (RUST). RUST is robust and outperforms other normalization techniques in the presence of heterogeneous noise. We illustrate how RUST can be used for identifying mRNA sequence features that affect ribosome footprint densities globally. We show that a few parameters extracted with RUST are sufficient for predicting experimental densities with high accuracy. Importantly the application of RUST to 30 publicly available Ribo-seq data sets revealed a substantial variation in sequence determinants of ribosome footprint frequencies, questioning the reliability of Ribo-seq as an accurate representation of local ribosome densities without prior quality control. This emphasizes our incomplete understanding of how protocol parameters affect ribosome footprint densities.

Synonymous co-variation across the E1/E2 gene junction of hepatitis C virus defines virion fitness

Hepatitis C virus is a positive-sense single-stranded RNA virus. The gene junction partitioning the viral glycoproteins E1 and E2 displays concurrent sequence evolution with the 3′-end of E1 highly conserved and the 5′-end of E2 highly heterogeneous. This gene junction is also believed to contain structured RNA elements, with a growing body of evidence suggesting that such structures can act as an additional level of viral replication and transcriptional control. We have previously used ultradeep pyrosequencing to analyze an amplicon library spanning the E1/E2 gene junction from a treatment naïve patient where samples were collected over 10 years of chronic HCV infection. During this timeframe maintenance of an in-frame insertion, recombination and humoral immune targeting of discrete virus sub-populations was reported. In the current study, we present evidence of epistatic evolution across the E1/E2 gene junction and observe the development of co-varying networks of codons set against a background of a complex virome with periodic shifts in population dominance. Overtime, the number of codons actively mutating decreases for all virus groupings. We identify strong synonymous co-variation between codon sites in a group of sequences harbouring a 3 bp in-frame insertion and propose that synonymous mutation acts to stabilize the RNA structural backbone.