Characteristics, causes and evolutionary consequences of male-biased mutation

Genic capture, sex linkage, and the heritability of fitness

Several sexual selection models predict that females will obtain indirect genetic benefits by preferentially mating with males that transmit high-quality genes to their offspring. However, despite widespread observations of additive population genetic variation for fitness as well as for male sexually selected traits, estimated fitness associations between fathers and offspring are often weak. Perhaps more puzzling, the strength of these associations differs drastically between species, leading many researchers to question the relevance of genetic benefits for processes of sexual selection. Here, I show that a species' sex chromosome system can strongly influence the genetic architecture of male and female fitness variation and, consequently, the heritability of fitness between fathers and their offspring. Indirect genetic benefits are reduced, and sexually antagonistic costs are pronounced, in species with X chromosomes relative to species with homomorphic sex chromosomes, environmental sex determination, or Z chromosomes. Data from the sexual selection literature are consistent with predictions of the models, though additional studies will be required to circumvent phylogenetic nonindependence between sex determination systems. This study strongly suggests that inferences about genetic benefits of female choice must be considered within a species-specific genomic context, and it has several implications for the evolution of female preferences as well as the genomic consequences of sex and sexual selection.

Comparative genomics based on massive parallel transcriptome sequencing reveals patterns of substitution and selection across 10 bird species

Next-generation sequencing technology provides an attractive means to obtain large-scale sequence data necessary for comparative genomic analysis. To analyse the patterns of mutation rate variation and selection intensity across the avian genome, we performed brain transcriptome sequencing using Roche 454 technology of 10 different non-model avian species. Contigs from de novo assemblies were aligned to the two available avian reference genomes, chicken and zebra finch. In total, we identified 6499 different genes across all 10 species, with approximately 1000 genes found in each full run per species. We found evidence for a higher mutation rate of the Z chromosome than of autosomes (male-biased mutation) and a negative correlation between the neutral substitution rate (d(S)) and chromosome size. Analyses of the mean d(N)/d(S) ratio (omega) of genes across chromosomes supported the Hill-Robertson effect (the effect of selection at linked loci) and point at stochastic problems with omega as an independent measure of selection. Overall, this study demonstrates the usefulness of next-generation sequencing for obtaining genomic resources for comparative genomic analysis of non-model organisms.

Female-to-male breeding ratio in modern humans-an analysis based on historical recombinations

Was the past genetic contribution of women and men to the current human population equal? Was polygyny (excess of breeding women) present among hominid lineages? We addressed these questions by measuring the ratio of population recombination rates between the X chromosome and the autosomes, rho(X)/rho(A). The X chromosome recombines only in female meiosis, whereas autosomes undergo crossovers in both sexes; thus, rho(X)/rho(A) reflects the female-to-male breeding ratio, beta. We estimated beta from rho(X)/rho(A) inferred from genomic diversity data and calibrated with recombination rates derived from pedigree data. For the HapMap populations, we obtained beta of 1.4 in the Yoruba from West Africa, 1.3 in Europeans, and 1.1 in East Asian samples. These values are consistent with a high prevalence of monogamy and limited polygyny in human populations. More mutations occur during male meiosis as compared to female meiosis at the rate ratio referred to as alpha. We show that at alpha not equal 1, the divergence rates and genetic diversities of the X chromosome relative to the autosomes are complex functions of both alpha and beta, making their independent estimation difficult. Because our estimator of beta does not require any knowledge of the mutation rates, our approach should allow us to dissociate the effects of alpha and beta on the genetic diversity and divergence rate ratios of the sex chromosomes to the autosomes.

A conceptual framework for the colonisation of urban areas: the blackbird Turdus merula as a case study

Despite increasing interest in urban ecology the factors limiting the colonisation of towns and cities by species from rural areas are poorly understood. This is largely due to the lack of a detailed conceptual framework for this urbanisation process, and of sufficient case studies. Here, we develop such a framework. This draws upon a wide range of ecological and evolutionary theory and the increasing number of studies of how the markedly divergent conditions in urban and rural areas influence the traits of urban populations and the structure of urban assemblages. We illustrate the importance of this framework by compiling a detailed case study of spatial and temporal variation in the urbanisation of the blackbird Turdus merula. Our framework identifies three separate stages in the urbanisation process: (i) arrival, (ii) adjustment, and (iii) spread. The rate of progress through each stage is influenced by environmental factors, especially human attitudes and socio-economic factors that determine the history of urban development and the quality of urban habitats, and by species' ecological and life-history traits. Some traits can positively influence progression through one stage, but delay progression through another. Rigorous assessment of the factors influencing urbanisation should thus ideally pay attention to the different stages. Urbanisation has some similarities to invasion of exotic species, but the two clearly differ. Invasion concerns geographic range expansion that is external to the species' original geographic range, whilst urbanisation typically relates to filling gaps within a species' original range. This process is exemplified by the blackbird which is now one of the commonest urban bird species throughout its Western Palearctic range. This is in stark contrast to the situation 150 years ago when the species was principally confined to forest. Blackbird urbanisation was first recorded in Germany in 1820, yet some European cities still lack urban blackbirds. This is especially so in the east, where urbanisation has spread more slowly than in the west. The timing of blackbird urbanisation exhibits a marked spatial pattern, with latitude and longitude explaining 76% of the variation. This strong spatial pattern contrasts with the weaker spatial pattern in timing of urbanisation exhibited by the woodpigeon Columba palumbus (with location explaining 39% of the variation), and with the very weak spatial pattern in timing of black-billed magpie Pica pica urbanisation (in which location explains 12% of the variation). Strong spatial patterns in timing of urbanisation are more compatible with the leap-frog urbanisation model, in which urban adapted or imprinted birds colonise other towns and cities, than with the independent urbanisation model, in which urban colonisation events occur independently of each other. Spatial patterns in isolation do not, however, confirm one particular model. Factors relating to the arrival and adjustment stages appear particularly likely to have influenced the timing of blackbird urbanisation. Spatial variation in the occurrence of urban populations and the timing of their establishment creates opportunities to assess the factors regulating urbanisation rates, and how the composition of urban assemblages develops as a result. These are major issues for urban ecology.

Timing of replication is a determinant of neutral substitution rates but does not explain slow Y chromosome evolution in rodents

Mutation rates, assayed as substitution rates of putatively neutral sites, are highly variable around mammalian genomes: There is heterogeneity between genes, between autosomes, and between X, Y, and autosomes. The differences between X, Y, and autosomes are typically assumed to reflect the greater number of cell divisions in the male germ-line. Such an effect can neither account for within-autosome differences nor does it predict the differences between X, Y, and autosome observed in rodents. It has recently been proposed that in primates, the time during S-phase when a gene is replicated is an important determinant of neutral rates of evolution. Here we ask 1) whether we can replicate this result in rodents, 2) whether different autosomes replicate on average at different times, and 3) whether this might explain differences in their substitution rates. Finally we ask 4) whether X, Y, and autosome replicate at different times and 5) whether any difference might explain why the number of replication events alone cannot explain their substitution rates. We find that, as in primates, autosomal intronic rates of evolution increase significantly during S-phase. Different autosomes do have different average replication times, and together with rearrangement, this is a significant predictor of between-autosome differences in substitution rate. Although we find that autosomal, X-, and Y-linked genes replicate at different times, it is paradoxical that the Y-linked genes replicate latest, and replicate more often, but are not especially fast evolving. These results support the hypothesis that replication timing is an important source of substitution rate heterogeneity.

Sex differences in intraventricular hemorrhage rates among very low birth weight newborns

BACKGROUND

The influence of male or female sex on newborn outcomes has been recognized for >30 years. Several studies have observed higher mortality and morbidity in males than in females. It is not clear how this sex difference is sustained in postnatal complications such as intraventricular hemorrhage (IVH), especially in very low birth weight (VLBW) newborns.

OBJECTIVE

This study examined possible sex-related differences in IVH rates among VLBW neonates.

METHODS

In a retrospective observational study conducted in Hospital Privado, Córdoba, Argentina, data from 332 consecutive VLBW newborns in a 12-year period were reviewed. Maternal factors, labor and delivery characteristics, and neonatal parameters, including the results of cranial ultrasound examination to detect IVH, were compared for males and females. Bivariate and multivariate logistic regression analyses were performed.

RESULTS

A total of 322 VLBW newborns were included, 168 males and 154 females. Compared with female neonates, male neonates had a higher risk of overall IVH (26.8% vs 9.7%; odds ratio [OR] = 3.4 [95% CI, 1.8-6.4]; P < 0.001) and for grades III or IV on the Papile scale (16.1% vs 1.9%; OR = 9.6 [95% CI, 2.9-32.5]; P < 0.001). In the multivariate logistic regression model, male sex sustained the association with a greater risk of IVH (OR = 6.8 [95% CI, 3.8-12.0]).

CONCLUSIONS

IVH was significantly associated with male sex in these VLBW newborns. Because other factors affect these differences, further research is required.

Effect of male age on sperm traits and sperm competition success in the guppy (Poecilia reticulata)

Deleterious mutations can accumulate in the germline with age, decreasing the genetic quality of sperm and imposing a cost on female fitness. If these mutations also affect sperm competition ability or sperm production, then females will benefit from polyandry as it incites sperm competition and, consequently, minimizes the mutational load in the offspring. We tested this hypothesis in the guppy (Poecilia reticulata), a species characterized by polyandry and intense sperm competition, by investigating whether age affects post-copulatory male traits and sperm competition success. Females did not discriminate between old and young males in a mate choice experiment. While old males produced longer and slower sperm with larger reserves of strippable sperm, compared to young males, artificial insemination did not reveal any effect of age on sperm competition success. Altogether, these results do not support the hypothesis that polyandry evolved in response to costs associated with mating with old males in the guppy.

Mosaicism and haemophilia

Mosaicism may affect the haemophilia phenotype. Well-known instances include chromosomal mosaicism due to aneuploidy and pseudo-mosaicism due to varying patterns of X-chromosome inactivation. Chromosomal mosaicism in a chimera is a potential source of phenotypic variation. Gene mosaicism is commonplace. Its pattern and effect depend on the stage of development at which a mutation occurs. Proven or possible genetic mosaicism is an important consideration when predicting the likelihood of transmission of haemophilia to a future generation.

Contrasting evolutionary dynamics between angiosperm and mammalian genomes

Continuing advances in genomics are revealing substantial differences between genomes of major eukaryotic lineages. Because most data (in terms of depth and phylogenetic breadth) are available for angiosperms and mammals, we explore differences between these groups and show that angiosperms have less highly compartmentalized and more diverse genomes than mammals. In considering the causes of these differences, four mechanisms are highlighted: polyploidy, recombination, retrotransposition and genome silencing, which have different modes and time scales of activity. Angiosperm genomes are evolutionarily more dynamic and labile, whereas mammalian genomes are more stable at both the sequence and chromosome level. We suggest that fundamentally different life strategies and development feedback on the genome exist, influencing dynamics and evolutionary trajectories at all levels from the gene to the genome.

Understanding what determines the frequency and pattern of human germline mutations

Surprising findings about human germline mutation have come from applying new technologies to detect rare mutations in germline DNA, from analysing DNA sequence divergence between humans and closely related species, and from investigating human polymorphic variation. In this Review we discuss how these approaches affect our current understanding of the roles of sex, age, mutation hot spots, germline selection and genomic factors in determining human nucleotide substitution mutation patterns and frequencies. To enhance our understanding of mutation and disease, more extensive molecular data on the human germ line with regard to mutation origin, DNA repair, epigenetic status and the effect of newly arisen mutations on gamete development are needed.

The different levels of genetic diversity in sex chromosomes and autosomes

Sex chromosomes and autosomes differ in their effective population size, mutation and demography, all of which affect the relative level of genetic diversity within the genome. Moreover, natural selection acts differentially on the two chromosomal categories, for example, because recessive mutations are directly exposed to selection on the single X chromosome of males. Recent genome analyses reveal a heterogeneous picture of the sex-chromosome-to-autosome diversity ratio in different organisms. Reduced X chromosome diversity has been interpreted to reflect demographic features such as bottlenecks and male-biased dispersal, whereas more equal diversity in sex chromosomes and autosomes has been explained by polygynous mating systems.

Evidence that replication-associated mutation alone does not explain between-chromosome differences in substitution rates

Since Haldane first noticed an excess of paternally derived mutations, it has been considered that most mutations derive from errors during germ line replication. Miyata et al. (1987) proposed that differences in the rate of neutral evolution on X, Y, and autosome can be employed to measure the extent of this male bias. This commonly applied method assumes replication to be the sole source of between-chromosome variation in substitution rates. We propose a simple test of this assumption: If true, estimates of the male bias should be independent of which two chromosomal classes are compared. Prior evidence from rodents suggested that this might not be true, but conclusions were limited by a lack of rat Y-linked sequence. We therefore sequenced two rat Y-linked bacterial artificial chromosomes and determined evolutionary rate by comparison with mouse. For estimation of rates we consider both introns and synonymous rates. Surprisingly, for both data sets the prediction of congruent estimates of alpha is strongly rejected. Indeed, some comparisons suggest a female bias with autosomes evolving faster than Y-linked sequence. We conclude that the method of Miyata et al. (1987) has the potential to provide incorrect estimates. Correcting the method requires understanding of the other causes of substitution that might differ between chromosomal classes. One possible cause is recombination-associated substitution bias for which we find some evidence. We note that if, as some suggest, this association is dominantly owing to male recombination, the high estimates of alpha seen in birds is to be expected as Z chromosomes recombine in males.

Mosaics and haemophilia

Some mosaic conditions may affect the haemophilia phenotype. Well-known instances include chromosomal mosaicism because of aneuploidy and pseudo-mosaicism because of varying patterns of X-chromosome inactivation. Chromosomal mosaicism in a chimera is a potential source of phenotypic variation. Gene mosaicism is commonplace. Its pattern and effect depend on the stage of development at which a mutation occurs. Proven or possible genetic mosaicism is an important consideration when predicting the likelihood of transmission of haemophilia to a future generation. A mosaic is an individual who has two or more cell lines, genetically different with regard to chromosomes or genes. As techniques improve and studies accumulate, mosaics are being found to be more common than hitherto believed. Some mosaic conditions may affect the phenotype of haemophilia in males and of the carrier state in females. Cells may be mosaic with regard to chromosomes, as in some instances of aneuploidy, and in chimeras, and in females owing to the pattern of X-chromosome inactivation. Cells may be mosaic with regard to new gene mutations. The pattern of mosaicism depends upon the stage in embryogenesis or in germ-cell formation in which the mutation arose.

Analysis of SINE and LINE repeat content of Y chromosomes in the platypus, Ornithorhynchus anatinus

Monotremes feature an extraordinary sex-chromosome system that consists of five X and five Y chromosomes in males. These sex chromosomes share homology with bird sex chromosomes but no homology with the therian X. The genome of a female platypus was recently completed, providing unique insights into sequence and gene content of autosomes and X chromosomes, but no Y-specific sequence has so far been analysed. Here we report the isolation, sequencing and analysis of ~700 kb of sequence of the non-recombining regions of Y2, Y3 and Y5, which revealed differences in base composition and repeat content between autosomes and sex chromosomes, and within the sex chromosomes themselves. This provides the first insights into repeat content of Y chromosomes in platypus, which overall show similar patterns of repeat composition to Y chromosomes in other species. Interestingly, we also observed differences between the various Y chromosomes, and in combination with timing and activity patterns we provide an approach that can be used to examine the evolutionary history of the platypus sex-chromosome chain.

Purging the genome with sexual selection: reducing mutation load through selection on males

Healthy males are likely to have higher mating success than unhealthy males because of differential expression of condition-dependent traits such as mate searching intensity, fighting ability, display vigor, and some types of exaggerated morphological characters. We therefore expect that most new mutations that are deleterious for overall fitness may also be deleterious for male mating success. From this perspective, sexual selection is not limited to influencing those genes directly involved in exaggerated morphological traits but rather affects most, if not all, genes in the genome. If true, sexual selection can be an important force acting to reduce the frequency of deleterious mutations and, as a result, mutation load. We review the literature and find various forms of indirect evidence that sexual selection helps to eliminate deleterious mutations. However, direct evidence is scant, and there are almost no data available to address a key issue: is selection in males stronger than selection in females? In addition, the total effect of sexual selection on mutation load is complicated by possible increases in mutation rate that may be attributable to sexual selection. Finally, sexual selection affects population fitness not only through mutation load but also through sexual conflict, making it difficult to empirically measure how sexual selection affects load. Several lines of enquiry are suggested to better fill large gaps in our understanding of sexual selection and its effect on genetic load.

The chicken (Gallus gallus) Z chromosome contains at least three nonlinear evolutionary strata

Birds have female heterogamety with Z and W sex chromosomes. These evolved from different autosomal precursor chromosomes than the mammalian X and Y. However, previous work has suggested that the pattern and process of sex chromosome evolution show many similarities across distantly related organisms. Here we show that stepwise restriction of recombination between the protosex chromosomes of birds has resulted in regions of the chicken Z chromosome showing discrete levels of divergence from W homologs (gametologs). The 12 genes analyzed fall into three levels of estimated divergence values, with the most recent divergence (d(S) = 0.18-0.21) displayed by 6 genes in a region on the Z chromosome corresponding to the interval 1-11 Mb of the assembled genome sequence. Another 4 genes show intermediate divergence (d(S) = 0.27-0.38) and are located in the interval 16-53 Mb. Two genes (at positions 42 and 50 Mb) with higher d(S) values are located proximal to the most distal of the 4 genes with intermediate divergence, suggesting an inversion event. The distribution of genes and their divergence indicate at least three evolutionary strata, with estimated times for cessation of recombination between Z and W of 132-150 (stratum 1), 71-99 (stratum 2), and 47-57 (stratum 3) million years ago. An inversion event, or some other form of intrachromosomal rearrangement, subsequent to the formation of strata 1 and 2 has scrambled the gene order to give rise to the nonlinear arrangement of evolutionary strata currently seen on the chicken Z chromosome. These observations suggest that the progressive restriction of recombination is an integral feature of sex chromosome evolution and occurs also in systems of female heterogamety.

Rates of nucleotide substitution in Cornaceae (Cornales)-Pattern of variation and underlying causal factors

Identifying causes of genetic divergence is a central goal in evolutionary biology. Although rates of nucleotide substitution vary among taxa and among genes, the causes of this variation tend to be poorly understood. In the present study, we examined the rate and pattern of molecular evolution for five DNA regions over a phylogeny of Cornus, the single genus of Cornaceae. To identify evolutionary mechanisms underlying the molecular variation, we employed Bayesian methods to estimate divergence times and to infer how absolute rates of synonymous and nonsynonymous substitutions and their ratios change over time. We found that the rates vary among genes, lineages, and through time, and differences in mutation rates, selection type and intensity, and possibly genetic drift all contributed to the variation of substitution rates observed among the major lineages of Cornus. We applied independent contrast analysis to explore whether speciation rates are linked to rates of molecular evolution. The results showed no relationships for individual genes, but suggested a possible localized link between species richness and rate of nonsynonymous nucleotide substitution for the combined cpDNA regions. Furthermore, we detected a positive correlation between rates of molecular evolution and morphological change in Cornus. This was particularly pronounced in the dwarf dogwood lineage, in which genome-wide acceleration in both molecular and morphological evolution has likely occurred.

CpG dinucleotides and the mutation rate of non-CpG DNA

The neutral mutation rate is equal to the base substitution rate when the latter is not affected by natural selection. Differences between these rates may reveal that factors such as natural selection, linkage, or a mutator locus are affecting a given sequence. We examined the neutral base substitution rate by measuring the sequence divergence of approximately 30,000 pairs of inactive orthologous L1 retrotransposon sequences interspersed throughout the human and chimpanzee genomes. In contrast to other studies, we related ortholog divergence to the time (age) that the L1 sequences resided in the genome prior to the chimpanzee and human speciation. As expected, the younger orthologs contained more hypermutable CpGs than the older ones because of their conversion to TpGs (and CpAs). Consequently, the younger orthologs accumulated more CpG mutations than the older ones during the approximately 5 million years since the human and chimpanzee lineages separated. But during this same time, the younger orthologs also accumulated more non-CpG mutations than the older ones. In fact, non-CpG and CpG mutations showed an almost perfect (R2 = 0.98) correlation for approximately 97% of the ortholog pairs. The correlation is independent of G + C content, recombination rate, and chromosomal location. Therefore, it likely reflects an intrinsic effect of CpGs, or mutations thereof, on non-CpG DNA rather than the joint manifestation of the chromosomal environment. The CpG effect is not uniform for all regions of non-CpG DNA. Therefore, the mutation rate of non-CpG DNA is contingent to varying extents on local CpG content. Aside from their implications for mutational mechanisms, these results indicate that a precise determination of a uniform genome-wide neutral mutation rate may not be attainable.

Human-macaque comparisons illuminate variation in neutral substitution rates

The evolutionary distance between human and macaque is particularly attractive for investigating neutral substitution rates, which were calculated as a function of a number of genomic parameters.

Background

The evolutionary distance between human and macaque is particularly attractive for investigating local variation in neutral substitution rates, because substitutions can be inferred more reliably than in comparisons with rodents and are less influenced by the effects of current and ancient diversity than in comparisons with closer primates. Here we investigate the human-macaque neutral substitution rate as a function of a number of genomic parameters.

Results

Using regression analyses we find that male mutation bias, male (but not female) recombination rate, distance to telomeres and substitution rates computed from orthologous regions in mouse-rat and dog-cow comparisons are prominent predictors of the neutral rate. Additionally, we demonstrate that the previously observed biphasic relationship between neutral rate and GC content can be accounted for by properly combining rates at CpG and non-CpG sites. Finally, we find the neutral rate to be negatively correlated with the densities of several classes of computationally predicted functional elements, and less so with the densities of certain classes of experimentally verified functional elements.

Conclusion

Our results suggest that while female recombination may be mainly responsible for driving evolution in GC content, male recombination may be mutagenic, and that other mutagenic mechanisms acting near telomeres, and mechanisms whose effects are shared across mammalian genomes, play significant roles. We also have evidence that the nonlinear increase in rates at high GC levels may be largely due to hyper-mutability of CpG dinucleotides. Finally, our results suggest that the performance of conservation-based prediction methods can be improved by accounting for neutral rates.

Sex: differences in mutation, recombination, selection, gene flow, and genetic drift

In many instances, there are large sex differences in mutation rates, recombination rates, selection, rates of gene flow, and genetic drift. Mutation rates are often higher in males, a difference that has been estimated both directly and indirectly. The higher male mutation rate appears related to the larger number of cell divisions in male lineages but mutation rates also appear gene- and organism-specific. When there is recombination in only one sex, it is always the homogametic sex. When there is recombination in both sexes, females often have higher recombination but there are many exceptions. There are a number of hypotheses to explain the sex differences in recombination. Sex-specific differences in selection may result in stable polymorphisms or for sex chromosomes, faster evolutionary change. In addition, sex-dependent selection may result in antagonistic pleiotropy or sexually antagonistic genes. There are many examples of sex-specific differences in gene flow (dispersal) and a number of adaptive explanations for these differences. The overall effective population size (genetic drift) is dominated by the lower sex-specific effective population size. The mean of the mutation, recombination, and gene flow rates over the two sexes can be used in a population genetics context unless there are sex-specific differences in selection or genetic drift. Sex-specific differences in these evolutionary factors appear to be unrelated to each other. The evolutionary explanations for sex-specific differences for each factor are multifaceted and, in addition, explanations may include chance, nonadaptive differences, or mechanistic, nonevolutionary factors.

Evolution of female preference for younger males

Previous theoretical work has suggested that females should prefer to mate with older males, as older males should have higher fitness than the average fitness of the cohort into which they were born. However, studies in humans and model organisms have shown that as males age, they accumulate deleterious mutations in their germ-line at an ever-increasing rate, thereby reducing the quality of genes passed on to the next generation. Thus, older males may produce relatively poor-quality offspring. To better understand how male age influences female mate preference and offspring quality, we used a genetic algorithm model to study the effect of age-related increases in male genetic load on female mate preference. When we incorporate age-related increases in mutation load in males into our model, we find that females evolve a preference for younger males. Females in this model could determine a male's age, but not his inherited genotype nor his mutation load. Nevertheless, females evolved age-preferences that led them to mate with males that had low mutation loads, but showed no preference for males with respect to their somatic quality. These results suggest that germ-line quality, rather than somatic quality, should be the focus of female preference in good genes models.

A macaque's-eye view of human insertions and deletions: differences in mechanisms

Insertions and deletions (indels) cause numerous genetic diseases and lead to pronounced evolutionary differences among genomes. The macaque sequences provide an opportunity to gain insights into the mechanisms generating these mutations on a genome-wide scale by establishing the polarity of indels occurring in the human lineage since its divergence from the chimpanzee. Here we apply novel regression techniques and multiscale analyses to demonstrate an extensive regional indel rate variation stemming from local fluctuations in divergence, GC content, male and female recombination rates, proximity to telomeres, and other genomic factors. We find that both replication and, surprisingly, recombination are significantly associated with the occurrence of small indels. Intriguingly, the relative inputs of replication versus recombination differ between insertions and deletions, thus the two types of mutations are likely guided in part by distinct mechanisms. Namely, insertions are more strongly associated with factors linked to recombination, while deletions are mostly associated with replication-related features. Indel as a term misleadingly groups the two types of mutations together by their effect on a sequence alignment. However, here we establish that the correct identification of a small gap as an insertion or a deletion (by use of an outgroup) is crucial to determining its mechanism of origin. In addition to providing novel insights into insertion and deletion mutagenesis, these results will assist in gap penalty modeling and eventually lead to more reliable genomic alignments.