Codon usage tabulated from the GenBank genetic sequence data

Prediction of gene expression under drought stress in spring wheat using codon usage pattern

Spring wheat (Triticum aestivum) is a staple food providing sources of essential proteins for human. In fact, gene expressions of wheat play an important role in growth and productivity that are affected by drought stress. The objective of this work focused on analysis gene feature on spring wheat represented by nucleotide and gene expressions under drought stress. It was found that the higher codon adaptation index was in both wheat root and L-galactono-1, 4-lactone dehydrogenase. It was also found that guanine and cytosine content were high (55.56%) in wheat root. Whereas, guanine and cytosine content were low (41.28%) in L-galactono-1, 4-lactone dehydrogenase. Moreover, the higher relative synonymous codon usage value was observed in codon CAA (1.20), GAA (1.33), GAT (1.00), and ATG (1.00) in wheat root and thus about 62.95% of the total variation in relative synonymous codon was explained by principal component analysis. Additionally, high averages frequency number of codon were (above 15.76) in Met, Lys, Ala, Gly, Phe, Asp, Glu, His, and Tyr; whereas, low averages were in remaining amino acids and majority (90%) of modified relative codon bias values was between 0.40 and 0.90. Shortly, calculations and analysis of codon usage pattern under drought stress would help for genetic engineering, molecular evolution, and gene prediction in wheat studies for developing varieties that associate with drought tolerance.

Investigating the role of site specific synonymous variation in disease association studies

Synonymous codon changes may not always be neutral indicating their significance in disease association studies, which is almost always overlooked. Synonymous substitutions may affect protein-folding rates leading to protein misfolding and aggregation. Genome wide analysis of 2301 mitochondrial genomes is performed to evaluate the significance of synonymous codons in disease association studies. The analysis revealed usage of rare codons at several sites in mitochondrial genes with rare codon usage higher for hydrophobic amino acids. The analysis suggests that variation data in association studies should be analyzed using site-specific codon usage values to infer the potential phenotypic impact of synonymous changes.

Biologically inspired intelligent decision making: a commentary on the use of artificial neural networks in bioinformatics

Artificial neural networks (ANNs) are a class of powerful machine learning models for classification and function approximation which have analogs in nature. An ANN learns to map stimuli to responses through repeated evaluation of exemplars of the mapping. This learning approach results in networks which are recognized for their noise tolerance and ability to generalize meaningful responses for novel stimuli. It is these properties of ANNs which make them appealing for applications to bioinformatics problems where interpretation of data may not always be obvious, and where the domain knowledge required for deductive techniques is incomplete or can cause a combinatorial explosion of rules. In this paper, we provide an introduction to artificial neural network theory and review some interesting recent applications to bioinformatics problems.

Nucleotide and phylogenetic analysis of human papillomavirus types 6 and 11 isolated from recurrent respiratory papillomatosis in Brazil

There are few studies about the distribution of natural molecular variants of low-risk HPVs. Our aim was to evaluate the E6 early gene variability among HPV-6 and HPV-11 isolates detected in recurrent respiratory papillomatosis (RRP) samples obtained in a cohort of Brazilian patients. We also performed a phylogenetic analysis in order to compare nucleotide sequences identified in our study with previously reported isolates from different anatomic sites (laryngeal papillomas, genital warts, cervical cancer and anal swabs) obtained from other parts of the world to determine the phylogenetic relationships of variants detected in Brazil. The complete coding region of the E6 gene of 25 samples was cloned and sequenced: 18 isolates of HPV-6 (72%) and 7 isolates of HPV-11 (28%). A total of four different HPV-6 genomic variants and two HPV-11 genomic variants was identified. It was not possible to correlate specific variants with disease severity. Phylogenetic trees for both HPV types were constructed enclosing both E6 sequences detected in our study and formerly published sequences. In both phylogenetic trees, the sequences from Brazil did not group together. We could not establish a geographical association between HPV-6 or HPV-11 variants, unlike HPV-16 and HPV-18.

Efficient production of HIV-1 virus-like particles from a mammalian expression vector requires the N-terminal capsid domain

It is now well accepted that the structural protein Pr55(Gag) is sufficient by itself to produce HIV-1 virus-like particles (VLPs). This polyprotein precursor contains different domains including matrix, capsid, SP1, nucleocapsid, SP2 and p6. In the present study, we wanted to determine by mutagenesis which region(s) is essential to the production of VLPs when Pr55(Gag) is inserted in a mammalian expression vector, which allows studying the protein of interest in the absence of other viral proteins. To do so, we first studied a minimal Pr55(Gag) sequence called Gag min that was used previously. We found that Gag min fails to produce VLPs when expressed in an expression vector instead of within a molecular clone. This failure occurs early in the cell at the assembly of viral proteins. We then generated a series of deletion and substitution mutants, and examined their ability to produce VLPs by combining biochemical and microscopic approaches. We demonstrate that the matrix region is not necessary, but that the efficiency of VLP production depends strongly on the presence of its basic region. Moreover, the presence of the N-terminal domain of capsid is required for VLP production when Gag is expressed alone. These findings, combined with previous observations indicating that HIV-1 Pr55(Gag)-derived VLPs act as potent stimulators of innate and acquired immunity, make the use of this strategy worth considering for vaccine development.

Codon usage roles in human papillomavirus

Human papillomavirus (HPV) genomes, similar to other virus genomes, frequently have a G + C content significantly different from their host species. The HPV genomes show a strong codon usage bias to 18 codons, with 14 showing T at the third position amongst degenerately encoded amino acids. The codon usage pattern in HPV genome plays an important role, which regulates low or non-translational expression of the viral capsid genes and results in very weak protein expression of oncogenes in a wide range of mammalian cells. Codon modification has been proved to be a powerful technology to overcome the translational blockage and weak expression of both HPV capsid genes and oncogenes in different expression systems. Furthermore, keratinocytes are the host cells of HPV infection; the codon usage in HPV capsid genes matches available aminoacyl-tRNAs in differentiated keratinocytes to modulate their protein expression. HPV DNA vaccines with codon optimization have been shown to have higher immunogenicity and induce both strong cellular and humoral responses in animal models, which may be a promising form of therapeutic HPV vaccines.

Computational identification of rare codons of Escherichia coli based on codon pairs preference

Background

Codon bias is believed to play an important role in the control of gene expression. In Escherichia coli, some rare codons, which can limit the expression level of exogenous protein, have been defined by gene engineering operations. Previous studies have confirmed the existence of codon pair's preference in many genomes, but the underlying cause of this bias has not been well established. Here we focus on the patterns of rarely-used synonymous codons. A novel method was introduced to identify the rare codons merely by codon pair bias in Escherichia coli.

Results

In Escherichia coli, we defined the "rare codon pairs" by calculating the frequency of occurrence of all codon pairs in coding sequences. Rare codons which are disliked in genes could make great contributions to forming rare codon pairs. Meanwhile our investigation showed that many of these rare codon pairs contain termination codons and the recognized sites of restriction enzymes. Furthermore, a new index (F_rare) was developed. Through comparison with the classical indices we found a significant negative correlation between F_rare and the indices which depend on reference datasets.

Conclusions

Our approach suggests that we can identify rare codons by studying the context in which a codon lies. Also, the frequency of rare codons (F_rare) could be a useful index of codon bias regardless of the lack of expression abundance information.

Preclinical evaluation of the immunogenicity of C-type HIV-1-based DNA and NYVAC vaccines in the Balb/C mouse model

As part of a European initiative (EuroVacc), we report the design, construction, and immunogenicity of two HIV-1 vaccine candidates based on a clade C virus strain (CN54) representing the current major epidemic in Asia and parts of Africa. Open reading frames encoding an artificial 160-kDa GagPolNef (GPN) polyprotein and the external glycoprotein gp120 were fully RNA and codon optimized. A DNA vaccine (DNA-GPN and DNA-gp120, referred to as DNA-C), and a replication-deficient vaccinia virus encoding both reading frames (NYVAC-C), were assessed regarding immunogenicity in Balb/C mice. The intramuscular administration of both plasmid DNA constructs, followed by two booster DNA immunizations, induced substantial T-cell responses against both antigens as well as Env-specific antibodies. Whereas low doses of NYVAC-C failed to induce specific CTL or antibodies, high doses generated cellular as well as humoral immune responses, but these did not reach the levels seen following DNA vaccination. The most potent immune responses were detectable using prime:boost protocols, regardless of whether DNA-C or NYVAC-C was used as the priming or boosting agent. These preclinical findings revealed the immunogenic response triggered by DNA-C and its enhancement by combining it with NYVAC-C, thus complementing the macaque preclinical and human phase I clinical studies of EuroVacc.

The complete nucleotide sequence of region 1 of the CFA/I fimbria! operon of human enterotoxigenicEscherichia coli

The production of the plasmid-encoded fimbrial antigen CFA/I of enterotoxigenic Escherichia coli requires two DNA regions: CFA/I region 1 and CFA/I region 2. These two regions are separated by about 40 kb on the wildtype plasmid. CFA/I region 1 contains the structural genes, whereas CFA/I region 2 contains a positive regulator. The first two genes (cfaA and cfaB) and the cfaD' sequence of region 1 have already been described. Here the total nucleotide sequence of region 1 is presented. Two new genes in region 1 are described, named cfaC and cfaE. The GC content of the genes in region 1 is 33.6% which is substantially lower than normally found in E. coli genes (50%). The codon usage also differs from the standard codons used in E. coli.

Better protective effects in rhesus macaques by combining systemic and mucosal application of a dual component vector vaccine after rectal SHIV89.6P challenge compared to systemic vaccination alone

In this study we investigated the efficacy of a multigenic DNA prime/modified vaccinia Ankara (MVA)boost vaccine approach, followed by mucosal challenge with highly pathogenic simian-human immunodeficiency virus (SHIV) 89.6P, using different routes for vaccine delivery. After three times of DNA priming (SIVmac239, GagPol, and SHIV 89.6P Env) one vaccine group of monkeys was immunized with MVA systemically via intramuscular (IM) and intradermal (ID) application, and in another vaccine group the MVA booster immunization comprised the IM, ID, and atraumatic oral route. Although all vaccinees became infected after intra-rectal challenge with SHIV 89.6P, substantial protection as indicated by lower peak and set point viral loads and unambiguous preservation of CD4 T cells could be achieved. As we could only transiently detect low levels of neutralizing antibodies in some vaccinees, these antibodies did not seem to add to the protection in the vaccinees. Our results indicate that both preventive multigenic DNA prime/MVA booster immunization strategies promote the control of virus replication and protect from disease progression. We also demonstrated that combining mucosal and systemic vaccination mediated better protective effects compared to systemic vaccination alone.

Differential CD4+ versus CD8+ T-cell responses elicited by different poxvirus-based human immunodeficiency virus type 1 vaccine candidates provide comparable efficacies in primates

Poxvirus vectors have proven to be highly effective for boosting immune responses in diverse vaccine settings. Recent reports reveal marked differences in the gene expression of human dendritic cells infected with two leading poxvirus-based human immunodeficiency virus (HIV) vaccine candidates, New York vaccinia virus (NYVAC) and modified vaccinia virus Ankara (MVA). To understand how complex genomic changes in these two vaccine vectors translate into antigen-specific systemic immune responses, we undertook a head-to-head vaccine immunogenicity and efficacy study in the pathogenic HIV type 1 (HIV-1) model of AIDS in Indian rhesus macaques. Differences in the immune responses in outbred animals were not distinguished by enzyme-linked immunospot assays, but differences were distinguished by multiparameter fluorescence-activated cell sorter analysis, revealing a difference between the number of animals with both CD4(+) and CD8(+) T-cell responses to vaccine inserts (MVA) and those that elicit a dominant CD4(+) T-cell response (NYVAC). Remarkably, vector-induced differences in CD4(+)/CD8(+) T-cell immune responses persisted for more than a year after challenge and even accompanied antigenic modulation throughout the control of chronic infection. Importantly, strong preexposure HIV-1/simian immunodeficiency virus-specific CD4(+) T-cell responses did not prove deleterious with respect to accelerated disease progression. In contrast, in this setting, animals with strong vaccine-induced polyfunctional CD4(+) T-cell responses showed efficacies similar to those with stronger CD8(+) T-cell responses.

Comparative analysis of predicted gene expression among deep-sea genomes

Deep-sea species live in an environment that is specifically characterized by extreme temperature and hydrostatic pressure. In this work, predicted highly expressed (PHX) genes are comparatively analyzed for six deep-sea microbes, which allows us to pinpoint the common highly expressed genes shared by them. The relationships between gene expression level and some basic properties such as genomic G + C content, optimal growth temperature (OGT), and environmental hydrostatic pressure of the six deep-sea species are also investigated. We find that the percentage of PHX genes out of a whole genome positively correlates to OGT for the deep-sea genomes, whereas such positive correlation seems not to exist between environmental hydrostatic pressure and percentage of PHX genes. Moreover, there exists a negative correlation between genomic G + C content and diversity of gene expression level for the deep-sea genomes, which is in sharp contrast to land-living microbes. We report the top 20 PHX genes for the six deep-sea genomes and find no common highly expressed genes shared by them except for ribosomal proteins, transcription factors, and translation factors. Our present work proffers a paradigm for studying the relationship between environmental factors and microbes' predicted gene expression level.

Gene expression level shapes the amino acid usages in Prochlorococcus marinus MED4

Prochlorococcus species are the first example of free-living bacteria with reduced genome. Codon and amino acid usages bias of Prochlorococcus marinus MED4 was investigated using all protein coding genes having length greater than or equal to 100 amino acids. Correspondence analysis on relative synonymous codon usage (RSCU) values shows that there is no such influence of translational selection in shaping the codon usage variation among the genes in this organism. However, amino acid usages were markedly different between the highly and lowly expressed genes in this organism and in particular, GC rich amino acids were found to occur significantly higher in highly expressed genes than the lowly expressed genes. Comparative analysis of the homologous genes of Synechococcus sp. WH8102 and Prochlorococcus marinus MED4 shows that amino acids conservation in highly expressed genes is significantly higher than lowly expressed genes. Based on our results we concluded that conservation of GC rich amino acids in the highly expressed genes to its ancestor is the major source of variation in amino acid usages in the organism.

The Most Deviated Codon Position in AT-rich Bacterial Genomes: A Function Related Analysis

We have performed systematic study on more than 120 archaeal and bacterial genomes. Based on the index proposed in the current paper, clear patterns are observed showing the relation between the base compositional deviation at three codon positions and the genomic GC content. For AT-rich genomes, the Most Deviated Codon Position (MDCP) is the 1st codon position, while for GC-rich genomes, MDCP appears at the 2nd or 3rd codon position alternatively. According to MDCP, the CDSs of a genome can be classified into two types: typical and atypical. In AT-rich genomes the typical represent the majority and account for about 3/4 of all the CDSs. Based on the functional classification of COG database, the two types of CDSs are examined. An apparent bias of distribution is observed that the CDSs with the function of 'information processing' are more likely to present in typical type.

Generation of a packaging cell line for prolonged large-scale production of high-titer HIV-1-based lentiviral vector

BACKGROUND

A stable packaging cell line facilitates large-scale lentivirus vector manufacture. However, it has been difficult to produce clinical-scale HIV-1-based lentiviral vectors using a packaging cell line, in part due to toxicity of packaging genes, and gene silencing that occurs during the long culture period necessary for sequential addition of packaging constructs.

METHODS

To avoid these problems, we developed a three-level cascade gene regulation system designed to remove tetracycline transactivator (tTA) from cytomegalovirus immediate early promoter (CMV)-controlled expression to reduce cytotoxicity from constitutive expression of tTA and leaky expression of packaging genes. We also performed a one-step integration of the three packaging plasmids to shorten the culture time for clonal selection.

RESULTS

Although leaky expression of p24 and vector production still occurred despite the three-level regulation system, little cytotoxicity was observed and producer cells could be expanded for large-scale production. Producer cells yielded remarkably stable vector production over a period greater than 11 days with the highest titer 3.5 x 10(7) transducing units (TU)/ml and p24 300 ng/ml, yielding 2.2 x 10(11) TU and 1.8 milligram (mg) p24 from one cell factory. No replication-competent lentivirus (RCL) was detected. Long-term analysis demonstrated that, although the cells are genetically stable, partial gene silencing occurs after 2-3 months in culture; however, the one-step construct integration allowed prolonged vector production before significant gene silencing. Concentrated vector resulted in 90% transduction in CD4+ lymphocytes at 20 TU per cell. CD34+ progenitor cells were transduced at 41-46% efficiency, and long-term initiating culture (LTC-IC) was transduced at 45-51%.

CONCLUSIONS

These results demonstrate for the first time HIV-1-based lentiviral vector production on the large scale using a packaging cell line.

Cloning and characterization of a phospholipase C-beta isoform from the sea urchin Lytechinus pictus

Calcium is a ubiquitous intracellular signaling molecule controlling a wide array of cellular processes including fertilization and egg activation. The mechanism for triggering intracellular Ca(2+) release in sea urchin eggs during fertilization is the generation of inositol-1,4,5-trisphosphate by phospholipase C (PLC) hydrolysis of phosphatidylinositol-4,5-bisphosphate. Of the five PLC isoforms identified in mammals (beta, gamma, delta, epsilon and zeta), only PLCgamma and PLCdelta have been detected in echinoderms. Here, we provide direct evidence of the presence of a PLCbeta isoform, named suPLCbeta, within sea urchin eggs. The coding sequence was cloned from eggs of Lytechinus pictus and determined to have the greatest degree of homology and identity with the mammalian PLCbeta4. The presence of suPLCbeta within the egg was verified using a specifically generated antibody. The majority of the enzyme is localized in the non-soluble fraction, presumably the plasma membrane of the unfertilized egg. This distribution remains unchanged 1 min postfertilization. Unlike PLCbeta4, suPLCbeta is activated by G protein betagamma subunits, and this activity is Ca(2+)-dependent. In contrast to all known PLCbeta enzymes, suPLCbeta is not activated by Galphaq-GTPgammaS subunit suggesting other protein regulators may be present in sea urchin eggs.

The gusBC genes of Escherichia coli encode a glucuronide transport system

Two genes, gusB and gusC, from a natural fecal isolate of Escherichia coli are shown to encode proteins responsible for transport of beta-glucuronides with synthetic [(14)C]phenyl-1-thio-beta-d-glucuronide as the substrate. These genes are located in the gus operon downstream of the gusA gene on the E. coli genome, and their expression is induced by a variety of beta-d-glucuronides. Measurements of transport in right-side-out subcellular vesicles show the system has the characteristics of secondary active transport energized by the respiration-generated proton motive force. When the genes were cloned together downstream of the tac operator-promoter in the plasmid pTTQ18 expression vector, transport activity was increased considerably with isopropylthiogalactopyranoside as the inducer. Amplified expression of the GusB and GusC proteins enabled visualization and identification by N-terminal sequencing of both proteins, which migrated at ca. 32 kDa and 44 kDa, respectively. Separate expression of the GusB protein showed that it is essential for glucuronide transport and is located in the inner membrane, while the GusC protein does not catalyze transport but assists in an as yet unknown manner and is located in the outer membrane. The output of glucuronides as waste by mammals and uptake for nutrition by gut bacteria or reabsorption by the mammalian host is discussed.

Evolutionary forces in shaping the codon and amino acid usages in Blochmannia floridanus

Endosymbiotic relationship has great effect on ecological system. Codon and amino acid usages bias of endosymbiotic bacteria Blochmannia floridanus (whose host is an ant Camponotus floridanus) was investigated using experimentally known genes of this organism. Correspondence Analysis on RSCU values show that there exists only one single explanatory major axis that is linked to the strand specific mutational biases. Majority of the genes have a tendency to concentrate on the leading strand, which may be related to the adaptive property related to the replication mechanisms. Amino acid usages were markedly different between the highly and lowly expressed genes in this organism and in particular, GC rich amino acids were found to occur significantly higher in highly expressed genes than the lowly expressed genes. Comparative analyses of the orthologous genes of Escherichia coli and Blochmannia floridanus show that highly expressed genes are significantly more conserved than lowly expressed genes. Based on our results we concluded that strand specific mutational bias is strongly operational in selecting the codon usage in this organism. Replicational-transcriptional selection can be invoked from the presence of majority of highly expressed genes in the leading strand. Conservation of GC rich amino acids in the highly expressed genes to its ancestor is the major source of variation in amino acid usages in the organism. Hydrophobicity of the genes is the second major source in differentiating the genes according to their amino acid usages in this organism.

Codon bias signatures, organization of microorganisms in codon space, and lifestyle

New and simple numerical criteria based on a codon adaptation index are applied to the complete genomic sequences of 80 Eubacteria and 16 Archaea, to infer weak and strong genome tendencies toward content bias, translational bias, and strand bias. These criteria can be applied to all microbial genomes, even those for which little biological information is known, and a codon bias signature, that is the collection of strong biases displayed by a genome, can be automatically derived. A codon bias space, where genomes are identified by their preferred codons, is proposed as a novel formal framework to interpret genomic relationships. Principal component analysis confirms that although GC content has a dominant effect on codon bias space, thermophilic and mesophilic species can be identified and separated by codon preferences. Two more examples concerning lifestyle are studied with linear discriminant analysis: suitable separating functions characterized by sets of preferred codons are provided to discriminate: translationally biased (hyper)thermophiles from mesophiles, and organisms with different respiratory characteristics, aerobic, anaerobic, facultative aerobic and facultative anaerobic. These results suggest that codon bias space might reflect the geometry of a prokaryotic "physiology space." Evolutionary perspectives are noted, numerical criteria and distances among organisms are validated on known cases, and various results and predictions are discussed both on methodological and biological grounds.

Synonymous Codon Usage inLactococcus lactis: Mutational Bias Versus Translational Selection

In this study codon usage bias of all experimentally known genes of Lactococcus lactis has been analyzed. Since Lactococcus lactis is an AT rich organism, it is expected to occur A and/or T at the third position of codons and detailed analysis of overall codon usage data indicates that A and/or T ending codons are predominant in this organism. However, multivariate statistical analyses based both on codon count and on relative synonymous codon usage (RSCU) detect a large number of genes, which are supposed to be highly expressed are clustered at one end of the first major axis, while majority of the putatively lowly expressed genes are clustered at the other end of the first major axis. It was observed that in the highly expressed genes C and T ending codons are significantly higher than the lowly expressed genes and also it was observed that C ending codons are predominant in the duets of highly expressed genes, whereas the T endings codons are abundant in the quartets. Abundance of C and T ending codons in the highly expressed genes suggest that, besides, compositional biases, translational selection are also operating in shaping the codon usage variation among the genes in this organism as observed in other compositionally skewed organisms. The second major axis generated by correspondence analysis on simple codon counts differentiates the genes into two distinct groups according to their hydrophobicity values, but the same analysis computed with relative synonymous codon usage values could not discriminate the genes according to the hydropathy values. This suggests that amino acid composition exerts constraints on codon usage in this organism. On the other hand the second major axis produced by correspondence analysis on RSCU values differentiates the genes into two groups according to the synonymous codon usage for cysteine residues (rarest amino acids in this organism), which is nothing but a artifactual effect induced by the RSCU values. Other factors such as length of the genes and the positions of the genes in the leading and lagging strand of replication have practically no influence in the codon usage variation among the genes in this organism.

Production of a fusion protein consisting of the enterotoxigenic Escherichia coli heat-labile toxin B subunit and a tuberculosis antigen in Arabidopsis thaliana

Transgenic plants are potentially safe and inexpensive vehicles to produce and mucosally deliver protective antigens. However, the application of this technology is limited by the poor response of the immune system to non-particulate, subunit vaccines. Co-delivery of therapeutic proteins with carrier proteins could increase the effectiveness of the antigen. This paper reports the ability of transgenic Arabidopsis thaliana plants to produce a fusion protein consisting of the B subunit of the Escherichia coli heat-labile enterotoxin and a 6 kDa tuberculosis antigen, the early secretory antigenic target ESAT-6. Both components of the fusion protein were detected using GM1-ganglioside-dependent enzyme-linked immunosorbant assay. This suggested the fusion protein retained both its native antigenicity and the ability to form pentamers.

Protein structure preference, tRNA copy number, and mRNA stem/loop content

From statistical analyses of protein sequences for humans and Escherichia coli we found that the messenger RNA segment of m-codons (for m=2 to 6) with average high tRNA copy number (TCN) (larger than approximately 10.5 for humans or approximately 1.95 for E. coli) preferably code for the alpha helix and that with low TCN (smaller than approximately 7.5 for humans or approximately 1.7 for E. coli) preferably code for coil. Between them there is an intermediate region without correlation to structure preference. For the beta strand the preference/ avoidance tendency is not obvious. All strong preference-modes of TCN for protein secondary structures have been deduced. The mutual interaction between two factors--protein secondary structural type and codon TCN--is tested by F distribution. A phenomenological model on the relation between structure preference and translational efficiency or accuracy is proposed. It is pointed out that the structure preference of codons is related to the distribution of mRNA stem/loop content in three TCN regions.

Oral immunogenicity of human papillomavirus-like particles expressed in potato

Human papillomavirus-like particles (HPV VLPs) have shown considerable promise as a parenteral vaccine for the prevention of cervical cancer and its precursor lesions. Parenteral vaccines are expensive to produce and deliver, however, and therefore are not optimal for use in resource-poor settings, where most cervical HPV disease occurs. Transgenic plants expressing recombinant vaccine immunogens offer an attractive and potentially inexpensive alternative to vaccination by injection. For example, edible plants can be grown locally and can be distributed easily without special training or equipment. To assess the feasibility of an HPV VLP-based edible vaccine, in this study we synthesized a plant codon-optimized version of the HPV type 11 (HPV11) L1 major capsid protein coding sequence and introduced it into tobacco and potato. We show that full-length L1 protein is expressed and localized in plant cell nuclei and that expression of L1 in plants is enhanced by removal of the carboxy-terminal nuclear localization signal sequence. We also show that plant-expressed L1 self-assembles into VLPs with immunological properties comparable to those of native HPV virions. Importantly, ingestion of transgenic L1 potato was associated with activation of an anti-VLP immune response in mice that was qualitatively similar to that induced by VLP parenteral administration, and this response was enhanced significantly by subsequent oral boosting with purified insect cell-derived VLPs. Thus, papillomavirus L1 protein can be expressed in transgenic plants to form immunologically functional VLPs, and ingestion of such material can activate potentially protective humoral immune responses.

Expression of the B subunit of Escherichia coli heat-labile enterotoxin as a fusion protein in transgenic tomato

Epitopes often require co-delivery with an adjuvant or targeting protein to enable recognition by the immune system. This paper reports the ability of transgenic tomato plants to express a fusion protein consisting of the B subunit of the Escherichia coli heat-labile enterotoxin (LTB) and an immunocontraceptive epitope. The fusion protein was found to assemble into pentamers, as evidenced by its ability to bind to gangliosides, and had an average expression level of 37.8 microg g(-1) in freeze-dried transgenic tissues. Processing of selected transgenic fruit resulted in a 16-fold increase in concentration of the antigen with minimal loss in detectable antigen. The species-specific nature of this epitope was shown by the inability of antibodies raised against non-target species to detect the LTB fusion protein. The immunocontraceptive ability of this vaccine will be tested in future pilot mice studies.

Analysis of nucleotide distribution in the genome of Streptomyces coelicolor A3(2) using the Z curve method

The nucleotide distribution of all 33,527 open reading frames (ORFs) (> or =300 bp) in the genome of Streptomyces coelicolor A3(2) has been analyzed using the Z curve method. Each ORF is mapped onto a point in a 9-dimensional space. To visualize the distribution of mapping points, the points are projected onto the principal plane based on principal component analysis. Consequently, the distribution pattern of the 33,527 points in the principal plane shows a flower-like shape, in which there are seven distinct regions. In addition to the central region, there are six petal-like regions around the center, one of which corresponds to 7172 coding sequences. The central region and the remaining five petal-like regions correspond to the intergenic sequences and out-of-frame non-coding ORFs, respectively. It is shown that selective pressure produces a remarkable bias of the G+C content among three codon positions, resulting in the interesting phenomenon observed. A similar phenomenon is also observed for other bacterial genomes with high genomic G+C content, such as Pseudomonas aeruginosa PA01 (G+C = 66.6%). However, for the genomes of Bacillus subtilis (G+C = 43.5%) and Clostridium perfringens (G+C = 28.6%), no similar phenomenon was observed. The finding presented here may be useful to improve the gene-finding algorithms for genomes with high G+C content. A set of supplementary materials including the plots displaying the base distribution patterns of ORFs in 12 prokaryotes is provided on the website http://tubic.tju.edu.cn/highGC/.

Sequence survey of the genome of the opportunistic microsporidian pathogen, Vittaforma corneae

The microsporidian Vittaforma corneae has been reported as a pathogen of the human stratum corneum, where it can cause keratitis, and is associated with systemic infections. In addition to this direct role as an infectious, etiologic agent of human disease, V. corneae has been used as a model organism for another microsporidian, Enterocytozoon bieneusi, a frequent and problematic pathogen of HIV-infected patients that, unlike V. corneae, is difficult to maintain and to study in vitro. Unfortunately, few molecular sequences are available for V. corneae. In this study, seventy-four genome survey sequences (GSS) were obtained from genomic DNA of spores of laboratory-cultured V. corneae. Approximately, 41 discontinuous kilobases of V. corneae were cloned and sequenced to generate these GSS. Putative identities were assigned to 44 of the V. corneae GSS based on BLASTX searches, representing 21 discrete proteins. Of these 21 deduced V. corneae proteins, only two had been reported previously from other microsporidia (until the recent report of the Encephalitozoon cuniculi genome). Two of the V. corneae proteins were of particular interest, reverse transcriptase and topoisomerase IV (parC). Since the existence of transposable elements in microsporidia is controversial, the presence of reverse transcriptase in V. corneae will contribute to resolution of this debate. The presence of topoisomerase IV was remarkable because this enzyme previously had been identified only from prokaryotes. The 74 GSS included 26.7 kilobases of unique sequences from which two statistics were generated: GC content and codon usage. The GC content of the unique GSS was 42%, lower than that of another microsporidian, E. cuniculi (48% for protein-encoding regions), and substantially higher than that predicted for a third microsporidian, Spraguea lophii (28%). A comparison using the Pearson correlation coefficient showed that codon usage in V. corneae was similar to that in the yeasts, Saccharomyces cerevisiae (r = 0.79) and Shizosaccharomyces pombe (r = 0.70), but was markedly dissimilar to E. cuniculi (r = 0.19).

Nitric oxide synthase sequences in the marine fish Stenotomus chrysops and the sea urchin Arbacia punctulata, and phylogenetic analysis of nitric oxide synthase calmodulin-binding domains

The phylogenetic distribution and structural diversity of the nitric oxide synthases (NOS) remain important and issues that are little understood. We present sequence information, as well as phylogenetic analysis, for three NOS cDNAs identified in two non-mammalian species: the vertebrate marine teleost fish Stenotomus chrysops (scup) and the invertebrate echinoderm Arbacia punctulata (sea urchin). Partial gene sequences containing the well-conserved calmodulin (CaM)-binding domain were amplified by RT-PCR. Identical 375-bp cDNAs were amplified from scup brain, heart, liver and spleen; this sequence shares 82% nucleic acid and 91% predicted amino acid identity with the corresponding region of human neuronal NOS. A 387-bp cDNA was amplified from sea urchin ovary and testes; this sequence shares 72% nucleic acid identity and 65% deduced amino acid identity with human neuronal NOS. A second cDNA of 381 bp was amplified from sea urchin ovary and it shares 66% nucleic acid and 57% deduced amino acid identity with the first sea urchin sequence. Together with earlier reports of neuronal and inducible NOS sequences in fish, these data indicate that multiple NOS isoforms exist in non-mammalian species. Phylogenetic analysis of these sequences confirms the conserved nature of NOS, particularly of the calmodulin-binding domains.

Gene expressivity is the main factor in dictating the codon usage variation among the genes in Pseudomonas aeruginosa

Codon usage biases of all DNA sequences (length greater than or equal to 300 bp) from the complete genome of Pseudomonas aeruginosa have been analyzed. As P. aeruginosa is a GC-rich organism, G and/or C are expected to predominate in their codons. Overall codon usage data analysis indicates that indeed codons ending in G and/or C are predominant in this organism. But multivariate statistical analysis indicates that there is a single major trend in the codon usage variation among the genes in this organism, which has a strong negative correlation with the expressivities of the genes. The majority of the lowly expressed genes are scattered towards the positive end of the major axis whereas the highly expressed genes are clustered towards the negative end. This is the first report where the prokaryotic organism having highly skewed base composition is dictated mainly by translational selection, though some other factors such as the lengths of the genes as well as the hydrophobicity of genes also influence the codon usage variation among the genes in this organism in a minor way.

Codon-usage based regulation of colicin K synthesis by the stress alarmone ppGpp

The molecular mechanism of the upregulation of Escherichia coli colicin K (Cka) synthesis during stress conditions was studied. Nutrient starvation experiments and the use of relA spoT mutant strains, IPTG-regulated overproduction of ppGpp and lacZ fusions revealed that the stringent response alarmone guanosine 3',5'-bispyrophosphate (ppGpp) is the main positive effector of Cka synthesis. Comparison of the amounts of protein produced (Western blotting) and specific mRNA (Northern blotting) before and after nutrient starvation demonstrated increases in Cka protein with unaltered specific mRNA levels, suggesting a post-transcriptional regulatory mechanism. Reporter (beta-galactosidase) assays using truncated cka of variable length fused to lacZ located the key regulatory region close to the 5' end of the cka mRNA. Closer analysis of this region indicated the presence of several rare codons, including the leucine-encoding codon CUA. Synonymous exchange of the rare codons with more frequently used ones abolished the regulatory effect of ppGpp. Supplementation of the strain with the plasmid CodonPlus carrying several rare tRNA genes yielded similar results, indicating that codon usage (in particular, the fifth codon for the amino acid leucine) and tRNA availability (i.e. tRNAleu) are the key elements of the regulatory function of ppGpp. We conclude that ppGpp regulates Cka synthesis via a novel post-transcriptional mechanism that is based on rare codon usage and variable cognate tRNA availability.

A screen for mutations in human homologues of mice exencephaly genes Tfap2alpha and Msx2 in patients with neural tube defects

BACKGROUND

Very little is known about the identity of genetic factors involved in the complex etiology of nonsyndromic neural tube defects (NTD). Potential susceptibility genes have emerged from the vast number of mutant mouse strains displaying NTD. Reasonable candidates are the human homologues of mice exencephaly genes Tfap2alpha and Msx2, which are expressed in the developing neural tube.

METHODS

A single-strand conformation analysis (SSCA) mutation screen of the coding sequences of TFAP2alpha and MSX2 was performed for 204 nonsyndromic NTD patients including cases of anencephaly (n = 10), encephalocele (n = 8), and spina bifida aperta, SBA (n = 183). A selected number of SBA patients was additionally tested for specific mutations in MTHFD, FRalpha, and PAX1 already shown to be related to NTD.

RESULTS

Two TFAP2alpha point mutations in individual SBA patients were silent on the amino acid level (C308C, T396T). On nucleic acid level, these mutations change evolutionary conserved codons and thus may influence mRNA processing and translation efficiency. One SBA patient displayed an exonic 9-bp deletion in MSX2 leading to a shortened and possibly less functional protein. None of these mutations was found in 222 controls. Seven polymorphisms detected in TFAP2alpha and MSX2 were equally distributed in patients and controls. Patients with combined heterozygosity of an exonic MSX2 and an intronic TFAP2alpha polymorphism were at a slightly increased risk of NTD (OR 1.71; 95% CI 0.57-5.39).

CONCLUSIONS

Although several new genetic variants were found in TFAP2 and MSX2, no statistically significant association was found between NTD cases and the new alleles or their combinations. Further studies are necessary to finally decide if these gene variants may have acted as susceptibility factors in our individual cases.

Rev-independent expression of synthetic gag-pol genes of human immunodeficiency virus type 1 and simian immunodeficiency virus: implications for the safety of lentiviral vectors

The safety of lentiviral vectors for clinical applications is still a major concern. The gag-pol expression plasmids and the lentiviral vectors used in previous studies contain homologous regions, which constitute a risk for recombination events. Synthetic gag-pol genes of human immunodeficiency virus type 1 (HIV-1) and simian immunodeficiency virus (SIV) were therefore constructed, in which the codon usage was optimized for expression in human cells without altering the amino acid sequences. The synthetic gag-pol genes allowed efficient expression of these genes in the absence of Rev and the 5' untranslated leader region. Both the HIV-1 and the SIV synthetic gag-pol expression plasmids could mediate transduction of an SIV vector into nondividing human cells with titers of about 10(6) transducing units/ml. Similar titers were obtained with a four-plasmid vector-packaging system based on HIV-1. Using a biological assay, homologous recombination events between the synthetic gag-pol expression plasmids and an SIV vector were undetectable and in comparison with a previously used gag-pol expression plasmid at least approximately 100-fold less frequent. By eliminating regions of homology and sequences involved in packaging, synthetic gag-pol genes should improve the safety profile of lentiviral vectors.

FOG acts as a repressor of red blood cell development in Xenopus

Members of the GATA family of zinc-finger transcription factors have critical roles in a variety of cell types. GATA-1, GATA-2 and GATA-3 are required for proliferation and differentiation of several hematopoietic lineages, whereas GATA-4, GATA-5 and GATA-6 activate cardiac and endoderm gene expression programs. Two GATA cofactors have recently been identified. Friend of GATA-1 (FOG-1) interacts with GATA-1 and is expressed principally in hematopoietic lineages, whereas FOG-2 is expressed predominantly in heart and brain. Although gene targeting experiments are consistent with an essential role for FOG-1 as an activator of GATA-1 function, reporter assays in transfected cells indicate that FOG-1 and FOG-2 can act as repressors. We have cloned a Xenopus laevis homologue of FOG that is structurally most similar to FOG-1, but is expressed predominantly in heart and brain, as well as the ventral blood island and adult spleen. Ectopic expression and explant assays demonstrate that FOG proteins can act as repressors in vivo, in part through interaction with the transcriptional co-repressor, C-terminal Binding Protein (CtBP). FOG may regulate the differentiation of red blood cells by modulating expression and activity of GATA-1 and GATA-2. We propose that the FOG proteins participate in the switch from progenitor proliferation to red blood cell maturation and differentiation.

The presence of the alternatively spliced A2 cassette in the vacuolar H+-ATPase subunit A prevents assembly of the V1 catalytic domain

Vacuolar ATPases (V-ATPases) are multisubunit enzymes that couple the hydrolysis of ATP to the transport of H+ across membranes, and thus acidify several intracellular compartments and some extracellular spaces. Despite the high degree of genetic and pharmacological homogeneity of V-ATPases, cells differentially modulate the lumenal pH of organelles and, in some cells, V-ATPases are selectively targetted to the plasma membrane. Although the mechanisms underlying such differences are not known, the subunit isoform composition of V-ATPases could contribute to altered assembly, targeting or activity. We previously identified an alternatively spliced variant of the chicken A subunit in which a 30 amino acid cassette (A1) containing the Walker consensus sequence for ATP binding is replaced by a 24 amino acid cassette (A2) that lacks this feature. We have examined the ability of chimeric yeast/chicken A subunits containing either the A1 or the A2 cassette to restore the V-ATPase activity of yeast that lack the A subunit. The A1-containing chimeric subunit, but not the chimera that contains the A2 cassette, partially restores the ability of the mutated yeast to grow at neutral pH. Both chimeric proteins are expressed, although at lower levels than the similarly transfected yeast A subunit. The A2-containing subunit fails to associate with the vacuolar membrane or support the assembly of V-ATPase complexes. Thus, the substitution of the A1 sequence by A2 not only removes the Walker nucleotide binding sequence but also compromises the ability of the A subunit to assemble with other V-ATPase subunits.

Identification of an intragenic ribosome binding site that affects expression of the uncB gene of the Escherichia coli proton-translocating ATPase (unc) operon

The uncB gene codes for the a subunit of the Fo proton channel sector of the Escherichia coli F1 Fo ATPase. Control of expression of uncB appears to be exerted at some step after translational initiation. Sequence analysis by the perceptron matrices (G. D. Stormo, T. D. Schneider, L. Gold, and A. Ehrenfeucht, Nucleic Acids Res. 10:2997-3011, 1982) identified a potential ribosome binding site within the uncB reading frame preceding a five-codon reading frame which is shifted one base relative to the uncB reading frame. Elimination of this binding site by mutagenesis resulted in a four- to fivefold increase in expression of an uncB'-'lacZ fusion gene containing most of uncB. Primer extension inhibition (toeprint) analysis to measure ribosome binding demonstrated that ribosomes could form an initiation complex at this alternative start site. Two fusions of lacZ to the alternative reading frame demonstrated that this site is recognized by ribosomes in vivo. The results suggest that expression of uncB is reduced by translational frameshifting and/or a translational false start at this site within the uncB reading frame.

The base contents of A, C, G or U for the three codon positions and the total coding sequences show positive correlation

The distribution of the occurrence frequencies of each of the four bases at the first, second and third codon positions and in the total coding sequences is analyzed by a graphic method. It is shown that for the coding sequences of 90 species, A has its largest frequency at the second codon position and the smallest one at the third position. C and U have their least frequencies at the first codon position, while G has its largest frequency at the first codon position. By this method, we also find that for each base, there is positive correlation between every two frequencies of the base in the first, second, third and the total coding sequences for 90 species. For each of the four bases, the correlation between the frequencies at the third codon position and that in the total coding sequences is more prominent than others. A statistical method is used to give a precise description of the correlation for the frequencies of every base and it is found that the conclusions drawn by the graphic method are consistent with that got by the statistical method.

Edible vaccine protects mice against Escherichia coli heat-labile enterotoxin (LT): potatoes expressing a synthetic LT-B gene

The authors have designed and constructed a plant-optimize synthetic gene encoding the Escherichia coli heat-labile enterotoxin B subunit (LT-B), for use in transgenic plants as an edible vaccine against enterotoxigenic E. coli. Expression of the synthetic LT-B gene in potato plants under the control of a constitutive promoter yielded increased accumulation of LT-B in leaves and tubers, as compared to the bacterial LT-B gene. The plant-derived LT-B assembled into native pentameric structures as evidenced by its ability to bind ganglioside. The authors demonstrated immunogenicity by feeding mice the raw tubers and comparing the anti-LT-B serum IgG and faecal IgA to that produced in mice gavaged with bacterial LT-B. Mice were fed three weekly doses of 5 g tuber tissue containing either 20 or 50 micrograms LT-B, or gavaged weekly with 5 micrograms of LT-B from recombinant E. coli. One week after the third dose, mice immunized with potato LT-B had higher levels of serum and mucosal anti-LT-B than those gavaged with bacterial LT-B. Mice were challenged by oral administration of 25 micrograms LT, and protection assessed by comparing the gut/carcass mass ratios. Although none of the mice were completely protected, the higher dose potato vaccine compared favourably with the bacterial vaccine. These findings show that an edible vaccine against E. coli LT-B is feasible.

Light regulation of Fed-1 mRNA requires an element in the 5' untranslated region and correlates with differential polyribosome association

Light regulation of Fed-1 mRNA abundance in the leaves of green plants is primarily a post-transcriptional process. Previously, we have shown that the Fed-1 mRNA light response requires an open reading frame, indicating that the light regulation of the mRNA depends on its concurrent translation. We now show that light-induced increases in Fed-1 mRNA abundance are associated with increases in polyribosome association that require both a functional AUG and a normal Fed-1 translational start context. We also present evidence that light regulation of Fed-1 mRNA levels requires more than efficient translation per se. Substitution of the efficiently translated tobacco mosaic virus Omega 5' untranslated region resulted in a loss of Fed-1 light regulation. In addition, we identified a CAT T repeat element located near the 5' terminus of the Fed-1 5' untranslated region that is essential for light regulation. We introduced two different mutations in the CAT T repeat element, but only one of these substitutions blocked the normal light effect on polyribosome association, whereas both altered dark-induced Fed-1 mRNA disappearance. The element may thus be important for Fed-1 mRNA stability rather than polyribosome loading. We propose a model in which Fed-1 mRNA is relatively stable when it is associated with polyribosomes in illuminated plants but in darkness is not polyribosome associated and is thus rapidly degraded by a process involving the CAT T repeat element.

Organization and transcriptional characterization of the cat1 gene cluster in Acinetobacter lwoffi K24

Previously, we have reported that two clustered cat genes from Acenitobacter lwoffi K24 had different arrangements, catB1C1A1 and catB2A2C2 (Kim, S.I., S.-H. Leem, J.-S. Choi, Y.H. Chung, S. Kim, Y.-M. Park, Y.K. Park, Y.N. Lee, and K.-S. Ha. 1997, J. Bacteriol. 179, 5226-5231). By further analysis of the organization of the cat1 gene cluster, we obtained a complete sequence of the catB1 gene, which encoded 40.8-kDa polypeptide containing 379 amino acids, and found a open reading frame (ORF) coding a putative regulatory protein in upstream region of catB1 on plasmid pCD1-1. This ORF encoded 34.2-kDa polypeptide containing 379 amino acids and had more than 40% identity with catR, LysR family regulatory protein of Pseudomonas putida. RT-PCR, Northern blot analysis and primer extension assay for transcriptional analysis of the cat1 gene cluster revealed that the catB1C1 genes were cotranscribed and the catA1 gene was independently transcribed.

An Integrated Sequence-Structure Database incorporating matching mRNA sequence, amino acid sequence and protein three-dimensional structure data

We have constructed a non-homologous database, termed the Integrated Sequence-Structure Database (ISSD) which comprises the coding sequences of genes, amino acid sequences of the corresponding proteins, their secondary structure and straight phi,psi angles assignments, and polypeptide backbone coordinates. Each protein entry in the database holds the alignment of nucleotide sequence, amino acid sequence and the PDB three-dimensional structure data. The nucleotide and amino acid sequences for each entry are selected on the basis of exact matches of the source organism and cell environment. The current version 1.0 of ISSD is available on the WWW at http://www.protein.bio.msu.su/issd/ and includes 107 non-homologous mammalian proteins, of which 80 are human proteins. The database has been used by us for the analysis of synonymous codon usage patterns in mRNA sequences showing their correlation with the three-dimensional structure features in the encoded proteins. Possible ISSD applications include optimisation of protein expression, improvement of the protein structure prediction accuracy, and analysis of evolutionary aspects of the nucleotide sequence-protein structure relationship.

Bacterial expression and purification of biologically active mouse c-Fos proteins by selective codon optimization

A simple strategy using selective codon optimization was devised to express mouse c-Fos protein in high levels in E. coli. Ten arginine codons located in the basic region were optimized to achieve high levels of protein expression. The c-Fos protein was purified to near homogeneity and was demonstrated to be biologically active by assaying its several biological activities.

Codon usage tabulated from the international DNA sequence databases

The codon usage in individual protein genes has been calculated using the nucleotide sequence obtained from the GenBank Genetic Sequence Database. Sum of the codon use of each organism has been also calculated. The data files can be obtained from anonymous ftp sites of DDBJ, DISC and EBI. The list of codon usage of genes in organisms was made searchable by name of organism through a web site. The compilation has been synchronized with a major release of GenBank.

Some molecular insights into schistosome evolution

Robust phylogenies based on molecular data for species within the genus Schistosoma have been generated in recent years. The considerable progress made in understanding the relationships between many of the 19 recognised species of Schistosoma is reviewed with particular attention being given to the detection and analysis of parasite variation as shown by studies on ribosomal RNA genes, mitochondrial DNA and RAPDs. For the most part, molecular phylogenies agree with observations based on morphological or life-history characteristics. It is clear that the parasites of man do not form a monophyletic group and that close relationships exist between parasites within species groups, especially in the S. haematobium group of species. The S. japonicum group appears to be the most divergent of the species groups and yet little DNA sequence variation has been observed between various isolates of S: japonicum. Some of the less studied schistosomes have yet to be examined at the molecular level and may prove to be interesting links between the species groups as has recently been shown with S. hippopotami. The power of molecular approaches for the analysis of schistosomes at the population and individual level is now apparent, especially for S. mansoni. Important questions remain concerning the maintenance of parasite diversity and how schistosomes respond to selection pressures imposed either during natural progression through the life-cycle or through drug treatment or vaccination. Gene discovery and gene mapping projects are leading to a better understanding of the schistosome genome and can be expected to contribute significantly to future comparative evolutionary studies.

Non-random usage of 'degenerate' codons is related to protein three-dimensional structure

We report an analysis of a novel sequence-structure database of mammalian proteins incorporating nucleotide sequences of the exon regions of their genes together with protein sequence and structural information. We find that synonymous codon families (i.e. coding the same residue) have non-random codon distribution frequencies between protein secondary structure types. Their structural preferences are related to the third, 'silent' nucleotide position in a codon. We also find that some synonymous codons show very different or even opposite structural preferences at the N- or C-termini of structure fragments, relative to those observed for their amino acid residues.