Gene flow and selection interact to promote adaptive divergence in regions of low recombination

Adaptation to new environments often occurs in the face of gene flow. Under these conditions, gene flow and recombination can impede adaptation by breaking down linkage disequilibrium between locally adapted alleles. Theory predicts that this decay can be halted or slowed if adaptive alleles are tightly linked in regions of low recombination, potentially favouring divergence and adaptive evolution in these regions over others. Here, we compiled a global genomic data set of over 1,300 individual threespine stickleback from 52 populations and compared the tendency for adaptive alleles to occur in regions of low recombination between populations that diverged with or without gene flow. In support of theory, we found that putatively adaptive alleles (F_ST and d_XY outliers) tend to occur more often in regions of low recombination in populations where divergent selection and gene flow have jointly occurred. This result remained significant when we employed different genomic window sizes, controlled for the effects of mutation rate and gene density, controlled for overall genetic differentiation, varied the genetic map used to estimate recombination and used a continuous (rather than discrete) measure of geographic distance as proxy for gene flow/shared ancestry. We argue that our study provides the first statistical evidence that the interaction of gene flow and selection biases divergence toward regions of low recombination.     

Recombination and the origin of sequence diversity in the DRB MHC class II locus in chamois (Rupicapra spp.)

We examined the evolutionary processes contributing to genetic diversity at the major histocompatibility complex (MHC) class II DRB locus in chamois (Rupicapra spp., subfamily Caprinae). We characterised the pattern of intragenic recombination (or homologous gene conversion) and quantified the amount of recombination in the genealogical history of the two chamois species, Pyrenean chamois (Rupicapra pyrenaica) and Alpine chamois (Rupicapra rupicapra). We found evidence for intragenic recombination, and the estimated amount of population recombination suggests that recombination has been a significant process in generating DRB allelic diversity in the genealogical history of the genus Rupicapra. Moreover, positive selection appears to act on the same peptide-binding residues in both analysed chamois species, but not in identical intensity. Recombination coupled with positive selection drives the rapid evolution at the peptide-binding sites in the MHC class II DRB gene. Many chamois MHC class II DRB alleles are thus much younger than previously assumed.     

Molecular population genetic analysis of the enn subdivision of group A streptococcal emm-like genes: horizontal gene transfer and restricted variation among enn genes

The group A streptococcal emm-like genes, which encode the cell-surface M and M-like proteins, are divided into distinct mrp, emm and enn subdivisions and are clustered together in a region of the chromosome called the vir regulon. In order to understand the mechanisms involved in the evolution of emm-like genes, a 180bp fragment of the 5 variable region of the enn gene was characterized in 31 strains for which emm sequences and multilocus enzyme electrophoretic profiles have been previously determined. The results demonstrate that nucleotide polymorphisms at the enn locus are generated predominantly by point mutations and short deletions or insertions, and that variation among enn and emm genes has arisen by similar mechanisms. However, diversity at the enn locus is restricted in comparison to the emm locus. Moreover, there is strong evidence for intragenic recombination at the enn locus and the pattern of distribution of emm and enn alleles among strains suggests that these genes may be independently acquired by horizontal transfer and recombination from distinct donor strains, thereby generating a mosaic structure for the vir regulon. The results add to a growing body of evidence that horizontal gene transfer has played a major role in the evolution of Streptococcus pyogenes vir regulons.     

Haplotypic Diversity Within the Ovine Calpastatin <Emphasis Type="Italic">(CAST)</Emphasis> Gene

Calpastatin (CAST) is a protein inhibitor that acts specifically on calpains and plays a regulatory role in postmortem beef tenderization and muscle proteolysis. Polymorphisms in the bovine CAST gene have been associated with meat tenderness, but little is known about how the ovine CAST gene may affect sheep meat quality traits. In this study, we selected two parts of the ovine CAST gene that have been previously reported to be polymorphic (region 1—part of intron 5 and exon 6, and region 2—part of intron 12), to investigate haplotype diversity across an extended region of the ovine gene. First, we developed a simple and efficient polymerase chain reaction-single-strand conformational polymorphism (PCR-SSCP) method for genotyping region 2, which allowed the detection of a novel allele as well as the three previously reported alleles. Next, we genotyped both regions 1 and 2 of the ovine CAST gene from a large number of sheep to determine the haplotypes present. Nine different haplotypes were found across this extended region of the ovine CAST gene and four haplotypes were identified that suggested historical recombination events within this gene. Haplotypes are typically more informative than single nucleotide polymorphisms (SNPs) for analyzing associations between genes and complex production traits, such as meat tenderness, but the potential for intragenic recombination within the ovine CAST may make finding associations challenging.     

The molecular mechanism underlying the “rare allele phenomenon” in a subspecific hybrid zone of the California field mouse,Peromyscus californicus

Natural hybrid zones are known to have unusually high levels of novel or otherwise rare electrophoretic variants (the rare allele phenomenon). These variant alleles are most likely the result either of high levels of unique mutations in hybrids or of intragenic recombination between divergent alleles from the parental populations. This study uses DNA sequence comparisons to determine which process has produced a rare allele of the 6-phosphogluconate dehydrogenase (6PGD) gene in a subspecific hybrid zone of the California field mouse (Peromyscus californicus).About 70% of the coding sequence of 6-PGD was cloned and sequenced from three alleles, including two widespread alleles and one rare allele unique to hybrid populations. Sequence comparisons among the three alleles reveal no patterns that would indicate that the variant was formed by intragenic recombination. Instead, the unique allele of 6-PGD studied seems to have developed by the accumulation of base substitutions, which supports the hypothesis of increased mutation rates in hybrids.     

Evolutionary history of an MHC gene in two leporid species: characterisation of <Emphasis Type="Italic">Mhc-DQA</Emphasis> in the European brown hare and comparison with the European rabbit

We surveyed the genetic diversity of the expressed major histocompatibility complex class II DQA locus in natural populations of European brown hares, Lepus europaeus, from Austria and Belgium (267 individuals in total). Based on cDNA sequences, we designed hare-specific primers to amplify the highly variable second exon of the DQA gene. Using cloning–sequencing methodology and capillary electrophoresis single-strand conformation polymorphism, we found ten alleles of the DQA exon 2 locus across these two European regions, of which eight are described for the first time. To search for signals of selection and recombination in the evolution of the DQA gene within the leporids, we augmented our sample with orthologous DQA alleles from the European rabbit, Oryctolagus cuniculus, in order to carry out a species level, species pairwise comparison. We found evidence of recombination in the history of the DQA sequences in leporids with some recombinant alleles bridging the species divide. In both species, selection on peptide binding site codons can be detected, though stronger for the rabbit. This result suggests that there may be a differential selection pressure in the deeper evolutionary history of these two species due to differences in several demographic and ecological traits likely subjecting them to differential selection by parasites. Finally, evolutionary relationships show a widespread and statistically significant intermingling of alleles from the two species. The many macroparasites shared between hares and rabbits may explain this pattern of trans-species polymorphism. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. The nucleotide sequence data reported in this paper have been submitted to Genbank and have been assigned the accession numbers FJ225335–FJ225346.     

Two-locus,fourth-order gene frequency moments: Implications for the variance of squared linkage disequilibrium and the variance of homozygosity

Identity coefficients are used to construct a sufficient set of equations to determine the fourth-order moments of gene frequencies for two linked loci. This allows the variance of the expected squared linkage disequilibrium to be found. It is shown that the coefficient of variation is generally greater than one and if the mutation rate is small, the standard deviation is more than four times the size of the mean. This demonstrates that squared linkage disequilibrium is a highly variable quantity. The variance of homozygosity for a gene which consists of two sites can also be obtained. Recombination between these sites increases the variance of homozygosity, suggesting that intragenic recombination significantly changes all the expected moments of gene frequencies if 4Nμ > 1.0 and r > μ (where N is the population size, μ is the mutation rate of the gene to neutral alleles, and r is the recombination rate between two sites within the gene).     

Direct selection of mutants influencing gene conversion in the yeast Schizosaccharomyces pombe

Summary In Schizosaccharomyces pombe, a suppressor-active mutation at the anticodon site of the tRNA _Ser ^UCA gene sup3 leads to opal (UGA)-specific suppression. Second-site mutations (rX) in sup3 inactivate the suppressor. The sup3-UGA, rX double mutants are genetically unstable in meiotic selfings, due to the intergenic transfer of information between sup3 and the unlinked genes sup9 and sup12 (Hofer et al. 1979; Munz and Leupold 1981; Munz et al. 1982). These three genes have considerable sequence homology over about 200 base pairs (Hottinger et al. 1982).Mutants showing a decrease or an increase of the meiotic instability at sup3 have been selected. One mutation (rec3-8) increases both the genetic instability and the frequency of intragenic recombination in sup3 by one order of magnitude. It has no effect on the stability of the nonsense alleles arg1-230 (UAA), ade6-704 and ura1-61 (UGA) or on the frequency of crossing-over between sup3 and the closely linked gene cdc8.The existence of a common genetic control over intragenic recombination and genetic instability at sup3 provides a direct way of selecting for rec mutants in homothallic haploid strains of S. pombe carrying a suppressor-inactive allele of sup3. It also supports the hypothesis that the instability of mutant alleles of this gene is due to chromosome mispairing at meiosis allowing sup3 to pair with sup9 or sup12 and then to undergo recombination by gene conversion restoring the suppressor-active allele sup3-UGA from the suppressor-inactive allele sup3-UGA,rX. Two mutations (rec2-5 and rec5-11) have no effect on intragenic recombination, but considerably reduce the meiotic instability of the sup3-UGA,rX alleles. They may suppress illegitimate pairing between sup3 and sup9 or sup12.     

Inter- and Intralocus Recombination Drive MHC Class IIB Gene Diversification in a Teleost, the Three-Spined Stickleback Gasterosteus aculeatus

The mutational mechanism underlying the striking diversity in MHC (major histocompatibility complex) genes in vertebrates is still controversial. In order to evaluate the role of inter- and intragenic recombination in MHC gene diversification, we examined patterns of nucleotide polymorphism across an exon/intron boundary in a sample of 31 MHC class IIB sequences of three-spined stickleback (Gasterosteus aculeatus). MHC class IIB genes of G. aculeatus were previously shown to be under diversifying (positive) selection in mate choice and pathogen selection experiments. Based on recoding of alignment gaps, complete intron 2 sequences were grouped into three clusters using maximum-parsimony analysis. Two of these groups had >90% bootstrap support and were tentatively assigned single locus status. Intron nucleotide diversity within and among loci was low (p-distance within and among groups = 0.016 and 0.019, respectively) and fourfold lower than the rate of silent mutations in exon 2, suggesting that noncoding regions are homogenized by frequent interlocus recombination. A substitution analysis using GENECONV revealed as many intergenic conversion events as intragenic ones. Recombination between loci may explain the occurrence of sequence variants that are particularly divergent, as is the case in three-spined stickleback, with nucleotide diversity attaining d_N = 0.39 (peptide-binding residues only). For both MHC class II loci we also estimated the amount of intragenic recombination as population rate (4N_er) under the coalescent and found it to be approximately three times higher compared to point mutations (Watterson estimate per gene, 4N_eμ). Nonindependence of molecular evolution across loci and frequent recombination suggest that MHC class II genes of bony fish may follow different evolutionary dynamics than those of mammals. Our finding of widespread recombination suggests that phylogenies of MHC genes should not be based on coding segments but rather on noncoding introns. [Reviewing Editor: Dr. Richard Kliman]     

Evidence for Recombination in the Microcystin Synthetase (<Emphasis Type="Italic">mcy</Emphasis>) Genes ofToxic Cyanobacteria <Emphasis Type="Italic">Microcystis</Emphasis><Emphasis Type="Bold">spp</Emphasis>

Recombination has been suggested to be an important factor for the genetic variation of bacterial genes, but few studies have dealt with intragenic recombination between the same or closely related species of cyanobacteria. Here we provide strong evidence for recombination in the microcystin synthetase (mcy) gene cluster of the toxic cyanobacteria Microcystis spp. This gene cluster contains 10 genes (mcyA to J) that encode a mixed polyketide synthase (PKS)/nonribosomal peptide synthetase (NRPS) complex. mcy gene sequences were determined for four selected regions (within mcyA, D, G, and J) within the mcy gene cluster from 1 Canadian and 10 Asian toxic Microcystis and compared with previously published mcy sequences. Split decomposition analysis indicated a reticulate phylogeny of mcyA, and several potential recombination tracts of mcyA were identified by the RDP analysis and a runs test implemented in GENECONV. In contrast, no recombination was detected in the mcyD, G, and J sequences. However, discrepancies among the four mcy gene genealogies were evident from the results of independent split decomposition analyses, which were further supported by incongruence length difference (ILD) tests. Taken together, these findings suggest that both intragenic and intergenic recombination within the mcy gene cluster contributes to the genetic diversity of the mcy genes of Microcystis spp.This article contains online supplementary material.     

UV-Induced mitotic recombination and its dependence on photoreactivation and liquid holding in the rad6-1 mutant of Saccharomyces cerevisiae

Summary Spontaneous and UV-induced mitotic recombination was compared in diploids homozygous for rad6-1 mutation and in the wild-type strain carrying heterozygous markers for detecting gene conversion (hom2-1, hom2-2) and crossing over (adel, ade2). Diploids homozygous for rad6-1 mutation were characterised by an elevated level of spontaneous and UV-induced mitotic recombination, particularly the intergenic events. Exposure of UV-irradiated strains to visible light resulted in an increased survival and decreased level of mitotic recombination. Liquid holding (LH) differentially affected frequency of mitotic intergenic and intragenic recombination in mutant and wild-type strains, being without any significant effect on cell survival. In a mutant strain intragenic recombination is significantly increased, intergenic only slightly. In the wild-type strain intragenic recombination is slightly decreased but intergenic is not changed by LH. Visible light applied after LH had no effect on survival and mitotic recombination in the wild type, while in the mutant strain photoreactivability of survival was fully preserved and accompanied by a decrease in the frequency of intragenic and intergenic recombination. The results suggest that metabolic pathways responsible for restoring cell survival are independent of or only partly overlapping with those concerning recombination events.     

Allelic dimorphism in a surface antigen gene of the malaria parasite Plasmodium falciparum

Merozoites of the malaria parasite Plasmodium falciparum carry surface proteins processed from a precursor termed p190 or p195. Polymorphism has been reported in this protein. Since the protein is a candidate for a malaria vaccine, it is important to understand the nature of this polymorphism. We have determined the complete nucleotide sequence of the p190 gene from the MAD20 strain (a Papua New Guinea isolate). Comparisons of the gene with that from other strains of P. falciparum allowed us to study the genetic basis of the antigen's polymorphism. The gene consists of sequences distributed in variable blocks, which are separated by conserved or semi-conserved sequences. Variable sequences occur both in regions that code for tripeptide repeats and in regions with no apparent repeats. Interestingly, according to the present data, variable sequences are not widely polymorphic but fall into two distinct types. We argue that the p190 protein is encoded by dimorphic alleles capable of limited genetic exchange and present evidence at the nucleotide level documenting intragenic recombination in Plasmodium.     

Footprints of intragenic recombination at HLA loci

 To evaluate the effect of balancing selection and intragenic recombination (or gene conversion) at six individual HLA loci, synonymous nucleotide diversity in different exon groups is examined within (π_w) and between (π_b) allelic lineages that may be defined by either serological or DNA sequence differences. Both π values are high in exons which encode for the peptide binding region (PBR) and tend to decrease in other exons. The value of π_w is significantly smaller than that of π_b in any exon of any locus. However, even π_w is much greater than nucleotide diversity at non-HLA loci. These observations provide additional strong evidence for the operation of balancing selection in PBR-encoding exons and its indirect effects on polymorphism at linked neighboring regions. It appears that allelic lineages have generally evolved in isolation but the linkage relationships within and between exons are incomplete throughout the long evolutionary history. To quantify intragenic recombination and account for the large discrepancy between the HLA and non-HLA diversity, a population genetics model is analyzed with special reference to the evolution of modern humans. The analysis suggests that the recombination rate between two sites 1000 base pairs apart is about 10^–5 per generation and that the effective size of human populations (equivalent roughly to the number of breeding individuals in a randomly mating population) has dropped from 10⁵ to 10⁴ in most of the Quaternary. One possibility for this reduction is discussed. Received: 11 August 1997 / Revised: 8 October 1997     

Variation of gene conversion and intragenic recombination frequencies in the genome of Ascobolus immersus.

Summary Eighty mutants in 17 ascospore character genes were studied for their conversion patterns. The correlation between conversion pattern and mutagenic origin, previously found in genes b1 and b2 was extended to all the genes studied. Aberrant 4:4 asci were found in most genes irrespective of their conversion frequency. From gene to gene, the conversion frequency showed an almost 100 times variation. The frequency of intragenic recombination also showed sharp variation from gene to gene. The mean conversion frequency and the maximal intragenic recombination frequency were shown to be highly correlated in 5 genes for which these 2 values are known. This correlation was extended to 12 other genes in other Ascomycetes: Saccharomyces cerevisiae, Schizosaccharomyces pombe, Neurospora, and Sordaria. From this study it is concluded that, 1) the probability of hybrid DNA formation undergoes considerable changes according to the region of the genome; 2) the intragenic recombination frequency primarily reflects the frequency of hybrid DNA formation rather than the physical length of the gene; 3) for a given physical distance on the DNA, a similar fraction of the gene conversion events lead to recombination in the 5 Ascomycetes.     

Natural selection and evolution of streptococcal virulence genes involved in tissue-specific adaptations

The molecular mechanisms underlying niche adaptation in bacteria are not fully understood. Primary infection by the pathogen group A streptococcus (GAS) takes place at either the throat or the skin of its human host, and GAS strains differ in tissue site preference. Many skin-tropic strains bind host plasminogen via the plasminogen-binding group A streptococcal M protein (PAM) present on the cell surface; inactivation of genes encoding either PAM or streptokinase (a plasminogen activator) leads to loss of virulence at the skin. Unlike PAM, which is present in only a subset of GAS strains, the gene encoding streptokinase (ska) is present in all GAS isolates. In this study, the evolution of the virulence genes known to be involved in skin infection was examined. Most genetic diversity within ska genes was localized to a region encoding the plasminogen-docking domain (beta-domain). The gene encoding PAM displayed strong linkage disequilibrium (P < 0.01) with a distinct phylogenetic cluster of the ska beta-domain-encoding region. Yet, ska alleles of distant taxa showed a history of intragenic recombination, and high intrinsic levels of recombination were found among GAS strains having different tissue tropisms. The data suggest that tissue-specific adaptations arise from epistatic coselection of bacterial virulence genes. Additional analysis of ska genes showed that approximately 4% of the codons underwent strong diversifying selection. Horizontal acquisition of one ska lineage from a commensal Streptococcus donor species was also evident. Together, the data suggest that new phenotypes can be acquired through interspecies recombination between orthologous genes, while constrained functions can be preserved; in this way, orthologous genes may provide a rich and ready source for new phenotypes and thereby play a facilitating role in the emergence of new niche adaptations in bacteria.     

Complementation pattern of lexB and recA mutations in Escherichia coli K12; mapping of tif-1, lexB and recA mutations

Summary Three lexB mutations, whose phenotypes have been previously characterized, are studied here in relation to a few recA mutations as to their complementation pattern and relative location.The restoration of resistance to UV-light and to X-rays in the hetero-allelic diploid bacteria was used as a test for dominance and complementation. The wild type allele was always dominant over the mutant allele. Only partial complementation was found between lexB and two recA alleles. There was no complementation between the recA alleles. All the data taken together strongly suggest that the complementations found are intragenic: lexB and recA mutations are in one gene.Mapping of lexB, recA and tif-1 mutations in relation to srl-1 and cysC by phage P1 transduction shows that lexB and the tif-1 mutations form a cluster proximal to srl-1 whereas recA mutations are located at the other extremity of the gene. Variability with temperature of cotransduction frequencies as well as their extended range of values prevent a meaningful calculation of the length of the recA gene.Our hypothesis is that the recA protein has two functional regions called A and B respectively defined at the genetical level by recA and lexB mutations and that it is, in vivo, an oligomeric protein forming a complex with the lexA protein. This complex is postulated to be multifunctional: recombination and control of exonuclease V are effected by the A region while the B region and lexA protein effect induced DNA repair and lysogenic induction.     

Control of recombination ability of vegetative cells in Saccharomyces cerevisiae

Summary In a diploid strain heteroallelic at the ade3 locus, the mitotic intragenic recombination frequency is enhanced ten fold when the cell population is starved for histidine (Hénaut et Luzzati, 1971). By studying simultaneous recombinational events at two independant loci, it is shown that the effect of histidine starvation is most simply explained in term of an increase in the frequency of cells capable of recombination. In these competent cells, intragenic recombination frequencies during mitosis are equal to those found during meiosis. However, the frequency of recombination between the gene and the centromere appears to be lower during mitosis than during meiosis.We believe that histidine starvation in ade3 strains stimulates chromosome pairing, and that there is no fundamental difference between mitotic and meiotic recombination.     

Transitions in the evolution of meiosis

Meiosis may have evolved gradually within the eukaryotes with the earliest forms having a one‐step meiosis. It has been speculated that the putative transition from a one‐step meiosis without recombination to one with recombination may have been stimulated by the invasion of Killer alleles. These imaginary selfish elements are considered to act prior to recombination. They prime for destruction (which occurs after cell division) the half of the cell on the opposite side of the meiotic spindle. Likewise the transition from one‐step to two‐step meiosis might have been stimulated by a subtly different sort of imaginary distorter allele, a SisterKiller. These are proposed to act after recombination. It has yet to be established that the presence of such distorter alleles could induce the transitions in question. To investigate these issues we have analysed the dynamics of a modifier (1) of recombination and (2) of the number of steps of meiosis, as they enter a population with one‐step meiosis. For the modifier of recombination, we find that invasion conditions are very broad and that persistence of Killer and modifier is likely through most parameter space, even when the recombination rate is low. However, if we allow a Killer element to mutate into one that is self‐tolerant, the modifier and the nonself‐tolerant alleles are typically both lost from the population. The modifier of the number of steps can invade if the SisterKiller acts at meiosis II. However, a SisterKiller acting at meiosis I, far from promoting the modifier’s spread, actually impedes it. In the former case the invasion is easiest if there is no recombination. The SisterKiller hypothesis therefore fails to provide a reasonable account of the evolution of two‐step meiosis with recombination. As before, the evolution of self‐tolerance on the part of the selfish element destroys the process. We conclude that the conditions under which SisterKillers promote the evolution of two‐step meiosis are very much more limited than originally considered. We also conclude that there is no universal agreement between ESS and modifier analyses of the same transitions.     

Positive selection and intragenic recombination contribute to high allelic diversity in effector genes of Mycosphaerella fijiensis,causal agent of the black leaf streak disease of banana

Previously, we have determined the nonhost‐mediated recognition of the MfAvr4 and MfEcp2 effector proteins from the banana pathogen Mycosphaerella fijiensis in tomato, by the cognate Cf‐4 and Cf‐Ecp2 resistance proteins, respectively. These two resistance proteins could thus mediate resistance against M. fijiensis if genetically transformed into banana (Musa spp.). However, disease resistance controlled by single dominant genes can be overcome by mutated effector alleles, whose products are not recognized by the cognate resistance proteins. Here, we surveyed the allelic variation within the MfAvr4, MfEcp2, MfEcp2‐2 and MfEcp2‐3 effector genes of M. fijiensis in a global population of the pathogen, and assayed its impact on recognition by the tomato Cf‐4 and Cf‐Ecp2 resistance proteins, respectively. We identified a large number of polymorphisms that could reflect a co‐evolutionary arms race between host and pathogen. The analysis of nucleotide substitution patterns suggests that both positive selection and intragenic recombination have shaped the evolution of M. fijiensis effectors. Clear differences in allelic diversity were observed between strains originating from South‐East Asia relative to strains from other banana‐producing continents, consistent with the hypothesis that M. fijiensis originated in the Asian‐Pacific region. Furthermore, transient co‐expression of the MfAvr4 effector alleles and the tomato Cf‐4 resistance gene, as well as of MfEcp2, MfEcp2‐2 and MfEcp2‐3 and the putative Cf‐Ecp2 resistance gene, indicated that effector alleles able to overcome these resistance genes are already present in natural populations of the pathogen, thus questioning the durability of resistance that can be provided by these genes in the field.