Initiation of translation in prokaryotes and eukaryotes

Stepwise modifications of genetic parts reinforce the secretory production of nattokinase in Bacillus subtilis

Nattokinase (NK) is an important serine‐protease with direct fibrinolytic activity involving the prevention of cardiovascular disease as an antithrombotic agent. Dozens of studies have focused on the characterization of intrinsic novel promoters and signal peptides to the secretory production of recombinant proteins in Bacillus subtilis. However, intrinsic genetic elements have several drawbacks, which cannot mediate the production of NK to the desired level. In this study, the genetic elements, which were used to overproduce the recombinant secretory NK, were rationally modified in B. subtilis in a stepwise manner. The first step was to select a suitable signal peptide for the highly efficient secretion of NK. By comparison of the secretory levels mediated by two different signal peptides, which were encoded by the genes of a minor extracellular protease epr (SP_epr) and cell‐wall associated protease wapA (SP_wapA), respectively, SP_wapA was verified as the superior secretory element. Second, P04, which was a synthetic promoter screened from an array of mutants based on the promoter cloned from the operon of a quorum‐sensing associated gene srfA (P_srfA), was paired to SP_wapA. The secretory level of NK was obviously augmented by the combination of these two genetic elements. Third, the cis‐acting element CodY‐binding sequence positioned at the 5′UTR was deleted (yielding P08), and thus the secretory level was significantly elevated. The activity of NK, which was defined as fibrinolytic units (FU), reached to a level of 270 FU ml⁻¹. Finally, the superior genetic element composed of P08 and SP_wapA was utilized to overproduce NK in the host B. subtilis WB800, which was able to produce the secretory NK at 292 FU ml⁻¹. The strategy established in this study can not only be used to overproduce NK in B. subtilis but also might be a promising pipeline to modify the genetic element for the synthetic secretory system.

Bioinformatics and Translation Elongation

Codon usage depends on mutation bias, tRNA-mediated selection, and the need for high efficiency and accuracy in translation. One codon in a synonymous codon family is often strongly over-used, especially in highly expressed genes, which often leads to a high dN/dS ratio because dS is very small. Many different codon usage indices have been proposed to measure codon usage and codon adaptation. Sense codon could be misread by release factors and stop codons misread by tRNAs, which also contribute to codon usage in rare cases. This chapter outlines the conceptual framework on codon evolution, illustrates codon-specific and gene-specific codon usage indices, and presents their applications. A new index for codon adaptation that accounts for background mutation bias (Index of Translation Elongation) is presented and contrasted with codon adaptation index (CAI) which does not consider background mutation bias. They are used to re-analyze data from a recent paper claiming that translation elongation efficiency matters little in protein production. The reanalysis disproves the claim.

TASEP modelling provides a parsimonious explanation for the ability of a single uORF to derepress translation during the integrated stress response

Translation initiation is the rate-limiting step of protein synthesis that is downregulated during the Integrated Stress Response (ISR). Previously, we demonstrated that most human mRNAs that are resistant to this inhibition possess translated upstream open reading frames (uORFs), and that in some cases a single uORF is sufficient for the resistance. Here we developed a computational model of Initiation Complexes Interference with Elongating Ribosomes (ICIER) to gain insight into the mechanism. We explored the relationship between the flux of scanning ribosomes upstream and downstream of a single uORF depending on uORF features. Paradoxically, our analysis predicts that reducing ribosome flux upstream of certain uORFs increases initiation downstream. The model supports the derepression of downstream translation as a general mechanism of uORF-mediated stress resistance. It predicts that stress resistance can be achieved with long slowly decoded uORFs that do not favor translation reinitiation and that start with initiators of low leakiness.

Ribosome Shunting, Polycistronic Translation, and Evasion of Antiviral Defenses in Plant Pararetroviruses and Beyond

Viruses have compact genomes and usually translate more than one protein from polycistronic RNAs using leaky scanning, frameshifting, stop codon suppression or reinitiation mechanisms. Viral (pre-)genomic RNAs often contain long 5′-leader sequences with short upstream open reading frames (uORFs) and secondary structure elements, which control both translation initiation and replication. In plants, viral RNA and DNA are targeted by RNA interference (RNAi) generating small RNAs that silence viral gene expression, while viral proteins are recognized by innate immunity and autophagy that restrict viral infection. In this review we focus on plant pararetroviruses of the family Caulimoviridae and describe the mechanisms of uORF- and secondary structure-driven ribosome shunting, leaky scanning and reinitiation after translation of short and long uORFs. We discuss conservation of these mechanisms in different genera of Caulimoviridae, including host genome-integrated endogenous viral elements, as well as in other viral families, and highlight a multipurpose use of the highly-structured leader sequence of plant pararetroviruses in regulation of translation, splicing, packaging, and reverse transcription of pregenomic RNA (pgRNA), and in evasion of RNAi. Furthermore, we illustrate how targeting of several host factors by a pararetroviral effector protein can lead to transactivation of viral polycistronic translation and concomitant suppression of antiviral defenses. Thus, activation of the plant protein kinase target of rapamycin (TOR) by the Cauliflower mosaic virus transactivator/viroplasmin (TAV) promotes reinitiation of translation after long ORFs on viral pgRNA and blocks antiviral autophagy and innate immunity responses, while interaction of TAV with the plant RNAi machinery interferes with antiviral silencing.

Genome scale analysis of Escherichia coli with a comprehensive prokaryotic sequence-based biophysical model of translation initiation and elongation

Translation initiation in prokaryotes is affected by the mRNA folding and interaction of the ribosome binding site with the ribosomal RNA. The elongation rate is affected, among other factors, by the local biophysical properties of the coding regions, the decoding rates of different codons, and the interactions among ribosomes. Currently, there is no comprehensive biophysical model of translation that enables the prediction of mRNA translation dynamics based only on the transcript sequence and while considering all of these fundamental aspects of translation. In this study, we provide, for the first time, a computational simulative biophysical model of both translation initiation and elongation with all aspects mentioned above. We demonstrate our model performance and advantages focusing on Escherichia coli genes. We further show that the model enables prediction of translation rate, protein levels, and ribosome densities. In addition, our model enables quantifying the effect of silent mutations on translation rate in different parts of the transcript, the relative effect of mutations on translation initiation and elongation, and the effect of mutations on ribosome traffic jams. Thus, unlike previous models, the proposed one provides comprehensive information, facilitating future research in disciplines such as molecular evolution, synthetic biology, and functional genomics. A toolkit to estimate translation dynamics of transcripts is available at: https://www.cs.tau.ac.il/∼tamirtul/transim.

Two I84 family protease inhibitors from Chinese razor clams Sinonovacula constricta expressed in response to environmental challenges

Protease inhibitors play critical roles in numerous biological processes including host defense in all multicellular organisms. Eighty three evolutionary families of protease inhibitors are currently accommodated in the MEROPS database and the I84 family currently consists of 3 novel serine protease inhibitors from the eastern oyster Crassostrea virginica. In this study, we identified 2 new I84 family members from the Chinese razor clam Sinonovacula constricta, scSI-1 and scSI-2, using cDNA cloning and sequencing. The scSI-1 cDNA consisted of 494 bp with a 282 bp ORF encoding a 93-amino acid polypeptide that was predicted to have a 19-amino acid signal peptide and a 74-residue mature protein with a calculated molecular mass of 8248.5 Da. The scSI-2 cDNA was 490 bp long with a 273 bp ORF encoding a 90-amino acid polypeptide that was predicted to have an 18-amino acid signal peptide and a 72-residue nature protein with a calculated molecular mass of 7528.4 Da. ScSI-1 and scSI-2 shared high sequence similarity with the 3 known members of I84 family and both expressed primarily in the clam digestive glands. Protease inhibitory activity in the clam plasma also exhibited the signature kinetic characteristics of the I84 members from the oyster. In addition, levels of scSI-1 and scSI-2 gene expression in digestive glands and the protease inhibitory activity in plasma elevated significantly in clams challenged by bacterial injections and Vibrio harveyi was more effective than Staphylococcus epidermidis in inducing the gene expression and plasma protease inhibitory activity. Moreover, drastic changes of salinity and temperature also caused significant changes in the gene expression and plasma activity. These results indicated that scSI-1 and scSI-2 represented 2 new members of the I84 family and they likely play a role in clam host defense against infections and in reactions against physiochemical stressors.

Complex regulation of the TaMyc1 gene expression in wheat grain synthesizing anthocyanin pigments

The wheat TaMyc1 gene encodes for transcriptional factor (TF) with bHLH domain. The gene is expressed in purple wheat grains and activates transcription of the anthocyanin biosynthesis structural genes. To reveal the features of TaMyc1 regulation in wheat pericarp transcription start sites (TSS) were identified by 5' RACE mean and translation efficiency was predicted by in silico methods. Three alternative transcript variants of TaMyc1 differing by 5' leader sequence only were identified in purple pericarp. The three transcripts are generated from distinct TATA boxes and thereby are differed by TSS. Two transcripts (TaMyc1a, -b) have identical initiation AUG codons that lead to the TaMYC1 regulatory protein with bHLH domain. However because of different stability of secondary structures predicted in 5' leader the two transcripts might be translated with different efficiency. The third transcript is assumed to be not effectively translated. qRT-PCR and colonies counting were applied to assess contribution each of the transcripts to total TaMyc1 gene transcription level. TaMyc1c has the lowest contribution (ca. 16%), whereas the others two transcripts contribute equally (ca. 42%) to total TaMyc1 expression level. The role of the tree mRNA isoforms transcribed in one tissue is discussed.

An early nonsense mutation facilitates the expression of a short isoform of CNGA3 by alternative translation initiation

The cyclic nucleotide-gated (CNG) channel - composed of CNGA3 and CNGB3 subunits - mediates the influx of cations in cone photoreceptors after light stimulation and thus is a key element in cone phototransduction. Mutations in CNGA3 and CNGB3 are associated with achromatopsia, a rare autosomal recessive retinal disorder. Here, we demonstrate that the presence of an early nonsense mutation in CNGA3 induces the usage of a downstream alternative translation initiation site giving rise to a short CNGA3 isoform. The expression of this short isoform was verified by Western blot analysis and DAB staining of HEK293 cells and cone photoreceptor-like 661W cells expressing CNGA3-GST fusion constructs. Functionality of the short isoform was confirmed by a cellular calcium influx assay. Furthermore, patients carrying an early nonsense mutation were analyzed for residual cone photoreceptor function in order to identify a potential role of the short isoform to modify the clinical outcome in achromatopsia patients. Yet the results suggest that the short isoform is not able to compensate for the loss of the long isoform leaving the biological role of this variant unclear.

Enhanced Solubilization of Class B Radical S-Adenosylmethionine Methylases by Improved Cobalamin Uptake in Escherichia coli

The methylation of unactivated carbon and phosphorus centers is a burgeoning area of biological chemistry, especially given that such reactions constitute key steps in the biosynthesis of numerous enzyme cofactors, antibiotics, and other natural products of clinical value. These kinetically challenging reactions are catalyzed exclusively by enzymes in the radical S-adenosylmethionine (SAM) superfamily and have been grouped into four classes (A-D). Class B radical SAM (RS) methylases require a cobalamin cofactor in addition to the [4Fe-4S] cluster that is characteristic of RS enzymes. However, their poor solubility upon overexpression and their generally poor turnover has hampered detailed in vitro studies of these enzymes. It has been suggested that improper folding, possibly caused by insufficient cobalamin during their overproduction in Escherichia coli, leads to formation of inclusion bodies. Herein, we report our efforts to improve the overproduction of class B RS methylases in a soluble form by engineering a strain of E. coli to take in more cobalamin. We cloned five genes ( btuC, btuE, btuD, btuF, and btuB) that encode proteins that are responsible for cobalamin uptake and transport in E. coli and co-expressed these genes with those that encode TsrM, Fom3, PhpK, and ThnK, four class B RS methylases that suffer from poor solubility during overproduction. This strategy markedly enhances the uptake of cobalamin into the cytoplasm and improves the solubility of the target enzymes significantly.

Functional genetic variant in the Kozak sequence of WW domain-containing oxidoreductase (WWOX) gene is associated with oral cancer risk

In Taiwan, oral cancer is the fourth leading cancer in males and is associated with exposure to environmental carcinogens. WW domain-containing oxidoreductase (WWOX), a tumor suppressor gene, is associated with the development of various cancers. We hypothesized that genetic variants of WWOX influence the susceptibility to oral cancer. Five polymorphisms of WWOX gene from 761 male patients with oral cancer and 1199 male cancer-free individuals were genotyped. We observed that individuals carrying the polymorphic allele of WWOX rs11545028 are more susceptible to oral cancer. Furthermore, patients with advanced-stage oral cancer were associated with a higher frequency of WWOX rs11545028 polymorphisms with the variant genotype TT than did patients with the wild-type gene. An additional integrated in silico analysis confirmed that rs11545028 affects WWOX expression, which significantly correlates with tumor expression and subsequently with tumor development and aggressiveness. In conclusion, genetic variants of WWOX contribute to the occurrence of oral cancer, and the findings regarding these biomarkers provided a prediction model for risk assessment.

The Regulation of Translation in Alphavirus-Infected Cells

Sindbis virus (SINV) contains an RNA genome of positive polarity with two open reading frames (ORFs). The first ORF is translated from the genomic RNA (gRNA), rendering the viral non-structural proteins, whereas the second ORF is translated from a subgenomic mRNA (sgRNA), which directs the synthesis of viral structural proteins. SINV infection strongly inhibits host cell translation through a variety of different mechanisms, including the phosphorylation of the eukaryotic initiation factor eIF2α and the redistribution of cellular proteins from the nucleus to the cytoplasm. A number of motifs have been identified in SINV sgRNA, including a hairpin downstream of the AUG initiation codon, which is involved in the translatability of the viral sgRNA when eIF2 is inactivated. Moreover, a 3′-UTR motif containing three stem-loop structures is involved in the enhancement of translation in insect cells, but not in mammalian cells. Accordingly, SINV sgRNA has evolved several structures to efficiently compete for the cellular translational machinery. Mechanistically, sgRNA translation involves scanning of the 5′-UTR following a non-canonical mode and without the requirement for several initiation factors. Indeed, sgRNA-directed polypeptide synthesis occurs even after eIF4G cleavage or inactivation of eIF4A by selective inhibitors. Remarkably, eIF2α phosphorylation does not hamper sgRNA translation during the late phase of SINV infection. SINV sgRNA thus constitutes a unique model of a capped viral mRNA that is efficiently translated in the absence of several canonical initiation factors. The present review will mainly focus in the non-canonical mechanism of translation of SINV sgRNA.

RNA Structure Design Improves Activity and Specificity of trans-Splicing-Triggered Cell Death in a Suicide Gene Therapy Approach

Spliceosome-mediated RNA trans-splicing enables correction or labeling of pre-mRNA, but therapeutic applications are hampered by issues related to the activity and target specificity of trans-splicing RNA (tsRNA). We employed computational RNA structure design to improve both on-target activity and specificity of tsRNA in a herpes simplex virus thymidine kinase/ganciclovir suicide gene therapy approach targeting alpha fetoprotein (AFP), a marker of hepatocellular carcinoma (HCC) or human papillomavirus type 16 (HPV-16) pre-mRNA. While unstructured, mismatched target binding domains significantly improved 3′ exon replacement (3’ER), 5′ exon replacement (5’ER) correlated with the thermodynamic stability of the tsRNA 3′ end. Alternative on-target trans-splicing was found to be a prevalent event. The specificity of trans-splicing with the intended target splice site was improved 10-fold by designing tsRNA that harbors secondary target binding domains shielding alternative on-target and blinding off-target splicing events. Such rationally designed suicide RNAs efficiently triggered death of HPV-16-transduced or hepatoblastoma-derived human tissue culture cells without evidence for off-target cell killing. Highest cell death activities were observed with novel dual-targeting tsRNAs programmed for trans-splicing toward AFP and a second HCC pre-mRNA biomarker. Our observations suggest trans-splicing represents a promising approach to suicide gene therapy.

Exploring Cytotoxic mRNAs as a Novel Class of Anti-cancer Biotherapeutics

New treatments to overcome the obstacles of conventional anti-cancer therapy are a permanent subject of investigation. One promising approach is the application of toxins linked to cell-specific ligands, so-called immunotoxins. Another attractive option is the employment of toxin-encoding plasmids. However, immunotoxins cause hepatoxicity, and DNA therapeutics, among other disadvantages, bear the risk of insertional mutagenesis. As an alternative, this study examined chemically modified mRNAs coding for diphtheria toxin, subtilase cytotoxin, and abrin-a for their ability to reduce cancer cell growth both in vitro and in vivo. The plant toxin abrin-a was the most promising candidate among the three tested toxins and was further investigated. Its expression was demonstrated by western blot. Experiments with firefly luciferase in reticulocyte lysates and co-transfection experiments with EGFP demonstrated the capability of abrin-a to inhibit protein synthesis. Its cytotoxic effect was quantified employing viability assays and propidium iodide staining. By studying caspase-3/7 activation, Annexin V-binding, and chromatin condensation with Hoechst33258 staining, apoptotic cell death could be confirmed. In mice, repeated intratumoral injections of complexed abrin-a mRNA resulted in a significant reduction (89%) of KB tumor size compared to a non-translatable control mRNA.

The gene fmt, encoding tRNAfMet-formyl transferase, is essential for normal growth of M. bovis, but not for viability

Mycobacterium tuberculosis is a major health threat, necessitating novel drug targets. Protein synthesis in bacteria uses initiator tRNA_i charged with formylated methionine residue. Deletion of the formylase gene, tRNA^fMet-formyl transferase (fmt), causes severe growth-retardation in E. coli and in S. pneumoniae, but not in P. aeruginosa or S. aureus. fmt was predicted to be essential in M. tuberculosis by transposon library analysis, but this was never formally tested in any mycobacteria. We performed a targeted deletion of fmt in M. smegmatis as well as Mtb-complex (M. bovis). In both cases, we created a mero-diploid strain, deleted the native gene by two-step allelic exchange or specialized-phage transduction, and then removed the complementing gene to create full deletion mutants. In M. smegmatis a full deletion strain could be easily created. In contrast, in M. bovis-BCG, a full deletion strain could only be created after incubation of 6 weeks, with a generation time ~2 times longer than for wt bacteria. Our results confirm the importance of this gene in pathogenic mycobacteria, but as the deletion mutant is viable, validity of fmt as a drug target remains unclear. Our results also refute the previous reports that fmt is essential in M. tuberculosis-complex.

Genomic and Transcriptomic Analysis Reveals Spliced Leader Trans-Splicing in Cryptomonads

Spliced leader trans-splicing (SLTS) is a poorly understood mechanism that is found in a diversity of eukaryotic lineages. In SLTS, a short RNA sequence is added near the 5′ ends of the transcripts of protein-coding genes by a modified spliceosomal reaction. Available data suggest that SLTS has evolved many times, and might be more likely to evolve in animals. That SLTS might be more likely to evolve in the context of the generally complex transcriptomes characteristic of animals suggests the possibility that SLTS functions in gene regulation or transcriptome diversification, however no general novel function for SLTS is known. Here, I report SLTS in a lineage of cellularly complex unicellular eukaryotes. Cryptomonads are a group of eukaryotic algae that acquired photosynthetic capacity by secondary endosymbiosis of a red alga, and that retain a reduced copy of the nucleus of the engulfed alga. I estimate that at least one-fifth of genes in the model cryptomonad Guillardia theta and its relative Hanusia phi undergo SLTS. I show that hundreds of genes in G. theta generate alternative transcripts by SLTS at alternative sites, however I find little evidence for alternative protein production by alternative SLTS splicing. Interestingly, I find no evidence for substantial operon structure in the G. theta genome, in contrast to previous findings in other lineages with SLTS. These results extend SLTS to another major group of eukaryotes, and heighten the mystery of the evolution of SLTS and its association with cellular and transcriptomic complexity.

Differential use of 3'CITEs by the subgenomic RNA of Pea enation mosaic virus 2

The genomic RNA (gRNA) of Pea enation mosaic virus 2 (PEMV2) is the template for p33 and -1 frameshift product p94. The PEMV2 subgenomic RNA (sgRNA) encodes two overlapping ORFs, p26 and p27, which are required for movement and stability of the gRNA. Efficient translation of p33 requires two of three 3' proximal cap-independent translation enhancers (3'CITEs): the kl-TSS, which binds ribosomes and engages in a long-distance interaction with the 5'end; and the adjacent eIF4E-binding PTE. Unlike the gRNA, all three 3'CITEs were required for efficient translation of the sgRNA, which included the ribosome-binding 3'TSS. A hairpin in the 5' proximal coding region of p26/p27 supported translation by the 3'CITEs by engaging in a long-distance RNA:RNA interaction with the kl-TSS. These results strongly suggest that the 5' ends of PEMV2 gRNA and sgRNA connect with the 3'UTR through similar long-distance interactions while having different requirements for 3'CITEs.

Switching Protein Localization by Site-Directed RNA Editing under Control of Light

Site directed RNA editing is an engineered tool for the posttranscriptional manipulation of RNA and proteins. Here, we demonstrate the inclusion of additional N- and C-terminal protein domains in an RNA editing-dependent manner to switch between protein isoforms in mammalian cell culture. By inclusion of localization signals, a switch of the subcellular protein localization was achieved. This included the shift from the cytoplasm to the outer-membrane, which typically is inaccessible at the protein-level. Furthermore, the strategy allows to implement photocaging to achieve spatiotemporal control of isoform switching. The strategy does not require substantial genetic engineering, and might well complement current optogenetic and optochemical approaches.

Chronicles of a reductase: Biochemistry, genetics and physio-pathological role of GSNOR

S-nitrosylation is a major redox posttranslational modification involved in cell signaling. The steady state concentration of S-nitrosylated proteins depends on the balance between the relative ability to generate nitric oxide (NO) via NO synthase and to reduce nitrosothiols by denitrosylases. Numerous works have been published in last decades regarding the role of NO and S-nitrosylation in the regulation of protein structure and function, and in driving cellular activities in vertebrates. Notwithstanding an increasing number of observations indicates that impairment of denitrosylation equally affects cellular homeostasis, there is still no report providing comprehensive knowledge on the impact that denitrosylation has on maintaining correct physiological processes and organ activities. Among denitrosylases, S-nitrosoglutathione reductase (GSNOR) represents the prototype enzyme to disclose how denitrosylation plays a crucial role in tuning NO-bioactivity and how much it deeply impacts on cell homeostasis and human patho-physiology. In this review we attempt to illustrate the history of GSNOR discovery and provide the evidence so far reported in support of GSNOR implications in development and human disease.

Ribosome signatures aid bacterial translation initiation site identification

BACKGROUND

While methods for annotation of genes are increasingly reliable, the exact identification of translation initiation sites remains a challenging problem. Since the N-termini of proteins often contain regulatory and targeting information, developing a robust method for start site identification is crucial. Ribosome profiling reads show distinct patterns of read length distributions around translation initiation sites. These patterns are typically lost in standard ribosome profiling analysis pipelines, when reads from footprints are adjusted to determine the specific codon being translated.

RESULTS

Utilising these signatures in combination with nucleotide sequence information, we build a model capable of predicting translation initiation sites and demonstrate its high accuracy using N-terminal proteomics. Applying this to prokaryotic translatomes, we re-annotate translation initiation sites and provide evidence of N-terminal truncations and extensions of previously annotated coding sequences. These re-annotations are supported by the presence of structural and sequence-based features next to N-terminal peptide evidence. Finally, our model identifies 61 novel genes previously undiscovered in the Salmonella enterica genome.

CONCLUSIONS

Signatures within ribosome profiling read length distributions can be used in combination with nucleotide sequence information to provide accurate genome-wide identification of translation initiation sites.

Measurements of translation initiation from all 64 codons in E. coli

Our understanding of translation underpins our capacity to engineer living systems. The canonical start codon (AUG) and a few near-cognates (GUG, UUG) are considered as the ‘start codons’ for translation initiation in Escherichia coli. Translation is typically not thought to initiate from the 61 remaining codons. Here, we quantified translation initiation of green fluorescent protein and nanoluciferase in E. coli from all 64 triplet codons and across a range of DNA copy number. We detected initiation of protein synthesis above measurement background for 47 codons. Translation from non-canonical start codons ranged from 0.007 to 3% relative to translation from AUG. Translation from 17 non-AUG codons exceeded the highest reported rates of non-cognate codon recognition. Translation initiation from non-canonical start codons may contribute to the synthesis of peptides in both natural and synthetic biological systems.

Direct Head-to-Head Evaluation of Recombinant Adeno-associated Viral Vectors Manufactured in Human versus Insect Cells

The major drawback of the Baculovirus/Sf9 system for recombinant adeno-associated viral (rAAV) manufacturing is that most of the Bac-derived rAAV vector serotypes, with few exceptions, demonstrate altered capsid compositions and lower biological potencies. Here, we describe a new insect cell-based production platform utilizing attenuated Kozak sequence and a leaky ribosome scanning to achieve a serotype-specific modulation of AAV capsid proteins stoichiometry. By way of example, rAAV5 and rAAV9 were produced and comprehensively characterized side by side with HEK293-derived vectors. A mass spectrometry analysis documented a 3-fold increase in both viral protein (VP)1 and VP2 capsid protein content compared with human cell-derived vectors. Furthermore, we conducted an extensive analysis of encapsidated single-stranded viral DNA using next-generation sequencing and show a 6-fold reduction in collaterally packaged contaminating DNA for rAAV5 produced in insect cells. Consequently, the re-designed rAAVs demonstrated significantly higher biological potencies, even in a comparison with HEK293-manufactured rAAVs mediating, in the case of rAAV5, 4-fold higher transduction of brain tissues in mice. Thus, the described system yields rAAV vectors of superior infectivity and higher genetic identity providing a scalable platform for good manufacturing practice (GMP)-grade vector production.

Role of Different CD40 Polymorphisms in Graves' Disease and Hashimoto's Thyroiditis

Genome-wide association studies have led to the discovery of several susceptibility genes related to autoimmune thyroid diseases (AITDs). However, controversial results have been reported regarding the role of single-nucleotide polymorphism (SNP) of CD40 in the disease susceptibility. The objective of this study was to identify the relationship of the polymorphisms of three sites of CD40 with the susceptibility to AITD in the Chinese population. We genotyped three polymorphisms of CD40: C/T -1 SNP, 58038T site of the third exon and C64610G site of the ninth exon in 196 GD cases, 121 HT cases and 122 control subjects. The three putative polymorphism sites were amplified by PCR for sequencing and analysis. The genotype frequencies of CD40 -1 C/C genotype and C allele were significantly higher in the GD group than those in normal control. For the C64610G polymorphism, the C/G genotype was significantly more frequent in HT group than in control group, and the G allele frequencies in the GD and HT group were both higher than those in control group. These results indicated that there exist different susceptibility loci for AITD within CD40, each contributing a different effect in the onset and development of AITDs.

Consensus architecture of promoters and transcription units in Escherichia coli: design principles for synthetic biology

Genetic information in genomes is ordered, arranged in such a way that it constitutes a code, the so-called cis regulatory code. The regulatory machinery of the cell, termed trans-factors, decodes and expresses this information. In this way, genomes maintain a potential repertoire of genetic programs, parts of which are executed depending on the presence of active regulators in each condition. These genetic programs, executed by the regulatory machinery, have functional units in the genome delimited by punctuation-like marks. In genetic terms, these informational phrases correspond to transcription units, which are the minimal genetic information expressed consistently from initiation to termination marks. Between the start and final punctuation marks, additional marks are present that are read by the transcriptional and translational machineries. In this work, we look at all the experimentally described and predicted genetic elements in the bacterium Escherichia coli K-12 MG1655 and define a comprehensive architectural organization of transcription units to reveal the natural genome-design and to guide the construction of synthetic genetic programs.

ROC'n'Ribo: Characterizing a Riboswitching Expression System by Modeling Single-Cell Data

RNA-engineered systems offer simple and versatile control over gene expression in many organisms. In particular, the design and implementation of riboswitches presents a unique opportunity to manipulate any reporter device in cis, executing tight temporal and spatial control at low metabolic costs. Assembled to higher order genetic circuits, such riboswitch-regulated devices may efficiently process logical operations. Here, we propose a hierarchical stochastic modeling approach to characterize an in silico repressor gate based on neomycin- and tetracycline-sensitive riboswitches. The model was calibrated on rich, transient in vivo single-cell data to account for cell-to-cell variability. To capture the effect of this variability on gate performance we employed the well-known ROC-analysis and derived a novel performance indicator for logic gates. Introduction of such a performance measure is necessary, since we aimed to assess the correct functionality of the gate at the single-cell level-a prerequisite for its further adaption to a genetic circuitry. Our results may be applied to other genetic devices to analyze their efficiency and ensure their correct performance in the light of cell-to-cell variability.

Auxin Signaling in Regulation of Plant Translation Reinitiation

The mRNA translation machinery directs protein production, and thus cell growth, according to prevailing cellular and environmental conditions. The target of rapamycin (TOR) signaling pathway-a major growth-related pathway-plays a pivotal role in optimizing protein synthesis in mammals, while its deregulation triggers uncontrolled cell proliferation and the development of severe diseases. In plants, several signaling pathways sensitive to environmental changes, hormones, and pathogens have been implicated in post-transcriptional control, and thus far phytohormones have attracted most attention as TOR upstream regulators in plants. Recent data have suggested that the coordinated actions of the phytohormone auxin, Rho-like small GTPases (ROPs) from plants, and TOR signaling contribute to translation regulation of mRNAs that harbor upstream open reading frames (uORFs) within their 5'-untranslated regions (5'-UTRs). This review will summarize recent advances in translational regulation of a specific set of uORF-containing mRNAs that encode regulatory proteins-transcription factors, protein kinases and other cellular controllers-and how their control can impact plant growth and development.

Time-Resolved Proteomics Extends Ribosome Profiling-Based Measurements of Protein Synthesis Dynamics

Ribosome profiling is a widespread tool for studying translational dynamics in human cells. Its central assumption is that ribosome footprint density on a transcript quantitatively reflects protein synthesis. Here, we test this assumption using pulsed-SILAC (pSILAC) high-accuracy targeted proteomics. We focus on multiple myeloma cells exposed to bortezomib, a first-line chemotherapy and proteasome inhibitor. In the absence of drug effects, we found that direct measurement of protein synthesis by pSILAC correlated well with indirect measurement of synthesis from ribosome footprint density. This correlation, however, broke down under bortezomib-induced stress. By developing a statistical model integrating longitudinal proteomic and mRNA-sequencing measurements, we found that proteomics could directly detect global alterations in translational rate caused by bortezomib; these changes are not detectable by ribosomal profiling alone. Further, by incorporating pSILAC data into a gene expression model, we predict cell-stress specific proteome remodeling events. These results demonstrate that pSILAC provides an important complement to ribosome profiling in measuring proteome dynamics.

Intracellular trafficking pathways of Cx43 gap junction channels

Gap Junction (GJ) channels, including the most common Connexin 43 (Cx43), have fundamental roles in excitable tissues by facilitating rapid transmission of action potentials between adjacent cells. For instance, synchronization during each heartbeat is regulated by these ion channels at the cardiomyocyte cell-cell border. Cx43 protein has a short half-life, and rapid synthesis and timely delivery of those proteins to particular subdomains are crucial for the cellular organization of gap junctions and maintenance of intracellular coupling. Impairment in gap junction trafficking contributes to dangerous complications in diseased hearts such as the arrhythmias of sudden cardiac death. Of recent interest are the protein-protein interactions with the Cx43 carboxy-terminus. These interactions have significant impact on the full length Cx43 lifecycle and also contribute to trafficking of Cx43 as well as possibly other functions. We are learning that many of the known non-canonical roles of Cx43 can be attributed to the recently identified six endogenous Cx43 truncated isoforms which are produced by internal translation. In general, alternative translation is a new leading edge for proteome expansion and therapeutic drug development. This review highlights recent mechanisms identified in the trafficking of gap junction channels, involvement of other proteins contributing to the delivery of channels to the cell-cell border, and understanding of possible roles of the newly discovered alternatively translated isoforms in Cx43 biology. This article is part of a Special Issue entitled: Gap Junction Proteins edited by Jean Claude Herve.

A silent mutation in human alpha-A crystallin gene in patients with age-related nuclear or cortical cataract

A cataract is a complex multifactorial disease that results from alterations in the cellular architecture, i.e. lens proteins. Genes associated with the development of lens include crystallin genes. Although crystallins are highly conserved proteins among vertebrates, a significant number of polymorphisms exist in human population. In this study, we screened for polymorphisms in crystallin alpha A (CRYAA) and alpha B (CRYAB) genes in 200 patients over 40 years of age, diagnosed with age-related cataract (ARC; nuclear and cortical cataracts). Genomic DNA was extracted from the peripheral blood. The coding regions of the CRYAA and CRYAB gene were amplified using polymerase chain reaction and subjected to restriction digestion. Restriction fragment length polymorphism (RFLP) was performed using known restriction enzymes for CRYAA and CRYAB genes. Denaturing high performance liquid chromatography and direct sequencing were performed to detect sequence variation in CRYAA gene. In silico analysis of secondary CRYAA mRNA structure was performed using CLC RNA Workbench. RFLP analysis did not show any changes in the restriction sites of CRYAA and CRYAB genes. In 6 patients (4 patients with nuclear cataract and 2 with cortical cataract), sequence analysis of the exon 1 in the CRYAA gene showed a silent single nucleotide polymorphism [D2D] (CRYAA: C to T transition). One of the patients with nuclear cataract was homozygous for this allele. The in silico analysis revealed that D2D mutation results in a compact CRYAA mRNA secondary structure, while the wild type CRYAA mRNA has a weak or loose secondary structure. D2D mutation in the CRYAA gene may be an additional risk factor for progression of ARC.

Translation regulation in plants: an interesting past, an exciting present and a promising future

Changes in gene expression are at the core of most biological processes, from cell differentiation to organ development, including the adaptation of the whole organism to the ever-changing environment. Although the central role of transcriptional regulation is solidly established and the general mechanisms involved in this type of regulation are relatively well understood, it is clear that regulation at a translational level also plays an essential role in modulating gene expression. Despite the large number of examples illustrating the critical role played by translational regulation in determining the expression levels of a gene, our understanding of the molecular mechanisms behind such types of regulation has been slow to emerge. With the recent development of high-throughput approaches to map and quantify different critical parameters affecting translation, such as RNA structure, protein-RNA interactions and ribosome occupancy at the genome level, a renewed enthusiasm toward studying translation regulation is warranted. The use of these new powerful technologies in well-established and uncharacterized translation-dependent processes holds the promise to decipher the likely complex and diverse, but also fascinating, mechanisms behind the regulation of translation.

PTENβ is an alternatively translated isoform of PTEN that regulates rDNA transcription

PTEN is a critical tumour suppressor that is frequently mutated in human cancer. We have previously identified a CUG initiated PTEN isoform designated PTENα, which functions in mitochondrial bioenergetics. Here we report the identification of another N-terminal extended PTEN isoform, designated PTENβ. PTENβ translation is initiated from an AUU codon upstream of and in-frame with the AUG initiation sequence for canonical PTEN. We show that the Kozak context and a downstream hairpin structure are critical for this alternative initiation. PTENβ localizes predominantly in the nucleolus, and physically associates with and dephosphorylates nucleolin, which is a multifunctional nucleolar phosphoprotein. Disruption of PTENβ alters rDNA transcription and promotes ribosomal biogenesis, and this effect can be reversed by re-introduction of PTENβ. Our data show that PTENβ regulates pre-rRNA synthesis and cellular proliferation. These results demonstrate the complexity of the PTEN protein family and the diversity of its functions.

High association of T1858-G1896 precore mutations with impaired base pairing and high hepatitis B virus DNA levels in HBeAg-negative chronically infected patients

The progression of liver disease in hepatitis B virus (HBV) infection is fostered by active virus replication. Mutations in the basal core promoter (BCP) and precore (PC) regions of the HBV genome are known to have an impact on viral replication. The aim of the present study was to assess the correlation of mutation profiles in the BCP and PC regions with the viral load in HBeAg-negative chronically infected patients. The HBV genotype, BCP/PC mutations, serum HBV DNA levels, and associated serological markers were analyzed in 92 HBeAg-negative chronically infected patients. Sequence analysis of the BCP and PC regions revealed variability of 19% and 24.1%, respectively. This variability was primarily associated with five critical positions (1753, 1762, 1764, 1896 and 1899). An elevated HBV viral load (>20,000 IU/ml) was classically correlated with F2-F4 liver fibrosis, elevated serum alanine aminotransferase levels, 1762/1764 and 1753 combination mutations, and surprisingly, with an 1858T-1896G double mutation that impairs base pairing at the base of the bulge in the ε encapsidation signal. An analysis of covariance confirmed the independent nature of the relationship between the 1858T-1896G double mutation and the HBV viral load. In conclusion, independently of conventional parameters, this study demonstrates that a high serum HBV DNA level was also associated with PC 1858-1896 mutations. These BCP/PC mutations may have important clinical implications as predictive factors for HBV DNA increase.

Identification of cvSI-3 and evidence for the wide distribution and active evolution of the I84 family of protease inhibitors in mollusks

Protease inhibitors are an extremely diverse group of proteins that control the proteolytic activities of proteases and play a crucial role in biological processes including host defenses. The I84 family of protease inhibitors in the MEROPS database currently consists of cvSI-1 and cvSI-2, two novel serine protease inhibitors purified and characterized from the eastern oyster Crassostrea virginica plasma and believed to play a role in host defense and disease resistance. In the present study, a third member of I84 family, named cvSI-3, was identified from C. virginica by cDNA cloning and sequencing. The full cvSI-3 cDNA was composed of 342 bp including a 255 bp open reading frame (ORF) that encodes an 84-amino acid peptide. The mature cvSI-3 molecule was predicted to have 68 amino acid residues after removal of a 16-amino acid signal peptide, with a calculated molecular mass of 7724.5 Da and a theoretical isoelectric point (pI) of 6.28. CvSI-3 amino acid sequence shared 41% identity with cvSI-2 and 37% identity with cvSI-1, which included 12 conserved cysteines. Quantitative real-time PCR determined that cvSI-3 gene expressed primarily in oyster digestive glands. Real-time PCR also detected that cvSI-1, cvSI-2 and cvSI-3 expression levels in digestive glands varied significantly, with cvSI-2 showing the highest expression level and cvSI-3 the lowest. Additionally, a significant correlation was detected between cvSI-2 and cvSI-3 mRNAs levels. Searches into sequence databases using cvSI-1, cvSI-2 and cvSI-3 as queries retrieved ESTs suggesting the possible existence of at least 9 more I84 family members in eastern oysters and of I84 family protease inhibitors in various bivalve and gastropod species. Moreover, orthologs of all C. virginica I84 family members or potential member genes were found to be present in the C. gigas genome, and their distributions among species provided important information about the evolution of the I84 family of protease inhibitors. It appears that the I84 family of protease inhibitors is widely distributed and actively evolving in the Phylum Mollusca.

The ins and outs of eukaryotic viruses: Knowledge base and ontology of a viral infection

Viruses are genetically diverse, infect a wide range of tissues and host cells and follow unique processes for replicating themselves. All these processes were investigated and indexed in ViralZone knowledge base. To facilitate standardizing data, a simple ontology of viral life-cycle terms was developed to provide a common vocabulary for annotating data sets. New terminology was developed to address unique viral replication cycle processes, and existing terminology was modified and adapted. The virus life-cycle is classically described by schematic pictures. Using this ontology, it can be represented by a combination of successive terms: “entry”, “latency”, “transcription”, “replication” and “exit”. Each of these parts is broken down into discrete steps. For example Zika virus “entry” is broken down in successive steps: “Attachment”, “Apoptotic mimicry”, “Viral endocytosis/ macropinocytosis”, “Fusion with host endosomal membrane”, “Viral factory”. To demonstrate the utility of a standard ontology for virus biology, this work was completed by annotating virus data in the ViralZone, UniProtKB and Gene Ontology databases.

Transductive learning as an alternative to translation initiation site identification

Background

The correct protein coding region identification is an important and latent problem in the molecular biology field. This problem becomes a challenge due to the lack of deep knowledge about the biological systems and unfamiliarity of conservative characteristics in the messenger RNA (mRNA). Therefore, it is fundamental to research for computational methods aiming to help the patterns discovery for identification of the Translation Initiation Sites (TIS). In the field of Bioinformatics, machine learning methods have been widely applied based on the inductive inference, as Inductive Support Vector Machine (ISVM). On the other hand, not so much attention has been given to transductive inference-based machine learning methods such as Transductive Support Vector Machine (TSVM). The transductive inference performs well for problems in which the amount of unlabeled sequences is considerably greater than the labeled ones. Similarly, the problem of predicting the TIS may take advantage of transductive methods due to the fact that the amount of new sequences grows rapidly with the progress of Genome Project that allows the study of new organisms. Consequently, this work aims to investigate the transductive learning towards TIS identification and compare the results with those obtained in inductive method.

Results

The transductive inference presents better results both in F-measure and in sensitivity in comparison with the inductive method for predicting the TIS. Additionally, it presents the least failure rate for identifying the TIS, presenting a smaller number of False Negatives (FN) than the ISVM. The ISVM and TSVM methods were validated with the molecules from the most representative organisms contained in the RefSeq database: Rattus norvegicus, Mus musculus, Homo sapiens, Drosophila melanogaster and Arabidopsis thaliana. The transductive method presented F-measure and sensitivity higher than 90% and also higher than the results obtained with ISVM. The ISVM and TSVM approaches were implemented in the TransduTIS tool, TransduTIS-I and TransduTIS-T respectively, available in a web interface. These approaches were compared with the TISHunter, TIS Miner, NetStart tools, presenting satisfactory results.

Conclusions

In relation to precision, the results are similar for the ISVM and TSVM classifiers. However, the results show that the application of TSVM approach ensured an improvement, specially for F-measure and sensitivity. Moreover, it was possible to identify a potential for the application of TSVM, which is for organisms in the initial study phase with few identified sequences in the databases.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1502-6) contains supplementary material, which is available to authorized users.

High level expression of recombinant human growth hormone in Escherichia coli: crucial role of translation initiation region

For high-throughput production of recombinant protein in Escherichia coli (E. coli), besides important parameters such as efficient vector with strong promoter and compatible host, other important issues including codon usage, rare codons, and GC content specially at N-terminal region should be considered. In the current study, the effect of decreasing the percentage of GC nucleotides and optimizing codon usage at N-terminal region of human growth hormone (hGH) cDNA on the level of its expression in E. coli were investigated. Mutation in cDNA of hGH was performed through site-directed mutagenesis using PCR. Then, the mutant genes were amplified and cloned into the expression vector, pET-28a. The new constructs were transformed into the BL21(DE3) strain of E. coli and chemically induced for hGH expression. At the final stage, expressed proteins were analyzed by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE), scanning gel densitometry, and western blot. SDS-PAGE scanning gel densitometry assay and western blot analysis revealed higher expression level of hGH by using the two new expressions constructs (mutant genes vectors with decreasing GC content and optimized-codon usage at N-terminal of cDNA) in comparison with wild gene expression vector. Obtained results demonstrated that decreasing the GC nucleotide content and optimization of codon usage at N-terminal of the hGH cDNA could significantly enhance the expression of the target protein in E. coli. Our results highlight the important role of both 5´ region of the heterologous genes in terms of codon usage and also GC content on non-host protein expression in E. coli.

Long Noncoding RNA Identification: Comparing Machine Learning Based Tools for Long Noncoding Transcripts Discrimination

Long noncoding RNA (lncRNA) is a kind of noncoding RNA with length more than 200 nucleotides, which aroused interest of people in recent years. Lots of studies have confirmed that human genome contains many thousands of lncRNAs which exert great influence over some critical regulators of cellular process. With the advent of high-throughput sequencing technologies, a great quantity of sequences is waiting for exploitation. Thus, many programs are developed to distinguish differences between coding and long noncoding transcripts. Different programs are generally designed to be utilised under different circumstances and it is sensible and practical to select an appropriate method according to a certain situation. In this review, several popular methods and their advantages, disadvantages, and application scopes are summarised to assist people in employing a suitable method and obtaining a more reliable result.

Engineering Translation in Mammalian Cell Factories to Increase Protein Yield: The Unexpected Use of Long Non-Coding SINEUP RNAs

Mammalian cells are an indispensable tool for the production of recombinant proteins in contexts where function depends on post-translational modifications. Among them, Chinese Hamster Ovary (CHO) cells are the primary factories for the production of therapeutic proteins, including monoclonal antibodies (MAbs). To improve expression and stability, several methodologies have been adopted, including methods based on media formulation, selective pressure and cell- or vector engineering. This review presents current approaches aimed at improving mammalian cell factories that are based on the enhancement of translation. Among well-established techniques (codon optimization and improvement of mRNA secondary structure), we describe SINEUPs, a family of antisense long non-coding RNAs that are able to increase translation of partially overlapping protein-coding mRNAs. By exploiting their modular structure, SINEUP molecules can be designed to target virtually any mRNA of interest, and thus to increase the production of secreted proteins. Thus, synthetic SINEUPs represent a new versatile tool to improve the production of secreted proteins in biomanufacturing processes.

Nucleotide diversity analysis highlights functionally important genomic regions

We analyzed functionality and relative distribution of genetic variants across the complete Oryza sativa genome, using the 40 million single nucleotide polymorphisms (SNPs) dataset from the 3,000 Rice Genomes Project (http://snp-seek.irri.org), the largest and highest density SNP collection for any higher plant. We have shown that the DNA-binding transcription factors (TFs) are the most conserved group of genes, whereas kinases and membrane-localized transporters are the most variable ones. TFs may be conserved because they belong to some of the most connected regulatory hubs that modulate transcription of vast downstream gene networks, whereas signaling kinases and transporters need to adapt rapidly to changing environmental conditions. In general, the observed profound patterns of nucleotide variability reveal functionally important genomic regions. As expected, nucleotide diversity is much higher in intergenic regions than within gene bodies (regions spanning gene models), and protein-coding sequences are more conserved than untranslated gene regions. We have observed a sharp decline in nucleotide diversity that begins at about 250 nucleotides upstream of the transcription start and reaches minimal diversity exactly at the transcription start. We found the transcription termination sites to have remarkably symmetrical patterns of SNP density, implying presence of functional sites near transcription termination. Also, nucleotide diversity was significantly lower near 3′ UTRs, the area rich with regulatory regions.

A Functional Polymorphism (rs10817938) in the XPA Promoter Region Is Associated with Poor Prognosis of Oral Squamous Cell Carcinoma in a Chinese Han Population

Single nucleotide polymorphisms of XPA gene have been studied in several cancers such as rs10817938, rs2808668. However, the role of XPA polymorphisms in patients with oral squamous cell carcinoma (OSCC) remains unclear. Thus, we analyzed the association of XPA polymorphisms with OSCC risk, clinicopathological characteristics and prognosis in the present study. TaqMan genotyping was used to evaluate the frequency of rs10817938, rs2808668 polymorphisms in OSCC patients. The prognostic significance of these polymorphisms was evaluated using Kaplan-Meier curves, Log-Rank analyses, and the Cox proportional hazard model. Luciferase reporter assay, RT-PCR and western blot were used to determine whether rs10817938 could influence transcription activity and XPA expression. The results showed that individuals carrying TC and CC genotypes had significantly greater risk of developing OSCC (OR = 1.42, 95% CI 1.04-1.93; OR = 2.75, 95% CI 1.32-5.71, respectively) when compared with wild-type TT genotype at rs10817938. OSCC patients with C allele at rs10817938 were more susceptible to lymph metastases, poor pathological differentiation and late TNM stage (OR = 1.67, 95% CI 1.17-2.37; OR = 1.64, 95% CI 1.18-2.28; OR = 1.54, 95% CI 1.11-2.14; respectively). A significant gene-environment interaction between smoking and CC genotype at rs10817938 was observed (COR = 3.60, 95% CI 1.20-10.9) and data also showed that OSCC patients with CC genotype and C allele had worse survival (p<0.001 for both). The T to C substitution at rs10817938 significantly decreased transcription activity of XPA gene, XPA mRNA and protein were also decreased in individuals with C allele at rs10817938. In addition, no significant association of rs2808668 polymorphism with OSCC risk, prognosis could be observed. In conclusion, the present study showed that XPA rs10817938 polymorphism is a functional SNP in vitro and in vivo and a biomarker for poor prognosis in OSCC patients.

A novel VHLα isoform inhibits Warburg effect via modulation of PKM splicing

Von Hippel-Lindau (VHL) is the most frequently mutated gene in clear cell renal carcinoma. Here, we identified a novel translational variant of VHL, termed VHLα, initiated from an alternative translational start site upstream and in frame with the ATG start codon. We showed that VHLα interacts with and regulates heterogeneous nuclear ribonucleoprotein A2B1 (hnRNPA2B1), which consequently modulates pyruvate kinase transcript splicing and reprograms cellular glucose metabolism. Our study demonstrated that a novel VHL isoform may function as a tumor suppressor through inhibiting the Warburg effect.

Regulation of Rev expression by the equine infectious anaemia virus tat-rev mRNA Kozak sequence and its potential influence on viral replication

Rev, an important accessory protein of equine infectious anaemia virus (EIAV), induces the nuclear export of incompletely spliced viral mRNAs. Rev is translated from the tat-rev mRNA through leaky scanning of the tat CUG. In this study, the function of the Kozak sequence at the beginning of the rev ORF was investigated. Deletion or attenuation of the Kozak sequence resulted in expression of an N-terminal 11  aa-truncated Rev in addition to WT Rev. Truncated Rev displayed weaker promotion of Gag expression and processing than WT Rev. Furthermore, EIAV rescued from an infectious molecular clone (pEIAVUK3) with Kozak attenuation exhibited decreased viral replication in host cells in vitro. These results provide a new understanding of the relationship between EIAV Rev expression and viral replication.