Contrasting population genetic structures using allozymes and the inversion polymorphism in Drosophila buzzatii

Second chromosome inversion and genotypic frequencies at seven allozyme loci, differentially associated with inversions, were determined in seven natural populations of Drosophila buzzatii. The patterns of variation of allozymes and the inversion polymorphisms were significantly different, indicating the role of adaptive differentiation for the latter. Moreover, the patterns of population structure varied among allozyme loci, suggesting the operation of diversifying selection for certain loci. Differentiation was negligible for Leucyl‐amino peptidase (Lap) and Peptidase‐2 (Pep‐2), low to moderate for Aldehyde oxidase (Aldox), Peptidase‐1 (Pep‐1) and Esterase‐1 (Est‐1) and high for Esterase‐2 (Est‐2) and Xanthine dehydrogenase (Xdh). Significant linkage disequilibria were detected between inversions and Aldox, Est‐1, Est‐2 and Xdh. Multiple regression analyses of inversion and allele frequencies on environmental variables revealed the existence of clines for inversions, Est‐1, Est‐2, Xdh and Aldox along altitudinal, latitudinal and/or climatic gradients. Tests using conditional allele frequencies showed that Est‐1 and Aldox clines could be accounted for by hitchhiking with inversions, whereas natural selection should be invoked to explain the clines observed for Est‐2 and Xdh.     

Evolution of putative barrier loci at an intermediate stage of speciation with gene flow in campions (Silene)

Understanding the origin of new species is a central goal in evolutionary biology. Diverging lineages often evolve highly heterogeneous patterns of genetic differentiation; however, the underlying mechanisms are not well understood. We investigated evolutionary processes governing genetic differentiation between the hybridizing campions Silene dioica (L.) Clairv. and S. latifolia Poiret. Demographic modelling indicated that the two species diverged with gene flow. The best‐supported scenario with heterogeneity in both migration rate and effective population size suggested that a small proportion of the loci evolved without gene flow. Differentiation (F _ST) and sequence divergence (d _XY) were correlated and both tended to peak in the middle of most linkage groups, consistent with reduced gene flow at highly differentiated loci. Highly differentiated loci further exhibited signatures of selection. In between‐species population pairs, isolation by distance was stronger for genomic regions with low between‐species differentiation than for highly differentiated regions that may contain barrier loci. Moreover, differentiation landscapes within and between species were only weakly correlated, suggesting that linked selection due to shared recombination and gene density landscapes is not the dominant determinant of genetic differentiation in these lineages. Instead, our results suggest that divergent selection shaped the genomic landscape of differentiation between the two Silene species, consistent with predictions for speciation in the face of gene flow.     

GENETIC HITCHHIKING AND THE DYNAMIC BUILDUP OF GENOMIC DIVERGENCE DURING SPECIATION WITH GENE FLOW

A major issue in evolutionary biology is explaining patterns of differentiation observed in population genomic data, as divergence can be due to both direct selection on a locus and genetic hitchhiking. “Divergence hitchhiking” (DH) theory postulates that divergent selection on a locus reduces gene flow at physically linked sites, facilitating the formation of localized clusters of tightly linked, diverged loci. “Genome hitchhiking” (GH) theory emphasizes genome‐wide effects of divergent selection. Past theoretical investigations of DH and GH focused on static snapshots of divergence. Here, we used simulations assessing a variety of strengths of selection, migration rates, population sizes, and mutation rates to investigate the relative importance of direct selection, GH, and DH in facilitating the dynamic buildup of genomic divergence as speciation proceeds through time. When divergently selected mutations were limiting, GH promoted divergence, but DH had little measurable effect. When populations were small and divergently selected mutations were common, DH enhanced the accumulation of weakly selected mutations, but this contributed little to reproductive isolation. In general, GH promoted reproductive isolation by reducing effective migration rates below that due to direct selection alone, and was important for genome‐wide “congealing” or “coupling” of differentiation (F_ST) across loci as speciation progressed.     

Genome‐wide association analyses reveal polygenic genomic architecture underlying divergent shell morphology in Spanish Littorina saxatilis ecotypes

Gene flow between diverging populations experiencing dissimilar ecological conditions can theoretically constrain adaptive evolution. To minimize the effect of gene flow, alleles underlying traits essential for local adaptation are predicted to be located in linked genome regions with reduced recombination. Local reduction in gene flow caused by selection is expected to produce elevated divergence in these regions. The highly divergent crab‐adapted and wave‐adapted ecotypes of the marine snail Littorina saxatilis present a model system to test these predictions. We used genome‐wide association (GWA) analysis of geometric morphometric shell traits associated with microgeographic divergence between the two L. saxatilis ecotypes within three separate sampling sites. A total of 477 snails that had individual geometric morphometric data and individual genotypes at 4,066 single nucleotide polymorphisms (SNPs) were analyzed using GWA methods that corrected for population structure among the three sites. This approach allowed dissection of the genomic architecture of shell shape divergence between ecotypes across a wide geographic range, spanning two glacial lineages. GWA revealed 216 quantitative trait loci (QTL) with shell size or shape differences between ecotypes, with most loci explaining a small proportion of phenotypic variation. We found that QTL were evenly distributed across 17 linkage groups, and exhibited elevated interchromosomal linkage, suggesting a genome‐wide response to divergent selection on shell shape between the two ecotypes. Shell shape trait‐associated loci showed partial overlap with previously identified outlier loci under divergent selection between the two ecotypes, supporting the hypothesis of diversifying selection on these genomic regions. These results suggest that divergence in shell shape between the crab‐adapted and wave‐adapted ecotypes is produced predominantly by a polygenic genomic architecture with positive linkage disequilibrium among loci of small effect.     

Tracking changes in chromosomal arrangements and their genetic content during adaptation

There is considerable evidence for an adaptive role of inversions, but how their genetic content evolves and affects the subsequent evolution of chromosomal polymorphism remains controversial. Here, we track how life‐history traits, chromosomal arrangements and 22 microsatellites, within and outside inversions, change in three replicated populations of Drosophila subobscura for 30 generations of laboratory evolution since founding from the wild. The dynamics of fitness‐related traits indicated adaptation to the new environment concomitant with directional evolution of chromosomal polymorphism. Evidence of selective changes in frequency of inversions was obtained for seven of 23 chromosomal arrangements, corroborating a role for inversions in adaptation. The evolution of linkage disequilibrium between some microsatellites and chromosomes suggested that adaptive changes in arrangements involved changes in their genetic content. Several microsatellite alleles increased in frequency more than expected by drift in targeted inversions in all replicate populations. In particular, there were signs of selection in the O₃₊₄ arrangement favouring a combination of alleles in two loci linked to the inversion and changing along with it, although the lack of linkage disequilibrium between these loci precludes epistatic selection. Seven other alleles increased in frequency within inversions more than expected by drift, but were not in linkage disequilibrium with them. Possibly these alleles were hitchhiking along with alleles under selection that were not specific to those inversions. Overall, the selection detected on the genetic content of inversions, despite limited coverage of the genome, suggests that genetic changes within inversions play an important role in adaptation.     

A long-term study on interactions between the Adh and αGpdh allozyme polymorphisms and the chromosomal inversion In(2L)t in a seminatural population of D. melanogaster

The Adh and αGpdh allozyme loci (both located on the second chromosome) showed considerable fluctuations in allele frequencies in a seminatural population of Drosophila melanogaster during 1972–97. Both long-term and short-term fluctuations were observed. The short-term fluctuations occurred within almost all years and comparison of allele frequencies between winters and summers showed significantly higher Adh^S (P < 0.001) and αGpdh^F (P < 0.01) allele frequencies in summers. Frequencies of these alleles were significantly positively correlated with environmental temperature, suggesting the adaptive significance of these allozyme polymorphisms. Frequency changes of the Odh locus (located on the third chromosome) showed no seasonal pattern and were not correlated with environmental temperature. Almost all short-term and long-term increases in Adh^S frequency were accompanied by a corresponding decrease in αGpdh^S frequency (r = –0.82, P < 0.001) and vice versa. Further analysis showed that gametic disequilibria between the Adh and αGpdh loci, which frequently occurred, were due to the presence of inversion In(2L)t located on the same chromosome arm and In(2L)t frequencies were positively correlated with environmental temperature. Gametic disequilibria between Adh and Odh and between Odh and αGpdh were hardly observed. Because In(2L)t is exclusively associated with the Adh^S/αGpdh^F allele combination, the observed correlated response in Adh/αGpdh allele frequencies is (at least partly) explained by hitchhiking effects with In(2L)t. This means that the adaptive value of the allozyme polymorphisms has been overestimated by ignoring In(2L)t polymorphism. Fluctuations in Adh allele frequencies are fully explained by selection on In(2L)t polymorphism, whereas we have shown that αGpdh frequency fluctuations are only partly explained by chromosomal hitchhiking, indicating the presence of selective differences among αGpdh genotypes in relation with temperature and independent of In(2L)t. Frequency fluctuations of αGpdh and In(2L)t are consistent with their latitudinal distributions, assuming that temperature is the main environmental factor varying with latitude that causes directly or indirectly these frequency distributions. However, the results of the tropical greenhouse population show no correlation of Adh (independent of In(2L)t) and Odh allele frequencies with environmental temperature, which may indicate that the latitudinal distribution in allele frequencies for these loci is not the result of selection on the F/S polymorphism in a direct way.     

Differential introgression from a sister species explains high FST outlier loci within a mussel species

Scanning genomes for loci with high levels of population differentiation has become a standard of population genetics. F_ST outlier loci are most often interpreted as signatures of local selection, but outliers might arise for many other reasons too often left unexplored. Here, we tried to identify further the history and genetic basis underlying strong differentiation at F_ST outlier loci in a marine mussel. A genome scan of genetic differentiation has been conducted between Atlantic and Mediterranean populations of Mytilus galloprovincialis. The differentiation was low overall (F_ST = 0.03), but seven loci (2%) were strong F_ST outliers. We then analysed DNA sequence polymorphism at two outlier loci. The genetic structure proved to be the consequence of differential introgression of alleles from the sister‐hybridizing species Mytilus edulis. Surprisingly, the Mediterranean population was the most introgressed at these two loci, although the contact zone between the two species is nowadays localized along the Atlantic coasts of France and the British Isles. A historical contact between M. edulis and Mediterranean M. galloprovincialis should have happened during glacial periods. It proved difficult to disentangle two hypotheses: (i) introgression was adaptive, implying edulis alleles have been favoured in Mediterranean populations, or (ii) the genetic architecture of the barrier to edulis gene flow is different between the two M. galloprovincialis backgrounds. Five of the seven outliers between M. galloprovincialis populations were also outliers between M. edulis and Atlantic M. galloprovincialis, which would support the latter hypothesis. Differential introgression across semi‐permeable barriers to gene flow is a neglected scenario to interpret outlying loci that may prove more widespread than anticipated.     

Local Patterns of Nucleotide Polymorphism Are Highly Variable in the Selfing Species Arabidopsis thaliana

Neighboring genes predictably share similar evolutionary histories to an extent delineated by recombination. This correlation should extend across multiple linked genes in a selfing species such as Arabidopsis thaliana due to its low effective recombination rate. To test this prediction, we performed a molecular population genetics analysis of nucleotide polymorphism and divergence in chromosomal regions surrounding four low-diversity loci. Three of these loci, At1g67140, At3g03700, and TERMINAL FLOWER1 (TFL1), have been previously implicated as targets of selection and we would predict stronger correlations in polymorphism between neighboring loci due to genetic hitchhiking around these loci. The remaining locus, At1g04300, was identified in a study of linkage disequilibrium surrounding the CRYPTOCHROME2 (CRY2) locus. Although we found broad valleys of reduced nucleotide variation around two of our focal genes, At1g67140 and At3g03700, all chromosomal regions exhibited extreme variation in the patterns of polymorphism and evolution between neighboring loci. Although three of our four regions contained potential targets of selection, application of the composite-likelihood-ratio test of selection in conjunction with a goodness-of-fit test supports the selection hypothesis only for the region containing At3g03700. The degree of discordance in evolutionary histories between linked loci within each region generally correlated with estimates of recombination and linkage disequilibrium for that region, with the exception of the region containing At1g04300. We discuss the implications of these data for future population genetics analyses and genomics studies in A. thaliana. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Linkage disequilibria and the site frequency spectra in the su(s) and su(w(a)) regions of the Drosophila melanogaster X chromosome

Over the last decade, surveys of DNA sequence variation in natural populations of several Drosophila species and other taxa have established that polymorphism is reduced in genomic regions characterized by low rates of crossing over per physical length. Parallel studies have also established that divergence between species is not reduced in these same genomic regions, thus eliminating explanations that rely on a correlation between the rates of mutation and crossing over. Several theoretical models (directional hitchhiking, background selection, and random environment) have been proposed as population genetic explanations. In this study samples from an African population (n = 50) and a European population (n = 51) were surveyed at the su(s) (1955 bp) and su(w(a)) (3213 bp) loci for DNA sequence polymorphism, utilizing a stratified SSCP/DNA sequencing protocol. These loci are located near the telomere of the X chromosome, in a region of reduced crossing over per physical length, and exhibit a significant reduction in DNA sequence polymorphism. Unlike most previously surveyed, these loci reveal substantial skews toward rare site frequencies, consistent with the predictions of directional hitchhiking and random environment models and inconsistent with the general predictions of the background selection model (or neutral theory). No evidence for excess geographic differentiation at these loci is observed. Although linkage disequilibrium is observed between closely linked sites within these loci, many recombination events in the genealogy of the sampled alleles can be inferred and the genomic scale of linkage disequilibrium, measured in base pairs between sites, is the same as that observed for loci in regions of normal crossing over. We conclude that gene conversion must be high in these regions of low crossing over.     

Evolution at the tip and base of the X chromosome in an African population of Drosophila melanogaster

Hitchhiking effects of advantageous mutations have been invoked to explain reduced polymorphism in regions of low crossing-over in Drosophila. Besides reducing DNA heterozygosity, hitchhiking effects should produce strong linkage disequilibrium and a frequency spectrum skewed toward an excess of rare polymorphisms (compared to the neutral expectation). We measured DNA polymorphism in a Zimbabwe population of D. melanogaster at three loci, yellow, achaete, and suppressor of forked, located in regions of reduced crossing-over. Similar to previously published surveys of these genomic regions in other populations, we observed low levels of nucleotide variability. However, the frequency spectrum was compatible with a neutral model, and there was abundant evidence for recombination in the history of the yellow and ac genes. Thus, some aspects of the data cannot be accounted for by a simple hitchhiking model. An alternative hypothesis, background selection, might be compatible with the observed patterns of linkage disequilibrium and the frequency spectrum. However, this model cannot account for the observed reduction in nucleotide heterozygosity. Thus, there is currently no satisfactory theoretical model for the data from the tip and base of the X chromosome in D. melanogaster.      

Natural selection at linked sites in humans

Theoretical and empirical work indicates that patterns of neutral polymorphism can be affected by linked, selected mutations. Under background selection, deleterious mutations removed from a population by purifying selection cause a reduction in linked neutral diversity. Under genetic hitchhiking, the rise in frequency and fixation of beneficial mutations also reduces the level of linked neutral polymorphism. Here we review the evidence that levels of neutral polymorphism in humans are affected by selection at linked sites. We then discuss four approaches for distinguishing between background selection and genetic hitchhiking based on (i) the relationship between polymorphism level and recombination rate for neutral loci with high mutation rates, (ii) relative levels of variation on the X chromosome and the autosomes, (iii) the frequency distribution of neutral polymorphisms, and (iv) population-specific patterns of genetic variation. Although the evidence for selection at linked sites in humans is clear, current methods and data do not allow us to clearly assess the relative importance of background selection and genetic hitchhiking in humans. These results contrast with those obtained for Drosophila, where the signals of positive selection are stronger.     

Footprints of selection in wild populations of Bicyclus anynana along a latitudinal cline

One of the major questions in ecology and evolutionary biology is how variation in the genome enables species to adapt to divergent environments. Here, we study footprints of thermal selection in candidate genes in six wild populations of the afrotropical butterfly Bicyclus anynana sampled along a c. 3000 km latitudinal cline. We sequenced coding regions of 31 selected genes with known functions in metabolism, pigment production, development and heat shock responses. These include genes for which we expect a priori a role in thermal adaptation and, thus, varying selection pressures along a latitudinal cline, and genes we do not expect to vary clinally and can be used as controls. We identified amino acid substitution polymorphisms in 13 genes and tested these for clinal variation by correlation analysis of allele frequencies with latitude. In addition, we used two F_ST‐based outlier methods to identify loci with higher population differentiation than expected under neutral evolution, while accounting for potentially confounding effects of population structure and demographic history. Two metabolic enzymes of the glycolytic pathway, UGP and Treh, showed clinal variation. The same loci showed elevated population differentiation and were identified as significant outliers. We found no evidence of clines in the pigmentation genes, heat shock proteins and developmental genes. However, we identified outlier loci in more localized parts of the range in the pigmentation genes yellow and black. We discuss that the observed clinal variation and elevated population divergence in UGP and Treh may reflect adaptation to a geographic thermal gradient.     

Genetic hitchhiking in a subdivided population of <Emphasis Type="Italic">Mytilus edulis</Emphasis>

Background  

Few models of genetic hitchhiking in subdivided populations have been developed and the rarity of empirical examples is even more striking. We here provide evidences of genetic hitchhiking in a subdivided population of the marine mussel Mytilus edulis. In the Bay of Biscay (France), a patch of M. edulis populations happens to be separated from its North Sea conspecifics by a wide region occupied only by the sister species M. galloprovincialis. Although genetic differentiation between the two M. edulis regions is largely non-significant at ten marker loci (average F_ST~0.007), a strong genetic differentiation is observed at a single locus (F_ST = 0.25). We validated the outlier status of this locus, and analysed DNA sequence polymorphism in order to identify the nature of the selection responsible for the unusual differentiation.     

Occasional recombination of a selfish X‐chromosome may permit its persistence at high frequencies in the wild

The sex‐ratio X‐chromosome (SR) is a selfish chromosome that promotes its own transmission to the next generation by destroying Y‐bearing sperm in the testes of carrier males. In some natural populations of the fly Drosophila neotestacea, up to 30% of the X‐chromosomes are SR chromosomes. To investigate the molecular evolutionary history and consequences of SR, we sequenced SR and standard (ST) males at 11 X‐linked loci that span the ST X‐chromosome and at seven arbitrarily chosen autosomal loci from a sample of D. neotestacea males from throughout the species range. We found that the evolutionary relationship between ST and SR varies among individual markers, but genetic differentiation between SR and ST is chromosome‐wide and likely due to large chromosomal inversions that suppress recombination. However, SR does not consist of a single multilocus haplotype: we find evidence for gene flow between ST and SR at every locus assayed. Furthermore, we do not find long‐distance linkage disequilibrium within SR chromosomes, suggesting that recombination occurs in females homozygous for SR. Finally, polymorphism on SR is reduced compared to that on ST, and loci displaying signatures of selection on ST do not show similar patterns on SR. Thus, even if selection is less effective on SR, our results suggest that gene flow with ST and recombination between SR chromosomes may prevent the accumulation of deleterious mutations and allow its long‐term persistence at relatively high frequencies.     

Divergent and linked selection shape patterns of genomic differentiation between European and North American Atlantic salmon (Salmo salar)

As populations diverge many processes can shape genomic patterns of differentiation. Regions of high differentiation can arise due to divergent selection acting on selected loci, genetic hitchhiking of nearby loci, or through repeated selection against deleterious alleles (linked background selection); this divergence may then be further elevated in regions of reduced recombination. Atlantic salmon (Salmo salar) from Europe and North America diverged >600,000 years ago and despite some evidence of secondary contact, the majority of genetic data indicate substantial divergence between lineages. This deep divergence with potential gene flow provides an opportunity to investigate the role of different mechanisms that shape the genomic landscape during early speciation. Here, using 184,295 single nucleotide polymorphisms (SNPs) and 80 populations, we investigate the genomic landscape of differentiation across the Atlantic Ocean with a focus on highly differentiated regions and the processes shaping them. We found evidence of high (mean F_ST = 0.26) and heterogeneous genomic differentiation between continents. Genomic regions associated with high trans‐Atlantic differentiation ranged in size from single loci (SNPs) within important genes to large regions (1–3 Mbp ) on four chromosomes (Ssa06, Ssa13, Ssa16 and Ssa19). These regions showed signatures consistent with selection, including high linkage disequilibrium, despite no significant reduction in recombination. Genes and functional enrichment of processes associated with differentiated regions may highlight continental differences in ocean navigation and parasite resistance. Our results provide insight into potential mechanisms underlying differences between continents, and evidence of near‐fixed and potentially adaptive trans‐Atlantic differences concurrent with a background of high genome‐wide differentiation supports subspecies designation in Atlantic salmon.     

Polymorphism pattern at a miniature inverted‐repeat transposable element locus downstream of the domestication gene Teosinte‐branched1 in wild and domesticated pearl millet

Unravelling the mechanisms involved in adaptation to understand plant morphological evolution is a challenging goal. For crop species, identification of molecular causal polymorphisms involved in domestication traits is central to this issue. Pearl millet, a domesticated grass mostly found in semi‐arid areas of Africa and India, is an interesting model to address this topic: the domesticated form shares common derived phenotypes with some other cereals such as a decreased ability to develop basal and axillary branches in comparison with the wild phenotype. Two recent studies have shown that the orthologue of the maize gene Teosinte‐Branched1 in pearl millet (PgTb1) was probably involved in branching evolution during domestication and that a miniature inverted‐repeat transposable element (MITE) of the Tuareg family was inserted in the 3′ untranslated region of PgTb1. For a set of 35 wild and domesticated populations, we compared the polymorphism patterns at this MITE and at microsatellite loci. The Tuareg insertion was nearly absent in the wild populations, whereas a strong longitudinal frequency cline was observed in the domesticated populations. The geographical pattern revealed by neutral microsatellite loci clearly demonstrated that isolation by distance does not account for the existence of this cline. However, comparison of population differentiation at the microsatellite and the MITE loci and analyses of the nucleotide polymorphism pattern in the downstream region of PgTb1 did not show evidence that the cline at the MITE locus has been shaped by selection, suggesting the implication of a neutral process. Alternative hypotheses are discussed.     

Fin whale MDH‐1 and MPI allozyme variation is not reflected in the corresponding DNA sequences

The appeal of genetic inference methods to assess population genetic structure and guide management efforts is grounded in the correlation between the genetic similarity and gene flow among populations. Effects of such gene flow are typically genomewide; however, some loci may appear as outliers, displaying above or below average genetic divergence relative to the genomewide level. Above average population, genetic divergence may be due to divergent selection as a result of local adaptation. Consequently, substantial efforts have been directed toward such outlying loci in order to identify traits subject to local adaptation. Here, we report the results of an investigation into the molecular basis of the substantial degree of genetic divergence previously reported at allozyme loci among North Atlantic fin whale (Balaenoptera physalus) populations. We sequenced the exons encoding for the two most divergent allozyme loci (MDH‐1 and MPI) and failed to detect any nonsynonymous substitutions. Following extensive error checking and analysis of additional bioinformatic and morphological data, we hypothesize that the observed allozyme polymorphisms may reflect phenotypic plasticity at the cellular level, perhaps as a response to nutritional stress. While such plasticity is intriguing in itself, and of fundamental evolutionary interest, our key finding is that the observed allozyme variation does not appear to be a result of genetic drift, migration, or selection on the MDH‐1 and MPI exons themselves, stressing the importance of interpreting allozyme data with caution. As for North Atlantic fin whale population structure, our findings support the low levels of differentiation found in previous analyses of DNA nucleotide loci.     

Outlier loci highlight the direction of introgression in oaks

Loci considered to be under selection are generally avoided in attempts to infer past demographic processes as they do not fit neutral model assumptions. However, opportunities to better reconstruct some aspects of past demography might thus be missed. Here we examined genetic differentiation between two sympatric European oak species with contrasting ecological dynamics (Quercus robur and Quercus petraea) with both outlier (i.e. loci possibly affected by divergent selection between species or by hitchhiking effects with genomic regions under selection) and nonoutlier loci. We sampled 855 individuals in six mixed forests in France and genotyped them with a set of 262 SNPs enriched with markers showing high interspecific differentiation, resulting in accurate species delimitation. We identified between 13 and 74 interspecific outlier loci, depending on the coalescent simulation models and parameters used. Greater genetic diversity was predicted in Q. petraea (a late‐successional species) than in Q. robur (an early successional species) as introgression should theoretically occur predominantly from the resident species to the invading species. Remarkably, this prediction was verified with outlier loci but not with nonoutlier loci. We suggest that the lower effective interspecific gene flow at loci showing high interspecific divergence has better preserved the signal of past asymmetric introgression towards Q. petraea caused by the species' contrasting dynamics. Using markers under selection to reconstruct past demographic processes could therefore have broader potential than generally recognized.     

Characterization of eight polymorphic microsatellite markers in the tarnished plant bug,Lygus lineolaris (Palisot de Beauvois)

A partial genomic library of the tarnished plant bug, Lygus lineolaris, enriched for microsatellite sequences was screened to identify marker loci. Eight polymorphic loci suitable for population genetic studies were identified by screening 192 field‐collected insects. The observed number of alleles ranged from four to 21 with an average of 12.25 (SE ± 1.94) while the effective number of alleles ranged from 1.23 to 11.05 with an average of 4.49 (SE ± 1.15). No linkage disequilibria or significant deviations from Hardy–Weinberg expectations were detected at any of the loci. Seven of the eight L. lineolaris microsatellite loci were transferable to Lygus hesperus.     

Sequence diversity patterns suggesting balancing selection in partially sex‐linked genes of the plant Silene latifolia are not generated by demographic history or gene flow

DNA sequence diversity in genes in the partially sex‐linked pseudoautosomal region (PAR) of the sex chromosomes of the plant Silene latifolia is higher than expected from within‐species diversity of other genes. This could be the footprint of sexually antagonistic (SA) alleles that are maintained by balancing selection in a PAR gene (or genes) and affect polymorphism in linked genome regions. SA selection is predicted to occur during sex chromosome evolution, but it is important to test whether the unexpectedly high sequence polymorphism could be explained without it, purely by the combined effects of partial linkage with the sex‐determining region and the population's demographic history, including possible introgression from Silene dioica. To test this, we applied approximate Bayesian computation‐based model choice to autosomal sequence diversity data, to find the most plausible scenario for the recent history of S. latifolia and then to estimate the posterior density of the most relevant parameters. We then used these densities to simulate variation to be expected at PAR genes. We conclude that an excess of variants at high frequencies at PAR genes should arise in S. latifolia populations only for genes with strong associations with fully sex‐linked genes, which requires closer linkage with the fully sex‐linked region than that estimated for the PAR genes where apparent deviations from neutrality were observed. These results support the need to invoke selection to explain the S. latifolia PAR gene diversity, and encourage further work to test the possibility of balancing selection due to sexual antagonism.