UAG readthrough in mammalian cells: Effect of upstream and downstream stop codon contexts reveal different signals

Background  

Translation termination is mediated through an interaction between the release factors eRF1 and eRF3 and the stop codon within its nucleotide context. Although it is well known that the nucleotide contexts both upstream and downstream of the stop codon, can modulate readthrough, little is known about the mechanisms involved.     

Analysis of Low Frequency Protein Truncating Stop-Codon Variants and Fasting Concentration of Growth Hormone

Background

The genetic background of Growth Hormone (GH) secretion is not well understood. Mutations giving rise to a stop codon have a high likelihood of affecting protein function.

Objectives

To analyze likely functional stop codon mutations that are associated with fasting plasma concentration of Growth Hormone.

Methods

We analyzed stop codon mutations in 5451 individuals in the Malmö Diet and Cancer study by genotyping the Illumina Exome Chip. To enrich for stop codon mutations with likely functional effects on protein function, we focused on those disrupting >80% of the predicted amino acid sequence, which were carried by ≥10 individuals. Such mutations were related to GH concentration, measured with a high sensitivity assay (hs-GH) and, if nominally significant, to GH related phenotypes, using linear regression analysis.

Results

Two stop codon mutations were associated with the fasting concentration of hs-GH. rs121909305 (NP_005370.1:p.R93*) [Minor Allele Frequency (MAF) = 0.8%] in the Myosin 1A gene (MYO1A) was associated with a 0.36 (95%CI, 0.04 to 0.54; p=0.02) increment of the standardized value of the natural logarithm of hs-GH per 1 minor allele and rs35699176 (NP_067040.1:p.Q100*) in the Zink Finger protein 77 gene (ZNF77) (MAF = 4.8%) was associated with a 0.12 (95%CI, 0.02 to 0.22; p = 0.02) increase of hs-GH. The mutated high hs-GH associated allele of MYO1A was related to lower BMI (β-coefficient, -0.22; p = 0.05), waist (β-coefficient, -0.22; p = 0.04), body fat percentage (β-coefficient, -0.23; p = 0.03) and with higher HDL (β-coefficient, 0.23; p = 0.04). The ZNF77 stop codon was associated with height (β-coefficient, 0.11; p = 0.02) but not with cardiometabolic risk factors.

Conclusion

We here suggest that a stop codon of MYO1A, disrupting 91% of the predicted amino acid sequence, is associated with higher hs-GH and GH-related traits suggesting that MYO1A is involved in GH metabolism and possibly body fat distribution. However, our results are preliminary and need replication in independent populations.     

Control of <Emphasis Type="Italic">gag-pol</Emphasis> gene expression in the <Emphasis Type="Italic">Candida albicans</Emphasis> retrotransposon Tca2

Background  

In the C. albicans retrotransposon Tca2, the gag and pol ORFs are separated by a UGA stop codon, 3' of which is a potential RNA pseudoknot. It is unclear how the Tca2 gag UGA codon is bypassed to allow pol expression. However, in other retroelements, translational readthrough of the gag stop codon can be directed by its flanking sequence, including a 3' pseudoknot.     

Versatile Dual Reporter Gene Systems for Investigating Stop Codon Readthrough in Plants

Background

Translation is most often terminated when a ribosome encounters the first in-frame stop codon (UAA, UAG or UGA) in an mRNA. However, many viruses (and some cellular mRNAs) contain “stop” codons that cause a proportion of ribosomes to terminate and others to incorporate an amino acid and continue to synthesize a “readthrough”, or C-terminally extended, protein. This dynamic redefinition of codon meaning is dependent on specific sequence context.

Methodology

We describe two versatile dual reporter systems which facilitate investigation of stop codon readthrough in vivo in intact plants, and identification of the amino acid incorporated at the decoded stop codon. The first is based on the reporter enzymes NAN and GUS for which sensitive fluorogenic and histochemical substrates are available; the second on GST and GFP.

Conclusions

We show that the NAN-GUS system can be used for direct in planta measurements of readthrough efficiency following transient expression of reporter constructs in leaves, and moreover, that the system is sufficiently sensitive to permit measurement of readthrough in stably transformed plants. We further show that the GST-GFP system can be used to affinity purify readthrough products for mass spectrometric analysis and provide the first definitive evidence that tyrosine alone is specified in vivo by a ‘leaky’ UAG codon, and tyrosine and tryptophan, respectively, at decoded UAA, and UGA codons in the Tobacco mosaic virus (TMV) readthrough context.     

The origin of a novel gene through overprinting in <Emphasis Type="Italic">Escherichia coli</Emphasis>

Background  

Overlapped genes originate by a) loss of a stop codon among contiguous genes coded in different frames; b) shift to an upstream initiation codon of one of the contiguous genes; or c) by overprinting, whereby a novel open reading frame originates through point mutation inside an existing gene. Although overlapped genes are common in viruses, it is not clear whether overprinting has led to new genes in prokaryotes.     

Comparative Analysis of Codon Usage Bias and Codon Context Patterns between Dipteran and Hymenopteran Sequenced Genomes

Background

Codon bias is a phenomenon of non-uniform usage of codons whereas codon context generally refers to sequential pair of codons in a gene. Although genome sequencing of multiple species of dipteran and hymenopteran insects have been completed only a few of these species have been analyzed for codon usage bias.

Methods and Principal Findings

Here, we use bioinformatics approaches to analyze codon usage bias and codon context patterns in a genome-wide manner among 15 dipteran and 7 hymenopteran insect species. Results show that GAA is the most frequent codon in the dipteran species whereas GAG is the most frequent codon in the hymenopteran species. Data reveals that codons ending with C or G are frequently used in the dipteran genomes whereas codons ending with A or T are frequently used in the hymenopteran genomes. Synonymous codon usage orders (SCUO) vary within genomes in a pattern that seems to be distinct for each species. Based on comparison of 30 one-to-one orthologous genes among 17 species, the fruit fly Drosophila willistoni shows the least codon usage bias whereas the honey bee (Apis mellifera) shows the highest bias. Analysis of codon context patterns of these insects shows that specific codons are frequently used as the 3′- and 5′-context of start and stop codons, respectively.

Conclusions

Codon bias pattern is distinct between dipteran and hymenopteran insects. While codon bias is favored by high GC content of dipteran genomes, high AT content of genes favors biased usage of synonymous codons in the hymenopteran insects. Also, codon context patterns vary among these species largely according to their phylogeny.     

Construction and eukaryotic expression of recombinant large hepatitis delta antigen

Background:

Hepatitis delta virus (HDV) is a subviral human pathogen that exploits host RNA editing activity to produce two essential forms of the sole viral protein, hepatitis delta antigen (HDAg). Editing at the amber/W site of HDV antigenomic RNA leads to the production of the large form (L-HDAg), which is required for RNA packaging.

Methods:

In this study, PCR-based site-directed mutagenesis by the overlap extension method was used to create the point mutation converting the small-HDAg (S-HDAg) stop codon to a tryptophan codon through three stages.

Results:

Sequencing confirmed the desirable mutation and integrity of the L-HDAg open reading frame. The amplicon was ligated into pcDNA3.1 and transfected to Huh7 and HEK 293 cell lines. Western blot analysis using enhanced chemiluminescence confirmed L-HDAg expression. The recombinant L-HDAg localized within the nuclei of cells as determined by immunofluorescence and confocal microscopy.

Conclusion:

Because L-HDAg requires extensive post-translational modifications, the recombinant protein expressed in a mammalian system might be fully functional and applicable as a tool in HDV molecular studies, as well as in future vaccine research.Key Words: Hepatitis Delta Virus, L-HDAg, SOEing-PCR     

Translation regulation of mammalian selenoproteins

Background

Interest in selenium research has considerably grown over the last decades owing to the association of selenium deficiencies with an increased risk of several human diseases, including cancers, cardiovascular disorders and infectious diseases. The discovery of a genetically encoded 21^st amino acid, selenocysteine, is a fascinating breakthrough in molecular biology as it is the first addition to the genetic code deciphered in the 1960s. Selenocysteine is a structural and functional analog of cysteine, where selenium replaces sulfur, and its presence is critical for the catalytic activity of selenoproteins.

Molecular mechanism of selenoprotein P synthesis

Background

Selenoprotein synthesis requires the reinterpretation of a UGA stop codon as one that encodes selenocysteine (Sec), a process that requires a set of dedicated translation factors. Among the mammalian selenoproteins, Selenoprotein P (SELENOP) is unique as it contains a selenocysteine-rich domain that requires multiple Sec incorporation events.

Detecting overlapping coding sequences in virus genomes

Background  

Detecting new coding sequences (CDSs) in viral genomes can be difficult for several reasons. The typically compact genomes often contain a number of overlapping coding and non-coding functional elements, which can result in unusual patterns of codon usage; conservation between related sequences can be difficult to interpret – especially within overlapping genes; and viruses often employ non-canonical translational mechanisms – e.g. frameshifting, stop codon read-through, leaky-scanning and internal ribosome entry sites – which can conceal potentially coding open reading frames (ORFs).     

Gene discovery and functional marker development for fragrance in sorghum (Sorghum bicolor (L.) Moench)

Key message

Sequence analysis and genetic mapping revealed that a 1,444 bp deletion causes a premature stop codon in SbBADH2 of sorghum IS19912. The non-function of SbBADH2 is responsible for fragrance in sorghum IS19912.

Abstract

2-acetyl-1-pyrroline (2AP) is a potent volatile compound causing fragrance in several plants and foods. Seeds of some varieties of rice, sorghum and soybean possess fragrance. The genes responsible for fragrance in rice and soybean are orthologs that correspond to betaine aldehyde dehydrogenase 2 (BADH2). Genotypes harboring fragrance in rice and soybean contain a premature stop codon in BADH2 which impairs the synthesis of full length functional BADH2 protein leading to the accumulation of 2AP. In this study, we reported an association between the BADH2 gene and fragrance in sorghum. An F₂ population of 187 plants developed from a cross between KU630 (non-fragrant) and IS19912 (fragrant) was used. Leaves of F₂ and F₃ progenies were evaluated for fragrance by organoleptic test, while seeds of F₂ plants were analyzed for 2AP. The tests consistently showed that the fragrance is controlled by a single recessive gene. Gene expression analysis of SbBADH1 and SbBADH2 in leaves of KU630 and IS19912 at various stages revealed that SbBADH1 and SbBADH2 were expressed in both accessions. Sequence comparison between KU630 and IS19912 revealed a continuous 1,444 bp deletion encompassing exon 12 to 15 of SbBADH2 in IS19912 which introduces a frameshift mutation and thus causes a premature stop codon. An indel marker was developed to detect polymorphism in SbBADH2. Bulk segregant and QTL analyses confirmed the association between SbBADH2 and fragrance.     

BglBricks: A flexible standard for biological part assembly

Background  

Standard biological parts, such as BioBricks™ parts, provide the foundation for a new engineering discipline that enables the design and construction of synthetic biological systems with a variety of applications in bioenergy, new materials, therapeutics, and environmental remediation. Although the original BioBricks™ assembly standard has found widespread use, it has several shortcomings that limit its range of potential applications. In particular, the system is not suitable for the construction of protein fusions due to an unfavorable scar sequence that encodes an in-frame stop codon.     

Autosomal recessive cerebellar ataxia caused by mutations in the PEX2 gene

Objective

To expand the spectrum of genetic causes of autosomal recessive cerebellar ataxia (ARCA).

Case report

Two brothers are described who developed progressive cerebellar ataxia at 3 1/2 and 18 years, respectively. After ruling out known common genetic causes of ARCA, analysis of blood peroxisomal markers strongly suggested a peroxisomal biogenesis disorder. Sequencing of candidate PEX genes revealed a homozygous c.865_866insA mutation in the PEX2 gene leading to a frameshift 17 codons upstream of the stop codon. PEX gene mutations usually result in a severe neurological phenotype (Zellweger spectrum disorders).

Conclusions

Genetic screening of PEX2 and other PEX genes involved in peroxisomal biogenesis is warranted in children and adults with ARCA.     

Genetic correction of splice site mutation in purified and enriched myoblasts isolated from mdx5cv mice

Background  

Duchenne Muscular Dystrophy (DMD) is an X-linked genetic disorder that results in the production of a dysfunctional form of the protein, dystrophin. The mdx^5cv mouse is a model of DMD in which a point mutation in exon 10 of the dystrophin gene creates an artificial splice site. As a result, a 53 base pair deletion of exon 10 occurs with a coincident creation of a frameshift and a premature stop codon. Using primary myoblasts from mdx^5cv mice, single-stranded DNA oligonucleotides were designed to correct this DNA mutation.     

Evolutionary history of selenocysteine incorporation from the perspective of SECIS binding proteins

Background  

The co-translational incorporation of selenocysteine into nascent polypeptides by recoding the UGA stop codon occurs in all domains of life. In eukaryotes, this event requires at least three specific factors: SECIS binding protein 2 (SBP2), a specific translation elongation factor (eEFSec), selenocysteinyl tRNA, and a cis -acting selenocysteine insertion sequence (SECIS) element in selenoprotein mRNAs. While the phylogenetic relationships of selenoprotein families and the evolution of selenocysteine usage are well documented, the evolutionary history of SECIS binding proteins has not been explored.     

Whole Exome Sequencing Identifies Recessive PKHD1 Mutations in a Chinese Twin Family with Caroli Disease

Background

Mutations in PKHD1 cause autosomal recessive Caroli disease, which is a rare congenital disorder involving cystic dilatation of the intrahepatic bile ducts. However, the mutational spectrum of PKHD1 and the phenotype-genotype correlations have not yet been fully established.

Methods

Whole exome sequencing (WES) was performed on one twin sample with Caroli disease from a Chinese family from Shandong province. Routine Sanger sequencing was used to validate the WES and to carry out segregation studies. We also described the PKHD1 mutation associated with the genotype-phenotype of this twin.

Results

A combination of WES and Sanger sequencing revealed the genetic defect to be a novel compound heterozygous genotype in PKHD1, including the missense mutation c.2507 T>C, predicted to cause a valine to alanine substitution at codon 836 (c.2507T>C, p.Val836Ala), and the nonsense mutation c.2341C>T, which is predicted to result in an arginine to stop codon at codon 781 (c.2341C>T, p.Arg781*). This compound heterozygous genotype co-segregates with the Caroli disease-affected pedigree members, but is absent in 200 normal chromosomes.

Conclusions

Our findings indicate exome sequencing can be useful in the diagnosis of Caroli disease patients and associate a compound heterozygous genotype in PKHD1 with Caroli disease, which further increases our understanding of the mutation spectrum of PKHD1 in association with Caroli disease.     

Stop codons in bacteria are not selectively equivalent

ABSTRACT: BACKGROUND: The evolution and genomic stop codon frequencies have not been rigorously studied with the exception of coding of non-canonical amino acids. Here we study the rate of evolution and frequency distribution of stop codons in bacterial genomes. RESULTS: We show that in bacteria stop codons evolve slower than synonymous sites, suggesting the action of weak negative selection. However, the frequency of stop codons relative to genomic nucleotide content indicated that this selection regime is not straightforward. The frequency of TAA and TGA stop codons is GC-content dependent, with TAA decreasing and TGA increasing with GC-content, while TAG frequency is independent of GC-content. Applying a formal, analytical model to these data we found that the relationship between stop codon frequencies and nucleotide content cannot be explained by mutational biases or selection on nucleotide content. However, with weak nucleotide content-dependent selection on TAG, -0.5 < Nes < 1.5, the model fits all of the data and recapitulates the relationship between TAG and nucleotide content. For biologically plausible rates of mutations we show that, in bacteria, TAG stop codon is universally associated with lower fitness, with TAA being the optimal for G-content < 16% while for G-content > 16% TGA has a higher fitness than TAG. CONCLUSIONS: Our data indicate that TAG codon is universally suboptimal in the bacterial lineage, such that TAA is likely to be the preferred stop codon for low GC content while the TGA is the preferred stop codon for high GC content. The optimization of stop codon usage may therefore be useful in genome engineering or gene expression optimization applications.ReviewersThis article was reviewed by Michail Gelfand, Arcady Mushegian and Shamil Sunyaev. For the full reviews, please go to the Reviewers' Comments section.     

Study of the HIV-2 Env Cytoplasmic Tail Variability and Its Impact on Tat,Rev and Nef

Background

The HIV-2 env’s 3’ end encodes the cytoplasmic tail (CT) of the Env protein. This genomic region also encodes the rev, Tat and Nef protein in overlapping reading frames. We studied the variability in the CT coding region in 46 clinical specimens and in 2 reference strains by sequencing and by culturing. The aims were to analyse the variability of Env CT and the evolution of proteins expressed from overlapping coding sequences.

Results

A 70% reduction of the length of the CT region affected the HIV-2 ROD and EHO strains in vitro due to a premature stop codon in the env gene. In clinical samples this wasn’t observed, but the CT length varied due to insertions and deletions. We noted 3 conserved and 3 variable regions in the CT. The conserved regions were those containing residues involved in Env endocytosis, the potential HIV-2 CT region implicated in the NF-kB activation and the potential end of the lentiviral lytic peptide one. The variable regions were the potential HIV-2 Kennedy region, the potential lentiviral lytic peptide two and the beginning of the potential lentiviral lytic peptide one. A very hydrophobic region was coded downstream of the premature stop codon observed in vitro, suggesting a membrane spanning region. Interestingly, the nucleotides that are responsible for the variability of the CT don’t impact rev and Nef. However, in the Kennedy-like coding region variability resulted only from nucleotide changes that impacted Env and Tat together.

Conclusion

The HIV-2 Env, Tat and Rev C-terminal part are subject to major length variations in both clinical samples and cultured strains. The HIV-2 Env CT contains variable and conserved regions. These regions don’t affect the rev and Nef amino acids composition which evolves independently. In contrast, Tat co-evolves with the Env CT.     

Live virus-free or die: coupling of antivirus immunity and programmed suicide or dormancy in prokaryotes

Background

mtRF1 is a vertebrate mitochondrial protein with an unknown function that arose from a duplication of the mitochondrial release factor mtRF1a. To elucidate the function of mtRF1, we determined the positions that are conserved among mtRF1 sequences but that are different in their mtRF1a paralogs. We subsequently modeled the 3D structure of mtRF1a and mtRF1 bound to the ribosome, highlighting the structural implications of these differences to derive a hypothesis for the function of mtRF1.

Results

Our model predicts, in agreement with the experimental data, that the 3D structure of mtRF1a allows it to recognize the stop codons UAA and UAG in the A-site of the ribosome. In contrast, we show that mtRF1 likely can only bind the ribosome when the A-site is devoid of mRNA. Furthermore, while mtRF1a will adopt its catalytic conformation, in which it functions as a peptidyl-tRNA hydrolase in the ribosome, only upon binding of a stop codon in the A-site, mtRF1 appears specifically adapted to assume this extended, peptidyl-tRNA hydrolyzing conformation in the absence of mRNA in the A-site.

Conclusions

We predict that mtRF1 specifically recognizes ribosomes with an empty A-site and is able to function as a peptidyl-tRNA hydrolase in those situations. Stalled ribosomes with empty A-sites that still contain a tRNA bound to a peptide chain can result from the translation of truncated, stop-codon less mRNAs. We hypothesize that mtRF1 recycles such stalled ribosomes, performing a function that is analogous to that of tmRNA in bacteria.

Reviewers

This article was reviewed by Dr. Eugene Koonin, Prof. Knud H. Nierhaus (nominated by Dr. Sarah Teichmann) and Dr. Shamil Sunyaev.