<title>Background</title> <p>TEs are widely distributed in eukaryotes, representing 50% of the human genome [##REF##11237011##1##], 15% of the Drosophila genome, and up to 70% in <italic>Zea mays </italic>[##UREF##0##2##]. Because of their capacity of transposition they are able to invade the genome and promote insertional mutations and chromosomal rearrangements. Recurrent mobility allows them to persist in spite of their harmful effects in the host [##REF##16888345##3##]. Most of the proposed models in population dynamic studies [##UREF##1##4##, ####REF##2559652##5##, ##REF##1334899##6##, ##REF##1334900##7##, ##REF##15466430##8####15466430##8##] suggest that TEs are able to invade the genome if their transposition rate is enough to balance out opposing forces as excision and selection against deleterious insertions and chromosomal arrangements. Yet, these models, often too general, do not consider that each element behaves depending on both its own characteristics and the history of the population to which it belongs. This challenge to standard reasoning is most relevant in colonizing populations [##REF##11486506##9##]. Several authors have suggested that bursts of transposition could be induced in colonization by the foreign, often stressful, environment faced by the founders of colonizing populations [##REF##12082134##10##,##REF##9391219##11##]. Moreover, colonizing populations are subjected to well documented founder, drift effects [##UREF##2##12##]. Both processes generate population instabilities that may incorporate new variables to the interpretation of TE occupancy profiles in colonizing populations. These considerations qualify the study of TEs in colonization as of prime interest to understanding their invasive dynamics and putative evolutionary role in populations.</p> <p>Colonization effects on TEs were studied in Drosophila species [##REF##11486506##9##,##REF##9391219##11##,##REF##12572617##13##] showing that this process plays an important role in the TE chromosomal distribution. In particular studies in colonizing populations of <italic>D. buzzatii </italic>showed a TE bimodal distribution with sites either highly occupied, in a few cases, or showing low insertion occupancy, in most cases. Molecular studies of TE copies from high and low occupied sites [##REF##17151248##14##] strongly indicated that the most reliable explanation of the observed bimodal distribution is that a founder effect followed by genetic drift occurred during the colonization process. These results notwithstanding, valid for <italic>D. buzzatii</italic>, cannot be generalized to other colonizing Drosophila species, with different genomic characteristics, and subjected to different environmental pressures.</p> <p><italic>D. subobscura</italic>, a Paleartic species belonging to the obscura group [##UREF##3##15##] and characterized by a rich inversion polymorphism [##UREF##4##16##], has colonized North and South America almost 30 years ago [##UREF##5##17##,##REF##16593967##18##]. It was found for the first time in Puerto Montt (Chile) in 1978 [##UREF##6##19##] and later near Port Townsend in Washington (USA) in 1982 [##UREF##7##20##]. Thereafter this species showed a rapid spread and adaptation to the new colonized environment in form of latitudinal clines for chromosomal polymorphism and body size that paralleled the Paleartic clines [##UREF##5##17##,##REF##16593967##18##,##UREF##8##21##]. Main after-colonization population effects were the presence of allelic lethal genes in different populations [##UREF##9##22##], the low genetic variability of mtDNA [##REF##16578796##23##,##REF##2299979##24##] and the reduction of microsatellite allele numbers [##REF##11319257##25##] compared to original founder populations. These are expected outcomes of the founder drift effect of colonization. However nothing is known of the impact of colonization on the TE chromosomal distribution in this species.</p> <p>Here we present the study of the distribution of two TEs, <italic>gypsy </italic>and <italic>bilbo</italic>, in original and colonizing populations of <italic>D. subobscura</italic>. Results show that TE frequency distribution differ between original and colonizing populations in a way that colonization, chromosomal inversion polymorphism and particular characteristics inherent to each element can provide a sufficient likely explanation. In this paper we particularly emphasize the importance of population structure and history to explain TE distribution in natural populations.</p>
<title>Methods</title> <title>Drosophila strains</title> <p>The control strain <italic>chcu </italic>carries the recessive markers cherry eyes and curled wings and is homokaryotypic for chromosomal arrangements A<sub>st</sub>, J<sub>st</sub>, U<sub>st</sub>, E<sub>st </sub>and O<sub>3+4</sub>. It is kept by mass-culturing to maintain its viability. <italic>In situ </italic>hybridization for insertions of <italic>bilbo </italic>and <italic>gypsy </italic>displayed high stability over generations in 19C, 46A, 46C, 73A, 81D, 84A, 96A for <italic>bilbo </italic>and in 7C and 52A for <italic>gypsy</italic>.</p> <p>The original population was sampled in Spring 2005 in Bordils (42.30°N, Girona, Spain). The colonizing populations were sampled in Spring 2004 in Davis (38.33°N, California, USA) and Bellingham (48.45°N, Washington, USA), and in Spring 2005 in Maipú (36.52°S, Argentina).</p> <title>Mating system (prior to "<italic>in situ</italic>" hybridization)</title> <p>Individual males of natural populations were crossed with virgin females of the control line <italic>chcu</italic>. Insertion profiles were analyzed in F<sub>1 </sub>female larval progeny to include the X chromosome. The TE insertion profile of each male was deduced by subtracting the TE insertion profile of the control line from that of the F1 larva.</p> <title><italic>In situ </italic>hybridization and DNA probes</title> <p>Polytene chromosome [##UREF##4##16##] squashes from salivary glands of third-instar larvae, prepared as described in [##REF##2185305##79##], were hybridized with digoxigenin labelled probes of <italic>bilbo </italic>and <italic>gypsy</italic>. The probes consisted of PCR fragments (2.6 and 2.8 kb long) which included the reverse transcriptase region. Prehybridization solutions and posthybridization washes were done following a protocol by Roche [##UREF##23##80##]. PCR reactions were carried out in a final volume of 25 μl, including 1× activity buffer (Ecogen), 1.6 mM MgCl<sub>2</sub>, 0.2 mM of each dNTP (Roche), 0.4 μM primer (Roche), 10–20 ng of genomic template DNA, and 0.04 units per μl of Taq polymerase (Ecotaq from Ecogen). Amplifications were run in a MJ Research Inc. thermocycler programmed as follows: 5 min preliminary denaturation at 94°, 30 cycles of 45 s at 94° (denaturation), 45 s at specific PCR annealing temperatures, 1.5 min at 72° (extension) and a final extension for 10 min at 72°. PCR products were gel purified with a Geneclean kit (BIO 101) and labelled using the random primer method. After hybridization signal development was done using an anti-digoxigenin antibody conjugated with alkaline phosphatase (Roche).</p> <p><italic>In situ </italic>hybridization is the more suitable method used in localization of TEs on chromosomal arms. However, the power of resolution of this technique allow us neither discriminate between closely neighbouring sites, nor between elements that diverge below 10%.</p> <title>Statistical analyses</title> <p>Statistical analyses were performed excluding centromeric and pericentromeric TEs insertions. The statistical software SPSS version 14.0 was used for most of the statistical data analyses.</p> <p>In cases of multiple testing, corrections were achieved measuring the significance of False Discovery Rates [##UREF##24##81##] through q values. To get the q-value we used the software QVALUE [##UREF##25##82##] on the p values obtained from the multiple test. When this test could not be applied, Bonferroni's correction was performed [##UREF##26##83##].</p>
<title>Results</title> <title>Chromosomal distribution of <italic>bilbo </italic>and <italic>gypsy</italic></title> <p>We analyzed the distribution of <italic>bilbo </italic>and <italic>gypsy </italic>in polytene chromosomes of <italic>D. subobscura</italic>. Fig. ##FIG##0##1## shows two examples of chromosomal distribution: <italic>bilbo </italic>in chromosome O and <italic>gypsy </italic>in chromosome U. A different distribution pattern is observed, in general, when we compare colonizing and original populations. Colonizing populations present insertion frequencies of <italic>bilbo </italic>and <italic>gypsy </italic>higher than those of the original population. In general the same distribution pattern is observed for the rest of chromosomes. Ten sites (7A, 16A, 20A, 45C, 58D, 74D 82A, 83C, 85A, 89C) show a <italic>bilbo </italic>insertion polymorphism greater than 32% in at least one colonizing population. <italic>Gypsy </italic>insertion frequencies are lower than those of <italic>bilbo </italic>with an occupancy of more than 10% in eight chromosomal sites (39D, 41C, 43D, 49D, 52D, 63C, 71B, 74D). Differences in occupancy profiles between original and colonizing populations are represented in Table ##TAB##0##1## that shows the distribution of the number of times that each site is occupied in the studied sample. Thus, in <italic>bilbo </italic>the occupancy frequency ranges from 1 to 51 times in colonizing populations and only from 1 to 19 in the original population. Although <italic>gypsy </italic>shows a low occupancy profile compared to <italic>bilbo </italic>(colonizing populations range: 1–15; original population range: 1–5), the occupancy rate of both TEs in colonizing populations is greater than in the original population. The highest <italic>bilbo </italic>insertion frequencies are observed, in decreasing order, in Bellingham, Maipú and Davis. In the original population of Bordils, the highest insertion frequency corresponds to one site observed 16 times.</p> <p>Table ##TAB##1##2## lists the means and variances of copy number for <italic>bilbo </italic>and <italic>gypsy </italic>per chromosome and haploid genome. The mean copy number of both TEs for the whole genome (HG) is always higher in colonizer populations than in the original one. The <italic>bilbo </italic>mean copy number differs greatly among chromosomes ranging from 2.58 copies in chromosome O from Bellingham to 0.55 in chromosome J from Davis. In fact chromosome J hosts the lowest number of <italic>bilbo </italic>in all populations. A different scenario is found for <italic>gypsy </italic>in which the A (X) chromosome contains the lowest number of insertions, in colonizer and original populations alike. However, among the autosomes J is the least occupied in all populations. Deviation from Poisson distribution was tested by chi-square goodness of fit tests (for details, see additional files ##SUPPL##0##1## and ##SUPPL##1##2##) pooling adjacent classes with low expected numbers. In colonizing populations <italic>bilbo </italic>distribution in each chromosome fits a Poisson distribution. <italic>Gypsy </italic>deviates from a Poisson distribution in E chromosome from Davis and Bellingham and in U chromosome from Maipú. When the whole genome is considered both TEs follow a Poisson distribution in the original population and deviate in all colonizing ones, except for <italic>bilbo </italic>in Bellingham and gypsy in Maipú. For this element the general trend in colonizing populations is a lower than expected number of genomes with a single copy and an excess of genomes with three or more copies (see Table ##TAB##0##1##). An alternative test was performed using dispersion coefficients (DC), which measure the ratio between the variance (V<sub>n</sub>) and the mean (m) (DC = V<sub>n</sub>/m, see table ##TAB##1##2##). DC of 1 indicates that TE distribution is Poisson, and DC &gt; 1 or DC &lt; 1 indicates contagious or repulsive distributions, respectively. When the haploid genome is considered, there is a general tendency towards DCs &gt; 1 for both elements in all populations except for <italic>gypsy </italic>in Maipú (these results are due to the greater effect of some chromosomes in the final result of the test).</p> <p>Because in some cases TE sites seem to be distributed in a contagious way (DC &gt; 1), linkage disequilibrium was computed for each pair of sites by way of 2 × 2 contingency tables [##REF##2558961##26##]. Linkage disequilibrium between TE sites could be responsible of the non-random distribution detected in some cases. The observed distribution of correlation coefficients between all paired sites was compared to the expected distribution in absence of linkage disequilibrium using Fisher's hypergeometric formula [##UREF##10##27##]. Figure ##FIG##1##2## depicts, as an example, the correlation coefficient distributions (pooled in intervals of 0.1) of <italic>bilbo </italic>in chromome E from Bordils and in chromosome O from Maipú; and of <italic>gypsy </italic>in chromosome E from Davis and Bellingham. Tests were significant in most cases where a deviation from Poisson distribution was observed. Moreover we also found significant results in some cases where departures from Poisson distribution were not detected (e.g <italic>bilbo </italic>on chromosome O of Maipú). The general trend is a defect of class (-0.09–0.00) and an excess of some positive correlation classes. This indicates that some sites tend to stay together, as indicated by a DC &gt; 1. This tendency was observed in all cases where deviations from Poisson distribution were observed, except for <italic>gypsy </italic>in chromosome U from Maipú where there is an overabundance in class (-0.09–0.00) and the DC is lower than 1.</p> <title>Copy number comparisons among chromosomes</title> <p>Montgomery et al [##REF##3032743##28##] proposed that selection against TE insertions would lead to a lower number of TE copies in chromosome X than in autosomes due to the stronger deleterious effect of recessive insertional mutations in the X chromosome of hemizygous males. In order to test this hypothesis we compared the copy number in the A (X) chromosome with that in autosomes. To estimate the expected number of insertions we multiply the relative proportion of chromatin of each chromosome by the number of total insertions in the population. The relative proportion of chromatin is that reported by Stumm-Zollinger and Goldschmidt [##UREF##11##29##] corrected by eliminating the dot chromosome, not included in our analyses. If TEs are randomly distributed, we expect a TE copy number per chromosome proportional to the amount of chromatin.</p> <p>Observed and expected proportions were compared by a G test [##UREF##12##30##] among all chromosomes (G<sub>a</sub>), between the A (X) chromosome and autosomes (G<sub>b</sub>), and among autosomes (G<sub>c</sub>), as indicated in Table ##TAB##2##3##. G<sub>b </sub>values were significant for <italic>gypsy </italic>in all populations, and for <italic>bilbo </italic>only in Maipú and Bordils. Because some differences may be due to high insertion sites, additional analyses were done after eliminating these sites. After elimination the significance was maintained for <italic>gypsy </italic>in all populations except in Maipú, and removed for <italic>bilbo</italic>. In general <italic>gypsy </italic>shows a low copy number in the A (X) chromosome compared to autosomes. However, this is not the rule for <italic>bilbo </italic>where Maipú and Bordils show a high copy number in A (X). Interestingly, those populations that display <italic>gypsy </italic>copy number differences between A (X) and autosomes, show also significant differences among autosomes (G<sub>c</sub>), specially in colonizer populations where chromosomes E and O show a higher copy number than expected.</p> <p>In general, copy number tend to be higher for <italic>bilbo </italic>in chromosome O and for <italic>gypsy </italic>in chromosome U in all colonizing populations, (except <italic>bilbo </italic>in Maipú), whereas in the original population the E chromosome hosts the highest proportion of <italic>gypsy </italic>and <italic>bilbo</italic>. In order to determine if chromosomal differences have the same tendency in colonizing populations, heterogeneity (H) tests were performed for comparisons among chromosomes, between A (X) and autosomes and among autosomes. Table ##TAB##2##3## shows that all cases were heterogeneous for both TEs. However when Maipú is excluded from the analyses and high insertion frequency sites are eliminated, Bellingham and Davis become homogeneous for <italic>bilbo </italic>(data not shown).</p> <title>Correlation studies between high frequency sites and chromosomal arrangements</title> <p>All five pairs of acrocentric chromosomes of <italic>D. subobscura </italic>are polymorphic for inversions. Frequencies of chromosomal arrangements show clinal variation correlated with latitude in Paleartic populations [##REF##14503625##31##,##UREF##13##32##] and clines that follow the same latitudinal gradient evolved in recent colonizing populations in both hemispheres of the Americas [##UREF##5##17##,##REF##16593967##18##]. These parallel observations across continents provided a natural experiment that supports the adaptative role of the chromosomal inversion polymorphism.</p> <p>Frequencies of chromosomal arrangements in the analyzed populations of this work are summarized in Table ##TAB##3##4##. Each arrangement is conventionally designed by the letter of the chromosome in which it occurs, followed by a combination of digits that identify the set of inversions included in it [##UREF##4##16##]. Arrangement frequencies are of the same order of magnitude as those previously reported, including the North-South latitudinal variation of most arrangements [##UREF##8##21##,##REF##14503625##31##,##UREF##13##32##]. However it is interesting to note that Maipú presents a higher O<sub>St </sub>frequency than expected according to its latitude.</p> <p>Some authors consider recombination as the main factor determining the chromosomal distribution of TEs [##REF##12032249##33##,##REF##11875027##34##], but see [##REF##11875027##34##]. The model of ectopic exchange, predicts a negative correlation between recombination rate and TE copy number if ectopic exchange is reduced in parallel with regular meiotic recombination rate [##REF##8878672##35##,##REF##2854088##36##]. Under this model, TEs are expected to be more abundant in regions of low recombination as inversions or inversion break-points. In these regions the probability of induction of deleterious rearrangements produced by unequal recombination between TEs, is low because most of the time inversions will be found in heterozygous state (recombination is suppressed inside). Experimental evidences [##REF##10411506##37##, ####UREF##14##38##, ##REF##9770505##39####9770505##39##] suggest that TEs are responsible of chromosomal inversions in natural populations of Diptera and are particularly abundant inside and near inversion break-points [##REF##1334899##6##,##REF##1334900##7##,##REF##8088526##40##].</p> <p>In order to know whether an association between high insertion sites and arrangements exist, we computed the product-moment correlation coefficient (r) for high-frequency sites (Table ##TAB##4##5##). We observed two <italic>bilbo </italic>sites of particular interest (67A and 89C) that show the highest correlation coefficients. The 67A site is located inside the breakpoint of arrangement E<sub>12 </sub>and is significantly associated with E<sub>1+2+9+12 </sub>in Davis (r = 0.64) and Maipú (r = 0.85). The 89C site is located near the break-point of O<sub>8 </sub>arrangement and is significantly correlated with arrangement O<sub>3+4+8 </sub>in Davis (r = 0.58) and Bellingham (r = 0.34), and only marginally (r = 0.26) in Maipú. Other instances of significant associations are not so easily explained because sites are external to inversion breakpoints. Thus, highly occupied 74D <italic>bilbo </italic>site is located outside of chromosomal inversions, yet, it is also significantly associated with E<sub>1+2+9+12 </sub>in Davis (r = 0,33) and Maipú (r = 0,24). This site is also highly occupied by <italic>gypsy </italic>but in this case associations are not significant. In other cases we observe associations of sites inside highly frequent inversions where the crossing-over is not reduced. This is the case of 11B <italic>bilbo </italic>site, for example, negatively associated to A<sub>2 </sub>arrangement in all populations except Bellingham but located inside it.</p>
<title>Discussion</title> <title><italic>Bilbo </italic>and <italic>gypsy </italic>distributions are different in original and colonizer populations</title> <p>Results show a clear differential TE distribution in original and colonizing populations. While in the original population most sites have low insertion frequencies, colonizing populations present some highly occupied sites, with frequencies higher than 50% for <italic>bilbo </italic>and close to 20% for <italic>gypsy</italic>. Interestingly, most of them are common to all populations. Mean copy number of both elements is higher in colonizing populations than in the original one due to the presence of these highly occupied sites.</p> <p>Low occupied sites would represent insertions occurred after colonization and/or copies from the original population whose frequency is decreasing in colonizing populations. An argument in favour of the former hypothesis is the existence of unique sites that would correspond to new transpositions (i.e. site 48D of gypsy), while the latter hypothesis explains the existence of low-occupancy original sites common to different populations (i.e. site 41A of gypsy or 85B of bilbo).</p> <p>High insertion frequency sites are most likely due to a founder event during the colonization process (the founder hypothesis), as previously reported in other Drosophila species [##REF##11486506##9##,##REF##9391219##11##,##REF##12572617##13##]. In <italic>D. buzzatii </italic>this hypothesis was also verified by molecular studies showing identical <italic>Osvaldo </italic>retrotransposon structures and flanking genomic sequences in high insertion frequency sites from different colonizing populations [##REF##17151248##14##].</p> <p>In this study two lines of evidence support the founder hypothesis. First, the two studied TEs belong to different subclasses, yet they show a similar population behaviour. Second, most highly occupied sites are located in colonizing population chromosomes, although some exceptions occur for <italic>bilbo </italic>whose insertion frequency exceeds 10% in 9 sites in the original population of Bordils. All these sites correspond also to high insertion sites in colonizing populations, except 90C site and 21A, which are, respectively, free of insertions or occupied at low frequency in America.</p> <p>The presence of high frequency sites in the original population could be a consequence of the transposition mechanism of <italic>bilbo</italic>, a LINE element. It has been shown that LINE elements (L1) make 5' truncated copies during their transposition mechanism indicating that 5' sequences are not absolutely necessary to insertion [##REF##17023124##41##, ####REF##2550903##42##, ##REF##2445384##43####2445384##43##]. In fact, the majority of the L1 copies present in mammalian genomes are 5' truncated with a length of not more than 1 kb [##REF##11237011##1##,##REF##12368240##44##]. We can think that selection against truncated, "dead-on-arrival" (DOA) copies should be weak because they are not transcribed, potentially immobile and shorter than full copies. Thus, deleted copies could persist in some genomic regions without being completely eliminated by natural selection. In fact, some Drosophila TE families (most of them LINE like elements) seem to be only marginally affected by purifying selection, reaching high insertion frequencies in euchromatin [##REF##12716993##45##].</p> <p>On the other hand, some of <italic>bilbo </italic>high frequency sites from Bordils could be explained by the dragging effect from the rich inversion polymorphism of <italic>D. subobscura</italic>. For example the 67A site located in the break-point of E<sub>12 </sub>arrangement presents highly significant correlations with this arrangement in 2 out of 3 colonizing populations. In Bordils, this correlation is not significant because of the lower frequency of this arrangement in this population. Arrangements of chromosome E cover approximately 75% of its length and it is not rare to find this kind of associations. In this chromosome another high insertion frequency site (74D) shows association with the same E<sub>1+2+9+12 </sub>chromosomal arrangement. This site corresponds to a heterochromatic telomeric site where it is not rare to find an accumulation of TE insertions. In fact <italic>gypsy </italic>is inserted also in this chromosomal site at occupation rates that range from 1.3 to 11.4%. Accumulation of TEs in heterochromatin is well documented in <italic>D. melanogaster </italic>where a significant excess of insertions were reported in heterochromatin, dot and Y chromosomes alike [##REF##15219154##46##, ####REF##7664614##47##, ##REF##8078581##48##, ##REF##12654931##49####12654931##49##].</p> <p>Seasonal fluctuations in population frequency of chromosomal rearrangements can modify recombination rates and associations between arrangements and genes. In <italic>D. subobscura </italic>no seasonal fluctuations were reported in some works [##REF##5584166##50##,##UREF##15##51##], but fluctuations and seasonal changes of associations between chromosomal inversions and allozymes were reported in others, specially in the O chromosome from original populations [##REF##17246183##52##,##REF##12778553##53##]. In the present case, we observe no associations between insertions and specific chromosomal arrangements in the original population, but we do detect this kind of associations in colonizing populations (where fluctuations were not studied). However, changes in associations between chromosomal arrangements and chromosomal sites do not follow a definite trend. As an example, the U<sub>ST </sub>arrangement, whose frequency has increased in all colonizing populations, shows a positive association to 43B and 45C sites but a negative one to 53A site in Maipu. This is a rather odd outcome since increase of rearrangement frequency is always expected to break down associations due to an increase of recombination rate. So, the likely explanation would be that fluctuations do not affect associations or at least not in the same way for every studied rearrangement polymorphism.</p> <p>On the other hand, we favor the general idea that the positive correlation between arrangements and TE copies is not due to an inversion effect but, most probably, to the founder event [##UREF##6##19##,##REF##11319257##25##,##REF##7995925##54##,##UREF##16##55##]. This could explain why arrangement E<sub>1+2+9+12</sub>and the 74D site, which is located outside of the inversion, show a positive association and also why an excess of classes including positive correlation coefficients between chromosomal sites was observed in some chromosomes like E. Genetic estimates suggest that the number of founders ranged from 10 to 150 [##REF##11319257##25##,##UREF##17##56##]. If some founders carried together this site and this arrangement, both will appear together in all populations because they are identical by descent. The founder hypothesis is favored by the fact that all correlations between sites and arrangements are significant only in colonizing populations. In the original population in spite of having correlation coefficients of 0.57 (in 57D) and 0.56 (in 83C) with E<sub>1+2+9 </sub>and O<sub>3+4+7 </sub>respectively, these are not significant. In fact, these two arrangements are currently decreasing in frequency in the Mediterranean populations and perhaps these combinations descend also from a few individuals. All these considerations suggest that most of the associations detected are due to a founder effect.</p> <p>The general rule, as reported in <italic>D. melanogaster </italic>[##REF##12716993##45##,##REF##1334919##57##], is that TEs are spread and have low insertion frequencies in euchromatin. In some cases, however, accumulations of TEs in some chromosomal sites have been reported, as in the 42B [##REF##2548982##58##], 87C [##UREF##18##59##] and 38 [##UREF##19##60##] regions, of <italic>D. melanogaster</italic>, and the 85D region of <italic>D. subobscura </italic>[##REF##2612293##61##,##UREF##20##62##], and even fixation has occurred, as in the 42C site in natural populations of <italic>D. simulans </italic>[##REF##8642613##63##]. Preferential insertion sites (hotspots) have been suggested for some Drosophila elements [##UREF##21##64##, ####REF##6266759##65##, ##REF##2832692##66####2832692##66##] and we cannot completely discard the possibility of an activation of transposition to specific hotspots during the colonization process. This hypothesis could be verified if a process affecting equally the two TEs studied occurred, as shown in <italic>D. melanogaster</italic>. In this species some proteins are involved in RNA-silencing mechanisms for retrotransposable elements repression [##REF##17277359##67##, ####REF##17135323##68##, ##REF##16809489##69####16809489##69##]. We cannot discard the existence of a similar mechanism in <italic>D. subobscura </italic>that was de-repressed as a consequence of the colonization process contributing simultaneously to an increase of transposition of different transposable elements.</p> <title>Factors affecting TEs distribution in <italic>D. subobscura</italic></title> <p>In Drosophila, TEs seem to be maintained in populations as the result of a balance between transposition and opposing forces that reduce their copy number. In this way selection can act either directly against deleterious insertions or indirectly against deleterious chromosomal rearrangements produced by ectopic recombination between TEs [##UREF##1##4##,##REF##2559652##5##,##REF##2854088##36##,##REF##17246144##70##]. In this work a test of selection against deleterious insertions was done by comparing copy numbers between X and autosomes, selection being more effective in the former than in the latter.</p> <p>For <italic>gypsy </italic>we observe a clear tendency to follow a selection model, except in Maipú. This result is in concordance with that observed in a natural population of <italic>D. melanogaster </italic>with this element [##REF##8082837##71##]. For <italic>bilbo </italic>the data do not fit a selection model against deleterious insertions; even in those cases where the test is significant, a higher copy number on A (X), compared to autosomes, is observed. A possible explanation of this result is that <italic>bilbo </italic>could have a differential transposition rate between X and autosomes. Some examples of transposition restricted to female or male <italic>D. melanogaster </italic>germ line have been reported [##UREF##22##72##,##REF##9230904##73##] and they should be taken into account when X and autosomes are compared. On the other hand, the discrepancies observed between the two elements may be accounted for by the different factors that control copy numbers in each of these elements. In <italic>D. melanogaster gypsy </italic>is a retrovirus [##REF##8108403##74##] submitted probably to a strong selection effect, its transposition depending on the presence of permissive alleles most likely segregating in natural populations. In <italic>D. subobscura </italic>this retrotransposon seems to be non infectious because current available copies have an apparently inactive <italic>env </italic>region [##REF##8600460##75##], but this does not discard the putative presence of alleles that control its transposition. On the other hand <italic>bilbo </italic>is a LINE element and could be submitted to a soft selection pressure due to its DOA transposition mechanism. Most of the copies are probably deleted and its deleterious capability by transposition is diminished. The model of selection against deleterious insertions has been questioned by some authors [##REF##3032743##28##,##REF##8078581##48##] because neither all ETs nor all populations had a lower insertion frequency on X chromosomes compared to autosomes. However in a later work [##REF##9440269##76##], where the authors reanalyze the data including more results from other species, selection against insertions is considered as the major mechanism of TE copy number control. On the other hand, values of selection coefficients against deleterious mutations could not be comparable to mutations associated to TE insertions. Moreover, deleterious effects of TEs can be species specific and populations may also sometimes suffer TE mobilizations that mask selection effects on TE distribution.</p> <p>In this work each element presents a different behavior probably due to their distinct transposition mechanisms. Moreover we should not forget that elements which are stable in some genome conditions could be unstable in others. Recently mobilized TEs and/or colonization events, in populations, could lead to a differential copy distribution between chromosomes, rendering the selection undetected. This could be the case of Maipú, a new colonizing Argentinian population, which shows a distribution pattern for <italic>gypsy </italic>and <italic>bilbo </italic>quite different from the other colonizing populations. In particular, some high insertion frequency sites are more represented, or even exclusive, in this population. It is possible that Maipú was established through a bottleneck of founder flies from Chile as a consequence of a secondary colonization. In this case, we cannot discard the existence of new transposition events in founders induced by the new environmental conditions encountered as previously proposed by other authors [##REF##12082134##10##,##REF##9391219##11##]. If this colonization occurred recently, as indicated by collecting records, selection has not had enough time to act, explaining the discrepancies in this population when comparing A (X) and autosome copy numbers in Table ##TAB##2##3## or when this population is included in heterogeneity tests. In addition if TEs are not at equilibrium, departures from random distribution across chromosomes could reflect the insertion pattern rather than the effect of natural selection.</p> <p>Another model proposed to explain the TE dynamics is the selection against deleterious arrangements produced by ectopic recombination between TEs. In <italic>D. subobscura </italic>accurate measures of recombination rate are not available and it is not possible to calculate a correlation between TE copies and recombination rates. This species has a rich inversion polymorphism in all chromosomes and recombination is reduced in heterokaryotypes. Under this model we expect accumulation of TEs in inverted segments, and in inversion break points or near them. In some cases arrangements include overlapping inverted fragments, often reaching frequencies higher than the standard arrangements, but in other cases, of low frequency arrangements, TE copy number is too low to allow statistical tests. Also recombination between non-overlapping inversions or inversion complexes may also be prevented [##REF##9652229##77##].</p> <p>We looked for accumulations of <italic>bilbo </italic>and <italic>gypsy </italic>in breakpoints of inversions but only one high insertion frequency <italic>bilbo </italic>site, 67A, coincides with an inversion breakpoint (E<sub>1+2+9+12</sub>). In another case the 89C high frequency site of <italic>bilbo </italic>is located near the inversion O<sub>8 </sub>and shows a significant correlation with O<sub>3+4+8 </sub>arrangement. This is in concordance with several unsuccessful attempts to localize <italic>in situ </italic>hybridization middle repeated sequences in <italic>D. subobscura </italic>inversions breakpoints [##REF##2612293##61##,##REF##7705632##78##]. These data notwithstanding, we cannot discard that other elements may be responsible of chromosomal inversion induction as reported in other Drosophila species [##REF##10411506##37##,##UREF##14##38##].</p>
<title>Conclusion</title> <p>We conclude that the differential distribution of <italic>bilbo </italic>and <italic>gypsy </italic>between original and colonizing <italic>D. subobscura </italic>populations, is mainly due to a founder effect occurred during the colonization process of this species. We have shown that both founder effect and inversion polymorphism contribute notably to an excess of positive correlations between site pairs. Moreover the two transposable elements show a different pattern of distribution in populations that might be due to their differences in transposition and copy number regulatory mechanisms. This paper is also an attempt to emphasize the importance of population structure and history to explain the TE chromosomal distribution. We highlight the fact that comparisons in TE copy number between X and autosomes have to be interpreted cautiously. Sometimes TEs mobilizations can mask the effect of selection on TE distribution.</p>
<title>Background</title> <p>Transposable elements (TEs) constitute a substantial amount of all eukaryotic genomes. They induce an important proportion of deleterious mutations by insertion into genes or gene regulatory regions. However, their mutational capabilities are not always adverse but can contribute to the genetic diversity and evolution of organisms. Knowledge of their distribution and activity in the genomes of populations under different environmental and demographic regimes, is important to understand their role in species evolution. In this work we study the chromosomal distribution of two TEs, <italic>gypsy </italic>and <italic>bilbo</italic>, in original and colonizing populations of <italic>Drosophila subobscura </italic>to reveal the putative effect of colonization on their insertion profile.</p> <title>Results</title> <p>Chromosomal frequency distribution of two TEs in one original and three colonizing populations of <italic>D. subobscura</italic>, is different. Whereas the original population shows a low insertion frequency in most TE sites, colonizing populations have a mixture of high (frequency ≥ 10%) and low insertion sites for both TEs. Most highly occupied sites are coincident among colonizing populations and some of them are correlated to chromosomal arrangements. Comparisons of TE copy number between the X chromosome and autosomes show that <italic>gypsy </italic>occupancy seems to be controlled by negative selection, but <italic>bilbo </italic>one does not.</p> <title>Conclusion</title> <p>These results are in accordance that TEs in <italic>Drosophila subobscura </italic>colonizing populations are submitted to a founder effect followed by genetic drift as a consequence of colonization. This would explain the high insertion frequencies of <italic>bilbo </italic>and <italic>gypsy </italic>in coincident sites of colonizing populations. High occupancy sites would represent insertion events prior to colonization. Sites of low frequency would be insertions that occurred after colonization and/or copies from the original population whose frequency is decreasing in colonizing populations. This work is a pioneer attempt to explain the chromosomal distribution of TEs in a colonizing species with high inversion polymorphism to reveal the putative effect of arrangements in TE insertion profiles. In general no associations between arrangements and TE have been found, except in a few cases where the association is very strong. Alternatively, founder drift effects, seem to play a leading role in TE genome distribution in colonizing populations.</p>
<title>Authors' contributions</title> <p>MPGG participated in the design, the chromosomal slides, some statistical analyses and the writing of the manuscript. BECS collected the Bordils population, performed most of the technical work, the reading of slides and some statistical analyses. JB collected Davis and Bellingham populations, supervised all arrangement readings and performed the data set analyses. LS collected the Bordils population, contributed to the design and thoroughly revised the manuscript AF participated in the design, directed the project, coordinated the data analyses, contributed to the writing of the manuscript and collected the Maipú and Bordils populations. All authors read and approved the final manuscript</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>This work was funded by projects BOS2003-05904-C02-01 and CGL2006-13423-C02-01/02 from the Ministerio de Educación y Ciencia, Spain, and project 2005SGR 00995 from Generalitat de Catalunya, Spain to the Grup de Biología Evolutiva (Principal Investigator: AF); and a doctoral fellowship from UAB to BECS. We acknowledge technical help by Montse Peiró and field assistance by Dr. E. Hasson, Dr. P. Fernández-Iriarte, Dr. F. García-Franco and D. Vela.</p>
<fig id="F1" position="float"><label>Figure 1</label><caption><p><bold>Distribution of bilbo and gypsy in chromosomes O and U, respectively, from colonizing (DA: Davis, BE: Bellingham, MA: Maipú) and original populations (BO: Bordils) of D. subobscura</bold>. Number of haploid genomes analyzed are given in parenthesis.</p></caption></fig> <fig id="F2" position="float"><label>Figure 2</label><caption><p><bold>Observed and expected frequency distributions of correlation coefficients between all pairs of sites in natural populations: A) </bold><italic>bilbo </italic>in chromosome E from Bordils and O from Maipú. <bold>B)</bold><italic>gypsy </italic>in chromosome E from Davis and Bellingham.</p></caption></fig>
<table-wrap id="T1" position="float"><label>Table 1</label><caption><p>Occupancy profiles of euchromatic sites in original and colonizing populations</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="left">TE</th><th align="center">Populations</th><th align="center" colspan="18">Occupancy profiles</th></tr></thead><tbody><tr><td/><td/><td/><td align="center"><bold>1</bold></td><td align="center"><bold>2</bold></td><td align="center"><bold>3</bold></td><td align="center"><bold>4</bold></td><td align="center"><bold>5</bold></td><td align="center"><bold>6</bold></td><td align="center"><bold>7</bold></td><td align="center"><bold>8</bold></td><td align="center"><bold>9</bold></td><td align="center"><bold>10</bold></td><td align="center"><bold>11</bold></td><td align="center"><bold>12</bold></td><td align="center"><bold>13</bold></td><td align="center"><bold>14</bold></td><td align="center"><bold>15</bold></td><td align="center"><bold>16</bold></td><td align="center"><bold>17–51</bold></td></tr><tr><td colspan="20"><hr/></td></tr><tr><td align="left"><bold>bilbo</bold></td><td align="center"><bold>Colonizing</bold></td><td align="left">DA</td><td align="center">7</td><td align="center">6</td><td align="center">9</td><td align="center">4</td><td align="center">1</td><td align="center">7</td><td align="center">0</td><td align="center">2</td><td align="center">2</td><td align="center">1</td><td align="center">5</td><td align="center">1</td><td align="center">1</td><td align="center">0</td><td align="center">1</td><td align="center">0</td><td align="center">7</td></tr><tr><td/><td/><td align="left">BE</td><td align="center">3</td><td align="center">5</td><td align="center">3</td><td align="center">2</td><td align="center">4</td><td align="center">2</td><td align="center">1</td><td align="center">2</td><td align="center">1</td><td align="center">2</td><td align="center">0</td><td align="center">0</td><td align="center">4</td><td align="center">1</td><td align="center">2</td><td align="center">2</td><td align="center">10</td></tr><tr><td/><td/><td align="left">MA</td><td align="center">8</td><td align="center">8</td><td align="center">4</td><td align="center">7</td><td align="center">4</td><td align="center">2</td><td align="center">2</td><td align="center">5</td><td align="center">2</td><td align="center">3</td><td align="center">2</td><td align="center">3</td><td align="center">2</td><td align="center">1</td><td align="center">2</td><td align="center">1</td><td align="center">12</td></tr><tr><td/><td align="center"><bold>Original</bold></td><td align="left">BO</td><td align="center">32</td><td align="center">19</td><td align="center">21</td><td align="center">11</td><td align="center">6</td><td align="center">2</td><td align="center">5</td><td align="center">1</td><td align="center">2</td><td align="center">0</td><td align="center">1</td><td align="center">0</td><td align="center">1</td><td align="center">0</td><td align="center">3</td><td align="center">1</td><td align="center">1</td></tr><tr><td colspan="20"><hr/></td></tr><tr><td align="left"><bold>gypsy</bold></td><td align="center"><bold>Colonizing</bold></td><td align="left">DA</td><td align="center">8</td><td align="center">4</td><td align="center">4</td><td align="center">0</td><td align="center">1</td><td align="center">2</td><td align="center">2</td><td align="center">2</td><td align="center">0</td><td align="center">0</td><td align="center">1</td><td align="center">0</td><td align="center">1</td><td align="center">1</td><td align="center">0</td><td/><td/></tr><tr><td/><td/><td align="left">BE</td><td align="center">11</td><td align="center">5</td><td align="center">1</td><td align="center">2</td><td align="center">2</td><td align="center">2</td><td align="center">1</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">1</td><td/><td/></tr><tr><td/><td/><td align="left">MA</td><td align="center">9</td><td align="center">10</td><td align="center">6</td><td align="center">3</td><td align="center">2</td><td align="center">1</td><td align="center">1</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">2</td><td align="center">0</td><td align="center">0</td><td/><td/></tr><tr><td/><td align="center"><bold>Original</bold></td><td align="left">BO</td><td align="center">10</td><td align="center">12</td><td align="center">5</td><td align="center">2</td><td align="center">1</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td/><td/></tr></tbody></table></table-wrap> <table-wrap id="T2" position="float"><label>Table 2</label><caption><p>Tests of the Poisson distribution of <italic>bilbo </italic>and <italic>gypsy </italic>per chromosome and haploid genome.</p></caption><table frame="hsides" rules="groups"><thead><tr><th/><th/><th align="center" colspan="10">Populations</th></tr><tr><th/><th/><th colspan="10"><hr/></th></tr><tr><th/><th/><th align="center" colspan="5">DA (76)</th><th align="center" colspan="5">BE (88)</th></tr><tr><th/><th/><th colspan="10"><hr/></th></tr><tr><th align="left">TE</th><th align="center">Ch.</th><th align="left">m</th><th align="center">V<sub>n</sub></th><th align="center">DC</th><th align="center"><bold>χ</bold><sup>2</sup></th><th align="center">df</th><th align="left">m</th><th align="center">V<sub>n</sub></th><th align="center">DC</th><th align="left"><bold>χ</bold><sup>2</sup></th><th align="left">df</th></tr></thead><tbody><tr><td/><td align="center"><bold>A</bold></td><td align="left">0.84</td><td align="center">0.72</td><td align="center">0.86</td><td align="center">0.61</td><td align="center">2</td><td align="left">1.07</td><td align="center">0.77</td><td align="center">0.73</td><td align="left">2.68</td><td align="left">2</td></tr><tr><td/><td align="center"><bold>J</bold></td><td align="left">0.55</td><td align="center">0.49</td><td align="center">0.89</td><td align="center">0.11</td><td align="center">1</td><td align="left">0.76</td><td align="center">0.67</td><td align="center">0.88</td><td align="left">1.61</td><td align="left">2</td></tr><tr><td/><td align="center"><bold>U</bold></td><td align="left">1.01</td><td align="center">0.87</td><td align="center">0.86</td><td align="center">5.74</td><td align="center">3</td><td align="left">1.11</td><td align="center">1.00</td><td align="center">0.90</td><td align="left">1.46</td><td align="left">3</td></tr><tr><td/><td align="center"><bold>E</bold></td><td align="left">1.38</td><td align="center">1.73</td><td align="center">1.25</td><td align="center">4.70</td><td align="center">4</td><td align="left">0.82</td><td align="center">0.84</td><td align="center">1.03</td><td align="left">0.69</td><td align="left">3</td></tr><tr><td/><td align="center"><bold>O</bold></td><td align="left">1.71</td><td align="center">1.67</td><td align="center">0.98</td><td align="center">3.47</td><td align="center">5</td><td align="left">2.58</td><td align="center">1.80</td><td align="center">0.70</td><td align="left">7.06</td><td align="left">4</td></tr><tr><td align="left"><bold>bilbo</bold></td><td align="center"><bold>HG</bold></td><td align="left">5.50</td><td align="center">9.40</td><td align="center">1.71**</td><td align="center">32.56**</td><td align="center">11</td><td align="left">6.34</td><td align="center">7.79</td><td align="center">1.23</td><td align="left">18.61</td><td align="left">11</td></tr><tr><td/><td/><td colspan="10"><hr/></td></tr><tr><td/><td/><td align="center" colspan="5"><bold>MA (81)</bold></td><td align="center" colspan="5"><bold>BO (81)</bold></td></tr><tr><td/><td/><td colspan="10"><hr/></td></tr><tr><td/><td/><td align="left"><bold>m</bold></td><td align="center"><bold>V<sub>n</sub></bold></td><td align="center"><bold>DC</bold></td><td align="center"><bold>χ</bold><sup>2</sup></td><td align="center"><bold>df</bold></td><td align="left"><bold>m</bold></td><td align="center"><bold>V<sub>n</sub></bold></td><td align="center"><bold>DC</bold></td><td align="left"><bold>χ</bold><sup>2</sup></td><td align="left"><bold>df</bold></td></tr><tr><td/><td/><td colspan="10"><hr/></td></tr><tr><td/><td align="center"><bold>A</bold></td><td align="left">1.75</td><td align="center">1.69</td><td align="center">0.96</td><td align="center">1.54</td><td align="center">5</td><td align="left">1.06</td><td align="center">0.71</td><td align="center">0.67</td><td align="left">8.11</td><td align="left">2</td></tr><tr><td/><td align="center"><bold>J</bold></td><td align="left">0.91</td><td align="center">0.93</td><td align="center">1.02</td><td align="center">0.25</td><td align="center">3</td><td align="left">0.75</td><td align="center">0.81</td><td align="center">1.08</td><td align="left">13.65*</td><td align="left">4</td></tr><tr><td/><td align="center"><bold>U</bold></td><td align="left">1.54</td><td align="center">1.28</td><td align="center">0.83</td><td align="center">4.15</td><td align="center">4</td><td align="left">0.84</td><td align="center">0.69</td><td align="center">0.82</td><td align="left">1.14</td><td align="left">2</td></tr><tr><td/><td align="center"><bold>E</bold></td><td align="left">1.80</td><td align="center">1.66</td><td align="center">0.92</td><td align="center">1.70</td><td align="center">4</td><td align="left">1.09</td><td align="center">1.70</td><td align="center">1.57**</td><td align="left">29.79**</td><td align="left">5</td></tr><tr><td/><td align="center"><bold>O</bold></td><td align="left">1.79</td><td align="center">1.99</td><td align="center">1.11</td><td align="center">6.53</td><td align="center">5</td><td align="left">0.96</td><td align="center">1.28</td><td align="center">1.34</td><td align="left">9.58</td><td align="left">3</td></tr><tr><td/><td align="center"><bold>HG</bold></td><td align="left">7.80</td><td align="center">19.51</td><td align="center">2.50**</td><td align="center">254.34**</td><td align="center">18</td><td align="left">4.70</td><td align="center">5.76</td><td align="center">1.22</td><td align="left">11.88</td><td align="left">10</td></tr><tr><td colspan="12"><hr/></td></tr><tr><td/><td/><td align="center" colspan="5"><bold>DA (70)</bold></td><td align="center" colspan="5"><bold>BE (84)</bold></td></tr><tr><td/><td/><td colspan="10"><hr/></td></tr><tr><td/><td/><td align="left"><bold>m</bold></td><td align="center"><bold>V<sub>n</sub></bold></td><td align="center"><bold>DC</bold></td><td align="center"><bold>χ</bold><sup>2</sup></td><td align="center"><bold>df</bold></td><td align="left"><bold>m</bold></td><td align="center"><bold>V<sub>n</sub></bold></td><td align="center"><bold>DC</bold></td><td align="left"><bold>χ</bold><sup>2</sup></td><td align="left"><bold>df</bold></td></tr><tr><td/><td/><td colspan="10"><hr/></td></tr><tr><td/><td align="center"><bold>A</bold></td><td align="left">0</td><td align="center">0</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="left">0.02</td><td align="center">0.02</td><td align="center">1.00</td><td align="left">-</td><td align="left">-</td></tr><tr><td/><td align="center"><bold>J</bold></td><td align="left">0</td><td align="center">0</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="left">0.02</td><td align="center">0.02</td><td align="center">1.00</td><td align="left">-</td><td align="left">-</td></tr><tr><td/><td align="center"><bold>U</bold></td><td align="left">0.98</td><td align="center">1.26</td><td align="center">1.28</td><td align="center">5.53</td><td align="center">3</td><td align="left">0.43</td><td align="center">0.49</td><td align="center">1.14</td><td align="left">1.76</td><td align="left">3</td></tr><tr><td/><td align="center"><bold>E</bold></td><td align="left">0.58</td><td align="center">0.97</td><td align="center">1.66**</td><td align="center">44.10**</td><td align="center">4</td><td align="left">0.33</td><td align="center">0.56</td><td align="center">1.69**</td><td align="left">38.10</td><td align="left">3</td></tr><tr><td/><td align="center"><bold>O</bold></td><td align="left">0.04</td><td align="center">0.04</td><td align="center">0.97</td><td align="center">-</td><td align="center">-</td><td align="left">0.09</td><td align="center">0.11</td><td align="center">1.17</td><td align="left">1.38</td><td align="left">1</td></tr><tr><td align="left"><bold>gypsy</bold></td><td align="center"><bold>HG</bold></td><td align="left">1.61</td><td align="center">2.44</td><td align="center">1.51*</td><td align="center">14.89</td><td align="center">5</td><td align="left">0.90</td><td align="center">1.60</td><td align="center">1.77**</td><td align="left">33.18**</td><td align="left">4</td></tr><tr><td/><td/><td colspan="10"><hr/></td></tr><tr><td/><td/><td align="center" colspan="5"><bold>MA (80)</bold></td><td align="center" colspan="5"><bold>BO (80)</bold></td></tr><tr><td/><td/><td colspan="10"><hr/></td></tr><tr><td/><td/><td align="left"><bold>m</bold></td><td align="center"><bold>V<sub>n</sub></bold></td><td align="center"><bold>DC</bold></td><td align="center"><bold>χ</bold><sup>2</sup></td><td align="center"><bold>df</bold></td><td align="left"><bold>m</bold></td><td align="center"><bold>V<sub>n</sub></bold></td><td align="center"><bold>DC</bold></td><td align="left"><bold>χ</bold><sup>2</sup></td><td align="left"><bold>df</bold></td></tr><tr><td/><td/><td colspan="10"><hr/></td></tr><tr><td/><td align="center"><bold>A</bold></td><td align="left">0.12</td><td align="center">0.16</td><td align="center">1.29</td><td align="center">4.46</td><td align="center">1</td><td align="left">0.02</td><td align="center">0.02</td><td align="center">1.00</td><td align="left">-</td><td align="left">-</td></tr><tr><td/><td align="center"><bold>J</bold></td><td align="left">0.14</td><td align="center">0.12</td><td align="center">0.87</td><td align="center">-</td><td align="center">-</td><td align="left">0.10</td><td align="center">0.12</td><td align="center">1.16</td><td align="left">1.26</td><td align="left">1</td></tr><tr><td/><td align="center"><bold>U</bold></td><td align="left">0.42</td><td align="center">0.27</td><td align="center">0.64</td><td align="center">8.49*</td><td align="center">1</td><td align="left">0.11</td><td align="center">0.10</td><td align="center">0.90</td><td align="left">-</td><td align="left">-</td></tr><tr><td/><td align="center"><bold>E</bold></td><td align="left">0.39</td><td align="center">0.32</td><td align="center">0.82</td><td align="center">1.43</td><td align="center">1</td><td align="left">0.30</td><td align="center">0.34</td><td align="center">1.13</td><td align="left">2.27</td><td align="left">1</td></tr><tr><td/><td align="center"><bold>O</bold></td><td align="left">0.27</td><td align="center">0.28</td><td align="center">1.01</td><td align="center">0.12</td><td align="center">1</td><td align="left">0.24</td><td align="center">0.23</td><td align="center">0.99</td><td align="left">0.00</td><td align="left">1</td></tr><tr><td/><td align="center"><bold>HG</bold></td><td align="left">1.35</td><td align="center">1.27</td><td align="center">0.94</td><td align="center">11.48</td><td align="center">4</td><td align="left">0.77</td><td align="center">0.88</td><td align="center">1.14</td><td align="left">3.44</td><td align="left">2</td></tr></tbody></table></table-wrap> <table-wrap id="T3" position="float"><label>Table 3</label><caption><p>Comparison of the proportion of <italic>gypsy </italic>and <italic>bilbo </italic>sites among chromosomes, autosomes and between chromosome A and autosomes</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="center">TE</th><th/><th align="center" colspan="4">gypsy</th><th align="center" colspan="4">bilbo</th></tr><tr><th colspan="6"><hr/></th><th colspan="4"><hr/></th></tr><tr><th align="center">Ch.</th><th align="center">P. chromat</th><th align="left">DA</th><th align="left">BE</th><th align="left">MA</th><th align="left">BO</th><th align="left">DA</th><th align="left">BE</th><th align="left">MA</th><th align="left">BO</th></tr></thead><tbody><tr><td align="center"><bold>A(X)</bold></td><td align="center">0.16</td><td align="left">0.00</td><td align="left">0.03</td><td align="left">0.09</td><td align="left">0.03</td><td align="left">0.15</td><td align="left">0.17</td><td align="left">0.22</td><td align="left">0.23</td></tr><tr><td align="center"><bold>J</bold></td><td align="center">0.20</td><td align="left">0.01</td><td align="left">0.03</td><td align="left">0.10</td><td align="left">0.13</td><td align="left">0.10</td><td align="left">0.12</td><td align="left">0.12</td><td align="left">0.16</td></tr><tr><td align="center"><bold>U</bold></td><td align="center">0.19</td><td align="left">0.61</td><td align="left">0.48</td><td align="left">0.31</td><td align="left">0.14</td><td align="left">0.18</td><td align="left">0.17</td><td align="left">0.20</td><td align="left">0.18</td></tr><tr><td align="center"><bold>E</bold></td><td align="center">0.20</td><td align="left">0.36</td><td align="left">0.37</td><td align="left">0.29</td><td align="left">0.39</td><td align="left">0.25</td><td align="left">0.13</td><td align="left">0.23</td><td align="left">0.23</td></tr><tr><td align="center"><bold>O</bold></td><td align="center">0.25</td><td align="left">0.03</td><td align="left">0.10</td><td align="left">0.20</td><td align="left">0.31</td><td align="left">0.31</td><td align="left">0.40</td><td align="left">0.23</td><td align="left">0.20</td></tr><tr><td colspan="10"><hr/></td></tr><tr><td/><td align="center"><bold>Df</bold></td><td/><td/><td/><td/><td/><td/><td/><td/></tr><tr><td align="center"><bold>Ga</bold></td><td align="center">4</td><td align="left">197.23** (115.82**)</td><td align="left">71.18** (61.31**)</td><td align="left">22.26** (5.08)</td><td align="left">21.40** (21.40**)</td><td align="left">36.98** (0.08)</td><td align="left">87.16** (16.98**)</td><td align="left">41.82** (4.80)</td><td align="left">15.67 (8.66)</td></tr><tr><td align="center"><bold>Gb</bold></td><td align="center">1</td><td align="left">40.48** (26.15**)</td><td align="left">15.24** (17.20**)</td><td align="left">4.63* (1.14)</td><td align="left">11.06**</td><td align="left">0.36 (4.8)</td><td align="left">0.08 (3.47)</td><td align="left">15.56** (2.9)</td><td align="left">9.69** (0.52)</td></tr><tr><td align="center"><bold>Gc</bold></td><td align="center">3</td><td align="left">156.75** (89.66**)</td><td align="left">55.93** (44.11**)</td><td align="left">17.63** 3.94</td><td align="left">10.34*</td><td align="left">36.61** (3.47)</td><td align="left">87.08** (13.50*)</td><td align="left">26.26** (1.87)</td><td align="left">5.97 (8.10*)</td></tr><tr><td colspan="10"><hr/></td></tr><tr><td align="center"><bold>Ga</bold></td><td/><td/><td/><td/><td/><td/><td/><td/><td/></tr><tr><td align="center">Total</td><td align="center">12</td><td align="center" colspan="4">290.67**(182.21**)</td><td align="center" colspan="4">165.96**(30.06**)</td></tr><tr><td align="center">Pooled</td><td align="center">4</td><td align="center" colspan="4">230.15*(110.30**)</td><td align="center" colspan="4">102.01**(1.8)</td></tr><tr><td align="center">H</td><td align="center">8</td><td align="center" colspan="4">60.52**(71.90**)</td><td align="center" colspan="4">63.95**(28.23**)</td></tr><tr><td align="center"><bold>Gb</bold></td><td/><td/><td/><td/><td/><td/><td/><td/><td/></tr><tr><td align="center">Total</td><td align="center">3</td><td align="center" colspan="4">60.36**(44.50**)</td><td align="center" colspan="4">16.00**(11.20**)</td></tr><tr><td align="center">Pooled</td><td align="center">1</td><td align="center" colspan="4">44.97**(23.78**)</td><td align="center" colspan="4">5.77*(1.09)</td></tr><tr><td align="center">H</td><td align="center">2</td><td align="center" colspan="4">15.38**(20.71**)</td><td align="center" colspan="4">10.24*(10.10*)</td></tr><tr><td align="center"><bold>Gc</bold></td><td/><td/><td/><td/><td/><td/><td/><td/><td/></tr><tr><td align="center">Total</td><td align="center">9</td><td align="center" colspan="4">230.32**(137.71**)</td><td align="center" colspan="4">149.95**(18.85)</td></tr><tr><td align="center">Pooled</td><td align="center">3</td><td align="center" colspan="4">185.18**(86.52**)</td><td align="center" colspan="4">96.25**(0.72)</td></tr><tr><td align="center">H</td><td align="center">6</td><td align="center" colspan="4">45.13**(51.19**)</td><td align="center" colspan="4">53.71**(18.12*)</td></tr></tbody></table></table-wrap> <table-wrap id="T4" position="float"><label>Table 4</label><caption><p>Frequencies of chromosomal arrangements in natural populations</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="left">Arrangements</th><th align="center" colspan="4">Populations</th></tr><tr><th colspan="1"><hr/></th><th colspan="4"><hr/></th></tr><tr><th/><th align="center">DA</th><th align="center">BE</th><th align="center">MA</th><th align="center">BO</th></tr></thead><tbody><tr><td align="left"><bold>Ast</bold></td><td align="center">0.57</td><td align="center">0.67</td><td align="center">0.55</td><td align="center">0.53</td></tr><tr><td align="left"><bold>A1</bold></td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">0.17</td></tr><tr><td align="left"><bold>A2</bold></td><td align="center">0.43</td><td align="center">0.33</td><td align="center">0.45</td><td align="center">0.30</td></tr><tr><td colspan="5"><hr/></td></tr><tr><td align="left"><bold>Jst</bold></td><td align="center">0.41</td><td align="center">0.45</td><td align="center">0.26</td><td align="center">0.40</td></tr><tr><td align="left"><bold>J1</bold></td><td align="center">0.59</td><td align="center">0.55</td><td align="center">0.74</td><td align="center">0.60</td></tr><tr><td colspan="5"><hr/></td></tr><tr><td align="left"><bold>Ust</bold></td><td align="center">0.32</td><td align="center">0.42</td><td align="center">0.50</td><td align="center">0.17</td></tr><tr><td align="left"><bold>U1</bold></td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">0.03</td></tr><tr><td align="left"><bold>U1+2</bold></td><td align="center">0.39</td><td align="center">0.40</td><td align="center">0.20</td><td align="center">0.72</td></tr><tr><td align="left"><bold>U1+2+8</bold></td><td align="center">0.29</td><td align="center">0.18</td><td align="center">0.30</td><td align="center">0.07</td></tr><tr><td align="left"><bold>U1+2+3</bold></td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">0.01</td></tr><tr><td colspan="5"><hr/></td></tr><tr><td align="left"><bold>Est</bold></td><td align="center">0.54</td><td align="center">0.74</td><td align="center">0.59</td><td align="center">0.53</td></tr><tr><td align="left"><bold>E8</bold></td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">0.01</td></tr><tr><td align="left"><bold>E1+2</bold></td><td align="center">0.07</td><td align="center">0.05</td><td align="center">-</td><td align="center">0.30</td></tr><tr><td align="left"><bold>E1+2+9</bold></td><td align="center">0.12</td><td align="center">0.01</td><td align="center">0.08</td><td align="center">0.03</td></tr><tr><td align="left"><bold>E1+2+9+12</bold></td><td align="center">0.05</td><td align="center">0.15</td><td align="center">0.23</td><td align="center">0.12</td></tr><tr><td align="left"><bold>E1+2+9+3</bold></td><td align="center">0.22</td><td align="center">0.05</td><td align="center">0.10</td><td align="center">0.01</td></tr><tr><td colspan="5"><hr/></td></tr><tr><td align="left"><bold>Ost</bold></td><td align="center">0.08</td><td align="center">0.23</td><td align="center">0.22</td><td align="center">0.22</td></tr><tr><td align="left"><bold>O2</bold></td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">0.01</td></tr><tr><td align="left"><bold>O5</bold></td><td align="center">0.01</td><td align="center">0.14</td><td align="center">-</td><td align="center">-</td></tr><tr><td align="left"><bold>O7</bold></td><td align="center">-</td><td align="center">-</td><td align="center">0.01</td><td align="center">0.04</td></tr><tr><td align="left"><bold>O3+4</bold></td><td align="center">0.20</td><td align="center">0.11</td><td align="center">0.30</td><td align="center">0.32</td></tr><tr><td align="left"><bold>O3+4+7</bold></td><td align="center">0.17</td><td align="center">0.01</td><td align="center">0.21</td><td align="center">0.04</td></tr><tr><td align="left"><bold>O3+4+2</bold></td><td align="center">0.34</td><td align="center">0.19</td><td align="center">0.22</td><td align="center">0.06</td></tr><tr><td align="left"><bold>O3+4+8</bold></td><td align="center">0.20</td><td align="center">0.32</td><td align="center">0.04</td><td align="center">0.30</td></tr><tr><td align="left"><bold>O3+4+23+2</bold></td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">0.01</td></tr></tbody></table></table-wrap> <table-wrap id="T5" position="float"><label>Table 5</label><caption><p>Correlation coefficients between chromosomal arrangements and high insertion frequency (HF) sites</p></caption><table frame="hsides" rules="groups"><thead><tr><th/><th/><th align="center" colspan="8">Populations</th></tr><tr><th/><th/><th colspan="8"><hr/></th></tr><tr><th/><th/><th align="center" colspan="2">DA</th><th align="center" colspan="2">BE</th><th align="center" colspan="2">MA</th><th align="center" colspan="2">BO</th></tr></thead><tbody><tr><td align="center"><bold>HF sites of bilbo</bold></td><td align="center"><bold>Arrang.</bold></td><td align="center"><bold>r</bold></td><td align="center"><bold>q-value</bold></td><td align="center"><bold>r</bold></td><td align="left"><bold>q-value</bold></td><td align="center"><bold>r</bold></td><td align="center"><bold>q-value</bold></td><td align="center"><bold>r</bold></td><td align="center"><bold>q-value</bold></td></tr><tr><td colspan="10"><hr/></td></tr><tr><td align="center">11B</td><td align="center">A2</td><td align="center">-0.13</td><td align="center">(0.42)</td><td align="center">0.35**</td><td align="left">(10<sup>-7</sup>)</td><td align="center">0.28*</td><td align="center">0.02</td><td align="center">-0.07</td><td align="center">(0.99)</td></tr><tr><td align="center">20A</td><td align="center">J1</td><td align="center">0.24</td><td align="center">(0.11)</td><td align="center">-</td><td align="left">-</td><td align="center">-0.06</td><td align="center">(0.39)</td><td align="center">-0.04</td><td align="center">(0.99)</td></tr><tr><td align="center">43B</td><td align="center">Ust</td><td align="center">-0.08</td><td align="center">(0.46)</td><td align="center">0.04</td><td align="left">0.44</td><td align="center">0.23*</td><td align="center">(0.04)</td><td align="center">0.24</td><td align="center">(0.99)</td></tr><tr><td align="center">45C</td><td align="center">Ust</td><td align="center">0.11</td><td align="center">(0.56)</td><td align="center">0.20</td><td align="left">(0.08)</td><td align="center">0.40**</td><td align="center">(3.10<sup>-7</sup>)</td><td align="center">0.05</td><td align="center">(0.99)</td></tr><tr><td/><td align="center">U1+2</td><td align="center">-0.13</td><td align="center">(0.59)</td><td align="center">0.22*</td><td align="left">(0.04)</td><td align="center">-0.19</td><td align="center">0.08</td><td align="center">0.02</td><td align="center">(0.99)</td></tr><tr><td align="center">45D</td><td align="center">Ust</td><td align="center">0.15</td><td align="center">(0.42)</td><td align="center">0.49**</td><td align="left">(3.10<sup>-7</sup>)</td><td align="center">--</td><td align="center">--</td><td align="center">--</td><td align="center">--</td></tr><tr><td/><td align="center">U1+2</td><td align="center">-0.02</td><td align="center">(0.86)</td><td align="center">-0.34**</td><td align="left">(2.10<sup>-3</sup>)</td><td align="center">--</td><td align="center">--</td><td align="center">--</td><td align="center">--</td></tr><tr><td align="center">53A</td><td align="center">Ust</td><td align="center">-0.20</td><td align="center">(0.15)</td><td align="center">-0.35**</td><td align="left">(3.10<sup>-7</sup>)</td><td align="center">-0.22*</td><td align="center">(0.04)</td><td align="center">-0.05</td><td align="center">(0.99)</td></tr><tr><td/><td align="center">U1+2</td><td align="center">0.26</td><td align="center">(0.11)</td><td align="center">0.05</td><td align="left">0.44</td><td align="center">0.18</td><td align="center">(0.08)</td><td align="center">0.07</td><td align="center">(0.99)</td></tr><tr><td colspan="10"><hr/></td></tr><tr><td align="center">57D</td><td align="center">E1+2+9</td><td align="center">-</td><td align="center">-</td><td align="center">0.25</td><td align="left">(0.12)</td><td align="center">-0.10</td><td align="center">(0.44)</td><td align="center">0.57</td><td align="center">(0.41)</td></tr><tr><td/><td align="center">E1+2+9+3</td><td align="center">-0.06</td><td align="center">(0.72)</td><td align="center">0.17</td><td align="left">(0.15)</td><td align="center">0.28*</td><td align="center">(0.04)</td><td align="center">-0.02</td><td align="center">(0.99)</td></tr><tr><td align="center">59C</td><td align="center">Est</td><td align="center">0.36**</td><td align="center">(5.10<sup>-3</sup>)</td><td align="center">-</td><td align="left">-</td><td align="center">-</td><td align="center">-</td><td align="center">-0.02</td><td align="center">(0.80)</td></tr><tr><td align="center">67A</td><td align="center">Est</td><td align="center">-0.15</td><td align="center">(0.42)</td><td align="center">-0.25*</td><td align="left">(0.04)</td><td align="center">-0.69**</td><td align="center">(4.10<sup>-8</sup>)</td><td align="center">-0.03</td><td align="center">(0.99)</td></tr><tr><td/><td align="center">E1+2+9+12</td><td align="center">0.64**</td><td align="center">(1.10<sup>-4</sup>)</td><td align="center">-</td><td align="left">-</td><td align="center">0.85**</td><td align="center">(4.10<sup>-8</sup>)</td><td align="center">0.22</td><td align="center">(0.51)</td></tr><tr><td align="center">74D</td><td align="center">E1+2+9+12</td><td align="center">0.33*</td><td align="center">(0.04)</td><td align="center">-</td><td align="left">-</td><td align="center">0.24*</td><td align="center">(0.04)</td><td align="center">0.13</td><td align="center">(0.99)</td></tr><tr><td colspan="10"><hr/></td></tr><tr><td align="center">82A</td><td align="center">Ost</td><td align="center">0.12</td><td align="center">(0.63)</td><td align="center">-0.10</td><td align="left">(0.28)</td><td align="center">0.29*</td><td align="center">(0.02)</td><td align="center">-0.16</td><td align="center">(0.99)</td></tr><tr><td align="center">83C</td><td align="center">O3+4+7</td><td align="center">0.00</td><td align="center">(0.72)</td><td align="center">0.09</td><td align="left">(0.08)</td><td align="center">-0.14</td><td align="center">(0.17)</td><td align="center">0.56</td><td align="center">(0.11)</td></tr><tr><td align="center">85A</td><td align="center">O3+4</td><td align="center">-0.14</td><td align="center">(0.46)</td><td align="center">0.27*</td><td align="left">(0.02)</td><td align="center">0.10</td><td align="center">(0.17)</td><td align="center">0.03</td><td align="center">(0.99)</td></tr><tr><td/><td align="center">O3+4+7</td><td align="center">-0.17</td><td align="center">(0.41)</td><td align="center">0.08</td><td align="left">0.53</td><td align="center">-0.24*</td><td align="center">(0.04)</td><td align="center">-0.04</td><td align="center">(0.99)</td></tr><tr><td align="center">89C</td><td align="center">O3+4</td><td align="center">-0.16</td><td align="center">0.42</td><td align="center">-0.25*</td><td align="left">(0.03)</td><td align="center">-0.21</td><td align="center">(0.06)</td><td align="center">-0.09</td><td align="center">(0.99)</td></tr><tr><td/><td align="center">O3+4+2</td><td align="center">-0.36**</td><td align="center">(5.10<sup>-3</sup>)</td><td align="center">-0.31**</td><td align="left">(7.10<sup>-3</sup>)</td><td align="center">0.09</td><td align="center">(0.22)</td><td align="center">0.07</td><td align="center">(0.99)</td></tr><tr><td/><td align="center">O3+4+8</td><td align="center">0.58**</td><td align="center">(1.10<sup>-4</sup>)</td><td align="center">0.34**</td><td align="left">(4.10<sup>-3</sup>)</td><td align="center">0.26</td><td align="center">(0.05)</td><td align="center">0.02</td><td align="center">(0.99)</td></tr><tr><td align="center">91B</td><td align="center">Ost</td><td align="center">0.28</td><td align="center">(0.16)</td><td align="center">0.02</td><td align="left">(0.12)</td><td align="center">0.32*</td><td align="center">(0.03)</td><td align="center">-0.06</td><td align="center">(0.99)</td></tr><tr><td/><td align="center">O5</td><td align="center">-0.03</td><td align="center">(0.86)</td><td align="center">0.28*</td><td align="left">(0.02)</td><td align="center">--</td><td align="center">--</td><td align="center">--</td><td align="center">--</td></tr><tr><td align="center">98D</td><td align="center">O3+4+2</td><td align="center">0.05</td><td align="center">(0.72</td><td align="center">0.12</td><td align="left">(0.18)</td><td align="center">0.32*</td><td align="center">(0.03)</td><td align="center">0.29</td><td align="center">(0.66)</td></tr><tr><td colspan="10"><hr/></td></tr><tr><td align="center"><bold>HF sites of gypsy</bold></td><td/><td/><td/><td/><td/><td/><td/><td/><td/></tr><tr><td align="center">41C</td><td align="center">U1+2</td><td align="center">0.35*</td><td/><td align="center">0.52**</td><td/><td align="center">0.25*</td><td/><td align="center">--</td><td/></tr><tr><td align="center">52D</td><td align="center">U1+2</td><td align="center">0.32*</td><td/><td align="center">0.28*</td><td/><td align="center">-0.12</td><td/><td align="center">--</td><td/></tr><tr><td align="center">63C</td><td align="center">Est</td><td align="center">--</td><td/><td align="center">--</td><td/><td align="center">-0.46**</td><td/><td align="center">--</td><td/></tr><tr><td/><td align="center">E1+2+9+12</td><td align="center">--</td><td/><td align="center">--</td><td/><td align="center">0.45**</td><td/><td align="center">--</td><td/></tr><tr><td align="center">74D</td><td align="center">E1+2+9+3</td><td align="center">0.31</td><td/><td align="center">-0.09</td><td/><td align="center">0.34</td><td/><td align="center">-0.04</td><td/></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p><bold>Poisson distribution: raw data of bilbo copy number per chromosome and population, P values and chi tests</bold>. A table of detailed tests of Poisson distribution of bilbo per chromosome and haploid genome.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p><bold>Poisson distribution: raw data of gypsy copy number per chromosome and population, P values and chi tests</bold>. A table of detailed tests of Poisson distribution of gypsy per chromosome and haploid genome.</p></caption></supplementary-material>
<table-wrap-foot><p>Population origin: Davis (DA), Bellingham (BE), Maipú (MA), Bordils (BO).</p><p>Occupancy profile: number of times that each site is occupied in populations.</p></table-wrap-foot> <table-wrap-foot><p>TE: Transposable elements; Ch: chromosome; HG: haploid genome; m: mean copy number; V<sub>n</sub>: variance of copy number; Numbers of haploid genomes are in parenthesis; DC: dispersion coefficient (V<sub>n</sub>/m); *P &lt; 0.05; **P &lt; 0.01. Bonferroni's correction was applied. See Table 1 for population origin</p></table-wrap-foot> <table-wrap-foot><p>P. chromat: Proportion of chromatin. Ga: Comparison of the proportion of TEs among chromosomes. Gb: Comparison of the proportion of TEs between chromosome A (X) and autosomes. Gc: Comparison of the proportion of TEs among autosomes; H: Heterogeneity test between colonizing populations; Df: Degrees of freedom; Pooled: Only colonizing populations; *P &lt; 0.05; **P &lt; 0.001. Bonferroni's correction was applied. Test values excluding high insertion frequency sites are in parenthesis.</p></table-wrap-foot> <table-wrap-foot><p>-: Arrangement absent</p></table-wrap-foot> <table-wrap-foot><p>Only high insertion sites showing correlation coefficient values either significant or higher than 0.20 at least in one population, are considered. Arrang: arrangement; r: Correlation coefficient; -: indicates cases where correlations cannot be computed because of low ETs copy number; --: indicates the lack of a site or an inversion in the population). See Table 1 for population origin. *P &lt; 0.05; **P &lt; 0.01. Q-value and Bonferroni corrections were applied to <italic>bilbo </italic>and <italic>gypsy </italic>respectively</p></table-wrap-foot>
<graphic xlink:href="1471-2148-8-234-1"/> <graphic xlink:href="1471-2148-8-234-2"/>
<media xlink:href="1471-2148-8-234-S1.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-234-S2.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Deininger", "Roy-Engel"], "given-names": ["PL", "AM"], "person-group": ["Craig NL, Craigie R, Gellert M, Lambowitz A"], "article-title": ["Mobile elements in animals and plant genomes"], "source": ["Mobile DNA II"], "year": ["2002"], "publisher-name": ["Washington D.C. , American Society for Microbiology"], "fpage": ["1074"], "lpage": ["1092"]}, {"surname": ["Charlesworth", "Charlesworth"], "given-names": ["B", "D"], "article-title": ["The population dynamics of transposable elements"], "source": ["Genet Res"], "year": ["1983"], "volume": ["42"], "fpage": ["1"], "lpage": ["27"]}, {"surname": ["Fontdevila"], "given-names": ["A"], "person-group": ["Hewitt GM, Johnston AWB, Young JPW"], "article-title": ["Colonizing Species of Drosophila"], "source": ["Molecular Techniques in Taxonomy"], "year": ["1991"], "volume": ["57"], "issue": ["H"], "publisher-name": ["Heidelberg , Springer-Verlag "], "fpage": ["249"], "lpage": ["269"]}, {"surname": ["Sturtevant"], "given-names": ["AH"], "article-title": ["The classification of the genus Drosophila with the description of nine new species"], "source": ["Univ Texas Publ"], "year": ["1942"], "volume": ["4213"], "fpage": ["5"]}, {"surname": ["Krimbas"], "given-names": ["CB"], "source": ["Drosophila subobscura."], "year": ["1993"], "publisher-name": ["Hamburg , Verlag Dr. Kovac"]}, {"surname": ["Ayala", "Serra", "Prevosti"], "given-names": ["FJ", "L", "A"], "article-title": ["A grand experiment in evolution: the Drosophila subobscura colonization of the Americas."], "source": ["Genome"], "year": ["1989"], "volume": ["31"], "fpage": ["246"], "lpage": ["255"]}, {"surname": ["Brncic", "Prevosti", "Budnik", "Moncl\u00fas", "Oca\u00f1a"], "given-names": ["D", "A", "M", "M", "J"], "article-title": ["Colonization of Drosophila subobscura in Chile I. First population and cytogenetic studies."], "source": ["Genetica"], "year": ["1981"], "volume": ["56"], "fpage": ["3"], "lpage": ["9"], "pub-id": ["10.1007/BF00126923"]}, {"surname": ["Beckenbach", "Prevosti", "Serra"], "given-names": ["AT", "A", "L"], "article-title": ["Colonization of North America by the European species Drosophila subobscura and D. ambigua"], "source": ["Am Midl Nat"], "year": ["1986"], "volume": ["115"], "fpage": ["10"], "lpage": ["18"], "pub-id": ["10.2307/2425832"]}, {"surname": ["Menozzi", "Krimbas"], "given-names": ["P", "CB"], "article-title": ["The inversion polymorphism of D. subobscura revisited: Synthetic maps of gene arrangement frequencies and their interpretation"], "source": ["J Evol Biol"], "year": ["1992"], "volume": ["5"], "fpage": ["625"], "lpage": ["641"], "pub-id": ["10.1046/j.1420-9101.1992.5040625.x"]}, {"surname": ["Mestres", "Bala\u00f1\u00e1", "Segarra", "Prevosti", "Serra"], "given-names": ["F", "J", "C", "A", "L"], "article-title": ["Colonization of America by Drosophila subobscura: Analysis of the O5 inversions from Europe and America and their implications for the colonizing process."], "source": ["Evolution"], "year": ["1992"], "volume": ["46"], "fpage": ["1564"], "lpage": ["1568"], "pub-id": ["10.2307/2409960"]}, {"surname": ["Fisher"], "given-names": ["RA"], "source": ["Statistical methods for research workers"], "year": ["1958"], "edition": ["13rd"], "publisher-name": ["New York , Hafner"]}, {"surname": ["Stumm-Zollinger", "Goldschmidt"], "given-names": ["E", "E"], "article-title": ["Geographical differentiation of inversion systems in Drosophila subobscura. "], "source": ["Evolution"], "year": ["1959"], "volume": ["13"], "fpage": ["89"], "lpage": ["98"], "pub-id": ["10.2307/2405946"]}, {"surname": ["Sokal", "Rohlf"], "given-names": ["RR", "FJ"], "source": ["Biometry"], "year": ["1995"], "edition": ["3rd"], "publisher-name": ["New York , Freeman"]}, {"surname": ["Balany\u00e0", "Sol\u00e9", "Oller", "Sperlich", "Serra"], "given-names": ["J", "E", "JM", "D", "L"], "article-title": ["Long-term changes in the chromosomal inversion polymorphism of Drosophila subobscura. II. European populations"], "source": ["J Zoolog Syst Evol Res"], "year": ["2004"], "volume": ["42"], "fpage": ["191"], "lpage": ["201"], "pub-id": ["10.1111/j.1439-0469.2004.00274.x"]}, {"surname": ["Lyttle", "Haymer"], "given-names": ["TW", "DS"], "article-title": ["The role of the transposable element hobo in the origin of the endemic inversions in wild populations of Drosophila melanogaster"], "source": ["Genetical Research"], "year": ["1992"], "volume": ["86"], "fpage": ["113"], "lpage": ["126"]}, {"surname": ["Prevosti"], "given-names": ["A"], "article-title": ["Chromosomal polymorphism in Drosophila subobscura populations from Barcelona (Spain)"], "source": ["Genet Res"], "year": ["1964"], "volume": ["5"], "fpage": ["27"], "lpage": ["38"]}, {"surname": ["Prevosti", "Serra", "Rib\u00f3", "Aguad\u00e9", "Sagarra", "Balany\u00e0", "Moncl\u00fas", "Garc\u00eda"], "given-names": ["A", "L", "G", "M", "E", "J", "M", "MP"], "article-title": ["The colonization of Drosophila subobscura in Chile. II. Clines in the chromosomal arrangements."], "source": ["Evolution"], "year": ["1985"], "volume": ["39"], "fpage": ["838"], "lpage": ["844"], "pub-id": ["10.2307/2408683"]}, {"surname": ["Mestres", "Pegueroles", "Prevosti", "Serra"], "given-names": ["F", "G", "A", "L"], "article-title": ["Colonization of America by Drosophila subobscura: lethal genes and the problem of the O5 inversion."], "source": ["Evolution"], "year": ["1990"], "volume": ["44"], "fpage": ["1823"], "lpage": ["1836"], "pub-id": ["10.2307/2409511"]}, {"surname": ["Maside", "Bartolom\u00e9", "Charlesworth"], "given-names": ["X", "C", "B"], "article-title": ["S-element insertions are associated with the evolution of the Hsp70 genes in Drosophila melanogaster."], "source": ["Current Biol"], "year": ["2002"], "volume": ["12"], "fpage": ["1686"], "lpage": ["1691"], "pub-id": ["10.1016/S0960-9822(02)01181-8"]}, {"surname": ["Bergman", "Quesneville", "Anxolab\u00e9h\u00e8re", "Ashburner"], "given-names": ["CM", "H", "D", "M"], "article-title": ["Recurrent insertion and duplication generate networks of transposable element sequences in the Drosophila melanogaster genome"], "source": ["Genome Biol"], "year": ["2006"], "volume": ["7"], "issue": ["11"]}, {"surname": ["Vivas", "Garc\u00eda-Planells", "Ru\u00edz", "Marfany", "Paricio", "Gonz\u00e1lez-Duarte", "de Frutos"], "given-names": ["MV", "J", "C", "G", "N", "R", "R"], "article-title": ["GEM, a cluster of repetitive sequences in the Drosophila subobscura genome"], "source": ["Genet Res"], "year": ["1999"], "volume": ["229"], "fpage": ["47"], "lpage": ["57"]}, {"surname": ["Belyaeva", "Ananiev", "Gvozdev"], "given-names": ["ES", "EV", "VA"], "article-title": ["Distribution of mobile dispersed genes (mdg-1 and mdg-3) in the chromosomes of Drosophila melanogaster ."], "source": ["Chromosoma"], "year": ["1984"], "volume": ["90"], "fpage": ["16"], "lpage": ["19"], "pub-id": ["10.1007/BF00352273"]}, {"surname": ["Bucheton", "Busseau", "Teninges"], "given-names": ["A", "I", "A"], "person-group": ["Craig NL, Craigie R, Gellert M, Lambowitz AM"], "article-title": ["I elements in Drosophila melanogaster"], "source": ["Mobile DNA II"], "year": ["2002"], "publisher-name": ["Washington DC , ASM press"], "fpage": ["796"], "lpage": ["813"]}, {"surname": ["Smidt"], "given-names": ["ER"], "person-group": ["Gr\u00fcnewald-Janho S, Keesey J, Leous M, van Miltenburg R, Schroeder C"], "article-title": ["A simplified and efficient protocol for non-radioactive in situ hybridization to polytene chromosomnes with a DIG-labeled DNA probe"], "source": ["Non-radioactive In situ hybridization Aplication Manual"], "year": ["1996"], "publisher-name": ["Mannheim , Roche"], "fpage": ["97"], "lpage": ["107"]}, {"surname": ["Benjamini", "Hochberg"], "given-names": ["Y", "Y"], "article-title": ["Controlling the false discovery rate. A practical and powerful approach to multiple testing."], "source": ["J R Stat Soc Ser B Stat Meth"], "year": ["1995"], "volume": ["57"], "fpage": ["289"], "lpage": ["300"]}, {"surname": ["Storey"], "given-names": ["JD"], "article-title": ["A direct approach to false discovery rates"], "source": ["J R Stat Soc Ser B Stat Meth"], "year": ["2002"], "volume": ["64"], "fpage": ["479"], "lpage": ["498"], "pub-id": ["10.1111/1467-9868.00346"]}, {"surname": ["Holm"], "given-names": ["S"], "article-title": ["A simple sequencially rejective multiple test procedure"], "source": ["Scand J Stat"], "year": ["1979"], "volume": ["6"], "fpage": ["65"], "lpage": ["70"]}]
BMC Evol Biol. 2008 Aug 14; 8:234