Figures
Abstract
Background
Although insertions and deletions (indels) account for a sizable portion of genetic changes within and among species, they have received little attention because they are difficult to type, are alignment dependent and their underlying mutational process is poorly understood. A fundamental question in this respect is whether insertions and deletions are governed by similar or different processes and, if so, what these differences are.
Methodology/Principal Findings
We use published resequencing data from Seattle SNPs and NIEHS human polymorphism databases to construct a genomewide data set of short polymorphic insertions and deletions in the human genome (n = 6228). We contrast these patterns of polymorphism with insertions and deletions fixed in the same regions since the divergence of human and chimpanzee (n = 10546). The macaque genome is used to resolve all indels into insertions and deletions. We find that the ratio of deletions to insertions is greater within humans than between human and chimpanzee. Deletions segregate at lower frequency in humans, providing evidence for deletions being under stronger purifying selection than insertions. The insertion and deletion rates correlate with several genomic features and we find evidence that both insertions and deletions are associated with point mutations. Finally, we find no evidence for a direct effect of the local recombination rate on the insertion and deletion rate.
Citation: Sjödin P, Bataillon T, Schierup MH (2010) Insertion and Deletion Processes in Recent Human History. PLoS ONE 5(1): e8650. https://doi.org/10.1371/journal.pone.0008650
Editor: Richard Cordaux, University of Poitiers, France
Received: August 13, 2009; Accepted: December 14, 2009; Published: January 19, 2010
Copyright: © 2010 Sjödin et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: PS was supported by a post doctoral grant (grant number CTS 06: 417, URL http://www.carltryggersstiftelse.se/) from Carl Trygger's Foundation. TB was supported by a Steno Fellowship (grant number 95092886, URL http://en.fi.dk/councils-commissions/the-danish-council-for-independent-research/) from the Danish Council for Natural Science (FNU). MHS was supported by a grant (grant number LSHC-CT-2005-018827, URL http://en.fi.dk/councils-commissions/the-danish-council-for-independent-research/) from the Danish Council for Technology and Production. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Although indels are less common than single nucleotide mutations, they generally account for the majority of differences between species [1]. Indels have also been implicated in many human diseases such as cystic fibrosis, fragile X syndrome and Huntingtons disease [2], [3] as well as in many cancers [4]. Very recently indels have been shown to influence the mutation rate of neighbouring genomic sequences [5]. In spite of this it is our view that, to date, indels have been overlooked in evolutionary studies. This can partly be attributed to difficulties of modelling indels, both because comparatively little is known about their origin and also because a model of indels has to deal with the length of an indel and not only its mutation rate. Indels are also more difficult to handle because they are alignment dependent [6], and indeed, have more often been treated as alignment noise rather than something biologically interesting (but see [7] for an exception).
When considering indels, a fundamental question that needs addressing is whether there is a need to treat insertions and deletions differently. This presents quite a challenge because correctly identifying an indel as either a deletion or an insertion event is very sensitive to alignment errors. It is also obvious that there are major differences in the evolutionary origin and dynamics of indels of small and large size. For instance, large scale indels caused by the proliferation and illegitimate recombination of transposable elements [8], [9] are clearly very different from short indels generated by polymerase slippage, as in micro-satellites [10]. Moreover, some authors suggest that deletions are more deleterious than insertions [11], [12], while others have argued that insertions may be deleterious because they increase the number of sites that can mutate into deleterious variation [13].
Another question to address when looking at indels from an evolutionary perspective is whether there is an association between indels and recombination rate. We believe that there are at least three reasons to expect such an association. The first is a consequence of recombination rates being defined as the number of crossing-over events per physical length unit and since indels by necessity affect the physical length of the sequence region they reside in, they affect this rate. The second reason is more mechanistic: given two homologous sequences the probability of initiating a recombination event in that specific region depends on the similarity of the sequences [14], [15]. In this respect, a length difference between two potentially recombining sequences may be even more likely to inhibit the initiation of a recombination event than a single nucleotide difference. We would, therefore, predict a negative correlation between recombination rates in a small region (i.e. one kb) and the heterozygosity of a polymorphic indel in the same region. The third reason is that the recombination process itself could be causing indel formations.
In this study we used polymorphic human data freely available from the University of Washington together with the human-chimpanzee-macaque alignment from the UCSC genome center. We estimate fine scale recombination rates in regions of the human genome covered by our polymorphic dataset. We then examine if the occurrence of polymorphic insertions or deletions affected the recombination in these regions. Also, by using the chimpanzee as outgroup to orient indels polymorphic in the selected human genome sequence, and the macaque genome sequence to orient fixed human and chimpanzee indels, we were able to investigate the dynamics of short (1 to 100 bp) insertions and deletions both on the intra- and the interspecific levels. Contrasting these levels, we gained information concerning differences between insertions and deletions with respect to their origin and genomic effect.
Results and Discussion
To compare polymorphic human indels with fixed human indels and chimpanzee indels (set-up shown in Figure 1), we scanned the Seattle SNPs and the NIEHS Environmental Genome Project for indels and SNPs. Then, for the homologous regions scanned for polymorphism, we compared the human, chimpanzee and macaque sequences in the 27 species multiple alignment available from the UCSC genome center to get fixed indel and single nucleotide differences between human and chimpanzee. In total we scanned 20.3 Mb of the human genome (Table 1).
Dating is approximate and drawing is not to scale. Thin blue lines denote species boundaries and thick black lines denote the gene genealogy underlying the sequence data used in that study.
Sample Frequency Distributions
The average sample frequency distributions of SNPs, insertions and deletions are shown in Table 2. In general, deletions have a slightly lower average frequency than insertions. This is also reflected in a significant Wilcoxon rank sum test: polymorphic insertions segregate at significantly higher frequencies than deletions. When comparing, separately, 1 bp insertions to 1 bp deletions and longer than 1 bp insertions to longer than 1 bp deletions, only the latter comparison is significant (Table 3). To get a more comprehensive view of these frequency differences, we also present the frequency spectrum of the derived variant of polymorphic insertions and deletions. We contrast these spectra with those obtained in the same genomic regions for: i) nonsynonymous SNPs, ii) synonymous SNPs and iii) SNPs in noncoding regions (Figure 2).
SNPs in noncoding regions always in grey. These are contrasted to synonymous SNPs (black) and nonsynonymous SNPs (white) in left histogram, to 1 bp insertions (black) and 1 bp deletions (black) in middle histogram and to longer than one bp insertions (black) and deletions (white) in right histogram.
Indel Counts, Ratios and McDonald-Kreitman-Like Tests
The number of insertions and deletions and their ratio are shown in Figure 3 for the three categories and for different lengths.
Deletion count (top), insertion count (middle) and DI-ratio (bottom) for different indel lengths (in bp) (white: polymorphic human, grey: fixed human, black: chimp). Note the difference in scale on the y-axis for indels of length 1 to 4 bp and those longer in the two top graphs. Length classes constructed such that the smallest count in any of the six possible indel categories was larger than 20 (except for indels of length 31 to 100 bp for which only the count of polymorphic human deletions and chimp insertions were larger than 20).
To test whether the distribution of polymorphic variation is different from fixed human variation we performed a - test on a 2×2 contingency table. This is in essence the same approach as used by McDonald and Kreitman [16] to contrast polymorphism and divergence in synonymous vs. non synonymous sites. We found a ratio of deletions to insertions (DI-ratio) of ∼2.4∶1 for polymorphic human indels, similar to that reported by Bhangale et al [17] and the Human Gene Mutation Database. This ratio is significantly higher than for fixed human indels (Figure 3 and Table 4). The ratio of non-synonymous to synonymous changes also mirrors this pattern (Table 4), which suggests that in these ratios, deletion events and non-synonymous changes are the more deleterious variants.
The two ratios were also marginally, but significantly, different between the fixed human and chimp categories. This significance disappears (data not shown), however, when the human and chimpanzee sequences from UCSC are compared directly (without dividing the human variation into polymorphic and fixed) as the chimp category also includes some variation that is still polymorphic within the chimpanzee population.
Genomic Features Affecting the Occurrence of Indels
Results of the linear model analysis of correlates of indels are shown in Table 5. The number of single nucleotide changes in a window strongly correlates with the number of indels and this correlation is exceptionally strong within the “polymorphic” and “fixed” categories. In other words, SNPs are strongly positively correlated with polymorphic indels but not so much to fixed human indels or chimp indels (and so on). The positive correlation between SNPs and polymorphic indels could partly be explained by variation in time to the most recent common ancestor of a window but this effect provides a poor explanation for why the number of indels on the chimp branch seem to be associated with single nucleotide substitutions on the chimp branch while fixed human indels are best explained by single nucleotide substitutions on the human branch. Furthermore, these associations seem to be weaker for AT to GC changes than other single nucleotide changes. This could be a consequence of there being fewer AT to GC changes than other changes, although the difference is only twofold (data not shown). Another strong explanatory factor in common for all indel categories is the presence of poly(A/T).
The number of indels in a window is in general negatively correlated with the percentage of sites in coding exons in a window, but positively correlated with the percentage sites in UTR. Of repetitive elements, only the percentage sites in SINEs is a significant explanatory variable in our linear model. This variable is negatively correlated to all categories of indels.
Although insertions and deletions are generally correlated with the same genomic features, there are some notable exceptions: the highly significant interaction term between indel type and GC-content and poly (A/T) for the polymorphic human and chimp data categories. Here the interpretation of the positive effect of the interaction terms between indel type and GC content and poly(A/T) is not straightforward since this indicates one of two possibilities: either i) a stronger positive correlation of GC content (and/or poly(A/T)) with insertions than with deletions, or ii) a weaker negative correlation of GC (and/or poly(A/T)) with insertions than with deletions. Closer inspection (by constructing separate models for insertions and for deletions) reveals that i) is the cause of the positive effect of “ID: poly AT” while ii) is the reason for the positive effect of “ID: GC” (data not shown). This also showed that the marginally significant effect of the interaction between indel type and logged recombination rate is due to a positive correlation between the number of polymorphic insertions and the logged recombination rate in our data.
Indels and Recombination
The effect of adding heterozygosity and the length effect of polymorphic insertions and deletions on a trained linear model (see Table S1 for specifications of the trained model) of logged recombination rate are shown in Table 6. A model using only windows with insertions does not show a significant improvement over the trained model but the recombination rate is significantly better predicted in the model fitted only on windows with polymorphic deletions. This could be a consequence of there being many more windows with deletions than windows with insertions. However, a model using all windows with indels (instead of treating insertions and deletions separately), was not significant (data not shown). In the deletion model, we find a negative correlation of length effect and recombination, as predicted. A positive effect of heterozygosity was also found although it was not predicted under our hypothesis.
Deletions Are More Deleterious Than Insertions
In many ways the difference between insertion and deletion counts mirrors the difference between synonymous and non-synonymous single nucleotide mutation counts (Table 4). The DI-ratio is much higher for polymorphic human indels than fixed human indels in our data (Figure 3 and Table 4). A sizable proportion of all polymorphisms in the human genome (SNPs, indels and CNVs) are expected to be weakly deleterious and not destined to become fixed. Accordingly, the ratio of non-synonymous to synonymous SNPs is expected to be much higher than the ratio for non-synonymous to synonymous fixed single nucleotide sites. Similarly, the DI-ratio for polymorphic human indels is much larger than for fixed human indels and can, therefore, be interpreted as stronger selection against deletions than insertions.
Additional evidence for deletions being more deleterious than insertions is provided by their sample frequency distributions: polymorphic deletions segregate at significantly lower frequencies than polymorphic insertions (Table 2 and 3). Stronger selection against deletions than insertions has also been suggested in several earlier studies [11], [12], [18], [19]. One explanation is based on deletions requiring two cut points, while an insertion only has one [11]. Briefly, if an important motif resides in sequence positions n1 to n2, an insertion of any length at position n1 to n2−1 will disrupt this motif. Deletions at these positions will also disrupt the motif but additionally, so will deletions of length ≥k at start position n1−k. This explanation to why deletions may be more deleterious than insertions predicts, first of all, that the detrimental effect of deletions is more length dependent than the detrimental effect of insertions. Our findings that the mean frequency of 1 bp insertions and 1 bp deletions are not significantly different while, despite considerably fewer data points, the mean frequency of deletions longer than 1 bp is significantly smaller than insertions longer than 1 bp (Table 3) thus provide support for this explanation. Moreover (this point is due to an insightful reviewer), assuming that the density of important motifs decreases with distance to coding regions, the difference in selection between the two indels types should be larger close to coding regions than far away. In our data this is reflected by a larger mean and variance of distance to the closest coding region (Figure 4), and also in a stronger negative correlation between sample frequency and distance to coding regions, for deletions than for insertions (Table 7).
Boxplots of distance to the closest coding region of polymorphic deletions (“poly del”, n = 4350), polymorphic insertions (“poly ins”, n = 1831), fixed human deletions (“hum del”, n = 5824) , fixed human insertions (“hum ins”, n = 4789) , chimp deletions (“pan del”, n = 7527) and chimp insertions (“pan ins”, n = 3929). Indels within coding regions were excluded.
Differences between the Origin of Insertions and the Origin of Deletions
If deletions are more deleterious than insertions this suggests that the DI-ratio also correlates with functional constraint across the genome. If this is true, the DI-ratio may provide an alternative method to locate functionally important intergenic regions. However, for this to work it has to be assumed that insertions occur at the same rate as deletions. In our study we found very little that differentiated insertions from deletions except that poly(A/T) and GC content were significantly more correlated with insertions than deletions. Both a high poly(A/T) and a high GC content result in an increased ‘repetitiveness’ and thus the propensity of polymerase slippage. One explanation of this difference is that insertions are relatively more likely to be the outcome of an indel event caused by polymerase slippage than other indel causing events. This is very much related to the dynamics of microsatellites as microsatellites by definition are low complexity regions. We do not wish to specifically discuss microsatellites here (see for instance [10] for a review). Suffice it to say that our data is consistent with the numerous reports that microsatellites tend to expand ( = insert repeats) when short but contract ( = delete repeats) when long and that although our data certainly contains some microsatellites, our filtering procedure should have excluded the majority of long (hyper mutable) microsatellites (see attached material “Data filtering and possible biases” for more details).
While some studies find no association between recombination and indels [5], [20], other authors have reported that insertions are associated with factors linked to recombination while deletions are mostly associated with replication-related features [21]. Our analysis provides little evidence to suggest that recombination has a differential role with respect to insertions and deletions. Although we find a marginally significant positive effect of recombination rate with polymorphic insertions, other expected patterns, assuming recombination promotes insertions but not deletions, such as a significant interaction term between indel type and number of AT to GC mutations (AT to GC mutations should be overrepresented in regions with high recombination rate [22]), are absent (Table 5). That recombination may be positively, but equally, associated with insertion and deletion formation is, however, indirectly supported by our linear model analysis. For instance, the significant positive correlation between UTR content and (both kinds of) indels is possibly a result of an increased recombination rate. Since recombination is known to be suppressed in coding regions, but elevated in flanking regions [23] this could be interpreted as indirect evidence for recombination and indels being positively correlated. This being said, the correlation is rather weak and only significant for fixed human and chimp indels. Moreover, our analysis suggests that any direct causal link between short indels and recombination can be neglected. Where an underestimate of recombination rate in windows with polymorphic indels is by construction a necessity, the effect of this is very weak (Table 6). This is likely due to a lack of statistical power as the variance of the recombination process is known to be very large [24], whereas the underestimate caused by the short indels we study is very small. In fact, given that the distance between SNPs with a polymorphic indel between them is on average shorter than the distance given in the databases (an indel in the Seattle SNPs and NIEHS databases is always represented by the long indel variant regardless of whether it is an insertion or a deletion), we underestimate the recombination rate. Given the average sizes of insertions and deletions in our data, we estimate roughly that the effect of this on our linear model of recombination rate should be smaller than 0.5% of the variance (see Text S1).
Based on a subset of the data analyzed here, Bhangale and coauthors [17] suggested that gene-conversion events between two regions may, compared to SNP differences, be particularly suppressed when there are indel differences. Their suggestion was an attempt to explain their interesting finding that while SNP diversity is significantly greater, indel diversity is significantly lower in repeat masked regions compared to other regions. Although a negative association between repetitive regions (more specifically SINEs) and indels exists in our data, the suggestion that recombination is greatly inhibited by short length differences between sequences seems to be more or less ruled out by our results (Table 6).
The Occurrence of Indels Is Associated with Single Nucleotide Changes
Our linear model of indels suggests that human SNPs are more closely associated with polymorphic human indels than with fixed human indels for which the strongest association is instead with fixed human single nucleotide differences. Likewise, chimp indels are more strongly correlated with single nucleotide changes on the chimp branch than with any such events on the human branch. In other words, indel events and single nucleotide changes are positively correlated. We interpret these correlations within data categories in three ways: either this is a demonstration of the recent finding that indels are mutagenic [5], or there are mutational hotspots along the genome affecting both single nucleotide changes and indel formation (these would then, similar to recombination hotspots [25], be transient over time), or, as suggested by a reviewer, it is an effect of a local variation in effective population size due to background selection or selective sweeps.
Underrepresentation of Indels on the X-Chromosome
Kvikstad and coauthors [21] noted that indels are strongly underrepresented on the X-chromosome, which they suggested could be evidence for indel formation by replication errors: the X-chromosome spends less time in a male background, and is thus less affected by cell divisions, than are autosomes. This interpretation is not corroborated by our results as only polymorphic indels are strongly underrepresented on the X-chromosome. Instead we propose that since the amount of polymorphism expected is 4Neμ, where Ne is the effective population size and μ is the mutation rate, the lower number of indels could be a consequence of the lower effective population size of the X-chromosome compared to autosomes. There are however several important differences between their study and ours that should be pointed out: their choice of window size is more than 2500 times larger than ours; they do not separate human indels into polymorphic and fixed; they do not consider indels in the chimpanzee lineage; they only study indels in interspersed ancestral repeats and their recombination rates are estimated on a much larger scale than ours (data from Kong et al's map [26]). It should also be mentioned that their data set is an order of magnitude larger than ours.
Conclusions
What we did find from using human polymorphism and human chimp divergence data was strong evidence for deletions being more deleterious than insertions, and contrary to earlier studies, we suggest that recombination does not play a differential role with respect to insertions and deletions. More generally, and although there definitely exist differences between insertions and deletions [20], [27], our analysis is more in line with earlier findings that these differences may be on a very local scale and can be ignored [28].
Over all, from the data analyzed during this study we believe that by studying the distribution of lengths between indels [7] and by comparing DI-ratios across genomic regions, indels provide unique information for quantifying the amount of functionally important material in different classes of non-coding regions.
Materials and Methods
Data Sets
Sequence data for this study was obtained from Seattle SNPs (http://pga.mbt.washington.edu/) and the NIEHS Environmental Genome Project (http://egp.gs.washington.edu/), both made available by the University of Washington. The Seattle SNPs data consists of a 9.9 Mb sequence divided among 301 candidate genes associated with inflammatory response, while the NIEHS dataset consists of 27.5 Mb in 616 genes thought to be involved in environmental response. For genes in both datasets, some parts were not resequenced; these are typically repeat masked regions. In total 21.0 Mb (56.3%) out of 37.4 Mb, were effectively resequenced. We also used the 27 species multiple alignment available from the UCSC genome center (http://genome.ucsc.edu/) to orient variation. After removing genes for which mapping of positions in the data from the University of Washington to positions in the UCSC alignment was problematic, we were left with 891 genes and 20.3 Mb of resequenced data (Table 1).
Derived and Ancestral States
Using the resequenced data from the University of Washington as a starting point, we retrieved the homologous human-chimpanzee-macaque alignment from the UCSC genome center. We categorized any variation within the data from the University of Washington as polymorphic. Differences between the human and chimp genomes in the human-chimpanzee-macaque alignment that 1) overlapped with the resequenced regions and 2) did not overlap with polymorphic variation were designated as fixed differences. The human-chimp alignment was used to orient human polymorphisms into ancestral and derived states by assuming that the chimp variant was ancestral. In the same way, fixed differences between human and chimp were oriented using the human-chimp-macaque alignment. Orientation is straightforward for segregating sites but slightly more complex for indels. A fixed indel difference can be classified into one of four categories; 1) an insertion on the human branch, 2) a deletion on the human branch, 3) an insertion on the chimpanzee branch, and 4) a deletion on the chimpanzee branch. For instance, if the chimp sequence is the long variant this implies that either a human deletion or a chimp insertion occurred and we need the macaque sequence to resolve this (see Figure 1).
By way of this orientation procedure, we categorized variation as: i) polymorphic human data; that is variation found within the Seattle SNPs or NIEHS databases, ii) fixed human data; events in the human lineage since the human-chimpanzee split not polymorphic in the Seattle SNPs or NIEHS databases, and iii) chimp data; events occurring in the chimpanzee lineage since the human-chimpanzee split (Figure 1). We use this denotation of the categories throughout.
Frequency of Derived Polymorphic Variation
For each SNP and indel polymorphic within the Seattle SNPs and NIEHS databases the genotypes of the individuals in the samples are provided. Restricting the data to variation with only two variants where one of them is the same as the variant found in the outgroup (see above), the frequency of the derived variant was calculated by dividing the count of the derived variant by the total number of identified variants. In obtaining the frequency spectrum a potential issue is the fact that samples from distinct population are pooled within the Seattle SNPs and NIEHS databases. Moreover, not all polymorphisms were systematically amplified for all individuals in a sample. We investigated the influence of this by binning frequencies in different ways but saw no qualitative differences and we settled for 10 frequency intervals.
Filtering
First, we corrected some obvious minor errors in the NIEHS and Seattle SNPs data. Alignment artifacts may be responsible for complex patterns of fixed segregating sites and indels overlapping SNPs and polymorphic indels. We conservatively removed such problematic regions in the following way: For SNPs, we required that at least 12 out of the 16 directly adjacent positions should be identical in the human chimp alignment. As the state of indels is much more sensitive to alignment problems, we used a more stringent qualification procedure for these. For polymorphic indels, we required that 28 out of 32 flanking positions were identical in the human-chimpanzee alignment and also that there were no gaps at these positions. We also required that either all human sites were gaps or that all chimp sites were gaps within the indel sites in the human-chimpanzee alignment. For fixed indels, we added the additional requirements that 20 out of 32 directly adjacent flanking positions were identical in the human-chimpanzee-macaque alignment and that there were no gaps at these positions. Again, no ambiguous information was allowed within the indel sites in the human-chimpanzee-macaque alignment. To avoid stronger selection against deletions than against insertions, the only states we considered within the indel sites were ‘gap’ or ‘not gap’ - we did not take into account whether nucleotides within the indel sites matched each other between species. While the complete data cleaning procedure leads to loss of data we believe that data are discarded in an unbiased way with respect to SNP-type and indel-type (see Text S2).
Estimating Recombination Rate
We used the SNPs provided in the ‘.prettybase.txt’ files provided in the bulk download of both databases with at most 10% missing data as input to the package LDhat [24] to estimate the recombination rate along the genes. Two different programs to estimate recombination rate are available in LDhat: inter [29] and rho [30]. The program rho fits a recombination hotspot process on top of a varying background recombination rate while the recombination hotspot process is not implemented in inter [30]. We ran inter with the recommended penalty of 5 and rho with default parameter values. Since there were no qualitative differences between the estimates from the two programs we report only results using inter.
Statistical Models
The LDhat program used an average recombination rate resolution of 386 bp. We therefore divided sequences into windows of length 386 bp, and calculated the recombination rate in each window as follows. Given that a specific standardized window is overlapped by n LDhat-windows with (overlapping) lengths L1, L2,…, Ln, and estimated recombination rates R1, R2,…, Rn, the average recombination rate of the standardized window was calculated as for i = 1,‥,n. Finally we discarded windows where the percentage of resequenced sites was lower than 90%.
Linear models of indels.
We constructed generalized linear models with the counts of polymorphic human, fixed human and chimp indels in each window as response variables (three models in total). The explanatory variables we used were the log transformed recombination rate, the GC content, CpG count, the number of bases that are part of a poly (A/T) tract (defined as a contiguous stretch of at least 4 bp with only A's or only T's), the percentage sites in coding regions, the percentage of sites in UTR (in Genbank the designated mRNA which is not coding) and the percentage of sites not scanned. The percentage of sites in repeat masked regions was divided into three categories depending on whether the region was masked due to the presence of SINEs, LINEs or something else. A chromosome X indicator variable and the physical distance of each window to the telomere/chromosome end were additionally used as explanatory variables. Finally, we included the number of single nucleotide differences in each window. These were divided into polymorphic human, fixed human and chimp and also according to whether they were a change of type AT→GC or something else. Instead of having separate models of insertions and deletions, we included an indel type indicator (ID, for insertion or deletion) as an explanatory variable to indicate whether the number of deletions or the number of insertions were to be predicted by the linear model. Each window was used twice: once to predict the number of insertions and once to predict the number of deletions. All interaction terms with this indicator and the explanatory variables above were included. In this way, if one of these variables affected insertions and deletions significantly differently, this should result in a significant interaction term between this variable and the indicator variable. As a consequence of the way we constructed this indicator variable, a positive coefficient of an interaction term between the indel indicator variable and another explanatory variable, X, means that X should have a larger coefficient to predict insertions than deletions while a negative value of the coefficient of the interaction term implies the opposite, that the coefficient of X should be larger for deletions than for insertions. Using a - test, each of the three GLM models were finally compared to a corresponding reduced model without the indicator variable.
Linear models of recombination.
When investigating how indels relate to the log-transformed recombination rate, all explanatory variables used in the prediction of the indel counts described above were used except for the recombination rate, variables related to fixed human variation and chimp variation. Polymorphic human indels were also treated separately. We trained a model using all windows without polymorphic indels (42214 windows). To do so, all two-way interactions were initially included and the model was subsequently reduced using the stepwise procedure implemented in R [31]. To investigate separately different aspects of polymorphic indels, this model was applied to 1): a data set consisting of all windows containing insertions (1642 windows), and 2): a data set consisting of all windows containing deletions (3837 windows). In both cases, the response variable was the log transformed recombination rate in a window minus the predicted value of the trained model for this window.
Since the average distance between two sequence positions is shorter in an alignment including a polymorphic indel, we expected a negative length effect of indels on recombination. In order to test for this, we used (1−p)l, where p is the frequency of the long variant and l is the length of the indel as explanatory variable (see Text S1). If a length difference between two sequences influences the probability of a recombination event, such an effect would correlate with how often the length variants occur together in an individual. Hence, p(1−p) was used as a variable to search for this type of effect which we label by heterozygosity (see Text S1). When there was more than one deletion, or more than one insertion per window, the maximum of the length effect and heterozygosity value was used. Finally, for the data set with deletions, the length effect and heterozygosity was calculated only from deletions and, likewise, deletions were ignored for the insertion data. All statistical computations described below were performed using the statistical computing language R [31].
Supporting Information
Text S1.
Indel recombination model. Derivation of the parameters ‘indel heterozygosity’ and ‘length effect’ used in our linear model of recombination.
https://doi.org/10.1371/journal.pone.0008650.s001
(0.09 MB DOC)
Text S2.
Data filtering and possible biases. Discussion of possible non-biological sources of our results.
https://doi.org/10.1371/journal.pone.0008650.s002
(0.27 MB DOC)
Table S1.
Trained recombination model. Table showing specifications of the trained recombination model.
https://doi.org/10.1371/journal.pone.0008650.s003
(0.08 MB DOC)
Acknowledgments
We are grateful to Gerton Lunter, Mathieu Emily and Julien Dutheil for statistical support and Claire Neesham for improving the English of this manuscript.
Author Contributions
Conceived and designed the experiments: PS TB MHS. Analyzed the data: PS TB MHS. Wrote the paper: PS TB MHS.
References
- 1. Britten RJ, Rowen L, Williams J, Cameron RA (2003) Majority of divergence between closely related DNA samples is due to indels. Proc Natl Acad Sci U S A 100: 4661–4665.
- 2. Ashley CT Jr, Warren ST (1995) Trinucleotide repeat expansion and human disease. Annu Rev Genet 29: 703–728.
- 3. Zielenski J, Tsui LC (1995) Cystic fibrosis: genotypic and phenotypic variations. Annu Rev Genet 29: 777–807.
- 4. Duval A, Hamelin R (2002) Mutations at coding repeat sequences in mismatch repair-deficient human cancers: toward a new concept of target genes for instability. Cancer Res 62: 2447–2454.
- 5. Tian D, Wang Q, Zhang P, Araki H, Yang S, et al. (2008) Single-nucleotide mutation rate increases close to insertions/deletions in eukaryotes. Nature 455: 105–108.
- 6. Lunter G, Rocco A, Mimouni N, Heger A, Caldeira A, et al. (2008) Uncertainty in homology inferences: Assessing and improving genomic sequence alignment. Genome Research 18: 298–309.
- 7. Lunter G, Ponting CP, Hein J (2006) Genome-wide identification of human functional DNA using a neutral indel model. PLoS Comput Biol 2: e5.
- 8. Devos KM, Brown JK, Bennetzen JL (2002) Genome size reduction through illegitimate recombination counteracts genome expansion in Arabidopsis. Genome Res 12: 1075–1079.
- 9. Kazazian HH Jr (2004) Mobile elements: drivers of genome evolution. Science 303: 1626–1632.
- 10. Ellegren H (2004) Microsatellites: simple sequences with complex evolution. Nat Rev Genet 5: 435–445.
- 11. Petrov DA (2002) Mutational equilibrium model of genome size evolution. Theor Popul Biol 61: 531–544.
- 12. Taylor MS, Kai C, Kawai J, Carninci P, Hayashizaki Y, et al. (2006) Heterotachy in mammalian promoter evolution. PLoS Genet 2: e30.
- 13.
Lynch M (2007) The Origins of Genome Architecture. Sunderland, MA, USA: Sinauer Associates.
- 14. Chen W, Jinks-Robertson S (1999) The role of the mismatch repair machinery in regulating mitotic and meiotic recombination between diverged sequences in yeast. Genetics 151: 1299–1313.
- 15. Opperman R, Emmanuel E, Levy AA (2004) The effect of sequence divergence on recombination between direct repeats in Arabidopsis. Genetics 168: 2207–2215.
- 16. McDonald JH, Kreitman M (1991) Adaptive protein evolution at the Adh locus in Drosophila. Nature 351: 652–654.
- 17. Bhangale TR, Rieder MJ, Livingston RJ, Nickerson DA (2005) Comprehensive identification and characterization of diallelic insertion-deletion polymorphisms in 330 human candidate genes. Hum Mol Genet 14: 59–69.
- 18. Kim J, He X, Sinha S (2009) Evolution of Regulatory Sequences in 12 Drosophila Species. PLoS Genet 5: e1000330.
- 19. Taylor MS, Ponting CP, Copley RR (2004) Occurrence and Consequences of Coding Sequence Insertions and Deletions in Mammalian Genomes. Genome Research 14: 555–566.
- 20. Messer PW, Arndt PF (2007) The Majority of Recent Short DNA Insertions in the Human Genome Are Tandem Duplications. Mol Biol Evol 24: 1190–1197.
- 21. Kvikstad EM, Tyekucheva S, Chiaromonte F, Makova KD (2007) A macaque's-eye view of human insertions and deletions: differences in mechanisms. PLoS Comput Biol 3: 1772–1782.
- 22. Galtier N, Piganeau G, Mouchiroud D, Duret L (2001) GC-content evolution in mammalian genomes: the biased gene conversion hypothesis. Genetics 159: 907–911.
- 23. Myers S, Bottolo L, Freeman C, McVean G, Donnelly P (2005) A fine-scale map of recombination rates and hotspots across the human genome. Science 310: 321–324.
- 24. McVean G, Awadalla P, Fearnhead P (2002) A coalescent-based method for detecting and estimating recombination from gene sequences. Genetics 160: 1231–1241.
- 25. Winckler W, Myers SR, Richter DJ, Onofrio RC, McDonald GJ, et al. (2005) Comparison of fine-scale recombination rates in humans and chimpanzees. Science 308: 107–111.
- 26. Kong A, Gudbjartsson DF, Sainz J, Jonsdottir GM, Gudjonsson SA, et al. (2002) A high-resolution recombination map of the human genome. Nat Genet 31: 241–247.
- 27. Tanay A, Siggia E (2008) Sequence context affects the rate of short insertions and deletions in flies and primates. Genome Biology 9: R37.
- 28. Ball EV, Stenson PD, Abeysinghe SS, Krawczak M, Cooper DN, et al. (2005) Microdeletions and microinsertions causing human genetic disease: common mechanisms of mutagenesis and the role of local DNA sequence complexity. Hum Mutat 26: 205–213.
- 29. McVean GA, Myers SR, Hunt S, Deloukas P, Bentley DR, et al. (2004) The fine-scale structure of recombination rate variation in the human genome. Science 304: 581–584.
- 30. Auton A, McVean G (2007) Recombination rate estimation in the presence of hotspots. Genome Res 17: 1219–1227.
- 31.
R_development_core_team (2007) R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing.