1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
|
\name{cai}
\alias{cai}
\title{Codon Adaptation Index}
\description{
The Codon Adaptation Index (Sharp and Li 1987) is the most popular
index of gene expressivity with about 1000 citations 20 years after its
publication. Its values range from 0 (low) to 1 (high). The implementation
here is intended to work exactly as in the program \code{codonW} written by
by John Peden during his PhD thesis under the supervision of P.M. Sharp.
}
\usage{
cai(seq, w, numcode = 1, zero.threshold = 0.0001, zero.to = 0.01)
}
\arguments{
\item{seq}{a coding sequence as a vector of single characters}
\item{w}{a vector for the relative adaptiveness of each codon}
\item{numcode}{the genetic code number as in \code{\link{translate}}}
\item{zero.threshold}{a value in \code{w} below this threshold is
considered as zero}
\item{zero.to}{a value considered as zero in \code{w} is forced to
this value. The default is from Bulmer (1988).}
}
\details{
Adapted from the documentation of the CAI function in the
program \code{codonW} writen by John Peden:
CAI is a measurement
of the relative adaptiveness of the codon usage of a gene towards the
codon usage of highly expressed genes. The relative adaptiveness (w) of
each codon is the ratio of the usage of each codon, to that of the most
abundant codon for the same amino acid. The CAI
index is defined as the geometric mean of these relative adaptiveness
values. Non-synonymous codons and termination codons (genetic code
dependent) are excluded. To aid computation, the CAI is calculated as
using a natural log summation, To prevent a codon having a relative
adaptiveness value of zero, which could result in a CAI of zero;
these codons have fitness of zero (<.0001) are adjusted to 0.01.
}
\value{
A single numerical value for the CAI.
}
\references{
Sharp, P.M., Li, W.-H. (1987) The codon adaptation index - a
measure of directional synonymous codon usage bias, and its
potential applications.
\emph{Nucleic Acids Research}, \bold{15}:1281-1295.
Bulmer, M. (1988).
Are codon usage patterns in unicellular organisms determined by
selection-mutation balance.
\emph{Journal of Evolutionary Biology}, \bold{1}:15-26.
Peden, J.F. (1999)
Analysis of codon usage.
PhD Thesis, University of Nottingham, UK.
The program \code{codonW} used here for comparison is available at
\url{http://codonw.sourceforge.net/} under a GPL licence.
\code{citation("seqinr")}.
}
\seealso{
\code{\link{caitab}} for some \code{w} values from \code{codonW}.
\code{\link{uco}} for codon usage tabulation.
}
\author{J.R. Lobry}
\examples{
#
# How to reproduce the results obtained with the C program codonW
# version 1.4.4 writen by John Peden. We use here the "input.dat"
# test file from codonW (Saccharomyces cerevisiae).
#
inputdatfile <- system.file("sequences/input.dat", package = "seqinr")
input <- read.fasta(file = inputdatfile) # read the FASTA file
#
# Import results obtained with codonW
#
scucofile <- system.file("sequences/scuco.txt", package = "seqinr")
scuco.res <- read.table(scucofile, header = TRUE) # read codonW result file
#
# Use w for Saccharomyces cerevisiae
#
data(caitab)
w <- caitab$sc
#
# Compute CAI and compare results:
#
cai.res <- sapply(input, cai, w = w)
plot(cai.res, scuco.res$CAI,
main = "Comparison of seqinR and codonW results",
xlab = "CAI from seqinR",
ylab = "CAI from codonW",
las = 1)
abline(c(0,1))
}
|
