Mitonuclear linkage disequilibrium in human populations

There is extensive evidence from model systems that disrupting associations between co-adapted mitochondrial and nuclear genotypes can lead to deleterious and even lethal consequences. While it is tempting to extrapolate from these observations and make inferences about the human-health effects of altering mitonuclear associations, the importance of such associations may vary greatly among species, depending on population genetics, demographic history and other factors. Remarkably, despite the extensive study of human population genetics, the statistical associations between nuclear and mitochondrial alleles remain largely uninvestigated. We analysed published population genomic data to test for signatures of historical selection to maintain mitonuclear associations, particularly those involving nuclear genes that encode mitochondrial-localized proteins (N-mt genes). We found that significant mitonuclear linkage disequilibrium (LD) exists throughout the human genome, but these associations were generally weak, which is consistent with the paucity of population genetic structure in humans. Although mitonuclear LD varied among genomic regions (with especially high levels on the X chromosome), N-mt genes were statistically indistinguishable from background levels, suggesting that selection on mitonuclear epistasis has not preferentially maintained associations involving this set of loci at a species-wide level. We discuss these findings in the context of the ongoing debate over mitochondrial replacement therapy.     

DOES LINKAGE DISEQUILIBRIUM GENERATE HETEROZYGOSITY‐FITNESS CORRELATIONS IN GREAT REED WARBLERS?

Heterozygosity-fitness correlations (HFCs) at noncoding genetic markers are commonly assumed to reflect fitness effects of heterozygosity at genomewide distributed genes in partially inbred populations. However, in populations with much linkage disequilibrium (LD), HFCs may arise also as a consequence of selection on fitness loci in the local chromosomal vicinity of the markers. Recent data suggest that relatively high levels of LD may prevail in many ecological situations. Consequently, LD may be an important factor, together with partial inbreeding, in causing HFCs in natural populations. In the present study, we evaluate whether LD can generate HFCs in a small and newly founded population of great reed warblers (Acrocephalus arundinaceus). For this purpose dyads of full siblings of which only one individual survived to adult age (i.e., returned to breed at the study area) were scored at 19 microsatellite loci, and at a gene region of hypothesized importance for survival, the major histocompatibility complex (MHC). By examining siblings, we controlled for variation in the inbreeding coefficient and thus excluded genome-wide fitness effects in our analyses. We found that recruited individuals had significantly higher multilocus heterozygosity (MLH), and mean d2 (a microsatellite-specific variable), than their nonrecruited siblings. There was a tendency for the survivors to have a more diverse MHC than the nonsurvivors. Single-locus analyses showed that the strength of the genotype-survival association was especially pronounced at four microsatellite loci. By using genotype data from the entire breeding population, we detected significant LD between five of 162 pairs of microsatellite loci after accounting for multiple tests. Our present finding of a significant within-family multilocus heterozygosity-survival association in a nonequilibrium population supports the view that LD generates HFCs in natural populations.     

Sexual antagonism and meiotic drive cause stable linkage disequilibrium and favour reduced recombination on the X chromosome

Sexual antagonism and meiotic drive are sex‐specific evolutionary forces with the potential to shape genomic architecture. Previous theory has found that pairing two sexually antagonistic loci or combining sexual antagonism with meiotic drive at linked autosomal loci augments genetic variation, produces stable linkage disequilibrium (LD) and favours reduced recombination. However, the influence of these two forces has not been examined on the X chromosome, which is thought to be enriched for sexual antagonism and meiotic drive. We investigate the evolution of the X chromosome under both sexual antagonism and meiotic drive with two models: in one, both loci experience sexual antagonism; in the other, we pair a meiotic drive locus with a sexually antagonistic locus. We find that LD arises between the two loci in both models, even when the two loci freely recombine in females and that driving haplotypes will be enriched for male‐beneficial alleles, further skewing sex ratios in these populations. We introduce a new measure of LD, , which accounts for population allele frequencies and is appropriate for instances where these are sex specific. Both models demonstrate that natural selection favours modifiers that reduce the recombination rate. These results inform observed patterns of congealment found on driving X chromosomes and have implications for patterns of natural variation and the evolution of recombination rates on the X chromosome.     

Significant evidence for linkage disequilibrium over a 5-cM region among Afrikaners

We explore the extent of deviations from Hardy-Weinberg equilibrium (HWE) at a marker locus and linkage disequilibrium (LD) between pairs of marker loci in the Afrikaner population of South Africa. DNA samples were used for genotyping of 23 loci on six chromosomes. The samples were collected from 91 healthy unrelated Afrikaner adults. Exact tests were used to determine evidence for deviations from HWE at a single marker locus or LD between pairs of marker loci. At the 0.05 level of significance, evidence was found for deviation from HWE at only one of the 23 loci. At the same level of significance, LD was found among 8 of the 34 intrachromosomal pairs of loci. On chromosome 21, there was evidence for LD (P = 0.02) between a pair of loci with a genetic distance of 5.51 cM. On chromosome 2, there was evidence for LD between a pair of loci with a genetic distance of 5.28 cM (P = 0.002) and a pair of loci with a genetic distance of 3.68 cM (P = 0.0004). Detailed analysis of LD for one locus pair indicated that only a few of all alleles participated in the LD and that strong LD was most often positive. Our findings indicate that Afrikaans-speaking Afrikaners represent one of those special populations deemed particularly suitable for disequilibrium mapping.     

Genomic breeding value prediction: methods and procedures

Animal breeding faces one of the most significant changes of the past decades - the implementation of genomic selection. Genomic selection uses dense marker maps to predict the breeding value of animals with reported accuracies that are up to 0.31 higher than those of pedigree indexes, without the need to phenotype the animals themselves, or close relatives thereof. The basic principle is that because of the high marker density, each quantitative trait loci (QTL) is in linkage disequilibrium (LD) with at least one nearby marker. The process involves putting a reference population together of animals with known phenotypes and genotypes to estimate the marker effects. Marker effects have been estimated with several different methods that generally aim at reducing the dimensions of the marker data. Nearly all reported models only included additive effects. Once the marker effects are estimated, breeding values of young selection candidates can be predicted with reported accuracies up to 0.85. Although results from simulation studies suggest that different models may yield more accurate genomic estimated breeding values (GEBVs) for different traits, depending on the underlying QTL distribution of the trait, there is so far only little evidence from studies based on real data to support this. The accuracy of genomic predictions strongly depends on characteristics of the reference populations, such as number of animals, number of markers, and the heritability of the recorded phenotype. Another important factor is the relationship between animals in the reference population and the evaluated animals. The breakup of LD between markers and QTL across generations advocates frequent re-estimation of marker effects to maintain the accuracy of GEBVs at an acceptable level. Therefore, at low frequencies of re-estimating marker effects, it becomes more important that the model that estimates the marker effects capitalizes on LD information that is persistent across generations.     

Nucleotide variation at Msn and Alas2, two genes flanking the centromere of the X chromosome in humans

The centromeric region of the X chromosome in humans experiences low rates of recombination over a considerable physical distance. In such a region, the effects of selection may extend to linked sites that are far away. To investigate the effects of this recombinational environment on patterns of nucleotide variability, we sequenced 4581 bp at Msn and 4697 bp at Alas2, two genes situated on either side of the X chromosome centromere, in a worldwide sample of 41 men, as well as in one common chimpanzee and one orangutan. To investigate patterns of linkage disequilibrium (LD) across the centromere, we also genotyped several informative sites from each gene in 120 men from sub-Saharan Africa. By studying X-linked loci in males, we were able to recover haplotypes and study long-range patterns of LD directly. Overall patterns of variability were remarkably similar at these two loci. Both loci exhibited (i) very low levels of nucleotide diversity (among the lowest seen in the human genome); (ii) a strong skew in the distribution of allele frequencies, with an excess of both very-low and very-high-frequency derived alleles in non-African populations; (iii) much less variation in the non-African than in the African samples; (iv) very high levels of population differentiation; and (v) complete LD among all sites within loci. We also observed significant LD between Msn and Alas2 in Africa, despite the fact that they are separated by approximately 10 Mb. These observations are difficult to reconcile with a simple demographic model but may be consistent with positive and/or purifying selection acting on loci within this large region of low recombination.     

Characterization of a divergent chromosome region in the willow warbler Phylloscopus trochilus using avian genomic resources

Genome scans have made it possible to find outlier markers thought to have been influenced by divergent selection in almost any wild population. However, the lack of genomic information in nonmodel species often makes it difficult to associate these markers with certain genes or chromosome regions. Furthermore, the extent of linkage disequilibrium (LD) in the genome will determine the density of markers required to identify the genes under selection. In this study, we investigated a chromosome region in the willow warbler Phylloscopus trochilus surrounding a single marker previously identified in a genome scan. We first located the marker in the assembled genome of another species, the zebra finch Taeniopygia guttata, and amplified surrounding sequences in Fennoscandian willow warblers. Within an investigated chromosome region of 7.3 Mb as mapped to the zebra finch genome, we observed elevated genetic differentiation between a southern and a northern population across a 2.5-Mb interval comprising numerous coding genes. Within the southern and northern populations, higher values of LD were mostly found between SNPs within the same locus, but extended across distantly situated loci when the analyses were restricted to sampling sites showing intermediate allele frequencies of southern and northern alleles. Our study shows that cross-species genome information is a useful resource to obtain candidate sequences adjacent to outlier markers in nonmodel species.     

Allele frequency changes due to hitch-hiking in genomic selection programs

Background

Genomic selection makes it possible to reduce pedigree-based inbreeding over best linear unbiased prediction (BLUP) by increasing emphasis on own rather than family information. However, pedigree inbreeding might not accurately reflect loss of genetic variation and the true level of inbreeding due to changes in allele frequencies and hitch-hiking. This study aimed at understanding the impact of using long-term genomic selection on changes in allele frequencies, genetic variation and level of inbreeding.

Methods

Selection was performed in simulated scenarios with a population of 400 animals for 25 consecutive generations. Six genetic models were considered with different heritabilities and numbers of QTL (quantitative trait loci) affecting the trait. Four selection criteria were used, including selection on own phenotype and on estimated breeding values (EBV) derived using phenotype-BLUP, genomic BLUP and Bayesian Lasso. Changes in allele frequencies at QTL, markers and linked neutral loci were investigated for the different selection criteria and different scenarios, along with the loss of favourable alleles and the rate of inbreeding measured by pedigree and runs of homozygosity.

Results

For each selection criterion, hitch-hiking in the vicinity of the QTL appeared more extensive when accuracy of selection was higher and the number of QTL was lower. When inbreeding was measured by pedigree information, selection on genomic BLUP EBV resulted in lower levels of inbreeding than selection on phenotype BLUP EBV, but this did not always apply when inbreeding was measured by runs of homozygosity. Compared to genomic BLUP, selection on EBV from Bayesian Lasso led to less genetic drift, reduced loss of favourable alleles and more effectively controlled the rate of both pedigree and genomic inbreeding in all simulated scenarios. In addition, selection on EBV from Bayesian Lasso showed a higher selection differential for mendelian sampling terms than selection on genomic BLUP EBV.

Conclusions

Neutral variation can be shaped to a great extent by the hitch-hiking effects associated with selection, rather than just by genetic drift. When implementing long-term genomic selection, strategies for genomic control of inbreeding are essential, due to a considerable hitch-hiking effect, regardless of the method that is used for prediction of EBV.     

Inferring the effects of demography and selection on Drosophila melanogaster populations from a chromosome-wide scan of DNA variation

Identifying regions of the Drosophila melanogaster genome that have been recent targets of positive Darwinian selection will provide evidence for adaptations that have helped this species to colonize temperate habitats. We have begun a search for such genomic regions by analyzing multiple loci (about 250) dispersed across the X chromosome in a putatively ancestral population from East Africa and a derived European population. For both populations we found evidence for past changes in population size. We estimated that a major bottleneck associated with the colonization of Europe occurred about 3,500-16,000 years ago. We also found that while this bottleneck can account for most of the reduction in variation observed in the European sample, there is a deficit of polymorphism in some genomic regions that cannot be explained by demography alone.     

Reduced X-linked rare polymorphism in males in comparison to females of Drosophila melanogaster

Natural selection is assumed to act more strongly on X-linked loci than on autosomal loci because the fitness effect of a recessive mutation on the X chromosome is fully expressed in hemizygous males. Therefore, selection is expected to fix or remove recessive mutations on the X chromosome more efficiently than those on autosomes. However, the assumption that hemizygosity of the X chromosome selectively accelerates changes in allele frequency has not been confirmed directly. To examine this assumption, we investigated current natural selection on X-linked chemoreceptor genes in a natural population of Drosophila melanogaster by comparing nucleotide diversity, linkage disequilibrium (LD), and departure from the neutrality in 4 chemoreceptor genes on 100 X chromosomes each from female and male flies. The general pattern of nucleotide diversity and LD for the genes investigated was similar in females and males. In contrast, males harbored significantly fewer rare polymorphisms defined as singletons and doubletons. When all the gene sequences were concatenated, Tajima's D showed a significant departure from the neutrality in both females and males, whereas Fu and Li's F* value revealed departure only in males. These results suggest that some rare polymorphisms on the X chromosome from females are recessively deleterious and are removed by stronger purifying selection when transferred to hemizygous males.     

Linkage disequilibrium and historical effective population size in the Thoroughbred horse

Many genomic methodologies rely on the presence and extent of linkage disequilibrium (LD) between markers and genetic variants underlying traits of interest, but the extent of LD in the horse has yet to be comprehensively characterized. In this study, we evaluate the extent and decay of LD in a sample of 817 Thoroughbreds. Horses were genotyped for over 50,000 single nucleotide polymorphism (SNP) markers across the genome, with 34,848 autosomal SNPs used in the final analysis. Linkage disequilibrium, as measured by the squared correlation coefficient (r(2)), was found to be relatively high between closely linked markers (>0.6 at 5 kb) and to extend over long distances, with average r(2) maintained above non-syntenic levels for single nucleotide polymorphisms (SNPs) up to 20 Mb apart. Using formulae which relate expected LD to effective population size (N(e)), and assuming a constant actual population size, N(e) was estimated to be 100 in our population. Values of historical N(e), calculated assuming linear population growth, suggested a decrease in N(e) since the distant past, reaching a minimum twenty generations ago, followed by a subsequent increase until the present time. The qualitative trends observed in N(e) can be rationalized by current knowledge of the history of the Thoroughbred breed, and inbreeding statistics obtained from published pedigree analyses are in agreement with observed values of N(e). Given the high LD observed and the small estimated N(e), genomic methodologies such as genomic selection could feasibly be applied to this population using the existing SNP marker set.     

Heterozygosity-fitness correlations within inbreeding classes: local or genome-wide effects?

Marker-based studies of inbreeding may lead to an enhanced understanding of inbreeding depression in natural populations, which is a major concern in conservation genetics. Correlations between marker heterozygosity and variation in fitness-associated traits—‘heterozygosity-fitness correlations’ (HFCs)—are of particular importance and have been widely applied in natural populations. In partially inbred populations, HFCs can be driven by selection against inbred individuals and thus reflect inbreeding depression. However, other explanations for HFCs also exist, such as functional effects of the markers per se or that the markers reveal selection on linked fitness genes due to extended linkage disequilibrium (LD) in the population. Accordingly, HFCs do not only arise in partially inbred populations, they may also occur within inbreeding classes such as families, i.e. in situations when there is no variation in the inbreeding coefficient. In this study we focus on the importance of LD for within-family HFCs, thereby aiming at enhancing our general understanding of HFCs. For non-coding markers, within-family HFCs have been proposed to be caused in two ways: either by ‘local effects’ at linked fitness genes in LD with the markers, or by ‘general effects’ due to a correlation between proportion of heterozygous markers (H _M ) and heterozygosity at genome-wide distributed fitness genes (H _GW ). To evaluate these contrasting hypotheses for within-family HFCs, we analysed simulated data sets of sexually reproducing populations with varying levels of LD. The results confirmed that segregation induces variation in both H _M and H _GW at a fixed level of inbreeding; as expected, the variation in H _M declined with increasing number of markers, whereas the variation in H _GW declined with decreasing LD. However, less intuitively, there was no positive correlation between the variation in H _M and H _GW within inbreeding classes when the local component of H _GW was accounted for (i.e. when the part of the chromosome in LD with the markers was excluded). This strongly suggests that within-family HFCs are not caused by general effects. Instead, our results support the idea that HFCs at a known level of inbreeding can be driven by local effects in populations with high to moderate LD. Note however that we define the local component of H _GW as the part of the chromosomes in LD with the markers. This implies that when LD is high, the local component will consist of a substantial part of the genome and thus provides a rather ‘genome-wide’ view. We caution against routinely interpreting positive HFCs as evidence of inbreeding depression and non-significant HFCs as lack thereof, especially when few markers are used.     

High-resolution mapping of quantitative trait loci for sternopleural bristle number in Drosophila melanogaster.

We have mapped quantitative trait loci (QTL) harboring naturally occurring allelic variation for Drosophila bristle number. Lines with high (H) and low (L) sternopleural bristle number were derived by artificial selection from a large base population. Isogenic H and L sublines were extracted from the selection lines, and populations of X and third chromosome H/L recombinant isogenic lines were constructed in the homozygous low line background. The polymorphic cytological locations of roo transposable elements provided a dense molecular marker map with an average intermarker distance of 4.5 cM. Two X chromosome and six chromosome 3 QTL affecting response to selection for sternopleural bristle number and three X chromosome and three chromosome 3 QTL affecting correlated response in abdominal bristle number were detected using a composite interval mapping method. The average effects of bristle number QTL were moderately large, and some had sex-specific effects. Epistasis between QTL affecting sternopleural bristle number was common, and interaction effects were large. Many of the intervals containing bristle number QTL coincided with those mapped in previous studies. However, resolution of bristle number QTL to the level of genetic loci is not trivial, because the genomic regions containing bristle number QTL often did not contain obvious candidate loci, and results of quantitative complementation tests to mutations at candidate loci affecting adult bristle number were ambiguous.     

Mutant alleles of small effect are primarily responsible for the loss of fitness with slow inbreeding in Drosophila melanogaster.

Multilocus simulation is used to identify genetic models that can account for the observed rates of inbreeding and fitness decline in laboratory populations of Drosophila melanogaster. The experimental populations were maintained under crowded conditions for approximately 200 generations at a harmonic mean population size of Nh approximately 65-70. With a simulated population size of N = 50, and a mean selective disadvantage of homozygotes at individual loci approximately 1-2% or less, it is demonstrated that the mean effective population size over a 200-generation period may be considerably greater than N, with a ratio matching the experimental estimate of Ne/Nh approximately 1.4. The buildup of associative overdominance at electrophoretic marker loci is largely responsible for the stability of gene frequencies and the observed reduction in the rate of inbreeding, with apparent selection coefficients in favor of the heterozygote at neutral marker loci increasing rapidly over the first N generations of inbreeding to values approximately 5-10%. The observed decline in fitness under competitive conditions in populations of size approximately 50 in D. melanogaster therefore primarily results from mutant alleles with mean effects on fitness as homozygotes of sm < or = 0.02. Models with deleterious recessive mutants at the background loci require that the mean selection coefficient against heterozygotes is at most hsm approximately 0.002, with a minimum mutation rate for a single Drosophila autosome 100 cM in length estimated to be in the range 0.05-0.25, assuming an exponential distribution of s. A typical chromosome would be expected to carry at least 100-200 such mutant alleles contributing to the decline in competitive fitness with slow inbreeding.     

Selection strength and hitchhiking around two anti-malarial resistance genes

Neutral mutations may hitchhike to high frequency when they are situated close to sites under positive selection, generating local reductions in genetic diversity. This process is thought to be an important determinant of levels of genomic variation in natural populations. The size of genome regions affected by genetic hitchhiking is expected to be dependent on the strength of selection, but there is little empirical data supporting this prediction. Here, we compare microsatellite variation around two drug resistance genes (chloroquine resistance transporter (pfcrt), chromosome 7, and dihydrofolate reductase (dhfr), chromosome 4) in malaria parasite populations exposed to strong (Thailand) or weak selection (Laos) by anti-malarial drugs. In each population, we examined the point mutations underlying resistance and length variation at 22 (chromosome 4) or 25 (chromosome 7) microsatellite markers across these chromosomes. All parasites from Thailand carried the K76T mutation in pfcrt conferring resistance to chloroquine (CQ) and 2-4 mutations in dhfr conferring resistance to pyrimethamine. By contrast, we found both wild-type and resistant alleles at both genes in Laos. There were dramatic differences in the extent of hitchhiking in the two countries. The size of genome regions affected was smaller in Laos than in Thailand. We observed significant reduction in variation relative to sensitive parasites for 34-64 kb (2-4 cM) in Laos on chromosome 4, compared with 98-137 kb (6-8 cM) in Thailand. Similarly, on chromosome 7, we observed reduced variation for 34-69 kb (2-4 cM) around pfcrt in Laos, but for 195-268 kb (11-16 cM) in Thailand. Reduction in genetic variation was also less extreme in Laos than in Thailand. Most loci were monomorphic in a 12 kb region surrounding both genes on resistant chromosomes from Thailand, whereas in Laos, even loci immediately proximal to selective sites showed some variation on resistant chromosomes. Finally, linkage disequilibrium (LD) decayed more rapidly around resistant pfcrt and dhfr alleles from Laos than from Thailand. These results demonstrate that different realizations of the same selective sweeps may vary considerably in size and shape, in a manner broadly consistent with selection history. From a practical perspective, genomic regions containing resistance genes may be most effectively located by genome-wide association in populations exposed to strong drug selection. However, the lower levels of LD surrounding resistance alleles in populations under weak selection may simplify identification of functional mutations.     

The effect of linkage and population size on inbreeding depression due to mutational load.

Using a stochastic model of a finite population in which there is mutation to partially recessive detrimental alleles at many loci, we study the effects of population size and linkage between the loci on the population mean fitness and inbreeding depression values. Although linkage between the selected loci decreases the amount of inbreeding depression, neither population size nor recombination rate have strong effects on these quantities, unless extremely small values are assumed. We also investigate how partial linkage between the loci that determine fitness affects the invasion of populations by alleles at a modifier locus that controls the selfing rate. In most of the cases studied, the direction of selection on modifiers was consistent with that found in our previous deterministic calculations. However, there was some evidence that linkage between the modifier locus and the selected loci makes outcrossing less likely to evolve; more losses of alleles promoting outcrossing occurred in runs with linkage than in runs with free recombination. We also studied the fate of neutral alleles introduced into populations carrying detrimental mutations. The times to loss of neutral alleles introduced at low frequency were shorter than those predicted for alleles in the absence of selected loci, taking into account the reduction of the effective population size due to inbreeding. Previous studies have been confined to outbreeding populations, and to alleles at frequencies close to one-half, and have found an effect in the opposite direction. It therefore appears that associations between neutral and selected loci may produce effects that differ according to the initial frequencies of the neutral alleles.     

Comparison of the physical and recombination maps of the mouse X chromosome

The locations of five random mouse genomic DNA markers and five cloned genes, including the genes for clotting factors VIII and IX (Cf-8 and Cf-9), Duchenne muscular dystrophy (Dmd), phosphoglycerate kinase-1 (Pgk-1), and alpha-galactosidase (Ags), on the mouse X chromosome were determined by in situ hybridization. The five random DNA markers provide new genetic loci with useful restriction fragment length polymorphisms between mouse strains and species, including one locus close to the centromeric region of the mouse X chromosome. The physical map and the recombination map of these loci on the X chromosome were compared. There was good agreement in the order of loci. Relative distances between loci were consistent along the X chromosome, with the exception of the telomeric end of the long arm, where the recombination fraction observed between loci closely associated on the physical map was higher than that between similarly spaced markers located in the proximal region of the X chromosome. These results are discussed in comparison to the human X-chromosome map.     

Assessing footprints of selection in commercial Atlantic salmon populations using microsatellite data

Relatively large rates of response to traits of economic importance have been observed in different selection experiments in salmon. Several QTL have been mapped in the salmon genome, explaining unprecedented levels of phenotypic variation. Owing to the relatively large selection intensity, individual loci may be indirectly selected, leaving molecular footprints of selection, together with increased inbreeding, as its likely relatives will share the selected loci. We used population differentiation and levels of linkage disequilibrium in chromosomes known to be harbouring QTL for body weight, infectious pancreatic necrosis resistance and infectious salmon anaemia resistance to assess the recent selection history at the genomic level in Atlantic salmon. The results clearly suggest that the marker SSA0343BSFU on chromosome 3 (body weight QTL) showed strong evidence of directional selection. It is intriguing that this marker is physically mapped to a region near the coding sequence of DVL2 , making it an ideal candidate gene to explain the rapid evolutionary response of this chromosome to selection for growth in Salmo salar. Weak evidence of diversifying selection was observed in the QTL associated with infectious pancreatic necrosis and infectious salmon anaemia resistance. Overall, this study showed that artificial selection has produced important changes in the Atlantic salmon genome, validating QTL in commercial salmon populations used for production purposes according to the recent selection history.     

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.     

A microsatellite variability screen for positive selection associated with the "out of Africa" habitat expansion of Drosophila melanogaster

We report a "hitchhiking mapping" study in D. melanogaster, which searches for genomic regions with reduced variability. The study's aim was to identify selective sweeps associated with the "out of Africa" habitat expansion. We scanned 103 microsatellites on chromosome 3 and 102 microsatellites on the X chromosome for reduced variability in non-African populations. When the chromosomes were analyzed separately, the number of loci with a significant reduction in variability only slightly exceeded the expectation under neutrality--six loci on the third chromosome and four loci on the X chromosome. However, non-African populations also have a more pronounced average loss in variability on the X chromosomes as compared to the third chromosome, which suggests the action of selection. Therefore, comparing the X chromosome to the autosome yields a higher number of significantly reduced loci. However, a more pronounced loss of variability on the X chromosome may be caused by demographic events rather than by natural selection. We therefore explored a range of demographic scenarios and found that some of these captured most, but not all aspects of our data. More theoretical work is needed to evaluate how demographic events might differentially affect X chromosomes and autosomes and to estimate the most likely scenario associated with the out of Africa expansion of D. melanogaster.