Comment on "Large Bottleneck Size in Cauliflower Mosaic Virus Populations during Host Plant Colonization" by Monsion et al. (2008).

Our research groups have previously published a study estimating the bottleneck size of Cauliflower mosaic virus (CaMV) during host plant colonization [1]. Two methods were used in that study, both based on the temporal evolution of neutral marker frequencies; one uses estimates of F ST, and the other tracks changes in marker frequency variance over time. Here we report that (i) the variance-based method was actually published by Felsenstein [2], a reference we were not aware of when publishing the original study; (ii) equation (4) in [1] is actually an approximation; and (iii) these methods rely on the assumption that the bottleneck size is constant, which can be relaxed by assuming it is a Poisson-distributed random variable.

Equation (4) in [1] reads N = p(1−p)/[Var(p’)−Var(p)], where N is the estimated bottleneck size, and p and p’ represent the frequency of a neutral marker before and after the bottleneck, respectively. The approximation in this equation concerns the numerator. Its exact formulation should be E(p(1 − p)) and not E(p)(1 − E(p)), which is the expression implied in [1] and used in the numerical application on CaMV (E denotes the expected value). Based on the definition of the variance, it can be shown that E(p(1 − p)) = E(p)(1 − E(p)) − Var(p), which is equivalent to the numerator of equation (14) in [2]. Thus, the approximation will give relatively accurate results whenever Var(p) is negligible relative to E(p)(1 − E(p)). This was indeed the case in the study by Monsion et al. [1], as can be seen in Table 1, which compares the N values published in [1] to those obtained when applying the exact expression. The results of that paper are thus not qualitatively affected and, consequently, its conclusions remain unchanged. This method was also used in another study [3], the results of which are also marginally affected quantitatively, but not at all qualitatively (Table 2). To illustrate the application of these methods in a situation in which the estimated bottleneck size is lower, we used data from Tromas et al. (Table 3) [4].

Table 1

Comparison of estimates of N using the approximate formula published in [1] (values in Table 2 of this publication), those using the exact expression, and those assuming that bottleneck size is a zero-truncated-Poisson-distributed random variable.

Marker	N v approximate	N v exact	n vP
VIT1	423.00	421.16	422.16
VIT3	298.45	296.31	297.31
VIT4	484.61	481.42	482.42

Table 2

Comparison of estimates of N published in [3] using the approximate formula (values in table S1 of [3]), those using the exact expression, and those assuming that bottleneck size is a zero-truncated-Poisson-distributed random variable.

Note that due to copying errors, the values reported in table S1 of [3] for leaves 5 and 21 were slightly wrong—we report here the correct values.

Leaf level	N v approximate	N v exact	n vP
# 5	8.78	8.63	9.79
# 16	127.70	116.98	117.99
# 21	15.45	13.90	14.99

Table 3

Estimates of N using data published in [4] for leaf level 5.

The authors of that paper had estimated a bottleneck size of 6 (see their figure 5B) using the F ST method. To obtain the reported estimates, Leaf 3 and Leaf 5 were considered as the levels before and after the bottleneck, respectively. From the histogram 4A in [4], we calculated a mean frequency of 0.622 for the marker in Leaf 3, and we used the variances in marker frequency reported for Leaf 3 and Leaf 5.

Leaf level	N v approximate	N v exact	n vP
# 5	6.91	6.83	8.02

Both Felsenstein’s formula [2] and its approximation assume that the bottleneck size N is constant (i.e., identical in all the sampled plants). However, we can make the more realistic assumption that N follows a zero-truncated Poisson (ZTP) distribution (i.e., N can vary among experimental replicates according to a Poisson distribution, but p’ cannot be measured when N = 0 and the plant is discarded, hence the zero-truncation):

Among the N genomes that go through the bottleneck, the number X bearing the neutral marker is drawn according to its pre-bottleneck frequency p from the binomial distribution: X~B(N,p). Thus, the marker frequency after the bottleneck is p’ = X/N. The variance of this marker frequency can be expressed as: Var(p’) = Var[E(p’|N,p)]+E[Var(p’|N,p)].

Because E(p’|N,p) = E(X|N,p)/N and Var(p’|N,p) = Var(X|N,p)/N , and because N and p are independent, we then get: Var(p’) = Var(p)+E[p(1−p)]×E(1/N).

Using the properties and ,

where is the exponential integral function and γ the Euler-Mascheroni constant, one can estimate the bottleneck size n by numerically solving the following equation:

We recommend using this estimation procedure when N values are below 20, i.e., when the relative estimation error is higher than 5%. The following lines of R code may be used to estimate n :

library(gsl)

EulerMaschCst<−(−digamma(1))

p<−c(…,.,…) # vector of initial frequencies

p_prime<−c(…,.,…) # vector of final frequencies

funestim<−function(n_vP) (mean(p)*(1−mean(p))−var(p))/(var(p_prime)−var(p)) −

(exp(n_vP)−1)/(expint_Ei(n_vP)−log(n_vP)−EulerMaschCst)

uniroot(funestim,interval = c(1e−6,600))