Experimental evolution across different thermal regimes yields genetic divergence in recombination fraction but no divergence in temperature associated plastic recombination

Phenotypic plasticity is pervasive in nature. One mechanism underlying the evolution and maintenance of such plasticity is environmental heterogeneity. Indeed, theory indicates that both spatial and temporal variation in the environment should favor the evolution of phenotypic plasticity under a variety of conditions. Cyclical environmental conditions have also been shown to yield evolved increases in recombination frequency. Here, we use a panel of replicated experimental evolution populations of D. melanogaster to test whether variable environments favor enhanced plasticity in recombination rate and/or increased recombination rate in response to temperature. In contrast to expectation, we find no evidence for either enhanced plasticity in recombination or increased rates of recombination in the variable environment lines. Our data confirm a role of temperature in mediating recombination fraction in D. melanogaster, and indicate that recombination is genetically and plastically depressed under lower temperatures. Our data further suggest that the genetic architectures underlying plastic recombination and population‐level variation in recombination rate are likely to be distinct.     

Estimation of Fine-Scale Recombination Intensity Variation in the <Emphasis Type="Italic">white–echinus</Emphasis> Interval of <Emphasis Type="Italic">D. melanogaster</Emphasis>

Accurate assessment of local recombination rate variation is crucial for understanding the recombination process and for determining the impact of natural selection on linked sites. In Drosophila, local recombination intensity has been estimated primarily by statistical approaches, by estimating the local slope of the relationship between the physical and genetic maps. However, these estimates are limited in resolution and, as a result, the physical scale at which recombination intensity varies in Drosophila is largely unknown. Although there is some evidence suggesting as much as a 40-fold variation in crossover rate at a local scale in D. pseudoobscura, little is known about the fine-scale structure of recombination rate variation in D. melanogaster. Here we experimentally examine the fine-scale distribution of crossover events in a 1.2-Mb region on the D. melanogaster X chromosome using a classic genetic mapping approach. Our results show that crossover frequency is significantly heterogeneous within this region, varying approximately 3.5-fold. Simulations suggest that this degree of heterogeneity is sufficient to affect levels of standing nucleotide diversity, although the magnitude of this effect is small. We recover no statistical association between empirical estimates of nucleotide diversity and recombination intensity, which is likely due to the limited number of loci sampled in our population genetic data set. However, codon bias is significantly negatively correlated with fine-scale recombination intensity estimates, as expected. Our results shed light on the relevant physical scale to consider in evolutionary analyses relating to recombination rate and highlight the motivations to increase the resolution of the recombination map in Drosophila.     

Analysis of a genetic hitchhiking model, and its application to DNA polymorphism data from Drosophila melanogaster

Begun and Aquadro have demonstrated that levels of nucleotide variation correlate with recombination rate among 20 gene regions from across the genome of Drosophila melanogaster. It has been suggested that this correlation results from genetic hitchhiking associated with the fixation of strongly selected mutants. The hitchhiking process can be described as a series of two-step events. The first step consists of a strongly selected substitution wiping out linked variation in a population; this is followed by a recovery period in which polymorphism can build up via neutral mutations and random genetic drift. Genetic hitchhiking has previously been modeled as a steady-state process driven by recurring selected substitutions. We show here that the characteristic parameter of this steady-state model is alpha v, the product of selection intensity (alpha = 2Ns) and the frequency of beneficial mutations v (where N is population size and s is the selective advantage of the favored allele). We also demonstrate that the steady-state model describes the hitchhiking process adequately, unless the recombination rate is very low. To estimate alpha v, we use the data of DNA sequence variation from 17 D. melanogaster loci from regions of intermediate to high recombination rates. We find that alpha v is likely to be > 1.3 x 10(-8). Additional data are needed to estimate this parameter more precisely. The estimation of alpha v is important, as this parameter determines the shape of the frequency distribution of strongly selected substitutions.      

No Evidence that Infection Alters Global Recombination Rate in House Mice

Recombination rate is a complex trait, with genetic and environmental factors shaping observed patterns of variation. Although recent studies have begun to unravel the genetic basis of recombination rate differences between organisms, less attention has focused on the environmental determinants of crossover rates. Here, we test the effect of one ubiquitous environmental pressure–bacterial infection–on global recombination frequency in mammals. We applied MLH1 mapping to assay global crossover rates in male mice infected with the pathogenic bacterium Borrelia burgdorferi, the causative agent of Lyme Disease, and uninfected control animals. Despite ample statistical power to identify biologically relevant differences between infected and uninfected animals, we find no evidence for a global recombination rate response to bacterial infection. Moreover, broad-scale patterns of crossover distribution, including the number of achiasmate bivalents, are not affected by infection status. Although pathogen exposure can plastically increase recombination in some species, our findings suggest that recombination rates in house mice may be resilient to at least some forms of infection stress. This negative result motivates future experiments with alternative house mouse pathogens to evaluate the generality of this conclusion.     

The genetic architecture of genome‐wide recombination rate variation in allopolyploid wheat revealed by nested association mapping

Recombination affects the fate of alleles in populations by imposing constraints on the reshuffling of genetic information. Understanding the genetic basis of these constraints is critical for manipulating the recombination process to improve the resolution of genetic mapping, and reducing the negative effects of linkage drag and deleterious genetic load in breeding. Using sequence‐based genotyping of a wheat nested association mapping (NAM) population of 2,100 recombinant inbred lines created by crossing 29 diverse lines, we mapped QTL affecting the distribution and frequency of 102 000 crossovers (CO). Genome‐wide recombination rate variation was mostly defined by rare alleles with small effects together explaining up to 48.6% of variation. Most QTL were additive and showed predominantly trans‐acting effects. The QTL affecting the proximal COs also acted additively without increasing the frequency of distal COs. We showed that the regions with decreased recombination carry more single nucleotide polymorphisms (SNPs) with possible deleterious effects than the regions with a high recombination rate. Therefore, our study offers insights into the genetic basis of recombination rate variation in wheat and its effect on the distribution of deleterious SNPs across the genome. The identified trans‐acting additive QTL can be utilized to manipulate CO frequency and distribution in the large polyploid wheat genome opening the possibility to improve the efficiency of gene pyramiding and reducing the deleterious genetic load in the low‐recombining pericentromeric regions of chromosomes.     

Regulation of genome stability by TEL1 and MEC1, yeast homologs of the mammalian ATM and ATR genes

DNA sequence surveys of Drosophila melanogaster populations show a strong positive correlation between the recombination rate experienced by a locus and its level of nucleotide polymorphism. In particular, surveys of the fourth chromosome gene ci(D) show greatly reduced levels of nucleotide variation; this observation was originally interpreted in terms of selective sweeps occurring on the nonrecombining fourth chromosome. Subsequent theoretical work has, however, uncovered several other selective processes that can reduce variation. In this study, we revisit the Drosophila fourth chromosome, investigating variation in 5-6 kb of the gene ankyrin in D. melanogaster and D. simulans. Silent nucleotide site diversity is approximately 5 x 10(-4) for both species, consistent with the previous observations of low variation at ci(D). Given the observed frequency spectra at ankyrin, coalescent simulations indicate that reduced diversity in the region is unlikely to be due to a selective sweep alone. We find evidence for recombinational exchange at this locus, and both species appear to be fixed for an insertion of the transposable element HB in an intron of ankyrin.     

The increase in the proportion of nervous animals bred for catatonia: The participation of central adrenoreceptors in catatonic reactions

Using a large amount of breeding material, the idea of D.K. Belyaev on the role of selection in the appearance of new behavioral and neuronal forms was confirmed. Experiments were performed using rats of the GC (genetics + catatonia) strain, which are prone to passive defensive reactions of cataleptic freezing. At the current breeding stage, elevation of the proportion of so-called ??nervous?? animals was demonstrated, both with respect to the expression of such reactions and their frequency. At this breeding stage, in the brains of GC rats, the mRNA levels of ??1A- and ??2A-adrenoreceptor genes were determined. A decrease of ??1A-adrenoreceptor gene expression in the midbrain and medulla oblongata, along with elevation of ??2A-adreno-receptor gene expression in the frontal cortex was observed. It was suggested that changes in the expression of ??-adrenoreceptor genes could be caused by an increase in the proportion of nervous animals and could contribute to the akinetic behavioral component in GC rats.     

Genome-Wide Fine-Scale Recombination Rate Variation in Drosophila melanogaster

Estimating fine-scale recombination maps of Drosophila from population genomic data is a challenging problem, in particular because of the high background recombination rate. In this paper, a new computational method is developed to address this challenge. Through an extensive simulation study, it is demonstrated that the method allows more accurate inference, and exhibits greater robustness to the effects of natural selection and noise, compared to a well-used previous method developed for studying fine-scale recombination rate variation in the human genome. As an application, a genome-wide analysis of genetic variation data is performed for two Drosophila melanogaster populations, one from North America (Raleigh, USA) and the other from Africa (Gikongoro, Rwanda). It is shown that fine-scale recombination rate variation is widespread throughout the D. melanogaster genome, across all chromosomes and in both populations. At the fine-scale, a conservative, systematic search for evidence of recombination hotspots suggests the existence of a handful of putative hotspots each with at least a tenfold increase in intensity over the background rate. A wavelet analysis is carried out to compare the estimated recombination maps in the two populations and to quantify the extent to which recombination rates are conserved. In general, similarity is observed at very broad scales, but substantial differences are seen at fine scales. The average recombination rate of the X chromosome appears to be higher than that of the autosomes in both populations, and this pattern is much more pronounced in the African population than the North American population. The correlation between various genomic features—including recombination rates, diversity, divergence, GC content, gene content, and sequence quality—is examined using the wavelet analysis, and it is shown that the most notable difference between D. melanogaster and humans is in the correlation between recombination and diversity.     

Kind granddaughters of angry grandmothers: The effect of domestication on vocalization in cross-bred silver foxes

The genetic basis of the effects of domestication has previously been examined in relation to morphological, physiological and behavioural traits, but not for vocalizations. According to Belyaev [Belyaev, D.K., 1979. Destabilizing selection as a factor in domestication. J. Hered. 70, 301-308], directional selection for tame behaviour toward humans resulted in domestication. This hypothesis has been confirmed experimentally on the farm-bred silver fox Vulpes vulpes population that has undergone 45 years of artificial selection for tameness and 35 years of selection for aggressiveness. These foxes, with their precisely known attitudes toward people, provide a means of examining vocal indicators of tameness and aggressiveness to establish the genetic basis for vocal production in canids. We examined vocalizations toward people in foxes selected for tameness and aggressiveness compared to those of three kinds of crosses: Hybrids (Tame × Aggressive), A-Backcrosses (Aggressive × Hybrid) and T-Backcrosses (Tame × Hybrid). We report the effects of selection for tameness on usage and structure of different vocalizations and suggest that vocal indicators for tameness and aggressiveness toward people are discrete phenotypic traits in silver foxes.     

Scrambling eggs: meiotic drive and the evolution of female recombination rates

Theories to explain the prevalence of sex and recombination have long been a central theme of evolutionary biology. Yet despite decades of attention dedicated to the evolution of sex and recombination, the widespread pattern of sex differences in the recombination rate is not well understood and has received relatively little theoretical attention. Here, we argue that female meiotic drivers--alleles that increase in frequency by exploiting the asymmetric cell division of oogenesis--present a potent selective pressure favoring the modification of the female recombination rate. Because recombination plays a central role in shaping patterns of variation within and among dyads, modifiers of the female recombination rate can function as potent suppressors or enhancers of female meiotic drive. We show that when female recombination modifiers are unlinked to female drivers, recombination modifiers that suppress harmful female drive can spread. By contrast, a recombination modifier tightly linked to a driver can increase in frequency by enhancing female drive. Our results predict that rapidly evolving female recombination rates, particularly around centromeres, should be a common outcome of meiotic drive. We discuss how selection to modify the efficacy of meiotic drive may contribute to commonly observed patterns of sex differences in recombination.     

DNA sequence variation and the recombinational landscape in Drosophila pseudoobscura: a study of the second chromosome.

The relationship between rates of recombination and DNA sequence polymorphism was analyzed for the second chromosome of Drosophila pseudoobscura. We constructed integrated genetic and physical maps of this chromosome using molecular markers at 10 loci spanning most of its physical length. The total length of the map was 128.2 cM, almost twice that of the homologous chromosome arm (3R) in D. melanogaster. There appears to be very little centromeric suppression of recombination, and rates of recombination are quite uniform across most of the chromosome. Levels of sequence variation (theta(W), based on the number of segregating sites) at seven loci (tropomyosin 1, Rhodopsin 3, Rhodopsin 1, bicoid, Xanthine dehydrogenase, Myosin light chain 1, and ribosomal protein 49) varied from 0.0036 to 0.0167. Generally consistent with earlier studies, the average estimate of theta(W) at total sites is 1.5-fold higher than that in D. melanogaster, while average theta(W) at silent sites is almost 3-fold higher. These estimates of variation were analyzed in the context of a background selection model under the same parameters of mutation rate and selection as have been proposed for D. melanogaster. It is likely that a significant fraction of the higher level of sequence variation in D. pseudoobscura can be explained by differences in regional rates of recombination rather than a larger species-level effective population size. However, the distribution of variation among synonymous, nonsynonymous, and noncoding sites appears to be quite different between the species, making direct comparisons of neutral variation, and hence inferences about effective population size, difficult. Tajima's D statistics for 6 out of the 7 loci surveyed are negative, suggesting that D. pseudoobscura may have experienced a rapid population expansion in the recent past or, alternatively, that slightly deleterious mutations constitute an important component of standing variation in this species.     

Effect of the major repeat sequence on mitotic recombination in Candida albicans

The major repeat sequence (MRS) is known to play a role in karyotypic variation in Candida albicans. The MRS affects karyotypic variation by expanding and contracting internal repeats, by altering the frequency of chromosome loss, and by serving as a hotspot for chromosome translocation. We proposed that the effects of the MRS on translocation could be better understood by examination of the effect of the MRS on a similar event, mitotic recombination between two chromosome homologs. We examined the frequency of mitotic recombination across an MRS of average size (approximately 50 kb) as well as the rate of recombination in a 325-kb stretch of DNA adjacent to the MRS. Our results indicate that mitotic recombination frequencies across the MRS were not enhanced compared to the frequencies measured across the 325-kb region adjacent to the MRS. Mitotic recombination events were found to occur throughout the 325-kb region analyzed as well as within the MRS itself. This analysis of mitotic recombination frequencies across a large portion of chromosome 5 is the first large-scale analysis of mitotic recombination done in C. albicans and indicates that mitotic recombination frequencies are similar to the rates found in Saccharomyces cerevisiae.     

Nucleotide variation and recombination along the fourth chromosome in Drosophila simulans

The fourth chromosome of Drosophila melanogaster and its sister species are believed to be nonrecombining and have been a model system for testing predictions of the effects of selection on linked, neutral variation. We recently examined nucleotide variation along the chromosome of D. melanogaster and revealed that a low average level of recombination could be associated with considerably high levels of nucleotide variation. In this report, we further investigate the variation along the fourth chromosome of D. simulans. We sequenced 12 gene regions evenly distributed along the fourth chromosome for a worldwide collection of 11 isofemale lines and 5 gene regions in a local population of 10 isofemale lines from South America. In contrast to predictions for regions of very low recombination, these data reveal that the variation levels in many gene regions, including an intron region of the ci gene, vary considerably along the fourth chromosome. Nucleotide diversity ranged from 0.0010 to 0.0074 in 9 gene regions interspersed with several regions of greatly reduced variation. Tests of recombination indicate that the recombination level is not as low as previously thought, likely an order of magnitude higher than that in D. melanogaster. Finally, estimates of the recombination parameters are shown to support a crossover-plus-conversion model.     

Analysis of biological features associated with meiotic recombination hot and cold spots in Saccharomyces cerevisiae

Meiotic recombination is not distributed uniformly throughout the genome. There are regions of high and low recombination rates called hot and cold spots, respectively. The recombination rate parallels the frequency of DNA double-strand breaks (DSBs) that initiate meiotic recombination. The aim is to identify biological features associated with DSB frequency. We constructed vectors representing various chromatin and sequence-based features for 1179 DSB hot spots and 1028 DSB cold spots. Using a feature selection approach, we have identified five features that distinguish hot from cold spots in Saccharomyces cerevisiae with high accuracy, namely the histone marks H3K4me3, H3K14ac, H3K36me3, and H3K79me3; and GC content. Previous studies have associated H3K4me3, H3K36me3, and GC content with areas of mitotic recombination. H3K14ac and H3K79me3 are novel predictions and thus represent good candidates for further experimental study. We also show nucleosome occupancy maps produced using next generation sequencing exhibit a bias at DSB hot spots and this bias is strong enough to obscure biologically relevant information. A computational approach using feature selection can productively be used to identify promising biological associations. H3K14ac and H3K79me3 are novel predictions of chromatin marks associated with meiotic DSBs. Next generation sequencing can exhibit a bias that is strong enough to lead to incorrect conclusions. Care must be taken when interpreting high throughput sequencing data where systematic biases have been documented.     

Hitchhiking under positive Darwinian selection

Positive selection can be inferred from its effect on linked neutral variation. In the restrictive case when there is no recombination, all linked variation is removed. If recombination is present but rare, both deterministic and stochastic models of positive selection show that linked variation hitchhikes to either low or high frequencies. While the frequency distribution of variation can be influenced by a number of evolutionary processes, an excess of derived variants at high frequency is a unique pattern produced by hitchhiking (derived refers to the nonancestral state as determined from an outgroup). We adopt a statistic, H, to measure an excess of high compared to intermediate frequency variants. Only a few high-frequency variants are needed to detect hitchhiking since not many are expected under neutrality. This is of particular utility in regions of low recombination where there is not much variation and in regions of normal or high recombination, where the hitchhiking effect can be limited to a small (<1 kb) region. Application of the H test to published surveys of Drosophila variation reveals an excess of high frequency variants that are likely to have been influenced by positive selection.     

Population genomics in bacteria: a case study of Staphylococcus aureus

We analyzed the genome-wide pattern of single nucleotide polymorphisms (SNPs) in a sample with 12 strains of Staphylococcus aureus. Population structure of S. aureus seems to be complex, and the 12 strains were divided into five groups, named A, B, C, D, and E. We conducted a detailed analysis of the topologies of gene genealogies across the genomes and observed a high rate and frequency of tree-shape switching, indicating extensive homologous recombination. Most of the detected recombination occurred in the ancestral population of A, B, and C, whereas there are a number of small regions that exhibit evidence for homologous recombination with a distinct related species. As such regions would contain a number of novel mutations, it is suggested that homologous recombination would play a crucial role to maintain genetic variation within species. In the A-B-C ancestral population, we found multiple lines of evidence that the coalescent pattern is very similar to what is expected in a panmictic population, suggesting that this population is suitable to apply the standard population genetic theories. Our analysis showed that homologous recombination caused a dramatic decay in linkage disequilibrium (LD) and there is almost no LD between SNPs with distance more than 10 kb. Coalescent simulations demonstrated that a high rate of homologous recombination-a relative rate of 0.6 to the mutation rate with an average tract length of about 10 kb-is required to produce patterns similar to those observed in the S. aureus genomes. Our results call for more research into the evolutionary role of homologous recombination in bacterial populations.     

Connecting recombination,nucleotide diversity,and species divergence in Drosophila

The association between recombination rate and nucleotide diversity provides compelling evidence for the action of natural selection across much of the Drosophila melanogaster genome. This conclusion is further supported by the lack of association between recombination rate and nucleotide divergence between species. However, studies of other species, including other Drosophila, have not always yielded the same results. Our recent study measured these parameters within the D. pseudoobscura species group using next-generation sequencing and high-throughput genotyping technologies. We documented fine-scale variation in crossover rate within D. pseudoobscura, and we observed that crossover variation was strongly associated with nucleotide diversity only when measured at a fine-scale. We also observed associations between crossover rate and sequence differences between D. pseudoobscura and its close relatives. These latter associations could have been driven in part by mutagenic effects associated with double-strand break repair, but we cannot exclude the possibility that it results primarily from shared ancestral polymorphisms. Overall, this work strongly underscores the importance of scale in testing for associations of recombination rate with other parameters, and it brings us one small step closer to understanding the role of natural selection and other evolutionary forces in shaping divergence among genomes.     

Long-term experimental evolution in Escherichia coli. V. Effects of recombination with immigrant genotypes on the rate of bacterial evolution

This study builds upon an earlier experiment that examined the dynamics of mean fitness in evolving populations of Escherichia coli in which mutations were the sole source of genetic variation. During thousands of generations in a constant environment, the rate of improvement in mean fitness of these asexual populations slowed considerably from an initially rapid pace. In this study, we sought to determine whether sexual recombination with novel genotypes would reaccelerate the rate of adaption in these populations. To that end, treatment populations were propagated for an additional 1000 generations in the same environment as their ancestors, but they were periodically allowed to mate with an immigrant pool of genetically distinct Hfr (high frequency recombination) donors. These donors could transfer genes to the resident populations by conjugation, but the donors themselves could not grow in the experimental environment. Control populations were propagated under identical conditions, but in the absence of sexual recombination with the donors. All twelve control populations retained the ancestral alleles at every locus that was scored. In contrast, the sexual recombination treatment yielded dramatic increases in genetic variation. Thus, there was a profound effect of recombination on the rate of genetic change. However, the increased genetic variation in the treatment populations had no significant effect on the rate of adaptive evolution, as measured by changes in mean fitness relative to a common competitor. We then considered three hypotheses that might reconcile these two outcomes: recombination pressure, hitchhiking of recombinant genotypes in association with beneficial mutations, and complex selection dynamics whereby certain genotypes may have a selective advantage only within a particular milieu of competitors. The estimated recombination rate was too low to explain the observed rate of genetic change, either alone or in combination with hitchhiking effects. However, we documented comple x ecological interactions among some recombinant genotypes, suggesting that our method for estimating fitness relative to a common competitor might have underestimated the rate of adaptive evolution in the treatment populations.     

Cytogenetic mapping with centromeric bacterial artificial chromosomes contigs shows that this recombination‐poor region comprises more than half of barley chromosome 3H

Genetic maps are based on the frequency of recombination and often show different positions of molecular markers in comparison to physical maps, particularly in the centromere that is generally poor in meiotic recombinations. To decipher the position and order of DNA sequences genetically mapped to the centromere of barley (Hordeum vulgare) chromosome 3H, fluorescence in situ hybridization with mitotic metaphase and meiotic pachytene chromosomes was performed with 70 genomic single‐copy probes derived from 65 fingerprinted bacterial artificial chromosomes (BAC) contigs genetically assigned to this recombination cold spot. The total physical distribution of the centromeric 5.5 cM bin of 3H comprises 58% of the mitotic metaphase chromosome length. Mitotic and meiotic chromatin of this recombination‐poor region is preferentially marked by a heterochromatin‐typical histone mark (H3K9me2), while recombination enriched subterminal chromosome regions are enriched in euchromatin‐typical histone marks (H3K4me2, H3K4me3, H3K27me3) suggesting that the meiotic recombination rate could be influenced by the chromatin landscape.