Unusual effect of clusters of rare arginine (AGG) codons on the expression of human interferon alpha 1 gene in Escherichia coli

Application of codon usage and context analysis in genes up- or down-regulated in neurodegeneration and cancer to combat comorbidities

Introduction

Neurodegeneration and cancer present in comorbidities with inverse effects due to the expression of genes and pathways acting in opposition. Identifying and studying the genes simultaneously up or downregulated during morbidities helps curb both ailments together.

Methods

This study examines four genes. Three of these (Amyloid Beta Precursor Protein (APP), Cyclin D1 (CCND1), and Cyclin E2 (CCNE2) are upregulated, and one protein phosphatase 2 phosphatase activator (PTPA) is simultaneously downregulated in both disorders. We investigated molecular patterns, codon usage, codon usage bias, nucleotide bias in the third codon position, preferred codons, preferred codon pairs, rare codons, and codon context.

Results

Parity analysis revealed that T is preferred over A, and G is preferred over C in the third codon position, suggesting composition plays no role in nucleotide bias in both the upregulated and downregulated gene sets and that mutational forces are stronger in upregulated gene sets than in downregulated ones. Transcript length influenced the overall %A composition and codon bias, and the codon AGG exerted the strongest influence on codon usage in both the upregulated and downregulated gene sets. Codons ending in G/C were preferred for 16 amino acids, and glutamic acid-, aspartic acid-, leucine-, valine-, and phenylalanine-initiated codon pairs were preferred in all genes. Codons CTA (Leu), GTA (Val), CAA (Gln), and CGT (Arg) were underrepresented in all examined genes.

Discussion

Using advanced gene editing tools such as CRISPR/Cas or any other gene augmentation technique, these recoded genes may be introduced into the human body to optimize gene expression levels to augment neurodegeneration and cancer therapeutic regimens simultaneously.

Accounting for Programmed Ribosomal Frameshifting in the Computation of Codon Usage Bias Indices

Experimental evidence shows that synonymous mutations can have important consequences on genetic fitness. Many organisms display codon usage bias (CUB), where synonymous codons that are translated into the same amino acid appear with distinct frequency. Within genomes, CUB is thought to arise from selection for translational efficiency and accuracy, termed the translational efficiency hypothesis (TEH). Indeed, CUB indices correlate with protein expression levels, which is widely interpreted as evidence for translational selection. However, these tests neglect -1 programmed ribosomal frameshifting (-1 PRF), an important translational disruption effect found across all organisms of the tree of life. Genes that contain -1 PRF signals should cost more to express than genes without. Thus, CUB indices that do not consider -1 PRF may overestimate genes’ true adaptation to translational efficiency and accuracy constraints. Here, we first investigate whether -1 PRF signals do indeed carry such translational cost. We then propose two corrections for CUB indices for genes containing -1 PRF signals. We retest the TEH in Saccharomyces cerevisiae under these corrections. We find that the correlation between corrected CUB index and protein expression remains intact for most levels of uniform -1 PRF efficiencies, and tends to increase when these efficiencies decline with protein expression. We conclude that the TEH is strengthened and that -1 PRF events constitute a promising and useful tool to examine the relationships between CUB and selection for translation efficiency and accuracy.

Asparaginase treatment side-effects may be due to genes with homopolymeric Asn codons (Review-Hypothesis)

The present treatment of childhood T-cell leukemias involves the systemic administration of prokary-otic L-asparaginase (ASNase), which depletes plasma Asparagine (Asn) and inhibits protein synthesis. The mechanism of therapeutic action of ASNase is poorly understood, as are the etiologies of the side-effects incurred by treatment. Protein expression from genes bearing Asn homopolymeric coding regions (N-hCR) may be particularly susceptible to Asn level fluctuation. In mammals, N-hCR are rare, short and conserved. In humans, misfunctions of genes encoding N-hCR are associated with a cluster of disorders that mimic ASNase therapy side-effects which include impaired glycemic control, dislipidemia, pancreatitis, compromised vascular integrity, and neurological dysfunction. This paper proposes that dysregulation of Asn homeostasis, potentially even by ASNase produced by the microbiome, may contribute to several clinically important syndromes by altering expression of N-hCR bearing genes. By altering amino acid abundance and modulating ribosome translocation rates at codon repeats, the microbiomic environment may contribute to genome decoding and to shaping the proteome. We suggest that impaired translation at poly Asn codons elevates diabetes risk and severity.

CBDB: the codon bias database

Background

In many genomes, a clear preference in the usage of particular codons exists. The mechanisms that induce codon biases remain an open question; studies have attributed codon usage to translational selection, mutational bias and drift. Furthermore, correlations between codon usage within host genomes and their viral pathogens have been observed for a myriad of host-virus systems. As such, numerous studies have investigated codon usage and codon bias in an effort to better understand how species evolve. Numerous metrics have been developed to identify biases in codon usage. In addition, a few data repositories of codon bias data are available, differing in the metrics reported as well as the number and taxonomy of strains examined.

Description

We have created a new web resource called the Codon Bias Database (CBDB) which provides information regarding the codon bias within the set of highly expressed genes for 300+ bacterial genomes. CBDB was developed to provide a resource for researchers investigating codon bias in bacteria, facilitating comparisons between strains and species. Furthermore, the site was created to serve those studying adaptation in phage; the genera selected for this first release of CBDB all have sequenced, annotated bacteriophages. The annotations and sequences for the highly expressed gene set are available for each strain in addition to the strain’s codon bias measurements.

Conclusions

Comparing species and strains provides a comprehensive look at how codon usage has been shaped over evolutionary time and can elucidate the putative mechanisms behind it. The Codon Bias Database provides a centralized repository of look-up tables and codon usage bias measures for a wide variety of genera, species and strains. Through our analysis of the variation in codon usage within the strains presently available, we find that most members of a genus have a codon composition most similar to other members of its genus, although not necessarily other members of its species.

Genetic code redundancy and its influence on the encoded polypeptides

The genetic code is said to be redundant in that the same amino acid residue can be encoded by multiple, so-called synonymous, codons. If all properties of synonymous codons were entirely equivalent, one would expect that they would be equally distributed along protein coding sequences. However, many studies over the last three decades have demonstrated that their distribution is not entirely random. It has been postulated that certain codons may be translated by the ribosome faster than others and thus their non-random distribution dictates how fast the ribosome moves along particular segments of the mRNA. The reasons behind such segmental variability in the rates of protein synthesis, and thus polypeptide emergence from the ribosome, have been explored by theoretical and experimental approaches. Predictions of the relative rates at which particular codons are translated and their impact on the nascent chain have not arrived at unequivocal conclusions. This is probably due, at least in part, to variation in the basis for classification of codons as “fast” or “slow”, as well as variability in the number and types of genes and proteins analyzed. Recent methodological advances have allowed nucleotide-resolution studies of ribosome residency times in entire transcriptomes, which confirm the non-uniform movement of ribosomes along mRNAs and shed light on the actual determinants of rate control. Moreover, experiments have begun to emerge that systematically examine the influence of variations in ribosomal movement and the fate of the emerging polypeptide chain.

Common polymorphisms in KCNJ5 [corrected] are associated with early-onset lone atrial fibrillation in Caucasians

OBJECTIVES

The aim of this study was to screen lone atrial fibrillation (AF) patients for mutations in the genes KCNJ2, KCNJ3 and KCNJ5, all encoding potassium channels. Furthermore, we wanted to replicate the prior association of two single-nucleotide polymorphisms (SNPs) in KCNJ5, C171T (rs6590357) and G810T (rs7118824), with lone AF in Han Chinese.

METHODS

We sequenced the coding region and splice site of KCNJ2, KCNJ3 and KCNJ5 in 187 early-onset lone-AF patients screening for mutations and counting SNP frequencies for the two noted SNPs in KCNJ5.

RESULTS

No mutations were found in KCNJ2, KCNJ3 or KCNJ5. Both genotype distribution and allele frequencies of the SNPs rs6590357 and rs7118824 significantly differed between the AF and control group (p(genotype) = 0.0067, p(allele) = 0.0021 and p(genotype) = 0.014, p(allele) = 0.0101, respectively). On allele level, the OR for lone AF for rs6590357 was 1.77 (95% CI 1.16-2.73, p = 0.009) and for rs7118824 it was 1.71 (95% CI 1.13-2.57, p = 0.01) in a model adjusted for age and gender.

CONCLUSIONS

Our findings indicate that rs6590357 and rs7118824 in KCNJ5 are associated with early-onset lone AF in Caucasians. No mutations were found in the exon or splice site of KCNJ2, KCNJ3 or KCNJ5.

Screening of KCNN3 in patients with early-onset lone atrial fibrillation

AIMS

The aim of this study was to screen KCNN3 encoding the small-conductance calcium-activated K+ channel (SK3) in lone atrial fibrillation patients. Atrial fibrillation (AF) is the most common cardiac arrhythmia. A genome-wide association study has recently associated an intronic single-nucleotide polymorphism (SNP) in KCNN3 with lone AF.

METHODS AND RESULTS

We sequenced the coding region and splice junctions of KCNN3 in 209 early-onset lone AF patients, screening for variations. A group of 208 healthy blood donors with normal ECGs and without cardiac symptoms were used as controls. All patients and controls were of Danish ethnicity. No mutations were found in the coding regions or splice sites of KCNN3. We found one known exonic synonymous SNP (rs1131820) in KCNN3 that was associated with AF. Both the genotype distribution and allele frequencies of SNP rs1131820 were significantly different between the AF cases and controls (PGenotype=0.047 and PAllele=0.027). Being a homozygous carrier of the major allele (GG) vs. the minor allele (AA) of rs1131820 was associated with an odds ratio of 2.85 (95% CI 1.13-7.18, P=0.026) for lone AF.

CONCLUSIONS

In this study of 209 young lone AF patients, we found no mutations in the exons or splice sites of KCNN3, but we found an association between the synonymous SNP rs1131820 in KCNN3 and lone AF.

Characterization of ovine Toll-like receptor 9 protein coding region, comparative analysis, detection of mutations and maedi visna infection

One of the major roles of innate immunity system is the recognition and the determination of the nature of the antigen. This ability is encompassed by specific receptors as Toll-like receptors (TLRs). TLR9 recognizes bacterial and viral CpG motifs, while their potent immunostimulation effect seems to be promising for lentiviral therapies. Recent studies, however, show the presence of a big polymorphism within the TLR genes and the linkage between substitutions and susceptibility to various infections. Moreover, different recognition ability seems to be utilized by different species and possibly breeds. In this study, we characterized the protein coding region of ovine TLR9 gene. By using comparative analysis of two closely related species and humans, we suggest, which characteristics of protein could be responsible for altered recognition. Furthermore, analyzing the presence of the substitutions, we show the intraspecies polymorphism and its possible implications, while attempting to define the association of discovered substitutions with the maedi visna infection.

Silent (synonymous) SNPs: should we care about them?

One of the surprising findings of the Human Genome Project was that single nucleotide polymorphisms (SNPs), which, by definition, have a minor allele frequency greater than 1%, occur at higher rates than previously suspected. When occurring in the gene coding regions, SNPs can be synonymous (i.e., not causing a change in the amino acid) or nonsynonymous (when the amino acid is altered). It has long been assumed that synonymous SNPs are inconsequential, as the primary sequence of the protein is retained. A number of studies have questioned this assumption over the last decade, showing that synonymous mutations are also under evolutionary pressure and they can be implicated in disease. More importantly, several of the mechanisms by which synonymous mutations alter the structure, function, and expression level of proteins are now being elucidated. Studies have demonstrated that synonymous polymorphisms can affect messenger RNA splicing, stability, and structure as well as protein folding. These changes can have a significant effect on the function of proteins, change cellular response to therapeutic targets, and often explain the different responses of individual patients to a certain medication.

Synonymous mutations and ribosome stalling can lead to altered folding pathways and distinct minima

How can we understand a case in which a given amino acid sequence folds into structurally and functionally distinct molecules? Synonymous single-nucleotide polymorphisms in the MDR1 (multidrug resistance 1 or ABCB1) gene involving frequent-to-rare codon substitutions lead to identical protein sequences. Remarkably, these alternative sequences give a protein product with similar but different structures and functions. Here, we propose that long-enough ribosomal pause time scales may lead to alternate folding pathways and distinct minima on the folding free energy surface. While the conformational and functional differences between the native and alternate states may be minor, the MDR1 case illustrates that the barriers may nevertheless constitute sufficiently high hurdles in physiological time scales, leading to kinetically trapped states with altered structures and functions. Different folding pathways leading to conformationally similar trapped states may be due to swapping of (fairly symmetric) segments. Domain swapping is more likely in the no-pause case in which the chain elongates and folds simultaneously; on the other hand, sufficiently long pause times between such segments may be expected to lessen the chances of swapping events. Here, we review the literature in this light.

In vivo introduction of unpreferred synonymous codons into the Drosophila Adh gene results in reduced levels of ADH protein

The evolution of codon bias, the unequal usage of synonymous codons, is thought to be due to natural selection for the use of preferred codons that match the most abundant species of isoaccepting tRNA, resulting in increased translational efficiency and accuracy. We examined this hypothesis by introducing 1, 6, and 10 unpreferred codons into the Drosophila alcohol dehydrogenase gene (Adh). We observed a significant decrease in ADH protein production with number of unpreferred codons, confirming the importance of natural selection as a mechanism leading to codon bias. We then used this empirical relationship to estimate the selection coefficient (s) against unpreferred synonymous mutations and found the value (s >or= 10(-5)) to be approximately one order of magnitude greater than previous estimates from population genetics theory. The observed differences in protein production appear to be too large to be consistent with current estimates of the strength of selection on synonymous sites in D. melanogaster.

Silent mutations affect in vivo protein folding in Escherichia coli

As an approach to investigate the molecular mechanism of in vivo protein folding and the role of translation kinetics on specific folding pathways, we made codon substitutions in the EgFABP1 (Echinococcus granulosus fatty acid binding protein1) gene that replaced five minor codons with their synonymous major ones. The altered region corresponds to a turn between two short alpha helices. One of the silent mutations of EgFABP1 markedly decreased the solubility of the protein when expressed in Escherichia coli. Expression of this protein also caused strong activation of a reporter gene designed to detect misfolded proteins, suggesting that the turn region seems to have special translation kinetic requirements that ensure proper folding of the protein. Our results highlight the importance of codon usage in the in vivo protein folding.

Depletion of free 30S ribosomal subunits in Escherichia coli by expression of RNA containing Shine-Dalgarno-like sequences

We have constructed synthetic coding sequences for the expression of poly(alpha,L-glutamic acid) (PLGA) as fusion proteins with dihydrofolate reductase (DHFR) in Escherichia coli. These PLGA coding sequences use both GAA and GAG codons for glutamic acid and contain sequence elements (5'-GAGGAGG-3') that resemble the consensus Shine-Dalgarno (SD) sequence found at translation initiation sites in bacterial mRNAs. An unusual feature of DHFR-PLGA expression is that accumulation of the protein is inversely related to the level of induction of its mRNA. Cellular protein synthesis was inhibited >95% by induction of constructs for either translatable or untranslatable PLGA RNAs. Induction of PLGA RNA resulted in the depletion of free 30S ribosomal subunits and the appearance of new complexes in the polyribosome region of the gradient. Unlike normal polyribosomes, these complexes were resistant to breakdown in the presence of puromycin. The novel complexes contained 16S rRNA, 23S rRNA, and PLGA RNA. We conclude that multiple noninitiator SD-like sequences in the PLGA RNA inhibit cellular protein synthesis by sequestering 30S small ribosomal subunits and 70S ribosomes in nonfunctional complexes on the PLGA mRNA.

Terminal proteins essential for the replication of linear plasmids and chromosomes in Streptomyces

Linear plasmids and chromosomes of the bacterial genus Streptomyces have proteins of unknown characteristics and function linked covalently to their 5' DNA termini. We purified protein attached to the end of the pSLA2 linear plasmid of Streptomyces rochei, determined the N-terminal amino acid sequence, and used this information to clone corresponding genes from a S. rochei cosmid library. Three separate terminal protein genes (here designated as tpgR1, tpgR2, and tpgR3), which map to the S. rochei chromosome and to 100-kb and 206-kb linear plasmids contained in S. rochei, were isolated and found to encode a family of similar but distinct 21-kD proteins. Using tpgR1 to probe a genomic DNA library of Streptomyces lividans ZX7, whose linear chromosome can undergo transition to a circular form, we isolated a S. lividans chromosomal gene (tpgL) that we found specifies a protein closely related to, and functionally interchangeable with, TpgR proteins for pSLA2 maintenance in S. lividans. Mutation of tpgL precluded propagation of the pSLA2 plasmid in a linear form and also prevented propagation of S. lividans cells that contain linear, but not circular, chromosomes, indicating a specific and essential role for tpg genes in linear DNA replication. Surprisingly, Tpg proteins were observed to contain a reverse transcriptase-like domain rather than sequences in common with proteins that attach covalently to the termini of linear DNA replicons.

NMR structure of the N-terminal J domain of murine polyomavirus T antigens. Implications for DnaJ-like domains and for mutations of T antigens

The NMR structure of the N-terminal, DnaJ-like domain of murine polyomavirus tumor antigens (PyJ) has been determined to high precision, with root mean square deviations to the mean structure of 0.38 A for backbone atoms and 0.94 A for all heavy atoms of ordered residues 5-41 and 50-69. PyJ possesses a three-helix fold, in which anti-parallel helices II and III are bridged by helix I, similar to the four-helix fold of the J domains of DnaJ and human DnaJ-1. PyJ differs significantly in the lengths of N terminus, helix I, and helix III. The universally conserved HPD motif appears to form a His-Pro C-cap of helix II. Helix I features a stabilizing Schellman C-cap that is probably conserved universally among J domains. On the helix II surface where positive charges of other J domains have been implicated in binding of hsp70s, PyJ contains glutamine residues. Nonetheless, chimeras that replace the J domain of DnaJ with PyJ function like wild-type DnaJ in promoting growth of Escherichia coli. This activity can be modulated by mutations of at least one of these glutamines. T antigen mutations reported to impair cellular transformation by the virus, presumably via interactions with PP2A, cluster in the hydrophobic folding core and at the extreme N terminus, remote from the HPD loop.

Ribosome traffic in E. coli and regulation of gene expression

The ribosome traffic during translation of E. coli coding sequences was simulated, assuming that the rate of translation of individual codons is limited by the cognate tRNA availability. Actual translation rates were taken from Solomovici et al. (J. theor. Biol. 185, 511-521, 1997). The mean translation rates of the 4271 sequences cover a broad, two-fold range, whereas the local rate of translation along messengers varies three-fold on average. The simulation allows one to sketch the ribosome traffic on the polysome, in particular by providing the extent of mRNA sequences uncovered between consecutive ribosomes and the time during which these sequences are exposed. These parameters may participate in the control of mRNA stability and transcriptional polarity. By averaging the translation rates in a 17-codon window, assumed to be the sequence covered by a translating ribosome, and sliding this window along a given coding sequence, the addresses KMAX and KMIN, and the times TMAX and TMIN of respectively the slowest and the fastest translated window were determined. It is shown that under the assumptions made, TMAX sets the number of proteins translated from a given mRNA molecule per unit time, in case the delay between consecutive translation starts is below TMAX. Both windows display two strong biases, one as expected on the usage of codon frequencies, and the other surprisingly on the occurrence of amino acids.