Codon usage bias and A+T content variation in human papillomavirus genomes

We analyzed the codon usage bias of eight open reading frames (ORFs) across up to 79 human papillomavirus (HPV) genotypes from three distinct phylogenetic groups. All eight ORFs across HPV genotypes show a strong codon usage bias, amongst degenerately encoded amino acids, toward 18 codons mainly with T at the 3rd position. For all 18 degenerately encoded amino acids, codon preferences amongst human and animal PV ORFs are significantly different from those averaged across mammalian genes. Across the HPV types, the L2 ORFs show the highest codon usage bias (73.2 +/- 1.6% and the E4 ORFs the lowest (51.1 +/- 0.5%), reflecting as similar bias in codon 3rd position A + T content (L2: 76.1 +/- 4.2%; E4: 58.6 +/- 4.5%). The E4 ORF, uniquely amongst the HPV ORFs, is G + C rich, while the other ORFs are A + T rich. Codon usage bias correlates positively with A + T content at the codon 3rd position in the E2, E6, L1 and L2 ORFs, but negatively in the E4 ORFs. A general conservation of preferred codon usage across human and non-human PV genotypes whether they originate from a same supergroup or not, together with observed difference between the preferred codon usage for HPV ORFs and for genes of the cells they infect, suggests that specific codon usage bias and A + T content variation may somehow increase the replicational fitness of HPVs in mammalian epithelial cells, and have practical implications for gene therapy of HPV infection. (C) 2003 Elsevier B.V. All rights reserved.

Exploring synonymous codon usage preferences of disulfide-bonded and non-disulfide bonded cysteines in the E-coli genome

High-quality data about protein structures and their gene sequences are essential to the understanding of the relationship between protein folding and protein coding sequences. Firstly we constructed the EcoPDB database, which is a high-quality database of Escherichia coli genes and their corresponding PDB structures. Based on EcoPDB, we presented a novel approach based on information theory to investigate the correlation between cysteine synonymous codon usages and local amino acids flanking cysteines, the correlation between cysteine synonymous codon usages and synonymous codon usages of local amino acids flanking cysteines, as well as the correlation between cysteine synonymous codon usages and the disulfide bonding states of cysteines in the E. coli genome. The results indicate that the nearest neighboring residues and their synonymous codons of the C-terminus have the greatest influence on the usages of the synonymous codons of cysteines and the usage of the synonymous codons has a specific correlation with the disulfide bond formation of cysteines in proteins. The correlations may result from the regulation mechanism of protein structures at gene sequence level and reflect the biological function restriction that cysteines pair to form disulfide bonds. The results may also be helpful in identifying residues that are important for synonymous codon selection of cysteines to introduce disulfide bridges in protein engineering and molecular biology. The approach presented in this paper can also be utilized as a complementary computational method and be applicable to analyse the synonymous codon usages in other model organisms. (c) 2005 Elsevier Ltd. All rights reserved.

Codon modified human papillomavirus type 16 E7 DNA vaccine enhances cytotoxic T-lymphocyte induction and anti-tumour activity

Polynucleotide immunisation with the E7 gene of human papillomavirus (HPV) type 16 induces only moderate levels of immune response, which may in part be due to limitation in E7 gene expression influenced by biased HPV codon usage. Here we compare for expression and immunogenicity polynucleotide expression plasmids encoding wild-type (pWE7) or synthetic codon optimised (pHE7) HPV16 E7 DNA. Cos-1 cells transfected with pHE7 expressed higher levels of E7 protein than similar cells transfected with pW7. C57BL/6 mice and F1 (C57X FVB) E7 transgenic mice immunised intradermally with E7 plasmids produced high levels of anti-E7 antibody. pHE7 induced a significantly stronger E7-specific cytotoxic T-lymphocyte response than pWE7 and 100% tumour protection in C57BL/6 mice, but neither vaccine induced CTL in partially E7 tolerant K14E7 transgenic mice. The data indicate that immunogenicity of an E7 polynucleotide vaccine can be enhanced by codon modification. However, this may be insufficient for priming E7 responses in animals with split tolerance to E7 as a consequence of expression of E7 in somatic cells. (C) 2002 Elsevier Science (USA).

Mitochondrial gene content, arrangement and composition compared in African and Asian schistosomes

Complete sequences were obtained for the coding portions of the mitochondrial (mt) genomes of Schistosoma mansoni (NMRI strain, Puerto Rico; 14415 bp), S. japonicum (Anhui strain, China; 14085 bp) and S. mekongi (Khong Island, Laos; 14072 bp). Each comprises 36 genes: 12 protein-encoding genes (cox1-3, nad1-6, nad4L, atp6 and cob); two ribosomal RNAs, rrnL (large subunit rRNA or 16S) and rrnS (small subunit rRNA or 12S); as well as 22 transfer RNA (tRNA) genes. The atp8 gene is absent. A large segment (9.6 kb) of the coding region (comprising 14 tRNAs, eight complete and two incomplete protein-encoding genes) for S. malayensis (Baling, Malaysian Peninsula) was also obtained. Each genome also possesses a long non-coding region that is divided into two parts (a small and a large non-coding region, the latter not fully sequenced in any species) by one or more tRNAs. The protein-encoding genes are similar in size, composition and codon usage in all species except for cox1 in S. mansoni (609 aa) and cox2 in S. mekongi (219 an), both of which are longer than homologues in other species. An unexpected finding in all the Schistosoma species was the presence of a leucine zipper motif in the nad4L gene. The gene order in S. mansoni is strikingly different from that seen in the S. japonicum group and other flatworms. There is a high level of identity (87-94% at both the nucleotide and amino acid levels) for all protein-encoding genes of S. mekongi and S. malayensis. The identity between genes of these two species and those of S. japonicum is less (56-83% for amino acids and 73-79 for nucleotides). The identity between the genes of S. mansoni and the Asian schistosomes is far less (33-66% for amino acids and 54-68% for nucleotides), an observation consistent with the known phylogenetic distance between S. mansoni and the other species. (C) 2001 Elsevier Science B.V. All rights reserved.

Whole genome analysis reveals a high incidence of non-optimal codons in secretory signal sequences of Escherichia coli

Translational pausing may occur due to a number of mechanisms, including the presence of non-optimal codons, and it is thought to play a role in the folding of specific polypeptide domains during translation and in the facilitation of signal peptide recognition during see-dependent protein targeting. In this whole genome analysis of Escherichia coli we have found that non-optimal codons in the signal peptide-encoding sequences of secretory genes are overrepresented relative to the mature portions of these genes; this is in addition to their overrepresentation in the 5'-regions of genes encoding non-secretory proteins. We also find increased non-optimal codon usage at the 3' ends of most E. coli genes, in both non-secretory and secretory sequences. Whereas presumptive translational pausing at the 5' and 3' ends of E. coli messenger RNAs may clearly have a general role in translation, we suggest that it also has a specific role in sec-dependent protein export, possibly in facilitating signal peptide recognition. This finding may have important implications for our understanding of how the majority of non-cytoplasmic proteins are targeted, a process that is essential to all biological cells. (C) 2004 Elsevier Inc. All rights reserved.

Investigation of recombinant Schistosoma japonicum paramyosin fragments for immunogenicity and vaccine efficacy in mice

Schistosoma japonicum paramyosin, a 97 kDa myofibrillar protein, is a recognized vaccine candidate against schistosomiasis. To improve its expression and to identify protective epitopic regions on paramyosin, the published Chinese Schistosoma japonicum paramyosin cDNA sequence was redesigned using Pichia codon usage and divided into four overlapping fragments (fragments 1, 2, 3, 4) of 747, 651, 669 and 678 bp, respectively. These gene fragments were synthesized and expressed in Pichia pastoris (fragments 2 and 3) or E. coli (fragments 1 and 4). The recombinant proteins were produced at high level and purified using a two-step process involving Ni-NTA affinity chromatography and gel filtration. BALB/c mice were immunized subcutaneously three times at 2-week-intervals with the purified proteins formulated in adjuvant Quil A. The protein fragments were highly immunogenic, inducing high, though variable, ELISA antibody titres, and each was shown to resemble native paramyosin in terms of its recognition by the anti-fragment antibodies in Western blotting. The immunized mice were subjected to cercarial challenge 2 weeks after the final injection and promising protective efficacy in terms of significant reductions in worm burdens, worm-pair numbers and liver eggs in the vaccinated mice resulted. There was no apparent correlation between the antibody titres generated and protective efficacy, as all fragments produced effective but similar levels of protection.

Incorporating a TEV cleavage site reduces the solubility of nine recombinant mouse proteins

Failure to express soluble proteins in bacteria is mainly attributed to the properties of the target protein itself, as well as the choice of the vector, the purification tag and the linker between the tag and protein, and codon usage. The expression of proteins with fusion tags to facilitate subsequent purification steps is a widely used procedure in the production of recombinant proteins. However, the additional residues can affect the properties of the protein; therefore, it is often desirable to remove the tag after purification. This is usually done by engineering a cleavage site between the tag and the encoded protein that is recognised by a site-specific protease, such as the one from tobacco etch virus (TEV). In this study, we investigated the effect of four different tags on the bacterial expression and solubility of nine mouse proteins. Two of the four engineered constructs contained hexahistidine tags with either a long or short linker. The other two constructs contained a TEV cleavage site engineered into the linker region. Our data show that inclusion of the TEV recognition site directly downstream of the recombination site of the Invitrogen Gateway vector resulted, in a loss of solubility of the nine mouse proteins. Our work suggests that one needs to be very careful when making modifications to expression vectors and combining different affinity and fusion tags and cleavage sites: (c) 2006 Elsevier Inc. All rights reserved.

Cardiac expression of the cystic fibrosis transmembrane conductance regulator involves novel Exon 1 usage to produce a unique amino-terminal protein

Cystic fibrosis is caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene, which encodes a chloride channel present in many cells. In cardiomyocytes, we report that multiple exon 1 usage and alternative splicing produces four CFTR transcripts, with different 5'-untranslated regions, CFTRTRAD-139, CFTR-1C/-1A, CFTR-1C, and CFTR-1B. CFTR transcripts containing the novel upstream exons (exons -1C, -1B, and -1A) represent more than 90% of cardiac expressed CFTR mRNA. Regulation of cardiac CFTR expression, in response to developmental and pathological stimuli, is exclusively due to the modulation of CFTR-1C and CFTR-1C/-1A expression. Upstream open reading frames have been identified in the 5'-untranslated regions of all CFTR transcripts that, in conjunction with adjacent stem-loop structures, modulate the efficiency of translation initiation at the AUG codon of the main CFTR coding region in CFTRTRAD-139 and CFTR-1C/-1A transcripts. Exon(-1A), only present in CFTR-1C/-1A transcripts, encodes an AUG codon that is in-frame with the main CFTR open reading frame, the efficient translation of which produces a novel CFTR protein isoform with a curtailed amino terminus. As the expression of this CFTR transcript parallels the spatial and temporal distribution of the cAMP-activated whole-cell current density in normal and diseased hearts, we suggest that CFTR-1C/-1A provides the molecular basis for the cardiac cAMP-activated chloride channel. Our findings provide further insight into the complex nature of in vivo CFTR expression, to which multiple mRNA transcripts, protein isoforms, and post-transcriptional regulatory mechanisms are now added.

Selection bias in gene extraction on the basis of microarray gene-expression data

In the context of cancer diagnosis and treatment, we consider the problem of constructing an accurate prediction rule on the basis of a relatively small number of tumor tissue samples of known type containing the expression data on very many (possibly thousands) genes. Recently, results have been presented in the literature suggesting that it is possible to construct a prediction rule from only a few genes such that it has a negligible prediction error rate. However, in these results the test error or the leave-one-out cross-validated error is calculated without allowance for the selection bias. There is no allowance because the rule is either tested on tissue samples that were used in the first instance to select the genes being used in the rule or because the cross-validation of the rule is not external to the selection process; that is, gene selection is not performed in training the rule at each stage of the cross-validation process. We describe how in practice the selection bias can be assessed and corrected for by either performing a cross-validation or applying the bootstrap external to the selection process. We recommend using 10-fold rather than leave-one-out cross-validation, and concerning the bootstrap, we suggest using the so-called. 632+ bootstrap error estimate designed to handle overfitted prediction rules. Using two published data sets, we demonstrate that when correction is made for the selection bias, the cross-validated error is no longer zero for a subset of only a few genes.

PCR bias toward the wild-type k-ras and p53 sequences: Implications for PCR detection of mutations and cancer diagnosis

PCR-based cancer diagnosis requires detection of rare mutations in k-ras, p53 or other genes. The assumption has been that mutant and wild-type sequences amplify with near equal efficiency, so that they are eventually present in proportions representative of the starting material. Work factor IX suggests that this assumption is invalid for one case of near-sequence identity To test the generality of this phenomenon and its relevance to cancer diagnosis, primers distant from point mutations in p53 and k-ras were used to amplify, wild-type and mutant sequences from these genes. A substantial bias against PCR amplification of mutants was observed for two regions of the p53 gene and one region of k-ras. For kras and p53, bias was observed when the wild-type and mutant sequences were amplified separately or when mixed in equal proportions before PCR. Bias was present with proofreading and non-proofreading polymerases. Mutant and wild-type segments of the factor V cystic fibrosis transmembrane conductance regulator and prothrombin genes were amplified and did not exhibit PCR bias. Therefore, the assumption of equal PCR efficiency for point mutant and wild-type sequences is invalid in several systems. Quantitative or diagnostic PCR will require validation for each locus, and enrichment strategies may be needed to optimize detection of mutants.

The novel mitochondrial gene arrangement of the cattle tick, Boophilus microplus: Fivefold tandem repetition of a coding region

We sequenced across all of the gene boundaries in the mitochondrial genome of the cattle tick, Boophilus microplus, to determine the arrangement of its genes. The mtDNA of B. microplus has a coding region, composed of tRNA(Glu) and 60 bp of the 3' end of ND1, that is repeated five times. Boophilus microplus is the first coelomate animal known to have more than two copies of a coding sequence. The mitochondrial genome of B, microplus has other unusual features, including (1) reduced T arms in tRNAs, (2) an AT bias in codon use, (3) two control regions that have evolved in concert, (4) three gene rearrangements, and (5) a stem-loop between tRNA(Gln) and tRNA(Phe). The short T arms and small control regions (CRs) of B. microplus and other ticks suggest strong selection for small genomes. Imprecise termination of replication beyond its origin, which can account for the evolution of tandem repeats of coding regions in other mitochondrial genomes, cannot explain the evolution of the fivefold repeated sequence in the mitochondrial genome of B. microplus. Instead, slipped-strand mispairing or recombination are the most plausible explanations for the evolution of these tandem repeats.

Mothers' mental illness and child behavior problems: Cause-effect association or observation bias

Objective: A number, of studies have consistently found that a mother's mental health (particularly her level of depression) is a strong predictor of mental health problems experienced by her child(ren). However, the validity of this finding is in doubt because the majority of these studies have relied on maternal reports as indicators of children's behavior. Method: This prospective, longitudinal study examines data an the mental health of the mother from prior to the birth of her child to when the child reaches 14 years of age. Child behavior is measured at 14 years of age using reports from mother and child. Mother and child responses are compared to provide an indication of the possible magnitude of maternal observation bias in the reporting of child behavior problems. Results: Anxious and/or depressed mothers tend to report more cases of child behavior problems than do their mentally healthy counterparts or children themselves. Differences between mothers and youths In reporting behavior problems appear to be related to the mothers' mental health. Conclusions: Current maternal mental health impairment appears to have a substantial effect on the reporting of child behavior problems by the mother, thereby raising questions about the validity of reports of child behavior by persons who are currently emotionally distressed.

Investigation of threat-related attentional bias in anxious children using the startle eyeblink modification paradigm

We used the startle eyeblink modification paradigm to investigate whether clinically anxious children, like high trait-anxious adults, display a bias in favour of threat words compared to neutral words. The present study included 16 clinically anxious children whose diagnostic status was determined using the parent version of a semistructured diagnostic interview as part of a larger childhood anxiety study. The children were presented with threat and neutral words fur 6 s each. A startle-eliciting auditory stimulus - a 100 dBA burst of white noise of 50 ms duration - was presented during the words at lead intervals of 60, 120, 240, or 3500 ms and during intertrial intervals. The overall pattern of startle eyeblink modification indicated inhibition at the 120 and 240 ms lead intervals and facilitation at the 3500 ms lead interval. startle-latency shortening during threat words at the :60 ms lead interval was larger than at other intervals, whereas there was no difference during neutral words. This result reflects an anxiety-related bias in favour of threat words occurring at a very early - and possibly preattentive stage - of information processing.

Size-based predation by kookaburras (Dacelo novaeguineae) on lizards (Eulamprus tympanum : Scincidae): what determines prey vulnerability?

Lizards and birds are both popular model organisms in behavioural ecology, but the interactions between them have attracted little study. Given the putative importance of birds as predators of diurnal Lizards, it is of considerable interest to know which traits (of lizards as well as birds) influence the outcome of a predatory attempt. We studied predation by giant terrestrial kingfishers (kookaburras, Dacelo novaeguineae: Alcedinidae) on heliothermic diurnal lizards (highland water skinks, Eulamprus tympanum: Scincidae), with particular reference to the role of prey (lizard) size. Our approach was twofold: to gather direct evidence (sizes of lizards consumed in the field, compared to those available) and indirect evidence rite-related shifts in lizard behaviour). We quantified the size structure of a natural population of skinks (determined by an extensive mark-recapture program), and compared it to the sizes of wild lizards taken by kookaburras (determined by analysis of prey remains left at the birds' nests,. Kookaburras showed size-based predation: they preyed mainly on small and medium-sized rather than large lizards in the field. However, the mechanism producing this bias remains elusive. It is not due to any distinctive behavioural attributes (locomotor ability, activity level, habitat usage) of the lizards of the size class disproportionately taken by the kookaburras. The greater vulnerability of subadult lizards may reflect subtle ontogenetic shifts in ecological and behavioural traits, but our data suggest that great caution is needed in inferring patterns of vulnerability to predation from indirect measures based on either the prey or the predator alone. Instead, we need direct observations on the interaction between the two.