Amélie Charles
email: acharles@audencia.com.

Olivier Darné
email: olivier.darne@univ-nantes.fr.

Lemna

The efficiency of the crude oil markets: Evidence from variance ratio tests

Keywords: Crude oil markets, market efficiency, random walk hypothesis, variance ratio test. JEL Classification: Q40, C14

This study examines the random walk hypothesis for the crude oil markets, using daily data over the period 1982-2008. The weak-form efficient market hypothesis for two oil crude markets (UK Brent and US West Texas Intermediate) is tested with non-parametric variance ratio tests developed by Wright (2000) and Belaire-Franch and Contreras (2004) as well as the wild-bootstrap variance ratio tests suggested by Kim (2006). We find that the Brent crude oil market is weak-form efficiency while the WTI crude oil market seems to be inefficiency on the 1994-2008 sub-period, suggesting that the deregulation have not improved the efficiency on the WTI crude oil market in the sense of making returns less predictable.

Introduction

Since the end of the 1990s oil prices have been steadily increasing, reflecting rising demand for crude oil, particularly from developing nations. Indeed, in the last decade the oil prices increased more than 700% (in average) and more than 80% (in average) between mid 2007 to mid 2008 due to rising demand, low spare capacity, weak dollar and geopolitical concerns (especially, tensions in Turkey, Nigeria and Iran). Oil prices have been very volatile, changing their trajectories and behavior with respect to the economic situation. Oil prices exhibit large upward or downward swings primarily caused by fluctuations in demand, extraction costs, and reserves (Pindyck, 1999).

Supply and demand remain the main factors determining oil prices. More precisely, oil demands depend on oil consumption by developed and developing countries, and oil supplies depend on geopolitical events, such as troubles between Venezuela and the US or Turkey and Kurdish Iraq or Iran and Israel, among others, as well as oil tank levels and the Organization of the Petroleum Exporting Countries (OPEC) decisions on adjusting production levels. However, oil investor behavior can be also a factor in the recent price behavior, especially increasingly speculative behavior of a more diverse set of investors, including hedge funds, pension funds, and investment banks. All these factors question on the issue of whether or not the crude oil markets are predictable and therefore efficient.

In this paper we analyze the efficiency of the crude oil markets. The literature on market efficiency and stock market predictability is vast, as researchers have been discussing this theme in depth for the past decades (see Fama, 1970; 1991; Fama and French, 1988; Lo and MacKinlay, 1988; among others). A capital market is considered as efficient if stock prices at any time fully reflect all available and relevant information. Therefore, given only past price and return data, the current price is the best predictor of the future price, and the price change or return is expected to be zero. Stock prices exhibit no serial dependencies, meaning that there are no patterns to asset prices. This implies that future price movements are determined entirely by information not contained in the price series. This is the essence of the weakform efficient market hypothesis [EMH], which implies a random walk. It is this random-walk implication of the weak-form EMH which is most commonly tested in the empirical literature [START_REF] Alvarez-Ramirez | Short-term predictability of crude oil markets: A detrended fluctuation analysis approach[END_REF] .

Recently, Tabak and Cajueiro (2007) We also investigate the EMH over two sub-periods in order to analyze the effects of the important structural change due to policy changes that attempted to increase the efficiency of the North American energy industry in 1993 [START_REF] Alvarez-Ramirez | Multifractal Hurst analysis of crude oil prices[END_REF] . Second, the weak form of the EMH is evaluated from an alternative method relative to the previous studies, namely the variance-ratio [VR] tests which are widely used in financial empirical studies (e.g., Lagoarde-Segot and Lucey, 2007; Kim and Shamsuddin, 2008) [START_REF] Belaire-Franch | Ranks and signs-based multiple variance ratio tests[END_REF] . More precisely, we adopt individual non-parametric VR tests suggested by Wright (2000) as well as its multiple version proposed by Belaire-Franch and Contreras (2004).

These VR tests are robust to heteroscedasticity and non-normality which are features displayed by the oil crude prices (Pindyck, 2004; Narayan and Narayan, 2007; Kang, Kang and Yoon, 2009), and powerful against fractionally integrated alternatives which are present in oil crude prices (Alvarez-Ramirez et al., 2002, 2008; Tabak and Cajueiro, 2007). We also apply the wild bootstrap VR tests suggested by Kim (2006) which are robust to heteroscedasticity and do not rely on asymptotic approximations.

The rest of this paper is organized as follows. Section 2 discusses the variance ratio tests. Section 3 summarizes the characteristics of the data on the Brent and WTI crude oil markets. Section 4 reports the empirical results. Section 5 concludes.

Variance ratio tests

Since the seminal work of Lo and MacKinlay (1988, 1989) and Poterba and Summers (1988), the standard variance ratio [VR] test or its improved modifications have been [START_REF] Alvarez-Ramirez | Multifractal Hurst analysis of crude oil prices[END_REF] See Serletis and Andreadis (2004) and Serletis and Rangel-Ruiz (2004) for a discussion on these policy changes. [START_REF] Belaire-Franch | Ranks and signs-based multiple variance ratio tests[END_REF] widely used for testing market efficiency [START_REF] Charles | Variance ratio tests of random walk: An overview[END_REF] .

The VR methodology consists of testing the random walk hypothesis [RWH] against stationary alternatives, by exploiting the fact that the variance of random walk increments is linear in all sampling intervals, i.e., the sample variance of k-period return (or k-period differences), y ty t-k , of the time series y t , is k times the sample variance of one-period return (or the first difference), y ty t-1 . The VR at lag k is then defined as the ratio between (1/k)th of the k-period return (or the kth difference) to the variance of the one-period return (or the first difference). Hence, for a random walk process, the variance computed at each individual lag interval k (k = 2, 3, . . . ) should be equal to unity.

In testing the null hypothesis of random walk, the VR test evaluates the hypothesis that a given time series or its first difference (or return), x t = y ty t-1 , is a collection of independent and identically distributed observations (i.i.d.) or that it follows a martingale difference sequence. Following Wright (2000), the VR statistic be written as

V R(x; k) = (T k) -1 T ∑ t=k (x t + • • • + x t-k+1 -k μ) 2 ÷ T -1 T ∑ t=1 (x t -μ) 2 (1)
where μ = T -1 ∑ T t=1 x t . If the stock return follows a random walk, the expected value of V R(x; k) should be equal to unity for all horizons k. If this ratio is less than one at long horizons, then we have indications of negative serial correlation (mean-reversion) and ratios greater than one at long horizons implies positive serial correlation (meanaversion). Lo and MacKinlay (1988) proposed the asymptotic distribution of V R(x; k) by assuming that k is fixed when T → ∞. They show that under the assumption of conditional heteroscedasticity, then under the null hypothesis that

V (k) = 1, the test statistic M(x; k) is given by 5 M(x; k) = V R(x; k) -1 φ * (k) 1/2 (2)
follows the standard normal distribution asymptotically, where Wright (2000) showed that rank-based tests display low size distortion, under conditional heteroscedasticity [START_REF] Davidson | The wild bootstrap, tamed at last[END_REF] .

φ * (k) = k-1 ∑ j=1 2(k -j) k 2 δ( j) δ( j) = T ∑ t= j+1 (x t -μ) 2 (x t-j -μ) 2 ÷    T ∑ t=1 (x t -μ) 2 2    2.
Given T observations of first differences of a variable, {x 1 , . . . , x T }, and let r(x) be the rank of x t among (x 1 , . . . , x T ). Under the null hypothesis that x t is generated from an i.i.d. sequence, r(x) is a random permutation of the numbers of 1, . . . , T with equal probability. Wright (2000) suggested the R 1 and R 2 statistics, defined as

R 1 (k) = (T k) -1 ∑ T t=k (r 1,t + • • • + r 1,t-k+1 ) 2 T -1 ∑ T t=k r 2 1,t -1 × φ(k) -1/2 (3) R 2 (k) = (T k) -1 ∑ T t=k (r 2,t + • • • + r 2,t-k+1 ) 2 T -1 ∑ T t=k r 2 2,t -1 × φ(k) -1/2 (4)
where the standardized ranks r 1,t and r 2,t are given by

r 1,t = r(x t ) -T +1 where φ(k) = 2(2k -1)(k -1)(3kT ) -1 ,
joint null hypothesis V (k i ) = 1 for i = 1, . . . , m against the alternative that V (k i ) = 1
for some holding period k i . The (heteroscedasticity-robust) test statistic can be written as The statistics are defined as

MV (x; k i ) = √ T max 1≤i≤m |M(x; k i )| (5)
CD (R 1 ) = max 1≤i≤m |R 1 (k i )| CD (R 2 ) = max 1≤i≤m |R 2 (k i )|
The rank-based procedures are exact under the i.i.d. assumption.

Kim (2006) test

Kim (2006) uses the wild bootstrap which is a resampling method that approximates the sampling distribution of the VR test statistic, and is applicable to data with unknown forms of conditional and unconditional heteroscedasticity (Mammen, 1993; Davidson and Flachaire, 2001; MacKinnon, 2002).

The wild bootstrap is applied to Chow-Denning, MV (x; k i ), VR test. The wild bootstrap test based on MV (x; k i ) can be conducted in three stages as below (i) Form a bootstrap sample of T observations x * t = η t x t (t = 1, . . . , T ) where η t is a random sequence with E(η) = 0 and E(η 2 ) = 1.

(ii) Calculate MV * = MV (x * ; k i ), the MV (x * ; k i ) statistic obtained from the bootstrap sample generated in stage (i).

(iii) Repeat (i) and (ii) sufficiently many, say m, times to form a bootstrap distribution of the test statistic {MV ( j)} m j=1 .

The bootstrap distribution {MV ( j)} m j=1 is used to approximate the sampling distribution of the MV (x; k i ) statistic. The p-value of the test can be obtained as the proportion of {MV ( j)} m j=1 greater than the MV (x; k i ) statistic calculated from the original data. The wild bootstrap version of Lo-MacKinlay M(k) test can be implemented in a similar manner as a two-tailed test, where we obtain M * = M(x * ; k) in stage (ii) and {M(x * ; k j ; j)} m j=1 in stage (iii). Conditionally on X t , X * t is a serially uncorrelated sequence with zero mean and variance X 2 t . As such, M * and MV * have the same asymptotic distributions as M 2 (k) and MV 2 (k i ), respectively. Since X * t is a serially uncorrelated sequence, wild bootstrapping approximates the sampling distributions under the null hypothesis, which is a desirable property for a bootstrap test.

To implement the wild bootstrap test, a specific form of η t should be chosen. Kim (2006) recommends using the standard normal distribution for η t since he reports that other choices provided qualitatively similar small sample results [START_REF] Fama | Efficient capital markets: A review of theory and empirical work[END_REF] .

Data and summary statistics

The data of the study consists of the daily closing spot prices for two oil crude markets: the US West Texas Intermediate (WTI) and the UK Brent. The data comes from Thomson Financial Datastream and is given in US dollar per barrel. The data spans from 1 June, 1982 to 4 July, 2008, namely 6809 observations. Figures 1 and2 Brent crude oils are light crude oils, with WTI lighter than Brent (the API9 gravity is around 39.6 degrees for WTI against around 38.3) and sweet crude oils, with Brent slightly less sweet than WTI (Brent contains 0.37% of sulfur against 0.24% for WTI).

Due to its very high quality, WTI crude oil is excellent for refining a larger portion of gasoline whereas Brent is ideal for making gasoline and middle distillates. The Brent crude oil is actually a combination of crude oil from 15 different oil fields in the Brent and Ninian systems located in the North Sea. Most WTI crude oil gets refined in the Midwest region of the US, more precisely at Cushing (Oklahoma), with some more refined within the Gulf Coast region. Brent crude oil is typically refined in Northwest Europe.

Table 1 presents summary statistics for the WTI and Brent crude oil returns calculated as the first differences in the logs of the spot prices. The two crude oil markets have approximately equal mean returns of about 0.02% per day, with the WTI returns actually marginally smaller than the Brent returns. All the returns are highly non-normal, i.e. showing evidence of significant negative skewness and excess kurtosis, as might be expected from daily returns. The Lagrange Multiplier test for the presence of the ARCH effect indicates clearly that all crude oil prices show strong conditional heteroscedasticity, which is a common feature of financial data.

Testing the efficient market hypothesis

We investigate the weak-form efficient market hypothesis by testing the random-walk hypothesis from VR tests. Table 2 reports Table 3 presents summary statistics for the two crude oil markets on the two subperiods. For the Brent market, the returns on the 1982-1993 sub-period are higher (-0.021% per day) than those on the second sub-period (-0.031% per day). In the contrast, the returns of the WTI market on the two sub-periods have approximately equal mean returns of about 0.03% per day. Note that the risk (in terms of the standard deviation) associated with the 1994-2008 sub-period for the two crude oil markets is higher (2.30% per day) than that associated with the 1982-1993 sub-period (around 1.0% per day).

On the whole, all the statistics increase over the two sub-periods, especially for the excess kurtosis and the conditional heteroscedasticity, suggesting an impact of the deregulation on the crude oil prices. Furthermore, the skewness is negative and high, implying an increase in the extreme losses for the investors. The kurtosis is very high for the second sub-period, implying that the variance of the crude oil prices is principally due to infrequent but extreme deviations.

Tables 4 and5 display the results of individual and multiple VR tests [START_REF] Fama | Permanent and temporary components of stock prices[END_REF] . For the Brent crude oil market, the VR tests show the efficiency on the two sub-periods. On the contrary, for the WTI crude oil market, the RWH is rejected on the 1994-2008

sub-period while it is not rejected for the 1982-1993 sub-period, suggesting that the deregulation have not improved the efficiency on the WTI crude oil market in the sense of making returns less predictable. Thus, the NAFTA deregulation does not seem to have a positive impact on the efficiency of the WTI crude oil market, contrarily to what was expected. Furthermore, the Brent crude oil market seems to be more weakform efficiency than the WTI crude oil market since the process of deregulation. This result can be explained by the fact that the WTI price can be affected by its storage tank capacity and its infrastructure logistics (Horsnell and Mabro, 1993). For example, high levels of stocks at Cushing (e.g., due to unsatisfactory supply-demand conditions) tend to depress the WTI price relative to the Brent price because the WTI crude oil market suffers from pipeline logistic, more precisely infrastructure bottlenecks. These results are in contradiction with those of Tabak and Cajueiro (2007), finding that the WTI crude oil prices seems to be more weak-form efficient than Brent prices.

Conclusion

This paper has explored weak efficiency in the Brent and WTI crude oil markets and the relationship between efficiency and changes in the regulation, using new variance The results suggest that the Brent crude oil prices is weak-form efficiency while the WTI crude oil market seems to be inefficiency on the 1994-2008 sub-period, suggesting that the deregulation has not improved the efficiency on the WTI crude oil market in the sense of making returns less predictable. The storage tank capacity and infrastructure logistics of the WTI crude oil market can explain the fact that the Brent crude oil market is more efficient than the WTI crude oil market.

  , Alvarez-Ramirez et al. (2008) and Maslyuk and Smyth (2009) investigated the efficiency of crude oil markets from timevarying long-range dependence, Hurst exponent dynamics (from detrended fluctuation analysis) and unit root tests, respectively. They found that these markets are weakform efficient. In this paper, we extend the examination of the weak-form of the EMH in the crude oil markets in two ways. First, this study is based on a more extensive sample. We study daily data for two crude oil markets, namely, the US West Texas Intermediate and the UK Brent, over the period June 1982 -July 2008.

  where M(x; k i ) is defined in (2). Recently, Belaire-Franch and Contreras (2004) proposed a multiple rank VR tests by substituting Wright's rank-based tests in the definition of the Chow-Denning (1993) multiple test procedure 7 .

  provide a graphical representation of these series. The WTI and Brent crude oils are traded on New York Mercantile Exchange (NYMEX) and InterContinental Exchange (ICE, London), respectively. The WTI and 8 Kim (2006) finds that his simulation results are largely insensitive to different choices of η t , including the two-point distribution of Mammen (1993), and the Rademacher distribution discussed in Davidson and Flachaire (2001).

  the individual and multiple non-parametric VR test results. The test statistics are displayed for the Wright (R 1 and R 2 ) and Belaire-Contreras (CD(R 1 ) and CD(R 2 )) tests while the p-values are reported for the Kim (M * and MV * ) tests. The holding periods (k i 's) considered are (2, 5, 10, 30). As advocated by Deo and Richardson (2003), we use relatively short holding periods when testing for the mean reversion using VR tests.

Figure 1 :Figure 2 :

 12 Figure 1: Daily Brent crude oil prices

  ratio tests. Theses tests, which are robust to heteroscedasticity and non-normality and powerful against fractionally integrated alternatives which are present in oil crude prices, are the non-parametric VR tests developed by Wright (2000) and Belaire-Franch and Contreras (2004) as well as the wild-bootstrap VR tests suggested by Kim (2006).

  1 Wright (2000) tests A well-known problem with the VR test is that the standard VR tests such as Lo and MacKinlay (1988) tests, which are based on asymptotic approximations, are biased

(severe size distortions and low power) and right-skewed in finite samples, resulting in misleading statistical inference.

Wright (2000)

proposed a non-parametric alternative to conventional asymptotic VR tests using ranks.

Wright's (2000)

tests have two advantages over Lo-MacKinlay test when sample size is relatively small: (i) as the rank (R 1 and R 2 ) tests have exact sampling distribution, there is no need to resort to asymptotic distribution approximation, and (ii) the tests may be more powerful than the conventional VR tests against a wide range of models displaying serial correlation, including fractionally integrated alternatives. The tests based on ranks are exact under the i.i.d. assumption. Moreover,

  and Φ -1 is the inverse of the standard normal cumulative distribution function. The R 1 and R 2 statistics follow the same exact sampling distribution. The critical values of these tests can be obtained by simulating their exact distributions.It is preferable to employ multiple VR tests rather than individual VR tests. Indeed, in practice, it is customary to examine the VR statistics for several holding periods (k) to test the RWH. The null is rejected if it is rejected for some k value. As stressed byChow and Denning (1993), this sequential procedure leads to size distortions. To overcome this problem, Chow and Denning (1993) provide a multiple VR test for the

	2.2 Belaire-Franch and Contreras (2004) tests

Table 1 :

 1 Statistical analysis of returnsNotes: The skewness and kurtosis statistics are standard-normally distributed under the null of normality distributed returns. ARCH(10) indicates the Lagrange multiplier test for conditional heteroscedasticity with 10 lags. LB(10) corresponds to the Ljung-Box test for serial autocorrelation with 10 lags. * and * * means significant at 1% and 5% level, respectively.

	Market Obs.	Mean	SD	Skewness Kurtosis LB(10) ARCH(10)
			(% daily) (% daily)
	Brent	6810 0.0209	2.20	-1.031 *	30.86 *	41.30 *	158.16 *
	WTI	6810 0.0203	2.30	-1.040 *	25.63 *	40.02 *	206.31 *

Table 2 :

 2 Results of individual and multiple VR tests -1982-2008.MV * 0.037 * * and * * Significant at the 5% and 10% level, respectively. We report the VR statistic for each test.

				k
			2	5	10	30
	Brent	R 1	0.903	-0.233	0.254	1.724
		R 2	0.836	-0.354	-0.356	1.111
		M *	0.494	0.578	0.255	0.776
		CD(R 1 )	1.270
		CD(R 2 )	0.852
		MV *	0.529
	WTI	R 1	-1.386	-1.749	-1.190	0.497
		R 2	-1.274	-2.064 *	-2.165 *	-0.499
		M *	0.717	0.071 * *	0.012 *	0.138
		CD(R 1 )	1.974 * *
		CD(R 2 )	2.236 * *

Table 3 :

 3 Statistical analysis of returnsNotes: The skewness and kurtosis statistics are standard-normally distributed under the null of normality distributed returns. ARCH(10) indicates the Lagrange multiplier test for conditional heteroscedasticity with 10 lags. LB(10) corresponds to the Ljung-Box test for serial autocorrelation with 10 lags.

	Market	Obs.	Mean	SD	Skewness Kurtosis LB(10) ARCH(10)
			(% daily) (% daily)
					Brent
	1982-1993 937	-0.021	0.90	-1.686 *	21.41 *	40.49 *	69.07 *
	1994-2008 3025	-0.031	2.30	-1.913 *	55.39 *	18.05 *	142.99 *
					WTI
	1982-1993 937	-0.033	1.00	-1.519 *	38.71 *	35.93 *	95.81 *
	1994-2008 3025	-0.030	2.30	-1.882 *	45.18 *	21.71 *	112.52 *

* and * * means significant at 1% and 5% level, respectively.

Table 4 :

 4 Results of individual and multiple VR tests -Brent.

				k
			2	5	10	30
	1982-1993	R 1	1.831 *	0.680	0.908	2.609 *
		R 2	1.920 *	0.617	0.394	2.065 *
		M *	0.274	0.912	0.402	0.896
		CD(R 1 )	1.876
		CD(R 2 )	1.962
		MV *	0.503
	1994-2008	R 1	-0.350	-1.043	-0.688	-0.226
		R 2	-0.564	-1.190	-0.988	-0.419
		M *	0.470	0.238	0.265	0.595
		CD(R 1 )	1.187
		CD(R 2 )	1.257
		MV *	0.517

Table 5 :

 5 Results of individual and multiple VR tests -WTI. 0.113 * and * * Significant at the 5% and 10% level, respectively. We report the VR statistic for each test.

				k
			2	5	10	30
	1982-1993	R 1	0.862	0.103	0.561	2.109 *
		R 2	0.601	-0.400	-0.445	1.376
		M *	0.940	0.171	0.054 *	0.509
		CD(R 1 )	0.959
		CD(R 2 )	0.686
		MV *	0.152
	1994-2008	R 1	-2.187 *	-2.346 *	-2.007 *	-1.253
		R 2	-1.811	-2.242 *	-2.303 *	-1.828 *
		M *	0.362	0.118	0.030 *	0.051 *
		CD(R 1 )	2.569 *
		CD(R 2 )	2.527

*

MV *

Note that if the random walk hypothesis is based on the theory of efficiency, the EMH does not imply that prices follow a random walk. Therefore, if prices do not follow a random walk, this does not imply inefficiency of the market. See[START_REF] Lo | A Non-Random Walk Down Wall Street[END_REF] for a discussion on random walk hypothesis and efficiency market hypothesis.

See Hoque, Kim and Pyun (2007) and Charles and Darné (2009) for a review.

[START_REF] Lo | Stock market prices do not follow random walk: Evidence from a simple specification test[END_REF] also propose a test statistic under the assumption of homoscedasticity. We focus only on VR statistic which is robust under heteroscedasticity since, as shown in Section 3, all the data display heteroscedasticity.

(T -1)(T + 1)/12r 2,t = Φ -1 r(x)T + 1[START_REF] Davidson | The wild bootstrap, tamed at last[END_REF] Wright (2000) also suggested sign-based tests.

Belaire-Franch and Contreras (2004) also suggested multiple sign-based VR tests based on Wright's sign-based tests but they showed that the rank-based tests are more powerful than the sign-based tests.

The American Petroleum Institute (API) gravity is a measure of how heavy or light a petroleum liquid is compared to water. If its API gravity is greater than

10, it is lighter and floats on water; if less than 10, it is heavier and sinks. API gravity is used to compare the relative densities of petroleum liquids.

We also studied the effect of the Gulf War of 1990 on the VR tests. The outcomes are similar to those on the 1982-1993 period.


Se ´bastien Leclercq

Richard Cordaux
email: richard.cordaux@univ-poitiers.fr

Selection-Driven Extinction Dynamics for Group II Introns in Enterobacteriales

Transposable elements (TEs) are one of the major driving forces of genome evolution, raising the question of the long-term dynamics underlying their evolutionary success. Some TEs were proposed to evolve under a pattern of periodic extinctionsrecolonizations, in which elements recurrently invade and quickly proliferate within their host genomes, then start to disappear until total extinction. Depending on the model, TE extinction is assumed to be driven by purifying selection against colonized host genomes (Sel-DE model) or by saturation of host genomes (Sat-DE model). Bacterial group II introns are suspected to follow an extinction-recolonization model of evolution, but whether they follow Sel-DE or Sat-DE dynamics is not known. Our analysis of almost 200 group II intron copies from 90 sequenced Enterobacteriales genomes confirms their extinction-recolonization dynamics: patchy element distributions among genera and even among strains within genera, acquisition of new group II introns through plasmids or other mobile genetic elements, and evidence for recent proliferations in some genomes. Distributions of recent and past proliferations and of their respective homing sites further provide strong support for the Sel-DE model, suggesting that group II introns are deleterious to their hosts. Overall, our observations emphasize the critical impact of host properties on TE dynamics.

Introduction

Transposable elements (TEs), which are mobilizing pieces of DNA, are widely distributed in eukaryote and prokaryote genomes, and they may represent substantial fractions of the genomes, as in Homo sapiens [START_REF] Lander | Initial sequencing and analysis of the human genome[END_REF] or Zea Mays [START_REF] Schnable | The B73 maize genome: complexity, diversity, and dynamics[END_REF]. TEs are major drivers of genome evolution, sometimes as actors of genetic innovation, or by creating genomic instability or genetic disorders [START_REF] Feschotte | DNA transposons and the evolution of eukaryotic genomes[END_REF][START_REF] Cordaux | The impact of retrotransposons on human genome evolution[END_REF][START_REF] Cerveau | Evolutionary dynamics and genomic impact of prokaryote transposable elements[END_REF]. Given their relative deleteriousness, their tremendous evolutionary success results from a complex interplay between transposition rate, fitness cost, and host effective population size [START_REF] Rouzic | Models of the population genetics of transposable elements[END_REF][START_REF] Rouzic | Long-term evolution of transposable elements[END_REF][START_REF] Goodier | Retrotransposons revisited: the restraint and rehabilitation of parasites[END_REF][START_REF] Yahara | Evolutionary maintenance of selfish homing endonuclease genes in the absence of horizontal transfer[END_REF][START_REF] Boutin | How does selfing affect the dynamics of selfish transposable elements?[END_REF]. In multicellular eukaryotes, the rather low effective population size (compared to unicellular organisms) induces a reduction of selection efficiency and enhanced genetic drift, so that slightly deleterious TE insertions may be maintained [START_REF] Rouzic | Long-term evolution of transposable elements[END_REF][START_REF] Lynch | The origins of genome complexity[END_REF]. Therefore, copies accumulate in nearly neutral genomic regions and undergo gradual degradation, which results in a large amount of old, non functional copies and few but still active elements [START_REF] Lander | Initial sequencing and analysis of the human genome[END_REF][START_REF] Luchetti | Molecular characterization, genomic distribution and evolutionary dynamics of Short INterspersed Elements in the termite genome[END_REF][START_REF] Castoe | Discovery of highly divergent repeat landscapes in snake genomes using high-throughput sequencing[END_REF].

By contrast, prokaryotes usually have much larger effective population sizes, and studies on insertion sequence (IS) TEs suggested that very few TE remnants were present within genomes while most active elements resulted from recent acquisitions and proliferations [START_REF] Lawrence | The evolution of insertion sequences within enteric bacteria[END_REF][START_REF] Touchon | Causes of insertion sequences abundance in prokaryotic genomes[END_REF][START_REF] Wagner | A survey of bacterial insertion sequences using IScan[END_REF]. To explain the evolutionary success of IS elements in prokaryotes, a periodic extinction-recolonization model has been proposed, in which bacterial genomes undergo recurrent TE acquisitions and proliferations, followed by rapid elimination of the resulting TE copies [START_REF] Wagner | Periodic extinctions of transposable elements in bacterial lineages: evidence from intragenomic variation in multiple genomes[END_REF]. This scenario has recently received direct empirical evidence, based on the analysis of the unusual IS fossil record buried in the genomes of the bacterial endosymbiont Wolbachia [START_REF] Cerveau | Short-and long-term evolutionary dynamics of bacterial insertion sequences: insights from Wolbachia endosymbionts[END_REF]. The model is driven by two antagonist properties of IS elements: their transposition rate and their deleteriousness to the host. After acquisition by an ISfree host through horizontal transfer, a novel active IS element proliferates in the host genome due to its high transposition rate. IS copies may then invade the population through multiplication of the colonized cell or IS transmission to neighboring cells through intra-population horizontal transfers. However, genomes harboring IS elements are rapidly removed from the bacterial population because of the deleterious effect of IS insertions, and the bacterial population finally goes back to its initial IS-free state. The second and third steps heavily depend on IS and host genetic properties, and host population sizes [START_REF] Cerveau | Evolutionary dynamics and genomic impact of prokaryote transposable elements[END_REF][START_REF] Cerveau | Short-and long-term evolutionary dynamics of bacterial insertion sequences: insights from Wolbachia endosymbionts[END_REF][START_REF] Nagy | Regulation of transposition in bacteria[END_REF][START_REF] Cordaux | Gene conversion maintains nonfunctional transposable elements in an obligate mutualistic endosymbiont[END_REF], such that most newly acquired IS are probably removed from the population even before having a chance to proliferate. In that way, a high rate of acquisition is critical to ensure IS persistence in prokaryotes [START_REF] Bichsel | The early phase of a bacterial insertion sequence infection[END_REF]. In prokaryotes, TE acquisition is thought to be mainly achieved through larger mobile genetic elements (MGEs), especially plasmids which carry numerous IS elements [START_REF] Leclercq | Do phages efficiently shuttle transposable elements among prokaryotes?[END_REF].

Such dynamics (accumulation in neutral genomic regions or periodic extinctions-recolonizations) are largely driven by the deleteriousness of TEs, and radically different dynamics are expected for TEs which are not constrained by host selective pressures. For instance, group I introns are genetic elements mainly found in housekeeping genes of eukaryote organelles and prokaryotes, which are able to self-splice from the transcript they are located in, making them silent at the translational level [START_REF] Raghavan | Group I introns and inteins: disparate origins but convergent parasitic strategies[END_REF]. Most of them harbor a homing endonuclease gene that allows them to duplicate from an occupied locus to a highly-related intron-free locus (mainly another allele of the same gene) through a process called homing [START_REF] Raghavan | Group I introns and inteins: disparate origins but convergent parasitic strategies[END_REF]. According to their very high site specificity and self-splicing ability, group I intron insertions are likely to be largely neutral, and a model of recurrent invasions and extinctions was proposed for group I introns in organelles and chromosomes of unicellular eukaryotes [START_REF] Goddard | Recurrent invasion and extinction of a selfish gene[END_REF][START_REF] Burt | Homing endonuclease genes: the rise and fall and rise again of a selfish element[END_REF]. In this model, a newly acquired element can spread in the population through homing during mating (when infected and free alleles are together in a unique cell) until fixation (i.e. no free allele available in the population). When the population is saturated, the element cannot home anymore and is slowly degraded through point mutations, to finally lead to static, splicing-only copies. These remnants are sometimes precisely removed, reconstructing a free homing site, which becomes available for future group I invasion.

Group II introns are TEs related to non-long terminal repeat retrotransposons, found in organellar and prokaryote genomes [START_REF] Dai | Compilation and analysis of group II intron insertions in bacterial genomes: evidence for retroelement behavior[END_REF], and the putative ancestors of spliceosomal introns in eukaryotes [START_REF] Martin | Introns and the origin of nucleus-cytosol compartmentalization[END_REF][START_REF] Koonin | Intron-dominated genomes of early ancestors of eukaryotes[END_REF]. Similar to group I introns, group II introns are catalytic RNAs which are able to self-splice from the transcript they are inserted in. The molecular splicing mechanism requires base-pairing interactions between three short motifs located on the intron RNA (EBS1 to 3, for exon-binding sites) and their complementary motifs on the transcript RNA (IBS1 to 3, for intron-binding sites) spanning positions from 212 to +1 relative to the intron insertion site [START_REF] Toro | Bacterial group II introns: not just splicing[END_REF][START_REF] Lambowitz | Group II Introns: Mobile Ribozymes that Invade DNA[END_REF]. Group II introns move via a target-primed reverse transcription mechanism called retrohoming, which is processed by the intron-encoded protein (IEP) [START_REF] Lambowitz | Group II Introns: Mobile Ribozymes that Invade DNA[END_REF][START_REF] Belfort | Mobile Introns: Pathways and Protein[END_REF]. The IEP possesses several catalytic domains necessary for intron mobility, such as reverse-transcriptase (RT) and maturase (X) domains, and sometimes an endonuclease (En) domain. When produced, the IEP binds to the intron RNA and helps in its efficient splicing using the maturase activity, resulting in a free ribonucleoprotein (RNP) particle composed of the intron ribozyme and the IEP. The RNP complex then recognizes a new insertion site through interactions between the IEP and a small number of specific nucleotides in the distal 5'-exon region of the target site, in addition to EBS-IBS binding [START_REF] Singh | Interaction of a group II intron ribonucleoprotein endonuclease with its DNA target site investigated by DNA footprinting and modification interference[END_REF][START_REF] Zhuang | EcI5, a group IIB intron with high retrohoming frequency: DNA target site recognition and use in gene targeting[END_REF][START_REF] Garcı ´a-Rodrı ´guez | Use of RmInt1, a group IIB intron lacking the intron-encoded protein endonuclease domain, in gene targeting[END_REF]. The intron RNA is then reverse-spliced at the integration position using these IBS-EBS interactions, and finally reverse-transcribed by the IEP [START_REF] Lambowitz | Group II Introns: Mobile Ribozymes that Invade DNA[END_REF]. Such a specific targeting ensures integration of group II introns primarily in intron-free alleles of the same gene, although some cases of non-specific retrotransposition have been reported, involving only (sometimes imperfect) IBS motifs [START_REF] Cousineau | Retrotransposition of a bacterial group II intron[END_REF][START_REF] Martı ´nez-Abarca | RecA-independent ectopic transposition in vivo of a bacterial group II intron[END_REF][START_REF] Ichiyanagi | Retrotransposition of the Ll.LtrB group II intron proceeds predominantly via reverse splicing into DNA targets[END_REF]. Some group II introns, called bacterial class C introns, also shift from this general pattern, as they harbor only two EBS (EBS1 and EBS3) and they specifically insert downstream of Rho-independent transcription terminators [START_REF] Lambowitz | Group II Introns: Mobile Ribozymes that Invade DNA[END_REF][START_REF] Robart | Insertion of group II intron retroelements after intrinsic transcriptional terminators[END_REF].

Group II introns are widespread in the bacterial kingdom and show a surprising diversity. They are separated into three ribozyme groups (IIA, IIB, and IIC) and nine ORF classes (A to F, CL1, CL2, ML) which may have diverged for several hundred million years [START_REF] Simon | A broadscale phylogenetic analysis of group II intron RNAs and intron-encoded reverse transcriptases[END_REF]. Moreover, studies conducted on natural populations/species from different bacterial groups revealed an extreme variability in group II intron abundance and diversity between bacterial strains [START_REF] Dai | The dispersal of five group II introns among natural populations of Escherichia coli[END_REF][START_REF] Ferna ´ndez-Lo ´pez | Dispersal and evolution of the Sinorhizobium meliloti group II RmInt1 intron in bacteria that interact with plants[END_REF][START_REF] Tourasse | Survey of group I and group II introns in 29 sequenced genomes of the Bacillus cereus group: insights into their spread and evolution[END_REF][START_REF] Leclercq | Remarkable abundance and evolution of mobile group II introns in Wolbachia bacterial endosymbionts[END_REF][START_REF] Biondi | Spread of the group II intron RmInt1 and its insertion sequence target sites in the plant endosymbiont Sinorhizobium meliloti[END_REF], suggesting recurrent extinctionrecolonization dynamics for these elements. According to group II intron properties (site specificity and self-splicing), we would expect a model of dynamics not constrained by host selection, i.e. rapid homing site saturation in the population, followed by slow sequence degradation and removal, similar to that of eukaryote group I introns. However, several observations suggest that group II introns may not be selectively neutral in bacteria. First, they are virtually never found in housekeeping genes [START_REF] Dai | Compilation and analysis of group II intron insertions in bacterial genomes: evidence for retroelement behavior[END_REF][START_REF] Robart | Group II intron retroelements: function and diversity[END_REF]. Rather, they are preferentially found inserted into regions non essential for the bacterial host, such as plasmids and other mobile genetic elements [START_REF] Tourasse | Survey of group I and group II introns in 29 sequenced genomes of the Bacillus cereus group: insights into their spread and evolution[END_REF][START_REF] Robart | Group II intron retroelements: function and diversity[END_REF][START_REF] Klein | Bacterial group II introns and their association with mobile genetic elements[END_REF][START_REF] Ichiyanagi | A bacterial group II intron favors retrotransposition into plasmid targets[END_REF][START_REF] Staddon | Conserved target for group II intron insertion in relaxase genes of conjugative elements of grampositive bacteria[END_REF]. In addition, some empirical evidence indicates very low efficiency to render viable exons after splicing [START_REF] Chen | Quantitative analysis of group II intron expression and splicing in Lactococcus lactis[END_REF][START_REF] Chillo ´n | Splicing of the Sinorhizobium meliloti RmInt1 group II intron provides evidence of retroelement behavior[END_REF]. Finally, group II introns are sources of genomic instability in some bacteria [START_REF] Leclercq | Remarkable abundance and evolution of mobile group II introns in Wolbachia bacterial endosymbionts[END_REF]. These observations thus suggest that bacterial group II intron dynamics may be governed by host selective pressures, similar to the extinction-recolonization model proposed for IS elements.

Applied to group II introns, these models can be summarized as proposed in Figure 1: Colonization always starts with an acquisition of a novel active element by one cell in an elementfree bacterial population, which then spreads in the host genome and in the population. A selection-driven extinction (Sel-DE) model then predicts that highly colonized genomes are removed from the population through purifying selection (Figure 1A), while a saturation-driven extinction (Sat-DE) model predicts a saturation of all available homing sites in the host population without elimination of highly colonized genomes (Figure 1B). Resulting copies are finally inactivated and degraded until their complete loss. The two models are expected to lead to critical differences in terms of observed group II intron distributions within genomes. Under the Sel-DE model, most genomes should be free of group II introns; and when present, elements should be mostly active and observed at very low copy number. Genomes saturated by active or degraded elements can be rarely observed. Under the Sat-DE model, we expect mainly degraded elements located in all potential homing sites, and probably shared between closely related strains. Few genomes at transient state with a variable number of active copies may also be observed. Finally a number of genomes should also be free of group II introns, depending on the probability of element acquisition.

To identify the evolutionary forces driving group II intron dynamics in bacteria, we analyzed their abundance and distribution in 90 complete genomes from Enterobacteriales (i.e. Escherichia coli and relatives) and compared our observations to expectations under the Sel-DE and Sat-DE models. First, we provide clear examples of plasmid-mediated acquisition and rapid, strainspecific proliferations, consistent with the extinction-recolonization behavior already suspected for group II introns. The distribution of recent and more ancient proliferations also support the Sel-DE model, in which group II intron extinctions are driven by removal of overloaded genomes from the population though purifying selection. Finally, our data highlight the impact of bacterial lifestyle on group II intron abundance and evolution.

Materials and Methods

Genome Sequences

Complete sequences of the 90 Enterobacteriales genomes available in GenBank as of September 14, 2009 were downloaded from the GenBank FTP site (ftp://ftp.ncbi.nih.gov/genomes/Bacteria/all. fna.tar.gz); they are listed in Table S1. Genomes were grouped into 18 clusters according to their genus name, except Shigella and Escherichia strains which belong to the same complex [START_REF] Touchon | Organised genome dynamics in the Escherichia coli species results in highly diverse adaptive paths[END_REF][START_REF] Zhou | Derivation of Escherichia coli O157:H7 from its O55:H7 precursor[END_REF][START_REF] Sahl | A comparative genomic analysis of diverse clonal types of enterotoxigenic Escherichia coli reveals pathovar-specific conservation[END_REF], and Erwinia and Pectobacterium strains which are also included in the same complex [START_REF] Samson | Transfer of Pectobacterium chrysanthemi (Burkholder et al. 1953) Brenner et al. 1973 and Brenneria paradisiaca to the genus Dickeya gen. nov. as Dickeya chrysanthemi comb. nov. and Dickeya paradisiaca comb. nov. and delineation of four novel species, Dickeya dadantii sp. nov., Dickeya dianthicola sp. nov., Dickeya dieffenbachiae sp. nov. and Dickeya zeae sp. nov[END_REF][START_REF] Kube | The genome of Erwinia tasmaniensis strain Et1/99, a non-pathogenic bacterium in the genus Erwinia[END_REF].

Group II Intron Detection

We downloaded the 397 reference sequences of putatively functional introns available in Zimmerly's group II intron database (update July 29, 2011), hereafter referred to as Zbase [START_REF] Candales | Database for bacterial group II introns[END_REF], and constructed nucleotide and protein libraries. The nucleotide library contained whole group II intron sequences while the protein library contained the amino acid residues encoded by their open reading frame (ORF). We complemented these datasets with sequences of full-length but putatively non-functional introns that were available in a previous Zbase version (update March 11, 2008). BLAST searches were then performed for nucleotide and protein intron libraries against Enterobacteriales genomes (e-value ,0.01), and results were filtered for overlapping detections to keep hits with the highest number of identical sites (nucleotides or amino acid residues). Detections shorter than 100 bp in length or less than 75% similarity to the reference sequence (for BLASTN only) were discarded to avoid false positives. Group II intron copies were categorized as full-length when their extremities could be unambiguously defined, and fragmented otherwise. For copies that did not cover the whole sequence of their intron reference, boundaries were searched by comparing flanking sequences of multiple instances of the same intron (when available), or by looking for an intron-free locus homologous to the insertion site in GenBank. A second round of BLASTN searches using all previously detected copies (full-length and fragmented copies) was performed to search for more divergent copies. Intron names were attributed according to the Zbase nomenclature for copies showing more than 10% nucleotide divergence compared to their reference intron sequence, and with the Zbase reference intron name otherwise.

Determination of Insertion Patterns and Orthology Relationships

Insertion patterns and orthology relationships were evaluated using the following comparison procedure. The immediate 3 kbp flanking sequences of each intron copy were searched against other genomes of the same genus/complex using BLASTN. When the 59 and 39 flanks aligned to directly adjacent regions in another genome, the copy was considered to result from retrohoming in the query genome. When 59 and/or 39 flanks of an intron copy aligned with the 59 and/or 39 flanks of an intron copy in another genome, copies were considered as orthologous, with a few exceptions (see below). For analyses of abundance of distinct Individuals with few copies ultimately disappear from the population through rampant purifying selection or eventually genetic drift, resulting in complete loss of the group II intron from the bacterial population (step 4). Alternatively, not proliferated copies may be inactivated by mutation (step 5), evolve as neutral genomic regions and eventually become fixed in the population (step 6). Residual inactivated copies then undergo slow degradation until complete loss from the bacterial genome (step 7). (B) Saturation-driven extinction (Sat-DE) model: the first two steps are identical to those of the Sel-DE model, but proliferation results in the saturation of all available homing sites in the population (step 3). Proliferated copies then start to degenerate (step 4) and are slowly removed from the population through random mutations and deletions (step 5). Bacterial populations are represented by sets of cells (rounded rectangles) harboring available homing sites (grey boxes) or homing sites occupied by active (green boxes) or inactive intron copies (dark red boxes). The starred state in each model is those predicted to be the most frequently observed when a group II intron copy is present. Smaller cartoon sizes (in A) reflect infrequent alternative. doi:10.1371/journal.pone.0052268.g001 elements, orthologous copies were counted as only one element. Flanking sequences of many group II intron copies were not readily identifiable using this procedure, mainly because they were parts of larger MGEs, such as integrated plasmids, genomic islands (GIs), or prophages. For these intron copies, the surrounding genomic region (up to 1 Mbp) was compared visually to closely related genomes, using the MaGe synteny browser [START_REF] Vallenet | MaGe: a microbial genome annotation system supported by synteny results[END_REF]. In most cases, this analysis revealed specific insertions of large genetic sequences (.10 Kbp), which were considered as integrated MGEs. During this analysis, we found that the SelC intron-carrying GI inserted independently at two different genomic loci in two E. coli lineages (CFT073/ED1a and S88/APEC_01). Thus, the group II intron copies carried by these two GIs were considered as distinct, while they were initially considered as orthologous with the 3 Kbp flanking sequence comparison (see above). Two other cases of shared but not vertically inherited GIs were detected: between the 55989 E. coli strain and the Shigella flexneri lineage, and between the two E. coli strains ED1a and 536. Group II intron copies inserted in these GIs were consistently counted as distinct copies. Insertion patterns, positions and orthology relationships are shown in Table S2 and Table S3.

Detection of Putative Available Homing Sites

Four major schemes were conducted to detect potential group II intron homing sites. For introns which retrohome into specific IS elements (E.c.I1, E.c.I3, E.c.I4, E.c.I9, and c-Ha.de.I1 subgroup 2), the 90 bp regions surrounding homing sites were searched within genome sequences using a BLASTN procedure. Hits with a evalue #0.0001, a similarity $90% with the reference 90 pb, and that overlapped with the precise integration site were counted as putative homing sites. For other already referenced introns (Kl.pn.I5, Di.ze.I1, and So.gl.I1), putative EBS motifs were retrieved from the secondary structure available on the current (http:// webapps2.ucalgary.ca/ ˜groupii/index.html) or previous (http:// www.fw.ucalgary.ca/group2introns/index.htm) versions of Zbase. Consensus IBS sequences were constructed by aligning the predicted EBS motifs to the 59 and 39 regions surrounding detected full-length copies. For Kl.pn.I5, sequence surrounding the Zbase copy was also added to the analysis, and all identical bases spanning from positions 225 to +10 relative to the insertion site were added to the consensus homing site sequence. Consensus sequences were then searched within genomes using a Perl script according to standard base-pairing rules with perfect match. For introns belonging to bacterial ORF class C (Di.ze.I1 and So.gl.I1), host genomes were searched for Rho-independent transcription terminators with the ARNold algorithm [START_REF] Naville | ARNold: a web tool for the prediction of Rho-independent transcription terminators[END_REF]. Putative IBS1/3 located 3' to a transcription terminator stem-loop were considered as potential homing sites. Finally, for c-Ha.de.I1 subgroup 1 for which EBS motifs are unknown, conserved bases flanking the 26 full-length copies were estimated using Weblogo 3.3 [START_REF] Crooks | WebLogo: a sequence logo generator[END_REF]. The estimated consensus sequence was then searched within the H. defensa genome using the same Perl script as above.

Computation of Intra-genus Mean Nucleotide Divergence

We estimated intra-genus mean nucleotide divergence for the Escherichia/Shigella complex and the Salmonella genus as a measure of genus diversity given by the sequenced strains. For each genus/ complex, we searched all genomes for the housekeeping genes used for MLST analyses, according to the MLST database (http://mlst.ucc.ie/): adk, fumC, gyrB, icd, mdh, purA and recA for Escherichia/Shigella genomes, and aroC, dnaN, hemD, hisD, purE, sucA and thrA for Salmonella genomes. Gene sequences were concatenated and aligned with the L-INS-i method of MAFFT [START_REF] Katoh | MAFFT version 5: improvement in accuracy of multiple sequence alignment[END_REF], and mean nucleotide divergence was computed for each genus/ complex with MEGA 5 [START_REF] Tamura | MEGA5: Molecular Evolutionary Genetics Analysis Using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods[END_REF] using third codon positions only. Genomes of E. fergusonii and S. enterica serovar arizonae can be considered as outgroups compared to other Escherichia/Shigella and Salmonella genomes respectively [START_REF] Sahl | A comparative genomic analysis of diverse clonal types of enterotoxigenic Escherichia coli reveals pathovar-specific conservation[END_REF][START_REF] Fricke | Comparative genomics of 28 Salmonella enterica isolates: evidence for CRISPRmediated adaptive sublineage evolution[END_REF], and were excluded from the analysis.

Results

General Abundance and Distribution of Group II Introns in Enterobacteriales

We searched 90 sequenced genomes of Enterobacteriales for group II intron elements and found a total of 198 copies longer than 100 bp (see Methods), 87 (44%) of which are full-length elements. We found that 132 (67%) copies are inserted in chromosomes and 66 (33%) in plasmids. Although it represents a density of 2.2 group II intron copies per genome on average, only 50 (56%) of the 90 analyzed genomes carry at least one group II intron (Figure 2A). Moreover, group II introns are not uniformly distributed among the different bacterial genera, and their abundance does not seem to be related to the number of sequenced genomes within genera (Figure 2B). Indeed, there was no significant relationship between the number of sequenced strains per genus/complex and the abundance of group II introns (Chi square test, P = 0.41). We searched for orthologous copies within genera/complexes to prevent counting several times copies that originated from a single ancestral insertion event, which would artificially inflate the overall abundance of group II introns in genera with more than one sequenced strain (see Materials and Methods). A total of 90 copies were found to be orthologous between at least two genomes within genera, reducing the total number of detections to 135 distinct intron copies (i.e. when orthologues are counted as a single copy; Table 1). Again, the relationship between the number of sequenced strains per genus/complex and the abundance of distinct copies was not significant (Chi square test, P = 0.23).

Among the four most sequenced genera (i.e. .5 sequenced strains), only the Escherichia/Shigella complex shows group II intron abundance exceeding two distinct copies (Figure 2B). The absence of group II intron in the six Buchnera aphidicola strains is expected, as these bacteria are ancient mutualistic endosymbionts with highly reduced genomes known to be completely devoid of TEs [START_REF] Toft | Evolutionary microbial genomics: insights into bacterial host adaptation[END_REF]. Concerning the Yersinia genus, 11 of the 12 sequenced strains belong to the very closely related and almost monomorphic Y. pestis and Y. pseudotuberculosis species [START_REF] Achtman | Yersinia pestis, the cause of plague, is a recently emerged clone of Yersinia pseudotuberculosis[END_REF][START_REF] Laukkanen-Ninios | Population structure of the Yersinia pseudotuberculosis complex according to multilocus sequence typing[END_REF]. These data therefore probably display only a small subset of the whole Yersinia diversity, which may account for the apparent paucity of group II intron in this genus. Sequenced Salmonella strains also represent a small subset of the whole genus diversity, as all strains (except S. enterica arizonae) belong to the S. enterica group I [START_REF] Fricke | Comparative genomics of 28 Salmonella enterica isolates: evidence for CRISPRmediated adaptive sublineage evolution[END_REF]. However, the sequenced Salmonella strains show a genetic diversity nearly equivalent to that of Escherichia/Shigella strains (mean nucleotide divergence on MLST genes of 0.027 and 0.038, respectively), and they have experienced several DNA acquisitions through plasmids and prophages since their divergence [START_REF] Fricke | Comparative genomics of 28 Salmonella enterica isolates: evidence for CRISPRmediated adaptive sublineage evolution[END_REF]. The complete lack of group II intron detection in the 16 Salmonella genomes may therefore be attributed to the erratic nature of extinctionrecolonization processes [START_REF] Bichsel | The early phase of a bacterial insertion sequence infection[END_REF].

Group II Intron Acquisition through MGE-mediated Horizontal Transfers

In both Sel-DE and Sat-DE models, horizontal transfers of group II intron copies between bacterial cells are essential to ensure periodic recolonizations. Such transfers are believed to occur via large MGEs such as plasmids or GIs [START_REF] Tourasse | Survey of group I and group II introns in 29 sequenced genomes of the Bacillus cereus group: insights into their spread and evolution[END_REF][START_REF] Klein | Bacterial group II introns and their association with mobile genetic elements[END_REF][START_REF] Ichiyanagi | A bacterial group II intron favors retrotransposition into plasmid targets[END_REF][START_REF] Staddon | Conserved target for group II intron insertion in relaxase genes of conjugative elements of grampositive bacteria[END_REF]. *: Two orthologous copies in the Escherichia/Shigella complex are full-length in some strains and fragmented in others; they are counted only as distinct full-length copies. doi:10.1371/journal.pone.0052268.t001

Genomic locations of group II introns in the Escherichia/Shigella complex support a high horizontal transfer potential, as 54 (81%) of the 67 distinct Escherichia/Shigella group II intron copies are located in free or integrated MGEs (Table S2). Although some intron copies may have inserted in MGEs after MGE acquisition by host bacterial cells, instead of being shuttled by them, we identified at least three highly probable cases of group II intron horizontal transfers via MGEs in Escherichia/Shigella strains. The first transfer concerns an E.c.I3 copy in E. coli strains ED1a and 536 (unique E.c.I3 copy in this strain), which is located at an orthologous position in the same GI (PAI1), but this GI is located at unrelated positions in the two bacterial strains. Similarly, E.c.I6 is present only once in E. coli strains CFT073, ED1a, S88, and APEC_O1, and is inserted at an orthologous site in a GI inserted at different loci in the various strains. Finally, the SHI-1 GI carries a E.c.I4 intron copy at the same position in the two Shigella flexneri 2a strains (the only full-length copy in these strains) and in E. coli 55989, but the GI is again inserted at different loci. These examples show group II intron acquisitions following chromosomal integration of the host MGE. In addition, other Enterobacteriales genomes provide us with strong support for chromosomal acquisition through another mechanism, namely MGE-to-chromosome retrohoming. For example, there are two full-length Kl.pn.I5 intron copies in the Klebsiella pneumoniae MGH genome: one located in a plasmid and the other one showing evidence of retrohoming in the chromosome. Both intron copies are identical at the nucleotide level, suggesting that the chromosomal copy was acquired from the plasmid through retrohoming (although a transfer from the chromosome to the plasmid cannot be formally ruled out). Dickeya genomes provide another example of MGE-to-chromosome intron acquisition. One of these genomes carries two Di.ze.I1 copies, which are identical and were acquired through retrohoming (as both homing sites are free in the other Dickeya genome). As no other Di.ze.I1 copy was detected in this genome, observed copies were most probably acquired from a non-integrated MGE which has since disappeared.

Homing Site Occupation of Newly Acquired, Potentially Active Group II Introns

Following acquisition, both models predict that incoming elements start to proliferate, sometimes quite rapidly for group II introns [START_REF] Leclercq | Remarkable abundance and evolution of mobile group II introns in Wolbachia bacterial endosymbionts[END_REF][START_REF] ´nez | Dispersion of the RmInt1 group II intron in the Sinorhizobium meliloti genome upon acquisition by conjugative transfer[END_REF]. Two main factors limit group II intron proliferation: the number of available homing sites and, in the Sel-DE model, the power of selection. Thus, genomes with saturated or almost saturated homing sites should be commonly observed under the Sat-DE hypothesis, while they should be very rarely observed under the Sel-DE hypothesis (except when the number of homing sites is low).

Our dataset contains 87 full-length copies from 14 group II introns distributed among 22 genomes (Table 2). For five of them (E.c.I5, Di.da.I1, P.l.I1, P.l.I2, and c-Ha.de.I2), only one distinct copy is detected and their homing site have never been described. We were thus unable to infer the number of potential homing sites for these introns. Secondly, when checking intron sequences with ORFinder (http://www.ncbi.nlm.nih.gov/ projects/gorf/) for the remaining elements, we found that E.c.I4 copies carried by E. coli SMS_3_5 and S. flexneri 2a 2457T, as well as the four E.c.I6 copies, exhibit a disrupted ORF. As all of them are present in free or integrated MGE (Table S2), they may have been acquired by the host genome already inactivated and unable to proliferate. Evaluating the proportion of occupation for these introns may thus not be relevant and was not conducted.

Among the remaining, potentially active group II intron copies, E.c.I2, E.c.I9, E.c.I3, and E.c.I4 copies detected in Escherichia/Shigella genomes are known to specifically target Rsh IS ( = ISEc125), Rsh IS at another position, IS679, and various IS3 elements, respectively [START_REF] Dai | The dispersal of five group II introns among natural populations of Escherichia coli[END_REF][START_REF] Han | Transposition of IntAs into the conserved regions of IS3 family elements[END_REF]. Both E.c.I2 and E.c.I9 are found only once in their respective genomes, in which they occupy only 20% and 14% of the potential homing sites, respectively (Table 2). E.c.I3 is also detected at low copy number per genome (2 and 1 copies), but no free IS679 was detected in colonized genomes (Table 2). On the contrary, E.c.I4 shows a large range of abundance, from 1 copy in three genomes up to 15 copies in the genome of IAI39 (Table 2). All copies show ,1% intra-strain nucleotide divergence and are strain-specific (except one shared by E24377A and 55989 strains), indicating that the E.c.I4 intron experienced very recent and independent proliferations in these E. coli strains. As expected, all copies are inserted at specific positions in IS3-like transposases, namely from IS629, IS911, ISEc16, and ISEc31 elements. Searches of these four potential host genes in E. coli/ Shigella genomes revealed that genomes of 55989 and E24377A strains are saturated or almost saturated, and 67% of the potential homing sites are occupied in ED1a (Table 2). By contrast, IAI39, UTI89, and UMN026 have less than 40% occupied homing sites, and only one copy is present in S. flexneri Sf301 despite 42 potentially available homing sites (Table 2).

Homing site occupation can also be estimated for potentially active copies of Kl.pn.I5 (in Klebsiella pneumoniae MGH), Di.ze.I1 (in Dickeya zeae Ech1591), So.gl.I1 (in Sodalis glossinidius 'morsitans'), and c-Ha.de.I1 (in Hamiltonella defensa 5AT). The two Kl.pn.I5 copies are identical and inserted in intergenic regions harboring well conserved IBS sequences, and comparison with the 225 to +10 flanking sequences of the reference Kl.pn.I5 (located in the K. pneumoniae plasmid pK245, cf. Zbase) extends the putative homing site to some additional bases (Figure 3A). When searching the MGH genome for available homing sites, no occurrence was detected, indicating an homing site saturation (Table 2). Di.ze.I1 and So.gl.I1 are bacterial class C introns and detected copies are located downstream of Rhoindependent transcription terminators, as expected (Figure 3B). Target regions also show well conserved IBS1 and IBS3 motifs, indicating that these detected copies resulted from retrohoming. The search for Rho-independent transcription terminators followed by consensus IBS1/3 sequences in D. zeae and S. glossinidius genomes yielded a homing site occupation of 11% (2/ 17) and 4% (3/73) for Di.ze.I1 and So.gl.I1 introns, respectively (Table 2). c-Ha.de.I1 in H. defensa is represented by 29 full-length copies, for which a previous analysis reported an average nucleotide divergence ,2% [START_REF] Degnan | Hamiltonella defensa, genome evolution of protective bacterial endosymbiont from pathogenic ancestors[END_REF]. Our inspection of c-Ha.de.I1 sequences revealed that the pool of copies is distributed into two subgroups separated by 6% average divergence. Subgroup 1 includes 26 copies showing an average nucleotide divergence of 0.1%. Subgroup 2 encompasses two identical copies and a third one with 2.5% divergence, largely due to a ,700 bp-long gene conversion event (involving a subgroup 1 copy) which artificially inflates the genetic distance with the two other subgroup 2 copies (Figure S1). The two subgroups thus probably result from two recent, independent proliferation events. Inspection of regions surrounding insertion sites revealed that c-Ha.de.I1 copies from subgroup 2 are all inserted in ISPlu15-like IS elements at a specific position. No intron-free occurrence of this IS was found in the chromosomal genome, but one degraded fragment including the homing site was detected in the plasmid (Table 2). By contrast, copies from subgroup 1 show a variety of genomic insertion sites. For example, one copy is inserted in another type of ISPlu15-like IS element (at the same position though), another copy is inserted in an IS427-like element, and three copies are inserted in other copies of c-Ha.de.I1 (at a specific position). Some conserved positions can however be detected within 225 to +15 intron flanking sequences, which were used to infer a consensus homing site (Figure 3C). According to this consensus, 17 homing sites (excluding uninterrupted intron copies) are still available in the H. defensa genome analyzed here (Table 2). Thus, c-Ha.de.I1 subgroup 1 occupies only 60% of its potential homing sites in H. defensa.

Overall, among the 16 cases in which proliferation of a potentially active group II intron can be estimated, only 5 exhibit a homing site occupation larger than 80%, and in all cases the number of potential homing sites was lower than 10. By contrast, none of the five genomes which harbored more than 10 potential homing sites for recently acquired group II introns are saturated. These results mainly support the Sel-DE model, in which genomes overloaded by group II intron copies are rapidly eliminated from bacterial populations.

Fate of Past Group II Intron Acquisitions and Proliferations

The section above provided homing site occupation for recently acquired and still potentially active group II introns. The fate of these proliferations on the long term is a critical difference between the Sat-DE and Sel-DE models (Figure 1). The Sat-DE model predicts that ancient proliferations may be visible for a while in bacterial genomes as degraded copies spanning all potential homing sites. On the contrary, the Sel-DE model predicts that ancient proliferations should not be observable because of the selection pressure acting on intron-loaded genomes. The Sel-DE model also predicts that intron copies may sometimes escape the selection filter when they are inactivated by mutation before proliferation.

One hundred eleven group II intron fragments were detected in our analysis, distributed in 55 distinct copies (Table 1). Most of them are only distantly related to a described active group II intron (Table S2), and their potential homing site cannot be inferred. We thus restricted our analysis to group II intron fragments related to full-length elements described in the previous section and from which we have information on putative homing sites. Two distinct E.c.I9 fragments were detected in 7 and 3 E. coli genomes, respectively. All these genomes also harbor several free Rsh IS copies, leading to a maximum homing site occupation of 33% (Table 3). Four distinct E.c.I3 fragments were detected, one on a E. coli E24377A plasmid and all others located on the Shigella virulence plasmid. Again, putative available homing sites were detected in host genomes, with 6 occurrences in E. coli E24377A and one in each S. boydii genome (Table 3). Four distinct E.c.I4 fragments were also detected in the Escherichia/Shigella genus.

Putative homing site occupation is #50% in all E. coli strains, and ,10% in all Shigella strains (Table 3). However, similar to inactivated full-length E.c.I6, all E.c.I4 fragments are located in large MGEs (Table S2), and we cannot exclude a scenario in which they have been acquired already inactivated and unable to proliferate. Two fragments of So.gl.I1 were detected in the S. glossinidius genome. When these copies were inactivated (that we assume before more recent copies were acquired), 76 homing sites still remained available (Table 3). Finally, four fragments of c-Ha.de.I1 were detected in the H. defensa genome, although three of them show 100% identity with copies from subgroup 1 and may not result from past proliferations. The last one is also related to subgroup 1 copies, but with 4% divergence. If we assume that its homing site was the same as those of subgroup 1 copies, this fragment left at least 98% (48/49) of its potential homing sites free before being inactivated (Table 3). In summary, despite numerous group II intron fragments detected in our dataset, no evidence of any large ancient proliferation that led to homing site saturation can be identified. On the contrary, detected fragments generally show no counterpart belonging to the same intron in the same genome, despite the presence of putative available homing sites. This suggests that these copies were inactivated before their proliferation, or that genomes in which proliferation occurred didn't survive the purifying selection filter. All these observations are more consistent with a Sel-DE model than with a Sat-DE model for explaining group II intron dynamics.

Discussion

Group II intron distribution in bacteria is known to be highly variable in abundance and diversity, suggesting periodic extinction-recolonization events [START_REF] Ferna ´ndez-Lo ´pez | Dispersal and evolution of the Sinorhizobium meliloti group II RmInt1 intron in bacteria that interact with plants[END_REF][START_REF] Tourasse | Survey of group I and group II introns in 29 sequenced genomes of the Bacillus cereus group: insights into their spread and evolution[END_REF][START_REF] Leclercq | Remarkable abundance and evolution of mobile group II introns in Wolbachia bacterial endosymbionts[END_REF]. However, whether extinctions are driven by purifying selection acting on group II intron proliferation (Sel-DE model) or by slow degradation of elements after genome saturation (Sat-DE model) was still unclear. Our analysis of group II introns in Enterobacteriales first confirms their heterogeneous distribution between and within bacterial species, and their tendency to be located in larger MGEs which may favor their horizontal transfer, consistent with an extinction-recolonization model. Next, introns subject to recent (or ongoing) and past proliferation events have not reached homing site saturation in genomes they colonized, except when the number of homing sites was very low. The Sat-De model assumes that there is no other barrier to proliferation than the number of available homing sites when a new group II intron enters a genome. On the contrary, the Sel-DE model assumes that the more loaded with group II introns copies the genomes are, the more they are counter-selected, even if not all homing sites are filled. Our observations thus better support group II intron dynamics constrained by selective forces limiting intra-genomic proliferations.

Deleteriousness of inserted elements is a prerequisite of the Sel-DE model, because it conditions their quick removal from the bacterial populations [START_REF] Wagner | Periodic extinctions of transposable elements in bacterial lineages: evidence from intragenomic variation in multiple genomes[END_REF][START_REF] Bichsel | The early phase of a bacterial insertion sequence infection[END_REF]. Individual TE copies may be deleterious when they insert within genes or regulatory regions [START_REF] Cerveau | Evolutionary dynamics and genomic impact of prokaryote transposable elements[END_REF][START_REF] ´nais | The impact of transposable elements on eukaryotic genomes: From genome size increase to genetic adaptation to stressful environments[END_REF]. In principle, group II introns inserted in transcribed regions have the ability to precisely self-splice from the RNA transcript, thereby allowing the correct expression of the genomic region they are inserted in [START_REF] Toro | Bacterial group II introns: not just splicing[END_REF]. Moreover, this study and others show that group II introns are generally located within MGEs or other non essential genomic regions [START_REF] Dai | Compilation and analysis of group II intron insertions in bacterial genomes: evidence for retroelement behavior[END_REF][START_REF] Tourasse | Survey of group I and group II introns in 29 sequenced genomes of the Bacillus cereus group: insights into their spread and evolution[END_REF][START_REF] Robart | Group II intron retroelements: function and diversity[END_REF][START_REF] Mohr | Mechanisms Used for Genomic Proliferation by Thermophilic Group II Introns[END_REF]. Sel-DE dynamics for these elements may therefore be somewhat unexpected, as intron insertions would be expected to have no or limited effect on bacterial fitness. Two explanations can resolve this apparent contradiction. First, group II introns may be found only in non-essential genomic regions because copies inserted in essential genes are so deleterious that they are quickly removed from the population and virtually undetectable. This would imply that the splicing efficiency of group II introns should be called into question, or at least their efficiency to render viable exons after splicing. Consistently, it was observed that introns located in an IS element and in a gene only essential for plasmid propagation were respectively very and moderately inefficient for reconstructing the interrupted exon after splicing [START_REF] Chen | Quantitative analysis of group II intron expression and splicing in Lactococcus lactis[END_REF][START_REF] Chillo ´n | Splicing of the Sinorhizobium meliloti RmInt1 group II intron provides evidence of retroelement behavior[END_REF]. On the other hand, an active group II intron was detected in the essential recA gene of Geobacillus kaustophilus, in which it seems not to alter host survival [START_REF] Chee | Housekeeping recA gene interrupted by group II intron in the thermophilic Geobacillus kaustophilus[END_REF]. This intron shows hallmarks of host adaptation (based on amelioration of nucleotide composition), and it may also have increased its exon-ligation efficiency to lower its impact on host fitness. The question of exon-ligation efficiency after splicing is thus still in great debate and clearly needs further experimental investigation.

Alternatively, group II intron deleterious effect may be caused by the accumulation of copies rather than by the deleteriousness of individual insertions. Repetitive elements, particularly TEs, are known to be factors of genomic instability because they are preferred targets of ectopic homologous recombination events [START_REF] ´nais | The impact of transposable elements on eukaryotic genomes: From genome size increase to genetic adaptation to stressful environments[END_REF][START_REF] Treangen | Genesis, effects and fates of repeats in prokaryotic genomes[END_REF][START_REF] Santoyo | Gene conversion and concerted evolution in bacterial genomes[END_REF], and this holds true for group II introns as well [START_REF] Leclercq | Remarkable abundance and evolution of mobile group II introns in Wolbachia bacterial endosymbionts[END_REF]. Although TE-mediated genomic rearrangements may occasionally result in genetic innovations [START_REF] Cerveau | Evolutionary dynamics and genomic impact of prokaryote transposable elements[END_REF][START_REF] Hua-Van | The struggle for life of the genome's selfish architects[END_REF], most events are structurally deleterious and subject to purifying selection. A large number of homologous TEs within a genome is thus by nature counterselected, as it mechanistically increases the probability of genomic rearrangements. Increasing the number of active copies also increases the overall activity of a given element. For group II introns, it would result in an accumulation of RNP particles wandering for available homing sites, increasing as much the probability of retrotransposition (non-specific insertions), sometimes in important genes.

The impact of negative selection on group II intron proliferation probably shaped the diversity of homing site targeting strategies. For instance, bacterial class C introns target downstream regions of transcription terminators, which means that they are never transcribed except when the transcription machinery miss the end checkpoint. This kind of specific targeting can thus be viewed as a self-regulation strategy to limit intron proliferation [START_REF] Dai | Compilation and analysis of group II intron insertions in bacterial genomes: evidence for retroelement behavior[END_REF][START_REF] Robart | Insertion of group II intron retroelements after intrinsic transcriptional terminators[END_REF]. Our data perfectly illustrate such a self-regulation process, with a very low number of Di.ze.I1 and So.gl.I1 detected in D. zeae and S. glossinidius genomes, respectively, despite numerous potential homing sites (Table 2). Another strategy widely used by group II introns is to target IS elements [START_REF] Dai | The dispersal of five group II introns among natural populations of Escherichia coli[END_REF][START_REF] Biondi | Spread of the group II intron RmInt1 and its insertion sequence target sites in the plant endosymbiont Sinorhizobium meliloti[END_REF][START_REF] Klein | Bacterial group II introns and their association with mobile genetic elements[END_REF]. This homing strategy provides several advantages, such as numerous putative homing sites, located in non essential genes, and which may help for horizontal transfers. Moreover, IS element abundance itself is limited by selection [START_REF] Wagner | Periodic extinctions of transposable elements in bacterial lineages: evidence from intragenomic variation in multiple genomes[END_REF], which means that intron proliferations are de facto controlled. Surprisingly, among the nine Escherichia/ Shigella genomes harboring the IS-targeting E.c.I4, only three show intron saturation (Table 2). Although it may be caused by intron inactivity (at least for genomes with only one full-length representative), it may also result from a cumulative effect of repeated elements on genome instability. Indeed, intron/IS tandems span .3 kb in length, while IS alone are generally ,1 kb-long. As recombination frequency increases with the size and similarity of the homologous region [START_REF] Petes | Recombination between repeated genes in microorganisms[END_REF][START_REF] Matic | Genetic barriers among bacteria[END_REF], IS interrupted by group II introns may be more prone than IS or introns alone to induce genomic rearrangements. Thus, group II intron proliferation within other repeated elements (such as IS) may be even more efficiently counter-selected than proliferation into less specific target regions.

Bacterial lifestyle may also play a critical role in group II intron proliferation. Reduction of effective population size of bacteria is generally linked with severe proliferations of IS elements, because of a weakening of the efficiency of purifying selection [START_REF] Moran | Genomic changes following host restriction in bacteria[END_REF][START_REF] Escobar-Pa ´ramo | Evidence for genetic drift in the diversification of a geographically isolated population of the hyperthermophilic archaeon Pyrococcus[END_REF][START_REF] Mira | The Neolithic revolution of bacterial genomes[END_REF]. This process is particularly striking in intracellular bacteria: genomes of facultative intracellular or recent obligate endosym-bionts tend to exhibit larger IS abundance compared to those of free-living species, which generally contain very few IS elements [START_REF] Toft | Evolutionary microbial genomics: insights into bacterial host adaptation[END_REF][START_REF] Moran | Genomic changes following host restriction in bacteria[END_REF]. Interestingly, two of the facultative endosymbionts in our dataset (H. defensa and S. glossinidius), have experienced recent group II intron proliferations. However, reduction of the efficiency of purifying selection only allows the proliferation of elements that were already present and functional in the genome of these bacteria. This is fairly well illustrated by the facultative intracellular bacterium Y. pestis, in which all IS families that have proliferated were present in the free-living ancestor Y. pseudotuberculosis [START_REF] Chain | Complete genome sequence of Yersinia pestis strains Antiqua and Nepal516: evidence of gene reduction in an emerging pathogen[END_REF]. Thus, at least one active group II intron should have been present in the genome of bacteria that shifted toward intracellularity to benefit from the reduction of effective population size and proliferate. This may explain why the three young and facultative intracellular species P. asymbiotica, S flexneri, and Y. pestis do not show any evidence for recent group II intron proliferation.

Intracellularity may also have a further impact on Sel-DE dynamics: such bacteria are confined within host cells, which reduces the opportunities for MGE exchange with other bacteria. For example, obligate mutualistic bacteria, which are transmitted essentially vertically and are clustered within specific cellular compartments, show virtually no DNA acquisition [START_REF] Toft | Evolutionary microbial genomics: insights into bacterial host adaptation[END_REF]. As these bacteria no longer exchange plasmids or phages, they are not subject to the first phase of the Sel-DE model, and are therefore ''protected'' against new group II intron colonizations. In our dataset, nine genomes are from obligate mutualistic endosymbionts (six Buchnera aphidicola strains, two candidatus Blochmania spp. strains, and one Wiggleworthia glossinidia strain). As expected, no group II intron was detected in their genomes, consistent with the view that they have eliminated copies of pre-mutualistic stage proliferations, and that their strict intracellular confinement mainly protects them against new invasions. Group I intron dynamics in organelles follow a saturationdriven extinction (Sat-DE) model, and it was proposed that this model could fit to other target-primed retroelements in other organisms [START_REF] Burt | Homing endonuclease genes: the rise and fall and rise again of a selfish element[END_REF]. On the contrary, our study indicates that group II introns evolve under a selection-driven extinction (Sel-DE) model initially proposed to describe the evolution of IS elements in bacterial genomes [START_REF] Wagner | Periodic extinctions of transposable elements in bacterial lineages: evidence from intragenomic variation in multiple genomes[END_REF]. It emphasizes that Sel-DE dynamics are mostly constrained by life history traits of their bacterial hosts rather than by TE intrinsic properties, although TE characteristics (such as specific or non-specific target insertion sites) may play a role on their relative abundance within genomes. This could explain why IS elements are generally much more common than group II introns in bacterial genomes. Effective population size as well as lateral gene transfers have also critical effects on these dynamics. Most bacteria show two of the prerequisites for Sel-DE dynamics, i.e. large population sizes and frequent genetic exchanges, contrary to eukaryotes which generally show smaller effective population sizes and lower rates of lateral gene transfer, and for which other TE evolutionary strategies were observed [START_REF] Hua-Van | The struggle for life of the genome's selfish architects[END_REF]. Yet, the Sel-DE model may also apply to some extent to eukaryote TEs such as mariner elements, which experience frequent horizontal transfers in their life cycles, followed by rapid proliferation and degeneration in host genomes [START_REF] Robertson | Evolution of DNA Transposons in Eukaryotes[END_REF][START_REF] Lampe | Recent horizontal transfer of mellifera subfamily mariner transposons into insect lineages representing four different orders shows that selection acts only during horizontal transfer[END_REF]. Whether Sel-DE dynamics apply to other eukaryote TEs is an open and essential question, which could be investigated for example by looking at TE dynamics in eukaryotes sharing properties similar to those of bacteria (e.g. unicellular free living eukaryotes).

Figure 1 .

 1 Figure 1. Two proposed extinction-recolonization dynamics of group II introns in a bacterial population. (A) Selection-driven extinction (Sel-DE) model: acquisition of a functional copy via a MGE-mediated horizontal transfer (step 1, red arrow), quick element proliferation in the host genome and in the population (step 2), and removal of bacterial individuals with proliferated copies through purifying selection (step 3).Individuals with few copies ultimately disappear from the population through rampant purifying selection or eventually genetic drift, resulting in complete loss of the group II intron from the bacterial population (step 4). Alternatively, not proliferated copies may be inactivated by mutation (step 5), evolve as neutral genomic regions and eventually become fixed in the population (step 6). Residual inactivated copies then undergo slow degradation until complete loss from the bacterial genome (step 7). (B) Saturation-driven extinction (Sat-DE) model: the first two steps are identical to those of the Sel-DE model, but proliferation results in the saturation of all available homing sites in the population (step 3). Proliferated copies then start to degenerate (step 4) and are slowly removed from the population through random mutations and deletions (step 5). Bacterial populations are represented by sets of cells (rounded rectangles) harboring available homing sites (grey boxes) or homing sites occupied by active (green boxes) or inactive intron copies (dark red boxes). The starred state in each model is those predicted to be the most frequently observed when a group II intron copy is present. Smaller cartoon sizes (in A) reflect infrequent alternative. doi:10.1371/journal.pone.0052268.g001

Figure 2 .

 2 Figure 2. Distribution of (A) 90 Enterobacteriales genomes (corresponding to 18 genera/complexes) according to their group II intron abundance, and (B) group II intron abundance according to sequenced strain number within genera/complexes. White dots: total number of detections considered; black dots: only distinct copies considered (see text). Names of genera/complexes with .20 intron copies and/or .5 sequenced strains are indicated. doi:10.1371/journal.pone.0052268.g002

Figure 3 .

 3 Figure 3. Inferred homing sites for Kl.pn.I5, Di.ze.I1, So.gl.I1, and c-Ha.de.I1 introns. (A) Sequence spanning from positions 225 to +10 relative to full-length Kl.pn.I5 copy insertion sites are displayed and compared to predicted EBS. Identical bases and those which fit the Watson-Crick (stars) or the Wobble (plus) base-pairing with EBS motifs were used to construct the consensus homing site. (B) Regions surrounding Di.ze.I1 and So.gl.I1 full-length copies are displayed and compared to predicted EBS1 and EBS3. Most represented nucleotides at a given position and which fit the Watson-Crick (stars) or the Wobble (plus) base-pairing with EBS motifs were used to construct the consensus homing site. Stem-loops belonging to Rho-independent transcription terminators inferred by ARNold [58] are colored in blue (stems) and red (loops). Underlined sequences are stem-loops not detected by ARNold. (C) Weblogo diagram of the 225 to +10 region surrounding the 26 full-length c-Ha.de.I1 copies belonging to subgroup 1. Bases which show relative conservation were used to construct the consensus homing site. doi:10.1371/journal.pone.0052268.g003

Table 1 .

 1 Distribution of the 198 detected group II intron copies among the 18 Enterobacteriales genera/complexes.

		Number of	Number of	Distinct	Full-length	Distinct
	Bacterial genus	strains	copies	copies	distinct copies	fragments
	Escherichia/Shigella	35	122	67	40*	27
	Hamiltonella	1	3 8	38	30	8
	Photorabdus	2	1 0	9	2	7
	Sodalis	1	8	8	3	5

Table 2 .

 2 Full length group II intron copies, number of potential homing sites available in their host genomes, and proportion of homing site occupation.

	Group II intron	Host genome	Copy number	Available homing sites	Proportion of occupation (%)
	E.c.I2	E. coli UMN026	1	4	20
	E.c.I3	E. coli ED1a	2	0	100
		E. coli 536	1	0	100
	E.c.I4	E. coli ED1a	2	1	67
		E. coli 55989	6	1	86
		E. coli E24377A	7	0	100
		E. coli IAI39	15	27	36
		E. coli SMS_3_5	1*	/	/
		E. coli UTI89	1	2	33
		E. coli UMN026	1	6	14
		S. flexneri Sf301

Table 3 .

 3 Group II intron fragments, number of potential homing sites available in their host genomes, and proportion of homing site occupation.

			Fragment	Available	Proportion of
	Group II intron	Host genome	number	homing sites*	occupation (%)
	E.c.I9	E. coli BL21	1	4	20
		E. coli BL21_DE3	1	4	20
		E. coli REL606	1	5	17
		E. coli IAI1	1	3	25
		E. coli SE11	1	6	14
		E. coli E24377A	1	3	25
		E. coli EDL933	1	2	33
		E. coli Sakai	1	2	33
		E. coli EC4115	1	2	33
		E. coli TW14359	1	2	33
	E.c.I3	E. coli E24377A	1	6	14
		S. sonnei Ss046	3	0	100
		S. boydii Sb227	3	1	75
		S. boydii SbCDC3083	3	1	75
		S. flexneri Sf301	3	0	100
		S. dysenteriae Sd197	3	0	100
	E.c.I4	E. coli BW2952	1	1	50
		E. coli DH10B	1	1	50
		E. coli MG1655	1	1	50
		E. coli W3110	1	4	50
		E. coli SMS_3_5	1	10	9
		E. coli UTI89	1	3	25
		E. coli UMN026	1	7	13
		S. sonnei Ss046	1	15	6
		S. boydiii SbCDC3083	1	57	2
		S. flexneri Sf301	1	43	2
		S. dysenteriae Sd197	1	25	4
	So.gl.I1	S. glossinidius 'morsitans'	2	76	3
	c-Ha.de.I1 (subgroup 1)	H. defensa 5AT	1	49	2

*: at the expected time of fragment acquisition (i.e. including homing sites currently occupied by full-length copies). doi:10.1371/journal.pone.0052268.t003

PLOS ONE | www.plosone.org 1 December 2012 | Volume 7 | Issue 12 | e52268

PLOS ONE | www.plosone.org

December 2012 | Volume 7 | Issue 12 | e52268

December 2012 | Volume | Issue 12 | e52268

Acknowledgments

We thank Cle ´ment Gilbert for comments on an earlier version of the manuscript.

This research was funded by a Young Investigator ATIP Award from the Centre National de la Recherche Scientifique (CNRS) and a European Research Council Starting Grant to RC. SL was supported by a postdoctoral fellowship from the CNRS. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.